前处理 — CAE术语解说

问得好。CAE中的"前处理"是指为了使求解器能够进行计算，对几何形状进行离散化并设置物理条件的一系列工作的总称。具体来说，包括几何模型的修复、网格生成、材料物性的分配、边界条件和荷载条件的设置、接触定义、分析类型的选择等。例如，在Ansys Mechanical中，从几何到模型的工作流程基本上都属于前处理。

这正是实务中最大的陷阱。设计用的CAD数据是以制造和外观为优先的，包含分析中不需要的微小面（例如R0.1mm的圆角）和不完全接合的"接触状态"的零件。求解器无法直接理解这些信息。例如，如果厚度0.5mm的加强筋和本体被建模为0.01mm的间隙，刚性会被严重低估，计算出的振动模式会与实际完全不同。前处理是把这种"设计模型"转换为"分析模型"的工程。

根据经验法则，在不评估应力集中的区域，通常删除应力评估目标部分板厚的10%以下的圆角或孔洞。在固有值分析中，对质量或刚性贡献小于总体1%的零件可以排除。重要的是这种判断取决于分析目的。如果要评估疲劳寿命，需要再现所有的微小缺口。Abaqus/CAE的"几何编辑"工具和ANSYS SpaceCllaim的"拖拉"功能专门用于这种形状简化（去特征化）。

主要有两种算法。一种是"推进锋法（Advancing Front）"，从边界向内部逐次生成三角形或四边形。另一种是"Delaunay三角剖分"，先生成点云，然后按照Delaunay条件（三角形的外接圆内不包含其他点）连接这些点。实用的求解器通常是这两种方法的组合。例如，Altair HyperMesh的2D网格划分器以Advancing Front为基础，但为了改善网格质量，也部分利用Delaunay重新划分。

数学本质在于是否能将任意的3D区域"拓扑上映射到长方体"。复杂的形状使这种映射变得非常困难。因此，六面体网格生成通常采用"多块法"。这是将形状分解为若干个六面体块，在每个块内生成结构化网格的方法。Siemens NX Nastran和MSC Apex的某些功能试图自动化这一过程，但通常仍需要用户进行块分割（分解）。相比之下，四面体是可以填充任意形状的单形，所以算法相对容易。

主要影响数值误差，特别是"刚度矩阵的条件数"。极端的纵横比（例如100:1）会导致单元刚度矩阵的最大特征值与最小特征值的比值变大，求解器在求解线性方程组时出现数值不稳定性。具体来说，迭代求解器（如PCG法）收敛变慢，直接求解器（如稀疏LU分解）的舍入误差也会放大。允许标准取决于分析类型，但对于静弹性分析，通常目标是纵横比小于10，偏斜角小于60度。Ansys的网格指标检查可以自动评估这些值。

有既定的最佳实践。首先是"明确分析目的"。要知道什么、精度如何。其次是"几何模型的评估和修复"。检查从CAD导入的错误（间隙、重叠面、微小边），删除分析不需要的特征。然后是"确定仿真方针"：利用对称性、选择线性还是非线性、接触定义方针。在完成这些之后，才进行"网格生成"。如果切出网格后才想到"这个零件可以是刚体"，就得全部重来。

这是物理现象和计算成本之间的权衡。"粘结"是最简单的，将接触面视为完全粘合（变位连续）。适用于初期阶段粗略观察配合部位应力。"摩擦"允许接触、分离和滑动，并考虑摩擦。对于螺栓接合部位的荷载传递或制动衬片分析是必需的，但计算变为非线性，可能出现收敛问题。"无摩擦接触"忽视摩擦，仅考虑分离和滑动，对许多非线性静分析来说是平衡良好的选择。在Ansys中，将"摩擦"的摩擦系数设为0就等同于"无摩擦"。

这是最危险的想法。盲目细化会导致单元数量膨胀，计算时间变得现实不可行，而且舍入误差的积累反而可能降低精度。实务中进行"收敛性研究"。例如，关注应力集中部位的最大应力值，逐步细化网格（如2mm、1mm、0.5mm），直到结果变化在允许范围内（例如5%以内），采用该网格大小。这个过程可以用Abaqus的"网格收敛性研究"功能或自编的Python脚本自动化。

取决于用途。集成环境（Ansys Workbench、SIMULIA Abaqus/CAE、COMSOL Multiphysics）工作流程协联顺畅，学习成本相对较低。特别是要与参数研究或优化相结合时优势明显。专用工具（Altair HyperMesh、ANSA）针对超复杂装配模型（如汽车完整车身）或特定行业规范（NASTRAN、LS-DYNA输入文件）具有高度专业化的网格控制功能。航空航天业界处理数百万单元的模型时，HyperMesh的批量网格划分功能非常有用。

设计理念完全不同。CATIA和SolidWorks是"基于历史的参数化CAD"，特征的创建顺序和尺寸约束以树形结构管理。这对设计变更很有利，但不利于编辑来自其他软件的导入模型。SpaceCllaim的中心是"直接建模"，可以直接推拉面和边改变形状。不受历史限制，对任何CAD源的模型都能快速修复（去特征化、实心变空心等）。这是一种专注于"处理形状"的CAE前处理专业工具。

可以，但适用范围有限。Gmsh在脚本化网格生成（.geo文件）上表现出色，对于研究用途或特定形状模式反复出现的模型非常强大。但是，通过GUI操作进行复杂装配模型的交互式修复功能上，不如商业工具。Salome作为集成平台，涵盖几何、网格和后处理，但不能期待商用软件那样高度的自动接触定义或网格质量批量检查功能。对于降低成本同时自动化特定的反复工作，或教育用途，十分有用。

几乎100%是网格质量问题。特别是在六面体单元或二阶单元（有中间节点）中容易出现。"负雅可比"表示单元形状极度扭曲，从等参数坐标到实坐标的映射不唯一（单元自己与自己相交）。具体原因是凹陷的六面体单元或中间节点距边界较远的"凸起"单元。对策是用网格质量检查工具确认"翘曲角"和"偏斜度"，对有问题的单元进行重新网格划分。在Abaqus中，也可以通过"网格→验证"工具放松雅可比检查的阈值（增加"容差"）作为暂时对策，但这无法根本解决问题。

首先要怀疑"初始接触状态"。零件间是否相互穿透？是否有意外的微小间隙？例如，在Ansys Mechanical中，可以用"接触工具"可视化初始接触状态。穿透会导致"刚体运动"，使求解器发散。间隙会引起接触的突然发生，造成"不连续"，难以收敛。对策是在CAD阶段设置适当的间隙，或用分析软件的"接触面偏移"功能调整。其次要检查"网格不匹配"。接触面之间的网格大小差异过大，会导致接触力的传递不稳定。接触面的网格大小应保持一致，这是原则。

前提是使用64位版本的软件。在此基础上，可以分解工作流。例如，HyperMesh中可以使用"轻量化"模式，将不显示的组件从内存中卸载。也可以将几何修复和网格生成分开进行。除了对整体装配体进行网格划分，还可以逐个零件（按子系统）进行网格划分，然后再结合，这是"子结构分析"的方法。如果物理内存不够，可在装有SSD的机器上设置虚拟内存页面文件（速度会下降）。根本上，工作站的RAM应尽可能大（目前128GB以上较为理想）。