24小时热门版块排行榜    

查看: 614  |  回复: 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 的主题更新
信息提示
请填处理意见