请教隧道开挖模拟前原始地应力场的求解问题

yan1984

在求解原始地应力后，怎么得到的节点力很小（大概5e-6），我看书上面的一般也是10的正五六次方啊。不知道问题出在哪里？请好心人指导一下。不胜感激！
命令流如下：/prep7
et,1,plane42
keyopt,1,3,2
mp,ex,1,1.0e9       !土层材料属性
mp,prxy,1,0.35
mp,dens,1,2000
tb,dp,1,,,
tbmodif,1,1,1.2e6
tbmodif,1,2,45
mp,ex,2,3.0e10       !初期支护材料属性 ，C30
mp,prxy,2,0.2
mp,dens,2,2700
r,2,0,0,1,0.01
k,1,-5.1345,3.6537
k,2,-6.0626,-3.4998
k,3,-4.7979,-4.7556
k,4,4.7979,-4.7556
k,5,6.0626,-3.4998
k,6,5.1345,3.6537
larc,1,2,5,8.06
larc,2,3,6,2
larc,3,4,6,15.25
larc,4,5,1,2
larc,5,6,2,8.06
larc,6,1,3,6.25
LOCAL,11,1,0,0,0, , , ,1,1,
LOCAL,12,1,1.57,-0.91,0, , , ,1,1,
LOCAL,16,1,-1.57,-0.91,0, , , ,1,1,
LOCAL,14,1,0,9.72,0, , , ,1,1,
LOCAL,13,1,-3.95,-2.7,0,,,,1,1,
LOCAL,15,1,3.95,-2.7,0,,,,1,1,
*do,i,1,6
csys,10+i
k,i+6,kx(i)+0.1,ky(i)
*enddo
CSYS,0,
CSDELE,11, ,1,  
CSDELE,12, ,1,  
CSDELE,13, ,1,  
CSDELE,14, ,1,  
CSDELE,15, ,1,  
CSDELE,16, ,1,
larc,7,8,11,8.16
larc,8,9,12,2.1
larc,9,10,12,15.35
larc,10,11,7,2.1
larc,11,12,8,8.16
larc,12,7,9,6.35
a,1,2,3,4,5,6
a,7,8,9,10,11,12
blc4,-20,-15,40,35
blc4,-30,-15,10,35
blc4,20,-15,10,35
blc4,-30,-30,60,15
aovl,1,2,3,4,5,6
nummrg,all,,,,low            !合并重复元素并保留低编号号码
numcmp,all
l,13,8
l,14,11
l,15,12
l,16,7
lsel,s,line,,26,29,3
lsel,a,line,,7
asbl,5,all
nummrg,all,,,,low            !合并重复元素并保留低编号号码
numcmp,all
lsel,s,line,,27,28,1
lsel,a,line,,11
asbl,7,all
nummrg,all,,,,low            !合并重复元素并保留低编号号码
numcmp,all
lsel,s,line,,8,10,1
lccat,all
lsel,s,line,,8,10,2
lesize,all,,,2
lsel,s,line,,9
lesize,9,,,8
lsel,s,line,,7,11,4
lesize,all,,,10
lsel,s,line,,14,16,2
lesize,all,,,10
lsel,s,line,,12,15,3
lsel,a,line,,13
lesize,all,,,12
lsel,s,line,,26,29,1
lesize,all,,,10,2
asel,s,area,,6,9,1
amesh,all
lsel,s,line,,2,4,2   !划分最里面的两个面
lesize,all,,,2
lsel,s,line,,3
lesize,3,,,8
lsel,s,line,,1,5,4
lesize,all,,,10
lsel,s,line,,6
lesize,all,,,12
asel,s,area,,1,5,4
amesh,all
lsel,s,line,,20,23,1  !划分外面的三个矩形
lesize,all,,,4,2
lsel,s,line,,17,18,1
lesize,all,,,10
lsel,s,line,,24,25,1
lesize,all,,,5,2
lsel,s,line,,20,22,2
lsel,a,line,,13
lccat,all
lsel,s,line,,19
lesize,all,,,20
asel,s,area,,2,4,1
amesh,all
lsel,s,lccat
ldele,all
allsel
nummrg,all,,,,low
numcmp,all
finil
/solu
antype,static
lsel,s,loc,x,-30
lsel,a,loc,x,30
nsll,s,1
d,all,ux,0
allsel
lsel,s,loc,y,-30
nsll,s,1
d,all,uy,0
allsel
acel,,10
deltim,0.1,0.05,0.2  !设定分析选项
autots,on
pred,on
lnsrch,on
nlgeom,on
nropt,full
cnvtol,f,,0.02,2,0.5
save
time,1
solve                   !求解原始地应力场
/post1
lsel,s,line,,7,12,1
nsll,s,1
nforce,all              !后处理得到开挖截面的节点力

[ 本帖最后由 yan1984 于 2007-10-26 12:59 编辑 ]