CN106971023B - 基于超椭圆曲线的轮盘异型孔结构设计方法 - Google Patents
基于超椭圆曲线的轮盘异型孔结构设计方法 Download PDFInfo
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000005457 optimization Methods 0.000 claims abstract description 21
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000005540 biological transmission Effects 0.000 abstract description 5
- 238000013178 mathematical model Methods 0.000 abstract description 4
- 238000012938 design process Methods 0.000 abstract 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000005498 polishing Methods 0.000 description 1
Images
Classifications
-
- 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/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical 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
技术领域
本发明涉及航空发动机涡轮转子部件的孔类结构优化设计领域,具体是一种基于超椭圆曲线的轮盘异型孔结构设计方法。
背景技术
航空发动机涡轮转子部件工作环境恶劣,其轮盘幅板上的螺栓孔、通气孔等孔类结构常因应力集中现象,产生孔边裂纹,成为轮盘失效的一个重要原因。传统措施如采用孔边倒角、抛光、增加结构重量等方法并不能从根本上解决该问题。
发明内容
本发明为了解决现有技术的问题,提供了一种基于超椭圆曲线的轮盘异型孔结构设计方法,按结构使用需求降低孔边应力,同时保证孔轮廓形状变化尽可能小,以此满足强度要求、装配要求并保证可靠传力。该方法所需设计变量较少,异型孔数学模型简洁。
本发明包括以下步骤:
步骤一、定义超椭圆异型孔曲线方程。
式中,a、b为超椭圆曲线的半轴值,m、n为指数,其取值不同,超椭圆异型孔形状不同,此发明中取m,n>2,此时异型孔轮廓介于方孔和圆孔之间,更加符合轮盘孔类结构的实际设计要求。
步骤二、采用超椭圆参数方程运用有限元分析软件进行参数化建模:
式中,θ为参数,sgn为符号判断函数:
步骤三、构建超椭圆异型孔多目标优化模型:
minf{f1(σmax),f2(m,n)}; (4)
其中,优化设计变量m、n为超椭圆曲线指数;σmax=σmax(m,n)为设计点λ(m,n)对应的超椭圆孔孔边最大主应力值,目标函数f1(σmax)代表孔边应力降低程度;目标函数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-原螺栓通孔(基圆)。
具体实施方式
下面结合附图对本发明作进一步说明。
下面结合具体实施方式对本发明作更进一步的说明。
首先,定义超椭圆异型孔曲线方程。
本发明中采用如下超椭圆曲线方程对轮盘孔类结构进行改进:
其中a、b为超椭圆曲线的半轴值,m、n为指数,其取值不同,超椭圆异型孔形状不同,此发明中取m,n>2,此时异型孔轮廓介于方孔和圆孔之间,更加符合轮盘孔类结构的实际设计要求。
按方程(2)运用有限元分析软件对超椭圆参数化方程进行建模,进一步对孔边应力情况进行分析:
式中,θ为参数,单位为角度,sgn为符号判断函数:
构建超椭圆异型孔多目标优化模型:
minf{f1(σmax),f2(m,n)}; (4)
其中,设计变量m、n为超椭圆曲线指数,它们代表了设计点λ(m,n);upper(=6~8)为m、n上界,lower(=2)为m、n下界;σmax=σmax(m,n)为设计点λ(m,n)对应的超椭圆孔边最大主应力值。目标函数f1(σmax)代表孔边应力降低率,可由设计者根据需要设定(一般下降率可为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.一种基于超椭圆曲线的轮盘异型孔结构设计方法,其特征在于包括以下步骤:
步骤一、定义超椭圆异型孔曲线方程:
式中,a、b为超椭圆曲线的半轴值,m、n为指数;
步骤二、采用超椭圆参数方程运用有限元分析软件进行参数化建模:
式中,θ为参数,sgn为符号判断函数:
步骤三、构建超椭圆异型孔多目标优化模型:
min f{f1(σmax),f2(m,n)}; (4)
其中,优化设计变量m、n为超椭圆曲线指数;σmax=σmax(m,n)为设计点λ(m,n)对应的超椭圆孔孔边最大主应力值,目标函数f1(σmax)代表孔边应力降低程度;目标函数f2(m,n)代表超椭圆异型孔轮廓变化;upper和lower分别为曲线指数的上、下界;
目标函数f2(m,n)计算方法如方程如下所示:
f2(m,n)=S1-S0; (6)
S1为超椭圆异型孔面积,S0为原圆孔面积;
步骤四、运用多目标优化方法对超椭圆异型孔优化模型进行寻优,最终得到异型孔轮廓优化结果。
2.根据权利要求1所述的一种基于超椭圆曲线的轮盘异型孔结构设计方法,其特征在于:步骤一中,m>2,n>2。
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)
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)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105046000A (zh) * | 2015-07-20 | 2015-11-11 | 南京航空航天大学 | 基于工程离散变量的两变量多圆弧异型孔均衡设计方法及模型 |
Family Cites Families (1)
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 |
-
2017
- 2017-03-01 CN CN201710117739.4A patent/CN106971023B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105046000A (zh) * | 2015-07-20 | 2015-11-11 | 南京航空航天大学 | 基于工程离散变量的两变量多圆弧异型孔均衡设计方法及模型 |
Non-Patent Citations (2)
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 |