# Soil-Structure Interaction of Flexible Temporary Trench Box: Parametric Studies Using 3D FE Modelling

Authors: Miah Alam; Omar Chaallal; Bertrand Galy
Journal: Modelling and Simulation in Engineering
 (2021)

---

## Abstract

This paper presents the results of two parametric finite-element studies that were carried out using the PLAXIS-3D finite element (FE) computer code. The following objectives and corresponding parameters were considered: (i) to evaluate the soil pressure on the steel trench box shield; the parameters studied were related to soil type and material, and the study considered till, dry sand, wet sand, and sensitive clay soil; (ii) to assess the effect of trench box material and geometry on earth pressure; the parameters studied were related to trench box material (steel versus aluminum) as well as geometry (plate thickness and strut diameter). These studies included simulation of two steel (or aluminum) trench box shields stacked upon each other to cover the total 6 m (20 ft) deep trench. A Mohr-Coulomb (MC) constitutive material model was chosen for FE analysis (FEA). The FEA results were compared to empirical apparent earth pressure diagrams for a sensitive clay. Comparisons showed that the parameters related to the soil and the trench box have a significant influence on earth pressures.


---

## Body

## 1. Introduction

Evaluation of horizontal earth pressure is of major importance when designing any soil retaining structure [1, 2]. In the last few years, numerous analytical studies have been carried out to better understand the behaviour of braced excavation protection systems [3]. The trench box shield is different from other shoring systems as it is intended not only to support the trench face but also to protect workers from cave-ins. Trench boxes can be assembled together to create a box system capable of shoring up trenches up to 11 m (35 ft) deep [4]. This system is relatively easy to transport, assemble, and install in narrow trenches compared to retaining structures built on site such as soldier pile walls. In order to keep the workers safe, shields should not be overloaded [5]. The earth pressure distribution on a shield depends on the shoring, the method of installation, and the type of soil. Many theoretical solutions to estimate earth pressure on a flexible temporary support have been published [6–8]. The Canadian Foundation Engineering Manual [9] recommends the use of pressure envelopes for excavation trench shields.Numerical FE and analytical methods helped understanding the behaviour of excavations and trenches with flexible retaining walls. Under a given set of soil conditions, the higher the wall rigidity, the greater was the apparent earth pressure [10]. Hashash and Whittle [11] studied shoring excavations in clay and showed that the earth pressure on a retaining structure (or wall) is a function of the flexibility of the retaining structure and the magnitude of shear deformations at depth. The arching effect in the ground in contact with the shoring also plays a significant role to estimate earth pressure on walls [3, 12–15]. In soft clays, earth pressure on supporting structures is a function that depends on the thickness of the clay layer in the trench bottom, the soil anisotropy, the flexibility of the support, and soil deformation [3].Horizontal and vertical strut spacing and wall stiffness also influence the earth pressure distribution on excavation support systems (ESS) and have been studied by Bryson and Zapata-Medina [16] for three types of clays. These researchers also defined a relative stiffness ratio that relates the stiffness of the soil to the stiffness of the ESS. This ratio depends on multiple factors, including Young’s modulus of the wall, Young’s modulus of the soil, moment of inertia per unit length of the wall, average vertical support spacing, average horizontal support spacing, excavation depth, total height of the wall, average unit weight of the soil, and undrained shear strength of the soil at the bottom of the excavation. Therefore, the interaction between soil mass and support structure is significant in the development of apparent earth pressure. A recently published paper by Alam et al. [17] showed that modern 3D FEM studies on clay trench SSI focused mainly on a depth greater than 10 m and on retaining wall or diaphragm wall. Not that many studies have been published on FEM of soil-structure interaction (SSI) with trench boxes in shallow depth trenches of clayey soil. Using 3D FEA and Abaqus software, LaBaw [18] evaluated earth pressure for a standard construction shoring (pneumatic shoring with steel panel) system, which provides an adequate and safe emergency rescue for first responders in medium clay of 2.4 m depth of trench excavation. Soil-structure interaction analysis is a powerful tool for analyzing, designing systems, and monitoring soil behaviour behind retaining wall under different types of loading in expansive, swelling, and clayey soil. Alam et al. [19] compared a full-scale experimental in situ test of 6 m depth with 3D FEM to investigate the earth pressure on a flexible temporary trench box shield in soft and sensitive clay soil using PLAXIS-3D software.This paper presents the influence of soil parameters and shield material type and geometry on the earth pressures on a flexible temporary trench box shield. An FE model using PLAXIS-3D software [20] shows a comprehensive picture of the influence of soil parameters on earth pressures on a flexible temporary trench box shield in till, dry sand, wet sand, and sensitive clay. These FEM analysis results are also compared to published empirical apparent earth pressure diagrams [7]. The influence of shield material type and geometry on the earth pressure in a trench excavation is also presented.


## 2. Finite-Element Analyses

Numerical analyses of ESS performance for recent deep excavation case histories have been done using three-dimensional FE analyses [21–23]. Ou et al. [24] showed that excavation-induced ground movements and soil-structure interactions of ESS are a three-dimensional (3D) phenomenon. FEA software may be a useful tool to obtain in-depth information on the effect of SSI on this type of trench box wall support [17]. To perform this 3D FEA and study the soil-structure interactions, the PLAXIS-3D software was used.


## 3. Parameters for FE

To simulate the trench and steel trench box shield, three sets of parameters were used in PLAXIS-3D: (i) soil model parameters; (ii) structural model parameters for the steel trench box, and (iii) interface parameters.


### 3.1. Soil Modelling

To simulate the trench, the soil model geometry ((7.5m25ft×12m40ft×15m50ft) was more than twice the size of the excavation contour (3m10ft×1.5m5ft×6m20ft) following recommendations by Hsieh and Ou [25] for shallow trench excavations (Figure 1(b)). The soil layers were modeled as follows: (i) fissured brown clay from 0 m to 0.6 m (2 ft) depth; (ii) plastic brown clay from 0.6 m (2 ft) to 2 m (6.56 ft) depth; and (iii) sensitive blue clay (Louiseville clay) from 2 m (6.56 ft) to 15 m (50 ft) depth, in accordance with the stratigraphic profile of the Louiseville soil [19]. Soil strength properties were selected carefully as they are a key factor in ESS performance [26]. The in situ value of undrained shear strength shown in Figure 1(a) was used in this FEA. On the other hand, till, dry sand, and wet sand layers were modelled using 0.0 to 15 m depths of the properties of each type of sand, as shown in Figure 1(c), with varying stiffness parameters along the depth.Figure 1
Geometric profile and excavation contour in PLAXIS-3D simulation of trench: (a) experimental values of undrained shear strength of sensitive blue clay (Louiseville clay, data from [33]); (b) PLAXIS-3D model profile for sensitive blue clay (number of soil elements: 23221; MC soil model); (c) PLAXIS-3D model profile for till and dry and wet sand (number of soil elements: 23221; MC soil model).
(a)(b)(c)


### 3.2. Structural Modelling

The stiffness of an ESS is a function of the flexural rigidity of the wall element, the structural stiffness of the support elements, the type of connections between the wall and the supports, and the vertical and horizontal spacing of the support system [16]. In the FE model, the two trench steel shields are stacked upon each other to cover the total depth (6 m (20 ft)) of excavation. The geometrical dimensions and mechanical properties are given in Tables 1 and 2, and the positions and locations of the struts used for simulation are shown in Figure 2. To account for the “sandwich-type” construction of the trench box plates, unit weight (γ) of the trench box plates was calculated from the real weight of the trench box used for experimental test, divided by its real volume. The real thickness of sandwich plates was the PLAXIS-3D input for the plate geometry. Four different thicknesses were considered in the parametric study: 200 mm (8 in.), 150 mm (6 in.), 100 mm (4 in.), and 75 mm (3 in.). For the struts, cross-sectional area (A) and moment of inertia (I) were calculated from the original strut diameter, while Young’s modulus (E) and unit weight (γ) were taken from the SCH80 chart. The vertical positions of the struts were kept different on both sides of the steel trench box as per the in situ experimental test at the Louiseville site. Note that 200 mm (8 in.) strut diameter and 100 mm (4 in.) shield plate were used for till, dry sand, wet sand, and sensitive clay soil analysis, whereas other dimensions were used for geometric influence FE analysis.Table 1
Plate properties (steel and aluminum) of trench box shield protection used for FE simulation in PLAXIS-3D.
Parameter/propertyPLAXIS nomenclatureSteel plateAluminum plate200 mm (8 in.)150 mm (6 in.)100 mm (4 in.)75 mm (3 in.)200 mm (8 in.)150 mm (6 in.)100 mm (4 in.)75 mm (3 in.)Material typeTypeElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicYoung’s modulusE1 (kN/m2)2×1082×1082×1082×1087×1077×1077×1077×107Unit weightγ (kN/m3)11.5911.8613.6715.702.683.224.294.99Thicknessd (m)0.200.150.100.080.200.150.100.08Poisson’s ratioν120.300.300.300.300.300.300.300.30Shear modulusG12 (kN/m2)8×1078×1078×1078×1072.7×1072.7×1072.7×1072.7×107G13 (kN/m2)8×1078×1078×1078×1072.7×1072.7×1072.7×1072.7×107G23 (kN/m2)8×1078×1078×1078×1072.7×1072.7×1072.7×1072.7×107One shield geometry;length=3 m (10 ft); height=3 m (10 ft); horizontal distance between struts=2.6416 m (104 in.).Table 2
Strut properties (steel and aluminum) of trench box shield protection used for FE modelling in PLAXIS-3D.
Parameter/propertyPLAXIS nomenclatureSteel strutAluminum strut200 mmØ (8 in. Ø)100 mmØ (4 in. Ø)50 mmØ (2 in. Ø)200 mmØ (8 in. Ø)100 mmØ (4 in. Ø)50 mmØ (2 in. Ø)Material typeTypeElasticElasticElasticElasticElasticElasticYoung’s modulusE1 or E2 (kN/m2)2×1082×1082×1087×1077×1077×107Unit weightγ (kN/m3)78.5078.5078.5026.9926.9926.99Cross-sectional areaA (m2)0.010.00280.00100.010.00280.001Moment of inertiaI2xx (m4)4.41×10−53.99×10−63.67×10−74.41×10−53.99×10−63.67×10−7I3yy (m4)4.41×10−53.99×10−63.67×10−74.41×10−53.99×10−63.67×10−7Poisson’s ratioν120.30.30.30.30.30.3200 mm strut (SCH80)Outer diameter strut, do (mm)2151126021511260Inner diameter strut, di (mm)19095481909548Area of strut (m2)8.24×10−32.83×10−39.64×10−48.24×10−32.83×10−39.64×10−4Iyy or Ixx for strut (m4)4.41×10−53.99×10−63.67×10−74.41×10−53.99×10−63.67×10−7One shield geometry;length=3 m (10 ft); height=3 m (10 ft); horizontal distance between struts=2.6416 m (104 in.).Figure 2
Elevation views and strut locations in the two assembled box shields (stacked upon each other) along the trench.In the PLAXIS-3D simulation, the steel trench box walls were modelled as “plate elements” and the struts as “beam elements.” “Hinge joints” were used between the two box shields (stacked upon each other) to represent the experimental setup (Figures3). The joints between struts and plates were also simulated as “hinge joints” as per the experimental setup (Figures 3). In the in situ experiment, when the 1st box is inserted into the soil and the upper box (2nd box) is pushed on top of the lower box (1st box) using an excavator, the soil provide the relevant lateral support for the bottom end of the lower trench box (1st box). Alam et al. [27] noticed in their experimental work that the bottom end part of the shield (approximately 300 mm) was embedded into the soil. Therefore, this part of the lower box acted more likely as a “fixed-end support” and was represented by a “fixed-end anchored element” in the FE model (Figure 3). The pressure measuring locations are shown in Figure 4. To analyze the nonlinear stress-strain behaviour of the soil, a Mohr-Coulomb (MC) constitutive material model was chosen in this study.Figure 3
Structural simulation view of the two assembled box shields (stacked upon each other) with protection in PLAXIS-3D (number of structural elements: 1050; MC soil model).Figure 4
Trench box supporting wall side views and pressure measuring locations from the top.


### 3.3. Generating the Mohr-Coulomb (MC) Model

The Mohr-Coulomb (MC) model is a simple and widely used linear elastic perfectly plastic model. Tables3 and 4 shows soil properties for the sensitive blue clay (Louiseville clay), till, dry sand, and wet sand model parameters used for FE modelling in PLAXIS-3D. The soil stiffness modulus (E) at different depths for sensitive blue clay (Louiseville clay) was calculated from the shear strength (Cu) using the equation proposed by R. Peck, The Selection of Soil Parameters for the Design ofFoundations. Second Nabor Carrillo Lecture, Mexican Societyfor Soil Mechanics, Mexico, 1975.:
(1)E=600∗Cu.Table 3
Soil properties for sensitive blue clay, till, and dry and wet sand model parameters used for FE modelling in PLAXIS-3D.
Parameter/propertyPLAXIS nomenclatureSensitive blue clayTillDry sand, medium densityWet sand, medium densityDepth of layerm0.0–150.0–150.0–150.0–15Material modelModelMCMCMCMCDrainage typeTypeUndrainedDrainedDrainedDrainedUnit wt. above phreatic levelγunsat (kN/m3)1621.519.519.5Unit wt. below phreatic levelγsat (kN/m3)17——19.5Maximum elastic/shear modulusEavg/Gmax (kN/m2)(Eavg as per Equation (1))(Gmax as per Equation (2))(Gmax as per Eq. (2))(Gmax as per Equation (2))K2max——705252CohesionCref′ (kN/m2)See Table4555Frictional angleϕ (° (degree))See Table4423535Dilatancy angleΨ (° (degree))0000Poisson’s ratioνur′0.30.30.30.3Interface strength reduction factorRinter0.50.650.650.65InitialK0 determinationK00.531.00.430.43Soil typeVery fine<2μm%741010102μm–50 μm%11131313Table 4
Parameters for different layers of sensitive blue clay (Louiseville clay, data from [34]).
Soil layerParameters in total stressesDrained conditionFissured brown clay (0 to 0.6 m)c′=5 kPa and ϕ′=28°DrainedPlastic brown clay (0.6 to 2 m)c=cu (Figure 1(a)) and ϕ=0°UndrainedSensitive blue clay (2 to 15 m)c=cu (Figure 1(a)) and ϕ=0°UndrainedUsing Equation (1), E was calculated for each layer, and the average values of E (Eavg) were taken as PLAXIS-3D input for the 1 m soil layers (further details are provided in [19]). This method to evaluate E is practical as an in situ shear strength testing and is very common in engineering practice, allowing thereby to practicing engineers to have access to those representative geotechnical properties. The total stress soil properties for sensitive blue clay (Louiseville clay) are presented in Table 3. For till, dry sand, and wet sand, which have different unit weights (γ) of sand and angles of friction (φ), these soil parameters for the MC model were used for parametric FE modelling, as shown in Table 3. For the different sandy soil stiffness moduli, the maximum shear modulus, Gmax, was calculated for each meter of depth using Equation (2), developed by Seed and Idriss [28]:
(2)Gmax=1000K2maxσm′0.5,where σm′ is effective mean stress=σv′+2σh′/3; σv′ is vertical effective stress; σh′ is horizontal effective stress; σh′=K0×σv′; and K0 is coefficient of earth pressure at rest. K2max is an empirical factor that varies with void ratio. The values of K2max considered in this study for the different soil types are presented in Table 3. Then, these calculated maximum shear modulus values (Gmax) for till, dry sand, and wet sand were used separately for each type of soil for drainage analysis in PLAXIS-3D simulation models.Soil stiffness parameters vary with the depth. Therefore, anisotropy is considered for the soil models, which are subdivided into 1 m deep layers below the first 2 meters, and the parametersE and Gmax are computed for each layer using Equations (1) and (2), respectively. Unloading/reloading stiffness parameter is given by Gmax and considered elastic, as in Karlsrud and Andresen’s [3] study. Further details related to the Louiseville clay soil model can be found in Alam et al. [19].


### 3.4. Interface Model

Trench box steel plate forms a smooth wall, but in construction practice, it is impossible to find a smooth wall that has no roughness [29]. Interfaces were added for proper modelling of SSI. They represent a thin zone of shearing material at the contact between a plate element and the surrounding soil. Negative and positive interfaces were added on either side of a plate surface (Figure 5(c)). The interaction between plate elements and soil elements was expressed by a strength reduction factor (Rinter), which relates the soil strength (friction angle and cohesion) to the interface strength (wall friction and adhesion). According to PLAXIS manual, recommended values for Rinter are 2/3 for fissured brown clay and <2/3 for plastic brown and sensitive blue soft clay [20] (Table 3). The interface was assigned a virtual thickness that varied with mesh type and in this case was 103×10−3 m. The interface element consisted of pairs of nodes that were compatible with a six-node triangular side of a soil element or a plate element [20].Figure 5
Soil and structural mesh elements and interfaces in PLAXIS-3D simulation.
(a)
Soil mesh(b)
Structural mesh of trench box inside the soil(c)
InterfacesFor till and dry sand analysis, water was not considered (water level below the bottom of the excavation). For the wet sand, water table was considered at 2.5 m below the top of soil. Inside the trench, water was considered pumped out (to allow work inside the trench). Therefore, water was present in the surrounding soil of excavation, creating thereby a horizontal pressure on the wall of trench box shield outside the trench in addition to soil pressure.


### 3.5. Mesh Generation

Medium coarse 10-noded elements were used Figures5(a) and 5(b). The mesh has a total of 24271 elements, i.e., 23221 soil elements and 1050 structural elements.


## 3.1. Soil Modelling

To simulate the trench, the soil model geometry ((7.5m25ft×12m40ft×15m50ft) was more than twice the size of the excavation contour (3m10ft×1.5m5ft×6m20ft) following recommendations by Hsieh and Ou [25] for shallow trench excavations (Figure 1(b)). The soil layers were modeled as follows: (i) fissured brown clay from 0 m to 0.6 m (2 ft) depth; (ii) plastic brown clay from 0.6 m (2 ft) to 2 m (6.56 ft) depth; and (iii) sensitive blue clay (Louiseville clay) from 2 m (6.56 ft) to 15 m (50 ft) depth, in accordance with the stratigraphic profile of the Louiseville soil [19]. Soil strength properties were selected carefully as they are a key factor in ESS performance [26]. The in situ value of undrained shear strength shown in Figure 1(a) was used in this FEA. On the other hand, till, dry sand, and wet sand layers were modelled using 0.0 to 15 m depths of the properties of each type of sand, as shown in Figure 1(c), with varying stiffness parameters along the depth.Figure 1
Geometric profile and excavation contour in PLAXIS-3D simulation of trench: (a) experimental values of undrained shear strength of sensitive blue clay (Louiseville clay, data from [33]); (b) PLAXIS-3D model profile for sensitive blue clay (number of soil elements: 23221; MC soil model); (c) PLAXIS-3D model profile for till and dry and wet sand (number of soil elements: 23221; MC soil model).
(a)(b)(c)


## 3.2. Structural Modelling

The stiffness of an ESS is a function of the flexural rigidity of the wall element, the structural stiffness of the support elements, the type of connections between the wall and the supports, and the vertical and horizontal spacing of the support system [16]. In the FE model, the two trench steel shields are stacked upon each other to cover the total depth (6 m (20 ft)) of excavation. The geometrical dimensions and mechanical properties are given in Tables 1 and 2, and the positions and locations of the struts used for simulation are shown in Figure 2. To account for the “sandwich-type” construction of the trench box plates, unit weight (γ) of the trench box plates was calculated from the real weight of the trench box used for experimental test, divided by its real volume. The real thickness of sandwich plates was the PLAXIS-3D input for the plate geometry. Four different thicknesses were considered in the parametric study: 200 mm (8 in.), 150 mm (6 in.), 100 mm (4 in.), and 75 mm (3 in.). For the struts, cross-sectional area (A) and moment of inertia (I) were calculated from the original strut diameter, while Young’s modulus (E) and unit weight (γ) were taken from the SCH80 chart. The vertical positions of the struts were kept different on both sides of the steel trench box as per the in situ experimental test at the Louiseville site. Note that 200 mm (8 in.) strut diameter and 100 mm (4 in.) shield plate were used for till, dry sand, wet sand, and sensitive clay soil analysis, whereas other dimensions were used for geometric influence FE analysis.Table 1
Plate properties (steel and aluminum) of trench box shield protection used for FE simulation in PLAXIS-3D.
Parameter/propertyPLAXIS nomenclatureSteel plateAluminum plate200 mm (8 in.)150 mm (6 in.)100 mm (4 in.)75 mm (3 in.)200 mm (8 in.)150 mm (6 in.)100 mm (4 in.)75 mm (3 in.)Material typeTypeElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicElastic, isotropicYoung’s modulusE1 (kN/m2)2×1082×1082×1082×1087×1077×1077×1077×107Unit weightγ (kN/m3)11.5911.8613.6715.702.683.224.294.99Thicknessd (m)0.200.150.100.080.200.150.100.08Poisson’s ratioν120.300.300.300.300.300.300.300.30Shear modulusG12 (kN/m2)8×1078×1078×1078×1072.7×1072.7×1072.7×1072.7×107G13 (kN/m2)8×1078×1078×1078×1072.7×1072.7×1072.7×1072.7×107G23 (kN/m2)8×1078×1078×1078×1072.7×1072.7×1072.7×1072.7×107One shield geometry;length=3 m (10 ft); height=3 m (10 ft); horizontal distance between struts=2.6416 m (104 in.).Table 2
Strut properties (steel and aluminum) of trench box shield protection used for FE modelling in PLAXIS-3D.
Parameter/propertyPLAXIS nomenclatureSteel strutAluminum strut200 mmØ (8 in. Ø)100 mmØ (4 in. Ø)50 mmØ (2 in. Ø)200 mmØ (8 in. Ø)100 mmØ (4 in. Ø)50 mmØ (2 in. Ø)Material typeTypeElasticElasticElasticElasticElasticElasticYoung’s modulusE1 or E2 (kN/m2)2×1082×1082×1087×1077×1077×107Unit weightγ (kN/m3)78.5078.5078.5026.9926.9926.99Cross-sectional areaA (m2)0.010.00280.00100.010.00280.001Moment of inertiaI2xx (m4)4.41×10−53.99×10−63.67×10−74.41×10−53.99×10−63.67×10−7I3yy (m4)4.41×10−53.99×10−63.67×10−74.41×10−53.99×10−63.67×10−7Poisson’s ratioν120.30.30.30.30.30.3200 mm strut (SCH80)Outer diameter strut, do (mm)2151126021511260Inner diameter strut, di (mm)19095481909548Area of strut (m2)8.24×10−32.83×10−39.64×10−48.24×10−32.83×10−39.64×10−4Iyy or Ixx for strut (m4)4.41×10−53.99×10−63.67×10−74.41×10−53.99×10−63.67×10−7One shield geometry;length=3 m (10 ft); height=3 m (10 ft); horizontal distance between struts=2.6416 m (104 in.).Figure 2
Elevation views and strut locations in the two assembled box shields (stacked upon each other) along the trench.In the PLAXIS-3D simulation, the steel trench box walls were modelled as “plate elements” and the struts as “beam elements.” “Hinge joints” were used between the two box shields (stacked upon each other) to represent the experimental setup (Figures3). The joints between struts and plates were also simulated as “hinge joints” as per the experimental setup (Figures 3). In the in situ experiment, when the 1st box is inserted into the soil and the upper box (2nd box) is pushed on top of the lower box (1st box) using an excavator, the soil provide the relevant lateral support for the bottom end of the lower trench box (1st box). Alam et al. [27] noticed in their experimental work that the bottom end part of the shield (approximately 300 mm) was embedded into the soil. Therefore, this part of the lower box acted more likely as a “fixed-end support” and was represented by a “fixed-end anchored element” in the FE model (Figure 3). The pressure measuring locations are shown in Figure 4. To analyze the nonlinear stress-strain behaviour of the soil, a Mohr-Coulomb (MC) constitutive material model was chosen in this study.Figure 3
Structural simulation view of the two assembled box shields (stacked upon each other) with protection in PLAXIS-3D (number of structural elements: 1050; MC soil model).Figure 4
Trench box supporting wall side views and pressure measuring locations from the top.


## 3.3. Generating the Mohr-Coulomb (MC) Model

The Mohr-Coulomb (MC) model is a simple and widely used linear elastic perfectly plastic model. Tables3 and 4 shows soil properties for the sensitive blue clay (Louiseville clay), till, dry sand, and wet sand model parameters used for FE modelling in PLAXIS-3D. The soil stiffness modulus (E) at different depths for sensitive blue clay (Louiseville clay) was calculated from the shear strength (Cu) using the equation proposed by R. Peck, The Selection of Soil Parameters for the Design ofFoundations. Second Nabor Carrillo Lecture, Mexican Societyfor Soil Mechanics, Mexico, 1975.:
(1)E=600∗Cu.Table 3
Soil properties for sensitive blue clay, till, and dry and wet sand model parameters used for FE modelling in PLAXIS-3D.
Parameter/propertyPLAXIS nomenclatureSensitive blue clayTillDry sand, medium densityWet sand, medium densityDepth of layerm0.0–150.0–150.0–150.0–15Material modelModelMCMCMCMCDrainage typeTypeUndrainedDrainedDrainedDrainedUnit wt. above phreatic levelγunsat (kN/m3)1621.519.519.5Unit wt. below phreatic levelγsat (kN/m3)17——19.5Maximum elastic/shear modulusEavg/Gmax (kN/m2)(Eavg as per Equation (1))(Gmax as per Equation (2))(Gmax as per Eq. (2))(Gmax as per Equation (2))K2max——705252CohesionCref′ (kN/m2)See Table4555Frictional angleϕ (° (degree))See Table4423535Dilatancy angleΨ (° (degree))0000Poisson’s ratioνur′0.30.30.30.3Interface strength reduction factorRinter0.50.650.650.65InitialK0 determinationK00.531.00.430.43Soil typeVery fine<2μm%741010102μm–50 μm%11131313Table 4
Parameters for different layers of sensitive blue clay (Louiseville clay, data from [34]).
Soil layerParameters in total stressesDrained conditionFissured brown clay (0 to 0.6 m)c′=5 kPa and ϕ′=28°DrainedPlastic brown clay (0.6 to 2 m)c=cu (Figure 1(a)) and ϕ=0°UndrainedSensitive blue clay (2 to 15 m)c=cu (Figure 1(a)) and ϕ=0°UndrainedUsing Equation (1), E was calculated for each layer, and the average values of E (Eavg) were taken as PLAXIS-3D input for the 1 m soil layers (further details are provided in [19]). This method to evaluate E is practical as an in situ shear strength testing and is very common in engineering practice, allowing thereby to practicing engineers to have access to those representative geotechnical properties. The total stress soil properties for sensitive blue clay (Louiseville clay) are presented in Table 3. For till, dry sand, and wet sand, which have different unit weights (γ) of sand and angles of friction (φ), these soil parameters for the MC model were used for parametric FE modelling, as shown in Table 3. For the different sandy soil stiffness moduli, the maximum shear modulus, Gmax, was calculated for each meter of depth using Equation (2), developed by Seed and Idriss [28]:
(2)Gmax=1000K2maxσm′0.5,where σm′ is effective mean stress=σv′+2σh′/3; σv′ is vertical effective stress; σh′ is horizontal effective stress; σh′=K0×σv′; and K0 is coefficient of earth pressure at rest. K2max is an empirical factor that varies with void ratio. The values of K2max considered in this study for the different soil types are presented in Table 3. Then, these calculated maximum shear modulus values (Gmax) for till, dry sand, and wet sand were used separately for each type of soil for drainage analysis in PLAXIS-3D simulation models.Soil stiffness parameters vary with the depth. Therefore, anisotropy is considered for the soil models, which are subdivided into 1 m deep layers below the first 2 meters, and the parametersE and Gmax are computed for each layer using Equations (1) and (2), respectively. Unloading/reloading stiffness parameter is given by Gmax and considered elastic, as in Karlsrud and Andresen’s [3] study. Further details related to the Louiseville clay soil model can be found in Alam et al. [19].


## 3.4. Interface Model

Trench box steel plate forms a smooth wall, but in construction practice, it is impossible to find a smooth wall that has no roughness [29]. Interfaces were added for proper modelling of SSI. They represent a thin zone of shearing material at the contact between a plate element and the surrounding soil. Negative and positive interfaces were added on either side of a plate surface (Figure 5(c)). The interaction between plate elements and soil elements was expressed by a strength reduction factor (Rinter), which relates the soil strength (friction angle and cohesion) to the interface strength (wall friction and adhesion). According to PLAXIS manual, recommended values for Rinter are 2/3 for fissured brown clay and <2/3 for plastic brown and sensitive blue soft clay [20] (Table 3). The interface was assigned a virtual thickness that varied with mesh type and in this case was 103×10−3 m. The interface element consisted of pairs of nodes that were compatible with a six-node triangular side of a soil element or a plate element [20].Figure 5
Soil and structural mesh elements and interfaces in PLAXIS-3D simulation.
(a)
Soil mesh(b)
Structural mesh of trench box inside the soil(c)
InterfacesFor till and dry sand analysis, water was not considered (water level below the bottom of the excavation). For the wet sand, water table was considered at 2.5 m below the top of soil. Inside the trench, water was considered pumped out (to allow work inside the trench). Therefore, water was present in the surrounding soil of excavation, creating thereby a horizontal pressure on the wall of trench box shield outside the trench in addition to soil pressure.


## 3.5. Mesh Generation

Medium coarse 10-noded elements were used Figures5(a) and 5(b). The mesh has a total of 24271 elements, i.e., 23221 soil elements and 1050 structural elements.


## 4. Construction Steps and Phases in PLAXIS-3D Simulation for Sensitive Blue Clay

Construction sequence has a significant impact on the lateral loads acting on a flexible retaining structure [30]. The model includes a step-by-step excavation sequence, including the structural element insertion into the excavation. The excavation sequence is summarized in Table 5, and construction and excavation steps are shown in Figure 6 for sensitive blue clay.Table 5
Construction sequence of FE modelling for sensitive blue clay in PLAXIS-3D.
PhaseIn the modelDrainagePhase 0Only soil elements (no structure)K0 procedurePhase 1Excavation (1st layer, depth: 0.6 m) fissured brown clayDrainedPhase 2Excavation (2nd layer: 0.6 to 2 m) brown clayPlastic: undrained B1 (soil model)Phase 3Adding structural elementPlastic: undrained B1 (soil model)Phase 4Excavation (3rd layer: 2 m to 6 m) sensitive blue clayPlastic: undrained B1 (soil model)1In the undrained (B) drainage type, modelling of undrained behaviour using effective parameters for stiffness and undrained shear strength parameters is allowed.Figure 6
PLAXIS-3D simulation of the construction steps.
(a)
Step 0 (initial phase)(b)
Step 1 (Phase 1)(c)
Step 2 (Phase 2)(d)
Step3 (Phase 3)(e)
Step 4 (Phase 4)


### 4.1. Calculations with DifferentK0

In practice,K0 for a drained soil is computed with Jaky’s [31] empirical expression. In an overconsolidated soil, K0 would be expected to be larger than the value given by Jaky’s formula. Hamouche et al. [32] evaluated the in situ coefficient of earth pressure at rest for Louiseville sensitive blue clay. Three values of K0 (0.55, 1, and 1.5) were selected to evaluate the effect on the soil pressure on the trench box as follows: (i) K0=0.55 represents a common value for normally consolidated eastern Canada clays (with ϕ′=27°); (ii) K0=1 represents a clay in undrained conditions (when ϕ=0); and (iii) K0=1.5 was determined experimentally at 6 m depth by Hamouche et al. [32] for Louiseville clay. Figure 7 shows the value of soil pressure computed with different K0 values using the PLAXIS model.Figure 7
Soil pressure at different depths of trench with differentK0 values.


## 4.1. Calculations with DifferentK0

In practice,K0 for a drained soil is computed with Jaky’s [31] empirical expression. In an overconsolidated soil, K0 would be expected to be larger than the value given by Jaky’s formula. Hamouche et al. [32] evaluated the in situ coefficient of earth pressure at rest for Louiseville sensitive blue clay. Three values of K0 (0.55, 1, and 1.5) were selected to evaluate the effect on the soil pressure on the trench box as follows: (i) K0=0.55 represents a common value for normally consolidated eastern Canada clays (with ϕ′=27°); (ii) K0=1 represents a clay in undrained conditions (when ϕ=0); and (iii) K0=1.5 was determined experimentally at 6 m depth by Hamouche et al. [32] for Louiseville clay. Figure 7 shows the value of soil pressure computed with different K0 values using the PLAXIS model.Figure 7
Soil pressure at different depths of trench with differentK0 values.


## 5. Construction Steps and Phases in PLAXIS-3D Simulation for Till, Dry Sand, and Wet Sand

To simulate till and dry sand, the water table was maintained below excavation level, and 2.5 m below the ground for wet sand, which was the lowest level recorded in experimental testing [19]. Then, the PLAXIS-3D simulation model was implemented to perform the analysis. However, when the cohesion value C=0 kPa was used, the PLAXIS-3D simulation showed that the excavation steps failed because the soil bodies failed. Hence, trial methods were used, increasing cohesion values in the range of 4–5 kPa. Thereafter, the excavation steps proceeded with the cohesion value of 5 kPa, and the model was run for analysis. However, the numerical analysis results showed that the soil pressure on the box shield was notably low, and at a certain depth, the soil ceased to be in an elastic state.When only the initial phase was used (meaning no excavation steps, installing the structure inside the soil straight away, and then removing the soil inside the trench box), the soil stresses were more reliable and comparable with the Terzaghi and Peck [7] apparent earth pressure curve for granular soil. It was observed that till created more stress on the shield than any other soil type. Dry sand yielded a higher soil pressure than wet sand except at 5.4 m depth, as shown in Figure 8. All the FEA results were compared with the Terzaghi and Peck [7] curves for till, dry and wet sand, and sensitive blue clay (Louiseville clay), as shown in Figure 8.Figure 8
Comparisons of soil pressure for different types of soil by numerical analysis with theoretical values in the trench.


## 6. Parametric Studies for Steel and Aluminum Trench Box for Plate and Strut Geometry

For this type of parametric study, the M-C model for sensitive blue clay (Louiseville clay) was chosen for analysis. The objective was to evaluate the influence of material type and geometry on the earth pressure on the temporary excavation trench box shield. Three scenarios were considered: (i) trench box materials (steel and aluminum), (ii) steel trench box with different plate thicknesses and strut diameters, and (iii) aluminum trench box with different plate thicknesses and strut diameters. A total of 16 trench box configurations were considered (4thicknesses×2materials×2strut diameters). Table 6 shows the simulation matrix, which was used for the PLAXIS-3D simulations.Table 6
Simulation matrix for parametric study of the trench box.
Wall thickness of the boxStrut diameterSteel boxAluminum box75 mm (3 in.)200 mm (8 in.)Ø and 50 mm (2 in.) Ø200 mm (8 in.)Ø and 50 mm (2 in.) Ø100 mm (4 in.)200 mm (8 in.)Ø and 100 mm (4 in.) Ø200 mm (8 in.)Ø and 100 mm (4 in.) Ø150 mm (6 in.)200 mm (8 in.)Ø and 100 mm (4 in.) Ø200 mm (8 in.)Ø and 100 mm (4 in.) Ø200 mm (8 in.)200 mm (8 in.)Ø and 100 mm (4 in.) Ø200 mm (8 in.)Ø and 100 mm (4 in.) Ø1in.=25 mm.


## 7. Overall Results and Comparisons

The Terzaghi and Peck [7] curves for till, dry sand, and wet sand were compared with the FEA result for homogeneous granular soils (Figure 8). A detailed comparison of the sensitive blue clay (Louiseville clay) FEA results to theoretical apparent earth pressure is provided in Alam et al. [19]. It is evident that when the trench box material type and geometry changed, the earth pressure on the temporary shield also changed for the same soil type. The results are summarized in Figure 8 for different types of soil and in Table 7 and Figures 9–12 for different materials and geometries. Based on these tables and figures, the following comparative observations can be made:
(1)
Till, dry sand, and wet sand versus theoretical predictions: the “Peck” curve was calculated from apparent earth pressure according to the Terzaghi and Peck [7] method (TPM). From the “Peck” curve, values for till, dry sand, and wet sand were overestimated in the upper part compared to the deeper part of the trench(2)
Steel versus aluminum trench box shield: under the same soil condition (Louiseville clay), a significant difference in earth pressures between the steel and aluminum temporary shields was observed. The maximum undrained total stresses decreased by 40.6% when shifting from a steel to an aluminum trench box shield with the same geometry of 100 mm (4 in.) plate thickness and 200 mm (8 in.) strut diameter(3)
Steel versus aluminum trench box shield plate thickness: reduction of shield plate thickness influenced soil pressure on both steel and aluminum trench box shields under the same soil condition (Louiseville clay). The FEA results showed that changing the plate thickness from 200 mm (8 in.) to 150 mm (6 in.), 150 mm (6 in.) to 100 mm (4 in.), and 100 mm (4 in.) to 75 mm (3 in.) resulted in a maximum decrease in soil pressure on the steel trench box of 5.5%, 9.2%, and 17%. These values reached 35.11%, 33.67%, and 13.5% for aluminum trench box for the same 200 mm (8 in.) strut diameter(4)
Steel versus aluminum trench box shield strut diameter (Louiseville clay): once again, a reduction in shield strut diameter influenced soil pressure on both steel and aluminum shields under the same soil condition. Furthermore, FEA results showed that changing the strut diameter from 200 mm (8 in.) to 100 mm (4 in.) gave a maximum of 1.12%, 1.9%, and 8.4% decrease in soil pressure on the steel trench box and 2.7%, 4.0%, and 9.5% on the aluminum trench box shield, respectively, for 200 mm (8 in.), 150 mm (6 in.), and 100 mm (4 in.) plate thickness for the same soil conditionTable 7
Soil pressure on steel and aluminum protection shield boxes with different strut diameters and plate thicknesses.
Steel trench box (kPa)Aluminum trench box (kPa)Depth from top (m)100 mm plate, 200 mmØ strut75 mm plate, 200 mm Ø strut100 mm plate, 200 mmØ strut75 mm plate, 50 mmØ strut150 mm plate, 200 mmØ strut150 mm plate, 100 mmØ strut200 mm plate, 200 mmØ strut200 mm plate, 100 mmØ strut100 mm plate, 200 mmØ strut75 mm plate, 200 mmØ strut100 mm plate, 100 mmØ strut75 mm plate, 50 mmØ strut150 mm plate, 200 mmØ strut150 mm plate, 100 mmØ strut200 mm plate, 200 mmØ strut200 mm plate, 100 mmØ strut0.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.00.93.04.54.96.52.02.30.80.95.805.396.786.244.64.93.22.82.439.233.234.829.237.436.938.738.428.9024.0426.2619.8133.632.536.635.93.951.038.344.933.456.255.159.558.830.2920.0927.3516.3446.744.853.952.55.479.465.572.763.681.280.683.883.459.5349.3758.4648.2273.972.979.879.0Figure 9
Soil pressure (Louiseville clay) on steel versus aluminum trench box shields with the same geometry.Figure 10
Soil pressure (Louiseville clay) on steel and aluminum trench box shields with different plate thicknesses.
(a)
Soil pressure on steel box shield of 200 mmØ struts(b)
Soil pressure on aluminum box shield of 200 mmØ struts(c)
Soil pressure on steel box shield of 100 mmØ struts(d)
Soil pressure on aluminum box shield of 100 mmØ strutsFigure 11
Soil pressure (Louiseville clay) on steel trench box shields with different plate thicknesses and strut diameters.
(a)
Change of strut diameter in steel box shield with 75 mm wall(b)
Change of strut diameter in steel box shield with 100 mm wall(c)
Change of strut diameter in steel box shield with 150 mm wall(d)
Change of strut diameter in steel box shield with 200 mm wallFigure 12
Soil pressure (Louiseville clay) on aluminum trench box shields with different plate thicknesses and strut diameters.
(a)
Change of strut diameter in aluminum box shield with 75 mm wall(b)
Change of strut diameter in aluminum box shield with 100 mm wall(c)
Change of strut diameter in aluminum box shield with 150 mm wall(d)
Change of strut diameter in aluminum box shield with 200 mm wallFrom the above discussions, it is apparent that soil pressure on the trench box shield depends not only on soil type but also on soil-structure stiffness. For the same material type and geometry, a trench box experiences different scales of pressure for different soil types. On the other hand, a given type of soil trench box experiences different scales of pressure for different types of material and geometry. Wall deformation shows some bulging as a result of flexure, and the level of bulging depends on the stiffness of the wall support system. A wall with a higher rigidity experiences less deformation and creates greater apparent earth pressure in a particular homogeneous soil. The results of these nonlinear parametric numerical models, in which only the wall (or trench box) stiffness (steel to aluminum) varied, show that the shape of the earth pressure distribution evolves from higher to lower values (i.e., 40.6% decrease of undrained total stresses) under a given set of soil conditions (i.e., the higher the wall rigidity, the greater the apparent earth pressure). Therefore, it is evident that soil pressure depends not only on soil stiffness but also on the stiffness of the trench box. Classical theories overestimate soil pressure in the upper part of the trench and may not be safe for the deeper part. When designing a trench box shield, the classical apparent earth pressure envelope may not be the only source of information. Adding in situ measurements and numerical data would be a better idea.


## 8. Conclusions

The parametric FE studies presented in this paper involve simulation of two trench box shields stacked upon each other to cover the total depth of the 6 m (20 ft) trench in sensitive clay and sandy soils, with the PLAXIS 3D FE computer code. The goal of these studies is to gain an insightful picture of the effect of (i) soil parameters on the steel trench box shield for till, dry sand, wet sand, and sensitive clay soil and (ii) the trench box material type (steel or aluminum) and geometry on the earth pressure. The FEM analysis results are also compared to published empirical apparent earth pressure diagrams.Based on the FEA results, the following observations can be made:(i)
Earth pressure on a steel trench box for till, dry sand, and wet sand seems to be underestimated by Terzaghi and Peck [7] earth pressure diagrams for the lower part of the trench box (below 4 m depth in the case of till and below 3 m depth in the case of dry and wet sand)(ii)
A significant difference in earth pressure was found between the steel and aluminum temporary shields. The maximum undrained total stresses decreased by 40.6% when shifting from a steel to an aluminum shield having the same geometry(iii)
FEA results showed that changing plate thickness from 100 mm (4 in.) to 75 mm (3 in.) resulted in a maximum 17% decrease in soil pressure on steel and a 33.67% decrease on aluminum trench box shields for the same type of soil(iv)
Changing the strut diameter from 200 mm (8 in.) to 100 mm (4 in.) resulted in a maximum 8.4% decrease in soil pressure on steel and a maximum 9.5% decrease on aluminum trench box shields for the same type of soil



---

*Source: 9949976-2021-09-21.xml*