gen zone radcyl p0 0 0 0 p1 7 0 0 p2 0 60 0 p3 0 0 7 dim 3 3 3 3 size 5 30 14 12 rat 1 1 1 1.2 group 围岩
gen zone cshell p0 0 0 0 p1 3 0 0 p2 0 60 0 p3 0 0 3 size 5 30 14 12 dim 2.7 2.7 2.7 2.7 group 衬砌 fill group 原土
gen zone reflect normal 0 0 1 ori 0 0 0
gen zone reflect normal 1 0 0 ori 0 0 0
gen zone brick p0 7 0 -7 p1 23 0 -7 p2 7 60 -7 p3 7 0 7 size 16 30 14 group 围岩
gen zone brick p0 -7 0 -23 p1 23 0 -23 p2 -7 60 -23 p3 -7 0 -7 size 30 30 16 group 围岩
gen zone brick p0 -7 0 7 p1 23 0 7 p2 -7 60 7 p3 -7 0 7.28 size 30 30 1 group 围岩
gen zone reflect normal 1 0 0 ori -7 0 0
gen zone brick p0 -37 0 7.28 p1 23 0 7.28 p2 -37 60 7.28 p3 -37 0 10.4 size 60 30 2 group 粘土
gen zone brick p0 -37 0 10.4 p1 23 0 10.4 p2 -37 60 10.4 p3 -37 0 14.4 size 60 30 2 group 素填土
;设置边界条件
fix x range x -36.99 -37.01
fix x range x 22.99 23.01
fix y range y -0.01 0.01
fix y range y 59.99 60.01
fix x range z -22.99 -23.01
fix y range z -22.99 -23.01
;设置重力
set grav 0 0 -9.81
;本构模型的确立以及材料的赋值
mod mohr
ini den 1.94e3 range group 素填土
ini den 1.91e3 range group 粘土
ini den 2.74e3 range group 围岩
ini den 2.74e3 range group 衬砌
ini den 2.74e3 range group 原土
pro bulk 6.67e6 shear 3.1e6 coh 1e4 fri 28 range group 素填土
pro bulk 1e7 shear 6e6 coh 6.7e4 fri 17.9 range group 粘土
pro bulk 3.2e8 shear 2.02e8 coh 3.26e7 fri 39.45 range group 围岩
pro bulk 3.2e8 shear 2.02e8 coh 3.26e7 fri 39.45 range group 衬砌
pro bulk 3.2e8 shear 2.02e8 coh 3.26e7 fri 39.45 range group 原土
hist unbal
step 1000
