Code covered by the BSD License  

Highlights from
Fast assembly of stiffness and matrices in finite element method

Be the first to rate this file! 46 Downloads (last 30 days) File Size: 26.9 KB File ID: #27826
image thumbnail

Fast assembly of stiffness and matrices in finite element method

by Jan Valdman

 

03 Jun 2010 (Updated 19 Mar 2013)

vectorization of FEM matrices assembly

| Watch this File

File Information
Description

This code demonstrates vectorization concepts from a 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)

Original technical report from 2011 is available at
http://www.mis.mpg.de/publications/other-series/tr/report-1111.html

To run the code, call

"start_assembly_P1_3D" or
"start_assembly_P1_2D"

in the Matlab window (for 3D or 2D cases)
 
Considering a given 3D or 2D triangulation, mass and stiffness matrices corresponding to linear nodal elements are generated in an efficient vectorized and fast way.

Recently, functionality was extended to generate mass and stiffness matrices for a (vector) linear elasticity problem.

To run the code, call

"start_assembly_P1_3D_elasticity" or
"start_assembly_P1_2D_elasticity"

MATLAB release MATLAB 7.9 (R2009b)
Tags for This File  
Everyone's Tags
fem(2), fem matrix assembly, linear elasticity(2), mass matrix, stiffness matrix, vectorization
Tags I've Applied
Add New Tags Please login to tag files.
Please login to add a comment or rating.
Updates
10 Jul 2010

2D functionality added

16 Sep 2010

added start_calculate_areas

17 Feb 2011

tested for both linux and windows

02 Mar 2011

report link add, extended to work for linear elasticity

19 Mar 2013

Paper information added

Contact us