| 查看: 2895 | 回复: 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朵）
闵磊子» 猜你喜欢
面上合作单位盖章
已经有9人回复
-
面上函评意见出来了，像什么等级？
已经有17人回复
-
为什么国自然不能直接公布
已经有5人回复
-
梦想
已经有7人回复
-
学科评审组评审是指会评吗？
已经有4人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有14人回复
-
小白求助 投论文要求的highlights应该如何写
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有12人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 请问下玻璃做纳米压痕有什么用？
已经有12人回复
闵磊子
银虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 455.2
- 散金: 50
- 帖子: 146
- 在线: 98.5小时
- 虫号: 3459945
- 注册: 2014-10-08
- 专业: 机械摩擦学与表面技术
8楼2015-04-01 23:01:45
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
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
4楼2013-03-22 11:29:18
20
投票:
