《有限元原理 结构矩阵分析(平面桁架 平面应力) 变分》由会员分享，可在线阅读，更多相关《有限元原理 结构矩阵分析(平面桁架 平面应力) 变分（23页珍藏版）》请在金锄头文库上搜索。
1、第二章第二章 结构矩阵分析结构矩阵分析由于有限元方法起源于力学中的结构分析，本章的作用是通过三个典型问题说明有限 元方法应用于结构分析时的一般步骤，并借此了解有限元方法的一些基本概念。2-1 平面桁架平面桁架 （直接法，结构矩阵分析中常用的力法，处理静定问题，位移法，可处理静定&静不定） 本节讨论的对象是图 2-1 所示的平面桁架。组成桁架的各杆为等截面直杆，外载荷 p 直接作用于杆的铰接点（结点）上。为简单起见不妨设各杆的截面积均为，材料的弹性 模量均为。我们可按下述步骤求得桁架的变形和内力。、结构的离散化 对结点及单元编号 取组成桁架的每根杆为一个单元(该问题本身为一离散结构的力学问题)，
2、以, , 加以编号；取杆的铰接点为结点，以、加以编号（总体结点序号） 。如图所 示，即：我们所讨论的桁架包括三个单元、三个结点。各单元（杆）仅在结点处连接。 、建立总体坐标系 并确定结点坐标和自由度 为了描述结构的平衡需要建立一个坐标系，称为总体坐标系，以区别于以后出现的 “局部坐标系” 。总体坐标系的选择原则上不受限制，但 希望使用方便。本节所选的总体坐标系示于图， 坐标原点与结点重合。以 u, v 分别表示沿 x, y 方向的 位移分量， p, q 分别表示力沿 x, y 轴的力分量（投影） 。在总体坐标系中各结点的坐标为：它们将作为程序的输入数据（几何参数） 。p 图 x p y 图x
3、u py v qijy v qx u pss图（x1, y1）=(0, 0 )、 （x2, y2）=(a, a )、 （x3, y3）=(a, 0 )每个结点有两个自由度，对结点、分别为若暂时不考虑支承约束条件，整个结构的结点自由度为、单元分析（建立结点力与结点位移之间的关系） 取一个一般性的单元，设它的两个结点在结构中的编号为 i, j （单元内部的结点序号） 。 由材料力学可知，杆的轴向刚度为 EA/L。其中为杆的长度：（）单元局部坐标系 现选取一典型单元对其进行单元分析，对所分析的单元按如下方式建立一个坐标系： 原点：与结点 i 重合, x轴：沿 i ,j 方向, y轴：与 x轴垂直。
4、如图 2-3 所示。这个坐标系只属于一个单元，故称为单元局部坐标系，不同单元的单元局 部坐标系一般是不相同的。在单元局部坐标系中可以规定： 结点自由度 ui vi T， uj vj T ； 单元结点自由度 u ui vi uj vj T 。 （）局部坐标系中的单元刚度矩阵 在外载荷作用下，结构发生变形，单元必受到来自结点的作用力。桁架中的杆只承受 轴向力，大小与杆的轴向伸长成正比在局部坐标系中这种特性可以得到清楚的表述（这一点也是引入局部坐标系的理由之一） 。 若以 pi , qi , pj , qj 分别表示结点 i, j 作用于单元的力在 x, y 轴上的投影，由 号单元的静力平衡有（图
5、2-3）有用矩阵的形式可以写成若引入单元广义力矢量： 则上式可缩写为 22 ijijyyxxLLLEAS)(ijjuuLEAp)(jiiuuLEAsp0ijqq(2-1-1) jjiijjiivuvuLEAqpqp0000010100000101 T jjiiqpqpru1， v1T 、 u2 ，v2T 、 u3， v3Tu1 v1 u2 v2 u3 v3 T其中称为局部坐标系中的单元刚度矩阵局部坐标系中的单元刚度矩阵,它只与杆的几个参数它只与杆的几个参数 E、A、L 有关，有关，与杆的方位无关与杆的方位无关。（）坐标变换 局部坐标系中的单元刚度矩阵公式简捷。但不同单元的局部坐标系一般不同，为
6、了研 究结构整体的平衡，必须将结点给单元的力以及相应的单元刚度矩阵转换到统一的坐标系 总体坐标系。在总体坐标系中 单元结点自由度 u ui vi uj vj T 结点给单元的力 r pi qi p j qj T 在图 2-3 中，x 轴与 x 轴的夹角为 结点的位移分量的坐标变换为单元的位移分量的坐标变换为或缩写为 类似，r 与 r 之间的转换关系为由于 ukr(2-1-2)(2-1-3) 0000010100000101LEAkLyyLxxijijsin,cos iiiiii vutvu vu cossinsincos jjjjjj vutvu vu cossinsincos jjiijji
7、ijjiivuvuvuvuvuvuTcossin00sincos0000cossin00sincos uuT(2-1-) rrT(2-1-) cossinsincost(2-1-6)是正交矩阵，因此也是正交矩阵。所以有将（2-1-4） 、 （2-1-5）代入（2-1-2）有从上式可得到其中称为单元在总体坐标系中的单元刚度矩阵。以后将会看到， （2-1-9）是一个具有普遍意义的公式。它表明，当单元的自由度由一 种形式换成另一种形式时，单元刚度矩阵只需进行一次相似变换。对于平面桁架单元，将 （2-1-3） 、 （2-1-6） 、 （2-1-7）代入（2-1-9）可得到更便于应用的单元刚度矩阵公式（
8、）具体结果 由（2-1-10）可求得各单元的刚度矩阵的具体形式如下： 单元：单元自由度 u1 v1 u2 v2 T ， 单元刚度矩阵为 单元：单元自由度 u v u v T ， 单元刚度矩阵为 tt 00T (2-1-7) T-1TT ukrTT ukukrTTT TTTkk(2-1-8)(2-1-9) 22222222sinsincossinsincossincoscossincoscossinsincossinsincossincoscossincoscosLEAk(2-1-10)o45 42 42 42 4242 42 42 4242 42 42 4242 42 42 421aEAko9
9、0 10100000101000002aEAk(2-1-12)(2-1-11)单元：单元自由度 u1 v1 u v T ， 单元刚度矩阵为请注意，单元刚度矩阵与单元自由度中位移分量的排列次序有关单元刚度矩阵与单元自由度中位移分量的排列次序有关。如果改动这种排列次序， 例如对号单元，将单元自由度次序由u1 v1 u2 v2 T改为 u2 v2 u1 v1 T，必然导致刚 度矩阵（2-1-11）元素位置的变动。 （）单元刚度矩阵的物理意义和特点 设平面桁架单元在总体坐标系中刚度矩阵的一般形式为由（2-1-8） ，当单元结点位移为1 0 0 0 T时，在单元各结点上施加的力刚好为单元刚度 矩阵中的第
10、一列：k11 k21 k31 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 To0 00000101000001013aEAk(2-1-13) 4443424
11、1343332312423222114131211kkkkkkkkkkkkkkkkk对图 2-4 中的各种情况，据平面力系的平衡条件应有这三个关系说明，k的四个行向量中只有一个线性独立（四个元素有三个约束方程） 。 从以上分析可以看出，一般的单元刚度矩阵均具备以下两个特征。 （对平面桁架单元而 言，从（2-1-10）也可以得出这些结论） （i）单元刚度矩阵是对称矩阵，）单元刚度矩阵是对称矩阵，这是线性系统互易定理的具体体现。由于对称性，对 行向量或列向量两者之一得到的结论，对另一个也适用。 (ii)单元刚度矩阵是奇异矩阵。单元刚度矩阵是奇异矩阵。它的行向量（或列向量）线性相关，具有零特征值，
12、detk=0。对平面桁架的单元刚度矩阵而言，它的四个行向量（或列向量）中只有一个线性 独立，而k有三个零特征值。这三个零特征值对应的特征向量相当于三种独立的刚体位移 模式：两个平移，一个旋转。这是我们在单元分析中不考虑位移约束条件的自然结果。 、总体刚度矩阵的组装 总体平衡方程 将图 2-1 所示的桁架中的支承约束以约束反力代替，如图 2-5 所示。下面来建立平 衡问题的有限元方程。 （） 结点平衡条件 作用于图 2-5 每个结点上的外载荷、支座反力以及来自单元的力应处于平衡。 以 pi(m)、qi(m) 示结点 i 作用于单元 m 的力在 x, y 轴上的投影，则单元 m 给结点 i 的力在
13、 x, y 轴上的投影应为 pi(m) 、qi(m)。 对结点：对结点：对结点；）（ 41000344231skyykxxkkkksijsijssss00)3( 1)1( 11)3( 1)1( 11 qqRppRYX00)2( 2)1( 2)2( 2)1( 2 qqppP00)3( 3)2( 33)3( 3)2( 3 qqRppY R3Y p R1X图R1Yk11k41k21k311jik12k42k22k321jik44jk24ik34k141k33k43k23jik131图可以合并成式（2-1-14）的右边为外载荷和支反力。左边则为单元给结点的力，它们是未知的，但可 以借助单元刚度矩阵以结点位移来表示。 （） 单元刚度矩阵的扩充 为了表示（2-1-14）左边的各个列向量，设想将每个单元的自由度扩充到与结构总体将每个单元的自
