10造船技术2008年第2期(总第282期)
内河大型船舶快速性能及结构力学特性综合优化方法研究
杨松林,杨昭龙,张火明
1
2
3
(1江苏科技大学,江苏镇江212003;2江苏大学,江苏镇江212013;3中国计量学院,浙江杭州310018)
传与分层分析思想所构造的一种基于敏感变量分段
提要建立了内河大型船舶快速性能及结构力学特性综合优化的数学模型,基于并行算法、遗传算法和分层思想,构造了一种基于敏感变量分段的分层并行遗传复合算法,并将其应用于求解此类综合优化计算问题。对于内河大型船舶快速性能及结构力学特性综合优化问题,进行了遗传算法或并行遗传算法及其分层并行遗传算法的大量优化计算,结果表明:该复合算法不但能有效地克服遗传算法的早熟问题,而且计算可靠、效率高;为内河大型船舶设计方案的综合评估及综合优化船型设计准备了前提条件。
关键词内河船快速性结构强度结构稳定性
最优设计遗传算法
中图分类号U662文献标识码B
的分层并行遗传复合算法,并介绍将其应用于内河
大型船舶快速性能及结构力学特性综合优化计算中的研究结果。
2分层并行遗传算法
2.1分层并行遗传算法概要
遗传算法是基于达尔文的生物遗传及进化原理和遗传学的随机交换理论发展起来的一种新的寻优方法;遗传算法与一般的优化方法相比,它具有全局性、快速性、适应性和鲁棒性等优点。近十多年的研究和应用实践表明,该算法具有解决复杂系统优化问题的非凡能力。但是,这种算法的运行过程及实现方法本质上仍是串行的。因此,在求解一些实际优化问题时,特别是非常复杂的系统优化问题时,它一般具有很大的群体规模,从而使得算法的遗传进化运算过程进展缓慢,导致计算耗时长、效率相对较低。应用实践表明,并行算法、遗传算法和分层分析思想相结合,形成一种合理寻优构造有望解决上述问题。考虑到遗传算法用于求解内河大型船舶快速性能及结构力学特性综合优化这类非常复杂的多目标、多约束和多变量工程优化设计问题所带来的特大种群规模和特别复杂的适应度函数,以及处理和计算时间过长或较大概率的早熟问题,本文采用的基于敏感变量分段的分层并行遗传算法优化计算分两步走:首先,采用遗传算法进行10个以上种群有限代数计算,寻优区间为敏感参数可行域各分段全面交叉与剩余变量可行域的组合;第二步是在遗传算法优化结果(6至12组最优结果)的小范围内,再采用遗传算法进一步完成最优点的计算。例如,某优化问题共涉及6个变量,其中2个敏感变量,而这两个敏感变量中,一个可行域等分为5个段,另一个分为6个段。由于两个敏感变量可行域分段全面交叉组合为30,因此,可得新的寻优区间等于30。一般种群数选定为新寻优区间数的整数倍,并且将其1引言
随着现代化生产的高速发展和科学技术的进步,特别是计算机软硬件技术的飞速发展,使得最优化理论及其应用日益受到科技人员的广泛重视。目前,它已渗透到科研与设计、生产与管理、商业与军事等各领域。优化方法通常可分为常规优化方法和现代优化方法两类。常规优化方法常常陷入局部最优而不能自拔,全局寻优能力差,一般不适合于多变量、多约束和多目标的复杂优化问题的求解;而现代优化方法的遗传算法及其复合方法,具有很强的全局寻优能力。近年来,遗传算法或模糊遗传等算法已被用于求解复杂船舶工程中的优化问题。但是,对于兼顾船舶主尺度受限条件下的船舶航行力学性能及结构力学特性综合优化这种必须合理处理大量错综复杂关系的工程优化问题,如果也采用常规遗传算法或模糊遗传算法完成寻优计算,那么,优化计算耗时将会很长;而如果将优化计算的代数降下来,那么优化计算过程中,早熟现象概率较高的问题又不可避免。因此,本文尝试论述将基于并行算法、遗
作者简介:杨松林(1956-),男,教授。
杨松林,等:内河大型船舶快速性能及结构力学特性综合优化方法研究11
平均分配给各区间。2.2遗传算法实现细节
(1)初始种群的确定。运用遗传算法进行优化计算时,如何确定初始种群是个很重要的问题,它和编码方案息息相关。在本文给出的算法中,假定一代种群的确定包括每个个体的适值、相对适值、累积适值、对应每维变量下限、对应每维变量上限和基因组共6项。
(2)编码。一般遗传算法中,染色体编码都采用0、1二进制串。考虑到优化数学模型变量较多,为节约编码、解码过程机时消耗,提高优化进程的速度,本算法的染色体采用浮点数据编码。
(3)适值(Fitness)。对于无约束优化问题,适值可取为该问题的优化目标值。但对于有约束的优化问题,适值就不一定是该问题的优化目标值,而其对应关系应是明显的。实际上,正是通过适值的合适选取,可将约束优化问题转化成无约束优化问题。约束优化问题的一个可行解应满足约束条件,而最优解不但要满足约束条件,还必须使得目标函数值最大或最小。这样, 染色体(Chromosome)!的 适值!应包括两方面信息,一是该 染色体!满足约束程度,一是该 染色体!优化目标函数值。
(4)选择。选择是实现进化的重要一环,许多学者提出了种种不同的选择方法,这里采用比较通用的轮盘赌法。
(5)交叉。遗传系统的一个重要的参数是杂交概率pc,本算法取为0.8,此概率给出预计要进行杂交的染色体个数是pc=0.85。
(6)变异。遗传系统的另一个重要系统参数是变异概率pm,它给出预计的变异位数,对杂交后的当前群体中的每个染色体和染色体中的每个位,即产生在区间[0,1]里的一随机浮点数r;如果r(7)进化停机准则。如果对遗传算法的搜索过程进行跟踪观察,我们会发现它将很快收敛于一稳定的最优解。这样,我们可以设定一个最大叠代次数,例如300,并认为进化到300代即可停止演化,亦即第300代中最好的染色体即是问题的最优解。在实际优化过程中,为提高计算速度,可每隔10代判断一下叠代是否收敛。方法是,计算出该代个体的适值的方差,当它小于给定值s(可取为1.0e5)时,即认为找到最优解,可提前停止搜索。3综合优化数学模型
内河大型船舶快速性能及结构力学特性综合优化,包括两大部分:即船舶快速性能和结构力学特性
的优化。
本文将快速性能指标和结构力学特性指标的加权和,作为综合优化目标函数,而将稳性和浮性等其他一些性能及船型主要要素(参数)的,作为优化约束条件,构成一种内河大型船舶性能和结构力学特性综合优化数学模型。具体描述如下。3.1设计变量
内河大型船舶快速性能及结构力学特性综合优化计算需要兼顾的关系错综复杂,涉及的船型参数等因素也很多。经分类与分层综合权衡和比较分析,选定以下12个参数作为主要设计变量:船长L,船宽B,吃水T,方形系数Cb,中横剖面系数Cm,设计水线面系数Cwp,浮心纵向位置Lcb,螺旋桨直径Dp,螺旋桨盘面比Ae/Ao,螺旋桨螺距比P/Dp,螺旋桨转速n,设计航速Vs。为描述简单起见,将设计变量用一向量X表示,即:
X={L,B,T,Cb,Cm,Cwp,Lcb,Dp,Ae/Ao,P/
Dp,n,Vs,其他变量}。3.2目标函数
设P(X)为总目标函数,Pp(X)为船舶快速性能优化目标函数,Ps(X)为船舶结构力学特性优化目标函数,那么有
P(X)=ApPp(X)+BpPs(X);Pp(X)=Csp;
Ps(X)=Bp1J(X)+Bp2D(X),式中AP∀0,B
p
T
∀0,Ap+Bp=1;
Bp1∀0,Bp2∀0,Bp1+Bp2=1;
Csp是压缩至(0-1)之间的船舶快速性指标C;
J(X)和D(X)分别是压缩至(0-1)之间的船舶结构静力特性和船舶结构动力特性指标。C=Pe/Vs 0 R H,
式中是排水量,Pe为有效功率, 0为螺旋桨敞水效率, H为船身效率, R为相对旋转效率。3.3约束条件3.3.1等式约束
(1)满足浮性条件,=LBTCb;
(2)螺旋桨的有效推力与船舶航行时遭受的阻力相等;2/3
3
12造船技术2008年第2期(总第282期)
(3)设计状态下主机供给螺旋桨的转矩必须等于螺旋桨吸收的转矩。
3.3.2不等式约束
(1)12个设计变量取值范围;
(2)初稳性高度条件:初稳性高GM>a;(3)极限横摇幅值衡准条件:!a#b∃;
(4)回转性衡准条件:相对回转直径Ds%非常复杂的工程优化问题。笔者采用如上所述基于敏感变量分段的分层并行遗传算法,并为上述数学模型编制了计算软件。由于篇幅所限,在此仅以双桨、排水量4.3万t大型中速内河运输船为计算和比较分析实例,该船优化设计变量范围如表1。
选取Ap=0.75,Bp=0.25。J(x)=(10000-0.28LBTCb)/3000。在一台P4(主频1.5GHz,内存512MB,硬盘40GB)计算机上完成了基于6000代遗传算法优化计算、基于1000代遗传算法的并行遗传算法的优化计算、基于敏感变量分段分层并行遗传算法的优化计算(第一层和第二层遗传计算分别为100代和300代),计算结果见表2。
Cwp
Lcb-2.00
Dp5.76.5
Ae/Ao0.50.9
P/Dp0.51.1
n129160
Vs2225
0.9624
0.9014
0.883
0.037
4优化计算及结果分析
4.1优化计算
内河大型船舶快速性能及结构力学特性综合优化是一个设计变量多、约束条件多和多目标权衡的
L
下限上限
219240
B2831
T8..7
Cb0.590.696
Cm0.940.99
表1设计变量范围
0.740.78
表2不同优化方法计算结果
各参数名或符号总目标函数值
L(m)B(m)T(m)CbCmCwpLcbDp(m)Ae/AoP/DpN(r/min)Vs(kn)1.0/Csp桨轴数桨类型桨叶数排水量(t)有效推力船体阻力螺旋桨转矩
主机供给螺旋桨转矩
敞水效率 0有效功率主机功率傅汝德数
遗传算法的结果(遗传计算6000代)
0.85074226.0429.029.1520.69550.94810.7418-1.4906.1320.51400.7882141.322.4769.882AU型443003.691654.681654.811847.781847.930.673519122.9827660.290.2455
并行遗传算法的结果(遗传计算1000代)
0.86350221.0029.069.4360.600.940.7436-1.8486.3100.50200.7936133.822.3771.432AU型443001.611618.681618.541879.171879.010.682918623.9226714.920.2473
分层并行遗传算法的结果(遗传计算100和300代)
0.87919219.728.719.5120.69580.950.7409-1.9126.3300.500.8120131.322.5072.362AU型442997.961627.211627.141920.861920.780.685318828.3326815.560.2494
[下转第46页]
46造船技术2008年第2期(总第282期)
入ACCESS内的一个数据库表内,已经达到整艘船生产设计数据同时能够检索的目的。如果多艘船所有分段零件需同时导入一个数据库的相同数据表内,则应保证两艘船分段名和零件代码不一样,否则系统将覆盖上一次导入的相同数据。多艘船生产设计零件表内容导入ACCESS内的数据表,应按上述步骤一一进行导入不同的数据表。在分段数据导入ACCESS内的数据表时,对导入数据进行检索分类是按照数据存储规则进行分类的。
据库,集合所有零散分段的加工与制造内容,方便管理人员是编制工作计划,使管理人员在信息处理、物料需求、人员需求等方面提前做足准备,不至于在某一工位出现人员窝工、物料等待等不利于流水生产的现象。
5参考文献
1翁德伟.造船成组技术.上海:上海交通大学出版社2中国船舶工业总公司造船生产设计指导组编.造船生产
设计.北京:国防工业出版社
3姚巍.VisualBasic数据库开发及工程实例.北京:人民邮
电出版社
4关清玉,陈宁.船舶生产设计日程管理系统研究.船舶,
2005,(5):24
4总结
本文对船体内场业务据流程进行了分析优化,
绘制出了船体内场业务流程图,同时设计了内场生产管理辅助软件,该软件能够将零件表内容转入数
[上接第12页]
由上表易得,三个等式约束的满足程度高达99.99%以上;同时,不等式约束均获得100%的满足,说明上述优化方法计算可靠性良好。4.2结果分析
根据表2,经过比较易得结论如下:
(1)基于遗传算法6000代的优化计算的总目标函数值为0.85074,较基于1000代遗传算法的(12个种群)并行遗传算法的优化计算的目标函数值0.86350低约1.48%,呈现出明显的早熟特征;可见并行算法明显占优。
(2)基于分层并行遗传算法优化计算的目标函数值等于0.87919,比基于1000代并行遗传算法的高1.82%、比6000代遗传算法的高3.34%。可见,分层并行遗传算法构造可以非常有效地提高求解复杂三多(多目标、多约束和多变量)工程优化问题的计算效率。
法不但能有效地克服遗传算法的早熟问题,而且计算效率高。为有效解决复杂的三多(多目标、多约束和多变量)工程优化设计问题,提供了一种可借鉴的途径。
6参考文献
1[美]Z.米凯利维茨,著.周家驹,译.演化程序&&&遗传算
法与数据编码的结合.北京:科学出版社,2000.
2刘勇,康立山,陈毓屏.非数值并行算法第二册遗传算
法.北京:科学出版社,1998.
3王健,王建华.标准遗传算法的研究进展.华东船舶工业
学院学报,2000,(3):28
4YangSonglin,ZhangHuoming,ZhuRenqing,etal.The
SynthesisOptimizationofShipNavigationPerformanceBasedonFuzzyGeneticAlgorithm.ThirdConferenceforNewShipandMarineTechnology.May2123,2002.Kobe,Japan.
5YangSonglin,ZhouYulong.Adesignmethodofhull
formparametersforlargescaleshipsoptimalintegrativenavigationperformance.2005,(6):
JournalofShipMechanics,
5结束语
综上所述,本文提出了一种基于敏感变量分段
的分层并行遗传算法,并应用于内河大型船舶快速性能及结构力学特性综合优化计算,结果表明:该算