Dynamic Analysis of Structures PDF
Dynamic Analysis of Structures PDF
@CivilMethod
Dynamic Analysis of
Structures
John T. Katsikadelis
Professor Emeritus of Structural Analysis
School of Civil Engineering
National Technical University of Athens
@CivilMethod
Academic Press is an imprint of Elsevier
125 London Wall, London EC2Y 5AS, United Kingdom
525 B Street, Suite 1650, San Diego, CA 92101, United States
50 Hampshire Street, 5th Floor, Cambridge, MA 02139, United States
The Boulevard, Langford Lane, Kidlington, Oxford OX5 1GB, United Kingdom
© 2020 Elsevier Inc. All rights reserved.
No part of this publication may be reproduced or transmitted in any form or by any means,
electronic or mechanical, including photocopying, recording, or any information storage and
retrieval system, without permission in writing from the publisher. Details on how to seek
permission, further information about the Publisher’s permissions policies and our
arrangements with organizations such as the Copyright Clearance Center and the Copyright
Licensing Agency, can be found at our website: [Link]/permissions.
This book and the individual contributions contained in it are protected under copyright by the
Publisher (other than as may be noted herein).
Notices
Knowledge and best practice in this field are constantly changing. As new research and experience
broaden our understanding, changes in research methods, professional practices, or medical
treatment may become necessary.
Practitioners and researchers must always rely on their own experience and knowledge in evaluating
and using any information, methods, compounds, or experiments described herein. In using
such information or methods they should be mindful of their own safety and the safety of others,
including parties for whom they have a professional responsibility.
To the fullest extent of the law, neither the Publisher nor the authors, contributors, or editors,
assume any liability for any injury and/or damage to persons or property as a matter of products
liability, negligence or otherwise, or from any use or operation of any methods, products,
instructions, or ideas contained in the material herein.
ISBN 978-0-12-818643-5
@CivilMethod
Dedication
To my wife Efi
for her loving patience and support,
to my children Stefan and Christina,
and to
my granddaughter Katharina
for patiently enduring and sharing the years of
preparation of the book with me.
@CivilMethod
Preface
The statement of the laws of motion by Newton 334 years ago (1686)a was a
milestone in the evolution of mechanics and modern engineering. The relation
between force (cause) and motion (effect) was quantified as a relation between
the linear momentum of the body and the force exerted on it. Thus, this relation
from a subject of philosophy up to that time turned out to be a valuable tool of
science for the study of the natural world. The subsequent developments in the
sciences were rapid. Astronomy, mathematics, mechanics of fluid and deform-
able bodies, and in general, mechanics of continuous media reached their peaks
in the centuries that followed, with immense applications to all engineering dis-
ciplines. Nevertheless, the laws of motion, which were stated as an axiom
(Axiomata sive Leges Motus) by Newton because, apparently, he could not jus-
tify their derivation, was a consequence of the discoveries of great scientists
who preceded him such as Galileo,b Kepler,c Hook, etc.
The implementation of the laws of motion leads to mathematical models
described by differential equations, ordinary or partial, whose solution effort
has given a great impetus to the development of mathematics. Unfortunately,
analytical solutions are limited to simple problems such as vibrations of discrete
systems with a few degrees of freedom; linear vibrations of beams, membranes,
plates, and shells with simple geometry; and simple support conditions made
from materials, mostly with a linear behavior. These solutions, while useful
for extracting qualitative conclusions about the dynamic response of structures,
are not capable of solving realistic problems in engineering, where the geometry
and loads are complicated while the response is generally nonlinear. Although it
has been shown that Newton’s law of motion is of an integer-order derivative,d
in recent years, the fractional derivatives have been proven more suitable for
modeling the actual structures. However, the use of fractional calculus has
not been employed in mathematical physics for three centuries because the
a. I. Newton, Philosophiae Naturalis Principia Mathematica, Royal Society Press, London, 1686.
b. J.T. Katsikadelis, Derivation of Newton’s law of motion using Galileo’s experimental data, Acta
Mech. 226 (2015) 3195–3204, doi:10.1007/s00707-015-1354.
c. J.T. Katsikadelis, Derivation of Newton’s law of motion from Kepler’s laws of planetary motion,
Arch. Appl. Mech. 88 (2018) (2017) 27–38, doi:10.1007/s00419-017-1245-x.
d. J.T. Katsikadelis, Is Newton’s law of motion really of integer differential form? Arch. Appl.
Mech. 89 (2019) 639–647, doi:10.1007/s00419-018-1486-3.
xv
@CivilMethod
xvi Preface
@CivilMethod
Preface xvii
particularly after the advent of computers. Therefore, many books have appeared
on this subject. However, these books provide a means only to understand the
response of simple and mostly unrealistic structures when subjected to dynamic
loads, especially to ground motion. The principles of dynamics are illustrated by
applying them to very simple models, which cannot describe actual structures
and therefore cannot be employed for dynamic analysis and design. We should
have in mind that the dynamic analysis of actual structures requires their model-
ing, the formulation of the governing equations of motion, their solution under
any dynamic load, and the physical interpretation of the results.
In the last 30 years, almost all seismic codes have encountered earthquake
ground motion as an effective dynamic load. The advent of computers in the
early 1960s encouraged engineers to develop methods of dynamic analysis
of structures, modeled first by the FEM and later by other advanced numerical
methods. Today, these methods constitute a powerful tool for dynamic engi-
neering analysis. Thanks to the availability of cheap computer power, every
engineer can use them. The essential ingredients of a book on Dynamics of
Structures for Civil Engineers should be:
(a) The basic concepts and principles of structural dynamics as they are
applied to particles as well as rigid and deformable bodies, enabling the
student or the engineer to formulate the equations of motion of any struc-
ture, no matter how complex, once the dynamic model has been adopted.
(b) Realistic modeling of actual structures under dynamic loads.
(c) Analysis of the dynamic response of the structure represented by its
model under any specified load. The analysis should include single- and
multiple-degree-of-freedom systems for linear and nonlinear response
under any dynamic excitation.
(d) Approximation of real structures using computational methods such as the
FEM, which replaces the actual structure (distributed parameter system)
with an approximate discrete system for which analysis methods can be
applied.
(e) Effective present-day numerical methods for dynamic analysis, including the
numerical solution of eigenvalue problems and the direct solution of the equa-
tions of motion of large systems, namely, systems with a large number of
degrees of freedom such as those resulting from the employed discretization.
Students attending a course on Dynamics of Structures should be exposed at
least to the above subjects. However, not all of them can be found in a single
book. Therefore, people interested in structural dynamics should refer to more
than one book in order to retrieve the required knowledge. Apparently, these
books cannot be used as integrated textbooks in the sense described above.
The student should be acquainted with different symbols and approaches, which
complicate the acquisition of knowledge, an approach that is, at least, educa-
tionally inappropriate.
@CivilMethod
xviii Preface
@CivilMethod
Preface xix
discussed. The derivation of the equivalent element nodal quantities, that is,
mass and stiffness matrices and forces, are derived by using the Lagrange equa-
tions instead of the principle of virtual works. Although the principle of vir-
tual works offers a handy tool for the derivation of these quantities, the use of
the Lagrange equations was preferred here. The reason is that the Lagrange
equations not only offer a straightforward method for the derivation of the
equivalent nodal quantities for all types of elements, especially in complex
systems with a nonlinear response, but they also allow understanding of their
physical significance. Chapter 12 studies the free vibrations of MDOF sys-
tems without and with damping. The linear eigenvalue problem is presented
from the mathematical point of view, aiming at drawing useful conclusions
about the eigenfrequencies and the mode shapes of the physical systems.
Chapter 13 presents the numerical methods for the computation of the eigen-
frequencies and mode shapes, especially for systems with a large number of
degrees of freedom such as those derived from the application of finite ele-
ments. Chapter 14 studies the forced vibrations of MDOF systems. A large
part of this chapter is devoted to the mode superposition method. It also dis-
cusses the use of Ritz vectors to reduce the degrees of freedom. Particular
emphasis is given to the response spectrum method. The response of linear
systems when they are subjected to a synchronous and an asynchronous
motion of the supports are also discussed. This chapter ends with the presen-
tation of the numerical methods, giving the respective pseudocodes for the
time integration of linear and nonlinear systems of equations of motion.
Chapter 15 discusses the approximation of multistory buildings by skeletal
structures and presents methods of formulating their equation of motion.
Finally, Chapter 16 discusses the response of seismically isolated buildings.
This chapter is introductory to the subject and aims primarily at understand-
ing the impact of base isolation on structures.
The book is supplemented by an appendix. Therein, the basic theory of rigid
body dynamics is presented for large and small displacements and the relevant
equations are derived, which are employed in the development of the material of
the book.
In closing, the author wishes to express his sincere thanks to his former stu-
dent and coworker Dr. A. J. Yiotis for carefully reading the manuscript as well
as for his suggestions, constructive recommendations and his overall contribu-
tion to minimizing the oversights in the text. Warm thanks also belong to Dr.
Nikos G. Babouskos, also a former student and coworker of the author, not only
for his careful reading of the manuscript and his apposite suggestions for
improvement of the book but also for his assistance in checking the computer
programs and in producing the numerical results of examples therein. Finally,
thanks belong to Dr. G. Dasios, professor of mathematics at the University of
Patras as well as to his former students, Dr. G. Tsiatas, associate professor of
mathematics at the University of Patras, and Dr. P. Tsopelas, associate professor
of mechanics at NTUA, for reading certain sections of the book and making
constructive suggestions.
@CivilMethod
xx Preface
J.T. Katsikadelis
Athens
April 2020
@CivilMethod
Chapter 1
1.1 Introduction
Apart from static loads, engineering structures may be subjected to dynamic
loads, that is, loads whose magnitude as well as direction of action and/or posi-
tion vary with time. The analysis of stresses and deflections developed in a
given structure undergoing dynamic loads is the fundamental objective of the
dynamic analysis of structures. Between static and dynamic analysis of struc-
tures, there exist two substantial differences:
(a) In static analysis, the loads are assumed time-invariant, and the resulting
response is unique, at least in linear theory. On the other hand, in dynamic
analysis the loads are time-varying and the deformations and stresses
depend on time, that is, at each instant the response of the structure is
different.
(b) In dynamics analysis, the material points of the structure change position
with the time, hence they have velocity and acceleration. Inasmuch as
the structure has a mass, inertial forces are produced due to the accelerations
of the material points. These inertial forces constitute an additional loading
that cannot be ignored. To make it tangible, we consider the cantilever beam
of Fig. 1.1.1a. The beam has a mass per unit length m and a flexural rigidity
EI , both assumed constant along the length, and it is subjected to the time-
(a)
(b)
(c)
FIG. 1.1.1 Vibrating cantilever beam.
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 5
The transverse deflection is a function not only of the spatial variable x but
also of time t, namely it is u ¼ u ðx, t Þ. As the element has a mass mdx, an inertial
€ arises, which, according to d’Alembert’s principle (see Section 1.5)
force m udx
opposes the motion, that is, if the positive transverse displacement u ðx, t Þ in the
beam is directed downward, the inertial force is directed upward (see Fig. 1.1.1b
and c). Similarly, due to angular acceleration ∂u€ðx, t Þ=∂x of the cross-section, an
inertial moment is also developed, which we may neglect [1]. Thus, referring to
Fig. 1.1.1b, we obtain the equation of dynamic equilibrium of the beam element
in the y direction as
€ ¼0
Q + Q + dQ + pðx, t Þdx m udx
or
∂Q ∂2 u
¼ pðx, t Þ + m 2 (1.1.1)
∂x ∂t
∂3 u
Q ¼ EI (1.1.2)
∂x 3
∂4 u ∂2 u
EI + m 2 ¼ pðx, t Þ (1.1.3)
∂x 4 ∂t
Eq. (1.1.3) is known as the equation of the dynamic equilibrium or the equa-
tion of motion of the vibrating beam. It is apparent that if we omit the inertial
term m∂2 u=∂t 2 in Eq. (1.1.3), we obtain the equation of the deflection of the
beam under static loading, that is,
d4u
EI ¼ pðx Þ (1.1.4)
dx 4
Fig. 1.1.1c shows the beam subjected to the inertial forces. These forces
resist the accelerations and they need to be accounted for in the solution. This
is the most important characteristic of the dynamic problem. Obviously, the
magnitude of the inertial forces depends on the magnitude of the acceleration.
When the produced accelerations are very small, as in the case of slow motion,
the inertial forces are very small too, and they can be neglected. In this case,
the time appears in the equation as a parameter and the response at any instant
@CivilMethod
6 PART I Single-degree-of-freedom systems
can be obtained by the static structural analysis, even though the load and the
response are time-varying. This response is pseudodynamic and is referred to
as quasistatic. The inertial forces appear in the equation of motion of the struc-
ture with the second derivatives of the displacements with respect to time.
Therefore, the equations that must be solved in dynamic analysis in order
to establish the deformations and stresses in the structure are differential
equations, contrary to static analysis where the governing equations are alge-
braic. For this reason, the solution procedure in dynamic analysis is essentially
different from that used in static analysis.
Dynamic loads can be classified into two great groups that characterize the
approach of evaluating the structural response: The deterministic dynamic loads
and the nondeterministic or random dynamic loads. In the first group are the
dynamic loads whose time variation is fully determined, regardless of the
complexity of their mathematical presentation. They are also referred to as pre-
scribed dynamic loads. They can be represented by an analytic or a generalized
function (Dirac or Heaviside) as well as numerically by a set of their values at
discrete time instances. The second group includes the loads, whose time var-
iation is not completely known but it can be defined in a statistical sense. In this
book, the dynamic response of structures only under deterministic loads is
studied.
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 7
5
2
H = me sin t
2.5
0
H(t)
–2.5
T T T
–5
0 2 4 6 8 10
t
FIG. 1.2.2 Harmonic loading due to an unbalanced rotating mass.
@CivilMethod
8 PART I Single-degree-of-freedom systems
structure and, as we will see later in this book (Chapter 6), it can be reduced to an
effective dynamic load if the accelerogram of the ground motion is known
(Fig. 1.2.4).
400
–400
0 10 20 30
FIG. 1.2.4 Effective dynamic load pðt Þ ¼ m u€g ðt Þ due to seismic ground motion.
If we examine static loading closer, we will see that even what we call static
loads are actually dynamic in nature. They are applied starting from a zero value
until the final prescribed value is reached within a time span. That is, they are
time-varying, thus dynamic. However, the duration of the application of the
static load is longer than the period of vibration of the structure. This produces
negligible accelerations and consequently the response under a “static load”
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 9
@CivilMethod
10 PART I Single-degree-of-freedom systems
Massless
Massless
columns
columns
(a) (b)
(c)
FIG. 1.3.1 Systems with one degree of freedom (SDOF).
Rigid
Rig
(a) (b)
FIG. 1.3.2 Systems with two degrees of freedom (2 DOF).
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 11
The lumped mass idealization provides a simple means of reducing the num-
ber of degrees of freedom. Fig. 1.3.3 represents the discrete model of a canti-
lever column, whose mass has been localized at three points. Neglecting the
axial deformation of the column and considering plane motion, the system
has six degrees of freedom, the three transnational ui ðt Þ and the three rotational,
fi ðt Þ. If the masses are fully concentrated so that their rotational inertia can be
ignored, the inertial moments Ii f€i are zero and the number of dynamic degrees
of freedom reduces to three. Obviously, the number of degrees of freedom
increases with the number of nodal points, where the mass of the structure is
lumped. As the number of points becomes infinitely large, the discretized struc-
ture approaches the continuous system (Fig. 1.3.4).
@CivilMethod
12 PART I Single-degree-of-freedom systems
Spring
Damper
Frictionless rollers
FIG. 1.4.1 Model of a SDOF system.
Center of mass
The forces applied to the body at time t are shown in the free body diagram
of Fig. 1.4.2. These are
(a) The external load pðt Þ
(b) The elastic force fS
(c) The damping force fD
(d) The inertial force fI .
The spring force fS depends on the displacement u ðt Þ and it is generally
expressed by a nonlinear function, fS ¼ fS ðu Þ. For linear response of the struc-
ture, the force fS is proportional to the displacement and is given by
fS ¼ ku (1.4.1)
where k is the constant that represents the spring stiffness coefficient, that is, the
force required to change the length of the spring by a unit. The force fS repre-
sents the elastic force of the structure that resists the motion and tends to bring
the body to its initial undeformed position.
The damping force fD also resists the motion. It represents the energy loss
due to internal or external dissipative forces. Damping forces are complex in
nature. Their exact expression in terms of the parameters of motion and of
the geometrical and material properties of the structure is complicated and dif-
ficult to determine. The simplest form of damping is linear viscous damping.
This produces damping forces, which are the easiest to handle mathematically
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 13
and provide analytical results for the response of a system close to the experi-
mental ones. The linear viscous damping mechanism is indicated by a dashpot,
as shown in Fig. 1.4.1. In viscous damping, the resisting force is proportional to
the velocity
fD ¼ cu_ (1.4.2)
where c is a constant that can be established experimentally. Inasmuch as the
work done by this force is converted to heat, the damping force is a nonconser-
vative force. It is the force that makes the amplitude of a vibrating
structure decay.
The inertial force fI depends on the mass m of the body and its acceleration
€ It also resists the motion. It is given by Newton’s second law of motionb
u.
fI ¼ m u€ (1.4.3)
A simple example of a structure that can be modeled as SDOF is the
one-story, one-bay frame of Fig. 1.4.3. It consists of two identical weightless
columns fixed on the ground and having height h, cross-sectional moment of
inertia Ι, and modulus of elasticity E. The cross-sectional moment of inertia
of the horizontal beam is assumed infinitely large. This means that the beam
behaves like a rigid body of mass m and hence the cross sections of the columns
at the roof level cannot rotate when the frame deforms. The frame is subjected to
an external horizontal force pðt Þ, as shown in Fig 1.4.3a, which forces the frame
to move. Neglecting the axial deformation of the beam and columns, an allow-
able assumption for frames, the only possible movement is the displacement
u ðt Þ at the roof level. The rotation of the beam as a rigid body is excluded
because this would cause a change in the length of columns.
(a)
(b) (c)
FIG. 1.4.3 Two-column shear frame.
b. Actually, this form of Newton’s law of motion is attributed to L. Euler, who defined it indepen-
dently as a mechanical principle [2, 3]. This law was recently derived from Kepler’s laws of plan-
etary motion [4].
@CivilMethod
14 PART I Single-degree-of-freedom systems
Referring to Fig. 1.4.3b, we see that the elastic forces are the shear forces Q
at the top cross-sections of the columns. These forces are given by the known
relation of statics
12EI
Q¼ u ðt Þ (1.4.4)
h3
The quantity 12EI =h 3 represents the translational stiffness of the column.
This is the force required to produce a unit relative displacement between
the end cross-sections of the column. These shear forces tend to restore the
frame to the undeformed position. Therefore, they play the role of the spring
in the SDOF model with a stiffness coefficient
12EI
k ¼2 (1.4.5)
h3
The inertial force is given by fI ¼ m u€ while the damping force by fD ¼ cu. _
Another convenient model to represent the single-story frame is shown in
Fig. 1.4.3c. It consists of a mass m placed at the top of a column with transla-
tional stiffness equal to the sum of the translational stiffness coefficients of the
columns of the frame. During the motion, the top cross-sections of columns
undergo only the translational displacement u ðt Þ. Models of this type are also
suitable to idealize multistory shear frames (see Fig. 1.4.4), in which the masses
are placed at the floor levels and the girders are assumed rigid.
FIG. 1.4.4 Four-story shear frame and its model without damping.
(a) (b)
FIG. 1.4.5 Two-story, two-bay shear frame and its model without damping.
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 15
Fig. 1.4.5a shows another two-story shear frame. The columns are assumed
weightless. Fig. 1.4.5b shows its dynamic model. The column 1-2 is represented
by a spring of stiffness k ¼ 12EI =h 3 . The stiffness coefficients k1 and k2 include
only the stiffness of the columns with heights h1 and h2 , respectively.
Given the dynamic model of the structure, the equation of motion of the sys-
tem is formulated. For the SDOF system, the equation of motion can be formu-
lated using Newton’s second law of motion as it is applied for the motion of a
particle
m u€ ¼ F (1.4.6)
where
F ¼ pðt Þ fS fD (1.4.7)
is the resultant of the external forces. Using Eqs. (1.4.1), (1.4.2), (1.4.7),
Eq. (1.4.6) is written
m u€ + cu_ + ku ¼ pðt Þ (1.4.8)
Eq. (1.4.8) is the equation of motion of the SDOF system. The equation of
motion represents the dynamic equilibrium of the system. It is an ordinary dif-
ferential equation of the second order with respect to the unknown variable u ðt Þ.
The solution of this equation yields the displacement as a function of time. For
MDOF systems, the number of equations of motion that must be formulated is
equal to the number of dynamic degrees of freedom. The use of Newton’s law of
motion is not always well suited to formulate the equations, especially for
MDOF systems or complex SDOF systems. It requires advanced knowledge
of the dynamics of the rigid and deformable body as well as mastering various
special methods. Generally, the equations of motion can be formulated using:
(a) d’Alembert’s principle or method of equilibrium of forces.
(b) Principle of virtual work.
(c) Hamilton’s principle.
(d) Lagrange’s equations.
These methods will be presented in the following and will be demonstrated by
appropriate examples. The acquaintance with the application of these methods
constitutes a fundamental presupposition for the analysis of the dynamic
response of structures.
@CivilMethod
16 PART I Single-degree-of-freedom systems
where F is the resultant of all external forces acting on the particle of mass m and
€ is its acceleration with respect to an inertial frame of reference.c If we consider
u
that the term m€ u is another force, known as inertial force, then Eq. (1.5.1) states
that the vector sum of all forces, external and inertial, is zero during the motion.
But this is the necessary and sufficient condition for the static equilibrium of the
particle. Thus, in a sense, the dynamic problem is reduced to a problem of statics
according to the following statement, known as d’Alembert’s principle.
The laws of static equilibrium can be applied also to a dynamic system with
respect to an inertial frame of reference if the inertial forces are considered as
applied forces on the system together with the actual external forces.
The motion of a rigid body of mass m with respect to an inertial frame of
reference X, Y ,Z is decomposed into a translational motion of its center of
mass, where the whole mass is considered to be concentrated, and a rotational
motion about it (Fig. 1.5.1).
FIG. 1.5.1 Rigid body moving with respect to the inertial. frame X,Y , Z .
c. In classical dynamics, an inertial frame of reference is a frame of reference in which a body with
zero force acting upon it is not accelerating; that is, the body is at rest or it is moving at a constant
velocity in a straight line [5].
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 17
respect to the center to mass, and H_ c is the rate of change of the angular momen-
tum Hc of the body with respect to the same point given as
ZZZ
H_ c¼ r r€rdV (1.5.3)
V
Path of P
FIG. 1.5.2 Plane body moving in the XY plane. The system of xy axes moves with P without
rotating.
@CivilMethod
18 PART I Single-degree-of-freedom systems
1 2 1
T ¼ m X_ c + Y_ c + Ic w2
2
(1.5.7)
2 2
(b) with respect to an arbitrary point P of the body (K€onig’s theorem)
1 2 1
T ¼ m X_ p + Y_ p + Ip w2 + m xc Y_ p yc X_ p w
2
(1.5.8)
2 2
We shall write now Eqs. (1.5.4a)–(1.5.4c) in terms of the displacement
vector. Apparently, the displacement vector from the beginning of the motion
is defined as
u ¼ Rðt Þ Rð0Þ ¼ u ðt Þi + v ðt Þj (1.5.9)
where
u ¼ X ðt Þ X ð0Þ, v ¼ Y ðt Þ Y ð0Þ (1.5.10)
Hence, X€ ¼ u,
€ Y€ ¼ v€. Moreover, if fðt Þ represents the change of the rota-
_ w_ ¼ f,
tion in the same time interval and set w ¼ f, € Eqs. (1.5.4a)–(1.5.4c) are
written in terms of displacements as
Fx ¼ m u€c (1.5.11a)
Fy ¼ m v€c (1.5.11b)
Mc ¼ Ic f€ (1.5.11c)
or in matrix form
€c
Fc ¼ m c U (1.5.12)
where
8 9 8 9 2 3
< Fx = < uc = m 0 0
Fc ¼ Fx , Uc ¼ vc , mc ¼ 4 0 m 0 5 (1.5.13)
: ; : ;
Mc f 0 0 Ic
are the force vector, the displacement vector, and the mass matrix of the body,
respectively.
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 19
€p
Fp ¼ m p U (1.5.15)
where
8 9
>
< Fx >
=
Fp ¼ F x (1.5.16a)
>
: >
;
Mp
8 9
< up >
> =
Up ¼ vp (1.5.16b)
: >
> ;
f
2 3
m 0 my c
mp ¼ 4 0 m mx c 5 (1.5.16c)
my c mx c Ic
Note that the mass matrix is not diagonal when the point of reference is not
the center of mass.
Finally, Eqs. (1.5.7), (1.5.8) are written as
1 1
T ¼ m u_ 2c + v_ 2c + Ic f_
2
2 2
(1.5.17)
1_T _c
¼ U mc U
2 c
1 1
T ¼ m u_ 2p + v_ 2p + Ip f_ + m xc u_ p yc v_ p f_
2
2 2
(1.5.18)
1_T _p
¼ U mp U
2 p
The set of equations with reference to point P can also be derived from the
set of equations with reference to point C by transforming the displacements
and the forces from point C to P (see Section 10.7).
Example 1.5.1 Equation of motion of an elastically supported body
Consider the rigid plate of constant thickness and total mass m shown in
Fig. E1.1a. The plate is hinged at O and elastically supported at A. Formulate
@CivilMethod
20 PART I Single-degree-of-freedom systems
the equation of motion of the system for small amplitude motion using the
method of equilibrium of forces.
Solution
The only possible motion of the plate is the rotation in its plane about the point
O. Hence, the system has one degree of freedom. The motion can be described
either by the rotation fðt Þ about O or by the translational displacement of
a point, for example, the displacement u ðt Þ of point B, which is related
to fðt Þ as
d2 1 b
fIx ¼ m ðCC 0 Þx ¼ m u€ (2c)
dt 2 2 a
d2 1
fIy ¼ m 2
ðCC 0 Þy ¼ m u€ (2d)
dt 2
u€
MIc ¼ IC f€ ¼ IC (2e)
a
The quantities ðCC 0 Þx and ðCC 0 Þy are the horizontal and the vertical dis-
placements of the center of mass C due to rotation, respectively. They are
obtained from Fig. E1.1b as
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 21
(a)
(b)
(c)
FIG. E1.1 Rigid plate in Example 1.5.1.
1b
ðCC 0 Þx ¼ ðOC Þf sin b ¼ u (3a)
2a
1
ðCC 0 Þy ¼ ðOC Þf cos b ¼ u (3b)
2
The equation of motion results from the dynamic equilibrium of moments
with respect to point O. Thus, we obtain
a 2a b a
W fS fIx fIy MIc + pðt Þa ¼ 0 (4)
2 3 2 2
@CivilMethod
22 PART I Single-degree-of-freedom systems
(a)
(b)
FIG. E1.2 Frame with a rigid column in Example 1.5.2.
Solution
The only possible motion of the system is the rotation fðt Þ of the column as a
rigid body about the hinged support at point A of its base. Because the rotation is
small, we have:
sin f f, cos f 1, f2 0
Hence
h h
u ¼ L sin f ¼ Lf, d ¼ sin f f
2 2
@CivilMethod
24 PART I Single-degree-of-freedom systems
h h
ð1 2Þ ¼ sin f + L cos f f + L,
2 2
h h
ð3 4Þ ¼ L sin f + cos f Lf + ,
2 2
h h
ð5 6Þ ¼ cos f
2 2
The forces applied on the column are shown in Fig. E1.2b. These are:
The elastic moment at the corner C MS ¼ ð1:5L
6EI
Þ2
4EI
d + 1:5L f
The elastic moment due to the rotational spring MR ¼ CR f ¼ EI
L f
MIA ¼ IA f€ ¼ mL €
2
The moment of inertia of the mass m 3 f
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 25
(a)
(b)
FIG. E1.3 System with two rigid bodies in Example 1.5.3.
Solution
As the bar AF is rigid, the only possible motion is its rotation about A. Hence,
the system has a SDOF. Its motion can be described either by the angle of rota-
tion fðt Þ about the hinge at A or by the transverse displacement of any point
along the axis of the bar. We choose the upward displacement u ðt Þ at point
B as the parameter of the motion. For small amplitude motion, the forces acting
on the system are shown in Fig. E1.3b. These are:
The elastic force fS at B: As it opposes the motion, it is directed downward
and is expressed as
fS ¼ ku (1)
The damping force fD at G: It is directed also downward and is expressed as
d d
fD ¼ c ðGG 0 Þ ¼ c ð1:625u Þ ¼ 1:625cu_ (2)
dt dt
The inertial force fIK and the inertial moment MIK at the center of mass K of
the bar due the distributed mass m are
d2
Þ
fIK ¼ ðm3L ðKK 0 Þ ¼ 0:75m u€ (3)
dt 2
MIK ¼ IK f€
@CivilMethod
26 PART I Single-degree-of-freedom systems
m ∗ u€ + c ∗ u_ + k ∗ u ¼ p ∗ ðt Þ (8)
where
m ∗ ¼ 6:833m, c ∗ ¼ 5:281c, k ∗ ¼ 2k, p ∗ ðt Þ ¼ 2L
pðt Þ (9)
The quantities defined by Eq. (9) are referred to as the generalized mass, the
generalized damping, the generalized stiffness, and the generalized load,
respectively.
Once the dynamic displacement u ðt Þ is established from the solution of
Eq. (8), the vertical reaction RA can be evaluated from the dynamic equilibrium
of forces in the direction of the y axis. This yields
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 27
(b)
(a)
FIG. E1.4 Single-story shear building in Example 1.5.4.
Solution
Taking into account that the structure is symmetric with respect to the x axis,
the columns are inextensible, and the load pðt Þ acts on the axis of symmetry, the
only possible motion of the plate is the horizontal displacement u ðt Þ in the
direction of the x axis. The SDOF model of the structure is shown in Fig. E1.4b.
The total mass of the system is due to the load of the plate and to half the
weight of the columns
5 10 20 + ð4 0:3 0:3 + 2 0:3 0:2Þ 2 24
m¼ ¼ 104:285
9:81
The stiffness of the system is equal to the sum of the translational stiffness
coefficients of all columns, which are given as
12EI i
ki ¼
hi3
where Ii is the moment of inertia of the cross-section of the i column with
respect to the y axis through its mass center and hi its height. Thus, we have
@CivilMethod
28 PART I Single-degree-of-freedom systems
Columns 30 30:
0:304
12 2:1 107
k3030 ¼ 12 ¼ 2657:8kN=m
43
Columns 30 20:
0:303 0:20
12 2:1 107
k3020 ¼ 12 ¼ 1771:9kN=m
43
Therefore the stiffness of the system is
k ¼ 4 2657:8 + 2 1771:9 ¼ 14175:0kN=m
The equation of motion results from the equilibrium of the forces shown in
Fig. E1.4b. This yields
fI fS + pðt Þ ¼ 0
or
m u€ + ku ¼ pðt Þ
Substituting the numerical values for m, k and the expression for pðt Þ, the
above equation of motion becomes
Rigid
Rigid
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 29
m€
u + ku ¼ pðt Þ (3)
where
u1 m1 0 k1 k1 p1 ðt Þ
u¼ , m¼ , k¼ , pðt Þ ¼
u2 0 m2 k1 k1 + k2 p2 ðt Þ
Solution
We choose O xy as the system of reference of the motion, whose origin coin-
cides with point O at the beginning of motion. Let xi , yi represent the coordi-
nates of the center of mass of the cross-section of i column and fi the angle
between its principal x axis and the x axis. The axes xy will be referred to
as the global axes of the system while the axes xy as the local axes of the
column.
Inasmuch as the axial deformation of columns is ignored, the only possible
motion of the plate is inside its plane, which can be determined by the two
translational displacements of a point and the rotation of the plate. We study
the motion of the plate with reference to point O and let U , V represent its
translational components with respect to the global axes xy, which are
assumed fixed in the plane, and W the rotation of the plate. As a consequence
of this motion, the cross-section of the i column at the level of the plate
undergoes the displacements u i , v i , wi , with respect to its base. These displace-
ments generate elastic forces X i , Y i , M i , which act on the plate. Thus,
we define the following vectors and matrices that will be used in the subse-
quent analysis.
@CivilMethod
30 PART I Single-degree-of-freedom systems
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 31
The transformation matrix for the vector quantities related to i column from the
global axes to the local axes is given as
2 3
cos fi sin fi 0
Ri ¼ 4 sin fi cos fi 0 5
0 0 1
Hence, the vectors are transformed from one system of axes to the other as
i
D i ¼ Ri D (1a)
i ¼ Ri T D i
D (1b)
i
FiS ¼ Ri F (2a)
S
i ¼ Ri T F i
F (2b)
S S
where
2 3
cos fi sin fi 0
R ¼ 4 sin fi
i T
cos fi 0 5
0 0 1
i 1 i T
is the transpose of R . Note that R
i
¼ R because Ri is orthonormal.
The elastic forces X , Y , M are related to the displacements u i , v i , wi by
i i i
12EI y i
Xi ¼ u (3a)
h3
12EI x i
Yi ¼ v (3b)
h3
GI t i
Mi ¼ w (3c)
h
where Ix ,Iy are the principal moments of inertia of the column cross-section and
It is the torsional constant, E and G are the material constants, and h is the
height of the column.
Setting
12EI y 12EI x GI t
i
k11 ¼ i
, k22 ¼ i
, k33 ¼ (4)
h3 h3 h
Eqs. (3a)–(3c) can be written in matrix form as
8 9 2 i 38 9
< Xi = k11 0 0 < u i =
Y i ¼ 4 0 k22 i
0 5 vi
: i; i : i;
M 0 0 k33 w
or
Fi ¼ k i D i (5)
@CivilMethod
32 PART I Single-degree-of-freedom systems
The matrix
2 i
3
k11 0 0
6 7
ki ¼ 4 0 i
k22 0 5
i
0 0 k22
or
F i D
i ¼ k i (6)
where
i ¼ Ri T k i R i
k (7)
is the stiffness matrix of the column in global axes, which becomes after per-
forming the matrix multiplications
2 i 3
k11 k12 0
i
6 7
i ¼ 6 ki ki 0 7
k (8)
4 21 22 5
0 0 k
i
33
where
9
k11 ¼ k11 >
i i
cos 2 fi + k22
i
sin 2 fi >
>
>
>
>
=
k22 ¼ k11
i i
sin 2 fi + k22
i
cos 2 fi
i (9)
k12 ¼ k21 ¼ k11 sin fi cos fi >
i i
k22
i >
>
>
>
>
;
i
k 33 ¼ k33
i
Inasmuch as the plate is rigid, the displacements ui , vi , wi of the i column
of point O. The geometrical rela-
depend on the plate displacements U , V , W
tions result from the following consideration.
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 33
(b) The rotation of the plate about O. Referring to Fig. E1.7 and observing that
¼ wi , we obtain
cos ai ¼ xi =ri , sin ai ¼ yi =ri , W
i
u r ¼ ri W sin ai ¼
yi W
i
v r ¼ ri W cosai ¼ xi W
Thus, we have
ui ¼ ui t + ui r ¼ U yi W (10a)
vi ¼ vi t + vi r ¼ V + xi W (10b)
wi ¼ W (10c)
The previous equations are written in matrix form as
8 9 2 38 9
< ui = 1 0 y i < U =
vi ¼ 4 0 1 xi 5 V (11)
: i; :;
w 0 0 1 W
or setting
2 3
1 0
yi
ei ¼ 4 0 1 xi 5
0 0 1
@CivilMethod
34 PART I Single-degree-of-freedom systems
X
K
P y ðt Þ Y ¼ m V€
i €
+ xc W (13b)
i¼1
K
X
ðt Þ
M xi Y yi X + M
i i €
i ¼ m xc V€ yc U€ + Io W (13c)
i¼1
or
X
K
ðt Þ
P
T i
ei F €
¼ ðec ÞT mðec ÞU (14)
i¼1
where
8 9
< P x ðt Þ >
> =
ðt Þ ¼ P y ðt Þ
P (15a)
: >
> ;
M ðt Þ
2 3
m 0 0
m ¼ 40 m 0 5 (15b)
0 0 Ic
Finally, using Eqs. (6), (12), we obtain the equation of motion
M € + K
U U ¼P
ðt Þ (16)
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 35
where
2 3
m 0 m yc
¼ ðec ÞT mec ¼ 4 0
M m m xc 5 (17a)
m yc m xc Io
X
K
¼
K ei
T T
Ri ki Ri ei (17b)
i¼1
T
d. The notation eT ¼ ðe1 Þ is employed.
@CivilMethod
36 PART I Single-degree-of-freedom systems
The off-diagonal terms cause coupling between the elastic force in one
direction and the displacement in another direction. For example, the element
k12 represents the force acting in the direction of the x axis when the plate
undergoes a unit displacement in the direction of the y axis. Similarly, the ele-
ment k31 represents the moment acting on the plate about the z axis, if the plate
undergoes a unit displacement, U ¼ 1, in the direction of the x axis. The elastic
center or center of resistance of the plate is defined as the point of the plate
where an applied force in any direction does not produce rotation. This implies
the vanishing of the elements k13 and k23 (hence also their symmetric k31 and
k32 ) in the stiffness matrix (23). This point can be established as follows.
The stiffness matrix K is transformed from point O to the sought elastic cen-
ter E ðxE , yE Þ according to Eq. (22a), if ec is replaced by eE . Namely
2 32 32 3
1 0 0 k11 k12 k13 1 0 yE
T 1
E ¼ eE
K eE
K ¼ 40 1 0 54 k21 k22 k23 54 0 1 xE 5
yE xE 1 k31 k32 k33 0 0 1
or after performing the matrix multiplications
2 E E E3
k11 k12 k13
6 7
6 E E E7
K ¼ 6 k21 k22 k23 7
E 6
7
4 5
E E E
k 31 k 32 k 33
2 3
k11 k12 k11 yE k12 xE + k13
6 7
6 7
¼ 6 k21 k22 k21 yE k22 xE + k23 7
4 5
k11 yE k21 xE + k31 k12 yE k22 xE + k32 2 2
k 13 yE + k 23 xE + k 33
The vanishing of the elements k13 and k23 yields
E E
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 37
Thus, the stiffness matrix with respect to the elastic center takes the form
2 E E 32 3
k11 k12 0 k11 k12 0
6 7
E ¼ 6 kE kE 0 76
K 7
4 21 22 54 k 21 k 22 0 5
0 0 k 13 yE + k 23 xE + k 33
0 0 k
E 2 2
33
2k12
tan 2q ¼ (26)
k 22 k11
The axes defined by angle q are referred to as the principal directions of stiff-
ness of the structure. The stiffness matrix becomes now diagonal
2 3
k^11 0 0
K^E ¼6 4 0 k^22 0 5
7
0 0 k^33
where
k^11 ¼ k11 cos 2 q + k22 sin 2 q k12 sin 2q (27a)
k^22 ¼ k11 sin 2 q + k22 cos 2 q + k12 cos 2q (27b)
k^33 ¼ k13 y2E + k23 x2E + k33 (27c)
@CivilMethod
38 PART I Single-degree-of-freedom systems
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 39
(a)
(b)
FIG. E1.8 System in Example 1.6.1.
Solution
Because the cable is inextensible, the displaced configuration of the system can
be specified either by the angle of rotation of one of the bars or by the transverse
displacement of a point on it. Thus, the system has only one degree of freedom.
If the upward transverse displacement u ðt Þ of point C is taken as the basic
parameter of the motion, then all other displacements can be expressed in terms
of it. Fig. E1.8b shows the deformed system with all forces applied to it.
The elastic forces fS1 and fS2 are due to the deformation of the springs k1 and
k2 . They are directed downward as they oppose the motion. The force fD is
due to the viscous damping mechanism and is directed upward as it also
opposes the motion. The inertia moments MIA , MIO , and MIE are due to the rota-
tion of the masses about A, O , and E, respectively. All forces are expressed in
terms of the single displacement u ðt Þ
fS1 ¼ k1 ðBB 0 Þ ¼ ku=2, fS2 ¼ k2 ðCC 0 Þ ¼ 2ku
d
fD ¼ c ðDD 0 Þ ¼ cu_
dt
ð2LÞ3 u€
m
MIA ¼ IA f€1 ¼ 2 u€
¼ 1:333mL
3 2L
ð1:5LÞ3 u€
m
MIE ¼ IE f€2 ¼ 2 u€
¼ 0:750mL
3 1:5L
ð0:8LÞ2 u€
MIO ¼ IO f€3 ¼ mL
2 u€
¼ 0:200mL
8 0:4L
@CivilMethod
40 PART I Single-degree-of-freedom systems
If point C is given a virtual displacement du, the forces ride the following
displacements
d ðCC 0 Þ ¼ du, dðBB 0 Þ ¼ du=2, dðDD0 Þ ¼ du
df1 ¼ du=2L, df2 ¼ du=1:5L, df3 ¼ du=0:4L
du ðx Þ ¼ xdf1 ¼ xdu=2L
The work done by the forces acting on the system due to the virtual displace-
ment should be set equal to zero, that is,
fS1 d ðBB 0 Þ fS2 dðCC 0 Þ fD dðDD 0 Þ MIA df1
Z L
(1)
MIE df2 MIO df3 + pðt Þdu ðx Þdx ¼ 0
0
Using the expressions for the forces and the displacements in terms of the
basic displacement derived previously, Eq. (1) yields
0:25ku 2ku cu_ 1:333mL
2 u=2L
€ 2 u=1:5L
0:750mL €
0:200mL €
2 u=0:4L + pðt ÞL=4du ¼ 0
or, inasmuch as du 6¼ 0, the expression within the square brackets should vanish.
This yields the equation of motion
m ∗ v€ + c ∗ v_ + k ∗ v ¼ p ∗ ðt Þ (2)
where
m ∗ ¼ 1:667mL, c ∗ ¼ c, k ∗ ¼ 2:25k, p ∗ ðt Þ ¼ 0:25
pðt ÞL
Example 1.6.2 Equation of motion of a rigid body assemblage
Formulate the equations of motion of the rigid body assemblage shown in-
Fig. E1.9a by using the principle of virtual displacements on the basis of small
amplitude motion.
Solution
Due to the spring k1 , the rigid bars can rotate independently from each other
about their hinged supports at A and F. Hence, the system has two degrees
of freedom. Its motion can be specified by the transverse downward displace-
ments u1 ðt Þ and u2 ðt Þ of points C and E, respectively. The forces applied to the
displaced system are shown in Fig. E1.9b. They are
The elastic force fS1 ¼ k1 ðCC 0 Þ ¼ k ðu2 u1 Þ
The elastic force fS2 ¼ k2 ðDD 0 Þ ¼ 4ku 2
The damping force fD ¼ c dtd ðBB 0 Þ ¼ c u_21
The inertial moment MIA ¼ IA f€1 ¼ IA 2a u€1 2
¼ 4ma
3 u €1
€
The inertial moment M ¼ IF f2 ¼ IF ¼
F u€2 2
8ma
u€2
I a 3
The system is given a virtual displacement pattern du1 and du2 corresponding
to the two degrees of freedom. The forces ride the following displacements
du1 du1
d ðBB 0 Þ ¼ , dðCC 0 Þ ¼ du1 , df1 ¼
2 2a
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 41
(a)
(b)
FIG. E1.9 System in Example 1.6.2.
du2
d ðDD 0 Þ ¼ 2du2 , d ðEE 0 Þ ¼ du2 , df2 ¼
a
According to the principle of virtual displacements, the work done by the
applied forces must be equal to zero, that is,
Inasmuch as the quantities du1 and du2 are arbitrary, Eq. (2) is valid
only if
2ma u_ 1
u€1 + c k ðu2 u1 Þ pðt Þ ¼ 0 (3a)
3 4
8ma
u€2 + k ð9u2 u1 Þ ¼ 0 (3b)
3
@CivilMethod
42 PART I Single-degree-of-freedom systems
Eqs. (3a), (3b) are the equations of motion of the system. In matrix form they
are written as
2 3 2 c 3( ) "
2ma #( ) ( )
6 3 0 7 u€1 0 u_ 1 k k u1 pðt Þ
4 44 5 ¼
5 u€2 +
8ma u_ 2
+
k 9k u2 0
(4)
0 0 0
3
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 43
d 2r
m F¼0 (1.7.1)
dt 2
We confine our attention to an interval of time during which the particle
moves from point 1 at t ¼ t1 to point 2 at t ¼ t2 . We consider now a varied path,
specified by rðt Þ + drðt Þ, adjacent to the actual one. We will refer to the quan-
tity drðt Þ ¼ dx ðt Þi + dy ðt Þj + dz ðt Þk as the variation of r. The only restriction is
that the two paths coincide at time t ¼ t1 and t ¼ t2 . This implies that the var-
iation dr ¼ drðt Þ vanishes at these instants, that is,
The first step to derive Hamilton’s principle is to take the inner product of
the left side of Eq. (1.7.1) with the vector dr and to integrate from time t1 to time
t2 . This gives
Z t2
d2r
m dr F dr dt ¼ 0 (1.7.3)
t1 dt 2
Integrating by parts the first term in the above integral and knowing that the
operator d acts like the differential operator [6], we obtain
Z
t2 Z t2
t2
d 2r dr dr dr
m 2 drdt ¼ m dr m d dt
t1 dt dt t1 t1 dt dt
@CivilMethod
44 PART I Single-degree-of-freedom systems
The term outside the integral is equal to zero because of Eq. (1.7.2). More-
over, we can write the integrand as
" #
dr dr 1 dr dr 1 dr 2 1 dr 2
m d ¼m d ¼m d ¼ d m ¼ dT
dt dt 2 dt dt 2 dt 2 dt
where
2
1 dr
T¼ m ¼ x_ ðt Þ2 + y_ ðt Þ2 + z_ ðt Þ2 (1.7.4)
2 dt
is the kinetic energy of the particle. Hence, the integral (1.7.3) takes the form
Z t2
ðdT + F drÞdt ¼ 0 (1.7.5)
t1
The variation dr is a virtual displacement that leads from the actual path to
the varied one. Hence the term F dr in Eq. (1.7.5) is the virtual work done
by the force Fðt Þ. Eq. (1.7.5) is a statement of Hamilton’s principle as it is
applied to a particle. This equation can be transformed into a more convenient
form if the force Fðt Þ is separated in its conservative and nonconservative
components, that is
Fðt Þ ¼ Fc ðt Þ + Fnc ðt Þ (1.7.6)
A potential function A ¼ Aðx, y, z, t Þ exists from which the conservative
force Fc ðt Þ is derived as its minus gradient
∂A ∂A ∂A
Fc ¼ i+ j+ k (1.7.7)
∂x ∂y ∂z
Hence
∂A ∂A ∂A
Fc dr ¼ dx + dy + dz
∂x ∂y ∂z
or
Fc dr ¼ dA (1.7.8)
Hence, Hamilton’s principle, Eq. (1.7.5), can be written as
Z t2 Z t2
d ðT AÞdt + dWnc dt ¼ 0 (1.7.9)
t1 t1
where
dWnc ¼ Fnc dr
represents the virtual work of the nonconservative force.
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 45
has a stationary value. In fact, it can be shown that this value is the minimum
value of the integral.
The derivation of Hamilton’s principle for a particle can be extended to
MDOF systems as well as to continuous systems. The potential energy usually
arises from the gravity field. However, it may also arise from other sources such
as electrical and magnetic fields. The strain energy U ðt Þ should be included as
an additional potential energy. Thus, we can write
Z t2 Z t2
d ðU T + AÞdt dWnc dt ¼ 0 (1.7.13)
t1 t1
@CivilMethod
46 PART I Single-degree-of-freedom systems
The next step is to remove the variation d u_ of the velocity u_ from Eq. (6).
This is achieved using integration by parts as follows:
Z t2 Z t2
du
_ udt
m ud _ ¼ _
m ud dt
t1 t1 dt
Z t2
d
¼ m u_ ðdu Þdt (7)
t1 dt
Z t2
_ tt21
¼ ½m udu €
m ududt
t1
In order that the integral in Eq. (9) is equal to zero for any time interval
½t1 , t2 , its integrand should vanish, that is,
½m u€ + cu_ + ku pðt Þdu ¼ 0
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 47
and taking into account that du1 ðt1 Þ ¼ du1 ðt2 Þ ¼ du2 ðt1 Þ ¼ du2 ðt2 Þ ¼ 0, we
obtain
Z t2 Z t2
dTdt ¼ ðm1 u€1 du1 + m2 u€2 du2 Þdt (3)
t1 t1
Moreover, it is
p
dWnc ¼ p1 ðt Þdu1 + p2 ðt Þdu2 and A ¼ 0 (4)
Introducing Eqs. (1), (3), (4) into Hamilton’s principle, Eq. (1.7.13), we
obtain Z t2
½k1 ðu1 u2 Þdu1 k1 ðu1 u2 Þdu2 + k2 u2 du2 + m1 u€1 du1
t1
+ m2 u€2 du2 p1 ðt Þdu1 p2 ðt Þdu2 dt ¼ 0
@CivilMethod
48 PART I Single-degree-of-freedom systems
or
Z t2
f½m1 u€1 + k1 ðu1 u2 Þ p1 ðt Þdu1 + ½m2 u€2 k1 u1 + ðk1 + k2 Þu2
t1
p2 ðt Þdu2 gdt ¼ 0 (5)
Because Eq. (5) is valid for any interval ½t1 , t2 , its integrand must be equal
to zero, that is,
½m1 u€1 + k1 ðu1 u2 Þ p1 ðt Þdu1 + ½m2 u€2 k1 u1 + ðk1 + k2 Þu2 p2 ðt Þdu2 ¼ 0
(6)
Inasmuch as the quantities du1 and du2 are arbitrary, Eq. (6) is valid only if
the quantities in the square brackets are equal to zero, that is,
m1 u€1 + k1 ðu1 u2 Þ p1 ðt Þ ¼ 0 (7a)
m2 u€2 k1 u1 + ðk1 + k2 Þu2 p2 ðt Þ ¼ 0 (7b)
which give the equations of motion
m1 u€1 + k1 u1 k1 u2 ¼ p1 ðt Þ (8a)
m2 u€2 k1 u1 + ðk1 + k2 Þu2 ¼ p2 ðt Þ (8b)
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 49
Solution
Inasmuch as the bars are assumed rigid, this system has only two degrees of
freedom. The displaced configuration of the system can be determined from
the two transverse displacements u1 ðt Þ and u2 ðt Þ of points B and C . Referring
to Fig. E1.11, we have
9
f1 ¼ u1 =4L =
f2 ¼ ðu2 u1 Þ=3L (1)
;
f3 ¼ u2 =3L
The displacements of points of application of the forces and the changes of
angles are expressed in terms of the basic quantities u1 and u2 as
9
EE 0 ¼ u1 =4, FF 0 ¼ u1 =2, GG 0 ¼ 3u1 =4 >>
=
HH 0 ¼ u1 + ðu2 u1 Þ=3, QQ 0 ¼ u2 =2
(2)
DfB ¼ f2 f1 ¼ ð4u2 7u1 Þ=12L >
>
;
DfC ¼ f3 + f2 ¼ ð2u2 u1 Þ=3L
The potential energy U due to the deformation of the springs is
1 1 1 1
U ¼ k1 ðEE 0 Þ + k2 ðQQ 0 Þ + k3 ðDfB Þ2 + k4 ðDfC Þ2
2 2
2 2 2 2
which by virtue of Eqs. (2) becomes
1 2 1 2 1 2
U¼ ku 1 + ku 2 + k ð4u2 7u1 Þ2 + k ð2u2 u1 Þ2
32 4 72 9
¼ 0:934ku 21 + 1:361ku 22 1:667u1 u2
Its variation is
dU ¼ k ð1:868u1 1:667u2 Þdu1 + k ð1:667u1 + 2:722u2 Þdu2 (3)
The kinetic energy consists of the kinetic energies T1 and T2 of the bars ΑΒ
and CD, and of the kinetic energy T3 of the rigid body S. Thus, we have
2
1 1 1 d 1
T ¼ IA f_ 1 + ID f_ 3 + m ðHH 0 Þ + IH f_ 2
2 2 2
(4)
2 2 2 dt 2
@CivilMethod
50 PART I Single-degree-of-freedom systems
where
ð4LÞ3 ð3LÞ3 L3
m ¼ mL,
IA ¼ m
, ID ¼ m
, IH ¼ m (5)
3 3 6
Introducing Eqs. (1), (2), (5) into Eq. (4) yields
2 1 1 1
T ¼ m u_ 21 + m u_ 22 + m ðu_ 2 + 2u_ 1 Þ2 + m ðu_ 2 u_ 1 Þ2
3 2 18 108
¼ 0:898m u_ 21 + 0:565m u_ 22 + 0:204u_ 1 u_ 2
we obtain
Z t2 Z t2
dTdt ¼ ½m ð1:796u€1 + 0:204u€2 Þdu1 + m ð0:204u€1 + 1:130u€2 Þdu2 dt
t1 t1
(6)
The nonconservative forces include the loading pðx, t Þ and the damping
forces. Their virtual work is expressed in terms of the basic quantities as follows:
Z 3L x
dWnc p
¼ pðx, t Þ 1 du2 dx
0 3L
Z 3L (7)
x x
¼ p f ðt Þ 1 du2 dx ¼ 1:5
pLf ðt Þdu2
0 L 3L
d d d
D
dWnc ¼ c1 ðFF 0 Þd ðFF 0 Þ c2 ðGG 0 ÞdðGG 0 Þ c3 ðDfB ÞdðDfB Þ
dt dt dt (8)
d
c4 ðDfC Þd ðDfC Þ
dt
Using Eq. (2) and taking into account that c1 ¼ c, c2 ¼ 3c, c3 ¼ c4 ¼ 2cL2 ,
we can write
D
dWnc ¼ cð2:840u_ 1 0:833u_ 2 Þdu1 + cð0:833u_ 1 1:111u_ 2 Þdu2
Hence, we have
dWnc ¼ cð2:395u_ 1 0:833u_ 2 Þdu1 + cð0:833u_ 1 1:111u_ 2 Þdu2 + 1:5
pLf ðt Þdu2
(9)
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 51
Finally, the potential A of the external conservative forces is due to the con-
stant axial force P . Hence it is
A ¼ P ðDD 0 Þ ¼ Pe
and
dA ¼ Pde (10)
The variation de is evaluated as follows.
Referring to Fig. E1.11, we have
e ¼ ðAD Þ ðAD0 Þ ¼ 10L 4L cos f1 3Lcos f2 3L cos f3
Therefore
de ¼ Lð4sin f1 df1 + 3 sin f2 df2 + 3sin f3 df3 Þ
(11)
¼ Lð4f1 df1 + 3f2 df2 + 3f3 df3 Þ
which is introduced into Eq. (10) to yield
7P P P 2P
dA ¼ u1 + u2 du1 + u1 u2 du2 (12)
12L 3L 3L 3L
Introducing the expressions for dU , dT , dWnc , and dA into Hamilton’s prin-
ciple, Eq. (1.7.13), we obtain the following equations of motion
7P
1:796m u€1 + 0:204m u€2 + 2:395cu_ 1 0:833cu_ 2 + 1:868k u1
12L
P
+ 1:667k + u2 ¼ 0
3L
P
0:204m u€1 + 1:130m u€2 0:833cu_ 1 + 1:111cu_ 2 + 1:667k + u1
3L
2P
+ 2:722k u2 ¼ 1:5
pLf ðt Þ
3L
or in the matrix form
" #( ) " #( )
1:796 0:204 u€1 2:395 0:833 u_ 1
m +c
0:204 1:130 u€2 0:833 1:111 u_ 2
" #( ) ( ) (13)
1:868 0:583l 1:667 + 0:333l u1 0
+k ¼
1:667 + 0:333l 2:722 0:667l u2 pLf ðt Þ
1:5
where l ¼ P=kL.
The elastic forces of the system are
fS1 ¼ k ð1:868 0:583lÞu1 + k ð1:667 + 0:333lÞu2
fS2 ¼ k ð1:667 + 0:333lÞu1 + k ð2:722 0:667lÞu2
@CivilMethod
52 PART I Single-degree-of-freedom systems
Z
1 L
∂u ðx, t Þ 2
T¼ m dx (1)
2 0 ∂t
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 53
The strain energy of the beam is obtained by integrating the strain energy
density over its volume V , namely
Z
1
U¼ sx ex dV (2)
2 V
From the beam theory we have
M ðx Þ sx ∂2 u ðx, t Þ
sx ¼ y, ex ¼ , M ðx Þ ¼ EI
I E ∂x 2
Substituting the previous equations into Eq. (2) and integrating over the
cross-section of the beam yield
Z 2 2
1 L ∂ u ðx, t Þ
U¼ EI dx (3)
2 0 ∂x 2
For the simplicity of the expressions, the differentiation with respect to time
t will be designated by an over-dot while that with respect to the spatial coor-
dinated x by a prime. Moreover, the arguments will be dropped for the same
reason. Hence, expressions (1) and (3) can be rewritten as
Z
1 L
T¼ m u_ 2 dx (4)
2 0
Z
1 L
EI ðu 00 Þ dx
2
U¼ (5)
2 0
Their variations are
Z L
dT ¼ _ udx
m ud _ (6)
0
Z L
dU ¼ EI u 00 du 00 dx (7)
0
@CivilMethod
54 PART I Single-degree-of-freedom systems
Introducing Eqs. (5), (9), (10) into Hamilton’s principle, Eq. (1.7.13), we
obtain
Z t2
Z L Z L Z L
EI u 0000 dudx _ udx
m ud _ pðx, t Þdudx dt ¼ 0 (11)
t1 0 0 0
Because Eq. (13) is valid for any interval ½t1 , t2 , the integrand must vanish,
namely
Z L
½EI u 0000 + m u€ pðx, t Þdudx ¼ 0 (14)
0
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 55
(a) (b)
FIG. 1.8.1 Simple (a) and double (b) pendulum.
@CivilMethod
56 PART I Single-degree-of-freedom systems
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 57
The variations associated with the kinetic energy and the potential energy
defined by Eqs. (1.8.2), (1.8.3), respectively, are of the form
∂T ∂T ∂T ∂T
dT ¼ dq1 + ⋯ + dqN + d q_ + ⋯ + d q_
∂q1 ∂qN ∂q_ 1 1 ∂q_ N N
∂A ∂A
dA ¼ dq1 + ⋯ + dqN
∂q1 ∂qN
Substituting these expressions into Eq. (1.8.5), integrating by parts the terms
including d q_ i and taking into account dq1 ¼ dq2 ¼ ⋯ ¼ dqN ¼ 0 at instants t1
and t2 , we obtain
Z t 2
∂T d ∂T ∂A ∂T d ∂T ∂A
dq1 + ⋯ + dqN dt ¼ 0
t1 ∂q1 dt ∂q_ 1 ∂q1 ∂qN dt ∂q_ N ∂qN
Because the time interval ½t1 , t2 as well as the virtual displacements dqi are
arbitrary, this previous equation results in the following equations
d ∂T ∂T ∂A
+ ¼ 0 ði ¼ 1, 2, …, N Þ (1.8.6)
dt ∂q_ i ∂qi ∂qi
which, using Eq. (1.8.4), become
d ∂T ∂T
¼ Qi ði ¼ 1, 2, …, N Þ (1.8.7)
dt ∂q_ i ∂qi
Eq. (1.8.6) or (1.8.7) are the Lagrange equations of motion.
When nonconservative forces act on the system in addition to the conserva-
tive forces, we can include them in Lagrange’s equations, if the work done by
the nonconservative forces riding the virtual displacements is expressed in
terms of the generalized forces, that is,
dWnc ¼ Q1 dq1 + Q2 dq2 + ⋯ + QN dqN (1.8.8)
Introducing Eq. (1.8.8) into Hamilton’s principle, Eq. (1.7.9), the Lagrange
equations (1.8.6) become
d ∂T ∂T ∂A
+ ¼ Qi ði ¼ 1, 2, …, N Þ (1.8.9)
dt ∂q_ i ∂qi ∂qi
The elastic force components, which are derivable from a potential U (strain
energy), can be also involved in Eq. (1.8.9). Noting that
U ¼ U ðq1 , q2 , …, qN Þ (1.8.10)
the associated variation is
∂U ∂U
dU ¼ dq1 + ⋯ + dqN
∂q1 ∂qN
Therefore, the components ∂U =∂qi express generalized elastic forces and
Lagrange’s equations become
@CivilMethod
58 PART I Single-degree-of-freedom systems
d ∂T ∂T ∂V
+ ¼ Qi ði ¼ 1, 2, …, N Þ (1.8.11)
dt ∂q_ i ∂qi ∂qi
where
V ¼U +A (1.8.12)
X
N
dW ¼ Qi dqi (1.8.13)
i¼1
From physical consideration, the work done by the two sets of forces is the
same. The only difference is that they are expressed in two different coordinate
systems. Therefore, we can write
X
N X
K
Qi dqi ¼ Fk dxk (1.8.15)
i¼1 k¼1
or in matrix form
QT dq ¼ FT dx (1.8.16)
where
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 59
@CivilMethod
60 PART I Single-degree-of-freedom systems
∂T
¼ m2 L1 L2 q_ 1 q_ 2 sin ðq1 q2 Þ
∂q2
∂A
¼ m2 gL2 sin q2
∂q2
Applying Eq. (1.8.6) for i ¼ 1, 2 and q1 ¼ q1 q2 ¼ q2 , we obtain the equa-
tions of motion of the double pendulum
ðm1 + m2 ÞL1 q€1 + m2 L2 q€2 cos a + q_ 2 sin a + ðm1 + m2 Þg sin q1 ¼ 0 (4a)
2
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 61
Solution
Because the rod is no more inextensional, the system has two degrees of freedom.
Its displaced configuration can be specified either by the orthogonal coordinates x
and y of the mass or by the angle of the q and the axial deformation of the rod.
The kinetic energy of the system is
1
T ¼ m x_ 2 + y_ 2 (1)
2
The potential energy of the external force (gravitational force) is
A ¼ mgy (2)
and the potential of the elastic force
1
U ¼ ke2 (3)
2
where k ¼ EA=L is the axial stiffness of the rod and e its elongation. The latter
is expressed in terms of x and y as
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
e ¼ x 2 + y2 L (4)
Introducing Eq. (4) in the expression for the axial stiffness, Eq. (3), yields
1 EA pffiffiffiffiffiffiffiffiffiffiffiffiffiffi
ffi 2
U¼ x 2 + y2 L (5)
2 L
Differentiating the energies, we obtain
!
d ∂T ∂A ∂U EA L
¼ m x€, ¼ 0, ¼ 1 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi x (6)
dt ∂x_ ∂x ∂x L x 2 + y2
!
d ∂T ∂A ∂U EA L
¼ m y€, ¼ mg, ¼ 1 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi y (7)
dt ∂y_ ∂y ∂y L x 2 + y2
Introducing Eqs. (6), (7) into Lagrange’s equations (1.8.11), we obtain the
equations of motion of the soft pendulum
!
EA L
m x€ + 1 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi x ¼ 0 (8a)
L x 2 + y2
@CivilMethod
62 PART I Single-degree-of-freedom systems
!
EA L
m x€ + 1 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi y ¼ mg (8b)
L x 2 + y2
1 1
_ 2 m y U_ W _
_ + m xc V W
T ¼ m U_ 2 + V 2 + Io W
c (1)
2 2
The potential energy U consists of the strain energy of all columns. For the i
column it is
1 h i i 2 i 2 i 2 i
Ui ¼ k11 u + k22i
v + k33i
w
2
or using matrix notation
i 8 9
k 0 0 < u i =
1 i i i
11 i 1 T
Ui ¼ u v w 0 k22 0 v i ¼ Di ki Di (2)
2 0 0 k i wi : ; 2
33
∂T
¼0 (7d)
∂V
d ∂T € m y U€
+ m xc V€
¼ Io W (7e)
_
dt ∂W c
∂T
∂W ¼0 (7f)
@CivilMethod
64 PART I Single-degree-of-freedom systems
∂U
¼ k11 U + k12 V + k13 W (8a)
∂U
∂U
¼ k21 U + k22 V + k23 W (8b)
∂V
∂U
∂W ¼ k 31 U + k 32 V + k 33 W (8c)
M € + K
U U ¼P
ðt Þ (9)
where
2 3 2 3
m 0 m yc k 11 k12 k13
6 7 6 7
¼6 0
M m mx c 7 ¼ 6 k21
K k22 k23 7
4 5, 4 5,
my c mx c Io k31 k32 k33
8 9
> P ðt Þ >
< x =
P ðt Þ ¼
P y ðt Þ
>
: >
;
y A P x ðt Þ + xA P y ðt Þ
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 65
where
dgj
aji ¼ (1.8.27)
dqi
If we assume that the constraints are frictionless, then no work is done by the
constraint forces Ri when they ride any virtual displacement dqi , that is,
X
K
Ri dqi ¼ 0 (1.8.28)
i¼1
@CivilMethod
66 PART I Single-degree-of-freedom systems
where we note that a separate equation is written for each of the constraints.
Next, we subtract the sum of equations of the form (1.8.29) from
Eq. (1.8.28) and, interchanging the order of summation, we obtain
!
XK X n
Ri lj aji dqi ¼ 0 (1.8.30)
i¼1 j¼1
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 67
Solution
The kinetic and the potential energies of the system are
1
T ¼ m x_ 2 + y_ 2
2
A ¼ mgy
U ¼0
Because the rod is rigid, the coordinates must satisfy the constraint equation
g1 ¼ x 2 + y 2 L2 ¼ 0
Differentiating the quantities T and A we obtain
d ∂T ∂T ∂A ∂g1
¼ m x€, ¼ 0, ¼ 0, ¼ 2x, Q1 ¼ px
dt ∂x_ ∂qi ∂x ∂x
d ∂T ∂T ∂A ∂g1
¼ m y€, ¼ 0, ¼ mg, ¼ 2y, Q2 ¼ py
dt ∂y_ ∂y ∂y ∂y
Applying Eq. (1.8.32) for q1 ¼ x and q2 ¼ y we obtain the equations of
motion
m x€ ¼ px + 2xl (1a)
m y€ + mg ¼ py + 2yl (1b)
x 2 + y2 l 2 ¼ 0 (1c)
Eqs. (1) must be solved for the three unknowns x, y, and l. It should be
noted that two of these equations are differential and one algebraic and therefore
special care is required for their solution. A convenient method is to differen-
tiate the constraint equation twice with respect to time and then to solve the
€ y€ and the parameter
resulting linear system of equations for the accelerations x,
l. For the problem at hand, we obtain
2T
x x€ + y y€ ¼ (2)
m
Eqs. (1a), (1b), (2) are combined and written in matrix form
2 38 9 8 9
m 0 2x < x€ = < px =
4 0 m 2y 5 y€ ¼ py mg (3)
: ; : ;
x y 0 l 2T =m
which are solved to yield
px py mg
L2 x€ + x x_ 2 + y_ 2 ¼ y 2 xy (4a)
m m
2 px 2 py mg
L y€ + y x_ + y_ ¼ xy + x
2 2
(4b)
m m
@CivilMethod
68 PART I Single-degree-of-freedom systems
xpx y py mg
L l¼
2
T (4c)
2 2
Eqs. (4a), (4b) are solved using techniques for nonlinear differential equa-
tions. Analytical solutions are in general out of the question. However, a numer-
ical solution is always feasible using the methods presented in Chapter 5. Once
the coordinates x ðt Þ, y ðt Þ and the Lagrange multiplier l have been established,
they are utilized in Eq. (1.8.31) to evaluate the constraint forces, which are the
components of the axial force of the rod. Thus, we have
∂g ∂g qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
Rx ¼ l ¼ 2lx, Rx ¼ l ¼ 2ly, S ¼ R2x + R2y ¼ 2lL (5)
∂x ∂y
for all virtual displacements dxi consistent with the constraints, which are
assumed workless and bilateral.
Inasmuch as the forces are conservative, they are derivable from a potential
function V ¼ V ðx1 , x2 , …, x3N Þ, V ¼ U + A, according to the relation
∂V
Fi ¼ (1.8.35)
∂xi
Using Eq. (1.8.35), Eq. (1.8.34) is written as
X
3N
∂V
dW ¼ dxi ¼ 0 (1.8.36)
i¼1
∂xi
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 69
Substituting the previous expression for dxi into Eq. (1.8.36), we obtain
3N X
X n
∂V ∂xi
dW ¼ dqj ¼ 0 (1.8.39)
i¼1 j¼1
∂xi ∂qj
Noting that
∂V X 3N
∂V ∂xi
¼ (1.8.40)
∂qj i¼1
∂xi ∂qj
Because dqj are assumed to be independent, the virtual work is zero only if
the coefficients of dqj are zero at the equilibrium condition, that is, if
∂V
¼ 0, j ¼ 1, 2, …, n (1.8.42)
∂qj 0
The subscript zero denotes that the derivatives refer to the equilibrium
position.
Let us expand now the potential energy function V ðq1 , q2 , …, qn Þ in a Tay-
lor series about the position of equilibrium
Xn n 2
∂V 1X n X
∂ V
V ¼ V0 + dqi + dqi dqj + ⋯ (1.8.43)
i¼1
∂qi 0 2 i¼1 j¼1 ∂qi ∂qj 0
We can arbitrarily set the potential energy at the reference position equal to
zero, that is,
V0 ¼ 0 (1.8.44)
@CivilMethod
70 PART I Single-degree-of-freedom systems
If we now assume that the displacements about the equilibrium position are
small, we can neglect terms of order higher than the second in Eq. (1.8.43).
Thus, using Eqs. (1.8.41), (1.8.44) the expression for the potential energy is sim-
plified as
n 2
1X n X
∂ V
V¼ qi qj (1.8.45)
2 i¼1 j¼1 ∂qi ∂qj 0
or setting
∂2 V
kij ¼ kji ¼ (1.8.46)
∂qi ∂qj 0
The quantities kij defined by Eq. (1.8.46) are the stiffness coefficients of the
system. Thus we see that the potential energy is expressed by a homogeneous
quadratic function of the generalized coordinates qi if small motions about the
position of equilibrium are examined.
Eq. (1.8.47) is written in matrix form
1
V ¼ qT kq (1.8.48)
2
where
8 9 2 3
>
> q1 >
> k11 k12 ⋯ k1n
< = 6 k21 k22
q2 ⋯ k2n 7
q¼ , k ¼6
4 ⋮ ⋮
7 (1.8.49)
>
> ⋮> > ⋱ ⋮ 5
: ;
qn kn1 kn2 ⋯ knn
The matrix k is called the stiffness matrix of the system.
The expression for the potential energy given in Eq. (1.8.47) is an example
of a quadratic form. For a system whose reference equilibrium configuration is
stable, the potential energy V is positive for all possible values of qi , except
q1 ¼ q2 ¼ … ¼ qn ¼ 0. In this case, the function V is referred to as positive def-
inite. This condition, however, puts restrictions on the allowable values of kij . It
is clear that all diagonal elements must be positive. The necessary and sufficient
condition that V be positive definite is that
2 3
k11 k12 ⋯ k1n
k k 6 k21 k22 ⋯ k2n 7
k11 > 0, 11 12 > 0, …, 6
4 ⋮ ⋮ ⋱ ⋮ 5>0
7 (1.8.50)
k21 k22
kn1 kn2 ⋯ knn
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 71
∂xk X n
∂xk
¼ q_ (1.8.52)
∂t j¼1
∂qj j
Introducing this expression into Eq. (1.8.51) we can write the kinetic energy
in the form
1X n X n
T¼ mij q_ i q_ j (1.8.54)
2 i¼1 j¼1
X
3N
∂xk ∂xk
mij ¼ mji ¼ mk (1.8.55)
k¼1
∂qi 0 ∂qj 0
The quantities mij defined by Eq. (1.8.55) are the inertia coefficients of the
system.
Eq. (1.8.54) is written in matrix form
1
T ¼ q_ T mq_ (1.8.56)
2
where
8 9 2 3
>
> q_ 1 >
> m11 m12 ⋯ m1n
< = 6 m21 m22
q_ 2 ⋯ m2n 7
q_ ¼ , m ¼6
4 ⋮
7 (1.8.57)
>
> ⋮> > ⋮ ⋱ ⋮ 5
: ;
q_ n mn1 mn2 ⋯ mnn
The matrix m is called the mass matrix of the system. The kinetic energy is a
positive definite quadratic function because it is the sum of positive quantities,
that is, the kinetic energies of the masses of the individual particles.
@CivilMethod
72 PART I Single-degree-of-freedom systems
∂V X n
¼ kij qj
∂qi j¼1
or in matrix form
m€
q + kq ¼ pðtÞ (1.8.59)
where p(t)¼Q
The matrices m and k are symmetric. It is an advantage of the Lagrange for-
mulation of the equations of motion that it preserves the symmetry of the coef-
ficient matrices for those cases where T and V are represented by quadratic
functions of the velocities and displacements, respectively.
where cij ¼ cji are the damping coefficients of the linear viscous damping.
Apparently, we can construct a quadratic function
1X n X n
R¼ cij q_ i q_ j (1.8.61)
2 i¼1 j¼1
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 73
which yields
∂R Xn
fDj ¼ ¼ cij q_ j (1.8.62)
∂q_ j i¼1
(a) (b)
FIG. 1.9.1 Influence of the gravity load.
@CivilMethod
74 PART I Single-degree-of-freedom systems
The elongation ust of the spring under its own weight will be
ust ¼ W =k ¼ constant (1.9.2)
Further we set
u ¼ ust + uðt Þ (1.9.3)
where uðt Þ represents the vertical displacement measured from the position of
the static equilibrium.
Differentiating Eq. (1.9.3) yields
_ u€ ¼ u€
u_ ¼ u, (1.9.4)
Using Eqs. (1.9.3), (1.9.4), the equation of motion (1.9.1) becomes
m u€
+ cu_ + ku st + k u ¼ pðt Þ + W
or using Eq. (1.9.2) we obtain
m u€
+ cu_ + k u ¼ pðt Þ (1.9.5)
The conclusion drawn from Eq. (1.9.5) states that, in the study of the
dynamic response of a system undergoing small displacements, the loads due
to gravity can be neglected. Of course, the total displacements will result as
the sum of the static plus dynamic displacements. That is, the superposition
principle is valid.
1.10 Problems
Problem P1.1 The plane square rigid body B of side length L and surface
mass density g is supported by two identical inclined columns having
cross-sectional moment of inertia I , modulus of elasticity E, and negligible
mass. Derive the equation of motion neglecting the axial deformation of
the columns (Fig. P1.1).
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 75
Problem P1.2 Consider the structure of Fig. P1.2a. The square plate of constant
thickness h ¼ a=10 and mass density g is supported at its center by a flexible
column having a circular cross-section with diameter d ¼ a=10, height a, and
material constants E, n. The plate is loaded by a force P acting in the plane
of the plate at point (Aða=8, a=6Þ and in the direction ∡x, P ¼ b ¼ 30° as
shown in Fig. P1.2b. Derive the equations of motion of the plate when the mass
of the column is neglected.
(a) (b)
FIG. P1.2 Structure in problem P1.2
Problem P1.3 The semicircular rigid plate of constant thickness and total mass
m is supported as shown in Fig. P1.3. Taking into account that the support at
point O is a hinge, formulate the equation of motion of the plate using (i)
the method of equilibrium of forces, (ii) the principle of virtual displacements,
and (iii) the method of the Lagrange equations.
Problem P.1.4 Consider the system shown in Fig. P1.4. The bars AD and EG
are rigid with masses m and m=3 , respectively. The mass at end D is concen-
trated. The elastic supports at points at B, E , and D are simulated by springs
with a stiffness k while the end G is supported by a viscous damper with a
damping coefficient c. The rod CE is weightless and rigid. Derive the equation
of motion using the principle of virtual displacements.
@CivilMethod
76 PART I Single-degree-of-freedom systems
Problem P1.5 Consider the system shown in Fig. P1.5. The mass m is sup-
ported at the top of the flexible and massless column 2 3, which is supported
on the ground by means of the rigid body 1 2 of mass 2a m. The support 1 is
elastically restrained by the rotational spring CR . Formulate the equation of
motion of the structure using CR ¼ EI =2a, m ¼ m=a.
Problem P1.6 Consider the two-story frame of Fig. P1.6. The columns 1 2,
10 20 , and the beam 3 30 are rigid while the columns 2 3, 20 30 , and the
beam 2 20 are massless and flexible with cross-sectional moment of inertia
I and modulus of elasticity E. The supports at 1 and 10 are elastically restrained
by rotational springs with a stiffness CR . Formulate the equation of motion of
the structure taking CR ¼ EI =2a and m ¼ m=a.
Problem P1.7 The system of Fig. P1.7 consists of the beam AB and the rigid
body S interconnected at B. The beam AB has a negligible mass, modulus of
elasticity E , and cross-sectional moment of inertia I . The beam is fixed at A
while the rigid body is elastically restrained at C by a rotational spring with
a stiffness CR ¼ EI =10L. The total mass m is uniformly distributed. The system
is loaded by the concentrated moment M ðt Þ at point B. Derive the equation of
motion of the system using Lagrange’s equations.
10 10
Problem P1.8 The frame of Fig. P1.8 consists of the rigid beam BD of total
mass m and the two massless and flexible columns AB and CD with a
cross-sectional moment of inertia I and modules of elasticity E. The two mass-
less cables FB and GD have cross-sectional area A and cannot undertake com-
pressive force. Derive the equation of motion of the structure taking
I =A ¼ a2 =25 and m ¼ m=5a.
Problem P1.9 Consider the two-story frame of Fig. P1.9. The columns of the
frame are rigid and have a surface mass density g. Their elastic support on the
ground is simulated by the rotational springs with a stiffness CR ¼ EI =10a.
The horizontal beams are massless and flexible with a cross-sectional moment
of inertia I and modulus of elasticity E. Derive the equation of motion when the
structure is subjected to the horizontal loads pðt Þ at the beam levels.
@CivilMethod
78 PART I Single-degree-of-freedom systems
Problem P1.10 The hinge O of the soft pendulum of Fig. P1.10 is elastically
restrained by the rotational spring with a stiffness CR ¼ EAL=10. The length of
the rod is L, its cross-sectional area A, and the modulus of elasticity E. Formu-
late the equation of motion of the pendulum.
Problem P1.11 The rigid bar AB of circular cross-section and mass density
m ¼ m=a is hinged at point A (Fig. P1.11). The cables DB,FB have cross-
sectional area A and modulus of elasticity E. They are assumed massless
and are prestressed so that they can undertake compressive forces. Formulate
the equation of motion of the structure taking into account that the load P is
removed suddenly at instant t ¼ 0. Evaluate the minimum prestressing force
of the cables DB, FB so that they can undertake compressive loads.
Problem P1.12 Consider the structure of Fig. P1.12. The column AC has a
circular cross-section and a mass per unit length m ¼ m=a; it is supported by
a spherical hinge on the ground and is kept in place by three elastic cables of
cross-sectional area A and modulus of elasticity E. The cables are assumed mass-
less and are prestressed so that they can undertake compressive force. Derive the
equation of motion of the structure when it is loaded by the horizontal force P ðt Þ
acting at the top of the column in the direction ∡x,P ¼ b (Fig. P1.12b).
(a) (b)
FIG. P1.12 Structure in problem P1.12
Problem P1.13 The silo of Fig. P1.13 is supported on its fundament by four
identical columns of a square cross-section. The silo is full of material of density
g. The ground yields elastically with a subgrade constant Ks . The silo and the
fundament are rigid. Derive the equation of motion of the structure when it is
loaded by the horizontal force P ðt Þ acting at the top of the silo in the direction
(a) (b)
FIG. P1.13 Silo on elastic subgrade. (a) Vertical section. (b) Plan form.
@CivilMethod
80 PART I Single-degree-of-freedom systems
∡x,P ¼ b (Fig. P1.13b) using the following data: Side of the columns a=4;
thickness of the bottom and walls of the silo a=8; density of the material of
the silo 1:5g; and soil constant Ks ¼ EI =1500a 3 .
Problem P1.14 Consider the one-story building of Fig. P1.14. The rigid plate is
an equilateral triangular with a side a and it is supported by three columns of
height a, rectangular cross-section a=10 a=20, and modulus of elasticity E.
The columns are fixed at both ends. Derive the equation of motion of the plate
when a horizontal force P ðt Þ acts at point Að0, a=5Þ in the direction ∡x, P ¼ b.
The dead weight of the plate is included in p (kN=m2 ).
Rigid plate
Problem P1.15 The two one-story buildings of Fig. P1.15 are connected with a
beam as shown in the figure. All columns have a square cross-section with a
moment of inertia Ic ¼ 2I . The connecting beam has a square cross-section with
moment of inertia Ib ¼ I . The structure is loaded by the horizontal force F(t) at the
level of the plates as shown in Fig. P1.15b. Formulate the equations of motion
using Lagrange’s equations. Assume: Torsion constant It ¼ 2:25d 4 =16, d ¼side
length of the square cross-section of the beam.
Rigid plate
(a)
beam
(b)
FIG. P1.15 Structure in problem P1.15. (a) vertical section, (b) plan form.
Problem P1.16 The system of Fig. P1.16 consists of the block of mass m1 ,
which can slide without friction on the inclined surface, and the pendulum of
@CivilMethod
General concepts and principles of structural dynamics Chapter 1 81
length L and mass m2 , which is pivoted at the center of mass of the block. The
rod of the pendulum has a cross-sectional area A and modulus of elasticity E.
Assuming plane motion, derive the equation of motion of the system taking
EA=L ¼ 5k and m1 ¼ 5m2 .
.
FIG. P1.16 System in problem P1.16
Problem P1.18 Consider the crane of Fig. P1.18. The horizontal beam is
assumed rigid. The column is flexible with a cross-sectional moment of inertia
I and the cable axially deformable with cross-sectional area A. The mass of the
cable and column is negligible. Derive the equation of motion of the system
when it is loaded by the horizontal force pðt Þ in the plane of the structure using
I =A ¼ a 2 =100 and a common modulus of elasticity E.
@CivilMethod
82 PART I Single-degree-of-freedom systems
@CivilMethod
Chapter 2
Single-degree-of-freedom
systems: Free vibrations
Chapter outline
2.1 Introduction 83 2.3.3 Overdamped system 96
2.2 Free undamped vibrations 83 2.4 Conservation of energy in an
2.3 Free damped vibrations 91 undamped system 97
2.3.1 Critically damped system 91 2.5 Problems 99
2.3.2 Underdamped system 92 References and further reading 103
2.1 Introduction
In this chapter, the free vibrations of a single-degree-of-freedom system
(SDOF) are studied, that is, its response when it is not subjected to any external
force, pðt Þ ¼ 0, but it is excited by an initial displacement and/or initial velocity.
The dynamic model of the system is shown in Fig. 1.4.1 and the equation of
motion (1.4.8) takes the form
m u€ + cu_ + ku ¼ 0 (2.1.1)
Eq. (2.1.1) is an ordinary linear homogeneous differential equation of the
second order with constant coefficients and its solution can be obtained using
known mathematical methods. Inasmuch as we are interested in the physical
response of the system described by this equation, it is advisable to analyze
the free vibration response in two stages, first for c ¼ 0 and then c 6¼ 0. In the
first case, we speak of free undamped vibrations while in the second case we
speak of free damped vibrations. Illustrative examples analyzing the free vibra-
tions of SDOF systems are presented. The pertinent bibliography with recom-
mended references for further study is also included.
u ¼ elt (2.2.2)
where l is an arbitrary constant to be determined. Substitution of Eq. (2.2.2)
into Eq. (2.2.1) gives
2
ml + k elt ¼ 0 (2.2.3)
ml2 + k ¼ 0 (2.2.4)
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 85
u_ 0
B¼ (2.2.12b)
w
and Eq. (2.2.10) becomes
u_ 0
u ðt Þ ¼ u0 cos wt + sin wt (2.2.13)
w
Obviously, it is u ðt Þ ¼ 0, when u0 ¼ u_ 0 ¼ 0. Hence, the system is set
to motion only if it is given an initial displacement and/or an initial
velocity.
We set
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
2ffi
u_ 0
r ¼ ðu 0 Þ2 + (2.2.14)
w
Inasmuch as it is
2 2
u0 u_ 0 u0 u_ 0
1, 1 and + ¼1 (2.2.16)
r rw r rw
we can set
u0 u_ 0
¼ cos q and ¼ sin q (2.2.17)
r rw
and Eq. (2.2.15) becomes
u ðt Þ ¼ rcos ðwt qÞ (2.2.18)
where
u_ 0
q ¼ tan 1 (2.2.19)
wu0
Eq. (2.2.18) states that the motion of the system is a harmonic vibration
with amplitude ju ðt Þjmax ¼ r, angular velocity w, and phase angle q. The
@CivilMethod
86 PART I Single-degree-of-freedom systems
The time required for the undamped system to complete one cycle of free
vibration is referred to as the natural period of vibration of the system, which
is denoted by T and measured in seconds. It is related to the natural frequency of
vibration through
2p
T¼ (2.2.20)
w
The inverse of the period
1
f¼ (2.2.21)
T
expresses the number of cycles that the system performs in 1 s. This is referred
to as the natural cyclic frequency. The unit of f is the hertz (Hz) (cycles per
second, cps) and it is related to w through
w
f¼ (2.2.22)
2p
The displacement versus time for a system with w ¼ 8s1 , u0 ¼ 0:05m, and
u_ 0 ¼ 1m=s is shown in Fig. 2.2.2.
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 87
0.2
(du/dt)0 T = 2p /w
0.15
0.1
0.05
u0
0
–0.05
-
–0.1
–0.15
T = 2p /w
–0.2
0 0.5 1 1.5 2
FIG. 2.2.2 Response of an undamped SDOF system.
(a) (b)
FIG. E2.1 One-story shear building in Example 2.2.1. (a) Vertical section. (b) Plan form.
@CivilMethod
88 PART I Single-degree-of-freedom systems
Solution
The motion of the slab is described by the two displacements U ,V of its
center O along the x,y axes, respectively, and its rotation W about O. The
stiffness and mass matrices of the structure can be established using
Eqs. (17a), (17b) for the single-story building in Example 1.5.6. However,
taking into account that the structure is symmetric with respect to both axes,
the three components of the motion are uncoupled and an ad hoc solution can
be readily obtained.
(i) The stiffness of the columns in the x and y directions are:
12EI 1y
kx1 ¼ ky1 ¼ ¼ 1360:8kN=m (1a)
h3
12EI 2y
kx2 ¼ ¼ 907:2kN=m (1b)
h3
12EI 2x
ky2 ¼ ¼ 403:2kN=m (1c)
h3
Hence the respective stiffnesses of the structure are
Kx ¼ 4kx1 + 2kx2 ¼ 7257:6kN=m (2a)
Ky ¼ 4ky1 + 2ky2 ¼ 6249:6kN=m (2b)
The torsional stiffness KW is equal to the moment produced by the elastic
forces of the columns for unit rotation of the slab. Referring to Fig. E2.2,
we have
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 89
fx i ¼ kxi yi W (4a)
fy i ¼ kyi xi W (4b)
GI it E
M i ¼ kWi W, kW i ¼ , G¼ (5)
h 2ð 1 + n Þ
The moment of the elastic forces of column i with respect O is
M i ¼ yi2 kxi + xi2 kyi + kWi W (6)
GI 1t GI 2
kW1 ¼ ¼ 1998:7kNm, kW 2 ¼ t ¼ 823:2kNm (8)
h h
x1 ¼ x6 ¼ 4:85m, x2 ¼ x5 ¼ 0, x3 ¼ x4 ¼ 4:85m
y1 ¼ y2 ¼ y3 ¼ 2:85m, y4 ¼ y5 ¼ y6 ¼ 2:85m
X
KW ¼ yi2 kxi + xi2 kyi + kWi ¼ 1:9663 105 kNm (9)
i¼1
The mass of the slab and its moment of inertia with respect to O are
10 53 + 103 5
Io ¼ 70=9:81 ¼ 3716:44kNs2 (11)
12
@CivilMethod
90 PART I Single-degree-of-freedom systems
rffiffiffiffiffiffi 9
Kx 1 >
w1 ¼ ¼ 4:510s > >
>
>
m >
>
rffiffiffiffiffiffi >
=
Ky 1
w2 ¼ ¼ 4:185s > (12)
m >
>
rffiffiffiffiffiffiffi >
>
>
KW 1 >
w3 ¼ ¼ 7:274s > ;
I0
Fig. E2.3 presents the displacements of the column at the upper right
corner (x1 ¼ 4:85m, y1 ¼ 2:85m)
0.2
u1(t)
0.15 v1(t)
w1(t)
0.1
0.05
–0.05
–0.1
–0.15
–0.2
0 1 2 3 4 5
t
FIG. E2.3 Displacements u ðt Þ,v ðt Þ and rotation wðt Þ of the top cross-section of column at
1 ¼ 4:85m, y1 ¼ 2:85m in Example 2.2.1.
@CivilMethod
x
Single-degree-of-freedom systems: Free vibrations Chapter 2 91
0.07
0.06 (du/dt)0
0.05
0.04
u(t)
0.03
u0
0.02
0.01
0
0 0.2 0.4 0.6 0.8 1
t
FIG. 2.3.1 Response of a system with critical damping.
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 93
or
l1,2 ¼ xw iwD (2.3.12)
where
qffiffiffiffiffiffiffiffiffiffiffiffi
wD ¼ w 1 x2 (2.3.13)
@CivilMethod
94 PART I Single-degree-of-freedom systems
dynamic response of a SDOF system for various values of the damping ratio and
u0 ¼ 1, u_ 0 ¼ 7ms1 , w ¼ 6s1 . From the plot of Eq. (2.3.19), we conclude that
the motion of an underdamped system is harmonic vibration with frequency
wD and period T ¼ 2p=wD , whose amplitude, however, decays exponentially
with time and is bounded by the envelops rexwt . Hence, wD is called the
damped frequency and the respective period the damped period of the under-
damped system.
1.5 T = 2p /wD
T T
1
re−xw t
0.5 u0
u(t)
–0.5
−re-xw t
–1
–1.5
–2
0 1 2 3 4 5
t
FIG. 2.3.2 Response of an underdamped SDOF system.
2
x=0
1.5 x=0.03
x=0.06
x=0.1
1
0.5
u(t)
–0.5
–1
–1.5
–2
0 1 2 3 4 5
t
FIG. 2.3.3 Response of an underdamped SDOF system for various values of the damping ratio.
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 95
pffiffiffiffiffiffiffi2ffi
2npx
¼ u ðt Þe 1x
from which we obtain
u ðt Þ 2npx
d ¼ ‘n ¼ pffiffiffiffiffiffiffiffiffiffiffiffi (2.3.24)
u ðt + nT Þ 1 x2
The quantity d defined by Eq. (2.3.24) is called the logarithmic decrement
and can be employed to determine the damping ratio x experimentally, when the
displacements u ðt Þ and u ðt + nT Þ between n consecutive cycles are measured.
Example 2.3.1 Free damped vibrations of a silo
Fig. E2.4 presents the idealization of a silo. It consists of two massless and inex-
tensional columns of cross-sectional area a a and a square rigid plate of mass
m. At time t ¼ 0, the silo is displaced horizontally by a force P. Then the force is
suddenly removed and the system starts to vibrate. For P ¼ 200 kN, h ¼ 5m,
a ¼ 0:3m, E ¼ 2.1 107 kN/m2, m ¼ 100kN=ms2 determine:
(i) The damping ratio x of the system, if the horizontal displacement is reduced
to u1 ¼ 0:1u0 after n ¼ 5 oscillations.
(ii) Determine the displacement at time t1 ¼ 2.
@CivilMethod
96 PART I Single-degree-of-freedom systems
Solution
(i) The only possible motion of the structure is its horizontal displacement
u ðt Þ. Hence, the system has a single degree of freedom. The stiffness k
of the structure is due to the relative displacement of the column ends,
which yield the stiffness
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 97
The arbitrary constants A,B are evaluated from the initial conditions
u ð0Þ ¼ u0 and u_ ð0Þ ¼ u_ 0 . Then, Eq. (2.3.29) becomes
u_ 0 + u0 xw
u ðt Þ ¼ exwt u0 cosh Wt + sinh Wt (2.3.30)
W
Eq. (2.3.30) has been plotted in Fig. 2.3.4 with u0 ¼ 1m, u_ 0 ¼ 10ms1 , and
w ¼ 6s1 . It becomes evident that the motion of the overdamped system is
nonoscillatory.
1.5
x =1
x =1.5
x =2
x =2.5
1
u(t)
0.5
0
0 0.5 1 1.5 2
t
FIG. 2.3.4 Response of an overdamped SDOF system for various values of the damping ratio.
@CivilMethod
98 PART I Single-degree-of-freedom systems
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 99
2.5 Problems
Problem P2.1 The structure of Fig. P2.1 consists of two identical rigid bars BA
and BC , both having line density m (mass/length). The support A is a hinge
while the support C is a simple support. The bracing rod DF has a
@CivilMethod
100 PART I Single-degree-of-freedom systems
Problem P2.4 A SDOF system of mass m and stiffness k performs free vibra-
tions. At the end of four complete cycles, the displacement is u ð0Þ=3. If the mass
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 101
Problem P2.6 The horizontal force P applied to the structure of Fig. P2.6 is
suddenly removed at t ¼ 0. Determine the motion of the structure when (i)
the cables are free of any pretension and (ii) have been prestressed to withstand
compression, and compute the minimum required pretension forces. In both
cases, the cables are assumed massless. Data: a ¼ 2m, P ¼ 100kN,
m ¼ 100 kNm1 s2 , E ¼ 2.1 107 kN/m2, I ¼ 880cm4 , A ¼ 5cm2
Problem P2.7 In the structure of Fig. P2.7, the rigid rod AB of circular cross-
section and line mass density m ¼ m=a is supported on the ground by a spher-
ical hinge and held in vertical position by three identical elastic massless cables
of cross-sectional area A and modulus of elasticity E. The cables are prestressed
@CivilMethod
102 PART I Single-degree-of-freedom systems
so that they can withstand compression. Their anchor points D,G,F form an
equilateral triangle. Determine the motion of the structure if the horizontal
force P applied at the top of the rod in the direction of the y axis is suddenly
removed at t ¼ 0. Evaluate the minimum values of the prestress forces, which
ensure the capability of the cables to withstand compression. Data: a ¼ 2m,
m ¼ 100 kN m1s2, A ¼ 5cm2 , E ¼ 2.1 108 kN/m2, P ¼ 100kN.
(a) (b)
FIG. P2.7 Structure in Problem P2.7. (A) Vertical view. (B) Plan form.
Problem P2.8 In the pendulum of Fig. P2.8, the rigid rod suspending the con-
centrated mass m has a line mass density m ¼ 2m=L. The hinge O is elastically
restrained by the rotational spring with stiffness CR ¼ kL2 =2. The rod is
displaced by an angle q0 from the vertical position and then is left to move.
Considering small displacements, derive the equation of motion and compute
its period.
Problem P2.9 The rigid silo of Fig. P2.9 is supported on its fundament by four
identical columns of square cross-sections with a side length a=8. The silo is full
of material with density g. The bottom and the walls of the silo have a thickness
a=8 and material density 1:5g. Compute the frequencies and the periods of the
structure.
@CivilMethod
Single-degree-of-freedom systems: Free vibrations Chapter 2 103
(a) (b)
FIG. P2.9 Silo on four columns. (a) Vertical section. (b) Plan form.
@CivilMethod
Chapter 3
Single-degree-of-freedom
systems: Forced vibrations
Chapter outline
3.1 Introduction 105 3.5.1 Rectangular pulse load 126
3.2 Response to harmonic loading 106 3.5.2 Triangular pulse load 128
3.2.1 Response of undamped 3.5.3 Asymmetrical triangular
systems to harmonic pulse load 131
loading 106 3.5.4 Response to piecewise
3.2.2 Response of damped linear loading 135
systems to harmonic 3.6 Response to a periodic loading 137
loading 110 3.6.1 Periodic loads 137
3.3 Response to arbitrary dynamic 3.6.2 Fourier series 138
loading—Duhamel’s integral 113 3.6.3 Response of the SDOF
3.3.1 Undamped vibrations 113 system to periodic
3.3.2 Damped vibrations 116 excitation 143
3.4 Analytical evaluation of the 3.7 Response to unit impulse 146
Duhamel integral-applications 117 3.7.1 The delta function or
3.4.1 Response to step Dirac’s delta function 146
function load 117 3.7.2 Response to unit impulse 148
3.4.2 Response to ramp 3.7.3 Response to arbitrary
function load 120 loading 151
3.4.3 Response to step function 3.7.4 The reciprocal theorem
load with finite in dynamics 151
rise time 121 3.8 Problems 152
3.5 Response to impulsive loads 125 References and further reading 157
3.1 Introduction
In this chapter, the forced vibrations of the SDOF system are studied. The
dynamic model of the system is shown in Fig. 1.4.1 and the motion of the system
is governed by Eq. (1.4.8), namely
m u€ + cu_ + ku ¼ pðt Þ (3.1.1)
where pðt Þ represents an arbitrary function of time. First, the response under a
harmonic load is examined. This type of loading is particularly important in the
dynamic analysis of structures because it allows understanding the major
differences between the static and dynamic response and identifying phenom-
ena such as resonance that are not conceived by the static consideration. More-
over, the harmonic load analysis allows studying the response of SDOF systems
under a general periodic load using the Fourier series representation of a peri-
odic load. Then, the response of a SDOF system under an arbitrary load is stud-
ied using Duhamel’s integral. Finally, the response to a unit load is discussed by
exploiting the properties of the Dirac delta function. The chapter ends by pre-
senting the dynamic reciprocal theorem. Throughout the chapter, illustrative
examples analyzing the forced vibrations of SDOF systems are presented.
The pertinent bibliography with recommended references for further study is
also included. The chapter is enriched with problems to be solved aiming at
better understanding the theoretical issues.
u ¼ uh + up (3.2.3)
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 107
@CivilMethod
108 PART I Single-degree-of-freedom systems
Obviously
ust ¼ p0 =k (3.2.11)
denotes the static displacement that would be produced by a load p0 , equal to
the amplitude of the harmonic excitation, if it were to be applied statically.
The time-dependent quantity
u ðt Þ 1
Rðt Þ ¼ ¼ b sin wt Þ
ð sin wt (3.2.12)
ust 1 b 2
is called the response ratio. It is dimensionless and expresses the number that
must multiply the static displacement at time t to obtain the respective dynamic
displacement. The response ratio provides a measure of the influence of the
dynamic loading.
The extreme value
D ¼ max jRðt Þj (3.2.13)
is referred to as the dynamic magnification factor (DMF). It is a very useful
quantity in dynamic analysis because, if it has been established for a given
loading, the extreme state of deformation and stress can be obtained by static
analysis.
Eq. (3.2.12) for b ¼ 1 takes the indeterminate form
0
Rðt Þ ¼
0
whose limit can be determined using the L’H^
opital rule. Thus, on the basis of
Eq. (3.2.6b) we obtain
b sin wt
sin wt
‘im Rðt Þ ¼ ‘im
b!1 b!1 1 b2
sin bwt b sin wt
¼ ‘im
b!1 1 b2
(3.2.14)
wt cos bwt sin wt
¼ ‘im
b!1 2b
wt cos wt sin wt
¼
2
From the latter relation, it is concluded that when b tends to 1, that is, when
the excitation frequency w of the harmonic force approaches the natural fre-
quency of the system, the dynamic displacement grows indefinitely with time,
although the amplitude of the harmonic loading is finite. This phenomenon is
called resonance. The growth of the amplitude of the displacement with time
due to resonance is shown in Fig. 3.2.1. The response is periodic with a period
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 109
2p=w. A measure of the growth rate can be obtained by taking the difference
of the amplitudes of two consecutive peaks.
The time peak occurs when
20
w t/2
15
2p /w 2p /w
10
5
R(t)
–5
p
–10 p
–15
–w t/2
–20
0 1 2 3 4 5
t
FIG. 3.2.1 Response ratio of an undamped system at resonance, D ¼ max jRðt Þj ! 1
when t ! 1 (w ¼ 7, b ¼ 1).
dRðt Þ w2 sin wt
¼ ¼0 (3.2.15)
dt 2
or
np
t¼ n ¼ 1, 2, … (3.2.16)
w
Hence the difference between consecutive peaks is
np p np
R + R ¼ p cos np
w w w (3.2.17)
¼ p
the
When the excitation force is of the cosine type, pðt Þ ¼ p0 cos wt,
employed procedure yields the particular solution
p0 1
up ¼
cos wt (3.2.18)
k 1 b2
and the general solution for zero initial conditions, u ð0Þ ¼ u_ ð0Þ ¼ 0, is
obtained as
p0 1
u ðt Þ ¼ b cos wt Þ
ð cos wt (3.2.19)
k 1 b2
@CivilMethod
110 PART I Single-degree-of-freedom systems
p0 1
+ 2
1 b2 sin wt 2xb cos wt (3.2.25a)
k 1 b2 + ð2xb Þ2
|fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}
Steady-state response
p0 1
+ 2
1 b2 2xb sin wt + 1 b2 cos wt (3.2.25b)
k 1 b2 + ð2xbÞ2
|fflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl{zfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflfflffl}
Steady-state response
It is clear that the first term in Eq. (3.2.25a) decays rapidly with time because
of the exponential term exwt , so its contribution becomes negligible after a
short time, which, of course, depends on the damping ratio. Thus, the response
of the system is governed by the second term. For this reason, we say that the
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 111
first term expresses the transient response while the second term expresses the
steady-state response of the system. This is shown in Fig. 3.2.2.
Referring to Eq. (3.2.25a) we see that the steady-state response can be writ-
ten in the form
qÞ
u ðt Þ ¼ rsin ðwt (3.2.26)
where
p0 h
2 2 2
i1=2
1 2xb
r¼ 1 b + ð2xbÞ , q ¼ tan (3.2.27)
k 1 b2
Therefore, the steady-state response of the underdamped system subjected
to a harmonic loading is an undamped free vibration.
The DMF is
r h 2 i1=2
D ¼ max jRðt Þj ¼ ¼ 1 b2 + ð2xbÞ2 (3.2.28)
p0 =k
We observe that D ¼ D ðbÞ. Consequently the maximum value of D is
obtained when
dD
¼0 (3.2.29)
db
This condition gives
qffiffiffiffiffiffiffiffiffiffiffiffiffiffi
b ¼ 1 2x 2 (3.2.30)
@CivilMethod
112 PART I Single-degree-of-freedom systems
6
x=0
5
x=0.1
4
D (b)
3
x=0.2
2 x=0.3
ξ=0.7
x=0.5
1
ξ=1.0
0
0 0.5 1 1.5 2 2.5 3
Frequency ratio, b
FIG. 3.2.3 Variation of the amplitude of the response ratio Dðb, x Þ.
and
1
Dmax ¼ pffiffiffiffiffiffiffiffiffiffiffiffi (3.2.31)
2x 1 x2
Apparently, the maximum value of the dynamic factor does not occur for
b ¼ 1. Nevertheless, for a small value of x it is b 1, for example, for
x ¼ 0:05 Eq. (3.2.30) gives b ¼ 0:9975. Fig. 3.2.3 shows the variation of D ver-
sus the frequency ratio b for different values of the damping ratio x. Note that
if x ¼ 0, Dmax becomes infinite.
4
1/2x
0
R(t)
–2
–1/2x
–4
–6
0 5 10 15
t
FIG. 3.2.4 Response ratio of an underdamped system at resonance (w ¼ 5, x ¼ 0:1, b ¼ 1,
D ¼ max jRðt Þj ! 1=2x when t ! 1).
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 113
Let us study now the response of the system when b ¼ 1, which is the
conventional value for resonance. For this value of b it is w ¼ w and
Eq. (3.2.25a) becomes
p0 cos wt
u ðt Þ ¼ exwt ðAsin wD t + B cos wD t Þ (3.2.32)
k 2x
If we assume zero initial conditions, u ð0Þ ¼ u_ ð0Þ ¼ 0, Eq. (3.2.32) gives
" ! #
1 xwt x
Rðt Þ ¼ e pffiffiffiffiffiffiffiffiffiffiffiffi sin wD t + cos wD t cos wt (3.2.33)
2x 1 x2
m u€ + ku ¼ pðt Þ (3.3.1)
1 pffiffiffiffiffiffiffiffiffi
u€ + w2 u ¼ pðt Þ, w ¼ k=m (3.3.2)
m
There are several methods to obtain the solution of Eq. (3.3.2). A con-
venient, straightforward, and rather simple method for solving ordinary dif-
ferential equations with constant coefficients is the Laplace transform
method [1, 3]. This method is based on the Laplace transform, which for a func-
tion u ðt Þ, t 0 is commonly denoted by L and defined as
Z 1
L½u ðt Þ ¼ U ðs Þ ¼ u ðt Þest dt (3.3.3)
0
@CivilMethod
114 PART I Single-degree-of-freedom systems
u ðt Þ ¼ L1 ½u ðt Þ (3.3.4)
1
L½u€ + w2 L½u ¼ L½pðt Þ (3.3.6)
m
s 1 1 1
U ðs Þ ¼ u ð 0Þ + 2 u_ ð0Þ + P ðs Þ (3.3.8)
s2 +w 2 s +w 2 m s + w2
2
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 115
or
1 s 1 s 1 1 s
u ðt Þ ¼ u ð0ÞL + u_ ð0ÞL + L P ðs Þ
s 2 + w2 s 2 + w2 m s 2 + w2
(3.3.9)
From the table of the Laplace transforms we obtain
s
L1 2 ¼ cos wt (3.3.10a)
s + w2
1 1 sin wt
L 2 2
¼ (3.3.10b)
s +w w
Now we focus our attention on the last term on the right side of Eq. (3.3.9).
Its inverse Laplace transform can be obtained using the convolution theorem.
The convolution of two functions f ðt Þ and g ðt Þ denoted by f ðt Þ∗ g ðt Þ or
ðf ∗ g Þt is defined as
Z t
f ðt Þ∗ g ðt Þ ¼ f ðt τÞg ðτÞdτ
0
Z t (3.3.11)
¼ g ðt τÞf ðτÞdτ
0
@CivilMethod
116 PART I Single-degree-of-freedom systems
It obvious that the first two terms in Eq. (3.3.14) express the contribution of
the initial conditions to the motion of the system. The third term, which
expresses the contribution of the external loading, is known as the Duhamel
integral for the undamped system.
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 117
s 1 l1 t
L1 ¼ l1 e l2 el2 t
ðs l 1 Þð s l 2 Þ l1 l2 (3.3.22b)
xwt
¼e cos wD t
Z t
P ðs Þ 1
L1 ¼ pðτÞexwðtτÞ sin wD ðt τÞdτ (3.3.22c)
ðs l 1 Þ ð s l 2 Þ wD 0
It is obvious that the first term in Eq. (3.3.23) expresses the contribution of
the initial conditions to the motion of the system. The last term represents the
Duhamel integral for the underdamped motion. It is obvious that Eq. (3.3.23)
for x ¼ 0 yields Eq. (3.3.14) as anticipated. The Duhamel integral can also be
derived using the method described in Section 3.7.3.
@CivilMethod
118 PART I Single-degree-of-freedom systems
is the Heaviside step function and p0 the magnitude of the load. Fig. 3.4.1 shows
the step function load applied at t ¼ 0.
For u ð0Þ ¼ u_ ð0Þ ¼ 0 the displacement is obtained from the Duhamel integral
in Eq. (3.3.14)
Z t Z t
1 p0
u ðt Þ ¼ p0 sin ½wðt τÞdτ ¼ sin ½wðt τÞd ½wðt τÞ
mw 0 mw2 0
p0
¼ ½ cos wðt τÞt0
mw2
or
p0
u ðt Þ ¼ ð1 cos wt Þ (3.4.3)
k
Taking into account that ust ¼ p0 =k represents the static displacement, the
response ratio is
u ðt Þ
R ðt Þ ¼ ¼ ð1 cos wt Þ (3.4.4)
ust
and the DMF
D ¼ max jRðt Þj ¼ 2 (3.4.5)
Eq. (3.4.5) shows that the suddenly applied load produces a maximum dis-
placement that is twice as large than the displacement that the load p0 would
produce if it were applied statically (slowly). This is an elementary but impor-
tant result that illustrates the difference between static and dynamic loading of a
structure.
When the damping is taken into account, the displacement for zero initial
conditions, u ð0Þ ¼ u_ ð0Þ ¼ 0, is obtained from Eq. (3.3.23)
Z t
p0
u ðt Þ ¼ exwðtτÞ sin ½wD ðt τÞdτ (3.4.6)
mwD 0
The evaluation of the Duhamel integral is more complicated. Nevertheless,
using MAPLE we obtain
" ! #
p0 x
u ðt Þ ¼ 1 cos wD t + pffiffiffiffiffiffiffiffiffiffiffiffi sin wD t exwt (3.4.7)
k 1 x2
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 119
2.5
x=0
2
x=0.05
x=0.15
1.5
1
R(t)
0.5
0.
–0.5
0 0.5 1 1.5 2 2.5 3 3.5 4
t/TD (TD=2p /w D)
FIG. 3.4.2 Response ratio of a system under step load.
Fig. 3.4.2 shows the plot of Eq. (3.4.8) for different values of the
damping ratio.
The extreme values of Rðt Þ occur when
" #
2
dR xwt ðwx Þ
¼e + wD sin wD t ¼ 0 (3.4.9)
dt wD
which gives
np
tn ¼ n ¼ 0,1, 2, … (3.4.10)
wD
@CivilMethod
120 PART I Single-degree-of-freedom systems
that is, the maximum displacement occurs at the first peak, where t1 ¼ p=wD .
This is shown in Fig. 3.4.2. Eq. (3.4.12) for x ¼ 0 gives D ¼ 2 as anticipated.
Fig. 3.4.4 shows the response ratio Rðt Þ ¼ u ðt Þ=ust , ust ¼ p0 =k, of the
undamped system to a ramp function load. We see that it oscillates about the
line p0 t=t1 .
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 121
30
25
20
R(t)
15
10
0
0 0.5 1 1.5 2 2.5 3
t/T
FIG. 3.4.4 Response ratio of a system under ramp function load (t1 ¼ T =10).
3.4.3 Response to step function load with finite rise time. Static load
The step function load with finite rise time is a constant load that, however, is
not applied suddenly, but rises linearly up to a value p0 within a time t1 , the rise
time, and thereafter remains constant, as shown in Fig. 3.4.5. Mathematically it
is defined by
(p
0
t t t1
pðt Þ ¼ t1 (3.4.16)
p0 t > t1
@CivilMethod
122 PART I Single-degree-of-freedom systems
u_ I ðt1 Þ p0
uII ðt Þ ¼ sin wt+ uI ðt1 Þcos wt+ ð1 cos wtÞ (3.4.18)
w k
p0
uII ðt Þ ¼ (3.4.20)
k
or
RII ðt Þ ¼ 1 (3.4.21)
which means that the motion in the constant load phase is not oscillatory but the
displacement is constant and equal to the static displacement.
Fig. 3.4.6 presents the response ratio Rðt Þ for different values of t1 .
We observe that for smaller values of t1 =T , the response is similar to that
of a step function load while for larger values, the response is similar to that of
a static load. Therefore, the loads in real structures should not be applied
suddenly but slowly rising, that is, in time much larger than the natural
period of the structure to avoid dynamic magnification effects. Due to this
property, the step function load with finite rise time is also called static
loading.
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 123
FIG. 3.4.6 Response ratio of a system under a step function load with finite rise time.
Solution
The mass of the structure, its stiffness in the x direction, and the corresponding
natural frequency were evaluated in Example 2.2.1
m ¼ 356:78kN m1 s2 , k ¼ Kx ¼ 7257:6kN=m, w ¼ w1 ¼ 4:51
@CivilMethod
124 PART I Single-degree-of-freedom systems
1.5
D = 1.265
0.5
R(t)
tmax = 0.563
–0.5
–1
–1.5
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
t
FIG. E3.2 Response ratio in Example 3.4.1.
Fig. E3.2 shows the plot of the response ratio Rðt Þ ¼ u ðt Þ=ðP=k Þ. The DMF
was found at D ¼ max jRðt Þj ¼ 1:1265. It was determined as the maximum
value of the array used to plot Rðt Þ. Thus, we have
umax ¼ Du st ¼ 1:1265 0:02067 ¼ 0:02328m (5)
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 125
The extreme values of the shear forces and bending moments are:
Columns 30 30
max Qx ¼ kx umax ¼ 1360:8 0:02328 ¼ 31:68kN
h
max Mx ¼ max Qx ¼ 31:68 5=2 ¼ 79:20kNm
2
Columns 30 20
max Qx ¼ kx umax ¼ 907:2 0:02328 ¼ 21:12kN
h
max Mx ¼ max Qx ¼ 21:12 5=2 ¼ 52:80kNm
2
t
@CivilMethod
FIG. 3.5.1 Time variation of an impulsive load.
126 PART I Single-degree-of-freedom systems
The response of the system to an impulsive load can be analyzed using the
methods presented in previous sections for the solution of the differential equa-
tion of motion under an arbitrary loading, that is, either by solving directly the
differential equation or by evaluating Duhamel’s integral. Another method to
obtain the response is to express the pulse as the superposition of two or more
simpler pulses for which the response solution is available or simple to deter-
mine (Fig. 3.5.2). Nevertheless, the analytical methods, especially for arbitrary
impulse loads, have lost their importance because of the development of effi-
cient numerical methods. For this reason, only two simple impulsive loads
are considered, the rectangular pulse load and the triangular pulse load.
p(t )
p(t ) p 0H (t ) if t t1
p0 p(t ) 0 if t t1
O t
t1
p(t )
p0 p(t ) p0 sin t if t t1
p(t ) 0 if t t1
O t
t1
FIG. 3.5.2 Impulsive loads represented by simple functions.
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 127
@CivilMethod
128 PART I Single-degree-of-freedom systems
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
u_ I ðt1 Þ 2
umax ¼ + ½ u I ðt1 Þ 2 (3.5.6)
w
Using Eqs. (3.5.4a), (3.5.4b) the previous relation becomes
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
p0 2p 2p 2p
umax ¼ 1 2cos t1 + cos 2 t1 + sin 2 t1
k T T T
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
ffi
p0 2p
¼ 2 1 cos t1
k T
or
2p0 t1
umax ¼ sin p (3.5.7)
k T
and the DMF D will be given by
umax t1 T
D¼ ¼ 2sin p , t1 (3.5.8)
p0 =k T 2
Eq. (3.5.8) shows that the maximum response depends only on the ratio
t1 =T . The plot of the function D ¼ D ðt1 =T Þ shown in Fig. 3.5.4 is referred
to as the displacement response spectrum or simply the response spectrum of
the impulsive load (see also Chapter 6). It is evident that the response spectrum
of a pulse load serves to determine the maximum response of the system under
this load without solving the differential equation of motion for the particular
pulse load.
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 129
u_ I ðt1 Þ
uII ðt Þ ¼ sin wt+ uI ðt1 Þcos wt, t¼ t t1 0 (3.5.11)
w
where
p0 sin wt1
u I ðt1 Þ ¼ coswt1 (3.5.12a)
k t1 w
p0 w cos wt1 1
u_ I ðt1 Þ ¼ + sin wt1 (3.5.12b)
k wt1 wt1
cos wt 1
u_ I ðt Þ ¼ w sin wt + ¼0 (3.5.13)
t1 t1
@CivilMethod
130 PART I Single-degree-of-freedom systems
2
3
t1
sin 2p
p0 6 T t1 7
uI ðt1 Þ ¼ 6 4 cos 2p 7 (3.5.16a)
k t1 T 5
2p
T
2
3
t1
cos 2p
p0 w 6
6 T t1 1 7 7
u_ I ðt1 Þ ¼ + sin 2p (3.5.16b)
k 4 t1 T t1 5
2p 2p
T T
The curve RII in Fig. 3.5.6 represents the function max uII =ust versus t1 =T .
From this figure, we conclude that for t1 =T 0:4 the maximum response D of
the system to the triangular pulse load occurs in Phase II while for t1 =T > 0:4 it
occurs in Phase I.
R = maxu /ust
D(0.4)=1.0513
D(t1/T)
R = maxu /ust
t1/T
FIG. 3.5.6 DMF for the triangular pulse load.
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 131
p0 sin wt
uI ðt Þ ¼ t , 0 t t1 (3.5.17)
k t1 w
u_ I ðt1 Þ
uII ðt Þ ¼ sin wet + uI ðt1 Þ cos wet
w
(3.5.18)
p0 e sin wet et
+ 1 cos wt + , 0 et ¼ t t1 < t2
k t2 w t2
u_ II ðt1 + t2 Þ
uIII ðt Þ ¼ sin wt+ uII ðt1 + t2 Þcos wt, t¼ t ðt1 + t2 Þ 0
w
(3.5.19)
@CivilMethod
132 PART I Single-degree-of-freedom systems
1 sin wt
Rðt Þ ¼ t , t t1
t1 w
1
Rðt Þ ¼ 1 + ½ sinwðt t1 Þ sin wt , t t1
wt1
sin wt t
Rðt Þ ¼ 1 cos wt + , t t1
wt1 t1
sinwt1 sinwðt t1 Þ
Rðt Þ ¼ cos wt, t t1
wt1
1 sin wt t1
Rðt Þ ¼ t , 0t
t1 w 2
2 1 t1 t1
Rðt Þ ¼ t1 t + 2sin w t sin wt , t t1
t1 w 2 2
2 t1
Rðt Þ ¼ sinwðt t1 Þ + 2sinw t sinwt , t t1
wt1 2
b sinwt
sin wt
Rðt Þ ¼
1 b2
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 133
y
k1 k2 k1
M (t ) x
O 5m
k1 k2 k1
10 m
(a) (b)
(c)
FIG. E3.3 One-story building under impulsive load.
Solution
Because the structure is symmetric with respect to both axes x and y, the only
possible motion of the slab due to the moment M ðt Þ is the rotation fðt Þ about its
center O. Hence the equation of motion of the slab is
where IO is the moment of inertia of the mass of the plate with respect to O and
Kf the torsional stiffness of the structure. These quantities have been computed
in Example 2.2.1. Thus we have
rffiffiffiffiffiffi
Kf
w¼ ¼ 7:2738 (4)
IO
2p
T¼ ¼ 0:869s (5)
w
t1 ¼ T ¼ 0:869s (6)
@CivilMethod
134 PART I Single-degree-of-freedom systems
f_ I ðt1 Þ
fII ðt Þ ¼ sin wðt t1 Þ + fI ðt1 Þcos wðt t1 Þ (13)
w
Eq. (12) for t ¼ t1 gives
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 135
–3
4 ×10
fmax = 3.594e-3
3
2
f(t)
0
tmax = 0.462
–1 t1
–2
0 0.5 1 1.5 2 2.5 3 3.5 4
t
FIG. E3.4 Response of the structure in Example 3.5.1.
The stress resultants of the columns are computed using Eqs. (4a), 4(b), (6)
of Example 2.2.1. Thus, for the upper right column we have
kx1 ¼ ky1 ¼ 1360:8kN=m, kf1 ¼ 1:9985 103 kNm
x1 ¼ 4:85m, y1 ¼ 2:35m
max Qx1 ¼ kx1 y1 fmax ¼ 11:49kN (16a)
max Qy1 ¼ ky1 x1 fmax ¼ 23:72kN (16b)
h
max Mx1 ¼ max Qx1 ¼ 28:73kNm (16c)
2
h
max My1 ¼ max Qy1 ¼ 59:30kNm (16d)
2
max Mf1 ¼ kf1 fmax ¼ 7:183kNm (16e)
@CivilMethod
136 PART I Single-degree-of-freedom systems
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 137
in which
1 xbi
Ai0 ¼ ai 2
k w
bi
Ai1 ¼ , Ai2 ¼ u ðti1 Þ Ai0 (3.5.25)
k
1
A3 ¼
i
u_ ðti1 Þ + xwAi2 Ai1
wD
Differentiating Eq. (3.5.24) with respect to time gives the velocity
u_ i ðt Þ ¼ Ai1 + exwt xwAi2 + wD Ai3 cos wD t xwAi3 + wD Ai2 sin wD t
(3.5.26)
The presented method is exact. However, numerical methods are more con-
venient to compute the response to a piecewise linear loading (see Chapter 4).
(a)
(b)
(c)
FIG. 3.6.1 Examples of periodic loads.
@CivilMethod
138 PART I Single-degree-of-freedom systems
Z T =2
1sin mw0 tdt ¼ 0 for each m (3.6.5b)
T =2
Z T =2
0 if m 6¼ n
cos mw0 t cos nw0 tdt ¼ (3.6.5c)
T =2 T=2 if m ¼ n
Z T =2
cos mw0 t sin nw0 tdt ¼ 0 for each m, n (3.6.5d)
T =2
Z T =2
0 if m 6¼ n
sin mw0 t sin nw0 tdt ¼ (3.6.5e)
T =2 T =2 if m ¼ n
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 139
Z T =2
1
a0 ¼ pðt Þdt (3.6.6a)
T T =2
Z T =2
2
an ¼ pðt Þcos nw0 tdt (3.6.6b)
T T =2
The series (3.6.2) represents the function pðt Þ, that is, converges to pðt Þ for
n ! 1, provided that it satisfies the following conditions, known as Dirichlet
conditions:
(a) The function pðt Þ has a finite number of discontinuities in one period.
(b) The function pðt Þ has a finite number of maxima and minima in one period.
(c) The function pðt Þ is absolutely integrable over a period, that is,
Z T =2
jpðt Þjdt ¼ k < 1 (3.6.7)
T =2
We shall say that that the function pðt Þ is piecewise continuous in the finite
interval ½T =2, T =2, if it satisfies conditions (a) and (b). At the points of dis-
continuity, for example, point t1 in Fig. 3.6.1c, the Fourier series converges to
the mean value
1
p t1 + p t1+ (3.6.8)
2
where p t1 and p t1+ are the left and right limits of pðt Þ at t1 .
In practice, the periodic function pðt Þ is approximated by a finite number of
terms of the Fourier series, that is, by a finite Fourier series.
Let
X
k
Sk ð t Þ ¼ a 0 + ðan cos nw0 t + bn sin nw0 t Þ (3.6.9)
n¼1
be the sum of the first k + 1 terms of the Fourier series, which will represent the
function pðt Þ in the interval ½T=2, T =2. Then we will have
X
k
pðt Þ ¼ a0 + ðan cosnw0 t + bn sin nw0 t Þ + ek ðt Þ (3.6.10)
n¼1
@CivilMethod
140 PART I Single-degree-of-freedom systems
where
ek ðt Þ ¼ pðt Þ Sk ðt Þ (3.6.11)
is the error between pðt Þ and its approximation. The mean square error Ek is
given by
Z
1 T =2
Ek ¼ ½ek ðt Þ2 dt
T T =2
Z " #2 (3.6.12)
1 T =2 Xk
¼ pðt Þ a0 ðan cos nw0 t + bn sin nw0 t Þ dt
T T =2 n¼1
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 141
(a) (b)
FIG. 3.6.2 Functions pðt Þ replacing pðt Þ in the interval ½t1 e, t1 + e.
@CivilMethod
142 PART I Single-degree-of-freedom systems
Z T
2 cos np
bn ¼ pðt Þ sin nw0 t ¼ p0 (2c)
T 0 np
Hence
p0 p0 X
N
cos np 1 p0 X
N
cos np
pðt Þ ¼ + 2 cos nw 0 t + sin nw0 t (3)
4 p n¼1 n2 p n¼1 n
Fig. E3.5 shows the graphical representation of the finite Fourier series in
Eq. (3) for various values of N with p0 ¼ 1, T ¼ 1. We observe that the conver-
gence is very slow at point t ¼ T =2 due to the Gibbs phenomenon.
N = 20 N = 100
1.2 1.2
1 1
0.8 0.8
0.6 0.6
p(t)
p(t)
0.4 0.4
0.2 0.2
0 0
–0.2 –0.2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
t t
N = 1000 N = 2000
1.2
1 1
0.8 0.8
0.6 0.6
p(t)
p(t)
0.4 0.4
0.2 0.2
0 0
–0.2 –0.2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
t t
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 143
Solution
The function pðt Þ within a period T ¼ 3t1 can be represented using the delta
function.a Thus we obtain
pðt Þ ¼ P ½d ðt T =3Þ d ðt 2T =3Þ (1)
Eqs. (3.6.6a)–(3.6.6c) give
Z
1 T
a0 ¼ pðt Þdt ¼ 0 (2a)
T 0
Z
2 T 2P 2np 4np
an ¼ pðt Þcos nw0 tdt ¼ cos cos (2b)
T 0 T 3 3
Z
2 T 2P 2np 4np
bn ¼ pðt Þ sin nw0 t ¼ sin sin (2c)
T 0 T 3 3
Inserting these values of the coefficients in Eq. (3.6.2) we obtain
(
1
2P X 2np 4np
pðt Þ ¼ cos cos cos nw0 t
T n¼1 3 3
)
2np 4np
+ sin sin sin nw0 t (3)
3 3
a. The delta function d ðt t0 Þ and its properties are discussed in Section 3.7.1.
@CivilMethod
144 PART I Single-degree-of-freedom systems
For the sine term bn sin nw0 t the steady-state response results from
Eq. (3.2.25a) by setting p0 ¼ bn , w ¼ nw0 , and b ¼ b n ð¼ nw0 =wÞ
bn 1
unsin ðt Þ ¼ 2 1 b2n sin nw0 t 2xbn cos nw0 t
k 1b 2
+ ð2xbn Þ 2
n
(3.6.16b)
while for the cosine term an cos nw0 t the steady-state response is obtained from
Eq. (3.2.25b) by setting p0 ¼ an , w ¼ nw0 , and b ¼ b n ð¼ nw0 =wÞ
an 1
uncos ðt Þ ¼ 2xbn sin nw0 t + 1 b 2n cos nw0 t
k 1 b 2 2 + ð2xb Þ2
n n
(3.6.16c)
The steady-state response of the damped system to periodic loading results
as the superposition of responses to individual terms of the Fourier series
X
n X
n
u ðt Þ ¼ u0 ðt Þ + unsin ðt Þ + uncos ðt Þ (3.6.17)
n¼1 n¼1
(3.6.19)
Theoretically, there is no transient response when x ¼ 0. However, for
small values of damping, which is the usual case in our structures, the
steady-state response can be obtained from Eqs. (3.6.16b), (3.6.16c) for
x 0. This yields
bn 1
unsin ðt Þ sin nw0 t (3.6.20a)
k 1 b2n
an 1
uncos ðt Þ cos nw0 t (3.6.20b)
k 1 b 2n
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 145
(a) (b)
FIG. E3.7 Response to periodic loading in Example 3.6.3.
@CivilMethod
146 PART I Single-degree-of-freedom systems
FIG. 3.7.1 Distribution of a concentrated impulsive force over the interval ½e, +e.
namely, the total force is equal to unity. If we skip the problem of determining
analytically the function pðt Þ, we may assume a priori a prescribed shape for
this function, for example,
k=2, jt j < 1=k
pk ðt Þ ¼ (3.7.2a)
0, jt j 1=k
or
k
pk ðt Þ ¼ (3.7.2b)
pð1 + k 2 t 2 Þ
where k is a positive number.
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 147
Fig. 3.7.2a and b show both functions pk defined in Eqs. (3.7.2a), (3.7.2b),
respectively. Moreover, they satisfy Eq. (3.7.1), which means that they are
equivalent to pðt Þ and can represent the actual force P.
(a) (b)
FIG. 3.7.2 Functions pk sufficiently concentrated for large values of k.
or by the relation
Z +1
dðt t1 Þf ðt Þ dt ¼ f ðt1 Þ (3.7.5)
1
@CivilMethod
148 PART I Single-degree-of-freedom systems
and
Z +1 Z e
d ðt Þdt ¼ d ðt Þdt ¼ 1 (3.7.6b)
1 e
or for e ! 0
mDu_ ¼ 1 (3.7.14)
where
Du_ ¼ u_ t1+ u_ t1 (3.7.15)
Hence the impulsive load produces an abrupt change (discontinuity) of the
t1 . If the system is at rest before the action of the impulsive load,
velocity at time
then it is u_ t1 ¼ 0 and Eq. (3.7.14) gives
1
u_ ðt1 Þ ¼ (3.7.16)
m
where t1 designates t1+ . Therefore, Eq. (3.7.16) presents the initial velocity
given to the system by the impulsive load. However, the displacement remains
continuous, which means that
Du ¼ u t1+ u t1 ¼ 0 (3.7.17)
If the elastic and damping forces are taken into account, then Eq. (3.7.12b)
is written
d u_ du
m ¼ pðt Þ c ku (3.7.18)
dt dt
which after integration over the interval ½e, +e gives
Z +e Z +e
+e +e
m ½u_ e ¼ pðt Þdt c½u e k u ðt Þdt (3.7.19)
e e
Applying the mean value theorem of integral calculus to the integral of elas-
tic force, we write Eq. (3.7.19) as
Z +e
+e +e
m ½u_ e ¼ pðt Þdt c½u e ku ðt ∗ Þ2e, e < t ∗ < e (3.7.20)
e
@CivilMethod
150 PART I Single-degree-of-freedom systems
h(t-t)
t =1
t
FIG. 3.7.3 Response to unit impulse.
1
u€ + 2xwu_ + w2 u ¼ d ðt τ Þ (3.7.23)
m
The solution of the above equation can be obtained from Eq. (3.3.23) for
zero initial conditions and pðt Þ ¼ dðt t1 Þ. Thus, we have
Z t
1
u ðt Þ ¼ dðτ t1 ÞexwðtτÞ sin ½wD ðt τÞdτ (3.7.24)
mwD 0
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 151
WI , II ¼ WII , I (3.7.28)
In dynamics of linear systems, the reciprocity is expressed by the dynamic
reciprocal theorem [10], known also as the dynamic Betti-Rayleigh theorem.
This theorem for SDOF systems reads
If two loadings pI ðt Þ and pII ðt Þ act separately on a linear dynamic system and
produce the responses uI ðt Þ and uII ðt Þ then the convolution CI , II of the loading
pI ðt Þ with the response uII ðt Þ is equal to the convolution CII , I of the loading pII ðt Þ
with the response uI ðt Þ
@CivilMethod
152 PART I Single-degree-of-freedom systems
CI , II ¼ CII , I (3.7.29)
or using definition (3.3.11) for the convolution, we may write
pI ðt Þ∗ uII ðt Þ ¼ pII ðt Þ∗ uI ðt Þ (3.7.30a)
or in integral form
Z t Z t
pI ðτÞuII ðt τÞdτ ¼ pII ðτÞuI ðt τÞdτ (3.7.30b)
0 0
3.8 Problems
Problem P3.1 A machine carrying a mass m0 is placed on the slab of the build-
ing of Fig. P3.1. The mass rotates eccentrically about the point ð2:5, 2:0Þ at a
distance s ¼ 1:0m with a frequency f ¼ 4Hz. The columns are massless, inex-
tensible, and fixed on the base while the slab is assumed uniform and rigid.
Determine the time history of the shear forces Qx , Qy , the bending moments
Mx , My , and the torsion moment Mw at the top of the columns. The material
constants are E ¼ 1:2 107 kN=m2 , n ¼ 0:2. The total load of the plate (dead
plus live) is q ¼ 12kN=m2 ; m0 ¼ m=5, where m is the total mass of the slab.
The acceleration of gravity is g ¼ 9:81m=s2 .
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 153
(a) (b)
FIG. P3.4 One-story building in problem P3.4
@CivilMethod
154 PART I Single-degree-of-freedom systems
the x direction. Plot the response ratio Rðt Þ of the structure and the stress resul-
tants at the base of the columns 1 and 3 using the following data: Height of col-
umns a ¼ 4m, cross-section of columns a=10 a=20, side length of the
triangular slab a, p0 ¼ 10kN, E ¼ 2:1 107 kN=m2 , load of the plate, including
the dead load, q ¼ 20kN=m2 , t1 ¼ p=2w, x ¼ 0:05, w ¼ 8s1 , acceleration of
gravity g ¼ 9:81m=s . 2
Problem P3.5 The structure of Fig. P3.5a consists of the rigid girder BC and
the two flexible columns AB and CD having a cross-sectional moment of inertia
I and a modulus of elasticity E. The cables FB and GC have cross-sectional area
A, cannot withstand compression, and are assumed massless. The structure is sub-
jected to the impulsive loads pðt Þ shown in Fig. P3.5b and c. Study the response of
the structure and determine the maximum error when the impulsive loads are
substituted by equivalent concentrated forces. Plot the function D ðt1 =T Þ for the
two load cases. For which value of the ratio t1 =T is the maximum error less than
2%? Data: a ¼ 1:5m, I ¼ 33,740cm4 (IPE450), E ¼ 2:1 108 kN=m2 ,
A ¼ 3cm2 , t1 ¼ 0:1, m ¼ 1:0kNm1 s2 =m, p10 ¼ 10kN. The value p20 is deter-
mined so that both loads have the same impulse.
(a)
(b) (c)
FIG. P3.5 Structure in problem P3.5
@CivilMethod
Single-degree-of-freedom systems: Forced vibrations Chapter 3 155
Problem P3.6 Consider the structure of Fig. P3.6. The rigid column AC of circu-
lar cross-section and mass per unit length m is supported by the three elastic cables
of cross-sectional area A and modulus of elasticity E. The support on the ground is
a spherical hinge. The cables have been prestressed so that they can withstand com-
pression. Three advertising panels are massless fixed at the top of the columns, as
shown in the figure. The structure is subjected to the wind blast load of Fig. P3.6c in
the y direction (see Fig. P3.6b). Determine the minimum prestress force of the
cables. The cables and the panels have negligible mass. Use the data a ¼ 5:0m,
p0 ¼ 4kN=m2 , m ¼ 0:5kNm1 s2 =m, E ¼ 2:1 108 kN=m2 , A ¼ 4cm2 ,
t1 ¼ 0:1s.
(a) (b)
(c)
FIG. P3.6 Structure in problem P3.6
@CivilMethod
156 PART I Single-degree-of-freedom systems
10
p
0
4
2
t1 t2
0
p1
–2
0 0.5 1 1.5 2
FIG. P3.7 Blast pressure in problem P3.7
Problem P3.8 Show that an impulsive load generated by the load pðt Þ
and acting over the interval ½t1 , t2 can be represented as pI ¼ ½H ðt t1 Þ
H ðt t2 Þpðt Þ, where H ðt ti Þ is the Heaviside step function. Write a
MATLAB program that constructs the impulsive load.
Problem P3.9 Determine the dynamic response of a SDOF system subjected to
the sine periodic loading of Fig. P3.9. Assume: T ¼ p=w, w ¼ 1:1w,
x ¼ 0, and
x ¼ 0:1.
@CivilMethod
Chapter 4
4.1 Introduction
The previous analysis shows that an analytical solution of the equation of
motion for an single-degree-of-freedom (SDOF) system is possible only if
the external force is described by a simple function. If the excitation force varies
arbitrarily with time or is given by a set of its values, an analytical solution is out
of the question. However, such problems can be tackled numerically by time
step integration methods for differential equations. The literature about these
methods is vast. Extensive chapters and whole books cover this subject. They
present the mathematical development of these methods, their computer imple-
mentation, and their accuracy, convergence, and stability. Several computer
packages include ready-to-use subroutines for the solution of the differential
equation of motion.
Some of these methods have been specially developed for the study of the
dynamic response of systems. A survey of these methods is given in [1, 2]. The
central difference method (CDM), Houbolt’s method, Wilson’s q-Method, and
Newmark’s method are the most well known among them [3, 4]. Nevertheless,
with the increase of cheap computer power, some of them have lost their impor-
tance while others have taken dominating places in the computational arena.
Further, adding Eqs. (4.2.1), (4.2.2) and neglecting the terms of order higher
than three, we obtain the following expression to approximate the second deriv-
ative of u ðt Þ at time t
u ðt + Dt Þ 2u ðt Þ + u ðt Dt Þ
u€ðt Þ (4.2.5)
Dt 2
Substitution of the derivatives u_ ðt Þ and u€ðt Þ into Eq. (4.1.1) with their
approximations (4.2.4) and (4.2.5) gives
u ðt + Dt Þ 2u ðt Þ + u ðt Dt Þ u ðt + Dt Þ u ðt Dt Þ
m +c + ku ðt Þ ¼ pðt Þ
Dt 2 2Dt
which is solved for u ðt + Dt Þ to yield
m
c 2m m c
+ u ð t + Dt Þ ¼ p ð t Þ k u ðt Þ u ðt Dt Þ
Dt 2 2Dt Dt 2 Dt 2 2Dt
(4.2.6)
K^ u ðt + Dt Þ ¼ P^ (4.2.7)
where
m c
K^ ¼ 2 + (4.2.8)
Dt 2Dt
2m m c
P^ ¼ pðt Þ k 2 u ðt Þ u ðt Dt Þ (4.2.9)
Dt Dt 2 2Dt
The quantities K^ and P^ are referred to as the effective stiffness and the
effective load, respectively. Obviously, Eq. (4.2.7) allows the evaluation of
the displacement at instant t + Dt, if the displacements at the two preceding
instants t and t Dt are known. Because u ð0Þ is known from the initial
@CivilMethod
162 PART I Single-degree-of-freedom systems
conditions, the procedure starts at t ¼ Dt. Obviously, this requires the value of
u ðDt Þ, which is unknown in the first instance, but it can be determined from
Eq. (4.2.2) for t ¼ 0. Thus neglecting terms of order higher than two, we have
1
u ðDt Þ u ð0Þ Dt u_ ð0Þ + Dt 2 u€ð0Þ (4.2.10)
2
In the above equation, the quantities u ð0Þ and u_ ð0Þ are known from the ini-
tial conditions while u€ð0Þ can be computed from the equation of motion,
Eq. (4.1.1), for t ¼ 0. Thus, we obtain
u€ð0Þ ¼ ½pð0Þ cu_ ð0Þ ku ð0Þ=m (4.2.11)
The stability of the CDM requires that time step Dt is less than a certain
critical value, that is, qffiffiffiffiffiffiffiffiffiffiffiffi
Dt Dtcr ¼ T 1 x2 =p (4.2.12)
where T is the period of the system (see Eq. 4.5.30). Otherwise, the procedure
“blows up” with time and the solution makes no sense. This is discussed in
Section 4.5.
Because T is usually a small number, Dt should be small, which implies
that a large number of time steps are required to solve the equation of motion.
This has been a major drawback of the method, especially in older times when
the computer capabilities in terms of memory and speed were restricted. This
fact has led researchers to develop integration methods in which the size of
the time step is not restricted by a critical value. Table 4.2.1 presents the
@CivilMethod
Numerical integration of the equation of motion Chapter 4 163
Solution
The solution is obtained using the program centr_diff_lin.m with Dt ¼ 0:01.
Fig. E4.1 gives the response of the SDOF system. Moreover, Fig. E4.2 shows
Displacement Velocity
0.2 0.5
0.1
u,t(t)
u(t)
0
0
–0.1 –0.5
0 5 10 0 5 10
t t
Acceleration Response ratio
4 2
2
1
u,tt (t)
R(t)
0
0
–2
–4 –1
0 5 10 0 5 10
t t
FIG. E4.1 Response of the SDOF system in Example 4.2.1.
@CivilMethod
164 PART I Single-degree-of-freedom systems
the displacement u ðt Þ as compared with the exact one together with the error
u ðt Þ uex ðt Þ. The exact solution was obtained by analytical evaluation of
Duhamel’s integral giving
4po ð1 2xwÞ
u ðt Þ ¼ e 10:5t
+ cos w D t + sin w D t e 1xwt
m ð1 4xw + 4w2 Þ 2wD
(1)
x 10 -4
0.2 1
u(t) computed
u(t) exact 0.8 u-uex
0.15 0.6
0.4
0.1 0.2
u(t)
0
0.05 –0.2
–0.4
0 –0.6
–0.8
-0.05 –1
0 2 4 6 8 10 0 2 4 6 8 10
t
FIG. E4.2 Computed solution and error in Example 4.2.1.
@CivilMethod
Numerical integration of the equation of motion Chapter 4 165
FIG. 4.3.1 Variation of the acceleration, velocity, and displacement in the average acceleration
method.
τ2
u ðt + τÞ ¼ τu_ ðt Þ + ½u€ðt Þ + u€ðt + Dt Þ + C2
4
@CivilMethod
166 PART I Single-degree-of-freedom systems
Setting
Du ¼ u ðt + Dt Þ u ðt Þ (4.3.6a)
Du_ ¼ u_ ðt + Dt Þ u_ ðt Þ (4.3.6b)
Du€ ¼ u€ðt + Dt Þ u€ðt Þ (4.3.6c)
Eqs. (4.3.4), (4.3.5) are written as
Dt
Du_ ¼ ½2u€ðt Þ + Du€ (4.3.7)
2
Dt 2
Du ¼ Dt u_ ðt Þ + ½2u€ðt Þ + Du€ (4.3.8)
4
€ we obtain
Solving Eq. (4.3.8) for Du,
4
Du€ ¼ ½Du Dt u_ ðt Þ 2u€ðt Þ (4.3.9)
Dt 2
2
Du_ ¼ Du 2u_ ðt Þ (4.3.10)
Dt
@CivilMethod
Numerical integration of the equation of motion Chapter 4 167
@CivilMethod
168 PART I Single-degree-of-freedom systems
As will be shown, contrary to the CDM, the stability of the AAM does not
demand any restriction on the size of the time step Dt. The time step, however,
is influenced by the accuracy of the method and its capability to describe an
oscillatory motion. Therefore, it must be small enough. The selection of Dt
equal to 1/10 of the period of the system or of the period of the excitation force
produces accurate results.
Adhering to the steps of Table 4.3.1, a computer program called av_acc_lin.m
has been written in MATLAB for the numerical integration of the equation of
motion. The program is available on this book’s companion website. It computes
the displacement u ðt Þ, the velocity u_ ðt Þ, the acceleration u€ðt Þ, and the response
ratio Rðt Þ ¼ u ðt Þ=ðpmax =k Þ and gives their graphical representation. Moreover, it
computes the dynamic magnification factor D ¼ max jRðt Þj and the time tmax it
occurs. The user of the program is responsible for providing the function of the
excitation force.
Solution
The solution was evaluated using the program av_acc_lin.m with Dt ¼ 0:01.
Fig. E4.3 gives the graphical representation of the displacement together with
the error u ðt Þ uex ðt Þ. The computed error by the AAM is almost double the
error of the CDM. Moreover, Fig. E4.4 shows the response of the system under
the harmonic load p ¼ 2 sin 5t. Obviously, this excitation produces resonance
(w ¼ w ¼ 5).
x 10−4
0.2 2
u(t) computed
u(t) exact 1.5 u-uex
0.15
1
0.5
0.1
u(t)
0.05
–0.5
–1
0
–1.5
−0.05 –2
0 2 4 6 8 10 0 2 4 6 8 10
t t
FIG. E4.3 Computed solution and error in Example 4.3.1.
@CivilMethod
Numerical integration of the equation of motion Chapter 4 169
Displacement Velocity
0.01 0.04
0.005 0.02
u,t(t)
u(t)
0 0
–0.005 –0.02
–0.01 –0.04
0 5 10 15 0 5 10 15
t t
Acceleration Response ratio
0.2 10
0.1 5
u,tt (t)
R(t)
0 0
–0.1 –5
–0.2 –10
0 5 10 15 0 5 10 15
t t
FIG. E4.4 Response of the system in Example 4.3.1 under harmonic load.
@CivilMethod
170 PART I Single-degree-of-freedom systems
it, we have
u€ ¼ q ðt Þ (4.4.3)
where q ðt Þ is a fictitious source, unknown in the first instance. Eq. (4.4.3) is the
analog equation of Eq. (4.4.1). It indicates that the solution of Eq. (4.4.1) can
be obtained by solving Eq. (4.4.3) with the initial conditions (4.4.2), if the q ðt Þ
is first established. This is achieved as follows.
Taking the Laplace transform of Eq. (4.4.3) we obtain
1 1 1
U ðs Þ ¼ u ð0Þ + 2 u_ ð0Þ + 2 Q ðs Þ (4.4.4)
s s s
where U ðs Þ,Q ðs Þ are the Laplace transforms of u ðt Þ,q ðt Þ, respectively. The
inverse Laplace transform of Eq. (4.4.4) gives the solution in integral from
Z t
u ðt Þ ¼ u ð0Þ + u_ ð0Þt + q ðτÞðt τÞdτ (4.4.5)
0
Thus, the IVP of Eqs. (4.4.1), (4.4.2) is transformed into the equivalent
Volterra integral equation for q ðt Þ.
Eq. (4.4.5) is solved numerically within a time interval ½0, T . The interval
½0, T is divided into N equal intervals, Dt ¼ h, h ¼ T =N , in which q ðt Þ is
assumed to vary according to a certain law, for example, constant, linear,
etc. In this analysis, q ðt Þ is assumed to be constant and equal to the mean value
in the interval h (Fig. 4.4.1). That is
qr1 + qr
qrm ¼ (4.4.6)
2
FIG. 4.4.1 Discretization of the interval ½0, T into N equal intervals h ¼ T=N .
(4.4.7)
@CivilMethod
Numerical integration of the equation of motion Chapter 4 171
where
h2
c1 ¼ (4.4.9)
2
The velocity is obtained by direct differentiation of Eq. (4.4.5) using
Leibnitz’ rule for integrals [12]. Thus, we have
Z t
u_ ðt Þ ¼ u_ ð0Þ + q ðτÞdτ (4.4.10)
0
Using the same discretization for the interval ½0, T to approximate the inte-
gral in Eq. (4.4.10), we have
X
n1
u_ n ¼ u_ 0 + c2 qrm + c2 qnm
r¼1 (4.4.11)
¼ u_ n1 + c2 qnm
where
c2 ¼ h (4.4.12)
P
n1
Solving Eq. (4.4.11) for qrm and substituting into Eq. (4.4.8) gives
r¼1
Because m 6¼ 0, the coefficient matrix in Eq. (4.4.17) is not singular for suf-
ficient small h and the system can be solved successively for n ¼ 1, 2, … to yield
the solution un and the derivatives u_ n , u€n ¼ qn at instant t ¼ nh T . For n ¼ 1,
the value q0 appears in the right side of Eq. (4.4.17). This quantity can be readily
obtained from Eq. (4.4.1) for t ¼ 0. Thus, we have
q0 ¼ ðp0 cu_ 0 ku 0 Þ=m (4.4.18)
Eq. (4.4.17) can be also written as
Un ¼ AUn1 + bpn , n ¼ 1, 2, …,N (4.4.19)
in which
8 9
< qn =
Un ¼ u_ n (4.4.20a)
: ;
un
2 31 2 3
m c k 0 0 0
6 1 7 6 1 7
6 c1 h 1 7 6 c1 0 1 7
6
A¼6 2 7 6 7 (4.4.20b)
7 6 2 7
4 1 5 4 1 5
c2 1 0 c2 1 0
2 2
2 31
m c k 8 9
6 1 7 <1=
6 c h 1 7
b¼66 2
1 7
7 :0; (4.4.20c)
4 1 5 0
c2 1 0
2
The recurrence formula (4.4.19) can be employed to construct the solution
algorithm. However, the solution procedure can be further simplified. Thus,
applying Eq. (4.4.19) for n ¼ 1, 2, … we have
U1 ¼ AUo + bp1
U2 ¼ AU1 + bp2
¼ AðAUo + bp1 Þ + bp2
(4.4.21)
¼ A2 Uo + Abp1 + bp2
⋯ ¼ ⋯⋯⋯⋯⋯⋯⋯⋯⋯
Un ¼ An Uo + An1 p1 + An2 p2 + ⋯A0 pn b
Obviously, the last of Eq. (4.4.21) gives the solution vector Un at instant
tn ¼ nh using only the known vector U0 at t ¼ 0. The matrix A and the vector
b are computed only once.
Table 4.4.1 presents the algorithm for the numerical implementation of
AEM in pseudocode-type notation so that the reader can write a computer code
in the language of his/her preference.
@CivilMethod
Numerical integration of the equation of motion Chapter 4 173
Example 4.4.1 Response of an SDOF system using the analog equation method
Determine the response the SDOF system in Example 4.2.1 using the AEM.
Solution
The solution is obtained using the program aem_lin.m with Dt ¼ 0:01.
Fig. E4.5 gives the graphical representation of the displacement together with
the error u ðt Þ uex ðt Þ. Moreover, Fig. E4.6 shows the response of the system
under the static load p ¼ t, if 0 t 1 and p ¼ 1, if 1 < t.
@CivilMethod
174 PART I Single-degree-of-freedom systems
x 10−4
0.2 2
u(t) computed
u-uex
u(t) exact 1.5
0.15
1
0.1 0.5
u(t)
0
0.05
–0.5
–1
0
–1.5
–0.05
0 2 4 6 8 10 –2
0 2 4 6 8 10
t
× 10
-4 Displacement × 10
-4 Velocity
6 10
4 5
du(t)
u(t)
2 0
0 –5
0 5 10 0 5 10
t t
× 10
-3 Acceleration Response ratio
2 1.5
0 1
ddu(t)
R(t)
–2 0.5
–4 0
0 5 10 0 5 10
t t
FIG. E4.6 Response of the system in Example 4.4.1 under static load.
@CivilMethod
Numerical integration of the equation of motion Chapter 4 175
@CivilMethod
176 PART I Single-degree-of-freedom systems
We observe that un + 1 can be computed if the two initial values u1 and uo
are known. Eq. (4.5.1) may also be written as
Eq. (4.5.3) is a difference equation of the second order [14] whose solution
can yield un + 1 using the recursive procedure (4.5.2). Because Eq. (4.5.3) is
linear with respect to un1 , un , and un + 1 , it is called a linear difference equa-
tion. If gn ¼ 0 the difference equation is called homogeneous while if gn 6¼ 0 it is
called nonhomogeneous. In general, an equation of the form
rn + k + a1 rn + k1 + a2 rn + k2 + … + ak rn ¼ 0
is also a solution for all n for which the difference equation is defined with
c1 ,c2 ,…, ck being arbitrary constants. Eq. (4.5.6) is the general solution of
the difference equation (4.5.4). The arbitrary constants c1 ,c2 ,…, ck are deter-
mined from the k initial conditions.
If two roots of the polynomial (4.5.5) are complex conjugate, say r1 ¼
a + ib and r2 ¼ a ib, then we can write them in exponential form
where
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
r¼ a2 + b 2 , q ¼ tan 1 ðb=aÞ (4.5.8)
@CivilMethod
Numerical integration of the equation of motion Chapter 4 177
@CivilMethod
Numerical integration of the equation of motion Chapter 4 179
‘im ð1 + eÞ1=e ¼ e
e!0
Using this relation and taking into account that n ¼ tn =h, we obtain for a
given tn
@CivilMethod
180 PART I Single-degree-of-freedom systems
1
un ¼ c1 e2tn + (4.5.18)
2
Evidently, the first term in Eq. (4.5.17) is the exact solution. The second
term is spurious (extraneous) and results from the fact that the first-order dif-
ferential equation is substituted by a second-order difference equation. The
application of the initial conditions would give c2 ¼ 0 if the computations were
exact. In practice, however, errors are introduced, which are mainly due to the
rounding of numbers or the inaccuracy of the starting value. Therefore, the con-
stant c2 is not exactly zero and consequently, a small error is introduced in each
integration step. This is magnified because it is multiplied by the factor
ð1Þn e2tn , which increases exponentially. Because the first term of the solution
(4.5.17) diminishes exponentially, the introduced error due to the spurious solu-
tion dominates the exact solution and leads to a totally wrong result. We
describe a method as unstable if the error increases exponentially with tn .
For the first order differential equations, the one-step integration methods do
not exhibit instability for small values of h. The multistep methods, however,
which lead to difference equations of order greater than one, introduce spurious
solutions and they may be unstable either for all values of h or for a certain
region of the values of h. In order to decide whether a multistep method is sta-
ble, we work as follows.
If the multistep method leads to a difference equation of order k, we find the
roots of the characteristic equation. If ri (i ¼ 1,2, …,kÞ are these roots, the gen-
eral solution will be
One of the solutions, say rn1 , will tend to the actual solution of the differential
equation. The remaining roots are spurious. We will say that a multistep method
is strongly stable if for h ! 0 the spurious roots satisfy the condition
Because we do not know which is the actual solution, the above condition
should apply to all roots ri . Apparently, this condition ensures that the error
diminishes as n increases. On the contrary, the error increases exponentially
if jri j > 1.
@CivilMethod
Numerical integration of the equation of motion Chapter 4 181
where
2 h 2 w2 1 hwx pn h 2
a1 ¼ , a2 ¼ , gn ¼ (4.5.21)
2ð1 + hwx Þ 1 + hwx m ð1 + hwx Þ
The characteristic equation of Eq. (4.5.20) is
r2 2a1 r + a2 ¼ 0 (4.5.22)
whose roots are
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
r1, 2 ¼ a1 a21 a2
The type of root depends on the sign of the discriminant D ¼ a21 a2 , which
may be written as
pffiffiffiffiffiffiffiffiffiffiffiffi pffiffiffiffiffiffiffiffiffiffiffiffi
s2 s 2 1 x2 s + 2 1 x2
D ðs Þ ¼ , s ¼ hw
4ð1 + sxÞ2
pffiffiffiffiffiffiffiffiffiffiffiffis > 0, the sign of the discriminant depends only on the factor
Because
s 2 1 x2 . Hence, we distinguish the following two cases
qffiffiffiffiffiffiffiffiffiffiffiffi pffiffiffiffiffiffiffiffiffiffiffiffi qffiffiffiffiffiffiffiffiffiffiffiffi
2 1 x2 T
ðiÞ Dðs Þ > 0 s > 2 1 x 2
or h > ¼ 1 x2 (4.5.23)
w p
qffiffiffiffiffiffiffiffiffiffiffiffi pffiffiffiffiffiffiffiffiffiffiffiffi qffiffiffiffiffiffiffiffiffiffiffiffi
2 1 x2 T
ðiiÞ Dðs Þ 0 s 2 1 x 2
or h ¼ 1 x 2 (4.5.24)
w p
In case (i), the characteristic equation has two real roots, r1 , r1 , The stability
of the solution requires that jr1 j < 1 and jr2 j < 1. But it can be shown that
jr2 j > 1, hence jr1 j > 1. Consequently, the solution is unstable in this case.
In case (ii), the characteristic equation has two complex conjugate roots
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
r1, 2 ¼ a1 i a2 a21 (4.5.25)
This gives
sffiffiffiffiffiffiffiffiffiffiffiffiffi
pffiffiffiffiffi 1 sx
j r1 j ¼ j r2 j ¼ a2 ¼ 1 (4.5.26)
1 + sx
@CivilMethod
182 PART I Single-degree-of-freedom systems
@CivilMethod
Numerical integration of the equation of motion Chapter 4 183
2 31
1 2xw w2 8 9
6 7 >
> 1>
6 h2 7 < > =
6 1 7
b ¼ 6 4 0 7 0 (4.5.36b)
6 7 >
> >
4 h 5 : > ;
0
1 0
2
Applying Eq. (4.5.35) for n ¼ 1, 2, … we have
U1 ¼ AUo + b
p1
U2 ¼ AU1 + b
p2
¼ AðAUo + b
p1 Þ + b
p2
(4.5.37)
¼ A2 Uo + Ab
p1 + b
p2
⋯ ¼ ⋯⋯⋯⋯⋯⋯⋯⋯⋯
Un ¼ An Uo + An1 p1 + An2 p2 + …A0 pn b
The matrix A is known as the amplification matrix. The stability of the
method requires that An is bounded. This is true if the spectral radius of A sat-
isfies the condition [15]
rðAÞ ¼ max fjr1 j, jr2 j, jr3 jg 1 (4.5.38)
where ri (i ¼ 1, 2, 3) are the eigenvalues of the matrix A.
Using a symbolic language (here MATLAB) we find
pffiffiffiffiffiffiffiffiffiffiffiffi pffiffiffiffiffiffiffiffiffiffiffiffi
4 s2 + 4s x2 1 4 s 2 4s x2 1
r1 ¼ , r2 ¼ , r3 ¼ 0 (4.5.39)
4 + 4sx + s 2 4 + 4sx + s 2
where s ¼ xw.
The type of the roots r1 , r2 depends on the sign of the discriminant
DðxÞ ¼ x2 1. Hence, we distinguish the following two cases
(i) If D ðx Þ > 0, both eigenvalues are real. It can be shown that jr1 j < 1, hence
jr2 j < 1. Therefore, the method is stable
(ii) If DðxÞ 0, the eigenvalues are complex conjugate
pffiffiffiffiffiffiffiffiffiffiffiffi
4 s 2 i4s 1 x2
r1,2 ¼
4 + 4sx + s 2
and
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
ð4 s 2 Þ2 + 16s2 1 x 2 ðs 2 + 4Þ2 16s 2 x2
j r1 j ¼ j r2 j ¼ ¼ 1 (4.5.40)
ð4 + 4sx + s 2 Þ ½ðs 2 + 4Þ + 4sx 2
The equality is valid for x ¼ 0.
@CivilMethod
184 PART I Single-degree-of-freedom systems
The conclusion is that the AAM is stable without imposing any constraint on
the size of the time step. We say in this case that the integration method is
unconditionally stable.
Essentially, the procedure based on the condition (4.5.38) to prove the sta-
bility of the AAM is not different from that presented in Section 4.5.3, where the
stability results from the response of the difference equation. This is shown in
what it follows.
We write Eq. (4.5.35) for tn ,tn1 , tn2
Un AUn1 ¼ b
pn
Un1 AUn2 ¼ b
pn1 (4.5.41)
Un2 AUn3 ¼ b
pn2
or in matrix form
8 9
2 3> Un > 8 9
I A 0 >
0 < >
= < pn =
4 0 I A 0 5 Un1
¼ b pn1 (4.5.42)
> Un2 > : ;
0 0 I A >
: >
; pn2
Un3
Eq. (4.5.42), beside the displacements un ,un1 ,un2 , un3 , contains the
velocities u_ n , u_ n1 , u_ n2 , u_ n3 and the accelerations u€n , u€n1 , u€n2 , u€n3 . Reor-
dering these equations and eliminating the velocities and accelerations yield the
equation
un 2a1 un1 + a2 un2 + a3 un3 ¼ c1 pn + c2 pn1 + c3 pn2 (4.5.43)
where
4 w2 h 2 4 4sx + w2 h 2
a1 ¼ , a2 ¼ , a3 ¼ 0 (4.5.44)
4a 4a
1 + xwh 1 xwh
c1 ¼ , c2 ¼ , c3 ¼ 0 (4.5.45)
a a
where
w2 h 2
a ¼ 1 + xwh +
4
Eq. (4.5.43) is a difference equation whose characteristic equation is
r3 2a1 r2 + a2 r ¼ 0 (4.5.46)
The roots of Eq. (4.5.46) are
pffiffiffiffiffiffiffiffiffiffiffiffi pffiffiffiffiffiffiffiffiffiffiffiffi
4 s 2 + 4s x 2 1 4 s 2 4s x2 1
r1 ¼ , r2 ¼ , r3 ¼ 0 (4.5.47)
4 + 4sx + s 2 4 + 4sx + s 2
that is, they are identical to the eigenvalues of the matrix A.
@CivilMethod
Numerical integration of the equation of motion Chapter 4 185
@CivilMethod
186 PART I Single-degree-of-freedom systems
or
un ¼ r n ðc1 sin nq + c2 cos nqÞ
(4.6.4)
n + c2 cos wt
¼ r n ðc1 sin wt nÞ
pffiffiffiffiffiffiffiffiffiffiffiffiffiffi
where r ¼ a 2 + b2 , q ¼ tan 1 ðb=a Þ, a ¼ Reðr1 Þ, b ¼ Imðr1 Þ, w ¼ q=h,
tn ¼ nh.
The corresponding exact solution is
qffiffiffiffiffiffiffiffiffiffiffiffi
un ¼ exwtn ðc1 sin wD tn + c2 cos wD tn Þ, wD ¼ 1 x2 , tn ¼ nh (4.6.5)
Comparison of Eqs. (4.6.4), (4.6.5) could show the accuracy of the numer-
ical scheme. Thus, the period elongation is
pffiffiffiffiffiffiffiffiffiffiffiffi
T T h 1 x2
PE ¼ ¼ 1 (4.6.6)
T q
For the amplitude decay, we can define an equivalent damping ratio x from
the relation
n
r n ¼ exwt ¼ exqn (4.6.7)
which gives
x ¼ ln r=q (4.6.8)
The difference Dx ¼ x x can be employed as a measure for the amplitude
decay. The dependence of the period elongation and amplitude decay on h=T is
shown in Figs. 4.6.1 and 4.6.2, respectively. Obviously, for small values of h=T
the scheme is accurate. Note that for x ¼ 0 it is jr2 j ¼ jr3 j ¼ r ¼ 1 and Eq. (4.6.8)
yields x ¼ 0. That is, there is no amplitude decay.
0.25
x=0
0.2 x = 0.1
Period elongation %
x = 0.2
0.15
0.1
0.05
0
0 0.05 0.1 0.15 0.2 0.25 0.3
h/T
FIG. 4.6.1 Period elongation versus h=T for different values x.
@CivilMethod
Numerical integration of the equation of motion Chapter 4 187
−3
x 10
2
x=0
0 x = 0.1
x = 0.2
−2
Δx
−4
−6
−8
0 0.05 0.1 0.15 0.2 0.25 0.3
h/T
FIG. 4.6.2 Amplitude decay Dx ¼ x x versus h=T for different values of x.
Solution
The characteristic equation of Eq. (1) is
2 + s 2 r3 5r2 + 4r 1 ¼ 0 (2)
pffiffiffi pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi1=3
b ¼ 720a + 108a 2 + 1000 + 12 3a 104a + 100 + 27a 2 (5)
The stability condition requires that jr1 j < 1 and jr2 j ¼ jr3 j < 1. This
condition is satisfied, as is shown in Fig. E4.7.
@CivilMethod
188 PART I Single-degree-of-freedom systems
1
|r1|
0.9
|r2|
0.8 |r3|
0.7
0.6
0.5
|r2|=|r3|
0.4
0.3
0.2
0.1
| r1| (spurious)
0
1 2 3 4 5 6 7 8 9 10
h/t
FIG. E4.7 Houbolt’s method. Roots of the characteristic equation.
4.7 Problems
Problem P4.1 Find the solution of the IVPs for the difference equations
un + 1 3un ¼ 5, u0 ¼ 1
un + 2 4un + 1 + 3un ¼ 2n , u0 ¼ 0, u1 ¼ 1
Hint: The particular solution of (ii) will be sought in the form gn ¼ 2n c,
where c is a constant.
Problem P4.2 A SDOF system is subjected to pulse the load
t
pðt Þ ¼ p0 1 1 + et=t1 , 0 t t1
t1
pðt Þ ¼ 0, t > t1
Study the response of the system using
(i) The central difference method.
(ii) The average acceleration method.
(iii) The analog equation method.
Compare the results with the exact solution.
Data: u ð0Þ ¼ u_ ð0Þ ¼ 0, m ¼ 50kNm1 s2 , x ¼ 0:05, w ¼ 5s1 , p0 ¼ 100kN,
t1 ¼ 0:5s and ttot ¼ 10s.
Problem P4.3 A SDOF system is subjected to the piecewise linear load pðt Þ
shown in Fig. P4.3. Use the three discussed numerical methods to establish
its response.
@CivilMethod
Numerical integration of the equation of motion Chapter 4 189
Compare the results with the exact solution. Data: u ð0Þ ¼ u_ ð0Þ ¼ 0,
m ¼ 50kN m1 s2 , x ¼ 0:05, w ¼ 5s1 , ti ¼ 0:01i, pi ¼ 20 1 + ð1Þi ði + 5Þ=
ði + 1Þ , i ¼ 0, 2, …,100.
@CivilMethod
Chapter 5
Nonlinear response:
Single-degree-of-freedom
systems
Chapter outline
5.1 Introduction 191 5.4 The analog equation
5.2 The central difference method 195 method 203
5.3 The average acceleration 5.5 Problems 211
method 197 References and further reading 215
5.1 Introduction
The equation of motion of a vibrating system expresses the equilibrium condi-
tion of all forces applied to the system, namely the external excitation force, the
inertial force, the damping force, and the elastic force. The equilibrium condi-
tion reads
fI ðt Þ + fD ðt Þ + fS ðt Þ ¼ pðt Þ (5.1.1)
The forces fI , fD , and fS depend on the physical properties of the system.
In the systems we analyzed, the physical properties are not time-dependent
and the dependence of these forces on the cause that produces them is linear,
that is,
fI ¼ m u€ (5.1.2a)
fD ¼ cu_ (5.1.2b)
fS ¼ ku (5.1.2c)
where m, c,k are constant quantities. Systems with such a physical response are
referred to as linear systems.
In general, however, Eq. (5.1.2a)-(5.1.2c) may be of the form
fI ¼ m ðt Þu€ (5.1.3a)
_ tÞ
fD ¼ fD ðu, u, (5.1.3b)
_ tÞ
fS ¼ fS ðu, u, (5.1.3c)
That is, the mass may vary with time and the damping and elastic forces may
be nonlinear functions of u, u, _ t Þ, fS ðu, u,
_ and t, that is, fD ðu, u, _ t Þ. In this case,
Εq. (5.1.1) takes the form
_ t Þ + fS ðu, u,
m ðt Þu€ + fD ðu, u, _ t Þ ¼ pðt Þ (5.1.4)
which is a nonlinear differential equation of the second order. A dynamic sys-
tem, whose response is described by Eq. (5.1.4), is referred to as nonlinear.
Although systems with variable mass are not unusual [1], our discussion will
be limited to systems with constant mass. Besides, the forces fD and fS will
be considered of the form fD ðu_ Þ and fD ðu Þ. Thus, Eq. (5.1.4) becomes
m u€ + fD ðu_ Þ + fS ðu Þ ¼ pðt Þ (5.1.5)
We distinguish two types of nonlinearity: the geometric nonlinearity, which
is due to large displacements implying large deformations of the structure, and
the material nonlinearity, which is due to nonlinear constitutive equations (e.g.,
hyperelastic or elastoplastic materials). Of course, both types of nonlinearity
can simultaneously characterize the response of a system.
The analytical solution of the nonlinear equations of motion is a difficult and
complicated mathematical problem. Exact solutions are available only for a few
cases and for differential equations of a specific form [2, 3]. The existing solu-
tions aim rather at a qualitative study of the response of the system described by
a nonlinear equation than at offering a computational means for practical ana-
lyses. The knowledge of the nonlinear response of the single-degree-of-freedom
(SDOF) systems comes from approximate methods, and mainly from numerical
methods. Therefore, the recourse to numerical methods to solve the nonlinear
equations of motion is inevitable. The step-by-step methods play a dominant
role. The Runge-Kutta methods, usually employed for the solution of nonlinear
equations, belong to these methods [4].
The dynamic response of nonlinear systems can be studied effectively by
demanding the fulfillment of equation motion (5.1.5) at discrete time instants
Dt apart by the use of the step-by-step integration methods we discussed in
Chapter 4. These methods as developed for nonlinear equations of motion
are presented directly below while for the analytical methods, the reader is
advised to look in the vast related literature [2, 3, 5].
Example 5.1.1 Systems with a geometrical nonlinearity
Derive the equation of motion of the system shown in Fig. E5.1. The supports at
A,B and the interconnection at C are hinges. The mass m at C is concentrated.
The system is set to motion by the initial conditions u0 , u_ 0 and/or the vertical
external force pðt Þ. The bars are assumed massless.
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 193
l l
C
A φ B
u(t )
E, A E, A
m
p(t )
fI
S fS S
p(t )
FIG. E5.1 SDOF system with nonlinear response.
Solution
During the motion, the force pðt Þ, the inertial force fI , and the elastic forces of
the bars are in equilibrium
The inertial force is given by the relation
fI ¼ m u€ (1)
The total elastic force is caused by the elongation of the bars and is given by
fS ¼ 2S sin f (2)
where
EA
S¼ d (3)
l
E is the modulus of elasticity of the material of the bar, A its cross-sectional area,
and d the elongation of the bars.
The elongation of the bars at time t is
pffiffiffiffiffiffiffiffiffiffiffiffiffiffi
d¼ l 2 + u2 l (4)
and
u
sin f ¼ pffiffiffiffiffiffiffiffiffiffiffiffiffiffi (5)
l 2 + u2
@CivilMethod
194 PART I Single-degree-of-freedom systems
Hence, the elastic force in Eq. (2) by virtue of Eqs. (3)–(5) is expressed as
2 3
u6 1 7
fS ðu Þ ¼ 2EA 41 qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi5 (6)
l 2
1 + ðu=l Þ
Obviously, the equation of motion is nonlinear. Fig. E5.2 shows the graph
of fS ðu Þ. Because d 2 fS =du 2 > 0 the curve is concave upward. In this system,
the slope kT increases continually, which implies that the elastic force
increases with increasing u. In this case, we say that the system exhibits hard-
ening, in contrast to other systems that exhibit softening. In the latter systems,
the curve fS ðu Þ is concave downward (d 2 fS =du 2 < 0), that is, the slope
kT decreases continually, which implies that elastic force decreases with
increasing u, for example, a system with stiffness fS ¼ 40u u 3 exhibits soft-
ening; see Fig. E5.3.
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 195
@CivilMethod
196 PART I Single-degree-of-freedom systems
where k ¼ sin 2 ðq0 =2Þ and T0 is the quarter of the period; sn represents the
sn-Jacobean elliptic function [7].
Solution
The response of the pendulum for l ¼ g, q0 ¼ 0:40p, and q_ 0 ¼ 0 is obtained using
the program centr_diff_nlin.m with Dt ¼ 0:01. It is shown in Fig. E5.5 as com-
pared with the exact one.
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 197
t t
(a) (b)
FIG. 5.3.1 Tangent and secant of (a) damping and (b) stiffness curve.
DfS ¼ fS ðt + Dt Þ fS ðt Þ
¼ fS ðu + Du Þ fS ðu Þ
(5.3.4)
df 1 d 2 fS
¼ S Du + ðDu Þ2 + ⋯
du 2 du 2
For small values of Dt, the quantities Du, Du_ are also small. Thus, neglect-
ing the nonlinear terms in Eqs. (5.3.3), (5.3.4), we obtain
DfD cT Du_ (5.3.5)
DfS kT Du (5.3.6)
Obviously, cT and kT express the slope of the tangent to the curves fD ðu_ Þ
and fS ðu Þ, respectively, at time t.
Referring to Fig. 5.3.1a, we have
fD ðt + Dt Þ ¼ fD ðt Þ + DfD
However, the exact value of DfD is
DfD ¼ Du_ tan fc
where fc is the angle of the secant. Hence, approximating DfD by Eq. (5.3.5)
introduces the error (see Fig. 5.3.1a)
eD ¼ ðcT tan fc ÞDu_ (5.3.7)
because cT is the slope of the tangent.
Similarly, the use of kT to approximate DfS introduces the error (see
Fig. 5.3.1b)
eS ¼ ðkT tanfk ÞDu (5.3.8)
The errors eD and eS cannot be avoided because u_ ðt + Dt Þ and u ðt + Dt Þ
are not known at instant t + Dt. However, as we will show, they can be kept
under a given bound, which specifies the accuracy of the solution procedure.
By virtue of Eqs. (5.3.5), (5.3.6), Eq. (5.3.2) is written in incremental form
mDu€ + cT Du_ + kT Du ¼ Dp (5.3.9)
The previous equation is of the form (4.3.11). Hence, the AAM is suitable to
solve it. Thus, using Eqs. (4.3.9), (4.3.10) to express Du€ and Du,
_ we obtain
2
Du_ ¼ Du 2u_ ðt Þ (5.3.10)
Dt
4
Du€ ¼ ½Du Dt u_ ðt Þ 2u€ðt Þ (5.3.11)
Dt 2
Then Eq. (5.3.9) becomes
2cT 4m 4m
kT + + 2 Du ¼ Dp + + 2cT u_ ðt Þ + 2m u€ðt Þ (5.3.12)
Dt Dt Dt
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 199
or
k ∗ Du ¼ Dp∗ (5.3.13)
where
2cT 4m
k ∗ ¼ kT + + 2 (5.3.14)
Dt Dt
and
4m
Dp∗ ¼ Dp + + 2cT u_ n + 2m u€n (5.3.15)
Dt
The value of Du obtained from Eq. (5.3.13) is used in Eq. (5.3.10) to eval-
_ Then we obtain
uate Du.
un + 1 ¼ un + Du (5.3.16a)
u_ n + 1 ¼ u_ n + Du_ (5.3.16b)
where
DfS DfD
ksec ¼ , csec ¼ (5.3.18)
Du Du_
and a is a small specified number defining the upper bound of the error, for
example, a ¼ 0:01. If it is
@CivilMethod
200 PART I Single-degree-of-freedom systems
When the damping force depends linearly on the velocity, then the response
is governed by Eq. (5.2.1). In this case, the previously presented incremental
method can be improved by employing an iterative procedure within each step,
which minimizes the error introduced by the tangent stiffness kT .
The starting point is Eq. (5.3.13), which we write as
k ∗ du ð1Þ ¼ Dp∗ (5.3.20)
The quantity du ð1Þ is the first approximation to Dun within the time step
from tn to tn + 1 . That is
Dunð1Þ ¼ du ð1Þ (5.3.21)
The index n denotes the number of the step that brings us from the displace-
ment un to un + 1 . Fig. 5.3.2 presents the graph of the function p∗ ðu Þ. Apparently,
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 201
the displacement du ð1Þ resulting from Eq. (5.3.20) decreases Dp∗ by dpð1Þ .
Hence, there is a remaining force
DF ð2Þ ¼ Dp∗ dpð1Þ (5.3.22)
which must be equilibrated. The change dpð1Þ is computed using Eq. (5.3.20), if
the tangential slope kT is replaced with the slope of the secant
" #
ð1Þ fS un + du ð1Þ fS ðun Þ 2c 4m
dp ¼ + + du ð1Þ
du ð1Þ Dt Dt 2
2c 4m
¼ fS un + du ð1Þ fS ðun Þ + + 2 du ð1Þ
Dt Dt
which by virtue of Eq. (5.3.14) becomes
dpð1Þ ¼ fS un + du ð1Þ fS ðun Þ + ðk ∗ kT Þdu ð1Þ (5.3.23)
The force DF ð3Þ produces the additional displacement du ð3Þ , which is com-
puted from the relation
k ∗ du ð3Þ ¼ DF ð3Þ (5.3.28)
Thus, the new approximation of Dun is
Dunð3Þ ¼ Dunð2Þ + du ð3Þ (5.3.29)
Consequently, for the i + 1 approximation it is
DF ði + 1Þ ¼ DF ðiÞ dpðiÞ (5.3.30)
dpðiÞ ¼ fS un + DunðiÞ fS un + Dunði1Þ + ðk ∗ kT Þdu ðiÞ (5.3.31)
@CivilMethod
202 PART I Single-degree-of-freedom systems
k ∗ du ði + 1Þ ¼ DF ði + 1Þ (5.3.32)
Dunði + 1Þ ¼ DunðiÞ + du ði + 1Þ (5.3.33)
Note that for i ¼ 1 it must be set Dunð0Þ ¼ 0 and DF ð1Þ ¼ Dp∗ .
The iteration procedure is terminated after I iterations, if
du ði + 1Þ
e (5.3.34)
DunðI Þ
where e is a specified small number. Then we assume that the convergence has
been achieved. The value Dun ¼ DunðI Þ is considered exact and it is used to
compute un + 1 and u_ n + 1 , u€n + 1 . Subsequently, the procedure continues to
the next step. The iterative procedure within the time step from tn to tn + 1 is
summarized in Table 5.3.2. This procedure is known as the modified
Newton-Raphson method.
5. du ði + 1Þ ¼ DRði + 1Þ =k ∗
6. If du ði + 1Þ =DunðI Þ > e set i ¼ i + 1 and go to B.2
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 203
solution procedure are presented in Table 5.4.1. Adhering to the steps of this
table, a computer program called aem_nlin.m has been written in MATLAB
for the numerical integration of the nonlinear equation of motion using the
AEM. The program is available on this book’s companion website. It computes
the displacement u ðt Þ, the velocity u_ ðt Þ, and the acceleration u€ðt Þ.
Example 5.4.1 The Duffing equation. AEM solution
Use the AEM to solve the IVP for the Duffing equation
u€ + 0:2u_ + u + u 3 ¼ pðt Þ (1)
u ð 0Þ ¼ 0 (2a)
u_ ð0Þ ¼ 1 (2b)
For
pðt Þ ¼ e0:1t ½ð0:01sin t 0:2cost sin t Þ 0:2ð0:1 sin t cos t Þ
+ sin t + e0:2t ð sin t Þ3
Eq. (1) admits an exact solution uexact ðt Þ ¼ e0:1t sin t.
Solution
The solution is obtained using the program aem_nlin.m with Dt ¼ 0:01. The
graph of the solution is shown in Fig. E5.6 as compared with the exact one.
1
Computed
0.8 Exact
error X103
0.6
0.4
0.2
−0.2
−0.4
−0.6
−0.8
0 5 10 15 20 25
t
FIG. E5.6 Solution u and error u uexact in Example 5.4.1.
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 205
which is solved for the data of the problem to yield ust ¼ 0:25m.
The external force is the constant weight of the body, that is, pðt Þ ¼ mg,
while the elastic force
fS ðust + u Þ (2)
where u ¼ u ðt Þ denotes the additional displacement due to the dynamic
response. Thus, the equation of motion becomes
m u€ + fS ðust + u Þ ¼ mg (3)
with initial conditions
u0 ¼ 0:05m, u_ 0 ¼ 0 (4)
It should be noted that due to the nonlinearity of the elastic force, the super-
position of the displacements does not apply. The computed response of the sys-
tem with Dt ¼ 0:01 is shown in Fig. E5.7.
(ii) In this case the equation of motion reads
m u€ + fS ðu Þ ¼ mg (5)
with initial conditions
u0 ¼ 0, u_ 0 ¼ 0 (6)
where u ¼ u ðt Þ denotes the total dynamic displacement from the undeformed
position. The computed response of the system with Dt ¼ 0:01 is shown in
Fig. E5.8.
(iii) In this case, the IVP becomes
m u€ + fS ðu Þ ¼ mg + 10sin Wt (7)
u0 ¼ 0, u_ 0 ¼ 0 (8)
The computed response with Dt ¼ 0:01 is shown in Fig. E5.9 for the ratio
pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
W=w ¼ 0:1, where w ¼ EA=lm . In all cases, the results are compared with
those obtained by the CDM
@CivilMethod
206 PART I Single-degree-of-freedom systems
u(t)
t
FIG. E5.7 Response in Example 5.4.2 (i).
u(t)
t
FIG. E5.8 Response in Example 5.4.2 (ii).
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 207
@CivilMethod
208 PART I Single-degree-of-freedom systems
(a) (b)
FIG. E5.11 System with a material nonlinearity in Example 5.4.4.
Solution
It is
fSy 69:68
uy ¼ ¼ ¼ 0:013m ¼ yield displacement
k 5360:00
um ¼ maximum displacement where the velocity changes sign
uR ¼ um uy ¼ remaining plastic deformation
The equation of motion is
m u€ + fS ðu Þ ¼ pðt Þ
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 209
@CivilMethod
TABLE E5.1 Numerical solution of elastoplastic system in
Example 5.4.4.—cont’d
u ðt Þ (cm)
Dt ¼ 0:01 Dt ¼ 0:001
t (s) AEM AAM AEM AAM
0.35 1.134 1.234 1.177 1.175
0.40 0.378 0.487 0.423 0.421
0.45 0.353 0.262 0.312 0.313
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 211
5.5 Problems
Problem P5.1 Show that the motion of the simple pendulum is governed by
the IVP
g
q€ + sin q ¼ 0, qð0Þ ¼ q0 , q_ ð0Þ ¼ q_ 0 (1)
l
Solve the equation of motion numerically when g=l ¼ 1, q0 ¼ 0:1p, q_ 0 ¼ 0
and compute the period T of the pendulum. Give the graphical representation
of the function T ¼ T ðq0 Þ, 0:1 q0 1 for the two time steps Dt ¼ 0:1 and
Dt ¼ 0:0001. Compare with the exact expression
sffiffiffi Z
l p=2 df sin ðq=2Þ
T ðq 0 Þ ¼ 4 pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi , k ¼ sin ðq0 =2Þ, sin f ¼ (2)
g 0 1 k 2 sin 2 f sin ðq0 =2Þ
@CivilMethod
212 PART I Single-degree-of-freedom systems
Problem P5.3 Study the response of the system shown in Fig. P5.3, when
(i) u ð0Þ ¼ 0:05m, u_ ð0Þ ¼ 0, pðt Þ ¼ 0. The initial displacement will be taken
from the position of the static equilibrium.
(ii) u ð0Þ ¼ u_ ð0Þ ¼ 0, pðt Þ ¼ mg ð1 t=t1 ÞH ðt1 t Þ:pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
(iii) u ð0Þ ¼ u_ ð0Þ ¼ 0, P ðt Þ ¼ 15sin wt,
where w ¼ kmin =m .
(iv) Plot of the dependence of the period T as a function of u0 , T ¼ T ðu0 Þ, if
u_ ð0Þ ¼ 0 and pðt Þ ¼ 0. The initial displacement will be taken from the
position of the static equilibrium.
Data: t1 ¼ 1s, ttot ¼ 5s, m ¼ 10kN m1 s2 , E ¼ 2:1 108 kN=m2 , A ¼ 1cm2 ,
L ¼ 4:0m, and g ¼ 9:81 ms2 . The cables are assumed massless to avoid the
sag due to self-weight [10–12].
Problem P5.4 Study the response of the system shown in Fig. P5.3 if the elastic
force is approximated by the first three terms of its Tailor series. Compare the
results with those in Problem P5.3 when pðt Þ ¼ mg ð1 t=t1 ÞH ðt1 t Þ and
u ð0Þ ¼ u_ ð0Þ ¼ 0.
Problem P5.5 The response of a nonlinear system is governed by the following
equation, known as the van der Pol equation
u€ m 1 u 2 u_ + u ¼ 0
Plot the solution u ðt Þ, if u ð0Þ ¼ 0, u_ ð0Þ ¼ 0:1 and for (i) m ¼ 0:2 (ii) m ¼ 1:2.
Problem P5.6 The horizontal beams of the frame in Fig. P5.6 have negligible
mass and are flexible while the shear walls having a uniform mass density with
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 213
specific weight g are assumed rigid. Their support on the ground is elastic and it
is expressed by the relation
1
MR ¼ CR f f2 (a)
4
The system is set in motion by the initial conditions u0 ¼ 5cm, u_ 0 ¼ 0. Plot
the dependence of the ratio T=T0 (T0 is the period of the structure resulting
when the nonlinear term in Eq. a is ignored) on the initial displacement and dis-
cuss the influence of the nonlinearity of the elastic supports on the response of
the structure. Data: a ¼ 5m, g ¼ 24kN=m3 , E¼2.1107kN/m2, CR ¼ EI =5a,
and cross-sectional dimensions of the beam 0:20 0:40m2 .
Problem P5.7 The buoy of Fig. P5.7 consists of two massless cones with a base
dimeter 2R and a height h. A concentrated mass m attached at the bottom of the
body keeps the buoy floating at the position shown in the figure. Study the
@CivilMethod
214 PART I Single-degree-of-freedom systems
vertical motion of the buoy if it is displaced vertically downward from the equi-
librium position by u0 . Data: m ¼ 10kNm1 s2 h ¼ 5m, R ¼ 2m,
u ð0Þ ¼ 0:30m, u_ ð0Þ ¼ 0, and specific weight of the liquid g ¼ 2kN=m3 .
Problem P5.8 The water tower of Fig. P5.8a is subjected to the blast load of
Fig. P5.8b. The response of the structuren is elastoplastic.
h The
iorestoring force in
the elastic branch is given by fS ¼ ku 1 + 1= 1 + ð10u Þ6 =2 (Fig. P5.8c).
The structure is modeled by a SDOF system. Study the response of the structure
in the interval of t 0 ¼ 1s. Data: m ¼ 50kNm1 s2 , k ¼ 2000kN=m, x ¼ 0:07,
p0 ¼ 25kN, t1 ¼ 0:1s, and yield displacement uy ¼ 0:1m.
(b)
(a)
{ }
(c)
FIG. P5.8 Water tower in problem P5.8.
@CivilMethod
Nonlinear response: Single-degree-of-freedom systems Chapter 5 215
@CivilMethod
Chapter 6
6.1 Introduction
Structural systems are often excited by the motion of their support. The
response of a structure to support excitation is dynamic even though no external
dynamic loads act on it. The seismic motion of the ground represents a typical
example of support excitation of structures. The study of the response of struc-
tures to earthquake-induced motion is a specific but very important subject of
structural dynamics. It is discussed in depth in books on earthquake engineering
as well as in books on structural dynamics, preparing engineers to design struc-
tures for earthquake-induced motion [1,2]. This book treats the dynamic
response of structures when the excitation force is known. Therefore, the dis-
cussion in this chapter is limited only to the study of the dynamic response
of the SDOF system due to support excitation. Besides, some basic concepts such
as the response spectrum concept, which facilitates the dynamic analysis of struc-
tures due to ground motion, are presented. The general problem of the support
excitation of structures will be examined later when the MDOF (multi-degree-
of-freedom) systems are studied. The transmission of vibrations from the
structure to the fundament and vice versa are also discussed. Illustrative examples
analyzing the response of SDOF systems due to ground motion are presented. The
pertinent bibliography with recommended references for further reading is also
included. The chapter is enriched with problems to be solved.
utot
m ug u
Fixed axis of reference
c m c
k k
(a) ug (t ) (b)
m
fI mutot fD cu
fS ku
(c)
FIG. 6.2.1 SDOF system subjected to ground motion (a), Dynamic model (b), Forces on the
free body (c).
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 219
max ug,tt=258.59
ug,tt (cm/s2)
t
FIG. 6.2.2 Accelerogram from Athens earthquake, Sept. 7, 1999 (Recorded by ITSAK, Α399-1.
V2, longitudinal component, max u€g ¼ 258:59cm=s2 ).
@CivilMethod
220 PART I Single-degree-of-freedom systems
Apparently, Eq. (6.2.7b) states that the deformation of the system due to
given ground acceleration u€g ðt Þ depends only on the natural frequency w, hence
on the natural period T ¼ 2p=w, and on the damping ratio x, that is,
u ¼ u ðt, T , x Þ. Consequently, two systems with the same natural period T
and the same damping ratio x will undergo the same displacement u ðt Þ under
the same ground motion, in spite of the fact that the two systems may have dif-
ferent masses or different stiffnesses.
The negative sign in the effective load pðt Þ ¼ m u€g ðt Þ affects only the
direction of the displacement and not its magnitude. In practice, this has little
significance inasmuch as the engineer is usually interested in the maximum
absolute value of u ðt Þ. Therefore, the sign can be omitted in this case. This
assumption allows us to write the Duhamel integral in the forma
1
u ðt Þ ¼ U ðt Þ (6.2.8)
wD
where
Z t
U ðt Þ ¼ u€g ðτÞexwðtτÞ sin wD ðt τÞdτ (6.2.9)
0
For a given ground motion and a fixed damping ratio, we can evaluate the
largest absolute value of the function U ðt Þ, hence of u ðt Þ, for an interval of
values of the natural period T of the damped SDOF system and plot the curves
u ðT , xÞ ¼ max t ju ðt, T , xÞj for discrete values of the damping ratio x. Fig. 6.2.3
shows the curves u ðT , xÞ ¼ U ðT , x Þ=wD , 0 < T 2 for different values of x,
when the ground motion is induced by the Athens earthquake, Sept. 7, 1999.
The respective accelerogram is shown in Fig. 6.2.2. The curves u ðT , x Þ were
obtained by direct solution of Eqs. (6.2.7a), (6.2.7b). The solution can be
obtained using either the analytic solution presented in Section 3.5.4 or numer-
ically using any of the methods presented in Chapter 4.
pffiffiffiffiffiffiffiffiffiffiffiffi
a. In Eqs. (6.2.8) and (6.2.9), we can set wD ¼ w 1 x 2 w, because in real structures the value of
the damping ratio is small (x ¼ 3% 15%, hence x2 ≪1) and the error due to this approximation is
much smaller than that due to the uncertainty of the determination of u€g ðt Þ.
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 221
u T
T (s)
FIG. 6.2.3 Deformation response spectrum u ðT , x Þ for Athens earthquake, Sept. 7, 1999.
@CivilMethod
222 PART I Single-degree-of-freedom systems
T (s)
FIG. 6.2.4 Response spectrum of the relative velocity u_ ðT , xÞ for the Athens Earthquake,
Sept. 7, 1999.
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 223
T (s)
FIG. 6.2.5 Response spectrum of the total acceleration u€tot ðT , xÞ for the Athens earthquake,
Sept. 7, 1999.
For small values of the damping ratio (say, 0 x 0:15), we may set x 2 0.
Thus, the last two equations become
Qo ¼ mwU ðt Þ (6.2.16)
max Qo ¼ mwU ðT , x Þ (6.2.17)
We observe that the quantity wU ðT , x Þ in Eq. (6.2.17) has dimensions of
acceleration. In earthquake engineering, this quantity is designated by
Spa ðT , x Þ and it is known as spectral pseudoacceleration. It is an important
quantity because it allows direct evaluation of the maximum elastic force (base
shear force) in the SDOF system from the graph of Spa ðT , xÞ.
The quantity Spv ðT , x Þ ¼ U ðT , x Þ has dimensions of velocity and it is
known as spectral pseudovelocity. Actually, the quantities Spv ðT , x Þ,
Spa ðT , x Þ are different from u_ ðT , xÞ and u€tot ðT , x Þ. Therefore, they should
not be confused. Nevertheless, it is Spa ðT , xÞ ¼ u€tot ðT , xÞ, if x ¼ 0. Indeed,
Eq. (6.2.3) for x ¼ 0 becomes
u€tot ¼ w2 u ¼ wU ðt Þ (6.2.18)
from which we obtain
u€tot ðT , x Þ ¼ Spa ðT , x Þ (6.2.19)
The deviation of the pseudoacceleration from the extreme value of the total
acceleration is small for small values of x (say 0 x 0:1). Thus we may set
Spa ðT , x Þ u€tot ðT , x Þ. This is shown in Fig. 6.2.6.
@CivilMethod
224 PART I Single-degree-of-freedom systems
(utot),tt(T,0)
Spa(T,0)
T (s)
FIG. 6.2.6 Response spectra Spa ðT , x Þ and u€tot ðT , x Þ, ðx ¼ 0Þ.
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 225
@CivilMethod
226 PART I Single-degree-of-freedom systems
Sd
(a)
(b)
Spa
(c)
FIG. 6.2.7 Response spectra for the Athens earthquake, Sept. 7, 1999 (x ¼ 0:1).
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 227
5m
z
2.5 m
x utot
1m P ug u
m c
10 m
k
x
(a) (b)
FIG. E6.1 Silo in Example 6.2.1
@CivilMethod
228 PART I Single-degree-of-freedom systems
Hence
c ¼ 2mwx ¼ 2 40:6 5:57 0:04265 ¼ 19:3kNm1 s
qffiffiffiffiffiffiffiffiffiffiffiffi pffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
wD ¼ w 1 x 2 ¼ 5:57 1 0:042652 ¼ 5:56s1
2p
to ¼ 6T ¼ 6 ¼ 6:78s
wD
t (s)
FIG. E6.2 Time history of the response ratio Rðt Þ in Example 6.2.1.
@CivilMethod
230 PART I Single-degree-of-freedom systems
and
12EI k
max Q ¼ 3
max u ¼ max uII ¼ 1:781kN (6a)
h 4
h
max M ¼ max Q ¼ 8:907kNm (6b)
2
(a) (b)
FIG. E6.3 SDOF system in Example 6.3.1
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 231
Solution
The function representing the ground motion is shown in Fig. E6.3b. The prob-
lem can be solved by the following two ways:
(a) Formulating the equation of motion in terms of the total displacement, that
is, Eq. (6.3.3). This yields
m u€tot + ku tot ¼ ku 0 (1)
This implies that an effective constant excitation force pðt Þ ¼ ku 0 is
suddenly applied. The solution of Eq. (1) is obtained from Eq. (3.4.3)
for p0 ¼ ku 0 . Hence
pffiffiffiffiffiffiffiffiffi
utot ðt Þ ¼ u0 ð1 cos wt Þ, w ¼ k=m (2)
We observe that
ðutot Þmax ¼ 2u0
The relative displacement is
u ðt Þ ¼ utot ug ¼ u0 cos wt (3)
which yields an elastic force
fS ¼ ku ¼ ku 0 cos wt (4)
and
max jfS j ¼ ku 0 (5)
(b) Formulating the equation of motion in terms of the relative displacement
u ðt Þ. In this case, we have
ug ð0Þ ¼ u0 , u_ g ð0Þ ¼ 0, u€g ð0Þ ¼ 0
Consequently
pðt Þ ¼ 0, u ð0Þ ¼ utot ð0Þ ug ð0Þ ¼ u0 , u_ ð0Þ ¼ u_ tot ð0Þ u_ g ð0Þ ¼ 0
and the equation of motion becomes
m u€ + ku ¼ 0 (4)
with initial conditions u ð0Þ ¼ u0 , u_ ð0Þ ¼ 0.
The solution of Eq. (4) is obtained from Eq. (2.2.13) as
u ðt Þ ¼ u0 cos wt (5)
which is identical with that given by Eq. (3).
Example 6.3.2 The supports of the columns of the one-story frame of Fig. E6.4
are subjected to the displacements ug1 ¼ uo sin ðwt qÞ and ug2 ¼ uo sin wt.
Determine the equation of motion of the structure and give the expressions
of the relative displacement u ðt Þ and the stress resultants Q ðt Þ, M ðt Þ at the
top cross-sections of the columns. The dead load of the rigid beam is included
@CivilMethod
232 PART I Single-degree-of-freedom systems
utot (t ) utot (t )
p
EI
h
35 × 70 40 × 80
ug 1 ug 2
L
FIG. E6.4 One-story frame in Example 6.3.2
in the load p. The material of the columns is reinforced concrete. Data: specific
weight of concrete g ¼ 24kN=m3 , u0 ¼ 0:02m, x ¼ 0:05, w ¼ 2:5rad=s,
u ð0Þ ¼ 0, u_ ð0Þ ¼ 0, L ¼ 15m, h ¼ 7m, E ¼ 2:1 107 kN=m2 , and
p ¼ 200kN=m. The mass of the columns is assumed lumped at their ends.
Solution
The system has one degree of freedom. The equation of motion with respect to
the total displacement utot ðt Þ is
m u€tot + c1 u_ tot u_ g1 + c2 u_ tot u_ g2 + k1 utot ug1 + k2 utot ug2 ¼ 0 (1)
or
m u€tot + cu_ tot + ku tot ¼ pðt Þ (2)
where
c ¼ c1 + c2 , k ¼ k1 + k2 (3)
p ¼ c1 u_ g1 + c2 u_ g2 + k1 ug1 + k2 ug2 (4)
The mass of the system is
200 15 + ð0:35 0:70 + 0:40 0:80Þ ð7=2Þ 24
m¼ ¼ 310:6kNm1 s2
9:81
The stiffness of the columns
12EI 1 12 2:1 107 0:70 0:353
k1 ¼ ¼ ¼ 1837:5kN=m
h3 12 73
3EI 2 3 2:1 107 0:80 0:403
k2 ¼ ¼ ¼ 783:7kN=m
h3 12 73
Hence the stiffness of the system is
k ¼ k1 + k2 ¼ 2621:2kN=m
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 233
t
FIG. E6.5 Total displacement for different values of the phase angle q in Example 6.3.2 (a ¼ 0:2).
@CivilMethod
234 PART I Single-degree-of-freedom systems
utot(t) (m)
t
FIG. E6.6 Total displacement for different values of the allocation coefficient a in Example 6.3.2
(q ¼ 3p=4).
max|utot(t)|
a
FIG. E6.7 Extreme value max|utot ðt Þ| versus the allocation coefficient a in Example 6.3.2
(x ¼ 0:05).
From the study of the numerical results, we may draw the following
conclusion:
The percentage allocation of the damping to the two columns does not affect
significantly the dynamic response of the structure. Consequently, an arbitrary
but reasonable allocation, for example, 0:4 < a < 0:6, allows treating practical
cases of asynchronous support excitations.
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 235
@CivilMethod
236 PART I Single-degree-of-freedom systems
examined when the motion is in the steady state phase. Hence, the vertical dis-
placement is given by Eq. (3.2.26), that is,
qÞ
u ðt Þ ¼ rsin ðwt (6.4.1)
where
p0 h 2 i 1
2
r¼ 1 b2 + ð2xb Þ2 (6.4.2a)
k
2xb
q ¼ tan 1
, b ¼ w=w (6.4.2b)
1 b2
The total force transmitted to the foundation is
f ¼ fS + fD
¼ ku + cu_
(6.4.3)
qÞ + crw cos ðwt
¼ krsin ðwt qÞ
q fÞ
¼ fT sin ðwt
where
qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
fT ¼ r k 2 + ðcwÞ2 (6.4.4)
1 cw
f ¼ tan (6.4.5)
k
Using Eqs. (6.4.2a), (6.4.2b) for the expression of r and c=k ¼ 2x=w, the
previous relations are written as
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
u
u 1 + ð2xbÞ2
fT ¼ p0 t 2 (6.4.6)
1 b2 + ð2xb Þ2
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 237
TR=fT /p0
@CivilMethod
238 PART I Single-degree-of-freedom systems
m w2 u0 1
u ðt Þ ¼ 2 2xb cos wt
1 b2 sin wt
k 1 b + ð2xbÞ
2 2
(6.4.10)
u0
¼ 2xb3 cos wt
1 b2 b 2 sin wt
2 2 2
1 b + ð2xbÞ
and the total displacement
utot ¼ ug + u
u0
+
¼ u0 sin wt 2
2xb3 cos wt
1 b2 b 2 sin wt
1 b2 + ð2xbÞ2
u0 nh i o
¼ 2
1 b2 + ð2xb Þ2 sin wt
2xb 3 cos wt
1 b2 + ð2xbÞ2
Þ
¼ uT sin ðwt
(6.4.11)
where
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
u
u 1 + ð2xbÞ2
uT ¼ u0 t 2 (6.4.12)
1 b2 + ð2xb Þ2
!
1 2xb3
¼ tan (6.4.13)
1 b 2 + ð2xb Þ2
From Eqs. (6.4.8), (6.4.12), we obtain
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
u
ðu€tot Þmax ðutot Þmax uT u 1 + ð2xbÞ2
TR ¼ ¼ ¼ ¼ t 2 (6.4.14)
u€g max ug max u0 1 b2 + ð2xb Þ2
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 239
W 30kN
m¼ ¼ ¼ 3:058kNm1 s2
g 9:81m=s2
rffiffiffiffiffi rffiffiffiffiffiffiffiffiffiffiffi
k 28
w¼ ¼ ¼ 3:026s1
m 3:058
The circular frequency of the harmonic motion and the maximum
acceleration are
w ¼ 2pf ¼ 94:248s1
u€g max ¼ w2 u0 ¼ 177:653 ¼ 18:109g g ¼ 9:81ms2
hence
w
b ¼ ¼ 31:147
w
and
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
u
u 1 + ð2xbÞ2
TR ¼ t 2 ¼ 0:0065
1 b 2 + ð2xbÞ2
For the same elastic layer, the stiffness k and the damping coefficient c
are unaltered. Consequently, the transmission ratio TR can be modified if
the ratio b is changed. But because w is prescribed, this ratio can be chan-
ged, only if the natural frequency w is changed. This is possible if the mass
of the system (instrument) is changed by an increment Dm.
Let m 0 ¼ m + Dm be the new mass and x 0 , w0 the new damping ratio and
the frequency of the system (instrument), respectively. Then we will have
rffiffiffiffiffiffi
w0 m
¼ (1)
w m0
c ¼ 2mwx ¼ 2m 0 w0 x0
@CivilMethod
240 PART I Single-degree-of-freedom systems
which yield
rffiffiffiffiffiffi
0 mw m
x ¼x 0 0 ¼x (2)
mw m0
Moreover, we have
rffiffiffiffiffiffi
0 w m0
b ¼ 0¼b (3)
w m
From Eqs. (2), (3) we obtain
x0 b0 ¼ xb ¼ 3:115
Consequently, it must be
vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi vffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
u u
u 1 + ð2x 0 b0 Þ
2
u 2
TR ¼ u
¼ u 1 +
ð2 3:115Þ < 0:0055
t 2 2
t 2 2
1 b0 + ð2x0 b 0 Þ 1 b0
2
+ ð2 3:115Þ2
which holds if
pffiffiffiffiffiffiffiffiffiffiffiffi
33:885 < b0 ¼ b m 0 =m or m 0 > 3:619kNm1 s2
Hence Dm 0:561kNm1 s2 .
6.5 Problems
Problem P6.1 The supports 1 and 2 of the structure shown in Fig. P6.1 are sub-
and ug2 ðt Þ ¼ u0 sin ðwt
jected to the motions ug1 ðt Þ ¼ u0 sin wt p=3Þ, respec-
tively. A plane square rigid body of side a and density g is attached to
node 3. The mass of the column and the beam is neglected. Determine the
response of the structure and the reactions at the supports as well as the forces
that produce the support excitations. Data: L ¼ 8m, a ¼ 2m, cross section of the
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 241
(ug)i
(ug)i
t
FIG. P6.2 Accelerogram in problem P6.2
Problem P6.3 The rigid vertical column AC (Fig. P6.3) of a circular cross-
section, line density m is supported on the ground by means of a spherical hinge
at A and three elastic cables of cross-sectional area A and modulus of elasticity E.
(a) (b)
FIG. P6.3 Structure in problem P6.3
@CivilMethod
242 PART I Single-degree-of-freedom systems
The cables have been prestressed so that they can undertake compressive forces.
The column carries three advertising panels of total mass 10ma, which are
arranged as in Fig. P6.3b. Their support on the column extends to a length
1:25a. The structure is subjected to ground motion in the y direction, whose accel-
erogram u€g ðt Þ is given in Problem P6.2. Determine the minimum prestressing
force of the cable GB using the results of Problem P6.2. The cables are assumed
massless. Data: a ¼ 5m, m ¼ 0:5kNm1 s2 =m, E ¼ 2:1 108 kN=m2 , and
A ¼ 4cm . 2
Problem P6.5 The vertical columns of the frame in Fig. P6.5 have specific
weight g b and are assumed rigid. The elastic support on the ground is simulated
by the rotational springs CR . The horizontal beams are flexible with cross-
sectional moment of inertia I and modulus of elasticity E while their mass
and axial deformation are assumed negligible. The structure is subjected to
of total duration ttot ¼ 3s. Give
the horizontal ground motion ug ðt Þ ¼ u0 sin wt
the graph of the response spectrum of the rotation of the structure and compute
the extreme values of the shear force and the bending moment of the beams.
Data: Cross-sectional area of the beams A ¼ a=10 a=10, CR ¼ EI =a,
E ¼ 2:1 107 kN=m, g b ¼ 24kN=m3 , u0 ¼ 0:03m, w ¼ 3, 5 and 7s1 .
@CivilMethod
Response to ground motion and vibration isolation Chapter 6 243
=∞
(a)
h
h
(b) L
@CivilMethod
244 PART I Single-degree-of-freedom systems
[4] F.B. Hildebrand, Advanced Calculus for Applications, Prentice Hall, Englewood Cliffs, NJ,
1962.
[5] R.R. Craig Jr., J. Andrew, A.J. Kurdila, Fundamentals of Structural Dynamics, second ed., John
Wiley, New Jersey, 2006.
[6] J.W. Leonard, Tension Structures, McGraw-Hill, New York, 1988.
[7] J.T. Katsikadelis, Finite deformation of cables under 3-D loading: an analytic solution, in:
D.E. Beskos, D.L. Karabalis, A.N. Kounadis (Eds.), Proc. of the 4th National Congress on Steel
Structures, Patras, May 24–25, vol. II, 2002, pp. 526–534.
@CivilMethod
Chapter 7
Damping in structures
Chapter outline
7.1 Introduction 245 7.6.1 Introduction 257
7.2 Loss of energy due to damping 246 7.6.2 The fractional derivative 258
7.3 Equivalent viscous damping 249 7.7 Measurement of damping 260
7.4 Hysteretic damping 250 7.7.1 Free vibration decay
7.5 Coulomb damping 252 method 261
7.5.1 Free vibrations with 7.7.2 Resonance amplitude
Coulomb damping 252 method 262
7.5.2 Forced vibrations with 7.7.3 Width of response curve
Coulomb damping 255 method 263
7.6 Damping modeling via 7.8 Problems 265
fractional derivatives 257 References and further reading 267
7.1 Introduction
In this chapter, the damping of structures is discussed. Damping appears in all
mechanical systems that perform vibrations. It is the dissipation of energy in a
vibrating structure. The type of energy into which the mechanical energy is
transformed depends on the system and the physical mechanism that causes
the dissipation. The energy is lost either in the form of heat or is radiated into
the environment. For example, the loss of energy in the form of heat is perceived
when an iron rod is subjected to alternating bending. The sound produced by a
body that is hit represents the loss of energy dissipated into the environment. In
the study of vibrations, we are interested in the damping related to the response
of the structure. Damping is due to different energy dissipation mechanisms act-
ing simultaneously. In spite of the age-long detailed studies on the damping of
structures, the understanding of damping mechanisms is quite primitive. A
well-known method to get rid of this problem is to use so-called viscous damp-
ing. This approach was first introduced by Rayleigh [1] via his famous dissipa-
tion function (see Section 1.8.4).
The loss of energy of a vibrating system reduces the amplitude of the free
vibration. When a system undergoing forced vibrations reaches the phase of
the steady-state response, the loss of energy is balanced by the energy input into
the system by the excitation force.
In vibrating systems, we distinguish different types of damping forces,
which may be due to the internal molecular friction, the sliding friction, or the
where T ¼ 2p= w the period of the vibration in the steady-state response. The
quantity WD depends on several factors such as temperature, frequency, or
amplitude of the vibration.
In this section, we consider the simplest form of energy loss, that is, the loss
due to viscous damping. The damping force, in this case, is given by the relation
_ Moreover, the displacement in the steady-state phase, due to the
fD ¼ cu.
t, is given by Eq. (3.2.26), namely
harmonic force pðt Þ ¼ p0 sin w
t qÞ
u ðt Þ ¼ rsin ðw (7.2.2)
where
p0 1
r¼ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi (7.2.3)
k
2 2
1 b + ð2xbÞ2
2xb
tan q ¼ (7.2.4)
1 b2
Consequently
u_ ðt Þ ¼ r t qÞ
w cos ðw (7.2.5)
The integral (7.2.1) yields
Z T
WD ¼ cu_ 2 dt
0
Z 2p=
w
(7.2.6)
w2 r 2
¼ c cos 2 ðw
t qÞdt
0
wr2
¼ pc
@CivilMethod
Damping in structures Chapter 7 247
Of particular
pffiffiffiffiffiffiffiffiffi interest is the
pffiffiffiffiffiffienergy
ffi loss at resonance. Then it is
¼ w ¼ k=m , c ¼ 2mxw ¼ 2x km and Eq. (7.2.6) becomes
w
WD ¼ 2xpkr2 (7.2.7)
Moreover, the energy input into the system is due to the work that the force
t produces in a complete oscillation. Namely
pðt Þ ¼ p0 sin w
Z T
Wp ¼ pðt Þdu
0
Z 2p=
w
(7.2.8)
¼ t qÞ sin w
w cos ðw
p0 r tdt
0
¼ p0 rp sin q
Taking into account that (see Eq. 3.2.27)
2xb
sin q ¼ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
2 ¼ 2xbrk=p0 (7.2.9)
1 b2 + ð2xbÞ2
@CivilMethod
248 PART I Single-degree-of-freedom systems
Thus, it must be
WI + WD + W S ¼ Wp (7.2.13)
which by virtue of Eqs. (7.2.11), (7.2.12) yields Eq. (7.2.10).
Eq. (7.2.5) yields
u_ ðt Þ2 ¼ w
2 r2 cos 2 ðw
t qÞ
¼w r 1 sin 2 ðw
2 2
t qÞ (7.2.14)
¼w 2 r2 u 2
hence
fD2 ¼ c2 u_ ðt Þ2
(7.2.15)
¼ c2 w
2 r2 u 2
which is readily transformed into
2
fD 2 u
+ ¼1 (7.2.16)
wr
c r
@CivilMethod
Damping in structures Chapter 7 249
wr2 2pc
pc w 2pð2mwx Þ
w
wD ¼ ¼ ¼ ¼ 4pxb (7.2.19)
kr2 =2 mw2 mw2
¼ 2xb (7.2.20)
wr2
WD ¼ ceq p (7.3.1)
@CivilMethod
250 PART I Single-degree-of-freedom systems
@CivilMethod
Damping in structures Chapter 7 251
t qh Þ
u ðt Þ ¼ rh sin ðw (7.4.7)
where
p0 1
rh ¼ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
ffi (7.4.8)
k 2
1 b2 + 2
tan qh ¼ (7.4.9)
1 b2
Working as in Section 7.2, we can determine the geometry of the hysteresis
loop for this type of damping. Thus, we obtain from Eq. (7.4.7)
u_ ðt Þ2 ¼ r2h w
2 cos 2 ðw
t qh Þ
2 2
¼ rh w 1 sin 2 ðw t qh Þ (7.4.10)
¼w 2 r2h u 2
Moreover, it is
2
k
fD2 ¼ ðf ku Þ2 ¼ u_ 2 (7.4.11)
w
which by virtue of (7.4.10) is written as
ðf ku Þ2 ¼ 2 k 2 r2h u 2 (7.4.12)
or
2
f ku 2 u
+ ¼1 (7.4.13)
krh rh
Therefore, in the plane u, f , Eq. (7.4.13) represents a rotated ellipse
(Fig. 7.4.1). The area of the ellipse is pkr2h and expresses the loss of energy
in a cycle. We observe that Eq. (7.4.13) does not involve the excitation fre-
quency. This implies that the hysteresis loop can be determined experimentally
using a low excitation frequency, that is, quasistatic, by plotting the load-
displacement curve.
For f ¼ 0, we obtain the abscissa of the ellipse on the u axis
d ¼ rh pffiffiffiffiffiffiffiffiffiffiffiffiffi (7.4.14)
1 + 2
which can be used to evaluate the damping coefficient . Hence
d
¼ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi (7.4.15)
r2h d 2
@CivilMethod
252 PART I Single-degree-of-freedom systems
while when the body is moving to the left (Fig. 7.5.1c) the equilibrium of forces
yields the equation of motion
m u€ + ku ¼ F (7.5.4)
@CivilMethod
Damping in structures Chapter 7 253
(a)
(b) (c)
FIG. 7.5.1 Motion with Coulomb friction.
@CivilMethod
254 PART I Single-degree-of-freedom systems
Eq. (7.5.11) holds until the velocity vanishes, namely t ¼ p=w. At that
instant, the body is at the extreme left position, where the displacement is
u ðp=wÞ ¼ ðu0 2F=k Þ. The body will now start moving to the right with ini-
tial conditions
p
p
F
u ¼ u0 2 , u_ ¼0 (7.5.12)
w k w
The motion is now described by Eq. (7.5.3), which has a particular solution
F
up ðt Þ ¼ (7.5.13)
k
and general solution
F
u ðt Þ ¼ A cos wt + B sin wt (7.5.14)
k
The initial conditions (7.5.12) give
3F
A ¼ u0 , B¼0 (7.5.15)
k
and Eq. (7.5.14) becomes
F F
u ðt Þ ¼ u 0 3 cos wt (7.5.16)
k k
@CivilMethod
Damping in structures Chapter 7 255
Eq. (7.5.16) holds until the body reaches the extreme right position, namely
until the instant t ¼ 2p=w. At that time, the body has completed a full oscillation
and the displacement is u ð2p=wÞ ¼ ðu0 4F=k Þ. This solution procedure con-
tinues to obtain the response of the next oscillations. The graphical representa-
tion of the displacement versus time is shown in Fig. 7.5.2. The curve was
obtained by numerical integration of the equation of motion (7.5.6) with
m ¼ 10kNm1 s2 , T ¼ 0:5s, F ¼ 23:685kN, u0 ¼ 0:55m, u_ 0 ¼ using the pro-
gram aem_nlin.m developed in Chapter 5. The numerical results coincide with
those obtained by the above-presented analytical solution. The motion is a
vibration with a period T ¼ 2p=w, which means that the Coulomb friction does
not affect the frequency or the period of vibration. The amplitude of vibration is
reduced in each cycle by 4F=k. A consequence of this is that the envelopes of
the curve are straight lines, unlike in the cases of viscous or hysteretic damping
where the envelopes are exponential functions. The motion of the system con-
tinues until the elastic force ku becomes smaller than the force F of the friction.
Until now, no difference was made between static friction Fs ¼ ms N and
dynamic friction Fd ¼ md N . The first occurs when the body is stationary and
the second when the body moves. Generally, it is md < ms , hence the dynamic
friction coefficient md will be used in the equation of motion while the static
friction coefficient ms is used for the control of the motion. In viscous or hys-
teretic damping, theoretically, the body does not stop moving because the
amplitude of the vibration reduces exponentially. Nevertheless, real structural
systems stop after a finite time. This is due to the fact that the Coulomb friction
coexists with other forms of damping and forces the moving systems to stop.
@CivilMethod
256 PART I Single-degree-of-freedom systems
WD ¼ 4Fr (7.5.23)
We can determine an equivalent coefficient of viscous damping by equating
the loss of energy (Eq. 7.5.23) with that of the viscous damping given by
Eq. (7.3.1). Namely,
wr2 ¼ 4Fr
ceq p (7.5.24)
which gives
4F
ceq ¼ (7.5.25)
p
wr
and an equivalent damping ratio
ceq 2F
xeq ¼ ¼ (7.5.26)
2mw pkrb
@CivilMethod
Damping in structures Chapter 7 257
Eq. (7.5.28) holds if 1 ð4F=pp0 Þ2 > 0, that is, when F=p0 < p=4. Obvi-
ously, for F=p0 > p=4 r becomes imaginary and this method for determining
an equivalent damping coefficient does not apply.
The phase angle results from Eq. (7.2.4) by setting x ¼ xeq and taking the
value of r from Eq. (7.5.28). Thus, we have
4F=pp0
tan q ¼ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi (7.5.29)
1 ð4F=pp0 Þ2
The plus sign is valid when b < 1 while the minus sign is valid when b > 1.