| 查看: 2599 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
zhoudeyuan银虫 (小有名气)
|
[求助]
纳米压痕过程中程序出错，哪位大神帮我解答下！
|
|
|
我的纳米压痕程序如下，运行结束后，牛顿原子层很不稳定，且log文件提示错误，不知道哪位大神能够帮我解答下，小弟万分感谢！ # 3d contact simulation(Sphere_Slid) units metal dimension 3 boundary p p p atom_style atomic neighbor 3.0 bin neigh_modify delay 5 newton on # create geometry region box block 0.0 40.0 0.0 40.0 0.0 50.0 units box create_box 4 box mass 1 29.0 mass 2 29.0 mass 3 29.0 mass 4 12.0 # atom regions region reg block 0.0 40.0 0.0 40.0 0.0 25.0 units box region reg1 block 3.0 37.0 3.0 37.0 3.0 25.0 units box region newtonian_layer block 8.0 32.0 8.0 32.0 8.0 25.0 units box region asperity sphere 20 20 30 5 units box lattice diamond 5.431 create_atoms 1 region reg units box lattice diamond 3.5667 create_atoms 2 region asperity units box # Tersoff Potentials pair_style tersoff pair_coeff * * SiC.tersoff Si Si Si C # define groups group reg region reg group reg1 region reg1 group newtonian_layer region newtonian_layer group thermostat_layer subtract reg1 newtonian_layer group boundary_layer_down subtract reg reg1 group asperity region asperity group integrater union reg asperity set group newtonian_layer type 1 set group thermostat_layer type 2 set group boundary_layer_down type 3 set group asperity type 4 # initial velocities compute newTemp thermostat_layer temp velocity thermostat_layer create 293.0 5812775 temp newTemp units box # fixed fix 1 integrater nve fix 2 boundary_layer_down setforce 0.0 0.0 0.0 velocity newtonian_layer set NULL NULL NULL units box velocity thermostat_layer create 293.0 5812775 temp newTemp units box fix 3 thermostat_layer temp/rescale 10 293.0 293.0 10 1 fix_modify 3 temp newTemp fix 4 asperity setforce 0.0 0.0 0.0 compute forcez newtonian_layer property/atom fz compute force newtonian_layer reduce sum c_forcez compute forcex newtonian_layer property/atom fx compute lateralforce newtonian_layer reduce sum c_forcex timestep 0.001 thermo 5 thermo_style custom step temp press pe ke etotal c_force c_lateralforce vol thermo_modify temp newTemp dump 1 all atom 50 dump.contact_block_slid fix 6 asperity move linear 0.0 0.0 -1.0 units box run 2000 fix 6 asperity move linear 0.0 0.0 1.0 units box run 2000 log文件出错： WARNING: Temperature for thermo pressure is not for group all (thermo.cpp:432) WARNING: One or more atoms are time integrated more than once (modify.cpp:216)， 我不知道问题出在哪儿了，希望高手给予指点，谢谢！ |
回复此楼» 本帖已获得的红花（最新10朵）
闵磊子» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有4人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问下玻璃做纳米压痕有什么用？
已经有12人回复
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
4楼2013-03-22 11:29:18
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
2楼2013-03-22 09:02:09
zhoudeyuan
银虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 485
- 红花: 2
- 帖子: 78
- 在线: 61.8小时
- 虫号: 2324272
- 注册: 2013-03-06
- 性别: GG
- 专业: 动力学与控制
3楼2013-03-22 09:21:32
yangwenbin99
铜虫 (小有名气)
- 应助: 15 (小学生)
- 金币: 545.9
- 散金: 215
- 红花: 10
- 帖子: 284
- 在线: 111.8小时
- 虫号: 2058623
- 注册: 2012-10-13
- 专业: 机械摩擦学与表面技术
5楼2013-03-22 14:20:16
10