CN113792467A - 高速铣削adc12铝合金刀具磨损预测方法 - Google Patents
高速铣削adc12铝合金刀具磨损预测方法 Download PDFInfo
- Publication number
- CN113792467A CN113792467A CN202111096328.4A CN202111096328A CN113792467A CN 113792467 A CN113792467 A CN 113792467A CN 202111096328 A CN202111096328 A CN 202111096328A CN 113792467 A CN113792467 A CN 113792467A
- Authority
- CN
- China
- Prior art keywords
- cutter
- wear
- model
- aluminum alloy
- workpiece
- 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.)
- Pending
Links
Images
Classifications
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- 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/08—Thermal analysis or thermal optimisation
-
- 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)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明提出一种高速铣削ADC12铝合金刀具磨损预测方法,采用有限元仿真技术,考虑材料的弹塑性本构关系,输出刀具磨损值情况,再基于刀具磨损几何关系,实现一种ADC12铝合金高速铣削过程刀具磨损的精准预测。与传统刀具磨损试验研究相比显著降低经济成本,操作简单,精度满足要求。预测方法通过对工件和刀具建模,考虑材料弹塑性本构关系,建立刀具‑工件摩擦类型、金属切削分离准则,可以结合Python语言输出有限元仿真不同切削参数下刀具磨损值,再通过移动节点位置确定铣刀磨损后几何形状及前、后刀面磨损值,预测高速铣削过程刀具磨损情况。
Description
技术领域
本发明属于高速铣削加工及模拟仿真技术领域,涉及一种高速铣削ADC12铝合金刀具磨损预测方法,尤其是一种通过仿真建模推导ADC12铝合金高速铣削过程刀具磨损的预测方法。
背景技术
由于汽车行业正朝着轻量化、低成本、节能减排的方向发展,ADC12铝合金因其良好的综合性能,被广泛用于制作汽车发动机缸体、缸盖、变速箱以达到减重的需求。这些零件因其较高的技术要求需要采用高速铣削进行加工,但由于ADC12铝合金材料特性,采用高速铣削易产生刀具磨损。刀具磨损会严重影响工件尺寸精度、表面质量及加工成本,是个需要重视的问题,若能准确地预测刀具磨损,将会提高加工效率,节省加工成本。
目前,关于刀具磨损的研究大多数是依靠大量的试验和经验性的刀具寿命公式进行计算,这些方法周期长、成本高、效率低,无法实现量化分析。近几十年,随着科技的进步,有限元方法备受关注,若可以使用有限元方法模拟刀具的磨损,预测刀具磨损及寿命将变得事半功倍。因此,对ADC12铝合金高速铣削刀具磨损的仿真进行研究具有重要意义与实用价值。
发明内容
考虑到现有技术存在的空白和不足,本发明的目的在于提供一种高速铣削ADC12铝合金刀具磨损预测方法,该方法采用有限元仿真技术,考虑材料的弹塑性本构关系,输出刀具磨损值情况,再基于刀具磨损几何关系,实现一种ADC12铝合金高速铣削过程刀具磨损的精准预测。与传统刀具磨损试验研究相比显著降低经济成本,操作简单,精度满足要求。
本发明采用的技术方案是一种ADC12铝合金高速铣削过程刀具磨损的预测方法,采用有限元仿真技术,预测方法通过对工件和刀具建模,考虑材料弹塑性本构关系,建立刀具-工件摩擦类型、金属切削分离准则,可以结合Python语言输出有限元仿真不同切削参数下刀具磨损值,再通过移动节点位置确定铣刀磨损后几何形状及前、后刀面磨损值,预测高速铣削过程刀具磨损情况。
本发明具体采用以下技术方案:
一种高速铣削ADC12铝合金刀具磨损预测方法,其特征在于,包括以下步骤:
步骤S1:通过有限元模拟,进行切屑成形过程仿真;并提取刀-屑间的接触应力、相对滑移速度及温度;
步骤S2:通过切削试验获得磨损率数据;
步骤S3:根据步骤S1和步骤S2的结果,构建磨损模型;
步骤S4:进行仿真,并计算Δt时间后刀具各节点的磨损率;
步骤S5:根据刀具各节点的磨损方向计算出刀具各节点的位移;
步骤S6:判断磨损率是否达到预设值,如果达到,则输出磨损模型,如果未达到,则对刀具形貌进行更新及网格重划分并回到步骤S1;
步骤S7:根据步骤S6输出的磨损模型进行高速铣削ADC12铝合金刀具磨损预测。
进一步地,在步骤S1中,通过有限元模拟,进行切屑成形过程仿真具体包括以下步骤:
步骤S11:测绘铣刀几何结构,基于此建立铣刀模型,导入ABAQUS软件中,将其设置为刚体;
步骤S12:建立高速铣削工件模型,设置为弹塑性体;
步骤S13:对模型进行网格划分,网格单元采用4节点非线性单元,并将刀具及工件接触切削区域的网格局部加密;
步骤S14:采用Johnson-Cook本构模型描述真实切削过程材料内流动应力与应变关系;
步骤S15:材料损伤判据采用Johnson-Cook失效模型,当材料应变超过失效应变后,材料失效,失效单元被删除,切屑开始形成;
步骤S16:定义接触性质,刀尖部位为粘结接触,后刀面与工件间为滑动接触;
步骤S17:定义边界条件,设置工件约束;设置刀具与工件之间接触条件及工件与环境之间热传导;
步骤S18:设置主轴转速、进给速度和切削深度;将计算设置为热力耦合模式,运行有限元模拟软件,进行连续性切屑成形仿真。
进一步地,在步骤S14中,ADC12铝合金的J-C本构模型表达式为:
式中,σ为Von-Mises流动应力;A为参考温度和参考应变速率下的屈服强度;B为应变强化系数;εp为等效塑性应变;n为应变硬化指数;C为应变速率硬化系数;为参考应变速率;Tmelt为材料熔点温度;T0为室温;m为加热软化指数。
进一步地,在步骤S15中,Johnson-Cook失效模型的表达式如式(2)所示:
进一步地,在步骤S16中,定义的接触性质为:
其中,μ为滑移区接触的摩擦系数,是固定值;τmax为极限剪切应力,值等于σy/3,σy为接触面的工件材料的等效屈服应力。
进一步地,在步骤S3中,所述磨损模型采用Usui模型,表达式为:
进一步地,在步骤S6中,当达到仿真所设定的时间时,停止循环,输出磨损模型。
与现有技术相比,本发明及其优选方案基于有限元仿真技术,针对ADC12铝合金高速铣削过程中刀具易产生粘结磨损的特点,提出了新的计算模型,使其可以利用ABAQUS软件及Python语言利用计算机的算力辅助进行刀具磨损预测,为选择合理切削参数提供参考,与磨损实验研究相比降低了经济成本,且操作简单。
附图说明
下面结合附图和具体实施方式对本发明进一步详细的说明:
图1为本发明实施例ADC12铝合金高速铣削刀具磨损预测方法流程示意图。
具体实施方式
为让本专利的特征和优点能更明显易懂,下文特举实施例,作详细说明如下:
如图1所示,本实施例以盘铣刀高速铣削平面为例,刀片采用金鹭公司生产的APKT1604PDFR-MA H01系列铣刀片,刀具主切削刃前角为0°,后角为11°,借助ABAQUS软件进行仿真。
第一步:首先确定选用工件材料ADC12铝合金的材料参数性能,密度:2670kg·m-3;硬度:74.1-80HB;弹性模量:76GPa;泊松比:0.33;热导率:92.6W·mol-1·K-1;屈服应力:170MPa;抗拉强度:230MPa。
第二步:将选择的铣刀片测绘,在ABAQUS软件中构建模型,导入刀具材料参数,并将其设置为刚体,在与工件接触的切削区域使用网格局部加密,采用4节点平面应变单元。模型工件为半圆弧形状,设置为弹塑性材料,网格采用四边形结构网格类型,在与刀具接触的切削区域使用网格局部加密。在有限元模型中创建工件材料ADC12,并根据上述材料性能参数定义,将材料属性赋予工件。
第三步:金属塑性行为采用Johnson-Cook本构进行模拟,在公式(1)中参数依次取为:参考温度和参考应变速率下的屈服强度A为240MPa;应变强化系数B为636MPa;应变速率硬化系数C为0.006998;加热软化指数m为1.8;应变硬化指数n为0.95;室温T0为20°。
第四步:材料损伤判据采用Johnson-Cook失效模型,式(2)中参数d1、d2、d3、d4、d5依次是,0.0939,0.00001,-29.4,1.5,0。
第五步:设置接触性质,依据公式(3),设置剪切摩擦系数λ为0.2;设置库伦摩擦系数μ为0.7。定义边界条件,将工件侧面及底面自由度严格约束,设刀具工件间热传导系数为6000kW/(m2·K-1);设工件与环境间热传导系数为18kW/(m2·K-1)。
第六步:设置刀具切削速度、进给量及切削深度,进行热力耦合计算,运行有限元模拟过程,得到仿真结果。
第七步:使用编写好的子程序从结果数据库中提取刀-屑间的接触应力、相对滑移速度及温度等相关参数。
第八步:针对ADC12铝合金高速铣削过程中出现的粘结磨损,磨损模型采用Usui模型,由切屑成形仿真得到正应力σn、滑动速度vs和温度T,以及试验得到的磨损率使用曲线拟合的方法得到Usui模型中的常数C、λ。
第九步:计算Δt时间后刀具各节点的磨损率,再根据刀具各节点的磨损方向计算出刀具各节点的位移,最后对刀具形貌进行更新及网格重划分,实现铝合金ADC12高速铣削过程刀具磨损预测,可循环以上步骤,直至刀具达到磨损标准或达到仿真所设定的时间。
本实施例基于有限元仿真技术,针对ADC12铝合金高速铣削过程中刀具易产生粘结磨损的特点,使用有限元仿真软件ABAQUS软件的二次开发技术进行刀具磨损预测,为选择合理切削参数提供参考,与磨损实验研究相比降低了经济成本,且操作简单。
本专利不局限于上述最佳实施方式,任何人在本专利的启示下都可以得出其它各种形式的高速铣削ADC12铝合金刀具磨损预测方法,凡依本发明申请专利范围所做的均等变化与修饰,皆应属本专利的涵盖范围。
Claims (7)
1.一种高速铣削ADC12铝合金刀具磨损预测方法,其特征在于,包括以下步骤:
步骤S1:通过有限元模拟,进行切屑成形过程仿真;并提取刀-屑间的接触应力、相对滑移速度及温度;
步骤S2:通过切削试验获得磨损率数据;
步骤S3:根据步骤S1和步骤S2的结果,构建磨损模型;
步骤S4:进行仿真,并计算Δt时间后刀具各节点的磨损率;
步骤S5:根据刀具各节点的磨损方向计算出刀具各节点的位移;
步骤S6:判断磨损率是否达到预设值,如果达到,则输出磨损模型,如果未达到,则对刀具形貌进行更新及网格重划分并回到步骤S1;
步骤S7:根据步骤S6输出的磨损模型进行高速铣削ADC12铝合金刀具磨损预测。
2.根据权利要求1所述的高速铣削ADC12铝合金刀具磨损预测方法,其特征在于:
在步骤S1中,通过有限元模拟,进行切屑成形过程仿真具体包括以下步骤:
步骤S11:测绘铣刀几何结构,基于此建立铣刀模型,导入ABAQUS软件中,将其设置为刚体;
步骤S12:建立高速铣削工件模型,设置为弹塑性体;
步骤S13:对模型进行网格划分,网格单元采用4节点非线性单元,并将刀具及工件接触切削区域的网格局部加密;
步骤S14:采用Johnson-Cook本构模型描述真实切削过程材料内流动应力与应变关系;
步骤S15:材料损伤判据采用Johnson-Cook失效模型,当材料应变超过失效应变后,材料失效,失效单元被删除,切屑开始形成;
步骤S16:定义接触性质,刀尖部位为粘结接触,后刀面与工件间为滑动接触;
步骤S17:定义边界条件,设置工件约束;设置刀具与工件之间接触条件及工件与环境之间热传导;
步骤S18:设置主轴转速、进给速度和切削深度;将计算设置为热力耦合模式,运行有限元模拟软件,进行连续性切屑成形仿真。
7.根据权利要求1所述的高速铣削ADC12铝合金刀具磨损预测方法,其特征在于:
在步骤S6中,当达到仿真所设定的时间时,停止循环,输出磨损模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111096328.4A CN113792467A (zh) | 2021-09-17 | 2021-09-17 | 高速铣削adc12铝合金刀具磨损预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111096328.4A CN113792467A (zh) | 2021-09-17 | 2021-09-17 | 高速铣削adc12铝合金刀具磨损预测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113792467A true CN113792467A (zh) | 2021-12-14 |
Family
ID=78878939
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111096328.4A Pending CN113792467A (zh) | 2021-09-17 | 2021-09-17 | 高速铣削adc12铝合金刀具磨损预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113792467A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114769692A (zh) * | 2022-04-19 | 2022-07-22 | 哈尔滨理工大学 | 高效铣刀后刀面摩擦磨损边界跨尺度识别方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102103646A (zh) * | 2010-12-14 | 2011-06-22 | 武汉理工大学 | 基于有限元技术和人工神经网络的精冲模具磨损预测方法 |
CN105643024A (zh) * | 2016-03-17 | 2016-06-08 | 哈尔滨理工大学 | 一种车削大螺距螺纹轴向分层切削方法、刀具磨损测试方法及其力热载荷计算方法 |
CN106407526A (zh) * | 2016-09-05 | 2017-02-15 | 大连理工大学 | 一种微铣削过程刀具后刀面磨损预测方法 |
CN110489931A (zh) * | 2019-09-18 | 2019-11-22 | 大连理工大学 | 一种基于切削过程仿真的薄壁微铣削力预测方法 |
CN110569586A (zh) * | 2019-08-28 | 2019-12-13 | 北京自动化控制设备研究所 | 行星滚柱丝杠副精密螺纹切削建模方法 |
-
2021
- 2021-09-17 CN CN202111096328.4A patent/CN113792467A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102103646A (zh) * | 2010-12-14 | 2011-06-22 | 武汉理工大学 | 基于有限元技术和人工神经网络的精冲模具磨损预测方法 |
CN105643024A (zh) * | 2016-03-17 | 2016-06-08 | 哈尔滨理工大学 | 一种车削大螺距螺纹轴向分层切削方法、刀具磨损测试方法及其力热载荷计算方法 |
CN106407526A (zh) * | 2016-09-05 | 2017-02-15 | 大连理工大学 | 一种微铣削过程刀具后刀面磨损预测方法 |
CN110569586A (zh) * | 2019-08-28 | 2019-12-13 | 北京自动化控制设备研究所 | 行星滚柱丝杠副精密螺纹切削建模方法 |
CN110489931A (zh) * | 2019-09-18 | 2019-11-22 | 大连理工大学 | 一种基于切削过程仿真的薄壁微铣削力预测方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114769692A (zh) * | 2022-04-19 | 2022-07-22 | 哈尔滨理工大学 | 高效铣刀后刀面摩擦磨损边界跨尺度识别方法 |
CN114769692B (zh) * | 2022-04-19 | 2023-01-03 | 哈尔滨理工大学 | 高效铣刀后刀面摩擦磨损边界跨尺度识别方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Özel et al. | Finite element modeling the influence of edge roundness on the stress and temperature fields induced by high-speed machining | |
Zanger et al. | Investigations on mechanisms of tool wear in machining of Ti-6Al-4V using FEM simulation | |
Özel et al. | Finite element method simulation of machining of AISI 1045 steel with a round edge cutting tool | |
CN105243195B (zh) | 一种微铣削镍基高温合金加工硬化的预测方法 | |
CN110597180B (zh) | 一种基于切削过程仿真的单晶铜微铣削力预测方法 | |
Denkena et al. | Inverse determination of constitutive equations and cutting force modelling for complex tools using Oxley's predictive machining theory | |
Wang et al. | A novel finite element method for the wear analysis of cemented carbide tool during high speed cutting Ti6Al4V process | |
Linjiang et al. | Inverse identification of constitutive parameters of Ti2AlNb intermetallic alloys based on cooperative particle swarm optimization | |
Ortiz-de-Zarate et al. | Experimental and FEM analysis of surface integrity when broaching Ti64 | |
Vishwakarma et al. | FEA modeling of material removal rate in electrical discharge machining of Al6063/SiC composites | |
Ren et al. | Inverse determination of improved constitutive equation for cutting titanium alloy Ti-6Al-4V based on finite element analysis | |
CN113792467A (zh) | 高速铣削adc12铝合金刀具磨损预测方法 | |
Yan et al. | Numerical simulation of finish hard turning for AISI H13 die steel | |
Tang et al. | Finite-element simulation of conventional and high-speed peripheral milling of hardened mold steel | |
Yue et al. | 3D FEM simulation of milling force in corner machining process | |
Long et al. | Finite element modeling of burr formation in orthogonal cutting | |
Cerutti | Numerical modelling and mechanical analysis of the machining of large aeronautical parts: Machining quality improvement | |
CN115544687A (zh) | 一种精密复杂刀具的切削性能预测方法 | |
Yang et al. | Study on mechanism of chip formation during high-speed milling of alloy cast iron | |
Huan et al. | Predicting cutting force and simulating material removal process in turning TiCp/Ti6Al4V composite | |
Akram et al. | A numerical investigation of effects of cutting velocity and feed rate on residual stresses in aluminum alloy Al-6061 | |
Zhang et al. | The effects of tool structure parameters on forming extrusion cutting (FEC) | |
Song et al. | Study on surface morphology and residual stress in inclined milling of titanium alloy TC11 | |
Necpal et al. | Assessment of cutting tool wear using a numerical FEM simulation model | |
Uthayakumar et al. | Precision machining of an aluminum alloy piston reinforced with a cast iron insert |
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 |