面向混合模态分层失效的复合材料梁离轴铺层角度稳健优选策略

《Journal of Materials Research and Technology》：Robust Selection of Off-Axis Ply Orientations for Enhancing Structural Load Capacity Against Mixed-Mode Delamination in Composite Beams

　　

分层（delamination）仍是碳纤维复合材料梁的主要失效限制，尤其是当既有界面缺陷处同时发生张开（opening）与滑移（sliding）时。本研究探讨离轴铺层角度（off-axis ply orientation）是否不仅能提升名义混合模态（mixed-mode）抗力，还能在材料、界面及制造参数变化时保持该提升效果。研究人员将AS4/PEEK混合模态弯曲（mixed-mode bending，MMB）试件的层合梁有限元模型与中面界面的双线性内聚区（cohesive-zone）描述相结合。混合模态损伤演化由Benzeggagh–Kenane关系控制，同时耦合COMSOL Multiphysics与MATLAB进行参数扫描、不确定性传播及稳健设计优选。数值模型复现了实验混合模态弯曲响应，峰值载荷误差为0.33%，初始刚度误差为2.98%，归一化均方根误差为2.26%。在初始角度研究中，将离轴角度从0°增大至60°使峰值载荷从233.0 N提升至276.8 N，并将最终分层扩展从34.8 mm减小至26.4 mm。这些变化的最大增量分别为18.8%和24.1%，对应于裂纹扩展的缩减。精细化搜索表明，名义最优值取决于模态混合比（mode mixity），从20% Mode-II贡献时的30°移至50%时的45°及80%时的60°。考虑变异性后，优选设计发生变化。尽管确定性分析倾向于45°构型，但47.5°层合板提供了259.6 N的第五百分位峰值载荷和30.0 mm的第九十五百分位分层扩展。结果表明，当偏离确定性最优值出现微小偏差时，响应变异性可在保持大部分名义承载能力的同时得以降低。所开发模型为混合模态断裂抗力及不确定性情形下的复合材料铺层选择提供了实用方法。

复合材料梁离轴铺层角度稳健优选：混合模态分层失效的结构抗力提升研究解读

一、研究背景与问题提出

碳纤维增强聚合物（carbon-fiber-reinforced polymer，CFRP）层合板因其低质量、高刚度与高强度特性，广泛用于翼梁、加筋条、桁条、肋、框架及载荷传递构件等梁类结构。层合复合材料可通过调整单层（ply）的取向与位置来定制结构响应，这是其相较于传统各向同性梁的主要优势。然而，铺层取向与厚度方向刚度的变化会产生层间应力场，使原本按全局刚度或面内强度准则表现优异的层合板，其界面可能反而易于分离。分层（delamination）是复合材料层合板中最难控制的失效模式之一，可在层合板内部发生而不在外部显现。层间分离可源于制造缺陷、富树脂区、孔隙、铺层波纹度、自由边区域、冲击损伤及几何不连续等。当界面开始张开或滑移时，弯曲与剪切刚度下降、内部应力场重新分布，进而可能触发局部屈曲、基体开裂与纤维断裂等后续失效机制。由于面内载荷主要由纤维承担，而平行于纤维方向的分离抗力主要取决于聚合物基体及纤维-基体界面，分层问题因而更为严峻。

复合材料梁中的裂纹尖端状态通常并非纯Mode-I或纯Mode-II。分离臂的面内弯曲产生类似铰链的张开元素，而横向剪切与差异纵向位移导致滑移。这些相对贡献随裂纹位置、子层合板刚度、载荷施加方式及试件变形而变化，因此混合模态断裂抗力不能通过独立考虑Mode-I与Mode-II来确定，需要采用相互作用关系来定义张开与滑移联合作用下的临界能量。Benzeggagh–Kenane（BK）准则是该领域最常采用的关系之一，已广泛用于内聚区与断裂力学分析。

混合模态弯曲（mixed-mode bending，MMB）试验提供了可控组合两种断裂模态的方式，但其测量响应不仅取决于名义杠杆比。初始裂纹长度、试件尺寸、加载位置、界面取向、子层合板刚度失配、裂纹前缘曲率及断裂过程区尺寸均可能影响表观抗力。在多向层合板中，裂纹附近铺层的纤维取向与各向异性刚度不同，裂纹扭转、非均匀桥接、迁移至相邻界面及层内开裂等行为可能伴随名义界面扩展。因此，为单向0°/0°界面标定的断裂性能不能不加区分地移植至θ/?θ、θ/90°等界面。铺层顺序在结构与局部尺度上均影响分层行为：通过层合板的A、B、D矩阵决定拉伸刚度、弯曲刚度、耦合行为及开裂臂的变形；在界面处，相邻铺层取向控制刚度失配并影响裂纹前方传递的法向与切向牵引力。具有相似全局弯曲刚度的层合板可能发展出不同的裂纹尖端场与不同的内聚损伤速率。

现有铺层优化研究大多优化结构全局量或面内失效量，分层通常被简化为约束或假定缺陷对刚度的影响，而非作为渐进混合模态萌生与扩展过程。同时，制造中的铺层角度误差、厚度变化、间隙、重叠、波纹度、孔隙率、固化条件及纤维体积分数变化均会影响制造层合板的响应；内聚强度与断裂能也存在实验离散性，载荷位置与模态混合比难以精确已知。因此，可靠性与稳健性方法在复合材料设计中日益受到重视。然而，现有可靠性优化研究虽纳入了预定义分层，但未将渐进混合模态裂纹作为核心优化响应。本研究面临的剩余问题因此不仅是识别单次模拟中产生最高载荷的角度——名义最大值可能出现在设计空间的狭窄区域内，且对铺层厚度、界面强度或断裂能的微小变化敏感；附近具有略低名义最大值的设计可能因提供更低响应方差与更保守载荷而更优。本研究通过将渐进混合模态内聚分析与确定性和不确定性感知的铺层优选相结合来解决该问题。

二、研究内容与主要结论

本研究针对基于[0/θ/?θ/90]_3s重复块的平衡对称层合板族，在广泛的离轴角度、模态比与初始裂纹长度范围内进行研究，随后在有利区域附近进行精细化角度搜索，最后引入材料、界面与制造变异性，识别确定性最优与稳健设计之间的差异。研究建立了AS4/PEEK梁的渐进混合模态分层分层模型并进行了验证；系统评估了离轴铺层角度、局部模态混合比与初始分层长度的影响；通过MATLAB–COMSOL工作流自动耦合非线性分析；比较了确定性与不确定性感知的界面角度建议；并从承载能力、裂纹扩展抑制、峰后响应与稳健性角度解释了所得设计。需要强调的是，该优化并非旨在改变材料的固有界面断裂韧性，而是通过调整子层合板的弯曲与扭转刚度（D_ij）来最小化到达裂纹前缘的裂纹驱动力，从而实现结构级分层抗力提升。

三、关键技术方法

研究人员采用的主要技术方法包括：其一，基于COMSOL Multiphysics的三维层合壳（layered-shell）有限元模型，通过Layered Material定义逐层赋予各单层厚度、取向与正交各向异性属性，中面界面采用双线性内聚区模型（cohesive zone model，CZM）描述，混合模态损伤演化由Benzeggagh–Kenane（BK）关系控制；其二，通过LiveLink for MATLAB实现COMSOL–MATLAB耦合工作流，由MATLAB控制参数研究、组织结果并计算不确定性，自动生成参数矩阵或拉丁超立方（Latin Hypercube Sampling，LHS）样本，调用非线性求解并提取响应数据；其三，采用位移控制的非线性静力分析并辅以延续扫描以追踪峰后软化分支；其四，不确定性传播采用n=50的LHS抽样，对临界断裂能（G_Ic±10%）、界面强度（±10%）、层合板纵向模量（E₁₁±5%）及铺层厚度（t_ply±3%）施加均匀分布容差，采用矩估计法计算第五百分位峰值载荷与第九十五百分位裂纹扩展作为稳健性指标。模型验证基于Reeder和Crews的实验基准数据（AS4/PEEK，η=0.50混合模态条件），材料参数直接取自文献独立实验数据，未进行曲线拟合。

四、研究结果

4.1 网格无关性与模型验证

四种沿分层路径的网格（25、50、100、200单元）研究表明，最大载荷从最粗网格的184.0 N增至极细网格的201.6 N，相对200单元模型的最大载荷误差分别为8.73%、3.27%与0.79%。100单元网格被选用于参数研究与优化，其在弹性斜率、最大载荷、峰后响应与最终裂纹扩展方面均准确捕捉了整体行为。内聚过程区长度（l_cz）估计为Mode-I约1.6 mm至混合模态约2.4 mm，100单元界面离散产生Δx≈0.35 mm的局部单元尺寸，每个活跃内聚过程区包含4–7个单元，满足裂纹扩展稳定性条件。模型验证结果显示，数值模型预测峰值载荷274.5 N与实验值275.4 N的相对误差为0.33%；初始刚度数值51.14 N/mm与实验52.71 N/mm相差2.98%；全载荷-位移响应的归一化均方根误差为2.26%。在η=0.20的较低模态比下额外基准验证的预测误差为2.11%。与辅助三维实体单元基准的数值比较显示全局载荷-位移行为与分层扩展速率的差异小于2.5%，验证了层合壳表示的适用性。

4.2 铺层角度对混合模态分层抗力的影响

初始角度扫描考虑五个离轴角度（0°、15°、30°、45°、60°）。将角度从0°增至60°使初始刚度从46.8 N/mm提升至52.4 N/mm；损伤萌生载荷从179.4 N增至221.4 N，对应内聚退化起始延迟23.4%。最大载荷分别为0°时233.0 N、15°时252.5 N、30°时276.8 N、45°时268.4 N、60°时244.1 N，最大值出现在30°，相对0°层合板提升18.8%。最终分层扩展从0°的34.8 mm降至30°的26.4 mm，缩减24.1%。局部Mode-II贡献随裂纹扩展的变化显示，30°层合板的局部模态混合比从0.4625逐渐增至0.5175，45°层合板从0.4940缓慢增至0.5060，60°层合板从0.5345逐渐降至0.4895，反映了离轴弯曲/扭转耦合引起的模态漂移。

4.3 模态混合比与初始裂纹长度的影响

精细化搜索表明，名义最优角度取决于模态混合比：在20% Mode-II贡献时最优为30°，在50%时为45°，在80%时为60°。初始裂纹长度通过归一化裂纹长度（a₀/2L）变化，影响裂纹尖端约束与模态分配，进而影响最优角度选择。

4.4 不确定性与稳健设计优选

考虑变异性后，优选设计发生变化。尽管确定性分析倾向于45°构型，但47.5°层合板提供了259.6 N的第五百分位峰值载荷和30.0 mm的第九十五百分位分层扩展。结果表明，当偏离确定性最优值出现微小偏差（2.5°）时，响应变异性可在保持大部分名义承载能力的同时得以降低。参数敏感性研究验证了绝对阈值载荷随界面韧性缩放，但优化趋势、局部模态相互作用及稳健构型（θ=47.5°）的优越性得以保持。

五、讨论与结论

研究讨论指出，离轴界面（θ/?θ）可能表现出纤维桥接、基体屈服等微观机制，理论上需要取向相关的内聚属性；本研究假设恒定基线界面并通过±10%扰动进行敏感性评估。AS4/PEEK体系的热塑性基体屈服与纤维桥接是额外耗散机制，实际中对峰后平台载荷与总断裂抗力（G_R）贡献更显著，这些微观机制所依赖的固有树脂性能与局部纤维-基体脱粘对所有铺层共享，因此是总断裂抗力的共同缩放因子。绝对峰后韧性被适度低估，但离轴界面的结构弯曲-延伸-扭转耦合与混合模态分配对于识别稳健角度（θ=47.5°）仍然有效。研究结论为：渐进混合模态分层抗力可通过调整子层合板刚度实现结构级提升；名义最优设计在参数变异下可能不再最优，考虑不确定性的稳健设计（47.5°）能在保持大部分名义承载能力的同时降低响应变异性；所开发的COMSOL–MATLAB耦合框架为混合模态断裂抗力与不确定性情形下的复合材料铺层选择提供了实用方法。