 Open Access
 Total Downloads : 43
 Authors : Ammar Jalil Almosawi , Anis A. MohamadAli , David A. M Jawad
 Paper ID : IJERTV8IS060386
 Volume & Issue : Volume 08, Issue 06 (June 2019)
 Published (First Online): 19062019
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Behaviour of Waler Retaining System on Sheet Piles using Finite Element Method
Ammar Jalil Almosawi, Anis A. MohamadAli, David A. M Jawad
University of Kufa, Iraq
Abstract: In this paper, a truss retaining system on a cantilever sheet piles in a 300 meter length berth was studied using two dimensional plain strain analysis. The analysis of the soil continuum under the effect of active and passive pressure in addition to the cyclic load effect was extended by using the finite element method with the aid of the software Autodesk FEM ROBOT 2019. The most common analyses, designs and properties of such load cases were reviewed. The soil continuum was modeled using 15 node triangular element and the stress concentration for all the nodes was analyzed. The behavior of truss retaining system under the effect of active, passive and cyclic loads are predicted in this study by using FEM ROBOT software 2019.
Keywords: – Sheet piles, truss retaining system, FEM, ROBOT program

INTRODUCTION
Sheet piles are of the most important retaining systems. Behavior of the sheet piles depends mainly on the active load, the passive load and the way of retaining system. A truss retaining system with a continuous waling beam connected to sheet pile subjected to several load types like active, passive and cyclic load was analyzed in the present work. Usually the waler beam is connected to the tie rods in which the waler beam behaves as continuous walling beam and the tie rods points represent the reactions of the continuous waler beam. In the present work a truss supports the waler beam in a cantilever sheet piles. However, the key of our motivation is to look for an approach that can be employed to avoid excavation problems during the tie rods installing along the port length.
Reference is made to several papers that have discussed various types of supporting systems and their construction, for example, Chang de Young (2007), Bogdan Rymsza (2008), Leila Eskandari (2011), and David Baxter, (2016) [4], Engineers EDGE (2018). A truss system is installed in xzplane with I sections properties of (W28X16). (Uchannels U25X18) were installed as bracing in xz, yz and xyplanes for the truss. The layout of the truss is as shown in Figure (1). The effect of using a truss as supporting system instead of tie rods was tested using the ROBOT software 2019 program in order to study the stress concentration behavior, bending moment behavior and the effect of truss supporting instead of tie rods.
Figure (1) Side view of truss system retaining structure
Figure (1) shows a side view section of a berth 300 meter long with three lines of drive piles along the port. Figure (2) shows top view of the port with the located waler beam in front of the supporting truss.
Figure (2) Top view of truss system retaining structure

SITE DESCRIPTION

Layout of the structure
A berth of 300 meter long has two main parts, land part (18.64m) length and sea part (15.66m) was retained by a 30m sheet pile (AZ 22 800) supported by waler beam (double back to back channels C300X90), the waler beam is supported by a truss with specific properties installed in two directions (xz plane) between the driven piles as shown in Figure (1) instead of tie rod system in order to retain the active load pressure from the land part and any other load types like crane load in addition to wave loads or cyclic loads due to ship impact from the sea side. The reinforced concrete deck provides a partial anchoring system throughout the shear connecters welded on the sheet piles tip connected to the deck reinforcement.

Load application methodology
The effect of cyclic load on the behavior of waler beam is shown in Table (1), at each cycle there were two cases, the loading case and the unloading case, and this may give an overall view for the behavior of the continuous waler beam in order to conclude the proper results. The cyclic load applied on the sheet piles comes from the repeated sea waves or from ships impact and was taken into account as a cyclic load effect in the software input data. The cyclic load is not constant it is transmitted as cumulative waves acting on the sheet piles, then transmitted directly to the continuous waler beam. The values of the cyclic load [P] were ranged from (10.8 347.7 kN/mÂ²/cycle) as the minimum cyclic load recorded in the port was (10 kN/mÂ²/cycle) and the maximum cyclic load at failure was (365 kN/mÂ²/cycle).
The problem has been considered as plane strain problem so in order to improve the behavior of the waler reaction points, under the effect of cyclic load pressure it is required to analyze the results from tables by using zero order Hermite finite element polynomial function for 5 nodes element depending on the cubic equation function for plane strain analysis of continuous waler beam under the effect of cyclic load in finite element method in order to investigate the displacement of waler beam under cyclic load effect. The equivalent axial displacement for all load cases as shown from the following Hermite FEM equations (11) to (1 3) [5]
() = 1 + 2 + 32 + 43 (1.1)
The constants 1, 2, 3 4 can be found by applying the boundary conditions;
5 5
= 5 3 (1.2)
0 0
=1 =1
(+0.1 ) + (0.1 )
=
2
(1.3)
The finite element simulation for the transverse deformations of each node will be:
() = [] (1.4)
Where [N] is the shape functions for each node given by:
[] = [ 1 2 3 4] (1.5)The nodal deformations of the continuous waler beam element are shown in the Figure (4). The shape functions of the nodal deformations are calculated from equations (1.6 To 1.7) are as follows;
1() =
2() =
3() =
233 +3
2
2
3 (1.6)
2
2
32 +2
2 (1.7)
2
2
(233 )
3 2
3 2
3 (1.8)
() =
(1.9)
4 2
And the matrix of equivalent nodal displacement for the continuous waler beam element;
1
= [2]
3
4
The finite element method and according to simple beam theory, plane sections of the beam remain plane after deformation and
hence the axial displacement (u) due to the transverse displacement (v) can be expressed as it is shown in Figure (3);
Figure (3) Transverse displacement for continuous beam element [1]
The relation between the axial and the transverse deformation is shown in equation (1.10):
=
(1.10)
Where y represents the distance from the neutral axis of the distorted continuous waler beam. While the axial strain is given by;
2
2
= = = [] (1.11)
2
If the nodal displacements of the element, 1, 2, 3 4 are known, the stress distribution in the element, can be found as
= (, ) = = [][] (1.12) Where (, ) represents the stress in the element at a point located at a specific distance from the origin (in horizontal direction from the left node) and from the neutral axis (in the vertical direction).
Where [R] is given by;
[R] = {(12 6) (6 4) (12 6) (6 2)} (1.13)3
Using the equations above with [D] = [E], the stiffness matrix can be found as; [6]
[()] = [][][] = [][]()
0
Where:
=
12 6
[ 6 4212 6
p>6 22
12 6
6 22]
12 6
6 42
: is to be determined as modulus of subgrade reaction in the contact area
= 2 Is the area moment of inertia of the cross section about the z axis. Noting that the nodal interpolation functions
N,(x) are the same as the firstorder Hermite polynomials.
Many approaches exists in the finite element method to determine the contact area and the [k] matrix, which include the Rayleigh Ritz approach, variational or method of weighted residual which depends on scaling down the powers of the equations in order to find the constants.
The theory of plane strain is applied for bodies that are long like long pipes as mentioned previously and whose geometry and loading do not vary significantly in the longitudinal direction. Therefore the analysis of sheet piling, dams, cylinders, and retaining walls use the assumption of plane strain. In plane strain distribution, it is assumed that = 0 at every cross section. Here, the dependent variables are assumed to be functions of only the x and y coordinates provided. We consider a cross section of the
body away from the ends, in this case the two dimensional stressstrain relations given by:
Where,
= {} and
= [] + Â°
= {}
(1.14)
Then,
[] =1 +
[1 0
1 0]
0 0 2
1
= []( Â°) = [] 1 2 {1}
0
1 0
[] = [(1 + )(1 2)
1 0
1 2]
0 0
2
The component of stress in the z direction will be nonzero and is given by:
= ( + )
and,
= = 0
Where i represents the number of points at the continuous waler beam, k is the higher order of the cubic function chosen for solving the problem and represent the displacement at the middle point of the continuous waler beam.
In order to recognize the optimum critical load that makes the waler at critical point, the cyclic load is divided into portions ranges
from 0.1 to 1, which represents the whole cyclic load acting on the sheet pile. [3]
Table (1) refers to the displacements at five nodes indicated at the critical section from the continuous waler beam with four spans specified. Six sets of cyclic load were applied on two types of sheet piles with (30m) penetration depths. The five nodes where the displacements calculations were located on the top face of the back to back channel waler beam. See figure (3)
254 11875 Mpa
1
2
LOADING CASE
CYCLIC LOAD
Range
10.8 kN/ mÂ²
Load set No.1
Safe
LOADING CASE
CYCLIC LOAD
Range
10.8 kN/ mÂ²
Load set No.1
Safe
3
4 5
Primary affected Low affected Medium affected
Loaded low Loaded medium
High load
High concentration Affected stress
Figure (4) Reaction points on the waler
Table (1) shows the cyclic load applied for set number one as an example for all the cyclic load applied. The behavior of the waler beam at the five nodes represents the displacements of the waler beam under the effect of applied cyclic load. Figure (4) shows the load displacement curves for six load cases applied on the walling beam.
Table (1) displacements calculation of waler beam under set 1 loading
Set No. Set 1
Cyclic load 108.2
kN/mÂ²
Displacements measured in (mm)
Load portion (%)
1
10.8
+0.1P
0.17
0.45
0.85
0.46
0.17
0.86
0.1P
0.15
0.43
0.87
0.44
0.16
2
21.6
+0.2P
0.34
0.47
0.96
0.48
0.36
0.95
0.2P
0.34
0.48
0.94
0.49
0.35
3
32.4
+0.3P
0.85
0.95
1.25
0.94
0.84
1.265
0.3P
0.91
0.94
1.28
0.95
0.93
4
43.2
+0.4P
1.22
1.45
1.64
1.46
1.23
1.64
0.4P
1.24
1.46
1.64
1.45
1.25
5
54.1
+0.5P
1.96
1.77
2.66
1.73
1.97
2.675
0.5P
2.00
1.76
2.69
1.71
1.99
6
64.9
+0.6P
2.51
1.97
3.11
1.96
2.52
3.12
0.6P
2.54
1.96
3.13
1.97
2.54
7
75.7
+0.7P
2.92
2.25
6.96
2.24
2.91
6.95
0.7P
2.94
2.22
6.94
2.26
2.94
8
86.5
+0.8P
3.11
3.45
8.60
3.56
3.10
8.61
0.8P
3.11
3.89
8.62
3.65
3.91
9
97.3
+0.9P
3.50
4.12
10.93
4.13
3.57
10.95
0.9P
3.51
4.33
10.98
4.12
3.52
10
108.2
+P
4.12
4.78
18.15
4.88
4.13
18.17
P
4.12
4.68
18.19
4.69
4.15
400
400
300
load kN/ mÂ²
200
300
load kN/ mÂ²
200
Set 1 loading case Set 2 loading case
Set 3 loading case
0
Set 1 loading case Set 2 loading case
Set 3 loading case
0
Set 4 loading case
Set 4 loading case
25
25
50
Displacement mm
50
Displacement mm
75
75
100
100
Set 5 loading case
Set 6 loading case
Set 5 loading case
St 6 loading case
Figure (4) Load displacements curves for loading cases


CONCLUSIONS
The effect of using truss retaining system instead of tie rod and deadman system for long length berths could be more effective to withstand the cyclic load due to ships impact or sea waves. The use of a steel structure truss as retaining system decreases the probability of soil slip during excavation works. There are many factors that need to be considered in the analysis of the wailing like soil properties, sheet pile penetration depth, active and passive load pressure as well as the cyclic load effect then comparing the results.
REFERENCES

John Wiley 5th edition 2014 by I M Smith Programming The Finite Element Method University of Manchester U.K

Singiresu S. Rao 4th edition, 2004 The Finite Element Method In Engineering Elsevier Science & Technology Books

David V Hutton, 1st edition 2004, Fundamentals of Finite Element Analysis Mc Graw Hill Higher Education

David Baxter and others, 9th edition 2016 Piling Handbook ArcelorMittal Design & Execution manual

Wang Lunyan, Qian DaXing, and Zhang Jianwe 2010, Study on PY Curve Method for Computing Laterally Loaded Piles under Horizontal Distribution Loads The International Conference on Computer Application and System Modeling, V5366369.