24小时热门版块排行榜    

查看: 538  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今年也是没消息就是没中么 +5 vito刘 2026-06-13 9/450 2026-06-14 11:31 by newfuzzy1
[考研] *O_o^ +3 jk31rb6vyz 2026-06-12 3/150 2026-06-14 09:08 by 8bmv3h4hvh
[找工作] *囍 +4 m584r2z4j3 2026-06-12 4/200 2026-06-14 09:06 by 8bmv3h4hvh
[论文投稿] ⊙▽⊙ +3 m584r2z4j3 2026-06-11 3/150 2026-06-14 08:46 by 8bmv3h4hvh
[硕博家园] T^T +7 m584r2z4j3 2026-06-11 8/400 2026-06-14 08:34 by 8bmv3h4hvh
[考博] →_→ +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-14 08:28 by 8bmv3h4hvh
[考博] T^T +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-14 08:26 by 8bmv3h4hvh
[硕博家园] O(∩_∩)O +6 jk31rb6vyz 2026-06-12 8/400 2026-06-14 08:10 by 8bmv3h4hvh
[考博] T^T +4 m584r2z4j3 2026-06-11 4/200 2026-06-14 08:08 by 8bmv3h4hvh
[公派出国] π_π +5 m584r2z4j3 2026-06-11 5/250 2026-06-14 08:08 by 8bmv3h4hvh
[基金申请] 函评 +3 布布和一二 2026-06-13 3/150 2026-06-13 23:35 by 芋圆香飘飘
[考研] T^T +6 m584r2z4j3 2026-06-11 6/300 2026-06-13 18:29 by lkmi35jdzt
[文学芳草园] 咸菜 +3 myrtle 2026-06-10 3/150 2026-06-13 18:29 by ahsoarli
[考博] π_π +4 m584r2z4j3 2026-06-11 4/200 2026-06-13 18:25 by lkmi35jdzt
[找工作] 售SCI文章,我:8.O.5.51O.5.4,科目齐全,可十急 +4 q1zg6x4lf3 2026-06-12 4/200 2026-06-13 08:53 by t4hz8376wy
[论文投稿] →_→ +5 m584r2z4j3 2026-06-11 5/250 2026-06-13 08:13 by t4hz8376wy
[考博] [急招]吉林大学招收26年9月入学考核制博士生1名,截止今晚24点前 +3 hehejin 2026-06-10 4/200 2026-06-12 18:24 by 冻雨天
[教师之家] AASSDDFF +3 m584r2z4j3 2026-06-11 3/150 2026-06-12 11:28 by zs48m4lz8g
[考博] 还有课题组有博士名额吗 +5 冻雨天 2026-06-07 6/300 2026-06-08 20:59 by RCS_tyro
[考博] 博士申请 +3 张城斐 2026-06-07 3/150 2026-06-08 09:15 by A大朋同学
信息提示
请填处理意见