huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析

Johnson-cook 损伤起始准则是延性损伤准则模型的一个特例,用于预测延性金属中孔洞的形核、生长和聚结导致的损伤起始。该模型假设损伤开始时的等效塑性应变是应力三轴性和应变率的函数。同时可以考虑温度的影响。

包含的材料参数有:

失效相关参数:d1-d5。

金属的熔化温度: θ melt

金属的转变温度: θ transition

参考应变率:ε0

当满足下列条件时,损伤初始化准则得以满足:

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的图1

等效塑性应变认为与应力三轴度和应变率相关联。

θ^是无量纲温度,表示为:

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的图2

其中,θ是当前温度,θ-melt是熔化温度,θ_transition是指转变温度,在该温度或低于该温度时,损伤应变εpl_D的表达式不存在温度依赖性。材料参数必须在转变温度或低于转变温度时测量。

损伤的发展可以公式化为:

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的图3

公式中分母表示单元失效对应的Johnson-cook等效塑性应变,公式为:

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的图4

分子表示为等效塑性应变增量,公式为:

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的图5

公式中可以看到,损伤随着塑性应变的增大不断累积,直至材料的失效,通过损伤变量进一步与晶体材料的屈服面或者弹性性能的退化可以实现材料弹-塑-损伤的耦合模拟,当不对其进行耦合时,可以用来判断材料的失效状态与相关参数的关系。

参考文献:《Crystal plasticity finite element modeling and simulation of diamond cutting of polycrystalline copper》编写对应的材料子程序。在huang晶体塑性程序的基础上,调用johnson-cookd损伤函数,编写过程中,需要自定义响应的状态变量,如等效塑性应变,等效塑性应变率,损伤变量,以及是否进行损伤单元的删除分析。其中等效塑性应变增量的计算,通过滑移系统的分切应力与对应滑移系统剪切应变的乘积绝对值之后与等效应力的比值获得。并最终实现损伤的表征,采用umat子程序进行编写。

umat子程序对应的材料参数为

弹性:C11=168400Mpa , C12=121400Mpa , 75400Mpa

塑性:n=10 ,γ_dot0=0.001 , h0=541500Mpa , Tau_s=109500Mpa , tau_0=60800Mpa

损伤:d1-d5 ( 0.54 , 4.59 , 3.03 , 0.014 , 1.12 ) ,

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的图6

可以看出损伤分布于强剪切带处,因为该区域具有最大的塑性应变。同时由于采用的损伤为考虑对弹塑性的影响,因此模型的计算效率与原始模型保持一致



θ melt=1083






ABAQUS二次开发

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的评论3条

  • 沉鱼海
    1
    想咨询一下大佬,我们用huang的子程序是对晶体的弹塑性力学变化进行分析,这时候如果考虑损伤,是不是应该侧重的是晶体层面的破坏,比如说沿晶断裂或者穿晶断裂,如果用J-C损伤模型,这是不是相当于用宏观的经验公式来代替实际的机理了?
  • JHChen
    0
    厉害了👍👍👍大佬,有偿求教程

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的相关视频课程

huang晶体塑性umat耦合Johnson-cook 损伤模型,实现晶体材料弹-塑-损伤模拟分析的相关案例教程

在运用JC本构模型的时候,不知如何控制损伤开始的位置,在学习总结之后分享出来,希望和大家一起进步。 JC本构模型包括塑性硬化段和损伤演化段 1 JC本构——塑性硬化段 方程: 式中:A,B,n,m 是控制塑性段硬化的材料参数,等号右侧第二个括号与第三个括号分别是应变率和温度对于塑性硬化段的影响。 关于第二个括号: \dot{\varepsilon}_{\mathbf{0}}:参考应变率,这个值越小
Johnson-Cook 材料模型及失效模型。一般用于描述大应变(large strains)、高应变率(high strain rates)、高温(high temperatures)环境下金属材料的强度极限以及失效过程。在Johnson-Cook强度模型中,屈服应力(yield stress)由应变、应变率以及温度决定。 屈服应力的表达式为: 其中,A,B,N,M是材料参数,epsilon_p
(原创,欢迎转载,转载请说明出处) 1 概述 本系列文章研究成熟的有限元理论基础及在商用有限元软件的实现方式,通过 (1) 基础理论 (2) 商软操作 (3) 自编程序 三者结合的方式将复杂繁琐的结构有限元理论通过简单直观的方式展现出来,同时深层次的学习有限元理论和商业软件的内部实现原理。 有限元的理论发展了几十年已经相当成熟,商用有限元软件同样也是采用这些成熟的有限元理论,只是在实际应用过程中,
一、前言 无人机的结构由机身、动力系统、飞行控制系统、链路系统、任务载荷等几个方面组成的。 无人机动力系统中的电机,俗称“马达”,是无人机的动力来源,无人机通过改变电机的转速来改变无人机的飞行状态。即改变每个电机的速度,使得无人机能够盘旋空中,上升或下降,或向各个方向移动。 电机又分为"无刷电机"和"有刷电机"。现在市面上大部分无人机都是用无刷电机,无刷电机是采用半导体开关器件来实现电子换向,具有
umat子程序只用于通用静力分析步,使用用户自定义材料,这里对横观各向同性材料本构模型进行构建。 对于横观各向同性的刚度矩阵为: 由矩阵可以知道,横观各向同性材料仅需5个独立常数,用横观各向同性的5个参数进行表示,进而可以构建出材料的本构方程。不懂去翻材料力学相关书籍。 这里就随便建立一个高20,半径为1的圆柱体对其进行施加一个单位为5的拉伸位移(较长方向拉伸)。 材料属性就随便找了一篇论文,按照
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    4 13