CN106294977B - 一种机器人铣削加工中工件装夹位置优划方法 - Google Patents

一种机器人铣削加工中工件装夹位置优划方法 Download PDF

Info

Publication number
CN106294977B
CN106294977B CN201610640732.6A CN201610640732A CN106294977B CN 106294977 B CN106294977 B CN 106294977B CN 201610640732 A CN201610640732 A CN 201610640732A CN 106294977 B CN106294977 B CN 106294977B
Authority
CN
China
Prior art keywords
milling
robot
workpiece
flutter
stability
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
Application number
CN201610640732.6A
Other languages
English (en)
Other versions
CN106294977A (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.)
University of Shanghai for Science and Technology
Original Assignee
University of Shanghai for Science and 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 University of Shanghai for Science and Technology filed Critical University of Shanghai for Science and Technology
Priority to CN201610640732.6A priority Critical patent/CN106294977B/zh
Publication of CN106294977A publication Critical patent/CN106294977A/zh
Application granted granted Critical
Publication of CN106294977B publication Critical patent/CN106294977B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Numerical Control (AREA)
  • Manipulator (AREA)

Abstract

本发明属于机器人铣削工艺优化技术领域,具体涉及一种机器人铣削加工中工件装夹位置优划方法。在末端装有电主轴的六轴工业机器人和固定的工作平台基础上,根据待加工工件及其夹具大小将工作平台适当地分为等面积的若干矩形区域作为备选装夹位置,并将夹具分别放置在这些备选区域的中心点位置。在每个区域内,分别对机器人铣削***进行模态实验,获取相应的模态参数,通过再生型与耦合型交叉颤振稳定性解析模型分析,计算得到颤振稳定性叶瓣图。选择稳定性叶瓣图中纵坐标最高的曲线对应的位置作为工件在工作台上的最佳安装位置。所述方法有效降低了机器人铣削过程中的颤振程度,提高了加工稳定性,从而保证加工零件质量与尺寸精度。

Description

一种机器人铣削加工中工件装夹位置优划方法
技术领域
本发明属于机器人铣削工艺优化技术领域,具体涉及一种机器人铣削加工中工件装夹位置优划方法。
背景技术
近年来,越来越多的工业机器人铣削***应用在汽车、航空制造等领域中塑料、铝制零件的加工生产中。相对于传统数控加工机床而言,使用工业机器人进行机械加工,具有成本低、自动化程度高、柔性好、易于集成传感器、外部驱动器等***设备的优点。尤其在大型复杂不规则零件的加工中,工业机器人切削***相对于数控机床,其占地面积与可加工工件规格相比明显减小;另外工业机器人能方便地实现切削加工工具在空间的各种位姿,避免刀具与零件复杂结构细节产生干涉。因此近年来,为适应多品种、小批量、现场加工的现代生产要求,通过在工业机器人末端安装电主轴和切削刀具,形成工业机器人切削加工***,使其具备类似于多轴数控加工机床的切削加工性能,在某些领域替代数控机床进行切削加工,已经成为工业机器人应用新的发展方向之一。
但在实际加工生产中,工业机器人由于刚度较差,在切削过程中,刀具会与被切削表面之间产生相对振动,这种切削振动会产生工件表面质量下降、刀具使用寿命减少、***可靠性降低等问题。在机器人铣削加工中,工件的初始工作位置对于颤振发生影响较大,当工件与机器人相对位置不同,机器人铣削时会产生不同振型。因此提出一种基于颤振稳定域分析的机器人铣削加工中工件装夹位置的优化方法,对提高机器人铣削加工稳定性具有非常重要的意义。
发明内容
本发明目的在于提供一种机器人铣削加工中工件装夹位置优划方法,通过对工件在工作台上装夹位置的合理规划,提高由工作平台和末端安装电主轴的六轴工业机器人组成的铣削加工***的切削稳定性。该方法可以有效降低机器人铣削过程中的颤振程度,保证良好的加工质量与精度。
为了达到上述目的,本发明的构思是:
首先在末端装有电主轴的六轴工业机器人和固定的工作平台基础上,根据待加工工件及其夹具大小将工作平台适当地分为等面积的若干矩形区域作为备选装夹位置,并将夹具分别放置在这些备选区域的中心点位置。在每个区域内,分别对机器人铣削***进行模态实验,获取相应的模态参数,通过再生型与耦合型交叉颤振稳定性解析模型分析,计算得到颤振稳定性叶瓣图。选择稳定性叶瓣图中纵坐标最高的曲线对应的位置作为工件在工作台上的最佳安装位置。
根据上述构思,本发明采用以下技术方案:
一种机器人铣削加工中工件装夹位置优划方法,基于颤振稳定域分析方法,包括如下步骤:
步骤1:根据待加工工件及其夹具与工作台的相对尺寸大小,确定工作台的装夹区域,从而将工作台分为等面积的若干备选装夹区域,并将夹具分别放置在这些备选区域的中心点位置;
步骤2:在每个区域内,利用力锤对机器人铣削***中的铣刀进行激励,获取铣刀在力锤激励下的加速度信号;
步骤3:根据力锤的激振力和铣刀在力锤激励下的加速度信号获得铣刀刀尖点处的位移频响函数,从位移频响函数获取机器人铣削***的模态参数,包括:模态频率、模态阻尼、模态质量;
步骤4:通过查询相关铣削参数表,确定铣削***中切削力参数及刀具参数;
步骤5:计算得到颤振稳定性叶瓣图。
利用频域法求解颤振稳定性解析模型,考虑交叉传递函数φxy和φyx,得到闭环动态铣削***的特征方程的特征值Λ:
Figure BDA0001071827500000021
式中:
a0=Φxx(iωcyy(iωc)(axxayy-axyayx)-
Φxy(iωcyx(iωc)(axxayy-axyayx)
a1=axxΦxx(iωc)+ayyΦyy(iωc)+axyΦyx(iωc)+ayxΦxy(iωc)
式中:Φxx(iωc)和Φyy(iωc)分别为X和Y方向上的直接传递函数;
Φxy(iωc)和Φyx(iωc)分别为X和Y方向上的交叉传递函数;
axx,ayy,axy和ayx为平均方向系数。
利用颤振理论中的频域解析法得到临界轴向切深为:
Figure BDA0001071827500000031
主轴转速n通过求得刀齿切削周期T得到:
Figure BDA0001071827500000032
式中,Ktc为径向切削力系数,ΛR和ΛI分别为***传递函数特征根的实部和虚部,k为叶瓣图中的叶瓣数,N为刀齿数;
利用轴向切深与主轴转速的计算公式,绘制出六个区域内的颤振稳定性叶瓣图;
步骤6:利用得到的稳定性叶瓣图,选择稳定性叶瓣图中纵坐标最高的曲线对应的位置作为待加工工件在工作台上的最佳安装位置。
与现有技术相比,本发明的有益效果是:
本方法得到机器人铣削加工时工件在工作台上的最佳安装位置,有效降低了机器人铣削过程中的颤振程度,从而大大降低机器人铣削加工过程中发生的颤振危害,提高了加工稳定性,从而保证加工零件质量与尺寸精度。
附图说明
图1为本发明的机器人铣削***平台示意图。
图2为本发明方法的流程示意图。
图3为本发明的工件在工作台上的备选装夹区域划分示意图。
图4为本发明的三个备选装夹位置下机器人铣削***的稳定性叶瓣图。
具体实施方式
以下结合附图详细叙述本发明的实施例。
如图1所示,机器人铣削颤振分析平台由机器人1、电主轴2、加速度传感器3、力锤4、振动采集仪5、计算机6、待加工工件7、夹具8、工作台9组成,所述电主轴2安装在机器人1的末端执行器上,所述加速度传感器3安装在电主轴2的末端,所述加速度传感器3、力锤4和计算机6分别与振动采集仪5连接。
如图2所示,一种机器人铣削加工中工件装夹位置优划方法,基于颤振稳定域分析方法,包括如下步骤:
第一步:如图3所示,在工作台9与六轴工业机器人1相对位置固定的情况下,工作台9尺寸为990mm×900mm,夹具8尺寸为300mm×400mm,待加工工件7尺寸为350mm×200mm。在工作台9可装夹区域内,将工作台9分为330mm×450mm的六个等面积的初始区域,作为备选待加工工件7的装夹位置。
第二步:将夹具8分别放置于初始划分的六个区域的中心点位置,如图3所示,在每个区域内,当铣削加工时,利用力锤4对于机器人1、电主轴2及刀具***进行X方向锤击,X方向采集加速度信号,Y方向锤击,Y方向采集加速度信号,X方向锤击,Y方向采集加速度信号,Y方向锤击,X方向采集加速度信号的模态试验。
第三步:利用加速度传感器3的信号与力锤4的激励信号,分别分析加速度频响函数,获得相应的模态频率、模态阻尼以及模态质量参数。
第四步:确定铣削***中切削力参数及刀具参数。通过查典型材料切削力系数表,确定相应材料切削力系数Ktc=800N/mm2。根据刀具型号确定刀齿数N=4。
第五步:计算轴向切深与主轴转速表达式,构建直观的稳定叶瓣图。利用
Figure BDA0001071827500000041
Figure BDA0001071827500000042
公式,计算得到六个区域内的稳定叶瓣图。
第六步:观察叶瓣图曲线,选择曲线纵坐标幅值最大的曲线所对应的装夹位置作为机器人铣削加工时工件最优位置。图4所示的是六个位置的叶瓣图曲线中轴向切深幅值最大、中间及最小的三条曲线。在稳定性叶瓣图中,曲线下方为稳定加工区域,曲线上方为颤振区。由此在六个备选装夹区域内,在位置3机器人铣削***的稳定加工区面积最大,位置5面积最小,位置2面积大小为中间值。所以位置3是这个机器人铣削***的最佳工件装夹位置。
所述的基于颤振稳定域分析的机器人铣削加工中工件装夹位置的优化方法,有效降低了机器人铣削加工中的颤振程度,从而保证了加工零件质量与尺寸精度。在实际应用中,本方法可以拓展至机器人钻削、磨削等***的颤振抑制中。

Claims (1)

1.一种机器人铣削加工中工件装夹位置优划方法,其特征在于,基于颤振稳定域分析方法,包括如下步骤:
步骤1:根据待加工工件(7)及其夹具(8)与工作台(9)的相对尺寸大小,确定工作台(9)的装夹区域,从而将工作台(9)分为等面积的若干备选装夹区域,并将夹具(8)分别放置在这些备选区域的中心点位置;
步骤2:在每个区域内,利用力锤(4)对机器人铣削***中的铣刀进行激励,获取铣刀在力锤(4)激励下的加速度信号;
步骤3:根据力锤(4)的激振力和铣刀在力锤(4)激励下的加速度信号获得铣刀刀尖点处的位移频响函数,从位移频响函数获取机器人铣削***的模态参数,包括:模态频率、模态阻尼、模态质量;
步骤4:通过查询相关铣削参数表,确定铣削***中切削力参数及刀具参数;
步骤5:计算得到颤振稳定性叶瓣图,
利用频域法求解颤振稳定性解析模型,考虑交叉传递函数φxy和φyx,得到闭环动态铣削***的特征方程的特征值Λ:
Figure FDA0002208343240000011
式中:
a0=Φxx(iωcyy(iωc)(axxayy-axyayx)-Φxy(iωcyx(iωc)(axxayy-axyayx)
a1=axxΦxx(iωc)+ayyΦyy(iωc)+axyΦyx(iωc)+ayxΦxy(iωc)
式中:Φxx(iωc)和Φyy(iωc)分别为X和Y方向上的直接传递函数;
Φxy(iωc)和Φyx(iωc)分别为X和Y方向上的交叉传递函数;
axx,ayy,axy和ayx为平均方向系数;
利用颤振理论中的频域解析法得到临界轴向切深为:
Figure FDA0002208343240000021
主轴转速n通过求得刀齿切削周期T得到:
Figure FDA0002208343240000022
式中,Ktc为径向切削力系数,ΛR和ΛI分别为***传递函数特征根的实部和虚部,k为叶瓣图中的叶瓣数,N为刀齿数;
利用轴向切深与主轴转速的计算公式,绘制出六个区域内的颤振稳定性叶瓣图;
步骤6:利用得到的稳定性叶瓣图,选择稳定性叶瓣图中纵坐标最高的曲线对应的位置作为待加工工件(7)在工作台(9)上的最佳安装位置。
CN201610640732.6A 2016-08-08 2016-08-08 一种机器人铣削加工中工件装夹位置优划方法 Active CN106294977B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610640732.6A CN106294977B (zh) 2016-08-08 2016-08-08 一种机器人铣削加工中工件装夹位置优划方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610640732.6A CN106294977B (zh) 2016-08-08 2016-08-08 一种机器人铣削加工中工件装夹位置优划方法

Publications (2)

Publication Number Publication Date
CN106294977A CN106294977A (zh) 2017-01-04
CN106294977B true CN106294977B (zh) 2020-01-21

Family

ID=57665901

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610640732.6A Active CN106294977B (zh) 2016-08-08 2016-08-08 一种机器人铣削加工中工件装夹位置优划方法

Country Status (1)

Country Link
CN (1) CN106294977B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107609311B (zh) * 2017-10-17 2020-02-18 西北工业大学 基于排屑力模型的深孔钻削深度优化方法
CN108519758B (zh) * 2018-03-21 2020-12-22 南京工业大学 一种避免干涉和颤振的多轴数控加工球头刀轴矢量规划方法
CN108638076B (zh) * 2018-07-25 2021-04-13 东北大学 一种六自由度串联机器人铣削加工三维稳定性预测方法
CN109093650B (zh) * 2018-09-25 2021-07-13 上海大学 一种机器人动态特性测定方法及***、装置
CN109375578B (zh) * 2018-10-31 2020-06-16 湖南工学院 一种机油泵壳体内深孔高效加工控制方法
CN117103280B (zh) * 2023-10-19 2023-12-22 中国长江电力股份有限公司 一种大型水轮机顶盖在位机器人减材加工方法及***

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102554326B (zh) * 2012-02-17 2014-02-26 上海交通大学 一种基于叶轮叶片动态刚度的铣削精加工方法
CN103823945A (zh) * 2014-03-13 2014-05-28 大连理工大学 一种平面切削过程的颤振稳定域建模方法
CN105678043A (zh) * 2014-11-17 2016-06-15 天津职业技术师范大学 一种考虑刚度时变的大切除率铣削颤振监测方法
CN104484506B (zh) * 2014-11-25 2017-05-24 东北大学 一种基于可靠性叶瓣图的车削加工颤振预测方法

Also Published As

Publication number Publication date
CN106294977A (zh) 2017-01-04

Similar Documents

Publication Publication Date Title
CN106294977B (zh) 一种机器人铣削加工中工件装夹位置优划方法
Li et al. Effect of the cutter path and the workpiece clamping position on the stability of the robotic milling system
Mejri et al. Dynamic characterization of machining robot and stability analysis
Scippa et al. FEM based cutting velocity selection for thin walled part machining
CN111390250B (zh) 弱刚度薄壁结构件及其加工方法和工位快换定位装夹装置
JP2005074568A (ja) 多軸加工機、ワークの加工方法
CN111002313B (zh) 一种机器人模态参数辨识与动态特性分析的方法
CN108339995B (zh) 一种可调角度的机器人末端执行器安装机构及其调整方法
CN113182567B (zh) 一种数控铣削刀具快速试切及切削参数优选方法
CN112276571B (zh) 一种斜孔加工方法
Qin et al. A novel approach for the acquisition of vibration signals of the end effector in robotic drilling
CN101859124B (zh) 数控龙门铣床五面体程序控制方法
WO2022149322A1 (ja) 工作機械、工具にかかる力の推定方法、および工具にかかる力の推定プログラム
Adetoro et al. Stability lobes prediction in thin wall machining
Grossi et al. A novel experimental-numerical approach to modeling machine tool dynamics for chatter stability prediction
JP5855840B2 (ja) 2主軸対向nc旋盤及び防振加工方法
Huynh et al. Report of robotic machining measurements using a stäubli tx200 robot: application to milling
Jayaweera et al. Robotic edge profiling of complex components
CN110877188A (zh) 一种飞机盒体零件的加工方法
CN206393209U (zh) 车铣一体模组
Yáñez-Valdez et al. Dynamic stability of a parallel kinematic machine
Zhang et al. A novel in-process machining deformation perception and control method
Ninevski et al. Instrumented Milling Process and Analysis for Tool Wear Measurement
Jayaweera et al. Measurement assisted automated robotic edge deburring of complex components
CN209793504U (zh) 多轴加工通用虎钳

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant