Abstract

The powerful genetic algorithm optimization technique is augmented with an innovative “domain-trimming” modification. The resulting adaptive, high-performance technique is called Genetic Algorithm with Domain-Trimming (GADT). As a proof of concept, the GADT is applied to a widely used benchmark problem. The 10-dimensional truss optimization benchmark problem has well documented global and local minima. The GADT is shown to outperform several published solutions. Subsequently, the GADT is deployed onto three-dimensional structural design optimization for offshore wind turbine supporting structures. The design problem involves complex least-weight topology as well as member size optimizations. The GADT is applied to two popular design alternatives: tripod and quadropod jackets. The two versions of the optimization problem are nonlinearly constrained where the objective function is the material weight of the supporting truss. The considered design variables are the truss members end node coordinates, as well as the cross-sectional areas of the truss members, whereas the constraints are the maximum stresses in members and the maximum displacements of the nodes. These constraints are managed via dynamically modified, nonstationary penalty functions. The structures are subject to gravity, wind, wave, and earthquake loading conditions. The results show that the GADT method is superior in finding best discovered optimal solutions.

1. Introduction and Background

Offshore structures offer presumably the greatest level of sophistication and difficulty in the analysis, design, and construction stages. This is due to multitude of complexities introduced by the offshore locations. The entailed unique conditions such as underwater currents, surface waves and wind, and seabed soil can often become problematic quickly. Such conditions introduce additional factors to the design challenge, such as extreme temperatures and rough weather conditions, marine growth, fatigue, corrosion, marine vessel impact, and interaction with aeroelastic and hydrodynamic loads. These conditions need to be thoughtfully taken into consideration due to the long and short term effects that can compromise the structure. The intricacy of the design and analysis process multiplies when earthquake loads are factored in. However, it is recommended to account for all loads simultaneously except for seismic loads, based on the potential occurrence combinations [1, 2]. Therefore, seismic effects are not typically considered to occur at the same time as other loads [3]. The considered combinations of different loads and their influence on the design of offshore structures are further elaborated on in a later section.

The design solution selected from a range of viable options often imposes additional restraints to the structural design as it shapes the structure’s response to the applied loads. The size, type, and configuration of the members supporting offshore structure, for instance, may produce a significant influence on marine growth size, which poses a consequent impact on hydrodynamic drag forces. The resulting effects of nonlinearities associated with fluid-structure interaction may be too significant for a thoughtless dismissal in the design phase [4]. Furthermore, the interaction between the soil, foundation, and the structure is crucial phenomenon to be studied for the development of enhanced structural system designs supporting offshore wind towers [510].

A common simplified approach to the complex process of offshore structural design is the use of uncoupled load effects [1113]. The justification of the above simplified method in fatigue analysis is attributable to the following: (a) material and geometric nonlinearities are not significant structural fatigue responses, (b) wind and wind-induced currents and waves can be assumed as stationary and independent processes, and (c) soil stress levels are moderate to presume linear behavior under wind and sea wave loading [14]. Employing such approach allows for the advantageous implementation of simple soil-structure interaction behavior.

Devising optimal designs that satisfy multiple performance criteria, such as minimizing cost and maximizing efficiency, is arguably one of the most influential factors in modern structural design. Member-level (local) optimization techniques are routinely used to identify design quantities for acceptable structural performance. Within this framework, several optimization techniques have been utilized in the past [15]. In the presence of multiple optima and nonsmooth constraints in the design variable space, it is difficult to obtain a set of optimum values using member-level optimization formulations. This shortcoming inspired researchers to explore using the relatively new and innovative evolutionary-based optimization techniques. For example, these methods have been used recently to address structural engineering optimization problems: Shuffled Complex Evolution Optimizer (SCEO), Ant Colony Optimization (ACO), the Genetic Algorithms (GA), and Particle Swarm Optimizer (PSO) (e.g., [1624]).

Many structural optimization problems involve problem-specific constraints applicable to the solutions limiting the feasible search space. In these types of problems, it is challenging to adapt traditional optimization techniques to handle constraints. One of the most popular constraint handling methods is by incorporating penalty functions due to the relative simplicity and ease of implementation.

Topology optimization was used by many researchers to generate alternative structural design concepts for benchmark wind turbine blades. In these studies, the focus is alternative structural layout for wind turbine blades with the aim of improving its design, minimizing weight, and, ultimately, wind energy cost reduction (e.g., [2535]).

The materials’ cost for offshore wind turbine supporting structures constitutes a considerable amount of overall cost. Potentially significant cost savings that result from optimized structural designs encourage the incorporation of efficient optimization techniques. However, determining the optimal structural shape and weight of the supporting system is not a trivial task. The complexity is increased due to the mathematical description of loads (aeroelastic [wind], hydrodynamic, and seismic), the many variables describing the geometry, the nonsmooth objective function, and the constraints that have to be satisfied. The above considerations have inspired the implementation of the efficient Genetic Algorithm with Domain-Trimming (GADT) optimization technique.

This study presents the implementation of the GADT to achieve superior optimization. The capabilities of the presented optimization tool are demonstrated on a well-established optimization benchmark problem known for being challenging with known global and multiple local minima. Finally, the GADT tool is used to solve two typical least-weight design problems for offshore wind turbine supporting structures through topology and member size optimization. The algorithm is coded using the MATLAB [36] commercial software package and its GA libraries, while the structural analysis and design for the different topology and loading configurations is performed using the SAP2000 [37] commercial software package. The linkage and automation of the structural analysis and design tasks are initiated from within the MATLAB code and established through the Open Application Programming Interface (OAPI) capabilities of the SAP2000 software.

2. Primary Loads on Offshore Structures

2.1. Hydrodynamic Loads (Sea Waves and Currents)

Hydrodynamics [the kinematics of the water particles] is a very sophisticated branch of applied science that has been studied extensively. The motion of salty ocean water is a result of several causes, such as tidal effects that initiate underwater currents, thermal gradient effects, surface-level winds, and topology of the seabed especially for shallow waters (where the water depth is lesser than half the characteristic wavelength). All of the mentioned which are site dependent and their individual influences and their interactions would potentially vary significantly.

Several wave theories exist which span in sophistication and acceptance. Among the widely used theories are the following: linear or Airy wave theory, Stokes 2nd-order and higher order theories, and stream-function and cnoidal wave theories [44, 45]. The range of applicability of the different wave theories can be found in [4]. The superposition of several Airy wavelets is routinely used to model the irregular sea states. The specifics of the site conditions can be captured through the different wavelength, amplitude, phase, and direction of the superimposed linear waves [4, 46].

When wave loading is applied to an offshore structure, by means of Airy’s wave theory and Morison’s equation, gravity loads are assumed to be simultaneous. The analysis often accounts for the various nonlinearities introduced by the nonlinear hydrodynamic force (drag and lift). Typical representation of wave loads is based on the three parameters: wave period, wave height, and mean water depth (Figure 1(a)). The velocity along the water depth is assigned a nonlinear profile to account for the appropriate experimentally observed conditions, as can be seen in Figure 1(b). Drag and inertial force profiles are also representative of realistic situations. Figure 1(c) schematically depicts the different wave and current actions and effects on an offshore pole with a circular cross-section of diameter ().

In such mathematical procedures, wave and current kinematic fields are often modeled using 5th-order Stokes wave theory. Forces on individual structural elements are then generated using Morison’s equation. These forces account for the inertial effects through the mass coefficient and the hydrodynamic drag via the coefficients [47]. The used values for this research are adopted from the American Petroleum Institute (API) standards [1, 2] and will be introduced in Section 4.1.

Specifics and more details of quantifying such forces can be found elsewhere [4, 4749]. Although the forces induced by ocean waves on platform are dynamic in nature, it is the common practice to design by static approaches. This is deemed a reasonable approximation especially for rigid offshore structures submerged in shallow to intermediate water levels. As the water depth increases and/or offshore structures become more compliant (flexible), the dynamic effects gain significance. In such cases, static loading conditions assumption will no longer be appropriate.

2.2. Current Loads

Waves induce orbital particle motion through the matter in which they travel; such orbits are typically closed but may undergo minor forward drift resulting from wind-surface interaction effects. So, in essence, currents are produced due to wave kinetic energy. If a current coincides with a wave propagation path or direction, the wave length is typically elongated [1, 48].

2.3. Wind Loads

Atmospheric temperature and pressure gradients provide wind with ample kinetic energy. As obstacles are placed in the path of this kinetic energy, some or all of it is converted into potential energy due to path deflection or obstruction and discontinuation. The generated potential energy is manifested through differential pressures on the obstacle body. The resulting pressures are influenced by many factors such as the obstacle shape, orientation in space, and contact area, as well as the wind speed and angle of approach. Wind-induced forces are highly dynamic in nature, but often times, for design purposes, it is deemed appropriate to represent them by equivalent static pressures.

Land-based structures approaches can be adopted for representing wind loads on offshore structures with proper adjustments accounting for open ocean surface roughness (lower category of surface roughness). The lesser roughness produces lower levels of turbulence than would be produced on land; this is translated into a slower rate of variation with height. That is, at any given height, the storm conditions are more severe and produce higher wind speeds in offshore locations as compared to onshore. An illustrative schematic is provided in Figure 2.

2.4. Earthquake Loads

As the ground shakes under a structure during a seismic event, inertial forces are generated within the structure. During that oscillatory vibrational movement, other forces are engaged as well, namely, dissipative forces (damping) and restoring forces. The interaction between all these forces to maintain dynamic equilibrium gives rise to complex structural response. Being an extreme event, seismic activity often involves high structural response demands where inelastic and nonlinear response is to be expected, thus adding to the already complex interactive loading situation in an offshore site. Luckily, significant seismic events are not very frequent and require relatively a long time to build up (store) strain energy. This warrants separate (nonsimultaneous) consideration of earthquake loads from other loads for offshore structures. At offshore sites, tidal waves may or may not require special attention depending on the specific site conditions. Tidal waves travel with great velocities at sea but with small amplitude (wave height) as large bodies of water are displaced during an earthquake. It is only as these waves approach the shoreline that they attain huge amplitudes due to the shallow depths of the seabed. These scenarios are considered on site-specific basis and require special seismic hazard studies prior to any design efforts. The authors are not aware of any specific tidal wave design codes for offshore structures.

When seismic loads are deemed significant to control the offshore structural design either partially or entirely, equivalent static representations are often utilized. This is very much similar to design considerations of on-land sites. If the offshore structure under consideration is irregular or complex in geometry, dynamic (linear or otherwise nonlinear) response needs to be properly modeled. The details of the used design code, design loads, design load combinations, and so forth are specified in Section 4.1 where a numerical example is introduced.

3. Genetic Algorithm with Domain-Trimming (GADT)

3.1. Development of the GADT Optimization Technique

In the master’s thesis of the third author [50], the GADT optimization algorithm was introduced by augmenting the GA with a domain-trimming (DT) to enhance its capacity of discovering global optima. Conventional GA optimization begins with a randomly generated population of solution alternatives (population size) using a uniform probability distribution; a combination of variables comprises each solution of the GA having its own fitness value. When the optimization is performed to find the least weight for a given problem, a function = total mass + penalty must be minimized. Low yielding solution alternatives would have better fitness while they maintain the problem constraints satisfied. A fitness value is allocated to each solution using (1) once is computed for all the solutions in the initial population. Solutions with higher than the population’s are eliminated by an assigned fitness value of zero as they are deemed unfit: where is the design variables vector.

The three basic operations of a GA—reproduction, crossover, and mutation—are the tools improving each population’s fitness across generations/iterations. During reproduction, better fit designs are selected, copied, and placed into a pool allowing each design to mate and reproduce. The selection method used in this study is the roulette wheel due to its simplicity and popularity. One prominent advantage of the GA optimization is its inherent resilience against converging to local optima rather than global optima. GA’s resilience can also be sometimes a potential shortcoming as its “low-resolution” becomes problematic. This is especially prominent when the population size is relatively limited. In such cases, the technique may converge to solutions that are in the vicinity of, but not necessarily at, the global minima. An explanation to such occurrence is that solution variables are first randomly selected from a predetermined domain, which do not change in nature but rather in location from one solution to the other, except for mutation, which occurs at small probability. Such issue is aggravated by poor choices of solution domains, for example, selecting a domain of 1, 1000 while the optimal is at a value of 2. While blending techniques aim at reducing this drawback, they may pose a distinct undesirable effect of refuting the strong traits of the parents.

This study features a new technique developed to tackle GA’s limitation, which involves first reducing the selected domain and then recommencing the GA algorithm to enhance the probability of selecting the optimal solution. To demonstrate, take the previous example of an optimum solution of 2 which is searched for in a domain of 1, 1000. When the same solution is searched for in a “trimmed” domain of 1, 10, the probability of being selected as the optimal solution is increased 200-fold, when the GA is reinitiated and all else maintained constant. Domain-trimming is carried out herein by reducing a percentage (10%) of the initial domain. A continuity of domain reduction and reinitiation is maintained until the global optimum solution is located.

Naturally, the possibility that the global optima may exist within the removed (trimmed) region of the solution space persists. In such cases, the GADT will converge to equivalent solution achieved by the typical GA algorithm, as trimming is queued only by the initial convergence of the GA process. Experimentation has demonstrated that the GADT technique will typically outperform conventional GA as will be shown herein with the offshore support structure examples.

3.2. Constraints Handling

Penalty functions are very useful in enabling the solution of constrained problems as unconstrained. The solutions that violate any number of constraints are penalized, which results in associating high objective function values to nonfeasible solutions. Of the two types of penalty functions, stationary (static) function generally exhibit an inferior performance to that of nonstationary (dynamic) function. A nonstationary penalty function is generally defined as follows:where is the original objective function of the constrained optimization problem and is a dynamically modified penalty, defined aswhere is the performance criteria of an individual member (often taken as the ratio of utilization to satisfaction of pertinent design code, or the demand-to-capacity ratio, etc.). Figure 3 pictorially lays out the flowchart for the general GADT algorithm as well as the specifics of the domain-trimming subalgorithm. The algorithm starts by generating a randomized set of potential population of “parent” solutions and then sets out to order them based on their fitness. For fit solutions, no actions are necessary at this point, while the unfit solutions are penalized such that their fitness function values will cause them to be eliminated from the pool. After convergence checking is performed (including maximum number of solution iterations), the GA operations are performed in order to create offsprings that carry the favorable traits of the parents. The top 10 most fit solutions are ordered, updated, and saved. Domain-trimming (DM) is performed once the top 100 most fit solutions are collected. In this step, the domains are reduced such that the least fit solutions are eliminated from the pool of feasible solutions. The domain-trimming was carefully executed in order to guarantee survival of the fittest solutions upon GA reinitiation, which utilizes this pool of solutions as the starting population for the next optimization cycle, and so on until convergence is achieved or any of the stopping criteria is reached, for example, maximum number of populations is reached.

While and are the penalty exponent and coefficient, respectively, they enable penalizing the optimization objective function value if exceeds unity. Both and will minimize the probability of unfit solutions reappearing in subsequent generations and trimmed solution domains by penalizing them with proportional severity to their excessive values above unity. In this investigation, several combinations were considered (Figure 4) and the used values of the penalty coefficient and exponent are 2 and 2.5, respectively.

3.3. GADT Verification Example: Cantilever 10-Bar Planar Truss Structure

The GADT technique is tested against a classical global optimization problem, the Cantilever 10-Bar Planar Truss Structure optimization problem, before using it to solve offshore problems. This 10-dimensional problem has been investigated by many researchers and has been well-established as a challenging optimization benchmark problem. It has a best discovered global minimum as well as multiple local minima.

Figure 5 shows an illustration of the benchmark structure tested. The length is 360 in (914.4 cm), the modulus of elasticity is 10,000 ksi (68,950 MPa), and the assumed material density is 0.1 lb/in3 (2767.990 kg/m3). The member’s deflection and stress restrictions are ±2.00 in (5.08 cm) and ±25.0 ksi (172.375 MPa), respectively. Cross-sectional areas for all structural members can assume any values between 0.1 in2 and 35 in2 (0.6452 cm2 and 225.806 cm2). All members are designed and optimized independently, which could yield in a unique cross-sectional area for each member (A1 to A10). This truss optimization problem has 10 design variables that are subject to 32 constraints (10 compression stresses, 10 tension stresses, and 12 displacements). The studied scenario is when = 150 kips (667.2 kN) and = 50 kips (222.4 kN).

Table 1 gives the best optimum computed solutions with the corresponding minimum weight for the considered scenario, compared against optimal designs presented by scholarly studies available in the literature. Note that it required nearly 4000 iterations to reach the best discovered solution with an initial population size of 500 solutions prior to any domain reductions.

The optimal solutions found by GADT satisfied all the problem constraints and Table 1 demonstrates how GADT’s performance is superior to preestablished solutions.

Figure 6 shows the convergence history for the GA and the GADT solutions. The GADT is distinctively outperforming the conventional GA, by converging to more optimal solution in lower number of iterations. Notice that GADT convergence plateaus after approximately 800 iterations after which any trimming has minimal effect on the GA, indicating that most or all unfit values have been eliminated.

4. Problem Formulation for Offshore Jacket Structure Optimization

The proposed optimization problem for the supporting wind turbine truss structure is highly sophisticated as it requires the determination of the optimum configuration of geometry, shape, and member sizes simultaneously. The mathematical expression of the problem is as follows:where is the design variable (member cross-sectional area), is the objective function (the structural weight), (, , and ) are the members’ end nodes coordinates, ρ is the material density, is the number of inequality constraints , is the number of the design variables, is the member length, and and are the upper and the lower bounds of the th variable. Joint displacements and truss member stresses based on the specified design code constitute the upper and lower bounds of the constraints given in (4).

Combining the penalty and objective functions produces the fitness function used in GADT:

To limit the spatial shapes of structure to practical and feasible potential solutions, joint coordinates are restricted to remain on the main (vertically inclined) legs but are free to move along them (see Figure 7). Thus, the design parameters and their combinations that would be explored are as follows: the slope of the main legs (), the members’ lengths (), diameters (), and wall thicknesses (WT). For each spatial configuration combined with members’ sizes, the total structural weight is calculated as the objective function to be minimized.

4.1. Numerical Illustrations

The efficiency of the GADT algorithm is demonstrated through two widely used offshore wind turbine support systems, namely, the quatropod and the tripod jacket trusses (Figure 8). The two design problems are subjected to the offshore loading conditions described earlier. A total mass of 3.4 × 105 kg represents the turbine mass assumed to be atop the towers. The supporting truss systems are analyzed for load combinations and designed for conditions specified according to the ASCE7-10 [51] design specifications and the 2010 AISC-LRFD design specifications for steel structures [52].

The design constraints included the traditional permissible member stresses and truss joint displacements prescribed in the adopted AISC-LRFD design standard. All structural elements are modeled as frame elements subjected to gravity as well as wind, wave, and seismic loads. The dead loads are accounted for by the member dimensions and density, while the turbine weight is superimposed at the appropriate location. Typical live loads are treated as movable and temporary, and, as such, a static equivalent uniform 50 psf (2.39 kPa) was used. A typical water depth for the type of jacket was adopted 110 ft (33.5 m). The wave loading application was carried out via a 5th-order Stokes theory representation. The design wave was designated as a 100-yr return event with a height of 26 ft (7.92 m) and the period of 8 sec. The horizontal velocity contours for the waves during the 100-yr design storm conditions are shown in Figure 9. The and values are adopted from the recommendations of the API document [1, 2] to be 0.65 and 1.6, respectively. Since the wave loads are represented in a static fashion, the load direction and incidence angle are varied and an envelope of the different scenarios was used for the analysis and design. Specifically, the wave loading was applied in the following configurations: 0°/180°, 45°/135°, and 90°/270°. The wind loading was combined with the wave loading due to their close coupling and associations. The design wind speed was chosen to be 100 mph (45 m/s) measured at 30 ft (10 m) above the Lowest Astronomical Tide (LAT). The wind speed variation with elevation is typically assumed to follow a power law relationship. Equation (6) is used for wind speed evaluations:where is velocity at height of feet (ft/sec); is velocity at height of 30 feet (ft/sec); is height above LAT (feet); = 0.125; shall not be less than .

Seismic loads were accounted for using the 2006 International Building Code [53] response spectrum constructed based on a 2,500-yr event corresponding to high seismicity. All specifics of the load case definition are in Figure 10. The earthquake loads were also applied in the same directions as the wind and wave loads in order to justify implementing them in load combinations as per the ASCE7-10 [51] design specifications. The ASCE7-10 strength deign load combinations are well-known and not listed herein for brevity.

5. Results and Discussion

The GADT algorithm is applied to the structures mentioned above. The effects of the GADT control parameters (population size, penalty coefficient, and exponent) on the optimum solutions are explored and the combination producing the best results is adopted but omitted for brevity.

5.1. Quatropod (64-Bar) Jacket Structure

The truss members are assigned to one of 14 member groups as the 64-bar 3D truss structure is symmetric about both principal axes ( and ). The 14 groups are as follows: 1 group for the slope of the main legs, 1 group for supports, 1 group for the rod, 1 group for platform beams, 5 groups for bracing members, and 5 groups for legs members. Figure 7 illustrated each of these groups. Therefore, the optimization problem has 33 independent design variables and 32 constraints (6 displacements, 13 tension stresses, and 13 compression stresses). Table 2 gives the best discovered solutions for the 64-bar 3D truss structure.

The GADT algorithm produced an optimum weight of 3,965,364 kgf after completing 500 searches as shown in Figure 11. It is worth mentioning that the optimal solutions found by the GADT meet all constraints and, thus, have no active constraints.

5.2. Tripod (48-Bar) Jacket Structure

Similar to the quatropod jacket, the tripod jacket truss members are assigned to one of the 14 member groups that are illustrated in Figure 7. Again, the optimization problem has 33 independent design variables and 32 constraints. Table 3 gives the best discovered solutions for the 48-bar 3D truss structure.

The GADT algorithm produced an optimum weight of 3,901,710 kgf after completing 1000 searches as shown in Figure 12. Similarly, the optimal solutions found by the GADT meet all constraints, and thus, have no active constraints.

5.3. Computational Efficiency versus Best Discovered Optima

The GADT is observed to outperform the conventional GA technique in terms of discovering optima consistently throughout the three problems presented in this paper. Typically, such improvements to optimal solutions come at a computational price. The domain-trimming and GA reinitiation are the main sources of additional computational cost. However, it can be reported that the GADT is still considered to be a rather computationally efficient technique. In comparison to the conventional GA, the GADT optimization algorithm offers improvements to the discovered optima, which far exceed the minute extra numerical burden. A head-to-head performance comparison on the 10-bar benchmark problem demonstrated the superiority of the GADT over conventional GA. A more thorough and conclusive comparison between the GADT versus conventional GA can be included in a future study.

6. Summary and Conclusions

The genetic algorithm optimization technique is augmented with a novel “domain-trimming” variation. The resulting superior optimization technique is called Genetic Algorithm with Domain-Trimming (GADT). Initially, the GADT performance is demonstrated by application to the widely used 10-dimensional truss benchmark optimization problem. The 10-bar truss has best discovered global and local minima that are published in the literature. The GADT is shown to be superior to several published solutions. Subsequently, the GADT is deployed onto the main focus of this work, which is three-dimensional structural design optimization for offshore wind turbine supporting structures. This design optimization problem is highly complex entailing least-weight topology and member size optimization. The GADT is applied to two popular alternatives for the offshore wind turbine support structures, namely, quadropod and tripod jacket trusses. This design optimization problem, with its two versions, is nonlinearly constrained where the objective function is the material weight of the supporting truss. The selected design variables are the truss members’ cross-sectional areas, as well as the truss members’ end node coordinates (which determine the members’ lengths). The maximum displacements of the nodes and the maximum stresses in members are the applied constraints. These constraints are managed via nonstationary, dynamically modified penalty functions. During analysis and design of the two jacket structures, they are subjected to gravity, wind, wave, and earthquake loading conditions. The results show that the (GADT) method is very efficient in finding the best discovered optimal solutions in both cases. The resulting savings in the optimally designed offshore support structures using the GADT far outweighs the necessary additional computational effort.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.