Seismic response of the trapezoidal alluvial hill located on a circular cavity: Incident SH-wave.

Document Type : Research Article


1 Department of Civil Engineering, Islamic Azad University, Zanjan, Iran.

2 Department of Civil Engineering, Zanjan Branch, Islamic Azad University, Zanjan, Iran.


A direct time-domain numerical approach named the half-plane boundary element method is proposed based on the half-space Green’s functions for seismic analysis of trapezoidal alluvial hill located on a circular cavity, subjected to propagating vertical incident SH-waves. To analyze the assumed model, a half plane time-domain Boundary Element Method (BEM) was used which can concentrate the meshes only around the boundary of desired features. First, the problem is decomposed into two parts including a pitted half-plane and a trapezoidal filled solid on the surface. Then, the influence coefficients of the matrices are obtained by applying the method to each part. By satisfying the boundary/continuity conditions on the interfaces, a coupled equation is formed to determine unknown boundary values in each time-step. After implementing the method in an advanced developed algorithm, its efficiency is investigated by solving some practical examples and compared with those of the published works. To complete the results, the sensitivity analysis was carried out to obtain the seismic response of hill by considering the key parameters including impedance and shape ratios. In the meantime, the effect of subsurface cavity on the amplification pattern of surface has been studied as well. The results showed that the impedance and shape ratios of the trapezoidal alluvial hill were very effective on the seismic response of surface. The results of the present study can be used by geotechnical engineers to completing and increasing the accuracy of existing codes around the subject of ground surface zonation in presence of different topographic features.


[1] مجموعه ی استانداردها و  آیین نامههای ساختمانی ایران. مرکز تحقیقات راه،  مسکن و شهرسازی.  شماره ی نشر:  ض-253.  آیین نامه ی طراحی ساختمان ها در برابر زلزله. استاندارد 28۰۰ )ویرایش 4(، 18-19.
[2] Aki, K., Larner, K.L., (1970). “Surface motion of a layered medium having an irregular interface due to incident plane SH-waves”. J Geophys Res. 75(5):933954.
[3] Amornwongpaibun, A., Lee V.W., (2013). “Scattering of anti-plane SH-waves by a semi-elliptical hill: deep hill”. Soil Dyn Earthq Eng. 52:126-137.
[4] Amornwongpaibun, A., Luo, H., Lee, V.W., (2015). “Scattering of anti-plane SH-waves by a shallow semielliptical hill with a concentric elliptical tunnel”. J Earthq Eng. 20(3):363-382.
[5] Boor, D.M., (1972). “A note on the effect of simple topography on seismic SH-waves”. Bull Seism Soc Am. 62(1):275-284.
[6] Boor, D.M., (1973). “The effect of simple topography on seismic waves: Implications for the accelerations recorded at pacoima dam, San Fernando Valley, California”. Bull Seism Soc Am. 63(5):1603-1609.
[7] Bouchon, M., (1973). “Effect of topography on sur-.face motion”. Bull Seism Soc Am. 63(3):615-632.
[8]Bard, P., (1982). “Diffracted waves and displacement field over two-dimensional elevated topographies”. Geophys J Int. 71(3):731-760.
[9] Brebbia, C.A., Dominguez, J., (1989). “Boundary elements, an Introductory Course”. Comp Mech Pub. Southampton, Boston.
[10] Benites, R., Aki, K., Yomigida, K., (1992). “Multiple scattering of SH-waves in 2-D media with many cavities”. Pure appl Geophys 138:353-390.
[11] Ba, Z., Liang, J., Zhang, Y., (2016). “Scattering and diffraction of plane SH-waves by periodically distributed canyons”. Earthq Eng Eng Vib. 15(2):325-339.
[12] Ba, Z., Liang, J., Zhang, Y., (2017). “Diffraction of SH-waves by topographic features in a layered transversely isotropic half-space”. Earthq Eng Eng Vib, 16(1):11-22.
[13] Cao, X.R., Song, T.S., Liu, D.K., (2001). “Scattering of plane SH-wave by a cylindrical hill of arbitrary shape”. Appl Math Mech. 22(9):1082- 1089.
[14] Chen, J.T., Lee, J.W., Wu, C.F., Chen, I.L., (2011). “SH-wave diffraction by a semi-circular hill revisited: A null-field boundary integral equation method using degenerate kernels”. Soil Dyn Earthq Eng. 31(5-6):729736.
[15] Chen, J.T., Lee, J., Shyu, W., (2012). “SH-wave scattering by a semi-elliptical hill using a null-field boundary integral equation method and a hybrid method”. Geophys J Int. 188:177-194.
[16] Datta, S.K., (1974). “Diffractionof SH-waves by an elliptic elastic cylinder”. Int J Sol Struct. 10(1):123.331
[17] Day, S.M., (1977). “Finite element analysis of seismic scattering problems” [Dissertation]. University of California, San Diego.
[18] Dominguez, J., Meise, T., (1991). “On the use of the BEM for wave propagation in infinite domains”. Eng Analy BE. 8(3):132-138.
[19] Dominguez, J., (1993). “Boundary elements in dynamics”. Comp Mech Pub, Southampton, Boston.
[20] Eringen, A.C. & Suhubi, E.S., (1975). “Elastodynamics”. Academic Press.
[21] Geli, L., Bard, P., Jullien, B., (1988). “The effect of topography on earthquake ground motion: A review and new results”. Bull Seism Soc Am. 78(1):42-63.
[22] Hadley, P.K., (1987). “Scattering of waves by inclusions in a nonhomogeneous elastic half-space solved by boundary element method” [Dissertation]. Princeton University, Princeton, United States. 
[23] Kamalian, M., Gatmiri, B., Sohrabi-Bidar, A., (2003). “On time-domain two-dimensional site response analysis of topographic structures by BEM”. J Seism Earthq Eng. 5(2):35-45.
[24] Kamalian, M., Jafari, M.K., Sohrabi-Bidar, A., Razmkhah, A., Gatmiri, B., (2006a). “Time-domain two-dimensional site response analysis of non-homogeneous topographic structures by a hybrid FE/BE method”. Soil Dyn Earthq Eng. 26(8):753-765.
[25] Kamalian, M., Jafari, M.K., Sohrabi-Bidar, A., Razmkhah, A., (2006b). “Shape effects on amplification potential of two-dimensional hills” [In Persian]. J Seismolog Earthq Eng. 8(2):59-70.
[26] Lee, V.W., (1977). “On the deformations near circular underground cavity subjected to incident plane SH-waves”. In: Conference Proceedings of application of Computer Methods in Engineering. Los Angeles: Uni of South California. 9:5162-.
[27] Lee, V.W., Hao, L., Jianwen, L., (2004). “Diffraction of anti-plane SH-waves by a semi-circular cylindrical hill with an inside concentric semi-circular tunnel”. Earthq Eng Eng Vib. 3(2):249-262.
[28]Liang, J.W., Luo, H., Lee, V.W., (2004). “Scattering of plane SH-waves by a circular-arc hill with a circular tunnel”. Acta Seism Sinica. 17(5):549-563.
[29] Lee, V.W., Luo, H., Liang, J.W., (2006). “Antiplane SH-wave diffraction by a semicircular cylindrical hill revisited: An improved analytic wave series solution”. J Eng Mech. 132(10):1106-1114.
[30] Lin, S., Qiu, F., Liu, D., (2010). “Scattering of SHwaves by a scalene triangular hill”. Earthq Eng Eng Vib. 9(1):23-38.
[31] Liu, G., Chen, H., Liu, D., Khoo, B.C., (2010). “Surface motion of a half-space with triangular and semicircular hills under incident SH-waves”. Bull Seism Soc Am. 100(3):1306-1319.
[32] Lee, V.W., Amornwongpaibun, A., (2013). “Scattering of anti-plane SH-waves by a semi-elliptical hill: Shallow hill”. Soil Dyn Earthq Eng, 52:116-125.
[33]  Moczo, P., (1989). “Finite difference technique forSH-waves in 2D media using irregular grids-applica tion to the seismic response problem”. Geophys J Int. .99(2):321-329
[34] Mogi, H., Kawakami, H., (2007). “Analysis of scattered waves on ground with irregular topography using the direct boundary element method and neumann series expansion”. Bull seism Soc Am. 97(4):1144-1157.
[35] MATLAB, (2020). “The language of technical computing”. V9.8 .. Natick, Massachusetts: The MathWorks Inc.
[36] Ohtsu, M., Uesugi, S., (1985). “Analysis of SHwave scattering in a half space and its applications to seismic responses of geological structures”. Eng Analy. 2(4):198-204.
[37] Panji, M., Kamalian, M., Asgari-Marnani, J., Jafari, M.K., (2013). “Transient analysis of wave propagation problems by half-plane BEM”. Geophys J Int. 194:1849-1865.
[38] Panji, M., Kamalian, M., Asgari-Marnani, J., Jafari, M.K., (2014). “Analysing seismic convex topographies by a half-plane time-domain BEM”. Geophys J Int. 197(1):591-607.
[39] Panji, M., Ansari, B., (2017). “Transient SH-wave scattering by the lined tunnels embedded in an elastic half-plane”. Eng Analy BE. 84:220-230.
[40] Panji, M., Mojtabazadeh-Hasanlouei, S., (2018). “Time-history responses on the surface by regularly distributed enormous embedded cavities: Incident SHwaves”. Earthq Sci. 31:1-17.
[41] Panji, M., Mojtabazadeh-Hasanlouei, S., Yasemi, F., (2020). “A half-plane time-domain BEM for SH-wave scattering by a subsurface inclusion”. Comp Geosci. 10.1016/j.cageo.2019.104342.
[42] Panji, M., Mojtabazadeh-Hasanlouei, S., (2020). “Transient response of irregular surface by periodically distributed semi-sine shaped valleys: Incident SHwaves”. J Earthq Tsu. 10.1142/S1793431120500050.
[43] Panji, M., Habibivand, M., (2020). “Seismic Analysis of Semi-Sine Shaped Alluvial Hills above Subsurface Circular Cavity”. Earthq Eng Eng Vib. Accepted.
[44] Qiu, F., Liu, D., (2005). “Antiplane response of isosceles triangular hill to iIncident SH-Waves”. Earthq Eng Eng Vib. 4(1): 37-46
[45] Ricker, N. (1953). “The form and laws of propagation of seismic wavelets”. Geophys. 18(1):10-40.
[46] Reinoso, E., Wrobel, L.C., Power, H., (1993), “Preliminary results of the modeling of the Mexico City valley with a two-dimensional boundary element method for the scattering of SH-waves”. Soil Dyn Earthq Eng. 12(8), 457-468.
[47] Sabina, F.J., Willis, J.R., (1975). “Scattering of SHwaves by a rough half-space of arbitrary slope”. Geophys J R Astr Soc. 42(2):685-703.
[48] Smith, W.D., (1975). “The application of finite element analysis to body wave propagation problems”. Geophys J Int. 42:747-768.
[49] Sills, L.B., (1978). “Scattering of horizontally-po- larized shear waves by surface irregularities”. Geophys .J Int. 54(2):319-348
[50] Sánchez-Sesma, F.J., Herrera, A., Aviles, J., (1982). “A boundary method for elastic wave diffraction: Application to scattering of SH-waves by surface irregularities”. Bull Seism Soc Am. 72(2):473-490.
[51] Sánchez-Sesma, F.J., Palencia, V.J., Luzon, F., (2002). “Estimation of local site effects during earthquakes: An overview”. ISET J Earthq. 39(3):167-193.
[52] Shyu, W., Teng, T., (2014). “Hybrid method com- bines transfinite interpolation with series expansion to simulate the anti-plane response of a surface irregular- ity”. J Mech. 4:349-360.
[53] Shyu, W., Teng, T., Chou, C., (2017). “Anti-plane response caused by interactions between a dike and the surrounding soil”. Soil Dyn Earthq Eng. 92:408-418
[54] Thambiratnam, D.P., Lee, S.L., (1990). “Scatter- ing of plane SH-waves by underground cavities”. Eng Struct. 12(3):215-221.
[55] Takemiya, H., Fujiwara, A., (1994). “SH-Wave scattering and propagation analysis at irregular sites by time-domain BEM”. Bull Seism Soc Am, 84(5), 1443- 1455.
[56] Tsaur, D., Chang, K., (2009). “Scattering and fo- cusing of SH-waves by a convex circular-arc topogra- phy”. Geophys J Int. 177(1):222-234.
[57] Tsaur, D., Chang, K., Hsu, M., (2018). “Ground motions around a deep semielliptic canyon with a hori- zontal edge subjected to incident plane SH-waves. J Seismol. 22(6):1579-1593.
[58] Virieux, J., (1984). “SH-wave propagation in het- erogeneous media: Velocity-stress finite-difference method”. Geophys. 49(2);1933-1957.
[59]Yuan, X., Men, F.L., (1992). “Scattering of plane SH-waves by a semi-cylindrical hill”. Earthq Eng Struct Dyn. 21(12):1091-1098.
[60] Yuan, X., Liao, Z., (1996). “Surface motion of a cylindrical hill of circular arc cross section for incident plane SH-waves”. Soil Dyn Earthq Eng,15(3):189-199.
[61] Yang, Z., Xu, H., Hei, B., Zhang, J., (2014). “An- tiplane response of two scalene triangular hills and a semi-cylindrical canyon by incident SH-waves”. Earthq Eng Eng Vib. 13(4): 569-581.
[62] Yang, Z., Song, Y., Li, X., Jiang, G., Yang, Y., (2019). “Scattering of plane SH-waves by an isosceles trapezoidal hill”. Wave Motion. 92, 102415.
[63] Zhou, H., Chen, X., (2006). “A new approach to simulate scattering of SH-waves by an irregular topog- raphy”. Geophys J Int. 164(2):449-459.