Vectorized Matlab codes for linear two-dimensional elasticity

Article type: Research Article

Authors: Koko, Jonas

Affiliations: LIMOS, Université Blaise Pascal, CNRS UMR 6158, and ISIMA, Campus des Cézeaux, BP 10125, 63173 Aubiére Cedex, France. E-mail: [email protected]

Abstract: A vectorized Matlab implementation for the linear finite element is provided for the two-dimensional linear elasticity with mixed boundary conditions. Vectorization means that there is no loop over triangles. Numerical experiments show that our implementation is more efficient than the standard implementation with a loop over all triangles.

Keywords: Finite element method, elasticity, Matlab

Journal: Scientific Programming, vol. 15, no. 3, pp. 157-172, 2007