《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》：Gradient-based initial postbuckling optimisation of imperfection-sensitive variable angle tow shells using an implicitly differentiated reduced-order model
编辑推荐：
轻量化复合材料壳体结构的后屈曲行为优化仍然是一个关键且很大程度上未解决的挑战，原因在于结构非线性、对不确定数据的敏感性以及通常大量的设计变量的综合作用。即使是很小振幅的几何缺陷也会显著改变后屈曲响应，因此必须在优化框架内明确考虑。基于Koiter渐近分析（Ko
轻量化复合材料壳体结构的后屈曲行为优化仍然是一个关键且很大程度上未解决的挑战，原因在于结构非线性、对不确定数据的敏感性以及通常大量的设计变量的综合作用。即使是很小振幅的几何缺陷也会显著改变后屈曲响应，因此必须在优化框架内明确考虑。基于Koiter渐近分析（Koiter's asymptotic analysis）的降阶模型（Reduced-Order Models, ROMs）为初始后屈曲阶段和缺陷敏感性提供了高效且准确的描述。然而，在涉及大量设计变量的梯度优化问题中，由于精确且高效的梯度评估的缺失，其系统应用仍受到阻碍。
研究人员通过引入一种专门针对多模态Koiter渐近分析与有限元固体壳离散化相结合的隐式算法微分策略，解决了上述局限性。所提方法应用于变角度铺丝（Variable Angle Tow, VAT）复合材料壳体结构的后屈曲优化，包括带开孔的构型。数值结果表明，与基于有限差分的方案相比，特别是存在缺陷时，所提解析梯度在稳健性和计算效率方面均有显著提升。
本研究发表于《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》，聚焦于复合材料壳体结构后屈曲优化这一长期存在的核心难题。壳体结构凭借其高比刚度与比强度，在航空航天、风电等领域具有巨大应用潜力，但其高柔度比导致结构极易发生几何非线性屈曲现象，且后屈曲响应复杂，存在多平衡分支与强模态耦合。更为关键的是，微小的几何缺陷、材料属性偏差或载荷不确定性都会严重改变后屈曲行为，使结构趋向某一特定平衡分支，形成所谓的"缺陷敏感性"。传统的线性化屈曲分析（Linearised Buckling Analysis）完全无法捕捉此类复杂现象，若用于设计将导致对承载能力的显著高估。尽管路径跟踪策略（Path-Following Strategies）和基于物理信息神经网络的方法被尝试用于后屈曲优化，前者计算成本过高以至于难以考虑缺陷敏感性，后者则局限于特定问题类别。基于Koiter弹性稳定理论的降阶模型凭借其独特优势——可在完美结构分析完成后引入缺陷进行高效敏感性分析——成为后屈曲分析的理想工具，但其在梯度优化中的瓶颈在于：针对大量设计变量，缺乏精确的梯度解析表达式，传统有限差分法需为每个设计变量重建降阶模型，计算量巨大且对扰动参数高度敏感。
研究人员首次提出了一种隐式算法微分策略，实现了Koiter降阶模型中所有组成成分的解析微分，包括线性弹性解、屈曲模态与载荷、以及二次修正项。该方法利用假设应力固体壳模型（Assumed-Stress Solid-Shell Model）采用Green-Lagrange应变度量，使得混合能量对模型变量呈三次多项式依赖，从而获得显著的代数简化。优化目标为改善初始后屈曲行为：降低缺陷敏感性、在可能时实现初始稳定行为并提升后屈曲承载能力；若不可行，则执行极限载荷最大化。研究通过多个数值算例验证了方法的有效性，包括带开孔的曲面壳体，结果表明解析梯度显著优于有限差分法。
研究采用的关键技术方法包括：（1）基于混合假设应力固体壳有限元离散化，采用8节点六面体单元，通过Hellinger-Reissner格式的混合应变能克服剪切、薄膜和体积自锁；（2）Koiter多模态渐近降阶模型构建，将配置表示为线性解、屈曲模态及二次修正的组合，通过Taylor展开获得降阶平衡方程；（3）隐式算法微分策略，对Koiter各能量贡献项进行精确微分，利用已分解的系数矩阵求解一系列线性系统获得解析梯度；（4）最坏工况缺陷识别算法，采用随机采样结合Gumbel极值分布拟合确定最坏几何缺陷组合；（5）全局收敛移动渐近线法（Globally Convergent Method of Moving Asymptotes, GCMMA）作为梯度优化求解器，配合多起点策略。样本来源为数值构造的复合材料壳体结构算例，无临床或生物样本队列。
**矩形板算例**：首先验证单层堆叠序列在单模态Koiter展开下的梯度精度。数值实验表明，单模态时解析梯度与有限差分结果一致。当考虑双模态并引入缺陷敏感性时，目标函数因缺陷引入产生微小振荡，有限差分梯度（尤其是小步长h时）出现严重非物理振荡，仅h=1时获得平滑但略欠准确的结果；而解析梯度不受影响。优化求解表明，解析梯度在5步内收敛至最优解，有限差分在小h时无法收敛。随后对八层对称堆叠序列±θ?/±θ??进行双变量优化，结果再次证实解析梯度在缺陷存在时的优势：有限差分近似因目标函数噪声而失效，解析梯度则保持稳健。
**变角度铺丝方板算例**：边长1m、厚20mm的单层方板，采用20×20网格，10×10网格描述VAT铺向角分布。优化目标为Koiter降阶空间中的后屈曲曲率b=B????/λ?。与Henrichsen等人的单模态Koiter优化方法对比，不同初始铺层下的最优解和平衡路径均显示出良好一致性，验证了方法在多模态条件下的正确性。
**曲面壳板算例**：曲率半径1000mm、厚度10mm的简支曲板，采用32×32固体壳有限元网格。研究涵盖从两变量（S1-SF）到六变量（S3-SF）直至变量角度铺丝大量设计变量（S1-VAT达338个变量，S2-VAT达676个变量）的多种配置。参数研究表明，目标函数高度非线性且非凸，存在多个局部最优，且分为不稳定区域（λ?≡λ???）和稳定区域（λ?=λ?）两部分。多起点随机初始化策略被采用。结果显示，解析梯度在多数情况下以更少的迭代步数和函数评估次数收敛至更优或相同解；有限差分在接近最优时因导数趋零而迭代次数激增。最为显著的是，S2-VAT案例的优化时间从有限差分的3.87小时骤降至解析方法的26分钟，效率提升近9倍。VAT技术相比直纤维层合板分别带来约24%（解析）和19%（数值）的性能提升，且能实现对称铺层下的稳定最优解。
**带开孔曲面壳板算例**：在曲板基础上引入三个椭圆形开孔（半轴25cm×12.5cm），增加几何复杂性。S1-SF和S3-SF案例表明，开孔导致不稳定区域扩大、设计载荷整体降低，且解析梯度始终获得优于数值梯度的最优解。S4-VAT案例仅采用解析梯度，结果显示VAT层合板实现了比SF配置显著更高的设计载荷（2.20 vs 1.15），且从初始不稳定响应优化为稳定后屈曲行为，Koiter降阶解与Riks全阶路径跟踪解吻合良好。
**风力发电机塔架算例**：高10m圆柱段加锥段的工业级壳体结构，承受顶部压缩和剪切载荷，采用平衡对称八层铺层±θ?/±θ?/±θ?/±θ??。从初始016铺层（λ?=0.3511）优化至±83/±50/±16/?5?，极限载荷提升30%至0.4570。模态收敛研究表明，考虑16个屈曲模态时Koiter分析与路径跟踪参考解高度一致，验证了方法的准确性和工程适用性。
在讨论部分，主要结论包括：所提隐式算法微分策略首次实现了Koiter降阶模型解析梯度的精确计算，克服了有限差分法在缺陷敏感优化中的根本局限；解析梯度不仅将优化时间降低一个数量级（如从3.87小时降至26分钟），更通过消除数值噪声和步长敏感性显著提升了优化鲁棒性，尤其在从非稳定区向稳定区过渡时表现关键；该方法使高维设计空间中复杂后屈曲行为的梯度优化成为可能，包括传统方法难以处理的变角度铺丝结构和带开孔几何；Koiter降阶模型在多达16个模态时仍保持良好精度，计算成本可控，为工程应用提供了实用工具。研究指出，该方法具有向拓扑优化扩展的潜力，这将是未来工作方向。
