扩展有限元方法和裂纹扩展

合集下载

扩展有限元法(XFEM)及其应用12

扩展有限元法(XFEM)及其应用12

• 1999年,以美国西北大学Belytschko教授为代表 的研究组首先提出了XFEM的思想[21],2000年, 他们正式提出了XFEM术语[22]。XFEM是迄今为 止求解不连续力学问题的最有效的数值方法,它 在标准有限元框架内研究问题,不需要对结构内 存在的几何或物理界面进行剖分,保留了CFEM 的所有优点。XFEM与CFEM的最根本区别在于所 使用的网格与结构内部的几何或物理界面无关, 从而克服了在诸如裂纹尖端等高应力和变形集中 区进行高密度网格剖分所带来的困难,当模拟裂 纹扩展时也无需对网格进行重新剖分。XFEM在 处理裂纹问题包括以下三个方面[23]:
• 2. 单位分解法(PUM) • 2.1 单位分解法的基本概念 • 1996年Melenk和Babuska[24]及Duarte和 Oden[25]先后提出了单位分解法(PUM), 其基本思想是任意函数ψ(x)都可以用域内一 组局部函数NI(x)ψ(x)表示,即 • ,, (1) • 其中,NI(x)为有限单元形状函数, 它形成一 个单位分解。 N ( x) 1 • , (2) • 基于此,可以对有限元形状函数根据需要 进行改进。
• 数值方法,如有限元、边界元、无单元法等,特别是有限 元法(FEM)已被广泛用于处理不连续问题。有限元法具 有其它数值方法无可比拟的优点,如适用于任意几何形状 和边界条件、材料和几何非线性问题、各向异性问题、容 易编程等,是数值分析裂纹问题的主要手段。这方面的工 作很多,无法一一列举。Oritz等[1]及Belytschko等[2]通过 使用多场变分原理,用可以横贯有限单元的“弱”(应变) 间断模拟剪切带。Dvorkin等[3]通过修改虚功原理表达式 考虑了“强”(位移)间断问题;Lotfi和Sheng[4]将HuWashizu变分原理推广至具有内部间断的物体中;通过考 虑软化本构律和界面上的面力-位移关系,Simo及其同 事[5,6]提出了分析强间断问题的统一框架,很多研究者 [7-12]将该法应用到变形局部化分析中。Borja[13]提出了 分析强间断问题的标准Galerkin公式,并证明它与假定改 进应变逼近等价。

《2024年重型卡车紧急制动制动鼓有限元分析及裂纹扩展模拟》范文

《2024年重型卡车紧急制动制动鼓有限元分析及裂纹扩展模拟》范文

《重型卡车紧急制动制动鼓有限元分析及裂纹扩展模拟》篇一一、引言随着现代物流业和交通运输业的快速发展,重型卡车作为主要的运输工具之一,其安全性能尤为重要。

紧急制动系统是重型卡车安全保障的关键部分,而制动鼓作为该系统的主要构件,其性能的优劣直接关系到车辆的安全性。

因此,对重型卡车紧急制动制动鼓的力学性能及裂纹扩展进行研究具有重要的工程实际意义。

本文采用有限元分析方法,对重型卡车紧急制动时的制动鼓进行应力分析,并模拟裂纹的扩展过程,以期为制动鼓的设计与优化提供理论依据。

二、有限元模型建立1. 模型简化与假设为便于计算和分析,需要对实际制动鼓进行一定的简化。

假设制动鼓为均质材料,忽略其内部结构细节,并假定材料属性在各方向上均匀一致。

2. 材料属性定义采用有限元分析软件,定义制动鼓的材料属性,包括弹性模量、泊松比、屈服强度等。

同时,考虑材料在紧急制动过程中的非线性行为。

3. 网格划分对制动鼓进行网格划分,确保在关键区域如应力集中处有足够的网格密度,以捕捉到更精细的应力分布。

三、紧急制动过程模拟1. 加载条件设定根据实际工况,设定紧急制动的加载条件,包括制动力矩、车轮转速等。

同时考虑制动过程中可能出现的热效应对材料性能的影响。

2. 应力分析通过有限元分析软件进行计算,得到制动鼓在紧急制动过程中的应力分布情况。

分析最大应力出现的位置,以及应力随时间的变化情况。

四、裂纹扩展模拟1. 裂纹初始化在有限元模型中预设裂纹,根据实际损伤情况设定裂纹的初始尺寸和位置。

2. 裂纹扩展准则采用合适的裂纹扩展准则,如能量释放率准则或应力强度因子准则,来描述裂纹的扩展行为。

3. 裂纹扩展模拟在有限元模型中逐步推进裂纹的扩展过程,观察裂纹的扩展路径、扩展速度以及扩展过程中应力的变化情况。

五、结果与讨论1. 应力分析结果通过有限元分析得到制动鼓的应力分布云图,分析最大应力出现的位置及原因。

讨论应急制动过程中可能出现的应力集中现象。

2. 裂纹扩展模拟结果模拟得到的裂纹扩展路径与实际损伤情况相比较,验证模拟的准确性。

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

裂纹扩展的扩展有限元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积分等参量,也可以模拟裂纹的开裂过程。

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.固定裂纹中,只有各向同性材料的裂纹尖端渐进场才被考虑。

基于xfem的垂直于双材料界面的裂纹扩展问题

基于xfem的垂直于双材料界面的裂纹扩展问题

文章主题:基于xfem的垂直于双材料界面的裂纹扩展问题在材料科学和工程领域,裂纹扩展问题一直是一个备受关注的研究课题。

特别是在双材料界面上的裂纹扩展问题,由于双材料特性的不均匀性和复杂性,增加了研究和分析的难度。

在本文中,我们将从基于xfem(扩展有限元法)的角度出发,探讨垂直于双材料界面的裂纹扩展问题,以期为这一领域的研究和实际应用提供新的思路和方法。

一、概述垂直于双材料界面的裂纹扩展问题是指在两种材料的交界面上,裂纹在垂直方向上的扩展行为。

这种情况下,裂纹扩展的受力和受约束条件都受到了双材料特性的影响,需要深入分析和研究。

传统的有限元法在模拟和分析这种问题时存在一定的局限性,而xfem则能够有效地刻画裂纹的扩展路径和受力情况,因此成为了研究这一问题的有力工具。

二、裂纹模型的建立在进行垂直于双材料界面的裂纹扩展问题建模时,需要考虑双材料界面的影响,分析裂纹在材料间传播的受力情况和速度。

利用xfem,可以方便地将裂纹扩展路径和扩展速度等参数纳入模型中,通过数值计算得到裂纹扩展的演化规律和裂纹尖端的受力情况。

这有助于更准确地理解和分析垂直双材料界面上的裂纹扩展问题。

三、影响因素分析垂直于双材料界面的裂纹扩展受到诸多因素的影响,其中包括材料性质、裂纹尺寸、应力状态等。

通过xfem的数值模拟,可以分析不同因素对裂纹扩展行为的影响程度,揭示裂纹扩展过程中的关键因素和规律。

这有助于为材料设计和工程应用提供更可靠的参考依据。

四、工程应用与展望垂直于双材料界面的裂纹扩展问题在工程应用中具有重要意义,例如在复合材料结构的设计和评估中。

通过对裂纹扩展行为的深入研究和分析,可以为工程实践提供更有效的裂纹控制和结构安全性评估方法。

未来,基于xfem的裂纹扩展问题研究还可以结合人工智能算法和大数据分析等方法,进一步提高模型精度和计算效率,拓展应用领域和深化理论研究。

个人观点与总结垂直于双材料界面的裂纹扩展问题是一个复杂而又具有挑战性的研究领域,需要运用先进的数值模拟方法和理论分析手段来探讨和解决。

航空器复合材料胶接接头设计(ABAQUS-XFEM)

航空器复合材料胶接接头设计(ABAQUS-XFEM)

摘要复合材料结构的连接形式主要分为胶接和机械连接,随着复合材料在航空航天领域的广泛应用,胶接因其在复合材料结构连接中的优良特性日益受到结构设计人员的青睐,具有连接效率高、结构轻、抗疲劳、密封性好等优点。

然而胶接设计也具有很大的挑战性,在结构强度计算中,胶接连接接头部位一般为危险部位,需要重点校核。

所以,对复合材料胶接接头的设计分析是十分必要的。

本选题利用成熟的有限元商用软件ABAQUS,使用XFEM(扩展有限元法)对胶层和复合材料层的应力场等进行分析。

通过分析计算这些应力,同时应用相应的失效准则,进而可预测初始裂纹的扩展与否及扩展的长度,为胶接接头设计的选择提供必要的依据。

在文章中,讨论了胶接长度、胶层厚度和初始裂纹的位置对裂纹扩展的影响。

通过对仿真结果的分析,提出了减小胶接长度和胶层厚度的观点,指出裂纹易于产生及扩展的区域,对胶接接头的设计进行了优化。

胶接接头的优化设计对拓宽复合材料在飞机结构上的应用范围,进一步减轻结构重量、提高疲劳性能和降低制造成本具有重要的工程使用价值。

关键词:复合材料板胶接接头扩展有限元裂纹扩展AbstractThe joint methods of composite structure contain cementing and mechanical connection.. With the use of composite in the field of aviation increased a lot in recent years for its high strength and lightness, the cementing is increasingly favored by the structure design staff for its excellent characteristics in the connection field of composite structure. The characteristics are high ligation efficiency, light structure, antifatigue and good sealing. However, glued design also has a great challenge. In the structural strength calculations, glued joints are generally connected to dangerous parts and need to focus on checking. Therefore, the design and analysis of composite bonded joint is very necessary.The topic use the sophisticated and commercial software -ABAQUS, in the field of finite element, and use XFEM ( extended finite element method ) as the foundation to analysis the stress field of bonding layers and composite layers. By analyzing and calculating these stresses, while applying the appropriate failure criterion, we can predict the initial crack extension and the length of the expansion. In this way, it can provide the necessary basis for the design of bonding joints. In the article, we discussed the impact of the bonding length, layer thickness and initial crack location on crack propagation. Through the analysis of simulation results, we presented two standpoints of reducing the length of bonding joint and the thickness of adhesive. Besides, we pointed the areas where cracks are easy to generate and expand. Optimal design of adhesive joints in composite materials has important engineering value to broaden the scope of application of the aircraft structure and further reduce the structural weight, improve the performance of fatigue and reduce manufacturing costs.Keywords:Composite plates, Adhesive joints, XFEM, Crack extension目录摘要 (I)Abstract ....................................................... I I 目录.......................................................... I II 第一章引言.. (1)1.1导言 (1)1.2胶接连接 (2)1.2.1 简介 (2)1.2.2胶接连接应当注意的问题 (3)1.2.3胶接连接研究现状 (3)1.3 胶接接头 (4)1.3.1胶接接头简介 (4)1.3.2胶接接头的基本形式 (5)1.3.3胶接接头的破坏模式 (6)1.3.4胶接接头处可能出现的裂纹及其影响 (7)第二章复合材料损伤和胶接连接的力学模型 (8)2.1导言 (8)2.2复合材料层板强度预测 (8)2.3复合材料和胶层断裂准则 (10)第三章利用ABAQUS建立复合材料胶接接的有限元模型 (13)3.1扩展有限元方法和工程软件ABAQUS简介 (13)3.1.1传统有限元方法 (13)3.1.2扩展有限元方法及基本原理 (14)3.1.3ABAQUS简介 (15)3.2利用ABAQUS建立复合材料板胶接模型的过程 (16)3.2.1几何模型的建立和约束条件 (16)3.2.2材料属性 (17)3.2.3定义接触 (19)3.2.4 对于XFEM定义 (19)第四章基于裂纹扩展分析的单面搭接接头设计 (21)4.1复合材料胶接接头在纵向载荷下的受力分析 (21)4.2不同搭接长度下胶接接头的裂纹扩展情况 (23)4.2.1搭接长度为15mm的情况 (23)4.2.2搭接长度为10mm的情况 (25)4.2.3搭接长度为20mm的情况 (26)4.2.4不同搭接长度下裂纹情况的对比及结论 (28)4.3不同胶层厚度下胶接接头的裂纹扩展情况 (29)4.3.1胶层厚度为0.1mm的情况 (29)4.3.2胶层厚度为0.2mm的情况 (31)4.3.3胶层厚度为0.3mm的情况 (33)4.3.4不同胶层厚度下裂纹情况的对比及结论 (34)带五章基于裂纹扩展的斜面搭接接头设计 (37)5.1斜面搭接接头在纵向载荷下的受力分析 (37)5.2不同裂纹位置下胶接接头的裂纹扩展情况 (38)5.2.1选取的三种不同裂纹位置 (39)5.2.2裂纹的扩展情况 (40)5.2.3三种情况对比及结论 (42)5.3单面搭接和斜面搭接情况的对比 (43)第六章全文总结及展望 (46)6.1全文总结 (46)6.2展望 (47)致谢辞 (49)参考文献 (50)第一章引言1.1导言复合材料作为一种新材料,在最近的半个多世纪中飞速发展,由于复合材料采用纤维加强结构,使得复合材料具有比重小、比强度和比模量大的特点,并且由于采用的是铺层结构,制造过程简单,容易成型。

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

裂纹扩展的扩展有限元(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积分等参量,也可以模拟裂纹的开裂过程。

基于扩展有限元的碳纤维复合材料裂纹扩展仿真

基于扩展有限元的碳纤维复合材料裂纹扩展仿真

基于扩展有限元的碳纤维复合材料裂纹扩展仿真韩少燕 门 静 韩海燕(西安交通大学城市学院,陕西 西安 710018)引言 碳纤维复合材料以其良好的力学性能被广泛的应用于汽车、航空航天等领域[1]。

碳纤维层合板在实际使用过程中容易受到冲击载荷产生大变形弯曲,导致局部产生应力集中与应变从而引起材料损伤,例如基体开裂、纤维断裂后或者层间分层等,材料损伤扩展会进一步导致力学性能降低,从而导致材料失效最终结构失效。

扩展有限元通过引入富集函数来修正传统有限元的近似位移函数,以描述间断界面,使间断的描述独立于有限元网格,避免了计算过程中的网格重构[2]。

本文采用扩展有限元法模拟了碳纤维复合材料层合板在弯曲载荷作用下的开裂过程,以预测材料抵抗外力损伤的性能。

1、扩展有限元 扩展有限元是以美国西北大学Belytschko 教授为首的研究组于1999年提出的一种求解不连续问题的数值方法,该方法可有效的求解强和弱不连续问题[2-3]。

扩展有限元的基本原理是基于单位分解法在传统有限元位移模式中加入特殊函数(加强函数),从而反应不连续性的存在,不同类型的不连续问题,只是加强函数不同而已。

1.1单位分解法单位分解法是Melenk 和Bubska 及Duarte 和Oden 于1996年先后提出的。

对于求解区域Ω,单位分解法用一些相互交叉的子域ΩI 来覆盖,每个子域都与一个函数()I ϕx 相联系。

函数()I ϕx 仅在ΩI 内非零,且满足单位分解条件()1I Iϕ=∑x (1)Duarte 和Oden 用K 阶移动最小二乘近似函数来构造单位分解,即1()()[()]mh k I iI i Ii b q ϕ==+∑∑ u x x u x (2) 其中:()i q x 可以是单项式基。

系数是未知量,可以通过Galerkin 法或配点法求解。

为了提高逼近精度,或满足对待定问题的特殊逼近要求,也可以包含其他一些形式的函数(称之为加强基函数)。

扩展有限元方法及应用综述_郭历伦

扩展有限元方法及应用综述_郭历伦
采用水平集方法后富集函数所采用的极坐标系由水平集函数定义见式14同时计算j积分时由水平集函数定义虚速度矢量14belytschko和blac以及mos等的研究均是基于准静态二维裂纹问题在提出xfem方法时并未涉及裂纹扩展因此仅采用有向距离函数就可以描述裂纹的几何特性无需求解式17sukumar等采用快速推进法研究了三维平面裂纹扩展问题并在后续工作中逐步研究了多个平面裂纹的扩展问题及三维非平面裂纹扩展问题stolarsko等10采用水平集方法研究二维裂纹扩展问题mos位错问题fig5twolevelsetfunctions1112将水平集方法应用于三维非平面裂纹问题中ventura等13提出了矢量水平集vectorlevelset14算法
[ 3~4]

。 对于任意函数 V( , 可得 x)
I=1
。 x) V( x) =V( x) ∑ (

设函数 VI( 为函数 u( 在子域 ΩI 内的局部近似函数 , 则函数 u( 在求解域的全局近似可取为 x) x) x)
1 单位分解函数
由于扩展有限元近似函数的基础是单位分解法 , 本节将简要介绍 单 位 分 解 法 。 单 位 分 解 法 使 用 一 些 , 。 在每个子域 ΩI 上定义一个仅在该子 域 以节点 x I 为中心的子域 Ω I 来覆盖整个求解区域 即 Ω ∪Ω I
I=1 N N () , () } 内非零的函数 并且要求它们满足单位分解条件 : x) =1。 则函数集 { I=1 称为属 于开 I x I x ∑ I( I=1 N N 覆盖 { ΩI} I=1 的单位分解函数
( , , ) I n s t i t u t e o f S s t e m s E n i n e e r i n C A E P, P. O. B o x 9 1 9 4 1 1, M i a n a n S i c h u a n 6 2 1 9 0 0, C h i n a - y g g g y g

XFEM实现裂纹扩展

XFEM实现裂纹扩展

---因为专注,所以卓越!
网格划分
焊缝在管道的上下起始位置,造成几何模型的急剧变化, 导致网格不容易划分,因此,使用专业的的前处理软件 ANSA进行网格划分,使得焊缝的网格密度大于其他位置 的网格密度。
初始裂纹在焊缝中的位置

---因为专注,所以卓越!
分析过程

---因为专注,所以卓越!
I型裂纹扩展过程的动画演示

---因为专注,所以卓越!
II型裂纹扩展过程的动画演示

---因为专注,所以卓越!
ABAQUS采用XFEM模 案例2 块实现压力容器的裂纹 过程的模拟,如果图所 示,压力容器与外部连 接的接口处存在初始微 裂纹,当容器内压力达 到一定程度,裂纹开始 启裂并扩展。 模型的建模与应用针对 工程实例,很好的展现 了XFEM强大的裂纹扩 展功能。
石油管道的裂纹扩展模拟
利用ABAQUS的XFEM方法实现石油管道的裂纹扩展,在 已知起始裂纹尺寸的情况下,根据外部载荷模拟裂纹的起 裂和扩展过程。 由于裂纹的尺寸较整体模型尺寸较小,因此采用用户子模 型的方法对局部进行更加细致的分析。
一、XFEM模块功能简述
ABAQUS V6.9及其以后的版本将拓展有限元方法引入到 其分析中,并增加了新的模块XFEM,该方法可以认为是 有限元方法处理不连续问题的革命性变革。这是第一个将 XFEM商用化的软件。 固体力学中存在两类典型的不连续问题,一类是因材料特 性突变引起的弱不连续问题,这类问题以双材料问题和夹 杂问题为代表,其复杂性由物理界面处的应变不连续性引 起;另一类是因物体内部几何突变引起的强不连续问题, 这类问题以裂纹问题为代表,其复杂性由几何界面处的位 移不连续性和端部的奇异性引起。物体内部物理界面的脱 粘或起裂,是上述两类问题的混合。

重力坝开裂过程扩展有限元数值模拟

重力坝开裂过程扩展有限元数值模拟

重力坝开裂过程扩展有限元数值模拟靳旭;董羽蕙【摘要】扩展有限元法(XFEM)是一种求解不连续问题的数值方法.它继承了常规有限元法(CFEM)的所有优点,在模拟裂纹扩展、界面、复杂流体等不连续问题时特别有效,近十多年得到了快速发展.介绍了XFEM的基本原理,给出了进行混凝土裂纹扩展分析的方法.利用XFEM模拟混凝土重力坝裂纹扩展,通过对比有、无裂纹情况下的重力坝应力分布,分析裂纹存在对重力坝应力场分布的影响;分析裂纹扩展受网格疏密程度的影响;计算在不同岩基弹性模量下裂纹的扩展方向.%Extended finite element method(XFEM)is a numerical solution for analyzing discontimuity problem . It inherited all the advantages of the conventional finite element method (CFEM) , in the simulation of crack extension , interface, complex fluid and other discontinuities are particularly effective , in the past decade it has been rapid development. The basic theory of XFEM in introduced and the method of analyzing concrete fracture is presented. The XFEM is utilized to simulate the crack propagation in concrete gravity dam. By the contrast of stress distribution under no crack and crack circumstance of gravity dam the discipline of stress field distribution is analyzed; It is also used for influence of mesh density to crack propagation and is calculated the crack propagation direction in batholith elastic modulus.【期刊名称】《科学技术与工程》【年(卷),期】2012(012)033【总页数】6页(P9100-9104,9109)【关键词】重力坝;扩展有限元法;裂纹扩展;网格疏密;弹性模量【作者】靳旭;董羽蕙【作者单位】昆明理工大学建筑工程学院,昆明650500;昆明理工大学建筑工程学院,昆明650500【正文语种】中文【中图分类】TV313;TV642.3实际工程中,无论采用多么严格的裂缝控制措施,混凝土结构仍然会带裂缝工作。

热腐蚀疲劳_裂纹扩展_有限元模拟_概述说明以及解释

热腐蚀疲劳_裂纹扩展_有限元模拟_概述说明以及解释

热腐蚀疲劳裂纹扩展有限元模拟概述说明以及解释1. 引言1.1 概述热腐蚀疲劳是金属材料在高温和腐蚀环境中长期使用过程中发生的一种疲劳现象。

在这种条件下,材料会受到复杂的力学、化学和热耦合作用,导致裂纹的产生和扩展。

因此,对于热腐蚀疲劳和裂纹扩展机制进行深入的研究具有重要意义。

1.2 文章结构本文将分为三个部分进行讨论,首先会对热腐蚀疲劳进行概述说明,包括其产生原因、影响因素等内容;然后将解释裂纹扩展机制,阐述不同环境条件下裂纹扩展行为的特点;最后介绍有限元模拟在热腐蚀疲劳中的应用,并探讨其优缺点以及未来发展方向。

1.3 目的本文旨在全面了解和阐述热腐蚀疲劳、裂纹扩展机制以及有限元模拟在该领域中的应用。

通过对已有文献和实验结果的综合分析,探讨热腐蚀疲劳对材料性能和工程结构寿命的影响,并为相关领域的进一步研究提供参考和展望。

2. 正文:2.1 热腐蚀疲劳概述说明:热腐蚀疲劳是指在高温和腐蚀介质中,材料发生的由于循环载荷引起的疲劳断裂现象。

它是一种常见的工程问题,广泛存在于航空、能源等领域。

热腐蚀疲劳具有复杂的损伤机制,包括氧化、腐蚀、应力诱导等多个因素相互作用所导致的材料失效。

2.2 裂纹扩展机制解释:裂纹扩展是热腐蚀疲劳过程中的一项重要现象。

在循环载荷作用下,由于应力集中或缺陷存在,材料表面产生微小裂纹。

随着时间和循环次数的增加,这些微小裂纹会逐渐扩展并相互连通,最终导致材料失效。

裂纹扩展受到多个因素影响,如应力水平、温度、腐蚀介质以及材料本身的力学性能等。

2.3 有限元模拟在热腐蚀疲劳中的应用:有限元模拟是一种常用的工程分析方法,可以应用于热腐蚀疲劳相关问题的研究。

通过建立适当的数学模型、物理方程和边界条件,有限元模拟可以在计算机上模拟材料在不同工况下的响应行为。

在热腐蚀疲劳中,有限元模拟可以用来预测材料的失效寿命、裂纹扩展速率以及评估其结构的可靠性。

有限元模拟在热腐蚀疲劳中的应用主要包括以下几个方面:首先,通过建立材料力学性能及裂纹扩展规律等方面的数学模型,可以使用有限元模拟来预测热腐蚀疲劳过程中裂纹扩展的发展趋势和速率。

扩展有限元法在疲劳裂纹扩展模拟中的应用

扩展有限元法在疲劳裂纹扩展模拟中的应用

2019年18期创新前沿科技创新与应用Technology Innovation and Application扩展有限元法在疲劳裂纹扩展模拟中的应用张芮晨(中国直升机设计研究所,天津300000)前言目前,有限元仿真大量应用于疲劳问题分析计算中。

然而常规有限元法采用连续函数作为形状(插值)函数,在单元内部形状函数和材料性能必须保证连续,进行裂纹分析时,裂纹面需与网格边界重合,网格节点需与间断面的尖端重合[1-2]。

如Abaqus 中针对裂纹扩展分析的两种技术:基于debond 裂纹扩展技术和基于cohesive 裂纹扩展技术。

两种技术均需提前指定裂纹扩展区域与方向,对于网格的划分限制较大,工作量大,效率低。

而在扩展有限元法(Extended Finite Element Method ,XFEM )中,网格与结构内部的几何或物理界面无关,两者间相互独立,裂纹尖端应力场的计算与裂纹面扩展计算相互独立,避免了在裂纹尖端进行高密度网格划分的问题和网格重新划分的工作[3-5]。

本文采用扩展有限元法对疲劳裂纹扩展进行仿真分析,达到无需指定裂纹扩展路径、无需进行高密网格划分也可模拟裂纹扩展的目的,节省建模所需要的时间,大大减少工作量,同时又提高了裂纹扩展模拟的准确度。

1扩展有限元法理论基础1.1单位分解法(PUM )扩展有限元通过在常规有限元位移模式中加进跳跃函数和渐进裂尖位移场函数等一些特殊的函数,来有效地模拟类似裂纹一类的强不连续问题,而它的理论基础既为单位分解法。

在单位分解法中,认为任何函数φ(x)都可以用域内局部函数H Γ(x)φ(x)表示,即其中,H Γ(x)为有限单元形状函数,它形成一个单位分解(2)以此为基础,对有限元形状函数进行改进。

重叠分片{θi }构成研究区域M 的一个覆盖,{φi }为覆盖上的一个单位分解。

在每一分片上,函数空间V i 为区域M 的局部逼近,总体试探空间V 为(3)总体空间V 不但具有局部空间V i 的逼近特性,又有单位分解φi 和局部空间V i 的光滑性,只要单位分解φi 足够光滑,就能构造出足够光滑的试探空间。

ANSYS16.0新增扩展有限元XFEM裂纹扩展仿真简介

ANSYS16.0新增扩展有限元XFEM裂纹扩展仿真简介

ANSYS16.0新增扩展有限元XFEM裂纹扩展仿真简介中国矿业大学, 师访, matmes@1 引言早在两年前,就听安世亚太的人说ANSYS15.0将加入XFEM,但结果令人失望。

左盼右盼,终于在ANSYS16.0中等来了扩展有限元(Extended Finite ElementMethod)XFEM 功能。

首先,对于不知道XFEM为何物的朋友们,建议看下这篇文献:断裂问题的扩展有限元法研究_茹忠亮_岩土力学_2011.pdf(834.69 KB, 下载次数: 0)欢迎联系我讨论关于XFEM断裂模拟的相关问题,QQ:15492217582 ANSYS16.0 XFEM简介ANSYS16.0发布时候关于结构分析的简介中并未提及XFEM,由此可见ANSYS中的XFEM功能也不会太让人满意,看了其帮助文档(ANSYS Mechanical APDL Fracture Analysis Guide.pdf(2.52 MB, 下载次数: 0),3.2节)后发现,事实确实如此。

(1)基于虚拟节点法,与Abaqus一致。

(2)仅支持线弹性材料。

(3)不支持裂尖增强,同样与Abaqus一致,这就导致裂纹尖端不能落在单元内部,只能位于单元边界上。

(4)必须事先给定初始裂纹,即不支持裂纹的自动萌生。

鸡肋的是,初始裂纹的定义居然要通过给定水平集值的方法来实现:XFDATA,LSM,ELEMNUM,NODENUM,PHI(5)支持粘聚裂纹。

(6)支持PLANE182(4节点四边形单元,用于2D平面问题分析)及SOLID185(8节点正方形单元,用于3D问题分析)这两种单元。

(7)仅支持准静态分析,不支持动态断裂。

(8)裂纹每次只能扩展一个单元长度。

(9)仅支持两个裂纹扩展准则:STTMAX,最大周向应力准则(Maximum circumferential stress criterion);PSMAX(Circumferential stress criterion)。

基于扩展有限元和循环内聚力模型的榫连结构微动疲劳裂纹萌生与扩展研究

基于扩展有限元和循环内聚力模型的榫连结构微动疲劳裂纹萌生与扩展研究

基于扩展有限元和循环内聚力模型的榫连结构微动疲劳裂纹萌生与扩展研究随着机械设备的广泛应用,微动疲劳问题越来越受到关注。

微动疲劳是指由于微小振动和相对运动而引起的表面疲劳现象。

榫连结构是一种常见的连接方式,但在实际应用中,榫连结构容易受到微动疲劳的影响,导致结构失效。

因此,研究榫连结构微动疲劳裂纹萌生与扩展具有重要的意义。

传统的微动疲劳研究方法主要是基于实验和有限元分析,但这些方法存在一定的局限性。

实验方法需要大量的时间和精力,而且结果受到许多因素的影响,如温度、湿度等。

有限元分析虽然可以模拟微动疲劳过程,但往往忽略了材料的微观结构和力学性能的差异。

因此,需要采用更加精确的方法来研究榫连结构微动疲劳裂纹萌生与扩展。

基于扩展有限元和循环内聚力模型,我们可以对榫连结构微动疲劳裂纹萌生与扩展进行更为精确的研究。

扩展有限元方法是一种数值模拟方法,可以模拟复杂几何形状和材料属性的微动疲劳过程。

循环内聚力模型是一种考虑材料微观结构和力学性能差异的模型,可以更好地描述微动疲劳过程中的裂纹萌生和扩展。

在本文中,我们首先对榫连结构微动疲劳裂纹萌生与扩展的理论基础进行了介绍,包括微动疲劳的基本概念、扩展有限元方法和循环内聚力模型的应用。

接着,我们建立了榫连结构微动疲劳的数值模型,并进行了参数分析和优化设计。

通过模拟和分析,我们发现材料的微观结构和力学性能、相对运动速度、表面状态等因素对微动疲劳裂纹萌生与扩展有显著影响。

此外,我们还发现扩展有限元和循环内聚力模型可以有效地模拟榫连结构微动疲劳裂纹萌生与扩展过程,并能够为实际应用提供指导。

例如,可以通过优化结构设计、选择合适的材料和表面处理方法来降低榫连结构的微动疲劳风险。

总之,基于扩展有限元和循环内聚力模型的榫连结构微动疲劳裂纹萌生与扩展研究具有重要的理论和实践意义。

通过这种方法,我们可以更好地了解微动疲劳过程,为实际应用提供更加精确的指导。

未来,我们将继续研究榫连结构微动疲劳裂纹萌生与扩展的其他影响因素,并探索更加有效的预防和修复方法。

Abaqus提取XFEM(扩展有限元)裂缝长度和缝隙面积(精品)

Abaqus提取XFEM(扩展有限元)裂缝长度和缝隙面积(精品)

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。

扩展有限元方法和裂纹扩展

扩展有限元方法和裂纹扩展

扩展有限元方法和裂纹扩展1.1扩展有限元方法(XFEM)基本理论1999 年,美国Northwestern University 的Belytschko 和Black 领导的研究小组提出了扩展有限元方法,为解决裂纹这类强不连续问题带来了曙光。

他们正式应用扩展有限元法(XFEM)这一专业术语是在2000年,截止到目前,扩展有限元法(XFEM)成为我们解决强不连续力学问题的最有效的数值计算方法,也成为计算断裂力学的重要分支。

XFEM在有限元的框架下进行求解,无需对构件内部的物理界面进行网格划分,具有常规有限元方法的所有优点。

它最明显的特点是用已知的特征函数作为形函数来使传统有限元的位移得到逼近,进而克服了在裂纹尖端和变形集中处进行高密度网络划分产生的困难,方便地模拟裂纹的任意路径,而且计算精度和效率得到了显著的提高⑹。

扩展有限元方法是将已知解析解的特征函数作为插值函数增强传统有限元的位移逼近,来使得单元内的真实位移特性得以体现,裂纹尖端和物理或儿何界面独立于有限元网格。

XFEM主要包括以下三部分内容:首先是不考虑构件的任何内部细节,按照构件的儿何外形尺寸生成有限元网格;其次,采用水平集方法跟踪裂纹的实际位置;根据已知解,改进影响区域的单元的形函数,来反映裂纹的扩展。

最后通过引入不连续位移模式来表示不连续儿何界面的演化。

因为改进的插值函数在单元内部具有单元分解的特性,其刚度矩阵的特点与常规有限元法的刚度矩阵特性保持一致。

单元分解法(Partition Of Unity Method)和水平集法(Level Set Method).节点扩展函数构成了扩展有限元法的基本理论,其中,单元分解法是通过引入加强函数讣算平面裂纹扩展问题,保证了XFEM的收敛性; 水平集法是跟踪裂纹的位置和模拟裂纹扩展的常用数值方法,任何内部儿何界面位置都可用它的零水平集函数来表示。

(1)单元分解法的基本思想是任意函数处工)都可以用子域内一组局部函数Ng(x)表示,满足如下等式:此丫)=工"应沁) ⑴I其中,它们满足单位分解条件:玄/%丫) = 1 M⑴是有限元法中的形函数,根I据上述理论,便可以根据需要对有限元的形函数进行改进。

飞行器机身结构的疲劳与断裂行为分析方法

飞行器机身结构的疲劳与断裂行为分析方法

飞行器机身结构的疲劳与断裂行为分析方法飞行器机身结构是整个飞行器系统中最关键的组成部分之一。

由于长期的运行和外界环境的影响,机身结构的疲劳与断裂问题可能会导致严重的安全事故。

因此,对于飞行器机身结构疲劳与断裂行为的分析方法具有重要的意义。

本文将介绍飞行器机身结构疲劳与断裂行为的分析方法,主要包括负载谱法、应力范围法、裂纹扩展法和有限元法。

1. 负载谱法负载谱法是一种通过对机身结构的负载进行统计和分析得出的疲劳寿命预测方法。

通过监测和记录机身结构在实际飞行中所受到的压力、力矩和振动等负载情况,建立负载谱。

然后,根据疲劳试验结果建立的S-N曲线,结合负载谱,可以推算出机身结构的疲劳寿命。

2. 应力范围法应力范围法是一种通过对机身结构中应力范围进行分析得出的疲劳寿命预测方法。

应力范围是指机身结构在一个完整的负荷周期内,应力的最大值和最小值之差。

通过对机身结构的应力范围进行测量和分析,对其进行统计性处理,可以得到疲劳寿命预测结果。

3. 裂纹扩展法裂纹扩展法是一种通过对机身结构中裂纹扩展的分析得出的疲劳寿命预测方法。

通过在机身结构中制造一定长度的初始裂纹,并以一定的载荷条件下进行疲劳试验,测量和记录裂纹扩展的速率,然后推算出机身结构的疲劳寿命。

4. 有限元法有限元法是一种通过使用数值计算方法对机身结构进行分析的方法。

该方法基于有限元理论,在计算机上建立机身结构的有限元模型,并通过应用适当的加载条件进行仿真分析。

通过分析机身结构在不同加载条件下的应力和变形分布,可以得出机身结构的疲劳寿命。

综上所述,飞行器机身结构的疲劳与断裂行为分析方法主要包括负载谱法、应力范围法、裂纹扩展法和有限元法。

这些方法可以帮助工程师预测和评估机身结构的疲劳性能,提前采取措施避免机身结构的疲劳破坏,确保飞行器的安全运行。

在实际应用中,可以根据具体情况结合使用这些方法,以获得更准确可靠的分析结果。

直缝焊管裂纹扩展过程有限元数值模拟

直缝焊管裂纹扩展过程有限元数值模拟
第3 4卷
第 2期








V o 1 . 3 4 N o . 2
A p r . 2 0 1 3
2 0 1 3年 4月
J O U R N A L O F T A I Y U A N U N I V E R S I T Y O F S C I E N C E A N D T E C HN O L O G Y
重要。
近, 但也不宜太远 , 因为式 ( 1 )已略去 r 的高阶项 。
设 裂纹 表 面两 点距 裂 纹尖 端 r 和r 2 两点 位 移分 别
为/ 3 和 , 计算
由于有 限元 方法 的发展 日趋 成 熟 , 运 用 有 限 元
理论 , 以非线性有限元 分析软件 A B A Q U S为平台 ,
作者简 介: 彭英 ( 1 9 7 9一) , 女, 讲师 , 主要研究 方向为裂 纹扩 展。
第3 4卷
第 2期








V o 1 . 3 4 N o . 2
A p r . 2 0 1 3
2 0 1 3年 4月
J O U R N A L O F T A I Y U A N U N I V E R S I T Y O F S C I E N C E A N D T E C H N O L O G Y
文章 编 号 : 1 6 7 3— 2 0 5 7 ( 2 0 1 3 ) 0 2—0 1 4 7—0 5
无 限大 功 能 梯 度 材 料 反 平面 裂 纹 应 力 场
胡 志新 , 张雪霞, 李 婵
( 太原科 技 大 学应 用科 学 学院 , 太原 0 3 0 0 2 4 )
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

扩展有限元方法和裂纹扩展
1.1 扩展有限元方法(XFEM )基本理论
1999年,美国Northwestern University 的Belytschko 和Black 领导的研究小
组提出了扩展有限元方法,为解决裂纹这类强不连续问题带来了曙光。

他们正式
应用扩展有限元法(XFEM )这一专业术语是在2000年,截止到目前,扩展有
限元法(XFEM )成为我们解决强不连续力学问题的最有效的数值计算方法,也
成为计算断裂力学的重要分支。

XFEM 在有限元的框架下进行求解,无需对构件
内部的物理界面进行网格划分,具有常规有限元方法的所有优点。

它最明显的特
点是用已知的特征函数作为形函数来使传统有限元的位移得到逼近,进而克服了
在裂纹尖端和变形集中处进行高密度网络划分产生的困难,方便地模拟裂纹的任
意路径,而且计算精度和效率得到了显著的提高[6]。

扩展有限元方法是将已知解析解的特征函数作为插值函数增强传统有限元
的位移逼近,来使得单元内的真实位移特性得以体现,裂纹尖端和物理或几何界
面独立于有限元网格。

XFEM 主要包括以下三部分内容:首先是不考虑构件的任
何内部细节,按照构件的几何外形尺寸生成有限元网格;其次,采用水平集方法
跟踪裂纹的实际位置;根据已知解,改进影响区域的单元的形函数,来反映裂纹
的扩展。

最后通过引入不连续位移模式来表示不连续几何界面的演化。

因为改进
的插值函数在单元内部具有单元分解的特性,其刚度矩阵的特点与常规有限元法
的刚度矩阵特性保持一致。

单元分解法(Partition Of Unity Method)和水平集法
(Level Set Method )、节点扩展函数构成了扩展有限元法的基本理论,其中,单
元分解法是通过引入加强函数计算平面裂纹扩展问题,保证了XFEM 的收敛性;
水平集法是跟踪裂纹的位置和模拟裂纹扩展的常用数值方法,任何内部几何界面
位置都可用它的零水平集函数来表示。

(1)单元分解法的基本思想是任意函数()x φ都可以用子域内一组局部函数
()()x x N I ϕ表示,满足如下等式:
()()()x x N x I
I ϕφ∑= (1)
其中,它们满足单位分解条件:f I I
åx ()=1 ()x N I 是有限元法中的形函数,根
据上述理论,便可以根据需要对有限元的形函数进行改进。

在XFEM 中,单元
分解的目的是进行数值积分,达到不引人额外的自由度的目的[7-8]。

(2)水平集法 使用水平集法来描述几何间断性。

在一般情形下,多用来追踪
界面的位置。

这里的界面可以是闭合曲线或者是计算域边界的曲线。

在扩展有限
元方法中,由于网格划分并不需要符合裂纹的几何性质,简化裂纹跟踪位置的关
键是对裂纹的几何描述。

而水平集法的强大优势是可以用于分析和计算界面运
动。

正好符合扩展有限元方法的要求即对于任意方向的裂纹增长不需要网格重
划。

裂纹的水平集函数常取下列符号的距离函数即:
()()
γγγψx x t x t x -±=∈min , (2) 若x 位于裂纹上方,上式符号为正,反之为负[9-10]。

(3)扩展有限元法的位移模式(节点扩展函数简介):
扩展有限元法的基本内容就是在含有位移不连续性的影响区域内通过一些
附加的加强函数来多常规有限元法的位移模式进行改进,进而对不连续的位移进
行描述。

为了实现断裂分析,裂纹尖端附近渐近函数用来模拟裂尖的应力奇异性,
间断函数则用来表示裂纹几何界面处的位移间断跳跃。

所以改进的位移逼近可以
写为:
()()()⎥⎦⎤⎢⎣⎡++=∑∑==411ααααI I I N
I I b x F x H u x N u (3)
其中,()x N I 为有限元形函数,I u 为节点位移向量的连续部分,I α为与
Heaviside 函数相关的节点扩展自由度向量,()x H 是裂纹面的间断跳跃函数,α
I b 是与弹性渐进裂尖函数有关的节点扩展自由度向量,()x F α是裂尖应力渐进函数。

等式右边第一项可用于构件中的所有节点,右边第二项单元节点形函数被裂纹完
全贯穿的情况,第三项裂尖停留在单元内部的节点。

引入位移扩展函数使得扩展有限元法确保裂尖可以被模拟为精确地停留在
单元内部,同时允许在相对粗糙的二维有限元网格上获得比较好的精度[11]。

1.2扩展有限元方法研究现状
1.1.1 有限元法研究裂纹扩展
有限元方法是目前来说应用非常广泛的数值计算方法之一,对弹性体进行模
拟裂纹扩展时,有限元网格的划分对整个计算过程最为重要,它不仅仅影响着计
算的精度和速率,而且还决定着此计算过程能否顺利进行。

模拟裂纹扩展问题的
关键是裂纹尖端处的应力场,所以裂纹尖端处的单元应当设计为奇异性单元,并
且在尖端领域内进行局部加密来提高计算精度,计算网格随着裂纹的扩展而不断
变化,用来保证奇异性单元和加密网格始终位于裂纹尖端[12] ,综上,网格的重新划分在有限元法模拟裂纹扩展过程中是非常关键的一个问题。

很多文献关于有限元方法模拟裂纹扩展过程的裂纹尖端处的网格生成进行了研究,比如杨庆生的《断裂过程的有限元模拟》中研究了自身适应的模拟裂纹扩展的有限元网格生成技术。

还有2003年陈永强和姚振汉在有限元技术的基础上,对表征非均匀材料且初始分布并不均匀的格子模型用重复多子域边界元法对非均匀材料在简单荷载作用下的破坏过程进行了相关研究和模拟。

诸如对裂纹扩展过程模拟的研究还有很多,但是有限元法是基于网格的数值方法,由于裂纹尖端的奇异性和强不连续性,使得网格划分变得困难,还是没有完全解决此类问题。

1.1.2扩展有限元方法研究裂纹扩展
扩展有限元方法(the extended finite method—XFEM)可以减少裂纹面网格划分带来的劣势。

该方法基于整体划分(单元分解)的概念,仍然属于传统有常规有限元方法的扩展。

整体划分的概念使扩展函数方便地进入到有限元当中,不连续性也可以通过引入额外自由度的扩展函数来解决,即可以确定裂纹的实际位置,跟踪裂纹的生长[13]。

当然,扩展有限元方法保留了传统有限元的框架,其刚度矩阵具有与常规有限单元一样的优点。

具有以下特点:
(1)单元划分的常规有限元扩展
(2)引入特殊的位移函数,扩展自由度实现不连续性
(3)不需要重新划分网格描述几何间断性
(4)所划分网格与构件的几何或物理界面无关
(5)可用于模拟任意性、求解相关路径裂纹的裂纹初始和扩展过程
(6)同时允许几何和材料的非线性的存在
(7)裂纹的生成与扩展路径可以完全基于计算的结果
(8)简化模型网格细化后的裂纹研究
(9)在裂尖应变奇异性改进方法的基础上,使裂纹分析的收敛速度得到提高
1.3断裂力学中裂纹相关内容
断裂力学的任务是:获得不同材料的断裂强度;预测物体在给定外力作用下是否将发生断裂,即建立材料断裂准则;研究载荷作用过程里的裂纹扩展规律;研究在复杂环境和应力同时作用下物体的断裂问题。

裂纹扩展有三种基本形式:张开型、滑开型以及撕开型。

在工程应用与实际
生产中,并非所有的断裂问题都可简化为上述三个基本模式之一,还有复合型断
裂是三种基本模式的组合,也是最常见的裂纹类型。

在断裂力学中,有最为基本
的理论即:能量释放率理论、裂纹尖端弹性应力场理论、J 积分理论[14]。

(1)能量释放率(G )理论:
R G G C =≥ (4)
其中R 是裂纹扩展的阻力,C G 是G 达到R 时的临界值。

该理论认为材料对
裂纹扩展的阻力等于弹性表面能与伴随裂纹产生的塑性应变功之和,成功的
解释了裂纹扩展现象,为线弹性断裂力学的发展奠定了基础。

(2)裂纹尖端应力场理论:
C K K = (5)
应力强度因子K 不仅仅适用于稳定的裂纹扩展,还适用于应力腐蚀和疲劳之
类的亚裂纹扩展。

(3)J 积分理论:
ds x
u T
wdy J ⎰Γ∂∂-= (6) 是一个回路积分,由Rice 提出,J 积分是一个与积分路径无关的常数,J 积分和 应力强度因子K 一样反映裂纹尖端的力学特性的。

1.4 本章内容小结
本章内容详细介绍了扩展有限元方法(XFEM )的基本理论和内容,包括单位分
解法、水平集法、引入节点扩展函数表征扩展有限元方法描述裂纹的位移模式,
以及扩展有限元方法基本特点和它的发展,以及与有限元方法模拟裂纹的对比,
最后还简要的介绍了断裂力学中裂纹的相关基本理论。

相关文档
最新文档