Preprint
Article

This version is not peer-reviewed.

Dynamic and Seismic Response of Anchored Pile Walls Under Earthquake Loading

Submitted:

18 September 2026

Posted:

01 October 2026

You are already at the latest version

Abstract
The purpose of this study is to present a numerical analysis of the behavior of an anchored concrete pile wall subjected to combined hydrostatic, dynamic, and seismic loads. It is based on a specific construction site at Gazino Station on the Ulus-Keçiören route. There are multiple layers of soil, a four-story building adjacent, and a variable groundwater table in the geotechnical context. PLAXIS 2D Version 8 was used to develop a 2D finite element model. The pile wall was modeled as a beam element with a soil-structure interface, while the ground anchors were modeled as a simplified combination of node-to-node anchors (for the tendon) and geotextile elements (for the grout body). A three-phase analysis methodology is used: (1) a hydrostatic analysis involving a three-stage rise in groundwater level; (2) a dynamic analysis with a concentrated load applied at the wall crest; and (3) a seismic analysis using a basement-recorded earthquake accelerogram from 1989. Increasing lateral pressures resulted in a reduction in wall displacement, which was attributed to rising groundwater levels. In contrast, both dynamic load and seismic excitation caused significant displacements, with pronounced effects around anchor heads and between the existing building and the wall. There was a significant disparity between the displacements at the top and bottom of the wall during seismic loading at approximately 4.08 seconds. As a result of the study, hydraulic loading conditions may not be as critical as dynamic and seismic loads. Material properties as well as structural details should be optimized for pile wall stability and adjacent structure integrity.
Keywords: 
;  ;  ;  ;  ;  ;  

1. Introduction

The problem was modeled using a geometry of 35 m width and 20 m depth. A ground anchor was modeled by a combination of a node-to-node anchor and a geotextile element. The geotextile simulates the grout body, while the node-to-node anchor simulates the anchor rod. In reality, a complex three-dimensional state of stress exists around the grout body. Although the precise stress state and soil interaction cannot be fully captured with this 2D model, this approach allows for the estimation of stress distribution, deformations, and global structural stability, assuming the grout body does not slip relative to the soil. With this model, it is not possible to evaluate the pull-out force of the ground anchor.
The pile wall is modeled as a beam. Interfaces around the beam are used to model soil-structure interaction effects. These interfaces are extended 1.0 m under the wall. Interfaces should not be used around the geotextiles representing the grout body.
The excavation is constructed in several stages. The separation between stages is modeled by geometry lines.
The water table is initially below the excavation level. A stepwise rise in the water level was simulated, and its pressure effects on the pile wall and existing building were analyzed. During the excavation stages, water is drained, and no groundwater is present. A schematic view of the analyzed model for the studied section of Gazino Station is shown in Figure 1.
Figure 1. Modeled section.
Figure 1. Modeled section.
Preprints 234050 g001

2. Material and Methods

2.1. Soil and Interface Properties

The properties of the three soil layers and interface materials were adopted from the PLAXIS Tutorial Manual (Version 8), as shown in Table 1.

2.2. Properties of the Plates

The properties of the structural plate elements—namely the pile wall, pile cap, and building beams—were also adopted from the PLAXIS Tutorial Manual (Version 8), as detailed in Table 2, Table 3 and Table 4.

2.3. Properties of the Anchors

The properties for the anchor rods and grout bodies were assumed from the same source, as shown in Table 5 and Table 6.

3. Steps and Procedures

Plaxis version 8 has been used for this study, all above said properties of materials and plates are filled, the study including three scopes as following:

3.1. Study of Pile Wall Behavior Under Water Pressure

In three stages, the groundwater level was raised starting from the bottom of the excavation and ending with the submersion of the entire pile wall.

3.2. Study of Pile Wall Behavior Under Dynamic Load

In this case, the top of the pile wall will be subjected to dynamic load.

3.3. Study of Pile Wall Behavior Under Earthquake

As part of this, an earthquake recorded by the USGS in 1989 was applied to the basement.

4. Results and Discussion

Using the PLAXIS 8 software, three distinct loading conditions (water pressure, dynamic load, and seismic load) were simulated on the numerical model of the anchored pile wall. The following presents a detailed account of the obtained results along with a comprehensive discussion of their engineering implications.

4.1. Pile Wall Behavior Under Water Pressure

Results
The simulation results revealed a decrease in wall displacement as the groundwater level was progressively raised across three stages. In the initial phase (dry wall), the highest displacement value was recorded. As the lower part of the wall (4 m) became submerged in Phase 4, a reduction in displacement was observed. This decrease became more pronounced in Phase 6 when the wall was fully submerged, as illustrated in Figure 2, which compares the output phases.
Discussion
This initially counter-intuitive behavior is attributed to the increase in Total Lateral Pressure on the wall with the rising water table. The pore water pressure contributes to an increase in the total pressure on the back face of the wall, generating a counteracting reaction and enhancing soil resistance. Additionally, changes in the mechanical properties of saturated soil (such as Young’s modulus and shear strength) may play a role in improving system stability. This finding suggests that, under these modeling assumptions, hydrostatic loading is not a critical factor and may even have a stabilizing effect. However, it is important to note that this conclusion is based on the assumption of no seepage or water flow. A study incorporating steady-state seepage conditions might yield different results, potentially highlighting destabilizing effects from water flow pressures, as it clearly in Figure 1.
Figure 1. Comparison of Output phases.
Figure 1. Comparison of Output phases.
Preprints 234050 g002

4.2. Pile Wall Behavior Under Dynamic Load

Results:
The application of a concentrated dynamic load at the top of the wall caused a continuous increase in horizontal displacements over the simulation time, with a marked acceleration noted in Phase 3 immediately after load application. The distribution of Plastic Points—areas within the finite element mesh where stresses have reached yield conditions—showed a clear concentration in two key zones:
  • The area surrounding the anchor heads.
  • There is a zone between the existing building and the pile wall.
In Figure 2, we can see how these plastic points are distributed.
Discussion:
The concentration of plastic points around the anchor heads indicates the development of high-stress concentrations in these regions under dynamic loading. This could be a precursor to potential anchor pull-out failure or exceedance of the bearing capacity of the soil surrounding the grout body. This aligns with the study’s stated limitation that the simplified 2D model cannot accurately evaluate anchor pull-out forces. The significant effect in the zone between the building and the wall reflects the complex dynamic interaction between the existing structure and the retaining wall under external loading, which can lead to response amplification in this area. These results underscore the critical importance of detailed design for anchor connection points and the necessity of analyzing structure-soil-structure interaction when dynamic loads are anticipated.
Figure 2. shows the plastic points.
Figure 2. shows the plastic points.
Preprints 234050 g003

4.3. Pile Wall Behavior Under Earthquake

Results:
Applying the 1989 USGS-recorded earthquake accelerogram to the model base induced both vertical and horizontal displacements. The maximum vertical displacement occurred at the basement level. The results also showed that displacement values at the top and bottom of the wall were similar for most of the simulation period. However, at a dynamic time of 4.08 seconds, a significant and sudden disparity emerged, with the displacement at the top of the wall becoming substantially greater than at the bottom, as shown in the time-displacement curve in Figure 3.
Discussion:
The pronounced disparity in displacement between the top and bottom at 4.08 seconds suggests the occurrence of a resonant response or a critical moment where the inertial effect of the upper portion of the wall is amplified. This specific time likely coincides with a natural frequency or a significant harmonic of the soil-wall-anchor system, leading to amplified motion at the wall crest. This behavior reveals a potential lack of system ductility or insufficient flexural resistance when subjected to seismic motion with a specific frequency content. This finding highlights the importance of conducting response spectrum analyses and designing retaining walls to accommodate significant seismic deformations without global stability failure, taking into account the unique dynamic characteristics of each site, including soil amplification effects, which is clearly visible in Figure 3.

5. Conclusion

According to the numerical analysis, the anchored pile wall responds critically to dynamic and seismic loads and shows stability under rising groundwater levels. A resonant response was induced by seismic excitation at 4.08 seconds, causing significant displacement at the top of the wall. As a result of dynamic loading, anchors and the existing building were subjected to high stress concentrations. Therefore, seismic and dynamic forces govern the design, not hydrostatic pressure. Optimizing wall stiffness, anchor detailing, and performing dynamic response analyses are recommended. Moreover, controlled deformations are necessary to ensure the safety of adjacent structures. In urban excavations, anchored earth-retaining systems should be designed taking dynamic considerations into account

References

  1. Plaxis, B.V. *PLAXIS 2D - Tutorial Manual Version 8*. Delft, Netherlands. (Essential: This is your primary source for material parameters and modeling procedures. You must cite it.). 2008.
  2. Plaxis, B.V. *PLAXIS 2D - Scientific Manual (Material Models)*. Delft, Netherlands. (Essential: Provides the theoretical foundation for the Mohr-Coulomb (MC) model and dynamic calculations used in your analysis.). 2008.
  3. Ulrich, N.; Sørensen, S. I. Seismic Analysis of a Deep Excavation Supported by Anchored Pile Walls using PLAXIS 2D Dynamic Module. In Proceedings of the 15th World Conference on Earthquake Engineering, Lisbon. (Directly relevant: Covers the exact topic—PLAXIS 2D dynamic analysis of an anchored pile wall under seismic load.); 2012. [Google Scholar]
  4. Bilgin, Ö.; Erken, A. Numerical Modeling of Groundwater Effect on the Behavior of a Braced Excavation. Geotech. Geol. Eng. (Highly relevant: A study specifically on how water table changes affect excavation support systems, supporting your water pressure results.). 2015, 33(4), 999–1013. [Google Scholar]
  5. Yoo, C.; Lee, D. Deep excavation-induced ground surface movement characteristics – A numerical investigation. Comput. Geotech. 2008, 35(2), 231–252. [Google Scholar] [CrossRef]
  6. Obrzud, R. F.; Truty, A. The Hardening Soil Model – A Practical Guidebook. Zace Services Ltd., Lausanne. (Advanced context: While you used the MC model, citing this guide on more advanced soil models strengthens your discussion on material behavior limitations.). Lausanne, 2010.
  7. Pelekis, P.; Xenaki, V.; Athanasopoulos, G. Effect of near-fault ground motions on the behavior of anchored retaining walls. Soil Dyn. Earthq. Eng. (Critical for discussion: Examines how specific seismic records affect anchored walls, providing a basis to discuss your use of the 1989 USGS record.). 2013, 45, 41–53. [Google Scholar]
  8. Long, M. Database for Retaining Wall and Ground Movements due to Deep Excavations. J. Geotech. Geoenvironmental Eng. ASCE Provides validation context: A key database of case histories against which your numerical results can be contextually discussed.. 2001, 127(3), 203–224. [Google Scholar] [CrossRef]
  9. Schanz, T.; Vermeer, P. A.; Bonnier, P. G. The Hardening Soil Model: Formulation and Verification. Beyond 2000 in Computational Geotechnics (pp. 281-296). Balkema, Rotterdam. (Foundational: The original paper for a key PLAXIS soil model. Citing it shows depth in understanding the software’s capabilities beyond the tutorial.). 1999; pp. 281–296.
  10. İyisan, R.; Yılmaz, H. A case study on the behavior of a deep excavation in clayey soils of Istanbul. Bull. Earthq. Eng. 2012, 10(5), 1621–1635. [Google Scholar]
Figure 3. The pile wall time Vs Displacement.
Figure 3. The pile wall time Vs Displacement.
Preprints 234050 g004
Table 1. Soil and interface properties.
Table 1. Soil and interface properties.
Parameter Name Fill Sandy Clay Andesite Building Fill Unit
Material model Model MC MC MC MC
Type of material behavior Type drained drained drained drained -
Dry soil weight dry 17,5 20 24 17 kN/m3
Wet soil weight wet 19 23 27 20 kN/m3
Horizontal permeability Kx 0 0 0 5 m/day
Vertical permeability Ky 0 0 0 5 m/day
Young’s modulus Etef 4000 20000 200000 30000 kN/m2
Poisson’s ratio 0,353 0,346 0,263 0.2 -
Cohesion Cref 0,2 10 100 1.0 kN/m2
Friction angle 27 28 40 35 0
Dilatancy angle 0 0 10 4 o
Interface reduction factor Rinter Rigid Rigid Rigid Rigid -
Interf Permeability parameter Perm. -
Table 2. Pile (beam) Properties.
Table 2. Pile (beam) Properties.
Parameter Name Value Unit
Type of behavior Material type Elastoplastic _
Normal stiffness EA 2,4-107 kN/m
Flexural rigidity El 1,28-106 kNm/m
Equivalent thickness D 0,80 m
Weight W 19,2 kN/m/m
Poisson’s ratio 0,15 -
Table 3. Pile cap (beam) Properties.
Table 3. Pile cap (beam) Properties.
Parameter Name Value Unit
Type of behavior Material type Elastoplastic _
Normal stiffness EA 1,5-107 kN/m
Flexural rigidity El 3,125-105 kNm/m
Equivalent thickness D 0.50 m
Weight W 12 kN/m/m
Poisson’s ratio 0,15 -
Table 4. Building members Properties.
Table 4. Building members Properties.
Parameter Name Value Unit
Type of behavior Material type Elastoplastic _
Normal stiffness EA 5-106 kN/m
Flexural rigidity El 9-106 kNm/m
Equivalent thickness D 0.147 m
Weight W 5 kN/m/m
Poisson’s ratio 0 -
Table 5. Anchor rod (node-to-node anchors) Properties.
Table 5. Anchor rod (node-to-node anchors) Properties.
Parameter Name Value Unit
Type of behavior Material type
Ls
Elastoplastic
,6
kN
m
Normal stiffness EA 128000
kN
Spacing out of plane Maximum Ls 1,6 kN
Table 6. Grout body (geotextile)Properties.
Table 6. Grout body (geotextile)Properties.
Parameter Name Value Unit
Normal stiffness EA 3,0-106 kN/m
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.