| 查看: 6856 | 回复: 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 。 。 。 |
回复此楼» 猜你喜欢
投稿Elsevier的Neoplasia杂志，到最后选publishing options时页面空白，不能完成投稿
已经有22人回复
-
申请26博士
已经有5人回复
-
职称评审没过，求安慰
已经有22人回复
-
垃圾破二本职称评审标准
已经有15人回复
-
EST投稿状态问题
已经有7人回复
-
毕业后当辅导员了，天天各种学生超烦
已经有4人回复
-
聘U V热熔胶研究人员
已经有10人回复
-
求助文献
已经有3人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有8人回复
-
三无产品还有机会吗
已经有6人回复
shiling1813
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 38
- 红花: 1
- 帖子: 128
- 在线: 88.2小时
- 虫号: 725356
- 注册: 2009-03-18
- 专业: 金属材料的微观结构
|
谢谢！基本解决了我的问题。看来我需要看看统计基础理论。 另外，再问个问题：从我的计算结果看，npt命令控压后，其它五个方向（包括动能项的）应力的输出结果可以接近零；而（不包括动能项的）应力的输出结果不能接近零。这是为什么？ （包括动能项的）应力的输出结果如下： # Time-averaged data for fix avttem # TimeStep Number-of-rows # Row c_ppptem 10000 6 1 70.5359 2 68.5396 3 71.3667 4 1.67524 5 0.703051 6 1.42924 20000 6 1 20.9834 2 66.5369 3 67.3381 4 4.24172e+07 5 -0.684267 6 0.38628 30000 6 1 27.1811 2 75.2281 3 67.9913 4 8.53567e+07 5 -0.400147 6 0.873144 （不包括动能项的）应力的输出结果： # Time-averaged data for fix avt # TimeStep Number-of-rows # Row c_ppp 10000 6 1 497893 2 500057 3 493104 4 3157.48 5 -2380.42 6 389.928 20000 6 1 947148 2 282209 3 263951 4 4.24208e+07 5 -1603.51 6 414.441 30000 6 1 1.14449e+06 2 183360 3 154346 4 8.53597e+07 5 3022.89 6 -75.536 。 。 。 |
5楼2017-08-11 13:35:21
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