有限元方法研究.ppt

有限元方法Finite Element Method 一.概 念 有限元法是求解数理方程的一种数值计算方法。它是将理论﹑计算和计算机软件有机的结合在一起的一种数值分析技术。 已经实现商业化： 一批专用或通用有限元分析软件，美国Nastran、德国的ASKA、英国的PAFEC、法国的SYSTUS、美国的ABAQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。 1. 工程技术领域内的第一类问题 可以归结为有限个已知单元体的组合。例如，材料力学中的连续梁、建筑结构框架和桁架结构，把这类问题称为离散系统。 2. 工程技术领域内的第二类问题 通常可以建立它们应遵循的基本方程，即微分方程和相应的边界条件。例如弹性力学问题，热传导问题，电磁场问题等。由于建立基本方程所研究的对象通常是无限小的单元，这类问题称为连续系统。 以热传导问题为例： 热传导问题的控制方程与换热边界条件： 初始温度场也可以是不均匀的，但各点温度值是已知的： 某种边界条件： 1.典型的结构问题 (1) 应力分析，包括衔架和框架分析，与孔、凸起式其他物体几何形状改变相关的应力集中问题。 (2)弯曲问题 (3)振动分析 2.非结构问题 (1)热传递 (2)流体流动，包括通过多孔材料的渗流 (3)电位或磁场的分布 有限元法使用广泛的主要原因是其一般性，能够求解任何形状、在任意边界条件下的物体，并能够得到变形分布和详细的应力状态。 三.有限元法的基本思路 将连续系统分割成有限个分区或单元，对每个单元提出一个近似解，再将所有单元按标准方法加以组合，从而形成原有系统的一个数值近似系统，也就是形成相应的数值模型。 在有限元方法中，控制方程按照单元的组合重写，成为齐次线性方程组，也就是通过求解物理上近似的有限元表示的控制方程，得到数学上的准确解。 有限元法在数学上用偏微分方程组表示变分分析，能够被应用于由偏微分方程组控制的任何现象的求解，如流体、热传导和电磁学问题等。 四.有限元法的计算步骤 对求解区域进行离散化，通常是剖分成一系列的三角形网络，或其他形状网格，并对所有的单元与节点进行编号； 根据原微分方程及其边界条件建立单元内部的泛函表达式 建立单元内部的插值函数 对单元进行变分计算 将各单元的方程进行总体合成，构成代数方程组 求解代数方程组，计算出各离散点上所求物理量的近似函数值 概括为：网格划分、单元分析、整体分析 1.网格划分 有限元法的基本做法是用有限个单元体的集合来代替原有的连续体。首先要对弹性体进行必要的简化，再将弹性体划分为有限个单元组成的离散体。单元之间通过节点相连接。由单元、节点、节点连线构成的集合称为网格。 离散后单元之间利用单元间的节点相互连接起来。单元节点的设置、性质、数目等应视问题的性质、描述变形形态的需要和计算进度而定。 用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理，则所获得的结果就与实际情况相符合。 通常把三维实体划分成四面体或六面体单元的实体网格，平面问题划分成三角形或四边形单元的面网格。 网格划分遵循的主要规则 互不重叠。如是三角形的网格，不能把一个三角形的定点取在另一个三角形的边上。 当边界是曲线时，则在划分时近似的用直线代替并取为网格的一边。网格形状和大小可任取，但长度不宜相差过大，以免影响计算精度。 需要把区域划分成内部单元和边界单元。 对所有的内部单元的节点编号位置可以任意，但是要遵循逆时针排列顺序。 规定边界单元只能有一条边在边界上，且在确定节点编号时将此边相对的节点作为本单元的第一个节点。 所有的单元与节点均须编号，单元中节点的编号应与周围节点的编号尽可能接近，同一单元中的节点编号相差不应太大。 单元编号要求先编内部单元，然后再编边界单元。 如果求解区域由两种或更多种的材料组成，则一个单元中只能包含一种材料，而且必须确定单元的材料种类。 网络自动生成及其计算机编程 利用计算机来实现自动生成的关键在于对单元作规格化处理：即找出不同单元之间的共同规律，对研究区域进行分层处理。 2.单元分析 根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等，找出单元节点力和节点位移的关系式，这是单元分析中的关键一步。此时需要应用弹性力