大家帮忙看看我这个程序有什么问题,打算模拟地震作用下,与断层正交的隧道衬砌在地震下的破坏情况,可是模拟的结果很不让人满意,衬砌的相对位移只有0.0几毫米,太小了,根本不符合实际啊,大家帮忙分析,看问题出在哪。
config dyn
set dyn off
gen zone radcyl size 6 5 6 6 ratio 1 1 1 1.1 p0 0 0 0 p1 0.5 0 0 p2 0 0.5 0 p3 0 0 0.5 group rock
gen zone radcyl size 6 3 6 6 ratio 1 1 1 1.1 p0 0 0.5 0 p1 0.5 0.5 0 p2 0 0.8 0 p3 0 0.5 0.5 group siol
gen zone radcyl size 6 5 6 6 ratio 1 1 1 1.1 p0 0 0.8 0 p1 0.5 0.8 0 p2 0 1.3 0 p3 0 0.8 0.5 group rock
gen zon cshell p0 0 0 0 p1 0.1 0 0 p2 0 1.3 0 p3 0 0 0.1 p8 0.08 0 0 p9 0 0 0.08 p10 0.08 1.3 0 p11 0 1.3 0.08 &
size 3 13 6 4 group shel
gen zone reflect dip 90 dd 90 orig 0 0 0
gen zone reflect dip 0 dd 90 orig 0 0 0
gen zone brick size 2 13 6 p0 0.5 0 -0.5 p1 0.6 0 -0.5 p2 0.5 1.3 -0.5 p3 0.5 0 0.5 group pmb
gen zone brick size 2 13 6 p0 -0.6 0 -0.5 p1 -0.5 0 -0.5 p2 -0.6 1.3 -0.5 p3 -0.6 0 0.5 group pmb
interface 1 face range plane norm 0 1 0 origin 0 0.5 0 dist 0.01
interface 2 face range plane norm 0 1 0 origin 0 0.8 0 dist 0.01
interface 3 face range plane norm 1 0 0 origin 0.5 0 0 dist 0.01
interface 4 face range plane norm 1 0 0 origin -0.5 0 0 dist 0.01
attach face cyl rad 0.09 0.11
model el range group rock
prop den 2300 bulk 5e8 shear 3e8 range group rock
model mo range group siol
prop den 2000 bulk 1e7 shear 1e7 fric 10 coh 15000 range group siol
model el range group shel
prop den 2500 bulk 1e9 shear 5e8 range group shel
model el range group pmb
prop den 50 bulk 2e5 shear 1e5 range group pmb
fix x range x -0.61 -0.59
fix x range x 0.59 0.61
fix y range y -0.01 0.01
fix y range y 1.29 1.31
set gravity 0 0 -10
ini xvel 0 yvel 0 zvel 0
ini xdisp 0 ydisp 0 zdisp 0
ini state 0
free x y z range z -0.51 -0.49
set dyn on
set large
table 1 read timeacc.txt
apply xacc=3.0 hist table 1 zacc=0 yacc=0 range z -0.51 -0.49
apply squiet dquiet nquiet range z -0.51 -0.49
set dyn damp rayl 1 20
set dyn dt 5.0e-6
set dyn multi on
his unbal
his dytime
his gp xdisp 0 0.65 0.09
his gp xvel 0 0.65 0.09
his gp xacc 0 0.65 0.09
his gp xdisp 0.09 0.65 0
his gp xvel 0.09 0.65 0
his gp xacc 0.09 0.65 0
his gp xdisp 0 0.65 -0.09
his gp xvel 0 0.65 -0.09
his gp xacc 0 0.65 -0.09
his gp xacc 0 0.65 -0.49
his gp xacc 0 0.65 0.49
plo add surf red
plo sho
solve age 8.0
save mxsy11.sav