Discussion of “On the Evaluation of Dynamic Stresses in Pipelines Using Limited Vibration Measurements and FEA in the Frequency Domain” (Moussa, W. A., and Abdelhamid, A. N., 1999, ASME J. Pressure Vessel Technol., 121, Aug., pp. 241–245)

Truong, , K. T.

doi:10.1115/1.1359528

A practical technique in determining the dynamic stresses in pipelines using finite element method is investigated by Moussa et al.. The authors propose to create an FEA model to excite the piping structure by a single concentrated harmonic force at its midspan with frequencies ranging from 0 to 160 Hz. The dynamic vibrating displacements are then determined as well as the von Mises stress at all elements. The displacement responses are assumed as if they were measured and are then used as input to excite the model. The von Mises stress in both cases (force input and motion input) are compared together and plotted for different tested sections. For the case of a single displacement used as the loading source, the maximum differences are 0.0 and 9.4 percent for Sections 1 and 12, respectively (Fig. 2 of the paper).

Here, the same pipeline structure is restudied using an analytical approach as well as a user-oriented piping stress analysis software, CAESAR II. The results show that both approaches are nearly identical, but differ seriously from the aforementioned work. The piping model is described as follows:

Total length:

240 in.

Modulus of elasticity:

30E6 psi

Outside diameter:

4.5 in.

Mass density:

0.2825

lb / {in}^{3}

Wall thickness:

0.74 in.

No. of elements:

22

Boundary conditions: Clamped at point 1 and hinged at point 23 (Fig. 1)

Natural frequencies (first 3 modes): 11.67, 37.87, and 79.02 Hz

1

View large Download slide

Clamped-hinged piping model

Basic Equations

Consider an undamped single degree of freedom SDOF system that is subjected to a harmonic force

P (t)

with amplitude

P_{0}

and a circular frequency

ω ¯

(Fig. 2). The equation of motion is given by

My ¨ + ky = P_{0} sin (ω \bar{t})

(1)

The solution of Eq. (¹) is

y (t) = A cos ω t + B sin ω t + \frac{P_{0}}{k} \frac{1}{1 - r^{2}} sin ω \bar{t}

(2)

where r is defined as the ratio of the circular frequency of the externally applied load to the natural circular frequency of the system; that is,

r = \frac{ω ¯}{ω} = \frac{f ¯}{f}

(3)

2

View large Download slide

Undamped SDOF system subjected to a harmonic load

The solution given by Eq. (²) is the superposition of the free vibration problem and the effect of the exciting force exposed by the last term of Eq. (²), which involves only the frequency of the harmonic load. For frequency response or harmonic analysis, only the steady-state response is considered, and Eq. (²) becomes

y (t) = \frac{P_{0}}{k} \frac{1}{1 - r^{2}} sin ω \bar{t}

(4)

The solution for maximum displacement from an unphased harmonic analysis is then

δ_{dyn} = \frac{P_{0}}{k} \frac{1}{(1 - r^{2})}

(5)

Let

\frac{P_{0}}{k} = δ_{static} \Rightarrow δ_{dyn} = \frac{δ_{static}}{(1 - r^{2})}

(6)

where k is the structural stiffness of the SDOF system.

Force Excitation at Midspan

Analytical Calculation

The beam deflection equations for the case, as shown in Fig. 3, are well known and written as follows:

at load, Δ_{x} = \frac{7 {PL}^{3}}{768 EI}

(7)

when x < L / 2, Δ_{x} = \frac{Px}{96 EI} (3 L^{2} - 5 x^{2})

(8)

when x > L / 2, Δ_{x} = \frac{P}{96 EI} {(x - L)}^{2} (11 x - 2 L)

(9)

For the force load case of 1000 lb and the excitation frequency of 5 Hz, at midspan (Node 12), the static deflection is calculated as per Eq. (⁷) and is equal to 0.267 in. The dynamic deflection is

δ_{dyn} = \frac{0.267}{1 - {(\frac{5}{11.675})}^{2}} = 0.32 in.

At Node 8, 84 in. from the clamped end

(x = 156 > L / 2),

the same procedure could be applied to calculate the dynamic deflection using Eq. (⁹) and is equal to 0.231 in. The foregoing procedure is then used to obtain the dynamic deflection profile, as shown in Table 1. In order to determine dynamic stresses, one need only apply the following relationships for bending moment, shear, and bending stress:

M = - EI \frac{d^{2} y}{{dx}^{2}}; V = \frac{dM}{dx}; σ = \frac{M}{Z}

Table

Displacement results due to force load case at midspan $(F = 1000$ lb)

3

View large Download slide

Concentrated load at midspan

Computer Simulation

A harmonic dynamic analysis of the foregoing system has been performed using CAESAR II. The harmonic load of 1000 lb at a frequency equal to 5 Hz was applied at Node 12 (midspan). The computer output is tabulated in Table 1 for comparison purpose. As one can see, the displacement results are nearly identical in both approaches.

The next step is to take displacement response at a certain location to excite the system in order to simulate the displacement excitation case and evaluate the dynamic stresses.

Displacement Excitation

Computer Simulation

In order to compare with the results of the paper, the displacement at Node 8 (equal to 0.2308 in. in CAESAR II from the force load case) was used to excite the system with the harmonic frequency of 5 Hz. A computer run was performed and the output results are shown in Table 2.

Table

Displacement results due to imposed displacement at Node 8

Analytical Calculation

Regarding the analytical approach for this loading case, one needs to evaluate the equivalent dynamic force applied at Node 8 to create a displacement of 0.2308 in. at this location. With trial and error processing using CAESAR II, a final force equal to 1220 lb applied at Node 8 will generate a identical result as the displacement load case. The problem now becomes how to evaluate the dynamic displacement with a concentrated load at any point. The beam deflection equations for Fig. 4 are written as follows:

at load, Δ_{x} = \frac{{Pa}^{2} b^{3}}{12 {EIL}^{3}} (3 L + a)

(10)

when x < a Δ_{x} = \frac{{Pb}^{2} x}{12 {EIL}^{3}} (3 {aL}^{2} - 2 {Lx}^{2} - {ax}^{2})

(11)

when x > a Δ_{x} = \frac{Pa}{12 {EIL}^{3}} {(L - x)}^{2} (3 L^{2} x - a^{2} x - 2 a^{2} L)

(12)

At load point, i.e., Node 8, the static deflection is calculated using Eq. (¹⁰)

Δ_{x} = \frac{1220 {(156)}^{2} {(84)}^{3} (3 \times 240 + 156)}{12 (30 E 6) (16.047) {(240)}^{3}} = 0.193 in.

and the dynamic deflection is

δ_{dyn} = \frac{0.193}{1 - {(5 / 11.675)}^{2}} = 0.236 in.

The deflection at other nodes are also calculated using Eqs. (¹¹) and (¹²), and the results are also shown in Table 2.

4

View large Download slide

Concentrated load at any point

Results

As we can see, both approaches (analytical and computer calculations), offer very similar results with the maximum error not exceeding 4.5 percent. We are going to evaluate the stresses induced by the force load and displacement load at 5 Hz. We use the same equation proposed in the paper

γ = |\frac{σ_{F} - σ_{D}}{σ_{F}}| \times 100

where

σ_{F}

and

σ_{D}

are the bending stresses in the tested sections for the force and the displacement excitation case, respectively. The calculation is shown in Table 3. For comparison purposes, the γ values of the paper at 5 Hz are reported in Table 3; as we can see, the γ values differ seriously. At the fixed end (Section 1), it is not 0 percent, but 16 percent and at Section 12, 33 percent in lieu of 9.5 percent.

Table

Comparison of γ values at frequency of 5 Hz

Until now, we used only one value of frequency of 5 Hz to establish the confidence level between the analytical and computer procedures. The same calculation was then repeated for the following cases using the piping computer code:

Run 1. Force load case for frequencies ranging from 5 to 145 Hz with increment of 20 Hz.

Run 2. Use Run 1’s displacement at Node 8 as input for displacement excitation case for the same frequency range.

Run 3. Use Run 1’s displacement at Nodes 8 and 16 as dual translational vibration measurements (TVM) input for the same frequency range.

For finite element analysis point of view, Runs 2 and 3 are treated as an excitation due to multi-base motion. The natural frequencies of the piping system is no longer the same as in Run 1. They are as follows:

Frequency (Hz)

Mode

Run 2

Run 3

1

23.2

70.3

2

75.2

102.2

3

115.0

145.5

4

172.3

210.8

After stress runs completed, the γ values were then calculated with the attempt to reproduce the plots of Figs. 2 and 4(a) in the paper. The outcome is completely different; the results are shown in Tables 4 and 5, respectively, rather than plotted. This is because, during the load steps execution, the imposed displacement piping system will go through its natural frequencies and there will be resonance situations where the stresses become excessively large. In the stress output, some frequency related results have been eliminated or shifted because of the aforementioned resonance reason. In light of this noticeable discrepancy, the precision of stress calculation of the paper is seriously questionable and definitely not reliable, and only one conclusion could be drawn: the dynamic stresses are only identical if, and only if, the motion measurement is applied exactly at the same force position. When the loading positions are unknown or different between the force and displacement excitation case, all boundary conditions, mode shapes, as well as the stiffness matrix change accordingly, and the errors are unpredictable.

Table

Spectrum of γ using displacement at Node 8 as displacement excitation

Table

Spectrum of γ using displacement at Nodes 8 and 16 simultaneously as displacement excitation

by ASME

Discussion of “On the Evaluation of Dynamic Stresses in Pipelines Using Limited Vibration Measurements and FEA in the Frequency Domain” (Moussa, W. A., and Abdelhamid, A. N., 1999, ASME J. Pressure Vessel Technol., 121, Aug., pp. 241–245)

Basic Equations

Force Excitation at Midspan

Analytical Calculation

Computer Simulation

Displacement Excitation

Computer Simulation

Analytical Calculation

Results

Get Email Alerts

Cited By

ASME Journals

ASME Conference Proceedings

ASME eBooks

Resources

Opportunities

Discussion of “On the Evaluation of Dynamic Stresses in Pipelines Using Limited Vibration Measurements and FEA in the Frequency Domain” (Moussa, W. A., and Abdelhamid, A. N., 1999, ASME J. Pressure Vessel Technol., 121, Aug., pp. 241–245)

Basic Equations

Force Excitation at Midspan

Analytical Calculation

Computer Simulation

Displacement Excitation

Computer Simulation

Analytical Calculation

Results

Get Email Alerts

Cited By

Related Articles

Related Proceedings Papers

Related Chapters

ASME Journals

ASME Conference Proceedings

ASME eBooks

Resources

Opportunities

This Feature Is Available To Subscribers Only