27.1.4 Section controls
Products: Abaqus/Standard Abaqus/Explicit Abaqus/CAE
References
Overview
Section controls in Abaqus/Standard:
choose the hourglass control formulation for most first-order elements with reduced integration;
define the distortion control for C3D10HS elements;
select the hourglass control scale factors for all elements with reduced integration; and
select the choice of element deletion and the value of maximum degradation for cohesive elements, connector elements, elements with plane stress formulations (plane stress, shell, continuum shell, and membrane elements) with constitutive behavior that includes damage evolution, any element that can be used with damage evolution models for ductile metals, and any element that can be used with the damage evolution law in a low-cycle fatigue analysis.
choose the hourglass control formulation or scale factors for all elements with reduced integration;
define the distortion control for solid elements;
select the scale factors for the drill stiffness of shell elements or deactivate the drill stiffness for small-strain shell elements S3RS and S4RS;
select an amplitude for ramping of any initial stresses in membrane elements;
select the kinematic formulation for hexahedron solid elements;
select the accuracy order of the formulation for solid and shell elements;
select the scale factors for linear and quadratic bulk viscosity parameters;
specify the size of the particle tracking box for discrete element method (DEM) analyses and smoothed particle hydrodynamic (SPH) analyses;
specify the formulation and additional control parameters for SPH analyses; and
select the choice of element deletion and the value of maximum degradation for elements with constitutive behavior that includes damage evolution.
Using section controls
In Abaqus/Standard section controls are used to select the enhanced hourglass control formulation for solid, shell, and membrane elements. This formulation provides improved coarse mesh accuracy with slightly higher computational cost and performs better for nonlinear material response at high strain levels when compared with the default total stiffness formulation. Section controls can also be used to select some element formulations that may be relevant for a subsequent Abaqus/Explicit analysis.
In Abaqus/Explicit the default formulations for solid, shell, and membrane elements have been chosen to perform satisfactorily on a wide class of quasi-static and explicit dynamic simulations. However, certain formulations give rise to some trade-off between accuracy and performance. Abaqus/Explicit provides section controls to modify these element formulations so that you can optimize these objectives for a specific application. Section controls can also be used in Abaqus/Explicit to specify scale factors for linear and quadratic bulk viscosity parameters. You can also control the initial stresses in membrane elements for applications such as airbags in crash simulations and introduce the initial stresses gradually based on an amplitude definition.
In addition, section controls are used to specify the maximum stiffness degradation and to choose the behavior upon complete failure of an element, once the material stiffness is fully degraded, including the removal of failed elements from the mesh. This functionality applies only to elements with a material definition that includes progressive damage (see “Progressive damage and failure,” Section 24.1.1; “Connector damage behavior,” Section 31.2.7; and “Defining the constitutive response of cohesive elements using a traction-separation description,” Section 32.5.6). In Abaqus/Standard this functionality is limited to
cohesive elements with a traction-separation constitutive response that includes damage evolution,
any element with a plane stress formulation that can be used with the damage evolution model for fiber-reinforced composites,
any element that can be used with the damage evolution models for ductile metals,
any element that can be used with the damage evolution law in a low-cycle fatigue analysis, and
connector elements with a constitutive response that includes damage evolution.
| Input File Usage: | Use the following option to specify a section controls definition: |
*SECTION CONTROLS, NAME=name This option is used in conjunction with one or more of the following options to associate the section control definition with an element section definition: *COHESIVE SECTION, CONTROLS=name *CONNECTOR SECTION, CONTROLS=name *DISCRETE SECTION, CONTROLS=name *EULERIAN SECTION, CONTROLS=name *MEMBRANE SECTION, CONTROLS=name *SHELL GENERAL SECTION, CONTROLS=name *SHELL SECTION, CONTROLS=name *SOLID SECTION, CONTROLS=name You can apply a single section control definition to several element section definitions. |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Methods for suppressing hourglass modes
The formulation for reduced-integration elements considers only the linearly varying part of the incremental displacement field in the element for the calculation of the increment of physical strain. The remaining part of the nodal incremental displacement field is the hourglass field and can be expressed in terms of hourglass modes.
Excitation of these modes may lead to severe mesh distortion, with no stresses resisting the deformation. Similarly, the formulation for element type C3D4H considers in the constraint equations only the constant part of the incremental pressure Lagrange multiplier field. The remaining part of the nodal incremental pressure Lagrange multiplier interpolation is comprised of hourglass modes.
Hourglass control attempts to minimize these problems without introducing excessive constraints on the element's physical response.
Several methods are available in Abaqus for suppressing the hourglass modes, as described below.
Integral viscoelastic approach in Abaqus/Explicit
The integral viscoelastic approach available in Abaqus/Explicit generates more resistance to hourglass forces early in the analysis step where sudden dynamic loading is more probable.
Let q be an hourglass mode magnitude and Q be the force (or moment) conjugate to q. The integral viscoelastic approach is defined as
, 
, and 
that you can define (by default, 
). The scale factors are dimensionless and relate to specific displacement degrees of freedom. For solid and membrane elements scales all hourglass stiffnesses. For shell elements scales the hourglass stiffnesses related to the in-plane displacement degrees of freedom, and scales the hourglass stiffnesses related to the rotational degrees of freedom. In addition, scales the hourglass stiffness related to the transverse displacement for small-strain shell elements.The integral viscoelastic form of hourglass control is available for all reduced-integration elements and is the default form in Abaqus/Explicit, except for elements modeled with hyperelastic, hyperfoam, and low-density foam materials. It is the most computationally intensive hourglass control method. It is not supported for Eulerian EC3D8R elements.
| Input File Usage: | *SECTION CONTROLS, NAME=name, HOURGLASS=RELAX STIFFNESS |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Kelvin viscoelastic approach in Abaqus/Explicit
The Kelvin-type viscoelastic approach available in Abaqus/Explicit is defined as
Three approaches are provided in Abaqus/Explicit for specifying Kelvin viscoelastic hourglass control.
Specifying the pure stiffness approach
The pure stiffness form of hourglass control is available for all reduced-integration elements and is recommended for both quasi-static and transient dynamic simulations.
| Input File Usage: | *SECTION CONTROLS, NAME=name, HOURGLASS=STIFFNESS |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Specifying the pure viscous approach
The pure viscous form of hourglass control is available only for solid and membrane elements with reduced integration and is the default form in Abaqus/Explicit for Eulerian EC3D8R elements. It is the most computationally efficient form of hourglass control and has been shown to be effective for high-rate dynamic simulations. However, the pure viscous method is not recommended for low frequency dynamic or quasi-static problems since continuous (static) loading in hourglass modes will result in excessive hourglass deformation due to the lack of any nominal stiffness.
| Input File Usage: | *SECTION CONTROLS, NAME=name, HOURGLASS=VISCOUS |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Specifying a combination of stiffness and viscous hourglass control
A linear combination of stiffness and viscous hourglass control is available only for solid and membrane elements with reduced integration. You can specify the blending weight factor 
(
) to scale the stiffness and viscous contributions. Specifying a weight factor equal to 0.0 or 1.0 results in the limiting cases of pure stiffness and pure viscous hourglass control, respectively. The default weight factor is 0.5.
| Input File Usage: | *SECTION CONTROLS, NAME=name, HOURGLASS=COMBINED, WEIGHT FACTOR= |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Total stiffness approach in Abaqus/Standard
The total stiffness approach available in Abaqus/Standard is the default hourglass control approach for all first-order, reduced-integration elements in Abaqus/Standard, except for elements modeled with hyperelastic, hyperfoam, or hysteresis materials. It is the only hourglass control approach available in Abaqus/Standard for S8R5, S9R5, and M3D9R elements and the only hourglass control approach available for the pressure Lagrange multiplier interpolation for C3D4H elements. Hourglass stiffness factors for first-order, reduced-integration elements depend on the shear modulus, while factors for C3D4H elements depend on the bulk modulus. A scale factor can be applied to these stiffness factors to increase or decrease the hourglass stiffness.
Let q be an hourglass mode magnitude and Q be the force (moment, pressure, or volumetric flux) conjugate to q. The total stiffness approach for hourglass control in membrane or solid elements or membrane hourglass control in shell elements is defined as
is a dimensionless scale factor (by default, 
); 
is an hourglass stiffness factor with units of stress; 
is the gradient interpolator used to define constant gradients in the element (
where the superscript P refers to an element node, the subscript refers to a direction, and 
is a material coordinate); and V is the element volume. Similarly, the hourglass control for the pressure Lagrange multiplier interpolation for C3D4H elements is defined as
is a dimensionless scale factor (by default, 
); 
is a volumetric gradient operator; and 
is an hourglass stiffness factor with units of stress for compressible hyperelastic and hyperfoam materials and units of stress compliance for all other materials. The total stiffness approach for bending hourglass control in shell elements is defined as 
is the scale factor (by default, 
), 
is the hourglass stiffness factor, t is the thickness of the shell element, and A is the area of the shell element.| Input File Usage: | *SECTION CONTROLS, NAME=name, HOURGLASS=STIFFNESS |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Default hourglass stiffness values
Normally the hourglass control stiffness is defined from the elasticity associated with the material. In most cases, the control stiffness of first-order, reduced-integration elements is based on a typical value of the initial shear modulus of the material, which may, for example, be given as part of the elastic material definition (“Linear elastic behavior,” Section 22.2.1). Similarly, hourglass control stiffness of the reduced-integration pressure and volumetric Lagrange multiplier interpolations of C3D4H elements is based on a typical value of the initial bulk modulus. For an isotropic elastic or hyperelastic material G is the shear modulus. For a nonisotropic elastic material average moduli are used to calculate the hourglass stiffness: for orthotropic elasticity defined by specifying the terms in the elastic stiffness matrix or for anisotropic elasticity
If the elastic moduli are dependent on temperature or field variables, the first value in the table is used. The default values for the stiffness factors are defined below.
For membrane or solid elements
is the section stiffness matrix.User-defined hourglass stiffness
When the initial shear modulus is not defined, you must define the hourglass stiffness parameters; an example is when user subroutine UMAT is used to describe the material behavior of elements with hourglassing modes. In some cases the default value provided for the hourglass control stiffness may not be suitable and you should define the value.
In some coupled pore fluid diffusion and stress analyses the prevailing pore pressure in the medium may approach the magnitude of the stiffness of the material skeleton, as measured by constitutive parameters such as the elastic modulus. These cases are expected in some partial saturation evaluations of the wetting of relatively compliant materials such as tissues or cloth. When reduced-integration or modified tetrahedral or triangular elements are used in such analyses, the default choice of the hourglass control stiffness parameter, which is based on a scaling of skeleton material constitutive parameters, may not be adequate to control hourglassing in the presence of large pore pressure fields. An appropriate hourglass control stiffness in these cases should scale with the expected magnitude of pore pressure changes over an element.
| Input File Usage: | Use the following option to specify nondefault values for the hourglass stiffness factors: |
*HOURGLASS STIFFNESS This option must immediately follow one of the following options: *MEMBRANE SECTION *SHELL GENERAL SECTION *SHELL SECTION *SOLID SECTION |
, 
, 
, drilling hourglass scaling factor for shells| Abaqus/CAE Usage: | Mesh module: Mesh |
Enhanced hourglass control approach in Abaqus/Standard and Abaqus/Explicit
The enhanced hourglass control approach available in both Abaqus/Standard and Abaqus/Explicit represents a refinement of the pure stiffness method in which the stiffness coefficients are based on the enhanced assumed strain method; no scale factor is required. It is the default hourglass control approach for hyperelastic, hyperfoam, and low-density foam materials in Abaqus/Explicit and for hyperelastic, hyperfoam, and hysteresis materials in Abaqus/Standard. This method gives more accurate displacement solutions for coarse meshes with linear elastic materials as compared to other hourglass control methods. It also provides increased resistance to hourglassing for nonlinear materials. Although generally beneficial, this may give overly stiff response in problems displaying plastic yielding under bending. In Abaqus/Explicit the enhanced hourglass method will generally predict a much better return to the original configuration for hyperelastic or hyperfoam materials when loading is removed.
The enhanced hourglass control approach is compatible between Abaqus/Standard and Abaqus/Explicit. It is recommended that enhanced hourglass control be used for both Abaqus/Standard and Abaqus/Explicit for all import analyses. See “Transferring results between Abaqus/Explicit and Abaqus/Standard,” Section 9.2.2.
The enhanced hourglass method is not supported for enriched elements (see “Modeling discontinuities as an enriched feature using the extended finite element method,” Section 10.7.1).
Specifying the enhanced hourglass control approach
The enhanced hourglass control method is available for first-order solid, membrane, and finite-strain shell elements with reduced integration. In Abaqus/Explicit it cannot be used for a hyperelastic or hyperfoam material when adaptive meshing is used on that domain (see the discussion below).
| Input File Usage: | *SECTION CONTROLS, NAME=name, HOURGLASS=ENHANCED |
Any scaling factors specified on the data line following this option will be ignored. |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Special considerations for hyperelastic and hyperfoam materials in an adaptive mesh domain in Abaqus/Explicit
The enhanced hourglass method cannot be used with elements modeled with hyperelastic or hyperfoam materials that are included in an adaptive mesh domain. Thus, if you decide to use hyperelastic or hyperfoam materials in an adaptive mesh domain, you must specify section controls to choose a different hourglass control approach. The use of adaptive meshing in domains modeled with finite-strain elastic materials is not recommended since better results are generally predicted using the enhanced hourglass method and, for solid elements, element distortion control (discussed below). Therefore, for these materials it is recommended that the analysis be run without adaptive meshing but with enhanced hourglass control.
Use in coupled pore pressure analysis
When first-order, reduced-integration, or modified tetrahedral or triangular elements are used in coupled pore fluid diffusion and stress analyses or coupled temperature–pore pressure analyses with enhanced hourglass control, the hourglass control stiffness, which is based on skeleton material constitutive parameters, may not be adequate to control hourglassing in the presence of large pore pressure fields. Since enhanced hourglass control does not allow you to change the hourglass control stiffness, it is recommended that total stiffness hourglass control be used in these cases with an appropriate hourglass control stiffness scaled with the expected magnitude of pore pressure changes over an element.
Controlling element distortion for crushable materials in Abaqus/Explicit
Many analyses with volumetrically compacting materials such as crushable foams see large compressive and shear deformations, especially when the crushable materials are used as energy absorbers between stiff or heavy components. The material behavior for crushable materials usually stiffens significantly under high compression. When a finer mesh is used, the stiffening behavior of the material model is enough to prevent excessive negative element volumes or other excessive distortion from occurring under high compressive loads. However, analyses may fail prematurely when the mesh is coarse relative to strain gradients and the amount of compression.
Abaqus/Explicit offers distortion control to prevent solid elements from inverting or distorting excessively for these cases. The constraint method used in Abaqus/Explicit prevents each node on an element from punching inward toward the center of the element past a point where the element would become non-convex. Constraints are enforced by using a penalty approach, and you can control the associated distortion length ratio.
Distortion control is available only for solid elements and cannot be used when the elements are included in an adaptive mesh domain. Distortion control is activated by default for elements modeled with hyperelastic, hyperfoam, or low-density foam materials. Using adaptive meshing in a domain modeled with hyperelastic or hyperfoam materials is not recommended since better results are generally predicted using the enhanced hourglass method in combination with element distortion control. However, if you decide to use hyperelastic or hyperfoam materials in an adaptive mesh domain, you must specify section controls to deactivate distortion control.
When element distortion control is used in combination with the enhanced hourglass method (default behavior for hyperelastic and hyperfoam materials), a small amount of viscous damping is added to the element formulation and the associated viscous energy dissipation is included in the output of artificial strain energy (ALLAE).
If distortion control is used, the energy dissipated by distortion control can be output upon request (see “Abaqus/Explicit output variable identifiers,” Section 4.2.2, for details). Although developed for analyses of energy absorbing, volumetrically compacting materials, distortion control can be used with any material model. However, care must be used in interpreting results since the distortion control constraints may inhibit legitimate deformation modes and lock up the mesh. Distortion control cannot prevent elements from being distorted due to temporal instabilities, hourglass instabilities, or physically unrealistic deformation.
| Input File Usage: | Use the following option to activate distortion control: |
*SECTION CONTROLS, NAME=name, DISTORTION CONTROL=YES Use the following option to deactivate distortion control: *SECTION CONTROLS, NAME=name, DISTORTION CONTROL=NO |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Controlling the distortion length ratio
By default, the constraint penalty forces are applied when the node moves to a point a small offset distance away from the actual plane of constraint. This appears to improve the robustness of the method and limits the reduction of time increment due to severe shortening of the element characteristic length. This offset distance is determined by the distortion length ratio times the initial element characteristic length. The default value of the distortion length ratio, r, is 0.1. You can change the distortion length ratio by specifying a value for r, 
.
| Input File Usage: | *SECTION CONTROLS, NAME=name, DISTORTION CONTROL=YES, LENGTH RATIO=r |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Selecting a scale factor for the drill stiffness in Abaqus/Explicit
A drill constraint acts to keep the element nodal rotations in the direction of the shell normal consistent with the average in-plane rotation of the element. Lack of such a constraint can lead to large rotations at these element nodes. Section controls can be used to select a scale factor for the default drill stiffness of an individual element set.
| Input File Usage: | Use the following options to specify a scale factor for the drill stiffness: |
*SECTION CONTROLS, NAME=name , , , , , , , scale factor for drill stiffness |
Drill constraint in small strain shell elements S3RS and S4RS in Abaqus/Explicit
The formulation of small strain shell elements S3RS and S4RS includes a drill constraint and does so by default. Alternatively, you can deactivate the drill constraint for these elements. The drill constraint is always active for the finite strain conventional shell elements such as S4R, but the default value of the drill stiffness can be scaled as mentioned above.
| Input File Usage: | Use the following option to activate the drill constraint (default): |
*SECTION CONTROLS, DRILL STIFFNESS=ON Use the following option to deactivate the drill constraint: *SECTION CONTROLS, DRILL STIFFNESS=OFF |
Ramping of initial stresses in membrane elements in Abaqus/Explicit
For applications such as airbags in crash simulations the initial strains (hence, the initial stresses) are introduced into the model through a reference configuration that is different from the initial configuration. Often the components that confine the airbag in the initial configuration are excluded from the numerical model causing motion of the airbag under initial stresses at the beginning of the analysis. Abaqus/Explicit provides a technique to introduce the initial stresses in the membrane elements gradually based on an amplitude definition. This amplitude must be defined with its value starting from zero and reaching a final value of one. The initial stresses will not be applied for the duration that the amplitude stays at zero.
| Input File Usage: | Use both of the following options: |
*AMPLITUDE, NAME=name *SECTION CONTROLS, RAMP INITIAL STRESS=name |
Defining the kinematic formulation for hexahedron solid elements
The default kinematic formulation for reduced-integration solid elements in Abaqus (and the only kinematic formulation available in Abaqus/Standard) is based on the uniform strain operator and the hourglass shape vectors. Details can be found in “Solid isoparametric quadrilaterals and hexahedra,” Section 3.2.4 of the Abaqus Theory Guide. These kinematic assumptions result in elements that pass the constant strain patch test for a general configuration and give zero strain under large rigid body rotation. However, the formulation is relatively expensive, especially in three dimensions.
Abaqus/Explicit offers two alternative kinematic formulations for the C3D8R solid element that can reduce the computational cost. The performance for each kinematic formulation on the patch test and under large rigid body rotation for various element configurations is summarized in Table 27.1.4–1. Suitable applications for each kinematic formulation are summarized in Table 27.1.4–2.
Table 27.1.4–1 Element performance for patch test and large rigid body rotations for various element configurations.
| Element configuration | Kinematic formulation type | |||
|---|---|---|---|---|
| Average strain | Orthogonal | Centroid | ||
| Satisfaction of the three-dimensional patch test | Parallelepiped | Yes | Yes | Yes |
| General | Yes | No | No | |
| Zero straining under rigid body rotation | Parallelepiped | Yes | Yes | Yes |
| General | Yes | Yes | No | |
Table 27.1.4–2 Different element formulations and their suitable applications. The default formulation is highlighted below.
| Kinematic formulation | Order of accuracy | Suitable applications |
|---|---|---|
| Average strain | Second-order | All; recommended for problems involving a large number of revolutions (>5). |
| Average strain | First-order | All; except those involving a large number of revolutions (>5). |
| Orthogonal | — | All; except those involving high confinement, very coarse meshes, or highly distorted elements. |
| Centroid | — | Problems with little rigid body rotation and reasonable mesh refinement. |
You can specify the kinematic formulation for 8-node brick elements.
Default formulation
The default average strain formulation of uniform strain and hourglass shape vectors is the only formulation available in Abaqus/Standard. This formulation is recommended for all problems and is particularly well suited for applications exhibiting high confinement, such as closed-die forming and bushing analyses.
| Input File Usage: | *SECTION CONTROLS, KINEMATIC SPLIT=AVERAGE STRAIN |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Orthogonal formulation in Abaqus/Explicit
A noticeable reduction in computational cost can be obtained by using the orthogonal formulation available in Abaqus/Explicit. This formulation is based on the centroidal strain operator and a slight modification to the hourglass shape vectors. The centroidal strain operator requires three times fewer floating point operations than the uniform strain operator. Elements formulated with an orthogonal kinematic split pass the patch test only for rectangular or parallelepiped element configurations. However, numerical experience has shown that the element converges on the exact solution for general element configurations as the mesh is refined. It also performs well for large rigid body motions.
This formulation provides a good balance between computational speed and accuracy. It is recommended for all analyses except those involving highly distorted elements, very coarse meshes, or high confinement. Suitable applications for this formulation include elastic drop testing.
| Input File Usage: | *SECTION CONTROLS, NAME=name, KINEMATIC SPLIT=ORTHOGONAL |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Centroid formulation in Abaqus/Explicit
The fastest formulation available in Abaqus/Explicit is specified by selecting the centroid formulation. The centroid formulation is based on the centroidal strain operator and the hourglass base vectors. Using the hourglass base vectors instead of the hourglass shape vectors reduces hourglass mode computations by a factor of three. However, the hourglass base vectors are not orthogonal to rigid body rotation for general element configurations, so that hourglass strain may be generated with large rigid body rotations with this formulation.
This formulation should be used only to improve computational performance on problems that have reasonable mesh refinement and no significant amount of rigid body rotation (e.g., transient flat rolling simulation).
| Input File Usage: | *SECTION CONTROLS, NAME=name, KINEMATIC SPLIT=CENTROID |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Choosing the order of accuracy in solid and shell element formulations
Abaqus/Standard offers only a second-order accurate formulation for all elements.
Abaqus/Explicit offers both first- and second-order accurate formulations for solid and shell elements. First-order accuracy is the default and yields sufficient accuracy for nearly all Abaqus/Explicit problems because of the inherently small time increment size. Second-order accuracy is usually required for analyses with components undergoing a large number of revolutions (>5). For three-dimensional solids the second-order accuracy formulation is available only with the default average strain kinematic formulation.
First-order accuracy
In Abaqus/Explicit the first-order accurate formulation for solid and shell elements is the default. This formulation is not available in Abaqus/Standard.
| Input File Usage: | *SECTION CONTROLS, NAME=name, SECOND ORDER ACCURACY=NO |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Second-order accuracy
The second-order accurate element formulation is appropriate for problems with a large number of revolutions (>5). This is the only formulation available in Abaqus/Standard. “Simulation of propeller rotation,” Section 2.3.15 of the Abaqus Benchmarks Guide, illustrates the performance of second-order accurate shell and solid elements in Abaqus/Explicit as they undergo about 100 revolutions.
| Input File Usage: | *SECTION CONTROLS, NAME=name, SECOND ORDER ACCURACY=YES |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Selecting scale factors for bulk viscosity in Abaqus/Explicit
Bulk viscosity introduces damping associated with volumetric straining. Its purpose is to improve the modeling of high-speed dynamic events. Abaqus/Explicit contains two forms of bulk viscosity, linear and quadratic, which can be defined for the whole model at each step of the analysis, as discussed in “Bulk viscosity” in “Explicit dynamic analysis,” Section 6.3.3. Section controls can be used to select scale factors for the linear and quadratic bulk viscosities of an individual element set.
The pressure term generated by bulk viscosity may introduce unexpected results in the volumetric response of highly compressible materials; therefore, it is recommended to suppress bulk viscosity for these materials by specifying scale factors equal to zero.
| Input File Usage: | Use the following options to specify scale factors for the linear and quadratic bulk viscosities: |
*SECTION CONTROLS, NAME=name , , , scale factor for linear bulk viscosity, scale factor for quadratic bulk viscosity |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Controlling element deletion and maximum degradation for materials with damage evolution
Abaqus offers a general capability for modeling progressive damage and failure of materials (“Progressive damage and failure,” Section 24.1.1). In Abaqus/Standard this capability is available only for cohesive elements, connector elements, elements with plane stress formulations (plane stress, shell, continuum shell, and membrane elements), any element that can be used with the damage evolution models for ductile metals, and any element that can be used with the damage evolution law in a low-cycle fatigue analysis. In Abaqus/Explicit this capability is available for all elements with progressive damage behavior except connector elements. Section controls are provided to specify the value of the maximum stiffness degradation, 
, and whether element deletion occurs when the degradation reaches this level. By default, an element is deleted when it is fully damaged (i.e., 
). The choice of element deletion also affects how the damage is applied; details can be found in the following sections:
| Input File Usage: | Use the following option to delete the element from the mesh: |
*SECTION CONTROLS, ELEMENT DELETION=YES Use the following option to keep the element in the computation: *SECTION CONTROLS, ELEMENT DELETION=NO Use the following option to specify *SECTION CONTROLS, MAX DEGRADATION= |
| Abaqus/CAE Usage: | Use the following option to control whether completely damaged elements remain in the computation: |
Mesh module: Mesh Use the following option to determine when an element is considered completely damaged: Mesh module: Mesh |
Using viscous regularization with cohesive elements, connector elements, and elements that can be used with the damage evolution models for ductile metals and fiber-reinforced composites in Abaqus/Standard
Material models exhibiting softening behavior and stiffness degradation often lead to severe convergence difficulties in implicit analysis programs, such as Abaqus/Standard. A common technique to overcome some of these convergence difficulties is the use of viscous regularization of the constitutive equations, which causes the tangent stiffness matrix of the softening material to be positive for sufficiently small time increments.
The traction-separation laws used to describe the constitutive behavior of cohesive elements can be regularized in Abaqus/Standard using viscosity, by permitting stresses to be outside the limits defined by the traction-separation law. The details of the regularization procedure are discussed in “Viscous regularization in Abaqus/Standard” in “Defining the constitutive response of cohesive elements using a traction-separation description,” Section 32.5.6. The same technique is also used to regularize the following:
damaged (softening) connector response (see “Connector damage behavior,” Section 31.2.7),
damaged response of elements with plane stress formulations when they are used with the damage model for fiber-reinforced materials (see “Viscous regularization” in “Damage evolution and element removal for fiber-reinforced composites,” Section 24.3.3), and
damage response of elements used with the damage model for ductile metals (see “Damage evolution and element removal for ductile metals,” Section 24.2.3).
| Input File Usage: | *SECTION CONTROLS, VISCOSITY= |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Using viscous damping with connector elements in Abaqus/Standard
Material failure in connector elements often causes convergence problems in Abaqus/Standard. To avoid such convergence problems, you can introduce viscous damping into the connector components by specifying the value of the damping coefficient as discussed in “Connector failure behavior,” Section 31.2.9. By default, no damping is included.
| Input File Usage: | *SECTION CONTROLS, VISCOSITY= |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Using section controls in an import analysis
The recommended procedure for doing import analysis is to specify the enhanced hourglass control formulation in the original analysis. Once the section controls have been specified in the original analysis, they cannot be modified in subsequent import analyses. This ensures that the enhanced hourglass control formulation is used in the original as well as import analyses. The default values for other section controls are usually appropriate and should not be changed. For further details on using section controls in an import analysis, see “Transferring results between Abaqus/Explicit and Abaqus/Standard,” Section 9.2.2.
Using section controls for flexion-torsion type connector
When the third axes of the two local coordinate systems for a flexion-torsion type connector are exactly aligned, a numerical singularity occurs that may lead to convergence difficulties. To avoid this, a small perturbation can be applied to the local coordinate system defined at the second connector node.
| Input File Usage: | *SECTION CONTROLS, PERTURBATION=small angle |
| Abaqus/CAE Usage: | You cannot specify a perturbation for flexion-torsion type connectors in Abaqus/CAE. |
Using section controls to define the particle tracking box for DEM and SPH particles
For discrete element method (DEM) analyses, a particle tracking box is established at the beginning of the analysis to define the rectangular region within which the particle search (finding all neighbors for all particles) is performed. A region that is 10% larger in all directions than the overall model initial dimensions and is centered at the geometric center of the model is used.
For smoothed particle hydrodynamic (SPH) analyses, all particles are tracked as the analysis progresses by default. For DEM analyses, particle tracking is based on the initially established tracking box by default. Alternatively, you can define a particle tracking box to define the region within which the particle search is performed.
You define a fixed size for the particle tracking box by specifying the coordinates of two opposite corners (lower left and upper right) of this box. As the analysis progresses, if a particle is outside this tracking box, it behaves like a free-flying point mass and does not contribute to the DEM or SPH calculations. If the particle reenters the box at a later stage, it is once again included in the calculations. If you want to track all of the particles during the analysis, you must ensure that the particle tracking box fully encompasses the domain through which the model moves; otherwise, you will lose tracking of the particle.
| Input File Usage: | Use the following option to specify a fixed size for the particle tracking box in a DEM analysis: |
*SECTION CONTROLS blank line blank line X, Y, and Z-coordinates (lower box corner) and X, Y, and Z-coordinates Use the following option to specify a fixed size for the particle tracking box in an SPH analysis: *SECTION CONTROLS first data line second data line X, Y, and Z-coordinates (lower box corner) and X, Y, and Z-coordinates |
| Abaqus/CAE Usage: | In Abaqus/CAE you can only specify section controls for SPH parameters in Abaqus/Explicit analyses involving the conversion of continuum elements to SPH particles. |
Using section controls for smoothed particle hydrodynamics (SPH)
In addition to controlling the size of the particle tracking box, you can control other aspects of the smoothed particle hydrodynamic (SPH) formulation implemented in Abaqus/Explicit.
Using section controls for specifying the SPH kernel
For a smoothed particle hydrodynamic analysis, you can choose the order of the kernel used for interpolation. For a list of references that discuss the various kernels that can be used, see “Smoothed particle hydrodynamics,” Section 15.2.1.
| Input File Usage: | Use one of the following options: |
*SECTION CONTROLS, KERNEL=CUBIC *SECTION CONTROLS, KERNEL=QUADRATIC *SECTION CONTROLS, KERNEL=QUINTIC |
| Abaqus/CAE Usage: | In Abaqus/CAE you can choose the order of the kernel used for interpolation only in Abaqus/Explicit analyses involving the conversion of continuum elements to SPH particles. |
Mesh module: Mesh |
Using section controls for specifying the SPH formulation
By default, the SPH kernels satisfy the zero-order completeness requirement. A first-order complete corrected (normalized) kernel is also available, which is sometimes referred in the literature as the normalized SPH (NSPH) method. In high-deformation solid mechanics analyses the use of the NSPH method may lead to more accurate results.
In the SPH methods, a mean velocity filtering coefficient can be used for the modified coordinate updates for particles. A nonzero value for this coefficient leads to the XSPH method, as discussed in “Smoothed particle hydrodynamics,” Section 15.2.1.
| Input File Usage: | Use one of the following options to specify the SPH formulation: |
*SECTION CONTROLS, SPH FORMULATION=CLASSICAL (default) *SECTION CONTROLS, SPH FORMULATION=NSPH *SECTION CONTROLS, SPH FORMULATION=XSPH |
| Abaqus/CAE Usage: | In Abaqus/CAE you can only specify section controls for SPH parameters in Abaqus/Explicit analyses involving the conversion of continuum elements to SPH particles. |
Using section controls for specifying SPH parameters
You can control the way the smoothing length is computed (see “Smoothed particle hydrodynamics,” Section 15.2.1). You can specify the smoothing length (units of length) for precise control of the radius of influence associated with a given particle. Alternatively, you can scale the default smoothing length by specifying a dimensionless smoothing length factor. By default, the smoothing length is kept constant throughout the analysis. You can specify a variable smoothing length that will increase or decrease during the analysis depending on the divergence of the velocity field, which is a measure of compressive or expansive behavior.
You can also specify the minimum number of particles within the sphere of influence for the given particle. If the total number of particles within the sphere of influence for the given particle is less than the specified minimum number of particles, the deformation gradient for this given particle is frozen, that is, unchanged between the previous and current time increment. In solid mechanics it means that the strain associated with this element will not be changed during the current time increment.
You can specify a mean velocity filtering coefficient that is used for the modified coordinate updates for particles using the XSPH method.
| Input File Usage: | Use the following option to specify SPH parameters: |
*SECTION CONTROLS first data line smoothing length, smoothing length factor, min number of neighboring particles, , mean velocity filtering coefficient Use one of the following options to define the smoothing length: *SECTION CONTROLS, SPH SMOOTHING LENGTH=CONSTANT (default) *SECTION CONTROLS, SPH SMOOTHING LENGTH=VARIABLE |
| Abaqus/CAE Usage: | In Abaqus/CAE you can only specify section controls for SPH parameters in Abaqus/Explicit analyses involving the conversion of continuum elements to SPH particles. |
Using section controls to convert continuum elements to particles
Reduced-integration continuum elements can convert to particles if a certain criterion is met, as discussed in “Finite element conversion to SPH particles,” Section 15.2.2. You can specify the number of particles per parent element to be generated. Several criteria to trigger the conversion are available.
| Input File Usage: | Use the following option to prevent finite elements from converting to particles: |
*SECTION CONTROLS, ELEMENT CONVERSION=NO (default) Use the following option to trigger the conversion of finite elements to particles: *SECTION CONTROLS, ELEMENT CONVERSION=YES Use the following option to trigger the conversion of finite elements to particles based on a uniform background grid: *SECTION CONTROLS, ELEMENT CONVERSION=BACKGROUND GRID |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Generating particles based on a uniform background grid is not supported in Abaqus/CAE. |
Specifying the number of particles generated
You specify the number of particles to be generated per isoparametric direction. The number of particles can range from 1 to 7.
| Input File Usage: | *SECTION CONTROLS, ELEMENT CONVERSION=YES first data line second data line third data line number of particles to be generated per isoparametric direction |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Specifying the background grid
You specify the spacing of the background grid and the name of an orientation definition to define a local coordinate system for the background grid.
| Input File Usage: | *SECTION CONTROLS, ELEMENT CONVERSION=BACKGROUND GRID first data line second data line third data line spacing of the background grid, name of an orientation definition |
| Abaqus/CAE Usage: | Generating particles based on a uniform background grid is not supported in Abaqus/CAE. |
Specifying the thickness of generated particles
The thickness of the particles is primarily used in resolving initial overclosures between the particles and the surfaces in the general contact. When particles are generated based on the uniform background method, you can specify the thickness of the generated particles to be either variable or uniform.
| Input File Usage: | Use one of the following options to define the thickness of the generated particles: |
*SECTION CONTROLS, PARTICLE THICKNESS=VARIABLE (default) *SECTION CONTROLS, PARTICLE THICKNESS=UNIFORM |
| Abaqus/CAE Usage: | Generating particles based on a uniform background grid is not supported in Abaqus/CAE. |
Specifying a time-based criterion
The time-based criterion is primarily intended as a modeling tool to allow all particles to convert from the defined finite element mesh at the same time.
| Input File Usage: | *SECTION CONTROLS, ELEMENT CONVERSION=YES, |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Specifying a strain-based criterion
The strain-based criterion is primarily intended for cases in which you want to use a progressive conversion approach. You specify the maximum principle strain (absolute value) when continuum elements are to convert to SPH particles.
| Input File Usage: | *SECTION CONTROLS, ELEMENT CONVERSION=YES, |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Specifying a stress-based criterion
Similar to the strain-based criterion, the stress-based criterion is primarily intended for cases in which you want to use a progressive conversion approach. You specify the maximum principle stress (absolute value) when continuum elements are to convert to SPH particles.
| Input File Usage: | *SECTION CONTROLS, ELEMENT CONVERSION=YES, |
| Abaqus/CAE Usage: | Mesh module: Mesh |
Specifying a user subroutine–based criterion
The user subroutine–based criterion allows you to implement a user-defined conversion criterion. You can control element conversion during the course of an Abaqus/Explicit analysis through any of the user subroutines that can actively modify state variables associated with a material point, such as VUSDFLD and VUMAT.
| Input File Usage: | Use the following option to trigger a user subroutine–based conversion criterion: |
*SECTION CONTROLS, ELEMENT CONVERSION=YES, |
| Abaqus/CAE Usage: | Specifying a user subroutine–based criterion for element conversion is not supported in Abaqus/CAE. |