24小时热门版块排行榜    

查看: 2886  |  回复: 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),
我不知道问题出在哪儿了,希望高手给予指点,谢谢!
回复此楼
在这个世界上,你必须成为你希望看到的改变!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhoudeyuan: 金币+5, ★★★很有帮助 2013-03-22 09:22:42
zhoudeyuan: 金币+2, ★★★很有帮助 2013-03-23 16:15:02
第一个问题手册里有说明。 第二个是对asperity的原子重复移动了, integrater组理包含了asperity, 所以fix        1 integrater nve
对asperity的原子多了一次位移。
2楼2013-03-22 09:02:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoudeyuan

银虫 (小有名气)

引用回帖:
2楼: Originally posted by ddll_slsl at 2013-03-22 09:02:09
第一个问题手册里有说明。 第二个是对asperity的原子重复移动了, integrater组理包含了asperity, 所以fix        1 integrater nve
对asperity的原子多了一次位移。

嗯,谢谢你!
在这个世界上,你必须成为你希望看到的改变!
3楼2013-03-22 09:21:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ddll_slsl

木虫 (正式写手)

【答案】应助回帖

第一个问题, 手册里说明了, 部分的thermo参数无法准确表示整个热压。 把对象设为all就可以了。
4楼2013-03-22 11:29:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yangwenbin99

铜虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
也不一定有错误吧,我昨天也出现差不多同样的问题,但是换了一个版本的执行程序就没报警错误了
5楼2013-03-22 14:20:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoudeyuan

银虫 (小有名气)

引用回帖:
4楼: Originally posted by ddll_slsl at 2013-03-22 11:29:18
第一个问题, 手册里说明了, 部分的thermo参数无法准确表示整个热压。 把对象设为all就可以了。

回复的迟了,谢谢你!

» 本帖已获得的红花(最新10朵)

在这个世界上,你必须成为你希望看到的改变!
6楼2013-03-23 16:13:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhoudeyuan

银虫 (小有名气)

引用回帖:
5楼: Originally posted by yangwenbin99 at 2013-03-22 14:20:16
也不一定有错误吧,我昨天也出现差不多同样的问题,但是换了一个版本的执行程序就没报警错误了

我是在win7下串行计算的,那我在换个机子试试,谢谢你!
在这个世界上,你必须成为你希望看到的改变!
7楼2013-03-23 16:14:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

闵磊子

银虫 (小有名气)

送红花一朵
请问如何把你这里的文件保存成lammps例子那样的格式啊,我直接复制,不会弄
8楼2015-04-01 23:01:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

PANGHS

木虫 (初入文坛)

送红花一朵
引用回帖:
6楼: Originally posted by zhoudeyuan at 2013-03-23 16:13:27
回复的迟了,谢谢你!...

请问楼主 是哪一个参数改为all 我也在做压痕
9楼2016-04-16 09:06:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhoudeyuan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 gy1nBQXYQJqL 2026-08-29 7/350 2026-08-30 06:46 by ZPa0EcMwuECS
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 G6APbkg8SA6w 2026-08-29 5/250 2026-08-30 02:40 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:45 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 01:34 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3 kynobel 2026-08-23 3/150 2026-08-29 21:57 by 余韵清
[基金申请] 面上意见出来了 +8 黄鸟于飞Chao 2026-08-29 16/800 2026-08-29 20:45 by 黄鸟于飞Chao
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +9 苏东坡二世 2026-08-23 9/450 2026-08-29 14:41 by hustersqt
[基金申请] 为什么资助数各大高校都创新高,自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了!? +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见