ABAQUS金属切削实例步骤

背景介绍：切削过程是一个很复杂的工艺过程，它不但涉及到弹性力学、塑性力 学、断裂力学，还有热力 学、摩擦学等。同时切削质量受到刀具形状、切屑流 动、温度分布、热流和刀具磨损等影响，切削表面的残余应力和残余应变严重影 响了工件的精度和疲劳寿命。利 用传统的解析方法，很难对切削机理进行定量 的分析和研究。 计算机技术的飞速发展使得利用有限元仿真方法来研究切削加工 过程以及各种参数之间的关系成为可能。近年来，有限元方法在切削工艺中的应 用表明， 切削工艺和切屑形成的有限元模拟对了解切削机理，提高切削质量是很 有帮助的。 这种有限元仿真方法适合于分析弹塑性大变形问题，包括分析与温度 相关的材料性能参数和很大的应变速率问题。ABAQUS 作为有限元的通用软件， 在处理这种高度非线性问题上体现了它独到的优势， 目前国际上对切削问题的研 究大都采用此软件， 因此， 下面针对 ABAQUS 的切削做一个入门的例子， 希望初 学 者能够尽快入门，当然要把切削做好，不单单是一个例子能够解决问题的，随着 深入的研究， 你会发现有很多因素影响切削的仿真的顺利进行，这个需要自己去 不 断探索， 在此本人权当抛砖引玉， 希望各位切削的大神们能够积极探讨起来， 让我们在切削仿真的探索上更加精确，更加完善。 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 切削参数：切削速度 300m/min，切削厚度 0.1mm，切削宽度 1mm 尺寸参数：本例作为入门例子，为了简化问题，假定刀具为解析刚体，因为在切 削过程中，一般我们更注重工件最终的切削质量，如应力场，温度场等，尤其是 残余应力场， 而如果是要进行刀具磨损或者涂层刀具失效的分析的话，那就要考 虑建立刀具为变形体来进行分析了。 工件就假定为一个长方形， 刀具设置前角 10°， 后角 6°， 具体尺寸见 INP 文件。 下面将切削过程按照 ABAQUS 的模块分别进行叙述，并对注意的问题作出相应的 解释。 ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 建模：建模过程其实没有什么好注意的，对于复杂的模型， 我一般用其他三维 软件导入进来，注意导入的时候尽量将格式转化为 IGES 格式，同时要把一些不 必要的东西去掉，比如一些尖角，圆角之类的，如果不是分析那 个部位的应力 集中的话就没必要导入它，如果导入，还要进行一些细化，大大降低了计算的效 率。我一般做的是二维切削，模型相对比较简单，所以一般都是直接在 ABAQUS 中进行建模。由于此处为刚体，要在 part 里面建立刚体参考点，而且注意不要 在装配模块建立参考点，因为有时候 ABAQUS 找不到装配模块相 应的参考点。 第二步定义材料参数，在这个阶段要注意的问题是要把参数定义全，尤其是你要考虑什么， 要计算什么，比如你要考虑热应力，那你就要定义热膨胀系数，你要考虑 刀具的温度场， 你就要定义刀具的密度，热导率，比热容等， 这部分最终要的是关于工件塑性的问题， 现在 切削大部分用的都是 J-C 模型，一方面它全面考虑了热 硬化，热软件等多种效应，另一方 面可以通过霍金普森实验比较容易确定， 因此广泛应用于现代切削加工的塑性模型中。 下面 针对这个模型做简单介绍。 J-C 模型 的式子我就不再输入了， 相信每一个做切削的人应该都 知道的，分为三个部分的乘积， 分别代表了热硬化，热软化和温升效应，运用这个模型的时 候，参数根据实验 输入即可，在定义的时候，有个参考应变率，这个一定记住是做实验时 候的参考应变率，所以有时候不同的参考应变率，拟合的参数A,B 等可能是不同的，因此， 在选择模型的时候一定要注意参考应变率和实验的参考应变率相对应， 这样你所得的参数才 是准确的。在定义材料参数时，要注意别忘了定义非弹性热系数，这个参 数将决定你的塑 性产热，如果你发现自己的切削温度场不对， 很大程度上是因为没有定义这个参数。 另外还 有一个非常重要的问题，单位一定要统一，不然你的计算 肯定会出错，记住 mm-t-MPa， 其他自己推一下就行。 第三部装配模块，这个没啥说的，就是把你自己工件的刀具的相对位置确定一下就行了。 先 通过 Instance part 建立工件和刀具的装配模型，然后通过translate instance调整工件和刀具 的位置，尽量保证工件和刀具充分接近， 但是不要产生干涉，工件和刀具如果如果太远，刀 具空走的路程会长，容易产生一些振动 什么的，如果太近，产生干涉，影响刀具和工件的 正确接触关系，因此要把握好工件和刀具的距离，其他的基本上都很简单。 第四部定义分析步， 目前大体在 ABAQUS 中做切削分析可以采用两种分析步， 一种是基于热 力耦合的分析步，这种分析步充分考虑了切削过程的热力相互耦合及 相应机制，另一种是 动态绝热分析，认为切削仿真是在短时间完成的， 在这么短的时间，热量还未来得及发生传 导，这种分析步一般进行锯齿切屑形状的分析， 容易产生锯齿的剪切带。如果是采用 ALE, 要定义相应的 ALE 设置。下面对 ALE 做相应简介，其实这涉及到 ABAQUS 中相应的算法，在 ABAQUS 中一 种是拉格朗日算法，一种是欧拉算法，拉格朗日算法的节点是基于实体单元 的， 也就是实体单元的节点的位置， 一般用于固体分析， 而由于算法本身的特点， 在遇到 大 变形时候很容易产生畸变， 欧拉算法是基于空间的， 一般用于流体，其相应的节点也是基于 空间的，网格划分后其在空间的位置是不变的，但是其边界条件难以定 义和琢磨，而 ALE 结合了拉格朗日算法和欧拉算法的优点， 一方面通过拉格朗日方法定义边界， 另一方面利用 欧拉方法解决网格畸变问题，因此目前在切削领域由 于很好的解决了问题而得到广泛的应 用。 关于 ALE 的具体设置可以参看 ABAQUS 经典例题 3 中的稳态切削的例子， 由于时间和空 间有限，有些东西我是说不完 的，希望大家能自己去探索。如果想输出切削力随时间的时 域图，在这里要选择参考点， 然后设置 set，然后在historyoutput里面输出参考点即 可。如 果想加快计算时间，可以采用质量放大，但是放大的尺度一定要把握好， 不然会造成结果失 真的，具体可以参照帮助文档关于质量放大的介绍。 第五部分定义接触，首先定义接触属性，这个要考虑自己具体的计算过程了， 比如考虑不考 虑摩擦等，一般对切削而言，需要设置相应的摩擦，采用硬接触的形式， 而且一定要定义 两个系数，一个是生热系数，一个是热量分配系数，生热系数主要是计算摩擦做的功有多大 的比例产生热，前面我们已经在材料里面定义了塑性产热 的转化系数，热分配系数主要是 摩擦最终产生的热量有多少流入工件， 有多少流入刀具。接下来就是定义接触了， 选择刀具 接触表面以及工件接触表面，可以建立相 应的接触关系，本例因为模型简