Matthijs Langelaar
Please Note
56 records found
1
Combining an Asynchronous Multi-Point and a Multi-Fidelity Infill Strategy for Surrogate-Based Optimisation
An Empirical Evaluation on Unconstrained Problems
This thesis introduces and evaluates an asynchronous MP MF infill strategy, benchmarked against eleven unconstrained MF numerical problems, using expected runtime (ERT).
The asynchronous MP single-fidelity strategy (16 ranks) reduced the geometric mean ERT by 72.7% compared with the baseline Efficient Global Optimisation with Expected Improvement. The single-point MF strategy, on the other hand, increased the ERT by 21.2%, degrading performance. The combined asynchronous MP MF strategy achieved a 68.2% reduction relative to the baseline.
These results show that the asynchronous MP strategy substantially improves performance. In contrast, the selected MF strategy proves detrimental, indicating that a revised MF strategy is required to yield further gains. ...
This thesis introduces and evaluates an asynchronous MP MF infill strategy, benchmarked against eleven unconstrained MF numerical problems, using expected runtime (ERT).
The asynchronous MP single-fidelity strategy (16 ranks) reduced the geometric mean ERT by 72.7% compared with the baseline Efficient Global Optimisation with Expected Improvement. The single-point MF strategy, on the other hand, increased the ERT by 21.2%, degrading performance. The combined asynchronous MP MF strategy achieved a 68.2% reduction relative to the baseline.
These results show that the asynchronous MP strategy substantially improves performance. In contrast, the selected MF strategy proves detrimental, indicating that a revised MF strategy is required to yield further gains.
Multidisciplinary Optimisation of Formula 1 Suspension Systems
A Conceptual Design Framework Integrating Kinematics, Structural Optimisation & Evaluation, Aerodynamics, and Vehicle Dynamics
The framework operates over a 49-dimensional design space of suspension pickup point coordinates, arm stiffness parameters, and a discrete architecture selector. A production run of 482 evaluations was completed in approximately 50 hours of simulation time, with 89.4% satisfying all 52 constraints. The resulting 10-solution, three-objective Pareto front reveals the trade-off structure between three performance objectives: tyre grip utilisation (exploitation of tyre load potential through camber and toe control), transient vehicle stability, and aerodynamic downforce. Two distinct solution families emerge: a kinematic-priority cluster that simultaneously optimises grip utilisation and vehicle stability, and an aerodynamic-priority cluster that prioritises downforce at the cost of kinematic precision. The dominant conflict is between kinematic performance and aerodynamic compliance; grip and stability are broadly complementary rather than conflicting. This is consistent with established suspension design theory and constitutes the primary design insight of the framework.
Structural optimisation is integrated via a Ground Structure Method (GSM) for upright topology optimisation, maintaining computational tractability by avoiding full finite element mesh regeneration at each iteration. The resulting structural superelements are directly coupled to the half-car compliance model, propagating stiffness information into the vehicle dynamics simulation. Domain expert knowledge is embedded through heuristics-based design partitioning and experience-based aerodynamic assessment, ensuring the optimiser explores only viable architectural configurations.
A Model-Based Systems Engineering (MBSE) development process is proposed in which the MDO execution can be structured and implemented in future work, providing a concrete foundation for transitioning the framework into a fully model-driven engineering environment.
The framework demonstrates that a 49-dimensional multidisciplinary suspension design problem is tractable as a single automated optimisation loop, suggesting that the conventional sequential, discipline-by-discipline approach at the conceptual phase is not the only viable option for problems of this complexity. The 50-hour unattended run completes a multidisciplinary concept exploration that would require multiple coordinated simulation campaigns under conventional manual workflows. ...
The framework operates over a 49-dimensional design space of suspension pickup point coordinates, arm stiffness parameters, and a discrete architecture selector. A production run of 482 evaluations was completed in approximately 50 hours of simulation time, with 89.4% satisfying all 52 constraints. The resulting 10-solution, three-objective Pareto front reveals the trade-off structure between three performance objectives: tyre grip utilisation (exploitation of tyre load potential through camber and toe control), transient vehicle stability, and aerodynamic downforce. Two distinct solution families emerge: a kinematic-priority cluster that simultaneously optimises grip utilisation and vehicle stability, and an aerodynamic-priority cluster that prioritises downforce at the cost of kinematic precision. The dominant conflict is between kinematic performance and aerodynamic compliance; grip and stability are broadly complementary rather than conflicting. This is consistent with established suspension design theory and constitutes the primary design insight of the framework.
Structural optimisation is integrated via a Ground Structure Method (GSM) for upright topology optimisation, maintaining computational tractability by avoiding full finite element mesh regeneration at each iteration. The resulting structural superelements are directly coupled to the half-car compliance model, propagating stiffness information into the vehicle dynamics simulation. Domain expert knowledge is embedded through heuristics-based design partitioning and experience-based aerodynamic assessment, ensuring the optimiser explores only viable architectural configurations.
A Model-Based Systems Engineering (MBSE) development process is proposed in which the MDO execution can be structured and implemented in future work, providing a concrete foundation for transitioning the framework into a fully model-driven engineering environment.
The framework demonstrates that a 49-dimensional multidisciplinary suspension design problem is tractable as a single automated optimisation loop, suggesting that the conventional sequential, discipline-by-discipline approach at the conceptual phase is not the only viable option for problems of this complexity. The 50-hour unattended run completes a multidisciplinary concept exploration that would require multiple coordinated simulation campaigns under conventional manual workflows.
Design, topology optimization, fabrication and testing of an adjustable compliant slit mechanism
Towards the future of X-ray astronomy
Vibroacoustic Optimisation of 3D Printed Loudspeaker Cabinets
Optimising Acoustic Performance Through Hybrid FEM/BEM Simulation of Structural Cabinet Dynamics and Sound Radiation
...
...
...
The current state-of-the-art in Topology Optimisation (TO) for design-dependent pressure-actuated CMs (PACMs) relies heavily on linear models. The determination of design-dependent pressure loads involves employing the Darcy method, which integrates Darcy's law with the drainage term to obtain the pressure field. Subsequently, the finite element method (FEM) is used to transform the pressure field into consistent nodal forces. However, it is crucial to acknowledge that these linear models are only valid for small displacements.
This thesis introduces a novel approach by incorporating nonlinearities into the solid mechanics of the TO process for PACMs in conjunction with the Darcy method. Additionally, this work incorporates nonlinearities into the solid mechanics of the TO process for PA multi-material compliant mechanisms, presenting another novel method.
Four nonlinearities in the solid mechanics of PA soft robots may occur, two of which are addressed in this thesis: geometric nonlinearities and a hyperelastic material model. Geometric nonlinearities arise from large deformations caused by high applied pressures. The Neo-Hookean material model is implemented to describe the low-stiffness material accurately.
The TO of pressure-actuated (PA) soft robots is simulated using COMSOL, a commercial software program for multi-physics simulation. This research presents a detailed comparison between theoretical predictions and practical outcomes as realised in COMSOL. Furthermore, this thesis includes a case study validating the successful implementation of the new method, covering a PA inverter, a PA compliant gripper, a PA member of the Pneumatic Networks, and a PA multi-material compliant gripper. The obtained results indicate limitations on the allowable applied pressure loads for the mechanisms, specifically in the case of the PA member of the Pneumatic Networks and a PA multi-material-compliant gripper. However, the PA inverter and PA compliant gripper validate the expectation that incorporating a hyperelastic material model yields significantly different results than the linear elastic material model. Moreover, the TO with the hyperelastic material model can predict displacements more accurately than the linear TO, as the differences between the displacements obtained from the TO and the analysis align more closely.
The Wang method is investigated to observe its influence on the range of the applied pressure loads during the TO of PA soft robots. The Wang method employs an interpolation technique that interpolates between linear and nonlinear theories. In this approach, void elements are described using linear theory, while solid elements are characterised by nonlinear theory. This interpolation method is developed to address distorted elements during large displacements. It effectively extended the range of applied loads during the TO of structures. However, it is found that this method does not influence the range of the applied load during the TO of CMs. ...
The current state-of-the-art in Topology Optimisation (TO) for design-dependent pressure-actuated CMs (PACMs) relies heavily on linear models. The determination of design-dependent pressure loads involves employing the Darcy method, which integrates Darcy's law with the drainage term to obtain the pressure field. Subsequently, the finite element method (FEM) is used to transform the pressure field into consistent nodal forces. However, it is crucial to acknowledge that these linear models are only valid for small displacements.
This thesis introduces a novel approach by incorporating nonlinearities into the solid mechanics of the TO process for PACMs in conjunction with the Darcy method. Additionally, this work incorporates nonlinearities into the solid mechanics of the TO process for PA multi-material compliant mechanisms, presenting another novel method.
Four nonlinearities in the solid mechanics of PA soft robots may occur, two of which are addressed in this thesis: geometric nonlinearities and a hyperelastic material model. Geometric nonlinearities arise from large deformations caused by high applied pressures. The Neo-Hookean material model is implemented to describe the low-stiffness material accurately.
The TO of pressure-actuated (PA) soft robots is simulated using COMSOL, a commercial software program for multi-physics simulation. This research presents a detailed comparison between theoretical predictions and practical outcomes as realised in COMSOL. Furthermore, this thesis includes a case study validating the successful implementation of the new method, covering a PA inverter, a PA compliant gripper, a PA member of the Pneumatic Networks, and a PA multi-material compliant gripper. The obtained results indicate limitations on the allowable applied pressure loads for the mechanisms, specifically in the case of the PA member of the Pneumatic Networks and a PA multi-material-compliant gripper. However, the PA inverter and PA compliant gripper validate the expectation that incorporating a hyperelastic material model yields significantly different results than the linear elastic material model. Moreover, the TO with the hyperelastic material model can predict displacements more accurately than the linear TO, as the differences between the displacements obtained from the TO and the analysis align more closely.
The Wang method is investigated to observe its influence on the range of the applied pressure loads during the TO of PA soft robots. The Wang method employs an interpolation technique that interpolates between linear and nonlinear theories. In this approach, void elements are described using linear theory, while solid elements are characterised by nonlinear theory. This interpolation method is developed to address distorted elements during large displacements. It effectively extended the range of applied loads during the TO of structures. However, it is found that this method does not influence the range of the applied load during the TO of CMs.
Unravelling gravel
Including stochastic behaviour of granular materials in design of bulk handling equipment
The structure and the position of the cut line and the connectors are optimized using the Method of Moving Asymptotes (MMA) method. A gradient based sensitivity analysis is used in the MMA. Afterwards, the influence of the cut line, the connector and the voids are observed individually. After optimizing the parts individually, the full optimization was performed, where the structure, the cut lines and the connectors with the voids were optimized. Furthermore, a parameter study was done to observe their influence on the final layout. The optimizer's behaviour was observed by looking at the optimization results and the parameter study. For example, how the optimizer tends to stack some connectors together to create a member of the structure or the influence of the voids.
With the approach presented, the main idea of optimizing a structure using topology optimization and simultaneously dividing it and optimizing the connector's position is obtained. However, the optimization has some limitations, as some assumptions and design considerations are not accurate, further research is needed to get accurate results.
...
The structure and the position of the cut line and the connectors are optimized using the Method of Moving Asymptotes (MMA) method. A gradient based sensitivity analysis is used in the MMA. Afterwards, the influence of the cut line, the connector and the voids are observed individually. After optimizing the parts individually, the full optimization was performed, where the structure, the cut lines and the connectors with the voids were optimized. Furthermore, a parameter study was done to observe their influence on the final layout. The optimizer's behaviour was observed by looking at the optimization results and the parameter study. For example, how the optimizer tends to stack some connectors together to create a member of the structure or the influence of the voids.
With the approach presented, the main idea of optimizing a structure using topology optimization and simultaneously dividing it and optimizing the connector's position is obtained. However, the optimization has some limitations, as some assumptions and design considerations are not accurate, further research is needed to get accurate results.
The goal for this research was therefore to research the effect of such a Topology Optimization based Generative Design approach on the design performance and experience. In order to do so, a robust and user-friendly TOP-GD tool was created. In this tool, multiple design solutions are explored quickly by implementing a batch-run setup that varies several chosen parameters, without needing to manually run several optimizations consecutively. Calculations are done with a simple TO script using coarse geometries, and without taking into account manufacturing methods yet. This asks for less demanding, detailed and complicated calculations than AI-based Generative Design tools currently offer, while at the same time moving from a single TO result to generating a range of candidate solutions. A lot of effort was put in the user-friendliness of the TOP-GD tool, enabling an easy workflow for the setup of design problems and a clear presentation of the results by means of a simple GUI.
The use of the TOP-GD tool in the design process was evaluated in an experiment, where it was compared with a more simple TO tool and a basic manual design approach using just pen and paper. This was done by giving the participants of the experiments three simple design assignments, that they had to carry out using each of the design approaches one by one. Evaluation of the approaches was done by comparing the design performance, and assessing the design experience with a survey and using Eye-tracking techniques.
The results of this experiment did not show enough evidence to conclude that the different design approaches had an effect on the design performance for the simple assignments executed during the experiment. However, the results of the survey show a clear positive impact of both the TO tools on the design experience, compared to manually designing. Furthermore, the TOP-GD tool has the largest positive impact on the design experience and its use in the design process is considered a big improvement, especially in quickly exploring new design directions and creating overview. This confirms the expectation that a Topology Optimization based Generative Design approach has a positive effect on the early stages of the design process. The differences found with Eye-tracking between the TO tools support this, although a more extensive experiment should be done to convincingly confirm this conclusion. ...
The goal for this research was therefore to research the effect of such a Topology Optimization based Generative Design approach on the design performance and experience. In order to do so, a robust and user-friendly TOP-GD tool was created. In this tool, multiple design solutions are explored quickly by implementing a batch-run setup that varies several chosen parameters, without needing to manually run several optimizations consecutively. Calculations are done with a simple TO script using coarse geometries, and without taking into account manufacturing methods yet. This asks for less demanding, detailed and complicated calculations than AI-based Generative Design tools currently offer, while at the same time moving from a single TO result to generating a range of candidate solutions. A lot of effort was put in the user-friendliness of the TOP-GD tool, enabling an easy workflow for the setup of design problems and a clear presentation of the results by means of a simple GUI.
The use of the TOP-GD tool in the design process was evaluated in an experiment, where it was compared with a more simple TO tool and a basic manual design approach using just pen and paper. This was done by giving the participants of the experiments three simple design assignments, that they had to carry out using each of the design approaches one by one. Evaluation of the approaches was done by comparing the design performance, and assessing the design experience with a survey and using Eye-tracking techniques.
The results of this experiment did not show enough evidence to conclude that the different design approaches had an effect on the design performance for the simple assignments executed during the experiment. However, the results of the survey show a clear positive impact of both the TO tools on the design experience, compared to manually designing. Furthermore, the TOP-GD tool has the largest positive impact on the design experience and its use in the design process is considered a big improvement, especially in quickly exploring new design directions and creating overview. This confirms the expectation that a Topology Optimization based Generative Design approach has a positive effect on the early stages of the design process. The differences found with Eye-tracking between the TO tools support this, although a more extensive experiment should be done to convincingly confirm this conclusion.
The company Delft Offshore Turbines (DOT) proposes to replace the conventional drivetrain in the top of the turbine with a hydraulic drivetrain, which has a lower mass-to-torque ratio. Although this concept shows promise, finding a low mass design that fulfils the infinite fatigue life requirement can be challenging. Using topology optimization to minimize mass while constraining the structural requirements could, therefore, be instrumental in the realization of this concept.
The design case by DOT can be classified as rotating machinery. In general rotating machinery is
commonly subjected to periodic loading that varies non-proportionally in time. This results in fluctuating stresses causing material fatigue. A consequence of the non-proportionality of loading is that the time response needs to be computed to evaluate for fatigue, which adds additional computation cost. However, another common aspect found in rotating machinery is that parts are cyclic symmetric, which allows for a potential reduction in computation cost.
In this thesis a method is presented to implement infinite fatigue life constraints into density based topology optimization for structures subjected to non-proportional loading, while cyclic symmetric properties are exploited to reduce computation cost. It was found that when the load case on a cyclic symmetric part adheres to certain conditions, a single static FE-analysis can provide multiple time steps for a quasi-static analysis. Decreasing the computational burden roughly proportional to the unique number of time steps obtained. The largest local variations in stress are estimated using a smooth min/max function and aggregated into a global constraint.
The method was tested on several numerical problems as well as applied to the DOT design case. The results showed that the method was able to properly constrain a global fatigue constraint while minimizing mass, achieving final designs that might not be trivial to find by hand.
...
The company Delft Offshore Turbines (DOT) proposes to replace the conventional drivetrain in the top of the turbine with a hydraulic drivetrain, which has a lower mass-to-torque ratio. Although this concept shows promise, finding a low mass design that fulfils the infinite fatigue life requirement can be challenging. Using topology optimization to minimize mass while constraining the structural requirements could, therefore, be instrumental in the realization of this concept.
The design case by DOT can be classified as rotating machinery. In general rotating machinery is
commonly subjected to periodic loading that varies non-proportionally in time. This results in fluctuating stresses causing material fatigue. A consequence of the non-proportionality of loading is that the time response needs to be computed to evaluate for fatigue, which adds additional computation cost. However, another common aspect found in rotating machinery is that parts are cyclic symmetric, which allows for a potential reduction in computation cost.
In this thesis a method is presented to implement infinite fatigue life constraints into density based topology optimization for structures subjected to non-proportional loading, while cyclic symmetric properties are exploited to reduce computation cost. It was found that when the load case on a cyclic symmetric part adheres to certain conditions, a single static FE-analysis can provide multiple time steps for a quasi-static analysis. Decreasing the computational burden roughly proportional to the unique number of time steps obtained. The largest local variations in stress are estimated using a smooth min/max function and aggregated into a global constraint.
The method was tested on several numerical problems as well as applied to the DOT design case. The results showed that the method was able to properly constrain a global fatigue constraint while minimizing mass, achieving final designs that might not be trivial to find by hand.
Silicon based power semiconductors have long been used as the standard in ‘semiconductor technology in power conversion applications’. Recent developments replaces the Silicon with Silicon Carbide as it results in superior performance of the power conversion applications. However, due to the increased performance, challenges regarding heat dissipation emerge and the lifetime of the power semiconductor packaging or power module is compromised. Since this leads to an increase power density, the cooling of the power module is becoming of more importance and the heat sink becomes an interesting component to optimize. The best performance of a heat sink can be obtained when the flow through the device is turbulent. Developing turbulent flow heat sinks by using topology optimization methods can significantly improve the cooling performance compared to the current designs. This work is thus aimed towards improving methods for topology optimization of turbulent flow cooling devices. However, this work focuses on turbulent flow topology optimization only and aims to improve the accuracy of current methods. It is important that the flow physics are accurate since the thermal energy transfer is dependent on the flow field. The current state-of-the-art method based on the 𝑘−𝜔 turbulence model developed by Dilgen et al. is investigated. A design domain is subdivided into elements since the finite element method (FEM) is used, such that an optimization algorithm is able to turn every element into either fluid or solid with the goal of finding the best performing structure. This density based approach, models the solid domain as a highly impermeably porous material. To inhibit flow in the solid domain a Darcy penalization is added to the momentum equation. Moreover, in the method by Dilgen et al. boundary conditions in the other turbulent fields are also enforced using a similar penalization approach. Weaknesses and errors in the density based method are investigated by comparing solutions to ones computed on a body fitted mesh. It has been found that the largest errors in the solution, by using the state-of-the-art method, appear at the solid/fluid interface in the design. In these regions the penalizations are not applied correctly for the desired boundary conditions. Therefore, in this work it is improved on by the enforcement of the boundary condition by using the Dilation method. The Dilation method focuses on the solid/fluid region where it shifts the boundary conditions for the specific dissipation rate (𝜔) and ensures it reaches the desired value at the solid/fluid interface. Secondly, severe flow leakage is found in the “porous” solid domains using the state-of-the-art method. Flow leakage is reduced by using an improved formulation of the maximum Darcy penalization in the solid domain. Finally, the improved approach is investigated in several topology optimization cases and compared to the state-of-the-art Dilgen method. It is shown that by using the new approach, different designs with a better accuracy can be obtained. In an extreme test case, the Dilgen method resulted in an infeasible design which disconnects the flow inlets from the outlets while the new and improved method resulted in a feasible design. ...
Silicon based power semiconductors have long been used as the standard in ‘semiconductor technology in power conversion applications’. Recent developments replaces the Silicon with Silicon Carbide as it results in superior performance of the power conversion applications. However, due to the increased performance, challenges regarding heat dissipation emerge and the lifetime of the power semiconductor packaging or power module is compromised. Since this leads to an increase power density, the cooling of the power module is becoming of more importance and the heat sink becomes an interesting component to optimize. The best performance of a heat sink can be obtained when the flow through the device is turbulent. Developing turbulent flow heat sinks by using topology optimization methods can significantly improve the cooling performance compared to the current designs. This work is thus aimed towards improving methods for topology optimization of turbulent flow cooling devices. However, this work focuses on turbulent flow topology optimization only and aims to improve the accuracy of current methods. It is important that the flow physics are accurate since the thermal energy transfer is dependent on the flow field. The current state-of-the-art method based on the 𝑘−𝜔 turbulence model developed by Dilgen et al. is investigated. A design domain is subdivided into elements since the finite element method (FEM) is used, such that an optimization algorithm is able to turn every element into either fluid or solid with the goal of finding the best performing structure. This density based approach, models the solid domain as a highly impermeably porous material. To inhibit flow in the solid domain a Darcy penalization is added to the momentum equation. Moreover, in the method by Dilgen et al. boundary conditions in the other turbulent fields are also enforced using a similar penalization approach. Weaknesses and errors in the density based method are investigated by comparing solutions to ones computed on a body fitted mesh. It has been found that the largest errors in the solution, by using the state-of-the-art method, appear at the solid/fluid interface in the design. In these regions the penalizations are not applied correctly for the desired boundary conditions. Therefore, in this work it is improved on by the enforcement of the boundary condition by using the Dilation method. The Dilation method focuses on the solid/fluid region where it shifts the boundary conditions for the specific dissipation rate (𝜔) and ensures it reaches the desired value at the solid/fluid interface. Secondly, severe flow leakage is found in the “porous” solid domains using the state-of-the-art method. Flow leakage is reduced by using an improved formulation of the maximum Darcy penalization in the solid domain. Finally, the improved approach is investigated in several topology optimization cases and compared to the state-of-the-art Dilgen method. It is shown that by using the new approach, different designs with a better accuracy can be obtained. In an extreme test case, the Dilgen method resulted in an infeasible design which disconnects the flow inlets from the outlets while the new and improved method resulted in a feasible design.
Efficient and Robust Topology Optimization of Compliant Mechanisms using Perturbed Geometrically Non-Linear Analysis
Paper Title: "Structural Topology Optimization using Bayesian-Enhanced Perturbed Non-Linear Analysis"