Fast Global Stiffness Matrix Assembly
function K=Fast_Matrix_Assembly(Elements)
Input:
---------
Elements: a structure contains Elements{i}.K and Elements{i}.DOFs (Stiffness and Degrees of freedoms of the ith element).
output:
---------
K: Sparse global stiffness matrix
References
-----------------
Cuvelier, François, Caroline Japhet, and Gilles Scarella. "An efficient way to perform the assembly of finite element matrices in Matlab and Octave." arXiv preprint arXiv:1305.3122 (2013).
Cite As
Ayad Al-Rumaithi (2021). Fast Global Stiffness Matrix Assembly (https://www.mathworks.com/matlabcentral/fileexchange/70447-fast-global-stiffness-matrix-assembly), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (11)
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxTags
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
hello
I am trying to assemble global stiffness matrix for CST element I want to ask if this file is useful for it plus, I don't understand how to use this m.file for my codes please help me
thanks.
It seems you didn’t add the effect of supports to the stiffness matrix. You should remove the rows and columns that corresponding to the supports degree of freedom
I tried this code but I always get the warning 'Matrix is singular to working precision'. I am unable to find where am I going wrong. Can you please give any suggestions?
Dear Anilcan Taner,
Elements{i}.K should be the stiffness matrix of ith element .Elements{i}.DOFs degrees of freedom vector of the ith element.
in your example it could be:
Elements{1}.K=KK{1};
Elements{1}.DOFs=[5 1 6];
Elements{2}.K=KK{2};
Elements{2}.DOFs=[6 1 2];
.....
K=Fast_Matrix_Assembly(Elements)
I am having problem with inputs for example;
i have 3x3 matrices named as KK{i} which i=1:9
but I am a bit lost with putting it to the code
Are you able to give an example of an input? I am trying to run this for a 3D frame model. Thanks in advance
Yes. It can be used with any type of elements
can this be used for 3-D elements?
Thank you
Very useful. Very fast