目前，国内对于6005A铝合金大型薄壁中空结构件的高速数控加工存在诸多问题，如数控加工效率较低、刀具损耗大、铣削速度较高时出现崩刃、打筋等，这些问题严重影响了企业生产效益。因此，研究6005A铝合金的高速铣削加工机理，以解决生产中出现的问题，具有重要的现实意义。

在铣削加工过程中，切削力决定切削热的产生，并影响刀具使用寿命、加工精度和加工表面质量，研究切削力的规律和计算方法有助于分析铣削机理。解析法和实验法是研究切削力的主要方法，但随着计算机技术的发展，有限元数值模拟方法为研究切削过程中的切削力提供了新的方法。国内外学者在二维、三维正交切削数值模拟方面做了很多的工作，取得了丰硕的成果，然而针对6005A铝合金的高速数控铣削加工的研究还未见报道。因此本课题采用三维斜角切削数值模拟的方法，对6005A铝合金高速铣削过程中的切削力变化规律进行探讨，为优选和优化高速铣削6005A铝合金结构件的刀具几何参数提供依据。

6005A铝合金材料特性

6005A合金是在6005合金的基础上发展的一种中等强度Al-Mg-Si系合金，具有优良的挤压性、焊接性和耐腐蚀性，广泛应用于地铁和巴士车体结构件、铁塔平台、管道设施等。6005A铝合金的主要化学成分和材料基本性能参数如表1和表2所示。

有限元模型建立

1 热力耦合有限元控制方程

在金属切削加工时，随着刀具的进给，材料应力不断增加。当应力超出弹性范围时，材料发生塑性变形。塑性变形功和刀屑之间的摩擦功转变成热而引起温度的升高，从而由于温度的升高产生热应变。在不考虑弹性蠕变的情况下，弹塑性区域的总应变分量dε 可以表示为：

dε=dεel+dεpl+dεT  （1）

式中，dεel 为弹性应变增量；dεpl为塑性应变增量；dεT为热应变增量。

基于虚功原理和大应变、大变形理论，在不考虑体力的情况下，建立热-力耦合弹塑性控制方程

Kep+ [KG]d˙=∫[Bε]T Depεt } dV−[Bε]T {R˙ εT} dV +˙F （2）

式中，[Kep] 为弹塑性刚度矩阵；[KG]为几何刚度矩阵；d为节点速度；{R˙ εT}为相对温度和等效应变率的变化率；˙F为外部载荷的变化率，包括实际载荷变化率和由变形体形状改变引起的载荷变化率。

2 几何模型

车辆铝合金6005A是易切削但容易产生加工缺陷的材料，一般采用2或3齿的螺旋立铣刀进行高速数控铣削加工。立铣刀是双切削刃刀具，圆柱面上的切削刃是主切削刃，端面上的切削刃是副切削刃，一般情况下主、副切削刃同时参与材料去除。根据文献，铣削加工过程可以看作是沿切削刃许多微小斜角切削过程的综合作用。因此，本课题采用斜角切削模型对切削力进行分析，更有助于铣削加工过程中刀具参数对切削力影响规律的理解，同时也简化了有限元数值计算模型。利用有限元软件本身所拥有的3D建模功能，分别建立刀具和工件的几何模型，并在装配模块下定义其装配关系。

3 材料模型

材料模型的正确定义是金属切削过程数值模拟成功的一个关键因素。有限元软件用一系列直线段连接给定的真实应力- 应变数据对来形成连续的、分段线性的塑性曲线，可以用任意多的数据对来逼近实际的材料行为，从而可得到近似于真实材料的塑性行为。图1为通过拉伸试验得到的6005A铝合金的工程应力应变曲线。

通过拉伸试验得到工程应力-应变数据，需要转换为真实应力-应变。真实应力-应变和工程应力-应变的转换关系为：

σ = σnom（1+εnom） （3）

ε = ln（1+εnom） （4）

式中，σ为真实应力； σnom为名义应力；ε为真实应变；εnom为名义应变。#p#副标题#e#

试验数据中得出的应变为总体应变，因此将总体应变分解为弹性应变和塑性应变分量。弹性应变等于真实应力与杨氏模量的比值，则塑性应变等于总体应变减去弹性应变。表3为经过数据整理计算后得到的6005A铝合金的真实应力-应变数据。在有限元软件中，输入材料的塑性参数时，初始真实应变必须为0。通过在ABAQUS有限元软件包Property模块中选择Mechanical→ Plasticity→Plastic，输入表3中所给出的6005A铝合金的真实塑性应力-应变数据，来描述实际的材料行为。

4 刀具切屑接触摩擦特性

金属切削过程即刀具与工件相互运动并相互作用，使切削层金属与工件母体金属分离开来。切削层金属在刀具刃口和前刀面的挤压、摩擦作用下发生剪切滑移变形和摩擦变形而形成切屑。前刀面和切屑的接触区域分为滑移区和黏结区。在黏结区，切削时产生的高温和高压使得切屑底层材料软化，切屑底层的金属材料黏嵌在前刀面上，实际上是金属内部的剪切滑移，刀屑接触点处的摩擦剪应力等于临界剪切流动应力。在滑移区服从库仑摩擦定律。刀屑摩擦特性表示为：

τ =μσn (σn < τpμ)

τp (σn τpμ) （ 5）

式中，τ为摩擦应力；τp为材料的最大剪切应力；σn为前刀面正应力；μ为摩擦系数。

但是滑移区的摩擦系数是要通过分析刀具前面的应力分布情况计算出来，因此本课题采取平均摩擦系数的方法，把摩擦系数设置为0.3，最大剪应力设为205MPa。

5 网格划分

工件和刀具均采用8节点六面体热力耦合线性减缩积分单元C3D8RT进行网格划分，采用增强沙漏控制。为了提高数值模拟精度和减少计算机时，只对切屑层和刀尖处进行了网格细化。因为刀具和工件的硬度相差很大，所以分析时视刀具为刚体，刀刃绝对锋利。工件共划分8280个单元，刀具共划分480个单元，有限元模型如图2所示，其中工件尺寸为5mm×1mm×1mm。

数值模拟结果及分析

利用ABAQUS有限元软件对车辆铝合金材料6005A高速铣削加工进行了数值模拟，模拟时刀具从图2所示位置开始切削，沿X方向做进给运动，边界条件如图2所示，约束工件底面6个自由度，刀具除X方向移动自由度外，其余全部受到约束。刀具材料选择硬质合金立铣刀，牌号为K20，弹性模量为600GPa，泊松比为0.25。模拟采用的切削参数: 进给速度为3m/min，切削深度为0.2mm，切削宽度为1mm。

在切削参数一定的情况下，模拟不同刀具前角和刃倾角情况下的斜角切削过程，分析刀具几何角度变化对切削力的影响。表4给出了2组模拟条件。在高速切削的条件下，刀具前角对切削过程的影响更大，因此在第一组的模拟条件中，变化刀具前角，后角和刃倾角为目前正在使用的刀具的参数一致。在第二组的模拟条件中，在高速铣削时，为了减小刀具——工件之间的摩擦，综合考虑可以将后角选为12°以上，本课题选择13°进行模拟，分析刃倾角对切削力的影响。