 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1850 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者nulixiaoda將贈(zèng)送您 15 個(gè)金幣 | ||
nulixiaoda新蟲 (初入文壇)
|
[求助]
lammps拉伸時(shí)變形特別大了還沒(méi)有拉斷 已有1人參與
|
|
|
我在用lammps進(jìn)行水合物的拉伸，但拉伸時(shí)三個(gè)尺寸都發(fā)成變形，而且水合物變形特別大了還是不會(huì)斷裂。導(dǎo)出來(lái)的應(yīng)力應(yīng)變曲線也是不規(guī)律。請(qǐng)問(wèn)各位大佬知不知道這是什么原因呢，小弟在這跪求了，謝謝各位！下面是我的in文件。 #水合物拉伸 #初始化 echo screen units real dimension 3 boundary p p p atom_style full #勢(shì)函數(shù) pair_style lj/cut/coul/long 14 kspace_style ewald/disp 1.0e-4 bond_style harmonic angle_style harmonic dihedral_style harmonic pair_modify mix arithmetic #定義文件 read_data "shw.data" #公共近鄰 neighbor 0.3 bin neigh_modify every 2 delay 10 check yes page 100000 velocity all create 273 12345 mom yes rot no #定義變量 計(jì)算值 #應(yīng)力 variable stress equal "-pxx/10000" #應(yīng)變 variable l_x equal lx variable lx0 equal ${l_x} variable strain equal (lx-v_lx0)/v_lx0 #馳豫 timestep 0.001 fix 1 all npt temp 243 243 0.1 z 50 50 0.1 y 50 50 0.1 drag 0.1 thermo 1000 thermo_style custom step lx ly lz press pxx pyy pzz pe temp etotal thermo_modify lost ignore run 20000 #拉伸 reset_timestep 0 unfix 1 fix 1 all npt temp 243 243 0.1 y 50 50 0.1 z 50 50 0.1 drag 0.1 fix 2 all deform 10 x erate 0.001 units box fix 3 all print 100 "${strain} ${stress}" file strain_stress.txt screen no thermo 1000 thermo_style custom step lx ly lz press pxx pyy pzz pe temp etotal thermo_modify lost ignore dump traj all custom 500 shw.xyz id type x y z dump_modify traj sort id timestep 0.001 run 200000 發(fā)自小木蟲Android客戶端 |
回復(fù)此樓» 猜你喜歡
308求調(diào)劑
已經(jīng)有7人回復(fù)
-
085600 材料與化工 295 求調(diào)劑
已經(jīng)有8人回復(fù)
-
求調(diào)劑材料專碩293
已經(jīng)有6人回復(fù)
-
歡迎211本科同學(xué)，過(guò)A區(qū)國(guó)家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組
已經(jīng)有37人回復(fù)
-
標(biāo)題:撿漏預(yù)警|08工科/09農(nóng)學(xué)調(diào)劑!英語(yǔ)要求低，過(guò)線即有機(jī)會(huì)!
已經(jīng)有10人回復(fù)
-
收調(diào)劑
已經(jīng)有7人回復(fù)
-
332材料求調(diào)劑
已經(jīng)有8人回復(fù)
-
調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿天津大學(xué)材料與化工275求調(diào)劑
已經(jīng)有17人回復(fù)
-
材料工程專碩，一志愿中國(guó)礦業(yè)大學(xué)，總分314，求調(diào)劑
已經(jīng)有3人回復(fù)
2樓2022-07-04 17:19:20
白馬踏平川
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 121.4
- 帖子: 37
- 在線: 3.1小時(shí)
- 蟲號(hào): 29139620
- 注冊(cè): 2022-03-14
- 專業(yè): 原子和分子物理
|
您好我想請(qǐng)教您幾個(gè)問(wèn)題，以下是我的in文件 #npt弛豫 units real timestep 1 #npt弛豫 fix NPT all npt temp 270 270 100 x 98 98 1000 y 98 98 1000 z 0 0 1000 thermo_style custom time atoms step temp press vol lx ly lz enthalpy pe ke etotal density pxx pyy pzz pxy pxz pyz thermo ${thermo} dump 1 all xyz ${thermo} npt_${name}.xyz dump_modify 1 element O H C fix pt all print 200 "${istep} ${etol} ${han} ${T0} ${V} ${PE} ${P0}" screen no file npt_${name}.dat run ${runmd} undump 1 unfix pt unfix NPT reset_timestep 0 #stress-strain variable tmpz equal "lz" variable Lz0 equal ${tmpz} variable strainz equal "(lz - v_Lz0)/v_Lz0" variable stressz equal "-pzz/10000" #z向拉伸 fix 2 all npt temp 270 270 100 x 98 98 1000 y 98 98 1000 fix 3 all deform 100 z erate 0.00000001 remap x units box thermo_style custom time atoms step temp press vol lx ly lz density pxx pyy pzz pxy pxz pyz thermo ${thermo} dump 1 all xyz ${thermo} tension_${name}.xyz dump_modify 1 element O H C #保存應(yīng)力應(yīng)變數(shù)據(jù) fix def3 all print 100 "${strainx} ${stressx} ${strainy} ${stressy} ${strainz} ${stressz}" screen no file stress-strain_${name}.dat run ${runten} 1.您所說(shuō)的x方向拉伸，在NPT時(shí)，使x方向的壓強(qiáng)調(diào)節(jié)為零，是在弛豫階段設(shè)定為0嗎？我是z方向拉伸，這么設(shè)置對(duì)嗎？但是我這么設(shè)置應(yīng)力依舊不是從0開(kāi)始的，我檢查了log文件發(fā)現(xiàn)我的模擬過(guò)程控不住壓，壓力波動(dòng)巨大，是否是這個(gè)原因那？應(yīng)該怎么解決？ 2.我的應(yīng)力應(yīng)變文件的屈服強(qiáng)度很大，有8GPa，這明顯是有問(wèn)題的，這是什么原因那？也是因?yàn)榭刂撇蛔毫?dǎo)致的嗎？ 非常希望能得到您的指點(diǎn)，感謝。 |
3樓2024-01-11 17:06:33
白馬踏平川
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 121.4
- 帖子: 37
- 在線: 3.1小時(shí)
- 蟲號(hào): 29139620
- 注冊(cè): 2022-03-14
- 專業(yè): 原子和分子物理
4樓2024-01-11 17:13:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 nulixiaoda 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0857環(huán)境調(diào)劑 +3 | 熠熠_11 2026-03-10 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 0703化學(xué)調(diào)劑 +3 | 三dd. 2026-03-10 | 3/150 |
|
|
[考研] 304求調(diào)劑（085602一志愿985） +8 | 化工人999 2026-03-09 | 8/400 |
|
|
|
[考研] 08工科 +5 | li李樂(lè)成 2026-03-06 | 5/250 |
|
|
|
[考研] 0817一志愿蘇州大學(xué)280復(fù)試調(diào)劑 +9 | kk扛 2026-03-07 | 9/450 |
|
|
|
[考研] 085600材料與化工，一志愿廣州985，求調(diào)劑 +15 | qqyyaill 2026-03-05 | 15/750 |
|
|
|
[考研] 一志愿南大化學(xué)339分求調(diào)劑，四六級(jí)已過(guò)，有比賽，有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | 鶴遨予卿 2026-03-09 | 3/150 |
|
|
|
[考研] 材料與化工，291分，求調(diào)劑 +13 | 咕嚕咕嚕123123 2026-03-04 | 15/750 |
|
|
|
[考研] 085701環(huán)境工程專碩求調(diào)劑 +5 | 幼稚園大班打手 2026-03-08 | 5/250 |
|
|
|
[考研] 化學(xué)工程求調(diào)劑 +12 | 化工人999 2026-03-04 | 12/600 |
|
|
|
[考研] 材料科學(xué)（0805）338 求調(diào)劑 +7 | xiaokang3286 2026-03-07 | 7/350 |
|
|
|
[考研] 322分 085600求調(diào)劑，有互聯(lián)網(wǎng)+國(guó)金及主持省級(jí)大創(chuàng)經(jīng)歷 +4 | 熊境喆 2026-03-05 | 4/200 |
|
|
|
[考研] 誠(chéng)求調(diào)劑，323分 有專利/科研/實(shí)習(xí)經(jīng)歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
|
[考研] 一志愿南京大學(xué)資源與環(huán)境專業(yè)283分求調(diào)劑 +8 | 追云?? 2026-03-04 | 8/400 |
|
|
|
[考研]
 085600材料與化工 292分 求調(diào)劑
+6
|
程晴之 2026-03-06 | 6/300 |
|
|
|
[考研] 材料考研339求調(diào)劑 +3 | Karry*^_^* 2026-03-04 | 3/150 |
|
|
|
[考研] 一志愿武漢理工大學(xué)-085602-總分296分-求調(diào)劑 +7 | 紫川葡柚 2026-03-04 | 7/350 |
|
|
|
[考研] 材料專碩346求調(diào)劑 +3 | 旺一下 2026-03-04 | 3/150 |
|
|
|
[考研] 281求調(diào)劑 +3 | 我是小小蔥蔥 2026-03-03 | 5/250 |
|