飞行器结构力学电子教案5-1 第五章 位移法.pptVIP

四、总刚度矩阵[K ]的组集 对于元素i j 的刚度矩阵， 遵照“按号就座”的原则，放置在总刚度矩阵[K ]中的{δ i } 和{δ j } 所在行和列的交叉点处。 注意：刚度系数在此叠加。 五、总刚度方程的求解 总刚度方程在没有引入位移约束条件之前，是不能求解的。这是因为，结构在没有引入位移约束条件之前是自由体，整体上存在3个(空间体6个)静力平衡条件，结点力列阵{P}中的项并不是完全独立的，或刚度矩阵[ K ]中的相关行和列有线性相关性，导致刚度矩阵的逆阵不存在，是奇异矩阵。 在结构分析中，采用引入边界条件的办法删去相关的行(列)。如果一个结构的刚体位移被限制，则与限制位移相对应的支反力自行保持全结构的平衡，也即这些支反力必然是其他外力的线性组合，而它所对应的行(列)，恰恰是线性相关的行列。位移边界条件处理之后，线性相关的行(列)就删去，处理之后的总刚度矩阵[ K ]变成可逆的，此时，总刚度方程才能求解。 对零位移边界条件，从总刚度矩阵中删去零位移所对应的行和列，这种方法称为删行删列法。 六、桁架的位移法求解 【例题1】试用位移法求图示桁架中3点的位移以及元素内力和支反力。设全部杆件的拉伸刚度均为EA。 解：以此例说明位移法的解题步骤。 1、结点编号、元素编号。 2、建立总体坐标系xoy 。 3、建立元素在总体坐标系xoy 中的刚度矩阵。 元素1-2，长度为L， 元素1-3，长度为 元素2-3，长度为L， 4、组集总刚度矩阵，建立总刚度方程。 按编号顺序列出结点位移和结点力列阵： 5、引入位移约束条件： 消除刚体位移。 简单做法：从总刚度方程中，删去等式两端零位移所在的行和列。 6、解引入位移约束条件后的总刚度方程，求出结点位移。 7、求支反力。 可以从总刚度方程中，求出支反力。 8、求元素内力，即杆轴力。 方法1：将已求出的结点位移代入到元素刚度方程中，可得到元素的结点力，然后合成得到杆的轴力。 对于元素i-j ： 注意：这里需要根据结点力的正负，来判断杆子是受拉、还是受压。 8、求元素内力，即杆轴力。 方法2：将总体坐标系中的结点位移，转化到元素局部坐标系中，在元素坐标系中计算杆的轴力。 对于元素i-j ： 并且：Nij 0，杆受拉； Nij 0，杆受压。 【例题2】用位移法求图示桁架中结点1、3的位移，以及元素内力，设全部杆件的拉伸刚度均为EA。 解： 1、结点编号、元素编号。 2、建立总体坐标系xoy ，考虑到对称性，取一半结构来分析。 3、建立元素刚度矩阵。 元素2-1， 元素2-3， 元素1-3， 4、组集总刚度矩阵，建立总刚度方程。 按编号顺序列出结点位移和结点力列阵： 考虑到位移约束条件： 因为要从总刚度方程中，删去等式两端零位移所在的行和列。因此，只需要针对 v1 和 v3 建立总刚度方程就可以了。 5、总刚度方程，求处结点位移。 6、求元素内力。 飞行器结构力学基础——电子教学教案 西北工业大学航空学院 航空结构工程系 第五章 位移法 Displacement Method of Structure Analysis 第一讲 位移法概述 杆元素与桁架的位移法求解 5.1 位移法概述 以结点位移(广义位移)作为基本未知量，写出由未知位移表示的应变，由物理方程写出仍由未知位移表示的应力表达式，最后由平衡条件解出所有的未知位移，这就是位移法的基本思路。 在计算机科学飞速发展的今天，适合于计算机应用的“有限元素法”正在逐步取代其他方法而成为飞行器结构分析方法的主流，并已发展为一门独立的新兴学科。本章所讨论的位移法，是以矩阵运算作为数学工具来处理结构位移计算的，故也称为矩阵位移法，它是有限元素法的基础。 矩阵位移法主要内容包括两个部分： (1)单元分析，即将结构分解为有限个较小的单元，进行所谓离散化。对于杆系结构，一般以一根杆件或杆件的一段作为一个单元，分析单元的内力与位移之间的关系，建立单元刚度矩阵。 (2)整体分析，即将各单元又集合成原来的结构，要求各单元满足原结构的变形协调条件和平衡条件，从而建立整个结构的刚度方程，以求解原结构的位移和内力。 在杆系结构中，若单元只受轴力作用，则称为杆元素，如桁架；若单元不仅受轴力，还受剪力和弯矩作用，则称为梁元素，如梁、刚架等。 由于杆元素和梁元素是最简单的元素，对这两个元素的分析，既有鲜明的物理意义，又能反映位移法的实质。所以，本章主要对杆元素和梁元素进行分析，并用于桁架和刚架的位移法求解。 5.2 杆元素与桁架位移法求解 本节将由最简单的杆元素和桁架开始，逐步介绍矩阵位移法的基本原理和计算过程 。 5.2 杆元素与桁