fbmjing 发表于 2006-10-2 09:34:09

是加载还是计算（收敛）方面出现问题了？（命令流）

最近在做一个联箱方面的题，主要是看它的管接头处的应力分布。思路如下：
先建模（由于是电厂过热器出口联箱，有很多孔，我只拿出其中一个的一小部分研究）
之后加载定义好的 时间－压力、时间－内表面温度、时间－外表面温度表格，
最后分网计算。
   问题就在求解上，求解老是出错，提示说”Solution not converged at time 600(load step 1 substep 1)“.应该是收敛方面出现问题了。(求解时间挺长的，大概要20分钟)
说实话，我对求解方面学习很浅，很多都不懂。自己的程序都不是很明白。请大家不吝赐教。


/Prep7
ET,1,SOLID98         !定义单元类型

TOFFST,273         

!定义材料参数
mp,c,1,1e-6
mp,dens,1,7.28e-6
mptemp,1,300,400,500,600,700
mpdata,kxx,1,1,0.0292,0.0291,0.029,0.0293,0.029
mpdata,alpx,1,1,11.7e-6,11.9e-6,12.2e-6,12.4e-6,12.6e-6
mpdata,ex,1,1,2.04e11,1.95e11,1.85e11,1.70e11,1.58e11
mpdata,prxy,1,1,0.3,0.3,0.3,0.3,0.3

cylind,0.31075,0.42725,-0.3,0.3,-9,9
wprota,0,0,90
cylind,0.014,0.021,0,0.537,0,360
vovlap,1,2
VDELE,3, , ,1
VDELE,5, , ,1
numcmp,all                  !压缩编号


!形成内表面集合
ASEL,S,AREA,,6,7,1,
ASEL,A,AREA,,11,11
ASEL,A,AREA,,14,16,1,
CM,IN_AREA,AREA

!形成外表面集合
ASEL,S,AREA,,12,13,1,
ASEL,A,AREA,,17,17
CM,OUT_AREA,AREA

!定义时间-压力表格，（表在附件里）
*DIM,TABPRESS,TABLE,53,1,1,TIME,PRES,
*tread,TABPRESS,p3,txt

!定义时间-内壁温度表格
*DIM,TABTEMPIN,TABLE,53,1,1,TIME,TEMP,
*tread,TABTEMPIN,t_in,txt

!定义时间-外壁温度表格
*DIM,TABTEMPOUT,TABLE,53,1,1,TIME,TEMP,
*tread,TABTEMPOUT,t_out,txt

!在内表面加载压力表格
SFA,IN_AREA,1,PRES, %TABPRESS%

!在内壁加载温度表格
DA,IN_AREA,TEMP, %TABTEMPIN% ,1

!在外壁加载温度表格
DA,OUT_AREA,TEMP, %TABTEMPOUT% ,1


!画网格
ALLSEL,ALL
VATT,1,1,,0,

VSEL,ALL
smrt,5
mshape,1,3d
mshkey,0
vmesh,all

ALLSEL,ALL

save
finish

!计算
/SOLU
ANTYPE,4
TRNOPT,FULL
LUMPM,0

solcontrol,on
autots,off
deltim,600,600,1200,on
time,31200

OUTPRE,basic,-100
outres,all,-100
kbc,0

solve

finish



所需的时间－压力、时间－内表面温度、时间－外表面温度表格都是以每600秒取一次值，可参见附件。

[ 本帖最后由 fbmjing 于 2006-10-2 09:27 编辑 ]

fanmingjian 发表于 2006-10-25 16:38:29

我也遇到同样的问题了，很郁闷，一直没有解决！

diaohuwang 发表于 2006-11-14 12:13:50

我最近正在作一个项目，前面建模等方面对我来说已经不是问题了（仅限2D),但是在求解的过程中，也经常遇到不收敛的问题，比较麻烦，我根据我的计算经验总结了几点，仅供参考！
1、有限元计算时，会经常遇到计算不收敛，出现这种情况的时候，按照下述两种方法，一般可以解决：
（1）网格剖分的问题
      可以采用加大网格的密度，我在计算中遇到网格过密，出现扭曲的单元过多而计算无法收敛，通过加大网格的密度，便可以计算收敛！
或者减小网格的密度，对局部单元过度扭曲的，采用局部优化的方法解决扭曲问题。
（2）计算参数的问题
      有时存在这样的情况，当某个计算参数过小的时候，会在计算中产生过大的位移，从而造成计算不收敛，此时可以采用修改计算参数的方法，使计算达到收敛。
2、改善计算收敛条件
   通过放宽程序自定义的收敛条件，可以解决计算不收敛的问题。我有时通过第2点即改善计算收敛条件，计算是收敛了，但是查看计算结果，总觉得不敢相信，我是从事稳定性分析的，每次作项目，都会遇到这问题，挺烦人的！这可能就是有限元的本质所在吧，建立一个漂亮的模型，也只是一个小小的开始，不等于万事大吉，对计算结果的分析和应用才是见真功夫的地方！
因为有感想，写了多一点，大家不见怪吧！

[ 本帖最后由 diaohuwang 于 2006-11-14 12:15 编辑 ]

dahaidahai 发表于 2006-11-15 17:00:12

做了下你的题目．．．．．

我把你的弹性模量改了－1e6～1e6．．．．．．．．．．．．．
页: [1]
查看完整版本: 是加载还是计算（收敛）方面出现问题了？（命令流）