Next Article in Journal
Probiotic Characteristics of Ligilactobacillus salivarius AS22 Isolated from Sheep Dung and Its Application in Corn-Fox Tail Millet Silage
Next Article in Special Issue
Experimental Evaluation of the Bending Behavior of a Drilled Shaft with Partial Casing under Lateral Loads
Previous Article in Journal
Biological Diabrotica Management and Monitoring of Metarhizium Diversity in Austrian Maize Fields Following Mass Application of the Entomopathogen Metarhizium brunneum
Previous Article in Special Issue
Study on Soil Displacement Fields around the Expanded Body of Drill-Expanded Concrete Piles Based on DIC Technique
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bearing Capacity of Single Pile-Friction Wheel Composite Foundation on Sand-over-Clay Deposit under V-H-M Combined Loadings

1
College of Civil Engineering, Hunan University, Changsha 410082, China
2
Key Laboratory of Building Safety and Energy Efficience of Ministry of Education, College of Civil Engineering, Hunan University, Changsha 410082, China
*
Author to whom correspondence should be addressed.
Submission received: 17 September 2021 / Revised: 3 October 2021 / Accepted: 7 October 2021 / Published: 12 October 2021
(This article belongs to the Special Issue Recent Progress on Advanced Foundation Engineering)

Abstract

:
This paper presents numerical modelling to investigate the bearing capacities and failure mechanisms of single pile-friction wheel composite foundation in sand-overlying-clay soil conditions under combined V-H-M (vertical-horizontal-moment) loadings. A series of detailed numerical models, with validations of centrifuge testing results, are generated to explore the potential factors influencing the bearing capacity of this composite system. Intensive parametric study is then performed to quantify the influences of the foundation geometry, soil properties, sand layer thickness, pre-vertical loading and lateral loading height on the failure envelopes in the V-H-M domain. Last but not least, an empirical design procedure is proposed based on a parametric study to predict the bearing capacity of this composite foundation under various loading conditions, which can provide guidance for its design and application.

1. Introduction

1.1. Concept

Large diameter monopile is often applied to support offshore structures such as wind turbines (OWTs) [1,2,3,4]. This is because of its sufficient theoretical basis for analysis design and mature manufacturing technology. However, the bearing capacities provided by monopiles are likely to not be suitable for the next generation of larger wind turbines in more severe environments [5,6,7,8,9]. For example, for OWTs with a full capacity of more than 10 MW, it would be costly and bring more challenges for the construction and installation procedure by merely increasing the diameter and embedment length of the monopile. Thus, a new type of single pile-friction wheel composite foundation has therefore been proposed, which combines the monopile and circular footing foundations. Compared to the traditional monopile foundation, this composite foundation provides higher bearing capacity and stiffness with less construction time and expense [10,11,12], as is shown in Figure 1. This composite foundation is initiated from the concept of pile caps and embedded retaining walls with stabilizing platforms, in which the addition of the circular wheel increases the shear stress and restoring moment against the lateral deflection [13,14,15]. The combined wheel can also be applied as an effective reinforcement for existing monopiles to provide additional bearing capacities and protections of the seabed from circum-pile scouring [15,16,17].
Sand-overlying-clay soil deposit is a complex, but commonly encountered, soil profile in petroleum regions and offshore wind farms, for example, in the Yellow Sea of Korea, North Sea, Gulf of Mexico, South China Sea, offshore India and Thailand [18,19]. Existing understandings only cover the failure mechanisms of conventional types of foundations in such sand-overlying-clay soil profiles [18,19,20,21,22]. For example, Zheng et al. (2018) [23] found a ‘punch-through’ failure pattern for strip footings in the sand-over-clay layered deposits, where the upper sand block is pushed into the underlying clay due to the vertical compression under V-H-M combined loadings. However, the bearing capacities and failure mechanisms of this composite foundation are still not well investigated in such soil conditions. It is also worth noticing that, in order to make full use of the bearing capacity provided by the friction wheel, the composite foundation is suggested to be installed in the regions where the strength of the seabed topsoil is great enough, such as pure sand, sand-over-clay and stiff clay-over-sand soil conditions [24,25,26]. Hence, a systematic study is badly needed for the application of composite foundations in sand-over-clay layered deposits.

1.2. Previous Work

A number of studies have been carried out on the bearing capacity of the composite foundation through centrifuge model tests, detailed three-dimensional finite element analyses and theoretical derivations. Most of these studies are only limited to the behaviors of these foundations in pure sand or clay [10,11,24,25,26,27,28,29,30,31,32,33,34,35,36].
For the composite system, the interaction effect between the monopile and the wheel has been mostly conducted in two cases: (i) the composite foundation with a wheel diameter (Dw) exceeds the pile embedment depth (L) (Dw/L > 1); and (ii) with a ratio of Dw/L = 0.5~1. When the wheel diameter is larger than the pile embedment depth as in case (i), Wang et al. (2018) [10] and Yang et al. (2019) [11] conducted centrifuge tests, in which it was found that the lateral bearing capacity of the composite system was actually smaller than the summation of the individual capacities of monopile and wheel. To quantify the reduction, Wang et al. (2018) [10] and Li et al. (2020) [35] proposed an add-up method with a reduction factor to estimate the lateral bearing capacity of the composite foundation. On the contrary, with a smaller wheel as in case (ii), a significant positive interaction effect was found by Lehane et al. (2014) [27] from their centrifuge tests, in which the bearing capacity of the composite foundation is greater than the sum of that provided by the individual monopile and wheel. In addition, Anastasopoulos and Theofilou (2016) [28] proposed that the increase of pile length L had a more pronounced effect on the bearing capacity of the composite foundation than the improvement through increasing wheel diameter Dw.
For the failure mechanism of the composite foundation, reported by Anastasopoulos and Theofilou (2016) [28], substantial soil yielding was captured at the edge of the wheel in front of the lateral loading direction. The increase of vertical load (V) on the foundation led to the mobilization of the wheel and the soil mass beneath. The load transfer mechanism was also found to be significantly different between two configurations, that is, coupled and decoupled composite systems, respectively. As suggested by Stone et al. (2018) [36], for a coupled system, the wheel rotated as a rigid body with the center on the pile axis; however, it was found that the uniform settlement of the wheel triggered by the vertical load played a dominant role in the failure mechanism for the decoupled arrangement. In order to assess the bearing capacity of the composite foundation, an empirical design procedure was proposed by adopting the concept of the soil wedge under the wheel in front of the composite foundation [34]. To investigate the bearing capacities of the composite foundation under V-H-M combined loadings, centrifuge tests have been carried out by El-Marassi (2011) [25] and Yang et al. (2019) [11]. It was concluded that the lateral load bearing capacity and the moment resistant capacity were improved by increasing the vertical load V on the composite foundation, and the maximum lateral load resistance occurred at V/Vult = 0.4~0.5 (Vult is the ultimate vertical loading capacity).
Unlike the previous studies, which mainly focus on the bearing behavior of the composite foundation in pure sand or clay, in this study detailed finite element analyses are carried out to study the bearing capacities and failure mechanisms of the composite foundation in sand-overlying-clay deposits under V-H-M combined loadings. Besides that, an empirical design method is innovatively proposed based on an intensive parametric study to quantify the failure envelopes considering the influences of foundation geometry (Dw/L and L/Dp: Dw is the wheel diameter; Dp is the wheel diameter; L is the pile embedment depth), soil layer distribution condition (Ts/L: Ts is the thickness of sand layer), soil properties (tan(φ)γsDp/sum: φ is the internal friction wheel of sand, γs is the buoyant unit weight of sand, sum is the undrained shear strength of clay at the soil surface) and the lateral loading height (e), which are essential for design management during application.

2. Materials and Methods

Detailed three-dimensional finite element analyses are conducted to examine the failure mechanism and the bearing capacity of the composite foundation in sand-overlying-clay using commercial software ABAQUS [37]. The numerical modeling helps to investigate the soil deformation and stress distribution around the composite foundation, and the failure surface in the limiting state.

2.1. Finite Element Model

As shown in Figure 2, the soil domain is set with the horizontal boundary width (horizontal distance from the wheel edge to the mesh boundary) as 5 Dw and the vertical boundary height (vertical distance from the pile end to the mesh boundary) as 2 L, which are large enough to ensure minimal boundary effects [30,31,32]. Three-dimensional 8-node linear brick elements (C3D8R) with reduced integration are adopted in the models. Relatively fine meshes are adopted around the foundation for a higher computational accuracy, and coarser mesh sizes are employed in the soil away from the foundation to optimizing the computer resources needed. Hinge supports are applied at the bottom of the model and roller conditions are set along the vertical sides of the soil domain. The vertical boundary in the front side of the model is fixed in the horizontal direction (y direction) and set free in the other two directions (x and z directions). In order to simplify the simulation, the analyses involve a “wished in place” pile installation before loading. Beam elements are used to model the tower structure with a linear elastic material property. A lumped mass is adopted to simulate the super-structure assembly, which applies the vertical pre-load by its self-weight [38,39]. The corresponding load can be extracted from the loading point as shown in Figure 2. The resultant lateral load (H) is applied at a certain height (e) on the tower above the mudline, which leads to the combined lateral load (H) and the designed moment (M = H × e) acting on the foundation. The height of the lateral load was determined according to the engineering practice [40]. The lateral load is applied in the form of continuous lateral displacement until the foundation reaches its ultimate bearing state, which is determined by the tangent intersection method [41,42]: the load corresponding to the intersection point of the lateral load-displacement curve is regarded as the ultimate bearing capacity of the foundation.

2.2. Constitutive Model and Soil Properties

In this study, a composite system used to support a wind turbine is considered, which consists of a hollow circular pipe pile with pile diameter Dp = 4 m, wall thickness t = 0.2 m and a solid wheel with the thickness Tw = 2 m. The wheel diameter Dw and pile embedment depth L are determined as ratios of Dw/L and L/Dp varying from 0 to 0.95 and 7 to 13, and the effects of wheel thickness (Tw) and height of wheel collar (Tcollar) are also investigated. The steel pile and wheel are modelled as an elastic-perfectly material with the Young’s modulus of 210 GPa, Poisson’s ratio ν of 0.17, and unit weight of 78.5 kN/m3. This assumption of the foundation has also been widely used in many other researches [11,28,31,32], and good predictions of the bearing performance have been obtained compared to that of the real structure. The sand is assumed to have the buoyant unit weight γs of 9.8 kN/m3, internal friction angle φ of 33°, Poisson’s ratio (υ) of 0.3, Young’s modulus Es of 40 MPa, dilation angle ψ of 3°, and cohesive force c of zero, which are obtained from the triaxial tests. The friction coefficient between the sand and foundation is measured as 0.3 from the lab tests as reference [32]. By adopting a perfect contact with Coulomb’s tangent friction model, interactions between the foundation and the sand are simulated allowing detachment and sliding. The buoyant unit weight of clay γc is taken as 8.5 kN/m3 [32]. A uniform stiffness ratio Ec/su = 500 (where Ec is the Young’s modulus of the clay, su is the undrained shear strength of clay) [43] is adopted throughout the soil domain and the Poisson’s ratio ν is set as 0.49 to model the fully undrained condition. The undrained shear strength of clay is described with a normalized soil strength gradient in a linearly increasing strength profile, as follows:
su = sum + kz,
where sum is the undrained shear strength at soil surface, k is the shear strength gradient with the submerged depth z and is taken as 1.75 kPa/m in the following analyses [4,44]. The shear strength of the nodal joint elements near the clay/pile interface is limited to αsu (α is the interface friction coefficient which is set as 0.3 for the composite foundation/clay interface). This method has been successfully adopted in simulating the interaction between pile and clay in small strain analyses, where good predictions were obtained [45,46]. The gravity is initially applied to the soil to establish geostatic stresses state in soils, and the coefficient of lateral earth pressure at rest is set as K0 = 1 − sinφ.

2.3. Numerical Model Validation

The generated numerical model is validated using the lateral load-displacement curves that are obtained from the centrifuge tests [27]. (Case 1: Dw/L = 0.5, V = 5.3 MN; Case 2: wheel only, V = 5.9 MN; Case 3: wheel only, V = 12.7 MN; Case 4: Dw/L = 0.5, V = 6.4 MN; Case 5: Dw/L = 0.5, V = 13.2 MN, where the other parameters are: Dp = 3.33 m, e/Dp = 9, Ts/L = 1; Group I, Table 1). As shown in Figure 3, reasonably close agreements can be found between the results obtained from finite element method (FEM) and the centrifuge tests. We can therefore be confident that the above developed numerical modeling method could reasonably closely replicate the bearing behavior of the composite foundation in sand-overlying-clay deposit.

3. Effect of Different Influencing Factors

3.1. Effects of Wheel Diameter and Pile Embedment Depth (Dw/L and L/Dp)

To investigate the effects of the wheel diameter and pile embedment depth on the lateral bearing performance, a group of numerical analyses is carried out with seven different wheel diameter ratios for the composite foundation (Dw/L = 0, 0.3, 0.4, 0.5, 0.6, 0.7 and 0.8, L/Dp = 7, 9, 11 and 13; Group I, Table 1). Figure 4a shows the normalized lateral capacities (Hult/(2πR3γs)) of the composite foundation with different wheel diameters. It is apparent that the lateral capacity and the wheel diameter show an almost linear relationship when Dw/L is within the range of 0 to 0.5. The slope of the curve experiences a rapid increase at the turning point of about Dw/L = 0.5, which shows a good agreement with the observation from the previous investigations [10,11,32]. The soil plastic strain distribution at the ultimate condition in terms of plastic strain within and around the composite foundation are also illustrated in Figure 4a. It can be seen that through comparing the composite foundation and the monopile, the main difference of the failure pattern lay on the area at the right side of the foundation. When Dw/L = 0, the plastic strain generates from the bed surface at about 3 Dp away from the right side of the pile, which propagates inwards to the right exterior side of the pile, forming a hill-type failure zone (f2). Besides, the failure zone f3 develops from the mudline to the interface of the sand and clay layer along the left pile body side. The failure zone f4 occurs on both sides of the pile in the clay layer due to the rotational movement of the pile. As for the composite foundation with Dw/L = 0.6, a semi ellipse-shape failure zone (f1) occurs near the right external side of the wheel edge due to the intrusion of the wheel during rotation. Failure zone (f2) at the intersection area of the pile and the wheel shows a relative smaller value compared to that of the monopile. The failure zone (f3) propagates along the left side of the pile, and eventually reaches the interface of the two soil layers. Similar shape of failure zone (f4) can be captured on the composite foundation as that on the monopile. It is evident that the failure pattern of the composite foundation is considerably different from that of the monopile due to the integration of the wheel. As the wheel diameter Dw/L increases from 0.6 to 0.8, the failure zone (f1) enlarges along the loading direction, indicating a greater bearing capacity.
Figure 4b shows the vertical passive earth pressure distribution under the wheel through presenting the normalized passive pressure (Pv/(γsDp)) versus normalized distance to the wheel center (dw/Dp) relationship. It is evident that the overall distribution of the vertical passive earth pressure is similar for the foundation with different wheel diameters, which indicates an increasing trend as the distance develops towards the wheel edge with the maximum value at about 0.75 Dw. It can be explained that the right intersection soil area between the pile and the wheel experiences soil densification due to the compression of the pile and the wheel, which generates a higher soil resistance vertically. The total vertical passive resistance generated under the wheel increases as the wheel diameter increases, which provides a larger bearing area.

3.2. Effects of Pre-Vertical Load and Lateral Loading Height (V/Vult and e/Dp)

Numerical modeling is carried out to examine the effect of the pre-vertical load and the lateral loading height, where the foundation geometry and the soil properties are listed in Group III of Table 1. The relationship of the normalized lateral bearing capacity and the pre-vertical load is plotted in Figure 5. It can be found that, with Dw/L = 0.6, the maximum lateral bearing capacity of the composite foundation occurs when V/Vult = 60% for different lateral loading heights. It can be concluded that the lateral loading height has little influence on the shape of the failure envelopes but only alter asymptotes.
Figure 5 also illustrates the failure mechanism of the composite foundation under the combined V-H-M loadings. The plastic strain area f1 is found to expand downwards as the pre-vertical load increases, which means more soils under the wheel are mobilized to provide the bearing capacity. In addition, with the elevation of the lateral loading height, plastic strain area f1 and f3 increase and show greater values. This behavior can be attributed to the increasing dominance of the rotational movement of the composite foundation as the lateral loading height increases, which results in a lower initial stiffness and lateral bearing capacity of the composite foundation.

3.3. Effects of Sand Layer Thickness and Soil Properties (sum, k and φ)

To study the effect of the soil layer thickness and the soil properties, a group of numerical analyses is conducted with Ts/L = 0.1~0.8, sum = 20 kPa and 40 kPa and k = 1.75 and 2 kPa/m (Group IV, Table 1). The influences of the internal friction angle on sand (φ) and Ts/L are presented in Figure 6, where the normalized lateral bearing capacity is plotted against Ts/L with different φ and sum. It can be found that the lateral bearing capacity of the composite foundation increases nonlinearly as Ts/L increases. The increasing trend becomes more evident with a greater value of φ, which is partly due to the P-Δ effect. From the failure mechanism of the composite foundation depicted in Figure 6a,b, it can be seen that, when Ts/L = 0.1 and φ = 20°, a “punching-through” failure pattern is developed in region (f1), which includes a rigid block near the pile beneath the footing, a passive wedge at lower right of the wheel in the sand layer, and a fan in the lower clay layer. Similar conclusions are also reported by Zheng et al. (2018) [23] in their investigation of the failure mechanism of the strip footing in sand-over-clay soil deposit. In addition, more evident plastic strain areas (f2) and (f3) at the interface of the two soil layers and the plastic strain areas (f4) under the wheel are captured when Ts/L = 0.1 and φ = 20° as compared to that of φ = 40°. When Ts/L increases to 0.8, a significant plastic strain area (f4) can be seen on the lower left region of the wheel when φ = 20°. Generally, increasing sum can significantly enhance the bearing capacity of the composite foundation when Ts/L is relatively small, while soil strength gradient k has limited influence on the bearing capacity, which can be concluded by the Hult/(2πR3γs)-Ts/L curves and the failure mechanisms.
Figure 7 demonstrates the lateral earth pressure distributions along the right side of the pile body under different Ts/L and φ by plotting normalized lateral earth pressure Pp/(2πs) versus depth below the mudline z/Dp. In terms of Ts/L = 0.1, there are no obvious differences between the earth pressure distributions when φ = 20° and 40°, with the same depth of maximum earth pressure zmax = 0.4 L, which is in accordance with the case when the composite foundation is in pure clay [32]. Besides that, no abrupt change is captured at Ts/L = 0.1, and it indicates that the bearing capacity is mainly controlled by the combined wheel rather than the soil resistance provided by the pile. For Ts/L = 0.4, the overall earth pressure distribution is significantly different from that of Ts/L = 0.1, with its zmax elevate to above 0.3 L. Apparently, a sharp decrease of soil resistance can be found at Ts/L = 0.4 as the depth descends. In addition, the depth of the maximum earth pressure raises with larger φ (φ = 20°: zmax/L = 0.28; φ = 40°: zmax/L = 0.19). As Ts/L increases from 0.4 to 0.8, for φ = 20°, the earth pressure along pile barely changes, while it experiences an obvious increase within the range of z/L = 0~0.55 when φ = 40°, which is consistent with the phenomenon shown in Figure 6.

3.4. Effects of Wheel Thickness and Height of Collar (Tw/Dp and Tcollar/Dp)

To explore the effect of wheel thickness Tw and the height of collar Tcollar, a group of analyses was conducted (Group V, Table 1). Figure 8 illustrates the relationship between ultimate lateral bearing capacities and normalized wheel thickness Tw/Dp with different Tcollar/Dp. It shows that Tcollar has a negligible influence on the capacity of composite system. As Tw/Dp changes from 0.5 to 1.25, the bearing capacity increases linearly by about 17%. This is because that additional vertical pressure is applied to the soil underneath the wheel when Tw is greater, and it can be validated by the plastic strain distribution displayed in Figure 8, which is similar to the situation in Figure 5 when vertical load is applied to the composite system. In addition, with the aim of cost-efficiency in the manufacture and installation of the composite foundation, a friction wheel with a greater thickness would not be suitable for the application of this type of foundation. Therefore, based on the above analysis, it can be confirmed that the Tw/Dp and Tcollar/Dp have minimal influences on the bearing capacity.

4. Failure Envelopes in the V-H-M Domain

4.1. Parametric Study

Parametric studies are carried out to quantify the influences of the soil properties (tan(φ)γsDp/sum), sand layer thickness (Ts/L) and foundation geometry (DwDp/L2) on the bearing capacity under combined V-H-M loadings. Figure 9 presents the influence of soil properties on the bearing capacities in the V-H-M domain by plotting failure envelopes with varying tan(φ)γsDp/sum. (DwDp/L2 = 0.062 and 0.08, Ts/L = 0.1 and 0.7, e/Dp = 15). Figure 9a shows that, the shape and value of failure envelopes are both related to the soil properties. From Figure 9b, it can be observed that the peak lateral capacity increases at an accelerating rate with the increase of tan(φ)γsDp/sum. Figure 10 shows the effect of the foundation geometry on the failure envelops in the V-H-M domain. The failure envelopes expand with the increase of DwDp/L2 under a constant upper sand thickness, which means the foundation geometry has little influence on the failure mechanism but affects the magnitude of the bearing capacity of the composite foundation. As can be seen in Figure 11, when Ts/L = 0.1, under a constant DwDp/L2, the lateral bearing capacity decreases continuously when the vertical load increases. However, as Ts/L exceeds 0.1, the peak bearing capacity can be observed in the failure envelopes. The normalized vertical load corresponding to the peak value grows as the thickness of the upper sand layer increases. For instance, it increases from 50 to 180 as Ts/L increases from 0.3 to 0.7 when DwDp/L2 = 0.062. This is because a thin layer of top sand may result in a different failure pattern like the “punch-through” failure under the combined V-H-M capacities.

4.2. Approximate Expression of the Bearing Capacity in the V-H-M Domain

Based on the parametric study, it can be concluded that both the sand layer thickness Ts/L and the soil properties tan(φ)γsDp/sum influence the shape of failure envelopes under the combined V-H-M loadings. To estimate the bearing capacities conveniently, the failure envelopes normalized by the corresponding H0 and Vult are shown in Figure 12 (H0 is the uniaxial horizontal capacity; Vult is the vertical capacity of the foundation), and consider the mobilization of the sand layer thickness and the soil properties. Inheriting the format of the empirical formula by El-Marassi (2011) for composite foundations in pure sand, Equation (2) is proposed to determine the shape of failure envelopes in the combined V-H-M domain. A, B, C, D, and E are coefficients related to Ts/L and tan(φ)γsDp/sum, which are derived in Table 2. (h0 is the dimensionless uniaxial horizontal capacity; v0 is the dimensionless vertical load)
H ult H 0 A × ( V V ult ) B × ( C D × ( V V ult ) E ) = 1
h 0 = a 1 + b 1 × D w 0 + ( c 1 + d 1 × D w 0 ) × e 0 + ( e 1 + f 1 × D w 0 + ( g 1 + h 1 × D w 0 ) × e 0 ) × φ 0
v 0 = a 2 + b 2 × D w 0 + ( c 2 + d 2 × D w 0 ) × φ 0
h 0 = H 0 2 π R 3 γ s ,   v 0 = V ult 2 π R 3 γ s ,   D w 0 = D w D p L 2 ,   φ 0 = tan ( φ ) γ s D p s u m ,   e 0 = e D p .

4.3. Empirical Design Procedure and Validation

Equations (2)–(4) can be applied based on geometry parameters as DwDp/L2 = 0.05~0.1 and soil properties with sum = 20~40 kPa, φ = 30°~40°, which is in line with engineering practice [47,48,49]. The equations are calculated with Ts/L = 0.1, 0.3, 0.5 and 0.7, which can also be used for other values of Ts/L through interpolating. The failure envelopes calculated by FEM and Equations (2)–(4) are compared in Figure 10b for an example, and good agreement is obtained. The maximum relative error of about 7% occurs when Ts/L = 0.7 and V/(2πR3γs) is relatively small, which means that the prediction of the equations can be considered conservative. Based on the above geometric dimensions and soil properties, a suggested design procedure is proposed: (i) determining the shape of the failure envelopes in the V-H-M domain with the assistance of Equation (2) and Table 2; (ii) calculating the absolute values of uniaxial lateral bearing capacity and vertical capacity using Equations (3) and (4) along with Table 3 and Table 4; (iii) obtaining the lateral and moment bearing capacities under a certain vertical load of the composite foundation.
To validate the suitability and accuracy of the proposed empirical design procedure, it is used to estimate the bearing capacity of composite foundations reported by other researchers. Firstly, the centrifuge test data reported by Lehane et al. (2014) [27] on a composite foundation in pure sand are compared with the predicted results using the proposed method in this study (DwDp/L2 = 0.048, e/Dp = 8, V/Vult = 0). The lateral bearing capacities from the empirical formula and test results are 59.23 and 63.14 in normalized form, respectively, which indicates an error of only 6.61%. Then, to further validate the accuracy, the laboratory testing results from the centrifuge tests carried out by Stone et al. (2018) [36] are also compared with the prediction results using the proposed empirical formula in this work (DwDp/L2 = 0.094, e/Dp = 8, V/Vult = 0). The normalized lateral bearing capacities are 143.64 and 181.89 for the empirical design procedure and the centrifuge test, respectively. It reflects an error of 21.03%. The difference is probably because in pure sand, the shape of the failure envelope is slightly different to that of Ts/L = 0.7 in certain situations, and the empirical design procedure for the bearing capacities of the composite foundation in sand-over-clay deposit may not be perfectly suitable for that in single soil layer conditions. Since very limited results are available for composite foundations in layered soil conditions, and the proposed method agrees well with the laboratory testing data in this study, further study with more testing should be carried out to validate the proposed method. Nevertheless, based on the above comparison with FEM results and the centrifuge testing data, it is believed that the proposed empirical formula could be employed with good accuracy.

5. Conclusions

In this study, numerical simulations are carried out on the composite foundation in sand-overlying-clay deposits under combined V-H-M loadings. Based on the parametric analyses, the bearing capacities and failure mechanisms are found to be closely related to the foundation geometry (Dw/L and L/Dp), the soil properties (tan(φ)γsDp/sum), the upper sand layer thickness (Ts/L), the lateral loading height (e/Dp = 8), and the pre-vertical load (V/Vult). The following conclusions can be drawn:
For the bearing capacities in the V-H-M domain, the foundation geometry has a more pronounced influence when the wheel diameter (Dw/L) exceeds 0.5. However, it makes little difference to the shape of the failure envelopes under combined V-H-M loadings. The sand layer thickness (Ts/L) influences the bearing capacity with a nonlinear manner and the trend is determined by the soil properties (tan(φ)γsDp/sum), which have non-negligible effects on the failure mechanism and the shape of failure envelopes in the V-H-M domain. Due to the P-Δ effect and the “punch-through” failure pattern, the vertical load (V/Vult) influences the lateral bearing capacity with an optimal value when the maximum capacity occurs. Through imposing an additional bending moment onto the composite foundation, the lateral loading height (e/Dp) shows a negative effect on the bearing capacity and the stability of the foundation. In the end, an empirical design approach is proposed with Equations (2)–(4) and Table 2, Table 3 and Table 4 for the prediction of the composite foundation in sand-overlying clay deposit under the V-H-M combined loadings.

Author Contributions

Y.W. (conceptualization, formal analysis, methodology, validation, writing-original draft preparation); X.Z. (formal analysis, Funding acquisition, supervision, writing-review & editing); J.H. (writing-review & editing; validation). All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (Grant No. 51578231 & 52178329).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Notations

αsoil/pile friction coefficient
Dwdiameter of friction wheel
Dw0dimensionless foundation geometry
Dpdiameter of pile
dwdistance to the foundation center
elateral loading height
e0normalized lateral loading height
EcYoung’s modulus of clay
EsYoung’s modulus of sand
Hlateral load
Hultlateral capacity
H0uniaxial horizontal capacity
h0dimensionless uniaxial horizontal capacity
HBDvertical boundary width
K0horizontal earth pressure coefficient
kgradient of strength with depth
Lembedment length of pile
Mmoment
Pwsoil passive resistance under the wheel
Pvvertical pressure under the wheel
Ppsoil passive resistance along the pile
Rradius of the pile
suundrained shear strength of clay
sumundrained shear strength of clay at the mudline
Twthickness of friction wheel
Tsthickness of sand layer
Tcollarthickness of wheel collar
v0dimensionless vertical load
VBDvertical boundary distance
Vvertical load
Vultvertical capacity
γceffective unit weight of clay
γseffective unit weight of sand
zdepth below the mudline alongside the pile
zmaxdepth of maximum earth pressure
φinternal friction angle of sand fill
φ0normalized internal friction angle of sand
ψdilation angle
υPoisson’s ratio
A, B, C, D, Ecoefficients of Equation (2)
a1, b1, c1, d1, e1, f1, g1, h1coefficients of Equation (3)
a2, b2, c2, d2coefficients of Equation (4)

References

  1. Lee, M.; Bae, K.; Lee, I.; Yoo, M. Cyclic p-y Curves of Monopiles in Dense Dry Sand Using Centrifuge Model Tests. Appl. Sci. 2019, 9, 1641. [Google Scholar] [CrossRef] [Green Version]
  2. Zhang, Z.; Guan, P.; Xu, J.; Wang, B.; Li, H.; Dong, Y. Horizontal Loading Performance of Offshore Wind Turbine Pile Foundation Based on DPP-BOTDA. Appl. Sci. 2020, 10, 492. [Google Scholar] [CrossRef] [Green Version]
  3. Plodpradit, P.; Dinh, V.; Kim, K. Coupled Analysis of Offshore Wind Turbine Jacket Structures with Pile-Soil-Structure Interaction Using FAST v8 and X-SEA. Appl. Sci. 2019, 9, 1633. [Google Scholar] [CrossRef] [Green Version]
  4. Zhang, Y.; Andersen, K.; Tedesco, G. Ultimate bearing capacity of laterally loaded piles in clay—Some practical considerations. Mar. Struct. 2016, 50, 260–275. [Google Scholar] [CrossRef]
  5. Jeong, Y.; Park, M.; Kim, J.; Song, S. Wave Force Characteristics of Large-Sized Offshore Wind Support Structures to Sea Levels and Wave Conditions. Appl. Sci. 2019, 9, 1855. [Google Scholar] [CrossRef] [Green Version]
  6. Pandit, R.; Kolios, A. SCADA Data-Based Support Vector Machine Wind Turbine Power Curve Uncertainty Estimation and Its Comparative Studies. Appl. Sci. 2020, 10, 8685. [Google Scholar] [CrossRef]
  7. Schafer, M.; Madabhushi, S.P.G. Uplift Resistance of Enlarged Base Pile Foundations. Indian Geotech. J. 2020, 50, 426–441. [Google Scholar] [CrossRef]
  8. Yang, Y.; Qiu, L. MPM simulation of uplift resistance of enlarged base piles in sand. Soils Found. 2020, 60, 1322–1330. [Google Scholar] [CrossRef]
  9. Xu, H.; Yue, Z.; Qian, Q. Predicting uplift resistance of deep piles with enlarged bases. Proc. Inst. Civil Eng. Geotech. Eng. 2009, 162, 225–238. [Google Scholar] [CrossRef]
  10. Wang, X.; Zeng, X.; Li, J.; Yang, X. Lateral bearing capacity of hybrid monopile-footing foundation for offshore wind turbines by centrifuge modelling. Ocean Eng. 2018, 148, 182–192. [Google Scholar] [CrossRef]
  11. Yang, X.; Zeng, X.; Wang, X.; Berrila, J.; Li, X. Performance and bearing behavior of monopile-footing foundations under lateral-moment loading for offshore wind turbines. Ocean Eng. 2019, 184, 159–172. [Google Scholar] [CrossRef]
  12. Wang, X.; Zeng, X.; Yang, X.; Li, J. Seismic response of offshore wind turbine with hybrid monopile foundation based on centrifuge modelling. Appl. Energy 2019, 235, 1335–1350. [Google Scholar] [CrossRef]
  13. Carder, D.; Watson, G.; Chandler, R.; Powrie, W. Long-term performance of an embedded retaining wall with a stabilizing base slab. Proc. Inst. Civil. Eng. Geotech. Eng. 1999, 137, 63–74. [Google Scholar] [CrossRef]
  14. Mokwa, R.; Duncan, J. Experimental evaluation of lateral-load resistance of pile caps. J. Geotech. Geoenviron. Eng. 2001, 127, 185–192. [Google Scholar] [CrossRef]
  15. Rollins, K.; Cole, R. Cyclic lateral load behavior of a pile cap and backfill. J. Geotech. Geoenviron. Eng. 2006, 132, 1143–1153. [Google Scholar] [CrossRef]
  16. Wang, Z.; Hu, R.; Leng, H.; Liu, H.; Bai, Y.; Lu, W. Deformation Analysis of Large Diameter Monopiles of Offshore Wind Turbines under Scour. Appl. Sci. 2020, 10, 7579. [Google Scholar] [CrossRef]
  17. Yang, B.; Wei, K.; Yang, W.; Li, Y.; Qin, B.; Ning, L. A Feasibility Study for Using Fishnet to Protect Offshore Wind Turbine Monopile Foundations from Damage by Scouring. Appl. Sci. 2019, 9, 5023. [Google Scholar] [CrossRef] [Green Version]
  18. Zou, X.; Hu, Y.; Hossain, M.; Zhou, M. Capacity of skirted foundations in sand-over-clay under combined V-H-M loading. Ocean Eng. 2018, 159, 201–218. [Google Scholar] [CrossRef]
  19. Yang, C.; Zhu, Z.; Xiao, Y. Bearing Capacity of Ring Foundations on Sand Overlying Clay. Appl. Sci. 2020, 10, 4675. [Google Scholar] [CrossRef]
  20. Jan, M.; Adeel, M.; Pervaiz, U.; Park, D. Estimation of Distribution Factor for Peak Penetration Resistance Prediction of Spudcan Foundations in Loose to Medium-Dense Sand Overlying Clay. Appl. Sci. 2020, 10, 8795. [Google Scholar] [CrossRef]
  21. Hu, P.; Cassidy, M.; Randolph, M. Bearing capacity on sand overlying clay: An analytical model for predicting post peak behaviour. Mar. Struct. 2018, 59, 94–104. [Google Scholar] [CrossRef]
  22. Hossain, M.; Hu, P.; Cassidy, M.; Menzies, D.; Wingate, A. Measured and calculated spudcan penetration profiles for case histories in sand-over-clay. Appl. Ocean Res. 2019, 82, 447–457. [Google Scholar] [CrossRef]
  23. Zheng, G.; Zhao, J.; Zhou, H.; Zhang, T. Ultimate Bearing Capacity of Strip Footings on Sand Overlying Clay under Inclined Loading. Comput. Geotech. 2018, 106, 266–273. [Google Scholar] [CrossRef]
  24. Stone, K.; Newson, T.; Sandon, J. An Investigation of the Performance of a ‘Hybrid’ Monopile-footing Foundation for Offshore Structures. In Proceedings of the 6th International Offshore Site Investigation and Geotechnics Conference: Confronting New Challenges and Sharing Knowledge, London, UK, 11–13 September 2007; pp. 11–13. [Google Scholar]
  25. El-Marassi, M. Investigation of Hybrid Monopile-Footing Foundation Systems Subjected to Combined Loading. Ph.D. Thesis, The University of Western Ontario, London, ON, Canada, 2011. [Google Scholar]
  26. Lehane, B.; Powrie, W.; Doherty, J. Centrifuge model tests on piled footings in clay for offshore wind turbines. In Proceedings of the International Conference in Physical Modelling in Geotechnics, Zurich, Switzerland, 28 June–1 July 2010; Balkema: Rotterdam, The Netherlands, 2010. [Google Scholar]
  27. Lehane, B.; Pedram, B.; Doherty, J.; Powrie, W. Improved performance of monopiles when combined with footings for tower foundations in sand. J. Geotech. Geoenviron. Eng. 2014, 140, 04014027. [Google Scholar] [CrossRef]
  28. Anastasopoulos, I.; Theofilou, M. Piled footing for offshore wind turbines: Environmental and seismic loading. Soil Dyn. Earthq. Eng. 2016, 80, 192–209. [Google Scholar] [CrossRef]
  29. Abdelkader, A.; Naggar, E.; Hesham, M. Piled footing System for Offshore Wind Turbine. Geotech. Geol. Eng. 2018, 36, 2921–2937. [Google Scholar] [CrossRef]
  30. Pedram, B. Behaviour of Hybrid Piled Footing Structures in Sands. Geotech. Geol. Eng. 2018, 36, 2273–2292. [Google Scholar] [CrossRef]
  31. Wang, X.; Li, J. Parametric study of hybrid monopile foundation for offshore wind turbines in cohesionless soil. Ocean Eng. 2020, 218, 108172. [Google Scholar] [CrossRef]
  32. Yang, X.; Zeng, X.; Wang, X.; Yu, H. Performance of monopile-footing foundations under lateral loading for offshore wind turbines. Appl. Ocean Res. 2018, 78, 14–24. [Google Scholar] [CrossRef]
  33. Trojnar, K. Multi scale studies of the new piled footings for offshore wind turbines. Ocean Eng. 2019, 192, 106506. [Google Scholar] [CrossRef]
  34. Trojnar, K. Simplified design of new hybrid monopile foundations for offshore wind turbines. Ocean Eng. 2021, 219, 108046. [Google Scholar] [CrossRef]
  35. Li, X.; Zeng, X.; Wang, X. Feasibility study of monopile-footing-bucket piled footing for offshore wind turbine. Ocean Eng. 2020, 204, 107276. [Google Scholar] [CrossRef]
  36. Stone, K.; Arshi, H.; Zdravkovic, L. Use of a bearing plate to enhance the lateral capacity of monopiles in sand. J. Geotech. Geoenviron. Eng. 2018, 144, 04018051. [Google Scholar] [CrossRef]
  37. Dassault Systèmes. Abaqus Analysis User’s Manual; Simulia Corp.: Santa Clara, CA, USA, 2016. [Google Scholar]
  38. Hao, E.; Liu, C. Evaluation and comparison of anti–impact performance to offshore wind turbine foundations: Monopile, tripod, and jacket. Ocean Eng. 2017, 130, 218–227. [Google Scholar] [CrossRef]
  39. Moulas, D.; Shafiee, M.; Mehmanparast, A. Damage analysis of ship collisions with offshore wind turbine foundations. Ocean Eng. 2017, 143, 149–162. [Google Scholar] [CrossRef]
  40. Malhotra, S. Selection, design and construction of offshore wind turbine foundations. In Wind Turbines; Al-Bahadly, I., Ed.; InTechOpen: London, UK, 2011. [Google Scholar]
  41. Al-Mhaidib, A.I. Experimental investigation of the behavior of pile groups in sand under different loading rates. Geotech. Geol. Eng. 2006, 24, 889–902. [Google Scholar] [CrossRef]
  42. Hung, L.; Kim, S. Evaluation of undrained bearing capacities of bucket foundations under combined loads. Mar. Georesour. Geotechnol. 2014, 32, 76–92. [Google Scholar] [CrossRef]
  43. Zhou, M.; Hossain, M.; Hu, Y.; Liu, H. Behaviour of ball penetrometer in uniform single-and double-layer clays. Géotechnique 2013, 63, 682–694. [Google Scholar] [CrossRef]
  44. Vulpe, C.; Gourvenec, S.; Leman, B.; Fung, K. Failure envelope approach for consolidated undrained capacity of shallow foundations. J. Geotech. Geoenviron. Eng. 2016, 142, 04016036. [Google Scholar] [CrossRef] [Green Version]
  45. Gourvenec, S.; Barnett, S. Undrained failure envelope for skirted foundations under general loading. Géotechnique 2011, 61, 263–270. [Google Scholar] [CrossRef]
  46. Gourvenec, S.; Mana, D. Undrained vertical bearing capacity factors for shallow foundations. Géotech. Lett. 2011, 1, 101–108. [Google Scholar] [CrossRef]
  47. Randolph, M.; Gourvenec, S.; Acosta-Martinez, H. Experimental Study of Uplift Resistance of Shallow Skirted Foundations in Clay under Transient and Sustained Concentric Loading. Géotechnique 2009, 59, 525–537. [Google Scholar]
  48. Mahmoud, G.; Javad, N. Bearing capacity of geosynthetic encased stone columns. Geotext. Geomembr. 2013, 38, 26–36. [Google Scholar]
  49. Jung, S.; Kim, S.; Patil, A.; Hung, L. Effect of monopile foundation modeling on the structural response of a 5-MW offshore wind turbine tower. Ocean Eng. 2015, 109, 479–488. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic of single pile-friction wheel composite foundation.
Figure 1. Schematic of single pile-friction wheel composite foundation.
Applsci 11 09446 g001
Figure 2. Geometry and mesh of the single pile-friction wheel composite foundation models.
Figure 2. Geometry and mesh of the single pile-friction wheel composite foundation models.
Applsci 11 09446 g002
Figure 3. Lateral load-deflection curves obtained from FEM and centrifuge tests.
Figure 3. Lateral load-deflection curves obtained from FEM and centrifuge tests.
Applsci 11 09446 g003
Figure 4. Effect of wheel diameter on the lateral capacity and the failure mechanism.
Figure 4. Effect of wheel diameter on the lateral capacity and the failure mechanism.
Applsci 11 09446 g004
Figure 5. Effect of the pre-vertical loading on the lateral capacity and the failure mechanism with different loading height.
Figure 5. Effect of the pre-vertical loading on the lateral capacity and the failure mechanism with different loading height.
Applsci 11 09446 g005
Figure 6. Effect of wheel diameter on the lateral capacity and the failure mechanism.
Figure 6. Effect of wheel diameter on the lateral capacity and the failure mechanism.
Applsci 11 09446 g006
Figure 7. Lateral soil resistance along the pile under different sand layer thicknesses.
Figure 7. Lateral soil resistance along the pile under different sand layer thicknesses.
Applsci 11 09446 g007
Figure 8. Effect of wheel thickness and height of wheel collar.
Figure 8. Effect of wheel thickness and height of wheel collar.
Applsci 11 09446 g008
Figure 9. Effects of soil properties (tan(φ)γsDp/sum).
Figure 9. Effects of soil properties (tan(φ)γsDp/sum).
Applsci 11 09446 g009
Figure 10. Effects of wheel diameter and pile embedment depth (DwDp/L2).
Figure 10. Effects of wheel diameter and pile embedment depth (DwDp/L2).
Applsci 11 09446 g010
Figure 11. Effects of sand layer thickness (Ts/L).
Figure 11. Effects of sand layer thickness (Ts/L).
Applsci 11 09446 g011
Figure 12. Typical FE results of the normalized bearing capacity in the V-H-M domain.
Figure 12. Typical FE results of the normalized bearing capacity in the V-H-M domain.
Applsci 11 09446 g012
Table 1. Summary of the numerical simulations.
Table 1. Summary of the numerical simulations.
AnalysisDw/LV/Vulte/DpTs/Lsumφ (sand)Tw/DpNotes
Group I0, 0.5, wheel5.9, 12.7, 6.4, 13.2 (MN)91033°0.8Numerical model validation [27]
Group II0~0.80100.530 kPa33°0.5Investigation for the effect of the wheel diameter
Group III0.60~95%6~160.530 kPa33°0.5Investigation for the effect of the combined pre-vertical load and moment
Group IV0.60100.1~0.820 kPa, 40 kPa
(k = 1.75 and 2 kPa)
20~40°0.5Investigation for the effect of thickness of sand layer
Group V0.60100.530 kPa33°0.5~1.25Investigation for the effect of wheel thickness
Table 2. Values of A, B, C, D and E.
Table 2. Values of A, B, C, D and E.
Ts/Ltan(φ)γsDp/sumABCDE
0.10.75411.177−0.5050.4393.784
0.10.84811.847−0.7900.20927.125
0.10.94912.200−0.7900.21045
0.11.0960.7391.913−0.5510.7744.49
0.30.75412.1720.5081.4110.194
0.30.84812.82.7253.6590.201
0.30.94911.850.3511.3071.55
0.31.09611.3550.3221.3142.552
0.50.75411.2560.9141.6530.891
0.50.84811.83114.90715.6650.084
0.50.94910.9610.3051.3003.593
0.51.09611.7181.5982.6231.621
0.70.7540.3020.0150.2033.4224.655
0.70.8480.6020.3210.1611.7944.572
0.70.9490.8740.8690.3141.4483.357
0.71.0961.0631.7090.8891.8623.471
Table 3. Values of parameters related to Ts/L in Equation (3) (h0).
Table 3. Values of parameters related to Ts/L in Equation (3) (h0).
Ts/La1b1c1d1e1f1g1h1
0.11.61230−0.21−2.78−1.22399.440.19−12.22
0.3−15.99457.220.53−11.67−16.99957.221.33−37.78
0.5−2.25164.440.07−2.220.09948.33−0.035−25.56
0.77.52−260.56−0.615−9.021506.670.54−45
Table 4. Values of parameters related to Ts/L in Equation (4) (v0).
Table 4. Values of parameters related to Ts/L in Equation (4) (v0).
Ts/La2b2c2d2
0.1−72.571876.67−11.97997.22
0.3−101.012206.11−21.932400
0.5−74.471263.89−3.823653.89
0.7−136.691598.3342.764313.89
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Y.; Zou, X.; Hu, J. Bearing Capacity of Single Pile-Friction Wheel Composite Foundation on Sand-over-Clay Deposit under V-H-M Combined Loadings. Appl. Sci. 2021, 11, 9446. https://0-doi-org.brum.beds.ac.uk/10.3390/app11209446

AMA Style

Wang Y, Zou X, Hu J. Bearing Capacity of Single Pile-Friction Wheel Composite Foundation on Sand-over-Clay Deposit under V-H-M Combined Loadings. Applied Sciences. 2021; 11(20):9446. https://0-doi-org.brum.beds.ac.uk/10.3390/app11209446

Chicago/Turabian Style

Wang, Yikang, Xinjun Zou, and Jianfeng Hu. 2021. "Bearing Capacity of Single Pile-Friction Wheel Composite Foundation on Sand-over-Clay Deposit under V-H-M Combined Loadings" Applied Sciences 11, no. 20: 9446. https://0-doi-org.brum.beds.ac.uk/10.3390/app11209446

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