1、 当弹性体的温度变化时，其体积将趋于膨胀和收缩，若外部的约束或内部的变形协调要求而使膨胀或收缩不能自由发生时，结构中就会出现附加的应力。这种因温度变化而引起的应力称为热应力，或温度应力。 忽略变温对材料性能的影响，为了求得温度应力，需要进行两方面的计算：（1）由问题的初始条件、边界条件，按热传导方程求解弹性体的温度场，而前后两个温度场之差就是弹性体的变温。（2）按热弹性力学的基本方程求解弹性体的温度应力。本章将对这两方面的计算进行简单的介绍。第六章第六章 温度应力问题的基本解法温度应力问题的基本解法6-4 6-4 按位移求解温度应力的平面问题按位移求解温度应力的平面问题6-3 6-3 温度场的
2、边界条件温度场的边界条件6-2 6-2 热传导微分方程热传导微分方程6-1 6-1 温度场和热传导的基本概念温度场和热传导的基本概念6-5 6-5 位移势函数的引用位移势函数的引用6-6 6-6 轴对称温度场平面热应力问题轴对称温度场平面热应力问题6-1 6-1 温度场和热传导的基本概念温度场和热传导的基本概念1.温度场：在任一瞬时，弹性体内所有各点的温度值的总体。用T表示。 不稳定温度场或非定常温度场：温度场的温度随时间而变化。 即 T=T（x，y，z，t） 稳定温度场或定常温度场：温度场的温度只是位置坐标的函数。 即 T=T（x，y，z） 平面温度场：温度场的温度只随平面内的两个位置坐标而
3、变。 即 T=T（x，y，t）2.等温面：在任一瞬时，连接温度场内温度相同各点的曲面。显然，沿着等温面，温度不变；沿着等温面的法线方向，温度的变化率最大。T+2TT+TTT-Txoy），（），（），（znCOSnTzTynCOSnTyTxnCOSnTxTnTn0nT3.温度梯度：沿等温面的法线方向，指向温度增大方向的矢量。用T T表示，其大小用 表示。其中n n为等温面的法线方向。温度梯度在各坐标轴的分量为0n取 为等温面法线方向且指向增温方向的单位矢量，则有T T（1）4.热流速度：在单位时间内通过等温面面积S 的热量。用 表示。 热流密度：通过等温面单位面积的热流速度。用 表示，则有 dt
4、dQqSdtdQq/其大小为SdtdQnq/0(2)5.热传导基本定理：热流密度与温度梯度成正比而方向相反。即 q(3)nTq由（1）和（3）可见，热流密度的大小可见，导热系数表示“在单位温度梯度下通过等温面单位面积的热流速度”。SnTdtdQ/称为导热系数。由（1）、（2）、（3）式得T T热流密度在坐标轴上的投影可见：热流密度在任一方向的分量，等于导热系数乘以温度在该方向的递减率。zTqyTqxTqzyx 热量平衡原理：在任意一段时间内，物体的任一微小部分所积蓄的热量，等于传入该微小部分的热量加上内部热源所供给的热量。6-2 6-2 热传导微分方程热传导微分方程dxxqqxxxqxyz 取
5、图示微小六面体dxdydz。假定该六面体的温度在dt时间内由T 升高到 。由温度所积蓄的热量是 ，其中 是物体的密度，C 是单位质量的物体升高一度时所需的热量比热容。dttTTdttTdxdydzC 在同一段时间dt内，由六面体左面传入热量qxdydzdt，由右面传出热量 。因此，传入的净热量为dydzdtdxxqqxx)(dxdydzdtxqxxTqx将 代入可见：dxdydzdtxT22dydzdxdtyT22由左右两面传入的净热量为：由上下两面传入的净热量为：dzdxdydtzT22由前后两面传入的净热量为：因此，传入六面体的总净热量为：简记为：dxdydzdtzTyTxT)(22222
6、2Tdxdydzdt2 假定物体内部有正热源供热，在单位时间、单位体积供热为W，则该热源在时间dt内所供热量为Wdxdydzdt。 根据热量平衡原理得：WdxdydzdtTdxdydzdtdxdydzdtxTc2cWTctT2化简后得：ca 记cWTatT2则这就是热传导微分方程。6-3 6-3 温度场的边值条件温度场的边值条件 初始条件： 边界条件分四种形式： 第一类边界条件 已知物体表面上任意一点在所有瞬时的温度，即 其中Ts 是物体表面温度。 第二类边界条件 已知物体表面上任意一点的法向热流密度，即 其中角码 s 表示“表面”，角码n 表示法向。),()(0zyxfTt)(tfTs)()
7、(tfqsn 为了能够求解热传导微分方程，从而求得温度场，必须已知物体在初瞬时的温度，即所谓初始条件；同时还必须已知初瞬时以后物体表面与周围介质之间热交换的规律，即所谓边界条件。初始条件和边界条件合称为初值条件。 第三类边界条件 已知物体边界上任意一点在所有瞬时的运流（对流）放热情况。按照热量的运流定理，在单位时间内从物体表面传向周围介质的热流密度，是和两者的温差成正比的，即 )()(essnTTqesTT其中Te是周围介质的温度； 称为运流放热系数，或简称热系数。 第四类边界条件 已知两物体完全接触，并以热传导方式进行热交换。即6-4 6-4 按位移求解温度应力的平面问题按位移求解温度应力的
8、平面问题 设弹性体内各点的温变为T。对于各向同性体，若不受约束，则弹性体内各点的微小长度，都将产生正应变 （ 是弹性体的膨胀系数），这样，弹性体内各点的形变分量为T0,xyzxyzzyxTTETETEyxzzxzyyzyxx)(1)(1)(1 但是，由于弹性体所受的外在约束以及体内各部分之间的相互约束，上述形变并不能自由发生，于是就产生了应力，即所谓温度应力。这个温度应力又将由于物体的弹性而引起附加的形变，如虎克定理所示。因此，弹性体总的形变分量是：对于平面应力的变温问题，上式简化为xyxyzxzxyzyzEEE)1(2)1(2)1(2xyxyxyyyxxETETE)1 (211这就是平面应力
9、问题热弹性力学的物理方程。将应力分量用形变分量和变温T表示的物理方程为：xyxyxyyyxxETEETEE)1 ( 21)(11)(122几何方程仍然为：yuxvyvxuxyyx,将几何方程代入物理方程，得用位移分量和变温T 表示的应力分量）（）（）（）（yuxvETExuyvETEyvxuExyyx12111122将上式代入不计体力的平衡微分方程00 xyyxxvyyxx简化得：这就是按位移求解温度应力平面应力问题的微分方程。 同理，将应力分量代入无面力的应力边界条件00sxysysyxsxlmml）（）（）（）（0121210121212222222222yTyxuxvyvxTyxvyux
10、u）（）（1）简化后得：这是按位移求解温度应力平面应力问题的应力边界条件。 位移边界条件仍然为：vvuuss ， 将式(1)、(2)与第二章2-8中式(1)、(2)对比，可见TmyuxvlxuyvmTlxvyumyvxulssss）（）（）（）（）（）（121121（2）代替了体力分量 X 及 Y ，而：则得到在平面应变条件下的相应方程。yTExTE11及11TEmTEl及代替了面力分量 及 。XY 对于温度应力的平面应变问题，只须将温度应力的平面应力问题的）（换成换成换成1112EE6-5 6-5 位移势函数的引用位移势函数的引用 由上一节知：在平面应力的情况下按位移求解温度应力问题时，须使
11、位移分量u 和v 满足微分方程：0121210121212222222222yTyxuxvyvxTyxvyuxu）（）（并在边界上满足位移边界条件和应力边界条件。实际求解时，宜分两步进行：（1）求出上述微分的任意一组特解，它只需满足微分方程，而不一定要满足边界条件。（2）不计变温T，求出微分方程的一组补充解，使它和特解叠加以后，能满足边界条件。 引用一个函数 ，将位移特解取为：yvxu,），（yx函数 称为位移势函数。以 和 分别作为u和v代入微分方程，简化后得：uvyTyxTx）（）（1122由于 和 都是常量，所以取：T）（ 12时， 满足微分方程。因此 ， 可以作为微分方程的一组特解。将
12、），（yx u v以及211Tyvxu,代入位移分量和变温T表示的应力分量表达式)()1 (21)(11)(122yuxvETExuyvETEyvxuExyyx可得相应位移特解的应力分量是：yxExEyExyyx22222111 设 ， 为位移的补充解，则 ， 需满足齐次微分方程：uu v v02121021212222222222yxuxvyvyxvyuxu相应于位移补充解的应力分量为（注意不计变温，即T=0）：)()1(2)(1)(122yuxvExuyvEyvxuExyyx总的应力分量是：, , zzzyyyxxx需满足应力边界条件。在应力边界问题中（没有位移边界条件），可以把相应于位移
13、补充解的应力分量直接用应力函数来表示，即其中的应力函数 可以按照应力边界条件的要求来选取。yxxyxyyx22222, 在平面应变条件下，将上述各方程中的）（换成换成换成1112EE这样总的位移分量是：, vvvuuu需满足位移边界条件。例1 图示矩形薄板中发生如下的变温：其中的T0 是常量。若 ，试求其温度应力。）（2201byTTbaxyoaabb解：位移势函数 所应满足的微分方程为）（）（220211byT22ByAy ）（）（220111222byTByA21212100bTBTA）（，）（比较两边系数，得代入上式，得取将A，B回代，得位移势函数于是相应于位移特解的应力分量为 为求补充
14、解，取 可得所需要的相应于位移补充解的应力分量：）（）（24201221byyT2cy0, 0),1 (220 xyyxbyTE0, 0,2222yxcyxyyx00)1 (2220 xyxyxyyyyxxxbyTEc因此，总的应力分量为0)( , 0)(0)( , 0)(byxybyyaxxyaxx边界条件要求显然，后三个条件是满足的；而第一个条件不能满足，但由于 ，可应用圣维南原理，把第一个条件变换为静力等效条件，即，在 的边界上， 的主矢量及主矩等于零：将baaxx0)(, 0)(ydydyaxbbxaxbbx)1 (2220byTEcx代入上式，求得于是矩形板的温度应力为：0322TE
15、c00)31(220 xyyxbyTE6-6 6-6 轴对称温度场平面热应力问题轴对称温度场平面热应力问题 对于圆形、圆环及圆筒等这类轴对称结构弹性体，若其变温也是轴对称的T=T（r），则可简化为轴对称温度场平面热应力问题。轴对称温度场平面热应力问题，宜采用极坐标求解。 不考虑体积力平面应力问题平衡方程02101rrrrrrrrrrr 在轴对称问题中得到简化，其第二式自然满足；而第一式成为0rrrr 几何方程简化为rudrdurrr, 物理方程简化为TETErrr)(1)(1 将应力用应变表示1)(11)(122TEETEErrr 将几何方程代入上式，然后将其代入平衡方程，得按位移求解轴对称热
16、应力的基本方程：或写成： drdTrudrdurdrudrrr)1 (1222drdTrudrdrdrdr)1 ()(1积分两次可得到轴对称问题位移分量：rBArTrdrrurar)1 (式中A，B为任意常数，积分下限取为a。由上式可得应力分量：0)1 ()1(1)1 ()1(1222222rrararTErBAETrdrrErBAETrdrrE其中常数A，B由边界条件确定。 在平面应变的情况下，只需在以上各式中将）（换成换成换成1112EE例2 设有一厚壁圆筒，内半径为a，外半径为b。从一均匀温度加热，内表面增温Ta ，外表面增温Tb，如图所示。试求筒内无热源，热流稳定后的热应力。abTaT
17、b 得无热源，热流稳定后的热传导微分方程为解：首先求温度场。由热传导微分方程cWzTyTxTtT)(22222202T对于轴对称温度场有 积分两次得：0)(10122drdTrdrdrTdrdrdrd）（或 由边界条件：求出A，B后回代，得温度场：BrAT lnbbraarTTTT）（）（baraTabrbTTbalnlnlnln积分后得TEababrbTTEabrbabrbTTEabrbabrbTTEbazbabar12ln1ln21211ln1ln1211lnln122222222222）（）（）（）（）（）（21112222222222222TTrdrabETrTrdrTrdrabarr
18、ETrdrTrdrabarrEbazrabarabar）（）（）（将T代入平面应变问题应力表达式练习6.1 图示矩形薄板中发生变温试求温度应力（假定a远大于b）xyoaabbbyTT2cos0解：byTT2cos1102取可解得所以byA2cos20214TbAbyTb2cos14022由此得01012cos1222022yxExEbyETyExyyx取2cy则0022 22 22 yxxcyxyyx所以0022cos 0 xyxyxyyyyxxxcbyET边界条件0, 0, 0axxybyxybyy显然满足由0dybbaxx即0422cos00cETbdycbyETbb得ETc0而边界条件0
19、ydyaxbbx恒成立。故002cos20 xyyxbyET练习6.2 已知半径为b的均质圆盘，置于等温刚性套箍内，圆盘和套箍由相同的材料制成，设圆盘按如下规律加热22011brTTT套箍温度则保持为常温T0，而由此温度所引起的应变可以忽略，试求距圆盘中心为r处的压应力值。解： 轴对称平板的平面应力问题的位移的表达式为 arcrcTrdrrur2011由于在中心处，即r=0处，u为有限值，因此c2=0。当r=b时，u=0，代入（a）即 bbcbbbbTT0421124201由（b）式得01141TTc将c1的值代入22122221221111111111rccEETTrdrrErccETrdrrErarar得220122013134134brTTEbrTTEr
