结构力学第8章.ppt

《结构力学第8章.ppt》由会员分享，可在线阅读，更多相关《结构力学第8章.ppt（23页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助结构力学第8章 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 矩阵位移法的基本思路是矩阵位移法的基本思路是:(1)先把结构离散成单元，进行单元分析先把结构离散成单元，进行单元分析,建立单元杆建立单元杆端力与杆端

2、位移之间的关系；端力与杆端位移之间的关系；(2)在单元分析的基础上，考虑结构的几何条件和平衡在单元分析的基础上，考虑结构的几何条件和平衡条件，将这些离散单元组合成原来的结构，进行整体分析，条件，将这些离散单元组合成原来的结构，进行整体分析，建立结构的结点力与结点位移之间的关系，即结构的总刚建立结构的结点力与结点位移之间的关系，即结构的总刚度方程，进而求解结构的结点位移和单元杆端力。度方程，进而求解结构的结点位移和单元杆端力。在从单元分析到整体分析的计算过程中，全部采用矩在从单元分析到整体分析的计算过程中，全部采用矩阵运算。阵运算。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计

3、算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 集成总刚度矩阵最常用的方法是直接刚度法，即由单集成总刚度矩阵最常用的方法是直接刚度法，即由单元刚度矩阵直接集成结构刚度矩阵，又可分为后处理法和元刚度矩阵直接集成结构刚度矩阵，又可分为后处理法和先处理法。先处理法。1.后处理法后处理法 (1)集成。集成。对所有单元不做边界条件处理，均采用自对所有单元不做边界条件处理，均采用自由式的单元刚度矩阵，按单元的结点编号将单元刚度矩阵由式的单元刚度矩阵，按单元的结点编号将单元刚度矩阵分为四个子块（阶数相同），逐块地将结点所对应的分为四个子块（阶数相同），逐块地将结点所对应的子块子块在结构的原

4、始刚度矩阵中对号入座，形成结构的原始刚度在结构的原始刚度矩阵中对号入座，形成结构的原始刚度矩阵。由于结点位移分量中包括了非自由结点的已知位移，矩阵。由于结点位移分量中包括了非自由结点的已知位移，原始刚度矩阵为奇异的，需进行边界条件处理，才能求解原始刚度矩阵为奇异的，需进行边界条件处理，才能求解自由结点位移。由于原始刚度矩阵的阶数较高，所以后处自由结点位移。由于原始刚度矩阵的阶数较高，所以后处理法的主要缺点是占用较多的计算机内存。理法的主要缺点是占用较多的计算机内存。二、总刚度矩阵的集成及约束处理二、总刚度矩阵的集成及约束处理烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟

5、台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 对于每个结点位移分量数相同的结构，原始刚度矩对于每个结点位移分量数相同的结构，原始刚度矩阵的阶数为结构的总结点数乘以结点位移分量的数目，阵的阶数为结构的总结点数乘以结点位移分量的数目，例如，每个结点位移分量数为例如，每个结点位移分量数为3 3的平面刚架，结构原始刚的平面刚架，结构原始刚度矩阵的阶数为度矩阵的阶数为3n3n。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 对于刚性支座，用划行划列法处理刚性支座，即直对于刚性支座，用划行划列法处理刚性支座，即

6、直接划去原始刚度方程中与零位移对应的行和列。这样做接划去原始刚度方程中与零位移对应的行和列。这样做有时要改变原方程的排列顺序有时要改变原方程的排列顺序,会给编程带来麻烦。为了会给编程带来麻烦。为了不改变原方程的排列顺序，同时又要引入边界条件不改变原方程的排列顺序，同时又要引入边界条件,采用采用“主一副零主一副零”法。法。（2）边界条件处理边界条件处理 设结点位移向量中第设结点位移向量中第r个位移等于零个位移等于零,即即 r=0，则在，则在结构的原始刚度矩阵结构的原始刚度矩阵k中的第中的第r行第行第r列中主对角元素列中主对角元素krr改改为为1其余元素改为零。同时将结点结点荷载列向量其余元素改为

7、零。同时将结点结点荷载列向量P中的第中的第r个分量也改为零。个分量也改为零。即即烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 对于支座位移等于给定值时，采用对于支座位移等于给定值时，采用“乘大数法乘大数法”。设。设结点位移向量中第结点位移向量中第r个位移等于个位移等于d0，在矩阵，在矩阵K与向量与向量P中，中，主对角元素主对角元素krr 改为改为Gkrr，将将Pr改为改为d0Gkrr，其中，其中G为一大为一大数通常取数通常取1081010。，烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟

8、台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 单元定位向量：单元定位向量：按单元连接结点编号顺序由结点未按单元连接结点编号顺序由结点未知位移编号组成的向量。知位移编号组成的向量。2.先处理法先处理法 (1)集成。集成。将单元刚度矩阵先按边界条件进行处理将单元刚度矩阵先按边界条件进行处理,然后按照单元连接结点的总位移编号将单元刚度矩阵的元然后按照单元连接结点的总位移编号将单元刚度矩阵的元素在结构的刚度矩阵中对号入座，形成总刚后即可进行求素在结构的刚度矩阵中对号入座，形成总刚后即可进行求解。上述过程可通过引入定位向量来实现。在单元定位向解。上述过程可通过引入定位向量来实现。在单元定位

9、向量中考虑边界条件，凡给定的结点位移分量，其位移总码量中考虑边界条件，凡给定的结点位移分量，其位移总码均编为零，与总码编为零相应的行、列元素在集成总刚时均编为零，与总码编为零相应的行、列元素在集成总刚时被屏弃在外。被屏弃在外。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 （2 2）边界条件处理。）边界条件处理。对于刚性支座，其位移总码均编对于刚性支座，其位移总码均编为零。对于支座位移等于给定值时，通常也将其位移总码为零。对于支座位移等于给定值时，通常也将其位移总码均编为零，将支座结点位移的影响转换成单元非结点荷

10、载均编为零，将支座结点位移的影响转换成单元非结点荷载,即，将支座结点位移转换成与该支座结点位移连接的各单即，将支座结点位移转换成与该支座结点位移连接的各单元在单元坐标系中的杆端位移，求出由此给定的杆端位移元在单元坐标系中的杆端位移，求出由此给定的杆端位移产生的单元固端力，然后转换成等效结点荷载。产生的单元固端力，然后转换成等效结点荷载。通常用主对角元素叠加法处理弹性支座。如果结构通常用主对角元素叠加法处理弹性支座。如果结构的第的第j个自由度是弹性约束，那么，把弹性支座的刚度系个自由度是弹性约束，那么，把弹性支座的刚度系数叠加到原始刚度矩阵主对角线的第数叠加到原始刚度矩阵主对角线的第j个元素上即

11、可得到个元素上即可得到经约束处理后的总刚度方程。经约束处理后的总刚度方程。3.弹性支座的处理弹性支座的处理烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 总刚度方程为整体结构的结点荷载与结点位移之间的总刚度方程为整体结构的结点荷载与结点位移之间的关系式，是结构应满足的平衡条件。无论何种结构，其总关系式，是结构应满足的平衡条件。无论何种结构，其总刚度方程都具有统一的形式刚度方程都具有统一的形式:4.总刚度方程和总刚度矩阵的性质与特点总刚度方程和总刚度矩阵的性质与特点 K=P 式中式中K为总刚度矩阵，为总刚度矩阵，为

12、结构的结点位移列向量，为结构的结点位移列向量，P为结点力列向量。为结点力列向量。总刚度矩阵总刚度矩阵K反应了整个结构的刚度，是描述结点反应了整个结构的刚度，是描述结点力与结点位移之间关系的系数矩阵。其矩阵的性质与力与结点位移之间关系的系数矩阵。其矩阵的性质与特点：特点：烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 (1)元素元素kij的物理意义为：当的物理意义为：当j=1=1而其他位移分量为而其他位移分量为零时产生在零时产生在i方向的杆端力。方向的杆端力。（2）主子块主子块Kii是由结点是由结点i的相关单元中与

13、结点的相关单元中与结点i相应的相应的主子块叠加而得。主子块叠加而得。（3）当当i、j为相关结点时，副子块为相关结点时，副子块Kij就等于连接就等于连接ij的杆的杆单元中相应的子块；若单元中相应的子块；若i、j不相关，则不相关，则Kij为零子块。为零子块。（4）总刚度矩阵为对称矩阵。总刚度矩阵为对称矩阵。（5）总刚度矩阵为稀疏带状矩阵。愈是大型结构，总刚度矩阵为稀疏带状矩阵。愈是大型结构，带状分布规律就愈明显。带状分布规律就愈明显。（6）总刚度矩阵主对角元素都大于零。通常是主对总刚度矩阵主对角元素都大于零。通常是主对角元素占优势的矩阵，因此，线形方程组的解有较好的稳角元素占优势的矩阵，因此，线形

14、方程组的解有较好的稳定性。定性。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助5.总刚度矩阵的最大半带宽总刚度矩阵的最大半带宽 总刚度矩阵的上三角部分，总刚度矩阵的上三角部分，从某行的主对角元素到该从某行的主对角元素到该行最末一个非零元素所具有的元素的个数称为该行的半带行最末一个非零元素所具有的元素的个数称为该行的半带宽。宽。各行半带宽的最大值称为总刚度矩阵的最大半带宽。各行半带宽的最大值称为总刚度矩阵的最大半带宽。对应于对应于后处理法后处理法，结构内部不存在组合结点时最大半，结构内部不存在组合结点时最大半带宽的

15、计算公式为带宽的计算公式为:d=(b+1)c，其中，其中b为单元为单元两端结点编码两端结点编码的最大差的最大差;c为结构中一个结点的位移分量数为结构中一个结点的位移分量数，显然，最大，显然，最大半带宽与结构的结点编码的顺序有关。通常应使相邻结点半带宽与结构的结点编码的顺序有关。通常应使相邻结点编码的最大差值为最小，即编码的最大差值为最小，即d 值为最小。值为最小。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 例如图示刚架，按图例如图示刚架，按图a 编码，编码，d=3(9+1)=30，而按，而按 b 图编码，图编

16、码，d=3(3+1)=12 。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 如果结构内部存在组合结点，并采用如果结构内部存在组合结点，并采用先处理法先处理法，则不，则不能用上述公式计算总刚度矩阵的最大半带宽，而应按照单能用上述公式计算总刚度矩阵的最大半带宽，而应按照单元编，利用单元定位向量求出总刚度矩阵的最大半带宽。元编，利用单元定位向量求出总刚度矩阵的最大半带宽。设用设用MAX表示单元表示单元(e)定位向量中的最大分量，定位向量中的最大分量，MIN表示表示单元单元(e)定位向量中的最小分量，则定位向量中的最小

17、分量，则d e=MAX-MIN+1总刚度矩阵的最大半带宽为：总刚度矩阵的最大半带宽为：d=MAX(d1d2dn)烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 （1）初学者易把单元的固端力与传统位移法中载常初学者易把单元的固端力与传统位移法中载常数混淆，造成求等效荷载时出错。单元的固端力是在固数混淆，造成求等效荷载时出错。单元的固端力是在固定单元的杆端其不能有任何位移时荷载作用下的杆端力定单元的杆端其不能有任何位移时荷载作用下的杆端力（即固端力）。（即固端力）。二、需要注意的几个问题二、需要注意的几个问题 （2）

18、在考虑轴向变形的单元刚度矩阵中剔除）在考虑轴向变形的单元刚度矩阵中剔除EA项，项，即得忽略轴向变形的单元刚度矩阵。即得忽略轴向变形的单元刚度矩阵。（3）为适应计算机计算、节省内存和机时，在对结）为适应计算机计算、节省内存和机时，在对结点编号时应力求使相关结点点编号时应力求使相关结点的最大差值为最小，以减小的最大差值为最小，以减小总刚度矩阵的带宽。总刚度矩阵的带宽。例如，对于梁式杆，不论连接该杆的结点是铰结点、例如，对于梁式杆，不论连接该杆的结点是铰结点、定向结点，均按两端固定梁计算固端力。定向结点，均按两端固定梁计算固端力。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟

19、台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 例例：图示梁用矩阵位移法求解时的基本未知量数图示梁用矩阵位移法求解时的基本未知量数目为多少？目为多少？解：解：基本未知量数目为基本未知量数目为2，即，即A点的竖向位移和转角。点的竖向位移和转角。三、例题三、例题烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助例：例：图示结构中单元图示结构中单元的定位向量为的定位向量为 。C.（0 0 1 3 2 4）T B.（2 3 4 0 0 1）T D.（3 2 4 0 0 1）T A.（0 0 1 2 3 4）T 解

20、：解：答案为答案为B。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助例：例：图示结构整体刚度矩阵图示结构整体刚度矩阵K中元素中元素k22等于（等于（）D.16EI/l A.28EI/3l B.12EI/l C.20EI/3l 解解：答案选答案选A。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 例：例：矩阵位移法中，结构的原始刚度方程是表示下矩阵位移法中，结构的原始刚度方程是表示下列两组量值之间的相互关系：（列两组量值之间的相互关系

21、：（）A杆端力与结点位移杆端力与结点位移 B杆端力与结点力杆端力与结点力 C结点力与结点位移结点力与结点位移 D结点位移与杆端力结点位移与杆端力 解：解：答案选答案选C。例：例：平面杆件结构用后处理法建立的原始刚度方程组平面杆件结构用后处理法建立的原始刚度方程组,（）A可求得全部结点位移可求得全部结点位移 B可求得可动结点的位移可求得可动结点的位移 C可求得支座结点位移可求得支座结点位移 D无法求得结点位移无法求得结点位移 解：解：答案选答案选D。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 例：例：图示结构若

22、考虑轴向变形，在未引入支撑条件时，图示结构若考虑轴向变形，在未引入支撑条件时，其整体刚度矩阵其整体刚度矩阵K是是_ 阶方阵。阶方阵。解解：答案为答案为2121。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 例：例：图示结构若只考虑弯曲变形，括号中的数字为结图示结构若只考虑弯曲变形，括号中的数字为结点位移分量编码，则其整体刚度矩阵中元素点位移分量编码，则其整体刚度矩阵中元素k11等于（等于（）.A.B.C.D.解：解：答案选答案选D。提示：提示：在不考率轴向变形时，在不考率轴向变形时，结点结点2和结点和结点3只有

23、水只有水平位移和转角，杆件平位移和转角，杆件12对对k11的贡献为的贡献为12（2EI）l 3，杆件杆件34对对k11的贡献为的贡献为12EI（l2）3。烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助 例例：用矩阵位移法计算图用矩阵位移法计算图a所示连续梁所示连续梁,并画并画M图，图，EI=常数。常数。q=12kN/m,l=6m。解：解：(1)建立坐标系、对单元和结点编号如图建立坐标系、对单元和结点编号如图b b，单，单元刚度矩阵元刚度矩阵 单元定位向量单元定位向量=（0 1）T，=（1 2）T，=（2 0）T

24、(2)将各单元刚度矩阵中的元素按单元定位向量在将各单元刚度矩阵中的元素按单元定位向量在K中中对号入座，得整体刚度矩阵对号入座，得整体刚度矩阵 烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助(3)连续梁的等效结点荷载连续梁的等效结点荷载(4)将整体刚度矩阵将整体刚度矩阵K和等效结点荷载和等效结点荷载P代入基本方程得代入基本方程得(5)解方程得解方程得 (6)求杆端力并绘制弯矩图如图所示求杆端力并绘制弯矩图如图所示c。M图（图（kNm）烟台大学烟台大学第四章第四章 静定结构的位移计算静定结构的位移计算 烟台大学烟台大学第第8章章 矩阵位移法矩阵位移法返回开篇帮助四、思考题四、思考题 1.单元刚度矩阵的物理意义及其性质与特点各是什么？单元刚度矩阵的物理意义及其性质与特点各是什么？2.单元定位向量是由什么组成？他的用处是什么？单元定位向量是由什么组成？他的用处是什么？3.刚架中有铰结点时应该怎样处理？刚架中有铰结点时应该怎样处理？