要列出方程,即方程(1.6),其中∫是“等效节点力”,用常 规的有限元方法,它是可以从边界值∫求出来的。目前,我们的 目的是给出口上组合刚度矩阵的算法,因此可以暂且不管上述 边界条件的差别。 我们将在下章证明,存在一个n阶实矩阵X,使 Uk+1=Xy (1.7) X称为转移矩阵。以(1.)式代人(1.6)得 K-A+KX-AX)yo=0 但是y0可以是任意向量,所以X满足方程 A TX2-KX+A-0 (1.8) 令 K A R A 0 RI 其中是单位阵,方程(1、6)k可以写成如下形式: J R R (1,10) 设A,9分别是矩阵X的特征值和特征向量,则有 X9=Ag (1.1L) 在(1.10)中取无=1,令孙=9,并且注意到(].7).(.I1)得 g R1 nR (1.12) g 因此)与j 分别恳矩阵束R1-λR2的广义特征值与广义特 征向量1},为了求矩阵￥,我们先解上述广义特征值问题,求 出它的全部2n个特征值与特征向量,然后丛中选出ⅹ的n个特 征值与特征向量