Initial Design Phase and Tender Designs of a Jacket Structure Converted into a Retrofitted Offshore Wind Turbine

Jackets are the most common structures in the Adriatic Sea for extracting natural gas. These structural typologies are suitable for relative low water depths and flat sandy sea floors. Most of them have been built in the last 50 years. When the underground source finishes, these structures should be moved to another location or removed if they have reached their design life. Nevertheless, another solution might be considered: change the future working life of these platforms by involving renewable energy and transforming them into offshore wind towers. The present research proposal aims to investigate the possibility of converting actual structures for gas extraction into offshore platforms for wind turbine towers. This simplified analysis is useful for initial design phases and tender design, or generally when available information is limited. The model proposed is a new simplified tool used to study the structural analysis of the jacket structure, developed and summarized in 10 steps, firstly adopted to study the behavior of the oil and gas structure and then for the retrofitted wind tower configuration.


Introduction
Policy support, technology advances, and a maturing supply chain are making offshore wind farms an increasingly viable option for renewable-based electricity generation, harnessing the more consistent and higher wind speeds available offshore instead the onshore solutions.Investment has picked up sharply in recent years, and with fewer restrictions on size and height with respect to their onshore counterparts, offshore wind turbines are becoming giants [1].The growth of offshore wind creates potential synergies with the offshore hydrocarbons sector; integration could bring benefits in terms of reduced costs and improved environmental performance and utilization of infrastructure [2,3].The possibility to electrify offshore oil and gas operations where there are wind farms nearby, or via floating turbines, reduces the need to run diesel or gas-fired generators on the platform and emissions of carbon dioxide (CO 2 ) and air pollutants [4].The scope is to find new uses for existing offshore infrastructure once it reaches the end of its operational life or when the underground source finishes, in ways that might aid energy transitions [5,6].The North Sea, a relatively mature oil and gas basin with a thriving renewable electricity industry, is already seeing some crossover between the sectors: some large oil and gas companies are major players in offshore wind; one former oil and gas company, Orsted in Denmark, has moved entirely to wind and other renewables [7,8].A major problem in offshore wind energy design is the estimation of fatigue life and fatigue loads in an offshore environment [9,10].The study of the best locations of offshore wind turbines is worth investigation as well [11].However, many studies have been devoted to how the dynamic behaviors of wind turbines couple with jacket and pile-supported foundations [12][13][14][15][16][17][18][19][20][21].
This work aims to provide a tool for a general structural check analysis of a different retrofitting activities for the same offshore lattice structure.The research should allow the structural designer to evaluate the environmental loads, the dynamics of the structure, and the effects of changing the structural behavior of the platform over the different cases of study.All these analyses should be available in the early design phase when the data available about the site's environment and seabed conditions are limited.Furthermore, the tools shown should skip the use of the sophisticated and expensive software analysis that requires detailed data and a high level of understanding of the underlying physics.It follows that the calculation procedure will be faster than the long run-times of the software, which is based on time domain simulations of a high number of different load cases.MATLAB and Strand7 have been used as software programs in this research.In terms of accuracy, the most important requirement is that the approach should provide conservative results and at the same time avoid underestimation of loads and deformations.The issue of this project is to find a useful tool that is able to represent the behavior of a generic offshore structure subjected to environmental actions.The structure under investigation can be one of the oil and gas platforms sustained by jacket structures with a long pile foundation.The authors are aware of the complex problem that is presented here.Therefore, the present research is limited, but in agreement with the assumptions considered in this study.Further investigations can be performed in a future work such as tubular joint design [22] and fatigue assessment [23,24].The influence of a different wind turbine such as vertical axis wind turbines [25] and wind turbines in other environments [26,27] is also worth investigating.

Theoretical Background
Offshore structures are complex systems subjected to severe dynamic loads due to wind, waves, currents, and mechanical loads.This makes the design of these structures a complex process that requires expertise from a wide range of fields.Foundations and substructures are specifically the key structural components that differentiate the design of wind towers from those of onshore ones.Due to the long service lifetime of offshore systems, the analysis of the performance of the substructure is a complex task that involves many steps including load analysis, dynamic analysis, evaluation of the fatigue life, as well as long-term deformations.This research presents a simplified methodology for the analysis of the ultimate and dynamics loads with a check on the most stressed pile of the different offshore retrofitted systems.The analysis was carried out such that the data required about the oil and gas platform, wind turbines, and the site are available in the early design phases.Thus, the methodology appears to obtain a conservative estimation.An apparent limitation in the research is the modeling of first the oil and gas platform and second of the offshore wind turbine called 6DOF with all the retrofitting configurations, which was realized in general form, and it found adaptation in different specific applications.The model proposed considers a discrete system of masses and stiffness, which is just a first rough representation of the real offshore lattice structure with pile foundation.Another inherent limitation in the methodology, which is a natural consequence of the approach taken, is the load estimation based on linear analysis and where they are applied.In the model proposed, all the environmental actions are considered at each node of the cantilever beam, while in the real world, all the steel elements are subjected to the wave and wind force.Another issue of the loads is that in this research is that they have not been calculated using partial safety factors regarding the load combination.The present analysis is performed only over one direction with the assumption that the wind and the wave act in the same direction at the same time.Fatigue state analysis is out of the scope of the present research; nevertheless, this is central to the structural analysis, the earthquake analysis, and the checking of the state of the structure like welding, bolts, and grouting, which for sake of simplicity's, have been bypassed.The simplified analysis proposed in this research can be used in the fatigue design and evaluation including the effects of welding, bolts, and the stiffness of tubular connections [28,29].The next stage of the current research study will address fatigue analysis using current (global rules) and advanced approaches (local rules) for the critical structural detail considering the wave and wind loads as indicated in the DNVGLstandards.
The equation of motion associated with a generic structure is a second order linear differential equation with constant coefficients, and it changes form in relation to the structural properties that have to be modeled.The general set of equations representing structural motion can be defined by the following matrix form: that represents a finite set of N ordinary, linear differential equations in N independent coordinates that contain the dominant features of the structural motion.These coordinates are the elements of the column vector: The loading vector p = p(t) is comprised of N loads p i , i = 1, 2, . . ., N where the load p i = p i (t) is located at the respective node i.The vector representation is: In this analysis, the mass matrix M is assumed to be a diagonal matrix with elements m i > 0, i = 1, 2, . . ., N, in which the element's subscript is its associated node point.It has been said that each generalized coordinate represents either the displacement or the rotation of a portion of the structure's mass at a specific node point.The mass lumping method is probably the most popular method of discretizing the supporting framework and the rigid body portions of an offshore structure.In this research, the mass lumping method will be applied where all the offshore mass structure is modeled as lumped at the top of a cantilever beam.The stiffness matrix K for an N degree of freedom structure is a symmetric matrix of N × N elements.The stiffness is defined as the force applied to the structure in order to produce a unitary displacement.The constant k ij in other terms is that force that is required at node i to counteract a unit elastic displacement ξ sj = 1 imposed at node j, under the condition that all displacements ξ si = 0 for i = j.If a displacement condition is applied sequentially to each node, then the net force at each node j can be obtained by superposition.For a structure with N degrees of freedom, the damping matrix C is defined as a symmetric array of N × N constants c ij .Damping is an influence within or upon an oscillatory system that has the effect of reducing, restricting, or preventing its oscillations.In physical systems, damping is the capacity of the system to dissipate the energy stored in the oscillation within itself without damage.In this analysis, the damping force q Di for the structural mode coordinate ξ is assumed to be a linear combination of the generalized coordinate velocities ξi , i = 1, 2, . . ., N: The damping matrix can be cast in several different specialized forms, each of which has the advantage of easily utilizing available experimental data to determine the elements c ij .One such form is Rayleigh damping in which C is proportional to the system's mass and also the system's stiffness: in which a 1 and a 2 are Rayleigh constants, and they are fixed for a given dynamic system.Rayleigh constants are determined using a standardized procedure [30].The latter equation is pre-and post-multiplied by the matrix of the shape functions as X T CX = a 1 X T KX + a 2 X T MX with the orthogonal properties X T KX = diag(ω n ) and X T MX = I.Introducing the nth modal damping factor ζ n , the last result becomes X T CX = diag(a 1 ω 2 n + a 2 ) = diag(2ζ n ω n ).After equating the nth diagonal terms above, it follows that ζ n = (ω n /2)a 1 + (0.5ω n )a 2 .This last result shows that, for an N degree of freedom system for which the N frequencies ω 1 , ω 2 , . . ., ω N are known, the two Rayleigh constants a 1 and a 2 are uniquely determined if any two of ζ n are specified.For instance, if ζ k and ζ m are specified, then the equation yields the following two simultaneous equations from which a 1 and a 2 can be calculated ζ k = (ω k /2)a 1 + (0.5ω k )a 2 and ζ m = (ω m /2)a 1 + (0.5ω m )a 2 .Then, the remaining N − 2 values of ζ n for all n = m can be determined.

Reference Model
Let us suppose that the oil and gas platform studied in this detection is part of a wide context of a gas field, in a water depth of 87 m.The platform sub-structure is a steel jacket, installed by lifting, supporting a three generic level deck for the drilling and extraction activities.The jacket configuration is represented by a four-legged tubular structure, and the jacket base dimensions are approximately 32 m by 20 m (excluding the sleeves).The jacket dimensions at the top interface are 8 m × 8 m.Jacket planes are located at the following elevations from top the bottom with respect to the still water level: +6 m, −8.5 m, −24.5 m, −43.5 m, −64.5 m.The leg diameter varies from 1200 mm by 40 mm at the top to 1600 mm by 50 mm at the bottom.For the bracing elements, the diameter varies from 600 mm by 10 mm-800 mm by 15 mm. Figure 1 is a possible representation of the jacket substructure, the subject of this analysis, realized using Strand7 and MATLAB software.All analysis steps necessary for the simplified design of an offshore structure have been done introducing a 10-step modeling approach: In order to test if the activity model conforms to the real response of the structure, two methods of analysis will be tested for the same structure.One has been developed using MATLAB software and the other one using the Strand7 software.If the results are coherent, the modeling could be said to be efficient and suitable for the representation of the behavior of a generic jacket steel structure.
Platform properties and environmental parameters are described below.The steel jacket is modeled with five lumped masses placed at different heights, one at each horizontal plane and clamped at the bottom.The generic lumped mass is the sum of all the tubular masses surrounding the center of the horizontal plane, in particular half above and half below the generic horizontal plane.A real complex steel jacket structure can be studied using an FE model in Strand7.All the masses have been calculated considering the specific weight of steel 7850 kg/m 3 multiplied by the generic cross-section and multiplied by the length of the tubular members.The diagonal mass matrix M of the 5DOF model expressed in kg developed in MATLAB: The stiffness matrix K has been determined along the direction x.Once having applied the definition of stiffness to all the system, it is possible to define (using MATLAB) the matrix K in N/m as: The sea state analysis has been simplified, and sample values have been chosen through Pierson-Moskowitz spectra.Such a distribution is an empirical equation for deep water waves that defines how the energy is distributed with respect to the wave frequency.It is based on the assumption of the fully-developed sea in which the wind blows constantly for a long time over an unlimited fetch area.Another wave spectrum (termed JONSWAP) might be taken into account.It is typical of the North Sea characterized by a limited fetch area.The aim of this project is not to reproduce a specific study of a particular case, but to find a generic tool to reproduce an easier and faster way expressed in terms of analysis cost of the dynamic behavior of a generic lattice structure.The quality of this model is also promoted by the flexibility, the adaptability of the frame, which represents a valid approach to consider different cases.The reader is free to reproduce the analysis using the JONSWAP spectra.
These data characterize the definition of the wave loads defined by Morison's equation against the structure.The following sea state conditions have been chosen as wave height H = 11.6 m, the wave period T = 12 s, the wavelength λ = 1.56T 2 = 226 m, the wave frequency ω = 2π/T = 0.5236 rad/s, and the salt water density ρ = 1031 kg/m 3 .
The free undamped motion of the structure is investigated solving the eigenvalue problem: Thus, the five undamped natural structural frequencies for free vibration in rad/s of the 5DOF model in rad/s are: Modal shape matrix X is also identified for each natural structural frequency ω n .Each modal vector is identified by ξn = [ ξ1n , ξ2n , . . . ,ξNn ] T = [1, ξ2n , . . ., ξNn ] T .Those vectors are normalized with respect to the mass matrix to form the new modal vectors x n = ξn /e n , where e n is a set of positive and real constants computed from ξT n M ξn = e 2 n .It follows that the modal shape matrix X of the 5DOF model is defined as the assembly of the normalized modal vectors x n : Associated with five undamped natural structural frequencies for free vibration ω n with n = 1, 2, . . ., 5 are the five corresponding vibration modes of the 5DOF model.
A benchmark test with Strand7 software is carried out for verification purposes using natural frequency analysis.Forty vibration modes have been considered, because the FE model is made of several 1D beam elements.The five vibration modes given by the MATLAB code are compared to the ones obtained via FE in Figure 2. The undamped natural frequencies in rad/s are: ω 1 = 14.33, ω 2 = 35.00,ω 3 = 57.93,ω 4 = 76.97, and ω 5 = 104.62.
It can be noted that the mode shapes obtained by MATLAB are close to the ones obtained by FE simulation, as well as the corresponding natural undamped frequencies.
Subsequently, damped natural frequencies should be computed because they better represent the structural behavior of real structures.Rayleigh damping is considered for this purpose.Damped frequencies are characterized by the two Rayleigh constants a 1 and a 2 , which are uniquely determined if any two of the modal damping factors ζ n are specified.Since the first few modes will dominate the motion, it is reasonable to choose Thus, the five damped natural structural frequencies for free vibration in rad/s of the 5DOF model are: ω d1 = 11.4355,ω d2 = 35.8437,ω d3 = 68.7227,ω d4 = 87.5516,ω d5 = 110.0883(11) The fourth step, of the present procedure, regards the response of the structure to the harmonic wave.Airy theory is considered, and it is assumed that the motion of the structure is much smaller than the motion of the wave; thus, Morison's equation can be applied.The inertia counterpart of the flow is represented through the inertia coefficient C M = 2, while the drag counterpart by the drag coefficient C D = 0.8.The structure has been modeled with four vertical legs N l = 4 plus two horizontal cross braces N c = 2, which are normal to the flow.The calculation of the wave forces follows the guidelines reported by classical books of offshore structural modeling [30].Wave loads for the present problem in kN are: The fifth step is related to the pile stress analysis at the foundation level.Global bending moment at the foundation level is considered due to all the forces of the model.The final stress on the pile is determined by the bending moment and the weight of the structure (Figure 3).Moreover, a wind force F wind = 410 kN is considered acting on the superstructure z wind = +7.5 m.Therefore, the force applied on the pile B-1 is F B1 = −8316.28kN.As the numerical benchmark, the same analysis is performed in the FE model.The maximum compression value of the stress in piles given by Strand7 is in the pile B1 and is equal to F B1 = −9922 kN, while the MATLAB calculation has given F B1 = −8316 kN.There is a high correlation between the two values; indeed, they differ by only of 16%.It follows that the 5DOF model proposed in this research is working properly.
The sixth step is based on the time domain solution.These five structural steady state displacements over a time of 30 s, or in other words, the response of the structure subjected to five cycles of wave loading, are depicted in Figure 4.The absolute value of the maximum peaks in m are: ξ 1 = 0.02813, ξ 2 = 0.027214, ξ 3 = 0.02356, ξ 4 = 0.017664, ξ 5 = 0.01073.These values are then compared with the results obtained with the 3σ approach that defines structural safety according to the stochastic approach.The simulation is verified via FE analysis, leading to ξ 1 = 0.0384, ξ 2 = 0.0348, ξ 3 = 0.0278, ξ 4 = 0.0187, ξ 5 = 0.0101.
The maximum displacements obtained with MATLAB and FE are of the same magnitude.The ratios between the displacement obtained with MATLAB and Strand7 are listed in Table 1.Small deviations can be observed.However, this is due to the simplifications introduced; nevertheless, as a preliminary design phase, the results can be considered worthy for proceeding in the present study.It follows that the model fully satisfies the requirements of this research.
The seventh step defines the wave transfer functions for the model.The transfer function G(ω) is a function that relates the wave height H of the incident wave to the load imparted to the structural component.It is generally defined in harmonic form: G(ω) = G 0 e jωt with j = √ −1, and G 0 is a complex number independent of time.The transfer function depends on the flow regimes and the structural component.Transfer functions are computed for all components and assembled for the structure.The transfer function can be applied only to linear analysis (as in the present case), and the effect of multiple wave excitation can be considered (which is not possible in the nonlinear regime).The calculation of the wave transfer function follows the indication provided in [30], and the plot of the modulus square of the five transfer functions |x T k G(p, ω)| 2 = |G( qk , ω)| 2 is given in Figure 5.

Spectral Density Function of the Generalized Force Component in Modal Coordinates
The eighth step is based on the definition of the spectral density function of the generalized force component in modal coordinates.It is necessary to determine first the Pierson-Moskowitz wave spectrum [30]: with B = 0.74(g/V) 4 where v is the wind speed in m/s at a height of 19.5 m above the still water level and ω is expressed in rad/s.Multiplying the transfer function of the load in modal coordinates by the Pierson-Moskowitz spectrum, it achieved the spectral density function of the generalized force component in modal coordinates (Figure 6).
is the modulus of the harmonic response function for the kth mode.Once the spectral density in terms of modal coordinates S(y k , ω) has been determined, is possible to define the spectral density in terms of physical coordinates S(ξ k , ω) = ∑ N n=1 x 2 kn S(y k , ω) applying the inverse of the modal analysis, as shown in [30] (Figure 7).The last step includes the stochastic approach design using the 3σ approach.By computing the area covered by the response displacement spectra, the variance of such displacement can be carried out σ 2 (ξ k ).The extreme limits of ξ k are ±3σ(ξ k ) with k = 1, 2, . . ., 5. If the static stresses and deflection of the members are within these extreme limit values, then the structure is assumed safe.Assuming a Gaussian process, there is only a 0.026% chance that each response exceeds the ±3σ(ξ k ) limits.

Jacket Supporting the Wind Turbine
The aim of the present paper is to find a useful tool able to model the behavior accurately of a new offshore structure that supports the wind turbine subjected to environmental actions.The structure under investigation presents the same jacket structure with a long pile foundation placed in the same position with the same characteristics.All analysis steps necessary for the simplified design for the new offshore structure have been introduced in the previous section.The introduced 10-step modeling approach will be considered also in the following.
The wind tower is made of steel with a length of 80 m.The nacelle mass is 240 tons.The rotor mass is 110 tons.The turbine diameter is 126 m (Figure 8).The wind tower has a truncated cone shape; at the bottom, the diameter is 4 m and the thickness is 0.18 m; at the top, the diameter is 3.5 m and the thickness is 0.10 m.At half tower height, the diameter is 3.75 m and the thickness is 0.14 m.The hub is at z = +90 m.The wind tower produces a rated power of 5 MW.It is attached to the substructure through a concrete (density 1600 kg/m 3 ) transition piece with 9.6 × 9.6 × 4 m 3 of volume.The whole structure will be treated as a discrete system with 5 + 1 lumped masses.
The stiffness matrix K has been determined along the x direction using the proper definition.Thus, matrix K in N/m is: The same sea state of the previous section has been considered.The free undamped motions of the present structure can be determined solving the corresponding eigenvalue problem.The six undamped natural frequencies computed using MATLAB and the FE model are listed in Table 2 in rad/s.The modal shape matrix X of the 6DOF model is defined as the assembly of the normalized modal vectors x n as: and its vibrational modes of shape are compared to the ones obtained by FE analysis in Figure 9.
As follows, all the frequencies computed via FE analysis in rad/s are: ω 1 = 2.74, ω 2 = 7.16, ω 3 = 26.13,ω 4 = 50, ω 5 = 63.81, and ω 6 = 80.39.The modes of shape obtained by MATLAB calculation and the ones obtained by Strand7 are quite similar, as well as the corresponding natural undamped frequencies.The correlation between the two approaches shows the quality of the model.
As done for the previous case, damped frequencies of the 6DOF model in rad/s are computed as:  The harmonic wave of the previous case is considered below.Regarding the wind load, the correspondent force is calculated at z = +90 m.The wind force is given by two contributions: one from the wind acting against the swept rotor area and the other from the wind acting against the tower.The force considered at the top node of the model p 1 = 3055 kN will take into account the wind on the swept rotor area and the wind on the half tower.The remaining half part of the wind on the tower will be add to the first wave load.The wind force is given by two contributions: one from the wind acting against the swept rotor area and the other from the wind acting against the tower.The force considered at the top node of the model p 1 will take into account the wind on the swept rotor area and the wind on the half tower.The remaining half part of the wind on the tower will be added to the first wave load.The counterpart wind on the tower is calculated by integrating the wind velocity over the tower height, and it is equal to f w = u(z)dz = 191.58N with u(z) = 1.5(z/3) 1/7 and with 0 ≤ z ≤ 90 m.The counterpart wind on the rotor is f rot = 0.5ρu rot A = 3054.9kN, where the wind velocity at the rotor hub is u rot = 20 m/s.It follows that the wind load applied at the top of the wind turbine is p 1 = f w /2 + f rot = 3055 kN.
The pile foundation check is carried out as previously described with reference to Figure 10.For this new configuration, the maximum compression value in the piles given in B-1 is equal to FB1 = −21,131 kN.The same pile in the 5DOF configuration withstood a compression equal to FB1 = −8316 kN.In the new configuration, the pile B-1 carries more than two and half-times the stress carried in the previous one.Supposing that the compression threshold pile limit is Fmax = −18,000 kN, it is quite clear that the new load produced by the 6DOF configuration is no longer bearable for the pile B-1.
The generic long steel pile foundation in this structure presents the following characteristics: diameter of 2.13 m, wall thickness of 50 mm, pile length of 96 m, of which 75 m are embedded.The corresponding axial pile capacity is Q r = 28,800 kN.In order to take into account the differences between the structure with the cantilever model, a safety factor equal to γ r = 1.60 has been chosen.It follows that the maximum bearing load for the pile is F max = Q r /γ r = 18,000 kN.Performing the static analysis with Strand7 of the 6DOF model, the maximum displacement for each node is: ξ 1 = 1.3893, ξ 2 = 0.1012, ξ 3 = 0.712, ξ 4 = 0.434, ξ 5 = 0.0244, ξ 6 = 0.0103.These values will be compared with the 3σ(ξ k ) approach below.Although the jacket substructure has the same configuration as the 5DOF model, the wave transfer functions are not identical, but similar to the previous case.This is due to the mode shape matrix X of the 6DOF model.It follows that the equations are the same as the 5DOF configuration, but the functions of the five wave transfer functions are slightly different.
In order to define the spectral density function of the generalized force component, it is requested to define first the wave and wind spectra.For the wind analysis, the Kaimal spectrum is used.The parameters chosen for the Kaimal spectrum are: where the mean wind velocity U mean = 6 m/s, the variance of the wind speed 81 (m/s) 2 with λ w = 6.7 and k w = 1.72 Weibull parameters of probability density of gamma function Γ, roughness parameter z 0 = 0.001, the integral length scale L u = 300( el 1 300 ) 0.46+0.074log(z 0 )= 321 m, and the air density ρ air = 1.225 kg/m 3 (Figure 11).Wave modeling is carried out using the Pierson-Moskowitz spectrum.Multiplying the wind transfer function of the load H u f = C T Aρ air U mean in the modal coordinates with the thrust coefficient C T = 0.8 by the Kaimal spectrum S u for the wind and the wave transfer functions |x T k G(p, ω)| 2 by the Pierson-Moskowitz S η , the spectral density function of the wind force in modal coordinates S f = H 2 u f • S u and the spectral density of the generalized wave force in modal coordinates S qi , respectively, are achieved (Figure 12).The spectral density of the wind load in modal coordinates i = 1 has values much higher than the wave loads; in particular, it tends to infinity when ω tends to zero.For the sake of simplicity, the spectral density function of the wind force in modal coordinates will be limited.The wind transfer function of the load H u f in modal coordinates, which converts the wind speed spectrum to the wind load spectrum, is based on the assumption that u(t) has the same value over the area, which in general is not true.
Once the spectral density function of the wind load and generalized wave component are defined in modal coordinates S f and S qi , the wind transfer function from the spectral density of wind force into spectral density in modal coordinates is carried out as: and the wave transfer functions from spectral density of the generalized wave force component into spectral density in modal coordinates |H k (ω)| 2 .The spectral density function in terms of modal coordinates S(y k , ω) is given by multiplying the spectral density function of the load in modal coordinates with the correspondent transfer function in modal coordinates.Let J 1 be the joint acceptance of the mode, which accounts for the correlation of wind loading along the rotor, considering the mode shape.
Once the spectral density in terms of modal coordinates S(y , ω) has been determined, it is possible to define the spectral density in terms of physical coordinates S(ξ k , ω) (Figure 13).The last 10-step methodology phase defines the variance of the structural displacements in physical coordinates σ 2 (ξ k ).The extreme limits of ξ k are ±3σ(ξ k ) with k = 1, 2, . . ., 6.If the static stresses and deflection of the members are within these extreme limit values, then the structure is assumed safe.Considering a Gaussian process, there is only a 0.026% chance that each response exceeds the ±3σ(ξ k ) limits.The extreme displacement values of the 6DOF configuration in m are: 3σ(ξ 1 ) = 1.3822, 3σ(ξ 2 ) = 0.3471, 3σ(ξ 3 ) = 0.2917, 3σ(ξ 4 ) = 0.2154, 3σ(ξ 5 ) = 0.1505, 3σ(ξ 6 ) = 0.0748.After comparing these results with the ones given above, it can be stated that the structure cannot be safe for the first node at the top of the wind tower, in particular ξ 1 = 1.40 m, whereas the limit is 3σ(ξ 1 ) = 1.38 m.The 6DOF model built from the adaptation and extension of the 5DOF model has to be modified in order to find the safety condition.In order to do that, hereafter, different retrofitting examples will be presented.

Improved Retrofitting Models
Retrofitting refers to the addition of new technology or features to older systems, extending the service life.In this case, different retrofitting typologies are illustrated.The retrofitting models proposed in this section are analyzed following the same methodology previously shown.The approach is suitable for tender design and early design steps of jacket structures.The aim is to identify a useful tool that is able to represent the behavior of the new offshore structure subjected to environmental actions.The system under investigation presents the same jacket and a big wind tower installed on top of it.The whole system will be analyzed with MATLAB calculation and supported by Strand7 software.The aim is to find the right intervention that will ensure the modified offshore structure without causing stability trouble for the system.This research proposes three main retrofitting activities that can be adopted for the offshore wind turbine.They are: crown pile, mooring lines, and 2-MW wind turbine.

Crown Pile Configuration
The crown pile solution represents the same configuration of the 6DOF model except for the piles at the base of the jacket.In addition to the four corner piles, here there are an extra four piles placed at the midspan of each side at z = −d connected each to the other piles with a rectangular ring (Figure 14).Regarding Figure 14 and all the following configurations, z = −d represents the sea bottom level when the vertical axis z is zero at the Still Water Line (SWL) pointing upward.The reference system shown in Figure 14 of the FE model considers the reference system at the sea bottom, so z = 0.However, all the present MATLAB models consider the z axis pointing upward with the sea bottom at z = −d.frequency ω n exists a modal vector ξ n , which can be computed solving the linear eigenvalue problem (see [30] p. 226).
The crown pile mode shapes are represented in Figure 15.
Wave and wind loads do not change with respect to the previous 6DOF case.
The pile foundation check leads to a maximum compression in the pile B-1 equal to F B1 = −17,178 kN lower than the pile compression limit F max = −18,000 kN.The reader can easily state that this first configuration helps the wind turbine in terms of stability.The maximum node displacements of the crown pile model in m are: ξ 1 = 1.3893, ξ 2 = 0.1011, ξ 3 = 0.711, ξ 4 = 0.433, ξ 5 = 0.0243, ξ 6 = 0.0101.These values will be compared with the results obtained with the 3σ(ξ k ) approach to ensure structural safety.
The five wave transfer function are rather similar to the ones obtained for the 6DOF.Multiplying the wind transfer function of the load H u f in modal coordinates by the Kaimal spectrum S u for the wind and the wave transfer functions |x T k G(p, ω)| 2 by the Pierson-Moskowitz S η , the spectral density function of the wind force in modal coordinates S f and the spectral density of the generalized wave force in modal coordinates S qi , respectively, are achieved.Such spectral densities are depicted in Figure 16.The response spectra (Figure 17) for the displacements in physical coordinates of the offshore wind turbine crown pile configuration have been determined.The area beneath these spectral density functions represents the variance of the generic displacement in physical coordinates σ 2 (ξ k ).The extreme limits of ξ k are ±3σ(ξ k ) with k = 1, 2, . . ., 6.If the static stresses and deflection of the members are within these extreme limit values, then the structure is assumed safe.Considering a Gaussian process, there is only a 0.026% chance that each response exceeds the ±3σ(ξ k ) limits.
Comparing these results with the ones given by Figure 15, the reader will easily state that the structure cannot be safe for the first node at the top of the wind tower, in particular ξ 1 = 1.3893 m and 3σ(ξ 1 ) = 1.3821.The crown pile configuration has obtained higher values of safety expressed by the 3σ(ξ) approach results, but this is not enough yet for the global stability.

Mooring Line Configuration
The mooring line configuration presents the same configuration of the 6DOF model, which in addition has a mooring line system that links the midspan second horizontal plane point of the jacket to the seabed, limiting the movement (Figure 18).The mooring line has been modeled as a steel truss element subjected to tension action given by the environmental loads.It has a diameter cross-section of 0.1 m and a length of 98 m, and it creates an angle of 60 • with the seabed.For this proposal, only the x direction will be considered, but for real cases, it could be extended and applied for each direction.
The stiffness matrix K has been determined over the x direction.Once having applied the definition of stiffness to all the system, it is possible to define the matrix K in N/m as: The sea state is unchanged with respect to the previous cases.
The response of the mooring line configuration to the harmonic waves has been considered using the same hypothesis for the 6DOF model.The wave loads and the wind load present the same values for the 6DOF system.A foundation check is performed for both the piles and the mooring line anchoring.The tension in the mooring line from Strand7 calculation is F ml = 427.5 kN.The calculation of the stress applied on the pile B-1 considers all the environmental loads p i with the horizontal and vertical component of F ml that are applied to the structure.The compression in the pile B-1 is equal to F B1 = −20,959 kN, which is still too high for the compression pile limit capacity.The response spectra (Figure 21) for the displacements in physical coordinates of the offshore wind turbine crown pile configuration have been determined.The area beneath these spectral density functions represents the variance of the generic displacement in physical coordinates σ 2 (ξ k ).The extreme limits of ξ k are ±3σ(ξ k ) with k = 1, 2, . . ., 6.If the static stresses and deflection of the members are within these extreme limit values, then the structure is assumed safe.Considering a Gaussian process, there is only a 0.026% chance that each response exceeds the ±3σ(ξ k ) limits.The extreme displacement values of the crown pile configuration in m are: 3σ(ξ 1 ) = 1.3965, 3σ(ξ 2 ) = 0.3768, 3σ(ξ 3 ) = 0.3204, 3σ(ξ 4 ) = 0.2459, 3σ(ξ 5 ) = 0.1724, 3σ(ξ 6 ) = 0.0995.
Comparing these results with the ones previously given, the reader will easily state that the structure cannot be safe for the first node at the top of the wind tower, in particular ξ 1 = 1.3854 m and 3σ(ξ 1 ) = 1.3965.The mooring configuration has obtained higher values of safety expressed by the 3σ(ξ) approach results, but it is not enough yet for the global stability.

The 2-MW Wind Turbine Configuration
The modeling of the 2-MW configuration is characterized by the same jacket sub-structure, but with a smaller wind tower.Tower length is 60 m.The nacelle mass is 57 tons.The rotor mass is 23 tons.The diameter is 66 m.The wind tower has a truncated cone shape.At the bottom, the diameter is 4 m and the thickness is 0.18 m.At the top, the diameter is 3.5 m and the thickness is 0.10 m.At half the tower height, the diameter is 3.75 m and the thickness is 0.14 m.The hub is at z = +70 m.The wind tower produces a rated power of 2 MW.As for the 6DOF model, it is attached to the sub-structure through the same concrete transition piece.The sketch of the 2-MW offshore wind turbine configuration is given in Figure 22.
The stiffness matrix K has been determined over the x direction.Once having applied the definition of stiffness to all the system, it is possible to define the matrix K in N/m as: The sea state remains unchanged with respect to the other configurations.
The six undamped natural structural frequencies for free vibration in rad/s are: ω 1 = 4.5039, ω 2 = 8.8834, ω 3 = 36.1127,ω 4 = 71.9487,ω 5 = 90.6417,ω 6 = 111.5328.They are also graphically represented in Figure 23.The response of the harmonic wave has been considered the same as the 6DOF model.The wind load due to the different wind towers will be lower in particular if it is applied at z = +70 m and is equal to P w = 838 kN.
The compression in the pile B1 is equal to F B1 = −13,288 kN, lower than the compression limit value of the pile.
The five wave transfer functions are similar to the ones obtained for the previous cases.
Multiplying the wind transfer function of the load H u f in modal coordinates by the Kaimal spectrum S u for the wind and the wave transfer functions |x T k G(p, ω)| 2 by the Pierson-Moskowitz S η , the spectral density function of the wind force in modal coordinates S f and the spectral density of the generalized wave force in modal coordinates S qi , respectively, are achieved (Figure 24).For the spectral density of the wind load in modal coordinates, the same consideration is made for the 6DOF case.
As done for the previous model, the response spectra for the displacements in physical coordinates of the offshore wind turbine crown pile configuration have been determined and are shown in Figure 25.The area beneath these spectral density functions represents the variance of the generic displacement in physical coordinates σ 2 (ξ k ).The extreme limits of ξ k are ±3σ(ξ k ) with k = 1, 2, . . ., 6.If the static stresses and deflection of the members are within these extreme limit values, then the structure is assumed safe.Considering a Gaussian process, there is only a 0.026% chance that each response exceeds the ±3σ(ξ k ) limits.
Comparing these results with the ones given above, the reader will easily state that the structure cannot be safe for the first node at the top of the wind tower, in particular ξ 1 = 0.2989m and 3σ(ξ 1 ) = 0.4208.The 2-MW configuration has obtained lower values of safety expressed by the 3σ(ξ) approach results and together with the pile check analysis completes the global stability of the 2-MW retrofitted configuration.

Conclusions
Due to the long service lifetime of offshore systems, the analysis of the performance of the substructure is a complex task that involves many steps including load analysis, dynamic analysis, evaluation of the fatigue life, as well as long-term deformations.
This work presents a simplified methodology for the analysis of the ultimate and dynamic loads with a check on the most stressed pile of the different offshore retrofitted systems.The analysis was carried out such that the data required about the oil and gas platform, wind turbines, and the site are available in the early design phases.Thus, the methodology appears to obtain a conservative estimation.It is obvious that the present study is not free from limitations, but it can be improved in further investigations.
The first limitation is due to the simplified modeling made of discrete masses, one for each jacket level.Nevertheless, high correlations have been shown with the example proposed by the references.
The present work proposes five retrofitted solutions: four related to the substructure operations and one regarding the wind turbine.They are respectively the crown pile configuration, long pile, mooring line, stirrups, and the 2-MW wind turbine configuration.The safety concept of the retrofitted models is expressed in terms of pile check analysis and the 3σ approach.
Nevertheless, the final values obtained by the different configuration typologies are slightly higher or lower than the threshold limit, and the check is slightly fulfilled or not; here, the goal is just to show a simple method to investigate an offshore structure.The present research can be considered as the starting point for a reliability analysis regarding decommissioned jacket platforms for wind energy generations in an offshore environment.

1 .
Structural and environmental parameters 2. Undamped motion of the structure 3. Damped frequencies' definition 4. Response to environmental loads 5. Maximum load pile applied 6.Time domain solution 7. Transfer function definition 8. Spectral density function of the generic load 9. Response spectrum 10. 3σ approach

Figure 2 .
Figure 2. Mode shape comparison between 5DOF and the FE model.

Figure 3 .
Figure 3. Pile stress analysis of the 5DOF model.

Figure 4 .
Figure 4. Time domain plot of the 5DOF model.

Figure 5 .
Figure 5. Wave transfer functions for the 5DOF model.

Figure 6 .
Figure 6.Spectral density functions of the generalized force component in modal coordinates for the 5DOF model.The ninth step is focused on the spectral density in terms of modal coordinates, and it is given by S(y k , ω) = |H k (ω)| 2 S( qk , ω) with k = 1, 2, . . ., 5. S( qk , ω) is the spectral density function of the generalized force component in modal coordinates, and |H k(ω)| = [(ω 2 k − ω 2 ) 2 + (2ζ K ω k ω) 2 ] −1/2is the modulus of the harmonic response function for the kth mode.Once the spectral density in terms of modal coordinates S(y k , ω) has been determined, is possible to define the spectral density in terms of physical coordinates S(ξ k , ω) = ∑ N n=1 x 2 kn S(y k , ω) applying the inverse of the modal analysis, as shown in[30] (Figure7).

Figure 7 .
Figure 7. Response spectra for the displacement in physical coordinates of the fixed-leg platform.

Figure 8 .
Figure 8. Wind tower model.The lumped masses of the jacket are calculated using the same procedure illustrated above.It follows that the mass matrix M of the 6DOF model in kg is:

Figure 9 .
Figure 9. Vibrational mode comparison between the 6DOF and FE models.

Figure 10 .
Figure 10.Pile stress analysis of the 6DOF model.

Figure 11 .
Figure 11.Wave and wind spectral density function.

Figure 12 .
Figure 12.Spectral density of the generic environmental load in modal coordinates of the 6DOF model.

Figure 13 .
Figure 13.Response spectra for the displacement in physical coordinates of the offshore wind turbine plot.

Figure 16 .
Figure 16.Crown pile spectral density of the generic environmental load in modal coordinates.

Figure 17 .
Figure 17.Response spectra for the displacements in physical coordinates of the offshore wind turbine crown pile configuration.

Figure 18 .
Figure 18.Mooring line configuration model detail.The mooring line masses are calculated following the same approach shown previously.Here, the mooring pile configuration mass matrix M in kg:
The damped frequencies of the 6DOF model in rad/s are: ω d1 = 3.3695, ω d2 = 8.2803, ω d3 = 33.0309,ω d4 = 64.0793,ω d5 = 80.2370, ω d6 = 95.8081.The maximum nodal displacements of the mooring model in m are: ξ 1 = 1.3854, ξ 2 = 0.0983, ξ 3 = 0.0685, ξ 4 = 0.0416, ξ 5 = 0.0235, ξ 6 = 0.0099.These values will be then compared with the results obtained with the 3σ(ξ k ) for checking structural safety.The five wave transfer functions are similar to the ones obtained for the previous case.Multiplying the wind transfer function of the load H u f in modal coordinates by the Kaimal spectrum S u for the wind and the wave transfer functions |x T k G(p, ω)| 2 by the Pierson-Moskowitz S η , the spectral density function of the wind force in modal coordinates S f and the spectral density of the generalized wave force in modal coordinates S qi , respectively, are achieved.These spectra are illustrated byFigure 20.

Figure 20 .
Figure 20.Mooring line spectral density of the generic environmental load in modal coordinates.

Figure 21 .
Figure 21.Response spectra for the displacements in physical coordinates of the offshore wind turbine mooring line configuration.

Figure 22 .
Figure 22.The 2-MW model detail.All the masses have been calculated considering the specific weight of the steel 7850 kg/m 3 multiplied by the generic cross-section and multiplied by the length of the tubular member.It follows that the mass matrix M of the 6DOF model in kg is:

Figure 24 .
Figure 24.The 2-MW configuration spectral density of the generic environmental load in modal coordinates.

Figure 25 .
Figure 25.Response spectra for the displacements in physical coordinates of the 2-MW offshore wind turbine configuration.