仿真 — CAE术语解说
理论与物理 — 基本概念、支配方程
理解CAE基本术语
「网格」就是简单地把模型切成细小部分吗？我听说网格越细密精度越高。
这是一个很大的误区。网格是计算的「基础栅格」，其质量决定了结果的好坏。仅仅细密是不够的，比如在应力集中部分应该细化，而在变形很小的部分则可以粗化的「自适应网格」才是基本方法。Ansys的用户手册推荐使用根据曲率和狭小间隙自动生成细网格的「Curvature」和「Proximity」功能。
「支配方程」这个词经常出现，CAE软件具体求解什么方程呢？
取决于分析领域。结构分析中，表示弹性体变位与应力关系的基本方程是
「收敛」是指计算结束吗？为什么迭代增加时会说「收敛性差」呢？
收敛是指反复迭代后解的变化量在允许误差（如残差范数小于1e-6）范围内，并不是指计算结束。非线性分析（接触、塑性）中，将荷载分成小步骤进行「增量分析」，如果某个步骤的迭代次数超过50次仍未收敛，就判定为「收敛性差」。这通常是由接触条件的急剧变化或材料不稳定行为引起的。
「边界条件」是将现实数值化吗？除了「完全约束」，还有其他常用的设置吗？
正是的，边界条件是将现实的约束和环境进行模型化。例如「对称条件」是当模型具有对称性时，在对称面上约束法向变位和面内旋转，这样可以将计算成本降低到1/2或1/4。「远场边界条件」在CFD中常用，用来设定进口速度或压力出口。Ansys Fluent中设为「velocity-inlet」或「pressure-outlet」。设置错误会产生非物理的流动，需要注意。
数值求解方法与实现 — FEM/CFD离散化、求解器设置
求解器背后发生的事
「直接法求解器」和「迭代法求解器」的区别是什么？软件自动选择的话不用在意吗？
绝对要注意。直接法（如PARDISO、MUMPS）消耗大量内存但数值上稳健，对于包含接触和不连续的非线性问题的每个步骤都有效。迭代法（如CG法、GMRES）在大规模问题中内存效率高，但如果刚度矩阵病态（条件数差），就可能不收敛。Abaqus/Standard默认使用直接法，但对于大规模线性弹性问题可选择迭代法求解器。在32GB内存的PC上求解百万节点问题时，直接法可能会因内存不足而失败。
有很多「单元类型」，一阶单元和二阶单元的实际区别是什么？二阶单元总是更优吗？
一阶单元（线性单元）仅在角点有节点，二阶单元在边的中点也有节点。要表达弯曲和应力集中，二阶单元的精度绝对更高。但计算成本也更高。经验法则是，要获得相同精度，一阶单元所需的节点数是二阶单元的约8倍（三维情况）。不过接触分析中二阶单元的接触面处理变得复杂，收敛性恶化。因此Abaqus/Explicit（动态分析）通常推荐使用一阶单元作为默认。
「离散化误差」和「舍入误差」有什么区别？哪个对结果影响更大？
离散化误差是将连续体近似为有限个单元所产生的本质误差，细化网格可以减小。舍入误差是计算机以有限位数表示浮点数所产生的误差。通常离散化误差要大好几个数量级。例如，用一阶单元计算悬臂梁末端挠度，相对理论值可能有10%以上的误差。而舍入误差在双精度计算中约为1e-16的数量级。不过当刚度矩阵条件数极差（如1e12以上）时，舍入误差会被放大产生「数值不稳定」，导致解没有意义。
「隐式」和「显式」求解的选择依据只是时间步长吗？
时间步长是重要区别，但本质是求解的稳定性。显式法（Abaqus/Explicit、LS-DYNA）条件稳定，时间步长
实践指南 — 工作流程、检查清单
如何实施可信度高的分析
「网格独立性检验」怎样高效地做？每次都改变网格计算3-4次太耗时了。
从一开始就均匀细化全部网格是低效的。实务中，先用「粗网格」把握整体行为，确定应力高的区域和关心的区域。再对该区域进行仅在该处细化的「局部加密」。使用Ansys Workbench的「Convergence」工具，可以自动将关心量（最大应力等）变化率控制在2%以内，每次将单元尺寸减半。最后比较中等网格和最细网格的结果，验证差值在允许范围内。这样2-3次计算就够了。
定义材料的「塑性」时，应力-应变曲线的输入点越多越好吗？能直接用实验数据吗？
直接使用实验数据会引起数值振荡。求解器对数据点间进行线性插值，输入点过多会导致斜率不连续变化，收敛性恶化。特别是屈服点附近数据要多，之后应该平缓变化地间隔。最稳定的方法是拟合应变硬化规律（如Swift规律
「验证」分析结果时，数值和实验必须完全一致吗？允许多少误差？
完全一致在现实中不可能，要根据目的设定允许误差。汽车业结构耐久性分析中，应力集中部最大应力±10%、整体变形模式和固有振动数±5%以内是常见目标。但这是在材料常数、边界条件、网格都设置正确基础上。实验本身也有误差（应变片贴装误差、测力传感器精度等）。重要的是趋势是否一致（改变参数时的变化趋势是否吻合）。ASME V&V 20和JSME S 012规格中规定了验证和确认指南。
分析报告中必须记载的「前提条件」具体是什么？
最少需要以下5项：1. 使用软件和版本（如Ansys Mechanical 2023 R2）、2. 单元类型和总节点/单元数（如二次四面体单元（SOLID187），约50万节点）、3. 材料模型和常数（如弹性：E=210GPa、ν=0.3，塑性：双线性硬化、屈服应力355MPa）、4. 边界条件和荷载条件详情（约束位置、荷载值和方向）、5. 收敛判定准则（如力的残差范数1.0e-5）。没有这些，第三者无法重现和评价结果。汽车业还从FMEA（故障模式影响分析）角度要求对预想外荷载工况的考察。
软件比较 — Ansys/Abaqus/COMSOL等
工具特性和选择指南
Ansys Mechanical和Abaqus/Standard都是通用结构分析求解器，决定性区别是什么？
核心数值求解哲学不同。Abaqus在非线性分析中非常强大，特别是「接触」与「材料非线性」的组合表现出色，因此在汽车橡胶衬套和螺栓连接部分析中市场占有率高。另一方面，Ansys Mechanical以「工作流集成」和「多物理耦合」见长。在Ansys Workbench上容易将流体（Fluent）与结构（Mechanical）耦合，可求解热应力和流固耦合（FSI）。HPC（高性能计算）许可证的灵活性也通常是Ansys更优。预算允许的话，根据用途使用不同工具或两者兼备是理想状态。
听说COMSOL Multiphysics是「基于方程」的，和Ansys、Abaqus根本上有什么区别？
COMSOL的核心是「弱形式」，用户可直接输入偏微分方程（PDE）。相比之下Ansys和Abaqus提供「物理接口」（结构、流体、电磁等），COMSOL通过「PDE组合」来实现。这使得学术上新的物理模型和复杂耦合现象（如电场诱导的化学反应及其热变形）能在一个环境中建模。不过在大规模三维结构分析和复杂接触问题上，相比Abaqus这样的专化求解器，计算效率倾向于下降。网格生成功能也比Ansys和Altair HyperMesh等显得生硬。
「开源CAE」免费能用，实用水平吗？比如CalculiX和Code_Aster。
特定领域已达到实用水平。CalculiX与Abaqus输入文件格式（.inp）兼容，可进行线性、非线性结构分析。但GUI很差，前处理、后处理需要独立工具（PrePoMax、Gmsh、Paraview）。Code_Aster由法国电力公司（EDF）开发，在破坏力学和土壤分析中有强大功能。两者都缺乏有偿软件那样的高级接触算法和求解器优化，计算时间往往较长。不过教育目的或作为开发自有求解器的基础极其有用。考虑到商用软件年许可费数百万元，性价比很高。
「专用求解器」比通用求解器优在什么地方？比如MSC Nastran或LS-DYNA。
针对特定问题，算法被调优到极限。MSC Nastran在航空航天业长期使用，线性静解、特征值分析（振动、失稳）的速度和精度，特别是大规模模型处理能力备受推崇。其输入形式（Bulk Data Deck）是业界事实标准。LS-DYNA是显式法的专家，在碰撞、爆炸、冲压成形等极端非线性动态问题上，接触算法和材料模型丰富度无与伦比。这些专用求解器内置许多特殊数值手段来处理通用求解器困难的「病态问题」。但许可证昂贵，学习曲线陡峭。
故障排除 — 常见错误和对策
计算停止时的处理方法
非线性静力分析出现「不收敛」错误。首先应该检查什么？
首先查看日志文件，确认「哪个自由度（U1、U2、U3、ROT等）不收敛」。变位不收敛则可能是边界条件不当（存在刚体运动）。也可能初始状态有问题。接触条件是否「突然」施加，检查接触面是否有微小初始间隙或重叠。Abaqus中打开「稳定化」或「自动步长」，Ansys中引入「弱弹簧」抑制初始数值不稳定。荷载过大时，将荷载减为1/10试验，逐步增加的「阶段加载」也有效。
出现「负雅可比」或「过度畸变」错误，原因是什么，怎样修复？
单元发生极端变形，数学映射不再成立的状态。主要原因有三：1. 网格品质差（扭斜度大、长宽比大）、2. 材料模型不当（例如大应变不考虑塑性，仅用弹性）、3. 荷载/变位条件不现实。对策是首先改善网格
相关主题
帮助
更详细
错误