量纲分析法建模.docVIP

§5 量 纲 分 析 法 建 模 量纲分析(Dimensional Analysis)是20世纪初提出的在物理领域中建立数学模型的一种方法，它在经验和实验的基础上利用物理定律的量纲齐次原则，确定各物理量之间的关系．本节在一个例子的引导下先介绍量纲齐次原则和著名的BuckinghamPi定理，然后用这个定理讨论一个力学问题的建模方法，并介绍量纲分析在物理模拟中的应用．最后给出一种简化模型的方法——无量纲化． 一、量纲齐次原则 许多物理量是有量纲的，有些物理量的量纲是基本的，另一些物理量的量纲则可 以由基本量纲根据其定义或某些物理定律推导出来．例如在研究动力学问题时常把长 度l、质量m和时间t的量纲作为基本量纲，记以相应的大写字母L，M和T．于 是速度v、加速度a的量纲可以按照其定义分别用LT 根据牛顿第二定律用质量和加速度量纲的乘积LMT 万有引力定律f=km1m2 r2-1和LT-2表示，力f的量纲则应-2表示．有些物理常数也有量纲，如fr2中的引力常数k，由 k=m1m2 -1可知其量纲应从力f、距离r和质量m的量纲求出，为LMT-2·L2·M-2=L3MT-2．通常，一个物理量q的量纲 记作[q]，于是上述各物理量的量纲为[l]=L，[m]=M，[t]=T，[v]=LT-1，[a]=LT-2，[f] =LMT-2， [k]= L3M-1T-2. 对于无量纲量α，我们记[α]=1(因为可视为[α]=L0M0T0)． 用数学公式表示一个物理定律时，等号两端必须保持量纲的一致，或称量纲齐次性(Dimensional Homogeneity)．量纲分析就是利用量纲齐次原则来寻求理量之间的关系 [6,20]．在叙述主要定理之前先看一个例子． 单摆运动 这是一个熟知的物理现象，质量为m的小球系在长度为l的线的一端，稍偏离平衡位置后小球在重力mg作用下(g为重力加速度)做往复摆动，忽略阻力．求摆动 第 36 页 周期t的表达式． 在这个问题中出现的物理量有t，m，l，g，设它们之间有关系式 其中α1，α2，α3是待定常数，λ是无量纲的比例系数．取(1)式的量纲表达式即 [t]=[m]α[l]α[g]α123将[t]=T，[m]=M，[l]=L，[g]=LT代入得 -2 按照量纲齐次原则应有 (3)的解为α1=0，α2=1/2，α3=-1/2，代人(1)式得t=λ (4)式与用力学规律得到的结果是一致的． 为了导出量纲分析建模的一般方法，将这个例子中各个变量之间的关系写作 lg （4） 进而假设(5)式形如 tmy1y2ly3gy4=π （6） 其中y1～y4是待定常数，π是无量纲常数．将t，m，l，g的量纲用基本量纲L，M，T 表示为[t]=LMT,[m]=LMT,[l]=LMT,[l]=LMT达式可写作(注意到[π ]=LMT) 00000101010010-2，则(6)的量纲表 即 Ly3+y4My2Ty1-2y4=LMT000 （7） 第 36 页 此方程组有一个基本解y=(y1,y2,y3,y4)T=(2,0,-1,1) （9） 代回(6)式得 t2l-1g=π （10） 而(5)式等价于F(π)= （11） (10)，(11)两式就是用量纲齐次原则从(5)式得到的结果．前面给出的(4)式只是它的特殊表达形式. 把从(5)式到(11)式的推导过程一般化，就是著名的Pi定理． 定理 设有m个物理量q1,q2, ,qm， 是与量纲单位的选取无关的物理定律*,X1,X2, ,Xn是基本量纲，n≤m． n q1,q2, ,qm的量纲可表为[qi]=∏ i=1Xiij,j=1,2,...,m (13) a 矩阵A={aij}n?m称量纲矩阵．若A的秩RankA=r (14) 设线性齐次方程组(y是m维向量) Ay=0 (15) T 的m-r个基本解为ys=(ys1,ys2, ,ysm),s=1,2, ,m-r (16) m 则πi=∏ j=1qjsj y 为m-r个相互独立的无量纲量．且 与(12)式等价．F表示一个未定的函数关系． [航船的阻力] 长l、吃水深度h的船以速度v航行，若不考虑风的影响，那么航船受的阻力f 除依 第 37 页 赖于船的诸变量l，h，v以外，还与水的参数——密度ρ、粘性系数μ，以及重力加速度g有关．下面用量纲分析方法确定阻力f和这些物理量之间的关系． 我们按照Pi定理中(12)～(18)式的步骤进行． 1．航船问题中涉及的物理量有：阻力f，船长l，吃水深度h，速度v，水的密度ρ，水的粘性系数μ，重力加速度g，要寻求的关系式记作 2 ．这是一个力学问题，基本量纲选为L，M，T．上述各物理量的量纲