声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3477|回复: 2

[CFD及热分析] 求教!热力学仿真中出现的问题!

[复制链接]
发表于 2010-12-1 08:56 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 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
回复
分享到:

使用道具 举报

发表于 2010-12-5 09:02 | 显示全部楼层
从提示信息上来看应该是在求解的过程中系数矩阵中出现小主元,最终导致结果没有收敛,主要原因可能是:约束不够,有刚体位移,或者是网格较差,网格扭曲较为严重,你可以从这两个方面找找原因

 楼主| 发表于 2010-12-5 20:22 | 显示全部楼层

问题查了下 还是约束不足的问题 已经基本解决啦 谢谢啊
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-30 02:14 , Processed in 0.057398 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表