关于如何在matlab中导入并翻译Hypemesh导出的大型刚度矩阵txt文本

Posted 2021-09-09 文宇肃然

    刚度矩阵是有限元计算中最重要的矩阵，它是一个典型的对称稀疏矩阵，并且非零元素全部集中

前言

在对角线附近，稀疏矩阵就是零元素远大于非零元素的数目，所谓大型刚度矩阵。这里指稀疏矩阵中非零元素的数目有500多万。而这个稀疏矩阵是125万x125万≈1e12。如果不用稀疏矩阵的格式存储，这对matlab甚至其他计算软件来说都将会是内存溢出。这两篇推送将说明如何在matlab里如何导入并翻译Hypemesh导出的大型刚度矩阵txt文本。

    作为第一部分是如何导入的问题。在matlab里可以导入txt文本和excel文本。

对于一个600多万行的文本，复制到excel和word里不现实，只会让软件崩掉。所以最好在txt进行批处理。用在主页面用Import可以手动导入一些较少数据的txt文本和excel文本。但是对于一个600多万行的文本，不建议使用手动导入，这将花费大约30min才能完成导入。这里建议使用textread函数导入。但是使用data_K=textread(‘data.txt’)导入的时候，文本中如果有字母和其他符号时，将会报错。这就需要我们要么在textread里增加格式说明要么在txt里将除数字外的不用的东西批量删掉。

       其实很简单，在txt中使用替换功能就行，将* $  DIMG  KAAX等挨个批量替换成空格。