- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Journal of Applied Mathematics
Volume 2013 (2013), Article ID 241482, 12 pages
Simplified Boundary Element Method for Kinematic Response of Single Piles in Two-Layer Soil
1Key Laboratory of Geotechnical and Underground Engineering of Ministry of Education, Tongji University, Shanghai 200092, China
2Department of Geotechnical Engineering, Tongji University, Shanghai 200092, China
3School of Civil Mining and Environmental Engineering, University of Wollongong, Wollongong, NSW 2522, Australia
Received 15 May 2013; Accepted 23 July 2013
Academic Editor: Ga Zhang
Copyright © 2013 Fayun Liang 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.
A simple approach is formulated to predict the elastic, kinematic pile bending during harmonic or transient excitation for a circular pile (rather than a simplified thin strip). The kinematic response of a pile embedded in two-layer soil is resolved in the frequency domain caused by the upward propagation of shear waves from the underlying bedrock. The simplified approach is generally valid to nonhomogeneous soil profiles, in light of the good comparison with the dynamic FE method and BDWF solution. It employs the soil-displacement-influence coefficients to consider the pile-soil interaction (resembling the spring constant in the BDWF) and provides conservative estimations of maximum kinematic bending moments at the soil-layer interface (with a sharper stiffness contrast). The accuracy of the approach may be improved by incorporating the interaction of soil into the soil-displacement-influence coefficients for such cases with .
Kinematic response is one of the key issues in seismic design of pile foundations [1–4], as the dynamic response of the pile structure largely differs from the response at the free field soil caused by seismic wave. In practice, the influence of inertial loading at the pile-head level has conventionally received sufficient attention by design engineers, but not the kinematic seismic response . The existing earthquake investigations [6, 7] and experimental studies  demonstrate that at the interface of two-layer soils with a sharp stiffness contrast [2, 3, 9, 10], or at a pile-head with fixed constraints [11, 12], large kinematic bending moments may be induced to inflict damage to the pile. This problem has attracted the attention researchers [13–15] and is highlighted in some advanced seismic codes .
The kinematic response in the pile-soil system has been analyzed by considering the effect of the passive pile using rigorous mechanical solutions [17, 18], numerical methods [19–25], and some simplified models [2, 3, 26–30]. Rigorous solutions and numerical methods for the kinematic pile bending are, however, not convenient to design purposes. Empirical formulas are developed for evaluating the bending moment at the pile-head or at the interface of two-layer soils [2, 3, 11]. Among them, a Beam-on-Dynamic-Winkler-Foundation (BDWF) formulation was used successfully in practice, although it was confined to the harmonic excitation at the pile head. The BDWF (or the Winkler model) is underpinned by a frequency-dependent impedance , in which the continuously distributed springs is empirically related to a dimensionless coefficient (used for a given pile-soil system, regardless of layers). Accurate selection of the value is not critical to calculating the pile-head deflection but is important to predicting the kinematic bending moment and shear force [19, 31]. An optimized is thus required to obtain correct kinematic pile bending at the interface of soil layers. On the other hand, a simplified boundary element formulation proposed by Poulos and Davis  offers good estimates of bending moment and shear force for static loading. It would be good to see its accuracy in predicting the kinematic seismic response of pile, especially at the sharp stiffness contrasts between adjacent soil layers.
In this paper, the simplified approach  is employed to evaluate the kinematic pile bending during harmonic or transient excitation, concerning piles in two-layer soil. The pile is modeled as a circular shape rather than a thin strip adopted previously , and the soil displacement is given by the Mindlin equation with corresponding elastic modulus. A nodal relative displacement is obtained by one-dimensional site response and is then imposed on the pile. The solution was compiled into a program operating in Matlab platform. For some typical cases, a comprehensive study on the kinematic seismic response during harmonic or transient excitation has been carried out, and the results are compared with available dynamic FE method and the BDWF solutions. The study sheds new light on the kinematic bending moment at the interface of two-layer soil and at the pile head and may facilitate the use of the simplified boundary element method to predict the kinematic seismic response of a single pile.
2. Simplified Analysis Procedure
2.1. Basic Assumptions for Pile and Soil Model
The one-dimensional model for a floating or end-bearing single pile embedded in a two-layer soil is shown in Figure 1. The circle pile is assumed as linearly hysteretic beam having a length , a diameter , a mass density , and a bending stiffness . The pile is discretized into segments of equal lengths , but for a length of for the top and the tip segments, respectively. Each segment is subjected to a uniformly distributed load over the semicircular area. The pile is head restrained (fixed head) or free to rotate (free head), and sits above a bedrock. The linearly hysteretic soil profile is characterized by an upper-layer of thickness and a shear wave velocity , which is underlain by a lower layer of thickness and shear wave velocity . The two layers have damping ratios and , mass densities and and Poisson’s ratios and . A shear wave propagates vertically through the free field soil, which induces the horizontal harmonic motion and horizontal displacements. The motion at the bedrock surface is expressed by the amplitude of either bedrock displacement or the bedrock acceleration .
2.2. Calculation of the Horizontal Displacement of the Pile
Kinematic response of a single pile is induced by the free-field soil displacement shown in Figure 1. The Mindlin hypothesis does not meet the needs of dynamic analysis. However, the Mindlin equation is still valid for calculating elastic displacement and stress fields caused by a dynamic loading, provided that the characteristic wavelength in the soil is sufficiently long in comparison with the horizontal distance across the zone of major influence [28, 33], as is noted for nonhomogeneous soil by Poulos and Davis . In the current, simplified BEM formulation, the soil displacement (due to the pile-soil interface pressure) is gained using the Mindlin solution , which is then added together with free-field soil displacement where is the soil-displacement-influence coefficient; is the vector of soil-pile interface pressure over the semi-circular area; and is the free-field soil displacement estimated using one-dimensional site response for vertically propagating shear waves through an unbounded medium .
The dynamic equilibrium under steady-state conditions for the pile may be written in the following form using the finite-difference method : where is the horizontal displacement of the pile, with the cap “·” indicating differentiation with time; is the vector of soil-pile interface pressure; is the matrix of finite difference coefficients; is the pile mass; and is the soil radiation damping. Here, the soil damping is the same as that of the simplified boundary element model , and where is the shear wave velocity of soil and is the density of soil.
The displacement compatibility between the pile and the adjacent soil offers . Taking the displacement as and substituting (1) into (2) result in the following: Equation (4), together with the pile-top and -bottom boundary conditions, leads to unknown displacements, which involve the pile nodes from 1 to (see Figure 1) and 4 additional nodes at the pile top and tip. Equation (4) is resolved in either frequency or time domain and may attain the required accuracy using 21 segments of the pile .
2.3. Soil-Displacement-Influence Coefficients
Poulos and Davis  obtained soil-displacement-influence coefficients by integrating the Mindlin equation over a rectangular plane and taking the pile as a thin rectangular vertical strip. Ideally, a circular pile (rather than a thin strip) should be used. In elastic, semi-infinite space, the force in the horizontal direction at a depth induces a displacement component , which at any other point is given by where is shear modulus of soil and is the Poisson ratio.
Secondly, (6) is reexpressed in a cylindrical coordinate, which is then integrated over the semi-circular pile surface with respect to the angle (see Figure 2) to gain the horizontal component of the displacement. The soil-displacement-influence coefficients, taking as weighted average of the integrated horizontal displacements, are deduced asEquation (7) is solved by numerical integration using an adaptive Lobatto rule .
2.4. Harmonic Response in the Frequency Domain
Wave produces horizontal harmonic motion in the free field, as shown in Figure 1. This is described by , and the associated free-field horizontal displacement is given by . One-dimensional site response analysis can be formulated as where is the mass density, is viscosity, and is displacement. Equation (8) is resolved in frequency domain analysis [37, 38], allowing the nodal relative displacement to be obtained.
The current dynamic analysis employs time domain analysis and frequency domain analysis [21, 35, 39]. Frequency domain methods are widely used to estimate the dynamic impendences of the pile head. In the strong seismic motions, time domain method (involving the Newton-Raphson iteration and the Newmark method ) is used to obtain the nonlinear results. Equation (4) was resolved in time domain . In contrast, to facilitate comparison with the rigorous FE method and BDWF model, (4) is resolved herein in the frequency domain by the following form: where ; is the excitation frequency; is the amplitude of pile displacement; is the amplitude of addition displacement in the free-field soil; and is the identity matrix.
A cut-off method [28, 33] is generally used to accommodate soil yield around the pile. If the pressure at the pile-soil interface exceeds the ultimate lateral pressure of soil, the excess pressure is redistributed to other segments through iteration until pressure at all pile nodes within the ultimate values. This study, however, will not consider this yield and will be confined to elastic analysis using the simplified approach for piles in two-layer soil.
3. Validation of Simplified Method
3.1. Comparison with Dynamic Finite-Element Solution
The proposed simple approach is compared with dynamic FE results concerning a free head pile embedded in a two-layer soil deposit . The pile-soil model is the same as that shown in Figure 1, except that the pile tip is extended into the underlying bedrock. The pile-soil system is featured by a ratio of soil layer thickness of 1, a soil density , a soil Poisson’s ratio , and a soil damping coefficient . The pile has a slenderness ratio of 20, a pile-to-soil stiffness ratio of 5000, and a pile density of 1.60 ( kg/m3).
A comparison between the simplified approach and FE solution  is presented in Figures 3(a) and 3(b), respectively, for the profiles of the pile deflection and bending moment amplitude at the natural frequency of soil deposit . A good agreement is evident. The current, simplified approach can reveal the kinematic bending moments at the interface of the two layers, despite the ~20% overestimation of the maximum bending moment (against FE result) of the pile in Case D.
Figure 4 shows the amplitude spectrum of maximum kinematic bending moment as a function of the frequency ratio . A good agreement is again observed between the simplified approach and the dynamic FE solution . Both indicate that the peak kinematic bending moment occurs at the inherent frequency of the soil.
3.2. Comparison with BDWF Formulation
The proposed approach for kinematic loading along the pile depth during the lateral ground movements is compared with the BDWF solution . The pile-soil system is the same as the case just discussed in the last section. To examine the sensitivity of the parameters, four groups of 12 cases (Table 1) were studied, by maintaining soil density , soil Poisson’s ratio , soil damping coefficient , and pile density .
As shown in Table 1, the BDWF method adopts an optimized to obtain kinematic pile bending at the interface of two-layer soil. In contrast, the current method uses the displacement-influence coefficients to consider the pile-soil interaction (resembling the spring constant in the Winkler model) and may incorporate the interaction of soil along the pile to improve the accuracy. Nevertheless, when the ratio of the shear wave velocities of two soil layers exceeds 3, a larger than 15% error (compared with the FE method) in maximum kinematic bending moment may be seen using the simplified method. This is discussed next concerning Case 12 for kinematic bending at the two-layer interface and at the pile head (at the natural frequency of soil deposit).
3.2.1. Kinematic Pile Bending
Figure 5 shows the amplitude distributions of kinematic pile bending and shear force from the simplified method and the BDWF solutions. The kinematic pile bending and shear force profiles are slightly sensitive to the value . At , the distribution profiles of pile moment and shear force agree with each other between the simplified results and the BDWF solution, although the maximum moment of the simplified method is 17.2% larger than the BDWF solution.
Figure 6 provides the corresponding amplitude spectrum of maximum kinematic pile bending moment in the two-layer soil, owing to variation in the frequency ratio , among the FE, the BDWF, and the simplified approaches. It shows a consistent trend of variation in bending moment among various approaches.
3.2.2. Kinematic Response of Pile Head
Kinematic responses are obtained in form of the ratio of the amplitudes of pile-head displacement over the excitation motion or the ratio of the head displacement over the free field surface displacement [17, 26, 39]. The responses for the free-head pile are plotted in Figure 7 for a spectrum of the frequency ratio . The good agreement of the factors and among the simplified approach, the FE method, and the BDWF solution has been attributed to the predominant effect of the free field soil displacement.
As for fixed-head piles, equally successful prediction is seen in Figure 8, concerning the kinematic response of pile to soil displacement among the simplified method, the rigorous boundary integral method (Fan et al. ), and the BDWF (Makris and Gazetas ). Nevertheless, the current approach offers slightly larger ratio for than the other two solutions at high frequencies.
4. Application under Seismic Excitation
4.1. Seismic Motion and Case Model
The harmonic steady state is rarely seen in a practical engineering design. Kinematic seismic response of a pile should be tailored to cater for the transient excitation  owing to earthquake shaking, as is noted in dynamic analysis involving nonlinear pile-soil interaction. A large soil resistance in certain depth may render the nonlinear kinematic pile bending (e.g., caused by SH wave) insignificant compared to pile-head inertial excitation. In addition, previous study  does not allow either a detrimental or a beneficial effect on kinematic pile bending to be concluded due to non-linear site response. Consequently, the validity of the simplified method is examined herein for elastic pile and soil.
The performance of the simplified approach for the transient response is examined for six typical seismic accelerograms (see Table 2), which include 4 actual records selected from ground motion database of the Pacific Earthquake Engineering Research Center (PEER)  and 2 artificial motions used in seismic design of a typical site in Shanghai (China). The acceleration time histories for seismic events are plotted in Figure 9 and the associated acceleration response spectra are provided in Figure 10.
4.2. Parametric Investigation
The parametric analysis is again conducted for the pile-soil system shown in Figure 1, with the following profile parameters: bedrock located at m, density of either soil layer = 1900 kg/m3, Poisson's ratio = 0.4, soil damping = 5%, pile density = 2500 kg/m3, and pile-to-soil stiffness ratio . The shear wave velocity is 100 m/s for the upper soil layer, 150 m/s for the lower layer, and 1000 m/s for the bedrock, respectively. The input signals by Sica et al.  were scaled in amplitude to a peak acceleration of 0.35 g . Figure 11 shows the comparison of kinematic pile bending between the simplified proposed approach and the BDWF formulation . A good agreement in the predicted bending moment is evident between the simplified method and the BDWF solution, but for the large difference in the peak bending moment around the layer interface and at a shear wave velocity ratio of 1/3.
Figure 12 shows the pile-diameter contrast on kinematic bending moment under the six input motions. The ratio of soil layer thickness is 1, the pile is 20 m in length and 0.6, 0.9, or 1.2 m in diameter, and the head is constrained against rotation (fixed head). The figure indicates that small-diameter piles accommodate more easily to seismically induced soil deformations than lager-diameter piles. Kinematic bending moments at the pile head and the interface of soil layers are nearly proportional to the diameters. This will reduce the safety of pile head, although it does not necessarily increase or reduce the seismic safety of the pile body (depending on the circumstances), as is discussed previously [10, 14].
Figure 13 presents the bending moment profiles owing to variation in the depths of the interface of the two layers under the six input motions. They were obtained for a fixed head pile (20 m in length, 0.6 m in diameter) embedded in a soil with a layer interface located at a depth of 5, 10, or 15 m, respectively. An increase in depth of the interface renders increase in the peak value of kinematic pile bending but has a negligible impact on the kinematic bending moment at the pile head. This preliminary analysis suggests the simplified approach has the potential in modeling kinematic seismic response of piles during the earthquake.
5. Conclusion and Discussions
A simple approach is formulated to predict the elastic, kinematic pile bending during harmonic or transient excitation. The approach employs a circular pile (rather than a simplified thin strip). The kinematic response of a pile embedded in two-layer soil is resolved in the frequency domain owing to specified soil displacement field. The simplified approach is generally valid to nonhomogeneous soil problems, in light of the good comparison with the dynamic FE method and BDWF solution. The main conclusions from the study are as follows.
The simplified method employs the soil-displacement-influence coefficients to consider the pile-soil interaction (resembling the spring constant in the BDWF). It provides conservative estimation of maximum kinematic bending moments at soil-layer interface (with a sharper stiffness contrast) despite an adequate accuracy in general. The accuracy may be improved by incorporating the interaction of soil into the soil-displacement-influence coefficients for such cases with .
The parametric studies during the seismic excitation show the impact of pile diameters and depths of the interface of two layers and demonstrate the simplified approach comparing well with published results in gaining kinematic pile bending during the earthquake.
The formulated simple method is intended for elastic soil, elastic pressure on pile-soil interface, and elastic seismic response of the free field soil, for which a judicious choice of elastic modulus of the soil is required. The method may be extended to elastoplastic case though introducing a similar technology as cut-off method [28, 33] in the current procedure.
The work reported herein was supported by the National Natural Science Foundation of China (Grants nos. 41172246, 90915011), the National Key Basic Research Program of China (Grant no. 2013CB036304), and the Fundamental Research Funds for the Central Universities, Tongji University. The above financial support is gratefully acknowledged.
- W. D. L. Finn, “A study of piles during earthquakes: issues of design and analysis,” Bulletin of Earthquake Engineering, vol. 3, no. 2, pp. 141–234, 2005.
- S. Nikolaou, G. Mylonakis, G. Gazetas, and T. Tazoh, “Kinematic pile bending during earthquakes: analysis and field measurements,” Geotechnique, vol. 51, no. 5, pp. 425–440, 2001.
- G. Mylonakis, “Simplified model for seismic pile bending at soil layer interfaces,” Soils and Foundations, vol. 41, no. 4, pp. 47–58, 2001.
- B. A. Bradley, M. Cubrinovski, and J. J. M. Haskell, “Probabilistic pseudo-static analysis of pile foundations in liquefiable soils,” Soil Dynamics and Earthquake Engineering, vol. 31, no. 10, pp. 1414–1425, 2011.
- J. P. Stewart, R. B. Seed, and G. L. Fenves, “Seismic soil-structure interaction in buildings. II: empirical findings,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 125, no. 1, pp. 38–48, 1999.
- T. Tazoh, K. Shimizu, and T. Wakahara, “Seismic observations and analysis of grouped piles,” in Dynamic Response of Pile Foundations-Experiment, Analysis and Observation (GSP 11), pp. 1–20, 1987.
- K. Tokimatsu, “Behaviour and design of pile foundations subjected to earthquakes,” in Proceedings of the 12th Asian Regional Conference on Soil Mechanics and Geotechnical Engineering,, vol. 1 and 2, Singapore, 2003.
- K. Tokimatsu, H. Suzuki, and M. Sato, “Effects of inertial and kinematic interaction on seismic behavior of pile with embedded foundation,” Soil Dynamics and Earthquake Engineering, vol. 25, no. 7-10, pp. 753–762, 2005.
- R. M. S. Maiorano, L. De Sanctis, S. Aversa, and A. Mandolini, “Kinematic response analysis of piled foundations under seismic excitation,” Canadian Geotechnical Journal, vol. 46, no. 5, pp. 571–584, 2009.
- F. Dezi, S. Carbonari, and G. Leoni, “Kinematic bending moments in pile foundations,” Soil Dynamics and Earthquake Engineering, vol. 30, no. 3, pp. 119–132, 2010.
- L. de Sanctis, R. M. S. Maiorano, and S. Aversa, “A method for assessing kinematic bending moments at the pile head,” Earthquake Engineering & Structural Dynamics, vol. 39, no. 10, pp. 1133–1154, 2010.
- R. Di Laora, G. Mylonakis, and A. Mandolini, “Pile-head kinematic bending in layered soil,” Earthquake Engineering & Structural Dynamics, vol. 42, pp. 319–337, 2013.
- S. Sica, G. Mylonakis, and A. L. Simonelli, “Strain effects on kinematic pile bending in layered soil,” Soil Dynamics and Earthquake Engineering, vol. 49, pp. 231–242, 2013.
- R. Di Laora, A. Mandolini, and G. Mylonakis, “Insight on kinematic bending of flexible piles in layered soil,” Soil Dynamics and Earthquake Engineering, vol. 43, pp. 309–322, 2012.
- M. Saitoh, “Fixed-head pile bending by kinematic interaction and criteria for its minimization at optimal pile radius,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 131, no. 10, pp. 1243–1251, 2005.
- R. Di Laora and A. Mandolini, “Some remarks about Eurocode and Italian code about piled foundations in seismic area,” in Presented at the ERTC-12 Workshop on Evaluation of EC8, Athens, Greece, 2011.
- A. M. Kaynia and E. Kausel, “Dynamic stiffness and seismic response of pile groups,” Research Report R82-03, Massachusetts Institute of Technology, Cambridge, Mass, USA, 1982.
- K. Fan, G. Gazetas, A. Kaynia, E. Kausel, and S. Ahmad, “Kinematic seismic response of single piles and pile groups,” Journal of Geotechnical Engineering, vol. 117, no. 12, pp. 1860–1879, 1991.
- M. Kavvadas and G. Gazetas, “Kinematic seismic response and bending of free-head piles in layered soil,” Geotechnique, vol. 43, no. 2, pp. 207–222, 1993.
- M. H. El Naggar and M. Novak, “Nonlinear analysis for dynamic lateral pile response,” Soil Dynamics and Earthquake Engineering, vol. 15, no. 4, pp. 233–244, 1996.
- G. Wu and W. D. L. Finn, “Dynamic nonlinear analysis of pile foundations using finite element method in the time domain,” Canadian Geotechnical Journal, vol. 34, no. 1, pp. 44–52, 1997.
- B. K. Maheshwari, K. Z. Truman, M. H. El Naggar, and P. L. Gould, “Three-dimensional nonlinear analysis for seismic soil-pile-structure interaction,” Soil Dynamics and Earthquake Engineering, vol. 24, no. 4, pp. 343–356, 2004.
- F. Dezi, S. Carbonari, and G. Leoni, “A model for the 3D kinematic interaction analysis of pile groups in layered soils,” Earthquake Engineering & Structural Dynamics, vol. 38, no. 11, pp. 1281–1305, 2009.
- M. A. Millán and J. Domínguez, “Simplified BEM/FEM model for dynamic analysis of structures on piles and pile groups in viscoelastic and poroelastic soils,” Engineering Analysis with Boundary Elements, vol. 33, no. 1, pp. 25–34, 2009.
- F. Dezi, S. Carbonari, and G. Leoni, “Static equivalent method for the kinematic interaction analysis of single piles,” Soil Dynamics and Earthquake Engineering, vol. 30, no. 8, pp. 679–690, 2010.
- N. Makris and G. Gazetas, “Dynamic pile-soil-pile interaction. Part II: lateral and seismic response,” Earthquake Engineering & Structural Dynamics, vol. 21, no. 2, pp. 145–162, 1992.
- G. Mylonakis, A. Nikolaou, and G. Gazetas, “Soil-pile-bridge seismic interaction: kinematic and inertial effects. Part I: soft soil,” Earthquake Engineering & Structural Dynamics, vol. 26, no. 3, pp. 337–359, 1997.
- A. Tabesh and H. G. Poulos, “Pseudostatic approach for seismic analysis of single piles,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 127, no. 9, pp. 757–765, 2001.
- T. Juirnarongrit and S. A. Ashford, “Soil-pile response to blast-induced lateral spreading. II: analysis and assessment of the p-y method,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 132, no. 2, pp. 163–172, 2006.
- F. Castelli and M. Maugeri, “Simplified approach for the seismic response of a pile foundation,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 135, no. 10, pp. 1440–1451, 2009.
- S. Sica, G. Mylonakis, and A. L. Simonelli, “Transient kinematic pile bending in two-layer soil,” Soil Dynamics and Earthquake Engineering, vol. 31, no. 7, pp. 891–905, 2011.
- H. G. Poulos and E. H. Davis, “Pile foundation analysis and design,” 1980.
- H. Elahi, M. Moradi, H. G. Poulos, and A. Ghalandarzadeh, “Pseudostatic approach for seismic analysis of pile group,” Computers and Geotechnics, vol. 37, no. 1-2, pp. 25–39, 2010.
- R. D. Mindlin, “Force at a point in the interior of a semi-infinite solid,” Journal of Applied Physics, vol. 7, no. 5, pp. 195–202, 1936.
- Y. M. A. Hashash, C. Phillips, and D. R. Groholski, “Recent advances in non-linear site response analysis,” in Proceedings of the 5th International Conference in Recent Advances in Geotechnical Eartqhuake Engineering and Soil Dynamics, 2010.
- L. F. Shampine, “Vectorized adaptive quadrature in MATLAB,” Journal of Computational and Applied Mathematics, vol. 211, no. 2, pp. 131–140, 2008.
- D. Park and Y. M. A. Hashash, “Soil damping formulation in nonlinear time domain site response analysis,” Journal of Earthquake Engineering, vol. 8, no. 2, pp. 249–274, 2004.
- G. A. Ordonez, “SHAKE2000: a computer program for the 1-D analysis of geotechnical earthquake engineering problems,” 2004.
- G. Wu and W. D. L. Finn, “Dynamic elastic analysis of pile foundations using finite element method in the frequency domain,” Canadian Geotechnical Journal, vol. 34, no. 1, pp. 34–43, 1997.
- R. W. Boulanger, C. J. Curras, B. L. Kutter, D. W. Wilson, and A. Abghari, “Seismic soil-pile-structure interaction experiments and analyses,” Journal of Geotechnical and Geoenvironmental Engineering, vol. 125, no. 9, pp. 750–759, 1999.
- P. G. Motion, “Pacific earthquake engineering research center,” Tech. Rep., University of California, Berkeley, Berkeley, Calif, USA, 2008.