哪位大神帮忙给看看,为什么我的模型中修正剑桥模型单元中无材料属性?万分感谢。能帮忙解决问题的给打赏。一下是模型部分;
config fluid
set fluid biot on
;第一步:建模
gen zone brick p0 0 0 0 p1 50 0 0 p2 0 50 0 p3 0 0 20 &
p4 50 50 0 p5 0 50 20 p6 50 0 20 p7 50 50 20 &
size 25 25 10 group 1
gen zone brick p0 0 0 20 p1 50 0 20 p2 0 50 20 p3 0 0 50 &
p4 50 50 20 p5 0 50 50 p6 50 0 50 p7 50 50 50 &
size 25 25 15 group 2
;第二步:设置边界条件
fix z ran z 0
fix x ran x 0
fix x ran x 50
fix y ran y 0
fix y ran y 50
;第三步:初始地应力计算
model mohr ;设置摩尔库伦本构模型
prop bulk 3.9e8 shear 8e7 coh 10e3 fri 35 ten 0 range group 2
prop bulk 3.9e8 shear 8e7 coh 25e3 fri 10 ten 0 range group 1
ini dens 1500;设置密度
model fl_iso
ini fden 1000 ften -5e5
ini biot_mod 4e9 range group 2
ini biot_mod 4e9 range group 1
prop perm 10e-7 biot_c 1 range group 2
prop perm 10e-13 biot_c 0.6 range group 1
ini szz -0.9e6 grad 0 0 20e3 ran z 0 30
ini syy -0.45e6 grad 0 0 10e3 ran z 0 30
ini sxx -0.45e6 grad 0 0 10e3 ran z 0 30
ini szz -3e5 grad 0 0 15e3 ran z 30 50
ini syy -1.5e5 grad 0 0 7.5e3 ran z 30 50
ini sxx -1.5e5 grad 0 0 7.5e3 ran z 30 50
ini pp 0.3e6 grad 0 0 -10e3 ran z 0 30
set grav 0 0 -10
fix pp range z -0.1 0.1
plot add cont pp out on ;在当前视图中绘出孔隙压力等值线阴影图
plot add flow plane ;在当前剖平面图内绘出流体流动向量
plot add cont szz ;绘出zz方向应力云图
plot add cont zdis ;绘出zz方向位移云图
plot show
set mech force 1e-3
solve ;按软件默认精度求解(自重作用下)初始地应力场
save 01.sav
;第四步:位移、速度清零
ini xvel 0 yvel 0 zvel 0
ini xdisp 0 ydisp 0 zdisp 0
config fluid
set fluid biot on
;第五步:施加恒荷载
app nstress -200e3 ran z 49.9 50.1 x 0 50 y 0 50
;第六步:流体模型
model fl_iso ;各项同性流体
prop perm 1e-13 range group 1 ;设置渗透系数
prop perm 1e-7 range group 2
;prop poro 0.7 ;设置孔隙率
;ini fmod 2.0e9 ; 初始化水的体积模量
ini ftens -1e-3 ; 初始化水的抗拉强度
ini fdensity 1e3 ; 初始化 水密度
;ini sat 0.7 ;饱和度
prop biot_c 0.8 ;比奥系数
;ini biot_mod 0.8e9 ;比奥模量
app pp 0.3e6 grad 0 0 -1.0e4 range z 0 30 x -0.1 0.1 ;透水边界
app pp 0.3e6 grad 0 0 -1.0e4 range z 0 30 x 49.9 50.1 ;透水边界
app pp 0.3e6 grad 0 0 -1.0e4 range z 0 30 y -0.1 0.1 ;透水边界
app pp 0.3e6 grad 0 0 -1.0e4 range z 0 30 y 49.9 50.1 ;透水边界
app pp 0.3e6 range z -0.1 0.1 ;透水边界
model cam-clay
prop shear 15e5 bulk_bound 20e6 range group 1
prop mm 0.85 lambda 0.054 kappa 0.0234 mpc 0.4e6 mp1 1.0 mv_l 3.32 range group 1
ini dens 1500 range group 1
;mohr-coulomb model 本构模型
model mohr ;设置摩尔库伦本构模型
prop bulk 35e6 shear 15e6 coh 10e3 fri 35 ten 0 range group 2
ini dens 2000 range group 2
def camclay_ini_p
pnt=zone_head
loop while pnt # null
OCR=1.0
s1=abs(z_sxx(pnt))
s2=abs(z_syy(pnt))
s3=abs(z_szz(pnt))
p0=(s1+s2+s3)/3.0-z_pp(pnt)
z_prop(pnt,'cam_cp')=p0
q0=sqrt(((s1-s2)^2+(s2-s3)^2+(s3-s1)^2)*0.5)
temp1=q0/(z_prop(pnt,'mm')*p0)
pc=p0*(1.0+temp1^2)*OCR
z_prop(pnt,'mpc')=pc
pnt=z_next(pnt)
endloop
end
camclay_ini_p
set gravity 0,0,-10
plot add cont pp out on ;在当前视图中绘出孔隙压力等值线阴影图
plot add flow plane ;在当前剖平面图内绘出流体流动向量
plot add cont szz ;绘出zz方向应力云图
plot add cont zdis ;绘出zz方向位移云图
plot show
hist id=1 gp zdis 0,0,0
hist id=2 gp zdis 25,25,25
hist id=3 gp zdis 0,0,50
hist id=4 gp zdis 50,0,50
hist id=5 gp zdis 50,50,50
hist id=6 gp zdis 25,25,50
hist id=7 gp pp 25 25 0
hist fltime
hist unbal
solve
save 02.sav