ProductsAbaqus/StandardAbaqus/CAE Defining contact pairsTo define a contact pair, you must indicate which pairs of surfaces may interact with one another or which surfaces may interact with themselves. Contact surfaces should extend far enough to include all regions that may come into contact during an analysis; however, including additional surface nodes and faces that never experience contact may result in significant extra computational cost (for example, extending a slave surface such that it includes many nodes that remain separated from the master surface throughout an analysis can significantly increase memory usage unless penalty contact enforcement is used). Every contact pair is assigned a contact formulation (either explicitly or by default) and must refer to an interaction property. Discussion of the various available contact formulations (based on whether the tracking approach assumes finite- or small-sliding—and whether the contact discretization is based on a node-to-surface or surface-to-surface approach) is provided in Contact formulations in Abaqus/Standard. Interaction property definitions are discussed in Assigning contact properties for contact pairs in Abaqus/Standard. Defining contact between two separate surfacesWhen a contact pair contains two surfaces, the two surfaces are not allowed to include any of the same nodes and you must choose which surface will be the slave and which will be the master. The selection of master and slave surfaces is discussed in detail in Choosing the master and slave roles in a two-surface contact pair. For simple contact pairs consisting of two deformable surfaces, the following basic guidelines can be used:
The finite-sliding, node-to-surface formulation is used by default (except in Abaqus/CAE, where the default is the finite-sliding, surface-to-surface formulation). Defining contact pairs using the finite-sliding, node-to-surface formulationA finite-sliding, node-to-surface formulation is available. Input File Usage CONTACT PAIR, INTERACTION=interaction_property_name slave_surface_name, master_surface_name You can also specify the contact discretization directly: CONTACT PAIR, INTERACTION=interaction_property_name, TYPE=NODE TO SURFACE slave_surface_name, master_surface_name Abaqus/CAE Usage Interaction module: Create Interaction: Surface-to-surface contact (Standard): select the master surface, click Surface or Node Region, select the slave surface, Interaction editor, Sliding formulation: Finite sliding, Discretization method: Node to surface, Contact interaction property: interaction_property_name Defining contact pairs using the finite-sliding, surface-to-surface formulationA node-based slave surface precludes the use of surface-to-surface discretization. Some contact capabilities are not available with the finite-sliding, surface-to-surface formulation, including crack propagation (see Crack propagation analysis). Input File Usage Use the following option to define contact constraints using the finite-sliding, surface-to-surface formulation: CONTACT PAIR, INTERACTION=interaction_property_name, TYPE=SURFACE TO SURFACE slave_surface_name, master_surface_name Abaqus/CAE Usage Interaction module: Create Interaction: Surface-to-surface contact (Standard): select the master surface, click Surface, select the slave surface, Interaction editor, Sliding formulation: Finite sliding, Discretization method: Surface to surface, Contact interaction property: interaction_property_name Defining contact pairs using the small-sliding, node-to-surface formulationThe small-sliding tracking approach uses node-to-surface discretization by default. For an explanation of when the small-sliding tracking approach is appropriate in an analysis, see Using the small-sliding tracking approach. Input File Usage CONTACT PAIR, INTERACTION=interaction_property_name, SMALL SLIDING slave_surface_name, master_surface_name You can also specify the contact discretization directly: CONTACT PAIR, INTERACTION=interaction_property_name, SMALL SLIDING, TYPE=NODE TO SURFACE slave_surface_name, master_surface_name Abaqus/CAE Usage Interaction module: Create Interaction: Surface-to-surface contact (Standard): select the master surface, click Surface or Node Region, select the slave surface, Interaction editor, Sliding formulation: Small sliding, Discretization method: Node to surface, Contact interaction property: interaction_property_name Defining contact pairs using the small-sliding, surface-to-surface formulationA node-based slave surface precludes the use of surface-to-surface discretization. Input File Usage CONTACT PAIR, INTERACTION=interaction_property_name, SMALL SLIDING, TYPE=SURFACE TO SURFACE slave_surface_name, master_surface_name Abaqus/CAE Usage Interaction module: Create Interaction: Surface-to-surface contact (Standard): select the master surface, click Surface, select the slave surface, Interaction editor, Sliding formulation: Small sliding, Discretization method: Surface to surface, Contact interaction property: interaction_property_name Using symmetric master-slave contact pairs to improve contact modelingFor node-to-surface contact it is possible for master surface nodes to penetrate the slave surface without resistance with the strict master-slave algorithm used by Abaqus/Standard. This penetration tends to occur if the master surface is more refined than the slave surface or a large contact pressure develops between soft bodies. Refining the slave surface mesh often minimizes the penetration of the master surface nodes. If the refinement technique does not work or is not practical, a symmetric master-slave method can be used if both surfaces are element-based surfaces with deformable or deformable-made-rigid parent elements. To use this method, define two contact pairs using the same two surfaces, but switch the roles of master and slave surface for the two contact pairs. This method causes Abaqus/Standard to treat each surface as a master surface and, thus, involves additional computational expense because contact searches must be conducted twice for the same contact pair. The increased accuracy provided by this method must be compared to the additional computational cost. All of the contact formulations are available for symmetric master-slave contact pairs, and can be applied using the same options discussed above. Input File Usage CONTACT PAIR, INTERACTION=interaction_property_name surface_1, surface_2 surface_2, surface_1 Abaqus/CAE Usage Interaction module: Create Interaction: Surface-to-surface contact (Standard): select the master surface, click Surface, select the slave surface Copy this interaction to a new interaction, and edit the new interaction. In the interaction editor, click Switch to reverse the master and slave surfaces. Limitations of symmetric master-slave contact pairsUsing symmetric master-slave contact pairs can lead to overconstraint problems when very stiff or hard contact conditions are enforced. See Contact constraint enforcement methods in Abaqus/Standard for a discussion of overconstraints and alternate constraint enforcement methods. For softened contact conditions, use of symmetric master-slave contact pairs will cause deviations from the specified pressure-versus-overclosure behavior, because both contact pairs contribute to the overall interface stress without accounting for one another. For example, symmetric master-slave contact pairs effectively double the overall contact stiffness if a linear pressure-overclosure relationship is specified. Likewise, use of symmetric master-slave contact pairs will cause deviations from the friction model if an optional shear stress limit is specified (see Using the optional shear stress limit), because the contact stresses observed by each contact pair will be approximately one-half of the total interface stress. Similarly, it can be difficult to interpret the results at the interface for symmetric master-slave contact pairs. In this case both surfaces at the interface act as slave surfaces, so each has contact constraint values associated with it. The constraint values that represent contact pressures are not independent of each other. Therefore, the constraint values reported in the data (.dat) and results (.fil) files represent only a part of the total interface pressure and have to be summed to obtain the total. In the output database, mechanical contact variables are reported at the nodes on both the master and slave surfaces per contact pair and not just the slave surface where constraints are formed. Consequently, two result sets are available per surface of a symmetric master-slave contact pair; once when a surface acts as a slave and once as a master. For nodal contact pressures the Visualization module of Abaqus/CAE only reports the maximum of the two pressure values associated with a node when the surface containing the node acts either as a master or as a slave surface. Even in this case, the contact pressures do not represent the true interface pressure. Apart from contact pressures, some contact output may be confusing with symmetric master-slave contact pairs. For example, Abaqus/Standard may report a positive opening distance on one side of a contact interface but zero opening distance (i.e., touching) on the opposite side of the interface. Typically this is caused by the shape or relative mesh refinement of the two surfaces. Defining self-contactDefine contact between a single surface and itself by specifying only a single surface or by specifying the same surface twice. The small-sliding tracking approach cannot be used with self-contact. Defining self-contact using node-to-surface discretizationAbaqus/Standard uses node-to-surface contact discretization by default for self-contact. Input File Usage Use either of the following options: CONTACT PAIR, INTERACTION=interaction_property_name surface_1, CONTACT PAIR, INTERACTION=interaction_property_name surface_1, surface_1 Abaqus/CAE Usage Interaction module: Create Interaction: Self-contact (Standard): select the surface Interaction editor, Discretization method: Node to surface, Contact interaction property: interaction_property_name or Interaction module: Create Interaction: Surface-to-surface contact (Standard): select the surface, click Surface, select the surface again Interaction editor, Sliding formulation: Finite sliding, Discretization method: Node to surface, Contact interaction property: interaction_property_name Defining self-contact using surface-to-surface discretizationSurface-to-surface discretization often leads to more accurate modeling of self-contact simulations. However, because the self-contact surface is acting as both a master and a slave, surface-to-surface discretization can sometimes significantly increase the solution cost. Input File Usage Use either of the following options: CONTACT PAIR, INTERACTION=interaction_property_name, TYPE=SURFACE TO SURFACE surface_1, CONTACT PAIR, INTERACTION=interaction_property_name, TYPE=SURFACE TO SURFACE surface_1, surface_1 Abaqus/CAE Usage Interaction module: Create Interaction: Self-contact (Standard): select the surface Interaction editor, Discretization method: Surface to surface, Contact interaction property: interaction_property_name or Interaction module: Create Interaction: Surface-to-surface contact (Standard): select the surface, click Surface, select the surface again Interaction editor, Sliding formulation: Finite sliding, Discretization method: Surface to surface, Contact interaction property: interaction_property_name Limitations of self-contactSelf-contact is valid only for mechanical surface interactions and is limited to finite sliding with element-based surfaces. A node of a self-contact surface can be both a slave node and a member of the master surface for two-dimensional self-contact using the surface-to-surface formulation and for all three-dimensional self-contact. In these cases the contact behavior is similar to symmetric master-slave contact pairs, and the issues discussed in Using symmetric master-slave contact pairs to improve contact modeling apply. Abaqus/Standard automatically applies some numerical softening to contact conditions in these cases, as discussed in Contact constraint enforcement methods in Abaqus/Standard. Direct enforcement of hard contact conditions is the default constraint enforcement method for two-dimensional self-contact using the node-to-surface formulation. In this case, each node adjacent to a vertex where a two-dimensional surface folds onto itself is automatically assigned a slave or master role during the analysis. Since contact constraints directly resist penetrations at nodes that act as slave nodes, there is some possibility of unresolved penetrations at nodes that only act as master nodes for two-dimensional self-contact using the node-to-surface formulation. Selecting surfaces used in contact pairsMethods for creating surfaces are discussed in Element-based surface definition, Node-based surface definition, and Analytical rigid surface definition; those sections discuss general restrictions for the various surface types. Considerations related to surface characteristics for various contact formulations are discussed in Contact formulations in Abaqus/Standard. Additional considerations for surfaces used in contact definitions are discussed below. Orientation considerations for shell-like surfacesAbaqus/Standard requires master contact surfaces to be single-sided for node-to-surface contact and for some surface-to-surface contact formulations (see Fundamental choices affecting the contact formulation for details). This requires that you consider the proper orientation for master surfaces defined on elements, such as shells and membranes, that have positive and negative directions. For node-to-surface contact the orientation of slave surface normals is irrelevant, but for surface-to-surface contact the orientation of single-sided slave surfaces is taken into consideration. Double-sided element-based surfaces are allowed for the default surface-to-surface contact formulations, although they are not always appropriate for cases with deep initial penetrations. If the master and slave surfaces are both double-sided, the positive or negative orientation of the contact normal direction will be chosen such as to minimize (or avoid) penetrations for each contact constraint. If either or both of the surfaces are single-sided, the positive or negative orientation of the contact normal direction will be determined from the single-sided surface normals rather than the relative positions of the surfaces. When the orientation of a contact surface is relevant to the contact formulation, you must consider the following aspects for surfaces on structural (beam and shell), membrane, truss, or rigid elements:
The following output from a data check analysis (see Abaqus/Standard and Abaqus/Explicit execution) can be useful in identifying incorrectly oriented master surfaces:
Surface connectivity restrictionsCertain connectivity restrictions apply to contact surfaces depending on the type of contact formulation. Surface connectivity restrictions for the various contact formulations are summarized in Table 1. As indicated in this table, the connectivity restrictions are sometimes different for master and slave surfaces. Self-contact surfaces act as both master and slave surfaces; therefore, if a restriction applies to either a master or slave surface, it also applies to self-contact. The potential connectivity restrictions referred to in Table 1 are described below:
Analytical rigid surfacesAnalytical rigid surfaces are often effective for efficiently modeling curved, rigid geometries, as discussed in Analytical rigid surface definition. For rare cases in which a very large number (thousands) of segments would be necessary to define an analytical rigid surface, better performance can be achieved with an element-based rigid surface (see Element-based surface definition). Three-dimensional beam and truss surfacesAbaqus/Standard cannot use three-dimensional beams or trusses to form a master surface because the elements do not have enough information to create unique surface normals. However, these elements can be used to define a slave surface. Two-dimensional beams and trusses can be used to form both master and slave surfaces. Edge-based surfacesEdge-based surfaces (Element-based surface definition) on three-dimensional shell elements cannot be used in a contact analysis in Abaqus/Standard. Limitations of node-based surfacesUse node-based surfaces with caution when the contact property definition includes user-defined softened contact properties or thermal or electrical interactions because the contact constitutive behavior (which relies on accurate calculation of contact pressure, heat flux, or electric current) will not be enforced correctly unless the precise surface area is associated with each node. For details, see Contact pressure-overclosure relationships, Thermal contact properties, or Electrical contact properties. Removing and reactivating contact pairsYou can temporarily remove contact pairs from a simulation, which may result in significant computational savings by eliminating unnecessary contact searches and updates of surface orientations during the simulation. Removal and reactivation of contact pairs is commonly used in complicated forming processes where multiple tools need to interact with the workpiece at different stages in the analysis. You cannot remove tied contact pairs from a simulation (see Defining tied contact in Abaqus/Standard). Removing contact pairsRemoval of contact pairs is a useful technique for uncoupling components of an assembly until they should be brought together (such as tooling in manufacturing process simulations). Significant computational expense may be saved by removing a contact pair and introducing it at the proper time, thus eliminating the need to monitor the contact conditions except when they are relevant. Input File Usage MODEL CHANGE, TYPE=CONTACT PAIR, REMOVE slave_surface, master_surface Repeat the data line as needed. Abaqus/CAE Usage Use one of the following options: Interaction module: Create Interaction: surface-to-surface contact or self-contact interaction editor: toggle off Active in this step Interaction module: interaction manager: select interaction, Removal of contact forces associated with closed contact pairsIf the surfaces are in contact when a contact pair is removed, Abaqus/Standard stores the corresponding contact forces (or heat fluxes if thermal interactions are present, or electrical currents if it is a coupled-thermal electrical analysis) for every node on each surface. Abaqus/Standard automatically ramps these forces (or heat fluxes or electrical currents) linearly down to zero magnitude during the removal step. Abaqus/Standard always removes the contact constraints for mechanical surface interactions instantaneously. Care must be taken in removing contact pairs in transient procedures. In transient heat transfer, fully coupled temperature-displacement, or fully coupled thermal-electrical-structural analysis if the fluxes are high and the step is long, this ramping down may have the effect of cooling down or heating up the rest of the body. In dynamic analysis if the forces are high and the step is long, kinetic energy can be imparted to the remaining portion of the model. This problem can be avoided by removing the contact pairs in a very short transient step prior to the rest of the analysis. This step can be done in a single increment. Using an allowable contact interference to deactivate contact pairsA contact pair with mechanical contact interactions can be deactivated during an analysis by assigning a very large allowable contact interference to the contact pairs (see Modeling contact interference fits in Abaqus/Standard). This method has the disadvantage of not reducing the computational cost of the analysis because the contact algorithm will still calculate the contact conditions for the contact pair in each increment. Reactivating contact pairsAll contact pairs that will be used in a simulation must be created at the start of the analysis; they cannot be created once the simulation has begun. However, contact pairs can be created, removed at the start of the analysis in the first step, and then reactivated at a later point during the simulation. In Abaqus/CAE you can create contact pairs in any step. If a contact pair is created in a step other than the initial step, Abaqus/CAE automatically deactivates the contact pair in the initial step and reactivates it in the step in which you created it. Input File Usage MODEL CHANGE, TYPE=CONTACT PAIR, ADD slave_surface, master_surface Repeat the data line as needed. Abaqus/CAE Usage Interaction module: Create Interaction: surface-to-surface contact or self-contact interaction editor: toggle on Active in this step Reactivating overclosed contact pairsWhen a contact pair is reactivated, the contact constraint becomes active immediately. In mechanical simulations it is possible for the surfaces of a contact pair to move such that they become overclosed while the contact pair is inactive. If this overclosure is too severe when the contact pair is reactivated, Abaqus/Standard may encounter convergence problems as it tries to enforce the suddenly activated contact constraint. To avoid such problems, you can specify a permissible interference value, v, for the contact pair that is larger than the overclosure for the contact pair. Abaqus/Standard will ramp v down to zero during the step. For details on specifying allowable interferences, see Modeling contact interference fits in Abaqus/Standard. OutputOutput variables associated with the interaction of contact pairs fall into two categories: nodal variables (sometimes called constraint variables) and whole surface variables. In addition, Abaqus outputs an array of diagnostic information associated with contact interactions, as discussed in Contact diagnostics in an Abaqus/Standard analysis. For more detailed discussions of variables associated with thermal, electrical, and pore fluid analyses, see the sections on the related contact properties in Contact Property Models. Nodal contact variablesNodal contact variables can be contoured on contact surfaces in the Visualization module of Abaqus/CAE. Nodal contact variables include contact pressure and force, frictional shear stress and force, relative tangential motion (slip) of the surfaces during contact, clearance between surfaces, heat or fluid flux per unit area, fluid pressure, and electrical current per unit area. Many of the nodal contact variables written to the output database (.odb) file are often available for all contact nodes, regardless of whether they act as slave or master nodes. Other nodal contact variables are available only at nodes acting as slave nodes. Most contact output to the data (.dat) file, results (.fil) file, and the utility subroutine GETVRMAVGATNODE is associated with individual constraints. For contact output to the output database (.odb) file, some filtering is applied to reduce contact output noise. The contact pressure distribution is of key interest in many Abaqus analyses. You can view the contact pressure on all contact surfaces except for analytical rigid surfaces and discrete rigid surfaces based on rigid-type elements (the latter restriction does not apply to general contact). You can view a contour plot of the contact pressure error indicator next to a contour plot of the contact pressure to gain perspective on local accuracy of the contact pressure solution in regions where the contact pressure solution is of interest (see Selection of error indicators influencing adaptive remeshing for further discussion of error indicator output). In some cases you may observe the contact pressure extending beyond the actual contact zone due to the following factors:
Due to these factors, trying to infer the contact force distribution from the contact stress distribution can be somewhat misleading. Instead, you can request nodal contact force output, which accurately represents the contact force distribution present in the analysis. Whole surface variablesWhole surface variables are attributes of an entire slave surface. Available as history output, these variables record the total force and moment due to contact pressure and frictional stress, the center of pressure and frictional stress (defined as the point closest to the centroid of the surface that lies on the line of action of the resultant force for which the resultant moment is minimal), and the total contact area (defined as the sum of all the facets where there is contact force). The last letter of each variable name (except the variable CAREA) denotes which contact force distribution on the surface is used to calculate the resultant:
For example, CFN is the total force due to contact pressure, CFS is the total force due to frictional stress, and CFT is the total force due to both contact pressure and frictional stress. Each total moment output variable will not necessarily equal the cross product of the respective center of force vector and resultant force vector. Forces acting on two different nodes of a surface may have components acting in opposite directions, such that these nodal force components generate a net moment but not a net force; therefore, the total moment may not arise entirely from the resultant force. The center of force output variables tend to be most meaningful when the surface nodal forces act in approximately the same direction. Requesting outputCertain contact variables must be requested as a group. For example, to output the clearance between surfaces (COPEN), you must request the variable CDISP (contact displacements). CDISP outputs both COPEN and CSLIP (tangential motion of the surfaces during contact). A complete listing of available contact pair variables and identifiers is given in Abaqus/Standard output variable identifiers. Output requests can be limited to individual contact pairs or portions of a slave surface. You can:
Instructions on forming these output requests are available in the following sections:
Differences for small-sliding and finite-sliding contactFor small-sliding contact problems the contact area is calculated in the input file preprocessor from the undeformed shape of the model; thus, it does not change throughout the analysis, and contact pressures for small-sliding contact are calculated according to this invariant contact area. This behavior is different from that in finite-sliding contact problems, where the contact area and contact pressures are calculated according to the deformed shape of the model. Output of tangential resultsAbaqus reports the values of tangential variables (frictional shear stress, viscous shear stress, and relative tangential motion) with respect to the local tangent directions defined on the surfaces. The local tangent directions CTANDIR1 and CTANDIR2 can be output by requesting the generic output variable CTANDIR. The definition of local tangent directions is explained in Local tangent directions on a surface. These directions do not always correspond to the global coordinate system, and they rotate with the contact pair in a geometrically nonlinear analysis. Abaqus/Standard calculates tangential results at each constraint point by taking the scalar product of the variable's vector and a local tangent direction, or , associated with the constraint point. The number at the end of a variable's name indicates whether the variable corresponds to the first or second local tangent direction. For example, CSHEAR1 is the frictional shear stress component in the first local tangent direction, while CSHEAR2 is the frictional shear stress component in the second local tangent direction. Definition of accumulated incremental relative motion (slip)Abaqus/Standard defines the incremental relative motion (also known as slip) as the scalar product of the incremental relative nodal displacement vector and a local tangent direction. The incremental relative nodal displacement vector measures the motion of a slave node relative to the motion of the master surface. The incremental slip is accumulated only when the slave node is contacting the master surface. The sums of all such incremental slips during the analysis are reported as CSLIP1 and CSLIP2. Details about the calculation of this quantity can be found in Small-sliding interaction between bodies, Finite-sliding interaction between deformable bodies, and Finite-sliding interaction between a deformable and a rigid body. Extending the range for which contact opening output is provided for gapsTo reduce computational costs, detailed computations to monitor potential points of interaction are avoided by default where surfaces are separated by a distance greater than the minimum gap distance at which contact forces (or thermal fluxes, etc.) may be transmitted. Therefore, contact opening (COPEN) output is typically not provided for finite-sliding contact where surfaces are opened by more than a small amount compared to surface facet dimensions. You can extend the range in which Abaqus/Standard provides contact opening output; COPEN will be provided up to gap distances equal to a specified “tracking thickness.” Using this control may increase computational cost due to extra contact tracking computations, especially if you specify a large tracking thickness value. Input File Usage SURFACE INTERACTION, TRACKING THICKNESS=value Abaqus/CAE Usage You cannot adjust the default tracking thickness in Abaqus/CAE. Output for axisymmetric modelsIn an axisymmetric analysis the total forces and moments transmitted between the contacting bodies as a result of contact pressure and frictional stress are computed in the same manner as in a two-dimensional analysis. Therefore, the component of the total forces along the r-axis is nonzero, and the components of the total moments include contributions from the total forces along the r-axis. Obtaining the “maximum torque” that can be transmitted about the z-axis in an axisymmetric analysisWhen modeling surface-based contact with axisymmetric elements (element types CAX and CGAX), Abaqus/Standard can calculate the maximum torque (output variable CTRQ) that can be transmitted about the z-axis. This capability is often of interest when modeling threaded connectors (see Axisymmetric analysis of a threaded connection). The maximum torque, T, is defined as where p is the pressure transmitted across the interface, r is the radius to a point on the interface, and s is the current distance along the interface in the r–z plane. This definition of “torque” effectively assumes a friction coefficient of unity. |