## Approximate Methods for Continuous Systems: Application to Lateral Vibrations of Beams

The systems we have considered so far have been relatively simple so we were able to find analytical solutions. This is usually possible in only the simplest situations. For more complicated systems, say ones which include

- non-uniform material properties and/or cross sections, or
- more complicated boundary conditions,

it is very difficult, if not impossible, to determine analytical solutions. In these situations we usually resort to numerical solutions. Sometimes however, other approximate techniques can also be useful. In particular, Rayleigh’s Method is well suited for use with continuous systems.

Recall that for discrete systems, Rayleigh’s Quotient was expressed as

(11.1)

where is the stiffness matrix, is the mass matrix and {} is a mode shape vector. To determine an approximate natural frequency, a trial mode shape vector {} is used in 11.1 to return an upper limit on the lowest natural frequency.

To apply Rayleigh’s Quotient to continuous systems, we recall that Rayleigh’s Quotient is based on the idea that the maximum kinetic energy and the maximum potential energy in a conservative system will be equal. This requires that the potential energy be zero in the equilibrium configuration.

when all of the masses in the system undergo simple simultaneous harmonic motion, as happens in the modal responses.

We can also apply this idea to continuous systems. If is some representative displacement function, the procedure will be to

- Estimate the fundamental mode shape .
- Assume simple simultaneous harmonic motion of the form
- Determine the maximum potential energy .
- Determine the maximum kinetic energy
- Set and solve for

It is important here that the mode shape assumed in step (a) *satisfies the (geometric) boundary conditions* of the problem. If this is the case then the frequency determined by this procedure will again be an upper bound on the lowest natural frequency. If the assumed mode shape doesn’t satisfy the (geometric) boundary conditions in the problem, then we no longer have this guarantee.

To apply Rayleigh’s Method to beams, we need to be able to determine the potential and kinetic energies associated with a beam vibrating in a given assumed mode shape.

### Potential Energy

Potential energy is stored in a beam as elastic (strain) energy. In most cases the energy stored due to bending is much larger than that stored due to shear or axial forces. As a result only the energy associated with bending is usually considered.

Consider an initially straight beam, deformed to some arbitrary configuration (limited to small motions) during its motion as shown in Figure 11.1{a}. Figure 11.1{b} shows an infinitesimal element of the beam during this motion. At this point the beam element has a radius of curvature~ which is related to the included angle by

(11.2)

For an elastic beam the bending moment is related to the deflection (curvature) by

(11.3)

Now, consider the infinitesimal work done on the beam element as it is deformed from its initially unloaded straight (equilibrium) configuration to one with an included angle as shown in Figure 11.1{b}. The work is

which is equal to the area under the – graph as shown. For an elastic beam, this work done on the beam is equal to the strain energy stored in the beam element, so that we get

Using equation 11.2 and equation 11.3 this becomes

or

(11.4)

For small deflections,

so 11.4 becomes

(11.5)

which gives the strain energy stored in the infintesimal beam element. To find the total strain energy stored, we integrate over the length of the beam

(11.6)

If we make the usual assumption that the motion of the beam is of the form

then

so equation 11.6 becomes

(11.7)

The maximum strain energy occurs when

which occurs at the same instant for every point on the beam. Using 11.7 we find that

(11.8)

Equation 11.8 represents the maximum strain energy stored in the beam for a given (assumed) mode shape . Note that since

Equation 11.8 can equivalently be written as

(11.9)

which may be easier to use in some cases.

### Kinetic Energy

Once again consider the infinitesimal beam element shown in Figure 11.1(b). The kinetic energy of that element is given by

where

is the mass of the infinitesimal element and its speed is given by (This considers only the translational kinetic energy and ignores any rotational energy, i.e. Euler beam theory.)

Assuming once again that the motion of the beam is of the form

(11.10)

Therefore the kinetic energy becomes

The total kinetic energy in the beam is then obtained by integration

This is a maximum when

which occurs simultaneously at all points on the beam. Therefore the maximum kinetic energy in the vibrating beam is

(11.11)

### Rayleigh’s Quotient For Lateral Vibrations of Beams

As we have seen previously, for a conservative system,

Using equation 11.8 and equation 11.11 this becomes

(11.12)

Equation 11.12 is Rayleigh’s Quotient for beam applications where, as in the multiple degree of freedom cases discussed previously, is an assumed mode shape for the beam. If the assumed mode shape satisfies the geometric boundary conditions in the system, the natural frequency obtained using equation 11.12 will be an upper bound on the lowest (fundamental) natural frequency.

#### Example

- Estimate the fundamental natural frequency of a uniform cantilever beam of length ( are all constants). Assume a mode shape of
- Repeat the calculation in part (1) but instead use an assumed shape function
which is the static equilibrium shape of the beam when it is subjected to a unit tip load. Compare these predictions to the true fundamental natural frequency

### Extensions to Rayleigh’s Method for Beams

For beams, Rayleigh’s Quotient, given by 11.12 was obtained by considering the ratio of the maximum potential and kinetic energies in the system. We can extend this idea to consider additional springs and masses in the system.

Consider the system shown in Figure 11.2. To determine the potential energy stored in the system for an assumed mode shape , we can include the potential energy stored in the springs as additional terms

where and are the deflections of the beam at locations and respectively. The maximum potential energy stored in the systems then becomes

(11.13)

Similarly, the extra kinetic energy in the beam due to the additional masses is given by

where

However, we have seen in equation 11.10 that

so that the speeds become

These will be maximum when so that the maximum additional kinetic energy is

As a result the total maximum kinetic energy in the system is given by

(11.14)

Equating 11.13 and equation 11.14 we see that

or

This can be generalized to a beam with masses and springs as

(11.15)

#### Example 1

A driveshaft for a marine propeller can be modeled as shown below.

- Estimate the fundamental natural frequency using Rayleigh’s quotient. Assume a mode shape of
and

*neglect the mass of the shaft*. - An alternative design for this propeller features a uniform shaft stiffened with an intermediate bearing (with stiffness ) as shown below.

Estimate the location for this bearing such that the natural frequency of the alternate design is the same as for the original design. Express your answer in terms of and .

Use the same assumed mode shape and once again neglect the mass of the shaft in your calculations.

#### Example 2

Two views of an airplane wing are shown in which the wing is modeled as a cantilever that has a varying stiffness and a wing tank of mass when filled with aviation gasoline.

- Determine an expression for the lowest natural frequency of free vertical vibrations of the wing and fuel tank. Assume a trial function
- If the natural frequency is not to be more than 25% greater when the tank is empty () than when it is full (mass is included), determine the maximum mass which can be added in comparison to the mass of the wing itself.