matlab中矩阵怎么敲_Ansys刚度(质量、阻尼)矩阵的提取(part 1)

在做数值计算分析中要验证正确性，刚度、质量和阻尼矩阵是有限元分析最终最重要的3个矩阵。本文只针对刚度矩阵，不过我认为质量和阻尼矩阵应该是同理。

王新敏老师的《ansys工程结构数值分析》中有具体的提取刚度矩阵的内容，书中P356页6.3.4节。

我采用的是python语言读取file.full文件内容，参考下面的连接

ANSYS刚度矩阵的提取与解析_长安CAE_新浪博客​blog.sina.com.cn

连接中的程序是可以运行的，就是在前面读取和输出文件的部分需要加点内容，可能是python的版本不同的原因，我的版本是python3.8

下面给出一个模型的提取结果，刚度矩阵是去除约束剩余的

采用的solid185单元，一面上的四个节点约束后刚度矩阵大小为12*12。

红框中的数值是矩阵对角线上的元素在矩阵具体数值(绿框)中的位置。蓝框中的数值是矩阵行数。