3D Solid FEM: N=ShapeFunctions(xi,eta,mu) - File Exchange

Code covered by the BSD License

Highlights from
3D Solid FEM

ElementMatrix.m
ExternalForces0D1D2D.m
N=ShapeFunctions(xi,eta,mu)
Plotmesh.m This function plots the mesh with or without the global nodenumbers. The
StrainMatrix.m Computes the six element strains, in case of the linear tetrahedral element, these strains are constant
StressRoutine.m
MainProgram.m
View all files

from 3D Solid FEM by Peter van Alem
3D FEM analysis using linear tetrahedral elements. Also a modal analysis is shown.

N=ShapeFunctions(xi,eta,mu)

function N=ShapeFunctions(xi,eta,mu)
% Isoparametric Shape Functions============================================
NN=[1-xi-eta-mu
    xi
    eta
    mu];
%==========================================================================     
N=[];
for i=1:numel(NN)
   Nsmall=[NN(i) 0     0
           0     NN(i) 0
           0     0     NN(i)];
   N=[N,Nsmall];
end

Highlights from 3D Solid FEM

Highlights from
3D Solid FEM