有限元方法的发展及应用文档格式.docx

有限元方法的发展及应用文档格式.docx

《有限元方法的发展及应用文档格式.docx》由会员分享，可在线阅读，更多相关《有限元方法的发展及应用文档格式.docx（13页珍藏版）》请在冰点文库上搜索。

同时，为了模拟裂纹的动态扩展过程，也需要不断地进行网格重构。

（3）在处理夹杂问题时，单元的边须位于夹杂与基体的界面处，即使对于网格自动化程度很高的二维问题这也很不容易，而三维问题则更复杂。

1.3有限元法的派生

有限元法作为数值方法中的基础方法，有其一定的使用范围，也由于一定的弊端决定了其不完全通用性。

在有限元方法基础上，发展出有其特殊使用范围的更精准的派生数值方法，下面介绍几种重要的数值方法。

1.3.1有限差分法

有限差分法（FDM，FiniteDifferenceMethod）已经发展的一些近似数值分析方法中，最初常用的是有限差分法，它可以处理一些相当困难的问题。

但对于几何形状复杂的边界条件，其解的精度受到限制，甚至发生困难。

作为60年代最重要的科技成就之一的有单元法。

在理论和工程应用上都得到迅速发展，几乎所有用经典力学解析方法难以解决的工程力学问题郁可以用有限元方法求解。

它将连续的求解域离散为一组有限个单元的组合体，解析地模拟或逼近求解区域。

由于单元能按各种不同的联结方式组合在一起，且单元本身又可有不同的几何形状，因此可以适应几何形状复杂的求解域。

有限元的另一特点是利用每一单元内假设的近似函数来表示全求解区域上待求的未知场函数。

单元内的近似函数由未知场函数在各个单元结点上数值以及插值函数表达，这就使未知场函数的结点值成为新的未知量，把一个连续的无限自由度问题变成离散的有限自由度问题，只要结点来知量解出，便可以确定单元组合体上的场函数。

随着单元数目的增加，近似解收敛于精确解。

但是有限元方法常常需要很大的存贮容量，甚至大得无法计算；

由于相邻界面上只能位移协调，对于奇异性问题（应力出现间断）的处理比较麻烦。

这是有限单元法的不足之处。

1.3.2边界元法

边界元法（BEM，BoundaryElementMethod）是在有限元法之后发展起来的一种较精确有效的工程数值分析方法。

与有限元法在连续体域内划分单元的基本思想不同，边界元法是在定义域的边界上划分单元，用满足控制方程的函数去逼近边界条件，通过对边界分元插值离散，化为代数方程组求解。

降低了问题的维数，可用较简单的单元准确地模拟边界形状，利用微分算子的解析的基本解作为边界积分方程的核函数，而具有解析与数值相结合的特点，通常具有较高的精度。

边界元法的主要缺点是它的应用范围以存在相应微分算子的基本解为前提，对于非均匀介质等问题难以应用，故其适用范围远不如有限元法广泛，而且通常由它建立的求解代数方程组的系数阵是非对称满阵，对解题规模产生较大限制。

上述两种数值方法的主要区别在于，边界元法是“边界”方法，而有限元法是“区域”方法，但都是针对连续介质而言，只能获得某一荷载或边界条件下的稳定解。

对于节理裂隙发育的岩体或颗粒散体的处理则要麻烦得多，更无法进行大变形、分离、回转及塌落过程的模拟。

这就使得人们去探索和寻求适合模拟节理岩体和颗粒散体运动变形特性的有效数值方法。

1.3.3离散元法

离散元法（DEM，DistinctElementMethod）是由CundallPA（1971）首先提出并应用于岩土体稳定性分析的一种数值分析方法。

它是一种动态的数值分析方法,可以用来模拟边坡岩体的非均质、不连续和大变形等特点,因而,也就成为目前较为流行的一种岩土体稳定性分析数值方法。

该方法在进行计算时,首先将边坡岩体划分为若干刚性块体（目前已可以考虑块体的弹性变形）,以牛顿第二运动定律为基础,结合不同本构关系,考虑块体受力后的运动及由此导致的受力状态和块体运动随时间的变化。

它允许块体间发生平动、转动,甚至脱离母体下落,结合CAD技术可以在计算机上形象地反应出边坡岩体中的应力场、位移及速度等力学参量的全程变化。

该方法对块状结构、层状破裂或一般碎裂结构岩体比较适合。

1.3.4广义有限元法

广义有限元方法（GFEM，GeneralizedFiniteMethod）是常规有限元方法在思想上的延伸，它基于单位分解方法,通过在结点处引入广义自由度，对结点自由度进行再次插值,从而提高有限元方法的逼近精度，或满足对特定问题的特殊逼近要求。

基于广义有限元方法对单元形状函数构造理论的深入研究，具有任意内部特征（空洞、夹杂、裂纹等）及外部特征（凹角、角点、棱边等）的复杂问题,都将在简单、且与区域无关的有限元网格上加以求解。

1.3.5扩展有限元法

扩展有限元（XFEM，ExtendedFiniteElementMethod）是在标准有限元方法的框架下，提出来的一种用于解决裂纹、孔洞、夹杂等间断问题的数值方法。

在有限元的近似函数中，增加能反映待求问题间断特性的附加函数项，采用水平集方法（LSM）描述间断面的几何特性及其移动规律。

扩展有限元方法与标准有限元方法相比，具有计算精度高、勿需网格重构等特点。

2有限元法的发展

有限元法是R.Courant于1943年首先提出的。

自从提出有限元概念以来,有限元理论及其应用得到了迅速发展。

过去不能解决或能解决但求解精度不高的问题，都得到了新的解决方案。

传统的FEM假设：

分析域是无限的；

材料是同质的，甚至在大部分的分析中认为材料是各向同性的;

对边界条件简化处理。

但实际问题往往是分析域有限、材料各向异性或边界条件难以确定等。

为解决这类问题,美国学者提出用GFEM（Gener-alizedFiniteElementMethod）解决分析域内含有大量孔洞特征的问题；

比利时学者提出用HSM（theHybridmetisSingularelementofMembraneplate）解决实际开裂问题。

在FEM应用领域不断扩展、求解精度不断提高的同时，FEM也从分析比较向优化设计方向发展。

印度Mahanty博士用ANSYS对拖拉机前桥进行优化设计，结果不但降低了约40%的前桥自重，还避免了在制造过程中的大量焊接工艺，降低了生产成本。

FEM在国内的应用也十分广泛。

自从我国成功开发了国内第一个通用有限元程序系统JIGFEX后，有限元法渗透到工程分析的各个领域中，从大型的三峡工程到微米级器件都采用FEM进行分析，在我国经济发展中拥有广阔的发展前景。

目前在进行大型复杂工程结构中的物理场分析时，为了估计并控制误差,常用基于后验误差估计的自适应有限元法。

基于后处理法计算误差,与传统算法不同，将网格自适应过程分成均匀化和变密度化2个迭代过程。

在均匀化迭代过程中，采用均匀网格尺寸对整体区域进行网格划分，以便得到一个合适的起始均匀网格；

在变密度化迭代过程中只进行网格的细化操作，并充分利用上一次迭代的结果，在单元所在的曲边三角形区域内部进行局部网格细化，保证了全局网格尺寸分布的合理性，使得不同尺寸的网格能光滑衔接，从而提高网格质量。

整个方案简单易行，稳定可靠，数次迭代即可快速收敛,生成的网格布局合理，质量高。

2.1有限元法的国内外研究现状

FEM作为求解数学物理问题的一种数值方法，已经历了50余年的发展。

20世纪50年代，它作为处理固体力学问题的方法出现。

1943年，Courant第一次提出单元概念。

1945~1955年,Argyris等人在结构矩阵分析方面取得了很大进展。

1956年，Turner、Clough等人把刚架位移法的思路推广应用于弹性力学平面问题。

1960年,Clough首先把解决弹性力学平面问题的方法称为“有限元法”,并描绘为“有限元法=RayleighRitz法+分片函数”。

几乎与此同时,我国数学家冯康也独立提出了类似方法。

FEM理论研究的重大进展，引起了数学界的高度重视。

自20世纪60年代以来,人们加强了对FEM数学基础的研究。

如大型线性方程组和特征值问题的数值方法、离散误差分析、解的收敛性和稳定性等。

FEM理论研究成果为其应用奠定了基础,计算机技术的发展为其提供了条件。

20世纪70年代以来，相继出现了一些通用的有限元分析（FEA:

FiniteElementAnalysis）系统,如SAP、ASKA、NASTRAN等，这些FEA系统可进行航空航天领域的结构强度、刚度分析,从而推动了FEM在工程中的实际应用。

20世纪80年代以来，随着工程工作站的出现和广泛应用,原来运行于大中型机上的FEA系统得以在其上运行，同时也出现了一批通用的FEA系统,如ANSYS-PC、NISA,SUPERSAP等。

20世纪90年代以来,随着微机性能的显著提高，大批FEA系统纷纷向微机移植，出现了基于Windows的微机版FEA系统。

经过半个多世纪的发展,FEM已从弹性力学平面问题扩展到空间问题、板壳问题;

从静力问题扩展到动力问题、稳定问题和波动问题;

从线性问题扩展到非线性问题；

从固体力学领域扩展到流体力学、传热学、电磁学等其他连续介质领域；

从单一物理场计算扩展到多物理场的耦合计算。

它经历了从低级到高级、从简单到复杂的发展过程,目前已成为工程计算最有效的办法之一。

2.2有限元法的网格化分发展

作为有限元走向工程应用枢纽的有限元网格划分，是有限元法的一个非常重要的研究领域，经历了40多年的发展历程。

有限元网格划分算法研究中的某些难点问题始终未能得到真正意义上的解决，它们的解决对工程问题具有重要的现实价值和理论意义。

有限元分析的基本过程可分为三个阶段：

有限元模型的建立（即前处理）、有限元解算、结果处理和评定（即后处理）。

根据经验,有限元分析各阶段所用的时间为：

40%-45%用于模型的前处理，50%-55%用于后处理，而分析计算只占5%左右；

更有指出有限元建模占有限元分析一半以上的工作量,甚至高达80%。

因此，有限元分析的前后处理一直都是有限元分析的瓶颈问题，严重地阻碍着有限元分析技术的应用和发展。

许多学者对有限元网格生成方法近30年的研究进行了概括和总结，对某些重要分支领域的研究进展方面也做出了贡献。

近年来,有限元网格生成方法研究有两个显著特点：

（1）经历了一个进化过程，一些方法的研究与应用出现停滞，而另外一些方法在不断地深入、完善和发展，成为适应性强、应用范围广泛的通用方法；

（2）领域和主题在不断扩展和深入，研究重点由二维平面问题转移到三维曲面和三维实体问题，从三角形、四面体网格自动生成转移到四边形、六面体网格自动生成。

3有限元法的应用

有限元法最初应用在求解结构的平面问题上,发展至今，已由二维问题扩展到三维问题、板壳问题，由静力学问题扩展到动力学问题、稳定性问题，由结构力学扩展到流体力学、电磁学、传热学等学科，由线性问题扩展到非线性问题，由弹性材料扩展到弹塑性、塑性、粘弹性、粘塑性和复合材料，从航空技术领域扩展到航天、土木建筑、机械制造、水利工程、造船、电子技术及原子能等，由单一物理场的求解扩展到多物理场的耦合，其应用的深度和广度都得到了极大的拓展。

3.1有限元法的应用过程

FEM应用于实际问题须经历以下过程，如图1所示。

应用

理论研究

否

是

计算结果

有限元建模

系统开发

图1.FEM的应用过程

（1）问题的数学描述。

对问题客观规律的数学描述（通常是微分方程及边界条件）是建立有限元方程的前提。

单元特性矩阵和整体有限元方程都是基于数学模型建立的。

常见的弹性力学基本方程、运动方程、热传导方程等都是对客观现象的数学描述。

（2）有限元方程的建立。

利用变分原理,通过离散、单元分析、整体分析等过程,建立数学模型的有限元方程，它通常是一组易于用数值方法求解的代数方程。

（3）算法研究。

有限元方程的计算量庞大，须有有效的算法来保证计算效率和精度，同时考虑对计算条件的要求。

如求解大型线性方程组的带宽法、波前法，求解大型特征值问题的分块Lanczos法等。

（4）程序开发。

数值计算依赖于计算机,因此求解算法需用相应的计算程序来实现。

（5）有限元建模。

对应于FEA系统的前处理（Pre-processing）。

它为数值计算提供所有原始输入数据（节点数据、单元数据和边界条件数据）。

因为模型形式直接决定计算精度和规模，且建模所需时间约占整个FEA的70%左右，所以建模质量和效率是FEA的关键。

图2列出了有限元建模中的关键技术。

图2有限元建模的关键技术

（6）数值计算。

对应于FEA系统的计算（Solving）。

它由一系列计算程序组成,计算程序又称求解器（solver）。

每个求解器完成特定类型的计算。

因此求解器越多,系统功能越强。

（7）结果处理。

对应于FEA系统的后处理（Post-processing）。

它对计算结果进行处理、显示、运算和列表等。

若按照

（1）~（7）过程,问题得以解决，则FEM应用结束;

反之,则需根据求解结果提出改进方案,循环执行（5）~（7）过程,直至问题解决或得到最佳设计。

对于一个全新的问题,必须从第一步开始。

而对已知的问题，可从第（5）步开始，即直接利用已有的FEA系统,建立有限元模型。

在实际应用中，绝大多数问题都属于第二类问题。

3.2有限元法的应用领域

FEM最早应用于固体力学领域,但由于其解决问题的有效性和实用性，很快推广应用于温度场、电磁场、流场、声场等连续介质领域。

目前FEM的应用领域主要包括：

（1）静力分析。

包括线性非线性静力分析。

线性静力分析研究线弹性结构的变形和应力，它是工程结构分析和设计中最基本的方法。

非线性结构静力分析主要研究外载作用下引起的非线性响应，其中非线性来源主要是材料非线性、几何非线性和边界条件非线性3大类。

（2）动力分析。

主要包括以下分析类型：

1）模态分析。

用于求解多自由度系统的模态参数。

为计算得到的计算机主板的前三阶振型。

2）瞬态响应分析。

求解在时域内结构承受随时间变化的载荷和速度作用时的动力响应。

3）简谐响应分析。

对简谐激励结构在其平衡位置的振动进行分析。

4）频谱响应分析和随机振动分析。

用于分析结构受已知频率激励时的最大响应。

5）屈曲和失稳分析。

分析考察结构的极限承载能力，研究结构总体或局部的稳定性，获得结构失稳形态和失稳路径。

6）自动接触分析。

用于接触边界定义和摩擦分析。

（3）失效和破坏分析。

包括断裂分析（线弹性断裂分析和弹塑性断裂分析）、裂纹萌生与扩展分析、跌落分析和疲劳失效分析。

（4）热传导分析。

包括稳态热传导分析、瞬态热传导分析、热辐射、强迫对流及温度的耦合分析。

（5）电磁场分析。

它用于对电磁场中电感、电容、磁通量密度、涡流、电场分布、磁力线分布、能量损失等物理量进行分析。

（6）声场分析。

它用来研究在含有流体介质中声波的传播问题,或分析浸在流体中的固体结构的动态特性。

（7）研究流体速度、压强、密度变化规律和粘滞流体的运动规律及粘滞流体中运动物体所受阻力及其它热力学性质。

（8）耦合场分析。

考虑两种或两种以上物理场的交叉作用和相互影响（耦合）。

3.3有限元法的应用的热点和前景

随着FEM研究的深入，过去不能解决或能解决但求解精度不高的问题，都得到了新的解决方案。

材料是同质的,甚至在大部分的分析中认为材料是各向同性的;

为解决这类问题，美国的heofanisStrouboulis&

LinZhang等人提出用GFEM（GeneralizedFiniteElementMethod）解决分析域内含有大量孔洞特征的问题；

比利时的NguyenDangHung和越南的TranThanhNgoc提出用HSM（theHybridmetisSingularelementofMembraneplate）解决实际开裂问题（结构尺寸有限,形状任意,边界条件复杂,材料特性任意）。

传统的有限元断裂力学技术（thefiniteelementfracturemechanicstechniques）在解决零件中出现裂缝这类问题时，需要在曲线型裂纹前缘附近的区域细分网格。

这样无论是从网格生成的角度看还是从求解的角度看，都需要花费大量的时间。

而且在循环加载的情况下产生的次裂纹将会使分析变得更加复杂。

为此，美国的DanielSPipkinsay&

SatyaNAtlurib提出了FEAM（FiniteElementAlternatingMethod）。

该方法在求解应力集中因子时，可在不牺牲精度的情况下节省时间，用它分析具有椭圆裂纹或部分椭圆裂纹的结构是很有用的。

此外，西班牙的OnateE和波兰的RojekJ将DEM（DiscreteElementMethod）和FEM结合解决地质力学中的动态分析问题;

瑞典的BirgerssonF和英国的FinnvedenS针对FEM在频域中的应用提出了SFEM（SpectralFiniteElementMethod）。

印度Mahanty博士用ANSYS对拖拉机前桥进行优化设计，结果不但降低了约40%的前桥自重,还避免了在制造过程中的大量焊接工艺，降低了生产成本。

20世纪80年代我国大连理工大学工程力学研究所开发成功了国内第一个通用有限元程序系统JIGFEX，并在1983年开发出了它的微机版JIGFEX-W。

目前,FEM已渗透到工程分析的各个领域,从大型的三峡工程到微米级器件都采用FEM进行分析。

它在我国经济发展中拥有广阔的发展前景。

FEM的研究热点目前表现在两个方面:

超收敛应力计算和有限元模型修正技术。