ABAQUS平台的扩展有限元方法模拟裂纹实现

ABAQUS平台的扩展有限元方法模拟裂纹实现

1.1 扩展有限元方法（XFEM）在ABAQUS上的实现

ABAQUS中XFEM的实现，两个步骤最为关键：

1、选择模型中可能出现的裂纹区域，将其单元设为具有扩展有限元性质的enrichment element.

2、其次重要的是选择恰当的破坏准则，使单元在达到给定的条件破坏，裂纹扩展。

在ABAQUS中模拟裂纹扩展的操作中，需要注意的是：

1、在Property模块，添加损伤演化参数、破坏法则、损伤稳定性参数

2、在Interaction模块，主菜单Special中创建XFEM的enrichment element

对于固定的裂纹模型，采用ABAQUS/STANDARD中使用奇异渐进函数。针对移动的裂纹问题，在XFEM中，有一种方法基于traction-separation cohesive behavior，即使用虚拟节点连续片段法进行移动裂纹建模，ABAQUS/STANDAR D 中用于计算脆性或韧性材料的裂纹初始化和扩展过程的模拟。另外一种cohesive segments method (粘性片段方法)可用于bulk material中的任意路径的裂纹初始化模拟扩展过程，由于裂纹扩展不依赖于单元边界，在XFEM中，裂纹每扩展一次需要通过一个完整单元，避免尖端应力奇异性。除此之外，ABAQUS为拥护提供了自定义子程序，来满足不同建模的需要。ABAQUS/STANDARD中的任意力学本构模型均可用来模拟扩展裂纹的力学特性。

由于XFEM采用的形函数在求解过程中，很容易造成逼近线性相关，极大的增加了收敛难度，到目前为止，能够实现扩展有限元的商业软件只有ABAQUS，但是ABAQUS为了减少求解难度，做了大量简化，因此用ABAQUS 扩展有限元模拟裂纹扩展时，有一些局限[16]：

1.扩展单元内不能同时存在两条裂纹，所以ABAQUS不能模拟分叉裂

纹；

2.在裂纹扩展分析过程中，每一个增量步的裂纹转角不允许超过90度；

3.自适应的网格是不被支持的；

4.固定裂纹中，只有各向同性材料的裂纹尖端渐进场才被考虑。

1.2 数值算例

1.2.1中心裂纹平板问题求解

给定一矩形板，长3m ，高6m,在板中心处有一初始直裂纹，并且受到拉伸应力的作用，左右两边施加载荷为1KN ，材料参数为：

在Property 功能模块中，我们采用最大主应力失效准（Maxps Damage ）则作为裂纹损伤开始的判据,损伤演化选取基于能量的、线性软化、混合模式的损伤演化规律，相关参数见上表。裂纹的扩展过程如下：

图1.1为裂纹第10步的扩展

QUANTITY

NAME

EXPRESSION Young's modulus E

210E9

Poisson's ratio

0.3 Max principle stress S

84.4E6 Fracture energy C G 1 42200 First direction C G 2 42200

Second direction

C G 3

42200

图1.2为裂纹第26步的扩展

图1.3为裂纹第56步扩展

图1.4为裂纹扩展的最终结果

Mises应力分布云图如下：

图1.5为Mises应力分布云图

图1.6为最大最小应力矢量分布云图

图1.7为初始裂纹左端点位移随时间的变化

1.2.2 边缘裂纹实体求解

给定一实体，长6m,高3m,宽1m,在边缘处有一初始直裂纹，前后两面施加载荷为10KN ，给定材料参数为：

同样，我们在Property 功能模块中，我们采用最大主应力失效准（Maxps Damage ）则作为裂纹损伤开始的判据,损伤演化选取基于能量的、线性软化、混合模式的损伤演化规律，相关参数见上表。计算结果如下：

图1.8为变形后的网格图

QUANTITY

NAME

EXPRESSION Young's modulus E

300E9

Poisson's ratio

0.2 Max principle stress S

84.5E6 Fracture energy C G 1 20580 First direction C G 2 20580

Second direction

C G 3

20580

图1.9为Mises应力分布云图

图1.10为变量的矢量图

图1.11为初始裂纹端点应力随时间的变化

图1.11为初始裂纹端点位移随时间的变化

图1.12为XFEM裂纹图像显示

1.3 本章小结

本章内容介绍了扩展有限元法在ABAQUS上实现的操作的关键步骤，需要注意的地方，也指出了ABAQUS为了实现XFEM的，进行了一些简化而造成了

某些局限。并且通过两个数值算例直观的反映了裂纹的扩展过程，包括裂纹随时间增量步长的显示分析结果的动画以及Mises应力的分布云图等其他参量的变化图。实现了本论文的最终目的即扩展有限元的ABAQUS实现。

ABAQUS平台的扩展有限元方法模拟裂纹实现

ABAQUS平台的扩展有限元方法模拟裂纹实现 1.1 扩展有限元方法（XFEM）在ABAQUS上的实现 ABAQUS中XFEM的实现，两个步骤最为关键： 1、选择模型中可能出现的裂纹区域，将其单元设为具有扩展有限元性质的enrichment element. 2、其次重要的是选择恰当的破坏准则，使单元在达到给定的条件破坏，裂纹扩展。 在ABAQUS中模拟裂纹扩展的操作中，需要注意的是： 1、在Property模块，添加损伤演化参数、破坏法则、损伤稳定性参数 2、在Interaction模块，主菜单Special中创建XFEM的enrichment element 对于固定的裂纹模型，采用ABAQUS/STANDARD中使用奇异渐进函数。针对移动的裂纹问题，在XFEM中，有一种方法基于traction-separation cohesive behavior，即使用虚拟节点连续片段法进行移动裂纹建模，ABAQUS/STANDAR D 中用于计算脆性或韧性材料的裂纹初始化和扩展过程的模拟。另外一种cohesive segments method (粘性片段方法)可用于bulk material中的任意路径的裂纹初始化模拟扩展过程，由于裂纹扩展不依赖于单元边界，在XFEM中，裂纹每扩展一次需要通过一个完整单元，避免尖端应力奇异性。除此之外，ABAQUS为拥护提供了自定义子程序，来满足不同建模的需要。ABAQUS/STANDARD中的任意力学本构模型均可用来模拟扩展裂纹的力学特性。 由于XFEM采用的形函数在求解过程中，很容易造成逼近线性相关，极大的增加了收敛难度，到目前为止，能够实现扩展有限元的商业软件只有ABAQUS，但是ABAQUS为了减少求解难度，做了大量简化，因此用ABAQUS 扩展有限元模拟裂纹扩展时，有一些局限[16]： 1.扩展单元内不能同时存在两条裂纹，所以ABAQUS不能模拟分叉裂 纹； 2.在裂纹扩展分析过程中，每一个增量步的裂纹转角不允许超过90度； 3.自适应的网格是不被支持的； 4.固定裂纹中，只有各向同性材料的裂纹尖端渐进场才被考虑。 1.2 数值算例

扩展有限元简介

扩展有限元 有限元是将一个物理实体模型离散成一组有限的相互连接的单元组合体, 该方法在考虑物体内部存在缺陷时间，单元边界与几何界面一致，会造成局部网格加密，其余区域稀疏的非均匀网格分布，在网格单元中最小的尺寸会增加计算成本，再者裂纹的扩展路径必须预先给定只能沿着单元边界发展。 1999年，美国西北大学Beleytachko 提出了扩展有限法，该方法是对传统有限元法进行了重大改进。扩展有限元法的核心思想是用扩充带有不连续性质的形函数来代表计算区域内的间断，在计算过程中，不连续场的描述完全独立于网格边界，在处理断裂问题有较好的优越性。利用扩展有限元，可以方便的模拟裂纹的任意路径，还可以模拟带有孔洞和夹杂的非均质材料。 扩展有限元是以标准有限元的理论为框架，保留传统有限元的优点，目前商业软件中如Abaqus 等都加入扩展有限元的分析模块。 扩展有限元以有限元为基本框架，主要针对不连续问题进行研究，相对于传统有限元方法，它克服了裂纹扩展问题的不足。其采用节点扩展函数，其中包括2个函数:裂纹尖端附近渐进函数表示裂纹尖端附近的应力奇异性;间断函数表示裂纹面处位移跳跃性。整体划分位移函数表示为 αααI =I I I =∑∑++=b x F a x H u x N x u N i )(])()[()('41 1 式中:)(x N I 为常用的节点位移函数；I u 为常规形状函数节点自由度，适用于模型中的所有节点；)(x H 为沿裂纹面间断跳跃函数；I a 为节点扩展自由度向量，这项只对形函数被裂纹切开的单元节点有效；)(x F α为裂纹尖端应力渐进函数；αI b 为节点扩展自由度向量，这项只对形函数被裂纹尖端切开的单元节点有效。 沿裂纹面间断跳跃函数)(x H 表达式为： otherwise n x x if x H 0)(11)(*≥-???-= 式中:x 为样本点；*x 距x 最近点；n 为单位外法线向量。 各向同性材料的裂纹尖端渐进函数)(x F α表达式为： ????? ?=2cos sin ,2sin sin ,2cos ,2sin )(θθθθθθαr r r r x F 裂纹尖端的渐进函数并不局限于各向同性弹性材料的裂纹建模。可用于弹塑性指数硬化材料，不同的裂纹尖端渐进函数的形式与裂纹位置、非线性材料变形程度有关。

abaqus扩展有限元(xfem)例子(裂缝发展) ()

Abaqus扩展有限元(XFEM)例子（裂缝发展） part模块中的操作： 1. 生成一个新的part，取名为plate，本part选取3D deformable solid extrusion类型（如图1） 2.通过Rectangle工具画出一长3，高6的矩形。考虑使用工具栏add-dimension和edit dimension来画出精确长度的模型。强烈建议此矩形的左上角坐标为（0，3），右下角坐标为（3，-3）（如图2） 3. 完成后拉伸此矩形，深度为1.（如图3） 4. 生成一个新的part，取名为crack，本part选取3D deformable shell extrusion类型（如图4）

5.生成一条线，此线的左端点坐标为（0，0.08），右端点坐标为（1.5，0.08） 6 . 完成后拉伸此线，深度为1.（如图6） 7.保存此模型为XFEMtutor（如图7），以后经常保存模型，不再累述。 8. 在part Plate中分别创建4个集合，分别为：all，bottom，top和fixZ，各部分的内容如图

8～11所示 Material模块中的操作： 1 创建材料elsa，其弹性参数为E=210GPa，泊松比为0.3（如图12） 2 最大主应力失效准则作为损伤起始的判据，最大主应力为84.4MPa（如图13）

3.损伤演化选取基于能量的、线性软化的、混合模式的指数损伤演化规律，有关参数为G1C= G2C= G3C=42200N/m，a=1.（如图14） 4.创建一个Solid Homogeneous 的section，名为solid（如图15），此section与材料elsa相

ABAQUS版本XFEM扩展元例子的详细图解

A B A Q U S6.9版本X F E M(扩展有限元)例子的详细图解 一、part模块中的操作： 二、 1.生成一个新的part，取名为plate，本part选取3Ddeformablesolidextrusion类型（如图1） 三、 2.通过Rectangle工具画出一长3，高6的矩形。考虑使用工具栏add-dimension和editdimension 来画出精确长度的模型。强烈建议此矩形的左上角坐标为（0，3），右下角坐标为（3，-3）（如图2） 四、 3.完成后拉伸此矩形，深度为1.（如图3） 五、图1，图2，图3， 4.生成一个新的part，取名为crack，本part选取3Ddeformableshellextrusion类型（如图4）

5.生成一条线，此线的左端点坐标为（0，0.08），右端点坐标为（1.5，0.08） 6.完成后拉伸此线，深度为1.（如图6） 7.保存此模型为XFEMtutor（如图7），以后经常保存模型，不再累述。 8.在partPlate中分别创建4个集合，分别为：all，bottom，top和fixZ，各部分的内容如图8～11所示。

二、Material模块中的操作： 1.创建材料elsa，其弹性参数为E=210GPa，泊松比为0.3（如图12） 最大主应力失效准则作为损伤起始的判据，最大主应力为84.4MPa（如图13）

损伤演化选取基于能量的、线性软化的、混合模式的指数损伤演化规律，有关参数为 G1C=G2C=G3C=42200N/m， =1.（如图14） 2.创建一个SolidHomogeneous的section，名为solid（如图15），此section与材料elsa相联（如图16），并将此section赋给platepart（也就是集合all）（如图17） 3.赋予材料取向，分别如图18～21所示。

abaqus中xfem扩展有限元教程

abaqus 中xfem扩展有限 元教程

part 模块中的操作: 1. 生成一个新的 part ，取名为 plate ，本 part 选取 3D deformable solid extrusion 类型 （如图1） 2. 通过Rectangle 工具画出一长 3，高6的矩形。考虑使用工具栏 add-dimension 和edit dimension 来画出精确长度的模型。强烈建议此矩形的左上角坐标为（ 0, 3），右下角坐标 为（3，-3）（如图2） 3.完成后拉伸此矩形，深度为 1.（如图3 ）3D '、2D Planw I ' Axisymmetric Tyre Options ” Di scr^te ri gi >1 f'■ Analytical ri

4.生成一个新的 part ，取名为 crack ，本 part 选取 3D deformable shell extrusion 类型 (如图4) 叩刊网扌 rr Ack M ud-el L iLg. Spa-j-e (*) 3D { ' 29 Pl war ( ) Ajci symmetri c Typ? @ H 栏 ￡oir.ahle： :；Di 5?r ?tc ari gi d Cj An>lytic41 rigid ■.. j Euler i an Opti QKS None available Q hl 迥 ⑥*1】 ■：.\ Wire (.Poiitt Base Feature Type Planar Ez trusi on Rezolution Swsep Xppr^MiTatt =it e Cine el 5.生成一条线，此线的左端点坐标为( 0, 0.08 ),右端点坐标为(1.5 , 0.08 )

裂纹扩展的扩展有限元(xfem)模拟实例详解

基于ABAQUS 扩展有限元的裂纹模拟 化工过程机械622080706010 李建 1 引言 1.1 ABAQUS 断裂力学问题模拟方法 在abaqus中求解断裂问题有两种方法（途径）：一种是基于经典断裂力学的模型；一种是基于损伤力学的模型。 断裂力学模型就是基于线弹性断裂力学及其基础上发展的弹塑性断裂力学等。如果不考虑裂纹的扩展，abaqus可采用seam型裂纹来分析(也可以不建seam，如notch型裂纹)，这就是基于断裂力学的方法。这种方法可以计算裂纹的应力强度因子，J积分及T-应力等。 损伤力学模型是指基于损伤力学发展而来的方法，单元在达到失效的条件后，刚度不断折减，并可能达到完全失效，最后形成断裂带。这两个模型是为解决不同的问题而提出来的，当然他们所处理的问题也有交叉的地方。 1.2 ABAQUS 裂纹扩展数值模拟方法 考虑模拟裂纹扩展，目前abaqus有两种技术：一种是基于debond的技术（包括VCCT）；一种是基于cohesive技术。 debond即节点松绑，或者称为节点释放，当满足一定得释放条件后（COD 等，目前abaqus提供了5种断裂准则），节点释放即裂纹扩展，采用这种方法时也可以计算出围线积分。 cohesive有人把它译为粘聚区模型，或带屈曲模型，多用于模拟film、裂纹扩展及复合材料层间开裂等。cohesive模型属于损伤力学模型，最先由Barenblatt 引入，使用拉伸-张开法则（traction-separation law）来模拟原子晶格的减聚力。这样就避免了裂纹尖端的奇异性。Cohesive 模型与有限元方法结合首先被用于混凝土计算和模拟，后来也被引入金属及复合材料。Cohesive界面单元要服从cohesive 分离法则，法则范围可包括粘塑性、粘弹性、破裂、纤维断裂、动力学失效及循环载荷失效等行为。 此外，从abaqus6.9版本开始还引入了扩展有限元法（XFEM），它既可以模拟静态裂纹，计算应力强度因子和J积分等参量，也可以模拟裂纹的开裂过程。被誉为最具有前途的裂纹数值模拟方法。本文将利用abaqus6.9版本中的扩展有限元法功能模拟常见的Ⅰ型裂纹的扩展。 2 Ⅰ型裂纹的扩展有限元分析 本文针对断裂力学中的平面Ⅰ型裂纹扩展问题用abaqus中的扩展有限元方法进行数值模拟，获得了裂纹扩展的整个过程，裂尖单元的应力变化曲线，以及裂纹尖端塑性区的形状。在此基础上绘制裂纹扩展的能量历史曲线变化趋势图。

ABAQUS中扩展有限元(XFEM)功能简介

ABAQUS中扩展有限元（XFEM）功能简介 扩展有限元（Extended Finite Element Method）是一种解决断裂力学问题的新的有限元方法，其理论最早于1999年，由美国西北大学的教授Belyschko和Black首次提出，主要是采用独立于网格剖分的思想解决有限元中的裂纹扩展问题，在保留传统有限元所有优点的同时，并不需要对结构内部存在的裂纹等缺陷进行网格划分。 ABAQUS基于在非线性方面的突出优势，在其6.9的版本中开始加入了扩展有限元功能，到6.13做了一些修正，加入了一些可以被CAE支持的关键字。目前为止，除了手动编程，能够实现扩展有限元常用的商业软件只有ABAQUS，今天，我们就来谈谈ABAQUS 中如何实现扩展有限元。 1. XFEM理论 在XFEM理论出现之前，所有对裂纹的静态模拟（断裂）都基本上是采用预留裂缝缺角，通过细化网格仿真裂缝的轮廓。而动态的模拟（损伤）基本上都是基于统计原理的Paris 方法。然而，断裂和损伤的结合问题却一直没有得到有效的解决，究其原因，在于断裂力学认可裂纹尖端的应力奇异现象（就是在靠近裂尖的区域应力值会变无穷大），并且尽可能的绕开这个区域。而损伤力学又没有办法回避这个问题（裂纹都是从尖端开裂的）。 从理论上讲，其实单元内部的位移函数（形函数）可以是任意形状的，但大多数的计算软件都采用了多项式或者插值多项式作为手段来描述单元内部的位移场，这是因为采用这种方法更加便于在编程中进行处理。但是这种方法的缺点就是，由于形函数的连续性，导致单元内部不可能存在间断。直到Belytschko提出采用水平集函数作为手段，其基本形式为 和 上面左边的等式描述了单元内裂缝的位置，右边的等式描述了裂尖的位置。与之对应的形函数便是

基于ABAQUS的渐开线齿轮齿根裂纹扩展仿真

基于AＢＡＱUＳ的渐开线齿轮齿根裂纹扩展仿真

————————————————————————————————作者: ————————————————————————————————日期: ?

基于ABＡQUＳ的渐开线齿轮齿根裂纹扩展仿真 齿轮传动是机械传动中最重要、应用最广泛的一种传动。齿轮传动的主要优点有:传动效率高，工作可靠，寿命长，传动比准确，结构紧凑。齿轮传动的失效一般发生在轮齿上,通常有齿面损伤和齿轮折断两种形式。齿轮折断一般发生在齿根部位,包括疲劳折断和过载折断。 为了提高齿轮的可靠性和使用寿命,有必要对齿轮根部的断裂现象进行研究。本文将从断裂力学角度出发,采用有限元的计算方法,研究齿根的断裂。 1 轮齿断裂分析 应力强度因子是描述裂纹尖端的一个参数,它与载荷大小以及几何有关,共有３种断裂模型（图１),在任何应力下的裂尖应力场为 ?图1 断裂模型 式中:r为距裂尖的距离；θ=arctan(x2/x1);KＩ为Ⅰ型(张开)裂纹应力强度因子；KⅡ为Ⅱ型（张开）应力强度因子。ＫⅢ为Ⅲ型(撕开)应力强度因子。 对于二维裂纹,假定KⅡ为0。

裂纹扩展方向根据条件аσθθ/аθ＝0或者γγθ=0，得到 为了计算二维情况下的积分，ABAQUＳ定义了围线围绕着裂尖由单元组成的环形域(图２)。 图2 裂纹尖端环形域 计算J积分时，围线外的节点处值为0，围线内的所有节点(裂纹 扩展方向)的值为ｌ,但外层单元的中间点除外,这些节点根据在单元中的位置被置于0和１之间。 裂纹扩展角度口可以参考裂纹平面计算，当裂纹扩展方向沿着初始裂纹方向时，θ＝0;当K1>0时,θ<0;当K1<0时,θ>0。裂纹扩展角度从q到n(图3)。

ABAQUS中Cohesive单元建模方法讲解

复合材料模型建模与分析 1. Cohesive单元建模方法 1.1 几何模型 使用内聚力模型（cohesive zone）模拟裂纹的产生和扩展，需要在预计产生裂纹的区域加入cohesive层。建立cohesive层的方法主要有： 方法一、建立完整的结构（如图1（a）所示），然后在上面切割出一个薄层来模拟cohesive 单元，用这种方法建立的cohesive单元与其他单元公用节点，并以此传递力和位移。 方法二、分别建立cohesive层和其他结构部件的实体模型，通过“tie”绑定约束，使得cohesive单元两侧的单元位移和应力协调，如图1（b）所示。 （a）cohesive单元与其他单元公用节点（b）独立的网格通过“tie”绑定 图1.建模方法 上述两种方法都可以用来模拟复合材料的分层失效，第一种方法划分网格比较复杂；第二种方法赋材料属性简单，划分网格也方便，但是装配及“tie”很繁琐；因此在实际建模中我们应根据实际结构选取较简单的方法。 1.2 材料属性 应用cohesive单元模拟复合材料失效，包括两种模型：一种是基于traction-separation 描述；另一种是基于连续体描述。其中基于traction-separation描述的方法应用更加广泛。 而在基于traction-separation描述的方法中，最常用的本构模型为图2所示的双线性本构模型。它给出了材料达到强度极限前的线弹性段和材料达到强度极限后的刚度线性降低软化阶段。注意图中纵坐标为应力，而横坐标为位移，因此线弹性段的斜率代表的实际是cohesive单元的刚度。曲线下的面积即为材料断裂时的能量释放率。因此在定义cohesive的力学性能时，实际就是要确定上述本构模型的具体形状：包括刚度、极限强度、以及临界断裂能量释放率，或者最终失效时单元的位移。常用的定义方法是给定上述参数中的前三项，也就确定了cohesive的本构模型。Cohesive单元可理解为一种准二维单元，可以将它看作被一个厚度隔开的两个面，这两个面分别和其他实体单元连接。Cohesive单元只考虑面外的力，包括法向的正应力以及XZ，YZ两个方向的剪应力。 下文对cohesive单元的参数进行阐述，并介绍参数的选择方法。

ABAQUS精选本FEM扩展元例子的详细图解

版本X F E M(扩展有限元)例子的详细图解 一、part模块中的操作： 二、 1.生成一个新的part，取名为plate，本part选取3Ddeformablesolidextrusion类型（如图1） 三、 2.通过Rectangle工具画出一长3，高6的矩形。考虑使用工具栏add-dimension和 editdimension来画出精确长度的模型。强烈建议此矩形的左上角坐标为（0，3），右下角坐标为（3，-3）（如图2） 四、 3.完成后拉伸此矩形，深度为1.（如图3） 五、图1，图2，图3， 4.生成一个新的part，取名为crack，本part选取3Ddeformableshellextrusion类型（如图4）

5.生成一条线，此线的左端点坐标为（0，），右端点坐标为（，） 6.完成后拉伸此线，深度为1.（如图6） 7.保存此模型为XFEMtutor（如图7），以后经常保存模型，不再累述。 8.在partPlate中分别创建4个集合，分别为：all，bottom，top和fixZ，各部分的内容如图8～11所示。

二、Material模块中的操作： 1.创建材料elsa，其弹性参数为E=210GPa，泊松比为（如图12） 最大主应力失效准则作为损伤起始的判据，最大主应力为（如图13）

损伤演化选取基于能量的、线性软化的、混合模式的指数损伤演化规律，有关参数为 G1C=G2C=G3C=42200N/m，=1.（如图14） 2.创建一个SolidHomogeneous的section，名为solid（如图15），此section与材料elsa相联（如图16），并将此section赋给platepart（也就是集合all）（如图17） 3.赋予材料取向，分别如图18～21所示。

abaqus裂纹模拟心得

abaqus裂纹模拟心得 baqus裂纹模拟心得（Contour Integral不是XFEM） 最近由于项目需要，做了一些裂纹相关的模拟，在此把一些心得体会贴到论坛上与大家分享，如有不当之处，欢迎大家指正！ 本帖主要侧重于介绍裂纹定义过程中各个选项的意义，具体的操作过程论坛里已经有高手做了很好的教程，至于断裂力学理论推荐大家看一下沈成康写的《断裂力学》一书。裂纹的定义和输出需要用到interaction模块和step模块： 一、Interaction模块 1.1 预制裂纹（步骤：菜单/special/crack/assign seam） 注意：并不是作裂纹分析都要定义seam，如果你的裂纹不是一条缝，而是一个缺口，则不需要assign seam，直接走下一步（定义裂纹）就行。 1.2 创建裂纹（步骤：菜单/special/crack/create，type：contour integral） —crack front：crack front是用来定义第一围线积分的区域，2D下我们可以选择包围裂尖点的面，3D则选择包围裂尖线的面；另外还有一种定义crack front的方法，就是直接选择裂尖点（2D）或裂尖线3D），用这个方法定义crack front不需要再定义下一步的crack tip/line，比较简便，两种方法算出的结果没有明显的差别，其实只是影响积分路线的问题，但是J积分值是路径无关的，看个人喜好吧 —crack tip/line：这个比较好理解就是裂尖点（2D）或线（3D），如果我们在上一步中用方法二定义crack front，这一步就直接跳过了 —crack extension direction（定义裂纹扩展方向）：这里定义的其实是一个虚拟的裂纹扩展方向，定义了这个参考方向后，我们才能通过输出的角度判断裂纹扩展方向，可以通过两种方法： o q vector：输入一个方向，用来作为计算裂纹的扩展方向的参考方向； o normal to crack plane：crack plane表示裂纹的对称面（当裂纹在一个平面内时，可能需要分开定义多个裂纹），这种方法下我们只需定义裂纹面的法线方向，通过（t表示裂纹尖端的切线）, 会在每个节点得出一个q方向（如下图）； o 注意：q的方向对输出的应力强度因子，J积分等都会有影响，一般情况下，q最好在裂纹平面内，且垂直于裂尖线的切线，否则算出的应力强度因子，J积分值等等在不同围线积分中会差别较大。 二、step模块 定义好了裂纹相关参数后，我们需要返回step模块定义输出变量： 步骤：菜单/output/history output requests/create，domain：crack，可以输出的值包括：J-integral，Ct-integral，stress intensity factor，T-stress —J-integral ：用于应变率无关材料的准静态分析过程，包括线弹性，非线性弹性，弹塑性材料（单调加载工况）的静态分析。J-integral的优点是和积分路径无关，从而可以避开尖端塑性区的

abaqus6.8和6.9的新功能

ABAQUS6.9新功能 一、扩展有限元（XFEM）的引入 我认为ABAQUS将扩展有限元引入是其最大亮点，也非常有市场前景。该方法可以认为是有限元方法处理不连续问题的革命性变革。据我所知，这是第一个将XFEM商用化的软件。我厂在生产过程中可能遇到的裂纹问题，夹杂问题，气孔，复合材料问题（复合材料的纤维相当于夹杂）都可以通过这种方法解决。虽然我对该方法的理论不胜了解，但是在此我想就我所知对该方法多做些介绍。+ u! _2 v# n& O& M* P 固体力学中存在两类典型的不连续问题，一类是因材料特性突变引起的弱不连续问题，这类问题以双材料问题和夹杂问题为代表，其复杂性由物理界面处的应变不连续性引起;另一类是因物体内部几何突变引起的强不连续问题，这类问题以裂纹问题为代表，其复杂性由几何界面处的位移不连续性和端部的奇异性引起.物体内部物理界面的脱粘或起裂，是上述两类问题的混合.另外，在复杂流体、复杂传热、物质微结构演化等复杂问题中，也存在许多不连续力学问题。* e4 o5 F2 K4 n& g n3 {4 U 数值方法，如有限元、边界元、无单元法等，一直是处理不连续问题的主要途径.有限元法具有其它数值方法无可比拟的优点，即适用于任意几何形状和边界条件、材料和几何非线性问题、各向异性问题、容易编程等，因而成为数值分析裂纹等不连续问题的主要手段。更为重要的是有限元方法的商业化程度和推广程度都很高，对于我们这样的企业已经有一大批熟悉有限元技术的工程师。因此通过有限元方法来解决这一问题是成本低效率高的途径。但是传统的有限元方法在处理裂纹，夹杂，空隙这些不连续问题时困难非常大。) t, Z- t8 ?2 J+ ~3 f: Q 常规有限元法(CFEM)采用连续函数作为形状(插值)函数，要求在单元内部形状函数连续且材料性能不能跳跃，在处理像裂纹这样的强不连续(位移不连续)问题时，必须将裂纹面设置为单元的边、裂尖设置为单元的结点、在裂尖附近的高应力区需要令人难以接受的网格密度，同时在模拟裂纹生长时还需要对网格进行重新剖分。现在绝大多数商业软件在模拟裂纹扩展问题时都需要预设裂纹的扩展方向，而且在裂纹扩展过程中不断的从新划分网格，效率极低甚至无能为力。在处理多裂纹问题时，其求解规模之大、网格剖分之难是不可想象的。处理夹杂问题时，要求单元的边必须位于夹杂与基体的界面处，即使对于网格自动化程度很高的二维问题这也不容易，更何况拓扑结构更复杂的三维问题。 1999年以来，在有限元框架内发展起来的扩展有限元法，以解决不连续问题为着眼点，对常规有限元法在求解裂纹问题时所遇到的困难提出了近乎完美的解决方案。' G, ]4 z! c/ i$ d2 T 1999年，以美国西北大学Belytschko教授为代表的研究组首先提出扩展有限元思想，20 00年，他们正式使用扩展有限元法(XFEM)这一术语即。XFEM是迄今为止求解不连续力学问题最有效的数值方法，它在标准有限元框架内研究问题，保留了CFEM的所有优点，但并不需要对结构内存在的几何或物理界面进行网格剖分。XFEM与CFEM的最根本区别在于所使用的网格与结构内部的几何或物理界面无关，从而克服了在诸如裂纹尖端等高应力和变形集中区进行高密度网格剖分所带来的困难，当模拟裂纹扩展时也无需对网格进行重新剖分。 也就是说在裂纹的扩展过程中裂纹可以穿透单元扩展。就其理论我们可以简单的理解为在单元内部有很多的潜在节点，当需要时这些节点被激活实现裂纹穿透单元扩展。在宣讲会上ABAQUS业务代表现场展示了其用XFEM完成的II裂纹扩展模拟，裂纹面沿70°方向穿透单元扩展 我有个师兄也在从事XFEM方面的研究工作。他在研究中遇到的主要问题是材料中有大小随机的夹杂的时候可能会出现一个单元将一个或多个夹杂完全包围在其中的情况。处理这种

abaqus有限元分析报告开裂梁要点

Abaqus梁的开裂模拟计算报告 1.问题描述 利用ABAQUS有限元软件分析如图1.1所示的钢筋混凝土梁的裂缝开展。参考文献Brena et al.（2003）得到梁的基本数据： 图1.1 Brena et al.（2003）中梁C尺寸 几何尺寸：跨度3000mm，截面宽203mm，高406mm的钢筋混凝土梁 由文献Chen et al. 2011得材料特性： 1.混凝土：抗压强度f c’=35.1MPa，抗拉强度f t= 2.721MPa，泊松比ν=0.2，弹性模量 E c=28020MPa； 2.钢筋：弹性模量为E c=200GPa，屈服强度f ys=f yc=440MPa，f yv=596MPa 3.混凝土垫块：弹性模量为E c=28020MPa，泊松比ν=0.2 2.建模过程 1)Part 打开ABAQUS使用功能模块，弹出窗口Create Part，参数为：Name：beam；Modeling Space：2D；Type：Deformable；Base Feature─Shell；Approximate size：2000。点击Continue 进入Sketch二维绘图区。由于该梁关于Y轴对称，建模的时候取沿X轴的一半作为模拟对象。 使用功能模块，分别键入独立点(0,0),(1600,0),(1600,406),(406,0),(0,0)并按下下方提 示区的Done，完成草图。 图2.1 beam 部件二维几何模型

相同的方法建立混凝土垫块： 图2.2 plate 部件二维几何模型 所选用的点有(0,0),(40,0),(40,10),(0,10) 受压区钢筋： 在选择钢筋的base feature的时候选择wire，即线模型。 图2.3 compression bar 部件二维几何模型 选取的点(0,0),(1575,0) 受拉区钢筋： 图2.4 tension bar 部件二维几何模型 选取的点(0,0),(1575,0) 箍筋： 图2.5 stirrup 部件二维几何模型 选取的点为(0,0),(0,330) 另外，此文里面为了作对比，部分的模型输入尺寸的时候为m，下面无特别说明尺寸都为mm。

第5章+在ABAQUS中构造有限元模型的若干问题-

第五章 构造有限元模型的若干问题 第三章所述的岩土工程中常用的本构模型可用于以下四种类型的单元中： 1. 平面应变单元 2. 广义平面应变单元 3. 轴对称单元 4. 三维单元 在应用平面应变单元时，只能用子午线作为线性的Drucker-Prager 模型。 考虑非线性变形的轴对称单元具有如下功能：初始变形是轴对称的，在进入非线性变形后，允许产生非轴对称变形，这个功能是十分有用的，比如桩基承受轴向力、水平力和弯矩时，它在r-z 平面上采用标准的插值函数，在与θ有关的方向上采用Fourier 级数插值函数，这类单元可以充分考虑不同方位上变形的差异，其计算结果更接近实际。 考虑流体在多孔介质中的作用时，需采用位移——孔隙压力耦合单元，耦合单元可应用于平面应变单元、轴对称/反对称单元、三维单元。此时，位移与孔隙压力的插值函数可以不同，孔隙压力通常采用线性插值已经足够，位移则可选择线性插值或二次插值。 有限元与无限元相结合来求解考虑这个问题是一个有效的方法。 5-1 广义平面单元 所谓广义平面单元是指位于两个受约束的平面之间的区域，这两个受约束的平面可以如同刚体一般绕轴相互转动，这个转动会导致沿厚度方向上的应变。假定这个应变在厚度方向上与位置无关，则两个平面的相对移动仅仅在厚度方向上引起应变。这样引起的应变以及一阶与二阶变量由下式定义： 设P 0（x 0，y 0）为受约束平面上的一个固定点，P 0到另一个受约束平面的距离为z u t ?+0，其中t 0是初始距离，z u ?是变形过程中产生的位移，z u ?为广义平面元节点自由度的一阶变量。对于用广义平面元进行连结的某一区域而言，只有一个广义平面元节点自由度。不同的连结区域可以有不同的广义平面应变节点，因此两个受约束的平面之间的这个区域可视为一个刚体，该区域中任何其它点（x ，y ）的纤维长度为： y x z Z x Y y u t t φ???φ??+?+=)()(000 （5-1-1） 其中 10)(|)(x x x φ?+φ?=φ? （5-1-2） 10)(|)(y y y φ?+φ?=φ? （5-1-3）

Abaqus裂纹模拟心得(Contour Integral不是XFEM)

Abaqus裂纹模拟心得（Contour Integral不是XFEM） 最近由于项目需要，做了一些裂纹相关的模拟，在此把一些心得体会贴到论坛上与大家分享，如有不当之处，欢迎大家指正！ 本帖主要侧重于介绍裂纹定义过程中各个选项的意义，具体的操作过程论坛里已经有高手做了很好的教程，至于断裂力学理论推荐大家看一下沈成康写的《断裂力学》一书。 裂纹的定义和输出需要用到interaction模块和step模块： 一、Interaction模块 1.1 预制裂纹（步骤：菜单/special/crack/assign seam） 注意：并不是作裂纹分析都要定义seam，如果你的裂纹不是一条缝，而是一个缺口，则不需要assign seam，直接走下一步（定义裂纹）就行。 1.2 创建裂纹（步骤：菜单/special/crack/create，type：contour integral） —crack front：crack front是用来定义第一围线积分的区域，2D下我们可以选择包围裂尖点的面，3D则选择包围裂尖线的面；另外还有一种定义crack front的方法，就是直接选择裂尖点（2D）或裂尖线3D），用这个方法定义crack front不需要再定义下一步的crack tip/line，比较简便，两种方法算出的结果没有明显的差别，其实只是影响积分路线的问题，但是J 积分值是路径无关的，看个人喜好吧 —crack tip/line：这个比较好理解就是裂尖点（2D）或线（3D），如果我们在上一步中用方法二定义crack front，这一步就直接跳过了 —crack extension direction（定义裂纹扩展方向）：这里定义的其实是一个虚拟的裂纹扩展方向，定义了这个参考方向后，我们才能通过输出的角度判断裂纹扩展方向，可以通过两种方法： o q vector：输入一个方向，用来作为计算裂纹的扩展方向的参考方向； o normal to crack plane：crack plane表示裂纹的对称面（当裂纹在一个平面内时，可能需要分开定义多个裂纹），这种方法下我们只需定义裂纹面的法线方向，通过（t表示裂纹尖端的切线）, 会在每个节点得出一个q方向（如下图）； o 注意：q的方向对输出的应力强度因子，J积分等都会有影响，一般情况下，q最好在裂纹平面内，且垂直于裂尖线的切线，否则算出的应力强度因子，J积分值等等在不同围线积分中会差别较大。

ABAQUS(Explict)混凝土开裂模型翻译

混凝土开裂模型 适用模块：Abaqus/Explicit Abaqus/CAE 参考 ●“Material library: overview,” Section 18.1.1 ●“Inelastic behavior,” Section 20.1.1 ●*BRITTLE CRACKING ●*BRITTLE FAILURE ●*BRITTLE SHEAR ●“Defining brittle cracking” in “Defining other mechanical models,” Section 12.9.4 of the Abaqus/CAE User's Manual 概述 Abaqus/Explicit模块中脆性断裂模型： ●提供一种通用模型来模拟包括梁单元、桁架单元、壳单元以及实体单元在内的所有单元 形式； ●也可以用来模拟诸如陶瓷及脆性岩石的其他材料； ●用于模拟受拉开裂占主导地位的材料本构行为； ●假设受压行为是线弹性的； ●必须与线弹性模型(“Linear elastic behavior,” Section 19.2.1)同时使用，它也 定义了材料开裂前的本构行为； ●用于模拟脆性行为占主导地位的本构关系是十分准确的，基于此，假设受压行为是线弹 性的是合理的； ●该模型主要是用于钢筋混凝土结构的分析，同时也适用于素混凝土； ●基于脆性失效准则，将失效单元删除； 关于失效单元删除的内容详见“A cracking model for concrete and other brittle materials,” Section 4.5.3 of the Abaqus Theory Manual. 关于ABAQUS中可用混凝土本构模型的相关讨论参见“Inelastic behavior,” Section 20.1.1。 钢筋 ABAQUS中，混凝土结构中的钢筋是通过指定Rebar单元实现的。Rebar单元是一维应变理论单元（杆单元），既可以单独定义，也可以镶嵌在有向曲面上。关于Rebar单元的相关讨论参见“Defining rebar as an element property,” Section 2.2.4。Rebar单元特别地用来模拟钢筋的弹塑性行为，并且可叠加在用于模拟素混凝土标准单元的网格上。基于这种模拟方法，混凝土的开裂行为与Rebar是没有关系的。混凝土和钢筋之间的相互作用，例如粘结滑移、销栓作用，可以引入拉神硬化（强化）的概念来近似模拟混凝土裂缝处荷载向钢筋转移的特点。 开裂 Abaqus/Explicit中使用弥散裂缝模型来表征混凝土脆性行为的非连续性。这种方法并不关注于单个宏观裂缝，相反地，只是独立地计算有限元模型质点处的本构关系。裂缝的存在对于计算的影响在于：裂缝的存在影响质点处的应力和材料刚度。 为简化本部分讨论内容，“开裂”一词实质上代表的是一个方向——所考虑单个质点处裂缝的方向，与其最相近的物理解释为：在一个质点附近出现一系列连续的微裂缝，其方向

abaqus扩展有限元例子

part模块中的操作： 1. 生成一个新的part，取名为plate，本part选取3D deformable solid extrusion类型（如图1） 2.通过Rectangle工具画出一长3，高6的矩形。考虑使用工具栏add-dimension和edit dimension来画出精确长度的模型。强烈建议此矩形的左上角坐标为（0，3），右下角坐标为（3，-3）（如图2） 3. 完成后拉伸此矩形，深度为1.（如图3） 4. 生成一个新的part，取名为crack，本part选取3D deformable shell extrusion类型（如图4） 2009-7-3 09:22:22 上传下载附件(30.84 KB) 5. 生成一条线，此线的左端点坐标为（0，0.08），右端点坐标为（1.5，0.08）

2009-7-3 09:22:23 上传下载附件(27.27 KB) 6 . 完成后拉伸此线，深度为1.（如图6） 7.保存此模型为XFEMtutor（如图7），以后经常保存模型，不再累述。 8. 在part Plate中分别创建4个集合，分别为：all，bottom，top和fixZ，各部分的内容如图8～11所示 Material模块中的操作： 1.创建材料elsa，其弹性参数为E=210GPa，泊松比为0.3（如图12）

2. 3.最大主应力失效准则作为损伤起始的判据，最大主应力为8 4.4MPa（如图13） 4. 损伤演化选取基于能量的、线性软化的、混合模式的指数损伤演化规律，有关参数为G1C= G2C= G3C=42200N/m，a=1.（如图14）

Abaqus裂纹扩展分析

定义初始粘合裂纹面 可能裂纹表面建模时采用采用主、从接触面来定义。在接触形式中，除了有限滑动、面对面形式以外，其他所有接触形式均可使用。预先定义的裂纹面在初始时应部分粘合，裂纹尖端因而可以被Abaqus/Standard显式识别。初始粘合裂纹面不能采用自接触形式。 定义初始状态（initial condition）以识别裂纹初始绑定部分。用户可以定义从接触面（slave surface）、主接触面（master surface）、以及用来识别从接触面初始部分粘结的节点。从接触面上没有粘结的部分表现为正常接触面。主接触面及从接触面均需要指明。 如果没有节点如上所述被定义，初始接触状态将被应用于整个接触对。这种情况下，不能识别出裂纹尖端，因而粘结面不能分开。 如果节点如上所述被定义，初始解除状态将被应用于从接触面上已定义的节点处。Abaqus/Standard将进行核对以确保所定义节点只包含从接触面上的节点。 *INITIAL CONDITIONS, TYPE=CONTACT 激活裂纹扩展能力（crack propagation capacibility） 裂纹扩展能力需要在STEP定义中被激活，以确保初始部分粘合的2个面有可能产生裂纹扩展。用户需要指明会产生裂纹扩展的面。 *DEBOND, SLA VE=slave_surface_name，MASTER=master_surface_name 多裂纹扩展 裂纹可以在一个或多个裂纹尖端处产生扩展。一个接触对可以在多个裂纹尖端处产生裂纹扩展。然而，对于给定的接触对只能拥有一个裂纹扩展准则（crack propagation criterion）。 定义开裂振幅曲线（debonding amplitude curve） 开裂产生后，通过从接触面节点及主接触面相应节点上大小相等方向相反的力产生面间牵引。当采用临界应力准则、临界裂纹开口位移准则、裂纹长度-时间破坏准则时，用户可以定义粘结面上某点产生开始时，上述力以何种方式降至零。 最佳振幅曲线的选择依赖于材料性质、具体加载情况、裂纹扩展准则。如果应力过快被消除，在裂纹尖端处应变的剧烈变化将会造成收敛困难。对于大变形问题，将会产生严重的网格畸形。对于采用率无关材料的问题，线性振幅曲线完全能够满足要求。对于采用率相关材料的问题，产生开裂时的应力应以较慢速度消耗，以避免收敛困难和网格畸形。在50%

Abaqus针对复合材料优势

四Abaqus在复合材料领域的优势 4.1 复合材料介绍 4.1.1 复合材料的应用 复合材料有许多特性： 1、制造工艺简单 2、比强度高，比刚度大 3、具有灵活的可设计性 4、耐腐蚀，对疲劳不敏感 5、热稳定性能、高温性能好 由于复合材料的上述优点，在航空航天、汽车、船舶等领域，都有广泛的应用。复合材料的大量应用对分析技术提出新的挑战。

4.1.2 复合材料的结构 复合材料是一种至少由两种材料混合而成的宏观材料，其中的一种材料被称作基体，其它的材料称作纤维。其中纤维可以包含很多不同的 形式：离散的宏观粒子，任意方向的短纤维，规则排列的纤维和织物。 1)单向纤维层合板----冲击分析

2)编织复合材料---- 挤压分析 3)蜂窝夹心复合材料----不可见冲击损伤分析

基体和纤维的存在形式以及材料属性对于复合材料的力学行为有着很大的影响。改变纤维和基体的属性目的就是在于生成一种复合材料具有如下性质： 1）低成本：原型，大规模生产，零件合并，维修，技术成熟。 2）期望的重量：轻重量，比重分配合理。 3）改进的强度和刚度：高强度/高刚度比。 4）改进的表面属性：良好的耐腐蚀性，表面抛光性好。 5）期望的热属性：较低的热传导性，热膨胀系数较低。 6）独特的电属性：具有较高的绝缘强度，无磁性。 7）空间适应性：大部件，特殊的几何构型。 4.1.4 复合材料的有限元模拟 根据不同的分析目的，可以采用不同的复合材料模拟技术： 1）微观模拟：将纤维和基体都分别模拟为可变形连续体。 2）宏观模拟：将复合材料模拟为一个正交各向异性体或是完全各向

异性体。 3）混合模拟：将复合材料模拟为一系列离散、可见的纤维层合板。4）离散纤维模拟：采用离散单元或是其它模拟工具进行模拟。5）子模型模拟：对于研究加强纤维周围点的应力集中问题比较有效。 微观模拟：纤维-基体的单胞模拟 混合模拟：层合板的混合模拟 Abaqus中复合材料的单元技术 Abaqus中复合材料的单元技术主要为三种：分层壳单元、分层实体单元以及实体壳单元。 分层壳单元 单元类型：S4, S3R；