24小时热门版块排行榜    

Znn3bq.jpeg
查看: 525  |  回复: 1

527751273

银虫 (正式写手)

[求助] 用LJ模拟裂纹扩展,可是结果出来都是从加力的地方断。。。。求解决啊。。。。

我想用LJ模拟裂纹扩展,可是结果出来都是从加力的地方断,而试验是在约束的地方开始断的,实在不懂啊。。

回复此楼

» 猜你喜欢

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

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

527751273

银虫 (正式写手)

in文件如下:
# 2d LJ crack simulation
dimension 2
boundary s s p
atom_style atomic
neighbor 0.3 bin
neigh_modify delay 3
# create geometry
lattice hex 0.93
region box block 0 10 0 40 -5 5
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
# LJ potentials
pair_style  lj/cut 2.5
pair_coeff  * * 1.0 1.0 2.5
# define groups
region 1 block 9.99 10 0 0.01 INF INF
group load1 region 1
region 2 block 9.99 10 23.99 24 INF INF
group load2 region 2
group load union load1 load2
group mobile subtract all load
region 3 block 0 0.01 15.99 16 INF INF
region 4 block 0 0.01 39.99 40 INF INF
group lower1 region 3
group lower2 region 4
group boundary union lower1 lower2
region  5 block INF 2 19.99 20.01 INF INF
region  6 block 8 INF 19.99 20.01 INF INF
group  lowernotch region 5
group  uppernotch region 6
group  notch union lowernotch uppernotch
delete_atoms  group notch
set  group load1 type 2
set  group load2 type 3
set  group lower1 type 4
set  group lower2 type 5
# initial velocities
velocity   load set -2.5e-3 0.0 0.0
minimize   1.0e-5 1.0e6 1000 1000
# fixes
fix  1 all nve
fix  2 load1 setforce 1000 NULL 0.0
fix  3 load2 setforce 5000 NULL 0.0
fix  4 boundary setforce 0.0 0.0 0.0
# run
timestep  0.00003
thermo  200
dump 1  all atom 500 dump8.lammpstrj
run  6000
2楼2011-08-24 21:24:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 527751273 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿085502,267分求调剂 +14 再忙也要吃饭啊 2026-04-08 15/750 2026-04-12 22:33 by fqwang
[教师之家] 山东双非院校考核超级无底线,领导幸灾乐祸,教师遭殃恐 +3 qut2026 2026-04-11 7/350 2026-04-12 20:24 by qut2026
[考研] 化学工程调剂289 +44 yang婷 2026-04-07 50/2500 2026-04-12 02:36 by 秋豆菜芽
[考研] 271求调剂 +20 2261744733 2026-04-11 22/1100 2026-04-11 23:14 by labixiaoqiao
[考研] 求调剂 +11 月@163.com 2026-04-07 13/650 2026-04-11 22:55 by BruceLiu320
[考研] 调剂 +6 青灯不负 2026-04-09 6/300 2026-04-11 20:35 by dongdian1
[考研] 085410 273分调剂 +4 X1999 2026-04-09 4/200 2026-04-11 13:05 by pies112
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 材料与化工调剂 10+11 下一站上岸@ 2026-04-10 36/1800 2026-04-11 10:26 by 89436494
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 284求调剂 +12 archer.. 2026-04-10 13/650 2026-04-11 08:44 by zhq0425
[材料工程] 材料调剂推荐 +8 蛋糕x2 2026-04-07 8/400 2026-04-10 23:13 by Ftglcn90
[考研] 085800 能源动力求调剂 +6 阿biu啊啊啊啊啊 2026-04-10 6/300 2026-04-10 15:03 by hemengdong
[考研] 一志愿211 0703化学 346分求调剂 +22 土豆er? 2026-04-09 23/1150 2026-04-10 10:58 by 高维春
[考研] 材料复试求调剂 +20 xhhdjdjsjks 2026-04-09 20/1000 2026-04-10 10:25 by 孙小小12457
[考研] 274求调剂 +5 山阿蔓 2026-04-07 5/250 2026-04-09 15:28 by 18828373951
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 22408 调剂材料 +7 我叫ez 2026-04-06 8/400 2026-04-07 17:12 by 蓝云思雨
[考研] 307求调剂 +3 所念及所望 2026-04-06 3/150 2026-04-06 17:30 by 土木硕士招生
信息提示
请填处理意见