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

回复此楼

» 猜你喜欢

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

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

» 抢金币啦！回帖就可以得到:
查看全部散金贴

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

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

tfl03

木虫 (职业作家)

应助: 24 (小学生)
金币: 3168.5
帖子: 4452
在线: 1812小时
虫号: 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