深海高强度钢环肋圆柱壳单位容积重量优化设计

余 俊，欧阳吕伟，李艳青，陆 波

(中国船舶科学研究中心，江苏 无锡 214082)

0 引 言

环肋圆柱壳是典型的深海耐压结构，在结构设计时，在一定容积下减轻结构重量，提高结构装载重量，是深海耐压结构设计常常关心的问题。高强度钢材料屈服强度、抗拉强度比普通钢要高，在满足同样的结构强度、使用容积的情况下，结构设计中采用高强度钢材料，可以减少钢材使用量，进而减小结构的单位容积重量，故高强度钢在深海耐压结构中被广泛应用。

为了减轻环肋圆柱壳重量，需要对其进行优化设计。国内外学者对环肋圆柱壳优化方法进行了广泛的研究。Liang 等[1]提出了在稳定性和屈服强度约束条件下，受静水压力的潜艇耐压壳体结构重量最轻的优化设计方案。由于该设计问题属于离散变量多模态约束最小值问题，作者采用回溯程序设计方法进行求解，获得的优化设计方案的重量比原有模型重量减轻6.65%。Bagheri 等[2]采用遗传算法对环肋圆柱壳进行多目标优化设计，以最大固有频率和最小结构重量为目标函数，采用权重法将多目标转化为单目标，研究了权重系数、加强筋结构形式和设计变量的取值范围对优化结果的影响。

国内丁海旭等[3]对潜艇耐压船体进行优化设计，以结构重量和满应力设计为目标，建立单目标非线性约束优化的目标规划模型，求解连续与离散混合变量的最优解，并进行规整，获得了关于满应力设计和结构重量较轻的耐压船体设计方案。杨卓懿等[4]对环肋圆柱壳进行了优化，建立了强度计算的响应面模型，并采用罚函数法进行优化，获得了结构重量最轻的优化解。王燕[5]对潜艇耐压结构优化设计方法进行探索，采用有限元法对系列潜艇环肋骨圆柱壳进行强度和稳定性分析，分别建立了分析的二次响应面模型和Kringing近似模型，然后以重量为目标，采用混合整型优化方法对环肋圆柱壳进行优化设计，并进行了对比。

目前国内外环肋圆柱壳结构优化研究均是针对具体模型的优化，主要研究优化算法和优化过程，优化结果只适用于具体的设计模型，缺少对高强度钢材料通用性的优化成果，以指导环肋圆柱壳结构设计。本文首先建立了600 MPa级高强度钢环肋圆柱壳结构优化计算方法。以多学科优化软件Isight为平台，联合利用多岛遗传算法（MIGA）和序列二次规划法（NLPQL）对高强度钢环肋圆柱壳单位容积重量进行优化，得到了最小环肋圆柱壳单位容积重量、最优的半径厚度比与最大工作压力的关系曲线和拟合公式。

1 环肋圆柱壳无量纲化的强度稳定性计算方法

适用于600 MPa级高强度钢环肋圆柱壳无量纲化的强度和稳定性计算方法，主要依据GJB/Z 21A-2001潜艇结构设计计算方法[6]和现代潜艇结构强度理论与试验[7]进行计算。计算时选用以下几个无量纲的变量、，，，，，。 其 中为 环 肋 圆柱壳耐压壳体中面半径；t 为耐压壳体厚度；为T型肋骨腹板厚度；为T型肋骨腹板高度；为T型肋骨的面板宽度；为T型肋骨面板厚度；和为无量纲化的结构参数，计算公式如下：

1.1 环肋圆柱壳无量纲化的强度计算方法

计算下列参数：

肋骨根部壳板内表面的纵向应力：

跨中壳板的中面周向应力：

肋骨平均应力：

1.2 环肋圆柱壳无量纲化的稳定性计算方法

相邻肋骨壳板欧拉失稳压力：

肋间壳板失稳临界压力：

肋骨型心具体壳板表面高度与厚度之比：

舱段间失稳临界压力：

2 优化模型和优化算法

2.1 设计变量

查阅相关资料，设计变量取值范围如下：

2.2 目标函数

单位容积重量为目标函数：

2.3 约束条件

1）强度约束

跨中壳板中面周向应力：

肋骨根部壳板内表面纵向应力：

肋骨平均应力：

2）稳定性约束

肋间壳板失稳临界压力：

舱段总体失稳临界压力：

3）几何约束

肋骨高度与壳板半径之比：

2.4 优化算法

该优化问题属于连续型变量约束优化问题，联合运用多岛遗传算法（MIGA）和序列二次规划法（NLPQL）进行优化，采用多岛遗传算法进行全局搜索，得到一个最优解，然后在该最优解基础上采用进一步局部寻优，从而得到满足约束条件单位容积重量最轻的设计方案。

多岛遗传算法是全局优化算法的一种，它能求解全局最优解，避免了在局部区域搜索；同时它只评价设计点，不计算任何函数的梯度，适合处理离散型和连续性设计变量。相比传统遗传算法，它将每一代的种群划分为几个“岛”，每个“岛”都进行传统遗传算法操作，岛与岛之间有个体发生“迁移”，这使其具有更优良的全局求解能力。

序列二次规划法（NLPQL）能够处理设计变量实型问题。NLPQL算法将目标函数以二阶泰勒级数展开，并把约束条件线性化，通过解二次规划得到下一个设计点。该方法能够有效地探索初始设计点周围局部区域，如果设计空间是连续单峰的形态，能够沿最快下降方向快速搜索。

2.5 优化过程

本文采用Isight软件对强度和稳定性计算模型进行集成，然后在Isight软件优化算法下进行驱动优化，图1为Isight对环肋圆柱壳进行确定性优化的流程图。

3 600 MPa 级高强度钢单位容积重量优化结果

表1列出了600 MPa级高强度钢环肋圆柱壳单位容积重量优化结果，从表1可以看出，时，在同一工作压力下，最小单位容积重量及对应化不大。

从统计学角度，一项试验中欲考察的变量称为因子或因素，因子被考察的范围称为试验范围，试验结果称为响应（Response）或输出（Output），在试验范围内，因素所取的值称为水平。常用的试验设计方法有全因子设计、正交设计、中心组合设计、拉丁方设计和最优拉丁方设计等。对于因子数目较大时，可以采用拉丁方设计或最优拉丁方设计[8]，拉丁方设计中每个因子的设计空间等分份，然后每个因子水平随机组合共确定个采样点（一个因子每个水平只研究一次）。而最优拉丁方法是在拉丁方法基础上，对每个因子组合方式进行优化，使设计点在设计空间均匀分布，如图2所示。

图1 环肋圆柱壳确定性优化流程图Fig.1 Deterministic optimization process of ring-stiffened cylindrical shell

表1 600 MPa级高强度钢环肋圆柱壳单位容积重量优化结果Tab.1 Weight to Volume Optimization results of the 600 MPa ring-stiffened cylindrical high-strength steel shell

本文采用最优拉丁方进行试验设计，取500个取样点，对每个采样点进行单位容积重量计算，然后进行各变量之间的相关性分析，得到单位容积重量和自变量的复相关系数，如图3所示。（）与自变量的复相关系数图如图4所示，与自变量的复相关系数图如图5所示，与自变量的复相关系数图如图6所示。图中纵坐标为因变量相关系数，横坐标为自变量。

图2 拉丁方设计（左）和最优拉丁方设计（右）Fig.2 Latin square design (left) and optimal Latin square design (right)

图3 与自变量的复相关系数Fig.3 Multiple correlation coefficient between and independent variable

图4 （）与自变量复相关系数Fig.4 Multiple correlation coefficient between ( ) and independent variable

图5 与 自变量复相关系数Fig.5 Multiple correlation coefficient between and independent variable

图6 与自变量复相关系数Fig.6 Multiple correlation coefficient between and independent variable

图7 600 MPa 级高强度钢环肋圆柱壳与关系Fig.7 Relation between and of the 600 MPa ring-stiffened cylindrical high-strength steel shell

与1接近，说明拟合精度较高，该公式可以用于600 MPa高强度钢环肋圆柱壳最小的单位容积重量计算。

图8 600 MPa级高强度钢环肋圆柱壳Fig.8 Relation betweenan of the 600 MPa ringstiffened cylindrical high-strength steel shell

4 结 语

1）本文依据GJB/Z 21A-2001潜艇结构设计计算方法建立了600 MPa级高强度钢环肋圆柱壳强度和稳定性无量纲的计算方法，该方法可用于环肋圆柱壳强度和稳定性校核和优化设计。

2）在上述计算方法的基础上对环肋圆柱壳单位体积重量进行优化，通过多岛遗传算法（MIGA）整体寻优和序列二次规划法（NLPQL）局部寻优，得到了在不同的圆柱壳长度半径比、不同的最大工作压力下，600 MPa级高强度钢环肋圆柱壳最小的单位容积重量和对应的最优的半径厚度比，在圆柱壳长度半径比为1，2，3，4情况下分别绘制了600 MPa级、高强度环肋圆柱壳最大工作压力与最小的单位容积重量关系曲线，最大工作压力与最优的半径厚度比关系曲线，并证明了最小的单位容积重量、最优的半径厚度比仅和材料、最大工作压力有关，和圆柱壳长度半径比关系不大。

3）根据上述优化结果，提出了计算600 MPa级高强度环肋圆柱壳最小的单位容积重量、最优的半径厚度比的拟合公式，该公式可以用于评估高强度钢环肋圆柱壳单位容积重量是否最优，可以计算单位容积重量最优情况下的半径厚度比，指导深海耐压环肋圆柱壳优化设计与评估。