有限元分析离散方法.doc

有限元分析离散方法 离散化方法：有限差分法，有限元法，有限体积法 所谓区域离散化， 实质上就是用一组有限个离散的点来代替原来连续的空间。实施过程是;把所计算的区域划分成许多互不重迭的子区域，确定每个子区域的节点位置及该节点所代表的控制容积。节点：需要求解的未知物理量的几何位置；控制容积: 应用控制方程或守恒定律的最小几何单位。一般把节点看成是控制容积的代表。控制容积和子区域并不总是重合的。在区域离散化过程开始时，由一系列与坐标轴相应的直线或曲线簇所划分出来的小区域称为子区域。网格是离散的基础，网格节点是离散化物理量的存储位置。 大家都知道，常用的离散化方法有：有限差分法，有限元法，有限体积法。 1.有限差分法是数值解法中最经典的方法。它是将求解区域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程（控制方程）的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组。这种方法发展比较早，比较成熟，较多用于求解双曲线和抛物线型问题。用它求解边界条件复杂、尤其是椭圆型问题不如有限元法或有限体积法方便。 2.有限元法是将一个连续的求解域任意分成适当形状的许多微小单元，并于各小单元分片构造插值函数，然后根据极值原理（变分或加权余量法），将问题的控制方程转化为所有单元上的有限元方程，把总体的极值作为各单元极值之和，即将局部单元总体合成，形成嵌入了指定边界条件的代数方程组，求解该方程组就得到各节点上待求的函数值。对椭圆型问题有更好的适应性。有限元法求解的速度较有限差分法和有限体积法慢，在商用CFD软件中应用并不广泛。目前的商用CFD软件中，FIDAP 采用的是有限元法。 3.有限体积法又称为控制体积法，是将计算区域划分为网格，并使每个网格点周围有一个互不重复的控制体积，将待解的微分方程对每个控制体积积分，从而得到一组离散方程。其中的未知数十网格节点上的因变量。子域法加离散，就是有限体积法的基本方法。就离散方法而言，有限体积法可视作有限元法和有限差分法的中间产物。 4.有限分析法：同有限差分法一样，用一系列网格线将区域离散，所不同的是每个节点与相邻8个邻点组成。在计算单元中把控制方程中的非线形项局部线形化，并对该单元上未知函数的变化型线作出假设，把所选定型线表达式中的系数和常数项用单元边界节点上未知的变量值来表示，这样该单元内的被求问题就转化为第一类边界条件下的一个定解问题，可以找出分析解；然后利用这一分析解，得出该单元中点及边界上8个邻点上未知值间的代数方程，此即为单元中点的离散方程。两种离散方法外节点法：节点在子域的四角，先定节点位置而计算相应的界面内节点法：节点在子域中心，子域与控制容积重合。 计算时先定界面后算出节点位置。 5.边界元法(Boudarv element method, BEM)上面四种方法都必须对整个区域作离散化处理，用分布在整个区域上的有限个节点上函数的近似值来代替连续问题的解。在边界元方法中应用格林函数公式，并通过选择适当的权函数把空间求解域匕的偏微分方程转换成为其边界上的积分方程，它把求解区中任一点的求解变量(如温度)与边界条件联系了起来。通过离散化处理，由积分方程导出边界节点上未知值的代数方程。解出边界上的未知值后就可以利用边界积分方程来获得内部任一点的被求函数之值。边界元法的最大优点是，可以使求解问题的空间维数降低一阶，从而使计算工作量及所需计算机容量大大减小。边界元法推广应用的一个最大限制是，需要已知所求解偏微分方程的格林函数基本解。虽然对不少偏微分方程这种基本解业已找出，但对Navier-Stoles方程这样的非线性偏微分方程，至今尚未找到其基本解。目前的一种处理方式是，把Navier-Stokes方程中的非线性项看作是扩散方程的源项并通过迭代的方式来求解，但一般只能获得Re较低情形的解。最近文献中采用高阶涡量一流函数方程的边界元方法，已使获得顶盖驱动流稳定解的Re高达10000. 本文来自: 博研石油论坛 详细出处参考/thread-11065-1-4-1.html 有限元法 finite element method 有限元法是一种高效能、常用的计算方法．有限元法在早期是以变分原理为基础发展起来的，所以它广泛地应用于以拉普拉斯方程和泊松方程所描述的各类物理场中（这类场与泛函的极值问题有着紧密的联系）。自从1969年以来，某些学者在流体力学中应用加权余数法中的迦辽金法(Galerkin)或最小二乘法等同样获得了有限元方程，因而有限元法可应用于以任何微分方程所描述的各类物理场中，而不再要求这类物理场和泛函的极值问题有所联系． 基本思想：由解给定的泊松方程化为求解泛函的极值问题。 方法运用的基本步骤： 步骤1：剖分: 将待解区域进行分割，离散成有限个元素的集合．元素（单元