结构稳定理论与设计-5修改.pptVIP

5.7 有限单元法（刚度矩阵法） 有限单元法 1. 先将构件划分为有限数量的单元（分段），将分段点的位移作为未知量； 2. 建立各单元两端的位移与内力之间的关系； 3. 利用各分段点处力的平衡和变形（位移）协调条件将各单元连接起来形成原构件。 各单元两端位移与内力之间的关系用矩阵表示，求解时可以得到足够精确的解。 单元刚度矩阵 整体刚度矩阵 5.7 有限单元法（刚度矩阵法） 结构力学中用直接刚度法建立了受弯构件内力与变形的关系式；对于受压（压弯）杆，建立单元刚度矩阵时，应考虑压力产生的附加内力（弯矩、剪力）的影响。 加载过程可视为: 第1步：加P引起轴向压缩且产生附加弯矩； 第2步：加q（杆端力），引起杆端变形。 5.7 有限单元法（刚度矩阵法） 只要变形很小且材料服从虎克定律，则其受力与加载的次序无关。 因此杆端位移矩阵[?]可考虑为第1步由P引起（只加P）；第2步由杆端力q引起（再加q），即[q]与[?]的关系是线形的。 单元刚度矩阵可用能量法求解。 1 单元刚度矩阵 根据能量法，因为在单元完成第2步时是平衡的，所以有内外虚功相等 有 因为 [q]= [k] [?]， [k]为单元刚度矩阵，故有 显然，若上式中的P=0，则相当于梁单元，即有： （1） 考虑压杆弯曲变形后压力所产生的附加内力时，有 [k]=（[k0] –P [k1]） [k1]：为当P=1时，由于P的作用而增加附加弯矩而引起杆端出现的附加反力与杆端位移的关系。即，杆端力q包括两部分，一部分是P=0时的杆端力；一部分是由P引起的附加弯矩所产生的附加反力，即有 [q]=（[q0] + [q1]） 矩阵式（1）中，每一列元素相当于当?i=1而其余杆端变形为零时的杆端力，故[k1]中的各列元素，亦应等于在压力P=1下?i=1而其余杆端变形为零时，由于附加弯矩而引起的杆端附加力。 假设挠曲线y为一个三次多项式（满足剪力为常数、无节间荷载时弯矩线形变化的条件），即 ？？ （2） 按图示坐标系，其边界条件为： 代入（2）式可解得待定常数： 即 因附加弯矩M1=Py = y（P=1），当仅?1=1而其余杆端变形为零时： 化为与其等效的分布荷载，有 按图示受力可求出杆端力，写 成矩阵的形式，即 q1 q2 q3 q4 当仅?2=1而其余杆端变形为零时： 等效分布荷载为： 写成矩阵的形式： 第2、3列同样按此方法导出，得： q1 q2 q3 q4 因[k]=（[k0] –P [k1]），即 对其加以简化，使有相同的量纲，即杆端力矩阵全部采用力的量纲，位移矩阵全部采用距离的量纲。 则因为 [q]=[k] [?]，有 以上即梁柱单元的单元刚度矩阵，由两部分组成。 第1部分是梁的刚度矩阵；第2部分系考虑轴压力的附加影响，亦称为初应力矩阵，也可用于轴向拉力的情况。 2 整体刚度矩阵 整体刚度矩阵可通过单元刚度矩阵与转换矩阵相乘而得。以两个单元的压杆为例，其节点位移用矩阵表示为： 节点外力用矩阵表示为