扩展有限元简介
扩展有限元法(XFEM)及其应用12
• 前 言 • 固体力学中存在两类典型的不连续问题,一类是 因材料特性突变引起的弱不连续问题,这类问题 以双材料问题和夹杂问题为代表,其复杂性是由 物理界面处的应变不连续性引起的;另一类是因 物体内部几何突变引起的强不连续问题,这类问 题以裂纹问题为代表,其复杂性是由几何界面处 的位移不连续性和端部的奇异性引起的。物体内 部物理界面的脱粘或起裂,是上述两类问题的混 合,也属于这里所讨论的范围。另外,在复杂流 体、复杂传热、物质微结构演化等复杂问题中, 也存在许多不连续力学问题。
• 在h型中,多项式次数p是固定 的(一般p≦2),逼近通过减 小网格尺寸予以实现,恰当插 值能产生一个满足连续性条件 的良好逼近;在p型中,局部逼 近通过增加多项式次数实现。
• 如果PUM中局部逼近空间Vi选作多项 式空间,PUM就可看作是经典h和p型 的推广,以这种方式构造的PUM空间 在逼近特性上与CFEM空间非常相似。 网格设计(即确定需要进行网格重构 的面积)和每个单元上多项式的次数 对h- p有限元方法性能影响很大,但它 们取决于对解局部行为的掌握。
• 可见,单位分解概念保证了XFEM的收敛, 基于此XFEM的逼近空间中增加了与问题相 关的特定函数;水平集法是XFEM中常用的 确定内部界面位置和跟踪其扩展的数值技 术,任何内部界面可用它的零水平集函数 表示。本章拟在第二和第三部分对单位分 解法和水平集法分别进行简要介绍;在第 四和第五部分重点介绍XFEM的基本思想、 实施步骤和若干应用实例;最后,初步展 望该领域尚需进一步研究的问题。
扩展有限元法在岩石断裂力学中的应用_王洪建
Abstract
Rock crack propagation cannot be simulated with conventional finite element method because it is a
wk.baidu.com
classic discontinuity problem. Extended finite element method ( XFEM ) can implement the mutual independence of computational grid and discontinuity. So there is no need to remesh the grid model for moving discontinuity simulation. This paper introduces the basic principles of XFEM and universal criterion for rock fracture mechanics. It utilizes XFEM to simulate I mode and II mode fracture in threepoint bending test as well as the doublecrack shear fracture test. All the obtained results indicate that using XFEM to simulate rock fracture cannot be restrained by mesh generation and the cracks propagate randomly according to the actual stress distribution. The whole process of micro crack generation and evolution until the final total destruction of rock can be directly observed,which have
扩展有限元培训
综述
扩展有限元(XFEM)是在标准有限元 方法的框架下,提出来的一种用于解决裂 纹、孔洞、夹杂等间断问题的数值方法。
虚拟节点-叠加于初始真实节点上-用于表示 断裂单元的间断性,如图所示。当单元保 持完整时,每一个虚拟节点被约束于相应 的真实节点上。当单元被裂纹切开时,断 裂单元被分成两部分。每一部分均由部分 真实节点和虚拟节点组成(与裂纹方向有 关)。每一个虚拟节点不再绑定与与其对 应的真实节点上,并可以独立移动。
百度文库
创建Section并把它赋予给断裂体
赋予材料取向
装配模块
网格模块
网格控制为:Hex型 structured单元类型为C3D8R,设置 plate各边的网格种子为8,26,36。
分析步模块
分析步控制
Load模块 施加边界条件与不均衡载荷
设置XFEM型裂纹
结果输出
实例介绍
Part模块操作 建立裂纹实体与裂纹位置 裂纹实体为3D deformable solid extrusion 裂纹位置为3D deformable shell extrusion
材料设置 弹性参数为E=210GPa,泊松比为0.3; 最大主应力失效准则作为损伤起始的判据; 损伤演化选取基于能量的、线性软化的、 混合模式的指数损伤演化规律; Damage Stabilization Cohesive
扩展有限元的基本知识1
制作时间:2014.12.12
有限元在处理间断问题的缺陷
有限元采用的是连续性的位移近似函数,对于裂 纹类强间断问题,为获得足够的计算精度,需要 对网格进行足够的细分,计算量极大。 采用拉格朗日法求解裂纹动态扩展、流固耦合、 局部剪切等特大变形问题时,有限元网格可能会 发生严重扭曲,使计算精度急剧下降甚至计算无 法继续,因此,需要不断的进行网格重构、计算 量极大,同时,也为了模拟裂纹的动态扩展过程, 也需要不断的进行网格重构。
ห้องสมุดไป่ตู้
单元分解函数
扩展有限元近似函数的基础是单元分解法。单元分 解法使用一些以节点 xi 为中心的子域 i 来覆盖整 个求解区域即 n
i
i 1
在每个子域 i上定义一个仅在该子域内非 零的函数 i (x) ,并且它们满足单位分解条 件: n 则函数集 {i (x)}in1 称为属于开 (x) 1
扩展有限元的概念
扩展有限元(XFEM):是在标准有限元方法的 框架下,提出来的一种用于解决裂纹、孔洞、夹 杂等间断问题的数值方法。在有限元的近似函数 中,增加能反映待求问题特性的附加函数项,采 用水平集法(LSM)描述间断面的几何特征及其 移动规律。
扩展有限元的优点
• 计算精度高 • 勿需网格重构
扩展有限元的提出
1999年,美国西北大belytschko 研 究组提出的扩展有限元。借助于对 研究问题的已有认识,在满足单位 分解的前提下,在位移近似函数中 增加更能反映实际间断特征的函数 项(称为富集函数)提高了计算精 度。采用水平集法(LSM)或快速 推进法(FMM)描述间断界面,使 间断的描述独立于有限元网格,避 免了计算过程中的重构。
扩展有限元法
扩展有限元法
#### 1.一维稳定性
有限元法可以用来解决一维稳定性问题,利用一维稳定性的基本原理,通过对实际结构的有限元单元网格进行分形集中或分散处理,以及对单元类型和形式进行选择,有效获取和分析各种参数,如材料参数、载荷参数、面宽参数以及梁中无弯度形状参数,从而得出有限元反应的刚度、弯矩、变形和应力的计算结果,进而估计和预测一维结构的承载能力、抗屈曲、变形或塑性扭曲变形。
#### 2.刚构结构
利用有限元法,可以进行刚构结构的力学分析,可以知晓构件在由外力所作用下,构件各点的位移、变形,以及构件各点的应力、应变,同时由于刚构结构个别构件之间也会存在约束关系,故必需考虑构件之间的相互影响,利用约束条件完成构件的组合,以此来估计稳定的系统的性能。
扩展有限元方法及应用综述_郭历伦
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 的单位分解函数
[ 3~4]
N
。 对于任意函数 V( , 可得 x)
I=1
。 x) Fra Baidu bibliotek( x) =V( x) ∑ (
I
设函数 VI( 为函数 u( 在子域 ΩI 内的局部近似函数 , 则函数 u( 在求解域的全局近似可取为 x) x) x)
采用水平集方法后富集函数所采用的极坐标系由水平集函数定义见式14同时计算j积分时由水平集函数定义虚速度矢量14belytschko和blac以及mos等的研究均是基于准静态二维裂纹问题在提出xfem方法时并未涉及裂纹扩展因此仅采用有向距离函数就可以描述裂纹的几何特性无需求解式17sukumar等采用快速推进法研究了三维平面裂纹扩展问题并在后续工作中逐步研究了多个平面裂纹的扩展问题及三维非平面裂纹扩展问题stolarsko等10采用水平集方法研究二维裂纹扩展问题mos位错问题fig5twolevelsetfunctions1112将水平集方法应用于三维非平面裂纹问题中ventura等13提出了矢量水平集vectorlevelset14算法
扩展有限元的基本原理
1 引言
有限元法是科学研究和工程分析中使用最多的一种数值分析方法,其数学理论基础和误差估计理论都 十分成熟。然而,有限元采用连续函数作为形函数,对于处理像裂纹这样的不连续问题时,需要将裂纹面 设置为单元的边、裂尖设置为单元的结点、在裂尖附近不连续体的奇异场内进行高密度网格划分以及在模 拟裂纹扩展时需要不断的进行网格的重新划分,使得有限元程序计算相当复杂,且效率极低[1]。边界元法 研究裂纹扩展方面有较成熟的应用,但它不便于处理非线性、多介质等复杂问题。无单元法将整个求解域 离散为独立的结点,无需将结点连成单元,因而在裂纹扩展数值模拟中得到了广泛的应用[2,3]。但是现有的 各种无单元法存在缺少坚实的理论基础和严格的数学证明; 计算时间长、 效率低; 存在一些未确定的参数, 如插值域的大小,背景积分域的大小等;解决复杂的工程和科学问题的研究不够; 没有成熟的商业软件包, 限制了其实际应用和推广等不足。 1999 年,以美国西北大学 Belytschko 教授为代表的研究组首先提出用扩展有限元(XFEM)来解决不 连续问题[4]。XFEM 是基于单位分解的方法(PUM)对单元的形函数加以改进,从而考虑所研究问题的不 连续、奇异性和边界层等特性。XFEM 所使用的网格与结构内部几何或物理界面无关,从而克服了裂纹尖 端等高应力和变形集中区网格划分的困难,使得模拟裂纹生长也无需对网格进行重新划分。自 XFEM 问世 以来,在国际上得到了很快的发展和广泛的应用。 在最初的 XFEM 中,位移模式中加进的是裂尖渐进位移场函数的主要项,且围绕裂尖加强结点的相应 系数是相互独立的,这样加强位移场并不是真实的裂尖附近的渐进位移场,因此,局部位移场的精度仍不 能令人满意,应力强度因子必须经过后处理才能求出。文[5][6]提出了一种改进的 XFEM,提高了局部位移 场的精度,且不需要经过后处理就可以直接求出应力强度因子,从而为分析裂纹扩展提供了方便。 文中首先简要介绍了改进扩展有限元的基本原理,然后详细分析了采用该方法分析不连续力学问题需 要注意的问题,包括:积分方案、裂纹闭合时的接触问题和裂纹扩展的分析。
扩展有限元方法和裂纹扩展
扩展有限元方法和裂纹扩展
1.1 扩展有限元方法(XFEM )基本理论
1999年,美国Northwestern University 的Belytschko 和Black 领导的研究小
组提出了扩展有限元方法,为解决裂纹这类强不连续问题带来了曙光。他们正式
应用扩展有限元法(XFEM )这一专业术语是在2000年,截止到目前,扩展有
限元法(XFEM )成为我们解决强不连续力学问题的最有效的数值计算方法,也
成为计算断裂力学的重要分支。XFEM 在有限元的框架下进行求解,无需对构件
内部的物理界面进行网格划分,具有常规有限元方法的所有优点。它最明显的特
点是用已知的特征函数作为形函数来使传统有限元的位移得到逼近,进而克服了
在裂纹尖端和变形集中处进行高密度网络划分产生的困难,方便地模拟裂纹的任
意路径,而且计算精度和效率得到了显著的提高[6]。
扩展有限元方法是将已知解析解的特征函数作为插值函数增强传统有限元
的位移逼近,来使得单元内的真实位移特性得以体现,裂纹尖端和物理或几何界
面独立于有限元网格。XFEM 主要包括以下三部分内容:首先是不考虑构件的任
何内部细节,按照构件的几何外形尺寸生成有限元网格;其次,采用水平集方法
跟踪裂纹的实际位置;根据已知解,改进影响区域的单元的形函数,来反映裂纹
的扩展。最后通过引入不连续位移模式来表示不连续几何界面的演化。因为改进
的插值函数在单元内部具有单元分解的特性,其刚度矩阵的特点与常规有限元法
的刚度矩阵特性保持一致。单元分解法(Partition Of Unity Method)和水平集法
有限元断裂模拟(动图)
Vcct就是虚拟裂纹闭合法,通过假设裂纹闭合所需的能量来解决这类问题。
扩展有限元法(XFEM)
扩展有限元法是解偏微分方程的一种数值方法
XFEM例子
XFEM与删除单元
谢谢!
有限元断裂模拟(动图)
断裂
工程上模拟断裂的方法有
1、自适应remeshing(最古老的方法)
Remeshing就是不断调整裂纹尖端网格密度,只有通过增加局部自由度,才能表征 裂尖行为。
2、单元删除法
单元删除法,原理是通过设定一定的失效准则,当单元达到准则时单元删除。
单元删除法劣势是不太准确,很依赖网格尺寸,模拟效果不好,除非特别对待。 3、扩展有限元法(XFEM)
扩展有限元,是目前世界上最流行的断裂模拟方法,并且效果理想。对于裂纹来 说,裂纹萌生模拟起来是困难的,所以xfem中裂尖前端使用粘聚区模型来表征其 萌生,这样其萌生和起裂方向就和实验比较吻合。
4、intபைடு நூலகம்relement法(cohesive 单元或接触行为)
Interelement, 就是在常规单元间加入特殊的单元,比如cohesive单元来表征裂纹
基于扩展有限元的应力强度因子的位移外推法
基于扩展有限元的应力强度因子的位移外推
法
扩展有限元法(XFEM)是有限元分析方法的一种扩展,主要用于模拟复杂断裂过程中的位移和应力分布。XFEM可以用于计算应力强度因子(SIF),这是一个用于描述裂缝尖端应力状态的参数。然而,计算SIF需要大量的计算资源和时间,并且对于大尺度的问题,计算结果也可能不精确。为了解决这个问题,研究人员提出了一种称为位移外推法的技术,它能够准确快速地计算裂纹尖端的位移场,并从中得出SIF的值。
位移外推法的基本思想是通过在裂纹尖端附近进行位移场外推来估计裂纹尖端的位移场。外推的基本思路是利用裂纹周围的位移场数据,通过合适的插值方法,在尖端附近构建一个位移场模型。理论上,如果外推的位移场能够准确地描述尖端附近的位移场,则可以通过该位移场计算裂纹尖端处的SIF。
在位移外推法中,首先需要确定一个合适的位移外推区域,通常该区域是以裂纹尖端为中心的一个圆盘。接下来,需要确定外推区域内的位移场边界条件,这是通过裂纹附近的位移场数据进行插值得到的。常用的插值方法有径向基函数插值、距离加权插值等。然后,根据插值得到的位移场构建一个位移模型,并将其应用于计算裂纹尖端处的SIF。
位移外推法的优点是可以在较小的计算域内获得高精度的位移场数据,并且能够通过这些数据快速准确地计算SIF。然而,该方法也存在一些限制。首先,位移外推法的精度依赖于外推区域的选择和位移场数据的精度。如果外推区域选择不合适或者位移场数据不准确,得到的SIF 结果可能不可靠。其次,位移外推法需要较为复杂的插值方法和位移模型,这需要一定的数值计算和编程技巧。
扩展有限元法的研究
1.1.2 Cohesive 裂纹 对于一个静态 cohesive 裂纹,Zi 和 Belytschko[18]提出使用符号函数(sign function)来改进所有三结点和六结点的包含裂尖的三角形单元,其中,为了保证 cohesive 裂纹的光滑闭合,强制令裂尖的应力投影与材料强度相等。 由于不需要每个有限元划分都是均匀的,Alfaiate 等[19]引入位移跳跃来描述 半脆性材料中的裂纹。对于半脆性材料的 cohesive 裂纹扩展问题,Mariani 和 Perego[20]在 XFEM 中引入高阶位移不连续性,研究了准脆性材料中准静态粘着 裂纹的传播, 其有效性通过 I 型和混合型数值试验给出了评定。 Patzak 和 Jirasek[21] 提出一种新的技术,用于计算半脆性材料狭窄损伤区域的局部应变,其中软化所 引起的局部化将由一个非局部损伤模型来描述, XFEM 使用与精确局部模式相 近的规则 Heaviside 函数来改进,使得位移逼近函数有更好的适应性。 Simone 等[22]提出一个用来计算描述含有不连续缺陷的应变软化连续体的计 算框架,并在连续体的本构方程中引入梯度项。在缺陷周围的区域使用基于 FE 形函数的单位分解的不连续插值函数来改进, 内部不连续面的边界条件也不使用 有限元划分引入,这就避免了规则化连续体模型的不真实损伤扩展。这种联合模 型允许分析整个失效过程,从细观裂纹的扩散到宏观裂纹的局部出现都能够分 析。De Borst 等[23]证明在混凝土断裂中所使用的裂纹模型,也可以通过合适的
基于ABAQUS平台的扩展有限元方法
基于ABAQUS平台的扩展有限元方法
扩展有限元方法(XFEM)是一种能够有效处理裂纹、接触、损伤等大
变形、大变位问题的计算方法。该方法扩展了传统有限元方法(FEM),
使其能够更准确地模拟物体的断裂行为。ABAQUS是一款常用的有限元分
析软件,提供了XFEM功能,可以在其平台上进行XFEM分析。
XFEM的主要思想是在有限元网格中引入额外的自由度,这些自由度
可以用来描述物体内部的裂纹、接触等特征。通过在网格中引入额外的基
函数,XFEM能够通过有限元分析获取到界面上的开裂和断裂行为,从而
更准确地预测物体的破坏。
XFEM在ABAQUS平台上的应用主要包括以下几个方面。
1.裂纹模拟:XFEM能够准确地模拟裂纹的扩展行为。在ABAQUS中,
用户可以通过定义裂纹路径和裂纹扩展准则,来模拟裂纹在不同加载条件
下的扩展过程。同时,用户还可以对裂纹的形状、长度、位置等进行控制,以得到更准确的结果。
2.接触分析:XFEM可以模拟接触问题,包括刚性接触和非线性接触。在ABAQUS中,用户可以通过定义接触面和接触行为,来模拟物体之间的
接触行为。XFEM能够考虑接触面的开裂和闭合,从而更准确地模拟接触
问题。
3.损伤模拟:XFEM可以模拟材料的损伤行为,包括塑性、弹塑性和
弹性损伤。在ABAQUS中,用户可以通过定义损伤模型和损伤准则,来模
拟材料的损伤行为。XFEM能够考虑材料中的裂纹行为,从而更准确地模
拟损伤问题。
4.多物理场耦合:XFEM可以模拟多个物理场的耦合问题,如固体力学和热传导、固体力学和流体力学等。在ABAQUS中,用户可以通过定义不同物理场的边界条件和耦合关系,来模拟多物理场耦合问题。XFEM能够考虑多物理场之间的相互作用,从而更准确地模拟多物理场问题。
扩展有限元简介
• 数值方法,如有限元、边界元、无网格法等,一直是处理 不连续问题的主要途径
1.1 有限元裂纹扩展模拟
• 裂纹扩展有限元模拟研究涉及三个问题:理论基础、扩展控 制参量及模拟方法。 • 理论基础直接影响有限元方程构成和具体实施的难易程度。
1.2Ⅰ、Ⅱ型复合裂纹扩展模拟
• 根据能量平衡理论,裂纹扩展的条件是裂纹体在裂纹扩展过程 中能释放出足够的弹性能,以供表面能增加的需要。定义裂纹 扩展单位长度(厚度为单位厚) 所释放出来的能量为裂纹扩展能 量释放率, 记为G。由于能量释放率G能描述在任何应力状况下 裂纹扩展时的基本量参数,所以在探讨断裂力学中应用广泛。 虚裂纹扩展法是一种能量方法, 是通过虚裂纹扩展这一技术求出 裂纹的能量释放率G, 该方法也叫能量释放率法。 • 在使用ANSYS软件建模时,通常将2D模型的裂纹尖端作为裂纹 的边缘,一般2D模型模拟裂纹扩展因试件长度和宽度方向的尺 寸远远大于厚度方向的尺寸,且边界条件满足平面应力问题的 条件,故可简化为平面应力问题求解。ANSYS软件采用KSCON命 令通过自动产生奇异单元来划分网格,并用KCALC命令来计算并 读取其应力强度因子值。
������, ������ = 0
(5)
水平集对裂纹的描述
• 定义裂纹面������(������) ⊂ ������2 可以被水平集曲面������(������, ������) ⊂ ������2 × ������ → ������确切地表达,其中 S ������ = ������ ∈ ������2 : ������ ������, ������ = 0 • 移动界面由������ ������, ������ 的演化方程确定������ ������, 0 给定 ������������ + ������ ������������ = 0 • F(x,t)为界面上点x在界面外法线方向的速度;水平集函数f 为符号距离函数 ������ ������, ������ = ± min ������ − ������������
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/1437402244
1.1概要
XFEM即扩展有限元方法,它在标准有限元框架内研究问题,保留了有限元方法的所有优点。扩展有限元法与有限元法最根本的区别在于所使用的网格与结构内部的几何或物理界面无关,从而克服了在诸如裂纹尖端等高应力和变形集中区进行高密度网格划分所带来的困难,在模拟裂纹扩展时也无需对网格进行重新划分。
有限元和扩展有限元的关系
有限元和扩展有限元的关系
有限元方法和扩展有限元方法都是求解力学问题的数值计算方法。在一些特定的工程和科学应用中,需要对非常复杂和不规则的几何形态的结构进行建模和分析,有限元方法和扩展有限元方法就可以很好地解决这些问题。
有限元方法将复杂的结构分解成许多小的、简单的有限元单元,然后通过对每个单元的力学特性进行分析,在整个结构上求解出力学问题。而扩展有限元方法则是对有限元方法的扩展,可以在有限元网格上增加额外的局部自由度,使其更加适用于分析复杂的断裂和损伤行为的问题。
因此,有限元方法可以看作是扩展有限元方法的基础。扩展有限元方法通过增加局部自由度,可以提高在裂纹或者材料失效的情况下的计算精度。它是在有限元方法的基础上发展而来的一种方法,除了更容易模拟复杂裂纹模式的形成和扩展外,还能够模拟在较大应变下的材料本构行为。
总之,有限元方法和扩展有限元方法都是有效的数值计算工具,可以帮助工程师和科学家对复杂的结构进行建模和分析,从而预测其力学行为。
扩展有限元matlab
扩展有限元matlab
有限元方法是一种常用的数值计算方法,它广泛应用于工程学、
物理学、生物学等领域中的复杂问题的建模和求解。在有限元分析中,使用离散化的方法将复杂的连续域分解成独立的小元素,并对每个小
元素进行有限差分计算。因此,有限元方法在工程学领域中,如结构
力学、热力学和流体力学等方面应用广泛。
Matlab是一个流行的数学计算环境,广泛应用于科学和工程领域。Matlab提供了一个强大的工具箱,支持有限元分析,称为PDE工具箱。PDE工具箱允许用户设置和求解的偏微分方程,并提供了大量的有限元网格,以允许用户执行自定义模拟。
为了扩展有限元Matlab工具,许多软件包已经开发出来以支持
特定领域的问题,如结构分析、流体力学、声学和电磁问题等。这些
扩展都提供了额外的工具,并通过提供针对特定应用的专业算法来改
进有限元分析的质量。
在使用Matlab进行有限元分析时,应遵循一般的数值方法要求。一些常见的数值问题包括稳定性问题、收敛性问题和模型不精确问题等。在处理这些问题时,必须仔细选择网格分辨率、求解器和数值方
法等参数。此外,还必须进行适当的验证和验证有限元模型。
因此,扩展有限元Matlab工具可以帮助使用者更好地应对这些
问题,提高有限元方法的可靠性和效率。同时,随着科学计算和工程
学领域中对模拟和优化的需求不断增长,有限元Matlab工具的发展和
改进将有助于更好地解决日益复杂的工程问题。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
扩展有限元
有限元是将一个物理实体模型离散成一组有限的相互连接的单元组合体, 该方法在考虑物体内部存在缺陷时间,单元边界与几何界面一致,会造成局部网格加密,其余区域稀疏的非均匀网格分布,在网格单元中最小的尺寸会增加计算成本,再者裂纹的扩展路径必须预先给定只能沿着单元边界发展。 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 裂纹尖端的渐进函数并不局限于各向同性弹性材料的裂纹建模。可用于弹塑性指数硬化材料,不同的裂纹尖端渐进函数的形式与裂纹位置、非线性材料变形程度有关。
扩展有限元方法(XFEM)可以缓解裂纹面网格划分带来的缺点,其基于整体划分的概念使扩展函数方便地插入到有限元中。间断性可以通过与额外自由度相关联的扩展函数来确定。同时扩展有限元方法保留了有限元框架及一些特性,如刚度矩阵的稀疏性及对称性等。因此以有限元为基本框架,主要针对不连续问题进行研究,相对于传统有限元方法,扩展有限元法提高了描述复杂位移场的能力,增加了对于演化的非连续边界进行跟踪的灵活性,避免了网格重划分,当考虑裂纹扩展时其计算工作量和准确性要明显优于有限元方法。