김경종
260 lines
20 KiB
Markdown
260 lines
20 KiB
Markdown
<!-- source-page: 1161 -->
|
||
|
||
# 12.2.7 ALE ADAPTIVE MESHING AND REMAPPING IN Abaqus/Standard
|
||
|
||
Products: Abaqus/Standard Abaqus/CAE
|
||
|
||
# References
|
||
|
||
• “Defining ALE adaptive mesh domains in Abaqus/Standard,” Section 12.2.6
|
||
• \*ADAPTIVE MESH
|
||
• \*ADAPTIVE MESH CONSTRAINT
|
||
• \*ADAPTIVE MESH CONTROLS
|
||
• “Customizing ALE adaptive meshing,” Section 14.14 of the Abaqus/CAE User’s Guide, in the HTML version of this guide
|
||
|
||
# Overview
|
||
|
||
ALE adaptive meshing consists of two fundamental tasks:
|
||
|
||
• creating a new mesh through a process called sweeping, and
|
||
• remapping solution variables from the old mesh to the new mesh with a process called advection.
|
||
|
||
You can control the process of mesh sweeping after which, if necessary, Abaqus/Standard will automatically perform advection. The default methods for creating a new mesh have been chosen carefully to work for acoustic analysis and for modeling the effects of ablation, or wear, of material. However, you may need to override the default choices to balance the robustness and efficiency of adaptive meshing or to extend the use of adaptive meshing for other types of applications.
|
||
|
||
Adaptive mesh smoothing is defined as part of a step definition. The adaptive meshing in Abaqus/Standard uses an operator split method wherein each analysis increment consists of a Lagrangian phase followed by an Eulerian phase. The Lagrangian phase is the typical Abaqus/Standard solution increment where neither mesh sweeps nor advection occur. Once the equilibrium equations have converged, mesh smoothing is applied. Following the adjustment of nodes through the mesh sweeping process, material point quantities are advected in an Eulerian phase to account for the revised meshing of the model in its current configuration. This operator split method is chosen to avoid unsymmetric Jacobian terms that would result when the advection and material straining occur simultaneously. Advection is not required for, and is not applied to, acoustic elements.
|
||
|
||
# The ALE adaptive mesh sweeping algorithm
|
||
|
||
Adaptive mesh smoothing is performed after the structural equilibrium equations have converged. The mesh smoothing equations are solved explicitly by sweeping iteratively over the adaptive mesh domain. During each mesh sweep, nodes in the domain are relocated—based on the positions of neighboring nodes obtained during the previous mesh sweep—to reduce element distortion. The new position, $\mathbf x _ { i + 1 }$ , of a node is obtained as
|
||
|
||
<!-- source-page: 1162 -->
|
||
|
||
$$
|
||
\mathbf {x} _ {i + 1} = \mathbf {X} + \mathbf {u} _ {i + 1} = \mathbf {N} ^ {N} \mathbf {x} _ {i} ^ {N},
|
||
$$
|
||
|
||
where is the original position of the node, $\mathbf { u } _ { i + 1 }$ is the nodal displacement, $\mathbf { x } _ { i } ^ { N }$ are the neighboring nodal positions obtained during the previous mesh sweep, and ${ \bf N } ^ { N }$ are weight functions obtained from one or a weighted mixture of the following methods. The displacements applied during sweeps are not associated with mechanical behavior.
|
||
|
||
# Original configuration projection
|
||
|
||
Original configuration projection is the default in Abaqus/Standard and determines the weight function from a least squares minimization procedure that minimizes node displacement in a projection of the mesh back to the original configuration. This method of smoothing affects only deformations of the mesh and not the original mesh.
|
||
|
||
# Volume smoothing
|
||
|
||
Volume smoothing determines the weight function by computing a volume-weighted average of the element centers in the elements surrounding the node. In Figure 12.2.7–1 the new position of node M is determined by a volume-weighted average of the positions of the element centers, C, of the four surrounding elements. The volume weighting will tend to push the node away from element center C1 and toward element center C3, thus reducing element distortion.
|
||
|
||
|
||
|
||
<details>
|
||
<summary>text_image</summary>
|
||
|
||
E
|
||
L3
|
||
E
|
||
C4
|
||
×
|
||
C3
|
||
L4
|
||
M
|
||
L2
|
||
×
|
||
C2
|
||
C1
|
||
E
|
||
L1
|
||
E
|
||
</details>
|
||
|
||
Figure 12.2.7–1 Relocation of a node during a mesh sweep.
|
||
|
||
Volume smoothing is supported in structured domains, where every node is surrounded by four elements in two dimensions or eight elements in three dimensions.
|
||
|
||
<!-- source-page: 1163 -->
|
||
|
||
# Combining smoothing methods
|
||
|
||
The default smoothing method in Abaqus/Standard is original configuration projection. To choose an alternate smoothing method or to combine the smoothing methods, you specify the weighting factor for each method. When more than one smoothing method is used, a node is relocated by computing a weighted average of the locations predicted by each chosen method. All weights must be zero or positive, and their sum must be nonzero. The weights are significant only in a relative sense; their values are normalized so that their sum is 1.0.
|
||
|
||
<table><tr><td>Input File Usage:</td><td>*ADAPTIVE MESH CONTROLS, NAME=nameoriginal configuration projection weight, volume smoothing weightFor example, the following option could be used to define an equal blend of original configuration projection and volume smoothing:*ADAPTIVE MESH CONTROLS, NAME=name0.5, 0.5</td></tr></table>
|
||
|
||
<table><tr><td>Abaqus/CAE Usage:</td><td>Step module: Other→ALE Adaptive Mesh Controls→Create: Name: name, Original configuration projection: original configuration projection weight, Volumetric: volume smoothing weight</td></tr></table>
|
||
|
||
# Geometric enhancements to the basic smoothing methods
|
||
|
||
The conventional forms of the basic smoothing methods may not perform well in highly distorted domains. You can use geometrically enhanced forms of the basic smoothing algorithms as a technique to mitigate distortion. These forms are heuristic and based on nodal locations only. Due to their heuristic nature, geometric enhancements may not always improve the mesh smoothing.
|
||
|
||
Input File Usage: Use the following option to apply geometric enhancements to the smoothing algorithm:
|
||
|
||
\*ADAPTIVE MESH CONTROLS, NAME=name, GEOMETRIC ENHANCEMENT=YES
|
||
|
||
Abaqus/CAE Usage: Step module: Other→ALE Adaptive Mesh Controls→Create: Name: name, toggle on Use enhanced algorithm based on evolving element geometry
|
||
|
||
# Application of the sweeping algorithm
|
||
|
||
The mesh smoothing process begins with the mesh in its current displaced equilibrium configuration. Nodes that have no displacement degrees of freedom, such as those connected to acoustic elements, are maintained at their most recent configuration. Mesh smoothing is then driven by distortions in the current configuration and by boundary constraints. These boundary constraints can be described directly through adaptive mesh constraints. In the case of structural-acoustic boundaries the structural mesh boundary provides a constraint that controls the smoothing of adjacent acoustic element regions.
|
||
|
||
When these boundary constraints are much larger than the characteristic element length in the adaptive mesh domain, significant geometric changes, such as the development of corners, can occur.
|
||
|
||
<!-- source-page: 1164 -->
|
||
|
||
To prevent such changes, the constraints are applied gradually over a series of “sub-increments” onto the domain boundary. The number of sub-increments used is determined on the basis of the magnitude of the maximum surface displacement and the characteristic element dimension.
|
||
|
||
The remaining nodes (nodes not driven by constraints) are identified as interior nodes, free surface nodes, edge nodes, or corner nodes. These nodes are treated as described in “Defining ALE adaptive mesh domains in Abaqus/Standard,” Section 12.2.6.
|
||
|
||
At the end of mesh sweeping the new geometry is checked to ensure that elements did not become severely distorted during mesh smoothing. Abaqus/Standard responds to severe distortion in different ways, depending on the elements and procedures used. When adaptive meshing is used with acoustic elements, the current analysis increment is repeated with a reduced time increment, followed by another adaptive mesh smoothing attempt. When adaptive meshing is used with other elements, severe distortion results in abandonment of mesh smoothing for that increment. In cases where adaptive mesh constraints are also defined, Abaqus/Standard aborts since the constraints cannot be satisfied.
|
||
|
||
# Controlling the frequency of ALE adaptive mesh smoothing
|
||
|
||
In most cases the frequency of adaptive meshing is the parameter that most affects the mesh quality. By default, mesh smoothing will be performed after each converged structural analysis increment. You can change the frequency of adaptive meshing, except when adaptive mesh constraints are defined.
|
||
|
||
Input File Usage: \*ADAPTIVE MESH, FREQUENCY=number of increments
|
||
|
||
Abaqus/CAE Usage: Step module: Other→ALE Adaptive Mesh Domain→Edit: toggle on Use the adaptive mesh domain below, Frequency: number of increments
|
||
|
||
# Controlling convergence of ALE adaptive mesh smoothing
|
||
|
||
The adaptive mesh smoothing equations are solved explicitly by sweeping iteratively over the adaptive mesh domain. During each mesh sweep, nodes in the domain are relocated based on the current positions of neighboring nodes to reduce element distortion.
|
||
|
||
Mesh smoothing is performed following the end of a converged increment. You can control the intensity of the mesh smoothing by defining the number of mesh sweeps required. When the displacements are large, more iterations are usually required. When used in acoustic analyses, more iterations are usually required when the volume of the elements in the acoustic domain decreases compared to the case when the volume increases during structural loading.
|
||
|
||
You can specify the number of mesh sweeps to be performed in each adaptive mesh increment. The default number of mesh sweeps is one.
|
||
|
||
By applying the mesh sweeping algorithm repeatedly, the mesh will converge; in other words, nodal positions are obtained that do not change with further mesh sweeping. However, it is usually not necessary to apply mesh smoothing until a converged mesh is obtained; the main objective is to reduce element distortion.
|
||
|
||
Input File Usage: \*ADAPTIVE MESH, MESH SWEEPS=number of sweeps
|
||
|
||
Abaqus/CAE Usage: Step module: Other→ALE Adaptive Mesh Domain→Edit: toggle on Use the adaptive mesh domain below, Remeshing sweeps per increment: number of sweeps
|
||
|
||
<!-- source-page: 1165 -->
|
||
|
||
# The ALE adaptive mesh advection algorithm
|
||
|
||
Abaqus/Standard applies an explicit method, based on the Lax-Wendroff method, to integrate the advection equation. The key principle of the Lax-Wendroff method is replacement of the time derivatives of the material point quantities with the spatial derivatives using the classical relationship between the material time derivative, the referential derivative, and the spatial derivative. The update scheme is second-order accurate and provides some upwinding. Nodal quantities are advected by first converting them to the material point quantities.
|
||
|
||
Advection of the material quantities will generally result in loss of equilibrium, for two main reasons. The first reason is the errors in the advection process itself. To minimize the errors in advection, Abaqus/Standard imposes restrictions on the magnitude of the advection velocity by requiring that the Courant number for every element in the adaptive domain be less than one. In cases where the Courant number is greater than one you will be informed and Abaqus/Standard will generate multiple advection passes per increment. The second reason for the loss of equilibrium is changes in the representation of the underlying material quantities by the changed mesh. For example, consider a region of the structure having some stress gradients spanned initially by two elements. After mesh smoothing, the same region might have more than two elements. This will lead to slightly different volume integration while computing the internal force even when there are no errors in advection.
|
||
|
||
These sources of error in equilibrium are significant only when the mesh is too coarse to provide a good solution and mesh smoothing is carried out with such small frequency that the mesh motion is larger than the average element size. In practical applications these errors are typically insignificant, the resulting loss of equilibrium is generally small, and the residuals generated by the loss of the equilibrium fall within the limits of the Abaqus/Standard convergence criterion. Any loss of equilibrium is not propagated since equilibrium will again be satisfied at the end of the Lagrangian phase of the next increment.
|
||
|
||
# Impact of advection on subsequent steps
|
||
|
||
To ensure that the results are output only for the configuration that satisfies equilibrium, Abaqus/Standard always outputs the results at the end of the Lagrangian phase. The Eulerian phase that follows the Lagrangian phase will leave the structure out of equilibrium for the next increment. This sequence has a consequence that after the last Eulerian phase is carried out at the end of the step, equilibrium will not be satisfied exactly at the beginning of the next step and the solution at the end of the step will differ slightly from the solution at the zero increment of the following step. Equilibrium can again be established by following the step that had adaptive meshing by a step that removes all the adaptive mesh domains and allows the structure to equilibrate. A one-increment step will usually suffice. This is particularly important when the following step is a perturbation procedure that uses the solution from the previous step as the base state.
|
||
|
||
Frequency steps that follow adaptive mesh steps will also be impacted, because element mass is not advected during mesh smoothing. This impact on the element mass can be significant, depending on the extent of adaptive mesh motion and change in element size due to mesh smoothing. Abaqus will provide a warning message in cases where adaptive meshing precedes a frequency step; you should evaluate
|
||
|
||
<!-- source-page: 1166 -->
|
||
|
||
the impact of your updated mesh configuration when interpreting results from a frequency step in these cases.
|
||
|
||
# Output
|
||
|
||
In adaptive meshing the integration point of an element will generally not refer to the same material point throughout the analysis. Contour plots of material variables will show correct spatial distribution, but history plots are not meaningful. The displacement of the nodes contains the material displacement as well as the displacement due to mesh motion. You can obtain measures of the volume lost due to adaptive mesh constraints with the partial model variable VOLC, which is useful when using adaptive mesh constraints to model ablation.
|
||
|
||
A summary of the adaptive meshing in each adaptive mesh domain is written to the message (.msg) file. This summary includes the total number of load increments over which the structural displacement is transferred to the fluid, the total number of mesh sweeps performed, the magnitude of the maximum displacement increment, and the node and degree of freedom at which the maximum displacement increment is measured. Warning messages are issued when geometric features change during mesh smoothing.
|
||
|
||
More detailed diagnostic output for adaptive mesh smoothing can be requested; see “The Abaqus/Standard message file” in “Output,” Section 4.1.1. This output provides the magnitude of the maximum displacement and the node and degree of freedom where the maximum displacement increment occurs during each mesh sweep. In addition, the nodes experiencing changes in geometric features are listed.
|
||
|
||
# Additional references
|
||
|
||
• Lax, P. D., and B. Wendroff, “Difference Schemes for Hyperbolic Equations with High-Order Accuracy,” Communications on Pure and Applied Mathematics, vol. 17, p. 381, 1964.
|
||
• Lax, P. D., and B. Wendroff, “Systems of Conservations Laws,” Communications on Pure and Applied Mathematics, vol. 13, pp. 217–237, 1960.
|
||
|
||
<!-- source-page: 1167 -->
|
||
|
||
# 12.3 Adaptive remeshing
|
||
|
||
• “Adaptive remeshing: overview,” Section 12.3.1
|
||
• “Selection of error indicators influencing adaptive remeshing,” Section 12.3.2
|
||
• “Solution-based mesh sizing,” Section 12.3.3
|
||
|
||
<!-- source-page: 1168 -->
|
||
|
||
<!-- source-page: 1169 -->
|
||
|
||
# 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:
|
||
|
||
• You identify regions of the model where you wish to apply one or more adaptive remeshing rules. A remeshing rule defines the step during which it will be applied, the error indicator output variables and targets for those error indicators, the sizing method, and any constraints on element size. See “What are remeshing rules?,” Section 17.13.1 of the Abaqus/CAE User’s Guide.
|
||
• You define a succession of analysis jobs, an “adaptivity process,” that will be run as Abaqus/CAE attempts to meet your remeshing rule targets. See “What is an adaptivity process?,” Section 19.3.1 of the Abaqus/CAE User’s Guide.
|
||
|
||
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 “Selection of error indicators influencing adaptive remeshing,” 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 19.3.2 of the Abaqus/CAE User’s Guide, for more details. Figure 12.3.1–1 shows the interaction of Abaqus products and files in this 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
|
||
|
||
<!-- source-page: 1170 -->
|
||
|
||
|
||
|
||
<details>
|
||
<summary>flowchart</summary>
|
||
|
||
```mermaid
|
||
graph TD
|
||
A["Create model"] --> B["Create mesh"]
|
||
B --> C["Create remeshing rules"]
|
||
C --> D["Create adaptivity process"]
|
||
D --> E["Submit adaptivity process"]
|
||
E --> F["Submit analysis job"]
|
||
F --> G{Are meshing rules satisfied?}
|
||
G -->|Yes| H["Adaptivity process is finished"]
|
||
G -->|No| I["Compute new element sizes"]
|
||
I --> J["Create mesh"]
|
||
H --> K["Review Results"]
|
||
J --> K
|
||
style A fill:#d4edda,stroke:#333
|
||
style K fill:#d4edda,stroke:#333
|
||
subgraph Setup Analysis
|
||
A
|
||
B
|
||
C
|
||
D
|
||
E
|
||
F
|
||
end
|
||
subgraph Monitor Progress
|
||
G
|
||
I
|
||
J
|
||
end
|
||
subgraph Perform Analysis
|
||
F
|
||
G
|
||
H
|
||
I
|
||
end
|
||
style G fill:#d4edda,stroke:#333
|
||
style H fill:#d4edda,stroke:#333
|
||
style I fill:#d4edda,stroke:#333
|
||
style J fill:#d4edda,stroke:#333
|
||
style K fill:#d4edda,stroke:#333
|
||
style_L["User Actions"] --> A
|
||
M["Automated Abaqus/CAE Actions"] --> H
|
||
```
|
||
</details>
|
||
|
||
Figure 12.3.1–1 User actions and automated Abaqus/CAE actions in the adaptive remeshing process.
|
||
|
||
• 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 Guide. The example includes a Python script that you can run from
|