| 2 | 1/1 | 返回列表 |
| 查看: 1640 | 回復(fù): 1 | |||
| 【懸賞金幣】回答本帖問題,作者小毛驢起飛將贈送您 5 個金幣 | |||
小毛驢起飛新蟲 (初入文壇)
|
[求助]
lammps新手想求助一下,fix deform算拉伸過程的應(yīng)力應(yīng)變曲線和fix ave/chunk命令相關(guān)
|
||
|
lammps新手想求助一下,本來是按照網(wǎng)上的案例想用fix deform算拉伸過程(x方向)的應(yīng)力應(yīng)變曲線,但是算出來pxx,pyy,pxx都是隨應(yīng)變增大而減小的趨勢,完全不符合, 然后就想用fix ave/chunk 輸出受力,in文件給出,initial.txt是有輸出fx,fy,fz都是接近于0,但是stress.txt這里就沒有輸出受力的數(shù)據(jù)了,想請大家指導(dǎo)討論下是哪里出錯了,感謝。。 #NVT ensemble for relaxation fix 1 all nvt temp 300 300 100 dump 1 all custom 500 nvt.xyz id type x y z vx vy vz thermo_style custom step temp press vol density etotal epair ebond eangle edihed pxx pyy pzz lx ly lz thermo 500 timestep 0.5 reset_timestep 0 run 50000 unfix 1 undump 1 variable tmp equal "lx" variable L0 equal ${tmp} variable strain equal "(lx - v_L0)/v_L0" variable stressx equal "pxx/10000" dump 1 all custom 500 traj.lammpstrj id type x y z ix iy iz vx vy vz fix 1 all npt temp 300 300 100 y 0 0 1000 z 0 0 1000 drag 2 fix 2 all deform 1 x erate 0.000005 units box remap x fix 3 all print 500 "${strain} ${stressx}" file strain_stress.txt screen no thermo_style custom step temp v_strain v_stressx pxx pyy pzz lx ly lz etotal epair ebond eangle edihed thermo 500 compute st all chunk/atom bin/1d x lower 40 units box fix 11 all ave/chunk 1 1 5 st fx fy fz ave one file stress.txt timestep 0.5 reset_timestep 0 run 20000 unfix 1 unfix 2 unfix 3 undump 1 unfix 11 |
新蟲 (初入文壇)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
一志愿天津大學(xué)材料與化工275求調(diào)劑
10+5
|
穿只靴子 2026-03-07 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 一志愿江南大學(xué)085701環(huán)境工程專碩總分287求調(diào)劑 +5 | 18266118446 2026-03-09 | 5/250 |
|
|
[考研] 288求調(diào)劑 +13 | 王曉陽- 2026-03-09 | 18/900 |
|
|
[考研] 哈工大材料324求調(diào)劑 +6 | 閆旭東 2026-03-10 | 8/400 |
|
|
[考研] 307求調(diào)劑 +6 | 超級伊昂大王 2026-03-10 | 6/300 |
|
|
[考研] 材料工程專碩調(diào)劑自薦信,初試323 +6 | 李白26 2026-03-07 | 6/300 |
|
|
[考研] 材料與化工求調(diào)劑 +11 | 與冬清寧 2026-03-07 | 14/700 |
|
|
[考研] 材料專碩調(diào)劑 +7 | 慕辰123 2026-03-05 | 10/500 |
|
|
[考研] 278求調(diào)劑 +8 | Gale1314 2026-03-06 | 8/400 |
|
|
[考研] 0817化學(xué)工程319求調(diào)劑 +7 | lv945 2026-03-08 | 9/450 |
|
|
[考研] 一志愿山東大學(xué),總分327,英語二79,有論文,有競賽,已過四六級 +3 | 木木目目1 2026-03-09 | 3/150 |
|
|
[考研] 293一志愿華東理工 0817化學(xué)工程與技術(shù) 調(diào)劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
[考研] 0817找調(diào)劑 +6 | kk扛 2026-03-08 | 6/300 |
|
|
[考研] 297求調(diào)劑 +3 | 胡達(dá)靈 2026-03-05 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | 拾柒12。 2026-03-08 | 3/150 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 322分 085600求調(diào)劑,有互聯(lián)網(wǎng)+國金及主持省級大創(chuàng)經(jīng)歷 +4 | 熊境喆 2026-03-05 | 4/200 |
|
|
[考研] 275求調(diào)劑 +4 | 大爆炸難民 2026-03-06 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | 曼殊2266 2026-03-05 | 4/200 |
|