Integrated design of cellular composites using a level-set topology optimization method

Design of shell-infill structures by a multiscale level set topology optimization method

This thesis aims at understanding and improving topology optimization techniques focusing on density-based level-set methods and geometrical nonlinearities. Central in this work are the numerical modeling of the mechanical response of a design and the consistency of the optimization process itself. Concerning the first topic, we investigate different means to improve the robustness of density-based numerical models including geometrical nonlinearities. The conventional approach (scaling the local material properties) can result in convergence problems due to excessive deformation in low-stiffness finite elements. To avoid excessive deformation, we combine the element connectivity parameterization method (adapting the connectivity between finite elements) with level-set-based topology optimization. Furthermore, we achieve greater robustness of analysis method using a second and improved approach called element deformation scaling. This approach eliminates the need for solving internal equilibrium equations (needed for the element connectivity parameterization method) via an explicit relation between the local internal and global external displacement field. The second focus of this thesis is on the optimization process of level-set-based topology optimization, and in particular its numerical consistency. We observe that signed-distance reinitialization of the level-set function affects the shape of a design in the optimization process. To minimize this effect, we propose a discrete level-set method that is based on an approximate Heaviside function and focusing on the implementation. Furthermore, we also propose a level-set-based topology optimization method using an exact Heaviside function and mathematical programming that effectively eliminates the need for reinitialization. We demonstrate that our density-based level-set method is closely related to conventional density-based topology optimization methods, while offering the advantage of more control over the geometrical complexity. On the other hand, we confirm that the dependence of the final result on the initial design which remains one of the big challenges for level-set-based topology optimization. The potential of the proposed level-set method is shown by applying it on problems with stress constraints and geometrical nonlinearities and performing manufacturing tolerant topology optimization. Finally, this thesis offers a review of level-set methods for structural topology optimization to identify and discuss the different approaches that are available in literature. We can distinguish between level-set methods by examination of their design parameterization, sensitivities, update procedures and regularization techniques. A level-set-based design parameterization offers the advantage of a crisp distinction between subdomains. For this reason, X-FEM approaches and conforming discretizations are an interesting option to retain the crisp nature of the level-set-based description of the design. Many level-set methods are combined with density-based numerical models and are, therefore, closely related to conventional density-based topology optimization methods. In particular, recently proposed projection methods have much in common with a level-set-based design description. The results of level-set-based TO methods often rely heavily on regularization techniques that introduce inconsistencies in the optimization process. Numerical consistency does not necessarily lead to the best search direction, but is essential to find a Karush-Kuhn-Tucker point of the discretized optimization problem.

SummaryThe parametric level set approach is an extension of the conventional level set methods for topology optimization. By parameterizing the level set function, level set methods can be directly coupled with mathematical programming to achieve better numerical robustness and computational efficiency. Moreover, the parametric level set scheme can not only inherit the primary advantages of the conventional level set methods, such as clear boundary representation and the flexibility in handling topological changes, but also alleviate some undesired features from the conventional level set methods, such as the need for reinitialization. However, in the existing radial basis function–based parametric level set method, it is difficult to identify the range of the design variables. Besides, the parametric level set evolution often struggles with large fluctuations during the optimization process. Those issues cause difficulties both in numerical stability and in material property mapping. In this paper, a cardinal basis function is constructed based on the radial basis function partition of unity collocation method to parameterize the level set function. The benefit of using cardinal basis function is that the range of the design variables can now be clearly specified as the value of the level set function. A distance regularization energy functional is also introduced, aiming to maintain the desired signed distance property during the level set evolution. With this desired feature, the level set evolution is stabilized against large fluctuations. In addition, the material properties mapped from the level set function to the finite element model can be more accurate.

A new multi-material level set topology optimization method with the length scale control capability

Level-set approaches are a family of domain classification techniques that rely on defining a scalar level-set function (LSF), then carrying out the classification based on the value of the function relative to one or more thresholds. Most continuum topology optimization formulations are at heart, a classification problem of the design domain into structural materials and void. As such, level-set approaches are gaining acceptance and popularity in structural topology optimization. In conventional level set approaches, finding an optimum LSF involves solution of a Hamilton-Jacobi system of partial differential equations with a large number of degrees of freedom, which in turn, cannot be accomplished without gradients information of the objective being optimized. A new approach is proposed in this paper where design variables are defined as the values of the LSF at knot points, then a Kriging model is used sto interpolate the LSF values within the rest of the domain so that classification into material or void can be performed. Perceived advantages of the Kriging-interpolated level-set (KLS) approach include alleviating the need for gradients of objectives and constraints, while maintaining a reasonable number of design variables that is independent from the mesh size. A hybrid genetic algorithm (GA) is then used for solving the optimization problem(s). An example problem of a short cantilever is studied under various settings of the KLS parameters in order to infer the best practice recommendations for tuning the approach. Capabilities of the approach are then further demonstrated by exploring its performance on several test problems.

This research is motivated by the need for control of flexural vibrations of lightweight plates. It addresses application of the level set method in optimal control of vibrations in plate-like structures. One of the most commonly practiced methods in control of vibration is to apply constrained layer damping patches to the surface of a structure. In order to consider the weight efficiency of the structure, the best shape and locations of the patches should be determined to achieve the optimum vibration suppression with lowest amount of damping patch. A novel topology optimization approach is proposed that is capable of finding the optimum shape and locations of the patches simultaneously. A 2D cantilever plate, undergoing flexural vibrations, will be considered. The optimal damping set will be found in the structure, such that the lowest modal energy in the fundamental vibration mode of the system is achieved. The proposed level set topology optimization method shows capability of determining the optimum damping set in structures accurately.

A topology optimization method for phononic crystals is developed for the design of sound barriers, using the level set approach. Given a frequency and an incident wave to the phononic crystals, an optimal shape of periodic inclusions is found by minimizing the norm of transmittance. In a sound field including scattering bodies, an acoustic wave can be refracted on the obstacle boundaries, which enables to control acoustic performance by taking the shape of inclusions as the design variables. In this research, we consider a layered structure which is composed of inclusions arranged periodically in horizontal direction while finite inclusions are distributed in vertical direction. Due to the periodicity of inclusions, a unit cell can be considered to analyze the wave propagation together with proper boundary conditions which are imposed on the left and right edges of the unit cell using the Bloch theorem. The boundary conditions for the lower and the upper boundaries of unit cell are described by impedance matrices, which represent the transmission of waves between the layered structure and the semi-infinite external media. A level set method is employed to describe the topology and the shape of inclusions. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. Through several numerical examples, the applicability of the proposed method is demonstrated.

Structural topology optimization.- On basic properties of Michell's structures.- Validation of numerical method by analytical benchmarks and verification of exact solutions by numerical methods.- Introduction to shape and topology optimization.- Homogenization method for shape and topology optimization.- Level set method for shape and topology optimization.- Compliance minimization of two-material elastic structures.- Some notes on topology optimization of vibrating continuum structures.- Topology optimization of vibrating continuum structures for maximum values of simple and multiple eigenfrequencies and frequency gaps.- On optimum shape design and periodicity of band-gape beam structures.- Topological design for minimum dynamic compliance of continuum structures subjected to forced vibration.- Topological design for minimum sound emission from bi-material structures subjected to forced vibration.- Discrete material optimization of vibrating laminated composite plates for minimum sound emission.- Optimization of diffusive transport problems.- Fluid topology optimization: stokes and Navier-stokes models.- Topology optimization of coupled multi-physics problems.- Topology optimization based on the extended finite element method.- Topology optimization of meso- and nano-scale problems.- Topology optimization under uncertainty.- A brief review of numerical methods of structural topology optimization.- The free material design in planar elasticity.

In recent years, additive manufacturing has been gaining ground among traditional manufacturing methods in many different applications. This is mainly due to its ease of prototyping, material savings and above all the highly complex geometries that can be achieved. Despite its many advantages, the use of additive manufacturing in electrical machines has been limited. The lack of maturity of the manufacturing processes and the increase in losses when applying additive approaches to some parts, mainly stators, have led little implementation in the field of electrical machines. The geometric freedom offered by additive manufacturing paves the way for the use of novel optimisation techniques, such as topology optimisation. Traditionally, topology optimisation has not been widely used because the complex geometries obtained by these algorithms are not easily produced by traditional manufacturing methods. However, thanks to additive manufacturing, topology optimisation has started to be used, mainly for mechanical applications and it has demonstrated its ability to reduce weight without compromising mechanical stability in several cases. Nevertheless, topology optimisation for physics other than mechanical has not been thoroughly studied, let alone for the simultaneous optimisation of different physics. In view of the scenario described, this thesis investigates the benefits that additive manufacturing and topology optimisation can bring to the design of electrical machines. For the additive manufacturing of electrical machines, a thorough literature review is carried out on the application of this manufacturing approach to different parts of electrical machines, active and non-active. Two main case studies are presented: in the first, the rotor, the shaft, the stator and the housing of an aerospace actuator are manufactured by L-PBF in FeCoV and the assembled actuator is tested. In the second case, a rotor is manufactured in FeSi by LP- DED. Finally, not considered an additive manufacturing case study itself in this thesis, but two electrical conductor prototypes are manufactured in CuCr1Zr and an additional geometry is manufactured in AlSi10Mg by L-PBF. With regard to topology optimisation of electrical machines, a description of the main topology optimisation methods used is given and their application to electrical machine components is presented. In a similar way to additive manufacturing, two case studies are analysed. First, a novel multiphysics - mechanical and electromagnetic - topology optimisation method is presented, in which the rotor and the stator of a permanent magnet motor are optimised simultaneously. This method is compared with two methods found in the literature, SIMP and on-off. Secondly, another topology optimisation method is described and applied to an electrical conductor model. The proposed algorithm involves a multiphysics - thermal and electromagnetic - hybrid parametric topology optimisation approach. Two reduced length geometries and one full length prototype are built via additive manufacturing. Finally, conclusions regarding additive manufacturing and topology optimisation of electrical machines are drawn from the work presented in this thesis. Future lines of work for additive manufacturing and topology optimisation of electrical machines are also presented.

In conventional level set methods, the slope of the level set function needs to be well controlled to maintain the numerical stability during the topology optimization process. One common solution is to regularize the level set function to be a signed distance function, which is usually achieved by periodically implementing the so called re-initialization scheme to force the level set function to gain the desired signed distance property. However, the re-initialization scheme will bring some unwanted drawbacks to the optimization process, such as zero level set drifting, time consuming etc. In addition, re-initialization is usually implemented outside the optimization loop, which will cause convergence issues. In this paper, a distance regularization functional is introduced to the structural topology optimization objective functional to ensure the signed distance property of the level set function near the structure boundaries. This functional can also keep the level set function to be constant-value at positions far away from the structural boundaries. The radial basis function (RBF) based parameterization technique together with the mathematical programming are utilized to improve the potential capability of handling multiple constraints for the topology optimization. The combination of these two techniques makes the level set based topology optimization be capable of handling complicated multi-constrained problems with higher numerical efficiency, leaving no compromise to multiple drawbacks. To demonstrate the validity of the proposed scheme, benchmark examples on minimum compliance structural optimization are employed. This type of problem is computed by the conventional level set method with the introduced distance regularization functional, the RBF based parametric level set and at last, the distance regularized RBF based parametric level set separately to demonstrate their differences.

Topological shape optimization of geometrically nonlinear structures using level set method

Topology optimization of natural convection using porous metal foam based on the adjoint lattice Boltzmann method and level set method

The aim of this paper is to provide a survey on the recent development in level set methods in inverse problems and optimal design. We give introductions on the general features of such problems involving geometries and on the general framework of the level set method. In subsequent parts we discuss shape sensitivity analysis and its relation to level set methods, various approaches on constructing optimization algorithms based on the level set approach, and special tools needed for the application of level set based optimization methods to ill-posed problems. Furthermore, we provide a review on numerical methods important in this context, and give an overview of applications treated with level set methods. Finally, we provide a discussion of the most challenging and interesting open problems in this field, that might be of interest for scientists who plan to start future research in this field.

This paper deals with simultaneous topology and shape optimization of elastic contact problems. The structural optimization problem for an elastic contact problem is formulated. Shape as well as topological derivatives formulae of the cost functional are provided using material derivative and asymptotic expansion methods, respectively. These derivatives are employed to formulate necessary optimality condition for simultaneous shape and topology optimization and to calculate a descent direction in numerical algorithm. Level set based numerical algorithm for the solution of this optimization problem is proposed. Numerical examples are provided and discussed.

The numerical method of stamp topological optimization taking into account fatigue strength is presented in the work. It is proposed to take into account the restrictions on the stress state in accordance with the curve of the dependence of the maximum stresses on the number of loading cycles in the ESO topological optimization method. An approach to the selection of the evolutionary coefficient with a step-by-step increase in the rejection coefficient is proposed when constructing an iterative scheme for the rejection of elements by the method of topological optimization. The calculation of the stamp optimal topology with a decrease in volume due to the removal and redistribution of material was carried out in the study. The new geometric model of the optimal topology stamp is based on the predicted distribution of elements with a minimum stress level. The verification calculation of the stress state of the stamp of optimal topology with an assessment of fatigue strength was carried out in the work. The numerical calculation was carried out using the finite element method in the Ansys software package. The minimized stamp volume decreased by 35% according to the calculation results. The results of the study can be further applied in the development of topological optimization methods and in the design of stamping tools of optimal topology.