7.4单元刚度矩阵组装及整体分析.docxVIP

7.4单元刚度矩阵组装及整体分析 741单刚组装形成总刚 根据全结构的平衡方程 1 . - /可知， 总体刚度矩阵是由单元刚度矩阵集合而成的 .如果一 个结构的计算模型分成盘个单元，那么总体刚度矩阵可由各个单元的刚度矩阵 ［］组装而成，即 ［K］是由每个单元的刚度矩阵［］的每个系数按其脚标编号“对号入座”叠加而成的 .这种叠加要 就必须要把它们转换到求在同一总体坐标系下进行.如果各单元的刚度矩阵是在单元局部坐标下建立的， 统一的结构（总体）坐标系.将总体坐标轴分别用 ■V7 ■■表 示，对某单元有 就必须要把它们转换到 式中,分别是局部坐标系和总体坐标系下的单元结点位移向量［T 式中, 分别是局部坐标系和总体坐标系下的单元结点位移向量 ［T］为坐标转换阵, 仅与两个坐标系的夹角有关，这样就有 .以后如不特别强调，总体坐标系下的各种物理参数 .以后如不特别强调，总体坐标系下的各种物理参数 均不加顶上的横杠 下面就通过简单的例子来说明如何形成总体刚度矩阵 .设有一个简单的平面结构， 选取6个结点，划 分为4个单元.单元及结点编号如图3-27所示.每个结点有两个自由度.总体刚度矩阵的组装过程可分为 下面几步： 4 5 6 图 7-27 (1)按单元局部编号顺序形成单元刚度矩阵 .图7-27中所示的单元③，结点的局部编号顺序为 形成的单元刚度矩阵以子矩阵的形式给出是 (2)将单元结点的局部编号换成总体编号, (2)将单元结点的局部编号 换成总体编号, 相应的把单元刚度矩阵中的子矩阵的下标也换 i j m3 = 冷 Kik 心冷冷 k 成总体编号.对下图3-27所示单元③的刚度矩阵转换成总体编号后为 ◎J ◎ J⑶ 吃) .单元③的各子矩阵将转换后的单元刚度矩阵的各子矩阵，投放到总体刚度矩阵的对应位置上 投放后情况如下: .单元③的各子矩阵 （4）将所有的单元都执行上述的 1 , 2, 3步，便可得到总体刚度矩阵，如式（3-9 ）.其中八H右上 角的上标：表示第：单元所累加上的子矩阵. [K] =心？ [K] = 心？ rl+3 x^l+5 护玄 jW-2-3 Kg p-3 K欝 % ” 2+3+4 (3-9) （5）从式（3-9 ）可看岀，总体刚度矩阵中的子矩阵 AB是单元刚度矩阵的子矩阵转换成总体编号后 ■ I [ 具有相同的下标.，一的那些子矩阵的累加.总体刚度矩阵第.行的非零子矩阵是由与结点.相联系的那 些单元的子矩阵向这行投放所构成的 742结点平衡方程 我们首先用结构力学方法建立结点平衡方程.连续介质用有限元法离散以后，取岀其中任意一个结点 ■，从环绕L_点各单元移置而来的结点载荷为 式中:表示对环绕结点丄的所有单元求和，环绕结点 』的各单元施加于结点[的结点力为 .因此，结点「的平衡方程可表示为 (3-10 ) 以［K］代入平衡方程，得到以结点位移表示的结点 「的平衡方程，对于每个结点，都可列岀平衡方程, 于是得到整个结构的平衡方程组如下： 式中，［K］为整体刚度矩阵，I： i为全部结点位移组成的向量， 「为全部结点载荷组成的向量 当然，如果各点的载荷向量也是在单元局部坐标下建立的，在合成以前，也应把它们转换到统一的 结构（总体）坐标系下，即 ｛即=［珈对 式中，｛Pi｝是总体坐标系下的结点载荷向量， ［』］为坐标转换阵. 743位移边界条件 在有限元法对结构进行整体分析时，建立了整体刚度矩阵 ［K］，也得到了结构的刚度平衡方程，即 .结构刚度方程的求解相当于总刚［K］求逆的过程. .结构刚度方程的求解相当于总刚 是对称、半正定的奇异矩阵，它的行列式值为零，不能立即求逆 .从物理意义看，在进行整体分析时，结 构是处于自由状态，在结点载荷 -的作用下，结构可以产生任意的刚体位移 .所以，在已知结点载荷 E2［的条件下，仍不能通过平衡方程惟一地解出结点位移 工.为了使问题可解，必须对结构加以足够 的位移约束，也就是应用位移边界条件.首先要通过施加适当的约束，消除结构的钢体位移，再根据问题 要求设定其他已知位移.所以，处理位移边界条件在有限元分析步骤中十分重要 约束的种类包括使某些自由度上位移为零， 4二Q,或给定其位移值，4 = 3,还有给定支承刚 度等，本书涉及前两种.处理约束的方法，常用的有删行删列法、分块法、置大数法和置“ 1”法等，下 面分别予以介绍. 1、删行删列法 若结构的某些结点位移值为零时（即与刚性支座连接点的位移），则可将总体刚度矩阵中相应的行 列、删行删列划掉，然后将矩阵压缩即可求解 .这种方法的优点是道理简单.如果删去的行列很多，则总 体刚度矩阵的阶数可大大缩小.通常用人工计算时常采用该方法.若用计算机算题，在程序编制上必带来 麻烦，因为刚度矩阵压缩以后，刚度矩阵中各元素的下标必全改变 .因而