第2章-有限元法的直接刚度法-1梁单元课件

知识点：直梁和平面刚架的直接刚度法重点：梁单元杆和刚架单元的自由度单元的坐标变换难点：直接刚度法的计算过程与物理意义第1章有限元法的直接刚度法知识点：第1章有限元法的直接刚度法1Ⅰ.关于梁和弯曲的概念受力特点：杆件在包含其轴线的纵向平面内，承受垂直于轴线的横向外力或外力偶作用。变形特点：直杆的轴线在变形后变为曲线。梁——以弯曲为主要变形的杆件称为梁。第二章有限元法的直接刚度法Ⅰ.关于梁和弯曲的概念受力特点：梁——以弯曲为2弯曲变形第二章有限元法的直接刚度法弯曲变形第二章有限元法的直接刚度法3工程实例F2F1第二章有限元法的直接刚度法工程实例F2F1第二章有限元法的直接刚度法4纵向对称面

对称弯曲——外力作用于梁的纵向对称面内，因而变形后梁的轴线(挠曲线)是在该纵对称面内的平面曲线。

非对称弯曲——梁不具有纵对称面(例如Z形截面梁)，因而挠曲线无与它对称的纵向平面；或梁虽有纵对称面但外力并不作用在纵对称面内，从而挠曲线不与梁的纵对称面一致。第二章有限元法的直接刚度法纵向对称面对称弯曲——外力作用于梁的纵向对称面内，因5本章讨论对称弯曲时梁的内力和应力。对称弯曲时和特定条件下的非对称弯曲时，梁的挠曲线与外力所在平面相重合，这种弯曲称为平面弯曲。第二章有限元法的直接刚度法本章讨论对称弯曲时梁的内力和应力。对称弯曲时和特定条6悬臂梁(2)梁的基本形式简支梁外伸梁第二章有限元法的直接刚度法悬臂梁(2)梁的基本形式简支梁外伸梁第二章有限元法的72.1直梁的有限元分析(a)直梁模型图2.1直梁(b)直梁的有限元模型第二章有限元法的直接刚度法以直梁为例来说明有限元法的直接刚度法。如图2.1(a)所示直梁，已知E、I、Z、M，AB=BC=CD=l，IAC=2l，ICD=l。2.1直梁的有限元分析图2.1直梁(b)直梁的有限元82.1.1划分单元两个节点之间的杆件构成一个单元，杆件结构的节点可按以下原则选取：1、杆件的交点一定要选为节点。2、阶梯形杆截面变化处一定要取为节点。3、支承点和自由端要取为节点。4、集中载荷作用处要取为节点。5、欲求位移的点要取为节点。6、单元长度不要相差太多。第二章有限元法的直接刚度法2.1.1划分单元两个节点之间的杆件构成一个单元，杆件结构的9按照杆件结构划分单元的原则，对图2.1(a)所示结构划分的单元如图2.1(b)所示图2.1(a)单元的节点位移(b)单元的节点力第二章有限元法的直接刚度法按照杆件结构划分单元的原则，对图2.1(a)所示结构划分的单102.1直梁的有限元分析任取一单元进行分析。根据材料力学的知识，梁单元上每个节点的节点位移分量有2个：挠度和转角，一般规定，向上为正，逆时针为正。写成列阵形式见式（2-1），表示节点的节点位移。（2-1）

图2.2(a)所示梁单元有、两个节点，共有4个节点位移分量：、、、，可用一个列阵表示，式（2-2）称为单元的节点位移列阵。（2-2）

第二章有限元法的直接刚度法2.1直梁的有限元分析任取一单元进行112.1直梁的有限元分析根据材料力学的知识，梁在外力作用下，横截面上的内力有2个：剪力、弯矩。所以，梁单元上每个节点的节点力有2个，用、来表示，规定：向上为正，逆时针为正。写成列阵形式见式（2-3），表示节点的节点力。

（2-3）图2.2(b)所示梁单元共有4个节点力分量：、、、，可用一个列阵表示，式（2-4）称为单元的节点力列阵。（2-4）

第二章有限元法的直接刚度法2.1直梁的有限元分析根据材料力学的知122.1直梁的有限元分析

梁单元上每个节点的节点载荷有2个：横向力和力偶，一般规定，向上为正，逆时针为正。写成列阵形式见式（2-5），表示节点的节点载荷。

（2-5）同理：

（2-6）第二章有限元法的直接刚度法2.1直梁的有限元分析梁单元上每个节点的节132.1直梁的有限元分析

节点力和节点载荷的区别：节点力是单元和节点之间的作用力，如果取整个结构为研究对象，节点力是内力；而节点载荷是结构在节点上所受到的外载荷或等效移置到节点上的外载荷。（2-7）第二章有限元法的直接刚度法2.1直梁的有限元分析

节点力和节点载荷142.1直梁的有限元分析根据材料力学的知识可知，在弹性范围和小变形的前提下，节点力和节点位移之间是线性关系。所以，单元的节点力和节点位移的关系可以表示为：（2-9）写成矩阵形式：（2-10）第二章有限元法的直接刚度法2.1直梁的有限元分析根据材料力学的知识可知152.1直梁的有限元分析简写为: （2-11）其中为单元节点力列阵，为单元节点位移列阵，称为单元刚度矩阵。单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵。单元刚度矩阵中各元素的物理意义：

(a)单元的节点位移图2.3单元刚度矩阵第1列元素的意义第二章有限元法的直接刚度法2.1直梁的有限元分析简写为: 第二章有162.1直梁的有限元分析

在点固定，令点有如图2.3(a)所示的位移，即有，，，。代入公式（2-10）中，得

（2-12）由式（2-12）可知，单元刚度矩阵中第一列元素的物理意义：为了使梁单元产生如图2.3(a)所示的位移，作用在单元节点上的节点力。第二章有限元法的直接刚度法2.1直梁的有限元分析在点固定，令172.1直梁的有限元分析的物理意义：单元第１个节点位移分量等于1，其它节点位移分量等于0时，对应的第1个节点力分量。的物理意义：单元第１个节点位移分量等于1，其它节点位移分量等于0时，对应的第2个节点力分量。的物理意义：单元第１个节点位移分量等于1，其它节点位移分量等于0时，对应的第3个节点力分量。的物理意义：单元第１个节点位移分量等于1，其它节点位移分量等于0时，对应的第4个节点力分量。

单元刚度矩阵中元素的物理意义：单元第个节点位移分量等于1，其它节点位移分量等于0时，对应的第个节点力分量。第二章有限元法的直接刚度法2.1直梁的有限元分析的物理意义：单元第１182.1直梁的有限元分析求单元刚度矩阵的第一列元素，由叠加原理，可得：（2-13）

其中，、为图2.3（b）所示单独作用所产生的位移，、为图2.3（b）所示单独作用所产生的位移。图2.3(b)节点i的节点力第二章有限元法的直接刚度法2.1直梁的有限元分析求单元刚度矩阵的第192.1直梁的有限元分析

可得到，，，（2-14）（2-15）

解方程（2-15）得：（2-16）第二章有限元法的直接刚度法2.1直梁的有限元分析第二章有限元法的直接刚度法20对梁单元分析受力，如图2.3（c）所示，列平衡方程（2-17）

图2.3(c)单元的节点力解方程（2-17）得

（2-18）

2.1直梁的有限元分析第二章有限元法的直接刚度法对梁单元分析受力，如图2.3（c）所示，列平衡方程212.1直梁的有限元分析单元刚度矩阵中第二列元素的物理意义是：，，，时，作用在单元节点上的节点力，如图2.4所示。求单元刚度矩阵的第二列元素，由叠加原理，可得：（2-19）解方程（2-19）得：（2-20）

图2.4单元刚度矩阵第2列元素的意义

第二章有限元法的直接刚度法2.1直梁的有限元分析单元刚度矩阵222.1直梁的有限元分析对梁单元分析受力，列平衡方程，解得：（2-21）同理，可求出单元刚度矩阵中的第三、四列元素，从而得到单元刚度矩阵。

（2-22）

第二章有限元法的直接刚度法2.1直梁的有限元分析对梁单元分析受力，列平衡方程232.1直梁的有限元分析从式（2-22）可以看出，单元刚度矩阵是一个对称矩阵，即。将单元刚度矩阵的公式，即式（2-22），应用于三个实际的梁单元，如图2.5所示，得到每个单元的节点力和节点位移的关系分别见式（2-23）、（2-24）和（2-25）。图2.5三个单元的受力图第二章有限元法的直接刚度法2.1直梁的有限元分析从式（2-22）可242.1直梁的有限元分析

（2-23）

（2-24）

（2-25）

第二章有限元法的直接刚度法2.1直梁的有限元分析第二章有限元法的直接刚度法252.1直梁的有限元分析单元的节点力和节点位移的关系，通常采用分块的方法表示，如2号单元的节点力和节点位移的关系见式（2-24），可表示为如下形式：（2-26）

第二章有限元法的直接刚度法2.1直梁的有限元分析单元的节点力262.1直梁的有限元分析2.1.3建立节点平衡方程式取图2.1(b)中各节点为研究对象，如图2..6所示，列平衡方程式：（2-27）图2.6各节点的受力图第二章有限元法的直接刚度法2.1直梁的有限元分析2.1.3建立节点平衡方程式第二章272.1直梁的有限元分析选取节点1为研究对象，分析受力，列平衡方程式（2-27），解得

（2-28）选取节点2为研究对象，分析受力，列平衡方程式（2-27），解得

（2-29）第二章有限元法的直接刚度法2.1直梁的有限元分析选取节点1为研究对象，分析282.1直梁的有限元分析同理，分别选取节点3、节点4为研究对象，分析受力，列平衡方程式（2-27），解得、和、，整个结构共得到8个平衡方程，称为有限元的基本方程，写成矩阵形式为：

（2-30）简记为 （2-31）第二章有限元法的直接刚度法2.1直梁的有限元分析同理，分别选取节点3、节点292.1直梁的有限元分析其中：——整个结构的节点载荷列阵（包括外载荷、约束反力）；——整个结构的节点位移列阵；——结构的整体刚度矩阵，又称总刚矩阵。整体刚度矩阵具有下列性质和特点：对称性：。奇异性：，不存在，这是因为尚未加入边界约束条件之前，整个系统可以作刚体运动，因而位移不是唯一的。只有加入边界约束条件后，约束了结构的刚体位移，才能使成为正定矩阵，从而得到位移的唯一解。稀疏性：整体刚度矩阵中非零的元素往往分布在对角线主元素的邻近，呈狭长的带状分布，这是因为任一节点只与围绕它的相连的单元发生联系，而其他单元的节点位移不会引起该节点处的节点力，所以整体刚度矩阵的每一行中会有大量的零元素。结构的节点数越多，整体刚度矩阵的这个特点就越明显。主对角线上的元素恒为正。第二章有限元法的直接刚度法2.1直梁的有限元分析其中：——整个结302.1直梁的有限元分析叠加法形成整体刚度矩阵的具体步骤如下：1）将单元刚度矩阵写成分块形式

（2-32）其中：——单元号；、——单元的两节点的编号；——单元在节点的节点力向量；——节点的节点位移向量；——单元上，节点单位位移在节点引起的节点力向量。1号单元： （2-33）2号单元： （2-34）3号单元： （2-35）第二章有限元法的直接刚度法2.1直梁的有限元分析叠加法形成整体刚度矩阵312.1直梁的有限元分析2）将整体刚度矩阵写成分块形式（分块矩阵的阶数等于结构的节点数）

（2-36）3）叠加形成整体刚度矩阵 （2-37）第二章有限元法的直接刚度法2.1直梁的有限元分析2）将整体刚度矩阵写成分块形式（分块矩322.1直梁的有限元分析2.1.4求解有限元基本方程分析图2.1所示梁结构的边界约束条件，可以得出结构的位移边界条件为：，载荷边界条件为：，，。将边界约束条件代入有限元基本方程（2-30），得

（2-39）

对式（2-39）求解，可以求得节点的未知位移值：将求得的节点未知位移值回代入式（2-38），可求得节点的未知载荷值：第二章有限元法的直接刚度法2.1直梁的有限元分析2.1.4求解有限元基本方程第二章332.1直梁的有限元分析2.1.5举例例：如图2.7（a）所示梁，已知、、、，求截面的转角及约束力。（a）直梁模型图2.7直梁第二章有限元法的直接刚度法2.1直梁的有限元分析2.1.5举例第二章有限元法的直342.1直梁的有限元分析解：（1）划分单元，对单元和节点进行编号，如图2.7(b)所示(b)直梁的有限元分析模型图2.7第二章有限元法的直接刚度法2.1直梁的有限元分析解：（1）划分单元，对单元和节点进行编352.1直梁的有限元分析（2）分析每一个单元，得到单元刚度矩阵，并写成分块形式。直梁单元刚度矩阵的公式前面已求得，见式（2-22），由此可以得出1号单元和2号单元的单元刚度矩阵和（2-40）

（2-41）第二章有限元法的直接刚度法2.1直梁的有限元分析（2）分析每一个单元，得到单元刚度矩阵362.1直梁的有限元分析（3）叠加法形成整体刚度矩阵。

（2-42）

第二章有限元法的直接刚度法2.1直梁的有限元分析（3）叠加法形成整体刚度矩阵。第二章372.1直梁的有限元分析（4）求解有限元基本方程。有限元基本方程为：（2-43）

第二章有限元法的直接刚度法2.1直梁的有限元分析（4）求解有限元基本方程。有限元基本方382.1直梁的有限元分析

分析结构的边界约束条件，位移边界条件为：，载荷边界条件为：。将结构的边界约束条件代入有限元基本方程（2-43），得（2-44）在上式中划掉已知位移所在的行和列，得

（2-45）解之，可得 （2-46）第二章有限元法的直接刚度法2.1直梁的有限元分析分析结构的边界约束条件392.1直梁的有限元分析将求得的回代入有限元基本方程（2-44），求得各节点的未知载荷，即约束力为：

（2-47）

第二章有限元法的直接刚度法2.1直梁的有限元分析将求得的402.1直梁的有限元分析2.1.6总结通过以上推导、分析和举例，我们可以总结出有限元求解问题的基本过程为：（1）结构的离散化：划分单元，对单元和节点进行编号。（2）单元分析，得到每一个单元的刚度矩阵，并将单元刚度矩阵写成分块形式。（3）求出结构的整体刚度矩阵，得到有限元基本方程。求出结构的整体刚度矩阵，有两种方法：①通过建立节点平衡方程式得到整体刚度矩阵。即对整个结构的每个节点进行受力分析，并列节点平衡方程式，得到整个结构的有限元基本方程，从而得到整体刚度矩阵。②通过叠加法得到整体刚度矩阵。即将单元刚度矩阵直接叠加形成整体刚度矩阵。

（4）分析结构的边界约束条件，求解有限元基本方程。将位移边界条件和载荷边界条件代入有限元基本方程，并降阶处理（划掉已知位移所在的行和列），求解降阶后的基本方程，求得节点的未知位移值，并将结果回代入没有降阶的有限元基本方程，求得节点的未知载荷值等。第二章有限元法的直接刚度法2.1直梁的有限元分析2.1.6总结第二章有限元法的直41知识点：直梁和平面刚架的直接刚度法重点：梁单元杆和刚架单元的自由度单元的坐标变换难点：直接刚度法的计算过程与物理意义第1章有限元法的直接刚度法知识点：第1章有限元法的直接刚度法42Ⅰ.关于梁和弯曲的概念受力特点：杆件在包含其轴线的纵向平面内，承受垂直于轴线的横向外力或外力偶作用。变形特点：直杆的轴线在变形后变为曲线。梁——以弯曲为主要变形的杆件称为梁。第二章有限元法的直接刚度法Ⅰ.关于梁和弯曲的概念受力特点：梁——以弯曲为43弯曲变形第二章有限元法的直接刚度法弯曲变形第二章有限元法的直接刚度法44工程实例F2F1第二章有限元法的直接刚度法工程实例F2F1第二章有限元法的直接刚度法45纵向对称面

对称弯曲——外力作用于梁的纵向对称面内，因而变形后梁的轴线(挠曲线)是在该纵对称面内的平面曲线。

非对称弯曲——梁不具有纵对称面(例如Z形截面梁)，因而挠曲线无与它对称的纵向平面；或梁虽有纵对称面但外力并不作用在纵对称面内，从而挠曲线不与梁的纵对称面一致。第二章有限元法的直接刚度法纵向对称面对称弯曲——外力作用于梁的纵向对称面内，因46本章讨论对称弯曲时梁的内力和应力。对称弯曲时和特定条件下的非对称弯曲时，梁的挠曲线与外力所在平面相重合，这种弯曲称为平面弯曲。第二章有限元法的直接刚度法本章讨论对称弯曲时梁的内力和应力。对称弯曲时和特定条47悬臂梁(2)梁的基本形式简支梁外伸梁第二章有限元法的直接刚度法悬臂梁(2)梁的基本形式简支梁外伸梁第二章有限元法的482.1直梁的有限元分析(a)直梁模型图2.1直梁(b)直梁的有限元模型第二章有限元法的直接刚度法以直梁为例来说明有限元法的直接刚度法。如图2.1(a)所示直梁，已知E、I、Z、M，AB=BC=CD=l，IAC=2l，ICD=l。2.1直梁的有限元分析图2.1直梁(b)直梁的有限元492.1.1划分单元两个节点之间的杆件构成一个单元，杆件结构的节点可按以下原则选取：1、杆件的交点一定要选为节点。2、阶梯形杆截面变化处一定要取为节点。3、支承点和自由端要取为节点。4、集中载荷作用处要取为节点。5、欲求位移的点要取为节点。6、单元长度不要相差太多。第二章有限元法的直接刚度法2.1.1划分单元两个节点之间的杆件构成一个单元，杆件结构的50按照杆件结构划分单元的原则，对图2.1(a)所示结构划分的单元如图2.1(b)所示图2.1(a)单元的节点位移(b)单元的节点力第二章有限元法的直接刚度法按照杆件结构划分单元的原则，对图2.1(a)所示结构划分的单512.1直梁的有限元分析任取一单元进行分析。根据材料力学的知识，梁单元上每个节点的节点位移分量有2个：挠度和转角，一般规定，向上为正，逆时针为正。写成列阵形式见式（2-1），表示节点的节点位移。（2-1）

图2.2(a)所示梁单元有、两个节点，共有4个节点位移分量：、、、，可用一个列阵表示，式（2-2）称为单元的节点位移列阵。（2-2）

第二章有限元法的直接刚度法2.1直梁的有限元分析任取一单元进行522.1直梁的有限元分析根据材料力学的知识，梁在外力作用下，横截面上的内力有2个：剪力、弯矩。所以，梁单元上每个节点的节点力有2个，用、来表示，规定：向上为正，逆时针为正。写成列阵形式见式（2-3），表示节点的节点力。

（2-3）图2.2(b)所示梁单元共有4个节点力分量：、、、，可用一个列阵表示，式（2-4）称为单元的节点力列阵。（2-4）

第二章有限元法的直接刚度法2.1直梁的有限元分析根据材料力学的知532.1直梁的有限元分析

梁单元上每个节点的节点载荷有2个：横向力和力偶，一般规定，向上为正，逆时针为正。写成列阵形式见式（2-5），表示节点的节点载荷。

（2-5）同理：

（2-6）第二章有限元法的直接刚度法2.1直梁的有限元分析梁单元上每个节点的节542.1直梁的有限元分析

节点力和节点载荷的区别：节点力是单元和节点之间的作用力，如果取整个结构为研究对象，节点力是内力；而节点载荷是结构在节点上所受到的外载荷或等效移置到节点上的外载荷。（2-7）第二章有限元法的直接刚度法2.1直梁的有限元分析

节点力和节点载荷552.1直梁的有限元分析根据材料力学的知识可知，在弹性范围和小变形的前提下，节点力和节点位移之间是线性关系。所以，单元的节点力和节点位移的关系可以表示为：（2-9）写成矩阵形式：（2-10）第二章有限元法的直接刚度法2.1直梁的有限元分析根据材料力学的知识可知562.1直梁的有限元分析简写为: （2-11）其中为单元节点力列阵，为单元节点位移列阵，称为单元刚度矩阵。单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵。单元刚度矩阵中各元素的物理意义：

(a)单元的节点位移图2.3单元刚度矩阵第1列元素的意义第二章有限元法的直接刚度法2.1直梁的有限元分析简写为: 第二章有572.1直梁的有限元分析

在点固定，令点有如图2.3(a)所示的位移，即有，，，。代入公式（2-10）中，得

（2-12）由式（2-12）可知，单元刚度矩阵中第一列元素的物理意义：为了使梁单元产生如图2.3(a)所示的位移，作用在单元节点上的节点力。第二章有限元法的直接刚度法2.1直梁的有限元分析在点固定，令582.1直梁的有限元分析的物理意义：单元第１个节点位移分量等于1，其它节点位移分量等于0时，对应的第1个节点力分量。的物理意义：单元第１个节点位移分量等于1，其它节点位移分量等于0时，对应的第2个节点力分量。的物理意义：单元第１个节点位移分量等于1，其它节点位移分量等于0时，对应的第3个节点力分量。的物理意义：单元第１个节点位移分量等于1，其它节点位移分量等于0时，对应的第4个节点力分量。

单元刚度矩阵中元素的物理意义：单元第个节点位移分量等于1，其它节点位移分量等于0时，对应的第个节点力分量。第二章有限元法的直接刚度法2.1直梁的有限元分析的物理意义：单元第１592.1直梁的有限元分析求单元刚度矩阵的第一列元素，由叠加原理，可得：（2-13）

其中，、为图2.3（b）所示单独作用所产生的位移，、为图2.3（b）所示单独作用所产生的位移。图2.3(b)节点i的节点力第二章有限元法的直接刚度法2.1直梁的有限元分析求单元刚度矩阵的第602.1直梁的有限元分析

可得到，，，（2-14）（2-15）

解方程（2-15）得：（2-16）第二章有限元法的直接刚度法2.1直梁的有限元分析第二章有限元法的直接刚度法61对梁单元分析受力，如图2.3（c）所示，列平衡方程（2-17）

图2.3(c)单元的节点力解方程（2-17）得

（2-18）

2.1直梁的有限元分析第二章有限元法的直接刚度法对梁单元分析受力，如图2.3（c）所示，列平衡方程622.1直梁的有限元分析单元刚度矩阵中第二列元素的物理意义是：，，，时，作用在单元节点上的节点力，如图2.4所示。求单元刚度矩阵的第二列元素，由叠加原理，可得：（2-19）解方程（2-19）得：（2-20）

图2.4单元刚度矩阵第2列元素的意义

第二章有限元法的直接刚度法2.1直梁的有限元分析单元刚度矩阵632.1直梁的有限元分析对梁单元分析受力，列平衡方程，解得：（2-21）同理，可求出单元刚度矩阵中的第三、四列元素，从而得到单元刚度矩阵。

（2-22）

第二章有限元法的直接刚度法2.1直梁的有限元分析对梁单元分析受力，列平衡方程642.1直梁的有限元分析从式（2-22）可以看出，单元刚度矩阵是一个对称矩阵，即。将单元刚度矩阵的公式，即式（2-22），应用于三个实际的梁单元，如图2.5所示，得到每个单元的节点力和节点位移的关系分别见式（2-23）、（2-24）和（2-25）。图2.5三个单元的受力图第二章有限元法的直接刚度法2.1直梁的有限元分析从式（2-22）可652.1直梁的有限元分析

（2-23）

（2-24）

（2-25）

第二章有限元法的直接刚度法2.1直梁的有限元分析第二章有限元法的直接刚度法662.1直梁的有限元分析单元的节点力和节点位移的关系，通常采用分块的方法表示，如2号单元的节点力和节点位移的关系见式（2-24），可表示为如下形式：（2-26）

第二章有限元法的直接刚度法2.1直梁的有限元分析单元的节点力672.1直梁的有限元分析2.1.3建立节点平衡方程式取图2.1(b)中各节点为研究对象，如图2..6所示，列平衡方程式：（2-27）图2.6各节点的受力图第二章有限元法的直接刚度法2.1直梁的有限元分析2.1.3建立节点平衡方程式第二章682.1直梁的有限元分析选取节点1为研究对象，分析受力，列平衡方程式（2-27），解得

（2-28）选取节点2为研究对象，分析受力，列平衡方程式（2-27），解得

（2-29）第二章有限元法的直接刚度法2.1直梁的有限元分析选取节点1为研究对象，分析692.1直梁的有限元分析同理，分别选取节点3、节点4为研究对象，分析受力，列平衡方程式（2-27），解得、和、，整个结构共得到8个平衡方程，称为有限元的基本方程，写成矩阵形式为：

（2-30）简记为 （2-31）第二章有限元法的直接刚度法2.1直梁的有限元分析同理，分别选取节点3、节点702.1直梁的有限元分析其中：——整个结构的节点载荷列阵（包括外载荷、约束反力）；——整个结构的节点位移列阵；——结构的整体刚度矩阵，又称总刚矩阵。整体刚度矩阵具有下列性质和特点：对称性：。奇异性：，不存在，这是因为尚未加入边界约束条件之前，整个系统可以作刚体运动，因而位移不是唯一的。只有加入边界约束条件后，约束了结构的刚体位移，才能使成为正定矩阵，从而得到位移的唯一解。稀疏性：整体刚度矩阵中非零的元素往往分布在对角线主元素的邻近，呈狭长的带状分布，这是因为任一节点只与围绕它的相连的单元发生联系，而其他单元的节点位移不会引起该节点处的节点力，所以整体刚度矩阵的每一行中会有大量的零元素。结构的节点数越多，整体刚度矩阵的这个特点就越明显。主对角线上的元素恒为正。第二章有限元法的直接刚度法2.1直梁的有限元分析其中：——整个结712.1直梁的有限元分析叠加法形成整体刚度矩阵的具体步骤如下：1）将单元刚度矩阵写成分块形式

（2-32）其中：——单元号；、——单元的两节点的编号；——单元在节点的节点力向量；——节点的节点位移向量；——单元上，节点单位位移在节点引起的节点力向量。1号单元： （2-33）2号单元： （2-34）3号单元： （2-35）第二章有限元法的直接刚度法2.1直梁的有限元分析叠加法形成整体刚度矩阵722.1直梁的有限元分析2）将整体刚度矩阵写成分块形式（分块矩阵的阶数等于结构的节点数）

（2-36）3）叠加形成整体刚度矩阵 （2-37）第二章有限元法的直接刚度法2.1直梁的有限元分析2）将整体刚度矩阵写成分块形式（分块矩732.1直梁的有限元分析2.1.4求解有限元基本方程分析图2.1所示梁结构的边界约束条件，可以得出结构的位移边界条件为：，载荷边界条件为：，，。将边界约束条件代入有限元基本方程（2-30），得