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:
PLAXIS
; pile wall
; soil-structure interaction
; water pressure
; dynamic load
; seismic analysis
; anchored retaining wall
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.

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
2.3. Properties of the Anchors
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.

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.

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
- 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.
- 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.
- 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]
- 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]
- Yoo, C.; Lee, D. Deep excavation-induced ground surface movement characteristics – A numerical investigation. Comput. Geotech. 2008, 35(2), 231–252. [Google Scholar] [CrossRef]
- 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.
- 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]
- 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]
- 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.
- İ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.

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.
| 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.
| 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.
| 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.
| 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.
| 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. |
© 2026 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
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.