利用LS/DYNA中的SPH法进行旋转刀具切削模拟

1、工程背景

传统有限元方法在切削仿真中应用已较为广泛,其在理论基础、边界处理、计算效率等方面发展均已较为成熟。但由于有限元方法是基于单元网格划分,在处理大变形问题时,网格极易发生畸变,导致计算崩溃。虽然有研究人员使用网格二次划分等方法应对该问题,但这些方法也往往导致计算时间大量延长等问题。针对这些情况,本文在前人研究基础上,采用一种SPH无网格法进行切削仿真。

SPH方法是基于插值理论的方法,是一种朗格朗日性质的无网格法。在对工程材料进行离散化处理时,它与有限元存在重大区别:对材料的离散化不使用单元,而是使用固定质量的质点或粒子,对某一系统的状态及其运动情况则通过这一系列彼此间相互作用的粒子来描述。质点的质量固定在坐标系上,所需的基本方程也是守恒方程和固体材料的本构方程。按照运算原理来说,SPH是把物理流场用具有一定流动速度的运动质点来描述,每个质点就是已知流场特性的插值点,这样整个问题的解可以通过质点的规则插值函数得到。SPH法不依赖网格,从而允许材料界面可以“天然”的存在,适用于解决工程材料在高速加载速率下的断裂等难题,能够较为简捷而精确地实现材料复杂的木构行为,因此如今已经广泛应用在水下爆炸仿真、高速碰撞等动态响应的数值模拟领域。

LS-DYNA是功能完备的大型商用计算软件,主要用来进行非线性显式分析。该软件适用于的求解类型包含以下几种:各种几何非线性问题、材料非线性与接触非线性问题、多种高度非线性的瞬态动力学问题。LS-DYNA提供的算法以Lagrange算法为主,但同时也拥有ALE算法和Euler算法。该软件主要以处理非线性动力分析为主,但同时也具有静力学分析的能力。LS-DYNA针对的工程领域以结构分析问题为主,但同时提供了热力学分析以及流-固耦合计算功能。

 2、模型介绍

模型分为两部分考虑:对于切削层区域,切削过程中会材料发生剧烈变形形成切屑,因此采用SPH粒子划分。很多时候工程问题中某些构件刚度远远大于其余部件,工作过程中可以认为是不发生弹性变形的。基于这一理论,ANSYS/LS-DYNA允许用户把有限元模型中刚度相对较大的构件定义为刚体,有利于大大减少显式分析所用时间。因此本文中刀具设置成刚性体,不考虑其在切削过程中的变形,采用有限元网格划分。具体模型如下图所示:

1.png
2.png

建立完几何模型之后需要进行相关材料本构的定义,本文刀具选用金刚石,并将其设置为刚体,不考虑骑在切削过程中的变形,具体材料参数如下图:

3.png

对于切削层的脆性材料,本文选用当前应用较为广泛的JH模型,具体参数如下:

4.png

除此之外,选用合适的接触算法也有助于提高计算精度,本文选用的是点-面接触,使用该接触类型时,接触节点与接触面之间会发生穿透,为解决这一问题,LS-DYNA提出的方案是定义contact表面和target表面,因此这种接触类型适用于一般情况下的两个物体表面之间的接触。LS-DYNA中这种接触类型主要包括NODES TO SURFACE,AUTOMATIC NODES TO SURFACE等。本文设置如下:

5.png

3、结果分析

6.png
7.png
8.png
9.png
10.png
12.png
13.png
15.png
16.png
18.png
19.png

21.png

对切削过程的模拟来说合理的流动应力模型是获得理想结果的关键。由图可知,在切削初始阶段,随着刀尖与工件接触,切削层材料产生很大的接触应力,形成一个局部高应力区,以工件与刀尖的接触点为中心,且与刀具的移动方向成一定角度。随着刀具的继续切入,等效应力场逐渐向外扩大,同时随着最大应力达到工件材料的屈服极限,切削层材料发生不可逆变形,开始沿着前刀面塑性流动而形成切屑。与实验方法相比,切削仿真的一大优点是可以获得切削过程中材料内部等效应力的变化情况,为切屑工件分离提供分析依据。

切削模拟SPH法
为TA投票

利用LS/DYNA中的SPH法进行旋转刀具切削模拟的评论2条

利用LS/DYNA中的SPH法进行旋转刀具切削模拟的相关案例教程

(I)节点穿透问题解决方案 摘录国外某论坛(亲测有效) You may also want to consider defining a contact thickness for your sph particles, as they have no associated contact thickness. (Use ‘SST’ parameter in the *CONTACT card),
1、背景 有限元方法作为数值计算的强大工具,计算结果精确且可重复,降低了试验成本,缩短了研发周期,但有限元方法在切削仿真时容易造成网格畸变,造成求解中断。 光滑粒子动力学(smoothed particle hydrodynamics,SPH)的基本思想是将连续体离散为相互作用的粒子,每个粒子具有密度、质量以及相关物理属性,粒子间运动遵循牛顿第二定律;其本质是一种拉格朗日方法,运用插值理论将宏观变
LS-DYNA中的显式SPH求解功能非常适合求解涉及超高速撞击、爆炸和其他瞬态事件等问题,但在涉及诸如涉水等较慢的流体流动仿真时仍需优化。在此基础之上,不可压缩SPH (ISPH)功能是专门为处理诸如涉水、电机冷却、齿轮润滑等大型不可压缩流体仿真而开发,它允许比通常的显式SPH仿真更大的时间步长,同时避免了对流体不可压缩性的妥协。与显式SPH和其他FVM方法相比,ISPH方法所需的仿真计算时间更少
传统的有限元方法在求解流固耦合问题时存在许多困难,而FEM-SPH(有限元-光滑粒子法)在求解流固耦合问题时可以完美解决这个问题,FEM-SPH耦合算法可以作为一种新的思路求解流固耦合问题。本案例中采用FE-SPH耦合算法有效地模拟了高突水问题对露天台阶的影响,露天台阶采用FEM有限元模型,高位水库采用SPH粒子,可以实现流体与固体系统的动态耦合分析。模拟结果较好地反映了突水对台阶的冲击造成的动力
LS-DYNA中,通过FEM-SPH方法可较好地模拟爆破飞石的效果,该方法可用于研究碎块的分布趋势,飞溅速度、位移等运动轨迹,对于损伤的块度、粒度分析则存在缺陷。 与自适应方法不同,炸药与部分岩石区域为单独的粒子算法,周围岩石为拉格朗日单元算法,能节省大量的计算时间。粒子与粒子之间的作用不需要额外定义接触,但是对于粒子的密度及粒子的个数需要合理控制,较多粒子数量可得到更好的爆破效果,同时会大幅提高
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    2 19