第13章矩阵位移法 学习目的和要求 矩阵位移法是以计算机为计算工具的现代化结构分析方法。基于该法的结构分 析程序在结构设计中得到了广泛的应用。因此,以计算机进行结构分析是本章的学 习目的。 矩阵位移法是以位移法为理论基础,以矩阵为表现形式,以计算机为为运算工 具的综合分析方法。引入矩阵运算的目的是式计算过程程序化,便于计算机自动化 处理尽管矩阵位移法运算模式呆板,过程繁杂,但这些正是计算机所需要的和十分 容易解决的。矩阵位移法的特点是用“机算”代替“手算”。因此,学习本章是既 要了解它与位移法的共同点,更要了解它的一些新手法和新思想。 本章的基本要求: 1.矩阵位移法包含两个基本环节:单元分析和整体分析。 2.在单元分析中,熟练掌握单元刚度矩阵和单元等效荷载的概念和 形成。熟练掌握已知结点位移求单元杆端力的计算方法。 3.在整体分析中,熟练掌握结构整体刚度矩阵元素的物理意义和集 成过程,熟练掌握结构综合结点荷载的集成过程。掌握单元 定位向量的建立,支撑条件的处理 1.自由式单元的单元刚度矩阵不要求背记,但要领会其物理意义 并会有它推出特殊单元的单元刚度矩阵。 学习内容 有限单元法的基本概念,结构离散化。 平面杆系结构的单元分析:局部坐标系下的单元刚度矩阵和整体坐标系下的单 元刚度矩阵 平面杆系结构的整体分析:结构整体刚度矩阵和结构整体刚度方程。 边界条件的处理,单元内力计算 矩阵位移法的计算步骤和应用举例
第 13 章 矩阵位移法 学习目的和要求 矩阵位移法是以计算机为计算工具的现代化结构分析方法。基于该法的结构分 析程序在结构设计中得到了广泛的应用。因此，以计算机进行结构分析是本章的学 习目的。 矩阵位移法是以位移法为理论基础，以矩阵为表现形式，以计算机为为运算工 具的综合分析方法。引入矩阵运算的目的是式计算过程程序化，便于计算机自动化 处理尽管矩阵位移法运算模式呆板，过程繁杂，但这些正是计算机所需要的和十分 容易解决的。矩阵位移法的特点是用“机算”代替“手算”。因此，学习本章是既 要了解它与位移法的共同点，更要了解它的一些新手法和新思想。 本章的基本要求： 1. 矩阵位移法包含两个基本环节：单元分析和整体分析。 2. 在单元分析中，熟练掌握单元刚度矩阵和单元等效荷载的概念和 形成。熟练掌握已知结点位移求单元杆端力的计算方法。 3. 在整体分析中，熟练掌握结构整体刚度矩阵元素的物理意义和集 成过程，熟练掌握结构综合结点荷载的集成过程。掌握单元 定位向量的建立，支撑条件的处理。 1. 自由式单元的单元刚度矩阵不要求背记，但要领会其物理意义， 并会有它推出特殊单元的单元刚度矩阵。 学习内容 有限单元法的基本概念，结构离散化。 平面杆系结构的单元分析：局部坐标系下的单元刚度矩阵和整体坐标系下的单 元刚度矩阵。 平面杆系结构的整体分析：结构整体刚度矩阵和结构整体刚度方程。 边界条件的处理，单元内力计算。 矩阵位移法的计算步骤和应用举例
§13.1矩阵位移法概述 1、矩阵位移法的基本思路 矩阵位移法是以位移法作为理论基础,以矩阵作为数学表达形式,以电子计算机作为计算 工具三位一体的分析方法。引入矩阵运算,使得公式排列紧凑,运算形式统一,便于计算过程 程序化,适宜于计算机进行自动化处理 矩阵位移法包含两个基本环节:单元分析和整体分析。 先将结构离散成有限个单元,按照单元的力学性质(物理关系),建立单元刚度方程,形 成单元刚度矩阵:然后在满足变形条件和平衡条件的前提下,将这些单元集合成整体,即由单 元刚度矩阵集成整体刚度矩阵,建立结构的位移法基本方程,进而求出结构的位移和内力。这 样,在一撤一搭的过程中,使一个复杂结构的计算问题转化为有限个简单单元的分析与集成问 2、单元划分 在杆件结构矩 11111111111I11 分析中,一般是把杆件 ② 的转折点、汇交点、边 界点、突变点或集中荷 载作用点等列为结点, 结点之间的杆件部分 作为单元。如图1(a) 所示。为了减少基本未 知量的数目,跨间集中 荷载作用点可不作为 结点,但要计算跨间荷 载的等效结点荷载:跨 间结点也可不作为结 点(如图1(b)所示), 但要推导相应的单元 刚度矩阵,编程序麻
§13.1 矩阵位移法概述 1、矩阵位移法的基本思路 矩阵位移法是以位移法作为理论基础，以矩阵作为数学表达形式，以电子计算机作为计算 工具三位一体的分析方法。引入矩阵运算，使得公式排列紧凑，运算形式统一，便于计算过程 程序化，适宜于计算机进行自动化处理。 矩阵位移法包含两个基本环节：单元分析和整体分析。 先将结构离散成有限个单元，按照单元的力学性质（物理关系），建立单元刚度方程，形 成单元刚度矩阵；然后在满足变形条件和平衡条件的前提下，将这些单元集合成整体，即由单 元刚度矩阵集成整体刚度矩阵，建立结构的位移法基本方程，进而求出结构的位移和内力。这 样，在一撤一搭的过程中，使一个复杂结构的计算问题转化为有限个简单单元的分析与集成问 题。 2、单元划分 在杆件结构矩阵 分析中，一般是把杆件 的转折点、汇交点、边 界点、突变点或集中荷 载作用点等列为结点， 结点之间的杆件部分 作为单元。如图 1（a） 所示。为了减少基本未 知量的数目，跨间集中 荷载作用点可不作为 结点，但要计算跨间荷 载的等效结点荷载；跨 间结点也可不作为结 点（如图 1（b）所示）， 但要推导相应的单元 刚度矩阵，编程序麻
§13.2单元分析——局部坐标系下 单元分析的目的是建立单元刚度方程,形成单元刚度矩阵 1、坐标系的选择 在矩阵位移法中采用两种坐标系:局部坐标 系和整体坐标系。 采用局部坐标系(以杆的轴线作为x轴如图 2),可直接由虎克定律、转角位移方程得到单元 刚度方程,导出的单元刚度矩阵具有最简单的形 取超静定的基本体系 在局部坐标系中,杆端力及杆端位移的正方 如图2所示。单元刚度方程可表示 4411 2 3 其中单元的杆端力列阵和杆端位移列阵为 可用力矩分配法画球61由载常数表画球求4p -{x,FM,兄,F1M 2、局部坐标系中的单元刚度矩阵: EA 0126E 12E6E bE 4E1 6El 2E 单元刚度矩阵为 EA EA 12E/ 6El 13-12 3 6El 2E1 6E/ 4E1 (例子112)
烦。 §13.2 单元分析——局部坐标系下 单元分析的目的是建立单元刚度方程，形成单元刚度矩阵。 1、坐标系的选择: 在矩阵位移法中采用两种坐标系：局部坐标 系和整体坐标系。 采用局部坐标系（以杆的轴线作为 轴如图 2），可直接由虎克定律、转角位移方程得到单元 刚度方程，导出的单元刚度矩阵具有最简单的形 式。 在局部坐标系中，杆端力及杆端位移的正方 向如图 2 所示。单元刚度方程可表示 为： 其中单元的杆端力列阵和杆端位移列阵为： ↓↓↓↓↓↓↓↓↓↓↓↓↓↓ q X1 可用力矩分配法画M求δ11 ↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓ q a a a 由载常数表画MP求Δ1P X1=1 ↓↓↓↓↓↓↓↓↓↓↓↓↓↓ q MP M 取超静定的基本体系 ↓↓↓↓↓↓↓↓↓↓↓↓↓↓ q X1 可用力矩分配法画M求δ11 ↓↓↓↓↓↓↓↓↓↓↓↓↓↓↓ q a a a 由载常数表画MP求Δ1P X1=1 ↓↓↓↓↓↓↓↓↓↓↓↓↓↓ q MP M 取超静定的基本体系 2、局部坐标系中的单元刚度矩阵： 单元刚度矩阵为： （例子 112）
3、单元刚度矩阵的特性 (1)单元刚度系数的意义:单刚中的每个元素称为单元刚度系数,代表由于单位杆端位移引 起的杆端力 如第讠行第j列元素代表当第j个杆端位移分量=1(其它位移分量为零)时引起的第i个杆端力 分量的值。 单刚中第j列元素代表当第j个杆端位移分量=1(其它位移分量为零)时引起的六个杆端力分量 的值 由图10-4可见,巧=1产生的单元变形及单元的杆端力与 =-1 生的单元变形及单元的杆端力相 由此得到:单元刚度矩阵的第二列元素变符号即第五列可= 元素,第一列元素变符号即第四列元素。第三列元素不变 符号即第六列元素,但要注意k36=k32,k6=2k。由 于单元刚度矩 阵是对称矩阵,所以,各行元素之间也具有类似的关系。 (2)单元刚度矩阵是对称矩阵:由反力互等定理可知,单元刚度矩阵是对称矩阵 (3)一般单元的单元刚度矩阵是奇异矩阵,不存在逆阵。因此上 由单元刚度 方程,如已知杆端位移可求出杆端力,且是唯一解。但如已知杄端力,则求不岀杆端位移,杆端 位移可能无解,可能无唯一解。 (4)可按杆端将单元刚度方程写成分块形式 〔(°[]° 4、特殊单元的单元刚度矩阵 忽略轴向变形时梁单元在局部坐标系中的单元刚度矩阵。 由k△+F=,求出△1 滑四“由M丽△+M叠加最后弯矩图。 力矩分配斗画求五 对称问要力分配烈 续梁单元的单元刚度矩阵 对称问题 反对称问 按位移法 题按力法 或力矩分 或无剪切 配法计算 分配法求
↓↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m ↓↓↓↓↓↓↓↓↓↓↓↓↓ F1P Δ 20kN/m ↓↓↓↓↓↓↓↓↓↓↓↓↓ F1P 20kN/m ↓↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m Δ=1 k11 由力矩分配法画MP求F1P 取无侧移的结构为位移法基本体系 由力矩分配法画MP求k11 ↓↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m ↓↓↓↓↓↓↓↓↓↓↓↓↓ F1P Δ 20kN/m ↓↓↓↓↓↓↓↓↓↓↓↓↓ F1P 20kN/m ↓↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m Δ=1 k11 由力矩分配法画MP求F1P 取无侧移的结构为位移法基本体系 由力矩分配法画MP求k11 P P/2 P/2 P/2 P/2 对称问题 反对称问题 P/2 P/2 对称问题 按位移法 或力矩分 配法计算 反对称问 题按力法 或无剪切 分配法求 P P/2 P/2 P/2 P/2 对称问题 反对称问题 P/2 P/2 对称问题 按位移法 或力矩分 配法计算 反对称问 题按力法 或无剪切 分配法求 3、单元刚度矩阵的特性 （1）单元刚度系数的意义： 单刚中的每个元素称为单元刚度系数，代表由于单位杆端位移引 起的杆端力。 如第 i 行第 j 列元素代表当第 j 个杆端位移分量=1（其它位移分量为零）时引起的第 i 个杆端力 分量的值。 单刚中第 j 列元素代表当第 j 个杆端位移分量=1（其它位移分量为零）时引起的六个杆端力分量 的值。 由图 10-4 可见， 产生的单元变形及单元的杆端力与 产生的单元变形及单元的杆端力相同。 由此得到：单元刚度矩阵的第二列元素变符号即第五列 元素，第一列元素变符号即第四列元素。第三列元素不变 符号即第六列元素，但要注意 , 。由 于单元刚度矩 阵是对称矩阵，所以，各行元素之间也具有类似的关系。 （2）单元刚度矩阵是对称矩阵：由反力互等定理可知， 由δ11X1+Δ1P=0，求出 X1， 由 M=M1X1+ MP叠加最后弯矩图。 单元刚度矩阵是对称矩阵。 （3）一般单元的单元刚度矩阵是奇异矩阵，不存在逆阵。因此上， 由单元刚度 方程，如已知杆端位移可求出杆端力，且是唯一解。但如已知杆端力，则求不出杆端位移，杆端 位移可能无解，可能无唯一解。 （4）可按杆端将单元刚度方程写成分块形式： 4、特殊单元的单元刚度矩阵 忽略轴向变形时梁单元在局部坐标系中的单元刚度矩阵。 由 k11Δ1+F1P=0，求出Δ1， 由 M=M1Δ1+ MP叠加最后弯矩图。 连续梁单元的单元刚度矩阵
30kN M图 位移法求R l+/法计第+ (KN.M) 3gl/8=30M 桁架单元在局部坐标系中的单元刚度矩阵 (刘老师:下边的图都不好了) §133单元分析一一整体坐标系下 选局部坐标系推导单元刚度矩阵方便且单元刚度矩阵的形式简单。但是,在一个复杂的结构中, 各单元的局部坐标系不尽相同,很不统一。为了进行整体分析,必须选一个统一的坐标系(称为整体 坐标系)。按这个统一的坐标系来建立各单元的刚度矩阵 1、单元坐标转换矩阵 (例子113)
i i 8m EI=∞ A B C D ↓↓↓↓↓↓↓↓↓ 10kN/m i i EI=∞ A B C D ↓↓↓↓↓↓↓↓↓ 10kN/m R 3ql/8=30kN R=30kN 4m 4m 30kNR 80 6 48 24 96 96 M图 (kN.M) 128 80 96 96 位移法求R 用剪力分 i 配法计算 i 8m EI=∞ A B C D ↓↓↓↓↓↓↓↓↓ 10kN/m i i 8m EI=∞ A B C D i i 8m EI=∞ i i 8m EI=∞ A B C D ↓↓↓↓↓↓↓↓↓ 10kN/m i i EI=∞ A B C D ↓↓↓↓↓↓↓↓↓ 10kN/m i i EI=∞ A B C D i i EI=∞ A B C D ↓↓↓↓↓↓↓↓↓ 10kN/m R 3ql/8=30kN R=30kN 4m 4m 30kNR 80 6 48 24 96 96 M图 (kN.M) 128 80 96 96 位移法求R 用剪力分 配法计算 桁架单元在局部坐标系中的单元刚度矩阵。 （刘老师：下边的 图都不好了） §13.3 单元分析——整体坐标系下 选局部坐标系推导单元刚度矩阵方便且单元刚度矩阵的形式简单。但是，在一个复杂的结构中， 各单元的局部坐标系不尽相同，很不统一。为了进行整体分析，必须选一个统一的坐标系（称为整体 坐标系)。按这个统一的坐标系来建立各 单元的刚度矩阵。 1、单元坐标转换矩阵： （例子 113）
整体坐标系如图1中的axy坐标系。 两种坐标系中单元的杆端力和杆端位移转换 关系 k-154 其中[刀称为单元坐标转换矩阵 cos a sin a 0 0 00 +∞|28E in a cos a 00 0 0 00 00 cos a sin a 0 0 0-sin a cos a 0 单元坐标转换矩阵[刁是正交矩阵 2、整体坐标系中的单元刚度矩阵:(例子114) 整体坐标系中的单元刚度矩阵的特性 将(F)=[7]F°,(△)0=[7]△ 整体坐标系中的单元刚度矩阵与局部坐标系中 代入(F}8=[k(△ 得:[z](P)·=[][7]4) 的单元刚度矩阵有类似的特性。 上式两边乘[r 另外还有,局部坐标系中的单元刚度矩阵,只 [T][](F)=[][k][](△) 与单元的几何形状、物理常数有关,而与单元的 (F)=[T][k][](△ 位置和方位无关。整体坐标系中的单元刚度矩 (F)=[k](△) k=[T]k][r] 阵,与单元的几何形状、物理常数及单元的方位 无有关 §13.4连续梁整体分析 整体分析的目的是建立整体刚度方程,导出整体刚度矩阵。具体作法有两种:一种是传统位 移法,另一种是单元集成法(即刚度集成法或直接刚度法)。下面以连续梁为例,用传统位移法 建立整体刚度方程,进而总结出单元集成法。 1、连续梁的整体刚度矩阵:(例子115 整体刚度方程是整体结构的结点力与结点位移之间的关系式,是通过考虑结构的变形连续条件和 平衡条件建立起来的。如右图结构有:
整体坐标系如图 1 中的 oxy 坐标系。 两种坐标系中单元的杆端力和杆端位移转换 关系： Pl/4 P P P P P l P 2 2 ( l) EI Pcr = μ=1 μ=1/2 l l/2 l/2 l EI EI EI EI Pl 48 3 3Pl/40 3Pl/40 EI Pl 960 11 3 7Pl/40 多余约束的存在，使结构的强度、刚度、稳定性都有所提高。 Pl/4 P P P P P l P 2 2 ( l) EI Pcr = μ=1 μ=1/2 l l/2 l/2 l EI EI EI l l/2 l/2 l EI EI EI EI Pl 48 3 EI Pl 48 3 3Pl/40 3Pl/40 EI Pl 960 11 3 EI Pl 960 11 3 7Pl/40 多余约束的存在，使结构的强度、刚度、稳定性都有所提高。 其中[T]称为单元坐标转换矩阵 单元坐标转换矩阵[T]是正交矩阵， 8m 6m ↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m I1 I2 I2 I1=2I2 53.3 53.3 106.7 ↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m I1 I2 I2 I1>>I2 ≈0 160 ≈0 ↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m I1 I2 I2 I1>I2 ≈0 160 ≈0 ↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m I1 I2 I2 ↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m I1 I2 I2 I1<<I2 106.7 106.7 53.3 ↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m I1 I2 I2 ↓↓↓↓↓↓↓↓↓↓↓↓ 20kN/m I1 I2 I2 I1=1.5I2 80 80 80 －35° －35° －35° +15° +15° +15° 40cm 60cm 8m 6m 94.2 N=－15.74 M&N×αEI －35° －35° －35° －35° －35° －35° +15° +15° +15° +15° +15° +15° 40cm 60cm 8m 6m 94.2 N=－15.74 M&N×αEI 2、整体坐标系中的单元刚度矩阵： （例子 114） 整体坐标系中的单元刚度矩阵的特性： 整体坐标系中的单元刚度矩阵与局部坐标系中 的单元刚度矩阵有类似的特性。 另外还有，局部坐标系中的单元刚度矩阵 ，只 与单元的几何形状、物理常数有关，而与单元的 位置和方位无关。整体坐标系中的单元刚度矩 阵，与单元的几何形状、物理常数及单元的方位 无有关 §13.4 连续梁整体分析 整体分析的目的是建立整体刚度方程，导出整体刚度矩阵。具体作法有两种：一种是传统位 移法，另一种是单元集成法（即刚度集成法或直接刚度法）。下面以连续梁为例，用传统位移法 建立整体刚度方程，进而总结出单元集成法。 1、 连续梁的整体刚度矩阵：（例子 115） 整体刚度方程是整体结构的结点力与结点位移之间的关系式，是通过考虑结构的变形连续条件和 平衡条件建立起来的。如右图结构有：
F1=K11△1+K12△2+K13△3 F2=K21A1+k22△2+K2343 F=k3141+k3242+k3343 {F}=({4 38 K11K1K [xFK21K22 K2 任何结构的整体刚度方程都具4, 有统一的形式: 4X [K△)={以 []是整体刚度矩阵,{4}结构的结点位移列向量,{结构的结点力列向量。 2、整体刚度矩阵的性质 ①[中的元素K称为整体刚度系数,它表示当第j个结点位移分量A=1(其他结点位移 分量为零)时所产生的第i个结点力F ②[]是对称矩阵,是稀疏带状矩阵。 ③引入支承条件之前是奇异矩阵,引入支承条件之后是非奇异矩阵,存在逆阵 3、直接刚度法:(例子116) 由变形连续条件知,结点发生单位位移,交与该结点的各单元的杆端也发生单位位移:由刚度系 数的物理意义知,单位杆端位移产生的杆端力是单元刚度矩阵中的元素,单位结点位移产生的结 点力是整体刚度矩阵中的元素:由平衡条件知交与某结点的各单元杆端力之和等于该结点的相应 结点力。故整体刚度矩阵中的元素是由对应的单元刚度矩阵中的元素叠加而成 综上所述,直接刚度法是根据单元的结点位移分量的局部码和总码之间的对应关系,由单元 刚度矩阵集成结构整体刚度矩阵。 8)9宽 KI 2 下一
Δ1 Δ2 Δ3 F1 F2 F3 Δ1 Δ2 Δ3 F1 F2 F3 Δ1=1 Δ1× K11 K21 K31 K12 K22 Δ K32 2=1 Δ2× K13 K23 K33 Δ3=1 Δ3× Δ1 Δ2 Δ3 F1 F2 F3 Δ1 Δ2 Δ3 1 Δ2 Δ3 F1 F2 F3 F1 F2 F3 Δ1=1 Δ1× K11 K21 K31 K11 K21 K31 K12 K22 K32 K12 K22 Δ K32 2=1 Δ2× K13 K23 K33 K13 K23 K33 Δ3=1 Δ3× Δ2=1 K12 K22 K32 1 2 =1 1 1 2 k12 1 k22 1 2 1 =1 2 k11 2 k21 K 2 22 k22 1 k11 2 k21 2 K32 Δ2=1 K12 K22 Δ K32 2=1 K12 K22 K32 K12 K22 K32 1 2 =1 1 1 2 k12 1 k22 1 2 1 =1 2 k11 2 k21 K 2 22 k22 1 k11 2 K22 k22 1 k11 2 k21 2 K32 3 31 1 32 2 33 3 2 21 1 22 2 23 3 1 11 1 12 2 13 3 = + + = + + = + + F K K K F K K K F K K K F = K Δ = 31 32 33 21 22 23 11 12 13 K K K K K K K K K K 3 31 1 32 2 33 3 2 21 1 22 2 23 3 1 11 1 12 2 13 3 = + + = + + = + + F K K K F K K K F K K K = 31 32 33 21 22 23 11 12 13 K K K K K K K K K K 任何结构的整体刚度方程都具 有统一的形式： [K]是整体刚度矩阵，{Δ}结构的结点位移列向量，{F}结构的结点力列向量。 2、整体刚度矩阵的性质 ： ①[K]中的元素 Kij称为整体刚度系数，它表示当第 j 个结点位移分量Δj=1（其他结点位移 分量为零）时所产生的第 i 个结点力 Fi。 ②[K]是对称矩阵，是稀疏带状矩阵。 ③引入支承条件之前是奇异矩阵，引入支承条件之后是非奇异矩阵，存在逆阵。 3、直接刚度法：（例子 116） 由变形连续条件知，结点发生单位位移，交与该结点的各单元的杆端也发生单位位移；由刚度系 数的物理意义知，单位杆端位移产生的杆端力是单元刚度矩阵中的元素，单位结点位移产生的结 点力是整体刚度矩阵中的元素；由平衡条件知交与某结点的各单元杆端力之和等于该结点的相应 结点力。故整体刚度矩阵中的元素是由对应的单元刚度矩阵中的元素叠加而成。 综上所述，直接刚度法是根据单元的结点位移分量的局部码和总码之间的对应关系，由单元 刚度矩阵集成结构整体刚度矩阵
4、单元定位向量: 集成整体刚度矩阵的关键,是确定单元刚度矩阵中的元素在整体刚度矩阵中的位置 首先要知道单元的结点位移分量的局部码和总码之间的对应关系,即单元定位向量。单元定位向 量是单元结点位移总码按局部码顺序排列而成的向量记为{A}。下图所示连续梁各单元的这种 对应关系及各单元定位向量如下表 BR (1) (2)(1 (2) 单元 局部码→总码 单元定位向量{4A 2 (1) 2 (2) 其次,要注意在单元刚度矩阵中,元素按局部码排列,在整体刚度矩阵中,元素按总码排列。所 以单元刚度矩阵中的元素在整体刚度矩阵中的定位原则是:将各单元的单刚的行列局部码(i)、 (j)换成对应的结点位移总码Ai、λj,按此行列总码将单刚元素送入总刚。即:→Ex2 直接刚度法的实施过程如下 ①将[置 ②将单元定位向量写在单元刚度矩阵的上方和右侧,确定出单元刚度矩阵 中的各元素在整体刚度矩阵[中的 位置并累加到[M ③对所有单元循环一边,最后得到整体刚度矩阵[ 2、单元定位向量 1)结点位移分量的统一编码一一总码 平面刚架中的一个结点可能有一个、两个或三个结点位移,在进行结点位移分量编码时,应考虑 每个结点的位移情况,对结构的所有结点位移分量进行统一编码。对每个结点的三个位移分量
2 1 1 2 2 1 3 A B C (1) (2) (1) (2) 2 1 1 2 2 1 3 A B C (1) (2) (1) (2) 单元 局部码→总码 单元定位向量{λ} e 1 2 θA （1） → 1 θB （2） → 2 {λ} = 1 1 2 θB （1） → 2 θC （2） → 3 {λ} = 2 2 3 单元 局部码→总码 单元定位向量{λ} e 1 2 θA （1） → 1 θB （2） → 2 {λ} = 1 1 2 {λ} = 1 1 2 θB （1） → 2 θC （2） → 3 {λ} = 2 2 3 {λ} = 2 2 3 4、单元定位向量 ： 集成整体刚度矩阵的关键，是确定单元刚度矩阵中的元素在整体刚度矩阵中的位置。 首先要知道单元的结点位移分量的局部码和总码之间的对应关系，即单元定位向量。单元定位向 量是单元结点位移总码按局部码顺序排列而成的向量记为{λ}。 下图所示连续梁各单元的这种 对应关系及各单元定位向量如下表。 其次，要注意在单元刚度矩阵中，元素按局部码排列，在整体刚度矩阵中，元素按总码排列。所 以单元刚度矩阵中的元素在整体刚度矩阵中的定位原则是：将各单元的单刚的行列局部码（i）、 （j）换成对应的结点位移总码λi、 λj，按此行列总码 将单刚元素送入总刚。即： 直接刚度法的实施过程如下： ①将[K]置零。 ②将单元定位向量 写在单元刚度矩阵 的上方和右侧，确定出单元刚度矩阵 中的各元素在整体刚度矩阵[K]中的 位置并累加到[K]。 ③对所有单元循环一边，最后得到整体刚度矩阵[K]。 2、单元定位向量 1）结点位移分量的统一编码——总码 平面刚架中的一个结点可能有一个、两个或三个结点位移，在进行结点位移分量编码时，应考虑 每个结点的位移情况，对结构的所有结点位移分量进行统一编码。对每个结点的三个位移分量
按照先x轴方向,再y轴方向后转动方向的顺序依次编码,编完一个结点再编下一个结点。对 已知为零的结点位移分量,其总码均编为0。如图(a) 忽略阻尼的振动规律 考虑阻尼的振动规律 结构的自振频率是结构的固有特性,与外因无关。 简谐荷载作用下有可能出现共振。 自由振动的振幅永不衰减。自由振动的振幅逐渐衰减。 共振时的振幅趋于无穷大。共振时的振幅较大但为有限值。 2)单元定位向量 图(a)所示各杆轴上的箭头表示各单元局部坐标系的轴的正方向,单元在始末两端的六个 位移分量的局部码(1)、(2)、…(6)也是按着先x轴方向,再y轴方向后转动放顺序依次编码的, 如图(b)所示。 单元定位向量仍然是由单元结点位移总码按局部码顺序排列而成的向量。各单元结 点位移分量局部码和总码之间的对应关系及单元定位向量如下表 单元① 单元② 局部码→总码 单元定位向量 局部码→总码单元定位向量 (2)→2 (2)→2 (3)→3 (3)→3 (4)→0 (4)→0 (5)→0 (5)→0 (6)→4 (6)→0 3、整刚集成过程:与连续梁整体刚度矩阵集成过程相同。(例子117) ①将[置零 ②将单元定位向量写在单元刚度矩阵的上方和右侧,确定出单元刚度矩阵中的各元素在整体 刚度矩阵[和中的位置并累加到[]。 ③对所有单元循环一边,最后得到整体刚度矩阵[和
忽略阻尼的振动规律 考虑阻尼的振动规律 结构的自振频率是结构的固有特性，与外因无关。 简谐荷载作用下有可能出现共振。 自由振动的振幅永不衰减。 自由振动的振幅逐渐衰减。 共振时的振幅趋于无穷大。 共振时的振幅较大但为有限值。 忽略阻尼的振动规律 考虑阻尼的振动规律 结构的自振频率是结构的固有特性，与外因无关。 简谐荷载作用下有可能出现共振。 自由振动的振幅永不衰减。 自由振动的振幅逐渐衰减。 共振时的振幅趋于无穷大。 共振时的振幅较大但为有限值。 按照先 x 轴方向，再 y 轴方向后转动 方向的顺序依次编码，编完一个结点再编下一个结点。对 于已知为零的结点位移分量，其总码均编为 0。如图（a）。 2）单元定位向量 图（a）所示各杆轴上的箭头表示各单元局部坐标系的 轴的正方向，单元在始末两端的六个 位移分量的局部码(1)、(2)、…(6)也是按着先 x 轴方向，再 y 轴方向后转动放顺序依次编码的， 如图(b)所示。 单元定位向量仍然是由单元结点位移总码按局部码顺序排列而成的向量。各单元结 点位移分量局部 码和总码之间的对应关系及单元定位向量如下表。 单元① 单元② 局部码→总码 单元定位向量 局部码→总码 单元定位向量 (1)→1 (2)→2 (3)→3 (4)→0 (5)→0 (6)→4 (1)→1 (2)→2 (3)→3 (4)→0 (5)→0 (6)→0 3、整刚集成过程：与连续梁整体刚度矩阵集成过程相同。（例子 117） ①将[K]置零。 ②将单元定位向量写在单元刚度矩阵的上方和右侧，确定出单元刚度矩阵中的各元素在整体 刚度矩阵[K]中的位置并累加到[K]。 ③对所有单元循环一边，最后得到整体刚度矩阵[K]
4、铰结点的处理:(例子118) 给结点位移分量进行统一编码时,考虑铰结点的特点。图11所示结构在铰结点处的两杆端结点 应看作半独立的两个结点(C1和C2)它们的线位移相同,编成同码,角位移不同,编成异码。 §13.6等效结点荷载 1、位移法基本方程: 1)整体刚度矩阵 前面讨论了结构的整体刚度矩阵,建立了整体刚度方程:F=[{△} 它表示由结点位移{△}推算结点力{F}的关系式。它只反映了结构的刚度性质,不涉及结构上的实际荷 并不是用以分析原结构的位移法方程 2)位移法基本方程 按位移法计算,就是将原结构分解为如图示位移法基本体系的两种状态: F3 原结构 状态 状态二 状态一只有荷载单独作用(结点位移为零)一一此时在基本结构中引起的结点约束力,记为{FP}。 状态二只有结点位移单独作用(荷载为零)一一此时在基本结构中引起的结点约束力为 (b)
4、铰结点的处理: （例子 118） 给结点位移分量进行统一编码时，考虑铰结点的特点。图 11 所示结构在铰结点处的两杆端结点 应看作半独立的两个结点（C1 和 C2）它们的线位移相同，编成同码，角位移不同，编成异码。 §13.6 等效结点荷载 1、位移法基本方程 ： 1）整体刚度矩阵 前面讨论了结构的整体刚度矩阵，建立了整体刚度方程：F}=[K]{Δ} ……（a） 它表示由结点位移{Δ}推算结点力{F}的关系式。它只反映了结构的刚度性质，不涉及结构上的实际荷载。 并不是用以分析原结构的位移法方程。 2）位移法基本方程 按位移法计算，就是将原结构分解为如图示位移法基本体系的两种状态： ↓↓↓↓↓↓↓↓↓ 原结构 F1 x y ① ② ③ ② ③ ③ ↓↓↓↓↓↓↓↓↓ ① ② ③ ② ③ ③ ① ② ③ ② ③ ③ F2 F3 F6 F5 F4 －F1 －F3 －F2 －F4 －F6 －F5 = + 状态一 状态二 状态一只有荷载单独作用（结点位移为零）——此时在基本结构中引起的结点约束力，记为{FP }。 状态二只有结点位移单独作用（荷载为零）——此时在基本结构中引起的结点约束力为 ……(b)