Abstract
In this study, the quantitative effect of the random distribution of the soil material properties to the probability density functions of the failure load and displacements is presented. A modified Cam-Clay failure criterion is embedded into a stochastic finite element numerical tool for this purpose. Various assumptions for the random distribution of the compressibility factor \(\kappa\), of the constitutive relation, the critical state line inclination c of the soil, and the permeability k have been tested and assessed with Latin hypercube sampling followed by Monte Carlo simulation. It is confirmed that both failure load and displacements follow Gaussian normal distribution despite the nonlinearity of the problem. Moreover, as the soil depth increases the mean value of failure load decreases and the failure displacement increases. Consequently, failure mechanism of clays can be determined in this work within an acceptable variability, taking into account the soil depth and nonlinear constitutive relations which in the analytical solutions is not feasible as it is assumed the Meyerhoff theory which considers the elastic halfspace.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
The limit load to shallow foundations and the corresponding displacements and material states is one of the most investigated areas in the field of Geomechanics. Kötter and Prandtl [26, 45] were the first to investigate the problem providing analytical solutions to sliding line equations and bearing capacity in the context of the infinite halfspace soils. [35, 60] investigated the bearing capacity of soils taking into account the friction angle, the soil depth and soil density. Subsequently, a number of papers have been published which investigate the problem [18, 29, 38, 58, 66]taking into account the failure mechanism and applying Mohr–Coulomb material yield criterion alongside with linear elastic analysis. In most cases 1D or 2D physical problems have been examined taking into account homogeneous or layered soils: [40, 47]. All studies presented were applied to the foundation design regulations by introducing the variables N and S which are for material influence and shape influence, respectively. The material related variables N are three respecting for the 3 addends of the total bearing capacity. \(N_q\), \(N_c\), \(N_{\gamma }\) are defined for the influence of possible vertical load in the lateral of the foundation, the cohesion of the soil and the footing dimensions alongside with the total weight of the soil, respectively. In a similar manner the shape variables \(S_q\), \(S_c\), \(S_{\gamma }\) are determined.( [31, 36, 37])
The uncertainty quantification of the bearing capacity of saturated porous media in relation with the input variability such as the Young Modulus or the permeability is investigated with the implementation of the stochastic finite element method [17, 44]. In the aforementioned method, the input variability can be simulated by using two alternative methods [8, 12, 20, 23, 30, 32, 33, 41, 52]. The first method consists of the nodal points to be assumed as random variables and deterministic shape function are implemented for the interpretation of the material spatial distribution [5, 14, 28, 58]. An alternative approach is the implementation of random field series expansion such as the spectral representation or the Karhunen–Loeve expansion or the spatial average method for providing the material input variability [1, 4, 13, 41]. In both approaches, the sampling can be pseudorandom and non-biased or importance sampling methods can be implemented such as the Latin Hypercube Sampling (LHS) [3, 53]. Then, for both methods, the standard Monte Carlo simulation can be implemented. In some cases there are studies that implement realizations of non-Gaussian random fields to the input material variables [43].
In the related studies the problems and methodology proposed refer to 1D-2D elastic halfspace theory with material criterion the Mohr–Coulomb yield stress model. In the present work a numerical simulation approach with a modified Cam-Clay material yield model proposed by [24] is adopted, which is an accurate and reliable quantitatively material for clayey soil simulation ( [63]).This material model with the implementation of the finite element method can be used in general 3D loading and geometric circumstances. The excessive computational cost which is required by the Crude Monte Carlo simulation method, can be reduced by efficient computational methods as suggested in [55, 56]. For the calculation of failure load, a recurrence relation algorithm which provides accurate results with a relative small amount of trials and with one initial guess trial is implemented. The algorithm is theoretically set, proved and compared with the classic intersection method. The goal of the present paper, which is an extension of previous work by the authors [51], is to quantify the uncertainty of the failure load and displacement alongside with the uncertainty quantification of the failure mechanism in relation to the variability of the input parameters, i.e., the spatial distribution of the material variables and the soil depth. Soil depth is defined as the distance between the upper layer of the soil and the impermeable and rigid bedrock. As the soil depth increases the total stiffness of the soil mass is reduced, leading to smaller values for failure loads. Moreover, with the increase of the soil depth the Meyerhoff spline can be expanded to a more generalized failure mechanism in comparison with localized failure surfaces or punching shear failure( [15, 19, 61, 62]).
The input variability refers to three material variables, namely the compressibility factor \(\kappa\), the permeability k and the critical state line inclination c of the soil. The possible spatial distributions of the material variables with respect to the depth Z of the soil domain are the constant, the linear variation and the random field process calculated through the Karhunen–Loeve series expansion with an exponential autocorrelation function. In the constant and linear variation with respect to depth, the truncated normal random variable [7, 49] is assumed at the nodal points. It should be noted that LHS importance sampling method is adopted for the construction of the input random vectors. The geometrical dimensions of the 3D soil domain and correlation lengths are parametrically studied and compared with respect to the corresponding solid problem without considering the pore pressure.
2 Formulation of the dynamic soil–pore–fluid interaction and the numerical solution
The Biot system of equations governs the physical behavior of porous media static or dynamic loading. If the dynamic load is of low frequency or if static loading problems are analyzed the aforementioned system is reduced leading to a notable decrease of the computational effort. The u-p formulation of Biot system of equations consists of soil-fluid momentum balance with the Darcian flow, the stress–strain law and the boundary conditions. In the present article static forces are applied to the clay soil domain consequently the u-p formulation is adopted.
The finite element discretization of the u-p formulation takes the form:
where
where \({\mathbf {k}}\) is the matrix of permeability n units \([length]^3 [time]/[mass]\), \({\mathbf {b}}\) is the loading vector divided by the total mixture density and Q is a factor influenced by the bulk moduli of fluid and soil skeleton. \(\mathbf {N^p}\) are the shape functions of pore pressure in matrix representation. \(\mathbf {M_1 , C_1 , K_1}\) are the standard mass, damping and stiffness matrices of the solid skeleton. Furthermore, \(\mathbf {Q_c , H, S}\) are the coupling, permeability and saturation matrices, respectively. Finally, \(\mathbf {F_1}\) corresponds to the equivalent forces due to the external loading. Numerical schemes, such as the Newmark direct integration method, are implemented to obtain the solution to the problem.
3 The yield stress model
3.1 Definition of plastic and bond strength envelope
The material model,which here is presented in a brief description, is a modified Cam-Clay type model. All the stresses in the following equations of Sect. 3.1 are effective stresses in the solid skeleton. The implemented model is a two surface model, namely the plastic yield envelope (PYE) which defines the elastic region and the bond strength envelope (BSE) which defines the region in which PYE can be located as PYE is always inside BSE. [9, 10, 22, 24, 63]. BSE is directly influenced by the structure of the microtiles of the clay. If a stress point lies in BSE, the structure degradation rate of the clay is maximum. Both envelopes have ellipsoidal shape as depicted in Fig. 1. The mathematical representation of BSE is given by the following:
where \(p_h\) is the hydrostatic component of the stress tensor, s is the deviatoric component of the stress tensor, c is the critical state line inclination and a is the half-size of the large diameter of the ellipse
Furthermore, PYE is described by the function
where \(p_L\) is the hydrostatic component of the centroid of PYE , \(\mathbf {s_L}\) is the deviatoric component of the centroid of PYE and \(\xi\) is the similarity factor
Finally, inside PYE the isotropic poroelasticity theory applies, where the bulk and shear moduli have constant ratio, with the assumption of constant Poisson ratio. The bulk modulus is given by:
where \(\nu\) is the specific volume of the soil.
The constitutive model is valid for clays, for static and dynamic loads, regardless of the overconsolidation ratio. Clays of friction angle between 17\(^o\)-30\(^o\) can be simulated with substantial accuracy. This range of friction angles corresponds to the majority of the natural clayey soils. Furthermore, it is proved to be numerically stable because the majority of the equations of the criterion are in a closed form. Furthermore, when highly consolidated clays are investigated it can be easily transformed to take into account a possible cementation of the soil, i.e., to bear tensile stresses and the BSE to move to the left. It should also be noted that the stresses and the strains transformations are made in order to have energy conjugate amounts by applying the numerical transformations used in Von Mises yield criterion. To this context, the following variables are introduced \(q_s=\sqrt{\frac{3}{2}\mathbf {s:s}}=\sqrt{3J_2}\) and \(e=\sqrt{\frac{2}{3}\mathbf {\epsilon _{dev}:\epsilon _{dev}}}\). The variable \(q_s\) stands for the deviatoric stress measure also known as Von Mises stress, \(J_2\) is the second invariant of the deviatoric part of the stress vector and e stands for deviatoric strain measure.
4 Random fields, the truncated normal variables and the Latin hypercube sampling
4.1 The Karhunen–Loeve series and the truncated normal distribution
The input material variability can be simulated either assuming that the nodal points are random variables with deterministic shape functions, or by implementing the Karhunen–Loeve series expansion for the construction or a stochastic process [17, 23, 27, 39, 42, 46, 64, 65]. In the present study, both approaches are adopted in order to evaluate each approach influence. For the first technique, which is proposed in [64] the random function f is estimated with the usage of shape functions \(N_i\) by
where \(N_0\) is the total number of shape functions and the \(f_i\) are the values of f in the nodal points which can be random variables following a probability density function (PDF). In the present work, \(N_i\) are linear functions and the \(f_i\) follow the truncated normal distribution [2, 6, 7, 49] which has a pdf described by the equation \(g_1(x)=\frac{\phi (X_0)}{{\sigma _d} (\Phi (B)-\Phi (A))}\), where \(\phi (X_0)\) and \(\Phi (X_0)\) are the standard normal probability and cumulative distribution function for \(X_0\), respectively, and A, B, \(X_0\) are the normalized coordinates of the subspace limits and x, respectively. Also, \(\sigma _d\) and the mean value of the normalization refers to the pdf of the random variable before truncation.
For the Karhunen–Loeve implementation, let \(H_1(\mathbf{x } , \omega )\) be a random field of mean \(\mu ({\mathbf {x}})\) based on a known autocovariance function \(C_h(\mathbf {x_1},\mathbf {x_2})=\sigma (\mathbf {x_1}) \sigma (\mathbf {x_2}) \rho (\mathbf {x_1},\mathbf {x_2})\), where \(\rho (\mathbf {x_1},\mathbf {x_2})\) is the correlation function and \(\sigma (\mathbf {x_1})\) is the standard deviation of \(\mathbf {x_1}\). Therefore, any realization \(H_1\) of the field with M number of eigenfunctions \(\phi _i\) with corresponding eigenvalues \(\lambda _i\) can be expanded as:
where \(\xi _i\) is a set of random variables of zero mean and covariance function \(COV(\xi _i ,\xi _j)=\delta _{ij}\). Finally, for a Gaussian random field, as implemented in the present study, the \(\xi _i\) functions are a set of standard normal random variables. This type of expansion is the most common because it is strong and robust. In the present work the Karhunen–Loeve expansion is applied by adopting an exponential autocovariance function which has an analytical solution of the Fredholm eigenfunction problem [23, 46]
where b is the correlation length. If the Fredholm equations cannot be solved analytically, and this occurs when the autocovariance function is more complicated, numerical methods can be applied such as the Galerkin method [20, 46]. In the present work, the input material variability is interpreted with both approaches discussed in this section. In the Karhunen–Loeve series expansion, the spatial variability is assumed only in vertical direction.
4.2 The Latin hypercube sampling
The selection of the samples of a random variable can be used by crude Monte Carlo simulation, using pseudorandom algorithms to take a large vector of random values, or variance reduction methods can be implemented such as the importance sampling [11] and the Latin Hypercube Sampling (LHS) [3, 53]. The implementation of LHS method saves a substantial computational effort in order to estimate mean value and coefficient of variation.
In the Latin Hypercube approach let X be a random vector (\(x_1\), \(x_2\),...\(x_n\)). The n dimensional LHS method is stated as follows: For each random variable \(x_i\), the interval [0,1] of the cumulative distribution function (CDF) into N equal subintervals is intersected. Then from each subinterval a random number is chosen and through the inverse CDF a sample of \(x_i\) is obtained. Once samples for all subintervals and all the random variables are acquired then the \(x_i\) vectors are randomly permuted and create the vector realization X. With this procedure, it is assured that to each possible row and each possible column in the nXn euclidean space exactly one sample is taken. The visualization of this feature is depicted in Fig. 2.
The advantage of this procedure is that a reduced amount of values in comparison with crude Monte Carlo approach are needed to integrate the PDF of the input and subsequently to estimate the variability of the output. Also, the subintervals in each dimension may not be equal thus taking into account possible asymmetries of the PDF of the input. Finally, this approach can be applied also to cross-correlated variables, i.e., when the correlation matrix is not diagonal thus it has a general use to the uncertainty quantification scientifical field. In the present work the variables are not cross-correlated and the LHS sampling is used to obtain the variables in the three dimensional space of the compressibility factor k, the critical state line inclination c and the permeability k. The above assumption is reasonable because from experimental data available in the literature the cross-correlation of \(\kappa\) and k can be set to zero due to the coefficient of variation of permeability can be found in some cases much higher than that of volume compressibility, while the correlation between \(\kappa\) and c can be neglected due to lack of experimental data that provide a correlation relation although indirect empirical expressions for the aforementioned variables have been conducted ( [21, 48]).
The material random variables expressed by the PDF \(g_1\) of Sect. 4.1 and the random fields realizations that occur from eq. (9) influence the finite element system of equations (1). The corresponding matrices \({\mathbf {C}}\), \({\mathbf {K}}\), and \({\mathbf {F}}\) are changing due to the randomness of the compressibility factor \(\kappa\), the critical state line inclination c and the permeability k. The selection of the samples follow the importance sampling method of Latin Hypercube Sampling Method.
5 An algorithm for the determination of failure load in ramp dynamic load function
In this section, an algorithm for the determination of failure load, when the dynamic loading function is the ramp loading function, is presented. Ramp loading time function is when the load-time relation is linear with zero initial load until time T. The aim of the algorithm is to find the failure load at the time of the end of the ramp loading. Assume the load function of the problem to be as depicted in Fig. 3.
The aim of this algorithm is to define the load factor \(\lambda ^*\) which causes failure of the continuum at exactly the time T. An initial guess of \(\lambda _1\) is taken leading to an initial time of failure \(t_1\) is obtained. Then, for each new trial \(\lambda _{n+1}\) if the load factor \(\lambda _n\) causes failure to the continuum is given by equation
Otherwise, if the load factor \(\lambda _n\) is not causing failure, then the maximum no failure factor \(\lambda _{max-no-failure}\) is obtained from all previous trials implemented for the calculation of \(\lambda _{n+1}\), as follows:
It should be noted that in practice this recurrence relation usually converges ”by the failure region” which means that only relation (11) is implemented. The difference between \(\lambda _{n+1}\) and \(\lambda _{n}\) is given by:
In eq. (13) it is obvious that as \(n \rightarrow \infty\), \(t_n \rightarrow T\) consequently, the left side of the equation tends to 0, thus the algorithm converges to the desired load factor \(\lambda ^*\).
In comparison with the typical bisection method, where one should guess initial values of failure, safety values \(\lambda _{1,fail}\) and \(\lambda _{1,no-failure}\) and then calculate the new load factor by the bisection of maximum safety factor and minimum failure factor, this algorithm provides less or equal number of trials for the same initial failure guess and convergence tolerance, as can be seen in Table 1. Moreover, it needs only one initial guess, which in high uncertainty problems is very useful in general for avoiding divergence of the solution. Finally, as proven by the numerical tests presented in Table 1 and in Sect. 6, the difference in failure load and failure displacement between the two algorithms is less than 1%. The absolute percentage difference is computed considering as exact solution the bisection algorithm with initial value of safety stress 100 KPa. The performance in terms of the computational time for the 100 deterministic analyses of a Monte Carlo simulation is also depicted in Table 1. The proposed algorithm can save up to 35% of the time demanded for a Monte Carlo analysis to be performed indicating a notable advantage of the aforementioned recurrence relation.
6 Numerical tests in bearing capacity estimation of cohesive soils with random linear and nonlinear material properties
6.1 Description of the problem
The proposed numerical simulation model, which is defined by the set of equations (1) is implemented to porous problems. The geometry of the problem is portrayed in Fig. 5. A uniform vertical load q is applied to the whole soil domain and the ultimate value \(q_{fail}\) alongside with the ultimate displacement \(u_{fail}\) of point A in Fig. 5 are the monitored output variables. The values for the depth are chosen \(h=20\), 40, 50 m. The finite element mesh is with 8 node hexahedral finite elements with linear shape functions for u and p, which provides high numerical accuracy [34, 59]. The length in X and Y directions are \(l_x=l_y=4h\). The stresses due to geostatic loading are directly imported as initial conditions with the relations \(\sigma _v=\gamma z\), \(\sigma _x=\sigma _y=0,85 \sigma _v\) as portrayed in stress point L of Fig. 1. The duration of the simulation in all cases is 1 day in order to have quasi static conditions and in each load case a time step of dt=0.001 day is implemented. The other deterministic properties of the soil are given in Table 2.
Here, \(\nu _0\) stands for the initial specific volume of the soil and \(\lambda _c\) stand for the inclination of Isotropic Compression Line for the respective virgin normally consolidated clay and is considered proportional to \(\kappa\). The following equations apply to boundary surfaces: \({\mathbf {u}}_{x}(z=h)={\mathbf {u}}_{y}(z=h)={\mathbf {u}}_{z}(z=h)={\mathbf {0}}\), and the lateral boundary surfaces are free of constraints. The input material uncertainty consist of the material variables of the compressibility factor \(\kappa\), the critical state line inclination c and the permeability k.
For the compressibility factor \(\kappa\), the spatial distribution along the depth is considered as linear denoted as \(\kappa _{\mathcal {L}}\), or constant denoted as \(\kappa _{\mathcal {C}}\). For \(\kappa _{\mathcal {L}}\), \(\kappa _{depth=0}=0,008686\) and the ratio \(R=\frac{\kappa _{depth=max}}{\kappa _{depth=0}}\) follows the truncated normal distribution. The mean value of the ratio is \(\mu _{R}=0,469\) and the corresponding CoV\(=\frac{\sigma _{R}}{\mu _{R}}\) is 0,25 , as a consequence \(\kappa _{z=max,mean}=0,004074\) and the CoV increases with depth. These values are adopted in order for the mean stiffness of the soil to correspond to a shear velocity of 200 \(\frac{m}{s}\). In Fig. 4, the linear spatial distribution for \(\kappa\) is depicted where the realizations of the compressibility factor are depicted for the values of R: \(\mu _{R}\), \(\mu _{R}+\sigma _{R}\), \(\mu _{R}-\sigma _{R}\). At this point, it is emphasized that that bulk and shear moduli are assumed proportional, since Poisson ratio is assumed constant, consequently \(\kappa\) is directly associated with the shear velocity. In the case of \(\kappa _{\mathcal {C}}\), the mean value of \(\kappa\) is \(\kappa _{\mu }=0,004074\) and the CoV is 0,25.
For the critical state inclination c the spatial distribution with respect to depth is assumed constant. Two possible cases for the absolute value are considered: the random variable analysis and the deterministic value. In the random variable analysis denoted as \(c_{\mathcal {R}}\), the friction angle \(\phi\) follows the truncated normal distribution PDF \(g_1\) of Sect. 6. The mean value is \(\mu _{\phi }=23^o\) while the standard deviation is \(\sigma _{\phi }=2^o\) and this satisfies the condition that the friction angle \(\phi\) to be within an acceptable closed space for physical clays [24]. Subsequently, values of \(\phi\) are generated and c is computed from \(c=\sqrt{\frac{2}{3}} \frac{6sin(\phi )}{3-sin(\phi )}\). The deterministic value denoted as \(c_{\mathcal {D}}\) is c=0,7336 for friction angle \(\mu _{\phi }=23^o\).
Two types of analyses are examined. The porous analyses, where the pore pressure of the soil domain is taken into consideration and the solid analyses, where the water flow is neglected. The solid analyses simulated, indicated with (\({\mathbf {S}}\)) are depicted in Table 3, associating linear (L) or constant (C) distribution for \(\kappa\) and deterministic (D) or random variable (R) cases for c. The porous analyses simulated are portrayed in Table 4, combining linear (L), constant (C), and random field (RF) distribution for \(\kappa\), deterministic (D), random variable case (R) and random field (RF) distribution for c and random field distribution (RF) for k.
In the random field processes, the mean values are assumed: \(\kappa _{mean}=0.008686\), \(c_{mean}\)=0,7336 and \(k_{mean}=10^{-8} \frac{m^3 s}{Mgr}\) [50, 57, 67]. The standard deviations are adopted: \(\sigma _{\kappa }=0.25\kappa _{mean}\), \(\sigma _{\phi }=2^o\) and \(\sigma _k=0.25 k_{mean}\). The autocorrelation function depicted in 10 is implemented in all stochastic processes. The values of correlation length are b=75 (\(k_{\mathcal {RF}{\mathcal {75}}}\)) and b=100 (\(k_{\mathcal {RF}{\mathcal {100}}}\) ). The \(\kappa _{\mathcal {L}}\) and \(\kappa _{\mathcal {C}}\) spatial distributions for \(\kappa\) as well as the random variable distributions for all material variables correspond to a random variable case analysis. In addition, for c a constant deterministic analysis is adopted. The random field (RF) distributions refer to the Karhunen–Loeve series expansion and realizations of the spatial stochastic process are computed through equation \(H_1\) of Sect. 6 with the usage of the autocovariance function 10.
The simulations are static, while the number of Fredholm eigenfunctions taken into account is eight. Failure is defined in the present work as when the first Gaussian Point is having softening behavior (i.e., plastic hardening modulus H0). Each Monte Carlo simulation was applied for 100 samples using the Latin Hypercube Sampling method, which were found sufficient in achieving convergence for the mean value and standard deviation of the monitored displacements as it is depicted in Fig. 6. Here, it should be noted that the cross-correlation in all material variables is set to zero and as a consequence the correlation matrix is diagonal.
6.2 Presentation of the results
6.2.1 Ultimate failure load and displacement
The results are depicted in Tables 5, 6, 7 and in Figs. 7, 8, 9, 10, 11, 12. In these tables the mean values \(q_{fail}\) and \(u_{fail}\) of point A in Fig. 5 are presented together with the standard deviation,the coefficient of variation (CoV), the maximum and minimum value of the Monte Carlo simulation.
When the pore pressure is ignored, larger mean failure displacements and smaller CoV are obtained when \(\kappa _{\mathcal {L}}\) is assumed compared to \(\kappa _{\mathcal {C}}\), while larger mean failure load and corresponding variability are obtained when \(\kappa _{\mathcal {C}}\) is assumed as can be observed from Table 5 for all depths considered. In the numerical tests performed, the largest CoV of the output for \(q_{fail}\) is almost half the input variability while the corresponding CoV of the output for \(u_{fail}\) is practically the same as the input variability of 0,25 (see \({\mathbf {S}}\)-\(\kappa _{\mathcal {C}}\)-\(c_{\mathcal {R}}\) for \(h=20\) m). The mean value of \(u_{fail}\) in \({\mathbf {S}}\)-\(\kappa _{\mathcal {L}}\)-\(c_{\mathcal {D}}\) for \(h=50\) m is 48 % larger than the corresponding value for \(h=20\) m while the mean value for \(q_{fail}\) in \({\mathbf {S}}\)-\(\kappa _{\mathcal {C}}\)-\(c_{\mathcal {D}}\) for \(h=20\) m is 20 % larger than the corresponding value for \(h=50\) m. Consequently, the critical spatial distribution of \(\kappa\) for the CoV of the output of both \(q_{fail}\) and \(u_{fail}\) is the \(\kappa _{\mathcal {C}}\) case. The PDF’s of solid analyses are depicted in Figs. 7 and 8 . This behavior in failure displacements can be explained by the fact that in the \(\kappa _{\mathcal {L}}\) case the upper layers of the soil, which are the most compressible, have low CoV of \(\kappa\) leading to low CoV for both displacements and strains, while in constant distribution along the depth the soil domain tends to be more homogeneous and stiff thus more Gauss Points have larger stiffness leading to larger failure loads.
The CoV of the failure load in porous analyses is slightly influenced by the change in the considered depth, unlike the output variability of the failure displacements, as is portrayed in Table 6. The largest output CoV of \(q_{fail}\), which is found at \(\mathbf {P-}\) \(\kappa _{\mathcal {C}}\)-\(c_{\mathcal {R}}\)-\(k_{\mathcal {RF}{\mathcal {75}}}\) for \(h=20\) m, is 47 % lower than the input variability, while the corresponding maximum CoV of \(u_{fail}\) is 26 % greater than the input variability of 0,25 and is located at \(\mathbf {P-}\) \(\kappa _{\mathcal {C}}\)-\(c_{\mathcal {R}}\)-\(k_{\mathcal {RF}{\mathcal {100}}}\) for \(h=20\) m. Thus, when considering the pore pressure in the soil domain a reduction of variability of the failure load occurs and in failure displacements when there is constant spatial distribution for \(\kappa\) a significant variability increase take place as is portrayed in Figs. 9 and 10 . This is in accordance with the results obtained in [17]. This way of behaving is attributed to the fact that when constant spatial variability of a material variable is assumed all Gauss Points are provided with larger uncertainty leading to greater CoV for strains and displacements. In addition, taking into account that the bulk modulus in porous problems as seen from equation 7, are generally smaller than the respective solid problems, it is concluded that smaller mean values of failure load are expected and smaller variability due to tensile failure of the first Gauss Point as it will be proven numerically in Sect. 6.2.2. and since the tensile stress is set to 0 there are limited acceptable positions for BSE leading to the aforementioned results.
The porous analysis with random field representations for all material variables is subsequently performed as a more general case, since it takes into consideration the spatial randomness of the material properties of the soil. In the case of porous random field analyses, the largest CoV of the output for \(q_{fail}\) is 3,5 times the input uncertainty, while for \(u_{fail}\) the largest output variability is 2,2 times the CoV of the input. The mean values of failure load in porous random field analyses are significantly smaller than the corresponding porous random field analyses with deterministic shape functions for \(\kappa\) and c, while the mean values for failure displacement are notably greater when the linear distribution over depth is assumed for \(\kappa\) as is depicted in Table 7 in \(\mathbf {P-}\) \(\kappa _{\mathcal {RF}}\)-\(c_{\mathcal {RF}}\)-\(k_{\mathcal {RF}{\mathcal {100}}}\) of \(h=50\) m and \(\mathbf {P-}\) \(\kappa _{\mathcal {L}}\)-\(c_{\mathcal {D}}\)-\(k_{\mathcal {RF}{\mathcal {100}}}\) of the same depth. In porous random field analyses the increase of correlation length decreases the CoV of the output in both \(q_{fail}\) and \(u_{fail}\) in depth 20 and 40 meters and reverses in depth 50 m.The opposite occurrence is present for the mean failure load. These behavior can be explained by the fact that the types of failure of Gauss points, when the proposed material yield model is applied, are two: the failure from the “wet” side which in Fig. 1 is the stress points laying in the right side of the intersection between BSE and the critical state line and the “dry” side which is the stress points laying in the left side of the same intersection. Considering the large uncertainty of c and triaxial loading of the Gauss points the failure stress states may be significantly different leading to larger uncertainty of the output. Consequently, the critical spatial distribution for smaller mean failure load, which is the unfavorable situation in the present work, is the Karhunen–Loeve random field representation for all three input variables, which in most cases it also provides the greater output uncertainty.The PDFs of \(\mathbf {P-}\) \(\kappa _{\mathcal {RF}}\)-\(c_{\mathcal {RF}}\)-\(k_{\mathcal {RF}}\) analyses are portrayed in Figs. 11 and 12 .
The results acquired by the analyses examined depict the effect of the input uncertainty of each material variable in porous failure problems. The compressibility factor \(\kappa\), influences the mean value and the variability of both monitored variables. This influence in CoV of the output is more notable when the distribution of \(\kappa\) with respect to depth is constant. This holds in both porous and non-porous problems. This behavior can be interpreted by the fact that \(\kappa\) is directly related with the bulk modulus, consequently has a significant influence on the strains, the displacements and the stiffness of the soil domain leading to an influence in the ultimate load.
The permeability k influences to a lesser extent failure load and failure displacement. The spatial variability of k has a slight sensitivity of the CoV of the output in most Monte Carlo simulations for \(\kappa\) and c with the exception for porous analyses with stochastic processes for all material variables. For the same porous consolidation problem with the same load and deterministic parameters, the output displacement and stress field are not influenced by the permeability since the pore pressures are fully dissipated, leading to the behavior observed by the numerical results.
In conclusion, the variability of critical state line inclination c of the material model appears to have a significant effect on the output variability when the stochastic process calculated from the Karhunen–Loeve series expansion is adopted, while greater mean values are obtained when the \(c_R\) case is implemented. If a deterministic value for c is chosen the output variability is negligible leading to a conclusion that the influence of this material variable is the most important in comparison with \(\kappa\) and k. This is explained as c is directly associated with the failure state of the Gauss Point, leading to influence the limit state of the soil mass.
In order to prove that the output displacement and failure load follows the truncated normal or the lognormal distribution, the histograms of three randomly selected Monte Carlo simulations for an output monitored variable with the corresponding distribution fitting are presented in Fig. 13. As proven by the diagrams, the empirical PDF estimated by the histograms can be modeled by the truncated normal PDF \(g_1\) described in Sect. 6 or the lognormal PDF, respectively. The limit values indicated in Tables 5, 6, 7 are accurate values to set the truncation of the PDF. In addition, for proving numerically the Gaussian nature of the output random vectors, the Kolmogorov–Smirnov test is used [16, 25, 54]. In all output probability density functions, the null hypothesis \(H_0\), at the 5\(\%\) significance level is accepted, as is portrayed in Table 8, where for the Monte Carlo simulation analyses portrayed in Fig. 13 the aforementioned test is implemented. The maximum absolute deviation of the theoretical and empirical cumulative distribution function (CDF) of the Monte Carlo simulations depicted in Fig. 13 are presented and compared to the critical difference for accepting \(H_0\). In all analyses, the critical value is greater than the maximum absolute difference of the CDFs consequently, the output random vectors follow a Gaussian type distribution. This comes in agreement with the previous research literature [17].
6.2.2 Ultimate stresses-strains and failure mechanism
The results are portrayed in Tables 9, 10, 11. In these tables the mean value, the CoV and the minimum values for the volumetric stresses denoted with \(p_h\), the deviatoric stresses denoted with \(q_s\), which is the Von Mises stress explained in Sect. 3.1. In addition, the mean value with the CoV of the volumetric strains denoted with \(\epsilon _{vol}\), deviatoric strains denoted with e are also depicted in Tables9, 10, 11. Moreover, the probability of first Gauss Point failure for all simulations is presented. Here, it should be noted that due to symmetry of the displacement field only the point with smaller coordinates is written. It is also pointed out that the volumetric strain at failure is tensile.
In solid analyses, larger CoV of the output and smaller minimum value of failure stress is obtained when \(c_R\) case is adopted as it can be seen from Table 9. In most cases, the aforementioned distribution also provides greater mean values. The same spatial distribution assumption for critical state line inclination provides larger variability of the strains at failure, which on average are on the proximity of 3-4 ‰. In general, from the results of the Monte Carlo simulations it can be concluded that the percentage plastic deviatoric strains are larger leading to the conclusion that the distortional failure is critical. Finally, in most cases the Gauss Point (2,11 , 2,11 , 12,11) is the most probable failure point with the exception of the \(\kappa _c\) distribution at \(h=40\) m.
In porous analyses with deterministic shape functions for \(\kappa\) and c larger stress variability is obtained at \(h=20\) m with deterministic value for c and at greater depths with \(c_R\) case. The change of the correlation length influences notably the output variability at \(h=20\) m as is portrayed in Table 10. The maximum output variability in failure strains occurs at \(h=50\) m and correlation length for k \(b=75\) m. From the Monte Carlo simulations implemented, it can be seen that when the \(\kappa _c\)-\(c_R\) distributions are assumed the volumetric failure is critical, while when \(\kappa _c\)-\(c_R\) distributions are adopted the deviatoric failure occurs. In conclusion, the most probable failure point is the (2,11 , 2,11 , 12,11) for \(h=20\) m, the (2,11 , 2,11 , 32,11) for \(h=40\) m and the (2,11 , 2,11 , 2,11) for \(h=50\) m and linear distribution for \(\kappa\), while (2,11 , 2,11 , 42,11) is the unfavorable Gauss Point for \(h=50\) m and constant distribution along the depth for the compressibility factor.
In porous analyses with random field representation for \(\kappa\), c, k larger output CoV in stresses is obtained in general when \(b=75\) m except from the deviatoric stresses in \(h=40\) and 50 m as depicted in Table 11. The volumetric strains and the deviatoric strains have in all depths critical correlation length \(b_{critical,vol}\)=75 m and \(b_{critical,dev}\)=100 m, respectively. On average, considering the numerical results obtained by the aforementioned analyses, at \(h=20\) and 50 m the volumetric failure is critical, while at depth 40 m the distortional failure occurs. Finally, the most probable failure points in 20 m are the (2,11 , 2,11 , 2,11) for \(b=75\) m and (2,11 , 2,11 , 12,11) for \(b=100\) m, while in depth 50 m is the (2,11 , 2,11 , 32,11) for both correlation lengths. If the soil depth is 40 meters, there are many points with equally high probability of failure thus leading to a more uncertain position of the failure mechanism.
7 Conclusions
In this paper, the uncertainty quantification of the failure of clayey soils taking into account the pore pressure–soil interaction with the implementation of the stochastic finite element method is presented. The aim of this work is to introduce accurate quantitative results on the failure load of 3D clayey soil domain in relation with the input variability of soil material variables. The proposed numerical simulation model holds for every possible assumptions for the geometry, the loading properties and the material distribution of the soil domain. In this context, a detailed finite element simulation, alongside with a complicated material constitutive model is used. It is evident that similar computational study can be conducted for shallow foundations and it is a future research work that will be performed by the authors.
The numerical results obtained indicate that the output failure load and displacement follow a Gaussian random distribution despite the material nonlinearity, which in failure phenomena is excessive. The randomness of material poroelasticity plays an important role in the output CoV for failure load and failure displacement, especially when it has a Karhunen–Loeve random field representation along the depth of the soil domain. The same thing applies for the critical state line inclination variable c. When the random Karhunen–Loeve field is assumed for all stochastic material variables in discussion the CoV of the output exceeds the corresponding variability of the input in half of the Monte Carlo simulations. This variability increase varies from 30% to 3,5 times larger. The porous variable of permeability influences to a lesser extent the uncertainty of the output material variables. In porous media problems the failure load is smaller compared to the corresponding solid problems.
The random field representations for \(\kappa\), c and k provide the maximum variability for failure stresses and strains. In porous random field analyses in depths \(h=20\), 50 m the volumetric failure occurs, while when the soil domain is 40 m deep the distortional failure is critical. When deterministic shape functions for \(\kappa\) and c are implemented, the \(\kappa _c\)-\(c_R\) distributions provide with larger probability volumetric failure of Gauss Points, while \(\kappa _L\)-\(c_D\) distributions are adopted in order to obtain with larger probability deviatoric failure. Finally, in most cases (2,11 , 2,11 , 12,11) is the critical Gauss Point of failure thus it can be the starting point of the failure mechanism through the Meyerhoff spline.
Availability of data and material
Upon communication with the corresponding author.
References
Ali A, Lyamin A, Huang J, Li J, Cassidy M, Sloan S (2017) Probabilistic stability assessment using adaptive limit analysis and random fields. Acta Geotech 12(4):937–948. https://doi.org/10.1007/s11440-016-0505-1
Ang AS, Tang W (1975) Probability concepts in engineering planning and design, vol 1. Wiley and sons, New Jersey
Olsson A, Sandberg G, Dahlblom O (2003) On latin hypercube sampling for structural reliability analysis. Struct Saf 25(1):47–68. https://doi.org/10.1016/S0167-4730(02)00039-5
Ashraf A, Soubra AH (2012) Probabilistic analysis of strip footings resting on a spatially random soil using subset simulation approach. Georisk Assess Manage Risk Eng Syst Geohazards 6(3):188–201. https://doi.org/10.1080/17499518.2012.678775
Assimaki D, Pecker A, Popescu R, Prevost J (2003) Effects of spatial variability of soil properties on surface ground motion. J Earthquake Eng 7:1–44. https://doi.org/10.1080/13632460309350472
Baecher G, Christian J (2003) Reliability and statistics in geotechnical engineering. pp 177–203, Wiley and sons, New Jersey
Barr DR, Sherrill ET (1999) Mean and variance of truncated normal distributions. Am Stat 53(4):357–361. https://doi.org/10.1080/00031305.1999.10474490
Blaheta R, Beres M, Domesova S (2016) A study of stochastic fem method for porous media flow problem. In: Proceedings of applied mathematics in engineering and reliability, Bris ISBN 978-1-138-02928-6
Borja R (1991) Cam-clay plasticity, part 2: Implicit integration of constitutive equation based on a nonlinear elastic stress predictor. Comput Methods Appl Mech Eng 88(2):225–240. https://doi.org/10.1016/0045-7825(91)90256-6
Borja R, Lee S (1990) Cam-clay plasticity, part 1: Implicit integration of elasto-plastic constitutive relations. Comput Methods Appl Mech Eng 78(1):49–72. https://doi.org/10.1016/0045-7825(90)90152-c
Bouhari A (2004) Adaptative monte carlo method, a variance reduction technique. Monte Carlo Methods Appl 10(1):1–24. https://doi.org/10.1515/156939604323091180
Brantson ET, Ju B, Wu D, Gyan PS (2018) Stochastic porous media modeling and high-resolution schemes for numerical simulation of subsurface immiscible fluid flow transport. Acta Geophys 66(3):243–266. https://doi.org/10.1007/s11600-018-0132-3
Chwała M (2019) Undrained bearing capacity of spatially random soil for rectangular footings. Soils Found 59:1508–1521. https://doi.org/10.1016/j.sandf.2019.07.005
Chwała M, Puła W (2020) Evaluation of shallow foundation bearing capacity in the case of a two-layered soil and spatial variability in soil strength parameters. PLoS One 15(4):e0231992. https://doi.org/10.1371/journal.pone.0231992
Das BM (2011) Principles of foundation engineering. Global Engineering- Christopher M Shortt
Dimitrova D, Kaishev V, Tan S (2019) Computing the Kolmogorov–Smirnov distribution when the underlying cdf is purely discrete, mixed or continuous. J Stat Softw 95(1):1–42
Fenton G, Griffiths D (2003) Bearing capacity prediction of spatially random c-\(\phi\) soils. Can Geotech J 40(1):54–65. https://doi.org/10.1016/j.probengmech.2005.06.003
Fu D, Zhang Y, Yan Y (2020) Bearing capacity of a side-rounded suction caisson foundation under general loading in clay. Comput Geotech 123:103543. https://doi.org/10.1016/j.compgeo.2020.103543
Ghalehjough BK, Akbulut S, Çelik S (2018) Effect of particle roundness and morphology on the shear failure mechanism of granular soil under strip footing. Acta Geotech Slov 15(1):43–53. https://doi.org/10.18690/actageotechslov.15.1.43-53.2018
Ghanem R, Spanos D (1991) Stochastic finite elements: a spectral approach. vol 1, pp 1–214. Springer, Berlin https://doi.org/10.1007/978-1-4612-3094-6
Huang J, Griffiths DV, Fenton GA (2010) Probabilistic analysis of coupled soil consolidation. J Geotech Geoenviron Eng ASCE 136(3):417–430. https://doi.org/10.1061/(ASCE)GT.1943-5606.0000238
Kalos A (2014) Investigation of the nonlinear time-dependent soil behavior. PhD Dissertation NTUA 1:193–236
Karhunen K (1947) Uber lineare methoden in der wahrscheinlichkeitsrechnung. Ann Acad Sci Fenn Ser A 37:1–79
Kavvadas M, Amorosi A (2000) A constitutive model for structured soils. Geotechnique 50(3):263–273. https://doi.org/10.1680/geot.2000.50.3.263
Kolmogorov A (1933) Sulla determinazione empirica di una legge di distribuzione. G Ist Ital Attuari 4:83–91
Kötter F (1903) Die bestimmung des drucks angekrummten, eineaufgabe aus der lehre vom erddruck. Sitzungsberichte derAkademie der Wissenschaften, pp 229–233. Berlin
Li C, Kiureghian AD (1993) Optimal discretization of random fields. J Eng Mech 119(6):1136–1154. https://doi.org/10.1061/(ASCE)0733-9399(1993)119:6(1136)
Li DQ, Qi XH, Cao ZJ, Tang XS, Zhou W, Phoon KK, Zhou CB (2015) Reliability analysis of strip footing considering spatially variable undrained shear strength that linearly increases with depth. Soils Found 55(3):866–880. https://doi.org/10.1016/j.sandf.2015.06.017
Li S, Yu J, Huang M, Leung G (2021) Upper bound analysis of rectangular surface footings on clay with linearly increasing strength. Comput Geotech 129:103896. https://doi.org/10.1016/j.compgeo.2020.103896
Liu W, Sun Q, Miao H, Li J (2015) Nonlinear stochastic seismic analysis of buried pipeline systems. Soil Dyn Earthq Eng 74:69–78. https://doi.org/10.1016/j.soildyn.2015.03.017
Martin C (2005) Exact bearing capacity calculations using the method of characteristics. Proceedings of the 11th International Conference IACMAG Graz, Austria
Matthies HG, Brenner CE, Butcher G, Soares CG (1997) Uncertainties in probabilistic numerical analysis of structures and solids- stochastic finite elements. Struct Saf 19(3):283–336. https://doi.org/10.1016/s0167-4730(97)00013-1
Meftah F, Dal-Pont S, Schrefler BA (2012) A three-dimensional staggered finite element approach for random parametric modeling of thermo-hygral coupled phenomena in porous media. Int J Numer Anal Meth Geomech 36:574–596. https://doi.org/10.1002/nag.1017
Melenk JM, Babuska I (1996) The partition of unity finite element method: basic theory and applications. Comput Methods Appl Mech Eng 139(1–4):289–314. https://doi.org/10.1016/s0045-7825(96)01087-0
Meyerhoff GG (1951) The ultimate bearing capacity of foundations. Geotechnique 2:301–332
Michalowski RL (1997) An estimate of the influence of soil weight on bearing capacity using limit analysis. Soils Found 37(4):57–64. https://doi.org/10.3208/sandf.37.4_57
Michalowski RL (2001) Upper-bound load estimates on square and rectangular footings. Geotechnique 51(9):787–798. https://doi.org/10.1680/geot.2001.51.9.787?journalCode=jgeot
Naderi E, Asakereh A, Dehghani M (2018) Bearing capacity of strip footing on clay slope reinforced with stone columns. Arab J Sci Eng 43:5559–5572. https://doi.org/10.1007/s13369-018-3231-1
Papadopoulos V, Giovanis D (2018) Stochastic finite element methods: an introduction. vol 1, pp 30–35 Springer, Berlin https://doi.org/10.1007/978-3-319-64528-5
Papadopoulou K, Gazetas G (2020) Shape effects on bearing capacity of footings on two-layered clay. Geotech Geol Eng 38:1347–1370. https://doi.org/10.1007/s10706-019-01095-6
Papadrakakis M, Papadopoulos V (1996) Robust and efficient methods for the stochastic finite element analysis using monte carlo simulation. Comput Methods Appl Mech Eng 134:325–340. https://doi.org/10.1016/0045-7825(95)00978-7
Peng X, Zhang L, Jeng D, Chenc L, Liao C, Yang H (2017) Effects of cross-correlated multiple spatially random soil properties on wave-induced oscillatory seabed response. Appl Ocean Res 62:57–69. https://doi.org/10.1016/j.apor.2016.11.004
Popescu R, Deodatis G, Prevost J (1998) Simulation of homogeneous nongaussian stochastic vector fields. Probab Eng Mech 13:1–13. https://doi.org/10.1016/s0266-8920(97)00001-5
Popescu R, Deodatis G, Nobahar A (2005) Effects of random heterogeneity of soil properties on bearing capacity. Probab Eng Mech 20:324–341. https://doi.org/10.1016/j.probengmech.2005.06.003
Prandtl L (1920) Uber die harte plastischer korper. nachrichten von derkonilichen gesellschaft der wissenschaften zu gottingen. Mathematisch-Physikalische Klasse au dem Jahre, pp 74–85. Berlin
Pryse S, Adhikari S (2017) Stochastic finite element response analysis using random eigenfunction expansion. Comput Struct 192:1–15. https://doi.org/10.1016/j.compstruc.2017.06.014
Rao P, Liu Y, Cui J (2015) Bearing capacity of strip footings on two-layered clay under combined loading. Comput Geotech 69:210–218. https://doi.org/10.1016/j.compgeo.2015.05.018
Reid D (2015) Estimating slope of critical state line from cone penetration test - an update. Can Geotech J 52(1):46–57. https://doi.org/10.1139/cgj-2014-0068
Robert CP (1995) Simulation of truncated normal variables. Stat Comput 5(2):121–125. https://doi.org/10.1007/BF00143942
Lewis RW, Schrefler BA (1988) The finite element method in the deformation and consolidation of porous media. vol 1, pp 1–508. Wiley and Sons, New Jersey https://doi.org/10.1137/1031039
Savvides A, Papadrakakis M (2020) A probabilistic assessment for porous consolidation of clays. SN Appl Sci 2:2115. https://doi.org/10.1007/s42452-020-03894-6
Sett K, Jeremic B (2007) Probabilistic elasto-plasticity: solution and verification in 1d. Acta Geotech 2(3):211–220. https://doi.org/10.1007/s11440-007-0037-9
Simoes J, Neves L, Antao A, Guerra N (2020) Reliability assessment of shallow foundations on undrained soils considering soil spatial variability. Comput Geotech 119:103369. https://doi.org/10.1016/j.compgeo.2019.103369
Smirnov N (1948) Table for estimating the goodness of fit of empirical distributions. Ann Math Stat 19(2):279–281. https://doi.org/10.1214/aoms/1177730256
Stavroulakis G, Giovanis D, Papadopoulos V, Papadrakakis M (2014a) A gpu domain decomposition solution for spectral stochastic finite element method. Comput Methods Appl Mech Eng 327:392–410. https://doi.org/10.1016/j.cma.2017.08.042
Stavroulakis G, Giovanis D, Papadopoulos V, Papadrakakis M (2014b) A new perspective on the solution of uncertainty quantification and reliability analysis of large-scale problems. Comput Methods Appl Mech Eng 276:627–658. https://doi.org/10.1016/j.cma.2014.03.009
Stickle MM, Yague A, Pastor M (2016) Free finite element approach for saturated porous media: consolidation. Math Prob Eng. https://doi.org/10.1155/2016/4256079
Sultana P, Dey AK (2019) Estimation of ultimate bearing capacity of footings on soft clay from plate load test data considering variability. Indian Geotech J 49:170–183. https://doi.org/10.1007/s40098-018-0311-9
Szabo B, Babuska I (2011) Intoduction to finite element analysis: formulation, verification and validation. Wiley Ser Comput Mech 1:1–382. https://doi.org/10.1002/9781119993834
Terzaghi KV (1966) Theoretical soil mechanics. Wiley and Sons, New Jersey
Vesic AS (1963) Bearing capacity of deep foundation in sand. Highw Res Rec Natl Acad Sci 39:112–153
Vesic AS (1973) Analysis of ultimate loads of shallow foundations. J Soil Mech Found Div ASCE 99(1):45–73
Vrakas A (2018) On the computational applicability of the modified cam-clay model on the “dry” side. Comput Geotech 94:214–230. https://doi.org/10.1016/j.compgeo.2017.09.013
Liu W, Belytschko T, Mani A (1986) Random fields finite element. Int J Numer Methods Eng 23:1831–1845. https://doi.org/10.1002/nme.1620231004
Yue Q, Yao J, Alfredo H, Spanos PD (2018) Efficient random field modeling of soil deposits properties. Soil Dyn Earthq Eng 108:1–12. https://doi.org/10.1016/j.soildyn.2018.01.036
Zhou H, Zheng G, Yin X, Jia R, Yang X (2018) The bearing capacity and failure mechanism of a vertically loaded strip footing placed on the top of slopes. Comput Geotech 94:12–21. https://doi.org/10.1016/j.compgeo.2017.08.009
Zienkiewicz OC, Chan AHC, Pastor M, Schrefler BA, Shiomi T (1999) Computational geomechanics with special reference to earthquake engineering. vol 1, pp 17–49. Wiley, Chichester
Acknowledgements
This work has been supported by the European Research Council Advanced Grant MASTER Mastering the computational challenges in numerical modeling and optimum design of CNT reinforced composites (ERC-2011-ADG-20110209). The author also acknowledges support from the Bodossaki Foundation. The authors state that there is no conflict of interest.
Open Access
This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Funding
As stated in the Acknowledgments section.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflicts of interest
As stated in the Acknowledgments section.
Code availability
Open source code MSolve in programming Language C# . Information to the link http://mgroup.ntua.gr/
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Savvides, A.A., Papadrakakis, M. A computational study on the uncertainty quantification of failure of clays with a modified Cam-Clay yield criterion. SN Appl. Sci. 3, 659 (2021). https://doi.org/10.1007/s42452-021-04631-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s42452-021-04631-3