Hair Settings PDF

ME 563
Mechanical Vibrations
Lecture #12
Multiple Degree of Freedom
Free Response + MATLAB
Free Response
We can solve for the homogeneous solution to a coupled set of

equations in a multiple degree of freedom linear system by:
- Identifying the initial conditions on all the states
- Assuming a solution of the form {x(t)}={A}est
What does this last assumption imply about the response?
Two DOF System

Consider the two degree of freedom system of equations:
If we make a solution of the form,

the single DOF case, we obtain:
, as we did for
Non-trivial solutions satisfy:
( Ms
+ 2Cs + 2K ) ( Ms + Cs + K ) (Cs + K ) = 0
2
Characteristic Equation
There are four solutions that satisfy the characteristic
equation and these solutions are expressed as follows when the
modal frequencies are underdamped:
s1=1+j1, s2=1-j1, s3=2+j2, and s4=2-j2

Real part: decay rate
Imaginary part: oscillation frequency
The solution to the homogeneous equation is then written as

follows (for the first two complex conjugate roots):
Free Response Form

The free response is usually written in the following form for
a multiple degree of freedom system:
Modal vector
(can be scaled)
Decaying co-sinusoid
(common to both degrees of freedom)
Four constants Four initial conditions are required.
Example Solution
For K=1 N/m, C=0.1 Ns/m, and M=1 kg, the solution is given
by:
MATLAB Solution
To obtain solutions for the free response in MATLAB, the
procedure explained before for a single DOF system is used:
Or
MATLAB Solution
We use the following formulation to define the outputs, {y}, if
we are only interested in the displacement variables:
Or:
{ y}21 = [C ] 24 {q}41 + [ D] 22 { f (t)}21
Free Response (Eigen-Analysis)
We can also solve the homogeneous equations of motion by:

- Identifying the initial conditions on all the states
- Identifying the modal frequencies, s, and vectors, {X}, using
eigen-analysis.
Do you remember what an eigenvalue problem looks like?
Eigenvalue Problem
Consider the case when we have no damping:
We can write this set of equations of motion as follows:
Now assume a solution of the form {x(t)}={A}est
10
Eigenvalue Problem
Consider the case when we have no damping:
This last equation is the standard eigenvalue problem.

Matlab solves this equation with the command eig(A):
[v,d]= eig(A);
11
Matlab eig
How do we interpret what Matlab gives us with this
command?
d is a diagonal matrix of -s1,22

v is a matrix of modal vectors
d=
2.6180
0
0 0.3820
v=
-0.8507 -0.5257
0.5257 -0.8507
How do we find the modal frequencies and vectors?
12
Matlab eig
First, we take the negative square root of the diagonal entries
of the d matrix providing the modal frequencies in rad/s:
-d =
0+1.6180i
0
0
0+0.6181i
Second, we can scale the columns of v to produce easy to

interpret modal vectors:
v=
1.0000 0.6180
-0.6180 1.0000
13
Matlab Form of Response

The form of the free response can then be written in the form:
In summary, the eigenvalue problem for an undamped system

is given by:
Where does the damping fit?
14
Damped Eigenvalue Problem

To obtain solutions for the free response in a damped system,
the state variable form of the equations of motion are used:
and then the eigenvalues and eigenvectors of the state matrix

are calculated using eig. This approach works because the
assumed solution {q}est is also used for the 1st order system:
[v,d]= eig(A);
15
Matlab eig
How do we interpret what Matlab gives us with this
command in the first order form of the state variable model?
d is a diagonal matrix of modal frequencies s1,2 and s3,4

v (1:2,1:4) is a matrix of modal vectors {X}1 and {X*}1
{X}2 and {X*}2
diag(d) =
-0.1309 + 1.6127i s3,4
-0.1309 - 1.6127i
-0.0191 + 0.6177i s1,2
-0.0191 - 0.6177i
15
Matlab eig
Eigenvectors are not the modal vectors of the 2nd order system
eigenvectors involve both position and velocity states. We
look to the position states to obtain the scaled modal vectors.
v=[
-0.0362 - 0.4457i
0.0224 + 0.2755i
0.7236
-0.4472 + 0.0000i
-0.0362 + 0.4457i
0.0224 - 0.2755i
0.7236
-0.4472 - 0.0000i
-0.4472 + 0.0000i
-0.7236
0.0085 - 0.2763i
0.0138 - 0.4470i
-0.4472 - 0.0000i
-0.7236
0.0085 + 0.2763i
0.0138 + 0.4470i
16
Damping in MDOF Systems

Note that for the previous example, the damping matrix [C]
corresponded to a proportionally viscously damped system:
For this type of damping, note that the modal vectors were
entirely real, i.e., the phase between DOFs was either 0 or 180o.
16

What if the system is not proportionally damped? The modal
vectors in this case are complex but the response is not:
For this type of damping, note that the modal vectors are
complex, so the phase between the DOFs is not 0 or 180o.
What does this mean?
17

What if the system is not proportionally damped? The modal
vectors in this case are complex but the response is not:
Imag
Mode #1
Mode #2
Real
2nd order EOMS

cannot be uncoupled
using real modal
vectors!

Hair Settings PDF

Uploaded by

Copyright:

Available Formats

Hair Settings PDF

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Hair Settings PDF

Uploaded by

Copyright:

Available Formats

ME 563

We can solve for the homogeneous solution to a coupled set of

Two DOF System

If we make a solution of the form,

Non-trivial solutions satisfy:

s1=1+j1, s2=1-j1, s3=2+j2, and s4=2-j2

Imaginary part: oscillation frequency

The solution to the homogeneous equation is then written as

Free Response Form

Four constants Four initial conditions are required.

{ y}21 = [C ] 24 {q}41 + [ D] 22 { f (t)}21

Free Response (Eigen-Analysis)

We can also solve the homogeneous equations of motion by:

We can write this set of equations of motion as follows:

Now assume a solution of the form {x(t)}={A}est

This last equation is the standard eigenvalue problem.

d is a diagonal matrix of -s1,22

How do we find the modal frequencies and vectors?

Second, we can scale the columns of v to produce easy to

Matlab Form of Response

In summary, the eigenvalue problem for an undamped system

Damped Eigenvalue Problem

and then the eigenvalues and eigenvectors of the state matrix

d is a diagonal matrix of modal frequencies s1,2 and s3,4

Damping in MDOF Systems

Damping in MDOF Systems

What does this mean?

Damping in MDOF Systems

2nd order EOMS

You might also like