The present study, which was carried out in collaboration with GKN Fokker, focuses on incorporating bird strike crashworthiness requirements within a multidisciplinary optimization (MDO) framework. During the preceding three-month internship in the same company, a pivotal contribution to this project was the development of an Abaqus interface for the Multidisciplinary Modeller, MDM, created within the Center of Competence in Design department. MDM is a Python/ParaPy-based automated generator of wings, moveables and flaps, starting from a set of user-specified parameters. The generation of ready-to-run input files thus lays the foundation for the subsequent optimization process, as any changes in materials or geometry can be easily accommodated.

The core objective of the research is to minimize the weight of an aircraft wing while taking into account additional requirements related to the extent of damage caused by bird strikes. Unfortunately, such events occur more frequently than one would be comfortable with, and stringent requirements are set in place to guarantee the safety of the passengers. Among these requirements, the aircraft must be capable of landing safely after such an event, being subject to loads associated with get-home conditions.

As a consequence, two critical constraints are formulated within the optimization framework, addressing the residual strength of the damaged front spar following a bird strike, coupled with a requirement based on a maximum penetration depth. The last constraint has also been included due to the rising popularity of the electric vertical take-off and landing aircraft, which not only fly at low altitudes, thus increasing the risk of bird strike, but may also contain battery packs in the leading edge, for instance, which can pose a significant risk if damaged. To tackle the complexity of this highly-dimensional optimization challenge, a methodology based on Bayesian optimization is proposed, employing surrogate models coupled with a preliminary variable ranking procedure.

The Kriging metamodel is identified as a suitable candidate, thanks to its error prediction capabilities, which are paramount in Bayesian optimization. A variance-based dimensionality reduction method is proposed, which makes use of an initial surrogate to estimate the main and interaction effects of the variables. The quantification of the significance of a variable is expressed as its percentage contribution to the total variance, thus allowing for an intuitive selection of the most important parameters. After the screening procedure is complete, the optimization procedure is carried out in the reduced design space, which uses the constrained expected improvement as an acquisition function. The proposed methodology is then applied on a case study problem, involving a five-bay metallic wing segment subject to the constraints aforementioned, involving 19 design variables representing the thicknesses of various components.

Remarkable weight savings have been achieved, the final result being 40\% lighter than the lightest feasible design among the initial data points. A significant dimensional reduction has also been attained for the maximum depth constraint, which is expected due to the local nature of the impact. Not only did the number of variables greatly decrease from 19 to just 3, but a considerable increase in the accuracy of the corresponding metamodel has also been registered, thanks to an increase in sampling density in the reduced space. However, the variable screening procedure revealed intricate interaction effects with respect to the residual strength of the front spar, emphasizing the nuanced complexity inherent in crashworthiness considerations. Nevertheless, a moderate dimensional reduction has been achieved for this constraint as well, reducing the number of variables to 8, thus proving the efficacy of the proposed variable screening procedure.

In conclusion, the utilization of Kriging models, variable ranking procedures, and Bayesian optimization collectively contributed to the success of achieving remarkable weight savings, proving the efficiency of the proposed methodology. Moreover, it has been shown that the integration of a residual strength requirement is necessary, as many cases were uncovered where no significant penetration occurred, although the application of the considered load case, which is not from critical to an undamaged wing, resulted in high stresses to the front spar of the damaged structure.