结构力学矩阵位移法演示文稿

结构力学矩阵位移法演示文稿现在是1页\一共有48页\编辑于星期四结构力学矩阵位移法现在是2页\一共有48页\编辑于星期四学习内容

有限单元法的基本概念，结构离散化。

平面杆系结构的单元分析：局部坐标系下的单元刚度矩阵和整体坐标系下的单元刚度矩阵。

平面杆系结构的整体分析：结构整体刚度矩阵和结构整体刚度方程。

边界条件的处理，单元内力计算。

利用对称性简化位移法计算。

矩阵位移法的计算步骤和应用举例。3现在是3页\一共有48页\编辑于星期四学习目的和要求

目的：矩阵位移法是以计算机为计算工具的现代化结构分析方法。基于该法的结构分析程序在结构设计中得到了广泛的应用。因此，以计算机进行结构分析是本章的学习目的。

矩阵位移法是以位移法为理论基础，以矩阵为表现形式，以计算机为运算工具的综合分析方法。引入矩阵运算的目的是使计算过程程序化，便于计算机自动化处理。尽管矩阵位移法运算模式呆板，过程繁杂，但这些正是计算机所需要的和十分容易解决的。矩阵位移法的特点是用“机算”代替“手算”。因此，学习本章是既要了解它与位移法的共同点，更要了解它的一些新手法和新思想。4现在是4页\一共有48页\编辑于星期四学习目的和要求

要求：矩阵位移法包含两个基本环节：单元分析和整体分析。

在单元分析中，熟练掌握单元刚度矩阵和单元等效荷载的概念和形成。熟练掌握已知结点位移求单元杆端力的计算方法。

在整体分析中，熟练掌握结构整体刚度矩阵元素的物理意义和集成过程，熟练掌握结构综合结点荷载的集成过程。掌握单元定位向量的建立，支撑条件的处理。

自由式单元的单元刚度矩阵不要求背记，但要领会其物理意义，并会有它推出特殊单元的单元刚度矩阵。5现在是5页\一共有48页\编辑于星期四矩阵位移法以传统的结构力学作为理论基础;

以矩阵作为数学表达形式;

以电子计算机作为计算手段三位一体的解决各种杆系结构受力、变形等计算的方法。第一节矩阵位移法概述采用矩阵进行运算，不仅公式紧凑，而且形式统一，便于使计算过程规格化和程序化。这些正是适应了电子计算机进行自动化计算的要求。6现在是6页\一共有48页\编辑于星期四第一节矩阵位移法概述结构力学传统方法与结构矩阵分析方法，二者同源而有别：

在原理上同源，在作法上有别前者在“手算”的年代形成，后者则着眼于“电算”，计算手段的不同，引起计算方法的差异。与传统的力法、位移法相对应，在结构矩阵分析中也有矩阵力法和矩阵位移法，或称柔度法与刚度法。矩阵位移法由于具有易于实现计算过程程序化的优点而广为流传。矩阵位移法是有限元法的雏形，因此结构矩阵分析有时也称为杆件结构的有限元法。在本章中将使用有限元法中的一些术语和提法。7现在是7页\一共有48页\编辑于星期四1、矩阵位移法的基本思路第一节矩阵位移法概述力法

需要选择基本体系和多余约束。所以较多地依赖于结构的具体情况，不宜实现计算机计算的自动化，但其优点是计算出的结果就是力。位移法

是先求结点位移，再换算成力，该法的计算自动化和通用性强，目前广为采用。结构结点力杆件杆端力杆件端点位移结构结点位移位移法力法位移法与力法之由于选取的基本未知量不同，因此计算次序不同a、方法的选择8现在是8页\一共有48页\编辑于星期四1、矩阵位移法的基本思路第一节矩阵位移法概述b、基本假设和基本原理线弹性、小变形。满足叠加原理、功能原理c、正负号规定（采用右手法则）杆端内力规定当与坐标轴正方向一致时为正；杆端位移和结点位移规定当与坐标轴正方向一致时为正。结点外力规定当与坐标轴正方向一致时为正；9现在是9页\一共有48页\编辑于星期四1、矩阵位移法的基本思路第一节矩阵位移法概述化整为零

（离散化、单元分析）集零为整

（结点力平衡、位移协调）

先把结构拆开，分解成若干个单元（在杆件结构中，一般把每个杆件取作一个单元），这个过程称作离散化。然后按单元力学性质对每个单元分析建立单元刚度方程，在满足变形条件和平衡条件的前提下，将这些单元集合成整体求解。在一分一合，先拆后搭的过程中，把复杂结构的计算问题转化为简单单元分析和集合问题。矩阵位移法的要点：10现在是10页\一共有48页\编辑于星期四2、单元划分第一节矩阵位移法概述将一个在荷载作用下的连续结构剖分成若干个各自独立的单元，原结构可以看成是由各单元在连接点（称结点）连接而成的体系——化整为零在杆件结构矩阵分析中，一般是把杆件的转折点、汇交点、边界点、突变点或集中荷载作用点等列为结点，结点之间的杆件部分作为单元。11现在是11页\一共有48页\编辑于星期四2、单元划分第一节矩阵位移法概述将一个在荷载作用下的连续结构剖分成若干个各自独立的单元，原结构可以看成是由各单元在连接点（称结点）连接而成的体系——化整为零为了减少基本未知量的数目，跨间集中荷载作用点可不作为结点，但要计算跨间荷载的等效结点荷载；跨间结点也可不作为结点，但要推导相应的单元刚度矩阵，编程序麻烦。12现在是12页\一共有48页\编辑于星期四第二节单元分析（局部坐标系下的单元分析）1、坐标系的选择:

在矩阵位移法中采用两种坐标系：

局部坐标系和整体坐标系。单元分析的目的是以结点位移为基本未知量，分析每个单元的结点力和结点位移及荷载之间的关系，即建立单元刚度方程，并用矩阵形式表示。整体坐标局部坐标FPxy13现在是13页\一共有48页\编辑于星期四第二节单元分析（局部坐标系下的单元分析）2、局部坐标系中的单元刚度矩阵采用局部坐标系（以杆的轴线作为x轴）时，杆端力及杆端位移的正方向以坐标轴正方向为正。杆端位移：杆端内力：杆件方向：

l14现在是14页\一共有48页\编辑于星期四第二节单元分析（局部坐标系下的单元分析）2、局部坐标系中的单元刚度矩阵局部坐标系下的单刚方程15现在是15页\一共有48页\编辑于星期四第二节单元分析（局部坐标系下的单元分析）2、局部坐标系中的单元刚度矩阵刚度系数的物理意义：单元刚度矩阵是杆端力与杆端位移之物理关系；矩阵的阶数与杆端位移分量数相等；表示引起的杆端力的大小。16现在是16页\一共有48页\编辑于星期四第二节单元分析（局部坐标系下的单元分析）3、局部坐标系中的单元刚度矩阵性质单刚一般具有奇异性：单刚具有对称性：由反力互等定理可知受力角度：存在刚体位移，杆端位移无法唯一确定数学角度：向量相关，矩阵不可逆，行列式为零。17现在是17页\一共有48页\编辑于星期四第二节单元分析（局部坐标系下的单元分析）3、局部坐标系中的单元刚度矩阵性质与单元刚度方程相应的正、反两类问题

为不平衡力系时，无解；为平衡力系时，有无穷多组解。

为任何值时，都有对应的唯一解，且总是平衡力系。将单元视为两端自由的杆件，直接加在自由端作为指定的杆端力。将单元视为两端有人为约束控制的杆件。控制附加约束加以指定。解的性质力学模型反问题正问题18现在是18页\一共有48页\编辑于星期四第三节单元分析（整体坐标系下的单元分析）1、单元坐标转换矩阵

同类单元在局部坐标系中具有相同的简洁形式。但不同单元在复杂结构中方位不尽相同，在整体分析中，为使分量叠加方便，需选择一个统一公共坐标系——整体坐标系。按整体坐标系来建立各单元的刚度矩阵。局部坐标系整体坐标系变换19现在是19页\一共有48页\编辑于星期四第三节单元分析（整体坐标系下的单元分析）1、单元坐标转换矩阵整体坐标系下的分量局部坐标系下的分量xy两种坐标系中单元的杆端位移转换关系为：用整体量表示局部量？20现在是20页\一共有48页\编辑于星期四第三节单元分析（整体坐标系下的单元分析）整体坐标系下的分量局部坐标系下的分量xy两种坐标系中单元的杆端力转换关系为：1、单元坐标转换矩阵用局部量表示整体量？21现在是21页\一共有48页\编辑于星期四第三节单元分析（整体坐标系下的单元分析）整体坐标系下的单刚与局部坐标系下的单刚性质相同2、整体坐标系中的单元刚度矩阵

22现在是22页\一共有48页\编辑于星期四第三节单元分析（整体坐标系下的单元分析）3、整体坐标系中的单元刚度矩阵的特性

整体坐标系中的单元刚度矩阵与局部坐标系中的单元刚度矩阵有类似的特性（对称、奇异）。另外，局部坐标系中的单元刚度矩阵只与单元的几何形状、物理常数有关；整体坐标系中的单元刚度矩阵不仅与单元的几何形状、物理常数有关，还与单元的位置和方位有关。

23现在是23页\一共有48页\编辑于星期四第四节整体分析整体刚度方程是整体结构的结点力与结点位移之间的关系式，是通过考虑结构的变形连续条件和平衡条件建立起来的。无论何种结构，其整体刚度方程都具有统一的形式：[K]是整体刚度矩阵；{Δ}结构的结点位移列向量；{P}结构的结点力列向量。利用结点位移协调和结点力平衡条件将各单元整合到一起,得到一个关于结点位移的线性代数方程——集零为整1、整体刚度矩阵的集成24现在是24页\一共有48页\编辑于星期四第四节整体分析1、整体刚度矩阵的集成

由变形连续条件知，结点发生单位位移，交于该结点的各单元的杆端也发生单位位移，由此产生的单元杆端力应是单元刚度矩阵中的元素。结点单位位移产生的结点力是整体刚度矩阵中的元素。由平衡条件知交于某结点的各单元杆端力之和等于该结点的相应结点力。故整体刚度矩阵中的元素应是由对应的单元刚度矩阵中的元素叠加而成。综上所述，整体刚度矩阵可以根据单元的结点位移分量的局部码和总码之间的对应关系，由单元刚度矩阵集成结构整体刚度矩阵。25现在是25页\一共有48页\编辑于星期四第四节整体分析FPxyFR1FR2FR3FR41(1,2)2(3,4)213(5,6)FPFR1FR2FR3FR41、整体刚度矩阵的集成离散26现在是26页\一共有48页\编辑于星期四第四节整体分析分别写出整体坐标系下单元刚度方程1、整体刚度矩阵的集成27现在是27页\一共有48页\编辑于星期四第四节整体分析将离散单元集合时应满足位移协调和平衡条件位移协调121(1,2)2(3,4)3(5,6)FP1、整体刚度矩阵的集成28现在是28页\一共有48页\编辑于星期四第四节整体分析将离散单元集合时应满足位移协调和平衡条件结点平衡121(1,2)2(3,4)3(5,6)FPFR1FR2FR3FR41、整体刚度矩阵的集成29现在是29页\一共有48页\编辑于星期四第四节整体分析1、整体刚度矩阵的集成集成整体刚度矩阵的关键，是确定单元刚度矩阵中的元素在整体刚度矩阵中的位置。首先要知道单元的结点位移分量的局部码和总码之间的对应关系，即单元定位向量。它是单元结点位移总码按局部码顺序排列而成的向量记为{λ}。单元结点分量按单元定位向量向整体刚度矩阵安装12(5,6)(1,2)(3,4)123FPxy30现在是30页\一共有48页\编辑于星期四第四节整体分析1、整体刚度矩阵的集成将单元刚度矩阵按单元定位向量扩展为单元贡献矩阵（换码扩阵）31现在是31页\一共有48页\编辑于星期四第四节整体分析1、整体刚度矩阵的集成按位叠加得整体刚度矩阵32现在是32页\一共有48页\编辑于星期四第四节整体分析1、整体刚度矩阵的集成按位叠加得整体刚度矩阵33现在是33页\一共有48页\编辑于星期四第四节整体分析根据集成法原理，计算过程分步进行：1、整体刚度矩阵的集成34现在是34页\一共有48页\编辑于星期四第四节整体分析1、整体刚度矩阵的集成整体刚度矩阵集成的两步可以合成一步完成：边换码，边叠加，对号入座(1)将[k]e

中元素重新在[K]e中定位

结点量值的两种编码：总码、局码；

量值分量在两种编码中的对应关系，由单元定位向量来表示；(2)将所有

[k]e

叠加生成

[K]

单刚元素在单元刚阵中按局码顺序排列；而在单元贡献阵中按总码顺序排列，即需要“换码重排座”35现在是35页\一共有48页\编辑于星期四第四节整体分析2、整体刚度矩阵系数的物理意义1.对称性3.稀疏性整体刚度矩阵是结点力与结点位移之间的物理关系；矩阵的阶数与结点位移分量数相等；Kij表示第j个结点位移分量为1时，第i个结点力分量的大小。3、整体刚度矩阵的性质2.奇异性36现在是36页\一共有48页\编辑于星期四第五节内力计算由前面的分析得到整体刚度矩阵和结点荷载列阵，形成关于结点位移的线性代数方程组；求解线性方程组可以得到结点位移列阵；利用结点位移列阵计算杆端内力。计算结点位移由定位向量e

提取单元杆端位移37现在是37页\一共有48页\编辑于星期四第五节内力计算FPxyl1l2EA=750l1

=4l2

=3FP

=40解:1.离散化局码12总码位移码单元编码(5,6)(1,2)(3,4)123FPxy例题：用矩阵位移法求结构内力。38现在是38页\一共有48页\编辑于星期四第五节内力计算FPxyl1l2例题：用矩阵位移法求结构内力。2.计算局部单刚3.求变换矩阵解:1.离散化39现在是39页\一共有48页\编辑于星期四第五节内力计算FPxyl1l2例题：用矩阵位移法求结构内力。2.计算局部单刚3.求变换矩阵解:1.离散化125612561234561234564.求单刚扩展矩阵40现在是40页\一共有48页\编辑于星期四第五节内力计算FPxyl1l2例题：用矩阵位移法求结构内力。2.计算局部单刚3.求变换矩阵解:1.离散化345634561234561234564.求单刚扩展矩阵41现在是41页\一共有48页\编辑于星期四第五节内力