udec初始应力计算不收敛，求助大家帮忙找下原因

qq846765019

以下是我的代码：
new
round 0.05
block 0,0 0,90.1 200,90.1 200,0
crack 0,30 200,30
crack 0,33.5 200,33.5
crack 0,58.5 200,58.5
crack 0,60.1 200,60.1

;1#
jregion id 1 0,0 0,30 200,30 200,0
jset 0,0 200,0 0,0 3,0 0,0 range jregion 1
jset 90,0 3,0 3,0 3,0 0,0 range jregion 1
jset 90,0 3,0 3,0 3,0 1.5,3 range jregion 1
;2#
jregion id 2 0,30 0,33.5 200,33.5 200,30
jset 0,0 200,0 0,0 1.75,0 0,30 range jregion 2
jset 90,0 1.75,0 1.75,0 1.75,0 0,30 range jregion 2
jset 90,0 1.75,0 1.75,0 1.75,0 0.8525,31.75 range jregion 2
;3#
jregion id 3 0,33.5 0,58.5 200,58.5 200,33.5
jset 0,0 200,0 0,0 2.5,0 0,33.5 range jregion 3
jset 90,0 2.5,0 2.5,0 2.5,0 0,33.5 range jregion 3
jset 90,0 2.5,0 2.5,0 2.5,0 1.25,36 range jregion 3
;4#
jregion id 4 0,58.5 0,60.1 200,60.1 200,58.5
jset 0,0 200,0 0,0 1.6,0 0,58.5 range jregion 4
jset 90,0 1.6,0 1.6,0 1.6,0 0,58.5 range jregion 4

;5#
jregion id 5 0,60.1 0,90.1 200,90.1 200,60.1
jset 0,0 200,0 0,0 3,0 0,60.1 range jregion 5
jset 90,0 3,0 3,0 3,0 0,60.1 range jregion 5
jset 90,0 3,0 3,0 3,0 1.5,63.1 range jregion 5

change cons=3
change mat=1 range jr 1
change mat=2 range jr 2
change mat=3 range jr 3
change mat=4 range jr 4
change mat=5 range jr 5

change jcons=2
change jmat=1 range jr 1
change jmat=2 range jr 2
change jmat=3 range jr 3
change jmat=4 range jr 4
change jmat=5 range jr 5

prop mat 1 dens=2566 bulk=0.26e11 shear=1.73e9 ten=2.5e6 coh=8.63e6 fric=35
prop mat 2 dens=1400 bulk=0.77e10 shear=0.86e9 ten=1.45e6 coh=0.86e6 fric=33
prop mat 3 dens=2561 bulk=1.87e10 shear=1.97e9 ten=4.43e6 coh=8.25e6 fric=33
prop mat 4 dens=1400 bulk=0.77e10 shear=0.86e9 ten=1.45e6 coh=0.86e6 fric=33
prop mat 5 dens=2571 bulk=1.82e10 shear=0.88e9 ten=5.66e6 coh=7.67e6 fric=32

prop jmat 1 jkn=2.2e9   jks=2.2e9   jfr=35  jcoh=2e6
prop jmat 2 jkn=1.2e9   jks=1.2e9   jfr=35  jcoh=3e6
prop jmat 3 jkn=1.35e9 jks=1.35e9 jfr=35  jcoh=2e6
prop jmat 4 jkn=1.2e9   jks=1.2e9   jfr=35  jcoh=3e6
prop jmat 5 jkn=3.2e9   jks=3.2e9   jfr=35  jcoh=3e6

set grav 0,-10
insitu stress -2.125e6,0,-2.125e7   szz -2.125e6
bound stress 0,0,-2.125e7 range        0,200        90.0 90.2
bound stress -2.125e6,0,0  range -0.1,0.1        0,90.1
bound stress -2.125e6,0,0  range 199.9,200.1  0,90.1
bound  xv 0  range    -0.1    0.1    0    90.1
bound  xv 0  range    199.9    200.1   0    90.1
bound  yv 0  range      0      200   -0.1   0.1
set ov 1
hist ydis 30 14.6
hist unbalance
so

ljm

好东西，谢谢楼主

我又不缺爱

ljm 发表于 2014-3-28 11:29
好东西，谢谢楼主

晕。。。你气死楼主了

长空皓月

初始应力方向