最近做计算遇到了几个问题,望大神解答。1、用切向位移边界条件不收敛;2、用应力边界条件做收敛速度特别慢
不知道其中的原因,是因为网格划分问题还是边界条件设置问题
new
ge zone brick p0=(0,0,0) p1=(600,0,0) p2=(0,600,0) p3=(0,0,500) size 60 60 50
plot
add su yellow
add axes black
show
model mohr
pro bu=8.93e9 sh=6147540984 fric=45 co=1.2e6 ten=9.49e6
set gravity 0,0,-9.81
in de=2600
fix x range x 599.9 600.1
fix y range y -0.9 1.1
fix y range y 599.9 600.1
fix y z range z -0.9 1.1
a sxx=-3e6 range x -0.9 1.1
set mec fo 5
so
位移边界条件如下:
set gravity 0 0 0
free x y z
fix y range y -0.1 0.1
fix y range y 599.9 600.1
fix x y z range z -0.1 0.1
free y range x -0.1 0.1 y -0.1 0.1 z -0.1 600
free y range x -0.1 0.1 y 599.9 600.1 z -0.1 600
free y range x 599.9 600.1 y -0.1 0.1 z -0.1 600
free y range x 599.9 600.1 y 599.9 600.1 z -0.1 600
free y range x -0.1 0.1 y -0.1 600.1 z -0.1 0.1
free y range x 599.9 600.1 y -0.1 600.1 z -0.1 0.1
a yvel=1E-5 range x -0.1 0.1
a yvel=1E-5 range x 599.9 600.1
step 100
set mec fo 5
so
print zone s range group dian1