Multi-domain optimization of cast iron components in wind turbines

The continuously rising demand for renewable energies leads to increased installations of wind turbines with higher power. While the current power-to-weight ratio of up to 20 metric tons of cast iron per megawatt is stagnating, cast iron components of modern wind turbines are facing new challenges in terms of weight, manufacturability, and castability. These challenges can be addressed by systematically using multi-domain optimization approaches to reduce component weight and increase local component utilization. In order to meet the requirements for modern cast iron components, this multi-domain approach must employ methods from casting simulation, micromechanical analysis, topology optimization, and strength assessment. Here, casting simulation is used to determine local microstructure descriptors, which are subsequently used in micromechanical shakedown analysis to estimate the local microstructure-dependent fatigue strength. In parallel to the fatigue strength estimation, topology optimization is performed iteratively in combination with a castability analysis. The component strength is evaluated using a strength assessment approach based on the previously determined local material properties in combination with the topology optimized component. In this study, the overall concept of the proposed multi-domain approach is presented and requirements for the application of such an approach are formulated. The use case of this study is a planet carrier of a wind turbine gearbox manufactured from austempered ductile cast iron ADI-GJS-1050‑6. For this use case, a weight reduction of 17% was achieved while maintaining the required stiffness, such that the microstructure variance along the component was significantly reduced. Furthermore, the potentials and limitations of the presented approach are outlined and discussed in the context of the design of heavy-section castings.


Introduction
The transition towards a carbon-neutral society, the stated aim of the European Union [1], requires ambitious expansion targets for the installed renewable energy production. As wind turbines have become a cornerstone of renewable energy production, the expansion targets require an increase in installed wind turbine power. This can be achieved through an accelerated installation of wind turbines and an increase in wind turbine-specific power. However, modern wind turbines still contain up to 20 tons of cast iron per installed megawatt due to heavy-section castings such as torque arm, main shaft, or planet carrier that are subjected to cyclic loading during operation [2]. Thus, accelerated installation and up-scaling create challenges in manufacturing and logistics due to the component weight. The multidomain lightweight design concepts enable the up-scaling of wind turbines to higher power density using weight-reduced components. The weight reduction in turn leads to an optimization of CO2 emissions over the entire life cycle of new wind turbines. Therefore, unused component potentials of heavy-section castings, such as those given by consideration of local material properties or new engineering materials, need to be exploited for component design process. This requires, firstly, the material properties to be numerically estimated based on limited microstructure information, and secondly, the local material properties to be considered in a component design concept.
For ductile cast iron, as used in wind turbines, microstructure-based modeling has been investigated using different approaches: Based on simplified unit-cell models, as used by Bonora et al. [3] to study damage mechanisms. More complex material models have been incorporated into advanced concepts, such as those used by Andriollo et al. [4,5]. Utilizing the results of advanced material modeling, the behavior of ductile cast iron under cyclic loading has been investigated experimentally and numerically at the microscale [6][7][8]. Thereby, numerical approaches, such as the shakedown theorem, have been developed to determine the fatigue strength of a given microstructure [6,9].
With respect to the consideration of local material properties in design concepts of cast iron components, Olofsson et al. [10,11] have presented first multi-domain approaches. Here, local material properties are considered along a cast iron component based on the casting simulation results. Thus, a strength assessment of the component has been performed using the estimated local material properties [11,12]. Progress has been made in other work in the field of topology optimization of castings. Li et al. took into account the manufacturing conditions in terms of gap dimensions and undercuts [13]. For this purpose, a quantity, virtual temperature, was defined, which allows specifying multidirectional demolding conditions. Approaches to stressbased topology optimization for castings were presented by Wang et al. [14]. Here, a comparison was made between a stress-based and a stiffness-based optimization. Recently, Lubenow et al. [15] used an optimization scheme of topology optimization, casting simulation, and local computeraided optimization to optimize a planet carrier.
However, multi-domain design concepts with numerical fatigue strength estimation have not yet been applied to heavy-section castings made of austempered ductile cast iron, such as the planet carrier. For such components, the local microstructure varies significantly such that a numerical estimation of the local fatigue strength can reduce experimental requirements. In combination with the localized strength assessment a component potential can be opened up.
This study presents an approach to consider a local strength assessment in the design of heavy-section castings made of austempered ductile cast iron for wind turbines using a planet carrier as a demonstrator. Potentials for weight reduction under certain boundary conditions given by operating requirements are determined. The material investigated, an austempered ductile cast iron, is described for this use case. The fatigue strength is estimated using shakedown analysis of the microstructure. Combined with the casting simulation, this allows for a high resolution of local fatigue strength within the casting. Subsequently, topology optimization of the planet carrier is performed and a modified strength assessment is used to determine local utilization rates. Finally, the optimization potential of the planet carrier is estimated based on topology optimization and local strength assessment results.

Material
The material investigated in this study is a heavy-section cast of austempered ductile cast iron grade ADI-GJS-1050-6. Sampling was done using cubes of side lengths 150, 200, and 250 mm (casting modules 2.5, 3.3, and 4.2 cm, respectively) cast by Buchholz & Cie. foundry in a sand mold and heat treated by ADI treatments. The material was subjected to a two-stage heat treatment: austenitizing at 900°C for 7 h and austempering at 300°C for 4 h. Specimens were taken from the center of each cube for material characterization as well as fatigue testing. Exemplary micrographs are shown in Fig. 1 for each casting modulus. The microstructure descriptors analyzed according to standard EN 945-4 using five micrographs of size 1.387 × 1.045 mm 2 at 100 × magnification are given in Table 1. Chemical analysis was performed using a melt sample during casting with the results given in Table 2. Results of uniaxial tensile testing of five specimens of type A with  Table 3. Additionally, results of fatigue testing (fatigue strength at a survival probability of 1%) according to DIN 50100 in a staircase testing procedure with a logarithmic increment of d log = 1.069 [17] are given in Table 3.
Fatigue testing was performed on a Roell Amsler HFP 422. The remaining specimens and runouts were used to determine the finite life regime. The results were evaluated using the software SAFD 5 [18].
Comparing the experimentally determined fatigue strengths, a maximum can be observed at a casting modulus of 3.3 cm. This fits with the microstructure descriptors, as casting modulus 2.5 cm shows a reduced nodularity and casting modulus 4.2 cm an increased size of graphite precipitates. A possible reason for the increased fatigue strength at a casting modulus of 3.3 cm can be seen in the required alloying for a successful heat treatment, which depends on the wall thickness [19]. With 11.7% and 16.7% change in local fatigue strength, respectively, a significant influence of local microstructure and, thus, local solidification and heat treatment can be identified.

Multi-domain optimization approach
The presented multi-domain optimization approach for cast iron components in wind turbines is characterized by the combination of casting simulation, topology optimization, micromechanical analysis, and strength assessment (see Fig. 2). Thus, it uses cross-scale simulation and allows an enhanced consideration of local material properties. Starting from an initial casting simulation, the microstructure of characteristic component regions is estimated and micromechanical analysis is performed for these regions. In parallel, a topology optimization in combination with castability analysis leads to a first optimized topology. In a subsequent optimization loop, the local material pa-rameters, determined by micromechanical analysis, are considered in an iteration of strength assessment, shape optimization, and casting simulation.

Casting simulation
Within the presented approach, the casting simulation is used to estimate local solidification behavior and determine the castability of the optimized component. In contrast to the classical approach of solidification time prediction using the casting modulus, defined by the volume V and the cooling surface S, the thermal modulus (or feedmod) is used. Based on Chvorinov's rule [20] the casting modulus is proportional to the solidification time: where c is a constant depending on cast alloy, mold material, and pouring temperature. Thus, the feedmod considers neighboring parts, auxiliary materials such as cooling chills, or transient heat flow conditions during cooling. Therefore, determining characteristic material regions in the investigated component under consideration of the feedmod allows an implicit consideration of solidification times.

Shakedown analysis
The micromechanical analysis, as a central element of the presented optimization approach, is based on the numerical determination of a shakedown state for a microstructure as a result of a casting simulation. A shakedown state can be reached under cyclic loading if, after initial loading cycles with plastic deformation, a purely elastic deformation occurs during loading. Weichert et al. [21] have investigated the shakedown state for composite materials. According to them, the shakedown state is governed by a time-invariant residual stress field , which is formed in the microstructure [21]. The residual stress field is caused by an accumulation of local plasticity and superimposes the purely elastic K Fig. 3 Exemplary RVE with graphite precipitates (dark) and ferrite matrix (gray) and applied boundary conditions for shakedown analysis stress field σ e due to external loading [21]. The locally acting stress can be defined as: The plastic dissipation is bounded once a shakedown state, characterized by purely linear elastic loading, is reached. This prevents the formation of any further plastic strains in the microstructure [8].
Using Melan's static theorem for shakedown [22], the shakedown state can be determined: "If there exists a factor oe > 1 and a time-independent residual stress field ¡ with R V 0ρ W D WρdV < 1, with compliance tensor D, such that for all loads B in the load domain the following inequality is satisfied f .λσ e +ρ/ Ä 0, then the structure will shake down under the given load domain." The FEM-based evaluation used in this study was implemented by Gebhardt [23], and yields the maximum admissible residual stress field under a given loading. Therefore, the following optimization problem, shown in non-discretized formulation, must be solved: Here, NV represents the number of load vertices that define the load domain, Ω the computational domain, and Γ the boundary of the computational domain.
In order to apply the FEM-based evaluation an applicable computation domain representing the microstructure of interest must be supplied, typically in the form of a representative volume element (RVE). Here, RVE define a computational domain that is sufficiently large to represent the characteristics of the local cast iron microstructure. For a detailed description of the RVE construction, the reader is referred to [9]. An exemplary RVE of size 1.387 × 1.045 mm 2 is given in Fig. 3.
Material parameters of ausferrite and graphite used for shakedown analysis are summarized in Table 4. Graphite and ausferrite were considered to be linear elastic-ideally plastic. The yield strength of graphite is artificially increased to ensure the significance of the ausferrite matrix during shakedown. No debonding is considered at the interface of graphite precipitates and ferrite matrix.

Topology optimization
Topology optimization is used to determine the general component layout within the presented optimization approach. Therefore, following the initial casting simulation, an FE model of the component is created and split into modifiable design and non-modifiable non-design spaces. To ensure the maximum potential of the topology optimization, design spaces should fill as much available space as possible. Finally, the topology optimization can be formulated as minimization of an objective function f depending on a design variable x. Here, gj are the inequality constraints, hk are the equality constraints, and x l i and x u i are the upper and lower explicit constraints, respectively. However, functions have several local minima or maxima and it is not immediately evident whether the global maximum has already been found [24].
While different approaches exist, here the SIMP (Solid Isotropic Material with Penalization) approach is used for topology optimization. In solving the optimization problem, the density in an individual element of the FE-model is continuously varied between a minimum value ρmin and the actual material density ρ 0 . If the ratio of the local element density ρe and the density of the material ρ 0 is formed, the relative density xe is obtained as a design variable: If xe takes the value 1, solid material is present in the component at this point. If the value is less than 1, the local material stiffness in the form of the Young's modulus is reduced accordingly: Here, p is the so-called penalty exponent and Ee is the local material stiffness. It usually covers an interval between 2 and 4 [25].

Strength assessment
The strength assessment of the cast iron components to be optimized in presented approach is based on the international standard IEC 61400 (IEC 61400-4 for wind turbine gearboxes) [26]. This specifies the essential design requirements for ensuring the structural integrity of wind turbines and their components. One requirement is that the strength assessment must be based on validated and recognized methods which consider, for example, the size effect on material strength.
The strength assessment in this optimization approach is performed according to the FKM guideline. The process of strength assessment according to FKM follows a defined procedure for determining material parameters and component design parameters up to the degree of utilization [27]. In this work, the process was automated and extended to include local material parameters.

Results and discussion
The concept of the presented optimization approach is demonstrated using a planet carrier made of ADI-GJS-1050-6. Therefore, first an initial topology optimization was performed and subsequently its effect on the local microstructure was analyzed. Finally, the effects of local fatigue strength on the utilization were analyzed.

Topology optimization of planet carrier
In order to show the maximum possible optimization potential, the initial geometry was reworked prior to topology optimization. This geometry extension is to ensure that the maximum possible weight reduction is achieved during topology optimization. Figure 4a shows a sectional view of the planet carrier. Here, the yellow part indicates the ini- tial geometry and the magenta part the volume extension for topology optimization. The outer diameter of the planet carrier was slightly enlarged compared to the initial geometry. Furthermore, the inside diameter was reduced to the diameter of the tip of the sun gear.
The design space available for topology optimization is visualized as blue regions in Fig. 4b. It is limited by the optimization restriction of constant local stiffness at the mounting for the planetary gears. The red area shows the non-design space, i.e. the part that remains unchanged after the topology optimization. Due to the increased design space, the mass of the planet carrier increases by 22%.
Boundary conditions and loadings considered during the topology optimization as well as for the determination of stresses for the strength assessment are visualized in Fig. 5.
Here, all degrees of freedom of the coupling surface were fixed. The maximum torsional load from the power production on normal turbulence design load case (DLC) 1.2 [28] was applied via circumferential forces acting on the planetary gear pins. The bearing surfaces of the planet carrier were only fixed in the radial direction. The topology optimization was carried out using constant material properties Concerning system stability, weight reduction was facilitated under the boundary condition of constant stiffness during power production load case DLC 1.2 and castability. The system stability requirement is intended to ensure that the system behavior remains unchanged after optimization for the power production load case. In this context, system stability was determined by evaluating allowable deformation and component stiffness. Therefore, seven different optimizations with different target volumes were carried out to determine the resulting change in torsional stiffness of the planet carrier for a given change in weight compared to Fig. 7 Change in torsional stiffness vs. change in weight of the planet carrier compared to the initial geometry for torsional loading Fig. 8 Result of the topology optimization of the planet carrier the initial geometry. Target volume vi specifies how much of the extended geometry volume remains after optimization. Figure 6 illustrates some results of topology optimization for different target volumes. Figure 7 depicts the relation between the change in torsional stiffness and weight of the planet carrier compared to the initial geometry for different target volumes. Positive values imply an improvement in torsional stiffness and an increase in the weight of the planet carrier. It can be seen that the torsional stiffness of the planet carrier can be maintained despite a reduction in mass. Moreover, compared to the initial geometry, an increase in torsional stiffness of about 16% can be achieved with the same mass. If the torsional stiffness is kept constant, a weight reduction of about 17% is achieved.
According to Fig. 7, a target volume of slightly less than 70% provides a torsional stiffness similar to the initial planet carrier design. For this reason, this geometry resulting from the topology optimization is chosen as the basis for the new design. In order to exploit the highest weightsaving potential, the planet carrier is reconstructed as close as possible to the optimized geometry. The first design of the planet carrier can be seen in Fig. 8. Compared to the initial geometry, a weight reduction of 17.6% is achieved.
As shown in Fig. 8, pockets were formed at the front and rear of the planet carrier, which protruded into the body. In top view, an x-shape structure results, caused by the fact that the planet carrier is predominantly subjected to torsion, so that mainly shear stresses have to be transmitted. Since the main stress direction is rotated by about 45°to the center axis, diagonal structures are formed. The structure described is smaller on the gearbox side, since a large proportion of the load is already distributed over the planetary gear pins. Noticeably, the pocket extends to the bearing seat and even forms a small opening. Comparing the influence of the target volume, only the intensity of the described characteristics varies. Figure 9 shows the deformation within the planet carrier and the planetary gear pins under the same loading for the initial planet carrier geometry and the 17.6% weightreduced planet carrier with the same torsional stiffness. As  indicated by the constant stiffness, no significant difference in local deformation can be observed. However, these values represent theoretical optimum values that may not be achieved in practice due to manufacturing-related design changes or cost considerations.

Analysis of topology optimized planet carrier with respect to manufacturability
As the planet carrier is made out of ADI-GJS-1050-6, component optimization must also consider the characteristics of the manufacturing route, e.g. heat treatment, besides the overall goal of weight reduction while maintaining system stability. Concerning the manufacturing route, the local distribution of feedmods and efficiency of the heat treatment were analyzed.
Since manufacturing of ADI depends on both initial raw material and heat treatment, the melt is adjusted to reliably set ADI microstructure in the critical component regions. Thus, to evaluate the component optimization with respect to the working material, the feedmod distribution and the phase fraction of ausferrite after heat treatment were analyzed. Figure 10 shows that the scatter of feedmods can be significantly reduced by reducing the maximum feedmod from 8.4 to 6.9 cm. A reduction of the feedmod range allows a better adjustment of the melt, especially regarding hardening agents and other alloying elements. The minimum of non-transformed austenite (after quenching), as shown in Fig. 11, increases from 2 to 98%. Thus, large non-ausferrite regions could be eliminated and a more successful ausferritization can be expected.

Strength assessment with local fatigue strength
Considering the complex adjustments required for successful ADI manufacturing, the importance of local material The strength curves as a function of feedmod are shown in Fig. 12. Here, the non-linearity of the fatigue strength with respect to the feedmod can be observed. This nonlinearity may be associated with the interdependence of graphite precipitates and ausferrite matrix. This may be related to the ADI manufacturing, which is governed by a heat treatment process window and special requirements regard- ing alloying of wall-thickness dependent elements [19]. Thus, optimal material parameters for the given use case may only be present at a specific feedmod, whereas lower feedmods may have passed the ideal process window while higher feedmods have not have reached the ideal process window. Due to limited experimental validation, the maximum feedmod of 5 cm was set as the bound of the validity range.
The locally resolved ultimate tensile strength as a function of the feedmod and the fatigue strength estimated using this ultimate tensile strength are shown in Fig. 13. The conversion of ultimate tensile strength into a fatigue strength (hereinafter referred to as FKM-fatigue strength) was carried out following the FKM guideline [27]: Highest ultimate tensile strength and fatigue strength can be observed in regions with lowest feedmods, e.g. close to the surface of the component and at the corners. However, this relation does not account for the observed non-linearity of fatigue strength with respect to the feedmod.
In Fig. 14, the locally resolved FKM-fatigue strength and fatigue strength according to shakedown analysis are compared within the critical notch. It can be seen that the regions of the highest fatigue strength can be found at a given offset from the surface of the component if fatigue strengths based on shakedown analysis are considered. Especially close to the surface a difference of 15% between FKM-fatigue strength and fatigue strength based on shakedown analysis can be observed.
The local utilization is shown in Fig. 15. For both cases, FKM-fatigue strength and fatigue strength based on shakedown analysis, the local utilization rates remain below 20%. This is related to the optimization requirements of the planet carrier, so local fatigue damage can be neglected for the given combination of considered loads and material. However, if the local utilization increases, the local fatigue strength within a component of austempered ductile cast iron becomes significantly more relevant.

Conclusion and outlook
This study demonstrates the consideration of local material properties in component design of heavy-section castings used in wind turbines. An improved understanding of the consideration of local material properties, especially for complex material systems such as austempered ductile cast iron, was achieved. Subsequent component optimization, resulted in significant weight reduction under the overall constraints of castability and system stability.
Heavy-section castings of austempered ductile iron grade ADI-GJS-1050-6 were characterized. Results of material characterization were used as calibration parameters for the presented multi-domain optimization approach, which was demonstrated using a planet carrier. Thereby, casting simulations have been used to compute local microstructure characteristics in critical regions of the planet carrier. For a set of supporting points of these microstructure characteristics, the shakedown analysis was used to numerically determine microstructure-dependent fatigue strengths. In parallel, a topology optimization of the planet carrier was performed and a significant weight reduction potential of 17% was identified while maintaining the required initial component stiffness. The effect of topology optimization on local material quality was determined by comparing feedmod scattering and expected local ausferrite phase fraction. Here, the topology optimization resulted in an improvement of both feedmod scattering and local ferrite phase fraction. Thus, it was shown that topology optimization reduces component weight and can increase the local material quality within the component for austempered ductile iron.
The local utilization was computed using fatigue strength estimates based on the ultimate tensile strength and fatigue strength obtained by shakedown analysis. Here, good agreement between both with a difference of 15% could be observed. It was shown that analysis of local fatigue strengths allows the derivation of required changes in manufacturing, e.g. the adjustment of the heat treatment process to achieve the highest fatigue strengths close to the surface. Overall, the presented multi-domain optimization approach resulted in an optimized planet carrier with a weight reduction of 17% while maintaining the required stiffness and reduced microstructure variance along the component. Due to the high strength of the used material ADI-GJS-1050-6, the strength assessment was neglectable compared to the required stiffness for the optimized planet carrier. Thus, the strength assessment was irrelevant to the demonstrated weight reduction. Nevertheless, this method can be applied to other heavy-section castings where consideration of local material properties during strength assessment may be advantageous for the component design.
Future work will continue to integrate local material properties into strength assessment concepts of heavy-section castings used in wind turbines. The applicability of the approach needs to be validated by investigating other geometries and other austempered ductile cast iron grades. The estimation of local fatigue strength using the shakedown approach will be extended to include a criterion to account for local permissible damage. To reduce safety factors in the future, the interaction of safety factors with local material properties and local material quality will be determined.
Funding This research has been funded by the Federal Ministry for Economic Affairs and Climate Action (BMWK) in the project LeKo-Guss WEA under grant number 0324279A. The authors acknowledge the support by Siemens PLM, agreement 60068580, for providing STAR-CCM+ licenses.
Funding Open Access funding enabled and organized by Projekt DEAL.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4. 0/.