我做的是地下洞室,想分析下塑性区的变形,但是塑性区显示的特别的小,用的是弹性求解计算的地应力,然后PWIPP模型计算的,洞室里面气压加6MPa,下面是程序,跪求高手指点!!!
set creep on
model mohr ran group 2 any group 6 any group 3 any group 4 any
prop bulk 14.667e9 shear 8.8e9 co 5.0e6 fri 30 ten 1.0e6 ran group 2 any group 6 any group 3 any group 4 any
model pwipp ran group 1 any group 5 any
prop bulk 26.67e9 shear 5.71e9 D_wipp 5.79e-36 act 50241.6 gas 8.3192 &
a_wipp 4.56 b_wipp 127 e_dot 5.39e-8 n_wipp 5.0 temp 300 &
tension 1e6 qvol=0.5 qdil=0.5 kshear=3.5e6 ran salt
set creep off
solv
def soften ; strain-softening law
while_stepping
soft_rep = soft_rep + 1
if soft_rep >= 10
soft_rep = 0
p_z = zone_head
if z_group(p_z) = 'salt'
loop while p_z # null
eplas = z_prop(p_z,'es_plastic')
if eplas # 0.0
rat_fac = 1.0 - eplas / 2e-4
if rat_fac > 0.0
z_prop(p_z,'kshear') = max_strength * rat_fac
else
z_prop(p_z,'kshear') = 0.0
endif
endif
p_z = z_next(p_z)
endLoop
endif
endif
end
set max_strength = 3.5, soft_rep=0
ini xdis 0 ydis 0 zdis 0
set creep on
m n range group 5 any group 6 any
apply nstress -6e6 range sphere center 0 0 0 rad 69 x -56.91 0 y 0 56.91 z -69 -39.02
apply nstress -6e6 range cyl end1 0 0 -37.02 end2 0 0 -39.02 rad 58.29
apply nstress -6e6 range cyl end1 0 0 -26.01 end2 0 0 -37.02 rad 62.86
apply nstress -6e6 range cyl end1 0 0 -15.01 end2 0 0 -26.01 rad 63.95
apply nstress -6e6 range cyl end1 0 0 -4 end2 0 0 -15.01 rad 63.95
apply nstress -6e6 range cyl end1 0 0 0 end2 0 0 -4 rad 63.18
apply nstress -6e6 range cyl end1 0 0 10.10 end2 0 0 0 rad 62.70
apply nstress -6e6 range cyl end1 0 0 20.02 end2 0 0 10.10 rad 60.93
apply nstress -6e6 range cyl end1 0 0 21.20 end2 0 0 20.20 rad 58.42
apply nstress -6e6 range cyl end1 0 0 33.31 end2 0 0 21.20 rad 58.13
apply nstress -6e6 range cyl end1 0 0 45.42 end2 0 0 33.31 rad 53.59
apply nstress -6e6 range cyl end1 0 0 57.53 end2 0 0 45.42 rad 47.81
apply nstress -6e6 range cyl end1 0 0 69.63 end2 0 0 57.53 rad 40.85
apply nstress -6e6 range cyl end1 0 0 71.13 end2 0 0 69.63 rad 32.78
apply nstress -6e6 range sphere center 0 0 57.53 rad 33.47 x -31.56 0 y 0 31.56 z 71.13 91
set large
set creep mindt 1e-16
set creep maxdt 1e3
set creep dt auto on
set creep lfob 5e-5
set creep ufob 1e-4
solv age 63.072e6
save try10_2.sav