Help Contents | Module Library
Spring-Mass System
An EngineersToolbox Calculation Module

Summary

The Spring-Mass System module calculates the natural frequencies, mode shapes, and the frequency response of a generic system of discrete springs, dampers, and masses.

 

Background

A frequency response analysis is a method used to compute structural response to steady-state harmonic excitation. Harmonic excitation is often encountered in engineering systems. For example, it is commonly produced by an unbalance in rotating machinery.

In a frequency response analysis, the excitation forces are explicitly defined in the frequency domain. The applied loads are implicitly sinusoidal in nature. In the simplest case, the loading is defined as having an amplitude at a specific frequency. The steady-state response occurs at the same frequency as the loading. If damping is present in the system, the response peak will not coincide with the peak loading. The shift in response is called a phase shift and is illustrated in Figure 1.

Figure 1. Phase shift between loading and steady-state response

The results obtained from a frequency response analysis include displacements, velocities, accelerations, and forces. The calculated values are complex numbers that can be represented by a magnitude and a phase angle, or by real and imaginary components.

In this module, the steady-state response of the system is computed at discrete excitation frequencies by solving a set of coupled matrix equations using complex algebra. The complex equation of motion for damped forced response is given by:

The steady-state harmonic response is:

Where {x(w)} is a complex displacement vector. Taking the first and second derivatives gives the associated velocity and acceleration vectors:

When the above expressions are substituted into the equation of motion, the following is obtained:

which simplifies to an equation of motion in the frequency domain:

This expression represents a system of equations with complex coefficients. The steady states response vector {x(w)} is computed by inserting the forcing frequency w and multiplying the inverse of the matrix on the left side of the equation by the vector of applied forces.

 

Input

The Spring-Mass System main input form is shown in Figure 2. The primary input sections are described in detail below.

Figure 2. Module input form

System Definition

The spring-mass system is defined using the numbering convention shown in Figure 3.  For each system with n input masses, there can be a maximum of n-1 springs and dampers.  Input values are entered using the ETBX spreadsheet panel shown in Figure 4. Adhere to the following rules when defining your dynamic system:

  • The spreadsheet line number will correspond to the numbering convention shown in Figure 3 and on the main input form.  A mass m2 will be entered on line 2, a spring rate k3 will be entered on line 3, and so on.
  • To constrain a mass, enter 'n' in the constraint column.
  • The number of mass input values defines the system order.  A blank cell in the mass column marks the end of input.
  • Mass input values must be positive real numbers.
  • Spring rate and damping input values must be non-negative real numbers.  Zero values are allowed.  Blank cells will be interpreted as a zero value.
  • The module will only recover output for degrees-of-freedom with a 'y' in the output column.
  • A maximum of 20 degrees-of-freedom are permitted.

Figure 3. Spring-Mass System

Figure 4. System definition input spreadsheet.

 

Force Input

The steady-state harmonic forcing function is defined in the frequency domain by the complex value:

The complex input force can be input using real and imaginary pairs, or in phasor form using the magnitude and phase angle.  ETBX input requirements for complex numbers are described in Table 1 below.

Real/Imaginary Input Format

Format: x + yi or x - yi

Description:

x is a real value representing the real part of the complex number

y is a real value representing the imaginary part of the complex number.

i is the imaginary unit equal to the square root of -1.  Note that ETBX does not recognize the variable j as the imaginary unit, although it is often used by engineers and physiscists.

Examples: 1.9+7.8i, 2.0-8.0i, 3i, -5.0i

Magnitude/Phase Input Format

Format: A | f

Description:

A is a real value representing the magnitude of the complex number.

f is a real value representing the phase angle, in degrees.  The phase angle corresponds to the counter-clockwise angle from the positive real axis, i.e., the value of f such that x = cos f and y = sin f, where x and y are the real and imaginary parts of the complex number.

A vertical bar character '|' is required as a separator between the magnitude and phase.

Examples: 9.0 | 180.0, 9.6 | 30.0, 1 | 30

Damping Input

The user can specify how numerical values entered in the damping column of the input spreadsheet will be interpreted.

If the user specifies direct damping, input values are interpreted as the viscous damping coefficients, Ci, of discrete dampers between degrees of freedom xi-1 and xi.  The form of the relationship is:

where Fi-1 and Fi are damping forces at degrees of freedom xi-1 and xi.

If the user specifies modal damping, input values are interpreted as the critical damping ratio, zi, used in the calculation of the generalized response for the ith mode in an uncoupled modal analysis:

Where, for the ith mode, di is the generalized displacement, mi is the generalized mass, ki is the generalized stiffness, wi is the undamped natural frequency, and fi is the generalized force.

Excitation Frequencies

A major consideration when you conduct a frequency response analysis is selecting the frequencies at which the solution is to be performed. The Excitation Frequencies input group provides two ways to define solution frequencies:

  • A list of discrete frequencies, separated by blanks spaces or commas.
  • A frequency range, defined by an initial frequency, a frequency increment, and the number of frequency increments.

Output Options

The important results obtained from a forced response analysis usually include displacements, velocities, accelerations, and forces. The Output Options group on the main panel lets you select the output type and format. The following types of output are available in either magnitude/phase or real/imaginary format:

  • Displacement
  • Velocity
  • Acceleration
  • Spring Force
  • Damper Force
  • Effective (Dynamic) Spring Rate
  • Impedance
  • Effective (Dynamic) Mass

 

Results

The Spring-Mass System module solves for natural frequencies, mode shapes, and the frequency response of a system of discrete springs, dampers, and masses.

Results Report

Figure 5 shows an ETBX results report generated by the Spring-Mass module. The report contains the following elements:

  • An input summary that documents the dynamic system that was analyzed
  • An eigenvalue summary that lists the calculated eigenvalues (l=wn2) and the undamped natural frequencies in radians per second (wn) and in cycles per second (fn).
  • A list of calculated eigenvectors for each mode
  • Frequency response calculations

Figure 6 shows an example of an eigenvalue summary for a dynamic system with 4 modes.

Frequency response output is in real/imaginary format or magnitude/phase format. Figure 7 shows frequency response output in magnitude/phase format. Remember that the module will only print output for degrees-of-freedom with a 'y' in the output column of the system input spreadsheet.

 

Figure 5. Module results report

 

Figure 6. Eigenvalue and eigenvector output

 

Figure 7. Frequency response magnitude/phase output

 

Frequency Response Plots

Frequency response output is plotted using the standard ETBX plot plot shown in Figure 8.

 

Figure 8. Frequency response plot

 

System Animations

The Spring-Mass module creates ETBX Animation Panels to display animated mode shapes and forced response shapes. Press the View System button in the System Definition input group to create a new animation panel.

Note that ETBX modules only pass the model and results data to the animation panel when the panel is created. The animation panel therefore only contains a 'snapshot' of the model and any associated results. After the animation panel is created it is NOT updated to reflect subsequent changes to the model or results calculations. This allows multiple models to be compared side-by-side, and avoids any confusion that might be caused by automatic updates. To display the latest updates to your model, simply create a new animation panel. For more details, please refer to the documentation on Animating Results.

 

Figure 9. Spring-Mass System Animation

 

References
Craig, R. R. (1981), Structural Dynamics, John Wiley & Sons (New York).
Thomson W.T., and Dahleh, M. D. (1997), Theory of Vibrations with Applications (5th Edition), Prentice Hall (New Jersey). 
Woods, R. L., Lawrence, K. L. (1997), Modeling and Simulation of Dynamic Systems, Prentice Hall (New Jersey).

Help Contents | Module Library
Copyright © 2001-2008 Engrasp, Inc.  All rights reserved.