DRAFT - 1

CHAPTER 4




Beam Element





4.1 Introduction

Beam element is a very versatile line-element, it has six degrees of freedom at each node,
which include, translations and rotations along the x, y, and z directions, respectively.
Figure 4.1 shows the positive directions of these displacements.


y
v
jy



jy
v
jx

i x

jz
j
jx



z v
jz




Figure 4.1 Beam Element with six degrees of freedom at each node



Beam element is employed to simulate a slender structure that has an uniform cross
section. The element is unsuitable for structures that have complex geometry, holes, and
points of stress concentration.

The stiffness constant of a beam element is derived by combining the stiffness constants
of a beam under pure bending, a truss element, and a torsion bar. Thus, a beam element
can represent a beam in bending, a truss element, and a torsion bar. In FEA its a
common practice to use beam elements to represent all or any of these three loads.

We will derive the element stiffness equation for a beam element by first deriving the
stiffness equation of a beam in bending, and then superimposing the stiffness of a truss
and a torsion bar element.
DRAFT 1
Chapter 4 Beam Element
4.1 Derivation of a Stiffness Equation for a Beam Element
Under Pure Bending

A beam, such as, a cantilever beam, under pure bending (without axial loads or torsional
loads), has two-degrees of freedom at any point, transverse deflection v and rotation , as
shown in Figure 4.2.
F





v


Figure 4.2 Cantilever Beam with its DOF, v and



A beam element has a total of four degrees of freedom, two at each node. Since there are
four degrees of freedom, the size of the stiffness matrix of a beam element has the size 4
x 4.

We will derive the stiffness matrix equation using a simple method, known as Stiffness
Influence Coefficient Method. In this procedure, a relationship between force and the
coefficients that influence stiffness is established. For a beam element, these coefficient
consist of: the modulus of elasticity, moment of inertia, and length of the element. For a
two-node beam element, there are two deflections and two rotations, namely, v
1
,
1
, v
2
,
and
2
. Force and influence coefficient relationship is established by setting each of the
four deflection values to unity, with the remaining deflection values equal to zero. The
procedure follows.

Consider a beam element, loaded in such a way that it has the deflection values: v
i
= 1,
i

= 0, v
j
= 0,
j
= 0



i j
v
i
,
i
v
j
,
j



Figure 4.3 Beam Element


The above deflections can be produced by a combination of load conditions, shown in
figure 4.4.
ME 273 Lecture Notes by R. B. Agarwal 4-2
DRAFT 1
Chapter 4 Beam Element

The following deflection relationships for loading of Figures 4.4 (a) and (b) can be found
in any Machine Design Handbook, and is given as,


v
max


v
max
= (FL
3
)/(3EI)
y
= - (FL
2
)/(2EI)
i L j x

F
(a)



y
M
i

L M
j
v
max
= - (ML
2
)/(2EI)
i j
v
max
x = (ML)/(EI)
(b)

Figure 4.4


Applying these relationships to the beam of Figure 4.3, we get,

1 = v
i
= (v
i
)
F
+ (v
i
)
M


1 = v
i
= (F
i
L
3
)/3 EI - (M
i
L
2
)/2EI (4.1)

and = 0 = ()
F
+ ()
M

0 = - (F
i
L
2
)/2EI + (M
i
L)/EI (4.2)

Solving Equations (4.1) and (4.2), we get,

F
i
= (12EI)/L
3
(A)

F
j
= - F
i
= -(12EI)/L
3
(B)

M
i
= (6EI)/L
2
(C)


From Figure 4. 4 (a) and (b),
ME 273 Lecture Notes by R. B. Agarwal 4-3
DRAFT 1
Chapter 4 Beam Element
M
j
= F
i
L - M
i

= (12EI)/L
2
= (6EI)/ L
2

= (6EI)/ L
2
(D)


Writing equations (A) through (D) in a matrix form we get,


F
i
(12EI)/L
3
1 (12EI)/ L
3
0 0 0 1


M
i
(6EI)/ L
2
1

(6EI)/ L
2
0 0 0 0
= =
F
j
-(12EI)/ L
3
1

-(12EI)/ L
3
0 0 0 0

M
j
(6EI)/ L
2
1

(6EI)/ L
2
0 0 0 0



Using a similar procedure and setting the following deflection values:

v
i
= 0,
i
= 1, v
j
= 0,
j
= 0, we get,



F
i
(6EI)/L
2
1 0 (6EI)/ L
2
0 0 0


M
i
(4EI)/ L 1

0

(4EI)/ L

0 0 1
= = (4.6)
F
j
-(6EI)/ L
2
1

0 -(6EI)/ L
2
0 0 0

M
j
(2EI)/ L

1 0 (2EI)/ L

0 0 0



Similarly, setting v
j
= 1 and ,
j
= 1, respectively, and keeping all other deflection values
to zero, we get the final matrix as,



F
i
(12EI)/L
3
(6EI)/ L
2
-(12EI)/ L
3
(6EI)/ L
2
1


M
i

(6EI)/ L2 (4EI)/ L

-(6EI)/ L
2
(2EI)/ L 1
= (4.7)
F
j

-(12EI)/ L
3
-(6EI)/ L
2
(12EI)/L
3
-(6EI)/ L
2
1

M
j
(6EI)/ L
2
(2EI)/ L -(6EI)/ L
2
(4EI)/ L

1


ME 273 Lecture Notes by R. B. Agarwal 4-4
DRAFT 1
Chapter 4 Beam Element



Note that, the first term on the RHS of the above equation is the stiffness matrix and the
second term is the deflection. In the case where deflections are other than unity, the
above equation will provide an element equation for a beam (in bending), which can be
written as,



F
i
(12EI)/L
3
(6EI)/ L
2
-(12EI)/ L
3
(6EI)/ L
2
v
i


M
i

(6EI)/ L2 (4EI)/ L

-(6EI)/ L
2
(2EI)/ L
i

= (4.7)
F
j

-(12EI)/ L
3
-(6EI)/ L
2
(12EI)/L
3
-(6EI)/ L
2
v
j


M
j
(6EI)/ L
2
(2EI)/ L -(6EI)/ L
2
(4EI)/ L

j





Where F
i
, M
i
, F
j
, M
j
are the loads corresponding to the deflections v
i
,
i
, v
j
,
j
.

Equation (4.7) is the equation of a beam element, which is under pure bending load (no
axial or torsion loads). The stiffness matrix is a 4 x 4, symmetric matrix. Using this
equation, we can solve problems in which several beam elements are connected in an uni-
axial direction. The assembly procedure is identical to the truss elements. However, if the
beam elements are oriented in more than one direction, we will have to first transform the
above equation (4.7) in to a global stiffness matrix equation (analogues to the procedure
used for truss elements).

For a beam element, transformation of a local stiffness matrix into a global equation
involves very complex trigonometric relations, and therefore, we will defer the
derivations at this time. However, Equation (4.7) can be used for solving a beam
problem, loaded under bending loads. In order to understand the application of this
equation, we will apply it to solve some statically indeterminate problems.











ME 273 Lecture Notes by R. B. Agarwal 4-5
DRAFT 1
Chapter 4 Beam Element

Example 1

For the beam shown, determine the displacements and slopes at the nodes, forces in each
element, and reactions at the supports.
5 ft 5 ft 100 lb
E = 1.4 x 10
6
psi, I = 2.4 in
4



K = 200 lb/in




Solution

The beam structure is descritized into three elements and 4-nodes, as shown.


[1] [2]
3
1 2

[3]

4


First, we will find the element stiffness matrix for each element, next we will assemble
the stiffness matrices, apply the boundary conditions, and finally, solve for node
deflection. Internal forces and reactions are calculated by back-substituting the
deflections in the structural equation.

[1]
Element 1
1 2

EI/L
3
= (1.4 x 10
6
) x (2.4)/(5x12)
3
= 15.55

The general equation of a stiffness matrix is given as,






ME 273 Lecture Notes by R. B. Agarwal 4-6
DRAFT 1
Chapter 4 Beam Element

12 6L -12 6L v
1


6L 4L
2
-6L 2 L
2

1

[K
e
]
(1)
= (EI/L
3
)
-12 -6L 12 -6L v
2


6L 2 L
2
-6L 4 L
2

2





[2]
Element 2 2 3



12 6L -12 6L v
2


6L 4L
2
-6L 2 L
2

2

[K
e
]
(1)
= (EI/L
3
)
-12 -6L 12 -6L v
3


6L 2 L
2
-6L 4 L
2

3






3
Element 3
[3]


4


[K
e
]
(3)
= K -K v
3

-K K v
4







ME 273 Lecture Notes by R. B. Agarwal 4-7
DRAFT 1
Chapter 4 Beam Element

To get the global stiffness matrix, we will use the same procedure used for assembling
truss element stiffness equations. In terms of E, L, and I the assembled global stiffness
matrix is,


v
1

1
v
2

2
v
3

3
v
4


v
1
12 6L -12 6L 0 0 0


1
4L
2
-6L 2 L
2
0 0 0

v
2
24 0 -12 6L 0
x (EI) /(L
3
)

2
8L
2
-6L 2L
2
0

v
3
12 +K

-6L - K

3
4L
2
0

v
4
SYMMETRY

K



Where K = (K) x [L
3
/ (EI)]

Our next step is to write the structural equation; however, we can reduce the size of the
stiffness matrix by applying the given boundary conditions:

v
1
=
1
= 0 node 1 is fixed

v
2
= 0

node 2 has no vertical deflection, but its free to rotate.

V
4
= 0 node 4 is fixed.


The reduced stiffness matrix is


8L
2
-6L 2L
2
K
G
= EI / (L
3
) -6L 12+K -6L
2L
2
-6L 4L
2


Substituting the values of E, L, and I the structural equation can be written as,

ME 273 Lecture Notes by R. B. Agarwal 4-8
DRAFT 1
Chapter 4 Beam Element

0 1152 -72 288
2
-100 = (15.55) -72 16.11 -72 v
3
0 288 -72 576
3




2
= - 0.0032 rad
Solving, we get v
3
= - 0.4412 in

3
= -0.0095 rad




ME 273 Lecture Notes by R. B. Agarwal 4-9
DRAFT 1
Chapter 4 Beam Element
4.2 Arbitrarily Oriented 2-D Beam Element

The stiffness equation for an arbitrarily oriented beam element can be derived with a
procedure similar to the truss element.

y d
2y

x
y

2

d
1y
d
1y


1
d
1x
x


d
1y
= d
1y
cos - d
1x
sin = d
1y
c - d
1x
s

d
2y
= d
1y
cos d
2x
sin = d
2y
c d
2x
s

and
1
=
1
,
2
=
2

Note: The underscored terms represent local coordinate values. Thus, x and y are local
coordinates and x and y are global coordinates.

The above equations can be written in a matrix form,

d
1x

d
1y
-s c 0 0 0 0 d
1y

1
= 0 0 1 0 0 0
1

d
2y
0 0 0 -s c 0 d
2x

2
0 0 0 0 0 1 d
2y


2


-s c 0 0 0 0
Let T

= 0 0 1 0 0 0
0 0 0 -s c 0

0 0 0 0 0 1 , the transformation matrix.



Thus, {d} = [T] {d}
Global
Local

Note that angle

is independent of the coordinate systems, and
1 =

1
,
2
=
2

ME 273 Lecture Notes by R. B. Agarwal 4-10
DRAFT 1
Chapter 4 Beam Element

As derived in the case of the truss element, relationship between local and global stiffness
matrices is given as

[k
g
] = [T] [k] [T]

Where, [k
g
] = Global stiffness matrix of an element
[T] = Transformation matrix
[k] = Local stiffness matrix of the element


Substituting the values of [T] and [k], we get the global equation of a beam element
oriented arbitrarily at an angle as,


12S
2
-12SC -6LS -12S
2
-12SC -6LS
12C
2
6LC 12SC -12C
2
6LC
k = EI/L 4L
2
6LS -6LC 2L
2

12S
2
-12SC 6LS
Symmetry 12C
2
4L
2



This is the equation of a beam element (without axial or torsional load, and oriented at an
angle .

Also, S = sin , C = cos in the above equation.


4.3 Beam Element with Combined Bending and Axial loads

First, we will derive the stiffness matrix in local coordinates and then convert it in to
global coordinates.

4.3.1 Stiffness matrix of a beam element with bending and axial loads in
local coordinates

The stiffness equation for the combined bending and axial load can be written by
superimposing the axial stiffness terms over the bending stiffness.

For axial loading, the structural equation is,


f
1x
1 -1 d
1x

= AE/L
3


f
2x
-1 1 d
2x


ME 273 Lecture Notes by R. B. Agarwal 4-11
DRAFT 1
Chapter 4 Beam Element

And for bending loading, the structural equation is,



f
1y
12 6L -12 6L d
1y


m
1
6L 4L
2
-6L 2L
2

1
= AE/L
3

f
2y
-12 -6L 12 -6L d
2y


m
2
6L 2L
2
-6L 4L
2

2





Therefore, the combined loading equation is


f
1x
C
1
0 0 - C
1
0 0 d
1x


f
1y
0 12 C
2
6C
2
L 0 -12 C
2
6C
2
L d
1y


m
1
0 6 C
2
L 4C
2
L
2
0 -6C
2
L 2C
2
L
2

1

=
f
2x
-C
1
0 0 C
1
0 0 d
2x


f
2y
0 -12 C
2
-6C
2
L 0 12 C
2
-6C
2
L d
2y


m
2
0 6 C
2
L 2C
2
L
2
0 -6C
2
L 4C
2
L
2

2






And, [
(
(
(
(
(
(
(

=
2
2 2
2
2 2
2 2 2 2
1 1
2
2 2
2
2 2
2 2 2 2
1 1
4 6 0 2 6 0
6 12 0 6 12 0
0 0 0 0
2 6 0 4 6 0
6 12 0 6 12 0
0 0 0 0
]
L C L C L C L C
L C C L C C
C C
L C L C L C L C
L C C L C C
C C
k


Where, C
1
= AE/L, and C
1
= EI/L
3
ME 273 Lecture Notes by R. B. Agarwal 4-12
DRAFT 1
Chapter 4 Beam Element


4.3.2 Transformation matrix for combined Bending and Axial loading.

For the axial loading, the relationship between the local and global coordinates was
derived earlier, as

S d C d d
S d C d
d d d
y x
x
y x
y x
x
2 2
2
^
1 1
1 1
1
^
sin cos
+ =
+ =
+ =


Also, for bending load, derived previously,



=
=
=
=
^
2
2 2
2
^
^
1
1 1
1
^
S d C d d
S d C d d
x y
y
x y
y


Therefore, the relationship for the combined bending and axial loading can be written as

(
(
(
(
(
(
(

2
2
2
1
1
1
^
2
^
2
^
2
^
1
^
1
^
1
1 0 0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 1 0 0
0 0 0 0
0 0 0 0

y
x
y
x
y
x
y
x
d
d
d
d
C S
S C
C S
S C
d
d
d
d

Or,
| |{ } d T d =
)
`

^
ME 273 Lecture Notes by R. B. Agarwal 4-13
DRAFT 1
Chapter 4 Beam Element
Where,
(
(
(
(
(
(
(

=
1 0 0 0 0 0
0 0 0 0
0 0 0 0
0 0 0 1 0 0
0 0 0 0
0 0 0 0
] [
C S
S C
C S
S C
T



4.3.3 2-D Beam Element Equation for Combined Loading Axial and
Bending at an Arbitrary Orientation

Substituting the values of and [T] into the equation | | , we get
(

^
K ] ][ [ ] [
^
T K T K
T
=

(
(
(
(
(
(
(
(
(
(
(
(

+
|
.
|

\
|
+

|
.
|

\
|
+
|
.
|

\
|
+

|
.
|

\
|

|
.
|

\
|
+
|
.
|

\
|
+
=
I
C
L
I
C
L
I
AS symmetry
S
L
I
CS
L
I
A S
L
I
AC
I C
L
I
S
L
I
I
C
L
I
C
L
I
AS CS
L
I
A C
L
I
C
L
I
AS
S
L
I
CS
L
I
A S
L
I
AC S
L
I
CS
L
I
A S
L
I
AC
L
E
K
4
6 12
6 12 12
2
6 6
4
6 12 12 6 12
6 12 12 6 12 12
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2
2


4.4 2-D Beam Element with combined loading Bending, Axial, and
Torsion ( = 0)

The torsional loads are m
1x
and m
2x
, and the corresponding deflections are,
x 1
and
x 2

The torsional structural equation is:



m
1x
1 -1
1x

= JG/L
m
1x
-1 1
2x

ME 273 Lecture Notes by R. B. Agarwal 4-14
DRAFT 1
Chapter 4 Beam Element



These terms can be superimposed on the stiffness equation derived previously for the
combined bending and axial loads.

d
y




y

3-D Beam Element: d
x


z

x

d
z



A 3-D beam element has 6 DOF at each node, and 12 DOF for each element. The
stiffness matrix can be derived by super-imposing the axial, bending, and torsion loadings
in the XY, XZ, and YZ planes. The equation is,























ME 273 Lecture Notes by R. B. Agarwal 4-15
DRAFT 1
Chapter 4 Beam Element


^
2
^
2
^
2
^
2
^
2
^
2
^
1
^
1
^
1
^
1
^
1
^
1 z y x z y x z y x z y x
d d d d d d
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(
(

=
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
GJ
L
GJ
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
AE
L
AE
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
GJ
L
GJ
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
EI
L
AE
L
AE
K
z z z z
y y y y
y y y y
z z z z
z z z z
y y y y
y y y y
z z z z
4
0 0 0
6
0
2
0 0 0
6
0
0
4
0
6
0 0 0
2
0
6
0 0
0 0 0 0 0 0 0 0 0 0
0
6
0
12
0 0 0
6
0
12
0 0
6
0 0 0
12
0
6
0 0 0
12
0
0 0 0 0 0 0 0 0 0 0
2
0 0 0
6
0
4
0 0 0
6
0
0
2
0
6
0 0 0
4
0
6
0 0
0 0 0 0 0 0 0 0 0 0
0
6
0
12
0 0 0
6
0
12
0 0
6
0 0 0
12
0
6
0 0 0
12
0
0 0 0 0 0 0 0 0 0 0
2 2
2 2
2 3 2 3
2 3 2 3
2 2
2 2
2 3 2 3
2 3 2 3


ME 273 Lecture Notes by R. B. Agarwal 4-16