TMD Intro
TMD Intro
TMD Intro
Tuned Mass
Damper
Systems
4.1
INTRODUCTION
A tuned mass damper (TMD) is a device consisting of a mass, a spring, and a
damper that is attached to a structure in order to reduce the dynamic
response of the structure. The frequency of the damper is tuned to a particular
structural frequency so that when that frequency is excited, the damper will
resonate out of phase with the structural motion. Energy is dissipated by the
damper inertia force acting on the structure. The TMD concept was first
applied by Frahm in 1909 (Frahm, 1909) to reduce the rolling motion of ships as
well as ship hull vibrations. A theory for the TMD was presented later in the
paper by Ormondroyd and Den Hartog (1928), followed by a detailed
discussion of optimal tuning and damping parameters in Den Hartogs book on
mechanical vibrations (1940). The initial theory was applicable for an undamped
SDOF system subjected to a sinusoidal force excitation. Extension of the theory
to damped SDOF systems has been investigated by numerous researchers.
Significant contributions were made by Randall et al. (1981), Warburton (1981,
1982), Warburton and Ayorinde (1980), and Tsai and Lin (1993).
This chapter starts with an introductory example of a TMD design and a
brief description of some of the implementations of tuned mass dampers in
building structures. A rigorous theory of tuned mass dampers for SDOF systems
subjected to harmonic force excitation and harmonic ground motion is discussed
next. Vari- ous cases, including an undamped TMD attached to an undamped
SDOF system, a damped TMD attached to an undamped SDOF system, and
a damped TMD attached to a damped SDOF system, are considered. Time
history responses for a
217
218 Chapter 4
range of SDOF systems connected to optimally tuned TMD and subjected to harmonic and seismic excitations are presented. The theory is then extended to MDOF
systems, where the TMD is used to dampen out the vibrations of a specific mode.
An assessment of the optimal placement locations of TMDs in building structures is
included. Numerous examples are provided to illustrate the level of control that can
be achieved with such passive devices for both harmonic and seismic excitations.
4.2
AN INTRODUCTORY EXAMPLE
In this section, the concept of the tuned mass damper is illustrated using the twomass system shown in Figure 4.1. Here, the subscript d refers to the tuned mass
damper; the structure is idealized as a single degree of freedom system.
Introducing the following notation
2
= k
-m
(4.1)
c = 2m
(4.2)
kd
2
d = ------md
(4.3)
cd = 2 d d m d
(4.4)
(4.5)
= - ---- m
u
Primary mass 1 + mu + 2u + u
m
kd
md
m
cd
c
u
(4.6)
FIGURE 4.1:
system.
SDOF-TMD
u + ud
Section 4.2
Tuned mass ud + 2 d d u d +
d u d
=
u
(4.7)
The purpose of adding the mass damper is to limit the motion of the structure
when it is subjected to a particular excitation. The design of the mass damper
involves specifying the mass md ,
kd, and damping coefficient cd . The
stiffness
optimal choice of these quantities is discussed in Section 4.4. In this example, the
near-optimal approximation for the frequency of the damper,
d =
(4.8)
is used to illustrate the design procedure. The stiffnesses for this frequency
combi- nation are related by
(4.9)
k d = mk
Equation (4.8) corresponds to tuning the damper to the fundamental period of the
structure.
Considering a periodic excitation,
p = p sint
(4.10)
u sin t + 1
u = u d sint + 1 +
2
ud =
(4.11)
(4.12)
where u and denote the displacement amplitude and phase shift, respectively.
The critical loading scenario is the resonant condition, = . The solution for this
case has the following form:
u =
--p
--
-k
m
------------------1
-------------------2
1 +1 ---- -- +
--- ----- 2
--1
----= u
2d
(4.13)
2d
(4.14)
2
1
tan1 = ------ +
-------m 2d
tan2 =
--2
(4.15)
(4.16)
220 Chapter 4
Note that the response of the tuned mass is 90 out of phase with the response of
the primary mass. This difference in phase produces the energy dissipation
contributed by the damper inertia force.
The response for no damper is given by
u =
p
-- ------
k 2
(4.17)
= ---
1
2
(4.18)
2
1
1 + ---- -- + ---- ---- 2
where
=
m
-----
(4.20)
2 d
Equation (4.20) shows the relative contribution of the damper parameters to the
total damping. Increasing the mass ratio magnifies the damping. However, since
the added mass also increases, there is a practical limit on m. Decreasing the
damping coefficient for the damper also increases the damping. Noting Eq.
(4.14), the rela- tive displacement also increases in this case, and just as for the
mass, there is a prac- tical limit on the relative motion of the damper. Selecting
the final design requires a compromise between these two constraints.
0 and we want to add a tuned mass damper such that the equivalent
damping ratio is 10%. Using Eq. (4.20), and setting e = 0.1, the following relation
between m and d is obtained:
m
2
1
----- 1 + ---- -- + --- ----- 2 =
2
2d
0.1
(1)
--1
---- u
2d
(2)
Section 4.2
(3)
u d
-----
-- 2
1+
---
2
u
10
u
------------2ue d u
m
Finally, taking u
=
(5)
2 0 . 1
m = - --- -- --- -- --- - = 0.02
10
(6)
(7)
kd = mk = 0.02
k
(8)
A description of some applications of tuned mass dampers to building structures is presented in the following section to provide additional background on
this type of device prior to entering into a detailed discussion of the underlying
theory.
222 Chapter 4
4.3
Floor beam
Direction of motion
Section 4.3
224 Chapter 4
Chiba Port Tower (completed in 1986) was the first tower in Japan to be equipped
with a TMD. Chiba Port Tower is a steel structure 125 m high weighing 1950
metric tons and having a rhombus-shaped plan with a side length of 15 m. The
first and second mode periods are 2.25 s and 0.51 s, respectively for the x
direction and 2.7 s and 0.57 s for the y direction. Damping for the fundamental
mode is estimated at 0.5%. Damping ratios proportional to frequencies were
assumed for the higher modes in the analysis. The purpose of the TMD is to
increase damping of the first mode for both the x and y directions. Figure 4.3
shows the damper system. Manu- factured by Mitsubishi Steel Manufacturing
Co., Ltd., the damper has mass ratios with respect to the modal mass of the first
mode of about 1/120 in the x direction and 1/80 in the y direction; periods in the
x and y directions of 2.24 s and 2.72 s, respectively; and a damper damping ratio
of 15%. The maximum relative displacement of the damper with respect to the tower is about 1 m in each direction.
Reductions of around 30 to 40% in the displacement of the top floor and 30% in
the peak bending moments are expected.
The early versions of TMDs employ complex mechanisms for the bearing and
damping elements, have relatively large masses, occupy considerable space, and are
quite expensive. Recent versions, such as the scheme shown in Figure 4.4, have
been designed to minimize these limitations. This scheme employs a multiassemblage of elastomeric rubber bearings, which function as shear springs, and bitumen
rubber compound (BRC) elements, which provide viscoelastic damping capability.
The device is compact in size, requires unsophisticated controls, is multidirectional,
and is easily assembled and modified. Figure 4.5 shows a full-scale damper being
subjected to dynamic excitation by a shaking table. An actual installation is contained in Figure 4.6.
FIGURE 4.3: Tuned mass damper for Chiba-Port Tower. (Courtesy of J. Connor.)
Section 4.3
BRC
Multi-stage
rubber bearings
Air brake
Limit switch
FIGURE 4.4: Tuned mass damper with spring and damper assemblage.
FIGURE 4.6: Tuned mass damperHuis Ten Bosch Tower, Nagasaki. (Courtesy
of J. Connor.)
226 Chapter 4
Floor beam
Direction of motion
L
sinL
T W d
Introducing these approximations transforms Eq. (4.21) to
(4.22)
Section 4.3
L
keq
8
md
ud
md
t=0
u + ud
W
md ud + -----L
= m u
d
(4.23)
W
k eq = ------d
-L
(4.25)
-L
(4.26)
--g
The simple pendulum tuned mass damper concept has a serious limitation.
Since the period depends on L, the required length for large Td may be greater than
the typical story height. For instance, the length for Td = 5 s is 6.2 meters whereas
the
story height is between 4 and 5 meters. This problem can be eliminated by resorting to
the scheme illustrated in Figure 4.9. The interior rigid link magnifies the
support motion for the pendulum and results in the following equilibrium
equation:
228 Chapter 4
md
u + u1
u + u1 + u d
pendulum.
W
d
+ ------ -- = 0
d
m d u + u1 + u d u
L
(4.27)
The rigid link moves in phase with the damper and has the same displacement
amplitude. Then, taking u1 = ud in Eq. (4.27) results in
m
Wd
mdu d + -------d =
2L
-----u
u
2
(4.28)
The equivalent stiffness is Wd/2L, and it follows that the effective length is equal to
2L. Each additional link increases the effective length by L. An example of a pendulum-type damper is described next.
Crystal Tower (Nagase and Hisatoku, 1990)
The tower, located in Osaka, Japan, is 157 m high and 28 m by 67 m in plan,
weighs 44000 metric tons, and has a fundamental period of approximately 4 s in
the north- south direction and 3 s in the east-west direction. A tuned pendulum
mass damper was included in the early phase of the design to decrease the windinduced motion of the building by about 50%. Six of the nine air cooling and
heating ice thermal storage tanks (each weighing 90 tons) are hung from the top
roof girders and used as a pendulum mass. Four tanks have a pendulum length
of 4 m and slide in the north-south direction; the other two tanks have a
pendulum length of about 3 m and slide in the east-west direction. Oil dampers
connected to the pendulums dissi- pate the pendulum energy. Figure 4.10 shows
the layout of the ice storage tanks that were used as damper masses. Views of the
actual building and one of the tanks are presented in Figure 4.11 on page 230.
The cost of this tuned mass damper system was around $350,000, less than 0.2%
of the construction cost.
Section 4.4
27.6 m
9.0 m
x
x
4 m67.2 m
Support frame
Suspension material
Ice storage
tank
Coil spring
Oil
damper
Stopper
Oil damper
Coil spring
230 Chapter 4
(4.29)
ma g +
mu + ku k du d =
p
(4.30)
where ag is the absolute ground acceleration and p is the force loading applied to
the primary mass.
Section 4.4
8
R
md
Floor
(a)
md
keq
(b)
kd
md
m
u + ug
ug
ud + u + ug
a g sin
t
(4.31)
232 Chapter 4
p = p sin t
(4.32)
u sin
t u
(4.33)
(4.34)
sin t
and substituting for these variables, the equilibrium equations are transformed to
md 2 + kd u d m d 2 u =
m d a
k d u d + m 2 + ku =
The solutions for u and u
by
(4.35)
g
ma g +
p
are given
(4.36)
1
p
d
u =
---- -
ma 1 + m 2
g
-------- d
----------------------- ---- ---- -----
k D1
k
D1
= --p--- -m---------
m
a g
m
------ --- - -- -----
k d D1
k d D 1
(4.37)
(4.38)
where
2
D1 = 1 1 d m
(4.39)
= =
d =
--- ---------------
k m
- =
- ------------------- ---
d
kd
md
(4.40)
(4.41)
Selecting the mass ratio and damper frequency ratio such that
2
1 +m = 0
d
(4.42)
Section 4.4
u =
-p
k
(4.43)
p
2
= -- --- +
m
a g
------ --- kd
kd
(4.44)
This choice isolates the primary mass from ground motion and reduces the
response due to external force to the pseudostatic value, p k. A typical range
for m is 0.01 to 0.1. Then the optimal damper frequency is very close to the
forcing frequency. The exact relationship follows from Eq. (4.42).
opt
------
----------
(4.45)
1+
m
We determine the corresponding damper stiffness with
kd
=
d
opt
m
opt
d
----2m
----m
= -----1- + m
(4.46)
a
= 1-----+----- --m
---- p--- g--
+
m k
2
(4.47)
The next level of complexity has damping included in the mass damper, as shown in
Figure 4.14. The equations of motion for this case are
m d u d + c d u d + k d u d + md u
=
md
ag
mu + ku cdud kdud = ma g +
p
(4.48)
(4.49)
234 Chapter 4
The inclusion of the damping terms in Eqs. (4.48) and (4.49) produces a phase
shift between the periodic excitation and the response. It is convenient to work
initially with the solution expressed in terms of complex quantities. We express
the excitation as
ag =
p = ge
(4.50)
it
(4.51)
p
e
where a
it
it
(4.52)
i
t
(4.53)
md
a g
ma
+ p
(4.54)
(4.55)
g
---kD 2
a g
2
---
m ------- 1 + mf
i2f1 + m
d
(4.56)
2
kD2
kd
k
md
m
cd
ug
u + ug
ud + u + ug
FIGURE 4.14: Undamped SDOF system coupled with a damped TMD system.
Section 4.4
2
ud = -p---------a g
- --- -m
----kD2 kD2
(4.57)
where
2
2 2
2
D2 = 1 f
m f
d
f = -----
+ i2f
1 1 + m
d
(4.58)
(4.59)
i
u = p H e 1 -a---g-m
----H
i 2
e
k
1
i
ud = p--H e
k
a
m
2
i 3
---- -----H e
(4.61)
where the H factors define the amplification of the pseudo-static responses, and the
s are the phase angles between the response and the excitation. The various H
and terms are as follows:
2
2 2
f + 2f d
H 1 = --------------------------------------------- -----------D2
2
2 2
1 + mf + 2f1 +
H 2 = m 2
d----------------------------------------------------------------------------------------------D2
H3 =
2
--- ------D2
H4 =
1
---- -----D2
(4.62)
(4.63)
(4.64)
(4.65)
D2 =
Also,
1
f
2 2 2
m f
1 = 1 3
+ 2 d f 1 1 +
m
(4.66)
(4.67)
236 Chapter 4
2 = 2 3
(4.68)
2
2f 1 1 + m
tan 3 = d-------------------------------------------------------------- 2
2
2
2
1 f m f
2
tan1 =
2 d f
---------------tan2 = f 2 2
(4.69)
(4.70)
(4.71)
2f1 + m
----------------------------------d
2
2
1 + mf
For most applications, the mass ratio is less than about 0.05. Then the amplification
factors for external loading H 1 and ground motion H 2 are essentially equal. A
similar conclusion applies for the phase shift. In what follows, the solution corresponding to ground motion is examined and the optimal values of the damper properties for this loading condition are established. An in-depth treatment of the
external forcing case is contained in Den Hartogs text (Den Hartog, 1940).
Figure 4.15 shows the variation of H2 with forcing frequency for specific values of damper mass m and frequency ratio f , and various values of the damper
damping ratio, d .
d = 0, there are two peaks with infinite amplitude
When
located on each side of = 1 . As d is increased, the peaks approach each other
and then merge into a single peak located at 1. The behavior of the
amplitudes suggests that there is an optimal value of d for a given damper
configuration ( md
and kd or, equivalently, m and f ). Another key observation is that all the curves
pass through two common points, P and Q. Since these curves correspond to different values of d , the location of P and Q must depend only on m and f .
Proceeding with this line of reasoning, the expression for H2 can be written as
2
H=
2
2
2 2
a a a +
a1 + d a 2
2
1
2
d
-------------------- =
-----------------------------2
2
-2
2
2 2
a
4 a3 a4 + d
a3 + da 4
(4.72)
H2
PQ
a2
---a
4
(4.74)
Section 4.4
30
m=
f=1
0.01
d
=1
=0
25
H2
20
15
10
0
0.8
0.85
0.9
0.95
p
1
1
D
p =w
1.05
1.1
1.15
Substituting for the a terms in Eq. (4.73), we obtain a quadratic equation for :
4
1 + mf 2 + 1----+ 0 .5 m
----- --- - ---- ---1+m
+f
0
(4.75)
The two positive roots 1 and 2 are the frequency ratios corresponding to points
P and Q. Similarly, Eq. (4.74) expands to
2 PQ
= --------------1----+----m
-----------------2
1 1 1 + m
(4.76)
238 Chapter 4
30
25
20
>
H2
15
d
d
Q
dlopt
10
0.85
p
1l opt
0.95
0.9
1
p
p
2lo
1.05
D
=w
1.1
1.15
pt
Requiring the amplitudes to be equal at P and Q is equivalent to the following condition on the roots:
2
1 1 + m=
1
(4.77)
1 1 + m
2
Then, substituting for 1 and 2 using Eq. (4.75), we obtain a relation between the
optimal tuning frequency and the mass ratio:
(4.78)
fopt =
d opt
----
---- ----0 .5
--- - ---m
---1+m
1
= f
opt
(4.79)
Section 4.4
12
opt
H2
opt
1 0.5m
------------------------1+m
(4.80)
(4.81)
=
1 +
-- ---- ----m
-----0.5m
The expression for the optimal damping at the optimal tuning frequency is
opt
m3 0.5m
-----------------------------------------------81 + m1 0.5m
(4.82)
Figures 4.17 through 4.20 show the variation of the optimal parameters with the
mass ratio, m.
The response of the damper is defined by Eq. (4.61). Specializing this equation for the optimal conditions leads to the plot of amplification versus mass ratio
contained in Figure 4.21. A comparison of the damper motion with respect to the
motion of the primary mass for optimal conditions is shown in Figure 4.22.
0.98
fopt
0.96
0.94
0.92
0.9
0.88 0
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
240 Chapter 4
1.1
fopt
1.05
0.95
1, 2lopt
p2 lopt
1 lopt
0.9
0.85
0.8 0
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
dlopt
0.12
0.1
0.08
0.06
0.04
0.02
0
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
Section 4.4
25
15
H2 l
opt
20
10
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
FIGURE 4.20: Maximum dynamic amplification factor for SDOF system (optimal
tuning and damping).
300
250
H4 l
opt
200
150
100
50
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
242 Chapter 4
20
fopt, popt,
18
dlopt
16
H
= 4
H2
10
u
12
u d
14
8
6
4
2
0
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
are plotted in Figure 4.23 and Figure 4.24. The response for no damper is also plotted in Figure 4.23. We observe that the effect of the damper is to limit the motion in
a frequency range centered on the natural frequency of the primary mass and
extending about 0.15. Outside of this range, the motion is not significantly influenced by the damper. d= 0.0
= 0.03
= 0.061 (op
= 0.1
d
timal)
damper
d
m = 0.01
fopt = 0.98
76
30
No
25
20
H2
15
10
0.8
0.85
0.9
0.95 p
1
1lopt
p D
=w
pt
p 1.05
2 lo
1.1
1.15
1.2
FIGURE 4.23: Response curves for amplitude of system with optimally tuned TMD.
Section 4.4
150
=
d=
d=
d=
d=
d
m=0
fopt =
0.0
0.03
0.06
0.1 1 (optimal)
0.2
.01
0.9876
H4
100
50
0.85
0.9
0.95
D
=w
1.05
1.1
1.15
(4.83)
1
H -- ---2
-----------------2H 2
1
(4.85)
opt
Figure 4.25 shows the variation of e with the mass ratio. A mass ratio of 0.02 is
equivalent to about 5% damping in the primary system.
244 Chapter 4
0.12
fopt, popt,
lopt
0.1
0.08
0.06
0.04
0.02
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
opt
(4.86)
Compute kd :
kd =
2
opt
(4.87)
d
d
Determine d
Compute cd :
mk
2
m = f
cd = 2 d
opt
dmd
= mf opt
2 d
opt
(4.88)
Section 4.4
2 opt
(1)
H4
---- ---------- 6
H-2
opt
(2)
Constraint Eq. (1) requires m 0.05. For constraint Eq. (2), we need to take
m 0.02 . Therefore, m 0.05 controls the design. The relevant parameters are
m = 0.05
f opt = 0.94
md = 0.05m
d = 0.94
kd = m f opt
k
opt
= 0.135
Then
2
= 0.044k
=
md u d + c d ud + k dud +
mda g
mdu
mu + cu + ku c d u d kd u d mag +
=
p
(4.89)
(4.90)
p
kd
md
m
c
cd
ug
u+
ug
FIGURE 4.26:
ud + u + ug
246 Chapter 4
Proceeding in the same way as for the undamped case, the solution due to
periodic excitation (both p and ug) is expressed in polar form:
i 5 a g m
i
p
= ---H
----------H 6 e 6
5 e
u
k
k
i 7 a g m
i
= p---H
----------H 8 e 8
7 e
ud
k
k
(4.91)
(4.92)
2 2
(4.93)
f + 2 f d
H5 = --------------------------------------------- -------------D3
2
2 2
(4.94)
1 + mf + 2f1 +
H6 = m 2
----------------------------------------------------------------d
-------------------------------D3
H7 =
----- --2-D-3
8 =
1+ 2
---- -- ---- --- ---- 2
H
----------
D3 =
(4.95)
(4.96)
D3
f 2 2 m + 1 2 f 2 2 4 d f 2 2
(4.97)
+ 4f 2 2 + d f1 2 1 + m
5 = 1 7
(4.98)
6 = 2 7
(4.99)
8 = 3 7
(4.100)
tan 7 =
f 2 2 + f 1 2 1 + m
d
2-------------------------------------------------------------------------------------------------- f 2 2 m + 1 2 f 2 2 4 d f
2
tan3 = 2
(4.101)
(4.102)
this case, these parameters also depend on . Numerical simulations can be applied
Section 4.4
247
to evaluate H5 and H7 for a range of , given the values for m, , f , and d . Starting with specific values for
and , plots of H5 versus can be generated for a
m
range of f and d . Each H5 plot has a peak value of H5 . The particular
combi- nation of f and d that corresponds to the lowest peak value of H5 is
taken as the
optimal state. Repeating this process for different values of m and produces the
behavioral data needed to design the damper system.
Figure 4.27 shows the variation of the maximum value of H5 for the optimal
state. The corresponding response of the damper is plotted in Figure 4.28. Adding
damping to the primary mass has an appreciable effect for small m. Noting
Eqs. (4.91) and (4.92), the ratio of damper displacement to primary mass displacement is given by
-
- --
---u
H7
(4.103)
Since is small, this ratio is essentially independent of . Figure 4.29 confirms this
statement. The optimal values of the frequency and damping ratios generated
through simulation are plotted in Figures 4.30 and 4.31. Lastly, using Eq. (4.85),
H5
can be converted to an equivalent damping ratio for the primary system.
opt
e =
-----------------2H 5
(4.104)
opt
1 0 .5 m
2
f = ----- ---- ----- --- - ---- ---- + 1 2
1
1+m
(4.105)
d =
------------------------------------------------- + 0.151
0.170
81 + m1 0.5m
2
+ 0.163+ 4.980 m
(4.106)
248 Chapter 4
40
=
=
=
=
35
30
0.0
0.01
0.02
0.05
0.1
H5l
opt
25
20
15
10
5
0
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
FIGURE 4.27: Maximum dynamic amplification factor for damped SDOF system.
300
=0.0
=0.01
=
=0.02
=0.05
0.1
250
H7l
opt
200
150
100
50
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
Section 4.4
20
= 0.0
= 0.01
= 0.02
18
16
= 0.05
= 0.1
12
10
u d
H7
=
uH5
14
8
6
4
2
0
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
0.98
0.96
0.94
0.0
0.01
0.02
0.05
0.1
fopt
0.92
0.9
0.88
0.86
0.84
0.82
0.8
0.01
0.02
FIGURE 4.30:
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
250 Chapter 4
0.2
0.18
0.16
0.14
d opt
0.12
0.1
0.08
= 0.0
= 0.01
= 0.02
= 0.05
= 0.1
0.06
0.04
0.02
0
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
0.12
0.1
0.08
= 0.0
= 0.01
= 0.02
= 0.05
= 0.1
0.06
0.04
0.02
0
0.01
0.02
0.03
0.04
0.05
m
0.06
0.07
0.08
0.09
0.1
Section 4.5
251
5 opt
(1)
H7
---- ---------- 6
H--5
opt
(2)
Using Figure 4.27, the required mass ratio for = 1.2 is m 0.03. The other optimal values are fopt = 0.965 and
= 0.105 . Then
opt
d
md = 0.03m
d = 0.955
kd =
mf opt
0.027k
In this case, there is a significant reduction in the damper mass required for this
set of motion constraints. The choice between including damping in the primary
system versus incorporating a tuned mass damper depends on the relative costs
and reli- ability of the two alternatives, and the nature of the structural
problem. A TMD system is generally more appropriate for upgrading an
existing structure where access to the structural elements is difficult.
4.5
252 Chapter 4
T = 0.49 s
= 0.02
m = 0.01
1
0.8
0.6
u (m)
0.4
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
0
10
20
With TMDWithout
TMD 30
Time (s)
40
50
60
ud (m)
-2
-4
-6
10
20
30
40
50
60
Time (s)
FIGURE 4.34:
Section 4.5
253
T = 0.49 s
= 0.02
m = 0.01
1
0.8
0.6
0.4
u (m)
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
Without TMD
10
With TMD 20
30
40
50
60
Time (s)
ud (m)
-2
-4
-6
FIGURE 4.36:
10
20
30
Time (s)
40
50
60
254 Chapter 4
1
0.8
0.6
0.4
u (m)
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
Without TMD
10
With TMD 20
30
40
50
60
Time (s)
ud (m)
-2
-4
-6
10
FIGURE 4.38:
20
30
Time (s)
40
50
60
Section 4.5
255
T = 5.35 s
= 0.02
m = 0.01
1
0.8
0.6
0.4
u (m)
0.2
0
-0.2
-0.4
-0.6
-0.8
Without TMD
With TMD
-1
0
10
20
30
40
50
60
Time (s)
ud (m)
-2
-4
-6
FIGURE 4.40:
10
20
30
Time (s)
40
50
60
256 Chapter 4
T = 5.35 s
= 0.02
m = 0.01
1
0.8
0.6
0.4
u (m)
0.2
0
-0.2
-0.4
-0.6
-0.8
-1
Without TMD
10
With TMD 20
30
40
50
60
Time (s)
ud (m)
-2
-4
-6
FIGURE 4.42:
10
20
30
Time (s)
40
50
60
Section 4.5
257
T = 5.35 s
= 0.02
m = 0.01
1
0.8
0.6
0.4
u (m)
0.2
0
-0.2
-0.4
-0.6
-0.8
Without TMD
With TMD
-1
0
10
20
30
40
50
60
Time (s)
ud (m)
-2
-4
-6
10
FIGURE 4.44:
20
30
Time (s)
40
50
60
258 Chapter 4
4.6
p2
k1
k2
kd
m1
m2
c1
md
c2
cd
u1 + ug
ug
u2 + ug
u2 + ug + ud
m1u1 + c 1u 1 + k1 u1 k2 u 2 u1 c 2 u 2 u1
p1 m 1 u g
=
m2 u2 + c 2 u2 u1 + k 2 u 2 u1 kdu d
= p2
c du d
m2ug
md ud + kdud + cd ud = md u 2 + ug
(4.107)
(4.108)
(4.109)
The key step is to combine Eqs. (4.107) and (4.108) and express the
resulting equation in a form similar to the SDOF case defined by Eq. (4.90). This
operation reduces the problem to an equivalent SDOF system, for which the
theory of Section 4.4 is applicable. The approach followed here is based on
transforming the original matrix equation to scalar modal equations.
Introducing matrix notation, Eqs. (4.107) and (4.108) are written as
0
p1 m 1 a g
MU + CU + KU =
p2
m2 a g
kdud +
c d ud
(4.110)
Section 4.6
U =
u1
(4.111)
u2
M=
m1
(4.112)
m2
K =
k2
C =
(4.113)
k1 + k2 k 2
k2
(4.114)
c1 + c2 c2
c2
c2
(4.115)
The modal vectors satisfy the following orthogonality relations [see Eq. (2.211)]:
T
(4.116)
j Ki = ij j j
Mi
j
Mj
T
(4.117)
2
= j
m j
= j
K j
c j = jTC j
(4.118)
(4.119)
j =
j1
j2
(4.120)
= K
(4.121)
260 Chapter 4
m j q j + c j q j +
kj q j =
j1 p 1 m1 a g
j = 1, 2
+ j2 p 2 m2 a g + kd u d + c du d
(4.122)
= j
--------2
(4.123)
m 1 q1 + c 1 q 1 + k1 q 1 = 11p1 + 12p2
(4.124)
m1 11 + m2 12 a g + 12 kd u d + cd ud
In general, u2 is obtained by superposing the modal contributions
u2 = 12q1 + 22 q2
(4.125)
However, when the external forcing frequency is close to , the first mode
1
response will dominate, and it is reasonable to assume
u2 12 q1
(4.126)
Solving for q1
q1 =
---1
----u
12
(4.127)
Section 4.6
1e
u 2 + c
1e
u2 + k1e u 2 = kd u d + cd ud
1e
(4.128)
1e m1e ag
1e
k1e =
(4.129)
--- m
1
------ 1
2
12
(4.130)
---- k 1
----2
12
1e
=
k1e
(4.131)
(4.132)
p
1 1 + 1 p
p = -----1------- ----------- ---2
2
1e
- --
--12- 12
1e = --------- m 1 11 +
m2 22
m 1
(4.133)
Equations (4.109) and (4.128) are similar in form to the SDOF equations
treated in the previous section. Both set of equations are compared next.
TMD equation
md ud + cd ud + kd ud = m d u a g
versus
(4.134)
md ud + cd ud + kd ud = m d u 2 ag
Primary mass equation
mu + cu + ku = cd u + kd ud + p mag
d
versus
m 1e u2 + c 1e u 2 +
k1e u 2
Taking
= c d u d + kd u d + p
u2 u
(4.135)
1e
1e m
1e a g
c 1e
c
1e
1e
1e
k 1e
k
(4.136)
262 Chapter 4
mu + cu + ku = cd u + k ud + p ma g
d
d
(4.137)
which differs from the corresponding SDOF equation by the factor . Therefore,
the solution for ground excitation generated earlier has to be modified to account
for the presence of .
The generalized solution is written in the same form as the SDOF case. We
need only to modify the terms associated with ag (i.e., H6, H8 and Their
expanded form is as follows:
2 2
H6 =
+ mf + 2f+ m
-------------------------------------------------------------------d
--------------------------------
(4.138)
D3
H8 =
(4.139)
2
2
1 + 1 + 2
-----------------------------------------------------------------D3
tan a2
tan a3
2f+
m
d
= ------------------------------------2
2
f +
-- m
2
= ------------ --- ---- ------------1 + 1 2
(4.140)
(4.141)
6 = a2 7
(4.142)
8 = a3 7
(4.143)
(4.144)
Given m and 1 , we find the tuning frequency and damper damping ratio using
Figures 4.30 and 4.31. The damper parameters are determined with
m d = m m
1e
d = f
opt
cd = 2
d
(4.145)
o
p
t
dm d
(4.146)
Section 4.6
1e
mM
1
1
----------------------------2
12
(4.148)
shows that we should select the TMD location to coincide with the maximum
amplitude of the mode shape that is being controlled. In this case, the first mode
is the target mode, and 12 is the maximum amplitude for 1.
This derivation can be readily generalized to allow for tuning on the i th
modal frequency. We write Eq. (4.127) as
qi ---1----u
i2 2
(4.149)
ie
k ie
ie
(4.150)
(4.151)
= m
ie
Given m
1
--- ----m i
2
i2
(4.152)
= 118.44
8 =
78.96
Requiring a fundamental mode damping ratio of 2%, and taking damping proportional to stiffness (C = K), the corresponding is
2
= 0.0064
= --------1 =
1 0 .
--- -0 2
--- -
264 Chapter 4
The mass, stiffness, and damping matrices for these design conditions are
M = 10
01
K = 197.39 78.96
78.96 78.96
C =
1.26 0.51
0.51 0.51
2 = 15.39 rad/s
1 = 0.5
1.0
2 =
1.0
0.5
1 M 1 = 1.25
k1 = 1 K 1 = 49.35
= 2T M2 = 1.25
k2 = 2
K 2
= 296.09
1 C1 = 0.32
c
= 2 C 2 = 1.90
2 =
= 0.02
1 = ------ 1
------2-1 m
c
------ -
= 0.049
------2-2 m
2
The optimal parameters for a TMD located at node 2, having a mass ratio of
0.01 and tuned to a specific mode, are as follows:
Mode 1: optimum location is node 2
fopt = 0.982
md = 0.0125
opt
= 0.062
kd = 0.4754
cd = 0.0096
md = 0.972
0.0125
opt
k d = 0.068
2.7974
cd = 0.0254
Section 4.6
265
This result is for the damper located at node 2. When located at node 1, the mass
and stiffness are reduced 75%.
The general case of a MDOF system with a tuned mass damper connected
to the nth degree of freedom is treated in a similar manner. Using the notation
defined previously, the jth modal equation can be expressed as
m j q j + c j qj + kj q j = p j + jn k d ud +
c d ud
where p
j = 1, 2, . . .
(4.153)
denotes the modal force due to ground motion and external forcing, and
qi ---1----u
in
(4.154)
m ie u n + c ie u n + kie
un =
where
ie
k ie
p ie + k d ud +
c d ud
i
--- M
----- =
2
in
---
----2
in
i M
= 2i
m
(4.155)
(4.156)
(4.157)
ie
ie
= kie
p = --ie
1
----p i
in
The remaining steps are the same as described previously. We specify m
(4.158)
(4.159)
and i ,
determine the optimal tuning and damping values with Figures 4.30 and 4.31, and
then compute md and d .
=
md
m m
=
ie
T
m
----- M
--i
2
i
(4.160)
in
d = f
opt i
(4.161)
The optimal mass damper for mode i is obtained by selecting n such that in is the
maximum element in i .
266 Chapter 4
y, u
P*
x*
EI constant
x
L
FIGURE E4.5a
Consider the simply supported beam shown in Figure E4.5a. The modal
shapes and frequencies for the case where the cross sectional properties are constant and the transverse shear deformation is negligible are
n x= sin n--- x
-L
E I n
= ----- -- ---- -4
n
m L
(1)
(2)
n = 12
We obtain a set of N equations in terms of N modal coordinates by
expressing the transverse displacement, u(x, t), as
(3)
u =
q i t j x
j=1
and substituting for u in the principle of virtual displacements specialized for negligible transverse shear deformation [see Eq. (2.157)],
L
M dx
(4)
b u
dx
Substituting for ,
dx2
(5)
Section 4.6
and taking
u = q j j
(6)
(7)
j = 12 . . . N
Lastly, we substitute for M and b in terms of and q and evaluate the integrals. The expressions for M and b are
(8)
M = EI = EI
qllxx
l=1
N
b = mu + bxt = m
lql + bxt
(9)
l=1
(10)
L
j k dx =
---0
jk
2
j
-- L
-- jxxkxxdx 4 ---jk 2
=
0
L
(11)
L m
= ---------2
-j
L--- 4
----
(12)
j
(13)
kj = EI
L 2
jx
= b sin ----
--0
dx L
(14)
(15)
268 Chapter 4
When the external loading consists of a concentrated force applied at the location x = x* (see Figure E4.5a), the corresponding modal load for the jth mode is
p
(16)
*
* sin
Pj
x
----------L
In this example, the force is considered to be due to a mass attached to the beam
as indicated in Figure E4.5b. The equations for the tuned mass and the force are
*
m d u + ud + kd u d + cd ud =
0
*
m d
u
+ ud = P*
(17)
(18)
u* + ud
md
kd
cd
u*
x*
FIGURE E4.5b
Suppose we want to control the ith modal response with a tuned mass damper
attached at x = x*. Taking j equal to i in Eqs. (12) and (13), the ith modal equation
has the form
m i q i + k i q i
sin
k du d + c du d
x
---- - *
----L
(19)
(20)
ie u *
ki
+ u* =
kd ud + c dud
(21)
Section 4.6
where
m
=
ie
--- m
----------------------- i --- 2
sin -- x
--- --*
L
(22)
The remaining steps utilize the results generated for the SDOF undamped structure
damped TMD system considered in Section 4.3. We use mie and kie as the mass
and stiffness parameters for the primary system.
To illustrate the procedure, consider the damper to be located at midspan,
and the first mode is to be controlled. Taking i = 1 and x* = L 2, the
corresponding parameters are
i x
= 1
sin -- ---- --*
L-
(23)
L
m
---------2
E
I
L
ie
= m
=
(24)
(25)
We specify the equivalent damping ratio, e , and determine the required mass ratio
from Figure 4.32. For example, taking e = 0.06 requires m = 0.03. The other
parameters corresponding to
= 0.03 follow from Figures 4.29, 4.30, and 4.31.
m
d = 0.965
f opt = ----1
d opt
= 0.105
d
---- - = 5
*
u
(26)
(27)
(28)
Using these parameters, the corresponding expression for the damper properties are
md =
0.03m
1
d =
kd =
0.9651
2
cd = d
m
2 d d m d
(29)
(31)
(30)
(32)
270 Chapter 4
Once m
and 1
(33)
E = 2 1011 N m2
The beam parameters
are
m 1
10000 kg
=
1 = 9.87 rad
s
(34)
(35)
The total mass of the girder is 20,000 kg. Adding 300 kg, which is just 1.5% of the
total mass, produces an effective damping ratio of 0.06 for the first mode response.
The mode shape for the second mode has a null point at x = L/2, and therefore
locating a tuned mass at this point would have no effect on the second modal
response. The optimal locations are x*
=
L 4 and x*
=
3L 4. Taking x* = L 4
and i = 2, we obtain
i x
sin -- ---- --* - = 1
LL
(36)
2e
k2e
= m
=
---------2
= k2 =
8EIL
(37)
- L
1 6 E
I
= -- --- ----- - -4 --
--
(38)
L
The procedure from here on is the same as before. We specify e and determine
the required mass ratio and then the frequency and damping parameters. It is of
interest to compare the damper properties corresponding to the same equivalent
damping
ratio. Taking e = 0.06, the damper properties for the example steel beam are
Section 4.7
md = 300 kg
kd =
(40)
(41)
435440 N
cd = m
(42)
2400 N s m
The required damper stiffness is an order of magnitude greater than the corresponding value for the first mode response.
CASE STUDIES: MDOF SYSTEMS
This section presents shear deformation profiles for the standard set of building
examples defined in Table 2-4. A single TMD is placed at the top floor and tuned
to either the first or second mode. The structures are subjected to harmonic
ground acceleration with a frequency equal to the fundamental frequency of the
buildings, as well as scaled versions of El Centro and Taft ground accelerations.
As expected, significant reduction in the response is observed for the harmonic
excitations (see Figures 4.46 through 4.49). The damper is generally less
effective for seismic excita- tion versus harmonic excitation (see Figures 4.50
through 4.61). Results for the low period structures show more influence of the
damper, which is to be expected since the response is primarily due to the first
mode. This data indicates that a TMD is not the optimal solution for controlling
the motion due to seismic excitation.
1
Building 1
Quadratic based Initial
H = 25 m
pm = 20,000 kg/m
s = 0.15
Sv = 1.2 m/s
1 = 2%
Harmonic TMD Mode 1
0.9
0.8
x
H
0.7
Normalized height
4.7
0.6
0.5
0.4
0.3
0.2
m = 0%
m = 1%
0.1
0
0.001
0.002
m = 2%
m = 5%
"{ *
0.003
0.004
0.005
0.006
0.007
0.008
Maximum shear deformation "{ (m/m)
0.009
0.01
272 Chapter 4
1
Building 2
Quadratic based Initial
H = 50 m
pm = 20,000 kg/m
s = 0.25
Sv = 1.2 m/s
1 = 2%
Harmonic
TMD Mode 1
0.9
0.8
Normalized height
x
H
0.7
0.6
0.5
0.4
0.3
0.2
m = 0%
m = 1%
m = 2%
m = 5%
0.1
0
"{ *
0
0.001
0.002
0.003
0.004
0.005
0.006
0.007
Maximum shear deformation "{ (m/m)
0.008
0.009
0.01
Normalized height
H
0.9
0.8
0.7
H = 100 m
pm = 20,000 kg/m
s = 0.40
Sv = 1.2 m/s
1 = 2%
Harmonic TMD Mode 1
0.6
0.5
0.4
0.3
0.2
m = 0%
m = 1%
0.1
0
0.001
"{*
0.002 0.003 0.004 0.005 0.006 0.007
Maximum shear deformation "{ (m/m)
m = 2%
m = 5%
0.008
0.009
0.01
Section 4.7
1
Building 4
Quadratic based Initial
H = 200 m
pm = 20,000 kg/m
s = 0.63
Sv = 1.2 m/s
1 = 2%
Harmonic TMD Mode 1
0.9
Normalized height
H
0.8
0.7
0.6
0.5
0.4
0.3
0.2
m = 0%
m = 1%
m = 2%
m = 5%
0.1
0
0.001
0.002
"{ *
0.003 0.004 0.005 0.006 0.007
Maximum shear deformation "{ (m/m)
0.008
0.009
0.01
Normalized height
H
0.9
0.8
0.7
H = 25 m
pm = 20,000 kg/m
s = 0.15
Sv = 1.2 m/s
1 = 2%
El Centro TMD Mode 1
0.6
0.5
0.4
0.3
0.2
m = 0%
m = 1%
0.1
0
0.001
"{ *
0.002 0.003 0.004 0.005 0.006 0.007
Maximum shear deformation "{ (m/m)
m = 2%
m = 5%
0.008
0.009
0.01
274 Chapter 4
1
Building 1
Quadratic based Initial
H = 25 m
pm = 20,000 kg/m
s = 0.15
Sv = 1.2 m/s
1 = 2%
Taft
TMD Mode 1
0.9
Normalized height
H
0.8
0.7
0.6
0.5
0.4
0.3
0.2
m = 0%
m = 1%
m = 2%
m = 5%
0.1
0
0.001
0.002
"{ *
0.003 0.004 0.005 0.006 0.007
Maximum shear deformation "{ (m/m)
0.008
0.009
0.01
Normalized height
H
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0.2
m = 0%
m = 1%
m = 2%
m = 5%
0.1
0
0.001
0.002
"{ *
0.003 0.004 0.005 0.006 0.007
Maximum shear deformation "{ (m/m)
0.008
0.009
0.01
Section 4.7
1
0.9
x Hheight
Normalized
0.8
0.7
0.6
0.5
0.4
0.3
m
m
m
m
0.2
0.1
0
0.001
0.002
0.003
0.004
"{ * 0.006
0.005
0.007
0.008
=
=
=
=
0%
1%
2%
5%
0.009
0.01
x Hheight
Normalized
1
0.9 Building 3
Quadratic based Initial
H = 100 m
0.8 pm = 20,000 kg/m
s = 0.40
0.7 Sv = 1.2 m/s
1 = 2%
El Centro
0.6 TMD Mode 1
0.5
0.4
0.3
m
m
m
m
0.2
0.1
0
"{*
0
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
=
=
=
=
0%
1%
2%
5%
0.009
0.01
276 Chapter 4
1
Building 3
Quadratic based Initial
H = 100 m
pm = 20,000 kg/m
s = 0.40
Sv = 1.2 m/s
1 = 2%
Taft
TMD Mode 1
0.9
x Hheight
Normalized
0.8
0.7
0.6
0.5
0.4
0.3
m
m
m
m
0.2
0.1
0
0.001
0.002
0.003
"{*
0.004 0.005
0.006
0.007
=
=
=
=
0.008
0%
1%
2%
5%
0.009
0.01
x Hheight
Normalized
0.2
0.1
0
"{*
0
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
=
=
=
=
0%
1%
2%
5%
0.009
0.01
Section 4.7
1
Building 3
Quadratic based Initial
H = 100 m
pm = 20,000 kg/m
s = 0.40
Sv = 1.2 m/s
1 = 2%
Taft
TMD Mode 2
0.9
x Hheight
Normalized
0.8
0.7
0.6
0.5
0.4
0.3
m
m
m
m
0.2
0.1
0
0.001
0.002
0.003
"{*
0.004 0.005
0.006
0.007
0.008
=
=
=
=
0%
1%
2%
5%
0.009
0.01
x Hheight
Normalized
0.8
0.7
0.6
0.5
0.4
0.3
m
m
m
m
0.2
0.1
0
"{*
0
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
=
=
=
=
0%
1%
2%
5%
0.009
0.01
278 Chapter 4
Building 4
Quadratic based Initial
H = 200 m
pm = 20,000 kg/m
s = 0.63
Sv = 1.2 m/s
1 = 2%
Taft
TMD Mode 1
0.9
x Hheight
Normalized
0.8
0.7
0.6
0.5
0.4
0.3
m
m
m
m
0.2
0.1
0
0.001
0.002
0.003
"{*
0.004 0.005
0.006
0.007
0.008
=
=
=
=
0%
1%
2%
5%
0.009
0.01
Building 4
Quadratic based Initial
H = 200 m
pm = 20,000 kg/m
s = 0.63
Sv = 1.2 m/s
1 = 2%
El Centro
TMD Mode 2
xH
Normalized
height
0.8
0.7
0.6
0.5
0.4
0.3
m
m
m
m
0.2
0.1
0
"{*
0
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
=
=
=
=
0%
1%
2%
5%
0.009
0.01
Problems 279
1
Building 4
Quadratic based Initial
H = 200 m
pm = 20,000 kg/m
s = 0.63
Sv = 1.2 m/s
1 = 2%
Taft
TMD Mode 2
0.9
x Hheight
Normalized
0.8
0.7
0.6
0.5
0.4
0.3
m
m
m
m
0.2
0.1
0
"{*
0
0.001
0.002
0.003
0.004
0.005
0.006
0.007
0.008
=
=
=
=
0%
1%
2%
5%
0.009
0.01
PROBLEMS
Problem 4.1
Verify Eqs. (4.13) through (4.17). Hint: Express p, u, and ud in complex form
p = p e it
u = ueit
ud = u e it
d
and solve Eqs. (4.6) and (4.7) for u and u d. Then take
u =
ud =
u e i1
u d e i1 + 2
= d =
280 Chapter 4
Problem 4.2
Refer to Eqs. (4.14) and (4.20). Express e as a function of m, , and u u
d.
Take = 0.05, and plot e versus
m
of the displacement ratio, u u d .
Problem 4.3
Figure 4.7 illustrates an active tuned mass damper configuration. The damper
can be modeled with the 2DOF system shown in Figure P4.3. The various terms are
as follows: us is the total displacement of the support attached to the floor beam;
Fa is the self-equilibrating force provided by the actuator; mdkdcd are parameters for the damper mass; ka and ma are parameters for the auxillary mass.
kd
ka
R
Fa
cd
md
ma
us + ud
us
us + ud + ua
FIGURE P4.3
(a) Derive the governing equation for md and ma . Also determine an
expression for the resultant force, R, that the system applies to the
floor beam.
(b) Consider ma to be several orders of magnitude smaller than md (e.g.,
ma = 0.01md ). Also take the actuator force to be a linear function of
the relative velocity of the damper mass.
Fa = c a u d
Specialize the equations for this case. How would you interpret the contribution
of the actuator force to the governing equation for the damper mass?
Problem 4.4
Design a pendulum damper system having a natural period of 6 seconds and
requiring less than 4 meters of vertical space.
Problems 281
Problem 4.5
L1
m1
m2
FIGURE P4.5
The pendulum shown in Figure P4.5 is connected to a second mass, which is
free to move horizontally. The connection between mass 1 and mass 2 carries only
shear. Derive an equation for the period of the compound pendulum and the length
of an equivalent simple pendulum. Assume the links are rigid.
Problem 4.6
Refer to Figure 4.12. Establish the equations of motion for the mass, md ,
con- sidering to be small. Verify that the equivalent stiffness is equal to Wd R.
Problem 4.7
Refer to Figure 4.15 and Eq. (4.76). Derive the corresponding expression
H
for P1 starting with Eq. (4.62) and using the same reasoning strategy.
Q
Considering
the mass ratio, m, to be less than 0.03, estimate the difference in
the optimal
values for the various parameters.
Problem 4.8
ranging from 0 to 0.2,
Generate plots of H1 versus for
m
d
f = 0.9876. Compare the results with the plots shown in Figure 4.23.
= 0.01, and
Problem 4.9
Consider a system composed of an undamped primary mass and a tuned
mass damper. The solution for periodic force excitation is given by [see Eqs.
(4.52) to (4.71)]
282 Chapter 4
u = ue it
(1)
ud = u d e it
(2)
u =
---H e i 1
(3)
k
p
---H
ud = e i3
k
(4)
3
2
2 2
f + 2f
d
-------------------- ------------H1 =
-----------------------D2
H3 =
2
--- ------D2
D2 =
1
f
(5)
(6)
2
2 2 2
m f
+ 2 d f 1 1 + m
(7)
The formulation for the optimal damper properties carried out in Section 4.3
was based on minimizing the peak value of H1 (actually H2 but H1 behaves in a
sim- ilar way) (i.e., on controlling the displacement of the primary mass).
Suppose the
design objective is to control the acceleration of the primary mass. Noting Eqs. (1)
and (3), the acceleration is given by
u = a = ae it
p
a = ---- --- 2 H e i +
--k
(8)
(9)
a = - ---- H e i1 +
m
(10)
where
2
H = H 1
(11)
shown in Figure 4.15, describe how you would establish the optimal values for
the various parameters, and also how you would design a tuned mass system
1
when H is
specified.
Problems 283
Problem 4.10
Design a TMD for a damped SDOF system having = 0.02. The design
motion constraints are
(a)
H
5 opt
10
H7
---- ----------- 5
H-5
opt
(b)
H5
opt
H7
--------- ------ 5
H-5
opt
(c) Repeat part (b), considering to be equal to 0.05.
Problem 4.11
This problem concerns the design of a tuned-mass damper for a damped
single degree of freedom system. The performance criteria are
eq =
0.1
u d u = 5
0.05
284 Chapter 4
(a) Determine the damper properties such that the equivalent damping
ratio for the fundamental mode is 0.16. Use the values of m, k, c
from Example 3.7. Assume stiffness proportional damping for c.
(b) Consider the tuned mass damper to be a pendulum attached to m5
(Figure P4.12). Determine md and L for the damper properties established in part (a).
u5
m5
L
k5, c 5
md
FIGURE
P4.12
(c) Repeat part (a) for the case where the mass damper is tuned for the
sec- ond mode rather than for the first mode, and the desired
equivalent modal damping ratio is 0.3. Use the same values of m, k, c
and assume stiffness proportional damping.
Problem 4.13
Consider a cantilever shear beam with the following properties:
= 50 m
= 20000 kg/m
DT
0 .6 x
--- - --- -- kN
8 10 1
H
5
Problems 285
L = 30 m
m = 1500 kg m
I = 1 10
m4
15 m
15 m
Constant EI
u1
FIGURE 4.15
Consider the simply supported beam shown in Figure P4.15. The beam has
a uniform weight of 15 kN/m and a concentrated weight at midspan of 100 kN.
The flexural rigidity is constant and equal to 200,000 kN-m2.
(a) Assume the first mode can be approximated by:
=
u
sin - --- x
u1 L