An Arc-Consistent Viscous-Spring Artificial Boundary for Numerical Analysis of Seismic Response of Underground Structures

Ye, Dan; Yin, Shangzhi; Wang, Yihong; Zuo, Taian

doi:https://doi.org/10.1155/2022/2791492

Shock and Vibration

On this page

Abstract Introduction Results and Discussion Conclusions Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Special Issue

Seismic Performance, Damage Mechanism and Seismic Design for Tunnel Engineering

View this Special Issue

Research Article | Open Access

Volume 2022 | Article ID 2791492 | https://doi.org/10.1155/2022/2791492

An Arc-Consistent Viscous-Spring Artificial Boundary for Numerical Analysis of Seismic Response of Underground Structures

Dan Ye,¹Shangzhi Yin,²Yihong Wang,³and Taian Zuo⁴

Academic Editor: Xuepeng Zhang

Received21 Oct 2021

Accepted20 Dec 2021

Published10 Jan 2022

Abstract

A new arc-consistent viscous-spring artificial boundary (ACVAB) was proposed by changing a traditional flat artificial boundary based on the theory of viscous-spring artificial boundaries. Through examples, the concept underpinning the establishment and specific setting of the boundary in the finite element software were described. Through comparison with other commonly used artificial boundaries in an example for near-field wave analysis using the two-dimensional (2D) half-space model, the reliability of the ACVAB was verified. Furthermore, the ACVAB was used in the numerical analysis of the effects of an earthquake on underground structures. The results were compared with the shaking table test results on underground structures. On this basis, the applicability of the ACVAB to a numerical model of seismic response of underground structures was evaluated. The results show that the boundary is superior to common viscous-spring boundaries in terms of accuracy and stability, and therefore, it can be used to evaluate radiation damping effects of seismic response of underground structures and is easier to use.

1. Introduction

Given increasing urbanisation, underground rail transit has developed rapidly. Numerical simulation, as one of the main methods to study the seismic performance of underground structures at present, can be used to reproduce dynamic characteristics of underground structures under seismic load. The infinite domain (infinite foundation) can be simulated by introducing an artificial boundary [1–8], so the artificial boundary theory is one of the important bases for seismic analysis and numerical model study of underground structures. In the numerical simulation of near-field wave motion, artificial boundary conditions directly affect the accuracy of the calculations and a reasonable artificial boundary can simulate the propagation process of waves in an infinite domain medium. An ideal artificial boundary should be stable, accurate, efficient, and easy to realise, but there are often contradictions among the above characteristics.

The generalised artificial boundaries are classified into global artificial boundaries and local artificial boundaries. Global artificial boundaries are mainly characterised by coupled motion of all boundary nodes in space and time and satisfaction of all field equations and mathematical and physical conditions in the infinite domain, which can accurately simulate the infinite domain [9–11]. In the analysis of large-scale complex wave motion, especially nonlinear generalised structures, global artificial boundaries impose an onerous computational burden. The boundary integral equation [12–14], thin-layer method [15–17], and exact Dirichlet-to-Neumann artificial boundaries [18–20] are categorised to global artificial boundary conditions. Local artificial boundaries allow approximate simulation of the infinite domain. Unlike global artificial boundaries, local artificial boundaries are mainly characterised by spatiotemporal decoupling [21], indicating that there is no need to solve simultaneous equations, which greatly reduce the amount of calculation required. At present, the main local artificial boundary conditions include a Sommerfeld boundary [22, 23], paraxial axial approximation artificial boundary (a Clayton–Engquist boundary) [24, 25], Higdon boundary [26–28], multi-transmitting boundary [29–32], viscous boundary [33, 34], and a viscous-spring boundary [35].

The artificial boundaries introduced above are mostly expressed in the form of differential equations of boundaries by an approximate treatment of wave equations. However, in the combination with the explicit finite element method, the discrete form is inconvenient to engineering design and there is no obvious advantage in its simulation accuracy. By referring to other studies [36–38], several boundaries, such as viscous boundary, viscous-spring boundary, consistent viscous-spring boundary elements, and remote artificial boundary, are discrete artificial boundary conditions established using approximate solutions to one-sided wave equations. They can be easily realised by the appropriate finite element analysis software.

Alterman et al. [36, 37] performed numerical simulation on near-field wave motion and compared it with concise analytical solutions of integral transformation based on the remote artificial boundary. The principle of the remote boundary is that for the problem of near-field wave motion, when the generalised structure is far enough from the artificial boundary, the artificial boundary problem can be ignored; that is, the artificial boundary can be set arbitrarily. This processing method of the remote artificial boundary has been still used to provide accurate solutions hitherto. However, this method runs contrary to the goal of high efficiency and is invalid in current engineering design and calculation. The viscous boundary [33, 34] was first studied by Lysmer et al., whose working principle [38, 39] is derived based on the assumption that there is no energy reflection at the boundary. It is found from the research that the viscous boundary does not permit attenuation of waves in the medium and then does not consider the stiffness recovery of the medium, which is likely to raise problems of structural drift and low accuracy. Deek and Randolph (1994) [35] were the first to propose the expression of the two-dimensional (2D) plane strain problem of the viscous-spring boundary. To cope with the problems of structural drift and low accuracy in the use of viscous boundaries, Liu (2006) [38, 40, 41] proposed the concept of consistent viscous-spring artificial boundary elements based on the cylindrical wave equation. A viscous-spring artificial boundary with good elastic recovery performance was established using discrete spring dampers at the boundary, which overcame the problem of low-frequency instability, indeed resulting in good stability [40–42]. At present, the boundary has been widely utilised in finite element analysis software; however, when solving the dynamic problem of the complex 2D semi-infinite domain, the section of the viscous-spring boundary is not directly introduced into the most finite element analysis software. For this reason, researchers need to calculate the parameters of spring dampers and apply spring dampers one by one, requiring significant time in preprocessing. Moreover, in the derivation of the viscous-spring artificial boundary, the cylindrical wave assumption is adopted, but the boundary is simplified to a flat boundary, which is inconsistent with reality.

The concentrated viscous-spring dynamic artificial boundary requires onerous preprocessing and is of low accuracy and inconsistent with the cylindrical wave assumption. In view of these problems, a new artificial boundary, namely a 2D arc-consistent viscous-spring artificial boundary (ACVAB) element, was established based on the theory of consistent viscous-spring artificial boundary and viscous-spring artificial boundary element [38, 40]. Section 1 introduces the concept and realisation method of the ACVAB. Section 2 introduces the methods of setting equivalent stiffness and damping. In Section 3, the example in near-field wave analysis using the 2D half-space model and example in seismic response of underground structures are described. Section 4 discusses the results of examples and verifies the accuracy of the ACVAB and its applicability to the analysis of the seismic response of underground structures. The conclusions are summarised in Section 5.

2. ACVAB

2.1. Description of the Equivalent Viscous-Spring Artificial Boundary

In the 2D arc-consistent viscous-spring artificial boundary element, a global stiffness matrix and a damping matrix of the model are derived, as shown in formulas (1) and (2) [35]. In the finite element analysis software, the consistent viscous-spring artificial boundary is realised by extending a layer of the same type of elements along the normal to the outer surface of the established model area, fixing nodes on the outer element surface, and setting springs and dampers. Based on the equivalence principle of the matrix, the consistent stiffness and damping matrixes of the element are transformed into an equivalent stiffness matrix and an equivalent damping matrix of the boundary element, thus meeting the requirements of resilience and damping on the boundary (Figure 1).where K_BT, K_BN, C_BT, and C_BN represent the tangential stiffness, normal spring stiffness, tangential damping coefficient, and normal damping coefficient, respectively.

2.2. Description of the ACVAB

When processing the artificial boundary in a semi-infinite domain, to simplify the model, previous researchers have mostly adopted the flat artificial boundary (Figure 1). In the derivation of the viscous-spring artificial boundary, the cylindrical wave assumption is used, so it is inconsistent with actual situations when simplifying the boundary to be flat. Therefore, a new viscous-spring boundary processing method, that is, the ACVAB, was proposed in this study.

In the specific setting, firstly, a finite element model of an arc boundary of radius R is established taking the point of action of the scattered wave source as the centre, so that the distance from the wave source to the arc boundary remains unchanged. Secondly, a layer of ordinary elements can be extended along the normal direction on the boundary of the established finite element analysis model and the outer boundary is fixed. An appropriate material coefficient is given through calculation to replace the original discrete spring dampers. The calculation and setting of the coefficient are described in the following section, and the model is displayed in Figure 2.

This method does not need to calculate and arrange spring dampers one by one, which greatly simplifies the preprocessing work. Meanwhile, because the radius of the model is fixed, the error caused by the approximation of the coefficient of material property is avoided, and the accuracy of the dynamic finite element analysis is thus improved.

3. Setting of Equivalent Stiffness and Damping of the ACVAB

The artificial boundary introduced in the analysis of dynamic interaction between structures and foundation is generally deduced based on the assumption that there is no energy reflected at the boundary, such as the viscous boundary [33]. Based on the assumption that 2D scattered waves are cylindrical, Deeks [35] deduced the artificial boundary condition as follows:where represents the coordinates of the artificial boundary in a polar coordinate system; and are the viscous damper and linear spring applied on the artificial boundary , respectively; denotes the shear wave velocity; and G and separately denote the shear modulus and mass density of the medium.

Based on the cylindrical wave equation, Liu et al. [38, 40–42] established a 2D viscous-spring artificial boundary, in which the spring coefficient K_B and damping coefficient C_B of the equivalent physical system separately are expressed as follows:

Tangential boundary:

Normal boundary:where and refer to the normal and tangential stiffness of the spring; R represents the distance from the wave source to the artificial boundary point; and denote the S-wave and P-wave velocities; and and denote the parameters of the tangential and normal viscous-spring artificial boundaries, respectively.

As shown in Figure 2, the artificial boundary element B consists of continuously and uniformly distributed springs and dampers. As illustrated in Figure 3, the displacement shape function of the two-node artificial boundary element can be expressed by the linear interpolation function. The springs distributed in the normal direction of element B and the corresponding displacement shape function are shown in Figure 4.

Assuming u and represent the tangential and normal displacements of the artificial boundary, the displacements of nodes i and j of the artificial boundary element are expressed as follows:

As shown in Figure 4, the normal displacement shape function of the artificial boundary element is expressed as follows:

The tangential displacement shape function is the same as the normal displacement shape function.

The geometry matrix of the artificial boundary element corresponding to formula (3) is shown as follows:

The elastic stiffness matrix of springs distributed on the viscous-spring artificial boundary is expressed as follows:

The stiffness matrix of the viscous-spring artificial boundary element can be derived based on the principle of virtual work as follows:

By substituting formulas (8) and (9) into formula (10), the stiffness matrix of the consistent viscous-spring artificial boundary element is obtained:

Through the same steps to form the stiffness matrix of the consistent viscous-spring boundary element, the damping matrix of the 2D arc-consistent viscous-spring boundary element can be attained:

Formulas (11) and (12) demonstrate the stiffness matrix and damping matrix of the 2D arc-consistent viscous-spring artificial boundary element. When the arc viscous-spring boundary is segmented into infinite elements with the minimum thickness (extremely dense meshing), each element approximates a rectangular element. A four-node rectangular element is used to simulate the 2D ACVAB. As shown in Figure 5, the displacement shape function of the element is described as follows:In formula (12), l and h represent the side length and the extended thickness of the element, respectively; i, j, k, and l represent the nodes of the element, of which i and j indicate the nodes of the aforementioned consistent viscous-spring artificial boundary: because nodes l and k of the element are fixed, the explicit expression of the stiffness matrix of plane strain can be found, thus giving the equivalent stiffness matrix of the element corresponding to two nodes i and j：where represents the mass density of the equivalent element; and indicate the S-wave and P-wave velocities of the medium of the equivalent element.

When the thickness of the artificial boundary zone is , the stiffness matrix is approximated as follows:

By comparing formulas (11) and (15),according to the following definitions of S-wave and P-wave velocities,the equivalent shear modulus and elastic modulus of the boundary element are obtained as follows:where , , and denote the equivalent shear modulus, equivalent elastic modulus, and equivalent Poisson’s ratio of the equivalent viscous-spring boundary element, respectively; h, R, and G refer to the thickness of the equivalent element, distance from source of the wave to the artificial boundary point, and shear modulus of the medium, respectively.

The damping matrix of the equivalent boundary element is formed using damping proportional to the stiffness, and then, let

Then,where and separately represent the scale factors related to tangential and normal stiffness. By substituting formulas (11) and (12) into formula (19), the following formula is obtained:

To realise the boundary element in a convenient manner, the damping scale coefficient of the equivalent element is set as the average of the parameters in two directions.

4. Example Verification

4.1. Example for Near-Field Wave Analysis Using the 2D Half-Space Model

To verify the accuracy of the ACVAB proposed in this study, an example of the near-field wave analysis using the 2D half-space model was compared with other commonly used artificial boundaries. When modelling the ACVAB, a semicircular near-field observation area with a radius of 50 m was taken, a layer of elements with the unit thickness was extended around the model, and the nodes on the outer layer were fixed. The medium of elastic semi-infinite domain was selected for research, and the density of the medium, Young’s modulus, and Poisson’s ratio are 2000 kg/m³, 20 MPa, and 0.3, respectively. Considering the Rayleigh damping of the medium, a = 0.616 and β = 0.000312. The coefficient of material property of the arc-consistent viscous-spring element is calculated by referring to formulas (19) and (21). Half of a sinusoidal pulse at the frequency of 1 Hz was applied at the centre, and surface point A at the axis of symmetry and surface point B at a burial depth of 10 m were taken as points for observation. The 2D half-space model for near-field wave analysis is shown in Figure 6. To verify the effects of the ACVAB, multiple models including the concentrated viscous boundary, ACVAB, concentrated viscous-spring boundary, fixed boundary, and remote boundary were also established for comparisons.

4.2. Example for Seismic Response of Underground Structures

To assess the applicability of the ACVAB proposed in this research to the seismic analysis of underground structures, a numerical model for seismic response of underground structures was built based on the ACVAB. The model results were compared with data from shaking table tests [43], thus allowing the analyst to determine whether the boundary is applicable to the analysis of radiation damping effects in the seismic response of underground structures.

4.2.1. Description of Project Background

The geomorphic unit of Xi’an Feitian Road Station belongs to the second- and third-grade loess tableland. The lithologic characteristics of the site strata are summarised in Tables 1 and 2. The table is arranged from top to bottom according to the order of soil layers from shallow to deep. The underground subway station is a reinforced concrete structure with a total height of 14.01 m and a total width of 19.2 m. The longitudinal spacing of the centre pillar is 9 m. The cross section of the centre pillar measures 0.8 m × 1.2 m. The depth of soil above the roof is 3.459 m. The density of concrete ρ is 2.5 g/cm³, its modulus of elasticity E is 35 GPa, and Poisson’s ratio is 0.15. A typical cross section is shown in Figure 7.

The station is located in Chang’an District of Xi’an City. The site type in this area is class II, the basic seismic intensity is 8 degrees, and the characteristic period of seismic response spectrum is 0.4 s.

4.2.2. Description of Calculation Model

Due to the length of Feitian Station, the three-dimensional dynamic interaction system of loess subway underground structure was considered as a two-dimensional plane strain problem in the seismic analysis by the time history analysis method. The typical section of loess site and subway underground structure was numerically simulated, and the model was established as shown in Figure 8.

The ACVAB proposed in this study is adopted as the artificial boundary of calculation model. TAFT waves are chosen as the incident seismic waves, and the seismic fluctuation input problem is transformed into a wave source problem to realise the seismic wave input [41]. The plastic damage model of concrete in ABAQUS finite element analysis software was adopted to simulate the mechanical behavior of the prototype subway station concrete. The CPE4R element was used for on-site soil, and the CPE4 element was used to model the subway station structure. The initial stress on the soil was calculated using the geostatic module in ABAQUS. To simulate the in situ stress on the soil and how it affects the adjacent underground structure, the model states of excavation, support, construction of underground structure, and backfilling were established to simulate the construction process of this subway station, and the stress state in the soil after construction was taken as the initial stress field for the later dynamic analysis.

5. Results and Discussion

5.1. Comparison and Discussion of Results of the Example for Near-Field Wave Analysis Using the 2D Half-Space Model

To validate the accuracy of the 2D ACVAB, the displacements of observation points on each boundary model were compared and analysed. The remote boundary with complex modelling, low calculation efficiency, and optimal simulation effects was taken as the standard solution to seek for the artificial boundary with accuracy approaching the remote boundary and high modelling efficiency and calculation efficiency. As shown in Figure 9, due to reflection-induced interference, the fixed boundary model differs greatly from the standard solution obtained by the remote boundary and cannot reflect the actual propagation of waves in the semi-infinite domain medium. In the model of the concentrated viscous-spring boundary, the reflection of waves is reduced, but there is a significant downward drift in the middle and rear sections, so it is considered that its accuracy cannot meet the requirements of engineering practice. For the models of viscous boundary, ACVAB, and remote boundary, after inputting pulse waves, energy rapidly decreases, showing only small fluctuations and being almost completely absorbed by the boundary element, which meets the accuracy requirements imposed by engineering practice. Through comparison, the advantages and disadvantages of the commonly used artificial boundaries and the ACVAB proposed in this study are summarised in Table 3.

(a)

(b)

To compare the accuracy of artificial boundaries, Table 4 is obtained by dividing the displacement of observation points of each boundary model by the displacement of observation points on the remote boundary (standard solutions). After processing the data in Table 4, the factors of accuracy improvement using the proposed ACVAB compared with other commonly used artificial boundaries can be determined, as listed in Table 5. By combining this with the data in Tables 4 and 5, the accuracy of the model of the ACVAB increases with the increase in the burial depth. For surface observation point A (used for monitoring displacement), the accuracy of the ACVAB is improved by 0.6 to 3.4% compared with other commonly used artificial boundaries. In terms of observation point B at the deep part, its accuracy increases by 3.1 to 16.7% in comparison with other commonly used artificial boundaries. From the perspective of the mean value, its accuracy increases by 3.5 to 15.4% on average. In conclusion, it is considered that the calculation accuracy and stability of the ACVAB are superior to those of the viscous boundary and concentrated viscous-spring boundary. The accuracy is closer to that of the remote boundary, and the setting method is simpler. It is easier to realise this in the finite element analysis software.

5.2. Comparison and Discussion of Results of Example for Seismic Response of Underground Structures

The numerical simulation results of loess subway stations based on the ACVAB were compared with shaking table test results [43]. Under the action of seismic waves, the acceleration response time history of the foundation in the model and corresponding Fourier spectra in the numerical simulation and shaking table test were compared, as shown in Figure 10. The comparison of acceleration time history of the structure and corresponding Fourier spectra is shown in Figure 11.

As illustrated in Figure 10, the waveforms and amplitudes of acceleration response time history of the foundation in the model and composition characteristics of Fourier spectra recorded in the numerical simulation and shaking table test are coincided. At a low acceleration, the numerical simulation and test results are consistent in terms of waveforms and amplitudes of acceleration response time history of the foundation in the model and composition characteristics of Fourier spectra; however, as the acceleration increases, the simulated accelerations of the foundation show certain deviation from the test result, with the error within 15%. This result indicates that the proposed ACVAB can be used to simulate the propagation of seismic waves in the soil medium of the semi-infinite domain, which well solves the transmission of ground motion energy from near-field to far-field areas.

As demonstrated in Figure 11, the waveforms and amplitudes of acceleration response time history of the structure in the model and composition characteristics of Fourier spectra recorded in the numerical simulation and shaking table test are similar. When the acceleration is low, consistent waveforms and amplitudes of acceleration response time history of the structure in the model and composition characteristics of Fourier spectra are found in the numerical simulation and test data. As the acceleration rises, the numerical simulation shows a certain error (within 25%) from the test results in terms of acceleration, which is related to the model of structural materials and accuracy of processing method for contact between the structure and foundation.

In summary, the numerical simulation based on the ACVAB effectively reduces boundary effects and shows high calculation accuracy and efficiency in the simulation of propagation of seismic waves in soil mass. This boundary can be used to ascertain radiation damping effects in the calculation of seismic response of structures.

6. Conclusions

The concept of a new ACVAB was proposed based on the theory of the consistent viscous-spring artificial boundary and viscous-spring artificial boundary element. The establishment ideas of the artificial boundary and the implementation method in the finite element analysis software were described, and its applicability was studied by analysing numerical examples.(1)Considering the error caused by the approximation in the setting of the concentrated viscous-spring boundary and the cumbersome preprocessing work, the new artificial boundary, namely the ACVAB, was proposed in this research. The operation of this element model is simpler than that of the concentrated viscous-spring boundary, and its high accuracy and good stability were proven by example.(2)Through the example of near-field wave analysis using the 2D half-space model, the ACVAB was compared with other commonly used artificial boundaries. The results demonstrate that the boundary with high accuracy and simple setting method is easy to realise in the finite element analysis software and improves both modelling efficiency and calculation speed.(3)The numerical simulation results pertaining to loess and subways stations therein based on the ACVAB were compared with shaking table test results. The results show that the ACVAB can reduce boundary effects and simulate the propagation of seismic waves in a soil medium forming a semi-infinite domain, which favourably solves the transmission of ground motion energy from near-field to far-field areas. It can be used to analyse radiation damping effects in the calculation of seismic response of underground structures.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Science and Technology Research Project of Chongqing Education Commission, China (KJQN201901603), the Natural Science Foundation of Chongqing, China (cstc2021jcyj-msxmX0963), the Humanities and Social Sciences Research Project of Chongqing Education Commission, China (20SKGH255), and 2021 Aging Scientific Research Project of Chongqing University of Education, China (KY202123C).

References

Z. Gao, M. Zhao, X. Du, and Z. Zhong, “A generalized response spectrum method for seismic response analysis of underground structure combined with viscous-spring artificial boundary,” Soil Dynamics and Earthquake Engineering, vol. 140, Article ID 106451, 2021.
View at: Publisher Site | Google Scholar
C. Song and J. P. Wolf, “Consistent infinitesimal finite element cell method for incompressible unbounded medium,” Communications in Numerical Methods in Engineering, vol. 13, no. 1, pp. 21–32, 1997.
View at: Publisher Site | Google Scholar
Y. Li, M. Zhao, C.-s. Xu, X.-l. Du, and Z. Li, “Earthquake input for finite element analysis of soil-structure interaction on rigid bedrock,” Tunnelling and Underground Space Technology, vol. 79, pp. 250–262, 2018.
View at: Publisher Site | Google Scholar
J. Song and S. C. Li, “Numerical analysis of viscous-spring artificial boundary element for underground explosion,” Applied Mechanics and Materials, vol. 353-356, pp. 3345–3349, 2013.
View at: Publisher Site | Google Scholar
G. Mangiante and S. Charles, “Absorbing boundary conditions for acoustic waves and Huygens principle,” Journal of the Acoustical Society of America, vol. 103, no. 5, p. 2971, 1998.
View at: Publisher Site | Google Scholar
Z. Ba, Q. Sang, and J. Liang, “Seismic analysis of a lined tunnel in a multi-layered TI saturated half-space due to qP1-and qSV-waves,” Tunnelling and Underground Space Technology, vol. 119, Article ID 104248, 2022.
View at: Publisher Site | Google Scholar
R. Zhang, D. Ye, J. Zhou, and D. Quan, “Seismic analysis method for underground structure in loess area based on the modified displacement-based method,” Applied Sciences, vol. 11, no. 23, p. 11245, 2021.
View at: Publisher Site | Google Scholar
D. Ye and S. Yin, “An arc consistent viscous-spring artificial boundary for numerical analysis of seismic response of underground structures,” Scientific Reports, 2021, https://www.researchsquare.com/article/rs-1049487/v1.
View at: Publisher Site | Google Scholar
Z.-h. Teng, “Exact boundary conditions for the initial value problem of convex conservation laws,” Journal of Computational Physics, vol. 229, no. 10, pp. 3792–3801, 2010.
View at: Publisher Site | Google Scholar
V. S. Ryaben'kii, S. V. Tsynkov, and V. I. Turchaninov, “Global discrete artificial boundary conditions for time-dependent wave propagation,” Journal of Computational Physics, vol. 174, no. 2, pp. 712–758, 2001.
View at: Publisher Site | Google Scholar
T. Liu, C. Zhao, C. Zhang, and F. Jin, “A semi-analytical artificial boundary for time-dependent elastic wave propagation in two-dimensional homogeneous half space,” Soil Dynamics and Earthquake Engineering, vol. 30, no. 11, pp. 1352–1360, 2010.
View at: Publisher Site | Google Scholar
G. Ih and G. Tonuc, “Boundary element analysis using artificial boundary node approach,” Communications in Numerical Methods in Engineering, vol. 16, pp. 769–776, 2000.
View at: Publisher Site | Google Scholar
G. E. Blandford, A. R. Ingraffea, and J. A. Liggett, “Two-dimensional stress intensity factor computations using the boundary element method,” International Journal for Numerical Methods in Engineering, vol. 17, no. 3, pp. 387–404, 1981.
View at: Publisher Site | Google Scholar
D. E. Beskos, “Boundary element methods in dynamic analysis: Part II (1986-1996),” Applied Mechanics Reviews, vol. 50, no. 3, pp. 149–197, 1997.
View at: Publisher Site | Google Scholar
J. Lysmer, “Lumped mass method for Rayleigh waves,” Bulletin of the Seismological Society of America, vol. 60, no. 1, pp. 89–104, 1970.
View at: Publisher Site | Google Scholar
J. Lysmer and G. Waas, “Shear waves in plane infinite structures,” Journal of the Engineering Mechanics Division, vol. 98, no. 1, pp. 85–105, 1972.
View at: Publisher Site | Google Scholar
E. Kausel, J. M. Roesset, and G. Waas, “Dynamic analysis of footings on Layered media,” Journal of Engineering Mechanics Division,ASCE, vol. 101, no. 5, pp. 679–693, 1975.
View at: Publisher Site | Google Scholar
D. Givoli, “Recent advances in the DtN FE method,” Archives of Computational Methods in Engineering, vol. 6, no. 2, pp. 71–116, 1999.
View at: Publisher Site | Google Scholar
L. Ting and M. J. Miksis, “Exact boundary conditions for scattering problems,” Journal of the Acoustical Society of America, vol. 80, no. 6, pp. 1825–1827, 1986.
View at: Publisher Site | Google Scholar
D. Givoli and D. Cohen, “Nonreflecting boundary conditions based on Kirchhoff-type Formulae,” Journal of Computational Physics, vol. 117, no. 1, pp. 102–113, 1995.
View at: Publisher Site | Google Scholar
L. C. Qiu and F. Jin, “Explicit time-domain radiation boundaries for analysis of wave propagation in unbounded media,” Journal of Tsinghua University, vol. 42, no. 11, pp. 153–1533, 2003, (in Chinese).
View at: Publisher Site | Google Scholar
R. Vichnevetsky and E. C. Pariser, “High order numerical sommerfeld boundary conditions Theory and experiments,” Computers & Mathematics with Applications, vol. 11, no. 1-3, pp. 67–78, 1985.
View at: Publisher Site | Google Scholar
C. G. Torre, “Uniqueness of solutions to the helically reduced wave equation with sommerfeld boundary conditions,” Journal of Mathematical Physics, vol. 47, no. 7, Article ID 073501, 2006.
View at: Publisher Site | Google Scholar
R. W and B. Clayton, “Engquist. Absorbing boundary conditions for wave-equation migration,” Geophysics, vol. 45, no. 5, pp. 895–904, 1980.
View at: Publisher Site | Google Scholar
R. W and B. Clayton, “Engquist.Absorbing boundary conditions for acoustic and elastic wave equations,” Bulletin of the Seismological Society of America, vol. 25, pp. 383–397, 1977.
View at: Publisher Site | Google Scholar
R. L. Higdon, “Absorbing boundary conditions for difference approximations to the multi-dimensional wave equation,” Mathematics of Computation, vol. 47, no. 176, pp. 437–459, 1986.
View at: Publisher Site | Google Scholar
R. L. Higdon, “Numerical absorbing boundary conditions for the wave equation,” Mathematics of Computation, vol. 49, no. 179, pp. 65–90, 1987.
View at: Publisher Site | Google Scholar
R. L. Higdon, “Absorbing boundary conditions for acoustic and elastic waves in stratified media,” Journal of Computational Physics, vol. 101, no. 2, pp. 386–418, 1992.
View at: Publisher Site | Google Scholar
J. P. Wolf and C. Song, “Doubly asymptotic multi-directional transmitting boundary for dynamic unbounded medium-structure-interaction analysis,” Earthquake Engineering & Structural Dynamics, vol. 24, no. 2, pp. 175–188, 1995.
View at: Publisher Site | Google Scholar
L. P. Jing, Z. P. Liao, and J. X. Zou, “A high-frequency instability mechanism in numerical realization of multi-transmitting formula,” Earthquake Engineering and Engineering Vibration, vol. 22, no. 1, pp. 7–13, 2002, (in Chinese).
View at: Publisher Site | Google Scholar
L. P. Jing, Z. Y. Wu, and J. X. Zou, “Stability analysis for numerical simulation of near-field wave motion,” Earthquake Engineering and Engineering Vibration, vol. 22, no. 2, pp. 17–21, 2002, (in Chinese).
View at: Publisher Site | Google Scholar
J. Lp, “Stability analysis of practical formula for multi-transmitting boundary,” Earthquake Engineering and Engineering Vibration, vol. 24, no. 4, pp. 20–24, 2004, (in Chinese).
View at: Publisher Site | Google Scholar
J. Lysmer and R. L. Kuhlemeyer, “Finite dynamic model for infinite media,” Journal of the EngineeringMechanics Division ASCE, vol. 95, no. EM4, pp. 869–877, 1969.
View at: Publisher Site | Google Scholar
W. White, S. Valliappan, and I. K. Lee, “Unified boundary for finite dynamic models,” Journal of the Engineering Mechanics Division on ASCE, vol. 103, no. EMS, pp. 949–964, 1977.
View at: Publisher Site | Google Scholar
A. J. Deeks and M. F. Randolph, “Axisymmetric time domain transmitting boundaries,” Journal of Enginecring Mechanics, vol. 120, no. 1, pp. 25–42, 1994.
View at: Publisher Site | Google Scholar
H. F. Sun, L. P. Jing, X. C. Meng, and N. W. Wang, “The selection of boundary conditions in dynamic problems by using ABAQUS,” Earthquake Engineering and Engineering Vibration, vol. 11, no. 3, pp. 71–76, 2011.
View at: Google Scholar
Z. S. Aterman and F. C. Karal, “Propagation of elastic waves in layered media by finite difference methods,” Bulletin of the Seismological Society of America, vol. 58, pp. 367–398, 1968.
View at: Publisher Site | Google Scholar
J. B. Liu, Y. Gu, and Y. X. Du, “Consistent viscous-spring artificial boundaries and viscous-spring boundary elements,” Chinese Journal of Geotechnical Engineering, vol. 129, pp. 1070–1075, 2006.
View at: Google Scholar
G. J. Wang and C. G. Zhao, “Viscus-Spring dynamical artificial boundary for saturated porous media,” World Earthquake Engineering, vol. 38, no. 5, pp. 607–611, 2006, (in Chinese).
View at: Publisher Site | Google Scholar
Y. Gu and J. B. Liu, “Du YX. 3D consistent viscous-spring artificial boundary and viscous- spring boundary element,” Engineering Mechanics, vol. 24, no. 12, pp. 31–37, 2007, (in Chinese).
View at: Publisher Site | Google Scholar
J. B. Liu and Y. D. Lv, “A direct method for analysis of dynamic soil-structure interaction,” Journal of Civil Engineering, vol. 31, no. 3, pp. 55–64, 1998, (in Chinese).
View at: Google Scholar
J. B. Liu, Y. X. Du, and Xl Du, “3D viscous-spring artificial boundary in time domain,” Earthquake engineering and vibration, vol. 5, no. 1, pp. 93–102, 2006, (in Chinese).
View at: Publisher Site | Google Scholar
D. Z. Quan, S. B. Chai, and S. Chen, “Experimental study on seismic response of subway station built in loess,” Latin American Journal of Solids and Structures, vol. 17, no. 5, pp. 1–17, 2020.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2022 Dan Ye et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

239

Downloads

392

Citations