<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Earth Sci.</journal-id>
<journal-title>Frontiers in Earth Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Earth Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-6463</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1057316</article-id>
<article-id pub-id-type="doi">10.3389/feart.2022.1057316</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Earth Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Seismic response analysis of slope sites exposed to obliquely incident P waves</article-title>
<alt-title alt-title-type="left-running-head">Jiao et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/feart.2022.1057316">10.3389/feart.2022.1057316</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Jiao</surname>
<given-names>Hongyun</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2026612/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Xie</surname>
<given-names>Junju</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1834759/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhao</surname>
<given-names>Mi</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Huang</surname>
<given-names>Jingqi</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Du</surname>
<given-names>Xiuli</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Juke</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1506502/overview"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Institute of Geophysics</institution>, <institution>China Earthquake Administration</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Key Laboratory of Urban Security and Disaster Engineering</institution>, <institution>Ministry of Education</institution>, <institution>Beijing University of Technology</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Beijing Key Laboratory of Urban Underground Space Engineering</institution>, <institution>School of Civil and Resource Engineering</institution>, <institution>University of Science and Technology Beijing</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1652730/overview">Yefei Ren</ext-link>, Institute of Engineering Mechanics, China Earthquake Administration, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2054057/overview">Zhang Yushan</ext-link>, China Earthquake Disaster Prevention Centre, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1840955/overview">Wujian Yan</ext-link>, Lanzhou Earthquake Research Institute, China Earthquake Administration, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Junju Xie, <email>xiejunjv05@mails.ucas.ac.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Structural Geology and Tectonics, a section of the journal Frontiers in Earth Science</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>04</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2022</year>
</pub-date>
<volume>10</volume>
<elocation-id>1057316</elocation-id>
<history>
<date date-type="received">
<day>29</day>
<month>09</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>05</day>
<month>12</month>
<year>2022</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Jiao, Xie, Zhao, Huang, Du and Wang.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Jiao, Xie, Zhao, Huang, Du and Wang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>This study proposes a seismic input method for layered slope sites exposed to obliquely-incident seismic waves which transforms the waves into equivalent nodal forces that act on the truncated boundary of a finite element model. The equivalent nodal forces at the left and right boundaries are obtained by combining the free field response of a one-dimensional layered model with a viscoelastic boundary. The equivalent nodal forces at the bottom boundary are obtained by combining the incident wave field with the viscoelastic boundary. This proposed seismic input method for slope sites exposed to obliquely incident seismic waves is implemented with the aid of MATLAB software; it is applied to the seismic response analysis of slope sites in the commercial finite element ABAQUS software. The calculation results are compared with the reference solutions obtained by using the extended model to verify the correctness of the established seismic input method. The proposed seismic input method is then employed to investigate the influencing factors of the seismic response of layered slope sites exposed to oblique incidence P waves. The results show that the angle of incidence, location of the interface between soft and hard rocks, and impedance ratio have significant effects on the seismic landslide.</p>
</abstract>
<kwd-group>
<kwd>seismic input method</kwd>
<kwd>obliquely incident P waves</kwd>
<kwd>layered slope site</kwd>
<kwd>influencing factors</kwd>
<kwd>seismic landslide</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Landslides are frequent during earthquakes in mountainous and hilly areas (<xref ref-type="bibr" rid="B34">Prestininzi and Romeo, 2000</xref>; <xref ref-type="bibr" rid="B9">Chigira et al., 2005</xref>; <xref ref-type="bibr" rid="B36">Sato et al., 2007</xref>; <xref ref-type="bibr" rid="B37">Semblat et al., 2011</xref>). Hillside topography can magnify seismic intensity and change seismic frequency content&#x2014;termed &#x201c;topographic effects&#x201d;. Landslide disasters caused by earthquakes have been a frequent subject of geological hazard research because of their wide distribution, considerable quantity, and great harm (<xref ref-type="bibr" rid="B6">Bird et al., 2004</xref>; <xref ref-type="bibr" rid="B30">Owen et al., 2008</xref>). Since the 1960s, scholars have begun identifying and analyzing the seismic response of slope sites exposed to seismic waves (<xref ref-type="bibr" rid="B8">Cavallin and Slejko, 1986</xref>; <xref ref-type="bibr" rid="B17">Jibson et al., 2000</xref>; <xref ref-type="bibr" rid="B31">Pareek and Arora, 2010</xref>; <xref ref-type="bibr" rid="B21">Li et al., 2022</xref>). Analysis methods have included landslide observation, model testing, and numerical simulations (<xref ref-type="bibr" rid="B19">Keeper, 1984</xref>). The finite element method (FEM) is the most common numerical method; it can effectively simulate the geometric and material non-linear characteristics of slope sites.</p>
<p>When analyzing the seismic response of slope sites based on FEM, it is necessary to introduce artificial boundary conditions at the truncated boundaries of the finite domain to simulate the radiation damping effect of the infinite domain on the finite domain. The artificial boundary appropriate to the characteristics and specific conditions of the research object must be deduced. Due to the irregular terrain of slope sites, poor geological conditions, and the existence of free surfaces, the topographic effect of slope site has attracted research attention (<xref ref-type="bibr" rid="B13">Gischig et al., 2015</xref>; <xref ref-type="bibr" rid="B33">Poursartip et al., 2017</xref>). It is necessary to deduce the appropriate artificial boundary according to the characteristics of the research object and the specific situation. In previous finite element analysis of the seismic response of slope sites, most scholars have assumed that the bottom boundary of the model is rigid, while the lateral side boundaries adopt roller boundaries (<xref ref-type="bibr" rid="B35">Rizzitano et al., 2014</xref>), viscous boundaries (<xref ref-type="bibr" rid="B27">Lysmer and Kuhlemeyer, 1969</xref>; <xref ref-type="bibr" rid="B3">Athanasopoulo et al., 1999</xref>; <xref ref-type="bibr" rid="B32">Pelekis, 2017</xref>), viscoelastic boundaries (<xref ref-type="bibr" rid="B11">Deeks and Randolph, 1994</xref>; <xref ref-type="bibr" rid="B24">Liu et al., 2006</xref>; <xref ref-type="bibr" rid="B28">Maleki and Khodakarami, 2017</xref>), transmissive boundaries (<xref ref-type="bibr" rid="B23">Liao and Wong, 1984</xref>), paraxial approximate boundaries (<xref ref-type="bibr" rid="B10">Clayton and Engquist, 1977</xref>), and infinite element boundaries (<xref ref-type="bibr" rid="B5">Bettess, 1977</xref>; <xref ref-type="bibr" rid="B2">Astley, 2000</xref>). Seismic input is completed by converting the seismic wave action into the equivalent forces applied to boundaries; however, this treatment method is not applicable in cases where the lower part of the site model is not bedrock.</p>
<p>Much recent research has been conducted on the seismic response analysis of slope sites with non-rigid bedrock at the bottom of the site model. Based on the viscous boundary, <xref ref-type="bibr" rid="B7">Bouckovalas and Papadimitriou (2005)</xref> presented numerical analyses for the seismic response of step-like ground slopes in uniform viscoelastic soil under vertically propagating SV seismic waves. <xref ref-type="bibr" rid="B1">Assimaki et al. (2005)</xref> employed the viscous boundary to address vertically incident seismic wave action and obtained a seismic input method suitable for layered slopes. Based on the viscous boundary, <xref ref-type="bibr" rid="B20">Lenti and Martino (2012)</xref> studied a landslide disaster on a stepped slope under the action of a vertically incident seismic wave. <xref ref-type="bibr" rid="B29">Nakamura (2012)</xref> adopted the energy-transmitting boundary to study the seismic input and established the seismic response analysis method for the layered slope. However, when the source is shallow or the site is far from the epicenter, the seismic input cannot be assumed to be vertically incident seismic waves but can be considered obliquely incident.</p>
<p>Scholars have thus conducted research on the seismic response of regular sites under the action of obliquely incident seismic waves. <xref ref-type="bibr" rid="B25">Liu and Lu (1998)</xref> proposed a seismic input method based on spring-buffer boundary conditions that convert seismic waves with arbitrary incident angles into equivalent nodal forces acting on boundary nodes. Based on the viscoelasticity boundary, <xref ref-type="bibr" rid="B16">Huang et al. (2017a)</xref> and <xref ref-type="bibr" rid="B15">Huang et al. (2017b)</xref> used the FEM to analyze the non-linear seismic response of tunnels with a normal fault ground subjected to obliquely incident seismic waves. <xref ref-type="bibr" rid="B4">Bazyar and Song (2017)</xref> expressed such a wave as the boundary condition applied to the near field by the proportional boundary FEM. Based on the one-dimensional time-domain FEM proposed by <xref ref-type="bibr" rid="B26">Liu and Wang (2007)</xref>, <xref ref-type="bibr" rid="B40">Zhao et al. (2013)</xref> and <xref ref-type="bibr" rid="B39">Zhao et al. (2017)</xref> proposed an improvement by establishing a site response analysis method and applying it to study the influence of the oblique incidence of ground motion on the seismic response of subway stations. However, previous seismic input methods have been established for seismic response analysis of regular sites under the action of obliquely incident seismic waves; this is not applicable to the seismic response of slope sites under the action of obliquely incident seismic waves. There are few studies on the seismic input method of slope sites exposed to obliquely incident seismic waves.</p>
<p>In this paper, obliquely incident seismic waves are transformed into equivalent nodal forces that act on the truncated boundary of the finite element model as the seismic input. Based on the viscoelastic artificial boundary, the free field responses of one-dimensional layered sites with different heights are taken as the seismic input for the left and right boundaries, and the incident wave field is used as the seismic input for the bottom. Then, the proposed seismic input method for slope sites under the action of obliquely incident seismic waves is implemented with the aid of MATLAB software and is applied to the seismic response analysis of slope sites in the commercial finite element ABAQUS software. Accordingly, the numerical simulation results of a two-dimensional layered slope site obtained using the established seismic input method are compared with the numerical results of an extended computational model to verify the accuracy of this method. Finally, with the aid of the proposed seismic input method, the influencing factors of the seismic response of layered slope sites under oblique incidence P waves are determined.</p>
</sec>
<sec id="s2">
<title>2 Seismic input mechanism for the slope site</title>
<sec id="s2-1">
<title>2.1 Governing equations</title>
<p>
<xref ref-type="fig" rid="F1">Figure 1A</xref> illustrates the schematic diagram of the seismic response analysis of a layered slope site subjected to P waves obliquely incident at an angle of <italic>&#x3b1;</italic>. When analyzing the seismic response of a two-dimensional layered slope site with the aid of the FEM, the finite computational domain needs to be cut off from the infinite ground. Artificial boundaries are usually used to simulate waves scattered by target structures and the non-reflecting wave effect of truncated infinite domains (<xref ref-type="bibr" rid="B12">Du et al., 2006</xref>; <xref ref-type="bibr" rid="B14">Huang et al., 2016</xref>; <xref ref-type="bibr" rid="B38">Zhao et al., 2019</xref>; <xref ref-type="bibr" rid="B22">Li et al., 2020</xref>). The stable and accurate viscoelasticity artificial boundary developed by <xref ref-type="bibr" rid="B12">Du et al. (2006)</xref> is adopted in this study. For a given boundary node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, <italic>y</italic>
<sub>
<italic>l</italic>
</sub>, <italic>z</italic>
<sub>
<italic>l</italic>
</sub>), one pair of dashpots and springs in the normal and tangential directions of the boundary plane are established respectively (<xref ref-type="fig" rid="F1">Figure 1B</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Model of a layered slope site under obliquely incident P waves. <bold>(A)</bold> Schematic diagram. <bold>(B)</bold> Proposed seismic input method.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g001.tif"/>
</fig>
<p>At a given boundary node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, <italic>y</italic>
<sub>
<italic>l</italic>
</sub>) of the 2-D finite element model, the parameters for the springs and dashpots are<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi mathvariant="italic">ln</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>G</mml:mi>
</mml:mrow>
<mml:mi>R</mml:mi>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>ln</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mfrac>
<mml:mrow>
<mml:mi>G</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>b</mml:mi>
<mml:mi>r</mml:mi>
</mml:msub>
<mml:mi>&#x3c1;</mml:mi>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>s</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>where <italic>K</italic>
<sub>
<italic>ln</italic>
</sub> and <italic>K</italic>
<sub>
<italic>ls</italic>
</sub> denote the normal and tangential spring stiffnesses at the boundary node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, <italic>y</italic>
<sub>
<italic>l</italic>
</sub>); <italic>C</italic>
<sub>
<italic>ln</italic>
</sub> and <italic>C</italic>
<sub>
<italic>ls</italic>
</sub> denote the normal and tangential damping coefficients at the boundary node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, <italic>y</italic>
<sub>
<italic>l</italic>
</sub>); <italic>A</italic>
<sub>
<italic>l</italic>
</sub> is half of the total length of all boundary elements containing the boundary node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, <italic>y</italic>
<sub>
<italic>l</italic>
</sub>); <italic>R</italic> is the distance between the scattering source and the artificial boundary node; <italic>&#x3bb;</italic>, <italic>G</italic>, and <italic>&#x3c1;</italic> are the Lam&#xe9; constant, shear modulus, and mass density of the ground, respectively; <italic>c</italic>
<sub>
<italic>s</italic>
</sub> and <italic>c</italic>
<sub>
<italic>p</italic>
</sub> represent the shear and compression wave velocities of the ground, respectively; <italic>a</italic>
<sub>
<italic>r</italic>
</sub> and <italic>b</italic>
<sub>
<italic>r</italic>
</sub> are the modified coefficients with good values of 0.8 and 1.1, respectively (<xref ref-type="bibr" rid="B12">Du et al., 2006</xref>).</p>
<p>The dynamic finite element equation for the finite domain has the form <disp-formula id="e2">
<mml:math id="m2">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>II</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>IB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>BI</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>BB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
<mml:mi mathvariant="normal">I</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi>II</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi>IB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi>BI</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi>BB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi mathvariant="normal">I</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mrow>
<mml:mi>II</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mrow>
<mml:mi>IB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mrow>
<mml:mi>BI</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mrow>
<mml:mi>BB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi mathvariant="normal">I</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold">I</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold">B</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where subscript <sub>B</sub> represents the boundary node and subscript <sub>I</sub> represents the internal node, <inline-formula id="inf1">
<mml:math id="m3">
<mml:mrow>
<mml:mi mathvariant="bold-italic">u</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf2">
<mml:math id="m4">
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf3">
<mml:math id="m5">
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf4">
<mml:math id="m6">
<mml:mrow>
<mml:mi mathvariant="bold-italic">F</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> represent displacement, velocity, displacement and force vectors, and <bold>
<italic>M</italic>
</bold>, <bold>C</bold> and <bold>K</bold> represent the mass matrix, damping matrix, and stiffness matrix.</p>
<p>The seismic waves can be separated into the free field motion of the infinite domain and scattering wave motion generated by all scattering sources. The free field motion in the infinite domain or engineering site is expressed by the superscript <sup>
<italic>f</italic>
</sup>, and the scattering wave motion is expressed by the superscript <sup>
<italic>s</italic>
</sup>. The displacement <inline-formula id="inf5">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, velocity <inline-formula id="inf6">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and force <inline-formula id="inf7">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> at a given boundary node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, <italic>y</italic>
<sub>
<italic>l</italic>
</sub>) of the 2-D finite element model can be written as<disp-formula id="e3">
<mml:math id="m10">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
<disp-formula id="e4">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<disp-formula id="e5">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mi mathvariant="bold-italic">f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
</mml:msubsup>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>where <inline-formula id="inf8">
<mml:math id="m13">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold">u</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf9">
<mml:math id="m14">
<mml:mrow>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi mathvariant="bold">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula id="inf10">
<mml:math id="m15">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold">F</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>T</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>. The load <inline-formula id="inf11">
<mml:math id="m16">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold">F</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> due to the scattering wave motion can be expressed as<disp-formula id="e6">
<mml:math id="m17">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mi mathvariant="bold-italic">s</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:msubsup>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>where <inline-formula id="inf12">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf13">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> represent the spring component and damping component of the stiffness matrix <inline-formula id="inf14">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">K</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and damping matrix <inline-formula id="inf15">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold">C</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> at an artificial boundary node <italic>l</italic> in the <italic>i</italic> direction, respectively.</p>
<p>Substituting Eqs <xref ref-type="disp-formula" rid="e3">3</xref>&#x2013;<xref ref-type="disp-formula" rid="e5">5</xref> into Eq. <xref ref-type="disp-formula" rid="e6">6</xref> after rearrangement, the following equation can be obtained:<disp-formula id="e7">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">A</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mi mathvariant="bold-italic">f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>where <inline-formula id="inf16">
<mml:math id="m23">
<mml:mrow>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mi>x</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mtext>&#x2002;</mml:mtext>
<mml:msubsup>
<mml:mi>f</mml:mi>
<mml:mi>y</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>T</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula> is the stress vector induced by free field motion.</p>
<p>Then, substituting Eq. <xref ref-type="disp-formula" rid="e7">7</xref> into Eq. <xref ref-type="disp-formula" rid="e2">2</xref>, the dynamic finite element equation is:<disp-formula id="e8">
<mml:math id="m24">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>II</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>IB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>BI</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">M</mml:mi>
<mml:mrow>
<mml:mi>BB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
<mml:mi mathvariant="normal">I</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#xa8;</mml:mo>
</mml:mover>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi>II</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi>RI</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi>BI</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mrow>
<mml:mi>BB</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi mathvariant="normal">I</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mrow>
<mml:mi>II</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mrow>
<mml:mi>IB</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mrow>
<mml:mi>BI</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mrow>
<mml:mi>BB</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi mathvariant="normal">I</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msub>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi mathvariant="normal">B</mml:mi>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="{" close="}" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mn>0</mml:mn>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msup>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold-italic">f</mml:mi>
</mml:msup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>where submatrices [<bold>
<italic>C</italic>
</bold>
<sub>B</sub>] and [<bold>
<italic>K</italic>
</bold>
<sub>B</sub>] are both diagonal. For the boundary node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, <italic>y</italic>
<sub>
<italic>l</italic>
</sub>), [<bold>C</bold>
<sub>B</sub>]<sub>
<italic>l</italic>
</sub> &#x3d; <bold>
<italic>C</italic>
</bold>
<sub>
<italic>l</italic>
</sub>, [<bold>
<italic>K</italic>
</bold>
<sub>B</sub>]<sub>
<italic>l</italic>
</sub> &#x3d; <bold>
<italic>K</italic>
</bold>
<sub>
<italic>l</italic>
</sub>. [<bold>
<italic>F</italic>
</bold>
<sup>
<italic>f</italic>
</sup>] represents the equivalent seismic forces acting on the artificial boundary. <inline-formula id="inf17">
<mml:math id="m25">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="normal">F</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the equivalent nodal force and has the form<disp-formula id="e9">
<mml:math id="m26">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">F</mml:mi>
<mml:mi mathvariant="bold-italic">l</mml:mi>
<mml:mi mathvariant="bold-italic">f</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">K</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">C</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mover accent="true">
<mml:mi mathvariant="bold-italic">u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">f</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>f</mml:mi>
</mml:msubsup>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>Equation <xref ref-type="disp-formula" rid="e9">9</xref> provides the method of converting seismic waves into equivalent nodal forces that act on artificial boundary nodes to realize the seismic wave input in the seismic response analysis of the finite element model. To calculate the equivalent node force at each boundary node, it is necessary to provide the free field velocity, displacement, and stress at the node.</p>
</sec>
<sec id="s2-2">
<title>2.2 Equivalent nodal forces for the slope site</title>
<p>For the incidence of the seismic waves, the wave motions can be transferred into equivalent node forces applied at the boundary nodes. However, due to the topographic effect of the slope site, the seismic input wave fields at the left, right, and bottom boundaries of the two-dimensional finite element model are different. Realization of the seismic wave input at the left and right artificial boundaries requires the free field response calculation for a one-dimensional layered site. Because of the different heights of the left and right sides of the slope terrain, one-dimensional FEM calculations are necessary (<xref ref-type="fig" rid="F2">Figure 2</xref>) (<xref ref-type="bibr" rid="B39">Zhao et al., 2017</xref>). The equivalent nodal forces at the left and right boundaries are then obtained by combining the calculation results of the one-dimensional FEM with the viscoelastic artificial boundary to complete the seismic input at the left and right boundaries.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>One-dimensional finite element model. <bold>(A)</bold> Left boundary. <bold>(B)</bold> Right boundary.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g002.tif"/>
</fig>
<p>For node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, <italic>y</italic>
<sub>
<italic>l</italic>
</sub>) on the left or right boundary, the equivalent nodal forces of obliquely incident P waves can be given as<disp-formula id="e10">
<mml:math id="m27">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi mathvariant="italic">ln</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi mathvariant="italic">ln</mml:mi>
</mml:msub>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>y</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>where <inline-formula id="inf18">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf19">
<mml:math id="m29">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula id="inf20">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the displacement, velocity, and internal stress of the one-dimensional finite element model at node <italic>l</italic> in the <italic>x</italic> direction respectively; <inline-formula id="inf21">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf22">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula id="inf23">
<mml:math id="m33">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the displacement, velocity, and internal stress of the one-dimensional finite element model at node <italic>l</italic> in the <italic>y</italic> direction, respectively; the boundary-dependent parameter <italic>&#x3b4;</italic> values are -1 for the left boundary and 1 for the right boundary.</p>
<p>For node <italic>l</italic> (<italic>x</italic>
<sub>
<italic>l</italic>
</sub>, 0) on the bottom boundary, the equivalent nodal forces of obliquely incident P waves can be given as<disp-formula id="e11">
<mml:math id="m34">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>x</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>l</mml:mi>
<mml:mi>s</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>y</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>K</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>ln</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mi>ln</mml:mi>
</mml:msub>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>where <inline-formula id="inf24">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf25">
<mml:math id="m36">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the displacement and velocity of the bottom boundary node in <italic>x</italic> the direction, respectively; <inline-formula id="inf26">
<mml:math id="m37">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf27">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the displacement and velocity of the bottom boundary node in the <italic>y</italic> direction, respectively; <inline-formula id="inf28">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf29">
<mml:math id="m40">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> are the stresses of the bottom boundary node in the <italic>x</italic> and <italic>y</italic> directions, respectively.</p>
<p>Since the P waves are obliquely incident at an angle of <italic>&#x3b1;</italic>, there is a time delay in the response of any node on the bottom boundary relative to that of the initial incident node in the finite element model. <inline-formula id="inf30">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf31">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> are the displacement and velocity at the initial incident node with an angle of <italic>&#x3b1;</italic>, respectively. Therefore, the displacement and velocity of each node at the bottom boundary are expressed as<disp-formula id="e12">
<mml:math id="m43">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>u</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>
<disp-formula id="e13">
<mml:math id="m44">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(13)</label>
</disp-formula>where <inline-formula id="inf32">
<mml:math id="m45">
<mml:mrow>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the time delay in the response of any node on the bottom boundary relative to that of the initial incident node.</p>
<p>The stresses along and perpendicular to the incident direction of the P wave are expressed as<disp-formula id="e14">
<mml:math id="m46">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>G</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(14)</label>
</disp-formula>
</p>
<p>The stresses <inline-formula id="inf33">
<mml:math id="m47">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf34">
<mml:math id="m48">
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> of the bottom boundary node in the <italic>x</italic> and <italic>y</italic> directions are<disp-formula id="e15">
<mml:math id="m49">
<mml:mrow>
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="left">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>G</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>f</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>&#x3c3;</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x22c5;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>sin</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>2</mml:mn>
<mml:mi>G</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>cos</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
<mml:msub>
<mml:mi>c</mml:mi>
<mml:mi>p</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:msub>
<mml:mover accent="true">
<mml:mi>u</mml:mi>
<mml:mo>&#x2d9;</mml:mo>
</mml:mover>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mo>&#x2206;</mml:mo>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>.</mml:mo>
</mml:mrow>
</mml:math>
<label>(15)</label>
</disp-formula>
</p>
</sec>
</sec>
<sec id="s3">
<title>3 Method validation</title>
<p>The seismic input method established in <xref ref-type="sec" rid="s2">Section 2</xref> is realized through programming with MATLAB software and is applied to the seismic response analysis of slope sites in the commercial finite element ABAQUS software. To verify the accuracy of the proposed seismic input method, this section simulates the free field response of a layered slope site under the action of obliquely incident P waves. The size of the finite element model of the two-dimensional slope site is 240&#xa0;m long, 70&#xa0;m high on the left side, and 100&#xa0;m high on the right side (<xref ref-type="fig" rid="F3">Figure 3</xref>). <italic>h</italic> represents the thickness of soft layer; in this section, <italic>h &#x3d;</italic> 30&#xa0;m. The reference solution is the calculation result of an extended computational model with a length of 4100&#xa0;m, a right-side height of 2100&#xa0;m, and a left-side height of 2070 m. The slope site is meshed by solid elements, with a mesh size of 1&#xa0;m. The material parameters of each layer of the slope site are shown in <xref ref-type="table" rid="T1">Table 1</xref> (<xref ref-type="bibr" rid="B18">JTG/T D70-2010, 2010</xref>). The form of the obliquely incident pulse wave with an incident angle of <italic>&#x3b1;</italic> acting on the slope site calculation model is shown in <xref ref-type="fig" rid="F4">Figure 4</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Calculation model for the layered slope site exposed to obliquely incident P waves.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g003.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Material constants of the soil.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Layer</th>
<th align="center">Density <italic>&#x3c1;</italic> (kg/m<sup>3</sup>)</th>
<th align="center">Elastic modulus <italic>E</italic> (Pa)</th>
<th align="center">Poisson&#x2019;s ratio <italic>&#x3c5;</italic>
</th>
<th align="center">Internal friction angle <italic>&#x3c6;</italic> (&#xb0;)</th>
<th align="center">Cohesion <italic>c</italic> (Pa)</th>
<th align="center">
<italic>C</italic>
<sub>
<italic>P</italic>
</sub> (m/s)</th>
<th align="center">
<italic>C</italic>
<sub>
<italic>SV</italic>
</sub> (m/s)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Upper</td>
<td align="center">1700</td>
<td align="center">1e9</td>
<td align="center">0.45</td>
<td align="center">20</td>
<td align="center">0.05e6</td>
<td align="center">1494</td>
<td align="center">450</td>
</tr>
<tr>
<td align="center">Lower</td>
<td align="center">2400</td>
<td align="center">17.7e9</td>
<td align="center">0.26</td>
<td align="center">45</td>
<td align="center">1.4e6</td>
<td align="center">3004</td>
<td align="center">1711</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Time history curve of the input seismic wave. <bold>(A)</bold> Displacement. <bold>(B)</bold> Velocity.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g004.tif"/>
</fig>
<p>In verifying the seismic input method for the layered slope site, the P wave is incident from the left corner of the model with angles of 0&#xb0;, 10&#xb0;, 20&#xb0;, and 30&#xb0;. The observation nodes are arranged at each node of the free surface within the range of 10&#xa0;m from the left side of the slope toe to 10&#xa0;m from the right side of the slope top (<xref ref-type="fig" rid="F3">Figure 3</xref>). Top node A and toe node B are the main observation nodes. The results of the acceleration, velocity, and displacement time history of the layered slope site model, calculated on the established seismic input method, are compared with the results of the extended model. The comparison results are shown in <xref ref-type="fig" rid="F5">Figure 5</xref>. The results obtained with the aid of the established seismic input method agree well with the reference solutions. The accuracy and applicability of the established seismic wave input method for slope sites are thus verified.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Comparison results at different input angles. <bold>(A)</bold> Comparison of the vertical acceleration results of node A at different input angles. <bold>(B)</bold> Comparison of the vertical velocity results of node A at different input angles. <bold>(C)</bold> Comparison of the vertical displacement results of node A at different input angles. <bold>(D)</bold> Comparison of the vertical acceleration results of node B at different input angles. <bold>(E)</bold> Comparison of the vertical velocity results of node B at different input angles. <bold>(F)</bold> Comparison of the vertical displacement results of node B at different input angles.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g005.tif"/>
</fig>
<p>The relative error of the peak displacement value is used to measure the accuracy of the calculation results; the calculation formula is <disp-formula id="e16">
<mml:math id="m50">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mi>r</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>r</mml:mi>
<mml:mn>0</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mi>max</mml:mi>
</mml:msub>
</mml:mfrac>
<mml:mo>,</mml:mo>
</mml:mrow>
</mml:math>
<label>(16)</label>
</disp-formula>where <italic>r</italic>
<sub>0</sub>(<italic>t</italic>) is the displacement reference solution, <italic>r</italic>(<italic>t</italic>) is the displacement calculation result obtained using the seismic input method proposed in this paper, &#x7c; &#x7c; represents the absolute value, and the subscript <sub>max</sub> is used to evaluate the maximum value.</p>
<p>The comparison of the results and relative errors of the vertical displacement peaks of each observation node on the free surface are shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. The change in the dynamic amplification effect of the displacement peak with topographic relief is negligible, and the maximum error of the displacement peak at four incident angles is less than 1%. Therefore, by comparing the calculated results of the small model of the slope site with the numerical results of the extended model, the high accuracy of the developed seismic input method is further demonstrated.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Comparison of the results and relative errors of the vertical displacement peaks of each observation node on the free surface. <bold>(A)</bold> Comparison of the results of the vertical displacement peaks. <bold>(B)</bold> Relative errors.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g006.tif"/>
</fig>
</sec>
<sec id="s4">
<title>4 Parameter study</title>
<p>This section discusses the effects of incident angles, soft-hard rock interface positions, and impedance ratios on the seismic responses of slope sites exposed to obliquely incident P waves. The seismic input method of the slope site established in <xref ref-type="sec" rid="s2">Section 2</xref> is employed.</p>
<p>The upper soft rock and lower hard rock of the slope site is modeled in ABAQUS software, and the Mohr-Coulomb elastic-plastic model is selected for the constitutive relationship. The model dimensions are shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. The Kobe wave (<xref ref-type="fig" rid="F7">Figure 7</xref>) is input by the seismic input method proposed in this paper to study the factors that influence seismic landslides.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Time histories of the Kobe wave for Kobe University. <bold>(A)</bold> Time history of displacement. <bold>(B)</bold> Time history of acceleration.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g007.tif"/>
</fig>
<sec id="s4-1">
<title>4.1 Input angle of the seismic wave</title>
<p>This subsection investigates the influence of the different incident angles (<italic>&#x3b1;</italic> &#x3d; 0&#xb0;, 15&#xb0;, 30&#xb0;, and 45&#xb0;) of P waves on the dynamic responses of the slope site. The site parameters are shown in <xref ref-type="table" rid="T1">Table 1</xref>, with the soft and hard rock interface 30&#xa0;m from the slope top. <xref ref-type="fig" rid="F8">Figure 8</xref> shows the plastic strain nephograms of the slope site at different timing exposed to obliquely incident P waves with different angles. The displacement nephograms of the slope site at <italic>t</italic> &#x3d; 20&#xa0;s are illustrated in <xref ref-type="fig" rid="F9">Figure 9A</xref>.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Plastic strain nephograms of the slope at different timing with different incident angles. <bold>(A)</bold> <italic>&#x3b1;</italic> &#x3d; 0&#xb0;. <bold>(B)</bold> <italic>&#x3b1;</italic> &#x3d; 15&#xb0;. <bold>(C)</bold> <italic>&#x3b1;</italic> &#x3d; 30&#xb0;. <bold>(D)</bold> <italic>&#x3b1;</italic> &#x3d; 45&#xb0;.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g008.tif"/>
</fig>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Displacement of the slope at different incident angles. <bold>(A)</bold> Displacement nephograms at <italic>t</italic> &#x3d; 20&#xa0;s. <bold>(B)</bold> Displacement time history comparison of the slider.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g009.tif"/>
</fig>
<p>As shown in <xref ref-type="fig" rid="F8">Figure 8</xref>, landslides occur in the slope sites of upper soft rock and lower hard rock exposed to obliquely incident P waves with different incident angles. Landslide surfaces appear on the slope bodies, and the blocks slide down along the landslide surfaces. With increasing incident angle, the maximum plastic strain value increases, and the plastic zone of the slope site enlarges. Meanwhile, the landslide surface starts earlier at larger incident angles. <xref ref-type="fig" rid="F9">Figure 9A</xref> shows that the volume and sliding displacement of the landslide mass increase with an increasing incident angle, and the sliding slope surface changes from steep to gentle. This means that the landslide hazard becomes higher with an increase in the oblique incidence angle of the P wave.</p>
<p>To better reflect the displacement of each slider at different incident angles, <xref ref-type="fig" rid="F9">Figure 9B</xref> shows the displacement time history comparison of a slider at such angles. As shown in <xref ref-type="fig" rid="F9">Figure 9B</xref>, the absolute value of the horizontal average displacement of the slider increases with an increase of the incident angle, while the absolute value of the vertical average displacement decreases. The resultant displacement of the sliding block first decreases and then increases with an increase of the incidence angle, and the corresponding deformation of the slope site gradually changes from elastic to plastic deformation.</p>
</sec>
<sec id="s4-2">
<title>4.2 Positions of the soft-hard rock interface</title>
<p>This subsection analyzes the effect of the soft and hard rock interface positions on the dynamic response of the slope site by adopting the slope site model of the upper soft rock and lower hard rock shown in <xref ref-type="fig" rid="F3">Figure 3</xref>. The distances between the soft&#x2013;hard rock interface and the slope top are <italic>h</italic> &#x3d; 30, 27, 24, and 21&#xa0;m, respectively. The smaller the distance between the soft-hard rock interface and the slope top, the thinner the soft rock overburden is above the hard rock. The site parameters are shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<p>P waves are obliquely incident at an angle of <italic>&#x3b1;</italic> &#x3d; 45&#xb0; from the lower left corner of the model. <xref ref-type="fig" rid="F10">Figure 10</xref> shows the plastic strain nephograms of the slope site under the action of obliquely incident P waves with different soft&#x2013;hard rock interface positions at different timings. <xref ref-type="fig" rid="F11">Figure 11A</xref> depicts the displacement nephograms of the slope site at <italic>t</italic> &#x3d; 20&#xa0;s. <xref ref-type="fig" rid="F10">Figures 10</xref>, <xref ref-type="fig" rid="F11">11A</xref> show that, for the slope site with upper soft rock and lower hard rock, the plastic zones first appear at the interfaces between soft rock and hard rock and then continuously develop upward to form landslide surfaces. A thicker soft-rock overburden tends to cause an earlier appearance of the plastic strain zone, faster sliding surface formation, a greater plastic zone and plastic deformation, and stronger sliding.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Plastic strain nephograms of the slope at different timings with different soft&#x2013;hard rock interface positions. <bold>(A)</bold> <italic>h</italic> &#x3d; 30&#xa0;m. <bold>(B)</bold> <italic>h</italic> &#x3d; 27&#xa0;m. <bold>(C)</bold> <italic>h</italic> &#x3d; 24&#xa0;m. <bold>(D)</bold> <italic>h</italic> &#x3d; 21&#xa0;m.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g010.tif"/>
</fig>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Displacement of the slope at different interface positions of soft and hard rock. <bold>(A)</bold> Displacement nephograms at <italic>t</italic> &#x3d; 20&#xa0;s. <bold>(B)</bold> Displacement time history comparison of the slider.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g011.tif"/>
</fig>
<p>In order to better reflect the influence of the soft&#x2013;hard rock interface position on the sliding displacement of the slider, <xref ref-type="fig" rid="F11">Figure 11B</xref> describes the displacement time history comparison results of the slider at different positions of the soft&#x2013;hard rock interface. Thus, the absolute value of the horizontal and vertical average displacement and the resultant displacement of the slider all increase as the distance between the soft&#x2013;hard rock interface to the top of the slope changes from near to far.</p>
<p>According to the comprehensive analysis in <xref ref-type="fig" rid="F10">Figures 10, 11</xref>, on slope sites of upper soft rock and lower hard rock, the smaller the thickness of the soft rock layer, the safer the site will be. When the interface between soft and hard rock is 21&#xa0;m from the slope top, the plastic strain zone appears late, and the plastic strain and sliding displacement of the slider are small. The thickness of soft rock is small enough to prevent landslides.</p>
</sec>
<sec id="s4-3">
<title>4.3 Impedance ratio</title>
<p>This subsection analyzes the dynamic responses of the layered slope site at different impedance ratios (<italic>&#x3b7;</italic> &#x3d; 0.2, 0.4, 0.6, 0.8, and 1.0) of the upper and lower layers. The size of the slope site finite element model is the same as in <xref ref-type="sec" rid="s4-1">Section 4.1</xref>; the site parameters are shown in <xref ref-type="table" rid="T2">Table 2</xref>. P waves are obliquely incident at an angle of <italic>&#x3b1;</italic> &#x3d; 15&#xb0; from the lower left corner of the model.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Material constants of the slope site.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Case</th>
<th align="center">Impedance ratio <italic>a</italic>
</th>
<th align="center">Layer</th>
<th align="center">Density <italic>&#x3c1;</italic> (kg/m<sup>3</sup>)</th>
<th align="center">Elastic modulus <italic>E</italic> (Pa)</th>
<th align="center">Poisson&#x2019;s ratio <italic>&#x3c5;</italic>
</th>
<th align="center">Internal friction angle <italic>&#x3c6;</italic> (&#xb0;)</th>
<th align="center">Cohesion <italic>c</italic> (Pa)</th>
<th align="center">
<italic>C</italic>
<sub>
<italic>P</italic>
</sub> (m/s)</th>
<th align="center">
<italic>C</italic>
<sub>
<italic>SV</italic>
</sub> (m/s)</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="2" align="center">Case 1</td>
<td rowspan="2" align="center">0.2</td>
<td align="center">Upper</td>
<td align="center">1700</td>
<td align="center">1.2e9</td>
<td align="center">0.42</td>
<td align="center">21</td>
<td align="center">0.08e6</td>
<td align="char" char=".">1342</td>
<td align="char" char=".">499</td>
</tr>
<tr>
<td valign="top" align="center">Lower</td>
<td valign="top" align="center">2500</td>
<td valign="top" align="center">20.2e9</td>
<td valign="top" align="center">0.25</td>
<td align="center">50</td>
<td align="center">1.5e6</td>
<td align="char" char=".">3114</td>
<td align="char" char=".">1798</td>
</tr>
<tr>
<td rowspan="2" align="center">Case 2</td>
<td rowspan="2" align="center">0.4</td>
<td valign="top" align="center">Upper</td>
<td valign="top" align="center">1700</td>
<td align="center">1.2e9</td>
<td align="center">0.42</td>
<td align="center">21</td>
<td align="center">0.08e6</td>
<td align="char" char=".">1342</td>
<td align="char" char=".">499</td>
</tr>
<tr>
<td align="center">Lower</td>
<td valign="top" align="center">2200</td>
<td valign="top" align="center">5.8e9</td>
<td valign="top" align="center">0.31</td>
<td align="center">37</td>
<td align="center">0.6e6</td>
<td align="char" char=".">1912</td>
<td align="char" char=".">1003</td>
</tr>
<tr>
<td rowspan="2" align="center">Case 3</td>
<td rowspan="2" align="center">0.6</td>
<td valign="top" align="center">Upper</td>
<td valign="top" align="center">1700</td>
<td align="center">1.2e9</td>
<td align="center">0.42</td>
<td align="center">21</td>
<td align="center">0.08e6</td>
<td align="char" char=".">1342</td>
<td align="char" char=".">499</td>
</tr>
<tr>
<td valign="top" align="center">Lower</td>
<td valign="top" align="center">2100</td>
<td valign="top" align="center">2.7e9</td>
<td valign="top" align="center">0.32</td>
<td align="center">31</td>
<td align="center">0.35e6</td>
<td align="char" char=".">1356</td>
<td align="char" char=".">698</td>
</tr>
<tr>
<td rowspan="2" align="center">Case 4</td>
<td rowspan="2" align="center">0.8</td>
<td valign="top" align="center">Upper</td>
<td valign="top" align="center">1700</td>
<td align="center">1.2e9</td>
<td align="center">0.42</td>
<td align="center">21</td>
<td align="center">0.08e6</td>
<td align="char" char=".">1342</td>
<td align="char" char=".">499</td>
</tr>
<tr>
<td valign="top" align="center">Lower</td>
<td valign="top" align="center">1900</td>
<td valign="top" align="center">1.68e9</td>
<td valign="top" align="center">0.37</td>
<td align="center">25</td>
<td align="center">0.16e6</td>
<td align="char" char=".">1251</td>
<td align="char" char=".">568</td>
</tr>
<tr>
<td rowspan="2" align="center">Case 5</td>
<td rowspan="2" align="center">1.0</td>
<td valign="top" align="center">Upper</td>
<td valign="top" align="center">1700</td>
<td align="center">1.2e9</td>
<td align="center">0.42</td>
<td align="center">21</td>
<td align="center">0.08e6</td>
<td align="char" char=".">1342</td>
<td align="char" char=".">499</td>
</tr>
<tr>
<td valign="top" align="center">Lower</td>
<td valign="top" align="center">1700</td>
<td valign="top" align="center">1.2e9</td>
<td valign="top" align="center">0.42</td>
<td align="center">21</td>
<td align="center">0.08e6</td>
<td align="char" char=".">1342</td>
<td align="char" char=".">499</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Plastic strain nephograms of the slope with different impedance ratios are given in <xref ref-type="fig" rid="F12">Figure 12</xref>. The greater the impedance ratio between the upper and lower media, the later the plastic strain zone appears, the slower the development speed, and the later the landslide surface will occur. The maximum plastic strain of the slope site decreases with increasing impedance ratio between the upper and lower media.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Plastic strain nephograms of the slope at different timing with different impedance ratios. <bold>(A)</bold> <italic>&#x3b7;</italic> &#x3d; 0.2. <bold>(B)</bold> <italic>&#x3b7;</italic> &#x3d; 0.4. <bold>(C)</bold> <italic>&#x3b7;</italic> &#x3d; 0.6. <bold>(D)</bold> <italic>&#x3b7;</italic> &#x3d; 0.8. <bold>(E)</bold> <italic>&#x3b7;</italic> &#x3d; 1.0.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g012.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F13">Figure 13A</xref> illustrates the displacement nephograms of the slope site with different impedance ratios at <italic>t</italic> &#x3d; 20&#xa0;s. The impedance ratio has no effect on the position of the landslide surface and the volume of the sliding block but has a significant effect on the sliding displacement of the sliding block.</p>
<fig id="F13" position="float">
<label>FIGURE 13</label>
<caption>
<p>Displacement of the slope at different impedance ratios. <bold>(A)</bold> Displacement nephograms at <italic>t</italic> &#x3d; 20&#xa0;s. <bold>(B)</bold> Displacement time history comparison of the slider.</p>
</caption>
<graphic xlink:href="feart-10-1057316-g013.tif"/>
</fig>
<p>To better reflect the influence of the impedance ratio between the upper and lower media on the sliding displacement of the slider, <xref ref-type="fig" rid="F13">Figure 13B</xref> describes the displacement time history comparison results of the slider at different impedance ratios. It can be concluded that the absolute value of the horizontal and vertical average displacement and the resultant displacement of the slider have the same change trend under different impedance ratios, and that they all decrease with the increase in impedance ratio. With this increase, the landslide hazard of the slope site decreases significantly.</p>
</sec>
</sec>
<sec id="s5">
<title>5 Summary and conclusion</title>
<p>Landslide seismic disasters have become a focus of geological hazard research because of their wide distribution, frequency, and great harm. Due to the topographic effect of slope sites, the seismic input method of regular sites is no longer applicable. This paper combines viscoelastic boundary with a seismic input method for slope sites under the action of obliquely incident P waves. Application is made with MATLAB programming software. The finite element model of the slope site is established in the application software ABAQUS, and the seismic mechanical behavior of the slope site is simulated by the Mohr-Coulomb model. Compared with the calculation results of the extended model, the established seismic input method is verified as correct. The established seismic input method is applied to analyze the parameters of landslide disasters in slope sites with upper soft rock and lower hard rock.</p>
<p>Under all calculation working conditions, the landslide hazard occurrence process of the slope site is basically the same, but the severity of the landslide hazard is different. Under the action of obliquely incident P waves, a small range of plastic deformation in upper soft rock slopes and lower hard rock slopes first occurs at the soft&#x2013;hard rock interface. Next, the plastic zone expands upward to form a small-scale plastic zone, and the maximum plastic strain value increases continuously; the input velocity time history is then close to the peak. Thereafter, the plastic zone expands upward. When the input seismic wave reaches peak value, an obvious slip surface is formed. The slope body slides downward along this, and the deformation of the slope top is obvious. After the peak value is reached, there is no obvious change in the size of the plastic zone, but the maximum plastic strain value still increases in a small range over time. The details are as follows.<list list-type="simple">
<list-item>
<p>(1) For the slope site with upper soft rock and lower hard rock, the plastic zone and plastic strain value increase as the incident angle of the P wave changes from 0&#xb0; to 45&#xb0;. The larger the incident angle of the P wave, the earlier the landslide surface will form, the larger the landslide volume and sliding displacement, and the more serious the landslide disaster will be.</p>
</list-item>
<list-item>
<p>(2) For slope sites with upper soft rock and lower hard rock, the plastic zones first appear at the soft&#x2013;hard rock interface and then continuously develop upward to form the landslide surfaces. The farther the soft&#x2013;hard interface is from the slope top, the earlier the plastic strain zone will appear, and the faster the slip surface forms. With the thickening of the soft rock overburden, the plastic zone area and plastic deformation increase, and the sliding displacement of the slider is more obvious.</p>
</list-item>
<list-item>
<p>(3) The greater the impedance ratio between the upper and lower media of the layered slope site, the later the plastic strain zone will appear, and the slower the plastic zone develops. At different impedance ratios, the location of the landslide surface and the volume of the sliding block are basically the same. The maximum plastic strain and sliding displacement of the sliding block decrease with the increase of impedance ratio.</p>
</list-item>
</list>
</p>
<p>Some limitations should be noted with this paper. As the type of seismic failure of hard rock slopes is different from that studied in this paper, the proposed method is not suitable for such slopes. These conclusions and findings are based on the simulation results of the calculation model of the soft rock slope site and slope site with upper soft rock and lower hard rock. Due to the complexity of the problems, the validity of these conclusions and findings for calculating a model of hard rock slope sites or soil slope sites need much more study. Moreover, the simulation of the material of the slope site is based on the Mohr-Coulomb model without considering the existence of cracks and structural planes in rock mass materials. In general, the limits mentioned above should be addressed in future work.</p>
<p>Finally, the proposed method only discusses seismic landslides under the action of P waves, and further research is needed on seismic landslides under the action of shear waves and Rayleigh waves. Follow-up study regarding these issues is needed.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>HJ: data curation, writing&#x2014;original draft preparation, and program. JX: conceptualization, writing-reviewing, and fund support. MZ: guide the establishment of method and writing-reviewing. JH: software and assist in the establishment of the method. XD: resources, supervision, writing&#x2014;review and editing. JW: visualization and data curation.</p>
</sec>
<sec id="s8">
<title>Funding</title>
<p>This work was supported by Postdoctoral Science Foundation of China (2022M722929), the Special Fund of the Institute of Geophysics, China Earthquake Administration (DQJB20B23), and Beijing Natural Science Foundation Program (JQ19029). The support is gratefully acknowledged.</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Assimaki</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Gazetas</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Kausel</surname>
<given-names>E.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Effects of local soil conditions on the topographic aggravation of seismic motion: Parametric investigation and recorded field evidence from the 1999 athens earthquake</article-title>. <source>Bull. Seismol. Soc. Am.</source> <volume>95</volume> (<issue>3</issue>), <fpage>1059</fpage>&#x2013;<lpage>1089</lpage>. <pub-id pub-id-type="doi">10.1785/0120040055</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Astley</surname>
<given-names>R. J.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Infinite elements for wave problems: A review of current formulations and an assessment of accuracy</article-title>. <source>Int. J. Numer. Methods Eng.</source> <volume>49</volume> (<issue>7</issue>), <fpage>951</fpage>&#x2013;<lpage>976</lpage>. <pub-id pub-id-type="doi">10.1002/1097-0207(20001110)49:7&#x3c;951::aid-nme989&#x3e;3.0.co;2-t</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Athanasopoulos</surname>
<given-names>G. A.</given-names>
</name>
<name>
<surname>Pelekis</surname>
<given-names>P. C.</given-names>
</name>
<name>
<surname>Leonidou</surname>
<given-names>E. A.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Effects of surface topography on seismic ground response in the Egion (Greece) 15 June 1995 earthquake</article-title>. <source>Soil Dyn. Earthq. Eng.</source> <volume>18</volume> (<issue>2</issue>), <fpage>135</fpage>&#x2013;<lpage>149</lpage>. <pub-id pub-id-type="doi">10.1016/s0267-7261(98)00041-4</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bazyar</surname>
<given-names>M. H.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Analysis of transient wave scattering and its applications to site response analysis using the scaled boundary finite-element method</article-title>. <source>Soil Dyn. Earthq. Eng.</source> <volume>98</volume>, <fpage>191</fpage>&#x2013;<lpage>205</lpage>. <pub-id pub-id-type="doi">10.1016/j.soildyn.2017.04.010</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bettess</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>1977</year>). <article-title>Infinite elements</article-title>. <source>Int. J. Numer. Methods Eng.</source> <volume>11</volume> (<issue>1</issue>), <fpage>53</fpage>&#x2013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1002/nme.1620110107</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bird</surname>
<given-names>J. F.</given-names>
</name>
<name>
<surname>Bommer</surname>
<given-names>J. J.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Earthquake losses due to ground failure</article-title>. <source>Eng. Geol.</source> <volume>75</volume> (<issue>2</issue>), <fpage>147</fpage>&#x2013;<lpage>179</lpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2004.05.006</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bouckovalas</surname>
<given-names>G. D.</given-names>
</name>
<name>
<surname>Papadimitriou</surname>
<given-names>A. G.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Numerical evaluation of slope topography effects on seismic ground motion</article-title>. <source>Soil Dyn. Earthq. Eng.</source> <volume>25</volume> (<issue>7-10</issue>), <fpage>547</fpage>&#x2013;<lpage>558</lpage>. <pub-id pub-id-type="doi">10.1016/j.soildyn.2004.11.008</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cavallin</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Slejko</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>1986</year>). <article-title>Statistical approach to landslide and soil liquefaction hazards in seismic areas</article-title>. <source>Geol. Appl. E Idrogeol.</source> <volume>21</volume> (<issue>2</issue>), <fpage>231</fpage>&#x2013;<lpage>236</lpage>.</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chigira</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yagi</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Geological and geomorphological characteristics of landslides triggered by the 2004 Mid Niigta prefecture earthquake in Japan</article-title>. <source>Eng. Geol.</source> <volume>82</volume> (<issue>4</issue>), <fpage>202</fpage>&#x2013;<lpage>221</lpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2005.10.006</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Clayton</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Engquist</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>1977</year>). <article-title>Absorbing boundary conditions for acoustic and elastic wave equations</article-title>. <source>Bull. Seismol. Soc. Am.</source> <volume>67</volume> (<issue>6</issue>), <fpage>1529</fpage>&#x2013;<lpage>1540</lpage>. <pub-id pub-id-type="doi">10.1785/bssa0670061529</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deeks</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Randolph</surname>
<given-names>M. F.</given-names>
</name>
</person-group> (<year>1994</year>). <article-title>Axisymmetric time-domain transmitting boundaries</article-title>. <source>J. Eng. Mech.</source> <volume>120</volume> (<issue>1</issue>), <fpage>25</fpage>&#x2013;<lpage>42</lpage>. <pub-id pub-id-type="doi">10.1061/(asce)0733-9399(1994)120:1(25)</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J. T.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>A stress artificial boundary in FEA for near-field wave problem</article-title>. <source>Chin. J. Theor. Appl. Mech.</source> <volume>38</volume>, <fpage>49</fpage>&#x2013;<lpage>56</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gischig</surname>
<given-names>V. S.</given-names>
</name>
<name>
<surname>Eberhardt</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Moore</surname>
<given-names>J. R.</given-names>
</name>
<name>
<surname>Hungr</surname>
<given-names>O.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>On the seismic response of deep-seated rock slope instabilities&#x2014;insights from numerical modeling</article-title>. <source>Eng. Geol.</source> <volume>193</volume>, <fpage>1</fpage>&#x2013;<lpage>18</lpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2015.04.003</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>J. Q.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Impact of incident angles of P waves on the dynamic responses of long lined tunnels</article-title>. <source>Earthq. Eng. Struct. Dyn.</source> <volume>45</volume> (<issue>15</issue>), <fpage>2435</fpage>&#x2013;<lpage>2454</lpage>. <pub-id pub-id-type="doi">10.1002/eqe.2772</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>J. Q.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2017b</year>). <article-title>Impact of incident angles of earthquake shear (S) waves on 3-D non-linear seismic responses of long lined tunnels</article-title>. <source>Eng. Geol.</source> <volume>222</volume>, <fpage>168</fpage>&#x2013;<lpage>185</lpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2017.03.017</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>J. Q.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
</person-group> (<year>2017a</year>). <article-title>Non-linear seismic responses of tunnels within normal fault ground under obliquely incident P waves</article-title>. <source>Tunn. Undergr. Space Technol.</source> <volume>61</volume>, <fpage>26</fpage>&#x2013;<lpage>39</lpage>. <pub-id pub-id-type="doi">10.1016/j.tust.2016.09.006</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jibson</surname>
<given-names>R. W.</given-names>
</name>
<name>
<surname>Harp</surname>
<given-names>E. L.</given-names>
</name>
<name>
<surname>Michael</surname>
<given-names>J. A.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>A method for producing digital probabilistic seismic landslide hazard maps</article-title>. <source>Eng. Geol.</source> <volume>58</volume> (<issue>3</issue>), <fpage>271</fpage>&#x2013;<lpage>289</lpage>. <pub-id pub-id-type="doi">10.1016/s0013-7952(00)00039-9</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>JTG/T D70-2010</collab>
</person-group> (<year>2010</year>). <article-title>Guidelines for Design of Highway Tunnel</article-title>. <source>The People&#x2019;s Republic of China: Industry Recommended Standards of the People&#x2019;s Republic of China</source>, <fpage>35</fpage>&#x2013;<lpage>48</lpage> <comment>(In Chinese)</comment>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Keeper</surname>
<given-names>D. K.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>Landslides caused by earthquakes</article-title>. <source>Geol. Soc. Am. Bull.</source> <volume>95</volume> (<issue>4</issue>), <fpage>406</fpage>. <pub-id pub-id-type="doi">10.1130/0016-7606(1984)95&#x3c;406:lcbe&#x3e;2.0.co;2</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lenti</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Martino</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>The interaction of seismic waves with step-like slopes and its influence on landslide movements</article-title>. <source>Eng. Geol.</source> <volume>126</volume>, <fpage>19</fpage>&#x2013;<lpage>36</lpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2011.12.002</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>A novel approach to probabilistic seismic landslide hazard mapping using Monte Carlo simulations</article-title>. <source>Eng. Geol.</source> <volume>301</volume>, <fpage>106616</fpage>. <pub-id pub-id-type="doi">10.1016/j.enggeo.2022.106616</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>H. F.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Accurate H-shaped absorbing boundary condition in frequency domain for scalar wave propagation in layered half-space</article-title>. <source>Int. J. Numer. Methods Eng.</source> <volume>121</volume>, <fpage>4268</fpage>&#x2013;<lpage>4291</lpage>. <pub-id pub-id-type="doi">10.1002/nme.6424</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liao</surname>
<given-names>Z. P.</given-names>
</name>
<name>
<surname>Wong</surname>
<given-names>H. L.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>A transmitting boundary for the numerical simulation of elastic wave propagation</article-title>. <source>Int. J. Soil Dyn. Earthq. Eng.</source> <volume>3</volume> (<issue>4</issue>), <fpage>174</fpage>&#x2013;<lpage>183</lpage>. <pub-id pub-id-type="doi">10.1016/0261-7277(84)90033-0</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>Y. X.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z. Y.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>3D viscous-spring artificial boundary in time domain</article-title>. <source>Earthq. Engin. Engin. Vib.</source> <volume>5</volume> (<issue>1</issue>), <fpage>93</fpage>&#x2013;<lpage>102</lpage>. <pub-id pub-id-type="doi">10.1007/s11803-006-0585-2</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>Y. D.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>A direct method for analysis of dynamic soil-structure interaction based on interface idea</article-title>. <source>Dev. Geotechnical Eng.</source> <volume>83</volume> (<issue>3</issue>), <fpage>261</fpage>&#x2013;<lpage>276</lpage>. <pub-id pub-id-type="doi">10.1016/s0165-1250(98)80018-7</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>A 1D Time-domain method for inplane wave motion of free field in layered media</article-title>. <source>Chin. J. Eng. Mech.</source> <volume>24</volume> (<issue>7</issue>), <fpage>16</fpage>&#x2013;<lpage>22</lpage>.</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lysmer</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Kuhlemeyer</surname>
<given-names>R. L.</given-names>
</name>
</person-group> (<year>1969</year>). <article-title>Finite dynamic model for infinite media</article-title>. <source>J. Engrg. Mech. Div.</source> <volume>95</volume> (<issue>4</issue>), <fpage>859</fpage>&#x2013;<lpage>877</lpage>. <pub-id pub-id-type="doi">10.1061/jmcea3.0001144</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Maleki</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Khodakarami</surname>
<given-names>M. I.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Feasibility analysis of using MetaSoil scatterers on the attenuation of seismic amplification in a site with triangular hill due to SV-waves</article-title>. <source>Soil Dyn. Earthq. Eng.</source> <volume>100</volume>, <fpage>169</fpage>&#x2013;<lpage>182</lpage>. <pub-id pub-id-type="doi">10.1016/j.soildyn.2017.05.036</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nakamura</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Two-dimensional energy transmitting boundary in the time domain</article-title>. <source>Earthquakes Struct.</source> <volume>3</volume> (<issue>2</issue>), <fpage>97</fpage>&#x2013;<lpage>115</lpage>. <pub-id pub-id-type="doi">10.12989/eas.2012.3.2.097</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Owen</surname>
<given-names>L. A.</given-names>
</name>
<name>
<surname>Kamp</surname>
<given-names>U.</given-names>
</name>
<name>
<surname>Khattak</surname>
<given-names>G. A.</given-names>
</name>
<name>
<surname>Harp</surname>
<given-names>E. L.</given-names>
</name>
<name>
<surname>Keefer</surname>
<given-names>D. K.</given-names>
</name>
<name>
<surname>Bauer</surname>
<given-names>M. A.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Landslides triggered by the 8 october 2005 kashmir earthquake</article-title>. <source>Geomorphology</source> <volume>94</volume> (<issue>1-2</issue>), <fpage>1</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1016/j.geomorph.2007.04.007</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pareek</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Arora</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Arora</surname>
<given-names>M. K.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Impact of seismic factors on landslide susceptibility zonation: A case study in part of Indian himalayas</article-title>. <source>Landslides</source> <volume>7</volume>, <fpage>191</fpage>&#x2013;<lpage>201</lpage>. <pub-id pub-id-type="doi">10.1007/s10346-009-0192-1</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pelekis</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Batilas</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Pefani</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Vlachakis</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Athanasopoulos</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Surface topography and site stratigraphy effects on the seismic response of a slope in the Achaia-Ilia (Greece) 2008 Mw6. 4 earthquake</article-title>. <source>Soil Dyn. Earthq. Eng.</source> <volume>100</volume>, <fpage>538</fpage>&#x2013;<lpage>554</lpage>. <pub-id pub-id-type="doi">10.1016/j.soildyn.2017.05.038</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Poursartip</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Fathi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Kallivokas</surname>
<given-names>L. F.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Seismic wave amplification by topographic features: A parametric study</article-title>. <source>Soil Dyn. Earthq. Eng.</source> <volume>92</volume>, <fpage>503</fpage>&#x2013;<lpage>527</lpage>. <pub-id pub-id-type="doi">10.1016/j.soildyn.2016.10.031</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Prestininzi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Romeo</surname>
<given-names>R. W.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Earthquake-induced ground failures in Italy</article-title>. <source>Eng. Geol.</source> <volume>58</volume> (<issue>3-4</issue>), <fpage>387</fpage>&#x2013;<lpage>397</lpage>. <pub-id pub-id-type="doi">10.1016/s0013-7952(00)00044-2</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rizzitano</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Cascone</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Biondi</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Coupling of topographic and stratigraphic effects on seismic response of slopes through 2D linear and equivalent linear analyses</article-title>. <source>Soil Dyn. Earthq. Eng.</source> <volume>67</volume>, <fpage>66</fpage>&#x2013;<lpage>84</lpage>. <pub-id pub-id-type="doi">10.1016/j.soildyn.2014.09.003</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sato</surname>
<given-names>P. H.</given-names>
</name>
<name>
<surname>Hasegawa</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Fujiwara</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Tobita</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Koarai</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Une</surname>
<given-names>H.</given-names>
</name>
<etal/>
</person-group> (<year>2007</year>). <article-title>Interpretation of landslide distribution triggered by the 2005 Northern Pakistan earthquake using SPOT 5 imagery</article-title>. <source>Landslides</source> <volume>4</volume> (<issue>2</issue>), <fpage>113</fpage>&#x2013;<lpage>122</lpage>. <pub-id pub-id-type="doi">10.1007/s10346-006-0069-5</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Semblat</surname>
<given-names>J. F.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Modeling seismic wave propagation and amplification in 1D/2D/3D linear and nonlinear unbounded media</article-title>. <source>Int. J. Geomech.</source> <volume>11</volume> (<issue>6</issue>), <fpage>440</fpage>&#x2013;<lpage>448</lpage>. <pub-id pub-id-type="doi">10.1061/(asce)gm.1943-5622.0000023</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>Z. D.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J. T.</given-names>
</name>
<name>
<surname>Zhong</surname>
<given-names>Z. L.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Response spectrum method for seismic soil-structure interaction analysis of underground structure</article-title>. <source>Bull. Earthq. Eng.</source> <volume>17</volume>, <fpage>5339</fpage>&#x2013;<lpage>5363</lpage>. <pub-id pub-id-type="doi">10.1007/s10518-019-00673-6</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>Z. D.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L. T.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>J. Q.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Obliquely incident earthquake input for soil-structure interaction in layered half space</article-title>. <source>Earthquakes Struct.</source> <volume>13</volume> (<issue>06</issue>), <fpage>573</fpage>&#x2013;<lpage>588</lpage>.</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yin</surname>
<given-names>H. Q.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J. B.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>1D finite element artificial boundary method for layered half space site response from obliquely incident earthquake</article-title>. <source>Earthquakes Struct.</source> <volume>1</volume> (<issue>9</issue>), <fpage>173</fpage>&#x2013;<lpage>194</lpage>. <pub-id pub-id-type="doi">10.12989/eas.2015.9.1.173</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>