有限元分析(FEA)软件

有限元分析(FEA)软件

   物理定律、偏微分方程和数值建模   有限元法   有限元分析软件 

有限元分析软件能够带来什么优势?

有限元分析(Finite Element Analysis,FEA)软件可以帮助企业减少在产品或者流程的设计、优化或控制环节中,原型测试的原型数量和测试次数。对于企业和研究机构来说,有限元仿真分析带来的不仅仅是成本的降低,更重要的是在激烈的市场竞争中赢得优势,为研发投入带来了更大的回报。正因如此,近年来,越来越多的企业将更多的研发资源投入到有限元分析中。

一旦建立了能够准确预测真实物理参数的有限元分析模型,工程师们就可以借助它来加强对物理现象的理解和认识,以大幅改进产品或过程的设计和运行。在此基础上,优化算法和自动控制的应用,可以进一步实现仅凭直觉完全无法达到的设计改进。目前的有限元分析软件大多已包含自动控制功能,并将这些功能嵌入数学和数值模型中,而优化算法也通常包含在求解过程中,下文将会详细介绍。

高保真模型的引入,可以帮助工程师们加深理解、激发灵感,带来全新的设计和方案。正是因为这个原因,对于面临着激烈竞争的企业来说,有限元分析是研发部门不可或缺的工具。近年来,有限元分析软件的使用越来越广泛,已经从大型企业以及工程师的培养机构,扩展到各行各业的中小型企业和涉及各个学科的研究型机构中。

深入解读有限元分析软件

基于数学模型表示的物理定律构成了有限元分析软件的基础。对于有限元分析来说,这些定律包括各项守恒定律、经典力学定律和电磁学定律。

通过使用有限元法(FEM)将数学模型离散化,可以得到相应的数值模型;随后求解离散方程,并对结果进行分析,这就是有限元分析 这一术语的含义。

通过数学语言对物理定律在空间和时间进行表述,即产生了偏微分方程(PDE)。偏微分方程的解用因变量表示,如结构位移、速度场、温度场和电势场,等等。解是基于自变量 x、y、z 和 t 在空间和时间尺度上进行描述的。

求解给定系统的偏微分方程,不仅可以帮助我们理解我们所研究的系统,还可以对其做出合理的预测。有限元分析主要用于理解、预测、优化以及控制产品或过程的设计和运行。

通过有限元分析软件生成的冷却风扇振动模型。 采用结构力学振型分析预测得到的冷却风扇叶片振动情况。 采用结构力学振型分析预测得到的冷却风扇叶片振动情况。

对于结构力学而言,物理描述基于力平衡定律以及将应力应变关联的本构关系,胡克定律 便是这种关系的典型示例。传统的有限元分析通常仅限于结构分析。

CFD 仿真示例。 CFD 仿真示例。 CFD 仿真示例。

流体流动、传热和传质中的描述基于动量、质量和能量守恒定律,其中的通量往往由对流和耗散或扩散组合而成。耗散和扩散的精确形式由某个本构关系给出,例如,牛顿流体的黏性应力、傅里叶热传导定律和菲克扩散定律。

微带贴片天线模型。
探针馈电的切角微带贴片天线在 GPS 频率范围附近的辐射图。箭头显示右旋圆极化,其分量几乎是径向的。近似球形的表面表示增益的大小。 探针馈电的切角微带贴片天线在 GPS 频率范围附近的辐射图。箭头显示右旋圆极化,其分量几乎是径向的。近似球形的表面表示增益的大小。

通过高斯定律可以将静电场与电荷量相关联,根据安培定律可以确定磁场与电流的关系,根据法拉第定律,可以使用麦克斯韦方程组分析动力学,比如时变磁场产生的感应电场。麦克斯韦自己的决定性贡献在于,他引入了时变电场产生的位移电流,从而完善了对电磁场波特性的描述。

MEMS 执行器仿真示例。

MEMS 执行器通过焦耳热效应使两个通电热臂(白色和黄色)发生热膨胀,热膨胀会产生结构位移,从而推动执行器工作。这是有限元分析中的多物理场典型案例。您可以观看以下视频,了解如何从头开始构建此模型:如何利用 COMSOL Multiphysics 建立并运行仿真

MEMS 执行器通过焦耳热效应使两个通电热臂(白色和黄色)发生热膨胀,热膨胀会产生结构位移,从而推动执行器工作。这是有限元分析中的多物理场典型案例。您可以观看以下视频,了解如何从头开始构建此模型:如何利用 COMSOL Multiphysics 建立并运行仿真

我们可以将结构力学、流体流动、传热、传质和电磁学中使用的定律结合起来,用于描述多种物理现象相互作用的系统,也就是通常所说的多物理场系统。例如,在描述微机电系统(MEMS)时,通常需要结合使用结构力学定律和电磁学定律。在描述流-固耦合(FSI)时,需要结合使用结构力学定律和流体流动定律。

数学模型和数值模型

系统的数学模型可以包含一个或多个偏微分方程来描述相关定律,同时包含边界条件和初始条件。边界条件用于对解和部分建模域(如表面、边或点)施加额外的条件,同一个模型可以使用多种不同的边界条件。初始条件用于定义系统在时变事件开始时的状态。

散热器的数学模型,突出显示了模型中使用的方程。 这是一个基于能量守恒建立的散热器数学模型,可以得到温度场 T 的热方程,即上图中的第一个方程。本构关系是傅里叶热传导定律,其中 k 表示导热系数。上述方程在散热器域(Ω)内部有效,其中 ρ 和 Cp 分别表示材料密度和热容。在散热器底座(∂Ω1)上,上图第二个方程中的温度(用 T 表示)被限制为 Thot,用来表示第一个边界条件。在除底座以外的其他所有边界(∂Ω2)上,每个点的热通量都与温度突变(突变到边界层外的环境温度)成正比。h 表示传热系数;n 表示向外指向边界表面的法矢,具有统一的长度。如果环境温度低于散热器的局部温度,则散热器会释放热量。 这是一个基于能量守恒建立的散热器数学模型,可以得到温度场 T 的热方程,即上图中的第一个方程。本构关系是傅里叶热传导定律,其中 k 表示导热系数。上述方程在散热器域(Ω)内部有效,其中 ρ 和 Cp 分别表示材料密度和热容。在散热器底座(∂Ω1)上,上图第二个方程中的温度(用 T 表示)被限制为 Thot,用来表示第一个边界条件。在除底座以外的其他所有边界(∂Ω2)上,每个点的热通量都与温度突变(突变到边界层外的环境温度)成正比。h 表示传热系数;n 表示向外指向边界表面的法矢,具有统一的长度。如果环境温度低于散热器的局部温度,则散热器会释放热量。

从物理学角度来看,边界条件和初始数据通常是模型的自然组成部分(例如,结构力学中的载荷和约束,流体流动分析中入口和出口的压力水平,以及静电学中的终端电势)。从数学的角度来看,边界条件和初始数据能够从无数的可能解中选出唯一的解。

适定的数学模型

一个正确定义的数学模型往往具有适定性。如果一个数学模型具有唯一解,并且这个解连续依赖于问题的数据(即源项、通量、约束值和初始值),则该模型是适定的。如果模型 适定,则会在数值模型中反映出来,并会在求解过程中出现问题。

“适定”可以认为是模型能否用于数值仿真(例如有限元分析)的最低要求。

从理论的角度来讲,通常很难确定现实中的非线性三维模型是否适定;鉴于此,用于基础分析的模型都经过了大量简化。通过这些简化模型得出的结论可以用来评估更贴近实际的模型的性能表现。即使是适定模型,也可能对模型数据的变化非常敏感,这类模型在本质上是病态或敏感模型。

通过现代数值方法求解偏微分方程,实现了数学应用的革命性突破。原因在于,只有在非常特殊的情况下(例如方程与简单几何的特定组合中)才能得到数学模型的解析解。尽管这些情况从理论上来说非常重要,但对工程师而言用处并不大。数值方法突破了这个限制,可以处理非线性问题和复杂的几何结构。虽然数值方法还存在其他方面的计算难题(请参见下面的“”部分),但对新模型和几何的适用性没有任何问题。

数值方法可以给出适定数学模型的近似解。大部分数值方法都以建模域和所描述因变量的离散为基础。有限差分法、有限体积法和有限元法是最常用的离散化方法。顾名思义,有限元法(FEM)用于进行有限元分析。

离散化的散热器模型。 有限元离散的散热器模型。底座的四面体有限元体网格构成三角形表面单元。翅片内部的棱柱单元构成翅片表面的矩形单元。 有限元离散的散热器模型。底座的四面体有限元体网格构成三角形表面单元。翅片内部的棱柱单元构成翅片表面的矩形单元。

对所描述系统的数学 模型进行离散化,可以得到对应的数值 模型,后者是前者的离散近似。使用数值模型代替数学模型会引入误差,这种误差称为截断误差

截断误差被定义为数值模型与数学模型的解之间的差值。如果数值模型稳定且一致,则截断误差接近于零,这是因为单元尺寸接近于零(即数值解收敛于数学模型的解)。截断误差会以一定的速度收敛;速度由精度阶数测得。如果精度阶数为正数,则说明模型具有一致性。

有限元法的出发点是数学模型的弱形式。通过引入试函数,并将偏微分方程与这些试函数相乘,然后在建模域中对其进行积分,即可从逐点偏微分方程得到这种形式(也称为强形式)。您可以选择将这一过程与分部积分法结合使用。每个试函数都必须保持积分关系。为了与偏微分方程一致,必须有无穷多的试函数,并且这些试函数必须具有普遍的适用性。如此一来,就必须保持无穷多个积分关系,而建模域中每个点的逐点偏微分方程也必须保持不变。

有限元法引入的试函数通过计算网格进行定义。每个计算单元或网格单元都有多个局部定义的试函数。此外,有限元法中还定义了形函数,作为其基本组成部分。这些形函数用来表示候选解。对于瞬态问题,有限元法往往只用于空间离散化。在这种情况下,也就是在有限元离散之后得到的方程组是常微分方程组(ODE)。这个方程组转而使用有限差分法或其他类似方法进行离散。

在这个示例中,温度的近似公式可以写为:

(1)

其中, 表示形函数, 表示未知权重。

这里,形函数只在空间上变化,而未知数只在时间上变化。在有限元法中,权重也称为自由度。在伽辽金有限元方法中,形函数与试函数具有相同的类型。然后,结合使用弱形式与边界条件为未知数建立有限代数方程组。

网格中的每个节点(用 j 表示)都有一个基或试函数 ,并根据以下公式给出方程:

(2)

每个节点 j 都只有几个形函数 i 与试函数 j 重叠;因此,上面的积分只能在共享同一节点 j 的单元上进行计算。请注意,如果节点 j 是一个边界节点,则边界贡献不为零。对未知数而言,这是一个线性常微分方程——至少在材料属性与温度无关的情况下如此。

根据一般的经验,如果偏微分方程是非线性的,则未知数问题通常也是非线性的。还请注意,这个方程组只包含通量边界条件,不包含指定温度条件,您必须将其作为附加约束条件进行添加。

在为不同的网格单元类型和方程选择试函数和形函数方面,有限元法提供了相当大的自由度。“非结构化”网格的参与单元在尺寸和形状上往往存在较大差异;这种网格的有限元公式使得该方法能以相对较低的成本应用于非常复杂的几何结构,自动生成非结构网格也更加容易一些。在大多数情况下,形函数和试函数为易于定义和积分的低阶多项式。此外,与弱形式的密切关系为该方法提供了坚实的数学基础,其中的数学理论已经发展得非常成熟。

一旦将数学模型离散化,便必须求解生成的数值模型中的方程。

使用有限元分析软件创建的散热器模型。 与从左向右流动的周围空气的温度相比,散热器内部的温度变化幅度很小。由于散热器材料(铝)的导热系数远大于空气(并且空气混合不均匀),因此能推测出散热器内部的温度变化不大。在该仿真中,散热器模型通过方程得到了扩展,包括散热器周围空气的能量守恒方程,以及体现质量守恒和动量守恒的流体流动方程。整个方程组通过有限元法进行离散和求解。 与从左向右流动的周围空气的温度相比,散热器内部的温度变化幅度很小。由于散热器材料(铝)的导热系数远大于空气(并且空气混合不均匀),因此能推测出散热器内部的温度变化不大。在该仿真中,散热器模型通过方程得到了扩展,包括散热器周围空气的能量守恒方程,以及体现质量守恒和动量守恒的流体流动方程。整个方程组通过有限元法进行离散和求解。

有限元分析的具体过程

由于历史原因,传统的有限元分析对象主要是基于结构力学的模型,只在较小范围内涉及传热。随着多物理场建模应用的日益广泛,有限元法在流体流动和电磁仿真中得到了普遍应用,“有限元分析”一词也逐渐被其他工程和科学领域所接受和认可。事实上,无论是在何种应用领域,有限元分析过程都是相同的。

下面简要概述有限元分析的主要工作流程,涵盖从几何到模型文档的各个过程。

几何

有限元分析要求模型几何“紧密相连”。计算机辅助设计(CAD)几何结构并不总是用于分析目的。举个例子,现实生活中的某个物体可以通过 CAD 绘图中一组松散连接的三维表面来描述;然而在有限元分析中,这些表面必须能够形成一个真实的体。

即使 CAD 绘图中的一组三维表面能够形成一个体,但很可能存在一些表面过于细长,而一些边对于几何尺寸而言又过短的情况。这样一来,这些不理想的几何特征上便会产生不符合要求的单元分布。

为了准备一个可用于有限元分析的 CAD 几何,通常需要对几何进行修复特征去除。修复操作可以修补几何中不“紧密相连”的部分,特征去除操作可以移除细长表面或合并多余的小边。

为了进行分析而创建、修复 CAD 几何以及去除特征的过程,通常是更大过程的一个环节——在有限元分析中,传统上称之为预处理

材料

数学模型中的本构关系涉及材料的物理属性,这些属性可能与建模变量(“因变量”)相关。例如,在热膨胀分析中,机械属性和热属性往往与温度相关。

空气导热系数的样本图。 空气导热系数随温度的变化情况。 空气导热系数随温度的变化情况。

在实际操作中,需要正确估计材料属性和参考点的有效间隔。除此之外,还必须为几何的各个部分指派不同属性的材料。

在传统的有限元分析中,定义和指派材料属性及材料属性函数的过程,通常也认为是预处理的一部分。

域设置、边界条件、载荷及约束

在结构力学中,通过为系统选择的材料、载荷及约束可以定义数学模型。一般情况下,设置材料、域方程、边界条件和初始条件,即可定义数学模型。

散热器有限元分析:相关方程。 通过选择固体传热接口,可自动定义散热器中的域方程。方程描述包含温度相关属性。通过显示这些方程,您可以了解模型背后的原理。 通过选择固体传热接口,可自动定义散热器中的域方程。方程描述包含温度相关属性。通过显示这些方程,您可以了解模型背后的原理。

这部分的分析涉及选择几何域、边界、边和点,以及为这些几何实体指派方程、载荷或约束。这些设置的定义过程通常也认为是传统有限元分析中预处理的一部分。

网格

几何、材料、域设置、边界条件、初始条件、荷载及约束的定义无需离散化即可进行。然而,在许多较旧的有限元分析软件中,仍然是基于离散模型来执行这些操作。

网格创建完成后,我们便得到一个完整的数值模型。不同的现象和分析需要使用不同的网格设置。例如,在波的传播问题(例如,结构力学中的弹性波建模,或者射频分析中的电磁波建模)中,最大单元的尺寸必须远小于波长才能正确求解问题。在流体流动中,可能需要边界层网格才能解析边界层,而单元雷诺数可以确定流体本体的单元尺寸。

边界层网格剖分示例模型。 通过适用于流体流动建模的边界层网格剖分方法,对固体表面网格进行细化。 通过适用于流体流动建模的边界层网格剖分方法,对固体表面网格进行细化。

在许多情况下,CAD 几何的不同部分必须单独进行网格剖分。各部分之间界面处的模型变量必须通过有限元分析软件进行匹配,这一操作可以通过连续性约束(即用于将不同部分的有限元离散相互关联的边界条件)来实现。由于这些条件可能具有非局部特性,它们通常被称为多点约束

在传统的有限元分析中,网格剖分被认为是难度最高的预处理任务之一;而在现代有限元分析软件中,初始网格可以在求解过程中自动改变,从而尽可能减小数值解的误差,这种方法称为自适应网格剖分

如果说创建网格是一项高难度的任务,那么在合理的计算时间内选择和设置求解器并求得方程的解(构成数值模型)便是一项更加艰巨的工作。这些困难与求解过程中面临的各项挑战紧密相关。

首先,使用代数方程的离散模型可能非常大。一个三维模型往往拥有数百万个自由度。基于有限元法的数值模型求解过程的中心环节,是求解大型线性代数方程组。非线性、参数化、特征值和瞬态问题则通过迭代法求解,该方法求解一系列大型线性方程组。

一般来说,大型线性方程组很难进行有效求解。尽管可以使用现成的黑箱法,但这种方法对于实际模型而言往往代价过高,相关示例包括基于 LU 因式分解法的直接求解器,或通用的迭代法

为了寻求一个成功且近乎最优的替代方案,必须利用基础系统的某种结构。对于多物理场问题,这种结构可能并不存在,或者难以识别。在这种情况下,一个有用的做法是,将问题按物理成分进行分解,从而使分解的结构能被有效利用。现代有限元分析软件使用几何或代数多重网格法,以加快线性系统的迭代求解过程。

有限元分析求解器的另一个问题来源是模型的非线性。牛顿方法使用局部导数信息来寻求更好的候选解;但只有在当前的估计解与真正的数值解足够接近 时,这种方法才可靠。实际上,对解的初始猜测值并不总是与真实值足够接近;在这样的情况下,使用牛顿方法通常无法正确求解。对问题进行不同的简化或松弛 可能会有所帮助。通过求解更为简单的相似问题来代替原始问题,可以得到候选解。例如,可以忽略某些非线性来得到一个容易求解的线性问题。为此,我们开发了分离式求解器和连续求解器。

有限元分析求解器遇到的第三类困难是,数值模型可能不稳定,或者由于其他原因,没有为数学模型提供较好的近似。与更完善、性能更优的数值模型相比,这些情况下的求解过程要显得困难得多。我们可能很难发现和理解引起这一问题的根本原因。在许多情况下,我们可以采用某种方式修改模型而不是通过设置求解器作为补救措施。更好的自适应网格通常是模型性能改进的重要组成部分。

总而言之,求解器的设置选项需要非常灵活;同时,还要能调用功能强大的各种方法。我们往往需要在稳定性和性能之间达成折衷。

通过扰乱数值模型来研究其灵敏度,始终是一个很好的做法。可以到达两个相关目的:检查模型是否具有数值上的稳定性,以及对当前有限元分析中的一些重要物理量的截断误差进行量化。有限元离散的截断误差通常(但不总是)决定着仿真中的误差。为此,可以将模型中的关键派生值与典型的网格单元尺寸进行比较,然后采用不同的网格(在理想情况下,采用与当前网格明显不同的其他两个网格,且它们彼此也存在较大差异)重复进行计算。

如果数值模型性能良好,便可以根据比较结果估计精度阶数。如果精度阶数为正数,则所研究的两个最精细网格之间的物理量的差值,可以作为该物理量的截断误差估计值。有时,我们很难创建出满足所有需求的多个网格;在这种情况下,需要使用两个截然不同的网格的比较结果作为替代方法。

如果差值较小,表明该数值模型表现良好,且这个物理量的截断误差也较小。但是,如果差值较大,则很难得出任何结论。这种差异可能是由于两个网格的不稳定性或精度低造成的。然而,细化网格的解可能是精确的,而粗化网格的解则可能精度较差。即使截断误差估计尚不确定,至少可以排除模型不稳定的风险(例如,使用与该细化网格相当但不同的其他网格)。其他离散参数和求解器设置(如容差)也应进行更改。如果对于所有“扰动”仿真来说,研究的物理量只发生了很小程度的变化,则说明该数值模型具有良好的稳定性。

结果

对数值求解器的计算结果进行的分析包括:研究建模场的三维绘图、横截面图(如 x-y 绘图)以及计算派生值,例如对体、表面或边求积分,或计算沿边或点的表达式的值。

在较旧的有限元分析软件中,必须先定义要分析的绘图和派生值,才能进行求解。如果遗漏某些关键定义,则意味着需要从头开始重新求解。因此,定义要在后处理 过程中分析的表达式和派生值,也认为是预处理的一部分。

现代有限元分析软件支持在计算出解之后,动态定义表达式和派生值。在这些软件中,表达式和派生值的定义是后处理的重要组成部分,用于对模型进行深入预测。

后处理示例:用于进一步分析的散热器三维绘图。 基于表达式和派生值生成三维绘图、表面图、x-y 绘图和表格值,然后进行分析,是正确的后处理操作。 基于表达式和派生值生成三维绘图、表面图、x-y 绘图和表格值,然后进行分析,是正确的后处理操作。
描述散热器周围压力损失的压力图。 翅片顶部(蓝色)和底部(绿色)沿流动方向的压力损失。由于流体必须通过横截面相对较大的散热器底座,因此底部的压力损失略大一些。 翅片顶部(蓝色)和底部(绿色)沿流动方向的压力损失。由于流体必须通过横截面相对较大的散热器底座,因此底部的压力损失略大一些。

后处理中的一项重要任务是估计数值解中的误差。如上所述,可以通过求解不同网格尺寸的数值模型方程,来估计数值解的收敛性,从而实现误差估计。

后处理的另一个重要部分是估计模型对不同数据(如材料属性、初始条件、边界条件、载荷、约束以及数学模型和数值模型所需的其他输入参数)的灵敏度。

自动生成模型文档

在运行仿真后,非常重要的一步是将输入数据和仿真结果汇总到报告中,并在其中记录特定的会话。现代有限元分析软件支持定义报告的结构,用户可以在其中选择要记录的输入和输出数据。系统可以自动生成此类报告,您可以将其另存为文档,在将来每次仿真时用作参考。

自动生成模型文档的示例报告。 散热器仿真报告的第一页。报告结构创建完成后,报告便会根据每次仿真结果自动更新,并能以不同的名称进行保存,以记录仿真信息。其中包含一个问题定义部分,记录域设置、边界设置、初始条件、网格、自由度数量等信息。结果包含派生值和模型文件中的绘图。 散热器仿真报告的第一页。报告结构创建完成后,报告便会根据每次仿真结果自动更新,并能以不同的名称进行保存,以记录仿真信息。其中包含一个问题定义部分,记录域设置、边界设置、初始条件、网格、自由度数量等信息。结果包含派生值和模型文件中的绘图。

有限元分析发展趋势

如上所述,有限元分析过程包含许多步骤。选择大量参数(这些参数用于控制求解过程)等许多细节操作已成功实现自动化,无需用户太多关注。现今的有限元分析软件与上一代产品相比,性能得到了显著提升,价格也明显降低,工程师和小企业也能购买使用。

然而,为了进一步发掘有限元分析的潜力,使其帮助人们将更好的工程设计变为现实,需要做的工作还有很多。算法和用户界面都在不断得到改进,对于在各自的特定应用领域使用有限元分析工具的工程师、设计者和研究人员,减轻他们的负担,使其不必花大量的时间和精力来研究计算方法的细节,是目前的一个重要趋势。新软件接口的开发工作正在进行中,希望能帮助有限元分析专业人员和应用专业人员一同构建专用的分析工具,使工程师能够专注于设计任务,而不必“时刻关注”不断变化的计算细节。

随着价格低廉的云计算资源成为现实,再加上安全的数据传输手段,设计项目中将引入越来越多的计算分析。数学建模和有限元分析软件已经在过去和现在取得了成功,下一代软件将实现质的飞跃。数值计算不仅能减少工程工作量,还能使分析更加精确,实现对从概念到生产的整个产品链提供有力支持。借助有限元分析软件进行数学建模,必能照亮未来发展之路!

使用有限元分析软件创建的太阳镜仿真。 有限元分析具有光明的前景。 有限元分析具有光明的前景。
发布日期:2015 年 1 月 13 日
上次修改日期:2017 年 2 月 21 日