水下非均匀环肋圆柱壳多学科设计优化

李静茹1，刘江2，黎胜*2 1海南大学机电工程学院，海南海口 570228

2大连理工大学工业装备结构分析国家重点实验室船舶工程学院，辽宁大连 116024

摘 要:［目的］水下结构对强度、质量和振动噪声等的要求严格，为提高结构整体性能，从多学科角度优化

水下结构设计具有重要意义。［方法］首先，采用改进的目标传递方法和代理模型优化水下结构振动声辐射设计，即采用改进的目标传递方法（将拉格朗日对偶原理运用到传统分级目标传递法中）实现对设计问题的分级和不同学科性能下的并行优化，并构建Kriging代理模型来解决水下结构振动声辐射计算耗时的问题；然后，以水下环肋圆柱壳（水下航行器）为例，针对轻量化、强度、振动与声辐射方面的要求，选取6个设计变量，以环肋圆柱壳的质量、应力、共振频率及其所处位置的声功率级作为目标函数，建立并行模型进行多学科设计优化。［结果］结果表明， Kriging代理模型可准确和实时预报空间任一点的响应；改进的目标传递法不仅获得了多学科优化结果，且与传统传递方法相比收敛性提高了30%。［结论］研究表明，基于所提方法

得到的多学科设计优化结果更理想且收敛性更好。关键词：非均匀环肋圆柱壳；改进的目标传递法；代理模型；多学科设计优化

中图分类号: U661.44文献标志码:A DOI：10.19693/j.issn.1673-3185.02111 Multidisciplinary design optimization of non-uniform stiffened cylindrical shells LI Jingru1, LIU Jiang2, LI Sheng*2

1 Mechanical and Electrical Engineering College, Hainan University, Haikou 570228, China

2 State Key Laboratory of Structural Analysis for Industrial Equipment, School of Naval Architecture and

Ocean Engineering, Dalian University of Technology, Dalian 116024, China

Abstract: ［Objectives］Underwater structures have strict requirements for intensity, mass, vibration, noise and so on. From a multidisciplinary perspective, realizing the design optimization of such structures to improve their overall performance is of great significance.［Methods ］The enhanced analytical target cascading method and surrogate model are employed for the design optimization of underwater structure vibration and acoustic radiation: an enhanced target cascading (ATC) method based on the Lagrange duality theory and traditional ATC method is established to achieve the classification of problems and parallel optimization of the performance of different disciplines; and the Kriging surrogate model is used to solve the time-consuming problem in the calculation of the vibration and acoustic radiation of underwater structures. Underwater nonuniform stiffened cylindrical shells (i.e., underwater unmanned vehicles) are selected for multidisciplinary design optimization (MDO): six design variables are selected, and the structural mass, stress components, resonant frequency and corresponding radiated power level are regarded as objectives for establishing the parallel optimization model when considering the vibration and acoustic radiation characteristics and lightweight requirements.［Results］The results show that the Kriging surrogate model can accurately predict the response in the design space, while the enhanced ATC method can obtain effective multidisciplinary optimization results with a 30% increase in convergence properties.［Conclusions ］This study shows that multidisciplinary design optimization based on the enhanced ATC method and surrogate model not only obtains effective results but also has better convergence properties.

Key words: non-uniform stiffened cylindrical shell；enhanced analytical target cascading method；surrogate model；multidisciplinary design optimization

0 引 言

舰艇声隐身总体设计是一个多学科设计优化问题，涉及了结构力学、振动、声学等学科。目前，常用的设计方法是串行设计方法，此方法人为割裂了各学科间的耦合关系，得到的设计结果会降低系统的总体性能。多学科设计优化 (multidisciplinary design optimization, MDO)是近几年发展起来的旨在提高复杂工程问题优化设计效率的一种分解策略。有别于传统的串行设计方法，多学科设计优化采用的是并行设计方法，不仅可提高效率，还可提高系统的整体性能。实际上，多学科设计优化是一种通过探索和利用工程系统中相互作用的协调机制来设计复杂产品及其子系统的方法[1] ，其主流算法是基于分解技术的多级优化算法[2] ，而分级目标传递法 (analytical target cascading, ATC ）就是典型算法之一，其收敛性已得到证明，可以收敛到全局最优点[3]。

多学科设计优化在水下结构振动声辐射领域的应用研究比较少。陈炉云等[4] 基于 iSIGHT 多学科优化设计软件，以结构模态频率和节点加速度为目标函数，对环肋圆柱壳进行了结构优化。黎胜[5]基于传统目标传递法进行了水下加筋板振动声辐射设计优化，证明了多学科设计优化方法应用的有效性。但是，相对于水下加筋板模型，实际设计中的结构都比较复杂，设计变量和目标函数会更多，从而大大增加了传统目标传递法的收敛难度。因此，本文采用改进的目标传递法[6]，以典型的水下环肋圆柱壳结构（例如水下航行器）为例，验证此方法在保证总体性能优化的同时具有更好的收敛性。

对于水下结构而言，动态性能的计算依赖于高精度的数值模型（例如有限元/边界元耦合方法），此类模型计算量大且耗时，无法完成整个过程的优化，而使用代理模型（例如 Kriging 代理模型）是解决此问题的可行方法。鉴于代理模型可以较好地实时预报水下相对复杂结构的振动声辐射性能，本文以 Kriging 代理模型作为优化过程中的分析模型来代替高精度数值模型。

舰艇总体设计对于轻量化、强度以及振动和声辐射等多个学科方面都有严格的要求。本文拟通过研究水下环肋圆柱壳的多学科设计优化方法，在满足多学科要求的前提下，为最优设计方法研究提供思路。首先，利用有限元/边界元方法（FEM/BEM ）计算样本点响应数据，建立 Kriging代理模型，以此作为设计优化过程中的目标响应分析模型；然后，将代理模型作为分析模型，分别利用传统和改进的目标传递方法对水下非均匀环肋圆柱壳进行设计优化；最后，通过结果对比，验证所用方法的有效性和优越性。

1 理 论

1.1改进的目标传递法目标传递法[7] 是一种基于分解技术的多层次、多学科优化策略，适用于大规模的复杂工程问题，具有响应和连接变量在不同层级之间按顺序传递及反馈的特点。不同于其他多学科优化算法（例如协同优化），目标传递法计算信息的传递发生在上、下两层级完全优化之后，所以显著减少了学科级优化的次数[8]。

为了增强分解后上、下层级之间参数变量的一致性，需要对目标传递法中的偏差项作一定的处理。通常，是在偏差项上乘以一个较大的正数，在相同容许误差下，偏差项本身会变得相对较小。但是，此处理方法也有一个弊端，即大权重系数会在最优解附近导致数值的不稳定性。因此，本文摒弃采用大权重系数，转而在目标传递法中设置拉格朗日乘子，利用对偶原理求解。

设置了拉格朗日乘子的改进目标传递法的数学表达式被分解为系统级和学科级问题来表达。

系统级问题：

式中： x˜ 为系统级的本地变量； ys,k ， Rs,k分别为系

sys

统级的第k个子系统的连接变量及响应值，其中，上标L表示学科级传递到系统级的变量，下标s表示子系统的变量； Rsys ， Tsys分别为系统响应和系统目标； λ R ， λ y分别为与第k个子系统响应偏差和连

k k接变量偏差相关的拉格朗日乘子；Csys为系统级问题的子问题集； c R和c y分别为与第k个子系统的响

k k

应偏差和连接变量偏差相关的惩罚因子； g和h分

别为不等式和等式约束条件；符号“◦”表示两向

量逐项相乘，例如：

(a1 , a2 , ··· , an ) ◦ (b1 , b2 , ··· , bn ) = (a1 b1 , a2 b2 , ··· , an bn )

式中： x˜s, ， ys, 和Rs, 分别为学科级第j个子系统的

j j j本地变量、连接变量及响应值，其中上标U表示从系统级传递到学科级的变量； λ R， λ y分别为与学

j j

科级第j个子系统响应偏差和连接变量偏差相关的拉格朗日乘子； c R和cj y分别为与学科级第j 个子

j系统响应偏差和连接变量偏差相关的惩罚因子。图 1所示为上述信息交换的流程。

1.2 Kriging 代理模型

利用数值模型求解水下非均匀环肋圆柱壳的振动声辐射问题所需的计算非常耗时，且计算量也过大，所以在设计优化过程中一般采用代理模

型来替代高精度的数值模型。

常用的代理模型构建方法包括多项式响应面法、径向基方法、人工神经网络方法和Kriging 方法等[10-13]。目前，Kriging 方法已成为用于多学科设计优化的一种具有代表性的代理模型。构建

Kriging 代理模型一般需要采取如下步骤：

1)确定设计变量，通过实验设计方法得到构建代理模型所需的样本点；

2)利用高精度数值模型，计算每个样本点相应的响应值；

3)根据样本点及其相应的响应值构建Kriging代理模型；

4)在设计空间内任选一点，比较分别经代理模型和高精度数值模型计算得到的结果，以验证代理模型的精度。

Kriging代理模型的一般数学表达式为

φˆ (x) = F (β:r , x) + zr (x) (7) r

式中： φˆ (x)为设计空间任一点x 的响应值，下标

r r为所研究的响应值的数量，r = 1, , q； F (β:r , x) =

··· β1,r f1 (x) + + β fp (x)，其中 f1 (x) , , fp (x)为事先

··· ···

p,r

选定的基函数， β:r为按照最优线性无偏估计所得到的参数；zr (x)为均值为 0的随机过程。为便于叙述，假定q = 1，则式（7）可变为

o(x) = [O (θ, x1 , x) , , O (θ, xm , x)] ···

(CB + HB) ps = GB vn (12)

式中： CB ， HB和GB为系数矩阵； ps为结构表面声压向量； vn为结构表面的速度向量，且vn = iωQs u（其中ω为角频率， Qs为结构表面法线方向矩阵），故式（12）代数方程可改写为

(CB + HB) ps = iωQs GB u = Gu (13)

式中，G为转换矩阵。通过结构和流体耦合的有限元方程，可以求出结构表面节点处的位移，进而利用边界元方法计算出结构表面声压向量ps ，则结构的辐射声功w

1

( ∗)

率为Π = Re ps v dS，其中， v ∗表示vn的共轭复

n n

2 S

数， Re表示取实部。

2算例分析

2.1环肋圆柱壳参数设置

图 2所示为建立的非均匀环肋圆柱壳有限元模型，其两端截面的最外圆周线固支。环肋圆柱壳的非均匀是指肋骨尺寸、肋距不等。本文研究中肋骨布置关于中间对称，共19 根，由左至右依次编号。其中，编号 8，9，10，11，12 的肋骨腹板高度最大，编号 13，14，15，5，6，7 的腹板高度次之，剩余编号的肋骨腹板高度最小。肋骨截面形状为 T型截面，肋距为中间小、两端大。为方便优化设计，将肋距写成等差数列形式来完成非均匀设定。但是，不同于文献[15]，本文仅对几种不同的离散情况进行了分析，即将3种肋骨腹板高度（h1，h2，h3）都作为连续变量来确定其最优解。图2 非均匀环肋圆柱壳有限元模型

Fig. 2 FEM model of non-uniform stiffened cylindrical shells

如上所述，非均匀环肋圆柱壳有限元模型长12 m ，直径 6m。3 种肋骨尺寸(单位：mm) 分别24 h1 14 h2 8 h3

× × ×

为：⊥ ⊥ ⊥

, , , h1 ⩾ h2 ⩾ h3。肋30 120 26 80 12 60

× × ×

距关于中间肋骨对称，并采用等差数列表示方法，可用2个参数代表：从中间编号为10 的肋骨开始，第1 档肋距l以及向两端依次增加的公差∆l。因此，距中间肋骨的第n档肋距可表示为： ln = l + (n 1) ∆l (n = 1, 2, , 9)。圆柱壳为钢材料： − × ···密度ρs = 7 850 kg/m3 ，弹性模量E = 2.1×1011 N/m2，泊松比υ = 0.3，屈服极限σs = 588 MPa。流体介质为水，其密度ρ = 1 025 kg/m3 ，声速c = 1400m /s。激励为幅值大小1N的垂向简谐力，作用位置位于距左端 6.3 m处。计算静水压力下的应力分量时，工作载荷Pj = 4.26 MPa。

非均匀环肋圆柱壳的设计变量取为6 个，分别是从中间肋骨开始的第1 档肋距l、肋距公差∆l、3 种肋骨腹板高度h1 ， h2 ， h3以及圆柱壳壳板厚度t。表 1所示为6个变量的取值范围。

考虑轻量化、强度以及振动和声辐射方面的要求，建立了6个目标函数，包括：环肋圆柱壳的质量M；整个肋距范围内的壳板横剖面正应力最大值σ1；壳板纵剖面正应力最大值σ2；肋骨横剖面正应力的最大值σf ；结构 1-2 阶共振频率f1；结构1-2阶共振频率处的声功率W。按照优化目标，希望在小于屈服极限的情况下最大化3个应力分量和 1-2 阶共振频率，以及最小化共振频率处的声功率级和圆柱壳质量。

根据目标传递法的要求，将此优化问题分为系统级和学科级。分解前，首先可以求出给各学科的最优响应，并将其作为理想目标值；然后，建立系统级和学科级的优化流程。系统级优化目标是协调各学科的一致性，以及最小化所有学科目标函数与理想目标间的偏差；学科级优化目标是最小化各学科与系统级传递过来的响应和连接变量之间的偏差值。值得注意的是，对于质量学科而言，连接变量只有后4个设计变量。图3 所示为非均匀环肋圆柱壳多学科优化的组织结构。图3 非均匀环肋圆柱壳多学科优化组织结构图Fig. 3 Flowchart of MDO in non-uniform stiffened cylindrical shells

2.2 建立 Kriging 代理模型

建立 Kriging 模型与选取的样本点有着重要关系。样本点数量、空间相关性影响着Kriging

模型的精度。科学合理地选取样本点是实验设计的主要内容。本文选用最优对称拉丁超立方体抽样方法来选择样本点。使用此方法选取的样本点数量一般为设计变量个数的4倍。当然，如果预测的响应在设计空间内变化剧烈，也可以多选择一些样本点来达到更好的全局近似。本文选取24个样本点，利用有限元/边界元方法计算24 种情况下的响应值，并建立了实时预报设计空间内任一点响应值的代理模型，再选取(0.1 ，0.11， 0.27，0.25，0.12，0.028) 和 (0.09，0.099，0.243，0.225， 0.108， 0.025 2) 这 2个非样本点（分别记为预测点1 和2）来检验所建代理模型的精度。表2示出了采用有限元耦合方法求得的各响应值与基于代理模型得到的响应值的对比结果。由表可见，基于代理模型的计算结果与常规数值模型方法所得结果最大相差3.5%，最小相差 0.084 2%。最大误差在允许范围内[16] ，这说明由24 个样本点建立的Kriging模型预测精度较高，在设计优化过程中可以代替数值模型作为目标响应值的分析模型。2.3基于改进目标传递法的优化结果

将所建的 Kriging 代理模型作为分析模型，根据 1.2节所述步骤求解。本文改进的目标传递法设置的容许误差为0.01。值得注意的是，为了尽量减少各优化目标因不同的绝对值及其取值范围带来的不同相对权重，避免结果“偏向”某一学科，本文将各目标值除以目标初始值进行了无量纲化，旨在一定程度上减少这种倾斜，然而仍然不能完全消除此影响[4]。

表 3为分别通过改进的目标传递法和传统的目标传递法得出的结果，即通过改进目标传递法得到的6个设计变量优化结果为 (0.1，0.111，0.270 4， 0.254，0.121， 0.022 4) ，记为结果1；通过传统的目标传递法得到的优化结果为 (0.082 5，0.09，0.217， 0.201，0.097， 0.026 6) ，记为结果2。由表可见：相较于初始值，两种目标传递法在6个目标函数上

都达到了更好的优化效果：应力分量、1-2阶共振频率更大，共振频率处的声功率级和圆柱壳质量更小；在整个肋距范围内，对于包括壳板横剖面正应力最大值σ1、壳板纵剖面正应力最大值σ2、肋骨横剖面正应力的最大值σf和结构 1-2 阶共振频率 f1这 4个目标函数，改进的目标传递法比传统的目标传递法都得到了较大的值，即优化效果更好；改进的目标传递方法对声功率级（共振频率处的声功率级） W和环肋圆柱壳质量M目标函数的优化结果比传统的目标传递法要大，即优化效果较差，但增大的最大幅度为声功率级目标，只有0.9%。当传统的目标传递法达到收敛后，计算第1级和第2级每个学科的响应以及相应连接变量的偏差，用向量的二范数度量。由于学科级有 6个不同的目标，因此可以得到6个偏差值(壳板横剖面正应力最大值偏差∆σ1、壳板纵剖面正应力最大值偏差∆σ2、肋骨横剖面正应力的最大值偏差∆σf、结构 1-2阶共振频率偏差∆ f1、共振频率处的声功率级偏差∆W、环肋圆柱壳质量偏差∆M ) ，如表 4 所示。由表可见，最小的偏差值为0.014 3，而改进的目标传递法得到的最小偏差值小于 0.01才会达到收敛标准。可见，利用拉格朗日对偶原理得到的改进目标传递法其收敛性强于传统的目标传递法。

综合表3 和表4的数据可知：改进的目标传递法在壳板横剖面正应力最大值σ1、壳板纵剖面正应力最大值σ2、肋骨横剖面正应力的最大值σf和结构 1-2 阶共振频率f1这 4个目标函数上与传统的目标传递法相比，达到了更好的优化效果，且其他2个目标函数的优化结果与传统的目

标传递法相差无几；在此基础上，改进的目标传递法比传统的目标传递法还具有更好的收敛性。因此，本文提出的利用拉格朗日对偶原理的改进目标传递法优于传统的目标传递法。

3 结 语

本文利用 Kriging 代理模型和改进的目标传递法完成了水下非均匀环肋圆柱壳的多学科优化。首先，根据水下航行器的轻量化、强度以及振动和声辐射方面的要求，选取了第1 档肋距、肋距公差、3种肋骨腹板高度、圆柱壳板厚度作为设计变量，质量、应力、共振频率、共振频率处的声功率级作为目标函数进行优化；然后，选用最优对称拉丁超立方体抽样方法选取样本点，利用有限元/边界元方法对样本点的响应进行数值计算，利用计算结果建立 Kriging 代理模型，并计算验证了其精度；最后，为了保证优化结果的收敛性，利用引入拉格朗日乘子的改进目标传递法，基于6个设计变量，对6个目标函数进行了优化设计。结果表明，代理模型预报的结果与常规数值模型所得结果最大相差3.5%，因此，采用 Kriging代理模型来代替高精度的数值计算模型有效地解决了计算量过大、计算时间过长的问题，且具有较高的精度；使用改进的目标传递法可以得到与传统的目标传递法基本相同的优化结果，且其收敛性相较于传统的目标传递法提高了30%。

水下非均匀环肋圆柱壳是水下航行器的经典简化模型，其性质和工况相比于实际的水下航行器结构比较简单，后续还可利用本文方法，对典型的实际水下航行器结构进行优化设计，以检验本文方法在实际生产中的适用性，以增强分析方法的稳定性，最终达到在实际生产中指导水下航行器设计的目的。

参考文献：

[1] SOBIESZCZANSKI-SOBIESKI J, HAFTKA R T. Multidisciplinary aerospace design optimization: survey of recent developments[J]. Structural Optimization, 1997, 14(1): 1–23.

[2] 姜哲, 崔维成.多学科设计优化算法比较及其在船舶和海洋平台设计上的应用[J]. 船舶力学, 2009, 13(1): 150–160.

JIANG Z, CUI W C. A comparison of multidisciplinary design optimization algorithms and their application to the design of ships and offshore platforms[J]. Journal of Ship Mechanics, 2009, 13(1): 150–160 (in Chinese). [3] MICHELENA N, PARK H, PAPALAMBROS P Y. Convergence properties of analytical target cascading[J]. AIAA Journal, 2003, 41(5): 897–905.

[4] 陈炉云, 王德禹. 基于 iSIGHT 的环肋圆柱壳动力学优

化设计 [J]. 中国舰船研究, 2007, 2(6): 1–3, 8. CHEN L Y, WANG D Y. Dynamic optimization design of ring-stiffened cylindrical shell based on iSIGHT[J]. Chinese Journal of Ship Research, 2007, 2(6): 1–3, 8 (in Chinese).

黎胜.基于目标传递和代理模型的水下结构振动声辐射设计优化 [J]. 振动与冲击, 2010, 29(10): 144–147.

LI S. Design optimization for underwater structural vibration and acoustic radiation based on target cascading and surrogate model[J]. Journal of Vibration and Shock, 2010, 29(10): 144–147 (in Chinese).

KIM H M, CHEN W, WIECEK M M. Lagrangian coordination for enhancing the convergence of analytical target cascading[J]. AIAA Journal, 2006, 44(10): 2197–2207.

[7] KIM H M, FELLOW R, MICHELENA N F, et al. Target cascading in optimal system design[J]. Journal of Mechanical Design, 2003, 125(3): 474–480.

[8] 吴蓓蓓, 黄海, 吴文瑞. ATC 与 CO 方法对比及其在卫星设计问题中的应用[J]. 计算机工程与设计, 2012, 33(6): 2455–2460.

WU B B, HUANG H, WU W R. Comparison of ATC with CO and its application in satellite design problem[J]. Computer Engineering and Design, 2012, 33(6): 2455–2460 (in Chinese). [9] BERTSEKAS D P. Nonlinear Programming[M]. 2nd ed. Belmont: Athena Scientific, 1999.

[10] SIMPSON T W, BOOKER A J, GHOSH D, et al. Approximation methods in multidisciplinary analysis and optimization: a panel discussion[J]. Structural and Multidisciplinary Optimization, 2004, 27(5): 302–313.

[11] JIN R, CHEN W, SIMPSON T W. Comparative studies of metamodelling techniques under multiple modelling criteria[J]. Structural and Multidisciplinary Optimization, 2001, 23(1): 1–13.

[12] SIMPSON T W, POPLINSKI J D, KOCH P N, et al. Metamodels for computer-based engineering design: survey and recommendations[J]. Engineering with Computers, 2001, 17(2): 129–150.

[13]曾会华, 余雄庆.基于代理模型的气动外形优化[J]. 航空计算技术, 2005, 35(4): 84–87.

ZENG H H, YU X Q. An approach to aerodynamic shape optimization using surrogate models[J]. Aeronautical Computer Technique, 2005, 35(4): 84–87 (in Chinese).

[14]王林.深海耐压结构型式及稳定性研究 [D]. 北京: 中国舰船研究院, 2011.

WANG L. Studies on deep-sea pressure structural type and stability[D]. Beijing: China Ship Research and Development Academy, 2011 (in Chinese).

[15]陈美霞, 徐鑫彤, 魏建辉, 等.非均匀圆柱壳振动及声辐射特性优化设计 [J]. 船舶力学, 2013, 17(1): 164–170. CHEN M X, XU X T, WEI J H, et al. Optimization design of vibro-acoustic behavior of non-uniform stiffened cylindrical shell[J]. Journal of Ship Mechanics, 2013, 17(1): 164–170 (in Chinese).

[16]郭明慧, 黎胜.运用代理模型方法预测潜艇结构模型的振动声辐射 [J]. 中国舰船研究, 2013, 8(6): 69–74. GUO M H, LI S. A surrogate model for structural vibration and acoustic radiation of underwater submarine structures[J]. Chinese Journal of Ship Research, 2013, 8(6): 69–74 (in Chinese). [5] [6]