CN103310046A - 基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法 - Google Patents
基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法 Download PDFInfo
- Publication number
- CN103310046A CN103310046A CN2013102109470A CN201310210947A CN103310046A CN 103310046 A CN103310046 A CN 103310046A CN 2013102109470 A CN2013102109470 A CN 2013102109470A CN 201310210947 A CN201310210947 A CN 201310210947A CN 103310046 A CN103310046 A CN 103310046A
- Authority
- CN
- China
- Prior art keywords
- optimization
- hyperelliptic
- design
- eta
- hyperelliptic curve
- 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
Links
Images
Landscapes
- Control Of Turbines (AREA)
- Turbine Rotor Nozzle Sealing (AREA)
Abstract
本发明公开了一种基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法,用于解决现有多弧段曲线形状优化设计方法设计周期长的技术问题。技术方案是采用超椭圆曲线边界条件对带辐板涡轮盘截面进行有限元建模、分析和灵敏度求解,由求得的灵敏度信息进行优化迭代,直至收敛到最终的设计结果。该方法对涡轮盘内腔的初始形状轮廓进行了有效的优化设计,设计周期短,并获得了较好的应力水平分布。优化后的最大应力为769.47MPa,比优化之前的1460.75MPa降低了47.32%。文献所述方法优化后的最大应力为842.7MPa,比优化之前的912MPa仅降低了7.6%。故本发明方法优化结果的应力集中水平低,明显优于背景技术。
Description
技术领域
本发明涉及一种带辐板涡轮盘内腔形状优化设计方法,特别涉及一种基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法。
背景技术
航空发动机涡轮盘为航空发动机中的关键零件之一,其体积和重量较大,在工作中受载荷包括涡轮盘高速旋转产生的自身离心力和叶片离心力,以及来自燃烧室高温、高压燃气的温度梯度载荷和叶片气动载荷等。因此涡轮盘在工作中所承受的载荷比较复杂,容易产生应力集中现象,降低其疲劳寿命。
文献“朱继宏,李军朔等,现代形状优化技术在航空发动机零部件设计中的应用。航空制造技术,2012(23/24)”公开了一种多圆弧曲线形状优化设计方法,实现了带辐板涡轮盘内腔的形状优化。文献公开的方法第一步使用常规的自由曲线进行优化迭代,这一步的结果对最终的形状优化结果影响并不大,因此这一步的迭代大大浪费了计算机机时,增长了优化设计周期。第二步采用多弧段曲线对自由曲线优化结果进行逼近,第三步采用多弧段的边界条件,重新定义圆弧设计变量,然后再次进行优化迭代,增加了整个优化过程中所需定义的设计变量个数,增加了优化过程的复杂性。优化后得到的多圆弧曲线,使用数学方程描述比较困难,需要知道每段圆弧曲线的起点、端点以及半径才能将其完整地描述。
文献所述方法的自由曲线初始设计的应力最大值为912MPa,多弧段曲线优化设计后的最大应力为842.7MPa,降幅为7.6%。
发明内容
为了克服现有多弧段曲线形状优化设计方法设计周期长的不足,本发明提供一种基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法。该方法采用超椭圆曲线边界条件对带辐板涡轮盘截面进行有限元建模、分析和灵敏度求解,由求得的灵敏度信息进行优化迭代,直至收敛到最终的设计结果。这种方法能够对涡轮盘内腔的初始形状轮廓进行有效的优化设计,设计周期短,并获得较好的应力水平分布设计,同时设计结果表述清晰。
本发明解决其技术问题所采用的技术方案是:一种基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法,其特点是包括以下步骤:
步骤一、定义超椭圆曲线的控制参数,构造超椭圆曲线。超椭圆曲线的解析式方程为:
其中,a、b和η均为正数,a和b中数值大的为超椭圆的长半轴,数值小的为超椭圆的短半轴,η为指数。与式(1)对应的参数方程为:
当a、b和η取不同的数值时,超椭圆的曲线具有不同的形状。
步骤二、采用部分超椭圆曲线定义优化模型的边界,选择开口宽度width作为另一个设计变量。取开口处两个点S、R的y坐标分别为yS、yR:
计算出开口处R点与超椭圆中心O的连线和长轴的夹角γ:
计算出超椭圆的中心坐标(x0,y0):
得到优化过程中所用的部分超椭圆曲线描述为:
步骤三、对于带辐板涡轮盘来说,还要选择开口外侧两处的倒圆角半径R1和R2作为设计变量。
步骤四、以部分超椭圆曲线为形状边界,构造形状优化的有限元模型:
find x二(x1,x2,…,xn)
min Φ(X) (7)
s.t.KU=F
其中,X为设计域上的形状变量向量;n为设计变量个数;Φ(X)为形状优化的目标函数;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点位移向量;Gj(X)为第j个约束函数;为第j个约束函数的上限;J为约束的数量。
步骤五、定义优化模型的设计变量,并对模型施加约束和边界载荷。
步骤六、用有限元软件Ansys将模型进行一次有限元分析;再通过结构优化平台Boss-Quattro进行优化灵敏度分析,求得目标函数和约束条件的灵敏度,选取梯度优化算法GCMMA进行优化设计,得到优化结果。
本发明的有益效果是:由于该方法采用超椭圆曲线边界条件对带辐板涡轮盘截面进行有限元建模、分析和灵敏度求解,由求得的灵敏度信息进行优化迭代,直至收敛到最终的设计结果。这种方法能够对涡轮盘内腔的初始形状轮廓进行有效的优化设计,设计周期短,并获得较好的应力水平分布设计,同时设计结果表述清晰。优化后的最大应力为769.47MPa,比优化之前的1460.75MPa降低了47.32%。
以下结合附图和实施例详细说明本发明。
附图说明
图1是超椭圆曲线在选取不同参数时的不同形状的示意图。
图2是带辐板涡轮盘截面模型尺寸形状以及受力和边界约束示意图
图3是带辐板涡轮盘内腔截面轮廓使用部分多圆弧曲线形状边界条件示意图。
图4是具体实施方式中内孔的多弧段设计结果图。
图5是具体实施例应用本发明方法的带辐板涡轮盘内腔形状优化设计结果。
具体实施方式
参照图1~5。以某发动机涡轮盘横截面的形状优化为例来说明本发明。涡轮绕X轴旋转。为简化模型,将涡轮叶片的离心载荷等效加载在涡轮盘径向外边界FG上。给定初始条件为涡轮叶片等效离心载荷为60MPa,涡轮转速为1400rad/s,涡轮盘温度为TAP=400℃,TFG=600℃,其间温度为线性分布。其整体结构虽不具有对称性,但考虑加工制造要求,要求涡轮盘内腔截面曲线关于O1O2对称,对称轴O1O2垂直于X轴。其中S、T、A、D、C、E、F分别与R、Q、P、I、J、H、G关于O1O2对称,T、A、D处的倒圆角半径值分别与Q、P、I处的倒圆角半径值相等,RI=5mm。该涡轮盘的材料为钛合金,其杨氏模量为1.6×105MPa,泊松比为0.3,密度为8.24×10-9ton/mm3,热膨胀系数为1.1×10-5。限制最左侧AB段和最右侧MN段辐板处X方向位移以及最下端AT段和PQ段的Y方向位移。设计涡轮盘内腔截面的形状,使得整个涡轮盘的质量不大于78kg。
步骤一、使用APDL语言自底向上建立涡轮盘截面的有限元模型:设定网格边长为2mm,自由划分网格。定义载荷:在涡轮盘截面Y方向最大处FG段节点上施加-60MPa(即方向为Y轴正方向的拉力)的压力为涡轮叶片的等效离心载荷,对涡轮整体施加1400rad/s的离心载荷,再对涡轮整体从AP处沿Y方向到FG处施加温度载荷,设定AP处温度为400℃,FG处温度为600℃,其间温度为线性分布。定义边界条件:限制最左侧AB段上的节点和最右侧MN段上的节点辐板处X方向位移以及最下端AT段上的节点和PQ段上的节点的Y方向位移。定义超椭圆曲线的控制参数,构造超椭圆曲线。超椭圆曲线的解析式方程为:
其中,a、b和η均为正数,a和b中数值大的为超椭圆的长半轴,数值小的为超椭圆的短半轴,η为指数。与式(1)对应的参数方程为:
当a、b和η取不同的数值时,超椭圆的曲线具有不同的形状,因此超椭圆曲线的优化设计方法实用性很强。
步骤二、采用部分(一般大于一半)超椭圆曲线定义优化模型的边界,因此还要选择开口宽度width作为另一个设计变量。取开口处两个点S、R的y坐标分别为yS、yR,它们都是定值:
计算出开口处R点与超椭圆中心O的连线和长轴的夹角γ:
计算出超椭圆的中心坐标(x0,y0):
步骤三、对于带辐板涡轮盘来说,还要选择开口外侧两处的倒圆角半径R1和R2作为设计变量。
步骤四、以部分超椭圆曲线为形状边界,构造形状优化的有限元模型:
find X=(x1,x2,...,xn)
min Φ(X) (7)
s.t.KU=F
其中,X为设计域上的形状变量向量;n为设计变量个数;Φ(X)为形状优化的目标函数;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点位移向量;Gj(X)为第j个约束函数;为第j个约束函数的上限;J为约束的数量。
步骤五、定义优化模型的设计变量,并对模型施加约束和边界载荷。
步骤六、用有限元软件Ansys将模型进行一次有限元分析;再通过结构优化平台Boss-Quattro进行优化灵敏度分析,求得目标函数和约束条件的灵敏度,选取梯度优化算法GCMMA(Globally Convergent Method of Moving Asymptotes)优化算法进行优化设计,得到优化结果。
由优化结果可以看出,应用本发明方法进行带辐板涡轮盘内腔的形状优化设计,得到了一个部分超椭圆曲线的设计结果。与文献中的方法相比,本发明方法所使用的设计变量减少了3个,而且不需要多次定义设计变量,优化过程简单,减少了形状优化的迭代时间。使用本发明方法优化后的最大应力为769.47MPa,比优化之前的1460.75MPa降低了47.32%。而使用文献中的方法,其自由曲线初始设计的应力最大值为912MPa,多弧段曲线优化设计后的最大应力为842.7MPa,降幅仅为7.6%。明显本发明所采用的方法优化结果的应力集中水平低,设计变量少,易于描述。
Claims (1)
1.一种基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法,其特征在于包括以下步骤:
步骤一、定义超椭圆曲线的控制参数,构造超椭圆曲线;超椭圆曲线的解析式方程为:
其中,a、b和η均为正数,a和b中数值大的为超椭圆的长半轴,数值小的为超椭圆的短半轴,η为指数;与式(1)对应的参数方程为:
当a、b和η取不同的数值时,超椭圆的曲线具有不同的形状;
步骤二、采用部分超椭圆曲线定义优化模型的边界,选择开口宽度width作为另一个设计变量;取开口处两个点S、R的y坐标分别为yS、yR:
计算出开口处R点与超椭圆中心O的连线和长轴的夹角γ:
计算出超椭圆的中心坐标(x0,y0):
得到优化过程中所用的部分超椭圆曲线描述为:
步骤三、对于带辐板涡轮盘来说,还要选择开口外侧两处的倒圆角半径R1和R2作为设计变量;
步骤四、以部分超椭圆曲线为形状边界,构造形状优化的有限元模型:
find X=(x1,x2,...,xn)
min Φ(X)
s.t.KU=F (7)
其中,X为设计域上的形状变量向量;n为设计变量个数;Φ(X)为形状优化的目标函数;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点位移向量;Gj(X)为第j个约束函数;为第j个约束函数的上限;J为约束的数量;
步骤五、定义优化模型的设计变量,并对模型施加约束和边界载荷;
步骤六、用有限元软件Ansys将模型进行一次有限元分析;再通过结构优化平台Boss-Quattro进行优化灵敏度分析,求得目标函数和约束条件的灵敏度,选取梯度优化算法GCMMA进行优化设计,得到优化结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310210947.0A CN103310046B (zh) | 2013-05-30 | 2013-05-30 | 基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310210947.0A CN103310046B (zh) | 2013-05-30 | 2013-05-30 | 基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103310046A true CN103310046A (zh) | 2013-09-18 |
CN103310046B CN103310046B (zh) | 2016-01-20 |
Family
ID=49135258
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310210947.0A Active CN103310046B (zh) | 2013-05-30 | 2013-05-30 | 基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103310046B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103790709A (zh) * | 2014-02-19 | 2014-05-14 | 中国航空动力机械研究所 | 涡轮轮盘 |
CN104268317A (zh) * | 2014-09-12 | 2015-01-07 | 西北工业大学 | 机械零件圆角结构形状优化方法 |
CN106089310A (zh) * | 2016-07-01 | 2016-11-09 | 中航空天发动机研究院有限公司 | 一种提高辐板承载能力的纤维增强双辐板涡轮盘 |
CN106202686A (zh) * | 2016-07-06 | 2016-12-07 | 中南大学 | 一种涡轮盘等温模锻预成形坯料的多目标设计方法 |
CN108664704A (zh) * | 2018-04-12 | 2018-10-16 | 南京航空航天大学 | 一种五参量多圆弧异型孔均衡优化方法 |
CN109033635A (zh) * | 2018-07-26 | 2018-12-18 | 中车青岛四方车辆研究所有限公司 | 一种s型辐板车轮优化设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101840452A (zh) * | 2010-04-22 | 2010-09-22 | 西北工业大学 | 带孔薄壁曲面结构的空间孔洞优化设计方法 |
CN102254066A (zh) * | 2011-06-30 | 2011-11-23 | 西北工业大学 | 带孔薄壁曲壳结构的曲面形状与孔洞形状协同优化设计方法 |
CN102368282A (zh) * | 2011-11-17 | 2012-03-07 | 西北工业大学 | 带孔薄壁曲面结构的孔洞隐式表达优化设计方法 |
-
2013
- 2013-05-30 CN CN201310210947.0A patent/CN103310046B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101840452A (zh) * | 2010-04-22 | 2010-09-22 | 西北工业大学 | 带孔薄壁曲面结构的空间孔洞优化设计方法 |
CN102254066A (zh) * | 2011-06-30 | 2011-11-23 | 西北工业大学 | 带孔薄壁曲壳结构的曲面形状与孔洞形状协同优化设计方法 |
CN102368282A (zh) * | 2011-11-17 | 2012-03-07 | 西北工业大学 | 带孔薄壁曲面结构的孔洞隐式表达优化设计方法 |
Non-Patent Citations (10)
Title |
---|
PEDERSEN N L: "Opminization of holes in plates for control of eigenfrequencies", 《STRUCT MULTIDISC OPTIM》 * |
WEIHONG ZHANG等: "A parametric mapping method for curve shape optimization on 3D panel structures", 《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING》 * |
周鸣等: "基于超椭圆函数的压力容器封头形状优化设计", 《中国计算力学大会’2010(CCCM2010)暨第八届南方计算力学学术会议(SCCM8)论文集》 * |
张乘齐等: "基于BOSS平台的热力耦合场涡轮盘优化设计", 《燃气涡轮试验与研究》 * |
张卫红等: "压力载荷下的结构拓扑_形状协同优化", 《航空学报》 * |
张正印: "关于超圆、超球等几何图形面、体积的计算", 《数学通报》 * |
张苗苗等: "带孔圆柱壳轴压屈曲与孔形优化设计", 《科学技术与工程》 * |
朱继宏等: "现代形状优化技术在航空发动机零部件设计中的应用", 《航空制造技术》 * |
芮井中: "具有最小应力集中的缺口形状优化方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》 * |
龚曙光: "形状优化数值方法的研究与应用", 《中国优秀博硕士学位论文全文数据库 (硕士) 基础科学辑》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103790709A (zh) * | 2014-02-19 | 2014-05-14 | 中国航空动力机械研究所 | 涡轮轮盘 |
CN104268317A (zh) * | 2014-09-12 | 2015-01-07 | 西北工业大学 | 机械零件圆角结构形状优化方法 |
CN106089310A (zh) * | 2016-07-01 | 2016-11-09 | 中航空天发动机研究院有限公司 | 一种提高辐板承载能力的纤维增强双辐板涡轮盘 |
CN106089310B (zh) * | 2016-07-01 | 2017-08-29 | 中航空天发动机研究院有限公司 | 一种提高辐板承载能力的纤维增强双辐板涡轮盘 |
CN106202686A (zh) * | 2016-07-06 | 2016-12-07 | 中南大学 | 一种涡轮盘等温模锻预成形坯料的多目标设计方法 |
CN106202686B (zh) * | 2016-07-06 | 2019-06-14 | 中南大学 | 一种涡轮盘等温模锻预成形坯料的多目标设计方法 |
CN108664704A (zh) * | 2018-04-12 | 2018-10-16 | 南京航空航天大学 | 一种五参量多圆弧异型孔均衡优化方法 |
CN108664704B (zh) * | 2018-04-12 | 2021-05-07 | 南京航空航天大学 | 一种五参量多圆弧异型孔均衡优化方法 |
CN109033635A (zh) * | 2018-07-26 | 2018-12-18 | 中车青岛四方车辆研究所有限公司 | 一种s型辐板车轮优化设计方法 |
CN109033635B (zh) * | 2018-07-26 | 2022-10-21 | 中车青岛四方车辆研究所有限公司 | 一种s型辐板车轮优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103310046B (zh) | 2016-01-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103310046B (zh) | 基于超椭圆曲线的带辐板涡轮盘内腔形状优化设计方法 | |
CN102799730B (zh) | 一种燃气轮机风扇叶片反扭过程的预估方法 | |
Krack et al. | Reliability optimization of friction-damped systems using nonlinear modes | |
Jutte et al. | Aeroelastic tailoring of the NASA common research model via novel material and structural configurations | |
CN105930620A (zh) | 一种叶轮结构的热耦合振动特性有限元分析方法 | |
Murthy et al. | Nonparametric stochastic modeling of uncertainty in rotordynamics—Part II: Applications | |
Diener et al. | Multi-Disciplinary optimization of a mixed-flow compressor impeller | |
Resor et al. | Aeroelastic instability of very large wind turbine blades | |
CN112765732A (zh) | 一种基于选区激光熔化工艺的航空叶片拓扑优化设计方法 | |
Faskhutdinov et al. | Topology optimization of a gas-turbine engine part | |
Waite et al. | The impact of blade loading and unsteady pressure bifurcations on low-pressure turbine flutter boundaries | |
CN111079326B (zh) | 二维各向异性网格单元度量张量场光滑化方法 | |
Muhamad et al. | Design improvement of steering knuckle component using shape optimization | |
Cheng et al. | Topology optimization of turbine disk considering maximum stress prediction and constraints | |
CN111079327B (zh) | 一种基于各向异性过滤技术的异形曲面加筋拓扑优化方法 | |
CN109815518A (zh) | 基于转动惯量约束的飞行器舵面设计方法 | |
Song et al. | Blade shape optimization of transonic axial flow fan in terms of sectional profiles and stacking line | |
Zhang et al. | Mistuning effects on aero-elastic stability of contra-rotating turbine blades | |
CN107092726A (zh) | 一种提高飞轮转子储能密度的结构设计方法及装置 | |
Martel et al. | Maximum mistuned vibration amplitude with aerodynamic damping | |
CN106777538A (zh) | 一种基于有限细胞描述的承载结构拓扑优化设计方法 | |
CN113486512A (zh) | 一种功能梯度变厚度叶片模型的颤振分析方法 | |
CN106156443A (zh) | 一种在冷却塔上布设加劲环和子午肋的方法及装置 | |
CN103279594B (zh) | 一种二维固-固声子晶体z模态带隙优化方法 | |
Wang et al. | Intercoupling aeroelastic stability topology optimization for the design of bend-twisted hollow blade in a turbomachinery cascade |
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 |