1. 复合材料仿真入门与ANSYS工具链

复合材料结构仿真与各向同性金属仿真最大的差别,在于材料本身。金属通常只需弹性模量、泊松比和密度三个常数即可描述线弹性行为,而一层碳纤维增强环氧树脂单向带,其纵向模量 E₁ 可达 130 GPa 以上,横向模量 E₂ 往往只有 8~10 GPa,面内剪切模量 G₁₂ 约 4~5 GPa,泊松比 ν₁₂ 约 0.3。这种"一个方向很强、垂直方向很弱"的特性,就是各向异性。仿真中如果把它当成各向同性材料处理,层合板的刚度和屈曲载荷会被严重高估,铺层顺序的优化也失去意义。

1.1 从单层到层合板:仿真的对象是什么

复合材料结构仿真的基本对象是"单层"(lamina)与由单层按不同角度堆叠而成的"层合板"(laminate)。单层在材料坐标系 1-2-3 下是横观各向同性的:1 为纤维方向,2 为面内垂直纤维方向,3 为厚度方向。工程上常用 9 个工程常数描述正交各向异性单层:E₁、E₂、E₃、ν₁₂、ν₁₃、ν₂₃、G₁₂、G₁₃、G₂₃。对于常见的二维壳单元建模,通常简化为平面应力状态,只需 E₁、E₂、ν₁₂、G₁₂ 四个参数,加上单层厚度和铺层角即可。

层合板中每一层的纤维方向与结构全局坐标系 X 轴之间的夹角称为铺层角,记为 0°、±45°、90° 等。铺层顺序用类似 [0/45/-45/90]s 的记号表示,下标 s 表示对称铺层。仿真中真正参与刚度计算的,是各层经过坐标变换后叠加得到的 ABD 矩阵:A 为面内刚度,B 为耦合刚度,D 为弯曲刚度。B 矩阵是否为零,直接决定层合板会不会在面内载荷下产生弯扭耦合。这正是复合材料仿真必须逐层建模、而不能用一层等效材料草率替代的原因。

1.2 ANSYS 工具链的分工

ANSYS 对复合材料仿真的支持由几个模块协同完成。Workbench 是流程集成平台,负责几何、材料、网格、求解和结果的统一管理;Engineering Data 用于录入各向异性工程常数;ACP(ANSYS Composite PrepPost)是专门的复合材料前后处理模块,负责铺层定义、铺层可视化、失效准则计算和结果提取;Mechanical 或 APDL 求解器负责实际有限元计算。典型分工是:几何在 SpaceClaim 或 Discovery 中准备,材料在 Engineering Data 中定义,网格在 Mechanical 中划分,铺层在 ACP 中定义并映射到网格上,求解后回到 ACP 或 Mechanical 查看层间应力与失效指标。

需要特别注意的是,ACP 中的铺层是定义在"壳网格"或"实体壳"上的,铺层信息通过坐标系和铺层角映射到单元。因此几何面的法向、参考方向、铺层角基准必须与设计意图一致,否则 0° 方向可能整体偏转,导致结果完全错误。

复合材料仿真工具链与数据流 几何准备 SpaceClaim / Discovery 材料定义 Engineering Data 网格划分 Mechanical ACP 铺层定义 铺层角 / 顺序 / 厚度 坐标系与参考方向 铺层可视化与检查 失效准则设置 求解器 Mechanical / APDL 层间应力计算 ABD 刚度组装 非线性 / 接触 结果判读(ACP Post) 层间剪应力 τ₁₃ / τ₂₃ · 纤维拉伸/压缩失效指数 · 基体失效指数 最大失效值 < 1 为安全 · 首层失效位置 · 铺层顺序优化依据 Tsai-Wu / Hashin / 最大应力准则 铺层迭代优化

1.3 典型仿真流程与文件管理

一个完整的复合材料仿真流程通常包含以下环节:几何清理与中面抽取、材料工程常数录入、壳网格划分、ACP 中建立铺层序列、定义参考坐标系与铺层角、设置失效准则、提交求解、在 ACP Post 中查看逐层结果。以一块 300 mm × 200 mm 的层合板为例,铺层设为 [0/45/-45/90]s,单层厚度 0.125 mm,共 8 层,总厚度 1 mm。材料采用 T300/环氧单向带,示例参数取 E₁ = 135 GPa、E₂ = 9.5 GPa、ν₁₂ = 0.3、G₁₂ = 4.8 GPa。一端固支,另一端施加 1000 N 面内拉伸载荷。求解后查看各层纤维方向应力 σ₁,若最大失效指数出现在 90° 层的基体拉伸模式,说明该层是薄弱环节,可考虑调整铺层顺序或增加 ±45° 层比例。

文件管理方面,Workbench 项目目录下会生成 .wbpj 工程文件、.mechdat 网格与模型数据、ACP 的 .acp 铺层文件以及求解结果 .rst 文件。建议按"项目名/几何/材料/网格/ACP/结果"分层存放,并在 ACP 中为每个铺层组命名,便于后续修改和追溯。铺层定义一旦改动,必须重新求解,不能只更新结果视图。

初学者最容易犯的错误是:在 ACP 中定义了铺层,却忘记在 Mechanical 中把壳单元的截面类型设为"复合材料"或未正确关联 ACP 设置,导致求解器仍按各向同性材料计算。提交前务必检查材料方向与铺层映射是否生效。

1.4 学习路径建议

建议按"材料常数 → 单层刚度 → 层合板 ABD → 铺层定义 → 失效判据 → 结果解读"的顺序推进。先用手算验证一个 [0/90]s 层合板的等效模量,再在 ANSYS 中复现,确认软件设置与理论一致后,再进入复杂几何与优化。这样能避免"软件会点、结果看不懂"的常见困境。