CN113868802B - 一种变壁厚约束下的空心叶片余量优化模型建立及求解方法 - Google Patents
一种变壁厚约束下的空心叶片余量优化模型建立及求解方法 Download PDFInfo
- Publication number
- CN113868802B CN113868802B CN202111169071.0A CN202111169071A CN113868802B CN 113868802 B CN113868802 B CN 113868802B CN 202111169071 A CN202111169071 A CN 202111169071A CN 113868802 B CN113868802 B CN 113868802B
- Authority
- CN
- China
- Prior art keywords
- blade
- curved surface
- allowance
- wall thickness
- processed
- 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
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P70/00—Climate change mitigation technologies in the production process for final industrial or consumer products
- Y02P70/10—Greenhouse gas [GHG] capture, material saving, heat recovery or other energy efficient measures, e.g. motor control, characterised by manufacturing processes, e.g. for rolling metal or metal working
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Electrical Discharge Machining, Electrochemical Machining, And Combined Machining (AREA)
- Numerical Control (AREA)
Abstract
本发明涉及一种变壁厚约束下的空心叶片余量优化模型建立及求解方法,属于制造加工领域。本发明结合空心叶片加工工艺特点,建立变壁厚公差约束、面轮廓度公差约束以及加工余量约束条件,以加工余量最均匀为优化目标,构建了新的余量优化建模方法及求解方案。本发明有效地优化了航空发动机空心叶片加工过程并一定程度上拯救报废零件,能够极大提高加工效率、降低加工成本,并且提高了航空发动机叶片精密制造的自适应程度。
Description
技术领域
本发明属于制造加工领域,涉及空心叶片余量优化建模及求解方法研究领域,具体地说是一种针对空心叶片变壁厚公差约束、面轮廓度公差约束以及加工余量约束下的余量优化模型建立及求解方法。
背景技术
自适应加工是当前制造领域的一个热门发展趋势,因为其依据待加工零件毛坯的实际状态自动地调整和改善加工方案的理念和特点,使其被广泛地运用到高端制造行业中去,成为提高复杂零部件加工质量、生产效率以及降低产品报废率的重要手段之一。余量优化技术是自适应加工中不可或缺的一环,其结合机内测量手段,可极大地提高产线的自适应程度。
航宇制造是制造业中的尖端领域,而航空发动机叶片更是航宇制造中的明珠。由于航空发动机叶片形状复杂、材料昂贵、生产周期长、加工精度要求高等特点,传统的余量优化方法已经不满足其高效、高精度生产的需求。特别是对于空心叶片,其毛坯受到超塑成形/扩散连接工艺水平的限制,往往变形较大,且空心结构复杂,传统余量优化方法仅以加工余量作为单一约束条件极易导致最终叶片部件报废。
为解决上述问题,提高航空发动机空心叶片制造的自适应程度,本发明结合空心叶片加工工艺特点,建立变壁厚公差约束、面轮廓度公差约束以及加工余量约束条件,以加工余量最均匀为优化目标,构建了新的余量优化建模方法及求解方案。
发明内容
要解决的技术问题
针对航空发动机空心叶片结构形状复杂,而传统的自适应加工余量优化过程仅考虑加工余量约束的优化结果极易导致空心叶片加工超差的问题,本发明旨在提供一种提高空心叶片余量优化结果质量的余量优化模型建立与求解方法。
技术方案
一种变壁厚约束下的空心叶片余量优化模型建立及求解方法,其特征在于步骤如下:
步骤1:测量空心叶片待加工面
将待加工空心叶片安装在数控机床工作台上,通过机内测量手段获取叶片待加工面的数据;所述的数据是叶片的待加工面外形离散测量数据点集以及空心叶片待加面区域的壁厚测量数据H={hi|i=1,2,…,NH},其中N和NH分别表示外形测量点和壁厚测量点的数量;
步骤2:构建空心叶片理论待加工面的余量优化模型
步骤2.1:参数化叶片理论待加工曲面
叶片待加工曲面的设计截面曲线族参数化表达为:
基于参数化的叶片设计截面曲线族,通过蒙面算法,可以进一步得到叶片理论待加工曲面的参数化表达:
其中u和v分别为叶片待加工曲面的U向和V向参数,u,v∈[0,1];di,j是参数曲面的控制顶点;m和n分别是参数曲面在U向和V向的控制顶点数量减1;Ni,3(u)和Nj,3(v)为参数曲面的基函数,3表示参数曲面为3次样条曲面;
步骤2.2:构建目标加工曲面
本发明中空心叶片的余量优化通过对叶片理论待加工曲面在空间中进行旋转和平移[R,T]实现,经过对理论待加工曲面在约束条件的约束下进行旋转、平移后得到目标加工曲面:
Sg(u,v)=Sd(u,v)·R·T (3)
其中R为绕叶片积叠轴的旋转矩阵,T为沿垂直叶片积叠轴平面的平移矩阵;
步骤2.3:构建3种约束条件
(1)加工余量约束条件F1:余量优化后的目标加工曲面Sg(u,v)必须位于叶片外形测量点集Pm的内侧,即处处都有加工余量;
(2)面轮廓度公差约束条件F2:余量优化后的目标加工曲面Sg(u,v)必须位于空心叶片待加工曲面设计的面轮廓度公差带之间;
(3)变壁厚公差约束条件F3:余量优化后的目标加工曲面Sg(u,v)必须满足叶片的设计壁厚要求,即Sg(u,v)必须位于空心叶片待加工曲面设计的壁厚公差带之间;
步骤2.4:构建优化目标函数
首先依据叶片外形测点集求解平均加工余量:
叶片余量优化的目标为加工余量最均匀:
步骤3:求解空心叶片余量优化模型
基于步骤2构建的3个约束条件和步骤3构建的优化目标函数,通过有效集方法active-set对余量优化模型进行求解。
步骤1中使用接触式测头和超声波壁厚测头对其待加工叶背曲面进行外形和壁厚测量。
有益效果
本发明提出的一种变壁厚约束下的空心叶片余量优化模型建立及求解方法,通过构建加工余量约束、面轮廓度公差约束以及变壁厚公差约束3个约束条件,有效地优化了航空发动机空心叶片加工过程并一定程度上拯救报废零件,能够极大提高加工效率、降低加工成本,并且提高了航空发动机叶片精密制造的自适应程度。
附图说明
附图仅用于示出具体实施例的目的,而并不认为是对本发明的限制,在整个附图中,相同的参考符号表示相同的部件。
图1说明示例叶片的叶背侧;
图2蒙面算法参数化叶片理论待加工曲面;
图3加工余量约束示意图;
图4面轮廓度公差约束示意图;
图5变壁厚公差约束示意图;
图6余量优化前加工余量分布;
图7余量优化后加工余量分布。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图和实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。此外,下面描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
以图1所示宽弦空心风扇叶片的叶背面为例,按照本发明内容中的步骤,参照附图2-7对本发明做出具体的说明,其具体实施方式如下:
步骤一:测量空心叶片待加工面
将航空发动机空心叶片安装在机床上,使用接触式测头和超声波壁厚测头对其待加工叶背曲面进行外形和壁厚测量,分别得到表1所示的外形点位数据和表2所示的壁厚点位数据。
表1叶片待加工曲面的外形测量数据
表2叶片待加工曲面的壁厚测量数据
步骤二:
步骤2.1:参数化叶片理论待加工曲面
如图2所示,将待加工曲面按照设计截面的位置取一系列的截面线族,设计截面曲线的参数化表达是已知的:
其中u表示设计截面曲线的U向参数,u∈[0,1];Ni,3(u)是参数曲线的基函数,3表示曲线的样条次数;为参数曲线的控制顶点;m是控制顶点数量减1;s表示设计截面线的数量,实施例中设计截面线数量为30。
然后通过蒙面算法,对叶片截面线族进行参数化曲面的拟合表达,得到实施例中叶片理论待加工曲面Sd(u,v)的参数化表达:
其中u和v分别为叶片待加工曲面的U向和V向参数,u,v∈[0,1];di,j是参数曲面的控制顶点;m和n分别是参数曲面在U向和V向的控制顶点数量减1;Ni,3(u)和Nj,3(v)为参数曲面的基函数,3表示参数曲面为3次样条曲面。实施例中叶片曲面是用双三次B样条曲面的形式进行参数化表达的。
步骤2.2:构建目标加工曲面
本发明中空心叶片的余量优化通过对叶片理论待加工曲面在空间中进行旋转和平移[R,T]实现,经过对理论待加工曲面在约束条件的约束下进行旋转、平移后得到目标加工曲面:
Sg(u,v)=Sd(u,v)·R·T (3)
其中R为绕叶片积叠轴的旋转矩阵,T为沿垂直叶片积叠轴平面的平移矩阵,在实施例中,依据待加工叶片的结构特点,设置如下旋转矩阵和平移矩阵:
T=[Δx,Δy,0,0]T (5)
α、Δx和Δy分别是Sd(u,v)绕积叠轴的旋转角度、沿X轴和Y轴的平移距离,在实施例中,α、Δx和Δy为余量优化的优化求解变量。
步骤2.3:构建3种约束条件
(1)加工余量约束条件F1:如图3所示,余量优化后的目标加工曲面Sg(u,v)必须位于叶片外形测量点集Pm的内侧,即处处都有加工余量。因此本发明中加工余量约束通过外形测量点集Pm中的点到Sg(u,v)上对应的最近点之间的有向距离大于0来表示,也就是点到点的矢量方向指向叶片的外侧;
(2)面轮廓度公差约束条件F2:如图4所示,余量优化后的目标加工曲面Sg(u,v)必须位于空心叶片待加工曲面设计的面轮廓度公差带之间。面轮廓度公差约束在本发明中定义成Sg(u,v)位于面轮廓度公差带上下两个公差带面之间。本发明为提升计算效率,采用离散方法来屁股Sg(u,v)是否位于叶片面轮廓度公差带内。首先在Sg(u,v)上取一定数量的离散点其中D表示离散点的数量。然后分别在上下公差带曲面上找到和对应的最近点和通过评估与之间的相对位置来评估Sg(u,v)与面轮廓度公差带之间的位置关系。即当Sg(u,v)满足面轮廓度公差约束时,预定有在Sg(u,v)外侧,在Sg(u,v)内侧:
(3)变壁厚公差约束条件F3:如图5所示,余量优化后的目标加工曲面Sg(u,v)必须满足叶片的设计壁厚要求,即Sg(u,v)必须位于空心叶片待加工曲面设计的壁厚公差带之间。图中和分别是壁厚测头在测量时的实际位置和测量法矢;为空心结构内表面上对应的点;和分别为叶片壁厚的上极限和下极限位置点;理想情况下,是从出发、沿方向的射线与Sg(u,v)的交点。但交点的理论解法复杂,本发明中的是Sg(u,v)上离的最近点。当Sg(u,v)满足变壁厚公差约束条件时,其必定在和之间。因此本发明中余量优化的变壁厚公差约束表达式为:
步骤2.4:构建优化目标函数
首先依据叶片外形测点集求解平均加工余量:
本发明叶片余量优化的目标为加工余量最均匀:
步骤三:求解空心叶片余量优化模型
基于步骤二构建的3个约束条件和步骤三构建的优化目标函数,通过有效集方法(active-set)对余量优化模型进行求解,最终解得优化变量如表3所示:
表3余量优化结果
经过余量优化后的目标加工曲面各处的加工余量分布如图7所示,可以看到加工余量分布情况较加工前分布更为均匀,且满足面轮廓度公差约束与变壁厚公差约束,有利于后续加工工艺优化。该实例证明了本发明的使用性和有效性。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明公开的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本发明的保护范围之内。
Claims (2)
1.一种变壁厚约束下的空心叶片余量优化模型建立及求解方法,其特征在于步骤如下:
步骤1:测量空心叶片待加工面
将待加工空心叶片安装在数控机床工作台上,通过机内测量手段获取叶片待加工面的数据;所述的数据是叶片的待加工面外形离散测量数据点集以及空心叶片待加面区域的壁厚测量数据H={hi|i=1,2,…,NH},其中N和NH分别表示外形测量点和壁厚测量点的数量;
步骤2:构建空心叶片理论待加工面的余量优化模型
步骤2.1:参数化叶片理论待加工曲面
叶片待加工曲面的设计截面曲线族参数化表达为:
基于参数化的叶片设计截面曲线族,通过蒙面算法,可以进一步得到叶片理论待加工曲面的参数化表达:
其中u和v分别为叶片待加工曲面的U向和V向参数,u,v∈[0,1];di,j是参数曲面的控制顶点;m和n分别是参数曲面在U向和V向的控制顶点数量减1;Ni,3(u)和Nj,3(v)为参数曲面的基函数,3表示参数曲面为3次样条曲面;
步骤2.2:构建目标加工曲面
本发明中空心叶片的余量优化通过对叶片理论待加工曲面在空间中进行旋转和平移[R,T]实现,经过对理论待加工曲面在约束条件的约束下进行旋转、平移后得到目标加工曲面:
Sg(u,v)=Sd(u,v)·R·T (3)
其中R为绕叶片积叠轴的旋转矩阵,T为沿垂直叶片积叠轴平面的平移矩阵;
步骤2.3:构建3种约束条件
(1)加工余量约束条件F1:余量优化后的目标加工曲面Sg(u,v)必须位于叶片外形测量点集Pm的内侧,即处处都有加工余量;
(2)面轮廓度公差约束条件F2:余量优化后的目标加工曲面Sg(u,v)必须位于空心叶片待加工曲面设计的面轮廓度公差带之间;
(3)变壁厚公差约束条件F3:余量优化后的目标加工曲面Sg(u,v)必须满足叶片的设计壁厚要求,即Sg(u,v)必须位于空心叶片待加工曲面设计的壁厚公差带之间;
步骤2.4:构建优化目标函数
首先依据叶片外形测点集求解平均加工余量:
叶片余量优化的目标为加工余量最均匀:
步骤3:求解空心叶片余量优化模型
基于步骤2构建的3个约束条件和步骤3构建的优化目标函数,通过有效集方法active-set对余量优化模型进行求解。
2.根据权利要求1所述的一种变壁厚约束下的空心叶片余量优化模型建立及求解方法,其特征在于步骤1中使用接触式测头和超声波壁厚测头对其待加工叶背曲面进行外形和壁厚测量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111169071.0A CN113868802B (zh) | 2021-10-04 | 2021-10-04 | 一种变壁厚约束下的空心叶片余量优化模型建立及求解方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111169071.0A CN113868802B (zh) | 2021-10-04 | 2021-10-04 | 一种变壁厚约束下的空心叶片余量优化模型建立及求解方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113868802A CN113868802A (zh) | 2021-12-31 |
CN113868802B true CN113868802B (zh) | 2022-09-06 |
Family
ID=79001664
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111169071.0A Active CN113868802B (zh) | 2021-10-04 | 2021-10-04 | 一种变壁厚约束下的空心叶片余量优化模型建立及求解方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113868802B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116991117B (zh) * | 2023-09-25 | 2024-01-05 | 南京航空航天大学 | 一种面向个性化零件加工的快速编程方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372291A (zh) * | 2016-08-29 | 2017-02-01 | 西北工业大学 | 一种公差约束下的叶片余量优化模型建立及求解方法 |
CN106897501A (zh) * | 2017-01-23 | 2017-06-27 | 西北工业大学 | 面向自适应加工中基于叶片类零件变形的定位优化方法 |
-
2021
- 2021-10-04 CN CN202111169071.0A patent/CN113868802B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106372291A (zh) * | 2016-08-29 | 2017-02-01 | 西北工业大学 | 一种公差约束下的叶片余量优化模型建立及求解方法 |
CN106897501A (zh) * | 2017-01-23 | 2017-06-27 | 西北工业大学 | 面向自适应加工中基于叶片类零件变形的定位优化方法 |
Non-Patent Citations (1)
Title |
---|
多目标约束的精锻叶片几何重构优化算法;冯亚洲等;《航空学报》;20180319(第07期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113868802A (zh) | 2021-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102982200B (zh) | 一种飞机框肋类钣金零件工艺模型设计方法 | |
CN105739440A (zh) | 一种宽弦空心风扇叶片的自适应加工方法 | |
CN103586737A (zh) | 一种叶片型面高精度数控铣加工的补偿方法 | |
CN113868802B (zh) | 一种变壁厚约束下的空心叶片余量优化模型建立及求解方法 | |
WO2019080539A1 (zh) | 一种复合材料风扇叶片金属加强边的制造方法 | |
CN107491034B (zh) | 一种基于协变场泛函的刀轴矢量优化方法 | |
US6340424B1 (en) | Manufacture of complexly shaped articles using an automated design technique | |
CN109693007B (zh) | 一种整体叶盘电解初成型加工轨迹参数的设计优化方法 | |
CN111266803A (zh) | 一种薄壁异型金属铸件的加工方法 | |
CN109590523B (zh) | 一种整体叶盘叶片扭转及弯曲变形的反向修正方法 | |
CN115146405B (zh) | 一种基于非刚性配准变形的薄壁零件模型重构方法 | |
CN114036668A (zh) | 一种基于离心叶轮中间面曲线的可变自由曲面叶片生成方法 | |
CN114692339A (zh) | 基于缘头补偿的整体叶盘叶片曲面重构方法 | |
CN109530768B (zh) | 一种叶片叶尖削边减薄区的加工方法 | |
CN115168986A (zh) | 一种叶片截面几何特征参数的获取方法 | |
CN112051803B (zh) | 一种基于空间平面法向量的小线段光顺方法 | |
CN114074170A (zh) | 一种航空发动机叶片精锻件终锻模的优化设计方法 | |
CN106001338A (zh) | 高温合金无余量叶片轧制方法 | |
Yilmaz et al. | A study of turbomachinery components machining and repairing methodologies | |
Khaimovich et al. | Practice in the application of specialized technologic process patterns of 5-axis machining operation of blade rings of blisks | |
CN116944395A (zh) | 一种非相切转接圆角精锻叶片及制造方法 | |
CN113029049B (zh) | 基于加权正负余量方差最小化算法的工件光学测量方法 | |
CN110421077B (zh) | 具有大角度薄壁扭曲窄深腔特征的构件的复合成形方法 | |
CN104438480A (zh) | 一种采用工程拼焊板实现复杂零件渐进成形的加工方法 | |
CN108363890B (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 |