第二章 弹性力学平面问题有限单元法.ppt

第二章 弹性力学平面问题的有限单元法§2.1 有限单元法的基本概念和基本思路 弹性力学问题的解析法是把弹性体看成由无限多个、无限小的微元体组合而成，通过对其中任一个微元体的分析，建立对弹性体内各点都适用的基本方程。然后用严格的数学方法求解这些方程并设法满足边界条件，从而 得到精确的解析解。 有限单元法则首先作了物理上的近似。它把连续弹性体近似看成由有限多个、有限大小、彼此只在有限个点相联接的单元组合而成，即近似把实际物体分割成有限个有限大小的多边形(对平面问题)、多面体(对空间问题)或杆件(结杆系结构)，这些多边形、多面体或杆件就称为有限单元(简称单元)。各单元之间只在结点处相联接，而单元周边的互相分离的。 现在所述的有限单元法(实为有限单元位移法)取全部结点的(广义)位移作为基本未知量。这样就把原来是无限多自由度的体系简化成有限多自由度的体系。以上过程称为连续体的有限单元离散化。 为了求得全部结点上的位移，并最终求得各单元内的位移，应变和应力，有限单元法是从任一个典型单元开始分析。有限单元法又作了一次数学上的近似：构造适当的位移函数(称为位移模式)，把单元内任一点的位移用单元的结点位移来表示。有限单元法认为每个单元内部符合弹性力学的基本假设，严格运用弹性力学的基本方程的几何方程和物理方程，从而建立起单元内任一点的应变、应力分别与单元结点位移的关系式，为最后当求出结点位移后，再顺利求得单元内的应变和应力作好了准备。单元的平衡条件是用虚功方程代替的，由此得到重要的单元刚度矩阵和单元刚度方程，即单元结点力和结点位移的关系式，为进一步的整体平衡分析作好准备。这个过程称为单元分析。 最后考虑整个体系的平衡。由于所有的单元上都已没有外载荷 (全部移置到结点上了)，全部结点的平衡就代表了整体的平衡。对所有的结点建立结点力与结点载荷 平衡的关系，然后利用单元刚度方程把所有的结点力都换成结点位移的表达式，就建立了结点位移与结点载荷的关系式——总刚度方程。这是一个代数方程组，未知量就是结点位移，引入已知的位移边界条件改造这个方程组后，就可借助计算机求解而得到所有结点的位移。有了结点位移，再利用单元分析中的位移函数和已准备好的几何方程和物理方程的演变式，可以按需要容易地求得任一单元内任一点位移、应变或应力。这个过程称为整体分析。 §2.2连续体的有限单元离散和计算简图 对于平面问题，最简单而且也是最常用的单元是三角形单元和矩形单元。例如， 图2-1所示为一个托架梁，可以当做弹性力学的平面问题处理。我们见将其划分位三角形网格，各单元再结点处用光滑的 平面较连接，每各单元所受的载荷也移置到结点上，成为结点载荷。 在位移为零的结点处，或者在位移很小可以忽略的结点处，设置相应的支座铰链。在图（a）中，固定边AB各点的位移均为零。所以在AB边上的各结点处应设置固定支座铰链。这样就得出托架得有限单元计算简图，如图（b）所示。 把连续体进行有限单元离散时，首先要考虑到选用哪一种形状的单元。这一选择取决于结构的几何形状、计算精度的要求及描述该问题所必须的独立空间坐标的数目。对于平面问题，通常采用直角三角形和矩形单元，特别是三角形单元比较适宜于模拟有曲线边界的物体或结构。一般说来，单元各边的比例不能相差太大。计算实践表明，单元各边的比例相差太大是影响计算精度的一个重要因素。故应避免取狭长的单元。对三角形单元而言。在划分网格时应尽量使所有单元接近于等边三角形。但通常为了适应结构的边界形状及单元由大到小的过渡，很难实现这一要求，不过应尽可能的满足。 在划分网格时，就整体来说，单元的大小（即网格的疏密）要根据精度要求和计算机的速度及容量来决定。一般讲来，单元越小，网格越密，计算结果越精确。但是，单元越多，要求计算机容量就越大。因此，单元划分多少合适，一方面要考虑计算精度的要求；；另一方面要根据计算机的条件，应在计算机的容量范围内来决定单元的大小和数量。原则是，在保证必要的计算精度条件下，单元应尽量取的少些。 在单元的排列上，应根据计算者的实践经验对所计算的对象进行判断，在应力剃度变化大的部位和重要的部位，单元应取小些，网格也划分的密些。反之，在应力变化平缓的部位和不重要的部位，单元可取大些，网格也就稀些 在划分单元还应考虑到，当计算对象的厚度或者其弹性性质有突变之处，除了在这些部位单元应取小之外，还应把突变线作为单元的边界线。如果结构受有集度突变的分布载荷或集中载荷时，在这些部位的单元同样应当取小些，并且在载荷突变处和集中力处应布置结点，以使应力的突变得到一定程度的反映。 总之，把连续体进行有限单元离散而成计算简图，是综合运用工程判断力的过程，在这个过程中，要决定单