 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 1859 | 回復(fù): 3 | ||
| 【懸賞金幣】回答本帖問題，作者nulixiaoda將贈送您 15 個(gè)金幣 | ||
nulixiaoda新蟲 (初入文壇)
|
[求助]
lammps拉伸時(shí)變形特別大了還沒有拉斷 已有1人參與
|
|
|
我在用lammps進(jìn)行水合物的拉伸，但拉伸時(shí)三個(gè)尺寸都發(fā)成變形，而且水合物變形特別大了還是不會斷裂。導(dǎo)出來的應(yīng)力應(yīng)變曲線也是不規(guī)律。請問各位大佬知不知道這是什么原因呢，小弟在這跪求了，謝謝各位！下面是我的in文件。 #水合物拉伸 #初始化 echo screen units real dimension 3 boundary p p p atom_style full #勢函數(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ù)此樓» 猜你喜歡
售SCI一區(qū)T0P文章，我:8.O.5.5.1.O.5.4,科目齊全,可+急
已經(jīng)有7人回復(fù)
-
（081700）化學(xué)工程與技術(shù)-298分求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中科大材料與化工。353分能調(diào)劑到哪？
已經(jīng)有7人回復(fù)
-
337一志愿華南理工0805材料求調(diào)劑
已經(jīng)有5人回復(fù)
-
伙伴們，祝我生日快樂吧
已經(jīng)有4人回復(fù)
-
301求調(diào)劑
已經(jīng)有3人回復(fù)
-
材料與化工304求B區(qū)調(diào)劑
已經(jīng)有4人回復(fù)
-
讀博申請
已經(jīng)有7人回復(fù)
-
一志愿江南大學(xué)085701環(huán)境工程專碩總分287求調(diào)劑
已經(jīng)有5人回復(fù)
2樓2022-07-04 17:19:20
白馬踏平川
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 121.4
- 帖子: 37
- 在線: 3.1小時(shí)
- 蟲號: 29139620
- 注冊: 2022-03-14
- 專業(yè): 原子和分子物理
|
您好我想請教您幾個(gè)問題，以下是我的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.您所說的x方向拉伸，在NPT時(shí)，使x方向的壓強(qiáng)調(diào)節(jié)為零，是在弛豫階段設(shè)定為0嗎？我是z方向拉伸，這么設(shè)置對嗎？但是我這么設(shè)置應(yīng)力依舊不是從0開始的，我檢查了log文件發(fā)現(xiàn)我的模擬過程控不住壓，壓力波動巨大，是否是這個(gè)原因那？應(yīng)該怎么解決？ 2.我的應(yīng)力應(yīng)變文件的屈服強(qiáng)度很大，有8GPa，這明顯是有問題的，這是什么原因那？也是因?yàn)榭刂撇蛔毫?dǎo)致的嗎？ 非常希望能得到您的指點(diǎn)，感謝。 |
3樓2024-01-11 17:06:33
白馬踏平川
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 121.4
- 帖子: 37
- 在線: 3.1小時(shí)
- 蟲號: 29139620
- 注冊: 2022-03-14
- 專業(yè): 原子和分子物理
4樓2024-01-11 17:13:36
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 nulixiaoda 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中科大材料與化工。353分能調(diào)劑到哪？ +5 | MayUxw1 2026-03-08 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 291分獸醫(yī)專碩 一志愿華中農(nóng) +3 | 種大大太陽 2026-03-05 | 3/150 |
|
|
|
[碩博家園] 深圳大學(xué)碩士招生（2026秋，傳感器方向，僅錄取第一志愿） +4 | xujiaoszu 2026-03-11 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +4 | 鶴遨予卿 2026-03-09 | 4/200 |
|
|
[考研] 086000生物與醫(yī)藥319分求調(diào)劑 +4 | Tolkien 2026-03-07 | 8/400 |
|
|
|
[考研] 收調(diào)劑 +7 | 調(diào)劑的考研學(xué)生 2026-03-10 | 7/350 |
|
|
|
[考研] 一志愿湖師大化學(xué)289求調(diào)劑 +5 | XMCMM3.14159 2026-03-10 | 5/250 |
|
|
|
[考研] 0856材料與化工309分求調(diào)劑 +4 | ZyZy…… 2026-03-10 | 4/200 |
|
|
|
[考研] 337求調(diào)劑 +3 | 睡醒，。 2026-03-09 | 3/150 |
|
|
|
[考研] 一志愿211 材料與化工 280求調(diào)劑 +3 | Sanity蔣 2026-03-08 | 3/150 |
|
|
|
[考研] 安徽農(nóng)業(yè)大學(xué)材料與化學(xué)學(xué)院0856材料招收調(diào)劑 +3 | akakk47 2026-03-05 | 9/450 |
|
|
|
[教師之家] 交大前校長王樹國：現(xiàn)在最先進(jìn)的科技并不在大學(xué)實(shí)驗(yàn)室，而是在企業(yè)研究院 +4 | zju2000 2026-03-08 | 6/300 |
|
|
|
[考研] 085701環(huán)境工程317分求調(diào)劑 +9 | 6汆尼9 2026-03-07 | 9/450 |
|
|
|
[考研] 一志愿武漢理工085601，初試301分，請問能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
|
[考研] 2026調(diào)劑】考試A區(qū)0703化學(xué)類323分 誠求接收 +3 | 卷柏卷柏 2026-03-05 | 4/200 |
|
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
|
[考研] 308求調(diào)劑 +7 | 倘若起風(fēng)了呢 2026-03-05 | 9/450 |
|
|
|
[考研] 306求調(diào)劑 +7 | Bahati 2026-03-05 | 7/350 |
|
|
|
[考研] 考研282分求調(diào)劑，接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|