24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 3647  |  回复: 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的回帖

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的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

优秀版主文献杰出贡献文献杰出贡献

【答案】应助回帖

fix   2 all npt temp 300.0 300.0 100.0 aniso 0.0 0.0 x 10.0 drag 1.0
x 为非周期,所以报错误
专业文献求助
5楼2013-07-19 17:41:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

天泊乐

金虫 (正式写手)

【答案】应助回帖

fix npt 命令不能用于非周期性边界,如要使用,涉及压强参数要全为NULL,也就是控不了压了
音乐you我!
6楼2013-08-07 11:39:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

largehoo

木虫 (小有名气)

【答案】应助回帖

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

chinaboy798

银虫 (小有名气)

【答案】应助回帖

这个in文件有几点问题,个人建议
1.非周期条件不能使用npt
2.为什么你要固定了Z轴两端的速度做能量最小化,而后在拉伸的时候是对所有原子施加npt系综做弛豫?
3.fix nvt已经有一个温度控制了,又加一个temp/rescale?
4.既然最后使用的fix deform来施行拉伸,之前施加Z轴两端速度为零的能量最小化有何意义?
8楼2014-11-11 14:45:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaohu199316

新虫 (初入文坛)

引用回帖:
1楼: Originally posted by 江流儿0305 at 2013-05-21 20:45:16
以下是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-- ...

请问,Al拉伸in文件还有吗?

发自小木虫Android客户端
9楼2017-04-14 13:08:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 江流儿0305 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +6 chenxrlkx 2026-04-05 6/300 2026-04-05 22:05 by imissbao
[考研] 求调剂求调剂 +8 121. 2026-04-02 8/400 2026-04-05 20:15 by lys0704
[考研] 296求调剂 +3 汪!?! 2026-04-05 4/200 2026-04-05 20:13 by 啵啵啵0119
[考研] 308求调剂 +3 终不似从前 2026-04-05 3/150 2026-04-05 20:07 by 啵啵啵0119
[考研] 329求调剂 +17 miaodesi 2026-04-02 20/1000 2026-04-05 18:33 by 蓝云思雨
[考研] 0860 求调剂 一志愿国科大 348 分 +3 WiiiP 2026-04-03 3/150 2026-04-05 17:43 by Ecowxq666!
[考研] 266分,一志愿电气工程,本科材料,求材料专业调剂 +8 哇呼哼呼哼 2026-04-02 9/450 2026-04-05 17:14 by lbsjt
[考研] 298分 070300求调剂 +15 zwen03 2026-04-02 15/750 2026-04-05 12:52 by Hdyxbekcb
[考研] 290求调剂085701 +10 1314捧花 2026-04-02 10/500 2026-04-05 10:19 by Sealedwind
[考研] 一志愿安徽某211 0703化学总分339求调剂 +6 晚风不晚 2026-04-04 6/300 2026-04-04 20:11 by dongzh2009
[考研] 333求调剂 +9 阿科逸 2026-03-31 9/450 2026-04-04 18:25 by macy2011
[考研] 277工科求调剂 +7 1915668 2026-04-04 7/350 2026-04-04 17:21 by 啊俊!
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 五邑大学土木工程招调剂生2026 +3 wyutj 2026-03-31 4/200 2026-04-03 18:21 by zengxj_7201
[考研] 交通运输考试264分求工科调剂 +4 jike777 2026-04-02 4/200 2026-04-02 21:53 by zllcz
[考研] 求调剂 302分初试 0854 +5 伶可乐 2026-04-02 5/250 2026-04-02 17:53 by 笔落锦州
[考研] 085900土木水利336分求调剂 +4 Zhangjiangj 2026-03-31 6/300 2026-04-02 11:40 by 1753564080
[考研] 372求调剂 +3 jj涌77 2026-04-02 3/150 2026-04-02 09:57 by olim
[考研] 302求调剂一志愿北航070300,本科郑大化学 +8 圣日耳曼条 2026-04-01 11/550 2026-04-02 07:40 by chemdavid
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
信息提示
请填处理意见