高速单体船航行性能与结构特性综合优化方法研究

《高速单体船航行性能与结构特性综合优化方法研究》由会员分享，可在线阅读，更多相关《高速单体船航行性能与结构特性综合优化方法研究（11页珍藏版）》请在装配图网上搜索。

1、高速单体船航行性能与结构特性综合优化方法研究杨松林，李密, 陈淑玲摘 要本文以典型高速单体船为基础，根据一些实船的资料和相关经验公式并结合高速单体船的结构特点，探讨一种将船舶航行性能与结构力学特征两者综合优化的方法。首先依次建立了基于 (VC+6.0)软件的高速单体船快速性优化、操纵性优化、耐波性优化的数学模型，并包括了这三个方面的优化变量、目标函数和约束条件。在选择算法方面，本文将遗传算法与混沌算法复合构造了一种遗传混沌算法，并利用C+语言对其编制了友好界面。大量的计算数据表明：遗传混沌算法不仅能克服遗传算法的不足，而且能减少计算、提高效率。其次，基于舰船通用设计规范（GJB 4000-20

2、00）对高速单体船舯横剖面进行结构优化设计。将舯剖面上影响其剖面模数的主要板厚、骨材型号等作为设计变量，以舯剖面纵向构件的单位长度重量最轻为目标函数，并根据相关规范提取相应的约束条件，采用遗传混沌算法对舯剖面的主要构件进行优化计算。最后，把船舶航行性能与结构力学特征综合优化的遗传混沌算法结果，与其它优化方法（例如，遗传算法、混沌算法）的优化结果进行比较。经比较可知，遗传混沌算法克服了遗传算法的早熟问题，其优化结果可靠且有效。综上，笔者采用并行遗传混沌算法，编制了计算软件。使用(VC+6.0)软件能方便的获取船舶航行性能分布函数中的每个最优设计变量，为船舶设计综合评估、船舶参数确定提供条件。关键

3、词：高速单体船；航行性能；结构力学特征；遗传混沌算法；综合优化1、 引言在工程领域中，我们可以利用经典的数学和相应的系统工程原则作为工具来完成分析、设计；用相关复杂数字来解释系统工程。但是，对于复杂的系统工程问题，尽管建立了分析、设计和制造模型，我们还将经常面对一些如相关环境等复杂因素和各种不确定性因素。这些因素涉及到一些限制指标和许多需要权衡的复杂条件，直接影响着设计的某些参数或进行评估的研究方法，这正是船舶设计的问题，也是复杂设计问题的典型代表。优化理论现在已广泛应用于科研，企业管理，军事事务等诸多领域。优化方法可分为常规优化方法和现代优化方法两种。前者通常会陷入局部最优，不适合多变量，多

4、约束，多目标函数的复杂优化。遗传算法，一种现代化的优化方法，具有很强的全局优化能力，可运用于船舶工程的复杂优化。但用遗传算法对多设计变量的函数进行优化时，削减的遗传算法一代很可能导致较大概率的早熟问题，而且会花费大量时间。因此，笔者提出了一种基于灵敏变量分段的分层并行遗传混沌算法，并成功地运用这一综合优化方法于船舶力学性能综合优化。2分层并行遗传混沌算法遗传算法是一种基于达尔文进化论和孟德尔遗传学理论的随机交换理论的新型优化方法。与常规优化方法相比，遗传算法具有许多优点，如计算速度快（快速性），全局搜索能力强（全局性），和良好的鲁棒性。大量研究和应用表明，遗传算法运用于复杂的综合优化效果显著。

5、然而，遗传算法是一个优胜劣汰的选择过程，所以会花大量的时间用来计算，在解决多设计变量优化时，效果不佳。并行遗传混沌算法使二者结合起来，加快运算速度。另一方面，遗传算法在处理含有多个目标函数和多约束条件的优化时，往往会出现早熟问题，带来特大种群规模和非常复杂的适应度函数。然而，运用分层并行遗传混沌算法，克服了这些缺点。基于灵敏变量分层并行遗传混沌算法，将以以下方式解决优化问题：1） 把每一个灵敏变量的值分成几个部分，并贯穿这些片段（灵敏变量是可变的，它的微小变化将引起目标函数值的重大变化）。然后用保留的变量值范围结合成片段，作为遗传算法的作业范围。最后，指定一个值到遗传算法的种群规模后，寻找恰当

6、的后代。例如，一个优化有12个设计变量，2个是灵敏变量。把其中一个变量值范围平分为6个片段，另一个平分为3个片段.在跨越这9个片段后由6 3我们可以得到18个片段组。最后，用其他的10个设计变量和每个片段组重组，便得到20个遗传算法运行范围。2） 在1）的基础上，挑选了几个用遗传算法算的好的结果，用混沌算法在前述优良结果的领域内进行再次运算，得到最优结果。21 遗传算法遗传算法是一种理想且高效的搜索(或优化)方法。考虑到遗传算法面对船舶航行性能综合优化这类非常复杂的多目标、多约束和多变量工程优化设计问题，带来的特大种群规模和特别复杂的适应度函数处理及计算时间过长和较大概率的早熟问题，本文采用的

7、优化计算分两步走：首先，采用遗传算法进行有限代数计算，并在遗传算法程序中嵌入辅助进化子函数，从每次计算后的种群中选择4个适应度相对较高的个体进行混沌搜索。接着，在遗传算法优化结果的小范围内采用二次载波混沌算法进一步完成最优值点的计算，作为改进后的遗传混沌算法，具体实现细节如下。第一步 遗传算法优化。遗传算法求解的组成为:编码方法在遗传算法中对每一个染色体进行适当的编码是必不可少的一步，算法中问题的构成以及遗传算子的使用都取决于编码方式。这些编码方法可以分为三大类：二进制编码方法、浮点数编码方法和符号编码方法Holland的简单遗传算法采用的是二进制编码，二进制编码理论分析方便，遗传操作与生物进

8、化相似，能在相同的范围内表示最多的模式，能够充分体现隐性的并行性。选择选择是实现进化的重要一环，许多学者提出了种种不同的选择方法，这里采用比较通用的轮盘赌法。交叉交叉运算使用单点交叉算子。在生物的自然进化过程中，两个同源染色体通过交配而重组，形成新的染色体，从而产生出新的个体或物种，在遗传算法中也使用交叉算子来产生新的个体或物种，在遗传算法中也使用交叉算子来产生新的个体。单点交叉算子是最常用和最基本的交叉操作算子。变异在生物的遗传和自然进化过程中，其细胞复制环节有可能会因为某些偶然因素的影响而产生一些复制差错，这样会导致生物的某些基因发生某些变异，从而产生新的染色体，表现出新的生物性状。在遗传

9、算法中也引入了变异算子来产生新的个体。基本位变异算子是最简单和最基本的变异操作算子。对于基本遗传算法中用二进制编码符号串所表示的个体，若需要进行变异操作的某一基因座上的原有基因值为0，则变异操作将该基因值变为1；反之亦然。适应度值在遗传算法中，适值是一个至关重要的概念。对于无约束优化问题，适值可取为该问题的优化目标值。但对于有约束的优化问题，适值就不一定是该问题的优化目标值。但其对应关系应是明显的。实际上，正是通过适值的合适选取，可将约束优化问题转化成无约束优化问题。我们知道，约束优化问题的一个可行解应满足约束条件，而最优解不但要满足约束条件，还必需使得目标函数值最大或最小。进化停机准则 如果

10、对遗传算法的搜索过程进行跟踪观察，我们会发现它将很快收敛于一个稳定的最优解。在实际优化过程中，为提高计算速度，可每隔10代判断一下迭代是否收敛，方法是计算出这代个体的适值的方差，当它小于给定值：(可取为1.0e-5)时，即认为找到最优解，可提前停止搜索。惩罚策略惩罚技术大概是用遗传算法解约束优化问题中最常用的技术。本质上它是通过惩罚不可行解将约束问题转化为无约束问题。在遗传算法中，惩罚技术用来在每代的种群中保持部分不可行解，使遗传算法可以从可行域和不可行域两边来达到最优解。惩罚策略的主要问题是如何设计一个惩罚函数，从而能有效地引导遗传搜索达到解空间的最好区域。2.2 混沌算法混沌算法求解的步骤

11、为:本文的混沌算法优化采用二次载波混沌算法。混沌优化是通过混沌变量实现的。混沌变量的产生有多种方法。这里选用应用较为广泛的Logistic映射，其方程如下： (2.1)式中：是控制参数，不难证明，当=4时，式(2.1)完全处于混沌状态，在数学上的解释就是通过式(2.1)的迭代可“随机” 产生(遍历)(0，1)之问除开不动点(即0.25，0.5，0.75)之外的所有数值。利用混沌对初值敏感的特性，赋予式(2.1)n个(0，1)之间的不同的初值(不动点除外)，即可得到n个轨迹不同的混沌变量。基于二次载波的混沌优化基本步骤为32： 初始化。k=0(k为首次载波迭代步数)，给予式(2.1)n个不同的初

12、值。(初始设计变量)，得到n个轨迹不同的混沌变量； 通过下式进行首次载波： (2.2)式中：和为常数，相当于平移和放大参数，其目的是将选定的n个混沌变量“映射”到优化变量中。显然，和分别为设计变量的下限与上限；如果自变量取值范围是开区间，则，其中和为很小的正数，可取为，且； 进行首次载波后的迭代搜索。先令表示(i =1，2，-，n)，相应的性能指标为，当前最优值点为，最优值为，并且，； 若经过若干次的搜索，都保持不变，则按下式进行二次载波 (2.3)式中：为二次载波后迭代步数，初始为0，为遍历区间很小的混沌变量，为调节常数； 用二次载波后的混沌变量继续迭代搜索，过程类似步骤； 如果满足终止条件

13、，则输出最优解及最优值，否则返回步骤。3 优化数学模型船舶力学性能综合优化包括两方面：航行性能和结构力学性能。笔者利用船舶的快速性、耐波性和操纵性作为分目标函数的航行性能；利用静态和动态性能的分目标函数的作为结构力学性能；这两个分目标函数构成总目标函数。不仅稳性、浮性和其他一些特性构成了约束条件，而且设计变量也构成约束条件。数学模型的详细描述如下：3.1 设计变量经综合分析现选择以下49个参数（包括35个关于舯剖面结构的参数）作为设计变量：船长L，船宽B，吃水T，纵向棱形系数Cp，舯横剖面系数，设计水线面系数，浮心纵向位置LCB(%LPP)，螺旋桨直径DP，盘面比AE/A0，螺距比P/DP，螺

14、旋桨转速n，设计航速Vs，半进水角，尾板相对浸湿面积At/Am等参数。同时，根据航行性能优化程序中的变量和总体设计已提前确定一些已知量，包括船舶的主尺度、舯剖面参数（甲板纵桁、船底纵桁、舷侧纵桁，支柱的位置和高度、中甲板位置）以及纵向构件的数量等。本文只好限制一部分横向构件和纵向构件的具体尺寸，以达到减少优化设计变量的目的。根据本文研究对象的结构特点，纵向构件的重量一般占整个舱段的35%-55%以上，进而部分增加了纵向构件的优化变量，减少横向构件和横舱壁的设计变量，并将桁材和肋板等具体参数设为已知量。同时为了减少结构分析和优化的规模，常把舯剖面上对总纵强度贡献较少的构件作为确定量，并在船舶的舯

15、横剖面处选择对总纵强度贡献较大的构件作为具体的优化变量,依次为：上甲板板厚、上甲板纵骨型号、舷侧板板厚、舷侧板纵骨型号、中甲板板厚、中甲板纵骨型号、舭列板厚度、舭列板纵骨型号、纵舱壁厚度、纵舱壁板纵骨型号、横舱壁厚度、横舱壁板纵骨型号、内底板厚度、内底板纵骨型号、船底板板厚、平板龙骨板厚、船底板纵骨型号、中桁材厚度、旁桁材厚度、上甲板纵骨间距、舷侧纵骨间距、中甲板纵骨间距、舱壁纵骨间距、舭板纵骨间距、内外底板纵骨间距。通常将它们定为设计变量，可采用下列向量形式表示为：X=x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x

16、15, x16, x17, x18, x19, x20, x21, x22, x23, x24, x25, x26, x27, x28, x29, x30, x31, x32, x33, x34, x35, x36, x37, x38, x39, x40, x41, x42, x43, x44, x45, x46, x47, x48, x493.2 目标函数P(X) 是总目标函数, PN(X) 是航行性能子目标函数, PS(X) 是结构特性子目标函数. 然而,P(X) =APPN(X) +BPPS(X) PN(X) =AP1Csp+ AP2Mv+ AP3Sv ；PS(X) =BP1J(X) +

17、BP2D(X)上式中： AP0, BP0, AP+BP=1; AP10, AP20, AP30, AP1+AP2+AP3=1; BP10, BP20, BP1+BP2=1. Csp, Mv 和 Sv 是在0,1内的相对标准化形式。快速性尺度C(Csp=PE/2/3Vs3oRH, 排水量, PE有效马力, o螺旋桨敞水系数, H船身效率, R相对旋转效率), 操纵性尺度M (M=pLVarL+pTVarT+pCVarC, VarL直线稳定系数, VarT旋转性系数, VarC流动变化系数, pL,pT and pC 重量数) 耐波性尺度S (S是俯仰角度和升沉幅度在0,1内的标准化形式). PS

18、(X) 是舯剖面纵向构件的单位长度重量式中：m-典型舱段横向强框架的数量； -材料密度；-典型舱段纵向构件的体积；-典型舱段第个横向强框架的体积；-典型舱段单个横舱壁的体积；-典型舱段的长度。3.3 约束条件3.3.1 等式约束条件(1)满足浮性条件，不变，LBTCB；(2)螺旋桨的推力与船舶航行时遭受的阻力平衡： 式中：为螺旋桨的桨轴数，t为推力减额分数，R为船之阻力。(3)设计状态下主机供给螺旋桨的转矩必须等于螺旋桨吸收的转矩： 式中：为主机功率，为桨轴转速，为水的质量密度。3.3.2. 不等式约束(1)设计变量的取值范围或上下限。(2)按海船稳性规范，初稳性高度GM，必须满足如下不等式约

19、束条件：GM h;(3)此外，设计的螺旋桨需满足空泡要求，按Keller公式计算; (4)回转无因次直径c。(5)结构的制约因素：构件许用应力如下表1所示。表.1 构件许用应力构件分类许用应力内外底板，舷侧板底板纵骨，舷侧纵骨横纵舱壁板，肋板甲板和平台结构平台纵骨，甲板纵骨，舷侧纵骨支柱注：为材料的屈服极限，4算例4.1优化计算 从数学模型可以看出，船舶动力学性能综合优化计算问题至少涉及49个设计变量，包括9个等式约束和5个不等式约束（若记入变量上下限约束，将有103个不等式约束）显然，这是一个非常复杂的工程优化问题。笔者采用基于灵敏变量的分层并行遗传混沌算法，编制了计算软件。 在此，以一艘排

20、水量为4200t的高速单体双桨船为例。该船设计变量范围如表2所示;笔者定义：AP=0.70, BP=0.30; AP1=0.50, AP2=0.30, AP3=0.20; 笔者运行了遗传算法8000代、并行遗传算法900代、混沌算法、并行混沌算法和并行遗传混沌算法。结果如表2所示表. 2 不同优化方法的计算结果项目混沌算法 遗传算法(8000 代)并行混沌算法 并行遗传算法(900 代)分层并行遗传混沌算法 (500 代)总目标函数值0.72740.75510.77750.82440.8469下限 /上限X1 (m)130/150141.44140.736139.104136.768134.4

21、64X2 (m)13/1514.46213.6914.49214.56414.468X3 (m)4 /5.54.88654.99454.99454.98254.933X10.4/0.540.407980.423240.405040.40980.4273X10.76/0. 860.85170.85780.83620.8450.8211X50.7/0.780.700640.741120.758080.770240.71936X-2/0-1.34-1.558-1.984-1.772-1.632X1 (m)2.4/33.1833.11223.1083.14283.4194X10.4/0.70.4729

22、0.41230.41950.580.4672X10.7/1.050.856450.95760.86030.83790.78085X1500/550531.3525.4510.1549.25509.2X1 (kn)30/3433.57233.80431.13232.63233.348X134/1210.35211.2885.6649.94411.312X140/0.180.133380.019980.060840.0360.08982X15(mm)7/97.8468.8548.2547.6567.942X16(mm)8/109.579.0488.0429.7649.75X171/91.1364.

23、3927.6082.4484.624X18(mm)44/7 /74.9635.5154.1174.5526.547X19(mm)12/1412.78212.4412.31413.17213.736X20(mm)10/1211.57611.02211.39210.29411.38X21(mm)12/1412.60212.36413.19413.49412.598X221/97.884.2882.8886.3685.896X23(mm)5/86.3387.2266.4466.6955.894X24(mm)5/75.1825.0245.745.0825.778X25(mm)5/65.4135.558

24、5.7725.9145.078X261/97.7762.7048.6323.048.872X27(mm)4/64.2485.2744.5864.8865.054X28(mm)7/87.2177.2587.0317.0157.984X29(mm)8/98.8018.7288.1668.28.894X30(mm)8/98.748.3268.0448.3718.968X311/94.5365.3445.4728.4962.144X32(mm)5/76.6726.4645.2326.6066.69X33(mm)5.5/76.42256.6136.20656.0616.6715X341/96.6565.

25、5843.8727.0087.864X35(mm)5.5/6.55.9266.2136.4225.8265.52X361/94.05658.7042.8963.064X37(mm)7.5/8.57.828.0987.9338.1717.629X381/91.1042.2886.0961.167.136X39(mm)9.5/1211.937511.8310.27510.287510.2825X40(mm)12/1512.49512.8112.21913.96513.494X411/91.9361.5924.884.924.28X42(mm)7/87.6777.6117.1497.9167.77X

26、43(mm)6/76.5636.296.9816.3586.306X44(m)0.7/10.96130.84460.95770.90730.8803X45(m)0.7/10.70750.89050.94840.96790.9256X46(m)1.25/1.3751.267881.374251.363751.373751.26625X47(m)0.5/0.70.62660.53040.60220.66120.5304X48(m)2/2.52.3482.45152.0322.35752.0125X49(m)0.43/0.540.450680.486210.510520.530870.47015排水

27、量 (t)4196.574191.264196.774185.454220.04有效推力TE (kN)1636.741772.681416.611645.441737.34阻力 (kN)1638.717791421.531647.861737.63Mp (kNm)776.218879.689669.523727.995860.467Md (kNm)777.146882.825671.845729.065860.614敞水效率o0.6801730.6599510.6613970.6865760.668634有效功率PE (kW)28299.430934.622764.727660.829807.

28、7主机功率 (kW)44299.649982.836695.342835.146535.2傅汝德数(Fr)0.4638520.4682240.4337360.45850.472558结果还包括其他20个参数，例如伴流分数、推力减额分数、初稳性高、摇周期、舵面积和湿面积等。 从结果中，我们可以看到，等式约束满意程度超过99.5 ，不等式约束则都全部满足。说明这种优化方法是可靠的。4.2 分析通过比较和分析这些不同优化方法的计算结果，得出以下4点结论：a、遗传算法进行8000代和并行遗传算法进行900代所得到的总目标函数值分别为0.7551和0.8244。前者低于后者8.4 ，表明遗传算法具有早熟

29、问题。b、混沌算法和并行混沌算法得到的总目标函数值分别为0.7274 和0.7775. 前者低于后者6.4%,表明并行混沌算法比混沌算法更有效。C、除b之外，我们还可以从并行混沌算法（基于灵敏变量分段）与混沌算法和并行混沌算法的总目标函数比较中，得到另一结论：基于灵敏变量片段的并行混沌算法更好。d、并行遗传混沌算法的总目标函数值为0.8469.比并行混沌算法高8.9%，比并行遗传算法高2.7%.这两组数据表明基于灵敏变量分段的并行遗传混沌算法是处理多目标，多约束，多变量的复杂工程优化的最好方法。5、结论本文提出了一种基于灵敏变量分段分层并行遗传混沌算法，并应用于高速单体船的动力学性能综合优化。

30、计算结果表明，这种解决方法的效率极高，为多目标，多约束，多变量工程优化设计提供了一个有效途径。大量的优化实例说明，该软件运行稳定，可靠。为高速单体船综合评估和船型参数的确定奠定了坚实的基础。参考文献：1 王健，王建华标准遗传算法的研究进展J华东船舶工业学院学报, 2000, (3):28342 刘如成，钉瑶君. 基于改进混沌算法的寻优策略J. 南京师大学报（自然科学版）, 2002,(4):3235.3 朱珉虎（译）. 高速方艉船型功率预报方法基于速度的再回归J. 中外船舶科技，2004,(4):3342.4 管义锋，俞铭华，赵建军，陈 练. 基于优化平台思想的舰船结构优化设计J. 中国造船,

31、 2004, (12)增刊:5256.5 吴剑国, 徐昌文. 船舶结构混合变量优化的遗传算法. 中国造船, 1998, (1):2126.6 中国船级社. 钢质海船入级与建造规范S. 北京:人民交通出版社,1996.7 俞铭华，谢祚水，吴建国.船舶舯剖面结构优化设计研究进展J. 华东船舶工业学院学报，1998, (9):3440.8 张火明，杨松林遗传算法在船舶航行性能优化中的惩罚策略J华东船舶工业学院学报, 2001, (4):4247.9 张火明，杨松林，朱仁庆，王志东船舶航行性能优化的模糊遗传算法J中国造船，2002，(3)：715.10 杨松林，陈淑玲，张火明. 船舶动力学性能及结构特

32、性的综合优化方法J. 中国造船, 2006,(2):3943.11 Yang Song-lin，Chen Shu-ling,Zhu Ren-qing,Wang Zhi-dong. Determination of Hull Form Parameters of Large-Scale Ship Based on Synthesis Optimization of Ship Navigation PerformanceJ. Journal of Ship Mechanics, V019 No6 Dec2005.12 杨松林, 奚炜, 张火明. 吃水受限的内河大型船舶航行性能综合优化方法研究J.

33、船海工程, 2006, (1):2330.On the overall optimization of navigation performance and structure characteristic for one kind of the high-speed shipsYang Song-lin，Li Mi, Chen Shu-ling(School of Naval Architecture and Ocean Engineering，Jiangsu University of Science and Technology，Zhenjiang Jiangsu 212003，Chin

34、a)AbstractBased on a kind of typical high speed monohull, the paper discusses an overall optimizations method of navigational performance and structure characteristic with the materials of some ships and some relevant experience expressions.At first, the mathematics model about overall optimization

35、of navigational performances for the high-speed monohull had built which based on a software (VC+6.0). The navigation performance goal function, the optimization variables and the constraint conditions make up of the mathematics model which includes ship fast, ship maneuverability and sea keeping. T

36、he author constructed an advanced algorithm that is composed of the chaos algorithm and the genetic algorithm. Then the software with friendly interface is programmed by C+ language. A large number of calculation results show that: the algorithm not only can conquer the prematurely problem of GA, bu

37、t also consume little calculation time and has higher efficiency. Secondly, the optimization computing of the midship section structure of the high-speed monohull with GA-chaos algorithm based on general specifications of naval ships (GJB 4000-2000) is introduced. The optimization model with some de

38、sign variables and extract some constraint conditions form General specifications for naval ships about longitudinal strength and buckling strength is established. The goal function is about minimal weight with the length of unit of the vertical component of midship section.At last, the data was gai

39、ned after this algorithm had applied to optimize overall the performance of navigational performance and structure characteristic by VC+，Quantities of computation results of different solving methods show that compared with other methods such as GA, chaos algorithm and their parallel algorithms, thi

40、s complex algorithm can overcome GAs premature convergence and is reliable and efficient.Above all, the software is programmed. Using this software can conveniently get the possibility distributing function of integrate optimum navigation performance of each design variables，and can provide the condition to integrate evaluation of the ship design project and the integrate decision-making of ship parameters. Keyword: high-speed monohull, navigational performance, structure characteristic, GA-chaos algorithm, integrate optimization11 / 11文档可自由编辑打印