Static and Dynamic Reliability Analysis of Laterally Loaded Pile Using Probability Density Function Method

: Pile foundation is one of the common foundation forms in marine geotechnical engineering, especially in wind power engineering. Its operation safety is seriously a ﬀ ected by many uncertainties, such as the randomness of ground motion in intensity and frequency. The stochastic reliability analysis method can better characterize these uncertainties in the evaluation of the safety performance of pile foundation. The probability density functions (PDFs) of stochastic systems are important prerequisites for reliability analysis. However, for geotechnical problems, the coupling between parametric and excitation randomness and the nonlinear mechanical properties of rock and soil make it very di ﬃ cult to obtain the associated PDFs. Instead, the probability density evolution method (PDEM) is introduced and is used to investigate the static and dynamic reliability of laterally loaded piles as an example of a geotechnical problem. Compared with Monte Carlo stochastic simulations, PDEM-based computing is shown to be highly e ﬃ cient when applied to the seismic design of pile in geotechnical engineering, and its calculation e ﬃ ciency is 20 times of the former for the seismic dynamic reliability of pile foundation. This study provides a new reference for the e ﬃ cient design and safety evaluation of o ﬀ shore pile foundation engineering based on static and dynamic reliability of multiple random factors. The results from the stochastic dynamic system of a laterally loaded pile subjected to seismic excitation show that the probability information is transported strongly by the nonlinear state evolution. They also indicate that the probability information of the seismic dynamic response of a pile cannot be easily determined from the excitation distribution. The assumption applied in traditional reliability analysis, namely, that the response and excitation have the same similar distribution, should be reconsidered. Different from the classical traditional reliability analysis, especially the random analysis based on Monte Carlo stochastic simulation, the main feature of this study is to obtain the probability density function of nonlinear static and dynamic pile foundation system. Combining with the new high-dimensional probability space dimension reduction sampling technology such as TS sampling, the calculation amount of stochastic analysis, especially the nonlinear seismic dynamic reliability analysis, can be effectively reduced. study, results


Introduction
In foundation engineering involving large piles, pile foundation is one of the common structural forms of offshore wind power. It is generally the case that the axial force is the main control load and buckling is the main failure mode [1,2]. However, because of lateral action from seismic excitation, wind, waves, and various other water and soil pressures, a single pile or a pile group is usually also subjected to huge lateral loading; examples of this are offshore platforms, transmission-line towers, high-piled wharfs, and the piers of long bridges. Therefore, there have been many studies of the static and dynamic responses of piles. These range from closed-form exact solutions of simplified problems [3,4] to numerical solutions of complicated nonlinear problems [5,6]. Meanwhile, many researchers have gradually come to realize the importance of uncertainty and randomness when analyzing the static and dynamic responses of a laterally loaded pile system. In pioneering work, Zhang et al. [7] systematically analyzed and summarized four main sources of analytical uncertainty and randomness for laterally loaded piles and proceeded to study the static reliability of piles considering the variability of the soil modulus. Aimed at soil-structure interaction, Barakat et al. [8] proposed a reliability and optimum analysis and design method for laterally loaded piles.
Recently, thanks to rapid computational developments, various analyses have been conducted of the linear and nonlinear stochastic responses of laterally loaded piles using finite-element (FE) and finite-difference approaches [9][10][11][12][13]. These studies considered multiple correlated sources of randomness in the static reliability analysis of a laterally loaded pile using the first-order and second-order reliability methods, the response-surface method, and Monte Carlo simulation, among others. However, for Monte Carlo simulation, the anomaly of stochastic convergence means that extensive computation is required [14,15]. Obviously, due to the limitation of the basic principle of Monte Carlo stochastic simulation, the static and dynamic reliability analysis of piles based on Monte Carlo method will inevitably face a huge amount of calculation. Furthermore, most other approaches to calculating reliability depend on the Gaussian distribution, which inevitably reduces the accuracy because the correlation coefficient is changed in the nonlinear transformation from others to Normality of probability distributions [16]. In essence, the stochastic analysis of static and dynamic response of geotechnical engineering is far less than that of other distribution types, whether the physical quantities are assumed to be random variables or random processes with normal distribution. This is because probability density function of normal distribution is determined by its mean and variance. Compared with other probability information, especially the probability density function, it is much less difficult to obtain the mean and variance, and the probability density function is also the most essential description of the objective physical world. Meanwhile, all the aforementioned studies involving reliability analyses of laterally loaded piles were aimed only at static or quasi-static analysis of piles subjected to seismic loading. Less attention has been paid to the nonlinear time-history analysis of the dynamic reliability of laterally loaded piles.
When a pile is subjected to either a wind or wave load or seismic excitation, analyzing the response of that pile is a dynamic problem; if the loads being considered are in any way uncertainty or stochastic, it becomes a stochastic dynamic problem. A previous study used time-history analysis and first-passage theory to investigate the linear dynamic response and reliability of a large offshore structure considering the variability of wave loading and earthquake excitation [17]. Many studies have also used the pseudo-static method to assess the dynamic reliability of piles considering the uncertainty of earthquake ground motion [18]. Unfortunately, none of those attempts succeeded in capturing the coupling between randomness and nonlinearity. Furthermore, for stochastic dynamic analysis, Monte Carlo simulation seems to be the most universal method for a stochastic system of arbitrary dimension. However, the already huge amount of computation will grow explosively as the dimensions of the random system are increased [15,[19][20][21], data dimensionality reduction is a common challenge for all stochastic uncertainty problems, and it is also a problem that this study needs to face. That is, by introducing a new sampling method of high-dimensional probability space, the computational complexity of stochastic static and dynamic analysis can be reduced.
Given the discussion above, the probability density evolution method is introduced herein [22]. It is used to assess the nonlinear static and dynamic responses of laterally loaded piles from the perspective of the probability density functions (PDF), which also provides a foundation for reliability analysis.

Generalized Probability Density Evolution Equation
A general static or dynamic system of laterally loaded piles subjected to static and earthquake loading can be expressed as where Θ is a vector of random parameters (e.g., the uncertainty of the soil properties or the randomness of the earthquake ground motion) whose role is to describe all the possible sources of randomness in the aforementioned static or dynamic system, ξ is a generalized time parameter representing the evolution direction of the system state, and u and t are the spatial location vector of the system and its situated time parameter, respectively.
The components of Equation (1) can be expressed as where n is the dimension of the system and combines with the physical quantity U r which attracts us. The switching rate of the system is the derivative of U r with respect to ξ: Here, it is assumed that the stochastic static or dynamic system of the laterally loaded pile is probability conserving and that the stochastic factors of the system are given entirely by Θ. The joint probability density of (U(ξ), Θ) is denoted as p UΘ (u, θ, ξ). It should be emphasized that if no new stochastic factors are added or removed in the process of system evolution, the probability of the stochastic system is conservative [23], namely, where Ω ξ is the value domain of U at time ξ and Ω θ is the value range of Θ. ξ 0 is treated as the initial time of system evolution. Equation (4) can be rewritten as Because of the arbitrariness of Ω ξ × Ω Θ , the generalized probability density evolution equation (PDEE) can be obtained from Equations (3) and (4): ∂p UΘ (u, θ, ξ) ∂u r = 0 (6) in which all the variables are the same as those in Equations (2)-(4). Equation (6) portrays the transmission regularity of the stochastic in general physical system; obviously, even the static or dynamic system of laterally loaded piles. Therefore, the PDF p U (u, ξ) of U(ξ) can be calculated by The aforementioned process of solving for the PDF of a physical quality is known as the probability density evolution method (PDEM) [23,24]. Generally speaking, for a laterally loaded pile whose engineering design is investigated from the perspective of reliability, there is only one response (e.g., displacement or internal force), whereupon Equation (6) can be reduced to a one-dimensional (1D) partial differential equation, namely, To solve the generalized PDEE, its initial conditions must be established here. For the 1D case, the initial conditions can be expressed as follows if the stochastic factor is without respect to the initial conditions of the pile: where u 0 is the deterministic initial value of U(ξ) at the time instant ξ 0 and δ(·) is the Dirac delta function. In general, although the nonlinearity of the soil properties makes it hugely challenging to obtain a closed-form exact solution of Equation (6) or Equation (8), the PDEE can be solved numerically as follows: (1) Conduct the probability division of the vector space Ω Θ of the stochastic parameters and then select representative sample point θ q (q = 1, 2, · · · , N sel ) in each subdivision, where N sel is the total number of subdomain divisions of the probability space and p q is the corresponding assigned probability measure of the subdivision region. By doing so, the stochastic response analysis of the laterally loaded pile can be translated into a series of deterministic calculations. The specific division technique (e.g., tangent sphere (TS) or number theory) can be chosen according to different probability spaces which combine with special stochastic parameters. This aspect is discussed in Sections 3.1 and 3.2 for static and dynamic problems, respectively; (2) Introduce the deterministic response into Equation (8) as a generalized velocity . U r , whereupon the PDF can be obtained by solving the PDEE by the finite-difference method.

Numerical Cases
This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

Static Reliability
In this section, the static reliability of a laterally loaded pile is investigated here (Figure 1). To compare the effectiveness and accuracy of the proposed method and solving techniques with those of existing research and the Monte Carlo method, a pile problem that has an exact solution is selected. Specifically, the lateral horizontal displacement of a pile under a force is introduced herein, the solution of which is known as the Hetényi equation [25,26]. Here, the diameter and length of the pile are d = 0.4 m and L = 5 m, respectively, and subject it to a horizontal force F and a moment M. It is assumed that these external influences as random variables with a correlation coefficient of ρ(F, M) = 0.5. The material properties of the pile are described by an elastic modulus of E p = 200 GPa and a section moment of inertia of I p = 3 × 10 −4 m 4 . The important parameter of the horizontal subgrade modulus k h is also selected as a random variable because of its variability in a free laterally loaded pile. obtain a closed-form exact solution of Equation (6) or Equation (8), the PDEE can be solved numerically as follows: (1) Conduct the probability division of the vector space Ω of the stochastic parameters and then select representative sample point ( = 1,2, ⋯ , ) in each subdivision, where is the total number of subdomain divisions of the probability space and is the corresponding assigned probability measure of the subdivision region. By doing so, the stochastic response analysis of the laterally loaded pile can be translated into a series of deterministic calculations. The specific division technique (e.g., tangent sphere (TS) or number theory) can be chosen according to different probability spaces which combine with special stochastic parameters. This aspect is discussed in Sections 3.1 and 3.2 for static and dynamic problems, respectively; (2) Introduce the deterministic response into Equation (8) as a generalized velocity , whereupon the PDF can be obtained by solving the PDEE by the finite-difference method.

Numerical Cases
This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusions that can be drawn.

Static Reliability
In this section, the static reliability of a laterally loaded pile is investigated here (Figure 1). To compare the effectiveness and accuracy of the proposed method and solving techniques with those of existing research and the Monte Carlo method, a pile problem that has an exact solution is selected. Specifically, the lateral horizontal displacement of a pile under a force is introduced herein, the solution of which is known as the Hetényi equation [25,26]. Here, the diameter and length of the pile are = 0.4m and = 5m, respectively, and subject it to a horizontal force and a moment . It is assumed that these external influences as random variables with a correlation coefficient of ( , ) = 0.5. The material properties of the pile are described by an elastic modulus of = 200GPa and a section moment of inertia of = 3 × 10 m . The important parameter of the horizontal subgrade modulus is also selected as a random variable because of its variability in a free laterally loaded pile.  All the probability information regarding the random variables is listed in Table 1. The horizontal displacement of the pile subjected to force F and moment M is expressed by Equations (10) and (11): where λ is a parameter expressed as Table 1. Probabilistic information on the random parameters for static reliability.

Parameter Distribution Mean COV Correlation Coefficient
Log-normal 10,000 0.2 Independent Therefore, the performance index of the pile is selected as y z=0 = y(0) . According to the numerical solution process of the PDEM, the first step is to select representative sample points by dividing the probability space which combined with the random variables F, M, and k h . In general, there are many random-variable sampling methods in the field of stochastic analysis. Considering sampling accuracy and efficiency, the TS point-selection technique is arguably the best one to combine with PDEM for a probability space of two or three random variables. The essence of this sampling technique is sphere packing: the probability space is represented as accurately as possible by a certain number of sampling points. Existing research suggests that there are no more than 12 equal-radius spheres tangent to a sphere in three-dimensional (3D) space. Therefore, the basic idea of TS sample points selection is that it facilitates easy construction of a coordinate representation of the centers of the TSs in 3D space. Suppose that the centers of the middle-layer spheres are located on the z = 0 plane and that the sphere-center coordinates of the other layers in the z direction can be expressed as follows according to the symmetry: where N z is the layer number of the TS sample technology in the z th layer. We assume that θ j ∈ [ −η, η ]( j = 1, 2, 3) are standardized random variables and that η is the point-selection boundary of the space of standard random variables. For instance, a reasonable range of η is 3.2-4.0 for a standard normal distribution. Here, the value of N z can be taken as int η ∆z + 1, where int(·) represents the maximum integer that is no more than the value in the brackets and ∆z = 2 where R is the radius of the TS and l is the number of TS layers around the sphere in the center in the z = 0 plane.
If the sphere-center coordinates of the z k -layer spheres are x i,k , y i,k , z i,k , then Coordinates (x i,0 , y i,0 , z k ) are those of the z = 0 plane and can be determined by where r i and ϕ are the radius and polar angle, respectively, in the polar coordinate system. The joint PDF in the 3D probability space formed by the three standardized random variables θ j , j = 1, 2, 3 is almost spherically symmetric and radially attenuated. Based on this situation, only those sampling points inside spheres of radius η are selected. In other words, the coordinates x i,k , y i,k , z i,k should meet the following requirement: Here, the coordinates are renumbered as x q , y q , z q , q = 1, 2, · · · , N sel . Distribute these sampling points around the z axis with fixed angle ψ = π 12 . to make the sampling more homogeneous, therefore: Accordingly, the assigned probability of the representative sampling points can be calculated as where V q is the volume enclosed by the 12 tangent planes and θ qi,j = θ qi + r η θ i,j (i = 1, 2, 3). The representative sampling points selected above are standardized to obey a standard normal distribution on the domain [−η, η]. However, the target probability distributions of the random variables F, M, and k h are Gumbel and log-normal. Therefore, the sampling points should be transformed from standard normal space into the corresponding target probability space. Because doing so will change the correlation coefficient, it uses the Nataf transformation to calculate a modified correlation coefficient in light of the following approximation [16]: where R coe is the ratio of the correlation coefficients before and after transformation, ρ Y i Y j is the correlation coefficient of the target distribution, and ρ X i X j is the original correlation coefficient between F and M, namely, ρ X i X j = ρ(F, M) = 0.5. Beyond that, however, the correlation coefficient can also be calculated numerically as where µ X and σ X are the mean and standard deviation, respectively, of random variable X. Herein, random variables X i and X j correspond to F and M, respectively, with ϕ 2 y i , y j ρ Y i Y j being the conditional cumulative distribution function (CDF). According to Equations (18) and (19), the correlation coefficients of the target probability distributions are 0.5154 and 0.5155, respectively. The results demonstrate that the above two categories methods are equally effective. Thus far, it has selected N sel = 527 sampling points in the 3D probability space formed by the random variables F, M, and k h according to the TS sampling technique in this study. By comparing the mean and variance of different sample points calculated by PDEM and Monte Carlo stochastic method, it can be found that when 527 sample points are selected, the calculation error accuracy of the random analysis can be controlled within 3%. Therefore, 527 sample points are selected for stochastic reliability analysis. Next, a series of deterministic displacements of the laterally loaded pile can be obtained using Equations (10) and (11). Introducing the displacements into the PDEE (Equation (8)) as the generalized velocity . U r , the PDF ( Figure 2a) and CDF (Figure 2b) of the displacement of the pile top can be obtained by solving the PDEE. In Monte Carlo stochastic simulation, the reliability (CDF) is achieved by calculating the percentage of the number of samples in safe and effective state to the total number of samples by different thresholds of performance (i.e., displacement) index. The corresponding probability density function is obtained by the Kernel Density Estimation of the total sample. Different from Monte Carlo stochastic simulation, PDEM is obtained by the probability density function of performance index distribution, while the corresponding reliability (CDF) is obtained by integrating the probability density function. The probability density function obtained by PDEM is based on the objective physical differential equation, so it has higher calculation accuracy and less calculation. The accuracy of Monte Carlo stochastic simulation is greatly affected by the calculation samples, especially the Kernel Density Estimation method, which is also affected by the selection of window function. The comparison of their accuracy can be characterized and found in Figure 2. To verify and validate the effectiveness and practicability of the PDEM and TS sampling technique, it also can calculate the PDF and CDF via kernel density estimation and Monte Carlo simulation. Herein, to compare to the previous research work [26], the number of Monte Carlo stochastic simulation also is 20,000. It should be emphasized that the accuracy and efficiency of the static reliability herein come from combining the PDEM and TS sampling technique organically. Even if with this sampling (TS) method, the Monte Carlo stochastic simulation is also impossible so effective. is obtained by the Kernel Density Estimation of the total sample. Different from Monte Carlo stochastic simulation, PDEM is obtained by the probability density function of performance index distribution, while the corresponding reliability (CDF) is obtained by integrating the probability density function. The probability density function obtained by PDEM is based on the objective physical differential equation, so it has higher calculation accuracy and less calculation. The accuracy of Monte Carlo stochastic simulation is greatly affected by the calculation samples, especially the Kernel Density Estimation method, which is also affected by the selection of window function. The comparison of their accuracy can be characterized and found in Figure 2. To verify and validate the effectiveness and practicability of the PDEM and TS sampling technique, it also can calculate the PDF and CDF via kernel density estimation and Monte Carlo simulation. Herein, to compare to the previous research work [26], the number of Monte Carlo stochastic simulation also is 20,000. It should be emphasized that the accuracy and efficiency of the static reliability herein come from combining the PDEM and TS sampling technique organically. Even if with this sampling (TS) method, the Monte Carlo stochastic simulation is also impossible so effective.

Seismic Dynamic Reliability
In general, a laterally loaded pile subjected to seismic excitation is a strongly nonlinear dynamic system because of the nonlinear properties of concrete (i.e., the pile) and soil under earthquake dynamic loading. For instance, the stress-strain hysteresis curve (Figure 3) of soil under earthquake loading demonstrates that its dynamic behavior is strongly nonlinear. The deformation characteristics differ considerably under different stress states, and the nonlinearity of the limit-state equation is also a crucial distinction. Obviously, the nonlinearity means that piles are extremely sensitive to seismic excitation; the seismic dynamic responses of identical laterally loaded piles will

Seismic Dynamic Reliability
In general, a laterally loaded pile subjected to seismic excitation is a strongly nonlinear dynamic system because of the nonlinear properties of concrete (i.e., the pile) and soil under earthquake dynamic loading. For instance, the stress-strain hysteresis curve (Figure 3) of soil under earthquake loading demonstrates that its dynamic behavior is strongly nonlinear. The deformation characteristics differ considerably under different stress states, and the nonlinearity of the limit-state equation is also a crucial distinction. Obviously, the nonlinearity means that piles are extremely sensitive to seismic excitation; the seismic dynamic responses of identical laterally loaded piles will not be exactly the same because of the stochastic nature of the earthquake ground motion. This variation is tremendous and has serious implications for seismic design and performance evaluation of laterally loaded piles. Therefore, the stochastic nature of seismic ground motion and the nonlinearity of material properties should be considered simultaneously in any system involving laterally loaded piles. However, the coupling between nonlinearity and randomness makes it a daunting challenge to solve for and extract the information regarding the stochastic seismic response (e.g., the PDF) of a nonlinear stochastic dynamic system. Nevertheless, the PDF of the seismic dynamic response is the foundation of any reliability calculation [27,28]. Generally speaking, only the PDF of a seismic response that obeys a Gaussian distribution can be defined by two parameters (mean and variance). Unfortunately, for most dynamic responses of a laterally loaded pile, the PDF is either not known a priori or is not Gaussian. In particular, even if the random parameters of the material properties and excitations are normally distributed, the response of the system might not have the same probability distribution as that of the initial stochastic source because of the nonlinear evolution of the system. Moreover, the PDF cannot be so easily available just only depend on the mean and variance. Hence, when analyzing the reliability of a system involving laterally loaded piles either analytically or numerically, it is essential to have access to the PDF. Here, this subsection investigates the stochastic dynamic responses of pile subjected to seismic excitation with random peak ground acceleration (PGA) and stochastic earthquake ground motion, respectively.

The Stochastic Dynamic Reliability of Pile under Seismic Excitation with Random PGA
In this section, we seek the PDF of a laterally loaded pile subjected to stochastic seismic excitation with a peak ground acceleration (PGA) given by a normally distributed random variable.
According to the theory of stochastic dynamics, the random seismic system of a laterally loaded pile can be expressed by the following stochastic differential equation: where ( ) = ( ), ( ), ⋯ , ( ) is the seismic response vector of the pile, and are functions representing the system properties (e.g., nonlinearity), the dynamic stress-strain curve of the soil under the deterministic seismic-acceleration time history is shown in Figure 3, which  However, the coupling between nonlinearity and randomness makes it a daunting challenge to solve for and extract the information regarding the stochastic seismic response (e.g., the PDF) of a nonlinear stochastic dynamic system. Nevertheless, the PDF of the seismic dynamic response is the foundation of any reliability calculation [27,28]. Generally speaking, only the PDF of a seismic response that obeys a Gaussian distribution can be defined by two parameters (mean and variance). Unfortunately, for most dynamic responses of a laterally loaded pile, the PDF is either not known a priori or is not Gaussian. In particular, even if the random parameters of the material properties and excitations are normally distributed, the response of the system might not have the same probability distribution as that of the initial stochastic source because of the nonlinear evolution of the system. Moreover, the PDF cannot be so easily available just only depend on the mean and variance. Hence, when analyzing the reliability of a system involving laterally loaded piles either analytically or numerically, it is essential to have access to the PDF. Here, this subsection investigates the stochastic dynamic responses of pile subjected to seismic excitation with random peak ground acceleration (PGA) and stochastic earthquake ground motion, respectively.

The Stochastic Dynamic Reliability of Pile under Seismic Excitation with Random PGA
In this section, we seek the PDF of a laterally loaded pile subjected to stochastic seismic excitation with a peak ground acceleration (PGA) given by a normally distributed random variable.
According to the theory of stochastic dynamics, the random seismic system of a laterally loaded pile can be expressed by the following stochastic differential equation: where X(t) = [X 1 (t), X 2 (t), · · · , X n (t)] T is the seismic response vector of the pile, f j and g j are functions representing the system properties (e.g., nonlinearity), the dynamic stress-strain curve of the soil under the deterministic seismic-acceleration time history is shown in Figure 3, which indicates that the soil goes into nonlinear state. Therefore, for nonlinear stochastic seismic dynamic response analysis, the coupling of randomness and nonlinearity is inevitable, which is also one of the difficulties in nonlinear dynamic reliability analysis. This also demonstrates that a system involving a laterally loaded pile subjected to earthquake excitation is strongly nonlinear. α is a random variable describing the variability of the PGA, and ϑ(t) is the time history of the deterministic seismic acceleration. Herein, it selects the strong-motion record (Figure 4) from the 1940 El Centro earthquake as the input excitation and modify its PGA by the standard normal random variable α. It is a well-known seismic record and is often selected for nonlinear dynamic response analysis. It can be found from the response spectrum ( Figure 4) that the energy is mainly distributed in the low frequency part of short period. Based on the Increment Dynamic Analysis (IDA) theory, the El Centro acceleration record is selected as the input excitation for the nonlinear dynamic response analysis of pile foundation. It is a well-known seismic record and is often selected for nonlinear dynamic response analysis. It can be found from the response spectrum ( Figure 4) that the energy is mainly distributed in the low frequency part of short period. Based on the Increment Dynamic Analysis (IDA) theory, the El Centro acceleration record is selected as the input excitation for the nonlinear dynamic response analysis of pile foundation. Suppose that ( ) (e.g., displacement, acceleration, moment or shear force) is the important response for the seismic design of a laterally loaded pile. The PDEE can be written as where ( = 1,2, ⋯ , ) are the discrete representative sampling points selected in the probability space. In this section, is treated as being related to a standard normal distribution and it selects roughly one thousand points (i.e., = 1,000 ). By doing so, the stochastic seismic dynamic response analysis of the pile is translated into a series of deterministic dynamic time-history analyses that we conduct by FE analysis in OpenSees. The FE model is shown in Figure 5, where the pile which made by concrete and steel is modeled by dispBeamColumn elements with concrete01 and hardening constitutive, and the soil is modeled by quad elements with nDMaterial constitutive. The corresponding parameters of each constitutive are listed in Table 2, where is the concrete compressive strength, is the concrete strain at maximum strength, , is the concrete crushing strength, is the concrete strain at crushing strength, is the elastic modulus, is the yield stress, and are the isotropic and kinematic hardening moduli, respectively, is the soil mass density, is the low-strain shear modulus, and is the bulk modulus. The soil-pile interface is very important for the dynamic history analysis. The interfacial shear slip between pile and soil,  Suppose that Z(t) (e.g., displacement, acceleration, moment or shear force) is the important response for the seismic design of a laterally loaded pile. The PDEE can be written as where α q (q = 1, 2, · · · , N sel ) are the discrete representative sampling points selected in the probability space. In this section, α is treated as being related to a standard normal distribution and it selects roughly one thousand points (i.e., N sel = 1000). By doing so, the stochastic seismic dynamic response analysis of the pile is translated into a series of deterministic dynamic time-history analyses that we conduct by FE analysis in OpenSees. The FE model is shown in Figure 5, where the pile which made by concrete and steel is modeled by dispBeamColumn elements with concrete01 and hardening constitutive, and the soil is modeled by quad elements with nDMaterial constitutive. The corresponding parameters of each constitutive are listed in Table 2, where f c is the concrete compressive strength, ε c is the concrete strain at maximum strength, f c,u is the concrete crushing strength, ε u is the concrete strain at crushing strength, E is the elastic modulus, σ s is the yield stress, H i and H k are the isotropic and kinematic hardening moduli, respectively, ρ 0 is the soil mass density, G is the low-strain shear modulus, and K is the bulk modulus. The soil-pile interface is very important for the dynamic history analysis. The interfacial shear slip between pile and soil, diameter of pile and non-tensile effect of soil should be considered in principle. A simplified method is adopted in this paper, namely, the translational movement of nodes of soil and pile at the same place is bonded together. The equal displacement boundary conditions are used for the left and right two boundaries to keep them synchronized, which can stimulate the simplified shear boundary, the belief that the dynamic behavior of soil is simple shear movement.

Concrete01
Harding nDMaterial For a series of deterministic nonlinear seismic dynamic time-history analyses by OpenSees (PDEM: 1000 times; Monte Carlo: 20,000 times), the nonlinear seismic dynamic response set of laterally loaded pile is achieved and introduced into the PDEE as the generalized velocity and Monte Carlo stochastic simulation to obtain the abundant probability information (e.g., PDF, mean and variance) of the stochastic dynamic system. The specific solution flowchart of PDEE is shown in Figure 6 as follows.
For a series of deterministic nonlinear seismic dynamic time-history analyses by OpenSees (PDEM: 1000 times; Monte Carlo: 20,000 times), the nonlinear seismic dynamic response set of laterally loaded pile is achieved and introduced into the PDEE as the generalized velocity and Monte Carlo stochastic simulation to obtain the abundant probability information (e.g., PDF, mean and variance) of the stochastic dynamic system. The specific solution flowchart of PDEE is shown in Figure 6 as follows.  To verify the efficiency and accuracy of the PDEM for a stochastic dynamic system, it also conducted 20,000 trials of Monte Carlo stochastic simulation. Figures 7 and 8 demonstrate the high accuracy of PDEM from the perspectives of the second-order statistics and the PDF, respectively, of the pile-top displacement.    The gap of several orders of magnitude between the mean and standard deviation of the displacement indicates the huge variability of the seismic dynamic displacement. This characteristic is also verified by the stochastic temporal fluctuations of the PDF of the pile-top displacement.
Meanwhile, the PDF evolution surface (Figure 9) of the displacement can also be obtained by the PDEM. Figure 9 shows that the PDF evolves as the water stream and rolling hills. It also reflects the fact that the seismic dynamic displacement of the pile fluctuates with time. This is because the statement is controlled by the nonlinear properties of the system, and the evolution of probability The gap of several orders of magnitude between the mean and standard deviation of the displacement indicates the huge variability of the seismic dynamic displacement. This characteristic is also verified by the stochastic temporal fluctuations of the PDF of the pile-top displacement.
Meanwhile, the PDF evolution surface (Figure 9) of the displacement can also be obtained by the PDEM. Figure 9 shows that the PDF evolves as the water stream and rolling hills. It also reflects the fact that the seismic dynamic displacement of the pile fluctuates with time. This is because the statement is controlled by the nonlinear properties of the system, and the evolution of probability information is transmitted by a series of samples. More importantly, the evolution driven by the coupling between nonlinearity and randomness causes the probability distribution of the displacement time history to deviate from the original distribution of the excitation. In other words, because of nonlinearity, the response to a normally distributed stochastic excitation may not be one with a Gaussian distribution. information is transmitted by a series of samples. More importantly, the evolution driven by the coupling between nonlinearity and randomness causes the probability distribution of the displacement time history to deviate from the original distribution of the excitation. In other words, because of nonlinearity, the response to a normally distributed stochastic excitation may not be one with a Gaussian distribution. Finally, the dynamic reliability ( Figure 10) of the laterally loaded pile subjected to seismic excitation is obtained using both the PDEM and Monte Carlo stochastic simulations. It can be seen from the figure that different performance indicators correspond to different safety assurance rates, that is, the so-called seismic dynamic reliability of pile foundation. The validity and high accuracy of PDEM method are verified by Monte Carlo method from the aspect of seismic dynamic reliability. Furthermore, to verify the efficiency of the PDEM, the computation times for static and dynamic reliability assessments are compared in Table 3 for the PDEM and Monte Carlo simulation. Table 3 demonstrates that the PDEM is significantly more computationally efficiency compared with the Finally, the dynamic reliability ( Figure 10) of the laterally loaded pile subjected to seismic excitation is obtained using both the PDEM and Monte Carlo stochastic simulations. It can be seen from the figure that different performance indicators correspond to different safety assurance rates, that is, the so-called seismic dynamic reliability of pile foundation. The validity and high accuracy of PDEM method are verified by Monte Carlo method from the aspect of seismic dynamic reliability. Furthermore, to verify the efficiency of the PDEM, the computation times for static and dynamic reliability assessments are compared in Table 3 for the PDEM and Monte Carlo simulation. Table 3 demonstrates that the PDEM is significantly more computationally efficiency compared with the classic Monte Carlo stochastic simulation. It should be noted that the Monte Carlo simulation method used in this paper is simple and plain without any other sampling techniques such as Latin hypercube. Here, the effect of different numbers (5000, 10,000, and 20,000) of simulations on results is compared with that of PDEM by which it can find the good comparison is observed when the simulation number of Monte Carlo is 20,000. Finally, the dynamic reliability ( Figure 10) of the laterally loaded pile subjected to seismic excitation is obtained using both the PDEM and Monte Carlo stochastic simulations. It can be seen from the figure that different performance indicators correspond to different safety assurance rates, that is, the so-called seismic dynamic reliability of pile foundation. The validity and high accuracy of PDEM method are verified by Monte Carlo method from the aspect of seismic dynamic reliability. Furthermore, to verify the efficiency of the PDEM, the computation times for static and dynamic reliability assessments are compared in Table 3 for the PDEM and Monte Carlo simulation. Table 3 demonstrates that the PDEM is significantly more computationally efficiency compared with the classic Monte Carlo stochastic simulation. It should be noted that the Monte Carlo simulation method used in this paper is simple and plain without any other sampling techniques such as Latin hypercube. Here, the effect of different numbers (5000, 10,000, and 20,000) of simulations on results is compared with that of PDEM by which it can find the good comparison is observed when the simulation number of Monte Carlo is 20,000.   In the previous section, we investigated and discussed the stochastic seismic dynamic response and dynamic reliability of pile under the seismic action with random PGA. However, for practical engineering application, the future earthquake in the engineering site of the pile cannot be accurately predicted. Moreover, the input seismic excitation may not purely have the randomness in PGA. It should have the stochastic characteristics in both intensity and frequency. Therefore, in this section, we will continue to discuss the stochastic seismic dynamic reliability of pile subjected to stochastic seismic ground motion. Here, the stochastic dynamic difference equation Equation (20) can be modified as follows: where ..
U g (Θ, t) is the stochastic seismic excitation and Θ is the random vector which describes the randomness in intensity and frequency. The detailed generation methodology of seismic ground motion time history samples can be found in the relevant references [29,30]. For the stochastic seismic dynamic analysis, there are 254 deterministic dynamic time-history analyses in OpenSees with the FE model in Figure 5, and the typical acceleration time history samples are shown in Figure 11.
It should be noted that the time histories set of the ground motion acceleration sample is the intensity frequency non-stationary random ground motion model of the same set system with corresponding given probability. It is determined based on the spectral representation of the corresponding site power spectral density function and random function. Its validity and rationality have been verified in a series of seismic dynamic evaluation of geotechnical engineering [30]. Moreover, the acceleration sample time histories in Figure 11 are only 3 out of 254 sample histories in the same set. Similarly, it also can obtain the one-dimensional PDEE Equation (8) of the key response quantity which impacts the seismic design. The CDF of displacement of laterally loaded pile subjected to stochastic seismic ground motion is shown in Figure 12.
where ( , ) is the stochastic seismic excitation and is the random vector which describes the randomness in intensity and frequency. The detailed generation methodology of seismic ground motion time history samples can be found in the relevant references [29,30]. For the stochastic seismic dynamic analysis, there are 254 deterministic dynamic time-history analyses in OpenSees with the FE model in Figure 5, and the typical acceleration time history samples are shown in Figure 11. It should be noted that the time histories set of the ground motion acceleration sample is the intensity frequency non-stationary random ground motion model of the same set system with corresponding given probability. It is determined based on the spectral representation of the corresponding site power spectral density function and random function. Its validity and rationality have been verified in a series of seismic dynamic evaluation of geotechnical engineering [30]. Moreover, the acceleration sample time histories in Figure 11 are only 3 out of 254 sample histories in the same set. Similarly, it also can obtain the one-dimensional PDEE Equation (8) of the key response quantity which impacts the seismic design. The CDF of displacement of laterally loaded pile subjected to stochastic seismic ground motion is shown in Figure 12.

Conclusions
In this paper, the static and dynamic (seismic) reliabilities of a laterally loaded pile are evaluated in light of the stochastic dynamics and a series of deterministic analyses. The innovation of this study is to introduce a new stochastic analysis method and high-dimensional spatial data reduction technology for obtaining the probability density function of multivariable stochastic static and seismic dynamic stochastic systems. The TS sampling technique is introduced to reduce the necessary sample numbers and calculation quantity. By combining TS sampling with the PDEM, the static reliability is obtained by integrating the PDF. The dynamic reliability is analyzed based on the PDEM and a series of deterministic time-history analyses using the FE technique. From comparisons with Monte Carlo stochastic simulation, the efficiency and accuracy of PDEM is demonstrated. The PDEM is clearly more efficient in assessing static and dynamic responses and reliability of a laterally loaded pile; when analyzing the static and dynamic reliabilities of a pile, it is roughly 38 and 20 times more

Conclusions
In this paper, the static and dynamic (seismic) reliabilities of a laterally loaded pile are evaluated in light of the stochastic dynamics and a series of deterministic analyses. The innovation of this study is to introduce a new stochastic analysis method and high-dimensional spatial data reduction technology for obtaining the probability density function of multivariable stochastic static and seismic dynamic stochastic systems. The TS sampling technique is introduced to reduce the necessary sample numbers and calculation quantity. By combining TS sampling with the PDEM, the static reliability is obtained by integrating the PDF. The dynamic reliability is analyzed based on the PDEM and a series of deterministic time-history analyses using the FE technique. From comparisons with Monte Carlo stochastic simulation, the efficiency and accuracy of PDEM is demonstrated. The PDEM is clearly more efficient in assessing static and dynamic responses and reliability of a laterally loaded pile; when analyzing the static and dynamic reliabilities of a pile, it is roughly 38 and 20 times more efficient, respectively, than Monte Carlo stochastic simulation.
The results from the stochastic dynamic system of a laterally loaded pile subjected to seismic excitation show that the probability information is transported strongly by the nonlinear state evolution. They also indicate that the probability information of the seismic dynamic response of a pile cannot be easily determined from the excitation distribution. The assumption applied in traditional reliability analysis, namely, that the response and excitation have the same similar distribution, should be reconsidered.
Different from the classical traditional reliability analysis, especially the random analysis based on Monte Carlo stochastic simulation, the main feature of this study is to obtain the probability density function of nonlinear static and dynamic pile foundation system. Combining with the new high-dimensional probability space dimension reduction sampling technology such as TS sampling, the calculation amount of stochastic analysis, especially the nonlinear seismic dynamic reliability analysis, can be effectively reduced.

Conflicts of Interest:
The authors declare no conflict of interest.