2. 有限元仿真基础回顾
在相关性分析中,有限元模型不是“算出一个结果”就结束,而是要与试验在同一物理量、同一位置、同一工况下逐点比对。因此,建模的每一步都要问一句:这个简化会不会改变我关心的响应?本节围绕单元类型、网格策略、边界条件与载荷施加四条主线,说明面向相关性分析的建模原则。
一、先明确“相关性目标”再建模
相关性分析通常关注三类量:位移/变形(如某测点挠度)、应变/应力(如应变片方向的正应变)、模态参数(固有频率与振型)。目标不同,建模精度分配完全不同。例如关注薄板弯曲振型时,壳单元的横向剪切与转动自由度必须保留;而只关注整体刚度时,可用更粗的网格。建议在建模前写下一张“目标—位置—工况”对照表,作为后续网格与边界处理的依据。
二、单元类型选择:维度与阶次
单元选择的核心是几何维度与插值阶次两个维度。几何维度上:实体单元(四面体、六面体)适合三维应力状态;壳/板单元适合一个方向尺寸远小于另两个方向的结构;梁/杆单元适合长细比大的构件。插值阶次上:线性单元(如C3D4、C3D8)节点少、刚度偏硬,弯曲问题易出现剪切自锁;二次单元(如C3D10、C3D20)对弯曲和应力梯度适应性更好,但计算量更大。
面向相关性分析,有两条经验值得记住:第一,优先使用与试验物理量匹配的单元——测应变就用能输出该方向应变的单元与积分点;第二,避免同一区域混用不兼容单元,例如壳与实体连接处若不做过渡处理,会引入虚假局部刚度,使比对曲线在连接附近出现无法解释的偏差。
三、网格策略:收敛性与测点对齐
网格不是越密越好,而是要在关注区域达到收敛。推荐做法是:先做一次网格收敛检查,将关注点响应(如最大挠度、某阶频率)随网格尺寸加密的变化画出来,当相邻两次加密的相对变化小于约2%时,可认为该区域已收敛。对于应力集中区,还需检查积分点应力是否稳定。
另一个常被忽视的要点是测点对齐。试验测点往往不在有限元节点上,若直接取最近节点值,会引入插值误差。正确做法是在模型上按试验坐标建立参考点,通过节点插值或建立刚性/耦合约束提取对应位置的响应。下表给出一个示例性的网格与测点处理对照(数值为教学假设):
| 关注量 | 建议单元 | 关注区网格尺寸 | 测点提取方式 |
|---|---|---|---|
| 整体挠度 | 二次壳/梁 | 约板厚的2~3倍 | 节点插值到试验坐标 |
| 局部应变 | 二次实体 | 约应变片长度的1/3 | 积分点外推至表面 |
| 低阶频率 | 二次实体/壳 | 全局均匀,圆角处加密 | 模态节点与测点匹配 |
四、边界条件与载荷施加原则
边界条件是相关性误差最常见的来源。试验中的夹具并非绝对刚性,而有限元中常简化为固支,这会高估约束刚度。面向相关性分析,应遵循三条原则:
- 可解释原则:每个约束都要能对应到试验中的真实支撑或夹具,不能为了“让模型能算”而随意添加约束。
- 最小约束原则:只约束消除刚体位移所必需的自由度,其余自由,避免过约束掩盖真实变形。
- 对称性一致原则:若试验利用对称面,模型也应在对称面施加对称边界;若试验未利用,模型不应擅自取半。
载荷施加同样要与试验一致。集中力应施加在真实加载点,并通过耦合或分布耦合避免单节点应力奇异;分布载荷(如压力、惯性载荷)应确认作用面与方向。对于螺栓预紧、接触等非线性因素,若试验中确实存在,应在模型中保留,否则相关性偏差会被误判为模型误差。
五、建模要点与相关性误差的对应关系
下面这张图把本节四条主线与它们最可能引起的相关性偏差对应起来,便于在比对结果异常时快速定位建模环节。
当比对曲线出现系统性偏差时,可按上图从“单元—网格—边界—载荷”四条线逐一排查:若整体频率普遍偏高,优先怀疑边界过约束;若局部应变峰值对网格敏感,优先检查网格收敛与积分点提取;若弯曲变形整体偏小,优先检查单元阶次与剪切自锁。把建模环节与偏差特征建立对应关系,才能让相关性分析从“看结果”变成“找原因”。