CN102996511B - 曲率连续的可控扩散叶型的cad辅助设计方法 - Google Patents
曲率连续的可控扩散叶型的cad辅助设计方法 Download PDFInfo
- Publication number
- CN102996511B CN102996511B CN201210539482.9A CN201210539482A CN102996511B CN 102996511 B CN102996511 B CN 102996511B CN 201210539482 A CN201210539482 A CN 201210539482A CN 102996511 B CN102996511 B CN 102996511B
- Authority
- CN
- China
- Prior art keywords
- point
- suction surface
- curvature
- 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.)
- Active
Links
Landscapes
- Moulds For Moulding Plastics Or The Like (AREA)
Abstract
曲率连续的可控扩散型叶型的CAD辅助设计方法属于压气机叶片设计领域,其特征在于,含有:分别利用压力面、吸力面与前缘小圆的切点求前缘点,用三次样条拟合,求出吸力面、压力面上各型值点的曲率,用三次贝塞尔曲线求出曲率连续的前缘的吸力面侧型线和压力面侧型线,用B样条曲线给定吸力面的曲率分布,求取B样条曲线的控制点,使曲率分布对应的曲率连续吸力面型线与原始吸力面型线偏差最小,求取优化后叶型尾缘小圆的位置。本发明与原始可控扩散叶型相比,冲角范围扩大2.5°,且在低损失区域内的总压损失降低。
Description
技术领域
本发明涉及压气机叶片的优化设计方法。
背景技术
压气机是航空发动机和燃气轮机的重要部件之一,叶型设计是压气机设计的关键环节。早期的压气机叶型采用系列叶型(如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,坐标分别记为 ...,其中P1 ps为所述压力面与前缘小圆的切点,为压力面与尾缘小圆的切点,
吸力面上均匀分布的N个型值点,记为P1 ss,P2 ss...,PN ss坐标分别为 ...,N=500,其中P1 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_LE+θps_LE)/2,其中:
θss_LE为吸力面与前缘的切点P1 ss在前缘小圆上的圆心角, 其中:
θps_LE为压力面与前缘的切点P1 ps在前缘小圆上的圆心角,
步骤(2.3),对所述吸力面和压力面的叶型坐标使用三弯矩法进行样条拟合,求出所述吸力面和压力面各型值点Pi ss、Pi ps处的型线斜率,用和(i=1,2,...,N)表示,
步骤(2.4),按下式求出所述吸力面各型值点Pi ss处的流向坐标i=1,2,...,N和中间型值点的曲率i=2,...,N-1,
同理,求出所述压力面各型值点处的流向坐标i=1,2,...,N以及曲率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),选择吸力面上的某一型值点作为控制点P2的初始值,的横坐标xj的取值满足1≤j≤N,令点的坐标为(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的方向取,坐标为计算中间型值点Pi LE_ss的曲率i=2,...,49,并计算与控制点P2处吸力面型线曲率C2的偏差
步骤(3.2.4),判断步骤(3.2.3)得到的error_C的值:
若error_C<0.1,则优化后的前缘吸力面侧型线满足要求,
若error_C≥0.1,则把控制点P′2选为吸力面上的下一个型值点
步骤(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):
(x′2,y′2)为优化后的曲率连续的前缘与吸力面的切点P′2的坐标,
步骤(5.5),使用步骤(5.4)得到的各离散点的坐标进行三次样条插值,得到与原始吸力面型线的有效型值点的横坐标对应的曲率连续的吸力面型线的纵坐标从而得到优化后的曲率连续的吸力面型线的型值点的坐标,同时求出曲率分布确定的曲率连续的吸力面型线与原始吸力面型线的偏差 为与优化后的控制点P′2对应的离散点的序号,
步骤(5.6),使用优化软件iSIGHT获得最优的曲率分布控制点:
设定:曲率分布曲线的控制点为CP1~CP7,坐标分别为:(sCP1,CCP1),(sCP2,CCP2),...,(sCP7,CCP7),其中:
为吸力面上对应于控制点P′2的序号为的离散点处的流向坐标
为吸力面上对应于控制点P′2的序号为的离散点处的曲率,
分别为原始吸力面型线倒数第二个型值点的流向坐标和曲率,
设定:优化变量为sCP2~sCP6,CCP2~CCP6,
设定:变量的约束为sCP1≤sCP2,sCP3,sCP4,sCP5,sCP6≤sCP7,以及-10≤CCP2,CCP3,CCP4,CCP5,CCP6≤10,
设定:目标函数为曲率连续的吸力面型线与原始吸力面型线的偏差error_ss最小,
设定:优化算法为序列二次规划方法NLPQL,
得到:error_ss最小的坐标为...,所对应的型值点...,
步骤(6),根据曲率连续的吸力面型线和原始压力面型线确定尾缘小圆以及尾缘与吸力面和压力面分别相切的切点:
步骤(6.1),确定下述三条直线:
a,曲率连续的吸力面型线的型值点和组成的直线,
b,原始压力面型线的型值点和组成的直线,
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样条表达的优化后曲率分布的比较:----·为原始曲率分布,——为优化后曲率分布,为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,...坐标分别记为 ...,其中P1 ps为所述压力面与前缘小圆的切点,为压力面与尾缘小圆的切点,
吸力面上均匀分布的500个型值点,记为P1 ss,...,坐标分别为 ...,其中P1 ss为所述吸力面与前缘小圆的切点,为所述吸力面与尾缘小圆的切点;
步骤(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_LE+θps_LE)/2,其中:
θss_LE为吸力面与前缘的切点P1 ss在前缘小圆上的圆心角, 其中:
θps_LE为压力面与前缘的切点P1 ps在前缘小圆上的圆心角,
步骤(2.3),对所述吸力面和压力面的叶型坐标使用三弯矩法进行样条拟合,求出所述吸力面和压力面各型值点Pi ss、Pi ps处的型线斜率,用和(i=1,2,...,500)表示,
步骤(2.4),按下式求出所述吸力面各型值点Pi ss处的流向坐标i=1,2,...,500和中间型值点的曲率i=2,...,499,
同理,求出所述压力面各型值点处的流向坐标i=1,2,...,500,以及曲率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,令点的坐标为(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的方向取,坐标为计算中间型值点Pi LE_ss的曲率i=2,...,49,并计算与控制点P2处吸力面型线曲率C2的偏差
步骤(3.2.4),判断步骤(3.2.3)得到的error_C的值:
若error_C<0.1,则优化后的前缘吸力面侧型线满足要求,
若error_C≥0.1,则把控制点P′2选为吸力面上的下一个型值点
步骤(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为新的前缘与吸力面的切点,并且满足在吸力面与前缘的切点处曲率的连续性,在本实施例中,为满足要求的控制点,则控制点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),按下式求出各离散点对应的曲线倾斜角
s′i为B样条曲线上的序号为i的离散点的流向坐标,k′2为优化后的曲率连续的前缘与吸力面的切点P′2的斜率,
步骤(5.4),按下式求出各离散点的横坐标和纵坐标(X′i,Y′i):
(x′2,y′2)为优化后的曲率连续的前缘与吸力面的切点P′2的坐标,
步骤(5.5),使用步骤(5.4)得到的各离散点的坐标进行三次样条插值,得到与原始吸力面型线的有效型值点的横坐标对应的曲率连续的吸力面型线的纵坐标从而得到优化后的曲率连续的吸力面型线的型值点的坐标,同时求出曲率分布确定的曲率连续的吸力面型线与原始吸力面型线的偏差 15≤i≤500,
步骤(5.6),使用优化软件iSIGHT获得最优的曲率分布控制点:
设定:曲率分布曲线的控制点为CP1~CP7,坐标分别为:(sCP1,CCP1),(sCP2,CCP2),...,(sCP7,CCP7),其中:
为吸力面上对应于控制点P′2的离散点处的流向坐标,
为吸力面上对应于控制点P′2的离散点处的曲率,
分别为原始吸力面型线倒数第二个型值点的流向坐标和曲率,
设定:优化变量为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最小的坐标为 ...,所对应的型值点 ...,
步骤(6),如图9,根据曲率连续的吸力面型线和原始压力面型线确定尾缘小圆以及尾缘与吸力面和压力面分别相切的切点:
步骤(6.1),确定下述三条直线:
a,曲率连续的吸力面型线的型值点和组成的直线,
b,原始压力面型线的型值点和组成的直线,
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,坐标分别记为 ...,其中P1 ps为所述压力面与前缘小圆的切点,为压力面与尾缘小圆的切点,
吸力面上均匀分布的N个型值点,记为P1 ss,P2 ss...,PN ss 坐标分别为 ...,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_LE+θps_LE)/2,其中:
θss_LE为吸力面与前缘的切点P1 ss在前缘小圆上的圆心角, 其中:
θps_LE为压力面与前缘的切点P1 ps在前缘小圆上的圆心角,
步骤(2.3),对所述吸力面和压力面的叶型坐标使用三弯矩法进行样条拟合,求出所述吸力面和压力面各型值点Pi ss、Pi ps处的型线斜率,用和(i=1,2,...,N)表示,
步骤(2.4),按下式求出所述吸力面各型值点Pi ss处的流向坐标i=1,2,...,N和中间型值点的曲率i=2,...,N-1,
同理,求出所述压力面各型值点处的流向坐标i=1,2,...,N以及曲率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),选择吸力面上的某一型值点作为控制点P2的初始值,的横坐标xj的取值满足1≤j≤N,令点的坐标为(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的曲率i=2,...,49,并计算与控制点P2处吸力面型线曲率C2的偏差
步骤(3.2.4),判断步骤(3.2.3)得到的error_C的值:
若error_C<0.1,则优化后的前缘吸力面侧型线满足要求,
若error_1≥0.1,则把控制点P′2选为吸力面上的下一个型值点
步骤(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):
(x′2,y′2)为优化后的曲率连续的前缘与吸力面的切点P′2的坐标,
步骤(5.5),使用步骤(5.4)得到的各离散点的坐标进行三次样条插值,得到与原始吸力面型线的有效型值点的横坐标对应的曲率连续的吸力面型线的纵坐标从而得到优化后的曲率连续的吸力面型线的型值点的坐标,同时求出曲率分布确定的曲率连续的吸力面型线与原始吸力面型线的偏差 为与优化后的控制点P′2对应的离散点的序号,
步骤(5.6),使用优化软件iSIGHT获得最优的曲率分布控制点:
设定:曲率分布曲线的控制点为CP1~CP7,坐标分别为:(sCP1,CCP1),(sCP2,CCP2),...,(sCP7,CCP7),其中:
为吸力面上对应于控制点P′2的序号为的离散点处的流向坐标
为吸力面上对应于控制点P′2的序号为的离散点处的曲率,
分别为原始吸力面型线倒数第二个型值点的流向坐标和曲率,
设定:优化变量为sCP2~sCP6,CCP2~CCP6,
设定:变量的约束为sCP1≤sCP2,sCP3,sCP4,sCP5,sCP6≤sCP7,以及-10≤CCP2,CCP3,CCP4,CCP5,CCP6≤10,
设定:目标函数为曲率连续的吸力面型线与原始吸力面型线的偏差error_ss最小,
设定:优化算法为序列二次规划方法NLPQL,
得到:error_ss最小的坐标为...,所对应的型值点...,
步骤(6),根据曲率连续的吸力面型线和原始压力面型线确定尾缘小圆以及尾缘与吸力面和压力面分别相切的切点:
步骤(6.1),确定下述三条直线:
a,曲率连续的吸力面型线的型值点和组成的直线,
b,原始压力面型线的型值点和组成的直线,
c,原始尾缘小圆的切线,切点为原始前缘小圆圆心与原始尾缘小圆圆心的连线的延长线与原始尾缘小圆的交点,
步骤(6.2),确定同时与步骤(6.1)中的三条直线a,b,c相切的圆即为尾缘小圆,
步骤(6.3),根据步骤(6.2)的结果计算出尾缘小圆的圆心坐标(xTE′,yTE′)和半径rTE′,
步骤(6.4),确定:
尾缘小圆与直线a的交点为尾缘与吸力面的切点,
尾缘小圆与直线b的交点为尾缘与压力面的切点。
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 CN102996511A (zh) | 2013-03-27 |
CN102996511B true 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 (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107178526A (zh) * | 2017-07-11 | 2017-09-19 | 上海电气燃气轮机有限公司 | 轴流压气机叶型设计方法及轴流压气机叶型 |
EP3613995A4 (en) * | 2017-04-17 | 2020-12-30 | Ihi Corporation | AXIAL FLUID MACHINE BLADE DESIGN PROCESS, AND BLADE |
Families Citing this family (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103398021A (zh) * | 2013-07-20 | 2013-11-20 | 大连理工大学 | 一种低转速高压比贝塞尔中弧线强扭曲叶片轴流叶轮 |
CN103527514A (zh) * | 2013-10-29 | 2014-01-22 | 中国北车集团大连机车研究所有限公司 | 机车用离心通风机叶轮 |
CN103541774B (zh) * | 2013-11-14 | 2015-06-17 | 上海汽轮机厂有限公司 | 涡轮叶片设计方法 |
CN104613157B (zh) * | 2015-01-26 | 2017-05-10 | 北京理工大学 | 基于贝塞尔曲线的液力变矩器二维叶片型线构造方法 |
CN106021681B (zh) * | 2016-05-13 | 2019-06-25 | 北京航空航天大学 | 一种压气机叶型优化方法及装置 |
CN106021782A (zh) * | 2016-05-31 | 2016-10-12 | 西北工业大学 | 基于中弧线的叶片前后缘拟合及截面线光滑重构方法 |
CN106151105B (zh) * | 2016-08-03 | 2018-07-06 | 华中科技大学 | 一种轴流叶型 |
CN106351878B (zh) * | 2016-09-28 | 2018-11-02 | 华中科技大学 | 一种轴流掠形叶片 |
CN109083798B (zh) * | 2017-06-13 | 2024-02-06 | 国网江苏省电力公司常州供电公司 | 流体发电装置 |
CN107657080A (zh) * | 2017-08-30 | 2018-02-02 | 中国科学院深海科学与工程研究所 | 一种基于自由曲面的深海照明灯设计方法及设计系统 |
CN107725477B (zh) * | 2017-10-10 | 2019-05-07 | 北京航空航天大学 | 一种优化吸力面波系抑制风扇激波噪声的前缘设计方法 |
CN108916113B (zh) * | 2018-06-13 | 2020-05-08 | 中国北方发动机研究所(天津) | 一种直纹面压气机叶轮叶片曲面的调整方法 |
CN110727995B (zh) * | 2018-06-29 | 2023-09-22 | 中国航发商用航空发动机有限责任公司 | 叶片形状的构造方法、叶片以及计算机设备 |
CN109571996B (zh) * | 2018-12-04 | 2021-05-04 | 北京宇航系统工程研究所 | 样条曲线回转体正置正交复合材料网格结构硬模成型装置 |
CN110593960B (zh) * | 2019-09-19 | 2021-03-23 | 西安交通大学 | 一种可用于弯扭掠的轴流透平机械叶片参数化方法 |
CN110566503B (zh) * | 2019-10-10 | 2020-05-15 | 珠海格力电器股份有限公司 | 离心风叶、离心风机、空调器和交通运输工具 |
CN111079239B (zh) * | 2019-12-19 | 2023-07-21 | 中国航空发动机研究院 | 一种仿生压气机叶栅造型方法 |
CN111400834A (zh) * | 2020-03-20 | 2020-07-10 | 国电联合动力技术有限公司 | 风力发电机组叶片翼型气动优化设计方法、模型及装置 |
CN111734678B (zh) * | 2020-06-24 | 2021-08-31 | 西北工业大学 | 一种压气机叶型非对称前缘设计方法 |
CN113221264B (zh) * | 2021-04-21 | 2022-09-27 | 山东双轮股份有限公司 | 海水淡化泵流道式导叶结构设计优化方法 |
CN113883093B (zh) * | 2021-09-22 | 2023-07-04 | 哈尔滨工业大学 | 一种低反力度压气机叶片设计方法、动叶及压气机 |
CN117763878B (zh) * | 2024-02-22 | 2024-05-07 | 清华大学 | 压气机叶型的确定方法、装置及存储介质 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201159202Y (zh) * | 2008-02-28 | 2008-12-03 | 大连海事大学 | 一种压气机动叶叶尖小翼 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5351637B2 (ja) * | 2009-07-09 | 2013-11-27 | 三菱重工業株式会社 | 翼体および回転機械 |
-
2012
- 2012-12-13 CN CN201210539482.9A patent/CN102996511B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201159202Y (zh) * | 2008-02-28 | 2008-12-03 | 大连海事大学 | 一种压气机动叶叶尖小翼 |
Non-Patent Citations (2)
Title |
---|
某离心叶轮的CFD分析和结果确认;李想等;《工程热物理学报》;20080831;第29卷(第8期);1292-1296 * |
马文生等.叶顶间隙对压气机性能的影响.《动力工程》.2007,第27卷(第6期),863-867. * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3613995A4 (en) * | 2017-04-17 | 2020-12-30 | Ihi Corporation | AXIAL FLUID MACHINE BLADE DESIGN PROCESS, AND BLADE |
CN107178526A (zh) * | 2017-07-11 | 2017-09-19 | 上海电气燃气轮机有限公司 | 轴流压气机叶型设计方法及轴流压气机叶型 |
CN107178526B (zh) * | 2017-07-11 | 2019-09-24 | 上海电气燃气轮机有限公司 | 轴流压气机叶型设计方法及轴流压气机叶型 |
Also Published As
Publication number | Publication date |
---|---|
CN102996511A (zh) | 2013-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102996511B (zh) | 曲率连续的可控扩散叶型的cad辅助设计方法 | |
CN108153998B (zh) | 离心鼓风机叶轮的全三维气动优化设计方法 | |
CN109779972B (zh) | 一种兼顾低雷诺数工况性能的压气机二维叶型优化方法 | |
CN106640210A (zh) | 一种径‑轴流式透平膨胀机叶轮叶片三维型线的设计方法 | |
Voß et al. | Automated multiobjective optimisation in axial compressor blade design | |
US11333167B2 (en) | Method of designing blade of axial flow fluid machine and blade | |
CN110727995B (zh) | 叶片形状的构造方法、叶片以及计算机设备 | |
CN112347579B (zh) | 一种压气机叶型设计方法及压气机叶型 | |
CN106089801B (zh) | 一种压气机叶片造型方法 | |
CN112576546B (zh) | 一种非等厚翼型轴流叶片的优化方法 | |
CN103425831A (zh) | 基于多弧段曲线的压力载荷下的结构拓扑-形状联合优化方法 | |
CN102587998B (zh) | 一种用于控制气流分离的叶片吸力面凹槽的设计方法 | |
CN101290636B (zh) | 高速列车车头优化设计方法 | |
CN105404743A (zh) | B样条与曲率光滑连续性结合的风力机翼型设计方法 | |
CN115952604A (zh) | 一种基于自由曲面叶片中性面的离心叶轮叶片生成方法 | |
CN104863799A (zh) | 一种利用贝塞尔函数曲线的风力机翼型设计方法 | |
CN109918778B (zh) | 一种霜冰条件下风力机钝尾缘翼型优化设计方法 | |
CN104317992A (zh) | 风力机翼型正设计方法及风力机翼型族 | |
Joly et al. | Full design of a highly loaded fan by multi-objective optimization of through-flow and high-fidelity aero-mechanical performances | |
Cheng et al. | A phased aerodynamic optimization method for compressors based on multi-degrees-of-freedom surface parameterization | |
CN109779971B (zh) | 基于曲率控制的高负荷压气机叶型径向积叠造型优化方法 | |
CN110245408B (zh) | 一种汽轮机单圆弧压力面叶型设计方法 | |
CN101229849A (zh) | 一种进气道斜切进口及其设计方法 | |
CN103670528B (zh) | 透平叶片的加载方法 | |
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 |