| 查看: 3696 | 回复: 11 | ||
[求助]
in文件运行失败（lost atoms）求助！！
|
||
|
刚刚接触lammps没多久，写个小的in文件试试手，运行中丢失原子，无法输出想要的结果~求指导！最好请大神运行下，给点建议！谢谢~ #3d contact simulation units metal dimension 3 boundary p p s atom_style atomic neighbor 3.0 bin neigh_modify delay 5 # create geometry region box block 0 80 0 80 0 100 units box create_box 3 box mass 1 64 mass 2 12 #atom region region boundary block INF INF INF INF INF 35 units box region coating block INF INF INF INF 35 40 units box region tool cone z 40 40 0 18 40 60 units box lattice fcc 3.6149 create_atoms 1 region boundary units box lattice diamond 3.57 create_atoms 2 region coating units box lattice diamond 3.57 create_atoms 3 region tool units box # Eam Potentials pair_style eam pair_coeff * * Ni_u3.eam # define groups group boundary region boundary group mobile subtract all boundary group tool region tool set group boundary type 1 set group mobile type 2 #initial velocities compute new mobile temp/partial 0 0 1 velocity mobile create 0.1 482748 temp new velocity tool set 0 0 -5 sum yes # fixes fix 1 all nve fix 2 boundary setforce 0 0 0 fix 3 mobile temp/rescale 100 0.1 0.1 0.02 1.0 fix_modify 3 temp new # Run timestep 0.01 thermo 100 thermo_modify temp new dump 1 all atom 500 dump.yahen.lammpstrj run 2000 |
回复此楼» 猜你喜欢
» 本主题相关价值贴推荐，对您同样有帮助:
求助，我安装auto CAD2007时，提示：错误 1402。无法打开主键：怎么解决呀！急~~~~
已经有8人回复
- 求助 molpro运行出错！急
已经有5人回复
- 求助！win7下运行vc6.0时出现msdev错误怎么办？？？
已经有1人回复
- 【求助】CASTEP计算出错，运行失败，请高手帮忙。
已经有11人回复
- 【求助】此程序没有错误，怎么运行不出来？
已经有6人回复
- 【求助】运行出错 虚拟内存
已经有10人回复
- 【求助】文件不能运行，这是什么问题呢
已经有10人回复
- 求助：电脑待机后无法启动
已经有11人回复
- 【求助】Fortran程序，如何让exe文件可以单独运行？
已经有6人回复
- 【求助】xshell的exe文件为什么不能运行？
已经有3人回复
- 【求助】程序运行后出现的错误 access violation
已经有15人回复
- 【求助】安装问题，安装后程序无法执行
已经有12人回复
- 【求助】安装ansys时，运行Keygen程序,但不会自动生成一个文件License.dat,为什么?
已经有12人回复
- 【求助】请问会不会出现在两台电脑上运行同一个模拟文件出现结果不同的情况？
已经有6人回复
- 【求助】vasp运行问题
已经有7人回复
- 【求助】运行错误，大家看看问题在哪？
已经有7人回复
- 【讨论】求助word 2007文件无法打开
已经有2人回复
王者111111
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 11039.4
- 散金: 1200
- 红花: 2
- 帖子: 1552
- 在线: 162.7小时
- 虫号: 2395717
- 注册: 2013-04-02
- 性别: GG
- 专业: 传热传质学
【答案】应助回帖
|
dimension 3 boundary s s p newton on units metal atom_style atomic neighbor 2.0 bin neigh_modify delay 5 lattice fcc 3.52 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 region box1 block 0 30 0 15 -3.0 3.0 side in units lattice region 1 sphere 15.0 20.0 0.0 2.0 side in units lattice region box2 block 0 30 15 22 -3.0 3.0 units lattice region box3 block 0 0.5 15 22 -3.0 3.0 region box4 block 29.5 30 15 22 -3.0 3.0 region box union 5 box1 1 box2 box3 box4 units lattice create_box 3 box mass 1 58.71 mass 2 58.71 mass 3 12 region 2 block INF INF INF 0.5 INF INF units lattice region 3 block INF 0.5 INF INF INF INF units lattice region 4 block 29.5 INF INF INF INF INF units lattice create_atoms 1 region box1 units lattice create_atoms 2 region 2 units lattice create_atoms 2 region 3 units lattice create_atoms 2 region 4 units lattice delete_atoms region box3 delete_atoms region box4 lattice none lattice diamond 3.57 orient x 1 0 0 orient y 0 1 0 orient z 0 0 1 create_atoms 3 region 1 units lattice pair_style eam pair_coeff * * Ni_u3.eam # define groups group sphere region 1 group lower region 2 group leftboundary region 3 group rightboundary region 4 group boundary union lower rightboundary leftboundary group mobile subtract all boundary set group mobile type 1 set group lower type 2 set group leftboundary type 2 set group rightboundary type 2 set group sphere type 3 # initial velocities compute new mobile temp velocity mobile create 0.1E-4 482748 temp new fix 1 all nve fix 2 boundary setforce 0.0 0.0 0.0 fix 3 all temp/rescale 100 0.1E-4 0.1E-4 0.1E-7 1.0 # run with indenter timestep 0.001 variable y equal "0.8*ylat - step*dt*1.6*ylat" fix 4 sphere move variable NULL v_y NULL NULL NULL NULL thermo 100 thermo_modify temp new compute stress mobile stress/atom pair compute displ mobile displace/atom compute centro mobile centro/atom fcc compute pe mobile pe/atom dump 1 all custom 100 dump.*.trajectory id type xs ys zs c_stress[1] c_stress[2] c_stress[3] c_centro c_pe c_displ[1] c_displ[2] c_displ[3] c_displ[4] run 3000 # run without indenter variable dy equal "step*dt*1.6*ylat-20.0" fix 6 sphere move variable NULL v_dy NULL NULL NULL NULL run 1000 我的源程序！ |
9楼2013-05-15 23:15:45
ddll_slsl
木虫 (正式写手)
- 应助: 101 (高中生)
- 金币: 2939.2
- 红花: 6
- 帖子: 670
- 在线: 255.3小时
- 虫号: 2121784
- 注册: 2012-11-12
- 专业: 动力学与控制
2楼2013-05-13 20:47:39
王者111111
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 11039.4
- 散金: 1200
- 红花: 2
- 帖子: 1552
- 在线: 162.7小时
- 虫号: 2395717
- 注册: 2013-04-02
- 性别: GG
- 专业: 传热传质学
3楼2013-05-14 09:22:51
4楼2013-05-14 12:46:50
5楼2013-05-14 12:47:25
王者111111
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 11039.4
- 散金: 1200
- 红花: 2
- 帖子: 1552
- 在线: 162.7小时
- 虫号: 2395717
- 注册: 2013-04-02
- 性别: GG
- 专业: 传热传质学
6楼2013-05-14 14:14:13
xmcyeyanping
金虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 687.6
- 散金: 7
- 红花: 1
- 帖子: 38
- 在线: 49小时
- 虫号: 1555886
- 注册: 2011-12-29
- 专业: 金属材料的微观结构
7楼2013-05-15 16:40:10
8楼2013-05-15 21:20:09
王者111111
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 11039.4
- 散金: 1200
- 红花: 2
- 帖子: 1552
- 在线: 162.7小时
- 虫号: 2395717
- 注册: 2013-04-02
- 性别: GG
- 专业: 传热传质学
10楼2013-05-15 23:16:27
