HS-4450: Multi-Objective Optimization of a Cantilever Ibeam using an Inclusion Matrix

Learn how to use the data created from a DOE and pass it to an optimization problem via an inclusion matrix which re-uses the data.

Before you begin, copy the model files used in this tutorial from <hst.zip>/HS-4450/ to your working directory.

The inclusion matrix feature passes an already existing set of data to the running process. In this tutorial, the data created from a DOE is passed to an optimization problem which reuses the data. This promotes efficient design exploration practices: an optimization using a direct solver call can still be done in combination with a DOE to study the system without any loss of data. This

Perform the Study Setup

  1. Start HyperStudy.
  2. Start a new study in the following ways:
    • From the menu bar, click File > New.
    • On the ribbon, click .
  3. In the Add Study dialog, enter a study name, select a location for the study, and click OK.
  4. Go to the Define Models step.
  5. Add a Parameterized File model.
    1. From the Directory, drag-and-drop the ibeam.tpl file into the work area.
      Figure 1.


    2. In the Solver input file column, enter ibeam.py.
      This is the name of the solver input file HyperStudy writes during the evaluation.
    3. In the Solver execution script column, select Python (py).
  6. Click Import Variables.
    Four input variables are imported from the ibeam.tpl resource file.
  7. Go to the Define Input Variables step.
  8. Review the input variable's lower and upper bounds ranges.

Perform Nominal Run

  1. Go to the Test Models step.
  2. Click Run Definition.
    An approaches/setup_1-def/ directory is created inside the study Directory. The approaches/setup_1-def/run__00001/m_1 directory contains the input file, which is the result of the nominal run.

Create and Evaluate Output Responses

  1. Go to the Define Output Responses step.
  2. Create the Iy output response for the y-axis moment of inertia.
    1. From the Directory, drag-and-drop the output.hstp file, located in approaches/setup_1-def/run__00001/m_1, into the work area.
    2. In the File Assistant dialog, set the Reading technology to Altair® HyperWorks® and click Next.
    3. Select Single item in a time series, then click Next.
    4. Define the following options, then click Next.
      • Set Type to Output.
      • Set Request to Iy.
      • Set Component to Value.
      Figure 2.


    5. Label the output response Iy.
    6. Set Expression to Maximum.
    7. Click Finish.
      Figure 3.


  3. Create four more output responses by repeating step 2, but change the Request assigned to each output according to Table 1:
    Table 1.
    Output Response Request
    Volume Vol
    IZ Iz
    Displacement d
    Frequency1 Freq
  4. Click Evaluate to extract the response values.

Run a Hammersley DOE Study

  1. Add a DOE.
    1. In the Explorer, right-click and select Add from the context menu.
      The Add dialog opens.
    2. From Select Type, choose DOE.
    3. For Definition from, select an approach.
    4. Select Setup and click OK.
  2. Go to the DOE 1 > Specifications step.
  3. In the work area, set the Mode to Hammersley.
  4. In the Settings tab, verify the Number of Runs is 17.
    Figure 4.


  5. Click Apply.
  6. Go to the DOE 1 > Evaluate step.
  7. Click Evaluate Tasks.
  8. Go to the DOE 1 > Post-Processing step.
  9. Click the Pareto Plot tab.
    In the options menu, ensure that Linear effects is enabled.
    A Pareto Plot shows the ranked influence of the input variables on the output response. For example, for the y-axis moment of the inertia, height has the largest influence and web thickness has the least. In contrast, for the z-axis moment of inertia, the flange length and flange thickness are the most influential variables. The size of the bar indicates the magnitude of the influence, and the hashed line’s slope indicates the sign of the effect: positive or negative. For example, increasing the height will increase Iy, but it will decrease displacement.
    Figure 5.


Run Optimization

  1. Add an Optimization.
    1. In the Explorer, right-click and select Add from the context menu.
    2. In the Add dialog, select Optimization.
    3. For Definition from, select Setup and click OK.
  2. Go to the Optimization > Definition > Define Output Responses step.
  3. Click the Objectives/Constraints - Goals tab.
  4. Apply an objective to the Volume output response.
    1. Click Add Goal.
    2. In the Apply On column, select Volume.
    3. In the Type column, select Minimize.
    Figure 6.


  5. Apply an objective to the Iy output response.
    1. Click Add Goal.
    2. In the Apply On column, select Iy.
    3. In the Type column, select Maximize.
    Figure 7.


  6. Go to the Optimization > Specifications step.
  7. In the work area, set the Mode to Global Response Search Method (GRSM).
  8. In the Settings tab, change the Number of Runs to 50.
  9. Click the More tab and define the following settings:
    1. Set Points per Iteration to 4.
    2. Set Use Inclusion Matrix to Without Initial.
    GRSM performs a global search, therefore the initial values of the variables are not important and do not have to be used within the optimization.
    Figure 8.


  10. Import run data from the DOE using an Inclusion Matrix.
    1. Click Edit Matrix > Inclusion Matrix from the top, right corner of the work area.
      Figure 9.


    2. In the Edit Inclusion Matrix dialog, click Import Values.
    3. In the Import Values dialog, select Approach evaluation data and click Next.
    4. Set the approach to DOE 1 and click Next.
    5. Click Finish.
    6. Review the imported run data and click OK.
  11. Click Apply.
  12. Go to the Evaluate step.
  13. Click Evaluate Tasks.
  14. Go to the Post-Processing step.
  15. Click the Optima tab.
  16. Observe the non-dominated front of designs.
    These points represent the trade-off between the objective of minimizing volume and maximizing the y-axis moment of inertia. In the plot, it is evident that as the moment of inertia increases, the volume increases as well. This curve represents the trade-off of the best available designs given the competing objective requirements.
    Figure 10.


  17. Click the Scatter tab to plot the objectives along the same axes shown in the Optima plot.
    This scatter plot shows all of the runs from the optimization. When comparing the scatter and optima plots, note that the optima plot contains only a subset of runs which are non-dominated. A dominated design is a design for which both objectives could be improved. A non-dominated design is one in which one objective may only be improved at the expense of another.
    Figure 11.