ProductsAbaqus/ExplicitAbaqus/CAE Contact surface weightingBoth the pure master-slave and the balanced master-slave contact algorithms are available in Abaqus/Explicit. By default, Abaqus/Explicit will decide which algorithm to use for any given contact pair based on the nature of the two surfaces forming the contact pair and whether kinematic or penalty enforcement of contact constraints is used. You can override the defaults in some cases. Default choices for the contact pair weightingAbaqus/Explicit uses the pure master-slave, kinematic contact algorithm, by default, in the following situations (the first surface in each situation listed is designated the master surface):
By default, Abaqus/Explicit uses the balanced master-slave, kinematic contact algorithm in the following situations:
If the penalty contact algorithm is specified, Abaqus/Explicit uses pure master-slave weighting, by default, in the following situations (the first surface in each situation listed is designated the master surface):
If the penalty contact algorithm is specified, Abaqus/Explicit chooses balanced master-slave weighting, by default, in the following situations:
Balanced master-slave weighting means that the corrections produced by both sets of contact calculations are weighted equally. Modifying the default choices for the contact pair weightingWhen the kinematic contact method is chosen, you can override the default contact pair weighting only when two separate deformable element-based surfaces are contacting each other, which corresponds to the last situation in each list for kinematic contact given in the previous section. The following aspects should be considered when deciding whether or not to override the default choice. First, the balanced master-slave contact algorithm requires more computational time, but it is typically more accurate. Second, when the densities differ by orders of magnitude, the less dense body should be a pure slave surface. Contact-induced noise can occur if a surface on a much denser body is at all weighted as a slave surface. Finally, to avoid significant penetration for hard contact, the surface with the finer mesh should not be the master surface in the pure master-slave contact pair. When the penalty contact method is chosen, you can choose to specify a pure master-slave weighting to reduce computational time. When two originally flat surfaces contact one another, a more uniform penetration distance distribution (and consequently pressure distribution) may result with pure master-slave weighting as compared to balanced master-slave weighting. This can be particularly evident if the mesh densities of the contacting surfaces differ significantly—with balanced weighting the contact penetrations will be smaller near the nodes of the coarsely meshed surface. However, balanced master-slave weighting provides better enforcement of contact constraints in most cases. You define a weighting factor, f, to specify the master-slave weighting. Set f=1.0 to designate the first surface in the contact pair as the master surface and the second surface as the slave surface. Set f=0.0 to designate the first surface in the contact pair as the slave surface and the second surface as the master surface. Specifying any value of f between 0 and 1.0 invokes the balanced master-slave contact algorithm. When f=0.5, which is the default for balanced master-slave contact pairs, Abaqus/Explicit weights each set of corrections equally. In contrast, Abaqus/Standard uses a pure master-slave contact algorithm; the slave surface must always be given first, as in the f=0.0 case above. Input File Usage CONTACT PAIR, WEIGHT=f Abaqus/CAE Usage Interaction module: interaction editor: Weighting factor Specify f Sliding formulationIn Abaqus/Explicit there are three approaches to account for the relative motion of the two surfaces forming a contact pair:
The small-sliding and infinitesimal-sliding formulations cannot be used for contact pairs using the penalty contact algorithm or involving self-contact or analytical rigid surfaces. Using the finite-sliding formulationThe finite-sliding formulation allows for arbitrary separation, sliding, and rotation of the surfaces. Abaqus/Explicit uses this formulation by default. Only the finite-sliding approach is available for self-contact or contact involving analytical rigid surfaces. Input File Usage CONTACT PAIR Abaqus/CAE Usage Interaction module: interaction editor: Sliding formulation: Finite sliding ExampleThe following input defines finite-sliding contact between the surfaces ASURF and BSURF, shown in Figure 1, with ASURF acting as the slave surface: SURFACE,NAME=ASURF ESETA, SURFACE,NAME=BSURF ESETB, CONTACT PAIR,INTERACTION=PAIR1, WEIGHT=0.0 ASURF, BSURF SURFACE INTERACTION,NAME=PAIR1 Figure 1. Contacting bodies.
In the example shown in Figure 1 slave node 101 may come into contact anywhere along the master surface BSURF. While in contact, it is constrained to slide along BSURF, irrespective of the orientation and deformation of this surface. This behavior is possible because Abaqus/Explicit tracks the position of node 101 relative to the master surface BSURF as the bodies deform. Figure 2 shows the possible evolution of the contact between node 101 and its master surface BSURF. Node 101 is in contact with the element face with end nodes 201 and 202 at time . The load transfer at this time occurs between node 101 and nodes 201 and 202 only. Later on, at time , node 101 may find itself in contact with the element face with end nodes 501 and 502. Then the load transfer will occur between node 101 and nodes 501 and 502. Figure 2. Trajectory of node 101 in finite-sliding contact.
Finite sliding in a geometrically linear analysisFinite-sliding simulations usually include nonlinear geometric effects because such simulations generally involve large deformations and large rotations. However, it is also possible to use the finite-sliding formulation in a geometrically linear analysis (see Geometric nonlinearity). The load transfer paths between the surfaces and the contact direction are updated in finite-sliding, geometrically linear analysis. This capability is useful for analyzing finite sliding between two stiff bodies that do not undergo large rotations. Using the small-sliding formulationFor a large class of contact problems the general tracking of the finite-sliding formulation is unnecessary, even though geometric nonlinearity must be considered. Abaqus/Explicit provides a small-sliding contact formulation for such problems. This formulation assumes that the surfaces may undergo arbitrarily large rotations but that a slave node will interact with the same local area of the master surface throughout the analysis. Contact pairs that use the small-sliding formulation must be defined in the first step of the simulation, although they may remain active after the first step. A large-displacement formulation (the default) should be used for the step in which the small-sliding contact formulation should be used. In a small-sliding analysis every slave node interacts with its own local tangent plane on the master surface (see Figure 3). The slave node is constrained not to penetrate this local tangent plane. Each local tangent plane, which is a line in two dimensions, is defined by an anchor point, , on the master surface and an orientation vector at the anchor point (see Figure 3). Figure 3. Definition of the anchor point and local tangent plane for node
103.
Having a local tangent plane for each slave node means that for the small-sliding formulation Abaqus/Explicit does not have to monitor slave nodes for possible contact along the entire master surface. Therefore, small-sliding contact is less expensive computationally than finite-sliding contact. The cost savings are most dramatic in three-dimensional contact problems. When the balanced master-slave contact algorithm is invoked with the small-sliding formulation, anchor points and tangent planes will be computed for both surfaces. Input File Usage Use both of the following options: STEP, NLGEOM=YES … CONTACT PAIR, SMALL SLIDING For example, the following options define small-sliding contact between the two bodies shown in Figure 1: STEP, NLGEOM=YES … SURFACE, NAME=ASURF ESETA, SURFACE, NAME=BSURF ESETB, CONTACT PAIR, SMALL SLIDING, WEIGHT=0.0 ASURF, BSURF Abaqus/CAE Usage Interaction module: interaction editor: Sliding formulation: Small sliding Step module: step editor: Nlgeom: On Anchor point and tangent plane definitionThe anchor point and the tangent plane orientation are chosen before the analysis starts using the initial configuration of the model. The anchor point and the tangent plane orientation remain fixed with respect to the master surface facet for all steps in which the contact pair is active. No contact constraints are enforced for slave nodes whose nearest point lies on the free perimeter of the master surface in the original configuration and that do not project onto any master surface facet. Abaqus/Explicit chooses the anchor point as the nearest point on the master surface. The orientation of the tangent plane is calculated by default from the normals at the master surface nodes, or you can specify it directly.
Local tangent plane rotationThe local tangent plane is always orthogonal to the contact direction. The contact direction is taken as the interpolated normal of the master surface at the anchor point, , or as the direction specified with a spatially varying clearance definition (see Specifying initial clearance values precisely). Once the contact direction has been defined, the orientation of the local tangent plane with respect to the master surface facet remains fixed. Because the small-sliding formulation considers nonlinear geometric effects, Abaqus/Explicit continuously updates the orientation of the local tangent plane to account for the rotation of the master surface facet. The position of the anchor point relative to the surrounding nodes on the master surface facet does not change as the master surface deforms. Load transferIn a small-sliding analysis the slave node will transfer load to the nodes of the master surface facet containing the anchor point, with the magnitude of the load transferred to each node weighted by its proximity to the anchor point. For example, in Figure 3 node 103 transmits load to both nodes 2 and 3 on the master surface. Thus, if node 103 impacts the local tangent plane, a larger share of the force would be transmitted to node 3 because it is closer to the anchor point . As a slave node slides along its local tangent plane, Abaqus/Explicit does not update the distribution of load transferred by a given slave node to its associated master surface nodes; the distribution is based solely on the position of the anchor point. This is unlike the small-sliding formulation in Abaqus/Standard, which does update the load distribution to the master surface nodes as sliding occurs, so that no net moment is associated with the contact forces acting on slave and master nodes per active contact constraint, regardless of the amount of sliding. Some net moment will be associated with the contact forces after sliding has occurred with the small-sliding formulation in Abaqus/Explicit. This net moment will not be significant if the sliding is truly small compared to element dimensions, but otherwise it can result in non-physical behavior and poor accounting of energy. Figure 7 shows the potential problem that arises if small sliding is used but the relative tangential motion of the surfaces is not “small.” Figure 7. Excessive sliding in a small-sliding contact analysis.
It shows the possible evolution of contact between slave node 101 in Figure 1 and its master surface BSURF. Using the unit normal vectors and , the anchor point was found for slave node 101; for the purposes of this example, assume that it lies at the midpoint of the 201–202 face. With this location of the local tangent plane for node 101 is parallel with the 201–202 face. The load transfer always occurs at the original anchor point between nodes 201 and 202, no matter how far node 101 has slid along the local tangent plane. Therefore, if node 101 moves as shown in Figure 7, it will continue to transmit load equally to nodes 201 and 202 when, in fact, it really slid off the mesh forming the master surface BSURF. What can be considered small slidingA contact pair in a small-sliding contact simulation should not grossly violate any of the assumptions or limitations outlined above. Adhere to the following guidelines:
Master surface refinement in small-sliding problemsThe basic guidelines for pure master-slave contact given previously in this section should still be followed in a small-sliding simulation. However, in a small-sliding simulation more thought must be given to the degree of refinement for the master surface. The smoothly varying master surface normal and the local tangent planes that are formed with it are crucial to the success of a small-sliding analysis. As has been mentioned previously, there are several methods that can be used to modify ; however, they only control the initial configuration of the local tangent planes. The deformation and rotation of the master surface can reorient the local tangent planes such that they become a poor representation of the master surface. Figure 8 shows an example where distortion of the master surface results in such a situation. Figure 8. Master surface deformation in a small-sliding contact analysis can
cause problems with the local tangent planes.
This problem can be minimized to some extent by using a more refined mesh on the master surface, thus providing more element faces to control the motion of the tangent planes. Excessive mesh refinement should not be necessary since only small sliding should occur. Using the infinitesimal-sliding formulationThe difference between the infinitesimal-sliding and small-sliding formulations is that the infinitesimal-sliding formulation ignores nonlinear geometric effects. To specify the infinitesimal-sliding formulation, you choose the small-sliding contact formulation and a small-displacement formulation for the analysis step. Infinitesimal sliding assumes that both the relative motions of the surfaces and the absolute motions of the model remain small. The orientations of the local tangent planes are not updated, and the load transfer paths and the weightings assigned to each master surface node remain constant during an infinitesimal-sliding simulation. Abaqus/CAE Usage Interaction module: interaction editor: Sliding formulation: Small sliding Step module: step editor: Nlgeom: Off Local tangent directions for contactLocal tangent directions for contact provide a reference frame for select contact pair output variables in Abaqus/Explicit (see About general contact in Abaqus/Explicit). These local tangent directions are separate from local coordinate systems associated with user subroutines VFRIC and VUINTER. Abaqus/Explicit establishes and updates the orientation of the first local contact tangent direction, , at slave nodes according to the logic described below for different contact formulation types within a contact pair. The orientation of the second local tangent direction, , is found as the cross product of the contact normal direction, with .
Contact tracking algorithmsA large portion of the computational cost associated with Abaqus/Explicit contact pairs derives from the algorithms used to track the relative motion between two contacting surfaces. There are two tracking approaches for the contact pair algorithm in Abaqus/Explicit, depending on the sliding formulation that is used: finite sliding and small/infinitesimal sliding. Finite-sliding trackingAbaqus/Explicit is designed to simulate highly nonlinear events or processes. Because it is possible for a node on one surface to contact any of the facets on the opposite surface, Abaqus/Explicit must use sophisticated search algorithms for tracking the motions of the surfaces. The contact search algorithm is designed to be robust, yet computationally efficient. This algorithm assumes that the incremental relative tangential motion between surfaces does not significantly exceed the dimensions of the master surface facets, but there is no limit to the overall relative motion between surfaces. It is rare for the incremental motion to exceed the facet size because of the small time increment used in explicit dynamic analyses. In cases involving relative surface velocities that exceed material wave speeds, it may be necessary to reduce the time increment. The contact search algorithm uses a global search at the beginning of each step, and a hierarchical global/local search algorithm is used for the other increments. The default contact search algorithm can handle the majority of typical contact situations. However, there are some situations that require special attention. We will consider a pure master-slave contact pair for discussion purposes. For a balanced master-slave contact pair, the contact search computations are performed twice for each contact pair. Global contact searchesA global search determines the globally nearest master surface facet for each slave node in a given contact pair. A bucket sorting algorithm is used to minimize the computational expense of these searches. A two-dimensional example, without consideration of “buckets,” is shown in Figure 9. Figure 9. Global search in two dimensions.
The global search computes the distance from node 50 to all of the master surface facets in the same bucket as node 50. It determines that the nearest facet on the master surface to node 50 is the facet of element 10. Node 100 is the node on this facet that is nearest to node 50, and it is designated the tracked master surface node. This search is conducted for each slave node, comparing each node against all of the facets on the master surface that are in the same bucket. By default, Abaqus/Explicit performs a global search every one hundred increments for two-surface contact pairs. The frequency of the global search can be manually adjusted, as discussed in Contact controls for contact pairs in Abaqus/Explicit. Despite the bucket sorting algorithm, global searches are computationally expensive: performing a global contact search in every increment will more than double the run time of many Abaqus/Explicit contact analyses. Local contact searchesAbaqus/Explicit uses a local contact search to track the motion of the surfaces during most increments of an analysis. In this approach a given slave node searches only the facets that are attached to the previously tracked master surface node. Abaqus/Explicit determines which adjacent facet is the nearest to the slave node. It then determines which node on that facet is the closest master surface node to the slave node and updates the tracked master surface node. If the closest master surface node is not the same as the previously tracked master surface node, Abaqus/Explicit performs another iteration of the local search. In the example shown in Figure 10, node 50 moves as shown during an increment. In the first iteration of the search Abaqus/Explicit finds that the master surface facet on element 10 is still the closest facet of those attached to node 100 but that node 101 is now the tracked master surface node. Because the previously tracked node was node 100, Abaqus/Explicit performs another iteration. In this second iteration a new element, element 11, is found to be the closest facet and the closest master surface node is 102. Another iteration is performed because the identity of the tracked master surface node changed. In the third iteration the identity of the tracked node does not change, so Abaqus/Explicit designates node 102 as the tracked master surface node for slave node 50. A local search is substantially less expensive computationally than a global search. A slightly more expensive local search algorithm can be employed in situations where contact is not being properly enforced; this alternate algorithm is discussed in Contact controls for contact pairs in Abaqus/Explicit. Figure 10. Local search in two dimensions.
Tracking approach for self-contact pairsAbaqus/Explicit uses similar contact searching methods for simulations with self-contact as for two-surface contact; however, more frequent global searches are often necessary for self-contact problems. By default, contact pairs with self-contact use a global contact search every four increments, compared to every 100 increments for two-surface contact pairs; the frequency of the global searches can be manually adjusted (see Contact controls for contact pairs in Abaqus/Explicit). If several facets that are unconnected to each other are found to be near a slave node during global tracking, global tracking automatically will be performed more frequently than the specified number of increments. Despite this precaution, the self-contact algorithm will be less robust if you specify a search frequency that is significantly lower than the default. Small-sliding (or infinitesimal-sliding) tracking approachWhen the small-sliding or infinitesimal-sliding contact approach is invoked (see Sliding formulation), Abaqus/Explicit performs a single global search at the beginning of the first step to determine the globally nearest master surface facet for each slave node in the given contact pair. Once the nearest facet has been determined, the nearest point on that facet defines the anchor point. Contact constraints will not be applied to slave nodes that do not project onto any master surface facet. No further tracking is performed during the step or for subsequent steps in which the contact pair remains active. This makes the small-sliding/infinitesimal-sliding contact approach less expensive computationally than the finite-sliding contact approach. The cost savings are most significant for three-dimensional contact problems. |