我最近在用直接耦合法算瞬态热应力,总是得不到一个满意的结果,很头痛。现
在发现了一个问题,就是我要先计算初始温度场,然后再加载初始温度场进行计
算,这是我算的初始温度场
然后我计算一分钟后的温度场就变成这样了
我在结算热应力之前用ldread读入了初始温度场的温度,在output窗口也看得到
它已经作为体载荷加载到各个结点上了
请各位大侠帮小弟看看这是什么个情况,小弟万分感激。
附命令流:
/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 |