网壳整体稳定分析命令流.docxVIP

网壳整体稳定分析命令流 !以一杆二单元网壳模型(Domel?2)分析命令流为例 FINISH /Clear /title,A nsysWorkforDomel?2 /Prep7 !定义基本儿何参数!初始化设计变量参数 Span=70!定义跨度为70米 Lf二0.25!定义矢跨比 f=Span*Lf!矢高 R=f/2+Spa n *Spa n/§/f!半径 Kn二8!默认网壳类型K8 Pi=3.14159!圆周率值 Nx=12!默认环向圈数12圈 DEAD=300!tH.荷载;数值 LIVE二500*30!活荷载数值，放人(17*0.80445*0.992)倍 ND=2!~根杆件划分单元 !定义单元类型和参数 MP,EX,l,2.060e+011 M^NUXYJ^.OOOe-OOl MP,DENS,l,7.850e+003 ET；1,BEAM188 SECTYPEZ1,BEAM,CTUBEZ/O SECOFFSET,CENT SECDATA,0.061500,0.063500,8 TYPE,1 MAT,1 SECNUM,1 !计算节点坐标位置,并定义关键点 DPha=Atn(Span/ySqrt(R*R-Span*Span/l))/Nx K,1,0,0,R *do,i,l/Nx *do,j,l,Kn*i x=R*sin(i*DPha)*cos((j-l)/i*2*Pi/Kn) y=R*sin(i*DPha)*sin((j-l)/i*2*Pi/Kn) z=R*cos(i*DPha) K/l+Kn*(i-l)*i/2+j/x/y/z *enddo *enddo !定义连接 *DO,i,l,K n *ENDDO !环向的单元连接 *DO,j,l,Nx *DO,i,l,Kn*j-l L/l+Kn*(j-l)*j/2+izl+Kn*(j-l)*j/2+i+l *ENDDO Lzl+Kn*(j-l)*j/2+l/l+Kn*(j-l)*j/2+Kn*j *ENDDO !径向的单元连接 *Do,j,l,Nx-l *Do,k,l,K n *lfj,eqj+l,the n *lf,k,eq,Kn/The n L/l+Kn*(j-l)*j/2+l/l+Kn*(j+l)*j/2+(k-l)*(j+l)+i *else L,l+Kn*(j-l)*j/2+(k-l)*j+i/l+Kn*(j+l)*j/2+(k-l)*(j+l)+l *Endif *Else Lzl+Kn*(j-l)*j/2+(k-l)*j+i/l+Kn*(j+l)*j/2+(k-l)*(j+l)+l *Endlf *Enddo Ll+Kn*(j-l)*j/2+(k-l)*j+U+Kn*(j+l)*j/2+(k-l)*(j+l)+i+l *Enddo *Enddo *Enddo LPLOT !定义边界约束和恒载 *Do/i,l,l+Kn*(Nx-l)*Nx/5+Kn*Nx *lf,i/lt/l+Kn*(Nx-l)*Nx/2+l/the n FKJ,fzz-DEAD*Pi*(SPAN/2)*(SPAN/2)/(l+Kn*(Nx-l)*Nx/2+Kn*Nx) *Else DK,i,UX DK,i,UY DK,I,UZ *Endif *Enddo !定义活荷载 FK,bfz,(?LIVE?DEAD)*P 产(SPAN/2)*(SPAN/2)/(l+Kn*(Nx-l)*Nx/2+Kn*Nx) *Doj,1,Nx-1 *Do/i/(j-l)*j*Kn/2+2z(j*j*Kn+4)/2 FKJ/fz/(-LIVE-DEAD)*Pi*(SPAN/2)*(SPAN/2)/(l+Kn*(Nx-l)*Nx/2+Kn*Nx) *Enddo *Enddo 戊I」分单元 lesize,alL,ND LPLOT Lmesh,all ACEL,0,0,9.8, EPLOT !(1)线性整体稳定分析 !Domel-2 FINISH /solu an type,0 pstres, on solve FINISH /solu an type,l bucoptja nb‘20 mxpa nd,20 outres,all,all solve finish !查看求解结果 /postl setjist set, 1,1 !(2)完善结构大位移几何非线性整体稳定分析 FINISH /solu antype,0!分析类型为静态分析 pstres,on!预应力开关 nlgeom,on!考虑大变形效应 arclen,on!激活弧长法 arctrm,U,0.5?uz!对弧长法求解进彳亍终止控制 neqit,15!容许的最人平衡迭代次数(15) nsubst,10!设置子步数(10) outres,all,all!输出结果类型及其频度 sol