CN113177280A - 球刀切削力系数标定方法 - Google Patents
球刀切削力系数标定方法 Download PDFInfo
- Publication number
- CN113177280A CN113177280A CN202110552105.8A CN202110552105A CN113177280A CN 113177280 A CN113177280 A CN 113177280A CN 202110552105 A CN202110552105 A CN 202110552105A CN 113177280 A CN113177280 A CN 113177280A
- Authority
- CN
- China
- Prior art keywords
- ball
- force coefficient
- end mill
- cutting force
- coefficient
- 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
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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Milling Processes (AREA)
Abstract
本发明公开了一种球刀切削力系数标定方法,包括如下步骤:1)构建球头铣刀的剪切力系数方程;2)对参考球头铣刀在不同加工参数下进行开槽试验,得到切削力数据;3)利用切削力数据得到参考球头铣刀旋转一周的平均切削力;4)利用平均切削力方程求解犁耕力系数和剪切力系数方程的常数,得到剪切力系数方程;5)根据待标定球头铣刀的几何参数和加工参数得到未知数il,利用所述剪切力系数方程得到该待标定球头铣刀的剪切力系数。本实施例的球刀切削力系数标定方法,采用平均切削力变系数法,仅用一个参考球头铣刀,即可辨识出犁耕力系数和剪切力系数方程;在球头铣刀及工件材料不变的情况下,利用该方法可以预测新的球头铣刀的切削力系数。
Description
技术领域
本发明涉及一种切削力系数计算方法,具体的为一种球刀切削力系数标定方法。
背景技术
切削力是影响零件表面加工性能和精度的重要因素之一。因此,高效精确的切削力预测模型一直是加工系统动态响应、颤振预测和工艺规划等研究的重点。可靠的铣削力预测不仅仅取决于准确的铣削力模型,还依赖于准确的切削力系数。目前,切削力系数的识别主要有两种方法:机械识别法和正交到斜角转换法。
机械识别法是通过使用特定的刀具进行切削试验来识别切削力系数。切削试验是在一定切削速度和进给速度的范围内进行。每次试验都需要测量切削力,计算出一个周期的平均力,并根据其平均力计算出切削力系数。Gradisek等人提出了一种通过平均槽切铣削力识别平头立铣刀切削常系数的方法。Wang等人利用傅里叶级数零频项,分别提出了用于球头铣刀和立铣刀的切削力系数识别方法。然而,上述方法没有考虑有效切削半径对球刀的剪切力系数的影响。另外,现有的机械识别法都只能求解出单一刀具的切削力系数,当刀具的几何形状发生改变时,需要重复所有的切削试验,这意味着需要消耗大量的时间和金钱成本。
正交到斜角转换法相比于机械识别法更加的灵活。当刀具的几何形状(螺旋角)发生改变时,该方法可以无需额外的切削实验来计算新刀具的切削力系数。通过一组不同切削速度和进给量的正交切削试验,计算出当前刀具材料下的正交数据库,即剪切应力、剪切角、摩擦角和犁耕力系数。然后,基于正交到斜角转换法,快速计算出具有不同螺旋角刀具的切削力系数。Lee和Altinas提出了通过正交车削试验来识别球头铣刀切削系数的方法。该方法中剪切力系数通过经典的斜角转换公式来计算,犁耕力系数通过将剪切力系数外推至切削厚度为零时获得。
机械识别法和正交到斜角转换法各有优缺点,尤其是面对多把刀具的切削力系数识别时,存在试验成本较高、物理参数测量不方便等问题,无法取得令人满意的结果。
发明内容
球头铣刀的切削力系数是变化的,现有技术中建立了基于轴向高度或轴向位置角的预测模型;然而,切削力系数变化的根本原因并不是由轴向高度或轴向位置角这两个因素导致的,其根本原因在于球头铣刀的切削刃倾角会发生变化,进而导致切削力系数发生改变。
有鉴于此,本发明的目的在于提供一种球刀切削力系数标定方法,能够满足球头铣刀的切削力系数的标定要求。
为达到上述目的,本发明提供如下技术方案:
一种球刀切削力系数标定方法,包括如下步骤:
1)构建球头铣刀的剪切力系数方程;
2)对参考球头铣刀在不同加工参数下进行开槽试验,得到切削力数据;
3)利用切削力数据得到参考球头铣刀旋转一周的平均切削力;
4)利用平均切削力方程求解犁耕力系数和剪切力系数方程的常数,得到以il为未知数的剪切力系数方程,且:
il=tan(ih)/tan(i0)
其中,i0为球头铣刀的螺旋角;ih为球头铣刀的局部螺旋角;
5)根据待标定球头铣刀的几何参数和加工参数得到未知数il,利用所述剪切力系数方程得到该待标定球头铣刀的剪切力系数。
进一步,所述步1)中,球头铣刀的剪切力系数方程的构建方法如下:
将球头铣刀沿其轴线分为Ne层,每层的局部螺旋角ih为:
ih=arctan(R(i)·tan(i0)/R)
其中,R(i)为第i层对应的球头铣刀的半径;R为球头铣刀的半径,i0为球头铣刀的螺旋角;
将球头铣刀的剪切力系数假定为局部螺旋角ih的三次多项式函数,得到球头铣刀的剪切力系数方程为:
其中,Ktc,Krc,Kac表示球刀剪切力系数;C0-C11为常数。
进一步,所述步骤4)中,球头铣刀旋转一周的平均切削力表示为:
其中,表示X方向上的平均切削力;表示Y方向上的平均切削力;表示Z方向上的平均切削力;N表示球头铣刀的刃数;表示局部螺旋角的上限;ψ表示切削刃微元点的径向位置角;[A]表示局部螺旋角与切削力的关系矩阵;[K]表示切削力系数矩阵;且:
其中,fc表示每齿进给量;Q、G、J、P、L和H均为无实际含义的变形参数;且:
[K]=[Kte Kre Kae C0 C4 C8 C1 C5 C9 C2 C6 C10 C3 C7 C11]T
其中,Kte,Kre,Kae表示球刀犁耕力系数。
进一步,由于[K]中包含15个未知数,因此所述步骤2)中,对参考球头铣刀在不同加工参数下进行至少5次开槽试验,在所述步骤4)中,利用最小二乘法辨识求解得到的犁耕力系数和剪切力系数方程的常数。
本发明的有益效果在于:
本发明的球刀切削力系数标定方法,采用平均切削力变系数法,仅用一个参考球头铣刀,即可辨识出犁耕力系数和得到以il为未知数的剪切力系数方程;由于未知数il与球头铣刀的局部螺旋角相关,也即本发明构建得到的剪切力系数方程表征了切削力系数随切削刃倾角变化而改变的原理;在球头铣刀及工件材料不变的情况下,利用该方法可以预测新的球头铣刀的切削力系数,从而显著减少了试验时间和成本,这对球头铣刀的设计和制造具有重要意义。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为球头铣刀切削力模型示意图;
图2为图1的仰视图。
其中,图1中,κ表示切削刃微元点的轴向侵入角,dz表示每一层微元的高度;图2中,ω表示旋转方向,θ表示切削刃位置角,φz表示螺旋滞后角,ψ表示切削刃微元点的径向位置角。
具体实施方式
下面结合附图和具体实施例对本发明作进一步说明,以使本领域的技术人员可以更好的理解本发明并能予以实施,但所举实施例不作为对本发明的限定。
在铣削加工中,由于球头铣刀的几何特性,局部螺旋角会导致切削刃倾角的变化。切削力系数由犁耕力系数(Kte,Kre,Kae)和剪切力系数(Ktc,Krc,Kac)组成。需要指出的是,局部螺旋角对犁耕力系数的影响可以忽略不计,因此,犁耕力系数可以视为常数,只需定量分析局部螺旋角与剪切力系数之间的关系。
具体的,本实施例的球刀切削力系数标定方法,包括如下步骤:
1)构建球头铣刀的剪切力系数方程。
具体的,本实施例的球头铣刀的剪切力系数方程的构建方法如下:
将球头铣刀沿其轴线分为Ne层,每层的局部螺旋角ih为:
ih=arctan(R(i)·tan(i0)/R)
其中,R(i)为第i层对应的球头铣刀的半径;R为球头铣刀的半径,i0为球头铣刀的螺旋角;
将球头铣刀的剪切力系数假定为局部螺旋角ih的三次多项式函数,得到球头铣刀的剪切力系数方程为:
其中,Ktc,Krc,Kac表示球刀剪切力系数;C0-C11为常数。
2)对参考球头铣刀在不同加工参数下进行开槽试验,得到切削力数据。具体的,加工参数包括:切削深度ap、主轴转速n、每齿进给量fp。
3)利用切削力数据得到参考球头铣刀旋转一周的平均切削力;
4)利用平均切削力方程求解犁耕力系数和剪切力系数方程的常数,得到以il为未知数的剪切力系数方程,且:
il=tan(ih)/tan(i0)
其中,i0为球头铣刀的螺旋角;ih为球头铣刀的局部螺旋角。
具体的,球头铣刀旋转一周的平均切削力表示为:
其中,表示X方向上的平均切削力;表示Y方向上的平均切削力;表示Z方向上的平均切削力;N表示球头铣刀的刃数;表示局部螺旋角的上限;ψ表示切削刃微元点的径向位置角;[A]表示局部螺旋角与切削力的关系矩阵;[K]表示切削力系数矩阵;且:
其中,fc表示每齿进给量;Q、G、J、P、L和H均为无实际含义的变形参数;且:
[K]=[Kte Kre Kae C0 C4 C8 C1 C5 C9 C2 C6 C10 C3 C7 C11]T
其中,Kte,Kre,Kae表示球刀犁耕力系数。
5)根据待标定球头铣刀的几何参数和加工参数得到未知数il,利用所述剪切力系数方程得到该待标定球头铣刀的剪切力系数。
优选的,由于[K]中包含15个未知数,对参考球头铣刀在不同加工参数下进行1次开槽试验,获得三个方程,因此所述步骤2)中,对参考球头铣刀在不同加工参数下进行至少5次开槽试验,即可从理论上计算出所有的未知量。然而,由于开槽过程中不可避免的误差,实验值和理论值可能会有所不同,因此,为了减小实验误差对该模型的影响,需要对参考球头铣刀在不同加工参数下进行至少5次开槽试验。本实施例对参考球头铣刀在不同加工参数下进行至少6次开槽试验,在所述步骤4)中,利用最小二乘法辨识求解得到的犁耕力系数和剪切力系数方程的常数。
本实施例的球刀切削力系数标定方法,采用平均切削力变系数法,仅用一个参考球头铣刀,即可辨识出犁耕力系数和得到以il为未知数的剪切力系数方程;由于未知数il与球头铣刀的局部螺旋角相关,也即本实施例构建得到的剪切力系数方程表征了切削力系数随切削刃倾角变化而改变的原理;在球头铣刀及工件材料不变的情况下,利用该方法可以预测新的球头铣刀的切削力系数,从而显著减少了试验时间和成本,这对球头铣刀的设计和制造具有重要意义。
以上所述实施例仅是为充分说明本发明而所举的较佳的实施例,本发明的保护范围不限于此。本技术领域的技术人员在本发明基础上所作的等同替代或变换,均在本发明的保护范围之内。本发明的保护范围以权利要求书为准。
Claims (4)
1.一种球刀切削力系数标定方法,其特征在于:包括如下步骤:
1)构建球头铣刀的剪切力系数方程;
2)对参考球头铣刀在不同加工参数下进行开槽试验,得到切削力数据;
3)利用切削力数据得到参考球头铣刀旋转一周的平均切削力;
4)利用平均切削力方程求解犁耕力系数和剪切力系数方程的常数,得到以il为未知数的剪切力系数方程,且:
il=tan(ih)/tan(i0)
其中,i0为球头铣刀的螺旋角;ih为球头铣刀的局部螺旋角;
5)根据待标定球头铣刀的几何参数和加工参数得到未知数il,利用所述剪切力系数方程得到该待标定球头铣刀的剪切力系数。
4.根据权利要求3所述的球刀切削力系数标定方法,其特征在于:由于[K]中包含15个未知数,因此所述步骤2)中,对参考球头铣刀在不同加工参数下进行至少5次开槽试验,在所述步骤4)中,利用最小二乘法辨识求解得到的犁耕力系数和剪切力系数方程的常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110552105.8A CN113177280B (zh) | 2021-05-20 | 2021-05-20 | 球刀切削力系数标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110552105.8A CN113177280B (zh) | 2021-05-20 | 2021-05-20 | 球刀切削力系数标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113177280A true CN113177280A (zh) | 2021-07-27 |
CN113177280B CN113177280B (zh) | 2022-10-25 |
Family
ID=76929402
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110552105.8A Active CN113177280B (zh) | 2021-05-20 | 2021-05-20 | 球刀切削力系数标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113177280B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113553727A (zh) * | 2021-08-12 | 2021-10-26 | 东北大学 | 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法 |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20010087976A (ko) * | 2000-03-09 | 2001-09-26 | 정융호 | 고속가공용 절삭부하 예측방법 |
US20040258495A1 (en) * | 2003-05-08 | 2004-12-23 | Yoshiaki Kakino | Cutting force detection method and machining control method and apparatus based on detected cutting force |
WO2012153157A2 (en) * | 2011-05-09 | 2012-11-15 | Pécsi Tudományegyetem | Method for optimizing cutting forces in a milling process and computer-controlled milling machine using the same method |
JP2013054431A (ja) * | 2011-09-01 | 2013-03-21 | Hiroshima Prefecture | 回転工具の実切削距離算出方法、実切削距離算出のためのコンピュータプログラム、切削力予測方法及び工具経路修正装置 |
US20140216170A1 (en) * | 2013-02-05 | 2014-08-07 | Georgia Tech Research Corporation | Systems And Methods For Monitoring Cutting Forces In Peripheral End Milling |
CN105069257A (zh) * | 2015-09-01 | 2015-11-18 | 华中科技大学 | 一种自由曲面微细铣削切削力建模方法 |
CN106808320A (zh) * | 2016-12-12 | 2017-06-09 | 华中科技大学 | 一种考虑刀具齿间角、螺旋角及偏心的铣削力预测方法 |
CN107944149A (zh) * | 2017-11-28 | 2018-04-20 | 清华大学 | 一种普通立铣刀圆弧段切削力系数辨识方法 |
CN107944176A (zh) * | 2017-12-06 | 2018-04-20 | 上海复合材料科技有限公司 | 一种球头铣刀多轴铣削钛合金铣削力预测方法 |
CN109656192A (zh) * | 2018-12-18 | 2019-04-19 | 华中科技大学 | 一种基于球头铣刀铣削力的加工优化方法 |
CN109732404A (zh) * | 2018-12-28 | 2019-05-10 | 哈尔滨理工大学 | 一种基于后刀面磨损的球头铣刀铣削力确定方法及系统 |
CN110270705A (zh) * | 2019-06-26 | 2019-09-24 | 东北大学 | 考虑刀具跳动影响的微细铣削加工切削力仿真预测方法 |
CN110405533A (zh) * | 2019-07-23 | 2019-11-05 | 重庆大学 | 一种刀具姿态优化方法 |
CN111459096A (zh) * | 2020-04-07 | 2020-07-28 | 西南交通大学 | 一种非均匀磨损圆弧头立铣刀的侧铣铣削力计算方法 |
-
2021
- 2021-05-20 CN CN202110552105.8A patent/CN113177280B/zh active Active
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20010087976A (ko) * | 2000-03-09 | 2001-09-26 | 정융호 | 고속가공용 절삭부하 예측방법 |
US20040258495A1 (en) * | 2003-05-08 | 2004-12-23 | Yoshiaki Kakino | Cutting force detection method and machining control method and apparatus based on detected cutting force |
WO2012153157A2 (en) * | 2011-05-09 | 2012-11-15 | Pécsi Tudományegyetem | Method for optimizing cutting forces in a milling process and computer-controlled milling machine using the same method |
JP2013054431A (ja) * | 2011-09-01 | 2013-03-21 | Hiroshima Prefecture | 回転工具の実切削距離算出方法、実切削距離算出のためのコンピュータプログラム、切削力予測方法及び工具経路修正装置 |
US20140216170A1 (en) * | 2013-02-05 | 2014-08-07 | Georgia Tech Research Corporation | Systems And Methods For Monitoring Cutting Forces In Peripheral End Milling |
CN105069257A (zh) * | 2015-09-01 | 2015-11-18 | 华中科技大学 | 一种自由曲面微细铣削切削力建模方法 |
CN106808320A (zh) * | 2016-12-12 | 2017-06-09 | 华中科技大学 | 一种考虑刀具齿间角、螺旋角及偏心的铣削力预测方法 |
CN107944149A (zh) * | 2017-11-28 | 2018-04-20 | 清华大学 | 一种普通立铣刀圆弧段切削力系数辨识方法 |
CN107944176A (zh) * | 2017-12-06 | 2018-04-20 | 上海复合材料科技有限公司 | 一种球头铣刀多轴铣削钛合金铣削力预测方法 |
CN109656192A (zh) * | 2018-12-18 | 2019-04-19 | 华中科技大学 | 一种基于球头铣刀铣削力的加工优化方法 |
CN109732404A (zh) * | 2018-12-28 | 2019-05-10 | 哈尔滨理工大学 | 一种基于后刀面磨损的球头铣刀铣削力确定方法及系统 |
CN110270705A (zh) * | 2019-06-26 | 2019-09-24 | 东北大学 | 考虑刀具跳动影响的微细铣削加工切削力仿真预测方法 |
CN110405533A (zh) * | 2019-07-23 | 2019-11-05 | 重庆大学 | 一种刀具姿态优化方法 |
CN111459096A (zh) * | 2020-04-07 | 2020-07-28 | 西南交通大学 | 一种非均匀磨损圆弧头立铣刀的侧铣铣削力计算方法 |
Non-Patent Citations (5)
Title |
---|
JEONG HOON KO等: "3D Ball-End Milling Force Model Using Instantaneous Cutting Force Coefficients", 《J. MANUF. SCI. ENG》 * |
ZENGYA ZHAO等: "Ball-end milling cutter design method towards the maximum material removal rate under surface roughness constraints", 《JOURNAL OF MANUFACTURING PROCESSES》 * |
丁云鹏: "多轴联动数控加工球头铣刀铣削力建模与仿真", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅰ辑》 * |
王海艳等: "基于灰色系统理论的球头铣刀螺旋铣孔切削参数优化", 《中国工程机械学报》 * |
郭一君等: "基于改进人工蜂群算法和极限学习机的刀具磨损监测", 《重庆大学学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113553727A (zh) * | 2021-08-12 | 2021-10-26 | 东北大学 | 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113177280B (zh) | 2022-10-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106808320B (zh) | 一种考虑刀具齿间角、螺旋角及偏心的铣削力预测方法 | |
Wojciechowski | The estimation of cutting forces and specific force coefficients during finishing ball end milling of inclined surfaces | |
CN105069257B (zh) | 一种自由曲面微细铣削切削力建模方法 | |
CN103646141B (zh) | 一种平底螺旋立铣刀正交车铣轴类零件的切削力建模方法 | |
CN101491844B (zh) | 圆周铣削过程铣削力系数与刀具径向偏心标定方法 | |
CN106682281A (zh) | 基于最大切削力的铣削瞬时切削力预测方法 | |
CN113177280B (zh) | 球刀切削力系数标定方法 | |
CN103258095A (zh) | 平底立铣刀通用铣削力建模方法 | |
CN105511397B (zh) | 统一犁切模型的通用铣削力建模方法 | |
CN111360582A (zh) | 一种刀具磨损状态辨识方法 | |
CN113158371A (zh) | 一种高速铣削加工动态切削力预测系统及参数优化方法 | |
Akyildiz | Evaluating of cutting forces in thread machining | |
CN107168245A (zh) | 一种考虑刃口效应的倒棱圆形车刀切削力精确预测方法 | |
CN107944149A (zh) | 一种普通立铣刀圆弧段切削力系数辨识方法 | |
Chen et al. | A unified analytical cutting force model for variable helix end mills | |
CN113400092B (zh) | 考虑材料堆积的金属切削力预测方法 | |
CN107025359B (zh) | 左、右切削刃的刀工界面切削行为差异性的计算与测试方法 | |
CN110253341B (zh) | 一种微铣削加工刀具跳动参数的快速识别方法 | |
CN113626953B (zh) | 高能效铣削加工误差动态分布特性识别方法 | |
Hu et al. | An effective thread milling force prediction model considering instantaneous cutting thickness based on the cylindrical thread milling simplified to side milling process | |
CN111553032B (zh) | 整体叶盘铣削温度预测方法、装置、设备及存储介质 | |
Zhang et al. | A cutting force model for a waved-edge end milling cutter | |
CN113553727A (zh) | 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法 | |
CN114036671A (zh) | 一种考虑主轴热误差和刀具跳动的切削力模型建模方法 | |
Yang et al. | Examination of surface milling performance of micro-textured ball milling cutter |
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 |