;开挖
new
;建模
gen zone brick &
p0 -8 0 -8 p1 16 0 -8 p2 -8 1 -8 p3 0 0 0 p4 16 1 -8 p5 0 1 0 p6 16 0 0 p7 16 1 0 &
size 8 1 4
pl sur
pause key
gen zone brick &
p0 4 0 0 p1 16 0 0 p2 4 1 0 p3 4 0 8 &
size 6 1 4
pause key
;开挖部分
gen zone brick &
p0 0 0 6 p1 4 0 6 p2 0 1 6 p3 0 0 8 &
size 2 1 1 group 1
pause key
gen zone brick &
p0 0 0 4 p1 4 0 4 p2 0 1 4 p3 0 0 6 &
size 2 1 1 group 2
pause key
gen zone brick &
p0 0 0 2 p1 4 0 2 p2 0 1 2 p3 0 0 4 &
size 2 1 1 group 3
pause key
p0 0 0 0 p1 4 0 0 p2 0 1 0 p3 0 0 2 &
size 2 1 1 group 4
pause key
;设置参数
; 定义模型材料参数
def modelMaterial
; 弹性模量(Pa)
E1=10e6
E2=20e6
E3=40e6
;E4=200e6
;
; 泊松比
poi1=0.3
poi2=0.3
poi3=0.3
;poi4=0.3
;
;粘结力(Pa)
coh1=20e3
coh2=22e3
coh3=32e3
;coh4=200e3
; 容重(N/m3)
weight1=16e3
weight2=17e3
weight3=19e3
;weight4=22e3
; 内摩擦角(度)
fri1=18
fri2=22
fri3=24
;fri4=30
;
;抗拉强度
ten1=4.3e6
ten2=5.645e6
grav0=-9.80 ; 重力加速度(N/kg)
dens1=-weight1/grav0 ; 密度(kg/m3)
dens2=-weight2/grav0
dens3=-weight3/grav0
;dens4=-weight4/grav0
K1=E1/(3*(1-2*poi1)) ; 体积模量(Pa)
K2=E2/(3*(1-2*poi2))
K3=E3/(3*(1-2*poi3))
;K4=E4/(3*(1-2*poi4))
G1=E1/(2*(1+poi1)) ; 剪切模量(Pa)
G2=E2/(2*(1+poi2))
G3=E3/(2*(1+poi3))
; G4=E4/(2*(1+poi4))
end
modelMaterial
model mohr
pro bulk K1 she G1 dens dens1 coh coh1 &
friction fri1 dilation 0. range z 6 8
pro bulk K2 she G2 dens dens2 coh coh2 &
friction fri2 dilation 0. range z 0 6
pro bulk K3 she G3 dens dens3 coh coh3 &
friction fri3 dilation 0. range z -8 0
fix z range z 7.9 8.1
fix x range x -7.1 -6.9
fix x range x 15.9 16.1
fix y
set gravity 0 0 -10
solve
save initialm.sav
ini xdis 0 ydis 0 zdis 0
model null range group 1
solve
model null range group 2
solve
model null range group 3
solve
model null range group 4
solve
save excavation.sav
solve fos file excfosm.sav associated
这个里面有错误还望指出来
