MANE 4240 & CIVL 4240
Introduction to Finite Elements
Prof. Suvranu De
Principles of minimum
potential energy and
Rayleigh-Ritz
Reading assignment:
Section 2.6 + Lecture notes
Summary:
Potential ener
%this script will solve reactions and displacements for a space truss.
d = zeros([3*nodes,1]);
F = zeros([3*nodes,1]);
disp(' ')
disp(' ')
disp('This part of the program will solve for gloabl reaction
%Practice with symbols
syms L E I A
k = E*I/L*[1,-1;-1,1];
j = L/L;
r = L^2;
f = L*L;
g = L^2/L;
syms b1 b2
A = sym('a',[2 2]);
A(1,1)=A(1,1)-A(1,1)+L;
i=1;
while i<=2
while j<=2
if A(i,j) = 2;
%Alexander Tadla
%CEE2330 Homework No. 7
%This program solves Problem 4.22 for Homework No. 7
%This program solves a displacements, slopes, member end forces, and reactions
tart script
clear
clc
%Defi
%Practice with symbols
syms L E I A
k = E*I/L*[1,-1;-1,1];
j = L/L;
r = L^2;
f = L*L;
g = L^2/L;
syms b1 b2
A = sym('O',[2 2]);
A(1,1)=A(1,1)-A(1,1)+L;
i=1;
while i<=2
j=1;
while j<=2
if size
%Alexander Tadla
%CEE2330 Homework No. 7
%This program solves Problem 4.22 for Homework No. 7
%This program determines the midspan deflection for the beam and the
%reactions.
tart script
clear
%Defi
%This program will create a global stiffness matrix
%And Calculate reactions and bar forces
%3D problem only with 3 DoF at each node
clear
disp(' ')
disp(' ')
0efine system
members = input('How many m
%Alexander Tadla
%CEE2330 Homework No. 7
%This program solves Problem 4.22 for Homework No. 7
%This program determines the midspan deflection for the beam and the
%reactions.
tart script
clear
%Defi
%Alexander Tadla
%CEE2330 Homework No. 7
%This program solves Problem 4.22 for Homework No. 7
%This program solves a displacements, slopes, member end forces, and reactions
tart script
clear
%Define G
%Alexander Tadla
%CEE2330 Homework No. 7
%This program solves Problem 4.34 for Homework No. 7
%This program solves a displacements, slopes, and reactions
clear
0efine givens
syms I
w=20000; %kN/m
L=6;
%This program will create a global stiffness matrix
%And Calculate reactions and bar forces
%3D problem only with 3 DoF at each node
clear
disp(' ')
disp(' ')
0efine system
members = input('How many m
%This program will solve reactions and displacements for a spring system,
%using the direct stiffness method, with known forces
clear
tart script
disp('Label the elements and nodes. Label the nodes wi
%this script will solve reactions and displacements for a space truss.
d = zeros([3*nodes,1]);
F = zeros([3*nodes,1]);
disp(' ')
disp(' ')
disp('This part of the program will solve for gloabl reaction
%this script will solve for member forces in a space truss
counter=1;
while counter<=members
cx = eval(['cos_' num2str(counter) '(1,1)']);
cy = eval(['cos_' num2str(counter) '(2,1)']);
cz = eval(
%this script will solve for member forces in a space truss
counter=1;
while counter<=members
cx = eval(['cos_' num2str(counter) '(1,1)']);
cy = eval(['cos_' num2str(counter) '(2,1)']);
cz = eval(
%CEE2330 Homework Set No. 3
%February 11, 2013
%This program will create a global stiffness matrix
%And Calculate reactions and bar forces
%2D problem only with 2 DoF at each node
clear
disp(' ')
disp
%this part of the program will find member forces
counter=1;
while counter<=members
c=cosd(theta_matrix(1,counter);
s=sind(theta_matrix(1,counter);
T = [c,s,0,0;-s,c,0,0;0,0,c,s;0,0,-s,c];
%this part of the program will find member forces
counter=1;
while counter<=members
c=cosd(theta_matrix(1,counter);
s=sind(theta_matrix(1,counter);
T = [c,s,0,0;-s,c,0,0;0,0,c,s;0,0,-s,c];
%CEE2330 Homework Set No. 3
%February 11, 2013
%This program will create a global stiffness matrix
%And Calculate reactions and bar forces
%2D problem only with 2 DoF at each node
clear
disp(' ')
disp
%Alexander Tadla
%CEE2330 Homework No. 6
%This program solves Number 3 for Homework No. 6
%This program solves a displacements, member end forces, and reactions for
%an indeterminate beam with a sprin
0ind reactions and find displacements of a 2D truss
displacement_matrix = zeros([2*nodes,1]);
force_matrix = zeros([2*nodes,1]);
disp(' ')
disp(' ')
disp('This part of the program will solve for gloab
%Alexander Tadla
%CEE2330 Homework No. 6
%This program solves Number 2 for Homework No. 6
%This program solves dispalcements, member end forces, and reactions for an
%indeterminate beam with fixed end
%Alexander Tadla
%CEE2330 Homework No. 6
%This program solves Number 1 for Homework No. 6
%This program solves a displacements, member end forces, and reactions for
%an indeterminate beam with a sprin