02-01平面问题的离散化.docVIP

第二章2009-11-12） 平面问题不仅在工程上具有实际意义，而且具有一定的典型性。工程实际问题中有许多问题可简化为平面应力问题或平面应变问题。如弯曲梁、水库大坝等等。使读者较容易地建立有限单元法的基本概念，掌握基本计算方法及步骤，以及使用中应注意的问题。本章将三角形常应变单元，介绍有限单元法的基本原理和方法。为进一步学习更加复杂问题的有限单元法打下基础。§2-1 平面问题的离散化 结构的离散化就是将任意形状的，受各种载荷和约束的连续体或结构划分成由有限个方位不同但形状相似，仅在节点处相连，仅在节点处受力和约束的单元的组合体。提取有限单元法分析所用的数据，为后续工作做准备。 有限单元计算所需用的数据包括： 节点的个数、单元个数、受力点数、约束点数； 节点坐标； 单元组成及其物理性质； 受力点号及受力值； 约束点号及约束类型。 实际上，结构的离散化过程，就是对真实结构的近似过程。为保证计算精度，要求离散体和真实结构应尽可能的达到几何近似和物理近似： 几何近似：即离散体的几何形状近似于真实结构。 物理近似：离散体内单元的物理特性(如受力、变形、材料等)近似于真实结构相应区域的物理性质。 这是结构离散化过程所必须遵循的总的原则。具体实施时，大致可分为如下四个步骤:①建立数学模型；②单元的划分；③载荷的移置；④约束的简化。下面将分别给予介绍。 建立数学模型 简化实际结构 保留几何形状、受力及约束的主要特征，而忽略次要的细节，使之既利于计算，而又不失其主要特性，满足一定的精度要求。例如，板带轧钢机机架，一般为一具有较大厚度，在上下中间部位受有轧制力的闭口框架。其在受力部位的厚度和受力分布是不规则的，另外，有时机架还受到横向力的作用。若是严格按实际结构分析则非常复杂。在一定精度要求下，可将其简化成如图(2-1a)所示的结构：不计厚度及受力沿厚度的变化和厚度较大等因素将其简化为平面应力问题；不计横向受力，只考虑轧制力，且将其看成作用于对称轴上的集中力。 在确定计算模型时，还要充分利用结构的对称性。在图(2-1a)中，机架的结构和受力均是对称的，所以其中的应力及位移也对称于该中间轴，且在对称轴上各点水平方向的位移为零。所以，可以取半个结构为计算简图。为保证其与原结构的同一性，在对称轴上各点加一附加约束，保证其水平方向不发生位移(图2-1b)。经过这样处理，可使同一问题采用相同单元数时所得解的精度提高；或满足同一精度要求时，减少单元总数，节省上机时数。 (a) (b) 图2-1结构与载荷对称于一个轴时的轧钢机机架的计算简图 图2-2结构与载荷对称于垂直的两个轴时的计算简图 当结构与载荷均对称于垂直的二个轴时(图2-2)，由于应力和位移也对称于该两轴，所以可以将结构沿此二对称轴分开，取1/4结构作为计算简图。这时，应加一些约束，以保证每个对称轴上各点在与该轴垂直的方向上的位移为零。在结构的厚度和材料不同时，还要在计算简图上画出不同厚度或不同材料的分界线。在以后划分单元时，应使单元不跨越此类分界线，以使每个单元内厚度和材料性质相同，便于处理(图2-3)。在分界线附近应将单元划得小一些，因为在材料或厚度突变的地方一般也伴随着应力的突变，应使单元能够反映这种变化。  (a) (b) 图2-3结构具有不同厚度或材料时单元划分简图 (a) 不同厚度的单元划分图；(b) 不同材料的单元划分图如果构件厚度是连续变化的，当梯度较小时，可用几个不同的平均厚度来代替；当厚度变化梯度较大时，就不能再按平面问题来计算了。应按空间问题处理。确定坐标系 本书系用右手坐标系。具体设置可由计算简图确定。但应尽可能使计算简图安置在XOY坐标系的第I象限，可以少出错误。 单元划分 即绘制网格图。在进行单元划分时，首先遇到的问题是单元类型的选择，表2-1列出了一些常用的单元。具体使用时应根据问题的具体类型、计算精度要求及手头可用程序等情况决定。 图2-4 三角形单元 对于平面问题，最简单最常用的是三节点三角形单元(如图2-4)，在平面应力问题中，它们是三角形板，而在平面应变问题中，则为三棱柱。三角形单元间通过节点以假想的铰链相互联接，它们只传递力，不传递力矩。当单元处于曲线边界时，用曲线的割来代替原曲线作为三角形的一条边。 为便于进行单元分析，在平面问题中，无论弹性体与材料厚度与材料如何变化，无论采用什么型式的单元，划分后的每个单元都看作是均匀连续的、完全弹性的、各向同性的等厚度平板。单元仅在节点处受力和约束。为达到上述目的，并保证计算的顺利进行及一定的精度，在单元划分时应注意以下几个问题： 在划分单元时，就整体来说，单元的大小(网格的疏密)应根据计算精度的要求和计算机的容量与速度来确