Me'Scopeves Application Note #7: Modal Analysis of A Mass-Spring-Damper System
Me'Scopeves Application Note #7: Modal Analysis of A Mass-Spring-Damper System
www.vibetech.com
1/14/02
BACKGROUND MATH
The time domain equation of motion for the mass-springdamper is represented by Newtons Second Law,
(1)
where:
M = mass value
C = damping coefficient
K = spring stiffness
&x&(t ) = acceleration
x& (t ) = velocity
x(t ) = displaceme nt
f (t ) = excitation force
Laplace Transforms
By taking Laplace transforms of the terms in Equation 1 and
setting initial conditions to zero, an equivalent frequency
domain equation of motion results,
[M s 2 + C s + K ] X(s ) = F(s )
Figure 1. Mass-Spring-Damper.
The purpose of this Application Note is to review the details
of modal analysis, to provide a better understanding of the
modal properties of all structures. The modal properties of
real world structures are analyzed using a multi-degree-offreedom (MDOF) dynamic model, whereas the model used
here is a single degree-of-freedom (SDOF) model. Nevertheless, the dynamics of MDOF structures are better understood
by analyzing the dynamics of this SDOF structure.
Modes are defined for structures, the dynamics of which can
be represented by linear ordinary differential equations like
Equation 1. The dynamic behavior of the mass-springdamper structure in Figure 1 is represented by a single (scalar) equation, Equation 1. (An MDOF structure is represented by multiple equations like Equation 1, which are written in matrix form.)
Because of the superposition property of linear systems, the
dynamics of an MDOF structure can be written as a summation of contributions due to each of its modes. Each mode
where:
(2)
Transfer Function
Equation 2 can be rewritten by simply dividing both sides by
the coefficients of the left-hand side.
1
F(s )
X(s ) =
2
Ms + C s + K
(3)
H (s ) =
Page 1 of 10
X(s )
1
=
2
F(s ) M s + C s + K
(4)
App Note #7
1/14/02
tions of the poles in the s-Plane have also been given some
other commonly used names, as shown in Figure 3.
Modal Parameters
The coordinates of the poles in the s-Plane are also modal
parameters. Rewriting Equation 4 in terms of its pole locations, or modal parameters,
1/ M
H (s ) = 2
2
s + 2 s +
0
0
where:
By inspection of Equation 4, it is clear that the Transfer Function goes to infinity for values on the s-Plane where its denominator is zero. It is also clear that as s goes to infinity,
the Transfer Function will go to zero.
The denominator is a second order polynomial in the S variable, called the characteristic polynomial. Since it is a second order, it has two roots (values of s) for which it will be
zero. These two roots of the denominator are called the poles
of the Transfer Function. Furthermore, the poles are complex
conjugates of one another. The poles therefore, are the locations in the s-plane where the Transfer Function has a value
of infinity. The poles are also called eigenvalues.
p 0 = 0 + j0 , p 0 = 0 j0
(5)
S-Plane Nomenclature
The real axis in the s-Plane is called the damping axis, and
the imaginary axis is called the frequency axis. The loca-
(6)
C
K
2
, 0 =
2M
M
0 = mod al damping coefficient
0 =
(7)
0 = 0 + 0
0 =
( 0 ) is written as,
C
0
=
0 2 M K
(8)
Page 2 of 10
App Note #7
1/14/02
FRF = H( j) = H(s ) s = j
=
X(s )
X( j)
=
F(s ) s = j F( j)
(9)
p 0 = 0 + j 0 & p 0 = 0 j 0 .
Equation 11 also says that the FRF for an SDOF is fully
represented by two poles and two residues. Furthermore,
since the residues are equal and the poles are complex conjugates of one another, the complete dynamics of an SDOF
0 ) ,
system is fully represented by a modal frequency ( j
modal damping
Mode Shapes
One final step is to represent the FRF in terms of mode
shapes instead of residues,
H ( j) =
{u 0 }2
1 {u 0 }2
2 j j p 0 j p 0
(13)
where:
H ( j) =
(1 / M )
( j) + 2 0 j + 0
2
(10)
Furthermore, using the poles of the Transfer Function, a partial fraction expansion can be performed on Equation 10
yielding,
R0
1 R0
H ( j) =
2 j j p 0 j p 0
where:
R 0 = 1 / 0 M
(14)
{u 0 } .
(11)
Page 3 of 10
App Note #7
1/14/02
Using Equation 11 to express the FRF in terms of modal parameters, the IRF is written,
1 R0
R0
h(t ) = FFT 1
2 j j p
j
p
0
0
or
h( t ) =
T
R 0 e p0t R 0e p 0 t
2j
where:
(16)
(17)
(sin(0t + 0 )) .
0 = angle of R 0
sinusoidal function
(15)
or
h(t ) = T R 0 e 0t (sin( 0 t + 0 ))
A NUMERICAL EXAMPLE
We will build a dynamic model for a mass-spring-damper
system using MEscopeVES. The Visual Modal, Modal Pro
or Visual SDM option is required to carry out all steps of this
exercise. If you have the Visual SDM option, you can add
mass, stiffness, and damping elements directly to the SDOF
model and generate its modal parameters. Otherwise, you
will have to enter the modal parameters into a Shape Table.
The mass, spring, & damper element values must have units
associated with them. Both English and Metric units are included below. Take your pick.
Execute File | Options in the Structure window, and select the following units.
Metric Units:
English Units:
Page 4 of 10
App Note #7
1/14/02
Click near the top Point to number it 1, and then near the
bottom Point to number it 2.
Metric Units:
1 N = 1 kg x 1(m/sec2), or
1 kg = 1 N-sec2/m
English Units:
This will add the mass element to the model and create a
new row in the mass properties spreadsheet.
On the Properties tab, change the Length & Width parameters to 3, and change their Points parameters to 2.
Drag the vertical blue bar to the left, and enter the
mass value (either 386.09 lbm or 1 kg) into the spreadsheet.
Page 5 of 10
App Note #7
1/14/02
Notice in Figure 9, that the measurement type is Unit Modal Mass mode shape, and that its magnitude value is 1.
, and click on
This will add a spring element to the model, and create a new
row in the spring properties spreadsheet.
Click on the Measurement Type cell heading. A dialog box will open. Select Residue Mode Shape from
the list, and click on OK.
Page 6 of 10
(m/N-sec)
App Note #7
1/14/02
(in/lbf-sec)
C = 2 0 M K = 2(0.05) (20) 2 = 2
R 0 = 1 / 0 M = 1 / 20 = 0.015915
, and
This will add the damper element to the model, and create a
new row in the damper properties spreadsheet, as shown
below.
{u 0 }M{u 0 } = 1 => A =
1
0
1
1
{u 0 } = M =
0 0
To scale the mode shape to unit modal masses,
The Shape Table will open with the damped mode in it.
Notice that the damped modal frequency is now slightly less
than 10 Hz. We can check this with the formula for damped
natural frequency.
0 =
If you have the Visual SDM option, you can add a damper to
the 3D model in the Structure window, and compute a new
damped mode of vibration. Equation 8 is the relationship
between Percent of Critical damping and the damper coefficient (C). Solving for this coefficient gives,
Page 7 of 10
C
2
rad
=
=
2M
2
sec
App Note #7
1/14/02
0 = 0 0
= ( 20 ) 2 2
rad
sec
= 9.9875 Hz
= 399
Stiffness Line
As the frequency approaches 0 Hz, Equation (4) shows that
the (displacement/force) FRF approaches a constant, equal
to the inverse of the stiffness. (This is also called the flexibility.) By letting S=0 in Equation (4), the flexibility is,
H ( 0)
Page 8 of 10
1
1
=
= 0.0002533
K 3947.8
App Note #7
1/14/02
Mass Line
As frequency becomes large, the FRF becomes dominated by
the 2 term in its denominator (characteristic polynomial).
(See Equation (10).) Therefore, at high frequencies the magnitude of the FRF can be approximated by,
H()
2
M
We can use the Tools | Math | Differentiate command twice
to multiply the (displacement/force) FRF by 2. Hence, its
magnitude for high frequencies will become a constant,
1
2 H()
M
Execute Modes | Modal Parameters in the (displacement/force) FRF Data Block window. The Curve Fitting panel will be displayed.
To check this,
Check the mass, damping & stiffness columns to display them in the spreadsheet, and click on OK.
Spread the edges of the Band cursor to enclose the modal peak at 10 Hz.
NOTE: In Application Note #16, Integration & Differentiation of FRFs, it is shown that multiplying a (displacement/force) FRF by 2 does not yield the correct (acceleration/force) FRF. Nevertheless, it is still useful for finding the
mass line of an SDOF system as done here.
The modal parameter estimates are displayed in the spreadsheet, as shown in Figure 16.
Page 9 of 10
App Note #7
1/14/02
(T R 0 e
0t
(e
0t
max . mag . = T R 0 e 0 t
t =0
= T R0
This is the product of the time length (T) of the Trace, and
the magnitude of the residue
R0 .
f =
1
= 1 Hz
T
Page 10 of 10