热应力说明书.docVIP

 PAGE \* MERGEFORMAT 11 热工及热应力基础 三级项目报告 组 员 指导教师 一、提出问题 某单程换热器的其中一根换热管和与其相连的两端管板结构，壳程介质为空气，管程介质为液体操作介质，换热管材料为不锈钢，膨胀系数为16.56×10-6°C-1,泊松比为0.3，弹性摸量为1.72×105MPa，热导率为15.1W/(m?°C)管板材料也为不锈钢，膨胀系数为17.79×10-6°C-1，泊松比为0.3，弹性摸量为1.73×105MPa，热导率为15.1W/(m?°C)。壳程气体温度为20°C，表面传热系数为3000W/m2°C，壳程压力为8.1MPa，管程液体温度为250°C，表面传热系数为426W/m2°C，管程压力为5.7MPa。 换热管内径0.01259m，管板厚度0.05m，换热管长度0.5m,部分管板长和宽度均为0.013m。 二、建立数学模型 此问题是一维稳态导热问题，温度场为常热性、无热源、一维稳态温度场。可简化为无限长圆筒壁的稳态导热问题。 三、温度场理论求解。 应采用圆柱坐标系下的拉普拉斯方程进行求解，由于 ?t?z=?t?θ=0 则导热方程式变为： d2tdr2+1rdtdr=0 即导热方程为： ddrrdtdr=0 （a） 边界条件为第一类边界条件： r=r1时,t=t1r=r2时,t=t2 （b） 式（a）、式（b）构成了求解此问题的数学模型，对式（a）积分两次，得： t=c1lnr+c2 （c） 将边界条件式（b）代入式（c），可求出积分常数c1和c2： c1=t2-t1lnr2r1 c2=t1-t2-t1lnlnr1lnr2r1 则温度场为： t=t1-t1-t2lnrr1lnr2r1 可见圆筒内的温度分布为对数曲线 r1=6.48mm t1=200° r2=9.53mm t2=25° r3=7.00mm t3=164.65° r4=8.00mm t4=104.11° r5=9.00mm t5=50.71° 四、温度场数值求解。 在工程中常遇到圆柱形或其他轴对称形状的物体，如机器的轴、轮、汽轮机和燃气轮机的转子，气缸、内燃机气阀、气缸套等，求解这类问题的零件温度场时，用圆柱坐标就比较方便。 图表示一个几何形状系沿轴对称，如果温度场也是对称的话，则t=tx,y，无内热源的稳定温度场微分方程为 ?2t?r2+1r?t?r+?2t?x2=0 （a） 这时划分网格的线族用 ri=i?ri=0,1,2,3??xi=j?xj=0,1,2,3?? （b） 式中?r和?x分别为r方向和x方向的步长。因为温度场是轴对称的，所以在r=0处，?t?r=0。于是，在讨论其温度分布时。可以取旋转体截面的一半来分析，如图所示。对轴对称，可视作绝热边界。然后对求解区域的任一内部节点ri,xi，列出相应的差分方程式。 αtαrri,xi≈ti+1,j-ti-1,j2?r?2t?r2ri,xi≈ti+1,j-2ti,j+ti-1,j?r2?2t?x2ri,xi≈ti,j+1-2ti,j+ti,j-1?x2 （c） 将式（c）的相应项代到式（a）中，就得到圆柱坐标的二维稳定导热差分方程式 ti+1,j-2ti,j+ti-1,j?r2+1riti+1,j-ti-1,j2?r+ti,j+1-2ti,j+ti,j-1?x2=0 根据式（b）的关系：ri=i?r,代入上式，并经过整理后，得到节点（i,j）的计算式 ti,j=121+?r?x21+12iti+1,j+1-12iti-1,j+?r?x2ti,j+1+?r?x2ti,j-1 如果取?r=?x，则上式可简化为 ti,j=141+12iti+1,j+1-12iti-1,j+ti,j-1+ti,j+1 计算结果： 温度场分布 五、热应力理论求解。 该问题为无限长圆柱体的热应力问题。对于无限长圆柱体或长圆柱体，其内径为a，外径为b，温度分布T=T（r），与轴向坐标无关。假如圆柱体两端是自由的，可先作平面应变来处理。 Ee=E1-ν2 νe=ν1-ν αe=αt1+ν 无限长圆柱体或长圆柱体热应力和唯一的解 σθ=αtE1-ν1r2r2+