​

​

==Numerical investigation of setup effects in the vertical holding capacity of finless torpedo anchors==

​

Guilherme K. Lopes<sup>a,*</sup>, José Renato M. de Sousa<sup>a</sup>, and Gilberto B. Ellwanger<sup>a</sup>

​

<sup>a </sup>Civil Engineering Department – PEC, Federal University of Rio de Janeiro, COPPE/UFRJ, 21945-970, Rio de Janeiro, RJ, Brazil.

​

*e-mail: [mailto:kronemberger.guilherme@laceo.coppe.ufrj.br kronemberger.guilherme@laceo.coppe.ufrj.br], Tel.: +55 24992272120

-->

​

==Abstract==

​

This article presents a numerical-based study on the vertical holding capacity of a typical finless torpedo anchor embedded in a pure cohesive isotropic sol considering setup effects. A strategy based on two axisymmetric nonlinear finite elements (FE) models is proposed. In these models, the soil is represented with continuous solid elements with both plane translational and pore pressure degrees of freedom, which allow the modeling of the soil's strength regain after the installation of the anchor, ''i. e.'', the soil setup. The anchor is also modeled with solid elements, and its interaction with the soil is addressed with contact elements. Several analyses are conducted to evaluate setup effects, and, additionally, a parametric study on different soil undrained shear strengths and permeability coefficients is performed. The results obtained show that the permeability of the soil and its undrained shear strength have an important role in the setup process, and an inverse relationship between consolidation time and the permeability coefficient of the soil is observed. Finally, an analytical approach is calibrated with a reasonable agreement to the numerical results.

​

'''Keywords''': Torpedo anchor, setup, finite elements, analytical models

​

==1. Introduction==

​

Dynamically installed anchors (DIAs) <span id='cite-_Ref31803537'></span>[[#_Ref31803537|[1]]], also known as gravity installed anchors (GIAs) <span id='cite-_Ref31803040'></span>[[#_Ref31803040|[2]]], are an attractive alternative for deep-water mooring systems due to their cost-effectiveness and capacity to withstand high vertical loads. Existing DIAs used in offshore engineering include deep penetration anchors (DPAs) <span id='cite-_Ref31817098'></span>[[#_Ref31817098|[3]]], torpedo anchors <span id='cite-_Ref31804075'></span>[[#_Ref31804075|[4]]], the OMNI-Max anchor <span id='cite-_Ref31804736'></span>[[#_Ref31804736|[5]]], and, more recently, an innovative lightweight gravity installed plate anchor (L-GIPLA) <span id='cite-_Ref31808849'></span>[[#_Ref31808849|[6]]]. Among these anchors, torpedo anchors have proven to be an outstanding alternative in Brazilian offshore fields.

​

As shown in <span id='cite-_Ref52370445'></span>[[#_Ref52370445|Fig. 1]], the torpedo anchor has a “rocket” shape with a long shaft made of carbon steel ballasted with lead and cast iron <span id='cite-_Ref31804075'></span>[[#_Ref31804075|[4]]]. Its structure typically comprises a loading point (or padeye) attached at the top, a conical tip, which is designed to improve the embedment depth of the anchor and may have up to 4 relatively small fins at the trailing edge <span id='cite-_Ref31810175'></span>[[#_Ref31810175|[7]]].

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image1.png|600px]] </div>

​

<div id="_Ref52370445" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 1 - Typical torpedo anchor: (a) conical tip and (b) padeye details ''<span id='cite-_Ref31804075'></span>[[#_Ref31804075|''[4]'']]''.''</div>

​

In general, a set of four to eight anchors can be transported to a location and deployed using a typical anchor-handling vessel (AHV). A schematic of the installation of a typical torpedo anchor is presented in <span id='cite-_Ref52370458'></span>[[#_Ref52370458|Fig.  2]]. The dropping anchor is connected to an installation line and an independent mooring line. The anchor is lowered to the predetermined drop height above the seabed and then released by disconnecting the installation line <span id='cite-_Ref31892302'></span>[[#_Ref31892302|[8]]]. The anchor falls freely through the water depth using its weight as driving energy before impacting and embedding into the seafloor sediments.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image2.png|384px]] </div>

​

<div id="_Ref52370458" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  2 – Typical installation procedure of a finless torpedo anchor.''</div>

​

Soil disturbance caused by the anchor installation induces an excess pore pressure around the anchor, which reduces its holding capacity. As the excess pore pressure dissipates over time, the holding capacity of the anchor progressively increases. This process is known as setup.

​

Setup is predominantly associated with an increase in the soil resistance around the anchor. Although the exact mechanism by which setup occurs is not entirely understood <span id='cite-_Ref502687747'></span>[[#_Ref502687747|[9]]]-<span id='cite-_Ref502153701'></span>[[#_Ref502153701|[11]]], the most important process is due to an increase in soil effective stresses with the dissipation of the excess pore pressure built up around the driven anchor and the gradual consolidation of the remolded soil. In addition, thixotropy and aging effects may also account for additional setup.

​

In general, the setup process can be divided into three main phases <span id='cite-_Ref502153693'></span>[[#_Ref502153693|[10]]]. <span id='cite-_Ref52286155'></span>[[#_Ref52286155|Fig.  3]] presents a typical setup curve, which indicates the variation of the anchor's normalized pullout (vertical) capacities with the logarithm (''log'') of time. In this figure, the normalized capacity is the ratio between the pullout capacity at time ''t'' and the capacity immediately after installing the anchor. During the first phase of setup, the rate of excess pore pressure dissipation is not linear with the ''log'' of time due to the highly disturbed state of the soil. According to Komurka and Wagner [10], in this phase, the effective stress increase and, therefore, the gain in soil strength is not well-understood and challenging to represent. In the second phase, the excess pore pressure dissipation rate becomes linear with the ''log'' of time, and the increase in the effective stresses follows the conventional consolidation theory. Theoretically, infinite time is required to dissipate the excess pore water pressure, but, in practice, after some time, the rate of dissipation is so slow that almost no strength regain is observed. This response corresponds to the third phase in Fig. 3. In this phase, primary consolidation is completed, but secondary compression continues and is independent of the effective stress [10].

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image3.png|600px]] </div>

​

<div id="_Ref52286155" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  3  – Example of a setup curve.''</div>

​

Setup mechanisms have an important role in predicting the ultimate holding capacities of DIAs. Richardson ''et al''. <span id='cite-_Ref31897166'></span>[[#_Ref31897166|[12]]] conducted several experimental tests, which indicated that 50% of the anchors’ holding capacities are achieved from 35 to 350 days after installation. Furthermore, 90% of the holding capacity is achieved between 2.4 to 24 years after installation but these values highly depend on the soil properties. Lieng ''et al. ''<span id='cite-_Ref31817098'></span>[[#_Ref31817098|[3]]]'' ''performed experimental tests and found that 70% of the holding capacity of a DPA is achieved 2 weeks after installation in the clay soil of Voring Bay in Norway. Radgahar ''et al.''<span id='cite-_Ref31899459'></span>[[#_Ref31899459|[13]]] and Raie and Tassoulas <span id='cite-_Ref31899465'></span>[[#_Ref31899465|[14]]] conducted finite element analyses of finless torpedo anchors. The results indicated that the time needed to achieve a degree of consolidation around 90% is highly inferior to the time needed for the complete consolidation of the soil. It is also emphasized that a significant regain of the anchor’s holding capacity was observed after its installation. Sabetamal ''et al''. <span id='cite-_Ref31900575'></span>[[#_Ref31900575|[15]]] focus on numerical analyses with the Arbitrary Lagrangian-Eulerian (ALE) approach to evaluate DIAs holding capacities simultaneously with pore pressure dissipation and soil consolidation with time. The authors concluded that most of the pullout capacity of the anchors should be available much earlier than the consolidation completion, as most of the excess pore pressure dissipates within a matter of days or weeks.

​

More recently, Zhao ''et al.'' <span id='cite-_Ref31810175'></span>[[#_Ref31810175|[7]]] presented the available numerical methods to analyze dynamically installed anchors and emphasized that an integrated and complete simulation of installation and pullout of a DIA is still a challenge. Regarding the setup analysis of finless torpedoes, a degree of consolidation of 90% for elements within a radial distance of 1 ''D<sub>a</sub>'' of the anchor tip, where ''D<sub>a</sub>'' is the anchor’s external diameter, was achieved within 51-82 days after the installation.

​

Despite its successful use in offshore applications, estimating a torpedo anchor holding capacity and predicting the stresses in its structure remains a challenge <span id='cite-_Ref31804075'></span>[[#_Ref31804075|[4]]]. Moreover, setup effects are not directly incorporated into the design of this structure. Sousa ''et al''. <span id='cite-_Ref31804075'></span>[[#_Ref31804075|[4]]] mention that, in practice, a typical torpedo anchor is first loaded approximately three months after its installation. However, the imposed loads are much lower than its holding capacity, as safety factors between 1.5 and 2.0 are employed in the anchor’s design. Nevertheless, there is no guarantee that this time is sufficient to achieve the required safety factor. On the other hand, a possible reduction in the time needed to connect the floating production system and the mooring line would be economically beneficial. Hence, there is still a gap of information concerning the setup process in the design of torpedo anchors, which is important to the offshore oil and gas industry.

​

In this work, a nonlinear axisymmetric finite element (FE) approach is proposed to predict the axial holding capacity of finless torpedo anchors considering setup effects. According to Medeiros Jr. <span id='cite-_Ref32417380'></span>[[#_Ref32417380|[16]]], finless torpedo anchors may be used to anchor flexible lines to avoid load transfer from flowlines to submarine equipment. Moreover, monobuoys and ships may also be anchored in shallow waters using this type of torpedo.

​

The proposed approach employs a FE model in which isoparametric solid elements are employed to represent both the soil and the anchor. These elements can describe the nonlinear physical behavior of the soil and, eventually, of the anchor. Large deformations are also addressed. Soil-anchor interaction is ensured by surface-to-surface contact elements placed on the external surface of the anchor and the surrounding soil. It is assumed that the anchor is “wished in place” and that soil initial stress conditions immediately after installation are obtained by applying the Cavity Expansion Method (CEM) in the proposed FE model. The initial stress conditions obtained are then input in another FE model, where the associated pullout capacity is determined. Both models share the same FE mesh to simplify the transference of the stress state to the holding capacity analysis. Moreover, this approach was implemented in Abaqus [[#_Ref534978548|[17]]], and a view of a typical FE mesh is presented in <span id='cite-_Ref52370499'></span>[[#_Ref52370499|Fig.  4]].

​

This approach was employed to analyze a finless torpedo anchor embedded in a purely cohesive isotropic soil. Different soil's undrained shear strengths and permeability coefficients were considered, and their impact on the anchor pullout capacity was assessed.<span id='cite-_Ref534978548'></span> Moreover, the pullout capacities obtained with the FE model were compared to those estimated using some existing analytical approaches, such as the ones of Skov and Denver <span id='cite-_Ref534917435'></span>[[#_Ref534917435|[18]]], Svinkin and Skov <span id='cite-_Ref534917465'></span>[[#_Ref534917465|[19]]], and ABS <span id='cite-_Ref51515640'></span>[[#_Ref51515640|[20]]]. Additionally, in this work, an analytical methodology based on the α-Method from API <span id='cite-_Ref502324791'></span>[[#_Ref502324791|[21]]]. is modified to account for setup effects. By relying on these comparisons, an analytical expression was calibrated to obtain the variation of the pullout capacity over time.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image4.png|510px]] </div>

​

<div id="_Ref52370499" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  4 – General view of the FE model.''</div>

​

==2. Finite Element Approach==

​

==2.1. Soil Modeling==

​

The soil is modeled in Abaqus <span id='cite-_Ref534978548'></span>[[#_Ref534978548|[17]]] by considering the medium as a multiphase material and adopting the effective stress principle to describe its behavior. The FE mesh is attached to the solid phase, and the liquid present in the interior of the soil can flow through this mesh. A continuity equation is, therefore, required for the liquid. This equation can be obtained by equating the rate of increase in liquid mass stored at a point to the rate of the mass of liquid flowing into this point in a time increment:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\int_{V}^{\, }\frac{1}{J}\frac{d}{dt}\left( J{\rho }_{w}{n}_{w}\right) dV=-\int_{S}^{\, }{\rho }_{w}{n}_{w}n\cdot {v}_{w}dS</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(1)

|}

​

​

where <math display="inline">J</math> is the ratio of the medium’s volume in the current configuration to its volume in the reference configuration (<math display="inline">J \overset{\underset{\mathrm{def}}{}}{=} \left| \frac{dV}{{dV}^{0}}\right| </math>), <math display="inline">{\rho }_{w}</math> is the density of the wetting liquid, <math display="inline">{n}_{w}</math> is the volume ratio of free wetting liquid at a point, <math display="inline">n</math>''' '''is the outward normal to the surface <math display="inline">S</math>, and <math display="inline">{v}_{w}</math> is the seepage velocity.

​

The constitutive behavior for pore fluid flow is assumed to be governed by Darcy’s law, which states that, under uniform conditions, the volumetric flow rate of the wetting liquid through a unit area of the medium is proportional to the negative of the piezometric head, ''i. e.'':

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>sn{v}_{w}=-k\cdot \frac{\partial \phi }{\partial x}</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(2)

|}

​

​

where <math display="inline">s</math> is the soil saturation, <math display="inline">n</math> is the porosity, <math display="inline">k</math>''' '''is the permeability of the medium, and <math display="inline">\phi</math>  is the piezometric head, defined as:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\phi \overset{\underset{\mathrm{def}}{}}{=}z+\frac{{u}_{w}}{g{\rho }_{w}}</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(3)

|}

​

​

where <math display="inline">z</math> is the elevation above a datum, <math display="inline">{u}_{w}</math> is the pore pressure, and <math display="inline">g</math> is the magnitude of the gravitational acceleration.

​

The soil is assumed to be a perfectly elastoplastic isotropic material with physical properties varying with depth. The constitutive model relies on simple bulk elasticity relationships (generalized Hook’s law) for the liquid and soil grains, combined with a constitutive theory for the soil skeleton.

​

The Mohr-Coulomb criterion was chosen to represent the failure of the soil. This criterion assumes that the maximum shear stress controls the failure, ''i. e.'':

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\tau ={c}^{'}-{\sigma }^{'}.tan\phi '</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='Eq_4'></span>(4)

|}

​

​

where <math display="inline">\tau</math>  is the shear stress, <math display="inline">{c}^{'}</math> is the effective cohesion of the soil, <math display="inline">{\sigma }^{'}</math> is the effective normal stress, and <math display="inline">{\phi }^{'}</math> is the effective internal friction angle of the soil.

​

An overview of the main dimensions of the soil’s mesh is shown in <span id='cite-_Ref52370499'></span>[[#_Ref52370499|Fig.  4]]. The proposed mesh is a cylinder with a base diameter of <math display="inline">20{D}_{a}</math>. The height of the cylinder is given by the sum of the penetration of the torpedo anchor, <math display="inline">{H}_{p}</math>, the length of the anchor, <math display="inline">{L}_{a}</math>, and the distance of the tip of the torpedo anchor to the bottom of the FE mesh, <math display="inline">{H}_{a}</math>.

​

The elements have dimensions varying between 0.10 m and 0.25 m in the regions where high plastic strains are expected to occur (close to the anchor) and between 0.25 m and 0.50 m in the regions far from the anchor, as shown in <span id='cite-_Ref52370499'></span>[[#_Ref52370499|Fig.  4]]. These dimensions were adopted after a mesh convergence study, where three different mesh densities (dense, intermediate, and coarse) were considered. The coarse mesh had about half the elements used in the intermediate mesh, while the number of elements in the dense mesh was twice the number in the intermediate mesh. As the intermediate and dense meshes did not present significant differences concerning the magnitude of the displacements and stresses, the intermediate mesh was adopted in this study.

​

The vertical and horizontal displacements of the cylinder were restrained at the nodes of its base. The displacements in the radial direction of the nodes associated with the outer wall of the cylinder were restrained. Moreover, once a plane of symmetry is found, the out-of-plane displacements are also restrained. Sousa ''et al''. <span id='cite-_Ref31804075'></span>[[#_Ref31804075|[4]]] performed several mesh tests to avoid any influence of the boundary conditions on the response of the anchor and concluded that a diameter of <math display="inline">20{D}_{a}</math> for the cylinder and a distance of 5.0m for <math display="inline">{H}_{a}</math> was enough to simulate an “infinite” media. Thus, these values were also considered in all analyses performed in this study.

​

==2.2. Anchor Modeling==

​

The finless torpedo anchor was modeled with isoparametric solid elements analogous to those used in soil representation but without the pore pressure degree of freedom. These elements, as shown before, are capable of considering both material and geometric nonlinearities.

​

It is worth mentioning that neither the padeye at the top of the anchor nor the mooring line were represented in the proposed model. Hence, the load from the mooring line was applied at a reference node placed 1.0 m above the anchor, and it was rigidly connected to the top of the anchor by rigid bars using the beam MPC (Multiple Point Constraint) in Abaqus <span id='cite-_Ref534978548'></span>[[#_Ref534978548|[17]]], as presented in <span id='cite-_Ref52370593'></span>[[#_Ref52370593|Fig.  5]].

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image5.png|216px]] </div>

​

<div id="_Ref52370593" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  5 – Detail of the point of load application.''</div>

​

Anchor loading was simulated by applying a concentrated force at the reference node following a simple ramp function. As the focus of this work is evaluating the undrained capacity of torpedo anchors, the load was applied at a very high rate. Hence, after a few trials, a value of 10 000 kN/s was adopted.

​

==2.3. Anchor–Soil Interaction==

​

In problems involving contact between two boundaries, one contact surface is usually defined as the master surface, and the other is the slave surface. These two surfaces together comprise the surface-to-surface contact pair.

​

The slave surface is constrained against penetrating the master surface, and, usually, the slave surface is defined as the softer of the two. In the proposed model, as the anchor is much stiffer than the surrounding soil, all master elements were placed on the outer wall of the anchor, and all slave elements were on the surrounding soil contact surface.

​

Furthermore, the impermeability of the anchor must also be considered. As the pore pressure degree of freedom is only active in the soil elements, the FE program automatically finds the contact surface impermeable; hence, no fluid flow occurs between the soil and the anchor.

​

The interaction between the soil and the anchor is simulated using a penalty-type contact approach, whose main parameter is the friction coefficient between the surfaces <span id='cite-_Ref503887814'></span>[[#_Ref503887814|[22]]]. Karlsrud <span id="cite-_Ref535247756"></span>[[#_Ref535247756|[23]]] states that, for piles, the axial holding capacity can be obtained by calculating the shear resistance over the pile shaft concerning time and effective radial stress. The author assumes that the radial effective stress over the pile can be multiplied by a factor, <math display="inline">{f}_{cons}</math>, typically between 0.2 and 0.4, thus:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>{f}_{cons}\cdot {{\sigma }^{'}}_{r}(U)={s}_{u}(U)</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='Eq_4_1'></span>(5)

|}

​

​

where <math display="inline">{{\sigma }^{'}}_{r}(U)</math> is the effective radial stress as a function of the soil consolidation ratio, and <math display="inline">{s}_{u}(U)</math> is the soil undrained shear strength as a function of the soil consolidation ratio.

​

During the analyses performed in this study, it was observed that when a factor <math display="inline">{f}_{cons}</math> between 0.2 and 0.4 was applied, the results obtained for the fully consolidated soil were much higher than the values considering the soil with intact undrained shear strength. Hence, the friction factor was calibrated against the undrained shear strength corresponding to the soil's 100% consolidation rate. It was assumed that the undrained shear strength of the fully consolidated soil is equal to the one obtained by the ''α''-method from API <span id='cite-_Ref502324791'></span>[[#_Ref502324791|[21]]], which will be described in section 3. Thus, the calibrated friction factor, <math display="inline">{f}_{calib}</math>, can be written as:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>{f}_{calib}=\frac{\alpha \cdot {s}_{u}}{{{\sigma }^{'}}_{r}\left( U=100\%\right) }</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(6)

|}

​

​

where <math display="inline">{s}_{u}</math> is the undrained shear strength of the soil, <math display="inline">{{\sigma }^{'}}_{r}\left( U=\right. </math><math>\left. 100\%\right)</math>  is the effective radial stress corresponding to the fully consolidated soil and <math display="inline">\alpha</math>  is a non-dimensional factor.

​

==2.4. Initial Stress State of the Soil==

​

An important aspect of the torpedo anchor’s analysis is the simulation of the initial stress state of the soil, ''i.e.'', the stresses in the soil before applying any structural load to the anchor. As the proposed FE model does not simulate the anchor penetration in the soil, stress changes in the soil surrounding the anchor were assumed to be close to those produced from the expansion of a cylindrical cavity.

​

The Cavity Expansion Method (CEM) assumes that the strains induced during the anchor installation arise from an ideal expansion of a cylindrical cavity. Randolph and Wroth <span id='cite-_Ref534824256'></span>[[#_Ref534824256|[24]]] present a solution based on the assumption of a cylindrical cavity in a perfect elastic-perfectly plastic (EP) soil model. The authors assumed axial symmetry and plane strain conditions, which implies that only radial displacement of soil particles occurs.

​

Hill <span id='cite-_Ref534824674'></span>[[#_Ref534824674|[25]]] and Gibson and Anderson <span id='cite-_Ref534824683'></span>[[#_Ref534824683|[26]]] demonstrate the expressions for the stresses around an expanded cavity. For a cavity expanded from zero radii to a radius <math display="inline">{r}_{0}</math>, the radial and circumferential stress changes within the plastic zone are respectively given by:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\Delta {\sigma }_{r}={s}_{u}\left[ 1+ln\left( \frac{{G}_{50}}{{s}_{u}}\right) \right]</math> 

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref52371212'></span><span id='_Ref52371242'></span><span id='Eq_6'></span>(7)

|}

​

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\Delta {\sigma }_{\theta }={s}_{u}\left[ -1+ln\left( \frac{{G}_{50}}{{s}_{u}}\right) \right]</math> 

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(8)

|}

​

​

where <math display="inline">\Delta {\sigma }_{r}</math> is the increment of the effective radial stress, <math display="inline">\Delta {\sigma }_{\theta }</math> is the increment of the effective circumferential stress, <math display="inline">{G}_{50}</math> is the shear modulus at 50% mobilization of soil strength.

​

The relation between <math display="inline">{G}_{50}</math> and <math display="inline">{s}_{u}</math> can be estimated by the following empirical expression <span id='cite-_Ref2592654'></span>[[#_Ref2592654|[27]]]:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\frac{{G}_{50}}{{s}_{u}}\approx \frac{{e}^{\left( \frac{137-IP}{23}\right) }}{{\left[ 1+ln\left( 1+\frac{{\left( OCR-1\right) }^{3.2}}{26}\right) \right] }^{0.8}}</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(9)

|}

​

​

where <math display="inline">IP</math> is the plasticity index of the soil and <math display="inline">OCR</math> is the overconsolidation ratio of the soil.

​

Furthermore, Randolph and Wroth <span id='cite-_Ref534824256'></span>[[#_Ref534824256|[24]]] estimate the excess pore pressure by assuming that the mean effective stress remains constant under undrained conditions. The distribution of the initial excess pore pressure can then be written as:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\Delta {u}_{0}=\left\{ \begin{matrix}2{s}_{u}ln\left( \frac{{r}_{p}}{r}\right) ,\quad \quad \quad \, \, {r}_{0}\leq r\leq {r}_{p}\\0,\quad \quad \quad \quad \quad \quad \quad \quad \quad r>{r}_{p}\end{matrix}\right.</math> 

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref52371231'></span><span id='Eq_9'></span>(10)

|}

​

​

where <math display="inline">\Delta {u}_{0}</math> is the distribution of the initial excess pore pressure, <math display="inline">{r}_{0}</math> is the outer radius of the anchor, and <math display="inline">{r}_{p}</math> is the plasticized radius of the soil, given by:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>{r}_{p}={r}_{0}\sqrt{\frac{{G}_{50}}{{s}_{u}}}</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(11)

|}

​

​

The installation of a closed-ended pile is, in theory, modeled by the expansion of a cylindrical cavity with an initial radius of zero. In contrast, numerical calculations must begin with a finite cavity radius to avoid infinite circumferential strains <span id='cite-_Ref502687747'></span>[[#_Ref502687747|[9]]]. Carter ''et al''. <span id='cite-_Ref2698172'></span>[[#_Ref2698172|[28]]] found that doubling the cavity radius is adequate for EP and modified Cam Clay models. Hence, expanding a cavity from <math display="inline">{a}_{0}</math> to <math display="inline">{2a}_{0}</math> can approximate the cavity expansion from <math display="inline">r=</math><math>0</math> to <math display="inline">{r}_{0}</math>. Thus, in all analyses conducted in this study, a relationship of <math display="inline">{a}_{0}=</math><math>0.5{r}_{0}</math> was considered.

​

==2.5. Solution Procedure==

​

A typical FE mesh to predict the load capacity and the structural behavior of a finless torpedo anchor considering setup effects must account for contact, geometric, and material nonlinearities. The solution procedure adopted in this study follows the steps presented by Lopes <span id='cite-_Ref32331392'></span>[[#_Ref32331392|[29]]] and is described here in detail.

​

Complete setup analysis is divided into four different step.:

* <u>Step 1:</u> A geostatic analysis was invoked to obtain the ''in situ ''conditions of the soil mass. In this step, the effective self-weight of the soil was applied using the “body-force” option in Abaqus <span id="cite-_Ref534978548"></span>[[#_Ref534978548|[17]]]. This option was defined to prescribe the loading per unit volume over a body.

* <u>Step 2:</u> The expansion of the cylindrical cavity was performed. Part of the soil corresponding to the initial cylindrical cavity radius was removed and replaced with a rigid bar. Then, the bar moved against the soil to simulate the expansion of a cavity, as shown in <span id="cite-_Ref52370778"></span>[[#_Ref52370778|Fig.  6]].

* <u>Step 3:</u> The second step results were imported to another FE model using the MAP SOLUTION command in Abaqus <span id="cite-_Ref534978548"></span>[[#_Ref534978548|[17]]], where a new undeformed finite element mesh was adopted. This procedure was essential to continue the analysis as the CEM generates high strain rates in the model. Finally, the anchor was positioned into the soil mass and the setup analysis was conducted.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image6.png|360px]] </div>

​

<div id="_Ref52370778" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  6 – Detail of the rigid bar and its imposed displacement.''</div>

* <u>Step 4:</u> In this step, the holding capacity analysis was subdivided into three load steps. In the first load step, a time step of 0.001 seconds was employed to simulate the conditions in the soil immediately after anchor driving. The second load step involved continuing consolidation analysis for lengths of time ranging from 1 day to 18 000 days (50 years approximately). In this loading step, the excess pore pressure dissipated as time progressed, and the soil consolidated. As the excess pore pressure dissipated, the effective radial stress in the soil adjoining the anchor increased with time. This increase in the effective radial stress resulted in a higher frictional resistance between the soil and the anchor. Finally, the third load step involved assessing the pullout strength by applying a concentrated force into the reference node positioned above the top of the anchor and monitoring the vertical motion.

The proposed models had, on average, 45 200 nodes, 14 850 elements, and 105 000 degrees of freedom. Hence, the total time needed to perform a complete setup analysis was about 39 hours, considering that the models were executed sequentially. All analyses were carried out in an Intel® Core™ i7-6500 CPU [mailto:@ @] 2.50GHz 2.60GHz machine with 16 GB of RAM.

​

<span id='_Ref52370739'></span>

​

==3. Analytical estimation of the vertical holding capacity of torpedo anchors==

​

The ultimate holding capacity of torpedo anchors can be calculated as the sum of three terms. The first term refers to the total lateral resistance provided by the friction between the soil and the lateral area of the anchor. The second term is associated with the top resistance, and the last term is the submerged anchor’s weight.

​

According to Sousa ''et al.'' <span id='cite-_Ref31804075'></span>[[#_Ref31804075|[4]]], to estimate the response of a torpedo anchor loaded vertically, the use of the expressions proposed by API <span id='cite-_Ref502324791'></span>[[#_Ref502324791|[21]]] proved to be adequate. The ultimate holding capacity, <math display="inline">{Q}_{d}</math>, can be obtained by the following equation:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>{Q}_{d}=\int_{{L}_{t}}^{\, }\alpha {\cdot s}_{u}\left( z\right) \cdot {P}_{s}(z)\cdot dz+</math><math>9{\cdot s}_{u}{\cdot A}_{p}+{W}_{a}</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(12)

|}

​

​

where <math display="inline">{P}_{s}(z)</math> is the perimeter of the anchor at a corresponding position z, <math display="inline">{L}_{t}</math> is the total length of the anchor, <math display="inline">{A}_{p}</math> is the top surface area of the anchor, <math display="inline">{W}_{a}</math> is the submerged weight of the anchor, and <math display="inline">\alpha</math>  is a non-dimensional factor, calculated as:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\alpha =\left\{ \begin{matrix}0.5{\left( \frac{{s}_{u}}{{p'}_{0}}\right) }^{-0.5},\, \, \left( \frac{{s}_{u}}{{p'}_{0}}\right) \leq 1.0\\0.5{\left( \frac{{s}_{u}}{{p'}_{0}}\right) }^{-0.25},\, \, \left( \frac{{s}_{u}}{{p'}_{0}}\right) >1.0\end{matrix}\right.</math> 

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(13)

|}

​

​

with the constraint that <math display="inline">\alpha \leq 1</math>, where <math display="inline">{p'}_{0}</math> is the effective overburden pressure.

​

However, the presented expressions do not incorporate the setup effects. By combining the calculus of the ultimate holding capacity using the ''α''-method with the CEM equations, the effective radial stress of consolidation <math display="inline">{p'}_{0}</math> can be written as:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>{{p}^{'}}_{0}\left( t\right) ={s}_{u}\left( 1+2\cdot ln\frac{{r}_{p}}{{r}_{0}}\right) -</math><math>U\cdot 2\cdot {s}_{u}\cdot ln\frac{{r}_{p}}{{r}_{0}}+{{\sigma }^{'}}_{v0}\cdot {K}_{0}</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(14)

|}

​

​

where <math display="inline">{{\sigma }^{'}}_{v0}</math> is the effective vertical stress ''in-situ'', <math display="inline">{K}_{0}</math> is the coefficient of earth pressure at rest and <math display="inline">U</math> is the consolidation ratio of the soil.

​

Moreover, some analytical expressions were proposed to estimate setup effects in cylindrical piles <span id='cite-_Ref502153693'></span>[[#_Ref502153693|[10]]]. However, these empirical expressions are limited to specific applications and vary according to the degree of complexity and coverage.

​

Skov and Denver <span id='cite-_Ref534917435'></span>[[#_Ref534917435|[18]]] presented an empirical relationship in which the setup process varies linearly with the ''log'' of time, ''i. e.'':

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\frac{{R}_{u}(t)}{{R}_{0}}=1+A\cdot log\left( \frac{t}{{t}_{0}}\right)</math> 

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref52446153'></span><span id='_Ref52376815'></span><span id='Eq_10'></span><span id='Eq_15'></span>(15)

|}

​

​

where <math display="inline">{R}_{u}(t)</math> is the axial capacity at time <math display="inline">t</math> after driving, <math display="inline">{R}_{0}</math> is the axial capacity at the time <math display="inline">{t}_{0}</math>, <math display="inline">A</math> is a constant that depends on the soil type, and <math display="inline">{t}_{0}</math> is an empirical value measured in days.

​

Furthermore, Svinkin and Skov <span id='cite-_Ref534917465'></span>[[#_Ref534917465|[19]]] derived a new formulation<span id='cite-Eq_10'></span> to quantify the setup effects on the holding capacity of piles installed in cohesive soils. The formula proposed assumes that <math display="inline">{t}_{0}</math> equals to 0.1, and is given by:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\frac{{R}_{u}(t)}{{R}_{EOD}}=1+B\left[ \mathrm{{log}_{10}}\,(t)+1\right]</math> 

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref57295891'></span><span id='_Ref52376819'></span><span id='Eq_14'></span>(16)

|}

​

​

where <math display="inline">{R}_{EOD}</math> is the axial capacity at the end of driving and <math display="inline">B</math> is a constant similar to <math display="inline">A</math> (<span id='cite-Eq_15'></span>[[#Eq_15|15]]).

​

Recently, ABS <span id='cite-_Ref51515640'></span>[[#_Ref51515640|[20]]] presented an empirical formulation that provides a preliminary assessment of the regain of dynamically installed piles capacity after driving. ABS <span id='cite-_Ref51515640'></span>[[#_Ref51515640|[20]]] states that the regain of the pile capacity <math display="inline">{F}_{v}</math> can be conservatively estimated with:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>\frac{{F}_{v}-{W}_{s}}{{F}_{max}}=1.1-\frac{1.08}{1+{\left( \frac{T}{6.5}\right) }^{0.42}}</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(17)

|}

​

​

where <math display="inline">{W}_{s}</math> is the submerged weight of the pile, <math display="inline">{F}_{max}</math> is the long term maximum pile capacity, and <math display="inline">T</math> is a non-dimensional time defined as:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>T=\frac{{c}_{h}t}{{d}^{2}}</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(18)

|}

​

​

where <math display="inline">{c}_{h}</math> is the horizontal coefficient of consolidation, <math display="inline">t</math> is the consolidation time after installation, and <math display="inline">d</math> is the pile shaft diameter.

​

==4. Parametric Study==

​

==4.1. Description==

​

In this work, several FE analyses were performed to study the effect of three parameters on the holding capacity of the anchor while considering setup effects. These parameters are the:

​

:1. Soil consolidation time;

​

:2. Undrained shear strength of the soil;

​

:3. Coefficient of permeability of the soil.

​

The geometry of the anchor is presented in <span id='cite-_Ref52370792'></span>[[#_Ref52370792|Fig.  7]]. The thickness of the anchor shaft and the conical tip are 65 mm and 50 mm, respectively. The material behavior of the anchor was assumed linear and elastic with Young’s modulus equal to 210 GPa and Poisson ratio of 0.3. The anchor was discretized using 8-noded plane axisymmetric elements with translational degrees of freedom. The submerged self-weight of the anchor was kept constant and equal to 500 kN.

​

It is worth mentioning that, in this work, the anchor was supposed to be installed perfectly vertically, while, in practice, inclinations concerning the horizontal plane are often observed. This assumption and the consideration of a finless torpedo anchor allow the use of the proposed axisymmetric model.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

'' [[Image:Draft_Kronemberger Lopes_857637204-image7.png|240px]] ''</div>

​

<div id="_Ref52370792" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  7 – Main dimensions of the analyzed finless torpedo anchor.''</div>

​

The soil was considered purely cohesive and isotropic with a Poisson ratio of 0.4 and Young’s modulus varying with depth according to the following expression:

​

{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 

|-

| 

{| style="margin:auto;width: 100%;"

|-

| <math>E\left( z\right) =550\cdot {s}_{u}(z)</math>

|}

|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19)

|}

​

​

The soil is normally consolidated ( <math display="inline">OCR=1</math>), with a plasticity index of 30%. In all FE analyses, the rupture of the soil was verified with the Mohr-Coulomb criterion, <span id='cite-Eq_4'></span>[[#Eq_4|4]], by assuming an internal friction angle of 0° and that the cohesion equals the undrained shear strength of the soil. The initial void ratio of the soil was set equal to 1.5, and its submerged self-weight was 5.5 kN/m³. Additionally, three different shear strength profiles with a fixed permeability coefficient were considered in this study, as presented in <span id='cite-_Ref32422685'></span>[[#_Ref32422685|Table 1]].

​

<div id="_Ref32422685" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Table 1 – Main parameters considered in the analyses.''</div>

​

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 

|-

|  style="border-top: 1pt solid black;text-align: center;"|Analysis Index

|  style="border-top: 1pt solid black;text-align: center;"|Soil Label

|  style="border-top: 1pt solid black;text-align: center;"|Shear strength profile (kPa)

|  style="border-top: 1pt solid black;text-align: center;"|Permeability coefficient (m/day)

|-

|  style="text-align: center;"|1

|  style="text-align: center;"|A

|  style="text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+1.24\cdot z\right)</math> 

|  style="text-align: center;"|<math>2.5\times {10}^{-5}</math>

|-

|  style="text-align: center;"|2

|  style="text-align: center;"|B

|  style="text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+1.86\cdot z\right)</math> 

|  style="text-align: center;"|<math>2.5\times {10}^{-5}</math>

|-

|  style="text-align: center;"|3

|  style="text-align: center;"|C

|  style="text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+2.48\cdot z\right)</math> 

|  style="text-align: center;"|<math>2.5\times {10}^{-5}</math>

|-

|  style="text-align: center;"|4

|  style="text-align: center;"|A

|  style="text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+1.24\cdot z\right)</math> 

|  style="text-align: center;"|<math>2.5\times {10}^{-4}</math>

|-

|  style="border-bottom: 1pt solid black;text-align: center;"|5

|  style="border-bottom: 1pt solid black;text-align: center;"|A

|  style="border-bottom: 1pt solid black;text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+1.24\cdot z\right)</math> 

|  style="border-bottom: 1pt solid black;text-align: center;"|<math>2.5\times {10}^{-6}</math>

|}

​

​

An important aspect that influences the holding capacity of torpedo anchors is its embedment depth. This depth is a function of the geometry and weight of the anchor, the installation characteristics, including the anchor drop height, and the various hydrodynamic and geotechnical parameters, such as the undrained shear strength [29, 30]. Therefore, torpedo anchors with different physical or geometric characteristics or embedded in soils with different properties would have different embedment depths if all other installation parameters were kept constant.

​

According to ABS <span id='cite-_Ref51515640'></span>[[#_Ref51515640|[20]]], the anchor embedment depth can be assessed using True`s method <span id='cite-_Ref32496303'></span>[[#_Ref32496303|[31]]], accounting for strain rate and drag effects. Therefore, a driving analysis was performed, based on the model presented by True <span id='cite-_Ref32496303'></span>[[#_Ref32496303|[31]]], to verify the embedment depth of the finless torpedo anchors. Normal and tangential drag coefficients were assumed equal to 1.51 and 0, respectively, while the soil sensitivity was set constant and equal to 3.1. Moreover, an empirical shear strain rate ( <math display="inline">{S}_{e}</math>) of 3.5, coefficients of shear strain <math display="inline">{C}_{e}</math> equal to 1.4kPa.s and <math display="inline">{C}_{0}</math> equal to 0.07, and an adhesion coefficient equal to 0.07 were assumed in all driving analyses. A bearing capacity coefficient of 12 for the anchor tip was also considered.

​

In this study, aiming to only access the effects of the undrained shear strength of the soil on the holding capacity of torpedo anchors, a unique embedment depth of 28.5 m was assumed in all cases. This particular value of embedment depth can be achieved by varying the drop height of the torpedo anchor in the driving analysis, as proved in the embedment analyses. For instance, the drop height required for soils A, B, and C were found to be, respectively, 8 m, 48 m, and 190 m.

​

==4.2. Results==

​

==4.2.1. Undrained shear strength effects==

​

<span id='cite-_Ref52282642'></span>[[#_Ref52282642|Fig. 8]] to <span id='cite-_Ref52282962'></span>[[#_Ref52282962|Fig. 19]] present plots of the effective radial stresses and excess pore pressures of soils A, B, and C (analyses 1, 2, and 3 from Table 1, respectively), for different instants of time after the installation of the anchor. These results refer to the final step before pullout analyses. The anchor was removed from these figures to evaluate the soil's response properly.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image8.png|600px]] </div>

​

<div id="_Ref52282642" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 8 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, immediately after installation (0 days) – Soil A, analysis 1.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image9.png|600px]] </div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 9 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 10 days after installation – Soil A, analysis 1.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image10.png|600px]] </div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 10 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 150 days after installation – Soil A, analysis 1.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image11.png|534px]] </div>

​

<div id="_Ref52282955" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 11 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 18 000 days after installation – Soil A, analysis 1.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image12.png|600px]] </div>

​

<div id="_Ref52282651" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 12 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, immediately after installation (0 days) – Soil B, analysis 2.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image13.png|600px]] </div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 13 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 10 days after installation – Soil B, analysis 2.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image14.png|600px]] </div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 14 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 150 days after installation – Soil B, analysis 2.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image15.png|600px]] </div>

​

<div id="_Ref52282960" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 15 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 18 000 days after installation – Soil B, analysis 2.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image16.png|600px]] </div>

​

<div id="_Ref52282662" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 16 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, immediately after installation (0 days) – Soil C, analysis 3.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image17.png|600px]] </div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 17 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 10 days after installation – Soil C, analysis 3.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image18.png|600px]] </div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 18 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 150 days after installation – Soil C, analysis 3.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image19.png|534px]] </div>

​

<div id="_Ref52282962" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 19 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 18 000 days after installation – Soil C, analysis 3.''</div>

​

<span id='cite-_Ref52282642'></span>[[#_Ref52282642|Fig. 8]] to <span id='cite-_Ref52282962'></span>[[#_Ref52282962|Fig. 19]] indicate that, as time passed, the initial excess pore pressure generated during the installation of the anchor dissipated, and the effective radial stress of the soil around the anchor increased, which characterizes the setup phenomenon.

​

Immediately after installation, the excess pore pressure generated was higher because the undrained soil strength increased, as can be observed in <span id='cite-_Ref52282642'></span>[[#_Ref52282642|Fig. 8]], <span id='cite-_Ref52282651'></span>[[#_Ref52282651|Fig. 12]], and <span id='cite-_Ref52282662'></span>[[#_Ref52282662|Fig. 16]]. Moreover, the magnitude of the radial effective stresses of the soil was higher for soils with a greater undrained shear strength profile (<span id='cite-_Ref52282955'></span>[[#_Ref52282955|Fig. 11]], <span id='cite-_Ref52282960'></span>[[#_Ref52282960|Fig. 15]], and <span id='cite-_Ref52282962'></span>[[#_Ref52282962|Fig. 19]]). This response was expected because <span id='cite-Eq_6'></span><span id="cite-Eq_9"></span>the CEM analytical formulations are directly related to the undrained shear strength of the soil, as indicated in Eqs.[[#Eq_6|7]] and [[#Eq_9|10]].

​

The plots presented in <span id='cite-_Ref52282642'></span>[[#_Ref52282642|Fig. 8]] to <span id='cite-_Ref52282962'></span>[[#_Ref52282962|Fig. 19]] indicate that the initial excess pore pressure was almost dissipated approximately 18 000 days after the anchor installation for all soils considered. Hence, even if the soil presents a higher magnitude for the initial excess pore pressure, the total time needed to achieve a fully consolidated scenario remains the same. It occurs mainly due to the permeability coefficient of the soil, confirming that this parameter plays an important role in consolidation analysis. Moreover, considering the same coefficient of permeability for all soils, it can be stated that the regain of the soil effective radial stresses was faster as the undrained shear strength increased.

​

Setup effects can be better observed in <span id='cite-_Ref52284960'></span>[[#_Ref52284960|Fig.  20]]. This figure shows the variation of the normalized displacement at the top of the anchor with the normalized applied load. The normalized displacement is defined as the ratio between the total displacement at the top of the anchor <math display="inline">d</math> and its shaft external diameter <math display="inline">{D}_{a}</math>. The normalized applied load is defined as the ratio between the total applied load at the reference node <math display="inline">{Q}_{t}</math> and the submerged weight of the anchor <math display="inline">{W}_{a}</math>. These curves refer to analysis 1 with the finless torpedo anchor embedded in soil A. In this figure, different curves indicate different consolidation times after the installation of the anchor. Additionally, the value obtained with the &#x03b1;-method is presented, which corresponds to a limit value, ''i. e.'', the pullout capacity after complete consolidation. Note that, as time passes, the solution converges to the ultimate value given by the analytical solution.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image20.png|456px]] </div>

​

<div id="_Ref52284960" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  20 – Holding capacity regain of the torpedo anchor in Soil A.''</div>

​

Setup curves referent to analyses 1 to 3, soils A, B, and C respectively, are presented in <span id='cite-_Ref52285415'></span>[[#_Ref52285415|Fig.  21]]. This figure indicates that the torpedo anchor’s holding capacity increased when considering soils with higher undrained shear strengths. Moreover, by evaluating the setup curves of soils A, B, and C, a higher ratio of recovery of the holding capacity when considering soils with higher undrained soil strengths can be noticed.

​

<span id='cite-_Ref52285415'></span>[[#_Ref52285415|Fig.  21]] also indicates that all setup curves obtained throughout this study converged to the respective analytical response of the ultimate holding capacity (''i. e.'', soil fully consolidated) obtained using the calibrated method based on API <span id='cite-_Ref502324791'></span>[[#_Ref502324791|[21]]].

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image21.png|600px]] </div>

​

<span id='_Ref52285415'></span><div id="_Ref52287206" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  21 – Setup curves obtained for Soils A, B, and C.''</div>

​

Phases 2 and 3 of the setup process (<span id='cite-_Ref52286155'></span>[[#_Ref52286155|Fig.  3]]) can be observed in <span id='cite-_Ref52285415'></span>[[#_Ref52285415|Fig.  21]]. However, the first phase of the setup could not be identified in the curves. It is important to highlight that the thixotropic effects were not considered in the analyses performed in this study. According to Bergset <span id='cite-_Ref502687747'></span>[[#_Ref502687747|[9]]], Komurka and Wagner <span id='cite-_Ref502153693'></span>[[#_Ref502153693|[10]]], and Richardson ''et al.'' <span id='cite-_Ref31897166'></span>[[#_Ref31897166|[12]]], thixotropic effects are more relevant in the first instants of time after installation, which corresponds to the first phase of setup. Thus, to correctly identify the first phase of setup, thixotropy effects need to be considered, which is a point that deserves future studies.

​

The transition between the second and the third phases of setup can be observed when the setup curve changes its inclination. For all the setup curves presented in <span id='cite-_Ref52285415'></span>[[#_Ref52285415|Fig.  21]], this transition occurs about 5 000 days after installation.

​

<span id='cite-_Ref34065797'></span>[[#_Ref34065797|Table 2]] summarizes the main results of <span id='cite-_Ref52285415'></span>[[#_Ref52285415|Fig.  21]]. This table indicates an increasing rate of the holding capacity during the first 15 days after the anchor’s installation in all cases analyzed. However, this rate behaves more linearly after 10 or 15 days after the anchor installation. For the different soils considered, the increase in the holding capacity regarding its initial value for the range of time between 15 and 100 days was approximately 11%.

​

<div id="_Ref34065797" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Table 2 – Summary of the normalized holding capacity obtained for soils A, B, and C.''</div>

​

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 

|-

|  style="border-top: 1pt solid black;text-align: center;"|Analysis

|  style="border-top: 1pt solid black;text-align: center;"|0 day

|  style="border-top: 1pt solid black;text-align: center;"|1 day

|  style="border-top: 1pt solid black;text-align: center;"|15 days

|  style="border-top: 1pt solid black;text-align: center;"|30 days

|  style="border-top: 1pt solid black;text-align: center;"|60 days

|  style="border-top: 1pt solid black;text-align: center;"|100 days

|-

|  style="text-align: center;"|1

|  style="text-align: center;"|1.00

|  style="text-align: center;"|1.12

|  style="text-align: center;"|1.30

|  style="text-align: center;"|1.35

|  style="text-align: center;"|1.39

|  style="text-align: center;"|1.42

|-

|  style="text-align: center;"|2

|  style="text-align: center;"|1.00

|  style="text-align: center;"|1.17

|  style="text-align: center;"|1.40

|  style="text-align: center;"|1.43

|  style="text-align: center;"|1.46

|  style="text-align: center;"|1.48

|-

|  style="border-bottom: 1pt solid black;text-align: center;"|3

|  style="border-bottom: 1pt solid black;text-align: center;"|1.00

|  style="border-bottom: 1pt solid black;text-align: center;"|1.20

|  style="border-bottom: 1pt solid black;text-align: center;"|1.45

|  style="border-bottom: 1pt solid black;text-align: center;"|1.51

|  style="border-bottom: 1pt solid black;text-align: center;"|1.56

|  style="border-bottom: 1pt solid black;text-align: center;"|1.59

|}

​

​

By comparing <span id='cite-_Ref34065797'></span>[[#_Ref34065797|Table 2]] and <span id='cite-_Ref52285415'></span>[[#_Ref52285415|Fig.  21]], it is possible to affirm that the regain rates of the holding capacity are higher in the initial moments. Immediately after, these rates start to decrease. It is important to emphasize that this conclusion is only valid because the same permeability coefficient was considered for all soils analyzed.

​

==4.2.2. Permeability coefficient effects==

​

In the investigation of the impact of the soil permeability coefficient in the ultimate holding capacity of the torpedo anchor, only the undrained shear strength profile of Soil A (<span id='cite-_Ref32422685'></span>[[#_Ref32422685|Table 1]]) was considered. All other properties of the model were the same described in section 4.1.

​

<span id='cite-_Ref52371035'></span>[[#_Ref52371035|Fig. 22]] shows the soil consolidation ratio considering a permeability coefficient of <math display="inline">2.5\times {10}^{-5}</math> m/day. The figure shows three different consolidation curves related to a specific point located along the anchor length. These points are also illustrated in this figure. It is emphasized that the points were chosen to ensure that the anchor geometry's boundaries do not influence the results.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image22.png|534px]] </div>

​

<div id="_Ref52371035" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 22 – Consolidation ratio over time in Soil A.''</div>

​

As can be observed in <span id='cite-_Ref52371035'></span>[[#_Ref52371035|Fig. 22]], the rate of excess pore pressure dissipation is slightly different at each of the points analyzed. However, all curves follow the same trend. Moreover, complete soil consolidation occurred practically at the same time. In this case, a time of 10 000 days is representative of this phenomenon.

​

[[#_img-23|Fig. 23]]<span id='cite-_Ref52371035'></span> compares the consolidation ratio along time for soil A, considering different permeability coefficients. These results refer to the point located at 7.5 m distant from the top of the anchor, as illustrated in [[#_Ref52371035|Fig. 22]]. It can be observed that the time needed to achieve a certain consolidation ratio of the soil is lower when the permeability coefficient increases. On the other hand, full consolidation time is higher when the permeability coefficient decreases.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image23.png|486px]] </div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  23 – Consolidation ratios for different permeability coefficients at 7.5m from the top of the anchor.''</div>

​

The total times needed for consolidation ratios of 0.1 to 0.99 are summarized in <span id='cite-_Ref34133953'></span>[[#_Ref34133953|Table 3]] for the different permeability coefficients considered. This table indicates that when the permeability coefficient increases an order of magnitude, the time needed to achieve a specific consolidation ratio decreases at the same ratio. Analogously, when decreasing one order of magnitude of the permeability coefficient, the total time required to achieve a certain consolidation ratio is increased by the same ratio.

​

<div id="_Ref34133953" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Table 3 – Time (in days) needed to achieve different consolidation ratios: 7.5 m above the top of the anchor.''</div>

​

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 

|-

|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<math display="inline">\mathit{\boldsymbol{k}}</math> (m/dia)

|  colspan='6'  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|Consolidation ratio, <math display="inline">\mathit{\boldsymbol{U}}</math>

|-

|  style="text-align: center;vertical-align: top;"|'''0.10'''

|  style="text-align: center;vertical-align: top;"|'''0.25'''

|  style="text-align: center;vertical-align: top;"|'''0.50'''

|  style="text-align: center;vertical-align: top;"|'''0.75'''

|  style="text-align: center;vertical-align: top;"|'''0.90'''

|  style="text-align: center;vertical-align: top;"|'''0.99'''

|-

|  style="text-align: center;vertical-align: top;"|2.50E-04

|  style="text-align: center;vertical-align: top;"|3.00E-02

|  style="text-align: center;vertical-align: top;"|3.30E-01

|  style="text-align: center;vertical-align: top;"|3.79E+00

|  style="text-align: center;vertical-align: top;"|6.00E+01

|  style="text-align: center;vertical-align: top;"|3.08E+02

|  style="text-align: center;vertical-align: top;"|1.41E+03

|-

|  style="text-align: center;vertical-align: top;"|2.50E-05

|  style="text-align: center;vertical-align: top;"|3.10E-01

|  style="text-align: center;vertical-align: top;"|3.26E+00

|  style="text-align: center;vertical-align: top;"|3.97E+01

|  style="text-align: center;vertical-align: top;"|6.23E+02

|  style="text-align: center;vertical-align: top;"|3.09E+03

|  style="text-align: center;vertical-align: top;"|1.52E+04

|-

|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.50E-06

|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3.09E+00

|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3.18E+01

|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4.04E+02

|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|6.32E+03

|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3.10E+04

|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1.12E+05

|}

​

Another important aspect that should be observed in <span id='cite-_Ref34133953'></span>[[#_Ref34133953|Table 3]] is the variation of consolidation ratio over time. It can be noted that, in the transition between the consolidation ratios chosen, a change of an order of magnitude into time scale occurs. The succeeding values of time are not precisely the previous value multiplied by a factor of 10, but they are relatively close to that.

​

==4.2.3. Comparison between FEM and analytical methodologies==

​

For this investigation, a unique permeability coefficient of <math display="inline">2.5\times {10}^{-5}</math> m/day was employed in all analyses. Moreover, different undrained shear strength profiles were considered, and the effects of these profiles on the holding capacity of the torpedo anchors were assessed.

​

<span id='cite-_Ref52371266'></span>[[#_Ref52371266|Fig.  24]] compares the radial stresses obtained with the FEM and the CEM analytical formulation, Eqs. <span id='cite-_Ref52371242'></span>[[#_Ref52371242|(7)]] and <span id='cite-_Ref52371231'></span>[[#_Ref52371231|(10)]]. In this figure, the radial stress was normalized by the undrained shear strength of the soil. Furthermore, the radial distance was normalized against the external radius of the anchor. In general, both methodologies are in good agreement throughout the radial distance within the soil. Although, a significant difference can be noticed in the region close to the anchor shaft in soil B.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image24.png|450px]] </div>

​

<div id="_Ref52371266" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  24 – Comparison between total radial stresses – Soils A and B.''</div>

​

Aiming at employing the analytical methodology based on the α-method from API <span id='cite-_Ref502324791'></span>[[#_Ref502324791|[21]]], CEM equations were used to obtain the effective radial stress over the external anchor wall along time and, thus, to calculate the values of the remaining parameters at each instant of time selected.

​

As mentioned before, the CEM is not suitable for estimating the stresses at the top and tip of the anchor. However, due to the lack of a better analytical formulation to estimate the stresses at these two regions, CEM equations were adopted in this study.

​

Input parameters of combined holding capacity <math display="inline">{R}_{0}</math> and <math display="inline">{R}_{EOD}</math>, used in the empirical formulations of Skov and Denver <span id='cite-_Ref534917435'></span>[[#_Ref534917435|[18]]], Eq. <span id='cite-_Ref52376815'></span>[[#_Ref52376815|(15)]], and Svinkin and Skov <span id='cite-_Ref534917465'></span>[[#_Ref534917465|[19]]], Eq. <span id='cite-_Ref52376819'></span>[[#_Ref52376819|(16)]], were obtained from the analyses using the FEM for time instants of 1 day and 0.01 day, respectively, after the installation of the anchor.

​

<span id='cite-_Ref52371311'></span>[[#_Ref52371311|Fig.  25]] to <span id='cite-_Ref52371314'></span>[[#_Ref52371314|Fig.  27]] present comparative plots, for different soil undrained shear strength profiles considering the FE and analytical predictions. In these figures, the holding capacity of the anchor was normalized by the submerged weight of the anchor (500 kN).

​

In general, all methodologies have similar behavior, except for the initial instants of time for the ABS <span id='cite-_Ref51515640'></span>[[#_Ref51515640|[20]]] formulation, which starts with a lower holding capacity value. Results obtained from FEM analyses estimated pullout capacities slightly lower for a consolidation time of 18 000 days when compared to the other methodologies. It is emphasized that for all soils analyzed, only the ABS <span id='cite-_Ref51515640'></span>[[#_Ref51515640|[20]]] formulation and the analytical methodology based on α-method of API <span id='cite-_Ref502324791'></span>[[#_Ref502324791|[21]]] did not present a predominant linear behavior.

​

Skov and Denver [17] empirical formulation, Eq. <span id='cite-_Ref52376815'></span>[[#_Ref52376815|(15)]], presented an excellent agreement for the first ten days after anchor installation. On the other hand, the analytical methodology based on the α-method of API <span id='cite-_Ref502324791'></span>[[#_Ref502324791|[21]]] showed better approaches for consolidation times over 100 days.

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image25.png|600px]] </div>

​

<div id="_Ref52371311" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  25 – Comparison between methodologies – Soil A.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image26.png|600px]] </div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig. 26 – Comparison between methodologies – Soil B.''</div>

​

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

 [[Image:Draft_Kronemberger Lopes_857637204-image27.png|600px]] </div>

​

<div id="_Ref52371314" class="center" style="width: auto; margin-left: auto; margin-right: auto;">

''Fig.  27 – Comparison between methodologies – Soil C.''</div>

​

The undrained shear strength of the soil has a major impact on the ABS <span id='cite-_Ref51515640'></span>[[#_Ref51515640|[20]]] formulation compared to the other methodologies. As shown in <span id='cite-_Ref52371311'></span>[[#_Ref52371311|Fig.  25]] to <span id='cite-_Ref52371314'></span>[[#_Ref52371314|Fig.  27]], the pullout capacity increases at a higher rate compared to other formulations. For soil C, the ABS <span id='cite-_Ref51515640'></span>[[#_Ref51515640|[20]]] formulation approaches the Skov and Denver <span id='cite-_Ref534917435'></span>[[#_Ref534917435|[18]]] formulation for a consolidation time of 18 000 days, while, for soils A and B, it approaches the other formulations.