ansys非线性分析不收敛,求大神解答?
麻烦各位帮我看看为什么计算结果不收敛,这个模型是在模拟拉索锚固区受力,很简单的模型,最后总是提醒:
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