| 7 | 1/1 | 返回列表 |
| 查看: 1286 | 回復(fù): 6 | ||
[求助]
vmd 已有1人參與
|
|
lammps 輸出的lammpsstrj導(dǎo)進(jìn)vmd中模型為什么會(huì)出現(xiàn)連續(xù)的移動(dòng)?我是做拉伸模擬 發(fā)自小木蟲Android客戶端 |
木蟲 (正式寫手)
木蟲 (正式寫手)
|
# LAMMPS cu units metal dimension 3 boundary s s s atom_style atomic neighbor 2.0 bin neigh_modify every 1 read_data data.cu pair_style hybrid eam/alloy airebo 3.0 1 1 lj/cut 10.0 pair_coeff * * eam/alloy Cu_zhou.eam.alloy NULL Cu pair_coeff * * airebo CH.airebo C NULL pair_coeff 1 2 lj/cut 0.019996 3.225 timestep 0.001 thermo 1000 group lower id <= 133 group force id >= 2539 group boundary union lower force group mobile subtract all lower force velocity mobile create 300 825577 dist gaussian velocity boundary set 0.0 0.0 0.0 units box fix 2 lower setforce 0.0 0.0 0.0 fix 3 force setforce 0.0 0.0 0.0 minimize 1.0e-8 1.0e-8 10000 10000 reset_timestep 0 compute 2 mobile temp/com compute 3 all pe/atom compute 4 all ke/atom compute 5 all coord/atom 3.0 compute s1 all stress/atom NULL compute strec1 mobile reduce sum c_s1[2] variable strecy equal c_strec1/(ly*lz*lx*10000) compute strec2 mobile reduce sum c_s1[1] variable strecx equal c_strec2/(ly*lz*lx*10000) thermo_style custom step c_2 pe ke vol v_strecx v_strecy lx ly lz dump 6 all custom 1000 lashen.lammpstrj id type x y z c_3 fix 10 all nvt temp 300 300 0.01 fix_modify 10 temp 2 label loop variable i loop 300 displace_atoms force move 0.0 0.05 0.0 units box displace_atoms lower move 0.0 -0.05 0.0 units box run 100000 next i jump in.cu loop |
木蟲 (正式寫手)
|
abel loop variable i loop 300 displace_atoms force move 0.0 0.05 0.0 units box displace_atoms lower move 0.0 -0.05 0.0 units box 這幾行放在了循環(huán)里,那么你的group force/lower的原子,在每個(gè)循環(huán)中都要平移,不是嗎? 樓主為什么不用fix deform拉升呢? http://lammps.sandia.gov/doc/fix_deform.html https://icme.hpc.msstate.edu/med ... anowire_Deformation |
送紅花一朵 |
謝謝,那應(yīng)該怎么修改一下in文件它才不會(huì)發(fā)生平移把displace atom這幾行移出label循環(huán)中?我之前用的是deform拉伸,每次拉伸都不在中間斷裂而且得到數(shù)值不準(zhǔn)確所以我就換了一種拉伸方法 發(fā)自小木蟲Android客戶端 |
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考! +3 | lemonzzn 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 0856 271求調(diào)劑 +3 | 展信悅_ 2026-03-21 | 3/150 |
|
|
[考研] 306求調(diào)劑 +4 | chuanzhu川燭 2026-03-18 | 4/200 |
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
[考研] 初始318分求調(diào)劑(有工作經(jīng)驗(yàn)) +3 | 1911236844 2026-03-17 | 3/150 |
|
|
[考研] 280求調(diào)劑 +7 | 咕嚕曉曉 2026-03-18 | 8/400 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠(chéng)心求調(diào)劑 +8 | 吃吃吃才有意義 2026-03-19 | 8/400 |
|
|
[考研] 295求調(diào)劑 +4 | 一志愿京區(qū)211 2026-03-18 | 6/300 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 334求調(diào)劑 +3 | 志存高遠(yuǎn)意在機(jī)?/a> 2026-03-16 | 3/150 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|