有限元-第1章.pdf

第二部分 有限单元法 第一章 有限单元法概述 § 1-1 引言 有限单元法作为固体力学的一种分析方法是在本世纪五十年代起源于航空工程中飞 机结构的矩阵分析方法。结构矩阵分析是把一个结构看成由许多元件互相连接而组合成的 集合体，通过对元件的受力分析，建立节点位移与节点力之间的关系，再将这些关系集合 起来形成结构方程组。根据选取的基本未知量是节点位移还是节点力，有位移法、力法和 混合法。结构矩阵分析的对象限于由杆、梁、受剪板等元件组成的结构。 1960 年Clough R. W. 等人将这种处理问题的方法推广用来解弹性力学的平面应力问 题，并第一次采用“有限单元法”这个术语。应用有限元法对任意连续体进行分析时，首 先将连续体划分成有限个单元，并在每个单元上指定有限个节点，认为相邻单元在节点处 相互连接构成一组单元的集合体。用以模拟或逼近原来的连续体。然后，由对单元的分析 和集合，得到描述该离散结构的代数方程组。 常规的结构矩阵分析法是将每个元件的力与位移之间的关系精确推导出来。而将连续 体离散为单元的有限单元法，是选定场函数的节点值，例如取节点位移作为基本未知量， 对于每个单元根据分区近似的思想，在单元内假设近似的位移插值函数，利用弹性力学的 变分原理建立节点力与节点位移之间的关系，得到一组以节点位移为未知量的方程组。有 限单元法是一种近似的数值方法，显然，如果单元满足问题的收敛性要求，那么随着求解 区域内单元数目的增加，单元尺寸的缩小，近似解将收敛于精确解。 从有限单元法所依据的变分原理来看，早期的有限元法大都依据最小位能原理，以位 移作为基本位知量。后来，有依据最小余能原理的有限元法，以内力作为基本未知量。再 后来则有许多某种形式的广义变分原理，同时将位移和内力作为独立的基本未知量。还有 将这些变分原理结合起来应用，例如在每个单元内用最小余能原理而对整个系统用最小位 能原理求近似解，这就是所谓杂交应力有限元法的基本思想。基于最小位能原理的有限元 位移法是用得最广的一种方法，本教材在第四章介绍混合杂交有限元法的基本概念和基本 理论外，其余各章都采用有限元法的基本理论和方法。 § 1-1 有限单元法进行结构分析的步骤及有限元列式 有限元位移法进行结构分析的过程概括起来可以分成以下几个步骤。 一、结构的离散化 所谓离散化，就是将连续体结构划分成有限个单元，并在单元的指定点设置节点，使 相邻单元的有关参数具有一定的连续性，构成单元的集合体，用来代替原来的结构。为了 使该离散的计算模型有效地逼近实际的连续体，就需要考虑合理地选择单元的形状和划分 网格的方案，确定单元和节点的数目、节点的位置和自由度等问题 。 1、单元的几何形状。选择单元的形状既要考虑单元的精度和计算的简便，又要考虑对边 界形状的适应性。如对二维问题，单元的形状有三角形、矩形、任意四边形……；对 于三维问题，有三棱锥的四面体、三棱柱和五面体、长方体……等。 2 、单元的划分和节点的布置。节点的布置与单元的划分是相互联系的。通常集中载荷的 作用点、分布载荷集度的突变点、支承点、边界的转折点等都应取作节点，如果物体 的厚度有突变或者物体由不同材料组成，在布置节点时，不要把厚度不同或材料不同 的区域划在同一个单元中。至于节点的多少和分布的疏密、也即单元的大小，要根据 对计算精度的要求和计算机容量等因数综合考虑，单元划分得越细计算精度越高，但 也使计算时间及准备工作增加。因此，在保证足够精度下力求减少单元的数目。故在 划分单元时，对应力变化急剧的区域单元网格应划分得细些，应力变化平缓的区域单 元可划分粗些。单元节点的设置应使单元的几何形状能够容易地由节点位置坐标来描 述。一些在单元的端点或角点上设置节点，还可在角点之间的边上增设节点，这样可 描述单元边界是曲线的情况，用二次或更高次的曲线来描述单元的几何特征。 3、节点自由度数的确定，节点的自由度数应能正确地模拟单元的力学特征和易于构造单 元的近似解，保证相邻单元在公共边界上的连续性。 综合起来，结构离散化应符合几何近似和物理近似的原则。有限元的计算模型在几何 上应近似于真实结构。各单元的物理特性和受力特性应与真实结构在单元所在区域的情况 近似。 二、选择位移模式 在对典