非线性屈曲的例子

这是我做的一个非线性屈曲的例子，是双套管的耗能支撑，请问接下来我如何做才能得到滞回曲线呢？
哪位大侠可以帮帮我啊

finish
/units,SI
/clear
/cwd,'E:\Ansys\PIPE'
/filname,PIPE0
/title,PIPE0
/prep7

ET,1,SOLID45     !外套管及连接段采用单元
ET,2,SOLID185     !芯材采用单元
ET,3,MESH200     !定义mesh单元
KEYOPT,3,1,6
KEYOPT,3,2,0

!*  定义材料属性
!*  Q235钢
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.06E11
MPDATA,PRXY,1,,0.3
TB,BKIN,1,1,2,1
TBTEMP,0
TBDATA,,235e6,0,,,,
MPDATA,DENS,1,,7800
!*  Q345钢
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,2.06E11
MPDATA,PRXY,2,,0.3
TB,BKIN,2,1,2,1
TBTEMP,0
TBDATA,,390e6,0,,,,
MPDATA,DENS,2,,7800

!定义尺寸参数
*set,R1,0.037      !定义钢管内径
*set,R2,0.045      !定义钢管外径
*set,R3,0.050      !定义钢管外径
*set,L1,3.0        !定义内钢管长度
*set,L2,2.4       !定义外钢管长度
*set,L3,0.5       !定义连接段长度

!!!!!!!!!!!!!!!!
!               建立核心段                   !
!!!!!!!!!!!!!!!!

CYLIND,R1,R2,-L1/2,-L1/2+L3
CYLIND,R1,R2,-L1/2+L3,L1/2-L3
CYLIND,R1,R2,L1/2-L3,L1/2
ALLSEL
NUMMRG,ALL, , , ,LOW
NUMCMP,ALL
!!!!!!!!!!!!!!!!
!               建立外包段                   !
!!!!!!!!!!!!!!!!

CYLIND,R2,R3,-L2/2,L2/2

!网格精度
*set,EDIV,5
/PREP7
!建立核心段
LSEL,S,LOC,Z,-L1/2+L3
!LESIZE,ALL,EDIV, , , , , , ,1
LESIZE,ALL, , ,EDIV, , , , ,1
LSEL,ALL
mshape,0
mshkey,2
TYPE,3
MAT,1
AMESH,2

TYPE,1
MAT,2
esize,0.02
VSWEEP,1

TYPE,2
MAT,1
esize,0.02
VSWEEP,2

TYPE,1
MAT,2
esize,0.02
VSWEEP,3

!建立外包钢管段
LSEL,S,LOC,Z,-L2/2
!!LESIZE,ALL,EDIV, , , , , , ,1
LESIZE,ALL, , ,EDIV, , , , ,1
LSEL,ALL
mshape,0
mshkey,2
TYPE,3
MAT,1
AMESH,17

TYPE,1
MAT,2
esize,0.02
VSWEEP,4

!建立连接单元 contac target
ET,4,CONTA174    !定义接触单元
ET,5,TARGE170    !定义目标单元
R,10
REAL,10
ALLSEL
ASEL,S,AREA,,8
ASEL,A,AREA,,9
type,4
NSLA,S,1
ESLN,S,0
esurf,all,
ALLSEL

R,10
REAL,10
ALLSEL
ASEL,S,AREA,,21
ASEL,A,AREA,,22
type,5
NSLA,S,1
ESLN,S,0
esurf,all,
ALLSEL

SAVE,'PIPE','DB'

 

!!*  静力求解
/SOLU
ALLSEL,ALL
ANTYPE,static
! NLGEOM,ON
LSEL,S,LOC,Z,L1/2
NSLL,S
D,ALL,UX,0
D,ALL,UY,0
D,ALL,UZ,0
allsel

LSEL,S,LOC,Z,-L1/2
NSLL,S
D,ALL,UX,0
D,ALL,UY,0
allsel

!施加外包钢管的约束
LSEL,S,LOC,Z,L2/2
LSEL,A,LOC,Z,-L2/2
LSEL,S,RADIUS,,R3
NSLL,S
D,ALL,ALL
allsel

!施加端部作用力
LSEL,S,LOC,Z,-L1/2
NSLL,S
F,ALL,FZ,1/16
ALLSEL
pstres,on
outres,all,all
ALLSEL,ALL
SOLVE
FINISH
!
!!*  特征值屈曲求解
/SOLU
ALLSEL,ALL
ANTYPE,BUCKLE
BUCOPT,LANB,10,0,0
MXPAND,10,0,0,1,0.001
SOLVE
FINISH

/post1
set,first
*get,freq,mode,1,freq
PLDISP,1
FINISH
/PREP7
UPGEOM,L1/1000/1.001,,,'PIPE0','rst' !施加1%的初始缺陷，以便进行几何非线性分析，upgeom只能在/prep7中使用
FINISH