CN110032748A - 适用于纵扭超声铣削的有限元仿真建模方法 - Google Patents
适用于纵扭超声铣削的有限元仿真建模方法 Download PDFInfo
- Publication number
- CN110032748A CN110032748A CN201811530530.1A CN201811530530A CN110032748A CN 110032748 A CN110032748 A CN 110032748A CN 201811530530 A CN201811530530 A CN 201811530530A CN 110032748 A CN110032748 A CN 110032748A
- Authority
- CN
- China
- Prior art keywords
- cutter
- workpiece
- parameter
- cutting
- milling
- 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
- 238000003801 milling Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005094 computer simulation Methods 0.000 title claims abstract description 8
- 238000005520 cutting process Methods 0.000 claims abstract description 43
- 239000000463 material Substances 0.000 claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 238000012545 processing Methods 0.000 claims abstract description 8
- 239000000758 substrate Substances 0.000 claims abstract description 8
- 238000002604 ultrasonography Methods 0.000 claims description 11
- 238000013499 data model Methods 0.000 claims description 5
- 239000002245 particle Substances 0.000 claims description 4
- 238000005516 engineering process Methods 0.000 claims description 3
- 241000208340 Araliaceae Species 0.000 claims description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 claims description 2
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 235000008434 ginseng Nutrition 0.000 claims description 2
- 230000035945 sensitivity Effects 0.000 claims description 2
- 238000005482 strain hardening Methods 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 3
- 239000011159 matrix material Substances 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000001133 acceleration Effects 0.000 description 2
- 230000002929 anti-fatigue Effects 0.000 description 2
- 238000003754 machining Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000000737 periodic effect Effects 0.000 description 2
- 239000000956 alloy Substances 0.000 description 1
- 229910045601 alloy Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000007514 turning Methods 0.000 description 1
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Numerical Control (AREA)
- Milling Processes (AREA)
Abstract
本发明公开了一种适用于纵扭超声铣削的有限元仿真建模方法,它包括以下步骤:设定铣刀几何参数、工件材料参数、铣刀材料参数、加工工艺参数以及超声振动参数;对铣削模型进行转换,建立纵扭超声铣削三维等效模型;建立刀具、工件、工件基体及装配三维模型;对刀具、工件、工件基体划分网格;实现刀具纵扭超声振动;设定刀具切削运动边界条件;提交计算,输出切削刃运动轨迹、切削力、切削温度的计算结果。本发明通过对铣削模型进行转换,成功实现了刀具作为纵扭超声振动载体的铣削三维有限元仿真,考虑了刀具螺旋角和切削刃半径,其模型和实际加工过程更为贴切,计算精度得以保证。同时本发明通过建立工件基体模型,成果实现了切削力的输出。
Description
技术领域:
本发明涉及一种超声振动加工建模方法,特别是涉及一种适用于纵扭超声铣削的有限元仿真建模方法。
背景技术:
超声振动尤其是二维纵扭超声振动技术由于其高频的振动冲击,以及分离特性,能够有效降低加工过程中的平均切削力和切削温度,在抗疲劳制造技术中有着重要的地位。而将超声技术与铣削加工方法相复合从而获得已加工表面压应力的主动控制,是目前机械制造领域抗疲劳制造技术研究的主要方法之一。
随着计算机硬件水平的不断提高,通过有限元仿真去模拟超声铣削加工不失为一种有效的方法,一方面可以通过完善材料物理参数,以及加工中存在的物理现象,去实现加工结果的预测,相对于解析模型有更高的精度;另一方面,为试验提供指导,减少了探索试验的次数和盲区,从而降低人力物力的消耗,提高了经济性。
目前高温合金铣削有限元仿真所存在问题:
(1)二维模型:忽略材料应变厚度的影响,忽略刀具螺旋角的影响,产生较大的误差。且无法实现刀具同时产生轴向和扭转方向超声振动的仿真。
(2)三维模型:模型复杂,计算耗时长,对计算机硬件有较高要求。
(3)纵扭超声铣削三维仿真:采用传统的三维模型,无法实现刀具同时产生轴向和扭转方向超声振动的仿真。
(4)超声振动铣削切削力仿真:对于刀具作为超声振动载体,由于超声振动会产生高频周期性的加速度,采用传统从刀具参考点作为切削力的输出质点不能准确反映实际切削力。
发明内容:
本发明所要解决的技术问题是:克服现有技术的不足,提供一种模型和实际加工过程更为贴切且计算精度得以保证的新型适用于纵扭超声铣削的有限元仿真建模方法。
本发明的技术方案是:
一种适用于纵扭超声铣削的有限元仿真建模方法,包括以下步骤:
(1)设定铣刀几何参数、工件材料参数、铣刀材料参数、加工工艺参数以及超声振动参数;所述铣刀几何参数包括刀具前角r0、刀具后角a0、铣刀直径D、刀具螺旋角b、切削刃半径re、铣刀齿数N;
所述工件材料参数包括材料密度、弹性模量、泊松比、热导率、比热容、非弹性热分数、膨胀系数、准静态下的屈服强度A,硬化模量B,应变率敏感系数C,温度敏感系数m,应变硬化指数n,材料的单元失效准则参数D1、D2、D3、D4和D5;
所述铣刀材料参数包括材料密度、弹性模量刀、泊松比、热导率、比热容;所述加工工艺参数包括每齿进给量fz、切削宽度ae、进给速度vf、主轴转速n;
所述超声振动参数包括纵向和扭转超声振动频率f1和f2、纵向和扭
(2)铣削模型等效转换:
利用所述步骤(1)中的铣刀几何参数和加工工艺参数,按照图1所述的铣 削模型转换几何关系,计算得到hmax、l
其中
(3)建立刀具、工件及装配三维模型:
利用所述步骤(1)中的铣刀几何参数和加工工艺参数,结合所述步骤2中的hmax和l,按照图2所示建立刀具模型和工件模型,并对刀具、工件划分网格,刀具-工件装配的相对位置如图3所示;
(4)在工件下方增加与工件尺寸相匹配的工件基体,施加刚体约束,建立参考点,作为切削力输出的质点,如图3所示。
(5)根据所述步骤(1)中的刀具材料参数和工件材料参数,分别赋予所述步骤(3)建立的刀具和工件模型材料属性;赋予所述步骤(4)建立的工件基体工件材料属性。
(6)赋予刀具纵扭超声振动:
利用所述步骤(1)中的超声振动参数和所述步骤(3)建立的刀具-工件装配模型,在刀具参考点RP位置同时施加纵向和扭转方向上的周期性位移边界条件,实现刀具纵扭超声振动,其中扭转超声振幅 A2和铣刀等效半径Rt'满足弧长计算公式,即A2=α×Rt'/180;α为一个扭转超声振动周期内刀具切削刃所经历圆心角的弧度数;
(7)赋予刀具切削运动:
根据所述步骤(1)中的加工工艺参数,赋予刀具切削运动;
(8)提交计算及后处理:
根据所述步骤(1)-(7)所建立的模型,提交计算,输出切削刃运动轨迹、切削力、切削温度等计算结果,如图4-6所示。
本发明的有益效果是:
1、本发明通过模型转换,在刀具上建立参考点RP,相对于传统的三维模型,由于刀具不再需要随主轴产生旋转运动,不再和扭转超声振动产生过约束,因而,可将纵向超声振动和扭转超声振动同时施加于刀具参考点位置。图中LUV和TUV分别代表纵向和扭转超声振动。在施加扭转超声振动时,需要考虑刀具的等效半径(如图3中Rt'),保证切削刃位置产生符合要求的扭转振幅。
2、本发明通过在工件下方增加基体,赋予刚体约束,建立参考点,作为切削力输出的质点,成果避免了当刀具作为超声振动载体时,由于刀具受超声振动作用所产生的周期性加速度,而导致的切削力的输出不准确现象。
3、本发明模型和实际加工过程更为贴切且计算精度得以保证,易于推广实施,具有良好的经济效益。
附图说明:
图1为本发明铣削模型转换几何关系示意图;
图2为本发明刀具和工件几何模型示意图;
图3为本发明纵扭超声铣削等效模型示意图;
图4为本发明切削刃运动轨迹计算结果示意图;
图5为本发明切削力计算结果示意图;
图6为本发明切削温度计算结果示意图。
具体实施方式:
实施例:参见图1-图6。适用于纵扭超声铣削的有限元仿真建模方法,设定以下参数:铣刀前角5°,后角12°,螺旋角30°,切削刃半径0.01mm,切削速度25m/min,每齿进给量0.008mm/z,径向切深0.15mm,轴向切深2mm,纵振振幅5μm,扭转振幅4μm,纵扭频率35kHz,纵扭相位差90°,刀具材料参数如表1所示,工件材料参数如表2和3所示,工件材料单元失效准则参数如表4所示。
表1刀具材料参数
表2工件材料参数
表3工件材料本构参数
表4工件材料单元失效准则参数
根据确定的参数,计算hmax、l,如下式:
式中,fz、ae、Rt、β、vf、n、N分别代表每齿进给量、切削宽度、铣刀半径、铣刀螺旋角、进给速度、主轴转速、铣刀齿数。
根据计算结果,以及设定的参数,建立刀具和工件模型如图2。其装配后的相对位置如图3所示。
在工件下方装配与工件尺寸相匹配的工件基体,赋予刚体约束,在其中心位置建立参考点,如图3所示。
将所设定的材料参数分别赋予工件和刀具材料属性。同时赋予工件基体工件材料属性。
根据所设定的超声振动参数,赋予刀具纵扭超声振动,如图3所示,在参考点位置分别施加纵向超声振动和扭转超声振动。
提交有限元分析软件ABAQUS进行计算,输出切削刃运动轨迹、切削力、切削温度的计算结果,如图4、5、6所示。
从图4可以看到,切削刃位置的幅值、相位差以及频率满足目标期望,这也进一步说明采用该方法对纵扭超声辅助加工进行仿真是有效的。
从图5可以看到,在纵扭超声振动切削中,切削力会随着超声振动产生周期性的“回零”波动,这主要是由于纵扭超声加工过程中刀具-工件分离特性造成的,证明了采用纵扭超声振动能够有效降低平均切削力。
从图6可以看到,在切削阶段,节点处的温度逐渐上升,当刀具距离节点最近时,节点温度达到最大值,同时切削温度会随着纵扭超声振动产生周期性的波动,刀具距离节点越近,温度波动约剧烈。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (1)
1.一种适用于纵扭超声铣削的有限元仿真建模方法,包括以下步骤:
(1)设定铣刀几何参数、工件材料参数、铣刀材料参数、加工工艺参数以及超声振动参数;所述铣刀几何参数包括刀具前角r0、刀具后角a0、铣刀直径D、刀具螺旋角b、切削刃半径re、铣刀齿数N;
所述工件材料参数包括材料密度、弹性模量、泊松比、热导率、比热容、非弹性热分数、膨胀系数、准静态下的屈服强度A,硬化模量B,应变率敏感系数C,温度敏感系数m,应变硬化指数n,材料的单元失效准则参数D1、D2、D3、D4和D5;
所述铣刀材料参数包括材料密度、弹性模量刀、泊松比、热导率、比热容;所述加工工艺参数包括每齿进给量fz、切削宽度ae、进给速度vf、主轴转速n;
所述超声振动参数包括纵向和扭转超声振动频率f1和f2、纵向和扭转超声振动幅度A1和A2、纵向和扭转超声相位差j;
(2)铣削模型等效转换:
利用所述步骤(1)中的铣刀几何参数和加工工艺参数,按照铣削模型转换几何关系,计算得到hmax、l其中
(3)建立刀具、工件及装配三维模型:
利用所述步骤(1)中的铣刀几何参数和加工工艺参数,结合所述步骤(2)中的hmax和l,按照刀具和工件几何模型对刀具、工件划分网格,刀具-工件装配的相对位置构成纵扭超声铣削等效模型;
(4)在工件下方增加与工件尺寸相匹配的工件基体,施加刚体约束,建立参考点,作为切削力输出的质点;
(5)根据所述步骤(1)中的刀具材料参数和工件材料参数,分别赋予所述步骤(3)建立的刀具和工件模型材料属性;赋予所述步骤(4)建立的工件基体工件材料属性;
(6)赋予刀具纵扭超声振动:
利用所述步骤(1)中的超声振动参数和所述步骤(3)建立的刀具-工件装配模型,在刀具参考点RP位置同时施加纵向和扭转方向上的周期性位移边界条件,实现刀具纵扭超声振动,其中扭转超声振幅A2和铣刀等效半径Rt'满足弧长计算公式,即A2=α×Rt'/180;α为一个扭转超声振动周期内刀具切削刃所经历圆心角的弧度数;
(7)赋予刀具切削运动:
根据所述步骤(1)中的加工工艺参数,赋予刀具切削运动;
(8)提交计算及后处理:
根据所述步骤(1)-(7)所建立的模型,提交计算,输出切削刃运动轨迹、切削力、切削温度等计算结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811530530.1A CN110032748B (zh) | 2018-12-14 | 2018-12-14 | 适用于纵扭超声铣削的有限元仿真建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811530530.1A CN110032748B (zh) | 2018-12-14 | 2018-12-14 | 适用于纵扭超声铣削的有限元仿真建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110032748A true CN110032748A (zh) | 2019-07-19 |
CN110032748B CN110032748B (zh) | 2023-03-31 |
Family
ID=67235345
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811530530.1A Active CN110032748B (zh) | 2018-12-14 | 2018-12-14 | 适用于纵扭超声铣削的有限元仿真建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110032748B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111736529A (zh) * | 2020-06-11 | 2020-10-02 | 南京理工大学 | 基于超声铣削的碳纤维增强复合材料纤维切削角改善方法 |
CN112496680A (zh) * | 2020-11-18 | 2021-03-16 | 北京卫星制造厂有限公司 | 一种高体分铝基碳化硅螺纹孔复合加工方法 |
CN113705057A (zh) * | 2021-09-07 | 2021-11-26 | 大连理工大学 | 一种基于有限元仿真的超声振动辅助铣削力建模方法 |
CN114309744A (zh) * | 2020-09-30 | 2022-04-12 | 南京理工大学 | 一种改善机器人铣边稳定性的方法 |
CN114812486A (zh) * | 2022-05-13 | 2022-07-29 | 武汉理工大学 | 一种加工工件表面粗糙度的获取方法、装置及电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070179661A1 (en) * | 2006-01-27 | 2007-08-02 | Hideaki Onozuka | Method and program for calculating maximum depth of cut without self-excited vibration of cutting tool |
US7933679B1 (en) * | 2007-10-23 | 2011-04-26 | Cessna Aircraft Company | Method for analyzing and optimizing a machining process |
CN104484515A (zh) * | 2014-12-02 | 2015-04-01 | 华中科技大学 | 一种基于有限元的钛合金变齿距铣削三维建模方法 |
CN105631072A (zh) * | 2014-11-03 | 2016-06-01 | 南京中车浦镇城轨车辆有限责任公司 | 一种铝合金粗加工刀具设计方法 |
-
2018
- 2018-12-14 CN CN201811530530.1A patent/CN110032748B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070179661A1 (en) * | 2006-01-27 | 2007-08-02 | Hideaki Onozuka | Method and program for calculating maximum depth of cut without self-excited vibration of cutting tool |
US7933679B1 (en) * | 2007-10-23 | 2011-04-26 | Cessna Aircraft Company | Method for analyzing and optimizing a machining process |
CN105631072A (zh) * | 2014-11-03 | 2016-06-01 | 南京中车浦镇城轨车辆有限责任公司 | 一种铝合金粗加工刀具设计方法 |
CN104484515A (zh) * | 2014-12-02 | 2015-04-01 | 华中科技大学 | 一种基于有限元的钛合金变齿距铣削三维建模方法 |
Non-Patent Citations (3)
Title |
---|
封志彬等: "钛合金超声扭振铣削强度优化设计仿真", 《声学技术》 * |
王明海等: "Ti6Al4V的超声振动铣削加工三维有限元仿真", 《塑性工程学报》 * |
罗海泉等: "TC4合金径向超声振动铣削切削力特征分析", 《南昌航空大学学报(自然科学版)》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111736529A (zh) * | 2020-06-11 | 2020-10-02 | 南京理工大学 | 基于超声铣削的碳纤维增强复合材料纤维切削角改善方法 |
CN114309744A (zh) * | 2020-09-30 | 2022-04-12 | 南京理工大学 | 一种改善机器人铣边稳定性的方法 |
CN114309744B (zh) * | 2020-09-30 | 2023-07-28 | 南京理工大学 | 一种改善机器人铣边稳定性的方法 |
CN112496680A (zh) * | 2020-11-18 | 2021-03-16 | 北京卫星制造厂有限公司 | 一种高体分铝基碳化硅螺纹孔复合加工方法 |
CN112496680B (zh) * | 2020-11-18 | 2022-07-29 | 北京卫星制造厂有限公司 | 一种高体分铝基碳化硅螺纹孔复合加工方法 |
CN113705057A (zh) * | 2021-09-07 | 2021-11-26 | 大连理工大学 | 一种基于有限元仿真的超声振动辅助铣削力建模方法 |
CN114812486A (zh) * | 2022-05-13 | 2022-07-29 | 武汉理工大学 | 一种加工工件表面粗糙度的获取方法、装置及电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN110032748B (zh) | 2023-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110032748A (zh) | 适用于纵扭超声铣削的有限元仿真建模方法 | |
Song et al. | Prediction of simultaneous dynamic stability limit of time–variable parameters system in thin-walled workpiece high-speed milling processes | |
Wang et al. | Milling stability analysis with considering process damping and mode shapes of in-process thin-walled workpiece | |
Luo et al. | Effect of cutting and vibration parameters on the cutting performance of 7075-T651 aluminum alloy by ultrasonic vibration | |
CN101758253A (zh) | 减小深孔加工镗杆自激振动的方法及其动力减振镗杆 | |
CN111291479A (zh) | 混联机床铣削稳定性预测方法 | |
Chen et al. | Iterative from error prediction for side-milling of thin-walled parts | |
Guo et al. | An integrated machine-process-controller model to predict milling surface topography considering vibration suppression | |
CN112222540A (zh) | 一种螺旋锥齿轮旋转超声磨齿加工装置的设计方法 | |
CN107315864A (zh) | 基于响应曲面法的超声辅助加工的尖形刀切削力预测模型建模方法 | |
Fomin | Determining undeformed chip thickness models in milling and its verification during wood processing | |
Yuan et al. | Experimental study on a two-dimensional ultrasonic vibration platform and milling of Ti2AlNb intermetallic alloy | |
Nagesh et al. | Machine tool design with preferentially asymmetrical structures to improve dynamics and productivity | |
Schulze et al. | Simulation approach for the prediction of surface deviations caused by process-machine-interaction during broaching | |
CN104008237A (zh) | 一种计算溢流环换能器阵声辐射特性的方法 | |
Cha et al. | Stability analysis for the crankshaft grinding machine subjected to a variable-position worktable | |
Zhang et al. | Influence factors of surface topography in micro-side milling | |
Li et al. | Surface profile and milling force prediction for milling thin-walled workpiece based on equivalent 3D undeformed chip thickness model | |
Meng et al. | Theoretical and numerical investigation of micro-textures fabrication by ultrasonic surface rolling process | |
CN204997078U (zh) | 一种适于立式外圆加工的并联式椭圆振动车削装置 | |
Hadraba et al. | Virtual twin of the multi-spindle lathe for the chatter time-domain analysis | |
CN108692806A (zh) | 固体中功率超声强度分布的测量方法 | |
Li et al. | Research on prediction and compensation strategy of milling deformation error of aitanium alloy integral blisk blade | |
CN203610779U (zh) | 一种非共振三维椭圆金刚石振动切削柔性装置 | |
Meng et al. | Investigation on stability of robotic rotary ultrasonic edge milling component with poor rigidity |
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 |