24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3720  |  回复: 7

朱雪刚521

铁杆木虫 (知名作家)

[求助] lammps运行出错:ERROR: Illegal fix nvt/npt/nph command (fix_nh.cpp:317)

整个结果如下。
LAMMPS (29 Jun 2012)
# 3d Lennard-Jones melt

units                metal
boundary        p p p

atom_style        atomic

variable        x equal 2.5                 #temperatue frome 2.5K

lattice                hcp 2.268                #short distance of Be hcp
Lattice spacing in x,y,z = 2.268 3.92829 3.70363

region                box block 0 10 0 10 0 10 #0 22.8560 0 22.8560 035.8320 frome MS creat 2000 atome

create_box        1 box
Created orthogonal box = (0 0 0) to (22.68 39.2829 37.0363)
  1 by 1 by 1 MPI processor grid

create_atoms         1 box
Created 4000 atoms

timestep        0.01                        #why others have 0.005

thermo                1000                        #why

#velocity        all create 3.0 87287

pair_style        tersoff
pair_coeff        * * Be.tersoff Be

neighbor        0.5 bin                        #why
neigh_modify        every 20 delay 0 check no #yes or no

#velocity        all create $x 825577 dist gaussian

fix                NVT all nvt $x $x 1.0 drag 0.2
fix                NVT all nvt 2.5 $x 1.0 drag 0.2
fix                NVT all nvt 2.5 2.5 1.0 drag 0.2
ERROR: Illegal fix nvt/npt/nph command (fix_nh.cpp:317)

加个reax?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

朱雪刚521

铁杆木虫 (知名作家)

跟编辑语言有关么?还是我的参数设置的不对
我的fix
fix                1 all nvt $x $x 1.0 drag 0.2

下边还有个fix暂时还没执行的,会不会也是出毛病呢
fix                 1 all npt $x 2000 4.00 xyz 0.0 0.0 6.0 drag 0.2
2楼2014-05-05 17:08:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
朱雪刚521: 金币+30, ★★★★★最佳答案, 谢谢,问题解决,但是我二楼列出的问题,表现出来了,加上temp后还是不行,加了couple也不行? 2014-05-06 14:49:53
月只蓝: 金币+3, 感谢指导! 2014-05-16 20:04:45
你要仔细看Manual,看看命令的格式,在fix nvt t命令当中,在nvt 字样之后,要有一个关键字的,你没有写这个关键字。fix npt也是一样。

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

3楼2014-05-06 09:43:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)

很早之前,大概两三年之前吧,你这个命令形式是对的,但是这个命令格式早就修改过了。
4楼2014-05-06 09:44:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wangshuai_xx

银虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
朱雪刚521: 金币+15, 谢谢你,问题解决了,但是在二楼的第二个命令没解决了,加了temp,不行,再在xyz前加上couple还是不行。 2014-05-06 14:51:39
月只蓝: 金币+2, 鼓励交流 2014-05-16 20:05:05
manual
Syntax:
fix ID group-ID style_name keyword value ...

你没有写keyword
temp or iso or aniso or tri or x or y or z or xy or yz or xz or couple or tchain or temp
Es muss sein
5楼2014-05-06 12:54:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

刘雪鹏

木虫 (小有名气)

少了keyword temp
6楼2014-05-07 09:28:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

朱雪刚521

铁杆木虫 (知名作家)

送红花一朵
引用回帖:
3楼: Originally posted by 老虎大王 at 2014-05-06 09:43:26
你要仔细看Manual,看看命令的格式,在fix nvt t命令当中,在nvt 字样之后,要有一个关键字的,你没有写这个关键字。fix npt也是一样。

再次求助:fix  和 msd  rdf有关,如下贴。
http://muchong.com/bbs/viewthread.php?tid=7413859
7楼2014-05-15 21:55:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tiancaisha

金虫 (正式写手)


引用回帖:
2楼: Originally posted by 朱雪刚521 at 2014-05-05 17:08:00
跟编辑语言有关么?还是我的参数设置的不对
我的fix
fix                1 all nvt $x $x 1.0 drag 0.2

下边还有个fix暂时还没执行的,会不会也是出毛病呢
fix                 1 all npt $x 2000 4.00 xyz 0.0 0.0 6.0 drag 0.2

亲 你这个是弛豫的意思吗?
8楼2015-12-25 14:05:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 朱雪刚521 的主题更新
信息提示
请填处理意见