以下是我的命令流，在面9上施加位移荷载，小于1时都能收敛，但是大于1时就无法收敛了，麻烦各位帮我看一下，本人感激不尽！
出现的警告主要有一下三种：
（1）There are 20 small equation solver pivot terms.
（2）Contact element 3610 (real ID 4) status changes abruptly from contact (with target element 3531) -> no-contact.
（3）One or more elements have become highly distorted. Excessive distortion of elements is usually a symptom indicating the need for corrective action elsewhere. Try incrementing the load more slowly (increase the number of substeps or decrease the time step size). You may need to improve your mesh to obtain elements with better aspect ratios. Also consider the behavior of materials, contact pairs, and/or constraint equations. If this message appears in the first iteration of first substep, be sure to perform element shape checking.
a=5
b=5
c=30
m=1
deep=1.5
r=0.5
t=0.015
weiyi=0.5
pi=3.14
angle=60
y=b/2-deep
duanju=1
!单元定义
/prep7
*AFUN,DEG
ET,1,SHELL281
sect,1,shell,,
secdata, 0.015,1,0.0,3
secoffset,MID
seccontrol,,,, , , ,
ET,2,SOLID45
!材料属性
!管道
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,E1
MPDATA,PRXY,1,,Po1
TB,MISO,1,1,3,0
TBTEMP,0
TBPT,,0.0022,462e6
TBPT,,0.040,517e6
TBPT,,0.044,517e6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7850
!砂粘土
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,8e6
MPDATA,PRXY,2,,0.37
TB,DP,2,,,
TBMODIF,1,1,5e6
TBMODIF,1,2,20
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,1900
!石灰岩
MPDATA,EX,3,,43.96e6
MPDATA,PRXY,3,,0.23
TB,DP,3,,,
TBMODIF,1,1,22.79e6
TBMODIF,1,2,44.1
TBMODIF,1,3,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,2710
!创建点
k,1,-a/2,-b/2,0
k,2,a/2,-b/2,0
k,3,a/2,b/2,0
k,4,-a/2,b/2,0
k,5,-a/2,-b/2,c/2-m/2+b/2/TAN(angle)
k,6,a/2,-b/2,c/2-m/2+b/2/TAN(angle)
k,7,a/2,b/2,c/2-m/2-b/2/TAN(angle)
k,8,-a/2,b/2,c/2-m/2-b/2/TAN(angle)
k,9,-a/2,-b/2,c/2+m/2+b/2/TAN(angle)
k,10,a/2,-b/2,c/2+m/2+b/2/TAN(angle)
k,11,a/2,b/2,c/2+m/2-b/2/TAN(angle)
k,12,-a/2,b/2,c/2+m/2-b/2/TAN(angle)
k,13,a/2,-b/2,c
k,14,-a/2,-b/2,c
k,15,-a/2,b/2,c
k,16,a/2,b/2,c
!创建面
a,1,2,3,4
a,4,1,5,8
a,2,3,7,6
a,3,4,8,7
a,5,6,2,1
a,5,6,7,8
a,12,11,10,9
a,11,16,13,10
a,11,12,15,16
a,14,15,12,9
a,13,14,15,16
a,9,10,13,14
!创建体
va,1,2,3,4,5,6
va,7,8,9,10,11,12
!挖空
cyl4,0,y,r,,,,c
FLST,2,2,6,ORDE,2
FITEM,2,1
FITEM,2,-2
VSBV,P51X, 3
!创建管道
cyl4,0,y,r-t,,r,,c
TYPE, 1
MAT, 1
REAL,
ESYS, 0
SECNUM,
!*
lesize,45,,,4,,,,,1
lesize,46,,,4,,,,,1
lesize,47,,,4,,,,,1
lesize,48,,,4,,,,,1
lesize,53,,,4,,,,,1
lesize,54,,,4,,,,,1
lesize,55,,,4,,,,,1
lesize,56,,,4,,,,,1
lesize,33,,,4,,,,,1
lesize,34,,,4,,,,,1
lesize,43,,,4,,,,,1
lesize,44,,,4,,,,,1
lesize,49,,,4,,,,,1
lesize,50,,,4,,,,,1
lesize,51,,,4,,,,,1
lesize,52,,,4,,,,,1
lesize,35,0.25,,,,,,,1
lesize,36,0.25,,,,,,,1
lesize,57,0.25,,,,,,,1
lesize,58,0.25,,,,,,,1
lesize,59,0.25,,,,,,,1
lesize,60,0.25,,,,,,,1
lesize,39,0.25,,,,,,,1
lesize,40,0.25,,,,,,,1
SMRT,6
AMESH,7
AMESH,11
!
TYPE, 2
MAT, 2
REAL,
ESYS, 0
SECNUM,
!*
lesize,29,,,4,,,,,1
lesize,30,,,4,,,,,1
lesize,31,,,4,,,,,1
lesize,32,,,4,,,,,1
lesize,25,,,4,,,,,1
lesize,26,,,4,,,,,1
lesize,27,,,4,,,,,1
lesize,28,,,4,,,,,1
lesize,18,0.25,,,,,,,1
lesize,21,0.25,,,,,,,1
lesize,22,0.25,,,,,,,1
lesize,24,0.25,,,,,,,1
lesize,17,0.25,,,,,,,1
lesize,19,0.25,,,,,,,1
lesize,20,0.25,,,,,,,1
lesize,23,0.25,,,,,,,1
VSWEEP,5
!*
TYPE, 2
MAT, 3
REAL,
ESYS, 0
SECNUM,
!*
lesize,1,0.25,,,,,,,1
lesize,2,0.25,,,,,,,1
lesize,3,0.25,,,,,,,1
lesize,4,0.25,,,,,,,1
lesize,5,0.25,,,,,,,1
lesize,7,0.25,,,,,,,1
lesize,8,0.25,,,,,,,1
lesize,10,0.25,,,,,,,1
VSWEEP,4
!创建接触对
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,3,0.3
MAT,3
MP,EMIS,3,7.88860905221e-031
R,3
REAL,3
ET,3,170
ET,4,174
R,3,,,0.1,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0
KEYOPT,4,10,2
KEYOPT,4,11,0
KEYOPT,4,12,0
KEYOPT,4,2,0
KEYOPT,3,5,0
! Generate the target surface
ASEL,S,,,7
CM,_TARGET,AREA
TYPE,3
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
NSLE,A,CT2
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,18
CM,_CONTACT,AREA
TYPE,4
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,3
ESEL,A,TYPE,,4
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,3
ESEL,A,TYPE,,4
ESEL,R,REAL,,3
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
!*
!*
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,3,0.3
MAT,3
MP,EMIS,3,7.88860905221e-031
R,4
REAL,4
ET,5,170
ET,6,174
R,4,,,0.1,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,6,4,0
KEYOPT,6,5,0
KEYOPT,6,7,0
KEYOPT,6,8,0
KEYOPT,6,9,0
KEYOPT,6,10,2
KEYOPT,6,11,0
KEYOPT,6,12,0
KEYOPT,6,2,0
KEYOPT,5,5,0
! Generate the target surface
ASEL,S,,,7
CM,_TARGET,AREA
TYPE,5
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
NSLE,A,CT2
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,21
CM,_CONTACT,AREA
TYPE,6
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,5
ESEL,A,TYPE,,6
ESEL,R,REAL,,4
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,5
ESEL,A,TYPE,,6
ESEL,R,REAL,,4
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
!*
!*
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,2,0.7
MAT,2
MP,EMIS,2,7.88860905221e-031
R,5
REAL,5
ET,7,170
ET,8,174
R,5,,,0.1,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,8,4,0
KEYOPT,8,5,0
KEYOPT,8,7,0
KEYOPT,8,8,0
KEYOPT,8,9,0
KEYOPT,8,10,2
KEYOPT,8,11,0
KEYOPT,8,12,0
KEYOPT,8,2,0
KEYOPT,7,5,0
! Generate the target surface
ASEL,S,,,11
CM,_TARGET,AREA
TYPE,7
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
NSLE,A,CT2
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,17
CM,_CONTACT,AREA
TYPE,8
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,7
ESEL,A,TYPE,,8
ESEL,R,REAL,,5
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,7
ESEL,A,TYPE,,8
ESEL,R,REAL,,5
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
!*
!*
/COM, CONTACT PAIR CREATION - START
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,2,0.7
MAT,2
MP,EMIS,2,7.88860905221e-031
R,6
REAL,6
ET,9,170
ET,10,174
R,6,,,0.1,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,10,4,0
KEYOPT,10,5,0
KEYOPT,10,7,0
KEYOPT,10,8,0
KEYOPT,10,9,0
KEYOPT,10,10,2
KEYOPT,10,11,0
KEYOPT,10,12,0
KEYOPT,10,2,0
KEYOPT,9,5,0
! Generate the target surface
ASEL,S,,,11
CM,_TARGET,AREA
TYPE,9
NSLA,S,1
ESLN,S,0
ESLL,U
ESEL,U,ENAME,,188,189
NSLE,A,CT2
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
ASEL,S,,,22
CM,_CONTACT,AREA
TYPE,10
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008)
ESURF
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,9
ESEL,A,TYPE,,10
ESEL,R,REAL,,6
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,9
ESEL,A,TYPE,,10
ESEL,R,REAL,,6
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/COM, CONTACT PAIR CREATION - END
FINISH
/SOL
ANTYPE,0 !静力分析
NLGEOM,1 !开启大变形
AUTOTS,1 !开启自动时间步长
TIME,1 !载荷步结束时间为1
LNSRCH,1 !开启线性搜索
OUTRES,ERASE
OUTRES,ALL,LAST !将最后一个子步写入结果文件
RESCONTRL,DEFINE,ALL,LAST,1 !重新启动时写入最后一个子步
SOLCONTROL,ON,0,
NSUBST,20,100,1 !子步50，最大子步5000最小子步50
NEQIT,50 !平衡迭代最大次数为50
PRED,ON !非线性分析中第一子步后打开预测器
NROPT,FULL, ,ON !完全牛顿拉夫逊法
LUMPM,1
EQSLV, , ,0,
PRECISION,0
MSAVE,0
PIVCHECK,1
SSTIF,ON
TOFFST,0,
CNVTOL, U, , 0.5, , 1
!*
DA,5,ALL,
DA,8,UX,
DA,10,UX
DA,9,UY,1
SOLVE
