| 查看: 4023 | 回复: 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步 |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有12人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 【求助】LAMMPS-运行报错
已经有9人回复
largehoo
木虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 2164.5
- 红花: 2
- 帖子: 150
- 在线: 207.5小时
- 虫号: 2855707
- 注册: 2013-12-07
- 专业: 爆炸与冲击动力学
7楼2014-07-12 15:40:05
2楼2013-05-21 21:06:40
3楼2013-05-22 13:47:03
4楼2013-07-17 09:30:28
