Fast FEM assembly: nodal elements
This code demonstrates vectorization concepts from the paper:
Talal Rahman and Jan Valdman: Fast MATLAB assembly of FEM matrices in 2D and 3D: nodal elements, Applied Mathematics and Computation 219, 7151–7158 (2013).
A link to the paper can be found at the author web page located at http://sites.google.com/site/janvaldman/publications
Please cite the paper if you find the code useful.
Available scripts from the paper:
"start_compute_areas"
"start_assembly_P1_2D"
"start_assembly_P1_3D"
"start_assembly_P1_2D_elasticity"
"start_assembly_P1_3D_elasticity"
Additional script (provided by Tien Dat Ngo from EFPL Lausanne):
"start_assembly_Q1_3D_elasticity"
Cite As
Jan Valdman (2021). Fast FEM assembly: nodal elements (https://www.mathworks.com/matlabcentral/fileexchange/27826-fast-fem-assembly-nodal-elements), MATLAB Central File Exchange. Retrieved .
Comments and Ratings (8)
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
Tags
Acknowledgements
Inspired: Element-based Solvers, Fast FEM assembly: edge elements, Evolutionary model for shape memory alloys, Hyper elasticity with a non-penetration condition , Implementation of C1 FEM, Simulation of von Kármán viscoelastic plates
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.
fast_fem_assembly
fast_fem_assembly/library_vectorization
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)
Is the function "phider" not uploaded ?
coeffs denotes a scalar coefficient in corresponding bilinear forms:
a(u,v) = (coeffs * grad u, grad v) (in the stiffness matrix function),
a(u,v) = (coeffs *u, v) (in the mass matrix function),
where u,v are P1 functions. Ignore it in the input of both functions to have coeffs=1.
It is straightforward to modify FEM assemblies to work with rectangular elements, but it is not done in this version. Read the paper and try it yourself :).
What's the 'coeffs' mean in function 'stifness_matrixP1' and 'mass_matrixP1'?
can this create a mass matrix for a rectangular (2 dimensional) mesh?