Accompany 发表于 2015-11-5 17:24

ABAQUS输出单元刚度矩阵

蔡和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。
页: [1]
查看完整版本: ABAQUS输出单元刚度矩阵