 Open Access
 Authors : Iberahin Jusoh
 Paper ID : IJERTV9IS050448
 Volume & Issue : Volume 09, Issue 05 (May 2020)
 Published (First Online): 30052020
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Stress Utilization of Jacket Structure Under Environmental Loading
Iberahin Jusoh
Mechanical Engineering Department, College of Engineering and Islamic Architecture Umm AlQura Univerrsity, Makkah, Saudi Arabia.
AbstractThis paper describes a procedure to identify the critical members in an offshore jacket structure in Malaysia waters by adopting the probabilistic approach. The hydrodynamic design loads on the jacket structure are imposed by considering the joint probability occurrence of the metocean phenomena (significant wave height and spectral peak period). This forms the basic calculation of the longterm response from which the extreme response can be estimated. Stochastic dynamic response analyses are performed in the time domain by stepwise integration. Necessary simulation lengths are established. Interaction ratios of the members are computed based on the API RP2AWSD design code to determine the level of stress utilization ratio of the jackets members. The inclined structural member located at elevation 21.3 m from the sea bed is found to have the highest stress utilization ratio of 0.2884. The results also confirm that the structure is strong by having a high reserve strength.
Keywords Stress utilization ratio; jacket structure; probabilistic analysis; environmental loading.

INTRODUCTION
The probabilistic approach is widely used in the design and analysis of offshore structures. The application of this approach has gained wide acceptance not only in the studies
response statistics under random wave load. The stress analysis is performed in the time domain to identify the critical members.

PROBABILISTIC ENVIRONMENTAL LOADING MODEL
The procedures used for designing offshore structures for extreme conditions, the environmental loads are based on a deterministic, regular wave. This deterministic fluid loading model is highly simplified of an abstract in nature. Studies made by [5] showed that the deterministic fluid loading model predicts fluid loads about 10% lower than those predicted by the probabilistic model. Therefore, the probabilistic model would best describe a random condition in the ocean environment.
The total wave force is calculated based on the extreme conditions, which produces the largest environmental load on the jacket. The sea condition: Hs = 5.2 m, Tp = 9.9s, = 3.3 and the current profile as illustrated in Fig. 1.
MLW
0.90 m/s
of fluid loading but also in the structural analysis due to its ability to capture the randomness of nature. Traditionally, the deterministic method which is the design wave approach is the simplest way of analyzing the fluid loading on a structure and the most wellestablished design practice. In this approach, all the environmental loadings are assumed to act as unidirectional. However, in reality, all the extreme environmental parameters will rarely occur simultaneously and cause the worst structural effect. Extensive researches have been carried out on the effects of the joint probability of environmental loading on the structural response [1], [2], [3],
10 m
0.80 m/s
[4]. The concept of joint probability has received increasing attention as the measured data set of the metocean has available and confidence in the hindcast model has grown. This approach is here employed for analysis fluid loading where the joint probability model for significant wave height and the spectral peak period is considered.In structural analysis, the spectral approach in the
73.3 m 0.40 m/s
Fig. 1: Current profile for storm condition
A Gaussian wave elevation (t), is first simulated as
k 1 k k k
k 1 k k k
(t) n a cos t
(1)
frequency domain is often employed compared to the time domain. It is the preferred method as it is less time consuming compared to the time domain approach. However, the wave
loads on complex ocean structures typically vary nonlinearly
The magnitude of the ak is related to the onesided JONSWAP power spectrum of the waves S(k) by
with the wave elevation. Thus, timedomain simulation remains one of the few general techniques for estimating
ak
2S k
(2)
where k is the circular frequency of a particular harmonic component and k is a random phase.
The applied force f(x,c) per unit length, at elevation z and leg location, is calculated from the Morison equation:
and response stochastic process is given by the dynamic equilibrium equation on the incremental form:
(4)
where M is a mass matrix, C is a damping matrix, K is a
f x, z RDCD ux, z uo ux, z uo R C ux, z
2
M
M
M
M
(3)
stiffness matrix, Ft is a vector of incremental hydrodynamic nodal point forces inline with the wave direction, xt is an
incremental structural displacement vector, x t is an
Absolute fluid velocity is used. The total forces are then found by integrating Morisons equation to the exact free surface. When this is above the mean water level, kinematics from linear wave theory is extended by constant stretching [6]. The Linear Airy Theory is used with empirical modification to take account of the constant stretching.
The current is assumed to be in the same direction as the wave loading. Adjustment to the current profile is made above the mean water level to take account of the wave and current interaction. The current profile is assumed to stretch above the mean wave line until it reaches the water surface as shown in Fig. 2. It is the most plausible since it approximately represents the convection of particles from crest to trough in the wave when that particle has its horizontal velocity increased by the presence of the current [7].
The force exerted by wind is assumed static. Wind loading is fixed at the highest elevation of the jacket structure at the incidence direction of 0degree for the simulation of 0 degree to 45degree. The input of the environmental parameters for the structural analysis is summarised in Table 1.
Mean Water Line
Fig. 2: Current Profile Stretching, the profile is stretched or compressed vertically but the current velocity at any proportion of the instantaneous depth d+ is constant.
TABLE 1: Input Parameters for Environmental Loading
Parameter
Value
Significant Wave Height
5.2 m
Associated Wave Period
9.9 s
Wave Spectrum
Jonswap
CD
0.7
CM
2.0
Wind Speed (Gust)
26.44 m/s
Current Profile
Bilinear Mudline 40cm/s
MWL 90cm/s

NON LINEAR DYNAMIC ANALYSIS
The longterm response sea state is specified by a model spectrum with significant wave height (Hs) and spectral peak period (Tp) as parameters. The relationship between the load
incremental structural velocity vector and x t is an incremental structural acceleration.
Sample time series of wave kinematics are simulated using independent random phases angles uniformly distributed between 0 and 2. The hydrodynamic force vector at each discrete time step is determined from the Morison equation containing the velocity and acceleration of the water particles.
The response process samples are obtained by a stepwise numerical integration procedure. A Newmark mehod with
= 0.25 and = 0.5 is employed, corresponding to the constant average acceleration algorithm. An iteration is required at each step as the force process depends on the response.

STRUCTURAL MODEL
The structure used in this study is a medium water depth jacket structure which is installed in the South China Sea at the water depth of 73.3m. It is a 4legged jacket that has 5 elevations that consist of horizontal, vertical and inclined bracing as shown in Fig. 3.
Fig. 3: Jacket Structure
The structure is made of steel tubular with the Modulus Young of 200 GN/m2 and a yield strength of 250 MPa. Each of the jacket legs is supported with a piledriven through the foundation and sealed with grout. The four legs are battered to achieve better stability against toppling.
The jacket has four elevations. The base, which is position on the seabed, has a dimension of 29.41 m Ã— 29.41 m from the plan view. The first elevation is 21.8 m from the mudline. It has a dimension of 25.05 m Ã— 25.05 m from the top view.
This elevation has the simplest structure as it only consists of horizontal Xbracing. The second and third elevations have the same structure which is the Xbracing and the Tbracing attached to one of its Xbracing. They are respectively 43.3 m and 61.8 m from the mudline with a dimension of 20.75 m Ã—
20.75 m and 17.05 m Ã— 17.05 m from the plan view. The top elevation which is the fourth elevation has the most complex structure as it consists mostly of Tbracing and Kbracing. It has a dimension of 14.05 m Ã— 14.05 m from the plan view at the elevation of 76.8 m from the mudline. The main jacket leg has the largest diameter, which is 1.1684 m and the wall thickness of 0.0127 m.
The jacket structure is modeled by using the finite element simulation software. Simplifications of the structure are made to minimize modeling error and to reduce the simulation time but retaining its essential features. The jacket frame is modeled by using the immersed pipe option. The element has six degrees of freedom, 3 translations in the nodal axis, x, y and z and 3 rotations about its nodal axis, x, y, z. The load of the top deck is evenly distributed on the jacket leg.

SIMULATION STUDY
The results of the probabilistic analysis are presented in the time series form. Unlike the spectral method, where the results are generated in a spectral form with the root means square value, the probabilistic method does not have a particular method in establishing its results. The results generated in the time domain probabilistic method show a degree of variability and randomness. Thus, this variability and randomness directly suggest the simulation lengths [8].
From simulations studies [8], it has been revealed that there is a roughly quadratic relation between maximum wave elevation, max and maximum wave force Fmax. This suggests that Fmax should show a greater scatter than max:
wave elevation generated is 7 sets of simulations with each set of 6 hours each. This makes the total of time which is needed to run a completed of simulation with the variance of 10% and error tolerance of 95% is 42 hours.

RESULTS AND DISCUSSION
The simulations are carried out in fourwave incidence directions; 0degree, 15degree, 30degree and 45degree. Due to the symmetrical features of the structure, the simulations are only performed up to 45 degree. The output generated from the simulations is axial stress for node I and J, bending stress for node I and J, and nodal displacement. These stresses are due to external and transferred forces exerted onto the structure. The horizontal forces cause horizontal shear load on the structure is called base shear and it is influenced by the wave height. In this study, the magnitude of the base shear on the structure is shown in Fig. 4. It clearly illustrated that the base shear increases quadratically with the increase of wave height.
VF max 2V max
(5)
A virtue of this estimate is that because (t) is Gaussian, analytical estimates of Vmax are available [9];
(6)
In the present case study, Tp = 9.9 s and = 3.3. Thus Tz =
7.73 s. Hence Vmax = 0.078. and VFmax = 0.156. These values agree well with the observed value; Vmax = 0.07 and VFmax= 0.14, [5].
From here, simulation lengths can be established. For example, 95% confidence intervals are found by multiplying our best (mean) estimate by 1 2V / n . Thus, if it seeks
(with 95% confidence) that the estimate lies within 10% on
Fig. 4. Variation of base shear to wave height increment for wave direction, 0degree
From the results of the simulation, the stress utilization ratio is calculated for every member for nodes I and J. The stress utilization ratio is a ratio to check the integrity of a member. If the ratio is low, then the reserve strength is high, it means that the structure is strong to withstand the loadings which are imposed on it and vice versa. Generally, if the ratio is less than one the structure is considered safe.
The definition of the stress utilization ratio which is used in this study is based on the criteria set by API Working Stress Design code [10]. Stress utilization ratio for a member under combined stress, combined axial compression/tension, and bending, can be defined as:
either side of the true value, we require
hence;
2V /
n =0.1 and
fa
0.6Fy

fb
Fb
1.0
(8)
required
20V 2
(7)
where fa and fb are output generate by simulations, while Fy
and Fb are allowable axial (tension or bending) stress and
required can apply to estimate the maximum wave elevation, force or response provided V reflects the relevant coefficient of variation (CoV) value. Therefore, the number of simulations required for this study based on the maximum
allowable bending stress.
The stress utilization ratios for critical members are tabulated in Table 2 to Table 5. In Table 2, the maxima value
for stress utilization ratio under the condition of axial and bending stress for 0degree of wave direction is 0.2864. This value decreases as the direction of the wave attack is increased from 0degree to 45degree. The stress utilization ratios are greater in the second and the base elevation. This is due to the force transfer by the members from the top to bottom. The base shear and overturning moment for the wave attack of 0degree are 393.56 kN and 127.89 kNm respectively. Further reporting on critical members stresses leading to structural failure may be found elsewhere [4].
From Tables 2 to 5, the results of the analysis show that the most critical element on the structure is element number
70. This is true for all three wavestructure incidence directions. The maximum stress utilization ratio for this element occurred when the wave is from 0degree direction and the lowest stress utilization ratio is when the wave from the 45degree direction as illustrated in Fig. 5.
TABLE 2: Stress utilization ratio for 5 most critical members in the wave direction of 0 degree
Rank
Node I
Node J
Element
Stress Utilisation
Element
Stress Utilisation
1
70
0.29
70
0.28
2
63
0.22
62
0.22
3
61
0.22
64
0.22
4
57
0.21
57
0.21
5
71
0.20
88
0.20
TABLE 3: Stress utilization ratio for 5 most critical members in the wave directionof 15 degree
Rank
Node I
Node J
Element
Stress Utilisation
Element
Stress Utilisation
1
70
0.29
70
0.28
2
69
0.22
64
0.22
3
63
0.22
62
0.22
4
61
0.22
69
0.20
5
57
0.21
57
0.20
TABLE 4: Stress utilization ratio for 5 most critical members in the wave direction of 30 degree
Rank
Node I
Node J
Element
Stress Utilisation
Element
Stress Utilisation
1
70
0.28
70
0.27
2
69
0.24
69
0.23
3
63
0.22
64
0.22
4
61
0.22
62
0.22
5
57
0.20
59
0.19
TABLE 5: Highest stress utilization ratio of 5 structural members (wave direction of 45 degree)
Rank
Node I
Node J
Element
Stress Utilisation
Element
Stress Utilisation
1
70
0.25
70
0.25
2
69
0.24
69
0.24
3
63
0.19
64
0.22
4
61
0.19
62
0.22
5
59
0.18
59
0.20
Fig. 5: Stress utilization ratio for critical element no. 70 versus wave incidence angle.
Stress Utilisation Ratio
Stress Utilisation Ratio
As a comparison to the results, an inplace analysis is carried out. In the inplace analysis, the only forces experienced by the structure due to its weight, i.e. due to gravity forces, as well as the buoyancy effect from the surrounding fluid are considered. The stress utilization ratio of the most critical member in place analysis is 0.2232. There is a decrease of 28.3% of the stress utilization ratio compared to the structure where the environmental loadings are imposed. From here, the author concludes that the structure is sturdy as the environmental loading only contribution 28.3% to the overall stress experienced by the structure. The comparison for the inplace analysis and the analysis executed with environmental load is illustrated in Fig. 6.
0.3500
0.3000
0.2500
0.2000
0.1500
0.1000
0.0500
0.0000
0.3500
0.3000
0.2500
0.2000
0.1500
0.1000
0.0500
0.0000
0
20
40
60
80
100
Strcutural Member
120
140 160
180
200
0
20
40
60
80
100
Strcutural Member
120
140 160
180
200
Analysis with Environmental Load In Place Analysis
Analysis with Environmental Load In Place Analysis
Fig. 6: Stress utilization ratio for the whole structure – In place analysis for wave direction of 0degree.


CONCLUSSIONS
A procedure for probabilistic analysis in the solution of dynamic timedomain analysis of marine structure has been outlined. The response statistic generated indicates that the simulations generated in the time domain can predict the results of the studies with the coefficient of variation, CoV of 5%.
Results obtained for the jacket platform, indicate that most of the critical members are situated at the second elevation from the mudline. The structural member that has the highest stress utilization ratio is element 70 which is an inclined member
between elevations 21.8 m and 43.3 m. It has a stress utilization ratio of 0.2884. This lowstress utilization ratio meant that the level of reserve strength is high and able to support all loads that may be encountered by the structure. Therefore, it can be concluded that the structure is robust.
REFERENCES

Madsen M.N., Nielsen J.B., Klinting P. and Knudsen J. (1988). A Design Load Method For Offshore Structures Based Upon the Joint Probability of Environmental Parameters. 7th International Conference on Offshore Mechanics and Arctic Engineering. February 712, 1988. USA: ASME, 7580.

BitnerGregersen E.M. and Haver S. (1991). Joint Environmental Model for Reliability Calculations. Proceedings of the 1st Int. Offshore and Polar Engineering Conf. 11th 16th August 1991. UK: ISOPE. 246253.

Kolios, A., Jusoh, I. and Brennan, F.P. (2009), Sensitivity analysis of partial safety factors to the design of offshore steel structures, Proceedings of the 7th International Probabilistic Workshop, The Netherlands, 24 – 25 November 2009

Jusoh, I., Okada, H., and Liu, H., (2000); Application of Probabilistic Method in Assessing the Reliability of Offshore Jacket Structure. Jurnal Teknologi, Universiti Teknologi Malaysia. Dec. 2000.

Hagemeijer P.M. (1989). A Comparison between a Deterministic and Probabilistic Fluid Loading Model for a Jacket Structure. 8th International Conference on Offshore Mechanics and Arctic Engineering. March 1923 1989. The Hague: OMAE. 8997.

Wheeler J.D. (1970). Methods for calculating Forces Produced by Irregular Waves. Journal of Petroleum Technology. 359367.

Gusmestad O.T. and Karunakaran D. (1990). Wave Current Interaction. In: Environmental Forces on Offshore Structures and their Prediction Vol. 26. London: Society for Underwater Technology. 81110.

Torhaug R. (1996). Extreme Response of Nonlinear Ocean Structures: Identification of Minimal Stochastic Wave Input for TimeDomain Simulation. PhD Thesis, Stanford University.

Winterstein S.R. and Torhaug R. (1993). Estimating Extreme Response of Jackup Structures from Limited TimeDomain Simulation. 12th International Conference on Offshore Mechanics and Arctic Engineering. United States of America: ASME, 251258

API RP2A WSD (America Petroleum Institute). (2010). Recommended Practice for Planning, Designing and Constructing Fixed Offshore PlatformsWorking Stress Design. Washington D.C.