盾构管片修正惯用法内力计算 ansys命令流.docVIP

!匀质圆环,在原程序上加了静水压力，故本命令流适合水土分算 fini /cle *dim,wxn,array,2000 !!定义名为wxn的数组，2000行1列，下同 *dim,wyn,array,2000 *dim,xn,array,2000 *dim,yn,array,2000 *dim,fxn,array,2000 *dim,fyn,array,2000 *dim,aa1,array,2000 *dim,bb1,array,2000 *dim,aa2,array,2000 *dim,bb2,array,2000 *dim,jx1,array,2000 *dim,jy1,array,2000 *dim,jx2,array,2000 *dim,jy2,array,2000 !!!!!!!!!!定义参数并赋值 nodesum=120 !!!!!!!!!!单元数 pi=3.1415927 height=0.35 !!!!!!环厚0.35米 width=1.0 !!!!!!环宽 area=height*width !!!!!!面积 inertia=width*height*height*height/12 !!!!!!惯性矩 emod=34500*1000000 !!!!!!弹性模量，按C50混凝土计 radius=2.925 !!!!!!计算半径 density=25*1000/10 !!!!!!材料密度 wradius=radius+0.5 !!!!!!弹簧节点所处半径 tankang=12*1000000*width !!!!!!弹簧刚度,MPa/m !!!!!!!!!!赋荷载值 ptop=200*1000*width !!!!!!顶板压力，单位按KN计 pltop=80*1000*width !!!!侧向压力 plbot=50*1000*width !!!!侧向沿竖向增加的侧压，即墙底水平压力为pltop+plbot pbot=140*1000*width !!!!底板压力 water=8.0*1000*width !!!!顶板处水位埋深 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! /prep7 et,1,beam3 !!!!!!定义单元类型,弹模,密度，实常数 mp,ex,1,emod mp,dens,1,density mp,prxy,1,0.2 r,1,area,inertia,height et,2,combin39 !!!!!定义弹簧单元 keyopt,2,1,0 !!!!0为沿与加载相同的路径卸载 keyopt,2,2,0 !!!!按照所定义的受压曲线加压力荷载 keyopt,2,3,0 !!!!沿节点X轴方向平动 keyopt,2,4,3 !!!!二维轴向单元 keyopt,2,6,0 !!!!基本单元输出 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! *do,i,1,nodesum,1 xn(i)=radius*cos(pi/2+(i-1)*(2*pi/nodesum)) !!!!!!定义管片节点I的x坐标 yn(i)=radius*sin(pi/2+(i-1)*(2*pi/nodesum)) !!!!!!定义管片节点I的y坐标 n,i,xn(i),yn(i) !!!!!!生成管片节点I *enddo *do,i,1,nodesum,1 wxn(i)=wradius*cos(pi/2+(i-1)*(2*pi/nodesum)) !!!!!!定义弹簧节点I的x坐标 wyn(i)=wradius*sin(pi/2+(i-1)*(2*pi/nodesum)) !!!!!!定义弹簧节点I的y坐标 n,i+nodesum,wxn(i),wyn(i) !!!!!!生成弹簧节点I *enddo !!!!!!!!!!!生成单元，类型为1，实常数为1 type,1 real,1 *do,i,1,nodesum-1,1 e,i,i+1 *enddo e,nodesum,1 !!!!!!!!!!!生成弹簧单元 len=radius*(2*pi/nodesum) wytankang=len*tankang r,2,-1,-wytankang,0,0,1,wytankang/1000000 !!!!!!! type,2 real,2 !!!!!!!!!!!!!!!! *do,i,