基于接口的huang显式晶体塑性计算与隐式计算结果对比

昨日展示了三维的对比情况,因此这里对比针对二维情况,考虑拉伸和剪切变形(即考虑不同的应力状态和单元类型)

模型包含500个晶粒,60000个单元,使用平面应变三节点单元(CPE3)

对比指标:等效应力分布,累计剪切应变分布,滑移系统当前强度分布

结果如下(默认左侧为显式结果,右侧为对应的隐式结果):

拉伸情况:

等效力分布;

基于接口的huang显式晶体塑性计算与隐式计算结果对比的图1

累计剪切应变分布:

基于接口的huang显式晶体塑性计算与隐式计算结果对比的图2

滑移系统当前强度分布:

基于接口的huang显式晶体塑性计算与隐式计算结果对比的图3

剪切情况:

等效力分布;

基于接口的huang显式晶体塑性计算与隐式计算结果对比的图4

累计剪切应变分布:

基于接口的huang显式晶体塑性计算与隐式计算结果对比的图5

滑移系统当前强度分布:

基于接口的huang显式晶体塑性计算与隐式计算结果对比的图6

可以看到使用隐式计算结果与显式计算结果几乎一致,然而显式的优势是显而易见的,尤其是在模拟高速冲击以及其他类似的接触问题

模拟多晶冲击的视频效果如下(隐式计算无法收敛,而显式可以轻松完成)

模拟效果可以在公众号查看

ABAQUS二次开发abaqus晶体塑性碰撞冲击显式动力学及仿真

基于接口的huang显式晶体塑性计算与隐式计算结果对比的评论1条

  • 晶体塑性有限元
    0
    欢迎关注新创建的微信订阅号,nzy17613019646,本订阅号将分享关于微细观塑性损伤等方面的精品文章,以及有限元建模分析的技巧和方法,重点关注GTN模型相关理论,晶体塑性相关理论,VPSC模型相关理论,以及abaqus,ansys相关的学习经验,心得。希望感兴趣的小伙伴共同交流学习,一起进步

基于接口的huang显式晶体塑性计算与隐式计算结果对比的相关案例教程

基于HUANG晶体塑性板材冲压成型模拟------案例十四 案例教学如下 1,分别建立板料半径80mm的1/4圆环,厚度0.65,夹具和冲头模型并装配 冲压的模型 2,分配材料属性:板材采用晶体塑性本构,夹具冲头采用纯弹性属性,并且在冲压过程形状中形状保持不变(约束成刚体) 3,建立接触条件,建立板料与上下夹具,以及冲头的接触 接触属性的建立 4,建立合适的约束条件,夹具完全固定,冲头施加Z方向1
黄永刚院士编写的单晶晶体塑性UMAT,主要用于在Abaqus有限元仿真中进行单晶及多晶晶体塑性变形的计算,是许多科研工作者学习晶体塑性模拟的教学资源。可以在其基础上对硬化模型进行修改,甚至引入损伤。 UMAT主要应用于隐式分析,而对于大变形接触问题,隐式分析往往计算效率较低。对于接触、碰撞、冲击等问题采用VUMAT往往具有更高的计算效率和收敛速度。本文旨在将Huang编写的UMAT改写为VUMAT
基于HUANG晶体塑性模型的三点弯曲案例 案例教学如下 1,建立几何模型:0.6*0.3*0.05(mm)的平板,半径0.05,长度0.3的下压棍和支撑并完成装配 三点弯曲模型图 2,分配材料属性:分别赋予板材晶体属性参数,棍子纯弹性参数,并假设下压辊和支撑的刚度远远大于平板,即变形过程中约束为刚体。 模型的材料属性分配 3,建立接触条件:板材与棍子支撑之间法向硬接触,切向摩擦系数设置为0.2 建
黄永刚原始晶体塑性具有良好的收敛性,以及高效的计算效率,在一般变形下无需修改,即可直接使用。然而一些特殊的工况,如切削,轧制,冲压等隐式存在收敛性问题。因此通常使用显示程序进行计算。但从头完成显式晶体塑性构造对于一般学者显然难度过高,一个简单的想法就是直接将现成的黄永刚隐式程序改成显式。abaqus里这是可以实现的。其基本的步骤是: 1,加入vumat接口程序(见附录abaqus官网有) 2,对n
基于凌志matlab脚本实现拉伸压缩,平面应变压缩过程中的织构演化 具体操作过程如下: 一,在for文件中加入输出内容如下: 代码: FILE1='C:\Users\Desktop\VPSC\cp-huang-al\texture2\temp.dat' OPEN(91,FILE=FILE1,STATUS='UNKNOWN',ACCESS='APPEND') 位置如图: 选择特定增量步格式化输出结果
影响力
粉丝
内容
获赞
收藏
    项目客服
    培训客服
    1 1