4. 仿真基础
在疲劳仿真与台架试验对标的工作中，有限元分析（FEA）是承载一切疲劳计算的“地基”。如果静力模型本身存在约束不当、接触定义错误或网格畸变，后续的疲劳寿命预测无论采用何种算法，都会与台架实测产生系统性偏差。因此，理解有限元分析的标准流程、熟悉疲劳分析工具链的分工，并掌握从静力结果向疲劳分析过渡的关键衔接点，是开展对标工作的前提。
一、有限元分析的标准流程
有限元分析并非“画网格—点求解”这么简单，它是一条从几何到结果的完整数据链。以结构疲劳评估为目标时，流程通常包含以下环节：
- 几何准备：清理CAD模型中的小圆角、螺纹孔、微小倒角等对疲劳热点不敏感但会急剧增加网格量的特征。对于焊缝区域，需保留焊趾几何或按结构应力法建立壳-实体耦合模型。
- 材料定义：静力分析至少需要弹性模量E、泊松比ν、密度ρ；若考虑塑性，还需真应力-真应变曲线。疲劳分析则额外需要S-N曲线或应变-寿命曲线。
- 网格划分：疲劳热点区域应采用二阶单元（如C3D10、SHELL181），网格尺寸建议不大于板厚的1/2；过渡区可适当放粗，但相邻单元尺寸比不宜超过2:1。
- 边界与载荷：约束需模拟台架实际装夹状态，避免过约束导致刚度虚高。载荷以台架实测谱或设计谱形式施加，静力分析中通常取单位载荷或极限载荷。
- 求解与后处理：静力求解得到位移、应力、应变场；疲劳分析则在此基础上叠加循环载荷信息，计算损伤与寿命。
需要强调的是，有限元模型的验证不能只看“是否收敛”。一个收敛的结果可能因为约束反力与台架实测差异过大而完全不可用。建议在静力阶段就提取关键测点的应变或位移，与台架静载标定数据对比，误差控制在10%以内再进入疲劳计算。
二、疲劳分析软件工具链
疲劳分析工具链通常由“通用有限元求解器 + 专业疲劳后处理软件 + 试验数据管理”三部分构成。常见组合如下：
| 环节 | 典型工具 | 核心作用 |
|---|---|---|
| 前处理 | HyperMesh、ANSA、Simcenter 3D | 几何清理、网格划分、连接与接触定义 |
| 静力求解 | Abaqus、Nastran、ANSYS | 输出应力/应变场、模态、频响结果 |
| 疲劳求解 | FE-SAFE、nCode DesignLife、Femfat | 载荷谱映射、损伤累积、寿命预测 |
| 试验对标 | MTS、Instron控制器 + 数据采集 | 实测应变、载荷、循环次数，用于修正仿真 |
工具链的衔接核心在于结果文件的兼容性。例如，Abaqus的.odb文件可直接被FE-SAFE读取，Nastran的.op2或.f06文件可导入nCode。若使用自定义脚本，则需确保应力张量、节点坐标、单元连接关系的映射关系正确。实际项目中，建议在静力分析阶段就统一单位制（如mm-N-MPa），避免疲劳软件因单位不一致导致寿命偏差几个数量级。
三、静力分析到疲劳分析的衔接要点
静力分析到疲劳分析的过渡，不是简单地把应力结果“拖进”疲劳软件，而是需要完成以下关键映射与转换：
- 载荷谱映射：静力分析通常只施加一个单位载荷或极限载荷，疲劳分析需要将实测载荷-时间历程按通道映射到有限元模型的加载点。对于多轴疲劳，还需定义各通道之间的相位关系。
- 应力组合：若静力分析采用线性叠加，疲劳软件需根据载荷谱的每个时间步，将单位载荷应力场乘以对应载荷系数后叠加，得到应力时间历程。对于非线性接触，则需在每个载荷步重新求解，计算量显著增加。
- 热点识别：静力结果中的高应力区不一定是疲劳危险点。疲劳分析需结合应力梯度、材料缺陷敏感性和载荷循环次数综合判断。例如，某点静应力虽高但循环次数极低，其疲劳损伤可能远小于一个中等应力但循环百万次的区域。
- 平均应力修正：静力分析得到的平均应力会显著影响疲劳寿命。在疲劳软件中需选择Goodman、Gerber或SWT等修正方法。对于焊接结构，通常采用结构应力法并配合主S-N曲线，此时平均应力已包含在曲线中，无需额外修正。
- 结果验证：完成疲劳计算后，应将预测寿命与台架试验的失效位置、循环次数进行对比。若偏差超过2倍，需回溯静力模型的约束、载荷施加方式或材料参数，而非直接调整疲劳算法参数。
四、从静力到疲劳的典型数据流
为了更直观地理解上述衔接关系，下图展示了一个典型的从有限元静力分析到疲劳寿命预测的数据流与工具链协作关系。图中左侧为静力分析阶段，右侧为疲劳分析阶段，中间通过结果映射与载荷谱叠加完成过渡。
上图清晰地展示了静力分析与疲劳分析之间的数据传递关系。左侧静力阶段输出的应力场是疲劳分析的“原材料”，中间的衔接要点决定了数据能否正确转换，右侧疲劳阶段则完成从载荷谱到寿命预测的最终计算。实际项目中，最容易出问题的环节往往不是疲劳算法本身，而是中间衔接处的单位制、坐标系、载荷通道映射和结果文件版本兼容性。
五、常见衔接错误与排查方向
根据工程经验，静力到疲劳衔接阶段的高频问题包括：
- 单位制不一致：静力模型用mm-N-MPa，疲劳软件默认m-Kg-Pa，导致寿命相差10⁶倍。排查方法是检查应力值量级是否合理（如钢材屈服约250MPa，若显示250e6则单位错误）。
- 载荷通道映射错误：台架有6个加载通道，疲劳软件只映射了3个，导致某些方向损伤被忽略。需核对通道编号、方向、量程与有限元加载点的一一对应关系。
- 坐标系不匹配：静力结果在全局坐标系下输出，疲劳软件按局部坐标系读取，导致应力张量方向错误。需在导出结果时明确坐标系定义，或在疲劳软件中统一转换。
- 平均应力修正方法不当：对焊接接头使用Goodman修正，而实际应使用结构应力法的主S-N曲线，导致寿命预测偏保守或偏危险。
掌握有限元分析流程、熟悉工具链分工、理解静力到疲劳的衔接要点，是开展疲劳仿真与台架试验对标的基础能力。只有在静力阶段就建立可信的模型，疲劳分析的结果才具有工程参考价值，后续的对标工作也才能有的放矢。