我写了一个分层土边坡的FLAC命令流,可是计算不出结果,运行了十几个小时还没结果,现在把它发出来,请各位大侠帮我看看毛病,谢谢了:) 4 -tC=>>wc
kaCN^yQ
(O+d6oT=Z2
new
$L= Dky7
;*****几何模型***** `w
6Qsah
impgrid Flacmodel_haitang.Flac3D \B4H0f
plot block group aPelt`
plot add axes red ym{@w3"S
set gravity 0,0,-9.81 Hq}g1?b
group soil range group 1 any nz?[
group rock range group 2 any o08g]a
;施加边界约束条件 A6=Z2i0w>X
fix x y z rang y -0.1 0.1 DKt98;
fix z range x 0 212 z -0.1 0.1 IVh5SS
fix x range x -0.1 0.1 z 0 130 #$BFTlm|
fix z range x 0 212 z 129.9 130.1 O z0-cM8t
fix x range x 211.9 212.1 z 0 130 Acw`ytV
;通过弹性求解生成初始应力场 u"qu!EY2
model elas xZ {6!=4!
pro bulk 1e10 she 1e10 range group soil w317]-n
pro bulk 1e10 she 1e10 range group rock i) e6U(H
solve fo 10
~frsgHW
;位移和速度归零 <CVX[R]U
ini xdisp 0 ydisp 0 zdisp 0 }3: mn
ini xvel 0 yvel 0 zvel 0 h7+"*fN
;定义材料特性 +=J$:/&U
model mohr x 4v:67_^
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6 & uNn1qV
dil 10 range group soil ysOf=~1
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6 & 29J|eBvxx
dil 15 range group rock @HbRfD/!
hist n=5 SEI0G_wk$
hist unbal +
-Rf@
plo hist 1 +1T>Ob;hk
set mech ratio 1.0e-6 G{: B'08
solve