24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

kite175

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 329.7
帖子: 22
在线: 33.4小时
虫号: 1886180

[交流] 有人做裂纹模拟的吗？看看我这个裂纹的in文件有什么问题

# 3d Morse crack simulation

units       metal
boundary s s p
atom_style atomic

# create geometry

lattice          fcc 4.05  orient x 1 0 0  orient y 0 1 0 orient z 0 0 1
region box block 0 60 0 40 -0.25 0.25
create_box 5 box
create_atoms 1 box

mass 1 1.0
mass 2 1.0
mass 3 1.0
mass 4 1.0
mass 5 1.0

# Morse potentials

pair_style morse 2.5
pair_coeff * * 1.0 1.0 2.5
neighbor 0.3 bin
neigh_modify every 5 delay 0 check yes

# define groups

region          1 block INF INF INF 2 INF INF
group lower region 1
region 2 block INF INF 38 INF INF INF
group upper region 2
group boundary union lower upper
group mobile subtract all boundary

region leftupper block INF 20 20 INF INF INF
region leftlower block INF 20 INF 20 INF INF

region       bonds block  INF 20  19.75 20.25 -0.25 0.25
delete_atoms region  bonds

group leftupper region leftupper
group leftlower region leftlower

group halfupper1 region halfupper1
group halfupper2 region halfupper2

set group leftupper type 3
set group leftlower type 4

set group lower type 1
set group upper type 2

neigh_modify exclude type 3 4
neigh_modify exclude type 1 2

# 能量最小化
min_style sd
minimize    1.0e-12 1.0e-12  10000 10000

# initial velocities 弛豫

compute       new mobile temp
velocity mobile create 0.001 887723 temp new
timestep 0.001
thermo 200
run 10000

velocity upper  set 0.0 0.01 0.0
velocity lower  set 0.0 -0.01 0.0
velocity mobile ramp vy 0.01 0.01 y 20 38 sum yes

# fixes
fix          1 all nve
fix          2 lower setforce  0.0 NULL 0.0
fix          3 upper setforce  0.0 NULL 0.0
# run
thermo 200
thermo_modify temp new

dump 1 all atom 1000 dump.crack.lammpstrj
run 300000

回复此楼

» 猜你喜欢

AI 太可怕了，写基金时，提出想法，直接生成的文字比自己想得深远，还有科学性已经有6人回复
有院领导为了换新车，用横向课题经费买了俩车已经有9人回复
酰胺脱乙酰基已经有13人回复
博士延得我，科研能力直往上蹿已经有8人回复
同年申请2项不同项目，第1个项目里不写第2个项目的信息，可以吗已经有4人回复
有时候真觉得大城市人没有县城人甚至个体户幸福已经有10人回复
天津大学招2026.09的博士生，欢迎大家推荐交流（博导是本人）已经有5人回复
遇见不省心的家人很难过已经有22人回复

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2012-08-27 09:17:48

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

tfl03

木虫 (职业作家)

应助: 24 (小学生)
金币: 3128.5
帖子: 4453
在线: 1812.2小时
虫号: 772313

★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zh1987hs: 金币+3, 谢谢 2012-10-01 11:02:11

# 3d Morse crack simulation

units       metal
boundary       s s p
atom_style       atomic

# create geometry

lattice             fcc 4.05  orient x 1 0 0  orient y 0 1 0 orient z 0 0 1
region             box block 0 60 0 40 -0.25 0.25
create_box       5 box
create_atoms       1 box

mass             1 1.0
mass             2 1.0
mass             3 1.0
mass             4 1.0
mass             5 1.0

# Morse potentials

pair_style       morse 2.5
pair_coeff       * * 1.0 1.0 2.5
neighbor       0.3 bin
neigh_modify       every 5 delay 0 check yes

# define groups

region             1 block INF INF INF 2 INF INF
group             lower region 1
region             2 block INF INF 38 INF INF INF
group             upper region 2
group             boundary union lower upper
group             mobile subtract all boundary

region             leftupper block INF 20 20 INF INF INF
region             leftlower block INF 20 INF 20 INF INF

region       bonds block  INF 20  19.75 20.25 -0.25 0.25
delete_atoms region  bonds

group             leftupper region leftupper
group             leftlower region leftlower

group             halfupper1 region halfupper1
group             halfupper2 region halfupper2

set             group leftupper type 3
set             group leftlower type 4

set             group lower type 1
set             group upper type 2

neigh_modify       exclude type 3 4
neigh_modify       exclude type 1 2

# 能量最小化
min_style sd
minimize    1.0e-12 1.0e-12  10000 10000

# initial velocities 弛豫

compute                new mobile temp
velocity       mobile create 0.001 887723 temp new
timestep       0.001
thermo             200
run 10000

velocity       upper  set 0.0 0.01 0.0
velocity       lower  set 0.0 -0.01 0.0
velocity       mobile ramp vy 0.01 0.01 y 20 38 sum yes

# fixes
fix          1 all nve
fix          2 lower setforce  0.0 NULL 0.0
fix          3 upper setforce  0.0 NULL 0.0
# run
thermo             200
thermo_modify       temp new

dump             1 all atom 1000 dump.crack.lammpstrj
run             300000