CN106971023B - 基于超椭圆曲线的轮盘异型孔结构设计方法 - Google Patents

基于超椭圆曲线的轮盘异型孔结构设计方法 Download PDF

Info

Publication number
CN106971023B
CN106971023B CN201710117739.4A CN201710117739A CN106971023B CN 106971023 B CN106971023 B CN 106971023B CN 201710117739 A CN201710117739 A CN 201710117739A CN 106971023 B CN106971023 B CN 106971023B
Authority
CN
China
Prior art keywords
hole
special
shaped hole
hyperelliptic
wheel disc
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
Application number
CN201710117739.4A
Other languages
English (en)
Other versions
CN106971023A (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.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
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 Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710117739.4A priority Critical patent/CN106971023B/zh
Publication of CN106971023A publication Critical patent/CN106971023A/zh
Application granted granted Critical
Publication of CN106971023B publication Critical patent/CN106971023B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)
  • Stereophonic System (AREA)

Abstract

本发明提供了一种基于超椭圆曲线的轮盘异型孔结构设计方法,用于解决目前轮盘孔结构的孔边应力集中问题,提高了设计结果的实用性,简化了异型孔结构的设计过程。技术方案是采用超椭圆曲线对轮盘孔类结构进行改进;构建一种超椭圆异型孔多目标优化模型,使超椭圆异型孔轮廓同时满足两个优化目标:1)孔边应力水平降低到设置水平2)孔形状变化尽可能小以满足传力要求;该模型可通过优化获得最佳设计。优化得到的超椭圆孔可按设计者需求将孔边最大应力降低10%~20%,且异型孔轮廓形状与先前圆孔的差异度较小,可满足原有的装配,保证传力可靠。此外,该方法所需设计变量较少,异型孔轮廓的数学模型简洁。

Description

基于超椭圆曲线的轮盘异型孔结构设计方法
技术领域
本发明涉及航空发动机涡轮转子部件的孔类结构优化设计领域,具体是一种基于超椭圆曲线的轮盘异型孔结构设计方法。
背景技术
航空发动机涡轮转子部件工作环境恶劣,其轮盘幅板上的螺栓孔、通气孔等孔类结构常因应力集中现象,产生孔边裂纹,成为轮盘失效的一个重要原因。传统措施如采用孔边倒角、抛光、增加结构重量等方法并不能从根本上解决该问题。
发明内容
本发明为了解决现有技术的问题,提供了一种基于超椭圆曲线的轮盘异型孔结构设计方法,按结构使用需求降低孔边应力,同时保证孔轮廓形状变化尽可能小,以此满足强度要求、装配要求并保证可靠传力。该方法所需设计变量较少,异型孔数学模型简洁。
本发明包括以下步骤:
步骤一、定义超椭圆异型孔曲线方程。
Figure BDA0001236080140000011
式中,a、b为超椭圆曲线的半轴值,m、n为指数,其取值不同,超椭圆异型孔形状不同,此发明中取m,n>2,此时异型孔轮廓介于方孔和圆孔之间,更加符合轮盘孔类结构的实际设计要求。
步骤二、采用超椭圆参数方程运用有限元分析软件进行参数化建模:
Figure BDA0001236080140000012
式中,θ为参数,sgn为符号判断函数:
Figure BDA0001236080140000013
步骤三、构建超椭圆异型孔多目标优化模型:
minf{f1max),f2(m,n)}; (4)
Figure BDA0001236080140000021
其中,优化设计变量m、n为超椭圆曲线指数;σmax=σmax(m,n)为设计点λ(m,n)对应的超椭圆孔孔边最大主应力值,目标函数f1max)代表孔边应力降低程度;目标函数f2(m,n)代表超椭圆异型孔轮廓变化;upper和lower分别为曲线指数的上、下界;
目标函数f2(m,n)计算方法如方程如下所示:
f2(m,n)=S1-S0。 (6)
S1为超椭圆异型孔面积,S0为原圆孔面积。
步骤四、运用多目标优化方法对超椭圆异型孔优化模型进行寻优,最终得到异型孔轮廓优化结果。
本发明有益效果在于:
1、运用本发明得到的超椭圆异型孔结构,在不增加结构重量和影响其他结构的基础上,按设计者需求应力目标,使孔边应力降低10%~20%。
2、本发明运用多目标优化方法,得到的异型孔在满足强度要求的前提下,轮廓尺寸改变最小,并保证了传力可靠,设计稳健性较好。
3、超椭圆曲线异型孔数学模型结构简洁,设计变量较少。
附图说明
图1是孔边应力降低10%时,单参量(m=n)超椭圆异型孔轮廓优化结果;
图2是孔边应力降低10%时,两参量(m≠n)超椭圆异型孔轮廓优化结果;
其中:1-原螺栓通孔(基圆)。
具体实施方式
下面结合附图对本发明作进一步说明。
下面结合具体实施方式对本发明作更进一步的说明。
首先,定义超椭圆异型孔曲线方程。
本发明中采用如下超椭圆曲线方程对轮盘孔类结构进行改进:
Figure BDA0001236080140000031
其中a、b为超椭圆曲线的半轴值,m、n为指数,其取值不同,超椭圆异型孔形状不同,此发明中取m,n>2,此时异型孔轮廓介于方孔和圆孔之间,更加符合轮盘孔类结构的实际设计要求。
按方程(2)运用有限元分析软件对超椭圆参数化方程进行建模,进一步对孔边应力情况进行分析:
Figure BDA0001236080140000032
式中,θ为参数,单位为角度,sgn为符号判断函数:
Figure BDA0001236080140000033
构建超椭圆异型孔多目标优化模型:
minf{f1max),f2(m,n)}; (4)
Figure BDA0001236080140000034
其中,设计变量m、n为超椭圆曲线指数,它们代表了设计点λ(m,n);upper(=6~8)为m、n上界,lower(=2)为m、n下界;σmax=σmax(m,n)为设计点λ(m,n)对应的超椭圆孔边最大主应力值。目标函数f1max)代表孔边应力降低率,可由设计者根据需要设定(一般下降率可为10%~20%);目标函数f2(m,n)代表超椭圆异型孔轮廓变化程度(异型度),其计算方法如方程(16)所示。
f2(m,n)=S1-S0; (6)
S1为超椭圆异型孔面积,S0为原圆孔面积。
例:运用多目标优化方法对超椭圆异型孔优化模型进行寻优,设a=b取为原螺栓通孔基圆半径,对某涡轮盘安装边螺栓孔进行优化,得到超椭圆异型孔轮廓优化结果如下:
1、当超椭圆优化模型为单参数时,异型孔轮廓优化结果如附图1所示,可降低孔边应力10%~20%,轮廓面积变化为5%~18%,图中虚线部分为原螺栓通孔1。
2、当超椭圆优化模型为双参数时,异型孔轮廓优化结果如附图2所示,可降低孔边应力10%~20%,轮廓面积变化为5%~17%,图中虚线部分为原螺栓通孔1。
本发明提出了一种采用超椭圆曲线对轮盘孔类结构进行改进的异型孔设计方法,使优化得到的超椭圆异型孔边最大应力可按结构使用要求降低10%~20%,且异型孔轮廓形状与原圆孔的差异度尽可能小,可满足原有的强度、装配和传力要求。此外,该方法所需设计变量较少,异型孔轮廓的数学模型简洁。
本发明具体应用途径很多,以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以作出若干改进,这些改进也应视为本发明的保护范围。

Claims (2)

1.一种基于超椭圆曲线的轮盘异型孔结构设计方法,其特征在于包括以下步骤:
步骤一、定义超椭圆异型孔曲线方程:
Figure FDA0002392842930000011
式中,a、b为超椭圆曲线的半轴值,m、n为指数;
步骤二、采用超椭圆参数方程运用有限元分析软件进行参数化建模:
Figure FDA0002392842930000012
式中,θ为参数,sgn为符号判断函数:
Figure FDA0002392842930000013
步骤三、构建超椭圆异型孔多目标优化模型:
min f{f1max),f2(m,n)}; (4)
Figure FDA0002392842930000014
其中,优化设计变量m、n为超椭圆曲线指数;σmax=σmax(m,n)为设计点λ(m,n)对应的超椭圆孔孔边最大主应力值,目标函数f1max)代表孔边应力降低程度;目标函数f2(m,n)代表超椭圆异型孔轮廓变化;upper和lower分别为曲线指数的上、下界;
目标函数f2(m,n)计算方法如方程如下所示:
f2(m,n)=S1-S0; (6)
S1为超椭圆异型孔面积,S0为原圆孔面积;
步骤四、运用多目标优化方法对超椭圆异型孔优化模型进行寻优,最终得到异型孔轮廓优化结果。
2.根据权利要求1所述的一种基于超椭圆曲线的轮盘异型孔结构设计方法,其特征在于:步骤一中,m>2,n>2。
CN201710117739.4A 2017-03-01 2017-03-01 基于超椭圆曲线的轮盘异型孔结构设计方法 Active CN106971023B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710117739.4A CN106971023B (zh) 2017-03-01 2017-03-01 基于超椭圆曲线的轮盘异型孔结构设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710117739.4A CN106971023B (zh) 2017-03-01 2017-03-01 基于超椭圆曲线的轮盘异型孔结构设计方法

Publications (2)

Publication Number Publication Date
CN106971023A CN106971023A (zh) 2017-07-21
CN106971023B true CN106971023B (zh) 2020-05-19

Family

ID=59328378

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710117739.4A Active CN106971023B (zh) 2017-03-01 2017-03-01 基于超椭圆曲线的轮盘异型孔结构设计方法

Country Status (1)

Country Link
CN (1) CN106971023B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108664704B (zh) * 2018-04-12 2021-05-07 南京航空航天大学 一种五参量多圆弧异型孔均衡优化方法
CN111353249B (zh) * 2020-03-02 2022-02-11 厦门大学 涡轮封严盘非圆通气孔集成设计优化方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105046000A (zh) * 2015-07-20 2015-11-11 南京航空航天大学 基于工程离散变量的两变量多圆弧异型孔均衡设计方法及模型

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7440569B2 (en) * 2003-07-28 2008-10-21 Microsoft Corporation Tate pairing techniques for use with hyperelliptic curves

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105046000A (zh) * 2015-07-20 2015-11-11 南京航空航天大学 基于工程离散变量的两变量多圆弧异型孔均衡设计方法及模型

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Memory efficient hyperelliptic curve point counting;Hendrik Hubrechts;《International Journal of Number Theory》;20111130;第07卷(第01期);第1-13页 *
基于粒子群神经网络的轮盘优化;殷艺云 等;《航空动力学报》;20070930;第22卷(第9期);第1578-1582页 *

Also Published As

Publication number Publication date
CN106971023A (zh) 2017-07-21

Similar Documents

Publication Publication Date Title
CN109460566B (zh) 一种风力机叶片内侧厚翼型的气动稳健优化设计方法
CN107391891B (zh) 一种基于模型融合方法的大展弦比机翼优化设计方法
US20170306970A1 (en) Extrapolation method of low rotational speed characteristic of compressor
CN106650125B (zh) 一种离心压气机叶轮的优化方法和系统
JP6708995B2 (ja) 軸流流体機械の翼の設計方法及び翼
CN110727995A (zh) 叶片形状的构造方法、叶片以及计算机设备
CN106971023B (zh) 基于超椭圆曲线的轮盘异型孔结构设计方法
CN105678015B (zh) 一种高超声速三维机翼的非概率可靠性气动结构耦合优化设计方法
CN106050739A (zh) 一种冷却风扇用高性能翼型
CN109408941B (zh) 基于数据挖掘和遗传算法的飞行器气动优化方法
CN109359360B (zh) 一种基于局部特征的结构应力处理方法
Yanhui et al. Performance improvement of optimization solutions by POD-based data mining
CN113569360A (zh) 一种风力机叶片抗颤振翼型簇设计方法
Buyukli et al. High-loaded compressor blisk-type impeller multidisciplinary optimization
CN106294908B (zh) 声衬设计方法
CN109409020B (zh) 一种空心风扇叶片的建模方法
CN112329196A (zh) 翼型的结构几何参数确定方法、装置、设备及存储介质
CN104376162B (zh) 叶片建模方法
CN106156409B (zh) 一种基于最小轮廓变化的两参量异型螺栓孔稳健设计方法
KR101811135B1 (ko) 최적 정밀도 결정을 통한 모델 단순화 장치 및 방법
CN111859817B (zh) 一种基于cfd软件激波模拟能力的飞行器气动选型方法
CN109241585B (zh) 一种高低压涡轮过渡流道型面反问题设计方法
CN108197397B (zh) 一种航空发动机紧固结合面动态性能的优化设计方法
CN106777495B (zh) 通过局部参数化修型控制旋翼翼型俯仰力矩的方法
CN116384007B (zh) 一种当表面存在镀层时镀前花键的设计方法

Legal Events

Date Code Title Description
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