| 查看: 6844 | 回复: 13 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
shiling1813新虫 (小有名气)
|
[求助]
lammps中的fix npt命令用后压强不为零的原因？已有1人参与
|
|
|
大家好！ 我最近在使用lammps计算金属的剪切模量，想在xy方向施加剪切应变，产生xy方向的剪切应力，且其它五个方向（x、y、z、xz、zy方向）的应力为零，从而计算xy方向的剪切模量。控压命令使用的是fix npt命令。 为什么计算出来的压强或应力在其它五个方向（x、y、z、xz、zy方向）不为零？ 我的in文件和应力计算结果如下： units metal boundary p p p atom_style atomic lattice fcc 1 read_data noreallrelax000lamyb.dat change_box all triclinic pair_style eam/alloy pair_coeff * * nialre_djp.eam.alloy Ni Al Re velocity all create 0.5 49280 rot yes dist gaussian thermo 1000 timestep 0.001 #fix 1 all nvt temp 0.5 0.5 0.1 fix 1 all npt temp 0.5 0.5 0.1 x 0.0 0.0 1 y 0.0 0.0 1 z 0.0 0.0 1 yz 0.0 0.0 1 xz 0.0 0.0 1 xy 0.0 0.0 1 couple none drag 1.0 compute spa all stress/atom NULL compute spb all stress/atom thermo_temp compute spc all stress/atom NULL virial compute ppptem zhongjian reduce sum c_spb[1] c_spb[2] c_spb[3] c_spb[4] c_spb[5] c_spb[6] compute ppp zhongjian reduce sum c_spc[1] c_spc[2] c_spc[3] c_spc[4] c_spc[5] c_spc[6] fix avt zhongjian ave/time 1 1000 10000 c_ppp fix avttem zhongjian ave/time 1 1000 10000 c_ppptem run 10000 unfix 1 fix 2 all npt temp 0.5 0.5 0.1 x 0.0 0.0 1 y 0.0 0.0 1 z 0.0 0.0 1 yz 0.0 0.0 1 xz 0.0 0.0 1 couple none drag 1.0 label loopa variable u loop 50 change_box all xy delta 0.1219364 units box run 10000 next u jump in.lamstressnpt loopa tmp1.rdf文件中的应力计算结果如下： # Time-averaged data for fix avt # TimeStep Number-of-rows # Row c_ppp 10000 6 1 -23328.7 2 -108325 3 -23041.5 4 -1066.45 5 7104.27 6 -11874 20000 6 1 -424350 2 35272.4 3 458826 4 -107510 5 -3254.15 6 5063.74 30000 6 1 568412 2 201561 3 -427249 4 2.29794e+06 5 -873.69 6 1600.17 40000 6 1 357465 2 258743 3 -154082 4 2.55284e+06 5 -2511.4 6 3019.35 50000 6 1 141415 2 221818 3 -8094.18 4 2.99189e+06 5 -5015.37 6 8539.99 。 。 。 |
回复此楼» 猜你喜欢
读博
已经有5人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有13人回复
-
博士申请都是内定的吗？
已经有6人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有5人回复
-
博士读完未来一定会好吗
已经有29人回复
-
投稿精细化工
已经有4人回复
-
高职单位投计算机相关的北核或SCI四区期刊推荐，求支招！
已经有4人回复
-
导师想让我从独立一作变成了共一第一
已经有9人回复
-
心脉受损
已经有5人回复
-
Springer期刊投稿求助
已经有4人回复
zhaijianhui
铁杆木虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 5793.8
- 散金: 1100
- 红花: 2
- 帖子: 953
- 在线: 663.1小时
- 虫号: 1030672
- 注册: 2010-05-28
- 专业: 金属材料的力学行为
13楼2017-09-06 06:29:35
2楼2017-08-10 09:50:59
shiling1813
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 38
- 红花: 1
- 帖子: 128
- 在线: 88.2小时
- 虫号: 725356
- 注册: 2009-03-18
- 专业: 金属材料的微观结构
送红花一朵|
谢谢！ 可是对时间统计平均也没看到结果应力控制为零啊？是我统计的原子数太少了吗？ 我从体系中选取了120个原子，统计了1000个时间步得到的应力结果如下： # Time-averaged data for fix avt # TimeStep Number-of-rows # Row c_ppp 10000 6 1 -23328.7 2 -108325 3 -23041.5 4 -1066.45 5 7104.27 6 -11874 20000 6 1 -424350 2 35272.4 3 458826 4 -107510 5 -3254.15 6 5063.74 30000 6 1 568412 2 201561 3 -427249 4 2.29794e+06 5 -873.69 6 1600.17 40000 6 1 357465 2 258743 3 -154082 4 2.55284e+06 5 -2511.4 6 3019.35 50000 6 1 141415 2 221818 3 -8094.18 4 2.99189e+06 5 -5015.37 6 8539.99 。 。 。 |
3楼2017-08-10 13:24:19
4楼2017-08-10 21:56:35