24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

江流儿0305

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 10.5
帖子: 4
在线: 1.6小时
虫号: 2364315
注册: 2013-03-20
性别: GG
专业: 计算机软件

[求助] 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步

回复此楼