MMS — CAE术语解说

好问题。MMS（Method of Manufactured Solutions，制造解法）是用于验证数值解法和代码的方法。首先，我们自己"制造"一个解。比如，二维热传导的温度分布可以设为

完全正确。通常的过程和MMS相反。定常热传导方程为

有几个重要原则。第一，解必须在解析意义上可以微分和积分。第二，解要反映问题的物理特性，比如非线性、对流项的主导性等。例如，验证湍流模型时，应该选择速度梯度较大的解。第三，边界条件要严格满足。比如，迪利克雷条件

比较差是第一步，但不够。关键是要检查"收敛率"。有规律地细化网格（例如，1mm, 0.5mm, 0.25mm），计算各网格下的误差范数（L2范数或H1范数）。用对数图来验证理论收敛率（比如线性元素的L2误差是

节点值差（逐点误差）也要看，但更严格的方法是用积分来评估两个函数的差。比如，L2范数是

可以，但需要高超的技巧。比如，考虑具有等向硬化的弹塑性材料。先制造一个位移场

建议从最简单的线性定常问题开始。步骤是：1) 选择一维泊松方程。2) 制造一个三次多项式解

没有直接的"MMS模式"，但可以用APDL或用户定义函数（UDF）实现。热分析中，可以用`BFE`命令逐单元定义体积热生成。更通用的方法是用`ANSYS UPFs`（用户可编程特性）在方程级别添加源项，开发自定义单元。不过这是高级用法。如果是为了验证，用开源FEM库（FEniCS、deal.II）可能更简单。

仅有小误差是不够的。最有力的证据是"收敛率"与理论值统计一致。比如，分4步细化网格，绘制误差图，其斜率应该等于预期的收敛阶（如2）。另外还要检查边界条件和守恒律（如总热流平衡）是否与网格粗度无关。ASME V&V 20规范中推荐代码验证要进行这样的系统收敛性测试。

非常合适。COMSOL本质上与MMS亲和性很强。在GUI中选"系数型PDE"界面，把制造解定义为"解析函数"。之后，在源项字段中直接输入将制造解微分得到的公式。比如，扩散方程的源项f可以写成 `-d(d(u_manufactured,x),x)-d(d(u_manufactured,y),y)` 。这比Abaqus的用户子程序（Fortran）或ANSYS的APDL简单得多。

标准的热和结构模块很难定义任意源项。主要有两个方法：一是用`DLOAD`或`DFLUX`等用户子程序逐单元定义源项；二是用`Abaqus/Standard`的`USER ELEMENT`（UEL）界面，从头创建单元矩阵和右端向量。后者灵活性极高但需要深入理解FEM，开发成本大。Dassault Systèmes的技术支持已公开了几个用UEL做验证的案例。

开源是MMS的主战场。Python库FEniCS和Firedrake都能直接把弱形式对应代码，制造解和源项的定义非常接近数学记号。比如，在变分形式中可以直接写 `f = -div(grad(u_manufactured))`。MFEM和deal.II的教程里有大量MMS收敛率验证的示例代码。研发现场通常用这些工具验证新的求解算法。

首先怀疑"边界条件设置错误"。确认数值计算中的边界条件与制造解完全一致。比如，诺依曼条件（梯度指定）的符号有没有错。其次怀疑源项实现错误。如果是自写代码，检查源项的积分（右端向量组装）有没有bug。商用软件的话，检查用户子程序中的坐标变换（全局坐标→自然坐标）和积分点权重是否正确。第三，制造解"多项式次数太高"。积分精度（高斯积分点数）可能不足，源项积分本身就有误差。

不一定是根本性错误，但肯定有问题。首先检查误差范数计算是否正确。特别是计算H1范数（梯度误差）时，数值解的梯度怎么算的（从节点值差分，还是单元内导数）。其次，网格是否足够细。初始粗网格可能还没进入渐近收敛区域，这会让斜率看起来偏低。最后，问题本身有没有奇性。制造解虽然光滑，但区域的形状（如重进型角）或材料常数的不连续可能破坏数值解的光滑性，从而降低收敛率。

不是无望，但初始条件和求解器设置要非常小心。非线性求解器（牛顿-拉夫逊法）初值离解太远会发散。MMS的优势是"已知解"，所以可以将初值设置得非常接近制造解。如果还是不收敛，这强烈提示雅可比行列式（切线刚度矩阵）的实现有误。另外，材料参数（如弹塑性模型的硬化系数）设得太极端会数值不稳定。安全的做法是先在线性弹性范围内验证，再逐步增强非线性性。