WES Wind Energy Science WESWind Energ. Sci. 2366-7451 Copernicus Publications Göttingen, Germany 10.5194/wes-4-23-2019A comparison study on jacket substructures for offshore wind turbines based on optimizationA comparison study on jacket substructures HäfeleJan j.haefele@isd.uni-hannover.de https://orcid.org/0000-0002-8896-4589 GebhardtCristian G. RolfesRaimund https://orcid.org/0000-0001-7714-3382 Leibniz Universität Hannover/ForWind, Institute of Structural Analysis, Appelstr. 9a, 30167 Hanover, Germany Jan Häfele (j.haefele@isd.uni-hannover.de)22January2019 4 1 2340 28August2018 17September2018 26December2018 6January2019 This work is licensed under the Creative Commons Attribution 4.0 International License. To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/This article is available from https://wes.copernicus.org/articles/4/23/2019/wes-4-23-2019.htmlThe full text article is available as a PDF file from https://wes.copernicus.org/articles/4/23/2019/wes-4-23-2019.pdf

The structural optimization problem of jacket substructures for offshore wind turbines is commonly regarded as a pure tube dimensioning problem, minimizing the entire mass of the structure. However, this approach goes along with the assumption that the given topology is fixed in any case. The present work contributes to the improvement of the state of the art by utilizing more detailed models for geometry, costs, and structural design code checks. They are assembled in an optimization scheme, in order to consider the jacket optimization problem from a different point of view that is closer to practical applications. The conventional mass objective function is replaced by a sum of various terms related to the cost of the structure. To address the issue of high demand of numerical capacity, a machine learning approach based on Gaussian process regression is applied to reduce numerical expenses and enhance the number of considered design load cases. The proposed approach is meant to provide decision guidance in the first phase of wind farm planning. A numerical example for a National Renewable Energy Laboratory (NREL) 5 MW turbine under FINO3 environmental conditions is computed by two effective optimization methods (sequential quadratic programming and an interior-point method), allowing for the estimation of characteristic design variables of a jacket substructure. In order to resolve the mixed-integer problem formulation, multiple subproblems with fixed-integer design variables are solved. The results show that three-legged jackets may be preferable to four-legged ones under the boundaries of this study. In addition, it is shown that mass-dependent cost functions can be easily improved by just considering the number of jacket legs to yield more reliable results.

 Introduction

The substructure contributes significantly to the total capital expenses of offshore wind turbines and thus to the levelized costs of offshore wind energy, which are still high compared to the onshore counterpart . Cost breakdowns show ratios of about 20 % such as depending on rated power, water depth, and what is regarded as capital expenses. In the face of wind farms with often more than 100 turbines, it is easily conceivable that a slight cost reduction can already render substantial economic advantages to prospective projects. Structural optimization is paramount because it provides the great opportunity to tap cost-saving potential with low economic effort. Technologically, it is expected that the jacket will supersede the mono-pile when reaching the imminent turbine generation or wind farm locations with intermediate water depths from about 40 to 60 m see, for instance,. According to current studies, there is an increasing market share of jackets . As it allows for many variants of structural design, the jacket structure is therefore a meaningful object of structural optimization approaches, which benefits massively from innovative design methods and tools .

It is state of the art in the field of jacket optimization to deal with optimal design in terms of a tube dimensioning problem, where the topology is fixed.

This work focuses on the problem of jacket optimization and disregards other substructure types. For a comprehensive overview of the structural optimization of wind turbine support structures, .

Structural design codes require the computation of time domain simulations to perform structural code checks for fatigue and ultimate limit state. As environmental conditions in offshore wind farm locations vary strongly, commonly thousands of simulations are necessary to cover the effect of varying wind and wave states for verification.

During conceptual design phases, the number of load cases is commonly reduced.

 Therefore, numerical limitations are a great issue in state-of-the-art jacket optimization approaches. In the literature, different approaches were presented to address this issue. proposed an optimization scheme based on a meta-heuristic genetic algorithm to guarantee global convergence. To increase the numerical efficiency, a reanalysis technique was applied. Later, an improved approach was illustrated , where the load calculation was decoupled from the actual tube dimensioning procedure and a simplified fatigue load set was applied. Similar approaches by and applied sequential quadratic or linear programming methods, respectively, with analytically derived gradients. Other optimization approaches using meta-heuristic algorithms were reported by and but without comprehensive load assumptions. The problem of discrete design variables was addressed by . presented a jacket optimization study, where different simulation codes were deployed to perform structural code checks. All mentioned works, except for the last one, represent tube sizing algorithms applied to the Offshore Code Comparison Collaboration Continuation (OC4) jacket substructure for the National Renewable Energy Laboratory (NREL) 5 MW reference turbine ,

It is worth mentioning that the Offshore Code Comparison Collaboration Continuation (OC4) jacket is actually a structurally reduced derivation of the so-called UpWind jacket , which was created to ease calculations within the verification efforts in the OC4 project. Therefore, it is not guaranteed that the OC4 jacket is an appropriate comparison object, as it does not incorporate details of tubular joints.

 where the initial structural topology is maintained even in the case of a strong tube diameter and wall thickness variations. Furthermore, it can be stated that all proposals share the entire mass of the jacket as an objective function to be minimized, which is meaningful in terms of tube sizing. Due to numerical limitations, the utilized load sets are altogether small, for instance with low numbers of production load cases or the omission of special extreme load events. These assumptions constitute drawbacks when considering jacket optimization as part of a decision process in early design stages, where basic properties like the numbers of legs or bays are more critical than the exact dimensions of each single tube. Therefore, an optimization scheme which addresses the early design phase is highly desirable to provide decision guidance for experienced designers. Proposals tackling this kind of problem were given by and , where technically oriented jacket models were proposed but lacking fatigue limit state checks in the first and detailed load assumptions in the second case. Based on the latter and with improved load assumptions, a hybrid jacket for offshore wind turbines with high rated power was designed . Due to innovative materials (the technology readiness level of such a structure is still low), this work lacked detailed cost assumptions. Another proposal for an integrated design approach was made by , considering varying bottom widths and soil properties. This work is meant as an approach for conceptual design phases. However, our conclusion on the state of the art is that an optimization approach without massive limitations is still missing.

This work is intended as a contribution to the improvement of the state of the art by considering jacket optimization in a different way. Compared to other works in this field, the focus is on

the incorporation of topological design variables in the optimization problem, while the dimensioning of tubes is characterized by global design variables;

more detailed cost assumptions;

more comprehensive load sets for fatigue and ultimate limit state structural design code checks;

a change in the exploitation of jacket optimization results. This work intends to consider jacket optimization as a part of the preliminary design phase because it is assumed that the (economically) most expensive mistakes in jacket design are made at this stage of the design process.

 A basis to address these points was given by , where appropriate geometry, cost, and structural code check models for fatigue and ultimate limit states were developed. In this study, these models are deployed within an optimization scheme to obtain optimal design solutions for jacket substructures. A more efficient or accurate method to solve the optimization problem is deliberately not provided in this study. The authors believe that there are numerous techniques presented in the literature that are able to solve the jacket optimization problem.

The paper is structured as follows. Sect.  describes the technical and mathematical problem statements. Both the objective and the constraints are presented and explained in Sect. . The optimization approach and methods to solve the problem are discussed in Sect. . Section  illustrates the application of the approach to a test problem, a comparison of jackets with different topologies, performed for an NREL 5 MW turbine under FINO3 environmental conditions. This section comprises a detailed setup of the problem and a discussion of the results. The work ends with a consideration of benefits and limitations (Sect. ) and conclusions (Sect. ).

 Problem statement

This paper presents a study on jacket substructures, based on optimization. The design of jackets is a complex task that requires profound expertise and experience. Therefore, it has to be clarified that this work does not provide a method replacing established design procedures. It is rather meant as guidance in early design phases, where it is desirable to define the basic topology and dimensions of the substructure. In industrial applications, this step is commonly highly dependent on the knowledge of experienced designers. Along with this statement, it has to be pointed out that the term “optimal solution” may indicate a solution that it is indeed optimal concerning the present problem formulation but not necessarily optimal in terms of a final design due to the following aspects.

Although the approach deploys more detailed assumptions on the modeling of costs and environmental conditions, compared to optimization approaches known from the literature, it still incorporates simplifications, mainly for the sake of numerical efficiency.

No sizing of each single tube is performed, for the same reason. This is a matter of subsequent design phases, and tube dimensioning approaches exist in the literature. Instead, tube dimensions are derived by global design variables.

The design of pile foundation and transition piece is not performed in this approach. The reason is that both are considered in models of the structure and the costs but are not impacted by the selected design variables.

Only fatigue and ultimate limit state are assumed to be design-driving constraints. Serviceability limit state, i.e., eigenfrequency constraints, is not regarded as design-driving in this work because the modal behavior of a wind turbine with jacket substructure is strongly dominated by the relatively soft tubular tower. In addition, a design leading to eigenfrequencies close to 1P or 3P excitation would probably fail due to high fatigue damage. Although the modal behavior is also impacted by the foundation, this is not significant here, as no foundation design is performed.

The overall goal of jacket optimization can be interpreted as a cost minimization problem involving certain design constraints. As stated before, it is assumed that the design-driving constraints of jackets are fatigue and extreme loads. In other words, a set of design variables for a parameterizable structure that minimizes its costs, Ctotal, is desirable, while fatigue and ultimate limit state constraints are satisfied; i.e., the maximal normalized tubular joint fatigue damage (among all tubular joints), hFLS, is less than or equal to 1,

All fatigue damage is normalized so that the lifetime fatigue damage corresponds to a value of 1.

and the extreme load utilization ratio (among all tubes), hULS, is less than or equal to 1.

The total expenses are defined as an objective function f(x), which depends on an array of design variables, x: f(x)=log⁡10Ctotal(x). In this equation, the cost value is logarithmized to obviate numerical issues. The constraints, h1(x) and h2(x), are formulated so as to match the requirements of mathematical problem statements; thus h1(x)=hFLS(x)-1,h2(x)=hULS(x)-1, depending also on the array of design variables, x.

Based on the technical problem statement, we define the mathematical problem statement in terms of a nonlinear program: min⁡f(x)such that xlb≤x≤xub,h1(x)≤0 and h2(x)≤0, where x is the array or vector of design variables, xlb and xub are the lower and upper boundaries, respectively, f(x) is the objective function, covering the costs related only to the substructure, and h1(x) and h2(x) are nonlinear constraints representing structural code checks for fatigue and ultimate limit state that are required to be satisfied for every design.

 Objective and constraints

This section illustrates the jacket model, which is the basis for the optimization study. Moreover, the models for costs and structural design code checks are described, which depict the objective and constraint functions, respectively. These models were elaborated on in a previous work .

Jacket modeling and design variables

In this work, it is assumed that a jacket substructure can be described by 20 parameters in total, of which 10 define topology, 7 tube dimensions, and 3 material properties. Topological parameters are the number of legs, NL, number of bays, NX (both integer variables), foot radius, Rfoot, head-to-foot radius ratio, ξ, jacket length, L, elevation of the transition piece over mean sea level, LMSL, lowermost segment height, LOSG, uppermost segment height, LTP, the ratio of two consecutive bay heights, q, and a boolean flag, xMB, determining whether the jacket has mud braces (horizontal tubes below the lowermost layer of K joints) or not. The topology of one example with four legs (NL=4), four bays (NX=4), and mud braces xMB=true is shown in Fig. . The tube sizing parameters are the leg diameter, DL, and six dependent parameters defining relations between tube diameters and wall thicknesses at the bottom and top of the structure: γb and γt are the leg radius to thickness ratios, βb and βt are the brace-to-leg diameter ratios, and τb and τt are the brace-to-leg thickness ratios, where the indices b and t indicate values at the bottom and the top of the jacket, respectively. Using dependent parameters is beneficial because structural code checks are valid for certain ranges of these dependent variables. Furthermore, for structural analysis, the material is assumed to be isotropic and can thus be described by a Young's modulus, E, a shear modulus, G, and density, ρ.

Jacket geometry model with variables characterizing the topology of the structure, shown exemplarily for a jacket with four legs, four bays, and mud braces. The ground layer is illustrated by the orange surface and the mean sea level and transition piece layers by the blue and gray surfaces, respectively.

To decrease the dimension of the problem, height measures related to the location of the wind farm (L, LMSL, LOSG, LTP) and the material parameters (E, G, ρ) are fixed. In addition, it is supposed that each design has mud braces (xMB=true). Although designs without mud braces are also imaginable, fixing this parameter is advantageous, as it is not continuous. The array of design variables therefore has a dimension of 12: x=NLNXRfootξqDLγbγtβbβtτbτtT. The number of design variables is not necessarily minimal, but, on the one hand, mathematically manageable and, on the other hand, meaningful from the technical point of view.

 Cost function (objective)

The total capital expenses, Ctotal, comprise several terms, Cj, expressed as the sum of so-called factors, cj, weighted by unit costs, aj:

Unit cost values are given in Sect. .

Ctotal(x)=∑Cj(x)=∑ajcj(x). A factor may be any property of the structure describing a cost contribution that can be expressed in terms of the design variables. A pure mass-dependent cost modeling approach, as used in most optimization approaches, would involve only one factor, while no unit cost value is required for weighting. However, a realistic cost assessment involves more than only the structural mass. For example, in the case of a structure with very lightweight tubes but many bays, it can be imagined that the manufacturing costs tend to be a cost-driving factor. To consider known, important impacts on jacket capital expenses, seven factors are incorporated, namely the following:

expenses for material, C1, depending on the mass, c1: c1(x)=2ρNLπDL2∑i=1NXβiτi2γi+τi24γi2Li2cos⁡2Φp+Ri+Ri+12sin⁡2ϑ2+xMBρNLπDL2βbτb2γb+τb24γb22R1sin⁡ϑ2+ρNLπDL2∑i=1NX12γi+14γi2Lm,icos⁡Φs+12γi+1+14γi+12Li-Lm,icos⁡Φs+ρNLπDL212γb+14γb2LOSGcos⁡Φs+ρNLπDL212γt+14γt2LTPcos⁡Φs;

expenses for fabrication, C2, depending on the entire volume of welds, c2: c2(x)=2NLπDL∑i=1NXβiDL2τi28γi2+t0DLτi22γi12sin⁡2ψ1,i+12+12sin⁡2ψ2,i+12+12sin⁡2ψ3,i+12+2xMBNLπDLβbDL2τb28γb2+t0DLτb22γb+NLπDL∑i=1NXDL2min⁡1γi2,1γi+128+DLt0min⁡1γi,1γi+122;

coating costs, C3, depending on the outer surface area of all tubes, c3: c3(x)=2NLπDL∑i=1NXβiLi2cos⁡2Φp+Ri+Ri+12sin⁡2ϑ2+xMBNLπDLβb2R1sin⁡ϑ2+NLπDLLcos⁡Φs;

costs for the transition piece, C4, proportional to the product of head radius and number of jacket legs, c4: c4(x)=NLRfootξ;

expenses for transport, C5, expressed by the mass-dependent factor, c5: c5(x)=c1(x);

and installation costs, C6, modeled by a factor only depending on the number of jacket legs, c6: c6(x)=NL.

 Fixed expenses, C7, are not dependent on any jacket parameter at all. Therefore, the factor, c7, simply takes c7(x)=1. In these equations, ϑ is the angle enclosed by two jacket legs: ϑ=2πNL. Bay heights, Li, intermediate bay heights, Lm,i, radii, Ri, and intermediate radii, Rm,i, are calculated by the following equations: Li=L-LOSG-LTP∑n=1NXqn-i,Lm,i=LiRiRi+Ri+1,Ri=Rfoot-tan⁡ΦsLOSG+∑n=1i-1Ln,Rm,i=Rfoot-tan⁡ΦsLOSG+∑n=1i-1Ln+Lm,i, with the spatial batter angle, Φs: Φs=arctan⁡Rfoot1-ξL. The interconnecting tube angles, ψ1,i, ψ2,i, and ψ3,i, are ψ1,i=π2-arctan⁡Rfoot1-ξsin⁡ϑ2cos⁡ΦpL-arctan⁡Lm,iRisin⁡ϑ2cos⁡Φp,ψ2,i=π2+arctan⁡Rfoot1-ξsin⁡ϑ2cos⁡ΦpL-arctan⁡Lm,iRisin⁡ϑ2cos⁡Φp, ψ3,i=2arctan⁡Lm,iRisin⁡ϑ2cos⁡Φp, with the planar batter angle, Φp: Φp=arctan⁡Rfoot1-ξsin⁡ϑ2L. γi, βi, and τi represent the ratios of leg radius-to-thickness, brace-to-leg diameter, and brace-to-leg thickness of the ith bay, respectively, obtained by linear stepwise interpolation and counted upwards.

The cost modeling is based on several simplifications and assumptions. The mass-proportional modeling of material costs, C1, is straightforward. Fabrication costs, C2, mainly arise from welding and grinding processes. Although the actual manufacturing processes are quite complex, the entire volume of welds can be regarded as a measure of the actual costs. Coating costs, C3, are quite easy to determine by the outer surface area of all tubes, i.e., the area to be coated. There may be synergy effects when coating larger areas, but these are neglected. The expenses for the (stellar-type) transition piece, C4, are assumed to be proportional to the head radius and the number of legs. There are more detailed approaches for this purpose, but no design of the transition piece is performed, which requires a simple approach. The determination of transport costs, C5, is very difficult. In this work, a mass-dependent approach was selected, which is, however, a large simplification. The mass dependence reflects that barges have a limited transport capacity, which is at least to some extent mass-dependent or dependent on factors partially related to mass (like the space on the deck of the barge covered by the jacket). Installation costs, C6, cover both the material and the manufacturing of the foundation and the installation at the wind farm location. In the case of a pile foundation, these costs are mainly governed by the number of piles, which is equal to the number of legs. The fixed expenses, C7, are not vital for the solution of the optimization problem but are required to shift the costs to more realistic values by covering expenses for cranes, scaffolds, and so forth.

 Structural code checks (constraints)

To check jacket designs – i.e., sets of design variables – for validity concerning fatigue and extreme load resistance, structural design code checks are performed. The standards DNV GL RP-C203 for fatigue and NORSOK N-004 for ultimate limit state checks are adopted. Both are widely accepted for practical applications and were used to design the UpWind and INNWIND.EU reference jackets.

Commonly, the numerical demand of structural code checks is one of the main problems in jacket optimization. To cover the characteristics of environmental impacts on wind turbines, representative loads are to be used for the load assessment. This involves numerous load simulations to consider all load combinations that might occur, particularly in the fatigue case, where the excitation is extrapolated for the entire turbine lifetime. As not only the number of load simulations but also the duration (in the case of time domain simulations) correlates to a high demand in numerical capacity, most approaches deploy very simple load assumptions like one design load case per iteration, as already discussed. Altogether, a high numerical effort is required. Utilizing simplified load assumptions like equivalent static loads, where the substructure decoupled from the overlying structure and all interactions are neglected, depicts, however, a massive simplification in the case of a wide range of design variables. By contrast, a pure simulation-based optimization is not applicable due to the aforementioned reasons.

To face this issue, a surrogate modeling approach based on Gaussian process regression (GPR) is deployed. It was shown previously that good regression results can be obtained by GPR for this purpose. In addition, the regression process relies on a mathematical process that can be interpreted easily and adapted to prior knowledge of the underlying physics. In the present case, the procedure is as follows: a load set with a defined number of design load cases is the basis for structural code checks. The size of the load sets and parameters of environmental and operational conditions are predetermined so as to represent the loads on the turbine adequately. With these load sets, numerical simulations are performed with the aero–hydro–servo–elastic simulation code FAST to obtain output data for the input space of the surrogate model.

FASTv8 was used for this study.

As this procedure requires much computational effort, the input space is limited to 200 jacket samples (excluding validation samples) in each case as a basis for both surrogate models (fatigue and ultimate limit state),

All parameters of these jacket samples are given in the publication where the surrogate modeling approach was reported .

 obtained by a Latin hypercube sampling of the input space. In both cases, the results are vectors of output variables, where each element corresponds to a row in the matrix of inputs, comprising parameters of the input space. Both (input matrix and output vector) build the training data. For each new sample, the corresponding output (result of a structural code check) is evaluated by GPR.

For the background theory of GPR, the reader is referred to , which is the standard reference in this field.

 The specific surrogate models for the considered test problems were derived in a previous work , which revealed that a Matérn 5/2 kernel function is well-suited for the present application.

Fatigue limit state

The evaluation of fatigue limit state code checks requires many simulations considering design load cases (DLCs) 1.2 and 6.4 production load cases according to IEC 61400-3 . Under defined conditions (5 MW turbine, 50 m water depth, FINO3 environmental conditions), the required number of design load cases with respect to uncertainty was analyzed in previous papers . In these papers, a load set with 2048 design load cases was gradually reduced to smaller load sets. A reduced load set with 128 design load cases turned out to be a good compromise between accuracy, as the uncertainty arising from the load set reduction is acceptable in this case, and numerical effort, which is significantly smaller compared to the initial load set; i.e., considering two X-joint positions, the standard deviation of fatigue damage increases by a factor of approximately 4 in the case of a 16-fold load set reduction (from 2048 to 128 design load cases). The actual fatigue assessment involves time domain simulations, an application of stress concentration factors according to to consider the amplification of stresses due to the geometry of tubular joints, rain flow cycle counting, and a lifetime prediction by S-N curves and linear damage accumulation.

It has to be stated that there are several ways to determine stress concentration factors for tubular joints. This is the approach proposed by the standard DNV GL RP-C203 .

The output value hFLS is the most critical fatigue damage among all damage values of the entire jacket (evaluated in eight circumferential points around each weld), normalized by the calculated damage at design lifetime. A design lifetime of 30 years is assumed, from which 25 years are the actual lifetime of the turbine and 5 years are added to consider malicious fatigue damage during the transport and installation process. Moreover, a partial safety factor of 1.25 is considered in the fatigue assessment.

 Ultimate limit state

The standard IEC 61400-3 requires several design load cases to perform structural code checks for the ultimate limit state. However, not every design load case is critical for the design of a jacket substructure. The relevant ones were analyzed and found to be DLC 1.3 (extreme turbulence during production), 1.6 (extreme sea state during production), 2.3 (grid loss fault during production), 6.1 (extreme sea state during idle), and 6.2 (extreme yaw error during idle) for a turbine with a rated power of 5 MW, under FINO3 environmental conditions and a water depth of 50 m. Extreme load parameters are derived by the block maximum method see, where the environmental data are divided into many segments featuring similarly distributed data. From this data set, the maximum values are extracted. Based on these maxima, return values (as required by IEC 61400-3) of environmental states are computed. To conduct the structural code checks for the ultimate limit state, time domain simulations are performed and evaluated with respect to the extreme load of the member, where the highest utilization ratio occurs. The result hULS is a value that approaches 0 in the case of infinite extreme load resistance and 1 in the case of equal resistance and loads, implying that values greater than 1 are related to designs not fulfilling the ultimate limit state code check. The procedure considers combined loads with axial tension, axial compression, and bending, with and without hydrostatic pressure, which may lead to failure modes like material yielding, overall column buckling, local buckling, or any combination of these. A global buckling check is not performed in this study, as it is known to be uncritical for jacket substructures .

 Optimization approach and solution methods

The optimization problem incorporates a mixed-integer formulation (due to discrete numbers of legs and bays of the jacket). In order to address this issue, the mixed-integer problem is transferred to multiple continuous problems by solving solutions with a fixed number of legs and bays. As only a few combinations of these discrete variables are regarded as realistic solutions for practical applications, this procedure leads to a very limited number of subproblems but eases the mathematical optimization process significantly. Furthermore, the optimization problem is generally non-convex; i.e., a local minimum in the feasible region satisfying the constraints is not necessarily a global solution. This is addressed by repeating the optimization with multiple starting points.

The development of new or improved optimization methods to solve the numerical optimization problem is not in the scope of this work because there are methods presented in the literature that are known to be suitable for this purpose. Meta-heuristic algorithms like genetic algorithms or particle swarm optimization are not considered in this work because they are known to be slow. With regard to efficiency and accuracy, two methods are regarded as the most powerful for optimization involving nonlinear constraints: sequential quadratic programming (SQP) and interior-point (IP) methods . SQP methods are known to be efficient, when the numbers of constraints and design variables are of the same order of magnitude. An advantage is that these methods usually converge better when the problem is badly scaled. In theory, IP methods have better convergence properties and often outperform SQP methods on large-scale or sparse problems. In this work, both approaches are used to solve the jacket optimization problem.

The function fmincon in MATLAB R2017b was used for this study.

They are outlined briefly in the following.

Sequential quadratic programming method

In principle, SQP can be seen as an adaption of Newton's method to nonlinear constrained optimization problems, computing the solution of the Karush–Kuhn–Tucker equations (necessary conditions for constrained problems). Here, a common approach is deployed, based on the works of , , and . In the first step, the Hessian of the so-called Lagrangian (a term incorporating the objective and the sum of all constraints weighted by Lagrange multipliers) is approximated by the Broyden–Fletcher–Goldfarb–Shanno method . In the next step, a quadratic programming subproblem is built, where the Lagrangian is approximated by a quadratic term and linearized constraints. This subproblem can be solved by any method able to solve quadratic programs. An active-set method described by is deployed for this task. The procedure is repeated until convergence is reached.

 Interior-point method

IP methods are barrier methods; i.e., the objective is approximated by a term that incorporates a barrier term, expressed by a sum of logarithmized slack variables. The actual problem itself, just like in SQP, is solved as a sequence of subproblems. In this work, an approach is deployed, which may switch between line search and trust region methods to approximated problem, depending of the success of each step. If the line search step fails, i.e., when the projected Hessian is not definitively positive, the algorithm performs a trusted region step, where the method of conjugate gradients is deployed. The algorithm is described in detail by .

 Jacket comparison study

In this section, the proposed approach is applied to find and compare optimal jacket designs for the NREL 5 MW reference turbine . The environmental conditions are adopted from measurements recorded at the research platform FINO3 in the German North Sea.

Reference turbine

The NREL 5 MW reference turbine, which was published almost 1 decade ago as a proposal to establish a standardized turbine for scientific purposes, is still an object of many studies in the literature dealing with intermediate- to high-power offshore wind applications. In fact, the market already provides turbines with 8 MW and aims for even higher ratings. Choosing this reference turbine is motivated by its excellent documentation and accessibility.

The rotor has a hub height of 90 m, and the rated wind speed is 11.4 m s-1, where the rotor speed is 12.1 min-1. This is equal to 1P and 3P excitations of 0.2 and 0.6 Hz, respectively. The critical first fore–aft and side–side bending eigenfrequencies of the entire structure are about 0.35 Hz and do not differ very much when considering only reasonable structural designs for the jacket because the modal behavior is strongly driven by the relatively soft tubular tower.

 Environmental conditions and design load sets

Due to excellent availability, the environmental data are derived from measurements taken from the offshore research platform FINO3, located in the German North Sea close to the wind farm “alpha ventus”. Compared to the environmental conditions documented in the UpWind design basis , the FINO3 measurements are much more comprehensive and allow for a better estimation of probability density functions as inputs for the determination of probabilistic loads . The probabilistic load set, which is based on probability density functions of environmental state parameters and reduced in size compared to full load sets used by industrial wind turbine designers, was described in recent studies . However, there are two drawbacks that have to be mentioned when using this data. First, the FINO3 platform was built at a location with quite a shallow water depth of 22 m, though the jacket is supposed to be an adequate substructure for water depths above 40 m and the design water depth in this study is 50 m. Nevertheless, this procedure was also performed in the UpWind project for the design of the OC4 jacket, where the K13 deep-water site was considered. Second, the soil properties of the Offshore Code Comparison Collaboration (OC3) are adopted to compute foundation inertias and stiffnesses, as these values are unknown for the FINO3 location. Moreover, it is assumed that the structural behavior of the OC4 jacket pile foundation is valid for all jacket designs, even with varying leg diameters and thicknesses.

 Boundaries of design variables and other parameters

The boundaries are chosen conservatively by means of quite narrow design variable ranges (see Table ), i.e., meaningful parameters that do not exhaust the possible range given by the structural code checks, in a realistic range around the values of the OC4 jacket . Only three- or four-legged structures with three, four, and five bays are regarded as valid solutions for this study. The fixed design variables are, if possible, adopted from the OC4 jacket, which can be seen as a kind of reference structure in this case. The material is steel (S355), with a Young's modulus of 210 GPa, a shear modulus of 81 GPa, and a density of 7850 kg m-3. According to , an S-N curve with an endurance stress limit of 52.63×106 N m-2 at ×107 cycles and slopes of 3 and 5 before and after endurance limit (curve T), respectively, is applied. The cost model parameters or unit costs, respectively, are adopted from the mean values given in and set to a1=1.0 kg-1 (material), a2=4.0×106 m-3 (fabrication), a3=1.0×102 m-2 (coating), a4=2.0×104 m-1 (transition piece), a5=2.0 kg-1 (transport), a6=2.0×105 (installation), and a7=1.0×105 (fixed). With these values, the cost function returns a dimensionless value, also interpretable as capital expenses in EUR.

Boundaries of jacket model parameters for design of experiments. Topological, tube sizing, and material parameters are separated into groups; single values mean that the corresponding value is held constant.

Parameter Description Lower boundary Upper boundary NL Number of legs 3 4 NX Number of bays 3 5 Rfoot Foot radius 6.792 m 12.735 m ξ Head-to-foot radius ratio 0.533 0.733 L Entire jacket length 70.0 m LMSL Transition piece elevation over mean sea level 20.0 m LOSG Lowest leg segment height 5.0 m LTP Transition piece segment height 4.0 m q Ratio of two consecutive bay heights 0.640 1.200 xMB Mud brace flag true (1) DL Leg diameter 0.960 m 1.440 m γb Leg radius-to-thickness ratio (bottom) 12.0 18.0 γt Leg radius-to-thickness ratio (top) 12.0 18.0 βb Brace-to-leg diameter ratio (bottom) 0.533 0.800 βt Brace-to-leg diameter ratio (top) 0.533 0.800 τb Brace-to-leg thickness ratio (bottom) 0.350 0.650 τt Brace-to-leg thickness ratio (top) 0.350 0.650 E Material Young's modulus 2.100×1011 N m-2 G Material shear modulus 8.077×1010 N m-2 ρ Material density 7.850×103 kg m-3 Results and discussion

To resolve the mixed-integer formulation of the optimization problem into continuous problems, six subproblems with three legs and three bays (NL=3, NX=3), three legs and four bays (NL=3, NX=4), three legs and five bays (NL=3, NX=5), four legs and three bays (NL=4, NX=3), four legs and four bays (NL=4, NX=4), and four legs and five bays (NL=4, NX=5) were solved using the SQP and IP methods. Therefore, multiple solutions are discussed and compared in the following. The optimization problem is non-convex; i.e., a local minimum in the feasible region satisfying the constraints is not necessarily a global solution. In theory, both algorithms converge from remote starting points. However, to guarantee global convergence to some extent, all six combinations of fixed-integer variables were solved using 100 randomly chosen starting points. Installation costs and fixed expenses were excluded from the objective function and included again after the optimization procedure because these terms do not have an effect on the individual optimization problems.

The values shown in the following include all cost terms. The exclusion is only performed during optimization.

Gradients were computed by finite differences. The optimization terminated, when the first-order optimality and feasibility measures were both less than 1×10-6. There was no limit to the maximum number of iterations.

The optimal solutions of all six subproblems do not depend on the starting point when using both optimization methods because there is only one array of optimal design variables in each case. The convergence behavior of both optimization methods is illustrated in Fig. , where the OC4 jacket with varying numbers of legs and bays was assumed as the starting point. This structure has a foot radius, Rfoot, of 8.79 m, a head-to-foot radius ratio, ξ, of 0.67, and a ratio of two consecutive bay heights, q, of 0.8. Moreover, it has a leg diameter, DL, of 1.2 m, and entirely constant tube dimensions from bottom to top, i.e., leg radius-to-thickness ratios, γb and γt, of 15, brace-to-leg diameter ratios, βb and βt, of 0.5, and brace-to-leg diameter ratios, τb and τt, of 0.5. The optimization process needed between 30 and 40 iterations using the SQP method and between 50 and 70 iterations using the IP method to converge. It is worth mentioning that the maximum constraint violation (feasibility) of the three-legged designs was higher at the beginning of the optimization process but converges stably. For the same reason, the four-legged designs have a higher improvement potential compared to the initial solution. The accuracy obtained by both methods is similar. The solutions are all feasible because they fulfill the Karush–Kuhn–Tucker conditions, and all constraint violations are around zero. Therefore, the optima are probably global optima for the given design variable boundaries.

Function and feasibility (maximum constraint violation) values during the optimization procedure of all six subproblems (blue line with circles: jacket with three legs and three bays; red line with triangles: jacket with three legs and four bays; brown line with diamonds: jacket with three legs and five bays; black line with pentagons: jacket with four legs and three bays; violet line with half-filled circles: jacket with four legs and four bays; green line with half-filled diamonds: jacket with four legs and five bays). The starting point (iteration “0”) is the OC4 jacket with a varying number of legs and bays in all cases. One iteration involves 11 evaluations of the objective function and the nonlinear constraints.

Topologies of optimal solutions x*. All images are displayed at the same scale. Line widths are not correlated to tube dimensions.

Optimal solutions of design variables x* obtained by the sequential quadratic programming method for fixed values of NL and NX.

Optimal solution x* NL 3 3 3 4 4 4 NX 3 4 5 3 4 5 Rfoot in m 12.735 12.735 12.735 10.894 10.459 10.549 ξ 0.533 0.533 0.533 0.533 0.533 0.533 q 0.937 0.941 0.936 0.813 0.809 0.977 DL in m 1.021 1.021 1.023 0.960 0.960 0.960 βb 0.800 0.800 0.800 0.800 0.799 0.787 βt 0.800 0.800 0.800 0.800 0.800 0.800 γb 12.000 12.000 12.000 12.680 12.259 12.000 γt 16.165 16.029 15.928 18.000 18.000 18.000 τb 0.513 0.505 0.493 0.497 0.493 0.478 τt 0.472 0.466 0.454 0.383 0.387 0.383 Overall mass in t 423 444 467 412 426 439 f(x*)=log⁡10Ctotal(x*) 6.452 6.472 6.493 6.487 6.500 6.514 h1(x*)=hFLS(x*)-1 1.172×10-10 3.966×10-11 1.151×10-10 1.450×10-10 -1.056×10-10 -1.721×10-10 h2(x*)=hULS(x*)-1 7.819×10-10 2.678×10-10 1.093×10-10 3.978×10-10 3.980×10-10 5.995×10-10

The optimal solutions obtained by the sequential quadratic programming method are illustrated in Table .

As the accuracy of the SQP and IP methods is similar here, only results obtained by the SQP method are shown in the following.

Additionally, the topologies of all optimal solutions are shown in Fig. . With respect to the constraints and assumptions of this study (5 MW turbine, 50 m water depth, given environmental conditions and cost parameters), jackets with three legs are beneficial in terms of capital expenses. The three-legged jacket with three bays (NL=3, NX=3) is the best solution, i.e., is related to the lowest total expenditures, among the considered jackets. The solutions show some interesting specialties. The foot radii, Rfoot, are at the upper boundaries in the case of the three-legged structures, while the head-to-foot radius ratios, ξ, are at the lower boundaries. Probably this arises from the combination of cost function and nonlinear constraints, where a large foot radius is quite beneficial because it generally provides a higher load capacity, while a small head radius is favorable due to lower transition piece costs. In the four-legged case, the foot radii are lower but still relatively high. In any case, it seems to be beneficial, when the ratio of two consecutive bay heights, q, is slightly below 1 (lower bays are higher than upper bays). Concerning tube dimensions, the leg diameters, DL, are relatively small, in the case of the four-legged jackets even at the lower boundary. The structural load capacity is established by high brace diameters (represented by design variables βb and βt, values at the bottom and top of the structures both at upper boundaries). The brace thicknesses, represented by τb and τt, show intermediate values in the range of design variables, while the values for τt are higher in the case of three-legged designs. Moreover, the structural resistance is strongly driven by the leg thicknesses. While the optimal values of γb are low in each case, implying high leg thicknesses at the jacket bottom, the values of γt are much higher. The impact of all design variables on the objective function is easier to understand when the sensitivities of cost model terms to variations in design variables are considered. In Fig. , each subplot shows the variation in the total costs, Ctotal, and the cost function terms C1 (proportional to C5), C2, C3, and C4 due to a 1 % one-at-a-time variation in each continuous design variable in three different phases of the optimization process (initial, intermediate, and final phase). The terms C6 and C7 are not impacted by any continuous design variable and therefore not considered. For instance, a 1 % increase in the foot radius, Rfoot, causes increasing material costs of ΔC1=0.14 %, evaluated for the initial design, but increasing material costs of ΔC1=0.26 %, evaluated for the optimal design. Therefore, the sensitivity of this cost term varies during the optimization process. In contrast, the variation in transition piece expenses does not change (which is reasonable because this term only depends linearly on the number of legs, NL, the foot radius, Rfoot, and the head-to-foot radius ratio, ξ). In general, Fig.  shows that there is no design variable with a strongly varying impact on any term of the cost function. It can also be concluded that tube sizing variables impact the costs much more strongly than topological variables, disregarding the number of legs and bays. Among the considered design variables, the leg diameter, DL, and leg radius-to-thickness ratios, γb and γt, are design-driving (together with the number of legs, NL) due to a significant impact both on the costs and on the structural code checks. In addition, an interesting specialty is featured by the cost term C4, which is only impacted by topological design variables, more precisely the foot radius, Rfoot, and the head-to-foot radius ratio, ξ. As a large foot radius, Rfoot, is needed to establish structural resistance, this cost term penalizes large head-to-foot radius ratios, ξ. For this reason, this value is at the lower boundary for all design solutions.

Variations in total costs, ΔCtotal, and cost function terms ΔC1 (material), ΔC2 (manufacturing), ΔC3 (coating), and ΔC4 (transition piece) due to 1 % one-at-a-time variations in design variables (subplots) in %. Derivatives were computed for the initial design (red bars), an intermediate design after 15 iterations (blue bars), and the optimal design (green bars) of the three-legged structure with three bays (NL=3, NX=3).

Regarding the costs of the jackets, the best solution with three legs and three bays is related to capital expenses of 106.452= 2 831 000. Altogether, this is a meaningful value and the designs are not far off from structural designs that are known from practical applications because it has already been reported in the literature that three-legged designs may be favorable in terms of costs and three-legged structures have already been built. However, the other solutions are more expensive but not completely off. As there is some uncertainty in the unit costs, the other jackets may also be reasonable designs with slightly different boundaries. A more detailed cost breakdown is given in Fig. , which shows the cost contributions of all six structures and where the actual cost savings come from. The lightest structure is the four-legged jacket with three bays, while the three-legged jacket with five bays is the heaviest one, which is illustrated by the expenses for material and transport according to the cost model used for this study. Nevertheless, the mass of all structures is quite similar. Other than expected, the jacket with the lowest expenditures for manufacturing is also the four-legged one with three bays and not the three-legged jacket with three bays, which has the least number of joints. The three-legged structures benefit – from the economic point of view – mainly from lower expenses for coating, the transition piece, and, most distinctly, installation costs. In total, these contributions add up to lower costs of the three-legged jackets, except for the one with five bays (106.493= 3 112 000), which is more expensive than the four-legged one with three bays (106.487= 3 069 000). The most expensive jackets are the four-legged ones with four (106.500= 3 162 000) and five (106.514= 3 266 000) bays, where the latter is about 15 % more expensive than the best solution among the six sub-solutions. A reasonable option may also be the jacket with three legs and four bays, which features a total cost value of 106.472= 2 965 000. In total, there is no jacket that is far too expensive compared to the others. It is indeed imaginable to find an appropriate application for each one.

Expenses comparison of optimal solutions of three-legged jacket with three bays (blue bars), three-legged jacket with four bays (red bars), three-legged jacket with five bays (brown bars), four-legged jacket with three bays (gray bars), four-legged jacket with four bays (violet bars), and four-legged jacket with five bays (green bars).

From the computational point of view, the optimization procedure based on surrogate models is very efficient. The numbers of iterations needed to find an optimal solution (from about 30 to 40 using the SQP method and from about 50 to 70 using the IP method) are related to computation times of about 15 to 30 min on a single core of a work station with an Intel Xeon E5-2687W v3 central processing unit and 64 GB random access memory. Compared to simulation-based approaches, this can be regarded as very fast. The number of iterations may be decreased, when using analytical gradients of the objective function because using finite differences is generally more prone to numerical errors but is not vital at this level of computational expenses. It has to be pointed out that the training data set of the surrogate models required 200×128= 25 600 time domain simulations in the fatigue and 200×10=2000 in the ultimate limit state case, thus 27 600 simulations in total, excluding validation samples. However, for the computation of the training data, a compute cluster was utilized, which allows for the computation of many design load cases in parallel. Therefore, the presented approach based on GPR allows for outsourcing computationally expensive simulations on high-performance clusters, while the closed-loop optimization, which cannot be parallelized completely, can be run on a workstation with lower computational capacity.

The question remains what happens when some cost terms are neglected. An associated question is how the approach performs compared to a pure mass-dependent one, which can be regarded as state of the art in jacket optimization. For this purpose, all unit costs except a1 were set to zero and the optimization procedure was repeated using the sequential quadratic programming method. The results, including optimal design variables and resulting values of objective and constraint functions, are shown in Table . Under these assumptions, the four-legged jackets are better (in terms of minimal mass) than the three-legged ones. Interestingly, similar design variables are obtained when comparing these values to the ones obtained by the more comprehensive cost model in Table , particularly in the case of the three-legged jackets. The resulting objective function values are, in comparison, similar to the material costs in Fig. . In other words, a pure mass-dependent cost function approach yields approximately proportional costs, when the installation costs (depending on the number of legs) are considered, and similar designs. The reason for this is that all cost terms C1…C5 depend in some way on the tube dimensions and the topology does not impact the costs to a great extent, as seen in Fig. . Indeed, the largest proportion of costs is purely mass-dependent, as the factors c1 and c5 are the mass of the structure. Therefore, the proposed cost model can lead to more accurate results, but a mass-dependent approach would be sufficient to draw the same conclusions.

Optimal solutions of design variables x* obtained by the sequential quadratic programming method for fixed values of NL and NX using a pure mass-dependent objective function.

Optimal solution x* NL 3 3 3 4 4 4 NX 3 4 5 3 4 5 Rfoot in m 12.735 12.735 12.735 12.735 12.735 12.735 ξ 0.533 0.533 0.533 0.533 0.533 0.533 q 1.062 0.987 0.936 1.200 1.200 1.178 DL in m 1.025 1.023 1.023 0.960 0.960 0.960 βb 0.800 0.800 0.800 0.730 0.757 0.800 βt 0.800 0.800 0.800 0.800 0.800 0.800 γb 12.000 12.000 12.000 13.194 13.318 12.000 γt 16.459 16.250 15.928 18.000 18.000 18.000 τb 0.509 0.502 0.493 0.510 0.470 0.443 τt 0.472 0.466 0.454 0.386 0.361 0.350 Overall mass in t 423 444 467 404 409 454 f(x*)=log⁡10Ctotal(x*) 5.627 5.647 5.669 5.606 5.612 5.657 h1(x*)=hFLS(x*)-1 -7.149×10-12 6.767×10-12 1.262×10-12 5.047×10-13 2.140×10-12 -1.017×10-12 h2(x*)=hULS(x*)-1 4.367×10-11 2.961×10-11 5.087×10-12 2.693×10-12 3.865×10-12 9.948×10-13 Benefits and limitations of the approach

With respect to the state of the art, the present approach can be regarded as the first one addressing the jacket optimization problem holistically, which incorporates four main improvements: a detailed geometry model with both topological and tube sizing design variables; an analytical cost model based on the main jacket cost contributions; sophisticated load assumptions and assessments; and a treatment of results that considers the optimization problem to be a methodology for early design stages. All these points lead to a better understanding how to address the multidisciplinary design optimization problem and to much more reliable results.

However, some drawbacks and limitations remain, which have to be considered when dealing with the results of this study. In general, the approach is easy to use, also in industrial applications, but needs some effort in implementation. Furthermore, the present study does not incorporate a completely reliability-based design procedure, which is not beyond the means when using Gaussian process regression to perform structural code checks. However, the question of how safety factors can be replaced by a meaningful probabilistic design is still a matter of research, and it is quite simple to advance the present approach to a robust one. In order to reduce the numerical cost (in particular concerning the number of time domain simulations needed to sample the input design space for surrogate modeling of structural code checks), the number of design variables is limited. The application of GPR as a machine learning approach to evaluate structural code checks performs in a numerically fast way but requires numerous time domain simulations to generate training and validation data sets. This is beneficial when dealing with numerically expensive studies (as in this case) but might lead to a numerical overhead when only considering one jacket design. Care has to be taken when transferring the results to designs with a more sophisticated geometry. Moreover, the parameterization of cost and structural code check models is site- and turbine-dependent. Therefore, the outcome of this study might not be directly transferable to other boundaries but requires recalculations. In particular, the utilized design standards and structural code checks are known to be conservative. The cost model also has shortcomings that must be mentioned. Some costs are affected by uncertain or indeterminable impacts. There is a number of examples. Transport and installation costs are strongly dependent on the availability of barges or vessels. The uncertainty in weather conditions can affect transport and installation costs. Furthermore, the design may be directly impacted if production facilities are not available. All these effects are not considered in the cost model.

In addition, it is important to highlight again that this study does not provide a detailed design methodology but an approach to obtain preliminary decision guidance at the earliest wind farm planning stage. This is actually not a limitation but has to be considered when dealing with the results of this study. There are indeed many studies known from the literature that address the tube dimensioning problem in a larger extension. However, these approaches assume that the structural topology is always optimal, even in the case of significant variations in tube dimensions. For instance, all optimal jackets have a larger bottom width than the OC4 jacket, while the design-driving leg diameters are relatively small. This indicates that topological design variables with minor impact on costs are useful factors to establish the structural resistance.

 Conclusions

The present work began by introducing four main points to be considered in order to improve the state of the art in the field of jacket optimization. The first one, the treatment of the jacket design problem in terms of a holistic topology and tube sizing problem instead of a pure tube dimensioning problem, was addressed by a 20-parameter jacket model, of which 12 parameters are design variables. The second, important point leads to the utilization of a more complex (compared to mass-dependent) but easy to handle cost model. In order to face the challenging task of numerically efficient structural code check evaluations, a machine learning approach based on Gaussian process regression was applied as the third point. On this basis, gradient-based optimization was deployed to find optimal design solutions. Lastly, optimization results were considered differently compared to approaches presented in the literature. It was pointed out that the solution is not supposed to be the final design but a very good starting point to find an initial solution for exact tube dimensioning.

The conclusions of this work are manifold. From the numerical point of view, surrogate modeling seems – as matters stand today – to be the most promising approach enabling us to address the computationally very expensive jacket optimization problem efficiently because other approaches in the literature go along with massive simplifications, mainly in load assumptions. The optimization methods that were used to find the optimal solution seem to be appropriate for the given problem, even in terms of finding a global optimum. The present paper does not provide improvements of state-of-the-art gradient-based optimization, but active-set SQP and IP methods both converge efficiently and accurately for the given problem.

From the application-oriented point of view, it can be stated that three-legged jackets with only three bays depict the best solution (in terms of costs) for offshore turbines with about 5 MW rated power in 50 m water depth, which confirms the results from other studies in the literature. Due to the cost model, the additional load-bearing capacity gained by the extra leg of a four-legged structure cannot compensate for the higher costs arising from several cost factors directly related to the number of legs. By contrast, it is instead beneficial to increase the tube dimensions and maintain the number of structural elements at a minimum level. It was shown that the same results were obtained when using a mass-dependent cost function, also considering the number of jacket legs.

With regard to turbines with a higher rated power or installations in deeper waters, the proposed methodology might lead to the result that the best jacket solution for this case looks completely different. Before this can be analyzed, simulation tools need to be improved to enable the consideration of nonlinear effects for rotors with a very large diameter and innovative control strategies.

This work is based on structural code checks computed and provided in Häfele et al. (2018a).

 Nomenclature DLC Design load case IP Interior-point method SQP Sequential quadratic programming method Φp Planar (two-dimensional) batter angle Φs Spatial (three-dimensional) batter angle βb Brace-to-leg diameter ratio at bottom (jacket model parameter) βi Brace-to-leg diameter ratio in the ith bay βt Brace-to-leg diameter ratio at top (jacket model parameter) γb Leg radius-to-thickness ratio at bottom (jacket model parameter) γi Leg radius-to-thickness ratio in the ith bay γt Leg radius-to-thickness ratio at top (jacket model parameter) ξ Head-to-foot radius ratio (jacket model parameter) ρ Material density (jacket model parameter) ϑ Angle enclosed by two jacket legs τb Brace-to-leg thickness ratio at bottom (jacket model parameter) τi Brace-to-leg thickness ratio in the ith bay τt Brace-to-leg thickness ratio at top (jacket model parameter) ψ1,i Lower brace-to-leg connection angle in the ith bay ψ2,i Upper brace-to-leg connection angle in the ith bay ψ3,i Brace-to-brace connection angle in the ith bay Cj Expenses related to jth cost factor Ctotal Total capital expenses DL Leg diameter (jacket model parameter) E Material Young's modulus (jacket model parameter) G Material shear modulus (jacket model parameter) L Overall jacket length (jacket model parameter) LMSL Transition piece elevation over mean sea level (jacket model parameter) LOSG Lowest leg segment height (jacket model parameter) LTP Transition piece segment height (jacket model parameter) Li ith jacket bay height Lm,i Distance between the lower layer of K joints and the layer of X joints of the ith bay NL Number of legs (jacket model parameter) NX Number of bays (jacket model parameter) RFoot Foot radius (jacket model parameter) Ri ith jacket bay radius at lower K joint layer Rm,i Radius of the ith X joint layer aj jth unit cost cj jth cost factor f Objective function value h1 First inequality constraint value h2 Second inequality constraint value hFLS Maximal normalized tubular joint fatigue damage hULS Maximal extreme load utilization ratio q Ratio of two consecutive bay heights (jacket model parameter) x Array of design variables xlb Array of lower boundaries xMB Mud brace flag (jacket model parameter) xub Array of upper boundaries

The authors declare that they have no conflict of interest.

Acknowledgements

This work was supported by the computer cluster, which is funded by the Leibniz Universität Hannover, the Lower Saxony Ministry of Science and Culture (MWK), and the German Research Foundation (DFG). The publication of this article was funded by the open-access fund of Leibniz Universität Hannover. Edited by: Lars Pilgaard Mikkelsen Reviewed by: Lars Einar S. Stieng and one anonymous referee

 References Agarwal, P. and Manuel, L.: Load Extrapolation Methods for Offshore Wind Turbines, in: Offshore Technology Conference, Houston, Texas, USA, 10.4043/21001-MS, 2010. AlHamaydeh, M., Barakat, S., and Nasif, O.: Optimization of Support Structures for Offshore Wind Turbines Using Genetic Algorithm with Domain-Trimming, Math. Probl. Eng., 2017, 1–14, 10.1155/2017/5978375, 2017. Biggs, M. C.: Constrained minimization using recursive quadratic programming, in: Towards global optimization, North-Holland, Cagliari, Italy, 341–349, 1975. BVGassociates: Offshore wind: Industry's journey to GBP 100/MWh – Cost breakdown and technology transition from 2013 to 2020, Tech. rep., 2013. Chew, K. H., Tai, K., Ng, E. Y. K., Muskulus, M., and Zwick, D.: Offshore Wind Turbine Jacket Substructure: A Comparison Study Between Four-Legged and Three-Legged Designs, Journal of Ocean and Wind Energy, 1, 74–81, 2014. Chew, K.-H., Tai, K., Ng, E., and Muskulus, M.: Optimization of Offshore Wind Turbine Support Structures Using an Analytical Gradient-based Method, Energy Proced., 80, 100–107, 10.1016/j.egypro.2015.11.412, 2015. Chew, K.-H., Tai, K., Ng, E., and Muskulus, M.: Analytical gradient-based optimization of offshore wind turbine substructures under fatigue and extreme loads, Mar. Struct., 47, 23–41, 10.1016/j.marstruc.2016.03.002, 2016. Damiani, R.: JacketSE: An Offshore Wind Turbine Jacket Sizing Tool, Technical Report NREL/TP-5000-65417, National Renewable Energy Laboratory, Golden, CO, USA, 2016. Damiani, R., Dykes, K., and Scott, G.: A comparison study of offshore wind support structures with monopiles and jackets for U.S. waters, J. Phys. Conf. Ser., 753, 092003, 10.1088/1742-6596/753/9/092003, 2016. DNV GL AS: Fatigue design of offshore steel structures, Recommended Practice DNVGL-RP-C203, 2016. Efthymiou, M.: Development of SCF formulae and generalised influence functions for use in fatigue analysis, in: Proceedings of the Conference on Recent Developments in Tubular Joints Technology, 4–5 October 1988, Surrey, UK, 1–13, 1988. Fischer, T., de Vries, W., and Schmidt, B.: Upwind Design Basis, Tech. rep., 2010. Fletcher, R.: Practical methods of optimization, Wiley, Chichester, UK, New York, USA, 2nd edn., 1987. Gill, P. E., Murray, W., and Wright, M. H.: Practical optimization, Academic Press, London, UK, New York, USA, 1981. Häfele, J. and Rolfes, R.: Approaching the ideal design of jacket substructures for offshore wind turbines with a Particle Swarm Optimization algorithm, in: Proceedings of the Twenty-sixth (2016) International Offshore and Polar Engineering Conference, 26 June–2 July 2016, Rhodes, Greece, 156–163, 2016. Häfele, J., Gebhardt, C. G., and Rolfes, R.: Innovative design of a hybrid-type jacket for 10MW turbines, Tech. Rep. INNWIND.EU Deliverable D4.3.5, 2016. Häfele, J., Hübler, C., Gebhardt, C. G., and Rolfes, R.: Efficient Fatigue Limit State Design Load Sets for Jacket Substructures Considering Probability Distributions of Environmental States, in: Proceedings of the Twenty-seventh (2017) International Offshore and Polar Engineering Conference, 25–30 June 2017, San Francisco, CA, 167–173, 2017. Häfele, J., Damiani, R. R., King, R. N., Gebhardt, C. G., and Rolfes, R.: A systematic approach to offshore wind turbine jacket predesign and optimization: geometry, cost, and surrogate structural code check models, Wind Energ. Sci., 3, 553–572, 10.5194/wes-3-553-2018, 2018a. Häfele, J., Hübler, C., Gebhardt, C. G., and Rolfes, R.: A comprehensive fatigue load set reduction study for offshore wind turbines with jacket substructures, Renew. Energ., 118, 99–112, 10.1016/j.renene.2017.10.097, 2018b. Han, S. P.: A globally convergent method for nonlinear programming, J. Optimiz. Theory. App., 22, 297–309, 10.1007/BF00932858, 1977. Hübler, C., Gebhardt, C. G., and Rolfes, R.: Development of a comprehensive database of scattering environmental conditions and simulation constraints for offshore wind turbines, Wind Energ. Sci., 2, 491–505, 10.5194/wes-2-491-2017, 2017. International Electrotechnical Commision: Wind turbines – Part 3: Design requirements for offshore wind turbines, Tech. Rep. IEC-61400-3:2009, 2009. Jonkman, J. and Musial, W.: Offshore Code Comparison Collaboration (OC3) for IEA Task 23 Offshore Wind Technology and Deployment, Technical Report NREL/TP-5000-48191, National Renewable Energy Laboratory, Golden, CO, USA, 2010. Jonkman, J., Butterfield, S., Musial, W., and Scott, G.: Definition of a 5-MW Reference Wind Turbine for Offshore System Development, Technical Report NREL/TP-500-38060, National Renewable Energy Laboratory, Golden, CO, USA, 2009. Kaveh, A. and Sabeti, S.: Structural optimization of jacket supporting structures for offshore wind turbines using colliding bodies optimization algorithm, Struct. Des. Tall Spec., 27, e1494, 10.1002/tal.1494, 2018. Mone, C., Hand, M., Bolinger, M., Rand, J., Heimiller, D., and Ho, J.: 2015 Cost of Wind Energy Review, Technical Report NREL/TP-6A20-66861, National Renewable Energy Laboratory, Golden, CO, USA, 2017. Muskulus, M. and Schafhirt, S.: Design Optimization of Wind Turbine Structures – A Review, Journal of Ocean and Wind Energy, 1, 12–22, 2014. National Wind Technology Center Information Portal: FAST v8, available at: https://nwtc.nrel.gov/FAST8 (last access: 28 August 2018), 2016. Nocedal, J. and Wright, S. J.: Numerical optimization, Springer series in operations research, Springer, New York, USA, 2nd edn., 2006. NORSOK: Design of steel structures, Standard N-004 Rev. 2, 2004. Oest, J., Sørensen, R., T. Overgaard, L. C., and Lund, E.: Structural optimization with fatigue and ultimate limit constraints of jacket structures for large offshore wind turbines, Struct. Multidiscip. O., 55, 779–793, 10.1007/s00158-016-1527-x, 2016. Oest, J., Sandal, K., Schafhirt, S., Stieng, L. E. S., and Muskulus, M.: On gradient-based optimization of jacket structures for offshore wind turbines, Wind Energy, 21, 953–967, 10.1002/we.2206, 2018. Popko, W., Vorpahl, F., Zuga, A., Kohlmeier, M., Jonkman, J., Robertson, A., Larsen, T. J., Yde, A., Saetertro, K., Okstad, K. M., Nichols, J., Nygaard, T. A., Gao, Z., Manolas, D., Kim, K., Yu, Q., Shi, W., Park, H., Vasquez-Rojas, A., Dubois, J., Kaufer, D., Thomassen, P., de Ruiter, M. J., Peeringa, J. M., Zhiwen, H., and von Waaden, H.: Offshore Code Comparison Collaboration Continuation (OC4), Phase I – Results of Coupled Simulations of an Offshore Wind Turbine with Jacket Support Structure, Journal of Ocean and Wind Energy, 1, 1–11, 2014. Powell, M. J. D.: The Convergence of Variable Metric Methods for Nonlinearly Constrained Optimization Calculations, Nonlinear Programming, 3, 27–63, 10.1016/B978-0-12-468660-1.50007-4, 1978a. Powell, M. J. D.: A fast algorithm for nonlinearly constrained optimization calculations, in: Numerical Analysis, Dundee, Scotland, 630, 144–157, 10.1007/BFb0067703, 1978b. Rasmussen, C. E. and Williams, C. K. I.: Gaussian processes for machine learning, MIT Press, Cambridge, MA, USA, 3rd print edn., 2008. Sandal, K., Latini, C., Zania, V., and Stolpe, M.: Integrated optimal design of jackets and foundations, Mar. Struct., 61, 398–418, 10.1016/j.marstruc.2018.06.012, 2018. Schafhirt, S., Zwick, D., and Muskulus, M.: Reanalysis of Jacket Support Structure for Computer-Aided Optimization of Offshore Wind Turbines with a Genetic Algorithm, Journal of Ocean and Wind Energy, 1, 209–216, 2014. Schafhirt, S., Zwick, D., and Muskulus, M.: Two-stage local optimization of lattice type support structures for offshore wind turbines, Ocean Eng., 117, 163–173, 10.1016/j.oceaneng.2016.03.035, 2016. Seidel, M.: Jacket substructures for the REpower 5M wind turbine, in: Conference Proceedings European Offshore Wind 2007, Berlin, Germany, 2007. Smith, A., Stehly, T., and Musial, W.: 2014-2015 Offshore Wind Technologies Market Report, Technical Report NREL/TP-5000-64283, National Renewable Energy Laboratory, Golden, CO, USA, 2015. Stolpe, M. and Sandal, K.: Structural optimization with several discrete design variables per part by outer approximation, Struct. Multidiscip. O., 57, 2061–2073, 10.1007/s00158-018-1941-3, 2018. The Crown Estate: Offshore Wind Cost Reduction – Pathways Study, Tech. rep., 2012. van Kuik, G. A. M., Peinke, J., Nijssen, R., Lekou, D., Mann, J., Sørensen, J. N., Ferreira, C., van Wingerden, J. W., Schlipf, D., Gebraad, P., Polinder, H., Abrahamsen, A., van Bussel, G. J. W., Sørensen, J. D., Tavner, P., Bottasso, C. L., Muskulus, M., Matha, D., Lindeboom, H. J., Degraer, S., Kramer, O., Lehnhoff, S., Sonnenschein, M., Sørensen, P. E., Künneke, R. W., Morthorst, P. E., and Skytte, K.: Long-term research challenges in wind energy – a research agenda by the European Academy of Wind Energy, Wind Energ. Sci., 1, 1–39, 10.5194/wes-1-1-2016, 2016. Vemula, N., DeVries, W., Fischer, T., Cordle, A., and Schmidt, B.: Design solution for the upwind reference offshore support structure, UpWind deliverable report D4.2.6, 2010. von Borstel, T.: Design Report Reference Jacket, Tech. Rep. INNWIND.EU deliverable report D4.3.1, 2013. Waltz, R. A., Morales, J. L., Nocedal, J., and Orban, D.: An interior algorithm for nonlinear optimization that combines line search and trust region steps, Math. Program., 107, 391–408, 10.1007/s10107-004-0560-5, 2006. Zwick, D. and Muskulus, M.: Simplified fatigue load assessment in offshore wind turbine structural analysis, Wind Energy, 19, 265–278, 10.1002/we.1831, 2016.