Abstract

This paper presents the results from numerical modeling of the keying process of plate anchors in three-layer soft-stiff-soft clay deposits. Three-dimensional large deformation finite-element analyses were carried out, and the results were firstly validated by the centrifuge test data and the previous numerical results. The soil flow mechanism during the keying process of plate anchors was examined, and a series of parametric studies were performed to investigate the factors affecting the rotation characteristics of plate anchors with an emphasis on the presence of the interbedded stiff soil layer. The results indicate that the loss of embedment depth of plate anchors decreases with the increase of the thickness of the first soil layer when the anchor is initially located at the middle of stiff soil layer. The flow velocity of soil around the anchor that is initially embedded at the first layer and adjacent to the underlying interbedded stiff soil layer is generally larger, resulting in a smaller embedment depth loss compared with the traditional normally consolidated soil layer. The interbedded stiff soil layer affects the keying process of plate anchors embedded 1.0B above and 2.0B below the interbedded stiff soil layer (B is the width of the square plate anchors). The increase of the strength of local soil around the plate anchors leads to the increase of the embedment depth loss, but the increase of the strength of soil slightly away from the plate anchors leads to the decrease of the embedment depth loss.

1. Introduction

Suction embedded plate anchor is a new type of foundation for marine platform structures. It has the advantages of accurate positioning, low cost, easy operation, and high vertical uplift bearing capacity. Thus, it has been extensively used in offshore engineering and has been used as a type of permanent anchor foundation in Angola [1].

When the plate anchor installation is completed by using the suction caisson, the anchor is generally embedded vertically in the seabed, as shown in Figure 1(a) (Zi is the initial embedment depth of the plate anchor, B is the width and t is the thickness of the plate anchor, and e is the anchor padeye eccentricity). Then, the plate anchor is rotated by pulling the mooring line that is attached to the plate anchor (Figures 1(b) and 1(c)), and the plate anchor is finally rotated from vertical to horizontal direction to serve as a foundation (Figure 1(d)). This process is termed “keying” process. During the keying process, the embedment depth of plate anchor will be greatly reduced, which may lead to the decrease of the bearing capacity of anchor. Therefore, the variation of embedment depth and the bearing capacity of plate anchors in the keying process are critical issues in the design of the plate anchor.

The bearing capacity of foundation is an important problem that has been widely studied (e.g., [2, 3]). Similarly, the pullout capacity of the plate anchors is also a critical problem in ocean engineering. Das [4] and Das and Puri [5] proposed a calculation program to estimate the ultimate pullout capacity of the plate anchors in clay by using the small-scale model test under 1 g condition. Merifield et al. [6, 7] presented rigorous numerical limit analysis for the ultimate capacity of horizontal, vertical, and inclined strip plate anchors in both homogeneous and inhomogeneous clay soils. Similarly, by using the finite-element limit analysis, the ultimate uplift capacity of multiplate strip anchors and strip anchors embedded adjacent to sloping ground in fully cohesive soil under undrained condition have also been studied [811]. These complex interactions between adjacent plates have also been presented in the seismic problems [12]. Zhang et al. [13] investigated the undrained behavior of vertical uplift plate anchors in remolded soft Ariake clay by laboratory tests and finite-element analyses.

Based on the field experiment and laboratory test, Wilde et al. [14] and Wilde [15] derived the load-displacement curve of the keying process of plate anchors and examined the relation between the ultimate pullout capacity of the plate anchor and the soil properties, disturbances, and strain rate. Gaudin et al. [1618] studied the suction caisson installation, recovery of the suction caisson, and the keying process of plate anchors through a series of centrifuge test. The experimental results were compared with that of the plate anchors directly embedded in the soil before the test, and the effects of the installation of the suction caisson, the recovery of the caisson, and the keying process of plate anchors were studied. O’Loughlin et al. [19] conducted centrifuge tests by using a semimodel to investigate the loss in plate anchor embedment during the keying process. The results show that the keying process will result in the decrease of the embedment depth of plate anchor, and anchor padeye eccentricity has great influence on the displacement of anchor. Song et al. [20, 21] investigated the pullout capacity and rotation behaviour of the plate anchors which were initially installed vertically in homogeneous clay based on model test under 1 and 100 g conditions and numerical analysis and proposed an influence factor on the embedment loss of the vertical pulling plate anchors. Yu et al. [22] discussed the influence of anchor padeye eccentricity, pullout angle, and shape of anchor on the embedment loss of plate anchors during keying by numerical simulation. It is concluded that the influence of anchor padeye eccentricity is the greatest and the shape of anchor is the least. Wang et al. [2325] analyzed the keying process of plate anchors embedded in the normally consolidated clays based on the remeshing and interpolation technique with small strain (RITSS), and they found that the losses of the embedment depth of the plate anchors were mainly dependent on the coupling effect of the anchors roughness coefficient, thickness of plate anchors, and anchor padeye eccentricity. Some researchers (e.g., [26]) also studied the effects of the angle of anchor-chain inclination on reducing the embedment depth loss of the plate anchors. Their results showed that with the decrease of the angle of anchor-chain inclination, the depth of the anchors being pulled out decreases linearly. Liu et al. [27] proposed a plasticity model of three-dimensional plate anchors with consideration of shank and flap. The effects of the flap on the embedment loss and pullout capacity, and the effects of loading eccentricity and pullout angle on diving behavior were investigated. Based on discrete-element method (DEM) simulations, Evans and Zhang [28] investigated the behavior of plate anchors during pullout from the view of microscale physical processes. The observations of the simulated contact force network and particle velocity during uplift were used to provide insight into anchor failure mechanisms. In addition, Singh and Ramaswamy [29, 30] performed model test of plate anchors under cyclic loadings and investigated the effects of the dynamic loading level, the frequency of cyclic loading, and the preloading level on the embedment depth loss. Rahimi et al. [31] performed a series of near-full-scale experiments and numerical models on plate anchors embedded in sand with and without geocell reinforcement, and the influence of geocell reinforcement on the uplift capacity of plate anchors was examined.

In summary, current studies on the keying process of plate anchors mainly focused on the homogeneous soil or the normally consolidated soil. However, the deep-sea soil layers are quite complex. In emerging provinces and fields, highly layered soils are prevalent and quite a number of stratified seabed profiles involve interbedded layers displaying strong variations in shear strength [32]. Therefore, this study utilized the three-dimensional Coupled Eulerian Lagrangian (CEL) large deformation finite element analysis to investigate the keying process of plate anchors on the more complicated layered soil, i.e., the three-layer soft-stiff-soft clay deposits. The flow of soil and the embedment depth loss of plate anchors in the keying process were analyzed, and the effects of soil strengths of different layers were investigated.

2. Validation of the Numerical Model

To validate the reliability of CEL method in predicting the interaction between plate anchors and soil, the loss of the embedment depth during keying process obtained by numerical calculation was compared with that from the existing centrifuge model test [20]. The numerical analyses were performed using the FE package ABAQUS/Explicit [33]. The feature of fixed Eulerian mesh and the contact between the Eulerian and Lagrangian materials allows for effective application of the CEL approach in large deformation geotechnical problems. The soil was assumed to be completely saturated and elastoplastic, and the Mohr–Coulomb model was adopted.

The centrifuge tests [20] were carried out by using the drum centrifuge which has a diameter of 1.2 m, a channel of radial depth of 0.2 m, and a height of 0.3 m to study the relation between the anchor padeye eccentricity of pulling force and the loss of the embedment depth. A uniform clay soil deposit was considered with an undrained shear strength, Su = 18 kPa, and a submerged unit weight of 9.5 kN/m3. The stainless-steel square plate anchor in the model test has a thickness of 2 mm, a length of 40 mm, and an eccentricity of 25 mm. The test was conducted at 100 g to model a prototype plate anchor of 4 m × 4 m × 0.2 m, and an eccentricity of 2.5 m. The submerged unit weight of the anchor is 17.0 kN/m3. The anchor was installed into the soil vertically and was pulled out vertically. The normalized initial embedment of the anchor is Zi/B = 3.

Following the previous studies, the loss of the embedment in this paper refers to the centroid displacement of plate anchor (ΔZ). The normalized pulling force is defined as F/SuA, where F is the vertical pulling force on the plate anchor, and A is the plate area. Figure 2(a) shows the variation of normalized pulling force with the normalized displacement of the anchor padeye (d/B) in this study. Also shown in the figure are the results from the model test [20] and the RITSS method [23]. It can be seen that the normalized pulling forces obtained from CEL technique are between the results of RITSS method and the results of model test. Figure 2(b) presents the variation of the inclination of the plate anchor (β) with the normalized loss of the embedment during keying process of plate anchors (ΔZ/B). The loss of the embedment depth of plate anchors from two types of numerical simulation is larger than the results from model tests. This may be related to the neglection in the numerical model of both the anchor shank and the embedded disturbances to soil of the plate anchor in the model test. In addition, the model anchor reached its stable state at an inclination of 10° from the horizontal direction. If the final orientation angle of the anchor is also taken as 10° in the numerical analyses, the loss of embedment would be much similar to the model test.

The second validation is for the model test results from O’Loughlin et al. [19], in which the keying process of the plate anchors in a normally consolidated clay was modeled by using the drum centrifuge. The undrained shear strength of the soil increases linearly with depth. The shear strength of soil on the surface of seabed is zero, and the rate of increase of shear strength is 0.7 kPa/m. The submerged unit weight of soil and plate anchor is 6.5 and 67 kN/m3, respectively. The width of the plate anchor is B = 3 m, the ratio of thickness to width is t/B = 0.067, and the eccentricity ratio of plate anchor e/B is 0.50. The normalized initial embedment of the anchor is Zi/B = 3. The loss of the embedment depth of plate anchors from CEL method and model test is shown in Figure 3. It can be seen that the numerical results in this study are generally close to the test, except for the stages when the keying process is nearly finished. Similar to the results of uniform soil, this is due to that the model anchor in the test reached its stable state with an inclination and not in the horizontal direction, resulting in some differences in the embedment loss compared to the numerical results.

3. Model Parameters

The soil in the model consists of three layers, i.e., the top and bottom layers of normally consolidated soil and the interbedded stiff layer of uniform soil (Figure 4). The undrained shear strength of the top and bottom soil layers increases linearly with depth as Su = (1 + kz) kPa, where k is the rate of increase of shear strength (kPa/m) and z is the depth (m). The undrained shear strength of the uniform stiff soil is Su = 40 kPa. The effective unit weight for soils is 6 kN/m3. Young’s modulus (E) is taken as E/Su = 500 and a Poisson’s ratio of 0.49 is adopted to give minimal volumetric strains, while maintaining numerical stability [32]. The plate anchor is a rigid body, and the submerged unit weight is 67.0 kN/m3.

The three-dimensional (3D) elements were used in the CEL analysis. For the parametric studies, only a half of the domain was involved in the analyses due to the symmetry inherent. The roughness factor of soil-anchor interface was taken as 0.3. The soil and plate anchor were discretized using Eulerian hexahedral elements of type EC3D8R and continuum elements of types C3D8R (hexahedron), respectively. The soil layers underlay an initially void layer to accommodate the possible soil heave resulting from keying process of the plate anchor. This void layer has no weight and strength, but it can provide an additional domain to accommodate the soil flow during the large deformations. The length, width, and thickness of the plate anchor are L = 4 m, B = 4 m, and t = 0.2 m, respectively. The numerical model of the plate anchor and soil profiles is shown in Figure 5. The soil domain is chosen as 10B in length and 5B in width, which is shown to be sufficiently large to avoid boundary effect. Similar domain sizes were considered by Zheng et al. [34] in spudcan continuous penetration analysis. For the selection of element size, routine preliminary simulations are carried out, varying the minimum element size around the plate anchor, hmin/B. Balancing the computational effort and accuracy of the computed results, and also complying with existing studies (e.g., Zheng et al. [34]), hmin/B = 0.025 is adopted.

4. Soil Flow Mechanism

The keying process of the plate anchors is a process combining translation and rotation and involves problems of highly nonlinear and large deformation. The soil generally flows faster around plate anchors during keying rotation process and the flow of soil has great influence on the loss of the embedment depth and bearing capacity of the plate anchors. In this section, the soil flow mechanism of plate anchors with various initial embedment depths in normally consolidated soil and three-layer soft-stiff-soft clay deposits was studied.

Figure 6 shows the soil flow around plate anchors in normally consolidated and layered soil for an initial embedment depth of Zi = 3.5B (i.e., 14 m). The plate anchor is thus located at the bottom of first soil layer and is close to the middle stiff soil layer for the three-layer deposits. It can be found that the flow of soil around the rotating plate anchor in the three-layer deposits is faster than the normally consolidated soil as a result of the presence of the underlying stiff layer. The backflow soils strike the plate of anchors, and a large velocity difference is formed on the upper and lower faces of the plate anchor. This accelerates the keying process of the plate anchors, which is conducive to reduce the loss of the embedment depth.

Figure 7 shows the soil flow around plate anchors in normally consolidated and layered soil for the initial embedment depth of Zi = 7.0B (i.e., 28 m) when the keying process is completed (i.e., the plate anchor is horizontal). In this case, the plate anchor is buried in the third soil layer (4 m away from the soil interface), the upward flow for layered soil is slower than that in the normally consolidated soil, which slightly reduces the loss of the embedment depth of plate anchors.

5. Parametric Study

5.1. Effect of the Thickness of Top Soil Layer

The influence of soil thickness of the top soil layer on the rotation characteristics of the plate anchor was analyzed. The thickness of the top soil layer H1 was selected as 8 m, 12 m, 16 m, 20 m, and 24 m, respectively. Considering that the thickness of the top soil layer may have different effects on plate anchors at different initial embedded depths, three typical initial embedment positions were selected, as shown in Figure 8. Among them, P1 is located at the interface of first and second soil layers, P2 is located at the middle of second soil layer, and P3 is located at the interface of second and third soil layers. In order to study the rotation characteristics of plate anchors at different positions, the variations of the inclination of the anchor plate with the normalized loss of the embedment depth (ΔZ/B) are examined.

Figure 9(a) shows the variation of plate inclinations with the loss of the embedment depth of plate anchors at position P2 for different thicknesses of the first soil layer (H1), and Figure 9(b) shows the ultimate loss of the embedment depth at the end of keying process of plate anchor (ΔZu). We can see the loss of the embedment depth decreases with the increase of the thickness of the first soil layer. The maximum embedment loss is about 1.10B and 0.75B for the thickness of the first soil layer of H1 = 8 m and H1 = 24 m, respectively. When the thickness of the first soil layer increases to 16 m, its influence on the rotation characteristics of the plate anchor is limited.

Figure 10 shows the variation of plate inclinations with the loss of the embedment depth of plate anchors at positions P1 and P3 during keying process for different H1 values. Observe that the effect of the thickness of the first soil layer on the rotation process of plate anchor is negligible when the plate anchor is initially buried at the interface of first and second soil layers or at the interface of second and third soil layers. In addition, the ultimate loss of the embedment depth is similar for the two positions (about 0.60B for the plate anchor at position P1 and 0.70B for the plate anchor at position P3).

5.2. Effect of Initial Embedment Depth

Figure 11 shows the variation of plate inclinations with the loss of the embedment depth for plate anchors initially buried in the first, second, and third soil layers, respectively. Observe that when the plate anchor is initially buried in the first soil layer, the embedment loss of the plate anchor increases with the increase of the initial embedment depth Zi. The embedment loss is minimum for Zi = 8 m (about 0.37B), and the maximum embedment loss is about 0.57B for Zi = 16 m. When the plate anchor is buried in the second soil layer, the embedment loss increases initially and then decreases with the increase of the initial embedment depth of plate anchors. When the plate anchor is located in the third soil layer, the influence of the initial embedment depth is negligible, with the maximum embedment loss of the plate anchor of about 0.7B.

Figure 12(a) shows the ultimate loss of embedment depth of plate anchors at different initial embedment depths in layered and normally consolidated soils. It can be seen that the influence of the interbedded uniform stiff soil layer on the embedment depth loss is significant for plate anchors embedded at 12 m∼32 m (i.e., from 1.0B above to 2.0B below the stiff soil layer). When the plate anchor is initially buried at other depths, the influence of the interbedded stiff soil layer is almost negligible. The embedment depth loss of plate anchors in the three-layer soft-stiff-soft and normally consolidated soil can be normalized aswhere and are the ultimate embedment depth loss of plate anchors in layered soil and normally consolidated soil, respectively. The influence of the interbedded stiff soil layer on the embedment loss of plate anchor is further shown in Figure 12(b). This result is useful for the further estimation on the variation of the uplift bearing capacity of plate anchors.

5.3. Effect of the Shear Strength of Layered Soil

The keying process of plate anchors initially embedded in the middle of interbedded uniform stiff soil layer was investigated for soil layers with different shear strengths. The first case is for different undrained shear strengths of the stiff layer Su2 = 30, 35, 40, 45, and 50 kPa and the strength gradient of the normally consolidated layers of 1.0 kPa/m. The other case is for Su2 = 40 kPa and the strength gradient of the normally consolidated layers of k = 0.5, 0.8, 1.0, 1.5 and 2.0 kPa/m, respectively. The results for different cases are shown in Figure 13.

For the first case (Figure 13(a)), the embedment depth loss of plate anchors increases significantly with the increase of the Su2. The minimum ultimate embedment depth loss is about 0.70B for Su2 = 30 kPa, and the maximum ultimate embedment depth loss is about 1.0B for Su2 = 50 kPa. For the second case (Figure 13(b)), the embedment depth loss of plate anchors decreases with the increase of the strength gradient (k) of the normally consolidated layer, and it is essentially in the range of 0.7B and 0.9B. Note that this is different from the observation of Wang et al. [23], in which they found that the embedment depth loss increases with the increase of k in normally consolidated soil. Therefore, it is inferred that the increase of the strength of local soil around the plate anchors will lead to the increase of the embedment depth loss, but the increase of the strength of soil slightly away from the plate anchors may lead to the decrease of the embedment depth loss.

6. Conclusions

In this study, the keying process of plate anchors in three-layer soft-stiff-soft clay deposits was investigated by using the CEL approach to perform the large deformation numerical simulation. The soil flow mechanism around the plate anchors and the influence of a series of parameters of layered soil on the embedment depth loss were examined. The results were compared with those from the traditional results of normally consolidated soil. The key conclusions are as follows:(1)The influence of the thickness of the top soil layer is significant on the keying process of plate anchors buried at the middle of the interbedded stiff soil layer. With the increase of the thickness of the top soil layer, the loss of embedment depth of plate anchors decreases. The influence of the thickness of the top soil layer on the rotation process of plate anchor is negligible when the plate anchor is initially buried at the interface of top and stiff soil layers or at the interface of stiff and bottom soil layers.(2)When the plate anchors are buried at the top soil layer and near the stiff soil layer, the surrounding soil flow velocity is larger and the embedment depth loss of plate anchors is smaller. Compared with the anchors buried in the case of normally consolidated soils, the interbedded stiff soil layer affects the keying process of plate anchors embedded in 1.0B above and 2.0B below the interbedded stiff soil layer.(3)The strength gradient of the normally consolidated soil layers and strength of the uniform stiff soil layer have significant effects on the keying process of plate anchors. With the increase of local soil strength around the plate anchor, the loss of embedment depth increases, but the loss of embedment depth decreases with the increase of strength gradient of the normally consolidated soils, which is different from previous observations on the normally consolidated soils.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research has been supported by the Fundamental Research Funding Scheme for the Central Universities in China (grant nos. 2018B42414, 2016B01114, and 2018B57014), the National Natural Science Foundation of China (grant nos. 41602280, 51808193, and 51778211), Jiangsu Province Postdoctoral Science Foundation (grant no. 2018K044B), and China Postdoctoral Science Foundation (grant no. 2016M601708). These supports are gratefully acknowledged.