有限元分析方法和材料断裂准则.docx

一、有限元模拟方法金属切削数值模拟常用到两种方法，欧拉方法和拉格朗日方法欧拉方法适合在一个 可以控制的体积内描述流体变形，这种方法的有限元网格描述的是空间域的，覆盖了可以 控制的体积在金属切削过程中，切屑形状的形成过程不是固定的，采用欧拉方法要不断 的调整网格来修改边界条件，因此用欧拉方法进行动态的切削过程模拟比较困难欧拉方 法适用于切削过程的稳态分析（即“Euler方法的模拟是在切削达到稳定状态后进行的”［2］） 仿真分析之前要通过实验的方法给定切屑的几何形状和剪切角［1］而拉格朗日方法是描述固体的方法，有限元网格由材料单元组成，这些网格依附在材 料上并且准确的描述了分析物体的几何形状，它们随着加工过程的变化而变化这种方法 在描述材料的无约束流动时是很方便的，有限元网格精确的描述了材料的变形情况实际 金属切削加工仿真中广泛采用的拉格朗日方法，它可以模拟从初始切削一直到稳态的过 程，能够预测切屑的形状和工件的残余应力等参数⑵但是用这种方法预定义分离准则和 切屑分离线来实现切屑和工件的分离，当物质发生大变形时常常使网格纠缠，轻则严重影 响了单元近似精度，重则使计算中止或者引起严重的局部变形［1］。

为了克服欧拉描述和拉格朗日描述各自的缺点，Noh和Hirt在研究有限差分法时提出 了 ALE（Arbitrary Lagrange-Euler）描述，后来又被 Hughes，liu 和 Belytschko 等人引入到有 限元中来其基本思想是：计算网格不再固定，也不依附于流体质点，而是可以相对于坐 标系做任意运动由于这种描述既包含Lagrange的观点，可应用于带自由液面的流动，也 包括了 Euler观点，克服了纯Lagrange方法常见的网格畸变不如意之处自20世纪80年 代中期以来，ALE描述己被广泛用来研究带自由液面的流体晃动问题、固体材料的大变形 问题、流固祸合问题等等金属的高速切削过程是一个大变形、高应变率的热力祸合过程， 正适合采用ALE方法采用ALE方法进行高速切削仿真克服了拉格朗日方法和欧拉方法需要预先定义分离 线、切屑和工件分离准则，假定切屑形状等缺点，避免了网格畸变以及网格再划分等问题， 使切屑和工件保持良好的接触，使计算易于收敛［1］囹二、材料断裂准则在金属切削成形有限元模拟中提出了多种切屑分离准则，这些准则可以分为两种类 型：物理准则和几何准则优点：几何分离准则需要预定义加工路径，在加工路径上判断刀尖与刀尖前单元节点的距离 变化来判断分离与否。

当两点的距离小于某个临界值时，刀尖前单元的节点被分成两个， 其中一个节点沿前刀面向上移动形成切屑，另一个保留在加工表面上形成己加工表面 ［1］［2］物理分离准则是基于刀尖前单元节点的应力、应变及应变能等物理量定义分离条件， 当单元中的该物理量的值超过给定材料的对应值时，单元节点就会分离［2］物理标准主要 是基于制定的一些物理量的值是否达到临界值而进行判断的，主要有基于等效塑性应变准 则、基于应变能密度准则、断裂应力准则等［5］）Carroll和Strenkowski使用了等效塑性应变作为物理分离准则的标准，在一些有限元 软件中该标准的演化得到了应用，ABAQUS/Explicit中的剪切失效准则（shear failure）就 是这样一种物理准则，它根据单元积分点处的等效塑性应变值是否到达预设值来判断材料 是否失效［1］缺点：采用物理准则使切削的有限元模拟更接近实际情况，但在实际的有限元模拟中，当刀 尖达到应该分离的节点时，该点的物理值并没有达到所给定的物理标准，导致切屑在该点 并没有分离［5］很难在有限软件中实现［3］几何标准的模型虽然简单，但是由于它不是基于切屑分离的物理条件，所以使用几何 标准很难找到一种通用的临界值，以适应切削加工中不同的材料以及不同的加工工艺［2］。

难以反映切屑分离过程中的力学和物理现象⑶为了克服物理准则和几何准则的缺点，又有学者提出几何-应力切屑分离准则该准则 以物理准则为主要判断依据，但当刀尖接近分离点并小于给定的几何值时，可以强迫结点 分离［2］标准如下［3］： f = Jb /b 2 + T卢2Johnson-Cook剪切失效准则提供了材料到达失效点时等效塑性应变的计算方法，特别 适用于金属的高应变率变形，因此本研究采用了 Johnson-Cook剪切失效准则作为刀屑分离 的物理准则，同时采用了 ALE方法进行刀屑分离［1］三、材料断裂模型在有限元模拟时，切屑形状的成形除了受刀屑分离准则的影响，同时受材料断裂模型 的影响，尤其是单元切屑和锯齿状切削的形成由材料的断裂模型来实现Elbestawi和 El-Wardanylsl应用断裂力学理论对高速切削中工件材料裂纹萌生和扩展方向进行了预测， 指出当工件材料自由表面的能量达到某一临界值时，裂纹开始产生，并沿着应变能密度最 小的方向不断扩展，进而造成材料的断裂目前用来进行高速切削有限元模拟的断裂准则有Rice&Tracy准则、Brozz准则、 Cockroft&Latham准则、MeClintock准则和Freudenthal准则。

大量研究者采用 Cockroft&Latham准则，该准则是从能量角度建立的，通过高温拉伸试验计算出断裂塑性 能，并同金属材料变形断裂所需的能量建立映射关系，将其作为判断金属材料延性断裂的 临界能量值考虑到高速切削的特点，塑性应变对工件材料失效断裂起重要影响，ABAQUS/Explicit 里定义材料属性时可在材料编辑（Edit Material）选用Johnson-Cook断裂应变模型，该模型（见 公式（3.7）、（3.8））提供了材料到达失效点时等效塑性应变的计算方法［1］本研究选用了 Johnson-Cook强度模型为材料的本构模型，选用了 Johnson-Cook剪切 失效准则为刀屑分离准则，同样采用了 Johnson-Cook断裂应变模型为材料的断裂模型［1］ABAQUS Explicit 里的 adaptive mesh采用 了 ALE技术剪切失效准则和 Johnson-Cook 断裂应变模型一起应用就能动态的判断材料的失效并达到工件与切屑的分离［1］参考文献［1］ 蒋志涛.高速金属铣削硕士学位论文,2009,2［2］ 金属切削过程有限元建模.［3］ 吴红兵，贾志欣，刘刚，毕运波，董辉跃.航空钛合金高速切削有限元建模.浙江大学学报（工学 版），2010,（44） 5.［4］ 成群林，柯映林，董辉跃，杨勇.高速硬加工中切屑成形的有限元模拟.浙江大学学报（工学版）， 2007 （41） 3 .［5］ 苏桂生.高速铣削中基于正交切削模型的有限元仿真.苏州大学硕士学位论文,2008,4。