24小时热门版块排行榜

返回列表

影月beatrice

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
帖子: 3
在线: 6.9小时
虫号: 3032055
注册: 2014-03-06
专业: 工程热力学

[求助] 关于fix shake以及compute命令的应用问题已有1人参与

我最近写了一个关于超临界水的in文件，里面用到了fix shake命令，但运行的时候一直报错ERROR: Invalid fix style，我对照manuel检查了命令格式，没有发现问题，禁了fix shake后运行又出现ERROR: Invalid compute style，都没查出问题在哪里，请各位大侠不吝赐教，小女子先谢过了
以下是我的in文件
units real
dimension 3
boundary p p p
atom_style full
pair_style lj/charmm/coul/charmm 7.5 8.0 8.0 10.0
bond_style harmonic
angle_style harmonic
read_data data.Mol
angle_coeff 1  55.0000 104.5200
bond_coeff  1  450.0000 0.9572
pair_coeff 1 1 0.1521 3.1507
pair_coeff 1 2 0.0836 1.7753
pair_coeff 2 2 0.0460 0.4000
region 1 block INF INF INF INF INF INF
velocity all create 640 100
velocity all scale 640
#group water type 1 2
fix 1 all shake 0.0001 20 0 b 1 a 1
fix 2 all npt temp 640 760 1000 iso 236.92 236.92 1000 drag 1
timestep 0.01
compute 1 all ke
compute 2 all pe
compute 3 all temp
compute 4 all dilatation/atom
compute 5 all displace/atom
thermo 500
thermo_style custom temp press density vol pe ke etotal enthalpy
run 1000

data文件如下
LAMMPS data file from msi2lmp v3.8 for Mol

   6 atoms
   4 bonds
   2 angles
   0 dihedrals
   0 impropers

2 atom types
1 bond types
1 angle types

-0.854007846    3.056692154 xlo xhi
-0.764937233    3.145762767 ylo yhi
-1.347689527    2.563010473 zlo zhi

Masses

1  15.999400
2 1.007970

Atoms

   1    1 1 -0.798200    0.121798883    0.387575451    0.043193656 0 1 0
   2    1 2  0.399100    0.308427222    1.054689176 -0.656528781 0 1 0
   3    1 2  0.399100    0.163119246 -0.433256524 -0.481969457 0 1 0
   4    2 1 -0.798200    1.973583684    2.001768480    1.988670964 0  -1 0
   5    2 2  0.399100    2.453643042    1.441265665    1.343893247 0  -1 0
   6    2 2  0.399100    1.587480847    2.690434354    1.408703210 0  -1 0

Bonds

   1 1    1    2
   2 1    1    3
   3 1    4    5
   4 1    4    6

Angles

   1 1    2    1    3
   2 1    5    4    6

回复此楼