1.3 惯性系、单位制与量纲分析

1.3 惯性系、单位制与量纲分析

本节摘要：方程写对只是半成品，数字要从合格的参考系、统一的单位里读出来才算数——惯性系是牛顿第二定律的舞台资格证，SI 单位制是全行业公认的刻度，量纲分析则是动笔前的免费体检。本节给出三件工具的用法与常见事故现场。

先决问题：刻度没校准，解剖全白做

对象定了（1.1）、法则定了（1.2），还剩最后一道术前核对：在哪个参考系里观察，用哪套单位记录，以及——有没有更便宜的办法在细算之前先掂出答案的量级。这一节承接前两节收尾本章，也是全册反复回访的工具间：第三章标约束反力时要默认惯性系，第七章在非惯性系里补惯性力，第八章的极限检验干脆就是量纲分析的近亲。

惯性系这个词听起来玄，判据其实朴素：牛顿第一定律在其中精确成立的参考系。地面算一个近似的——地球自转带来的偏差在多数工程问题里可忽略，但摆钟的长周期漂移、远程火炮的落点偏移就藏不住它了；地心参考系更好，日心参考系更好上加好。绝对严格的惯性系并不存在，工程的做法是按精度需求选"够格"的舞台：教室里算滑块用地面系毫无问题，设计傅科摆就必须把地球自转请回方程。

一、单位制：统一刻度与工程换算

SI 制把力学刻度归到三个基本量：米、千克、秒，力的单位牛顿是导出单位——使一千克质量产生一米每二次方秒加速度的力。工程现场还游荡着一堆历史遗留单位：千克力（kgf）、工程大气压（atm）、毫米汞柱，与 SI 制的换算系数必须形成肌肉记忆，因为单位事故从不挑时间。当年一颗火星气候探测器因地面软件用了英制单位而航天器按公制解算，数亿美元的使命在一瞬间变成流星——单位核对不是官僚流程，是保命动作。

# 常用力学单位的 SI 换算自查表
kmh_to_ms = lambda v: v / 3.6        # 千米每小时 转米每秒
print("72 km/h =", kmh_to_ms(72), "m/s")          # 20.0
kgf = 9.80665                       # 1 千克力定义为标准重力下的重量
print("5.0 kgf =", 5.0 * kgf, "N")                # 49.03325
atm = 101325.0                      # 1 标准大气压 Pa
print("2.0 atm =", 2.0 * atm / 1000.0, "kPa")     # 202.65
mmHg = 133.322                      # 1 毫米汞柱 Pa
print("760 mmHg =", 760 * mmHg, "Pa")             # 101324.72（与 1 atm 仅差舍入）
# 典型事故自查：把 kgf 当成 N 代入方程，力会放大 9.8 倍
print("误把 10 kgf 当 10 N，误差倍数 =", kgf)      # 9.80665
​

这张小表的正确用法是"每天扫一眼"而不是"考试前背一夜"。尤其最后两行：kgf 与 N 差将近十倍，是工程单位事故里出现频率最高的一对。

二、量纲分析：不花钱的预检

每个物理量都由基本量纲拼成：力的量纲是质量乘加速度，记作 L M T 的组合 L M T⁻²。量纲齐次性是一条铁律——任何正确的物理等式，两端量纲必须相同。这条铁律反过来用，就成了 π 定理：若一个问题涉及 n 个物理量、含 k 个基本量纲，则物理规律必可写成 n 减 k 个无量纲数之间的关系。很多场合它直接把未知函数的形态框死，例如单摆小摆幅周期只可能与摆长 l、重力加速度 g 有关，量纲上凑遍所有组合，唯一的选择是 T 正比于 l/g 的平方根的二倍π倍以内的常数——精确公式里有 2π，但"平方根"这个骨架由量纲分析白送。

量纲分析给的是骨架不是血肉：它告诉你 T 依赖 sqrt(l/g)，但给不出前面的系数。系数要靠实验或完整求解。便宜之处在于：一旦骨架错了（比如有人声称周期正比于摆长本身），量纲体检当场判死，省下全部计算。

# 单摆周期：量纲骨架的数值验证
import math
g = 9.8
for L in [0.25, 1.0, 4.0]:
    T = 2 * math.pi * math.sqrt(L / g)
    print(f"L={L:5.2f} m  T={T:.4f} s")
# 输出：
# L= 0.25 m  T=1.0030 s
# L= 1.00 m  T=2.0061 s
# L= 4.00 m  T=4.0122 s
# 摆长翻两番，周期恰好翻一番——与 sqrt(L) 骨架完全一致

# 无量纲数预检：简支梁最大挠度公式的量纲核对
# 挠度 f = C * F * L^3 / (E * I)，C 为无量纲系数（取决于支承与加载方式）
F, L = 1000.0, 2.0
E, I = 2.1e11, 8.33e-6     # 钢的弹性模量 Pa，惯性矩 m^4
f = F * L**3 / (E * I)
print("系数 C 取 1/48 时挠度 =", f / 48, "m")     # 0.0003968...
​

第二段代码顺带演示无量纲化的另一次要好处：E 乘 I 与 F 乘 L 三次方相除后剩下长度量纲，等式两端量纲自动闭合——凡是量纲闭不上的"公式"，连检算的资格都没有。

三、无量纲数：给病例贴标签

π 定理最常用的产出，是一批有名字的无量纲数——雷诺数、弗劳德数、马赫数。它们不是装饰品，而是给力学病例贴的分类标签：两个尺寸悬殊的系统，只要这批无量纲数彼此相等，力学上就是"同一种病"，实验室里的小模型因此能代替大原型做试验。船模水池把模型速度按弗劳德数相似律缩放，风洞把飞机模型按雷诺数相似律吹风，背后的账全部由量纲分析预先算好。

雷诺数是惯性项与黏性项之比，弗劳德数是惯性项与重力项之比。下面用一段数值演算演示"同一种病"怎么从标签里读出来：同一个船模在模型尺度与原型尺度下，要满足弗劳德数相等，速度该缩多少、雷诺数又差多远。

# 船模试验的相似律：弗劳德数相等决定拖曳速度，雷诺数揭示黏性失真
import math
L_model, L_ship = 2.0, 120.0          # 船模与原型的水线长，m
g = 9.8
Fr = 0.30                              # 设计弗劳德数（兴波阻力敏感区）
# 弗劳德数 Fr = v / sqrt(g*L)，两边 Fr 相等 → v 与 sqrt(L) 成比例
v_model = Fr * math.sqrt(g * L_model)
v_ship  = Fr * math.sqrt(g * L_ship)
print("船模拖速 =", round(v_model, 3), "m/s")     # 1.328
print("原型航速 =", round(v_ship, 3), "m/s")      # 10.288
rho, nu = 1000.0, 1.0e-6              # 水的密度与运动黏度
Re_model = v_model * L_model / nu
Re_ship  = v_ship  * L_ship  / nu
print("船模雷诺数 =", Re_model)                    # 2656313.2...
print("原型雷诺数 =", Re_ship)                     # 1234542830.4...
print("黏性相似失真倍数 =", Re_ship / Re_model)    # 464.758...
​

输出里那四百多倍的雷诺数差距，就是模型试验绕不开的"黏性失真"：兴波阻力按弗劳德律相似，摩擦阻力却因雷诺数悬殊不能直接换算，工程上要靠经验摩擦线分项修正。一个量纲分析同时说清了模型试验"能信什么、信到什么程度"——这正是它比一沓盲算公式值钱的地方。

四、实操清单与事故现场

⚠️ 常见坑：量纲正确不等于公式正确。2π 与 1/48 这类系数、以及函数形式里的加减号，量纲分析一概管不着——它只能判死刑，不能发通行证。

要点复盘：

• 惯性系有资格审核：第一定律在其中成立才算数，地面是近似成员，精度需求越高舞台越要挑；
• 单位事故是高频事故：kgf 与 N 差九倍八，火星探测器的教训值得刻在每张草稿纸角上；
• 量纲齐次是铁律：两端量纲不闭合的等式必然有错，闭合也只是及格；
• π 定理压缩问题：n 个量 k 个基本量纲，规律落在 n 减 k 个无量纲数上，骨架白送、系数另付。

器械间到此清点完毕：对象、法则、刻度三件套全部就位。下一章开进运动学——先把"怎么动"描述清楚，原因分析的手术留到再下一章动刀。