Next Article in Journal
Effect of Graphene Nanoplatelet on the Carbonation Depth of Concrete under Changing Climate Conditions
Next Article in Special Issue
Seismic Damage Probability Assessment of Isolated Girder Bridges Based on Performance under Near-Field Earthquakes
Previous Article in Journal
Development of an Analytical Model to Determine the Heat Fluxes to a Structural Element Due to a Travelling Fire
Previous Article in Special Issue
Time-Dependent Damage Estimation of a High-Rise Steel Building Equipped with Buckling-Restrained Brace under a Series of Earthquakes and Winds
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Shaking Table Tests to Validate Inelastic Seismic Analysis Method Applicable to Nuclear Metal Components

1
Korea Atomic Energy Research Institute, Daejeon 34057, Korea
2
SMARTEK E&C, Daejeon 34013, Korea
3
Department of Nuclear Engineering, Sejong University, Seoul 05006, Korea
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(19), 9264; https://doi.org/10.3390/app11199264
Submission received: 2 September 2021 / Revised: 23 September 2021 / Accepted: 29 September 2021 / Published: 6 October 2021
(This article belongs to the Special Issue Seismic Assessment and Design of Structures)

Abstract

:
The main purpose of this study is to perform shaking table tests to validate the inelastic seismic analysis method applicable to pressure-retaining metal components in nuclear power plants (NPPs). To do this, the test mockup was designed and fabricated to be able to describe the hot leg surge line nozzle with a piping system, which is known to be one of the seismically fragile components in nuclear steam supply systems (NSSS). The used input motions are the displacement time histories corresponding to the design floor response spectrum at an elevation of 136 ft in the in-structure building in NPPs. Two earthquake levels are used in this study. One is the design-basis safe shutdown earthquake level (SSE, PGA = 0.3 g) and the other is the beyond-design-basis earthquake level (BDBE, PGA = 0.6 g), which is linearly scaled from the SSE level. To measure the inelastic strain responses, five strain gauges were attached at the expected critical locations in the target nozzle, and three accelerometers were installed at the shaking table and piping system to measure the dynamic responses. From the results of the shaking table tests, it was found that the plastic strain response at the target nozzle and the acceleration response at the piping system were not amplified by as much as two times the input earthquake level because the plastic behavior in the piping system significantly contributed to energy dissipation during the seismic events. To simulate the test results, elastoplastic seismic analyses with the well-known Chaboche kinematic hardening model and the Voce isotropic hardening model for Type 316 stainless steel were carried out, and the results of the principal strain and the acceleration responses were compared with the test results. From the comparison, it was found that the inelastic seismic analysis method can give very reasonable results when the earthquake level is large enough to invoke plastic behavior in nuclear metal components.

1. Introduction

In the current seismic design of nuclear power plants (NPPs), the peak ground acceleration (PGA) level of safe shutdown earthquake (SSE) is generally taken to be 0.3 g as a design-basis earthquake. Recently, there have been attempts to increase the PGA level above 0.3 g to enhance the seismic safety of NPPs. Furthermore, after the Fukushima NPP accident, the consideration of the beyond-design-basis earthquake (BDBE) became one of the big issues at the design stage, especially for safety-related nuclear facility components having small seismic margins. In fact, it is true that increasing the value of the SSE level is a big burden for the seismic design of NPPs. There have been many research and development efforts to resolve BDBE issues for the seismic design of NPPs. As one of the hardware approaches to accommodate large earthquakes, seismic isolation design has been studied by many countries for a whole nuclear island building isolation [1,2] or an individual facility component isolation [3,4]. Additionally, seismic energy absorbers such as tuned mass and dampers (TMDs) have been studied for application to seismic fragile piping systems or any other safety-related facilities [5].
Recently, as one of the approaches against large earthquakes, studies on the new seismic design acceptance criteria for large earthquakes have been actively conducted for the application of nuclear metal components [6,7,8,9]. The main concept of this design criteria provides strain-based design limits to protect the pressure-retaining nuclear metal components against an inelastic strain-induced failure mode. Actually, current design criteria in nuclear codes and standards such as ASME BPVC III [10] and RCC-MR [11] provide the stress-based design limits for the elastically calculated stress values. Therefore, they may not provide adequate limits for actual seismic failure modes in cases when the reversing dynamic plastic responses are significant.
In the previous studies on the new strain-based seismic design criteria, the feasibility of this approach was investigated for the nuclear pipe and nozzle. The most important thing when applying this approach is to calculate accurate inelastic strain values using inelastic seismic analysis [12,13]. Then, seismic shaking table tests are inevitably required in order to validate the inelastic seismic analysis method. There has been research on seismic shaking table tests for the piping system [14,15].
In this study, the test mockup for simulating the hot leg surge line nozzle, which is known as one of the fragile components in nuclear steam supply systems (NSSS), is designed and tested on a shaking table for two earthquake levels of PGA = 0.3 and 0.6 g. For validation of the inelastic seismic analysis method, the finite element model using the ANSYS program [16] is established with material constitutive equations of the Chaboche kinematic hardening model and the Voce isotropic hardening model for Type 316 stainless steel. Through the investigation of the test results and the comparison of seismic time history responses such as the principal strain and accelerations between the measured data and the inelastic seismic analysis results, the inelastic seismic analysis method for large earthquakes was validated and confirmed to be useful for the seismic design of nuclear metal components.

2. Design of Seismic Test Mockup

2.1. Configuration and Dimensions

In this study, the seismic test mockup is designed to be able to simulate the actual nuclear component installed in the in-structure building. This will be consistent with the purpose of validating the strain-based seismic design method developed in previous studies [9,13]. To do this, the metal component nozzle, which is known as one of the seismically fragile components in nuclear steam supply systems (NSSS), was selected as the reference for the test mockup. Specifically, the hot leg surge line nozzle connecting the hot-leg pipe and the surge line piping system, which is installed at the elevation of 136 ft in the in-structure building, is referenced in this test mockup design.
The main design concept of the test mockup is to configure the piping layout, resulting in the maximum accumulated plastic strain at the safe-end region connecting the nozzle and the piping system in the same condition as the actual hot-leg surge line nozzle in NSSS. Figure 1 presents the designed overall configuration and seismic directional axis of the test mockup. As shown in the figure, the test mockup consists of a target nozzle, horizontal and vertical pipes, one elbow, and two added masses. The total height of the test mockup is 1132 mm, and the total mass is about 158 kg. Table 1 reveals the summary of the design parameter values of the test mockup.
The used material of the test mockup is Type 316 stainless steel, which has an elastic modulus of 190 GPa, a density of 7970 kg/m3, and a Poison’s ratio of 0.27 at room temperature.

2.2. Analysis Model and Dynamic Characteristics of Test Mockup

To derive an appropriate plastic strain in the nozzle part, enough to validate the inelastic seismic analysis method, the test mockup is designed to have dominant dynamic characteristics resonant in the peak spectral frequency band of the seismic input motions. Figure 2 presents the target required response spectrum, representing those at the elevation of 136 ft in the in-structure building (ISB) of NSSS. As shown in the figure, the peak spectral frequency band is 9 to 12 Hz for the EW (east–west) direction, 5 to 20 Hz for the NS (north–south) direction, and 15 to 21 Hz for the V (vertical) direction.
The dynamic characteristics of the test mockup are investigated by the modal analysis by using the ANSYS [16] commercial finite element program. Figure 3 shows the finite element seismic analysis model of the test mockup. As shown in the model, the element type of the nozzle is a three-dimensional solid element of SOLID185. For the piping system, the three-dimensional pipe elements of PIPE289 and ELBOW290 are used for the straight run section and for the elbow, respectively. To rigidly couple the solid element of the nozzle and the line element of the pipe, the element type of MPC184 is used. The added masses in the test mockup are modeled with a structural mass element of MASS21, including the effects of the rotary moments of inertia.
Related to the boundary conditions, it is assumed that the seismic input motions are applied to the nozzle bottom surface as the design response spectrum is specified at the hot-leg piping system. Then, the modal analysis is carried out with the fixed condition of the nozzle bottom surface.
Table 2 presents the summary results of the modal analysis. For the horizontal directions, the dominant fundamental frequency is 11.74 Hz for the NS direction and 11.99 Hz for the EW direction. Additionally, the dominant natural frequency in the vertical direction is 27.63 Hz, slightly apart from the vertical target peak response frequency band. From the modal analysis, it is confirmed that the dynamic characteristics of the test mockup will meet the intended design goal, deriving the sufficiently inelastic strains at the nozzle.
Figure 4 presents the results of the mode shapes for dominant natural frequencies. From the results, it is expected that Mode 1 and Mode 3 will contribute to the torsional seismic response at the nozzle, and Mode 2, Mode 4, and Mode 5 will contribute to the bending seismic responses at the nozzle.
The deadweight effect of the added masses attached to the pipe was examined and included in the modal analysis, but it was found that it has minimal effects on the dynamic characteristics of the test mockup.

3. Seismic Shaking Table Tests

3.1. Performance of Shaking Table Test Facility

Table 3 shows the performance information of the MTS shaking table (2209 × 2173 mm) used in tests. As shown in the table, the test facility can perform the maximum excitation level of 4.9 g for the horizontal direction and 6.0 g for the vertical direction with the loaded maximum of 2000 kg on the table. The frequency range for excitation is 0.8 to 100 Hz.

3.2. Description of Test System and Sensors

As shown in Figure 5, the nozzle bottom flange is connected to the adapter flange by bolts, and then the adapter flange is installed on the shaking table by the bolted connection. Then, it can be assumed that the bolted connection of the test mockup with the shaking table is strong enough to assure the fixed boundary condition, as intended. In order to apply appropriate seismic inertia load on the nozzle, two added masses are attached to the pipes. These masses, Added Mass-1 and Added Mass-2, are composed of eight and four plates combined by bolts, respectively.
For the measurement of seismic strain responses, five strain gauges are attached to the nozzle safe end region where the maximum plastic strain responses are expected. They are almost equally spaced in the circumferential direction at the location of the safe-end region, as shown in Figure 6a. Figure 6b shows the pattern of the general-purpose foil strain gauge type used in this test. Table 4 presents the specification of the used strain gauge.
As shown in Figure 5, three accelerometers are installed at the test system (one on the shaking table and two on the piping system) to measure the seismic acceleration responses. Since the seismic inelastic strain responses at the nozzle depend on what the seismic responses are at the piping system, it is important to investigate the acceleration responses. The accelerometer installed on the shaking table is to measure the input seismic motions needed for the inelastic seismic analyses. Table 5 presents the specifications of the used accelerometers.

3.3. Results of Shaking Table Tests

3.3.1. Strain Time History Seismic Responses

From the measured strain data (εa, εb, and εc) obtained from the triaxial 0°–45°–90° rosette strain gauges, as conceptually described in Figure 7, the maximum principal strain (εmax), the minimum principal strain (εmin), the maximum shear strain (τmax), and the directional angle of the principal strain (θ) can be calculated using the equations, as follows [17]:
ε max = 1 2 ε a + ε c + 2 ε a ε b 2 + ε b ε c 2
ε min = 1 2 ε a + ε c 2 ε a ε b 2 + ε b ε c 2
τ max = 2 ε a ε b 2 + ε b ε c 2
θ = 1 2   tan 1 2 ε b ε a ε c ε a ε c
In the above equation, the angle, θ, is the angle of the maximum principal strain to the εa axis when εa > εc or the angle of the minimum principal strain to the εa axis when εa < εc.
The accurate coordinates of the five strain gauges attached at the nozzle region and their corresponding node numbers in the seismic analysis model of Figure 3 are listed in Table 6.
Figure 8 and Figure 9 present the test results of the principal strain time history responses for PGA = 0.3 and 0.6 g, respectively. These are calculated from the measured strains of εa, εb, and εc using Equations (1) and (2). As shown in the figures, the strain response time histories of the maximum principal strain and the minimum principal strain are almost symmetric. This means that the seismic strain responses at the nozzle have fully reversing characteristics with almost the same tension and compression behavior. Among the strain gauges, we can see that SG-5 shows the most significant seismic strain responses.

3.3.2. Acceleration Time History Seismic Responses

Figure 10, Figure 11 and Figure 12 present the measured acceleration time history responses for PGA = 0.3 g. As shown in Figure 10, measured at the shaking table, the maximum values are 1.74 g for EW, 1.52 g for NS, and 1.92 g for the vertical direction. These measured data will be used for the validation of the inelastic seismic analysis as input motions.
Table 7 presents the test results of the zero-period acceleration (ZPA) obtained from the shaking table tests of two scale input motions (PGA = 0.3 g and PGA = 0.6 g) and the ZPA ratio of the two cases for each exciting direction. As shown in the table, the ZPA values of target input motions for PGA = 0.6 g are linearly set to be two times those of PGA = 0.3 g. Then, the ZPA ratio measured at the shaking table is almost 2.0, the same as the target ZPA ratio. However, the ZPA ratio values measured at the end of the pipe are much less than 2.0. These results mean that significant plastic behavior might be occurring in the piping system during the PGA = 0.6 g scale seismic shaking table test. Since this plastic behavior in the piping system has a role in energy dissipation in the seismic responses, the acceleration responses do not increase linearly two times even if the input motions increase two times. From these results, it can be seen that in the case of an earthquake level large enough to exceed the design level, a more accurate seismic response can be obtained using the inelastic seismic analysis method, and the effect of reducing the seismic response can be obtained. This will be discussed in the results of the inelastic analysis in the section below.

4. Validation of Inelastic Seismic Analysis

To validate the inelastic seismic analysis method, which can be used for strain-based seismic design, inelastic seismic time history analyses were performed, and their results were compared with those of the tests.

4.1. Analysis Modeling

4.1.1. Dynamic Characteristics

The used finite element model for the test mockup is shown in Figure 3. To confirm the dynamic characteristics of the analysis model, resonance searching tests were carried out, with the random input motions having the frequency range of 1.0 to 100.0 Hz. Figure 13 presents the frequency response functions measured by Accelerometer-3 (shown in Figure 5).
As shown in Figure 13, the measured dominant resonance frequencies are about 11.75 Hz for the 1st mode and 28 Hz for the 2nd mode. Table 8 presents the comparison of the resonance frequencies between tests and analyses.
As shown in Table 8, the dominant modes for EW and NS directions are almost the same in the tests but slightly different in the modal analysis results. In the test results, the frequency of Mode 1 is the same in both directions because it is not easy to separate such a closed mode in the test. However, it is confirmed that the used seismic analysis model describes the dynamic characteristics of the test mockup in good agreement.

4.1.2. Structural Damping Value

To be used for inelastic seismic analyses, the actual structural damping value corresponding to the test mockup is required. As shown in Figure 13, the test mockup is considered to be a very light damped system with dominantly well-separated modes. Therefore, the well-known half-power bandwidth method [18] can be used to identify the structural damping ratio. Figure 14 illustrates the concept for the dominant first natural frequency of the EW direction. In this method, it is assumed that half the total power of dissipation in this mode occurs in the frequency band between f1 and f2, where f1 and f2 are the frequencies corresponding to an amplitude of f c / 2 .
As illustrated in Figure 14, the critical damping ratio, ξ can be approximately determined by the following relationship;
ξ = f 1 f 2 2 f c
From Equation (5), the obtained critical damping values based on the first natural frequency mode for each direction are 0.5% for EW, 0.4% for NS, and 0.8% for the vertical direction. Then, the averaged value of 0.57% for the three directions is used in this study.

4.1.3. Inelastic Material Model

For inelastic seismic analyses, the well-known constitutive equations of Chaboche’s kinematic hardening model [19,20] is used as follows:
α ˙ i j     =     k = 1 3 2 3 C k ε ˙ i j p         γ k α i j k p ˙
where α ˙ i j and p ˙ indicate the revolution of back stress and an accumulated plastic strain, respectively. Additionally, Ck and γk (k = 1~3) are material constants to be used in the ANSYS program.
For the isotropic hardening model, the inelastic Voce model [21] is used as follows:
R ˙     =     b   Q         R   p ˙
where R ˙ indicates the revolution of drag stress. Additionally, b and Q are material constants.
Table 9 presents the used material constants required in Equations (6) and (7) for Type 316 stainless steel [22].

4.1.4. Seismic Input Motions

Figure 15 presents the displacement time history input motions used for the inelastic seismic analyses, which are measured at the shaking table. The correlation coefficients are 0.015 between EW and NS, 0.0015 between NS and V and 0.0003 between V and EW, which are much less than the criteria value of 0.16 required for the seismic input motions [23]. Therefore, since the independence of the input motions is guaranteed, the inelastic seismic analyses can be performed by applying them simultaneously.
The time interval (Δt) used in the inelastic seismic analysis is 3.91 ms; then, the cut-off frequency (1/2Δt), 128 Hz, can sufficiently cover the seismic cut-off frequency of 33 Hz in NPP seismic design. The total analysis duration is 32 s. Figure 16 presents the test response spectrum (TRS) calculated from the measured shaking table motions of Figure 15. Compared with the target required response spectrum in Figure 2, it is confirmed that the TRS envelops the target response spectrum.

4.2. Results of Inelastic Seismic Analyses

4.2.1. Validation of Strain Time History Responses

Figure 17 presents the inelastic seismic time history analysis results, representing the accumulated equivalent plastic strain distribution at the end time of PGA = 0.3 g in the nozzle. As shown in the results, we can see that the maximum seismic strain responses occur in the region of the nozzle’s safe end, as expected in the test mockup design. This result is considered to be dominantly caused by the bending motions of the piping system. In addition, the frequencies of Mode 1 (11.74 Hz, torsional) and Mode 2 (11.99, bending) are in a closed mode. It is judged that these closed bending and torsional modes strongly influence the determination of the location where the maximum strain response occurs.
If we look more closely in Figure 17b, which shows the inelastic strain distribution in the cross-section where the maximum strain occurs, the largest inelastic strain occurs along the (EW, NS) direction. From the sectional view, it is expected that the largest seismic responses will occur at the SG-5 location among the strain gauges in tests.
Figure 18 presents the material hysteretic responses at the location of the maximum accumulated equivalent plastic strain in Figure 17.
To compare the inelastic seismic analysis results with the test results, the maximum and minimum principal strains are calculated at each node, corresponding to the locations of the strain gauges, as shown in Figure 17b. Figure 19 and Figure 20 present the analysis results of the principal strains for PGA = 0.3 and 0.6 g, respectively. As expected from the inelastic seismic analysis, the largest strain responses occur at SG-5.
When compared, the seismic strain time history responses of Figure 19 and Figure 20 with the test results of Figure 8 and Figure 9, the overall strain wave shapes and amplitudes are similar. Specifically, we can see that the largest strain responses occur at the same location of SG-5 in the tests and analyses, as expected in Figure 17. In general, the maximum accumulated equivalent plastic strain, which can be a failure mode in strain-based seismic design criteria, occurs at the location where the maximum value of the time history response occurs, but this is not always the case.
Table 10 presents the comparison results of the maximum principal strain values of the time history responses between the tests and inelastic seismic analyses for the case of PGA = 0.3 g and PGA = 0.6 g. As shown in the table, the results of the inelastic strain responses are in good agreement with the test results.

4.2.2. Validation of Acceleration Time History Responses

The investigation of the acceleration time history responses at the pipe is important in pointing out that the seismic strain responses in the nozzle depend on the piping’s seismic behavior. Figure 21 and Figure 22 present the acceleration time histories at the locations of Accelerometer-2 and Accelerometer-3 (see Figure 5). As shown in the figures, the overall response waveforms are very similar to the test results of Figure 11 and Figure 12. The largest acceleration response occurs in the NS (Z) direction. This will dominantly excite Mode 1 and Mode 3, invoking torsional responses at the nozzle.
Table 11 presents the comparison of ZPA values between the tests and analyses at the pipe end (Acceleration-3). The ZPA value corresponds to the maximum amplitude in the acceleration time history responses. As shown in the table, the analysis results reveal slightly larger values compared with the test results, but we can see that overall response characteristics are in good agreement.
Figure 23 presents the comparison results of the response spectrum, calculated from the acceleration time history responses at the pipe end (Accelerometer-3). As shown in the figure, the results obtained from the inelastic seismic analysis reveal good agreement with those of the tests.

5. Conclusions

In this study, shaking table tests are performed to validate the inelastic seismic analysis method applicable to the nuclear metal components. To do this, a test mockup, which can simulate the hot leg surge line nozzle, known as one of the nuclear seismic fragile components in NSSS, was designed and tested on the shaking table with seismic input motions corresponding to design-basis earthquake (PGA = 0.3 g) and beyond-design-basis earthquake (PGA = 0.6 g) levels. To validate the inelastic seismic analysis method, detailed comparisons of seismic responses between the tests and inelastic seismic analyses were carried out, especially for the seismic responses of the principal strains at the nozzle and the accelerations at the pipe. From this study, some meaningful conclusions have been derived, as follows:
  • The seismic responses obtained from the inelastic seismic time history analyses with an accurate inelastic material model using Chaboche’s kinematic hardening model and the Voce isotropic hardening model for Type 316 stainless steel are in good agreement with those of the seismic shaking table tests.
  • From the comparison of the seismic strain time history responses at the nozzle between the tests and inelastic analyses, the location of the maximum strain responses from the inelastic analyses was found to be almost the same as the locations in the test results.
  • The structure damping value for the piping systems, recommended in the US NRC RG 1.61 [24], is 4% for the SSE level. However, in this study, the piping system made of Type 316 stainless steel revealed a much lighter damping value of about 0.57%.
  • For earthquakes large enough to result in inelastic behavior at the nuclear metal components, the inelastic seismic analysis is useful for reducing the seismic responses by energy dissipation due to hysteretic damping.
  • From the validation results obtained from the test mockup simulating the actual nuclear component in this study, it is assured that the inelastic seismic analysis method can be used for the seismic design of nuclear metal components in large earthquake scenarios such as the beyond-design-basis earthquake, which can cause significant plastic deformations.

Author Contributions

Conceptualization, G.-H.K., S.-W.A., J.-K.H. and J.-S.K.; methodology, G.-H.K., S.-W.A. and J.-S.K.; validation, G.-H.K.; formal analysis, G.-H.K.; investigation, S.-W.A. and J.-S.K.; data curation, S.-W.A.; writing—original draft preparation, G.-H.K.; writing—review and editing, S.-W.A., J.-K.H. and J.-S.K.; funding acquisition, G.-H.K., J.-K.H. and J.-S.K. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Ministry of Trade, Industry, and Energy through KETEP (Korea Institute of Energy Technology Evaluation Planning; No. 20181510102380).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This study was supported by the Ministry of Trade, Industry, and Energy through KETEP (Korea Institute of Energy Technology Evaluation Planning; No. 20181510102380).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Seismic Isolation Systems for Nuclear Installations; IAEA-TECDOC-1905; IAEA: Vienna, Austria, 2020.
  2. Koshida, H.; Yasaka, A. Vibration tests and earthquake observation results of base-isolated building. Seism. Shock. Vib. Isol. 1989, 181, 135–140. [Google Scholar]
  3. Koo, G.H.; Jung, J.Y.; Lee, J.H.; Shin, T.M. Development of Small-Sized Lead Inserted Laminated Rubber Bearing for Nuclear Component Seismic Isolation. Energies 2020, 13, 3193. [Google Scholar] [CrossRef]
  4. Kostarev, V.; Petrenko, A.; Vasilev, P. A New Method for Essential Reduction of Seismic and External Loads on NPP’s Structures, Systems and Components. In Proceedings of the 17th International Conference on Structural Mechanics in Reactor Technology, Prague, Czech Republic, 17–22 August 2003; p. 3216. [Google Scholar]
  5. Kwag, S.Y.; Kwag, J.S.; Lee, H.H.; Oh, J.H.; Koo, G.H. Enhancement in the Seismic Performance of a Nuclear Piping System using Multiple Tuned Mass Dampers. Energies 2019, 12, 2077. [Google Scholar] [CrossRef] [Green Version]
  6. Alternative Rules for Level D Service Limits of Class 1, 2, and 3 Piping Systems, Section III, Division 1, ASME Code Case (Record No. 13-1438); ASME: New York, NY, USA, 2019.
  7. Masaki, M.; Akihiro, O.; Tomoyoshi, W.; Izumi, N.; Tadahiro, S.; Masaki, S. Seismic Qualification of Piping Systems by Detailed Inelastic Response Analysis: Part 1—A Code Case for Piping Seismic Evaluation Based on Detailed Inelastic Response Analyses. In Proceedings of the ASME 2017 PVP Conference, Waikoloa, HI, USA, 16–20 July 2017; Volume 8. [Google Scholar]
  8. Akihiro, O.; Tadahiro, S.; Masaki, M.; Izumi, N.; Tomoyoshi, W.; Masaki, S. Seismic Qualification of Piping Systems by Detailed Inelastic Response Analysis: Part 2—A Guideline for Piping Seismic Inelastic Response Analysis. In Proceedings of the ASME 2017 PVP Conference, Waikoloa, HI, USA, 16–20 July 2017; Volume 8. [Google Scholar]
  9. Koo, G.H.; Kim, J.S.; Kim, Y.J. Feasibility Study on Strain-Based Seismic Design Criteria for Nuclear Components. Energies 2020, 13, 4435. [Google Scholar] [CrossRef]
  10. ASME. BPVC Section III, Rules for Construction of Nuclear Facility Components, Mandatory Appendix XIII, Design Based on Stress Analysis; ASME: New York, NY, USA, 2021. [Google Scholar]
  11. AFCEN. RCC-M: Design and Construction Rules for Mechanical Components of PWR Nuclear Islands; AFCEN: Lyon, France, 2017. [Google Scholar]
  12. Kim, J.Y.; Lee, J.M.; Park, J.G.; Kim, J.S.; Cho, M.K.; Ahn, S.W.; Koo, G.H.; Lee, B.H.; Huh, N.S.; Kim, Y.J.; et al. Round robin analysis to investigate sensitivity of analysis results to finite element elastic-plastic analysis variables for nuclear safety class 1 components under severe seismic load. Nucl. Eng. Technol. 2021. [Google Scholar] [CrossRef]
  13. Koo, G.H.; Kwag, S.; Nam, H.S. Study on Inelastic Strain-Based Seismic Fragility Analysis for Nuclear Metal Components. Energies 2021, 14, 3269. [Google Scholar] [CrossRef]
  14. Ravikiran, A.; Dubey, P.N.; Agrawal, M.K.; Reddy, G.R.; Singh, R.K.; Vaze, K.K. Experimental and Numerical Studies of Ratcheting in a Pressurized Piping System under Seismic Load. J. Press. Vessel. Technol. 2015, 137, 3. [Google Scholar] [CrossRef]
  15. Shirai, E.; Eto, K.; Umemoto, A.; Yoshii, T.; Kondo, M.; Shimizu, H.; Tai, K. Inelastic Seismic Test of the Small Bore Piping and Support System: Part 2—Static Failure Test for Piping Support Equipment. In Proceedings of the ASME 2008 PVP Conference, Chicago, IL, USA, 27–31 July 2008; Volume 8. [Google Scholar]
  16. ANSYS. ANSYS Mechanical APDL Release 15.0; ANSYS, Inc.: Canonsburg, PA, USA, 2014. [Google Scholar]
  17. Young, W.C.; Budynas, R.G.; Sadegh, A.M. Roark’s Formulas for Stress and Strain, 8th ed.; McGraw Hill: New York, NY, USA, 2012. [Google Scholar]
  18. Harris, C.M. Shock and Vibration Handbook; McGraw-Hill Book Company: New York, NY, USA, 1988. [Google Scholar]
  19. Chaboche, J.L.; Rousselier, G. On the plastic and viscoplastic constitutive equations—Part II: Application of internal variable concepts to the 316 stainless steel. J. Press. Vessel. Technol. 1983, 105, 159–164. [Google Scholar] [CrossRef]
  20. Chaboche, J.L. Constitutive equations for cyclic plasticity and cyclic viscoplasticity. Int. J. Plast. 1989, 5, 247–302. [Google Scholar] [CrossRef]
  21. Voce, E. A Practical Strain hardening Function. Metallurgia 1955, 51, 219–226. [Google Scholar]
  22. Koo, G.H.; Yoon, J.H. Inelastic Material Models of Type 316H for Elevated Temperature Design of Advanced High Temperature Reactors. Energies 2020, 13, 4548. [Google Scholar] [CrossRef]
  23. Standard Review Plan, 3.7.1 Seismic Design Parameters; NUREG-0800; U.S. Nuclear Regulatory Commission: Washington, DC, USA, 2012.
  24. Damping Values for Seismic Design of Nuclear Power Plants; Regulatory Guide 1.61; U.S. Nuclear Regulatory Commission: Washington, DC, USA, 2007.
Figure 1. Overall configuration of the test mockup and seismic directional axis.
Figure 1. Overall configuration of the test mockup and seismic directional axis.
Applsci 11 09264 g001
Figure 2. Target required response spectrum corresponding to PGA = 0.3 g.
Figure 2. Target required response spectrum corresponding to PGA = 0.3 g.
Applsci 11 09264 g002
Figure 3. Finite element analysis model: (a) total model; (b) zoomed nozzle model.
Figure 3. Finite element analysis model: (a) total model; (b) zoomed nozzle model.
Applsci 11 09264 g003
Figure 4. Mode shapes for test mockup: (a) Mode 1 (11.74 Hz); (b) Mode 2 (11.99 Hz); (c) Mode 3 (25.89 Hz); (d) Mode 4 (27.63 Hz); (e) Mode 5 (101.61 Hz).
Figure 4. Mode shapes for test mockup: (a) Mode 1 (11.74 Hz); (b) Mode 2 (11.99 Hz); (c) Mode 3 (25.89 Hz); (d) Mode 4 (27.63 Hz); (e) Mode 5 (101.61 Hz).
Applsci 11 09264 g004
Figure 5. Photo of test mockup with sensors on the shaking table.
Figure 5. Photo of test mockup with sensors on the shaking table.
Applsci 11 09264 g005
Figure 6. Location and pattern of the strain gauges attached on the test mockup: (a) photo of attached strain gauge; (b) strain gauge rosette pattern.
Figure 6. Location and pattern of the strain gauges attached on the test mockup: (a) photo of attached strain gauge; (b) strain gauge rosette pattern.
Applsci 11 09264 g006
Figure 7. Concept of principal strains from the triaxial 0°–45°–90° rosette gauge pattern.
Figure 7. Concept of principal strains from the triaxial 0°–45°–90° rosette gauge pattern.
Applsci 11 09264 g007
Figure 8. Test results of principal strain time history responses for PGA = 0.3 g: (a) SG-1; (b) SG-2; (c) SG-3; (d) SG-4; (e) SG-5.
Figure 8. Test results of principal strain time history responses for PGA = 0.3 g: (a) SG-1; (b) SG-2; (c) SG-3; (d) SG-4; (e) SG-5.
Applsci 11 09264 g008
Figure 9. Test results of principal strain time history responses for PGA = 0.6 g: (a) SG-1; (b) SG-2; (c) SG-3; (d) SG-4; (e) SG-5.
Figure 9. Test results of principal strain time history responses for PGA = 0.6 g: (a) SG-1; (b) SG-2; (c) SG-3; (d) SG-4; (e) SG-5.
Applsci 11 09264 g009
Figure 10. Test results of acceleration responses at the shaking table for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Figure 10. Test results of acceleration responses at the shaking table for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Applsci 11 09264 g010
Figure 11. Test results of acceleration responses at Accelerometer-2 for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Figure 11. Test results of acceleration responses at Accelerometer-2 for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Applsci 11 09264 g011
Figure 12. Test results of acceleration responses at Accelerometer-3 for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Figure 12. Test results of acceleration responses at Accelerometer-3 for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Applsci 11 09264 g012
Figure 13. Resonance test results for each direction: (a) EW (X); (b) NS (Z); (c) vertical (Y).
Figure 13. Resonance test results for each direction: (a) EW (X); (b) NS (Z); (c) vertical (Y).
Applsci 11 09264 g013
Figure 14. Spectrum illustrating the half-power bandwidth method for damping estimation (1st EW mode).
Figure 14. Spectrum illustrating the half-power bandwidth method for damping estimation (1st EW mode).
Applsci 11 09264 g014
Figure 15. Displacement input motions used for analyses for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Figure 15. Displacement input motions used for analyses for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Applsci 11 09264 g015
Figure 16. Test response spectrum of shaking table motions.
Figure 16. Test response spectrum of shaking table motions.
Applsci 11 09264 g016
Figure 17. Inelastic seismic analysis results of accumulated equivalent plastic strain distributions at the nozzle (PGA = 0.3 g): (a) isometric view; (b) sectional view at the maximum region and node points corresponding to locations of strain gauges attached in the tests.
Figure 17. Inelastic seismic analysis results of accumulated equivalent plastic strain distributions at the nozzle (PGA = 0.3 g): (a) isometric view; (b) sectional view at the maximum region and node points corresponding to locations of strain gauges attached in the tests.
Applsci 11 09264 g017
Figure 18. Hysteretic behavior of the material at the location of the maximum accumulated equivalent plastic strain for PGA = 0.6 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Figure 18. Hysteretic behavior of the material at the location of the maximum accumulated equivalent plastic strain for PGA = 0.6 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Applsci 11 09264 g018
Figure 19. Inelastic seismic analysis results of principal strain time history responses for PGA = 0.3 g: (a) SG-1; (b) SG-2; (c) SG-3; (d) SG-4; (e) SG-5.
Figure 19. Inelastic seismic analysis results of principal strain time history responses for PGA = 0.3 g: (a) SG-1; (b) SG-2; (c) SG-3; (d) SG-4; (e) SG-5.
Applsci 11 09264 g019
Figure 20. Inelastic seismic analysis results of principal strain time history responses for PGA = 0.6 g: (a) SG-1; (b) SG-2; (c) SG-3; (d) SG-4; (e) SG-5.
Figure 20. Inelastic seismic analysis results of principal strain time history responses for PGA = 0.6 g: (a) SG-1; (b) SG-2; (c) SG-3; (d) SG-4; (e) SG-5.
Applsci 11 09264 g020aApplsci 11 09264 g020b
Figure 21. Inelastic seismic analysis results of acceleration responses at Accelerometer-2 for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Figure 21. Inelastic seismic analysis results of acceleration responses at Accelerometer-2 for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Applsci 11 09264 g021
Figure 22. Inelastic seismic analysis results of acceleration responses at Accelerometer-3 for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Figure 22. Inelastic seismic analysis results of acceleration responses at Accelerometer-3 for PGA = 0.3 g: (a) EW (X); (b) NS (Z); (c) V (Y).
Applsci 11 09264 g022
Figure 23. Comparison results of the response spectrum at the pipe end (Acceleromter-3) for PGA = 0.3 g.
Figure 23. Comparison results of the response spectrum at the pipe end (Acceleromter-3) for PGA = 0.3 g.
Applsci 11 09264 g023
Table 1. Summary of test mockup design parameters.
Table 1. Summary of test mockup design parameters.
Design ParametersDesign Value
Pipe Outer Diameter (mm)89.1
Pipe Thickness (mm)7.6
Vertical Pipe Length (mm)850
Horizontal Pipe Length (mm)550
Elbow Curvature (mm)114.3
Added Mass-1 (kg)80
Added Mass-2 (kg)40
Total Height of Nozzle (mm)167.35
Nozzle Outer Diameter (mm)133.34
Table 2. Summaries of modal analysis.
Table 2. Summaries of modal analysis.
Mode No.Frequencies
(Hz)
Modal Participation FactorsEffective Mass
EWNSVEWNSV
111.740.0010.070.000.00101.460.00
211.99−8.740.004.3576.300.0018.96
325.890.005.470.000.0029.940.00
427.637.740.005.0559.840.0025.54
5101.610.000.070.000.005.810.00
Table 3. Shaking table performance.
Table 3. Shaking table performance.
ItemPerformance
Max. Loading (kg)2000
Table Size (mm)2209 × 2173
Excitation Axes6 DOF (Translation 3 axes, Rotational 3 Axes)
Max. Displacement (mm)Hor. (X, Z) = (±125, ±140), Ver. (Y) = ±110
Max. Accel. for bare table (g)Hor. (X, Z) = 17, Ver. (Y) = 21
Max. Accel. with max load (g)Hor. (X, Z) = 4.9, Ver. (Y) = 6.0
Frequency Range (Hz)0.8~100
Excitation MechanismElectro-hydraulic Servo, 3 Variable Control
Control Software354.20/MTS
Table 4. Specifications of the used strain gauge.
Table 4. Specifications of the used strain gauge.
ItemSpecifications
Model NameKFGS-1-120-D17-11L15M3S (KYOWA)
Gauge Length (mm)1
Gauge Resistance (Ω)120 ± 0.7%
Gauge PatternTriaxial 0°–45°–90° Rosette, Round base
Adoptable Thermal Expansion11.7 × 10−6/°C
Applicable AdhesiveCC-33A, EP-340
Table 5. Specifications of the used accelerometers.
Table 5. Specifications of the used accelerometers.
ItemSpecifications
Model NameType 8396A (KISTLER)
Sensing TypeMEMS Variable Capacitance, Silicon Sensing
Element
Measuring Range (g)50
Measuring AxisTriaxial
Measuring Freq. Range (Hz)0.5~5000
Operating Temperature (°C)−55~125
Table 6. Coordinates of the attached strain gauges.
Table 6. Coordinates of the attached strain gauges.
Strain Gauge IDX (EW)
(mm)
Y (V)
(mm)
Z (NS)
(mm)
Corresponding
Node (1)
SG-126.2156.036.019118
SG-240.7156.018.119125
SG-344.5156.00.0137
SG-440.7156.0−18.1851
SG-526.2156.0−36.0858
The origin global coordinates (X = 0, Y = 0, Z = 0) are at the center of the nozzle bottom surface. (1) The detailed node locations are described in Figure 17(b) below.
Table 7. Test results of ZPA responses.
Table 7. Test results of ZPA responses.
DirectionsPGATarget Input
Motions
Shaking Table
(Accelerometer-1)
End of Pipe
(Accelerometer-3)
ZPA (g)ZPA Ratio *ZPA (g)ZPA RatioZPA (g)ZPA Ratio
EW (X)0.3 g1.182.01.711.955.551.57
0.6 g2.353.338.76
NS (Z)0.3 g0.862.01.472.129.261.37
0.6 g1.753.1312.70
V (Y)0.3 g1.332.01.861.956.231.58
0.6 g2.653.629.88
* The ZPA ratio represents the ratio of ZPA responses corresponding to PGA = 0.6 g and PGA = 0.3 g.
Table 8. Comparison of resonance frequencies between tests and analyses.
Table 8. Comparison of resonance frequencies between tests and analyses.
DirectionsMode 1 (Hz)Mode 2 (Hz)
EW (X)Test11.7528.75
Analysis11.9927.63
NS (Z)Test11.7527.38
Analysis11.7425.89
V (Y)Test11.5628.56
Analysis11.9927.63
Table 9. Material constants for inelastic material models.
Table 9. Material constants for inelastic material models.
Materialσyo
× 106 (Pa)
E
× 109
(Pa)
C1
× 109
C2
× 109
C3
× 109
γ1
× 103
γ2
× 103
γ3bQ
× 106
Type 316SS13519012020.210.671.01.01.045.085
Table 10. Comparison results of principal strain responses.
Table 10. Comparison results of principal strain responses.
Strain Gauge IDTests (%)Inelastic Seismic Analysis (%)
PGA = 0.3 gPGA = 0.6 gPGA = 0.3 gPGA = 0.6 g
Max(εmax)Min(εmin)Max(εmax)Min(εmin)Max(εmax)Min(εmin)Max(εmax)Min(εmin)
SG-10.139−0.1440.185−0.2020.101−0.1390.184−0.247
SG-20.125−0.1190.212−0.2110.123−0.1200.199−0.197
SG-30.139−0.1580.249−0.2620.119−0.0990.255−0.189
SG-40.148−0.1610.217−0.2390.127-0.1280.273−0.194
SG-50.152−0.1680.228−0.2300.165−0.1670.290−0.216
Table 11. Comparison results of ZPA responses at the pipe end (Accelerometer-3).
Table 11. Comparison results of ZPA responses at the pipe end (Accelerometer-3).
DirectionsPGA
(g)
Tests
(g)
Inelastic Seismic
Analysis (g)
EW0.35.556.17
0.68.769.90
NS0.39.2612.68
0.612.7016.80
V0.36.238.71
0.69.8812.60
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Koo, G.-H.; Ahn, S.-W.; Hwang, J.-K.; Kim, J.-S. Shaking Table Tests to Validate Inelastic Seismic Analysis Method Applicable to Nuclear Metal Components. Appl. Sci. 2021, 11, 9264. https://doi.org/10.3390/app11199264

AMA Style

Koo G-H, Ahn S-W, Hwang J-K, Kim J-S. Shaking Table Tests to Validate Inelastic Seismic Analysis Method Applicable to Nuclear Metal Components. Applied Sciences. 2021; 11(19):9264. https://doi.org/10.3390/app11199264

Chicago/Turabian Style

Koo, Gyeong-Hoi, Sang-Won Ahn, Jong-Keun Hwang, and Jong-Sung Kim. 2021. "Shaking Table Tests to Validate Inelastic Seismic Analysis Method Applicable to Nuclear Metal Components" Applied Sciences 11, no. 19: 9264. https://doi.org/10.3390/app11199264

APA Style

Koo, G. -H., Ahn, S. -W., Hwang, J. -K., & Kim, J. -S. (2021). Shaking Table Tests to Validate Inelastic Seismic Analysis Method Applicable to Nuclear Metal Components. Applied Sciences, 11(19), 9264. https://doi.org/10.3390/app11199264

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop