Scribd
Upload
138 views

MATLAB Codes For Finite Element Analysis

This MATLAB code performs a finite element analysis on a 2D frame structure with 3 elements. It defines material properties, generates node coordinates and element connectivities, applies forces, forms the stiffness matrix, applies boundary conditions, solves for displacements, and outputs/plots the deformed mesh.

Uploaded by

tekili
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
138 views

MATLAB Codes For Finite Element Analysis

This MATLAB code performs a finite element analysis on a 2D frame structure with 3 elements. It defines material properties, generates node coordinates and element connectivities, applies forces, forms the stiffness matrix, applies boundary conditions, solves for displacements, and outputs/plots the deformed mesh.

Uploaded by

tekili
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

% MATLAB codes for Finite Element Analysis

% problem10.m
% antonio ferreira 2008

% clear memory
clear all

% E; modulus of elasticity
% I: second moment of area
% L: length of bar
E=210000; A=100; I=2e8; EA=E*A; EI=E*I;

% generation of coordinates and connectivities


numberElements=3;
p1=3000*(1+cos(pi/4));
nodeCoordinates=[0 3000;3000 3000;p1 0;p1+3000 0];
xx=nodeCoordinates;
for i=1:numberElements;
elementNodes(i,1)=i;
elementNodes(i,2)=i+1;
end
numberNodes=size(nodeCoordinates,1);
xx=nodeCoordinates(:,1);
yy=nodeCoordinates(:,2);

% for structure:
% displacements: displacement vector
% force : force vector
% stiffness: stiffness matrix
% GDof: global number of degrees of freedom
GDof=3*numberNodes;
U=zeros(GDof,1);
force=zeros(GDof,1);

%force vector
force(6)=-10000;
force(7)=-10000;
force(10)=-5e6;
force(11)=5e6;

% stiffness matrix
[stiffness]=...
formStiffness2Dframe(GDof,numberElements,...
elementNodes,numberNodes,xx,yy,EI,EA);

% boundary conditions and solution


prescribedDof=[1 4 5 8 9 12]';

% solution
displacements=solution(GDof,prescribedDof,stiffness,force);

% output displacements/reactions
outputDisplacementsReactions(displacements,stiffness,...
GDof,prescribedDof)
% drawing undeformed and deformed meshes
U=displacements;
clf
drawingMesh(nodeCoordinates+500*[U(1:numberNodes)...
U(numberNodes+1:2*numberNodes)],elementNodes,'L2','k.-');
drawingMesh(nodeCoordinates,elementNodes,'L2','k--');

You might also like