12.3.1 Adaptive remeshing: overview

ABAQUS/CAE provides an automated process to remesh your model adaptively. The goal of the adaptive remeshing process is to approach or reach targets on selected error indicators for a specified model and its accompanying load history. See Adaptivity techniques, Section 12.1.1, for a comparison of this process to other ABAQUS adaptivity methods.

Overview

The following steps are required to incorporate adaptive remeshing into your ABAQUS/CAE model:

Based on these remeshing rules and your adaptivity process definition, ABAQUS/CAE iteratively:
  • executes an ABAQUS/Standard analysis, which will write selected error indicator output variables based on your remeshing rule settings (See Error indicators, Section 12.3.2),

  • uses the error indicator variables in a sizing function to compute element sizes for a new mesh, respecting any size constraints you might specify (See Solution-based mesh sizing, Section 12.3.3), and

  • generates a new mesh in the regions specified, based on the computed element sizes. The neighboring regions will also be remeshed.

These iterations continue until either:
  • all remeshing rule targets are satisfied, or

  • a maximum number of remesh iterations is reached.

See When will my mesh adaptivity stop iterating?, Section 18.3.2 of the ABAQUS/CAE User's Manual, for more details. Figure 12.3.1–1 shows the interaction of ABAQUS products and files in this process.

Figure 12.3.1–1 User actions and automated ABAQUS/CAE actions in the adaptive remeshing process.

Typical applications

Adaptive remeshing can improve the quality of your simulation results. Adaptive remeshing can be helpful when:

  • you are unsure how refined a mesh needs to be to reach a particular level of accuracy or how coarse the mesh can be without unacceptably impacting solution accuracy;

  • it is difficult to design an adequately refined mesh near a region of interest, such as near a stress riser; or

  • you do not know a location of interest, such as with formation of a plastic zone, a priori.

An example of using adaptive remeshing to study the thermal and stress behavior of a bolted vessel is provided in Thermal-stress analysis of a reactor pressure vessel bolted closure, Section 5.1.6 of the ABAQUS Example Problems Manual. The example includes a Python script that you can run from ABAQUS/CAE to create the model and the remeshing rules. A second script allows you to submit the adaptivity process and to view the changing mesh as ABAQUS/CAE computes new element sizes.

Example: stress riser

Figure 12.3.1–2 shows how adaptive remeshing generates a high-quality mesh for a typical notched specimen subjected to axial loading.

Figure 12.3.1–2 Stress riser mesh before and after refinement.

Figure 12.3.1–3 shows the effect of these mesh changes on solution accuracy in comparison to the effect of uniform mesh refinement on solution accuracy. Adaptive mesh refinement is much more efficient than uniform mesh refinement at reducing solution error.

Figure 12.3.1–3 Comparison of adaptive remeshing to uniform mesh refinement based on boundary seeding.

Example: plastic hinge

This example, a doubly-notched specimen axially strained until a plastic hinge or band forms, is used to demonstrate how adaptive remeshing will focus a mesh on a plastic hinge. It illustrates the value of adaptive remeshing in cases where the region of interest may not be known a priori. Figure 12.3.1–4 shows the specimen and the region of active yielding. Figure 12.3.1–5 shows the original mesh and the adapted mesh after three adaptive remeshing iterations.

Figure 12.3.1–4 Region of active yielding in a doubly-notched specimen.

Figure 12.3.1–5 Mesh of doubly-notched specimen before and after adaptive remeshing.

Preparing your model for adaptive remeshing

You use ABAQUS/CAE to do the following when performing adaptive remeshing:

  • create the model and specify the boundary conditions and loading history,

  • create remeshing rules,

  • create an adaptivity process, and

  • start and monitor the progress of the adaptivity process.

Creating the model

You do not have to consider adaptive remeshing when you create the model and specify the boundary conditions and loading history; however, before using adaptive remeshing you must do the following:

  • create the geometry of the model—you cannot use an orphan mesh part—and

  • provide an initial, nominal, mesh. This mesh can be fairly coarse. Providing an extremely coarse mesh, however, can result in more adaptive remesh iterations due to the poor quality of early remesh iteration error indicator calculations. You can, in typical cases, define a reasonable initial mesh by using the default part instance mesh seeding in ABAQUS/CAE.

Creating a remeshing rule

You create and configure a remeshing rule using the Mesh module in ABAQUS/CAE. See Creating a remeshing rule, Section 17.19.1 of the ABAQUS/CAE User's Manual, for details on defining remesh rules. Refer to Error indicators, Section 12.3.2, and Solution-based mesh sizing, Section 12.3.3, for details on the methods used to determine revised mesh size distributions.

ABAQUS/CAE Usage: 

Mesh module: AdaptivityRemeshing RuleCreate


Creating an adaptivity process

You create and configure an adaptivity process using the Job module in ABAQUS/CAE. When you create an adaptivity process, you can specify the maximum number of remesh iterations to be performed and set various system resource parameters. See Creating, editing, and manipulating jobs, Section 18.5 of the ABAQUS/CAE User's Manual, for details.

ABAQUS/CAE Usage: 

Job module: AdaptivityCreate


Performing adaptive remeshing with a provisional analysis

In some cases you will want to determine an adequate mesh for your model prior to conducting a fully detailed analysis, which might include many steps and complex behavior. A “provisional” analysis can often be used, along with adaptive remeshing, to efficiently determine a good mesh for a model. The provisional analysis may include various simplifications of your fully detailed analyis, such as

  • replacing your steps with a single linear perturbation step with loading that adequately reflects your more general loading cases,

  • removing plasticity and other material nonlinearities, and

  • disabling geometric nonlinearity.

The provisional analysis approach may result in a mesh that is not ideally suited to your ultimate choice of loading. However, the cost for obtaining a mesh from a provisional model may be significantly lower than the case where your adaptivity process considers all of the complexity in the fully detailed analysis, and you may find the refined mesh adequate for use in a variety of analysis situations.

Special considerations

In general, the ABAQUS adaptive remeshing process iterates automatically toward a better quality mesh; however, you should be aware of certain considerations.

Singularities

Stress singularities frequently result from geometric abstractions, such as reentrant corners and contact of a sharp edge in elastic materials, and from point loads or abruptly ended distributed load regions. In these situations the stress field near the singularity is unbounded, and no amount of mesh refinement will enable resolution of the correct solution. If you apply the adaptive remeshing process to regions of your model that include singularities, the process will drive elements near the singularity to very small sizes in an unbounded fashion. The end result may be unacceptably expensive analyses.

You can prevent excessively expensive analyses of models with singularities using the following techniques:

  • Exclude the region of the singularity from consideration in the remeshing process. You exclude a region by partitioning the model and assigning remeshing rules only to regions away from the singularity.

  • Apply a minimum element size constraint in the remeshing rule. ABAQUS/CAE does assign a minimum element size by default, which is a fraction of the default part instance mesh seed. You can modify this constraint to achieve a quality solution near the singularity while avoiding an excessively refined mesh. You can also use the remeshing rule to control the rate at which ABAQUS/CAE refines the size of the elements. Element size constraints may prevent an adaptivity process from achieving specified error indicator targets.

Convergence issues

Figure 12.3.1–6 shows a typical history of an error indicator and the computational cost, in ABAQUS/Standard, versus remesh iteration.

Figure 12.3.1–6 Error indicator and computational cost versus iteration for a model with a 25% error indicator target.

The example in Figure 12.3.1–6 shows a desirable convergence profile. The solution error indicator decreases monotonically and quickly to the desired 25% error indicator target. Accompanying this error indicator decrease is a moderate increase in computational cost, measured either in model degrees of freedom or time in ABAQUS/Standard. Certain situations can interfere with this desirable convergence profile, as follows:

  • If your initial mesh is too coarse, the error indicator variables may be of insufficient quality to result in a mesh that is sufficiently improved in the next iteration. The adaptive remeshing process typically creates a high-quality mesh eventually even if the initial mesh is quite coarse. However, some mesh iterations can be avoided with a reasonably refined initial mesh.

  • Minimum element size constraints that you specify when creating the remeshing rule can prevent the mesh from achieving sufficient refinement (in the extreme case of singularities this will always be the case) to satisfy your error indicator targets. You may be able to satisfy your targets by relaxing these constraints; for example, by decreasing the minimum element size. For more information, see What are remeshing rules?, Section 17.11.1 of the ABAQUS/CAE User's Manual.

  • Singularities can cause unbounded focusing of the mesh and an inability of the adaptive remeshing process to achieve error targets. As described in “Singularities,” above, you can control the computational cost by specifying a minimum element size constraint. In any case where a singularity does exist within a remeshing rule region, you may see poor convergence in error indicator results.

Continuing a stopped adaptive remeshing process

The adaptive remeshing process is designed to be automatic—ABAQUS/CAE performs a sequence of analyses as it continues to refine your mesh. However, there are occasions where the process will stop and you will want to continue adaptive remeshing from your most recent mesh:

  • when you want to change remeshing rules for later remesh iterations, or

  • when the adaptive remesh process fails to complete due to machine resource problems.

You can continue the adaptive remeshing process by resubmitting an existing adaptivity process, creating and submitting a new adaptivity process, or performing manual remeshing. See Manually resizing and remeshing, Section 17.19.6 of the ABAQUS/CAE User's Manual.

Limitations

Adaptive remeshing requires the use of ABAQUS/CAE, and only ABAQUS/Standard procedures are supported. Other specific limitations also apply.

Element types

ABAQUS/CAE can perform adaptive remeshing only with elements of the following shapes (see Which mesh controls can I use with adaptive remeshing?, Section 17.11.2 of the ABAQUS/CAE User's Manual):

  • Planar continuum triangles and quadrilaterals

  • Shell triangles and quadrilaterals

  • Tetrahedrals