1. Introduction
The study of earthquake ground motions and associated earthquake risks plays an important role in the sustainable development of earthquake-prone countries. Site response analyses are usually performed to analyze ground amplification and to develop the design response spectra [1].
The northern part of Vietnam is located in the plate boundary between the South China plate and the Indochina plate. The plate boundary is the Red River fault zone that strikes in the NW–SE direction and extends over a length of 1000 km between the Gulf of Tonkin and Tibet [2]. It caused large scale displacement and major tectonic features in Vietnam, such as rifting and folding. Hanoi is located to the North and the capital of Vietnam lies in the center of the triangular basin of the Red River. The center of the city is placed on the south bank of the river. Most of the deltaic land lies on both sides of the Red River and its tributaries. Hanoi is also affected by many adjacent zones such as the Lai Chau–Dien Bien (LC–DB), Son La, Ma River, Lo River, and Chay River faults (Figure 1) [3].
A total of 1645 earthquakes with a local magnitude (ML) of 3.0 and higher were recorded in Vietnam from 114 to 2003 [4]. There were three earthquakes recorded at Hanoi in 1276, 1278, and 1285 with ML=5.5–6.0. In 1983, several moderate earthquakes happened in the northern provinces of Vietnam. However, even with a distance of approximately 300 km from these epicenters, the people of Hanoi experienced a strong vibration and a number of buildings were damaged. It revealed that amplification of seismic ground motion through soft soil caused a high amplification. Recently, earthquakes occurred in the neighboring provinces of Hanoi such as Dien Bien in 2001, Ha Giang in 2005, Ninh Binh in 2005, and Hoa Binh in 2010 with a ML of 5.3, 4.7, 3.1, and 4.0, respectively.
Site-specific ground response analysis is a usual practice to predict the local site condition effects. Shylamoni et al. [5] performed the site-specific one-dimensional (1-D) equivalent linear (EQL) site response analysis for the nuclear power plants in Japan. Grasso and Maugeri [6] carried out the seismic response analysis on the Italian site conditions to develop a seismic microzoning map for Catania City. Naik and Choudhury [7] conducted the 1-D EQL site response analysis for Goa City in India using the data of the Bhuj earthquake 2001. Recently, Aaqib et al. [8] performed 1-D nonlinear site response analyses for the Korean sites to investigate the effect of strength correction on the seismic site coefficients. Mahmood et al. [9] implemented the 1-D EQL site response analysis for sites in Islamabad, Pakistan.
Due to the recent earthquakes as well as the urban infrastructure development, the impact of the earthquake in Hanoi has also been considered and studied by performing site-specific seismic response analyses [10,11,12,13]. Nhung et al. [10] and Alexandr and Chi [11] assessed the hazard and possibility of liquefaction of soil layers caused by earthquakes in the urban area of the city. Wen et al. [13] concluded that the soil deposit in the Hanoi region may increase earthquake damage because of the amplification of ground motion after performing site response analysis from the microtremor survey. The current seismic design code of Vietnam, TCVN 9386:2012 [14], was constructed based on Eurocode 8 (EC-8) [15] with consideration of the Vietnamese seismic map. The above-mentioned studies revealed that there is limited literature available on the site-specific response analysis in Vietnam. The relative difference between the design spectra based on the site-specific amplification factors and the code-based spectra has not been investigated yet. Moreover, all the aforesaid studies did not consider a range of input motions to account for the variability in frequency contents of ground motions.
The purpose of this study is to perform site-specific nonlinear (NL) and EQL response analyses based on the profiles from different locations of Central Hanoi City (Vietnam), using a range of input ground motions. The obtained spectra shapes are compared with those of the Vietnamese seismic design code, TCVN 9386:2012, and differences between NL and EQL methods are thoroughly discussed. 2. An Insight of the Vietnamese Seismic Design Code TCVN 9386: 2012
The TCVN 9386: 2012 was developed based on EC-8. However, there were some significant changes [4] to conform with the local Vietnamese conditions. Firstly, the seismic map has been developed from the data obtained from Vietnam Institute of Geophysics (VIG) [16] as shown in Figure 2. Secondly, the importance factors and importance classes in TCVN 9386: 2012 are slightly different from those of EC-8. In TCVN 9386: 2012, there were five importance classes with importance factors ranged from 0.75 to the value corresponding to the possible maximum acceleration. Note that there is no required earthquake resistance for the final importance class (i.e., class IV) in the Vietnamese design code, whereas the EC-8 is classified into four importance classes with importance factors ranged from 0.8 to 1.4. The importance classes and importance factors specified by EC-8 and TCVN 9386: 2012 are listed in Table 1 and Table 2, respectively. Thirdly, as the strongest earthquake regions have a magnitude greater than 5.5 on the local magnitude scale and lack data to establish the spectral curves, the TCVN 9386: 2012 used only the type 1 elastic response spectra in two types of EC-8 as the design response spectrum for all regions in Vietnam. Figure 3a shows the detailed type 1 elastic response spectra of EC-8 or the design response spectrum in TCVN 9386: 2012.
The horizontal elastic response spectrumSe(T)is determined by the following expressions:
0≤T≤TB: Se(T)=ag⋅S⋅[1+TTB⋅(η⋅2.5−1)]
TB≤T≤TC: Se(T)=ag⋅S⋅η⋅2.5
TC≤T≤TD: Se(T)=ag⋅S⋅η⋅2.5⋅[TCT]
TB≤T≤TC: Se(T)=ag⋅S⋅η⋅2.5
where:
-
Se(T): is the elastic response spectrum.
-
T: is the vibration period of a linear single-degree-of-freedom system.
-
ag: is the design ground acceleration on type A ground (ag=γI⋅agR).
-
TB: is the lower limit of the period of the constant spectral acceleration branch.
-
TC: is the upper limit of the period of the constant spectral acceleration branch.
-
TD: is the value defining the beginning of the constant displacement response range of the spectrum.
-
S: is the soil factor.
-
η: is the damping correction factor with a reference value of η = 1 for viscous damping ξ = 5%
η=10/(5+ξ)≥0.55 in TCVN 9386: 2012
η=10/(5+ξ)≥0.55 in EC−8
For the five soil types, A, B, C, D, and E, the values of parameters S, TB, TC, and TD are given in Table 3, and normalized shapes by ag for 5% damping are shown in Figure 3b.
3. Site Response Analysis Methodology
DEEPSOIL v7.0 [17] was used to perform equivalent linear and nonlinear 1-D site response analyses. The most widely used pressure-dependent hyperbolic, the Modified Kodner–Zelasko (MKZ) model [18], implemented in DEEPSOIL, was used. The dynamic curves proposed by Darendeli [19] were adopted. The variables required to generate the nonlinear curves for each layer are the coefficient of lateral earth pressure (K0), plasticity index (PI), number and frequency of cycles (N), loading frequency (f), and over consolidation ratio (OCR). Because of the unavailability of site-specific index properties, PI for the sand and clay was assumed as 0% and 15% [20,21,22,23,24,25], respectively, whereas OCR was taken as 1. N and f were set to 10 and 1, respectively, as recommended by Darendeli [19], Ding et al. [22], and Harmon et al. [23]. K0 was calculated as 0.5 using the Jacky [26] equation based on the representative friction angle of 30 degrees, as specified in the laboratory investigation report. Small strain damping (Dmin) was also modeled using the functions of Darendeli [19], which predicts generally larger damping values near the ground surface, whereas it keeps on decreasing along with the depth. Figure 4 shows the site response analysis flowchart in this study. A detailed discussion about soil profiles, input ground motions is presented in sections below.
3.1. Soil Profiles
A total of three districts in Hanoi were selected for the analyses. Two soil profiles were selected from Cau Giay District (District 1), one from Hoang Mai District (District 2), and three soil profiles were selected from Hai Ba Trung District (District 3). The typical soil stratigraphy includes fill, sand, clay, and gravel. Due to the lack of shear wave velocity (VS) profile measurements, the in situ standard penetration test blow count (SPT-N) measurements were used to evaluate the VS profiles. The correlation between SPT-N and VS proposed by Marto et al. [27] was used to determine the VS profiles. This correlation was also employed by Trung et al. [28] for Vietnamese soil conditions, expressed as:
Vs=77.13N0.377
The depth to bedrock of the profiles range from 33 m to 75 m. Natural periods of the sites range from 0.64 to 1.39 s. The time-average shear wave velocity to 30 m depth (VS30) of the sites range from 160 m/s to 202 m/s. Sites are classified as types C and D according to the Vietnamese design code. The bedrock is considered at the point where the VS reached to 800 m/s, according to the bedrock definition in TCVN 9386: 2012. The detail of soil profiles is shown in Table 4 and Figure 5.
3.2. Input Ground Motions
To ensure an appropriate representation of ground motion intensities for the districts in Hanoi, earthquake motions representing the seismicity of the region were selected. Due to the lack of earthquake records available in Hanoi, ground motions recorded in other regions with a range of magnitude similar to the ancient earthquakes in Vietnam were selected. The ground motions recorded from earthquakes having a variety of rupture mechanisms were selected to account for the effect of the fault mechanism. The recorded ground motions were selected from the NGA-west2 database (https://ngawest2.berkeley.edu/). The ground motions with the following features were selected:
- A range of moment magnitudes (Mw) 5.5–7.0.
- Rupture distance (Rrup) 0.92–85.17 km.
- VS30 of 760–1500 m/s.
A total of 15 motions with peak ground acceleration (PGA) ranging from 0.04 g to 0.25 g were selected from six earthquake events, as shown in Table 5. The ground motions were scaled to the representative shaking intensity levels of each analyzed district, in which PGA levels at Cau Giay, Hoang Mai, and Hai Ba Trung Districts are scaled to 0.1032 g, 0.1001 g, and 0.0959 g, respectively. Furthermore, the mean response spectrum of the input ground motions matched closely to the target spectrum of site class A of the Vietnamese design code [14], as shown in Figure 6.
4. Spectral Shapes The soil nonlinearity can be characterized either through the equivalent linear or the fully nonlinear approach. In this section, the normalized spectral shapes from both the NL and EQL approaches are compared with the normalized design spectrum from TCVN 9386: 2012 to evaluate the applicability of the design code.
Figure 7 and Figure 8 compare the spectral shapes from the NL and EQL analyses with the design spectrum for the profiles of both site classes C and D. Districts 1 and 3 consist of site class C profiles. The mean normalized spectrum calculated in this study is higher than the normalized design spectrum for the period range from 0.6 to 2 s in the case of both the NL and EQL approaches. District 2 and District 3 consist of site class D profiles. The normalized spectral shape is lower than the design spectrum across the whole period range in the case of the NL method, whereas it matches the design spectrum well in the case of the EQL method. The underestimation of calculated spectral shape in the case of the NL approach can be attributed to the high levels of nonlinearity induced because of the softness of site class D profiles, which could lead to a certain amount of deamplification. It should be noted that most building periods are greater than 0.3 s [29].
Figure 9 shows a comparison of the NL and EQL analysis response spectra and the design spectrum. As can be seen from the figure, the mean normalized response spectrum in the case of the EQL analysis is higher than that of the NL analysis in both the site classes, especially near the natural site periods. The NL analysis predicts smaller spectral accelerations from 0.2 s to 1.6 s in case of profiles of site class C, whereas the period range extends up to 2.4 s in case of site class D profiles. This difference is attributed to the persistent softening of properties used in the nonlinear analysis which is not the case in the EQL analysis, where constant and strain compatible properties are used. These results are in line with the study of Rathje and Kottke [30].
5. Maximum Shear Strain and PGA 5.1. Maximum Shear Strain
Figure 10a–f show the maximum shear strain results obtained from the NL and EQL analysis for soil profiles. Overall, the NL and EQL results are closely matched. The difference between the EQL and NL increases as the higher strains are approached, which is in line with the study of Kim et al. [31]. The maximum shear strain of soil profiles calculated in this study is less than 0.5%, excluding soil profile P6 (site class D) which reaches up to 1% and 2.5% for the NL and EQL approaches, respectively. In each soil profile, both the NL and EQL approaches share the same depth where the maximum shear strain is encountered, except soil profile P3. The difference in maximum shear strain depths between the two approaches in profile P3 can be attributed to the presence of a significant reversal of shear wave velocity in between 31 m to 40 m.
5.2. Maximum PGA
Figure 11a–f show the NL and EQL maximum PGA for soil profiles. The difference in maximum PGA of NL and EQL at the surface is significant. The difference can be reached by 60%, i.e., in the case of soil profile P6 which can be attributed to the softness of the profile leading to an increase in the shear strains, inducing more nonlinearity and resulting in the PGA difference at the surface. In the case of the NL approach, a certain amount of deamplification at the surface is observed in the case of soil profile P3 and soil profile P5. This can be attributed to the shear wave reversals present throughout the mentioned soil profiles. However, in the case of the EQL approach, the amount of deamplification is negligible as the surface is approached. This can be demonstrated by the fact that shear wave reversals result in the increased nonlinearity and a certain amount of deamplification occurs while performing the NL approach which is not the case in the EQL approach.
6. Conclusions This study investigates the applicability of the Vietnamese seismic design code TCVN 9386:2012 using one-dimensional equivalent linear (EQL) and nonlinear (NL) site response approaches and discusses the relative differences between the two methods. A total of six soil profiles from three districts in Hanoi (Vietnam) and fifteen ground motions compatible with the rock design spectrum of the Vietnam design code were used as input data to conduct 1-D site response analyses. The following conclusions are achieved based on numerical analyses. (1) The TCVN 9386:2012 underestimates the response amplification for soil profiles falling to site class C in the range of period from 0.6 to 2 s and it overestimates the response across the whole period range for site class D (except in the case of the EQL method). It is therefore suggested that the seismic design code of Vietnam needs to be modified for both stiff profiles (i.e., site class C) and soft profiles (i.e., site class D).
(2) It is well-known that the NL approach is appropriate to capture the actual response rather than the EQL approach, which is recommended to be used only in the case of shear strain level of less than 0.4% [31,32]. In this study, the spectra shapes from the NL analyses are higher than those of the EQL analyses, especially in soft soil (i.e., site class D). Overall, the shear strain and PGA obtained from the EQL method can adequately approximate the NL below the threshold shear strain levels. The difference between the NL and EQL results are greater for soft soil. For the PGA, the difference occurs at the ground surface, the PGA of the EQL analysis is higher than that of the NL. Based on the results of the NL analyses, the maximum shear strain is less than 0.5%, except soil profile P6 with 1%, and PGA is lower than 0.16 g.
(3) The NL results revealed the occurrence of deamplification in the soft soils of Hanoi. It is recommended that the site-specific response analysis should be considered to perform the seismic design and response of above structures in site class D.
Importance Class | Buildings | Recommended Importance Factor, γI |
---|---|---|
I | Building of minor importance for public safety, e.g., agricultural building, etc. | 0.8 |
II | Ordinary building, not belonging in the other categories | 1.0 |
III | Building whose seismic resistance is of importance in view of the consequences associated with a collapse, e.g., schools, assembly halls, cultural institutions, etc. | 1.2 |
IV | Building whose integrity during earthquakes is of vital importance of civil protection, e.g., hospitals, fire stations, power plants, etc. | 1.4 |
Importance Class | Buildings | Recommended Importance Factor, γI |
---|---|---|
Special | The building has special importance, not allowing any damage by the earthquake. | Designed with maximum possible acceleration |
I | The building has vital importance for the protection of the community, its function must not be interrupted during the earthquake. | 1.25 |
II | The building has an importance in preventing earthquake consequences, if collapsed causing great loss of life and property. | 1.0 |
III | The building does not belong to importance classes: special, I, II, and IV. | |
IV | The building has minor importance for the safety of human life. | 0.75 |
Ground Type | S | TB (s) | TC (s) | TD (s) |
---|---|---|---|---|
A | 1.0 | 0.15 | 0.4 | 2.0 |
B | 1.2 | 0.15 | 0.5 | 2.0 |
C | 1.15 | 0.2 | 0.8 | 2.0 |
D | 1.35 | 0.2 | 0.8 | 2.0 |
E | 1.4 | 0.15 | 0.5 | 2.0 |
Soil Profile | Soil Class (According to EC-8 and TCVN 9386: 2012) | Location | Notation |
---|---|---|---|
P1 | C | Cau Giay District (District 1) | 1-C |
P2 | C | Cau Giay District | 1-C |
P3 | D | Hoang Mai District (District 2) | 2-D |
P4 | C | Hai Ba Trung District (District 3) | 3-C |
P5 | D | Hai Ba Trung District | 3-D |
P6 | D | Hai Ba Trung District | 3-D |
No | Earthquake Name | Year | Station | Mag. (MW) | Mechanism | Rrup (Km) | VS30 (m/s) |
---|---|---|---|---|---|---|---|
1 | San Fernando | 1971 | Pasadena—Old Seismo Lab | 6.61 | Reverse | 21.5 | 969.07 |
2 | Whittier Narrows-01 | 1987 | Pasadena—CIT Kresge Lab | 5.99 | Reverse | 18.12 | 969.07 |
3 | Loma Prieta | 1989 | Piedmont Jr High School Grounds | 6.93 | Reverse Oblique | 73 | 895.36 |
4 | Loma Prieta | 1989 | Point Bonita | 6.93 | Reverse Oblique | 83.45 | 1315.92 |
5 | Loma Prieta | 1989 | SF—Pacific Heights | 6.93 | Reverse Oblique | 75.96 | 1249.86 |
6 | Loma Prieta | 1989 | So. San Francisco_Sierra Pt. | 6.93 | Reverse Oblique | 63.15 | 1020.62 |
7 | Northridge-01 | 1994 | LA—Wonderland Ave | 6.99 | Reverse | 20.29 | 1222.52 |
8 | Northridge-01 | 1994 | Vasquez Rocks Park | 6.99 | Reverse | 23.64 | 996.43 |
9 | Chi-Chi_Taiwan-05 | 1999 | TTN042 | 6.20 | Reverse | 85.17 | 845.34 |
10 | Umbria-03_Italy | 1984 | Gubbio | 5.60 | Normal | 15.72 | 922 |
11 | Kobe_Japan | 1995 | Kobe University | 6.90 | Strike slip | 0.92 | 1043 |
12 | Chi-Chi_Taiwan-05 | 1999 | HWA002 | 6.20 | Reverse | 45.03 | 789.18 |
Author Contributions
Conceptualization and methodology, M.A. and D.P.; formal analysis, resources, and writing-original draft, V.-Q.N.; review and editing, D.-D.N. and N.-V.L. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by a National Research Foundation of Korea (NRF) grant funded by the Korean government (MEST), grant number NRF-2019R1A2C1011323.
Acknowledgments
The authors would like to express their thanks to all participants.
Conflicts of Interest
The authors declare no conflicts of interest.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2020. This work is licensed under http://creativecommons.org/licenses/by/3.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
A series of one-dimensional (1-D) site response analyses were performed using the nonlinear (NL) and equivalent linear (EQL) approaches to assess the applicability of the Vietnamese earthquake-resistance design code TCVN 9386: 2012. Six soil profiles were selected from three districts in Hanoi (Vietnam). A number of ground motions compatible with the rock design spectrum were used as input for carrying out analyses. The results highlight that the calculated response is higher than the design spectrum for site class C and lower for site class D. The normalized response spectra of the EQL approach results are higher than those of the NL approach. Moreover, the peak ground accelerations at the surface from EQL analyses are greater than those of the NL method because the latter generates a higher amount of nonlinearity. The results from the NL approach also illustrate that the deamplification phenomenon occurs in the soft soils of the Hanoi region (e.g., soil profile P3 and P5 of site class D). Additionally, the shear strains calculated from the NL method are closely matched with those from the EQL method, the difference between them increasing with a decrease in soil stiffness.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer