An Optimization Algorithm for Exponential Curve Model of Single Pile Bearing Capacity

The ultimate bearing capacity of single pile is very important to engineering safety, so correctly predicting its value becomes an important part of engineering safety. Based on the traditional exponential curve model, a parameter optimization algorithm of the exponential curve model of single pile bearing capacity, which combines the golden section method and the linear least square method, is proposed. In order to verify the reliability of the proposed optimization algorithm, the measured data of the building engineering in the literature were optimized and calculated. Through comparison, it is found that the optimization algorithm is closer to the measured value than the traditional exponential curve model algorithm, which can better guide the engineering practice, and verify the effectiveness and superiority of the proposed parameter optimization method.


Introduction
The bearing capacity of piles can be determined based on conventional semi-empirical equations. In this regard, the Meyerhof equation is still well-respected. In addition, there are some simple correlations between the bearing capacity of piles and in-situ tests (like cone penetration test, CPT or standard penetration test, SPT); however, some studies suggest that the aforementioned correlations overestimate the bearing capacity. Furthermore, some other studies recommend the use of dynamic equations, which are based on the pile and hammer properties, to estimate the bearing capacity of piles; however, as stated by Milad et al. (Milad et al. 2015), there are many assumption and simplification in estimating the bearing capacity using the aforementioned equations. (Ehsan Momeni et al. 2020; HAN Jiwei et al. 2020).
Reasonable and correct prediction of vertical ultimate bearing capacity of single pile is the most important safety premise in pile foundation engineering. The existing models are widely used because of the advantages of exponential curve model. In the use of exponential curve model, the reasonable solution of its unknown parameters is the key problem for prediction. The existing methods mainly use the traditional intelligent optimization algorithm or gradient optimization algorithm to complete this task (Peng et al. 2010). In engineering, we find that we can explore more optimized algorithm, so we propose an optimization algorithm combining golden section method and linear least square algorithm. (Julong et al. 2012;Xing et al. 2016).

Problem Description
In order to predict the ultimate bearing capacity, the exponential curve model is used in this paper (Deling 2003;Nehdi 2015).
In it, S represents stable settlement (mm), Q is the test load corresponding to s ( kN), Unknown parameters in the model are A,B and theoretical ultimate bearing capacity Q u (kN) at s ? ?
It is assumed that the stable settlement of n group S i (i = 1,2,… n) and its corresponding test load Q i (i = 1,2,…n) are obtained by test pile test, then the unknown parameter can be obtained by optimizing that nonlinear least squares cost function (Jiang et al. 2016).
3 Exponential Curve Fitting Method

Separable Least Squares Problem
Genetic algorithm may also be used Genetic algorithms can also be used. The nonlinear least squares cost function defined by (2) can be solved by standard nonlinear least squares optimization method, such as Gauss-Newton method or Marquardt method. Genetic algorithm, differential evolution and other intelligent optimization methods can also be used to solve this problem. However (Yan et al. 2015;Zhang et al. 2015), the above methods do not take into account the characteristics of the unknown parameters in the model (1). For the convenience of analysis, order: Then the formula (1) and (2) are converted into: And then it's easy to know, If B is known quantity, The nonlinear least squares problem defined by (5) will be simplified to a linear least squares problem. In addition, if let: Another two unknown parameters in the formula (4): Through the analysis above, we can see that the nonlinear least squares problem defined by (5) is variable separable. For this kind of nonlinear least squares problem, the nonlinear parameters and linear parameters can be divided into two groups estimated separately. The nonlinear parameter B can be estimated by nonlinear optimization techniques. In this paper, golden section method as a commonly used method in univariate optimization problems is used. Q u and A1 are linear parameters, them can be obtained through formula (7).

Golden Section Method
The main idea of the golden section method is to find the optimal solution by narrowing the area of the optimal solution step by step. In the process of optimization, only the value of the finite cost function needs to be calculated. To calculate the objective function value required by the golden section method (Arnau et al. 2012), let B take a particular value, First, the linear parameters Q u and A1 corresponding to B is calculated by formula (7), the objective function J corresponding to B, Q u and A1 is calculated by the formula (5).
For one dimension optimization problem: min f ðxÞ; x 2 ½a; b Order constant h = 0.618, the first step in the Golden Section, let x 1 ¼ a þ ð1 À hÞ Â ðb À aÞ, The second step in the Golden Section is to compare the size of f ðx 1 Þ and f ðx 2 Þ. If f ðx 1 Þ f ðx 2 Þ, the interval of the optimal solution is reduced to [a, x 2 ], this situation order b ¼ x 2 , the interval of the optimal solution can still be recorded as [a, b], at the same time, let x 2 ¼ x 1 , f ðx 2 Þ ¼ f ðx 1 Þ, x 1 ¼ a þ ð1 À hÞ Â ðb À aÞ, calculate the value of f ðx 1 Þ; If f ðx 1 Þ ¼ f ðx 2 Þ, the interval of the optimal solution is reduced to ½x 1 ; b, at this point, a ¼ x 1 , the interval of the optimal solution can still be recorded as ½a; b, at the same time, let The second step is repeated several times until the length of the interval ½a; b is less than a predetermined constant e, so that the optimal solution x b ¼ ða þ bÞ=2 can be obtained.
After each iteration, the golden section reduces the length of the search interval to h times that of the previous one (Xue et al. 2015;Sun et al. 2015).

Algorithm Flow
In summary, the golden section-least square procedure of estimating the unknown parameters a\b and qu in the exponential curve model defined by the formula (1) is as follows: (1) Sets the constant e required to end the golden section search, order constant h ¼ 0:618; Sets the range of values of the parameter B ½a; b.
(4) If b À a ! e, return to the step (3); Otherwise, let the optimal solution B Ã ¼ ða þ bÞ=2, the linear parameter and corresponding to B* are calculated by formula(3), (4), the objective function J* corresponding to, and can be calculated by formula (5). Let A * = A* 1/Q* u, to find out the optimal solution B * , A * , Q* u and the optimal objective function J*.

Example of Practical Calculation
To verify the performance of the above algorithm, Taking into account the static load test data of Wuhan Power Market Building provided by reference (Deling 2003), Set the constant e ¼ 1 Â 10 À4 required by the golden section method, Set the range of parameter B½0; 20 (Fig. 1). According to the static load test data of Wuhan power market building provided by literature (Deling 2003), the predicted value is calculated by golden section least square method according to the algorithm flow of 3.3, and compared with the prediction data in literature (Deling 2003) and Literature. The detailed results are shown in Table 1.
The model fitting accuracy in Table 1 is defined as: where Q m i corresponds to Q i , The load value calculated by the model is KN, Q is the average value of measured load Q i (i = 1,2,… n) ( Fig. 2 and Table 2). Fig. 1 The picture of single pile

Conclusion
According to the traditional exponential model curve fitting method, an exponential model curve fitting optimization algorithm is proposed, and its theoretical prediction of the ultimate load of a single pile is deduced. The calculation results for typical problems show that this method has the advantages of small amount of calculation and high model accuracy. From the data, the fitting precision of curve fitting optimization algorithm of exponential model is higher than that of traditional method. Through comparison, it can be seen that the exponential model curve fitting optimization algorithm is more safe and reliable to predict the limit load value, and can meet the need of high precision in prediction or design.