留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

2019年  第34卷  第4期

动力传输
刷式密封摩擦生热温度场数值计算及试验
吴施志, 江平, 力宁, 潘君, 王智加
2019, 34(4): 737-743. doi: 10.13224/j.cnki.jasp.2019.04.001
摘要:
在摩擦热热源条件下建立了CFD多孔介质数值计算模型,对某型航空发动机刷式密封摩擦生热温度场进行计算分析,开展了刷式密封装置全工况条件下摩擦生热试验研究,采用红外成像技术实现了刷式密封动态温度场实时监测。根据试验结果对摩擦生热温度场计算方法进行了修正,引入了刚度修正系数,并对刚度修正系数进行了确定及验证。总结出经过试验验证的航空发动机刷式密封摩擦生热温度场计算分析及试验方法,结果表明:与全工况试验结果相比,计算误差值从48.15%减少到10.67%。
基于ACMPE、ISSL-Isomap和GWO-SVM的行星齿轮箱故障诊断
戚晓利, 王振亚, 吴保林, 叶绪丹, 潘紫微
2019, 34(4): 744-755. doi: 10.13224/j.cnki.jasp.2019.04.002
摘要:
针对从行星齿轮箱非线性、非平稳振动信号特征提取困难的问题,提出了一种基于自适应复合多尺度排列熵(ACMPE)、改进监督型自组织增量学习神经网络界标点等度规映射(ISSL-Isomap)和灰狼群优化支持向量机(GWO-SVM)相结合的行星齿轮箱故障诊断方法。利用ACMPE从复杂域提取振动信号的故障特征,构建高维故障特征集;采用ISSL-Isomap方法对高维故障特征集进行维数约简,提取出低维、敏感故障特征;应用GWO -SVM分类器对低维故障特征进行模式识别,判断故障类型。行星齿轮箱故障诊断实验结果分析表明:与多尺度排列熵(MPE)、复合多尺度排列熵(CMPE)等特征提取方法相比,ACMPE方法在分类效果和识别精度上更具优势;与局部切空间排列(LTSA)、等度规映射(Isomap)、加权Isomap(W-Isomap)、监督Isomap(S-Isomap)和监督型自组织增量学习神经网络界标点Isomap(SSL-Isomap)等降维方法进行比较,ISSL-Isomap方法降维效果最佳;所提方法的故障识别率达到100%,具有一定优越性。
一种基于Deep-GBM的航空发动机中介轴承故障诊断方法
田晶, 李有儒, 艾延廷
2019, 34(4): 756-763. doi: 10.13224/j.cnki.jasp.2019.04.003
摘要:
针对航空发动机中介轴承故障信号难于识别的特点,提出了一种深度梯度提升模型(Deep-GBM)对振动信号特征进行逐层学习以提高分类模型的准确率。开展某型航空发动机中介轴承故障模拟实验,并采用经验模式分解(EMD)方法对采集的振动信号进行分解,提取内蕴模式函数(IMF)分量非线性动力学参数样本熵作为原始故障特征。采用Deep-GBM对中介轴承内环故障、内环和滚动体综合故障、正常、滚棒剥落、滚棒划伤五种不同状态进行识别。实验结果表明,所提出的Deep-GBM故障诊断准确率达到87%,相对于传统的机器学习模型准确率最高提升了28%,并具有良好的泛化能力。
基于MED和LMD的自动倾斜器轴承广义Shannon熵谱分析
张先辉, 李新民, 金小强
2019, 34(4): 764-771. doi: 10.13224/j.cnki.jasp.2019.04.004
摘要:
针对轴承信号微弱故障特征易被强背景噪声淹没的问题,提出采用最小熵反褶积,通过逆滤波器最优化设计,对目标信号进行降噪处理,其峭度值提高了约3.8倍,增强了信号的微弱故障特征;针对非平稳非线性信号频率成分复杂难以解调的问题,提出采用局部均值分解(LMD)和峭度-相关系数筛选准则,其可对非平稳非线性信号进行自适应分解和最优重构,提高了信号的信噪比;针对信号耦合调制及边频突出的问题,通过引入广义Shannon熵进行包络谱带内降噪处理,信号一阶故障特征调制频率与故障特征频率的幅度比降低了24%~43%。通过实验室信号及某型直升机自动倾斜器轴承故障诊断地面试验的分析结果验证了该方法的合理性和可行性。
基于强抗噪威格纳威利分析的滚动轴承故障诊断
王宏超, 杜文辽
2019, 34(4): 772-777. doi: 10.13224/j.cnki.jasp.2019.04.005
摘要:
为了解决威格纳威利谱(Wigner-Vile spectrum, WVS)时频分析方法对受强背景噪声影响下的滚动轴承冲击性故障信号特征提取难的问题,根据滚动轴承发生故障时呈现出的循环平稳特征,将基于2阶循环统计量的循环谱密度算法(cyclic spectral density, CSD)与WVS相结合,提出基于CSD的WVS分析方法,即循环谱密度威格纳威利谱(CSDWVS)时频分析方法。经仿真及实验验证,相对传统WVS分析方法,该方法能有效提取出强背景噪声影响下的滚动轴承内圈故障特征频率为51.9 Hz及外圈故障特征频率为32.1 Hz。
球轴承内圈剥落缺陷双冲击特征动力学建模
罗茂林, 郭瑜, 伍星
2019, 34(4): 778-786. doi: 10.13224/j.cnki.jasp.2019.04.006
摘要:
基于Hertz接触理论,以内圈滚道表面存在局部剥落的球轴承为研究对象,对球轴承局部剥落故障所激起的非线性激励机理进行研究,分析球轴承双冲击现象激励机理。将滚动体与球轴承内外圈非线性接触特性纳入考虑,提出时变位移激励与时变接触力激励相耦合的球轴承局部剥落双冲击现象动力学模型。研究中采用Runge-Kutta数值积分法求解该二阶非线性动力学微分方程,并对故障球轴承仿真信号进行双冲击时间间隔特征验证。在不同转速和内圈剥落尺寸条件下开展仿真与实验研究,通过仿真和实测信号的双冲击时间间隔与理论双冲击时间间隔对比分析,仿真信号的双冲击时间间隔与理论值对比的误差结果皆小于2%,与实测信号的双冲击时间间隔值对比的误差结果皆小于12%;同时,仿真信号与实测信号具有很好的相似性,从而验证了该动力学模型的有效性。
一种关于静压气体轴承节流孔系数的计算方法
张建波, 焦春晓, 邹冬林, 塔娜, 饶柱石
2019, 34(4): 787-795. doi: 10.13224/j.cnki.jasp.2019.04.007
摘要:
基于层流边界层方程的分离变量算法和雷诺方程的解析算法,提出了一种关于单节流孔静压气体止推轴承的节流孔系数的计算方法。该方法通过比较层流边界层方程计算获得的气体轴承的质量流量和雷诺方程计算获得的质量流量计算获得了节流孔系数。将计算获得的节流孔系数和节流孔系数为常数0.8代入单节流孔气体止推轴承的雷诺方程中,计算获得的承载力与分离变量算法求解层流边界层方程获得的承载力进行对比,可以发现,相对于采用节流孔系数为0.8来说, 采用该计算的节流孔系数求解雷诺方程的承载力与分离变量算法求解获得的承载力结果精度最大提高了8%。从而验证了该计算节流孔系数方法的正确性。
弧线齿面齿轮齿面几何设计及径向齿宽特性
盛伟, 王利霞, 冯占荣, 王文轩, 陈勇, 赵玉龙
2019, 34(4): 796-803. doi: 10.13224/j.cnki.jasp.2019.04.008
摘要:
为适应航空领域对面齿轮传动的设计要求,综合弧线齿圆柱齿轮、面齿轮传动的各项优点,提出了弧线齿面齿轮齿面设计并研究了其径向齿宽特性。以端面渐开线齿廓弧线齿圆柱齿轮作为产形齿轮通过包络成形建立面齿轮齿面数学模型;结合齿面方程在MATLAB中进行齿面可视化分析,并通过MATLAB软件的数值计算和CATIA软件的复杂曲面造型实现面齿轮的参数化建模;确立内端根切、外端顶尖对径向有效齿宽的限制条件,计算面齿轮的径向有效齿宽,分析刀具齿轮的齿线圆弧半径和坐标系位置参数对径向有效齿宽的影响。研究表明:当增大齿线的圆弧半径,面齿轮的内径减小,外径、径向有效齿宽均增大;当增大坐标系位置参数,面齿轮的内径、外径均增大,径向有效齿宽减小。
气动热力学与总体设计
真实进气条件下发动机气动稳定性计算方法
谢业平, 刘永泉, 潘宝军
2019, 34(4): 804-812. doi: 10.13224/j.cnki.jasp.2019.04.009
摘要:
基于平行压气机理论,采用非定常、二维、无黏的积分型欧拉方程,以蒸汽吸入试验中产生的温度和压力组合畸变为边界条件,开展真实进气条件下温度和压力组合畸变对发动机影响的计算方法研究。在保证畸变强度、高温区范围、低压区范围等特征不变的前提下,对畸变图谱进行规则化等效转换,同时协调温度和压力畸变图谱,使两者可以进行相互叠加计算。结果表明:当周向子发动机数量大于2时,子发动机数量对计算结果影响不大;组合畸变对风扇、压气机压比影响规律与单独的压力畸变相同;与试验的推力损失对比显示,对畸变图谱进行等效规则化转换,可以用于真实组合畸变的计算。
某型航空发动机起动故障试验
姜新瑞, 王赫, 唐震
2019, 34(4): 813-820. doi: 10.13224/j.cnki.jasp.2019.04.010
摘要:
为解决某型航空发动机冷热态起动不兼容问题,基于发动机室内试车台,对起动过程进行研究,分析了起动故障原因。研究了起动机功率、起动机脱开转速及起动供油逻辑对起动性能的影响,通过多种调整措施,对起动过程进行了优化,并试验验证了措施的有效性。结果表明:起动机功率提高10%,起动时间缩短5%;起动机脱开转速提高2%,起动时间缩短18%~25%;调整优化起动供油逻辑,增加起动前中期供油量,减少后期供油量,可有效改善发动机的起动性能。
吸气式飞行器进气道唇口三维激波/激波干扰
姜宝森, 张亮, 李俊红, 施家桐
2019, 34(4): 821-828. doi: 10.13224/j.cnki.jasp.2019.04.011
摘要:
对吸气式飞行器进气道唇口处三维曲面激波/弓形激波干扰流场进行数值模拟,利用典型三维气动干扰试验对采用的数值计算方法进行验证。利用拼接网格技术及逆距离加权插值方法获得入口处流场的非守恒变量,作为激波干扰研究的入口边界条件。数值模拟表明,唇口处激波干扰流动的三维效应十分显著,曲面激波与弓形激波产生斜交,尽管唇口前缘半径很小,但Edney提出的6类激波干扰类型可能沿唇口展向方向同时存在;第Ⅲ和Ⅳ类激波/激波干扰的诱导使得唇口热流分布异常严酷;激波相交处形成斜向“伤疤”状局部高热流条带,峰值热流达到参考热流的4~6倍,可能引起唇口结构的局部烧蚀或破坏。
对转螺旋桨流场气动干扰数值模拟
史文博, 李杰
2019, 34(4): 829-837. doi: 10.13224/j.cnki.jasp.2019.04.012
摘要:
基于结构网格动态面搭接技术,通过求解三维非定常雷诺平均纳维-斯托克斯(RANS)方程,对单独螺旋桨(SRP)滑流流场进行了数值模拟,结果与风洞试验数据吻合良好,在此基础上深入研究了对转螺旋桨(CRP)的非定常滑流效应及前后桨之间的气动干扰现象。计算结果表明:所采用的非定常RANS方法能够良好刻画对转螺旋桨滑流流场的发展变化特征,适用于分析前后桨气动干扰特性。前后桨桨尖涡的复杂相互作用是导致对转螺旋桨气动干扰的直接原因,后桨桨尖涡对于前桨桨尖涡具有一定耗散作用。在气动干扰影响下,对转螺旋桨气动力会产生周期性波动,波动周期与螺旋桨桨叶数相关, 6×6对转螺旋桨在1个旋转周期内产生了12次波动。气流穿过对转螺旋桨会发生两次加速,因此在相同工况下,对转螺旋桨的拉力系数和功率系数相比于单独螺旋桨分别增大约1倍,效率提升约1.5%。
悬停状态下旋翼噪声声源反演的正则化方法
赵庆贺, 李晓东
2019, 34(4): 838-846. doi: 10.13224/j.cnki.jasp.2019.04.013
摘要:
根据旋翼辐射的气动噪声求解气动声学反问题,给定合理的旋翼表面载荷并结合声类比方法得到观测点的信号,以此进行表面载荷反演。对建立的大型不适定性方程组采用奇异值分解结合正则化方法进行求解,为克服反演结果存在的数值振荡,采用多项式函数来表示未知量,并引入库塔条件作为先验信息。结果表明:理想条件下,该方法得到的旋翼表面的压力分布更接近输入值,尤其在载荷分布较大的叶尖区域,考虑干扰噪声影响后,该方法得到的展向截面吸力峰位置的最大相对误差从15%减小到4%。
结构、强度、振动
基于最差件等寿命曲线的桨轴高低周复合疲劳寿命分析
章胜, 张勇, 李锦红
2019, 34(4): 847-851. doi: 10.13224/j.cnki.jasp.2019.04.014
摘要:
以Goodman曲线为基础,构造了考虑应力集中系数和疲劳极限强度分散系数影响的最差件等寿命曲线,然后基于最差件等寿命曲线建立了桨轴高低周复合疲劳寿命分析方法,用于预测桨轴的安全寿命,并给出了预测方法的适用范围为低周103~105次循环。采用此方法对两种载荷下的某型航空发动机桨轴进行了高低周复合疲劳寿命预测,并与试验结果进行了对比。结果表明:基于最差件等寿命曲线的桨轴高低周复合疲劳寿命分析方法成功预测一组安全寿命低于低周1 000次循环的载荷,且另一组载荷预测结果也与试验结果相吻合,预测方法是可行的。
三维四向C/C复合材料高温氧化环境强度预测
陈波, 温卫东, 崔海涛, 张宏建, 徐颖
2019, 34(4): 852-863. doi: 10.13224/j.cnki.jasp.2019.04.015
摘要:
在细观尺度下,基于逐渐损伤理论,建立了一种三维四向C/C复合材料高温强度预测模型。模型考虑了纤维束挤压后的截面形状、单胞的周期性以及纤维束和基体的脱黏等因素,引入考虑温度的三维Hashin失效准则进行单元的失效判定,预测了三维四向C/C复合材料室温和有防氧化保护700 ℃的拉伸强度。为了将模型发展到高温氧化环境,建立了考虑氧化速率的纤维束高温氧化环境力学性能退化模型,结合纤维束和单向板力学性能等价性原理,实现了无防氧化保护下三维四向C/C复合材料700 ℃拉伸强度的预测。研究了切边加工对三维四向复合材料强度的影响,建立了考虑切边宽度的切边三维四向复合材料强度预测模型,预测了有、无防氧化保护切边宽度为18 mm的三维四向C/C复合材料拉伸强度。结果表明:对非切边试验件在室温、有防氧化涂层700 ℃和无防氧化涂层700 ℃的预测误差分别为5.51%、7.20%和7.13%,拉伸过程的应力-应变曲线与试验结果吻合度较好;对切边试验件在室温和有防氧化涂层700 ℃的预测误差分别为0.88%和4.53%。多种类的算例表明预测模型合理、可靠。
热塑性定型剂对干铺丝-RTM工艺及力学性能影响
李晓涵, 李勇, 还大军, 刘丽, 汪愿
2019, 34(4): 864-874. doi: 10.13224/j.cnki.jasp.2019.04.016
摘要:
为了研究热塑性定型剂的添加对干铺丝-树脂传递模塑(RTM)工艺的影响,采用PA-6作为热塑性定型剂,通过T-剥离确定定型剂的质量分数,结合试验和工艺仿真,研究了定型剂对预制体面内渗透率的影响,并对其充模过程进行工艺仿真和分析,在此基础上研究了定型剂对环氧复合材料制品层间性能及韧性的影响。渗透率试验和仿真结果表明:加入定型剂后,纤维预制体的沿纤维和垂直纤维的面内渗透率分别降低了55.5%和52.0%;工艺仿真结果与试验结果一致,不同渗透长度下的注胶时间偏差均小于10%。力学性能试验结果表明:加入热塑性定型剂可显著提升环氧树脂复合材料的韧性,Ⅰ型和Ⅱ型层间断裂韧性分别提高了89.3%和54.1%,冲击韧性提高了18.7%,复合材料的层间剪切强度降低了23.6%。
环境氧化物对EB-PVD热障涂层残余应力的影响
陈智, 贾文斌, 赵凯, 黄红梅, 方磊
2019, 34(4): 875-884. doi: 10.13224/j.cnki.jasp.2019.04.017
摘要:
通过建立EB-PVD涂层的柱状晶模型,系统地研究了CMAS沉积和渗透对残余应力的影响。结果表明:CMAS的沉积和渗透产生了很大的面内、面外拉应力,导致垂直、水平裂纹的萌生;随着裂纹扩展,最终涂层发生剥落和分层失效。该失效主要发生在TC层表面下方、CMAS渗透交界处、陶瓷/黏结(TC/BC)界面附近。
Z-pin增强复合材料开孔层合板压缩性能
孙一凡, 李勇, 还大军, 刘洪全, 杨潇
2019, 34(4): 885-893. doi: 10.13224/j.cnki.jasp.2019.04.018
摘要:
针对复合材料开孔层合板孔边应力集中而导致承载能力降低现象,采用Z-pin增强技术提高性能。重点了研究开孔层合板的压缩性能,制备了不同Z-pin植入体积分数的开孔层合板试样并进行了压缩试验,在试验基础上,应用有限元分析方法对Z-pin增强开孔层合板进行建模分析,基于渐进损伤分析法和内聚单元法模拟了复合材料基体以及Z-pin的增强效果。试验结果表明:Z-pin的植入可以显著提高开孔层合板的抗压性能,且随着Z-pin植入体积分数的提高,压缩性能不断上升,压缩强度最高可以提高23.06%;只有沿载荷方向,位于开孔两侧的Z-pin可以提高结构的承载能力;在Z-pin桥联力的作用下,降低了损伤的扩展速度,但并没有改变层合板的最终失效模式。模拟结果表明:压缩损伤从孔边向两侧扩展,位于层合板中心的铺层最先出现损伤,其中第6层与第10层的增强单元最先发生破坏。模拟计算结果与试验结果得到了较好的吻合,该有限元模型对于Z-pin增强复合材料开孔层合板的适用性得到了验证。
燃烧、传热、传质
OH与CH2O双组分平面激光诱导荧光对旋流燃烧室火焰结构与脉动特征的研究
严浩, 张少华, 余西龙, 李飞, 林鑫
2019, 34(4): 894-907. doi: 10.13224/j.cnki.jasp.2019.04.019
摘要:
介绍了利用平面激光诱导荧光(PLIF)技术对航空发动机的旋流燃烧室模型在贫燃状态下工作特性的研究。通过对OH与CH2O双组分进行同步PLIF测量,获得了不同工况下燃烧室反应区以及预热区的瞬态结构信息。应用本征正交分解(POD)方法对OH PLIF的图像进行处理,得到了旋流火焰的主要脉动模态,并通过扩展本征正交分解(EPOD)方法计算出了相应POD模态的CH2O荧光信号分布。实验结果表明:随着燃烧室热功率的增大,火焰的整体结构、脉动模式均出现了明显的变化。在火焰高度增加的同时,轴向不稳定性逐渐增强,涡核旋进(PVC)的脉动特征相对减弱。在较大的热功率下,在燃烧室的外回流区(ERZ)出现未燃烧的燃料。
超声速内流道摩擦阻力分析及减阻技术研究
王帅, 何国强, 秦飞, 魏祥庚
2019, 34(4): 908-919. doi: 10.13224/j.cnki.jasp.2019.04.020
摘要:
针对带有后向台阶的等截面受限空间,通过三维数值模拟开展了超声速内流道摩擦阻力分析及减阻技术研究。分析对比了飞行马赫数为5、6、6.5及7对应的燃烧室入口条件下相同质量氢气喷注、燃烧对壁面摩擦阻力的影响机制以及不同喷注压力对喷孔下游壁面剪应力的影响。研究结果表明,同等质量的氢气,低速喷注优于高速喷注(507、50.7 kPa喷注压力分别得到10%、5%左右的减阻效果)。近壁区燃烧得到接近70%的减阻效果;气流经过突扩结构之后,壁面剪应力呈现规律地不均匀变化,最大差异达100%;剪应力与密度变化趋势基本吻合。因此,发动机内流道减阻的关键在于营造近壁区低密度场;稳定、有效的减阻区域发生在靠后方的位置,但由于流动掺混、燃料的燃烧消耗,减阻效果沿流向逐渐减弱。
基于CH自由基发光的旋流火焰放热率时空特性研究
王宽亮, 李飞, 刘德峰, 余西龙
2019, 34(4): 920-928. doi: 10.13224/j.cnki.jasp.2019.04.021
摘要:
对不同工况下CH4/air旋流火焰的放热率在时间上的热声振荡现象和空间的三维形态转变两方面进行了研究。在燃烧形态转变方面,由于旋流火焰的复杂流场分布特性,采用基于化学自发光的三维计算层析技术(3D-CTC),测量了雷诺数从5 000到20 000的三个工况下旋流燃烧的CH*发光三维火焰结构。以此表征放热率的三维分布,实现对旋流火焰放热空间形态的测量。该诊断方法通过对旋流火焰发光在8个视角下的二维成像,结合层析重建算法得到其三维CH*分布信息。为验证重建保真度,将重建后结果二维可视化与高速摄影下的二维时均结果进行对比,结果表明重建误差在5%以内。研究中,分析了不同雷诺数下放热率的空间变化规律,结果显示所有实验工况下放热率的垂直于喷嘴方向的变化程度比沿喷嘴轴向的要剧烈;而随着雷诺数增加,最大的放热区表现出了明显的向后推进趋势。在旋流燃烧的热声振荡方面,利用CH*的二维高速摄影,对旋流燃烧的放热率不稳定性进行研究,发现放热率的振荡频率随着雷诺数的增大逐渐增加。
火箭发动机
液体火箭发动机非稳态燃烧过程对其稳定性的影响
汪广旭, 谭永华, 陈建华, 陈宏玉
2019, 34(4): 929-936. doi: 10.13224/j.cnki.jasp.2019.04.022
摘要:
基于模型分析了由非稳态燃烧过程引起的不稳定热释放这一激励源项对稳定性的影响,探讨了相应的激励机理。采用小扰动分析方法,将一维层流预混火焰燃烧模型方程分解为稳态和非稳态两类控制方程,建立了线性增长率中燃烧释热项的计算方法,获得了温度及释热量的空间分布,探讨了初始温度及比热比变化对其影响规律。结果表明:采用高能推进剂的燃烧系统具有相对较低的稳定性,而对于同一种推进剂,比热比的增加会导致稳定性对频率的敏感程度增加,且具有明显的频率选择性。
旋流燃烧室对固体燃料冲压发动机药柱表面传热以及燃速的影响
李唯暄, 陈雄, 周长省, MUSA Omer
2019, 34(4): 937-949. doi: 10.13224/j.cnki.jasp.2019.04.023
摘要:
为研究旋流燃烧室对固体燃料冲压发动机(SFRJ)药柱表面传热以及燃速的影响,以高密度聚乙烯(HDPE)为燃料,对旋流和无旋工况下的固体燃料冲压发动机进行了连管实验研究,并且编制了二维轴对称湍流燃烧仿真程序,采用流固耦合传热的方法以及非定常时间推进方式,对实验工况进行了数值模拟。结果表明:①药柱表面热流密度对燃速有显著影响,在回流区与附着点处,药柱表面的对流换热能力要明显优于再发展区;②在旋流工况下,在离心力与切向速度的作用下,使热解产物在药柱表面附近区域停留时间更长,有助于热解产物的充分反应,并且明显增强药柱表面对流换热能力,与无旋工况相比,提高幅度可达100%,并且在旋流工况下发动机可更快建立自持燃烧;③通过实验研究发现,旋流的引入提高发动机的燃速有积极作用,增幅可达26%,但会导致固体燃料冲压发动机补燃室压强出现周期性振荡。
叶轮机械
双调和方程在航空离心泵叶型优化设计的应用
刘显为, 李华聪, 史新兴, 符江锋
2019, 34(4): 950-960. doi: 10.13224/j.cnki.jasp.2019.04.024
摘要:
为了减少构建代理模型的计算量,提高优化效率,基于双调和方程模型和混合人工鱼群算法对航空离心泵叶型进行优化改型设计。在Matlab平台下分别采用5点4次Bezier曲线和线性函数控制叶型圆周角分布和积叠变化规律,通过软件UG和Fluent的联合批处理方法对15组设计结果进行内流场数值仿真。在曲面插值获得计算空间边界条件的基础上,采用中心差分格式对双调和偏微分方程进行数值求解,建立超曲面性能代理模型。以效率最高为目标函数基于人工鱼群算法对设计变量进行全局寻优。结果表明:基于双调和方程的超曲面代理模型能够保证试验值与预测值完全一致,而优化后的离心叶轮流动中的尾迹效应被明显减弱,水力效率比原型泵提高5.4%。