Holooly Plus Logo

Question 8.1: A helicopter manufacturer is developing a code of practice f......

A helicopter manufacturer is developing a code of practice for the installation of hydraulic service pipes, and wishes to comply with a ‘minimum standard of integrity’ vibration standard, which defines, among other requirements, a sinusoidal acceleration level of ± 5.0 g from 50 to 500 Hz, at the pipe’s attachments to the aircraft structure. One investigation considers straight lengths of steel pipe, installed with spacing, L, between support centers, as shown in Fig. 8.1(a). Only the fundamental vibration mode of the pipe is considered, and this is taken as the exact first bending mode of a simply supported, uniform beam. It may be assumed that a separate investigation has shown the response in higher order modes to be negligible.

Note: The system of units used in this example is the ‘British’ lbf inch system, in which the acceleration due to gravity, g, is 386 in./s², and mass is expressed in lb \mathrm{in}^{-1}s², a unit sometimes known as the ‘mug’.

The properties of the pipe are as follows:
E = Young’s modulus for the material = 30 × 10^{6}   lbf/in.²
D = outer diameter = 0.3125 in.
d =inner diameter = 0.2625 in.
I = second moment of area of cross-section, = (\pi / 64) (D^{4} – d^{4}) = 0.2350 \times 10^{-3} \mathrm{in.}^{4}
μ = total mass of pipe per inch = 0.0214 × 10^{–3} mug/in. (including contained fluid).
The non-dimensional viscous damping coefficient, γ, can be taken as 0.02 of critical.

Find

(a) The vertical single-peak displacement at the center of the pipe span, relative to the supporting structure, if L = 17 in.;

(b) The maximum oscillatory stress in the pipe.

8.1
Step-by-Step
The 'Blue Check Mark' means that this solution was answered by an expert.
Learn more on how do we answer questions.

Part (a)
From Table 8.1, for the first normal bending mode of a simply supported uniform beam:

\omega_{1} = \beta^{2}_{1} \sqrt{\frac{EI}{\mu}}     \mathrm{and}     \beta_{1}L = \pi \\ \mathrm{therefore:} \\ \omega_{1} = \frac{\pi^{2}}{L^{2}} \sqrt{\frac{EI}{\mu}}                (A)

where ω_{1} is the natural frequency in rad/s, and E,  I,  L,  and  μ are given above.

Also from Table 8.1, the shape of the first mode, with i = 1, is given by y_{1} = \sin \beta_{1}x.

Since \beta_{1} = \pi/L,   \mathrm{then}  y_{1} = \sin (\pi x / L).   \mathrm{Defining}   y_{1} = y / y_{c} where y is the actual displacement at distance x along the pipe and y_{c} is the maximum displacement at the half-span position, then:

y = y_{c} \sin \left(\frac{\pi x }{L}\right)               0 < x < L                    (B)

which is sketched in Fig. 8.1 at (b). The slope, dy / dx, and the curvature, d²y / dx², are given by differentiating Eq. (B) twice with respect to x:

\frac{dy}{dx} = y_{c} \frac{\pi}{L} \cos \left(\frac{\pi x }{L} \right) ;               0 < x < L                    (C)\\ \frac{d^{2}y}{dx^{2}} = y_{c} \frac{\pi^{2}}{L^{2}} \left[- \sin \left(\frac{\pi x }{L} \right) \right]               0 < x < L                   (D)

These are sketched in Fig. 8.1 at (c) and (d).

We can now derive an equivalent system to describe the vibration of the pipe, lumping the mass, stiffness and damping at the center of the pipe, where the displacement relative to the base is y_{c}. With the base fixed, the kinetic energy, T, in the pipe is

T = \frac{1}{2} \mu \int_{0}^{L}{\dot{y}^{2}}dx                   (E)

Now differentiating Eq. (B) with respect to time:

\dot{y} = \dot{y}_{c} \sin \frac{\pi x}{L}                   (F)

which gives the vertical velocity at every point on the pipe.

Therefore:

T = \frac{1}{2} \mu \dot{y}_{c}^2 \int_{0}^{L}{\sin^{2} \frac{\pi x}{L}}dx = \frac{1}{4} \mu L \dot{y}_{c}^{2}                   (G)

This must be equal to the kinetic energy of the equivalent mass, \underline{m}, which is

T = \frac{1}{2} \underline{m} \dot{y}_{c}^{2}                     (H) \\ \mathrm{Thus,} \\ T = \frac{1}{4} \mu L \dot{y}_{c}^{2} = \frac{1}{2} \underline{m} \dot{y}_{c}^{2}  \\  \mathrm{and  so}  \\  \underline{m} = \frac{1}{2} \mu L                     (I)

The equivalent stiffness, \underline{k}, could be found in a similar way by equating the potential energy U in the pipe to that of the equivalent spring stiffness, \underline{k}, but a short cut is to use the fact that we already know the first natural frequency, ω_{1}, from Eq. (A):

\omega _{1} = \frac{\pi ^{2}}{L^{2}} \sqrt{\frac{EI}{\mu } } , \\  \mathrm{and  since}   \underline{k} = \underline{m} \omega ^{2}_{1} ,  \mathrm{then:} \\  \underline{k} = \underline{m} \omega ^{2}_{1} = \frac{1}{2} EI \frac{\pi ^{4}}{L^{3}}                    (J)

The pipe with damping and base motion can now be represented by the lumped model shown at Fig. 8.1(e). This is seen to be the same as that of the system shown in Fig. 2.3, and the equation of motion is given by Eq. (2.4). In the present notation this is

{m} \ddot{y} +{c} \dot{y} + {k} y = – {m} \ddot{x}                             (2.4)

\underline{m} \ddot{y}_{c} + \underline{c} \dot{y}_{c} + \underline{k} y_{c} = – \underline{m} \ddot{x}_{B}                             (K)

where \ddot{x}_{B} is the base acceleration input, in this case the acceleration at both pipe attachments.
The modulus of the frequency response function \left|y_{c}\right| / \left|x_{B}\right| of the equivalent system represented by Eq. (K) was given as Eq. (4.26):

\frac{\left|y\right|}{\left|x\right|} = \frac{\Omega^{2}}{\sqrt{(1 – \Omega^{2})^{2} + (2\gamma \Omega )^{2}} }                           (4.26) \\ \frac{\left|y_{c}\right|}{\left|x_{B}\right|} = \frac{\Omega^{2}}{\sqrt{(1 – \Omega^{2})^{2} + (2\gamma \Omega )^{2}} }                           (L)\\ \mathrm{where} \\ \Omega = \frac{\omega }{\omega _{1}} = \frac{f}{f_{1}};

and

ω is the excitation frequency in rad/s, ω_{1} the undamped natural frequency in rad/s, f the excitation frequency in Hz, f_{1} the undamped natural frequency in Hz and γ the viscous damping coefficient.
Also, from Eq. (4.32), but using the present notation:

\frac{\left|y\right| }{\left|\dot{x} \right| } = \frac{1}{\omega } \cdot \left(\frac{\left|y\right| }{\left|x\right| } \right) , \frac{\left|y\right| }{\left|\ddot{x} \right| } = \frac{1}{\omega^{2}} \cdot \left(\frac{\left|y\right| }{\left|x\right| } \right) , \frac{\left|\dot{y} \right| }{\left|x \right| } =\omega \cdot \left(\frac{\left|y\right| }{\left|x\right| } \right) , \frac{\left|\ddot{y} \right| }{\left|x \right| } =\omega^{2} \cdot \left(\frac{\left|y\right| }{\left|x\right| } \right)                          (4.32)

\frac{\left|y_{c}\right|}{\left|\ddot{x} _{B}\right|} = \frac{1}{\omega ^{2}} \cdot \left(\frac{\left|y_{c}\right|}{\left|x_{B}\right|}\right)                          (M)

From Eqs (L) and (M):

\frac{\left|y_{c}\right|}{\left|\ddot{x} _{B}\right|} = \frac{\Omega ^{2}}{\omega ^{2}\sqrt{(1 – \Omega^{2})^{2} + (2\gamma \Omega )^{2}} } = \frac{1}{\omega_{1}^{2}\sqrt{(1 – \Omega^{2})^{2} + (2\gamma \Omega )^{2}}}                 (N)

Equation (N) now gives the magnitude of the displacement, \left|y_{c}\right| at the center of the pipe in terms of the base acceleration magnitude, \left|\ddot{x} _{B}\right|. The maximum value of \left|y_{c}\right| will occur, for practical purposes, when the excitation frequency coincides with the natural frequency, i.e. when Ω = 1 . Therefore:

\left|y_{c}\right| _{MAX} = \frac{\left|\ddot{x} _{B}\right|}{2\gamma \omega^{2}_{1}}                   (O)

The natural frequency ω_{1}, in rad/s, is given by Eq. (A). Inserting numerical values:

\omega_{1} = \frac{\pi^{2}}{L^{2}} \sqrt{\frac{EI}{\mu}} = \frac{\pi^{2}}{17^{2}} \sqrt{\frac{(30 \times 10^{6}) \times (0.2350 \times 10^{-3})}{(0.0214 \times 10^{-3})}} \\ = 619.8   rad/s   (or   f_{1} = 98.6   Hz).

Referring to the vibration specification, the required input at this frequency is ± 5g, so \left|\ddot{x}_{B}\right| = (5 × 386) in. / s² = 1930 in./s² . Substituting these values and γ = 0.02 into Eq. (O), the required single-peak displacement amplitude at the center of the pipe is

\left|y_{c}\right|_{MAX} = \frac{\left|\ddot{x}_{B}\right|}{2\gamma \omega^{2}_{1}} = \frac{1930}{2 \times 0.02 \times (619.8)^{2}} = 0.126  in.

Part (b)
From Eq. (D), the curvature of the pipe is given by:

\frac{d^{2}y}{dx^{2}} = y_{c} \frac{\pi^{2}}{L^{2}} \left[- \sin \left(\frac{\pi x}{L} \right) \right].

The curvature at the mid point, (d^{2}y / dx^{2})_{c} , where x = L/2, is

\left(\frac{d^{2}y}{dx^{2}}\right)_{c} = y_{c} \frac{\pi^{2}}{L^{2}} \left(- \sin \frac{\pi}{2} \right) = – y_{c} \frac{\pi^{2}}{L^{2}}                        (P)

The corresponding bending moment is M_{c} = EI (d^{2}y / dx^{2})_{c} and the corresponding stress, s_{c}, at the center of the span and distance D/2 from the bending neutral axis, is

s_{c} = M_{c} \frac{D}{2I} = \frac{DE}{2} \left(\frac{d^{2}y}{dx^{2}}\right)_{c} = – \frac{\pi^{2}DE}{2L^{2}} y_{c}                        (Q)

The negative sign in Eq. (Q) is of no importance, since it is arbitrary whether tensile or compressive stress is taken as positive, so we can simply say that

\left|s_{c}\right|_{MAX} = \frac{\pi^{2}DE}{2L^{2}} \left|y_{c}\right|_{MAX}                        (R)

Inserting numerical values into Eq. (R), D = 0.3125 in.; E = 30 × 10^{6}   lbf/in.², L = 17 in. and \left|y_{c}\right|_{MAX} = 0.126   in., gives

\left|s_{c}\right|_{MAX} = 20100   lbf/in.²
Table 8.1 Natural Frequencies and Mode Shapes for Uniform Beams in Bending
End Conditions Characteristic equation and roots \beta_{i} L A B C D
Simply-supported (Pinned-pinned) \sin \beta_{i}L = 0 \\ \beta_{1}L = \pi \\ \beta_{2}L = 2\pi \\ \beta_{3}L = 3\pi \\ \beta_{4}L = 4\pi 1 0 0 0
Free-free \cos \beta_{i}L \cdot \cosh \beta_{i}L = 1 \\ \beta_{1}L = 4.73004 \\ \beta_{2}L = 7.85321 \\ \beta_{3}L = 10.99561 \\ \beta_{4}L = 14.13717 1 \frac{\sin \beta_{i}L – \sinh \beta_{i}L}{\cosh \beta_{i}L – \cos \beta_{i}L} 1 \frac{\sin \beta_{i}L – \sinh \beta_{i}L}{\cosh \beta_{i}L – \cos \beta_{i}L}
Fixed-fixed \cos \beta_{i}L \cdot \cosh \beta_{i}L = 1 \\ \beta_{1}L = 4.73004 \\ \beta_{2}L = 7.85321 \\ \beta_{3}L = 10.99561 \\ \beta_{4}L = 14.13717 -1 -\frac{\sinh \beta_{i}L – \sin \beta_{i}L}{\cos \beta_{i}L – \cosh \beta_{i}L} 1 \frac{\sinh \beta_{i}L – \sin \beta_{i}L}{\cos \beta_{i}L – \cosh \beta_{i}L}
Cantilever (Fixed-free) (x is measured from the fixed end) \cos \beta_{i}L \cdot \cosh \beta_{i}L = -1 \\ \beta_{1}L = 1.87510 \\ \beta_{2}L = 4.69409 \\ \beta_{3}L = 7.85475 \\ \beta_{4}L = 10.99554 1 -\frac{\sin \beta_{i}L + \sinh \beta_{i}L}{\cos \beta_{i}L + \cosh \beta_{i}L} -1 \frac{\sin \beta_{i}L + \sinh \beta_{i}L}{\cos \beta_{i}L + \cosh \beta_{i}L}
Fixed-pinned (x is measured from the fixed end) \tan \beta_{i}L – \tanh \beta_{i}L = 0 \\ \beta_{1}L = 3.92660 \\ \beta_{2}L = 7.06858 \\ \beta_{3}L = 10.21017 \\ \beta_{4}L = 13.35177 1 -\frac{\sin \beta_{i}L – \sinh \beta_{i}L}{\cos \beta_{i}L – \cosh \beta_{i}L} -1 \frac{\sin \beta_{i}L – \sinh \beta_{i}L}{\cos \beta_{i}L – \cosh \beta_{i}L}
8.1
2.3

Related Answered Questions