MATLAB Functions For Finite Element Analysis
MATLAB Functions For Finite Element Analysis
Introduction
4. Applying the boundary conditions – like supports and applied loads and dis-
placements. In this book this step will be performed manually.
5. Solving the equations – this will be done by partitioning the global stiffness
matrix and then solving the resulting equations using Gaussian elimination.
In this book the partitioning process will be performed manually while the
solution part will be performed using MATLAB with Gaussian elimination.
6. Post-processing – to obtain additional information like the reactions and el-
ement forces and stresses. In this book this step will be performed using
MATLAB.
It is seen from the above steps that the solution process involves using a com-
bination of MATLAB and some limited manual operations. The manual operations
employed are very simple dealing only with discretization (step 1), applying bound-
ary conditions (step 4) and partitioning the global stiffness matrix (part of step 5). It
can be seen that all the tedious, lengthy and repetitive calculations will be performed
using MATLAB.
1.2
MATLAB Functions for Finite Element Analysis
The CD-ROM accompanying this book includes 84 MATLAB functions (M-files)
specifically written by the author to be used for finite element analysis with this book.
They comprise what may be called the MATLAB Finite Element Toolbox. These
functions have been tested with version 7 of MATLAB and should work with any
later versions. The following is a listing of all the functions available on the CD-ROM.
The reader can refer to each chapter for specific usage details.
SpringElementStiffness(k)
SpringAssemble(K, k, i, j)
SpringElementForces(k, u)
LinearBarElementStiffness(E, A, L)
LinearBarAssemble(K, k, i, j)
LinearBarElementForces(k, u)
LinearBarElementStresses(k, u, A)
QuadraticBarElementStiffness(E, A, L)
QuadraticBarAssemble(K, k, i, j, m)
QuadraticBarElementForces(k, u)
QuadraticBarElementStresses(k, u, A)
PlaneTrussElementLength(x1 , y1 , x2 , y2 )
PlaneTrussElementStiffness(E, A, L, theta)
PlaneTrussAssemble(K, k, i, j)
PlaneTrussElementForce(E, A, L, theta, u)