无量纲化的方法

第10卷第3期　

2008年6月

安顺学院学报

JOURNALOFANSHUNUNIVERSITY

Vol110　No13Jun12008

理科教学研究

无量纲化的方法

刘锋　贾多杰　李晓礼　席国柱　吉永林

(西北师范大学物理与电子工程学院,甘肃　兰州730070)

摘　要:文章在详细叙述量纲Ⅱ定理后,的式子,并详细讨论“回代”方法;(其实就是做无量纲化处理)之后的方程形式,关键词:量纲;单位制;量纲π定理;中图分类号:O303　　-9507(2008)03-0078-03

　　引言:

,在作理论运算和数值计算时,往往必须做无量纲化(其实就是用所求的体系的几个主要的特征量作为相应物理量的单位)处理,这样做可以使理论运算简单,数值计算方便,物理方程转换为特定的数学方程时,便于数学处理。另外,对那些无法严格求解只能做近似计算的方程往往要做小参数展开,如参数展开摄动法[1],此时小参数才能有准确的“比较小”的意义。针对目前量子力学参考书提及的薛定谔方程的无量纲化变换[2]及自然单位法[3-4],文章做了比较。将物理方程转化为无量纲方程,通常的方法是令方程中的物理常数为1,得无量纲方程后进行求解。在所求结果中用一些“关系式”回代。这样所得结果与原方程的解是一致的,避免了计算过程量纲的出现,以上做法的依据是量纲Ⅱ定理。

(一)量纲

基本单位一旦确立,某种物理量的量度单位就由它们与基本量的关系式导出,通过基本量度单位表示的导出量度单位的表达式称为这个物理量的量纲式。量纲式可以用符号写成公式的形式,设x1,x2,…xm,是所选用单位制中的m个基本单位(在我所举的例子中以符号L表示长度单位,T表示时间单位,M表示质量单位),剧[P]代表导出量P的量纲式,如果有

[p]=[x1]1[x1]

a

a

2

…[xm]am

(1)

则指数(a1,a2,…am)称为物理量p的“量纲”。一个物理量的量纲与单位制的选择有关,所以在没有选择单位制之前谈量纲是没有意义的。

量纲可以看成是某个“矢量空间”的矢量。对(1)式两端取自然对数则有:

ln[p]=a1ln[x1]+a2ln[x2]+…+amln[xm]

(2)

由于各物理量以一定的关系式联系着,所以我们可以取其中的一些独立的物理量作为“基本量”并给它们规定一个“基本量度单位”,其它的物理量的量度单位将以确定的形式来导出。我们把基本量所采用的量度单位叫做基本量度单位,其它的物理量的单位称为导出单位。按照此种方法构成

[5]

的一套单位,构成一定的“单位制”。在不同的学科中我们

我们把ln[x1],ln[x2],…,ln[xm]看作是m维空间的“正交基矢”,则(a1,a2,…am)就是矢量ln[p]在各个基矢上的投影。则物理量p的“量纲”可以记作:

ln[p]→(a1,a2,…am)(二)量纲π定理

量纲π定理:假设某个物理问题中涉及n个物理量p1,

p2,…pn。他们量纲为:

ln[pi]→(a1i,a2i,…ami)(i=1,2,…,n)　　(4)

通常取不同的物理量作为基本量。例如:在力学问题中,在物理学研究中取长度,时间和质量作为基本量,而在工程技术上则以长度,时间和力作为基本量[6]而且基本量的个数也可以不一样。

收稿日期:2008-03-14

而我们所选单位制中有m个基本量(选取不唯一),则由此可以组成(n-m)个无量纲的量pm+1,pm+2,…,pn-m

,

作者简介:刘锋(1982106～),西北师范大学物理与电子工程学院研究生。研究方向:量子场论。

・78・

安顺学院学报　2008年第3期

物理理p1,p2,…,pn,存在的关系式:

f(p1p2,…,pn)=0　　(5)

0　-2　-2　0　1　-2　-1　01　-1　-0　0　0　-0　-0　00　-0　-1　1　0　-1　-1　1

假设正好是前m个物理量是我们取的基本量,则还有n

-m个物理量的量纲式可由那m个基本量表示出来.例如:

[pm+j]=[p1]

x

1j[

的秩是3。所以针对这个物理问题有3个线性独立的物理量

(但这个选取不唯一)现在根据π定理选取m,k,g作为基本

p2]

x

2j

…[pm]

x

mj

(j=1,2,…,n-m)　　(6)物理量,(他们的独立性可由他们量纲分量组成的行列式:

0　-2　-2

(7)

可以写作:

ln[pm+j]=x1jln[p1]+x2jln[p2]+…+xmjln[pm]

1　-1　-0≠0来证明)0　-0　-1

其写成分量的形式为:

a11　a12　…　a1ma21　a22　…　a2m

x1jx2j

=

a1,m+ja2,m+j

第二步:寻找剩余物理量对应的无量纲量。根据(7)式则有:

ln[l]=x1ln[m]+x2ln[k]+x3ln[g]

(13)

…　…　…　…

am1　am2　…　amm

2x1j

2x2j

…

xmj

…

am,m+j

根据(8)式则有:

0　-2　-21　-1　-00　-0　-x1x23

0(14)

则pm+j这个物理量可以得到一个无量纲量

pm+j=p1

′

p2

…pm2xmjpm+j　　(9)

=0　(10)

(5)式写成相应的无量纲的形式:f(m),pm+1,pm+2,…,pm+(n-m)

所以从(5)式量纲的式子化成无量纲式子(10)式,把我们选择的基本物理量换成1,的无量纲的量p即可,式子((三)举例

x2=-1

x3=1

所以可以得到一个无量纲的量:l′=m-x1k-x2g-x3l=

mg

(16)

,下面举一例来说明:例如我们要处理这样一个问题:将一质量为m的质点悬挂在劲度系数为k的弹簧下端,l为弹簧的原长,g为重力加速度,将质点从未伸长的弹簧个下端由静止释放,求弹簧长度x随时间6的变化(以弹簧上端点为坐标原点,竖直向下

为正方向)。牛顿运动方程为:m=2k(x-2

dt

mg　　(11)

2

同理可以得到另四个无量纲量:

x′=t′=

mg2m2

(16)(17)(18)(19)

l)+

初始条件为:

t=0时x=l和

=0　　(12)dt

221

=gdt′dtm2

2221

g2=dt′dt2

我们如何使(11)式化成无量纲的式子,当然初始条件也要变换。

第一步:选取本问题中的基本物理量。先写出本问题所涉及的所有物理量的量纲(我将采用国际单位制)

mT

k

g

l

t

22dt-2

dt-1

x

由(17)式可得

为时间量纲即为本问题的特征时间。k2

第三步:根据原来带量纲的函数关系,写出无量的关系

2

式。

把(11)式中在第一步中选取的基本物理量m,k,g换成

l。在第二步中找到无量纲的量p替换掉其相应的物理量p′,

010

-2-2001

100

则(11)式改写成无量纲的式子[对应于第(10)式]为:

001

2)+1-l′2=-(x′dt′

(20)

ML

10

01

01

01

初始条件(12)式为〔根据(16),(18)式〕

x′=

==l′mgmg

(21a)(21b)

我们可以把以上六个物理量看成是三维空间中的矢量,所以其基本量最多可以选三个(由线性代数中的最大线性无关向量知识可以得到)又因为矩阵:

2=　g21=0dt′dtm2

解(20)式满足初始条件(21)式可得:x′=l′+1-cost′

(22)

第四步:回代

・79・

安顺学院学报　2008年第3期

把(22)式根据(15),(16),(17)化成有量纲的结果:

2=+1-cos()mgmgm2即为:x=l+

(1-k

2cos()

m2

第一种方法:在自然单位制下即令:

(23)

h=u=ω=1

(29)

则无量纲的方程为:

22

)=2E′)+x′

dx

解出无量纲的方程的结果后用:-(30)

(24)

经检验(24)式是(11)式满足初始条们:(12)式的解。第五步:验证(通常可以省略)

为了看得更清楚我下面来比较(11)式和(20)式。把(20)式中带“′’的量按照(15),(16),(19)式换成不带“′”的量得到:

-1)+1=-(-2gmgmgdt

2

(31)x′=()2x

μω(32)E′=

ωh

(31)利(32)这两个式子回代即可回到原来单位制下的结果。

h

2

(25)

为特征长度

即为:

2m2=-k(x-l)+mgdt

(26)式利(11)式一样。

(26)

ω为特征能量h

第二种方法:(变量代换):令η()xxμω

h

:

(33a)(33b)(33b)

-22

)=2E′)(34)+x′

dx

(30)式和(34)式一样都是无量纲的式子。可见无量

综上所述:我们进行无量纲化,首先选取基本物理量,找其余物理量对应的无量纲的量(带“′”的),的物理量替换掉(所选基本物理量换成l,应的“′”换掉)即可,按照找寻带“′”的量所用[的(15),

(16),(17),(18),(19][类似于上例中的(22)式变成(23)式]。

(四)薛定谔方程,在自然单位中,怎么应用量纲Ⅱ定理

纲化其实就是重新选择合适的单位制。

参考文献:

金・摄动方法[M]1北京:高等教育山版社119931[1]顾德氵

[2]苏汝铿・量子力学[M]1上海:复旦大学出版社119971[3]曾谨言・量子力学(上册)[M]1北京:科学山版社119901

[4]曾谨言,钱伯初・量子力学专题分析(上)[M]1上海:高等教育出版社。

选用自然单位制其实就是用所求的体系的几个主要的特征量作为相应物理量的单位,在具体计算中,令相应的物理量为l。例如:一维线性谐振子的势能函数为:

ω2x22

薛定谔方程:V(x)=

2

2

(27)

19901

[5]赵凯华・定性与半定量物理学[M]1北京:高等教育山版社119941[6](苏)л1и1・谢多夫著・力学中的相似方法与量纲理论[M]1北京:科学

22

ωx

22dx

(28)

出版社,19821

[7]周世勋・量子力学教程[M]1北京:高等教育山版社119791

Methodofnormalization

LiuFengJiaDuoLiXiaoliXiGuozhuYonglin

(CollegeofPhysicsandElectrpnicEnginering,NorthwestNormalUniversity,Lanzhou730070,Gansu,China)

Abstract:Wepresentthedetailedcontentofthetheoremofdimensionless,Anexampleisgiventoshowhowtotyrnadimensionformulaintodimentionless1Inadditon,themethodhowtogetthedimensionformula,isstat2ed.FinallywemaketheSchrodingerEquationwithaconfiningpotentialdimensionlessby-usingthemethodwestatedpreviously1

Keywrods:Dimension;SI;Theoremofdimensionless;Dimentionless

(责任编辑:王德红)

・80・