0.引言
要想无量纲化一个方程,必须先确定方程涉及到的基本量纲。
Buckingham
π
\pi
π定理表明,
无
量
纲
参
数
个
数
=
有
参
数
个
数
−
基
本
量
纲
无量纲参数个数=有参数个数-基本量纲
无量纲参数个数=有参数个数−基本量纲
1.例子
下面用单摆问题作为例子来说明:
m
l
θ
¨
+
c
l
θ
˙
+
m
g
sin
θ
=
F
0
cos
ω
t
(1)
m l \ddot{\theta}+c l \dot{\theta}+m g \sin \theta=F_{0} \cos \omega t\tag{1}
mlθ¨+clθ˙+mgsinθ=F0cosωt(1)
基本量纲是时间(
t
t
t),长度(
l
l
l)和质量(
m
m
m)。1==角
θ
\theta
θ用弧度表示是无量纲的,因为弧度是两个长度的比值。==因此有量纲的量总共有6个
m
,
l
,
c
,
g
,
F
0
m, l, c, g, F_{0}
m,l,c,g,F0和
ω
\omega
ω。根据
π
\pi
π定理,应该有不超过
6
−
3
=
3
6-3=3
6−3=3个无量纲参数。
由于方程的每一项量纲必须相同,由式(1)的第一项可知,每一项的量纲都是
m
l
/
t
2
m l / t^{2}
ml/t2。由此可以反算
c
c
c的量纲是
m
/
t
m / t
m/t。
定义自然频率为
ω
0
=
g
/
l
\omega_{0}=\sqrt{g / l}
ω0=g/l,则方程可以化为
θ
¨
+
c
m
θ
˙
+
(
ω
0
2
)
sin
θ
=
F
0
m
l
cos
ω
t
(2)
\ddot{\theta}+\frac{c}{m} \dot{\theta}+\left(\omega_{0}^{2}\right) \sin \theta=\frac{F_{0}}{m l} \cos \omega t\tag{2}
θ¨+mcθ˙+(ω02)sinθ=mlF0cosωt(2)各参数的单位为
[
c
m
]
=
1
/
t
,
[
ω
0
]
=
1
/
t
,
[
F
0
m
l
]
=
1
/
t
2
,
[
ω
]
=
1
/
t
\left[\frac{c}{m}\right]=1 / t, \quad\left[\omega_{0}\right]=1 / t, \quad\left[\frac{F_{0}}{m l}\right]=1 / t^{2}, \quad[\omega]=1 / t
[mc]=1/t,[ω0]=1/t,[mlF0]=1/t2,[ω]=1/t
现在有4个有量纲参数和一个基本量纲
t
t
t。因此需要一个时间尺度,把时间
t
t
t也无量纲化。不失一般性,令时间尺度为
t
∗
t_{*}
t∗,则无量纲时间为
τ
=
t
/
t
∗
\tau=t / t_{*}
τ=t/t∗。
将以下公式
d
/
d
t
=
(
d
τ
/
d
t
)
d
/
d
τ
=
t
∗
−
1
d
/
d
τ
d / d t=(d \tau / d t) d / d \tau=t_{*}^{-1} d / d \tau
d/dt=(dτ/dt)d/dτ=t∗−1d/dτ
d 2 / d t 2 = t ∗ − 2 d 2 / d τ 2 d^{2} / d t^{2}=t_{*}^{-2} d^{2} / d \tau^{2} d2/dt2=t∗−2d2/dτ2代入式(2)中,重新整理,得 d 2 θ d τ 2 + c t ∗ m d θ d τ + ω 0 2 t ∗ 2 sin θ = F 0 t ∗ 2 m l cos ( ω t ∗ τ ) (3) \frac{d^{2} \theta}{d \tau^{2}}+\frac{c t_{*}}{m} \frac{d \theta}{d \tau}+\omega_{0}^{2} t_{*}^{2} \sin \theta=\frac{F_{0} t_{*}^{2}}{m l} \cos \left(\omega t_{*} \tau\right)\tag{3} dτ2d2θ+mct∗dτdθ+ω02t∗2sinθ=mlF0t∗2cos(ωt∗τ)(3)
在式(3)的基础上选择合适的 t ∗ t_{*} t∗可以将(3)继续化为 d 2 θ d τ 2 + α d θ d τ + sin θ = γ cos β τ \frac{d^{2} \theta}{d \tau^{2}}+\alpha \frac{d \theta}{d \tau}+\sin \theta=\gamma \cos \beta \tau dτ2d2θ+αdτdθ+sinθ=γcosβτ或 α d 2 θ d τ 2 + d θ d τ + sin θ = γ cos β τ \alpha \frac{d^{2} \theta}{d \tau^{2}}+\frac{d \theta}{d \tau}+\sin \theta=\gamma \cos \beta \tau αdτ2d2θ+dτdθ+sinθ=γcosβτ的无量纲形式。(式(3)和式(4)的参数定义不同)
牛不是基本单位,是导出单位。由 F = m a F=ma F=ma可知, N = 1 k g ⋅ m / s 2 \mathrm{N}=1 \mathrm{kg} \cdot \mathrm{m} / \mathrm{s}^{2} N=1kg⋅m/s2 ↩︎