24小时热门版块排行榜    

查看: 2625  |  回复: 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朵)

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

在这个世界上,你必须成为你希望看到的改变!
已阅   回复此楼   关注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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 9/450 2026-03-17 10:29 by xujiaoszu
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +5 大火山小火山 2026-03-16 7/350 2026-03-17 10:26 by peike
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 304求调剂 +5 素年祭语 2026-03-15 5/250 2026-03-16 17:00 by 我的船我的海
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 招收0805(材料)调剂 +3 18595523086 2026-03-13 3/150 2026-03-14 00:33 by 123%、
[考研] 求调剂(材料与化工327) +4 爱吃香菜啦 2026-03-11 4/200 2026-03-13 22:11 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 301求调剂 +6 Liyouyumairs 2026-03-11 6/300 2026-03-13 20:11 by JourneyLucky
[考研] 310求调剂 +3 【上上签】 2026-03-11 3/150 2026-03-13 16:16 by JourneyLucky
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 308求调剂 +3 是Lupa啊 2026-03-12 3/150 2026-03-13 14:30 by 求调剂zz
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
信息提示
请填处理意见