Fibre-reinforced laminated composites are constructed layer-by-layer, enabling ease of directional stiffness tailoring. Their vast design space is typically explored using two-steps. First, the optimum stiffness for given loads is conceptualised using continuous optimisation of lamination parameters (LPs). Then, discrete optimisation determines a fibre stacking sequence (SS) that closely matches these LPs. While fibre angles are conventionally limited to 45°multiples, finer increments (e.g., ±15°) can enable lighter structures. However, existing SS design methods do not scale well with this increased problem dimensionality. To overcome this challenge, we propose LP2SS, a novel methodology utilising fast Fourier transforms (FFT) and a branch-and-bound optimiser. By treating LPs as a signal, FFTs identify the number of fibre layers oriented at different angles, akin to estimating the magnitude of different frequencies within a signal. This fibre angle distribution guides the branch-and-bound optimiser, enabling efficient SS design with accurate LP matching, while satisfying empirical design rules. The ingenious use of FFTs is key to LP2SS's performance, achieving solutions within tenths of a second, compared to minutes required by state-of-the-art methods. Validated on established benchmarks and a newly proposed comprehensive test set, LP2SS marks a significant advancement in the optimal design of large-scale laminated composite structures.

Laser heating is the most common method in thermoplastic Automated Fiber Placement (AFP) due to its precision and speed, but it poses safety and cost challenges. The HUMM3 pulsed flashlamp offers a promising broadband, programmable, and relatively safer alternative. This study investigates the thermal response and deconsolidation behavior of unidirectional carbon fiber/LM-PAEK (CF/LM-PAEK™) composite tape during HUMM3 heating in AFP. A static setup replicating AFP conditions was used to investigate the thermal response of the composite tape as a function of the programmable parameters: voltage, pulse width and frequency. Deconsolidation of the tape under HUMM3 heating was assessed from micrographs and surface topography, quantified by thickness change, void content, waviness, and roughness. Results showed that voltage is the dominant parameter influencing the thermal response, whereas pulse width and frequency showed no significant effects when total energy input was held constant. The deconsolidation analysis revealed a strong correlation of thickness change, void content, waviness, and surface roughness with local temperature. These changes are believed to be driven by polymer softening, fiber decompaction, internal gas pressure buildup, and local thermal gradients across the tape width. Compared to laser heating, HUMM3 produces less severe deconsolidation, likely due to its broadband UV–VIS–NIR spectrum enabling partial matrix absorption, reducing temperature gradients across the tape width.

Coordinating and executing an MSc program is accompanied by various challenges. To name a few: a) students joining the program means have diverse backgrounds which means varying levels of prior knowledge, and b) in a program that offers specialization into different profiles, it is important to ensure that every student attains the required breadth in foundation as well, which the program mandates. At the Aerospace Structures and Materials (ASM) Department in the Faculty of Aerospace Engineering at TU Delft, we are developing a core curriculum that encompasses the first education quarter (8 weeks, 15 ECTS), which aims to successfully tackle these challenges in order to train up future aerospace engineers with solid technical foundation, critical thinking ability and collaborative working spirit. The module starts with a “soft-landing” for the students and is followed by three content sprints that cover the essentials in the three principal domains of ASM – Materials, Manufacturing and Modelling. Engineering and consulting assignments have been embedded that tie the sprints together and offer inspiration and motivation for the students. In parallel, students get trained in transversal skills. A student facing dashboard will provide comprehensive feedback on learning progress and outcome. The dashboard will additionally support the transition towards a portfolio-based assessment. In this paper, we outline key features of the core curriculum along with the process that was employed for its development.

This work proposes a methodology for the characterisation of complex pore features in unidirectional composite prepregs, and provides insights into the interaction between fibre architecture and pores. The method showcased allows to compare spatial distributions at a three-dimensional level, highlighting in the tape analysed a significant correspondence between regions of elevated tortuosity and increased pore fractions. Regions associated with highly tortuous meandering fibres exhibit a pronounced association with porosity located both in the bulk and at the tape surface, suggesting a strong interaction between non-collective fibre displacement and the probability of pore location. Furthermore, our study quantifies the length scale of feature propagation, shedding light on the spatial extent of microstructural pore occurrence within the composite. These findings have significant implications from a characterisation perspective to aid modelling approaches and manufacturing processes for high-performance composite prepregs tapes.

Understanding the microstructural variability in unidirectional composite prepreg tapes is relevant to investigating mechanisms of tape microstructure formation, their impact on its processability and the mechanical performance of the final composite part. It has been shown that three-dimensional microstructural variability at the single-fibre level can be resolved by X-ray microcomputed tomography (XCT). However, to define a representative microstructural fingerprint of a given tape, investigations at the required small voxel size lead to limited volumes of observation, which might not be representative. This research aims to extend these findings via a multiscale approach, considering scales of observations, from microscopic (single fibre) up to mesoscopic (dimension of tape) length scale, to generate further insight into the microstructural organisation of thermoplastic prepreg tapes. By exploring the ability of XCT imaging for carbon fibre-reinforced thermoplastic composites at different voxel sizes, the work aims to identify the limitations of the use of different scales of observations to capture features of microstructures and their propagation from micro- to mesoscale level. While structure tensor analysis appeared to correctly capture misaligned regions in XCT images with small voxel size (1/10 of the fibre diameter), the method proved ineffective for larger voxel size images (1/2 of the fibre diameter).

This study investigates a filament-wound tube model incorporating fiber undulation from the filament winding process. The model was analyzed using the finite element method in the linear regime, then compared with the shell model and radial crushing experiment. Results showed that the solid model predicts the radial compression stiffness with a higher level of accuracy than the shell model due to the inclusion of the fiber undulation feature. This model is a first step towards the development of a composite pressure vessel model where fiber undulation is more frequent, and also for predicting failure initiation and damage propagation.

The high void content in laser-assisted fiber placement (LAFP)-manufactured thermoplastic (TP) hinders industrial adoption, with tape deconsolidation being a critical yet understudied factor. To address this research gap, this study provides an in-depth investigation into the deconsolidation mechanisms of TP tapes during the LAFP heating phase. A series of comparative experiments were conducted to systematically evaluate the effects of tape residual stress states (fiber-matrix combined, fiber-dominated, and near stress-free) and heating methods (laser vs. oven heating) on deconsolidation behavior. Deformation along the width, voids, thickness variations and surface roughness were identified as key factors to characterize deconsolidation behavior and to elucidate its underlying mechanisms. The results reveal that the matrix residual stress plays a dominant role in exacerbating deformation along the width, nonuniformity in thickness and intralaminar voids. Additionally, fiber decompaction—induced by the recovery of elastic deformation—contributes to surface deformation by generating voids near the surfaces. Furthermore, laser-heated tapes exhibit more pronounced intralaminar voids and higher surface roughness than oven-heated counterparts, underscoring the influence of heating rate on the release of residual stress. This study advances the understanding of deconsolidation mechanisms during the LAFP heating phase, and provides recommendations for optimizing the manufacturing of LAFP-grade TP tapes.

The present study introduces an automated multidisciplinary optimization (MDO) workflow that, for the first time, couples an explicit dynamic bird strike analysis with a post-impact static stress check. This joint problem is solved during preliminary wing sizing by integrating batch Bayesian optimization on Kriging surrogates with a variance-based variable screening procedure. The optimization problem comprises 19 thickness design variables and two highly non-linear constraints, imposing a maximum leading edge penetration and a maximum post-impact front spar stress while minimizing wing mass. The workflow is demonstrated on a five-bay metallic wing segment, yielding a 43% weight saving over the best-performing design during initial data generation while respecting CS 25.631 crashworthiness limits. Results demonstrated substantial computational savings by variable screening and highlighted the necessity of the stress constraint, as designs satisfying only the penetration depth requirement could still experience critical post-impact stress levels.

Manufacturing variations in the automated fiber placement (AFP) process are one of the causes of gaps and overlaps. These manufacturing variations can be due to robot inaccuracy, tow lateral movement on the roller, tow width variation or tow compaction. An experimental setup was built to measure and investigate these various sources of manufacturing variations and their relative contributions to gap and overlap defects. This setup consisted of a commercial AFP head instrumented with additional sensors. Among all the measured sources of variations, lateral movement of the tow on the compaction roller was the biggest contributor to gaps and overlaps. The distributions of these sources of variations were fit with probability density functions. Random samples from these fits were used to simulate adjacent tows and predict the occurrence of gap and overlap defects. The distribution of predicted gaps correlated closely with the distribution of experimentally measured gaps. Thus, this approach of using statistical information about the sources of manufacturing variations to predict the frequency and magnitude of defects in a layup was validated.

This work studies the effect of compaction of tapes that have been heated using a laser during automated fibre placement. The deconsolidation has been shown to have a significant effect on the surface roughness, degree of effective intimate contact and void content. This work investigates whether after compaction the as-received (i.e., before heating) properties are obtained, or whether the deconsolidationcompaction cycle has an influence on the final tape quality. First, a lab test set-up is designed and manufactured to mimic the real manufacturing conditions. Next, the set-up is used to study the influence of the placement speed and pressure on the final quality. The results show that the effects of deconsolidation are mostly reversed, but the final tape still has slightly worse qualities. This effect will have to be taken into account for accurate modelling of the laser-assisted automated fibre placement process.

Double-curved composite structures that are manufactured via automated fiber placement, such as pressure vessels, can take advantage of tow steering to reduce weight. This design freedom comes with the cost of adding internal normal stresses to the tow, possibly leading to wrinkles or pull up. The present work investigates tow pull up both experimentally and analytically and details a correlation between tow pull up and the minimum critical steering radius, where the material tackiness and the modelled plate’s length are found to be the most influential parameters. The experimental determination of the material tackiness is the next step to improve the predictive capabilities of the proposed model.

In-situ Automated Fiber Placement (AFP) manufacturing of thermoplastic prepreg tapes has the potential to provide a fast and cost-effective manufacturing solution for large composite structures. However, it is prone to several defects, especially gaps and overlaps. One of the primary causes for this is the tape width deformation during placement. Current literature is not enough to understand this mechanism completely and the conventionally considered tape width deformation mechanism i.e., transverse squeeze flow has been suggested to be incorrect for AFP. Therefore, the research objective was to experimentally investigate the width deformation mechanism and the influence of processing parameters for thermoplastic prepreg tapes using humm3® as the heating device. The results show that the tape width deformation takes place both in the heating and consolidation phase of the process and involves spreading of the fiber-resin mixture. Additionally, the conformable roller has an influence on the post-processed tape cross-section profile. The surface roughness and the influence of the processing parameters indicate the role of temperature distribution as the influencing factor. Therefore, these results can be used to accurately model the tape width deformation to mitigate the problem of gaps and overlaps.

Manufacturing variations in the AFP process are one of the causes of gaps and overlaps. These manufacturing variations can be due to robot inaccuracy, tape lateral movement on the roller, tape width variation or tape compaction. These manufacturing variations result in incorrect position or incorrect geometry of the laid tape. An experimental setup was built to measure and investigate the causes of these manufacturing variations and their relative contributions to gap and overlap defects. This setup consisted of an instrumented AFP head. A laser tracker measured the achieved trajectory of the AFP head. A camera measured the lateral movement of the tape on the roller. Laser line scanners measured tape width before and after layup. Experimental results show that the 99^th percentile absolute deviations for each of the four measured sources vary from 0.119 mm to 0.534 mm as compared to the specified tape width of 6.35 mm. Among all the measured sources of variations, lateral movement of the tape on the compaction roller was the biggest contributor to gaps and overlaps.

This work aims to improve the flexural behaviour of unidirectional fibre-reinforced laminates by means of coupling an optimization procedure for quasi-isotropic configurations with the design space opened by dispersed-ply orientations. The design approach consists of finding suitable alternatives to traditional laminates (with fibre orientations limited to 0°, ±45^∘, and 90°), while maintaining their stiffness characteristics. This strategy isolates the interlaminar response as the objective function that is optimized to improve their flexural behaviour. To this end, a modified Ant Colony Optimization was implemented and geared towards optimizing the interlaminar stress profile, allowing plies at every possible 5° orientation, with the ultimate goal of delaying delamination. To validate the approach, a traditional reference laminate and derived fully dispersed designs were experimentally tested. The correlated responses show that it was not possible to improve flexural resistance. However, the typical flexural brittleness of laminates can be modified into a pseudo-ductile behaviour.

Anisotropic materials formed by living organisms possess remarkable mechanical properties due to their intricate microstructure and directional freedom. In contrast, human-made materials face challenges in achieving similar levels of directionality due to material and manufacturability constraints. To overcome these limitations, an approach using 3D printing of self-assembling thermotropic liquid crystal polymers (LCPs) is presented. Their high stiffness and strength is granted by nematic domains aligning during the extrusion process. Here, a remarkably wide range of Young's modulus from 3 to 40 GPa is obtained during by utilizing directionality of the nematic flow during the printing process. By determining a relationship between stiffness, nozzle diameter, and line width, a design space where shaping and mechanical performance can be combined is identified. The ability to print LCPs with on-the-fly width changes to accommodate arbitrary spatially varying directions is demonstrated. This unlocks the possibility to manufacture exquisite patterns inspired by fluid dynamics with steep curvature variations. Utilizing the synergy between this path-planning method and LCPs, functional objects with stiffness and curvature gradients can be 3D-printed, offering potential applications in lightweight sustainable structures embedding crack-mitigation strategies. This method also opens avenues for studying and replicating intricate patterns observed in nature, such as wood or turbulent flow using 3D printing.

For pick-and-place processes to become widely implemented in industry a consistent and acceptable product quality needs to be achieved. One important quality criterion is the fiber angle deviations in the reinforcement. Handling a reinforcement will subject it to forces due to e.g. gravity and accelerations. These forces can result in in-plane shear and subsequently in fiber angle deviations. The current work looks at predicting and preventing in-plane shear induced fiber angle deviations by studying the positioning of pick-up points. In the state of the art the positioning of individual pick-up points is typically either not discussed or is based on the mould where the fabric is to be draped on – not on the effect of the handling on the fabric. The relationship between the positioning of the pick-up points and the behavior of the fabric should however also be considered. Finite element simulations validated through experimental work will be used to study the influence of pick-up point location on the in-plane shear strain for a bi-axial Non Crimp Fabric [NCF]. In [1] tolerances have been set for the fiber angle deviations, additionally the relationship between in-plane shear strain and fiber angle deviations has been demonstrated for the specific NCF. These results are used in the present work to evaluate the results from the simulations. The current work will demonstrate that it is possible to control the in-plane shear induced fiber angle deviations by varying the position and number of the pick-up points. Additionally, It will show whether simulations for the positioning of pick-up points on large reinforcements can be simplified by looking at one instance of a repeating pattern. The paper will provide a framework for the determination of positioning of pick-up points while predicting and preventing in-plane shear induced fiber angle deviations.

Finding new ways to evaluate the variability of microstructures, and its effect on macroscopic properties such as permeability and mechanical performance [1,2] is of increasing interest in the composite field. The variability of microstructural features at a three-dimensional level is not fully understood and its effect on macroscale properties is not well established, and mostly analyzed at a phenomenological level [3]. We introduced in recent work a method based on X-ray Computed Tomography for the threedimensional reconstruction of the fibrous microstructure of unidirectional tapes at a single fibre resolution [4]. A schematic of the workflow is represented in Figure 1. Three descriptors are introduced in the work to describe increasing level of complexity in the microstructural organization, from a single fiber path level with differential tortuosity, to group behavior with collective motion, to fibre network connectivity with length of contact. These descriptors and their interdependence highlight local effects like edge-core segregation in microstructural characteristics. However, in order to achieve a more complete definition of the unidirectional tape domain, understanding of matrix-based features and its interrelation with fiber architecture descriptors is needed. In this work, we expand the methodology of Gomarasca et al. [4], to account for matrix-based phenomena such as tape boundary variability, and void formation and morphology. This will be showcased on a unidirectional composite tape including both fiber-based and matrix-based analysis. These methods enable advanced characterization and modelling of microstructural formation and evolution during composite manufacturing.

This work presents a new aeroelastic optimisation framework for the preliminary design of variable stiffness composite wing structures. The framework is constructed by sequentially and iteratively solving two sub-problems: aeroelastic tailoring and lay-up retrieval, using gradient-based algorithms with full-analytical sensitivities provided. During aeroelastic tailoring, the wing mass is minimised by optimising the lamination parameters and thickness of wing laminates together with wing jig twist distribution. The load cases cover not only static loads, but also the critical gust loads that are identified across the entire flight envelop at every iteration of optimisation. Further, a cruise shape constraint is included in addition to other aerostructural constraints, so that the optimal aircraft performance can be ensured. During lay-up retrieval, the manufacturable stacking sequence is retrieved according to the optimal lamination parameters with the consideration of minimal steering radius constraint. Moreover, to fix the possible constraint violations caused by lay-up retrieval, a correction strategy is incorporated to tighten the violated constraints for repeating aeroelastic tailoring. Finally, several case studies on the design of NASA common research model wing are carried out and investigated. The results indicate that the critical gust loads and cruise shape constraint have a large influence on the design of tow-steered composite wing structures, which therefore demonstrate the usefulness and benefits of the proposed optimisation framework.

A computationally-efficient strength optimization method tailoring novel composite laminates using lamination parameters is developed. The method adopts a global p-norm approach to aggregate local failure indices into a global failure index, based on the Tsai-Wu failure criterion. For design purposes, the novel composite laminates are characterized via lamination parameters that can subsequently be transformed into locally variable fiber orientations in an existing three-step optimization method. An elliptical formulation of the conservative failure envelope is applied to represent the Tsai-Wu criterion in terms of lamination parameters. A lamination-parameter-based two-level approximation for the global failure index is derived, which guarantees the anticipated conservativeness and convexity in a gradient-based optimization framework. Numerical results show that the computational efficiency of the proposed strength optimization method improves remarkably with a proper value of p, compared to the existing local-based min-max method. The method is also shown to be robust and generate converged optimum designs even in the presence of stress concentrations and singularities.

Despite their importance in benchmarking numerical simulations, buckling tests still feature compromises between component-level and high-fidelity large-scale tests. For example, compression-induced buckling tests cannot capture the through-thickness or span-wise stress gradients in wing skins. Consequently, the results obtained often require careful interpretation and conservative considerations before applying to a structure. Alternatively, a system-level large-scale test can be used, yet at considerably increased time and expense. There has been little progression towards capturing system-level behaviour in a simplified test. Herein, for buckling behaviour assessment, a three-point bending test is used, which is quick, simple to implement, and cost-effective compared to existing conventional methods. The proposed method relies on subjecting a panel with auxiliary stiffeners to bending to introduce compression-induced buckling in the skin. The three-point bend test is used, because it provides readily controllable loading and boundary conditions. The location of the neutral plane can be tailored via design of the stiffeners, thus allowing for control of the through-thickness stress gradient induced in the skin. This method is applicable to buckling of stiffened structures subject to bending (e.g., aircraft wingboxes). Numerical models are used to explore the limits of the proposed method and comparing it against traditional coupon and full-scale structural level tests. The test method is experimentally demonstrated for capturing the buckling behaviour of a thermoplastic composite panel made via automated fibre placement. The proposed approach is shown to reliably capture the buckling behaviour of a large-scale test using a simpler and more time and cost-efficient setup than conventional methods.