热应力分析例题.docVIP

例1：有一截面为圆环形的输暖管道，如图7—17内外管道半径分别为200mm、800mm，管道内水的半径为80℃，管外表层温度为10℃，求管道内的热应力分布(假设管道内充满水)。 材料参数： 弹性模量Z—120GPa 泊松比v＝0.3 线膨胀系数α＝1.3×10^-6m/m℃ 导热系数k=l.2w／m℃ 分析：该问题属于轴对称问题。在进行有限元计算时，沿管道横截面取宽度50mm的矩形截面(如图7—18为计算模型，首先采用间接法进行对其进行热应力分析，然后再采用直接法进行分析，最后对二者求解结果进行比较分析。 间接法求解： Finish /clear /filname,thermal stresses in a long cylinder-indirect solution /prep7 Et,1,plane55,,,1 !定义轴对称单元 Mp,kxx,1,1.2 Rectng,0.2,0.8,0,0.05 Type,1 Lsel,s,line,,1,3,2 Lesize,all,,,10 Lsel,s,line,,2,4,2 Lesize,all,,,2 Amesh,1 Finish /solu !热传导求解 Antype,static Lsel,s,line,,4 Nsll,s,1 D,all,temp,80 Lsel,s,line,,2 Nsll,s,1 D,all,temp,10 Allsel Outpr,basic,all Solve finish /post1 Plnsol,temp Finish /prep7 !重新进入前处理，进行热应力耦合分析 Etchg,tts !转化单元类型热单元55为结构单元42 Keyopt,1,3,1 !定义单元关键字选项3为1（轴对称） Keyopt,1,6,1 !定义单元关键字选项6为1（无表面输出） Mp,ex,1,120e9 !输入材料结构性能参数 Mp,alpx,,1.3e-6 Mp,nuxy,1,0.3 Lsel,s,line,,3 Nsll,a,1 Cp,7,uy,all allsel Cp,8,ux,1,14,24 Cp,9,ux,2,12,13 Finish /solu !热应力求解 Antype,static D,1,uy,,,11 !施加位移约束 Ldread,temp,,,,,,rth Solve Finish /post1 Plnsol,s,x Plnsol,s,y Plnsol,s,z Finish 直接法求解： Finish /clear /filname,thermal stresses in a long cylinder-direct solution /prep7 Et,1,plane13,,,1 !定义热应力耦合单元 Keyopt,1,1,4 !定义单元关键字选项1=4（ux,uy,temp,az自由度） Keyopt,1,3,1 !定义单元关键字选项3=1（轴对称） Mp,ex,1,120e9 ！输入材料性能参数 Mp,nuxy,1,0.3 Mp,alpx,,1.3e-6 Mp,kxx,1,1.2 Rectng,0.2,0.8,0,0.05 !创建几何模型 Type,1 !进行网格划分 Lsel,s,line,,1,3,2 Lesize,all,,,10 Lsel,s,line,,2,4,2 Lesize,all,,,2 Amesh,1 Lsel,s,line,,3 Nsll,s,1 Cp,7,uy,all !耦合节点 Allsel Cp,8,ux,1,14,24 Cp,9,ux,2,12,13 Finish /solu Antype,static !热应力耦合求解 Lsel,s,line,,4 Nsll,s,1 D,all,temp,80 !施加温度载荷 Lsel,s,line,,2 Nsll,s,1 D,all,temp,10 !施加温度载荷 Allsel D,1,uy,,,11 !施加位移约束 Outpr,basic,all solve Finish /post1 Plnsol,temp !输出温度场 Plnsol,s,x !输出应力场 Plnsol,s,y Plnsol,s,z finish