涡轮叶片是涡轮机械的主要组成部分, 是航空发动机最重要的零件之一。新一代航空发动机进气口温度高达1 600~1 900℃, 单晶+复合气膜冷却技术是涡轮叶片的主流结构形式, 采用定向凝固、激光脉冲加工、涂层、超精密磨削等制造工艺, 工序复杂, 变形难以控制[1]。目前, 三坐标测量是涡轮叶片外形尺寸的主要检测方式, 测量精度高、灵活性强, 是目前最可信的涡轮叶片数字化检测方法[2]。
采样点规划是三坐标测量的重要组成部分, 也是自由曲线曲面测量的难点。很多学者对自由曲面曲线采样规划方法做了大量研究, 出现了随机采样[3]、曲率/曲率矩加权采样法[4-6]、细分采样法[7-10]、基于面片曲率的采样法[11]、基于面片大小的采样法[12]、混合采样法[13]、删点采样法[14-15]、二次采样法[16]、弦公差采样法[17-18]等。几乎所有研究都是针对轮廓度误差、弯扭度误差或逆向建模等, 很少研究叶片采样点对叶型参数的影响。而涡轮叶片特征参数是有明确的物理意义的, 是叶片设计制造过程中必须考虑的因素。
涡轮叶片的设计方法可分为2类:①基于控制点, 即按照气动载荷与性能要求, 通过求解微分方程获得叶盆、叶背上关键控制点的坐标, 根据这些控制点, 由分段圆弧、多项式曲线、Bezier曲线或B样条曲线等插值生成叶型线; ②基于叶型特征参数, 由中弧线控制参数及壁厚控制参数生成叶片的中弧线模型及厚度模型, 通过圆心在中弧线上的内切圆构成包络线来生成叶型线。基于叶型特征参数的涡轮叶片设计, 可描述叶片变形与特征参数之间的关系, 能够可靠地设计出高性能的叶型, 越来越受到研究者的重视, 结合叶型参数化设计及计算流体动力学理论出现很多涡轮机械设计软件, 如TURBODesign-1, PropCad, FINE/Design3DTM, BladePro, BladeCAD, AXCAD, T4T等[19]。
涡轮叶片叶型特征参数与其气动性能直接相关, 特征参数检测是叶片尺寸检测的重要内容之一。而三坐标测量点对计算涡轮叶片叶型特征参数的影响尚处于空白。因此, 研究不同三坐标测量点对涡轮叶片叶型特征参数检测精度的影响是非常有意义的。
1 叶型特征参数涡轮叶片叶型特征参数分为中弧线控制参数:内弦长、安装角、最大挠度、最大挠度位置、前缘几何进气角、后缘几何出气角; 厚度控制参数:最大厚度、最大厚度位置、前缘半径、后缘半径、前缘楔形角、后缘楔形角。如图 1所示, 其中L为弦长, r1为前缘半径, r2为后缘半径, rmax为最大厚度, Fmax为最大挠度, Pmax为最大挠度位置, β为安装角, α1为前缘进气角, α2为后缘出气角, γ1为前缘楔形角, γ2为后缘楔形角。
叶型特征参数提取[20]分为以下几步:
1) 去除噪声点:
a) 手动删除明显噪声点;
b) 删除凸侧噪声:设Pi-1, Pi, Pi+1为3个相邻测量点, ki-1, ki, ki+1分别为点Pi-1, Pi, Pi+1的估算曲率, Vi-1, Vi, Vi+1分别为点Pi-1, Pi, Pi+1指向曲率圆心的矢量。如果ki满足(1)式则Pi为凸侧噪声点, 删除Pi, 并重新计算ki-1, ki+1;
(1) |
c) 删除凹侧噪声:满足(2)式, 则Pi为凹侧坏点, 删除Pi, 并重新计算Vi-1, Vi+1;
(2) |
2) 弦线的提取:通过凸包算法, 以叶型测量点凸包多边形最长边为弦线, 其长度为弦长, 其与发动机轴线夹角为安装角。
3) 前后缘半径的提取:通过最小二乘圆弧拟合算法, 选取前(后)缘测量点, 记为(xi, yi) i=1, 2, …, n, 构建矩阵A,B,T, 则前(后)缘圆心为(c1, c2), 前(后)缘半径为
(3) |
(4) |
(5) |
4) 中弧线的提取:通过叶背(盆)点, 使用三次B样条插值算法, 以前后缘与叶背(盆)分界点作为端点约束, 以分界点切矢作为端点切矢约束, 构建叶背(盆)曲线, 然后采用等半径法求中弧线[21]。在中弧线提取过程中可得到最大厚度、最大厚度位置、最大挠度、最大挠度位置、前缘进气角、后缘出气角、前缘楔形角、后缘楔形角。
结合上述叶型参数提取算法, 采用VC++6.0以及VTK(visualization tool kit)可视化工具包, 开发了涡轮叶片尺寸测量软件TBGeoInspect, 界面如图 2所示。
2 自由曲线采样算法自由曲线采样方法可以分为3类:①删点法:是在密集均匀采样点的基础上, 以法矢、切矢或相邻采样点连接直线的夹角为判定条件删除不必要的点; ②增点法:在采样关键点的基础上, 以最小采样密度、弦公差或拟合误差为采样准则, 迭代增加采样点;③直接法:直接对曲线按一定准则计算采样点。删点法主要依据曲率相关准则删除点, 与曲率采样结果近似, 对叶型的采样效果并不好[22]; 增点法需要计算叶型关键点, 而且逐步迭代增加采样点, 计算过程复杂。直接法直接针对曲线形状, 几何意义明确, 计算简单, 通过调整采样准则可方便的调整采样点的分布, 在自由曲面测量中应用较广。杠杆平衡原理[5]在自由曲线直接采样中有广泛的应用, 如图 3所示, 杠杆平衡条件满足
(6) |
式中, mi是点i和点i+1之间曲线段的“质量”。
根据mi的不同计算方式, 可实现不同的采样方法:
1) 均匀采样(包含等参数采样和等弧长采样):
等参数采样以曲线段的参数增量作为曲线段质量, 如(7)式所示
(7) |
式中, ti和ti+1分别是点i和点i+1的曲线参数。
等弧长采样以曲线段的弧长作为曲线段质量, 可由参数曲线的弧长微分公式表示, 如(8)式所示
(8) |
式中,
2) 弦公差采样:
以曲线段的弦高作为曲线段质量, 弦高可由(9)式[17]表示
(9) |
式中, k(ti)为参数ti曲线的曲率, li为点i和点i+1间弧线的长度。
3) 曲率采样:
以曲线段的曲率积分作为曲线段质量, 如(10)式所示
(10) |
式中, k(t)为参数t处曲线的曲率。
4) 加权曲率采样[5]:
以曲线段的加权归一化曲率作为曲线段质量, 如(11)式所示
(11) |
式中, kmax和kmin分别是曲线C(t)的最大曲率和最小曲率, λ为权重系数。
5) 曲率-弧长采样[23]:
以曲线段曲率比重和弧长比重的加权和作为曲线段质量, 如(12)式所示
(12) |
式中, t为曲线参数, t∈[0, 1];λ为权重系数。
加权曲率采样和曲率-弧长采样算法中, 通过调整权重系数λ的大小, 可调整曲率和弧长的比重。λ值越大, 越接近于等弧长采样, λ值越小, 越接近于曲率采样。
3 仿真讨论以某叶型截面线为例, 分别采用均匀采样(等弧长)、弦公差采样、曲率采样、加权曲率采样(λ=0.03)和曲率-弧长采样(λ=0.5)规划测量点, 算例中的权重系数λ是以拟合误差为目标函数的优选值, 采样点分布如图 4所示:
选取前缘半径、后缘半径、弦长、最大厚度、最大挠度点作为叶型的主要特征参数。采用上述5种采样算法分别规划20~200个测量点, 计算特征参数并与理论值比较, 结果如表 1所示, 其中, 最大挠度点误差以最大挠度点的波动半径表示。对测量点进行三次B样条插值并与理论曲线对比, 插值误差如图 5所示。
mm | |||||
采样 方法 | 前缘半 径误差 | 后缘半 径误差 | 弦长 误差 | 最大厚 度误差 | 最大挠度 点误差 |
均匀 采样 | -0.022~ 0.505 | -0.278 ~0.005 | -0.422~ 0.002 | 0 ~0.002 | 2.728 3 |
曲率 采样 | -0.005~ 0.072 | -0.01~ 0.044 | -0.013~ 0.05 | -0.001~ 0.027 | 0.915 0 |
弦公差 采样 | -0.055~ 0.089 | -0.025~ 0.001 | -0.005~ 0.034 | -0.001~ 0.019 | 0.575 6 |
加权曲 率采样 | -0.026~ 0.119 | -0.01~ 0.093 | -0.001~ 0.047 | -0.001~ 0.020 | 0.115 8 |
曲率弧 长采样 | -0.024~ 0.063 | 0.01~ 0.016 | 0~ 0.05 | -0.001~ 0.003 | 0.074 6 |
通过表 1可以看出, 5种采样算法对最大厚度的计算精度都比较高; 均匀采样由于在前后缘点分布密度不够, 导致前后缘识别不是很准确, 从而影响前后缘半径及弦长的计算精度; 曲率采样由于大部分采样点都分布在前后缘, 对前后缘半径的计算精度最高; 曲率-弧长采样对5个关键叶型特征参数的计算精度都比较高, 而且对最大挠度点的计算精度明显高于其他采样方法。从图 5a)最大插值曲线误差和图 5b)平均插值曲线误差可以看出, 采样点数小于100时, 曲率-弧长采样点的曲线插值误差小于加权曲率采样点和弦公差采样点, 明显小于均匀采样点和曲率采样点; 曲率-弧长采样算法能用更少的点表征叶型曲线的曲率形状信息。因此, 曲率-弧长采样更适合叶片采样点规划。
4 测量实验以涡轮导叶精铸件为例, 如图 6所示, 取A-A截面, 分别采用均匀采样、曲率采样、弦公差采样、曲率-弧长采样、加权曲率采样规划50, 100, 150, 200个采样点。通过专用夹具固定导叶, 使用三坐标测量机(型号型号:GLOBAL D STATUS 121510)触测采样点。对采样点的测量值去噪声处理后进行三次B样条插值, 并导入到TBGeoInspect软件中, 计算A-A截面的叶型参数, 结果如表 2所示。
mm | |||||
采样 方法 | 前缘 半径 | 后缘 半径 | 弦长 | 最大 厚度 | 最大挠度 点误差 |
均匀 采样 | 3.941~ 4.746 | 0.623~ 1.079 | 80.286~ 80.479 | 15.377~ 15.397 | 0.894 |
曲率 采样 | 4.148~ 4.787 | 0.847~ 0.954 | 80.364~ 80.385 | 15.355~ 15.389 | 0.972 |
弦公差 采样 | 4.005~ 4.824 | 0.931~ 0.987 | 80.346~ 80.379 | 15.358~ 15.409 | 0.619 |
加权曲率 采样 | 4.053~ 4.771 | 0.870~ 0.924 | 80.324~ 80.388 | 15.372~ 15.397 | 0.645 |
曲率弧长 采样 | 4.600~ 4.771 | 0.871~ 0.945 | 80.371~ 80.386 | 15.371~ 15.400 | 0.371 |
图 7中, A, B, C, D分别为前缘半径波动范围、后缘半径波动范围、弦长波动范围和最大壁厚波动范围; E为最大挠度点波动半径, 用所求最大挠度点最小外接圆的半径表示。通过表 2和图 7a)~e)可以看出, 弦长和最大厚度的波动误差最小; 前缘半径的波动误差大于后缘半径的波动误差; 曲率-弧长采样计算叶型参数的波动误差明显小于其它4种采样算法。
5 结论本文通过仿真与实验, 研究了均匀采样、曲率采样、弦公差采样、加权曲率采样、曲率-弧长采样点对叶型特征参数及轮廓误差计算精度的影响, 可得出如下结论:
1) 最大厚度和弦长的计算精度受采样点分布的影响较小;
2) 前缘半径受采样点分布的影响大于后缘半径;
3) 曲率-弧长采样对叶型特征参数计算的误差小于均匀采样、曲率采样、弦公差采样和加权曲率采样;
4) 采样点数相同时, 曲率-弧长采样点的插值误差小于加权曲率采样和弦公差采样, 明显小于均匀采样和曲率采样。
[1] | QIU F, BU K, SONG J H, et al. Dimensional Control of Nickel-Based Single Crystal Turbine Blade Investment Casting by Process Control Optimization[J]. International Journal of Metalcasting, 2018, 12(3): 469-479. DOI:10.1007/s40962-017-0180-5 |
[2] |
蔺小军, 郭研, 吴广, 等. 等高法叶片型面CMM测量数据处理算法研究[J]. 仪器仪表学报, 2013, 34(11): 2442-2450.
LIN Xiaojun, GUO Yan, WU Guang, et al. CMM Measuring Data Processing Algorithms for Blades about the Contour Measurment[J]. Chinese Journal of Scientific Instrument, 2013, 34(11): 2442-2450. (in Chinese) |
[3] | LEE G, MOU J, SHEN Y. Sampling Strategy Design for Dimensional Measurement of Geometric Features Using Coordinate Measuring Machine[J]. International Journal of Machine Tools & Manufacture, 1997, 37(7): 917-934. |
[4] |
何改云, 季彩霞, 郭龙真, 等. 基于等矩策略的自由曲面测量方法误差分析[J]. 天津大学学报, 2016, 49(2): 128-133.
HE Gaiyun, JI Caixia, GUO Longzhen, et al. Analysis of Measurement Method Error in Freeform Surface Based on Equal Moment Strategy[J]. Journal of Tianjin University, 2016, 49(2): 128-133. (in Chinese) |
[5] | LI S Z. Adaotive Sampling and Mesh Generation[J]. Computer-Aided Design, 1995, 27(3): 235-240. DOI:10.1016/0010-4485(95)95872-C |
[6] | HE G Y, JIA H Y, GUO L Z, et al. Adaptive Sampling Strategy for Free-Form Surface Based on CAD Model[J]. Automatic Manufacturing Systems, 2012, 1(2): 541-544. |
[7] | YU M R, ZHANG Y J, LI Y L, et al. Adaptive Sampling Method for Inspection Planning on CMM for Free-Form Surfaces[J]. International Journal of Advanced Manufacturing Technology, 2013, 67(9/10/11/12): 1967-1975. |
[8] | AINSWORTH I, RISTIC M, BRUJIC D. CAD-Based Measurement Path Planning for Free-Form Shapes Using Contact Probes[J]. International Journal of Advanced Manufacturing Technology, 2000, 16(1): 23-31. |
[9] | REN M J, KONG L B, SUN L J, et al. A Curve Network Sampling Strategy for Measurement of Freeform Surfaces on Coordinate Measuring Machines[J]. IEEE Trans on Instrumentation and Measurement, 2017, 66(11): 3032-3043. DOI:10.1109/TIM.2017.2717283 |
[10] | OBEIDAT S M, RAMAN S. An Intelligent Sampling Method for Inspecting Free-Form Surfaces[J]. International Journal of Advanced Manufacturing Technology, 2009, 40(11/12): 1125-1136. |
[11] | RAJAMOHAN G, SHUNMUGAM M S, SAMUEL G L. Practical Measurement Stratigies for Verification of Freeform Surfaces Using Coordinate Measuring Machines[J]. Metrology and Measurement Systems, 2011, 18(2): 209-222. |
[12] | ELKOTT D F, VELDHUIS S C. Isoparametric Line Sampling for the Inspection Planning of Sculptured Surfaces[J]. Computer-Aided Design, 2005, 37(2): 189-200. DOI:10.1016/j.cad.2004.06.006 |
[13] | ELKOTT D F, ELMARAGHY H A, Elmaraghy W H. Automatic Sampling for CMM Inspection Planning of Free-Form Surfaces[J]. International Journal of Production Research, 2002, 40(11): 2653-2676. DOI:10.1080/00207540210133435 |
[14] |
龚玉玲, 平雪良, 周保珍, 等. 基于矢量相似淘汰算法的CMM测量路径规划及仿真检验[J]. 工具技术, 2008(1): 114-118.
GONG Yuling, PING Xueliang, ZHOU Baozhen, et al. Planning and Simulation of CMM Measuring Path Based on Vector Resemble Selection Algorithm[J]. Tool Engineering, 2008(1): 114-118. (in Chinese) DOI:10.3969/j.issn.1000-7008.2008.01.037 |
[15] |
周保珍. 曲面测量方法研究及测量仿真[D]. 无锡: 江南大学, 2008 ZHOU Baozhen. Study on Measuring Method of Surface and Measuring Simulation[D]. Wuxi, Jiangnan University, 2008(in Chinese) http://cdmd.cnki.com.cn/Article/CDMD-10295-2009013791.htm |
[16] |
潘金川. 整体叶轮的CMM测量规划与仿真技术研究[D]. 南京: 南京航空航天大学, 2014 PAN Jinchuan. Research on Planning and Simulation of CMM Measuring the Integral Impeller[D]. Nanjing, Nanjing University of Aeronautics and Astronautics, 2014(in Chinese) http://cdmd.cnki.com.cn/Article/CDMD-10287-1014061281.htm |
[17] | JIANG R S, WANG W H, ZHANG D H, et al. A Practical Sampling Method for Profile Measurement of Complex Blades[J]. Measurement, 2016, 81: 57-65. DOI:10.1016/j.measurement.2015.11.039 |
[18] |
程云勇, 朱真真, 李维亮, 等. 面向叶片弯扭变形分析的测量采样方法[J]. 计算机集成制造系统, 2013, 19(6): 1329-1335.
CHENG Yunong, ZHU Zhenzhen, LI Weiliang, et al. Measurement Sampling Method Oriented to Blade Bending and Torsion Deformation Analysis[J]. Computer Integrated Manufacturing Systems, 2013, 19(6): 1329-1335. (in Chinese) |
[19] | KOINI G N, SARAKINOS S S, NIKOLOS I K. A Software Tool for Parametric Design of Turbomachinery Blades[J]. Advances in Engineering Software, 2009, 40(1): 41-51. DOI:10.1016/j.advengsoft.2008.03.008 |
[20] |
陈志强. 基于测量数据的叶片截面特征参数提取技术研究[D]. 西安: 西北工业大学, 2007 CHEN Zhiqiang. Research of Cross-Sectional Feature Parameter Extraction of Blade Based on Measurement Data[D]. Xi'an, Northwestern Polytechnical University, 2007(in Chinese) http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y1034501 |
[21] |
张丹, 张卫红, 万敏, 等. 基于位移场仿真与特征参数提取的精铸模具型面逆向设计方法[J]. 航空学报, 2006, 27(3): 509-514.
ZHANG Dan, ZHANG Weihong, WAN Min, et al. Reversing Design Methodology of the Die Profile in Investment Casting Based on the Simulation of Displacement Field and Identification of Featured Parameters[J]. ACTA Aeronautica et Astronautica Sinica, 2006, 27(3): 509-514. (in Chinese) DOI:10.3321/j.issn:1000-6893.2006.03.031 |
[22] |
张现东, 卜昆, 程云勇, 等. 基于CMM的叶片测点采样分布方法研究[J]. 锻压技术, 2012, 37(4): 131-135.
ZHANG Xiandong, BU Kun, CHENG Yunyong, et al. Research on CMM-Based Measuring Points Sampling and Distribution Method of Blades[J]. Forging & Stamping Technology, 2012, 37(4): 131-135. (in Chinese) DOI:10.3969/j.issn.1000-3940.2012.04.032 |
[23] | PAGANI L, SCOTT P J. Curvature Based Sampling of Curves and Surfaces[J]. Computer Aided Geometric Design, 2018, 59: 32-48. DOI:10.1016/j.cagd.2017.11.004 |