CN110738728A - 基于线性组合过渡算法的叶片修复模型重建方法 - Google Patents

基于线性组合过渡算法的叶片修复模型重建方法 Download PDF

Info

Publication number
CN110738728A
CN110738728A CN201910954558.6A CN201910954558A CN110738728A CN 110738728 A CN110738728 A CN 110738728A CN 201910954558 A CN201910954558 A CN 201910954558A CN 110738728 A CN110738728 A CN 110738728A
Authority
CN
China
Prior art keywords
curve
transition
blade
linear combination
section
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
CN201910954558.6A
Other languages
English (en)
Other versions
CN110738728B (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.)
Northwest University of Technology
Original Assignee
Northwest University of Technology
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 Northwest University of Technology filed Critical Northwest University of Technology
Priority to CN201910954558.6A priority Critical patent/CN110738728B/zh
Publication of CN110738728A publication Critical patent/CN110738728A/zh
Application granted granted Critical
Publication of CN110738728B publication Critical patent/CN110738728B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Numerical Control (AREA)
  • Image Generation (AREA)

Abstract

本发明公开了一种基于线性组合过渡算法的叶片修复模型重建方法,用于解决现有叶片修复模型重建方法过渡区域连续性差的技术问题。技术方案是通过重节点插入技术将B样条基曲线退化成Bezier曲线;再通过Bezier曲线参数的一致性,结合基曲线线性组合过渡算法,实现B样条基曲线的光滑过渡;然后将所有过渡截面曲线实现高阶连续的过渡;最后放样得到重构叶片模型,该模型过渡区域与其他区域可以实现Gn连续,并可通过调控参数进行形状调控。

Description

基于线性组合过渡算法的叶片修复模型重建方法
技术领域
本发明涉及一种叶片修复模型重建方法,特别涉及一种基于线性组合过渡算法的叶片修复模型重建方法。
背景技术
文献“Repair of defective 3D blade model based on deformation ofadjacent non-defective cross-sectional curve,The International Journal ofAdvanced Manufacturing Technology.2018,95(5-8):3045-3055”公开了一种叶片修复模型重建方法。该方法针对非接触测量修复精度不高的问题,提出了一种接触测量法来进行叶片修复,并提出了基于参考截面曲线与缺陷截面曲线进行递推变形的修复算法,为了实现曲线的变形还提出了一种参考曲线上所有控制点的迭代移动分布算法,较好地重构了叶片的修复模型。但是文献所述方法没有考虑曲面过渡的连续性问题,如果重建模型过渡区域连续性过低,则加工中容易出现台阶问题。
发明内容
为了克服现有叶片修复模型重建方法过渡区域连续性差的不足,本发明提供一种基于线性组合过渡算法的叶片修复模型重建方法。该方法通过重节点插入技术将B样条基曲线退化成Bezier曲线;再通过Bezier曲线参数的一致性,结合基曲线线性组合过渡算法,实现B样条基曲线的光滑过渡;然后将所有过渡截面曲线实现高阶连续的过渡;最后放样得到重构叶片模型,该模型过渡区域与其他区域可以实现Gn连续,并可通过调控参数进行形状调控。
本发明解决其技术问题所采用的技术方案:一种基于线性组合过渡算法的叶片修复模型重建方法,其特点是包括以下步骤:
步骤一、将增材加工后的叶片进行三坐标测量,得到测量点集并进行测量数据预处理。
步骤二、将测量点集拟合成曲线,并从设计模型上端提取截面曲线。
步骤三、测量截面曲线与设计截面曲线沿叶高方向进行过渡。将测量截面曲线与设计截面曲线沿叶高方向重新拟合成待过渡基曲线。通过重节点插入方法将基曲线待过渡端退化成Bezier曲线。利用基曲线线性组合过渡算法对基曲线进行过渡,得到过渡曲线。
γ(t)=f(t)γ1(t)+(1-f(t))γ2(t),t∈[0,1]
式中,t是曲线参数,γ(t)为过渡曲线,γ1(t),γ2(t)分别为基曲线1和基曲线2,f(t)表示线性组合方法的过渡函数
Figure BDA0002226842910000021
式中,μ表示调控参数,t表示曲线参数,n表示曲线次数。
步骤四、将过渡曲线与设计前后缘进行桥接,得到完整截面曲线,重建叶片模型。
步骤五、判断重建模型是否符合性能要求,如果不符合,则调整参数,重新建模。
本发明的有益效果是:该方法通过重节点插入技术将B样条基曲线退化成Bezier曲线;再通过Bezier曲线参数的一致性,结合基曲线线性组合过渡算法,实现B样条基曲线的光滑过渡;然后将所有过渡截面曲线实现高阶连续的过渡;最后放样得到重构叶片模型,该模型过渡区域与其他区域可以实现Gn连续,并可通过调控参数进行形状调控。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明基于线性组合过渡算法的叶片修复模型重建方法的流程图。
图2是本发明方法中线性组合过渡算法高阶连续的证明图。
具体实施方式
参照图1-2。本发明基于线性组合过渡算法的叶片修复模型重建方法具体步骤如下:
步骤一、将增材加工后的叶片进行三坐标测量,得到测量点集并进行测量数据预处理。
步骤二、将测量点集拟合成曲线,并从设计模型上端提取截面曲线。
步骤三、测量截面曲线与设计截面曲线沿叶高方向进行过渡。将测量截面曲线与设计截面曲线沿叶高方向重新拟合成待过渡基曲线。通过重节点插入方法将基曲线待过渡端退化成Bezier曲线。利用基曲线线性组合过渡算法对基曲线进行过渡,得到过渡曲线。
γ(t)=f(t)γ1(t)+(1-f(t))γ2(t),t∈[0,1]
式中,t是曲线参数,γ(t)为过渡曲线,γ1(t),γ2(t)分别为基曲线1和基曲线2,f(t)表示线性组合方法的过渡函数
Figure BDA0002226842910000031
式中μ表示调控参数,t表示曲线参数,n表示曲线次数。
步骤四、将过渡曲线与设计前后缘进行桥接,得到完整截面曲线,重建叶片模型。
步骤五、判断重建模型是否符合性能要求,如果不符合,则调整参数,重新建模。
应用实施例。
步骤一、获得测量点集。
将增材加工后的叶片通过三坐标测量机进行等截面测量,获得测量数据点,并对数据点进行处理。
步骤二、拟合测量截面曲线与抽取设计截面曲线。
将处理过的测量数据按截面拟合成曲线,再从设计模型上端抽取等V参数线,本实验从设计模型V=0.6到V=1.0区域抽取等V线10条。其中V为曲面参数方向。
步骤三、将测量曲线与设计曲线进行过渡。
将测量截面曲线与设计截面曲线沿叶高方向离散为30个点,再沿横截面方向分别拟合成曲线,然后计算出两条B样条曲线待过渡端的最后一段的控制顶点与节点矢量,接着插入重节点,使之达到满重复度k,k为B样条曲线的次数。本实验的曲线次数为3,即k=3,则插入两重节点即可将B样条曲线最后一段退化成Bezier曲线。然后通过线性组合方程r(t)=f(t)r1(t)+(1-f(t))r2(t),t∈[0,1]计算过渡曲线。其中r(t)为过渡曲线,r1(t)与r2(t)分别为退化的Bezier曲线。表达式为
Figure BDA0002226842910000034
Figure BDA0002226842910000032
式中Jk,i(t)是伯恩斯坦基函数,
Figure BDA0002226842910000033
di为控制点。得到过渡曲线点后,将之前离散的测量截面数据点,设计截面曲线点与过渡曲线点沿横截面向拟合成曲线。最后将拟合的曲线统一离散成30个点,沿叶高方向拟合成叶片横截面曲线。
步骤四、将重建截面曲线与设计前后缘进行桥接,生成完整模型。
从设计前后缘上抽取与重建截面曲线等高的截面曲线,再将重建叶片截面曲线与设计前后缘截面曲线进行桥接,得到完整的叶片截面曲线。将截面曲线放样得到重建叶片模型。
步骤五、判断重建模型是否满足要求
重建模型如果满足叶片的性能要求,或者重建模型的性能与设计性能偏差在允许范围内,则结束建模。如果不满足性能要求,则调整参数,重复步骤三到四,直至模型满足要求。
图2表明采用传统的叶片修复建模曲面过渡手段只能达到G1连续,而本发明的曲面过渡方法可以达到Gn连续,重建模型质量更好。

Claims (1)

1.一种基于线性组合过渡算法的叶片修复模型重建方法,其特征在于包括以下步骤:
步骤一、将增材加工后的叶片进行三坐标测量,得到测量点集并进行测量数据预处理;
步骤二、将测量点集拟合成曲线,并从设计模型上端提取截面曲线;
步骤三、测量截面曲线与设计截面曲线沿叶高方向进行过渡;将测量截面曲线与设计截面曲线沿叶高方向重新拟合成待过渡基曲线;通过重节点插入方法将基曲线待过渡端退化成Bezier曲线;利用基曲线线性组合过渡算法对基曲线进行过渡,得到过渡曲线;
γ(t)=f(t)γ1(t)+(1-f(t))γ2(t),t∈[0,1]
式中,t是曲线参数,γ(t)为过渡曲线,γ1(t),γ2(t)分别为基曲线1和基曲线2,f(t)表示线性组合方法的过渡函数
Figure FDA0002226842900000011
式中,μ表示调控参数,t表示曲线参数,n表示曲线次数;
步骤四、将过渡曲线与设计前后缘进行桥接,得到完整截面曲线,重建叶片模型;
步骤五、判断重建模型是否符合性能要求,如果不符合,则调整参数,重新建模。
CN201910954558.6A 2019-10-09 2019-10-09 基于线性组合过渡算法的叶片修复模型重建方法 Active CN110738728B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910954558.6A CN110738728B (zh) 2019-10-09 2019-10-09 基于线性组合过渡算法的叶片修复模型重建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910954558.6A CN110738728B (zh) 2019-10-09 2019-10-09 基于线性组合过渡算法的叶片修复模型重建方法

Publications (2)

Publication Number Publication Date
CN110738728A true CN110738728A (zh) 2020-01-31
CN110738728B CN110738728B (zh) 2022-05-03

Family

ID=69268493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910954558.6A Active CN110738728B (zh) 2019-10-09 2019-10-09 基于线性组合过渡算法的叶片修复模型重建方法

Country Status (1)

Country Link
CN (1) CN110738728B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113714731A (zh) * 2021-11-01 2021-11-30 中国航发沈阳黎明航空发动机有限责任公司 一种基于公差约束的叶片表面振纹修复方法
CN116663205A (zh) * 2023-08-01 2023-08-29 武汉华锋惠众科技有限公司 一种复杂工艺补充面截面线构造方法和转换方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103093065A (zh) * 2013-02-20 2013-05-08 天津工业大学 一种航空发动机叶片叶尖缺失部位三维模型重建方法
CN103372713A (zh) * 2013-02-20 2013-10-30 天津工业大学 一种航空发动机叶片损伤的修复方法
CN103486996A (zh) * 2013-08-14 2014-01-01 西北工业大学 未知cad模型的航空发动机叶片型面测量方法
EP3242177A1 (en) * 2016-04-29 2017-11-08 Rolls-Royce plc Adaptive repair method for aerofoil blades
CN108229081A (zh) * 2018-03-27 2018-06-29 华中科技大学 一种组合曲面的几何修复的方法
CN108544181A (zh) * 2018-03-27 2018-09-18 西北工业大学 一种整体叶盘损伤叶片的修复方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103093065A (zh) * 2013-02-20 2013-05-08 天津工业大学 一种航空发动机叶片叶尖缺失部位三维模型重建方法
CN103372713A (zh) * 2013-02-20 2013-10-30 天津工业大学 一种航空发动机叶片损伤的修复方法
CN103486996A (zh) * 2013-08-14 2014-01-01 西北工业大学 未知cad模型的航空发动机叶片型面测量方法
EP3242177A1 (en) * 2016-04-29 2017-11-08 Rolls-Royce plc Adaptive repair method for aerofoil blades
CN108229081A (zh) * 2018-03-27 2018-06-29 华中科技大学 一种组合曲面的几何修复的方法
CN108544181A (zh) * 2018-03-27 2018-09-18 西北工业大学 一种整体叶盘损伤叶片的修复方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
YING ZHANG等: "Aero engine and blade repair Path Generation Method for Aero-Engine Free-Form Surface Blade in Laser Solid Forming", 《IEEE ACCESS》 *
聂兆伟等: "航空发动机叶片自适应修复目标曲面重构", 《计算机集成制造系统》 *
董艇舰等: "一种受损叶片数字化3D模型重建方法研究", 《机械设计与制造》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113714731A (zh) * 2021-11-01 2021-11-30 中国航发沈阳黎明航空发动机有限责任公司 一种基于公差约束的叶片表面振纹修复方法
CN113714731B (zh) * 2021-11-01 2021-12-31 中国航发沈阳黎明航空发动机有限责任公司 一种基于公差约束的叶片表面振纹修复方法
CN116663205A (zh) * 2023-08-01 2023-08-29 武汉华锋惠众科技有限公司 一种复杂工艺补充面截面线构造方法和转换方法
CN116663205B (zh) * 2023-08-01 2023-09-22 武汉华锋惠众科技有限公司 一种复杂工艺补充面截面线构造方法和转换方法

Also Published As

Publication number Publication date
CN110738728B (zh) 2022-05-03

Similar Documents

Publication Publication Date Title
CN108549325B (zh) 一种自由曲面弧长参数曲线加工轨迹生成方法
CN110738728B (zh) 基于线性组合过渡算法的叶片修复模型重建方法
US8805908B2 (en) Approximation processing method and approximation processing device
WO2021203711A1 (zh) 一种基于几何重建模型的等几何分析方法
CN108599170B (zh) 一种适用于交直流系统的潮流获取方法
US7142206B1 (en) Shared N-patch edges
JP2015001789A (ja) 曲面生成装置、曲面生成プログラム及び曲面生成方法
Li et al. An alternative method for constructing interpolatory subdivision from approximating subdivision
CN116720268A (zh) 一种周期性描述的叶片叶型全局光顺重建方法
WO2010097076A1 (de) Verfahren zum automatisierten erzeugen eines flächenmodells eines objekts
Liu et al. Optimal design of preform shape based on EFA-FEM-GA integrated methodology
CN109308382B (zh) 一种针对奇异点附近接触应力的分析方法
CN114626157A (zh) 一种基于数字化模型的实体轮廓曲线抽取方法
CN113658077A (zh) 一种基于曲率的分区域双边海量点云降噪方法
CN107798732B (zh) 一种自由曲面形态控制方法
Ma et al. Subdivision surface fitting from a dense triangle mesh
CN107330977B (zh) 一种基于Shepard插值的根号三细分飞艇囊体曲面重构方法
Ahmad et al. Curvature comparison of bézier curve, ball curve and trigonometric curve in preserving the positivity of real data
Wu et al. Segmentation and Iterative Algorithm for Repair of Turbine Blade
CN112698623B (zh) 一种用于多轴轮廓应用场合中的多轴轮廓控制方法
Pérez Rubio et al. Reconstruction of a Steam Turbine Blade Using Piecewise Bernstein Polynomials and Transfinite Interpolation
Chen et al. A high-precision T-spline blade modeling method with 2D and 3D knot adaptive refinement
CN117850345A (zh) 基于t样条骨骼线的曲面分区域刀轨生成方法
CN117934766A (zh) 一种样条曲面的网格重参数化与质量评估方法
Sun et al. 3D B-Spline surface reconstruction based on the contour information

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