3. 汽车结构有限元分析基础
在 Altair OptiStruct 中做汽车结构轻量化,第一步不是急着画网格,而是先想清楚“我要回答什么问题”。同一个白车身模型,问“扭转刚度够不够”和问“一阶扭转模态是多少”,对应的分析类型、边界条件、输出量完全不同。本节把 OptiStruct 常用的四类结构分析——线性静力(刚度/强度)、模态、线性屈曲——放在一起讲,重点说明它们各自解决什么问题、边界条件怎么加、求解器参数怎么设、结果怎么看,以及精度控制的关键点。
一、四类分析类型:各自回答什么问题
OptiStruct 的分析类型由 ANALYSIS 卡片和求解序列(如 STATICS、MODES、BUCKLING)共同决定。汽车结构中最常用的是:
- 线性静力分析(STATICS):在给定载荷下求位移和应力。刚度看位移,强度看应力。适用于小变形、材料线弹性、载荷缓慢加载的工况,如白车身弯曲/扭转刚度、悬架支架强度。
- 模态分析(MODES):求结构固有频率和振型,不加载荷,只加约束。用于判断整车/部件是否与激励源(发动机怠速、路面、电机)发生共振。
- 线性屈曲分析(BUCKLING):求临界屈曲载荷因子。用于薄壁梁、地板、电池包下箱体等受压部件的失稳风险评估。
- 强度校核:严格说不是独立分析类型,而是静力分析结果的判读方式——把 von Mises 应力与材料屈服强度对比,得到安全系数。
四类分析共享同一套网格和材料,但边界条件与输出请求不同。下面这张图给出从模型到结果判读的完整关系。
二、边界条件与载荷工况设置
边界条件决定“结构怎么被固定”,载荷决定“结构受什么力”。两者共同构成工况(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% 以内,可认为网格已收敛。对于应力集中区域,可局部加密后单独评估。
.out 文件中的警告和错误信息。