复合材料 — CAE术语解释

分类: 术语集 | 2026-01-15

CAE visualization for composite material - technical simulation diagram

复合材料

🧑‍🎓

老师，复合材料不仅仅是CFRP吧？从CAE分析的角度来看，应该如何分类？

复合材料的理论基础

复合材料的基本概念和支配方程

🧑‍🎓

复合材料的"异向性"具体是指什么特性的差异？与金属等向性材料的根本区别是什么？

🎓

根本区别在于材料常数的数量。等向性材料只需要杨氏模量和泊松比这2个参数就能描述弹性，而具有一般正交异向性的复合材料层板需要9个独立的弹性常数。例如，碳纤维增强塑料（CFRP）单向材在纤维方向的杨氏模量

$$ E_1 $$

约为140GPa，垂直方向

$$ E_2 $$

约为10GPa，相差10倍以上。这种方向依赖性既是设计的自由度来源，也是分析的难点所在。

🧑‍🎓

9个弹性常数具体是什么？它们如何通过实验确定？

🎓

正交异向性的本构关系表示为

$$ \{\sigma\} = [C]\{\epsilon\} $$

，刚度矩阵[C]的独立分量为

$$ E_1, E_2, E_3, \nu_{12}, \nu_{13}, \nu_{23}, G_{12}, G_{13}, G_{23} $$

这9个。这些参数通过ASTM D3039（拉伸试验）和ASTM D3518（面内剪切试验）等标准规程的试验来确定。在实际工程中，这些值通常可以在制造商目录中找到，例如Torayca T800S的数据表。

🧑‍🎓

单层和层板的区别是什么？在CAE中通常处理哪一个？

🎓

单层是指具有单一纤维方向的一层（相当于一张预浸料），而层板是指这些层的积层结构。在CAE中，首先定义各单层的性质，然后输入积层顺序（例如[0/90/±45]s）和各层的厚度，再用古典积层理论（CLT）计算层板的等效刚性。在实际设计中，如飞机机翼蒙皮等常由数十层组成的层板构成。

🧑‍🎓

古典积层理论（CLT）是如何计算层板刚性的？不是简单平均吧？

🎓

完全正确，这涉及厚度方向的积分。各单层的刚度矩阵[Q]经过其纤维方向角

$$ \theta $$

的坐标变换得到

$$ [\bar{Q}] $$

，然后用它来计算层板的面内刚度[A]、耦合刚度[B]和弯曲刚度[D]。用公式表示为

$$ A_{ij} = \sum_{k=1}^{N} (\bar{Q}_{ij})_k (z_k - z_{k-1}) $$

，其中

$$ z_k $$

是各层的位置。在非对称积层中，[B]不为零，导致拉伸和弯曲耦合，这是一个重要的特征。

复合材料的数值计算方法

FEM中的复合材料建模

🧑‍🎓

用FEM建模复合材料时，应该选用实心单元还是壳单元？判断标准是什么？

🎓

基本判断标准是厚度与跨度的比值。对于飞机或风电叶片等薄壳结构，厚度通常小于跨度的1/10，这种情况下应使用壳单元（如Abaqus的S4R、Ansys的Shell181）。而当需要评估粘结层或厚度方向的应力集中时，需要用连续体壳单元或实心单元（各层单独建模），但这会导致单元数量爆炸性增长，因此通常的做法是将分析范围局限于所关注的区域。

🧑‍🎓

用壳单元定义积层的"积层壳"功能是什么？需要输入什么数据？

🎓

积层壳中需要定义各层（单层）的材料、厚度、纤维方向角和积层顺序。Ansys Composite PrepPost（ACP）或Abaqus/CAE中的复合材料铺层功能提供表格形式的输入方式。具体来说，例如第1层：材料="T800S/Epoxy"、厚度=0.2mm、方向=0度，以此类推。求解器内部使用CLT计算[A]、[B]、[D]矩阵，并将其融入单元刚度。

🧑‍🎓

单元坐标系与材料坐标系的区别是什么？0度方向是相对于哪个坐标系的？

🎓

这是非常重要的细节。单元坐标系由单元的几何形状决定（例如壳单元的面内方向），是局部参考系。材料坐标系（或材料主方向）定义纤维的方向。纤维方向角（例如0度）定义为"从单元坐标系第1轴到材料坐标系第1轴（纤维方向）的旋转角"。在Ansys中以"SHELL181单元的X轴"为参考，在Abaqus中以"壳层平面的1方向"为参考。此设置错误会导致刚度应用在完全错误的方向上。

🧑‍🎓

划分网格时，应该以复合材料的哪个特征尺寸为基准？

🎓

对复合材料来说，面内的特征长度比层厚度（通常0.1～0.3mm）更重要。特别是在螺栓孔周围或形状急剧变化的位置，需要细化网格。一个参考指标是"剪切滞后长度"（shear lag length），其定义为

$$ \lambda^{-1} = \sqrt{\frac{G_{13}}{E_1} \cdot \frac{t}{h}} $$

，表示荷载有效传递的长度范围。其中

$$ t $$

是单层厚度，

$$ h $$

是总厚度。用这个长度的1/2至1/3作为网格尺寸，可以对层间应力的评估获得合理的精度。

复合材料的实务应用

分析工作流程和验证

🧑‍🎓

进行复合材料结构的线性静力分析时，具体的工作流程步骤是什么？

🎓

实务工作流程如下所示。

1. **几何模型简化**：仅提取肋板或加强筋作为面，创建中面模型。 2. **材料属性定义**：输入目录值（例如Hexcel IM7/8552的

$$ E_1=171 GPa, E_2=9.08 GPa, G_{12}=5.29 GPa, \nu_{12}=0.32 $$

）。同时准备强度值（

$$ X_t, X_c, Y_t, Y_c, S $$

）用于破坏判定。 3. **积层定义**：将设计图纸上的积层表（[45/0/-45/90]s等）转录到软件的积层定义工具中。 4. **网格生成**：曲率变化大的区域细化网格（单元尺寸≤5mm），平坦区域粗化（约20mm）。 5. **边界条件和荷载**：不过度简化实际支持条件（例如用面约束而非点支持）。 6. **求解**：用线性求解器执行。 7. **后处理**：评估位移、各层应力和破坏指标（详见后文）。

🧑‍🎓

为了信任分析结果，"验证"具体要做什么？仅看变形形状是不够的吧？

🎓

完全同意。至少应该进行以下3项验证。

1. **质量和体积检查**：软件计算的层板总质量是否与设计值（层数×层厚×密度×面积）在几个百分点以内一致？ 2. **反力检查**：所有反力的合力是否与施加的荷载平衡（通常误差<1e-3）？ 3. **简单情况验证**：例如，对单向材简支梁在自由端施加荷载，其端部挠度是否与梁理论

$$ \delta = PL^3/(3E_1I) $$

的计算值相符？其中

$$ I $$

是截面二次矩。这个验证可以确认材料常数和纤维方向角的设置是否正确。

🧑‍🎓

复合材料的破坏如何判定？不能像金属那样用单一等效应力来评估吧？

🎓

复合材料没有通用的破坏准则，需要选择不同的破坏模式。应用最广泛的是Tsai-Wu准则，它的形式为

$$ F_i \sigma_i + F_{ij} \sigma_i \sigma_j = 1 $$

，是张量形式的公式，可以考虑拉伸和压缩强度的差异。此外还有按模式判定破坏的Puck准则和Hashin准则。在实务中，使用Ansys ACP的"破坏准则"或Abaqus的"Hashin损伤"等软件内置功能，输出各层的破坏指标（0至1的值，1表示破坏），从而识别最危险的层及其破坏模式。

🧑‍🎓

层间脱层（分层）的分析需要特殊的技术吗？

🎓

是的，普通线性静力分析无法评估，需要采用基于断裂力学的方法或内聚区模型（CZM）。实用做法是在层间插入极薄的内聚单元，定义其界面强度和破坏能（例如Mode I的

$$ G_{Ic} = 0.2 N/mm $$

）。Abaqus使用内聚单元，Ansys使用INTER202单元等。这属于非线性分析，计算成本较高。在初期检查中，可以简单地核查层间剪切应力（

$$ \tau_{13}, \tau_{23} $$

）是否超过许可值（例如40MPa）。

复合材料软件对比

主要CAE软件功能对比

🧑‍🎓

用Ansys、Abaqus、MSC Nastran进行复合材料分析，各自的优势和特色功能是什么？

🎓

这是个很好的问题。软件选择往往取决于项目阶段。

- **Ansys（ACP + Mechanical）**：前处理工具"ACP"功能强大。可以直接从CAD定义积层区域，同时可视化纤维方向的定义。还支持制造性评估（悬垂仿真），在设计初期阶段经常使用。 - **Abaqus/CAE**：善于非线性分析（接触、大变形、破坏）与复合材料的耦合。通过VUMAT（用户定义材料子程序）可以在研发现场实现自定义的复合材破坏准则。 - **MSC Nastran**：在航空航天工业中具有历史标准地位。利用SOL 200可以进行复合材积层的优化（层厚和纤维方向角作为设计变量）。大规模线性分析（特征值、频率响应）计算速度快、稳定性好。

🧑‍🎓

"悬垂仿真"具体是计算什么？分析结果如何应用？

🎓

这是对预浸料铺贴到曲面金型时的仿真，会出现材料伸展、纤维方向改变、起皱等现象。Ansys ACP和PAM-FORM等专用软件进行计算。输出包括设计纤维方向的偏差（扭曲角）和纤维伸长率。这些与制造允许值（扭曲角通常在±5度以内）比较，判断设计是否可制造，或决定是否需要修改金型形状。这是连接设计和制造的关键工序。

🧑‍🎓

如果不用专业复合材料前后处理工具（如ACP），用通用前后处理器进行复合材分析，限制在哪里？

🎓

主要在"效率"和"表现力"两方面受限。用通用工具定义数十层的积层，必须逐一手动输入各单元的属性，实际上不可行。而且后处理中提取"特定层的纤维方向应力"等操作也很繁琐。进行多个设计方案的快速比较也很困难。专用工具可以统一管理积层表，自动将结果映射到各层，还能生成报告，大大提高工作效率。除了小规模检研外，专用工具的引进是必需的。

🧑‍🎓

开源FEM软件（CalculiX、Code_Aster）能进行复合材料分析吗？与商业软件有什么差距？

🎓

核心功能（异向性材料的线性静力分析）是可以的。CalculiX和Code_Aster都支持壳单元的积层定义。但差距明显，主要有3点。

1. **图形界面和前处理**：没有商业软件那样直观的积层定义界面，需要手工编写输入文件（.inp、.comm）或通过脚本生成。 2. **预验证的破坏准则**：Tsai-Wu等标准破坏准则不总是预装，用户需要自己实现（用户子程序）。 3. **技术支持和保证**：在涉及飞机等安全系数的领域，需要分析结果的保证和技术支持，商业软件是必须的。研究、教学或概念验证等级可以考虑开源软件。

复合材料故障排除

常见错误和解决方案

🧑‍🎓

复合材料分析中出现"刚度矩阵非正定"的错误，原因和解决方法是什么？

🎓

这是最常见的错误之一，通常源于材料常数输入错误或物理上不可能的值。具体原因包括：

1. **杨氏模量或剪切模量为负值**：简单的输入错误。 2. **泊松比不满足异向性稳定条件**：对正交异向性材料，应满足

$$ \nu_{12}^2 < E_1/E_2, \quad \nu_{13}^2 < E_1/E_3, \quad \nu_{23}^2 < E_2/E_3 $$

等条件。例如当

$$ E_1=140GPa, E_2=10GPa $$

时，

$$ \nu $$