第二章 连续介质力学的基本定律
在第一章中,我们仅考察了连续介质运动的运动学描述,而没有考虑到引起运动和变形的因素。本章我们将引入应力等概念,并给出连续介质力学的基本定律:质量守恒定律、动量平衡定律、动量矩平衡定律、能量守恒定律及熵不等式。
2.1 应力矢量与应力张量
在物体的运动中,物体的两部分之间或物体与其外界间的力学作用是通过力来描述的。在连续介质力学中我们主要研究三种类型的力:(1)一个物体的两部分之间的接触力;(2)由外界作用于物体边界上的接触力;(3)由外界作用于物体内部点的非接触力(如重力、离心力等) 。在另一方面,由于(1)(2)型的力总是通过某一接触面发生作用的,因此通常把作用于单位接触面积上的接触力称为表面力,或简称面力;由于(3)型力作用于物体整个体积内所含的物质点,因此通常把它称为体积力,或简称体力。
在连续介质力学中重要的公理之一就是关于接触力形式的柯西假设。柯西假设在运动过程中的时刻t 对于任何物质坐标X 和与之对应的接触面S 上的单位法矢量n ,表面力的存在形式为
t =t (X , t , n ) (2.101) 通常,我们规定t =t (X , t , n )指向接触面S 的外法向时为正,反之为负(见图2.1). 现在不管在X 和S 面与S' 面的曲率相差多少。
为了研究物体内部的力学状态,我们把一物体用一假想平面S 截断成两部分A 和B ,如图2.3所示。此时S 面就是A 和B 相互作用的接触面,B 部分对A 部分一点的作用,便可以用A 部分截面上的表面力t n 来表征,我们称之为应力矢量。反过来,考虑A 部分对B 部分作用,按照牛顿的作用与反作用定律可得应力矢量t -n 。它与t n 作用于同一平面上的同一点处,并且大小相等,方向相反。即
t n =-t n (2.102) 对于物体内部的一点P ,通过它可以有无穷多个方向的截面,而对于不同方向的截面,应力矢量也就不同,这种复杂情况只有引进应力张量的概念才能充分地加以描述。为了刻画一点的应力状态,设想在一点P 的附近任意给定一个单位法矢量为
n =(cos α1, cos α2, cos α3, )
=(e 1⋅n , e 2⋅n , e 3⋅n ) (2.103) 的平截面。相应地,过P 点沿活动标架作三个坐标平面。于是它们在物体内截得一个微小四面体,如图2.4所示。在这个微小四面体的每一个面上,都受有物体的其余部分给它的作用力,不妨设在ABC 上受到的作用力为t ∆A ,在PBC ,PCA 与PAB 上的作用力分别为-t 1∆A 1、-t 2∆A 2与-t 3∆A 3,其中∆A 与∆A i 分别为各微小平面的面积,作用于微小四面体ABCP 上单位质量的体力为b 。
现在假设对物体的任何部分,特别是对微小四面体ABCP 而言,动量的变化率与作用的合力成正比。虽然这是个很自然且牛顿第二定律更强的新假设(因为牛顿第二定律只适用于整个物体) ,然而,它却不能用实验直接验证,因为不可
能做内部表面接触力的直接测定,这种力的存在与大小只能由其它量的观测推知。描述一点是应力张量,描述通过一点的某一截面是应力矢量。
对于微小四面体ABCP ,柯西定律给出 t ∆A -t 1∆A 1-t 2∆A 2-t 3∆A 3+ρb ∆V =t ∆A -t i ∆A i +ρb ∆V
1
=t ∆A -t i ∆A cos αi +ρbh ∆V
3
=tma =ρ∆Va
1
=ρh ∆Aa (2.104)
3
其中ρ为物体的密度,h 为P 点到ABC 面的距离,并且考虑到微小四面体的体积.
1∆V =h ∆A (2.105)
3
2.104式也可写成
11
t -t cos α+ρbh =ρha (2.106) i i
33
当微小四面体体积趋于零时,即∆A →0,∆h →0,则有
t =t i cos αi (2.107) 考虑到2.103式,并令
t i =T i 1e 1+T i 2e 2+T i 3e 3
=T ij e i (2.108) 则式2.107可写成
t =cos αi t i =(n ⋅e i )(T ij e j )
=n ⋅(T ij e i e j )=n ⋅T
=(T ij e i e j )⋅n =T T ⋅n (2.109)
=t i cos αi =(T ij e j )(e i ⋅n )
当T 对称时,则
t =n ⋅T =T ⋅n (2.110) 其中
T =T ij e i e j (2.111) 称为应力张量,其矩阵形式为
⎡T 11T 12T 13⎤
⎥ (2.112) T T T [T ]=⎢212223⎢⎥
⎢⎣T 31T 32T 33⎥⎦ 如果物体中一点处的应力张量已知,那么由式2.112可以得到通过该点的任
何截面上的应力矢量,因此应力张量完全地刻画了物体中一点的应力状态。 由A i 面上的应力矢量t i 的定义可知,t i =t i (X , t ),而由式2.108知
T ij =T ij (X , t ),因此式2.109变为
t (X , t , n )=n ⋅T (X , t ) (2.113) 上式就是柯西假设的具体形式,常称之为柯西基本定理。 下面我们研究应力张量T 的各分量的力学意义。考虑到
T ij =e i ⋅T ⋅e j =t i ⋅e j
故知,T ij 代表作用于e i 方向截面上的应力矢量t i 在e j 方向上的分量,如图2.5所示。
我们从图2.5看到,应力张量T 的对角线元素T ij (i =j )位于所作用平面的法线方向内,故称之为法向应力分量;应力张量T 的非对角线元素T ij (i ≠j )位于所作用的平面内,故称为剪切应力分量。
2.2 质量守恒定律
物质无论经过怎样形式运动,其总质量是不变的,这就是古典连续介质力学中的最重要规律之一—质量守恒定律。下面我们研究质量守恒定律的数学表达式。
设ρ为物体的密度,dV 表示物质点的体积,由于在运动过程中质量保持不变,所以
D
(ρdV )=0 (2.201) Dt
展开有
D
(dV )+D ρdV =0 (2.202) ρDt Dt
又由式
∂v D
(dV )=i dV =(divv )dV (2.203)
Dt ∂x i
于是式2.202可写成
D ρ∂v
+ρi =0 (2.204) Dt ∂x i
其不变性形式为
D ρ
+ρdivv =0 (2.205) Dt
其中
D ρ∂ρ∂ρ=+v (2.206) i Dt ∂t ∂x i
∂ρ
+v ⋅∇ρ =∂t
把上式代入式2.204,则得
∂ρ∂(ρv i ) +=0 (2.207) ∂t ∂x i
其不变性形式为 ∂ρ
+div (ρv )=0注明v 是张量,ρ只是一个函数,既不是矢量,又不是张量 ∂t
(2.208)
式2.205和式2.208就是质量守恒定律的数学表达式质量守恒方程,在连续介质力学中常称为连续性方程。
在正交曲线坐标系中,利用式:H i =
g i ⋅g j ,连续性方程可写为
∂ρ1
[∂1(ρv 1H 2H 3)+∂2(ρv 2H 1H 3)+∂3(ρv 3H 1H 2)]=0 (2.209) +
∂t H 1H 2H 3
在直角坐标系中,连续性方程为
∂ρ∂(ρv x )∂(ρv y )∂(ρv z ) +++=0 (2.210) ∂t ∂x ∂y ∂z
在柱面坐标系中,利用第第一部分二章式2.13.03,连续性方程为
∂ρ1∂(ρrv r )1∂(ρv θ)∂(ρv z )+++=0 (2.211)
∂t r ∂r r ∂θ∂z
在球面坐标系中,利用第一部分二章式式2.13.04,连续性方程为
∂ρ1∂ρr 2v r 1∂(ρsin θv θ)1∂(ρv ϕ) +++=0 (2.212) ∂t r 2∂r r sin θ∂θr sin θ∂ϕ
连续性方程也可用物质描述法表示。在这种情况下质量定恒定律要求
⎰V 0ρ(X , t 0)dV 0=⎰V ρ(x , t )dV (2.213)
()
其中V 是物质在现时刻所占据的体积,而V 0是物质在时刻t 0所占据的体积。于是 ⎰V 0ρ(X , t 0)dV 0=⎰V 0ρ[x (X , t ), t ]JdV 0
=⎰V 0ρ(X , t )JdV 0 (2.214) 因为这个关系式对任意体积V 0都必须成立,故得
ρ0=ρJ (2.215) 它表示ρJ 与时间无关,即
ρJ =const (2.216) 这就是物质形式的连续性方程。
2.3 动量平衡定律
欧拉把下列关系作为在连续介质中普遍成立的一般性原理:
Dm
=f (2.301) Dt
它称为欧拉第一运动定律。上式说明任意物体具有的动量的变化率等于作用于该物体上的合力f 。
设所研究物体在其体积V 上受有连续分布的体力和在其体积的边界面S 上连续分布的接触力f c ,因此物体上所受合力为
f =f b +f c (2.302) 其中
f b =⎰V ρbdV (2.303) f c =⎰S tdS (2.304) 物体的动量为
m =⎰V ρvdV (2.305) =⎰V ρ
Dx
dV Dt
于是将式2.302和式2.305代入式2.301则
⎰V ρadV =⎰S tdS +⎰V ρbdV (2.306)
D 2x
其中a =表示x 点的加速度。由式2.109,可将上式改写为
Dt 2
⎰S n ⋅TdS +⎰V ρbdV =⎰V ρadV (2.307)
利用高斯公式
⎰S n ⋅TdS =⎰V ∇⋅TdV (2.308) 则得
⎰S ∇⋅TdV +⎰V ρbdV =⎰V ρadV (2.309) 即
⎰V (∇⋅T +ρb -ρa )dV =0 (2.310) 考虑到V 的任意性,则
∇⋅T +ρb -ρa =0 (2.311) 即
divT +ρb =ρa (2.312) 需要指出的是,这里的散度是对于空间坐标的。上式称为柯西第一运动定律。其指标形式为
T ji ; i +ρb i =ρa i (2.313) 展开得
∂T ∂T ∂T
11+21+31+ρb 1=ρa 1 (2.314)
∂x 1∂x 2∂x 3∂T ∂T ∂T
12+22+32+ρb 2=ρa 2 (2.315)
∂x 1∂x 2∂x 3∂T ∂T ∂T
13+23+33+ρb 3=ρa 3 (2.316)
∂x 1∂x 2∂x 3
特别地,在静止的情况下,物体的加速度为零,则式2.313化为
divT +ρb =0 (2.317) 在弹性力学中,上式称为平衡方程。
在柱面坐标系中,利用第一部分第二章2.13.4.d 可得上式化为
∂T rr 1∂T θr ∂T zr T rr -T θθ
++++ρb r =0 (2.318)
∂r r ∂θ∂z r ∂T 1∂T θθ∂T z θT r θ-T θr
+++ρb θ=0 (2.319) r θ+
∂r r ∂θ∂z r ∂T 1∂T θz ∂T zz T rz
+++ρb z =0 (2.320) rz +
∂r r ∂θ∂z r
在球面坐标系中,利用第一部分第二章2.13.4.e ,则2.317式可化为
∂T rr 1∂T θr 1∂T ϕr 1 +++(2T rr +cot θT θr -T θθ-T ϕϕ)+ρb r =0 (2.321)
∂r r ∂θr sin θ∂ϕr ∂T 1∂T θθ1∂T θθ1 r θ+++2T r θ+T θr +cot (T θθ-T ϕϕ)+ρb θ=0 (2.322)
∂r r ∂θr sin θ∂ϕr
[]
∂T r ϕ1∂T θϕ1∂T ϕϕ1
+++[2T r ϕ+T ϕr +cot θ(T θϕ-T ϕθ)]+ρb ϕ=0(2.323) ∂r r ∂θr sin θ∂ϕr
2.4 动量矩平衡定律
对于任意物体下列关系式成立:
DM x 0
=l x 0 (2.401) Dt
其中M x 0表示物体绕x 0点的动量矩,l x 0表示作用于物体上的力对x 0点的合力矩。上式称为欧拉第二运动定律。
设作用于物体上的力矩只是由体力和接触力引起的, 故其合力矩为
l x 0=⎰V ρ(x -x 0)⨯bdV +⎰S (x -x 0)⨯td S (2.402) 而物体的动量矩为 M x 0=⎰V ρ(x -x 0)⨯
Dx
dV (2.403) Dt
将式2.402和式2.403代入式2.401,并考虑到
D Dx
ρx -x ⨯dV (2.404) ()V 0⎰Dt Dt D (x -x 0)Dx D 2x Dx D
⨯dV +⎰V ρ(x -x 0)⨯2dV +⎰V (x -x 0)⨯ =⎰V ρ(ρdV ) Dt Dt Dt Dt Dt
Dx Dx D 2x Dx D =⎰V ρ⨯dV 张量本身叉乘是0+⎰V ρ(x -x 0)⨯2dV +⎰V (x -x 0)⨯(ρdV )质量守恒0
Dt Dt Dt D t Dt
D 2x
=⎰V ρ(x -x 0)⨯2dV (2.405)
Dt
可得
⎰V ρ(x -x 0)⨯adV =⎰V ρ(x -x 0)⨯bdV +⎰S (x -x 0)⨯td S (2.406)
D 2x
其中a =表示x 点的加速度。考虑到式2.110和高斯公式,则
Dt 2
⎰V ρ(x -x 0)⨯bdV +⎰S (x -x 0)⨯td S -⎰V ρ(x -x 0)⨯ad V 可知⎰S (x -x 0)⨯td S =⎰V (x -x 0)⨯(n T )d S
=⎰V ρ(x -x 0)⨯bdV +⎰S (x -x 0)⨯(n ⋅T )dS -⎰V ρ(x -x 0)⨯adV =⎰V ρ(x -x 0)⨯(b -a )dV +⎰S n ⋅T ⨯(x -x 0)dS 混合积互换 =⎰V ρ(x -x 0)⨯(b -a )+∇⋅⎡⎣T ⨯(x -x 0)⎤⎦dV 积分定理
=⎰V ρεljk (x l -x 0l )⨯(b j -a j )e k +∂i ⋅⎡⎣T ij ⨯(x l -x 0l )εjlk e k ⎤⎦dV 张量运算 =⎰V εljk e k {ρ(x l -x 0l )⨯(b j -a j )+∂i [T ij (x l -x 0l )]}dV
=⎰V εljk e k [ρ(x l -x 0l )⨯(b j -a j )+T ij ; i (x l -x 0l )+T ij ∂i (x l -x 0l )]dV
{}
{}
⎤=⎰V εljk e k ⎡⎣(x l -x 0l )⨯(T ij ; i +ρb j -ρa j )+T i j δil ⎦dV 根据平衡方程,红色部分为0
=⎰V εljk T ij δil e k dV
=
⎰
V
εijk T ij e k dV
=0 (2.407) 考虑到体积V 的任意性,得
εijk T ij =0 (2.408) 因此,T ij 必须对称张量,即
T ij =T ji (2.409) 或
T
T =T (2.410) 上式叫做柯西第二运动定律。柯西第二运动定律限定应力张量为对称张量,其中只有六个独立分量。
2.5 能量守恒定律
在连续介质中,如果只研究力学量的影响,而不考虑热学效应,那么连续介质的能量守恒定律可以直接由运动方程导出。首先,将运动方程
Dv
∇⋅T +ρb =ρ (2.501)
Dt
点乘速度矢量v
Dv
v ⋅(∇⋅T )+ρv ⋅b =ρv ⋅ (2.502)
Dt
在体积V 上积分
Dv
=V v ⋅(∇⋅T )dV +⎰V ρv ⋅bdV (2.503) ⎰V ρv ⋅
Dt ⎰
考虑到
Dv D ⎛1⎫
⎰V ρv ⋅=⎰V ρ v ⋅v ⎪dV
Dt Dt ⎝2⎭
D ⎛ρ1D ⎫
=⎰V v ⋅vdV -v ⋅v (ρdV )质量守恒0 ⎪⎰V
D t ⎝22Dt ⎭D 1
ρv ⋅vdV =V
Dt ⎰2D 12
ρv dV =V ⎰Dt 2
DK
= (2.504)
Dt
1
上式表示在体积V 中的总动能K =⎰V ρv 2dV 的时间变化率。另外,考虑到
2
v ⋅(∇⋅T )=v j T ij ; i =v j T ij
()
, i
-v j ; i T ij
=∇⋅(T ⋅v )-(∇v ):T
=∇⋅(T ⋅v )-(D +W ):T
=∇⋅(T ⋅v )-D :T -W :T 反对陈与对称双点乘是0
=∇⋅(T ⋅v )-D :T (2.505) 这里利用了反称张量W 与对称张量T 之间的双重点积为零的性质。 把式2.504和式2.505代回到式2.503中去,则得
DK
+V D :TdV =⎰V ∇⋅(T ⋅v )dV +⎰V ρv ⋅bdV (2.506)
Dt ⎰
运用高斯公式把上式右边第一体积分化为面积分,并利用柯西假设t =t =n ⋅T ,则
⎰V ∇⋅(T ⋅v )dV =⎰S n ⋅(T ⋅v )dS 添加取掉()无影响
=⎰S t ⋅vdS (2.507)
将上式代入式2.506,于是我们得到在纯力学作用下的能量方程
DK
+⎰V D :TdV =⎰S t ⋅vdS +⎰V ρb ⋅vdV 其中D 是速度梯度的对称部分 Dt (2.508)
其中方程左边两项分别表示连续介质的动能和内能(应力生热) 的时间变化率,右边两项分别表示接触力和体力所做的功率。若令U 表示内能,则能量方程5.508也可简洁地写成
DK DU += (2.509) Dt Dt Dt 其中表示接触力和体力的功率,记号表示这个量不一定能写成某个函数
Dt
的全微分形式。
如果同时考虑机械能和非机械能,那么就必须用能量守恒定律的一般形式。能量守恒定律的一般形式可以表述为:动能加上内能对时间的变化率等于总功率加上在单位时间内供给物体的各种其它形式的能量。这些能量包括热能、化学能、电磁能等等。本书只考虑机械能和热能,于是能量守恒定律就化为著名的热力学第一定律的形式。
对于热力连续介质(thermomechanical continua)来说,通常把内能的时间变化率写成
DU D
= V ρudV ⎰Dt Dt
Du D dV +⎰V u =⎰V ρ(ρdV )是0 Dt Dt Du
dV (2.510) =⎰V ρDt
其中u 称为比内能,表示每单位质量的内能密度。另外,我们定义矢量f 为在单位时间内每单位面积的热通量,函数q 为在单位时间内每单位质量的热辐射量,于是物体总热量的增量变化率为
Q
=-⎰S f ⋅ndS +⎰V ρqdV (2.511) Dt
其中n 为物体表面的外法向,热通量矢量f 由傅立叶定律给出,即
f =k ∇T (2.512)
这里k 为热传导系数,T 为温度。
于是热力连续介质的能量方程可以写成
DK DU +=+ (2.513) Dt Dt Dt Dt 或写成积分形式
D 1Du
ρv ⋅vdV +ρ V ⎰V Dt dV =⎰S t ⋅vdS +⎰V ρv ⋅bdV -⎰S f ⋅ndS +⎰V ρqdV Dt ⎰2
(2.514) 把上式右边面积分化为体积分后再移到左端,则有
⎡1D (v ⋅v )Du ⎤
⎰V ⎢ρ+ρ∇⋅(T ⋅v )+ρv ⋅b -∇⋅f +ρq ⎤dV 高斯公式 ⎥dV =⎰V ⎡⎣⎦Dt Dt ⎦⎣2(2.515)
由于体积V 是任意的,故有
D ⎛v ⋅v 1⎫1 +u ⎪=∇⋅(T ⋅v )+v ⋅b -∇⋅f +q (2.516)
Dt ⎝2ρ⎭ρ
利用式2.505,则上式化为
Dv Du 11
v ⋅+=[D :T +v ⋅(∇⋅T )]+v ⋅b -∇⋅f +q (2.517)
Dt Dt ρρ
整理得
Du 111⎛Dv ⎫ =D :T -∇⋅f +q +v ⋅ ∇⋅T +ρb -ρ⎪平衡方程0
Dt ρρρ⎝D t ⎭(2.518)
考虑到运动方程成立,则有
Du 11
=D :T -∇⋅f +q (2.519)
Dt ρρ或
Du 11∂f i
=D ij T ij -+q (2.520)
Dt ρρ∂x i
上式表示物体内能的时间变化率等于应力功率和吸收的热量之和。 式2.513、式2.514、和式2.519都是能量守恒定律的表现形式。
2.6 状态方程熵定律
完整地表征一个热力学统称做是对这个系统状态的描述。用来描述这个状态的物理量称状态参数。状态参数随着时间变化表征一个热力学过程。但是,在一般情况下,这些状态参数并不全是独立的,它们之间存在着某种关系。这种关系就称为状态方程。如果某个状态参数可以通过其它几个状态参数表出,则称它为状态函数。
现在,我们考虑一个均匀的热力学系统,它处于平衡状态,即在没有外界影响的条件下,系统的各部分在长时间内不发生任何变化。描述这样一个热力学系统的状态参数为:几何参数V(体积) 、力学参数p(压力) 及热力学参数T(温度) 。联系这三个量的关系的状态方程可写成
F (p , V , T )=0 (2.601)
这里需要指出的是,对于一定的物质来说,状态方程是普遍适用的,也就是说,构成热力学系统的物质一经选定,状态方程的具体形式也就确定了。 例如对于完全气体而言,状态方程的具体形式可写成
m pV =R 0T (2.602)
M
其中m 为气体的质量,M 为分子量,R 0是克分子气体常数。
在上一节我们曾叙述过热力学第一定律,它公设机械能和热能可以互相转换,但是,只根据热力学第一定律还不能判定这种转换过程是否可逆。事实上,所有的真实过程都是不可逆的,但可逆过程却是一个非常有用的假设,因为在许多情况下,能量耗损是可以忽略不计的。可逆性判据由热力学第二定律给出。
热力学第二定律公设存在两个独立状态函数:绝对温度T 和熵S 。它们有如下性质:绝对温度T 为一正量,它仅仅是经验温度θ(即我们通常见到的温度) 的函数,熵S 和体积V 一样,是一个广延量,而温度是与熵相对应的强度量,正如压强是与体积相对应的强度量一样。一个物体的强度量代表物质的内在性质,与物体的质量大小无关,而一个物体的广延量则可分解为物体上各个子部分上的广延量之和。因此,一连续介质的总熵S 可写成下列形式:
S =⎰V ρsdV (2.603) 这里s 表示连续介质中的熵密度,即每单位质量中的熵。
一个系统的熵既可由于与外界相互作用而发生改变,也可由于系统内部发生变化而改变,因此 ds =ds (e )+ds (i ) (2.604) 这里ds 是熵密度的增量, ds (e )是由于与外部相互作用而引起的熵密度增量。ds (i )是由于系统内部发生变化而引起的熵密度的增量。ds (i )决不能为负值。它在可逆过程中为零,在不可逆过程中为正,即
ds (i )>0 (不可逆过程) (2.605) ds (i )=0 (可逆过程) (2.606) 在可逆过程中,如果令dq (R )表示供给系统的每单位质量的热量,则ds (e )可表示为 ds (e )=
dq (R )
T
按照热力学第二定律,在连续介质所占据的物理空间中总熵的时间变率不小于通过连续介质表面流入的熵与连续体内部源产生的熵之和。在数学上,这个熵原理可以以积分形式表示为
f ⋅n d
(2.608) ⎰V ρsdV ≥⎰V ρedV -⎰S
dt T
称之为克劳修斯—杜姆不等式,其中e 为单位质量中的局部熵源。上式中的等号成立时表示可逆过程,不等号成立时代表不可逆过程。 利用质量守恒定律
d ds d
ρdV +⎰V S (ρdV ) ⎰V ρsdV =⎰V
dt dt dt
ds
ρdV =⎰V
dt
(可逆过程) (2.607)
和高斯公式 f ⋅n ⎛f ⎫ ⎰V dS =⎰V ∇⋅ ⎪dV T ⎝T ⎭
考虑到体积V 的任意性,则由式2.608可得克劳修斯—杜姆不等式的微分形式 ds 1⎛f ⎫ -e -∇⋅ ⎪≥0 (2.609) dt ρ⎝T ⎭
2.7 主应力最大剪应力
t =n ⋅T 表示物体中一点周围不同方向上的应力矢量公式,当应力张量已知时,在给定的任何一个方向n 上的应力矢量就由t =n ⋅T 给出。下面,我们将要讨论的问题是,对于某给定点来说,在什么方向上法向应力T n 取驻值。这个问题归结为在n 为单位矢量的条件下,即 2222 n =n ⋅n =n 1+n 2+n 3
=n k n k =1 (2.701) 时,求T n 的条件极值问题。运用大家所熟知的拉格朗日乘子法,有 ∂T n ∂f -λ=0 (2.702) ∂n i ∂n i
其中f 为约束条件
f (n )=n ⋅n -1=-1n k n k =0 (2.703) 考虑到T ij =T ji ,则由式2.110可得
T n =n ⋅t =n ⋅T ⋅n
=(n k e k )⋅(T pq e p e q )⋅(n l e l )
=δkp n k T pq n l δql
=n p T pq n q (2.704) 将上式代入式2.702,则 ∂T n ∂f -λ ∂n i ∂n i
∂(n p T pq n q )-λ∂(n k n k -1) ∂n i ∂n i
∂n p ∂n q ∂n k =T n +n T -2λn k pq q p pq ∂n i ∂n i ∂n i
=δpi T pq n q +n p T pq δqi -2λδki n k
=2(T iq n q -λn i )=0 (2.705) =-或写成不变性形式,即
T ⋅n =λn (2.706) 或
(T -λI )⋅n =0 (2.707) 写成展开形式,则为
(T 11-λ)n 1+T 12n 2+T 13n 3=0
T 21n 1+(T 22-λ)n 2+T 23n 3=0
T 31n 1+T 32n 2+(T 33-λ)n 3=0 (2.708) 上列方程中n 具有非零解的充分必要条件是它的系数行列式为零,即 T ij -λδij =0 (2.709) 或
32 λ-I 1λ+I 2λ-I 3=0 (2.710)
其中
I 1=T 11+T 22+T 33=T ii =trT (2.711) T 11T 12T 11T 13T 22T 23++ I 2= T 21T 22T 31T 33T 32T 33
1 =(T ii T jj -T ij ) 2
12 =(trT )-trT 2 (2.712) 2
T 11T 12T 13[]
I 3=T 12T 22T 23=T ij =det T (2.713)
T 31T 32T 33
这里I 1,I 2,I 3是应力张量T 的三个主不变量,分别称为第一、第二、第三应力
123不变量。方程的解λ1,λ2,λ3为特征值,n ,n ,n 为特征矢量。其中若
λi ≠λj ,则n i ⊥n j 。
i i i 事实上,在n 方向上法向应力值就是n 所对应的特征值。将式2.706与n 点
乘,得 i i i i i λi =λi n ⋅n =n ⋅T ⋅n =t i ⋅n (2.714)
i i 则λi 就是n 方向上的应力,称为主应力,而n 称为主方向,主方向所确定的平面
称为主平面。
i i j j 若n 和n 不两个不同的主方向(i ≠j ),则在n 面上n 方向的剪应力T ij 为 i j i j T ij =n ⋅T ⋅n =λi n n =0 (2.715)
123故主应力平面上的剪应力为零。若以(n ,n ,n ) 为坐标单位基矢量,并令
λi =T i ,则应力张量矩阵具有下列形式:
⎡T 100⎤⎥ (2.716) 0T 0 [T ]=⎢
2⎢⎥⎢⎣00T 3⎥⎦
即 i i T =Tn i n (2.717)
现在我们来讨论最大剪切应力问题。为了计算方便,不妨将坐标系选取在主方向上,即取(e 1,e 2,e 3) 为主方向。设n 是通过物体内一点的某一平面的单位法向矢量,则
n =n 1e 1+n 2e 2+n 3e 3=n k e k (2.718) 作用于该平面的应力矢量分量为 t =n ⋅T =n k e k ⋅Te i e i (2.719) n k T δki e i =n i Te i (2.720) 在该平面上的法向应力为
T n =t ⋅n =n i ⋅Te i i ⋅n j e j =n i n j T i δij =n i 2T i (2.721) 若以T S 表示该平面的总剪应力的大小(如图2.6) ,则 22 T S =t -T n (2.722)
即
22 T S =t ⋅t -T n 2 =n i T i e i ⋅n j T j e j -(n i n i T i ) 2
=n i n i T i 2-n i n i T i (2.723) =n i 2T i 2-T n 2
2 我们仍运用拉格朗日乘子法计算T S 的驻值,考虑到n 为单位矢量,令
f (n )=n ⋅n -1=n i n i -1=0 (2.724) ()2则 ∂T S 2∂f +λ=0 (2.725) ∂n k ∂n k
其中 ∂n ∂T S 2∂T n 2=2n T -2T i i n ∂n k ∂n k ∂n k
⎛∂n i ⎫2 =2nT δ-2T 2nT ⎪ i ik n i ∂n k ⎭⎝
=2n i T i 2δik -4T n (n i T i δik ) =2n k T i 2-4n k T n T i =2n k (T i 2-2TT i n ) (2.726) ∂f ∂(n i n i -1) =∂n k ∂n k
∂n i =2n =2n i δik =2n k (2.727) i ∂n k
于是 2n k (T i 2-2TT i n +λ)=0 (2.729) 即
2n 1T 12-2T n T 1+λ=0 (
2n (T
2n (T 2
3-2T n T 3
利用条件2.724,则方程组2.730显然有一组解
n 1=±1,n 2=n 3=0
n 2=±1,n 1=n 3=0
n 3=±1,n 1=n 2=0 (2.731)
2但是这组解所确定的平面就是主平面,而在主平面上T S =0,这不是我们所要
求的解。 2223-2T n T 2)+λ)=0 +λ)=0 (2.730)
假定在式2.730中n 1≠0,n 2≠0,n 3=0,则
2 T 1-2T n T 1+λ=0
2 T 2-2T n T 2+λ=0 (2.732) 将上列两式相减,则有
T 12-T 22-2T n (T 1-T 2)=0 (2.733) 故得 T 1-T 2 (2.734) 222把它代入式2.721中并与n 1+n 2=1联立,则可解得 11 n 1=±,n 2=±,n 3=0 (2.735) 2
这时n 方向与主方向e 2,e 3成45度角。
同样,若设n 1≠0,n 2=0,n 3≠0和n 1=0,n 2≠0,n 3≠0则对应的n 值分别为 11 n 1=±,n 2=0,n 3=± (2.736) 22
和
11n =0n =±n =± 1,2,3 (2.737) 22
考虑到上列三组驻值,则
1(e 1+e 2)时,T S =±1(T 1-T 2) (2.738) 当n =±22
1(e 3+e 1)时,T S =±1(T 3-T 1) (2.739) 当n =±22
1(e 2+e 1)时,T S =±1(T 2-T 3) (2.740) 当n =±22
因此,剪切应力的最大值由下列三个值中的最大值给出 T 1-T 2T 3-T 1T 2-T 3 ,, (2.741) 222
或 (T )-(T n )min (T S )max =n max (2.742) 2
T n =
第二章 连续介质力学的基本定律
在第一章中,我们仅考察了连续介质运动的运动学描述,而没有考虑到引起运动和变形的因素。本章我们将引入应力等概念,并给出连续介质力学的基本定律:质量守恒定律、动量平衡定律、动量矩平衡定律、能量守恒定律及熵不等式。
2.1 应力矢量与应力张量
在物体的运动中,物体的两部分之间或物体与其外界间的力学作用是通过力来描述的。在连续介质力学中我们主要研究三种类型的力:(1)一个物体的两部分之间的接触力;(2)由外界作用于物体边界上的接触力;(3)由外界作用于物体内部点的非接触力(如重力、离心力等) 。在另一方面,由于(1)(2)型的力总是通过某一接触面发生作用的,因此通常把作用于单位接触面积上的接触力称为表面力,或简称面力;由于(3)型力作用于物体整个体积内所含的物质点,因此通常把它称为体积力,或简称体力。
在连续介质力学中重要的公理之一就是关于接触力形式的柯西假设。柯西假设在运动过程中的时刻t 对于任何物质坐标X 和与之对应的接触面S 上的单位法矢量n ,表面力的存在形式为
t =t (X , t , n ) (2.101) 通常,我们规定t =t (X , t , n )指向接触面S 的外法向时为正,反之为负(见图2.1). 现在不管在X 和S 面与S' 面的曲率相差多少。
为了研究物体内部的力学状态,我们把一物体用一假想平面S 截断成两部分A 和B ,如图2.3所示。此时S 面就是A 和B 相互作用的接触面,B 部分对A 部分一点的作用,便可以用A 部分截面上的表面力t n 来表征,我们称之为应力矢量。反过来,考虑A 部分对B 部分作用,按照牛顿的作用与反作用定律可得应力矢量t -n 。它与t n 作用于同一平面上的同一点处,并且大小相等,方向相反。即
t n =-t n (2.102) 对于物体内部的一点P ,通过它可以有无穷多个方向的截面,而对于不同方向的截面,应力矢量也就不同,这种复杂情况只有引进应力张量的概念才能充分地加以描述。为了刻画一点的应力状态,设想在一点P 的附近任意给定一个单位法矢量为
n =(cos α1, cos α2, cos α3, )
=(e 1⋅n , e 2⋅n , e 3⋅n ) (2.103) 的平截面。相应地,过P 点沿活动标架作三个坐标平面。于是它们在物体内截得一个微小四面体,如图2.4所示。在这个微小四面体的每一个面上,都受有物体的其余部分给它的作用力,不妨设在ABC 上受到的作用力为t ∆A ,在PBC ,PCA 与PAB 上的作用力分别为-t 1∆A 1、-t 2∆A 2与-t 3∆A 3,其中∆A 与∆A i 分别为各微小平面的面积,作用于微小四面体ABCP 上单位质量的体力为b 。
现在假设对物体的任何部分,特别是对微小四面体ABCP 而言,动量的变化率与作用的合力成正比。虽然这是个很自然且牛顿第二定律更强的新假设(因为牛顿第二定律只适用于整个物体) ,然而,它却不能用实验直接验证,因为不可
能做内部表面接触力的直接测定,这种力的存在与大小只能由其它量的观测推知。描述一点是应力张量,描述通过一点的某一截面是应力矢量。
对于微小四面体ABCP ,柯西定律给出 t ∆A -t 1∆A 1-t 2∆A 2-t 3∆A 3+ρb ∆V =t ∆A -t i ∆A i +ρb ∆V
1
=t ∆A -t i ∆A cos αi +ρbh ∆V
3
=tma =ρ∆Va
1
=ρh ∆Aa (2.104)
3
其中ρ为物体的密度,h 为P 点到ABC 面的距离,并且考虑到微小四面体的体积.
1∆V =h ∆A (2.105)
3
2.104式也可写成
11
t -t cos α+ρbh =ρha (2.106) i i
33
当微小四面体体积趋于零时,即∆A →0,∆h →0,则有
t =t i cos αi (2.107) 考虑到2.103式,并令
t i =T i 1e 1+T i 2e 2+T i 3e 3
=T ij e i (2.108) 则式2.107可写成
t =cos αi t i =(n ⋅e i )(T ij e j )
=n ⋅(T ij e i e j )=n ⋅T
=(T ij e i e j )⋅n =T T ⋅n (2.109)
=t i cos αi =(T ij e j )(e i ⋅n )
当T 对称时,则
t =n ⋅T =T ⋅n (2.110) 其中
T =T ij e i e j (2.111) 称为应力张量,其矩阵形式为
⎡T 11T 12T 13⎤
⎥ (2.112) T T T [T ]=⎢212223⎢⎥
⎢⎣T 31T 32T 33⎥⎦ 如果物体中一点处的应力张量已知,那么由式2.112可以得到通过该点的任
何截面上的应力矢量,因此应力张量完全地刻画了物体中一点的应力状态。 由A i 面上的应力矢量t i 的定义可知,t i =t i (X , t ),而由式2.108知
T ij =T ij (X , t ),因此式2.109变为
t (X , t , n )=n ⋅T (X , t ) (2.113) 上式就是柯西假设的具体形式,常称之为柯西基本定理。 下面我们研究应力张量T 的各分量的力学意义。考虑到
T ij =e i ⋅T ⋅e j =t i ⋅e j
故知,T ij 代表作用于e i 方向截面上的应力矢量t i 在e j 方向上的分量,如图2.5所示。
我们从图2.5看到,应力张量T 的对角线元素T ij (i =j )位于所作用平面的法线方向内,故称之为法向应力分量;应力张量T 的非对角线元素T ij (i ≠j )位于所作用的平面内,故称为剪切应力分量。
2.2 质量守恒定律
物质无论经过怎样形式运动,其总质量是不变的,这就是古典连续介质力学中的最重要规律之一—质量守恒定律。下面我们研究质量守恒定律的数学表达式。
设ρ为物体的密度,dV 表示物质点的体积,由于在运动过程中质量保持不变,所以
D
(ρdV )=0 (2.201) Dt
展开有
D
(dV )+D ρdV =0 (2.202) ρDt Dt
又由式
∂v D
(dV )=i dV =(divv )dV (2.203)
Dt ∂x i
于是式2.202可写成
D ρ∂v
+ρi =0 (2.204) Dt ∂x i
其不变性形式为
D ρ
+ρdivv =0 (2.205) Dt
其中
D ρ∂ρ∂ρ=+v (2.206) i Dt ∂t ∂x i
∂ρ
+v ⋅∇ρ =∂t
把上式代入式2.204,则得
∂ρ∂(ρv i ) +=0 (2.207) ∂t ∂x i
其不变性形式为 ∂ρ
+div (ρv )=0注明v 是张量,ρ只是一个函数,既不是矢量,又不是张量 ∂t
(2.208)
式2.205和式2.208就是质量守恒定律的数学表达式质量守恒方程,在连续介质力学中常称为连续性方程。
在正交曲线坐标系中,利用式:H i =
g i ⋅g j ,连续性方程可写为
∂ρ1
[∂1(ρv 1H 2H 3)+∂2(ρv 2H 1H 3)+∂3(ρv 3H 1H 2)]=0 (2.209) +
∂t H 1H 2H 3
在直角坐标系中,连续性方程为
∂ρ∂(ρv x )∂(ρv y )∂(ρv z ) +++=0 (2.210) ∂t ∂x ∂y ∂z
在柱面坐标系中,利用第第一部分二章式2.13.03,连续性方程为
∂ρ1∂(ρrv r )1∂(ρv θ)∂(ρv z )+++=0 (2.211)
∂t r ∂r r ∂θ∂z
在球面坐标系中,利用第一部分二章式式2.13.04,连续性方程为
∂ρ1∂ρr 2v r 1∂(ρsin θv θ)1∂(ρv ϕ) +++=0 (2.212) ∂t r 2∂r r sin θ∂θr sin θ∂ϕ
连续性方程也可用物质描述法表示。在这种情况下质量定恒定律要求
⎰V 0ρ(X , t 0)dV 0=⎰V ρ(x , t )dV (2.213)
()
其中V 是物质在现时刻所占据的体积,而V 0是物质在时刻t 0所占据的体积。于是 ⎰V 0ρ(X , t 0)dV 0=⎰V 0ρ[x (X , t ), t ]JdV 0
=⎰V 0ρ(X , t )JdV 0 (2.214) 因为这个关系式对任意体积V 0都必须成立,故得
ρ0=ρJ (2.215) 它表示ρJ 与时间无关,即
ρJ =const (2.216) 这就是物质形式的连续性方程。
2.3 动量平衡定律
欧拉把下列关系作为在连续介质中普遍成立的一般性原理:
Dm
=f (2.301) Dt
它称为欧拉第一运动定律。上式说明任意物体具有的动量的变化率等于作用于该物体上的合力f 。
设所研究物体在其体积V 上受有连续分布的体力和在其体积的边界面S 上连续分布的接触力f c ,因此物体上所受合力为
f =f b +f c (2.302) 其中
f b =⎰V ρbdV (2.303) f c =⎰S tdS (2.304) 物体的动量为
m =⎰V ρvdV (2.305) =⎰V ρ
Dx
dV Dt
于是将式2.302和式2.305代入式2.301则
⎰V ρadV =⎰S tdS +⎰V ρbdV (2.306)
D 2x
其中a =表示x 点的加速度。由式2.109,可将上式改写为
Dt 2
⎰S n ⋅TdS +⎰V ρbdV =⎰V ρadV (2.307)
利用高斯公式
⎰S n ⋅TdS =⎰V ∇⋅TdV (2.308) 则得
⎰S ∇⋅TdV +⎰V ρbdV =⎰V ρadV (2.309) 即
⎰V (∇⋅T +ρb -ρa )dV =0 (2.310) 考虑到V 的任意性,则
∇⋅T +ρb -ρa =0 (2.311) 即
divT +ρb =ρa (2.312) 需要指出的是,这里的散度是对于空间坐标的。上式称为柯西第一运动定律。其指标形式为
T ji ; i +ρb i =ρa i (2.313) 展开得
∂T ∂T ∂T
11+21+31+ρb 1=ρa 1 (2.314)
∂x 1∂x 2∂x 3∂T ∂T ∂T
12+22+32+ρb 2=ρa 2 (2.315)
∂x 1∂x 2∂x 3∂T ∂T ∂T
13+23+33+ρb 3=ρa 3 (2.316)
∂x 1∂x 2∂x 3
特别地,在静止的情况下,物体的加速度为零,则式2.313化为
divT +ρb =0 (2.317) 在弹性力学中,上式称为平衡方程。
在柱面坐标系中,利用第一部分第二章2.13.4.d 可得上式化为
∂T rr 1∂T θr ∂T zr T rr -T θθ
++++ρb r =0 (2.318)
∂r r ∂θ∂z r ∂T 1∂T θθ∂T z θT r θ-T θr
+++ρb θ=0 (2.319) r θ+
∂r r ∂θ∂z r ∂T 1∂T θz ∂T zz T rz
+++ρb z =0 (2.320) rz +
∂r r ∂θ∂z r
在球面坐标系中,利用第一部分第二章2.13.4.e ,则2.317式可化为
∂T rr 1∂T θr 1∂T ϕr 1 +++(2T rr +cot θT θr -T θθ-T ϕϕ)+ρb r =0 (2.321)
∂r r ∂θr sin θ∂ϕr ∂T 1∂T θθ1∂T θθ1 r θ+++2T r θ+T θr +cot (T θθ-T ϕϕ)+ρb θ=0 (2.322)
∂r r ∂θr sin θ∂ϕr
[]
∂T r ϕ1∂T θϕ1∂T ϕϕ1
+++[2T r ϕ+T ϕr +cot θ(T θϕ-T ϕθ)]+ρb ϕ=0(2.323) ∂r r ∂θr sin θ∂ϕr
2.4 动量矩平衡定律
对于任意物体下列关系式成立:
DM x 0
=l x 0 (2.401) Dt
其中M x 0表示物体绕x 0点的动量矩,l x 0表示作用于物体上的力对x 0点的合力矩。上式称为欧拉第二运动定律。
设作用于物体上的力矩只是由体力和接触力引起的, 故其合力矩为
l x 0=⎰V ρ(x -x 0)⨯bdV +⎰S (x -x 0)⨯td S (2.402) 而物体的动量矩为 M x 0=⎰V ρ(x -x 0)⨯
Dx
dV (2.403) Dt
将式2.402和式2.403代入式2.401,并考虑到
D Dx
ρx -x ⨯dV (2.404) ()V 0⎰Dt Dt D (x -x 0)Dx D 2x Dx D
⨯dV +⎰V ρ(x -x 0)⨯2dV +⎰V (x -x 0)⨯ =⎰V ρ(ρdV ) Dt Dt Dt Dt Dt
Dx Dx D 2x Dx D =⎰V ρ⨯dV 张量本身叉乘是0+⎰V ρ(x -x 0)⨯2dV +⎰V (x -x 0)⨯(ρdV )质量守恒0
Dt Dt Dt D t Dt
D 2x
=⎰V ρ(x -x 0)⨯2dV (2.405)
Dt
可得
⎰V ρ(x -x 0)⨯adV =⎰V ρ(x -x 0)⨯bdV +⎰S (x -x 0)⨯td S (2.406)
D 2x
其中a =表示x 点的加速度。考虑到式2.110和高斯公式,则
Dt 2
⎰V ρ(x -x 0)⨯bdV +⎰S (x -x 0)⨯td S -⎰V ρ(x -x 0)⨯ad V 可知⎰S (x -x 0)⨯td S =⎰V (x -x 0)⨯(n T )d S
=⎰V ρ(x -x 0)⨯bdV +⎰S (x -x 0)⨯(n ⋅T )dS -⎰V ρ(x -x 0)⨯adV =⎰V ρ(x -x 0)⨯(b -a )dV +⎰S n ⋅T ⨯(x -x 0)dS 混合积互换 =⎰V ρ(x -x 0)⨯(b -a )+∇⋅⎡⎣T ⨯(x -x 0)⎤⎦dV 积分定理
=⎰V ρεljk (x l -x 0l )⨯(b j -a j )e k +∂i ⋅⎡⎣T ij ⨯(x l -x 0l )εjlk e k ⎤⎦dV 张量运算 =⎰V εljk e k {ρ(x l -x 0l )⨯(b j -a j )+∂i [T ij (x l -x 0l )]}dV
=⎰V εljk e k [ρ(x l -x 0l )⨯(b j -a j )+T ij ; i (x l -x 0l )+T ij ∂i (x l -x 0l )]dV
{}
{}
⎤=⎰V εljk e k ⎡⎣(x l -x 0l )⨯(T ij ; i +ρb j -ρa j )+T i j δil ⎦dV 根据平衡方程,红色部分为0
=⎰V εljk T ij δil e k dV
=
⎰
V
εijk T ij e k dV
=0 (2.407) 考虑到体积V 的任意性,得
εijk T ij =0 (2.408) 因此,T ij 必须对称张量,即
T ij =T ji (2.409) 或
T
T =T (2.410) 上式叫做柯西第二运动定律。柯西第二运动定律限定应力张量为对称张量,其中只有六个独立分量。
2.5 能量守恒定律
在连续介质中,如果只研究力学量的影响,而不考虑热学效应,那么连续介质的能量守恒定律可以直接由运动方程导出。首先,将运动方程
Dv
∇⋅T +ρb =ρ (2.501)
Dt
点乘速度矢量v
Dv
v ⋅(∇⋅T )+ρv ⋅b =ρv ⋅ (2.502)
Dt
在体积V 上积分
Dv
=V v ⋅(∇⋅T )dV +⎰V ρv ⋅bdV (2.503) ⎰V ρv ⋅
Dt ⎰
考虑到
Dv D ⎛1⎫
⎰V ρv ⋅=⎰V ρ v ⋅v ⎪dV
Dt Dt ⎝2⎭
D ⎛ρ1D ⎫
=⎰V v ⋅vdV -v ⋅v (ρdV )质量守恒0 ⎪⎰V
D t ⎝22Dt ⎭D 1
ρv ⋅vdV =V
Dt ⎰2D 12
ρv dV =V ⎰Dt 2
DK
= (2.504)
Dt
1
上式表示在体积V 中的总动能K =⎰V ρv 2dV 的时间变化率。另外,考虑到
2
v ⋅(∇⋅T )=v j T ij ; i =v j T ij
()
, i
-v j ; i T ij
=∇⋅(T ⋅v )-(∇v ):T
=∇⋅(T ⋅v )-(D +W ):T
=∇⋅(T ⋅v )-D :T -W :T 反对陈与对称双点乘是0
=∇⋅(T ⋅v )-D :T (2.505) 这里利用了反称张量W 与对称张量T 之间的双重点积为零的性质。 把式2.504和式2.505代回到式2.503中去,则得
DK
+V D :TdV =⎰V ∇⋅(T ⋅v )dV +⎰V ρv ⋅bdV (2.506)
Dt ⎰
运用高斯公式把上式右边第一体积分化为面积分,并利用柯西假设t =t =n ⋅T ,则
⎰V ∇⋅(T ⋅v )dV =⎰S n ⋅(T ⋅v )dS 添加取掉()无影响
=⎰S t ⋅vdS (2.507)
将上式代入式2.506,于是我们得到在纯力学作用下的能量方程
DK
+⎰V D :TdV =⎰S t ⋅vdS +⎰V ρb ⋅vdV 其中D 是速度梯度的对称部分 Dt (2.508)
其中方程左边两项分别表示连续介质的动能和内能(应力生热) 的时间变化率,右边两项分别表示接触力和体力所做的功率。若令U 表示内能,则能量方程5.508也可简洁地写成
DK DU += (2.509) Dt Dt Dt 其中表示接触力和体力的功率,记号表示这个量不一定能写成某个函数
Dt
的全微分形式。
如果同时考虑机械能和非机械能,那么就必须用能量守恒定律的一般形式。能量守恒定律的一般形式可以表述为:动能加上内能对时间的变化率等于总功率加上在单位时间内供给物体的各种其它形式的能量。这些能量包括热能、化学能、电磁能等等。本书只考虑机械能和热能,于是能量守恒定律就化为著名的热力学第一定律的形式。
对于热力连续介质(thermomechanical continua)来说,通常把内能的时间变化率写成
DU D
= V ρudV ⎰Dt Dt
Du D dV +⎰V u =⎰V ρ(ρdV )是0 Dt Dt Du
dV (2.510) =⎰V ρDt
其中u 称为比内能,表示每单位质量的内能密度。另外,我们定义矢量f 为在单位时间内每单位面积的热通量,函数q 为在单位时间内每单位质量的热辐射量,于是物体总热量的增量变化率为
Q
=-⎰S f ⋅ndS +⎰V ρqdV (2.511) Dt
其中n 为物体表面的外法向,热通量矢量f 由傅立叶定律给出,即
f =k ∇T (2.512)
这里k 为热传导系数,T 为温度。
于是热力连续介质的能量方程可以写成
DK DU +=+ (2.513) Dt Dt Dt Dt 或写成积分形式
D 1Du
ρv ⋅vdV +ρ V ⎰V Dt dV =⎰S t ⋅vdS +⎰V ρv ⋅bdV -⎰S f ⋅ndS +⎰V ρqdV Dt ⎰2
(2.514) 把上式右边面积分化为体积分后再移到左端,则有
⎡1D (v ⋅v )Du ⎤
⎰V ⎢ρ+ρ∇⋅(T ⋅v )+ρv ⋅b -∇⋅f +ρq ⎤dV 高斯公式 ⎥dV =⎰V ⎡⎣⎦Dt Dt ⎦⎣2(2.515)
由于体积V 是任意的,故有
D ⎛v ⋅v 1⎫1 +u ⎪=∇⋅(T ⋅v )+v ⋅b -∇⋅f +q (2.516)
Dt ⎝2ρ⎭ρ
利用式2.505,则上式化为
Dv Du 11
v ⋅+=[D :T +v ⋅(∇⋅T )]+v ⋅b -∇⋅f +q (2.517)
Dt Dt ρρ
整理得
Du 111⎛Dv ⎫ =D :T -∇⋅f +q +v ⋅ ∇⋅T +ρb -ρ⎪平衡方程0
Dt ρρρ⎝D t ⎭(2.518)
考虑到运动方程成立,则有
Du 11
=D :T -∇⋅f +q (2.519)
Dt ρρ或
Du 11∂f i
=D ij T ij -+q (2.520)
Dt ρρ∂x i
上式表示物体内能的时间变化率等于应力功率和吸收的热量之和。 式2.513、式2.514、和式2.519都是能量守恒定律的表现形式。
2.6 状态方程熵定律
完整地表征一个热力学统称做是对这个系统状态的描述。用来描述这个状态的物理量称状态参数。状态参数随着时间变化表征一个热力学过程。但是,在一般情况下,这些状态参数并不全是独立的,它们之间存在着某种关系。这种关系就称为状态方程。如果某个状态参数可以通过其它几个状态参数表出,则称它为状态函数。
现在,我们考虑一个均匀的热力学系统,它处于平衡状态,即在没有外界影响的条件下,系统的各部分在长时间内不发生任何变化。描述这样一个热力学系统的状态参数为:几何参数V(体积) 、力学参数p(压力) 及热力学参数T(温度) 。联系这三个量的关系的状态方程可写成
F (p , V , T )=0 (2.601)
这里需要指出的是,对于一定的物质来说,状态方程是普遍适用的,也就是说,构成热力学系统的物质一经选定,状态方程的具体形式也就确定了。 例如对于完全气体而言,状态方程的具体形式可写成
m pV =R 0T (2.602)
M
其中m 为气体的质量,M 为分子量,R 0是克分子气体常数。
在上一节我们曾叙述过热力学第一定律,它公设机械能和热能可以互相转换,但是,只根据热力学第一定律还不能判定这种转换过程是否可逆。事实上,所有的真实过程都是不可逆的,但可逆过程却是一个非常有用的假设,因为在许多情况下,能量耗损是可以忽略不计的。可逆性判据由热力学第二定律给出。
热力学第二定律公设存在两个独立状态函数:绝对温度T 和熵S 。它们有如下性质:绝对温度T 为一正量,它仅仅是经验温度θ(即我们通常见到的温度) 的函数,熵S 和体积V 一样,是一个广延量,而温度是与熵相对应的强度量,正如压强是与体积相对应的强度量一样。一个物体的强度量代表物质的内在性质,与物体的质量大小无关,而一个物体的广延量则可分解为物体上各个子部分上的广延量之和。因此,一连续介质的总熵S 可写成下列形式:
S =⎰V ρsdV (2.603) 这里s 表示连续介质中的熵密度,即每单位质量中的熵。
一个系统的熵既可由于与外界相互作用而发生改变,也可由于系统内部发生变化而改变,因此 ds =ds (e )+ds (i ) (2.604) 这里ds 是熵密度的增量, ds (e )是由于与外部相互作用而引起的熵密度增量。ds (i )是由于系统内部发生变化而引起的熵密度的增量。ds (i )决不能为负值。它在可逆过程中为零,在不可逆过程中为正,即
ds (i )>0 (不可逆过程) (2.605) ds (i )=0 (可逆过程) (2.606) 在可逆过程中,如果令dq (R )表示供给系统的每单位质量的热量,则ds (e )可表示为 ds (e )=
dq (R )
T
按照热力学第二定律,在连续介质所占据的物理空间中总熵的时间变率不小于通过连续介质表面流入的熵与连续体内部源产生的熵之和。在数学上,这个熵原理可以以积分形式表示为
f ⋅n d
(2.608) ⎰V ρsdV ≥⎰V ρedV -⎰S
dt T
称之为克劳修斯—杜姆不等式,其中e 为单位质量中的局部熵源。上式中的等号成立时表示可逆过程,不等号成立时代表不可逆过程。 利用质量守恒定律
d ds d
ρdV +⎰V S (ρdV ) ⎰V ρsdV =⎰V
dt dt dt
ds
ρdV =⎰V
dt
(可逆过程) (2.607)
和高斯公式 f ⋅n ⎛f ⎫ ⎰V dS =⎰V ∇⋅ ⎪dV T ⎝T ⎭
考虑到体积V 的任意性,则由式2.608可得克劳修斯—杜姆不等式的微分形式 ds 1⎛f ⎫ -e -∇⋅ ⎪≥0 (2.609) dt ρ⎝T ⎭
2.7 主应力最大剪应力
t =n ⋅T 表示物体中一点周围不同方向上的应力矢量公式,当应力张量已知时,在给定的任何一个方向n 上的应力矢量就由t =n ⋅T 给出。下面,我们将要讨论的问题是,对于某给定点来说,在什么方向上法向应力T n 取驻值。这个问题归结为在n 为单位矢量的条件下,即 2222 n =n ⋅n =n 1+n 2+n 3
=n k n k =1 (2.701) 时,求T n 的条件极值问题。运用大家所熟知的拉格朗日乘子法,有 ∂T n ∂f -λ=0 (2.702) ∂n i ∂n i
其中f 为约束条件
f (n )=n ⋅n -1=-1n k n k =0 (2.703) 考虑到T ij =T ji ,则由式2.110可得
T n =n ⋅t =n ⋅T ⋅n
=(n k e k )⋅(T pq e p e q )⋅(n l e l )
=δkp n k T pq n l δql
=n p T pq n q (2.704) 将上式代入式2.702,则 ∂T n ∂f -λ ∂n i ∂n i
∂(n p T pq n q )-λ∂(n k n k -1) ∂n i ∂n i
∂n p ∂n q ∂n k =T n +n T -2λn k pq q p pq ∂n i ∂n i ∂n i
=δpi T pq n q +n p T pq δqi -2λδki n k
=2(T iq n q -λn i )=0 (2.705) =-或写成不变性形式,即
T ⋅n =λn (2.706) 或
(T -λI )⋅n =0 (2.707) 写成展开形式,则为
(T 11-λ)n 1+T 12n 2+T 13n 3=0
T 21n 1+(T 22-λ)n 2+T 23n 3=0
T 31n 1+T 32n 2+(T 33-λ)n 3=0 (2.708) 上列方程中n 具有非零解的充分必要条件是它的系数行列式为零,即 T ij -λδij =0 (2.709) 或
32 λ-I 1λ+I 2λ-I 3=0 (2.710)
其中
I 1=T 11+T 22+T 33=T ii =trT (2.711) T 11T 12T 11T 13T 22T 23++ I 2= T 21T 22T 31T 33T 32T 33
1 =(T ii T jj -T ij ) 2
12 =(trT )-trT 2 (2.712) 2
T 11T 12T 13[]
I 3=T 12T 22T 23=T ij =det T (2.713)
T 31T 32T 33
这里I 1,I 2,I 3是应力张量T 的三个主不变量,分别称为第一、第二、第三应力
123不变量。方程的解λ1,λ2,λ3为特征值,n ,n ,n 为特征矢量。其中若
λi ≠λj ,则n i ⊥n j 。
i i i 事实上,在n 方向上法向应力值就是n 所对应的特征值。将式2.706与n 点
乘,得 i i i i i λi =λi n ⋅n =n ⋅T ⋅n =t i ⋅n (2.714)
i i 则λi 就是n 方向上的应力,称为主应力,而n 称为主方向,主方向所确定的平面
称为主平面。
i i j j 若n 和n 不两个不同的主方向(i ≠j ),则在n 面上n 方向的剪应力T ij 为 i j i j T ij =n ⋅T ⋅n =λi n n =0 (2.715)
123故主应力平面上的剪应力为零。若以(n ,n ,n ) 为坐标单位基矢量,并令
λi =T i ,则应力张量矩阵具有下列形式:
⎡T 100⎤⎥ (2.716) 0T 0 [T ]=⎢
2⎢⎥⎢⎣00T 3⎥⎦
即 i i T =Tn i n (2.717)
现在我们来讨论最大剪切应力问题。为了计算方便,不妨将坐标系选取在主方向上,即取(e 1,e 2,e 3) 为主方向。设n 是通过物体内一点的某一平面的单位法向矢量,则
n =n 1e 1+n 2e 2+n 3e 3=n k e k (2.718) 作用于该平面的应力矢量分量为 t =n ⋅T =n k e k ⋅Te i e i (2.719) n k T δki e i =n i Te i (2.720) 在该平面上的法向应力为
T n =t ⋅n =n i ⋅Te i i ⋅n j e j =n i n j T i δij =n i 2T i (2.721) 若以T S 表示该平面的总剪应力的大小(如图2.6) ,则 22 T S =t -T n (2.722)
即
22 T S =t ⋅t -T n 2 =n i T i e i ⋅n j T j e j -(n i n i T i ) 2
=n i n i T i 2-n i n i T i (2.723) =n i 2T i 2-T n 2
2 我们仍运用拉格朗日乘子法计算T S 的驻值,考虑到n 为单位矢量,令
f (n )=n ⋅n -1=n i n i -1=0 (2.724) ()2则 ∂T S 2∂f +λ=0 (2.725) ∂n k ∂n k
其中 ∂n ∂T S 2∂T n 2=2n T -2T i i n ∂n k ∂n k ∂n k
⎛∂n i ⎫2 =2nT δ-2T 2nT ⎪ i ik n i ∂n k ⎭⎝
=2n i T i 2δik -4T n (n i T i δik ) =2n k T i 2-4n k T n T i =2n k (T i 2-2TT i n ) (2.726) ∂f ∂(n i n i -1) =∂n k ∂n k
∂n i =2n =2n i δik =2n k (2.727) i ∂n k
于是 2n k (T i 2-2TT i n +λ)=0 (2.729) 即
2n 1T 12-2T n T 1+λ=0 (
2n (T
2n (T 2
3-2T n T 3
利用条件2.724,则方程组2.730显然有一组解
n 1=±1,n 2=n 3=0
n 2=±1,n 1=n 3=0
n 3=±1,n 1=n 2=0 (2.731)
2但是这组解所确定的平面就是主平面,而在主平面上T S =0,这不是我们所要
求的解。 2223-2T n T 2)+λ)=0 +λ)=0 (2.730)
假定在式2.730中n 1≠0,n 2≠0,n 3=0,则
2 T 1-2T n T 1+λ=0
2 T 2-2T n T 2+λ=0 (2.732) 将上列两式相减,则有
T 12-T 22-2T n (T 1-T 2)=0 (2.733) 故得 T 1-T 2 (2.734) 222把它代入式2.721中并与n 1+n 2=1联立,则可解得 11 n 1=±,n 2=±,n 3=0 (2.735) 2
这时n 方向与主方向e 2,e 3成45度角。
同样,若设n 1≠0,n 2=0,n 3≠0和n 1=0,n 2≠0,n 3≠0则对应的n 值分别为 11 n 1=±,n 2=0,n 3=± (2.736) 22
和
11n =0n =±n =± 1,2,3 (2.737) 22
考虑到上列三组驻值,则
1(e 1+e 2)时,T S =±1(T 1-T 2) (2.738) 当n =±22
1(e 3+e 1)时,T S =±1(T 3-T 1) (2.739) 当n =±22
1(e 2+e 1)时,T S =±1(T 2-T 3) (2.740) 当n =±22
因此,剪切应力的最大值由下列三个值中的最大值给出 T 1-T 2T 3-T 1T 2-T 3 ,, (2.741) 222
或 (T )-(T n )min (T S )max =n max (2.742) 2
T n =