抛砖引玉,欢迎大家讨论。
new
; w.gang.89@foxmail.com
gen zone brick size 10 10 10 edge 0.1
def para
; concrete C35
e_1 = 31e9
p_1 = 0.2
c_1 = 3.9e6
f_1 = 61.8
t_1 = 3.0e6
d_1 = 2500
s_1 = e_1 / (2.0 * (1.0 + p_1 ))
b_1 = e_1 / (3.0 * (1.0 - p_1 * 2.0))
end
@para
model mohr
property bulk @b_1 shear @s_1 cohesion @c_1 friction @f_1 tension @t_1
fix z range union z 0 z 0.1
initial zvel 1e-8 range z 0
initial zvel -1e-8 range z 0.1
history id 1 gp zdisp 0.05 0.05 0.1
history id 2 zone szz 0.05 0.05 0.1
history id 3 zone szz 0.05 0.0 0.05
plot clear
plot add zone colorby state
plot add zcon szz
plot add con zdisp
plot add history -3 vs -1
plot show
step 6000