3. 汽车结构有限元分析基础

在 Altair OptiStruct 中做汽车结构轻量化,第一步不是急着画网格,而是先想清楚“我要回答什么问题”。同一个白车身模型,问“扭转刚度够不够”和问“一阶扭转模态是多少”,对应的分析类型、边界条件、输出量完全不同。本节把 OptiStruct 常用的四类结构分析——线性静力(刚度/强度)、模态、线性屈曲——放在一起讲,重点说明它们各自解决什么问题、边界条件怎么加、求解器参数怎么设、结果怎么看,以及精度控制的关键点。

一、四类分析类型:各自回答什么问题

OptiStruct 的分析类型由 ANALYSIS 卡片和求解序列(如 STATICS、MODES、BUCKLING)共同决定。汽车结构中最常用的是:

  • 线性静力分析(STATICS):在给定载荷下求位移和应力。刚度看位移,强度看应力。适用于小变形、材料线弹性、载荷缓慢加载的工况,如白车身弯曲/扭转刚度、悬架支架强度。
  • 模态分析(MODES):求结构固有频率和振型,不加载荷,只加约束。用于判断整车/部件是否与激励源(发动机怠速、路面、电机)发生共振。
  • 线性屈曲分析(BUCKLING):求临界屈曲载荷因子。用于薄壁梁、地板、电池包下箱体等受压部件的失稳风险评估。
  • 强度校核:严格说不是独立分析类型,而是静力分析结果的判读方式——把 von Mises 应力与材料屈服强度对比,得到安全系数。

四类分析共享同一套网格和材料,但边界条件与输出请求不同。下面这张图给出从模型到结果判读的完整关系。

OptiStruct 结构分析:类型—边界—求解—判读关系图 共享模型输入 网格 + 材料 + 属性 单位制统一 (mm-t-N-MPa) 线性静力 STATICS 刚度(位移) / 强度(应力) 模态 MODES 固有频率 / 振型 线性屈曲 BUCKLING 临界载荷因子 λ 强度校核 von Mises vs 屈服强度 边界条件与载荷 静力:约束 + 集中/分布力 模态:仅约束,无外载 屈曲:约束 + 参考载荷 SPC 约束 / RBE2 刚性连接 RBE3 柔性加载 / FORCE GRAV 重力 / PLOAD4 压力 工况组合 LOAD 卡片 约束不足 → 刚体位移 过约束 → 虚假应力 求解器参数 METHOD(模态提取) EIGRL / EIGRA PARAM,AUTOSPC PARAM,COUPMASS 结果判读 位移云图 应力云图 频率列表 振型动画 屈曲因子 λ 安全系数 能量误差 约束反力 网格收敛性 单元质量检查 应力奇异识别 模态参与因子 屈曲模态形状 结果文件 .h3d HyperView 后处理 精度控制主线:单位制统一 → 网格质量(雅可比、翘曲、长宽比)→ 连接建模(焊点/粘胶/螺栓)→ 边界条件合理性 → 求解器参数 → 结果收敛性验证 常见陷阱:约束点选在非节点、RBE2 跨度过大、模态分析误加载荷、屈曲分析未加参考载荷、应力集中处网格过粗

二、边界条件与载荷工况设置

边界条件决定“结构怎么被固定”,载荷决定“结构受什么力”。两者共同构成工况(Subcase)。在 OptiStruct 中,约束用 SPC 卡片,集中力用 FORCE,压力用 PLOAD4,重力用 GRAV,工况用 SUBCASE 组织。

以白车身扭转刚度分析为例,典型设置是:在后悬架左右安装点约束全部 6 个自由度(123456),在前悬架一侧施加垂向力、另一侧施加反向力,形成绕纵轴的扭矩。扭转刚度计算公式为 K = T / θ,其中 T 为施加扭矩,θ 为前后轴间相对扭转角。位移结果从 .h3d 文件读取,取加载点垂向位移差除以轮距得到扭转角。

模态分析的边界条件要特别注意:自由模态不加任何约束,前 6 阶为刚体模态,频率接近 0;约束模态按实际安装状态约束,用于评价装车状态下的动态特性。汽车行业常说的“白车身一阶扭转模态”通常指自由-自由状态下的第一阶弹性模态,目标值一般在 40 Hz 以上(具体因车型而异)。

屈曲分析需要先施加参考载荷,求解器计算临界载荷因子 λ,实际临界载荷 = λ × 参考载荷。若 λ < 1,说明参考载荷已超过临界值,结构会失稳。薄壁件在屈曲分析前建议先做一次静力分析,确认应力水平合理。

三、求解器参数与结果解读

OptiStruct 的求解控制通过 PARAM 卡片和求解序列实现。模态分析常用 EIGRL 卡片指定提取频率范围或阶数,例如提取 0~200 Hz 内的前 20 阶。静力分析中 PARAM,AUTOSPC 可自动约束无刚度自由度,避免奇异;PARAM,COUPMASS 用于考虑耦合质量矩阵,对高精度模态分析有意义。

结果解读要抓住三个层次:

  • 位移:看最大值位置和分布趋势,判断刚度薄弱区。注意位移是相对量,与约束位置强相关。
  • 应力:看 von Mises 应力云图,重点关注高应力区和应力梯度。安全系数 = 材料屈服强度 / 最大 von Mises 应力。但要注意应力奇异——尖角、集中载荷点附近的应力会随网格加密无限增大,这类位置的应力值不可直接用于强度判断。
  • 频率与振型:看频率列表和振型动画,判断模态是否密集、振型是否局部。模态参与因子可帮助识别对整体响应贡献大的模态。
结果文件默认输出为 .h3d,用 HyperView 打开后可查看云图、动画和曲线。建议同时输出 .out 文件,里面包含求解信息、警告和错误,是排障的第一手资料。

四、分析精度控制要点

精度控制贯穿建模到后处理全过程,核心要点如下:

控制环节关键要点典型判据/做法
单位制长度、质量、力、时间单位必须自洽汽车常用 mm-t-N-MPa-s,密度 7.85e-9 t/mm³
网格质量雅可比、翘曲度、长宽比、内角雅可比 > 0.6,翘曲 < 5°,长宽比 < 5:1
连接建模焊点、粘胶、螺栓的等效方式焊点用 CWELD/ACM,粘胶用实体或壳偏置
边界条件约束位置与实际安装一致避免约束在非节点、避免过约束
载荷施加加载方式与工况匹配RBE3 分散加载,避免单点集中
收敛性网格加密后结果变化幅度位移变化 < 5%,应力变化 < 10%

实际项目中,建议对关键工况做一次网格收敛性检查:将全局网格尺寸缩小 20%~30%,重新求解,对比位移和应力变化。若位移变化在 5% 以内、应力变化在 10% 以内,可认为网格已收敛。对于应力集中区域,可局部加密后单独评估。

模态分析中若忘记删除静力载荷,OptiStruct 会忽略载荷但仍可能产生警告;屈曲分析中若参考载荷过小,λ 会很大,但实际临界载荷仍等于 λ × 参考载荷,不要误以为结构很安全。所有分析完成后,务必检查 .out 文件中的警告和错误信息。