西安交通大学有限元分析版第二章

1、第二章 结构矩阵分析由于有限元方法起源于力学中的结构分析，本章的作用是通过三个典型问题说明有限元方法应用于结构分析时的一般步骤，并借此了解有限元方法的一些基本概念。§2-1平面桁架（直接法，结构矩阵分析中常用的力法，处理静定问题，位移法，可处理静定&静不定）本节讨论的对象是图2-1所示的平面桁架。组成桁架的各杆为等截面直杆，外载荷p直接作用于杆的铰接点（结点）上。为简单起见不妨设各杆的截面积均为，材料的弹性模量均为。我们可按下述步骤求得桁架的变形和内力。p图x u ppy v q图、 结构的离散化对结点及单元编号取组成桁架的每根杆为一个单元(该问题本身为一离散结构的力学问题)

2、，以, , 加以编号；取杆的铰接点为结点，以、加以编号（总体结点序号）。如图所示，即：我们所讨论的桁架包括三个单元、三个结点。各单元（杆）仅在结点处连接。、 建立总体坐标系并确定结点坐标和自由度为了描述结构的平衡需要建立一个坐标系，称为总体坐标系，以区别于以后出现的“局部坐标系”。总体坐标系的选择原则上不受限制，但希望使用方便。本节所选的总体坐标系示于图，坐标原点与结点重合。以u, v 分别表示沿 x, y方向的位移分量， p, q分别表示力沿 x, y 轴的力分量（投影）。在总体坐标系中各结点的坐标为：x u py v qijy v qx u pss图（x1, y1）=(0, 0 )、（x2

3、, y2）=(a, a )、（x3, y3）=(a, 0 )它们将作为程序的输入数据（几何参数）。每个结点有两个自由度，对结点、分别为u1， v1T 、 u2 ，v2T 、 u3， v3T若暂时不考虑支承约束条件，整个结构的结点自由度为u1 v1 u2 v2u3 v3T、 单元分析（建立结点力与结点位移之间的关系）取一个一般性的单元，设它的两个结点在结构中的编号为i, j （单元内部的结点序号）。由材料力学可知，杆的轴向刚度为EA/L。其中为杆的长度：（）单元局部坐标系现选取一典型单元对其进行单元分析，对所分析的单元按如下方式建立一个坐标系：原点：与结点i重合,x轴：沿 i ,j方向,y轴：与

4、x轴垂直。如图2-3所示。这个坐标系只属于一个单元，故称为单元局部坐标系，不同单元的单元局部坐标系一般是不相同的。在单元局部坐标系中可以规定：结点自由度 ui vi T， uj vj T ；单元结点自由度 u ui vi uj vj T。（）局部坐标系中的单元刚度矩阵在外载荷作用下，结构发生变形，单元必受到来自结点的作用力。桁架中的杆只承受轴向力，大小与杆的轴向伸长成正比在局部坐标系中这种特性可以得到清楚的表述（这一点也是引入局部坐标系的理由之一）。若以pi , qi , pj , qj 分别表示结点i, j作用于单元的力在 x, y 轴上的投影，由号单元的静力平衡有（图2-3）有(2-1-1

5、)用矩阵的形式可以写成若引入单元广义力矢量： 则上式可缩写为(2-1-2)其中(2-1-3)称为局部坐标系中的单元刚度矩阵,它只与杆的几个参数E、A、L有关，与杆的方位无关。（）坐标变换局部坐标系中的单元刚度矩阵公式简捷。但不同单元的局部坐标系一般不同，为了研究结构整体的平衡，必须将结点给单元的力以及相应的单元刚度矩阵转换到统一的坐标系总体坐标系。在总体坐标系中单元结点自由度 u ui vi uj vj T结点给单元的力r pi qi p j qj T在图2-3中，x 轴与x轴的夹角为 结点的位移分量的坐标变换为单元的位移分量的坐标变换为(2-1-)或缩写为类似，r 与 r 之间的转换关系为(

6、2-1-)由于(2-1-6)是正交矩阵，因此(2-1-7)也是正交矩阵。所以有将（2-1-4）、（2-1-5）代入（2-1-2）有从上式可得到(2-1-8)其中(2-1-9)称为单元在总体坐标系中的单元刚度矩阵。以后将会看到，（2-1-9）是一个具有普遍意义的公式。它表明，当单元的自由度由一种形式换成另一种形式时，单元刚度矩阵只需进行一次相似变换。对于平面桁架单元，将（2-1-3）、（2-1-6）、（2-1-7）代入（2-1-9）可得到更便于应用的单元刚度矩阵公式(2-1-10)（）具体结果由（2-1-10）可求得各单元的刚度矩阵的具体形式如下：单元：单元自由度u1 v1 u2 v2 T ，单

7、元刚度矩阵为 (2-1-11)单元：单元自由度u v u v T ，单元刚度矩阵为(2-1-12)单元：单元自由度u1 v1 u v T ，单元刚度矩阵为(2-1-13)请注意，单元刚度矩阵与单元自由度中位移分量的排列次序有关。如果改动这种排列次序，例如对号单元，将单元自由度次序由u1 v1 u2 v2 T改为 u2 v2 u1 v1 T，必然导致刚度矩阵（2-1-11）元素位置的变动。（）单元刚度矩阵的物理意义和特点设平面桁架单元在总体坐标系中刚度矩阵的一般形式为由（2-1-8），当单元结点位移为1 0 0 0 T时，在单元各结点上施加的力刚好为单元刚度矩阵中的第一列：k11 k21 k31

8、 k41 T。对k的其他各列也可做出类似的解释。即单元刚度矩阵的每一列相当于一组特定位移下的结点力，如表2-1所示。由图2-4可以获得更为直观的理解。表2-1平面桁架单元刚度矩阵的物理意义单元结点位移作用于单元的结点力1 0 0 0 Tk11 k21 k31 k41 T0 1 0 0 Tk12 k22 k32 k42 T0 0 1 0 Tk13 k23 k33 k43 T0 0 0 1 Tk14 k24 k34 k44 Tk11k41k21k311jik12k42k22k321jik44jk24ik34k141k33k43k23jik131图对图2-4中的各种情况，据平面力系的平衡条件应有这三

9、个关系说明，k的四个行向量中只有一个线性独立（四个元素有三个约束方程）。从以上分析可以看出，一般的单元刚度矩阵均具备以下两个特征。（对平面桁架单元而言，从（2-1-10）也可以得出这些结论）（i）单元刚度矩阵是对称矩阵，这是线性系统互易定理的具体体现。由于对称性，对行向量或列向量两者之一得到的结论，对另一个也适用。(ii)单元刚度矩阵是奇异矩阵。它的行向量（或列向量）线性相关，具有零特征值，detk=0。对平面桁架的单元刚度矩阵而言，它的四个行向量（或列向量）中只有一个线性独立，而k有三个零特征值。这三个零特征值对应的特征向量相当于三种独立的刚体位移模式：两个平移，一个旋转。这是我们在单元分析

10、中不考虑位移约束条件的自然结果。、总体刚度矩阵的组装总体平衡方程将图2-1所示的桁架中的支承约束以约束反力代替，如图2-5所示。下面来建立平衡问题的有限元方程。（） 结点平衡条件作用于图2-5每个结点上的外载荷、支座反力以及来自单元的力应处于平衡。R3YpR1X图R1Y以pi(m)、qi(m)示结点i作用于单元m的力在x, y轴上的投影，则单元m给结点 i的力在x, y轴上的投影应为pi(m)、qi(m)。对结点：对结点：对结点；可以合并成(2-1-1)式（2-1-14）的右边为外载荷和支反力。左边则为单元给结点的力，它们是未知的，但可以借助单元刚度矩阵以结点位移来表示。（） 单元刚度矩阵的扩

11、充为了表示（2-1-14）左边的各个列向量，设想将每个单元的自由度扩充到与结构总体自由度相同（本例为），并在单元刚度矩阵中补充零元素，由（2-1-11）、（2-1-12）、（2-1-13）和（2-1-8）可以用结点位移表示（2-1-14）左边的各列向量。由单元(2-1-15)由单元(2-1-16)由单元(2-1-17)（） 组装总体刚度矩阵将（2-1-15）、（2-1-16）、（2-1-17）代入（2-1-14）得到或(2-1-18)(2-1-19)其中上式称为没有考虑位移约束条件情况下的总体刚度矩阵（求和对所有单元进行）。对本节所分析的平面桁架有(2-1-20)要形成总体平衡方程（2-1-1

12、8），只需组装出它的右端项和总体刚度矩阵k就足够了，这只是同一件事的两种不同的提法而已。（2-1-19）用来“书写“组装总刚度矩阵的过程是简单而明了的。但事实上不能这样做，将所有km补充零元素扩充为km极大地浪费了宝贵的存贮空间，这些零元素仅起到使km的元素在总刚度矩阵中就位的作用。实际上采用的是另一种方法。如果已选定各结点位移在结构总体自由度中的排列次序为u1 v1 u2 v2 u3 v3 T，对每个单元在形成单元刚度矩阵的同时还形成了个定位数组LM，它将指出单元自由度的各分量在总体自由度中的序号，如下表：单元号单元自由度LM(1)LM(2)LM(3)LM(4)u1 v1 u2 v2T123

13、4u2 v2 u3 v3 T3456u1 v1 u3 v3T1256单元刚度矩阵中第s行第t列的元素kst加到总刚度矩阵的第LM(s)行LM(t)列即可。这一组装总体刚度矩阵的方法被形象的称为“对号入座”。从（2-1-20）可以看出总刚阵k是奇异阵，它的六个行向量（或列向量）中只有三个线性独立。这是尚未考虑位移约束条件，结构的刚体位移未受到限制的必然结果。（） 引入位移约束条件(2-1-21)由图2-1，平面桁架的位移约束条件为 u1=v1=v3=0 代入方程（2-1-18），得到它显然可以分成两个方程组(2-1-22)和(2-1-23)(2-1-24)（2-1-22）就是平常所说的平衡问题的

14、有限元方程，一般把它写成K为考虑了约束条件（2-1-21）后的总体刚度矩阵。U为非约束自由度位移向量，F为作用于这些自由度上的载荷向量。矩阵K相当从矩阵k中划去了第、行和第、列。这说明（2-1-21）这种零位移约束条件，可以通过对总刚度矩阵划行划列的方法来实现。至于其他形式的位移约束条件的实现方法，后面将有专门的一章加以讨论。实际上，在组装总刚度矩阵的过程中，只要不组装应划去的行和列即可直接得到K和F，即方程（2-1-24）。至于方程（2-1-23），尽管求得结点位移后可以由它求得支承反力，但是有限元分析中一般不形成这个方程。如确有必要求支座反力，将另找其他途径。、解有限元方程对于平衡问题，归

15、结为解一个代数方程。本例的自由度较少，从（2-1-22）可以方便地求得非约束自由度的位移被约束自由度根据（2-1-21）可直接赋零、 求单元内力桁架单元的内力只有一个轴力。由（2-1-1）S的大小和正负与 pj相同。由（2-1-2）和（2-1-4）取出它的第三行即得到(2-1-24)对单元、求得内力于下表单元号单元结点位移内力（以拉为正） 0 0 T 0 0T 0 0 0 0 T0由于本节所讨论的桁架是一个十分简单的静定桁架，用理论力学的知识即可得到各杆的内力，结果相同。但是，若桁架为静不定桁架且杆的数目有上千个，那么本节所讨论的方法原则上不会遇到任何困难，我们要解决的课题将转为如何管理有关的

16、大量数据和如何解一个数千阶的代数方程组这样一些技术问题。§2-2平面框架本节将讨论另一个典型问题平面框架，示于图。构成它的元件是梁。为讨论简单，不妨设外载荷只作用于梁的刚结点上。我们用同上一节类似的步骤求得框架的变形和内力。、 结构的离散化取每根梁为一个单元，编号为、。取梁之间的刚结点为结点，以、对结点进行编号。、 总体坐标系结点坐标和自由度2ax,u,paz,M4231（图2）y,v,qPx、y轴在框架平面内，原点与结点重合，z轴与框架平面垂直，如图所示。u, v分别表示各结点沿x, y 方向的位移，为绕z轴的转角。p, q;分别表示力在x, y轴上的投影，M为绕z轴之力偶矩。不难

17、定出框架四个结点的坐标为暂时不考虑支承约束，整个结构的结点总自由度为(2-2-1)考虑、结点处的约束条件结构的非约束结点自由度为(2-2-2)x,u,py,v,qz,Mij（图）、 单元分析取一个一般性的单元进行分析，设它的两个结点在结构中的编号为i、j。见图。（）建立单元局部坐标系原点：与结点i重合；x轴：沿i， j方向；y轴：在x y平面内，与x轴垂直；z轴：与z轴一致。单元局部坐标系中，各结点自由度为单元结点自由度可表示为（）局部坐标系中的单元刚度矩阵根据单元刚度矩阵的物理意义，它的元素相当单元结点位移在六种位移模式情况下的结点力。如果我们约定：（i）单元内弯曲刚度EI为常数； (ii)

18、 单元只在端点受到外力； (iii) 不计剪切变形。则u将是x的线性函数，v是x的三次函数， 。 利用材料力学知识不难求得梁的变形曲线和结点力，如表所示。设单元长度为L，截面积为。图中实线为现实构形，虚线为参照构形。表平面梁单元的位移场和结点力序号结点位移和单元位移场结点作用于单元的力ij1ijji1radij1i1ji1radj由此求得平面梁单元在局部坐标系中的单元刚度矩阵k(2-2-3)它是一个对称阵；也是一个奇阵。在单元分析中没有考虑位移约束条件，允许单元作刚体运动。我们还可以看到，局部坐标系为单元分析提供了方便。（）坐标变换为了讨论结构的平衡，需要把局部坐标系中形成的单元刚度矩阵转换到

19、总体坐标系。若图中x轴与x的夹角为 ，则两个坐标系中位移分量（以结点i为例）的转换关系为(2-2-)每个单元有两个结点，单元结点自由度的转换关系为(2-2-) 利用与2-1节中类似的步骤，可推出单元刚度矩阵的变换公式为(2-2-6)形式上与(2-1-9)完全相同。但请注意，矩阵t的定义不同，因而T的具体表达式也不同。利用（2-2-3）和（2-2-6）可求得单元、在总体坐标系中的单元刚度矩阵k1、k2、k3。、 总体刚度矩阵和载荷向量的组装本节侧重讨论在组装刚度矩阵的过程中实现划行划列，即由k直接组装K的方法。对每个单元，根据（2-2-1）和(2-2-2)可形成一个数组LM，如下表所示：单元号单

20、元自由度LM(1)LM(2)LM(3)LM(4)LM(5)LM(6)000123123456456000每个LM数组有六个元素，分别与单元自由度中的各位移分量对应。当这个位移分量根据（2-2-1）被约束时（例如u1, v1, 1, u4, v4, 4）对应的LM元素取零；若这个位移属于非约束自由度，对应的LM元素取该位移在（2-2-1）中的序号，利用LM数组即可实现由k直接组装K。对于k的第s行第t列的元素kst，若LM(s)和LM(t)中至少一个为零，则对该元素不予组装（相当于划行划列）；否则将kst迭加入总刚度矩阵K的第LM(s)行LM(t)第列。非约束自由度中只有u2的自由度上作用外载荷

21、，故总体载荷向量为(2-2-7)最后得到有限元方程、 有限元方程（2-2-7）得到非约束自由度结点位移，被约束自由度上的位移，据（2-2-1）直接赋零。、 由结点位移求单元内力。例如：轴力、剪力、弯矩，对空间梁单元还可求扭矩。§2-3 平面应力问题常应变三角形qqt图图2-8为一边长为a、厚度为t的正方形薄板。其中AB边固定，BC、CD边自由，AD边作用均布压力q。对这一问题，有限元分析的步骤是：、将ABCD划分（离散）为个三角形（单元），编号。各单元仅在顶点（结点）铰接，结点编号19。建立坐标系后，不难定出各结点的坐标（xi, yi）。、单元分析任取一个一般性的单元，如图所示。三个

22、结点的编号为i, j, k。结点位移为( ui, vi )、( uj, vj ) 、( uk, vk )单元结点位移为（） 假定单元内位移场u, v 是x, y的一次函数y, v, qijk图（）x, u, p(2-3-1)为待定常数，在结点处应有可解出：其中当 i, j, k的位置为逆时针排列时，2恒正，且等于三角形单元面积的两倍。将这些结果代入（2-3-1）有类似可得到可以合并成(2-3-2)（）单元的应变、应力利用（2-3-1）不难求得(2-3-3)其中(2-3-4)在假定单元内位移场u、v是 x, y的一次函数的前提下，单元内的应变和应力将是常数，故这种单元又称为常应变三角元。图（10）xyijk0(a)xyijkF1xF3xF2xF3yF1yF2yxxy0y(b)（）为了在单元内构成均匀应力场，必须在单元的各边施加均布载荷，它们的合力一定作用在各边的中点，如图2-10(a)所示。再将各边上的合力平分到这边的两点结点，由图2-10 (b)不难得出类似可求得 qi、pj、qj、pk、qk 并可合并写成根据单元刚度