2. 招商局重工(江苏)有限公司,江苏 南通 226000;
3. 华中科技大学,湖北 武汉 423201
2. China Merchants Heavy Industry (Jiangsu) Co., Ltd., Nantong 226000, China;
3. Huazhong University of Science and Technology, Wuhan 423201, China
螺旋桨设计及优化一直是研究船舶快速化、经济性的紧迫任务之一。但由于螺旋桨参数众多、建模繁琐、水动力计算周期长等问题[1],长期限制船舶设计初期对螺旋桨的高效优化设计。建立一套包括螺旋桨几何变形和重构、水动力性能预测近似模型以及多目标优化的螺旋桨构型优化设计方法[2],对促进螺旋桨高效优化设计具有重要意义。
螺旋桨的快速建模、高效优化始终是国内外学者关注的热点研究。Ramin等[3]采用Non-Gradinet-based Algorithm和梯度法分别对弦长、厚度以及螺距沿半径的分布进行数值优化;Florian等[4]基于B-spline对螺距、拱度、厚度、弦长、曲率和倾角沿半径分布进行拟合,并以最大效率和最小脉动压力为目标进行了螺旋桨优化。叶礼裕[5]结合定常面元法和克里格近似理论建立螺旋桨参数与性能特征的预测模型,并利用粒子群算法对螺旋桨参数进行了多目标优化设计。王超[6]用试验设计方法为神经网络近似模型的建立提供足量且高质量的训练样本,利用遗传算法可以获得更广解集范围更广,而且耗时间大幅度缩减。马丹萍[2]借助于ISIGHT优化平台,集成螺旋桨的FFD方法和CFD水动力数值评估,提出了一套从构型-数值求解-优化的螺旋桨设计优化方法。
针对目前螺旋桨研究中存在的螺旋桨模型构造周期长、水动力性能评估耗时多的问题,本文结合基于NURBS基函数的自由变形算法(NURBS-based free-from deformation,NFFD)、多输出高斯近似模型以及NSGA-Ⅱ多目标优化算法,构建一套包括螺旋桨变形重构—水动力性能快速预测—多目标优化的螺旋桨高效自动优化方法,即以NFFD技术实现螺旋桨三维模型的变形与重构,采用有限元数值仿真得出样本螺旋桨的水动力性能,基于多输出高斯近似理论建立螺旋桨水动力性能预测近似模型,结合近似模型和NSGA-Ⅱ算法对KP505桨进行多目标优化设计,验证该方法的可行性。
1 NFFD方法及螺旋桨变形控制基于Sederberg等[7]提出的自由曲面变形(free form deformation,FFD)方法,研究人员将FFD技术扩展为很多不同形式[8-9],基本包括:构建参数控制体;将变形物嵌入参数控制体;参数控制体变形;计算变形物形变后的几何坐标[10]。基于FFD方法,Lamousin等[8]提出NFFD算法,将顶点加权与非均匀框架结合,引入权因子和节点作为变量,增加了变形操作的自由度,从而达到更好的变形。马丹萍[2]、张人会[10]、李冬琴[11]、WU[12]、马明生[13]、唐静[14]已将NFFD技术应用于船型与螺旋桨变形控制与重构,并在船型与螺旋桨优化设计方面取得了丰硕成果。
螺旋桨上任意一点绝对坐标可以用局部坐标系
Q=Q0+tT+uU+vV。 | (1) |
式中:
Pi,j,k=Q0+i/l⋅T+j/m⋅U+k/n⋅V。 | (2) |
式中:l,m,n分别表示局部坐标在
更进一步,螺旋桨上任一点绝对坐标可以用局部坐标表示为:
Q(t,u,v)=l∑i=0m∑j=0n∑k=0Pi,j,kBil(t)Bjm(u)Bkn(v)。 | (3) |
式中:(v,u,t)为局部坐标;Pi,j,k为参数化控制节点;
Bil(t)=l!i!(l−i)!ti(l−t)(l−i)。 | (4) |
NFFD方法利用NURBS基函数替换Bernstein函数,螺旋桨任意一点Q坐标可表示为:
Q(t,u,v)=l∑i=0m∑j=0n∑k=0Pi,j,kWi,j,kBil(t)Bjm(u)Bil(v)l∑i=0m∑j=0n∑k=0Wi,j,kBil(t)Bjm(u)Bil(v)。 | (5) |
式中,Wi,j,k为参数化控制节点Pi,j,k的权重因子。
在式(5)基础上采用deBoor-Cox递推定义[15]和Deboor算法[16]推导得出NURBS体三次矩阵表示方法:
Q(t,u,v)=ˆVMk(ˆTMiPi,j,kWi,j,kMTjˆUT)ˆVMk(ˆTMiWi,j,kMTjˆUT)。 | (6) |
式中:
最后将螺旋桨嵌入以上参数控制体,并对采样点进行归一化处理,利用牛顿法逆向搜索非线性方程组的最优解,即可得到螺旋桨采样点局部坐标。
在建立螺旋桨与参数体的相互关系后,保持每个采样点在参数体中相应的局部坐标(t,u,v)不变,通过改变参数化控制节点Pi,j,k或其权重因子Wi,j,k即可得到新的带权控制点矩阵,代入式(6)中即可求解得出螺旋桨变形后的局部坐标。
以KP505桨为例,利用Matlab实现了螺旋桨几何形变与重构,初始输入螺旋桨表面离散点集的坐标值,并选取0.2R,0.4R,0.6R,0.7R和0.9R处的叶宽、厚度、螺距比、曲率、倾角、最大厚度位置共30个参数作为螺旋桨形状控制设计变量,输出螺旋桨变形后的表面离散点集的坐标值。KP505桨形状控制效果和网格模型如图1所示。
![]() |
图 1 KP505桨形状控制效果和网格模型 Fig. 1 KP505 propeller shape control effect and mesh model |
近似模型是基于不同模型函数对设计变量和设计目标组成的样本数据进行近似拟合,具有计算成本低、函数关系明确等特点[17]。常用的近似模型方法包括:响应面法(RSM)、克里格法(Kringing)和人工神经网络技术[18]。高斯近似是基于Kringing和贝叶斯定理的一种机器学习方法,对少量样本的高阶非线性问题具有良好的适应性[19-22]。
给定训练样本:
Slearn={(xi,yi),i=1,⋯,Nlearn}, |
其中,任意一组样本
ˆyi≈f(xi), |
高斯近似模型为了保证预测值
ˆyi=f(xi)+ε, |
其中,
含有高斯噪声的预测值
cov(ˆy)=Kn(X,X)+σ2nI。 |
式中:
多输出高斯近似模型协方差矩阵K表示为:
Kn=Kf⊗Kx, |
式中,
训练样本的输出向量
[yf(xtest)]∼N([00],[K(X,X)+σ2nIk(X,xtest)k(xtest,X)k(xtest,xtest)])。 |
式中,
f(xtest)|X,y,xtest∼N(ˉf(xtest),σ2(xtest)), |
其中预测均值
ˉf(xtest)=k(X,xtest)T[K(X,X)+σ2nI]−1y, |
σ2(xtest)=k(xtest,xtest)−k(X,xtest)T×[K(X,X)+σ2nI]−1k(X,xtest)。 |
本研究中的输入向量
{X={xai,xbi,xci,xdi,xei,xfi⋯⋯},Y={KT,10KQ,η0,Pmin}。 |
式中:a~f分别表示某一剖面处的叶宽、厚度、螺距比、曲率、倾角、最大厚度位置;i=0.2R,0.4R,0.6R,0.7R和0.9R,KT,10KQ,η0以及Pmin分别表示螺旋桨推力系数、扭矩系数、效率以及最小空化压力。
2.2 敞水性能数值仿真以STAR-CCM+13.0流体动力学软件对KP505桨的敞水性能进行仿真计算。利用NFFD建立的原始螺旋桨曲面生成实体模型,构建非结构化六面体网格,采用k-ωSST模型,并基于MRF方法求解雷诺平均Navier-Stokes方程计算不同进速下的螺旋桨性能。
以进速系数J=0.893为例,在计算单元达到26000000时,螺旋桨的端涡和根涡现象如图2所示,即在此时实现了网格收敛。
![]() |
图 2 高雷诺数下螺旋桨的端涡和根涡现象 Fig. 2 End vortex and root vortex phenomenon of propeller at high Reynolds number |
图3为进速系数J=0.2~1.1时,螺旋桨推力系统、仿真计算结果与试验结果对比,仿真与试验所得推力系数KT、扭矩系数10KQ、效率η0的均方差分别为1.63%,2.41%,0.63%,说明以该软件和网格数量计算的螺旋桨敞水性能精度可信。
![]() |
图 3 螺旋桨仿真计算结果与试验结果对比 Fig. 3 Comparison between the simulation calculation results of the propeller and the test results |
为给高斯近似模型提供初始样本,以0.2R,0.4R,0.6R,0.7R和0.9R处的叶宽、厚度、螺距、曲率、倾角、最大厚度位置共30个参数作为螺旋桨形状控制设计变量,采用NFFD技术生成100个螺旋桨计算样本,并计算进速系数J=0.893时KT,10KQ,η0以及Pmin。
2.3 近似模型的验证对100个初始样本采用交叉验证方法对训练后的多输出高斯代理模型进行评估。通过交叉验证得到的100个测试样本预测误差如图4所示。
![]() |
图 4 高斯近似模型的交叉验证预测误差 Fig. 4 Cross-validation prediction error of Gaussian approximation model |
结果表明,各测试样本KT,10KQ,η0以及Pmin预测值与数值模拟结果的相对误差在3%以内,精度满足要求。训练后的多输出高斯代理模型可用于多目标优化中目标函数值的评估。
3 螺旋桨多目标优化 3.1 NSGA-Ⅱ多目标优化NSGA-Ⅱ是由Kalyanmoy等[23]针对NSGA提出的一种改进进化算法,该算法具有多目标适度简化、最优解拓展、求解精度高、高鲁棒性等优点,尤其针对少量样本、多变量多目标优化问题就较快的求解速度。
螺旋桨敞水性能多目标优化问题可以描述为:
{目标:{maxη0,minKQ,约束:{KT⩾KTinitial,KQ⩽KQinitial,Pmin⩾Pmin-initial,xmin⩽xa/b/c/d/e/f⩽xmax,0.4⩽Ae/Ad⩽1.4。 |
式中:叶宽
利用以100个样本建立的多输出高斯近似模型计算螺旋桨敞水性能,并以NSGA-Ⅱ算法进行螺旋桨的优化设计,将种群交叉遗传概率设置为0.92,变异概率设置为0.08。经过1000代的演化,可得到Pareto非支配最优解集。
初始种群和Pareto解集的扭矩系数和效率如图5所示。可以看出最优解是在原始设计空间之外的,其中Pareto前沿上A点为综合最优设计,表示Pareto前沿中的最优综合性能点,同时兼顾了扭矩系数和效率。
![]() |
图 5 初始种群和Pareto解集扭矩系数和效率 Fig. 5 Torque coefficient and efficiency of initial population and Pareto solution set |
选取优化结果中A点进行根据优化结果分析,原模型的效率提高了2.6%、扭矩降低了3.3%。此外,盘比降低了2.5%,但螺旋桨模型的空化特性并未恶化,因为推进器旋转区域的最小压力增加了7%。优化方案4桨叶r/R=0.7处剖面轮廓如图6所示。可以看出,经过优化,叶片的厚度和宽度发生了明显的变化,在相对半径0.7处变小,但倾角有所增加。
![]() |
图 6 优化方案4桨叶r/R=0.7处剖面轮廓 Fig. 6 Section profile at r/R=0.7 of the optimized scheme 4 blade |
桨叶宽度、剖面厚度、倾斜角沿径向的分布如图7所示。可以看出,经过优化,桨叶厚度、宽度、倾斜角在相对半径r/R=0.7处明显变小,而曲率在此处所有增加。
![]() |
图 7 KP505模型桨叶宽度、厚度和倾斜角沿半径分布 Fig. 7 KP505 model blade width, thickness and inclination angle distribution along the radius |
优化前后螺距比和曲率沿径向的分布如图8所示。优化前后螺旋桨的螺距比趋势相似,而且各半径处螺距比较对相近。此外,优化结果表明,螺旋桨螺距比在靠近桨叶根部较大。从推进设计的角度来看,螺旋桨螺距比在靠近桨叶根部的逐渐增大可以增加推力,从而提高效率。但是,通过增加螺距比,需要提高螺旋桨轴功率以及螺旋桨流动区域的最小压力。优化前后螺旋桨曲率沿半径分布特征相似,优化后的螺旋桨桨叶根部曲率变为负值,而且优化后的曲率变小,有利于降低螺旋桨噪声。
![]() |
图 8 KP505模型桨叶螺距比和曲率沿半径分布 Fig. 8 KP505 model blade pitch ratio and curvature distribution along the radius |
1)以有限元仿真结果为样本,以螺旋桨参数沿半径的分布为输入,基于多输出高斯近似模型建立螺旋桨水动力性能预测模型。该模型与数值仿真结果得出的螺旋桨推力系数和效率误差在3%以内,可以较为准确预测螺旋桨水动力性能。
2)应用本文提出的螺旋桨高效自动优化方法,开展KP505桨的优化设计,实现了对KP505桨降低扭矩系数、提高效率的多目标优化,形成了Pareto解集,获得了理论上的最佳设计方案。
3)根据对优化结果的分析,0.65R~0.95R处螺距比和0.4R~0.6R处倾斜角对KP505桨的效率有较大影响。
综上所述,本文方法可以有效提高螺旋桨优化的效率和自动化程度,并能得到最佳设计方案。在后续工作中,将进一步探讨应用该方法对考虑船后流场响应的螺旋桨优化设计的有效性。
[1] |
韩用波, 董郑庆, 吕江, 等. 优化理论在螺旋桨水动力设计中的应用[J]. 中国造船, 2019, 60(1): 52-59. HAN Yongbo, DONG Zhengqing, LV Jiang, et al. Application of optimization theory in propeller hydrodynamic design[J]. China Shipbuilding, 2019, 60(1): 52-59. |
[2] |
马丹萍, 李胜忠, 赵峰. 基于FFD几何变形方法的螺旋桨水动性能优化设计[C]//年中国造船工程学会优秀学术论文集, 2017, 39–47. MA Danping, LI Shengzhong, ZHAO Feng. Optimal design of propeller hydrodynamic performance based on FFD geometric deformation method [C]//Outstanding Academic Papers of Chinese Society of Shipbuilding Engineering, 2017, 39–47. |
[3] |
RAMIN T, KARI M. Hydrodynamic optimization of marine-propeller using gradient and non-gradient-based algorithms[J]. Acta Poly technical Hungarica, 2013(10): 221-237. |
[4] |
FLORIAN V, RICKARD J, RICKARD E B. On surrogate mehods in propeller optimization[J]. Ocean Engineering, 2014, 214-227. |
[5] |
叶礼裕, 王超, 孙文林, 等. 近似模型方法在螺旋桨优化设计中的应用[J]. 上海交通大学学报, 2016, 50(8): 1173-1179. DOI:10.16183/j.cnki.jsjtu.2016.08.006 |
[6] |
王超, 韩康, 孙聪, 等. 船用螺旋桨优化设计与参数分析[J]. 华中科技大学学报(自然科学版), 2020, 48(4): 97-102. DOI:10.13245/j.hust.200418 |
[7] |
SEDERBERG T W , PARRY S R. Free-form deformation of solid geometric models[C]// Proceedings of the 13 thannual conference on Computer Graphics and Interactive Techniques, 1986, 1151–1160.
|
[8] |
LAMOUSIN H J, WAGGENSPACK N N. NURBS-based free-form deformations[J]. IEEE Computer Graphics and Applications, 1994(6): 59-65. |
[9] |
NOBLE R A, CLAPWORTHY G J. Direct manipulation of surfaces using NURBS-based free-form deformations[C]// IEEE International Conference on Information Visualization, 1999, 1230–1239.
|
[10] |
张人会, 张书玮, 杨军虎, 等. 基于自由曲面变形方法的离心泵三维[J]. 工程热物理学报, 2018, 39(1): 98-103. |
[11] |
李冬琴, 李国焕, 戴晶晶, 等. 基于NFFD算法的船体几何变形技术[J]. 船舶工程, 2018, 40(6): 17-22. DOI:10.13788/j.cnki.cbgc.2018.06.017 |
[12] |
WU J W, LIU X Y, WAN D C. Multi-objective Hydrodynamic Optimization of Ship Hull Based on Approximation Model [C]//26th International Ocean and Polar Engineering Conference. 2016.
|
[13] |
马明生, 唐静, 李彬, 等. NFFD控制点分布对气动外形优化的影响[J]. 北京航空航天大学学报, 2016, 42(3): 1-9. MA Mingsheng, TANG Jing, LI Bin, et al. Influence of NFFD control point distribution on aerodynamic shape optimization[J]. Journal of Beijing University of Aeronautics and Astronautics, 2016, 42(3): 1-9. DOI:10.13700/j.bh.1001-5965.2016.0277 |
[14] |
唐静, 邓有奇. 飞翼气动优化中参数化和网格变形技术[J]. 航空学报, 2015, 36(5): 1480-1490. TANG Jing, DENG Youqi. Parameterization and mesh deformation technology in aerodynamic optimization of flying wings[J]. Chinese Journal of Aeronautics and Astronautics, 2015, 36(5): 1480-1490. |
[15] |
施法中. 计算机辅助几何设计与非均匀有理B样条[M]. 北京: 高等教育出版社, 2013.
|
[16] |
CHOI B K, YOO W S, LEE C S. Matrix representation for NURBS curves and B-spline curves[J]. Computer aided Design, 1990, 22(4): 235-240. DOI:10.1016/0010-4485(90)90052-E |
[17] |
REZA T M, MOZHGAN A, AZAR S, B. Solvinga multi-objective job shop scheduling problem with sequence-dependent setup times by a Paretoarchive PSO combined with genetic operators and VNS[J]. The International Journal of Advanced Manufacturing Technology volume, 2011, 53: 733-750. DOI:10.1007/s00170-010-2847-4 |
[18] |
高亮德. 基于多输出高斯过程回归的离心泵叶轮优化设计研究[D]. 兰州: 兰州理工大学, 2021.
|
[19] |
БурнаевЕ. В, ЗайцевА. А. Суррогатноемоделированиеразноточныхданныхвслучаевыборокбольшогоразмера[J]. Информационныепроцессы. 2015, 15: 97–109.
|
[20] |
LIU X, ZHU Q, LU H. Modeling Multi response surfaces for airfoil design with multiple-output-Gaussian-process regression[J]. Journal of Aircraft, 2014, 51(3): 740-747. DOI:10.2514/1.C032465 |
[21] |
ZHANG R, ZHAO X. Inverse method of centrifugal pump blade based on Gaussian process regression[J]. Mathematical Problems in Engineering, 2020, 4605625. |
[22] |
吴宽展, 刘学军, 吕宏强. 基于多输出高斯过程的超临界翼型优化[J]. 空气动力学学报, 2015, 33(6): 728-732. WU Kuanzhan, LIU Xuejun, LV Hongqiang. Optimization of supercritical airfoil based on multi-output Gaussian process[J]. Journal of Aerodynamics, 2015, 33(6): 728-732. DOI:10.7638/kqdlxxb-2014.0059 |
[23] |
DEB K, PRATAP A, AGARWAL S, et al. A fast and elitist multiobjective genetic algorithm: NSGA-II[J]. IEEE Transactions on Evolutionary Computation, 2002, 6(2): 182-197. DOI:10.1109/4235.996017 |