[attachment=55214] [attachment=55214] [attachment=55214] 打算模拟与断层正交的隧道在地震作用下的反应，可是开始计算时系统就出错，说有一个节点刚度为零，这是什么原因啊，大家帮忙看看，谢谢
config dyn
gen zone radcyl size 6 20 6 10 ratio 1 1 1 1.1 p0 0 0 0 p1 15 0 0 p2 0 20 0 p3 0 0 15 group rock1
gen zone radcyl size 6 2 6 10 ratio 1 1 1 1.1 p0 0 20 0 p1 15 20 0 p2 0 22 0 p3 0 20 15 group siol
gen zone radcyl size 6 20 6 10 ratio 1 1 1 1.1 p0 0 22 0 p1 15 22 0 p2 0 42 0 p3 0 22 15 group rock2
gen zon cshell p0 0 0 0 p1 3.0 0 0 p2 0 42 0 p3 0 0 3.0 p8 2.5 0 0 p9 0 0 2.5 p10 2.5 42 0 p11 0 42 2.5 &
size 4 42 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
attach face y 19.9 20.1
attach face y 21.9 22.1
attach face cyl rad 2.9 3.1
prop den 2500 bulk 3e10 shear 5e9 range group rock1
mo el
prop den 2000 bulk 1e8 shear 3e7 range group siol
model mohr
prop den 2500 bulk 3e10 shear 5e9 range group rock2
mo el
prop den 2600 bulk 1e11 shear 1e10 range group shel
mo el
fix x range x -15.1 -14.9
fix x range x 14.9 15.1
fix y range y -0.1 0.1
fix y range y 41.9 42.1
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 -15.1 -14.9
table 1 read timeacc.txt
apply xacc=1.0 hist table 1 yacc=0 zacc=0 range z -15.1 -14.9
apply squiet dquiet nquiet range z -15.1 -14.9
set dyn damp local 1.57
set dyn dt 1.0e-4
set dyn multi on
his unbal
his dytime
plo add surf red
plo sho
solve
[ 此贴被wpnjut在2008-09-03 21:31重新编辑 ]
