A Multi-Mesh Strategy for Continuum Topology Optimization under Correlated Uncertainties

A new fuzzy strategy for size and topology optimization of truss structures

A globally convergent method to accelerate topology optimization using on-the-fly model reduction

The center of mass and the inertia tensor are important parameters characterizing the kinetic properties of many structural systems, e.g., the rotor, spacecraft, and aircraft. In these occasions, designs not consistent with the requirements on the center of mass and inertia tensor cannot be accepted. In this paper, mathematical models are established to carry out the topology optimization of continua considering the mass and inertia characteristics. Two numerical examples, one for minimizing the static compliance considering locations of the center of mass and the other for maximizing the fundamental eigenfrequency considering both locations of the center of mass and values of the inertia products, are presented to demonstrate that the requirements on mass and inertia characteristics could have a profound influence on the final optimized layouts in structural topology optimization problems. On the other hand, the mass and inertia characteristics can also be tuned by topology optimization methods. In the second numerical problem, the smooth minimum operator is used to approximate the fundamental eigenfrequency with the first 4 eigenfrequencies so that the possible non-differentiability of the repeated eigenfrequencies can be avoided. Due to its extreme importance in engineering practice, the work of this paper can be useful to carry out designs in accordance with the requirements on the mass and inertia characteristics under the framework of traditional density-based structural topology optimization methods.

This paper presents an efficient and comprehensive MATLAB code to solve two-dimensional structural topology optimization problems, including minimum mean compliance, compliant mechanism synthesis and multi-load compliance problems. The Unsmooth Variational Topology Optimization (UNVARTOP) method, developed by the authors in a previous work, is used in the topology optimization code, based on the finite element method (FEM), to compute the sensitivity and update the topology. The paper also includes instructions to improve the bisection algorithm, modify the computation of the Lagrangian multiplier by using an Augmented Lagrangian to impose the constraint, implement heat conduction problems and extend the code to three-dimensional topology optimization problems. The code, intended for students and newcomers in topology optimization, is included as an appendix (Appendix A) and it can be downloaded from https://github.com/DanielYago together with supplementary material.

In this paper, topology optimization (TO) problems of structures subjected to design-dependent loads are formulated in the context of a special density-based TO approach wherein a Non Uniform Rational Basis Spline (NURBS) hyper-surface is used to represent the topological descriptor. Unlike classical density-based TO approaches, the NURBS-density-based method allows representing the pseudo-density field through a purely geometric computer-aided design compatible entity. In this context, the problem is formulated in the most general case of inhomogeneous Neumann-Dirichlet boundary conditions. Moreover, a thorough study of the penalty function of the design-dependent loads is carried out to investigate its effect on the optimized topologies and overcome the singularity effect related to the zones characterized by low values of the pseudo-density field. Finally, a wide campaign of sensitivity analyses is conducted to investigate the influence of the integer parameters of the NURBS entity, of the combination of design-depedent loads and inhomogeneous Neumann-Dirichlet boundary conditions, and of the concentrated load on the optimized topology. The effectiveness of the approach is tested on 2D and 3D benchmark problems.

With the rapid advancement of computing technology, mathematical modeling and computer simulations have become common and standard methods for design of complex systems. To exercise these models intelligently and eliminate the burden of manual iteration, manipulating inputs and reviewing outputs, optimization strategies are applied in a simulation-based design environment. In reality, the knowledge that engineers have about the design problem is imperfect and incomplete, uncertainty always exists in simulation-based design. The conventional deterministic optimization techniques do not consider the impact of such uncertainties. If the system responses driving the design decisions are very sensitive to these uncertainties, serious problems might confront the designer. Performance of the actual system may differ greatly from expected optimum performance. Again, if the design solution is near one or more constraint boundaries, then even slight uncertainties or changes in the operating environment could produce failed, unsafe designs, and/or result in substantial performance degradation. In recent years, the concept of reliability and robust design has become very popular. Probabilistic design analysis and statistical modeling allow the identification of designs that qualify as not only feasible, but as consistently feasible in the face of uncertainty. Reliability based topology optimization (RBTO) extends the reliability notion to the area of structural topology optimization. Although reliability based design optimization (RBDO) is a well established research area, very few work has been presented in the area of RBTO. In this study, we give an overview of the various works done in the field of RBTO. Ever since Bendsoe and Kikuchi introduced the topology optimization using a homogenization method, many topology optimization methods have been developed for the both linear and nonlinear structures. Traditional topology optimization methods drive the topology of a structure to an optimum design based on the constraint on mass. Kharmanda et al. introduced reliability constraints in deterministic topology optimization problem. They proposed a heuristic strategy that aims to reduce mass while improving the reliability level of the structure without greatly increasing its weight. But the limit state function used by them was not based on failure criteria for the structure. Their formulation considered uncertainty with respect to ∗Design Automation Laboratory †Advisor geometrical dimension and applied load only. Also their reliability analysis seems to be independent of the boundary and loading condition, so their results showed similar values for the uncertain variables for different structures. Tovar et al. have developed the Hybrid Cellular Automaton (HCA) method for structural synthesis of continuum material where the state of each cell is defined by both density and strain energy. In Agarwal, a decoupled RBDO approach is employed such that the topology optimization is separate from the reliability analysis. Patel et al. showed the use of RBTO using the gradient free Hybrid Cellular Automata (HCA) method. Their formulation incorporates uncertainty with respect to material property also. They considered limit state function based on failure modes on the output displacements. This study aims at giving a comparative study on these methods. It points out the advantages and disadvantages of the methods, and their significance in the field of RBTO. As large numbers of design variables are associated with continuum topology optimization problems, RBTO methods are inherently computationally expensive because of additional required system analysis associate with RBDO. Based on the observations, this study presents an efficient method for incorporating reliability analysis in topology optimization.

Robust topology optimization of multi-material structures under load uncertainty using the alternating active-phase method

The increasing structural flexibility of large aircraft leads to significant aeroelastic effects. More efficient topology optimization techniques are required for the design to further take advantage of aeroelasticity and obtain lightweight structures. This paper proposes a moving boundary meshfree topology optimization that combines the Galerkin method of weighted residuals and non-uniform rational B-splines (NURBS). The solution domain is described by the control points of NURBS and its property is calculated adaptively with an integration subtraction technique. The minimal compliance is searched for using the globally convergent method of moving asymptotes (GCMMA) by designing the locations of control points as subject to volume and flux constraints. The method is first applied to a typical two-dimensional design example with symmetric boundary conditions. The results show that the shape constraints can be conveniently applied, and smoother boundaries are obtained with fewer parameters. Then, a three-dimensional wing structure with asymmetric boundary conditions is optimized. A three-dimensional flight load that combines the high-order-panel and meshfree methods is employed to calculate the elastic loads and update asymmetric external loads during the optimization process. The designed wing satisfies engineering requirements and the presented method can solve the practical topology optimization problems of three-dimensional structures.

We generalize stochastic mathematical programs with equilibrium constraints (SMPEC) introduced by Patriksson and Wynter (Ref. 1) to allow for the inclusion of joint upper-level constraints and to change the continuity assumptions with respect to the uncertainty parameters assumed before by measurability assumptions. For this problem, we prove the existence of solutions. We discuss also algorithmic aspects of the problem, in particular the construction of an inexact penalty function for the SMPEC problem. We apply the theory to the problem of structural topology optimization.

A method is developed to systematically remove and reintroduce low density elements from and into the finite element mesh on which the structural topology optimization problem is defined. The material density field which defines the topology and the local ‘stiffness’ of the structure is optimally distributed via non‐linear programming techniques. To prevent elements from having zero stiffness, an arbitrarily small lower bound on the material density is typically imposed to ensure that the global stiffness matrix does not become singular. While this approach works well for most minimum compliance problems, the presence of low density elements can cause computational problems, particularly in structures that exhibit geometric non‐linearities, e.g. in compliant mechanisms. To resolve this problem, a systematic approach for removing and reintroducing low density elements is presented, and the substantial performance improvements both in design and computational efficiency of the method over current methods are discussed. Several structures and compliant mechanisms are designed to demonstrate the method. Copyright © 2003 John Wiley & Sons, Ltd.

In this paper, design-dependent pressure loads are addressed for topology optimization problems of multiphase material structures. The modeling of unidirectional load distribution with respect to the material density is studied and mathematical formulations of exponential form are established for the first time to model unidirectional pressure loads on the moving surface and contact surface of fixed shape, respectively. In both load cases, sensitivity analysis scheme of the structural compliance is developed to show the effects of the design-dependent loads. By means of 2D and 3D numerical examples, comparative studies are performed in detail to show the effects of different load cases upon the optimal configurations. It is seen that the proposed design procedure is reliable to deal with the design-dependent load cases.

An approach for the concurrent homogenization-based microstructure type and topology optimization problem

A robust topological derivative-based multi-material optimization approach: Optimality condition and computational algorithm

Isogeometric analysis (IGA) connects CAE with standard CAD design tools and has been utilized for a decade in structural optimization. FEA has a number of drawbacks, including a lack of linkage between analytic and geometric models, a shortage of continuity between components, and a lack of accuracy for high-order elements. Isogeometric analysis effectively overcomes these obstacles. This paper proposes the isogeometric topological optimization mathematical model depending on the asymmetric constitutive relationship of tension and compression microscopic materials. We provide a new density-based topology optimization formulation in which the design domain is limited to the B-spline space. A rectangular domain with a tensor-product is embedded in a design domain that is arbitrarily shaped. The density field is represented using B-splines. The intrinsic filter provided by the B-spline representation for topology optimization is governed by the degrees of B-spline and the knot spans. This B-spline effectively removes numerical artifacts and controls the minimum feature-length in; when the B-spline basis functions span many analysis elements. Benchmark problems in structural topology optimization are computed through Evolutionary Structural Optimization (ESO) algorithm. Comparisons of the optimal structure topologies using finite element analysis and isogeometric collocation methods are analyzed. The influence of different tension-compression stiffness ratios on the final structure topologies is listed and analyzed. The numerical results show the reliability and performance of the proposed mathematical formulation and asymmetry of natural structure.

Robust topology and discrete fiber orientation optimization under principal material uncertainty