1. 复合材料损伤建模工程导论

复合材料层合板由纤维与基体两种材料在细观尺度上组合而成，宏观上表现为明显的各向异性。以常见的碳纤维/环氧单向带为例，沿纤维方向的模量 E₁ 约为 130～150 GPa，而垂直纤维方向的模量 E₂ 通常只有 8～12 GPa，两者相差十倍以上；面内剪切模量 G₁₂ 约为 4～6 GPa，泊松比 ν₁₂ 约 0.25～0.35。这种"方向决定刚度"的特性意味着：同一块层合板，铺层角度不同，其承载能力、失效模式与损伤演化路径都会显著改变。损伤建模工程的核心任务，就是把这种各向异性与损伤的耦合关系，转化为有限元可以稳定求解、工程师可以判读的数值模型。

各向异性与损伤的耦合特征

复合材料的损伤不是单一裂纹的扩展，而是多种机制并存并相互耦合：基体开裂、纤维-基体界面脱粘、层间分层、纤维断裂。它们在不同应力状态下被激活，并随载荷增加彼此触发。例如横向拉伸时，基体开裂先出现，裂纹沿纤维方向扩展；当裂纹密度饱和后，局部应力集中会诱发分层，分层又改变层间载荷传递，反过来加剧相邻层的基体开裂。这种"损伤引发损伤"的链式关系，使复合材料失效呈现明显的渐进性和非线性。

更关键的是，损伤会改变材料的各向异性程度。初始时材料在纤维方向很强、横向很弱；基体开裂后，横向刚度进一步退化，而纤维方向刚度几乎不变。因此，损伤变量不能用一个标量描述，而必须按方向、按机制分别定义。工程上常用一组内部状态变量来表征：纤维拉伸损伤 d_ft、纤维压缩损伤 d_fc、基体拉伸损伤 d_mt、基体压缩损伤 d_mc。每个变量在 0～1 之间演化，0 表示完好，1 表示完全失效。刚度退化通常写成如下形式（以面内为例）：

E1_d = (1 - d_ft) * E1
E2_d = (1 - d_mt) * E2
G12_d = (1 - d_mt) * (1 - d_ft) * G12

这里 d_ft 与 d_mt 的乘积项体现了纤维与基体损伤对剪切刚度的共同削弱。需要强调的是，上述公式是工程简化表达，实际本构中还需考虑损伤变量的演化阈值、正则化与能量耗散，避免网格依赖和数值不稳定。

有限元建模在工程中的作用

有限元在复合材料损伤分析中承担三类角色。第一是应力状态预测：通过层合板壳单元或实体单元，计算各铺层在给定载荷下的面内应力与层间应力，为损伤判据提供输入。第二是损伤演化模拟：在单元积分点上引入损伤变量，按判据判断是否起裂，再按演化律更新刚度，实现渐进失效分析。第三是承载能力与失效模式评估：通过载荷-位移曲线的拐点、损伤分布云图和最终失效路径，判断结构是纤维主导失效还是基体/分层主导失效。

与金属结构不同，复合材料建模对网格、铺层定义和边界条件更敏感。层间应力对网格尺寸高度依赖，分层分析往往需要内聚力单元或虚拟裂纹闭合技术；损伤演化又要求特征长度正则化，否则结果会随网格加密而失真。因此，建模工程不是"画网格、加约束、点求解"，而是要在建模前明确分析目标、损伤机制和可接受的简化程度。

典型损伤问题与建模目标

工程中常见的复合材料损伤问题可归纳为四类：开孔层合板的挤压与拉伸失效、低速冲击后的分层与基体损伤、螺栓连接的孔边应力集中与拉脱、以及加筋壁板的屈曲后损伤扩展。不同问题对应不同的建模目标。开孔拉伸关注孔边纤维断裂与基体开裂的竞争；低速冲击关注层间分层面积与冲击后剩余压缩强度；螺栓连接关注孔边挤压损伤与层间剥离；加筋壁板则关注屈曲波形与筋条脱粘的耦合。

建模目标一旦确定，模型粒度也随之确定。若只关心整体刚度退化，可采用壳单元加渐进损伤本构；若关心分层位置与面积，则必须在层间布置内聚力单元；若关心冲击凹坑与纤维断裂，则需要实体单元加接触与失效删除。目标不同，计算成本可能相差一个数量级，因此"先定目标、再定模型"是复合材料损伤建模的基本工程原则。

课程项目对象与交付物定义

本课程以一个贯穿始终的项目对象为载体：含中心开孔的碳纤维/环氧层合板，铺层顺序为 [0/±45/90]ₛ，共 8 层，单层厚度 0.125 mm，总厚度 1 mm，开孔直径 10 mm，板尺寸 150 mm × 100 mm。材料参数采用上文给出的典型值，作为教学示例参数，实际工程中应以材料数据表或试验标定为准。加载方式为两端施加位移载荷，边界条件为一端固支、一端沿纤维方向拉伸。

项目交付物包括四项。第一，有限元模型文件：包含几何、网格、铺层定义、材料本构与边界条件，网格在孔边加密，层间布置内聚力单元。第二，损伤演化结果：输出载荷-位移曲线、各损伤变量的云图与孔边损伤起始载荷。第三，结果判读报告：说明失效模式是纤维主导还是基体主导，给出损伤起始与最终失效的载荷区间，并讨论网格敏感性。第四，建模决策记录：记录单元类型选择、特征长度取值、损伤演化参数与正则化方式，便于复现与审查。