CN107862104A - 一种环形喷嘴直线‑弧线叶片型线设计方法 - Google Patents
一种环形喷嘴直线‑弧线叶片型线设计方法 Download PDFInfo
- Publication number
- CN107862104A CN107862104A CN201710889452.3A CN201710889452A CN107862104A CN 107862104 A CN107862104 A CN 107862104A CN 201710889452 A CN201710889452 A CN 201710889452A CN 107862104 A CN107862104 A CN 107862104A
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtr
- mtd
- mtable
- 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
Classifications
-
- 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
-
- 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
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)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种用于径‑轴流式透平膨胀机环形喷嘴直线‑弧线叶片型线设计的方法,该方法包括方法将叶片分为吸力面和压力面两个部分,每个部分的叶片型线又分为直线段和弧线段,通过编程利用数学表达式对直线段和弧线段进行坐标求解。本发明采用Bezier曲线求解弧线段坐标以实现直线段和弧线段的光滑过渡,从而可以很好地保证喷嘴叶片型线的几何曲率连续性,避免更多的流动损失;二是通过编程将叶片型线分为直线段和弧线段并分别用数学式来表示,设计过程更加简便、快捷,方便后续对喷嘴叶片型线的加工制造、数值模拟和优化设计等工作过程。
Description
技术领域
本发明涉及一种环形喷嘴叶片型线设计方法,特别是涉及一种 用于径-轴流式透平膨胀机环形喷嘴直线-弧线叶片型线设计的方法。
背景技术
透平膨胀机是低温***中主要的冷量来源,也是现代空分装置 和低温液化***中的核心部件,其性能对整套装置的经济性与安全 性具有重要意义。透平膨胀机的主通流部分包括环形喷嘴、径-轴流 式工作轮和扩压器,其中环形喷嘴是透平膨胀机发生能量的主要流 道之一,环形喷嘴叶片型线的设计对于透平膨胀机性能具有重要影 响。直线-弧线叶片是环形喷嘴叶片型式的一种,其叶片型线采用简 单的直线和弧线组成,具有设计简单、加工方便等优点。传统方法 中弧线段采用圆弧,圆弧参数则根据经验参数选取。这样的几何造 型方法,主要有以下两方面的缺点:一是通过经验参数法设计弧线 段,很难保证直线段和弧线段的光滑过渡,从而不可避免地造成更 多的流动损失;二是通过经验作图法选取直线段和弧线段的比例, 修改型线过程繁琐复杂,不利于后续对叶片型线的加工制造、数值 模拟和优化设计等工作过程。
发明内容
针对背景技术中的问题,本发明的主要目的在于提供一种用于 径-轴流式透平膨胀机环形喷嘴直线-弧线叶片型线设计的方法,
为实现本发明目的而采用的技术方案是这样的,一种环形喷嘴 直线-弧线叶片型线设计方法,其特征在于,包括以下步骤:
1〕读入环形喷嘴叶片基本参数,包括进口直径R1、出口直径R2、 出口叶片角α2、叶片数目NB、喉部宽度b及叶片高度h;
2〕以环形喷嘴叶片俯视面(从叶片上方俯视)圆心为原点,选 定二维x-y正交坐标系,
自定义参数下标:“1”表示与进口相关,“2”表示与出口相关, “h”表示与根部相关(hub)线,“s”表示与顶部相关,“PS”表示 与压力面相关,“SS”表示与吸力面相关;
自定义压力面直线段选取比例lPS、吸力面直线段选取比例lSS、 叶片前缘顶点选取比例lp及吸力面Bezier曲线插值点选取比例lc,SS等型线设计参数,以及直线段输出点个数NL和弧线段输出点个数NC;
3〕根据基本参数和设计参数,求解单元流道中心线Lmid方程:
再通过直线平移法则计算得到压力面基准线LPS和吸力面基准线LSS:
4〕求解压力面基准线LPS与x和y轴的交点坐标分别为:
基于基准线和叶片数目,对压力面基准线LPS进行阵列,得到与吸力 面基准线LSS相交的压力面基准线LPS,a:
并求解两条直线交点A坐标:
以及两条直线与喷嘴外径圆的交点B和C的坐标:
5〕依据压力面和吸力面直线段选取比例lPS和lSS,分别计算压 力面和吸力面基准线直线段和弧线段分割点M和N的坐标:
6〕依据叶片前缘顶点选取比例lp,计算压力面叶片型线和吸力 面型线前缘交点D坐标,首先计算辅助线LPSa,SS方程:
再联立求解辅助线LPSa,SS和辅助外径圆O1的方程得前缘交点D坐标:
7〕依据分割点M和N及前缘交点D,基于Bezier曲线方程,求 解压力面根部型线和吸力面根部型线的曲线段,具体步骤分为:
7-1〕分别求解Bezier曲线的插值点坐标
7-1-1〕求解吸力面叶片型线Bezier曲线段插值点P1、P2、P3和P4坐标,P1选为N点,P4选为D点,即:
P2和P3则通过辅助切线Lt与直线段LSS联立求解求得。辅助切线Lt方程为:
则辅助线Lt与吸力面基准线LSS的交点E的坐标:
根据E点坐标及吸力面Bezier曲线插值点选取比例lc,SS(理论值范 围为0~1)求解得插值点P2和P3的坐标:
7-1-2〕求解压力面叶片型线Bezier曲线段插值点Q1、Q2和Q3坐标,Q1选为M点,Q3选为D点,即:
Q2则选取为辅助线Lt与压力面基准线LPS,a的交点:
7-2〕基于Bezier曲线方程,依据Bezier曲线段插值点坐标和 Bezier曲线段输出点个数NC分别求解得到压力面和吸力面弧线段 输出点坐标:
7-2-1〕求解吸力面弧线段输出点坐标
构建向量UXSS=[xP1,xP2,xP3,xP4]和向量UYSS=[yP1,yP2,yP3,yP4],则第 I个点,采用三次方Bezier公式,其对应的计算向量BI:
t=(I-1)/(NC-1)
I=1,...,NL,进而可以求得第I个点坐标为:
7-2-2〕求解压力面弧线段输出点坐标
构建向量UXPS=[xQ1,xQ2,xQ3]和向量UYPS=[yQ1,yQ2,yQ3],采用二次方 Bezier公式,则第J个点,采用二次方Bezier公式,其对应的计算 向量BJ:
t=(J-1)/(NC-1)
J=1,...,NL,进而可以求得第J个点坐标为:
8〕依据吸力面和压力面基准线方程和直线段输出点个数NL,求 解直线段输出点坐标:
对于吸力面叶片型线直线段第K个点坐标为:
K=1,...,NC
对于压力面叶片型线直线段第L个点坐标:
L=1,...,NC
9〕结合吸力面和压力面直线段输出点坐标和Bezier曲线段输 出点坐标,吸力面和压力面叶片型线输出点坐标具体可表示为:
吸力面:
压力面:
10〕建立三维坐标系,以步骤8〕中构造型线所在平面为x-y 平面,以垂直于x-y平面的直线并过圆心O的直线为z轴,则步骤8 中构造的叶片型线为根部型线,其三维点坐标可以表示为(xh,yh,0), 由输入参数叶片高度h,进而可以计算得到与根部型线对应的顶部型 线为:
从而完成吸力面和压力面三维叶片的型线设计。
值得说明的是,本方法将叶片分为吸力面和压力面两个部分, 每个部分的叶片型线又分为直线段和弧线段,通过编程利用数学表 达式对直线段和弧线段进行坐标求解。本发明采用Bezier曲线求解 弧线段坐标以实现直线段和弧线段的光滑过渡,从而可以很好地保 证喷嘴叶片型线的几何曲率连续性,避免更多的流动损失;二是通 过编程将叶片型线分为直线段和弧线段并分别用数学式来表示,设 计过程更加简便、快捷,方便后续对喷嘴叶片型线的加工制造、数 值模拟和优化设计等工作过程。
附图说明
图1为本发明环形喷嘴直线-弧线叶片型线设计方法的流程示意 图;
图2为本发明环形喷嘴叶片型线设计二维坐标系示意意图;
图3为本发明环形喷嘴叶片型线设计过程示意图;
图4为实施例中计算建立的环形喷嘴直线-弧线叶片三维整体结 构图;
图5为三种不同吸力面直线段选取比例lss计算所得的喷嘴叶片 型线图。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本 发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思 想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换 和变更,均应包括在本发明的保护范围内。
本实施例以某低压空分流程配套的透平膨胀机为实施例,透平 膨胀机设计参数及叶轮主要结构尺寸如表1,且为轴向出气。
表1透平膨胀机设计参数及叶轮主要结构尺寸
按以下步骤进行设计(参数中,用“1”表示进口,“2”表示出 口,“h”表示hub曲线,即轮盘基线,“s”表示shroud曲线,即轮 盖基线或叶顶线):
1〕读入环形喷嘴基本参数,包括进口直径R1=39mm、出口直径 R2=29mm、出口叶片角α2=16°、叶片数目NB=25、喉部宽度b=1.8mm 及叶片高度h=2.7mm;
2〕以环形喷嘴圆心为原点,选定二维x-y坐标系,自定义压力 面直线段选取比例lPS=0.95、吸力面直线段选取比例lSS=0.5、叶片 前缘顶点选取比例lp=0.9及吸力面Bezier曲线插值点选取比例 lc,SS=0.6等型线设计参数,以及直线段输出点个数NL=51和弧线段输 出点个数NC=101;
3〕根据基本参数和设计参数,求解单元流道中心线Lmid方程:
再通过直线平移法则计算得到压力面基准线LPS和吸力面基准线 LSS:
4〕求解压力面基准线LPS与x和y轴的交点坐标分别为:
基于基准线和叶片数目,对压力面基准线LPS进行阵列,得到与 吸力面基准线LSS相交的压力面基准线LPS,a:
并求解两条直线交点A坐标:
以及两条直线与喷嘴外径圆的交点B和C的坐标:
5〕依据压力面和吸力面直线段选取比例lPS和lSS,分别计算压 力面和吸力面基准线直线段和弧线段分割点M和N的坐标:
6〕依据叶片前缘顶点选取比例lp,计算压力面叶片型线和吸力 面型线前缘交点D坐标,首先计算辅助线LPSa,SS方程:
再联立求解辅助线LPS,SS和辅助外径圆O1的方程得前缘交点D坐 标:
7〕依据分割点M和N及前缘交点D,基于Bezier曲线方程,求 解压力面根部型线和吸力面根部型线的曲线段,具体步骤分为:
7-1〕分别求解Bezier曲线的插值点坐标
7-1-1〕求解吸力面叶片型线Bezier曲线段插值点P1、P2、P3和P4坐标,P1选为N点,P4选为D点,即:
P2和P3则通过辅助切线Lt与直线段LSS联立求解求得。辅助切线 Lt方程为:
则辅助线Lt与吸力面基准线LSS的交点E的坐标:
根据E点坐标及吸力面Bezier曲线插值点选取比例lc,SS求解得 插值点P2和P3的坐标:
7-1-2〕求解压力面叶片型线Bezier曲线段插值点Q1、Q2和Q3坐标,Q1选为M点,Q3选为D点,即:
Q2则选取为辅助线Lt与压力面基准线LPS,a的交点:
7-2〕基于Bezier曲线方程,依据Bezier曲线段插值点坐标和 Bezier曲线段输出点个数NC分别求解得到压力面和吸力面弧线段 输出点坐标:
7-2-1〕求解吸力面弧线段输出点坐标
构建向量UXSS=[31.786,32.861,35.680,37.081]和向量 UYSS=[-6.451,-10.199,-11.533,-10.758],则第I个点,采用三次 方Bezier公式,其对应的计算向量BI:
t=(I-1)/(NC-1)
进而可以求得第I个点坐标为:
以第51个点为例,求得计算向量B51为:
t=(51-1)/(101-1)=0.5
求得第51个点坐标为:
7-2-2〕求解压力面弧线段输出点坐标
构建向量UXPS=[37.077,37.464,37.081]和向量 UYPS=[-9.886,-10.546,-10.758],采用二次方Bezier公式,则第J 个点,采用二次方Bezier公式,其对应的计算向量BJ:
t=(J-1)/(NC-1)
进而可以求得第J个点坐标为:
以第51个点为例,求得计算向量B51为:
t=(51-1)/(101-1)=0.5
求得第51个点坐标为:
8〕依据吸力面和压力面基准线方程和直线段输出点个数NL,求 解直线段输出点坐标。对于吸力面叶片型线直线段第K个点坐标为:
以第21个点为例,求得坐标为:
对于压力面叶片型线直线段第L个点坐标:
以第31个点为例,求得坐标为:
9〕结合吸力面和压力面直线段输出点坐标和Bezier曲线段输 出点坐标,吸力面和压力面叶片型线输出点坐标具体可表示为:
吸力面:
压力面:
10〕如图3所示,建立三维坐标系,以步骤8〕中构造型线所在 平面为x-y平面,以垂直于x-y平面的直线并过圆心O的直线为z 轴,则步骤8中构造的叶片型线为根部型线,其三维点坐标可以表 示为(xh,yh,0),由输入参数叶片高度h,进而可以计算得到与根部型 线对应的顶部型线为:
从而完成吸力面和压力面三维叶片的型线设计。
如图4所示为通过上述计算建立的环形喷嘴直线-弧线叶片三维 整体结构图。当修改步骤1中的结构参数或者步骤2中的设计参数 时,重复步骤3-步骤10可以计算得到新的叶片结构,从而方便快捷 地得到一系列具有不同参数的喷嘴叶片结构,如图5为三种不同吸 力面直线段选取比例lSS计算所得的喷嘴叶片型线图,从而有利于后 续工作中研究不同参数对膨胀流动和透平膨胀机性能的影响。
Claims (1)
1.一种环形喷嘴直线-弧线叶片型线设计方法,其特征在于,包括以下步骤:
1〕读入所述环形喷嘴叶片基本参数,包括所述进口直径R1、出口直径R2、出口叶片角α2、叶片数目NB、喉部宽度b及叶片高度h。
2〕以环形喷嘴叶片俯视面圆心为原点,选定二维x-y正交坐标系,
自定义参数下标:“1”表示与进口相关,“2”表示与出口相关,“h”表示与根部相关(hub)线,“s”表示与顶部相关,“PS”表示与压力面相关,“SS”表示与吸力面相关;
自定义压力面直线段选取比例lPS、吸力面直线段选取比例lSS、叶片前缘顶点选取比例lp及吸力面Bezier曲线插值点选取比例lc,SS等型线设计参数,以及直线段输出点个数NL和弧线段输出点个数NC;
3〕根据基本参数和设计参数,求解单元流道中心线Lmid方程:
再通过直线平移法则计算得到压力面基准线LPS和吸力面基准线LSS:
4〕求解压力面基准线LPS与x和y轴的交点坐标分别为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>d</mi>
<mi>x</mi>
<mo>;</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<mn>0</mn>
<mo>;</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
基于基准线和叶片数目,对压力面基准线LPS进行阵列,得到与吸力面基准线LSS相交的压力面基准线LPS,a:
并求解两条直线交点A坐标:
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>A</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>A</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>:</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>:</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
以及两条直线与喷嘴外径圆的交点B和C的坐标:
<mrow>
<mi>B</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>B</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>B</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>:</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>O</mi>
<mo>:</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>C</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>C</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>C</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>:</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>O</mi>
<mo>:</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<msubsup>
<mi>R</mi>
<mn>1</mn>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
5〕依据压力面和吸力面直线段选取比例lPS和lSS,分别计算压力面和吸力面基准线直线段和弧线段分割点M和N的坐标:
<mrow>
<mi>M</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>M</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>M</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>B</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>B</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>N</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>C</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>C</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
6〕依据叶片前缘顶点选取比例lp,计算压力面叶片型线和吸力面型线前缘交点D坐标,首先计算辅助线LPSa,SS方程:
再联立求解辅助线LPSa,SS和辅助外径圆O1的方程得前缘交点D坐标:
<mrow>
<mi>D</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
<mi>a</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>:</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>A</mi>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>d</mi>
<mo>,</mo>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>A</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>O</mi>
<mn>1</mn>
</msub>
<mo>:</mo>
<msup>
<mi>x</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>y</mi>
<mn>2</mn>
</msup>
<mo>=</mo>
<msup>
<mrow>
<mo>(</mo>
<mn>0.99</mn>
<msub>
<mi>R</mi>
<mn>1</mn>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
7〕依据分割点M和N及前缘交点D,基于Bezier曲线方程,求解压力面根部型线和吸力面根部型线的曲线段,具体步骤分为:
7-1〕分别求解Bezier曲线的插值点坐标。
7-1-1〕求解吸力面叶片型线Bezier曲线段插值点P1、P2、P3和P4坐标,P1选为N点,P4选为D点,即:
<mrow>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<msub>
<mi>P</mi>
<mn>1</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mi>P</mi>
<mn>4</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<msub>
<mi>P</mi>
<mn>4</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<msub>
<mi>P</mi>
<mn>4</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<msub>
<mi>P</mi>
<mn>4</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
P2和P3则通过辅助切线Lt与直线段LSS联立求解求得。辅助切线Lt方程为:
则辅助线Lt与吸力面基准线LSS的交点E的坐标:
根据E点坐标及吸力面Bezier曲线插值点选取比例lc,SS求解得插值点P2和P3的坐标:
<mrow>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mi>E</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<msub>
<mi>P</mi>
<mn>2</mn>
</msub>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>E</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>E</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<msub>
<mi>P</mi>
<mn>3</mn>
</msub>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>E</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
7-1-2〕求解压力面叶片型线Bezier曲线段插值点Q1、Q2和Q3坐标,Q1选为M点,Q3选为D点,即:
<mrow>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>M</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<msub>
<mi>Q</mi>
<mn>1</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<msub>
<mi>Q</mi>
<mn>3</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mi>Q</mi>
<mn>3</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<msub>
<mi>Q</mi>
<mn>3</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<msub>
<mi>Q</mi>
<mn>3</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<msub>
<mi>Q</mi>
<mn>3</mn>
</msub>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Q2则选取为辅助线Lt与压力面基准线LPS,a的交点:
<mrow>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<msub>
<mi>Q</mi>
<mn>2</mn>
</msub>
</msub>
<mo>)</mo>
</mrow>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mi>t</mi>
</msub>
<mo>:</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>D</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>L</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>:</mo>
<mi>y</mi>
<mo>-</mo>
<msub>
<mi>y</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>P</mi>
<mi>S</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
7-2〕基于Bezier曲线方程,依据Bezier曲线段插值点坐标和Bezier曲线段输出点个数NC分别求解得到压力面和吸力面弧线段输出点坐标:
7-2-1〕求解吸力面弧线段输出点坐标
构建向量UXSS=[xP1,xP2,xP3,xP4]和向量UYSS=[yP1,yP2,yP3,yP4],则第I个点,采用三次方Bezier公式,其对应的计算向量BI:
<mrow>
<msub>
<mi>B</mi>
<mi>I</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>3</mn>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>3</mn>
<mi>t</mi>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>3</mn>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>t</mi>
<mn>3</mn>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
t=(I-1)/(NC-1)
I=1,...,NL,进而可以求得第I个点坐标为:
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>I</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>I</mi>
</msub>
<mo>)</mo>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>I</mi>
</msub>
<mo>=</mo>
<msub>
<mi>UX</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>B</mi>
<mi>I</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>I</mi>
</msub>
<mo>=</mo>
<msub>
<mi>UY</mi>
<mrow>
<mi>S</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>B</mi>
<mi>I</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
7-2-2〕求解压力面弧线段输出点坐标
构建向量UXPS=[xQ1,xQ2,xQ3]和向量UYPS=[yQ1,yQ2,yQ3],采用二次方Bezier公式,则第J个点,采用二次方Bezier公式,其对应的计算向量BJ:
<mrow>
<msub>
<mi>B</mi>
<mi>J</mi>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<mi>t</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>t</mi>
<mn>2</mn>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
t=(J-1)/(NC-1)
J=1,...,NL,进而可以求得第J个点坐标为:
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>J</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>J</mi>
</msub>
<mo>)</mo>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>J</mi>
</msub>
<mo>=</mo>
<msub>
<mi>UX</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>B</mi>
<mi>J</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>J</mi>
</msub>
<mo>=</mo>
<msub>
<mi>UY</mi>
<mrow>
<mi>P</mi>
<mi>S</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>B</mi>
<mi>J</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
8〕依据吸力面和压力面基准线方程和直线段输出点个数NL,求解直线段输出点坐标:
对于吸力面叶片型线直线段第K个点坐标为:
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>K</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>K</mi>
</msub>
<mo>)</mo>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>K</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>K</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mi>t</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>K</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mi>t</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>N</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
K=1,...,NC
对于压力面叶片型线直线段第L个点坐标:
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>L</mi>
</msub>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>L</mi>
</msub>
<mo>)</mo>
<mo>:</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mi>L</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mi>t</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>x</mi>
<mi>M</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>A</mi>
</msub>
<mo>+</mo>
<mi>t</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>y</mi>
<mi>M</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
L=1,...,NC
9〕结合吸力面和压力面直线段输出点坐标和Bezier曲线段输出点坐标,吸力面和压力面叶片型线输出点坐标具体可表示为:
10〕建立三维坐标系,以步骤8〕中构造型线所在平面为x-y平面,以垂直于x-y平面的直线并过圆心O的直线为z轴,则步骤8中构造的叶片型线为根部型线,其三维点坐标可以表示为(xh,yh,0),由输入参数叶片高度h,进而可以计算得到与根部型线对应的顶部型线为:
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<msub>
<mi>x</mi>
<mi>h</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>y</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<msub>
<mi>y</mi>
<mi>h</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>z</mi>
<mi>s</mi>
</msub>
<mo>=</mo>
<mi>h</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
从而完成吸力面和压力面三维叶片的型线设计。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710889452.3A CN107862104B (zh) | 2017-09-27 | 2017-09-27 | 一种环形喷嘴直线-弧线叶片型线设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710889452.3A CN107862104B (zh) | 2017-09-27 | 2017-09-27 | 一种环形喷嘴直线-弧线叶片型线设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107862104A true CN107862104A (zh) | 2018-03-30 |
CN107862104B CN107862104B (zh) | 2020-09-29 |
Family
ID=61698731
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710889452.3A Active CN107862104B (zh) | 2017-09-27 | 2017-09-27 | 一种环形喷嘴直线-弧线叶片型线设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107862104B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114962322A (zh) * | 2022-05-24 | 2022-08-30 | 威灵(芜湖)电机制造有限公司 | 叶轮结构、风机及家用电器 |
CN115329506A (zh) * | 2022-10-14 | 2022-11-11 | 中国空气动力研究与发展中心空天技术研究所 | 用于自由射流试验的动力装置安装架的设计方法及安装架 |
CN117235937A (zh) * | 2023-11-14 | 2023-12-15 | 广东顺威精密塑料股份有限公司 | 一种多翼离心风机叶轮的周向弯曲叶片的设计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011243028A (ja) * | 2010-05-19 | 2011-12-01 | Ihi Corp | ターボ機械の翼形状設計方法およびプログラム |
CN104376144A (zh) * | 2014-08-06 | 2015-02-25 | 如皋透平叶片制造有限公司 | 一种透平叶片型线的设计方法 |
CN105373074A (zh) * | 2015-12-04 | 2016-03-02 | 哈尔滨汽轮机厂有限责任公司 | 基于型线调整的汽轮机叶片加工方法 |
CN106640210A (zh) * | 2016-09-14 | 2017-05-10 | 西安交通大学苏州研究院 | 一种径‑轴流式透平膨胀机叶轮叶片三维型线的设计方法 |
CN106682287A (zh) * | 2016-12-08 | 2017-05-17 | 西北工业大学 | 一种涡扇发动机内外涵道型线设计方法 |
-
2017
- 2017-09-27 CN CN201710889452.3A patent/CN107862104B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011243028A (ja) * | 2010-05-19 | 2011-12-01 | Ihi Corp | ターボ機械の翼形状設計方法およびプログラム |
CN104376144A (zh) * | 2014-08-06 | 2015-02-25 | 如皋透平叶片制造有限公司 | 一种透平叶片型线的设计方法 |
CN105373074A (zh) * | 2015-12-04 | 2016-03-02 | 哈尔滨汽轮机厂有限责任公司 | 基于型线调整的汽轮机叶片加工方法 |
CN106640210A (zh) * | 2016-09-14 | 2017-05-10 | 西安交通大学苏州研究院 | 一种径‑轴流式透平膨胀机叶轮叶片三维型线的设计方法 |
CN106682287A (zh) * | 2016-12-08 | 2017-05-17 | 西北工业大学 | 一种涡扇发动机内外涵道型线设计方法 |
Non-Patent Citations (3)
Title |
---|
GEORGIA N. KOINI等: "A software tool for parametric design of turbomachinery blades", 《ADVANCES IN ENGINEERING SOFTWARE》 * |
SHUANGTAO CHEN等: "Effect of impeller blade profile on the cryogenic two-phase turbo expander performance", 《APPLIED THERMAL ENGINEERING》 * |
刘城等: "基于贝塞尔曲线的液力变矩器三维叶片造型方法", 《机械工程学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114962322A (zh) * | 2022-05-24 | 2022-08-30 | 威灵(芜湖)电机制造有限公司 | 叶轮结构、风机及家用电器 |
CN114962322B (zh) * | 2022-05-24 | 2023-09-26 | 威灵(芜湖)电机制造有限公司 | 叶轮结构、风机及家用电器 |
CN115329506A (zh) * | 2022-10-14 | 2022-11-11 | 中国空气动力研究与发展中心空天技术研究所 | 用于自由射流试验的动力装置安装架的设计方法及安装架 |
CN115329506B (zh) * | 2022-10-14 | 2022-12-16 | 中国空气动力研究与发展中心空天技术研究所 | 用于自由射流试验的动力装置安装架的设计方法及安装架 |
CN117235937A (zh) * | 2023-11-14 | 2023-12-15 | 广东顺威精密塑料股份有限公司 | 一种多翼离心风机叶轮的周向弯曲叶片的设计方法 |
CN117235937B (zh) * | 2023-11-14 | 2024-03-26 | 广东顺威精密塑料股份有限公司 | 一种多翼离心风机叶轮的周向弯曲叶片的设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107862104B (zh) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107862104A (zh) | 一种环形喷嘴直线‑弧线叶片型线设计方法 | |
CN104907617B (zh) | 基于分区域切削的离心压缩机叶轮五坐标铣削方法 | |
CN104537191B (zh) | 一种三维钢筋模型示意图的绘制方法 | |
CN106640210A (zh) | 一种径‑轴流式透平膨胀机叶轮叶片三维型线的设计方法 | |
CN104166374B (zh) | 一种不可展直纹曲面五轴侧铣双参数刀位优化方法 | |
CN107288690B (zh) | 一种基于曲线参数化的涡轮叶片气膜孔形位参数修正方法 | |
CN102054099B (zh) | 一种贯流式锥形空间导叶三维密合设计方法和装置 | |
CN103994195B (zh) | 基于儒科夫斯基型线的液力变矩器叶片造型方法 | |
CN106502202A (zh) | 一种球头铣刀与导向叶片接触区域的半解析建模方法 | |
CN109570591A (zh) | 离心叶轮切削加工方法及装置和离心叶轮加工设备 | |
Li et al. | Influence of Coanda surface curvature on performance of bladeless fan | |
CN114036668B (zh) | 一种基于离心叶轮中间面曲线的可变自由曲面叶片生成方法 | |
CN112001131B (zh) | 一种改进风力机尾流叠加计算方法 | |
CN104392053A (zh) | 一种蒙皮滚弯零件截面曲率分析方法 | |
CN108876911A (zh) | 一种超声速流道三维流场计算方法及装置 | |
CN103902786B (zh) | 一种涡轮气冷动叶片伸根段外形参数化设计方法 | |
CN106844963B (zh) | 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法 | |
CN106780742A (zh) | 一种用于蜗壳结构的三维布筋的方法 | |
CN106089806A (zh) | 一种减小扩压器分离损失的端壁处理方法 | |
CN105401982B (zh) | 半转速核电汽轮机用末级动叶片叶身结构 | |
CN105179322A (zh) | 叶根开设等宽直线槽的压气机静子叶栅 | |
CN108399291A (zh) | 一种用于叶轮机械有叶部件内部流道的网格划分方法 | |
CN104573207B (zh) | 基于保形变换的液力变矩器直纹叶片三维造型方法 | |
CN103413019A (zh) | 一种不规则轮廓双曲率外形面整体壁板离散方法 | |
CN103870626A (zh) | 一种径轴流式透平膨胀机叶轮子午面型线设计及校核方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
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 |