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]