** **

By Franck Delcroix

Program Manager – HyperStudy

**Multi-Objective Optimization in the Design Cycle**

Design optimization seeks a design that minimizes the objective function while satisfying design constraints. Real-world design problems are usually characterized by the presence of many conflicting objectives. Therefore, it is natural to look at the engineering design problems as multi-objective optimization problems. For example, we may want to maximize range and payload mass while trying to minimize manufacturing costs for an airplane design.

Figure 1. Multiple Objectives : A design dilemma

So, in practice, the designer often needs to consider several design criteria or objective functions simultaneously. In most cases, trade-offs exist among the objectives, where improvement in one objective cannot be achieved without deteriorating another. As a consequence, it is rare for a multi-objective optimization problem to produce a single optimal solution; instead, a set of equally valid solutions will be proposed.

In this article, we will demonstrate the design process of a control arm with multiple objectives using Altair HyperWorks.

** **

**Introduction to Multi-Objective Optimization**

Design variables are a set of independent parameters that can be changed to improve system performance. The objective function measures the quality of a design. Constraints are design requirements, and they must be satisfied by the optimal design. Both the objective and the constraints are functions of design variables. So, the optimization problem reads: Find the design variable values that minimize (maximize) the objective function(s) and satisfy the design constraints.

A multi-objective optimization (MOO) problem is formulated as follows:

minimize f(x)={f_{1}(x),f_{2}(x),…,f_{n}(x)}

subject to g_{j}(x) < 0

where x is the design variables, f_{1}(x), …, f_{n}(x) are multiple objective functions and g_{j}(x) are the constraints.

When dealing with multiple objectives (f_{1}, f_{2}, …) the natural “better than”/“worse than” relationships between alternative solutions no longer hold, and the relationship among individuals is described as a “domination” relationship, leading to Pareto-Optimal set instead of a single optimal design. Let’s introduce these two key definitions.

Domination

** **

- For any two solutions x
_{1}and x_{2}, x_{1}is said to dominate x_{2}if these conditions hold:- x
_{1}is not worse than x_{2}in all objectives - x
_{1}is strictly better than x_{2}in at least one objective

- x
- If one of the above conditions does not hold x
_{1}does not dominate x_{2}

** **

Pareto-Optimal set

- Non-dominated set

The set of all solutions that are not dominated by any other solution in the design space.

- Global Pareto optimal set

There exists no other solution in the entire design space that dominates any member of the set.

Navigation along the Pareto front allows engineers to do trade-off analyses and pick the best design for various scenarios.

From these definitions, in MOO, the approach is to search the design space for a set of Pareto optimal solutions, from which the designer can choose the final design. Pareto optimality is defined as a set where every element is a possible solution for which no other solutions can be better in *all* design objectives. A solution in the Pareto optimal set cannot be declared as “better” than others in the set without using other information (such as preference) to rank the competing objectives. For the two-dimensional case, the Pareto front is a curve that clearly illustrates the tradeoff between the two objectives (see figure 2).

* *

Figure 2. Pareto Front and non-domination illustration

** **

**Handling Multi-Objective Optimization Problems with HyperStudy**

Altair’s HyperStudy provides a number of algorithms (SQP, ARSM, MFD, GA, SORA) to cover a wide range of optimization problems: constrained optimization, gradient-based methods vs. heuristic approaches, continuous and discrete variables, adaptive response surface based method, reliability-based design optimization.

Multi-objective optimization problems cannot be solved using the usual design optimization methods for a scalar objective function, and specific algorithms are needed instead. HyperStudy’s offering for MOO includes:

- Weighted sum approach
- MOGA (Multi-Objective Genetic Algorithm)
- GMMO (Gradient based Method for Multi-Objective Optimization)

** Weighted Sum** is an approach for treating a MOO problem as a single objective optimization problem by summing up the weighted individual objectives. HyperStudy provides a user-friendly GUI to perform this task:

* min f(x)={f _{1}(x),f_{2}(x),…,f_{n}(x)} *

*à*

*min Σ ω*

_{i}.f_{i}Since the problem is no longer a MOO problem, any standard optimizer of HyperStudy (GA, ARSM, MFD, SQP) can be used to solve it.

**Multi-objective genetic algorithm**, unlike the weighted sum approach, produces a set of Pareto-optimal solutions (non-dominated solutions). The MOGA implementation in HyperStudy is based on the standard feature of GA and includes:

- Non-dominated classification strategy
- Crowding distance evaluation to help create a good distribution of the points on the Pareto front.
- Storage of all the global non-dominated points
- Flexible termination criterion for better run control

**Gradient-based method for multi-objective optimization** also produces a set of Pareto-optimal solutions (non-dominated solutions). The GMMO implementation in HyperStudy includes:

- A proprietary method that extends a typical gradient-based algorithm to multi-objective formulation
- A method constructed for the orderly and efficient exploration of the Pareto-front
- A method where the user can control the number of analyses to be performed

**Multi-Objective Optimization of an Arm**

This example illustrates the tradeoff challenges faced by the designer of a typical mechanical component. The arm shown below is clamped at one end and is under an axial loading on the other end (see figure 4).

The model has been meshed and modeled in HyperMesh. Linear static analysis is performed using RADIOSS, the HyperWorks finite element solver,.

Figure 3. Arm model

Figure 4. Boundary and loading conditions

The problem is stated as follows: minimize the mass and minimize the maximum displacement, but respect an upper limit on the stress level throughout all the part.

The current (nominal) design gives a volume of 1.7667E06 mm^{3}, a maximum displacement of 1.41 mm and a maximum stress of 195.29 MPa.

The designer can change the following properties of the arm: six shape variables for the overall length and height of the part, three shape variables for the three radii (see figure 5 below).

Figure 5. Shapes defined on the model

The first step is to create the shape variables using HyperMesh. HyperMesh includes HyperMorph, a tool for doing parametric mesh-based shape changes (morphing). In the morphing process, domains are created around the FE model features. Handles are defined on the edges of these domains. By moving a handle, nodes in the associated domain also will move. Each of these movements then can be saved as a shape variable and can be exported directly to HyperStudy, the HyperWorks solver-neutral design exploration, study and optimization tool.

In HyperStudy, the study is set up by identifying these variables (nine shape variables), setting up the simulations to run (RADIOSS), and extracting the responses (maximum displacement, volume, maximum stress).

Figure 6. Table of all the design variables

After defining the problem, a nominal run is executed within HyperStudy. Once the user makes sure that nominal run is executed correctly, the next task is to pick a study method (i.e. DOE, optimization, stochastics). In this case, an optimization based on an approximation is chosen.

Figure 7. Study flow and solver execution in HyperStudy

The study starts with an attempt to reduce the number of variables. A ffactional Factorial (16 runs) DOE (design of experiment), using the nine variables is realized. The main effects plots (figure 8 ) show that three design variables (the three radii are of smaller impact on the responses than the other six variables and thus can be screened out.

Only six design variables are then kept to perform:

- A full factorial DOE (64 runs) and an Hammersley DOE (86 runs), to be used as input matrix for the approximation
- A complementary Latin HyperCube DOE (20 runs) to be used for approximation validation

Moving Least Square Method (MLSM) approximations are built and checked for all responses using 64+86 runs as the input matrix and 20 runs as the validation matrix.

A multi-objective optimization is performed using MOGA. The target is to minimize both the volume and the maximum displacement, while satisfying a constraint on the maximum stress in the part.

Figure 9. MOO setup in HyperStudy

After running 25 iterations of MOGA, with a total of 1,878 analyses on the meta-model, the results of the last iteration as shown below are obtained (figure 10). Note that, contrary to standard optimization, an MOO will produce, for a given iteration, a set of designs (a Pareto front). The front displayed below is the 25^{th} iteration and is made of 116 non-dominated points.

Figure 10. Pareto front obtained at iteration 25 of a MOGA optimization

In this plot, X-axis is the volume and Y-axis is the maximum displacement: The tradeoff is based on selecting the best compromise between the two competing objectives.

The black lines on the graph below show where the nominal point (initial design) is situated.

Figure 11. Initial design vs. Pareto Front

Notice also that MOGA produces a front where all points are admissible (maximum stress < 300 MPa).

For example in the graph below, we can see that, of the 116 non-dominated points of iteration, 25 are all giving a max_stress response under 297.97 MPa (x-axis)

Figure 12. Pareto front for the space (volume ; max_stress)

**Conclusions**

Handling multiple objectives is a constant dilemma for designers. Often, the “most important” objective is kept and other objectives are transformed into constraints. But this does not reflect the tradeoff among all possible designs. Implementation of MOO methods in HyperStudy allows for extended exploration of solution and do tradeoffs.

This study illustrates the use of HyperWorks to search for an optimal arm design such that the volume (mass) is minimized and the maximum displacement is minimized, while meeting stress constraints. In this study, HyperStudy automatically updates the design using HyperMorph for shape variables, runs RADIOSS as the FE solver and uses a specific algorithm for multi-objective optimization. With this process, a set of solutions is proposed to the designer to select the best compromise.