关于全局刚度矩阵计算的问题:如何得到类示例的全局刚度矩阵?类中的例子: 3332 31 12 13 11 23 21 22 3331 32 1311 12 2321 22 同理,可得到第二、三、四单元对应的单元刚度矩阵。第二个和第四个单元与第一个单元类似,得到的刚度矩阵在相应的节点处相等。每个单元的刚度矩阵展开为1212的方阵并相加。展开时刚度矩阵计算公式,除了对应的i个双行双列上的9个子矩阵,其余元素全部为0。所以整体刚度矩阵1112 13 2122 23 24 25 3132 33 35 36 4244 45 5354 55 56 6365 66 12以12为例,即第一元素的刚度矩阵 K22 与第二和第三元素的刚度矩阵之和。计算结果为:2222 22 22刚度矩阵计算公式,将该矩阵放入整体刚度矩阵的第二双行第二双列。然后通过(1,3,7,8,10,12行均为0)和上式的整体刚度矩阵,我们得到1112 13 2122 23 24 25 3132 33 35 36 4244 45 5354 55 56 6365 66 12 12 1011 12 1, 3, 7, 8, 10, 12 对计算没有影响,所以考虑到在位移矩阵中,第 1,3,7,8,10,12 行都是 0,第一个, 3、7、8、10、12的行和列对结果没有影响,所以可以简化为两个66方阵的乘积,即取2、4、5、6、9 , 11 扩展的全局刚度矩阵 由行和列组成的方阵,位移矩阵舍弃了几个受约束的自由度。以全局刚度矩阵的前几项为例: 1011 12 省略1、3、7、8、10、12行和列得到的全局刚度矩阵为:








