| 查看: 670 | 回复: 2 | ||
[求助]
如何输出单原子应力数据 已有1人参与
|
|
原本我是想输出单原子应力数据，所以在dump命令最后加了v_Press，但是无法运行，所以请大神帮我看看哪里有问题？？万分感谢！！！！ # tensile testing boundary s s p units metal atom_style atomic # construct [001] [110] [111] lattice orientation lattice fcc 4.065 #lattice fcc 4.065 origin 0 0 0 orient x 1 1 0 orient y -1 1 0 orient z 0 0 1 spacing 1 1 1 #lattice fcc 4.065 origin 0 0 0 orient x 1 -1 0 orient y 1 1 -2 orient z 1 1 1 spacing 1 1 1 #region mybox block -1 1 -1 1 -1 1 region mybox cylinder z 0 0 7 -28 28 create_box 1 mybox create_atoms 1 region mybox # potentials pair_style eam/alloy pair_coeff * * Au.lammps.eam Au neighbor 1.0 bin neigh_modify delay 1 check yes timestep 0.001 variable N equal step variable pote equal pe variable T equal temp variable V equal vol variable Press equal press variable Etotal equal etotal variable H equal etotal+press*vol*0.625E-6 variable vol equal volume # define computes #compute myTemp all temp #compute disp all displace/atom compute pot all pe/atom compute ke all ke/atom compute s all stress/atom virial compute s1 all stress/atom compute p all reduce sum c_s1[1] c_s1[2] c_s1[3] variable Press equal -(c_p[1]+c_p[2]+c_p[3])/(3*vol) thermo_style custom step temp etotal ke pe lz ly lx press pxx pyy pzz dump 21 all custom 1000 AuAu*.txt id type x y z c_pot c_ke v_Press min_style sd minimize 1.0e-6 1.0e-6 1000 1000 # initialvelocities velocity all create 273 887723 dist gaussian units box thermo 1000 fix 1 all nve fix 2 all temp/berendsen 273 273 1.0 run 20000 dump 5 all custom 1000 AuAu*.txt id type x y z c_pot c_ke v_Press unfix 1 unfix 2 |
回复此楼» 猜你喜欢
PRB投稿，Acknowledgment sent to author状态十天了
已经有3人回复
-
PVA溶解
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有285人回复
-
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
2楼2016-06-09 08:29:54
guozhijing
银虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 254.1
- 帖子: 44
- 在线: 36.8小时
- 虫号: 4435776
- 注册: 2016-02-25
- 性别: MM
- 专业: 高分子材料结构与性能
3楼2016-09-18 19:21:24
投票:
10