层状岩质边坡,命令如下,但最终的solve fos计算不出来,求大神看看哪里出了问题
config dyn
set dyn off
model elas
prop bulk 28.95e9 shear 21.71e9 density 2500 range group a1
prop bulk 29.47e9 shear 22.11e9 density 2500 range group a2
prop bulk 14.11e9 shear 5.88e9 density 2250 range group a3
prop bulk 31.73e9 shear 23.79e9 density 2690 range group a4
prop bulk 62.00e9 shear 23.96e9 density 2700 range group B
prop bulk 65.00e10 shear 27.08e10 density 7500 range group C
gen separate group a1
int 1 wrap a1 a2
gen separate group a2
int 2 wrap a2 a3
gen separate group a3
int 3 wrap a3 a4
int 1 prop ks 2e9 kn 2e9 cohesion 50e4 friction 40
int 2 prop ks 2e9 kn 2e9 cohesion 50e4 friction 40
int 3 prop ks 2e9 kn 2e9 cohesion 50e4 friction 40
set grav 0,0,-10
fix x y z range z=-5
fix x range x=-30
fix x range x=10
fix y
hist unbal
solve
set dyn on
free x y z ran z -5
apply nquiet squiet dquiet ran z -5
table 1 read 20080512.txt
apply xacc 0.01 hist table 1 range z -5
apply ff
set dyn time=0
hist reset
hist dytime
set dyn damp local 0.0628
solve age=4
set dyn off
free range z -5
mod null range x=-40,-30
mod null range x=10,20
mod null range y=-10,0
mod null range y=20,30
fix x y z range z -5
fix x range x -30.1 -29.9 z -4.9 10
fix x range x 9.9 10.1 z -4.9 0
fix y range z -4.9 10
model mohr
apply rem ff
prop bulk 28.95e9 cohesion 60e4 shear 21.71e9 friction 30 density 2500 range group a1
prop bulk 29.47e9 cohesion 60e4 shear 22.11e9 friction 30 density 2500 range group a2
prop bulk 14.11e9 cohesion 20e3 shear 5.88e9 friction 21.0 density 2250 range group a3
prop bulk 31.73e9 cohesion 80e4 shear 23.79e9 friction 33.8 density 2690 range group a4
prop bulk 62.00e9 cohesion 100e4 shear 23.96e9 friction 42.6 density 2700 range group B
prop bulk 65.00e10 cohesion 100e5 shear 27.08e10 friction 45 density 7500 range group C
solve fos