CN102339351B - 离心压气机叶轮轮盘截面形状优化方法 - Google Patents
离心压气机叶轮轮盘截面形状优化方法 Download PDFInfo
- Publication number
- CN102339351B CN102339351B CN 201110324654 CN201110324654A CN102339351B CN 102339351 B CN102339351 B CN 102339351B CN 201110324654 CN201110324654 CN 201110324654 CN 201110324654 A CN201110324654 A CN 201110324654A CN 102339351 B CN102339351 B CN 102339351B
- Authority
- CN
- China
- Prior art keywords
- sigma
- impeller
- centrifugal compressor
- design
- compressor impeller
- 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
Landscapes
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
离心压气机叶轮轮盘截面形状优化方法,对离心叶轮轮盘截面形状采用了参数化建模的方式,以最少的设计参数描述叶轮轮盘截面形状,并将神经网络近似模型和优化设计方法结合起来对叶轮轮盘截面形状进行优化设计,在保证叶轮重量不增加的情况下有效降低叶轮的等效应力达7.4%以上,提高其工作可靠性。本设计方法大幅降低叶轮强度设计周期,设计结果可靠,节约了设计成本,具有很高的社会效益及推广价值。
Description
技术领域
本发明属于机械设计领域,设计一种离心压气机叶轮轮盘截面形状优化设计的方法。
背景技术
离心压气机在国民经济与国防建设各部门中发挥着举足轻重的作用,广泛应用于航空、航天、能源、交通、运输、化工及冶金等部门。随着现代规模工业、尖端科技的发展和社会对环境保护的呼声越来越高,高负荷、大流量、高效率、低噪声以及小体积成为离心压气机的发展方向。离心叶轮作为压缩机中的关键部件,随着技术的发展,离心叶轮的设计同时向着转速更高,重量更轻的方向发展,但是随着转速的提高,重量的减轻,离心力越来越大,叶片与轮盘、轮盖之间的相互影响和相互作用也越来越大,因此如何设计出高强度的叶轮是离心压缩机设计的难点。
研究(Valakos IM,Ntipteni MS,Nikolos IK,Structuraloptimization of a centrifugal impeller using differentialevolution in CATIA environment.Operational Research-AnInternational Journal,2007,17(2):185-212.和Verstraete.T,Alsalihi Z,Van den Braembussche RA.Multidisciplinaryoptimization of a radial compressor for micro gas turbineapplication[C].Proceedings of ASME Turbo Expo.2007,GT-2007-27484.)表明,叶轮变形与叶片形状和叶轮轮盘截面形状有密切的关系。一旦气动设计确定了叶片形状,那么叶轮变形和应力分布情况主要由轮盘截面形状决定。随着有限元仿真技术的发展,整体离心压气机叶轮设计完成后即可进行三维实体模型的建立以及强度应力的分析。通过对整体叶轮应力分布的分析来修改叶轮截面形状,降低叶轮的应力达到强度设计要求。以往的离心叶轮轮盘截面形状设计没有固定准则,对设计者的设计经验要求较高。而且叶轮强度设计往往需要较长的设计周期,设计效率低。
发明内容
本发明的目的在于提供一种将优化设计方法引入叶轮轮盘截面的设计过程中,能够大幅减少设计时间,提高设计结果的可靠性的离心压气机叶轮轮盘截面形状优化方法。该方法基于对离心压气机叶轮轮盘截面形状的优化,实现在不增加叶轮重量的条件下降低整体叶轮在工作转速下的等效应力大小,提高离心叶轮工作安全可靠性能。
为达到上述目的,本发明的方法包括以下步骤:
1)首先根据离心压气机叶轮轮盘截面模型数据以及压气机叶片数据对离心叶轮进行实体建模;
2)然后以叶轮出口背盘的厚度值、轮盘截面凹点深度和径向尺寸、圆台距轮盘右端面距离大小、轮盘加强肋条的径向尺寸和圆台半径等特征参数对轮盘进行参数化设计,将叶轮截面的形状尺寸用对应的设计参数表示,变化一个参数值,即改变与它相关的尺寸,通过调整参数来修改和控制几何形状,实现叶轮论盘截面的造型;
3)根据均匀设计方法对轮盘设计参数在一定范围内(10%-15%)变化并根据实验设计方法进行样本空间分布设计,从而设计出多个不同轮盘截面参数条件下的离心叶轮。然后进行实体建模,并计算叶轮体积V;
4)对多个不同截面形状下的离心叶轮在工作转速条件下进行有限元分析,并提取叶轮的最大等效应力和叶轮质量数据;
5)将不同的轮盘截面设计参数和对应的叶轮最大等效应力作为样本来训练神经网络近似模型,在参数变化范围内任意生成若干个检验神经网络精度的样本,对检验样本的等效应力采用有限元分析和神经网络预测并将结果比较,对比结果误差低于1%则训练的神经网络满足精度要求。若高于1%则将检验样本加入训练样本继续训练,直到满足精度要求。将训练好的神经网络近似模型作为优化的评价函数;
6)最后采用优化算法结合神经网络近似模型以离心叶轮最大等效应力为目标函数,以叶轮重量为算法约束条件,优化设计出离心叶轮等效应力最小的叶轮轮盘截面形状。
对离心叶轮轮盘截面的参数化设计时,轮盘背面选取了叶轮出口背盘的厚度值T、轮盘截面凹点深度D和径向尺寸R2、圆台距轮盘右端面距离大小Z、轮盘加强肋条的径向尺寸R1和圆台半径R3等特征参数,通过改变这些特征参数的大小,能够设计出不同形状的离心叶轮轮盘。
根据均匀设计方法对轮盘设计参数在一定范围(10%-15%)变化并进行样本空间分布设计,从而设计出多个不同轮盘截面参数条件下的离心叶轮;采用试验设计(DOE)方法在设计空间安排训练样本点,使训练样本在设计空间中均匀分布,在保证拟合精度的基础上可以最大限度的减少样本数,节约样本计算时间。均匀试验设计的试验点代表性较正交试验设计的试验点强得多,且能大幅度地减少试验次数。因而能够大大减少神经网络所需的训练样本数。均匀试验设计所使用的基本工具为均匀设计表,用Un(mk)表示,表中U是均匀设计表代号,n表示需要安排的试验次数,m表示每个因素的水平数,k表示该表最多安排的因素数。
对这些采用均匀设计后形成的不同截面形状下的离心叶轮在工作转速条件下进行有限元分析,通过三阶应力矩阵方程可以计算出对应的主应力σ0的分量(σ1,σ2,σ3)。
σ0=(σ1,σ2,σ3),
再根据第四强度理论计算各节点的等效应力为:
或者
其中σx、σy、σz为x,y,z方向上的正应力,σxy,σzy,σxz为xy,yz,zx平面上的切应力。
叶轮的质量为M=ρ*V,V为叶轮体积,ρ为材料密度。
在完成所有样本的有限元分析过程后,进行神经网络近似模型的建立。训练神经网络的输入值X为均匀设计表对应的样本叶轮设计参数,输出值D为样本叶轮的最大等效应力S和叶轮质量M。网络的输入矩阵X与输出矩阵D关系表达式为:
其中x为N维输入向量,μi为第i个隐节点的中心值,||||为欧式泛数,σi为基宽度,m为感受器单元的个数。
通过多个样本的训练过程,可以求出预测误差低于1%的对应的网络权值矩阵w,从而完成对神经网络近似模型的训练,可作为优化算法的评价函数来使用。
最后对叶轮截面几何结构的优化来降低叶轮应力,保证叶轮的安全可靠性。整体叶轮的最大应力Smax(Xi)为目标函数,算法约束条件为叶轮质量Mimpeller(Xi)不增加,优化问题目标是要降低应力,由下面的表达式描述:
优化过程完成后,等效应力最小的叶轮轮盘截面即为最优轮盘截面形状。
本发明离心压气机叶轮轮盘截面形状优化设计方法至少具有以下优点:本发明对离心叶轮轮盘截面形状采用了参数化建模的方式,以最少的设计参数描述叶轮轮盘截面形状,并将神经网络近似模型和优化设计方法结合起来对叶轮轮盘截面形状进行优化设计,在保证叶轮重量不增加的情况下有效降低叶轮的等效应力达7.4%以上,提高其工作可靠性。本设计方法大幅降低叶轮强度设计周期,设计结果可靠,节约了设计成本,具有很高的社会效益及推广价值。
附图说明
图1为某离心压气机叶轮实体模型图;
图2为某离心压气机叶轮轮盘截面参数化建模图;
图3为三层RBF径向基神经网络结构图;
图4为离心压气机叶轮轮盘截面优化前后应力分布对比图。
表1为某离心压气机轮盘截面参数均匀设计样本空间分布及对应整体叶轮应力和质量结果。
表2为神经网络训练后对随机样本的预测与有限元建模与分析结果比较。
表3为对离心压气机轮盘截面进行优化前后叶轮各参数比较。
具体实施方式
下面以某高速离心压气机叶轮的轮盘截面设计为例,结合附图,对离心压气机叶轮轮盘截面形状优化设计方法做详细说明:
1.叶轮实体模型建立
首先根据需要加工的离心压气机叶轮轮盘截面模型数据以及离心压气机叶轮叶片数据对离心压气机叶轮进行实体建模(图1);
2.叶轮实体模型的参数化
如图2,提取反映叶轮轮盘截面形状特征的关键数据为参数对轮盘进行参数化设计,对离心叶轮轮盘截面的参数化设计时,轮毂型线采用三次Bezier曲线进行拟和,轮毂型线由叶片气动设计过程确定,在强度设计过程中保持不变。轮盘背面选取了叶轮出口背盘的厚度值T、轮盘截面凹点深度D和径向尺寸R2、圆台距轮盘右端面距离大小Z、轮盘加强肋条的径向尺寸R1和圆台半径R3等特征参数,通过改变这些特征参数的大小,设计出不同形状的离心叶轮轮盘。
3.样本空间分布设计
根据均匀设计方法对轮盘设计参数在一定范围10%-15%变化,采用试验设计(DOE)方法在设计空间安排训练样本点,试验设计所使用的基本工具为均匀设计表,用Un(mk)表示,表中U是均匀设计表代号,n表示需要安排的试验次数,m表示每个因素的水平数,k表示该表最多安排的因素数。从而设计出多个不同轮盘截面参数条件下的离心叶轮,表1列出了6个设计参数在其变法范围内根据均匀设计表U21(216)均匀设计方法得到的所有组合;
4.对所有设计叶轮进行有限元分析
对这些不同截面形状下的离心叶轮在工作转速条件下进行有限元分析,并提取叶轮的最大等效应力数据和叶轮质量,其过程如下:
对这些不同截面形状下的离心压气机叶轮在工作转速条件下进行有限元分析,通过三阶应力矩阵方程计算出对应的主应力σ0的分量(σ1,σ2,σ3);
σ0=(σ1,σ2,σ3)
再根据第四强度理论计算各节点的最大等效应力为:
或者
其中σx、σy、σz为x,y,z方向上的正应力,σxy,σzy,σxz为xy,yz,zx平面上的切应力。
叶轮的质量为M=ρ*V,V为叶轮体积,ρ为材料密度。
结果如表1。
5.建立神经网络近似模型,将训练精度达到要求的近似模型作为优化算法的评价函数
将不同的轮盘截面设计参数Xi(R1,R2,R3,T,Z,D),i=1,...21和对应的叶轮最大等效应力Si,i=1,...21和质量Mi,i=1,...21作为样本来训练神经网络近似模型。图3为含一个输入、输出和隐含层的三层神经网络结构。表2为最后训练完成的神经网络对检验样本的预测结果与有限元分析结果的比较,神经网络对整体叶轮应力预测的最大相对误差为0.95%,可见RBF网络能够较好的对叶轮的强度性能进行预测。这样训练精度达到要求的神经网络就可以作为优化算法的评价函数来使用。
6.以离心叶轮最大应力为目标函数对叶轮截面进行优化设计
通过对叶轮截面几何结构的优化来降低叶轮应力,保证叶轮的安全可靠性。整体叶轮的最大应力Smax(Xi)为目标函数,算法约束条件为叶轮质量Mimpeller(Xi)不增加,优化问题目标是要降低应力,优化问题由下面的表达式描述:
优化过程完成后,等效应力最小的叶轮轮盘截面即为最优轮盘截面形状。
图4给出了目标函数收敛过程。可以看到,目标函数在算法迭代至80代以后基本收敛。优化后整个叶轮的最大应力从284MPa下降到263MPa,减小约7.4%叶轮的质量略有减小(图4,表3)。
表一
表二
表三
以上所述仅为本发明的一种实施方式,不是全部或唯一的实施方式,本领域普通技术人员通过阅读本发明书而对本发明技术方案采取的任何等效的交换,均为本发明的权利要求所涵盖。
Claims (3)
1.离心压气机叶轮轮盘截面形状优化方法,其特征在于:包括以下步骤:
1)首先根据需要加工的离心压气机叶轮轮盘截面模型数据以及离心压气机叶轮叶片数据对离心压气机叶轮进行实体建模;
2)然后以离心压气机叶轮出口背盘的厚度值、轮盘截面凹点深度和径向尺寸、圆台距轮盘右端面距离大小、轮盘加强肋条的径向尺寸和圆台半径特征参数对离心压气机叶轮轮盘截面进行参数化设计,将叶轮截面的形状尺寸用对应的设计参数表示,变化一个参数值,即改变与它相关的尺寸,通过调整参数来修改和控制几何形状,实现叶轮轮盘截面的造型;
3)根据均匀设计方法对离心压气机叶轮轮盘设计参数在10%-15%范围内变化进行样本空间分布设计,设计出多个不同轮盘截面参数条件下的离心压气机叶轮,然后进行实体建模,并计算离心压气机叶轮叶轮体积V;
4)对多个不同截面形状下的离心压气机叶轮在工作转速条件下进行有限元分析,并提取各个离心压气机叶轮的最大等效应力和离心压气机叶轮质量数据;
所述的对多个不同截面形状下的离心压气机叶轮在工作转速条件下进行有限元分析,通过三阶应力矩阵方程计算出对应的主应力σ0的分量(σ1,σ2,σ3);
σ0=(σ1,σ2,σ3)
再根据第四强度理论计算各节点的最大等效应力为:
或者
其中σx、σy、σz为x,y,z方向上的正应力,σxy,σzy,σxz为xy,yz,zx平面上的切应力;
叶轮的质量为M=p*V,V为叶轮体积,p为材料密度;
5)将不同的离心压气机叶轮轮盘截面设计参数和对应的离心压气机叶轮叶轮最大等效应力作为训练样本来训练神经网络近似模型,在参数变化范围任意生成若干个检验神经网络精度的样本,对检验样本的等效应力采用有限元分析和神经网络预测并将结果比较,对比结果误差低于1%则训练的神经网络满足精度要求,若高于1%则将检验样本加入训练样本继续训练,直到满足精度要求,将训练好的神经网络近似模型作为优化的评价函数;
所述的训练神经网络的输入值X为均匀设计表对应的样本叶轮设计参数,输出值D为样本叶轮的最大等效应力S和叶轮重量M,网络的输入矩阵X与输出矩阵D关系表达式为:
其中x为N维输入向量,μi为第i个隐节点的中心值,|| ||为欧式泛数,σi为基宽度,m为感受器单元的个数;
通过多个样本的训练过程求出预测误差低于1%的对应的网络权值矩阵w,训练好神经网络近似模型作为优化算法的评价函数;
6)最后采用优化算法结合神经网络近似模型以离心压气机叶轮最大等效应力为目标函数,以离心压气机叶轮重量为算法约束条件,得到离心压气机叶轮等效应力最小的离心压气机叶轮轮盘截面形状。
2.如权利要求1所述的离心压气机叶轮轮盘截面形状优化方法,其特征在于:所述步骤3)设计参数在10%-15%范围内变化,采用试验设计(DOE)方法在设计空间安排训练样本点,试验设计方法所使用的基本工具为均匀设计表,用Un(mk)表示,表中U是均匀设计表代号,n表示需要安排的试验次数,m表示每个因素的水平数,k表示该表最多安排的因素数。
3.如权利要求1所述的离心压气机叶轮轮盘截面形状优化方法,其特征在于:所述步骤6)以整体叶轮的最大应力Smax(Xi)为目标函数,算法约束条件为叶轮质量Mimpeller(Xi),优化问题目标是要降低等效应力,由下面的表达式描述:
优化过程完成后,等效应力最小的叶轮轮盘截面即为最优轮盘截面形状。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110324654 CN102339351B (zh) | 2011-10-24 | 2011-10-24 | 离心压气机叶轮轮盘截面形状优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110324654 CN102339351B (zh) | 2011-10-24 | 2011-10-24 | 离心压气机叶轮轮盘截面形状优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102339351A CN102339351A (zh) | 2012-02-01 |
CN102339351B true CN102339351B (zh) | 2013-07-10 |
Family
ID=45515079
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110324654 Active CN102339351B (zh) | 2011-10-24 | 2011-10-24 | 离心压气机叶轮轮盘截面形状优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102339351B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111720356A (zh) * | 2020-06-28 | 2020-09-29 | 辽宁福鞍燃气轮机有限公司 | 一种半开式离心压气机轴向推力平衡方法 |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR3009104B1 (fr) * | 2013-07-29 | 2015-08-21 | Snecma | Procede de modelisation d'une pale d'une helice non-carenee |
FR3023635B1 (fr) * | 2014-07-10 | 2018-05-25 | Safran Aircraft Engines | Procede de modelisation d'une baignoire d'une aube |
CN104331553B (zh) * | 2014-10-29 | 2017-07-11 | 浙江大学 | 考虑缺陷的大型透平膨胀机叶轮叶片结构优化设计方法 |
CN104573178B (zh) * | 2014-12-02 | 2017-05-31 | 中国航空动力机械研究所 | 一种整体叶轮强度有限元计算方法 |
CN106289893B (zh) * | 2015-05-25 | 2019-01-25 | 西安航空动力股份有限公司 | 一种轮盘盘心孔疲劳试验件的制造方法 |
CN106021773B (zh) * | 2016-05-30 | 2020-06-16 | 宁波至高点工业设计有限公司 | 一种针对叶扇自动适配算法的优化方法 |
CN106777482A (zh) * | 2016-11-18 | 2017-05-31 | 西北工业大学 | 一种基于网格参数化的结构多学科设计优化方法 |
CN106650125B (zh) * | 2016-12-28 | 2019-12-06 | 湘潭大学 | 一种离心压气机叶轮的优化方法和系统 |
CN107679270A (zh) * | 2017-08-28 | 2018-02-09 | 西北工业大学 | 离心压气机优化设计方法及系统 |
CN112214846A (zh) * | 2020-09-09 | 2021-01-12 | 浙江意动科技股份有限公司 | 一种减少应力集中的方法 |
CN113536474A (zh) * | 2021-06-04 | 2021-10-22 | 四川大学 | 一种离心式风机叶轮参数化设计系统及方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101363450A (zh) * | 2008-09-02 | 2009-02-11 | 西安交通大学 | 一种多翼式离心风机的叶轮结构 |
EP2098704A2 (en) * | 2008-03-05 | 2009-09-09 | United Technologies Corporation | Variable area fan nozzle fan flutter management system |
-
2011
- 2011-10-24 CN CN 201110324654 patent/CN102339351B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2098704A2 (en) * | 2008-03-05 | 2009-09-09 | United Technologies Corporation | Variable area fan nozzle fan flutter management system |
CN101363450A (zh) * | 2008-09-02 | 2009-02-11 | 西安交通大学 | 一种多翼式离心风机的叶轮结构 |
Non-Patent Citations (3)
Title |
---|
席光,王志恒,王尚锦.《叶轮机械气动优化设计中的近似模型方法及其应用》.《西安交通大学学报》.2007,第41卷(第2期),第125-135、184页. * |
王宏亮,席光.《离心压气机叶轮多学科优化设计方法研究》.《工程热物理学报》.2010,第31卷(第11期),第1835-1838页. * |
王宏亮,席光.《离心叶轮几何形变对气动性能的影响》.《西安交通大学学报》.2009,第43卷(第5期),第46-50页. * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111720356A (zh) * | 2020-06-28 | 2020-09-29 | 辽宁福鞍燃气轮机有限公司 | 一种半开式离心压气机轴向推力平衡方法 |
Also Published As
Publication number | Publication date |
---|---|
CN102339351A (zh) | 2012-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102339351B (zh) | 离心压气机叶轮轮盘截面形状优化方法 | |
Dhert et al. | Aerodynamic shape optimization of wind turbine blades using a Reynolds‐averaged Navier–Stokes model and an adjoint method | |
Li et al. | Deep learning regression-based stratified probabilistic combined cycle fatigue damage evaluation for turbine bladed disks | |
Guida et al. | SPH–Lagrangian study of bird impact on leading edge wing | |
Mengistu et al. | Aerodynamic optimization of turbomachinery blades using evolutionary methods and ANN-based surrogate models | |
Toal et al. | Efficient multipoint aerodynamic design optimization via cokriging | |
Fei et al. | Vectorial surrogate modeling method for multi-objective reliability design | |
Javed et al. | Optimization of a centrifugal compressor impeller for robustness to manufacturing uncertainties | |
CN102789539B (zh) | 一种飞行器升力面结构优化设计方法 | |
Karsh et al. | A hybrid stochastic sensitivity analysis for low-frequency vibration and low-velocity impact of functionally graded plates | |
CN104866652A (zh) | 一种基于abaqus的喷丸强化变形的有限元模拟方法 | |
Han et al. | Surrogate-based aerodynamic shape optimization with application to wind turbine airfoils | |
CN106934074A (zh) | 全局最优涡扇发动机进气道减噪设计方法 | |
Fei et al. | Distributed collaborative extremum response surface method for mechanical dynamic assembly reliability analysis | |
Li et al. | Designable mechanical properties of modified body-centered cubic lattice materials | |
Vesting | Marine propeller optimisation-strategy and algorithm development | |
Caprace et al. | Incorporating High-Fidelity Aerostructural Analyses in Wind Turbine Rotor Optimization | |
Xu et al. | A study of Nash-EGO algorithm for aerodynamic shape design optimizations | |
Díaz-Casás et al. | Wind turbine design through evolutionary algorithms based on surrogate CFD methods | |
CN103870614A (zh) | 一种结构概率优化设计方法 | |
Esfahlani et al. | Hexagonal honeycomb cell optimisation by way of meta-model techniques | |
Garcia et al. | CFD based wing shape optimization through gradient-based method | |
Yang et al. | Surrogate Model-Based Acoustic Optimisation of Jet Nozzle Exit Geometry | |
Changli et al. | Surrogate model-based optimization for the headform design of an axisymmetric body | |
Yadav et al. | Optimization of car rim using Optistruct |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |