麻烦各位帮我看看为什么计算结果不收敛，这个模型是在模拟拉索锚固区受力，很简单的模型，最后总是提醒：

 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.  

不知如何解决，命令流如下，需要导入截面，无法上传还请见谅。

FINI

/CLEAR,START

/AUX15

~SATIN,'small','sat',,ALL,0   !从CAD输出的.sat文件导入小截面small

/NOPR 

/PREP7

!0. 定义变量(单位采用N-mm制)------------------------------

!注：不设置密度DENS，即不考虑重力影响

EMCI=2.1e5         !锚杯铸钢ZG20Mn弹性模量 ELASTICITY MODULOUS OF CAST 

IRON

EMP=1.5e5          !热铸锌铜合金弹性模量 ELASTICITY MODULOUS OF PLATINE

!1. 定义单元、材料属性及截面---------------------------

ET,1,SOLID187       !铸体、锚杯

       

MP,EX,1,EMCI        !锚杯材料特性，考虑材料非线性，采用多线性等向强化模

型MISO

MP,PRXY,1,0.3

TB,MISO,1,,3

TBPT,,0.001905,400

TBPT,,0.003,401

TBPT,,0.021,625

MP,EX,2,EMP         !铸体材料特性，考虑材料非线性，采用多线性等向强化模

型MISO

MP,PRXY,2,0.27

TB,MISO,2,,3

TBPT,,0.002667,400

TBPT,,0.003,401

TBPT,,0.021,625

R,1                  !锚杯实常数

R,2                  !铸体实常数

!2.1 铸体几何模型建立----------------------------

AGLUE,ALL

LGLUE,ALL

NUMMRG,ALL

NUMCMP,ALL

ASEL,S,LOC,Z,0

CM,DY1,AREA           !定义组件“定位1”

PCIRC,36.5

ASEL,S,LOC,Z,0

ASEL,R,LOC,Y,0

ASEL,R,LOC,X,0

CM,LK1,AREA           !定义组件“轮廓1”

ALLSEL,ALL

ASBA,LK1,DY1,0,DELETE,KEEP

VEXT,ALL,,,,,435,2.386767,2.386767

VATT,2,2,1            !分别表示材料号、实常数号、单元类型号

!2.2 锚杯几何模型建立----------------------------

NUMSTR,KP,1001

NUMSTR,LINE,1001

NUMSTR,AREA,1001

NUMSTR,VOLU,1001

CSYS,0

K, 1001  ,0,83.000,495.000,

K, 1002  ,0,83.000,511.000,

K, 1003  ,0,83.000,516.000,

K, 1004  ,0,87.000,520.000,

K, 1005  ,0,120.000,520.000,

K, 1006  ,0,125.000,515.000,

K, 1007  ,0,125.000,511.000,

K, 1008  ,0,125.000,495.000,

K, 1009  ,0,125.000,436.288,

K, 1010  ,0,79.306,115.643,

K, 1011  ,0,118.905,70.000,

K, 1012  ,0,125.000,70.000,

K, 1013  ,0,125.000,2.000,

K, 1014  ,0,123.000,.000,

L,260,1001

L,1001,1002

L,1002,1003

L,1003,1004

L,1004,1005

L,1005,1006

L,1006,1007

L,1007,1008

L,1008,1009

L,1009,1010

L,1010,1011

L,1011,1012

L,1012,1013

L,1013,1014

L,1014,129

LSEL,S,LINE,,1001,1015

LSEL,A,LINE,,745

AL,ALL

VROTAT,1001,,,,,,33,138,,4      !面绕着旋转轴生成体

ALLS

NUMMRG,ALL

VSEL,S,,,1001,1004

VATT,1,1,1                      !分别表示材料号、实常数号、单元类型号

!3. 网格划分----------------------------

LSEL,S,LOC,Z,1,434

CSWPLA,11,1

LSEL,R,LOC,X,-60,60

CSYS,0

LESIZE,ALL,,,40

ALLS

LSEL,S,,,1,1000

LSEL,U,LOC,Z,100,300

LSEL,R,LENGTH,,1,20

LESIZE,ALL,,,1  

ALLS

LSEL,S,,,1,1000

LSEL,U,LOC,Z,100,300

LSEL,U,LENGTH,,1,20

LESIZE,ALL,,,12

ALLS

MSHAPE,1,3D

MSHKEY,0

SMRTSIZE,3

VMESH,ALL

!*GET,AA,ELEM,0,COUNT

!*STATUS,AA                     !查询单元数量

!4. 施加荷载及边界条件----------------------------

ALLS

NSEL,S,LOC,Z,0

CSWPLA,11,1

NSEL,R,LOC,X,100,10000

CSYS,0

D,ALL,ALL

ALLS

LSEL,S,LOC,Z,0

LSEL,R,LENGTH,,1,10

NSLL,S,1

F,ALL,FZ,-10000

ALLSEL,ALL

FINI

/SOLU

CNVTOL,U,,0.05,2              !位移的收敛准则为0.05

CNVTOL,F,,0.01,2,1            !位移的收敛准则为0.01

ANTYPE,0                      !定义分析类型：静态分析

NLGEOM,ON                     !几何非线性开关

TIME,1                        !荷载步结束时的时间值

NSUBST,10                     !荷载步的子步数

OUTRES,ALL,ALL

AUTOTS,ON                     !自动时间步

NEQIT,50                      !平衡迭代的最大次数

SOLVE

FINISH

!5. 后处理-------------------------

/POST1

SET,LAST                      !查看某个时刻的计算结果

PLDISP,2                      !显示结构变形图

PLESOL,S,X,0                  !显示单元X方向应力

PLESOL,S,Y,0                  !显示单元Y方向应力

PLESOL,S,Z,0                  !显示单元Z方向应力

FINI