matlab提取刚度矩阵,ANSYS刚度矩阵的提取与解析

微信公众号：CAE技术分享

结构有限元最后的求解过程总是归结到求解一个大型矩阵方程Ax=b，对于一些情况还需要考虑质量矩阵M和阻尼矩阵C。有限元程序在组装完所有单元的刚度矩阵后，考虑模型所施加的约束和载荷，最终将刚度矩阵进行一些处理，例如乘大数法，变成Ax=b的形式，其中A是刚度矩阵，b是节点载荷，x为待求的节点位移，A和b全为已知量。

基本上各类有限元软件均能够提取模型的刚度矩阵，此次针对刚度矩阵的提取与解析做一个例子，采用的软件是ANSYS经典。

在ANSYS中建立一个简单的模型，划分网格后共12个节点，定义材料参数，施加约束和载荷后求解。有限元模型如下所示。

图1

待求解结束后，会在工作目录下生成一个后缀为full的文件，之后即可进行刚度矩阵的提取。

通过主菜单，如下所示。

图2

选择Matrix后，弹出如下所示的界面。