Research Article  Open Access
Luís V. S. Sagrilo, José Renato M. de Sousa, Edison C. P. Lima, Elisabeth C. Porto, Jane V. V. Fernandes, "A Study on the Holding Capacity Safety Factors for Torpedo Anchors", Journal of Applied Mathematics, vol. 2012, Article ID 102618, 18 pages, 2012. https://doi.org/10.1155/2012/102618
A Study on the Holding Capacity Safety Factors for Torpedo Anchors
Abstract
The use of powerful numerical tools based on the finiteelement method has been improving the prediction of the holding capacity of fixed anchors employed by the offshore oil industry. One of the main achievements of these tools is the reduction of the uncertainty related to the holding capacity calculation of these anchors. Therefore, it is also possible to reduce the values of the associated design safety factors, which have been calibrated relying on models with higher uncertainty, without impairing the original level of structural safety. This paper presents a study on the calibration of reliabilitybased safety factors for the design of torpedo anchors considering the statistical model uncertainty evaluated using results from experimental tests and their correspondent finiteelementbased numerical predictions. Both working stress design (WSD) and load and resistance factors design (LRFD) design methodologies are investigated. Considering the WSD design methodology, the single safety is considerably lower than the value typically employed in the design of torpedo anchors. Moreover, a LRFD design code format for torpedo anchors is more appropriate since it leads to designs having lessscattered safety levels around the target value.
1. Introduction
Torpedo anchors, as shown in Figure 1, are a costeffective and practical solution for anchoring taut mooring lines of floating units as shown in Figure 2. This anchor has a “rocket” shape and is embedded into the seabed by free fall using its own weight as driven energy. The structure of a typical torpedo anchor usually comprises four different components [1, 2]: a padeye that connects the first chain segment of the mooring line to the anchor; a ballasted shaft of carbon steel; four flukes, which are approximately rectangular (torpedo anchors with no flukes may also be employed); a conical tip that is designed to help the embedment of the anchor. Typical torpedo anchors weight varies from 350 kN up to 980 kN.
(a)
(b)
These anchors have been first employed to hold flexible risers just after their touchdown points (TDPs) in order to avoid that the loads generated by the floating unit movements were transferred to subsea equipments, such as wet Christmas trees [1]. In 2003 a high holding capacity torpedo anchor was designed to sustain the loads imposed by the mooring lines of a floating production, storage, and offloading (FPSO) [3] and, since then, various other floating production systems have been anchored using the torpedo concept [4]. Furthermore, recently, Henriques Jr. et al. [4] presented a new torpedo anchor with a total weight of 1200 kN and larger flukes that will be employed to moor FPSOs in ultradeep waters. These characteristics, altogether, significantly increase its hold capacity when compared to the traditional anchors with 980 kN in weight.
Due to their “rocket” shape, the load position and inclination with respect to their flukes, traditional methods, such as conventional beamcolumn analysis procedures (otherwise known as the PY method) or the simple formulae proposed by API guidelines [5] for predicting the axial capacity of slender piles, often cannot be directly employed to estimate the holding capacity of these anchors.
However, numerical procedures, based on the finite element (FE) method, have recently been proposed by Brandão et al. [2] and de Sousa et al. [6]), and experimental tests performed by PETROBRAS have demonstrated their good accuracy [7]. An important outcome of these FEbased models is the lower uncertainty in their predictions when compared to traditional methods. Therefore, it is possible to reduce the design safety factors values that have been calibrated relying on prediction models with higher uncertainty, without impairing the original level of structural safety.
There is no specific guideline for the design of torpedo anchors. Therefore, it is common practice to follow the indications from API RP 2A [5] for piles foundations to support offshore structures. This guideline prescribes that the safety factor associated with these piles, according to a WSD design methodology, under operating environmental conditions and during producing operations is equal to 2.0. It seems that this relatively large safety factor is mainly related to the high prediction model uncertainty due to the large number of pile capacity prediction models available at the moment of the calibration process.
This paper presents a study on the calibration of safety factors for the ultimate limit state design of torpedo anchors. One important aspect of this study is the availability of results from holding capacity tests of six torpedo anchors installed in Campos Basin, offshore Brazil [7]. These results made it possible to assess the model uncertainty statistics associated with the FEbased model proposed by de Sousa et al. [6].
Both working stress design (WSD) and load and resistance factors design (LRFD) methodologies are investigated. Concerning the traditional WSD methodology, it is shown that, for the same actual design safety level of the traditional offshore piles, its single safety factor can be significantly lowered. However, the use of WSD design methodology results in designs with very scattered safety levels depending on the ratio of the functional and environmental load actions. Aiming at overcoming this drawback, a LRFD methodology calibration is also investigated in this paper. The results show that the structural safety levels of LRFDbased designs are more uniform than the WSDbased ones.
In the following sections, firstly, a brief review of the FE model proposed by de Sousa et al. [6] is presented followed by a comparison between the theoretical and experimental holding capacities. After that, the reliability baseddesign study is presented, and, finally, the main conclusions of this work are stated.
2. FiniteElement (FE) Model
2.1. Overview
Recently, de Sousa et al. [6] have proposed a FE model devoted to predict the undrained load capacity of torpedo anchors. This model employs isoparametric solid finite elements to model both the soil medium and the anchor. These elements are capable of representing the physical nonlinear behavior of the soil medium and of the anchor itself. Large deformations may also be considered. The soilanchor interaction is ensured by surfacetosurface contact elements placed on the external surface of the anchor and the surrounding soil medium. A general overview of this model is presented in Figure 3. Next, its main aspects are briefly reviewed.
(a)
(b)
2.2. Soil Modeling
2.2.1. Constitutive Matrix
The soil is assumed to be a perfectly elastoplastic isotropic material with physical properties variable with depth. Hence, the relation between stresses and strains is given in the following form:
The constitutive elastoplastic matrix, , is given by [9] as follows:
The elastic total stress constitutive matrix can be split into two parts, as follows [9]:
In order to represent the nonlinear material behavior of the soil, the DruckerPrager model was chosen. This model approximates the irregular hexagon of the MohrCoulomb failure surface by a circle, and, consequently, the DruckerPrager yield function is a cylindrical cone. The yield function, , and plastic potential function, , for this model are given by [10] as follows: where , , , and are stress state parameters. These parameters depend on the adopted approximation to the MohrCoulomb hexagon and are functions of the internal friction angle, , the dilation angle, , and the cohesion, , of the soil. Figure 4 presents three possible approximations. In the proposed model, the circumscribed DruckerPrager approximation is used. It is assumed that the DruckerPrager circle touches the MohrCoulomb hexagon at a Lode angle of +30°, which leads to the parameters [9] as follows:
It is worth mentioning that the rules for designing deep offshore foundations, such as the API RP 2A [5], typically assume that loading occurs rapidly enough so that no drainage and hence no dissipation of excess pore pressures occur. This leads to an undrained condition and, as a consequence, the undrained shear strength of the soil, , may replace the cohesion, , in (2.5).
On the other hand, if the rate of loading is slow enough such that no excess pore pressures are developed and sufficient time has elapsed since any previous application of stresses such that all excess pore pressures have been dissipated, a drained analysis would have to be performed. Lieng et al. [11], however, indicate that the undrained loading capacity of a deep penetrating anchor, such as a torpedo anchor, is quite close to the drained capacity, and this is also assumed in this work. Nevertheless, this aspect has yet to be validated for torpedo anchors.
2.2.2. FE Mesh Characteristics
The undrained response of a torpedo anchor embedded in clay is a classical elastoplastic problem that requires no locking for incompressible materials and reasonable bending behavior in order to obtain acceptable answers. Moreover, as high plastic strains are expected in the soil, the elements may become highly distorted, and insensitivity to these distortions is demanded. Finally, coarse meshes are often required for computational efficiency, and elements capable of providing sufficient accuracy are also needed. Aiming at addressing all these requirements, a possible choice, according to Wriggers and Korelc [12], is the use of enhanced strain hexahedrical and prismatic isoparametric finite elements with three degrees of freedom at each node: translations in directions , , and . These elements avoid both bend and volumetric locking by adding 13 internal degrees of freedom to the element [13, 14].
An outline of the main dimensions of the soil mesh is shown in Figure 3. The proposed mesh is a cylinder with a base diameter of , where is the pile diameter. The height of the cylinder is given by the sum of the penetration of the torpedo anchor, , the length of the anchor, , and the distance of the tip of the torpedo anchor to the bottom of the FE mesh, .
Each “slice” of the cylinder has its own physical properties, and, consequently, variable strength, density, longitudinal and transverse moduli, and Poisson coefficients may be assumed.
The vertical displacement of the cylinder is restrained at the nodes associated with its base. The displacements in and directions (see Figure 3) of nodes in the outer wall of the cylinder are restrained. Once a plane of symmetry is found, the outofplane displacements are also restrained. In order to avoid any influence of these boundary conditions on the response of the anchor, several mesh tests were performed. In these tests, the cylinder diameter and the distance from the tip of the anchor to the bottom of the mesh () were varied. It was observed that a diameter of for the cylinder and a distance of 5.0 m for was enough to simulate an “infinite” media.
2.3. Anchor Modeling
The torpedo anchor is modeled with eight node isoparametric solid elements analogous to those used in soil representation. These elements, as pointed out before, are capable of considering both material and geometrical nonlinearities. A typical FE mesh of a torpedo anchor is shown in Figure 5.
(a)
(b)
In order to properly represent the bending of the flukes and the shaft and, consequently, avoid shear locking problems, the numerical integration of each element is performed with the enhanced strain method.
It is worth mentioning that neither the padeye at the top of the anchor nor the mooring line is represented in the proposed model. The load from the mooring line is applied at the gravity center of the padeye, where a node is placed and rigidly connected to the top of the anchor by rigid bars, as presented in Figure 6.
2.4. AnchorSoil Interaction
The model proposed by de Sousa et al. [6] employs contact elements that allow relative large displacement and separation between the surfaces in contact. These elements are placed on the external surface of the torpedo anchor and on the surrounding soil. A contact detection algorithm based on the pinball technique and contact forces evaluated with the augmented Lagrangian method are employed in the formulation of these elements.
Relative sliding between the anchor and the soil is allowed when the shear stresses along the outer wall of the anchor exceed the values estimated with the MohrCoulomb friction model, which is given by [5] as follows:
The friction angle between the soil and the anchor wall, , is stated as [5]
The adhesion factor, , may be calculated by the following formula [5]: where
2.5. Solution Procedure and Implementation
The FE analysis is divided in three different steps. In the first one, the initial stress state, that is, the stresses in the soil prior to the imposition of any kind of structural load to the anchor, is generated. The model proposed by de Sousa et al. [6] does not simulate the anchor penetration in the soil, and it is assumed, by hypothesis, that the anchor is loaded with the soil stress state undisturbed by the installation process. The anchor is thereby supposed to be “wished in place,” and this undisturbed stress state is imposed to the soil in the first step of the analysis.
However, the undrained shear strength of the soil and its stress state are deeply affected by the dynamic installation process and the setup time for reconsolidation or recovery of shear strength after installation [15]. The direct consideration of this modified stress state and undrained shear strength is quite complex and, therefore, the hypothesis of an undisturbed stress state, as employed in this work, is commonly assumed in the design of different types of offshore anchors such as suction anchors or other types of dynamic penetrating anchors [16–18]. However, the validity of this hypothesis is subject of intense debate. Lieng et al. [11], relying on FE simulations, postulate that reconsolidation after installation of a typical dynamic anchor would occur in less than six months, and, as a consequence, a holding capacity determined from undrained analyses with an adhesion factor close to 1.0 is expected to be valid. Richardson et al. [19], relying on centrifuge tests and theoretical models, indicate that 50% of the operational capacity of a dynamic anchor is attained 35 up to 350 days after its installation. According to the same authors, 90% of the operational capacity is reached 2.4 up to 24 years after the installation of the anchor. These values are extremely dependent on the horizontal coefficient of consolidation, and, therefore, authors recommend that dynamic anchors should be installed in soils that consolidation proceeds quickly. Furthermore, Richardson et al. [19] also indicate that values higher than 1.0 for the adhesion factor are reached in the long term leading to the conclusion that the undrained holding capacity of the anchor would increase with time, and this capacity would even be higher than the one estimated with an undisturbed stress state. Mirza [20] state that the time required by offshore piles to reach their fulloperational capacity varies between 1 and 2 years, and Jeanjean [21] indicates that 90% of the operational capacity of suction anchors installed in soft Gulf of Mexico clays is attained in 90 days or less.
A typical torpedo anchor is firstly loaded approximately three months after its installation, and, therefore, a great part of its full operational capacity has already been reached. Moreover, the initial imposed loads are normally much lower than its actual load capacity, as its design relies on the maximum load obtained in mooring analyses that consider both intact and damaged (one mooring line broken) conditions in the mooring system and also several environmental load cases. These load cases encompass, for instance, extreme waves and currents. Hence, as the objective of the proposed FE model is to evaluate the undrained operational holding capacity of the anchor, it seems to be reasonable to assume that full reconsolidation is achieved before the maximum operational load is imposed to the anchor. Nevertheless, in order to improve the proposed design method, more experimental and also theoretical studies on this aspect have to be performed.
In the second step of the analysis, the self weight of the anchor is imposed. The analysis is restarted with the stress state from the first step, and gravity acts on the anchor. Finally, in the third step, the total load is applied with adaptable increments since, as the surrounding soil progressively fails, the stability of the solution is affected, and lower load increments are required. Typically, this last load step starts with an initial load increment of 5% of the total load applied, and as the analysis progresses, it can be reduced to 1% of the total load.
At each load step, convergence is achieved if the Euclidian norm (L2) of the residual forces and moments is less than 0.1% of the absolute value of the total forces and moments applied.
The model was implemented in a program called ESTACAS. This software generates FE meshes to be analyzed with ANSYS [22] program, where the following finite elements are employed: SOLID185 in order to model the soil and the torpedo anchor; CONTA174 and TARGE170 are used to simulate the contact between the soil and the anchor.
3. Comparison with Experimental Tests
Aiming at validating the torpedo anchor geotechnical and structural design methodology, PETROBRAS conducted six fullscale tests in Campos Basin, offshore Brazil [7]. These six fullscale tests were divided into two sets.
In the first set of tests, three torpedoes (type 1), each one weighing 350 kN, were installed at a location with a water depth of 500 m. At this location, soil was characterized as being pure clay after cone penetration tests (CPTs) had been performed. The torpedoes embedment depths varied from 8 m up to 11 m. Besides, the torpedoes were pulled out with load inclinations of about 40° [7]. The second set of tests consisted of pulling out three torpedoes (type 2) weighting 430 kN each. The water depth at this second location is 150 m, and the soil comprises both clay and sand layers. The torpedoes final embedment depths, in this case, varied from 6 m to 8 m, and the load was imposed with inclinations varying from 50° to 67°.
Table 1 presents the relation between the experimental measured pullout loads and those estimated with the model proposed by de Sousa et al. [6].

Despite all uncertainties related to the model and the geotechnical parameters involved, Table 1 indicates that the pullout loads estimated with the FE model agree quite well to the measured loads. Hence, these results will serve as a basis to the reliabilitybased design that is presented next.
4. ReliabilityBased Design
Modern probabilistic design codes aim at designing new structures with a specified target probability of failure. This can be achieved by fully probabilistic or deterministic codes [23]. The latter is more common in engineering design. In these codes the target probability of failure is not stated but it is implicitly achieved by means of safety factors. These safety factors are calibrated by using standard structural reliability methods [23, 24], in which all random variables (and their uncertainties) are represented by means of probability distributions.
The two main deterministic design methodologies for ultimate limit state (ULS) are the working stress design (WSD) and load and resistance factor design (LRFD). Focusing specifically on the case of torpedo anchors for mooring systems (see Figure 6) under the WSD methodology, the design equation can be written as where the characteristic functional load effect, , is associated with the initial mooring line pretension applied to the anchor, and the characteristic environmental load effect, , corresponds to the total tension acting on the anchor, when environmental loads are applied, minus the functional one. According to this classification, one can observe that the functional load effect has smaller uncertainties than the environmental one.
In the case of the LRFD methodology, the design equation can be written as where and are the partial safety factors associated with functional and environmental load effects, respectively.
The safety factors of (4.1) and (4.2) are obtained by means of a calibration process. The main idea of the calibration process is to obtain safety factors that guarantee a certain target failure probability, , for the design. The target probability of failure may be stated as
The calibration process is undertaken by means of reliability analyses considering various design cases with different ratios and also different statistical descriptions of loadings and anchor resistances that are expected to be found in practice. Then, the safety factors can be defined, for instance, for LRFD methodology (similar procedure applies to the WSD methodology) through an optimization process to solve the following problem: where is the reliability index of a given torpedo anchor design case that can be obtained by a reliability analysis method, for example, firstorder reliability method (FORM), secondorder reliability method (SORM), or Monte Carlo simulation [23, 24]. In order to pursue this kind of analysis, firstly, the limit state function must be defined. For the present study, it can be stated as where is the vector containing the random variables associated with the torpedo anchor holding capacity , functional load , and environmental load . The variable can be represented by where is the predicted holding capacity evaluated by a given methodology of analysis, and is the modeling uncertainty that represents the statistical bias between the measured and predicted pilebearing capacity. Therefore, some experimental data is necessary to statistically define this random variable. In fact, takes into account all uncertainties regarding the soil parameters, model representation, and so forth.
In this work, all reliability analyses have been performed with Monte Carlo simulation approach which is applied to solve the following multidimensional integral: where is the failure domain, that is, , and is the joint probability distribution of the random variables for the case under consideration.
4.1. Safety Factor Calibration
4.1.1. Random Variables Modeling
Tests with six torpedo anchors installed in clay and claysand soils have recently been performed [7]. The measured torpedo anchor holding capacities has been divided by their corresponding predicted numerical estimation, , according to the model presented before, to determine the model uncertainty, , statistics. The values obtained are presented in Table 1. The mean, , and standard deviation, , for according to Table 1 data correspond to 0.993 and 0.118, respectively. In order to properly take into account the limited number of experiments () in the reliability analysis, it is necessary to apply some statistical modeling technique for reduced sample data as explained in Ditlevsen and Madsen [25]. Assuming that the variable has a previous lognormal distribution, as indicated in Fenton [8], this variable should be modeled by the following Student’s distribution to account for the number of experiments: where
Then, considering Table 1 data, the mean and standard deviation of this Student’s distribution becomes
It is important to notice that, due to the limited data sample, is 64% higher than the sample standard deviation, . These values also show that the FE model for the torpedoholding capacity has a lower uncertainty when compared, for instance, with prediction methods reported in an APIsupported study for axial bearing capacity of single piles in clays [8]. In this latter case, the mean and standard deviation for the modeling uncertainty variable has been, respectively, 1.04 and 0.34.
In the reliability analysis of marine structures, the functional loading is usually modeled by a normal distribution with a coefficient of variation (CoV or ) in the order of 0.05–0.10. In this work, a CoV () of 0.07 is assumed. The environmental load effect is usually represented by its annual extreme value distribution which, in most of the cases, corresponds to a Gumbel distribution. Considering the inherent random variability of the environmental parameters related to waves, wind and current, and also modeling uncertainties in the estimation of this complex load effect [26], it is reasonable to assume a CoV () in the order of 0.10–0.30 for this distribution.
4.1.2. Characteristic Design Parameters
The characteristic design parameters used in (4.1) and (4.2) must be clearly defined in any calibration process or design code. Different values for these parameters are also related to different sets of safety factors.
In the present study, the characteristic value of the torpedo anchorbearing capacity, , is defined as that predicted by the FEbased model previously described [6]. The characteristic value of the functional load effect, , has been assumed as its estimated mean value, that is, . In the case of the environmental load effect, following the traditional design practice in offshore engineering, its characteristic value has been defined as the most probable value of the 100yr extreme distribution. This value corresponds to that related to 99% fractile of the annual extreme cumulative distribution [27].
In the calibration process, a broader range of location water depths for the floating units has been assumed by considering the ratio to be between 1.0 and 5.5. The lower ratio is associated with deeper waters in which the functional load at the anchor can be of the same order as the environmental one. The converse applies to the higher ratio. These ratios account for actual floating systems installed in water depths ranging approximately from 200 m to 2500 m [28].
4.1.3. Target Safety Index
In order to obtain the target safety index for the calibration process, firstly, all cases investigated have been designed according to the API guidelines [5] considering an intact mooring system, that is, using the WSD methodology with a safety factor equals to 2.0. Secondly, the reliability index of each case has been evaluated considering the uncertainty modeling mean and standard deviation parameters as 1.04 and 0.34, respectively, taken from Fenton [8]. Finally, the target reliability index has been taken as the average of all these reliability indexes. Considering specifically cases for , the target reliability index, , has been found to be approximately 2.9. Figure 7 illustrates the results that have been obtained. In summary, the target reliability index assumed in this work is more or less the same that is implicitly assumed when a single pile in clay is designed to support an axial loading.
4.1.4. WSD Calibration
The optimized WSD safety factor, , for the holding capacity design of torpedo anchors (see (4.1)), considering the previously discussed FE model uncertainty, has been found to be 1.45. Figure 8 shows a comparison between the target reliability index and the reliability indexes of all design cases investigated considering the optimized safety factor. When compared to the safety factor of 2.0, which is usually employed by the offshore oil industry, the modeling uncertainty reduction obtained with the FE model causes a significant reduction in the safety factor.
Comparing Figures 7 and 8, it is important to notice that a larger model uncertainty also results in more scattered reliability indexes among the design cases.
4.1.5. LRFD Calibration
As it is well known in the public literature [29], the use of a single safety factor in the WSD methodology leads to more scattered designs concerning their safety levels. The LRFD methodology, which attributes partial safety factors for different loading sources (see (4.2)), is an alternative to reduce this scatter.
In order to obtain lessscattered reliability indices for design cases considered in the calibration process, two sets of safety factors have been identified. The LRFD methodology design check for the holding capacity of torpedo anchors will be read as where
The two sets of safety factors in (4.12), and , have been obtained by applying the optimization process (4.4) twice: firstly, for all design cases in the range and, after that, for the other design cases in the range .
A comparison between the target reliability index and the reliability indexes of all cases investigated considering the LRFD design, (4.4), is shown in Figure 9.
For both WSD and LRFD design methodologies, the average of safety index values has been the target reliability index ; however, their coefficients of variation have been, respectively, 4.5% and 2.5%. This shows that LRFD methodology gives designs with safety levels that are slightly less scattered around the target value established for the calibration process.
5. Final Remarks
This paper has presented a reliabilitybased design study for torpedo anchors considering that their ultimate bearing capacity is predicted numerically by a nonlinear finite element model that represents the soil medium and the anchor itself. By considering the results of six fullscale experimental pull out tests, it was observed that the numerical model is practically unbiased and presents a relatively low level of statistical uncertainty. Using the model uncertainty statistics in reliabilitybased calibration process, it was shown that the same level of structural safety implied in the traditional design of offshore piles can be achieved with the use of lower safety factors. For instance, considering the WSD design methodology, the single safety factor has a significant reduction dropping from 2.0 to 1.45.
This calibration study has also shown that a LRFD design code format for torpedo anchors is the most appropriate one since it leads to designs having lessscattered safety levels around the target value. In this case, different safety factors are proposed for the environmental and functional parcels of the total load. These factors are also lower than the one proposed by API RP 2A [5].
An important point to be noticed is that it is very likely that the safety factors presented in this paper can become even lower when more experimental tests are performed since a great part of the model uncertainty considered in the present work comes from the very limited number of experiments available. More experiments are thereby required in order to estimate these new safety factors.
Nomenclature
:  Purely elastic total stress constitutive matrix 
:  Effective stress constitutive matrix 
:  Elastoplastic constitutive matrix of the soil 
:  Constitutive matrix of the fluid 
:  Cohesion of the soil 
:  Model uncertainty 
:  Diameter of the torpedo anchor including its flukes 
:  Environmental load effect applied to the anchor 
:  Yield function 
:  Limit state function 
:  Distance of the tip of the torpedo anchor to the bottom of the FE mesh 
:  Length of the torpedo anchor 
:  Penetration of the torpedo anchor into the soil (embedment depth) 
:  First invariant of the stress tensor 
:  Second invariant of the stress tensor 
k = :  Stress state parameters related to the yield function 
:  Coefficient of lateral earth pressure 
:  Functional load effect applied to the anchor 
LRFD:  Load and resistance factors design 
= :  Stress state parameters related to the plastic potential function 
:  Total number of design cases considered in the calibration process 
:  Plastic potential function 
:  Effective overburden pressure 
:  Target probability of failure 
:  Holding capacity of the anchor 
:  Undrained shear strength of the soil 
:  Safety factor 
WSD:  Working stress design 
:  Random variables vector 
:  Depth below mudline. 
Greek
:  Incremental total strain vector 
:  Incremental total stress vector 
:  Adhesion factor 
:  Reliability index of th case designed using the current set of safety factors 
:  Reliability index associated with 
:  Friction angle between the soil and the torpedo anchor wall 
:  Coefficient of variation 
:  Internal friction angle of the soil 
:  Standard normal cumulative distribution 
, :  Partial safety factors applied to functional and environmental load effects, respectively 
:  Mean 
:  Stress state in the soil 
:  Standard deviation 
:  Equivalent failure stress 
:  Limiting shear stress 
:  Dilation angle of the soil. 
Subscripts
:  Model uncertainty 
DP:  Related to the DruckerPrager criterion 
:  Environmental load effect 
:  Related to characteristic (or nominal) values 
:  Functional load effect 
:  Related to model predictions 
:  Target value. 
Acknowledgment
The authors would like to thank PETROBRAS Research Center for allowing them to publish the present paper.
References
 C. J. Medeiros Jr., “Low cost anchor system for flexible risers in deep waters,” in Offshore Technology Conference, pp. 1333–1337, Houston, Tex, USA, May 2002. View at: Google Scholar
 F. E. N. Brandão, C. C. D. Henriques, J. B. De Araújo, O. G. G. Ferreira, and C. Dos Santos Amaral, “Albacora Leste Field Development—FPSO P50 mooring system concept and installation,” in 38th Offshore Technology Conference, pp. 1443–1452, Houston, Tex, USA, May 2006. View at: Google Scholar
 J. B. De Araujo, R. D. Machado, and C. J. De Medeiros Junior, “High holding power torpedo pile—results for the first long term application,” in 23rd International Conference on Offshore Mechanics and Arctic Engineering, pp. 417–421, Vancouver, Canada, June 2004. View at: Google Scholar
 P. R. D. Henriques Jr., D. Foppa, E. C. Porto, and C. J. Medeiros Jr., “A New Torpedo Pile Conception for High Mooring Loads and Application in a Floating Production Unity in the PreSalt Area,” in Rio Oil & Gas Expo and Conference, Rio de Janeiro, Brazil, 2010. View at: Google Scholar
 API, Recommended Practice for Planning, Designing and Constructing Fixed Offshore Platforms – Working Stress Design (RP 2AWSD), American Petroleum Institute, 20th edition, 2005.
 J. R. M. de Sousa, C. S. de Aguiar, G. B. Ellwanger, E. C. Porto, D. Foppa, and C. J. de Medeiros, “Undrained load capacity of torpedo anchors embedded in cohesive soils,” Journal of Offshore Mechanics and Arctic Engineering, vol. 133, no. 2, Article ID 021102, 2010. View at: Publisher Site  Google Scholar
 E. C. Porto, C. S. Amaral, D. Foppa, J. V. V. Fernandes, and D. C. Corrêa, “Pullout Tests and Holding Capacity Evaluation of Torpedo Piles,” Tech. Rep. 004/2009, PETROBRAS, Rio de Janeiro, Brazil, 2009. View at: Google Scholar
 G. A. Fenton, “Probabilistic Methods in Geotechnical Engineering,” Short Course Notes, 1997, http://www.engmath.dal.ca/risk/publications.html. View at: Google Scholar
 D. M. Potts and L. Zdravkovic, Finite Element Analysis in Geotechnical Engineering – Theory, Thomas Telford, 1st edition, 1999.
 W. F. Chen and G. Y. Baladi, Soil Plasticity: Theory and Implementation, Elsevier, Amsterdam, The Netherlands, 1st edition, 1985.
 J. T. Lieng, A. Kavli, F. Hove, and T. I. Tjelta, “Deep penetrating anchor: further development, optimization and capacity verification,” in 10th International Offshore and Polar Engineering Conference, pp. 410–416, June 2000. View at: Google Scholar
 P. Wriggers and J. Korelc, “On enhanced strain methods for small and finite deformations of solids,” Computational Mechanics, vol. 18, no. 6, pp. 413–428, 1996. View at: Google Scholar
 J. C. Simo and F. Armero, “Geometrically nonlinear enhanced strain mixed methods and the method of incompatible modes,” International Journal for Numerical Methods in Engineering, vol. 33, no. 7, pp. 1413–1449, 1992. View at: Google Scholar
 J. C. Simo, F. Armero, and R. L. Taylor, “Improved versions of assumed enhanced strain trilinear elements for 3D finite deformation problems,” Computer Methods in Applied Mechanics and Engineering, vol. 110, no. 34, pp. 359–386, 1993. View at: Google Scholar
 H. Sturm and L. Andresen, “Large deformation analysis of the installation of dynamic anchor,” in 7th European Conference on Numerical Methods in Geotechnical Engineering (NUMGE '10), pp. 255–260, 2010. View at: Google Scholar
 L. Zdravkovic, D. M. Potts, and R. J. Jardine, “A parametric study of the pullout capacity of bucket foundations in soft clay,” Geotechnique, vol. 51, no. 1, pp. 55–67, 2001. View at: Publisher Site  Google Scholar
 C. D. O’Loughlin, M. F. Randolph, and M. Richardson, “Experimental and Theoretical Studies of Deep Penetrating Anchors,” in Offshore Technology Conference, Houston, Tex, USA, 2004. View at: Google Scholar
 S. Narasimha Rao, K. Hema Latha, B. Pallavi, and S. Surendran, “Studies on pullout capacity of anchors in marine clays for mooring systems,” Applied Ocean Research, vol. 28, no. 2, pp. 103–111, 2006. View at: Publisher Site  Google Scholar
 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 Geoenvironmental Engineering, vol. 135, no. 4, pp. 487–496, 2009. View at: Publisher Site  Google Scholar
 U. A. A. Mirza, “Pile shortterm capacity in clays,” in 9th International Offshore and Polar Engineering Conference, pp. 693–699, Brest, France, June 1999. View at: Google Scholar
 P. Jeanjean, “Setup characteristics of suction anchors for soft gulf of Mexico clays: experience from field installation and retrieval,” in 38th Offshore Technology Conference, Houston, Tex, USA, 2006. View at: Google Scholar
 ANSYS, ANSYS Help (version 12) – Mechanical APDL, Ansys Inc, 2009.
 H. Madsen, S. Krenk, and N. C. Lind, Methods of Structural Safety, Dover, EnglewoodCliffs, NJ, USA, 2006.
 R. E. Melchers, Structural Reliability Analysis and Prediction, John Willey and Sons, Chichester, UK, 2nd edition, 1999.
 O. Ditlevsen and H. O. Madsen, Structural Reliability Methods, John Wiley and Sons, New York, NY, USA, 1996.
 A. O. VazquezHernandez, G. B. Ellwanger, and L. V. S. Sagrilo, “Reliabilitybased comparative study for mooring lines design criteria,” Applied Ocean Research, vol. 28, no. 6, pp. 398–406, 2006. View at: Publisher Site  Google Scholar
 A. H. S. Ang and W. H. Tang, Probability Concepts in Engineering Planning and Design, vol. 2, John Willey and Sons, New York, NY, USA, 1984.
 A. O. VazquezHernandez, Reliabilitybased safety factors calibration for the design of mooring lines, D.Sc. thesis, Department of Civil Engineering, COPPE/UFRJ, Brazil, 2004.
 National Cooperative Highway Research Program (NCHRP), “Load and Resistance Factor Design (LRFD) for Deep Foundations,” Report 507, Edited by S. G. Paikowsky, 2004. View at: Google Scholar
Copyright
Copyright © 2012 Luís V. S. Sagrilo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.