Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Seismic Design Criteria of Refinery Piping Systems: June 2014

Download as pdf or txt
Download as pdf or txt
You are on page 1of 23

See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/263131870

Seismic Design Criteria of Refinery Piping Systems

Article · June 2014

CITATIONS READS

12 5,710

3 authors:

Fabrizio Paolacci Md Shahin Reza


Università Degli Studi Roma Tre Università degli Studi di Trento
166 PUBLICATIONS   1,149 CITATIONS    27 PUBLICATIONS   167 CITATIONS   

SEE PROFILE SEE PROFILE

O. S. Bursi
Università degli Studi di Trento
276 PUBLICATIONS   2,591 CITATIONS   

SEE PROFILE

Some of the authors of this publication are also working on these related projects:

SEISMIC PERFORMANCE OF MULTI-COMPONENT SYSTEMS IN SPECIAL RISK INDUSTRIAL FACILITIES (SPIF) View project

1st International Workshop on Risk and Resilience of Industrial installations Against Natural Threats and Mitigation Strategies View project

All content following this page was uploaded by Fabrizio Paolacci on 25 April 2016.

The user has requested enhancement of the downloaded file.


Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

SEISMIC DESIGN CRITERIA OF REFINERY PIPING SYSTEMS


Fabrizio Paolacci1, Md Shahin Reza2, Oreste S. Bursi2

1
University Roma Tre, Department of Engineering
Via Vito Volterra 62, 00146, Rome, Italy
paolacci@uniroma3.it
2
University of Trento, Department of Civil, Environmental and Mechanical Engineering
Via Mesiano 77, 38121, Trento, Italy
md.shahin.reza@unitn.it, oreste.bursi.@unitn.it

Keywords: Piping system, Seismic design, Dynamic analysis

Abstract

Piping systems are of paramount importance to many industries, e.g. refinery, oil & gas and
petrochemical plants, where they are often employed to transport raw and refined materials,
e.g. oil and gas, from one point to another connecting several components, such as tanks, dis-
tillations columns, furnaces and pressure vessels. Recent studies showed that such structures
are quite vulnerable under earthquakes, and reported many damages of piping systems and
their components which led to catastrophic accidents. Therefore, particular attention must be
paid to safeguard them against seismic events. However, there exists a clear lack of proper
seismic design guidelines for piping systems, whereas researchers have confirmed the inade-
quacy of available standards. As a result, in-depth research is required to understand better the
seismic design and analysis criteria of piping systems. Along this line, this paper studies the
main issues related to the analysis and design of refinery piping systems through a realistic
case study. Initially, seismic analysis and component design methods of refinery piping sys-
tems is analysed. A review of the current design approaches according to European (EN13480
- 3) and American (ASME B31.1 & B31.3) standards is illustrated by using a typical refinery
piping system as a case study. The analysis permitted to identify the limits of design standards
and some critical aspects of the problem, e.g. dynamic interaction between pipes and rack,
correct definition of response factors and strain versus stress approach. Finally, a discussion
on the main results of an experimental activity carried out on flange joints under strong cyclic
loading is presented, which allowed to highlight some limitations of current standards in
terms of conservativeness.
1 INTRODUCTION
The piping systems typically found in a refinery complex contain various components and
support structures and operate in a broad range of working environments. Some common
components usually used in piping systems include straight pipes, elbows, Tee-Joints, various
valves, flanged joints, pressure vessels, tanks, strainers and reducers. Depending on the nature
of the working fluids, piping systems are designed to work over a wide range of temperature
and pressure. A typical piping system is presented in Fig. 1.

Figure 1 Typical piping layout Figure 2 : Breakage of a piping flanged connection

Currently, both European and American codes are available for the design of piping sys-
tems in seismic-prone zones. The main European contribution is chiefly represented by the
standard EN13480, dedicated to metalling piping systems [1]. The Eurocode 8 - Part 4, the
European seismic code for industrial components, is also devoted to pipelines, but only of
above-ground type, which differs from metallic piping system for many aspects, and then use-
less.
American experience on piping system is instead very rich, especially in terms of design
standardization and seismic design calculation, and the long list of standards and codes avail-
able it is a clear demonstration of it. The main standard is represented by ASME B31.3 [2],
but many other contributions and guidelines are also available [3, 4].
The seismic analysis of a piping system involves several basic steps that allow defining the
proper seismic action, the suitable numerical model and analysis method and the verification
format to be used. The European (EN13480:3-2002) and the American standards (ASME
B31.3-2006) for piping systems differ for several of these aspects. Thus, in order to under-
stand these differences and the consequences on the seismic response evaluation, in the fol-
lowing both European and American standards are applied to a representative case study.
It is necessary to stress that the American Standard does not contain explicit indications on
the seismic analysis of piping systems, but rather refers to the American standard for seismic
analysis of structure ASCE7-05, which includes all the required prescriptions. On the contrary
the European Standard contains an entire Annex (A) dedicated to the dynamic and seismic
analysis of piping systems, but does not contains explicit quantification of the seismic action.
At this end the Eurocode 8 (prEN1998:1 2004) should be used.
Despite this articulated framework of standards and codes, that could help engineers to cor-
rectly design piping systems under seismic action, recent earthquakes showed a quite high
vulnerability of these structures, where damage ranges from the simple failure of joints to the
failure of supporting structures [5-8]. For example, the failure of a bolted flange connection is
shown in Fig. 2.
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

Consequences can be characterized by several degrees of severity, depending on the mate-


rial delivered by pipes. For dangerous liquids or gases, even a simple failure of a joint can
represent the trigger of a significant accidental chain, with severe consequences both for the
environment and human lives. Unfortunately, few contributions in the literature are available,
to clarify the seismic requirements that piping systems have to comply with. In addition, as-
pects like action and structural modelling have not yet been treated in a satisfactorily manner.
Moreover, current American and European standards do not contain enough rules and details
for a proper seismic analysis and design of piping systems.
Along these lines, some problems relevant to seismic analysis criteria of piping systems are
addressed in this paper. In a greater detail, several aspects that characterise the problem are
treated: 1) modelling of pipes and pipe-racks; 2) selection of the analysis method; 3) defini-
tion of the seismic action; 4) dynamic analysis of the system; 5) stress analysis of pipes; 6)
definition of the ultimate capacity of pipes and joints between pipes.
To this end a representative case study of an actual piping system is analysed. The re-
quirements of American and European codes were compared and important aspects were
highlighted, like: a) dynamic coupling between pipes and pipe-rack, often erroneously ne-
glected; b) definition of proper restraint conditions between pipes and support structures and
between adjacent piping systems; c) evaluation of the ultimate capacity of pipes and joints
necessary for a correct design of a structure, as suggested by modern approaches like Perfor-
mance-based Engineering.
About the aforementioned item c) an experimental campaign is has been undertaken at the
University of Trento in order to characterise the cyclic behaviour of flanged joints between
pipes, particularly important to avoid leakage of dangerous substances during a severe seismic
event. Some information about the testing results and conclusions are herein illustrated and
commented, in the light of the current codes for piping systems.

2 DESCRIPTION OF A CASE STUDY


The piping system here analysed belongs to a refinery, whose plan view is shown in Fig.3.
The supporting steel structure is composed by seven transverse moment resisting frames
placed every 6 m, made of commercial HEA/B steel profiles. In the longitudinal direction it
behaves as truss structure, which is reinforced with 6 braces (see Fig. 4). Horizontal bracings
are also installed to avoid excessive relative displacements between the pipe supports.
The piping system presents a typical piping layout with pipes having different diameters.
To simplify the analysis, only the structural contribution of 8’’ pipes has been considered,
whose layout is shown in Fig. 5. The remaining pipes are considered only as weight. Several
flange elbows are present within the piping system.
z  

x  
y  
(1)   (2)   (3)   (4)   (5)   (6)   (7)  

Figure 3 – Plan view of the refinery Figure 4 – Case Study

The fluids contained in the pipes are several, but essentially Amine, cooling water and high
to medium pressure steam. The vertical loads corresponding to the weight of the pipes, insula-
tion and fluid are considered as uniformly distributed equal to 12 kN/m. The main characteris-
tics of the piping system are the following:

• Structural steel S-275 JR according to EN 10025 (2005)


• Pipe steel - ASTM A106 Grade B
• pipes with diameter of 8’’
• Pressure of the pipes: 0.5÷5 MPa
• Temperature range 47 °C ÷360 °C.
• Seismic category of Importance Ip=1.5, PGA=0.24 g, Soil conditions: D

First floor Second floor Third floor


Figure 5 - Pipes Layout
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

3 MAIN ISSUES IN SEISMIC DESIGN OF PIPING SYSTEMS


The seismic design of a piping system entails many issues. As stated above, they are essen-
tially related to the modelling of the structure, to a correct definition of the seismic action, to a
proper analysis method to be applied, and finally, to a appropriate design method. In addition,
nowadays, the Load and Resistance Factor Design method (LRFD) is certainly the standard
method for designing structures and the allowable stress design method is by now abandoned
because often considered too much conservative. Unfortunately, this latter is still the current
approach for designing a piping system. This represents a limit to be overcome. But to the
purpose, it is firstly necessary to understand the limit of such an approach; and the more direct
way is to apply the codes prescriptions to the case study, trying to individuate limits and
drawbacks. In doing this, it seems proper to compare European and American experiences,
applying both the EN13480:3 and ASME B31.3 to the case study described above along with
a comparison of the results regarding all the treated aspects.

3.1 Definition of the numerical model


A synthetic scheme of what European and American standard prescribe for a correct nu-
merical modelling of a piping system is reported in Table 1, according to the seismic condi-
tions defined in chapter 3.3.
The table clearly shows that the suggested numerical model for the seismic analysis is al-
ways elastic both for EN13480 and ASME B31.3. This choice comes certainly from the old
way to evaluate the safety level of a structure: the allowable stress method, still diffused in
designing of piping systems. Usually, only the piping system is modelled, using the support-
ing structure to evaluate the seismic action at pipes level (e.g., in-structure spectra). The sup-
porting structure (e.g. pipe-tack) is treated as elastic too. The assumption of elastic behavior
would not be a strong limitation if a correct value of the behavior factor were adopted. Some
comments on this aspect will be provided afterwards.
Another important aspect, often related to the assumption on the numerical model is the
analysis method to be adopted. This ranges from the very simple equivalent static method to
the time-history analysis. This aspect will be treated in detail in the next section.
A key point in modelling a piping system is the possibility to neglect the interaction (static
and dynamic) between the pipes and the supporting structure. EN13480 does not provide any
indication, whereas ASME B31.3, by means of ASCE-07, prescribes a crude rule based on the
ratio WR between the weights of pipes and supporting structure. In particular, if WR <
25% the interaction can be excluded and the piping system can be considered as a non-
building structure, loaded by a seismic action coming from the supporting structure at pipes
level.
This rule has been recently analysed by several authors. For example in [9] the rule has
been studied using time-history analysis. From the results and discussion the author conclud-
ed that in some cases this decoupling rule could produce gross errors in the evaluation of the
dynamic behavior of piping systems. In particular, it seems that in dynamic assessment of
such systems, in addition to the primary-secondary system weight ratio criteria, attention
should be paid to other aspects as “end conditions of pipes”, “relative stiffness of supporting
structure to piping system” and “relative stiffness of pipes to pipe-supports”, even though on-
ly partial conclusions where reached by the authors, that suggested more investigations on this
matter. In the literature many works have been dedicated to the problem of dynamic coupling
between primary and secondary systems [10-12], but the results are difficult to be extended to
the case of piping systems for several reasons, but mainly because in case of piping systems
the secondary system is composed by more sub-systems (pipes with different diameters, dif-
ferent end conditions, different supports, etc..). Therefore, assumption of a single secondary
system could lead to gross errors in response prediction of piping systems.

Table 1 – Code prescriptions for numerical modelling as well as seismic conditions

Code Model Analysis Seismic Dynamic Relative Pipe modelling


type condition interaction motion
EN13480:3 Elastic (1) Equivalent Static SSE (2) NO YES Beam elements
Modal (G,IS4) OBE (FF and SIF)
Time-History
ASME Elastic Static SSE YES/NO YES Beam elements
B31.3 Modal (FF and SIF)(3)
Time-History
(1) elastic calculation shall be used although some part might be exhibit plastic deformations (p. 4.1), (2) SSE=safe shutdown earth-
quake, DBE= Occasional operation condition, (3) FF=Flexibility Factor, SIF= Stress Intensification Factor, (4) G=ground motion spectra,
IS=In-structure spectra.

Concerning the case study, applying the weight ratio rule, dynamic interaction has been
considered, because WR > 25%. Moreover, for comparison, this rule has also been considered
in applying the European standard. In order to evaluate the effectiveness of the weight rule,
the behavior of the piping system without pipe-rack has also been analysed.
The pipes are usually connected to the pipe-rack by mechanical supports, often flexible.
They are usually designed to accommodate thermal and pressure effects, avoiding excessive
stress in the pipes. Unfortunately, no indications on how to model them are provided by Eu-
ropean and American standards. Nevertheless, this is an important aspect that have to be
treated with particular attention because can cause important changes in dynamic behavior of
the system. Dissipative supports can be used instead of traditional ones in order to introduce
artificially more damping and reduce the response of the pipes [13].
The analysed case present a quite stiff support systems, modelled as elastic spring in the
transverse direction (Y), leaving free the relative displacements in longitudinal direction (X)
and using fix restraints conditions in vertical direction. Moreover, as usual, all the rotations
between pipe and pipe-rack have been unrestrained.
Another relevant aspect about modelling of piping systems is the adoption of a proper
model for pipes and fittings (elbows, tee-joints, nozzles, etc.). At this regard, usually beam
elements with hollow section are used for straight pipes. The fittings are also modelled using
beam elements, but modifying the stiffness for the effect of geometry. At this regards both
European and American codes define a flexibility factor (k >1) using which the moment of
inertia of the pipe is reduced. In addition, to take into account the stress concentration effect,
the Stress Intensification Factor (SIF) is used to increase the stress calculated using the beam
theory. The value of k and SIF calculated according to EN13480 and ASME B31.3 are very
similar (see Table 2).
In presence of high pressure condition the stiffness of the pipe could increase. To account
for this effect the flexibility factor k is reduced. Unfortunately, only ASME B313.3 provides
an explicit expression of the pressure reducing factor. For the analyzed case and for a pressure
of 5 MPa this factor halves the flexibility (see Table 3).
Alternatively, it is possible to use shell elements to model fittings [14]. This approach is
appropriate to account for ovalization of the section and stiffening pressure effect. For this
reasons this model has been used for the case study and a comparison with beam model has
been carried out (Fig. 6). We can anticipate that the numerical simulations have shown a simi-
lar behavior of both the models and the reliability of modified beam element, at least for
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

standard fittings, like the one here analyzed. Therefore, in what follows only the results of the
more refined one will be shown.

Table 2 Flexibility factor and Stress Intensification Factor for unflanged elbows

Standard k SIF h
In plane Out-of-plane

EN13480:3
-----

ASME B31.3
(*)

(*) ASME B31.3 allows to use in-plane SIF also for out-of-plane bending.

Table 3 – Pressure factor for k and SIF

code SF SSF
EN13480 ---- ----

ASME B31.3

(*) P= pressure, E=elastic modulus of steel

A last but not less important aspect regards the boundary conditions of the pipes. In fact,
because a piping system is realized by hundreds of miles of pipes, the analysis involves nec-
essary a limited part of the structure. Consequently proper boundary conditions have to be ac-
curately adopted to simulate the remaining part of the structure. Also for this delicate aspect
no indications are provided by the European and American standards. As already shown in
literature, the correct boundary conditions may cause important modifications of the dynamics
of the piping system but the correct choice depends on the single case [9].
Figure 6 – Shell FEM for an elbow

Because one of the aims of this work is to compare European and American standards, it
has been decided to limit the possible cases adopting as restraint conditions of the pipe ends
only hinges. In this respect, one can look at Fig. 7.

Trasversal  pipe    
Simulated  by  restrain  
conditions  

 
Fig. 7 – Boundary condition in one if the pipe ends

The model of the piping system and the support structure has been built in the general pur-
pose software MIDAS Gen [15]. It is shown in Fig. 8.

 
1   2  
z   3   4  
x   y   5   6   7  
B  

A  
   
Figure 8 – The FE model of the piping system

3.2 Seismic actions and analysis methods


Both European and American standards assume the following two types of analysis man-
datory for the pipes:

• Movements due to inertia effects


• Differential movement of the supports (within the supporting structure or between adja-
cent pipe-racks)

The first type of analysis is essentially related to the effects of the absolute acceleration on
the pipe mass. The second one is due to the relative movements between two supports, within
the supporting structure or belonging to adjacent structures. Often the relevant effects are due
to the displacement effect rather than acceleration effects.
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

Concerning the case study, the entire model here considered (pipe + pipe-rack) allow iden-
tifying both the effects. On the contrary the model without considering the supporting struc-
ture, here also considered, allows identifying only the inertia effect of the pipes, unless a
multi-support excitation would be used.
The seismic action for pipe-racks is usually represented by design response spectra or ac-
celerograms (natural records or synthetic accelerograms). For the analysis of pipes only, “In-
structure” spectra or “filtered response spectra” are used.

Design response spectrum method


The design spectra are the main representation of a seismic action and usually are defined
by the seismic codes in terms of hazard conditions of the site, the level of dissipation capabil-
ity of the supporting structure and pipes (response or behavior factor), the right level of damp-
ing to be employed, and the level of structure reliability to impose, identified by the
importance factor.
Concerning the support structure, hazard conditions apart, the damping usually adopted is
equal to 5%, as suggested by Eurocode 8 and 3 for steel structures, whereas the behavior fac-
tor q depends on the type of structure used for the pipe-rack. While for building-type structure
this aspect has been well identified and quantified, for structures like pipe-racks, that may of-
ten be considered as non-building structure (ASCE-07:2005), the problem may be quite dif-
ferent.
The current American and European seismic codes provide a q factor for steel racks equal to 3
½ and 4 respectively (Table 4). This choice probably derives from the hypothesis of no-
coupling between the rack (primary system) and the pipes (secondary system). In fact, usually
the level of dynamic coupling between pipes and rack can be neglected. But in other some
cases its influence cannot be excluded a priori [9].
In order to compare the results of EN13480 and ASME B31.3 the same elastic spectrum, has
been here adopted defined according to prEN1998:1, and modified using the behavior factor
evaluated using Incremental Dynamic Analyses [16].
The resultant forces and displacements from bi-directional analysis are typically obtained
by the square root sum of square of the response in each direction, or by applying the so-
called 100-30 rule.

Figure 9 First two vibration modes of the pipe-rack model


In Table 5, the main dynamic characteristics of the system are shown. As mentioned before,
the modal analysis on the entire system allows to highlight the important role of the pipes in
realizing the structural coupling between the several frames of the pipe-rack. For example in
Fig. 9 and 10 the vibration modes of the rack with and without the pipes is shown. They are
quite different especially in terms of excited mass. For example, the period of the first mode
of the rack with and without pipes is very similar, whereas the excited mass is higher in the
first case, showing the coupling effect of the transverse frames due to the pipes.

Figure 10 First three vibration modes of the pipe plus pipe-rack model

Table 4 – Behaviour factor, damping and Importance Factor

Code Behavior factor response factor Rp Damping Importance F.


R (pipe-rack) (pipe) (%)
EN13480:3 (EC8) 4 Not indicated 5 1.5
ASME B31.3 3¼ 6÷12 5 1.5

Table 5 – First three vibration modes of the entire piping system

Case Mode 1 Mode 2 Mode 3


f (Hz) T(sec) MPS(%) f (Hz) T(sec) MPS(%) f (Hz) T(sec) MPS(%)

TD-PR (*) 2.17 0.459 28.29 2.33 0.456 47.62 2.93 0.440 5.68
TD- 2.15 0.456 67.61 2.33 0.429 8.96 2.93 0.353 7.64
PRP(*)
LD- 2.81 0.355 53.27 5.15 0.194 36.35 --- --- ---
PRP(*)
(*) TD-PR=Transversal direction – Pipe-rack, TD-PPR=Trans. dir. – Pipe-rack+pipes, TD-PR=Longitudinal dir. – Pipe-rack+pipes

In-structure spectra
The in-structure spectra allow to define the seismic action for the single pipe at several
floors of the pipe-rack. Both European and American codes give their explicit expressions.
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

They are defined as the spectrum acceleration multiplied by the amplification factor AF
shown in Fig. 11. The expression of the amplification factor expressed by EN13480 (EC8)
and ASME B31.3 (ASCE-07) differ for the presence of the period ratio Ta/T1, presents only
in the seismic European Standard which takes into account the dynamic interaction between
pipes (Ta) and pipe-rack (T1); see Table 6.

Amplification Factor AF
6
5.5 Ta/T1=0.01
Ta/T1=0.1
Ta/T1=0.222
5 Ta/T1=0.5
AF_E ( z , 0.01 ) Ta/T1=1.0
AF_E ( z , 0.1 ) Ta/T1=2.0
AF_E ( z , 0.222) 4 ASCE 07
AF_E ( z , 0.5 )
AF_E ( z , 1 )
AF_E ( z , 2 ) 3
AF_A( z)

1
1
0 0.2 0.4 0.6 0.8
0 z 1

Figure 61 Amplification factor of in-structure spectra v/s Ta/T1

In case of pipes, because in general any yielding phenomenon in the pipes is avoided, the
dissipation capability is generally restricted only to the supporting structure and the relative
response modification factor depends on its structural configuration. In some cases dissipation
phenomena may have place also in the pipes and a specific behavior factor have to be defined
if the elastic analysis with response spectra is used.
The behavior factor provided by the codes, especially by the American one, seems to be
overestimated. For example, ASCE-07 prescribes the use of a behavior factor 6 or 12 accord-
ing to the deformability of the material used. In some cases this hypothesis may not be totally
true.
For example Okeil and Tung in 1996 [17] using an idealized piping system have suggested
a closed-form formula that provides the reduction factor q as function of the ductility of the
support structure, µ and the piping frequency to seismic action frequency ratio. The reduction
q increases with the ductility, µ. Moreover, the stiffer is the piping system (low values of f)
the higher is q. In any case q is never greater than 3-4. These results are in contrast with the
provisions of ASCE-07.
Other authors reached similar conclusions as well. For example in [18] the authors investi-
gated on the dynamic response of a piping system on a rack with gap and friction. By using
dynamic harmonic analysis on a simple system they found that in presence of nonlinearity
(friction between pipes and rack) the reduction of the acceleration, with respect to the elastic
case, could be of the order of 2-3, value suggested also by the American Lifelines Alliance [3]
and FEMA 450.
Instead of using static method employing the seismic force of Table 6, a dynamic approach is
also possible, using generated floor response spectra, whose shape is based on typical Soil-
Structure interaction theory. For industrial piping systems few contributions have been found
in literature [6-8, 19].
For the case study of Fig. 4, an example of floor spectra generated by using time-history
analyses by the 1940 N-S El-Centro record is shown in Fig. 12, where is clear the filtering
effect of the pipe-rack. The maximum amplification is now restricted to the range 2-3 Hz. The
labels Frame 1, 3 and 5 correspond to the labels of the frames indicated in Fig. 6. It is also
possible to note the variability of the maximum amplification effect within the structure. The
maximum acceleration peak is found at frame 3; this because more mass is placed there with
respect to the other frames.

Table 6 In-Structure spectra definitions

Code Expression note


ASME B31.3 ap=2.5 (Tab. 13.6-1)
0.4ap SDS Wp
(ASCE-07) Fp =Ip AF SDS=Maximum spectral acceleration
Rp
z Wp=Weight of the pipe
AF= 1+2 z=height where the pipe is placed
H
H=total height of the pipe-rack
Ip=Importance factor
Rp=response factor
EN13480:3 (EC8)
ag SWp
Fa =Ip  AF ag=Peak Ground Acceleration
Rp
S=soil factor
z
3 1+
H -­‐ 1 Ta=fundamental period of the pipes
AF=
T 2 2 T1=fundamental period of the pipe-rack
1+ 1-­‐ Ta
1

Because the natural frequencies of the piping systems are in the range 15-40 Hz the ampli-
fication effect due to the inertia is very limited. In fact, the maximum acceleration applied to
the pipes corresponds more or less to the maximum acceleration at floor in which the pipes
are placed, some peaks a part, at frequency 7 and 10 Hz.
 
8 5

7 4.5

6
4 Frame  5
Frame  5
Acceleration  (g)

3.5 Frame  3
Acceleration  (g)

5 Frame  3
3 Frame  1
4 Frame  1 2.5
Envelope
3 2

2 1.5

1 1

0.5
0
0
0 5 10 15 20 0 5 10 15 20
Frequency  (Hz)
  Frequency  (Hz)
 
Figure 72 - Case study - Spatial variability of the floor Figure 83 - Case study – Mean floor spectra for sev-
spectra for the El-Centro earthquake eral frames and Envelope Floor Spectrum

 
A possible representation of floor spectra for the case-study is shown in Fig.13 as the enve-
lope of the mean spectra at several frames. This allows to account for the spatial variability of
the action and its frequency content. This has been used to get the response of the pipes ana-
lyzed individually.

Time history analysis


A time history seismic analysis is rarely used for the seismic design or retrofit of piping
systems. Often it is used to generate facility specific response spectra, or as a research tool, to
study in detail the behavior of a component or system.
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

Nowadays, the scientific community has widely accepted the use of natural records to re-
produce a real input, for several reasons. For many engineering application, the purpose of
selection and scaling of real earthquake is to fit the code design spectrum considering the
seismological and geological parameters of the specific site. To comply with the seismic
codes set of accelerograms, regardless its type, should basically match the following criteria:

• minimum of 3 accelerograms should be used;


• the mean of the zero period spectral response acceleration values (calculated from the indi-
vidual time histories) should not be smaller than the value of ag× S for the site in question
(S is the soil factor, ag is the Peak Ground Acceleration);
• in the range of periods between 0.2 T1 and 2 T1, where T1 is the fundamental period of the
structure in the direction where the accelerogram will be applied, no value of the mean 5%
damping elastic spectrum, calculated from all time histories, should be <90% of the corre-
sponding value of the 5% damping elastic response spectrum.

To help engineers in selecting a proper set of records, some tools have been proposed in
the literature. The most recent is REXEL proposed by Iervolino et al [20].
More rarely artificial or synthetic accelerograms are used. The first ones are generated to
match a target response spectrum, the second ones from seismological source models and ac-
counting for path and site effects. For artificial signals, even though, it is possible to obtain
acceleration time series that are almost completely compatible with the elastic design spec-
trum, the generated accelerograms often have an excessive number of cycles of strong motion,
and consequently have unrealistically high-energy content. For this reason are less used than
natural records. Moreover, to generate synthetic accelerograms there is the need to define a
specific earthquake scenario in terms of magnitude, rupture mechanism, in addition to geolog-
ical conditions and location of the site. Generally, most of this information is not often availa-
ble, particularly when using seismic design codes. This representation of a seismic input is
rarely used.
According to the above considerations, for the analysed case study, it has been decided to
use natural records and select 7 accelerograms compatible with the spectra (EC8, type 1, soil
D) of Fig. 14 (see Table 7 and 8). For the case study the vertical component has been consid-
ered as negligible, so only bi-directional motion has been adopted. The records were taken
from Pacific Earthquake Engineering Research center database (http://peer.berkeley.edu), us-
ing the following hazard and compatibility conditions (Table 7). Fig. 14 shows the elastic
spectrum of the selected natural records, in comparison with the target spectrum.

Table 7 – Hazard and compatibility conditions for natural record selection

Magnitude Soil conditions Distance PGA Lower bound Upper bound


Mw (Km) (g)
6÷7 D 0÷20 0.35 -10% +30%
Fig. 14. Elastic spectra of the set of natural accelerograms;

Table 8 - List of the selected natural records

Epicenter PGA_X PGA_Y Site


Earthquake Name Date Mw
[km] [m/s2] [m/s2] Class
Montenegro 24/05/79 6.2 20 0.56 0.54 B
Erzincan 17/12/39 6.6 13 3.81 5.03 B
Ano Liosia 07/09/99 6 18 1.09 0.84 B
Montenegro 24/05/79 6.2 20 0.56 0.54 B
Campano-Lucan 23/11/80 6.9 16 1.53 1.72 B
Umbria Marche 26/08/97 6 11 5.14 4.54 B
Montenegro 24/05/79 6.9 16 3.68 3.56 B

Differential movement effects


No specific indications are provided by the codes for the differential movement between
adjacent supports or structures. Only ASCE-07 provides a simplified criterion based on the
elastic analysis of the pipe-rack. Chapter 13.3.2.1 suggests evaluating the relative displace-
ments between two connection points within the structure and at the same level using the dif-
ferential movement for each vibration modes combined using proper modal combination rules.
Unlike the cases of differential movements between adjacent pipe-racks connected by
pipes, this effect on pipes within the structure is usually neglected. As observed above, the
complete model, pipes + pipe-rack, allows to automatically accounting for this effect.

3.3 Piping stress analysis and checks


One of the fundamental steps for the qualification of a pipe system is the fulfillment of
some limits of the pipe stress or strain, for a given working condition. For a seismic action,
usually two working conditions are considered:

• Operating condition or design basis earthquake condition: (OBE) is that earthquake


which, considering the regional, local geology and seismology, could be reasonably
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

expected to affect the site during the operating life of the plant. It is the earthquake
during which the operating conditions of the plant can be still assured.
• Safe shutdown earthquake condition: (SSE) is the maximum ground motion for which
some critical components of the plant must be designed to remain functional.

A rational definition of the OBE and SSE seismic conditions, according to the current codes,
is provided in [21].

In order to evaluate the safety level stress-based or strain-based approach can be used. The
first approach intends to evaluate the maximum stress in the pipes and the calculation is usu-
ally based on elastic analysis of the structure.
While stress based approach for pipelines is acceptable for a material with a well defined
yield point and with a well defined yield ductility and strength, this design criteria becomes
invalid when the stress in pipelines exceeds the limit under some displacement control loads,
such as earthquakes and landslides [22]. In this case, strain based approach provides the de-
sign rule where the strain in the pipeline is allowed to exceed the specified yield strain pro-
vided that the safe operation can be ensured under displacement load. This method allows
selected extensions to the stress-based design possibilities to take advantage of steel’s well-
known ability to deform plastically, but remain a stable structure. Codes and Standards are
available for the strain based design approach, see [22] for reference. With the strain-based
approach maximum strain in the pipes is calculated and compared with specified strain-limits
related to limit states usually identified with buckling or ovalization of the pipes. Unfortunate-
ly, this approach needs of the calculation of the seismic response in the non-linear range. This
is one of the reasons because the stress approach in more used.
Only EN13480 contains explicit indications for calculating the pipe stresses limits, consid-
ering OBE and SSE conditions, whereas ASME B31.3 indicates only occasional load condi-
tions that can be identified with the OBE condition.
For the verification of the pipes against earthquake, the allowable stress approach is usual-
ly adopted and will be also adopted in the following.

Load combinations and stress calculation


The response to seismic and other loads (sustained, thermal, pressure, etc..), have to be
properly combined. European and American code prescribe similar combinations. In this re-
spect the seismic load prescribed by the code (prEN1998:1 and ASCE-07) can be considered
as an exceptional seismic action. Under this condition, usually Load and Resistance Factor
design (LRFD) approach is adopted. If the allowable stress approach is used, the seismic ac-
tion has to be reduced, as usual, of a certain safety factor, typically 1.4 (see ASCE-07).
Because the allowable stress approach is widely diffused in piping design in the following
section it will be used for the calculation of the stresses and the verification of the safety level.
In doing this, the ASCE-07 combination formula will be adopted.
ASME B31.3 does not provide an explicit equation for calculating the longitudinal stress,
whereas EN13480 provides at point 12.3.3 the formula to evaluate the longitudinal stresses
due to sustained, occasional and exceptional loads (e.g. the earthquake). This comes from the
beam theory and includes internal pressure P and the resultant moments due to the sustained
loads (MA) and the earthquake (MB):

pD M + MB
σ= + 0.75 × SFI A
4t Z (1)
where SIF  is the stress intensification factor.
Because ASME B31.1 provides a similar expression of the stress, in the following this will be
also adopted for ASME B31.3.
ASME B31.3 at point 319.4.4 indicates the way to calculate the resultant moments, similar to
that of EN13480:3. In particular:

MR = (SIF y ) (
× M y 2 + SIFz × M z )
2

(2)

This definition has also been adopted for calculating the resultant moments using EN13480.
Studies have shown that the present standards for piping system design under seismic loads
are over conservative and modifications have been proposed to relax this over conservatism
[23, 24]. For example, [24] has applied the present codes to calculate the stress limit on piping
systems and experimental results showed significant discrepancy from the real behavior. This
aspect will be treated in further studies.

Definition of the allowable stress


EN13480 includes par. 12.3.3., which is dedicated to the definition of the allowable stress
calculation. They are defined according to the conditions previously recalled: Operating Basic
Earthquake (OBE) and Safe Shutdown Earthquake (SSE). The basic inequality to be respected
is:

σ ≤ kf h (3)

where

σ  = maximum stress in the pipes due to the sustained and seismic loads.  
k  = 1.2 for design basic earthquake  
k  = 1.8 for safe shutdown earthquake
fh  = is the basic allowable stress given by code

ASME B31.3 does not provide explicit differentiation between OBE and SSE but refers to a
Design Earthquake that can be identified with the OBE condition [21]. It is prescribed that the
maximum stress cannot be greater than 1.33 times the basic allowable stress for pipes, indi-
cated in the Appendix A of the same code.

3.4 Calculation of the seismic response of the case study


The results for each analysis method applied to the Case Study presented in previous sec-
tions are reported in Table 9 in terms of moments along the local axes y and z of the pipe. The
resultant moment, MR, of the single moments along local axes y and z, calculated according
to EN13480 and ASME B31.3 are also shown.
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

Table 9 – Results relevant to the seismic analysis of the case study

The maximum moment is found near the left edge of the rack (bay 2), even if similar val-
ues are also obtained within bay 6 and 7. In addition, the maximum stress level of the pipe in
the same points has been calculated according to Eq. (1), obtaining a mean value in case of
non-linear time history analysis of about 86 MPa against the allowable stress indicated in Ta-
ble 10.
The above calculations confirm that, using allowable stress approach and then comparing
the capacity with the demand (OBE), the system is highly overdesigned. The attempt to modi-
fy Eq. (1) is not enough to fully exploit the plastic properties of pipes. In this sense the strain-
based approach is more promising, especially if the Performance-Based approach would be
used to design new piping systems or to assess the seismic behavior of existing ones [25].
This is the main reason that induced researchers to explore the plastic behavior of pipes
and pipe-joints using the way of experimental investigation. More details can be found for
example in [21] where the results of an experimental campaign on non-standard bolted pipe
flange joints (see Fig. 15 and 16) are presented and discussed. In particular, from the above
Case Study, it was found that the leakage loads were well above the allowable design loads
and loads demanded by the earthquake.

Table 10 – Pipes allowable stress (MPa) for the analysed case study

Code OBE SSE


EN13480 165.60 248.40
ASME B31.3 183.50 ---

See, for example Fig. 16, where the moment-rotation diagram of a non-standard bolted
flange joint for an 8” pipe, de-signed according to the Eurocode EN1993-1-8 [26] is shown.
One may clearly note a safe margin between leakage moment and allowable moment.
Fig. 15 Experimental Setup for cyclic tests on elbows Fig. 16 Moment-rotation cycle of a flange joint after [21]
[21],

Fig. 17 shows the maximum relative displacement calculated according to chapter 3.2,
which include or not the presence of pipes in the support structure. It can be noted that the
presence of pipes reduces the displacements, especially in the central part of the pipe-way,
whereas it is influenced by the boundary conditions in case the presence of the pipe would be
considered.

Table 11 Maximum relative displacements (m)

bay 1 2 3 4 5 6 7
SS RSM 0.030 0.070 0.089 0.086 0.072 0.054 0.039
pipes T-H 0.037 0.053 0.003 0.024 0.026 0.026 0.020
RSM 0.034 0.078 0.086 0.082 0.002 0.074 0.130
SS
T-H 0.041 0.043 0.003 0.021 0.021 0.026 0.050
SS= support structure, SS pipe= support structure + pipes

In Table 11, the maximum relative displacements for both the DA-RSA and DA-NTHA
methods are also reported. The influence of pipes in the relative displacement of the pipes is
evident. Moreover, the displacements calculated using DA-NTHA analysis method are always
the lesser, especially in the central part of the pipe (bay 2-6).
This corresponds to a certain degree of conservatism inherent in the combination of the
several vibration modal responses using SSRS rule.

0.1
0.08 w  pipes
0.06 w/o  pipes
0.04
0.02
0

2 3 4 5
6
7
Fig. 17 Maximum relative displacement for response spectrum method
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

4 CONCLUSIONS
The seismic analysis of piping systems is very different from the analysis of like-building
structures for which a long experience allowed to have enough provisions for obtaining a reli-
able design against earthquakes. Therefore, the present contribution tried to clarify all the de-
sign steps of this type of structures and to identify all aspects that need to be clarified. This
was made through a representative case study analysed according to both European
(EN13480:3) and American standards (ASME B31.3) devoted to piping systems. The com-
parison between EN13480 and ASME B31.3 standards showed that the two codes provide
similar indications for the evaluation of the seismic response of pipelines supported by a pipe-
rack. From the analysis of the results the following conclusions, considered valid for both the
standards, can be drawn.
ASCE-07 provides a simple rule to establish when the dynamic coupling between pipes
and supporting structure need to be considered. It is based only on a weight ratio. Actually,
other ingredients should be included; for example the vibration periods range of the pipe. In
the analysed case the single pipe has a limited weight and therefore short vibration periods.
This means that, as already shown, the dynamic coupling between pipes and the pipe-rack is
limited; whereas the relative displacements between the pipe supports would represent a more
relevant effect. The pipes provide a sort of static coupling between the transverse frames. In
fact, the models with and without the pipe stiffness contribution have more or less the same
period, whereas the excited mass is differently distributed between the several vibration
modes. Therefore, the consequence of considering the structural contribution of the pipes is
the introduction of a sort of a rigid floor effect that couples the horizontal movements of the
nodes of the same floor. For the above-mentioned reasons an extensive analysis on this par-
ticular but important aspect is recommended.
The in-structure spectra suggested by prEN1998:1 and ASCE-07 differs for a term that de-
pends on the pipe-structure frequency ratio, included only in the European code. This aspect
is strictly related to the previous observation on the dynamic coupling. The exclusion of the
dynamic interaction could introduce a high error in the evaluation of the amplification factor
of in-structure spectra. Because usually the frequency of the piping system is relatively high
with respect to that of the supporting structure, the amplification effect would be relatively
low. In any case, to cover all the possible cases, the use of the European formula for in-
structure spectra is highly recommended. The result of this investigation suggested at point 1,
should also help to clarify this important aspect toward the direction of simplified but more
realistic formulas for in-structure action on pipes.
In the analysis both beam element and shell elements were used to model elbows. The
comparison showed the reliability of the more simple elements based on the beam theory, at
least for the simple case of elbows. For more complicated situations, the use of finite shell
elements should be considered.
The behavior factor q is usually indicated by standards. In particular the values provided
for piping systems by ASCE-07 seems to be too high, i.e., q=6-12. This is in contrast with the
idea of avoiding yielding phenomena in the pipes and to maintain operational conditions also
in case of strong seismic events. Certainly, the contribution of the supports in providing dissi-
pation capability should be taken into account; however on the basis of several investigations
found in the literature, this exclude the possibility of having so high values of q. For this rea-
son the next step will be an investigation on the reduction factor q also considering the non-
linear behaviour of pipes. This could be performed using non-linear beam elements for the
straight pipes, and plastic hinges for the simulation of the elbows in the plastic range. In doing
so, time-history analysis would be used, even if static non-linear analysis is nowadays very
diffused, but difficult to be applied to piping systems for the presence of a uniformly distrib-
uted mass of pipes.
The usual way of designing piping systems is based on the allowable stress approach. This
means that the structure is considered elastic. Consequently only an operating basis earth-
quake condition (OBE) can be taken into account for design. Conversely, the modern ap-
proach to the seismic design of structures is to differentiate serviceability from ultimate limit
states. This latter condition could be represented, for example, by the Safe Shutdown Earth-
quake (SSE). The problem is that a proper definition of the deformation capability of the
pipes and fittings (e.g., ovalization) in the plastic range would be necessary, and on this aspect
the research has not reached solid conclusions. But the advantages would be several and first
of all the design optimisation. For example in the analysed case, we are very far from the
moment capacity of the bolted flanges located in the pipeline. This high level of conservatism
seems to be in contrast with the modern performance based-design approach, for which a cer-
tain level of yielding in the structure is admitted according to a specific performance. For in-
stance, it would be possible to accept in pipes a certain level of yielding on the condition that
leakage would not occur, and the consequence, for example, to accept a reduction of the
thickness of the flange. The experimental campaign, briefly described in the paper, which has
been performed at the University of Trento in order to evaluate the cyclic behavior of flanged
joints, provided useful information for the design of the flanged joints in the more proper
way. Moreover, it should allow one to link the capacity and the demand for several limit
states.

ACKNOWNLEDGMENTS
This work was carried out with a financial grant from the Research Fund for Coal and Steel of
the European Commission, within the INDUSE project: “Structural Safety of Industrial Steel
Tanks, Pressure Vessels and Piping Systems Under Seismic Loading”, Grant No. RFSR-CT-
2009-00022.

REFERENCES
[1] EN13480-3- Metallic industrial piping - Part 3: Design and calculation, June, 2002
[2] ASME Code for Pressure Piping, B31 - ASME B31.3 (2006)
[3] American Lifelines Alliance - Seismic Design and Retrofit of Piping Systems, 2002.
[4] Guidelines for the Seismic Design of Oil and Gas Piping Systems, Committee on Gas and
Liquid Fuel Lifelines of the ASCE Technical Council on Lifeline Earthquake Engineer-
ing, 1984
[5] Suzuki K., 2006 Earthquake Damage to Industrial Facilities and Development of Seismic
and Vibration Control Technology – Based on Experience from the 1995 Kobe (Hanshin-
Awaji) Earthquake –Journal of Disaster ResearchVol.1 No.2, 2006, Journal of System,
Design and Dynamics, N. 2, 1, 2008. 1-11
[6] Halil Sezen, M.ASCE; and Andrew S. Whittaker, M.ASCE, Seismic Performance of In-
dustrial Facilities Affected by the 1999 Turkey Earthquake, Journal of Performance of
Constructed Facilities, Vol. 20, No. 1, February 1, 2006. ©ASCE, ISSN 0887-
3828/2006/1-28–36/$25.00.
Fabrizio Paolacci, Md Shahin Reza and Oreste S. Bursi

[7] Sassan Eshghi and Mehran Seyed Razzaghi, The Behavior of Special Structures During
the Bam Earthquake of 26 December 2003, JSEE: Special Issue on Bam Earthquake,
197-207.
[8] 2007 Niigata Chuetsu-Oki Japan Earthquake, Reconnaissance Report,
http://www.grmcat.com/images/Niigata-Chuetsu-Oki-Japan-Report.pdf
[9] Azizpour O, Hosseisni M. (2009). A verification of ASCE Recommended Guidelines for
seismic evaluation and design of combination structures in petrochemical facilities, J. of
Applied Sciences.
[10] Asfura, A. and Kiureghian, A. D. (1986), Floor response spectrum method for seismic
analysis of multiply supported secondary systems. Earthquake Engineering & Structural
Dynamics, 14: 245–265. doi: 10.1002/eqe.4290140206
[11] Burdisso, R. A. and Singh, M. P. (1987). Multiply supported secondary systems part I:
Response spectrum analysis. Earthquake Engineering & Structural Dynamics, 15: 53–72.
doi: 10.1002/eqe.4290150105
[12] Saudy A., Ghobarah A., Aziz T.S., A modified CCFS approach for the seismic analysis
of multiply supported MDOF secondary systems, Nuclear Engineering and Design, Vol-
ume 133, Issue 2, March 1992, Pages 183-197, ISSN 0029-5493, DOI: 10.1016/0029-
5493(92)90179-Y
[13] Bakre S.V., Jangid R.S., Reddy G.R., Optimum X-plate dampers for seismic response
control of piping systems, International Journal of Pressure Vessels and Piping, Volume
83, Issue 9, September 2006, Pages 672-685
[14] DeGrassi, G. and Hofmayer, C. (2005). “Seismic analysis of simplified piping systems
for the NUPEC ultimate strength piping test program,” NUREG/CR-6889, by
Brookhaven National Laboratory for the US, Nuclear Regulatory Commission, December,
2005.
[15] MIDAS [2010] - Gen 2010 (v1.1), Midas Information Technology Co. Ltd., 2010
[16] Vamvatsikos, D. and Cornell, C. A. (2002), Incremental dynamic analysis. Earthquake
Engng. Struct. Dyn., 31: 491–514. doi: 10.1002/eqe.141
[17] Okeil A., Tung C., (1996), Effects of ductility on seismic response of piping systems and
their implication on design and qualification, Nuclear Engineering and Design,166: 69-83
[18] Kobayashi H., Ochi Y., (1989), Dynamic response of piping system on rack structure
with gaps and friction, Nuclear Engineering and Design,111: 341-350
[19] Ajaya K. Gupta, Jing-Wen Jaw, A new instructure response spectrum (IRS) method for
multiply connected secondary systems with coupling effects, Nuclear Engineering and
Design, Volume 96, Issue 1, September 1986, Pages 63-80, ISSN 0029-5493, DOI:
10.1016/0029-5493(86)90162-7.
[20] Iervolino I., Galasso C., Cosenza E. (2009). REXEL: computer aided record selection for
code-based seismic structural analysis. Bulletin of Earthquake Engineering, 8:339-362.
DOI 10.1007/s10518-009-9146-1
[21] Reza M.S., Bursi O.S., Paolacci F, Kumar A, (2014), Performance of non-standard bolted
flange joints in industrial piping systems subjected to seismic loading, Journal of Loss
Prevention in Process Industry, DOI: 10.1016/j.jlp.2014.05.011
[22] Liu B., Liu X.J. and Zhang Hong, Strain-based design criteria of pipelines, Journal of
Loss Prevention in the Process Industries 22 (2009) 884-888.
[23] Touboul F., Sollogoub P., Blay N., Seismic behaviour of piping systems with and without
defects: experimental and numerical evaluations, Nuclear Engineering and Design 192
(1999) 243–260.
[24] Touboul F., Blay N., Sollogoub P., Chapuliot S., Enhanced seismic criteria for piping,
Nuclear Engineering and Design 236 (2006) 1–9.
[25] CSA-Z662 (2007). Canadian standards association. Oil and gas pipeline systems. Exdale,
Ontario, Canada.
[26] EN 1993-1-8: Eurocode 3: Design of steel structures - Part 1-8: Design of joints (2005)

View publication stats

You might also like