新型仿生螺旋桨优化设计

doi:10.16183/j.cnki.jsjtu.2022.174

新型仿生螺旋桨优化设计

吴春晓, 卢雨^,, 刘社文, 顾朱浩, 邵思雨, 邵武, 李闯

大连海事大学船舶与海洋工程学院,辽宁大连 116026

Optimization Design of New Bionic Propeller

WU Chunxiao, LU Yu^,, LIU Shewen, GU Zhuhao, SHAO Siyu, SHAO Wu, LI Chuang

College of Naval Architecture and Ocean Engineering, Dalian Maritime University, Dalian 116026, Liaoning, China

通讯作者: 卢雨,副教授;E-mail:luyu@dlmu.edu.cn.

责任编辑: 王一凡

收稿日期: 2022-05-20 修回日期: 2022-07-27 接受日期: 2022-08-18

基金资助:

国家自然科学基金面上项目(52171293)
国家自然科学基金青年项目(51809029)
大连市支持高层次人才创新创业项目(2020RQ009)
大连市科技创新基金重点学科重大项目(2020JJ25CY016)
辽宁省自然科学基金博士科研启动计划项目(2019-BS-025)
中央高校基本科研业务费专项资金资助(3132022112)

Received: 2022-05-20 Revised: 2022-07-27 Accepted: 2022-08-18

作者简介 About authors

吴春晓(1998-),硕士生.主要从事船舶推进水动力学研究.

摘要

对导边变形的仿生螺旋桨的水动力性能优化设计提出一种新颖的方法.基于仿生学原理与参数化建模的方法,将座头鲸前鳍的凹凸结构应用于螺旋桨导边,即在螺旋桨导边迎流区依据指数衰减曲线和正弦函数曲线将标准光滑导边进行类似座头鲸鳍突起结构状的凹凸变形,获得的导边凹凸的仿生螺旋桨.分别对指数衰减型仿生桨与正弦函数型仿生桨进行水动力性能、空泡性能以及噪声性能数值模拟.选出其中性能较优的螺旋桨,然后将基于仿真设计(SBD)技术引入新型仿生螺旋桨优化设计中,以控制导边变形的指数衰减曲线形状的参数为优化设计变量,以母型桨的转矩作为约束条件,选取敞水效率为目标函数,采用Sobol与T-Search优化算法,构建基于指数衰减曲线的仿生螺旋桨优化研究系统.研究结果表明:将座头鲸前鳍的凹凸结构应用于螺旋桨导边对螺旋桨的空泡性能与噪声性能有所提高,但对于螺旋桨敞水性能的提高不显著,验证了本研究所建立的仿生桨水动力性能优化设计方法是有效可靠的,为仿生螺旋桨的性能数值计算及构型优化设计提供了一定的理论依据和技术指导.

关键词： 仿生螺旋桨; 指数衰减曲线; 正弦函数曲线; 基于仿真设计(SBD)技术; 优化算法

Abstract

A novel method for optimal design of hydrodynamic performance of bionic propeller with a deformable leading edge is proposed. Based on the bionics principle and method of parameterized modeling, the fore-fin concave-convex structure of humpback whales is applied to the propeller leading edge, the leading edge in the propeller to meet flow region according to the exponential decay curve and the standard sine curve smooth leading edge for similar humpback fins protuberant structure of concave and convex deformation, and the leading edge of concave and convex bionic propeller. The hydrodynamic performance, the cavitation performance, and the noise performance of the exponential decay bionic propeller and the sinusoidal function bionic propeller were simulated respectively. The propeller with a better performance is selected, and the simulation based design (SBD) technology is introduced into the optimization design of the new bionic propeller. The parameters controlling the shape of the exponential attenuation curve of the guide edge deformation are taken as optimization design variables, the torque of the parent propeller is taken as the constraint condition, the open water efficiency is selected as the objective function, and the optimization algorithm of Sobol and T-Search is adopted. A bionic propeller optimization system based on the exponential decay curve is constructed. The results show that the application of the concave and convex structure of the humpback whale fore-fin to the guide edge of the propeller improves the cavitation performance and noise performance of the propeller, but the improvement of the open water performance of the propeller is not particularly significant. It is verified that the hydrodynamic performance optimization design method of the bionic propeller established in this paper is effective and reliable, which provides a certain theoretical basis and technical guidance for the performance numerical calculation and configuration optimization design of the bionic propeller.

Keywords： bionic propeller; exponential decay curve; sine function curve; simulation based design (SBD) technology; optimization

PDF (10709KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

吴春晓, 卢雨, 刘社文, 顾朱浩, 邵思雨, 邵武, 李闯. 新型仿生螺旋桨优化设计[J]. 上海交通大学学报, 2023, 57(11): 1421-1431 doi:10.16183/j.cnki.jsjtu.2022.174

WU Chunxiao, LU Yu, LIU Shewen, GU Zhuhao, SHAO Siyu, SHAO Wu, LI Chuang. Optimization Design of New Bionic Propeller[J]. Journal of Shanghai Jiaotong University, 2023, 57(11): 1421-1431 doi:10.16183/j.cnki.jsjtu.2022.174

仿生学是将自然界中动物的某些功能性特征应用到具体的工程研究中, 众多学者在螺旋桨的外形设计上根据仿生学的原理对模拟生物某些外形的桨叶形状进行优化.在研究座头鲸前鳍凹凸结构的基础上,国内外学者对其进行了研究.Fish等^[1]在20世纪末提出导缘凹凸这一概念,并且在风洞中完成了仿生鳍的试验,将仿生机翼和无仿生结构标准翼的空气力性能进行比较,证明凹凸导缘可以改善生物机翼的空气力性能.凹凸前缘结构可以得到提升升力系数、降低阻力系数的效果,同时有助于改善气动力的失速特性.Pedro等^[2]对仿生凹凸结构进行计算流体动力学(CFD)数值计算,探讨凹凸结构在低雷诺数下对机翼性能的影响.Nierop等^[3]采用CFD数值计算得到了带有凹凸结构翼型的表面压力.Sandberg等^[4]在NACA0012机翼的噪声性能预测中采用了直接数字模拟(DNS),根据流场的变化,发现流场的跨度关联性被仿生锯齿削弱了,涡流被分解成若干小块,加速了涡流的扩散,导致了噪声的降低.Chong等^[5]把粒子图像测速仪(PIV)技术引入到流场观测中,结果表明,随着后缘锯齿的安装,机翼的声压级有所下降,而且增加锯齿的宽度和高度会使降噪效果更好.Leon等^[6]研究了带有后缘锯齿的NACA0018机翼产生的宽带噪声.

基于仿真设计(Simulation Based Design,SBD)的概念是由意大利INSEAN水池的Campana教授提出的,旨在将数值模拟技术、全局优化算法及几何重构技术集成起来,从而形成一种由优化目标驱动形状改变的新模式^[7-8].SBD理论通过利用CFD技术对设定的优化目标(仿生螺旋桨的一项或多项性能)进行数值模拟,然后利用全局优化算法和几何重构技术对构型(仿生螺旋桨)设计空间进行探索寻优,从而得到满足给定约束条件的性能最优的构型.这种理论将大大促进仿生螺旋桨设计新方法的研究.

提出了基于指数衰减曲线及正弦函数曲线的导边变形的两种新型仿生螺旋桨.基于仿生学原理与曲线参数化建模,将凹凸结构应用于船舶螺旋桨导边.即在螺旋桨导边迎流区分别依据指数衰减曲线和正弦函数曲线将标准光滑导边进行类似座头鲸鳍突起结构状的凹凸变形,获得导边凹凸的仿生螺旋桨.分别对指数衰减型仿生桨与正弦函数型仿生桨进行水动力性能、空泡性能以及噪声性能数值模拟.选出其中性能较优的螺旋桨,然后将SBD技术引入新型仿生螺旋桨优化设计中.以控制导边变形的指数衰减曲线形状的参数为优化设计变量,以母型桨的转矩作为约束条件,以敞水效率为目标函数,采用Sobol与T-Search优化算法,构建基于指数衰减曲线的仿生螺旋桨优化研究系统.

1 仿生螺旋桨半参数化模型建立

1.1 母型桨模型建立

采用的母型桨是AU-5-50,螺旋桨主要参数^[9]:直径D=0.25 m;桨叶数L=5;盘面比A_E/A₀=0.5;螺距比P_mean/D=1;纵倾角α=10°.根据AU-5-50的数据进行建模,模型如图1所示.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 母型螺旋桨AU-5-50

Fig.1 Parent propeller AU-5-50

1.2 指数衰减曲线型与正弦函数型仿生桨模型

建立通过研究座头鲸鳍外形图像发现突起结构并非以往研究中采用的简单正弦函数曲线,其突起结构部分变化与指数衰减曲线形式存在一定相似性.座头鲸鳍自根部起第一个结节凸起幅度较高,自中部起产生又一个比较明显的凸起之后跟随几个幅度相差较小的凸起直至鳍梢处结束.测试发现指数衰减曲线以正弦函数或余弦函数等规律前进能够很好地拟合座头鲸鳍的前缘突起变化结构.指数衰减曲线方程如下:

(1)N(t)=N₀e^-^λt

式中:N(t)为数量;N₀为原始数量;λ为正数为衰减量的衰减常数;t为时间.

由于座头鲸鳍存在两个较为明显的突起结构,故变形曲线选择主要由两条指数衰减曲线组合,定义变形曲线如下:

(2)

\begin{array}{l} X (t) = - A_{n} e^{- E_{n} t} s i n (T_{n} t) \\ Z (t) = K t + Z_{n} \end{array}\}

式中:A_n、E_n、T_n分别代表振幅、衰减和周期;X(t)为不同时刻下曲线的值;Z(t)为方向控制曲线的高度,定义指数衰减曲线影响导边的程度;K为常数;Z_n代表起始高度.根据导边的范围以及座头鲸鳍前缘凹凸形式定义Z方向上设置两条指数衰减曲线;应用F-spline对两条指数衰减曲线进行连接,确保曲线光顺性,最终将F-spline与指数衰减曲线结合得到整体变形曲线如图2所示.

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 指数衰减周期拟合曲线示意图

Fig.2 Schematic diagram of exponential decay period fitting curve

将得到的变形曲线施加于三维螺旋桨AU-5-50导边,对其导缘进行参数化曲线控制变形,绘制出指数衰减曲线型仿生桨如图3所示.

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 具有指数衰减曲线的仿生螺旋桨

Fig.3 Bionic propeller with exponential decay curve

为进行与指数衰减曲线型桨进行对比,建立了导边变形曲线为正弦函数的仿生螺旋桨,定义变形曲线如下:

(3)

\begin{array}{l} X (t) = s i n (T_{n} t) \\ Z (t) = K t + Z_{n} \end{array}\}

正弦函数型仿生桨模型如图4所示.

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 正弦函数曲线仿生螺旋桨

Fig.4 Bionic propeller of sine function curve

2 仿生螺旋桨数值预报方法

2.1 数值模拟设置

进行螺旋桨真实的旋转环境的模拟,创建计算域:旋转域与静止域.旋转域直径为1.25D, 长度为1.1D,静止域直径为5D,静止域速度进口至螺旋桨距离为3D,静止域压力出口至螺旋桨距离为6D.静止域边界条件设定:静止域左边设为速度进口,静止域右边设为压力出口,轴向方向设为对称平面边界条件,旋转域与静止域创建交界面.螺旋桨计算域如图5所示.

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 螺旋桨计算域

Fig.5 Calculation domain of propeller

网格化分旋转域选择切割体网格与棱柱层网格,静止域选择切割体网格.螺旋桨桨盘面处网格与轴向切面网格如图6所示.

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 螺旋桨网格示意图

Fig.6 Diagram of propeller meshing

2.2 数值方法

在进行水动力性能计算时,物理模型选择K-Epsilon 模型;在进行空泡性能计算时,选择延迟分离涡模拟(DDES)模型和用于多相相互作用的Schnerr-Sauer模型^[10⇓-12].Schnerr-Sauer空化模型基于Rayleigh-Plesset (RP) 简化方程并忽略气泡增长加速的影响、黏性效应和表面张力效应,可用于对单组分材料和多组分材料的气泡增长率和塌陷率进行比例缩放.蒸汽相质量分数的输运方程如下:

(4)

\frac{\partial (ρ w_{v})}{\partial t}

\frac{\partial (ρ u_{j} w_{v})}{\partial x_{j}}

=m_e-m_c

式中:ρ为混合介质密度;m_e为气泡增加质量;m_c为气泡减少质量;w_v为蒸汽的质量分数;u_j为速度分量;x_j为空间分量.空泡生成速度如下:

(5)

v_{r}^{2}

\frac{2}{3} \frac{p_{s a t} - p}{ρ_{1}}

式中:p_sat为饱和压力,对应于气泡表面处的温度;p为周围液体的压力;ρ₁为液体密度^[13].

FW-H积分方程^[14]在螺旋桨的噪声预测等领域得到了广泛的应用,FW-H方程是将连续性方程和动量方程精确地重新整理为不均匀波方程形式的结果.即使在积分表面位于非线性流体区域中的情况下,FW-H方程的结果同样精确.它根据自由空间格林函数来计算观察器位置x处的声压.对于区域或一组表面中的流体辐射到静止介质中的压力,FW-H方程^[15]的实质是从非线性流场中提取出等价声源, 将流场计算域声场计算解耦分开计算, 使得问题得到简化.FW-H方程如下:

(6)

\begin{array}{l} P' (x, t) = P'_{T} (x, t) + P'_{L} (x, t) + P'_{Q} (x, t) \\ P'_{T} (x, t) = \frac{1}{4 π} \frac{\partial}{\partial t} \int_{S} [\frac{Q}{r (1 - M r)}] d S \\ P'_{L} (x, t) = \frac{1}{4 π} (- \frac{\partial}{\partial x_{i}} \int_{S} [\frac{L_{i}}{r (1 - M r)}] d S) \\ P'_{Q} (x, t) = \\ \frac{1}{4 π} \frac{\partial^{2}}{\partial x_{i} \partial x_{j}} \int_{V} [\frac{T_{i j}}{r (1 - M r)}] d V \end{array}\}

(7)

\begin{array}{l} Q = ρ_{0} U_{i} n_{i} \\ U_{i} = (1 - \frac{ρ}{ρ_{0}}) v_{i} + \frac{ρ u_{i}}{ρ_{0}} \\ L_{i} = P_{i j} n_{i} + ρ u_{i} (u_{n} - v_{n}) \end{array}\}

式中:P'_T(x, t)、P'_L(x, t)、P'_Q(x, t)分别代表单极子、偶极子和四极子源项.单极子是由固体表面有规律排开流体引起的;偶极子是由流体施加在固体表面所引起的非定常负荷作用;四极子是由于湍流应力张量所引起.r为源点到观察器的距离;Mr为源点朝向观察器的马赫数;u_i表示i方向上的流体速度分量;u_n为与表面垂直的流体速度分量;v_i表示i方向上的表面速度分量;v_n为与表面垂直的表面速度分量;n_i为表面法向矢量;ρ₀为远场密度;P_ij为压缩应力张量;T_ij为Lighthill应力张量.

2.3 数值验证

由于网格的数量没有固定的标准,且网格是在数值模拟计算中最容易受到人为因素干扰的一个关键环节,因而网格无关性验证就显得尤为重要.进行网格无关性验证时使用网格数量分别为146万、202万、277万、392万的网格进行验证.不同网格数量的网格计算的螺旋桨在进速系数J=0.4下的转矩系数10K_Q、敞水效率η₀、推进系数K_T的CFD计算结果与试验结果(exp)如图7所示.

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 计算结果与试验结果的比较(J=0.4)

Fig.7 Comparison of calculation results and experiment results (J=0.4)

通过图7的计算结果,得出在4种不同网格数量下的10K_Q、η₀、K_T与母型桨的误差对比.根据比较结果可知,网格数量为202万时已经达到了网格独立性的要求,另外网格数量为277万和392万的10K_Q、η₀、K_T的相对差值分别为0.095%、0.24%、0.29%.采用网格数量为277万的网格能够得到可靠的螺旋桨水动力性能结果,同时也满足了尽可能减少所需要消耗机时的要求.因此本次模拟中使用网格数量为277万的网格进行计算.

使用该网格验证AU-5-50型螺旋桨,如果该型螺旋桨用此套网格数值模拟的数据与其试验值相差较小,说明用于该桨数值模拟的网格是合理的且具有可靠性的.图8所示为CFD数值计算AU-5-50型螺旋桨水动力结果误差对比,由图8可见,10K_Q、η₀、K_T与试验值误差均在5%以内,说明用此套网格计算的结果是可信的.

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 计算结果与试验结果的比较

Fig.8 Comparison of calculation results and experiment results

3 仿生螺旋桨性能对比分析

3.1 水动力性能对比

建立了母型桨(parent)、指数衰减曲线型仿生螺旋桨(decay)、正弦函数曲线型仿生螺旋桨(sin).并对其进行了CFD数值模拟^[16-17],得到了全进速下仿生螺旋桨与母型桨水动力性能对比如图9所示.指数衰减曲线型仿生桨在J=0.4~0.7时,η₀有明显增大,但10K_Q、K_T有所降低.

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 指数衰减曲线型仿生螺旋桨与母型螺旋桨水动力结果的比较

Fig.9 Comparison of hydrodynamic results of exponential decay curve type bionic propeller and parent propeller

图10为正弦函数曲线型螺旋桨与母型桨水动力性能对比,正弦函数型螺旋桨在J=0.5~0.7时,η₀有所提升,在其他进速下的水动力性能效果不是很明显.由此可以看出,以AU-5-50为母型的仿生螺旋桨对水动力性能的改善不是特别明显.

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10 正弦函数曲线型仿生螺旋桨与母型螺旋桨水动力性能比较

Fig.10 Comparison of hydrodynamic results of sine function curve type bionic propeller and parent propeller

图11为指数衰减曲线型、正弦函数曲线型仿生螺旋桨与母型桨在J=0.5时的叶面表面压力p对比图.在凹凸结构的影响下,仿生螺旋桨叶背压力分布规律发生了变化.母型桨的负压最小值为-30 078 Pa,正压最大值为14 966 Pa;指数衰减曲线型的负压最小值为-24 762 Pa,正压最大值为16 210 Pa;正弦函数曲线型负压的最小值为-28 710 Pa,正压最大值为15 784 Pa.对比发现仿生螺旋桨的负压值与正压值均有所增大,表明在凹凸结构的影响下导边压力过渡的较为平缓,叶背表面压力分布更加均匀,这也使得仿生螺旋桨敞水效率在一定程度上得到了提升并有可能改善螺旋桨的空泡性能.

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11 母型螺旋桨、指数衰减曲线型和正弦函数曲线型的叶背表面压力对比图

Fig.11 Comparison of surface pressure of parent propeller blade, exponential decay curve type, and sine function curve type

螺旋桨桨后涡流场能够体现螺旋桨在流场内的涡流状态及能量分布,有助于对螺旋桨敞水性能进行分析.通过对比分析各个螺旋桨桨后涡流场,能够在一定程度上分析凹凸结节结构对螺旋桨性能的影响.在J=0.5时,对比分析各个仿生桨与母型桨轴向涡量(Ω)图.母型桨的最大涡量值为45.9 s^-1,最小涡量为-85.2 s^-1;指数衰减曲线型的最大涡量值为45.4 s^-1,最小涡量为-84.4 s^-1;正弦函数曲线型的最大涡量值为45.2 s^-1,最小涡量为-84.2 s^-1.图12所示为母型螺旋桨、指数衰减曲线型和正弦函数曲线型仿生桨的涡量对比.由图可见,指数衰减曲线型、正弦函数曲线型涡量值的绝对值均低于母型桨.说明凹凸结构可以将涡分解,加速了涡的扩散,能够有效提升螺旋桨的水动力性能.

图12

新窗口打开| 下载原图ZIP| 生成PPT

图12 母型螺旋桨、指数衰减曲线型和正弦函数曲线型仿生桨的涡量对比图

Fig.12 Comparison of axial vorticity of parent propeller blade, exponential decay curve type, and sine function curve type

3.2 空泡性能对比

在进行各个仿生螺旋桨水动力性能数值模拟之后,对各个仿生螺旋桨以及母型桨进行了空泡性能的数值模拟,对空泡性能的计算工况取空泡数为2.024,空化压力在温度为23.2 ℃时为2 818 Pa.图13所示为母型桨、指数衰减曲线型、正弦函数曲线型在J=1时的桨叶空泡体积分数φ的对比.母型桨φ的最大值为0.962,指数衰减曲线型φ的最大值为0.928,正弦函数曲线型φ的最大值为0.842.由图可见,仿生桨的空化面积与φ均略小于母型桨,说明仿生螺旋桨对母型桨的空泡性能有所提升.而优化指数衰减曲线型螺旋桨在凹凸结节处可能会出现绝对压力低于水的饱和蒸气压的情况,并在此处发生空化现象.整体来说,仿生螺旋桨对螺旋桨的空化性能有所提升.

图13

新窗口打开| 下载原图ZIP| 生成PPT

图13 母型螺旋桨、指数衰减曲线型和正弦函数曲线型仿生桨的空泡体积分数对比图

Fig.13 Comparison of volume fraction of vapor of parent propeller blade, exponential decay curve type, and sine function curve type

3.3 噪声性能对比

在进行各个仿生螺旋桨空泡性能数值模拟之后,对各个仿生螺旋桨以及母型桨进行了噪声性能的数值模拟,仿生桨与母型桨模型特征点位置噪声总声级如表1所示,母型螺旋桨、指数衰减曲线型和正弦函数曲线型仿生桨的声压级对比如图14所示,纵轴SPL表示特征点位置处不同频率f下的噪声声压级.通过表1与图14,综合分析指数衰减函数型、正弦函数曲线型仿生螺旋桨与母型桨在特征点(桨上方一倍直径处)处的声压频谱曲线,可以看出:相较于母型桨,各个仿生螺旋桨在特征点位置的噪声总声级都有所下降,其中指数衰减函数曲线型仿生桨减小幅度较大, 相比母型桨降低了1.43%,正弦函数曲线型仿生桨相比母型桨降低了1.35%.表明了仿生结构可以将涡分解,加速了涡的扩散,能够有效降低螺旋桨的噪声,提高水下航行器的隐身性能.将仿生结构在螺旋桨上进行合理的布置与运用,能有效改善其噪声性能.

表1 仿生桨与母型桨模型特征点位置噪声总声级

Tab.1 Total noise level of model feature points of bionic propeller and prototype propeller

模型	特征点位置噪声总声级/dB	总声级降低比例/%
母型桨	147.32	—
指数衰减函数曲线型	145.21	1.43
正弦函数曲线型	145.32	1.35

新窗口打开| 下载CSV

图14

新窗口打开| 下载原图ZIP| 生成PPT

图14 母型螺旋桨、指数衰减曲线型和正弦函数曲线型仿生桨的声压级对比图

Fig.14 Comparison of sound pressure spectrum curve of parent propeller blade, exponential decay curve type, and sine function curve type

4 仿生螺旋桨水动力性能优化

4.1 构建优化设计系统

通过对比上文建立的指数衰减函数型仿生螺旋桨与正弦函数型仿生螺旋桨,发现指数衰减函数型仿生螺旋桨整体性能较优.参照文献[18],本文以指数衰减型仿生桨进行水动力性能优化设计系统的建立.图15为仿生螺旋桨优化设计系统的整个流程.

图15

新窗口打开| 下载原图ZIP| 生成PPT

图15 优化设计系统

Fig.15 Optimization design system

根据指数衰减曲线对螺旋桨进行几何半参数化建模.因此将曲线的振幅、衰减和周期作为设计优化变量,通过指数衰减曲线的变化达到对螺旋桨几何模型的改变.选取敞水效率为目标函数,采用Sobol与T-Search优化算法^[19],从而形成一种由目标驱动设计的形状构型设计新模式.

4.2 灵敏度分析与最优解的选取

将曲线振幅A₀₁/A₀₃、衰减E₀₁/E₀₃和周期T₀₁/T₀₃作为设计优化变量,优化目标为仿生桨的η₀,对其进行全局目标优化.优化过程为先使用Sobol算法对其进行全局均匀搜索,找到较优的个体,然后以此为基础用T-Search算法对其进行目标优化.以此进行迭代,从而得到最优桨型方案.Sobol搜寻了100个个体,A₀₁/A₀₃、E₀₁/E₀₃和T₀₁/T₀₃与η₀、K_T和10K_Q的相关性如图16所示,对其进行了线性拟合发现,A₀₁的变化对η₀、K_T和10K_Q影响较小,A₀₃的变化对η₀、K_T和10K_Q影响较大,E₀₁/E₀₃的变化对η₀、K_T和10K_Q影响不是特别显著,周期T₀₁的变化对η₀、K_T和10K_Q影响不是特别显著,而周期T₀₃的变化对η₀、K_T和10K_Q影响较为显著.由Sobol优化可知,A₀₃和T₀₃对η₀、K_T和10K_Q影响较为显著,在接下来使用T-Search优化算法进行优化中可以着重考虑.

图16

新窗口打开| 下载原图ZIP| 生成PPT

图16 仿生螺旋桨衰减曲线振幅A₀₁/A₀₃、衰减E₀₁/E₀₃和周期T₀₁/T₀₃的相关性

Fig.16 Correlation of attenuation curve amplitude A₀₁/A₀₃, attenuation E₀₁/E₀₃ and period T₀₁/T₀₃ of bionic propeller

在T-Search优化中选取对仿生桨影响较大的参数进行优化,图17为敞水效率的收敛曲线以及最优解(optimized)的选取.

图17

新窗口打开| 下载原图ZIP| 生成PPT

图17 敞水效率收敛曲线及最优解的选取

Fig.17 Convergence curve of open water efficiency and selection of optimal solution

对最优解进行全进速下的敞水性能进行计算,如图18所示.从对比图中可以看出:仿生桨的敞水效率在J=0.1~0.8下高于母型桨,但是可能由于优化的个体较少,以及计算资源的限制,优化的效果不是很显著.

图18

新窗口打开| 下载原图ZIP| 生成PPT

图18 优化仿生螺旋桨与母型螺旋桨水动力结果的比较

Fig.18 Comparison of hydrodynamic results of optimized bionic propeller and parent propeller

5 结论

利用半参数化建模的方式建立指数衰减曲线型仿生桨以及正弦曲线型仿生桨的模型,然后对其进行水动力性能、空泡性能以及噪声性能的数值模拟.对比分析发现仿生螺旋桨的空泡性能与噪声性能优于母型桨.选取指数衰减函数型仿生桨对其进行水动力性能优化设计研究.将螺旋桨的几何重构技术、水动性能评估技术、最优化技术结合起来,实现整个设计过程的自动化.首先采用半参数化建模方法完成螺旋桨的几何建模,基于数值模拟技术进行螺旋桨水动力性能.将曲线的振幅、衰减和周期作为设计优化变量,以母型桨的转矩作为约束条件,选取敞水效率为目标函数,以水动力性能优化目标函数,借助全局智能优化算法完成仿生螺旋桨优化设计系统的构建.得到以下结论:

(1) 将座头鲸前鳍的凹凸结构应用于螺旋桨导边对螺旋桨的空泡性能与噪声性能有所提高.

(2) 优化设计结果表明螺旋桨的敞水效率提高不是特别显著,可以尝试对仿生螺旋桨噪声性能的优化设计研究.

(3) 将SBD技术应用于仿生螺旋桨性能优化设计系统,证明了该技术在螺旋桨性能设计领域的可行性.

该研究为仿生螺旋桨的性能数值计算及构型优化设计提供了一定的参考和借鉴作用.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

FISH

F E

, BATTLE

J M

Hydrodynamic design of the humpback whale flipper

[J]. Journal of Morphology, 1995, 225(1): 51-60.

PMID:7650744 [本文引用: 1]

The humpback whale (Megaptera novaeangliae) is reported to use its elongate pectoral flippers during swimming maneuvers. The morphology of the flipper from a 9.02-m whale was evaluated with regard to this hydrodynamic function. The flipper had a wing-like, high aspect ratio planform. Rounded tubercles were regularly interspersed along the flipper's leading edge. The flipper was cut into 71 2.5-cm cross-sections and photographed. Except for sections near the distal tip, flipper sections were symmetrical with no camber. Flipper sections had a blunt, rounded leading edge and a highly tapered trailing edge. Placement of the maximum thickness placement for each cross-section varied from 49% of chord at the tip to 19% at mid-span. Section thickness ratio averaged 0.23 with a range of 0.20-0.28. The humpback whale flipper had a cross-sectional design typical of manufactured aerodynamic foils for lift generation. The morphology and placement of leading edge tubercles suggest that they function as enhanced lift devices to control flow over the flipper and maintain lift at high angles of attack. The morphology of the humpback whale flipper suggests that it is adapted for high maneuverability associated with the whale's unique feeding behavior.

[2]

PEDRO

, KOBAYASHI

Numerical study of stall delay on humpback whale flippers

[C]// 46th AIAA Aerospace Sciences Meeting and Exhibit. Reno, Nevada, USA: American Institute of Aeronautics and Astronautics, 2008: 1-8.