lianglcy2008 发表于 2012-7-25 17:09:39

ansys计算地下隧洞衬砌实例问题

我做的地下隧洞衬砌实例,beam3和combin14的组合使用
遇到的问题是衬砌全部加弹簧,计算出来的竟是全部受拉,也就是说要去掉全部弹簧,和实际情况不符合,请高手帮忙看看啊,以下是命令流
/CWD,E:\LINER\               !改变工作目录

fini
/cle
/filname,HV TUNNEL LINER
/titile,analysis of HV TUNNEL LINER
!*create,liner,mac

!!!!!!!!!!!!!!!!!!!!!!!!!!!!定义单元类型,实常数

/prep7


et,1,beam3
et,2,mesh200
et,3,combin14

r,1,0.3,0.00225,0.3
r,2,2246694632

mp,ex,1,24.87e9                                 !!!!!!!!!混凝土材料
mp,prxy,1,0.25
mp,dens,1,2400

k,1,-1.9,0
k,2,-1.9,5.11
k,3,0,7.01
k,4,1.9,5.11
k,5,1.9,0
k,6,-2.9,0
k,7,-2.9,5.11
k,8,0,8.01
k,9,2.9,5.11
k,10,2.9,0

l,1,2
larc,2,3,4,1.9
larc,3,4,2,1.9
l,4,5
l,6,7
larc,7,8,9,2.9
larc,8,9,7,2.9
l,9,10

lsel,s,line,,1,4
latt,1,1,1                              !!!!!!!!!!!材料号,实常数,单元类型号

lsel,s,line,,5,8
latt,,,2

alls
lesize,all,,,20,,,,,
lmesh,all
/eshape,1

dk,1,ux,0
dk,1,uy,0
dk,5,ux,0
dk,5,uy,0

nsel,s,node,,82,162,1
d,all,all

type,3            !!!!!!!!!!指定弹簧单元
real,2

*do,I,82,162,1
e,I,I-81
*enddo

SFCUM,PRES,ADD   !!!!!设置力为叠加
FCUM,ADD

!!!!!!!!!!!!!!!!!!!!!!!!荷载
!!!!!!!!!!!!!!!!!!!!!!!!!顶拱
rrp=166.20*1000      !!!!!!!顶拱岩石压力Pa            

lsel,s,line,,2,3
esll,s
sfbeam,all,1,pres,rrp,rrp      

rwp=25.11*1000         !!!!!!!顶拱水压力Pa            

lsel,s,line,,2,3
esll,s
sfbeam,all,1,pres,rwp,rwp

rgp=200*1000      !!!!!!!顶拱灌浆压力Pa            

lsel,s,line,,2,3
esll,s
sfbeam,all,1,pres,rgp,rgp

!!!!!!!!!!!!!!!!!!!!!!!边墙
wrp=42.6*1000                  !!!!!!!!!!!!!!边墙岩石压力Pa

alls
*do,I,20,1,-1
sfbeam,I,1,pres,wrp,wrp
*ENDDO
*do,I,61,80,1
sfbeam,I,1,pres,wrp,wrp
*ENDDO

w1=50.11*1000                  !!!!!!!!!!!!!!拱座水压力Pa
w2=68.74*1000                  !!!!!!!!!!!!!!墙底水压力Pa
d=(w2-w1)/20
*do,I,20,1,-1
sfbeam,I,1,pres,(20-(I-1))*d+w1,(20-I)*d+w1
*ENDDO

*do,I,61,80,1
sfbeam,I,1,pres,(I-61)*d+w1,(I-60)*d+w1
*ENDDO

alls
dtran

csys,0
alls
acel,,9.8,

/solu
antype,static
NROPT, FULL,,                       !采用全牛顿-拉普森法进行求解
alls
solve
finish
页: [1]
查看完整版本: ansys计算地下隧洞衬砌实例问题