Analytical Model and Back-Analysis for Pile-Soil System Behavior under Axial Loading

.e interaction mechanism between piles and soils is very complicated. .e load transfer function is generally nonlinear and is affected by factors such as pile side roughness, soil characteristics, section depth, and displacement..erefore, it is difficult to solve the pile-soil system based on load transfer function. .is paper presents a new method to study the soil-pile interaction problem with respect to axial loads. First, the shapes of the axial force-displacement curves at different depths and the displacement distribution curves along pile axis at different pile-top displacements were analyzed. A simple exponential function was taken as relationship model to express the relationship curves between two distribution functions of axial force and displacement along pile shaft obtained by using the geometric drawing method. Second, a new analytical model of the pile-soil system was established based on the basic differential equations for pile-soil load transfer theory and the relationship model and was used to derive the mathematical expressions on the distribution functions of the axial force, the lateral friction, and the displacement along pile shaft and the load transfer function of pile-side. We wrote the MATLAB program for the analytical model to analyze the influence laws of the parameters u andm on the pile-soil system characteristics. .ird, the back-analysis method and steps of the pile-soil system characteristics were proposed according to the analytical model. .e back-analysis results were in good agreement with the experimental results for the examples. .e analysis model provides an effective way for the accurate design of piles under axial loading.


Introduction
e pile foundation has the characteristics of high bearing capacity, good stability, small settlement, simple construction, and consumption of low material consumption. It has been widely used in construction engineering, marine engineering, industrial engineering, and other fields. In the study of pile-soil system under vertical load, there are actually two major problems: (a) to estimate the load-displacement curve of pile top by pile-soil parameters and (b) to investigate the system characteristics of pile-soil from piletop load-displacement curve and soil parameters. e former is a forward-analysis problem, and traditional theories such as Poulos elastic theory method [1], Seed and Reese load transfer method [2], and shear displacement method [3,4] are typical methods to solve this problem. e latter is a back-analysis problem and the researchers are fewer. e load transfer method was first proposed by Seed and Reese in 1955 [2]. Since the method is clear in concept and easy to understand, it has been immediately responded by scholars and engineers and has been rapidly developed [5][6][7][8][9][10][11][12][13][14]. e basic concept of the method is to assume that the pile-side and the pile-bottom are connected with the soil by springs. e constitutive relationships of the springs represent the load transfer relationship between piles and soils, and it is called the piles-soil load transfer functions. An element body is intercepted at a certain section along the pile axis (see Figure 1). According to the static balance and deformation compatibility of the element body, the corresponding differential equations can be established. Solving the differential equations depends on the pile-soil load transfer functions. erefore, many scholars have conducted extensive research on transfer functions and proposed various pile-side load function forms. Seed and Reese proposed hyperbolic model as the transfer function [2]. Kezdi proposed an exponential function relating skin friction to strain for the case of piles embedded in sand based on the model tests carried out in the laboratory or in the field [6]. Vijayvergiya [7] used the parabolic model as pile-soil load transfer function to analyze the axial behavior of piles. Armaleh and Desai [8] used the Ramberg-Osgood model to describe the pile side-soil transfer curve and the model provides a general method and includes other special functions, such as the hyperbola. According to the shear displacement theory that Randolph and Wroth proposed, the pile-soil load transfer curve in uniform soil is a linear function [4]. Sato [13] used a linear elastic perfect plastic model to calculate the bearing capacity of pile in multilayered soil. Currently, bilinear model [9] and trilinear model [9,11] were used to describe the pile side-soil load transfer curve and soften [9,11] and harden [9] phenomenon. e methods of solving the differential equation vary with the transfer function forms. Generally, for linear or exponential transfer functions, the results of elementary function expressions [9,11,13] or series form expressions [6] can be obtained by the analytic method solving the differential equations. For complex nonlinear pile-soil load transfer function, it is very difficult to solve differential equations using analytic method, which is generally solved by numerical methods, such as displacement coordination method [5,12,14,15], finite element method [8], and finite difference method [16]. e back-analysis method of pile and soil characteristics is often used to determine the geotechnical parameters. Danziger et al. [17] carried out plenty of back analyses for field data of offshore piles driven in calcareous sands using an improved soil model developed by Simons and the back analyses applied in a one-dimensional finite element to calculate the relative displacements between pile and soil at various points along the pile. Xiao and Yang [18] used Chin's plot method to back-analyze a large number of driven pile test data to obtain the two key design parameters, skin friction coefficient, and bearing coefficient. McCarron [19] executed the back-analysis of pile tests using the empirical one-dimensional load transfer method and used the model to research the cyclic capacity of the pile under the same conditions. e complex pile-soil lateral load transfer function is not only related to the nature of the soil, the nature of the pile, the lateral conditions of the pile, the pile forming method, the groundwater, the soil pressure on the pile-side, the distribution of the soil layer, and other factors, but also the function of two independent variables, the depth, and relative displacement between pile and soil at a cross section.
In pile foundation engineering, the designers pay close attention to the moments, shear forces, or axial force and deformation of piles in different depth piles so that it is properly reinforced. It must be obtained through special experiments or complex numerical analysis. None of the above forward-or back-analysis can accurately grasp the characteristics of the pile-soil system, which is quite different from the actual engineering properties. is paper focuses on the vertical forces of pile foundation such as bridge foundations or foundations of bracing systems.
Hence, this paper constructs a new analytical model of the pile-soil system based on the pile-soil load transfer theory and the back-analysis methods. It is a fact that the pile foundation is a concealed project, and most of the available test data which can be directly obtained are load-displacement curve and stratum survey data. It is beneficial to the pile foundation design and will be welcomed by engineering design technicians if the load-displacement curve of the pile top can be used to find the characteristics of the pile-soil system, such as the axial force distribution curve, the side resistance distribution curves, and the load transfer curves.

Differential Equation of Pile-Soil
Load Transfer e load transfer method was proposed by Seed and Reese [2]. e pile-soil system can be simplified as shown in  . e pile-side springs indicate the relationship between the pile-side resistance τ (z) and the shear displacement s (z), which are the pile-side transfer function (τ − s); the pile-bottom spring indicates the relationship between the pile-bottom resistance Q b and the pile-bottom displacement s b , which is the pile-bottom load transfer function (Q b − s b ). Figure 1(b) is a microelement body of pile and regardless of the weight of the pile. According to the static balance condition of the element body, we can obtain where z is the depth of the element body, Q (z) is the force distribution function along the pile axis, τ (z) is the skin friction of the pile-soil contact surface, and C is the perimeter of the pile section. By calculating the amount of elastic elongation produced by the element body, it is obtained that where E p is the elastic modulus of the pile, s (z) is the displacement distribution function along the pile axis, and A is the cross-sectional area of the pile. By taking the derivative of equation (2) and substituting equation (1) into it, the basic differential equation of the load transfer method can be obtained: According to equation (3), the displacement distribution function can be solved only if the pile transfer function τ (s, z) is correctly determined. Since the transfer function is very complicated, it is difficult to express it with a mathematical function. In this paper, through the relationship between the distribution function of axial force and the distribution function of displacement along the pile axis, the differential equations are further solved. Finally, a new mathematical model of the pile-soil system is obtained, which provides a new way for the accurate design of the pile.

The Relationship between Two Distribution Functions of Force and Displacement along Pile Axis
Xu et al. [20] used the geometric drawing method to analyze the geometric relationship curves, Q (z) − s (z), between the axial force distribution function along pile shaft, Q (z), and the displacement distribution function along pile shaft, s (z), for the vertical pile buried into soft soil and subjected to axial load. e test has confirmed that the displacement distribution function, s (z), decreases gradually along pile shaft [21]. e three curves in the s-z coordinate system of Figure 2 expressed the law and their corresponding pile-top displacements, s 0 , are s 01 , s 02 , s 03 , respectively, and the pilebottom displacements, s b , are s b1 , s b2 , s b3 , respectively. e test found that the pile-top load-displacement curves, Q 0 − s 0 , are similar to the models of the hyperbolic or exponential function [23,24]. Xu et al. synthesized the characteristics of two models and proposed a universally applicable power function model [22]. e relationship between the axial force, Q z , of the pile section at depth z and the displacement, s z , is similar to the curve of Q 0 − s 0 [6,24]. e initial stiffness and ultimate axial force of the Q z − s z curves decrease gradually with depth increasing. e five curves in the s − Q coordinate system of Figure 2 reflect the rules and their corresponding depths are 0, z 1 , z 2 , z 3 , z b , where 0 is the pile-top, z b is the pile-bottom, and z 1 > z 2 > z 3 are the pile sections of different depth. e Q z − s z curve at z � 0 is actually the pile-top load-displacement curve, Q 0 − s 0 . e Q z − s z curve at z � z b is actually the load transfer function Q b − s b of pile-bottom. e model of Q z − s z curve was proposed by Xu et al. [22] as follows: where K zi is the initial stiffness of the axial force-displacement curve at the depth z, Q zm is the asymptotic value of the ultimate axial force at the depth z, and n z is the exponent. Afterwards, we draw vertical lines upwards from any point of displacement distribution curve s(z) and intersect with the corresponding Q z − s z curves at depth z corresponding to the same depth. Connect these intersections to get the rough shape of the corresponding Q(z) − s(z) curve corresponding to the same pile-top displacement. e curves a, b, and c in Figure 1 are Q(z) − s(z) curves when the displacement of pile top is s 01 , s 02 , and s 03 . It can be seen that the functions Q(z) and s(z) are related to the displacement of pile-top, s 0 .  [20].

Mathematical Problems in Engineering
Xu et al. gave the mathematical expression of Q(z) − s(z) curve [20]. When the pile-top displacement is s 0 , where η and b are parameters related to the displacements of the pile-top and the pile-bottom, η > 0, 0 ≤ b ≤ s b ; u is the exponent, u > 0; and ϕ (z) is the influence function controlling the curve shape related to pile-top displacement s 0 , generally not less than 1 and monotone increasing. According to the boundary conditions of the pile-top and pile-bottom, it can be obtained from equation (5) that Solve equation (6) to obtain the parameters η and b.

Mathematical Model of Pile-Soil System and Its Analysis
Substituting equation (5) into equation (2), According to the displacement boundary condition of pile top, we obtain through solving the upper differential equation: where Φ is the original function of the influence function. Let Obviously, G (z) ≥ 0. From equations (9) and (10), e first and second derivatives of equation (11) can be obtained: where ϕ′(z) is the first derivative of the influence function ϕ (z). By substituting equations (12) and (9) into equations (2) and (3), respectively, the axial force distribution function, Q, and pile skin friction distribution function, τ, can be obtained as follows: Xu et al. suggested a specific form of the influence function [20]: where L is the pile length, r is the influence factor which is related to the displacement of the pile top, and m is the exponent. According to the condition that the influence function is an increasing function, it is easy to obtain a parameter range of (m According to equation (16), the original function of the influence function on z is Substituting the influence function (16) into equation (7), the specific expression of parameters b and η are obtained as follows: Substituting the original function (17) into equation (10), the specific expression of the G function is Substituting equations (16) and (20) into equations (11), (14), and (15), the displacement distribution function and the axial force distribution function along the pile axis and the frictional resistance distribution function are

Mathematical Problems in Engineering
From equations (16) to (23), it can be seen that Φ, ϕ, parameters b, η, and functions G, s, Q, and τ are all related to pile-top displacement s 0 .
Because the pile-top is the position with the largest displacement along pile axis and there is always the shearslip state between pile and soil, the skin friction at pile-top can be approximated as a constant, τ 0 . If the pile-top is at the same level as the ground and not loaded on ground, the lateral pressure of the pile-top can be approximated as 0. It can be considered that the lateral friction at pile-top is approximately equal to the cohesion between the pile and soil. erefore, when pile-top displacement, s 0 , is sufficiently large, the influence factor can be obtained from equation (23) according to the boundary conditions τ(z, s 0 )| z�0 � τ 0 , as shown in the following equation: e above analysis shows that u and m are pile-soil characteristic parameters related to pile-top displacement. To find the distribution functions of the displacement, axial force, and skin friction, the parameters b, η, and r must be obtained at first. In fact, by solving the equations of equations (18), (19), and (24), the values of b, η, and r can be obtained. Due to the complexity of above equations, the iterative algorithm can be used to obtain approximate values.
According to the force boundary conditions of the piletop, it is obtained from equation (5) that Substituting the above equation into equation (24) and simplifying, we can obtain that us, the following iterative equations can be established to obtain approximations of r and b, according to equations (18) and (26): where j is the number of iterations, r j represents the value of r at the jth iteration, b j represents the value of b at the jth iteration, and the initial value of the iteration b 0 � 0. e iteration continues until |b j − b j− 1 | ≤ δ; δ is a small positive real number, preferably 10 − 8 . At depth z k , the pile-side load transfer function τ z − s z can be expressed by the parametric equation composed of equations (21) and (23) taking the displacement of the piletop as the parameter: When the pile-top displacement is small, or the pile is exceptionally long or flexible, where the displacement zero point is yet to be transmitted to the pile-bottom, it is easy to (18) and (19). e depth of the displacement zero point (z 0 ) along the pile axis meets that s(z, s 0 )| z�z 0 � 0. From equation (21), we can obtain the depth of the displacement zero point as follows: Because z 0 > 0, then 1 − u > 0. Combined with the condition of equation (5), 0 < u < 1.
When the displacement is transmitted to the pile-bottom and the displacement s b is generated, the axial force of the pile-bottom is determined by the load transfer function of the pile-bottom, can be assumed to be the same model as the pile-top load-displacement curve model. According to equation (4), we obtain where K 0i and K bi are the initial stiffness of the load-displacement curve at the pile-top and the pile-bottom, respectively; Q 0m and Q bm are the asymptotic values of the ultimate load at the pile-top and the pile-bottom; and n 0 and n b are exponents. ese parameters are generally obtained by a field test data with fitting method. If the pile-bottom displacement is s b , the whole pile has a rigid body displacement s b . e pile-bottom displacement s b can be obtained by an iterative method, and the iteration equation is where s (n) b is the pile-bottom displacement at the nth iteration. Generally, the iteration initial value where δ is a small positive real number, preferably 10 − 8 .

Forward Analysis Program (FAP).
e question of analysis is knowing the geometric, physical, and mechanical parameters of pile and soil, the pile-top load-displacement curve, the pile-bottom load transfer function, and the values u and m, seeking the distribution curves of the displacement, the axial force, and the skin friction along pile axis.

Influence of the Parameter u on Pile-Soil
Characteristics. For the super-long pile with equal section, let the pile diameter d � 0.35 m, the pile-top displacement s 0 � 5 mm, the pile-top load Q 0 � 521.94 kN, the fixed parameter m � 5, the parameter u vary from 0 to 1. rough calculations, the influence of the parameter u on the characteristics of the pile soil is shown in Figures 4 and 5. It can be seen from the figures that the influence of the parameter u on the pile-soil characteristics is as follows.
(a) e depth of displacement zero point (or axial force zero point) increases as u increases. (b) When u ≤ 0.5, the pile end friction has a concentrated phenomenon. e pile end resistance is relatively large, and the upper side resistance is relatively small, representing the soil above soft and below hard. When u is small, it means that the pilebottom is in the base layer, and pile-bottom axial force is mainly borne by pile end resistance. (c) When u > 0.5, the lateral frictional resistance is distributed in an arch shape. When u � 0.6∼0.7, the distribution of lateral resistance is similar to that of equal section piles in uniform soil. When u � 0.7∼0.9, it forms an asymmetrical normal distribution and shows that the soil has certain characteristics of upper hard and lower soft.

Combined Control of Pile-Soil Characteristics with
Parameters u and m. e above analysis shows that u and m are closely related to the physical and mechanical properties of pile and soil. e parameters u and m always jointly control the pile-soil characteristics. When m is small (e.g., less than 0.5), u and m are equally important in the degree of control; when m is large (e.g., greater than 1), the pile-soil characteristics are mainly controlled by u, and that m only plays a minor role in the adjustment of the distribution pattern.

Back-Analysis of Pile-Soil System
e above analysis indicated that parameters u and m determine the distribution patterns of pile axis displacement, axial force, and pile side frictional resistance. If the pile-soil characteristic parameters (u, m) can be obtained using the 6 Mathematical Problems in Engineering  (21) and (22), and drawing Calculate s (z), Q (z) from equations (21) and (22), and drawing Calculate τ(z) from equation (23), and drawing Calculate τ (z) -s (z) from equation (28), and drawing End Mathematical Problems in Engineering back-analysis method through known test data, such as loaddisplacement curve of pile-top and pile-bottom and survey data of soil, we can solve the distribution curve of axial displacement, axial force, and pile side frictional resistance according to the analytical model established above, which is very beneficial to the design of pile.

Diversity of Back-Analysis Results.
In general, the forward-analysis results are unique; that is, for a particular pilesoil system we can obtained only the pile-top Q 0 − s 0 curve. Otherwise, there can be several pile-soil systems corresponding to a pile-top Q 0 − s 0 curve. In other words, the results of the back-analysis are not unique. e correspondence between the parameters (u, m) and the pile-soil system is complex and difficult to express with mathematical formulas, so it is impossible to find a unique result through mathematical relationships. In order to obtain the unique pile-soil system, other known information in the pile-soil system may be added, such as the distribution curve of frictional resistance, the distribution curve of axial force, the distribution curve of displacement, and the pilesoil load transfer curve of a certain soil layer. e optimal pile-soil system can be found by means of optimization methods through the known information.

Optimization Objective Function.
e best parameter set (u, m) can be found through the optimization method. Assume that the distribution curve of the displacement, the axial force, or the shear force along the pile axis has been measured by test under the pile-top displacement s 0 , and the test value (y j ) on the j th section depth (z j ) is {(z j , y j ): (j � 1, 2, . . ., k)}, where k is pile section all number. Set a parameter group (u, m), to calculate the theoretical distribution curve of the displacement, the axial force, or the shear force along the pile axis under the pile-top displacement s 0 through mathematical model of pile-soil system, and the theory value (y j ) on the j th section depth (z j ) is {(z j y j ): (j � 1, 2, . . ., k)}. In order to verify whether the setting parameter group (u, m) is optimal, the theoretical distribution curve must be the most consistent with the experimental distribution curve. e coefficient of determination is used as the optimization objective function [25] as follows: where R 2 is the coefficient of determination, representing the degree of conformity between the test curve and the theoretical curve, and the closer R 2 to 1, the better the degree of conformity between the two curves; SST stands for the total sum of squared deviations; SSE represents the sum of

Back-Analysis and Optimization
Steps. In the backanalysis, first fix an m value, such as m � 2, then let u go from 0 to 1 by increasing step size Δu, and calculate the objective function R 2 for each increment step. When the value of R 2 is the closest to 1, the value of u is the optimal value. en fix u as the optimal value, increase m from 0 to 10, and find the optimal value of m when R 2 is the closest to 1. e precision of optimal parameter set (u, m) is related to increasing step size. Generally, the small step length corresponds to high precision but slow calculation speed. e calculation steps are as follows.
(a) Selection of optimizable object: one of the distribution curves of the frictional resistance, the axial force and the displacement distribution curve, or the pile-soil load transfer curve can be selected as the optimization object. e back-analysis results are compared with the measured results. (b) Optimization of the parameter u: fixing m value, u is incremented from 0 to 1, and the incremental step is Δu. e theoretical value of optimization object can be calculated using FAP. By comparing the theoretical value with the experimental value, the variation rule of the objective function R 2 with the parameter u can obtain, and the optimal parameter u 1 is found. (c) Optimization of the parameter m: fixing u � u 1 , m is incremented from 0 to 10, and the incremental step is Δm. e theoretical value of optimization object can be calculated using FAP. By comparing the theoretical value with the experimental value, the variation rule of the objective function R 2 with the parameter m can obtain, and the optimal parameter m 1 is found. (d) Optimization of the parameter group (u, m) again: fix m � m 1 and repeat from step (b) to get the better parameter u 2 . en fix u � u 2 and repeat from step (c) to get the better parameter m 2 . is process can be repeated until the best parameter group (u, m) is obtained.

Example of Back-Analysis for Tension Pile
Example 1. Knowing that a cross-section pile diameter d � 500 mm and length L � 10 m, the mechanical parameters of the pile are elastic mode E p � 2.7e7 kPa and Poisson's ratio μ p � 0.18. e basic mechanical parameters of the cohesive soil are elastic modulus Es � 30000 kPa, Poisson's ratio μ s � 0.38, cohesion c � 10 kPa, internal friction angle φ � 20°, and pile-soil friction coefficient f � 0.4. e uplift load-displacement curve of pile-top is obtained by finite element calculation, as shown in Figure 8 [26].
e pile-top load-displacement curve is fitted by using equation (30), and the values of the parameters in the model are obtained that Q 0m � 348.7 kN, n 0 � 2, K 0i � 96.5 kN/mm, and the model equation According to the above back-analysis steps, the parameter group (u, m) is optimized as follows.
(a) When the pile-top load is 243.4 kN, the axial force distribution is selected as the optimization object. e theory results of the back-analysis are compared with the axial force obtained by the finite element calculation along pile axis (see "o" in Figure 9). (b) Preliminary optimization of u. Make m � 2, u increase from 0 to 1, the incremental step is 0.1, and the change of the objective function R 2 with u is shown in Figure 10(a). It can be seen that u � 0.7 is the optimal value. (c) Fine step optimization of u. Make m � 2, u increase from 0.6 to 0.8, the incremental step is 0.01. We can obtain optimal values that u � 0.67 and R 2 � 0.999382. (d) Preliminary optimization of m. Make u � 0.67, m increase from 0 to 10, the incremental step is 0.5, and the change of the objective function R 2 with m is shown in Figure 10(b). It can be seen that m � 2 is the optimal value. (e) Fine step optimization of m. Make u � 0.68, m increase from 1.5 to 2.5, and the incremental step is 0.01. We can obtain optimal values that m � 2.02 and R 2 � 0.999383.
Finally, the optimal parameter group of the pile-soil system is (u, m) � (0.67, 2.02). Figure 9 shows comparison between two results of back-analysis and FEM.

Mathematical Problems in Engineering
As the distribution pattern of pile lateral friction changes dynamically during the uplifting process, the parameter group (u, m) also changes with the pile-top uplift displacement. According to the influence of the parameter group (u, m) on the pile-soil characteristics introduced in Section 5.2, it can be concluded for tension piles that u is a decreasing function and m is an increasing function. Based on a lot of trial compute and analysis, a coincidental result showed that the variation law of the parameter group (u, m) with pile-top displacement is related to the model of pile-top load-displacement as follows:  .4974) can be obtained. By substituting above parameter groups into the FAP, the axial force distribution curves and the frictional resistance distribution curves can be obtained (see Figures 11 and 12 respectively).
According to the parametric equation (28), the curves of pile-soil load transfer can be obtained at the different depth (see Figure 13). It can be seen that the maximum friction increases with the increase of the section depth except the pile-end and softening appears in the shallow layer. ese laws are in good agreement with the experimental results of Zhao et al. [27]. It shows that the results of the back-analysis are reliable.

Back-Analysis Example for Compression Pile
Example 2. (see [28]). Two engineering test piles for the pile foundation of coal storage bunker were designed that the diameter is 1 m and depth is 25 to 27 m. e no. 2 test pile is selected as the back-analysis calculation. e test pile no. 2 uses cast in situ bored forming method, and its average diameter d � 1.1 m and the length L � 27 m. e top to down soil layers are filled soil, clayey soil, clay soil, clayey soil, coarse sand, clayey soil, coarse sand, clayey soil, coarse sand, clay soil, coarse sand, and hard clay. e pile-bottom passes through a coarse sand layer and is supported on a hard clay layer. e test load-displacement curves of pile-top and pilebottom are shown in Figure 14.
e measured value of ultimate frictional resistance τ m of pile side is shown in Figure 15. (b) Selection of the optimization object: take the pileside ultimate friction as the optimization object when s 0 � 50 mm. Due to the heterogeneous soil layer, the ultimate friction distribution is very discrete and it may cause the optimization distortion. erefore, the polymorphic fitting is first performed for the test ultimate friction τ m , the fitting curve is shown in Figure 15, and the expression is as follows (correlation coefficient 0.98783): e fitting curve is taken as optimization object to be compared with the back-analysis results. (c) Optimization: first, the initial values are estimated that Q bm � 1600 kN and τ 0 � 0. Second, with the method of Example 1 for optimization, u � 0.41 and m � 0.13 are obtained. ird, Q bm is optimized to get Q bm � 1786 kN. Fourthly, τ 0 is optimized to get τ 0 � 0. e above optimization process can be repeated 2 to 3 times, and the final objective function value of R 2 � 0.994946. e maximum resistance of the pilebottom is 1710 kN by Yukang and Qiang-hua [28]. e dashed line in Figure 14 is the load-displacement curve of pile-bottom which was obtained by backanalysis. It can be seen that the error between the back-analysis result and the measured value is small. (d) Considering the dynamic variation law of the lateral resistance, which is along with the displacement of the pile top, is different between the compression pile and the uplift pile Considering that the dynamic change of lateral resistance of the compressive pile with the pile-top displacement is different from that of the tension pile [26]; therefore, the variation rule of parameters (u, m) with the pile-top displacement is also different from that of drawing pile.
According to the influence rule of parameters (u, m) on the pile-soil characteristics described in Section 5.2, it can be concluded that for the compression pile, u is a decreasing function, but the change magnitude is smaller than that of equation (35), and m basically stays the same. e functions of the parameters (u, m) are determined as follows after a lot of trial calculation:   where u c and m c are constants, which are related to the physical and mechanical properties of the pile and soil and are independent of the pile-top displacement. e optimizational parameters u � 0.41 and m � 0.13 are substituted into equations (38) and (39) to obtain u c � 0.36 and m c � 0.13. en we calculate the axial force distribution, the friction distribution, and the load transfer curves at the different pile-top displacement. Figure 16 is the friction distribution curves along pile axis. e pile-top displacement and pressure forces corresponding to every curve from the left to the right are (s 0 , Q 0 ) � (1.183, 1600), (2.876, 3200.3), (5.617, 4800.4), (11.404, 6400.3), and (50, 7923.6) (mm, kPa), respectively. Symbol "○" in Figure 16 is the measured ultimate friction. Figure 17 is the load transfer curves (τ z − s z ) at different depth. It can be seen that the back-analysis results reflect the basic law of load transfer correctly.

Conclusion
e piles acting as axial force are common in foundation engineering, such as high-rise building foundation, bridge foundation, and antifloating pile. e pile-soil load transfer theory is widely used in pile foundation analysis. At present, there are two main algorithms of the load transfer analysis based on the theory of the analytical method and deformation coordination method [1,4,6]. ese two algorithms need to establish the load transfer function (τ z − s z ) of the pile-side and the load-displacement function (Q b − s b ) of the pile-bottom. Due to the complexity of pile-soil interaction, load transfer functions are usually nonlinear and vary with depth. erefore, it is difficult to solve the basic differential equations based on load transfer theory. erefore, the analytical method is generally difficult to implement for solving basic differential equations based on load transfer theory. In order to overcome this problem, a new algorithm is proposed in this paper, that is, the basic differential equations of load transfer are solved by constructing a relationship model between the distribution functions of axial force and displacement along the pile axis, and the analytical results of the force and the displacement on pile-soil system are obtained.
Establishing the relationship model is the key in the new algorithm. After analyzing the shapes of the axial loaddisplacement curves at different depth and the displacement distribution curves along the pile axis under different piletop displacements, the relationship curves between the distribution functions of axial force and displacement along pile shaft are obtained by using geometric mapping method. e relationship curves are simulated by a simple exponential function, that is, the relationship model between two distribution functions of axial force and displacement. e relational model was applied for the first time to solve the differential equations of load transfer, and then the analytical model of the force and deformation analysis of the pile-soil system was obtained. rough the analytical model, the axis displacement distribution function, axial force distribution function, pile side friction distribution function, and pile side load transfer function can be obtained. e MATLAB program of the analytical model is written.
e program is used to analyze the influence of parameters u and m on the characteristics of pile-soil system. e analysis shows that the parameters u and m are closely related to the physical and mechanical properties of the pile and soil and control the shape of lateral friction distribution curves. When m is small (e.g., less than 0.5), u and m are equally important in controlling the shapes of the  distribution curves; when m is large (e.g., greater than 1), the shapes of the distribution curves are mainly controlled by u, and m only plays an auxiliary effect. If the parameter m is unchanged (m � 5), when u ≤ 0.5, the soil layer shows the characteristics of upper-soft and lower-hard, and the smaller u is, the softer the upper is and the harder the lower is; when u � 0.6∼0.7, the soil layer showed uniform characteristics; when u � 0.7∼0.9, the soil layer exhibits upper-hard and lower-soft characteristics. If the parameter u is fixed (u � 0.6), when m is small, the soil layer exhibits upper-hard and lower-soft characteristics; with the gradual increase of m, the influence degree of m on the upper-hard and lower-soft characteristics gradually decreases until it is mainly controlled by parameter u.
According to the analytical model, a back-analysis method for the characteristics of pile and soil system is proposed. e distribution curves of the axial force and the displacement along pile axis and the pile-soil load transfer curve can be obtained from the pile-top load-displacement curve and the stratum survey data using the back-analysis method. e examples illustrate that the back-analysis results are consistent with the experimental results. e analysis finds that the parameters u and m are related to the pile-top load-displacement curve, and then the empirical formulas for the variation of the parameters u and m with the pile-top displacement under tension and compression are given.
When using the analytical model in this paper, it is important to determine the values of parameters u and m. In engineering design, the values of parameters u and m can be obtained by back-analysis method. In order to facilitate the use of this analytical model by engineers, the authors will further study the correlations between the parameters u and m and the mechanical parameters of the soil and pile.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest regarding the publication of this paper.