对称内罚间断伽勒金（SIPG）法：面向复合材料层合板的高效三维高精度应力分析

《Composites Communications》：Efficient 3D analysis of composite laminates using high-aspect-ratio discontinuous Galerkin finite elements

本文针对复合材料层合板自由边分层失效预测中，传统共形有限元（Conforming FEM）计算成本高、高宽厚比模型难以划分高质量网格的问题，开展了基于对称内罚间断伽勒金（SIPG）法的三维分析研究。结果表明，该方法在高长宽比六面体网格上仍能保持最优L2收敛速率（位移p+1，应力p），并能以少一个数量级的自由度精确复现层间应力分布，为层合板高效仿真与优化提供了新途径。

　　

想象一下，你正在设计下一代航空发动机的叶片或是超跑的车身，这些结构大量使用了碳纤维复合材料层合板。这类材料强度高、重量轻，但有一个“阿喀琉斯之踵”——自由边分层。当相邻铺层的纤维方向不同时，材料属性的不匹配会在自由边附近产生复杂的层间应力（interlaminar stress），这往往是分层失效的元凶。传统的经典层合板理论（CLT）无法捕捉这种三维应力状态，工程师们通常求助于三维有限元分析（FEM）。

然而，这里又碰到了老大难问题：为了精确模拟沿厚度方向的应力梯度，往往需要划分非常细密的网格，导致计算自由度（DOFs）爆炸，耗时极长；而如果为了控制成本使用较粗的高长宽比（高Aspect Ratio）六面体网格，传统的共形有限元（Conforming FEM）又容易出现数值锁定或精度下降。如何在保证精度的前提下，大幅降低计算成本？这正是本研究要解决的核心痛点。

为了解决这一问题，研究人员提出了采用对称内罚间断伽勒金（Symmetric Interior Penalty Discontinuous Galerkin, SIPG）法，结合高长宽比六面体网格，来开展复合材料层合板的三维力学分析。该研究通过一系列严格的数值基准测试，验证了该方法的收敛性、稳定性以及对层间应力（尤其是自由边效应）的预测能力。研究结果表明，SIPG法即使在网格长宽比高达约38.6的情况下，依然能保持数值稳定，并以远少于传统共形有限元的自由度，准确计算出关键的层间应力分布。该论文发表在《Composites Communications》。

作者开展研究用到几个主要关键技术方法包括：1. 对称内罚间断伽勒金（SIPG）离散格式，用于线弹性问题的弱形式求解，通过惩罚项保证单元间连续性；2. 结构化高长宽比六面体网格生成技术，结合光滑映射函数聚类节点以解析陡峭应力梯度；3. 分层（Manufactured）解法，构造带界面跳跃条件的解析位移场以验证收敛率；4. 自由边层合板基准模型（如交叉铺层[0/90/90/0]和角铺层[+45/-45/-45/+45]），并与已有的三维分层共形有限元（Hierarchic FEM）结果对标。

误差度量与数值积分

研究定义了L₂范数下的位移误差（‖e_u‖_L2(Ω)）和应力误差（‖e_σ‖_L2(Ω)，以及材料界面上的层间牵引连续误差（‖e_t‖_L2(ΓI)）。对于解析算例，体积和界面积分通过单元面和单元内的高斯积分点近似；对于无解析解的案例（如悬臂梁），则通过嵌套网格策略，将粗网格的解与高精度细密参考解（由形状正则立方体网格计算）进行对比。

悬臂梁（各向同性材料）

通过1000×50×50 mm的悬臂梁模型，评估了SIPG法在高长宽比网格下的基准精度和鲁棒性。结果显示：在各项同性细化下，位移误差收敛斜率为p+1（p为多项式阶次），应力误差斜率为p，符合理论最优速率。在仅沿长度和深度方向细化（各向异性细化，宽度方向保持1个单元）时，由于弯曲问题的主要变形梯度位于x-z平面，该方法表现出更陡的表观收敛斜率（位移约2.7，应力约1.2），且即使长宽比最终增至22，求解仍保持稳定且条件良好。

制造解：分段各向同性材料

引入带材料界面（z=0.5）的两层各向同性层合板，每层具有不同的弹性模量（E₁=1, E₂=2）。通过构造连续的位移场但具有跳跃的应变/应力梯度的解析解，验证了SIPG处理材料间断的能力。结果显示：L₂位移误差收敛率为p+1，应力误差为p；界面上的层间牵引连续误差表现出与多项式阶次奇偶性相关的超收敛行为（奇数阶p+1，偶数阶p），这与高斯点处的导数误差超收敛特性一致。

制造解：分段正交各向异性材料

将上述模型扩展至正交各向异性层合板（性能见表1），同样通过构造解析解验证。结果表明，即便两层具有方向性不同的刚度耦合（如不同的E₁₁, E₂₂和G₁₂），SIPG法依然保持了最优的L₂收敛速率（位移p+1，应力p），且界面牵引连续误差的奇偶依赖性依然成立，证明了该方法在处理多层正交各向异性系统及高长宽比单元时的稳健性。

交叉铺层层合板（Cross-ply laminate）基准

复现了Ullah等人的基准算例：[0/90/90/0]对称层合板，尺寸2l×2l×4h（l=8h），受y方向宏观拉伸应变ε₀=10^-6。利用几何与载荷对称性取1/8子模型，采用4×4×4结构化六面体网格，沿自由边、铺层界面及对称面聚类节点（最大长宽比38.6），多项式阶次p=3。结果显示：沿路径分布的层间法向应力σ_zz和层间剪切应力σ_zx与参考分层共形有限元结果高度吻合，准确捕捉了自由边处的应力集中与符号反转。该模型仅用12,288 DOFs，在苹果M2笔记本上仅需4.82秒，而对比方法用了134,655 DOFs和154.02秒，显示了极高的计算效率。

角铺层层合板（Angle-ply laminate）基准

进一步测试了更具挑战性的[+45/-45/-45/+45]对称角铺层板。由于仅存一个对称面，取1/2子模型，采用8×8×4网格，并采用各向异性的多项式阶次（p_x=3, p_y=3, p_z=1）。结果显示：沿路径的σ_xy、σ_yy和σ_yz均与参考结果紧密贴合，证明了SIPG法能有效处理角度铺层带来的强面内耦合效应。此案例使用24,576 DOFs，计算耗时6.39秒。

结论

研究表明，对称内罚间断伽勒金（SIPG）法结合高长宽比六面体网格，是进行复合材料层合板高效三维分析的有力框架。它能在极其粗化的网格和拉伸单元下，保持最优的L₂收敛速率，并精确预测自由边效应引发的关键层间应力分布，计算成本较传统共形高阶有限元降低约一个数量级。该方法不仅解决了厚方向网格划分的计算瓶颈，还为后续开展层合板优化设计、自适应hp加密策略以及更复杂的工程结构分析提供了坚实的数值基础。