← 返回知识卡片总目录

DS1号 · 自学教材系列 | 流体力学

第 8 章 相似原理与量纲分析

造一座大坝、一架飞机,不可能先做 1:1 的原型试错。工程师的办法是:做一个按比例缩小的小模型做实验,再把结果「换算」到原型上。凭什么小模型能预测大工程?答案是「相似原理」和「量纲分析」——它们用几个无量纲数,把不同尺度的流动统一起来。

🎯 本章学习目标

  1. 理解量纲与单位的区别,以及基本量纲与导出量纲。
  2. 理解量纲一致性原则,并会用它检查公式是否出错。
  3. 掌握π 定理(白金汉定理)的核心思想:把物理关系改写成无量纲数的关系。
  4. 认识几个最重要的相似准则数(Re、Fr、Eu 等)及其物理意义。
  5. 理解「模型与原型相似」需要满足什么条件,以及为什么难以完全满足。
  6. 会用白金汉 π 定理完整推导无量纲关系(如绕流阻力、泵相似三定律),并理解每个准则数的物理出身。

8.1量纲与单位:为什么要区分

先说清一对容易混淆的概念。

量纲是「种类」,单位是「尺子」。同一个量纲可以用不同单位量,但量纲本身不变。

基本量纲与导出量纲

力学里通常选三个基本量纲:

其他一切力学量都可以由这三者组合出来,称为导出量纲。例如:

常见物理量的量纲
物理量符号量纲(用 M, L, T 表示)
速度vL·T⁻¹
加速度aL·T⁻²
力FM·L·T⁻²
压强(应力)pM·L⁻¹·T⁻²
密度ρM·L⁻³
动力粘度μM·L⁻¹·T⁻¹

用方括号表示「取量纲」:比如 [v] = LT⁻¹、[p] = ML⁻¹T⁻²。记住这些基本写法,后面量纲分析全靠它。

🧠 一句话记忆 量纲是「种类」,单位是「尺子」。基本量纲 M、L、T,其余都是它们的组合。

量纲的「代数」:物理量的身份证

量纲不只是一个标签,它自带一套严格的运算规则,像代数一样可算:[x·y] = [x]·[y]、[x/y] = [x]/[y]、[xⁿ] = [x]ⁿ,而能相加的量必须量纲相同([x+y] = [x] = [y])。任何导出量的量纲都能唯一地「分解」成 M、L、T 的幂次组合——如同合数分解成质数,分解方式只有一个。速度 LT⁻¹、力 MLT⁻²、压强 ML⁻¹T⁻²、密度 ML⁻³……这些幂次组合就是每个物理量的「身份证号」。

为什么偏偏选 M、L、T 当基本量纲?这是人为「选坐标系」的问题,就像平面上可以任选方向当 x 轴。选 M、L、T,力就成了导出量([F] = MLT⁻²,来自 F = ma);老的工程单位制反过来选 F、L、T 当基本量纲,质量就变成导出量——两套坐标系完全等价,只是写法不同。国际单位制 SI 共 7 个基本单位,力学只用到 kg、m、s 三个,其余全是导出单位。一个重要推论由此而来:「基本」与「导出」是相对的,量纲体系的本质是那组幂次组合,而不是哪三个量被选中。

量纲运算直接产生单位换算。动力粘度 μ 的量纲是 ML⁻¹T⁻¹,SI 单位 Pa·s;由 1 Pa = 1 kg/(m·s²) 得 1 Pa·s = 1 kg/(m·s)。CGS 制把粘度单位叫「泊」(P):1 P = 1 g/(cm·s) = 0.1 Pa·s。再看运动粘度 ν = μ/ρ:[ν] = [μ]/[ρ] = (ML⁻¹T⁻¹)/(ML⁻³) = L²T⁻¹——量纲运算自动告诉你 ν 的 SI 单位是 m²/s,CGS 制叫「斯」(St),1 St = 1 cm²/s = 10⁻⁴ m²/s。整套粘度单位体系,全部从量纲代数里长出来。

基本量纲组合出导出量纲 M、L、T 三块积木按幂次组合,得到速度、力、压强、密度等导出量的量纲。 M 质量 L 长度 T 时间 速度 LT⁻¹ 力 MLT⁻² 压强 ML⁻¹T⁻² 密度 ML⁻³ 幂次组合唯一——量纲是物理量的「质因数分解」
图 8-1三个基本量纲像积木:不同幂次组合拼出全部导出量纲。正因为分解唯一,量纲才能当作物理量的「身份证」来检查方程。

📐 例题 1 从定义推量纲,再换算「厘泊 / 厘斯」

(a) 由牛顿内摩擦定律 τ = μ·(du/dy) 推导 [μ];(b) 推导运动粘度 ν = μ/ρ 的量纲与 SI 单位;(c) 20 ℃ 水的 μ = 1.0 mPa·s、ρ ≈ 1000 kg/m³,把 μ 与 ν 换算成 CGS 制的厘泊(cP)与厘斯(cSt)。

思路 量纲从定义式出发逐项运算;单位换算跟着量纲走,一步一验。

解 (a) [τ] = ML⁻¹T⁻²(压强类);[du/dy] = (LT⁻¹)/L = T⁻¹;故 [μ] = [τ]/[du/dy] = (ML⁻¹T⁻²)/(T⁻¹) = ML⁻¹T⁻¹。 (b) [ν] = [μ]/[ρ] = (ML⁻¹T⁻¹)/(ML⁻³) = L²T⁻¹,SI 单位 m²/s。 (c) 1 mPa·s = 10⁻³ Pa·s = 10⁻² P = 1 cP → μ = 1.0 厘泊; ν = μ/ρ = 1.0×10⁻³/1000 = 1.0×10⁻⁶ m²/s = 0.01 cm²/s = 1.0 厘斯。 交叉核对:1 cSt = 10⁻² St = 10⁻² cm²/s = 10⁻⁶ m²/s ✓。

答:[μ] = ML⁻¹T⁻¹、[ν] = L²T⁻¹;水的 μ = 1.0 cP、ν = 1.0 cSt。 ⚠ 易错点:1 Pa·s = 10 P(不是 1 P);「厘泊」配动力粘度、「厘斯」配运动粘度,两个「厘」不可互换——润滑油牌号(如 VG 46 指 40 ℃ 运动粘度约 46 cSt)用的是厘斯,拿厘泊去对牌号表必错。

工程警示:单位混用的 3.27 亿美元

1999 年,NASA 的「火星气候探测器」抵达火星时因轨道过低坠入大气层烧毁。事故调查结论令人瞠目:制造方提供的推进数据用英制单位(磅力),地面软件却按国际单位(牛顿)解读——1 lbf = 4.448 N,一个 4.45 倍的换算因子贯穿整段航程,轨道修正越修越偏,最终任务全损(探测器造价约 3.27 亿美元)。这次事故成了「量纲与单位必须严管」的最著名教材:现代工程规范因此普遍要求全流程使用统一单位制(SI),任何引入的实验数据、手册数据都要在交接处标注单位并复核量纲。审图与复核公式时先做一遍量纲检查,是最廉价也最有效的第一道防火墙。

辨析:无量纲量的量纲是「1」 雷诺数、马赫数、弧度、应变、阻力系数——这些量的量纲是 M⁰L⁰T⁰,即无量纲。无量纲量有三个宝贵性质:① 数值不随单位制改变(SI 与英制里同一个 Re 数值相同);② 可直接比较大小(Re = 2000 就是 2000,没有「米」「英尺」之争);③ 跨越具体尺度、连接模型与原型。量纲分析的终点,就是把一切有量纲量消化成无量纲数——这正是下一节的主题。

工程视角:图纸上的「单位陷阱」清单

量纲与单位的区分不是学究气——工程图纸与老设备铭牌上布满「历史遗留单位」,每年都在制造事故与返工。三对高频陷阱值得背下来:① 工程单位制残留——老压力表标 kgf/cm²(1 kgf = 9.81 N),1 kgf/cm² ≈ 0.098 MPa,与 bar(0.1 MPa)只差 2%,读错不至于出事、校准却必须较真;通风与过滤器行业至今爱用 mmH₂O(1 mmH₂O ≈ 9.81 Pa)。② 热量与功—— kcal 与 J 并存(1 kcal = 4187 J),制冷行业还用「冷吨」(1 冷吨 ≈ 3.52 kW);③ 动力粘度与运动粘度混标——油品样本几乎一律给运动粘度(cSt),传热计算却要动力粘度(Pa·s),中间隔一个密度。职业习惯只有一条:拿到任何数,先问「单位是什么」,再问「数值是多少」——顺序反了,量纲检查就形同虚设。

8.2量纲一致性原则

量纲分析的第一条铁律:量纲一致性原则。

量纲一致性原则 任何正确的物理方程,等号两边、以及相加的各项,量纲必须相同。你不能把「长度」和「速度」加起来,也不能让「力」等于「质量」。

这条原则有两个重要用途:

举个例子:伯努利方程的每一项 p、ρv²/2、ρgz,量纲都是 ML⁻¹T⁻²(压强),所以它们能相加,方程成立。如果哪天你写成了「p + ρv² = 常数」,量纲就不一致,一定是错的。

🧠 一句话记忆 方程两边量纲必须一致——「长度不能加速度」。量纲不对,公式必错。

为什么量纲必须一致:加法是「同类项合并」

「3 米 + 2 秒 = ?」——这个式子读不出任何结果,因为加法的本质是同类项的合并计数:相加的量必须度量同一件事,否则「和」没有定义。物理定律是对自然的陈述,它不应依赖你选米还是英尺、克还是千克;把单位换一下,方程两边必须等比例缩放,这只有在「两边量纲一致、每一项量纲一致」时才可能成立。量纲一致性因此是物理定律「单位制无关性」的代价:接受了普适性,就必须接受齐次性。

数学上这类方程叫量纲齐次方程。齐次性对运算有保护:求导与积分不破坏量纲一致性(d/dx 让量纲除以 [x],∫ 让量纲乘回 [x]),所以对物理方程放心地微分、积分都是安全的。但有一个硬约束:超越函数的参数必须是无量纲量。sin(x)、ex 的泰勒展开是 1 + x + x²/2 + …,第一项「1」是无量纲的,若 x 带量纲,x 与 1 相加立刻违反齐次性。于是看到 sin(ωt)、e−t/τ(τ 有时间量纲)就放心,看到 sin(t)、et(t 带时间量纲)就该报警——「超越函数只吃无量纲量」,是极实用的检查口诀。

📐 例题 2 量纲齐次性检查:三连问

(a) 同学抄公式时把泊肃叶流量写成 Q = πΔp·d³/(128μL),用量纲检查判断对错;(b) 验证孔口出流公式 v = √(2gh) 的量纲;(c) 有人说「量纲检查全过的公式一定对」,请举一个反例。

思路 把每个量的量纲代入候选公式,检查等号两边是否一致;第三问构造「量纲一致但数值错误」的公式。

解 (a) 左边 [Q] = L³T⁻¹。右边:[Δp·d³/(μL)] = (ML⁻¹T⁻²·L³)/(ML⁻¹T⁻¹·L) = (ML²T⁻²)/(M·T⁻¹) = L²T⁻¹ ≠ L³T⁻¹,量纲差一个长度——公式必错,正确形式应为 d⁴(分子变成 ML³T⁻²,除以 μL 后得 L³T⁻¹ ✓)。 (b) [2gh] = (LT⁻²)·L = L²T⁻²,开方得 LT⁻¹,与 [v] 一致 ✓(数字 2 与 √2 均无量纲)。 (c) 反例:把层流的 λ = 64/Re 抄成 λ = 32/Re——64 与 32 都是无量纲数,量纲检查完全通过,但数值是错的。量纲检查只能排除「结构性错误」,排除不了「系数级错误」。

答:(a) d³ 错,应为 d⁴;(b) 量纲正确;(c) 量纲对 ≠ 公式对。 ⚠ 易错点:(a) 中 μL 在分母上是「μ 乘 L」整体作除数,别把 L 漏乘;(b) 开方运算把量纲也开方(√(L²T⁻²) = LT⁻¹)——很多人开方时只开数字不开量纲;(c) 提醒:过不了量纲关的公式一定错,过了量纲关的公式还要过数值关。

工程视角:警惕「带量纲常数」的经验公式

审图复核的口头禅是「先查量纲、再查数值」,但有一类公式要格外当心:经验公式。明渠水力学广泛使用的曼宁公式 v = (1/n)·R2/3·√J(R 水力半径,J 水力坡度)简洁好用,但它的 1/n 带着量纲——按 SI 使用时 [n] = T·L⁻¹/3,换到英制手册同一个 n 数值就不通用,必须按手册规定整体换算。机械工程里常见的「功率 kW = T·n/9549」(T 单位 N·m、n 单位 r/min)同理:9549 = 60×1000/(2π) 这个常数把「r/min → rad/s」与「W → kW」两步单位换算打包在了一起——常数里浓缩了整套单位约定,输入单位一变常数就失效。判别办法:看公式里有没有「数值随单位制而定的常数」;有,就必须连同单位约定一起抄写,绝不能只抄数字。这类公式不属于量纲齐次方程,量纲一致性检查对它们只能「按约定单位验、不能跨单位验」。

辨析:量纲一致 ≠ 一定对,量纲不一致 = 一定错 把两个方向分开记:结构错误(量纲不一致)一票否决,量纲检查是高灵敏度的「排除法」;系数错误、符号错误、漏项错误都发生在量纲一致的前提之下,要靠代入数值、极端情形检验(如令 v→0、g→0 看结果合不合理)来抓。工程上常用的「代入特例检验」与量纲检查配合,就是公式自查的完整工具箱。

把量纲检查内化成三个快速动作,可以在十秒内给可疑公式做「体检」:① 单位回代——把每一项连同单位代进去算一遍,左右两边单位必须对得上;② 极限检验——让某个参数趋于零或无穷,公式行为是否符合直觉(例如阻力公式在 v→0 时必须归零);③ 特例对表——代入一个已知标准工况,与手册或经验值比对量级。三个动作任何一个报警,公式就该回去重推。这套动作对验算别人的公式、检查自己的推导、乃至审视 AI 生成或网上抄来的公式,同样有效。

8.3π 定理:把问题「压缩」成无量纲数

量纲分析最强大的工具是π 定理(也称白金汉 π 定理,因美国物理学家白金汉命名)。它的思想极其精妙,先看结论:

π 定理 如果一个物理问题涉及 n 个物理量,而这些量又只包含 r 个基本量纲,那么这个问题可以改写成由 (n − r) 个无量纲数(π 数) 组成的关系。

为什么这很强大?

它把问题的「自由度」大大压缩了。原来要研究 n 个变量的关系(比如 5 个变量),量纲分析告诉你:真正独立的只有 (n−r) 个无量纲组合(比如 2 个)。也就是说,决定现象的往往不是各个量本身,而是它们的几个无量纲组合。这既减少了实验次数,又揭示了现象背后的本质结构。

回想第 6 章的雷诺数——它就是 π 定理的一个现成例子。管流阻力本来取决于 ρ、v、d、μ 等一堆量,但量纲分析把它们压缩成「阻力系数 λ 只依赖于雷诺数 Re」这一个无量纲关系。莫迪图本质就是「λ = f(Re, ε/d)」这个无量纲关系的图像。你看,量纲分析早就悄悄在你身边了。

怎么用(思想层面)

  1. 列出影响问题的所有物理量(n 个)。
  2. 确定基本量纲个数(通常 r = 3:M、L、T)。
  3. 选出 r 个「重复变量」(要包含所有基本量纲),和其余变量逐个组合,消去量纲,得到 (n−r) 个无量纲 π 数。
  4. 把问题写成 π₁ = f(π₂, π₃, …) 的形式。

具体代数细节这里不展开,但核心思想要抓住:量纲分析不能给你完整的公式,但能告诉你「哪些量必须以什么组合出现」,从而把实验从「大海捞针」变成「有的放矢」。

🧠 一句话记忆 π 定理:n 个量、r 个基本量纲 → 只剩 (n−r) 个无量纲数决定一切。变量多不怕,无量纲组合才是关键。

把 π 定理用到底:白金汉方法的完整步骤

原版「思想层面」四步还不够过瘾,这里把可操作的完整流程走一遍,并说清每一步的「为什么」:

  1. 列全变量(最脆弱的一步):写出影响现象的全部有量纲量(n 个)。漏了关键变量,后面全盘皆错;混入无关变量,会多出无意义的 π 数。这一步靠物理洞察,不靠数学。
  2. 数出基本量纲(r 个,力学通常 3)。
  3. 选 r 个重复变量:要求它们合起来覆盖全部基本量纲,彼此又组合不出无量纲数(ρ、v、L 是黄金组合)。重复变量是搭积木的底座,其余变量轮流与它们组合。
  4. 逐个求解 π 数:π = (重复变量)a·(目标量),令总量纲为零,解出指数——本质是解一个 r 元线性方程组。
  5. 写出无量纲关系 π₁ = f(π₂, π₃, …),函数形式留给实验。

两个常被忽略的细节:① π 数的组合不唯一——π₁、π₂ 的幂次组合、乘积仍是合法 π 数(把 Re 换成 1/Re、Fr 换成 Fr² 都不算错),选哪一组看物理意义与行业习惯,不是数学规定的;② r 的严格说法是量纲矩阵的秩:若变量表里混进了本身就是无量纲的量(比如把 Re 也当变量列进去),有效量纲数会减少、π 数相应增多。工程实践按「r = 3」起步,再检查有无量纲相关的变量即可。

📐 例题 3 π 定理实战:推导绕流阻力的表达式

某物体(特征长度 L、迎流面积 A ∝ L²)以速度 v 在密度 ρ、动力粘度 μ 的流体中运动,受阻力 FD;重力与表面张力可忽略。用 π 定理推导 FD 的表达式,并说明它如何变成工程上熟悉的「阻力系数」形式。

思路 变量 FD、ρ、v、L、μ 共 5 个(n = 5),基本量纲 3 个(r = 3)→ 应得 2 个 π 数。选 ρ、v、L 为重复变量,逐个解指数。

解 ① π₁ = ρavbLcFD。量纲:Ma+1L−3a+b+c+1T−b−2 = 1。解:M: a+1=0;T: −b−2=0;L: −3a+b+c+1=0 → a = −1、b = −2、c = −2。 π₁ = FD/(ρv²L²)。 ② π₂ = ρavbLcμ。量纲:Ma+1L−3a+b+c−1T−b−1 = 1。解:a = −1、b = −1、c = −1。 π₂ = μ/(ρvL) = 1/Re。 ③ 无量纲关系:FD/(ρv²L²) = f(μ/(ρvL)),即 FD = ρv²L²·f(1/Re)。 ④ 工程习惯用迎流面积 A 替代 L²、把系数 ½ 吸收进「阻力系数」: FD = CD(Re)·(1/2)ρv²A。 这正是第 10 章绕流阻力要用的公式——量纲分析提前把它「猜」了出来,实验只剩测定 CD 随 Re 变化这一条曲线。

答:FD/(ρv²L²) = f(1/Re),工程形式 FD = CD(Re)·½ρv²A。 ⚠ 易错点:① 解指数时逐条写 M、L、T 三个方程再联立,别跳步——T 方程里 FD 贡献的是 T⁻²、μ 贡献的是 T⁻¹,容易混;② π₂ 取倒数写成正 Re 更合乎习惯,两者等价;③ 若问题发生在自由液面(兴波重要),变量必须加 g,于是多出 π₃ = v²/(gL) = Fr²——漏掉 g 就漏掉整个重力世界,模型实验将无法预测兴波阻力。

工程视角:无量纲化让实验「一劳永逸」

风洞、水池实验的数据为什么一律整理成无量纲曲线?因为只有无量纲数能跨越尺度。同一副翼型(不涉任何型号),在 30 m/s 的小风洞与 100 m/s 的大风洞里各测一组升力:直接画「升力 vs 风速」,两组点各奔东西;换成「升力系数 CL vs 攻角 α(按 Re 分层)」,两组点立刻落到同一族曲线上。模型与原型、实验与设计之间,靠无量纲数传递信息;手册里的翼型数据、管道阻力图、泵与风机的特性曲线,全部以无量纲形式出版,正是这个道理。反过来说,整理实验数据的第一课永远是:先把数据除以合适的尺度变成无量纲数,再去找规律。

辨析:量纲分析的能与不能 π 定理能告诉你「关系里该有几个无量纲数、它们由哪些量组合」,也能把五变量实验压缩成一条按 Re 分层的曲线;但它不能给出函数形式——f 是常数、是幂函数还是带突变,全靠实验或更深的理论;它也不能替你判断哪些变量该进入问题——变量清单是物理判断,清单错了,π 定理只是忠实地把错误的清单压缩成错误的框架。量纲分析的定位是「伟大的组织者」:它不生产物理,但让物理各就各位。

工程视角:泵与风机的「相似三定律」

π 定理最有威力的工业应用之一,是同一族几何相似的泵与风机之间的性能换算。以泵为对象列出变量:流量 Q(L³T⁻¹)、扬程能量头 gH(L²T⁻²)、轴功率 P(ML²T⁻³)、叶轮直径 D、转速 n(T⁻¹)、流体密度 ρ——共 6 个量、3 个基本量纲,π 定理给出 3 个无量纲数。取 ρ、n、D 为重复变量解出:

流量系数 πQ = Q/(nD³)  扬程系数 πH = gH/(n²D²)  功率系数 πP = P/(ρn³D⁵)

几何相似、工况相似的两台泵,这三个 π 必相等。把它变形,立即得到流体机械行业的「相似三定律」:

Q ∝ nD³  H ∝ n²D²  P ∝ ρn³D⁵

最常用的特例是同一台泵(D 不变)改变转速:流量随转速一次方、扬程随转速平方、功率随转速三次方。这正是第 7 章「变频泵比关阀门省能」的定量内核:转速降到 80%,功率只剩 0.8³ ≈ 51%;而若维持额定转速、用阀门把流量压到 80%,泵的功率几乎没降,被「省」掉的流量全部变成阀门的节流损失。

这组定律的工业价值在于「一次实验,全系列受益」:在模型泵上测出一组性能曲线,靠三条定律即可换算出同一系列(几何相似)任何叶轮直径、任何转速、甚至换一种密度流体(如输送烟气代替空气)的性能。泵与风机样本手册里「一张无量纲曲线覆盖整个系列」、水利装置用空气实验代替水实验,背后全是这几行量纲代数。旋转机械行业还把 πQ、πH 组合成一个「比转速」概念,用单一数字给泵「定性格」(离心式、混流式还是轴流式),这又是无量纲数的组织威力。

📐 例题 4 变频调速的账:转速降 20%,功率省一半

某离心泵在额定转速 n₁ = 1450 r/min 下运行:流量 Q₁ = 100 m³/h、扬程 H₁ = 32 m、轴功率 P₁ = 11 kW。因工艺需求流量降至 80%,改用变频调速(降转速)而非阀门节流。求新转速及对应的流量、扬程与轴功率;若年运行 8000 h,估算与「额定转速+阀门节流」相比的节能量(阀门节流时轴功率近似不变)。

思路 同一台泵,D 不变,用相似三定律的「变转速特例」:Q ∝ n、H ∝ n²、P ∝ n³。

解 ① 流量降到 80% → 转速等比:n₂ = 0.8×1450 = 1160 r/min。 ② H₂ = H₁×(0.8)² = 32×0.64 = 20.5 m(扬程自然降到系统所需水平,无需阀门「憋」出阻力)。 ③ P₂ = P₁×(0.8)³ = 11×0.512 = 5.6 kW。 ④ 节能:阀门节流时轴功率仍约 11 kW,变频后仅 5.6 kW,节省约 5.4 kW;按年运行 8000 h 计,年省电约 5.4×8000 = 4.3×10⁴ kW·h。

答:n₂ = 1160 r/min;Q₂ = 80 m³/h、H₂ = 20.5 m、P₂ ≈ 5.6 kW;年省电约 4.3 万 kW·h。 ⚠ 易错点:① 是转速的三次方不是平方——若误用平方(0.8²=0.64),P₂ 会算成 7.0 kW,节能账直接少算近三分之一;② 三定律的适用前提是「工况相似」(转速变化时流量、扬程沿同一条相似工况线移动),大幅偏离工况点时效率本身会变,需按泵特性曲线修正;③ 阀门节流「功率近似不变」是工程近似——实际会略降,但与 0.8³ 的降幅相比完全是两个量级。

8.4相似准则:几个关键的无量纲数

在流体力学里,有几个无量纲数反复出现,每个都代表「两种力量的比值」。它们是判断流动相似性的「身份证」。

流体力学重要相似准则数
名称符号与定义物理意义(力的比值)典型应用
雷诺数Re = ρvL/μ惯性力 / 粘性力管流、绕流的流态判断
弗劳德数Fr = v/√(gL)惯性力 / 重力明渠、船舶兴波、水坝溢流
欧拉数Eu = Δp/(ρv²)压力 / 惯性力压强场的相似
马赫数Ma = v/c流速 / 声速(惯性/压缩性)可压缩流、高速飞行
韦伯数We = ρv²L/σ惯性力 / 表面张力液滴、气泡、毛细流

这些数为什么叫「相似准则」?因为:两个流动如果几何相似,且对应的相似准则数相等,那么它们的流动就「相似」(速度场、压强场形状相同)。这就是模型实验的理论基础——只要小模型和大原型的关键准则数相等,小模型里发生的事,就能按比例「映射」到大原型上。

抓住本质 每个相似准则数都是「某种力」与「惯性力」(或另一种力)的比值。当某个准则数起主导作用时,说明那种力对流动的影响最重要。选择哪个准则数,取决于「哪个力主导这个流动」。

相似准则数的「出身」:从力的量纲比说起

表格里 Re、Fr、Eu 为什么恰好是「力的比值」?这不是巧合,而是把流体运动基本方程(N–S 方程)各项直接量纲化的结果。定常流动中,单位体积流体所受各种力的量纲分别是:惯性力 ~ ρv²/L、粘性力 ~ μv/L²、重力 ~ ρg、压力 ~ Δp/L。两两相除,比值自动消去尺度 L:

惯性/粘性 = (ρv²/L)/(μv/L²) = ρvL/μ = Re  惯性/重力 = (ρv²/L)/(ρg) = v²/(gL) = Fr²  压力/惯性 = (Δp/L)/(ρv²/L) = Δp/(ρv²) = Eu

看清楚了:每个准则数都是 N–S 方程里两项的相对权重。某个比值极大,说明那种力相对惯性可以忽略(或反之);比值在 1 附近,两种力势均力敌、谁也不能丢。所谓「判断哪个力主导」,本质是给方程做量级体检:逐项估算比值,把远小于 1 的项扔掉,留下的项决定流动的基本面貌。这一思想在流体力学里有个正式名字——尺度分析,它也是从 N–S 方程走到边界层方程(第 9 章)的那座桥。

准则数=两种力的比值 以惯性力为参照,分别与粘性力、重力、压力相比,得到雷诺数、弗劳德数平方、欧拉数。 惯性力 ~ ρv²/L 惯性力 ~ ρv²/L 惯性力 ~ ρv²/L Re ρvL/μ Fr² v²/(gL) Eu Δp/(ρv²) 粘性力 ~ μv/L² 重力 ~ ρg 压力 ~ Δp/L
图 8-2相似准则数的出身:以惯性力为「天平」的参照臂,与粘性力之比是 Re、与重力之比是 Fr²、与压力之比是 Eu——每个准则数都在称量 N–S 方程中两项的相对权重。

📐 例题 5 给三个流动「验身份证」:谁主导?

估算下列流动的主导相似准则(水的 ν 取 1.0×10⁻⁶ m²/s,海水取 1.2×10⁻⁶ m²/s):① 桥墩绕流——淡水,v = 2 m/s,墩宽 L = 10 m;② 船舶巡航——海水,v = 12 m/s,船长 L = 100 m;③ 高空巡航——10 km 高度声速 c ≈ 300 m/s,飞行速度为 0.82c。

思路 逐个计算准则数并比较量级:自由液面看 Fr,水下与绕流看 Re,高速可压缩看 Ma。

解 ① 桥墩:Re = vL/ν = 2×10/1.0×10⁻⁶ = 2×10⁷(强湍流,粘性只在近壁薄层起作用);Fr = v/√(gL) = 2/√(9.81×10) = 2/9.90 = 0.20。主导准则:Re(涡脱落、阻力由它控制)。 ② 船:Re = 12×100/1.2×10⁻⁶ = 1×10⁹;Fr = 12/√(9.81×100) = 12/31.3 = 0.38。两者都很大,但兴波阻力只认 Fr、且最难用理论修正——行业共识是模型实验保 Fr(详见 8.5)。 ③ 巡航:Ma = 0.82 > 0.3,压缩性不可忽略,主导准则为 Ma;此时 Re 仍有 10⁷ 量级,但流动的「性格」已由压缩性决定——第 11 章的主题。

答:① Re 主导;② 保 Fr 为主(Re 同样巨大,见 8.5 的取舍);③ Ma 主导。 ⚠ 易错点:「两个准则数都很大」不等于「都不重要」——准则数是两种力的比值,比值大说明被除的那个力可忽略、主导权在另一种力手里;② 也说明真实工程常是多力并存、抓主放次,这正是下一节模型实验困境的根源。

工程视角:风洞里的 Re 缺口

全尺寸风洞全世界屈指可数,大多数实验风洞的尺寸与风速决定了模型 Re 只有 10⁵~10⁶ 量级,比真实飞行器的 10⁷~10⁸ 差一个数量级以上。这个「Re 缺口」不只是数字差异:Re 直接控制边界层的转捩与分离位置(第 9 章),低 Re 模型上气流可能在本该湍流的地方保持层流,失速特性、最小阻力随之失真。工程师的三件应对:① 变流量补偿——加大风速把 Re 往上凑(受 Ma 限制,不能过头);② 增压、低温风洞——提高气体密度、降低粘度(低温氮气风洞能把 Re 提高 1~2 个数量级,代价高昂);③ 人工转捩——在模型表面贴粗糙带强制提前转捩,让边界层状态与原型对齐。三件手段共同的思想:让最关键的准则数对上,其余用修正与经验兜底。

把 Ma 与 We 也放上「天平」

同一套「力的比值」套路可以补全表格里剩下的准则数。弹性力(流体抵抗压缩的能力)~ ρc²(c 为声速,ρc² 正是体积模量),于是惯性/弹性 = ρv²/(ρc²) = (v/c)² = Ma²——马赫数平方称量的是惯性力与「流体被压缩的反抗力」的强弱。经验阈值 Ma > 0.3 的来历就在这里:小马赫数下密度变化 Δρ/ρ ≈ Ma²/2,Ma = 0.3 时约 4.5%,流体几乎「感觉不到」自己被压缩,按不可压缩处理误差可忽略;而巡航状态 Ma = 0.82 时这个量已达约三分之一(精确等熵关系给出的还要更高),气体的「性格」已彻底改变——那是第 11 章的世界。

表面张力(单位体积流体所受)~ σ/L²(σ 为表面张力系数),惯性/表面张力 = (ρv²/L)/(σ/L²) = ρv²L/σ = We。日常生活里 We 往往很小:一杯水表面平静光洁,是张力说了算;而高速小尺度流动里 We 急剧增大——雨滴以每秒几十米撞上机身、喷嘴把液体撕成雾滴、气泡在液体中受剪切变形破碎,全由 We 主导,超过临界值液滴就碎。选准则数的顺序建议:先问「有没有自由液面或界面」——没有界面,σ 根本不出场;有界面且尺度小、速度高,We 必查。

顺带补一个工程高频的「编外准则数」:空化数 σ = (p∞ − pv)/(½ρv²),即远场压力与饱和蒸气压之差同动压的比值——它本质是 Eu 的变体(把参考压差换成 p∞ − pv),决定液流中最低压点会不会「局部沸腾」出气泡(空化)。泵进口、螺旋桨叶背、泄洪洞的高流速区都怕它:σ 一旦跌破临界值,气泡溃灭的微射流会在金属表面啃出蜂窝状麻面(空蚀),几周就能报废一个叶轮。「按主导力挑准则数」的思想,在这里又一次出场。

8.5模型实验:相似条件的困境

理论很美,但真做模型实验时会撞上一堵墙:很难让所有相似准则数同时相等。

为什么难以同时满足?

举个具体例子:做一个船的缩比模型,在水池里做阻力实验。要让流动相似,需要弗劳德数 Fr 相等(因为兴波阻力主要由重力决定)和雷诺数 Re 相等(摩擦阻力由粘性决定)。但矛盾来了:

两个要求互相矛盾——用同一种流体(同 ν)根本不可能同时满足。这就是著名的「模型实验困境」。

📐 例题 6 把「困境」算出来:等 Re 要 66.7 m/s

某水下航行器原型长 Lp = 50 m、航速 vp = 8 m/s,在海水(νp = 1.2×10⁻⁶ m²/s)中工作。现做 K = 10 的缩比模型(Lm = 5 m),在淡水(νm = 1.0×10⁻⁶ m²/s)中拖曳实验。(a) 若要求 Re 相等,模型拖速多少?(b) 该拖速下模型的 Fr 是多少?与原型的 Fr 相比如何?

思路 由 Rem = Rep 解出 vm,再算两边的 Fr 对比,把「矛盾」变成看得见的数字。

解 (a) Re 相等:vmLm/νm = vpLp/νp,解得 vm = vp·(Lp/Lm)·(νm/νp) = 8×10×(1.0/1.2) = 66.7 m/s。 (b) 原型:Frp = 8/√(9.81×50) = 8/22.1 = 0.36; 模型:Frm = 66.7/√(9.81×5) = 66.7/7.00 = 9.5。 相差约 26 倍——模型被拖着以近乎失控的速度冲过水池,兴起的波纹与原型毫无相似可言。

答:等 Re 需拖速约 66.7 m/s,此时 Fr ≈ 9.5,是原型(0.36)的约 26 倍——常规水池里 Re 与 Fr 根本无法同时满足。 ⚠ 易错点:① K = 10 时 Lp/Lm = 10,代入公式别把分子分母颠倒;② 淡水粘度略低(ν 更小)只能「帮点小忙」,救不了数量级;③ 若反过来先满足 Fr 相等(vm = vp/√K = 2.53 m/s),则模型 Rem = 2.53×5/1.0×10⁻⁶ ≈ 1.3×10⁷,只有原型(3.3×10⁸)的约 1/26——两边都过不好,这就是困境的定量面貌。

工程师怎么办?

现实的做法是抓住主要矛盾、放弃次要:

工程智慧 模型实验几乎永远是「部分相似」——只保证主导力的准则数相等,其余靠经验修正。所以模型实验既是科学,也带着工程师的经验判断。

模型实验与原型换算 缩比船模型在水池中实验,保证关键相似准则数相等,即可把模型测得的阻力换算到原型。 缩比模型(小) 保证 Fr(或 Re)相等 原型(大) 按相似律换算阻力 →
图 8-3用缩比模型做实验,保证关键相似准则数相等,再把结果按相似律换算到原型。

破局的第三条路:自模区(自动模化区)

「抓主要矛盾」和「换介质」之外还有第三条路,它藏在阻力曲线的形状里。大量实验数据显示:许多无量纲系数在准则数足够大之后趋于常数——完全粗糙区的 λ 只看 ε/d(第 7 章)、钝体绕流的 CD 在 Re 超过 10³~10⁴ 后几乎不变(第 10 章)。既然系数不再随 Re 变化,那么只要模型与原型的 Re 都落在这个「平台」里,两者便自动相似,不必数值相等——这个区间就叫自模区(自动模化区)。它把「准则数相等」的苛刻要求松绑成「准则数同处平台区」,工程价值巨大:拖曳水池、建筑风洞里的许多常规实验,正是靠自模区才做得起、做得快。

自模区示意 无量纲阻力系数随雷诺数先剧烈变化,超过临界值后进入平台区;模型与原型的Re都落在平台上即自动相似。 雷诺数 Re(对数坐标,增大 →) 阻力系数 ζ 或 C_D 临界 Re 过渡区:ζ 随 Re 变化剧烈 自模区(平台):ζ ≈ 常数 模型与原型的 Re 都落在平台上 即自动相似,无需 Re 相等 模型(低 Re) 原型(高 Re)
图 8-4自模区:无量纲系数在 Re 足够大后进入平台。模型与原型的 Re 不必相等,只要都落在平台上,系数自动相同——这是「准则数相等」在工程上的重要松绑。

当然,自模区有边界:一旦流动发生「质变」——如边界层提前分离、转捩位置移动、涡脱落模态切换——曲线会跳出新台阶,平台失效。所以使用自模区之前必须确认:模型与原型没有跨越任何临界点。这既靠已有实验数据积累,也靠工程师的判断——「模型实验带着经验色彩」在这里又一次体现。

📐 例题 7 船模阻力的 Fr 换算:K³ 规则全链条

某船原型长 100 m、巡航速度 12 m/s。做 K = 25 的缩比模型(长 4 m)在水池中拖曳,保持 Fr 相等,测得模型总阻力 Fm = 18 N。忽略 Re 修正,求模型拖速并估算原型阻力。

思路 Fr 相等定拖速;阻力按惯性力相似换算 F ∝ ρv²L²,同种水密度相同,故 Fp/Fm = K³。

解 ① 拖速:Fr 相等 → vm = vp/√K = 12/√25 = 2.4 m/s。 验证:Frp = 12/√(9.81×100) = 12/31.3 = 0.383;Frm = 2.4/√(9.81×4) = 2.4/6.26 = 0.383 ✓。 ② 换算比:vp²/vm² = K(速度贡献一次)、Lp²/Lm² = K²(面积贡献两次)→ Fp/Fm = K³ = 25³ = 15625。 ③ 原型阻力:Fp = 18×15625 = 281250 N ≈ 281 kN(按 g = 9.81 折算约 28.7 吨的等效拖曳力)。

答:模型拖速 2.4 m/s;原型阻力约 281 kN。 ⚠ 易错点:① 是 K³ 不是 K²——速度平方贡献一个 K、长度平方贡献 K²,合起来三次方;② 这是「兴波主导、忽略粘性修正」的乐观估计,真实做法要把阻力拆成摩擦(对模型与原型分别按等价平板公式计算并修正)与剩余(按 K³ 缩放)两部分;③ 千牛换「吨力」要除以 g:281250/9.81 ≈ 28.7 吨,不是 281 吨。

衔接:从这一章到整本书 相似原理与量纲分析是全书的「通用语法」:第 7 章的莫迪图(λ 只随 Re、ε/d 变)是 π 定理的成品;第 9 章将把「尺度分析」用于 N–S 方程,导出边界层厚度 δ/L ~ Re^(-1/2) 这样的量纲结论;第 10 章的 CD 曲线、第 11 章的马赫数,全是本章思想在不同力组合下的再现。学完本章回头看:流体力学并不研究「每一个具体的流动」,它研究的是无量纲数之间的关系。

🧠 一句话记忆 模型相似 = 几何相似 + 关键准则数相等;但多个准则数常互相矛盾,只能「抓主要矛盾、做部分相似」。

本章小结

  1. 量纲是物理量的种类,单位是度量标准;基本量纲 M、L、T,其余量纲皆由其组合。
  2. 量纲一致性原则:正确方程的每项量纲必须相同,可用于检查公式对错。
  3. π 定理:n 个量、r 个基本量纲 → 可化为 (n−r) 个无量纲数的关系,把问题压缩、揭示本质结构。
  4. 重要相似准则数:Re(惯性/粘性)、Fr(惯性/重力)、Eu(压力/惯性)、Ma(流速/声速)等,每个都是两种力的比值。
  5. 模型实验靠「几何相似 + 准则数相等」实现相似,但多准则数常矛盾,工程上只能做部分相似并辅以经验修正。
  6. 相似准则数是 N–S 方程中两种力的比值:Re 称量惯性/粘性,Fr² 称量惯性/重力,Eu 称量压力/惯性,Ma² 称量惯性/弹性——判断主导力,就是给方程做量级体检。
  7. 自模区与 K³ 换算是模型实验的两大工程利器:前者把「准则数相等」松绑为「同处平台」,后者把模型测量换算回原型。

读完本章,你应该能……

思考题

1. 量纲和单位有什么区别?为什么量纲分析「与单位制无关」?

提示:量纲是种类(长度、时间…),单位是具体尺子(米、秒…)。量纲分析处理的是量纲之间的关系,无量纲数不受单位制影响,所以普适。

2. 用 π 定理解释:为什么莫迪图里 λ 只用 Re 和 ε/d 两个无量纲数就能描述?

提示:管流阻力涉及 ρ、v、d、μ、ε、Δp 等多个量,量纲分析把它们压缩成无量纲组合 λ=Δp/(ρv²L/d)、Re=ρvd/μ、ε/d,所以只需这几个无量纲数即可描述。

3. 弗劳德数和雷诺数分别代表什么力的比值?为什么它们不能同时在一个缩比船模上满足?

提示:Fr=惯性/重力,Re=惯性/粘性。缩比模型要求 Fr 相等则速度要缩小,要求 Re 相等则速度要放大,方向相反、同种流体无法同时满足。

4. 为什么说「模型实验既是科学也是经验」?

提示:理论(相似准则)提供骨架,但现实中只能部分相似,未满足的相似带来的偏差要靠经验修正系数弥补,需要工程师判断哪个力主导、如何修正。

下列哪个说法最能体现「量纲齐次性是必要不充分条件」的含义?

量纲正确只保证公式「有资格对」,不保证系数与指数也对——例如 λ = 32/Re 量纲完全正确(无量纲=无量纲),但对管流层流却是错的(正确为 64/Re)。所以量纲检查用于「排错」,数值正确性靠代入检验或实验。

什么是自模区(自动模化区)?它为什么能松绑「准则数相等」的要求?

当无量纲系数在准则数足够大后趋于常数(如钝体 CD 在高 Re 的平台),只要模型与原型的 Re 都落在平台上,系数自动相同,无需数值相等。前提是两者都不跨越任何临界点(无分离/转捩状态的突变)。

为什么像 P(kW) = T·n/9549 这样的「9549 型经验公式」要慎用?

9549 这个数字浓缩了 r/min→rad/s 与 W→kW 两层单位换算(9549 ≈ 60×1000/2π),公式只对「kW、N·m、r/min」这套特定单位成立。一旦换单位(如用 lb·ft),9549 必须重算。量纲齐次的物理公式换算系数后依然成立,而这类经验公式的常数与单位死死绑定——这是二者最实用的区别。