论坛 › C04：ANSYS--流体软件（Fluen\CFX\Icepark) › 请教：关于计算热应力的问题

leeppon 发表于 2010-8-18 21:01:44

请教：关于计算热应力的问题

我最近在用直接耦合法算瞬态热应力，总是得不到一个满意的结果，很头痛。现
在发现了一个问题，就是我要先计算初始温度场，然后再加载初始温度场进行计
算，这是我算的初始温度场
http://att.zju88.org/files/133388.ced2d23e.jpg
然后我计算一分钟后的温度场就变成这样了
http://att.zju88.org/files/133389.d46fd83c.jpg
我在结算热应力之前用ldread读入了初始温度场的温度，在output窗口也看得到
它已经作为体载荷加载到各个结点上了
http://att.zju88.org/files/133390.13d856eb.jpg
http://att.zju88.org/files/133391.99ae4cf0.jpg
请各位大侠帮小弟看看这是什么个情况，小弟万分感激。
附命令流：
/clear,start
finish
/FILNAME,warm_250
/TITTLE,THERMAL_STRESS_ANALYSIS_OF_ROTOR_DIRECT__warm_250
/PREP7
/UNITS,MKS!使用温度为摄氏度，其他均为国际标准单位的单位系统
!以下省略建模信息

mat,1!定义材料属性
MPTEMP,,,,,,,,
MPTEMP,1,100
MPTEMP,2,200
MPTEMP,3,300
MPTEMP,4,400
MPTEMP,5,500
MPTEMP,6,600
!定义比热容

MPDATA,C,1,,479.53
MPDATA,C,1,,523.66
MPDATA,C,1,,563.98
MPDATA,C,1,,607.48
MPDATA,C,1,,668.94
MPDATA,C,1,,760.9
!定义导热系数

MPDATA,KXX,1,,38.9
MPDATA,KXX,1,,38.1
MPDATA,KXX,1,,33.9
MPDATA,KXX,1,,33.1
MPDATA,KXX,1,,30.1
MPDATA,KXX,1,,26.4

!定义密度
MPDE,DENS,1
MPDATA,DENS,1,,7750
MPDATA,DENS,1,,7750
MPDATA,DENS,1,,7750
MPDATA,DENS,1,,7750
MPDATA,DENS,1,,7750
MPDATA,DENS,1,,7750
ET,1,PLANE55,2,,1

!定义热膨胀系数
MPTEMP,,,,,,,,
MPTEMP,1,100
MPTEMP,2,200
MPTEMP,3,300
MPTEMP,4,400
MPTEMP,5,500
MPTEMP,6,600
UIMP,1,REFT,,,22
MPDE,ALPX,1
MPDATA,ALPX,1,,11.5e-6
MPDATA,ALPX,1,,12e-6
MPDATA,ALPX,1,,12.5e-6
MPDATA,ALPX,1,,13.4e-6
MPDATA,ALPX,1,,13.8e-6
MPDATA,ALPX,1,,14e-6

!定义弹性模量和泊松比
MPTEMP,,,,,,,,
MPTEMP,1,100
MPTEMP,2,200
MPTEMP,3,300
MPTEMP,4,400
MPTEMP,5,500
MPTEMP,6,600
MPDE,EX,1
MPDE,NUXY,1
MPDATA,EX,1,,2.14e11
MPDATA,EX,1,,2.06e11
MPDATA,EX,1,,1.96e11
MPDATA,EX,1,,1.87e11
MPDATA,EX,1,,1.78e11
MPDATA,EX,1,,1.74e11
MPDATA,NUXY,1,,0.33
MPDATA,NUXY,1,,0.33
MPDATA,NUXY,1,,0.33
MPDATA,NUXY,1,,0.33
MPDATA,NUXY,1,,0.33
MPDATA,NUXY,1,,0.33

et,1,1plane13!选择13单元
keyopt,1,1,4
keyopt,1,3,1

mshkey,1!画网格
smrtsize,2
aesize,1,0.02
amesh,all

/solu
!计算初始温度场
ANTYPE,TRANS,NEW
TRNOPT,FULL

DL,8,,TEMP,70
DL,149,,TEMP,70
dl,160,,temp,250
tref,20
time,1e-6
timint,off

!以下省略定义表边界条件，不是我想省略，主要是太长了...都是类似这样的命
令
!汽封
SFL,110,CONV,%stage14_seal_conv%,,%stage14_seal_temp%
!级前
SFL,219,CONV,%stage14_front_conv%,,%stage14_front_temp%
SFL,111,CONV,%stage14_front_conv%,,%stage14_front_temp%

ldread,temp,,,1e-6,,warm_250_direct_intital_temp,rst!加载初始温度场
timint,on!打开积分

!施加角速度和角加速度
OMEGA,,%omega%
DOMEGA,,%angle_acceleration%

autots,on!打开自动时间步长
!TIME,timeunit
time,i60
DELTIME,1,10,60,on
!NSUBST,1
KBC,1
outpr,basic,all
outres,basic,all
allsel,all
solve

minxiang 发表于 2010-8-19 10:37:45

你的问题是什么？

yzzyinyin 发表于 2010-8-19 11:16:11

稳态分析完成后，进行瞬态分析时，应该把以前施加的初始温度载荷删除。

leeppon 发表于 2010-8-22 13:53:54

2# minxiang 现在计算的热应力结果的数量级和预想的结果大一个数量级，只计算转速的话，结果比较合理，但是只计算热应力以及热应力和离心力一块计算的话，结果就很诡异，找了好久的问题，网格划分没有警告，材料属性检查好几遍也没错误，所以就不知道问题出在哪里了...

leeppon 发表于 2010-8-22 13:57:18

3# yzzyinyin 不好意思，前几天有点忙忘记上论坛了
谢谢你的回复，我在初始温度场有三个：
DL,8,,TEMP,70
DL,149,,TEMP,70
dl,160,,temp,250

其中前两个地方在整个过程中都是70度，所以就没删除
最后一个在加载之前已经删除掉了

leeppon 发表于 2010-8-24 09:08:15

请大侠指点啊，在err文件里也只发现这四个错误信息：
*** WARNING ***                         CP =       4.259   TIME= 16:02:53
Element 1 references undefined MURX or BH table for material 1.   
这个是我用13单元，打开timint的时候也把mag打开了，但是并没有用到关于磁的东西，所以这个应该可以忽略，如果我算之前把timint，off，mag，就会没有这个警告了
*** WARNING ***                         CP =       4.774   TIME= 16:02:57
Material property ALPX of material 1 of element 10857 is evaluated at a
temperature of 20, which is below the supplied temperature range.      
Temperature range checking terminates.
我定义的材料的属性是随时间变化的，20，100，200，300，400，500
*** WARNING ***                         CP =       8.190   TIME= 16:03:05
There are 1 small equation solver pivot terms.May be a transitional   
effect caused by the Newton-Raphson logic.   
这个不知道什么意思
*** WARNING ***   SUPPRESSED MESSAGE    CP =      35.100   TIME= 16:03:41
The calculated reference FORCE CONVERGENCE VALUE = 4.941050718E-04 is   
less than the default minimum of 1 (MINREF field of CNVTOL command).   
Check results carefully.   
这个应该是关于收敛的，也不知道怎么解决

请大侠帮帮忙，最近快崩溃了，检查了好几遍了，都不知道自己错在哪里...主要是应力值偏大而且应力分布也不理想

yzzyinyin 发表于 2010-8-25 08:50:46

3# yzzyinyin 不好意思，前几天有点忙忘记上论坛了
谢谢你的回复，我在初始温度场有三个：
DL,8,,TEMP,70
DL,149,,TEMP,70
dl,160,,temp,250

其中前两个地方在整个过程中都是70度，所以就没删除
最后一个在 ...
leeppon 发表于 2010-8-22 13:57 http://forum.simwe.com/images/common/back.gif

前两个也应该删除，如果是施加的均匀温度或者参考温度，才可以不删除。你试试。

aganspring 发表于 2010-8-29 21:55:01

我也是做这个的，也出现了同样的问题，我还没有先进行稳态热分析，看了你写的之后有所启发。楼主有什么进展没？交流下啊。我qq379483295 6# leeppon

leeppon 发表于 2010-9-2 15:19:34

7# yzzyinyin 我现在用间接耦合法，没有删除那两个温度载荷，现在得到的结果还比较满意，谢谢您的帮忙~~

leeppon 发表于 2010-9-2 15:21:26

8# aganspring 我现在用直接耦合法算依旧得不到满意的结果（当使用均匀的初始温度场的时候结果是比较满意的），现在改用了间接耦合法，算出了比较满意的结果

查看完整版本: 请教：关于计算热应力的问题