断裂失效!LS-DYNA材料及结构Gissmo失效模型仿真分析如何做?

文摘   2024-11-30 09:01   广东  

点击文尾阅读原文观看

导读:随着经济全球化的快速发展,产品的开发周期越来越短,特别是汽车行业。之前一辆的仿真轮次有几轮,一轮有一两个月,现在轮次大幅度压缩,每轮也极限缩短,有些产品甚至压缩至一个月之内。如何在这么卷的时间内做好CAE,做对CAE,如何高效的做CAE,成为我们从事CAE的头等问题。


某整车正面碰撞图

CAE在各行各业的进程中应用愈加广泛,尤其是汽车行业,如主机厂、零部件企业、新能源企业等。随着CAE技术的快速发展,CAE在产品开中的应用场景也越来越多,传统的一些仿真模拟方法的精度不能满足要求,如碰撞仿真、冲击、挤压、动态球击等工况,在极限载荷下结构件容易出现失效断裂,这种断裂在有些情况下是可以接受的,也是必须的(如前正碰撞过程副车架也车身连接后螺栓需要剪断,进而保护乘员安全);但很多情况是不允许出现断裂失效,否则会出现安全问题。

图2 某正碰过程中副车架螺栓断裂图

图3 某新能源动力电池边框挤压失效图

全球钢铁行业工程研究显示,新的先进高强度钢与现代工程设计和计算机辅助工程(CAE)工具的结合使用可以减少25%的车辆零部件。根据不同的碰撞载荷要求,新型的钢材发挥着重要作用。如通常使用双相或TRIP钢来确保在前部和后部压溃区结构中吸收高能量,而其他高强度钢则用于抵御侧面冲击和车顶的压溃,进而实现性能及重量的平衡。
图4 钢材的性能发展演变图

一、当前主流的失效类型

如何在当前更加准确的模拟出断裂情形,方法有很多种,如发展出来的各种断裂失效模型(如JC失效模型、GISSMO失效模型、MMC模型等)也在各行各业的产品开发发挥着重要作用。

目前常用的断裂失效准则既包括:
① 最大主应力失效
② 最大等效塑性应变失效等简单常值失效判据
③ 成形极限曲线 ( Forming Limit Curve,FLC) 或成形极限图(Form-ing Limit Diagram,FLD)
④ Johnson-Cook损伤模型 ( JC 损伤模型)
⑤ 损伤起始与演化准则( Damage Initiation and Evolution Criteria,DIEC)
⑥ 修正的 Mohr-Coulomb 模型 ( Modified Mohr-Coulomb,MMC)
⑦ CrachFEM模型

⑧ Gissmo失效模型等。

其中使用相对比较广泛的为1)、2)、4)和8)等,其中1)和2)为常应力或应变失效,适用于恒定载荷加载;而Johnson-Cook 断裂准则以线性方式计算损伤积累,对于高速冲击或碰撞精度相对一般;Gissmo失效模型考虑了非线性损伤累积方式和等效应力与损伤之间的耦合,适合复杂工况下的断裂失效表征。当 GISSMO 损伤模型中的断裂极限曲线取为MMC 或 JC 模型断裂极限曲线,损伤累积方式为线性,且不考虑材料软化、失稳时,GISSMO 损伤模型自动退化为 MMC 或 JC 损伤模型。

在Lsdyna中对于金属材料常见的断裂失效模型包括:

(a)固定临界值:如最大/最小压力、最大/最小应力、最大等效失效应变(也称为常应变失效)/应力、最大体积应变、最大剪切应变、成形极限图(FLD)等;

(b)失效模型:Johnson-Cook、Gissmo、DIEM、MMC等,如JC模型可考虑温度、应力三轴度、断裂应变、应变率等效应;Gissmo可考虑应力三轴度、断裂应变、洛德角、网格尺寸等之间的关系。

二、GISSMO失效模型简介 
Gissmo一种唯象损伤力学模型,以非线性损伤累积的方式描述材料从变形到破坏的整个过程,而不用追究损伤的物理背景和材料内部微观结构变化。由于同时考虑了材料在不同应力状态下的失效应变以及应变路径对失效的影响和非线性损伤积累方式,适用于超高强钢板材在复杂工况下的断裂行为表征。
Gissmo涉及的相关理论非常多,涉及的相关参数也非常丰富,如应力三轴度、失效应变、不稳定性应变、损伤累积指数、应力退化指数、软化指数、中面失效、应变率、网格尺寸等。

Gissmo失效模型是基于材料开发的一种失效模型,材料的力学性能测试及转换至关重要,如测试得到的工程应力应变如何转换成真实应力应变,以及塑性应变等;同时得到的数据如何通过硬化公式外推等等这些都是最基础,也是最重要的。

以下是工程应力应变与真实应力应变的转换公式,如下所示。

1、工程与真实应力的关系

(1)工程应力

指在拉伸或压缩试验过程中,施加在样件的载荷除以样件的初始横截面积:

式中σe为工程应力,A0为样件初始截面积。

(2)真实应力

指在拉伸或压缩试验过程中,施加在样件的载荷除以实际截面积。

式中σt为真实应力,A为样件当前实际截面积。

(3)真实应力与工程应力的关系

根据样件拉伸过程中塑性阶段的体积不变假设,即

2、工程与真实应变的关系

1)工程应变

指样件标距范围内的伸长量除以标距的初始值:

式中εe为工程应变,∆l为样件的长度变化量。

5 GISSMO失效模型主要参数图

三、GISSMO失效模型的关键参数
图6 某样件准静态拉伸失效对比图

GISSMO主要的失效相关方程及参数有以下:

(1)损伤的定义

金属韧性材料的断裂过程实质上是微孔洞的成核、长大、聚合及开裂的过程,此时的损伤累积D会对材料的本构关系产生重要影响,其定义如下:

图7 损伤定义图

假设总面积S对应的应力为σ∗,有效承载面积S ̂的应力为σ,那么:

σ∗≤σ,进一步可得:σ∗=σ1- D)

由上式可知,在样件拉伸过程中,样件有效截面积S不断减小;当有效承载面积S0时,此时的D1,即损伤达到1时,样件断裂。

(2)损伤增量公式的定义

GISSMO模型中,损伤的表达式如下所示:

损伤增量的表达式如下:

对损伤增量进行积分后,可得总的损伤:

在软件中采用的损伤增量表达式如下:

上式中表示损伤增量,表示不同应力状态下的失效应变,n为损伤累积指数,表示真实的等效塑性应变增量。
当损伤D1时,材料失效,单元删除。

(3)不稳定性增量公式的定义

当材料开始出现颈缩时,表现出材料的不稳定性,采用类似损伤增量的关系式来表征材料的不稳定性,如下式所示。即当稳定性因子F1时,损伤与应力开始耦合,常见的材料不稳定性主要表现为扩散颈缩和局部颈缩。

上式中表示稳定性变量增量,表示不同应力状态下的材料出现不稳定性变形时的局部等效塑性应变,n为损伤累积指数,表示真实的等效塑性应变增量。

(4)应力耦合公式的定义

随着外力的继续加载,损伤不断增大,而宏观应力不断减小,当损伤为1时,材料失效,裂纹产生,损伤与应力耦合的关系图如下所示。从图中可以看出当材料达到临界应变后,即F=1,此时的应变即为GISSMO中的ECRIT,此时的损伤将会储存到DCRIT中;在此之后,损伤不断增大,宏观应力表现出明显的下降,刚度开始退化,当损伤达到1时,材料断裂,宏观应力为0,宏观应力与损伤的耦合表达式如下所示。

图8 损伤与应力耦合关系图

9 损伤与应力耦合结果图
图10 GISSMO中LCREGD结果图

三、Gissmo失效模型工程案例应用实操

GISSMO涉及的相关参数非常多,其主要的三个公式如上所示,分别为损伤增量公式、不稳定性公式及损伤与应力耦合公式。笔者通过总结梳理得到Gissmo专题课程《LS-DYNA材料及结构Gissmo失效模型工程案例应用实操48讲(附学习模型)》原创独家首发仿真秀。感兴趣的小伙伴请关注本套课程,课程目录初步如下,请识别下方二维码可以试看。
限时特价(限20名)
券后价1599元价值:1999 
课程可随时回放,可开具发票,奖学金、直播加餐
讲师提供vip群答疑和模型下载,开具电子发票

 扫码立即试看

1 本课适合人群:

(1)学习仿真工程师;

(2)Ls-dyna软件的学习者和使用者;

(3)碰撞、强度及材料等相关仿真工程师;

(4)理工科院校学生及老师;

(5)对失效仿真感兴趣的朋友;

(6)汽车及机械等相关行业的从业者。

2 你会得到什么:

(1)掌握Ls-dyna常用材料的特性及使用方法;

(2)掌握金属材料力学性能的方法和技巧

(3)掌握当前主流失效模型类型

(4)掌握MAT_ADD_EROSION关键字的使用方法和技巧

(5)掌握GISSMO失效模型相关基础

(6)掌握GISSMO失效模型的理论基础

(7)掌握金属材料拉伸数据处理及拟合方法、技巧和流程

(8)掌握GISSMO失效模型理论之应力三轴度

(9)掌握GISSMO失效模型理论之罗德角

(10)掌握GISSMO失效模型理论之应变率

(11)掌握GISSMO失效模型理论之损伤累积计算等

(12)掌握GISSMO失效各个关键字调用方法

(13)掌握GISSMO失效各个关键字含义及使用方法

(14)掌握GISSMO失效标定技巧及分析流程

(15)掌握GISSMO失效模型相关案例和经验积累

(16)为订阅用户提供VIP群、课程资料、不定期加餐答疑、高薪内推和电子发票等服务。

注:购买课程的可以加入专属学习交流群,同时根据需要可提供相关的学习模型,提供的学习模型仅供购买者自己使用,本课程不提供课件!整理的相关模型及干货学习资料如下。参考:笔记总结、Lsdyna帮助文档及相关文献等。
课程配套资料

点击“阅读原文”了解更多

纵横CAE
技术交流,知识共享。
 最新文章