量纲分析法
第三节 量纲分析法
量纲分析是20世纪初提出的, 在物理领域中建立数学模型的一种方法,它是在经验和实验的基础上, 利用物理定律的量纲齐次原则,确定各物理量之间的关系。
3.1 量纲齐次原则与Pi 定理
许多物理量是有量纲的,有些物理量的量纲是基本的,另一些物理量的量纲则可以由基本量纲根据其定义或某些物理定律推导出来。例如在动力学中,把长度l , 质量m 和时间t 的量纲作为基本量纲,记为
[l ]=L , [m ]=M , [t ]=T ; 而速度v ,力f 的量纲可表示为[v ]=LT -1, [f ]=MLT -2.
在国际单位制中,有7个基本量:长度、质量、时间、电流、温度、光强度和物质的量,它们的量纲分别为L 、M 、T 、I 、Θ、J 、和N ;称为基本量纲。任一个物理量q 的量纲都可以表成基本量纲的幂次之积,
[q ]=L αM βT γI δΘεN ξJ η
量纲齐次性原则:用数学公式表示一个物理定律时,等式两端必须保持量纲一致。 量纲分析就是在保证量纲一致的原则下,分析和探求物理量之间关系;先看一个具体的例子,再给出量纲分析的一般方法。
例3—1: 单摆运动,质量为m 的小球系在长度为l 的线的一端,线的另一端固定,小球偏离平衡位置后,在重力
mg 作用下做往复摆动,忽略阻力,求摆动周期t 的表达式。
解:在这个问题中有关的物理量有t , m , l , g 设它们之间有关系式
α1α2α3
t =λm 1l g
---------------(3.1)
其中α, α2, α3为待定常数,入为无量纲的比例系数,取(3.1)式的量纲表达式有
[t ]=[m ]α[l ]α[g ]α
1
2
3
整理得:T =M
α1
L α2+α3T -2α3 --------------(3.2)
由量纲齐次原则应有
⎧α1=0
⎪
⎨α2+α3=0 ---------------(3.3)
⎪-2α=1
3⎩
解得:α1=0,
α2=
1
, 2
α3=-
1
, 代入(3.1)得 t =λ2
l -------(3.4) g
(3.4)式与单摆的周期公式是一致的
下面我们给出用于量纲分析建模的 Buckingham Pi定理,
定理:设n 个物理量x 1, x 2, , x n 之间存在一个函数关系
f (x 1, x 2, , x n )=0 --------------(3.5)
[x 1] [x m ]为基本量纲,m ≤n 。x 1的量纲可表示为
[x i ]=∏[x j ]
j =1
m
αij
i =1, 2, , n
矩阵A =(αij ) n ⨯m 称为量纲距阵,若RankA =r , 则(3.5)式与下式等价,
F (π1π2 πn -r ) =0
其中F 为一个未定的函数关系,πs 为无量纲量,且πs 可表示为
πs =∏x i
i =1
n
βi (s )
-----------(3.6)
T
(s ) (s )
而β(s ) =(β1(s ) β2 βn ) 为线性齐次方程组A
β=0的基本解向量.
利用Pi 定理建模,关键是确定与该问题相关的几个基本量纲的无量纲量
π1π2 πn -r ,
作为量纲分析法的应用,下面我们介绍航船阻力的建模.
3.2 航船的阻力,
长l 吃水深度h 的船以速度v 航行,若不考虑风的影响,航船受到的阻力f 除依赖于船的诸变量l 、h 、v 以外,还与水的参数——密度P ,粘性系数μ, 以及重力加速度g 有关。
我们利用pi 定理分析f 和上述物理之间的关系 1. 航船问题中涉及的物理量及其量纲为
⎧⎪⎪⎪⎪⎨⎪⎪⎪⎪⎩
[f ]=LMT -2[l ]=L [h ]=L [v ]=LT -1[ρ]=L -3M [μ]=L -1MT -1[g ]=LT -2
我们要寻求的关系式为
ϕ(f , l , h , v , ρ, μ, g ) =0 ---------------(3.7)
这些物理量中涉及到的基本量纲为L 、M 、T 2. 写出量纲距阵
⎡1111-3-11⎤⎢⎥T
110 A =1000⎢⎥
⎢⎣-200-10-1-2⎥⎦
rank A =3
L
M T
3. 解齐次线性方程组A T β=0 , 可得 n -r a n k A =4 个基本解向量
⎧⎪⎪
⎨
⎪⎪⎩
β(1) =(0、1、-1、0、0、0、0) T
β(2) =(0、1、0、-2、0、0、1) T
(3) T
β=(0、1、0、1、1、-1、0)
T
β(4) =(1、-2、0、-2、-1、0、0)
由(3.6)式,可给出4个无量量纲
⎧⎪⎪
⎨
⎪⎪⎩
π1=lh -1
π2=lv -2g
-------------------(3.8)
π3=lv ρμ-1
π4=fl -2v -2ρ-1
由Pi 定理,(3.7)等价于下列方程
Φ(π1, π2, π3, π4) =0 -------------------(3.9)
这里Φ是未定的函数
由(3. 8), (3. 9)可得阻力f 的显式表达式,
f =l 2v 2ρψ(π1, π2, π3) ------------------(3.10)
其中ψ是由(3. 9)可得到的未知的函数关系, 在力学上,
v lg
称为Froude 数,记为
Fr ;
lv ρ
μ
称为Reynold 数,记为Re , 因此(3.10)又可写为
f =l 2v 2ρψ(, Fr , Re) ------------------(3.11)
4. 下面我们利用物理模拟进一步确定航船在水中的阻力。
设:f 、l 、h 、v 、ρ、μ、g 和f '、l '、h '、v '、ρ'、μ'、g '分别表示模型和原型中的各物理量, 由(3.11)有
f =l v
22
ρψ(
l v lv ρ, , ) h μlg
f '=l '2v '2ρ'ψ(
l 'v 'l 'v 'ρ', , ) h 'μ'l g l l '
=h h '
v lg =
v ''g '
lv ρ
=
l 'v 'ρ'
-------------------μ'
当无量纲量 (3.12)
μ
f '⎛l 'v '⎫ρ'
成立时 , 可得 --------------------= ⎪
f ⎝lv ⎭ρ
(3.13)
则此时由模型船的阻力f ,及其它的l 、v 、ρ、l '、v '、ρ';可确定原型船的阻力
2
f '.
下面, 我们讨论一下(3.12)成立的条件, 如果在实验中采用跟实际同样的水质, 则
p '=p , μ'=μ; 又g =g '
故可得 : (3.14)
要使得(3.14)成立 , 必有l =l ', h =h ';也即模型船与原型船一样大,这显然排除了物理模拟的可行性。若考虑选用不同的水质,s ⋅t 化为 (3.15)
由(3.15)可得
l l '
=h h 'v l =v 'l 'v l '
= ---------------------------v 'l
μ'≠μ, 仍设 ρ'=ρ 则(3.12)式
l l '=h h 'v l =v 'l '
3
2
v l 'μ
---------------------=⋅v 'l μ'
μ⎛l ⎫
= ⎪ , μ'⎝l '⎭
μ',显然无法找到如此小的粘性系数的液体 若按1:20的比例,μ=0⋅011
实际上的一种近似处理方法是,在一定条件下Re 数的影响很小,这样可近似得到, f ≈l v
2
2
ρψ(
l v , ) h lg
类似地分析,只要
l l '=h h '
3
v l
即有 =
''v l
f '⎛l '⎫
= ⎪ ----------------(3.16) f ⎝l ⎭
由(3.16)式就容易确定原型船的阻力f '
3.3 无量纲化 抛射问题
下面我们通过一个例子, 介绍如何使用无量纲化方法简化模型。
抛射问题:在某星球表面以初速度v 竖直向上发射火箭,记星球半径为r ,星球表面重力加速度为g ,忽略阻力,讨论发射高度x 随时间T 的变化规律.
模型建立:设x 轴竖直向上,t =0 时 x =0,火箭和星球质量分别记为m 1和m 2,由牛顿第二定律和万有引力定律可得:
=-k m 1 x
以x (0) =0,
m 1m 2
-----------------------(3.17)
(x +r ) 2
(0) =-g 代入(3.17)可得 km 2=gr 2 x
故得如下初值问题
⎧⎪⎪
⎨
⎪⎪⎩
r 2g
=-x
(x +y ) 2
x (0) =0 ----------------------(3.18) (0) =v x
(3.18)式的解可以表为 x =(t 、r 、v 、g ) 也即发射高度是以 r 、v 、g 为参数的t 的函数, 下面我们采用无量纲化方法化简方程(3.18),
显然抛射问题中的基本量纲为 L , T 而[x ]=L [t ]=T
;
[r ]=L [v ]=LT -1[g ]=LT -2
所谓无量纲化是指,对(3.18)式中的x 和t 分别构造且有相同的参数组合x c 和t c ,使得新变量
=
x x 0
=
t t 0
为无量纲量,其中 x c ,
t c 称为特征尺度或参考尺度;把方程(3.18)化为 对 的微分方程,即可简化模型,如何寻找特征尺度,这里我们以t c 为例,首先写出参数r 、v 、g 的量纲距阵A
1⎤⎡11 A =⎢⎥
0-1-2⎣⎦
t 的量纲向量为 (0, 1) T 记为 : β0
求解线性方程组 A β=
β0 得通解: β=(1, -1, 0) T +k (1, -2, 1) T
k =-1得
任取k ,即得到一种特征尺度,例如 k =0 得 t c =rv -1;
t c =vg -1; k =-
1-12-1
得t c =rg 同理可得x 的几种特征尺度r , v g 等 2
以下,我们利用不同的x c 和t c 化简(3.18) 1. 令 x c =r ,
t c =rv -1; 则=, =-1
d v 2d 2 =v =x 由 x ,
d r d 2
代入(3.18)可得:
2
1v =-ε, ε=
rg (+1) 2
(0) =0 ----------(0) =1⎧
⎪⎪
⎨
⎪⎪⎩
(3.19)
(3.19)式的解可表为 =(, ε) , 含一个独立参数且ε为无量纲量. 2. 令 x c =r , t c =
rg -1 , 类似地可将(3.18)化为 :
⎧⎪⎪⎨⎪⎪⎩
1
,
(+1) 2
(0) =0
2v (0) =, ε=rg =-2
-1
---------------(3.20)
3. 令 x 0=v g
, t 0=vg -1; 可将(3.19)化为
⎧⎪
⎪
⎨⎪⎪⎩
=-ε1v 2
, ε=2
rg (ε+1)
(0) =0(0) =1 ---------------(3.21)
按照现有科技能力, v
∴ε
=-1⎧在(3.21)中令ε=0,则有
⎪ ------------------------(3.22)
⎨(0) =0
⎪⎩(0) =1
2
+ , (3.22) 的解为: () =-2
代回原变量 得: x (t ) =-
12
gt +vt , ------------------------(3.23) 2
(3.23)式恰为假定火箭运动过程中所受星球引力 不变的运动方程。
小结:无量纲化是用数学工具研究物理问题时常用的方法,恰当地选择特征尺度不仅
可以减少参数的个数,而且可以帮助人们决定舍弃哪些次要因素
思考题: 在(3.19)和(3.20)式中, 令ε=0会得到什么结果?