CN103268430A - 基于机床刀具动刚度测量的铣削工艺参数优化方法 - Google Patents
基于机床刀具动刚度测量的铣削工艺参数优化方法 Download PDFInfo
- Publication number
- CN103268430A CN103268430A CN201310244572XA CN201310244572A CN103268430A CN 103268430 A CN103268430 A CN 103268430A CN 201310244572X A CN201310244572X A CN 201310244572XA CN 201310244572 A CN201310244572 A CN 201310244572A CN 103268430 A CN103268430 A CN 103268430A
- Authority
- CN
- China
- Prior art keywords
- milling
- lathe
- cutter
- tooling system
- dynamic stiffness
- 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
Links
Images
Landscapes
- Automatic Control Of Machine Tools (AREA)
Abstract
本发明涉及一种基于铣床刀具动刚度测量的铣削加工工艺参数优化方法,其步骤是:1.机床铣刀系统的动刚度测量计算,1)建立机床铣刀系统动力学模型,2)机床铣刀静刚度测量计算,3)机床铣刀系统固有频率和阻尼比测量计算;2.铣削加工工艺优化。本发明在测量计算铣床刀具系统的静刚度和动刚度的基础上,通过测量计算机床刀具在X和Y方向上的静刚度与动刚度,并结合切削过程中铣削力获得铣削加工过程中刀具的变形量,利用遗传算法优化的铣削加工工艺参数,该方法易操作,实用性强,对实际铣削加工提高加工精度和加工效率有重要意义。
Description
技术领域
本发明涉及一种铣削加工工艺参数优化方法,尤其是一种基于机床刀具动刚度测量的铣削工艺参数优化方法。
背景技术
数控铣削加工过程中,机床、刀具、工件共同组成一个切削过程动力学系统,切削力、加工变形、振动等动力学现象已成为决定数控铣削加工质量与加工效率的关键因素。因此,通过对切削过程中的动力学系统的测量计算来优化铣削加工工艺参数对提供加工质量与加工效率有重要意义。以往一般是根据磨削颤振的稳定性曲线,确定磨削加工中的工艺参数,以使磨削过程在不发生颤振的前提下,最高效率的完成磨削任务。然而,该方法在实际工程确定磨削颤振的稳定性曲线较为复杂,难以在被测对象安装各种测量装置,不适于机床操作者掌握。而实际铣削加工过程中铣削力具有一定的交变频率,因此,本发明专利提出一种基于机床刀具系统动刚度测量的外圆磨削工艺参数优化方法,该方法通过利用压力传感器和电涡流位移传感器测量计算机床刀具系统的静刚度和动刚度,并结合切削过程中铣削力获得铣削加工过程中刀具的变形量,利用遗传算法优化的铣削加工工艺参数。该方法易操作,有很强的实用性,对实际铣削提高加工精度和加工效率有重要意义。
发明内容
本发明是要提供一种基于铣床刀具动刚度测量计算的铣削加工工艺参数优化方法,来获得最小刀具切削变形的优化铣削加工工艺参数,以此有效地提高铣削加工效率和加工质量。
为实现上述目的,本发明的技术方案是:一种基于铣床刀具动刚度测量的铣削加工工艺参数优化方法,其步骤是:
一 .机床铣刀系统的动刚度测量计算
1)建立机床铣刀系统动力学模型
机床-刀具系统简化为X和Y相互垂直方向上的二自由度振动系统,如下式所示:
经拉普拉斯转换得到:
由上式(2)可知,机床-刀具系统的动刚度可表示为:
2)机床铣刀静刚度测量计算
利用实验方法求刀具系统静刚度,刀具通过刀柄固定在主轴上,使用压力传感器在刀具的刀尖位置上施加一定力,同时电涡流位移传感器测量出刀具的刀尖变形量,通过数据采集卡和计算机记录压力信号和电涡流位移传感器的位移信号,利用下式计算出刀具静刚度:
3)机床铣刀系统固有频率和阻尼比测量计算
使用阶跃响应法来计算刀具系统的固有频率和阻尼比,步骤如下:首先利用拉力传感器通过细绳对刀具施加一定的拉力作为系统输入信号;然后使用剪刀把细绳瞬间剪断,通过电涡流位移传感器采集刀具另一端的位移振动的输出信号,最后通过阶跃响应频响函数计算出刀具系统的固有频率和阻尼比;
通过上述静刚度测量和上述阶跃响应法可测量出机床-刀具系统X方向上的静刚度、固有频率和阻尼比,进而求出X方向的动刚度;用相同方法测量出机床-刀具系统Y方向上的静刚度、固有频率和阻尼比,求出Y方向的动刚度;
二 .铣削加工工艺优化
根据铣削加工原理,得到作用在整个铣刀上在X和Y方向上的瞬时切削力,如下式:
(5)
通过上述加工工艺参数可确定铣刀X方向和Y方向的铣削力值,结合机床-刀具系统的动刚度测量计算式(3),计算出铣削刀具在X方向和Y方向上的动态变形量,如式(6)所示:
(7)
正规化后,各铣削工艺参数作为个体的表现形式为:
(8)
本发明的有益效果:
该发明专利在测量计算铣床刀具系统的静刚度和动刚度的基础上,通过测量计算机床刀具在X和Y方向上的静刚度与动刚度,并结合切削过程中铣削力获得铣削加工过程中刀具的变形量,利用遗传算法优化的铣削加工工艺参数,该方法易操作,实用性强,对实际铣削加工提高加工精度和加工效率有重要意义。
附图说明
图1是铣削刀具系统动力学模型图;
图2是铣削刀具静刚度测量示意图;
图3是铣削刀具动刚度测量示意图;
图4 是阶跃响应的刀尖位移振动信号图。
具体实施方式
本发明的基于铣床刀具动刚度测量的铣削加工工艺参数优化方法,其步骤是:
一 机床铣刀系统的动刚度测量计算
1)机床铣刀系统动力学模型
如图1所示,机床-刀具系统可简化为X和Y相互垂直方向上的二自由度振动系统,如下式所示:
经拉普拉斯转换可得:
由上式(2)可知,机床-刀具系统的动刚度可表示为:
由上式(3)可知,分别计算机床-刀具系统X方向和Y方向的静刚度、固有频率和阻尼比,可计算出各方向的动刚度。因为X方向和Y方向的动刚度计算方法相同,所以本发明专利以X方向为例进行计算说明。
2)机床铣刀静刚度测量计算
如图2所示,利用实验方法求刀具系统静刚度,刀具3通过刀柄2固定在主轴1上,使用压力传感器4在刀具3的刀尖位置上施加一定力,同时电涡流位移传感器5测量出刀具3的刀尖变形量,通过数据采集卡6和计算机7记录压力信号和电涡流位移传感器的位移信号,利用下式计算出刀具静刚度:
3)机床铣刀系统固有频率和阻尼比测量计算
如图3所示,本发明专利使用阶跃响应法来计算刀具系统的固有频率和阻尼比,步骤如下:首先利用拉力传感器9通过细绳8对刀具3施加一定的拉力作为系统输入信号;然后使用剪刀把细绳瞬间剪断,通过电涡流位移传感器5采集刀具另一端的位移振动的输出信号,其测量结果如图4所示;最后通过阶跃响应频响函数计算出刀具系统的固有频率和阻尼比。
通过上述静刚度测量和上述阶跃响应法可测量出机床-刀具系统X方向上的静刚度、固有频率和阻尼比,进而求出X方向的动刚度;用相同方法也可以测量出机床-刀具系统Y方向上的静刚度、固有频率和阻尼比,求出Y方向的动刚度。
二 铣削加工工艺优化
根据铣削加工原理,可得到作用在整个铣刀上在X和Y方向上的瞬时切削力,如下式:
通过上述加工工艺参数可确定铣刀X方向和Y方向的铣削力大小,结合机床-刀具系统的动刚度测量计算式(3),可计算出在一定切削工艺参数的情况下的铣削刀具在X方向和Y方向上的动态变形量,如式(6)所示:
专利利用遗传算法对此铣削工艺参数优化,设主轴转速、切削宽度、工件进给速度、刀具切削深度及每齿切削量的约束条件为、、、、,对应于任一组参数,根据公式(6)可计算出铣削刀具的总变形,铣削材料去除率为,建立目标优化函数和:
正规化后,各铣削工艺参数作为个体的表现形式为:
应用案例
该案例中所使用机床为沈阳机床厂生产的三轴立式加工中心VMC1165B,刀具的直径为10mm,材质为硬质合金刚,刀的齿数为2。因为在机床XY平面上的X和Y方向的动刚度计算方法相同,本案例以X方向动刚度计算为例。
1)铣床刀具静刚度测量计算
按照图2的方法测量刀具的静刚度。利用公式(4)可计算其静刚度:
2)铣床刀具系统固有频率和阻尼比测量计算
3)铣床刀具动刚度计算结果
根据公式(7),可计算动刚度如下:
4)铣削工艺参数优化计算
设铣削各工艺参数的约束条件为:主轴转速1000r/min 8000r/min、切削宽度0.1mm 10mm、工件进给速度100mm/min 500mm/min、刀具切削深度0.1mm 5mm及每齿切削量0.01mm 0.1mm。建立目标优化函数:
正规化后,各铣削工艺参数作为个体的表现形式为:
Claims (1)
1.一种基于铣床刀具动刚度测量的铣削加工工艺参数优化方法,其特征在于,其步骤是:
一 .机床铣刀系统的动刚度测量计算
1)建立机床铣刀系统动力学模型
机床-刀具系统简化为X和Y相互垂直方向上的二自由度振动系统,如下式所示:
经拉普拉斯转换得到:
由上式(2)得知,机床-刀具系统的动刚度表示为:
2)机床铣刀静刚度测量计算
利用实验方法求刀具系统静刚度,刀具(3)通过刀柄(2)固定在主轴(1)上,使用压力传感器(4)在刀具(3)的刀尖位置上施加一定力,同时电涡流位移传感器(5)测量出刀具(3)的刀尖变形量,通过数据采集卡(6)和计算机(7)记录压力信号和电涡流位移传感器的位移信号,利用下式计算出刀具静刚度:
3)机床铣刀系统固有频率和阻尼比测量计算
使用阶跃响应法来计算刀具系统的固有频率和阻尼比,步骤如下:首先利用拉力传感器(9)通过细绳(8)对刀具(3)施加一定的拉力作为系统输入信号;然后使用剪刀把细绳瞬间剪断,通过电涡流位移传感器(5)采集刀具另一端的位移振动的输出信号,最后通过阶跃响应频响函数计算出刀具系统的固有频率和阻尼比;
通过上述静刚度测量和上述阶跃响应法可测量出机床-刀具系统X方向上的静刚度、固有频率和阻尼比,进而求出X方向的动刚度;用相同方法测量出机床-刀具系统Y方向上的静刚度、固有频率和阻尼比,求出Y方向的动刚度;
二. 铣削加工工艺优化
根据铣削加工原理,得到作用在整个铣刀上在X和Y方向上的瞬时切削力,如下式:
通过上述加工工艺参数确定铣刀X方向和Y方向的铣削力值,结合机床-刀具系统的动刚度测量计算式(3),计算出铣削刀具在X方向和Y方向上的动态变形量,如式(6)所示:
(7)
正规化后,各铣削工艺参数作为个体的表现形式为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310244572.XA CN103268430B (zh) | 2013-06-20 | 2013-06-20 | 基于机床刀具动刚度测量的铣削工艺参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310244572.XA CN103268430B (zh) | 2013-06-20 | 2013-06-20 | 基于机床刀具动刚度测量的铣削工艺参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103268430A true CN103268430A (zh) | 2013-08-28 |
CN103268430B CN103268430B (zh) | 2016-08-10 |
Family
ID=49012058
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310244572.XA Expired - Fee Related CN103268430B (zh) | 2013-06-20 | 2013-06-20 | 基于机床刀具动刚度测量的铣削工艺参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103268430B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104834271A (zh) * | 2015-06-03 | 2015-08-12 | 上海理工大学 | 基于动刚度评价的直接进给轴伺服参数优化方法 |
CN105699035A (zh) * | 2016-01-28 | 2016-06-22 | 哈尔滨理工大学 | 一种车削大螺距梯形外螺纹的振动测试方法 |
CN106002486A (zh) * | 2016-07-27 | 2016-10-12 | 大连理工大学 | 基于微铣床主传动系统功率的微铣削力测量方法 |
CN109406189A (zh) * | 2018-12-17 | 2019-03-01 | 武汉科技大学 | 一种五轴机床的刚度场半解析构建方法 |
CN110044625A (zh) * | 2019-04-30 | 2019-07-23 | 东华大学 | 一种飞轮轴承动刚度的测量方法 |
CN110103078A (zh) * | 2019-06-11 | 2019-08-09 | 贵州理工学院 | 一种基于普通车床的工艺系统的动刚度测定装置 |
CN114137907A (zh) * | 2021-11-30 | 2022-03-04 | 清华大学 | 数控机床振动幅值预测的建模方法、装置、设备及介质 |
CN114988679A (zh) * | 2022-08-08 | 2022-09-02 | 启东市云鹏玻璃机械有限公司 | 一种玻璃切割机床智能走刀控制方法 |
CN117991718A (zh) * | 2022-10-18 | 2024-05-07 | 力劲精密机械(昆山)有限公司 | 一种卧式加工中心多轴协同状态监测方法和系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102248209A (zh) * | 2011-07-01 | 2011-11-23 | 黑龙江科技学院 | 薄壁复杂曲面工件铣削加工时机床极限稳定工艺参数的确定方法 |
CN102592035A (zh) * | 2012-03-20 | 2012-07-18 | 北京航空航天大学 | 一种车铣复合切削加工表面粗糙度及表面形貌仿真预测方法 |
-
2013
- 2013-06-20 CN CN201310244572.XA patent/CN103268430B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102248209A (zh) * | 2011-07-01 | 2011-11-23 | 黑龙江科技学院 | 薄壁复杂曲面工件铣削加工时机床极限稳定工艺参数的确定方法 |
CN102592035A (zh) * | 2012-03-20 | 2012-07-18 | 北京航空航天大学 | 一种车铣复合切削加工表面粗糙度及表面形貌仿真预测方法 |
Non-Patent Citations (3)
Title |
---|
FRANC CUS,JOZE BALIC,JANEZ KOPAC: "determination and optimization of parameters in milling by GA approach", 《PROCEEDINGS OF COBEM》, 14 October 2003 (2003-10-14) * |
李忠群: "复杂切削条件高速铣削加工动力学建模、仿真与切削参数优化研究", 《北京航空航天大学博士学位论文》, 15 December 2008 (2008-12-15) * |
李琦 等: "基于遗传算法的切削工艺参数优化", 《机械制造》, vol. 42, no. 474, 29 February 2004 (2004-02-29) * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104834271A (zh) * | 2015-06-03 | 2015-08-12 | 上海理工大学 | 基于动刚度评价的直接进给轴伺服参数优化方法 |
CN104834271B (zh) * | 2015-06-03 | 2017-08-25 | 上海理工大学 | 基于动刚度评价的直接进给轴伺服参数优化方法 |
CN105699035A (zh) * | 2016-01-28 | 2016-06-22 | 哈尔滨理工大学 | 一种车削大螺距梯形外螺纹的振动测试方法 |
CN105699035B (zh) * | 2016-01-28 | 2018-01-12 | 哈尔滨理工大学 | 一种车削大螺距梯形外螺纹的振动测试方法 |
CN106002486A (zh) * | 2016-07-27 | 2016-10-12 | 大连理工大学 | 基于微铣床主传动系统功率的微铣削力测量方法 |
CN109406189A (zh) * | 2018-12-17 | 2019-03-01 | 武汉科技大学 | 一种五轴机床的刚度场半解析构建方法 |
CN110044625A (zh) * | 2019-04-30 | 2019-07-23 | 东华大学 | 一种飞轮轴承动刚度的测量方法 |
CN110103078A (zh) * | 2019-06-11 | 2019-08-09 | 贵州理工学院 | 一种基于普通车床的工艺系统的动刚度测定装置 |
CN114137907A (zh) * | 2021-11-30 | 2022-03-04 | 清华大学 | 数控机床振动幅值预测的建模方法、装置、设备及介质 |
CN114988679A (zh) * | 2022-08-08 | 2022-09-02 | 启东市云鹏玻璃机械有限公司 | 一种玻璃切割机床智能走刀控制方法 |
CN114988679B (zh) * | 2022-08-08 | 2022-10-25 | 启东市云鹏玻璃机械有限公司 | 一种玻璃切割机床智能走刀控制方法 |
CN117991718A (zh) * | 2022-10-18 | 2024-05-07 | 力劲精密机械(昆山)有限公司 | 一种卧式加工中心多轴协同状态监测方法和系统 |
CN117991718B (zh) * | 2022-10-18 | 2024-09-10 | 力劲精密机械(昆山)有限公司 | 一种卧式加工中心多轴协同状态监测方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103268430B (zh) | 2016-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103268430A (zh) | 基于机床刀具动刚度测量的铣削工艺参数优化方法 | |
CN102929210B (zh) | 基于特征的数控加工过程控制和优化系统及方法 | |
CN102248209B (zh) | 薄壁复杂曲面工件铣削加工时机床极限稳定工艺参数的确定方法 | |
CN107457609B (zh) | 基于刚度变化的铣削颤振抑制方法及铣削颤振优化系统 | |
CN203894596U (zh) | 数控机床加工状态多参数在线主动监控系统 | |
CN104517033B (zh) | 一种面向能量效率的数控加工工艺参数多目标优化方法 | |
CN103252680B (zh) | 一种铣削加工颤振主动控制系统及其方法 | |
CN107368665B (zh) | 高进给车削外螺纹件时变动力学模型构建方法 | |
CN101905340B (zh) | 一种高速铣削稳定性快速判定方法 | |
CN103433807B (zh) | 一种铣削力模型工艺参数的优化方法 | |
CN101930223A (zh) | 基于难加工金属数控加工工艺的智能筛选系统 | |
CN108405941A (zh) | 航空发动机叶片叶身型面高效精密铣削加工方法 | |
CN105750570B (zh) | 一种电主轴铣削颤振时延主动控制方法及其系统 | |
CN107145662B (zh) | 一种介观尺度车削变形预测方法 | |
Bagci | Monitoring and analysis of MRR-based feedrate optimization approach and effects of cutting conditions using acoustic sound pressure level in free-form surface milling | |
CN102081376A (zh) | 一种基于指令序列优化的加工负荷控制系统 | |
CN104898568A (zh) | 基于刚度辨识的数控机床进给系统控制参数优化方法 | |
CN103273387A (zh) | 基于刚度测量的外圆磨削工艺参数优化方法 | |
Zhuang et al. | On cutting parameters selection for plunge milling of heat-resistant-super-alloys based on precise cutting geometry | |
CN101870073A (zh) | 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法 | |
CN108647413B (zh) | 一种微细表面位置误差与稳定性综合预测方法 | |
CN105242634A (zh) | 一种数控机床加工极限切削深度预测方法 | |
CN104656562B (zh) | 一种柔性结构件铣削加工的工艺参数优化方法 | |
CN103926872A (zh) | 一种基于切削激励的数控机床激励方法 | |
Sampath et al. | Modeling and prediction of cutting noise in the face-milling process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160810 Termination date: 20190620 |
|
CF01 | Termination of patent right due to non-payment of annual fee |