什么是CAE？ — 定义·优势·应用案例详解

本质上是在用计算机求解描述物理现象的"支配方程"。例如在结构分析中，是用来求解材料变形和应力的平衡方程。具体来说，考虑物体内部微小立方体（单元）所受的力，其平衡式为

因为对于复杂形状，边界条件太复杂，解析解根本不存在。比如丰田在开发新车时，发动机缸体在多个螺栓拧紧力作用下的应力分布，用手工计算绝对不可能。这就是为什么要把区域分割成小的三角形或四边形（有限单元），在每个单元上近似求解方程，然后把所有单元组合起来，形成数万到数百万的联立方程组来求解，这就是FEM（有限元法）的核心。

完全正确，误差是必然存在的。主要有"离散化误差"和"数值误差"两种。通过细化网格（h-细化）或提高单元阶数（p-细化），解会逐渐收敛到真实解。在实务中，当把网格大小缩小一半后，结果变化只有1~2%时，就认为"网格收敛"了。例如最大应力是350MPa，细化网格后变成353MPa，工程上就认为已经收敛了。

每个单元的刚度矩阵$$ [k_e] $$组装成全局刚度矩阵$$ [K] $$。应用边界条件后，求解联立方程

不是。采用直接法（例如LDLT分解）和迭代法（例如共轭梯度法）两种方法灵活应用。直接法稳定但内存占用大，迭代法适合大规模问题但收敛性是个问题。Ansys Mechanical默认使用稀疏矩阵直接求解器"Sparse"，内存不足时推荐用迭代求解器"PCG"。近年来GPU加速求解器也在增多，使用NVIDIA的CUDA库，某些线性求解器可以加速10倍以上。

要进行"增量分析"，逐步增加荷载或位移，每个增量步用Newton-Raphson法等迭代方法达到平衡。公式是，$$ [K_T]_i \{\Delta u\} = \{F\}_{ext} - \{F\}_{int}_i $$

绝对不行。首先要进行"几何检查"。在导入时查看是否有面缺失、微小边或重叠面。其次是"单位系统的一致性"。CAD用mm建模，材料数据用MPa，密度用ton/mm³吗？常见的错误是，长度用mm建模，但把钢的杨氏模量写成210而不是210,000 MPa，这样结果就会相差1000倍。

要仔细观察实物。用螺栓连接的法兰不是"完全固定"，而是"只在螺栓孔周围限制位移，考虑摩擦力"。结构物放在地面上时，可能应该用"滚动支座"，只在竖直方向限制位移。盲目使用完全固定会高估刚性，导致应力评估过低，有很大风险。ISO 10791-7等机床试验规范中有接近实物的边界条件设置示例。

那还不够。至少要通过以下检查清单：1) 反力平衡（施加荷载总和与反力总和相等吗），2) 网格依赖性验证（细化网格后结果大幅变化吗），3) 变形模式物理合理性（变形看起来不合理吗）。还要检查应力集中部位的单元形状（宽高比、倾斜角）。JIS B 9951（结构分析的验证和妥当性确认指南）很有参考价值。

价格虽然不同，但哲学和强项不同。Ansys是"针对各领域优化的专用求解器套件"。结构用Mechanical、流体用Fluent/CFX、电磁场用Maxwell，各自都是最强的独立求解器，通过System Coupling进行耦合。Abaqus分为Standard（隐式）和Explicit两种，特别是在非线性接触问题（如轮胎接地）和破坏分析上很强。COMSOL Multiphysics的卖点是"在一个环境中组建任意物理场"，还能直接编写自定义PDE（偏微分方程）。

历史因素和对汽车主要非线性问题的适应性。Abaqus一开始就专注非线性结构分析，对橡胶材料的超弹性模型、车身零件的大变形接触分析很有优势。Explicit求解器在碰撞分析上是业界标准。丰田、本田等长期使用Abaqus，积累了大量内部材料库和分析手册。但随着电动化，Ansys Maxwell（电机分析）和Fluent（电池热流体）的采用也在增加。

取决于用途。CalculiX输入格式和Abaqus类似，简单线性静力学或一般非线性分析还能用。但商用软件的高级接触算法（如面对面接触的稳定性）、丰富的材料模型（复合材料、超弹性、蠕变）和技术支持都没有。用在试制前的最终验证风险太高，但在设计初期概念验证或学生学习非常适合。最大优点是开源，能完全追踪计算过程，对教育非常有利。

这是非线性分析中最常见的错误。原因和对策是：1) 网格质量：特别是变形大的区域，单元是否有极度歪斜。宽高比大于20是危险信号。2) 材料模型：出现大变形是否用了线性弹性材料。橡胶需要超弹性模型。3) 荷载和边界条件：是否施加了现实中不可能的大荷载或约束。首先试着把荷载改为1/10再运行，看看能否收敛。在Abaqus中可以调整"稳定化"或"自动增量"参数。

主要原因有三个。第一是"初始穿透"。分析开始时网格就重叠了，求解器无法修正就会崩溃。Ansys有"Adjust to Touch"功能可以微调。第二是"接触面定义不足"。从属面网格太粗，可能被漏掉。第三是"刚体运动"。只有接触约束的零件在其他方向没有支撑会运动。检查所有零件是否有适当的约束。

很可能"刚体模式"混入了。结构不充分约束（悬空）时，不产生变形的刚体运动（平移和旋转）会作为频率为0Hz的模式出现。数值误差导致这个模式被计算为很低的非零频率，混入变形模式中。对策是施加约束来抑制所有刚体运动。常用的是三点支持（1点固定XYZ，另1点固定YZ，再1点固定Z）。查看结果时，首先检查最低阶模式的振动数是否接近0，变形是否像刚体运动。