1. 复合材料结构分析概述与Abaqus工作流
复合材料结构分析的核心难点,并不在于“会不会点按钮”,而在于是否理解材料本身的各向异性、层合板的铺层构造,以及这些构造如何决定分析类型与建模路径。Abaqus 提供了从材料定义、截面指派、铺层建模到求解与后处理的完整链路,本章先建立整体认识,再进入具体操作。
一、复合材料为什么不能按各向同性材料处理
金属材料通常只需弹性模量 E 与泊松比 ν 即可描述面内刚度,而单层复合材料(lamina)在纤维方向与垂直纤维方向上的刚度差异可达一个数量级。以典型碳纤维/环氧单层板为例,其工程常数大致为:E₁ = 135 GPa,E₂ = 9 GPa,G₁₂ = 5 GPa,ν₁₂ = 0.30。可以看到 E₁ 约为 E₂ 的 15 倍,这种强烈的方向依赖性意味着:如果按各向同性材料建模,层合板的弯曲刚度、屈曲载荷与失效位置都会出现明显偏差。
因此,复合材料单层在 Abaqus 中通常采用 工程常数(Engineering Constants) 定义正交各向异性,或采用 Lamina 类型直接输入 E₁、E₂、ν₁₂、G₁₂、G₁₃、G₂₃。前者适合三维实体单元,后者常用于壳单元与连续壳单元。需要强调的是,1、2、3 方向是材料坐标系:1 为纤维方向,2 为面内垂直纤维方向,3 为厚度方向。后续所有铺层角度、应力输出与失效判据都基于这一约定,全文与图示必须保持一致。
二、层合板结构与铺层定义
层合板(laminate)由多个单层按指定角度堆叠而成。铺层顺序通常写作 [0/±45/90]s 这类形式,其中 s 表示对称铺层。铺层角度是相对参考轴(一般为全局 1 轴)旋转的角度,正负号遵循右手定则。层合板的关键参数包括:单层厚度、铺层角、铺层顺序、总厚度与对称性。
在 Abaqus/CAE 中,层合板通过 Composite Layup 定义,可对壳单元或连续壳单元逐层指定材料、厚度与角度。对于壳单元,Abaqus 还提供 Section Points 用于输出沿厚度方向的应力分布;对于连续壳单元,则通过厚度方向积分点直接获得三维应力状态。理解这一点,对后续判断分层、纤维拉伸/压缩失效至关重要。
三、分析类型选择
复合材料结构分析常见类型包括:静力分析、屈曲分析、模态分析、瞬态动力学分析与显式动力学分析。选择依据主要看载荷性质、是否涉及接触与失效、以及是否需要捕捉分层扩展。
- 静力分析(Static, General):适用于大多数层合板强度与刚度校核,可结合 Hashin、Tsai-Wu 等失效准则输出损伤起始。
- 线性摄动屈曲(Buckle):用于薄壁层合板受压屈曲载荷预测,需先施加参考载荷再提取特征值。
- 频率分析(Frequency):获取固有频率与振型,判断结构是否与激励频率耦合。
- 显式动力学(Explicit):适合冲击、低速撞击与分层扩展,常配合内聚力单元(Cohesive)模拟层间开裂。
对于层间分层问题,Abaqus 提供 Cohesive Surface 与 Cohesive Element 两种途径。前者基于接触面定义,适合壳单元层合板;后者需在铺层之间插入零厚度或薄层单元,适合连续壳或实体单元。选择时需权衡建模工作量与计算成本。
四、Abaqus/CAE 复合材料建模模块总览
Abaqus/CAE 中与复合材料直接相关的模块包括:Property 模块中的 Composite Layup、Section 指派、Material 中的 Lamina/Engineering Constants,以及 Mesh 模块中的壳单元与连续壳单元设置。典型工作流如下:
- 在 Part 模块创建几何,壳结构可直接建面,实体结构建三维体。
- 在 Property 模块定义材料:选择 Lamina 或 Engineering Constants,输入 E₁、E₂、ν₁₂、G₁₂ 等参数。
- 创建 Composite Layup,逐层指定材料、厚度与角度,并指派到对应区域。
- 在 Mesh 模块划分网格,壳单元选用 S4R 或 SC8R,注意厚度方向积分点数量。
- 在 Step 模块选择分析类型,设置输出变量如 S、E、DAMAGEFT 等。
- 提交 Job 求解,在 Visualization 模块查看应力云图、失效指数与变形。
五、建模前的关键检查项
在进入具体分析前,建议确认以下内容:材料方向是否与几何参考轴一致;铺层角度正负号是否按右手定则定义;壳单元厚度方向积分点是否足够捕捉弯曲应力;输出请求是否包含层间应力与失效判据;边界条件是否与真实约束等效。这些检查项直接决定后续结果是否可信。