连续体结构分析的有限元方法—平面问题有限元分析.ppt

4节点矩形单元的线性应变和应力 由单元的位移表达式可知，4节点矩形单元的位移在x，y方向呈线性变化，所以称为双线性位移模式，正因为在单元的边界x=±a和y=±b上，位移是按线性变化的，且相邻单元公共节点上有共同的节点位移值，可保证两个相邻单元在其公共边界上的位移是连续的，这种单元的位移模式是完备(completeness)和协调(compatibility)的，它的应变和应力为一次线性变化，因而比3节点常应变单元精度高。 4.4 三角形单元与矩形单元计算精度的比较 第4章 连续体结构分析的有限元方法 杆梁结构由于有自然的连接关系，可以将其进行自然的离散，而连续体则不同，它的内部没有自然的连接节点，必须完全通过人工划分节点的方法进行离散，因此成为逼近性离散；有限元方法的真正意义或价值在于它成功地处理了连续体(场)问题。人们公认的有限元方法的开拓者之一Courant，就是在1943年处理连续体问题时使用了三角形区域的分片连续函数和最小势能原理；而“有限单元”的名称，是Clough于1960年在处理平面连续体问题时正式提出的。 本章将就连续体问题的有限元方法进行全面的讨论和研究。 4.1 连续体结构分析的基本力学原理 连续体问题的三大类变量 连续体问题的三大类方程及边界条件 4.2 平面问题的3节点三角形单元 平面问题3节点单元具有几何特征简单、描述能力强的特点，是平面问题有限元分析中最基础的单元，也是最重要的单元之一。 单元的几何和节点描述 该单元共有6个节点位移自由度(DOF)，将所有节点上的位移组成一个列阵，记作qe；同样，将所有节点上的各个力也组成一个列阵，记作Pe，那么 利用函数插值、几何方程、物理方程以及势能计算公式，可以将单元的所有力学参量(u(x,y) ，ε(x, y) ，σ (x, y) 和Πe)用节点位移列阵及相关的插值函数来表示，下面进行具体推导。 单元位移场的表达 平面3节点三角形单元，每一个节点有两个位移，因此共有6个节点位移，考虑到简单性、完备性、连续性及待定系数的唯一确定性原则，分别选取单元中各个方向的位移模式为： 由节点条件，在(x=xi，y=yi)处，有 将上式代入可求解出其中的待定系数，即 上式中， 上式中的符号（1，2，3）表示下标轮换，如1→2，2→3，3→1同时更换。 重写位移函数，以节点位移的形式表示，有 单元应变场的表达 单元应力场的表达 单元的势能的表达 单元的刚度方程 常用的等效节点外载列阵 （平面3节点三角形单元） 单元在承受非节点载荷时，如在边线上承受一个分布载荷，这时应根据外力功的计算公式来获得节点载荷的等效值，常见的平面问题3节点三角形单元的节点等效外载荷列阵如下表所示。 平面3节点三角形单元的位移坐标变换问题 由于该单元的节点位移是以整体坐标系中的x方向位移u和y方向方向位移v来定义的，所以没有坐标变换问题。 平面3节点三角形单元的常系数应变和应力 由于该单元的位移场为线性关系式，且系数ai ， bi ， ci只与三个节点的坐标位置（xi，yi）相关，是常系数，因而求出的单元的B(x, y) 和S(x, y) 都为常系数矩阵，不随x、y变化，由相应公式可知，单元内任意一点的应变和应力都为常数，因此，3节点三角形单元称为常应变(应力)CST单元(constant strain triangle)。在实际使用过程中，对于应变梯度较大（也即应力梯度比较大）的区域，单元划分应适当加密，否则将不能反映应变(应力)的真实变化情况，从而导致较大的误差。 4.3 平面问题的4节点矩形单元 单元的几何和节点描述 单元位移场的表达 单元应变场的表达 单元应力场的表达 单元势能的表达