| 查看: 588 | 回复: 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 |
回复此楼» 猜你喜欢
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有32人回复
-
物理学I论文润色/翻译怎么收费?
已经有139人回复
-
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有78人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
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
