## Advanced Dynamics and Vibrations: Review of single and multi-degree of freedom (mdof) systems

### Introduction and Review

We are interested in small motions about some equilibrium positions. For SDOF systems, is usually measured from this static equilibrium position.

Let for static equilibrium. Therefore:

Let . Therefore:

All SDOF systems can be reduced to this form (assuming viscous damping):

Where means equivalent. To determine the nature of the motion, first consider this linearized equation of free vibration (Note: May have to linearize).

Assume (homogeneous), where is constant and is to be determined.

If are real and negative, then as . If either of or are real and positive, then the solution increases with time.

If are complex, then they are complex conjugates and the nature of the motion depends on the real part. If the real part is negative, as . If positive, as . If the real part is zero, the motion is oscillatory and this is the borderline case between stability and instability.

Consider all cases:

- . The roots have a negative real part and is stable regardless of . If is larger than , there will be an oscillatory part, but as . This is called asymptotically stable.
- . The roots are imaginary, and the system is borderline stable.
- . The roots have a positive real part, or are complex conjugates. Therefore, the system is unstable.
- One root has a positive real part. Therefore the system is unstable. The damping cannot overcome the negative spring.
- . The roots have a positive real part, and therefore the system is unstable. Negative damping feeds energy into the system.

Consider the solution for free damped vibration:

Drop , but remember that ,, are equivalent. Thus:

Again, assume that . Thus:

The critical value is then when:

Alternatively:

Define . Then:

Thus:

Where is the damping ratio (system parameter).

Assuming :

And using :

Using initial conditions , , we get the following:

We can write this as:

Where , , and .

In some problems, we have to find the equilibrium (static) configurations and then determine the motion about these positions.

#### Example

Find the motion about the equilibrium configuration.

slides on a wire with viscous friction .

##### Solution

We must find static equilibrium configurations and determine equation of motion for small perturbations about these static configurations.

Set , and let . Therefore:

Therefore. the static configurations are given by:

Assuming , there are 3 positions.

To find the motion, linearize (1) about . Set , where is a small angle.

Put the above into the original equation of motion to get:

Note that .

Using the fact that the term is small compared to 1, it becomes:

The above equation looks like the standard form except for a constant term that is zero for all conditions.

Now consider all three positions:

1. .

Therefore, this is stable if as , but unstable if .

2..

Therefore, this is always stable.

3..

Therefore, this is always unstable, regardless of damping.

###### Example

1., , . Therefore:

2., . Therefore:

### Equivalent spring constants

One of the components we need in these equations of motion is the spring constant . We can often find this for a system using static techniques. This is easily generalized for MDOF systems.

For SDOF systems, we can imagine the static response of the system using one of the approaches:

1.Apply a unit force (moment) to the mass (inertia) in the positive direction of motion then calculate the displacement that occurs (flexibility approach).

2.Apply a unit displacement (rotation) to the mass (inertia) in the positive direction of motion then calculate the force (moment) required to maintain it (stiffness approach).

It depends on the situation to determine the best one to use. Generally, for a series situation, use flexibility, while for parallel, use stiffness.

#### Springs in Series

Apply unit load, and calculate .

#### Springs in Parallel

Apply unit deflection, and calculate load.

#### Axially Loaded Bar

Therefore:

Therefore:

#### Inclined Axial Spring

Apply a in the direction of motion.

Therefore, force in spring when extended is , and thus:

#### Example 1

Stiffness

Apply unit deflection, and calculate total load.

Therefore:

If , then:

#### Example 2

Flexibility

Apply unit load, and calculate total deflection.

, therefore:

Therefore :

Thus:

If , then:

#### Equivalent Spring Constant Equations

1. axial springs in parallel

2. axial springs in series

3.Springs in parallel and series

4.Inclined axial spring

5.Rotating bar with spring support

6.Rigid bar supported on two springs

7.Rigid bar supported on three springs

8.Axially loaded bar

Where is the cross-sectional area, and is the elastic modulus.

9.Axially loaded tapered bar

10.Axially helical spring

Where is the active number of turns, and is the elastic shear modulus.

11.Torsion of a uniform shaft

Where is the torsional constant of cross section ().

12.Torsion of a tapered circular shaft

13.Spiral torsional spring

Where is Young’s modulus. is the moment of inertia of cross-sectional area, and is the total length of the spiral.

14.Cantilever bean, end load

15.Simply supported bean, load at midspan

16.Simply supported bean, load anywhere between supports

17.Fixed-fixed beam, load at midspan

18.Fixed-fixed beam, off-center load

19.Propped cantilever, load at midspan

20.Propped cantilever, load at free end

#### Table 2-1

###### Fixed-fixed beam*

###### Fixed-pinned beam with overhang*

###### Fixed-pinned beam with overhang (P at x = l + a)*

###### Pinned-pinned beam with overhang*

###### Pinned-pinned beam with overhang (P at x = l + a)*

###### Fixed-fixed beam with lateral displacement

* – Axial extensions due to axial end constraints considered negligible

#### Torsional Oscillations

Let . From the strength of the material:

Where is the modulus of rigidity, is the angle of twist/unit length and is the polar second moment of inertia. Therefore:

For soild circular cross sections:

### Center of Percussion and Vibration of a Compound Pendulum

The concept of the center of percussion is usually discussed using the analysis of the changes in angular momentum and linear momentum associated with a body impacting a rotating object and essentially determining the point at which there is no change in acceleration at the pivot point due to this impact. (The usual example is a ball hitting a baseball or cricket bat.) For the compound pendulum shown the natural frequency is where is the moment of inertia of the body about the axis of rotation at . Instead of having a distributed mass, let the mass of the pendulum be concentrated at a single point located at a distance from the pivot point . The natural frequency of this simple pendulum is . If the natural frequencies of these two pendulums are the same, then the point at which is concentrated is called the center of oscillation and is identical to the center of percussion. If instead, we consider rotation about the center of oscillation (percussion) then the original pivot at becomes the center of percussion (oscillation).

Early automobile manufacturers did not realize this and as a result, cars had much more uncomfortable rides because impacts at one set of wheels (front or rear) were felt at the other. One of the first automobile designs to incorporate this idea was the very futuristic Chysler Airflow (circa 1930) which was years ahead in its design in a number of aspects. (So far ahead that it was a commercial failure.) For the illustration of the automobile shown the relationship between the position of the center of gravity and the positions of the front and rear wheels is given by where is the radius of gyration of the vehicle about its center of gravity and is the wheelbase ( and are the horizontal distances from the front and rear wheels to the center of gravity). From a vibration perspective if a car is designed with a mass distribution to meet these criteria (almost all recent automobiles are) then we can examine what happens at one set of wheels individually for certain situations (e.g. hitting a bump with the front wheels only).

Using the diagram below show that the condition for the center of percussion to be at the front wheels with the pivot at the rear wheels that and that this is also true for the pivot point moving to the front wheels with the impact at the rear ones.

Where is the moment of inertia about and is the radius of gyration about :

Therefore:

Therefore:

#### Equivalent Pendulum

Therefore choose:

#### Example

For application to the auto above (wheelcase). For rotation about :

Therefore:

Therefore:

For rotation about :

Therefore:

#### SDPF damped system

We know that for any SDOF damped system the equation for free vibration is:

and the system parameters are the natural frequency:

and the damping ratio:

We can determine the damping ratio from how rapidly the vibration decays. To do this we use the so-called logarithmic decrement. Consider the damped response:

which for looks like:

These peaks occur approximately when:

Therefore if the first occurs at then the next time it is at:

Therefore:

and this is the same ratio for , , etc.

Therefore this creates the logarithmic decrement equation:

If we take n cycles:

Therefore:

Therefore if we measure we can easily determine the damping ratio .

If

0.05 | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 | |

1.37 | 1.87 | 3.61 | 7.21 | 15.5 | 37.6 |

It does not take a large to quickly damp out vibrations.

### Steady State Harmonic Response

Therefore:

Therefore find solution as homogeneous solution and particular integral.

###### a) Homogeneous Solution:

###### b) Particular integral:

Try into

Therefore:

Note:

Therefore:

Solving the above equations for C & D gives:

when and :

However, if is after some time the homogeneous part will decay and the steady state solution becomes:

If we take the special case of no clamping, the homogeneous equation is:

and the particular integral is

Therefore:

Therefore:

Therefore the total solution is:

Note that as there is no damping present that the homogeneous solution does not as time increases.

and we combine them to get:

Where:

This result is extremely useful in practical situations especially for controlling noise and vibration.

While this tells us the dynamic motion of the system, there is still the static deflection of the springs to be remembered. These are also related situations which are obtained by reconsideration of the fundamental equation of motion.

From the view of real situations, these are other quantities that are of interest. One of these is the force transmitted to the supporting structure.

Where is:

Therefore:

and the magnitude of using the fact that is:

A second practical situation accuse when the structure is vibrating and we wish to isolate a machine from the vibration (base excitation).

So the equation of motion is:

The steady state solution comes from the assumption from :

and the solution is:

A third potential situation is essentially the original forced situation where the magnitude of the forcing function is due to rotating imbalance:

Where past of the mass , at an eccentricity causes the dynamic force.

Equations of motion:

This looks like the previous forced case where is replaced with so that:

and reduced to non-dimensional form as:

so that the steady state solution is:

where:

The frequency response for all these cases are summarized on the following page

#### Single Degree of Freedom Systems – Useful Relations

The frequency response curves are helpful in visualizing where the specific system you are analyzing lies. It must be remembered that there is a phase difference between excitation and response with the change from being nearly in phase to being nearly 180 degrees out of phase happens near resonance (). For the undamped case, this means the response is undefined at and the response becomes:

The results for the steady state forced single degree of freedom systems are usually present in the form of graphs showing various non-dimensional parameters. Because of the technical importance of these results, it is necessary to understand the differences between these representations. Two of the most useful are given by the expressions for displacement of the mass:

And the displacement of the main mass in the case of rotating imbalance:

a) Show that (B) can be derived from (A) for the special case in which the forcing function is due to rotating imbalance

b) For the undamped case determine the result of (A) and (B) both below and above resonance

a) Special cases for rotating imbalance:

If :

Note: includes the imbalanced mass

b) for is:

There are two important physical phenomena that are easily illustrated using the forced undamped SDOF system as we approach resonance. The entire solution:

Consider the case

,

. From the first of these and

Therefore at t = 0:

If , where is small compared to or , then:

Known , therefore:

Therefore:

This expression can be considered as representing a vibration with period that has a slowly varying amplitude equal to , ( has a much longer period than )

While the original solution to the forced undamped differential equation is not valid at , it is of value to consider what happens at 0 for the displacement. Using L’Hospital’s rule for the following limit:

The amplitude of vibration increases indefinitely with time and the infinite amplitude is in the steady state. This means that to operate the system it is important to move through resonance as quickly as possible so that the amplitude does not “build up”.

### Energy Considerations in Vibrations

For many dynamic/vibration systems energy considerations are useful especially in modelling. This becomes evident when the complexity of the system increases as the formulation of the models will usually involve velocities rather than accelerations and the kinematics for general motion are easier to describe. Specifically for vibrations we can appreciate the exchange of energy occurring between the kinetic and potential components with the damping providing the energy dissipation. This occurs in the consideration of the Lagrange formulation of dynamical systems

If we choose the static equilibrium position as the datum for the energy calculation we must be aware of the energy involved in the static equilibrium deflection. There will be instances in which the static condition does not have a deflection then the energy calculation must include consideration of the position of the mass of the system.

**One DOF Spring Oscillator**

For a single DOF undamped system we can calculate the natural frequency from energy considerations as we know the motion is Simple Harmonic Motion (SHM). In addition it is possible to use energy as an approach to finding approximate solutions for the natural frequency without a large computational effort. For these systems, termed conservative, the sum of the kinetic and potential energy for free vibration is a constant.

T (kinetic) + U (potential) = C (constant)

In addition the energy varies over time from being all potential to all kinetic so that

Under the assumed SHM

Consider the SDOF system from an energy perspective where a gain * *is measured from the static equilibrium configuration.

Potential energy is the work done by F in moving through (energy in spring) :

While the change in potential energy due to position of is . Therefore, the total potential energy is , assuming :

Therefore,

We can extend this idea to consider approximating the effect on the natural frequency when the mass of the spring is included. Assume that the spring of length has a uniform mass over its length. Then at any position on the spring the velocity is and the kinetic energy of an infinitesimal mass is and the kinetic energy of the spring is:

If we now make an assumption of how the spring “moves” we can calculate e.g. assume the spring displaces linearly with position y so that at

Now using :

If instead we assumed the “spring” was a simply supported beam of mass and length and that the dynamic shape is a half-sine then at any position the velocity is: and the mass/unit length is

Then using the same approach:

### Complex Representation for Vibrations

When dealing with vibrating equations it is often more efficient to use complex numbers as they allow the phase relationships to be easily handled. To do this use the complex number:

is called the complex conjugate

Alternatively we write as:

is called a unit phasor. For fractional expressions we can rationalize the denominator to separate real & imaginary parts:

Consider the steady-state forced SDOF system:

and we want to imaginary part of the solution

Try , where is complex, therefore:

Note: , is defined as the mechanical impedence.

Rationalizing the denominator gives:

as

Then the imaginary part of the solution is:

We can now re-interpret equation as:

Each of the force terms above can be interpreted in their respective phases.

Many real automotive suspension systems use the so-called Macpherson strut which is a combined spring and damper along with a coil spring as shown in the diagram. The system can be represented by the diagram where the forcing function comes from road surface. (This system is sometimes called a 1 & DOF system).

a) Using the mechanical impedance to show that the displacement and .

b) For the model shown, determine the transmissibility of the system.

c) For the frequency ratio , and for . Calculate for the transmissibility and compare the values to the case in which .

a)

Using

But (3) is of the form:

Therefore:

If :

, if :

If we set ,

Therefore,

Where

Where

b) The force transmitted is

c) For the case :

For the case :

The added spring reduces the transmissivity in both cases, however it is most effective at higher ratio.

If

## Extra Details

If :

Where , , and . Thus:

Where:

If , :

### Aerodynamics Excitation

A relatively common excitation which occurs in diverse situations from subsea pipelines to slender tall structures is due to vortex shedding. This shedding creates a differential pressure perpendicular to the fluid flow. For lower Reynolds numbers (Re 300) the vortex shedding frequency is estimated to the Strouhal number by

where is the shedding frequency, is the velocity of the external flow and is a characteristic dimension (e.g. diameter of the building or the pipeline). Assuming that the lift force per unit length is harmonic it can be written as

where lift coefficient, density of the fluid (mass), and

For the case of a bluff body, the total force to the flow is modelled as

where is the frontal area to the flow. Because of the complexity of the flow there can be components of alternating forces at 1/2 the frequency estimated by the Strouhal number as well as at twice the .

It should be noted that the shape of the body can also have a major influence on the effective force magnitude. Catastrophic examples include the Tacoma Narrows Bridge and the destruction of powerlines that change their X-sectional shape due to icing.

For most of these lower Reynold’s number flows

3. In order to produce a relatively uniform acoustic field to the audience in an open stadium, the sound system is often suspended above the center of the playing surface. (This allows the distance from the source to the audience to be approximately the same for all fans.) This is shown diagrammatically, with typical dimensions in Figure 1. Especially in northern climates, there is the potential for the system to be covered in ice and the usual porous surfaces are replaced by an effectively solid “bluff” body that can be excited as a result of wind induced vortices being shed from the body. Alternate shedding from the upper and lower surfaces cause a pressure differential and therefore a vertical excitation of the body. The frequency f of this excitation is related to the wind velocity v and a representative vertical dimension D (taken as the total height of the system in this case, Figure 2) through the Strouhal number S given by

(S = 0.22 for this outdoor situation). The concern is that the frequency of vortex shedding could match the natural frequency and result in large even catastrophic oscillations of the sound system (google Tacoma Narrows Bridge).Add block

(a) For a cable of length L,

therefore for the 4 cables, the for the vertical motion

wind speed:

(b) For low damping ratio:

Therefore:

(c) For the steady state SDOF damped system at

For resonance at the vortex shedding frequency:

now using:

at Edmonton altitude

Therefore:

Therefore total movement is just over 1m.

## VIBRATION ANALYSIS AND TESTING OF THE COMMONWEALTH GAMES STADIUM SOUND SYSTEM

D.J. Wilson and M.G. Faulkner

Department of Mechanical Engineering

University Of Alberta, Edmonton, Alberta, Canada

ABSTRACT — The possible effects of wind—induced oscillations of the Commonwealth Games Stadium sound system were analyzed before construction. Tests after installation confirmed the original estimates of the system’s damping and natural frequencies.

INTRODUCTION

The Commonwealth Games Stadium sound system is one of the most sophisticated centrally suspended systems ever installed in an open—air athletic stadium. The 6350 kg (14,000 lb) speaker system is suspended above the centre of the playing field by four 110 m (360 ft) bridge strand cables. Figure 1 shows schematically the system and supporting cables. The central location was chosen because of the acoustical advantages of a frontal sound field and the lower power necessary to provide adequate sound pressure levels for all the spectators.

With the structure located 30 m (100 ft) above field level there existed the possibility of wind induced vibrations. While wind tunnel model studies are possible to estimate the magnitude of these vibrations, the relatively simple configuration lent itself to using computations alone. Calculations of the resonant frequencies of the system were combined with the methods outlined in Commentary B of the National Building Code of Canada [1] for estimating wind excitation forces.

The final design chosen has the speaker cluster enclosed by a 2.5 m (8 ft) radius cylinder 3 m (10 ft) high with a solid dome on the top and an expanded metal mesh on the sides and bottom. In addition struts (shown in Figure 2) were added for stability and to increase certain of the resonant frequencies.

DESIGN CONDITIONS AND ANALYSIS

The idealized configuration of the speaker cluster shown in Figure 2 was used for design purposes. From a structural viewpoint this meant a annular cylinder was suspended as shown by four 25 mm (1 inch) diameter bridge strand cables. Because of the symmetry, this structure was assumed to vibrate in any of six modes of vibration (see Figure 3) independently.

From an aerodynamic viewpoint the wind velocities for a 30 year return period were chosen. This is recommended by the National Building Code (NBC) because of the possible hazard to occupants of the Stadium. Based on data in Supplement #1 of the NBC the maximum one hour average wind speed for Edmonton is 90 kph (56 mph) for a 30 year return period. The most severe condition occurs when the holes in the expanded metal mesh covering the structure are open to the wind. The site exposure was taken as urban and an attempt made to account for the influence Of the spectator stands on the wind patterns.

Natural Frequencies of Vibration

The calculated natural frequencies of vibration for each of the six modes are shown in Figure 3. Initial estimates put the level of uncertainty for the translational modes at 10% while the pitching and torsional modes had levels of 25%. The latter uncertainty was due to the difficulty of estimating the mass moment of inertia for the system.

After initial calculations, the compression struts connecting the two easterly cables and the two westerly ones were included to provide additional torsional stiffness. Addition of these struts approximately doubled the frequencies for the torsional, pitching and rocking modes.

Wind Loads

The calculation of the wind loads was done using the procedures recommended in the NBC Commentary B. Two contributions to the wind load, gust buffeting and vortex shedding were considered.

Using the 30 year return wind speed and the wind exposure factors in the NBC gust buffeting loads were computed. Table I shows the parameters related to peak gust buffeting loads. For all wind directions the fluctuation frequency due to gust buffeting is far removed from the resonant frequency resulting in small gust buffeting loads.

Preliminary estimates indicated that vortex shedding would cause the largest deflections and cable loads. For a fixed geometric shape the maximum deflection is related to the aerodynamic and vibration parameters according to the proportionality;

Max. resonant deflection is

where

lift coefficient

damping ratio

Strouhal number

As seen from the above equation the damping ratio is a most important factor. In the initial design this was taken as 0.01 for vertical oscillation. This damping was assumed to be from aerodynamic sources (0.003) and an equivalent viscous damping of 0.007 from the structure. The values of the lift coefficients and Strouhal number were estimated by a literature survey. The lift coefficient evaluated in this manner should represent an upper limit since limited water channel tests of a sound system model showed no evidence of vortex shedding.

Of the six modes of vibration only vertical motion and a pitching motion about the east-west axis (see axis in Figure 3) had natural frequencies low enough to produce critical airport windspeeds for resonant vibration of less than 160 kph (100mph). Table II shows critical windspeeds equivalent maximum static loads and deflections caused by vortex shedding.

FULL SCALE TESTING

The lack of information on various system parameters caused considerable uncertainty in the design analysis. The estimates of the cable support stiffness, and the damping ratio made the resonant frequency and maximum amplitude calculations suspect. To confirm the previous estimates of these parameters a field test was carried out to directly measure the damping ratio and to confirm the previous estimates for certain of the natural frequencies.

Test Procedure

The damping ratio and natural frequencies were evaluated by using a motion picture camera to record artificially excited vibrations. The film was then projected on a motion analyzer enabling the amplitude, period of oscillation and decay rate to be measured.

The final filming was done with a Cinema Products CP-16 camera with a crystal controlled framing speed and a zoom lens. It was estimated that in conjunction with the motion analyzer the precision and resolution of amplitudes would be 10 mm. The crystal speed control maintained the base accuracy with a variation of only a few parts per million which is negligible for the present application.

Excitation of the system was provided by normally pulling on wires attached to the base of the sound system. One was located under the mass center and pulling on it produced nearly vertical oscillation alone. The other wire was attached at the northern edge of the system. Pulling on this wire produced a combination of pitching about the east-west axis and vertical oscillation. Amplitudes of 300 mm were easily obtained by the manual excitation. The camera recorded approximately 1 minute of decaying vibration after the excitation had ceased.

Data evaluation

The amplitude and period of decaying vibration was measured using a rear projection analyzer. A moveable cross-hair cursor for the projected image allowed a frame by frame analysis of the motion. Time periods between successive maxima were used to evaluate the natural frequency.

The damping ratio was determined by a semi-log plot of the relative amplitudes. The log decrement was evaluated using a best fit straight line (See Figure 4). The fit to a straight line verifies the assumption of equivalent viscous damping. The log decrement is related to the damping ratio by

Using the average of two tests the damping ratio for vertical motion was found to be

with an estimated uncertainty of 10%.

For the pitching mode, a beat formed between the vertical and pitching modes made it impossible to measure the damping in the pitching mode. The natural frequency however could be evaluated. A comparison between the measured and predicted values is shown in Table III.

CONCLUDING REMARKS

The results of the initial design analysis and the results of field testing shown in Cables i, ii, and iii show good agreement. The damping ratio originally estimated at 0.01 compares to the sum of measured values of 0.0054 and calculated aerodynamic damping of 0.0031 (total of 0.0085). The measured and calculated natural frequencies for the vertical and pitching modes are within 10-15% of those measured. Of the six modes of vibration these two are the only ones which could occur at airport wind speeds of under 160 kph.

The system has been in place nearly two years and while some deflections have been observed these have not approached the maximum predicted leading one to believe that estimates of the increased loads and deflections were quite conservative.

ACKNOWLEDGEMENTS

The authors would like to thank Ragan, Bell, McManus Consultants Ltd. for suggesting our involvement in this project. In particular we would like to mention Stan Ragan and Dan O’Brien.

REFERENCES

[1] National Building Code of Canada 1977 NRCC No. 15555 including Supplement No. 1 NRCC No. 15556 and Supplement No. 4 NRCC No. 15558

Gust Buffeting at 90 KPH Design Hourly Airport Windspeed (30 Year Return)

Wind Direction | Max 5 sec Gust on System (KPH)MPH | Total System Damping | Max Gust Load kN | Max Deflection mm | System Natural Frequency Hz | Fluctuation Frequency Hz |

North or South | (132)82 | 0.0073 | 11.7 | 4.0 | 1.04 | 0.24 |

East of West | (85)53 | 0.0058 | 4.8 | 2.5 | 2.68 | 0.68 |

Table I

Vibration Induced by Vortex Shedding

Direction of Motion | Natural Frequency Hz | Critical Airport Windspeed, North or South Direction (KPH)MPH | Critical Airport Windspeed, East or West Direction (KPH)MPH | Critical Windspeed at System (KPH)MPH | Maximum Equivalent Static Load kN | Half Amplitude mm |

Vertical Z | 0.48 | (63)39 | (135)84 | (53)33 | 38 | 795 |

Pitching about X-(West) Axis | 0.53 | (69)43 | (150)93 | (58)36 | unknown | unknown |

Table II

Comparison of Measured and Computed Resonant Frequencies

Direction of Motion | Computed , Hz | Measured , Hz | |

Vertical Translation along Z axis | 0.42 | 0.48 | 1.14 |

Pitchings: Rotation about X-West axis | 0.49 | 0.53 | 1.08 |

Table III

### Transient Vibrations

For dynamical systems excited by non-periodic forces, displacements, accelerations, etc., the response over the initial few periods of the system is termed transient response. This is important for many diverse situations including packaging of equipment, protecting people from automobile impacts, and loading of structures due to blasts and earthquakes. The modeling of these situations is often difficult and may require specific detail of the excitation as the response can vary dramatically depending on the detail of the input.

For a SDOF system the solution requires the general solution of and becomes difficult for complicated

###### Example

Consider first the undamped case for a step input:

###### Solution

Assume

and the total solution:

If , then:

While this case is straightforward, the solution is more difficult with damping or for different forcing functions.

###### Example

The response can be determined in a general form by considering as being made up from a series of short impulses. A general short impulse at some time is shown. Consider what happens(response) if __only__ this impulse is applied to the system:

If is vanishingly small the impulse can be considered as changing the velocity of using impulses/momentum relations so that from initial (just before) to final condition:

Therefore:

This impulses will result in a free vibration with initial conditions:

If this is the only impulse then the response is:

with

Therefore:

Where is the time since the impulse.

Therefore the response at due to the impulse at is :

The response due to all the impulses from to is the superposition of all the impulses’ response.

Therefore:

Assuming i.e. the mass was stationary prior to being applied.

For the step input

If the initial conditions are , then:

#### Example

For a package being dropped or for an airplane landing, a SDOF model can be utilized.

For the undamped case, as it strikes, and , and the forcing function is (constant). NOTE: is not measured from the static equilibrium configuration. Using the general result:

Therefore, the maximum acceleration is:

With:

In forensic studies of automobile collisions, the concern is the acceleration of a passenger’s head and the potential for internal damage. Assume the body/skull is modelled as a SDOF system.

Biomechanical studies gives the spinal stiffness . Assume the person weighs 160 lbs and falls 6 inches because of being unrestrained. Thus:

Therefore:

If we add a seat cushion to the seat (cushion stiffness is ). What is the change in maximum acceleration?

The effective stiffness is now:

Thus:

Therefore:

The cushion reduces the acceleration drastically.

For many situations, the support of the system is subjected to the motion specified by its displacement, velocity, or acceleration.

Therefore, for base acceleration:

Using the relative displacement :

The convolution integral for this case is:

The same approach can be used for a base excitation with velocity.

##### Examples of Dynamic Responses of Undamped SDF Systems to Different Pulse Forces

Note that the static solution is shown in the dashed lines.

###### Rectangular Pulse Force

###### Half Cycle Sine Pulse Force

###### Triangular Pulse Force

## Shock Response Spectrum

Previously we developed a technique to find the response of a damped spring mass system to an arbitrary excitation. When the duration of the pulse is small compared to the natural period (), the excitation is called a shock. It has been found that the concept of the shock response spectrum is useful for these situations.

The Shock Response Spectrum (SRS) is a plot of the maximum peak response of a SDOF oscillator as a function of the natural period of the oscillator. The maximum of the peaks (called maxima) represents only a single point on the time response curve.

To illustrate the SRS concept consider a SDOF undamped system with no initial motion (). The SRS is

If F(t) is a rectangular pulse of length

then during the pulse

where and after

where . Often the results are determined for the two sections (during & after) individually then combined as the primary response (during) and the residual (after pulse) response.

During the pulse

Where for a maximum:

therefore:

therefore up to , the maximum is the value at the time considered as it is still increasing therefore:

where after the maximum is:

These results are called the primary values as they correspond to the result during the pulse.

after

where

and we find the maxima by differentiating

at

therefore:

This value for must be inserted into the original expression to find the maximum displacement. This will require and in terms of and .

therefore consider the triangle

[ as tan is negative], sin is negativetherefore:

now

but therefore:

Residual Response

The total SRS can now be determined

If damping is added then the magnitude can be reduced considerably. However, the calculation of the SRS must be done numerically

For cases which the pulses variation is less than one-half the period , the overall maximum occurs after the pulse in the “free” vibration phase. As the pulse duration compared to the period becomes smaller we can approximate the response where it is the magnitude of the impulse(not the details of the shape) that matters. In the development of the theory, the response for a damped system was:

Where is the time since the impulse and for for the undamped situation

Where is the impulse, I , and the maximum displacement is:

Consider three different impulses:

1) rectangular 2) half sine 3) triangular

1)

Therefore, compared to the static response due to

2)

Therefore,

3)

If the effect of damping is included then:

and for we assume the maximum occurs at (or near) where is 1.

Therefore (time at ) is of a so-called damped period:

Therefore:

and for :

and we can multiply the undamped values by :

0.05 | 0.10 | 0.15 | 0.20 | 0.25 | 0.3 | |

0.921 | 0.843 | 0.764 | 0.686 | 0.607 | 0.529 |

### Vibration and Applications of MDOF Systems

For more than a SDOF system the analysis requires a more generalized approach. However, overall it doesn’t matter if there are 2 or 22 DOFs. Consider first a special 2 DOF free vibration system:

The simplest approach is to look for Simultaneous Simple Harmonic Motion (SSHM). That is:

Which is a solution (not the general one) if:

This can be true if and only if (iff), . Therefore:

There are two soloutions for .

Therefore:

As a result, there are 2 frequencies at which our assumption is true. NOTE: It turns out that the general solution can be determined from these , , and the ratio of the amplitudes between the two masses during each of the two SSHMs. Therefore:

Or:

Therefore, for either or :

If we put into either we get:

And if we put into either we get:

Thus, we define the mode shapes corresponding to each ratio:

Therefore, we have a mode shape corresponding to each of the natural frequencies. and are “normalized” vectors where the first component is set arbitrarily to unity. The general solution to the free vibration problem is then:

There are 4 constants , , , for the 4 initial conditions for , , , . If the initial conditions are selected then the solution will include both modes in general but only one mode for certain cases. For example, when , , , and , mode 1!

Consider the following initial conditions:

At :

Therefore:

Therefore:

From (1) and (2) (+):

From (1) and (2) (-):

Therefore:

Using identities:

These may be interpreted as a higher frequency oscillation at the average of the two natural frequencies with a variable amplitude, given by a lower frequency given by the difference in natural frequencies.

If and are close to each relative to their magnitude then , the motion becomes:

### Stiffness and Flexibility Influence Coefficients

We can use the Newtonian approach to find the stiffness matrix and equations of motion. However, we can also use influence coefficients for both the mass and stiffness matrices.

Consider the definition of stiffness and flexibility coefficients:

## S.I.C

is the force (moment) required at coordinate to maintain a unit linear (angular) displacement at coordinate with all other coordinates held fixed at zero

## F.I.C

is the linear (angular) displacement at coordinate due to a unit force (moment) applied at coordinate with all other coordinates free to move.

NOTE: All forces (moments) and displacements (rotations) are taken to be applied in the positive direction. The sign of the resulting forces (moments) or displacements (rotations) are determined to be positive or negative.

Also note that:

## Example

###### Stiffness

Therefore:

Therefore:

###### Flexibilty

Therefore:

Thus:

Therefore:

As an example of MDOF, consider a general planar situation which can be specialized to different industrial situations. Many of these problems can be initially modelled as SDOF ones, and then MDOF analysis will allow interpretations of the simplifying assumptions made.

NOTE: While the body shown is a rectangular shape, it can be any rigid body with a mass and moment of inertia about an axis to the plane through the center of gravity. The horizontal stiffness components, designated by are most often just the lateral stiffness of the vertically mounted spring. For real springs, the lateral spring may be hard to find in the literature. For coil springs, it is often in the range of times the longitudinal value.

In matrix form:

Part 2

Now find the stiffness matrix using influence coefficients.

First consider a unit displacement in the direction only:

For the direction unit displacement:

Assuming all summations for direction, direction and moment equal to :

Therefore:

This yields the same stiffness matrix.

As long as we use the motions about the center of gravity then the inertia(mass) matrix is just the mass on the mass moment of inertia about the center of gravity. If we use the coordinates to decide the motion that is not related to the absolute motion of the center of gravity then the stiffness matrix may not be symmetric and the mass matrix will not be diagonal.

#### Inertia Influence Coefficients

The elements of the mass matrix, can be determined through the use of inertia influence coefficients. They are defined as the set of impulses applied at the points (representing the coordinate directions) respectively to produce a unit velocity at point and zero velocity at every other point.

NOTE: If denotes an angular coordinate then represents an angular velocity

Thus for a multi-degree of freedom system, the total impulses at points is:

To find the

- Assume that a set of impulses are applied at all points as to produce a unit velocity at point only(with zero velocity in all the after dierections(coordinates)). By definition the set of impulses denote the influence inertia coefficients
- Repeat the procedure for each point .

###### Example

First apply unit velocity to only:

Apply unit velocity(angular) to only

Therefore:

Now find the stiffness matrix

Apply unit deflection to only:

Therefore

Apply unit deflection to only:

Therefore:

### Applications of MDOF systems

The concept of modes and their corresponding natural frequencies is useful in understanding the behavior of lightly damped systems. Before adding the complication of damping to our models, we will consider several applications to industrial situations.

###### Example

The model of the automobile shown is used to predict the low frequency modes of vibrations as these are responsible for the low frequency noise in the passenger compartment. Knowing these models can assist the designer in controlling them. In addition, it helps the designer in predicting the response of the frame to various inputs from road irregularities. As an initial estimate, the three modes of vibration of the vehicle may be considered as uncoupled. For the specific case below, the assumption of uncoupled modes is to be compared to the coupled model.

a) For the longitudinal model of the vehicle, specialize the equations of motion for the case of

b) If any of the modes are uncoupled determine an expression of the natural frequency. For the coupled modes derive the equations of motion for free vibration starting from Newton’s laws.

c) Calculate the natural frequencies and mode shapes for the two coupled modes of vibration. The weight of the vehicle is and the radius of gyration on the vehicle about is .

###### Solution

a)

b)

For horizontal motion:

Therefore:

& are the same as .

c)

Therefore:

Assume :

[DIAGRAM]

If we “uncoupled” the modes

Assume vertical node is

Rotational Mode

If the horizontal stiffness, , were the same as the vertical then the third frequency, , would be the same as the vertical.

The assumption of uncoupling gives essentially the same answer.

The illustration shows the detail of the proposed isolation system for the reciprocating compressor and motor considered previously and modelled as a single dof system. The actual system does not have the center of gravity symmetrically located on the supporting structure and isolation springs. In situations like these the industrial installation guides recommend that the isolation springs be selected so that the vertical static deflection of all springs is the same. The reason stated is that this will tend to reduce the frequency of the “rocking motions”. This exercise is to investigate this idea as well as the differences found for the isolation of the system compared with the single degree approximation used earlier. More details of the installation are given below. The stiffnesses shown are representative of both springs at the ends of the supports with the lateral stiffnesses assumed to be equal to the vertical values.

a) If the vertical static deflection of the spring is the same what is the ratio of and ?

For static equilibrium

therefore:

(b) Under this assumption what happens to the equations of motion?

as the vertical mode is uncoupled from the other two

becomes

this uncoupling solves the natural frequencies

(c) Noting that the original had an operating speed of 1750 RPM and the original transmissibility was 1/15, calculate for this installation assuming

therefore:

therefore:

therefore:

(d) Assuming SSHM

therefore:

Divide by to get

therefore:

(e)

The rotational mode can be calculated starting from a FBD

therefore:

(f) Do the calculation using a program

(g) From d)

from e)

from (f)

The results show that because is small compared to and that there is little difference. The “rotational” () is very much higher than the other two and closer to the excitation frequency of 29.2 Hz (1750 RPM).

To lower it, the spacing of the supports would have to be reduced. For example if it were brought to 1/2 of original (5.1′ instead of a total of 10.4′) then the rotational natural frequency would be

It is obvious from this calculation that the rule of thumb should be that the spacing should be no more than twice the radius of gyration of the total mass!

A commercial front loading washing machine is being used by a spa located adjacent to a medical clinic in a strip mall. During the spin cycle of the washing machine the fluctuating force due to the imbalance of the load causes the floor to vibrate and induces unacceptable levels of noise in the clinic. To reduce this disturbance, an engineering consultant recommends using spring isolation between the washing machine and the floor to reduce the dynamic force transmitted to 20% of the exciting force. The specifications of the washing machine are summarized below. Assume only vertical vibration and neglect damping.

Net weight | 1049 lbs |

Maximum load | 255lbs |

Maximum Dynamic Force | 212 lbs |

Spin Speed (Frequency of Dynamic Force) | 950 RPM (15.8 Hz) |

Cylinder Diameter | 27.6 inch |

a) Calculate the appropriate spring constant for each of four springs used to mount the washing machine to meet the specification for the transmissibility and select them from the enclosed table.

b) This washing machine can also be run at a spin speed of 750 RPM. Calculate the Dynamic Force at this RPM and then calculate the force transmitted to the supporting structure compared with the force transmitted at 950 RPM.

c) Using the size of the rotating wash cylinder estimate the unbalanced weight in the maximum load.

MountNo. | MountConstant () | Load at2-in. Deflection () |

SLF-401 | 50 | 100 |

SLF-402 | 70 | 140 |

SLF-403 | 90 | 180 |

SLF-404 | 115 | 230 |

SLF-405 | 150 | 300 |

SLF-406 | 200 | 400 |

SLF-407 | 305 | 610 |

SLF-408 | 365 | 730 |

SLF-409 | 500 | 1,000 |

SLF-410 | 650 | 1,300 |

SLF-411 | 740 | 1,880 |

SLF-412 | 1,150 | 2,300 |

SLF-413 | 1,500 | 3,000 |

SLF-414 | 2,000 | 4,000 |

SLF-415 | 2,680 | 5,360 |

SLF-416 | 3,500 | 7,000 |

SLF-417 | 4,750 | 9,500 |

SLF-418 | 6,250 | 12,500 |

SLF-419 | 8,200 | 16,400 |

SLF-420 | 11,250 | 22,500 |

SLF-421 | 15,000 | 30,000 |

SLF-422 | 20,000 | 40,000 |

a)

Therefore,

Therefore, each spring has

Therefore, SLF 412 and 413 would not meet the transmissibility spec.

b)

Therefore slightly higher than at 950 RPM

c)

Perhaps a wet hand towel out of balance

The commercial dryer was considered earlier under the assumption that the vibration was only a vertical mode. However the excitation die to rotating imbalance would cause excitation both horizontally and vertically that would excite coupled modes as well. As a result it is necessary to consider all the planar modes of vibration to investigate if they are near the excitation frequency. The diagram of the dryer is shown below. The lateral stiffness of the supporting springs is estimated to be 0.8 of the vertical stiffness. The excitation frequency is 950 rpm.

a) Specialize the equations of motion for the case shown in the above diagram

b) Determine if any of the modes are uncoupled. For the coupled modes derive the equations of motion for free vibration starting from Newton’s laws.

c) Calculate the natural frequencies and mode shapes for the dryer.

d) If the supports for the dryer were redesigned to make loading easier with h = 0 (as illustrated) what are the natural frequencies compared to the original design. How do they compare with the excitation frequency?

a)

b) Therefore ,

The vertical mode is uncoupled

Therefore

Which matches the reduction of the general case where has been replaced by and by

c) Calculate the natural frequencies and mode shapes for the dryer

The vertical mode is as previously

which is well below the excitation frequency of 950RPM (15.8Hz)

Assume that

Therefore

Mode Shapes

*** QuickLaTeX cannot compile formula: \[(4\ell-mp^{2})X+4\ellh\Theta = 0\] *** Error message: Undefined control sequence \ellh. leading text: \[(4\ell-mp^{2})X+4\ellh

Therefore

*** QuickLaTeX cannot compile formula: \[\begin{split} \frac{X}{\Theta} &= \frac{-4\ellh}{(4\ell-mp^{2})} \\&= \frac{-h}{1-\frac{m}{4\ell}p^{2}} \end{split}\] *** Error message: Undefined control sequence \ellh. leading text: ...rac{-h}{1-\frac{m}{4\ell}p^{2}} \end{split}

d)

now all natural frequencies are all below the excitation frequency of 158Hz (950RPM)

Semi-Definite Systems – Rigid Body Modes

So far we have dealt with systems that have positive definite stiffness and mass matrices. This means the eigenvalues are positive & real. We can then define an orthonormal basic system. However consider the example

but the stiffness matrix is not positive definitive. It is singular. This is because one “natural frequency” is zero.

is called a zeroth mode

Therefore

but since is positive definite , therefore ,

and the system is simply moving off in one direction, let’s suppose we set:

This is the constraint equation:

or

Where this expression says the momentum of the 3 masses must be zero. This essentially says this is a 2 DOF system.

Note: If we tried to find a flexibility matrix it is not defined.

We can write the constraint in a compound form as:

We can now solve the eigenvalue problem using

so that

For the example above set ,

A second example of a 2DOF system is the vibration absorber as it can provide quite spectacular results when applied correctly. We will consider it first as a 2DOF undamped system, then include viscous damping. The closed form solution for the case of damping is much more involved.

#### Undamped Vibration Absorber

For S.S. solution, assume:

Under the steady state assumption, the equations of motion become:

Therefore, from :

And then from :

Therefore:

For use in vibration absorber applications, we define:

Then divide numerator and denominator by , .

For the 2DOF system shown with the forcing appiled to the mass , the response of the masses and is:

It is seen that from that if or , then the mass is stationary. This leads to the idea that the motion of can theoretically be reduced to zero if we employ , at the operating frequency . For this application. the subsystem , is called a vibration absorber and at the operating frequency:

For this application it is useful to modify and to reflect what we can call – the original system , , the absorber system , using the definitions:

Using those definitions show that and become:

And when :

This shows that when the subsystem , is timed so that is motionless as the force in the spring exactly opposes the excitation force .

The natural frequencies of the combined system are found from setting the denominator of the responses , to zero. Therefore, the natural frequencies are found is replaced by and the expression becomes:

Therefore:

If , this becomes:

If instead of the mass being excited by the force , the exication is through the base as , what is the difference in the response of and of masses and ?

The mass ratio, , essentially determines the spread of the two resonances:

()

Therefore, the equations are identical except that is replaced by . The response of and for the vibration absorber is:

Again when :

Therefore, so that at any time:

And the net dynamic force on is zero!

To solve the system with damping it is more advantageous to use complex numbers:

Or:

There are a number of special cases of the general forced damped analysis that are useful.

#### Damped Forced 2DOF System

Consider a general damped system that can be specialized for special applications.

Thus:

Therefore:

And we wish to find and :

Where . Thus:

There is a lot hidden in these symbols so we can consider some special cases.

##### 1. Untuned Viscous Vibration Absorber

###### (Houdaille damper or viscous Lanchester damper)

If we look at our general solution then:

Therefore:

Put :

Therefore:

We are really interested in the amplitudes. Use the result that if:

We wish to write .

Thus:

Where . Thus:

We wish to find the optimum ratio which minimizes this amplitude;

Multiply by . Therefore:

Therefore:

Therefore, the minimum amplitude for optimum occurs when . To find the value of , use this value of .

The amplitude at this value of is independent of !

Multiply by . Then:

Now evaluate this for .

P2

Note that:

When :

When :

When :

When :

When :

###### Example

A second very practical application is the viciously damped vibration absorber.

###### Solution

From the general solution:

Therefore:

If we set :

In order to apply this to various situations, it is useful to non-dimensionalize these relationships similar to that of the undamped vibration absorber.

Set:

NOTE:

Therefore:

###### CHECK

- If :
- If , then

Note that only when

If :

As a result somewhere between and , will give an optimal result.

Note that all the curves intersect at and . That is these points are independent of the damping. If we calculate their location, our problem is essentially solved as we must find the curve that passes through them with a horizontal tangent. Also by changing the “tuning” the two points can be shifted up and down the curve. We can do this until the two points have the same height and the horizontal tangent through one of them.

To find the two points where the amplitude is independent of the damping and this is independent when :

This leads to:

The negative gives while the positive gives .

The two answers for the forcing frequency ratio are functions of and but not of the damping .

We now wish to adjust so that & are equal. Since their magnitude is independent of the damping, chose then the magnitude:

For and (the two roots), we want the same magnitude:

Therefore:

Now it is not necessary to solve for and if we remember the sum of the roots in a quadratic in the negative coefficient of the middle term.

Therefore:

Therefore:

Amplitudes of the main mass for various values of absorber damping. The absorber is twenty times as small as the main machine and is tuned to the same frequency. All curves pass through the fixed points and .

Resonance curves for the motion of the main mass fitted with the most favorably tuned vibration-absorber system of one-forth of the size of the main machine.

(a) Peak amplitudes of the main mass as a function of the ratio for various absorbers attached to the main mass. (b) Peak relative amplitudes between the masses and for various absorbers. (c) Damping constants required for the most favorable operation of the absorber, , for obtaining the results of (a) and (b).

Curve 1 for the most favorably tuned and damped absorber; curve 2 for the most favorably damped absorber tuned to the frequency of the main system; curve 3 of the most favorably damped viscous Lanchester damper; curve 4 for the most favorably damped Coulomb Lanchester damper.

and using one of the roots:

while the optimum damping is given by:

### MTS Tuned MASS Damper Systems

The purpose of a Tuned Mass Damper is to reduce wind excited periodic motion by increasing the apparent damping of structures with low inherent damping and predominant fundamental vibration mode response to wind excitation. MTS Tuned Mass Damper systems are operated essentially as classical Tuned Mass Dampers but are unique in their adaptability, long-term stability, and size. They are designed to minimize energy consumption and offer the structural designer a powerful new means of reducing and controlling dynamic motion.

Although the primary building structure would not depend on the TMD operation for integrity, operation during high wind conditions will help minimize the possibilities of minor non-structural cosmetic damage to building interior components, such as partitions and walls. The TMD system will reduce building motion(induced by moderate and high winds) to approximately one-half of what the building motion would have been without the system operating. If an extremely high gust of wind should impact the building and TMD system can automatically shift its damping coefficient so that the mass does not overstroke. Whenever the system is operating, electronic and mechanical interlock devices ensure that the system is always in damping mode, attempting to reduce building motion, System shut down will occur if an undesirable operating condition should ever develop.

#### Understanding the Tuned Mass Damper

Understanding the tuned mass damper begins with the basic concept of mechanical vibration: some object is simply moving back and forth in a regular manner. in all cases some physical mass, having inertia, is moving back and forth against a resistance, or restoring force. These two basic ingredients: inertia and restoring force are necessary to all vibrations of common physical objects; in fact,. the presence of both guarantees that vibrations can take place.

One of the simplest possible examples of vibration is a weight, hanging under gravity from a simple coiled spring (Figure 1). Pulled down slightly from the rest position and then released, the weight oscillates up and down. Another example, (Figure 2) is the same weight mounted on the elastic post. Pulled to one side and released, the weight will oscillate back and forth laterally.

Everything which has both inertia and sporing qualities can thus be expected to exhibit vibration tendencies if once started or excited in some way. Structures such as tall building are exposed to gusty winds and as a result they swap back and forth, or vibrate. Some structures vibrate slowly(i.e., at low frequencies) while others tend to vibrate more rapidly (i.e., at high frequencies).

###### Example

Consider a multi-story building(Figure 3) having a structural steel frame. Such a building fits all the conditions for vibration. and it will in fact tend to sway very slightly in the wind, moving slowly, requiring perhaps several seconds to accomplish a full swing over and back.

Tall steel freamed building occasionally vibrate enough to disturb their occupants, or at least to draw their attention to the fact that the building is moving or swaying. Actually, buildings oscillate to a much lesser extent then such common structures as aircraft wings and fuselages; yet, because, psychologically, buildings are”suppose” to stand still, occupant objection to building oscillation is often much greater than it is to aircraft vibratory motion. As a result, tuned mass damped systems may be used to lessen building vibration when it occurs.

Other methods may be used to reduce building vibration, but tend to be more costly, or are less certain and reliable. Two possible alternatives are to increase the building mass while maintaining its natural vibration frequency approximately constant, or to increase building damping by adding special energy-absorbing devices or using masonry filler walls.

To achieve the same degree of motion reduction as TMD would require increasing the building weight and lateral stiffness(the structure’s resistance to lateral sway, measured in pounds per inch for example) by a factor of 1 2/3 to as much as 3. The added weight requires heavier structural members and foundations. and the added lateral stiffness requires heavier structure, all of which tend to be very costly. It’s inefficient to invest in very costly large increases in building weight and structure when these increases are unnecessary for safety and are done in order to increase occupant comfort. Methods which reduce building motion by increasing damping are more attractive.

Damping resistance to an vibration is simply any force which tends to dissipate or drain away the energy of that vibration. Damping devices are thus of a type tending to produce friction or bring other energy-absorbing mechanisms into place. For example, the shock absorbers on an automobile dissipate oscillatory energy by hydraulic means. Some tall buildings possess natural damping in the form of masonry cladding walls (the Empire State Building is an example). The inherent friction in masonry construction (or between masonry and steel) is often an effective oscillation damper. Modern steel-and-glass skyscrapers, on the other hand, possess relatively little natural damping.

The tuned mass damper represents an alternative to the damping devices mentioned above. TMD devices have been successfully used to arrest vibrations of power line cables, bridge, and highway-sign structural members, chimneys, tall towers, and antennas, and are used extensively in the automotive and aircraft industry. A TMD does not really act as a damper through energy drainage; instead, it performs in a “reactive” way, temporarily “bouncing” away from the mass it controls. The principles of the tuned mass damper as discussed in the following paragraphs.

###### Example

Consider the single spring-mass oscillator depicted in Figure 1. The tuned mass damper consists basically of another mass ( the secondary mass) attached by a spring to the first () as suggested in Figure 4. The secondary mass is “tuned” to the first in the following sense: The stiffness of spring and the value of the mass are adjusted so that the natural vibration frequency of mass when is held rigidly still, is the same as that of mass when is absent entirely. Clearly, the sense of the word “tuned” in this context is that one vibration frequency(that of the secondary mass alone) is brought into coordination with another frequency (that of the first mass alone).

The effect of the presence and tuning of the secondary mass of the TMD is usually dramatic upon the action of the first mass. It tends to greatly reduce or suppress (i.e., “damp”) the vibratory motion of the primary mass. The exact conditions of this action require, and are susceptible to, extensive mathematical analysis.

Suppose a primary mass is excited by some repetitive force that recurs at somewhere near the natural frequency of the primary system. The mass responds by oscillating at large amplitudes, i.e.,”resonating”. Again, suppose a randomly varying force, repetitive but not just occurring at one single frequency, excites the primary mass in a “broad spectrum” of way (i.e., at a mixture of oscillations of many frequencies simultaneously). In either of these cases, the primary mass oscillates manly at its own natural frequency; in the randomly excited case it will also take on random amplitudes.

The action of a tuned secondary mass in either of the above cases is to “take upon itself” the energy of the initial motion of the primary mass and to drastically reduce the response of the latter. This means that any exciting force which acts upon the primary mass must begin over and over again to build up the response of the latter from a low point, since each time that the response of this mass has been even temporarily excited, the TMD reacts against it and suppresses it. As a result there is much less occasion for the primary mass ever to be excited to a condition of excessive displacement, as it might be without the TMD.

The amount of secondary mass needed to accomplish the described effect can be variable, but it has been found that even if it is as small as on percent of the primary mass, it can still produce rather dramatic effects.

When a large, spatially-distributed primary mass consists of a building, like a skyscraper, the entire building may be represented mathematically as a single spring-mass vibrating system, considering not its entire mass as the primary oscillatory element, but only some equivalent(generalized) mass as the necessary element to be dealt with. This is to say that, since part of the building(for example, at the foundation) moves very little, while other parts(say at the top) move much more, the net effective moving mass has a value somewhere between the full mass and zero. This generalized mass also has a location or height coordinate – usually near the top of the building. The magnitude of the generalized mass may then typically range from down to or even less of the total mass of the building.

The actual amount of generalized mass to be considered depends upon just how the building vibrates – in technical terms. the modes in which it sways. Figure 5 depicts the amount of building sway at maximum travel in typically the first mode. The first mode of tall building sway consists of almost pure “lean” of the building all to one side, then all to the opposite side. Torsional or building twist modes can also occur. In these, the building rotates about an axis more or less coinciding with a straight, vertical line running through the geometric center of each story. Modes higher than first and second can occur in tall buildings but in almost practical cases responses in modes other than the first mode are less important.

In connection with a tall building, a TMD will be aimed principally at suppressing building sway or twist in “early” modes, say first in sway and first in twist. A possibility, if there is space, is to place inside , as suggested in Figure 6. This is precisely the kind of situation which exists in a tall building, where a floor may contain the TMD ().

Typically, TMD installations would be a high upper floor (refer to Figure 5) to work against the first swat mode(and twist, if important). The best choice of TMD placement is near the location of the expected highest amplitude of the mode to be suppressed.

The displacement of the secondary mass of the TMD will, when the device is working best, be greater by a factor of three or more depending on TMD design than the displacement of the primary mass. Thus, adequate space and other provisions must be provided for such motion. A second point is that some local damping provided to mass “broadens” its ability to effectively respond to excitations not precisely at the frequency of tuning. The amount of this local TMD damping also controls the relative displacement of the TMD mass . The relative displacement may be reduced by increasing TMD local damping, at the expense of somewhat reduced TMD effectiveness.

The net result of the TMD installation is that the primary mass will respond to wind as if it were damped to an amount at least twice as much as its actual inherent damping. For exmaple, if actual building damping is , a TMD with mass of the primary (effective) mass can increase total building damping in the first mode to about while a TMD with mass of the primary mass can increase total damping in the first mode to over .

For one of the first TMD systems in a tall office building it was estimated that, to achieve the same reduction of building wind motion, the cost used a TMD system was less than one third of the cost of a conventional approach involving increased building mass and lateral stiffness.

#### System Description

###### Principle of Operation:

The MTS Tuned Mass Damper Systems use a moderate force servohydraulic actuator to control the position of a passively resonant spring-mass system. The servoactuator is sized and electronically programmed so that the exact damping and resonant frequency can be electronically set to assure linearity, long term stability, and optimum system characteristics for the particular input motion experienced.

The Tuned Mass Damper systems combine the damper and the servo-hydraulic drive into one servo valve/actuator combination. The system uses substantial flow rates of oil, but requires only a moderate supply pressure minimizing energy consumption while permitting sufficient active actuator forcing capability for a