ansys结构分析不收敛，请高手指点，附详细命令流

a130676

先热分析，求解完后直接在命令栏输入结构分析命令。热分析结果文件不是rth，是rst，之后结构分析不收敛，请高手指点
FINI                                       !
/CLE
/FILNAME,1
/TITLE,refenxi
/UNITS,MKS
/PREP7
ET,1,PLANE77,,,1                                              !定义单元类型并设置轴对称
ET,2,PLANE183,,,1
MPTEMP,1,0,50,100,200,400,600                                 !定义材料表格数据温度点
MPTEMP,7,700
MPDATA,C,1,1,610,680,783,826,967,1072                        !定义材料表格数据比热容
MPDATA,C,1,7,1079
MPTEMP,,,,,,,,                                                !定义材料表格导热系数
MPTEMP,1,20,100,200,400,600,700
MPDATA,KXX,1,1,1.02,1.14,1.37,1.52,1.97,2.79
MPTEMP,,,,,,,,                                                !定义材料表格热膨胀系数
MPTEMP,1,0,100,200,300,400,500
MPTEMP,7,600,700
MPDATA,ALPX,1,1,2.6e-7,4.35e-7,5.12e-7,5.98e-7,5.14e-7,5.01e-7
MPDATA,ALPX,1,7,4.93e-7,
UIMP,1,REFT,,,20                                              !定义材料热膨胀参考温度
MPTEMP,,,,,,,,
MPTEMP,1,20,200,300,700                            !定义材料表格弹性模量和泊松比
MPDATA,EX,1,1,6.69e10,6.78e10,7.4e10,7.4e10,
MPDATA,PRXY,1,1,0.14,0.14,0.14,0.14
MP,DENS,1,1923                                                !定义材料密度
TB,PRONY,1,0,4,SHEAR                                          !定义粘性参数（剪切）viscosity parameters（shear）
TBDATA,1,0.0081426,0.10927,0.017958,2.7854,0.017962,3.123
TBDATA,7,0.017962,3.1237
TB,SHIFT,1,0,1,2                                              !定义shift function
TBDATA,1,1000,4963.3
BLC4,0,0,0.00004,0.015                                        !开始建模
/PNUM,LINE,1                                                  !显示所有线号
APLOT                                                         !显示面
/PNUM,AREA,1                                                  !显示面号
!开始划分网格
MAT,1                                                         
TYPE,1
ESIZE,0.00001
AMESH,1
!进入求解器
/SOL
ANTYPE,4                                                     !激活瞬态分析
NLGEOM,1                                                     !激活大变形
TRNOPT,FULL                                                  !采用完全法
!定义边界条件
NSEL,ALL                                                     !选择X轴对称约束节点
NSEL,S,LOC,Y,0                                                                              
DSYM,SYMM,Y  
NSEL,ALL                                                   
NSEL,S,LOC,X,0
DSYM,SYMM  
NSEL,ALL
NSEL,S,LOC,X,0.00004
NSEL,R,LOC,Y,0,0.005
D,ALL,TEMP,600                                               !施加温度载荷
NSEL,ALL                                                     !选择加载温度载荷节点            
NSEL,S,LOC,X,0.00004
NSEL,R,LOC,Y,0.005,0.015
SF,ALL,CONV,7.5,20                                           !施加对流载荷
TREF,20                                                      !定义参考温度
TUNIF,20                                                     !定义整体温度
KBC,1                                                        !加载方式为阶跃载荷
DELTIM,0.01                                                  !最小时间步长0.001s
TIME,200                                                     !最终时间200s
NSUBST,400,800,200                                           !求解设置子步、最大最小子步数
OUTRES,,ALL                                                  !设置结果输出
AUTOTS,OFF                                                   !关闭自动时间步长
NSEL,ALL
SOLVE
!后处理
/POST1  
/EFACET,1   
PLNSOL, TEMP,, 0,1.0                                         !显示节点温度
求解完后直接在命令栏输入以下命令进行结构分析，为什么不收敛？将迭代次数改到1000还不行，麻烦高手帮忙看看
/title,jiegoufenxi
/prep7
nsel,all
lsclear,all
ddel,all
etchg,tts
keyopt,1,3,1
keyopt,1,6,1
ANTYPE,4                                                     !激活瞬态分析
NLGEOM,1                                                     !激活大变形
TRNOPT,FULL                                                  !采用完全法
ldread,temp,,,,,'1','rst'
NSEL,ALL                                                     !选择X轴对称约束节点
NSEL,S,LOC,Y,0                                                                              
DSYM,SYMM,Y  
NSEL,ALL                                                   
NSEL,S,LOC,X,0
DSYM,SYMM  
nsel,s,loc,y,0
d,all,uy,0
nsel,s,loc,y,0.015
d,all,uy,12
tunif,20
autots,on
deltim,0.1,0.01,1
time,400
kbc,0
NEQIT,50
/solu
nsel,all
solve

htbbzzg

可能有几个问题：
1  ET,2,PLANE183,,,1    ! 在热分析部分定义了结构单元,造成结果文件扩展名为 .rst，删除该命令，结果文件就是 .rth 了。
2  不知是什么材料，导热系数太小了，是非金属材料吗？
3  热载荷为点载荷，会造成局部温度梯度过大，应力分析时会造成局部相对变形过大，单元高度扭曲。
4  结构载荷有致命问题：
    nsel,s,loc,y,0.015
    d,all,uy,12
    !  Uy 值太大了 - 圆柱长度只有 0.015，位移载荷是其  12/0.015 = 800 倍，更会造成变形后单元高度扭曲
可能还有其它问题，自己再试试。

a130676

htbbzzg 发表于 2014-10-27 11:35
可能有几个问题：
1  ET,2,PLANE183,,,1    ! 在热分析部分定义了结构单元,造成结果文件扩展名为 .rst，删 ...

受教了，谢谢。我去调调试试

369

{:{39}:}

a130676

htbbzzg 发表于 2014-10-27 11:35
可能有几个问题：
1  ET,2,PLANE183,,,1    ! 在热分析部分定义了结构单元,造成结果文件扩展名为 .rst，删 ...

nsel,s,loc,y,0.015   d,all,uy,12 将位移12改为0.012，deltim 0.01,0.00001,1后在164.72444s时不收敛，此时位移0.00492，是位移太大了吗

htbbzzg

    d,all,uy,12 改为 d,all,uy,0.012，与圆柱长度 0.015 相比，仍然有点大，但是计算结果位移 0.00492 (是 Uy 吗？) 好像还可以。
    不过，前面提出过单点温度载荷造成温度梯度过大的问题，不知你是如何处理的？过大的温度梯度会造成局部相对变形较大，也可能造成不收敛的。
    另外，不收敛时，仍可以到后处理看一下已算出的结果，其中不收敛的载荷步的结果以 999999 标记。看一下不收敛的结果，特别是应力结果，也有助于判断不收敛的原因。

a130676

htbbzzg 发表于 2014-10-30 18:15
d,all,uy,12 改为 d,all,uy,0.012，与圆柱长度 0.015 相比，仍然有点大，但是计算结果位移 0.00492 (是 ...

在加热范围相对总体坐标系固定不变且拉伸时仍在加热的情况下除了通过节点加载外还有其他加载方式吗？局部变形过大甚至扭曲这个缺陷还没找到解决办法，请再指点指点。另外程序因为不收敛而停止求解，怎么才能查看已求解的那部分时长的结果？

htbbzzg

    一般的载荷，除了节点加载外，也可以加载到单元或实体，要看具体是什么载荷。
    计算不收敛，一般会把最后不收敛的结果存为结果组 999999，可以到一般后处理去查看 (包括已完成并存储的子步的结果)。