我写了一个分层土边坡的FLAC命令流,可是计算不出结果,运行了十几个小时还没结果,现在把它发出来,请各位大侠帮我看看毛病,谢谢了:) 02tt.0go
ohbU~R3{U
@ta:9wZ
new :%z#s
;*****几何模型***** Lk!m1J5
impgrid Flacmodel_haitang.Flac3D \FUMfo^
plot block group ]
s^7c
plot add axes red <(@Z#%O9)
set gravity 0,0,-9.81 i\_LLXc
group soil range group 1 any suzK)rJ9i
group rock range group 2 any kia[d984w
;施加边界约束条件 gD51N()s,
fix x y z rang y -0.1 0.1 R[14scV
fix z range x 0 212 z -0.1 0.1 H~TuQ
fix x range x -0.1 0.1 z 0 130 pSb tm74
fix z range x 0 212 z 129.9 130.1 'pT13RFD
fix x range x 211.9 212.1 z 0 130 b*(K;`9)B
;通过弹性求解生成初始应力场 Yj'9|4%+|
model elas 2WDe34
pro bulk 1e10 she 1e10 range group soil /* qx5$~
pro bulk 1e10 she 1e10 range group rock ZY8w1:'
solve fo 10 &l0K~7)b
;位移和速度归零 t=X=",)f
ini xdisp 0 ydisp 0 zdisp 0 h=S7Z:IaM
ini xvel 0 yvel 0 zvel 0 I eJI-lo
;定义材料特性 >|c?ZqW
model mohr \gA<yz-;N
pro bulk 9.1e8 she 6.5e8 co 2.10e4 fric 28 ten 1e6 & '#yqw%
dil 10 range group soil >DUTmJxv
pro bulk 32.5e8 she 35.1e8 co 1.19e6 fric 45 ten 1e6 & n
7i5A:
dil 15 range group rock UOFb.FRP>
hist n=5 mI5J]hk
hist unbal ;:_AOb31N
plo hist 1 1a/C(4_k
set mech ratio 1.0e-6 2Mk;r*FT
solve