zhuweihong2004 发表于 2009-5-16 17:55

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

我提取了的是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-010.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 发表于 2009-5-16 20:15

回复 楼主 zhuweihong2004 的帖子

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

zhuweihong2004 发表于 2009-5-17 22:42

回复 楼主 zhuweihong2004 的帖子

问题已解决，丢失了一个节点的（约束点的信息）——这个方法提取的刚度矩阵信息不全，建议使用/debug命令
页: [1]
查看完整版本: 帮我看看提取的刚度与质量矩阵