[go: up one dir, main page]

CN108388747A - 安装角误差叶片的多通道周向类正弦分布样本实现方法 - Google Patents

安装角误差叶片的多通道周向类正弦分布样本实现方法 Download PDF

Info

Publication number
CN108388747A
CN108388747A CN201810198451.9A CN201810198451A CN108388747A CN 108388747 A CN108388747 A CN 108388747A CN 201810198451 A CN201810198451 A CN 201810198451A CN 108388747 A CN108388747 A CN 108388747A
Authority
CN
China
Prior art keywords
error
blade
fixed angles
distribution
circumferential
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.)
Granted
Application number
CN201810198451.9A
Other languages
English (en)
Other versions
CN108388747B (zh
Inventor
郑似玉
滕金芳
羌晓青
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Jiao Tong University
Original Assignee
Shanghai Jiao Tong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shanghai Jiao Tong University filed Critical Shanghai Jiao Tong University
Priority to CN201810198451.9A priority Critical patent/CN108388747B/zh
Publication of CN108388747A publication Critical patent/CN108388747A/zh
Application granted granted Critical
Publication of CN108388747B publication Critical patent/CN108388747B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

一种安装角误差叶片的多通道周向类正弦分布样本实现方法,通过设置叶片安装角公差范围及高斯分布概率密度函数,随机生成安装角误差并得到多通道模拟样本,经过重排后得到周向类正弦分布样本。本发明将正弦曲线与叶片安装角误差的周向变化趋势相结合,以减小安装角误差造成的压气机性能的恶化程度,在适应叶轮机械高性能的发展要求的同时控制压气机叶片的加工成本。

Description

安装角误差叶片的多通道周向类正弦分布样本实现方法
技术领域
本发明涉及的是一种叶轮机械领域的技术,具体是一种安装角误差叶片的多通道周向类正弦分布样本实现方法。
背景技术
在航空发动机轴流高压压气机中,叶片是实现能量转换的零件,也是数量最多、加工最复杂的零件,叶片的加工质量对目前高负荷设计的压气机效率和气动稳定性起着至关重要的作用。然而由于叶片加工的复杂性,实际加工得到的几何极易偏离设计参数产生误差,尤其是压气机叶片的安装角误差。大量的安装角误差叶片将会导致压气机气动性能的急剧下降,尤其是喘振裕度的剧烈减小。
在实际加工中,通过尽可能地缩小规定公差范围,可以有效地减小叶片加工误差的大小,从而减少性能损失,但这种方式会导致加工成本的显著增加。若能在叶片误差范围不变的前提下,通过其他方法减少压气机的性能损失,将给压气机叶片的制造领域带来重大的意义。
发明内容
本发明针对现有技术存在的上述不足,提出一种安装角误差叶片的多通道周向类正弦分布样本实现方法,将正弦曲线与叶片安装角误差的周向变化趋势相结合,以减小安装角误差造成的压气机性能的恶化程度,在适应叶轮机械高性能的发展要求的同时控制压气机叶片的加工成本。
本发明是通过以下技术方案实现的:
本发明通过设置叶片安装角公差范围及高斯分布概率密度函数(PDF),随机生成安装角误差并得到多通道模拟样本,经过重排后得到周向类正弦分布样本。
所述的叶片安装角公差范围是指:参考航标规定的各级准则,确定压气机叶片的安装角公差范围采用但不限于[-1.5°,+1.5°]。
所述的高斯分布概率密度函数,在选定的公差范围下,建立其中:μ为几何误差的平均值,取为0,σ为标准差,由公差范围而定;使几何偏差的分布满足:偏离原型越远,其概率密度越低,且到达公差边界时概率密度接近于0的状态。
所述的安装角误差,通过将公差范围进行合理划分,构建出多个小区间,基于概率密度函数所对应的累积分布函数计算各个小区间的发生概率;确定所需叶片制造的总基数,由得到的各个小区间的发生概率值分别计算各个小区间内的实际叶片数,然后在这些小区间内采用随机函数构建一组随机数,随机数的数量为该区间内叶片的数量,这些随机数即为各个叶片的安装角误差。
所述的多通道模拟样本,通过确定压气机叶片的计算周期(相当于确定计算的通道数)。将所有叶片的安装角误差混合后乱序随机排列并分组,每组的叶片数即为计算通道数,由此得到大量的安装角误差叶片的随机组合,即为安装角误差叶片的多通道模拟样本。
所述的周向类正弦分布方式是指:通过重新分配多通道模拟样本中各个叶片的安装角误差以正弦曲线为参照对象进行分布重排,选定正弦曲线方程为y=Asin[α(x-1)],其中x和y分别为叶片周向排布序号和安装角误差值,振幅A为最初选定的公差范围边界,α则是取决于压气机叶片总数和计算周期的相关系数,则此方程即为基于更新误差叶片的周向分布顺序而得到的周向类正弦分布样本。
技术效果
与现有技术相比,本发明与实际的叶片加工和安装情况相吻合,为研究叶片的安装角误差对压气机性能的影响规律提供了更真实的数据结果。同时,这种沿周向类正弦的连续性安装角误差叶片分布方式,可以有效减小相邻叶片间安装角误差大梯度变化的情况,降低相邻通道流场不均匀变化以及相邻叶片负载不匹配情况的发生概率,使得相邻通道间流场变化的程度极为细微,与设计状态更加一致,从而减少气动性能相对于原型的恶化程度。
附图说明
图1为本发明压气机转子示意图;
图2为实施例的某压气机转子叶片安装角定义及其产生偏差时的几何示意图;
图中:γ表示安装角,c表示弦长,数字为安装角误差大小示例;
图3为实施例在某一公差范围上的高斯分布概率密度函数(PDF)曲线;
图4为本发明的某一多通道模拟样本及其对应的类正弦分布模式示意图。
具体实施方式
如图1所示,本实施例具体包括以下步骤:
步骤1.叶片安装角公差范围及PDF的选取
1.1首先,对于本压气机叶片,选定其安装角公差范围为[-1.5°,+1.5°]。
1.2在选定的公差范围下,建立合理的高斯分布概率密度函数,使几何分布满足:偏离原型越远,其概率密度越低,且到达公差边界时概率密度接近于0的状态,因此本例的高斯分布概率密度函数中,标准差σ取为0.5,则函数图像如图3所示。
步骤2.安装角误差的随机生成
2.1将公差范围[-1.5°,+1.5°]进行合理划分,构建出多个小区间,此实例中划分为[-1.5°,-0.5°],[-0.5°,+0.5°],[+0.5°,+1.5°]三个区间进行后续计算。
2.2基于步骤1中的概率密度函数所对应的累积分布函数,计算得到这三个小区间的发生概率分别为0.16,0.68,0.16。
2.3同时,确定所需叶片制造的总基数,在此例中,为方便说明,假设叶片制造总数为100片,则由上述概率值,得到三个区间内的实际叶片数分别为16,68,16。
2.4随后,采用随机函数,在这三个小区间内,分别构建一组相应数量的随机数,即在[-1.5°,-0.5°]和[+0.5°,+1.5°]区间内分别生成16个随机数,[-0.5°,+0.5°]区间内则生成68个随机数。这些随机数即代表着各个叶片的安装角误差大小。
步骤3.多通道模拟样本的生成
3.1确定压气机叶片的计算周期(相当于确定计算的通道数)。在本实施例中,转子叶片总数为64,根据实际的数值计算能力,本实例中选择1/4周期为计算周期,即计算样本的通道数选为16通道。
3.2将步骤2中得到的100个随机误差叶片混合并乱序,进行随机排列并分组,每组叶片数为计算通道数,因此为16片,由此可以得到至少6组的安装角误差叶片随机组合,这些组合即为安装角误差叶片的多通道模拟样本。
3.3若需要更多的几何样本,只需增加叶片总基数,并重复2.3,2.4,及3.2的步骤即可。
步骤4.周向类正弦分布方式的实现
4.1在步骤3得到的多通道模拟样本中,沿周向的安装角误差叶片均为随机分布,现重新分配其安装的周向排列顺序,将各个叶片的安装角误差大小以正弦曲线为参照对象进行分布。选定正弦曲线方程为y=Asin[α(x-1)],振幅A为最初选定的公差范围边界,因此本实例中为1.5,α则是取决于压气机叶片总数(64)和其计算周期(1/4周期)的相关系数,因此本实例中为π/8,使得在这16通道的模拟中,叶片误差的分布至少满足正弦函数的1个周期(如图4所示)。
4.2以方程为基准(见图4黑色曲线)更新误差叶片的周向分布顺序而得到的多通道样本即为类正弦分布样本。选取3.2中的某一随机几何样本为例(见图4红色空心柱状图),将其每个叶片的顺序按照各自的误差大小进行重新排列,得到图4黑色实心柱状图,其沿周向的变化趋势基本符合上述正弦曲线,即为此随机样本的类正弦分布样本,因此将此种分布方法称为“类正弦分布法”。
本发明的效果评价指标
本发明的步骤4对最终效果给出了最大贡献,使得总体性能的恶化程度显著减少。在本实施例中,以等熵效率、喘振裕度及相对总压损失系数为性能指标进行本方法的评估,对比类正弦分布模式、随机分布模式与原型的性能结果,如表1所示。由表1可见,相对于原型,安装角误差造成了叶片性能的恶化,但是,通过将叶片控制为周向类正弦的分布模式,三种总体性能指标相对于随机分布的结果均有极为显著的提升。
表1各性能指标对比
与现有技术相比,本发明将正弦曲线与叶片安装角误差的周向变化趋势相结合,极大地减少了压气机性能的恶化程度。在压气机制造工业中,若将本发明应用于实际叶片的加工安装过程,即,在实际叶片安装时,施加类正弦的周向分布方法,那么在不缩减公差范围甚至轻微增大公差范围的情况下,可以有效减少安装角误差对压气机性能的恶化作用,缩减制造成本,也满足性能要求。
上述具体实施可由本领域技术人员在不背离本发明原理和宗旨的前提下以不同的方式对其进行局部调整,本发明的保护范围以权利要求书为准且不由上述具体实施所限,在其范围内的各个实现方案均受本发明之约束。

Claims (5)

1.一种安装角误差叶片的多通道周向类正弦分布样本实现方法,其特征在于,通过设置叶片安装角公差范围及高斯分布概率密度函数,随机生成安装角误差并得到多通道模拟样本,经过重排后得到周向类正弦分布样本。
2.根据权利要求1所述的方法,其特征是,所述的高斯分布概率密度函数,在选定的公差范围下,建立其中:μ为几何误差的平均值,取为0,σ为标准差。
3.根据权利要求1所述的方法,其特征是,所述的安装角误差,通过将公差范围进行合理划分,构建出多个小区间,基于概率密度函数所对应的累积分布函数计算各个小区间的发生概率;确定所需叶片制造的总基数,由得到的各个小区间的发生概率值分别计算各个小区间内的实际叶片数,然后在这些小区间内采用随机函数构建一组随机数,随机数的数量为该区间内叶片的数量,这些随机数即为各个叶片的安装角误差。
4.根据权利要求1所述的方法,其特征是,所述的多通道模拟样本,通过确定压气机叶片的计算周期;将所有叶片的安装角误差混合后乱序随机排列并分组,每组的叶片数即为计算通道数,由此得到大量的安装角误差叶片的随机组合,即为安装角误差叶片的多通道模拟样本。
5.根据权利要求1所述的方法,其特征是,所述的周向类正弦分布方式是指:通过重新分配多通道模拟样本中各个叶片的安装角误差以正弦曲线为参照对象进行分布重排,选定正弦曲线方程为y=Asin[α(x-1)],其中x和y分别为叶片周向排布序号和安装角误差值,振幅A为最初选定的公差范围边界,α则是取决于压气机叶片总数和计算周期的相关系数,则此方程即为基准更新误差叶片的周向分布顺序而得到的周向类正弦分布样本。
CN201810198451.9A 2018-03-12 2018-03-12 安装角误差叶片的多通道周向类正弦分布样本实现方法 Active CN108388747B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810198451.9A CN108388747B (zh) 2018-03-12 2018-03-12 安装角误差叶片的多通道周向类正弦分布样本实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810198451.9A CN108388747B (zh) 2018-03-12 2018-03-12 安装角误差叶片的多通道周向类正弦分布样本实现方法

Publications (2)

Publication Number Publication Date
CN108388747A true CN108388747A (zh) 2018-08-10
CN108388747B CN108388747B (zh) 2019-08-02

Family

ID=63067376

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810198451.9A Active CN108388747B (zh) 2018-03-12 2018-03-12 安装角误差叶片的多通道周向类正弦分布样本实现方法

Country Status (1)

Country Link
CN (1) CN108388747B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110059414A (zh) * 2019-04-22 2019-07-26 北京理工大学 一种直接控制通道的二维叶片造型方法
CN111241653A (zh) * 2019-04-19 2020-06-05 上海交通大学 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103649528A (zh) * 2011-05-19 2014-03-19 米塔科技有限公司 风力涡轮机和风力涡轮机偏航角控制的方法
US8718975B2 (en) * 2009-03-02 2014-05-06 Rolls-Royce, Plc Surface profile evaluation
CN107489651A (zh) * 2017-10-10 2017-12-19 北京航空航天大学 一种基于二次函数的可抑制风扇激波噪声的叶型优化方法
CN108804853A (zh) * 2018-06-28 2018-11-13 东北大学 基于变截面梁的弹性支承下扭形凸肩叶片动力学建模方法
CN109340042A (zh) * 2018-11-19 2019-02-15 内蒙古工业大学 涡轮机叶片及其设计方法和太阳能烟囱发电系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8718975B2 (en) * 2009-03-02 2014-05-06 Rolls-Royce, Plc Surface profile evaluation
CN103649528A (zh) * 2011-05-19 2014-03-19 米塔科技有限公司 风力涡轮机和风力涡轮机偏航角控制的方法
CN107489651A (zh) * 2017-10-10 2017-12-19 北京航空航天大学 一种基于二次函数的可抑制风扇激波噪声的叶型优化方法
CN108804853A (zh) * 2018-06-28 2018-11-13 东北大学 基于变截面梁的弹性支承下扭形凸肩叶片动力学建模方法
CN109340042A (zh) * 2018-11-19 2019-02-15 内蒙古工业大学 涡轮机叶片及其设计方法和太阳能烟囱发电系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郑似玉 等: "叶片加工超差对高压压气机性能影响和敏感性分析", 《机械工程学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111241653A (zh) * 2019-04-19 2020-06-05 上海交通大学 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法
CN111241653B (zh) * 2019-04-19 2024-03-19 上海交通大学 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法
CN110059414A (zh) * 2019-04-22 2019-07-26 北京理工大学 一种直接控制通道的二维叶片造型方法

Also Published As

Publication number Publication date
CN108388747B (zh) 2019-08-02

Similar Documents

Publication Publication Date Title
CN112528575B (zh) 一种船舶燃气轮机多级轴流压气机多维度协同设计方法
RU2666933C1 (ru) Деталь или узел газотурбинного двигателя и соответствующий газотурбинный двигатель
JP2009144716A (ja) ターボ機械用の多段タービンを設計する方法
CN109598081B (zh) 基于数据降维及多二维流面的径流式透平气动优化方法
CN109356666A (zh) 一种轴流式涡轮大小叶片组合叶栅的叶型设计方法
US20160245299A1 (en) Turbomachine part with a non-axisymmetric surface
CN105465037B (zh) 一种双吸离心泵叶轮的水力优化方法以及装置
CN108108528B (zh) 一种分轴式航改燃机动力涡轮一维匹配设计方法
CN108388747B (zh) 安装角误差叶片的多通道周向类正弦分布样本实现方法
CN112287580A (zh) 一种基于全三维数值模拟的轴流压气机喘振边界计算方法
CN117764771A (zh) 基于k-shape聚类海上风电出力预测修正方法系统
CN113094964A (zh) 一种生成叶片加工坐标的方法和装置
Goryachkin et al. Three-stage low pressure compressor modernization by means of optimization methods
CN107956710A (zh) 基于级间匹配效应的立式多级泵叶轮水力设计方法
Zheng et al. Impact of nonuniform stagger angle distribution on high-pressure compressor rotor performance
US9482237B1 (en) Method of designing a multi-stage turbomachine compressor
Lyes et al. Experimental evaluation of the high-to-low speed transformation process for a highly loaded core compressor stage
CN113361029B (zh) 枞树形榫接结构配合间隙的计算方法、装置、设备及介质
CN111241653A (zh) 基于蒙特卡洛法的多维度叶片误差性能影响仿真实现方法
Wang et al. Adjoint aerodynamic design optimization for blades in multi-stage turbomachines: part ii—validation and application
Reitz et al. Performance analysis of deteriorated high pressure compressor blades
CN114493211A (zh) 一种内燃机车的排放量评估方法及装置
CN114282323A (zh) 一种涡轮叶片层板冷却结构的流量分配预测方法
Shao et al. Aerodynamic optimization of the radial inflow turbine for a 100kW-class micro gas turbine based on metamodel-semi-assisted method
Zheng et al. The impact of uncertain stagger angle variation on high-pressure compressor rotor performance

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