Design of periodic elastoplastic energy dissipating microstructures

The design of periodic elastoplastic microstructures for maximum energy dissipation is carried out using topology optimization. While the topology optimization of elastic microstructures has been performed in numerous studies, microstructural design considering inelastic behavior is relatively untouched due to a number of reasons which are addressed in this study. An RVE-based multiscale model is employed for computational homogenization with periodic boundary constraints, satisfying the Hill-Mandel principle. The plastic anisotropy which may be prevalent in materials fabricated through additive manufacturing processes is considered by modeling the constitutive behavior at the microscale with Hoffman plasticity. Discretization is done using enhanced assumed strain elements to avoid locking from incompressible plastic flow under plane strain conditions and a Lagrange multiplier approach is used to enforce periodic boundary constraints in the discrete system. The design problem is formulated using a density-based parameterization in conjunction with a SIMP-like material interpolation scheme. Attention is devoted to issues such as dependence on initial design and enforcement of microstructural connectivity, and a number of optimized microstructural designs are obtained under different prescribed deformation modes.

     

Topology optimization of continuum structures with local and global stress constraints

Topology structural optimization problems have been usually stated in terms of a maximum stiffness (minimum compliance) approach. The objective of this type of approach is to distribute a given amount of material in a certain domain, so that the stiffness of the resulting structure is maximized (that is, the compliance, or energy of deformation, is minimized) for a given load case. Thus, the material mass is restricted to a predefined percentage of the maximum possible mass, while no stress or displacement constraints are taken into account. This paper presents a different strategy to deal with topology optimization: a minimum weight with stress constraints Finite Element formulation for the topology optimization of continuum structures. We propose two different approaches in order to take into account stress constraints in the optimization formulation. The local approach of the stress constraints imposes stress constraints at predefined points of the domain (i.e. at the central point of each element). On the contrary, the global approach only imposes one global constraint that gathers the effect of all the local constraints by means of a certain so-called aggregation function. Finally, some application examples are solved with both formulations in order to compare the obtained solutions.     

Layout and material gradation in topology optimization of functionally graded structures: a global–local approach

By means of continuous topology optimization, this paper discusses the influence of material gradation and layout in the overall stiffness behavior of functionally graded structures. The formulation is associated to symmetry and pattern repetition constraints, including material gradation effects at both global and local levels. For instance, constraints associated with pattern repetition are applied by considering material gradation either on the global structure or locally over the specific pattern. By means of pattern repetition, we recover previous results in the literature which were obtained using homogenization and optimization of cellular materials.     

Constrained isogeometric design optimization of lattice structures on curved surfaces: computation of design velocity field

This paper presents an immersed boundary approach for level set topology optimization considering stress constraints. A constraint agglomeration technique is used to combine the local stress constraints into one global constraint. The structural response is predicted by the eXtended Finite Element Method. A Heaviside enrichment strategy is used to model strong and weak discontinuities with great ease of implementation. This work focuses on low-order finite elements, which given their simplicity are the most popular choice of interpolation for topology optimization problems. The predicted stresses strongly depend on the intersection configuration of the elements and are prone to significant errors. Robust computation of stresses, regardless of the interface position, is essential for reliable stress constraint prediction and sensitivities. This study adopts a recently proposed fictitious domain approach for penalization of displacement gradients across element faces surrounding the material interface. In addition, a novel XFEM informed stabilization scheme is proposed for robust computation of stresses. Through numerical studies the penalized spatial gradients combined with the stabilization scheme is shown to improve prediction of stresses along the material interface. The proposed approach is applied to the benchmark topology optimization problem of an L-shaped beam in two and three dimensions using material-void and material-material problem setups. Linear and hyperelastic materials are considered. The stress constraints are shown to be efficient in eliminating regions with high stress concentration in all scenarios considered.     

Microstructural topology optimization of structural-acoustic coupled systems for minimizing sound pressure level

The aim of this paper is to present a microstructural topology optimization methodology for the structural-acoustic coupled system. In the structural-acoustic system, the structure is considered to be a thin composite plate composed of periodic uniform microstructures. The discrete design variables are used in the microstructural topology optimization, and the constitutive matrix is interpolated by the power-law scheme at the micro scale. The equivalent macro material properties of the microstructure are computed through the homogenization method. The design objective is to minimize the sound pressure level (SPL) in an interior acoustic medium. The sensitivities of the SPL with respect to design variables are derived. The bi-directional evolutionary structural optimization (BESO) method is extended to solve the structural-acoustic coupled optimization problem to find the optimal material distribution of the microstructure. Numerical examples of a hexahedral box and an automobile passenger compartment are given to demonstrate the efficiency of the presented microstructural topology optimization method.     

On an alternative approach to stress constraints relaxation in topology optimization

The paper deals with the imposition of local stress constraints in topology optimization. The aim of the work is to analyze the performances of an alternative methodology to the ε-relaxation introduced in Cheng and Guo (Struct Optim 13:258–266, 1997), which handles the well-known stress singularity problem. The proposed methodology consists in introducing, in the SIMP law used to apply stress constraints, suitable penalty exponents that are different from those that interpolate stiffness parameters. The approach is similar to the classical one because its main effect is to produce a relaxation of the stress constraints, but it is different in terms of convergence features. The technique is compared with the classical one in the context of stress-constrained minimum-weight topology optimization. Firstly, the problem is studied in a modified truss design framework, where the arising of the singularity phenomenon can be easily shown analytically. Afterwards, the analysis is extended to its natural context of topology bidimensional problems.     

Concurrent topology optimization of multiscale structures with multiple porous materials under random field loading uncertainty

Aim of this work is the synthesis of auxetic structures using a topology optimization approach for micropolar (or Cosserat) materials. A distributed compliant mechanism design problem is formulated, adopting a SIMP–like model to approximate the constitutive parameters of 2D micropolar bodies. The robustness of the proposed approach is assessed through numerical examples concerning the optimal design of structures that can expand perpendicularly to an applied tensile stress. The influence of the material characteristic length on the optimal layouts is investigated. Depending on the inherent flexural stiffness of micropolar solids, truss–like solutions typical of Cauchy solids are replaced by curved beam–like material distributions. No homogenization technique is implemented, since the proposed design approach applies to elements made of microstructured material with prescribed properties and not to the material itself.     

On the influence of local and global stress constraint and filtering radius on the design of hinge-free compliant mechanisms

Distributed compliant mechanisms are components that use elastic strain to obtain a desired kinematic behavior. Compliant mechanisms obtained via topology optimization using the standard approach of minimizing/maximizing the output displacement with a spring at the output port, representing the stiffness of the external medium, usually contain one-node connected hinges. Those hinges are undesired since an ideal compliant mechanism should be a continuous part. This work compares the use of two strategies for stress constrained problems: local and global stress constraints, and analyses their influence in eliminating the one-node connected hinges. Also, the influence of spatial filtering in eliminating the hinges is studied. An Augmented Lagrangian formulation is used to couple the objective function and constraints, and the resulting optimization problem is solved by using an algorithm based on the classical optimality criteria approach. Two compliant mechanisms problems are studied by varying the stress limit and filtering radius. It is observed that a proper combination of filtering radius and stress limit can eliminate one-node connected hinges.     

Topology optimization of cellular materials with periodic microstructure under stress constraints

Material design is a critical development area for industries dealing with lightweight construction. Trying to respond to these industrial needs topology optimization has been extended from structural optimization to the design of material microstructures to improve overall structural performance. Traditional formulations based on compliance and volume control result in stiffness-oriented optimal designs. However, strength-oriented designs are crucial in engineering practice. Topology optimization with stress control has been applied mainly to (macro) structures, but here it is applied to material microstructure design. Here, in the context of density-based topology optimization, well-established techniques and analyses are used to address known difficulties of stress control in optimization problems. A convergence analysis is performed and a density filtering technique is used to minimize the risk of results inaccuracy due to coarser finite element meshes associated with highly non-linear stress behavior. A stress-constraint relaxation technique (qp-approach) is applied to overcome the singularity phenomenon. Parallel computing is used to minimize the impact of the local nature of the stress constraints and the finite difference design sensitivities on the overall computational cost of the problem. Finally, several examples test the developed model showing its inherent difficulties.

     

Stress-based design of thermal structures via topology optimization

The design of thermal structures in the aerospace industry, including exhaust structures on embedded engine aircraft and hypersonic thermal protection systems, poses a number of complex design challenges. These challenges are particularly well addressed by the material layout capabilities of structural topology optimization; however, no topology optimization methods are readily available with the necessary thermoelastic considerations for these problems. This is due in large part to the emphasis on cases of maximum stiffness design for structures subjected to externally applied mechanical loads in the majority of topology optimization applications. In addition, while limited work in the literature has investigated thermoelastic topology optimization, a direct treatment of thermal stresses is not well documented. Such a treatment is critical in the design of thermal structures where excessive thermal stresses are a primary failure mode. In this paper, we present a method for the topology optimization of structures with combined mechanical and thermoelastic (temperature) loads that are subject to stress constraints. We present the necessary steps needed to address both the design-dependent thermal loads and accommodate the challenges of stress-based design criteria. A relaxation technique is utilized to remove the singularity phenomenon in stresses and the large number of stress constraints is handled using a scaled aggregation technique that has been shown previously to satisfy prescribed stress limits in mechanical problems. Finally, the stress-based thermoelastic formulation is applied to two numerical example problems to demonstrate its effectiveness.     

Topology optimization of frame structures—joint penalty and material selection

This paper deals with joint penalization and material selection in frame topology optimization. The models used in this study are frame structures with flexible joints. The problem considered is to find the frame design which fulfills a stiffness requirement at the lowest structural weight. To support topological change of joints, each joint is modelled as a set of subelements. A set of design variables are applied to each beam and joint subelement. Two kinds of design variables are used. One of these variables is an area-type design variable used to control the global element size and support a topology change. The other variables are length ratio variables controlling the cross section of beams and internal stiffness properties of the joints. This paper presents two extensions to classical frame topology optimization. Firstly, penalization of structural joints is presented. This introduces the possibility of finding a topology with less complexity in terms of the number of beam connections. Secondly, a material interpolation scheme is introduced to support mixed material design.     

Optimum material design of minimum structural compliance under seepage constraint

In filtration and chemical engineering industry the load carrying capacity and seepage performances are very important for a successful filter design. We study a two-scale structural design optimization problem to minimize structural compliance under given seepage flow rate and material porosity constraints. Structural size, shape and topology are given because of other functional requirements. Structural material used is macro homogeneous porous material with periodic microstructure and is to be designed. Since structural compliance and seepage performances in macro-scale are implicit functions of material microstructural topology, it becomes a two-scale design optimization problem. The cross scale sensitivities are derived by the adjoint method. A new volume preserving nonlinear density filter is proposed which makes the process of optimization iteration more stable. The optimization problem is solved by GCMMA. Examples under the equality constraints of different seepage flow rate are presented to illustrate the effectiveness of two-scale design optimization formulation and solution approach.     

Optimal design of periodic functionally graded composites with prescribed properties

The computational design of a composite where the properties of its constituents change gradually within a unit cell can be successfully achieved by means of a material design method that combines topology optimization with homogenization. This is an iterative numerical method, which leads to changes in the composite material unit cell until desired properties (or performance) are obtained. Such method has been applied to several types of materials in the last few years. In this work, the objective is to extend the material design method to obtain functionally graded material architectures, i.e. materials that are graded at the local level (e.g. microstructural level). Consistent with this goal, a continuum distribution of the design variable inside the finite element domain is considered to represent a fully continuous material variation during the design process. Thus the topology optimization naturally leads to a smoothly graded material system. To illustrate the theoretical and numerical approaches, numerical examples are provided. The homogenization method is verified by considering one-dimensional material gradation profiles for which analytical solutions for the effective elastic properties are available. The verification of the homogenization method is extended to two dimensions considering a trigonometric material gradation, and a material variation with discontinuous derivatives. These are also used as benchmark examples to verify the optimization method for functionally graded material cell design. Finally the influence of material gradation on extreme materials is investigated, which includes materials with near-zero shear modulus, and materials with negative Poisson’s ratio.     

Structural topology optimization for frequency response problem using model reduction schemes

This study uses model reduction (MR) schemes such as the mode superposition (MS), Ritz vector (RV), and quasi-static Ritz vector (QSRV) methods, which reduce the size of the dynamic stiffness matrix of dynamic structures, to calculate dynamic responses and sensitivity values with adequate efficiency and accuracy for topology optimization in the frequency domain. The calculation of structural responses to dynamic excitation using the framework of the finite element (FE) procedure usually requires a significant amount of computation time; that is mainly attributable to repeated inversions of dynamic stiffness matrices depending on time or frequency intervals, which hastens the dissemination of the MR schemes in the analysis. However, using well-established MR schemes in topology optimization has not been prevalent. Therefore, this study conducted a comprehensive investigation to highlight the drawbacks and advantages of these MR schemes for topology optimization. In the results, the MS method, which generates reduction bases by considering some of the lowest eigenmodes, can lose the accuracy in both approximated structural responses and sensitivity values due to locally vibrating eigenmodes and higher mode truncation in the solid isotropic material with penalization (SIMP) approach. In addition, the RV and QSRV methods, which generate reduction bases by considering the external force, mass, and stiffness matrices of a structure, can be used as alterative model reduction schemes for stable optimization. Through several analysis and design examples, the efficiency and reliability of the model reduction schemes for topology optimization are compared and validated.     

Solving stress constrained problems in topology and material optimization

This article is a continuation of the paper Ko?vara and Stingl (Struct Multidisc Optim 33(4?C5):323?C335, 2007). The aim is to describe numerical techniques for the solution of topology and material optimization problems with local stress constraints. In particular, we consider the topology optimization (variable thickness sheet or ??free sizing??) and the free material optimization problems. We will present an efficient algorithm for solving large scale instances of these problems. Examples will demonstrate the efficiency of the algorithm and the importance of the local stress constraints. In particular, we will argue that in certain topology optimization problems, the addition of stress constraints must necessarily lead not only to the change of optimal topology but also optimal geometry. Contrary to that, in material optimization problems the stress singularity is treated by the change in the optimal material properties.     

Reliability-based topology optimization of continuum structures subject to local stress constraints

Topology optimization of continuum structures is a challenging problem to solve, when stress constraints are considered for every finite element in the mesh. Difficulties are compounding in the reliability-based formulation, since a probabilistic problem needs to be solved for each stress constraint. This paper proposes a methodology to solve reliability-based topology optimization problems of continuum domains with stress constraints and uncertainties in magnitude of applied loads considering the whole set of local stress constrains, without using aggregation techniques. Probabilistic constraints are handled via a first-order approach, where the principle of superposition is used to alleviate the computational burden associated with inner optimization problems. Augmented Lagrangian method is used to solve the outer problem, where all stress constraints are included in the augmented Lagrangian function; hence sensitivity analysis may be performed only for the augmented Lagrangian function, instead of for each stress constraint. Two example problems are addressed, for which crisp black and white topologies are obtained. The proposed methodology is shown to be accurate by checking reliability indices of final topologies with Monte Carlo Simulation.     

Topology optimization of couple-stress material structures

Conventional topology optimization is concerned with the structures modeled by classical theory of mechanics. Since it does not consider the effects of the microstructures of materials, the classical theory can not reveal the size effect due to material’s heterogeneity. Couple-stress theory, which takes account of the microscopic properties of the material, is capable of describing the size effect in deformations. The purpose of this paper is to investigate the formulation for topology optimization of couple-stress material structures. The artificial material density of each element is chosen as design variable. Based on the basic idea of SIMP (Solid Isotropic Material with Penalization) method, the effective material stiffness matrix of couple-stress material is related to the artificial density by power law with penalty. The structural analysis is implemented by finite element method for couple-stress materials, and a 4-noded quadrilateral couple-stress element is formulated in which C ₁ continuity requirement is relaxed. Some typical problems are solved and the optimal results based on the couple-stress theory are compared with the conventional ones. It is found that the optimal topologies of couple-stress continuum show remarkable size effect.     

Local and global Pareto dominance applied to optimal design and material selection of composite structures

The optimal design of hybrid composite structures considering sizing, topology and material selection is addressed in a multi-objective optimization framework. The proposed algorithm, denoted by Multi-objective Hierarchical Genetic Algorithm (MOHGA), searches for the Pareto-optimal front enforcing population diversity by using a hierarchical genetic structure based on co-evolution of multi-populations. An age structured population is used to store the ranked solutions aiming to obtain the Pareto front. A self-adaptive genetic search incorporating Pareto dominance and elitism is presented. Two concepts of dominance are used: the first one denoted by local non-dominance is implemented at the isolation stage of populations and the second one called global non-dominance is considered at age structured population. The age control emulates the human life cycle and enables to apply the species conservation paradigm. A new mating and offspring selection mechanisms considering age control and dominance are adopted in crossover operator applied to age-structured population. Application to hybrid composite structures requiring the compromise between minimum strain energy and minimum weight is presented. The structural integrity is checked for stress, buckling and displacement constraints considered in the multi-objective optimization. The design variables are ply angles and ply thicknesses of shell laminates, the cross section dimensions of beam stiffeners and the variables associated with the material distribution at laminate level and structure level. The properties of the proposed approach are discussed in detail.     

Non-probabilistic reliability-based topology optimization with multidimensional parallelepiped convex model

In this paper, a new non-probabilistic reliability-based topology optimization (NRBTO) method is proposed to account for interval uncertainties considering parametric correlations. Firstly, a reliability index is defined based on a newly developed multidimensional parallelepiped (MP) convex model, and the reliability-based topology optimization problem is formulated to optimize the topology of the structure, to minimize material volume under displacement constraints. Secondly, an efficient decoupling scheme is applied to transform the double-loop NRBTO into a sequential optimization process, using the sequential optimization & reliability assessment (SORA) method associated with the performance measurement approach (PMA). Thirdly, the adjoint variable method is used to obtain the sensitivity information for both uncertain and design variables, and a gradient-based algorithm is employed to solve the optimization problem. Finally, typical numerical examples are used to demonstrate the effectiveness of the proposed topology optimization method.