ansys实例——混凝土加弹簧

!定义相关参量*DIM,MK,ARRAY,1,4,,,,!定义水平抗力系数比例常数m的矩阵*SET,MK(1,1,1),6000!第四层土的m值*SET,MK(1,2,1),5000!第三层土的m值*SET,MK(1,3,1),4000!第二层土的m值*SET,MK(1,4,1),2000!第一层土的m值*DIM,FQ,ARRAY,1,4,,,,!定义桩周单位面积摩擦力*SET,FQ(1,1,1),40.53!中性点以下的单位面积正摩擦力*SET,FQ(1,2,1),-33.78!中性点以上第三层的单位面积负摩擦力*SET,FQ(1,3,1),-16.14*SET,FQ(1,4,1),-4.8*SET,QX,100!水平荷载*SET,QZ,100*SET,b0,1*SET,h0,0.5!进入前处理器/prep7et,1,solid45!桩身单元et,2,combin14!弹簧单元et,3,surf154!表面效应单元keyopt,1,1,0keyopt,2,3,0!作为纵向三维弹簧keyopt,3,4,1!不包含中间节点r,1mp,ex,1,3.25e10!桩的材料属性mp,nuxy,1,0.3mp,dens,1,2700!建立整个桩模型的节点n,1n,6,1fill,1,6ngen,6,6,1,6,1,,,0.2ngen,21,36,1,36,1,,0.5/view,1,1,1,1!以正等侧显示/ang,1/REP,FASTGPLOT!建立桩的几何实体模型k,1k,2,1k,3,1,,1k,4,,,1k,5,,10,k,6,1,10,k,7,1,10,1k,8,,10,1/PNUM,kp,1kplotv,1,2,3,4,5,6,7,8lsel,s,loc,y,0lsel,a,loc,y,10lesize,all,,,5allsellsel,s,loc,y,5lesize,all,,,20vmesh,1nummrg,all!压缩重复元素numcmp,all!合并重复元素!建立第一组弹簧模型,共7*16=102个弹簧nsel,s,loc,x,1nsel,r,loc,y,-0.1,8.1nplotngen,2,1000,all,,,2type,2*do,i,1,17,1*if,i,le,4,then!判断属于哪一层土kk=1*elseif,i,le,8,thenkk=2*elseif,i,le,12,thenkk=3*elseif,i,le,17,thenkk=4*endif*do,j,6,36,6*get,vy,node,j+(i-1)*36,loc,y!获取节点的Y坐标值r,kk+1,mk(1,kk,1)*abs(vy-8)*b0*h0/6!根据Y值计算水平抗力系数，即弹簧的刚度real,kk+1e,j+(i-1)*36,1000+j+(i-1)*36*enddo*enddoallsel!最大单元号602!第二组弹簧单元,共7*16=102个弹簧allselnsel,s,loc,x,0nsel,r,loc,y,-0.1,8.1nsel,r,loc,z,-0.1,1.1nplotngen,2,2000,all,,,-2type,2*do,i,1,17,1*if,i,le,4,thenkk=1*elseif,i,le,8,thenkk=2*elseif,i,le,12,thenkk=3*elseif,i,le,17,thenkk=4*endif*do,j,1,31,6*get,vy,node,j+(i-1)*36,loc,yr,kk+10,mk(1,kk,1)*abs(vy-8)*b0*h0/6real,kk+10e,j+(i-1)*36,2000+j+(i-1)*36*enddo*enddoallsel!最大单元号704!第三组弹簧单元,共7*16=102个弹簧local,11,0,,,,,90csys,0allselnsel,s,loc,z,1nsel,r,loc,y,-0.1,8.1nsel,r,loc,x,-0.1,1.1nplotngen,2,3000,all,,,,,2type,2*do,i,1,17,1*if,i,le,4,thenkk=1*elseif,i,le,8,thenkk=2*elseif,i,le,12,thenkk=3*elseif,i,le,17,thenkk=4*endif*do,j,31,36,1*get,vy,node,j+(i-1)*36,loc,yr,kk+20,mk(1,kk,1)*abs(vy-8)*b0*h0/6real,kk+20esys,0e,j+(i-1)*36,3000+j+(i-1)*36*enddo*enddoallsel!最大单元号码806!第四组弹簧单元,共7*16=102个弹簧allselnsel,s,loc,z,0nsel,r,loc,y,-0.1,8.1nsel,r,loc,x,-0.1,1.1nplotngen,2,4000,all,,,,,-2type,2*do,i,1,17,1*if,i,le,4,thenkk=1*elseif,i,le,8,thenkk=2*elseif,i,le,12,thenkk=3*elseif,i,le,17,thenkk=4*endif*do,j,1,6,1*get,vy,node,j+(i-1)*36,loc,yr,kk+20,mk(1,kk,1)*abs(vy-8)*b0*h0/6real,kk+20esys,0e,j+(i-1)*36,4000+j+(i-1)*36*enddo*enddoallsel!最大单元号码为908!底层弹簧allselcsys,0ngen,2,5000,1,36,1,,-1type,2r,1000,mk(1,1,1)*8*b0*h0/6*20real,1000*do,i,1,36,1e,i,5000+i*enddoallselgplot!第一组表面效应单元allselr,1001real,1001local,13,0,,,,90type,3*do,i,1,16,1*do,j,6,30,6nsel,s,node,,j+(i-1)*36nsel,a,node,,j+6+(i-1)*36nsel,a,node,,j+42+(i-1)*36nsel,a,node,,j+36+(i-1)*36esys,13esurf*enddo*enddo!第二组表面效应单元allsellocal,13,0,,,,90type,3*do,i,1,16,1*do,j,1,29,6nsel,s,node,,j+(i-1)*36nsel,a,node,,j+6+(i-1)*36nsel,a,node,,j+42+(i-1)*36nsel,a,node,,j+36+(i-1)*36esys,13esurf*enddo*enddo!第三组表面效应单元allsellocal,14,0,,,,90type,3*do,i,1,16,1*do,j,1,5,1nsel,s,node,,j+(i-1)*36nsel,a,node,,j+1+(i-1)*36nsel,a,node,,j+37+(i-1)*36nsel,a,node,,j+36+(i-1)*36esys,14esurf*enddo*enddocsys,0!第四组表面效应单元allsellocal,15,0,,,,90type,3*do,i,1,16,1*do,j,31,35,1nsel,s,node,,j+(i-1)*36nsel,a,node,,j+1+(i-1)*36nsel,a,node,,j+37+(i-1)*36nsel,a,node,,j+36+(i-1)*36esys,15esurf*enddo*enddosavefinish!进入求解器/solucsys,0allselnsel,s,loc,x,3nsel,a,loc,x,-2nsel,a,loc,z,3nsel,a,loc,z,-2d,all,all!约束边界条件nsel,a,loc,y,-1d,all,allnsel,a,loc,y,0d,all,uxd,all,uzACEL,0,9.81,0,!施加重力加速度!桩顶表面平均分布载荷allselnsel,s,loc,y,10nsel,r,loc,x,-0.1,1.1nsel,r,loc,z,-0.1,1.1nplotesln,ssfe,all,6,pres,,5000!桩顶表面分布荷载nsel,r,loc,x,1f,all,fx,QX!水平载荷nsel,s,loc,y,10nsel,r,loc,z,1f,all,fz,QZ!水平载荷allselgplot!第四层土的表面正摩擦力csys,0allselnsel,s,loc,y,0,1.5esln,sesel,r,type,,3sfe,all,2,pres,,FQ(1,1,1)/PSF,PRES,TANx,2,0,1gplot!第三层土的表面负摩擦力csys,0allselnsel,s,loc,y,2.5,3.5esln,sesel,r,type,,3sfe,all,2,pres,,FQ(1,2,1)/PSF,PRES,TANx,2,0,1gplot!第二层土的表面负摩擦力csys,0allselnsel,s,loc,y,4.5,5.5esln,sesel,r,type,,3sfe,all,2,pres,,FQ(1,3,1)/PSF,PRES,TANx,2,0,1gplot!第一层土的表面负摩擦力csys,0allselnsel,s,loc,y,6.5,7.5esln,sesel,r,type,,3sfe,all,2,pres,,FQ(1,4,1)/PSF,PRES,TANx,2,0,1gplotallselgplottime,1autot,onnlgeom,onnropt,fullcsys,0solvefinish/POST1!PLNSOL,U,Y,1,1!PLNSOL,S,Z,1,1PATH,path1,2,30,20,PPATH,P51X,1PATH,STAT/PBC,PATH,1!/REPLOT/PBC,PATH,0AVPRIN,0,0,PDEF,SY,S,Y,AVG/PBC,PATH,,0PRANGE,1,1,10,YG,PLPATH,SY