CN102996511A - 曲率连续的可控扩散叶型的cad辅助设计方法 - Google Patents

曲率连续的可控扩散叶型的cad辅助设计方法 Download PDF

Info

Publication number
CN102996511A
CN102996511A CN2012105394829A CN201210539482A CN102996511A CN 102996511 A CN102996511 A CN 102996511A CN 2012105394829 A CN2012105394829 A CN 2012105394829A CN 201210539482 A CN201210539482 A CN 201210539482A CN 102996511 A CN102996511 A CN 102996511A
Authority
CN
China
Prior art keywords
point
curvature
suction surface
leading edge
coordinate
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
CN2012105394829A
Other languages
English (en)
Other versions
CN102996511B (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201210539482.9A priority Critical patent/CN102996511B/zh
Publication of CN102996511A publication Critical patent/CN102996511A/zh
Application granted granted Critical
Publication of CN102996511B publication Critical patent/CN102996511B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Moulds For Moulding Plastics Or The Like (AREA)

Abstract

曲率连续的可控扩散型叶型的CAD辅助设计方法属于压气机叶片设计领域,其特征在于,含有:分别利用压力面、吸力面与前缘小圆的切点求前缘点,用三次样条拟合,求出吸力面、压力面上各型值点的曲率,用三次贝塞尔曲线求出曲率连续的前缘的吸力面侧型线和压力面侧型线,用B样条曲线给定吸力面的曲率分布,求取B样条曲线的控制点,使曲率分布对应的曲率连续吸力面型线与原始吸力面型线偏差最小,求取优化后叶型尾缘小圆的位置。本发明与原始可控扩散叶型相比,冲角范围扩大2.5°,且在低损失区域内的总压损失降低。

Description

曲率连续的可控扩散叶型的CAD辅助设计方法
技术领域
本发明涉及压气机叶片的优化设计方法。
背景技术
压气机是航空发动机和燃气轮机的重要部件之一,叶型设计是压气机设计的关键环节。早期的压气机叶型采用系列叶型(如NACA-64系列,C-4系列等)。上世纪五六十年代,可控扩散叶型开始在压气机中得到广泛应用,其与传统系列叶型相比,能够更好地控制扩散、抑制边界层的分离,在高马赫数下降低激波损失,因而能够减小型面损失、扩大攻角范围,且在多级环境下易于匹配。而在压气机级负荷不断提高的趋势下,发展高负荷、高效率、高可靠性的压气机对叶型的损失和攻角范围提出了更高的要求。
目前的可控扩散叶型(如图1),其型面曲率不连续(具体表现为:①压力面或吸力面型线与前缘的切点处曲率不连续,如图2,②压力面或吸力面型线本身曲率不连续,如图3),而曲率对与边界层的发展,特别是压气机叶片吸力面侧强逆压梯度下的边界层发展有明显影响。叶型曲率连续有利于改善吸力面边界层发展,抑制边界层分离,进而降低叶型损失,扩大叶型攻角范围。
一些研究者开展了控制叶型曲率的叶型设计方法的研究。Korakianitis等的论文“Design ofhigh-efficiency turbomachinery blades for energy conversion devices with the three-dimensionalprescribed surface curvature distribution blade design(CIRCLE)method”提出了一种预置曲率的叶型优化设计方法,其采用贝塞尔曲线给定叶片中部型线(即压力面和吸力面型线)的曲率沿轴向的分布,叶片前缘采用常规前缘形状(如圆弧),通过在前缘和中部型线之间添加一段过渡型线(其由一抛物线表达的构造线叠加一厚度分布来获得)来保证前缘与中部型线之间曲率光滑过渡。Sommer等的论文“Curvature driven two-dimensional multi-objective optimizationof compressor blade sections”提出了一种基于曲率的叶型优化设计方法,其将从前缘点到型线与尾缘切点处的型线(即前缘吸力面侧型线加上吸力面型线或是前缘压力面侧型线加上压力面型线)的曲率分布采用一个B样条曲线定义,以叶型损失为目标函数,使用遗传算法对B样条曲线的控制点坐标进行优化来对叶型进行优化设计。Korakianitis的方法由于采用了常规前缘,其前缘与中部型线的过渡区域曲率的梯度仍然较大,另外,该方法的优化过程并没有考虑叶型结构强度的要求。而Sommer的方法将前缘型线与中部型线的曲率分布统一使用一个B样条曲线来表示,由于前缘曲率与中部型线的曲率差别很大,因此Sommer的方法获得的叶型形状与原始叶型形状不可避免地存在较大差别,尽管该方法在优化时将叶片角度和最小截面积等结构气动参数与原始值的偏差进行了一定的限制。宋寅等的论文“曲率连续的压气机叶片前缘设计及其对气动性能的影响”提出使用三次贝塞尔曲线定义前缘吸力面侧和压力面侧型线,但是没有公开设计参数或是贝塞尔曲线控制点的选择方法。
考虑到高负荷压气机对于级匹配和结构强度的苛刻要求,应当发展一种叶型优化设计方法,其具有以下特性:①叶型前缘与吸力面型线切点处曲率连续,且切点附近曲率变化平缓,②吸力面型线曲率连续,③为保证对叶型优化设计后不改变级匹配和结构强度特性,优化后叶型的中部型线与原始型线具有较小的偏差。
发明内容
现代航空发动机或重型燃气轮机的压气机级负荷不断提高,对叶型的损失水平和攻角范围提出了更高的要求。已有的可控扩散叶型(如图1)其型面曲率不连续(具体表现为:①压力面或吸力面型线与前缘的切点处曲率不连续,如图2,②压力面或吸力面型线本身曲率不连续,如图3),而曲率对与边界层的发展,特别是压气机叶片吸力面侧强逆压梯度下的边界层发展有明显影响。叶型曲率连续有利于改善吸力面边界层发展,抑制边界层分离,进而降低叶型损失,扩大叶型攻角范围。为了提高现有可控扩散叶型的性能,本发明提出了曲率连续的可控扩散叶型的设计技术。
本发明的目的在于对可控扩散叶型进行优化设计,使其叶片前缘与吸力面、压力面的切点处曲率连续以及吸力面型线的曲率连续,进而提高叶型的性能(降低叶型损失、扩大叶型攻角范围)。
本发明的特征在于,是在计算机中依次按以下步骤实现的:
步骤(1),计算机初始化:
输入可控扩散叶片,以下简称叶片的以下参数:
前缘的圆心坐标(xLE,yLE)和半径rLE
尾缘的圆心坐标(xTE,yTE)和半径rTE
压力面上均匀分布的N个型值点,记为P1 ps,P2 ps...PN ps,N=500,坐标分别记为
Figure BDA00002579712400023
Figure BDA00002579712400024
...,其中P1 ps为所述压力面与前缘小圆的切点,
Figure BDA00002579712400026
为压力面与尾缘小圆的切点,
吸力面上均匀分布的N个型值点,记为P1 ss,P2 ss...,PN ss坐标分别为
Figure BDA00002579712400029
Figure BDA000025797124000210
...,
Figure BDA000025797124000211
N=500,其中P1 ss为所述吸力面与前缘小圆的切点,
Figure BDA000025797124000212
为所述吸力面与尾缘小圆的切点;
步骤(2),按以下步骤计算前缘点P1的坐标和该前缘点P1对应的圆心角θ1
步骤(2.1),在所述前缘的圆心(xLE,yLE)建立一个直角坐标系,横轴x轴向右为正,
步骤(2.2),按下式计算所述前缘点P1(x1,y1)对应的圆心角θ1以及P1点的斜率k1
x1=xLE+rLE·cosθ1,y1=yLE+rLE·cosθ1,k1=-tanθ1
θ1=(θss_LEps_LE)/2,其中:
θss_LE为吸力面与前缘的切点P1 ss在前缘小圆上的圆心角, θ ss _ LE = arctan ( y 1 ss - y LE x 1 ss - x LE ) + π , 其中:
θps_LE为压力面与前缘的切点P1 ps在前缘小圆上的圆心角, θ ps _ LE = arctan ( y 1 ps - y LE x 1 ps - x LE ) + 2 π ,
步骤(2.3),对所述吸力面和压力面的叶型坐标使用三弯矩法进行样条拟合,求出所述吸力面和压力面各型值点Pi ss、Pi ps处的型线斜率,用
Figure BDA00002579712400033
Figure BDA00002579712400034
(i=1,2,...,N)表示,
步骤(2.4),按下式求出所述吸力面各型值点Pi ss处的流向坐标
Figure BDA00002579712400035
i=1,2,...,N和中间型值点的曲率
Figure BDA00002579712400036
i=2,...,N-1,
s i ss = 0 , i = 1 , s i - 1 ss + ( x i ss - x i - 1 ss ) 2 + ( y i ss - y i - 1 ss ) 2 i = 2,3 , . . . , N ,
C i ss = ( dx · ddy - ddx · dy ) / [ ( dx 2 + dy 2 ) 1.5 ] , 其中:
dx = 1 2 · ( x i + 1 ss - x i ss s i + 1 ss - s i ss + x i ss - x i - 1 ss s i ss - s i - 1 ss ) ,
dy = 1 2 · ( y i + 1 ss - y i ss s i + 1 ss - s i ss + y i ss - y i - 1 ss s i ss - s i - 1 ss ) ,
ddx = 2 s i + 1 ss - s i - 1 ss · ( x i + 1 ss - x i ss s i + 1 ss - s i ss - x i ss - x i - 1 ss s i ss - s i - 1 ss ) ,
ddy = 2 s i + 1 ss - s i - 1 ss · ( y i + 1 ss - y i ss s i + 1 ss - s i ss - y i ss - y i - 1 ss s i ss - s i - 1 ss ) ,
同理,求出所述压力面各型值点处的流向坐标
Figure BDA00002579712400041
i=1,2,...,N以及曲率
Figure BDA00002579712400042
i=2,...,N-1;
步骤(3),使用三次贝塞尔曲线按以下步骤定义新的前缘吸力面侧型线,所述型线是指前缘与中部型线之间添加的一段过渡型线,以便使叶片前缘分别与吸力面、压力面的切点处曲率连续,以及吸力面侧型线的曲率连续:
步骤(3.1),定义三次贝塞尔曲线的四个控制点P1、P2、P3和P4,其中除了P1是所述前缘点外:
P2是原始的吸力面型线上也是优化后的前缘吸力面侧型线上的待定点,
P3是位于前缘点P1的切线和点P2的切线的相交点P0(x0,y0)这一个点和点P1的连线上的一个待定点,
P4是位于前缘点P1的切线和点P2的切线的相交点P0(x0,y0)这一个点和点P2的连线上的一个待定点,
步骤(3.2),按下述逐点迭代的步骤确定控制点P2的位置坐标:
步骤(3.2.1),选择吸力面上的某一型值点
Figure BDA00002579712400043
作为控制点P2的初始值,
Figure BDA00002579712400044
的横坐标xj的取值满足
Figure BDA00002579712400045
1≤j≤N,令点
Figure BDA00002579712400046
的坐标为(x2,y2),x2=xj,斜率为k2,曲率为C2
步骤(3.2.2),设定一组控制点P3的坐标控制系数a1和控制点P4的控制系数a2,0<a1<1,0<a2<1,其中:
控制点P3的坐标:x3=a1·x1+(1-a1)·x0,y3=a1·y1+(1-a1)·y0
控制点P4的坐标:x4=a2·x2+(1-a2)·x0,y4=a2·y2+(1-a2)·y0
步骤(323),按步骤(322)得到的四个控制点生成三次贝塞尔曲线,并将该贝塞尔曲线插值为50个点,记为Pi LE_ss,i=1,2,...,50,序号沿P1向P2的方向取,坐标为
Figure BDA00002579712400047
计算中间型值点Pi LE_ss的曲率
Figure BDA00002579712400048
i=2,...,49,并计算
Figure BDA00002579712400049
与控制点P2处吸力面型线曲率C2的偏差 error _ C = | C 49 LE _ ss - C 2 | ,
步骤(3.2.4),判断步骤(3.2.3)得到的error_C的值:
若error_C<0.1,则优化后的前缘吸力面侧型线满足要求,
若error_C≥0.1,则把控制点P′2选为吸力面上的下一个型值点
Figure BDA00002579712400051
步骤(3.2.5),重复步骤(3.2.1)~步骤(3.2.4),直到error_C<0.1为止,得到了一条由P1、P′3、P′4、P′2四个控制点确定的三次贝塞尔曲线,重新离散后的型值点为Pi LE_ss′(i=1,3,...,50),共50个,控制点P′2为新的前缘与吸力面的切点,并且满足在吸力面与前缘的切点处曲率的连续性;
步骤(4),按步骤(3)所述的方法求出压力面与前缘切点处曲率连续的前缘压力面侧型线,离散后的型值点为Pi LE_ps(i=1,3,...,50),共50个;
步骤(5),以步骤(3)得到的优化后的前缘与吸力面的切点P′2为起点,按如下步骤求取曲率连续的吸力面型线:
步骤(5.1),设定:使用7个控制点的3次B样条曲线来表达曲率连续的吸力面型线的曲率分布C(s),该B样条曲线的节点矢量取为U=[0,0,0,0,0.25,0.5,0.75,1,1,1,1]T
步骤(5.2),把所述曲率分布曲线C(s)离散为N个点,用(s′1,C′1),(s′2,C′2),…,(s′N,C′N)表示,N=500,i=1,2,...,N,C′i为曲率,
步骤(5.3),按下式求出各离散点对应的曲线倾斜角
s′i为B样条曲线上的序号为i的离散点的流向坐标,k′2为优化后的曲率连续的前缘与吸力面的切点P′2的斜率,
步骤(5.4),按下式求出各离散点的横坐标和纵坐标(X′i,Y′i):
Figure BDA00002579712400054
Figure BDA00002579712400055
(x′2,y′2)为优化后的曲率连续的前缘与吸力面的切点P′2的坐标,
步骤(5.5),使用步骤(5.4)得到的各离散点的坐标进行三次样条插值,得到与原始吸力面型线的有效型值点的横坐标对应的曲率连续的吸力面型线的纵坐标
Figure BDA00002579712400061
从而得到优化后的曲率连续的吸力面型线的型值点的坐标,同时求出曲率分布确定的曲率连续的吸力面型线与原始吸力面型线的偏差 error _ ss = max ( | y i ss ′ - y i ss | ) ,
Figure BDA00002579712400063
为与优化后的控制点P′2对应的离散点的序号,
步骤(5.6),使用优化软件iSIGHT获得最优的曲率分布控制点:
设定:曲率分布曲线的控制点为CP1~CP7,坐标分别为:(sCP1,CCP1),(sCP2,CCP2),...,(sCP7,CCP7),其中:
s CP 1 = s i P 2 ′ ss , C CP 1 = C i P 2 ′ ss , s CP 7 = S N - 1 ss , C CP 7 = C N - 1 ss ,
Figure BDA00002579712400069
为吸力面上对应于控制点P′2的序号为
Figure BDA000025797124000610
的离散点处的流向坐标
Figure BDA000025797124000611
Figure BDA000025797124000612
为吸力面上对应于控制点P′2的序号为
Figure BDA000025797124000613
的离散点处的曲率,
Figure BDA000025797124000614
分别为原始吸力面型线倒数第二个型值点的流向坐标和曲率,
设定:优化变量为sCP2~sCP6,CCP2~CCP6
设定:变量的约束为sCP1≤sCP2,sCP3,sCP4,sCP5,sCP6≤sCP7,以及-10≤CCP2,CCP3,CCP4,CCP5,CCP6≤10,
设定:目标函数为曲率连续的吸力面型线与原始吸力面型线的偏差error_ss最小,
设定:优化算法为序列二次规划方法NLPQL,
得到:error_ss最小的坐标为
Figure BDA000025797124000616
...,
Figure BDA000025797124000617
所对应的型值点
Figure BDA000025797124000618
...,
Figure BDA000025797124000619
步骤(6),根据曲率连续的吸力面型线和原始压力面型线确定尾缘小圆以及尾缘与吸力面和压力面分别相切的切点:
步骤(6.1),确定下述三条直线:
a,曲率连续的吸力面型线的型值点
Figure BDA000025797124000620
Figure BDA000025797124000621
组成的直线,
b,原始压力面型线的型值点
Figure BDA000025797124000622
Figure BDA000025797124000623
组成的直线,
c,原始尾缘小圆的切线,切点为原始前缘小圆圆心与原始尾缘小圆圆心的连线的延长线与原始尾缘小圆的交点,
步骤(6.2),确定同时与步骤(6.1)中的三条直线a,b,c相切的圆即为尾缘小圆,
步骤(6.3),根据步骤(6.2)的结果计算出尾缘小圆的圆心坐标(xTE′,yTE′)和半径rTE′,
步骤(6.4),确定:
尾缘小圆与直线a的交点为尾缘与吸力面的切点,
尾缘小圆与直线b的交点为尾缘与压力面的切点。
本发明的优点在于,优化后的曲率连续的可控扩散叶型与原始可控扩散叶型相比,冲角范围(即总压损失小于2倍的最小总压损失的冲角范围)扩大,且在低损失区域内的总压损失降低。
附图说明
图1为可控扩散叶型的示意图,1为前缘(为一圆弧),记为LE;2为吸力面,记为ss;3为压力面,记为ps;4为尾缘(为一圆弧),记为TE;5为吸力面与前缘小圆的切点;6为压力面与前缘小圆的切点;7为吸力面与尾缘小圆的切点;8为压力面与尾缘小圆的切点。
图2为图1所示的叶型前缘的吸力面侧型线的曲率沿型线长度方向的分布。
图3为图1所示的叶型吸力面型线的曲率沿型线长度方向的分布。
图4为叶片前缘点及其圆心角的示意图,9为前缘点。
图5为曲率连续的前缘吸力面侧型线的定义方式的示意图。
图6为前缘吸力面侧型线优化前后曲率比较:——为原始叶型前缘的曲率,----·为优化后曲率连续前缘的曲率。
图7为原始叶型吸力面曲率分布与用B样条表达的优化后曲率分布的比较:----·为原始曲率分布,——为优化后曲率分布,
Figure BDA00002579712400071
为B样条的控制点。
图8为iSIGHT优化求取B样条控制点的过程的示意图。
图9为曲率连续叶型的尾缘小圆的示意图。
图10为优化设计前后的叶型的比较:——为原始的可控扩散叶型,----·为优化后的曲率连续的可控扩散叶型。
图11为优化设计前后叶型的冲角-损失特性:——为原始叶型的冲角-损失特性,----·为优化后曲率连续叶型的的冲角-损失特性,…………为2倍的最小总压损失。
图12为本发明的程序流程框图。
具体实施方式
本发明的思路如下:
1.根据吸力面与前缘小圆的切点、压力面与前缘小圆的切点求取前缘点,
2.将吸力面与压力面进行三次样条拟合,求出各型值点的曲线斜率,
3.使用三次贝塞尔曲线给定曲率连续的前缘吸力面侧型线和前缘压力面侧型线,其控制点通过给定的初值选择和迭代方法获得,
5.使用B样条曲线给定吸力面的曲率分布,求取B样条曲线的控制点,使得该曲率分布曲线对应的曲率连续吸力面型线与原始吸力面型线的偏差最小(即获得原始吸力面型线的曲率连续拟合),
6.根据曲率连续的吸力面型线和原始压力面型线,结合原始尾缘小圆,求取优化后叶型尾缘小圆的位置。
为了更好地说明本发明的目的和优点,下面结合附图和实施例对本发明的内容作进一步说明。
本实施例根据发明内容所述的设计方法对一可控扩散叶型进行优化设计,并用计算流体力学方法验证其作用效果。
步骤(1),计算机初始化:
输入可控扩散叶片,以下简称叶片的以下参数:
前缘的圆心坐标(xLE,yLE)和半径rLE
尾缘的圆心坐标(xTE,yTE)和半径rTE
压力面上均匀分布的500个型值点,记为P1 ps,P2 ps,...
Figure BDA00002579712400082
坐标分别记为
Figure BDA00002579712400083
Figure BDA00002579712400084
...,
Figure BDA00002579712400085
其中P1 ps为所述压力面与前缘小圆的切点,
Figure BDA00002579712400086
为压力面与尾缘小圆的切点,
吸力面上均匀分布的500个型值点,记为P1 ss
Figure BDA00002579712400087
...,
Figure BDA00002579712400088
坐标分别为
Figure BDA000025797124000810
...,其中P1 ss为所述吸力面与前缘小圆的切点,
Figure BDA000025797124000812
为所述吸力面与尾缘小圆的切点;
步骤(2),如图4,按以下步骤计算前缘点P1的坐标和该前缘点P1对应的圆心角θ1
步骤(2.1),在所述前缘的圆心(xLE,yLE)建立一个直角坐标系,横轴x轴向右为正,
步骤(2.2),按下式计算所述前缘点P1(x1,y1)对应的圆心角θ1以及P1点的斜率k1
x1=xLE+rLE·cosθ1,y1=yLE+rLE·cosθ1,k1=-tanθ1
θ1=(θss_LEps_LE)/2,其中:
θss_LE为吸力面与前缘的切点P1 ss在前缘小圆上的圆心角, θ ss _ LE = arctan ( y 1 ss - y LE x 1 ss - x LE ) + π , 其中:
θps_LE为压力面与前缘的切点P1 ps在前缘小圆上的圆心角, θ ps _ LE = arctan ( y 1 ps - y LE x 1 PS - x LE ) + 2 π ,
步骤(2.3),对所述吸力面和压力面的叶型坐标使用三弯矩法进行样条拟合,求出所述吸力面和压力面各型值点Pi ss、Pi ps处的型线斜率,用
Figure BDA00002579712400093
(i=1,2,...,500)表示,
步骤(2.4),按下式求出所述吸力面各型值点Pi ss处的流向坐标
Figure BDA00002579712400095
i=1,2,...,500和中间型值点的曲率
Figure BDA00002579712400096
i=2,...,499,
s i ss = 0 , i = 1 , s i - 1 ss + ( x i ss - x i - 1 ss ) 2 + ( y i ss - y i - 1 ss ) 2 , i = 2,3 , . . . , 500 ,
C i ss = ( dx · ddy - ddx · dy ) / [ ( dx 2 + dy 2 ) 1.5 ] , 其中:
dx = 1 2 · ( x i + 1 ss - x i ss s i + 1 ss - s i ss + x i ss - x i - 1 ss s i ss - s i - 1 ss ) ,
dy = 1 2 · ( y i + 1 ss - y i ss s i + 1 ss - s i ss + y i ss - y i - 1 ss s i ss - s i - 1 ss ) ,
ddx = 2 s i + 1 ss - s i - 1 ss · ( x i + 1 ss - x i ss s i + 1 ss - s i ss - x i ss - x i - 1 ss s i ss - s i - 1 ss ) ,
ddy = 2 s i + 1 ss - s i - 1 ss · ( y i + 1 ss - y i ss s i + 1 ss - s i ss - y i ss - y i - 1 ss s i ss - s i - 1 ss ) ,
同理,求出所述压力面各型值点处的流向坐标
Figure BDA000025797124000913
i=1,2,...,500,以及曲率
Figure BDA000025797124000914
i=2,...,499;
步骤(3),如图5,使用三次贝塞尔曲线按以下步骤定义新的前缘吸力面侧型线,所述型线是指前缘与中部型线之间添加的一段过渡型线,以便使叶片前缘分别与吸力面、压力面的切点处曲率连续,以及吸力面侧型线的曲率连续:
步骤(3.1),定义三次贝塞尔曲线的四个控制点P1、P2、P3和P4,其中除了P1是所述前缘点外:
P2是原始的吸力面型线上也是优化后的前缘吸力面侧型线上的待定点,
P3是位于前缘点P1的切线和点P2的切线的相交点P0(x0,y0)这一个点和点P1的连线上的一个待定点,
P4是位于前缘点P1的切线和点P2的切线的相交点P0(x0,y0)这一个点和点P2的连线上的一个待定点,
步骤(3.2),按下述逐点迭代的步骤确定控制点P2的位置坐标:
步骤(3.2.1),选择吸力面上的某一型值点作为控制点P2的初始值,的横坐标xj的取值满足1≤j≤500,令点
Figure BDA00002579712400104
的坐标为(x2,y2),x2=xj,斜率为k2,曲率为C2
步骤(3.2.2),设定一组控制点P3的坐标控制系数a1和控制点P4的控制系数a2,其中:
控制点P3的坐标:x3=a1·x1+(1-a1)·x0,y3=a1·y1+(1-a1)·y0
控制点P4的坐标:x4=a2·x2+(1-a2)·x0,y4=a2·y2+(1-a2)·y0
在本实施例中,在0<a1<1,0<a2<1的范围内取a1=0.75,a2=0.1,
步骤(323),按步骤(322)得到的四个控制点生成三次贝塞尔曲线,并将该贝塞尔曲线插值为50个点,记为Pi LE_ss,i=1,2,...,50,序号沿P1向P2的方向取,坐标为
Figure BDA00002579712400105
计算中间型值点Pi LE_ss的曲率
Figure BDA00002579712400106
i=2,...,49,并计算
Figure BDA00002579712400107
与控制点P2处吸力面型线曲率C2的偏差 error _ C = | C 49 LE _ ss - C 2 | ,
步骤(3.2.4),判断步骤(3.2.3)得到的error_C的值:
若error_C<0.1,则优化后的前缘吸力面侧型线满足要求,
若error_C≥0.1,则把控制点P′2选为吸力面上的下一个型值点
Figure BDA00002579712400109
步骤(3.2.5),重复步骤(3.2.1)~步骤(3.2.4),直到error_C<0.1为止,得到了一条由P1、P′3、P′4、P′2四个控制点确定的三次贝塞尔曲线,重新离散后的型值点为Pi LE_ss′(i=1,3,...,50),共50个,控制点P′2为新的前缘与吸力面的切点,并且满足在吸力面与前缘的切点处曲率的连续性,在本实施例中,
Figure BDA00002579712400111
为满足要求的控制点,则控制点P′2即为
步骤(4),按步骤(3)所述的方法求出压力面与前缘切点处曲率连续的前缘压力面侧型线,离散后的型值点为Pi LE_ps(i=1,3,...,50),共50个;
步骤(5),以步骤(3)得到的优化后的前缘与吸力面的切点P′2为起点,按如下步骤求取曲率连续的吸力面型线:
步骤(5.1),设定:如图7,使用7个控制点的3次B样条曲线来表达曲率连续的吸力面型线的曲率分布C(s),该B样条曲线的节点矢量取为U=[0,0,0,0,0.25,0.5,0.75,1,1,1,1]T
步骤(5.2),把所述曲率分布曲线C(s)离散为500个点,用(s′1,C′1),(s′2,C′2),...,(s′500,C′500)表示,i=1,2,...,500,C′i为曲率,
步骤(5.3),按下式求出各离散点对应的曲线倾斜角
Figure BDA00002579712400113
Figure BDA00002579712400114
s′i为B样条曲线上的序号为i的离散点的流向坐标,k′2为优化后的曲率连续的前缘与吸力面的切点P′2的斜率,
步骤(5.4),按下式求出各离散点的横坐标和纵坐标(X′i,Y′i):
Figure BDA00002579712400115
Figure BDA00002579712400116
(x′2,y′2)为优化后的曲率连续的前缘与吸力面的切点P′2的坐标,
步骤(5.5),使用步骤(5.4)得到的各离散点的坐标进行三次样条插值,得到与原始吸力面型线的有效型值点的横坐标对应的曲率连续的吸力面型线的纵坐标
Figure BDA00002579712400117
从而得到优化后的曲率连续的吸力面型线的型值点的坐标,同时求出曲率分布确定的曲率连续的吸力面型线与原始吸力面型线的偏差 error _ ss = max ( | y i ss ′ - y i ss | ) , 15≤i≤500,
步骤(5.6),使用优化软件iSIGHT获得最优的曲率分布控制点:
设定:曲率分布曲线的控制点为CP1~CP7,坐标分别为:(sCP1,CCP1),(sCP2,CCP2),...,(sCP7,CCP7),其中:
s CP 1 = s 15 ss , C CP 1 = C 15 ss , s CP 7 = s 499 ss , C CP 7 = C 499 ss ,
Figure BDA00002579712400126
为吸力面上对应于控制点P′2的离散点
Figure BDA00002579712400127
处的流向坐标,
Figure BDA00002579712400128
为吸力面上对应于控制点P′2的离散点
Figure BDA00002579712400129
处的曲率,
Figure BDA000025797124001210
分别为原始吸力面型线倒数第二个型值点的流向坐标和曲率,
设定:优化变量为sCP2~sCP6,CCP2~CCP6
设定:变量的约束为sCP1≤sCP2,sCP3,sCP4,sCP5,sCP6≤sCP7,以及-10≤CCP2,CCP3,CCP4,CCP5,CCP6≤10,
设定:目标函数为曲率连续的吸力面型线与原始吸力面型线的偏差error_ss最小,
设定:优化算法为序列二次规划方法NLPQL,
本实施例iSIGHT进行28步迭代后得到:error_ss最小的坐标为
Figure BDA000025797124001212
Figure BDA000025797124001213
...,
Figure BDA000025797124001214
所对应的型值点
Figure BDA000025797124001215
Figure BDA000025797124001216
...,
步骤(6),如图9,根据曲率连续的吸力面型线和原始压力面型线确定尾缘小圆以及尾缘与吸力面和压力面分别相切的切点:
步骤(6.1),确定下述三条直线:
a,曲率连续的吸力面型线的型值点
Figure BDA000025797124001218
组成的直线,
b,原始压力面型线的型值点
Figure BDA000025797124001220
Figure BDA000025797124001221
组成的直线,
c,原始尾缘小圆的切线,切点为原始前缘小圆圆心与原始尾缘小圆圆心的连线的延长线与原始尾缘小圆的交点,
步骤(6.2),确定同时与步骤(6.1)中的三条直线a,b,c相切的圆即为尾缘小圆,
步骤(6.3),根据步骤(6.2)的结果计算出尾缘小圆的圆心坐标(xTE′,yTE′)和半径rTE′,
步骤(6.4),确定:
尾缘小圆与直线a的交点为尾缘与吸力面的切点,
尾缘小圆与直线b的交点为尾缘与压力面的切点。
至此,如图10,求出了优化后的叶型全部几何信息(包括前缘吸力面侧型线,前缘压力面侧型线,吸力面型线,压力面型线以及尾缘型线),优化后的叶型吸力面侧曲率连续。
为验证上述曲率连续叶型优化设计技术的效果,使用计算流体力学软件CFX对原始可控扩散叶型和曲率连续的可控扩散叶型的叶型冲角-损失特性进行计算,计算结果如图11,曲率连续的可控扩散叶型与原始可控扩散叶型相比,冲角范围(即总压损失小于2倍的最小总压损失的冲角范围)扩大了2.5°,且在低损失区域内的总压损失降低。

Claims (1)

1.曲率连续的可控扩散叶型的CAD辅助设计方法,其特征在于,是在计算机中依次按以下步骤实现的:
步骤(1),计算机初始化:
输入可控扩散叶片,以下简称叶片的以下参数:
前缘的圆心坐标(xLE,yLE)和半径rLE
尾缘的圆心坐标(xTE,yTE)和半径rTE
压力面上均匀分布的N个型值点,记为P1 ps,P2 ps,...PN ps,N=500,坐标分别记为
Figure FDA00002579712300013
Figure FDA00002579712300014
...,
Figure FDA00002579712300015
其中P1 ps为所述压力面与前缘小圆的切点,
Figure FDA00002579712300016
为压力面与尾缘小圆的切点,
吸力面上均匀分布的N个型值点,记为P1 ss,P2 ss...,PN ss 坐标分别为
Figure FDA00002579712300019
Figure FDA000025797123000110
...,
Figure FDA000025797123000111
N=500,其中P1 ss为所述吸力面与前缘小圆的切点,PN ss为所述吸力面与尾缘小圆的切点;
步骤(2),按以下步骤计算前缘点P1的坐标和该前缘点P1对应的圆心角θ1
步骤(2.1),在所述前缘的圆心(xLE,yLE)建立一个直角坐标系,横轴x轴向右为正,
步骤(2.2),按下式计算所述前缘点P1(x1,y1)对应的圆心角θ1以及P1点的斜率k1
x1=xLE+rLE·cosθ1,y1=yLE+rLE·cosθ1,k1=-tanθ1
θ1=(θss_LEps_LE)/2,其中:
θss_LE为吸力面与前缘的切点P1 ss在前缘小圆上的圆心角, θ ss _ LE = arctan ( y 1 ss - y LE x 1 ss - x LE ) + π , 其中:
θps_LE为压力面与前缘的切点P1 ps在前缘小圆上的圆心角, θ ps _ LE = arctan ( y 1 ps - y LE x 1 ps - x LE ) + 2 π ,
步骤(2.3),对所述吸力面和压力面的叶型坐标使用三弯矩法进行样条拟合,求出所述吸力面和压力面各型值点Pi ss、Pi ps处的型线斜率,用
Figure FDA000025797123000115
Figure FDA000025797123000116
(i=1,2,...,N)表示,
步骤(2.4),按下式求出所述吸力面各型值点Pi ss处的流向坐标
Figure FDA000025797123000117
i=1,2,...,N和中间型值点的曲率
Figure FDA00002579712300021
i=2,...,N-1,
s i ss = 0 , i = 1 , s i - 1 ss + ( x i ss - x i - 1 ss ) 2 + ( y i ss - y i - 1 ss ) 2 , i = 2,3 , . . . , N ,
C i ss = ( dx · ddy - ddx · dy ) / [ ( dx 2 + dy 2 ) 1.5 ] , 其中:
dx = 1 2 · ( x i + 1 ss - x i ss s i + 1 ss - s i ss + x i ss - x i - 1 ss s i ss - s i - 1 ss ) ,
dy = 1 2 · ( y i + 1 ss - y i ss s i + 1 ss - s i ss + y i ss - y i - 1 ss s i ss - s i - 1 ss ) ,
ddx = 2 s i + 1 ss - s i - 1 ss · ( x i + 1 ss - x i ss s i + 1 ss - s i ss - x i ss - x i - 1 ss s i ss - s i - 1 ss ) ,
ddy = 2 s i + 1 ss - s i - 1 ss · ( y i + 1 ss - y i ss s i + 1 ss - s i ss - y i ss - y i - 1 ss s i ss - s i - 1 ss ) ,
同理,求出所述压力面各型值点处的流向坐标
Figure FDA00002579712300028
i=1,2,...,N以及曲率
Figure FDA00002579712300029
i=2,...,N-1;
步骤(3),使用三次贝塞尔曲线按以下步骤定义新的前缘吸力面侧型线,所述型线是指前缘与中部型线之间添加的一段过渡型线,以便使叶片前缘分别与吸力面、压力面的切点处曲率连续,以及吸力面侧型线的曲率连续:
步骤(3.1),定义三次贝塞尔曲线的四个控制点P1、P2、P3和P4,其中除了P1是所述前缘点外:
P2是原始的吸力面型线上也是优化后的前缘吸力面侧型线上的待定点,
P3是位于前缘点P1的切线和点P2的切线的相交点P0(x0,y0)这一个点和点P1的连线上的一个待定点,
P4是位于前缘点P1的切线和点P2的切线的相交点P0(x0,y0)这一个点和点P2的连线上的一个待定点,
步骤(3.2),按下述逐点迭代的步骤确定控制点P2的位置坐标:
步骤(3.2.1),选择吸力面上的某一型值点
Figure FDA000025797123000210
作为控制点P2的初始值,
Figure FDA000025797123000211
的横坐标xj的取值满足
Figure FDA000025797123000212
1≤j≤N,令点
Figure FDA000025797123000213
的坐标为(x2,y2),x2=xj,斜率为k2,曲率为C2
步骤(3.2.2),设定一组控制点P3的坐标控制系数a1和控制点P4的控制系数a2,0<a1<1,0<a2<1,其中:
控制点P3的坐标:x3=a1·x1+(1-a1)·x0,y3=a1·y1+(1-a1)·y0
控制点P4的坐标:x4=a2·x2+(1-a2)·x0,y4=a2·y2+(1-a2)·y0
步骤(3.2.3),按步骤(3.2.2)得到的四个控制点生成三次贝塞尔曲线,并将该贝塞尔曲线插值为50个点,记为Pi LE_ss,i=1,2,...,50,序号沿P1向P2的方向取,坐标为计算中间型值点Pi LE_ss的曲率
Figure FDA00002579712300032
i=2,...,49,并计算
Figure FDA00002579712300033
与控制点P2处吸力面型线曲率C2的偏差 error _ C = | C 49 LE _ ss - C 2 | ,
步骤(3.2.4),判断步骤(3.2.3)得到的error_C的值:
若error_C<0.1,则优化后的前缘吸力面侧型线满足要求,
若error_1≥0.1,则把控制点P′2选为吸力面上的下一个型值点
Figure FDA00002579712300035
步骤(3.2.5),重复步骤(3.2.1)~步骤(3.2.4),直到error_C<0.1为止,得到了一条由P1、P′3、P′4、P′2四个控制点确定的三次贝塞尔曲线,重新离散后的型值点为Pi LE_ss′(i=1,3,...,50),共50个,控制点P′2为新的前缘与吸力面的切点,并且满足在吸力面与前缘的切点处曲率的连续性;
步骤(4),按步骤(3)所述的方法求出压力面与前缘切点处曲率连续的前缘压力面侧型线,离散后的型值点为Pi LE_ps(i=1,3,...,50),共50个;
步骤(5),以步骤(3)得到的优化后的前缘与吸力面的切点P′2为起点,按如下步骤求取曲率连续的吸力面型线:
步骤(5.1),设定:使用7个控制点的3次B样条曲线来表达曲率连续的吸力面型线的曲率分布C(s),该B样条曲线的节点矢量取为U=[0,0,0,0,0.25,0.5,0.75,1,1,1,1]T
步骤(5.2),把所述曲率分布曲线C(s)离散为N个点,用(s′1,C′1),(s′2,C′2),…,(s′N,C′N)表示,N=500,i=1,2,...,N,C′i为曲率,
步骤(5.3),按下式求出各离散点对应的曲线倾斜角
Figure FDA00002579712300036
Figure FDA00002579712300041
s′i为B样条曲线上的序号为i的离散点的流向坐标,k′2为优化后的曲率连续的前缘与吸力面的切点P′2的斜率,
步骤(5.4),按下式求出各离散点的横坐标和纵坐标(X′i,Y′i):
Figure FDA00002579712300042
Figure FDA00002579712300043
(x′2,y′2)为优化后的曲率连续的前缘与吸力面的切点P′2的坐标,
步骤(5.5),使用步骤(5.4)得到的各离散点的坐标进行三次样条插值,得到与原始吸力面型线的有效型值点的横坐标对应的曲率连续的吸力面型线的纵坐标
Figure FDA00002579712300044
从而得到优化后的曲率连续的吸力面型线的型值点的坐标,同时求出曲率分布确定的曲率连续的吸力面型线与原始吸力面型线的偏差 error _ ss = max ( | y i ss ′ - y i ss | ) ,
Figure FDA00002579712300046
Figure FDA00002579712300047
为与优化后的控制点P′2对应的离散点的序号,
步骤(5.6),使用优化软件iSIGHT获得最优的曲率分布控制点:
设定:曲率分布曲线的控制点为CP1~CP7,坐标分别为:(sCP1,CCP1),(sCP2,CCP2),...,(sCP7,CCP7),其中:
s CP 1 = s i P 2 ′ ss , C CP 1 = C i P 2 ′ ss , s CP 7 = S N - 1 ss , C CP 7 = C N - 1 ss ,
Figure FDA000025797123000412
为吸力面上对应于控制点P′2的序号为
Figure FDA000025797123000413
的离散点处的流向坐标
Figure FDA000025797123000414
Figure FDA000025797123000415
为吸力面上对应于控制点P′2的序号为
Figure FDA000025797123000416
的离散点处的曲率,
Figure FDA000025797123000418
分别为原始吸力面型线倒数第二个型值点的流向坐标和曲率,
设定:优化变量为sCP2~sCP6,CCP2~CCP6
设定:变量的约束为sCP1≤sCP2,sCP3,sCP4,sCP5,sCP6≤sCP7,以及-10≤CCP2,CCP3,CCP4,CCP5,CCP6≤10,
设定:目标函数为曲率连续的吸力面型线与原始吸力面型线的偏差error_ss最小,
设定:优化算法为序列二次规划方法NLPQL,
得到:error_ss最小的坐标为...,
Figure FDA00002579712300052
所对应的型值点
Figure FDA00002579712300053
...,
Figure FDA00002579712300054
步骤(6),根据曲率连续的吸力面型线和原始压力面型线确定尾缘小圆以及尾缘与吸力面和压力面分别相切的切点:
步骤(6.1),确定下述三条直线:
a,曲率连续的吸力面型线的型值点
Figure FDA00002579712300055
组成的直线,
b,原始压力面型线的型值点
Figure FDA00002579712300057
Figure FDA00002579712300058
组成的直线,
c,原始尾缘小圆的切线,切点为原始前缘小圆圆心与原始尾缘小圆圆心的连线的延长线与原始尾缘小圆的交点,
步骤(6.2),确定同时与步骤(6.1)中的三条直线a,b,c相切的圆即为尾缘小圆,
步骤(6.3),根据步骤(6.2)的结果计算出尾缘小圆的圆心坐标(xTE′,yTE′)和半径rTE′,
步骤(6.4),确定:
尾缘小圆与直线a的交点为尾缘与吸力面的切点,
尾缘小圆与直线b的交点为尾缘与压力面的切点。
CN201210539482.9A 2012-12-13 2012-12-13 曲率连续的可控扩散叶型的cad辅助设计方法 Active CN102996511B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210539482.9A CN102996511B (zh) 2012-12-13 2012-12-13 曲率连续的可控扩散叶型的cad辅助设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210539482.9A CN102996511B (zh) 2012-12-13 2012-12-13 曲率连续的可控扩散叶型的cad辅助设计方法

Publications (2)

Publication Number Publication Date
CN102996511A true CN102996511A (zh) 2013-03-27
CN102996511B CN102996511B (zh) 2015-04-29

Family

ID=47925492

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210539482.9A Active CN102996511B (zh) 2012-12-13 2012-12-13 曲率连续的可控扩散叶型的cad辅助设计方法

Country Status (1)

Country Link
CN (1) CN102996511B (zh)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103398021A (zh) * 2013-07-20 2013-11-20 大连理工大学 一种低转速高压比贝塞尔中弧线强扭曲叶片轴流叶轮
CN103527514A (zh) * 2013-10-29 2014-01-22 中国北车集团大连机车研究所有限公司 机车用离心通风机叶轮
CN103541774A (zh) * 2013-11-14 2014-01-29 上海汽轮机厂有限公司 涡轮叶片设计方法
CN104613157A (zh) * 2015-01-26 2015-05-13 北京理工大学 基于贝塞尔曲线的液力变矩器二维叶片型线构造方法
CN106021681A (zh) * 2016-05-13 2016-10-12 北京航空航天大学 一种压气机叶型优化方法及装置
CN106021782A (zh) * 2016-05-31 2016-10-12 西北工业大学 基于中弧线的叶片前后缘拟合及截面线光滑重构方法
CN106151105A (zh) * 2016-08-03 2016-11-23 华中科技大学 一种轴流叶型
CN106351878A (zh) * 2016-09-28 2017-01-25 华中科技大学 一种轴流掠形叶片
CN107657080A (zh) * 2017-08-30 2018-02-02 中国科学院深海科学与工程研究所 一种基于自由曲面的深海照明灯设计方法及设计系统
CN107725477A (zh) * 2017-10-10 2018-02-23 北京航空航天大学 一种优化吸力面波系抑制风扇激波噪声的前缘设计方法
CN108916113A (zh) * 2018-06-13 2018-11-30 中国北方发动机研究所(天津) 一种直纹面压气机叶轮叶片曲面的调整方法
CN109083798A (zh) * 2017-06-13 2018-12-25 国网江苏省电力公司常州供电公司 流体发电装置
CN109571996A (zh) * 2018-12-04 2019-04-05 北京宇航系统工程研究所 样条曲线回转体正置正交复合材料网格结构硬模成型装置
CN110566503A (zh) * 2019-10-10 2019-12-13 珠海格力电器股份有限公司 离心风叶、离心风机、空调器和交通运输工具
CN110593960A (zh) * 2019-09-19 2019-12-20 西安交通大学 一种可用于弯扭掠的轴流透平机械叶片参数化方法
CN110727995A (zh) * 2018-06-29 2020-01-24 中国航发商用航空发动机有限责任公司 叶片形状的构造方法、叶片以及计算机设备
CN111079239A (zh) * 2019-12-19 2020-04-28 中国航空发动机研究院 一种仿生压气机叶栅造型方法
CN111400834A (zh) * 2020-03-20 2020-07-10 国电联合动力技术有限公司 风力发电机组叶片翼型气动优化设计方法、模型及装置
CN111734678A (zh) * 2020-06-24 2020-10-02 西北工业大学 一种压气机叶型非对称前缘设计方法
CN113221264A (zh) * 2021-04-21 2021-08-06 山东双轮股份有限公司 海水淡化泵流道式导叶结构设计优化方法
CN113883093A (zh) * 2021-09-22 2022-01-04 哈尔滨工业大学 一种低反力度压气机叶片设计方法、动叶及压气机
US11333167B2 (en) 2017-04-17 2022-05-17 Ihi Corporation Method of designing blade of axial flow fluid machine and blade
CN117763878A (zh) * 2024-02-22 2024-03-26 清华大学 压气机叶型的确定方法、装置及存储介质

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107178526B (zh) * 2017-07-11 2019-09-24 上海电气燃气轮机有限公司 轴流压气机叶型设计方法及轴流压气机叶型

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201159202Y (zh) * 2008-02-28 2008-12-03 大连海事大学 一种压气机动叶叶尖小翼
JP2011017290A (ja) * 2009-07-09 2011-01-27 Mitsubishi Heavy Ind Ltd 翼体および回転機械

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201159202Y (zh) * 2008-02-28 2008-12-03 大连海事大学 一种压气机动叶叶尖小翼
JP2011017290A (ja) * 2009-07-09 2011-01-27 Mitsubishi Heavy Ind Ltd 翼体および回転機械

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李想等: "某离心叶轮的CFD分析和结果确认", 《工程热物理学报》 *
马文生等: "叶顶间隙对压气机性能的影响", 《动力工程》 *

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103398021A (zh) * 2013-07-20 2013-11-20 大连理工大学 一种低转速高压比贝塞尔中弧线强扭曲叶片轴流叶轮
CN103527514A (zh) * 2013-10-29 2014-01-22 中国北车集团大连机车研究所有限公司 机车用离心通风机叶轮
CN103541774A (zh) * 2013-11-14 2014-01-29 上海汽轮机厂有限公司 涡轮叶片设计方法
CN103541774B (zh) * 2013-11-14 2015-06-17 上海汽轮机厂有限公司 涡轮叶片设计方法
CN104613157A (zh) * 2015-01-26 2015-05-13 北京理工大学 基于贝塞尔曲线的液力变矩器二维叶片型线构造方法
CN106021681B (zh) * 2016-05-13 2019-06-25 北京航空航天大学 一种压气机叶型优化方法及装置
CN106021681A (zh) * 2016-05-13 2016-10-12 北京航空航天大学 一种压气机叶型优化方法及装置
CN106021782A (zh) * 2016-05-31 2016-10-12 西北工业大学 基于中弧线的叶片前后缘拟合及截面线光滑重构方法
CN106151105A (zh) * 2016-08-03 2016-11-23 华中科技大学 一种轴流叶型
CN106151105B (zh) * 2016-08-03 2018-07-06 华中科技大学 一种轴流叶型
CN106351878A (zh) * 2016-09-28 2017-01-25 华中科技大学 一种轴流掠形叶片
US11333167B2 (en) 2017-04-17 2022-05-17 Ihi Corporation Method of designing blade of axial flow fluid machine and blade
CN109083798A (zh) * 2017-06-13 2018-12-25 国网江苏省电力公司常州供电公司 流体发电装置
CN109083798B (zh) * 2017-06-13 2024-02-06 国网江苏省电力公司常州供电公司 流体发电装置
CN107657080A (zh) * 2017-08-30 2018-02-02 中国科学院深海科学与工程研究所 一种基于自由曲面的深海照明灯设计方法及设计系统
CN107725477A (zh) * 2017-10-10 2018-02-23 北京航空航天大学 一种优化吸力面波系抑制风扇激波噪声的前缘设计方法
CN108916113A (zh) * 2018-06-13 2018-11-30 中国北方发动机研究所(天津) 一种直纹面压气机叶轮叶片曲面的调整方法
CN108916113B (zh) * 2018-06-13 2020-05-08 中国北方发动机研究所(天津) 一种直纹面压气机叶轮叶片曲面的调整方法
CN110727995B (zh) * 2018-06-29 2023-09-22 中国航发商用航空发动机有限责任公司 叶片形状的构造方法、叶片以及计算机设备
CN110727995A (zh) * 2018-06-29 2020-01-24 中国航发商用航空发动机有限责任公司 叶片形状的构造方法、叶片以及计算机设备
CN109571996B (zh) * 2018-12-04 2021-05-04 北京宇航系统工程研究所 样条曲线回转体正置正交复合材料网格结构硬模成型装置
CN109571996A (zh) * 2018-12-04 2019-04-05 北京宇航系统工程研究所 样条曲线回转体正置正交复合材料网格结构硬模成型装置
CN110593960A (zh) * 2019-09-19 2019-12-20 西安交通大学 一种可用于弯扭掠的轴流透平机械叶片参数化方法
CN110566503A (zh) * 2019-10-10 2019-12-13 珠海格力电器股份有限公司 离心风叶、离心风机、空调器和交通运输工具
CN110566503B (zh) * 2019-10-10 2020-05-15 珠海格力电器股份有限公司 离心风叶、离心风机、空调器和交通运输工具
CN111079239A (zh) * 2019-12-19 2020-04-28 中国航空发动机研究院 一种仿生压气机叶栅造型方法
CN111079239B (zh) * 2019-12-19 2023-07-21 中国航空发动机研究院 一种仿生压气机叶栅造型方法
CN111400834A (zh) * 2020-03-20 2020-07-10 国电联合动力技术有限公司 风力发电机组叶片翼型气动优化设计方法、模型及装置
CN111734678A (zh) * 2020-06-24 2020-10-02 西北工业大学 一种压气机叶型非对称前缘设计方法
CN113221264A (zh) * 2021-04-21 2021-08-06 山东双轮股份有限公司 海水淡化泵流道式导叶结构设计优化方法
CN113221264B (zh) * 2021-04-21 2022-09-27 山东双轮股份有限公司 海水淡化泵流道式导叶结构设计优化方法
CN113883093A (zh) * 2021-09-22 2022-01-04 哈尔滨工业大学 一种低反力度压气机叶片设计方法、动叶及压气机
CN117763878A (zh) * 2024-02-22 2024-03-26 清华大学 压气机叶型的确定方法、装置及存储介质
CN117763878B (zh) * 2024-02-22 2024-05-07 清华大学 压气机叶型的确定方法、装置及存储介质

Also Published As

Publication number Publication date
CN102996511B (zh) 2015-04-29

Similar Documents

Publication Publication Date Title
CN102996511B (zh) 曲率连续的可控扩散叶型的cad辅助设计方法
CN108153998B (zh) 离心鼓风机叶轮的全三维气动优化设计方法
EP3613995B1 (en) Method for designing blade of axial-flow fluid machine and blade
CN112347579B (zh) 一种压气机叶型设计方法及压气机叶型
CN106640210A (zh) 一种径‑轴流式透平膨胀机叶轮叶片三维型线的设计方法
CN110727995B (zh) 叶片形状的构造方法、叶片以及计算机设备
CN109145506B (zh) 一种高气动性能低噪声水平风力机外侧翼型的优化设计方法
CN105404743B (zh) B样条与曲率光滑连续性结合的风力机翼型设计方法
CN110059414A (zh) 一种直接控制通道的二维叶片造型方法
CN112576546B (zh) 一种非等厚翼型轴流叶片的优化方法
CN102587998B (zh) 一种用于控制气流分离的叶片吸力面凹槽的设计方法
CN101290636B (zh) 高速列车车头优化设计方法
CN114722518A (zh) 一种涡轮基本叶型参数化设计方法
CN104863799A (zh) 一种利用贝塞尔函数曲线的风力机翼型设计方法
CN109918778B (zh) 一种霜冰条件下风力机钝尾缘翼型优化设计方法
CN114186513A (zh) 一种具有反s型前缘的轴流压气机叶片造型设计方法
CN112685852B (zh) 一种轴流压气机负荷定制化气动优化方法
CN114036668A (zh) 一种基于离心叶轮中间面曲线的可变自由曲面叶片生成方法
CN112373673B (zh) 改善双凸起翼段性能的前缘双凸起结构的流动控制方法
CN109779971B (zh) 基于曲率控制的高负荷压气机叶型径向积叠造型优化方法
CN110245408B (zh) 一种汽轮机单圆弧压力面叶型设计方法
CN117329147A (zh) 一种叶轮机一体化结构及其设计方法
CN113569498A (zh) 一种轴流压气机端部弯曲静叶片设计方法
Sun et al. Optimization design of IGV profile in centrifugal compressor
Senoo et al. Three-dimensional design method for long blades of steam turbines using fourth-degree nurbs surface

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