基于塑性损伤理论的金属裂纹扩展模拟.docxVIP

基于损伤塑性理论的延性金属裂纹扩展模拟

江金锋1孙秦1

(1.西北工业大学航空学院，西安，710072)

摘要：为了研究延性金属裂纹扩展特点，本文基于塑性损伤理论，在ABAQUS软件平台上，

编写了材料子程序VUMAT，对一中间带孔金属板进行了裂纹萌发及扩展数值计算，模拟出

了结构中裂纹的出现及扩展过程，并提取了结构的载荷位移曲线。

关键词：损伤塑性;延性;ABAQUS;VUMAT

Ductilefracturesimulationofmetalusingdamageplasticity

theory

JIANGJinfeng1SUNQin1

(1.SchoolofAeronautics,NorthwesternPolytechnicalUniversity,Xi’an,710072)

ABSTRACT:Inordertosimulateprogressivefractureofmetal,damageplasticitytheoryistaken

intoaccount.FracturealgorithmVUAMTisthencodedusingplasticitytheoryandincorporatedin

thecommercialsoftwareABAQUS.WiththeVUMATsubroutine,numericalfracturesimulation

ofametalplatewithaholeiscarriedout,andPVcurveisextractedfromtheresult.

Keywords:damageplasticity;ductility;ABAQUS;VUMAT

材料的断裂通常被认为是工程实际中是设计师们的首选材料。

最危险的结构材料失效模式。材料一旦发生众所周知，金属材料普遍具有良好的延

断裂，后果将是不堪设想的。当今工程应用展性，因此金属材料的断裂也是一个包含了

领域中，复合材料得到越来越广泛的运用，弹塑性变形的过程，这种断裂过程被称为

但是金属材料凭借着其长时间积累下来的

“延性断裂”[1]。“延性断裂”一般是指以

使用经验在实际设计中仍占有重要地位。特微孔洞生长方式的断裂形式。但是广义上来

别是在结构中的关键受力部位的设计上，由讲，“延性断裂”可以理解成与材料的塑性

于复合材料的运用还不是很成熟，金属往往变形相关-材料中发生断裂的部位往往表现

出很大的塑性变形和良好的延展性。一个世

纪以来，广大学者在研究断裂现象的过程中

1有限元前处理[2]

1.1几何模型、边界条件及载荷作出了各种各样的努力，相继提出了应力强

度因子，能量释放率，J积分理论以及连续

本算例选取一三维带孔矩形板作为计

介质损伤力学。应力强度因子主要解决的是

算模型，具体细节如图1所示线弹性断裂问题，J积分主要运用在小变形

假设下。而断裂问题是一个高度非线性问

题，因此就必须采用损伤力学的方法来处理

金属材料断裂问题。由于塑性变形本身是不

可逆的，因此可以在塑性变形和损伤累积之

模型尺寸为100?30?5(mm)，中间圆孔

间建立关系，这种建模方法也被称为累积的

直径D?2mm，并且在板的左端面固定6个

应变损伤方法。那么，损伤累积过程就可以

方向的自由度，在板的右端面沿x方向施加

采取对塑性变形进行加权积分的形式进行，

5mm的位移载荷。

如下式所示：

D?

?????

????

c

f,,T,?d

pD

c

0

1.2材料参数

式中D为损伤变量，

?为断裂应变阈

c

值，f??,?,T,???为场变量函数，

d?为等

p

效塑性应变增量，

D为损伤阈值。当损伤

c

值D达到阈值D，认为材料断裂。本文算

c

例选取最简单的线性损伤演化方程，如下式

材料选取任意一种铝合金，在自定义材所示：

D?

?

???

c

d

?p?D

pc

0

?

c

料的参数中输入相应的弹性模量、泊松比以

及初始屈服时的屈服应力。

在ABAQUS软件平台上，将损伤演化

1.3有限元网格

方程写入用户子程序VUMAT，针对一金属

带孔板进行了无初始缺陷的断裂过程模拟。

整体模型选用C3D8R单元，采用扫掠

方式划分网格。孔边为危险区域，因此在

XY面内设定较小单元尺寸0.5mm，沿Z方

向整个模型设定统一的单元尺寸0.5mm。整

?

p

?；

体网格图如图2所示，一共包含20340