18.6.2 Configuring a topology optimization task

The Optimization module provides a variety of settings that allow you to configure a topology optimization task. The configuration settings depend on whether you are configuring an optimization task for a general topology optimization or a condition-based topology optimization. The following topics are covered:

Configuring a general topology optimization task

A general topology optimization is a flexible, sensitivity-based optimization that allows you to apply a range of constraints and objective functions to your model. You use the optimization task editor to customize various aspects of a general topology optimization. To locate the editor, select TaskEditoptimization task name from the main menu bar. To specify a general topology optimization, select the Advanced tab and choose General optimization (sensitivity-based).

The following topics are covered:

Configuring basic settings

To configure basic settings:

  1. In the optimization task editor, click the Basic tab.

  2. Choose whether to freeze load or boundary condition regions.

    It is recommended that you freeze regions to which prescribed conditions are applied because you do not want these regions to be removed during the optimization. Freezing these regions stabilizes the optimization and often leads to a significantly lower number of iterations.

Configuring the density settings

To configure density settings:

  1. In the optimization task editor, click the Density tab.

  2. Select the Density update strategy.

    This setting controls the rate at which the Optimization module updates the relative material density of design elements during the optimization. In most cases you should accept the default setting (Normal). However, if the design responses are very sensitive and you have problems fulfilling the constraints, you may need a more conservative rate that requires more optimization iterations.

  3. Do either of the following to specify the relative density of each element during the initial optimization iteration:

    • Select Optimization product default to allow the Optimization module to determine the initial density. If the material volume is selected as a constraint, the Optimization module calculates the initial density such that the volume constraint is fulfilled exactly. If the material volume is selected as an objective function, each element has an initial relative density of 50%.

    • Select Specify and enter a value (0.0 < initial density ≤ 1.0). You should use this option only if volume is selected as an objective function and not as a constraint and if you know, prior to the optimization, that setting the initial density to a larger or smaller value will fulfill other constraints; for example, displacement constraints. You can use a value greater that 0.5 in conjunction with volume constraints to stabilize nonlinear or contact problems and to improve the convergence behavior.

  4. Enter the Minimum density, the Maximum density, and the Maximum change per design cycle.

    The minimum density must be greater than 0.0, and the maximum density must be less than or equal to 1.0. Changing the density bounds is not recommended, in particular the upper bound. You may need to increase the lower bound if the default value leads to a nearly singular stiffness matrix.

    Numerical experiments indicate that a value of 0.25 (default) is acceptable for the maximum change in density. A lower limit in the change of density, such as 0.1, is recommended for complicated design responses and optimization formulations. However, a lower limit often leads to a higher number of optimization iterations.

Configuring the perturbation settings

To configure the perturbation settings:

  1. In the optimization task editor, click the Perturbation tab.

  2. Enter the number of eigenmodes to track. The default value is five, which means that the Optimization module tracks the five lowest eigenfrequencies.

    In some cases many local low frequency eigenmodes appear during the optimization iterations, which leads to a high number of modes to track and degrades performance. You can avoid tracking a high number of modes by choosing the lower bound of the eigenfrequencies to be 25% of the eigenfrequency of interest in the first optimization iteration.

    Mode tracking is not required if your design response will use the Kreisselmaier-Steinhauser formulation to evaluate the eigenfrequencies. Your Abaqus model must include an output request for at least the number of eigenfrequencies you are tracking.

  3. Select the region over which the Optimization module should track the eigenmodes.

    By default, the Optimization module tracks the eigenmodes of all the nodes in the model, which can degrade performance if you have a large model. You can improve the performance by tracking the eigenmodes over only a selected region; for example, over selected surfaces of your model or over points where lumped or rigid masses are attached.

Configuring convergence options

To configure convergence options:

  1. In the optimization task editor, click the Convergence tab.

  2. Specify the Convergence Criteria. The following options allow you to specify the convergence criteria for a general topology optimization:

    Specifying when to start checking for convergence

    You can specify the iteration during which the Optimization module will begin to check the two convergence criteria. The optimization will always continue at least until this value has been reached. The default value is 4.

    Specifying which convergence criterion to check

    You can specify whether the optimization should end when either of the convergence criterion has been fulfilled or both of the criteria have been fulfilled. The default value is that both criteria must be fulfilled.

    Convergence based on the change in optimization function

    You can specify that the optimization will end based on the change in the objective function from one iteration to the next. The default value is 0.001.

    Convergence based on the change in element densities

    Element density is the design variable for a topology optimization. You can specify that the optimization will end based on the average change in the element density from one iteration to the next. The default value is 0.005.

Configuring advanced options

To configure advanced options:

  1. In the optimization task editor, click the Advanced tab.

  2. Select the General optimization algorithm.

  3. Choose whether to Delete soft elements in region.

    During the topology optimization process, the Optimization module distributes a given mass within the design area while it tries to satisfy the constraints and optimize the objective. At the end of the optimization, the structure contains hard (filled) and soft (void) elements. The soft elements have a negligible influence on the stiffness of the structure; but they are still relevant for the number of degrees of freedom of the structure and, hence, influence the speed of the optimization process. The Delete soft elements option allows you to select a region from which soft elements that have only soft neighboring elements will be removed. The deleted elements are reactivated if needed; for example, if the force flow changes during the optimization.

    Choosing to delete soft elements can help Abaqus converge on a solution because those elements would otherwise degenerate or collapse and is recommended when you are optimizing a nonlinear model. In addition. selecting a Conservative density update strategy and a small change in density per design cycle will improve the accuracy of the results. See Configuring the density settings” in “Configuring a topology optimization task, Section 18.6.2, for more information.

  4. If you chose to delete soft elements, you can prevent isolated soft elements from being removed by choosing to delete only soft elements that have neighboring soft elements. You can define a neighboring element as being within the radius specified by the Average edge length (default) or specified by a value that you enter. If the element edge length varies considerably within the mesh, the radius calculated from the average edge length can be misleading.

  5. If you chose to delete soft elements, you can select the method that the Optimization module will use to delete elements:

    Favor continuity (Standard)

    Choose Favor continuity (Standard) and enter a Relative material density threshold to check for continuity before deleting soft elements. If the optimized model contains an “island” of hard elements that are separated from the rest of the model by soft elements, the Optimization module does not remove the soft elements. In addition, the Optimization module retains soft elements that are preventing hard elements from moving with respect to each other; for example, hard elements that share a common edge but not a common face. An element is considered “soft” if its relative material density is less than the threshold value, and the Optimization module removes it from the analysis.

    Favor continuity (Aggressive)

    Choose Favor continuity (Aggressive) and enter a Relative material density threshold to remove soft elements regardless of continuity. An element is considered “soft” if its relative material density is less than the threshold value, and the Optimization module removes it from the analysis.

    Maximum shear strain

    Choose Maximum shear strain and enter a Maximum shear strain threshold. The Optimization module removes elements from the analysis that have a shear strain larger than the threshold.

    Minimum principal strain

    Choose Minimum principal strain and enter a Minimum principal strain threshold. The Optimization module removes elements that have a principal strain lower than the threshold.

    Maximum elastoplastic strain

    Choose Maximum elastoplastic strain and enter a Maximum elastoplastic strain threshold. The Optimization module removes elements that have an elastoplastic strain larger than the threshold.

    Volume compression

    Choose Volume compression and enter a Relative volume compression. The Optimization module removes elements that are compressing and have a relative volume that is lower than the threshold. The relative volume is defined as , where is the deformed element volume and is the original element volume.

    You should choose Volume compression if your model uses shell or membrane elements or if your model is experiencing large deformations.

    Note:  The soft delete method that you select is dependent on the material behavior and the element type, and you may have to experiment to determine the best method and its threshold value. The file, TOSCA.OUT, contains information about the elements that are being removed and will help you determine the best soft delete method and threshold value. The Favor continuity methods provide a default Relative material density threshold of 0.05. In contrast, the strain and volume methods do not provide a default threshold because the appropriate value depends on your model; for example, on the properties of the materials.

  6. Choose the Material interpolation technique and the Penalty factor.

    Optimization generates hard elements with a density close to one or void elements with a density close to zero. Topology optimization introduces elements with a density between one and zero, and the material interpolation technique calculates the relationship between density and stiffness for these intermediate elements. The SIMP (solid isotropic material with penalization) interpolation scheme defines an exponential relationship between the density and the stiffness of an element and is suitable for static problems. The penalty factor should be greater than 1, and numerical experiments indicate that the default value of 3 produces good results. The RAMP (rational approximation of material properties) interpolation scheme is suitable for dynamic problems. The penalty factor should be greater than 0, and numerical experiments indicate that the default value of 3 produces good results.

    By default, the Optimization module selects the SIMP interpolation scheme for static problems and the RAMP interpolation scheme if at least one dynamic load case appears in your model.

Configuring a condition-based topology optimization task

A condition-based topology optimization uses a strain energy objective function and a volume constraint. You use the optimization task editor to customize various aspects of the condition-based topology optimization. To locate the editor, select TaskEditoptimization task name from the main menu bar. To specify a condition-based topology optimization, select the Advanced tab and choose Condition-based optimization.

The following topics are covered:

Configuring basic settings

To configure basic settings:

  1. In the optimization task editor, click the Basic tab.

  2. Choose whether to freeze load or boundary condition regions.

    It is recommended that you freeze regions to which prescribed conditions are applied because you do not want these regions to be removed during the optimization. Freezing these regions stabilizes the optimization and often leads to a significantly lower number of iterations.

Configuring advanced options

To configure advanced options:

  1. In the optimization task editor, click the Advanced tab.

  2. Select the Condition-based optimization algorithm.

  3. Choose whether to Delete soft elements in region.

    During the topology optimization process, the Optimization module distributes a given mass within the design area while it tries to satisfy the constraints and optimize the objective. At the end of the optimization, the structure contains hard (filled) and soft (void) elements. The soft elements have a negligible influence on the stiffness of the structure; but they are still relevant for the number of degrees of freedom of the structure and, hence, influence the speed of the optimization process. The Delete soft elements option allows you to select a region from which soft elements that have only soft neighboring elements will be removed. The deleted elements are reactivated if needed; for example, if the force flow changes during the optimization.

  4. If you chose to delete soft elements, you can prevent isolated soft elements from being removed by choosing to delete only soft elements that have neighboring soft elements. You can define a neighboring element as being within the radius specified by the Average edge length (default) or specified by a value that you enter. If the element edge length varies considerably within the mesh, the radius calculated from the average edge length can be misleading.

  5. If you chose to delete soft elements, you can choose the method that the Optimization module will use to delete elements:

    • Choose Standard deletion to check for continuity before deleting soft elements. If the optimized model contains an “island” of hard elements that are separated from the rest of the model by soft elements, the Optimization module does not remove the soft elements. In addition, the Optimization module retains soft elements that are preventing hard elements from moving with respect to each other; for example, hard elements that share a common edge but not a common face.

    • Choose Aggressive deletion to delete soft elements regardless of continuity.

  6. If desired, enter the value of the Relative material density threshold. An element is considered “soft” if its relative material density that is less than this value, and the Optimization module removes it from the analysis.

  7. Select the rate at which the Optimization module will modify the element properties during a topology optimization. You can select the rate (Very small, Small, Moderate, Medium, or Large) and allow the Optimization module to calculate the number of design cycles required to meet this rate.

    Alternatively, you can select Dynamic and enter the maximum number of design cycles. The minimum number of design cycles is 10, and the default value is 15. A reduction in the number of design cycles can lead to undesired effects in the optimization. Although the resulting structures have the same stiffness (the sum of the strain energy is almost equal for the different results), changing the optimization speed can cause a different truss configuration in the solution.

  8. Select the volume deleted after the first cycle. You can enter a percentage or an absolute value.

    By default, the Optimization module removes 5% of the optimization region volume in the first iteration. In some cases increasing this starting value will accelerate the optimization without influencing the solution, especially for models where relatively low stresses are present in large areas. Conversely, the Optimization module may remove too many elements in the first iteration if the starting value is too high, leading to a failure in the optimization or a coarse structure.