24小时热门版块排行榜

返回列表

张牧童

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 3392.8
帖子: 559
在线: 79.3小时
虫号: 2306377
注册: 2013-02-28
专业: 燃烧学

[求助] 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

回复此楼

» 猜你喜欢

天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有6人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助，我安装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人回复

1楼2013-05-13 20:31:34

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖支持 ( 显示支持度最高的前 50 名 )

王者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
我的源程序！