连续介质力学几个定律汇总情况

实用标准第二章 连续介质力学的基本定律在第一章中，我们仅考察了连续介质运动的运动学描述，而没有考虑到引起运动和变形的因素。本章我们将引入应力等概念，并给出连续介质力学的基本定律：质量守恒定律、动量平衡定律、动量矩平衡定律、能量守恒定律及熵不等式。2.1 应力矢量与应力张量在物体的运动中，物体的两部分之间或物体与其外界间的力学作用是通过力来描述的。在连续介质力学中我们主要研究三种类型的力：(1)一个物体的两部分之间的接触力；(2)由外界作用于物体边界上的接触力；(3)由外界作用于物体内部点的非接触力(如重力、离心力等)。在另一方面，由于(1)(2)型的力总是通过某一接触面发生作用的，因此通常把作用于单位接触面积上的接触力称为表面力，或简称面力；由于(3)型力作用于物体整个体积内所含的物质点，因此通常把它称为体积力，或简称体力。在连续介质力学中重要的公理之一就是关于接触力形式的柯西假设。柯西假设在运动过程中的时刻t对于任何物质坐标X和与之对应的接触面S上的单位法矢量n，表面力的存在形式为ttX,t,n(2.101)通常，我们规定ttX,t,n指向接触面S的外法向时为正，反之为负(见图2.1).现在不管在X和S面与S'面的曲率相差多少。为了研究物体内部的力学状态，我们把一物体用一假想平面S截断成两部分A和B，如图2.3所示。此时S面就是A和B相互作用的接触面，B部分对A部分一点的作用，便可以用A部分截面上的表面力tn来表征，我们称之为应力矢量。反过来，考虑A部分对B部分作用，按照牛顿的作用与反作用定律可得应力矢量tn。它与tn作用于同一平面上的同一点处，并且大小相等，方向相反。即tn tn (2.102)对于物体内部的一点P，通过它可以有无穷多个方向的截面，而对于不同方向的截面，应力矢量也就不同，这种复杂情况只有引进应力张量的概念才能充分地加以描述。为了刻画一点的应力状态，设想在一点 P的附近任意给定一个单位法矢量为文档实用标准ncos1,cos2,cos3,e1n,e2n,e3n(2.103)的平截面。相应地，过 P点沿活动标架作三个坐标平面。于是它们在物体内截得一个微小四面体，如图 2.4所示。在这个微小四面体的每一个面上，都受有物体的其余部分给它的作用力，不妨设在 ABC上受到的作用力为 tA，在PBC，PCA与PAB上的作用力分别为t1A1、t2A2与t3A3，其中A与Ai分别为各微小平面的面积，作用于微小四面体ABCP上单位质量的体力为b。现在假设对物体的任何部分，特别是对微小四面体 ABCP而言，动量的变化率与作用的合力成正比。虽然这是个很自然且牛顿第二定律更强的新假设 (因为牛顿第二定律只适用于整个物体 )，然而，它却不能用实验直接验证，因为不可能做内部表面接触力的直接测定，这种力的存在与大小只能由其它量的观测推知。描述一点是应力张量，描述通过一点的某一截面是应力矢量。对于微小四面体ABCP，柯西定律给出tAt1A1t2A2t3A3bVtAtiAibVtAtiAcos1VibhtmaVa31hAa(2.104)3其中 为物体的密度，h为P点到ABC面的距离，并且考虑到微小四面体的体积.V1hA(2.105)32.104式也可写成tticos11(2.106)ibhha33当微小四面体体积趋于零时，即A0，h0，则有tticosi(2.107)考虑到2.103式，并令tiTi1e1Ti2e2Ti3e3Tijei(2.108)则式2.107可写成tcositineiTijejnTijeiejnTticosiTijejeinTijeiejnTTn(2.109)文档实用标准当T对称时，则tnTTn(2.110)其中TTijeiej(2.111)称为应力张量，其矩阵形式为T11T12T13TT21T22T23(2.112)T31T32T33如果物体中一点处的应力张量已知，那么由式2.112可以得到通过该点的任何截面上的应力矢量，因此应力张量完全地刻画了物体中一点的应力状态。由Ai面上的应力矢量ti的定义可知，titiX,t，而由式2.108知TijTijX,t，因此式2.109变为tX,t,nnTX,t(2.113)上式就是柯西假设的具体形式，常称之为柯西基本定理。下面我们研究应力张量T的各分量的力学意义。考虑到TijeiTejtiej故知，Tij代表作用于ei方向截面上的应力矢量ti在ej方向上的分量，如图2.5所示。我们从图2.5看到，应力张量T的对角线元素Tijij位于所作用平面的法线方向内，故称之为法向应力分量；应力张量T的非对角线元素Tijij位于所作用的平面内，故称为剪切应力分量。2.2 质量守恒定律物质无论经过怎样形式运动，其总质量是不变的，这就是古典连续介质力学中的最重要规律之一—质量守恒定律。下面我们研究质量守恒定律的数学表达式。设为物体的密度，dV表示物质点的体积，由于在运动过程中质量保持不变，所以D(2.201)dV0Dt展开有DdVDdV0(2.202)DtDt文档实用标准又由式DdVvidVdivvdV(2.203)Dtxi于是式2.202可写成Dvi0(2.204)Dtxi其不变性形式为Ddivv0(2.205)Dt其中Dvi(2.206)Dtxitvt把上式代入式2.204，则得vi0(2.207)t xi其不变性形式为div v 0注明v是张量， 只是一个函数，既不是矢量，又不是张量t(2.208)式2.205和式2.208就是质量守恒定律的数学表达式质量守恒方程，在连续介质力学中常称为连续性方程。在正交曲线坐标系中，利用式:Higigj，连续性方程可写为11v1H2H32v2H1H33v3H1H20(2.209)tH1H2H3在直角坐标系中，连续性方程为vxvyvz0(2.210)txyz在柱面坐标系中，利用第第一部分二章式2.13.03，连续性方程为1rvr1vvz0(2.211)trrrz在球面坐标系中，利用第一部分二章式式2.13.04，连续性方程为1r2vr1sinv1v(2.212)tr2rrsinrsin0连续性方程也可用物质描述法表示。在这种情况下质量定恒定律要求V0X,t0dV0Vx,tdV(2.213)文档实用标准其中V是物质在现时刻所占据的体积，而V0是物质在时刻t0所占据的体积。于是V0X,t0dV0V0xX,t,tJdV0V0X,tJdV0(2.214)因为这个关系式对任意体积V0都必须成立，故得0J(2.215)它表示J与时间无关，即Jconst(2.216)这就是物质形式的连续性方程。2.3 动量平衡定律欧拉把下列关系作为在连续介质中普遍成立的一般性原理：Dm(2.301)fDt它称为欧拉第一运动定律。上式说明任意物体具有的动量的变化率等于作用于该物体上的合力f。设所研究物体在其体积 V上受有连续分布的体力和在其体积的边界面 S上连续分布的接触力 fc，因此物体上所受合力为f fb fc (2.302)其中fbVbdV(2.303)fcStdS(2.304)物体的动量为m V vdV (2.305)DxdVDt于是将式2.302和式2.305代入式2.301则VadVStdSVbdV(2.306)D2x表示x点的加速度。由式2.109，可将上式改写为其中a2DtSnTdSVbdVVadV(2.307)利用高斯公式SnTdSVTdV(2.308)则得文档实用标准S TdV V bdV V adV (2.309)即V T b adV 0 (2.310)考虑到V的任意性，则Tba0(2.311)即divTba(2.312)需要指出的是，这里的散度是对于空间坐标的。上式称为柯西第一运动定律。其指标形式为Tji;ibiai(2.313)展开得T11T21T31b1a1(2.314)x1x2x3T12T22T32b2a2(2.315)x1x2x3T13T23T33b3a3(2.316)x1x2x3特别地，在静止的情况下，物体的加速度为零，则式2.313化为divTb0(2.317)在弹性力学中，上式称为平衡方程。在柱面坐标系中，利用第一部分第二章2.13.4.d可得上式化为Trr1TrTzrTrrTbr0(2.318)rrzrTr1TTzTrTrb0(2.319)rrzrTrz1TzTzzTrzbz0(2.320)rrzr在球面坐标系中，利用第一部分第二章2.13.4.e，则2.317式可化为Trr1Tr1Tr1TcotTrTTbr0(2.321)rrrsinrTr1T1T12TrTrcotTTb0(2.322)rrrsinrTr1T1T12cotTTb0(2.323)rrrsinrTrTr2.4 动量矩平衡定律文档实用标准对于任意物体下列关系式成立：DMx0lx0(2.401)Dt其中Mx0表示物体绕 x0点的动量矩，lx0表示作用于物体上的力对x0点的合力矩。上式称为欧拉第二运动定律。设作用于物体上的力矩只是由体力和接触力引起的 ,故其合力矩为lx0Vxx0bdVSxx0tdS(2.402)而物体的动量矩为Mx0xDxdV(2.403)Vx0Dt将式2.402和式2.403代入式2.401，并考虑到Dxx0Dx(2.404)VdVDtDtD2xDxx0Dxxx0dVVxx0DxDVDtdVVDt2dVDtDtDtVDxDxdV张量本身叉乘是0Vxx0D22xdVVxx0DxDdV质量守恒0DtDtDtDtDtVxx0D22xdV(2.405)Dt可得Vxx0adVVxx0bdVSxx0tdS(2.406)其中aD2x表示x点的加速度。考虑到式2.110和高斯公式，则Dt2V x x0 bdV S x x0 tdS V x x0 adV可知 S x x0 tdS V x x0 ngTdSVxx0bdVSxx0nTdSVxx0adVVxx0badV混合积互换SnTxx0dSVxx0baTxx0dV积分定理Vljkxlx0lbjajekiTijxlx0ljlkekdV张量运算Vljkekxlx0lbjajiTijxlx0ldVVljkekxlx0lbjajTij;ixlx0lTijixlx0ldVV ljkek xl x0l Tij;i bj aj Tijil dV根据平衡方程，红色部分为 0文档实用标准VljkTijilekdVVijkTijekdV0(2.407)考虑到体积V的任意性，得ijkTij0(2.408)因此，Tij必须对称张量，即TijTji(2.409)或TTT(2.410)上式叫做柯西第二运动定律。柯西第二运动定律限定应力张量为对称张量，其中只有六个独立分量。2.5 能量守恒定律在连续介质中，如果只研究力学量的影响，而不考虑热学效应，那么连续介质的能量守恒定律可以直接由运动方程导出。首先，将运动方程TbDv(2.501)Dt点乘速度矢量vvTvbvDv(2.502)Dt在体积V上积分DvVvTdVvbdV(2.503)VvVDt考虑到VvDvVD1vvdVDtDt2VD2vvdVV1vvDdV质量守恒0Dt2DtDV1vvdVDt2DV1v2dVDt2DK(2.504)Dt1上式表示在体积V中的总动能K2的时间变化率。另外，考虑到V2vdV文档实用标准vTvjTij;ivjTij,ivj;iTijTvv:TTvDW:TvD:TW:T反对陈与对称双点乘是0TvD:T(2.505)这里利用了反称张量W与对称张量T之间的双重点积为零的性质。把式2.504和式2.505代回到式2.503中去，则得DKVD:TdVVTvdVVvbdV(2.506)Dt运用高斯公式把上式右边第一体积分化为面积分，并利用柯西假设ttnT，则VTvdVSn添加取掉无影响TvdSStvdS(2.507)将上式代入式2.506，于是我们得到在纯力学作用下的能量方程DKVD:TdVStvdSVbvdV其中D是速度梯度的对称部分Dt(2.508)其中方程左边两项分别表示连续介质的动能和内能 (应力生热)的时间变化率，右边两项分别表示接触力和体力所做的功率。若令 U表示内能，则能量方程5.508也可简洁地写成DK DU DWDt Dt DtDW

(2.509)其中Dt表示接触力和体力的功率，记号D表示这个量不一定能写成某个函数的全微分形式。如果同时考虑机械能和非机械能，那么就必须用能量守恒定律的一般形式。能量守恒定律的一般形式可以表述为：动能加上内能对时间的变化率等于总功率加上在单位时间内供给物体的各种其它形式的能量。这些能量包括热能、化学能、电磁能等等。本书只考虑机械能和热能，于是能量守恒定律就化为著名的热力学第一定律的形式。对于热力连续介质(thermomechanical continua)来说，通常把内能的时间变化率写成DUDudVDtVDtDuVuD是0VdVdVDtDt文档实用标准DudV(2.510)VDt其中u称为比内能，表示每单位质量的内能密度。另外，我们定义矢量f为在单位时间内每单位面积的热通量，函数q为在单位时间内每单位质量的热辐射量，于是物体总热量的增量变化率为DQndSVqdV(2.511)SfDt其中n为物体表面的外法向，热通量矢量f由傅立叶定律给出，即fkT(2.512)这里k为热传导系数，T为温度。于是热力连续介质的能量方程可以写成DK DU DW DQDt Dt Dt Dt或写成积分形式

(2.513)D V1 vvdVDt 2

V DudV StvdS V vbdV Sf ndS V qdVDt(2.514)把上式右边面积分化为体积分后再移到左端，则有V1DvvDudVVDt2Dt(2.515)

Tv vb f qdV高斯公式由于体积V是任意的，故有Dvv1Tvvb1fq(2.516)Dtu2利用式2.505，则上式化为vDvDu1D:TvTvb1fq(2.517)DtDt整理得Du1D:T1fq1vTbDv平衡方程0DtDt(2.518)考虑到运动方程成立，则有Du11fq(2.519)DtD:T或Du1DijTij1fiq(2.520)Dtxi文档实用标准上式表示物体内能的时间变化率等于应力功率和吸收的热量之和 。式2.513、式2.514、和式2.519都是能量守恒定律的表现形式。2.6 状态方程熵定律完整地表征一个热力学统称做是对这个系统状态的描述。用来描述这个状态的物理量称状态参数。状态参数随着时间变化表征一个热力学过程。但是，在一般情况下，这些状态参数并不全是独立的，它们之间存在着某种关系。这种关系就称为状态方程。如果某个状态参数可以通过其它几个状态参数表出，则称它为状态函数。现在，我们考虑一个均匀的热力学系统，它处于平衡状态，即在没有外界影响的条件下，系统的各部分在长时间内不发生任何变化。描述这样一个热力学系统的状态参数为：几何参数 V(体积)、力学参数p(压力)及热力学参数T(温度)。联系这三个量的关系的状态方程可写成Fp,V,T0(2.601)这里需要指出的是，对于一定的物质来说，状态方程是普遍适用的，也就是说，构成热力学系统的物质一经选定，状态方程的具体形式也就确定了。例如对于完全气体而言，状态方程的具体形式可写成m(2.602)pVR0TM其中m为气体的质量，M为分子量，R0是克分子气体常数。在上一节我们曾叙述过热力学第一定律，它公设机械能和热能可以互相转换，但是，只根据热力学第一定律还不能判定这种转换过程是否可逆。事实上，所有的真实过程都是不可逆的，但可逆过程却是一个非常有用的假设，因为在许多情况下，能量耗损是可以忽略不计的。可逆性判据由热力学第二定律给出。热力学第二定律公设存在两个独立状态函数：绝对温度 T和熵S。它们有如下性质：绝对温度T为一正量，它仅仅是经验温度(即我们通常见到的温度)的函数，熵S和体积V一样，是一个广延量，而温度是与熵相对应的强度量，正如压强是与体积相对应的强度量一样。一个物体的强度量代表物质的内在性质，与物体的质量大小无关，而一个物体的广延量则可分解为物体上各个子部分上的广延量之和。因此，一连续介质的总熵S可写成下列形式：SVsdV(2.603)这里s表示连续介质中的熵密度，即每单位质量中的熵。文档实用标准一个系统的熵既可由于与外界相互作用而发生改变，也可由于系统内部发生变化而改变，因此dsdsedsi(2.604)这里ds是熵密度的增量,dse是由于与外部相互作用而引起的熵密度增量。dsi是由于系统内部发生变化而引起的熵密度的增量。dsi决不能为负值。它在可逆过程中为零，在不可逆过程中为正，即dsi0(不可逆过程)(2.605)dsi0(可逆过程)(2.606)在可逆过程中，如果令dqR表示供给系统的每单位质量的热量，则dse可表示为dsedqR(可逆过程)(2.607)T按照热力学第二定律，在连续介质所占据的物理空间中总熵的时间变率不小于通过连续介质表面流入的熵与连续体内部源产生的熵之和。在数学上，这个熵原理可以以积分形式表示为dfn(2.608)dtVsdVVedVSTdS称之为克劳修斯—杜姆不等式，其中e为单位质量中的局部熵源。上式中的等号成立时表示可逆过程，不等号成立时代表不可逆过程。利用质量守恒定律dVsdVVdsdVVSddVdtdtdtdsdVVdt和高斯公式V f ndST

V

dVT考虑到体积V的任意性，则由式2.608可得克劳修斯—杜姆不等式的微分形式ds1f(2.609)e0dtT2.7 主应力最大剪应力nT表示物体中一点周围不同方向上的应力矢量公式，当应力张量已知时，在给定的任何一个方向n上的应力矢量就由tnT给出。下面，我们将要文档实用标准讨论的问题是，对于某给定点来说，在什么方向上法向应力 Tn取驻值。这个问题归结为在n为单位矢量的条件下，即n2nnn12n22n32nknk1(2.701)时，求Tn的条件极值问题。运用大家所熟知的拉格朗日乘子法，有Tnf0(2.702)nini其中f为约束条件fnnn11kk0(2.703)nn考虑到TijTji，则由式2.110可得TnntnTnnkekTpqepeqnlelkpnkTpqnlqlnpTpqnq(2.704)将上式代入式2.702，则TnfninininpTpqnqninknk1npTpqnqnpTpqnqnknknini2nipiTpqnqnpTpqqi2kink2Tiqnqni0(2.705)或写成不变性形式，即Tnn(2.706)或TIn0(2.707)写成展开形式，则为T11n1T12n2T13n30T21n1T22n2T23n30T31n1T32n2T33n30(2.708)上列方程中n具有非零解的充分必要条件是它的系数行列式为零，即Tijij0(2.709)或3I12I2I30(2.710)其中文档实用标准I1T11T22T33TiitrT(2.711)I2T11T12T11T13T22T23T21T22T31T33T32T331TiiTjjTij21trT2trT2(2.712)2T11T12T13I3T12T22T23TijdetT(2.713)T31T32T33这里I1，I2，I3是应力张量T的三个主不变量，分别称为第一、第二、第三应力不变量。方程的解1，2，3为特征值，n1，n2，n3为特征矢量。其中若i j，则ni nj。事实上，在ni方向上法向应力值就是ni所对应的特征值。将式2.706与ni点乘，得i inini niTni tini (2.714)则i就是ni方向上的应力，称为主应力，而ni称为主方向，主方向所确定的平面称为主平面。若ni和nj不两个不同的主方向ij，则在ni面上nj方向的剪应力Tij为TijniTnjininj0(2.715)故主应力平面上的剪应力为零。若以(n1，n2，n3)为坐标单位基矢量，并令Ti，则应力张量矩阵具有下列形式：T100T0T20(2.716)00T3即T Tniini (2.717)现在我们来讨论最大剪切应力问题。为了计算方便，不妨将坐标系选取在主方向上，即取(e1，e2，e3)为主方向。设n是通过物体内一点的某一平面的单位法向矢量，则nn1e1n2e2n3e3nkek(2.718)作用于该平面的应力矢量分量为tnTneTee(2.719)kkiiinTenTei(2.720)k