《第十章_温度应力》由会员分享，可在线阅读，更多相关《第十章_温度应力（35页珍藏版）》请在金锄头文库上搜索。
1、第十章_温度应力第十章 温度应力前面几章讨论了弹性体在外荷载或边界位移（基础沉降）作用下引起的应力及变形，本章将讨论弹性体由于温度的改变而引起的应力与变形。当弹性体的温度改变时，它的每一部分都将由于温度的升高或降低而趋于膨胀或收缩，但是，由于弹性体所受到的外在约束，以及弹性体自身各个部分之间的相互约束，使得弹性体的膨胀或收缩并不能自由地发生，于是就产生了应力，这种由于温度改变所产生的应力称之为温度应力。温度应力是混凝土结构产生裂缝的主要原因之一，温度裂缝对结构造成严重的危害，对结构的安全性与耐久性产生重要影响。为了确定弹性体内的温度应力，需进行两方面的计算：（1）按照热传导理论，根据弹性体的热
2、学性质、内部热源、初始条件和边界条件，计算弹性体内各点在各瞬时的温度，即确定温度场，前后两个时刻的温度场之差就是弹性体的温度改变。（2）按照“热弹性力学”，根据弹性体的温度改变来求出体内各点的温度应力，即确定应力场。关于温度场的计算，在传热学中已有详细讨论，本章仅介绍传热学的基本方程及初始、边界条件。本章将主要讨论温度应力的计算，暂不考虑荷载或边界位移对应力的影响，对于实际既有荷载作用、边界位移作用以及温度改变作用的弹性力学问题，可分别单独计算单个作用因素下的应力结果，然后应用叠加原理，将单个作用因素下的应力结果叠加起来，即可得到实际问题的解答。 10-1温度场和热传导方程（1）温度场一般情况
3、下，在热传导的过程中，弹性体内各点的温度随位置和时间而变化，因而温度T 是位置坐标和时间t 的函数：(,)T T x y z t = (10.1)在任一瞬时，所有各点的温度值的总体称为温度场。一个温度场，如果它的温度随时间而变化，如式(10.1)所示，就称为不稳定温度场；如果温度场不随时间而变化，即有：0T t?=?，就称为稳定温度场。在稳定温度场中，温度只是位置坐标的函数，即：(,)T T x y z = （a ）如果稳定温度场只随着两个位置坐标而变化，即平面稳定温度场为：(,)T T x y = （b ）温度场的等温线可以表示为：(,)T x y C = （c ）式中：C 为常数，平面温度
4、场的等温线如图10-1所示。（2）温度梯度沿着等温面（或等温线，见如图10-1），温度不变；但沿着其他方向，温度都有变化，沿着等温面（或等温线）的法线方向，温度的变化率最大，这个最大的变化率可由温度梯度来表示：T T T T i j k x y z?=+?r r r （d ）或：(T T T i jx y)?=+?r r平面温度场 （e ）温度梯度为向量场，方向总是指向温度增加的一方（见图10-1）。 图10-1平面温度场的等温线 （3）热传导定律在弹性体内部，热量从弹性体的一个部位传到另一个部位，热量的传导符合下列热传导定律：q T =?r(10.2)式中：q 为热流密度，表示单位时间内通过
5、单位等温面面积的热量，它的单位是热量 长度r-2 时间-1，为热传导系数，与弹性体的材料有关。(10.2)式表明热流密度与温度梯度成正比，方程中的负号表示热量流动的方向与温度梯度相反，即热量总是从高温向低温方向流动。（4）牛顿冷却定律弹性体与其周边介质的热交换符合下列牛顿冷却定律： (ns sq T =?)a T (10.3)式中：s T 为弹性体表面温度，为弹性体周边介质温度；a T n s q 为弹性体外法向表面流出的热流密度（这里假定s a T T ）；为热交换系数。将热传导定律(10.2)式用于弹性体表面，有：nssTq n ?=?，所以牛顿冷却定律又可表示为： ()s a sTT T
6、 n ?=? (10.4) （5）热传导微分方程热传导微分方程的建立，是基于如下的热学基本原理：在任意一段时间内，物体的任一微小部分由于温度升高所积蓄的热量，等于传入该微小部分的热量加上内部热源所供给的热量。取直角坐标系并取微小六面体，如图10-2所示。假定该六面体的温度在时间内由T 升高到dxdydz dt dt t T T ?+。由于温度升高了dt t T?，它所积蓄的热量是dt tT dxdydz c ?，其中是物体的密度；是比热系数，也就是单位质量的物体温度升高一度时所需要的热量。c图10-2在同一时间内，由六面体左面传入热量，由右面传出dt dydzdt q x dydzdt dx
7、xq q xx )(?+。因此，传入的净热量为dxdydzdt x q x ?，并由热传导定律（10.2）式，可见其等于dxdydzdt xT22?。同样，由上下两面及前后两面传入的净热量分别为dxdydzdt y T 22?及dxdydzdt z T22?。这样，传入六面体的净热量总共是dxdydzdt z Ty T xT (222222?+?+?，或简写为。 Tdxdydzdt 2?设该六面体的内部有热源，其强度为W （在单位时间、单位体积内供给的热量），则该热源在时间内所供给的热量为Wdxdydzdt 。在这里，供热的热源作为正的热源，例如金属通电时发热，混泥土硬化时发热，水分结冰时发热
8、，等等；吸热的热源作为负的热源，例如水分蒸发时吸热，冰粒溶解时吸热，等等。dt 于是，根据热量平衡原则有：Wdxdydzdt Tdxdydzdt dt tTdxdydzc +?=?2 （f ） 除以dxdydzdt c ，移项后，即得热传导微分方程：c W T c t T =?2(10.5) 对于稳定温度场，有：02=+?WT (10.6)当无热源时：，则：0),(=z y x W 02=?T (10.7)即稳定的无热源温度场满足Laplace 方程。10-2温度场的定解条件为了能够求解热传导微分方程，从而求得温度场，必须已知物体在初瞬时得温度分布，即初始条件；同时还必须已知初瞬时以后弹性体表
9、面与周围介质之间进行热交换的规律，即边界条件。初始条件和边界条件合称为定解条件。（1） 初始条件00(,)(,)t T x y z t T x y z = (10.8) （2） 边界条件第一类边界条件：已知弹性体表面S 上任意一点的瞬时温度，即(,)()S S T x y z t T t = (10.9)式中：为弹性体表面的温度，在最简单的情况下：)(t T S 常数=)(t T S 。 第二类边界条件：已知弹性体表面S 上任意一点的法向热流密度，即(,)(,)n S S q x y z t q x y z t = (10.10)式中：为已知函数，将 (10-2)式用于弹性体表面，上式也可写为
10、：(,)S q x y z t (,)(,)S ST x y z t q x y z t n ?=? (10.11)当为绝热边界时，无热量流过表面S ，此时(,)0S q x y z t =，所以绝热边界条件为：(,)0ST x y z t n ?= （10-12）第三类边界条件：已知弹性体周边介质的温度为，弹性体与其周边介质的热交换符合牛顿冷却定律：a T (,)(,)(,)n a S S ST x y z t q x y z t T x y z t T n ?=?=? （10.13） 上式可以改写为：(,)(,)a ST x y z t T x y z t T n ?+?=? （10.14
11、）10-3 热弹性力学方程现在讨论温度应力计算的第二部分，即根据弹性体内的已知温度改变来决定体内的温度应力。以下仅考虑温度改变作用下的热弹性基本方程。（1） 平衡方程平衡方程仅从静力学的观点出发，研究弹性体内任一微体的平衡，与弹性体温度的改变无关。 与2-4节的推导相同，不考虑体力作用下，平衡方程为：x xy xz yx y yz zx zy z x y z ?= （10.15）000?（2） 几何方程几何方程纯粹从几何的角度出发研究弹性体内位移与应变之间的关系，并不考虑引起位移与应变的原因，与弹性体温度的改变无关，与2-5节的推导相同，几何方程为：（10.16）（3） 物理方程设弹性体内各点
12、的温变为，为了与温度值T 区分，这里将温变用符号（Temperature variation ）表示，即后一瞬时的温度减去前一瞬时的温度。规定升温为正，降温为负。由于温变，弹性体内各点的微小长度，如果不受约束，将发生正应变v T v T v T v T ，其中是弹性体的热膨胀系数，它的单位是温度-1。在各向同性体中，系数为常数，不随方向而变，所以这种正应变在所有各个方向都相同，因而也就不伴随着任何剪应变。弹性体内各点由于温度产生的应变分量为：v z y x T =，0=xy zx yz （g ）但是，由于弹性体所受的外在约束以及体内各部分之间的相互约束，上述的应变并不能自由发生，于是就产生了应
13、力，即所谓的温度应力。这个温度应力又将由于物体的弹性而引起附加的应变，连同式（g ）所示的应变，总的应变分量是：?+=+=+=+?=+?=+?=xy xy zxzx yzyz v y x z z v x z y y v z y x x E E ET ET ET E)1(2)1(2)1(2)(1)(1)(1 （10.17）（4） 边界条件在应力边界S 上，不考虑外荷载情况下，应力边界条件为：（10.18）000x xy xz yx y yz zx zy z S l m n ?=? 在位移边界上，有已知的位移u S (,Tu v w )时，位移边界条件为：u S u u v w w ?=?v （1
14、0.19）10-4 温度应力的等效荷载法以平面问题为例说明温度应力的等效荷载法。假定在如图10-3(a)所示的等厚度薄板及坐标系中，没有体力和面力的作用，但是有温变的作用，而这个温变只是v T v T x 和的函数，不随y z 而变化，且有0=z ，0=yz ，0=zx 。这个问题为温度平面应力问题。 (a) (b)图10-3 (a)温度平面应力问题：等厚度薄板，无体力和面力的作用，有温变的作用v T (b) 平面应力问题：等厚度薄板，有体力和面力的作用。 温度平面应力问题的物理方程为:1()1()x x y y y x xy xy T ET E G ?=?+v v ?=?+?=?（10.20）将几何方程（3.2）代入上式，得到以位移表示的应力分量为：2222E E ()()111E E ()()111E E()2(1+)2(1+)v v x x yv y y xxy xy E T E T u
