You do not have permission to edit this page, for the following reason:

You are not allowed to execute the action you have requested.


You can view and copy the source of this page.

x
 
1
<!-- metadata commented in wiki content
2
3
4
==Numerical investigation of setup effects in the vertical holding capacity of finless torpedo anchors==
5
6
Guilherme K. Lopes<sup>a,*</sup>, José Renato M. de Sousa<sup>a</sup>, and Gilberto B. Ellwanger<sup>a</sup>
7
8
<sup>a </sup>Civil Engineering Department – PEC, Federal University of Rio de Janeiro, COPPE/UFRJ, 21945-970, Rio de Janeiro, RJ, Brazil.
9
10
*e-mail: [mailto:kronemberger.guilherme@laceo.coppe.ufrj.br kronemberger.guilherme@laceo.coppe.ufrj.br], Tel.: +55 24992272120
11
-->
12
13
==Abstract==
14
15
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.
16
17
'''Keywords''': Torpedo anchor, setup, finite elements, analytical models
18
19
==1. Introduction==
20
21
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.
22
23
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]]].
24
25
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
26
 [[Image:Draft_Kronemberger Lopes_857637204-image1.png|600px]] </div>
27
28
<div id="_Ref52370445" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
29
''Fig. 1 - Typical torpedo anchor: (a) conical tip and (b) padeye details ''<span id='cite-_Ref31804075'></span>[[#_Ref31804075|''[4]'']]''.''</div>
30
31
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.
32
33
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
34
 [[Image:Draft_Kronemberger Lopes_857637204-image2.png|384px]] </div>
35
36
<div id="_Ref52370458" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
37
''Fig.  2 – Typical installation procedure of a finless torpedo anchor.''</div>
38
39
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.
40
41
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.
42
43
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].
44
45
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
46
 [[Image:Draft_Kronemberger Lopes_857637204-image3.png|600px]] </div>
47
48
<div id="_Ref52286155" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
49
''Fig.  3  – Example of a setup curve.''</div>
50
51
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.
52
53
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.
54
55
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.
56
57
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.
58
59
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]].
60
61
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.
62
63
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
64
 [[Image:Draft_Kronemberger Lopes_857637204-image4.png|510px]] </div>
65
66
<div id="_Ref52370499" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
67
''Fig.  4 – General view of the FE model.''</div>
68
69
==2. Finite Element Approach==
70
71
==2.1. Soil Modeling==
72
73
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:
74
75
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
76
|-
77
| 
78
{| style="margin:auto;width: 100%;"
79
|-
80
| <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>
81
|}
82
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(1)
83
|}
84
85
86
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.
87
88
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.'':
89
90
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
91
|-
92
| 
93
{| style="margin:auto;width: 100%;"
94
|-
95
| <math>sn{v}_{w}=-k\cdot \frac{\partial \phi }{\partial x}</math>
96
|}
97
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(2)
98
|}
99
100
101
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:
102
103
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
104
|-
105
| 
106
{| style="margin:auto;width: 100%;"
107
|-
108
| <math>\phi \overset{\underset{\mathrm{def}}{}}{=}z+\frac{{u}_{w}}{g{\rho }_{w}}</math>
109
|}
110
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(3)
111
|}
112
113
114
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.
115
116
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.
117
118
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.'':
119
120
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
121
|-
122
| 
123
{| style="margin:auto;width: 100%;"
124
|-
125
| <math>\tau ={c}^{'}-{\sigma }^{'}.tan\phi '</math>
126
|}
127
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='Eq_4'></span>(4)
128
|}
129
130
131
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.
132
133
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>.
134
135
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.
136
137
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.
138
139
==2.2. Anchor Modeling==
140
141
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.
142
143
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]].
144
145
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
146
 [[Image:Draft_Kronemberger Lopes_857637204-image5.png|216px]] </div>
147
148
<div id="_Ref52370593" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
149
''Fig.  5 – Detail of the point of load application.''</div>
150
151
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.
152
153
==2.3. Anchor–Soil Interaction==
154
155
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.
156
157
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.
158
159
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.
160
161
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:
162
163
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
164
|-
165
| 
166
{| style="margin:auto;width: 100%;"
167
|-
168
| <math>{f}_{cons}\cdot {{\sigma }^{'}}_{r}(U)={s}_{u}(U)</math>
169
|}
170
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='Eq_4_1'></span>(5)
171
|}
172
173
174
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.
175
176
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:
177
178
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
179
|-
180
| 
181
{| style="margin:auto;width: 100%;"
182
|-
183
| <math>{f}_{calib}=\frac{\alpha \cdot {s}_{u}}{{{\sigma }^{'}}_{r}\left( U=100\%\right) }</math>
184
|}
185
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(6)
186
|}
187
188
189
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.
190
191
==2.4. Initial Stress State of the Soil==
192
193
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.
194
195
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.
196
197
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:
198
199
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
200
|-
201
| 
202
{| style="margin:auto;width: 100%;"
203
|-
204
| <math>\Delta {\sigma }_{r}={s}_{u}\left[ 1+ln\left( \frac{{G}_{50}}{{s}_{u}}\right) \right]</math> 
205
|}
206
|  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)
207
|}
208
209
210
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
211
|-
212
| 
213
{| style="margin:auto;width: 100%;"
214
|-
215
| <math>\Delta {\sigma }_{\theta }={s}_{u}\left[ -1+ln\left( \frac{{G}_{50}}{{s}_{u}}\right) \right]</math> 
216
|}
217
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(8)
218
|}
219
220
221
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.
222
223
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]]]:
224
225
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
226
|-
227
| 
228
{| style="margin:auto;width: 100%;"
229
|-
230
| <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>
231
|}
232
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(9)
233
|}
234
235
236
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.
237
238
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:
239
240
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
241
|-
242
| 
243
{| style="margin:auto;width: 100%;"
244
|-
245
| <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> 
246
|}
247
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|<span id='_Ref52371231'></span><span id='Eq_9'></span>(10)
248
|}
249
250
251
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:
252
253
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
254
|-
255
| 
256
{| style="margin:auto;width: 100%;"
257
|-
258
| <math>{r}_{p}={r}_{0}\sqrt{\frac{{G}_{50}}{{s}_{u}}}</math>
259
|}
260
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(11)
261
|}
262
263
264
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.
265
266
==2.5. Solution Procedure==
267
268
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.
269
270
Complete setup analysis is divided into four different step.:
271
* <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.
272
* <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]].
273
* <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.
274
275
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
276
 [[Image:Draft_Kronemberger Lopes_857637204-image6.png|360px]] </div>
277
278
<div id="_Ref52370778" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
279
''Fig.  6 – Detail of the rigid bar and its imposed displacement.''</div>
280
* <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.
281
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.
282
283
<span id='_Ref52370739'></span>
284
285
==3. Analytical estimation of the vertical holding capacity of torpedo anchors==
286
287
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.
288
289
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:
290
291
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
292
|-
293
| 
294
{| style="margin:auto;width: 100%;"
295
|-
296
| <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>
297
|}
298
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(12)
299
|}
300
301
302
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:
303
304
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
305
|-
306
| 
307
{| style="margin:auto;width: 100%;"
308
|-
309
| <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> 
310
|}
311
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(13)
312
|}
313
314
315
with the constraint that <math display="inline">\alpha \leq 1</math>, where <math display="inline">{p'}_{0}</math> is the effective overburden pressure.
316
317
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:
318
319
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
320
|-
321
| 
322
{| style="margin:auto;width: 100%;"
323
|-
324
| <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>
325
|}
326
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(14)
327
|}
328
329
330
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.
331
332
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.
333
334
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.'':
335
336
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
337
|-
338
| 
339
{| style="margin:auto;width: 100%;"
340
|-
341
| <math>\frac{{R}_{u}(t)}{{R}_{0}}=1+A\cdot log\left( \frac{t}{{t}_{0}}\right)</math> 
342
|}
343
|  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)
344
|}
345
346
347
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.
348
349
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:
350
351
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
352
|-
353
| 
354
{| style="margin:auto;width: 100%;"
355
|-
356
| <math>\frac{{R}_{u}(t)}{{R}_{EOD}}=1+B\left[ \mathrm{{log}_{10}}\,(t)+1\right]</math> 
357
|}
358
|  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)
359
|}
360
361
362
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]]).
363
364
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:
365
366
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
367
|-
368
| 
369
{| style="margin:auto;width: 100%;"
370
|-
371
| <math>\frac{{F}_{v}-{W}_{s}}{{F}_{max}}=1.1-\frac{1.08}{1+{\left( \frac{T}{6.5}\right) }^{0.42}}</math>
372
|}
373
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(17)
374
|}
375
376
377
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:
378
379
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
380
|-
381
| 
382
{| style="margin:auto;width: 100%;"
383
|-
384
| <math>T=\frac{{c}_{h}t}{{d}^{2}}</math>
385
|}
386
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(18)
387
|}
388
389
390
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.
391
392
==4. Parametric Study==
393
394
==4.1. Description==
395
396
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:
397
398
:1. Soil consolidation time;
399
400
:2. Undrained shear strength of the soil;
401
402
:3. Coefficient of permeability of the soil.
403
404
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.
405
406
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.
407
408
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
409
'' [[Image:Draft_Kronemberger Lopes_857637204-image7.png|240px]] ''</div>
410
411
<div id="_Ref52370792" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
412
''Fig.  7 – Main dimensions of the analyzed finless torpedo anchor.''</div>
413
414
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:
415
416
{| class="formulaSCP" style="width: 100%;margin: 1em auto 0.1em auto;width: 100%;text-align: center;" 
417
|-
418
| 
419
{| style="margin:auto;width: 100%;"
420
|-
421
| <math>E\left( z\right) =550\cdot {s}_{u}(z)</math>
422
|}
423
|  style="text-align: right;width: 5px;text-align: right;white-space: nowrap;"|(19)
424
|}
425
426
427
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]].
428
429
<div id="_Ref32422685" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
430
''Table 1 – Main parameters considered in the analyses.''</div>
431
432
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
433
|-
434
|  style="border-top: 1pt solid black;text-align: center;"|Analysis Index
435
|  style="border-top: 1pt solid black;text-align: center;"|Soil Label
436
|  style="border-top: 1pt solid black;text-align: center;"|Shear strength profile (kPa)
437
|  style="border-top: 1pt solid black;text-align: center;"|Permeability coefficient (m/day)
438
|-
439
|  style="text-align: center;"|1
440
|  style="text-align: center;"|A
441
|  style="text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+1.24\cdot z\right)</math> 
442
|  style="text-align: center;"|<math>2.5\times {10}^{-5}</math>
443
|-
444
|  style="text-align: center;"|2
445
|  style="text-align: center;"|B
446
|  style="text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+1.86\cdot z\right)</math> 
447
|  style="text-align: center;"|<math>2.5\times {10}^{-5}</math>
448
|-
449
|  style="text-align: center;"|3
450
|  style="text-align: center;"|C
451
|  style="text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+2.48\cdot z\right)</math> 
452
|  style="text-align: center;"|<math>2.5\times {10}^{-5}</math>
453
|-
454
|  style="text-align: center;"|4
455
|  style="text-align: center;"|A
456
|  style="text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+1.24\cdot z\right)</math> 
457
|  style="text-align: center;"|<math>2.5\times {10}^{-4}</math>
458
|-
459
|  style="border-bottom: 1pt solid black;text-align: center;"|5
460
|  style="border-bottom: 1pt solid black;text-align: center;"|A
461
|  style="border-bottom: 1pt solid black;text-align: center;"|<math>{s}_{u}\left( z\right) =\left( 3+1.24\cdot z\right)</math> 
462
|  style="border-bottom: 1pt solid black;text-align: center;"|<math>2.5\times {10}^{-6}</math>
463
|}
464
465
466
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.
467
468
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.
469
470
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.
471
472
==4.2. Results==
473
474
==4.2.1. Undrained shear strength effects==
475
476
<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.
477
478
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
479
 [[Image:Draft_Kronemberger Lopes_857637204-image8.png|600px]] </div>
480
481
<div id="_Ref52282642" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
482
''Fig. 8 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, immediately after installation (0 days) – Soil A, analysis 1.''</div>
483
484
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
485
 [[Image:Draft_Kronemberger Lopes_857637204-image9.png|600px]] </div>
486
487
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
488
''Fig. 9 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 10 days after installation – Soil A, analysis 1.''</div>
489
490
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
491
 [[Image:Draft_Kronemberger Lopes_857637204-image10.png|600px]] </div>
492
493
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
494
''Fig. 10 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 150 days after installation – Soil A, analysis 1.''</div>
495
496
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
497
 [[Image:Draft_Kronemberger Lopes_857637204-image11.png|534px]] </div>
498
499
<div id="_Ref52282955" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
500
''Fig. 11 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 18 000 days after installation – Soil A, analysis 1.''</div>
501
502
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
503
 [[Image:Draft_Kronemberger Lopes_857637204-image12.png|600px]] </div>
504
505
<div id="_Ref52282651" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
506
''Fig. 12 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, immediately after installation (0 days) – Soil B, analysis 2.''</div>
507
508
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
509
 [[Image:Draft_Kronemberger Lopes_857637204-image13.png|600px]] </div>
510
511
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
512
''Fig. 13 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 10 days after installation – Soil B, analysis 2.''</div>
513
514
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
515
 [[Image:Draft_Kronemberger Lopes_857637204-image14.png|600px]] </div>
516
517
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
518
''Fig. 14 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 150 days after installation – Soil B, analysis 2.''</div>
519
520
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
521
 [[Image:Draft_Kronemberger Lopes_857637204-image15.png|600px]] </div>
522
523
<div id="_Ref52282960" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
524
''Fig. 15 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 18 000 days after installation – Soil B, analysis 2.''</div>
525
526
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
527
 [[Image:Draft_Kronemberger Lopes_857637204-image16.png|600px]] </div>
528
529
<div id="_Ref52282662" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
530
''Fig. 16 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, immediately after installation (0 days) – Soil C, analysis 3.''</div>
531
532
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
533
 [[Image:Draft_Kronemberger Lopes_857637204-image17.png|600px]] </div>
534
535
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
536
''Fig. 17 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 10 days after installation – Soil C, analysis 3.''</div>
537
538
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
539
 [[Image:Draft_Kronemberger Lopes_857637204-image18.png|600px]] </div>
540
541
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
542
''Fig. 18 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 150 days after installation – Soil C, analysis 3.''</div>
543
544
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
545
 [[Image:Draft_Kronemberger Lopes_857637204-image19.png|534px]] </div>
546
547
<div id="_Ref52282962" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
548
''Fig. 19 – Effective radial stress (S11) and excess pore pressure (POR) in kPa, 18 000 days after installation – Soil C, analysis 3.''</div>
549
550
<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.
551
552
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]].
553
554
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.
555
556
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.
557
558
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
559
 [[Image:Draft_Kronemberger Lopes_857637204-image20.png|456px]] </div>
560
561
<div id="_Ref52284960" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
562
''Fig.  20 – Holding capacity regain of the torpedo anchor in Soil A.''</div>
563
564
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.
565
566
<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]]].
567
568
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
569
 [[Image:Draft_Kronemberger Lopes_857637204-image21.png|600px]] </div>
570
571
<span id='_Ref52285415'></span><div id="_Ref52287206" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
572
''Fig.  21 – Setup curves obtained for Soils A, B, and C.''</div>
573
574
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.
575
576
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.
577
578
<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%.
579
580
<div id="_Ref34065797" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
581
''Table 2 – Summary of the normalized holding capacity obtained for soils A, B, and C.''</div>
582
583
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
584
|-
585
|  style="border-top: 1pt solid black;text-align: center;"|Analysis
586
|  style="border-top: 1pt solid black;text-align: center;"|0 day
587
|  style="border-top: 1pt solid black;text-align: center;"|1 day
588
|  style="border-top: 1pt solid black;text-align: center;"|15 days
589
|  style="border-top: 1pt solid black;text-align: center;"|30 days
590
|  style="border-top: 1pt solid black;text-align: center;"|60 days
591
|  style="border-top: 1pt solid black;text-align: center;"|100 days
592
|-
593
|  style="text-align: center;"|1
594
|  style="text-align: center;"|1.00
595
|  style="text-align: center;"|1.12
596
|  style="text-align: center;"|1.30
597
|  style="text-align: center;"|1.35
598
|  style="text-align: center;"|1.39
599
|  style="text-align: center;"|1.42
600
|-
601
|  style="text-align: center;"|2
602
|  style="text-align: center;"|1.00
603
|  style="text-align: center;"|1.17
604
|  style="text-align: center;"|1.40
605
|  style="text-align: center;"|1.43
606
|  style="text-align: center;"|1.46
607
|  style="text-align: center;"|1.48
608
|-
609
|  style="border-bottom: 1pt solid black;text-align: center;"|3
610
|  style="border-bottom: 1pt solid black;text-align: center;"|1.00
611
|  style="border-bottom: 1pt solid black;text-align: center;"|1.20
612
|  style="border-bottom: 1pt solid black;text-align: center;"|1.45
613
|  style="border-bottom: 1pt solid black;text-align: center;"|1.51
614
|  style="border-bottom: 1pt solid black;text-align: center;"|1.56
615
|  style="border-bottom: 1pt solid black;text-align: center;"|1.59
616
|}
617
618
619
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.
620
621
==4.2.2. Permeability coefficient effects==
622
623
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.
624
625
<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.
626
627
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
628
 [[Image:Draft_Kronemberger Lopes_857637204-image22.png|534px]] </div>
629
630
<div id="_Ref52371035" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
631
''Fig. 22 – Consolidation ratio over time in Soil A.''</div>
632
633
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.
634
635
[[#_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.
636
637
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
638
 [[Image:Draft_Kronemberger Lopes_857637204-image23.png|486px]] </div>
639
640
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
641
''Fig.  23 – Consolidation ratios for different permeability coefficients at 7.5m from the top of the anchor.''</div>
642
643
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.
644
645
<div id="_Ref34133953" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
646
''Table 3 – Time (in days) needed to achieve different consolidation ratios: 7.5 m above the top of the anchor.''</div>
647
648
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
649
|-
650
|  rowspan='2' style="border-top: 1pt solid black;text-align: center;"|<math display="inline">\mathit{\boldsymbol{k}}</math> (m/dia)
651
|  colspan='6'  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|Consolidation ratio, <math display="inline">\mathit{\boldsymbol{U}}</math>
652
|-
653
|  style="text-align: center;vertical-align: top;"|'''0.10'''
654
|  style="text-align: center;vertical-align: top;"|'''0.25'''
655
|  style="text-align: center;vertical-align: top;"|'''0.50'''
656
|  style="text-align: center;vertical-align: top;"|'''0.75'''
657
|  style="text-align: center;vertical-align: top;"|'''0.90'''
658
|  style="text-align: center;vertical-align: top;"|'''0.99'''
659
|-
660
|  style="text-align: center;vertical-align: top;"|2.50E-04
661
|  style="text-align: center;vertical-align: top;"|3.00E-02
662
|  style="text-align: center;vertical-align: top;"|3.30E-01
663
|  style="text-align: center;vertical-align: top;"|3.79E+00
664
|  style="text-align: center;vertical-align: top;"|6.00E+01
665
|  style="text-align: center;vertical-align: top;"|3.08E+02
666
|  style="text-align: center;vertical-align: top;"|1.41E+03
667
|-
668
|  style="text-align: center;vertical-align: top;"|2.50E-05
669
|  style="text-align: center;vertical-align: top;"|3.10E-01
670
|  style="text-align: center;vertical-align: top;"|3.26E+00
671
|  style="text-align: center;vertical-align: top;"|3.97E+01
672
|  style="text-align: center;vertical-align: top;"|6.23E+02
673
|  style="text-align: center;vertical-align: top;"|3.09E+03
674
|  style="text-align: center;vertical-align: top;"|1.52E+04
675
|-
676
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.50E-06
677
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3.09E+00
678
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3.18E+01
679
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4.04E+02
680
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|6.32E+03
681
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|3.10E+04
682
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1.12E+05
683
|}
684
685
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.
686
687
==4.2.3. Comparison between FEM and analytical methodologies==
688
689
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.
690
691
<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.
692
693
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
694
 [[Image:Draft_Kronemberger Lopes_857637204-image24.png|450px]] </div>
695
696
<div id="_Ref52371266" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
697
''Fig.  24 – Comparison between total radial stresses – Soils A and B.''</div>
698
699
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.
700
701
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.
702
703
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.
704
705
<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).
706
707
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.
708
709
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.
710
711
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
712
 [[Image:Draft_Kronemberger Lopes_857637204-image25.png|600px]] </div>
713
714
<div id="_Ref52371311" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
715
''Fig.  25 – Comparison between methodologies – Soil A.''</div>
716
717
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
718
 [[Image:Draft_Kronemberger Lopes_857637204-image26.png|600px]] </div>
719
720
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
721
''Fig. 26 – Comparison between methodologies – Soil B.''</div>
722
723
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
724
 [[Image:Draft_Kronemberger Lopes_857637204-image27.png|600px]] </div>
725
726
<div id="_Ref52371314" class="center" style="width: auto; margin-left: auto; margin-right: auto;">
727
''Fig.  27 – Comparison between methodologies – Soil C.''</div>
728
729
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.
730
731
<span id='_GoBack'></span>Finally, it is emphasized that the FE analyses presented a setup curve with a predominantly linear behavior with the ''log'' of time. Among the empirical methodologies, the Svinkin and Skov <span id='cite-_Ref534917465'></span>[[#_Ref534917465|[19]]] empirical formulation, Eq. <span id='cite-_Ref52376819'></span>[[#_Ref52376819|(16)]], presented a better match to those obtained with the FE approach. This indicates that a reasonable approach to obtain the setup curve is to conduct a pair of analyses with the FE approach, ''e. g.'', for consolidation times of zero (''i. e.'', immediately after the installation of the anchor) and several days after the installation (typically, 5 000 or 10 000 days). Then, the coefficient ''B'', Eq. <span id='cite-_Ref52376819'></span>[[#_Ref52376819|(16)]], can be calibrated using these two points of the curve. In this study, a value of 0.2 for the coefficient ''B'' was found adequate to simulate the setup for all soils considered.
732
733
==5. Conclusions==
734
735
This study presented a numerical approach based on axisymmetric FE models to evaluate the behavior of a finless torpedo anchor submitted to axial loads at different instants of time after its installation.
736
737
The proposed approach was employed in a parametric study aiming to evaluate the setup effect on the holding capacity of the finless torpedo anchor considering different undrained shear strength profiles, elasticity moduli, and permeability coefficients of the soil.
738
739
The parametric study evidenced the importance of the soil undrained shear strength and the permeability coefficient in estimating the holding capacity of the torpedo anchor and in the total time of setup. Soils with higher undrained shear strengths present states of effective stress with higher magnitudes and, consequently, at the end of the setup process, higher pullout capacities. It was also identified that higher undrained shear strengths are associated with higher rates of increase of soil consolidation at the first instants after installation. Regarding the anchor's ultimate capacity increase, a higher ratio was observed between the initial and ultimate capacities considering soils with higher undrained shear strengths.
740
741
When varying the permeability coefficient of the soil, an inverse proportional relation between this coefficient and the consolidation time was observed. For instance, if the permeability coefficient increases an order of magnitude, the consolidation time decreases in the same proportion, and vice-versa.
742
743
The empirical formulations presented in this study proved their potential to estimate the holding capacity of the finless torpedo anchors with good agreement to the FE approach. However, each of the empirical formulations fits better for different consolidation times after installation. Skov and Denver <span id="cite-_Ref534917435"></span>[[#_Ref534917435|[18]]] formulation presented an excellent agreement for the first ten days after installation. ABS <span id="cite-_Ref51515640"></span>[[#_Ref51515640|[20]]] formulation and API <span id="cite-_Ref502324791"></span>[[#_Ref502324791|[21]]] adapted formulation agreed better for consolidation times over 100 days after the anchor installation. The Svinkin and Skov <span id="cite-_Ref534917465"></span>[[#_Ref534917465|[19]]] formulation agreed well with the FE results. Moreover, with just a pair of numerical analyses, ''i.e''., immediately after the anchor installation and the full consolidated soil, for instance, it is possible to calibrate the Svinkin and Skov <span id="cite-_Ref534917465"></span>[[#_Ref534917465|[19]]] formulation and obtain a reasonable estimation of the anchor setup process with a relatively low computational effort.
744
745
Despite the insights obtained in the analyses of finless torpedo anchors, which may help design fixed points to anchor flexible lines, monobuoys or floating production systems in shallow waters, some important aspects were not addressed in this work. It is important to emphasize that the FE analyses did not consider thixotropic effects, and further investigations over this phenomenon are recommended. Thixotropy may rule the soil regain in the first instants of time after installation. Moreover, the effect of the setup when inclined loads are imposed on the anchor is another aspect that deserves future investigation. However, Sousa ''et al.'' [4] showed that, for inclinations between 40 degrees and 90 degrees (vertical), the holding capacity of the anchor can be obtained by simply projecting the pullout load in the considered direction. Finally, this study concentrated on finless torpedo anchors, but several torpedo anchors have fins that may also affect the soil stress state during installation. The effect of the fins has also to be considered in future studies.
746
747
==6. Acknowledgments==
748
749
Authors also thank the National Council for Scientific and Technological Development (CNPq) for the research grants PQ-311084/2019-2 and PQ-306357/2018-6, the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior- Brasil (CAPES), Finance Code 001, and the Fundação Carlos Chagas de Amparo à Pesquisa do Estado do Rio de Janeiro (FAPERJ).
750
751
==REFERENCES==
752
753
<span id='_Ref31803537'></span><span id='_Ref2530944'></span>[[#cite-_Ref2530944|[1].]] Y. H. Kim, M. S. Hossain, D. Wang and M. F. Randolph. “Numerical investigation of dynamic installation of torpedo anchors in clay”. Ocean Engineering 108 (2015) 820-832.
754
755
<span id='_Ref31803040'></span>[[#cite-_Ref31803040|[2].]] Y. Zhao and H. Liu. “Toward a quick evaluation of the performance of gravity installed anchors in clay: Penetration and keying”. Applied Ocean Research 69 (2017) 148-159.
756
757
<span id='_Ref31817098'></span>[[#cite-_Ref31817098|[3].]] J. T. Lieng, A. Kavli, F. Hove and T. I. Tjelta, “Installation of Two Prototype Deep Penetrating Anchors at the Gjoa Field in the North Sea”, in: Proceedings of the Offshore Technology Conference, Houston, Texas, USA, May 3-6, 2010.
758
759
<span id='_Ref31804075'></span>[[#cite-_Ref31804075|[4].]] J. R. M. de Sousa, C. S. de Aguiar, G. B. Ellwanger, E. C. Porto, D. Foppa and C. J. Medeiros Jr. “Undrained Load Capacity of Torpedo Anchors Embedded in Cohesive Soils”. Journal of Offshore Mechanics and Arctic Engineering 133 (2) (2011).
760
761
<span id='_Ref31804736'></span>[[#cite-_Ref31804736|[5].]] C. Han, J. Liu, Y. Zhang and W. Zhao. “An innovative booster for dynamic installation of OMNI-Max anchors in clay: Physical modeling”. Ocean Engineering 171 (2019) 345-360.
762
763
<span id='_Ref31808849'></span>[[#cite-_Ref31808849|[6].]] Y. Tong, C. Han and J. Liu. “An innovative lightweight gravity installed plate anchor and its keying properties in clay”. Applied Ocean Research 94 (2020) 101974.
764
765
<span id='_Ref31810175'></span>[[#cite-_Ref31810175|[7].]] Y. Zhao, Y. Kim, M. S. Hossain, M. Nazem, J. Liu and Y. Hu. “Numerical advancements on the analysis of dynamically installed anchors”. Ocean Engineering 182 (2019) 343-359.
766
767
<span id='_Ref31892302'></span>[[#cite-_Ref31892302|[8].]] M. S. Hossain, Y. Kim and C. Gaudin. “Experimental investigation of installation and pullout of dynamically penetrating anchors in clay and silt”. Journal of Geotechnical and Environmental Engineering 140 (7) (2014) 04014026.
768
769
<span id='_Ref502687747'></span>[[#cite-_Ref502687747|[9].]] K. Bergset. “Radial Consolidation of Pore Pressure Induced by Pile Driving”. Master Thesis, Norwegian University of Science and Technology, Trondheim, Norway (2015).
770
771
<span id='_Ref502153693'></span>[[#cite-_Ref502153693|[10].]] V. E. Komurka and A. B. Wagner. “Estimating Soil/Pile Set-up”. Final Report, University of Wisconsin-Madison (2003).
772
773
<span id='_Ref502153701'></span>[[#cite-_Ref502153701|[11].]] Simulia. “Analysis of Driven Pile Setup with Abaqus/Standard”. Abaqus Technology Brief (2007).
774
775
<span id='_Ref31897166'></span><span id='_Ref2531416'></span>[[#cite-_Ref2531416|[12].]] M. D. Richardson, C. D. O’Loughlin, M. F. Randolph and C. Gaudin. “Setup following installation of dynamic anchors in normally consolidated clay”. Journal of Geotechnical and Environmental Engineering 135 (4) (2009) 487-496.
776
777
<span id='_Ref31899459'></span>[[#cite-_Ref31899459|[13].]] E. Radgahar, M. S. Raie and N. Motamani. “Simulation of torpedo-shaped (without fins) for semisubmersible offshore platforms under tensile force”. Journal of Applied Environmental and Biological Sciences (5) (2015) 107-111.
778
779
<span id='_Ref31899465'></span>[[#cite-_Ref31899465|[14].]] M. S. Raie and J. L. Tassoulas. “Simulation of torpedo anchor set-up”. Marine structures 49 (2016) 138-147.
780
781
<span id='_Ref31900575'></span>[[#cite-_Ref31900575|[15].]] H. Sabetamal, J. P. Carter, M. Nazem and S. W. Sloan. “Coupled analysis of dynamically penetrating anchors”. Computers and Geotechnics 77 (2016) 24-44.
782
783
<span id='_Ref32417380'></span>[[#cite-_Ref32417380|[16].]] C. J. Medeiros Jr. “Low cost anchor system for flexible risers in deep Waters”. Offshore Technology Conference, Houston, Texas, U.S.A., 6-9 May 2002.
784
785
<span id='_Ref534978548'></span>[[#cite-_Ref534978548|[17].]] Abaqus/CAE v. 6.13-1. © Dassault Systèmes Simulia Corp., 2013.
786
787
<span id='_Ref534917435'></span>[[#cite-_Ref534917435|[18].]] R. Skov and H. Denver. “Time-dependence of bearing capacity of piles”, in: Proceedings of the 3rd International Conference on Application of Stress-Waves Theory to Piles, Ottawa, May 25-27, 1988.
788
789
<span id='_Ref534917465'></span>[[#cite-_Ref534917465|[19].]] M. R. Svinkin and R. Skov. “Set-up effects of cohesive soils in pile capacity”, in: Proceedings of the 6th International Conference on Application of Stress Waves to Piles, Sao Paulo, Brazil, September 11-13, 2000.
790
791
<span id='_Ref51515640'></span>[[#cite-_Ref51515640|[20].]] ABS. “Guidance notes on design and installation of dynamically installed piles”, American Bureau of Shipping, USA, 2017.
792
793
<span id='_Ref502324791'></span>[[#cite-_Ref502324791|[21].]] API. “Recommended Practice for Planning, Designing and Constructing Fixed Offshore Platforms – Working Stress Design (RP 2A-WSD)”, 20th ed., American Petroleum Institute, USA, 2005.
794
795
<span id='_Ref503887814'></span>[[#cite-_Ref503887814|[22].]] S. Helwany, “Applied Soil Mechanics with ABAQUS Applications”, John Wiley & Sons, Inc., Hoboken, New Jersey, 2007.
796
797
<span id='_Ref535247756'></span>[[#cite-_Ref535247756|[23].]] K. Karlsrud, “Prediction of load-displacement behavior and capacity of axially loaded piles in clay based on analyses and interpretation of pile load test results”. Doctoral thesis at NTNU, Trodheim, 2012.
798
799
<span id='_Ref534824256'></span>[[#cite-_Ref534824256|[24].]] M. F. Randolph and C. P. Wroth. “Anaytical Solution for the Consolidation around a Driven Pile”. International Journal for Numerical and Analytical Methods in Geomechanics, vol. 3, (1979) 217-229.
800
801
<span id='_Ref534824674'></span>[[#cite-_Ref534824674|[25].]] R. Hill. “The mathematical theory of plasticity”. Oxford University Press, 1950.
802
803
<span id='_Ref534824683'></span>[[#cite-_Ref534824683|[26].]] R. Gibson and W. Anderson, “In situ measurement of soil properties with the pressumeter”. Civil engineering and public works review, vol. 56 (1961) 615-618.
804
805
<span id='_Ref2592654'></span>[[#cite-_Ref2592654|[27].]] J. M. Keaveny and J. K. Mitchell. “Strength of fine-grained soils using the piezocone”. Use of In-Situ Tests in Geotechnical Engineering (GSP 6), ASCE, Reston/VA (1986) 668-699.
806
807
<span id='_Ref2698172'></span>[[#cite-_Ref2698172|[28].]] J. P. Carter, M. F. Randolph and C. P. Wroth. “Stress and pore pressure changes in clay during and after the expansion of a cylindrical cavity”. International Journal for Numerical and Analytical Methods in Geomechanics (3) (1979) 305-322.
808
809
<span id='_Ref32331392'></span>[[#cite-_Ref32331392|[29].]] G. K. Lopes. “Análise numérica da capacidade de carga de âncoras torpedo considerando efeitos de setup”. Master thesis at Federal University of Rio de Janeiro, Rio de Janeiro, Brazil, 2019.
810
811
<span id='_Ref32496298'></span>[[#cite-_Ref32496298|[30].]] C. D. O’Loughlin, M. F. Randolph and M. Richardson. “Experimental and theoretical studies of deep penetration anchors”. Offshore Technology Conference, Houston, Texas, U.S.A., 3-6 May 2004.
812
813
<span id='_Ref32496303'></span>[[#cite-_Ref32496303|[31].]] D. G. True. “Undrained vertical penetration into ocean bottom soils”. Ph.D. Thesis, University of California, Berkeley, California, 1976.
814

Return to Lopes et al 2020a.

Back to Top