Research Article  Open Access
Design of Jetty Piles Using Artificial Neural Networks
Abstract
To overcome the complication of jetty pile design process, artificial neural networks (ANN) are adopted. To generate the training samples for training ANN, finite element (FE) analysis was performed 50 times for 50 different design cases. The trained ANN was verified with another FE analysis case and then used as a structural analyzer. The multilayer neural network (MBPNN) with two hidden layers was used for ANN. The framework of MBPNN was defined as the input with the lateral forces on the jetty structure and the type of piles and the output with the stress ratio of the piles. The results from the MBPNN agree well with those from FE analysis. Particularly for more complex modes with hundreds of different design cases, the MBPNN would possibly substitute parametric studies with FE analysis saving design time and cost.
1. Introduction
Mooring dolphins are usually constructed when it would be impractical to extend the shore to provide access points to moor vessels. A typical mooring dolphin consists of a platform and several piles supporting the platform, which is socalled jetty. The vertical or battered piles are driven into the seabed. In design practice, deciding whether and where to use the vertical or battered piles is important issue. In the practical design process, the arrangement, the number, and the inclination of the piles are tentatively decided based on previous design experiences and then confirmed through finite element (FE) analysis. Therefore, building and analyzing lots of FE models adopting trial and error process are needed to find the optimum design.
Many researches have been performed to help designers to make decisions. An experimental study showed that the pile group effect is an important factor to resist horizontal loads [1]. For cyclic lateral loading, a zigzag arrangement shows higher resistance than an inline arrangement. Also it was shown that as the pile center distance increases, the stresses on the front piles decrease, while those on the rear piles increase [2]. When the center distance between piles becomes more than 3~5 times of pile diameter, the group effect decreases so that each pile can be considered as a single pile when the distance reaches 6 times the pile diameter [3]. The battered piles are commonly considered to resist lateral loads solely while the vertical piles resist gravity loads only. This traditional design assumption would make the design process easier but also it usually results in overestimated design. Moreover, it is well known that the vertical piles can also resist bending moments from the lateral loads. Through empirical studies, the py method has been proposed and developed by Kondner [4], Reese et al. [5], Scott [6], and Norris [7] to help the design of jetty structure. Though it is still a commonly adopted method, some concepts of the method are based on oversimplified or improper assumptions, especially in the effects of actual soil parameters after pile driving [8].
To overcome the complication of jetty pile design originated from mutual interaction among a number of design parameters, artificial neural networks (ANN) have been introduced in geotechnical engineering [9, 10]. This technique has also been applied successfully in static and dynamic pile systems [11, 12]. Kim et al. [8] predicted the lateral behavior of single and group piles using ANN and compared the results from ANN with the model test results. In this paper, as a suggestive solution of difficulties and cumbersome processes in building and analyzing lots of FE models, the ANN is adopted. Possibility of substituting ANN as a jetty structure analyzer for FE analysis is examined.
2. Methodology: Application of ANN as a Structure Analyzer
The jetty design process involves, as mentioned above, searching for the optimum pile pattern which results in the most effective pile usage within feasible design region. The internal forces of the piles of jetty structure subjected to horizontal mooring load vary unexpectedly depending on the inclination of the piles and deployment pattern of piles. Therefore, developing ANN, the input data to ANN are decided as horizontal load exerting on the jetty platform and the information of jetty piles, such as arrangement and inclination of piles, and the output results as the stress ratios of piles to confirm the feasibility of design candidate. Whole concept of methodology adopted in this paper is summarized in Figure 1. As shown in Figure 1, the trained ANN is used as a structural analyzer in this research, placing FE analysis. Firstly, the training samples are generated using FE analysis for various design conditions. It is important that the training samples should be generated from various available design conditions so that the trained ANN may predict adequately when it encounters real new design data. Also the number of training samples should be large enough to avoid overfitting. In this research, total of fifty design cases with different loading conditions and pile patterns are considered for generating training samples through FE analysis.
To construct the ANN architecture with predefined input and output layers, type of ANN, the number of hidden layers and neurons in each hidden layer, and type of transfer function for each layer should be determined. So in this research, because of the complexity of the problem, multilayer backpropagation neural network (MBPNN) with two hidden layers, shown in Figure 2, is adopted to tackle the problem. For the transfer function tangent sigmoid function and pure linear function are adopted for hidden layers and output layer, respectively, since the stress ratio, output from MBPNN, could be either compressive or tensional value. Each neuron of the hidden and output layers has bias and the neurons in one layer are interconnected with the neurons before and after the layer through weights.
Though the number of hidden layers of the MBPNN is determined as two, the performance of the MBPNN will vary depending on the number of neurons in hidden layers. Regarding the number of neurons in each hidden layer, however, there is no general rule to determine. Thus, in this study several neural network architectures with different number of neurons are examined for the best performance and generalization to new data based on the Kfold crossvalidation method [13–15]. The performance and generalization of the MBPNN are summarized as an average of root mean squared error (RMSE) from the Kfold crossvalidation. After fixing the architecture of MBPNN, training process is conducted to find the optimum values of the biases and weights using all fifty training samples. Finally, the trained MBPNN is used as a structural analyzer to produce the stress ratio of each pile for real design conditions.
3. Description of Jetty Structure for Analysis
In this study, a mooring dolphin which was designed for a real project is used. Variations of the pile layout which had been proposed from the early design stage were also considered. The main purpose of the project was to design and build a liquefied natural gas (LNG) terminal at a port area so that the gas product would be transmitted from floating storage and regasification unit (FSRU) to natural gas network onshore by pipelines. The dimension of the platform is 16 m in length, 10 m in width, and 2 m in thickness. The platform is made of reinforced concrete and its piles are made of steel. The plan and elevation view of the testbed mooring dolphin are shown in Figure 3.
3.1. Materials
The material properties of C35/45 concrete for the platform are shown in Table 1. In the latest European standard BS EN 2061 [16], the strength classes are classified using cylinder strength as well as a cube strength. S355 European standard steel is used for most of the structural members [17]. The material properties of S355 steel are shown in Table 2.


3.2. Load Conditions
The expected largest FSRU at the mooring dolphin has a capacity of 266,000 m^{3} and the largest LNG carrier has a capacity of 177,400 m^{3}. Maximum mooring force is calculated as 3750 kN. Dead loads are listed in Table 3. All permanent structural members as well as nonstructural members have been considered as dead loads on the structure. Nonstructural member (appurtenance) includes quick release hook (QRH), fender, handrail, and grating. Detail appurtenance loads are shown in Table 4. Pedestrian live load of 4.0 kN/m^{2} is assumed. The maximum wave height varies between 2.5 and 2.75 m during a year, but about 60% days of a year wave height is less than 0.5 m. The mean () wave period varies between 2.5 sec and 7 sec. Measurements about 4 km offshore indicate that the typical astronomical velocities are in the order of 0.5 m/s. The largest sea water current speed is 0.7 m/sec at the project site. The wind speed at the location is considered as 18 m/s. The maximum wind speed with 100 years of return period is 32.2 m/sec.


For the load combination, BS63492 [18] is adopted as shown: where , are partial factors, is prestressing, is leading variable action, and is combination factor.
3.3. Soil Conditions
The location of the virtual fixity points was computed by various methods: Chang’s method, AASHTO, Hansen’s method, and Lpile method. The pile penetration depth under the maximum tensile force was also computed by the Japanese bridge construction standard (2002), API recommended practice 2AWSD, AASHTO, and Broms’ analysis method. Based on those methods, it turned out that the penetration depth of 5 m into the bedrock would provide a fixed boundary condition at the bedrock level.
3.4. Design Configuration of Jetty Pile Pattern
Ten different configurations of jetty pile pattern are considered in this study depending on whether vertical or battered, if inclined, the batter direction, and the number of piles (Figure 4). The combination of ten different configurations and five different mooring forces (70%, 80%, 90%, 100%, and 110% of original mooring force) produced 50 FE models and they were analyzed to compute the stress ratios of piles. Among FE models, Patterns 1 and 2 are shown in Figure 5.
(a) Pattern 1 (proposed design)
(b) Pattern 2
(c) Pattern 3
(d) Pattern 4
(e) Pattern 5
(f) Pattern 6
(g) Pattern 7
(h) Pattern 8
(i) Pattern 9
(j) Pattern 10
(a)
(b)
4. Preliminary FE Analysis
After analyzing 50 FE models of mooring dolphins, the ratios of maximum stress to allowable stress of each pile were obtained under given loading condition as shown in Figure 6. Here, a pile pattern with the less number of required piles as well as with the smaller stress ratio is considered as the improved.
(a) Pattern 1 (proposed design)
(b) Pattern 2
(c) Pattern 3
(d) Pattern 4
(e) Pattern 5
(f) Pattern 6
(g) Pattern 7
(h) Pattern 8
(i) Pattern 9
(j) Pattern 10
4.1. Comparison of Battered and Vertical Piles
In Pattern 1, the absolute values of the stress ratio are all less than unity except Pile 5. Pile 5 shows the maximum compressive stress and Pile 8 shows the maximum tensile stress. Pattern 2, whose piles are all vertical, shows compressive stress and most of the stress ratios are greater than unity.
4.2. Patterns 3 and 4
When Pile 5 is absent (Pattern 3), the tensile stress of Pile 8 becomes bigger than that of the proposed design (Pattern 1). In Pattern 3, the compressive force on Pile 5 in Pattern 1 redistributes to the adjacent piles. Considering all stress ratios and the number of piles, Pattern 3 can be considered as more improved design than the proposed one. Pile 11 of Pattern 4 is compressive within 90% of applied force but it turns to be tensile when the mooring force is equal to or more than 100%. In this case, the absence of Pile 8 causes a rapid stress change in Pile 11 and design of reinforcing bars in concrete platform is difficult; therefore Pattern 4 shall be avoided.
4.3. Patterns 5, 6, and 7
When it is compared to Pattern 1, Pile 4 of Pattern 5 remains tensile. However, Piles 4 and 8 of Pattern 1 change to be compressive in Pattern 5, while Piles 5 and 11 become tensile. The use of vertical piles in Pattern 5 makes stress sign changed. All absolute values of stress ratio of Pattern 5 are less than unity. In this point of view, Pattern 5 is more effective than the proposed design. With absence of Pile 5, Pattern 6 shows slightly higher stress than Pattern 5. Similar to Pattern 4, the absence of Pile 8 causes rapid stress change in Pile 11 of Pattern 7. Pattern 7 shows higher compressive stress than Pattern 5.
4.4. Patterns 8, 9, and 10
The stress ratios of Patterns 8, 9, and 10 are similar to those of Patterns 5, 6, and 7 but they are not compatible for horizontal load direction change.
5. Architecture and Training of MBPNN
5.1. Training Samples
The first neuron of the input layer is assigned for load condition, and the other input neurons take the information of piles. To distinguish the battered and the vertical piles, the numbers “1” and “2” are assigned as neuron input values for each pile. To the location where pile is absent, the number “0” is assigned to the corresponding input neuron. For the value of mooring force corresponding to the value of input neuron 1 which is very big compared with the values of the other input neurons might lead to a failure of MBPNN training, the mooring force is normalized to the mooring force of 3750 kN. Table 5 shows the example of input values and corresponding pile locations.

In this study, MBPNN, with two hidden layers, utilizing backpropagation process is used. To obtain the training samples, five load cases—70%, 80%, 90%, 100%, and 110% of mooring force of 3750 kN—are applied to 10 different pile patterns of jetty structures. The combination of 10 jetty pile patterns and 5 load cases makes 50 training cases in total. Table 6 shows an example of input and target output of training samples.
(a)  
 
(b)  

5.2. Construction of MBPNN Architecture and Training
Since the performance and generalization of MBPNN to new design data will vary depending on the number of neurons in hidden layers, four different topologies of MBPNN with different number of neurons in hidden layers are examined: 13 (input neurons)15(1st hidden layer)15 (2nd hidden layer)12 (output layers), 13101012, 1371512, and (4) 13151012. The number of neurons in hidden layers of the first MBPNN model (13151512) is greater than that of input or output layers and vice versa in the second model (13101012). In this paper, the Kfold crossvalidation method is used to assess the generalization of model and to select the best architecture of MBPNN. For the Kfold crossvalidation the fifty training samples are randomly divided into 10 subsets, that is, 10fold crossvalidation. In the Kfold crossvalidation one subset is assigned as validation data set and the other nine subsets as training data set. Each MBPNN model is trained using the training data set and RMSE is computed for the validation data set. This procedure continues 10 times changing validation data set and training data set. Finally, a MBPNN model with the least averaged RMSE is selected as the best model.
The neural network toolbox provided by commercial program MATLAB was used to construct and for training of MBPNN models. The LevenbergMarquardt method with backpropagation process was adopted for the optimization algorithm for training. The LevenbergMarquardt method is considered to be effective for the complicated MBPNN for the fastest training [19, 20]. The objective function in the neural network training is defined as the minimization of mean squared error as shown: where total number of training pattern, = th target, = th input, and = weights and biases of neural network.
In the LevenbergMarquardt method, the optimum weights and biases are searched using where th weights and biases, = Jacobian matrix, , and .
The successful performance of the LevenbergMarquardt method depends on the choice of . Where the gradient is small, the search movement should be large so that the slow convergence is avoided. However, it should be small for the steeper gradient region. The initial value of was assumed as 0.001, and the increase and decrease factor of were assigned as 10 and 0.1, respectively. Thus, during the training, will take the value of , where is zero or natural number. The training process starts with the initial value and at the second step the objective function (i.e., mean squared error) is computed with the previous value of and . If both of these values do not result in good performance, a new value of is adopted for the next step. In the following steps, the best value is searched among “previous value,” “,” and “.” If a previous value results in reduction of the objective function, the value is not changed.
Table 7 summarizes the averaged RMSE of each MBPNN model from the Kfold crossvalidation. Interestingly the first model (13151512) which is the most complex one among the four models does not show the least averaged RMSE, rather than the second model of which the number of neurons in hidden layers is between that of input and output layer which shows the least. From the Kfold crossvalidation results, the second MBPNN model with 10 neurons for each hidden layer is selected as the best architecture.

After fixing the architecture of the MBPNN, the MBPNN was trained again with all the training samples. The training process was completed at 147 epochs, as shown in Figure 7, and terminated at performance goal of 10^{−6} which was set as one of the termination conditions. Figure 8 shows the gradient changes of the problem surface and changes of the value during the training process. If a value gives the reduction of objective function, the value is kept for the following steps. However, if it does not result in good performance, it is updated using the increase or decrease factor. To find the optimum point, the value is continuously updated during the training process and with the change of the value the searching direction and gradient of problem surface are changed. As long as the gradient is large enough to improve the training performance, the value is unchanged; however, when the gradient gets smaller and the training performance gets worse, which means the surface of the objective function becomes flat, the value is updated. Comparing the training performance graph and the graph, it is observed that the value starts from the initial value of 0.001 and is updated at epochs 2, 7, 16, 42, and 82. In Figure 8 with the constant values the performance and gradient become flat; however, the training performance and gradient are improved dramatically at those epochs.
5.3. Verification of Training
Cases 3, 15, 27, 39, and 50 were selected for verification of MBPNN training. Table 8 shows the comparison of 5 known targets of training cases and the simulated results from the trained MBPNN. The values in the last row are root mean squared error (RMSE) between the known targets and the simulated values. At the locations of no pile, the target value is 0 and the simulation results also show similar value and RMSE is very small. This ensures that the MBPNN is well trained.
 
RMSE: root mean squared error. ^{†}No pile at this location. 
6. Design with Trained MBPNN
Through the trained MBPNN, the stress ratios of jetty piles were obtained under different loading conditions which were not included in the training samples. The feasibility of the MBPNN was verified by comparing the results from FE model and the MBPNN. Table 9 shows the stress ratios computed by FE analysis and the MBPNN. The results from the MBPNN are very close to the FE analysis results. The RMSE is also very small regardless of the pile patterns.
 
RMSE: root mean squared error. ^{†}No pile at this location. 
7. Conclusions
In this paper, the application of MBPNN as a structural analyzer for jetty structures is explored. The framework of MBPNN is defined as the input with the lateral forces on the jetty structure and the type of piles and the output with the stress ratios of the piles. For the highly complex jetty pile patterns the results from the MBPNN show very good agreement with those from FE analysis. With the more training samples and the expansion of input parameters for jetty structure design, the MBPNN shows possibility to replace the repetitive and timeconsuming FE analysis. Although only 50 cases have been modeled for study purpose in this paper, the merit of MBPNN would be clearer as the number of cases increases.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
References
 J. S. Bae, T. Y. Kim, and S. H. Kim, “The lateral behavior characteristics of group piles in sand ground,” Journal of Korea Society of Civil Engineers C, vol. 20, no. 6, pp. 517–524, 2000. View at: Google Scholar
 J. S. Bae and S. H. Kim, “The behavior of group piles for pile arrangement under repeated lateral loads,” Journal of Korea Society of Civil Engineers, vol. 23, no. 4, pp. 231–239, 2003. View at: Google Scholar
 J. Y. Paek, J. Y. Cho, S. S. Jeong, and T. J. Hwang, “Shaft group efficiency of friction pile groups in deep soft clay,” Journal of Korea Society of Civil Engineers, vol. 32, no. 2, pp. 49–60, 2012. View at: Google Scholar
 R. L. Kondner, “Hyperbolic stressstrain response: cohesive soils,” Journal of the Soil Mechanics and Foundations Division, vol. 89, no. 1, pp. 115–143, 1963. View at: Google Scholar
 L. C. Reese, W. R. Cox, and F. D. Koop, “Analysis of laterally loaded piles in sand,” in Proceedings of the 6th Offshore Technology Conference, pp. 473–483, Dallas, Tex, USA, 1974. View at: Google Scholar
 R. F. Scott, “Analysis of centrifuge pile tests: simulation of pile driving,” Research Report OSAPR Project 13, American Petroleum Institute, 1980. View at: Google Scholar
 G. Norris, “Theoretically based BEF laterally loaded pile analysis,” in Proceedings of the 3rd International Conference on Numerical Methods in Offshore Piling, pp. 361–386, Navtes, 1986. View at: Google Scholar
 B. T. Kim, S. Kim, and S. H. Lee, “Prediction of lateral behavior of single and group piles using artificial neural networks,” KSCE Journal of Civil Engineering, vol. 5, no. 2, pp. 185–198, 2001. View at: Google Scholar
 W. T. Chan, Y. K. Chow, and L. F. Liu, “Neural network: an alternative to pile driving formulas,” Computers and Geotechnics, vol. 17, no. 2, pp. 135–156, 1995. View at: Publisher Site  Google Scholar
 G. W. Willis, C. Yao, R. Zhao, and D. Penumadu, “Stressstrain modeling of sands using artificial neural networks,” Journal of Geotechnical Engineering, vol. 121, no. 5, pp. 429–435, 1995. View at: Publisher Site  Google Scholar
 I. Lee and J. Lee, “Prediction of pile bearing capacity using artificial neural networks,” Computers and Geotechnics, vol. 18, no. 3, pp. 189–200, 1996. View at: Publisher Site  Google Scholar
 C. I. Teh, K. S. Wong, A. T. C. Goh, and S. Jaritngam, “Prediction of pile capacity using neural networks,” Journal of Computing in Civil Engineering, vol. 11, no. 2, pp. 129–138, 1997. View at: Publisher Site  Google Scholar
 S. Arlot and A. Celisse, “A survey of crossvalidation procedures for model selection,” Statistics Surveys, vol. 4, pp. 40–79, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 T. Fushiki, “Estimation of prediction error by using $K$fold crossvalidation,” Statistics and Computing, vol. 21, no. 2, pp. 137–146, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 R. Hornung, C. Bernau, C. Truntzer, T. Stadler, and A. Boulesteix, “Full versus incomplete crossvalidation: measuring the impact of imperfect separation between training and test sets in prediction error estimation,” Tech. Rep. 159, Department of Statistics, University of Munich, 2014. View at: Google Scholar
 CEN/TC 104, Concrete—Part 1: Specification, performance, production and conformity, EN 2061, 2000.
 CEN/TC 104, “Structural steel part 1: general technical delivery conditions,” Tech. Rep. EN 100251, 2004. View at: Google Scholar
 British Standard Institution, “Maritime works—part 2: Code of practice for the design of quay walls, jetties and dolphins,” BS 63492, 2010. View at: Google Scholar
 K. Levenberg, “A method for the solution of certain nonlinear problems in least squares,” Quarterly of Applied Mathematics, vol. 2, pp. 164–168, 1944. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 D. W. Marquardt, “An algorithm for leastsquares estimation of nonlinear parameters,” SIAM Journal on Applied Mathematics, vol. 11, no. 2, pp. 431–441, 1963. View at: Google Scholar  MathSciNet
Copyright
Copyright © 2014 Yongjei Lee 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.