第2章：有限单元法及解题思路2.ppt

确定单元刚度矩阵的方法主要有直接刚度法和能量原理法。对于像三角形这样稍复杂的单元是较难采用直接刚度法确定单元刚度矩阵的，这时可采用弹性力学或材料力学中的虚功原理法。 虚功原理的含义是：当结构受载荷（外力）作用处于平衡状态时，在任意给出的结点虚位移下，外力（结点力）[F]所做的虚功等于内力（应力）所做的虚功，即?AF=?A? 3.7 结构刚度矩阵的快速建立法 前面已经通过平面桁架结构介绍了结构刚度方程的建立过程，即通过单元刚度方程写出每一个结点分力表达式，再根据变形协调条件、静力平衡条件写出各结点力学平衡方程，最后得出结构刚度方程。这是一个极其繁琐的过程，对于结点数较多的结构分析问题，其辛苦程度也是无法忍受的。事实上，通过单元特性分析，建立了单元刚度矩阵[K]e；通过单元载荷移置，建立了结点载荷列阵[F]e。在此基础上，可按照一定的方式，将单元的结点位移列阵[?]e组合成结构的结点位移列阵[?]；将单元的结点载荷列阵[F]e组合成结构的结点载荷列阵[F]；将单元刚度矩阵[K]e组合成结构刚度矩阵[K]；从而直接快速建立结构刚度方程：[F]=[K][?]，并且这种操作方法具有规律性，容易实现计算机编程。 上述3个方面的组合过程统称为集合。本节将再通过一个连续体结构的简单例子说明这一集合过程，并介绍由单元子刚阵Krs直接形成总体刚度矩阵[K]的方法。 ?在相互连接的公共结点处，诸单元的结点位移必须相等，即必须满足变形协调条件。对于下图所示结构，在公共结点i处的位移必须满足变形协调条件，即： ? ? ? ? i 因此结点位移不需要按单元来划分，结点i处的位移可写为?i=[ui,vi]T，结构结点位移列阵[?]只需按结点号递增的顺序直接写出，即[?]=[u1,v1,u2,v2,u3,v3,…] T； 在相互连接的公共结点处，诸单元对结点的作用力（即诸单元结点力的反力）与作用在该结点上的外载荷Ri之间必须满足静力平衡条件。对于下图所示的结构，在公共结点i处有： 3.7.1 集合的基本原则 0.5P 0.5P ? ? 1 4 3 2 图3－7受拉薄板模型 1.结构的结点位移列阵 根据公共结点处的变形协调条件，不同单元在公共结点处的位移相等。因此结构的结点位移列阵只需按结点编号顺序排列，即有： 3.7.2 结构刚度方程的快速建立实例 现通过一个简单例子来说明集合的过程。图3－7为一块受拉薄板，一端固定在两个铰链上，另一端作用有两个各为0.5P的集中载荷。为简单起见，将结构离散为两个单元，共4个结点，坐标系及结点编号如图示 2.结构的结点载荷列阵 在形成结构的结点载荷列阵时，需要注意以下3点： (1)当存在非结点载荷时，首先要进行单元载荷等效移置；在公共结点i处，需将与其有关的诸单元移置后的等效结点载荷Rie按分量进行叠加，得出结点i处的结点载荷Ri。 (2) 与结构的结点位移列阵?相对应，结构的结点载荷列阵R亦按照总体结点编号顺序排列。 (3)在位移型有限元法中，约束是通过限制结点位移来体现。因此不求解约束反力的情况下，在结点载荷列阵中不必考虑约束反力的作用。 综上所述，图3－7所示结构的结点载荷列阵为： 3.结构刚度方程 单元刚度矩阵的阶数与单元的自由度数相同，因结构刚度矩阵K是由单元刚度矩阵Ke集合而成，故结构刚度矩阵的阶数亦与结构的自由度相同。该结构有4个结点，每个结点有2个自由度，因此结构刚度矩阵是一个8×8的方阵。 下面介绍由单元子刚Krs矩阵形成结构刚度矩阵的两种常用方法。 ?按单元形成结构刚度矩阵 先将存放结构刚度矩阵的数组充零，然后从第1个单元开始，计算单元刚度矩阵[K]e,并将[K]e的每个元素存放到结构刚度矩阵 0.5P 0.5P ? ? 1 4 3 2 图3－7受拉薄板模型 的相应位置上。若在送入某单元刚度元素时某子刚阵位置上已有数值，则需要叠加上去。当依次作完最后一个单元时，便形成了结构刚度矩阵[K]。 3.8 结构刚度矩阵的性质 结构刚度矩阵是由单元刚度矩阵集合而成。因此它具有单元刚度矩阵的某些性质，如对称性、奇异性等。在进行有限元程序设计时，利用这些性质可以减少计算工作量，并节省计算机的存储量。 （1）结构刚度矩阵是一个对称方阵 因为单元刚度矩阵是对称方阵，因此，由单元刚度矩阵依次叠加而成的结构刚度矩阵必然也是对称方阵。在进行有限元程序设计时，利用这一性质，可以只计算及存储结构刚度矩阵的上三角或下三角，从而大大减少了计算工作量及对计算机存储要求。 （2）结构刚度矩阵是一个奇异矩阵 从物理上讲，在建立结构刚度矩阵的过程中，并没有对结构施加约束，因而没有消除结构