2. 复合材料铺层基础与HyperMesh表达

复合材料结构分析与金属结构最大的不同,在于“方向”本身就是结构的一部分。同一块板,纤维铺在哪、朝哪个方向、按什么顺序叠起来,都会直接改变刚度、强度和失效位置。因此在 HyperMesh 里做复合材料前处理,第一步不是画网格,而是把“单层—层合板—铺层顺序—坐标系”这条表达链建立正确。只要这条链中任意一环定义错了,后面的求解结果即使算得再顺,也没有工程意义。

单层:复合材料的最小表达单元

单层(lamina/ply)是复合材料铺层的最小单元,通常由纤维和基体组成,在有限元中一般被简化为正交各向异性材料。它需要两类信息:一是材料属性,二是几何属性。材料属性包括 E1、E2、G12、ν12 等工程常数,其中 1 方向沿纤维方向,2 方向垂直于纤维方向;几何属性则是单层厚度 t,常见预浸料单层厚度在 0.1~0.25 mm 之间,例如 0.125 mm、0.2 mm 都是教学中常用的示例值。

在 HyperMesh 中,单层并不直接等同于一个实体单元,而是通过属性(property)和铺层定义来表达。典型做法是:先建立二维壳网格,再给壳单元赋予 PSHELL 或 PCOMP 类属性,其中 PCOMP 用于定义层合板,每一层通过 MID(材料号)、T(厚度)、THETA(铺层角)来描述。也就是说,单层在 HyperMesh 里是“材料 + 厚度 + 角度”的组合,而不是单独的一个几何面。

层合板:单层按顺序叠加的整体

层合板(laminate)由多层单层按一定顺序和角度叠合而成。它的力学行为不是各层简单相加,而是通过经典层合板理论中的 A、B、D 矩阵耦合起来。A 是面内刚度,D 是弯曲刚度,B 是拉弯耦合刚度。铺层顺序对称时 B 矩阵接近零,可以显著减少固化后的翘曲;铺层顺序不对称时,即使总厚度相同,弯曲和面内变形也会明显不同。

在 HyperMesh 中表达层合板,核心是铺层表(laminate layup)。例如一个四层板可以写成 [0/90]s,展开为 0°、90°、90°、0°,共 4 层;若单层厚度 0.2 mm,则总厚度 0.8 mm。再如 [0/±45/90]s,展开为 0°、+45°、−45°、90°、90°、−45°、+45°、0°,共 8 层,总厚度 1.6 mm。这些角度和顺序必须与设计图一致,不能随意调换。

纤维方向与铺层角

纤维方向是单层刚度的主方向,通常记为 1 方向。铺层角 θ 是纤维方向与参考轴之间的夹角,逆时针为正。工程上常以结构的主受力方向或全局 X 轴作为参考轴。例如 0° 表示纤维沿参考轴,90° 表示垂直,+45° 和 −45° 表示斜向。铺层角定义必须与材料坐标系、单元坐标系一起看,否则同一个 45° 在不同单元上可能指向完全不同的方向。

材料坐标系与单元坐标系的关系

这是 HyperMesh 复合材料前处理中最容易出错的地方。材料坐标系(material coordinate system)定义单层材料的主方向,即 1-2 方向;单元坐标系(element coordinate system)是壳单元自身的局部坐标,用于定义法向和面内方向。铺层角 θ 通常是在单元坐标系面内,从单元参考轴旋转到材料 1 方向的角度。因此,同一个 θ 值,如果单元坐标系不同,实际纤维方向也会不同。

在 HyperMesh 中,壳单元默认的单元坐标系与网格生成方向、节点编号顺序有关。对于规则平板,若网格沿全局 X 方向排列,单元 x 轴通常与全局 X 一致;但对于曲面、过渡区或映射网格,单元坐标系可能发生偏转。此时必须通过材料方向(material orientation)或铺层方向(ply orientation)显式指定参考方向,例如使用 vector、curve 或 coordinate system 来定义 0° 基准。

常见错误:只设置了铺层角,却没有检查单元坐标系。结果在平板上看似正确,到了曲面或圆角区域,纤维方向发生偏转,导致局部刚度异常。排查时应打开单元坐标系显示,并逐区检查 0° 参考方向是否与设计一致。

HyperMesh 中的典型表达流程

以一个 300 mm × 200 mm 的平板为例,铺层为 [0/45/−45/90]s,单层厚度 0.125 mm,材料为碳纤维/环氧,E1=135 GPa,E2=10 GPa,G12=5 GPa,ν12=0.3。操作上可以按以下顺序进行:

  1. 创建材料:在 HyperMesh 的 Material 面板中新建 MAT8,输入 E1、E2、G12、ν12 和密度等参数。
  2. 创建铺层属性:新建 PCOMPP 或 PCOMP,按铺层表逐层输入 MID、T、THETA,注意对称铺层可用 SYM 选项简化。
  3. 划分壳网格:用 2D 壳单元划分平板,控制单元尺寸约 5 mm,保证每层方向变化区域有足够单元。
  4. 指定材料方向:在需要时创建局部坐标系或方向矢量,作为 0° 参考轴,并赋给对应单元或属性。
  5. 检查单元坐标系:打开 element system 显示,确认各区域单元 x 轴方向与预期一致。
  6. 输出求解文件:导出为 OptiStruct 或 Nastran 格式,检查 PCOMP 卡片中的层顺序、角度和厚度是否正确。

完成上述设置后,可以通过查看导出的 PCOMP 卡片来验证。例如铺层 [0/45/−45/90]s 应展开为 0、45、−45、90、90、−45、45、0 共 8 层,每层厚度 0.125 mm,总厚度 1.0 mm。若卡片中层数、角度或厚度与设计不符,说明铺层表或对称设置有问题,应回到 HyperMesh 中修正,而不是在求解器里强行改卡片。

单层 → 层合板 → 铺层顺序 → 坐标系关系 单层 Lamina 材料:E1, E2, G12, ν12 厚度:t = 0.125 mm 方向:1 沿纤维,2 垂直 HyperMesh:MAT8 + 厚度 层合板 Laminate 铺层表:[0/45/−45/90]s 展开:8 层,总厚 1.0 mm A / B / D 刚度矩阵 HyperMesh:PCOMP 属性 铺层顺序 0° → 45° → −45° → 90° 90° → −45° → 45° → 0° 对称 → B ≈ 0 顺序影响翘曲与耦合 材料坐标系与单元坐标系的关系 单元坐标系 (Element CS) x_e y_e 由网格生成方向决定,可能偏转 材料坐标系 (Material CS) 1 (纤维) 2 由铺层角 θ 和参考方向定义 θ 旋转 逆时针为正 检查要点:单元坐标系显示 → 0° 参考方向 → 铺层角 θ → PCOMP 卡片层序与厚度

从上面的关系可以看出,HyperMesh 中的复合材料表达不是孤立的材料卡片,而是一条从单层到层合板、从铺层角到坐标系的完整链条。实际建模时,建议每完成一个区域就检查一次单元坐标系和铺层方向,尤其是曲面过渡区、开孔区和加强区。只有把“方向”这件事在网格层面表达清楚,后续的刚度计算、失效判断和优化才有可靠基础。