4.3.3 计算结果分析

首先施加自重应力场，命令输入如下：

ef moduli_data

s_mod=y_mod/(2.0*(1.0+p_ratio))

b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))

end

set y_mod 1.6e12

set p_ratio 0.2

moduli_data

mo el

prop bulk b_mod shear s_mod

fix x range x -100.1 -99.9

fix x range x 149.9 150.1

fix z range z -60.1 -59.9

fix y range y -0.1 0.1

fix y range y 1.9 2.1

ini density 2000

set gravity 0 0 -9.8

set large

;apply szz-0.045e6 range z 11.9 12.1

;释放0

solve

save地震地应力.sav

加载地震荷载

;围岩(5级)

def moduli_data11

s_mod=y_mod/(2.0*(1.0+p_ratio))

b_mod=y_mod/(3.0*(1.0-2.0*p_ratio))

end

set y_mod le9

set p_ratio 0.45

moduli_data11

;围岩(2级)

def moduli_data111

s_mod1=y_mod1/(2.0*(1.0+p_ratio1))

b_mod1=y_mod1/(3.0*(1.0-2.0*p_ratio1))

end

set y_mod1 3.3e10

set p_ratio1 0.2

moduli_data111

mo mohr

prop bulk b_mod1 shear s_mod1 fric 6 ......     （共4983字）    [阅读本文]>>