24小时热门版块排行榜    

查看: 4018  |  回复: 8
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

江流儿0305

新虫 (初入文坛)

[求助] lammps实验拉伸Al纳米线的in文件的错误提示好久找不出? 已有2人参与

以下是in文件一直提示错误ERROR: Illegal fix nvt/npt/nph command (fix_nh.cpp:317)
还有WARNING: Resetting reneighboring criteria during minimization (min.cpp:173)求高手。
# 3d EAM tensile simulation---------三维EAM拉伸仿真

#units         metal        #设定原子模型
dimension     3        #几何维度

boundary     s s p        #边界设定
atom_style     atomic

# create geometry---------------创建几何模型

lattice        fcc 0.4056 orient x 1 1 -2 orient y -1 1 0 orient z 1 1 1 origin 0.0 0.0 0.5  #晶格结构,晶格常数x,y,z坐标的晶向
region        box block 0 4 0 4 0 20    #模拟几何区域设定
create_box 1 box        #创建一个box
create_atoms 1 box         #在box中创建1一种原子
mass        1 26.981    #定义原子质量为26.981

# EAM potentials----------------势函数

pair_style    eam        #定义势函数
pair_coeff    * * Al_jnp.eam    #势所需要的参数在此文件里

neighbor    2.0 bin        #截断半径
neigh_modify    delay 1 check yes #邻近原子列表更新速度

# define groups----------------定义组

region         1 block INF INF INF INF INF 3         #定义一个叫1的区域
group        lower region 1                #此区域里的原子叫lower
region        2 block  INF INF INF INF 17 INF        #定义一个叫2的区域
group        upper region 2                #此区域里的原子叫upper
group        boundary union lower upper            #定义lower+upper = boundary
group        mobile subtract all boundary        #定义mobile=all- boundary

#initial velocities-------------设置原子初速度

velocity    lower set 0.0 0.0 0.0    #设置left区域原子初速度为0
velocity    upper set 0.0 0.0 0.3   #设置right区域z方向原子速度为0.3
velocity    mobile ramp vz 0.0 0.3 z 3 17 sum yes   

min_style sd                     #能量最小模式   
minimize 1.0e-6 1.0e-6 1000 10000        #驰豫

fix   1 lower setforce 0.0 NULL NULL     #固定下边的原子
fix   2 all npt temp 300.0 300.0 100.0 aniso 0.0 0.0 x 10.0 drag 1.0

thermo 200                 #每200步输出一次结果
timestep    0.003                       #时间步长0.003ps
run  5000                 #运行5000步

unfix2
fix   3 all nvt  
fix   4 all temp/rescale 300 1.0 1.0 5 1.0
fix   5 all deform 1000 z erate 0.002 units box  #设置沿z方向拉伸,应变率为0.002(1/ps)

dump 1 all atom 500 dump.tensile       #输出结果到dump.tensile
min_style sd                     #能量最小模式   
minimize 1.0e-6 1.0e-6 1000 10000        #驰豫
run  100000                            #运行10000步
回复此楼
倚天照海花无数,高山流水心自知。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

largehoo

木虫 (小有名气)

【答案】应助回帖

lattice 的单位是埃
7楼2014-07-12 15:40:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

jiaoyixiong

荣誉版主 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
ben_ladeng: 金币+1, 谢谢版主 2013-05-21 22:58:32
Illegal fix nvt/npt/nph command


错误提示很明显
2楼2013-05-21 21:06:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

江流儿0305

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by jiaoyixiong at 2013-05-21 21:06:40
Illegal fix nvt/npt/nph command


错误提示很明显

对啊,但就是找不出错误,求教!
倚天照海花无数,高山流水心自知。
3楼2013-05-22 13:47:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

韩欢欢

金虫 (初入文坛)

引用回帖:
3楼: Originally posted by 江流儿0305 at 2013-05-22 13:47:03
对啊,但就是找不出错误,求教!...

请问你的问题解决了吗?我也遇到了这个问题,不知道怎么解决。我在拉伸的时候用的是fix  1 nvt 300.0 300.0 100.0 drag 0.2 chain no 是想在300k恒温拉伸。运行也是这个错误提示,请教下怎么修改
Demi
4楼2013-07-17 09:30:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +10 家与远方 2026-08-10 15/750 2026-08-14 02:08 by 绵羊哥哥
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +9 Tide man 2026-08-10 10/500 2026-08-13 20:33 by 温淼淼淼
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] filecode +14 documentary 2026-08-10 16/800 2026-08-13 16:39 by 润生71
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 好奇怪的filecode +5 布布和一二 2026-08-08 6/300 2026-08-12 16:11 by 云上清扬
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
信息提示
请填处理意见