求教!热力学仿真中出现的问题！
本帖最后由 zou 于 2010-12-1 08:57 编辑做的是一个芯片的热力学仿真模型
想观察的结果是芯片在温度持续大范围变化下所受应力，变形，应变及应力最大点的剪切应力随时间变化等情况
代码如下
前面建模部分省略了 太多太长了 如果有需要的告诉我 我再放上去
哦 忘记把问题抛出来了~
首先是求解过程中出现的三个错误或警告信息（出现了N多）：
1．
WARNING:
There are 1 small equation solver pivot terms.2．
WARNING:The unconverged solution (identified as time 1800 substep 999999) is output for analysis debug purposes. Results should not be used for any other purpose.3. ERROR:Solution not converged at time 105 (load step 1 substep 1). Run terminated.
然后是解决结果不能顺利查看
命令如下：
1。查看总体变形分布图 PLNSOL,U,SUM,0,1
2。查看等效应力分布图 PLBSOL,S,EQV,0,1
3。查看等效应变分布图 PLNSOL,EPTO,EQV,0,1
以上三条命令显示结果均没有变化
4。进入时间历程后处理ANSOL,2,1103，S,XY,SXY_2
XVAR,1
PLVAR,2,
提示没有数据点被存储在这一存储单元步中。
脚本如下：（只放入划分至求解部分）
!进行自由网格划分
MSHAPE,0,2D
MSHKEY,0
!*
FLST,5,7,5,ORDE,4
FITEM,5,1
FITEM,5,-5
FITEM,5,13
FITEM,5,20
CM,_Y,AREA
ASEL, , , ,P51X
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
FINISH
SAVE
/SOLU
!设定为瞬态分析
ANTYPE,4
!设定为大变形分析
NLGEOM,1
!施加对称边界条件
LSEL,S,LOC,X,0
DL,ALL, ,SYMM
NSEL,S,LOC,X,0
NSEL,R,LOC,Y,0
D,ALL,,0,,,,UY,,,,,
ALLSEL,ALL
TREF,25
!设定摄氏温度与开尔文温度之间的差值
TOFFSET,273
!设定只有一个时间很小的载荷步，作为初始条件
TUNIF,25
TIME,1E-5
!设定非线性分析的收敛值
CNVTOL,F,,,,1e-5 ----很多帖子介绍上面的错误2可以用这句约束来解决，可是加了以后仍然会报错
SOLVE
!设定子步数为20
NSUBST,6
!输出每个子步的所有值
OUTRES,ALL,ALL
!升温阶段
TUNIF,125
TIME,5*60
SOLVE
!保温阶段
/SOLU ----求解部分主要参考了一本教材，但是这一句原来是没有的，可不加上去后面的TIME,SOLVE指令就提示是不合法的指令了，感觉整个时间历程被分别执行了，而不是一个整体过程了，个人感觉是这样
TIME,30*60
SOLVE
!降温阶段
TUNIF,-55
TIME,35*60
SOLVE
!保温阶段
TIME,60*60
SOLVE
………………
(中间均是升降温之间的变化过程)
!降温阶段
TUNIF,-55
TIME,215*60
SOLVE
!保温阶段
TIME,240*60
SOLVE
!升温阶段
TUNIF,125
TIME,245*60
SOLVE
!保温阶段
TIME,270*60
SOLVE
!回到初始温度阶段
TUNIF,25
TIME,275*60
SOLVE
SAVE
FINISH
从提示信息上来看应该是在求解的过程中系数矩阵中出现小主元，最终导致结果没有收敛，主要原因可能是：约束不够，有刚体位移，或者是网格较差，网格扭曲较为严重，你可以从这两个方面找找原因
{:{23}:}
问题查了下 还是约束不足的问题 已经基本解决啦 谢谢啊
页:
[1]