盾构隧道荷载结构法
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
finish
/clear
*afun,deg
E_c=3.45e7
u_c=0.25
t_c=1.5
r_c=t_c**3/12
density_c=2.5
r0=22.15
k_s=20000
g=9.8
eta=0.7
r_c=r_c*eta
p1=1472.62
p2=p1+density_c*g*3.14*t_c
q1=1336.13
q2=1667.88
/prep7
et,1,beam3
mp,ex,1,E_c
mp,prxy,1,u_c
mp,dens,1,density_c
r,1,t_c,r_c,t_c
csys,1
*do,i,1,37,1
n,i,r0,5*(i-1)+90
*enddo
*do,i,1,36,1
e,i,i+1
*enddo
/pnum,node,1
nplot
csys,0
*do,i,1,36,1
*get,y1,node,i,loc,y
*get,y2,node,i+1,loc,y
f1=q1+(q2-q1)*(r0-y1)/(2*r0)
f2=q1+(q2-q1)*(r0-y2)/(2*r0)
fn1=f1*(y1-y2)/2+(f2-f1)*(y1-y2)/2*1/3
fn2=f1*(y1-y2)/2+(f2-f1)*(y1-y2)/2*2/3
fcum,add
f,i,fx,fn1
f,i+1,fx,fn2
*enddo
*do,i,1,18,1
*get,x1,node,i,loc,x
*get,x2,node,i+1,loc,x
fn1=p1*(x2-x1)/2
fn2=p1*(x2-x1)/2
fcum,add
f,i,fy,fn1
f,i+1,fy,fn2
*enddo
*do,i,19,36,1
*get,x1,node,i,loc,x
*get,x2,node,i+1,loc,x
fn1=p2*(x2-x1)/2
fn2=p2*(x2-x1)/2
fcum,add
f,i,fy,fn1
f,i+1,fy,fn2
*enddo
csys,1
*do,i,10,27
*get,len,elem,i,leng
fcum,add
psprng,i,tran,0.5*len*k_s,1
psprng,i+1,tran,0.5*len*k_s,1
*enddo
csys,0
d,1,ux,0
d,1,rotz,0
d,37,ux,0
d,37,rotz,0
finish
/solu
acel,,g
solve
finish
/post1
pldisp,1
etable,m1,smisc,6
etable,m2,smisc,12
plls,m1,m2,-1
etable,n1,smisc,1
etable,n2,smisc,7
plls,n1,n2,-1
etable,q1,smisc,2
etable,q2,smisc,8
plls,q1,q2,-1