Investigation of Pullout Capacity Characteristics of Suction Anchors Under Inclined Loads in Layered Soil

1. Introduction

Suction anchors typically refer to cylindrical structures that are closed at the top and open at the bottom, commonly constructed of steel or concrete [1,2]. Its diameter is typically 3–7 meters, with wall thickness ranging from 25 to 75 millimeters, and the L/D is typically between 1 and 10 [3].Suction anchors are widely used in coastal floating platforms, offshore wind turbine foundations, and deepwater projects due to their rapid installation, low cost, high load-bearing capacity, and minimal disturbance to seabed soil layers [4,5]. During the installation of suction anchors, the anchors are initially driven into the seabed soil using their own weight. Subsequently, water is pumped out through the suction holes at the top of the anchor, which creates a pressure differential between the interior and exterior of the anchor. This pressure differential forces the anchor to continue driving into the soil. When the top cover of the suction anchor makes contact with the seabed soil, the anchor is considered fully penetrated [6,7].

In actual engineering environments, offshore structures inevitably experience the effects of waves, wind, and currents [8,9]. Since the load must be transmitted to the suction anchor via the suspension chain and mooring chain, the suction anchor is continuously subjected to inclined loads [10]. The bearing mechanism of suction anchors subjected to inclined pullout loads differs from that of suction anchor foundations bearing vertical loads at the top. The loading angle and padeye location alter the failure mode of suction anchors [11]. Keaveny et al. [12]conducted prototype tests on a series of suction anchors installed in clay, and applied mooring forces to the padeyes of the anchors below the mud surface. The test results indicate that when the mooring force is applied at the midpoint below the mud surface of the suction anchor, the corresponding ultimate bearing capacity is twice that when applied at the mud surface. Newlin [13] studied that when mooring forces are applied below the optimal padeye position, the top of the suction anchor tilts away from the mooring force direction upon failure. In this case, the failure mode is not pullout failure but rotation within the soil. The suction anchor can obtain its maximum bearing capacity only when the anchor undergoes translational displacement. Therefore, the loading angle and padeye are critical factors affecting suction anchors.

To explore the ultimate bearing capacity of suction anchors under inclined loads, numerous studies have been conducted through theoretical analysis, numerical simulation, and physical testing. Current research primarily focuses on suction anchors in homogeneous soil. Wang et al. [14] conducted load-controlled model tests on the ultimate bearing capacity of suction anchors in soft clay, and analyzed the effects of loading direction on failure modes and ultimate bearing capacity. The test results indicate that the bearing capacity of suction anchors is influenced by the loading direction. The reverse bearing capacity of the soil layer below the anchor base is a key factor affecting the ultimate bearing capacity of suction anchors. Zhang et al. [15] investigated the effect of cutting grooves on the ultimate bearing capacity of suction anchors in NC through numerical simulation methods. Kim et al. [16] investigated the behavior of suction anchor foundations installed in cohesionless soil under horizontal loading through numerical simulation. They conducted parametric studies on different length to diameter ratios and padeye position of the suction anchors. The horizontal bearing capacity of the suction anchors was estimated, and soil stress distribution was analyzed under loading applied at the optimal padeye position. Existing studies on the bearing capacity of suction anchors have primarily focused on single homogeneous soil layers, which lack consideration of the actual layered characteristics of soil masses [17,18]. However, under the influence of geological movements and sedimentation processes, the actual seabed typically exhibits natural layering or anisotropic properties along the depth direction. The characteristics of such layered soil differ significantly from those of the current single-layer homogeneous soil layers studies [19,20]. In natural layered soils, a typical layered soil structure occurs where stronger soil layers overlie weaker layers [21]. For example, in the Bohai oil and gas fields off of China, the seabed exhibits a layered structure with upper- stiff-lower-soft layers depending on factors such as estuarine sediment transport patterns and tidal currents [22,23]. Additionally, in areas such as the Gulf of Guinea, the Bass Strait off southeastern Australia, and the Joa oil field in the North Sea, the seabed also contains hard soil layers with significantly greater strength than the underlying soil [24,25]. Furthermore, existing design guidelines for piles in layered soils primarily employ the equivalent lateral friction coefficient of the piles and the p–y curve method for lateral friction resistance. This approach is unsuitable for predicting the pullout resistance of suction anchors in layered soils. Existing guidelines, such as those from the American Petroleum Institute (API) [26] and DNV [27], typically assume homogeneous and uniform soil properties. Applying these existing guidelines to predict bearing capacity under layered soil conditions may be overly conservative. Therefore, it is important to fully consider the influence of layered soil properties to gain a deeper understanding of the pullout resistance capacity of suction anchors.

Based on the above-mentioned research status, this study aims to investigate the ultimate bearing capacity of suction anchors in upper- stiff-lower-soft layered clay. To analyze the bearing capacity of suction anchors under inclined loads, numerical simulation methods were employed. Considering variations in clay strength (S_um/kD) and layer thickness ratios (T_h/L、T_c/L、T_n/L), the ultimate bearing capacity of suction anchors subjected to inclined loads in layered clay (lightly overconsolidated clay (LOC), normally consolidated clay (NC), low-strength normally consolidated clay (LNC)) was investigated. The primary research content includes: (1) Determined the optimal centerline loading depth Z_cl，opt in layered clay; (2) Obtained the corresponding uniaxial ultimate bearing capacities H_ult and V_ult; (3) Quantified the suction anchor VH failure envelope under inclined loading and derived an analytical expression for its lower bound solution.

2. Numerical Modeling

2.1. Geometry and Mesh

By referring the studies of Kim et al. [28] and Fu et al. [29], the numerical simulation model selected the suction anchor diameter D=5m, depth of embedment L=7.5m, with L/D=1.5 and thickness t =50mm. The selected suction anchor dimensions are commonly adopted in engineering practice [30,31]. As shown in Figure 1(a), the suction anchor is installed in layered clay. As the stiffness of suction anchors and soil differs significantly, suction anchors are typically treated as rigid bodies to improve computational efficiency. This method is widely applied in offshore geotechnical engineering [28,29,32]. The suction anchor is modeled using rigid materials with a density of ρ = 6800 kg/m³, elastic modulus E = 210 GPa, and Poisson’s ratio v = 0.3. An inclined load is applied by setting a reference point on the suction anchor.

In this study, a three-dimensional finite element model was constructed using the commercial software ABAQUS 2020. The suction anchor and soil mass were modeled using three-dimensional eight-node solid elements (C3D8H). Due to load symmetry, only half of the soil mass was modeled to improve computational efficiency. The minimum mesh size for the model is 0.02D. Based on existing research and sensitivity analysis of mesh sizes, the selected minimum mesh size is reasonable [33,34]. The length and height of the soil are B_H=10D and B_V=3L, respectively, with the suction anchor positioned at the center of the soil model. The model dimensions are selected sufficiently to eliminate the influence of boundary effects [18,35]. The relevant constraints on the soil are shown in Figure 1(b).

2.2. Geometry and Mesh

The clay in this study is treated as an elastoplastic material, adhering to the Tresca yield criterion and associated flow rules. The Poisson’s ratio of the clay is v=0.49, which is used to simulate undrained conditions where volume remains constant. The undrained shear strength of the clay is satisfied:

$$S_u=S_{um}+kz$$

(1)

Where, $S_{um}$ represents the undrained shear strength of the soil at the mud surface, $k$ represents the gradient of clay strength with depth, and $z$ represents the depth below the mud surface. The soil elastic modulus is taken as E=500S_u, and the normalized parameter S_um/kD is used to evaluate the influence of soil properties on bearing capacity. Referencing relevant studies and geological exploration reports, S_um/kD = 0, 0.4, 0.67, 1.0 [36,37,38,39], the relevant parameters of the clay are shown in Table 1. The numerical simulation tests for suction anchors are outlined in Table 2.

To investigate the bearing capacity of suction anchors in upper- stiff-lower-soft layered clay under inclined loading, the soil in this study was divided into three layers. The first layer consists of LOC with S_um/kD=1.0, and layer thickness ratios T_h/L=0.5, 0.75 and 1.25. The second layer comprised either NC or LOC, with S_um/kD values of 0.4, 0.67, and 1.0, and layer thickness ratios T_c/L =0.5, 0.75 and 1.25. The third layer consists of LNC, with S_um/kD = 0, and a soil layer thickness ratio T_n/L=(3L-T_h-T_c)/L. The relevant parameters were obtained by referring to related research and geological exploration reports. Regarding the bearing capacity of suction anchors in clay, due to the low permeability coefficient of clay, numerous studies have assumed the soil as undrained conditions[34,40,41,42].

2.3. Details of Analysis and Loading Methods

The finite element modeling process consists of three steps:

(1) The initial geostress equilibrium stage: This section performs geostress equilibrium analysis of the soil mass only, and the suction anchor is treated as void space. In the initial geostress equilibrium stage, the anchor-soil interface contact is not activated.

(2) The contact activation stage: After the soil stress equilibrium is established, the anchor-soil interface contact is activated. This means the penetration process of the suction anchor is disregarded, and the suction anchor is considered to be in the “in-place” state.

(3) The load application stage: A reference point RP is established at the suction anchor to apply the inclined load. The inclined pull-out load is applied under displacement control, with a maximum allowable displacement of 0.2D. The selection of the maximum displacement references the research by Supachawarote[32].

The tangent method was employed to determine the uniaxial bearing load capacity, by referring to the study of Long et al. [43], for which no distinct yield value was observed when the load-displacement curve reached its maximum displacement. When the suction anchor moves horizontally or vertically without rotation, the uniaxial bearing capacity is equivalent to the uniaxial ultimate bearing capacity H_ult and V_ult. The complete VH failure envelopes for the suction anchor can be obtained through sweep analysis, reverse sweep analysis, and probe analysis, which has been widely applied[32].

As shown in Figure 2, when the suction anchor moves horizontally under inclined loading without rotation, the inclined load can be simplified into a horizontal force H, a vertical force V, and a constraining moment M. By applying the principle of moment equilibrium, the optimal padeye depth Z_pad,opt at the corresponding angle and the optimal centerline loading depth Z_cl,opt can be obtained, which satisfy the equation:

$$Z_{cl,opt}={\displaystyle \frac{M}{H}}$$

(2)

$$Z_{cl,opt}=Z_{pad,opt}+{\displaystyle \frac{D}{2}}\times \tan \theta$$

(3)

3. Model Validation

To validate the accuracy of the constructed numerical model, the numerical simulation results of suction anchors under inclined loading were compared with existing experimental results, as shown in Figure 3. It can be seen from Figure 3 (a), for NC with S_u=1.25z and L/D=1.5, the numerical simulation results show excellent consistency with test results from Supachawarote [32] and experimental results from the Norwegian Geotechnical Institute (NGI), Offshore Technology Research Center (OTRC), and Center for Offshore Foundation Systems (COFS). The VH envelope obtained from the numerical model closely matches their findings, which indicates that the numerical model developed in this study can effectively reproduce the pullout behavior of suction anchors in clay.

Furthermore, the accuracy of the constructed numerical simulation model was further validated by referencing centrifuge test data and finite element modeling results for non-cohesive soils from Kim et al. [17,28].The Loads were applied at the top of the suction anchor and below the mud surface at Z_pad=2/3L, with cu=7 kPa, internal friction angle φ=32°, and dilation angle ψ=8°. The selection of relevant parameters strictly followed those used in the centrifuge tests and finite element model. The corresponding comparison results are shown in Figure 3(b), which shows that the constructed numerical model can effectively characterize the bearing capacity of the suction anchor.

4. Results of Numerical Analyses

4.1. Optimum Centerline Loading Depth Z_cl，opt

The optimal centerline loading depth Z_cl,opt is the key factor determining whether the suction anchor undergoes translational or rotational motion under inclined loading. When the inclined load vector applied to the suction anchors padeye passes through Z_cl,opt, the suction anchors move horizontally under the inclined load without rotation. At this point, the inclined load corresponds to the maximum pullout force at the given mooring angle. When the inclined load vector applied to the anchor padeye deviates from Z_cl,opt, the suction anchor will rotate forward or backward during the pullout process, correspondingly reducing the pullout resistance of the suction anchor. Therefore, Z_cl,opt is an indicator reflecting the pullout resistance capacity and failure mode of the suction anchor under corresponding soil parameter conditions. In engineering design, the location of Z_cl,opt must be prioritized for determination.

Figure 4 shows the influence of different soil mechanical parameters and layer thickness on the optimal centerline loading depth Z_cl,opt for layered clay. It can be seen that the LOC thickness ratio T_h/L is the key factor controlling Z_cl,opt. Taking S_um/kD = 0.4 as an example, when T_h/L = 0.5, half of the suction anchor is embedded in the LOC layer, with Z_cl,opt = 0.65. With the increase in T_h/L to 0.75 and 1.25, the upper portion of the suction anchor contributes more lateral friction resistance due to the increased thickness of the LOC layer. To balance the moments, the suction anchors Z_cl,opt shift downward, with corresponding Z_cl,opt being 0.67 and 0.71, respectively. When T_h/L > 1, the suction anchor is fully embedded in the LOC layer. The values of Z_cl,opt remain nearly unchanged for different S_um/kD and T_c/L, as the high strength of the LOC soil dominates the pullout resistance. When T_h/L = 0.5 and T_c/L = 0.5, the suction anchor end is located between the NC layer and the LNC layer interface. It can be observed that the variation of Z_cl,opt at this point is still influenced by the higher-strength clay within the suction anchor, which means that the impact of high-strength clay on Z_cl,opt is significantly greater than that of low-strength clay.

4.2. Uniaxial Ultimate Bearing Capacity (H_ult and V_ult)

The uniaxial ultimate bearing capacity (H_ult and V_ult) represents the pullout resistance of the suction anchors when subjected to horizontal or vertical pullout forces without rotation. In numerical simulations, the uniaxial ultimate bearing capacity equals the integral of the ultimate lateral friction resistance along the anchorage depth. Based on the above principles, the effects of suction anchor relative positions (T_h/L and T_c/L) and soil strength (S_um/kD) on uniaxial ultimate bearing capacity were quantified. The uniaxial ultimate bearing capacities H_ult and V_ult of suction anchors under different conditions are shown in Table 3, Table 4 and Table 5.

As shown in Table 3, Table 4 and Table 5, the LOC thickness ratio T_h/L is the core factor determining the uniaxial ultimate bearing capacity. Specifically, the undrained shear strength of LOC is significantly higher than that of NC and LNC. As the embedment length of the suction anchor into the LOC increases, the combined effect of lateral friction resistance and end resistance at the anchor-soil interface becomes more pronounced. Taking S_um/kD = 0.67 and T_c/L = 0.5 as examples, when T_h/L increases from 0.5 to 1.25, the corresponding horizontal ultimate bearing capacity H_ult of the suction anchor rises from 4182.01 kN to 5096.61 kN, representing an increase of approximately 21.9%. The vertical ultimate bearing capacity V_ult increases from 4457.01 kN to 6566.56 kN, representing an increase of approximately 47.3%. This supports the findings of Zhou et al. [18,43], which concluded that the pullout resistance of suction anchors in layered soils is positively correlated with the thickness of the anchor embedded in the higher undrained strength clay layer. As the anchor embedment depth in the higher undrained strength clay layer increases, both lateral friction resistance and end resistance enhance during anchor-soil interface interaction under pullout forces.

The corresponding displacement contour map (Figure 5) also shows that as T_h/L increases from 0.5 to 1.25, the suction anchors experience passive earth pressure on the external wall under horizontal loading, which results in the formation of a passive earth pressure zone. The extent of this passive earth pressure zone is primarily influenced by the length L of the suction anchor. When the suction anchor length remains constant, the horizontal ultimate bearing capacity is primarily influenced by the strength of the clay surrounding the anchor.

When subjected to vertical loads, the pullout resistance of suction anchors is primarily composed of the anchor-soil lateral friction resistance and end resistance. With the increase in the thickness ratio T_h/L of the LOC soil layer, Figure 5 shows that the end resistance influence zone expands. This means the contribution of end resistance to the vertical bearing capacity of the suction anchor increases, and consequently, the corresponding vertical ultimate bearing capacity V_ult increases. The vertical ultimate bearing capacity is primarily influenced by the thickness ratios T_h/L and T_c/L.

The layered clay makes the soil strength along the clay surface exhibit non-uniformity. The increase in S_um/kD enhances the uniaxial bearing capacity of layered clay. Taking T_h/L=0.5 and T_c/L=0.5 as examples, when S_um/kD increases from 0.4 to 1.0, H_ult rises from 3990.1 kN to 5017.5 kN, an increase of approximately 25.7%; V_ult increases from 4195.2 kN to 5341.32 kN, an increase of approximately 27.3%. The higher the S_um/kD , the greater the overall shear strength of the layered soil, and both the friction resistance and end resistance at the anchor-soil interface will increase simultaneously. The soil strength S_um/kD is positively correlated with bearing capacity; the higher the soil strength, the greater the corresponding bearing resistance. The thickness of the LNC soil layer has no significant effect on H_ult and V_ult. The displacement vector contour plots for layered soil subjected to horizontal/vertical loads are shown in Figure 6.

It can be observed that when the thickness ratios T_h/L and T_c/L remain unchanged, the corresponding horizontal/vertical load-displacement vector plots show little variation. This indicates that at this point, the bearing capacity of the suction anchor is primarily influenced by the friction resistance at the anchor-soil interface.

The thickness ratio T_c/L of the NC soil layer has little effect on the horizontal ultimate bearing capacity H_ult but significantly affects the vertical ultimate bearing capacity V_ult. For example, when S_um/kD=1.0 and T_h/L=0.5, as T_c/L increases from 0.5 to 1.25, H_ult changes little, rising only from 5017.5 kN to 5085.19 kN. Meanwhile, V_ult increases from 5341.32 kN to 7211.5 kN, representing an increase of approximately 35.0%. Combining Figure 5 and 6 with the studies by Supachawarote [32]: Vertical bearing capacity relies more on the accumulated effect of friction resistance at the anchor-soil interface. As the strength of the clay in the suction anchorage zone increases, the friction resistance becomes more pronounced. In contrast, horizontal bearing capacity is dominated by the distribution pattern of lateral resistance. The influence of T_c/L on H_ult is smaller than that of T_h/L and S_um/kD, which leads to a relatively minor overall change in H_ult.

When T_h/L > 1, which corresponded to the suction anchor being fully embedded in the LOC layer, the horizontal ultimate bearing capacity H_ult tended toward stability. Meanwhile, the vertical ultimate bearing capacity V_ult still shows a slight increase, which reflects that the vertical bearing capacity is more complex affected by the coupling of multiple soil layers. The lower soil layers (NC, LNC) still contribute to the end resistance.

4.3. The VH Failure Envelopes

The VH failure envelopes for layered soil are shown in Figure 7. By dimensionless processing of horizontal and vertical bearing capacities using H₀/γ’D and V₀/γ’D, respectively, the effects of parameters such as suction anchor dimensions and soil weight are eliminated.

As shown in Figure 7, with the increase in the LOC layer thickness ratio T_h/L (from 0.5 to 1.25), the horizontal bearing capacity of the VH failure envelopes significantly increases under the same vertical load. This indicates that the total area of the VH failure envelopes expands, which means the corresponding combined uniaxial bearing capacity improves. The LOC soil layer exhibits relatively high strength. As the thickness ratio T_h/L increases, the contact area between the suction anchor and the high-strength clay expands. This allows the anchor-soil interface to provide more uniform and substantial lateral friction resistance. Consequently, when horizontal loads increase, the attenuation of vertical loads becomes more gradual. The research findings align with the conclusions of Zhao et al. [44], indicating that the embedment length of suction anchors in high-strength soil layers is positively correlated with the combined bearing capacity.

The VH failure envelopes in layered clay exhibit asymmetry, which indicates that layered clay has a distinct influence on the horizontal and vertical bearing capacities of suction anchors. The horizontal friction resistance of the suction anchors dominates when horizontal loads dominate, with the thickness of the LOC layer determining the horizontal bearing capacity. When vertical loads dominate, the friction resistance at the anchor-soil interface and the end resistance play key roles, while the vertical bearing capacity is influenced by both the LOC and the underlying NC and LNC layers.

To visually represent the relative shape of the VH failure envelope, the results can be normalized based on the uniaxial bearing capacity H₀ and V₀. Extensive research indicates that the only lower bound solution of the VH failure envelope can be observed through relevant processing, which shows that the normalized VH failure envelope is independent of soil layer thickness and soil properties[45,46]. This conservative solution method allows the yield envelope in the VH envelope to be obtained simply by scaling H₀ and V₀. Based on previous research findings and regression analysis, the lower bound of the VH failure envelope, as shown in Figure 8, and can be expressed as:

$${\displaystyle \frac{H_0}{H_{ult}}}={\left\{1-{\left({\displaystyle \frac{V_0}{V_{ult}}}\right)}^3\right\}}^{1/0.85}$$

(4)

Where, ${\displaystyle \frac{V_0}{V_{ult}}}$ represents the normalized vertical bearing capacity; ${\displaystyle \frac{H_0}{H_{ult}}}$ represents the normalized horizontal bearing capacity. The exponent 1/0.85 indicates that vertical bearing capacity reduces the capacity of horizontal bearing capacity. Meanwhile, the ultimate bearing capacity of suction anchors installed in layered clay can be obtained by multiplying the corresponding ultimate bearing capacities H_ult and V_ult by Equation (4).

5. Discussion

In this paper, the bearing capacity of suction anchors subjected to inclined loads in upper- stiff-lower-soft layered clay is systematically investigated. Numerical simulation results indicate that the clay strength S_um/kD and the thickness ratio T/L are the core factors determining the bearing capacity of layered clay. To investigate the difference in bearing capacity of suction anchors in layered clay and homogeneous clay under inclined loading, the same numerical model shown in Figure 1 was employed to analyze the pullout bearing capacity of suction anchors for different clay strengths S_um/kD (Clay A to Clay D) as listed in Table 1. The optimum centerline loading depth, uniaxial ultimate load capacity, and VH failure envelopes under corresponding operating conditions were investigated.

It can be seen from Figure 9 that for homogeneous clay with varying S_um/kD, as the clay strength increases, Z_cl,opt decreases from 0.71 to 0.69. The increase in S_um/kD indicates a higher proportion of shear strength near the mud surface, shifting the center of gravity of the anchor shear strength upward. According to the principle of moment equilibrium, to prevent anchor rotation, Z_cl,opt shifts slightly upward. Overall, Z_cl,opt had little change.

Compared with the results in Section 4.1, it can be observed that the homogeneous clay Zcl,opt is solely influenced by S_um/kD, which has a more stable variation pattern. For layered clays, Z_cl,opt is primarily affected by the LOC layer, while the influence of the NC and LNC layers on Zcl,opt is masked by the high-strength LOC layer. This indicates that layered soils differ from homogeneous soils in their single correlation of ‘S_um/kD-Z_cl,opt,’ with the thickness of high-strength soil layers being the core variable controlling Z_cl,opt. By comparing the findings of Fu et al. [29] on homogeneous clay Z_cl,opt, the accuracy of the constructed numerical model can be further validated.

In combination with Figure 10, it can be observed that as S_um/kD increases, the uniaxial ultimate bearing capacity grows synchronously. However, the rate of increase varies: when S_um/kD rises from 0.4 to 1.0, the horizontal ultimate bearing capacity H_ult increases from 2915.65 kN to 5095.8 kN (an increase of approximately 74.8%), while the vertical ultimate bearing capacity V_ult increases from 3976.02 kN to 7317.5 kN (an increase of approximately 84.0%). The increase in H_ult is more pronounced. This difference stems from the distinct mechanisms governing bearing capacity. Horizontal bearing capacity relies on the uniform mobilization of lateral friction resistance throughout the entire anchor length; as S_um/kD increases, lateral friction resistance at all anchor depths rises synchronously. In contrast, vertical bearing capacity is simultaneously influenced by end resistance, this also results in a higher increment of vertical ultimate bearing capacity than that of horizontal ultimate bearing capacity.

In comparison with Section 4.2, it is evident that H_ult in layered clay is largely unaffected by the thicknesses of NC and LNC layers. When T_c/L increases from 0.5 to 1.25, H_ult remains nearly constant, being controlled solely by the LOC layer thickness T_h/L; whereas V_ult is influenced by both T_h/L and T_c/L. When T_c/L increases to 1.25, V_ult increases by approximately 35%. This contrasts with the behavior in homogeneous clay, where both H_ult and V_ult are governed by a single S_um/kD. Horizontal bearing capacity is more sensitive to the upper high-strength soil layer, while vertical bearing capacity is influenced by the coupled effects of multiple soil layers.

The VH failure envelopes for homogeneous clay are shown in Figure 11. By comparing with Section 4.3, it can be observed that the VH envelopes of layered clay exhibits distinct differences from that of homogeneous clay. Specifically, when horizontal bearing capacity is high, the envelope retains relatively high vertical bearing capacity, with the LOC layer providing significant lateral friction resistance. Conversely, when vertical bearing capacity is high, the envelope corresponds to a decrease in bearing capacity, where the anchor tip may be located in the NC or LNC layer, resulting in reduced end resistance.

A comprehensive analysis reveals significant differences in the pullout resistance of suction anchors in homogeneous clay and layered clay. While extensive research has focused on the bearing capacity of suction anchors in homogeneous clay, it is essential to investigate their pullout resistance in layered clay.

6. Limitations

This study conducted a series of pullout bearing capacity analyses for suction anchors with a commonly used L/D=1.5 in engineering applications, specifically within upper- stiff-lower-soft clay. It systematically investigated the effects of upper- stiff-lower-soft layered clay on the optimal centerline loading depth Z_cl,opt, uniaxial ultimate bearing capacity H_ult and V_ult, and VH failure envelopes of suction anchors. The findings provide theoretical basis and practical guidance for optimizing the design and safety assessment of suction anchor foundations in layered clay. However, certain limitations remain:

(1) This paper investigates the pullout bearing capacity of suction anchors subjected to inclined loads in layered clay with upper- stiff-lower-soft. In practical engineering, layered clay exhibits diverse characteristics such as upper- stiff-lower-soft, upper-soft-lower hard, soft-hard-soft, and hard-soft-hard etc. Therefore, based on existing research, a systematic investigation into the bearing capacity of suction anchors in layered clay is necessary.

(2) This paper investigates the pullout bearing capacity of suction anchors with a L/D =1.5 in layered clay. In practical engineering applications, the L/D of suction anchors varies according to project design (ranging from L/D=1 to 10). Therefore, further research is required to study the pullout bearing capacity in layered clay under a wide range of L/D.

(3) This paper only systematically analyzes the pullout capacity of suction anchors and does not examine in detail how pullout capacity varies with different padeyes. Further research is needed to analyze the pullout capacity of suction anchors at different padeyes under inclined loads.

7. Conclusion

In this paper, finite element numerical simulations were conducted to investigate the pullout capacity of suction anchors under inclined loading. The study examined the effects of varying soil parameters S_um/kD and layer thickness ratios T_h/L and T_c/L on the optimum centerline loading depth Z_cl,opt, uniaxial ultimate bearing capacity (H_ult and V_ult), and VH failure envelope of suction anchors. The main conclusions are as follows.

(1) The layer thickness ratio Th/L of the LOC is the core factor controlling Z_cl,opt. When the suction anchor partially embeds into the LOC layer, as the LOC layer thickness increases, the upper portion of the suction anchor contributes more lateral friction resistance. To balance the moments, the suction anchors Z_cl,opt shifts downward. When the suction anchor is fully embedded in the LOC soil layer, Z_cl,opt depends solely on the parameters of the LOC soil mass.

(2) The thickness ratio T_h/L of the LOC soil layer is the core factor determining the uniaxial ultimate bearing capacity, as the undrained shear strength of the LOC layer is significantly higher than that of the NC layer.With the larger the length of the suction anchor embedded into the LOC layer, the combined effect of lateral friction resistance and end resistance at the anchor-soil interface becomes more pronounced. The layered soil makes the soil strength along the clay surface exhibit non-uniformity. The increase in S_um/kD enhances the uniaxial bearing capacity of layered clay. The thickness of the NC layer T_c/L, has little effect on the horizontal ultimate bearing capacity H_ult but significantly influences the vertical ultimate bearing capacity V_ult.

(3) The strength of the LOC soil is relatively high. As the thickness ratio T_h/L increases, the contact area between the suction anchor and the high-strength clay expands, which provides more uniform and greater lateral friction at the anchor-soil interface. When the horizontal bearing capacity increases, the decrease in vertical bearing capacity becomes more gradual. The lower bound solution for the VH failure envelope of the corresponding layered soil has been obtained.

(4) A comparative analysis of the pullout bearing capacity of suction anchors under inclined loads in layered clay and homogeneous clay reveals that soil layering significantly affects the bearing capacity of suction anchors. In homogeneous clay, the pullout bearing capacity of suction anchors is solely influenced by the soil strength S_um/kD. In contrast, the pullout bearing capacity in layered clay is affected by factors including the soil strength S_um/kD and the soil layer thickness ratio T/L.