高速单体船航行性能与结构特性综合优化方法研究杨松林，李密, 陈淑玲摘 要本文以典型高速单体船为基础，根据一些实船的资料和相关经验公式并结合高速单体船的结构特点，探讨一种将船舶航行性能与结构力学特征两者综合优化的方法。首先依次建立了基于 (VC+6.0)软件的高速单体船快速性优化、操纵性优化、耐波性优化的数学模型，并包括了这三个方面的优化变量、目标函数和约束条件。在选择算法方面，本文将遗传算法与混沌算法复合构造了一种遗传混沌算法，并利用C+语言对其编制了友好界面。大量的计算数据表明：遗传混沌算法不仅能克服遗传算法的不足，而且能减少计算、提高效率。其次，基于舰船通用设计规范（GJB 4000-2000）对高速单体船舯横剖面进行结构优化设计。将舯剖面上影响其剖面模数的主要板厚、骨材型号等作为设计变量，以舯剖面纵向构件的单位长度重量最轻为目标函数，并根据相关规范提取相应的约束条件，采用遗传混沌算法对舯剖面的主要构件进行优化计算。最后，把船舶航行性能与结构力学特征综合优化的遗传混沌算法结果，与其它优化方法（例如，遗传算法、混沌算法）的优化结果进行比较。经比较可知，遗传混沌算法克服了遗传算法的早熟问题，其优化结果可靠且有效。综上，笔者采用并行遗传混沌算法，编制了计算软件。使用(VC+6.0)软件能方便的获取船舶航行性能分布函数中的每个最优设计变量，为船舶设计综合评估、船舶参数确定提供条件。关键词：高速单体船；航行性能；结构力学特征；遗传混沌算法；综合优化1、 引言在工程领域中，我们可以利用经典的数学和相应的系统工程原则作为工具来完成分析、设计；用相关复杂数字来解释系统工程。但是，对于复杂的系统工程问题，尽管建立了分析、设计和制造模型，我们还将经常面对一些如相关环境等复杂因素和各种不确定性因素。这些因素涉及到一些限制指标和许多需要权衡的复杂条件，直接影响着设计的某些参数或进行评估的研究方法，这正是船舶设计的问题，也是复杂设计问题的典型代表。优化理论现在已广泛应用于科研，企业管理，军事事务等诸多领域。优化方法可分为常规优化方法和现代优化方法两种。前者通常会陷入局部最优，不适合多变量，多约束，多目标函数的复杂优化。遗传算法，一种现代化的优化方法，具有很强的全局优化能力，可运用于船舶工程的复杂优化。但用遗传算法对多设计变量的函数进行优化时，削减的遗传算法一代很可能导致较大概率的早熟问题，而且会花费大量时间。因此，笔者提出了一种基于灵敏变量分段的分层并行遗传混沌算法，并成功地运用这一综合优化方法于船舶力学性能综合优化。2分层并行遗传混沌算法遗传算法是一种基于达尔文进化论和孟德尔遗传学理论的随机交换理论的新型优化方法。与常规优化方法相比，遗传算法具有许多优点，如计算速度快（快速性），全局搜索能力强（全局性），和良好的鲁棒性。大量研究和应用表明，遗传算法运用于复杂的综合优化效果显著。然而，遗传算法是一个优胜劣汰的选择过程，所以会花大量的时间用来计算，在解决多设计变量优化时，效果不佳。并行遗传混沌算法使二者结合起来，加快运算速度。另一方面，遗传算法在处理含有多个目标函数和多约束条件的优化时，往往会出现早熟问题，带来特大种群规模和非常复杂的适应度函数。然而，运用分层并行遗传混沌算法，克服了这些缺点。基于灵敏变量分层并行遗传混沌算法，将以以下方式解决优化问题：1） 把每一个灵敏变量的值分成几个部分，并贯穿这些片段（灵敏变量是可变的，它的微小变化将引起目标函数值的重大变化）。然后用保留的变量值范围结合成片段，作为遗传算法的作业范围。最后，指定一个值到遗传算法的种群规模后，寻找恰当的后代。例如，一个优化有12个设计变量，2个是灵敏变量。把其中一个变量值范围平分为6个片段，另一个平分为3个片段.在跨越这9个片段后由6 3我们可以得到18个片段组。最后，用其他的10个设计变量和每个片段组重组，便得到20个遗传算法运行范围。2） 在1）的基础上，挑选了几个用遗传算法算的好的结果，用混沌算法在前述优良结果的领域内进行再次运算，得到最优结果。21 遗传算法遗传算法是一种理想且高效的搜索(或优化)方法。考虑到遗传算法面对船舶航行性能综合优化这类非常复杂的多目标、多约束和多变量工程优化设计问题，带来的特大种群规模和特别复杂的适应度函数处理及计算时间过长和较大概率的早熟问题，本文采用的优化计算分两步走：首先，采用遗传算法进行有限代数计算，并在遗传算法程序中嵌入辅助进化子函数，从每次计算后的种群中选择4个适应度相对较高的个体进行混沌搜索。接着，在遗传算法优化结果的小范围内采用二次载波混沌算法进一步完成最优值点的计算，作为改进后的遗传混沌算法，具体实现细节如下。第一步 遗传算法优化。遗传算法求解的组成为:编码方法在遗传算法中对每一个染色体进行适当的编码是必不可少的一步，算法中问题的构成以及遗传算子的使用都取决于编码方式。这些编码方法可以分为三大类：二进制编码方法、浮点数编码方法和符号编码方法Holland的简单遗传算法采用的是二进制编码，二进制编码理论分析方便，遗传操作与生物进化相似，能在相同的范围内表示最多的模式，能够充分体现隐性的并行性。选择选择是实现进化的重要一环，许多学者提出了种种不同的选择方法，这里采用比较通用的轮盘赌法。交叉交叉运算使用单点交叉算子。在生物的自然进化过程中，两个同源染色体通过交配而重组，形成新的染色体，从而产生出新的个体或物种，在遗传算法中也使用交叉算子来产生新的个体或物种，在遗传算法中也使用交叉算子来产生新的个体。单点交叉算子是最常用和最基本的交叉操作算子。变异在生物的遗传和自然进化过程中，其细胞复制环节有可能会因为某些偶然因素的影响而产生一些复制差错，这样会导致生物的某些基因发生某些变异，从而产生新的染色体，表现出新的生物性状。在遗传算法中也引入了变异算子来产生新的个体。基本位变异算子是最简单和最基本的变异操作算子。对于基本遗传算法中用二进制编码符号串所表示的个体，若需要进行变异操作的某一基因座上的原有基因值为0，则变异操作将该基因值变为1；反之亦然。适应度值在遗传算法中，适值是一个至关重要的概念。对于无约束优化问题，适值可取为该问题的优化目标值。但对于有约束的优化问题，适值就不一定是该问题的优化目标值。但其对应关系应是明显的。实际上，正是通过适值的合适选取，可将约束优化问题转化成无约束优化问题。我们知道，约束优化问题的一个可行解应满足约束条件，而最优解不但要满足约束条件，还必需使得目标函数值最大或最小。进化停机准则 如果对遗传算法的搜索过程进行跟踪观察，我们会发现它将很快收敛于一个稳定的最优解。在实际优化过程中，为提高计算速度，可每隔10代判断一下迭代是否收敛，方法是计算出这代个体的适值的方差，当它小于给定值：(可取为1.0e-5)时，即认为找到最优解，可提前停止搜索。惩罚策略惩罚技术大概是用遗传算法解约束优化问题中最常用的技术。本质上它是通过惩罚不可行解将约束问题转化为无约束问题。在遗传算法中，惩罚技术用来在每代的种群中保持部分不可行解，使遗传算法可以从可行域和不可行域两边来达到最优解。惩罚策略的主要问题是如何设计一个惩罚函数，从而能有效地引导遗传搜索达到解空间的最好区域。2.2 混沌算法混沌算法求解的步骤为:本文的混沌算法优化采用二次载波混沌算法。混沌优化是通过混沌变量实现的。混沌变量的产生有多种方法。这里选用应用较为广泛的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个不同的初值。(初始设计变量)，得到n个轨迹不同的混沌变量； 通过下式进行首次载波： (2.2)式中：和为常数，相当于平移和放大参数，其目的是将选定的n个混沌变量“映射”到优化变量中。显然，和分别为设计变量的下限与上限；如果自变量取值范围是开区间，则，其中和为很小的正数，可取为，且； 进行首次载波后的迭代搜索。先令表示(i =1，2，-，n)，相应的性能指标为，当前最优值点为，最优值为，并且，； 若经过若干次的搜索，都保持不变，则按下式进行二次载波 (2.3)式中：为二次载波后迭代步数，初始为0，为遍历区间很小的混沌变量，为调节常数； 用二次载波后的混沌变量继续迭代搜索，过程类似步骤； 如果满足终止条件，则输出最优解及最优值，否则返回步骤。3 优化数学模型船舶力学性能综合优化包括两方面：航行性能和结构力学性能。笔者利用船舶的快速性、耐波性和操纵性作为分目标函数的航行性能；利用静态和动态性能的分目标函数的作为结构力学性能；这两个分目标函数构成总目标函数。不仅稳性、浮性和其他一些特性构成了约束条件，而且设计变量也构成约束条件。数学模型的详细描述如下：3.1 设计变量经综合分析现选择以下49个参数（包括35个关于舯剖面结构的参数）作为设计变量：船长L，船宽B，吃水T，纵向棱形系数Cp，舯横剖面系数，设计水线面系数，浮心纵向位置LCB(%LPP)，螺旋桨直径DP，盘面比AE/A0，螺距比P/DP，螺旋桨转速n，设计航速Vs，半进水角，尾板相对浸湿面积At/Am等参数。同时，根据航行性能优化程序中的变量和总体设计已提前确定一些已知量，包括船舶的主尺度、舯剖面参数（甲板纵桁、船底纵桁、舷侧纵桁，支柱的位置和高度、中甲板位置）以及纵向构件的数量等。本文只好限制一部分横向构件和纵向构件的具体尺寸，以达到减少优化设计变量的目的。根据本文研究对象的结构特点，纵向构件的重量一般占整个舱段的35%-55%以上，进而部分增加了纵向构件的优化变量，减少横向构件和横舱壁的设计变量，并将桁材和肋板等具体参数设为已知量。同时为了减少结构分析和优化的规模，常把舯剖面上对总纵强度贡献较少的构件作为确定量，并在船舶的舯横剖面处选择对总纵强度贡献较大的构件作为具体的优化变量,依次为：上甲板板厚、上甲板纵骨型号、舷侧板板厚、舷侧板纵骨型号、中甲板板厚、中甲板纵骨型号、舭列板厚度、舭列板纵骨型号、纵舱壁厚度、纵舱壁板纵骨型号、横舱壁厚度、横舱壁板纵骨型号、内底板厚度、内底板纵骨型号、船底板板厚、平板龙骨板厚、船底板纵骨型号、中桁材厚度、旁桁材厚度、上甲板纵骨间距、舷侧纵骨间距、中甲板纵骨间距、舱壁纵骨间距、舭板纵骨间距、内外底板纵骨间距。通常将它们定为设计变量，可采用下列向量形式表示为：X=x1, x2, x3, x4, x5, x6, x7, x8, x9, x10, x11, x12, x13, x14, x15, 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) 是结构特性子