在做数值计算分析中要验证正确性,刚度、质量和阻尼矩阵是有限元分析最终最重要的3个矩阵。本文只针对刚度矩阵,不过我认为质量和阻尼矩阵应该是同理。
王新敏老师的《ansys工程结构数值分析》中有具体的提取刚度矩阵的内容,书中P356页6.3.4节。
我采用的是python语言读取file.full文件内容,参考下面的连接
ANSYS刚度矩阵的提取与解析_长安CAE_新浪博客blog.sina.com.cn连接中的程序是可以运行的,就是在前面读取和输出文件的部分需要加点内容,可能是python的版本不同的原因,我的版本是python3.8
下面给出一个模型的提取结果,刚度矩阵是去除约束剩余的
采用的solid185单元,一面上的四个节点约束后刚度矩阵大小为12*12。
红框中的数值是矩阵对角线上的元素在矩阵具体数值(绿框)中的位置。蓝框中的数值是矩阵行数。