Introduction To Matrix Methods For Structural Analysis: Definitions Elements
Introduction To Matrix Methods For Structural Analysis: Definitions Elements
DEFINITIONS
Truss Element: A member subjected to axial forces only. Shown here is the type of
force and the associated axial deformation assigned on a truss element.
Beam Element: A member that resists transverse forces and moments. Shown here
are the shear and moment forces and associated transverse and rotational
displacements assigned on a beam element.
Frame Element: A member that resists axial forces in addition to moments and
transverse shear. Shown here are the forces and associated axial, transverse and
rotational displacements assigned on a frame element.
Nodes: Nodes are locations in the structure where elements are connected.
Node
Coordinates
Global Coordinates: Identify (with a unique number assigned) each DOF of the
structure in the global coordinates.
11
21
12
22
Definition 1: The stiffness coefficient Kij is the force at coordinate i when a unit
displacement is imposed at coordinate j, and the displacements at all other coordinates are
constrained to zero (fixed).
Consider the spring in the figure. The spring constant, k=10 lb/in is the force required to
cause a unit displacement: f = kx . Here, a force of 10 lb placed on node 2, causes a unit
displacement at node 2. At the same time, a reaction of 10 lbs develops at node 1 to keep
the system in static equilibrium and to ensure that the displacement of node 1 is zero.
If we consider Definition 1, we see that we defined K22 as a force at coordinate 2 when a unit
displacement is imposed at coordinate 2. Therefore, K22 has a value of 10. We also defined
K12 as the force that develops at coordinate 1 when a unit displacement is imposed at
coordinate 2. Therefore, K12 has the value of 10. These values are inserted into the
Stiffness Matrix for the element as follows:
Ke =
11
21
10
10
As seen, node 2 was used to fill-in column two of the stiffness matrix. Now, to fill-in column
one, we will work with the release of node 1 while fixing node 2.
A 10-lb force is applied at node 1, imposing a unit displacement on node 1. This force is
equivalent to K11. The reaction force of -10 lbs that develops at Node 1 is equivalent to K21.
Completing the remainder of the stiffness matrix, we obtain:
10
Ke =
10
10
10
Note that the stiffness of the bar is 10 k/in. The stiffness matrix, K, can be rewritten as:
Ke = 10
1 1
1
1
In general, the stiffness matrix for a single spring element of stiffness k can be written as:
1 1
Ke =
1
1
Now assume that the spring is loaded as shown below. We can solve for the unknown
displacement and force using the relationship {f}=[K]{x}, where [K] is given above.
200
F2
10 10
x1
10 10
The structural stiffness matrix is of size (4x4) since we are dealing with four degrees of
freedom. It is developed one column at a time. To find the first column:
Namely:
SF@node1: K11 10 20 = 0;
S@node2: K21 + 20 = 0;
SF@node3: K31 + 10 = 0;
K11 = 30
K21 = - 20
K31 = -10
7
30
20
10
0
Repeat the calculations after fixing all nodes except node 2 to find the second column.
Namely:
SF@1: K12 + 20 = 0;
K12 = -20
SF@2: K22 30 40 = 0; K22 = 50
SF@4: K42 + 30 = 0;
K32 = -30
Repeat the calculations after fixing all nodes except node 3 to find the third column.
Namely:
K13 = -10
K23 = 0
K33 = 10
K43 = 0
Repeat the calculations after fixing all nodes except node 4 to find the fourth column.
30 20
20 50
10
0
0
30
10
0
10
0
Repeat the calculations after fixing all nodes except node 4 to find the fourth column.
Namely:
K14 =
K24 =
K34 =
K44 =
0
-30
0
30
30
20
10
0
20 10
0
50
0
30
0
10
0
30
0
30
Now assume that the spring is loaded as shown below. We can solve for the unknown
displacement and force using the relationship {f}=[K]{x}, where [K] is given above.
Solution
Known:
x3 = 0
x4 = 0
F1 = 40 lb
F2 = 100 lb
Unknown:
40
30 20
20 50
100
=
10
0
3
4
0
30
x1 = ?
x2 = ?
F3 = ?
F4 = ?
10
0
1
0
30 2
10
0
0
0
30
0
This matrix equation cannot be solved as is stands. It must divided into parts where only
forces or displacements are the only unknowns. Also, since x3 and x4 are equal to zero,
columns 3 and 4 can be eliminated from the equations.
40
30
=
100
20
20 1
50 2
Once all the displacements are known, they can be used to solve for the remaining forces.
10
0
3.636
3
=
0
30 3.454
4
10
Here we will solve the same three-spring system by using the energy expressions.
where xi is the displacement of every node i. The potential energy P can be expressed as:
P=F1 x1+F2 x2+F3 x3+F4 x4
The total energy E in the system at any time is E(x)= U-P. In view of the equations above,
this relationship can be written as:
where
k2
k1 0
k1 + k 2
k 2 k 2 + k3 0
k 3
K=
k1
k1
0
0
k3
0
k3
0
Note: In order to understand how to work with matrices Please expand Equation 1 using the
K matrix above and show that it produces the energy of the structure U-P.
After the loading F is applied to the system, the system deforms. In this case, each DOF i
displaces an amount xi. We can solve for these displacements by solving for the vector x
that minimizes the energy of the system. Taking the derivative of the energy with the
respect to x and setting this derivative to zero:
we arrive at:
/ x [ xTKx xTF]=0
Kx=F
11
The solution of this equation for x gives the displacement that each DOF has undergone to
help the system achieve equilibrium. In other words, under loading, the system will
deform to a configuration that minimizes the total energy of the system. This is its
equilibrium position.
Note: Another way to arrive at the K x=F is to minimize the energy wrt x1, x2, x3, and x4.
Then write the four resulting equations in matrix form.
12
1) Repeat the lecture on energy minimization. Solve for all with your MATLAB (do not turn
this in)
2) Given that k1=10 and k2=30lb/in for the spring below.
Create the stiffness matrix of the structure by freeing one degree-of-freedom at a
time while constraining all others to zero.
Create the stiffness matrix for the structure by writing the kinetic and potential
energy expressions for the springs and minimizing the energy.
Apply boundary conditions and loading (assume a force F1= -30 is applied on node
one.) Solve for the equilibrium position.
Plot the energy of the system in MATLAB (note that you must apply boundary
conditions and loading before you solve for the energy.)
Show graphically that the vector x that minimizes the energy of the system is the
same as the one that solves the equations Kx=F
Use a MATLAB optimization function to find the minimum of the energy. Show that
minimization of energy gives you the same displacements as the solution of
Newtons Equations Kx=f.