我写了一个分层土边坡的FLAC命令流,可是计算不出结果,运行了十几个小时还没结果,现在把它发出来,请各位大侠帮我看看毛病,谢谢了:) q;*'V9#
Kg`P@
S|af?IW
new
En)Ptz#0
;*****几何模型***** hV:++g
impgrid Flacmodel_haitang.Flac3D %J Jp/I
plot block group &XCP@@T
plot add axes red P+BGCc%);B
set gravity 0,0,-9.81 uQ|LkL%<^
group soil range group 1 any Ve${g`7&
group rock range group 2 any Y>atJ
;施加边界约束条件 P^!g0K
fix x y z rang y -0.1 0.1 JI
cm$
fix z range x 0 212 z -0.1 0.1 O(VWJ@EHn
fix x range x -0.1 0.1 z 0 130 cQ'x]u_
fix z range x 0 212 z 129.9 130.1 Y% JE})
fix x range x 211.9 212.1 z 0 130 %|D)U>o{
;通过弹性求解生成初始应力场 =CO) Q2
model elas vKWi?}1
pro bulk 1e10 she 1e10 range group soil ^4NRmlb
pro bulk 1e10 she 1e10 range group rock vxOnv8(
solve fo 10 \GQRpJ#h1
;位移和速度归零 R !yh0y}Z
ini xdisp 0 ydisp 0 zdisp 0 UBJYs{zz
ini xvel 0 yvel 0 zvel 0 {~w( pAx
;定义材料特性 P&=YLL<W
model mohr &*Xrh7K2e
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6 & OQ4Pk/-'
dil 10 range group soil QO0@Ax\b
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6 & y5F"JjQAa
dil 15 range group rock ?zQW9e
hist n=5 v 2p
hist unbal !3iZa*
plo hist 1 #d*0
)w
set mech ratio 1.0e-6 crOSr/I$
solve