CN113048086B - 一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法 - Google Patents
一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法 Download PDFInfo
- Publication number
- CN113048086B CN113048086B CN202110291616.9A CN202110291616A CN113048086B CN 113048086 B CN113048086 B CN 113048086B CN 202110291616 A CN202110291616 A CN 202110291616A CN 113048086 B CN113048086 B CN 113048086B
- Authority
- CN
- China
- Prior art keywords
- fan
- noise
- unequal
- distance
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D25/00—Pumping installations or systems
- F04D25/02—Units comprising pumps and their driving means
- F04D25/08—Units comprising pumps and their driving means the working fluid being air, e.g. for ventilation
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D17/00—Radial-flow pumps, e.g. centrifugal pumps; Helico-centrifugal pumps
- F04D17/08—Centrifugal pumps
- F04D17/16—Centrifugal pumps for displacing without appreciable compression
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/26—Rotors specially for elastic fluids
- F04D29/28—Rotors specially for elastic fluids for centrifugal or helico-centrifugal pumps for radial-flow or helico-centrifugal pumps
- F04D29/281—Rotors specially for elastic fluids for centrifugal or helico-centrifugal pumps for radial-flow or helico-centrifugal pumps for fans or blowers
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/26—Rotors specially for elastic fluids
- F04D29/28—Rotors specially for elastic fluids for centrifugal or helico-centrifugal pumps for radial-flow or helico-centrifugal pumps
- F04D29/30—Vanes
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F04—POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
- F04D—NON-POSITIVE-DISPLACEMENT PUMPS
- F04D29/00—Details, component parts, or accessories
- F04D29/66—Combating cavitation, whirls, noise, vibration or the like; Balancing
- F04D29/661—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps
- F04D29/666—Combating cavitation, whirls, noise, vibration or the like; Balancing especially adapted for elastic fluid pumps by means of rotor construction or layout, e.g. unequal distribution of blades or vanes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明公开了一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法,包括:参数化设计不等距离心风扇三维初始模型;数值分析风扇流场域,获取原始风扇气动噪声级大小及风扇性能参数;以风扇相邻叶片间距为设计变量,最优拉丁超立方方法进行样本点试验设计;选用优化分析集成平台,批量模拟分析风扇流场,风扇气动噪声与风扇性能为响应输出,采集样本数据;拟合样本分析结果,构建径向基神经网络近似模型;以风扇噪声级最小为优化目标,静压及流量为约束函数,基于近似模型,运用存档微遗传算法(AMGA)对冷却风扇相邻叶片间距进行优化设计,通过选取帕雷托(Pareto)解集最优解,得到最优不等距低噪声离心风扇优化设计方案。
Description
技术领域
本发明属于一种低噪声不等距离心风扇优化设计领域,尤其涉及的是一种基于径向基神经网络模型低噪声不等距小型风冷内燃机离心风扇优化设计方法。
背景技术
噪声是人类环境重要污染之一,长期处于高分贝噪声环境中,不仅会降低工作效率,影响情绪,还会严重影响身体机能。随着科学技术的发展很人类生活水平的提高,噪声控制问题越来越受到人们的重视,不同国家、地区和行业都制定了严格的噪声限值标准来限制噪声污染。
离心风扇是小型风冷内燃机冷却系统必不可少的一个重要部件。由飞轮带动离心风扇,通过风扇旋转运动驱动冷空气克服系统阻力,使冷空气经过导风罩、导流板等引流装置,高速吹过气缸外壁及机体散热片表面,带走散出的热量实现冷却。离心风扇高速旋转运动会产生机械振动噪声和气动噪声,气动噪声是风扇驱动空气而产生的气体扰动,主要包括旋转噪声和涡流噪声。旋转噪声是由于风扇旋转叶片周期性地打击空气质点,使空气形成周期性压力脉动而产生的一系列离散的峰值频率噪声;涡流噪声是由于风扇旋转时使周围的空气产生涡流在分裂时空气发生的扰动而激发出的噪声,是连续的宽频带噪声。
离心风扇的气动噪声相比机械噪声要显著很多,因此降噪的主要目标是降低风扇的气动噪声,特别是降低对人体危害最大的旋转噪声的大小。影响离心风扇气动噪声大小的主要因素有扇叶结构及叶片布局方式,在不改变叶片几何参数的情况下,可以通过合理安排叶片布局,改变风扇辐射声场,达到降低噪声的目的。
小型风冷内燃机上大多数采用前弯叶片及叶片等高、间距均匀的布局形式,风扇转动时,叶片通过频率(BPF)及其谐波处会产生较强的离散峰值频率噪声,这些峰值组成了风扇的旋转噪声。采用不等间距风叶周向角布置,通过分散基频与谐波峰值,可以减少这一周期激励效果,能够降低风扇基频处的噪声峰值,有效降低风扇旋转噪声,从而达到降低风扇总体噪声级的目的,如CN201410379874-一种柴油机用的冷却风扇,CN200610088548-用于气冷式内燃机多叶式风扇等已证明该方法的有效性。
不等距叶片的设计在降低风扇噪声的同时还应该满足两个条件,一是叶片不等距布置后应该仍能保证气动性能要求,或者变化不大;二是叶片布置方式改变后,仍然要能够满足风扇运转时的动平衡要求。目前不等距风扇叶片布置设计通常采用分布公式或试验方法,研制周期较长,研发费用高,且有很大的随机性,设计方法与理论研究的还比较少。最常用的不等距风扇分布方式是按公式θ′=θ+Asinnθ进行布置,n为等分组数,(当n≥2时风扇为自平衡状态,可以满足动平衡要求),A为相位调制量。选择合适的A值与分组数n可以降低谐波振幅,进而降低风扇的总体噪声,实际上由于风扇结构参数与运转工况的影响,合适的A值与分组数n很难确定,选择不当不仅不能降低风扇噪声,甚至还会破坏风扇的气动性能。
随着计算机技术的发展,特别是计算流体动力学(CFD)、计算空气动声学(CAA)、多目标优化技术等的发展,为优化叶片布置,降低风扇噪声,设计高性能不等距低噪声离心风扇,提供了有效方法与手段。
发明内容
本发明的目的是降低离心风扇噪声,改进不等距离心风扇设计方法上的不足,提供一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法。通过三维软件设计不等距离心风扇参数化三维初始模型;利用计算流体力学(CFD)及计算气动声学(CAA)技术,对风扇流场进行数值模拟分析,获取原始风扇气动噪声级大小及风扇性能参数;以风扇相邻叶片间距为设计变量,采用最优拉丁超立方方法进行样本点试验设计,选用ISIGHT优化分析集成平台,批量模拟分析风扇流场,输出风扇气动噪声级大小及风扇流量、静压等性能样本数据;依据样本结果数据,构建径向基神经网络近似模型,建立离心风扇不等距布局与风扇性能及噪声级大小的量化关系,以风扇噪声声压级最小为优化目标,静压及流量为约束函数,运用存档微遗传算法(AMGA)对离心风扇相邻叶片间距进行优化设计,通过选取Pareto解集最优解,得到最佳不等距低噪声离心风扇设计方案。
本发明集成了参数化建模技术、计算流体动力学(CFD)、计算空气动声学(CAA)、径向基神经网络模型、存档微遗传优化搜索技术等,能够快速、准确的建立离心风扇噪声及性能及与不等距叶片布局的关系,优化设计相邻叶片不等距角度,满足风扇性能的同时,最大限度的降低离心风扇噪声的大小。
本发明采取的技术方案是:
步骤1,确定优化方案,初始设计不等距离心风扇三维参数化数字模型。根据冷却风扇叶片总个数Z,叶片均分为N(N≥2)组,分组内相邻叶片间距初步设计按等差数列进行布置,相邻叶片最大变化范围不超过等距布局时相邻叶片角度值的±20%,过大的角度偏差会造成风扇性能的下降。
步骤2,采用ANSYS WORKBENCH软件,通过其FLUENT子模块,数值模拟不等距冷却风扇工作流场,获取风扇流场性能参数与噪声仿真数据。风扇噪声模拟分析采用大涡模拟(LES)与威廉姆-霍金斯方程(FW-H)声学模型,求解瞬态风扇压力脉动,经傅里叶变换(FFT),输出风扇噪声级大小,并进行试验验证。
步骤3,以相邻叶片间距为设计变量,风扇性能参数与噪声大小为输出响应,最优拉丁超立方方法进行样本点试验设计。为了得到较精确的近似模型,样本点至少大于3n个(n为分组内叶片个数)。
步骤4,基于ISIGHT软件搭建风扇多目标优化设计平台,集成风扇三维数字化设计软件CATIA与流场数值模拟软件ANSYS WORKBENCH,选取离心风扇分组内叶片不等间距为设计变量,创建设计风扇输入变量与风扇性能及气动噪声级大小的输出结果映射关系,批量模拟分析风扇流场,获取风扇噪声级大小,采集近似模型所需样本数据。
步骤5,根据样本数据,建立径向基神经网络近似模型,并进行误差分析。径向基函数是一种将样本点与待测点之间的欧氏距离当作自变量的函数,径向基函数近似模型是把径向基函数当作基函数,然后通过线性拟合构建而来的近似模型。基本表达式如下
式中:wi是基函数的线性叠加权系数;Hi(r)是径向基函数;ns是径向基个数;r是样本点与待测点之间的欧式距离。
径向基神经网络近似模型采用插值方法,以样本点为中心,通过选用不同的基函数及线性叠加权系数wi,能够最终得到在预测域的响应的近似值,样本点处误差为零,具有将多维问题转化为一维问题的能力,并且转化后预测点与样本点之间的欧氏距离是一维问题自变量,拟合复杂非线性模型能力强、鲁棒性好、学习速度快、容错功能强等优点。
近似模型构建完成后必须对近似模型的精度进行验证,以保证近似模型的可信度。常用的误差评价指标有:平均绝对值误差(简称AMAE)、最大绝对值误差(简称MAE)、均方根误差(简称RMSE)与决定系数R2,各误差评价指标计算公式如下:
式中:ntest——样本点数目;
yi——第i个样本点的响应值;
步骤6,运用存档微遗传算法(AMGA),以风扇噪声级最小为优化目标,静压及流量为约束函数,对不等距冷却风扇叶片间距进行优化分析,获得Pareto非劣解的解集。在AMGA微遗传算法进化过程中,通过Pareto最优关系将群体中的每个个体同目标函数的向量进行比较,最终将每个个体分为多个Pareto前沿层,运用Pareto优越性评价法对每层的个体进行评价,其中在Pareto同层中拥有更大拥挤距离的个体更优异,进而获得Pareto非劣解的解集。
步骤7,选取Pareto最优解,获取最优不等距风扇结构模型,对优化结果进行试验验证。
本发明的有益效果:
(1)采用径向基神经网络模型建立离心风扇不等距叶片布置与风扇性能及噪声之间近似模型,能够精确的反映出叶片布置与性能及噪声之间的关系,为不等距低噪声离心风扇设计提供一定的理论依据。
(2)采用径向基神经网络模型近似模型进行不等距风扇设计,为改善离心风扇性能,降低离心风扇噪声提供一种可行的设计方法。
(3)采用径向基神经网络模型近似模型进行不等距风扇设计,克服了目前不等距风扇设计的局限性,避免了大量的试算与试验工作,能缩短设计周期,减少研发成本。
附图说明
图1为本发明一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法流程图;
图2为小型发动机冷却风扇原机等距叶片布局几何结构示意图;
图3为小型发动机冷却风扇径向基神经网络模型误差适应度分析示意图;
图4为存档微遗传算法(AMGA)计算流程图;
图5为小型发动机冷却风扇优化设计几何结构示意图;
具体实施方式
下面结合附图以及具体实施对本发明作进一步的说明,但本发明的保护范围并不限于此。
一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法,以一台小型风冷汽油机不等距离心风扇低噪声设计为例,具体实施流程图如图1所示,具体实施步骤包括:
步骤1,利用CATIA三维软件初始设计不等距离心风扇三维参数化数字模型,建立不等距离心风扇气动噪声数值分析流场域三维模型。原等距离心风扇如图2所示,包括叶盘,27个叶片及定位柱等结构。不等距叶片初始布局将叶片均分为3组,相邻叶片间距参数分布定义为ω1至ω8,分组内以相邻叶片角度最大变化值不超过等距布局时±10%为原则,以等差递减方式不等距布置,初始不等距离心风扇叶间距分别定义为ω1=17°、ω2=16°、ω3=15°、ω4=14°、ω5=13°、ω6=12°、ω7=11°、ω8=10°。绘制小型风冷内燃机导风罩、飞轮、主轴箱、缸盖、导风板、启动杯等离心风扇主要关联部件及试验测试风筒三维模型,根据各关联部件在整机中的位置关系进行装配,提取由离心风扇、离心风扇关联部件与测试风筒共同包围的区域,建立初始不等距离心风扇气动噪声分析流场域三维模型,并保存成*.Catpart文件。
步骤2,采用ANSYS Workbench有限元综合分析软件,将不等距离心风扇气动噪声数值分析流场域三维模型文件导入到Geometry子模块,然后通过Mesh子模块对流场域模型进行网格划分,网格划分后转入到Fluent子模块。在该模块下选取瞬态分析模式,大涡模拟(LES)湍流求解模型;开启廉姆-霍金斯方程(FW-H)声学模型;设定离心风扇旋转区域为动网格形式,设置压力进出口边界条件,风扇转速设定为发动机额定转速3600转/分钟;以进口流量、进出口静压、风扇扭矩系数为参数化输出结果数据;求解方法采用PISO方法,二阶迎风式,数值模拟初始不等距设计离心风扇工作流场,获取离心风扇流场进口流量、进出口静压、风扇扭矩系数等性能参数及瞬态风扇压力脉动数据;仿真计算结果不等距初始设计风扇静压Pst=204.0Pa、流量Q=0.090782Kg/s、风扇效率ηst=16.87%;同方法亦可模拟计算原等距风扇风扇静压Pst=205.8Pa、流量Q=0.091726Kg/s、风扇效率ηst=16.99%;在结果分析选项中选择傅里叶变换(FFT),按照往复式内燃机声压法声功率级的测定标准,设置噪声测试观测点,对瞬态风扇压力脉动结果进行换算,得到不等距初始风扇观测点产生的气动脉冲噪声级大小最大值为76.8dB,原等距风扇噪声级在观测点出最大值为78.1dB,并将模拟分析结果保存成*.wbpj文件格式。快速成型制作初始不等距离心风扇,经风扇性能与噪声综合试验平台验证,在转速3600转/分钟发动机额定转速下倒拖试验中,实测噪声观察点位置采用不等距初始设计风扇整机噪声最大值78.3dB,采用原等距风扇整机噪声最大值79.8dB,去除发动机机械噪声等影响因素,风扇噪声模拟分析精度满足要求。
步骤3,基于ISIGHT软件搭建不等距离心风扇优化集成平台,选用最优拉丁超立方方法进行样本点试验设计,批量求解近似模型样本点数据。在ISIGHT软件中选择试验设计(DOE)模块,依次添加CATIA、ANSYS Workbench、Caculator等应用组件,集成CATIA、ANSYSWorkbench等应用程序。利用CATIA组件解析不等距离心风扇流场域三维模型*.Catpart文件,将不等距风扇离心风扇分组内叶片间距设置为输入变量,相邻叶片间距变量从大到小,分别用ω1至ω8表示;由ANSYS Workbench组件解析不等距离心风扇气动噪声分析*.wbpj文件,将流场模拟分析结果数据风扇流量、进出口静压、风扇扭矩系数、噪声级大小等设置为输出变量;将CATIA中不等距离心风扇流场域三维模型存档文件*.Catpart指定为ANSYSWorkbench中风扇气动噪声流场分析的几何模型输入文件,建立设计不等距离心风扇输入变量与风扇气动噪声分析输出结果之间的映射关系。在试验设计(DOE)模块中,选用最优拉丁超立方方法进行样本点试验设计,样本点采集个数为30个,相邻叶片间距参数ω1至ω8为设计因子,初始不等间距为基准线,因子以基准线±20%变化,试验设计样本矩阵如表1所示;由流场模拟分析结果风扇流量、进出口静压、风扇扭矩、噪声级大小经Caculator应用组件计算处理后的结果数据风扇流量Q、风扇静压Pst、风扇效率ηst、风扇噪声级dB等变量为响应参数;建立设计风扇输入变量与风扇性能及气动噪声级大小的输出结果之间映射关系,批量集成分析风扇流场,获取不同试验设计样本矩阵点风扇噪声级大小与风扇性能参数。
表1.试验设计样本矩阵
步骤4,将样本点输出结果数据如表2所示保存成文本格式,在ISIGHT集成软件中选择优化分析(Optimization)模块,将近似模型Approximation组件添加到该组件中,选择近似模型建模方法为径向基神经网络(RBF)近似模型,载入样本点数据,选择径向基函数,以待测点与样本点之间的欧几里德距离为自变量,以样本点为中心进行插值,建立径向基近似模型。近似模型构建完成后必须对近似模型的精度进行验证,以保证近似模型的可信度,分析表明风扇静压Pst、流量Q、风扇噪声级dB、风扇效率ηst等响应参数,平均绝对值误差(简称AMAE)分别为0.05701、0.0585、0.03374、0.05662;最大绝对值误差(简称MAE)分别为0.17641、0.15535、0.10183、0.17029;均方根误差(简称RMSE)分别为0.06754、0.06812、0.04163、0.06729;决定系数R2分别为0.94154、0.94299、0.97697、0.94274,四项误差指标均满足要求,近似模型响应适应度如图3所示,可以看出样本点响应值均在对角线附近,拟合精度较高。
表2.样本点结果数据
步骤5,径向基模型创建之后,在优化分析(Optimization)模块中,选取如图4所示的存档微遗传算法(AMGA),以风扇噪声级dB最小、流量Q最大、风扇效率ηst最大为目标,扇静压Pst不低于200Pa为限制条件,对不等距冷却风扇叶片间距ω1至ω8进行寻优探索分析。设置交叉概率为0.9,变异概率0.05,初始种群大小40,父种群大小40,档案种群大小为5000,最大评价次数为5000,AMGA算法的计算流程如图5所示,经5001次循环探索后最终得到优化目标函数的Pareto最优解,如表3所示。
表3.Pareto最优解集
步骤6,综合考虑选取表3中循环数为4702的一组Pareto最优解,生成最终冷却风扇三维结构优化设计模型,如图5所示,该处风扇噪声级大小近似模型输出数值为73.3dB,较优化前下降了3.5dB,较原等距风扇噪声级下降4.8dB,风扇流量、静压、功率等性能参数基本维持不变,进一步通过快速成型制作不等距离心风扇,经风扇性能测试平台测试,在转速3600转/分钟发动机额定转速下倒拖试验中,实测噪声观察点位置不等距初始设计风扇整机噪声最大值74.2dB,风扇静压Pst=201.2Pa、流量Q=0.0906Kg/s、风扇静压效率ηst=16.12%,与近似模型优化结果基本一致;在转速3600转/分钟额定工况下,通过发动机台架试验测试,离心风扇优化前小型风冷汽油机整机在额定转速额定载荷下噪声级92.3dB,采用不等距叶片布局优化后实测整机噪声值90.2dB,降低了2.1dB,且发动机机油温度基本保持不变,油耗与功率也有所改善,说明了通过本发明方法,在不改变风扇性气动性能的前提下,通过合理叶片不等距布置可以减小发动机整体噪声。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明,它们并非用以限制本发明的保护范围,凡未脱离本发明技术所创的等效方式或变更均应包含在本发明的保护范围之内。
Claims (9)
1.一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法,其特征在于,包括如下步骤:
S1,初始设计不等距离心风扇三维参数化数字模型;
S2,基于计算流体与计算声学方法数值模拟不等距离心风扇工作流场,获取风扇流场性能参数与噪声级大小;
S3,以相邻叶片间距为设计变量,风扇性能参数与噪声级大小为输出响应,最优拉丁超立方方法进行不等距叶片布局样本点设计;
S4,通过集成优化分析平台ISIGHT,根据样本点数据,建立径向基神经网络近似模型;
S5,运用存档微遗传算法AMGA,以风扇噪声级最小为优化目标,静压及流量等风扇性能参数为约束函数,对不等距冷却风扇叶片间距进行优化分析,获得Pareto最优解集;
S6,选取Pareto最优解,获取最优不等距风扇结构模型;
所述S1的具体实现包括:
S1.1,设计等距离心风扇,包括叶盘、27个叶片及定位柱;
S1.2,设计不等距叶片初始布局:将叶片均分为3组,相邻叶片间距参数分布定义为ω1至ω8,分组内以相邻叶片角度最大变化值不超过等距布局时±10%为原则,以等差递减方式不等距布置,初始不等距离心风扇叶间距分别定义为ω1=17°、ω2=16°、ω3=15°、ω4=14°、ω5=13°、ω6=12°、ω7=11°、ω8=10°;
S1.3,绘制小型风冷内燃机导风罩、飞轮、主轴箱、缸盖、导风板、启动杯等离心风扇主要关联部件及试验测试风筒三维模型,根据各关联部件在整机中的位置关系进行装配,提取由离心风扇、离心风扇关联部件与测试风筒共同包围的区域,建立初始不等距离心风扇气动噪声分析流场域三维模型。
2.根据权利要求1所述的一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法,其特征在于,所述S2的具体实现包括:
S2.1,采用ANSYS Workbench有限元综合分析软件,将不等距离心风扇气动噪声数值分析流场域三维模型文件导入到Geometry子模块,然后通过Mesh子模块对流场域模型进行网格划分,网格划分后转入到Fluent子模块,在该模块下选取瞬态分析模式,大涡模拟(LES)湍流求解模型;开启廉姆-霍金斯方程(FW-H)声学模型;设定离心风扇旋转区域为动网格形式,设置压力进出口边界条件,风扇转速设定为发动机额定转速3600转/分钟;以进口流量、进出口静压、风扇扭矩系数为参数化输出结果数据;求解方法采用PISO方法,二阶迎风式,数值模拟初始不等距设计离心风扇工作流场,获取离心风扇流场进口流量、进出口静压、风扇扭矩系数等性能参数及瞬态风扇压力脉动数据;仿真计算结果不等距初始设计风扇静压Pst=204.0Pa、流量Q=0.090782Kg/s、风扇效率ηst=16.87%;
S2.2,采用上述方法模拟计算原等距风扇风扇静压Pst=205.8Pa、流量Q=0.091726Kg/s、风扇效率ηst=16.99%;
S2.3,在结果分析选项中选择傅里叶变换(FFT),按照往复式内燃机声压法声功率级的测定标准,设置噪声测试观测点,对瞬态风扇压力脉动结果进行换算,得到不等距初始风扇观测点产生的气动脉冲噪声级大小最大值为76.8dB,原等距风扇噪声级在观测点出最大值为78.1dB,将模拟分析结果保存成*.wbpj文件格式。
3.根据权利要求1所述的一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法,其特征在于,所述S3的具体实现包括:
S3.1,在ISIGHT软件中选择试验设计模块,依次添加CATIA、ANSYS Workbench、Caculator应用组件,以集成CATIA、ANSYS Workbench应用程序;利用CATIA组件解析不等距离心风扇流场域三维模型*.Catpart文件,将不等距风扇离心风扇分组内叶片间距设置为输入变量,相邻叶片间距变量从大到小,分别用ω1至ω8表示;
S3.2由ANSYS Workbench组件解析不等距离心风扇气动噪声分析*.wbpj文件,将流场模拟分析结果数据风扇流量、进出口静压、风扇扭矩系数、噪声级大小等设置为输出变量;将CATIA中不等距离心风扇流场域三维模型存档文件*.Catpart指定为ANSYS Workbench中风扇气动噪声流场分析的几何模型输入文件,建立设计不等距离心风扇输入变量与风扇气动噪声分析输出结果之间的映射关系;
S3.3在试验设计模块中,采用最优拉丁超立方方法进行样本点试验设计,样本点采集个数为30个,相邻叶片间距参数ω1至ω8为设计因子,初始不等间距为基准线,因子以基准线±20%变化,设计的样本点如下表1所示:
表1
5.根据权利要求4所述的一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法,其特征在于,所述建立的径向基模型满足以下条件:
风扇静压Pst、流量Q、风扇噪声级dB、风扇效率ηst,四个参数的平均绝对值误差分别为0.05701、0.0585、0.03374、0.05662;最大绝对值误差分别为0.17641、0.15535、0.10183、0.17029;均方根误差分别为0.06754、0.06812、0.04163、0.06729;决定系数R2分别为0.94154、0.94299、0.97697、0.94274。
6.根据权利要求4所述的一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法,其特征在于,所述S5的具体实现包括:
在优化分析模块中,选取存档微遗传算法,以风扇噪声级dB最小、流量Q最大、风扇效率ηst最大为目标,扇静压Pst不低于200Pa为限制条件,对不等距冷却风扇叶片间距ω1至ω8进行寻优探索分析:设置交叉概率为0.9,变异概率0.05,初始种群大小40,父种群大小40,档案种群大小为5000,最大评价次数为5000,经5001次循环探索后最终得到优化目标函数的Pareto最优解集。
8.根据权利要求7所述的一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法,其特征在于,所述S6中的最优解为权利要求7中对应的循环次数在4702时对应的参数值。
9.根据权利要求8所述的一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法制作的不等距离心风扇,其特征在于,在转速3600转/分钟发动机额定转速下,不等距初始设计风扇整机噪声最大值74.2dB,风扇静压Pst=201.2Pa、流量Q=0.0906Kg/s、风扇静压效率ηst=16.12%,整机噪声值90.2dB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110291616.9A CN113048086B (zh) | 2021-03-18 | 2021-03-18 | 一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110291616.9A CN113048086B (zh) | 2021-03-18 | 2021-03-18 | 一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113048086A CN113048086A (zh) | 2021-06-29 |
CN113048086B true CN113048086B (zh) | 2022-05-20 |
Family
ID=76513657
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110291616.9A Active CN113048086B (zh) | 2021-03-18 | 2021-03-18 | 一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113048086B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113268818A (zh) * | 2021-07-19 | 2021-08-17 | 中国空气动力研究与发展中心计算空气动力研究所 | 基于生成拓扑映射的气动全局优化方法、存储介质及终端 |
CN113806977B (zh) * | 2021-08-16 | 2024-03-22 | 江铃汽车股份有限公司 | 汽车冷却风扇的噪声振动分析方法及系统 |
CN115099115B (zh) * | 2022-08-24 | 2022-12-30 | 国网江西省电力有限公司电力科学研究院 | 一种户内变电站通风降噪深度强化学习设计方法 |
CN115544815B (zh) * | 2022-11-30 | 2023-03-21 | 苏州浪潮智能科技有限公司 | 一种风扇模型的生成方法和装置 |
CN117078030B (zh) * | 2023-07-12 | 2024-05-03 | 贵州大学 | 一种基于车速预测的燃料电池公交车能量管理方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1692962A1 (fr) * | 2005-02-18 | 2006-08-23 | Faco S.A. | Sèche-cheveux à confort acoustique amelioré |
CN102072196A (zh) * | 2011-02-28 | 2011-05-25 | 上海日用-友捷汽车电气有限公司 | 冷却风扇的低噪声风叶 |
CN202659563U (zh) * | 2012-06-08 | 2013-01-09 | 上海日用-友捷汽车电气有限公司 | 高效率低噪声风叶 |
CN105117564A (zh) * | 2015-09-18 | 2015-12-02 | 中国人民解放军海军工程大学 | 一种前置定子周向非对称布置的泵喷推进器水力模型及其设计方法 |
CN109460629A (zh) * | 2018-08-30 | 2019-03-12 | 华南理工大学 | 一种基于近似模型法的冷却风扇性能优化方法 |
CN111046504A (zh) * | 2019-11-27 | 2020-04-21 | 江苏大学 | 一种适用于仿生叶片泵几何参数的低噪声优化方法 |
CN111188793A (zh) * | 2020-01-17 | 2020-05-22 | 湘潭大学 | 一种离心压气机叶轮分流叶片周向角设计方法及离心压气机叶轮 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101671946B1 (ko) * | 2014-12-04 | 2016-11-16 | 한국생산기술연구원 | 부등피치 재생 블로워 및 이의 최적화 설계 방법 |
-
2021
- 2021-03-18 CN CN202110291616.9A patent/CN113048086B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1692962A1 (fr) * | 2005-02-18 | 2006-08-23 | Faco S.A. | Sèche-cheveux à confort acoustique amelioré |
CN102072196A (zh) * | 2011-02-28 | 2011-05-25 | 上海日用-友捷汽车电气有限公司 | 冷却风扇的低噪声风叶 |
CN202659563U (zh) * | 2012-06-08 | 2013-01-09 | 上海日用-友捷汽车电气有限公司 | 高效率低噪声风叶 |
CN105117564A (zh) * | 2015-09-18 | 2015-12-02 | 中国人民解放军海军工程大学 | 一种前置定子周向非对称布置的泵喷推进器水力模型及其设计方法 |
CN109460629A (zh) * | 2018-08-30 | 2019-03-12 | 华南理工大学 | 一种基于近似模型法的冷却风扇性能优化方法 |
CN111046504A (zh) * | 2019-11-27 | 2020-04-21 | 江苏大学 | 一种适用于仿生叶片泵几何参数的低噪声优化方法 |
CN111188793A (zh) * | 2020-01-17 | 2020-05-22 | 湘潭大学 | 一种离心压气机叶轮分流叶片周向角设计方法及离心压气机叶轮 |
Also Published As
Publication number | Publication date |
---|---|
CN113048086A (zh) | 2021-06-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113048086B (zh) | 一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法 | |
CN112949224A (zh) | 一种基于响应面模型的小型风冷内燃机冷却风扇优化设计方法 | |
RU2559718C2 (ru) | Система и способ настройки/калибровки семейств ступеней турбомашин | |
CN106777783B (zh) | 一种航空发动机叶片裂纹预测方法 | |
CN111159941B (zh) | 一种用于汽车液力变矩器内流场瞬态数值模拟的方法 | |
Kersken et al. | Time-Linearized and Time-Accurate 3D RANS Methods for Aeroelastic Analysis in Turbomachinery | |
CN112417773B (zh) | 多级轴流膨胀机的多学科优化设计方法、装置及设备 | |
CN107122512B (zh) | 液环泵非稳态气体流场及吸入压缩性能的简化计算方法 | |
CN111859730A (zh) | 一种燃料电池离心压缩机转子构形优化设计方法 | |
CN114117849A (zh) | 低压涡轮叶/盘转子的叶冠阻尼减振分析方法 | |
Zhuge et al. | Development of an advanced turbocharger simulation method for cycle simulation of turbocharged internal combustion engines | |
CN117454721B (zh) | 基于数字仿真实验的风电场尾流叠加效应评估方法和介质 | |
CN106844956A (zh) | 一种航空发动机涡轮盘结构可靠性分析方法 | |
CN106599422A (zh) | 叶片泵转子系统的振动仿真分析方法及其装置 | |
Lengyel et al. | Design of a counter rotating fan-an aircraft engine technology to reduce noise and CO2-emissions | |
CN112270139A (zh) | 基于母型库的燃料电池离心压缩机气动优化设计方法 | |
Lainé et al. | Blackbox optimization for aircraft engine blades with contact interfaces | |
CN112182802A (zh) | 计入随机失谐的叶轮机械气动弹性优化设计实现方法 | |
CN106503375A (zh) | 一种基于cn群理论确定汽轮机转子临界转速的方法及系统 | |
CN116595874A (zh) | 叶轮机械性能预测模型参数优化方法及装置、存储介质 | |
CN116702359A (zh) | 脱硫循环泵关键部件结构型线优化方法和装置 | |
Fruth et al. | Influence of the Blade Count Ratio on Aerodynamic Forcing: Part II—High Pressure Transonic Turbine | |
CN113761800A (zh) | 基于临界转速对应的轴系动力学参数模型缩比设计方法 | |
Fuhrer et al. | Multi-Objective Numerical Optimization of Radial Turbines | |
Saini et al. | The effect of blade count on body force model performance for axial fans |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |