Piping leakage can occur at T-joint, elbows, valves, or nozzles in nuclear power plants and nonnuclear power plants such as petrochemical plants when subjected to extreme loads and such leakage of piping systems can also lead to fire or explosion. For example, leakage of sodium, toxic gases, or nitrogen in hospitals can cause man-made hazards. The primary objective of this research is to understand the vibration effects due to classical/nonclassical damping with building-piping systems under extreme loads. The current evaluation employed finite-element analysis to calculate the effects of the responses of classically and nonclassically damped building-piping systems. Classical and nonclassical damping matrices for a coupled primary-secondary system were developed based on the Rayleigh equation. A total of 10 selected ground motions were applied to single degree of freedom (SDOF) primary-SDOF secondary (2-DOF coupled) systems in which the ratios of the natural frequencies between the primary and secondary systems ranged between 0.9 and 1.1. It revealed that the vibration effect of nonclassical damping was significant where the natural frequencies of the two systems were nearly tuned. For piping-material nonlinearity, the effects of nonclassical damping on the result forces of piping systems were not significantly different from those of classical damping.
1. Introduction
The vibration of a piping system is significantly influenced by the building in which it is installed, owing to the fact that, typically for such systems, the secondary structure is coupled with the primary structure. The individual components of piping systems and their overall structural integrity are critical to risk management and the continued safe operation of critical facilities such as hospitals and nuclear power plants during extreme loads such as blasts or earthquakes, most especially so as to avoid potential shutdowns. As seen in the recent sequence of events at the Fukushima Daiichi Nuclear power plant in 2011, keeping the core of the reactor cool by long-duration circulation of coolants is a necessity, even after a plant has been shut down due to an earthquake. A piping-system or individual-component failure can result in the loss of the capacity to circulate coolant, which, in turn, can increase the risk of the release of radioactivity. Furthermore, leakage of liquid coolant such as sodium in fast breeder reactor can trigger the internal fires and leakage of oxygen or other gases initiated from piping systems in a hospital can result in fire hazards or explosions [1]. For a typical medium-sized hospital, nonstructural components account for fully 85% of the total construction cost (structural components: 15%), of which portion mechanical, electrical, and piping systems make up 35%, as illustrated in Figure 1 [2]. It has been reported that nonstructural damage, and its results, is the cause of approximately $6.3 billion worth of economic loss [3]. Consequently, many engineers have recognized the significance of nonstructural components and the importance of ensuring that they remain operational, functional, and safe after an earthquake.
Construction costs of medium-sized hospital building.
In the past, calculating the structural response of nonclassically damped systems using uncoupled primary and secondary modal properties was more convenient than calculating the response using completely coupled primary-secondary systems, for example, the studies of the uncoupled primary and secondary systems conducted to evaluate the response of nonclassically damped systems using response spectrum at the base of a primary structure [4, 5]. The assumptions made (e.g., representation of ground motion using Gaussian process models) in developing their formulations, however, led to inaccurate outputs in the low- or high-frequency range. Gupta [6] pointed out that the calculated vibrations of a coupled primary-secondary system were more comprehensive and accurate than those obtained from uncoupled systems. Moreover, Gupta [7] demonstrated that the effect of nonclassical damping is significant when uncoupled systems are (nearly) tuned and when the modal mass ratios are sufficiently small. Furthermore, in recent years, many researchers have recognized the need to control the excessive behavior under seismic ground motions, in order to reduce the seismic damage or system drift of structural and nonstructural components [8–13].
Consequently, this paper reports the results of an evaluation of the vibrations of classical and nonclassical damping for coupled primary-secondary systems including material nonlinearity. Finite-element building-piping models derived from the Open System for Earthquake Engineering Simulation (OpenSees) [14] were used to carry out the Rayleigh equation to develop classical and nonclassical damping matrices for a 2-DOF coupled primary-secondary system. A total of 10 ground motions originally proposed by the ATC-63 project (publication: FEMA P695 [15]) as extreme loading conditions were selected for this study. A numerical model was developed for direct-integration time-history analyses of each ground motion case. The piping models included the effects of piping nonlinearity due to piping components’ nonlinear behaviors. The effect of the mass ratio between building and piping systems also was studied.
2. Analysis of Primary and Secondary Systems
Piping systems installed in typical hospitals or office buildings are comprised of main-pipe lines, piping branches, and supports. An equivalent 2-degree of freedom (DOF) coupled (SDOF primary-SDOF secondary) system is shown in Figure 2. The diagram represents the piping subsystem in a building in which the mass of the frame is assumed to be concentrated at the floor levels or main piping lines, and the floor beams or branches are assumed to be rigid, with the story stiffness provided by the flexure of the columns or piping supports. The typical values of the modal mass ratios in actual building-piping systems are on the order of 0.0001 or lower [16].
SDOF building-SDOF piping system.
2.1. Equation of Motion for Coupled System
The equation of motion for an N-DOF coupled primary-secondary system can be written as follows:(1)MU¨+CU˙+KU=-Mu¨gt,where M is the mass matrix, C is the damping matrix, K is the stiffness matrix, U is the displacement vector with respect to the fixed base, and ug is the ground acceleration.
These matrices and vectors can be expressed in terms of the matrices and vectors of the primary and secondary systems:(2)M=Mp00Ms,C=Cp+CpsCpsCspCs,K=Kp+KpsKpsKspKs,U=UpUs,where those matrices and vectors are denoted by subscripts p and s, respectively, and the matrices Kps and Cps are the stiffness and damping contributions of the secondary system to the connected degrees of freedom in the primary system.
Direct integration of the equation of motion into the time-history analysis enables calculation of the system responses at discrete, usually equally apportioned intervals of time. As seen in (1), displacement, velocity, and acceleration are the basic computation parameters for determination of structural response. The integration algorithms are based on appropriate expressions selected for their relation of the response parameters at the given intervals of time to their historic values.
If the acceleration function in terms of force is represented by a series of straight lines, precise formulations for integration of the equation of motion can be developed in a linear system. With due consideration of the time-dependent acceleration u¨(τ), the variation of acceleration with time can reasonably be assumed to be linear (shown by the dashed line in Figure 3) when a sufficiently small interval of time Δt is used. With the origin location at time τ, u¨(τ) can be written as follows:(3)u¨τ=u¨n+u¨n+1-u¨ntn+1-tnτ.Using this equation’s acceleration function in terms of force, the equation of motion for an N-DOF coupled primary-secondary system can be solved.
Linear representation of acceleration u¨(τ).
2.2. Classical and Nonclassical Damping Matrices
The response of the primary system is dominated by the first-mode frequency of the 2-DOF coupled system, and the response of the secondary system is dominated by the second-mode frequency. Therefore, a classical damping matrix for a 2-DOF coupled system can be evaluated using the Rayleigh equation:(4)Ccl=αMp00Ms+βKp+Ks-Ks-KsKs,where(5)α=2ωiωjωj2-ωi2ωjξp-ωiξs,β=2ωj2-ωi2ωjξs-ωiξp.The transformed damping matrix C¯ for the coupled system is(6)C¯=ϕTCϕ.Using (4) to (6), the transformed damping matrix Ccl¯ for a 2-DOF coupled system is always classical, which is to say, the off-diagonal terms are zeros. The nonclassical damping matrix for a 2-DOF coupled system also can be evaluated using the Rayleigh equation(7)Cp=αpMp000+βpKp000,Cs=αs000Ms+βsKs-Ks-KsKs,Cnc=Cp+Cs,where αp and βs are calculated as(8)αp=2ωiωjωj2-ωi2ωjξp-ωiξp,βp=2ωj2-ωi2ωjξp-ωiξp,αs=2ωiωjωj2-ωi2ωjξs-ωiξs,βs=2ωj2-ωi2ωjξs-ωiξs.Since the damping ratio between the primary system and the secondary system is not the same, the transformed damping matrix Cnc¯ is nonclassical, or, in other words, the off-diagonal terms are nonzero.
3. Ground Motion and Response Spectra
In order to evaluate classically and nonclassically damped coupled primary-secondary systems, a total of 10 ground motions originally proposed by the ATC-63 project, published as FEMA P695 [15], were selected in this study. The records correspond to earthquakes resulting from a variety of fault mechanisms including normal, reverse, and strike-slip fault. The local condition for each record was considered to be either site class C or D, and no more than two records were selected for a single earthquake. Table 1 lists the detailed characteristics of the ground motion records. As is indicated, the data were recorded in the US, Turkey, Iran, and Taiwan. The earthquake magnitude was greater than 6.0 Mw in every case, and the range of peak ground accelerations (PGA) was 0.210 g–0.822 g. To investigate the effect of the frequency tunings between two systems, linear and nonlinear time-history analyses were conducted using building-piping models subjected to the 10 different ground motions. For that purpose, the ground motions were first normalized to a given value of PGA (see Figure 4), so that the effect of the tuning for the different ground motions could be indicated more clearly. Figure 5 plots the response spectra of the total ground motions.
Details of selected ground motions.
Number
Date
Event
Location
Magnitude(Mw)
Dt (s)
1
08/171999
Kocaeli
Arcelik,Turkey
7.5
0.01
2
10/161999
Hector Mine
Hector,USA
7.1
0.01
3
10/151979
Imperial Valley
Delta,USA
6.5
0.005
4
08/171999
Kocaeli
Duzce,Turkey
7.5
0.01
5
09/201999
Chi-Chi
Tcu045,Taiwan
7.6
0.02
6
06/201990
Manjil
Manjil,Iran
7.4
0.005
7
01/171994
Northridge
Beverly Hills,USA
6.7
0.01
8
04/251992
Cape Mendocino
Rio Dell Overpass Ff,USA
7.0
0.005
9
10/181989
Loma Prieta
Gilroy Array,USA
6.9
0.005
10
11/121999
Duzce
Bolu,Turkey
7.1
0.01
Normalized ground motions.
Response spectra.
4. Classical and Nonclassical Damping in Coupled System
Linear and nonlinear direct-integration time-history analyses were performed using SDOF primary-SDOF secondary (2-DOF coupled) finite-element models subjected to each normalized ground motion. The different ratios of the natural frequencies or masses between the primary and secondary systems were considered.
4.1. Response Using Classical or Nonclassical Damping Matrix
For the uncoupled primary and secondary systems, the masses Mp = 90.7 kg and Ms = 0.00907 kg, respectively, were considered. The damping ratios for the primary and secondary systems were ξp = 5% and ξs = 2%, respectively. In general, the 5% damping ratio represents concrete material, and the 2% damping ratio, steel material. To represent the respective stiffness of the primary and secondary systems, Kp = 10,205 N/mm and Ks = 0.827–1.235 N/mm were applied, and the frequencies were nearly tuned with fp = 2.717 Hz and fs = 2.978–2.988 Hz. To further illustrate the significance of nonclassical damping, a parametric study was conducted using SDOF primary-SDOF secondary (2-DOF coupled) systems and varying the tuning between their frequencies within the 0.9–1.1 range. Figure 6 shows the differences in the secondary-system result forces obtained from the nonclassical and classical damping matrices for different frequency-ratio values: specifically, the ratios of secondary to primary system frequency for each ground motion. As seen in the figure, the effect of nonclassical damping was significant only in the systems that had nearly tuned modes. Table 2 lists the result forces for the secondary systems subjected to ground motions.
Result forces for secondary system in case of fs/fp=1.0.
Number
Force (N)
1
552.92
2
448.71
3
1117.07
4
424.80
5
758.64
6
341.06
7
520.78
8
350.02
9
470.27
10
926.48
Differences in scale of 2-DOF coupled system.
4.2. Responses of Two Systems Including Material Nonlinearity
A nonlinear finite-element model for a piping system was developed using OpenSees [17] and a hysteretic material model. The hysteretic behavior of the piping component was characterized by the basic parameters: yield strength Fy, initial stiffness k1, and hardening stiffness k2, as shown in Figure 7, and elasticbeamcolumn element given in Table 3 was also used in OpenSees. Typically, the material property of the piping systems derived from the yield strength was defined by its initial stiffness k1, listed in Table 4. Figure 8 plots the result forces of the vibrations for classically or nonclassically damped coupled primary-secondary systems including piping-material nonlinearity. The PGA for the first ground motion, as multiplied by 0.1, was selected for this study. As indicated in the figure, the result forces in the range of fs/fp<1.05, between the classically and nonclassically damped systems, did not show significant differences using the nonlinear piping system. Even if the fs/fp was greater than 1.00, the result forces between the classically and nonclassically damped systems again showed no significant differences. This was due to the fact that the piping system still remained in the elastic range (see Figure 7). It should be noted that the displacements of the piping system in the nonlinear range showed differences due to the nature of the nonlinear material behavior, even when the result forces with respect to the piping nonlinearity were very close.
Definition of the elasticbeamcolumn element parameters in OpenSees (after Mazzoni et al., 2006 [17]).
Parameters
Definitions
$eleTag
Unique element object tag
$iNode $jNode
End nodes
$A
Cross-sectional area of element
$E
Young’s modulus
$G
Shear modulus
$J
Torsional moment of inertia of cross section
$Iz
Second moment of area about the local z-axis
$Iy
Second moment of area about the local y-axis
$transfTag
Identifier for previously defined coordinate-transformation object
The parameters of coupled building-piping systems.
Parameter
fs/fp=0.9
fs/fp=1.05
Fy
42 N
57 N
k1
826.61 N/m
1125.11 N/m
k2
8.26 N/m
11.25 N/m
Nonlinear material behavior for piping system.
fs/fp = 0.9
fs/fp = 1.05
Linear versus nonlinear responses of piping system.
4.3. Effect of Mass Ratios for Nonlinear Piping Systems
Gupta (1999) earlier found that the effect of nonclassical damping was significant when the modal mass ratios of coupled primary-secondary systems were sufficiently small and those systems remained in the linear elastic range. Figure 9 shows the differences between nonclassical damped systems and classically damped systems for a linear building-piping system as well as a linear building system and a nonlinear piping system. The effect of nonclassical damping was significant when the mass ratios between the primary-secondary systems were on the order of 0.0001 or lower. However, if the piping system was in the nonlinear range, the differences in the result forces were significantly reduced. Next, the vibration effect of between classically/nonclassically damped coupled primary-secondary systems subjected to ground motions will be discussed.
Difference between linear and nonlinear models with respect to mass ratio.
5. Vibration Attenuation of Building-Piping System under Seismic Ground Motions
For vibration control of the building-piping system, the relative effectiveness of the classical and nonclassical damping system was considered in this study. In particular, five simulations corresponding to the different frequency ratio between the systems were carried out to measure the acceleration time histories at the top floor of the system. A vibration comparison between classical and nonclassical damped systems conducted in OpenSees was described in Figure 10. Also, it showed the quantitative comparison of the acceleration capacities of the systems. It can be seen that there was a similar manner in all of the acceleration time histories: nonclassical damped building-piping system was more conservative than classical coupled building-piping system, in terms of reducing the floor accelerations through the ground acceleration using the 1999 Chi-Chi earthquake. Furthermore, the nonclassical damped system showed the vibration attenuation about 20% with 0.95, 0.99, and 1.05 frequency ratios. In terms of frequency ratios—0.99 and 1.01, however, the vibration attenuation was significantly influenced by nonzero off-diagonal term, resulting in a decrease of more than 50%. Consequently, the vibration of the coupled building-piping system was sensitive to the mechanical properties such as mass and stiffness ratio of the primary and secondary system.
Acceleration time histories with different frequencies in classically and nonclassically damped coupled building-piping systems.
fs/fp = 0.95
fs/fp = 0.99
fs/fp = 1.00
fs/fp = 1.01
fs/fp = 1.05
6. Conclusions
The present study evaluated the effects of the seismic responses of classically and nonclassically damped coupled primary-secondary piping systems, including nonlinearity behavior. A finite-element building-piping model using OpenSees was employed to develop, using the Rayleigh equation, classical and nonclassical damping matrices for a 2-DOF coupled primary-secondary system. A total of 10 selected ground motions were applied to SDOF primary-SDOF secondary (2-DOF coupled) systems in which the ratios of the natural frequencies between the primary and secondary systems ranged between 0.9 and 1.1. Linear and nonlinear direct-integration time-history analyses were performed for each ground motion case. It was shown that the effect of nonclassical damping was significant in systems that have nearly tuned modes. For piping-material nonlinearity, the effects of nonclassical damping on the result forces of the piping system were not significantly different from those of classical damping. In fact, they were insignificant, even though the modal mass ratios were sufficiently small. However, it should be noted that, even when the result forces are very close, piping-system displacements in the nonlinear range can show differences due to the nature of nonlinear material behavior. Additionally, in the simulations using classically and nonclassically damped systems, it was interesting to find that nonzero off-diagonal term in nonclassical damping matrices was very effective to the vibration attenuation, especially, coupled building-piping system with frequency range of 0.95 to 1.05. Also, the probabilistic risk assessment of coupled building-piping systems as the analysis of the vibration attenuation must be achieved, in order to evaluate the probability of failure with respect to material uncertainty and ground motion uncertainty using classical and nonclassical damping matrices.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work was supported by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MEST) (no. 2012-0008762).
RyuY. H.GuptaA.JuB. S.JungW. Y.A framework for evaluating the probability of fire hazard due to leakage of liquid coolant in fast breeder reactorProceedings of the 23rd International Conference on Structural Mechanics in Reactor Technology Conference (SMiRT '15)August 2015Manchester, UKFEMA 5772007Washington, DC, USAFederal Emergency Management AgencyKircherC. A.It makes dollars and sense to improve nonstructural system performanceProceedings of the Seminar on Seismic Design, Performance, and Retrofit of Nonstructural Components in Critical FacilitiesOctober 2003Newport Beach, Calif, USA109120ATC-29-2SinghM. P.ChuS. L.Stochastic considerations in seismic analysis of structures19764329530710.1002/eqe.42900403092-s2.0-0016871298VillavardeR.NewmarkN. M.Seismic response of light attachments to buildings1980469Urbana, Ill, USAUniversity of IllinoisGuptaA. K.1992Boca Raton, Fla, USACRC PressGuptaA.Significance of nonclassical damping in coupled system analysisProceedings of the 15th Structural Mechanics in Reactor Technology Conference (SMiRT '99)August 1999Seoul, Republic of KoreaParkK.KimD.YangD.DaekiJ.HaI.KimS.A comparison study of conventional construction methods and outrigger damper system for the compensation of differential column shortening in high-rise buildings201010431732410.1007/bf032158402-s2.0-80053289275AghajanianS.BaghiH.AminiF.SamaniM. Z.Optimal control of steel structures by improved particle swarm201414222323010.1007/s13296-014-2003-32-s2.0-84904296746OhS.-H.SongS.-H.LeeS.-H.KimH.-J.Seismic response of base isolating systems with U-shaped hysteretic dampers201212228529810.1007/s13296-012-2011-02-s2.0-84866000152ShinD.-H.KimH.-J.Probabilistic assessment of structural seismic performance influenced by the characteristics of hysteretic energy dissipating devices201514469771010.1007/s13296-014-1202-22-s2.0-84921047384MaC. F.ZhangY. H.TanP.ZhouF. L.Seismic response of base-isolated high-rise buildings under fully nonstationary excitation201420141140146910.1155/2014/4014692-s2.0-84904135140FotiD.On the seismic response of protected and unprotected middle-rise steel frames in far-field and near-field areas201420141139387010.1155/2014/393870OpenSeesOpen System for Earthquake Engineering Simulation2006, http://opensees.berkeley.edu/FEMA P6952009Washington, DC, USAFederal Emergency Management AgencyBoseM. K.2001Raleigh, NC, USANorth Carolina State UniversityMazzoniS.McKennaF.ScottM. H.FenvesG. L.OpenSees command language manual2006, http://opensees.berkeley.edu/