Next Article in Journal
Approaches to Assessing the Strategic Sustainability of High-Risk Offshore Oil and Gas Projects
Next Article in Special Issue
Study on Reactive Power Compensation Strategies for Long Distance Submarine Cables Considering Electrothermal Coordination
Previous Article in Journal
Multi-Step-Ahead Forecasting of Wave Conditions Based on a Physics-Based Machine Learning (PBML) Model for Marine Operations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Department of Geotechnical Engineering, College of Civil Engineering, Tongji University, Shanghai 200092, China
2
Key Laboratory of Geotechnical and Underground Engineering of the Ministry of Education, Tongji University, Shanghai 200092, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2020, 8(12), 994; https://doi.org/10.3390/jmse8120994
Submission received: 17 November 2020 / Revised: 1 December 2020 / Accepted: 3 December 2020 / Published: 5 December 2020

Abstract

:
Pile foundation is one of the common foundation forms in marine geotechnical engineering, especially in wind power engineering. Its operation safety is seriously affected 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 difficult 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 efficient when applied to the seismic design of pile in geotechnical engineering, and its calculation efficiency is 20 times of the former for the seismic dynamic reliability of pile foundation. This study provides a new reference for the efficient design and safety evaluation of offshore pile foundation engineering based on static and dynamic reliability of multiple random factors.

1. 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.

2. 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
U = Γ ( Θ , ξ , u , t ) ,
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
U r = Γ r ( Θ , ξ , u , t ) ( r = 1 , 2 , , n ) ,
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 ξ :
U ˙ r ( ξ ) = U r ξ ( r = 1 , 2 , , n ) .
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,
d d ξ Ω ξ × Ω Θ p U Θ ( u , θ , ξ ) d u d θ = 0 ,
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
d d ξ Ω ξ × Ω Θ p U Θ ( u , θ , ξ ) d u d θ = d d ξ Ω ξ 0 × Ω Θ p U Θ ( u , θ , ξ ) | J | d u d θ = Ω ξ 0 × Ω Θ p U Θ ( p U Θ ξ + r = 1 n U ˙ r p U Θ u r ) | J | d u d θ . = Ω ξ × Ω Θ p U Θ ( p U Θ ξ + r = 1 n U ˙ r p U Θ u r ) d u d θ
Because of the arbitrariness of Ω ξ × Ω Θ , the generalized probability density evolution equation (PDEE) can be obtained from Equations (3) and (4):
p U Θ ( u , θ , ξ ) ξ + r = 1 n p U Θ ( u , θ , ξ ) u r = 0
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
p U ( u , ξ ) = Ω Θ p U Θ ( u , θ , ξ ) d θ .
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,
p U r Θ ( u r , θ , ξ ) ξ + U ˙ r p U r Θ ( u r , θ , ξ ) u r = 0 .
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:
p U Θ ( u , θ , ξ ) ξ = 0 = δ ( u u 0 ) p Θ ( θ ) ,
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 s e l ) in each subdivision, where N s e l 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 Section 3.1 and Section 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.

3. 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.

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 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.
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):
y ( z ) = 2 F λ k h d { sinh ( λ L ) cos ( λ L ) cosh [ λ ( L z ) ] sin ( λ L ) cosh ( λ z ) cos [ λ ( L z ) ] sinh 2 ( λ L ) sin 2 ( λ L ) } 2 M λ 2 k h d { sinh ( λ L ) [ cos ( λ L ) cos h ( λ ( L z ) ) sin ( λ z ) sinh ( λ ( L z ) ) cos ( λ z ) ] sinh 2 ( λ L ) sin 2 ( λ L ) + sin ( λ L ) [ sinh ( λ z ) cos h ( λ ( L z ) ) sin ( λ z ) cosh ( λ z ) cos ( λ ( L z ) ) ] sinh 2 ( λ L ) sin 2 ( λ L ) }
where λ is a parameter expressed as
λ = ( k h d 4 E p I p ) 1 4 .
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:
z k = k Δ z ( k = 0 , ± 1 , , ± N ) ,
where N z is the layer number of the TS sample technology in the z t h 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 i n t ( η Δ z ) + 1 , where i n t ( · ) represents the maximum integer that is no more than the value in the brackets and z = 2 6 3 R , 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
{ x i , k = x i , 0 + 1 2 [ 1 ( 1 ) k ] x y i , k = y i , 0 + 1 2 [ 1 ( 1 ) k ] y , ( k = 0 , ± 1 , ± 2 , , ± N z ) ,
where x = 2 3 R sin ψ 3 and x = 2 3 R cos ψ 3 , with ψ = π 12 . Coordinates ( x i , 0 , y i , 0 , z k ) are those of the z = 0 plane and can be determined by
x i , 0 = r i cos φ , y i , 0 = r i sin φ ,
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:
x i , k 2 + y i , k 2 + z i , k 2 η 2 .
Here, the coordinates are renumbered as ( x q , y q , z q ) , q = 1 , 2 , , N s e l . Distribute these sampling points around the z axis with fixed angle ψ = π 12 . to make the sampling more homogeneous, therefore:
{ θ q 1 = x q θ q 2 = y q cos ψ + z q sin ψ θ q 3 = z q cos ψ y q sin ψ .
Accordingly, the assigned probability of the representative sampling points can be calculated as
P q = V q N s e l j = 1 N s e l p ( θ q 1 , j , θ q 2 , j , θ q 3 , j ) , j = 1 , 2 , , N s e l ,
where V q is the volume enclosed by the 12 tangent planes and θ q i , j = θ q i + ( 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]:
R c o e = ρ Y i Y j ρ X i X j = 1.064 + ( 0.069 ) ρ X i X j + 0.005 ρ X i X j 2 .
where R c o e 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
ρ X i X j = + + x i μ X i σ X i x j μ X j σ X j φ 2 ( y i , y j ρ Y i Y j ) d y i d y j ,
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 s e l = 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.

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:
d d t X j ( t ) = f j [ X ( t ) , t ] + g j [ X ( t ) , t ] α ϑ ( t ) ,
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.
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
p Z ( z , α q , t ) t + Z ˙ p Z ( z , α q , t ) z = 0 ,
where α q ( q = 1 , 2 , , N s e l ) 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 s e l = 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.
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:
d d t X j ( t ) = f j [ X ( t ) , t ] + g j [ X ( t ) , t ] U ¨ g ( Θ , t ) ,
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.

4. 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.

Author Contributions

M.X. and Y.H. are responsible for the preparation of this study, M.X. is responsible for performance of all the computations, the collection of data, and all the passages in this manuscript, while M.X. and Y.H. revised the paper. This study was supervised by Y.H. The results were discussed and the conclusions drafted jointly by all authors. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Fund of China, grant number 41902274 and the National Key R&D Program of China, grant number 2017YFC1501304.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vogt, N.; Vogt, S.; Kellner, C. Buckling of slender piles in soft soils. Bautechnik 2009, 86, 98–112. [Google Scholar] [CrossRef]
  2. Yang, J.; Ye, J.Q. Dynamic elastic local buckling of piles under impact loads. Struct. Eng. Mech. 2002, 13, 543–556. [Google Scholar] [CrossRef]
  3. Liang, F.; Li, Y.; Li, L.; Wang, J. Analytical solution for laterally loaded long piles based on Fourier-Laplace integral. Appl. Math. Model. 2004, 38, 5198–5216. [Google Scholar] [CrossRef]
  4. Bolton, M.D. An alternative mechanism of pile failure in liquefiable deposits during earthquakes. Geotechnique 2004, 54, 203–213. [Google Scholar]
  5. Mehta, A.; Pachpor, P.D. Behavior of laterally loaded pile. Int. J. Eng. Sci. 2010, 2, 77–86. [Google Scholar]
  6. Han, F.; Salgado, R.; Prezzi, M. Nonlinear analyses of laterally loaded piles-A semi-analytical approach. Comput. Geotech. 2015, 70, 116–129. [Google Scholar] [CrossRef] [Green Version]
  7. Zhang, J.; Zhang, L.M.; Tang, W.H. Reliability based design of pile foundations considering both parameter and model uncertainties. J. Geoeng. 2009, 4, 119–127. [Google Scholar]
  8. Barakat, S.A.; Malkawi, A.I.H.; Ra’ed, H.T. Reliability-based optimization of laterally loaded piles. Struct. Saf. 1999, 21, 45–64. [Google Scholar] [CrossRef]
  9. Tandjiria, V.; Teh, C.I.; Low, B.K. Reliability analysis of laterally loaded piles using response surface methods. Struct. Saf. 2000, 22, 335–355. [Google Scholar] [CrossRef]
  10. Kim, S.B.; Yoon, G.L.; Yi, J.H.; Lee, J.H. Reliability analysis of laterally loaded piles for an offshore wind turbine support structure using response surface methodology. Wind Struct. 2015, 21, 597–607. [Google Scholar] [CrossRef]
  11. Haldar, S.; Babu, G.L.S. Design of laterally loaded piles in clays based on cone penetration test data: A reliability-based approach. Geotechnique 2009, 59, 593–607. [Google Scholar] [CrossRef]
  12. Low, B.K.; Tang, W.H. Efficient spreadsheet algorithm for first-order reliability method. J Eng. Mech. 2007, 133, 1378–1387. [Google Scholar] [CrossRef]
  13. Rad, M.M. Reliability based analysis and optimum design of laterally loaded piles. Period. Polytech. Civ. Eng. 2017, 61, 491–497. [Google Scholar]
  14. Schuëller, G.I. A state-of-the-art report on computational stochastic mechanics. Probab. Eng. Mech. 1997, 12, 197–321. [Google Scholar] [CrossRef]
  15. Engelund, S.; Rackwitz, R. A benchmark study on importance sampling techniques in structural reliability. Struct. Saf. 1993, 12, 255–276. [Google Scholar] [CrossRef]
  16. Ditlevsen, O.; Madsen, H.O. Structural Reliability Methods; Wiley: New York, NY, USA, 1996; pp. 331–335. [Google Scholar]
  17. Kawano, K.; Venkataramana, K. Dynamic response and reliability analysis of large offshore structures. Comput. Methods Appl. Mech. Eng. 1999, 168, 255–272. [Google Scholar] [CrossRef]
  18. Unsever, Y.S.; Matsumoto, T.; Esashi, K.; Kobayashi, S. Behaviour of model pile foundations under dynamic loads in saturated sand. Bull. Earthq. Eng. 2017, 15, 1355–1373. [Google Scholar] [CrossRef]
  19. Chen, J.B.; Li, J. Stochastic seismic response analysis of structures exhibiting high nonlinearity. Comput. Struct. 2010, 88, 395–412. [Google Scholar] [CrossRef]
  20. Zhang, Y.; Lam, J.S.L. Reliability analysis of offshore structures within a time varying environment. Stoch. Environ. Res. Risk. Assess. 2015, 29, 1615–1636. [Google Scholar] [CrossRef]
  21. Chen, Q.; Wang, C.; Juang, C.H. Probabilistic and spatial assessment of liquefaction-induced settlements through multiscale random field models. Eng. Geol. 2016, 211, 135–149. [Google Scholar] [CrossRef]
  22. Li, J.; Chen, J.B. The principle of preservation of probability and the generalized density evolution equation. Struct. Saf. 2008, 30, 65–77. [Google Scholar] [CrossRef]
  23. Li, J.; Chen, J.B. Probability density evolution method for dynamic response analysis of structures with uncertain parameters. Comput. Mech. 2004, 34, 400–409. [Google Scholar] [CrossRef]
  24. Chen, J.B.; Li, J. Dynamic response and reliability analysis of non-linear stochastic structures. Probab. Eng. Mech. 2005, 20, 33–44. [Google Scholar] [CrossRef]
  25. Hetényi, M. Beams on Elastic Foundation: Theory with Applications in the Fields of Civil and Mechanical Engineering; University of Michigan Press: Ann Arbor, MI, USA, 1952; pp. 52–55. [Google Scholar]
  26. Chan, C.L.; Low, B.K. Probabilistic analysis of laterally loaded piles using response surface and neural network approaches. Comput. Geotech. 2012, 43, 101–110. [Google Scholar] [CrossRef]
  27. Wang, L.; Hwang, J.H.; Juang, C.H.; Atamturktur, S. Reliability-based design of rock slopes-A new perspective on design robustness. Eng. Geol. 2013, 154, 56–63. [Google Scholar] [CrossRef]
  28. Juang, C.H.; Ching, J.; Luo, Z.; Ku, C.S. New models for probability of liquefaction using standard penetration tests based on an updated database of case histories. Eng. Geol. 2012, 133, 85–93. [Google Scholar] [CrossRef]
  29. Liu, Z.J.; Liu, W.; Peng, Y.B. Random function based spectral representation of stationary and non-stationary stochastic processes. Probab. Eng. Mech. 2016, 45, 115–126. [Google Scholar] [CrossRef]
  30. Huang, Y.; Xiong, M. Probability density evolution method for seismic liquefaction performance analysis of earth dam. Earthq. Eng. Struct. Dyn. 2017, 46, 925–943. [Google Scholar] [CrossRef]
Figure 1. Laterally loaded pile.
Figure 1. Laterally loaded pile.
Jmse 08 00994 g001
Figure 2. (a) Probability density function (PDF) of horizontal displacement of pile top. (b) Cumulative distribution function (CDF; i.e., reliability) of horizontal displacement of pile top.
Figure 2. (a) Probability density function (PDF) of horizontal displacement of pile top. (b) Cumulative distribution function (CDF; i.e., reliability) of horizontal displacement of pile top.
Jmse 08 00994 g002
Figure 3. Stress–strain curve of soil under seismic loading.
Figure 3. Stress–strain curve of soil under seismic loading.
Jmse 08 00994 g003
Figure 4. Seismic record (upper) and spectrum (lower) of the 1940 El Centro earthquake (NS).
Figure 4. Seismic record (upper) and spectrum (lower) of the 1940 El Centro earthquake (NS).
Jmse 08 00994 g004
Figure 5. Finite-element (FE) model of laterally loaded pile.
Figure 5. Finite-element (FE) model of laterally loaded pile.
Jmse 08 00994 g005
Figure 6. Computational process of PDEM.
Figure 6. Computational process of PDEM.
Jmse 08 00994 g006
Figure 7. Mean and variance time histories of pile displacement.
Figure 7. Mean and variance time histories of pile displacement.
Jmse 08 00994 g007
Figure 8. Comparing between probability density evolution method (PDEM) and Monte Carlo of PDFs of displacement at different times.
Figure 8. Comparing between probability density evolution method (PDEM) and Monte Carlo of PDFs of displacement at different times.
Jmse 08 00994 g008
Figure 9. PDFs evolution surface of displacement.
Figure 9. PDFs evolution surface of displacement.
Jmse 08 00994 g009
Figure 10. CDF of displacement of laterally loaded pile subjected to seismic excitation with random peak ground acceleration (PGA).
Figure 10. CDF of displacement of laterally loaded pile subjected to seismic excitation with random peak ground acceleration (PGA).
Jmse 08 00994 g010
Figure 11. Typical time history samples of the seismic excitation input.
Figure 11. Typical time history samples of the seismic excitation input.
Jmse 08 00994 g011
Figure 12. CDF of displacement of laterally loaded pile subjected to stochastic seismic ground motion.
Figure 12. CDF of displacement of laterally loaded pile subjected to stochastic seismic ground motion.
Jmse 08 00994 g012
Table 1. Probabilistic information on the random parameters for static reliability.
Table 1. Probabilistic information on the random parameters for static reliability.
ParameterDistributionMeanCOVCorrelation Coefficient
F ( kN ) Gumbel2400.10.5
M ( kN · m ) Gumbel800.1
k h ( kN / m 3 ) Log-normal10,0000.2Independent
Table 2. Parametric properties of the constitutive.
Table 2. Parametric properties of the constitutive.
Concrete01HardingnDMaterial
f c ( kPa ) −34485.6 E ( kPa ) 2.0 × 108 ρ ( kg / m 3 ) 2.0 × 103
ε c −0.004 σ s ( kPa ) 248,200.0 G ( kPa ) 54,550
f c , u ( kPa ) −20691.4 H i ( kPa ) 0.0 K ( kPa ) 1.6 × 105
ε u −0.014 H k ( kPa ) 1.6129 × 106--
Table 3. Comparison of computation amounts and times between PDEM and Monte Carlo simulation.
Table 3. Comparison of computation amounts and times between PDEM and Monte Carlo simulation.
Analysis TypePDEMMonte Carlo SimulationRatio
Samples/Time (h)Samples/Time (h)
Static52720,00038
Dynamic1000/5.6 h20,000/166.67 h20
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xiong, M.; Huang, Y. Static and Dynamic Reliability Analysis of Laterally Loaded Pile Using Probability Density Function Method. J. Mar. Sci. Eng. 2020, 8, 994. https://doi.org/10.3390/jmse8120994

AMA Style

Xiong M, Huang Y. Static and Dynamic Reliability Analysis of Laterally Loaded Pile Using Probability Density Function Method. Journal of Marine Science and Engineering. 2020; 8(12):994. https://doi.org/10.3390/jmse8120994

Chicago/Turabian Style

Xiong, Min, and Yu Huang. 2020. "Static and Dynamic Reliability Analysis of Laterally Loaded Pile Using Probability Density Function Method" Journal of Marine Science and Engineering 8, no. 12: 994. https://doi.org/10.3390/jmse8120994

APA Style

Xiong, M., & Huang, Y. (2020). Static and Dynamic Reliability Analysis of Laterally Loaded Pile Using Probability Density Function Method. Journal of Marine Science and Engineering, 8(12), 994. https://doi.org/10.3390/jmse8120994

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop