Improvements in the treatment of stress constraints in structural topology optimization problems

Topology optimization of continuum structures is a relatively new branch of the structural optimization field. Since the basic principles were first proposed by Bendsøe and Kikuchi in 1988, most of the work has been dedicated to the so-called maximum stiffness (or minimum compliance) formulations. However, since a few years different approaches have been proposed in terms of minimum weight with stress (and/or displacement) constraints.These formulations give rise to more complex mathematical programming problems, since a large number of highly non-linear (local) constraints must be taken into account. In an attempt to reduce the computational requirements, in this paper, we propose different alternatives to consider stress constraints and some ideas about the numerical implementation of these algorithms. Finally, we present some application examples.     

Handling equality constraints in evolutionary optimization

Over the last few decades several methods have been proposed for handling functional constraints while solving optimization problems using evolutionary algorithms (EAs). However, the presence of equality constraints makes the feasible space very small compared to the entire search space. As a consequence, the handling of equality constraints has long been a difficult issue for evolutionary optimization methods. This paper presents a Hybrid Evolutionary Algorithm (HEA) for solving optimization problems with both equality and inequality constraints. In HEA, we propose a new local search technique with special emphasis on equality constraints. The basic concept of the new technique is to reach a point on the equality constraint from the current position of an individual solution, and then explore on the constraint landscape. We believe this new concept will influence the future research direction for constrained optimization using population based algorithms. The proposed algorithm is tested on a set of standard benchmark problems. The results show that the proposed technique works very well on those benchmark problems.     

An active set algorithm for nonlinear optimization with polyhedral constraints

A polyhedral active set algorithm PASA is developed for solving a nonlinear optimization problem whose feasible set is a polyhedron. Phase one of the algorithm is the gradient projection method, while phase two is any algorithm for solving a linearly constrained optimization problem. Rules are provided for branching between the two phases. Global convergence to a stationary point is established, while asymptotically PASA performs only phase two when either a nondegeneracy assumption holds, or the active constraints are linearly independent and a strong second-order sufficient optimality condition holds.     

An algorithm for solving global optimization problems with nonlinear constraints

In this paper we propose an algorithm using only the values of the objective function and constraints for solving one-dimensional global optimization problems where both the objective function and constraints are Lipschitzean and nonlinear. The constrained problem is reduced to an unconstrained one by the index scheme. To solve the reduced problem a new method with local tuning on the behavior of the objective function and constraints over different sectors of the search region is proposed. Sufficient conditions of global convergence are established. We also present results of some numerical experiments.     

On various aspects of evolutionary structural optimization for problems with stiffness constraints

An evolutionary structural optimization (ESO) method for problems with stiffness constraints which is capable of performing simultaneous shape and topology optimization has been recently presented. This paper discusses various aspects of this method such as influences of the element removal ratio, the mesh size and the element type on optimal designs.     

A continuum of solutions to variational inequality with nonlinear constraints: Existence and simplicial algorithm

In this paper, we obtain an existence theorem of a connected set of solutions to a nonlinear variational inequality with explicit nonlinear constraints. This result follows in a constructive way by designing a simplicial algorithm. The algorithm operates on a triangulation of the unbounded regions and generates a piecewise linear path of parametrized stationary points. Each point on the path is an approximate solution.     

A supermemory gradient projection algorithm for optimization problem with nonlinear constraints

In this paper, by extending concept of the supermemory gradient method for unconstrained optimization problems, we present a supermemory gradient projection algorithm for nonlinear programming with nonlinear constraints. Under some suitable conditions we prove its global convergence.     

Robust topology optimization for multi-material structures under interval uncertainty

In this paper, we propose an efficient method to design robust multi-material structures under interval loading uncertainty. The objective of this study is to minimize the structural compliance of linear elastic structures. First, the loading uncertainty can be decomposed into two unit forces in the horizontal and vertical directions based on the orthogonal decomposition, which separates the uncertainty into the calculation coefficients of structural compliance that are not related to the finite element analysis. In this manner, the time-consuming procedure, namely, the nested double-loop optimization, can be avoided. Second, the uncertainty problem can be transformed into an augmented deterministic problem by means of uniform sampling, which exploits the coefficients related to interval variables. Finally, an efficient sensitivity analysis method is explicitly developed. Thus, the robust topology optimization (RTO) problem considering interval uncertainty can be solved by combining orthogonal decomposition with uniform sampling (ODUS). In order to eliminate the influence of numerical units when comparing the optimal results to deterministic and RTO solutions, the relative uncertainty related to interval objective function is employed to characterize the structural robustness. Several multi-material structure optimization cases are provided to demonstrate the feasibility and efficiency of the proposed method, where the magnitude uncertainty, directional uncertainty, and combined uncertainty are investigated.     

An efficient sensitivity computation strategy for the evolutionary structural optimization (ESO) of continuum structures subjected to self-weight loads

This work presents a modified version of the evolutionary structural optimization procedure for topology optimization of continuum structures subjected to self-weight forces. Here we present an extension of this procedure to deal with maximum stiffness topology optimization of structures when different combinations of body forces and fixed loads are applied. Body forces depend on the density distribution over the design domain. Therefore, the value and direction of the loading are coupled to the shape of the structure and they change as the material layout of the structure is modified in the course of the optimization process. It will be shown that the traditional calculation of the sensitivity number used in the ESO procedure does not lead to the optimum solution. Therefore, it is necessary to correct the computation of the element sensitivity numbers in order to achieve the optimum design. This paper proposes an original correction factor to compute the sensitivities and enhance the convergence of the algorithm. The procedure has been implemented into a general optimization software and tested in several numerical applications and benchmark examples to illustrate and validate the approach, and satisfactorily applied to the solution of 2D, 3D and shell structures, considering self-weight load conditions. Solutions obtained with this method compare favourably with the results derived using the SIMP interpolation scheme.     

Optimality and duality theory for stochastic optimization problems with nonlinear dominance constraints

We consider a new class of optimization problems involving stochastic dominance constraints of second order. We develop a new splitting approach to these models, optimality conditions and duality theory. These results are used to construct special decomposition methods.This research was supported by the NSF awards DMS-0303545 and DMS-0303728.Key words.Stochastic programming – stochastic ordering – semi-infinite optimized – decomposition     

Robust truss topology optimization via semidefinite programming with complementarity constraints: a difference-of-convex programming approach

The robust truss topology optimization against the uncertain static external load can be formulated as mixed-integer semidefinite programming. Although a global optimal solution can be computed with a branch-and-bound method, it is very time-consuming. This paper presents an alternative formulation, semidefinite programming with complementarity constraints, and proposes an efficient heuristic. The proposed method is based upon the concave–convex procedure for difference-of-convex programming. It is shown that the method can often find a practically reasonable truss design within the computational cost of solving some dozen of convex optimization subproblems.     

Theoretical justification of interior point algorithms for solving optimization problems with nonlinear constraints

A family of interior point algorithms is considered. These algorithms can be used for solving mathematical programming problems with nonlinear inequality constraints. Some weighted Euclidean norms are applied to finding a descent direction for improving the solution. These norms vary with iterations. A theoretical justification of the algorithms with some assumptions (including the nonsingularity of the problem) is presented.     

Non-monotone trust-region algorithms for nonlinear optimization subject to convex constraints

This paper presents two new trust-region methods for solving nonlinear optimization problems over convex feasible domains. These methods are distinguished by the fact that they do not enforce strict monotonicity of the objective function values at successive iterates. The algorithms are proved to be convergent to critical points of the problem from any starting point. Extensive numerical experiments show that this approach is competitive with the LANCELOT package.     

Progressive genetic algorithm for solution of optimization problems with nonlinear equality and inequality constraints

A new approach, identified as progressive genetic algorithm (PGA), is proposed for the solutions of optimization problems with nonlinear equality and inequality constraints. Based on genetic algorithms (GAs) and iteration method, PGA divides the optimization process into two steps; iteration and search steps. In the iteration step, the constraints of the original problem are linearized using truncated Taylor series expansion, yielding an approximate problem with linearized constraints. In the search step, GA is applied to the problem with linearized constraints for the local optimal solution. The final solution is obtained from a progressive iterative process. Application of the proposed method to two simple examples is given to demonstrate the algorithm.     

Global optimization of nonlinear least-squares problems by branch-and-bound and optimality constraints

We study a simple, yet unconventional approach to the global optimization of unconstrained nonlinear least-squares problems. Non-convexity of the sum of least-squares objective in parameter estimation problems may often lead to the presence of multiple local minima. Here, we focus on the spatial branch-and-bound algorithm for global optimization and experiment with one of its implementations, BARON (Sahinidis in J. Glob. Optim. 8(2):201–205, 1996), to solve parameter estimation problems. Through the explicit use of first-order optimality conditions, we are able to significantly expedite convergence to global optimality by strengthening the relaxation of the lower-bounding problem that forms a crucial part of the spatial branch-and-bound technique. We analyze the results obtained from 69 test cases taken from the statistics literature and discuss the successes and limitations of the proposed idea. In addition, we discuss software implementation for the automation of our strategy.     

Matrix factorizations in optimization of nonlinear functions subject to linear constraints

Several ways of implementing methods for solving nonlinear optimization problems involving linear inequality and equality constraints using numerically stable matrix factorizations are described. The methods considered all follow an active constraint set approach and include quadratic programming, variable metric, and modified Newton methods.Part of this work was performed while the author was a visitor at Stanford University. This research was supported in part by the National Science Foundation under Grant GJ 36472 and in part by the Atomic Energy Commission Contract No. AT(04-3)-326PA30.     

New differentiated consensus problem with topology cost constraints

A new differentiated consensus problem is studied. The problem is, given a system with multiple classes, consensus is targeted for each class and the consensus values can be different among the classes. Specifically, differentiated consensus is studied in a distributed stochastic network of nodes (or agents), where tasks assigned with different priorities are serviced. The network is assumed to have a switching topology and involves noises, delays in measurements, and topology cost constraints. The goal is to reach a balanced load (i.e., consensus) across the network and, at the same time, to satisfy the topology cost constraint, both for each priority class. A new control protocol is proposed, with which the network resources are allocated in a randomized way with a probability assigned to each priority class. It is shown that the control protocol meets the topology cost constraint and can be used to reach an approximate consensus for each of the priority classes in the network.     

An active set strategy for solving optimization problems with up to 200,000,000 nonlinear constraints

Numerical test results are presented for solving smooth nonlinear programming problems with a large number of constraints, but a moderate number of variables. The active set method proceeds from a given bound for the maximum number of expected active constraints at an optimal solution, which must be less than the total number of constraints. A quadratic programming subproblem is generated with a reduced number of linear constraints from the so-called working set, which is internally changed from one iterate to the next. Only for active constraints, i.e., a certain subset of the working set, new gradient values must be computed. The line search is adapted to avoid too many active constraints which do not fit into the working set. The active set strategy is an extension of an algorithm described earlier by the author together with a rigorous convergence proof. Numerical results for some simple academic test problems show that nonlinear programs with up to 200,000,000 nonlinear constraints are efficiently solved on a standard PC.     

Extreme point programming with nonlinear constraints

This paper summarizes previous results obtained by the authors on methods of solving extreme point mathematical programming problems with linear constraints. It is also shown how these results can be extended to yield an algorithm for solving extreme point mathematical programming problems with nonlinear constraints. Numerical examples to illustrate the algorithms are included.     

A NURBS-based Multi-Material Interpolation (N-MMI) for isogeometric topology optimization of structures

In this paper, the main intention is to propose a new Multi-material Isogeometric Topology Optimization (M-ITO) method for the optimization of multiple materials distribution, where an improved Multi-Material Interpolation model is developed using Non-Uniform Rational B-splines (NURBS), namely the “NURBS-based Multi-Material Interpolation (N-MMI)”. In the N-MMI model, three key components are involved: (1) multiple Fields of Design Variables (DVFs): NURBS basis functions with control design variables are applied to construct DVFs with the sufficient smoothness and continuity; (2) multiple Fields of Topology Variables (TVFs): each TVF is expressed by a combination of all DVFs to present the layout of a distinct material in the design domain; (3) Multi-material interpolation: the material property at each point is equal to the summation of all TVFs interpolated with constitutive elastic properties. DVFs and TVFs are in the decoupled expression and optimized in a serial evolving mechanism. This feature can ensure the constraint functions are separate and linear with respect to TVFs, which can be beneficial to lower the complexity of numerical computations and eliminate numerical troubles in the multi-material optimization. Two kinds of multi-material topology optimization problems are discussed, i.e., one with multiple volume constraints and the other with the total mass constraint. Finally, several numerical examples in 2D and 3D are provided to demonstrate the effectiveness of the M-ITO method.