下载文档原格式
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所示。
三、划分网格：网格控制为：Hex型structured（如图22），单元类型为C3D8R（如图23）设置plate各边的网格种子为8，26，36（如图24），各边种子的个数不能改变（如图25）四、装配模块：选中plate和crack两个part，分别生成2个实体（如图26），生成一个参考点，参考点的坐标为（1.5，-3，0）（如图27，28）。
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.固定裂纹中，只有各向同性材料的裂纹尖端渐进场才被考虑。
ABAQUS中扩展有限元（XFEM）功能简介扩展有限元（Extended Finite Element Method）是一种解决断裂力学问题的新的有限元方法，其理论最早于1999年，由美国西北大学的教授Belyschko和Black首次提出，主要是采用独立于网格剖分的思想解决有限元中的裂纹扩展问题，在保留传统有限元所有优点的同时，并不需要对结构内部存在的裂纹等缺陷进行网格划分。
ABAQUS基于在非线性方面的突出优势，在其6.9的版本中开始加入了扩展有限元功能，到6.13做了一些修正，加入了一些可以被CAE支持的关键字。
目前为止，除了手动编程，能够实现扩展有限元常用的商业软件只有ABAQUS，今天，我们就来谈谈ABAQUS 中如何实现扩展有限元。
1. XFEM理论在XFEM理论出现之前，所有对裂纹的静态模拟（断裂）都基本上是采用预留裂缝缺角，通过细化网格仿真裂缝的轮廓。
而动态的模拟（损伤）基本上都是基于统计原理的Paris 方法。
然而，断裂和损伤的结合问题却一直没有得到有效的解决，究其原因，在于断裂力学认可裂纹尖端的应力奇异现象（就是在靠近裂尖的区域应力值会变无穷大），并且尽可能的绕开这个区域。
而损伤力学又没有办法回避这个问题（裂纹都是从尖端开裂的）。
从理论上讲，其实单元内部的位移函数（形函数）可以是任意形状的，但大多数的计算软件都采用了多项式或者插值多项式作为手段来描述单元内部的位移场，这是因为采用这种方法更加便于在编程中进行处理。
但是这种方法的缺点就是，由于形函数的连续性，导致单元内部不可能存在间断。
直到Belytschko提出采用水平集函数作为手段，其基本形式为和上面左边的等式描述了单元内裂缝的位置，右边的等式描述了裂尖的位置。
与之对应的形函数便是和其中H(x)是阶跃函数。
想要了解更深的内容，大家可以参考《Extended Finite Element Method》和庄老师的《扩展有限单元法》这种扩充形函数能够描述单元内位移场在裂缝两边的跳跃性，同时，由于裂缝存在于单元内部，其扩展独立与其他单元，使得计算变得高效。
扩展有限元的ABAQUS实现常规有限元方法（CEFM）和其他数值方法相比，具有一些无法比拟的优点，但仍存在一些缺陷。
比如在解决类似裂纹这样的强不连续问题，由于裂纹尖端处的应力奇异性，导致计算量巨大而且精度不高。
然而扩展有限元方法（extended finite element method，XFEM）的出现，和常规有限元方法相比具有显著的优势，使得我们可以在裂尖和应力、变形集中处划分高密度的网格，也可以方便的模拟裂纹的扩展，使计算量不那么巨大，保留了常规有限元法的所有优点。
因此，扩展有限元得到了快速发展和应用，而且在裂纹的扩展研究中要的意义。
本文开展对扩展有限元方法和裂纹问题的研究，并且基于限元ABAQUS平台，对扩展有限元方法针对裂纹扩展问题进行模拟实现。
关键词：扩展有限元方法，裂纹扩展，ABAQUS第一章绪论1.1 引言21世纪以来，计算机硬件和数值仿真的快速发展以及工业工程实践与科学研究中存在的大量运算需求，世界上涌现出一批大型科研运算及科学模拟软件，能够极大的简化运算问题以及计算机模拟实验，使我们能够更加方便地研究虚拟工程及相关科学问题。
有限元方法的出现为数值分析方法的研究带来了新的曙光，力学学科本来就是连接理工学科的桥梁，计算力学是目前力学发展的一个重要分支。
有限元软件则是我们到达工程科学领域彼岸的非常重要的工具和桥梁之一。
ABAQUS软件是世界上最强大的大型有限元计算分析软件之一，具有不同种类的单元类型、材料类型和不同的分析过程，拥有很好的计算功能和模拟性能。
ABAQUS软件不但可以进行一种部件和复杂物理场的分析，而且可以处理多系统的部件分析；不仅可以分析简单的线弹性问题，还可以处理复杂的非线性组合问题等，相比其它软件具有无可比拟的优势[1,2]。
固体力学中存在的两类不连续问题之一则是因为物体内部几何结构突变引起的强不连续问题，裂纹问题就是这类问题的代表。
由于几何界面处的位移不连续性和裂纹尖端的应力奇异性使得这类问题的处理变得比较复杂。
基于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积分等参量，也可以模拟裂纹的开裂过程。
版本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所示。
三、划分网格：网格控制为：Hex型structured（如图22），单元类型为C3D8R（如图23）设置plate各边的网格种子为8，26，36（如图24），各边种子的个数不能改变（如图25）四、装配模块：选中plate和crack两个part，分别生成2个实体（如图26），生成一个参考点，参考点的坐标为（，-3，0）（如图27，28）。
基于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积分等参量，也可以模拟裂纹的开裂过程。
在Abaqus中进行疲劳裂纹扩展模拟通常需要使用ABAQUS/Standard或ABAQUS/Explicit这两个分析模块。
ABAQUS提供了丰富的工具和元素来模拟疲劳裂纹扩展，以下是一个基本的步骤：1. 建模：-使用ABAQUS/CAE（图形用户界面）或ABAQUS脚本语言（Python）创建模型。
确保模型包含准确的几何形状和边界条件。
2. 网格划分：-确保模型的网格划分足够细致，特别是在裂纹尖端区域。
使用ABAQUS 提供的适当类型的网格元素，如二维或三维等元素。
3. 材料定义：-定义材料的力学性质和断裂参数。
在疲劳分析中，通常需要使用合适的疲劳材料参数。
4. 加载和约束：-定义加载和约束条件。
对于疲劳裂纹扩展，通常使用周期性的加载。
加载可以是压力、力、位移等。
5. 疲劳裂纹增长：-使用ABAQUS的断裂力学（XFEM）方法来模拟裂纹的扩展。
你可以使用ABAQUS/Standard的XFEM方法来处理裂纹尖端的应力集中。
6. 结果输出：-设置合适的输出请求以获得关于裂纹扩展和结构响应的信息。
这可能包括应力、应变、位移、裂纹长度等。
7. 迭代分析：-如果需要模拟多个加载循环的疲劳裂纹扩展，你可能需要使用ABAQUS/Standard的循环加载功能，或者通过ABAQUS/Explicit进行显式动态疲劳分析。
8. 后处理：-使用ABAQUS/CAE或Python脚本进行后处理，绘制结果图形，分析裂纹扩展速率等。
请注意，这仅仅是一个基本的指南。
实际应用中，还需要考虑更多因素，如裂纹尖端应力场的准确建模、裂纹扩展准则的选择等。
确保在模拟前仔细阅读ABAQUS文档，并根据具体问题和标准进行模拟设置。
基于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积分等参量，也可以模拟裂纹的开裂过程。
abaqus XFEM 方ń樯abaqus XFEM 方法介绍aspan class= 裂纹扩展分析体验热分析中的热物性参数:材料密度，热导率，比热容，电阻率，弹性模量，融化潜热的焓、泊松比、散热系数。
在ABAQUS模型中，需用3D的deformable、shell、exctrusion方式建立一个初始裂纹，长短适宜，初始裂纹要从开始起裂的点设置。
由于计算方法目前还不稳定，参数要适当调整。
设置网格划分参数的时候，应该对称设置，否则网格不对称。
断裂应力的大小要和断裂能量的设置相应，能量太大、太小导致不易收敛，断裂区域的网格要规则，各个方向尺寸要差不多，整个厚度方向单元数量一致，且越少越好，即使裂纹起始点两侧单元未参与裂纹，也要尽可能均匀规则，裂纹扩展的区域不能被PARTION开，应该是一体的。
冲击动载荷时，载荷步时间应尽可能小，maxps damage应力应大于ductile damage应力。
初始裂纹不能在单元界限扩展，否则导致不收敛。
裂纹可在两种弹性金属界面上。
适用于弹性材料、两种弹性材料界面裂纹和幂硬化材料。
我建立的弹塑性材料模型不容易收敛，把塑性去掉后反而容易收敛。
建立ductile manage模型时，需要材料的塑性行为，但必须同时有traction manage模型，否则就提示某些单元的fracture damage模型未能建立,。
traction manage模型和ductile manage模型中的damage evolution中的类型和数值要相同。
初始屈服应力和断裂应力不能差太多。
动力学的冲击裂纹分析，不能有塑性(当然也不能有ductile damage)，可能是因为increment time 太大了,(或者是先分析一下不带塑性材料的，通过了才分析有塑性材料的,) 断裂能量要和断裂应力相适应，否则可能因为能量太高，还未达到，但应力达到了，导致裂纹不开裂、计算不收敛。
能量太低，很容易就满足了能量开裂原则，但应力未达到，裂纹不能开裂、计算不收敛。
基于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 积分等参量，也可以模拟裂纹的开裂过程。
Abaqus提取XFEM(扩展有限元)裂缝长度【壹讲壹插件】2015-7-20作者：星辰-北极星Abaqus提取XFEM(扩展有限元)裂缝长度 (1)第一部分：Abaqus 扩展有限元方法XFEM (2)1.1概要 (2)1.2这些你有注意到吗？ (2)1.3 圆孔内压裂缝模拟实例 (2)1.3.1 部件建立 (2)1.3.2 材料性质定义(part1) (2)1.3.3 分析步定义 (3)1.3.4 参数输出 (3)1.3.5 接触模块定义Crack (3)1.3.6边界条件定义 (4)1.3.7 网格划分 (4)1.3.8初始地应力施加 (4)1.3.9 计算结果： (4)第二部分：扩展有限元裂缝长度求解 (5)2.1 概要 (5)2.2 基本求解思路： (5)第三部分：星辰-北极星插件介绍：POLARIS-XFEMCreckGeo2D (6)3.1 概要 (6)3.2 插件的主要功能 (6)3.3 使用注意事项 (6)3.4 插件使用简介 (7)3.4.1 打开插件 (7)3.4.2 数据获取 (7)3.4.3 裂缝信息获取 (8)3.4.4 示例 (8)第一部分：Abaqus 扩展有限元方法XFEM文章转自：/908754116/blog/14374022441.1概要XFEM即扩展有限元方法，它在标准有限元框架内研究问题，保留了有限元方法的所有优点。
扩展有限元法与有限元法最根本的区别在于所使用的网格与结构内部的几何或物理界面无关，从而克服了在诸如裂纹尖端等高应力和变形集中区进行高密度网格划分所带来的困难，在模拟裂纹扩展时也无需对网格进行重新划分。
如果要正常地使用它，我们首先要了解Abaqus中的扩展有限元方法有哪些特别，它在理论上做了哪些简化等，帮助文档进行了很好的讲解：《Abaqus Analysis User's Manual》10.7.1 Modeling discontinuities as an enriched feature using the extended finite element method。
基于扩展有限元（XFEM ）裂纹扩展总结通过四个算例总结了用ABAQUS 计算裂纹扩展应用情况。
算例1基于XFEM 使用虚拟裂缝闭合技术结合Cohesive 单元，实现混凝土基体断裂和钢筋混凝土界面脱层的混合失效模式；算例2基于XFEM 以VCCT 准则判断裂缝的开裂扩展，研究了偏荷载作用下不同配筋率对裂缝扩展方向的影响，并对比了考虑钢筋与混凝土粘结滑移与不考虑粘结滑移的裂缝扩展情况；算例3则是以粘聚力模型判断裂缝扩展,研究了裂缝扩展情况；算例4对比了Cohesive 和VCCT 两种开裂准则下钢筋混凝土（纵、箍筋组合）的裂缝扩展情况。
扩展有限元基本原理扩展有限元法（XFEM ）是在单位分解法的基础上对常规有限元位移逼近函数进行改进加强，引入附加函数。
以二维裂纹（图1）为例，对于裂纹贯穿单元，采用Heaviside 函数来描述裂纹两侧的不连续性；对于裂尖单元，采用裂尖渐进函数来反映裂纹尖端应力的奇异性。
扩展有限元的位移逼近为：()()()()()()∑∑∑∑∈=∈∈⎥⎦⎤⎢⎣⎡++=K i i i i J i i i I i i xfem b x F x N a x H x N u x N x u 41ααα （1）式中，I 为所有节点集合，()x N i 为节点i 的形函数，i u 为节点i 的标准自由度，J 为裂纹贯穿单元节点集合（图1中圆圈所示节点），K 为裂尖单元节点集合（图1中方形所示节点），()x H 和()x F α分别为Heaviside 形函数和裂尖渐进函数，i a 和αi b 为相应节点自由度。
图1 扩展有限元中的富集节点描述裂纹面不连续性的Heaviside 形函数可表示为 ()⎩⎨⎧−≥•=*otherwise 10n )x -(x if 1x H （2）式中，*x 为点x 到裂纹面最近处的投影，n 为*x 点处的单位外法线向量（如图2所示）。
可以看出，节点位于裂纹面上侧时()1=x H ，节点位于裂纹面下侧时()1−=x H ，Heaviside 形函数能较好的描述裂纹面两侧的不连续性。
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）最大主应力失效准则作为损伤起始的判据，最大主应力为84.4MPa（如图13）损伤演化选取基于能量的、线性软化的、混合模式的指数损伤演化规律，有关参数为G1C= G2C= G3C=42200N/m， =1.（如图14）2.创建一个Solid Homogeneous 的section，名为solid（如图15），此section与材料elsa相联（如图16），并将此section赋给plate part（也就是集合all）（如图17）3.赋予材料取向，分别如图18～21所示。
三、划分网格：网格控制为：Hex型structured（如图22），单元类型为C3D8R（如图23）设置plate各边的网格种子为8，26，36（如图24），各边种子的个数不能改变（如图25）四、装配模块：选中plate和crack两个part，分别生成2个实体（如图26），生成一个参考点，参考点的坐标为（1.5，-3，0）（如图27，28）。
ABAQUS XFEM 柱面与平面滑动接触中裂纹的扩展xfem(step by step)
1.part部分：
plate模型为2D 变形壳体尺寸为0.5x0.2M(因为后面采用的是MPa，所以这里采用的是M，为了单位统一) 。
用partition命令将plate切割成装配图上面一样（尺寸0.24x0.06）Rigid模型为2D解析刚体尺寸为R0.06圆的下面部分的2/3半圆（看装配模型就知道）。
注意要在上面设置一个参考点，在load部分加载位移边界用。
wire模型2D变形线尺寸为0.01位置坐标（0.25,0.2），（0.25,0.19）线两端点的坐标，大家懂的。
添加一个装配图对照模型。
2.property部分这一部分尽量多用图片
主要设置了一下几部分的材料属性，用的材料为Q345-steel（因为刚好找到了它的应力应变数据）1、密度2、弹性属性3、塑性属性4、损失准则
5.interaction部分
1.首先创建contact controls，创建步棸菜单栏interaction->contact controls->create,参数的设置见下图：
2.裂纹的设置，创建步棸菜单栏special->crack->create->xfem ,区域的选择见下图
3其他包括相互作用的一些属性见下图
6.load部分
位移边界条件的加载如下图
7.Mesh部分
种子布置见下图
网格控制主要采取的是structured
8.Job就ok了，祝大家好运。
也同时感谢论坛里面帮助过我的朋友。
基于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积分等参量，也可以模拟裂纹的开裂过程。
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）3D '、2D Planw I ' AxisymmetricTyre Options” Di scr^te ri gi >1f'■ Analytical ri<id0 Euler i anBaseShapeSolidC? Shello灿电0 FeintApproKimat色size： 20Cancel3.完成后拉伸此矩形，深度为 1.（如图3）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 dCj An>lytic41 rigid ■.. j Euler i an Opti QKS None available Q hl 迥⑥*1】■：.\ Wire (.PoiittBase FeatureTypePlanar Ez trusi on Rezolution SwsepXppr^MiTatt =it eCine el5.生成一条线，此线的左端点坐标为( 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）最大主应力失效准则作为损伤起始的判据，最大主应力为84.4MPa （如图13）2.创建一个 Solid Homogeneous 的 section ，名为 solid （如图 15），此 section 与材料elsa 相联（如图16），并将此section 赋给plate part（也就是集合 all ）（如图17）T K I KFB ¥dit ■时” iiU■ 1； Tl h M 丿；N 心* D SJLIfTirit CT ci fi vzxi ablas ：损伤演化选取基于能量的、 线性软化的、混合模式的指数损伤演化规律,有关参数为G1 C=G2C= G3C=42200N/m =1.（如图 14）［恥Hiwl□罚也4 ■ - Ja#ih-Ka«i!杠ini-M*l| .111出帝HI錨Z扁玳刃3.赋予材料取向，分别如图18〜21所示。
WE3-c tf込：=4节-T px li』F弓二"・二鑑 日.】M - !-•■»«■■■;吕•g 応一■■■=■」！，■■■二 」ij WM L V -P工 f融.n d・Hr-m口」-g ・r m p二』wwlr-n 存翌 myu:i?E --&M-:L4」$F占報仪腿B旺u（如图22）,单元类型为 C3D8R （如图23）设置plate各边的网格种子为 8，26，36 （如图24），各边种子的个数不能改变（如图三、划分网格:窗 fiH<l -' :T网格控制为：Hex 型structured 25)四、装配模块：选中plate和crack两个part，分别生成2个实体（如图26）,生成一个参考点，参考点的坐标为（1.5，-3，0）（如图27, 28 ）。
将参考点的名字改为db （如图29，30）o生成集合bdisp，此集合包含db这个点。
和:呻祜•单i■刖■Triili Fblll ER-«ik ZitxIrMLZiLa卜'ir Todfrl ■ist.-baar :z to U,Z XETEMt ator .匚' “ rodZd Prt•呼Hr七三：--3c te^Ti s-T-^e 1 to 匚yrEHtTa+cr. c TS-'1F-bs h -bb - E kwa bM>・■<•'¥■4「■厂I v . r —*yr ■抑■“实用标准文档-F eatures (2)Hats c sys~ 1The model database 743& el emeuts have The model d^t^base The jnodel d^tabas^has been saved ta been geneia tud on hes be^n to hn 守 been saved to五、In teraction 模块中创建约束方程 ce_bot ,如图31所示。
、Bw px pp m pDelete.,. Ehh 卯 匚hildrenQuary Sh^w Parents Stej Fi 创Kisi T inn AT.EIdit..・六.Step模块中的操作:建立个static的分析步，如图32所示。
打开几何非线性开关，如图33所示。
设置增量步，如图34所示。
打开解控制，如图35和图36所示，并进行有关参数的设置，如图37和图38所示。
设置场输出变量，如图39所示。
请求一个新的历史输出变量，如图40所示。
5t*ti SUT L C Ci:«ic, 仙*rd.I■ »>■''■ in h 讥 5■■IXiiMEtlrVal Snlul inr int F *1 i ： FdlitPT育 314 曲MG"'：■ Ft a^hfll I CT Mi pr IWL IQS I.I QO' ■>»■!■ 71『u*av..4«ir tv - 'U IHLT ^«tw—diEi »®d dtiFnLLiTvii 1 i *r i7EW»LiLrrV：Ti•昌©D LIIA SIR -IIPL Hi i-l Hhi-IFI mHiu-jAjia^niF * s'气 I IMP ■■»•*■ ■!■ Ja lAq -ai H ■■匕 JiTflJ? IZI* 茅 F>dd 汕即诫 EtaiiHtiiS y) 1 土 n-«lp1L!： l«I'L«r Si d'ln 日皿丄 A 也 A4申胖刊 klhH [otirii£liixnCiAac HE I Liu吝i CnwLralsL 3C.L kCL ki-LXL^LJ-ULJC-vi<ri.iii!» 慟CoUMCl■ uLtiLily J«Ji antl l! ■ Pdel 。
时命星花 FyMhisJrj -ni 兰L*皿加 ■ JHiis-ld Output K^uaerfts l^n-direE Y.3uJ- St-: 4,，F 孚h E ■! Id Oul - ptlk Kr sI F -Cutp^i-L Stb-Lic『Dy Output Jtil G:5OiiPrut Viri^Mts*) t J.4：t £CDO 11t.&!□% C_.' Fr^S-iliftct ad ^B XV A-'JE ? All IjCirzlie[呵h J ■徑Fzelr 血甘皿I 】) r■ 昭Nyx 朝but 1■空gpxt 嗖 pffi- -dur- P 寸■, "「E J皿 汇叱 | U KF. 5 STmEI ：.U □£• C±*£>Aildd i n. tKis. sl4>pivp A_jto S'w'f trei Y& C :iuL*Et«r AzmiRumd 龙j*d 凸鬧二护■山电Bo.luf.fl-J 孰C ) :2 iTiilidH ptBffl[■] ?MFf □「T WY 口 ①回j F Al 11 IT 昨| Fhar i >«1□ vawii m St«i<J I4«lThi 41s 1. L "代 4cr h uttcI IflE-ilil bl.l! L= BEAY r ：i« ku I ：I .M*d*l*j 口右F Tt-匠 ViterL dl =5 (1JvllHL F £丹号Um Q 甬"FM LL" I-非 Aw«-#.¥1 r十 *r 、© ii.nf— ERad-uiB 曰号£t«p -BOH-UT I BOH-UT IM Edie 191 Jt urp O UT 国亡町店:vlH L <JT T Q irfc pw * 耳导■■召时 E ： V SU * rUH"、 *% <Jk«Ll. 2*. Udl 丄屮tM-tiiH. f-sihlt七Load 模块中的操作：……E t A.t 1 FT SFILEIC . LutpX iwiift tana-EJ'h -s fHW.-T ：耳T 也节 r 1tr-eyi. 5 *1 * :Wi IJ]- □ 5tr*sS*t•□ sotlrwHJ 巩沖1狀•■虱,帕】ici曲】训枇1 Q 1-[Hl Tr.r -JIV J'E. - li i n- s*u-<>口a:pjA V«c x dblvaih<n lul b«lM«d»l«.ul«L 匚 Jdl.*» iAldm设置4个位移边界条件，位移值、边界条件的名称、类型、作用载荷步、作用区域等分别 如图41〜44所示。
