?——内摩擦角
P0——上覆荷载
?——土的容重
c——土的粘聚力
三、管片内力有限元计算单元选择
1、beam3单元
输入参数：
（1）实常数（Real Constants） AREA – 横截面积 IZZ – 横截面惯性矩 HEIGHT – 梁高
（2）材料参数（Material Properties） EX——弹模 DENS——密度 PRXY——泊松比 其它选项采用默认即可 2、Link10单元
（1）Real Constants
AREA – 横截面积，取单位面积1。 （2）Material Properties
EX——弹模，需满足E*A=地层抗力刚度 其它选项采用默认即可 3、link1单元
（1）Real Constants
AREA – 横截面积，取单位面积1。 （2）Material Properties
EX——弹模，为了不影响计算结果取足够小，如1。
四、结点力计算
1、均布荷载时
中间结点：
Fn?p(xn?1?xn?1)/2
两端结点：
F1?p(x2?x1)/2 FN?p(xN?xN?1)/2
2、梯形荷载时
中间结点：
1Fn??(pn?1?2pn)(yn?yn?1)?(2pn?pn?1)(yn?1?yn)?
6两端结点：
1F1??(2p1?p2)(y2?y1)?
61FN??(pN?1?2pN)(yN?yN?1)?
6
五、圆形结构内力计算
1、计算参数
R=4.5m；B=1.5m；H=0.5m(厚)； Pw1=0.1MPa;Pe1=0.14MPa; Pw2=0.19MPa;Pe2=0.05MPa; Q1=0.14MPa Q2=0.23MPa
2、在自重情况下的结构内力分布
计算命令流为： /prep7 et,1,beam3 et,2,link10
keyopt,2,3,1 !*****只受压****** R,1,0.75,0.015625,0.5, , , , R,2,1, ,
!定义混凝土衬砌材料 mp,ex,1,2.76e10 mp,prxy,1,0.2 mp,dens,1,2600 mp,ex,2,5e7 mp,prxy,2,0.01 !建立几何模型 k,1,0,0,0 circle,1,4.5
!设置划分单元大小为0.5m/个 Lsel,all
LESIZE,all,0.5, , , , , , ,1 Type,1 Mat,1
Real,1 Lmesh,all
!复制建立地层弹簧 CSYS,1 wpcsys,-1 nsel,s,loc,x,4.4,4.6 NGEN,2,500,all, , ,1 , , ,1, allsel,all *do,i,1,60 type,2 real,2 mat,2 e,i,i+500 *enddo
nsel,s,loc,x,5.4,5.6 d,all,all nsel,s,,,32
D,all, , , , , ,UX, , , , , allsel,all csys,0 wpcsys,-1 /sol
ACEL,0,10,0, !********施加重力命令********* solve /post1
ETABLE, ,SMISC, 1 ETABLE, ,SMISC, 7 ETABLE, ,SMISC, 2 ETABLE, ,SMISC, 8 ETABLE, ,SMISC, 6