利用ANSYS/LS-DYNA的SPH-FEM耦合拉伸模拟

1、背景

    有限元方法作为数值计算的强大工具,计算结果精确且可重复,降低了试验成本,缩短了研发周期,但有限元方法在切削仿真时容易造成网格畸变,造成求解中断。

    光滑粒子动力学(smoothed particle hydrodynamics,SPH)的基本思想是将连续体离散为相互作用的粒子,每个粒子具有密度、质量以及相关物理属性,粒子间运动遵循牛顿第二定律;其本质是一种拉格朗日方法,运用插值理论将宏观变量(如压力、密度以及温度等)一系列无序点的值通过微分形式转换成积分运算。SPH法采用粒子划分,不依赖于网格,具有很好的自适应性,可以避免网格畸变,适合切削引起的大变形问题。然而SPH法对每个粒子进行计算时,需要搜索影响区域内近邻的粒子信息、粒子物理量计算和搜索信息都比较费时,因此计算效率比普通的有限元法低,对于三维模型占用计算机资源较大。

    针对SPH与FEM的各自特点,为提高计算效率并消除网格畸变,采用SPH与FEM耦合的方法解决切削数值模拟问题。在变形大的区域采用SPH,避免FEM的网格畸变过大造成计算困难。在变形小的区域采用FEM,以提高计算效率。SPH与FEM耦合算法分为固定耦合算法和自适应耦合算法。固定耦合算法在计算之前就已确定SPH区域和FEM区域。自适应耦合算法则在计算之前都是FEM网格,在计算过程中自动地将大变形的有限元网格单元转换为光滑粒子,并按SPH法计算物理量。

    基于以上考量,本文运用ANSYS/LS-DYNA进行了SPH-FEM耦合算法的拉伸试验模拟。

2、模型设置

    分析模型如下图所示,拉伸件两端采用壳单元,中间段采用SPH粒子法划分。粒子与壳单元接触段采用tie功能进行绑定,以实现FEM与SPH之间的耦合计算。

1.png

由于采用了耦合算法,还需要对壳单元和SPH粒子进行相关的设置,具体内容如下:

2.png
3.png

    对于模型的材料设置,考虑到模型的形状,断裂破坏肯定会发生在中间粒子区域,而模型的两端壳单元区域属于加载区域,不会发生破坏,也不是本次模拟的关心区域,因此为了进一步提高求解效率和节约求解资源,模型将壳单元区域赋予刚体材料模型,即不考虑模型两端的变形情况。粒子区域的具体材料参数如下图所示:

4.png

    为模拟拉伸工况,本次模拟中将模型的一端壳单元的自由度全部约束,使其成为固定端,在另一端壳单元采用线性位移加载,加载曲线如下图所示:

5.png

    除此之外,还需要设置相关的输出,计算终止时间等内容,在此不进行一一赘述。模型摄制完成之后即可导出K文件,利用ANSYS/LS-DYNA求解器进行求解。

3、结果分析

6.png
7.png
8.png
9.png
10.png
11.png
12.png
13.png

14.png

    以上为拉伸件的塑性应变随时间的分布图,可以看出断裂发生在预期位置,证明了采用SPH-FEM耦合方法进行联合仿真是可行的。SPH-FEM耦合的方法,吸收了FEM法计算效率高和SPH法模拟大变形能力强的优点,可以为大变形的材料仿真如切削等提供一种高效、准确的途径。

SPH-FEM耦合切削仿真
为TA投票

利用ANSYS/LS-DYNA的SPH-FEM耦合拉伸模拟的评论2条

利用ANSYS/LS-DYNA的SPH-FEM耦合拉伸模拟的相关案例教程

高速冲击涉及材料的大变形、破碎和飞溅等现象,应用基于网格的方法对其进行数值模拟存在困难。拉格朗日网格方法会遭遇单元畸变而计算终止。光滑粒子流体动力学(SPH)作为一种无网格、拉格朗日粒子法,能克服基于网格的方法的缺陷。SPH在处理大变形方面较有限元法(FEM)等拉格朗日网格方法有优势,计算精度和效率都不及FEM,并且SPH的边界处理不如FEM方便。 基于此,发展了将SPH与FEM进行耦合的方法,有
LS-DYNA中,通过FEM-SPH方法可较好地模拟爆破飞石的效果,该方法可用于研究碎块的分布趋势,飞溅速度、位移等运动轨迹,对于损伤的块度、粒度分析则存在缺陷。 与自适应方法不同,炸药与部分岩石区域为单独的粒子算法,周围岩石为拉格朗日单元算法,能节省大量的计算时间。粒子与粒子之间的作用不需要额外定义接触,但是对于粒子的密度及粒子的个数需要合理控制,较多粒子数量可得到更好的爆破效果,同时会大幅提高
岩石在爆破荷载作用下表现的性质十分复杂,在实际的钻爆开挖中容易引发安全事故。埋入岩体中的炸药发生化学爆炸后,岩石将在爆破近区发生压剪破坏形成粉碎区,在爆破远区发生张拉破坏形成裂隙区。 由于在冲击荷载下存在有限元大变形的问题,因此SPH-FEM耦合法常被用于模拟爆炸、侵彻等问题。下面给出SPH-FEM模拟爆破的结果,并附有k文件供参考学习。 炮孔附近的岩石及炸药采用SPH粒子,远端采用有限元。FEM
传统的有限元方法在求解流固耦合问题时存在许多困难,而FEM-SPH(有限元-光滑粒子法)在求解流固耦合问题时可以完美解决这个问题,FEM-SPH耦合算法可以作为一种新的思路求解流固耦合问题。本案例中采用FE-SPH耦合算法有效地模拟了高突水问题对露天台阶的影响,露天台阶采用FEM有限元模型,高位水库采用SPH粒子,可以实现流体与固体系统的动态耦合分析。模拟结果较好地反映了突水对台阶的冲击造成的动力
内容简介 该案例以药柱在混凝土内部爆炸为例,讲解如何采用SPH_FEM耦合算法实现药柱爆炸对混凝土损伤的数值模拟。该案例主要内容如下: (1)如何建立SPH_FEM爆炸模型, (2)SPH相关控制关键字如何设置, (3)如何实现SPH和FEM之间的耦合, (4)如何控制不同药柱的起爆时间, (5)如何查看混凝土的损伤参数。
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    2 11