| 查看: 4116 | 回复: 8 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
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步 |
回复此楼» 猜你喜欢
申博发邮件
已经有5人回复
-
各位大神，目前国内有哪些比较好用的逆合成软件？
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
-
西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后
已经有5人回复
-
国社科系统bug了，是不是要放榜了？
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】LAMMPS-运行报错
已经有9人回复
4楼2013-07-17 09:30:28
2楼2013-05-21 21:06:40
3楼2013-05-22 13:47:03
w76990480
荣誉版主 (文坛精英)
super-emuch
- 应助: 47 (小学生)
- 贵宾: 13.729
- 金币: 109765
- 散金: 29386
- 红花: 337
- 沙发: 50
- 帖子: 23156
- 在线: 3166.8小时
- 虫号: 942733
- 注册: 2010-01-15
- 专业: 固体力学
- 管辖: 文献求助
5楼2013-07-19 17:41:05
