| 查看: 1114 | 回复: 1 | ||
[求助]
使用read_data后，如何建立模型group？
|
|
小弟最近在学习lammps软件进行分子模拟，自己建了一个polycrystal模型，想对模型进行group，但是运行后报错，应该是region语句的格式问题，但是刚开始学习不知道该怎么改，请虫友们指点，附件上是in文件（in.4txt），data文件（fetwin.lmp）和势函数（用的是铝的势函数，随便找了一个，应该不是势函数的问题） 这是in文件 # ------------------------ INITIALIZATION ---------------------------- units metal dimension 3 boundary p p p atom_style atomic # variable latparam equal 4.05 # ----------------------- ATOM DEFINITION ---------------------------- # lattice fcc ${latparam} # region whole block 0 10 0 10 0 10 # create_box 1 whole # lattice fcc ${latparam} orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 # create_atoms 1 region whole read_data fetwin.data region up block -INF INF -INF INF 80 INF region low block -INF INF -INF INF -INF 20 group up region up group low region low # ------------------------ FORCE FIELDS ------------------------------ pair_style eam/alloy pair_coeff * * Al99.eam.alloy Al # ------------------------- SETTINGS --------------------------------- compute csym all centro/atom fcc compute peratom all pe/atom compute stress all stress/atom NULL ###################################### reset_timestep 0 timestep 0.001 velocity all create 300 12345 mom yes rot no fix 1 all npt temp 300 300 1 iso 0 0 1 drag 1 # Set thermo output thermo 10 thermo_style custom step lx ly lz press pxx pyy pzz pe temp run 500 unfix 1 # Store final cell length for strain calculations variable tmp equal "lx" variable L0 equal ${tmp} print "Initial Length, L0: ${L0}" 报错的内容：： PS C:\Users\xxpc-office\Desktop\tensile-twin> lmp_mpi -in in.4txt LAMMPS (11 Aug 2017) OMP_NUM_THREADS environment is not set. Defaulting to 1 thread. (../comm.cpp:90) using 1 OpenMP thread(s) per MPI task Reading data file ... orthogonal box = (0 0 0) to (100 100 100) 1 by 1 by 1 MPI processor grid reading atoms ... 84355 atoms ERROR: Expected floating point parameter in input script or data file (../region_block.cpp:38) Last command: region up block -INF INF -INF INF 80 INF data文件开头内容：： 84355 atoms 1 atom types 0.00000000 100.00000000 xlo xhi 0.00000000 100.00000000 ylo yhi 0.00000000 100.00000000 zlo zhi Masses 1 55.84500000 # Fe Atoms # atomic 1 1 25.11560200 83.72774500 76.47610300 2 1 28.27335500 85.32570400 75.71484200 3 1 22.03587500 80.44872200 74.03231600 4 1 21.11270900 82.67894400 74.88441700 5 1 17.95495500 81.08098500 75.64567900 6 1 23.57573800 82.08823300 75.25421000 7 1 24.27046200 84.27690300 74.12315500 8 1 28.96807900 87.51437300 74.58378700 9 1 30.58597000 87.47282000 72.60063100 10 1 27.42821600 85.87486200 73.36189300 11 1 25.81032600 85.91641400 75.34504900 12 1 32.12583300 89.11233200 73.82252500 13 1 32.82055700 91.30100100 72.69147000 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : Al99.eam.alloy
- 附件 2 : fetwin.lmp
2018-07-03 10:24:26, 762.16 K
2018-07-03 10:24:28, 5.55 M
» 猜你喜欢
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
为什么网上很多人说本周 12号出结果
已经有3人回复
-
静等基金结果
已经有16人回复
» 本主题相关商家推荐: (我也要在这里推广)
2楼2018-07-03 10:30:58
100