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.
3.1. 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
and
, 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
. The material properties of the pile are described by an elastic modulus of
and a section moment of inertia of
. 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
and moment
is expressed by Equations (10) and (11):
where
is a parameter expressed as
Therefore, the performance index of the pile is selected as .
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
,
, and
. 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
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
is the layer number of the TS sample technology in the
layer. We assume that
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
can be taken as
, where
represents the maximum integer that is no more than the value in the brackets and
, where
is the radius of the TS and l is the number of TS layers around the sphere in the center in the
plane.
If the sphere-center coordinates of the
-layer spheres are
, then
where
and
, with
. Coordinates
are those of the
plane and can be determined by
where
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
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
should meet the following requirement:
Here, the coordinates are renumbered as
,
. Distribute these sampling points around the
axis with fixed angle
. to make the sampling more homogeneous, therefore:
Accordingly, the assigned probability of the representative sampling points can be calculated as
where
is the volume enclosed by the 12 tangent planes and
.
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
,
, and
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
is the ratio of the correlation coefficients before and after transformation,
is the correlation coefficient of the target distribution, and
is the original correlation coefficient between
and
, namely,
. Beyond that, however, the correlation coefficient can also be calculated numerically as
where
and
are the mean and standard deviation, respectively, of random variable
. Herein, random variables
and
correspond to
and
, respectively, with
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
sampling points in the 3D probability space formed by the random variables
,
, and
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
, 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.
3.2. 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.
3.2.1. 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 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
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.
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
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.,
). 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, 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.
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.
Figure 7 and
Figure 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 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 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.
3.2.2. The Stochastic Dynamic Reliability of Pile Subjected to Stochastic Seismic Ground Motion
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
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.