03a Direct Stiffness Method Spring and Bar
03a Direct Stiffness Method Spring and Bar
ME751
Mechanical IV/II
Chapter 3
Direct Stiffness Method: Discrete
Finite Elements (Spring/Bar)
3.1 INTRODUCTION
The primary characteristics of a finite element are embodied in the
element stiffness matrix.
2
FINITE ELEMENT METHOD
3.2 SPRING AND BAR ELEMENTS
3.2.1 LINEAR SPRING AS FINITE ELEMENT
A linear elastic spring is a mechanical device capable of supporting
axial loading only and constructed such that, over a reasonable
operating range (meaning extension or compression beyond
undeformed length), the elongation or contraction of the spring is
directly proportional to the applied axial load.
The constant of proportionality between deformation and load is
referred to as the spring constant, spring rate, or spring stiffness,
generally denoted as k, and has units of force per unit length.
3
FINITE ELEMENT METHOD
As an elastic spring supports axial
loading only, we select an element
coordinate system (also known as a
local coordinate system) as an x axis
oriented along the length of the
spring, as shown
The ends of the spring are the nodes and the nodal displacements are
denoted by u1 and u2 and are shown in the positive sense.
5
FINITE ELEMENT METHOD
where
8
FINITE ELEMENT METHOD
System Assembly in Global Coordinates
Derivation of the element stiffness matrix for a spring element was
based on equilibrium conditions. The same procedure can be applied
to a connected system of spring elements by writing the equilibrium
equation for each node.
10
FINITE ELEMENT METHOD
Assuming the system of two spring elements to be in equilibrium, we
examine free-body diagrams of the springs individually and express
the equilibrium conditions for each spring as
The compatibility conditions state the physical fact that the springs
are connected at node 2, remain connected at node 2 after
deformation, and hence, must have the same nodal displacement at
node 2. FINITE ELEMENT METHOD
11
Substituting Equations (3.9) into Equations (3.8), we obtain
12
FINITE ELEMENT METHOD
The addition of Equations (3.11) and (3.12) yields
13
FINITE ELEMENT METHOD
which is of the form [K]{U} = {F}, similar to Equation (3.5).
However, Equation (3.15) represents the equations governing the
system composed of two connected spring elements. By direct
consideration of the equilibrium conditions, we obtain the system
stiffness matrix or global stiffness matrix [K] as
16
FINITE ELEMENT METHOD
EXAMPLE 3.1
Given: For the spring system shown in Figure E3.1 ,
k1 = 100 N/mm, k2 = 200 N/mm, k3 = 100 N/mm P = 500 N
Find: (a) the global stiffness matrix
(b) displacements of nodes 2 and 3
(c) the reaction forces at nodes 1 and 4.
Figure E3.1
17
FINITE ELEMENT METHOD
EXAMPLE 3.2
Three rigid bodies, 2, 3, and 4, are connected by four springs, as
shown in the Figure E3.2. A horizontal force of 1,000 N is applied
on Body 4. Find the displacements of the three bodies and the
forces (tensile/compressive) in the springs. What is the reaction at
the wall? Assume the bodies can undergo only translation in the
horizontal direction. The spring constants (N/mm) are k1 = 400,
k2 = 500, k3 = 500, k4 = 300.
Figure E3.2
18
FINITE ELEMENT METHOD
3.2 ELASTIC BAR AS FINITE ELEMENT
While the linear elastic spring serves to introduce the concept of the
stiffness matrix, the usefulness of such an element in finite element
analysis is rather limited. Certainly, springs are used in machinery in
many cases and the availability of a finite element representation of
a linear spring is quite useful in such cases. The spring element is
also often used to represent the elastic nature of supports for more
complicated systems.
19
FINITE ELEMENT METHOD
Consider an elastic bar of length L to which is affixed a uniaxial
coordinate system x with its origin arbitrarily placed at the left end.
This is the element coordinate system or reference frame.
20
FINITE ELEMENT METHOD
Using Equation (3.17), we obtain the equivalent spring constant of an
elastic bar as
21
FINITE ELEMENT METHOD
The associated axial force is then given by
Equation (3.21) is now used to relate the applied nodal forces f1 and
f2 to the nodal displacements u1 and u2 as
22
FINITE ELEMENT METHOD
EXAMPLE 3.3
For the bar assemblages shown in Figure E3.3, determine the
nodal displacements, the stresses in each element and the
reactions. Use the direct stiffness method for these problems.
Figure E3.3
23
FINITE ELEMENT METHOD