ANSYS非线性分析

跨次元

各位i大神帮忙看一下吧，非线性分析始终做不出结果来，因为是初学者对于非线性分析和弧长法不是很懂，拜托了！

跨次元

! ******************************************************
!!!!!!!!!!!三级飞行段弹体参数化建模强度分析模型
!!!!!!!!!!!对①材料参数；②几何参数；③载荷大小；
!!!!!!!!!!!④网格尺寸进行了参数化。
! ******************************************************
finish
/clear
/uis,msgpop,3                          !控制运算时只显示错误，不显示警告
/rgb,index,100,100,100,0                 !!!!!!!调换为ANSYS白色界面背景
/rgb,index,80,80,80,13
/rgb,index,60,60,60,14
/rgb,index,0,0,0,15
/uis,msgpop,3  
/FILNAME,cylindricity,1                    !不显示警告信息
/prep7
!******************************************************
!!!!!!!!!!!材料参数化
!******************************************************
!!!!材料力学参数化
E1=71e9
!FYX=15E7
V1=0.3
dens1=2670
!******************************************************
!!!!几何模型参数化
!******************************************************
!!!!!!定义各段特征尺寸
T=0.001                 !弹体厚度
lh=0.72                  !圆柱壳长度
rc=0.5                       !圆柱壳半径
R_laser=0.1              !光斑半径

!******************************************************
!!!!有限元网格尺寸参数化
!******************************************************
esize1=0.02                   !圆柱壳网格尺寸

!******************************************************
!!!!载荷大小参数化
!******************************************************

!******************************************************
!!!!!!!!!!!设置单元，并定义单元属性
!******************************************************
/prep7
*afun,deg                         !指定角度为度数单位
et,1,shell181                      !指定单元类型181
keyopt,1,8,2
!******************************************************
!!!!定义各向同性材料
!******************************************************
MP,EX,1,E1
MP,PRXY,1,V1
MP,dens,1,dens1
sect,1,shell,,
secdata,0.0015,1,0.0,5
secoffset,MID
seccontrol,,,,,,,
!******************************************************
!!!!!建立模型
!******************************************************
!!!!!建立模型
!******************************************************
!!建立关键点
K,1
K,2,rc
K,3,rc,lh
K,4,,lh

L,2,3

arotat,all,,,,,,1,4,360,4
asel,all
aglue,all
save
wpave,, lh /2
Asel,all
cm,a1,area
cswpla,11,1
csys,11
wpcsys,-1,0
csys,0
asel,all
allsel
aglue,all
save
wpcsys,-1,0
asbw,all

wpcsys,-1,0
csys,0
save

!******************************************************
!!!!!!划分网格
!******************************************************

!!!为关键点划分网格
allsel
mat,1              !指定材料的类型
type,1              !指定单元的类型
real,1              !指定实常数

allsel
lsel,all
lesize,all, esize1,,,,,,,1                       !定义单元的边长尺寸
allsel

asel,s,all
aesize,all,esize1
mshkey,1
mshape,0
amesh,all
allsel,all
save

wpcsys,-1,0
csys,0
!!!!!!施加力学边界条件
!******************************************************
/solu
ANTYPE,0
Pstres,on
allsel
nsel,s,loc,y,0               
!施加约束
d,all,all,0
!!由下向上施加力载
!!由下向上施加力载
allsel
nsel,s,loc,y,lh            !施加约束
d,all,ux,0
d,all,uz,0
d,all,rotx,0
d,all,roty,0
d,all,rotz,0
allsel
Sfl,6,Pres,60000
Sfl,8,Pres,60000
Sfl,10,Pres,60000
Sfl,12,Pres,60000
Allsel
solve
finish
/solu
Antype,buckle                        !  指定屈曲分析
Bucopt,lanb,1
Mxpand,1                            !  模态扩展
Outres,all,all
Solve
finish
!线性屈曲分析完成
/post1
*get,load1,mode,1,freq
pi=4*atan(1)
bucklingload=load1*pi*(0.5015**2-0.4985**2)*60000
nsort,u,sum
*get,dmax,sort,,max
*set,fscale,0.00015/ dmax
finish
!接着进行非线性屈曲分析
/solu
upcoord, fscale,off
!upcoord,0.2,off  !指定结构缺陷因子为0.05
antype,0   !选择静态分析选项
eqslv,pcg,1E-4  !选择pcg预条件共轭梯度迭代方程求解器,精度1e-4
lnsrch,1        !打开线性搜索技术
neqit,15   !指定平衡迭代次数最大为15，在默认情况下为25次太耗时间
nropt,full  !应用newton-raphson选项
nlgeom,on  !打开大变形选项
sfscale,pres,load1*0.5*60000 !放大结构上施加的载荷为预测值的0.7倍
time,1          !指定时间为1
autots,1         !打开自动时间步长
nsubst,200   !指定子步数为100
kbc,0        !逐渐递增载荷
outres,,1    !指定每1个子步输出一次结果
lswrite,1    !写入第1载荷步文件
lsread,1
lssolve,1,1        !读入1载荷步并求解
autots,1
nsubst,50000
Sfl,6,Pres, load1*1.2*60000
Sfl,8,Pres, load1*1.2*60000
Sfl,10,Pres, load1*1.2*60000
Sfl,12,Pres, load1*1.2*60000  !在模型顶端施加压力
ARCLEN,ON,               ! 打开弧长法，并设置最大乘子系数为1
lswrite,2            !写入第2载荷步文件
lsread,2
!ARCTRM,L
lssolve,2,2        !读入2载荷步并求解
finish
/post26
!*get,nmax,node,,num,max  !选取最大节点数
rforce,2,1,f,y
*do,i,2,41    !计算反作用力之和,其中变量2表示反作用力总和
rforce,3,i,f,y
add,2,3,2
*enddo
*do,i,1518,1557    !计算反作用力之和,其中变量2表示反作用力总和
rforce,4,i,f,y
add,2,4,2
*enddo
*do,i,2998,3037   !计算反作用力之和,其中变量2表示反作用力总和
rforce,5,i,f,y
add,2,5,2
*enddo
*do,i,4478,4516    !计算反作用力之和,其中变量2表示反作用力总和
rforce,6,i,f,y
add,2,6,2
*enddo
prod,5,4,,,,,,0.001
!反作用力计算完毕

nsol,4,3075,u,y  !提取节点18209的y向位移赋予变量4
prod,5,4,,,,,,-1  !改变位移值的符号，因为求得的位移值为负改变符号后变为正值以利于结果处理
!prvar,2,3,4 !列表显示变量2,3,4
/axlab,x,deflection(M)  !定义x轴名称
/axlab,y, TOTAL LOAD(N)  !定义y轴名称
!/xrange,0,0.5E-3          !指定x轴范围
!/yrange,-0.5*bucklingload,1.5*bucklingload  !指定y轴范围
/grid,1 !绘制x、y 栅格
xvar,5                    !将变量4指定为x轴变量
plvar,2                  !绘制变量2
finish