一种压气机叶片中弧线计算方法技术

一种压气机叶片中弧线计算方法，属于燃气轮机技术领域。解决了现有计算方法迭代次数多，且在压气机叶片前缘、尾缘附近往往误差较大，结果不理想的技术问题。本发明专利技术的技术要点：直接在背弧上求得一点E，使得该点的横坐标为C、D两点横坐标的平均值。然后在法线上寻找一点O，使得O到点A、点E的距离相等。以O为圆心，||OA||为半径作圆，根据该圆与背弧的交点个数判断点O是否为内切圆圆心。如若不是，以新的C、D点构成的圆弧为搜索范围继续寻找。每循环一次都会将目标范围缩小，一般只需要迭代3‑5次即可找到中弧线上的点，缩短了计算时间，提高了计算效率。

【技术实现步骤摘要】
一种压气机叶片中弧线计算方法
本专利技术涉及一种叶片中弧线计算方法，具体涉及一种压气机叶片中弧线计算方法，属于燃气轮机

技术介绍
叶片是燃气轮机发动机的重要零件，它在高温、高压、高转速状态下运转，决定了其工作环境的恶劣。为满足发动机高性能、高可靠性及长寿命的要求，叶片必须具有精确的尺寸、准确的形状。尤其是在逆向工程领域，根据叶片的测量结果反求出叶片截面的各项设计参数是必须的，其中进口几何角、出口几何角、前缘角、尾缘角、最大挠度、最大厚度等参数的求解都需要首先知道叶片截面的中弧线。因此，中弧线作为提取叶片截面参数的依据，对叶型质量检测及叶片造型都具有非常重要的影响。目前已有的提取叶片中弧线的方法有等半径法、等角度法。等半径法和等角度法类似，都是应用内切圆的圆心到叶盆、叶背上的距离相等原理，利用叶盆或叶背上的一点法线与另一型线有一交点，若两点的切线平行，两点的中点就是内切圆圆心，若不平行，则过中点做另一交点处切线的垂线，判断垂线与另一型线的交点和法线与另一型线的交点之间的距离小于给定值时，认为求得中弧线上的点，否则继续重复这个过程。这种计算方法迭代次数多，且在压气机叶片前缘、尾缘附近往往误差较大，结果不理想。现有技术中：申请人为上海电气电站设备有限公司，申请号为CN201010148017.3，专利技术创造名称为叶轮机械叶片中弧线计算方法的专利技术专利申请，其公开了如下步骤：(1)读入叶片型线离散数据，将叶片型线划分为内弧、背弧、前缘和尾缘四部分；(2)将内弧、背弧用三次样条曲线拟合；(3)在背弧上选取一点A，求出背弧在点A的法线以及该法线与内弧的交点C；(4)以线段AC为直径设置圆O2，并求出圆O2与内弧的另一个交点D，如点C与点D重合，则圆O2的圆心即为中弧线上的点，此时转入步骤(8)，如点C与点D不重合，则转入步骤(5)；(5)按照某一步长值inc，将内弧上点C和点D之间曲线细分成n-1段，记这些分点为E2、E3、…、En-1，并用点E1表示点C，点En表示点D；(6)过E1作内弧的法线，计算该法线与AC的交点Q1的坐标，如果线段E1Q1和线段AQ1的长度差值小于预先设定的误差error，则认为Q1就是所要求的中弧线上的点，如果E1不能满足要求，则依次计算E2至En各点，直至得到满足误差要求的中弧线上的点，并转入步骤(8)；(7)如果一直计算到En都没有得到满足误差要求的中弧线上的点，则减小步长inc，重复步骤(5)和(6)，直到求得达到误差要求的中弧线上的点；(8)对背弧上所有的原有数据点，按照步骤(3)到步骤(7)求中弧线上的点；(9)将求得的所有中弧线上的离散点用三次样条曲线拟合，检验每一个点的斜率及曲率是否存在突然跳跃，将存在突然跳跃的点舍去，剩余的点就构成了中弧线。其依然存在迭代次数多，计算时间长，计算效率低的问题。
技术实现思路
在下文中给出了关于本专利技术的简要概述，以便提供关于本专利技术的某些方面的基本理解。应当理解，这个概述并不是关于本专利技术的穷举性概述。它并不是意图确定本专利技术的关键或重要部分，也不是意图限定本专利技术的范围。其目的仅仅是以简化的形式给出某些概念，以此作为稍后论述的更详细描述的前序。鉴于此，为了克服现有计算方法迭代次数多，且在压气机叶片前缘、尾缘附近往往误差较大，结果不理想的技术问题，本专利技术提供了一种压气机叶片中弧线计算方法。本专利技术的一种压气机叶片中弧线计算方法，具体步骤为：步骤一，读入叶片截面离散数据，划分为叶盆、叶背、前缘圆弧、尾缘圆弧四部分；步骤二，应用最小二乘法分别拟合叶盆曲线fy1、叶背曲线fy2；步骤三，叶盆曲线fy1上选取一点A，求fy1在A点处法线ff1，与叶背曲线fy2交于点B，以线段AB中点O’为圆心，||AB||为直径作圆；步骤四，若该圆与fy2交于一点C且B点与C点重合，则该圆为内弧、背弧共有的内切圆，O’为内切圆圆心，即中弧线上的一点；步骤五，若该圆与fy2交于两点C、D，且两点之间的距离||CD||<delta，其中delta为设置的阈值，认为该圆为内弧、背弧共有的内切圆，O’为内切圆圆心，即中弧线上的一点；步骤六，若该圆与fy2交于两点C、D，且两点之间的距离||CD||>delta，在fy2上找一点E，使得该点的横坐标为C、D两点横坐标的平均值；在法线ff1上寻找一点O，使得O到点A、点E的距离相等，即||OA||＝||OE||；以O为圆心，||OA||为半径作圆；步骤七，判断该圆与fy2的交点个数，重复步骤四到步骤六，直到找到内弧、背弧共有的内切圆，其圆心O即是中弧线上的一点；步骤八，对内弧上所有点，按照步骤四到步骤七的方法，基于叶盆曲线、叶背曲线获得中弧线上所有的点；步骤九，将求得的中弧线上所有的点用三次样条曲线拟合，检验每一个点斜率、曲率是否存在突然跳跃，若有舍去突跃点，结束。进一步地：步骤二中叶盆曲线fy1、叶背曲线fy2的拟合采用移动的最小二乘法，即每次拟合仅取几个连续的离散点，既要保证曲线拟合的精度，又要确保A点处法线ff1与叶背曲线fy2的交点落在叶背曲线fy2上。进一步地：步骤三中叶盆曲线fy1在A点处的内切圆圆心落在A点处法线ff1上。进一步地：步骤四中以线段AB中点O’为圆心，线段AB的长度为直径作圆，该圆是叶盆曲线fy1在A点处的内切圆，且与fy2至少有一个交点。进一步地：步骤六中以O为圆心，||OA||为半径作圆时，该圆与叶背曲线fy2有一个或两个交点，点E必为其中一个。进一步地：以O为圆心，||OA||为半径作圆，该圆与叶背曲线fy2有两个交点(C’、D’)时，再次寻找叶背曲线fy2上一点E’时，C、D、C’、D’四点按横坐标排序，应以中间两点横坐标的平均值作为E’的横坐标。如此重复，不断缩小寻找范围。在上述方法中，也可将A点选在叶背曲线fy2上，采用相同的方法获得叶盆曲线上的C、D、E点，确定叶盆、叶背共有的内切圆，进而获得内切圆的圆心，即中弧线上的点。有益效果：本专利技术通过叶盆上某一点A的法线，与叶背相交于点B，以AB为直径作圆，交叶背于两点C、D，通过C、D两点横坐标平均值找到叶背曲线上点E，在法线上寻找与点A、E距离相等的点的方法寻找内切圆圆心。一般循环迭代3-5次即可找到中弧线上的点，缩短了计算时间，提高了计算效率，该方法简单、快捷，可满足计算精度和速度的要求。附图说明图1为本专利技术所述一种压气机叶片中弧线计算方法流程图；图2为叶盆、叶背共有内切圆圆心计算原理示意图；图3为压气机叶片截面中弧线计算结果示意图。具体实施方式在下文中将结合附图对本专利技术的示范性实施例进行描述。为了清楚和简明起见，在说明书中并未描述实际实施方式的所有特征。然而，应该了解，在开发任何这种实际实施例的过程中必须做出很多特定于实施方式的决定，以便实现开发人员的具体目标，例如，符合与系统及业务相关的那些限制条件，并且这些限制条件可能会随着实施方式的不同而有所改变。此外，还应该了解，虽然开发工作有可能是非常复杂和费时的，但对得益于本专利技术公开内容的本领域技术人员来说，这种开发工作仅仅是例行的任务。在此，还需要说明的一点是，为了避免因不必要的细节而模糊了本专利技术，在附图中仅仅示出了与根据本专利技术的方案密切相关的装置结构和/或处理步骤，而省略本文档来自技高网...

【技术保护点】
1.一种压气机叶片中弧线计算方法，其特征在于：具体步骤如下，步骤一，读入叶片截面离散数据，划分为叶盆、叶背、前缘圆弧、尾缘圆弧四部分；步骤二，应用最小二乘法分别拟合叶盆曲线fy1、叶背曲线fy2；步骤三，叶盆曲线fy1上选取一点A，求fy1在A点处法线ff1，与叶背曲线fy2交于点B，以线段AB中点O’为圆心，||AB||为直径作圆；步骤四，若该圆与fy2交于一点C且B点与C点重合，则该圆为内弧、背弧共有的内切圆，O’为内切圆圆心，即中弧线上的一点；步骤五，若该圆与fy2交于两点C、D，且两点之间的距离||CD||<delta，其中delta为设置的阈值，认为该圆为内弧、背弧共有的内切圆，O’为内切圆圆心，即中弧线上的一点；步骤六，若该圆与fy2交于两点C、D，且两点之间的距离||CD||>delta，在fy2上找一点E，使得该点的横坐标为C、D两点横坐标的平均值；在法线ff1上寻找一点O，使得O到点A、点E的距离相等，即||OA||＝||OE||；以O为圆心，||OA||为半径作圆；步骤七，判断该圆与fy2的交点个数，重复步骤四到步骤六，直到找到内弧、背弧共有的内切圆，其圆心O即是中弧线上的一点；步骤八，对内弧上所有点，按照步骤四到步骤七的方法，基于叶盆曲线、叶背曲线获得中弧线上所有的点；步骤九，将求得的中弧线上所有的点用三次样条曲线拟合，检验每一个点斜率、曲率是否存在突然跳跃，若有舍去突跃点，结束。...

【技术特征摘要】
1.一种压气机叶片中弧线计算方法，其特征在于：具体步骤如下，步骤一，读入叶片截面离散数据，划分为叶盆、叶背、前缘圆弧、尾缘圆弧四部分；步骤二，应用最小二乘法分别拟合叶盆曲线fy1、叶背曲线fy2；步骤三，叶盆曲线fy1上选取一点A，求fy1在A点处法线ff1，与叶背曲线fy2交于点B，以线段AB中点O’为圆心，||AB||为直径作圆；步骤四，若该圆与fy2交于一点C且B点与C点重合，则该圆为内弧、背弧共有的内切圆，O’为内切圆圆心，即中弧线上的一点；步骤五，若该圆与fy2交于两点C、D，且两点之间的距离||CD||<delta，其中delta为设置的阈值，认为该圆为内弧、背弧共有的内切圆，O’为内切圆圆心，即中弧线上的一点；步骤六，若该圆与fy2交于两点C、D，且两点之间的距离||CD||>delta，在fy2上找一点E，使得该点的横坐标为C、D两点横坐标的平均值；在法线ff1上寻找一点O，使得O到点A、点E的距离相等，即||OA||＝||OE||；以O为圆心，||OA||为半径作圆；步骤七，判断该圆与fy2的交点个数，重复步骤四到步骤六，直到找到内弧、背弧共有的内切圆，其圆心O即是中弧线上的一点；步骤八，对内弧上所有点，按照步骤四到步骤七的方法，基于叶盆曲线、叶背曲线获得中弧线上所有的点...

【专利技术属性】
技术研发人员：刘海旭，冯永志，张春梅，郑智文，孙涛，邓祥泉，李翔宇，王鑫，
申请(专利权)人：哈尔滨电气股份有限公司，
类型：发明
国别省市：黑龙江,23