ANSYS热应力分析实例.docx
《ANSYS热应力分析实例.docx》由会员分享，可在线阅读，更多相关《ANSYS热应力分析实例.docx（13页珍藏版）》请在冰点文库上搜索。
ANSYS热应力分析实例
ANSYS热应力分析实例
当一个结构加热或冷却时，会发生膨胀或收缩。
如果结构各部分之间膨胀收缩程度不同，和结构的膨胀、收缩受到限制，就会产生热应力。
热应力分析的分类
ANSY號供三种进行热应力分析的方法：
在结构应力分析中直接定义节点的温度。
如果所以节点的温度已知，贝冋以
通过命令直接定义节点温度。
节点温度在应力分析中作为体载荷，而不是节点自由度
间接法。
首先进行热分析，然后将求得的节点温度作为体载荷施加在结构应力分析中。
直接法。
使用具有温度和位移自由度的耦合单元，同时得到热分析和结构应
力分析的结果。
如果节点温度已知，适合第一种方法。
但节点温度一般是不知道的。
对于大
多数问题，推荐使用第二种方法一间接法。
因为这种方法可以使用所有热分析的功能和结构分析的功能。
如果热分析是瞬态的，只需要找出温度梯度最大的时间点，并将此时间点的节点温度作为荷载施加到结构应力分析中去。
如果热和结构
的耦合是双向的，即热分析影响结构应力分析，同时结构变形又会影响热分析（如大变形、接触等），则可以使用第三种直接法一使用耦合单元。
此外只有第三种方法可以考虑其他分析领域（电磁、流体等）对热和结构的影响。
间接法进行热应力分析的步骤
首先进行热分析。
可以使用热分析的所有功能，包括传导、对流、辐射和表面效应单元等，进行稳态或瞬态热分析。
但要注意划分单元时要充分考虑结构分
析的要求。
例如，在有可能有应力集中的地方的网格要密一些。
如果进行瞬态分
析，在后处理中要找出热梯度最大的时间点或载荷步。
热单元
结构单元
LINK32
LINK1
LINK33
LINK8
PLANE35
PLANE2
PLANE55
PLANE42
SHELL57
SHELL63
PLANE67
PLANE42
LINK68
LINK8
SOLID79
SOLID45
MASS71
MASS21
PLANE75
PLANE25
PLANE77
PLANE82
PLANE78
PLANE83
PLANE87
PLANE92
PLANE90
PLANE95
SHELL157
SHELL63
重新进入前处理，将热单元转换为相应的结构单元，表7-1是热单元与结构
单元的对应表。
可以使用菜单进行转换:
MainMenu>Preprocessor>ElementType>SwitchElementTyp，选择Thermalt
Structual。
但要注意设定相应的单元选项。
例如热单元的轴对称不能自动转换到结构单元中，需要手工设置一下。
在命令流中，可将原热单元的编号重新定义为结构单元，并设置相应的单元选项。
设置结构分析中的材料属性（包括热膨胀系数）以及前处理细节，如节点耦合、约束方程等。
读入热分析中的节点温度，
GUI:
Solution>LoadApply>Temperature>FromThermalAnalysis输入或选择热分析的结果文件名*.rth。
如果热分析是瞬态的，则还需要输入热梯度最大时的时间点或载荷步。
节点温度是作为体载荷施加的，可通过Utility
Menu>List>Load>BodyLoad>Onallnode列表输出。
设置参考温度，MainMenu>Solution>LoadSetting>RefereneeTemp
进行求解、后处理。
间接法热应力分析实例
问题描述
iwri™niMritd
fiM
3
中
国软件网5OFF6.匚口刚
热流体在代有冷却栅的管道里流动，如图为其轴对称截面图。
管道及冷却栅
的材料均为不锈钢，导热系数为hr-in-oF,弹性模量为28E6lb/in2泊松比为。
管
内压力为1000Ib/in2,管内流体温度为450oF,对流系数为1Btu/hr-in2-oF,外界流体温度为70oF，对流系数为Btu/hr-in2-oF。
求温度及应力分布。
菜单操作过程
设置分析标题
1、选择“UtilityMenu>File>ChangeTitle输入Indirectthermal-stressAnalysisofacoolingfin。
2、选择“UtilityMenu>File>ChangeFilename输入”IPE_FIN
进入热分析，定义热单元和热材料属性
1、选择“MainMenu>Preprocessor>ElementType>Add/Edit/Delete选择”PLANE55设定单元选项为轴对称。
2、设定导热系数：
选择“MainMenu>Preprocessor>MaterialPorps>Material
Models”点击Thermal，Conductivity，Isotropic，输入。
创建模型
1、创建八个关键点，选择“MainMenu>Preprocessor>Creat>Keypoints>On
ActiveCS，”关键点的坐标如下：
编号
1
2
3
4
5
6
7
8
X
5
6
12
12
6
6
5
5
Y
0
0
0
1
1
2、组成三个面：
选择“MainMenu>Preprocessor>Creat>Area>Arbitrary>Throuth
Kps”，由1,2,5,8组成面1;由2,3,4,5组成面2;由8,5,6,7组成面3。
3、设定单元尺寸，并划分网格：
“MainMenu>Preprocessor>Meshtool,设定
globalsize为，选择AREAMapped.Mesh，点击Pickall。
施加荷载
1、选择“UtilityMenu>Select>Entities>Nodes>Bylocation>XcoordinatesFrom
Full，”输入5,点击OK选择管内壁节点；
2、在管内壁节点上施加对流边界条件：
选择“Main
Menu>Solution>Apply>Convection>Onnodes，点击Pick,all，输入对流换热系数1,流体环境温度450。
3、选择“UtilityMenu>Select>Entities>Nodes>Bylocation>XcoordinatesFromFull，”
输入6,12，点击Apply；
4、选择“UtilityMenu>Select>Entities>Nodes>Bylocation>Ycoordinat，sReselect，”输入，1，点击Apply；
5、选择“UtilityMenu>Select>Entities>Nodes>Bylocation>Ycoordinat，sAlsoselect，”输入12，点击OK
6、在管外边界上施加对流边界条件：
选择“Main
Menu>Solution>Apply>Coivection>Onnodes，'点击Pick,all，输入对流换热系数，流体环境温度70。
求解
1、选择“UtilityMenu>Select>SelectEverything”
2、选择“MainMenu>Solution>SolveCurrentLS”
后处理
1显示温度分布：
选择“MainMenu>GeneralPostproc>PlotResult>Nodal
Solution>
Temperature。
重新进入前处理，改变单元，定义结构材料
1选择“MainMenu>Preprocessor>ElementType>SwitchElemTyp选择”
ThermaltoStructure0
2、选择“MainMenu>Preprocessor>ElementType>Add/Edit/Delete点击”
Option，将结构单元设置为轴对称。
3、选择“MainMenu>Preprocessor>MaterialPorps>MaterialModels输入材”
料的EX为28E6,PRXY为，ALPX为。
定义对称边界条件
1选择“UtityMenu>Select>Entities>Nodes>Bylocation>YcoordinatesFrom
Full，”
输入0，点击Apply;
2、选择“UtilityMenu>Select>Entities>Nodes>Bylocation>丫coordinat，sAlso
select，”输入1，点击Apply；
3、选择“MainMenu>Solution>Apply>Displacement>Symmetry.OnNodes点
击PickAll,选择Yaxis点击OK；
施加管内壁压力
1、选择“UtilityMenu>Select>Entities>Nodes>Bylocation>XcoordinatesFromFull，”输入5,点击OK
2、选择“MainMenu>Solution>Apply>Pressure>Onnode,点击”PickAll,输入1000。
设置参考温度
1、选择“UtilityMenu>Select>SelectEverything”
2、选择“MainMenu>Solution>_oads-Seting>RefereneeTemp输入70。
读入热分析结果
1、选择“MainMenu>Solution>Apply>Temperature>FromThermalAnalysis>选择。
求解
选择“MainMenu>Solution>SolveCurrentLS”
后处理
选择“MainMenu>GeneralPostpro>PlotResult>NodalSolution>Stress>VonMises”。
显示等效应力。
等效的命令流方法
/filename,pipe_fin
/TITLE,Thermal-StressAnalysisofacoolingfin
/prep7!
进入前处理
et,1,plane55!
定义热单元
keyopt,1,3,1!
定义轴对称
mp,kxx,1,!
定义导热系数
k,1,5!
建模
k,2,6
k,3,12
k,4,12,
k,5,6,
k,6,6,1
k,7,5,1
k,8,5,
a,1,2,5,8
a,2,3,4,5
a,8,5,6,7
esize，定义网格尺寸
amesh,all!
划分网格
eplot
finish
/solu!
热分析求解
nsel,s,loc,x,5选择内表面节点
sf,all,conv,1,450施加对流边界条件
nsel,r,loc,y,,1
nsel,a,loc,x,12
sf,all,conv,,70施加对流边界条件
nsel,all
/pse,conv,hcoef,1
nplot
solve!
求解生成文件
finish
/postl
pinsol,temp!
得到温度场分布
finish
/prep7!
重新进入前处理
etchg,tts!
将热单元转换为结构单元
plane42
keyopt,1,3,1!
定义轴对称特性mp,ex,1,28e6定义弹性模量
mp,nuxy,1,!
定义泊松比
mp,alpx,1,!
定义热膨胀系数
finish
/solu!
进入结构分析求解
nsel,s,loc,y,0选择对称边界
nsel,a,loc,y,1
dsym,symm,y!
定义对称条件
nsel,s,loc,x,5选择内表面
sf,all,pres,1000施加压力边界条件
nsel,all
/pbc,all,1
/psf,pres,,1
nplot
tref,70!
设定参考温度
ldread,temp,,,,,,rth!
读入节点温度
/pbc,all,0
/psf,pres,,0分布
/pbf,temp,,1
eplot
solve!
求解
finish
/post1,plnsol,s,eqv!
得到等效应力
finish
直接法热应力分析实例
问题描述
两个同心圆管之间有一个小间隙，内管中突然流入一种热流体，求经过
钟后外管表面的温度。
已知条件：
热膨胀系数：
5E-41/oF
泊松比：
导热系数：
10W/
密度：
7880Kg/m3
比热：
500J/
外管外半径：
m
外管内半径：
m
内管外半径：
内管内半径：
流体温度：
300oC
流体与内管内壁对流系数：
300W/
内、外管接触热导：
oC
命令流方法
/filename,contact_thermal
/title,contact_thermalexample
/prep7
ux,uy,temp
et,1,13,4,,1!
选择直接耦合单元PLANE13单元自由度为
!
定义为轴对称
et,2,48!
定义结构接触单元
keyopt,2,1,1!
设定接触单元的相应选项
keyopt,2,7,1
r,2,2e11,0,,,,!
定义接触单元实常数
mp,ex,1,2e11!
定义管材结构及热属性
mp,alpx,1,5e-5
mp,kxx,1,10
mp,dens,1,7880
mp,c,1,500
rect,,,0,!
建模
rect,,,0,
amesh,all
nsel,s,loc,x,!
将内管内壁的X方向位移及温度耦合
cp,1,ux,all
cp,2,temp,all
nsel,s,loc,x,!
将内管外壁的X方向位移及温度耦合
cp,3,ux,all
cp,4,temp,all
nsel,,x,!
将外管内壁的X方向位移及温度耦合
cp,5,ux,all
cp,6,temp,all
nsel,s,loc,x,!
将外管外壁的X方向位移及温度耦合
cp,7,ux,all
cp,8,temp,all
nsel,s,loc,y,!
将内管顶部节点的丫方向位移及温度耦合
nsel,r,loc,x,0,
cp,9,uy,all
nsel,s,loc,y,!
将外管顶部节点的丫方向位移及温度耦合
nsel,r,loc,x,,
cp,10,uy,all
nsel,s,loc,x,!
创建接触单元
cm,cont,node
nsel,s,loc,x,
cm,targ,node
type,2
real,2
gcgen,cont,targ,3
/solu
antype,trans!
瞬态分析
tunif,20!
初始平均温度
tref,20!
参考温度
sfl,4,conv,300,,300!
内管内壁对流边界
sfl,6,conv,10,,20!
外管外壁对流边界
nsel,s,loc,y,0!
约束所有底边单元的丫向位移
d,all,uy,0
time,180!
载荷步时间
deltime,10,5,15!
定义时间步长
outres,all,all
kbc,1
autots,on!
自动时间步长
allsel
solve!
求解
/postl
pinsol,temp!
显示温度分布
pinsol,s,eqv!
显示等效应力
升级会员 