梯度理论的应用

基于应变梯度理论的微切削毛刺仿真研究

原创 工具技术 工具技术 2019-12-30

作为复杂微小零件的主要加工方式，微铣削技术具有加工效率高、精度高、范围广的优点，是微细加工领域的研究热点，具有很大的发展空间。由于微铣削所加工的零件特征尺寸小，其表面形成的毛刺尺寸相对较大，会严重影响零件的使用性能和寿命。因此对微铣削过程中产生的毛刺进行深入研究具有重要意义和应用价值。

近年来，国内外学者对微铣削过程中产生的毛刺进行了研究。1973年，美国学者L.K.Gillespie等最早提出了毛刺这一概念，并且指出加工时工件材料发生剧烈的塑性变形导致了毛刺产生；Schueler G.M.等基于毛刺产生机理，将毛刺分为泊松毛刺、翻转毛刺、撕裂毛刺和剪断毛刺四种类型。微铣削是一个非常复杂的过程，根据刀具参数、切削参数以及所采用工件材料的不同，加工过程中所产生的毛刺会有不同的形态；Deng W.J.等通过建立基于热力耦合的二维平面应变正交切削模型对毛刺的形成进行研究，得到出口毛刺的影响因素有工件出口角度、刀具前角、每齿进给量以及切削刃圆弧半径；Guo Y.B.等通过建立三维有限元模型对钻削进行仿真，发现毛刺的形成过程需经历开始、发展、工件材料弯曲以及毛刺最终形成四个阶段；Özel T.等分别用未涂层刀具和CBN涂层刀具对Ti-6Al-4V材料进行试验和有限元仿真，发现CBN涂层刀具的加工表面质量更好，毛刺尺寸更小，刀具磨损量更小。

目前，对微铣削毛刺的研究主要使用有限元仿真法和试验法，少数学者则通过研究毛刺产生机理，建立毛刺尺寸的数学模型。与铣削相比，微铣削不仅尺寸上减小，而且具有尺寸效应等特性。上述对微铣削毛刺进行的有限元仿真研究均使用传统的Johnson-Cook材料本构模型，未考虑材料的尺寸效应，因此本文应用基于位错机制的应变梯度理论修正的Johnson-Cook材料本构模型进行二维正交切削仿真，分析工件材料尺寸效应对出口毛刺的影响。

1  基于应变梯度理论的无氧铜本构模型

金属材料进行有限元仿真分析的基础是建立合适的材料本构模型。因为金属材料切削是一个大d塑性变形、高应变率和大温度变化的复杂过程，所以与之相对的材料本构方程也较为复杂。许多学者通过大量的试验与仿真对比，发现由Johnson和Cook提出的材料本构模型能够很好地模拟金属切削过程，即

式中，σJC为材料流动应力；ε为应变；为应变率；0为参考应变率(1s-1)；T为温度；T0为室温(20℃)；Tmelt为材料熔点；A、B、n、C、m为材料本构参数，可由试验测得。

无氧铜的J-C本构模型参数见表1。无氧铜的其他相关参数见表2。

表1  J-C材料模型参数

上述Johnson-Cook材料本构模型不能反映微切削过程中的材料尺寸效应，应用基于位错机制的应变梯度理论对传统的Johnson-Cook材料本构模型进行修正，得到反映工件材料尺寸效应的本构关系式为

式中，G为材料的剪切模量；b为Burgers矢量；μ为修正系数；α为材料系数；L为主剪切变形区长度。

从式(2)可知，修正的JC本构方程流动应力的求解主要取决于主剪切变形区长度L的求解。在微切削中，根据切削深度与刀具圆弧半径比值的不同，会出现形成切屑和不形成切屑两种状态。当有切屑形成时，主剪切变形区长度L为

式中，h为未变形切削厚度（μm）；φ为剪切角（°）。

通过研究发现剪切角的表达式为

式中，t为常量，刀具材料为硬质合金，工件材料为无氧铜，此时t取值为47°；β为刀具和工件接触的摩擦角（°）；αe为刀具前角（°）。

当没有切屑形成时，主剪切变形区长度L为刀具和工件接触的弧长，其表达式为

式中，R为刀具切削刃圆弧半径（μm）。

上述修正的材料Johnson-Cook本构方程同时具有传统Johnson-Cook本构方程和能够描述材料尺寸效应的特点，因此适用于微切削有限元仿真。在切削仿真时，需要设定材料的断裂准则来模拟切屑形成。断裂时的等效塑性应变为

式中，d1-d5为无氧铜失效常数，可由试验测得，分别为0.54、4.89、-3.03、0.014、1.12；εf为断裂时等效塑性应变；σ*为应力三轴度。

在切削仿真过程中，刀具使工件不断发生塑性变形，当塑性应变积累到一定值时，该部分材料失效，单元从网格中删除，从而实现切屑与工件的分离。该过程表达式为

式中，Δε为等效塑性应变增量；D为失效参数，当D达到1时，单元失效。

2  微切削有限元仿真模型的建立

采用ABAQUS软件对微切削进行建模和仿真。刀具材料选用硬质合金，刀具前角为5°，后角为10°，切削刃钝圆半径为4μm，工件材料选用无氧铜，仿真中刀具材料和工件材料属性见表3。

表3  工件和刀具材料属性

由于本次仿真不考虑刀具磨损情况，因此将刀具设置为刚体。为减少仿真时间，切削区域内的网格划分较细密，非切削区域内的网格划分较稀疏，划分后的网格情况见图1。

图1  二维切削网格划分模型

3  材料本构模型的二次开发与验证

ABAQUS有限元软件中有很多材料本构模型，如考虑应变梯度，需要用户自编vumat子程序。使用Fortran语言根据上述修正的材料本构模型编写子程序vumat，利用ABAQUS调用调试成功的vumat子程序进行仿真。分别使用软件自带的Johnson-Cook材料本构模型和修正的Johnson-Cook材料本构模型进行切削仿真，通过对比仿真结果验证修正的本构模型是否正确。仿真时，切削深度为4μm，切削速度为45m/min，切削条件为干切削。图2为不考虑材料尺寸效应的应力云分布，图3为考虑材料尺寸效应的应力云分布。

图2  不考虑材料尺寸效应的应力云分布

图3  考虑材料尺寸效应的应力云分布

通过比较图2和图3的仿真应力云图可以得出：二者仿真所得的应力场分布大致相同，由于考虑材料尺寸效应的原因，图3的最大应力值（708.0MPa）大于图2的最大应力值（474.0MPa）。为了进一步研究材料尺寸效应对剪切区内最大应力的影响，分别在切削厚度为0.5μm、2μm、4μm、10μm、50μm、100μm的条件下进行切削仿真，测量各切削厚度下同一位置的最大应力值（见图4）。

由图4可以看出：考虑材料尺寸效应时，随着切削厚度的减小，剪切区内最大应力值明显增加；不考虑材料尺寸效应时，随着切削厚度的变化，最大应力值基本上保持不变。同时还可以看出：随着切削厚度的增加，两种模型最大应力的差值不断减小。说明随着切削厚度的增加，考虑材料尺寸效应的本构模型逐渐向传统的本构模型靠近。因此，考虑材料尺寸效应的本构模型既适用于微切削仿真，又适用于传统切削仿真。

图4  切削厚度与最大应力的关系

4  二维正交切削出口毛刺仿真

为了验证材料尺寸效应对出口毛刺的影响，对微切削过程进行二维正交切削仿真。切削深度为4μm，切削速度为45m/min，切削条件为干切削。图5为不考虑材料尺寸效应的仿真云图，图6为考虑材料尺寸效应的仿真云图。

图5  不考虑材料尺寸效应的仿真云图

通过图5和图6可以看出：考虑材料尺寸效应时的出口毛刺高度（10.576μm）大于不考虑材料尺寸效应时的出口毛刺高度（8.285μm）；考虑材料尺寸效应时的出口毛刺宽度（14.013μm）大于不考虑材料尺寸效应时的毛刺宽度（7.909μm）。

图6  考虑材料尺寸效应的仿真云图

为了进一步研究材料尺寸效应对出口毛刺高度和宽度的影响，分别在切削厚度为0.5μm、2μm、4μm、7μm、10μm的条件下进行切削仿真，分别测量出口毛刺高度值（见图7）和宽度值（见图8）。

由图7和图8可以得出：随着切削厚度的增加，出口毛刺高度和厚度均增加；考虑材料尺寸效应时的出口毛刺高度和宽度分别大于不考虑材料尺寸效应时的出口毛刺高度和宽度。

图7  切削厚度与出口毛刺高度的关系

图8  切削厚度与出口毛刺宽度的关系

为了分析上述两种情况出口毛刺高度和宽度不同的原因，在切削速度为45m/min、干切削条件下，分别测量切削厚度为0.5μm、2μm、4μm、7μm、10μm时的水平切削力（见图9）。由图可知，随着切削厚度的增加，水平方向上的切削力会增加；并且在同一切削厚度条件下，考虑材料尺寸效应时的切削力大于不考虑材料尺寸效应时的切削力。这是因为在考虑材料尺寸效应时，出口处负剪切区域形成更早，从而导致出口毛刺高度和宽度更大。

图9  切削厚度与水平切削力的关系

小结

（1）基于材料应变梯度理论建立的本构模型既可以描述切削过程中大应变、高应变率、高温度的特性，又可以描述微切削过程中材料尺寸效应的特性，因此该模型适用于微切削仿真。

（2）考虑材料尺寸效应时，随着切削厚度的增加，剪切变形区内的最大应力值会减小；当切削厚度增加到一定值时，最大应力值和使用传统本构模型仿真的最大应力值趋于相同。

（3）在同一切削厚度条件下，考虑材料尺寸效应的出口毛刺高度值和宽度值均大于不考虑材料尺寸效应的对应值。这主要是因为在考虑材料尺寸效应的条件下，水平方向上的切削力更大，导致出口处塑性变形严重，更早地形成负剪切区，从而产生更大的出口毛刺高度和宽度。

1、首先,先将链接里的vs2010打开安装文件,点击setup,显示出安装界面,点击next。

2、紧接着点击我接受许可,点击next。

3、选择安装目录,选择完成后点击install,等待导入自动导入,导入完成后,点击下一步,然后点击finish即可。

4、首先将上面的intel fortran安装包打开,点击开始安装,d出下面界面时点击下一步,然后选中我接受许可,点击下一步。

5、在下面d出的窗口选中使用许可文件,然后点击下一步,然后点击browse将安装包里的.lic文件导入,点击下一步。

6、这样，就将fortran软件集成到Visual Studio上,然后就可以在Visual Studio上进行编写fortran程序了。

安装gfortran：sudo apt-get install gfortran

其他fortran编译器像intel fortran，g77等也按同样的方法安装，没有源就自己下载个包安装。terminal下输入命令 g77 或者 gfortran 或ifort查看安装是否成功。

linux下的编译器一般没有界面，用记事本之类的编写文件，然后编译就好。如果是gfortran的话在命令行下gfortran g.f(Fortran文件名)；如果是intel fortran的话 ifort g.f.前提是你要在命令行下先进入你文件的文件夹下，以便编译器找到文件。

---