帮我看看提取的刚度与质量矩阵

zhuweihong2004

我提取了的是4个单元5节点的梁模型的刚度与质量矩阵，看不明白是怎么回事，高手指点下，怎么看：
ROW     1 MATRIX     1
    84000000.       0.0000000       0.0000000       0.0000000       0.0000000   
    0.0000000       0.0000000       0.0000000       0.0000000      -84000000.   
    0.0000000       0.0000000   
ROW     1 MATRIX     2
   0.26000000       0.0000000       0.0000000       0.0000000       0.0000000   
    0.0000000       0.0000000       0.0000000       0.0000000      0.13000000   
    0.0000000       0.0000000   
ROW     2 MATRIX     1
    0.0000000       33667.200      -8416.8000       0.0000000       0.0000000   
    0.0000000       0.0000000       0.0000000       0.0000000       0.0000000   
   -33667.200      -8416.8000   
ROW     2 MATRIX     2
    0.0000000      0.28974555     -0.20429874E-01   0.0000000       0.0000000   
    0.0000000       0.0000000       0.0000000       0.0000000       0.0000000   
   0.10025445      0.12070126E-01
ROW     3 MATRIX     1
    0.0000000      -8416.8000       2805.6000       0.0000000       0.0000000   
    0.0000000       0.0000000       0.0000000       0.0000000       0.0000000   
    8416.8000       1402.8000   
ROW     3 MATRIX     2
    0.0000000     -0.20429874E-01  0.18580113E-02   0.0000000       0.0000000   
    0.0000000       0.0000000       0.0000000       0.0000000       0.0000000   
  -0.12070126E-01 -0.13930742E-02
ROW     4 MATRIX     1
    0.0000000       0.0000000       0.0000000      0.16800000E+09   0.0000000   
    0.0000000      -84000000.       0.0000000       0.0000000       0.0000000   
    0.0000000       0.0000000   
ROW     4 MATRIX     2
    0.0000000       0.0000000       0.0000000      0.52000000       0.0000000   
    0.0000000      0.13000000       0.0000000       0.0000000       0.0000000   
    0.0000000       0.0000000   
ROW     5 MATRIX     1
    0.0000000       0.0000000       0.0000000       0.0000000       67334.400   
    0.0000000       0.0000000      -33667.200       8416.8000       0.0000000   
    0.0000000       0.0000000   
ROW     5 MATRIX     2
    0.0000000       0.0000000       0.0000000       0.0000000      0.57949110   
    0.0000000       0.0000000      0.10025445     -0.12070126E-01   0.0000000   
    0.0000000       0.0000000   
ROW     6 MATRIX     1
    0.0000000       0.0000000       0.0000000       0.0000000       0.0000000   
    5611.2000       0.0000000      -8416.8000       1402.8000       0.0000000   
    0.0000000       0.0000000   
ROW     6 MATRIX     2
    0.0000000       0.0000000       0.0000000       0.0000000       0.0000000   
   0.37160225E-02   0.0000000      0.12070126E-01 -0.13930742E-02   0.0000000   
    0.0000000       0.0000000   
ROW     7 MATRIX     1
    0.0000000       0.0000000       0.0000000      -84000000.       0.0000000   
    0.0000000      0.16800000E+09   0.0000000       0.0000000      -84000000.   
    0.0000000       0.0000000   
ROW     7 MATRIX     2
    0.0000000       0.0000000       0.0000000      0.13000000       0.0000000   
    0.0000000      0.52000000       0.0000000       0.0000000      0.13000000   
    0.0000000       0.0000000   
ROW     8 MATRIX     1
    0.0000000       0.0000000       0.0000000       0.0000000      -33667.200   
   -8416.8000       0.0000000       67334.400       0.0000000       0.0000000   
   -33667.200       8416.8000   
ROW     8 MATRIX     2
    0.0000000       0.0000000       0.0000000       0.0000000      0.10025445   
   0.12070126E-01   0.0000000      0.57949110       0.0000000       0.0000000   
   0.10025445     -0.12070126E-01
ROW     9 MATRIX     1
    0.0000000       0.0000000       0.0000000       0.0000000       8416.8000   
    1402.8000       0.0000000       0.0000000       5611.2000       0.0000000   
   -8416.8000       1402.8000   
ROW     9 MATRIX     2
    0.0000000       0.0000000       0.0000000       0.0000000     -0.12070126E-01
  -0.13930742E-02   0.0000000       0.0000000      0.37160225E-02   0.0000000   
   0.12070126E-01 -0.13930742E-02
ROW    10 MATRIX     1
   -84000000.       0.0000000       0.0000000       0.0000000       0.0000000   
    0.0000000      -84000000.       0.0000000       0.0000000      0.16800000E+09
    0.0000000       0.0000000   
ROW    10 MATRIX     2
   0.13000000       0.0000000       0.0000000       0.0000000       0.0000000   
    0.0000000      0.13000000       0.0000000       0.0000000      0.52000000   
    0.0000000       0.0000000   
ROW    11 MATRIX     1
    0.0000000      -33667.200       8416.8000       0.0000000       0.0000000   
    0.0000000       0.0000000      -33667.200      -8416.8000       0.0000000   
    67334.400       0.0000000   
ROW    11 MATRIX     2
    0.0000000      0.10025445     -0.12070126E-01   0.0000000       0.0000000   
    0.0000000       0.0000000      0.10025445      0.12070126E-01   0.0000000   
   0.57949110       0.0000000   
ROW    12 MATRIX     1
    0.0000000      -8416.8000       1402.8000       0.0000000       0.0000000   
    0.0000000       0.0000000       8416.8000       1402.8000       0.0000000   
    0.0000000       5611.2000   
ROW    12 MATRIX     2
    0.0000000      0.12070126E-01 -0.13930742E-02   0.0000000       0.0000000   
    0.0000000       0.0000000     -0.12070126E-01 -0.13930742E-02   0.0000000   
    0.0000000      0.37160225E-02

[ 本帖最后由 zhuweihong2004 于 2009-5-16 19:39 编辑 ]

zhuweihong2004

上面貌似是12X12的两个矩阵，我不知道怎么能分别组成刚度和质量矩阵，还有个问题，我的是一个悬臂梁4单元五节点的结构，怎么会是12X12呢，如果是刚度矩阵式12X12，那么我的节点只有5个，得到的位移会是10X1或者15X1，那么k*u=f左边就不能相乘，为什么

zhuweihong2004

问题已解决，丢失了一个节点的（约束点的信息）——这个方法提取的刚度矩阵信息不全，建议使用/debug命令