力学仿真(结构有限元分析)是CAE领域最成熟的应用方向。从线性到非线性，聊聊力学仿真中的核心概念、常见误区与工程实践经验。
Young woman scientist using tablet pc in her hands. Mixed media
1. 力学仿真的本质
有限元分析的核心思想其实不难理解：把连续的求解域离散成有限个单元，在每个单元内假设一个近似的位移场，然后通过变分原理推导出每个单元的刚度矩阵，最后组装成全局方程组求解。这个过程说起来简单，但里面每一个环节都可能引入误差。
我刚学有限元的时候，把软件当黑盒用——输入几何、材料、载荷，点一下”求解”，出来一个云图，就觉得完事了。后来因为一个项目的结果和实验差了40%，才开始认真去翻有限元法的教材，搞清楚每个环节在做什么。这次”回头补理论”的经历让我养成了一个习惯：用任何一个仿真工具之前，先搞清楚它背后的数值方法大概是怎么回事。
2. 线性静力学：大多数人的起点
线性静力学分析假设：材料是线弹性的、变形是小变形(几何线性)、边界条件不随变形改变(接触除外)。在这三个假设都满足的前提下，有限元方程的求解是直接的——只需要解一次线性方程组，不需要迭代。
但实际工程问题经常突破这些假设。比如橡胶类材料，应力-应变关系明显非线性;薄壁结构在大载荷下会出现大变形;螺栓连接处的接触状态会随着加载过程改变。遇到这些情况，线性静力学分析的结果就只能作为粗略参考，不能用于最终设计决策。我的经验是：在做一个新类型的结构分析时，先跑一个线性静力学作为baseline，再用非线性分析看两者的差异有多大——如果差异在5%以内，线性结果可以接受;如果超过15%，必须以非线性结果为依据。
3. 非线性分析的三大来源
非线性来源有三个：材料非线性(塑性、超弹性、蠕变等)、几何非线性(大变形、屈曲)、接触非线性(接触状态改变)。一个实际工程问题可能同时包含这三种非线性，这时候求解就变成了迭代过程——需要用牛顿-拉夫逊方法或者改进的求解算法，逐步逼近正确解。
接触非线性的处理是有限元分析里最棘手的部分之一。我在做一个螺栓连接分析时，初始接触状态设置不当，导致计算收敛极为困难。后来学到一个技巧：先用一个很小的载荷做一个”接触初始化”分析步，让接触对正确建立接触关系，再逐步加载到目标载荷——这个分步加载的策略，显著改善了收敛性。
4. 动力学分析：不只是静力学加了个时间轴
很多做静力学分析的人，第一次做动力学分析的时候会以为”只是加了个时间轴”。实际上动力学分析的方程组里多了惯性项和阻尼项，求解难度和研究目标都和静力学有本质区别。
模态分析(Modal Analysis)是动力学分析的基础——它告诉你结构的固有频率和振型。如果外载荷的频率接近某个固有频率，就会发生共振，振幅可能超出设计预期。我在做一个悬臂结构的设计时，模态分析显示第二阶固有频率在工作转速范围内，后来增加了刚度把固有频率移出了工作范围，才解决了振动问题。
谐响应分析(Harmonic Response)用来研究结构在正弦载荷下的稳态响应;瞬态动力学分析(Transient Dynamics)则用来研究任意载荷历程下的动态响应。两者的选择取决于载荷的性质：如果载荷是周期性且频率固定的，谐响应分析足够;如果载荷是冲击、爆炸这类瞬态载荷，必须用瞬态分析。
5. 结果验证：仿真不是算出来就完事
仿真结果的可信度验证，是很多工程师忽略的环节。我在团队里推动的一个做法是：每个重要的仿真项目，都必须有验证环节——要么是和实验数据对比，要么是和解析解(简单几何下有理论解的那种)对比，实在没有对比条件的，至少要做网格无关性验证。
一个具体的验证案例：一个悬臂梁的弯曲问题，我用ABAQUS做了精细化仿真，最大弯曲应力结果是285 MPa。然后用材料力学里的悬臂梁公式手算了一遍，结果是292 MPa，两者差异不到3%。这个对比让我对仿真模型的可靠性有了信心。如果没有这一步验证，我始终会对那个285 MPa的结果心存疑虑。
