HB矩阵转换程序

leap

ansys输出的刚度阵是Harwell-Boeing形式，但HB矩阵对人们并不是很熟悉，于是从网上下载了HB矩阵的matlab转换程序，与大家分享。
另外，做了一个简单的板的算例（算例见《MATLAB有限元分析与应用》第240页），发现ansys输出的刚度矩阵并不对，不知道大家有没有出现过这样的问题。
ansys输出的HB矩阵如下
Stiffness matrix from ANSYS FULL file dumped into Harwell-Boeing format         
            73             9            32            32             0
RSA                        8             8            32             0
(I14)           (I14)           (d25.15)            (d25.15)            
             1
             7
            13
            18
            23
            27
            30
            32
            33
             1
             3
             5
             6
             7
             8
             2
             4
             5
             6
             7
             8
             3
             5
             6
             7
             8
             4
             5
             6
             7
             8
             5
             6
             7
             8
             6
             7
             8
             7
             8
             8
    0.476923076923069D+07
    0.100000000000008D+07
   -0.137499999999996D+07
   -0.721153846153845D+05
   -0.150961538461542D+07
   -0.937500000000000D+06
    0.476923076923069D+07
   -0.274999999999992D+07
    0.721153846153845D+05
    0.500000000000039D+06
   -0.937500000000000D+06
   -0.150961538461542D+07
    0.476923076923069D+07
   -0.150961538461542D+07
    0.937500000000000D+06
   -0.137499999999996D+07
    0.721153846153845D+05
    0.476923076923069D+07
    0.937500000000000D+06
   -0.150961538461542D+07
   -0.721153846153845D+05
    0.500000000000039D+06
    0.238461538461535D+07
   -0.937500000000000D+06
    0.500000000000039D+06
   -0.721153846153845D+05
    0.238461538461535D+07
    0.721153846153845D+05
   -0.137499999999996D+07
    0.238461538461535D+07
    0.937500000000000D+06
    0.238461538461535D+07
   有做过类似问题的请赐教！