单层网壳整体稳定性分析-ANSYS命令流.docx
《单层网壳整体稳定性分析-ANSYS命令流.docx》由会员分享，可在线阅读，更多相关《单层网壳整体稳定性分析-ANSYS命令流.docx（8页珍藏版）》请在冰豆网上搜索。
单层网壳整体稳定性分析-ANSYS命令流
网壳整体稳定分析命令流
第8页
!
以一杆二单元网壳模型（Dome1-2）分析命令流为例
FINISH
/Clear
/title,AnsysWorkforDome1-2
/Prep7
!
定义基本几何参数!
初始化设计变量参数
Span=70!
定义跨度为70米
Lf=0.25!
定义矢跨比
f=Span*Lf!
矢高
R=f/2+Span*Span/8/f!
半径
Kn=8!
默认网壳类型K8
Pi=3.14159!
圆周率值
Nx=12!
默认环向圈数12圈
DEAD=300!
恒荷载数值
LIVE=500*30!
活荷载数值，放大（17*0.80445*0.992）倍
ND=2!
一根杆件划分单元
!
定义单元类型和参数
MP,EX,1,2.060e+011
MP,NUXY,1,3.000e-001
MP,DENS,1,7.850e+003
ET,1,BEAM188
SECTYPE,1,BEAM,CTUBE,,0
SECOFFSET,CENT
SECDATA,0.061500,0.063500,8
TYPE,1
MAT,1
SECNUM,1
!
计算节点坐标位置,并定义关键点
DPha=Atn（Span/2/Sqrt（R*R-Span*Span/4））/Nx
K,1,0,0,R
*do,i,1,Nx
*do,j,1,Kn*i
x=R*sin（i*DPha）*cos（（j-1）/i*2*Pi/Kn）
y=R*sin（i*DPha）*sin（（j-1）/i*2*Pi/Kn）
z=R*cos（i*DPha）
K,1+Kn*（i-1）*i/2+j,x,y,z
*enddo
*enddo
!
定义连接
*DO,i,1,Kn
L,1,1+i
*ENDDO
!
环向的单元连接
*DO,j,1,Nx
*DO,i,1,Kn*j-1
L,1+Kn*（j-1）*j/2+i,1+Kn*（j-1）*j/2+i+1
*ENDDO
L,1+Kn*（j-1）*j/2+1,1+Kn*（j-1）*j/2+Kn*j
*ENDDO
!
径向的单元连接
*Do,j,1,Nx-1
*Do,k,1,Kn
*Do,i,1,j+1
*If,i,eq,j+1,then
*If,k,eq,Kn,Then
L,1+Kn*（j-1）*j/2+1,1+Kn*（j+1）*j/2+（k-1）*（j+1）+i
*else
L,1+Kn*（j-1）*j/2+（k-1）*j+i,1+Kn*（j+1）*j/2+（k-1）*（j+1）+I
*Endif
*Else
L,1+Kn*（j-1）*j/2+（k-1）*j+i,1+Kn*（j+1）*j/2+（k-1）*（j+1）+I
*EndIf
*Enddo
*Do,i,1,j
L,1+Kn*（j-1）*j/2+（k-1）*j+i,1+Kn*（j+1）*j/2+（k-1）*（j+1）+i+1
*Enddo
*Enddo
*Enddo
LPLOT
!
定义边界约束和恒载
*Do,i,1,1+Kn*（Nx-1）*Nx/2+Kn*Nx
*If,i,lt,1+Kn*（Nx-1）*Nx/2+1,then
FK,i,fz,-DEAD*Pi*（SPAN/2）*（SPAN/2）/（1+Kn*（Nx-1）*Nx/2+Kn*Nx）
*Else
DK,i,UX
DK,i,UY
DK,I,UZ
*Endif
*Enddo
!
定义活荷载
FK,1,fz,（-LIVE-DEAD）*Pi*（SPAN/2）*（SPAN/2）/（1+Kn*（Nx-1）*Nx/2+Kn*Nx）
*Do,j,1,Nx-1
*Do,i,（j-1）*j*Kn/2+2,（j*j*Kn+4）/2
FK,i,fz,（-LIVE-DEAD）*Pi*（SPAN/2）*（SPAN/2）/（1+Kn*（Nx-1）*Nx/2+Kn*Nx）
*Enddo
*Enddo
!
划分单元
lesize,all,,,ND
LPLOT
Lmesh,all
ACEL,0,0,9.8,
EPLOT
!
（1）线性整体稳定分析
!
Dome1-2
FINISH
/solu
antype,0
pstres,on
solve
FINISH
/solu
antype,1
bucopt,lanb,20
mxpand,20
outres,all,all
solve
finish
!
查看求解结果
/post1
set,list
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!
输出结果类型及其频度
solve
finish
/post26
nsol,2,29,u,z,uz_29
prod,3,2,,,,,,-1000
prod,4,1,,,,,,30
/axlab,x,Uz（mm）
/axlab,y,K
xvar,3
plvar,4
!
（3）考虑初始缺陷的大位移几何非线性整体稳定分析
finish
/solu
nlgeom,off
arclen,off
antype,0
pstres,on
solve
finish
/solu
antype,buckle
bucopt,lanb,5
mxpand,5
solve
finish
/prep7
UPGEOM,-0.2333,LAST,LAST,'AnsysWorkforDome','rst',''
finish
/solu
antype,0
nlgeom,on
arclen,on
arctrm,U,1,,uz!
对弧长法求解进行终止控制
neqit,15!
容许的最大平衡迭代次数（15）
nsubst,20!
设置子步数（20）
outres,all,all
solve
finish
/post26
nsol,2,69,u,z,uz_69!
33,67,69,121
prod,3,2,,,,,,-1000
prod,4,1,,,,,,30
/axlab,x,Uz（mm）
/axlab,y,K
xvar,3
plvar,4
!
（4）弹塑性整体稳定分析
finish
/solu
antype,0
pstres,on
nlgeom,off
arclen,off
solve
finish
/solu
antype,buckle
bucopt,lanb,5
mxpand,5
solve
finish
/prep7
UPGEOM,-0.2333,LAST,LAST,'AnsysWorkforDome','rst',''
finish
/prep7
TB,BKIN,1
TBDATA,1,235E6,7.9E8!
Tbplot,bkin,1
FINISH
/solu
antype,0
nlgeom,on
arclen,on
arctrm,u,1,,uz
neqit,20!
容许的最大平衡迭代次数（20）
nsubst,500!
设置子步数（500）
outres,all,all
solve
finish
/post26
nsol,2,69,u,z,uz_69!
33,67,69,71,121
prod,3,2,,,,,,-1000
prod,4,1,,,,,,30
/axlab,x,Uz（mm）
/axlab,y,K
xvar,3
plvar,4
升级会员 