| 7 | 1/1 | 返回列表 |
| 查看: 1287 | 回復: 6 | ||
[求助]
vmd 已有1人參與
|
|
lammps 輸出的lammpsstrj導進vmd中模型為什么會出現(xiàn)連續(xù)的移動?我是做拉伸模擬 發(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的原子,在每個循環(huán)中都要平移,不是嗎? 樓主為什么不用fix deform拉升呢? http://lammps.sandia.gov/doc/fix_deform.html https://icme.hpc.msstate.edu/med ... anowire_Deformation |
送紅花一朵 |
謝謝,那應該怎么修改一下in文件它才不會發(fā)生平移把displace atom這幾行移出label循環(huán)中?我之前用的是deform拉伸,每次拉伸都不在中間斷裂而且得到數(shù)值不準確所以我就換了一種拉伸方法 發(fā)自小木蟲Android客戶端 |
| 7 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 26考研一志愿中國石油大學(華東)305分求調(diào)劑 +6 | 嘉年新程 2026-03-15 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿重慶大學085700資源與環(huán)境總分308求調(diào)劑 +7 | 墨墨漠 2026-03-20 | 7/350 |
|
|
[基金申請]
學校已經(jīng)提交到NSFC,還能修改嗎?
40+4
|
babangida 2026-03-19 | 9/450 |
|
|
[考研] 316求調(diào)劑 +6 | 梁茜雯 2026-03-19 | 6/300 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 一志愿中國石油大學(華東) 本科齊魯工業(yè)大學 +3 | 石能偉 2026-03-17 | 3/150 |
|
|
[考研] 華東師范大學-071000生物學-293分-求調(diào)劑 +3 | 研究生何瑤明 2026-03-18 | 3/150 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 311求調(diào)劑 +5 | 冬十三 2026-03-18 | 5/250 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 一志愿蘇州大學材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
[考研] 一志愿中南化學(0703)總分337求調(diào)劑 +8 | niko- 2026-03-19 | 9/450 |
|
|
[考研] 一志愿中國海洋大學,生物學,301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑 +7 | 為學666 2026-03-16 | 7/350 |
|
|
[考研] 一志愿福大288有機化學,求調(diào)劑 +3 | 小木蟲200408204 2026-03-18 | 3/150 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|