Matrix Analysis of The Beam
Matrix Analysis of The Beam
TE 52S
EARTHQUAKE ENGINEERING
PROJECT #2
Submitted by:
John Peter L. Garay
BSCE3-LDP
Submitted to:
Engr. Robert Lucena
Professor
TABLE OF CONTENTS
I. Introduction
II. Formation of Local Stiffness Matrix
III. Formation of Global Stiffness Matrix
IV. Computation of Displacements on Nodes
V. Computation of Moments at Beam Supports
VI. Computation of the Reactions
I. INTRODUCTION
This project is mainly focuses on a structure where an axial forces, torsion, and beam bent are
present and we should compute for the values of displacements on each node, moments at
supports and the reaction.
To start the project, we must need to make first a visual figure of the given structure so that we
can easily identify the forces acting on it.
Figure 1: A sample figure of the given structure with forces on different nodes.
The sample figure above is a part of a 40 nodes with 39 member structure. Each node has a force
where we can put label on it by abiding on conditions instructed to us using our student I.D.
number.
The formation of the local stiffness matrix per member of the structure comes from the
bisymmetrical framework element stiffness matrix and can be solved with variables listed below:
By that, you can now compute and assemble for the local stiffness matrix on each member. This
is a 12x12 matrix.
One of the hardest part and it can consume more time on this project is the assembling of the
global stiffness matrix. The global stiffness matrix is just a combination of the entire local
stiffness matrix on each member that we computed (shown on my computations spread sheets)
where we need to be careful on putting values onto the matrix.
Figure 4: Global stiffness matrix from my Microsoft Excel spread sheet
Now, I will tackle on how I’ve computed the displacement on each node using the global
stiffness matrix and the Microsoft Excel. First, we need to rearrange the global stiffness matrix. I
omitted all of the out-of-plane shear and bending forces. After that, we need to identify the forces
that have a displacement. Rearranged the global stiffness matrix with four parts which are the
Kff, Ksf, Kfs, and Kss. The Kff part of the rearranged matrix is the one we need to solve for the
displacement on the nodes.
Kff Kfs
Ksf Kss
Figure 5: A representation of a rearranged global stiffness matrix
After the rearrangement of the global stiffness matrix, we now have the Kff which we need to
solve the displacement on nodes. In order to solve for that, we need to get the inverse of the Kff
matrix. Using a Microsoft Excel, just use the MINVERSE function to get the inverse of it.
By using another function on Microsoft Excel, we should use the MMULT to multiply the
inverse of the Kff on the forces acting on nodes.
{=MMULT(I7:K9,A9:K10)}
{=MINVERSE(I7:K9)}
Figure 6: Examples of the functions used in Microsoft Excel
FORCES DISPLACE
ACTING ON MENTS ON
EACH EACH
NODE NODE
The computation of the moments on beam supports is quite similar with the computation of the
displacement. But the first thing to do is that we need to have the values of the displacement on
each node so that we can proceed on the computation of the moments.
Some of the beam supports has already a value of the moment but also some supports do not has
already of it so we compute for it. In order to get the values, rearranged again the global stiffness
matrix and do not remove the forces (Mz) with their corresponding degree of freedom.
Using an MMULT function again, we can now use it to multiply the matrix of the forces with
unknown moments and to their corresponding displacement.
{=MMULT(I7:K9,A9:K10)}
Figure 8: Examples of the functions used in Microsoft Excel
CORRESPONDING
MOMENTS ON
DISPLACEMENT
BEAM SUPPORTS
OF THE FORCES
In order to compute for the reaction, we need to go back once again on the rearranged global
stiffness matrix. Earlier, we used the Kff part of it to compute for the displacements. Now, we
will use the Ksf part of it in order to obtain the reactions on each support and nodes.
Using once again the MMULT function in the Microsoft Excel to multiply the Ksf into the
obtained displacement.
{=MMULT(I7:K9,A9:K10)}
Figure 10: Examples of the functions used in Microsoft Excel
KSF
COMPUTED REACTIONS
DISPLACEM AT
ENTS SUPPORTS
Figure 11: A visual representation on how we compute for the reaction at supports.