CN109773593B - 一种基于余量约束条件下的磨削方法 - Google Patents
一种基于余量约束条件下的磨削方法 Download PDFInfo
- Publication number
- CN109773593B CN109773593B CN201910070955.7A CN201910070955A CN109773593B CN 109773593 B CN109773593 B CN 109773593B CN 201910070955 A CN201910070955 A CN 201910070955A CN 109773593 B CN109773593 B CN 109773593B
- Authority
- CN
- China
- Prior art keywords
- model
- allowance
- grinding
- theoretical model
- coordinate system
- 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
Landscapes
- Numerical Control (AREA)
- Constituent Portions Of Griding Lathes, Driving, Sensing And Control (AREA)
- Grinding And Polishing Of Tertiary Curved Surfaces And Surfaces With Complex Shapes (AREA)
Abstract
本发明公开了一种基于余量约束条件下的磨削方法,其包括如下步骤:步骤一,获取零部件的三维测量数据,根据三维测量数据逆向重构零部件的重构模型;步骤二,初始匹配,得到重构模型到理论模型的初始变换矩阵;步骤三,精确匹配,通过迭代运算完成重构模型和理论模型的精确匹配,得到零部件的加工基准和加工余量;步骤四,由步骤三得到的加工基准和加工余量设定磨削参数,对零部件进行磨削加工。其能够有效确定零部件的加工基准和加工余量,有效解决无基准叶片类零部件的数控砂带自适应磨削工艺难的问题。
Description
技术领域
本发明涉及叶片类零件的磨削加工,具体涉及一种基于余量约束条件下的磨削方法。
附图说明
无基准叶片类零件其叶身型面及榫头缘板是一般未经过机械加工,故在实际数控砂带自适应磨削过程中,无法以其榫头等部位确定加工基准。此类无基准叶片类零件后续磨削加工,此前一直是通过人力的方式进行加工,而人力加工的方式不仅效率低下,一致性差,良品率低,无法满足此类零件的加工需求,而且劳动环境恶劣,劳动强度大,不利于从业人员的身体健康。近年来,数控砂带自适应磨削加工此类无基准叶片零件是发展的趋势,而在自适应磨削过程中,加工基准的确定是影响此类零件加工精度的重要因素。同时,零部件的加工余量是磨抛加工的重要依据,通过在零部件的理论模型上规划磨抛加工轨迹,获得刀触点坐标,求取零部件理论模型上待加工区域的刀触点沿法失方向到重构模型的距离,即为零部件的型面的加工余量,因此,重构模型与理论模型的匹配精度将直接影响加工基准和加工余量的确定。
发明内容
本发明的目的是提供一种基于余量约束条件下的磨削方法,其能够有效确定零部件的加工基准和加工余量,有效解决无基准叶片类零部件的数控砂带自适应磨削工艺难的问题。
本发明所述的基于余量约束条件下的磨削方法,其包括如下步骤:
步骤一,获取零部件的三维测量数据,根据三维测量数据逆向重构零部件的重构模型,重构模型设为曲面Q(u1,v1),理论模型设为曲面p(u,v);
步骤二,初始匹配,在重构模型上选取三个数据点m1、m2和m3,理论模型上选取三个数据点n1、n2和n3,通过m1、m2和m3在重构模型上建立以m1为坐标原点的坐标系通过n1、n2和n3在理论模型上建立以n1为坐标原点的坐标系得到重构模型到理论模型的平移向量为:
式中,ω1为三个数据点m1、m2和m3的形心,ω1=(m1+m2+m3)/3,ω2为三个数据点n1、n2和n3的形心,ω2=(n1+n2+n3)/3,R为重构模型到理论模型的初始变换矩阵;
步骤三,精确匹配,通过迭代运算完成重构模型和理论模型的精确匹配,得到零部件的加工基准和加工余量;
步骤四,由步骤三得到的加工基准和加工余量设定磨削参数,对零部件进行磨削加工。
进一步,所述步骤二中的坐标系∑1的基矢为:
坐标系∑2的基矢为:
进一步,所述步骤二中的初始变换矩阵R的计算公式为:
进一步,所述步骤三中的迭代运算为:首先建立目标函数E(g),设定迭代终止精度εe;
通过求解非线性方程,得到第k次迭代的刚体变换矩阵gk:
式中,α、β和γ为迭代前、后重构模型坐标系∑1的基矢之间的夹角,Δx、Δy和Δz为迭代前、后重构模型坐标系∑1的基矢平移量;
将计算结果代入目标函数,计算迭代收敛准则εiter:
若εiter≤εe,则结束迭代运算,若εiter>εe,则增加迭代次数,继续迭代运算。
本发明与现有技术相比具有如下有益效果:
1、本发明采用初始匹配和精确匹配相结合的方法来确定无基准零部件的加工基准和加工余量,保证了无基准零部件的磨削加工。
2、本发采用“三点法”实现重构模型和理论模型的初始匹配,便于操作,可行性高,提高了初始匹配的精度。
3、本发明通过建立特定的优化目标函数,保证了精确匹配的精度,进而更好地确定了零部件的加工基准和加工余量,保证了零部件的磨削精度。
附图说明
图1是本发明的流程示意图。
具体实施方式
下面结合附图和具体实施例对本发明作详细说明。
船用整体螺旋桨成型工艺一般是采用毛坯浇注,没有进行机械加工,故无加工基准,同时由于船用螺旋桨尺寸大、体积重量大,且需要多道工序完成加工,本实施例选用机器人夹持磨头完成船用螺旋桨的后续自适应磨抛加工,参见图1,所示的基于余量约束条件下的磨削方法,其包括如下步骤。
步骤一,选用海克斯康GLOBAL SILVER CLASSIC SR 05.07.05型号接触式三坐标测量机对其进行检测。根据检测获取的螺旋桨点云数据,利用NX8.5三维造型软件,采用三次非均匀有理B样条曲线曲面算法对螺旋桨叶身型面进行模型重构,重构模型设为曲面Q(u1,v1),理论模型设为曲面p(u,v)。
步骤二,将船用螺旋桨的重构模型与理论模型进行初始匹配,在重构模型上选取三个数据点m1、m2和m3,理论模型上选取三个数据点n1、n2和n3,通过m1、m2和m3在重构模型上建立以m1为坐标原点的坐标系通过n1、n2和n3在理论模型上建立以n1为坐标原点的坐标系所述坐标系∑1的基矢为:
坐标系∑2的基矢为:
由于两个初始重构模型和理论模型的形心应该重合,得到重构模型到理论模型的平移向量为:式中,ω1为三个数据点m1、m2和m3的形心,ω1=(m1+m2+m3)/3,ω2为三个数据点n1、n2和n3的形心,ω2=(n1+n2+n3)/3,R为重构模型到理论模型的初始变换矩阵。
步骤三,精确匹配,建立非线性最小二乘目标函数E(g),设定迭代终止精度εe;所述目标函数为:
为了求得ming,采用迭代运算来加快全局最优收敛速度,设经过初始匹配得到的变换矩阵为g0,g0=R,当迭代次数k=0,计算目标函数
通过求解非线性方程,得到第k次迭代的刚体变换矩阵gk:
式中,α、β和γ为迭代前、后重构模型坐标系∑1的基矢之间的夹角,Δx、Δy和Δz为迭代前、后重构模型坐标系∑1的基矢平移量;
将计算结果代入目标函数,计算迭代收敛准则εiter:
若εiter≤εe,则结束迭代运算,实现了重构模型与理论模型的精确匹配,从而确定船用螺旋桨的加工基准和加工余量。
若εiter>εe,则增加迭代次数,继续迭代运算。
步骤四,由步骤三得到的加工基准和加工余量分布信息,生成机器人自适应磨削加工程序,对船用螺旋桨型面局部超差区域进行磨削加工,再进行精磨抛光,完成船用螺旋桨的磨抛加工。
Claims (3)
1.一种基于余量约束条件下的磨削方法,其特征在于包括如下步骤:
步骤一,获取零部件的三维测量数据,根据三维测量数据逆向重构零部件的重构模型,重构模型设为曲面Q(u1,v1),理论模型设为曲面p(u,v);
步骤二,初始匹配,在重构模型上选取三个数据点m1、m2和m3,理论模型上选取三个数据点n1、n2和n3,通过m1、m2和m3在重构模型上建立以m1为坐标原点的坐标系通过n1、n2和n3在理论模型上建立以n1为坐标原点的坐标系得到重构模型到理论模型的平移向量为:
式中,ω1为三个数据点m1、m2和m3的形心,ω1=(m1+m2+m3)/3,ω2为三个数据点n1、n2和n3的形心,ω2=(n1+n2+n3)/3,R为重构模型到理论模型的初始变换矩阵;
步骤三,精确匹配,通过迭代运算完成重构模型和理论模型的精确匹配,得到零部件的加工基准和加工余量;
所述迭代运算为:首先建立目标函数E(g),设定迭代终止精度εe;
通过求解非线性方程,得到第k次迭代的刚体变换矩阵gk:
式中,α、β和γ为迭代前、后重构模型坐标系∑1的基矢之间的夹角,Δx、Δy和Δz为迭代前、后重构模型坐标系∑1的基矢平移量;
将计算结果代入目标函数,计算迭代收敛准则εiter:
若εiter≤εe,则结束迭代运算,若εiter>εe,则增加迭代次数,继续迭代运算;
步骤四,由步骤三得到的加工基准和加工余量设定磨削参数,对零部件进行磨削加工。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910070955.7A CN109773593B (zh) | 2019-01-25 | 2019-01-25 | 一种基于余量约束条件下的磨削方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910070955.7A CN109773593B (zh) | 2019-01-25 | 2019-01-25 | 一种基于余量约束条件下的磨削方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109773593A CN109773593A (zh) | 2019-05-21 |
CN109773593B true CN109773593B (zh) | 2021-01-22 |
Family
ID=66502451
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910070955.7A Active CN109773593B (zh) | 2019-01-25 | 2019-01-25 | 一种基于余量约束条件下的磨削方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109773593B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111300208A (zh) * | 2019-12-19 | 2020-06-19 | 重庆大学 | 一种复杂曲面机器人砂带变进给自适应磨削方法及设备 |
CN111618654B (zh) * | 2020-06-04 | 2022-07-08 | 清华大学 | 一种叶片加工工序的自适应调控方法、装置、及电子设备 |
CN115625595B (zh) * | 2022-12-21 | 2023-03-17 | 太原理工大学 | 一种定量控制的叶片抛磨方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102183205A (zh) * | 2011-01-19 | 2011-09-14 | 北京航空航天大学 | 一种大型零部件最佳装配位姿匹配的方法 |
CN104616311A (zh) * | 2015-02-13 | 2015-05-13 | 重庆大学 | 基于改进icp算法的损伤零部件精确配准方法 |
CN106407683B (zh) * | 2016-09-19 | 2019-01-15 | 上海理工大学 | 基于磨削去除率模型的切入磨削工艺参数优化方法 |
CN108279643B (zh) * | 2018-01-29 | 2020-08-07 | 西南交通大学 | 一种基于测量点与自适应差分进化算法的工件姿态调整方法 |
CN108917600A (zh) * | 2018-07-05 | 2018-11-30 | 西安交通大学 | 一种基于三维反求的复杂结构毛坯加工定位方法 |
CN109117560B (zh) * | 2018-08-17 | 2020-03-31 | 武汉理工大学 | 基于mbd的汽车典型机加工零件三维工艺设计方法及平台 |
-
2019
- 2019-01-25 CN CN201910070955.7A patent/CN109773593B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109773593A (zh) | 2019-05-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Calleja et al. | Highly accurate 5-axis flank CNC machining with conical tools | |
CN109773593B (zh) | 一种基于余量约束条件下的磨削方法 | |
CN108279643B (zh) | 一种基于测量点与自适应差分进化算法的工件姿态调整方法 | |
CN110434671A (zh) | 一种基于特征测量的铸造件表面机加工轨迹校准方法 | |
CN110069041B (zh) | 一种基于在机测量的工件加工方法及系统 | |
CN101000499A (zh) | 基于多传感器集成测量的仿形加工方法及系统 | |
CN110103071B (zh) | 一种变形复杂零件的数字化寻位加工方法 | |
CN102357666B (zh) | 自由曲面叶片平底刀三坐标端铣加工方法 | |
Warkentin et al. | Computer aided 5-axis machining | |
CN113276130B (zh) | 一种基于点云切片的自由曲面喷涂路径规划方法及系统 | |
Chong et al. | A trajectory planning method on error compensation of residual height for aero-engine blades of robotic belt grinding | |
CN111830900B (zh) | 一种机器人磨抛加工整体叶盘的无干涉刀具轨迹生成方法 | |
Zhang et al. | 3D curvature grinding path planning based on point cloud data | |
CN105643395A (zh) | 一种光学自由曲面磨削成型方法 | |
WO2022021479A1 (zh) | 一种基于 knn-icp 算法的自由曲面定位方法 | |
Xiao et al. | A novel trajectory planning method based on reverse compensation of profile error for robotic belt grinding of blisk | |
CN104392030B (zh) | 一种基于stl三维模型的曲面加工方法 | |
CN110227994B (zh) | 一种高阶切触包络式磨抛加工方法及其应用 | |
Hauth et al. | Constant cusp toolpath generation in configuration space based on offset curves | |
CN109035238B (zh) | 一种面向自由曲面零件的加工余量离线分析方法 | |
Chu et al. | Spline-constrained tool-path planning in five-axis flank machining of ruled surfaces | |
Feiyan et al. | A high-efficiency generation method of integral impeller channel tool path based on parametric domain template trajectory mapping | |
Zhang et al. | Robot automation grinding process for nuclear reactor coolant pump based on reverse engineering | |
Wu et al. | Rigid shape matching for 3-D robotic grinding measurement with applications to blades | |
CN110221575A (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 |