 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1853 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問(wèn)題，作者nulixiaoda將贈(zèng)送您 15 個(gè)金幣 | ||
nulixiaoda新蟲(chóng) (初入文壇)
|
[求助]
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ā)自小木蟲(chóng)Android客戶端 |
回復(fù)此樓» 猜你喜歡
282分材料專業(yè)求調(diào)劑院校
已經(jīng)有22人回復(fù)
-
售SCI一區(qū)T0P文章，我:8.O.5.5.1.O.5.4,科目齊全,可+急
已經(jīng)有5人回復(fù)
-
0857 資源與環(huán)境 285分
已經(jīng)有3人回復(fù)
-
標(biāo)題:撿漏預(yù)警|08工科/09農(nóng)學(xué)調(diào)劑!英語(yǔ)要求低，過(guò)線即有機(jī)會(huì)!
已經(jīng)有12人回復(fù)
-
一志愿安徽大學(xué)材料工程專碩313分，求調(diào)劑的學(xué)校
已經(jīng)有8人回復(fù)
-
313求調(diào)劑
已經(jīng)有4人回復(fù)
-
083000環(huán)境科學(xué)與工程調(diào)劑
已經(jīng)有8人回復(fù)
-
材料調(diào)劑，307分
已經(jīng)有12人回復(fù)
-
070300化學(xué) 280 一志愿太原理工 求調(diào)劑
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)材料與化工296分求調(diào)劑
已經(jīng)有12人回復(fù)
2樓2022-07-04 17:19:20
白馬踏平川
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 121.4
- 帖子: 37
- 在線: 3.1小時(shí)
- 蟲(chóng)號(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
白馬踏平川
新蟲(chóng) (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 121.4
- 帖子: 37
- 在線: 3.1小時(shí)
- 蟲(chóng)號(hào): 29139620
- 注冊(cè): 2022-03-14
- 專業(yè): 原子和分子物理
4樓2024-01-11 17:13:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 nulixiaoda 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 標(biāo)題:撿漏預(yù)警|08工科/09農(nóng)學(xué)調(diào)劑!英語(yǔ)要求低，過(guò)線即有機(jī)會(huì)! +7 | 馬超放煙花 2026-03-07 | 12/600 |
|
|---|---|---|---|---|
|
|
[考研] 320材料與化工，求調(diào)劑 +12 | 鶴遨予卿 2026-03-04 | 14/700 |
|
|
|
[考研] 0703求調(diào)劑 +3 | jtyq001 2026-03-10 | 3/150 |
|
|
|
[考研] 一志愿江南大學(xué)085701環(huán)境工程專碩總分287求調(diào)劑 +3 | 18266118446 2026-03-09 | 3/150 |
|
|
[考研] 一志愿南大化學(xué)339分求調(diào)劑，四六級(jí)已過(guò)，有比賽，有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
|
[考研] 313求調(diào)劑 +4 | Yyt楊1 2026-03-07 | 5/250 |
|
|
|
[考研] 0856求調(diào)劑 +3 | squirtle11 2026-03-07 | 3/150 |
|
|
|
[考研] 一志愿蘇大材料工程專碩293求調(diào)劑 +5 | 鋼鐵大炮 2026-03-04 | 5/250 |
|
|
|
[考研] 一志愿武理314求調(diào)劑 +4 | （￣～￣;） 2026-03-08 | 5/250 |
|
|
|
[考研] 材料調(diào)劑 +7 | ounce. 2026-03-05 | 15/750 |
|
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
|
[考研] 267化工調(diào)劑求助 +7 | 聰少OZ 2026-03-04 | 7/350 |
|
|
|
[考研] 材料277分求調(diào)劑 +13 | 飯飯星球 2026-03-04 | 14/700 |
|
|
|
[考研]
 材料工程調(diào)劑
+6
|
zbcm_zbcm 2026-03-05 | 6/300 |
|
|
|
[考研] 287求調(diào)劑 +3 | 看看我. 2026-03-05 | 6/300 |
|
|
|
[考博] 2026申博自薦 六級(jí)440電催化方向 +4 | 櫻落成影花成雙 2026-03-05 | 4/200 |
|
|
|
[考研] 化工282求調(diào)劑一志愿211 +5 | NA0912 2026-03-05 | 6/300 |
|
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
|
[考研] 能動(dòng)297求調(diào)劑，本科川大 +4 | 邵11 2026-03-04 | 4/200 |
|
|
|
[考研] 一志愿西交大材料學(xué)碩 346 求調(diào)劑 +3 | zju51 2026-03-04 | 3/150 |
|