This section presents an overview of the contact analysis capabilities in Abaqus.
Abaqus provides more than one approach for defining contact. Abaqus/Standard includes the following approaches for defining contact:
general contact;
contact pairs; and
contact elements.
general contact; and
contact pairs.
The remainder of this section is organized as follows:
first, discuss common aspects of the surface-based contact-definition approaches (i.e., contact pairs and general contact);
next, provide an overview of the contact definition approaches in Abaqus/Standard and the contact definition approaches in Abaqus/Explicit;
finally, discuss compatibility between the contact algorithms in Abaqus/Standard and Abaqus/Explicit.
A contact simulation using contact pairs or general contact is defined by specifying:
surface definitions for the bodies that could potentially be in contact;
the surfaces that interact with one another (the contact interactions);
any nondefault surface properties to be considered in the contact interactions;
the mechanical and thermal contact property models, such as the pressure-overclosure relationship, the friction coefficient, or the contact conduction coefficient;
any nondefault aspects of the contact formulation; and
any algorithmic contact controls for the analysis.
Surfaces can be defined at the beginning of a simulation or upon restart as part of the model definition (see “Surfaces: overview,” Section 2.3.1). Abaqus has four classifications of contact surfaces:
element-based deformable and rigid surfaces (“Element-based surface definition,” Section 2.3.2);
node-based deformable and rigid surfaces (“Node-based surface definition,” Section 2.3.3);
analytical rigid surfaces (“Analytical rigid surface definition,” Section 2.3.4); and
Eulerian material surfaces for Abaqus/Explicit (“Eulerian surface definition,” Section 2.3.5).
When the general contact algorithm is used, Abaqus also provides a default all-inclusive, automatically defined surface that includes all element-based surface facets (in Abaqus/Standard and in Abaqus/Explicit), all analytical rigid surfaces (in Abaqus/Explicit only), and all Eulerian materials (in Abaqus/Explicit only) in the model.
Contact interactions for contact pairs and general contact are defined by specifying surface pairings and self-contact surfaces. General contact interactions typically are defined by specifying self-contact for the default surface, which allows an easy, yet powerful, definition of contact. (Self-contact for a surface that spans multiple bodies implies self-contact for each body as well as contact between the bodies.)
At least one surface in an interaction must be a non-node-based surface, and at least one surface in an interaction must be a non-analytical rigid surface. Additional restrictions and guidelines for contact surfaces are discussed for each contact definition approach. The definition of contact pairs is discussed in detail in “Defining contact pairs in Abaqus/Standard,” Section 35.3.1, and “Defining contact pairs in Abaqus/Explicit,” Section 35.5.1. The definition of general contact interactions is discussed in detail in “Defining general contact interactions in Abaqus/Standard,” Section 35.2.1, and “Defining general contact interactions in Abaqus/Explicit,” Section 35.4.1.
Nondefault surface properties (such as thickness and, in some cases, offset) can be defined for particular surfaces in a contact model. In addition, you can control which edges of a surface will be included in the general contact domain in Abaqus/Explicit. Surface properties for contact pairs are discussed in “Assigning surface properties for contact pairs in Abaqus/Standard,” Section 35.3.2, and “Assigning surface properties for contact pairs in Abaqus/Explicit,” Section 35.5.2. Surface properties for general contact are discussed in “Surface properties for general contact in Abaqus/Standard,” Section 35.2.2, and “Assigning surface properties for general contact in Abaqus/Explicit,” Section 35.4.2.
Contact interactions in a model can refer to a contact property definition, in much the same way that elements refer to an element property definition. By default, the surfaces interact (have constraints) only in the normal direction to resist penetration. The other mechanical contact interaction models available depend on the contact algorithm and whether Abaqus/Standard or Abaqus/Explicit is used (see “Mechanical contact properties: overview,” Section 36.1.1). Some of the available models are:
softened contact (“Contact pressure-overclosure relationships,” Section 36.1.2, and “Frictional behavior,” Section 36.1.5);
contact damping (“Contact damping,” Section 36.1.3, and “Frictional behavior,” Section 36.1.5);
friction (“Frictional behavior,” Section 36.1.5);
a user-defined constitutive model for surface interactions (“User-defined interfacial constitutive behavior,” Section 36.1.6); and
spot welds bonding two surfaces together until the welds fail (“Breakable bonds,” Section 36.1.9).
Contact interaction models are defined as model data except for contact pairs in Abaqus/Explicit, in which case they are defined as history data. Information on assigning contact properties to contact pairs can be found in “Assigning contact properties for contact pairs in Abaqus/Standard,” Section 35.3.3, and “Assigning contact properties for contact pairs in Abaqus/Explicit,” Section 35.5.3. Information on assigning contact properties to general contact interactions can be found in “Contact properties for general contact in Abaqus/Standard,” Section 35.2.3, and “Assigning contact properties for general contact in Abaqus/Explicit,” Section 35.4.3.
The default algorithmic controls for contact analyses are usually sufficient, but you can adjust numerical controls for some special cases. For example, depending on the contact algorithm used, the numerical controls for the contact formulation, the master and slave roles for the contact surfaces, and the sliding formulation are provided. Information on contact formulations and numerical methods used by the contact algorithms is provided in “Contact formulations in Abaqus/Standard,” Section 37.1.1, and “Contact formulations for contact pairs in Abaqus/Explicit,” Section 37.2.2. The available numerical controls for the various contact algorithms are discussed in “Numerical controls for general contact in Abaqus/Standard,” Section 35.2.6; “Adjusting contact controls in Abaqus/Standard,” Section 35.3.6; “Contact controls for general contact in Abaqus/Explicit,” Section 35.4.5; and “Contact controls for contact pairs in Abaqus/Explicit,” Section 35.5.5.
Abaqus/Standard provides the following approaches for defining contact interactions: general contact, contact pairs, and contact elements. Contact pairs and general contact both use surfaces to define contact; comparisons of these approaches are provided later in this section. Contact elements are provided for certain interactions that cannot be modeled with either general contact or contact pairs; however, it is generally recommended to use general contact or contact pairs if possible.
Contact pairs and general contact combine to provide the following capabilities in Abaqus/Standard:
Contact between two deformable bodies. The structures can be either two- or three-dimensional, and they can undergo either small or finite sliding. Examples of such problems include the assembly of a cylinder head gasket and the slipping between the two components of a threaded connector.
Contact between a rigid surface and a deformable body. The structures can be either two- or three-dimensional, and they can undergo either small or finite sliding. Examples of such problems include metal forming simulations and analyses of rubber seals being compressed between two components.
Finite-sliding self-contact of a single deformable body. An example of such a problem is a complex rubber seal that folds over on itself.
Small-sliding or finite-sliding interaction between a set of points and a rigid surface. These models can be either two- or three-dimensional. An example of this type of problem is the pull-in of an underwater cable that is resting on the seabed, with the seabed modeled as a rigid surface.
Contact between a set of points and a deformable surface. These models can be either two- or three-dimensional. An example of this class of contact problem is the design of a bearing where one of the bearing surfaces is modeled with substructures.
Problems where two separate surfaces need to be “tied” together so that there is no relative motion between them. This modeling technique allows for joining dissimilar meshes.
Coupled thermal-mechanical interaction between deformable bodies with finite relative motion. The analysis of a disc brake is an example of such a problem.
Coupled thermal-electrical-structural interaction between deformable bodies with finite relative motion. An example of this type of problem is the analysis of resistance spot welding.
Coupled pore fluid-mechanical interaction between bodies. An example of this type of problem is the analysis of the interfaces between layered soil material at a waste disposal site.
For most contact problems you have a choice of whether to define contact interactions using general contact or contact pairs. In Abaqus/Standard the distinction between general contact and contact pairs lies primarily in the user interface, the default numerical settings, and the available options. The general contact and contact pair implementations share many underlying algorithms.
The contact interaction domain, contact properties, and surface attributes are specified independently for general contact, offering a more flexible way to add detail incrementally to a model. The simple interface for specifying general contact allows for a highly automated contact definition; however, it is also possible to define contact with the general contact interface to mimic traditional contact pairs. Conversely, specifying self-contact of a surface spanning multiple bodies with the contact pair user interface (if the surface-to-surface formulation is used) mimics the highly automated approach often used for general contact.
In Abaqus/Standard, traditional pairwise specifications of contact interactions will often result in more efficient or robust analyses as compared to an all-inclusive self-contact approach to defining contact. Therefore, there is often a trade-off between ease of defining contact and analysis performance. Abaqus/CAE provides a contact detection tool that greatly simplifies the process of creating traditional contact pairs for Abaqus/Standard (see “Understanding contact and constraint detection,” Section 15.6 of the Abaqus/CAE User's Manual).
Differences in default settings for general contact and contact pairs in Abaqus/Standard include the following:
Contact formulation: General contact uses the finite-sliding, surface-to-surface formulation supplemented by the finite-sliding, edge-to-surface formulation. Contact pairs use the finite-sliding, node-to-surface formulation by default except when the contact detection tool in Abaqus/CAE is used to create the contact pairs, in which case the finite-sliding, surface-to-surface formulation is used by default. See “Contact formulations in Abaqus/Standard,” Section 37.1.1, for a discussion of contact formulations.
Treatment of shell thickness and offset: General contact automatically accounts for thicknesses and offsets associated with shell-like surfaces. Contact pairs that use the finite-sliding, node-to-surface formulation do not account for shell thicknesses and offsets. See “Surface properties for general contact in Abaqus/Standard,” Section 35.2.2, and “Assigning surface properties for contact pairs in Abaqus/Standard,” Section 35.3.2, for discussions of surface properties for contact in Abaqus/Standard.
Contact constraint enforcement: General contact uses the penalty method to enforce the contact constraints by default. Contact pairs that use the finite-sliding, node-to-surface formulation use a Lagrange multiplier method to enforce contact constraints by default in most cases. See “Contact constraint enforcement methods in Abaqus/Standard,” Section 37.1.2, for a discussion of contact constraint enforcement methods.
Treatment of initial overclosures: General contact eliminates initial overclosures with strain-free adjustments by default. Contact pairs treat initial overclosures as interference fits to be resolved in the first increment of the analysis by default. See “Controlling initial contact status in Abaqus/Standard,” Section 35.2.4; “Modeling contact interference fits in Abaqus/Standard,” Section 35.3.4; and “Adjusting initial surface positions and specifying initial clearances in Abaqus/Standard contact pairs,” Section 35.3.5; for more information on contact initialization in Abaqus/Standard.
Master-slave assignments: General contact automatically assigns pure master and slave roles for most contact interactions and automatically assigns balanced master-slave roles to other contact interactions. The user must assign master and slave roles for most contact pairs. See “Numerical controls for general contact in Abaqus/Standard,” Section 35.2.6, and “Choosing the master and slave roles in a two-surface contact pair” in “Contact formulations in Abaqus/Standard,” Section 37.1.1, for discussions of master and slave roles for contacting surfaces.
The following capabilities are available only for contact pairs in Abaqus/Standard (they are not available for general contact in Abaqus/Standard):
Contact involving analytical rigid surfaces or rigid surfaces defined with user subroutine RSURFU (however, element-based rigid surfaces can be included in either general contact or contact pairs).
Contact involving node-based surfaces or surfaces on three-dimensional beam elements.
Small-sliding contact and tied contact.
The finite-sliding, node-to-surface contact formulation.
Debonding and cohesive contact behavior.
Surface interactions in analyses without displacement degrees of freedom, such as pure heat transfer.
Pressure-penetration loading.
Local definitions of some numerical contact controls.
Symmetric model generation.
A single analysis can include general contact and contact pair definitions. For example, you may choose to model contact interactions involving analytical rigid surfaces with contact pairs and other contact interactions with general contact. General contact automatically avoids processing contact interactions that are treated by contact pairs.
Surface-based contact methods associated with general contact and contact pairs cannot be used for certain classes of problems. Abaqus/Standard provides a library of contact elements for these problems. Examples of such problems are:
Contact interaction between two pipelines or tubes modeled with pipe, beam, or truss elements where one pipe lies inside the other (such as a J-tube pull in offshore piping installation) or the pipes lie next to each other (available in both two and three dimensions; see “Tube-to-tube contact elements,” Section 39.3.1).
Contact between two nodes along a fixed direction in space. An example of such a problem is the interaction of a piping system with its supports (see “Gap contact elements,” Section 39.2.1).
Simulations using axisymmetric elements with asymmetric deformations, CAXAn and SAXAn elements. See “Contact modeling if asymmetric-axisymmetric elements are present,” Section 35.3.10, for details.
Heat transfer analyses where the heat flow is one-dimensional. An example of such a problem is the heat flow in a piping system that is discontinuous. The thermal interaction in this problem is one-dimensional, so no surfaces can be defined (see “Gap contact elements,” Section 39.2.1).
The steps required for defining a contact simulation using contact elements are similar to those needed when defining a surface-based contact simulation:
create the contact elements or slide lines;
assign element section properties to the contact elements;
associate sets of contact elements with the slide lines if applicable; and
define the contact property models for the contact elements.
Abaqus/Explicit provides two algorithms for modeling contact interactions. The general (“automatic”) contact algorithm allows very simple definitions of contact with very few restrictions on the types of surfaces involved (see “Defining general contact in Abaqus/Explicit,” Section 35.4). The contact pair algorithm has more restrictions on the types of surfaces involved and often requires more careful definition of contact; however, it allows for some interaction behaviors that currently are not available with the general contact algorithm (see “Defining contact pairs in Abaqus/Explicit,” Section 35.5). The general contact and contact pairs algoirthms in Abaqus/Explicit differ by more than the user interface; in general they use completely separate implementations with many key differences in the designs of the numerical algorithms.
The two contact algorithms combine to provide the following capabilities in Abaqus/Explicit:
Contact between rigid and/or deformable bodies.
Contact of a body with itself.
Finite-sliding or small-sliding contact.
Contact with eroding bodies (due to element failure). A node-based surface must be used to model the eroding body if contact pairs are used. General contact allows element-based surfaces to be defined on eroding bodies, so contact between any number of eroding bodies can be modeled.
General constitutive models for the contact behavior, including user-defined models through user subroutines, relating constraint pressure and shear traction to penetration distance and relative tangential motion.
Thermal interaction at the surface of a body; for example, conductive heat transfer.
Contact between Eulerian material and Lagrangian bodies.
A friction coefficient defined in terms of average surface temperature and/or field variables.
Contact definitions are not entirely automatic with the general contact algorithm but are greatly simplified. The generality of this algorithm is primarily in the relaxed restrictions on the surfaces that can be used in contact. The general contact algorithm in Abaqus/Explicit allows the following (none of which are allowed with the contact pair algorithm in Abaqus/Explicit):
A surface can span unattached bodies.
More than two surface facets can share a common edge (allowing “T-intersections” in shells, for example).
A surface can include deformable and rigid regions; furthermore, the rigid regions need not be from the same rigid body.
A surface can have mixed parent element types; for example, adjacent surface facets can be on shell and solid elements.
A surface can be based on combinations of surfaces of the same type.
An element-based surface can be defined on the interior of solid bodies for use in modeling erosion due to element failure.
A surface can be defined on the exterior of an Eulerian material instance (see “Eulerian surface definition,” Section 2.3.5).
The general contact algorithm can enforce edge-to-edge contact for geometric feature edges, perimeter edges of structural elements, and edges defined by beam and truss elements, unlike the contact pair algorithm.
The general contact algorithm is the only option for enforcing contact between Eulerian materials and Lagrangian bodies (see “Interactions” in “Eulerian analysis,” Section 14.1.1).
The general contact algorithm eliminates problematic, nonphysical “bull-nose” extensions that may arise at shell surface perimeters in the contact pair algorithm.
With the general contact algorithm each slave node can see contact with multiple facets per increment; with the contact pair algorithm each slave node can see contact with only one facet per increment unless multiple surface pairings are specified. Likewise, each contact edge can see contact with multiple edges per increment when the general contact algorithm is used.
The general contact algorithm has some built-in smoothing for element-based surfaces that can be beneficial for modeling contact near corners.
The general contact algorithm, unlike the contact pair algorithm, removes contact faces and contact edges from the contact domain and, if an interior surface is defined, activates newly exposed surface faces as elements fail. Thus, element-based surfaces can be used to describe eroding solids. This allows contact between multiple eroding solids to be modeled since a node-based surface does not need to be defined on the eroding solid.
Contact state information (such as the proper contact normal orientation for double-sided surfaces) is transferred across step boundaries in the general contact algorithm even if the contact domain is modified; in the contact pair algorithm, contact state information is transferred across step boundaries only for contact pairs with no modifications.
The contact interaction domain, contact properties, and surface attributes are specified independently for the general contact algorithm, offering a more flexible way to add detail incrementally to a model.
The general contact algorithm does not place any restrictions on the domain decomposition for domain level parallelization (see “Parallel execution in Abaqus/Explicit,” Section 3.5.3).
The general contact algorithm in Abaqus/Explicit has been developed to minimize the need for algorithmic controls.
Although the general contact algorithm is more powerful and allows for simpler contact definitions, the contact pair algorithm must be used in certain cases where more specialized contact features are desired. The following features are available in Abaqus/Explicit only when the contact pair algorithm is used:
Two-dimensional surfaces
Kinematically enforced contact (see “Contact constraint enforcement methods in Abaqus/Explicit,” Section 37.2.3; the general contact algorithm uses only penalty enforcement)
Small-sliding contact (see “Contact formulations for contact pairs in Abaqus/Explicit,” Section 37.2.2)
Exponential and no separation contact pressure-overclosure models
Breakable bonds, such as spot welds (however, mesh-independent spot welds can be used with either contact algorithm; see “Mesh-independent fasteners,” Section 34.3.4)
The two contact algorithms can be used together in the same Abaqus/Explicit analysis. The general contact algorithm automatically avoids processing interactions that are treated by the contact pair algorithm.
There are fundamental differences in the mechanical contact algorithms in Abaqus/Standard and Abaqus/Explicit even though the input syntax is similar. The main differences are the following:
Contact pair and general contact definitions in Abaqus/Standard are model definition data (although contact pairs can be removed for a portion of the analysis and added back to the model in a later step of the analysis, as discussed in “Removing and reactivating contact pairs” in “Defining contact pairs in Abaqus/Standard,” Section 35.3.1). In the contact pair algorithm in Abaqus/Explicit contact constraints are history definition data (see “Defining a model in Abaqus,” Section 1.3.1); in the general contact algorithm in Abaqus/Explicit contact definitions can be either model or history data.
Abaqus/Standard typically uses a pure master-slave relationship for the contact constraints; whereas Abaqus/Explicit typically uses balanced master-slave contact by default. This difference is primarily due to overconstraint issues unique to Abaqus/Standard.
The contact formulations in Abaqus/Standard and Abaqus/Explicit differ in many respects due to different convergence, performance, and numerical requirements:
Abaqus/Standard provides surface-to-surface and edge-to-surface formulations, which Abaqus/Explicit does not;
Abaqus/Explicit provides an edge-to-edge formulation, which Abaqus/Standard does not;
Abaqus/Standard and Abaqus/Explicit both provide node-to-surface formulations, but some details associated with surface smoothing, etc. differ in the respective implementations.
The constraint enforcement methods in Abaqus/Standard and Abaqus/Explicit differ in some respects. For example, both analysis codes provide penalty constraint methods, but the default penalty stiffnesses differ (this is primarily due to the effect of the penalty stiffness on the stable time increment for Abaqus/Explicit).
The small-sliding contact capability in Abaqus/Standard transfers the load to the master nodes according to the current position of the slave node, but the small-sliding contact capability in Abaqus/Explicit always transfers the load through the anchor point due to a numerical limitation associated with the implementation.
Abaqus/Explicit can account for the thickness and midsurface offset of shells and membranes in the contact penetration calculations (although in some cases changes in the thickness upon deformation are not accounted for in the contact calculations). Abaqus/Standard cannot account for the thickness and offset of shells and membranes when using the finite-sliding, node-to-surface contact formulation (but can account for the original thickness and offset in all other contact formulations).