ABAQUS输出单元刚度矩阵

Accompany

蔡和Lin的文章中，一道关键的推导步骤中，需要设定整体体系广义质量的对角和为单位值，并以此缩放刚度和质量矩阵，即可得到缩减模态等效体系。
此处虽然验证存在此关键一步，但并未理解为何这样做，或者说，是基于什么让他们想到这样处理就可以恰好推导成功呢？
难不成真的是存在“偶为之”的因素？

输出单元刚度矩阵的方法!
1.用命令：*ELEMENT MATRIX OUTPUT
只设定Required parameter：ELSET的话，由于结果文件（*.fil）是二进制文件，用文本编辑器打开是一堆我们看不明白的乱码，所以有必要设置一下文件格式。

2.用命令：*FILE FORMAT设定文件格式
将文件格式设为ASCII，那么结果文件我们就能够看明白了。

3.命令：*ELEMENT MATRIX OUTPUT的其他参数设定：
首先FREQUENCY要设定，否则会生成每一步的结果文件，这个文件会很大，举例来说，1800个单元，共50步，全部输出的单元刚度矩阵有600M，这是一个很大的数据，所以要把不需要的删除。并且在打开这个600M大的文件时发现，UntraEdit比记事本和写字板的打开速度快得多!假如共50步，我只想要最后一步的单元刚度矩阵，那么可以设定
FREQUENCY＝50，设定FREQUENCY＝0表示不输出，默认为每一步都输出。
其次可以设定结果文件名，设定时要使用OUTPUT FILE=USER DEFINED这个选项，将FILE NAME的值设定为想要的文件名，注意不要扩展名。
例如：
*Elset, Elset=All, Generate
1,   1800,     1
*File Format,Ascii
*Element Matrix Output,Elset=All,File Name=abc,Frequency=50,Output File=User Defined,Stiffness=Yes

计算完后可以在工作目录下找到文件名为ABC.mtx的单元刚度矩阵的结果文件，结果数据为十进制,大小仅为12M。