为什么一个简单的组装刚度阵需要这么长时间？

yuanwucai

大家请看我用matlab编写的有限元程序中组装刚度阵的一段：

其中的M是单元刚度阵在总刚中的位置。用matlab的profiler测试发现，第86行组装单刚所花的时间竟然占整个有限元程序的83%!

对于这一点我非常不理解，为什么一个简单的对号入座会花这么长时间？是不是我的程序存在问题？

请各位高手指点！！！

feynmand

简单说一点自己的感觉：
可能你这里M是变化的，这样对gk赋值的话可能就是在一个循环中改变了矩阵的大小，从而导致这一步的开销非常大。这里需要注意内存的使用。

yuanwucai

谢谢！正是如此，gK的大小是不断增大的。

caoer

接下来可能会更慢，用Matlab编fem不是一个好的选择，尤其是3D的

xyz999

试试在程序开始的地方，把gk初始化为一个比较大的矩阵，比如用gk=zeros(500,500)。

如果总刚矩阵很大，可以试试下面的语句（加在程序开始处）

1. 
   
2. dofpnd = 3; %每一个节点的自由度数目
   
3. dof0   = n*dofpnd; % n是节点的个数
   
4. nn_z = floor(dof0*dof0*0.1); % 这儿假设刚度矩阵是稀疏矩阵。0.1是我使用的值，你可以设的大一点或小一点
   
5. gk = spalloc(dof0, dof0,nn_z);

复制代码

如果matlab报错（memory不够），或还是很慢，只好建议你使用其它语言了。