Research Article  Open Access
Vibration Effects of Nonclassically Damped BuildingPiping Systems Subjected to Extreme Loads
Abstract
Piping leakage can occur at Tjoint, 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 manmade hazards. The primary objective of this research is to understand the vibration effects due to classical/nonclassical damping with buildingpiping systems under extreme loads. The current evaluation employed finiteelement analysis to calculate the effects of the responses of classically and nonclassically damped buildingpiping systems. Classical and nonclassical damping matrices for a coupled primarysecondary system were developed based on the Rayleigh equation. A total of 10 selected ground motions were applied to single degree of freedom (SDOF) primarySDOF secondary (2DOF 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 pipingmaterial 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 longduration circulation of coolants is a necessity, even after a plant has been shut down due to an earthquake. A pipingsystem or individualcomponent 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 mediumsized 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.
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 primarysecondary 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 highfrequency range. Gupta [6] pointed out that the calculated vibrations of a coupled primarysecondary 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 primarysecondary systems including material nonlinearity. Finiteelement buildingpiping 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 2DOF coupled primarysecondary system. A total of 10 ground motions originally proposed by the ATC63 project (publication: FEMA P695 [15]) as extreme loading conditions were selected for this study. A numerical model was developed for directintegration timehistory 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 mainpipe lines, piping branches, and supports. An equivalent 2degree of freedom (DOF) coupled (SDOF primarySDOF 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 buildingpiping systems are on the order of 0.0001 or lower [16].
2.1. Equation of Motion for Coupled System
The equation of motion for an NDOF coupled primarysecondary system can be written as follows:where is the mass matrix, is the damping matrix, is the stiffness matrix, is the displacement vector with respect to the fixed base, and is the ground acceleration.
These matrices and vectors can be expressed in terms of the matrices and vectors of the primary and secondary systems:where those matrices and vectors are denoted by subscripts and , respectively, and the matrices and 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 timehistory 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 timedependent acceleration , 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 is used. With the origin location at time , can be written as follows:Using this equation’s acceleration function in terms of force, the equation of motion for an NDOF coupled primarysecondary system can be solved.
2.2. Classical and Nonclassical Damping Matrices
The response of the primary system is dominated by the firstmode frequency of the 2DOF coupled system, and the response of the secondary system is dominated by the secondmode frequency. Therefore, a classical damping matrix for a 2DOF coupled system can be evaluated using the Rayleigh equation:whereThe transformed damping matrix for the coupled system isUsing (4) to (6), the transformed damping matrix for a 2DOF coupled system is always classical, which is to say, the offdiagonal terms are zeros. The nonclassical damping matrix for a 2DOF coupled system also can be evaluated using the Rayleigh equationwhere and are calculated asSince the damping ratio between the primary system and the secondary system is not the same, the transformed damping matrix is nonclassical, or, in other words, the offdiagonal terms are nonzero.
3. Ground Motion and Response Spectra
In order to evaluate classically and nonclassically damped coupled primarysecondary systems, a total of 10 ground motions originally proposed by the ATC63 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 strikeslip 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 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 timehistory analyses were conducted using buildingpiping 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.

4. Classical and Nonclassical Damping in Coupled System
Linear and nonlinear directintegration timehistory analyses were performed using SDOF primarySDOF secondary (2DOF coupled) finiteelement 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 = 90.7 kg and = 0.00907 kg, respectively, were considered. The damping ratios for the primary and secondary systems were = 5% and = 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, = 10,205 N/mm and = 0.827–1.235 N/mm were applied, and the frequencies were nearly tuned with = 2.717 Hz and = 2.978–2.988 Hz. To further illustrate the significance of nonclassical damping, a parametric study was conducted using SDOF primarySDOF secondary (2DOF coupled) systems and varying the tuning between their frequencies within the 0.9–1.1 range. Figure 6 shows the differences in the secondarysystem result forces obtained from the nonclassical and classical damping matrices for different frequencyratio 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.

4.2. Responses of Two Systems Including Material Nonlinearity
A nonlinear finiteelement 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 , initial stiffness , and hardening stiffness , 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 , listed in Table 4. Figure 8 plots the result forces of the vibrations for classically or nonclassically damped coupled primarysecondary systems including pipingmaterial 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 , between the classically and nonclassically damped systems, did not show significant differences using the nonlinear piping system. Even if the 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.


(a) = 0.9
(b) = 1.05
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 primarysecondary 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 buildingpiping 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 primarysecondary 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 primarysecondary systems subjected to ground motions will be discussed.
5. Vibration Attenuation of BuildingPiping System under Seismic Ground Motions
For vibration control of the buildingpiping 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 buildingpiping system was more conservative than classical coupled buildingpiping system, in terms of reducing the floor accelerations through the ground acceleration using the 1999 ChiChi 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 offdiagonal term, resulting in a decrease of more than 50%. Consequently, the vibration of the coupled buildingpiping system was sensitive to the mechanical properties such as mass and stiffness ratio of the primary and secondary system.
(a) = 0.95
(b) = 0.99
(c) = 1.00
(d) = 1.01
(e) = 1.05
6. Conclusions
The present study evaluated the effects of the seismic responses of classically and nonclassically damped coupled primarysecondary piping systems, including nonlinearity behavior. A finiteelement buildingpiping model using OpenSees was employed to develop, using the Rayleigh equation, classical and nonclassical damping matrices for a 2DOF coupled primarysecondary system. A total of 10 selected ground motions were applied to SDOF primarySDOF secondary (2DOF 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 directintegration timehistory 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 pipingmaterial 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, pipingsystem 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 offdiagonal term in nonclassical damping matrices was very effective to the vibration attenuation, especially, coupled buildingpiping system with frequency range of 0.95 to 1.05. Also, the probabilistic risk assessment of coupled buildingpiping 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. 20120008762).
References
 Y. H. Ryu, A. Gupta, B. S. Ju, and W. Y. Jung, “A framework for evaluating the probability of fire hazard due to leakage of liquid coolant in fast breeder reactor,” in Proceedings of the 23rd International Conference on Structural Mechanics in Reactor Technology Conference (SMiRT '15), Manchester, UK, August 2015. View at: Google Scholar
 FEMA 577, Risk Management Series Design Guide for Improving Hospital Safety in Earthquakes, Floods, and High Winds, Federal Emergency Management Agency, Washington, DC, USA, 2007.
 C. A. Kircher, “It makes dollars and sense to improve nonstructural system performance,” in Proceedings of the Seminar on Seismic Design, Performance, and Retrofit of Nonstructural Components in Critical Facilities, ATC292, pp. 109–120, Newport Beach, Calif, USA, October 2003. View at: Google Scholar
 M. P. Singh and S. L. Chu, “Stochastic considerations in seismic analysis of structures,” Earthquake Engineering and Structural Dynamics, vol. 4, no. 3, pp. 295–307, 1976. View at: Publisher Site  Google Scholar
 R. Villavarde and N. M. Newmark, “Seismic response of light attachments to buildings,” Structural Research Series no. 469, University of Illinois, Urbana, Ill, USA, 1980. View at: Google Scholar
 A. K. Gupta, Response Spectrum Method in Seismic Analysis and Design of Structures, CRC Press, Boca Raton, Fla, USA, 1992.
 A. Gupta, “Significance of nonclassical damping in coupled system analysis,” in Proceedings of the 15th Structural Mechanics in Reactor Technology Conference (SMiRT '99), Seoul, Republic of Korea, August 1999. View at: Google Scholar
 K. Park, D. Kim, D. Yang, J. Daeki, I. Ha, and S. Kim, “A comparison study of conventional construction methods and outrigger damper system for the compensation of differential column shortening in highrise buildings,” International Journal of Steel Structures, vol. 10, no. 4, pp. 317–324, 2010. View at: Publisher Site  Google Scholar
 S. Aghajanian, H. Baghi, F. Amini, and M. Z. Samani, “Optimal control of steel structures by improved particle swarm,” International Journal of Steel Structures, vol. 14, no. 2, pp. 223–230, 2014. View at: Publisher Site  Google Scholar
 S.H. Oh, S.H. Song, S.H. Lee, and H.J. Kim, “Seismic response of base isolating systems with Ushaped hysteretic dampers,” International Journal of Steel Structures, vol. 12, no. 2, pp. 285–298, 2012. View at: Publisher Site  Google Scholar
 D.H. Shin and H.J. Kim, “Probabilistic assessment of structural seismic performance influenced by the characteristics of hysteretic energy dissipating devices,” International Journal of Steel Structures, vol. 14, no. 4, pp. 697–710, 2015. View at: Publisher Site  Google Scholar
 C. F. Ma, Y. H. Zhang, P. Tan, and F. L. Zhou, “Seismic response of baseisolated highrise buildings under fully nonstationary excitation,” Shock and Vibration, vol. 2014, Article ID 401469, 11 pages, 2014. View at: Publisher Site  Google Scholar
 D. Foti, “On the seismic response of protected and unprotected middlerise steel frames in farfield and nearfield areas,” Shock and Vibration, vol. 2014, Article ID 393870, 11 pages, 2014. View at: Publisher Site  Google Scholar
 OpenSees, “Open System for Earthquake Engineering Simulation,” 2006, http://opensees.berkeley.edu/. View at: Google Scholar
 FEMA P695, Quantification of Building Seismic Performance Factors, Federal Emergency Management Agency, Washington, DC, USA, 2009.
 M. K. Bose, Nonclassical damping properties and modal correlation coefficient for dynamic analysis of structures [Ph.D. thesis], North Carolina State University, Raleigh, NC, USA, 2001.
 S. Mazzoni, F. McKenna, M. H. Scott, and G. L. Fenves, “OpenSees command language manual,” 2006, http://opensees.berkeley.edu/. View at: Google Scholar
Copyright
Copyright © 2016 YongHee Ryu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.