Seismic Design Criteria of Refinery Piping Systems: June 2014
Seismic Design Criteria of Refinery Piping Systems: June 2014
Seismic Design Criteria of Refinery Piping Systems: June 2014
net/publication/263131870
CITATIONS READS
12 5,710
3 authors:
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.
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
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.
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
x
y
(1)
(2)
(3)
(4)
(5)
(6)
(7)
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:
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.
code SF SSF
EN13480 ---- ----
ASME B31.3
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
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.
Figure 10 First three vibration modes of the pipe plus pipe-rack model
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
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.
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.
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:
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.
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.
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.
σ ≤ 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.
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
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.
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)