Research article 22 Jan 2019
Research article  22 Jan 2019
A comparison study on jacket substructures for offshore wind turbines based on optimization
 Leibniz Universität Hannover/ForWind, Institute of Structural Analysis, Appelstr. 9a, 30167 Hanover, Germany
 Leibniz Universität Hannover/ForWind, Institute of Structural Analysis, Appelstr. 9a, 30167 Hanover, Germany
Correspondence: Jan Häfele (j.haefele@isd.unihannover.de)
Hide author detailsCorrespondence: Jan Häfele (j.haefele@isd.unihannover.de)
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 interiorpoint method), allowing for the estimation of characteristic design variables of a jacket substructure. In order to resolve the mixedinteger problem formulation, multiple subproblems with fixedinteger design variables are solved. The results show that threelegged jackets may be preferable to fourlegged ones under the boundaries of this study. In addition, it is shown that massdependent cost functions can be easily improved by just considering the number of jacket legs to yield more reliable results.
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 (Mone et al., 2017). Cost breakdowns show ratios of about 20 % (such as The Crown Estate, 2012; BVGassociates, 2013) 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 costsaving potential with low economic effort. Technologically, it is expected that the jacket will supersede the monopile when reaching the imminent turbine generation or wind farm locations with intermediate water depths from about 40 to 60 m (see, for instance, Seidel, 2007; Damiani et al., 2016). According to current studies, there is an increasing market share of jackets (Smith et al., 2015). 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 (van Kuik et al., 2016).
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.^{1} 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.^{2} Therefore, numerical limitations are a great issue in stateoftheart jacket optimization approaches. In the literature, different approaches were presented to address this issue. Schafhirt et al. (2014) proposed an optimization scheme based on a metaheuristic genetic algorithm to guarantee global convergence. To increase the numerical efficiency, a reanalysis technique was applied. Later, an improved approach was illustrated (Schafhirt et al., 2016), where the load calculation was decoupled from the actual tube dimensioning procedure and a simplified fatigue load set (Zwick and Muskulus, 2016) was applied. Similar approaches by Chew et al. (2015, 2016) and Oest et al. (2016) applied sequential quadratic or linear programming methods, respectively, with analytically derived gradients. Other optimization approaches using metaheuristic algorithms were reported by AlHamaydeh et al. (2017) and Kaveh and Sabeti (2018) but without comprehensive load assumptions. The problem of discrete design variables was addressed by Stolpe and Sandal (2018). Oest et al. (2018) 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 (Popko et al., 2014) for the National Renewable Energy Laboratory (NREL) 5 MW reference turbine (Jonkman et al., 2009),^{3} 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 Damiani (2016) and Häfele and Rolfes (2016), 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 (Häfele et al., 2016). 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 Sandal et al. (2018), 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 Häfele et al. (2018a), 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. 2 describes the technical and mathematical problem statements. Both the objective and the constraints are presented and explained in Sect. 3. The optimization approach and methods to solve the problem are discussed in Sect. 4. Section 5 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. 6) and conclusions (Sect. 7).
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 designdriving constraints. Serviceability limit state, i.e., eigenfrequency constraints, is not regarded as designdriving 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 designdriving constraints of jackets are fatigue and extreme loads. In other words, a set of design variables for a parameterizable structure that minimizes its costs, C_{total}, is desirable, while fatigue and ultimate limit state constraints are satisfied; i.e., the maximal normalized tubular joint fatigue damage (among all tubular joints), h_{FLS}, is less than or equal to 1,^{4} and the extreme load utilization ratio (among all tubes), h_{ULS}, 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:
In this equation, the cost value is logarithmized to obviate numerical issues. The constraints, h_{1}(x) and h_{2}(x), are formulated so as to match the requirements of mathematical problem statements; thus
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:
where x is the array or vector of design variables, x_{lb} and x_{ub} are the lower and upper boundaries, respectively, f(x) is the objective function, covering the costs related only to the substructure, and h_{1}(x) and h_{2}(x) are nonlinear constraints representing structural code checks for fatigue and ultimate limit state that are required to be satisfied for every design.
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 (Häfele et al., 2018a).
3.1 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, N_{L}, number of bays, N_{X} (both integer variables), foot radius, R_{foot}, headtofoot radius ratio, ξ, jacket length, L, elevation of the transition piece over mean sea level, L_{MSL}, lowermost segment height, L_{OSG}, uppermost segment height, L_{TP}, the ratio of two consecutive bay heights, q, and a boolean flag, x_{MB}, 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 (N_{L}=4), four bays (N_{X}=4), and mud braces x_{MB}=true is shown in Fig. 1. The tube sizing parameters are the leg diameter, D_{L}, 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 bracetoleg diameter ratios, and τ_{b} and τ_{t} are the bracetoleg 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, ρ.
To decrease the dimension of the problem, height measures related to the location of the wind farm (L, L_{MSL}, L_{OSG}, L_{TP}) and the material parameters (E, G, ρ) are fixed. In addition, it is supposed that each design has mud braces (x_{MB}=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:
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.
3.2 Cost function (objective)
The total capital expenses, C_{total}, comprise several terms, C_{j}, expressed as the sum of socalled factors, c_{j}, weighted by unit costs, a_{j}:^{5}
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 massdependent 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 costdriving factor. To consider known, important impacts on jacket capital expenses, seven factors are incorporated, namely the following:

expenses for material, C_{1}, depending on the mass, c_{1}:
$$\begin{array}{ll}{\displaystyle}{c}_{\mathrm{1}}\left(\mathit{x}\right)=& {\displaystyle}\mathrm{2}\mathit{\rho}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}^{\mathrm{2}}\sum _{i=\mathrm{1}}^{{N}_{X}}\left(\left({\displaystyle \frac{{\mathit{\beta}}_{i}{\mathit{\tau}}_{i}}{\mathrm{2}{\mathit{\gamma}}_{i}}}+{\displaystyle \frac{{\mathit{\tau}}_{i}^{\mathrm{2}}}{\mathrm{4}{\mathit{\gamma}}_{i}^{\mathrm{2}}}}\right)\right.\\ {\displaystyle}& {\displaystyle}\left.\sqrt{{\displaystyle \frac{{L}_{i}^{\mathrm{2}}}{{\mathrm{cos}}^{\mathrm{2}}\left({\mathrm{\Phi}}_{\mathrm{p}}\right)}}+{\left({R}_{i}+{R}_{i+\mathrm{1}}\right)}^{\mathrm{2}}{\mathrm{sin}}^{\mathrm{2}}\left({\displaystyle \frac{\mathit{\vartheta}}{\mathrm{2}}}\right)}\right)\\ {\displaystyle}& {\displaystyle}+{x}_{\mathrm{MB}}\mathit{\rho}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}^{\mathrm{2}}\left({\displaystyle \frac{{\mathit{\beta}}_{\mathrm{b}}{\mathit{\tau}}_{\mathrm{b}}}{\mathrm{2}{\mathit{\gamma}}_{\mathrm{b}}}}+{\displaystyle \frac{{\mathit{\tau}}_{\mathrm{b}}^{\mathrm{2}}}{\mathrm{4}{\mathit{\gamma}}_{\mathrm{b}}^{\mathrm{2}}}}\right)\mathrm{2}{R}_{\mathrm{1}}\mathrm{sin}\left({\displaystyle \frac{\mathit{\vartheta}}{\mathrm{2}}}\right)\\ {\displaystyle}& {\displaystyle}+\mathit{\rho}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}^{\mathrm{2}}\sum _{i=\mathrm{1}}^{{N}_{X}}\left(\left({\displaystyle \frac{\mathrm{1}}{\mathrm{2}{\mathit{\gamma}}_{i}}}+{\displaystyle \frac{\mathrm{1}}{\mathrm{4}{\mathit{\gamma}}_{i}^{\mathrm{2}}}}\right){\displaystyle \frac{{L}_{m,i}}{\mathrm{cos}\left({\mathrm{\Phi}}_{s}\right)}}\right.\\ {\displaystyle}& {\displaystyle}\left.+\left({\displaystyle \frac{\mathrm{1}}{\mathrm{2}{\mathit{\gamma}}_{i+\mathrm{1}}}}+{\displaystyle \frac{\mathrm{1}}{\mathrm{4}{\mathit{\gamma}}_{i+\mathrm{1}}^{\mathrm{2}}}}\right){\displaystyle \frac{\left({L}_{i}{L}_{m,i}\right)}{\mathrm{cos}\left({\mathrm{\Phi}}_{s}\right)}}\right)\\ {\displaystyle}& {\displaystyle}+\mathit{\rho}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}^{\mathrm{2}}\left({\displaystyle \frac{\mathrm{1}}{\mathrm{2}{\mathit{\gamma}}_{\mathrm{b}}}}+{\displaystyle \frac{\mathrm{1}}{\mathrm{4}{\mathit{\gamma}}_{\mathrm{b}}^{\mathrm{2}}}}\right){\displaystyle \frac{{L}_{\mathrm{OSG}}}{\mathrm{cos}\left({\mathrm{\Phi}}_{s}\right)}}\\ \text{(6)}& {\displaystyle}& {\displaystyle}+\mathit{\rho}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}^{\mathrm{2}}\left({\displaystyle \frac{\mathrm{1}}{\mathrm{2}{\mathit{\gamma}}_{\mathrm{t}}}}+{\displaystyle \frac{\mathrm{1}}{\mathrm{4}{\mathit{\gamma}}_{\mathrm{t}}^{\mathrm{2}}}}\right){\displaystyle \frac{{L}_{\mathrm{TP}}}{\mathrm{cos}\left({\mathrm{\Phi}}_{s}\right)}}\text{;}\end{array}$$ 
expenses for fabrication, C_{2}, depending on the entire volume of welds, c_{2}:
$$\begin{array}{ll}{\displaystyle}{c}_{\mathrm{2}}\left(\mathit{x}\right)=& {\displaystyle}\mathrm{2}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}\sum _{i=\mathrm{1}}^{{N}_{X}}\left({\mathit{\beta}}_{i}\left({\displaystyle \frac{{D}_{\mathrm{L}}^{\mathrm{2}}{\mathit{\tau}}_{i}^{\mathrm{2}}}{\mathrm{8}{\mathit{\gamma}}_{i}^{\mathrm{2}}}}+{\displaystyle \frac{{t}_{\mathrm{0}}{D}_{\mathrm{L}}{\mathit{\tau}}_{i}}{\mathrm{2}\sqrt{\mathrm{2}}{\mathit{\gamma}}_{i}}}\right)\right.\\ {\displaystyle}& {\displaystyle}\left(\sqrt{{\displaystyle \frac{\mathrm{1}}{\mathrm{2}{\mathrm{sin}}^{\mathrm{2}}\left({\mathit{\psi}}_{\mathrm{1},i}\right)}}+{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}}+\sqrt{{\displaystyle \frac{\mathrm{1}}{\mathrm{2}{\mathrm{sin}}^{\mathrm{2}}\left({\mathit{\psi}}_{\mathrm{2},i}\right)}}+{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}}\right.\\ {\displaystyle}& {\displaystyle}\left.\left.+\sqrt{{\displaystyle \frac{\mathrm{1}}{\mathrm{2}{\mathrm{sin}}^{\mathrm{2}}\left({\mathit{\psi}}_{\mathrm{3},i}\right)}}+{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}}\right)\right)\\ {\displaystyle}& {\displaystyle}+\mathrm{2}{x}_{\mathrm{MB}}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}{\mathit{\beta}}_{\mathrm{b}}\left({\displaystyle \frac{{D}_{\mathrm{L}}^{\mathrm{2}}{\mathit{\tau}}_{\mathrm{b}}^{\mathrm{2}}}{\mathrm{8}{\mathit{\gamma}}_{\mathrm{b}}^{\mathrm{2}}}}+{\displaystyle \frac{{t}_{\mathrm{0}}{D}_{\mathrm{L}}{\mathit{\tau}}_{\mathrm{b}}}{\mathrm{2}\sqrt{\mathrm{2}}{\mathit{\gamma}}_{\mathrm{b}}}}\right)\\ {\displaystyle}& {\displaystyle}+{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}\sum _{i=\mathrm{1}}^{{N}_{X}}\left({\displaystyle \frac{{D}_{\mathrm{L}}^{\mathrm{2}}min\left(\frac{\mathrm{1}}{{\mathit{\gamma}}_{i}^{\mathrm{2}}},\frac{\mathrm{1}}{{\mathit{\gamma}}_{i+\mathrm{1}}^{\mathrm{2}}}\right)}{\mathrm{8}}}\right.\\ \text{(7)}& {\displaystyle}& {\displaystyle}\left.+{\displaystyle \frac{{D}_{\mathrm{L}}{t}_{\mathrm{0}}min\left(\frac{\mathrm{1}}{{\mathit{\gamma}}_{i}},\frac{\mathrm{1}}{{\mathit{\gamma}}_{i+\mathrm{1}}}\right)}{\mathrm{2}\sqrt{\mathrm{2}}}}\right)\text{;}\end{array}$$ 
coating costs, C_{3}, depending on the outer surface area of all tubes, c_{3}:
$$\begin{array}{ll}{\displaystyle}{c}_{\mathrm{3}}\left(\mathit{x}\right)=& {\displaystyle}\mathrm{2}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}\sum _{i=\mathrm{1}}^{{N}_{X}}\\ {\displaystyle}& {\displaystyle}\left({\mathit{\beta}}_{i}\sqrt{{\displaystyle \frac{{L}_{i}^{\mathrm{2}}}{{\mathrm{cos}}^{\mathrm{2}}\left({\mathrm{\Phi}}_{\mathrm{p}}\right)}}+{\left({R}_{i}+{R}_{i+\mathrm{1}}\right)}^{\mathrm{2}}{\mathrm{sin}}^{\mathrm{2}}\left({\displaystyle \frac{\mathit{\vartheta}}{\mathrm{2}}}\right)}\right)\\ {\displaystyle}& {\displaystyle}+{x}_{\mathrm{MB}}{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}{\mathit{\beta}}_{\mathrm{b}}\left(\mathrm{2}{R}_{\mathrm{1}}\mathrm{sin}\left({\displaystyle \frac{\mathit{\vartheta}}{\mathrm{2}}}\right)\right)\\ \text{(8)}& {\displaystyle}& {\displaystyle}+{N}_{\mathrm{L}}\mathit{\pi}{D}_{\mathrm{L}}{\displaystyle \frac{L}{\mathrm{cos}\left({\mathrm{\Phi}}_{s}\right)}}\text{;}\end{array}$$ 
costs for the transition piece, C_{4}, proportional to the product of head radius and number of jacket legs, c_{4}:
$$\begin{array}{}\text{(9)}& {c}_{\mathrm{4}}\left(\mathit{x}\right)={N}_{\mathrm{L}}{R}_{\mathrm{foot}}\mathit{\xi}\text{;}\end{array}$$ 
expenses for transport, C_{5}, expressed by the massdependent factor, c_{5}:
$$\begin{array}{}\text{(10)}& {c}_{\mathrm{5}}\left(\mathit{x}\right)={c}_{\mathrm{1}}\left(\mathit{x}\right)\text{;}\end{array}$$ 
and installation costs, C_{6}, modeled by a factor only depending on the number of jacket legs, c_{6}:
$$\begin{array}{}\text{(11)}& {c}_{\mathrm{6}}\left(\mathit{x}\right)={N}_{\mathrm{L}}\text{.}\phantom{\rule{0.33em}{0ex}}\end{array}$$
Fixed expenses, C_{7}, are not dependent on any jacket parameter at all. Therefore, the factor, c_{7}, simply takes
In these equations, ϑ is the angle enclosed by two jacket legs:
Bay heights, L_{i}, intermediate bay heights, L_{m,i}, radii, R_{i}, and intermediate radii, R_{m,i}, are calculated by the following equations:
with the spatial batter angle, Φ_{s}:
The interconnecting tube angles, ψ_{1,i}, ψ_{2,i}, and ψ_{3,i}, are
with the planar batter angle, Φ_{p}:
γ_{i}, β_{i}, and τ_{i} represent the ratios of leg radiustothickness, bracetoleg diameter, and bracetoleg 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 massproportional modeling of material costs, C_{1}, is straightforward. Fabrication costs, C_{2}, 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, C_{3}, 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 (stellartype) transition piece, C_{4}, 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, C_{5}, is very difficult. In this work, a massdependent 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 massdependent or dependent on factors partially related to mass (like the space on the deck of the barge covered by the jacket). Installation costs, C_{6}, 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, C_{7}, 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.
3.3 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 RPC203 (DNV GL AS, 2016) for fatigue and NORSOK N004 (NORSOK, 2004) for ultimate limit state checks are adopted. Both are widely accepted for practical applications and were used to design the UpWind (Vemula et al., 2010) and INNWIND.EU (von Borstel, 2013) 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 simulationbased 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 (Häfele et al., 2018a) 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.^{6} 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),^{7} 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.^{8} The specific surrogate models for the considered test problems were derived in a previous work (Häfele et al., 2018a), which revealed that a Matérn 5∕2 kernel function is wellsuited for the present application.
3.3.1 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 614003 (International Electrotechnical Commision, 2009). 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 (Häfele et al., 2017, 2018b). 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 Xjoint positions, the standard deviation of fatigue damage increases by a factor of approximately 4 in the case of a 16fold 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 Efthymiou (1988) to consider the amplification of stresses due to the geometry of tubular joints, rain flow cycle counting, and a lifetime prediction by SN curves and linear damage accumulation.^{9} The output value h_{FLS} 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.
3.3.2 Ultimate limit state
The standard IEC 614003 (International Electrotechnical Commision, 2009) 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 Agarwal and Manuel, 2010), 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 614003) 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 h_{ULS} 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 (Oest et al., 2016).
The optimization problem incorporates a mixedinteger formulation (due to discrete numbers of legs and bays of the jacket). In order to address this issue, the mixedinteger 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 nonconvex; 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. Metaheuristic 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 interiorpoint (IP) methods (Nocedal and Wright, 2006). 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 largescale or sparse problems. In this work, both approaches are used to solve the jacket optimization problem.^{10} They are outlined briefly in the following.
4.1 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 Biggs (1975), Han (1977), and Powell (1978a, b). In the first step, the Hessian of the socalled 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 (Fletcher, 1987). 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 activeset method described by Gill et al. (1981) is deployed for this task. The procedure is repeated until convergence is reached.
4.2 Interiorpoint 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 Waltz et al. (2006).
In this section, the proposed approach is applied to find and compare optimal jacket designs for the NREL 5 MW reference turbine (Jonkman et al., 2009). The environmental conditions are adopted from measurements recorded at the research platform FINO3 in the German North Sea.
5.1 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 highpower 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.
5.2 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 (Fischer et al., 2010), 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 (Hübler et al., 2017). 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 (Häfele et al., 2017, 2018b). 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 deepwater site was considered. Second, the soil properties of the Offshore Code Comparison Collaboration (OC3) (Jonkman and Musial, 2010) 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.
5.3 Boundaries of design variables and other parameters
The boundaries are chosen conservatively by means of quite narrow design variable ranges (see Table 1), 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 (Popko et al., 2014). Only three or fourlegged 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 DNV GL AS (2016), an SN curve with an endurance stress limit of 52.63×10^{6} N m^{−2} at ×10^{7} 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 Häfele et al. (2018a) and set to a_{1}=1.0 kg^{−1} (material), ${a}_{\mathrm{2}}=\mathrm{4.0}\times {\mathrm{10}}^{\mathrm{6}}$ m^{−3} (fabrication), ${a}_{\mathrm{3}}=\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{2}}$ m^{−2} (coating), ${a}_{\mathrm{4}}=\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{4}}$ m^{−1} (transition piece), a_{5}=2.0 kg^{−1} (transport), ${a}_{\mathrm{6}}=\mathrm{2.0}\times {\mathrm{10}}^{\mathrm{5}}$ (installation), and ${a}_{\mathrm{7}}=\mathrm{1.0}\times {\mathrm{10}}^{\mathrm{5}}$ (fixed). With these values, the cost function returns a dimensionless value, also interpretable as capital expenses in EUR.
5.4 Results and discussion
To resolve the mixedinteger formulation of the optimization problem into continuous problems, six subproblems with three legs and three bays (N_{L}=3, N_{X}=3), three legs and four bays (N_{L}=3, N_{X}=4), three legs and five bays (N_{L}=3, N_{X}=5), four legs and three bays (N_{L}=4, N_{X}=3), four legs and four bays (N_{L}=4, N_{X}=4), and four legs and five bays (N_{L}=4, N_{X}=5) were solved using the SQP and IP methods. Therefore, multiple solutions are discussed and compared in the following. The optimization problem is nonconvex; 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 fixedinteger 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.^{11} Gradients were computed by finite differences. The optimization terminated, when the firstorder optimality and feasibility measures were both less than $\mathrm{1}\times {\mathrm{10}}^{\mathrm{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. 2, where the OC4 jacket with varying numbers of legs and bays was assumed as the starting point. This structure has a foot radius, R_{foot}, of 8.79 m, a headtofoot radius ratio, ξ, of 0.67, and a ratio of two consecutive bay heights, q, of 0.8. Moreover, it has a leg diameter, D_{L}, of 1.2 m, and entirely constant tube dimensions from bottom to top, i.e., leg radiustothickness ratios, γ_{b} and γ_{t}, of 15, bracetoleg diameter ratios, β_{b} and β_{t}, of 0.5, and bracetoleg 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 threelegged designs was higher at the beginning of the optimization process but converges stably. For the same reason, the fourlegged 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.
The optimal solutions obtained by the sequential quadratic programming method are illustrated in Table 2.^{12} Additionally, the topologies of all optimal solutions are shown in Fig. 3. 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 threelegged jacket with three bays (N_{L}=3, N_{X}=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, R_{foot}, are at the upper boundaries in the case of the threelegged structures, while the headtofoot 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 fourlegged 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, D_{L}, are relatively small, in the case of the fourlegged 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 threelegged 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. 4, each subplot shows the variation in the total costs, C_{total}, and the cost function terms C_{1} (proportional to C_{5}), C_{2}, C_{3}, and C_{4} due to a 1 % oneatatime variation in each continuous design variable in three different phases of the optimization process (initial, intermediate, and final phase). The terms C_{6} and C_{7} are not impacted by any continuous design variable and therefore not considered. For instance, a 1 % increase in the foot radius, R_{foot}, causes increasing material costs of ΔC_{1}=0.14 %, evaluated for the initial design, but increasing material costs of ΔC_{1}=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, N_{L}, the foot radius, R_{foot}, and the headtofoot radius ratio, ξ). In general, Fig. 4 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, D_{L}, and leg radiustothickness ratios, γ_{b} and γ_{t}, are designdriving (together with the number of legs, N_{L}) 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 C_{4}, which is only impacted by topological design variables, more precisely the foot radius, R_{foot}, and the headtofoot radius ratio, ξ. As a large foot radius, R_{foot}, is needed to establish structural resistance, this cost term penalizes large headtofoot radius ratios, ξ. For this reason, this value is at the lower boundary for all design solutions.
Regarding the costs of the jackets, the best solution with three legs and three bays is related to capital expenses of 10^{6.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 threelegged designs may be favorable in terms of costs (Chew et al., 2014) and threelegged 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. 5, which shows the cost contributions of all six structures and where the actual cost savings come from. The lightest structure is the fourlegged jacket with three bays, while the threelegged 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 fourlegged one with three bays and not the threelegged jacket with three bays, which has the least number of joints. The threelegged 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 threelegged jackets, except for the one with five bays (10^{6.493}= 3 112 000), which is more expensive than the fourlegged one with three bays (10^{6.487}= 3 069 000). The most expensive jackets are the fourlegged ones with four (10^{6.500}= 3 162 000) and five (10^{6.514}= 3 266 000) bays, where the latter is about 15 % more expensive than the best solution among the six subsolutions. A reasonable option may also be the jacket with three legs and four bays, which features a total cost value of 10^{6.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.
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 E52687W v3 central processing unit and 64 GB random access memory. Compared to simulationbased 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 $\mathrm{200}\times \mathrm{128}=$ 25 600 time domain simulations in the fatigue and $\mathrm{200}\times \mathrm{10}=\mathrm{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 highperformance clusters, while the closedloop 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 massdependent one, which can be regarded as state of the art in jacket optimization. For this purpose, all unit costs except a_{1} 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 3. Under these assumptions, the fourlegged jackets are better (in terms of minimal mass) than the threelegged ones. Interestingly, similar design variables are obtained when comparing these values to the ones obtained by the more comprehensive cost model in Table 2, particularly in the case of the threelegged jackets. The resulting objective function values are, in comparison, similar to the material costs in Fig. 5. In other words, a pure massdependent 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 C_{1}…C_{5} depend in some way on the tube dimensions and the topology does not impact the costs to a great extent, as seen in Fig. 4. Indeed, the largest proportion of costs is purely massdependent, as the factors c_{1} and c_{5} are the mass of the structure. Therefore, the proposed cost model can lead to more accurate results, but a massdependent approach would be sufficient to draw the same conclusions.
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 reliabilitybased 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 turbinedependent. 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 designdriving leg diameters are relatively small. This indicates that topological design variables with minor impact on costs are useful factors to establish the structural resistance.
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 20parameter jacket model, of which 12 parameters are design variables. The second, important point leads to the utilization of a more complex (compared to massdependent) 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, gradientbased 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 stateoftheart gradientbased optimization, but activeset SQP and IP methods both converge efficiently and accurately for the given problem.
From the applicationoriented point of view, it can be stated that threelegged 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 loadbearing capacity gained by the extra leg of a fourlegged 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 massdependent 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).
DLC  Design load case 
IP  Interiorpoint method 
SQP  Sequential quadratic programming method 
Φ_{p}  Planar (twodimensional) batter angle 
Φ_{s}  Spatial (threedimensional) batter angle 
β_{b}  Bracetoleg diameter ratio at bottom (jacket model parameter) 
β_{i}  Bracetoleg diameter ratio in the ith bay 
β_{t}  Bracetoleg diameter ratio at top (jacket model parameter) 
γ_{b}  Leg radiustothickness ratio at bottom (jacket model parameter) 
γ_{i}  Leg radiustothickness ratio in the ith bay 
γ_{t}  Leg radiustothickness ratio at top (jacket model parameter) 
ξ  Headtofoot radius ratio (jacket model parameter) 
ρ  Material density (jacket model parameter) 
ϑ  Angle enclosed by two jacket legs 
τ_{b}  Bracetoleg thickness ratio at bottom (jacket model parameter) 
τ_{i}  Bracetoleg thickness ratio in the ith bay 
τ_{t}  Bracetoleg thickness ratio at top (jacket model parameter) 
ψ_{1,i}  Lower bracetoleg connection angle in the ith bay 
ψ_{2,i}  Upper bracetoleg connection angle in the ith bay 
ψ_{3,i}  Bracetobrace connection angle in the ith bay 
C_{j}  Expenses related to jth cost factor 
C_{total}  Total capital expenses 
D_{L}  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) 
L_{MSL}  Transition piece elevation over mean sea level (jacket model parameter) 
L_{OSG}  Lowest leg segment height (jacket model parameter) 
L_{TP}  Transition piece segment height (jacket model parameter) 
L_{i}  ith jacket bay height 
L_{m,i}  Distance between the lower layer of K joints and the layer of X joints of the ith bay 
N_{L}  Number of legs (jacket model parameter) 
N_{X}  Number of bays (jacket model parameter) 
R_{Foot}  Foot radius (jacket model parameter) 
R_{i}  ith jacket bay radius at lower K joint layer 
R_{m,i}  Radius of the ith X joint layer 
a_{j}  jth unit cost 
c_{j}  jth cost factor 
f  Objective function value 
h_{1}  First inequality constraint value 
h_{2}  Second inequality constraint value 
h_{FLS}  Maximal normalized tubular joint fatigue damage 
h_{ULS}  Maximal extreme load utilization ratio 
q  Ratio of two consecutive bay heights (jacket model parameter) 
x  Array of design variables 
x_{lb}  Array of lower boundaries 
x_{MB}  Mud brace flag (jacket model parameter) 
x_{ub}  Array of upper boundaries 
The authors declare that they have no conflict of interest.
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 openaccess
fund of Leibniz Universität Hannover.
Edited by: Lars Pilgaard Mikkelsen
Reviewed by: Lars Einar S. Stieng and one anonymous referee
Agarwal, P. and Manuel, L.: Load Extrapolation Methods for Offshore Wind Turbines, in: Offshore Technology Conference, Houston, Texas, USA, https://doi.org/10.4043/21001MS, 2010. a
AlHamaydeh, M., Barakat, S., and Nasif, O.: Optimization of Support Structures for Offshore Wind Turbines Using Genetic Algorithm with DomainTrimming, Math. Probl. Eng., 2017, 1–14, https://doi.org/10.1155/2017/5978375, 2017. a
Biggs, M. C.: Constrained minimization using recursive quadratic programming, in: Towards global optimization, NorthHolland, Cagliari, Italy, 341–349, 1975. a
BVGassociates: Offshore wind: Industry's journey to GBP 100/MWh – Cost breakdown and technology transition from 2013 to 2020, Tech. rep., 2013. a
Chew, K. H., Tai, K., Ng, E. Y. K., Muskulus, M., and Zwick, D.: Offshore Wind Turbine Jacket Substructure: A Comparison Study Between FourLegged and ThreeLegged Designs, Journal of Ocean and Wind Energy, 1, 74–81, 2014. a
Chew, K.H., Tai, K., Ng, E., and Muskulus, M.: Optimization of Offshore Wind Turbine Support Structures Using an Analytical Gradientbased Method, Energy Proced., 80, 100–107, https://doi.org/10.1016/j.egypro.2015.11.412, 2015. a
Chew, K.H., Tai, K., Ng, E., and Muskulus, M.: Analytical gradientbased optimization of offshore wind turbine substructures under fatigue and extreme loads, Mar. Struct., 47, 23–41, https://doi.org/10.1016/j.marstruc.2016.03.002, 2016. a
Damiani, R.: JacketSE: An Offshore Wind Turbine Jacket Sizing Tool, Technical Report NREL/TP500065417, National Renewable Energy Laboratory, Golden, CO, USA, 2016. a
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, https://doi.org/10.1088/17426596/753/9/092003, 2016. a
DNV GL AS: Fatigue design of offshore steel structures, Recommended Practice DNVGLRPC203, 2016. a, b, c
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. a
Fischer, T., de Vries, W., and Schmidt, B.: Upwind Design Basis, Tech. rep., 2010. a
Fletcher, R.: Practical methods of optimization, Wiley, Chichester, UK, New York, USA, 2nd edn., 1987. a
Gill, P. E., Murray, W., and Wright, M. H.: Practical optimization, Academic Press, London, UK, New York, USA, 1981. a
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 Twentysixth (2016) International Offshore and Polar Engineering Conference, 26 June–2 July 2016, Rhodes, Greece, 156–163, 2016. a
Häfele, J., Gebhardt, C. G., and Rolfes, R.: Innovative design of a hybridtype jacket for 10MW turbines, Tech. Rep. INNWIND.EU Deliverable D4.3.5, 2016. a
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 Twentyseventh (2017) International Offshore and Polar Engineering Conference, 25–30 June 2017, San Francisco, CA, 167–173, 2017. a, b
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, https://doi.org/10.5194/wes35532018, 2018a. a, b, c, d, e, f
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, https://doi.org/10.1016/j.renene.2017.10.097, 2018b. a, b
Han, S. P.: A globally convergent method for nonlinear programming, J. Optimiz. Theory. App., 22, 297–309, https://doi.org/10.1007/BF00932858, 1977. a
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, https://doi.org/10.5194/wes24912017, 2017. a
International Electrotechnical Commision: Wind turbines – Part 3: Design requirements for offshore wind turbines, Tech. Rep. IEC614003:2009, 2009. a, b
Jonkman, J. and Musial, W.: Offshore Code Comparison Collaboration (OC3) for IEA Task 23 Offshore Wind Technology and Deployment, Technical Report NREL/TP500048191, National Renewable Energy Laboratory, Golden, CO, USA, 2010. a
Jonkman, J., Butterfield, S., Musial, W., and Scott, G.: Definition of a 5MW Reference Wind Turbine for Offshore System Development, Technical Report NREL/TP50038060, National Renewable Energy Laboratory, Golden, CO, USA, 2009. a, b
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, https://doi.org/10.1002/tal.1494, 2018. a
Mone, C., Hand, M., Bolinger, M., Rand, J., Heimiller, D., and Ho, J.: 2015 Cost of Wind Energy Review, Technical Report NREL/TP6A2066861, National Renewable Energy Laboratory, Golden, CO, USA, 2017. a
Muskulus, M. and Schafhirt, S.: Design Optimization of Wind Turbine Structures – A Review, Journal of Ocean and Wind Energy, 1, 12–22, 2014. a
National Wind Technology Center Information Portal: FAST v8, available at: https://nwtc.nrel.gov/FAST8 (last access: 28 August 2018), 2016. a
Nocedal, J. and Wright, S. J.: Numerical optimization, Springer series in operations research, Springer, New York, USA, 2nd edn., 2006. a
NORSOK: Design of steel structures, Standard N004 Rev. 2, 2004. a
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, https://doi.org/10.1007/s001580161527x, 2016. a, b
Oest, J., Sandal, K., Schafhirt, S., Stieng, L. E. S., and Muskulus, M.: On gradientbased optimization of jacket structures for offshore wind turbines, Wind Energy, 21, 953–967, https://doi.org/10.1002/we.2206, 2018. a
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., VasquezRojas, 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. a, b
Powell, M. J. D.: The Convergence of Variable Metric Methods for Nonlinearly Constrained Optimization Calculations, Nonlinear Programming, 3, 27–63, https://doi.org/10.1016/B9780124686601.500074, 1978a. a
Powell, M. J. D.: A fast algorithm for nonlinearly constrained optimization calculations, in: Numerical Analysis, Dundee, Scotland, 630, 144–157, https://doi.org/10.1007/BFb0067703, 1978b. a
Rasmussen, C. E. and Williams, C. K. I.: Gaussian processes for machine learning, MIT Press, Cambridge, MA, USA, 3rd print edn., 2008. a
Sandal, K., Latini, C., Zania, V., and Stolpe, M.: Integrated optimal design of jackets and foundations, Mar. Struct., 61, 398–418, https://doi.org/10.1016/j.marstruc.2018.06.012, 2018. a
Schafhirt, S., Zwick, D., and Muskulus, M.: Reanalysis of Jacket Support Structure for ComputerAided Optimization of Offshore Wind Turbines with a Genetic Algorithm, Journal of Ocean and Wind Energy, 1, 209–216, 2014. a
Schafhirt, S., Zwick, D., and Muskulus, M.: Twostage local optimization of lattice type support structures for offshore wind turbines, Ocean Eng., 117, 163–173, https://doi.org/10.1016/j.oceaneng.2016.03.035, 2016. a
Seidel, M.: Jacket substructures for the REpower 5M wind turbine, in: Conference Proceedings European Offshore Wind 2007, Berlin, Germany, 2007. a
Smith, A., Stehly, T., and Musial, W.: 20142015 Offshore Wind Technologies Market Report, Technical Report NREL/TP500064283, National Renewable Energy Laboratory, Golden, CO, USA, 2015. a
Stolpe, M. and Sandal, K.: Structural optimization with several discrete design variables per part by outer approximation, Struct. Multidiscip. O., 57, 2061–2073, https://doi.org/10.1007/s0015801819413, 2018. a
The Crown Estate: Offshore Wind Cost Reduction – Pathways Study, Tech. rep., 2012. a
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.: Longterm research challenges in wind energy – a research agenda by the European Academy of Wind Energy, Wind Energ. Sci., 1, 1–39, https://doi.org/10.5194/wes112016, 2016. a
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. a, b
von Borstel, T.: Design Report Reference Jacket, Tech. Rep. INNWIND.EU deliverable report D4.3.1, 2013. a
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, https://doi.org/10.1007/s1010700405605, 2006. a
Zwick, D. and Muskulus, M.: Simplified fatigue load assessment in offshore wind turbine structural analysis, Wind Energy, 19, 265–278, https://doi.org/10.1002/we.1831, 2016. a
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, Muskulus and Schafhirt (2014).
During conceptual design phases, the number of load cases is commonly reduced.
It is worth mentioning that the Offshore Code Comparison Collaboration Continuation (OC4) jacket is actually a structurally reduced derivation of the socalled UpWind jacket (Vemula et al., 2010), 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.
All fatigue damage is normalized so that the lifetime fatigue damage corresponds to a value of 1.
FASTv8 (National Wind Technology Center Information Portal, 2016) was used for this study.
All parameters of these jacket samples are given in the publication where the surrogate modeling approach was reported (Häfele et al., 2018a).
For the background theory of GPR, the reader is referred to Rasmussen and Williams (2008), which is the standard reference in this field.
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 RPC203 (DNV GL AS, 2016).
The function fmincon in MATLAB R2017b was used for this study.
The values shown in the following include all cost terms. The exclusion is only performed during optimization.
As the accuracy of the SQP and IP methods is similar here, only results obtained by the SQP method are shown in the following.