3. 物理模型简化误差

在仿真误差的整条传递链上，物理模型简化误差位于最上游，也最难通过加密网格或收紧求解器残差来消除。数值误差可以靠更细的离散和更严的收敛判据压缩，而模型误差一旦引入，后续所有计算都只是在“忠实求解一个被简化过的物理问题”。因此，评估模型误差的关键不是问“算得准不准”，而是先问“我求解的方程、材料关系和边界条件，是否足以支撑我要回答的工程问题”。

一、三类简化误差的来源与识别

控制方程选择误差来自对物理机制的取舍。例如可压缩流动中，若马赫数约 0.25、温度变化不大，用不可压 N-S 方程通常可接受；但当局部出现激波或强热耦合时，忽略密度变化与能量方程会直接丢失压力—密度耦合。结构分析中，把三维实体简化为梁或壳，本质是引入 Kirchhoff 或 Mindlin 假设，当厚度与跨度之比超过约 1/10 时，剪切变形不可忽略，薄壳理论就会产生系统性偏差。识别方法是做量级分析：用特征数（Re、Ma、Pe、Biot 数、长细比）判断被忽略项相对主导项是否真的小。

本构关系近似误差来自材料模型的降阶。线弹性假设忽略塑性、蠕变与损伤；理想弹塑性忽略硬化；Drucker-Prager 或 Mohr-Coulomb 用于岩土时对围压相关性的描述有限；湍流模型中 k-ε 在强分离、强旋流区会高估湍动能，而 LES 的亚格子模型又依赖网格尺度。这类误差往往表现为“趋势对、幅值偏”，且随载荷水平非线性放大。

边界条件理想化误差最容易被低估。把支座写成完全固支，实际可能是带柔度的弹性约束；把入口写成均匀速度，实际来流有边界层与湍流度；把外表面当作绝热，实际存在对流与辐射。理想化相当于人为改变了系统的刚度、阻尼或能量交换路径，误差会直接进入响应幅值与频率。

二、影响评估：用“逐项放松”量化模型误差

评估模型误差最实用的做法是逐项放松简化假设，观察目标响应量的变化幅度。以某承压管道三通为例，目标量是峰值等效应力，工程判据为相对偏差不超过 5%。基准模型采用线弹性、内压 6 MPa、两端固支。第一步把固支改为等效弹性约束（依据法兰刚度估算，轴向刚度约 2×10⁶ N/mm），峰值应力下降约 8%；第二步把线弹性改为双线性等向硬化（屈服 235 MPa、切线模量 2 GPa），在局部应力集中区峰值应力重分布，下降约 12%；第三步加入温度载荷 80 ℃ 引起的热应力，峰值上升约 15%。三项叠加后与基准相差超过 20%，说明原简化模型不足以支撑该判据，必须至少保留弹塑性与弹性约束。

这里的关键是把“模型误差”变成可比较的数值：固定网格与求解器设置，只改物理假设，差值即模型误差的近似估计。若放松后响应变化小于判据，简化可接受；若大于判据，则需升级模型或修正判据本身。

三、控制策略与建模检查清单

控制模型误差的核心思路是“按需建模、留有余量、记录假设”。可参考以下检查清单逐项确认：