广义有限元方法

✍ dations ◷ 2025-08-28 01:05:54 #连续介质力学,数值微分方程,偏微分方程

广义有限元方法(英语:Generalized Finite Element Method,GFEM),是基于传统有限元方法(FEM)的一种数值分析方法,能够更加优雅、准确地解决材料交界面和断裂力学等非连续问题。这类问题在划分单元网格时,应用传统有限元方法往往需要使网格边界与非连续界面尽可能重合,带来极大的网格依赖(mesh-dependence)。对于复杂边界问题及形状变化问题(例如拓扑优化(英语:Topology optimization)和裂纹扩展等),传统有限元需要花费大量的时间建立单元网格以及重建网格,同时还需要引入计算量处理裂尖端(crack tips)来避免奇异刚性矩阵。广义有限元正是在这样的背景下诞生,它对传统有限元的形函数(shape functions)进行扩展,在包含非连续的单元网格内使用强化方程(enrichment functions),来准确模拟非连续的物理场,这种方法又可称为扩展有限元(extended finite element method,XFEM)。

非连续性,在有限元方法中的定义,指物理场量或其梯度在极小的区间内发生了极大的改变,其中“极小的区间”是相对于观察定义域而言,“物理场量“包括位移场(displacement field)、温度场(temperature field)、势场(potential field)等。物理场量的非连续往往被称为强非连续(strong discontinuities),其梯度的非连续被称为弱非连续(weak discontinuities),前者的表现方式包括断裂力学中裂缝两边的位移差,后者可举例为复合材料受力时在材料交界面的位移梯度差或是位移转折(kink)。

非连续问题的建模方法,可从基本原理上分为两类。第一类方法依赖于多项式近似空间(polynomaial approximation spaces)和单元网格与非连续界面的重合程度。其中最高效的多项式近似空间构造方法即为传统有限元方法,同时还有一些无网格法(英语:meshfree method)给多项式近似空间的构造带来不同的可行性。然而该类方法在处理复杂边界问题时,需要花费大量的时间建立单元网格,对于形状变化问题(例如拓扑优化和裂纹扩展等),求解过程中重建网格不可避免,对于断裂力学还需要引入计算量处理裂尖端来避免奇异刚性矩阵,给求解过程带来极大不便,求解问题被严重限制。即便如此,在广义有限元未被推出时代,传统有限元被广泛用在分析非连续问题上,尤其在土木工程中对河床、片岩、叶理、节理、裂纹对管道及地基等基础设施建造的影响分析中起到重要作用。

随着计算机硬件的提升,工程学上对力学建模的精确程度要求逐渐提高,又由于传统有限元分析非连续问题带来的种种复杂性,广义有限元应运而生。这就属于非连续问题的第二种建模方法——强化有限元多项式近似空间,使得非连续的建模可以独立于单元网格。这种方法在传统有限元的公式中,加入一项特殊的函数,称为强化项(enrichment),该强化项包含解中非连续性的信息,即物理场或其梯度的在非连续界面两端的差值,详细原理参见下一章节。该概念基于有限元单位分解法(英语:partition of unity)(partition of unity method),最初于1995年J.M. Melenk的博士论文中被提出,随后Melenk与他的博士导师、有限元领域的地基人物Ivo Babuška(英语:Ivo Babuška)于1996年整理发表。同一时期,另一所学校的博士生Duarte也在毕业论文中发表了相似理论,称为hp云方法(hp cloud method)。随后,广义有限元在断裂力学中的应用立即被深入研究,由Belytschko及其同事延伸成扩展有限元方法,由于最终被证明广义有限元与扩展有限元其实是同样的一种方法,该两个名词可互换使用。广义有限元的重要特性即在于解决了传统有限元中求解过程与单元网格的依赖性,网格的划分可以更加自由,同时,诸如裂纹扩展、拓扑优化、凝固过程(solidification)等涉及到非连续界面不断发生改变的复杂的问题,该方法都能更有效率、更准确地计算。

在广义有限元的基本公式中,第一项与传统有限元完全一致,第二项是其对传统有限元近似空间的扩展,即强化项,此处为了简化,只使用了一项强化项作为说明。

其中 I {\displaystyle I} 为定义域内所有节点的集合; N i {\displaystyle N_{i}} N i {\displaystyle N_{i}^{*}} 是传统有限元中的形函数,他们大部分情况下会相等(这里区别表示是由于,在后期一些高阶广义有限元的发展中,一些研究者曾尝试在这两项中使用不同阶数的形函数);系数 u i {\displaystyle u_{i}} 是标准有限元的节点所对应的自由度(degrees of freedom); a i {\displaystyle a_{i}} 是强化项引入的额外节点自由度,对应强化项作用的强化节点(enriched nodes),即是包含非连续的单元网格的节点。


ψ ( x ) {\displaystyle \psi (x)} 作为该公式中最重要的部分,是结合对解当中非连续性的预了解而建立的强化函数(enrichment function)。 ψ ( x ) {\displaystyle \psi (x)} N i ( x ) {\displaystyle N_{i}^{*}(x)} 的乘积,即运用了传统有限元拉格朗日形函数的单位分解性质(partition of unity property):

使 ψ ( x ) {\displaystyle \psi (x)} 仅仅作用在包含非连续界面的网格单元中。


下面解释最初的广义有限元概念中,使用的强化函数 ψ ( x ) {\displaystyle \psi (x)}

其中, ϕ ( x ) {\displaystyle \phi (x)} 被称为Level set函数,用于描述非连续界面的位置:

Ω {\displaystyle \Omega ^{-}} Ω + {\displaystyle \Omega ^{+}} 是定义域 Ω {\displaystyle \Omega } 被非连续界面分割出的子域。当中 ϕ d ( x ) {\displaystyle \phi _{d}(x)} 是每个点 x {\displaystyle x} 到非连续界面 Γ {\displaystyle \Gamma } 的最短距离:


对于一个一纬杆的弱非连续问题,如果使用1个单元,并给每个节点添加两个强化函数,用上述基本公式构建的广义有限元近似空间,其解 u h {\displaystyle u^{h}} 可用矩阵形式表示为:


虽然广义有限元的网格独立性很大程度上提高了非连续的建模,可它基本公示中的强化项却引入了新的问题。

Abaqus是采用了广义有限元算法的商业有限元软件。

相关

  • 干支纪年法干支是天干与地支的合称,由两者经一定的组合方式搭配成六十对,为一个周期,循环往复,称为六十甲子或六十花甲。古代中国用以记录年、月、日、时。日本、朝鲜、越南、琉球等汉字文
  • 主楼南开大学主楼1963年建成,2003年又进行了大规模的装修和翻新。主楼的设计沿中轴对称,气势雄伟,广场中间巍然挺立的周恩来总理雕像高大而伟岸。 2009年3月,入选天津十大标志性建筑
  • 臧伯平臧伯平(1913年6月-2005年10月4日),出生于河北省唐县,教育家,曾于1978年2月至1979年1月担任南开大学校长。臧伯平出生于1913年6月,曾就读于保定第二师范、北京民国大学国文系,早年即
  • 悉达医学悉达医学(泰米尔语:சித்த மருத்துவம் 或 தமிழ் மருத்துவம்,英语:Siddha medicine),起源自印度南部泰米尔纳德邦的传统医学,为古印度医学三大传统之
  • 蛛网纹蛛网纹(Spider Veins)又名蜘蛛网状静脉曲张,是指表肤底层的静脉血管扩张,导致血液沉积并出现弯弯曲曲的情况。令皮肤表面呈现出底层静脉的颜色,如红、蓝、绿、紫等颜色,外观好像蜘
  • NHK鸟取放送局NHK鸟取放送局,是日本放送协会位于鸟取县鸟取市的地方放送局,也是负责主管当地事务的放送局。
  • 约翰·赫顿·爱德华兹约翰·赫顿·爱德华兹(英语:John Hilton Edwards,1928年3月26日-2007年10月11日)是英国医药遗传学家,曾经在1946年至1949年间攻读剑桥大学药剂学系及动物学系,之后转读密德塞斯医院
  • 米罗斯拉瓦乡 (雅西县)坐标:47°09′N 27°31′E / 47.150°N 27.517°E / 47.150; 27.517米罗斯拉瓦乡(罗马尼亚语:Comuna Miroslava, Iași),是罗马尼亚的乡份,位于该国东北部,由雅西县负责管辖,面积83
  • 佩内洛普·菲兹杰拉德佩内洛普·菲兹杰拉德 (1916年12月17日–2000年4月28日) 布克奖获奖者,是英国的小说家、诗人,作家和传记作者。 2008年, 泰晤士报称她为"自1945年以来的50个最伟大的英国作家"
  • 王平 (1955年)王平(1955年4月-),男,陕西清涧人,中国人民解放军中将。早年长期任职于沈阳军区陆军第三十九集团军,后曾任陆军第十六集团军政治部主任、副政委。2006年12月,任南京军区陆军第一集团