24小时热门版块排行榜    

查看: 670  |  回复: 0

li_xue_song

新虫 (小有名气)

[求助] 为什么模拟的裂纹拉不开?

1,谢谢您来解决问题!
2,问题:建立了含裂纹铝模型,y方向上设置速度,但是拉伸的时候拉不开是什么原因?
3,in文件如下所示:
## fourth simulation ##
units metal
boundary p s p
dimension 3

lattice fcc 4.0503
region box block 0 20 0 20 0 2
create_box 1 box
create_atoms 1 box

neighbor 1.5 bin
neigh_modify every 2 delay 10

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

pair_style eam/fs
pair_coeff * * Al_mm.eam.fs Al

mass * 26.981

fix 1 all nve

velocity all create 300 12345 rot yes
#velocity all scale 300
reset_timestep 0
timestep 0.01
thermo 100
thermo_style custom step temp pe ke press
dump 1 all atom 50 dump.xue4_1
log log.xue4_1
run 500

reset_timestep 0

region crack block INF 5 9 11 INF INF
group crack region crack
delete_atoms group crack

#fix 2 lower setforce 0 0 0
velocity upper set 0 1 0
velocity lower set 0 -1 0
#velocity mobile ramp vy -1 1 y 2 18 sum yes

variable strain equal (ly-81)/81

thermo_style custom step temp pe ke ly pxx pyy pzz v_strain
dump 2 all atom 100 dump.xue4_2
log log.xue4_2
#if "$(temp)>1000" then quit

run 5000
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 li_xue_song 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 是这周出结果还是下周出结果? +3 yuleib84 2026-08-11 3/150 2026-08-11 21:53 by jnhyjjm
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 F4zcb7xRCEJA 2026-08-11 3/150 2026-08-11 21:51 by pMlf9MbdMHyc
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +17 Lanmanbaby 2026-08-09 27/1350 2026-08-11 21:12 by klxdzh
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 yPv6ndDtiRxQ 2026-08-11 3/150 2026-08-11 21:09 by thgx1vqMqrLw
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 yPv6ndDtiRxQ 2026-08-11 3/150 2026-08-11 20:51 by thgx1vqMqrLw
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 yPv6ndDtiRxQ 2026-08-10 7/350 2026-08-11 20:39 by thgx1vqMqrLw
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +3 majunge000 2026-08-11 4/200 2026-08-11 20:13 by lch2012
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 yPv6ndDtiRxQ 2026-08-10 8/400 2026-08-11 19:56 by thgx1vqMqrLw
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 GpvwGNVshCOp 2026-08-10 7/350 2026-08-11 19:39 by thgx1vqMqrLw
[基金申请] 听说今天filecode变了 +26 布布和一二 2026-08-06 49/2450 2026-08-11 13:18 by WH3796
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 17/850 2026-08-11 10:50 by 要乐观耀哥
[基金申请] 基金中了 +15 laoda193707 2026-08-06 15/750 2026-08-11 00:11 by jiafei2190
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见