Numerical Methods in Civil Engineering

Numerical Methods in Civil Engineering

Evaluation of ground motion amplification across empty valleys subjected to inclined SV waves using boundary element method

Document Type : Research

Authors
Ehsan Nafici 1
Behrouz Gatmiri 2
1 Graduated Student in Geotechnical Engineering, Department of Civil Engineering, College of Engineering, University of Tehran, 16 Azar St., Enghelab St., Tehran, Iran,
2 Full Professor, Department of Civil Engineering, College of Engineering, University of Tehran,
10.61186/NMCE.2303.1031
Abstract
Many recorded earthquakes in last recent decades have shown that the damage distribution and seismic response of an area depends on Source characteristics, pathway of seismic waves and also on local conditions of the site; these conditions can change the amplitude of displacement and frequency properties of ground movement. Local site effects include geotechnical and topographical effects. Unfortunately, most modern seismic codes consider geotechnical effects, and those who consider topographical effects, focus on simple geometries but other parameters such as wave characteristics are not taken into account; the main reason is lack of quantitative amplification predicting methods and practical information in topography constructions.2D seismic site effects in empty valleys are evaluated in present study. Numerical analysis is carried out by HYBRID software. HYBRID is a nonlinear, two phase, code for solving 2D problems of propagation of the waves; this software combines F.E.M in near field and B.E.M in far field. Empty valleys are modeled, with different shapes including rectangle, triangle and trapezoid and with depth ratios equal to 1, 0.8, 0.6, 0.4 and 0.2. Oblique SV waves of Ricker type are considered as seismic solicitation and its angle of incidence varies from -90º to 90º every 10º. It is shown that the incidence angle is able to change the maximum amplification and critical incidence angle of a half-space is changed by topography.
Keywords
topographical effects
HYBRID
2D site effect
oblique SV waves
empty valleys

Subjects

[1] Pagliaroli, A., Pergalani, F., Ciancimino, A., Chiaradonna, A., Compagnoni, M. A., De Silva, F., Foti, S., Giallini, S., Lanzo, G., Lombardi, F., Luzi, L., Macerola, L., Nocentini, M., Pizzi, A., Tallini, M., & Teramo, C. (2019). Site response analyses for complex geological and morphological conditions: relevant case-histories from 3rd level seismic microzonation in Central Italy. Bulletin of Earthquake Engineering, 18(12), 5741–5777.
[2] Çetin, K. Ö., Papadimitriou, A. G., Altun, S., Pelekis, P., Unutmaz, B., Rovithis, E., Akgün, M., Klimis, N., Gündoğan, A. A., Ziotopoulou, K., Sezer, A., Kıncal, C., Ilgaç, M., Can, G., Çakır, E., Söylemez, B., Al-Suhaily, A., Elsaid, A., Zarzour, M., . . . Mylonakis, G. (2021). The role of site effects on elevated seismic demands and corollary structural damage during the October 30, 2020, M7.0 Samos Island (Aegean Sea) Earthquake. Bulletin of Earthquake Engineering, 20(14), 7763–7792.
[3] Hartzell, S. H., Carver, D., & King, K. W. (1994). Initial investigation of site and topographic effects at Robinwood Ridge, California. Bulletin of the Seismological Society of America, 84(5), 1336–1349.
[4] Bouchon, M., & Barker, J. S. (1996). Seismic response of a hill: The example of Tarzana, California. Bulletin of the Seismological Society of America, 86(1A), 66–72.
[5] Spudich, P., Hellweg, M., & Lee, W. H. K. (1996). Directional topographic site response at Tarzana observed in aftershocks of the 1994 Northridge, California, earthquake: Implications for mainshock motions. Bulletin of the Seismological Society of America, 86(1B), S193–S208.
[6] Athanasopoulos, G., Pelekis, P., & Leonidou, E. (1999). Effects of surface topography on seismic ground response in the Egion (Greece) 15 June 1995 earthquake. Soil Dynamics and Earthquake Engineering, 18(2), 135–149.
[7] Bouckovalas, GD, Kouretzis, G. Review of soil and topography effects in the September 7, 1999 Athens (Greece) earthquake. In: Proceedings of the fourth international conference on recent advances in geotechnical earthquake engineering and soil dynamics and symposium in honor of Professor WD Liam Finn. San Diego, California; 2001.
[8] Trifunac, M. D., & Hudson, D. E. (1971). Analysis of the Pacoima dam accelerogram—San Fernando, California, earthquake of 1971. Bulletin of the Seismological Society of America, 61(5), 1393–1411.
[9] Assimaki, D. (2005). Effects of Local Soil Conditions on the Topographic Aggravation of Seismic Motion: Parametric Investigation and Recorded Field Evidence from the 1999 Athens Earthquake. Bulletin of the Seismological Society of America, 95(3), 1059–1089.
[10] Géli, L., Bard, P., & Jullien, B. (1988). The effect of topography on earthquake ground motion: A review and new results. Bulletin of the Seismological Society of America, 78(1), 42–63.
[11] Ameri, G., Massa, M., Bindi, D., D’Alema, E., Gorini, A., Luzi, L., Marzorati, S., Pacor, F., Paolucci, R., Puglia, R., & Smerzini, C. (2009). The 6 April 2009 MW 6.3 L’Aquila (Central Italy) earthquake: strong-motion observations. Seismological Research Letters, 80(6), 951–966.
[12] Biondi, G., & Maugeri, M. (2005). Seismic response analysis of Monte Po hill (Catania). In WIT transactions on state-of-the-art in science and engineering.
[13] Fukushima, Y., Irikura, K., Uetake, T., & Hisashi, M. (2000). Characteristics of Observed Peak Amplitude for Strong Ground Motion from the 1995 Hyogoken Nanbu (Kobe) Earthquake. Bulletin of the Seismological Society of America, 90(3), 545–565.
[14] HB.Ozmen, An Investigation on Soil Amplification through Site Factors Used in Seismic Design Codes, Advances in Civil Engineering Volume 2023, Article ID 6858371, 11 pages
[15] Bard, P. (1997). Local effects of strong ground motion: Basic physical phenomena and estimation methods for microzoning studies.
[16] Lacave, C & Bard, Pierre-Yves & Koller, Martin. (1999). Microzonation: Techniques and examples. Block 15: Naturgefahren-Erdbebenrisiko
[17] Sánchez-Sesma, F. J., Palencia, V. J., & Luzón, F. (2004). Estimation of local site effects during earthquakes: An overview. From Seismic Source to Structural Response: Contributions of Professor Mihailo D. Trifunac, 44-70.
[18] Giallini, S., Pizzi, A., Pagliaroli, A., Moscatelli, M., Vignaroli, G., Sirianni, P., Mancini, M., & Laurenzano, G. (2020). Evaluation of complex site effects through experimental methods and numerical modelling: The case history of Arquata del Tronto, central Italy. Engineering Geology, 272, 105646. https://doi.org/10.1016/j.enggeo.2020.105646
[19] Borcherdt, R. D. (1970). Effects of local geology on ground motion near San Francisco Bay. Bulletin of the Seismological Society of America, 60(1), 29–61. https://doi.org/10.1785/bssa0600010029
[20] Bard, P. Y., & Riepl-Thomas, J. (2000). Wave propagation in complex geological structures and their effects on strong ground motion. Wave motion in earthquake engineering, 37-95.
[21] Andrews, D. J. (1986). Objective determination of source parameters and similarity of earthquakes of different size. Earthquake source mechanics, 37, 259-267..
[22] Castro, R. R., Anderson, J. G., & Singh, S. K. (1990). Site response, attenuation and source spectra of S waves along the Guerrero, Mexico, subduction zone. Bulletin of the Seismological Society of America, 80(6A), 1481-1503.
[23] Boatwright, J., Seekins, L. C., Fumal, T. E., Liu, H. P., & Mueller, C. S. (1991). Ground motion amplification in the Marina District. Bulletin of the Seismological Society of America, 81(5), 1980-1997.
[24] Hartzell, S. H. (1992). Site response estimation from earthquake data. Bulletin of the Seismological Society of America, 82(6), 2308-2327.
[25] Nogoshi, M., & Igarashi, T. (1970). On the Propagation Characteristics of Microtremor. Journal of the Seismological Society of Japan, 23, 264-280.
[26] Nakamura, Y. (1989). A method for dynamic characteristics estimation of subsurface using microtremor on the ground surface. Railway Technical Research Institute, Quarterly Reports, 30(1).
[27] Lermo, J., & Chávez-García, F. J. (1993). Site effect evaluation using spectral ratios with only one station. Bulletin of the seismological society of America, 83(5), 1574-1594.
[28] Sylvette, B. C., Cécile, C., Pierre-Yves, B., Fabrice, C., Peter, M., Jozef, K., & Fäh, D. (2006). H/V ratio: a tool for site effects evaluation. Results from 1-D noise simulations. Geophysical Journal International, 167(2), 827-837.
[29] Bonnefoy-Claudet, S., Köhler, A., Cornou, C., Wathelet, M., & Bard, P. Y. (2008). Effects of Love waves on microtremor H/V ratio. Bulletin of the Seismological Society of America, 98(1), 288-300.
[30] Field, E. H., & Jacob, K. H. (1995). A comparison and test of various site-response estimation techniques, including three that are not reference-site dependent. Bulletin of the seismological society of America, 85(4), 1127-1143.
[31] Guéguen, P., Cornou, C., Garambois, S., & Banton, J. (2007). On the limitation of the H/V spectral ratio using seismic noise as an exploration tool: application to the Grenoble valley (France), a small apex ratio basin. Pure and applied geophysics, 164, 115-134.
[32] Poursartip, B., Fathi, A., & Tassoulas, J. L. (2020). Large-scale simulation of seismic wave motion: A review. Soil Dynamics and Earthquake Engineering, 129, 105909.
[33] Lamb, H. (1904). On the Propagation of Tremors over the Surface of an Elastic Solid. Proceedings of the royal society of London, 72, 128-130.
[34] Hisada, Y. (1994). An efficient method for computing Green's functions for a layered half-space with sources and receivers at close depths. Bulletin of the Seismological Society of America, 84(5), 1456-1472.
[35] Hisada, Y. (1995). An efficient method for computing Green's functions for a layered half-space with sources and receivers at close depths (Part 2). Bulletin of the Seismological Society of America, 85(4), 1080-1093.
[36] Garvin WW. Exact transient solution of the buried line source problem. Proc R Soc Lond Ser A Math Phys Sci 1956;234(1199):528–41.
[37] Thomson, W. T. (1950). Transmission of elastic waves through a stratified solid medium. Journal of applied Physics, 21(2), 89-93.
[38] Haskell, N. A. (1953). The dispersion of surface waves on multilayered media*. Bulletin of the Seismological Society of America, 43(1), 17–34.
[39] Haskell, N. A. (1960). Crustal reflection of plane SH waves. Journal of Geophysical Research, 65(12), 4147-4150.
[40] Haskell, N. A. (1962). Crustal reflection of plane P and SV waves. Journal of Geophysical Research, 67(12), 4751-4768.
[41] Huang, M., Cegla, F., & Lan, B. (2023). Stiffness matrix method for modelling wave propagation in arbitrary multilayers. International Journal of Engineering Science, 190, 103888.
[42] Poursartip, B., Fathi, A., & Kallivokas, L. F. (2017). Seismic wave amplification by topographic features: A parametric study. Soil Dynamics and Earthquake Engineering, 92, 503-527.
[43] Alterman, Z., & Karal Jr, F. C. (1968). Propagation of elastic waves in layered media by finite difference methods. Bulletin of the Seismological Society of America, 58(1), 367-398.
[44] Hill, N. R., & Levander, A. R. (1984). Resonances of low-velocity layers with lateral variations. Bulletin of the Seismological Society of America, 74(2), 521-537.
[45] Moczo, P. (1989). Finite-difference technique for SH-waves in 2-D media using irregular grids—application to the seismic response problem. Geophysical Journal International, 99(2), 321-329.
[46] Frankel, A., & Vidale, J. (1992). A three-dimensional simulation of seismic waves in the Santa Clara Valley, California, from a Loma Prieta aftershock. Bulletin of the Seismological Society of America, 82(5), 2045-2074.
[47] Aki, K., & Larner, K. L. (1970). Surface motion of a layered medium having an irregular interface due to incident plane SH waves. Journal of geophysical research, 75(5), 933-954.
[48] Bard, P. Y., & Bouchon, M. (1980). The seismic response of sediment-filled valleys. Part 1. The case of incident SH waves. Bulletin of the seismological society of America, 70(4), 1263-1286.
[49] Bard, P. Y., & Bouchon, M. (1980). The seismic response of sediment-filled valleys. Part 2. The case of incident P and SV waves. Bulletin of the Seismological Society of America, 70(5), 1921-1941.
[50] Geli, L., Bard, P. Y., & Jullien, B. (1988). The effect of topography on earthquake ground motion: a review and new results. Bulletin of the Seismological Society of America, 78(1), 42-63.
[51] Bouchon, M. (1973). Effect of topography on surface motion. Bulletin of the Seismological Society of America, 63(2), 615–632.
[52] Bard, P. Y. (1982). Diffracted waves and displacement field over two-dimensional elevated topographies. Geophysical Journal International, 71(3), 731-760.
[53] Zhao, C., & Valliappan, S. (1993). Seismic wave scattering effects under different canyon topographic and geological conditions. Soil Dynamics and Earthquake Engineering, 12(3), 129-143.
[54] Zhao, C., Valliappan, S., & Wang, Y. C. (1992). A numerical model for wave scattering problems in infinite media due to p‐and sv‐wave incidences. International Journal for Numerical Methods in Engineering, 33(8), 1661-1682.
[55] Bravo, M. A., & Sanchez-Sesma, F. J. (1990). Seismic response of alluvial valleys for incident P, SV and Rayleigh waves. Soil Dynamics and Earthquake Engineering, 9(1), 16-19.
[56] Moczo, P., Bystrický, E., Kristek, J., Carcione, J. M., & Bouchon, M. (1997). Hybrid modeling of P-SV seismic motion at inhomogeneous viscoelastic topographic structures. Bulletin of the seismological Society of America, 87(5), 1305-1323.
[57] Zahradník, J., & Moczo, P. (1996). Hybrid seismic modeling based on discrete-wave number and finite-difference methods. Seismic Waves in Laterally Inhomogeneous Media: Part 1, 21-38.
[58] Ohtsuki, A., & Harumi, K. (1983). Effect of topography and subsurface inhomogeneities on seismic SV waves. Earthquake Engineering & Structural Dynamics, 11(4), 441-462.
[59] Mossessian, T. K., & Dravinski, M. (1987). Application of a hybrid method for scattering of P, SV, and Rayleigh waves by near-surface irregularities. Bulletin of the Seismological Society of America, 77(5), 1784-1803.
[60] Mossessian, T. K., & Dravinski, M. (1992). A hybrid approach for scattering of elastic waves by three-dimensional irregularities of arbitrary shape. Journal of Physics of the Earth, 40(1), 241-261.
[61] Khair, K. R., Datta, S. K., & Shah, A. H. (1989). Amplification of obliquely incident seismic waves by cylindrical alluvial valleys of arbitrary cross-sectional shape. Part I. Incident P and SV waves. Bulletin of the Seismological Society of America, 79(3), 610-630.
[62] Khair, K. R., Datta, S. K., & Shah, A. H. (1991). Amplification of obliquely incident seismic waves by cylindrical alluvial valley of arbitrary cross-sectional shape. Part II. Incident SH and Rayleigh waves. Bulletin of the Seismological Society of America, 81(2), 346-357.
[63] Gatmiri, B., & Kamalian, M. (2002), "Two-dimensional transient wave propagation in an elastic saturated porous media by a hybrid FE/BE method," fifth European conference of numerical methods in geotechnical, paris,
[64] Gatmiri, B., & Kamalian, M. (2002). On the fundamental solution of dynamic poroelastic boundary integral equations in the time domain. International Journal of Geomechanics, 2(4), 381-398.
[65] Gatmiri, B., Arson, C., & Nguyen, K. V. (2008). Seismic site effects by an optimized 2D BE/FE method I. Theory, numerical optimization and application to topographical irregularities. Soil dynamics and earthquake engineering, 28(8), 632-645.
[66] Nguyen, K. V., & Gatmiri, B. (2007). Evaluation of seismic ground motion induced by topographic irregularity. Soil dynamics and earthquake engineering, 27(2), 183-188.
[67] Eshraghi, H., & Dravinski, M. (1989). Scattering of plane harmonic SH, SV, P and Rayleigh waves by non‐axisymmetric three‐dimensional canyons: A wave function expansion approach. Earthquake engineering & structural dynamics, 18(7), 983-998.
[68] Van Nguyen, K. (2005). Étude des effets de site dus aux conditions topographiques et géotechniques par une méthode hybride éléments finis/éléments frontières (Doctoral dissertation, Ecole des Ponts ParisTech).
[69] Gatmiri, B., Nguyen, K. V., & Dehghan, K. (2007). Seismic response of slopes subjected to incident SV wave by an improved boundary element approach. International journal for numerical and analytical methods in geomechanics, 31(10), 1183-1195.
[70] Gatmiri, B., Maghoul, P., & Arson, C. (2009). Site-specific spectral response of seismic movement due to geometrical and geotechnical characteristics of sites. Soil Dynamics and Earthquake Engineering, 29(1), 51-70.
[71] Razmkhah, A., Kamalian, M., & Sadrolddini, S. M. A. (2008, October). Time Domain Modeling of Topographic Effects on the Seismic Response of Slopes. In Proceedings of the 12th International Conference of International Association for Computer Methods and Advances in Geomechanics, Goa, India (pp. 2940-2947).
[72] Le Pense, S., Gatmiri, B., & Maghoul, P. (2011, July). Influence of soil properties and geometrical characteristics of sediment-filled valleys on earthquake response spectra. In 8th International Conference on Structural Dynamics (EURODYN 2011 (pp. 130-136).
[73] Gatmiri, B., & Foroutan, T. (2012). New criteria on the filling ratio and impedance ratio effects in seismic response evaluation of the partial filled alluvial valleys. Soil Dynamics and Earthquake Engineering, 41, 89-101.
[74] Kamalian, M., Gatmiri, B., Sohrabi‐Bidar, A., & Khalaj, A. (2007). Amplification pattern of 2D semi‐sine‐shaped valleys subjected to vertically propagating incident waves. Communications in Numerical Methods in Engineering, 23(9), 871-887.
[75] Gatmiri, B., & Amini-baneh, D. (2014). Impact of geometrical and mechanical characteristics on the spectral response of sediment-filled valleys. Soil Dynamics and Earthquake Engineering, 67, 233-250.
[76] Gatmiri, B., Amini-Baneh, D., Dorostkar, O., & Vakili, M. R. (2013). Practical recommendations of spectral response analysis in non-curved alluvial valleys using hybrid FE/BE method. Journal of Multiscale Modelling, 5(02), 1350006.

  • Receive Date 11 September 2023
  • Revise Date 10 October 2023
  • Accept Date 20 January 2024