量纲分析与CAE模拟
量纲分析与CAE的理论基础
为什么量纲分析对CAE必不可少
老师，量纲分析与CAE有什么关系？我在大学物理中听过"单位要统一"这样的话，但在模拟现场真的会用到吗？
简单说就是模拟结果的合理性检查必不可少。只需检查输出的应力是否为Pa、热流密度是否为W/m²，仅这一项就能找出输入错误造成的一半漏洞。
真的有那么多单位错误吗？
多得很。著名的例子是1999年NASA火星气候探测器。洛克希德公司用磅力（lbf）输出数据，NASA用牛顿（N）来读取进行轨道计算。结果探测器冲入火星大气销毁。损失3.28亿美元，仅仅因为一个单位制错误。
3亿美元…！CAE中也会有同样的事吗？
而且量纲分析还有另一个强大的用途。通过Buckingham π定理将参数无量纲化，就能推导出实验的尺度变换规律。风洞试验1/10模型的结果如何外推到实物、CFD的计算域如何缩小来减少网格数量等，这些都是每天都在用的技术。
量纲的基础与量纲齐次原理
那先从基础开始吧。"量纲"具体指什么？
物理量的"种类"的标签。SI单位中有7个基本量纲。CAE中特别重要的是这4个：
- 质量 M（kg）
- 长度 L（m）
- 时间 T（s）
- 温度 Θ（K）
所有物理量都可以表示为这些基本量纲的幂的乘积。例如力：
应力是力除以面积：
量纲齐次原理（principle of dimensional homogeneity）：物理上有意义的方程，等号两边的量纲必须相等。如果不相等，这个方程就是错的。
也就是说FEM的控制方程 $[K]\{u\} = \{F\}$ 中，左边（刚度乘变位移）与右边（力）的量纲必须相等？
正是这样。验证一下。刚度的量纲是 $[K] = \mathrm{M \, T^{-2}}$（N/m = kg/s²），位移是 $[u] = \mathrm{L}$（m）。所以左边是 $\mathrm{M \, L \, T^{-2}}$ = 力的量纲。匹配了。如果这个搭不上，说明输入数据的单位制混淆了。
Buckingham π定理
你前面提到的Buckingham π定理，具体是怎样的定理？
简单说就是"物理问题所依赖的独立无量纲参数个数"的定理。
也就是说原物理问题可以完全用 $p$ 个无量纲数来描述。变量个数减少了，所以参数研究的工作量大幅下降。
具体例子看看。比如管内流的压力损失怎么分析？
好问题。管内流的压力损失 $\Delta p$ 相关的物理量列出来：
- $\Delta p$（压力损失）: $\mathrm{M \, L^{-1} \, T^{-2}}$
- $V$（流速）: $\mathrm{L \, T^{-1}}$
- $D$（管径）: $\mathrm{L}$
- $L$（管长）: $\mathrm{L}$
- $\rho$（密度）: $\mathrm{M \, L^{-3}}$
- $\mu$（粘度）: $\mathrm{M \, L^{-1} \, T^{-1}}$
- $\varepsilon$（表面粗糙度）: $\mathrm{L}$
$n = 7$ 个变量，基本量纲 $k = 3$（M, L, T），所以 $p = 7 - 3 = 4$ 个π组存在。
管的压力损失就表示为 $f = \phi(Re, L/D, \varepsilon/D)$。原来7变量的问题整理成4个无量纲参数了。CFD参数研究与其逐个尝试所有组合不同，只需改变Re、$L/D$ 和 $\varepsilon/D$ 就行。
也就是说原来7变量每个取5个值要 $5^7 = 78,125$ 个工况，但4变量只需 $5^4 = 625$ 个工况。成本降低100倍以上…！
这就是量纲分析的真正威力 — 同时降低计算成本和理清物理本质的参数。
主要无量纲数的推导
除了Reynolds数以外，CAE中还有哪些常用的无量纲数？推导方法也告诉我。
CAE中常见的无量纲数，加上推导过程来看3个。
推导：Navier-Stokes方程用代表长度 $L$、代表速度 $V$ 进行无量纲化。对流项 $\rho(\mathbf{v} \cdot \nabla)\mathbf{v}$ 的量级是 $\rho V^2/L$，粘性项 $\mu \nabla^2 \mathbf{v}$ 的量级是 $\mu V/L^2$。两者的比是：
实际上，管内流中 $Re < 2300$ 是层流，$Re > 4000$ 是湍流。CFD中要不要用湍流模型（$k$-$\varepsilon$、SST $k$-$\omega$ 等）就以此判断。
推导：能量方程无量纲化。壁面热流 $q = h(T_w - T_\infty)$ 与流体内导热 $q = k_f \Delta T / L$ 的比就是 $Nu = hL/k_f$。$Nu = 1$ 意思是"没有对流效果也一样"，$Nu = 100$ 意思是"对流把热传递快了100倍"。
也就是说Ansys Fluent的共轭传热模拟中，计算出壁面的Nusselt数就能看"对流冷却效果"？
对。而且Nu数有很多经验关联式。比如湍流管内流就有Dittus-Boelert公式：
CFD结果与这个关联式差很大，那就要考虑网格不足或湍流模型选择有问题了。这就是用量纲分析进行合理性检查。
这里 $L_c$ 是特征长度（=体积/表面积），$k_s$ 是固体的热导率。$Bi < 0.1$ 时固体内温度分布几乎均匀，不用做3D传热分析，用集总参数模型（牛顿冷却律）就够。
实际工作中怎么用Bi数？
比如铝制电子零件（$k_s \approx 200$ W/(m·K)）自然对流冷却。零件10mm、$h \approx 10$ W/(m²·K) 的话：
Bi数非常小，部件内温度梯度忽略。不用网格化做热分析。反之，塑料壳体（$k_s \approx 0.2$ W/(m·K)）的话 $Bi = 0.5$，需要3D分析才行。
分析前算算Bi数就能决定要不要做3D。太实用了…！
相似律与尺度变换
相似律的三个条件
"相似律"这个词经常听到，但没真正理解。是不是风洞试验中缩小模型用的那种？
对。相似律（相似性）就是"缩小模型的结果正确地放大到实物的条件"。有3个相似条件：
- 几何相似：形状相似（全部尺寸按比例 $\lambda_L$ 对应）
- 运动学相似：速度场相似（速度比 $\lambda_V$ 在所有点都一样）
- 动力学相似：作用力的比成立（有关的全无量纲数都相等）
第3条"动力学相似"听起来最难。全部无量纲数都相等不可能吧？
确实一般不可能。船舶的模型试验就无法同时让Reynolds数（粘性主导）和Froude数（重力波主导）相等。所以实际上采用"只让主导的无量纲数相等"的部分相似。
尺度变换规律的推导例
尺度变换规律的推导过程看一个例子。
风洞试验为例。实车（全长 $L_p = 5$ m）的1/5模型（$L_m = 1$ m）试验。Reynolds数要相等的情况：
同样空气（$\nu_m = \nu_p$）的话：
实车100 km/h，模型就需500 km/h（约Mach 0.4）风速。可压缩性影响开始显现，实际用增压风洞把密度提高，或者放弃完全Reynolds数相等而用修正。
缩放也有限制啊。阻力系数和升力系数怎么处理？
Reynolds数相等的话，阻力系数 $C_D$ 在模型和实物中相同：
模型中测得 $C_D = 0.30$，实车也是 $C_D = 0.30$。实车的阻力是 $F_{D,p} = \frac{1}{2} \rho V_p^2 A_p \times 0.30$ 来计算。这就是相似律的威力。
CAE中的模型缩放
这个思想能应用到CFD模拟吗？
当然能。特别在这些场景很有效：
- 网格独立性检查的效率化：用无量纲坐标整理结果，不同网格密度的比较容易
- 计算域缩放：利用对称性、周期性。无量纲解相等的话，计算域缩小1/4也能捕捉同样物理
- 参数研究的规划：无量纲参数空间上设计试验计划法（DOE），需要的工况数最少
- 与实验对比：用无量纲数整理，不同尺度、流体的实验数据可以直接比较
结构分析也能这样用吗？
能用。比如弹性问题，把应力无量纲化：
无量纲应力只与弹性模量 $E$、荷载参数 $P/(EL^2)$、泊松比 $\nu$、几何比有关。材料和尺寸变了，这些无量纲参数相等就能得到同样的无量纲应力分布。FEM结果能一般化。
量纲分析与CAE的实际应用
单位制的选择与统一
理论懂了。实际CAE建模时单位制怎么管理？
最开始就定规则。CAE中一般用的单位制有3种：
| 单位制 | 长度 | 质量 | 时间 | 力 | 应力 | 密度 | 主要应用 |
|---|---|---|---|---|---|---|---|
| SI | m | kg | s | N | Pa | kg/m³ | 学术、大尺度 |
| mm-tonne-s | mm | tonne | s | N | MPa | tonne/mm³ | 结构分析（最常用） |
| mm-kg-ms | mm | kg | ms | kN | GPa | kg/mm³ | 碰撞、显式求解 |
结构分析中最常见是 mm-tonne-s系。CAD一般mm做的，长度直接用mm。但密度要小心。铁的密度是 $7.85 \times 10^{-9}$ tonne/mm³。用 $7850$ kg/m³ 直接输入是典型的单位制错误。
密度的数值差10倍！ 这样就会出bug…。
所以量纲分析就是武器。结果应力是 $10^{12}$ Pa（太帕量级），对比钢的拉伸强度400 MPa，一秒内就能看出桁数错了。期望值的粗估用手算，这是最强大的调试工具。
合理性检查步骤
具体的合理性检查步骤怎么做？
结果出来了，5步必须做：
- 单位数值检查：应力是Pa还是MPa，变位移是m还是mm。数值的数量级符合期望吗
- 量纲检查：$[K]\{u\} = \{F\}$ 各项量纲一致吗。反力总和与外力平衡吗
- 守恒律检查：能量守恒、质量守恒、动量守恒成立吗。CFD的话残差够小吗
- 无量纲数检查：Re数、Nu数、CFL数在物理合理范围内吗
- 粗估对比：梁论、平板近似、Moody线图等手算的粗估值能比较吗。差超1数量级就要找原因
第5条的"粗估"，我听资深工程师说过"感觉"，就这个吗？
就是这个。比如悬臂梁的自由端挠度用 $\delta = PL^3/(3EI)$ 粗算。FEM结果是这个值的10倍，怀疑约束条件错了；是1/1000，怀疑荷载输入有问题。"FEM结果对不对用手算来判定" — 这是专业CAE工程师的标准做法。
不同工具的单位制设置
各个求解器中单位制怎么设置，有什么要注意的地方？
求解器各有各的特点，汇总一下：
| 求解器 | 单位制处理 | 注意事项 |
|---|---|---|
| Abaqus | 无单位制概念（用户责任统一） | 输入文件（.inp）没有单位头，文档必须写。密度、弹性率、荷载都要整合，用户全责 |
| Ansys Mechanical | Workbench可选单位制 | GUI设定的单位制与APDL命令数值一致性。MP材料库固定SI |
| Nastran | 无单位制概念（用户责任） | PARAM, WTMASS质量变换系数设定可能。BDF文件注释中单位写入是最佳实践 |
| OpenFOAM | SI固定 | 全m-kg-s。CADmm取入的几何一定要缩放（scaleMesh） |
| COMSOL | GUI上管理单位制 | 方程视图中可查次元。不同单位的分量组合时自动变换确认 |
| LS-DYNA | 无单位制概念（用户责任） | 碰撞分析mm-kg-ms系实质标准。重力加速度用 $9.81 \times 10^{-3}$ mm/ms² 常忘 |
软件中的量纲管理
各求解器的单位制处理
求解器单位处理这么分散，自动检查的功能没吗？
Abaqus、Nastran、LS-DYNA这种"单位自由"的求解器自动检查不了。全部用户責任。模型作成強勧。感：
| 物理量 | 记号 | 采用单位 | 输入值 | 量纲检查 |
|---|---|---|---|---|
| 长度 | $L$ | mm | — | $\mathrm{L}$ |
| 质量 | $m$ | tonne | — | $\mathrm{M}$ |
| 时间 | $t$ | s | — | $\mathrm{T}$ |
| 弹性模量 | $E$ | MPa | $2.1 \times 10^5$ | $\mathrm{M \, L^{-1} \, T^{-2}} ✓$ |
| 密度 | $\rho$ | tonne/mm³ | $7.85 \times 10^{-9}$ | $\mathrm{M \, L^{-3}} ✓$ |
| 荷载 | $F$ | N | $1000$ | $\mathrm{M \, L \, T^{-2}} ✓$ |
| 应力（输出） | $\sigma$ | MPa | — | $\mathrm{M \, L^{-1} \, T^{-2}} ✓$ |
这个表格一开始就做好，团队所有人都用同样单位制。
无量纲化功能
求解器有无量纲化的功能吗？
有些求解器有组建无量纲化功能：
- COMSOL：方程视图显示各项量纲。不一致就警告。非常方便
- Ansys Fluent：边界条件能用无量纲数（Re、Ma等）设定。Reference Values基准值设定
- OpenFOAM：dimensionedScalar型保持量纲信息。量纲不一致时
特OpenFOAM的量纲检查機能優秀。例：
dimensionedScalar nu("nu", dimViscosity, 1e-6); // m²/s
dimensionedScalar rho("rho", dimDensity, 1000); // kg/m³
// nu * rho dimViscosity * dimDensity = 动粘度 × 密度 = 粘度的次元
// nu + rho 、次元不一致！编译时就量纲检查，太棒了。Abaqus也有这种功能就好…。
Abaqus、Python自分作手。.inp、材料参数的桁数、荷重応力的期待値比较自動化企業多。量纲分析「工具頼」「自分的物理感覚判定」力鍛重要。
先进主题
分数量纲与分形分析
量纲分析有什么最新研究？
有趣的领域不少。先说分数阶微积分（fractional calculus）。古典量纲分析是整数量纲（L¹、T²等），但分数阶导数用上就能自然地模化粘弹性材料的记忆效应和异常扩散。比如Creep行为：
$\alpha = 1$ 通常的弾性体、$\alpha = 0$ 完全粘性体。中间的 $\alpha$ 粘弹性表现。 [t^α] 的量纲 $\mathrm{T}^\alpha$ 、無次元数的構成複雑。
机器学习与量纲分析的融合
最近流行的机器学习与量纲分析有关系吗？
大有关系。最近注目的Dimensionally Consistent Neural Networks。神经网络输入無次元数限定、出力無次元量、物理法则整合模型構築。
例、100案例的CFD数据空力係数予測作場合：
- 素朴的方法：速度、密度、粘度、长度、角度的入力 → 单位変精度崩壊
- 量纲分析：Re、Ma、迎角α入力 → 単位依存汎化性能
PINN（Physics-Informed Neural Networks）、损失函数次元整合性的制約入、学習数据少物理的妥当予測。
200年以上前的理論最前沿AI技術活躍。基础大事…。
量纲分析与CAE的问题排查
由单位制错误引起的典型错误
先生也单位制错过吗？
年轻时几次了。典型的方式总结一下：
| 现象 | 原因 | 识别方法（量纲分析） | 对策 |
|---|---|---|---|
| 应力 $10^9$ 倍大 | 长度mm输入，密度用kg/m³ | 应力期望值（MPa量级）对比，桁位异常 | 单位系检查表确认密度变换 |
| 固有振动数 $10^3$ 倍高 | 质量 $10^6$ 倍小（kg/m³→tonne/mm³变换忘） | $f \propto \sqrt{K/M}$，M小则f大 | 总质量检查（手算对比） |
| 热分析温度发散 | 热导率单位W/(m·K)直接输入mm系 | 计算Bi数,检查合理性 | $k$转换成mW/(mm·K) |
| CFD流速异常 | 动粘度输入错（SI值输入mm系） | 逆算Re数,检查流场合理 | 从Re数逆算输入参数,对比 |
| 碰撞分析部件穿透 | LS-DYNA时间单位错（s与ms） | 动能粗估对比 | 检查接触力时间序列，$F \cdot \Delta t$量纲检查 |
固有振数那个，超常见啊。密度一个桁错就振动数1000倍…。
一番多。対策簡単。模型的総質量手算比较。体積×密度的概算質量求解器出力的質量一致、少密度入力正。「」呼、的企業的第一项目。
量纲检验检查表
最后实务使检查清单作。
、分析模型的作成時結果確認時使检查清单总结：
的检查清单、印刷貼。量纲分析地味見、仿真的品質根本支技術。
的通。高価求解器使、入力数据的単位間違出。"Garbage in, garbage out" 防最初的防壁次元分析。物理的基本大切。
帮助
详细
错误