第3章_基于虚位移原理的变分法(1).pptVIP

* 在前述有限元的直接刚度法里，我们采用材料力学中的叠加原理来求解直梁和刚架的单元刚度矩阵。但对于复杂的二维平面甚至三维立体的力学问题，就无法采用直接刚度法来求解单元刚度矩阵，此时我们可以采用基于虚位移原理的变分法。 1、背景 2、基本概念 2.1 虚位移原理（虚功原理） 在材料力学中，虚功原理表述为：在虚位移中，外力所做虚功等于内力在相应的虚变形上所做的虚功（虚应变能）。 2.1 虚位移原理（虚功原理） 在材料力学中，虚功原理表述为：在虚位移中，外力所做虚功等于内力在相应的虚变形上所做的虚功（虚应变能）。 2.2 虚位移和变分 在材料力学中，我们规定实位移是v，对应的虚位移就是v*。但严格来讲，虚位移应该用变分符号表示为δv。变分与微分有相似的地方，都具有无限小变化的意思。但严格来讲，二者是两个完全不同的数学概念。 3、应用基于虚位移原理的变分法求解直梁的单元刚度矩阵 3.1 位移函数 假设图示直梁单元的位移函数v(x)（插值函数）为： 已知单元在两个节点的边界条件： 实际上对于直梁单元，v(x) 就是在其上各点的挠度f 。 根据材料力学可知，截面转角和挠度之间的关系是 。 将上述v(x)的表达式（位移函数）代入四个边界条件得到如下四个方程式。 将所有的ai代入前述v(x)的表达式得到插值形式的位移函数。 解得 写成矩阵形式 式中 上述位移函数可进一步缩写成 式中 ； 3.2 形状函数 上式中的Ni（i = 1, 2, 3, 4)叫做形状函数，有时简称为形函数。在梁单元中,它表示一个两端固定的梁只产生一个单位位移时梁弯曲成的形状（如下图所示）。其性质将在下一章讲解。 3.3 用虚位移原理求单元刚度矩阵 在平衡状态下，梁单元有挠曲线v(x)，其节点位移分别为fi、θi、fj、θj 。给该单元的节点任意虚位移δfi 、 、 、 ，由此在挠曲线各点上产生相应的虚位移 和虚转角 （即虚功原理中的虚变形）。根据虚功原理得到以下方程式。 由挠曲线的近似微分方程，弯矩可写成 因为虚位移具有和实位移相同的性质， 因此以上方程也适用于虚位移。故 上述虚功原理方程式变为 式中单元节点虚位移列阵 单元节点力列阵 位移函数 其变分为 将以上两式代入 得到 将上式与 相比较，得到单元刚度矩阵的表达式。