This thesis examines the buckling behavior of sandwich composite cylindrical shells, which are integral to the primary structure of launch vehicles. The central aim is to develop a scaling methodology based on nondimensional parameters, enabling large-scale composite structures to be scaled down to laboratory-size models while maintaining their buckling response.

The study begins with an overview of existing analytical, numerical, and experimental techniques for shell buckling analysis, alongside a review of current scaling methods in structural mechanics. Building on insights from the literature, a scaling framework is introduced by reformulating classical buckling equations in a nondimensional form. This approach allows the scaling laws to be directly derived from the equations' components but requires a comprehensive and adaptable nondimensional formulation of the structural behavior.

To establish this formulation, the work extends the nondimensional framework by first incorporating the effects of transverse shear deformations, relevant when the shells are reduced in size. Then the framework is extended to included the theoretical impact of imperfection sensitivity, tackled by incorporating a trigonometric imperfection model into the nondimensional framework.

With these extensions in place, a systematic scaling methodology is proposed. Two distinct strategies are developed: one that directly scales sandwich composite shells while preserving their structural characteristics, and another that substitutes them with equivalent composite laminate shells. While the first strategy offers greater theoretical accuracy, it is constrained by practical limitations in manufacturing scaled thicknesses. The second strategy, although more feasible for experimental implementation, introduces new complexities in ensuring equivalence between different structural configurations.

The proposed methodologies are validated through comparisons between analytical predictions, numerical simulations, and experimental results. Scaled laboratory models, produced using the more practical laminate-based strategy, show an 8% discrepancy in nondimensional buckling loads in theory compared to full-scale counterparts. However, experimental observations reveal a larger deviation of approximately 22%, underscoring the limitations of current imperfection modeling and the need for refinement in the scaling approach.

The thesis concludes by emphasizing the contributions made to nondimensional scaling methods for composite structures and highlights the importance of further experimental validation. In particular, improved modeling of imperfections and expanded laboratory testing are recommended to bridge the gap between theoretical predictions and real-world behavior. This work lays a solid foundation for the development of scaled testing protocols, advancing the design and verification of large-scale sandwich composite shell structures in aerospace applications.

This thesis project focuses on the evaluation of the Extended Finite Element Model (XFEM) for modelling skin-stringer separation in thermoplastic composite stiffened panels. The study uses test results from a reference study in the literature on a PEKK-FC carbon composite panel with three stringers joined to the skin using a short-fibre reinforced butt joint. The modelling strategy was developed using Double Cantilever Beam (DCB) specimens. The XFEM-based panel models were then developed using two different damage initiation criteria: Quadratic Stress Criterion (QUADS) and Maximum Principal Stress Criterion (MAXPS). These models were compared with the Virtual Crack Closure Technique (VCCT) model from the reference study. The project also addresses the challenges posed by the crack tip impingement problem in the XFEM-based model with MAXPS and tests different strategies to overcome it. A proof of concept is provided for using the UDMGINI subroutine to address these challenges, setting the stage for future studies and algorithm development.

This dissertation explores the development of a framework for wing aeroelastic tailoring in the preliminary design phase of aircraft, focusing on weight reduction and aerodynamic efficiency through the use of composite materials. The framework utilizes low-fidelity, fast computational tools to handle numerous load cases, modeling the wing structure with FEM shell elements and aerodynamics with the doublet lattice method (DLM). Composite material properties are optimized using lamination parameters, with structural continuity enforced through blending constraints.

To address the limitations of low-fidelity methods, CFD-based corrections are applied to the DLM for accurate aerodynamic load representation, and a surrogate model is introduced for aerodynamic drag estimation. This model approximates wing deflection efficiently, reducing the need for costly high-fidelity simulations.

The framework, implemented in Nastran, is applied to industrial cases focusing on the impact of blending constraints, dynamic load cases, and multi-objective optimizations balancing weight and drag. Blending constraints result in smoother transitions between wing sections, ultimately producing a light final design, even though there is an initial 5% increase in structural weight due to a reduced design space. Dynamic loads, incorporated through the Equivalent Static Load (ESL) method, have minimal impact on the overall design.

Multi-objective optimization reveals trade-offs between structural weight and aerodynamic drag. Weight minimization favors structures with wash-out behavior, leading to higher drag due to increased angle of attack. Conversely, drag minimization results in stiffer wings with less wash-out. The framework's results align closely with full-order CFD models, validating its effectiveness.

Offshore installation vessels are growing in size in order that they are capable of handling the next-generation 15 and 20MW wind turbines. At the same time, it is important to gain a better understanding of their structural behaviour using the finite element (FE) method. One approach that is often used is the full vessel FE analysis, which studies the vessel in its entirety. However, the drawbacks are that a coarse mesh must be adopted to have an acceptable simulation time and it takes months to create a 3D model of a vessel. Another approach is the partial vessel FE analysis, which considers only a section of the vessel. This research aims to develop a methodology for the partial vessel FE analysis of an offshore jack-up installation vessel.

The studied section of the vessel is between the four legs, which is where the highest bending moment is expected and the wind turbine towers are located on deck. The methodology was developed to analyse the vessel in still water and waves separately. The boundary conditions and hull girder load adjustments of the partial vessel model were described, which cause the vessel section to reflect the behaviour of the full vessel. The procedure to apply the loads on the model was also given. The finite element analyses were performed and the stress results compared to a simple analytical solution of the vessel based on Euler-Bernoulli beam theory.

The methodology of the partial vessel FE analysis was successfully applied to the offshore jack-up installation vessel. For the vessel in still water, the stress results were in agreement with the analytical solution. The results were accurate in the middle region for approximately 50% of the partial vessel model and thus it is advised to have a larger model for future use. The effect of torsion could not be fully justified as it is not captured in the analytical solution, but resulted in a significant improvement. For the vessel in waves, the stress results also showed a good agreement with the analytical solution. The accuracy depended on the equivalent design wave approach, which assumed a regular travelling wave. Lastly, it was proven that the presence of payload only leads to local changes in the stress field, but not in the global behaviour of the vessel.

The issue of space debris in Low Earth Orbit is a growing concern that requires comprehensive solutions. This includes preventing further pollution, removing existing debris, and designing resilient spacecraft that can withstand impacts. This thesis focuses on the improvement of spacecraft shielding structures to provide effective protection against hypervelocity impacts. The Smoothed-Particle Hydrodynamics method was used for the simulation of hypervelocity impacts using LS-DYNA. The method is preferred for the simulation of impacts at high velocity and was used to reproduce two studies from the past. The simulation was then performed on a double plate system, with results compared to data from experiments provided by Airbus Defence and Space GmbH, under co-supervison from ESA. The simulation accurately predicted the hole diameter in the 1st plate with less than 2% error, while the damage zone in the second plate showed considerable variance. The simulations consistently overpredict the damage zone, leading to conservative results.

The rudder is one of the few aircraft primary composite structures that experience relatively low external aerodynamic loading during regular flight conditions. This renders the rudder shell the perfect candidate to study the possibility of introducing sustainable materials with lower mechanical properties in primary structures. However, high internal loads are introduced in the shells during actuator jamming events, which occur with low frequency. Locally redesigning the rudder spar as a torsion box to take up the jamming loads and alleviate the shell loading could enable the safe usage of green materials in the rudder shells. The aim of this thesis is to study the design feasibility of three rudder torsion-box spar concepts to separate the aerodynamic and jamming loads in the rudder and provide a preliminary optimised structural design based on a parametric study of the torsion-box spar design variables and static FEM simulations.

This thesis aims to study the design and optimization of composite wing components for aircraft, focusing on the pivotal goal of reducing weight while preserving and enhancing structural performance. This is one of the goals of the COST Action joined in the span of the thesis: CA18203 - Optimising Design for Inspection (ODIN). Within this COST Action, a comparison between different Finite Element modeling predictions of the response of this structure subjected to 4-point loading has been presented
by Bisagni et al. [1] before this thesis. The thesis, which corresponds to the optimization of this same structure, is the continuity of this work. The foundation of this thesis is built upon a literature study that investigates the optimization of composite wing structures. It explores diverse optimization formulations, including design variables, objectives, and constraints. Different methodologies for optimizing wing structures, like multi-objective and probabilistic methods, are studied. Various algorithms used to tackle these challenges are presented, shedding light on their advantages and drawbacks. Once a solid knowledge of wing structure optimization is built, preliminary analyses (study of coupon geometries, failure criteria, and the buckling behavior of various materials and geometries) form the basis for the upcoming optimization strategy. These preliminary investigations refine our understanding and guide the selection of optimal constraints, such as failure criteria. Based on these preliminary analyses and the literature study, the optimization strategy has been established. A new configuration with cut-outs is developed to explore new innovative designs. To solve the optimization problem, Genetic Algorithms, inspired by real-life behaviors and processes, emerge as the optimization algorithm of choice. These algorithms contribute to the development of design solutions that are better than traditional designs. The objective is to achieve weight reduction in the structure while considering its structural performance. Consequently, the fitness function, employed to assess candidate designs within the Genetic Algorithm, is based on three primary factors: weight, stiffness, and buckling behavior. Weight carries the most significant part of the fitness function, accounting for approximately 80% of the total value, whereas the remaining two factors contribute approximately 10% each. To evaluate these aspects, weight is determined through basic calculations, while stiffness and buckling behavior are assessed by simulating a 4-point bending test using Abaqus.
The implementation of the Genetic Algorithm, adapted to our special case, allows us to showcase its effectiveness in achieving optimal designs. These optimal designs achieve remarkable weight reductions while maintaining great structural performance. Following 20 generations, each comprising 10 potential candidates, one of the two termination criteria is triggered: there is no improvement in the best solution for eight consecutive generations. The optimized solution obtained is 16.89 kg lighter compared to the baseline configuration studied by Bisagni et al. [1]. This represents a significant 25% reduction in structure weight. Concerning the structural performance, only a reduction of 1.5% and 0.25% is observed for the stiffness and buckling performance, respectively, which is acceptable regarding the weight reduction. Nonetheless, achieving such a reduction in wing weight may not be realistic as it was specifically tailored for the test case of this thesis, other loading scenarios need to be studied to assess its applicability.
These outcomes underscore the potential of optimization techniques in aerospace engineering, paving the way for the development of lightweight and high-performing composite wing components.

The exponential growth of the aviation sector combined with environmental and energy supply challenges have led companies to innovate to reduce aircraft fuel consumption. Among the many areas currently under research, this thesis investigates the possibility of making aircraft wings lighter thanks to flight load alleviation systems and more accurate modelling methods.

The development process of an aircraft involves using low-fidelity aerodynamic models at the early stage of the design process to rapidly compute the loads acting on the airframe, and to evaluate the efficiency of wing control surfaces. These models are, however, limited to linear flow conditions, and transonic shock or flow separation cannot be simulated with such methods. This requires important safety factors and leads to generally heavier designs. Computational Fluid Dynamic with Reynolds-Averaged Navier-Stokes (CFD-RANS) analysis is capable of better aerodynamic predictions, but the computational time required for such simulations is too long to be efficiently included in the sizing process of the airframe. The approach proposed in this thesis aims to combine the accuracy of CFD with fast linear loads estimation. This is achieved by deriving reduced-order models (ROM) of the aircraft control surfaces and manoeuvre loads from rigid CFD analysis to improve the accuracy of faster but lower-fidelity results where needed. These fast aerodynamic models for the control surfaces also allow rapid control optimisation to evaluate their load alleviation potential.

The thesis starts by introducing and validating the unsteady and non-linear models with 2D examples. Then, it covers the application of these models to a flexible 3D wing. The models are validated against high-fidelity steady and dynamic Fluid-Structure Interaction simulations and show good agreement with a 5% to 10% error margin in loads and deformations in most of the cases. Finally, a wingbox sizing optimization is performed with active load alleviation. Choosing to either use the linear or the non-linear aileron model for the GLA alone leads to a 2.5% difference in the wingbox structural weight.

A building block pyramid is designed used to evaluate a composite aircraft structure, combining tests and analyses, using stiffened panels, single-stringer specimens and coupons. Failure in these structures can be caused by postbuckling-induced skin-stringer separation, which is complex, involving matrix cracks, fibre bridges, and delamination migrations. Standardized tests covering postbuckling-induced skin-stringer separation are lacking. First, a material characterization is performed on the coupon level. These properties are applied to a stiffened panel model to identify critical postbuckling regions. Then a single-stringer specimen is designed that combines the material complexities of the coupon level and the geometrical complexities of the panel by mimicking the postbuckling shape. Specifically a seven-point bending configuration to study bending-induced separation and a four-point twisting configuration for skin twisting-induced separation. The guidelines for modelling and testing can assist in standardising these test methods.
Using composite materials in aircraft structures can reduce weight compared to conventional metals. However, utilizing more of the material's load-carrying capabilities can further reduce weight.

Thermoplastic composites enable new manufacturing techniques such as conduction welding to make the aviation industry more sustainable, while at the same time, provide great benefits to cost-efficient high-volume production. One of the benefits of welding is that it reduces the amount of mechanical fasteners required. Fastener-free joining also poses new challenges, because the performance of these highly loaded structural joints relies heavily on the performance of the thermoplastic polymer matrix. Furthermore, there is currently not much understanding of the mechanisms involved in thermoplastic welded joint failure, and the numerical and experimental methodologies, originally developed and validated on thermoset composites, have not yet been fully assessed for thermoplastic composites. On top of that, the process conditions to manufacture these new structures may have a significant influence on the mechanical performance of the material and can thus play an important role in the design of thermoplastic composite structures.
The objective of this research is to analyse matrix dominated failure of thermoplastic composites and conduction welded joints and to develop both experimental and numerical methodologies to support the design of thermoplastic composites structures. The research addresses important linkages between the three main pillars of Manufacturing, Experimental and Numerical analysis....

The state-of-the-art design approach of composite structures in Formula One (F1) rely on expensive full-scale crash tests and simple models that cannot describe the underlying physical phenomena during crushing. Although the models can describe the crush forces and resulting accelerations, their predictive capabilities are poor, and several iterations are usually required to arrive at a satisfactory design, causing the teams money and resources. Reductions in costs could be possible by reducing the number of full-scale experiments and increasing the amount of analysis.

This thesis aims to evaluate the capabilities of a commercially available material model at two distinct modelling scales by conducting physical and virtual experiments on two different specimen geometries; tubular coupons and a new open-section channel specimen. Firstly, the new open-section specimen with a geometry representing the typical F1 crash structures was designed, manufactured and tested. The newly proposed specimen showed controlled failure modes, and displayed typical characteristics of composite structures under crushing loading. Secondly, the two specimen geometries were numerically analyzed on the macro-scale (laminate-level) and meso-scale (ply-level) using a commercially available material model implemented in the Abaqus finite element software. Mesh-dependency of the macro-scale approach was discussed and numerical aspects required to perform crush simulations on the meso-scale were identified and discussed. Lastly, calibration of the numerical models was performed to understand the effects of model inputs on the output responses, to evaluate the validity of the input parameters, and to assess two different calibration approaches; the first method is based on macro-scale modeling, while the second method is based on surrogate-based inverse parameter identification at the meso-scale.

The numerical models could not reproduce the experimental crush responses of either geometry when using experimentally characterized material data. Following the second calibration approach, the qualitative and quantitative crush responses of the tubular specimen could be accurately described, but the channel specimen showed no improvements in their response. The macro-scale modelling approach with the current material model was judged to be impractical for large-scale analysis, and the meso-scale approach showed need for further developments.

The safety of a commercial aircraft is a factor to consider from the early stages of its design to guarantee that this one's structure can protect the occupants inside in a low-speed crash. In an aircraft's primary structure such as a fuselage, where the occupants and cargo should be protected, the design regulations imposed by the EASA have been applied for decades to conventional metallic designs. These designs had aluminum as the main constituent for being a lightweight material. Nevertheless, this project improves a component's design from a fully thermoplastic fuselage. While two commercial aircraft (A350, and B787) have already implemented some composite parts in their design to aim for sustainable and lightweight components, thermosets were the composites used in primary structures mainly. Composite materials do not present the plasticity metals do, which means that the regulations from EASA can no longer apply in the A350 and B787 hybrid designs. With the introduction of composite materials in these aircraft, organizations like the FAA and EASA had to write some new guidelines for these aircraft structures. These enable the protection of the passengers and cargo from the loads exerted on the aircraft during a survivable crash. Therefore, the new regulations are the ones taken into account for this thesis project, since these can be applied to similar composite structural designs, such as the fully thermoplastic fuselage section, whose numerical model is inherited from the Clean Sky STUNNING project.

What is wanted for the fuselage section's crashworthiness is for its structure to absorb as much energy as possible from the crash. All while minimizing the peaks of force that this one experiences. That way, the loads that arrive at the cabin of the aircraft are reduced and the structure becomes safer for a low-speed crash of under 30 ft/s (about 10 m/s), meaning that the passengers inside have a higher chance of survival and it is less probable they suffer from severe injuries. That is the aim of this project, focusing only on the structure's behavior caused by the changes in the fuselage struts.

The numerical analyses performed in LS-DYNA increase single components’ structural complexity until optimizing the layup of a tube to improve its crashworthiness behavior. It is after this optimization that the component is introduced as a strut in the STUNNING fuselage section. To further improve its global crushing in low-speed crushing conditions, the struts (or energy absorbers) keep the square geometry cross-section from the tube exercise, as well as the optimized layup. However, these are changed in size to improve the fuselage’s crushing behavior with little to no increase in mass. This study considers a fuselage section that neglects the passengers’ mass, their cargo, and the upper part of the fuselage airframe in the section's model. Considering only the cargo mass for six passengers on the model causes a change in the fuselage’s final design crashworthiness, which proves that a more representative crushing conditions setup should be considered in future studies to better predict the structure’s crushing behavior numerically.

A novel twist morphing concept is explored, that exploits the buckling instabilities of slender spar webs integrated into a wing structure to control the twisting response to external loads. The main novelty of the concept is in controlling the effective shear stiffness of the post-buckled slender spar webs through external variable constraints acting on the out-of-plane buckling deformations. A methodology for the design and analysis of these novel morphing structures is proposed and implemented in the design of a wing box structure of promising twist morphing capabilities. The overall design process is structured into a multilevel process of increased complexity. In the first level, the morphing structure is simplified to a wing box with slender spar webs. With the objective of maximizing the morphing twists that can be achieved under the action of an external quasi-static torque, the wing box design space is explored in terms of its cross-sectional dimensions and the material assigned to the slender spar webs. In the second level, the morphing structure is expanded to include both the wing box and the external devices required to implement the adaptive constraints acting on the slender spar webs' out-of-plane buckling deformations. At this level, the objective is to design adaptive constraining devices that maximize the twist morphing capabilities, for which the influence of the constraining devices over the twisting response and their effectiveness in restraining the slender spar webs' buckling deformations become the main concerns. After an extensive design process, a design solution for the adaptive constraining devices is proposed, for which thorough analyses on the twist morphing capabilities are performed. In addition, the sensitivity of the twist morphing capabilities to the slender spar webs' geometrical imperfections is investigated.

The goal of this thesis is to build an artificial neural network(ANN) surrogate model, that predicts the crashworthiness performance of a structure. The structure used in this research is a thermoplastic fibre-reinforced composite aircraft subfloor section which is part of the SmarT multifUNctioNal and INteGrated thermoplastic fuselage (STUNNING) fuselage demonstrator section. The structure will be analysed numerically, by making use of the finite element software LS-Dyna by Ansys. Within this research, a surrogate model is defined as a model that has been fed data from a computationally expensive model, yet is computationally considerably more inexpensive to run to obtain similar results as the expensive model. The main methodology to do so in this research is to build a two-dimensional design space, which is sampled using optimal latin hypercube sampling. By using an existing model for the STUNNING fuselage subfloor section and by automatic generation of geometry changes of the structure through a Python script, input files can be created according to the sampled design space. The simulations are then run on a high performance cluster. Once the simulations are run the results can be fed to an artificial neural network that uses this data to predict crashworthiness performance of the structure within the design space.

To build towards the goal of the thesis, a small comparative study has been performed on material models MAT054, MAT058, and MAT 261 in LS-Dyna. These material models are based on different failure criteria, and thus have different performance in coupon simulations. Using these simulations as a basis, a simple design space was created, sampled, and modelled so that a first ANN surrogate model could be built that could predict the failure displacement and the ultimate load of a composite coupon in tension. It was found that for a small dataset (<10 samples), the neural network could easily fit to the data, and could predict accurately on samples within the design space that were withheld from the network during training.

Then, the STUNNING subfloor section model was introduced and a number of models that section this structure are investigated. The main purpose of this investigation is to find out if the model can be reduced in size without changing the resulting failure mode. When a satisfactory reduced model was found, the research carried on by doing a sensitivity analysis on a proposed two-dimensional design space, and similarly a sensitivity analysis on the parallelization settings. From this study, a second failure mode of the model was discovered and the effect of the parallelization settings was quantified. From here the limits of the design space could be formed and sampled. Finally, different variations of ANN surrogate models were defined, trained, and tested. In the end, a root mean square discrepancy of around 10\% was achieved during testing of the network. This means that the methodology could prove useful in design applications, depending on the design stage and needed accuracy.

In recent decades, fibre-reinforced polymer composite structures have seen increased aerospace applications. A typical composite structural element widely used in primary load-bearing aircraft components is the T-joint, which is formed when two orthogonal laminates are joined in a T-shaped arrangement to enable out-of-plane load transfer. Under tensile loading, the structure tends to fail by the separation of the horizontal and vertical laminates at the joining interface. The incorporation of an interleaf, i.e. a thin layer of material with good interfacial bonding properties, may prevent or delay the failure event. A potential interleaf material is the carbon nanotube (CNT). In the literature, its addition as a reinforcing phase in interleaves has been shown to improve fracture toughness in standard coupon tests, but the technique has yet to be tested on a structural element level.

This thesis research investigated the effect of CNT interleaves on the mechanical behaviour of composite T-joints. Manufacturing processes were developed to produce test samples consistent in quality. T-joint samples with CNT-modified interleaves were manufactured alongside neat samples with zero CNT content. The joints were formed using custom-built tools and co-cured in an autoclave. Tensile tests were carried out and the T-joints were loaded until failure. Load-displacement data were acquired and further supported by results from a Digital Image Correlation (DIC) system.

Data analysis suggested that CNT interleaves modified the mechanical behaviour of T-joints by delaying the failure of the structure. While their stiffness and ultimate strength were largely unaltered, the CNT-modified T-joints were able to undergo higher displacements and absorb higher strain energy before failure compared to the neat samples. The delay of failure was enabled by extensive crack development within CNT interleaves preceding the final rupture. In comparison, the neat samples failed abruptly after very limited damage evolution. Observation of the crack surfaces revealed evidence of tortuous crack paths through CNT interleaves. The finding suggests that CNT interleaves are capable of improving the failure resistance of the T-joint structure in a targeted fashion.

The collapse of composite stiffened panels in compression is complex and involves the interactions of many different failure modes. However, when the panel is allowed to enter the post-buckling regime, skin-stringer separation due to the interaction of the post-buckling deformations with the skin-stringer interface is often the critical failure mode. The accurate prediction of skin-stringer separation in these types of panels is therefore crucial for their design. Due to the high computational costs associated with modelling damage in a large multistringer panel and the manufacturing and testing costs associated with testing such a panel, single stringer specimens that accurately represent the behaviour of critical regions in these panels have been proposed in literature. In this thesis, transversely loaded single stringer specimens are designed and verified and are used to predict skin-stringer separation in the critical regions of a specific multistringer panel. A specific composite stiffened panel is considered that contains a skin of a tape material and four co-cured stringers of a fabric material. At the intersections of the skin and stringers, resin-rich noodle regions are created. A buckling analysis is first performed on the panel to obtain its post-buckling deformations. Based on these, two types of regions where skin-stringer separation may occur are identified; a mode I dominated skin-stringer separation at the minima of the buckling waves and a combined mode II and mode III skin-stringer separation at the inflection points of the buckling waves. A seven point bending specimen is designed based on the panel deformations near the former critical region, while a four point twisting specimen is designed based on the deformations near the latter region. Finally, detailed damage models are set up for the single stringer specimens to predict skin-stringer separation and the models are verified with the multistringer panel. The analyses showed that a mode I dominated skin-stringer separation is critical for this panel. The detailed analysis of the four point twisting specimen showed that no damage occurred at the inflection point. To predict the skin-stringer separation at the minimum, a verified seven point bending specimen has successfully been obtained. To do so, it was important to accurately simulate the deformation of the skin at the stringer flange edge, since this proved to be driving for the initiation and propagation of skin-stringer separation. Next to that, a trilinear cohesive law was used to capture the R-curve effect in the tape/fabric interfaces. Finally, a first assessment of the noodle region behaviour showed that damage in this region will probably happen at the location of maximum bending. Recommendations for future research have been provided.

Launch vehicle structures are commonly composed of cylindrical and conical shells, which are inherently sensitive to buckling. The axial compression experienced during launch can consequently be a sizing load case, so it is important to understand the axial buckling behavior of these shells. Experimental testing is an essential part of studying this phenomenon because manufacturing imperfections can cause large discrepancies between theory and reality. In addition, conical shells have been researched less frequently than cylindrical shells, such that their behavior is less well understood. Experimental testing of launch vehicle structures is difficult and expensive due to their large size, hence it is preferred to test reduced-scale shells, representative of the full-scale ones. In this thesis, a scaling methodology is developed which allows designing representative reduced-scale conical shells for full-scale composite conical shells buckling in axial compression. The conical shells are assumed to have a symmetric, balanced layup with negligible flexural anisotropy. The scaling methodology is developed using the nondimensional governing equations, obtained through Nemeth's procedure, which allows to directly use the coefficients of the equations as scaling parameters. It also provides a framework to not only compare the buckling load of the shells of different sizes, but also the displacement upon buckling, the deformation shape, and the radial displacement. The methodology is set up such that the reduced-scale design parameters are determined sequentially. The buckling behavior of the two shells is compared using a semi-analytical approach, linear eigenvalue, and implicit dynamic finite element analyses. The eigenmode imperfection sensitivity is also evaluated. The methodology is successfully applied to isotropic, cross-ply, quasi-isotropic, and sandwich conical shells. The prediction accuracy is mainly affected by not being able to simultaneously satisfy all scaling parameters, by non-negligible flexural anisotropy and transverse shear, and by differences in imperfection sensitivity between the full-scale and reduced-scale shells. In any case, accurate results are obtained for the considered shells. The radial displacement is most difficult to predict, which is attributed to the membrane prebuckling assumption and neglecting the presence of imperfections. Finally, it is observed that larger eigenmode imperfections affect the accuracy, but they do not cause the methodology to fail. For future work, it is recommended to validate the methodology through experimental testing.

This thesis aims to provide more insight into the crashworthiness behaviour of a composite aircraft fuselage. This is achieved by studying both analytically and numerically the crushing behaviour of composite energy absorbers. The analytical model, which is based on energy dissipation rates, generalises previously derived analytical models to study a wider variety of structural components. It is found that the analytical model can give an estimation of the mean crushing load of square, C-shaped, and corrugated specimens. Next to the analytical study, numerical models are created to gain insight into the behaviour of the studied shapes when subjected to crash loading. Another part of the numerical study investigates the effectiveness of different material models to capture the complex failure mechanism of the composite crash absorber. In LS-Dyna, the MAT054, MAT058, and MAR262 material models are used, which are all able to reproduce the results obtained from reference tests. The ability of the material models to recreate the crushing phenomenon is thought to originate from the opportunity LS-Dyna offers to degrade material properties in the crush frond, which incorporates the formation of cracks and delaminations. In the final stage of this research, a case study is performed, where different absorbers are introduced into a simplified digital twin of the Next Generation Multifunctional Fuselage Demonstration as developed by the STUNNING project. Here, it is found that the fuselage's energy absorption can be increased by introducing energy absorbers. All the simulated fuselage sections show similar behaviour in the first crash stages, characterised by the flattening of the lower section of the frame. Subsequently, the behaviour of the absorber and the surrounding structure highly depends on the absorber's integration. Finally, it is found that by combining the numerical results of a baseline fuselage section and the analytical model for the energy absorbers, one can estimate the energy absorption of the augmented structure with a discrepancy of less than 20$\%$. This proves that the suggested analytical-numerical method can aid engineers during the preliminary design for crashworthiness of a thermoplastic composite aircraft fuselage.

Predicting the critical buckling load of cylindrical shells with circular cutouts subjected to uniform axial compression is an important part of the structural design in the aerospace industry as buckling significantly reduces the load-carrying capability of the structure. A cutout constitutes a major disruption in the shell geometry, and therefore it should be expected that it has a significant effect on the sustainable buckling load. An analytical solution for estimating the buckling load of isotropic and quasi-isotropic composite cylindrical shells with circular cutouts is developed to assess changes made to the geometry and the material during the preliminary design phase quickly. The Ritz method is employed to minimize the total potential energy of an ideal shell that contains a central opening in order to predict a linear buckling load. Finite element simulations are conducted to verify the accuracy of the analytical solution. In addition, they are used to investigate the evolution of buckling modes, the effects of initial geometric imperfections, as well as the shell failure mode. The nondimensional curvature parameter α can be used to categorize the buckling behavior of cylindrical shells and is a function of the cutout radius, the shell radius, and the shell thickness. A small cutout has virtually no influence on the buckling load compared to a pristine shell and the displacement pattern at buckling is global. The buckling load decreases rapidly for moderately large cutouts where the stability loss is the result of a local buckling mode that immediately leads to global buckling. Cylindrical shells with large cutouts are again relatively insensitive to an increase of the cutout size, but the buckling load is greatly reduced relative to a shell without a cutout. Large openings also feature a stable local buckling mode where substantial lateral prebuckling displacements emerge before the structure buckles globally. While the analytical procedure theoretically should not capture the onset of global buckling independent of local buckling, it follows numerical trends for cutouts of moderate and large size regardless. Therefore, it may be used during preliminary design to estimate the impact of changes made to the shell geometry and material. Local buckling is caused by high compressive stresses next to the cutout and, in some cases, large lateral prebuckling displacements. The detrimental effect of the stress field may be partially relieved in composite cylindrical shells by reducing the amount of axial bending stresses that occur. Hence, the chosen stacking sequence can have a significant influence on the buckling load of shells with moderate and large openings.

The COVID-19 pandemic had a significant impact on the aviation industry with more than 60% reduction in the passenger traffic in the year 2020 compared to year 2019. The passenger traffic which is now expected to reach the 2019 level only around the year 2024 will still continue to grow but at a lower pace compared to the pre pandemic levels. The environmental issues which were a major concern for the aviation industry before pandemic will still be relevant. The objective of achieving an environmentally friendly zero emission aircraft will not be met only with alternative fuels and propulsion concepts but also require advanced material technologies and novel designs. A promising technology having the potential to improve the performance of an aircraft by improving the structural efficiency is the application of aeroelastic tailoring with the help of composite materials. However, incorporating aeroelastic tailoring with composite materials in the design process is not a trivial task. In the traditional design process, knowledge about the design increases, while the design freedom decreases as one goes from conceptual to preliminary and finally to the detailed design. For conventional designs, the lack of knowledge during the initial stages is compensated through empirical knowledge. However, the lack of such empirical knowledge for novel design and advanced technology, results in the need for increased physics-based knowledge during the initial design process. In the research presented in this dissertation, the focus was on increasing knowledge in the early stages of the aeroelastic design process of a composite wing. As current state of art in aeroelastic tailoring does not include critical gust and fatigue loads, this thesis is focused on including critical gust loads and fatigue loading requirements in the preliminary aeroelastic optimization framework…