CONTA177
3D Line-to-Surface
Contact
CONTA177 Element Description
Use CONTA177 to represent contact and sliding between 3D target surfaces and a deformable line segment, defined by this element. The element is applicable to 3D beam-to-surface, 3D shell edge-to-surface, and 3D beam-to-beam (or edge-to-edge) structural contact analyses. It supports pair-based contact and general contact.
This element is located on the surfaces of 3D beam or pipe elements with or without midside nodes (such as BEAM188, BEAM189, PIPE288, PIPE289, and ELBOW290). It can also be located on feature edges of 3D solid elements and perimeter edges of 3D shell elements, with or without midside nodes (such as SHELL181 and SHELL281). Contact occurs when the element surface penetrates an associated target surface.
For pair-based contact, the target surface is defined by the 3D target element type, TARGE170. In the case of general contact, the target surface can be defined by CONTA174 elements (for deformable surfaces), CONTA177 elements (for 3D beams and 3D edges), or TARGE170 elements (for rigid bodies only).
Coulomb friction, shear stress friction, user-defined friction
with the USERFRIC
subroutine, and user-defined
contact interaction with the USERINTER
subroutine
are allowed. This element also allows separation of bonded contact
to simulate interface delamination.
See CONTA177 in the Mechanical APDL Theory Reference for more details about this element.
CONTA177 Input Data
The geometry and node locations are shown in Figure 177.1: CONTA177 Geometry. The element is defined by two nodes (if the underlying beam or shell element does not have a midside node) or three nodes (if the underlying beam or shell element has a midside node). The element x axis is along the I-J line of the element. Correct node ordering of the contact element is critical for proper detection of contact. The nodes must be ordered in a sequence that defines a continuous line.
Four different scenarios can be modeled by CONTA177:
Contact between one beam (or edge) and the surface of a solid or shell
Internal contact where one beam (or pipe) slides inside another hollow beam (or pipe); see Figure 177.2: Beam Sliding Inside a Hollow Beam
External contact between two beams (or edges) that lie next to each other and are roughly parallel; see Figure 177.3: Parallel Beams in Contact
External or internal contact between two beams (or edges) that cross; see Figure 177.4: Crossing Beams in Contact
KEYOPT(3) controls which of the above scenarios are allowed for the element type, and also controls the contact model used (force-based or traction-based):
Use KEYOPT(3) = 0 for the first three scenarios. The contact condition is only checked at contact nodes. The program reports contact force (contact force-based model).
Use KEYOPT(3) = 1 for the first three scenarios. The contact condition is only checked at contact nodes. The program reports contact pressure (contact traction-based model).
Use KEYOPT(3) = 2 for all scenarios. The contact condition is only checked at contact nodes for the first three scenarios, and on an intersection along the beams for the fourth scenario. The program reports contact pressure (contact traction-based model).
Use KEYOPT(3) = 3 for the fourth scenario. The contact condition is only checked on an intersection along the beams. The program reports contact pressure (contact traction-based model).
Use KEYOPT(3) = 4 for the fourth scenario. The contact condition is only checked on an intersection along the beams. The program reports contact force (contact force-based model).
The units for certain real constants (FKN, FKT, TNOP) and postprocessing items (PRES, TAUR, TAUS, SFRIC, and so on) vary by a factor of AREA, depending on whether the contact force-based model (KEYOPT(3) = 0 or 4) or the contact traction-based model (KEYOPT(3) = 1, 2, or 3) is specified. See the real constant table and output definitions table for details. For more information, see KEYOPT(3) in the Contact Technology Guide.
Pair-Based Contact versus General Contact
There are two methods to define a contact interaction: the pair-based contact definition and the general contact definition. Both contact definitions can exist in the same model. CONTA177 can be used in either type of contact definition.
The pair-based contact definition is usually more efficient and more robust than the general contact definition; it supports more options and specific contact features.
Pair-Based Contact
In a pair-based contact definition, the 3D line contact elements (CONTA177) are associated with 3D target segment elements (TARGE170) via a shared real constant set. The program looks for contact only between contact and target surfaces with the same real constant set ID (which is greater than zero). The material ID associated with the contact element is used to specify interaction properties (such as friction coefficient) defined via the MP or TB command.
If more than one target surface will make contact with the same boundary of line elements, you must define several contact elements that share the same geometry but relate to separate targets (targets which have different real constant numbers). Alternatively, you can combine several target surfaces into one (that is, multiple targets sharing the same real constant numbers). See Identifying Contact Pairs in the Contact Technology Guide for more information.
For rigid-flexible and flexible-flexible contact, one of the deformable surfaces (beam or shell edge) must be represented by a contact surface. See Designating Contact and Target Surfaces in the Contact Technology Guide for more information.
See Generating Contact Elements in the Contact Technology Guide for more information on generating elements automatically via the ESURF command.
General Contact
In a general contact definition, the general contact surfaces are generated automatically by the GCGEN command based on physical parts and geometric shapes in the model. The program overlays contact surface elements (CONTA174) on 3D deformable bodies (on both lower- and higher-order elements); 3D contact line elements (CONTA177) on 3D beams, on feature edges of 3D deformable bodies, and on the perimeter edges of shell structures; and vertex-to-surface elements (CONTA175) on convex corners of 3D solid bodies and/or shell structures. The general contact definition may also contain target elements (TARGE170) overlaid on the surfaces of standalone rigid bodies.
GCGEN automatically assigns section IDs and element type IDs for each general contact surface. As a result, each general contact surface consists of contact or target elements that are easily identified by a unique section ID number. The real constant ID and material ID are always set to zero for contact and target elements in the general contact definition.
The program looks for contact interaction among all surfaces and within each surface. You can further control contact interactions between specific surfaces that could potentially be in contact via the GCDEF command. The material ID and real constant ID input on GCDEF identify interface properties (defined via MP or TB) and contact control parameters (defined via the R command) for a specific contact interaction. Unlike a pair-based contact definition, the contact and target elements in the general contact definition are not associated with these material and real constant ID numbers.
If both pair-based contact and general contact are defined in a model, the pair-based contact definitions are preserved, and the general contact definition automatically excludes overlapping interactions wherever pair-based contact exists.
Some element key options are not used or are set automatically for general contact. See the individual KEYOPT descriptions in "CONTA177 Input Summary" for details.
Contact Radius and Target Radius
To model beam-to-beam or edge-to-edge contact, the contact/target surface is assumed to be the surface of a cylinder. For a general beam cross section other than the circular section, you need to estimate an equivalent circular radius. Follow these guidelines to define the equivalent radius:
Determine the smallest cross section along the beam axis.
Determine the largest circle embedded in that cross section.
By default, the program models contact between exterior surfaces of two cylindrical beams for both pair-based and general contact. However, defining contact and target radii differs for the two contact methods.
In a pair-based contact definition, the associated TARGE170 target segment elements are either LINE or PARA segment types. Use the first real constant, R1, to define the equivalent radius on the target side, and use the second real constant, R2, to define the equivalent radius on the contact side. To specify internal beam-to-beam contact (a beam sliding inside a hollow beam, or pipe sliding inside another pipe), set KEYOPT(9) = 1 for the associated TARGE170 element type and input R1 as the inner radius of the outer beam (see Figure 177.2: Beam Sliding Inside a Hollow Beam).
In a general contact definition, the equivalent beam radius is specified via SECTYPE and SECDATA commands, as shown:
SECTYPE,SECID
,CONTACT,RADIUS ! Set Type = CONTACT and Subtype = RADIUS for user-defined contact radius SECDATA,VAL1
,VAL2
,VAL3
!VAL1
= equivalent outer radius (external beam-to-beam contact) !VAL2
= equivalent inner radius; also set VAL3 = 1 for internal beam-to-beam
To model internal beam-to-beam contact in a general contact definition,
specify VAL3
= 1 and input VAL2
equal to
the inner radius of the outer beam.
For the case of internal contact, the inner beam should usually be considered the contact surface and the outer beam should be the target surface. The inner beam can be considered as the target surface only when the inner beam is much stiffer than the outer beam.
Contact is detected when two circular beams touch or overlap each other. The non-penetration condition for beams with a circular cross section can be defined as follows.
For internal contact:
and for external contact:
where rc and rt are the radii of the cross sections of the beams on the contact and target sides, respectively; and d is the minimal distance between the two beams which also determines the contact normal direction (see Figure 177.4: Crossing Beams in Contact).[1] Contact occurs for negative values of g.
For beam-to-beam contact modeled with either pair-based contact or general contact, if the contact radius and/or target radius are not defined, the program automatically calculates the equivalent radius for each individual contact/target element based on the associated geometry of underlying elements. As a result, the equivalent radius may vary within a contact pair or within a general contact surface.
For rigid targets, the program cannot compute the target radius since underlying elements do not exist. In this case you must explicitly specify the target radius.
Friction
CONTA177 supports isotropic and orthotropic Coulomb friction. For isotropic friction, specify a single coefficient of friction, MU, using either TB,FRIC (recommended) or MP. For orthotropic friction, specify two coefficients of friction, MU1 and MU2, in two principal directions using TB,FRIC. (See Contact Friction in the Material Reference for more information.)
For isotropic friction, local element coordinates based on the nodal connectivity are used as principal directions. In the case of two crossing beams in contact, the first principal direction is defined by 1/2(t1 + t2 ). The first vector, t1 , points from the first contact node to the second contact node, and the second vector, t2 , points from the first target node to the second target node. In all other cases, the first principal direction points from node I to node J, and the second principal direction is defined by taking a cross product of the first principal direction and the contact normal.
For orthotropic friction, the principal directions are determined as follows. The global coordinate system is used by default, or you can define a local element coordinate system via the ESYS command. The first principal direction is defined by projecting the first direction of the chosen coordinate system onto the contact element. The second principal direction is defined by taking a cross product of the first principal direction and the contact normal. These directions also follow the rigid body rotation of the contact element to correctly model the directional dependence of friction. Be careful to choose the coordinate system (global or local) so that only the z direction parallels (or is nearly parallel to) the surface normal.
If you want to set the coordinate directions for isotropic friction (to the global Cartesian system or another system via ESYS), you can define orthotropic friction and set MU1 = MU2.
To define a coefficient of friction for isotropic or orthotropic friction that is dependent on temperature, time, normal pressure, sliding distance, or sliding relative velocity, use the TBFIELD command with TB,FRIC. See Contact Friction in the Material Reference for more information.
To implement a user-defined friction model, issue TB,FRIC with
TBOPT
= USER to specify friction properties and write a
USERFRIC
subroutine to calculate friction forces. (See Writing Your Own Friction Law (USERFRIC
) in the Contact Technology Guide for more
information. See also the Guide to User-Programmable Features in the Programmer's Reference for a
detailed description of USERFRIC
.)
Other Input
The contact interaction subroutine USERINTER
is available for
user-defined interface interactions, including interactions in the normal and tangential
directions. See Defining Your Own Contact Interaction (USERINTER
) in the Contact Technology Guide for more information on how to use this feature. See also the
Guide to User-Programmable Features in the Programmer's Reference for a detailed description of
the USERINTER
subroutine.
To model proper momentum transfer and energy balance between contact and target surfaces, impact constraints should be used in transient dynamic analysis. See the description of KEYOPT(7) below and the contact element discussion in the Mechanical APDL Theory Reference for details.
To model separation of bonded contact with KEYOPT(12) = 2, 3, 4, 5, or 6, use TB,CZM. See Debonding in the Contact Technology Guide for more information.
See the Contact Technology Guide for a detailed discussion on contact and using the contact elements. 3D Beam-to-Beam Contact (Pair-Based) and Line-to-Surface Contact (Pair-Based) discuss CONTA177 specifically, including the use of real constants and KEYOPTs.
The following table summarizes the element input. Element Input gives a general description of element input.
CONTA177 Input Summary
- Nodes
I, J, (K)
- Degrees of Freedom
UX, UY, UZ - Real Constants
R1, R2, FKN, FTOLN, ICONT, PINB, (Blank), (Blank), TAUMAX, CNOF, FKOP, FKT, COHE, (Blank), (Blank), (Blank), (Blank), (Blank), PZER, CZER, FACT, DC, SLTO, TNOP, TOLS, (Blank), (Blank), (Blank), COR, STRM FDMN, FDMT, (Blank), (Blank), TBND, (Blank) (Blank), (Blank), (Blank), (Blank), (Blank), (Blank), (Blank), (Blank), BSRL, KSYM, TFOR, TEND See Table 177.1: CONTA177 Real Constants for descriptions of the real constants. - Material Properties
TB command: See Element Support for Material Models for this element. MP command: MU, DMPR, DMPS - Special Features
- KEYOPTs
Presented below is a list of KEYOPTS available for this element. Included are links to sections in the Contact Technology Guide where more information is available on a particular topic.
- KEYOPT(1)
Selects degrees of freedom. Currently, the default (UX, UY, UZ) is the only valid option:
- 0 --
UX, UY, UZ
- KEYOPT(2)
Contact algorithm:
- 0 --
Augmented Lagrangian (default)
- 1 --
Penalty function
- 2 --
Multipoint constraint (MPC); see Multipoint Constraints and Assemblies in the Contact Technology Guide for more information
- 3 --
Lagrange multiplier on contact normal and penalty on tangent
- 4 --
Pure Lagrange multiplier on contact normal and tangent
For general contact, GCGEN sets KEYOPT(2) = 1 (penalty function) automatically.
- KEYOPT(3)
Contact model:
- 0 --
Exclude crossing beam-to-beam contact (contact force-based model) (default)
- 1 --
Exclude crossing beam-to-beam contact (contact traction-based model)
- 2 --
Include all scenarios: beam/edge to surface contact, parallel beam-to-beam contact, and crossing beam-to-beam contact (contact traction-based model)
- 3 --
Crossing beam-to-beam contact (contact traction-based model)
- 4 --
Crossing beam-to-beam contact (contact force-based model)
- KEYOPT(4)
Standard KEYOPT(4) Usage (line-to-surface contact only)
Location of contact detection point:
- 0 --
On nodal point - normal to target surface
- 3 --
On nodal point - normal to target surface (standard projection-based method)
- 4 --
On nodal point - normal from contact surface (dual shape function projection-based method)
Note: Some restrictions apply when the surface-projection-based method (KEYOPT(4) = 3 or 4) is defined. See KEYOPT(4) for more information.
Alternate KEYOPT(4) Usage with Surface-Based Constraints
To define a surface-based constraint, set KEYOPT(4) as follows:
- 0 --
Rigid surface constraint. The rigid surface constraint can be based on either the MPC approach (KEYOPT(2) = 2) or the Lagrange multiplier method (KEYOPT(2) = 3).
- 1 --
Force-distributed constraint. The force-distributed constraint can be based on either the MPC approach (KEYOPT(2) = 2) or the Lagrange multiplier method (KEYOPT(2) = 3).
- 3 --
Coupling constraint
For more information, see Surface-Based Constraints in the Contact Technology Guide.
- KEYOPT(5)
CNOF/ICONT Automated adjustment:
- 0 --
No automated adjustment
- 1 --
Close gap with auto CNOF
- 2 --
Reduce penetration with auto CNOF
- 3 --
Close gap/reduce penetration with auto CNOF
- 4 --
Auto ICONT
- KEYOPT(6)
Normal contact stiffness variation (used to enhance stiffness updating when KEYOPT(10) ≠ 1):
- 0 --
Use default range for stiffness updating
- 1 --
Make a nominal refinement to the allowable stiffness range
- 2 --
Make an aggressive refinement to the allowable stiffness range
- 3 --
Use an exponential pressure-penetration relationship
- KEYOPT(7)
Element level time incrementation control / impact constraints:
- 0 --
No control
- 1 --
Automatic bisection of increment
- 2 --
Change in contact predictions are made to maintain a reasonable time/load increment
- 3 --
Change in contact predictions made to achieve the minimum time/load increment whenever a change in contact status occurs
- 4 --
Use impact constraints for standard or rough contact (KEYOPT(12) = 0 or 1) in a transient dynamic analysis with automatic adjustment of time increment
Note: KEYOPT(7) = 4 is not supported for contact elements used in a general contact definition.
- KEYOPT(8)
Symmetric contact behavior:
- 0 --
Both symmetric pairs are active. However, each pair has its own contact characteristics.
- 1 --
Both symmetric pairs are active and have the same contact characteristics.
- 2 --
The program internally selects which asymmetric contact pair is used at the solution stage (used only when symmetric contact is defined). However, the contact stiffness of the active contact pair is influenced by the underlying element stiffness of the inactive pair.
- 3 --
The program internally selects which asymmetric contact pair is used at the solution stage (used only when symmetric contact is defined). The contact characteristics of the active contact pair are completely independent of the inactive pair.
KEYOPT(8) settings are ignored for asymmetric contact pairs and rigid-to-rigid contact pairs.
KEYOPT(8) is ignored for contact elements used in a general contact definition. Instead, issue GCDEF,AUTO to enable auto-asymmetric pairing logic.
- KEYOPT(9)
Effect of initial penetration or gap:
- 0 --
Include both initial geometrical penetration or gap and offset
- 1 --
Exclude both initial geometrical penetration or gap and offset
- 2 --
Include both initial geometrical penetration or gap and offset, but with ramped effects
- 3 --
Include offset only (exclude initial geometrical penetration or gap)
- 4 --
Include offset only (exclude initial geometrical penetration or gap), but with ramped effects
- 5 --
Include offset only (exclude initial geometrical penetration or gap) regardless of the initial contact status (near-field or closed)
- 6 --
Include offset only (exclude initial geometrical penetration or gap), but with ramped effects regardless of the initial contact status (near-field or closed)
The effects of KEYOPT(9) are dependent on settings for other KEYOPTs. The indicated initial gap effect is considered only if KEYOPT(12) = 4 or 5. See the discussion on using KEYOPT(9) in the Contact Technology Guide for more information.
KEYOPT(9) is not supported for contact elements used in a general contact definition. Instead, use TB,INTER and TBDATA,,
C1
to specify the effect of initial penetration or gap. If TBDATA,,C1
is not specified, the default for general contact is to exclude initial penetration/gap and offset. For more information, see Interaction Options for General Contact Definitions in the Material Reference.- KEYOPT(10)
Contact Stiffness Update:
- 0 --
Each iteration based on the current mean stress of underlying elements. The actual elastic slip never exceeds the maximum allowable limit (SLTO) during the entire solution.
- 1 --
Each load step if FKN is redefined during the load step.
- 2 --
Each iteration based on the current mean stress of underlying elements. The actual elastic slip does not exceed the maximum allowable limit (SLTO) within a substep.
For general contact, GCGEN sets KEYOPT(10) = 0 automatically.
- KEYOPT(11)
Shell/beam thickness effect (line-to-surface contact only):
- 0 --
Exclude
- 1 --
Include
Note: KEYOPT(11) = 1 is not valid when the underlying elements are part of a pre-integrated shell section (SECTYPE,,GENS).
Note: In the case of general contact, GCGEN sets KEYOPT(11) = 1 automatically for beam-to-surface contact.
- KEYOPT(12)
Behavior of contact surface:
- 0 --
Standard
- 1 --
Rough
- 2 --
No separation (sliding permitted)
- 3 --
Bonded
- 4 --
No separation (always)
- 5 --
Bonded (always)
- 6 --
Bonded (initial contact)
When KEYOPT(12) = 5 or 6 is used with the MPC algorithm to model surface-based constraints, the KEYOPT(12) setting will have an affect on the local coordinate system of the contact element nodes. See Specifying a Local Coordinate System in the Contact Technology Guide for more information.
KEYOPT(12) is not supported for contact elements used in a general contact definition. Instead, issue TB,INTER with the appropriate
TBOPT
label to specify the behavior at the contact surface. For more information, see Interaction Options for General Contact Definitions in the Material Reference.- KEYOPT(13)
Tangential contact stiffness variation only for frictional contact.
- 0 --
Use the default range for stiffness updating.
- 1 --
Make an aggressive refinement to the allowable stiffness range.
- KEYOPT(14)
Number of target segments interacting with each contact detection point:
- 0 --
Only one target segment
- 1 --
Up to four target segments
- 2 --
Up to eight target segments
In the case of general contact, GCGEN sets KEYOPT(14) = 1 automatically for beam-to-beam contact.
- KEYOPT(15)
Effect of contact stabilization damping:
- 0 --
Damping is activated only in the first load step (default).
- 1 --
Deactivate automatic damping.
- 2 --
Damping is activated for all load steps.
- 3 --
Damping is activated at all times regardless of the contact status of previous substeps.
Note: Normal stabilization damping is only applied to the contact element when the current contact status of the contact detection point is near-field. When KEYOPT(15) = 0, 1, or 2, normal stabilization damping is not applied in the current substep if any contact detection point has a closed status. However, when KEYOPT(15) = 3, normal stabilization damping is always applied as long as the current contact status of the contact detection point is near-field. Tangential stabilization damping is automatically activated when normal damping is activated. Tangential damping can also be applied independent of normal damping for sliding contact. See Applying Contact Stabilization Damping in the Contact Technology Guide for more information.
- KEYOPT(18)
Sliding behavior:
- 0 --
Finite sliding (default). The contacting interface can undergo separation, relative large sliding, and arbitrary rotation.
- 1 --
Small sliding. The contacting interface can undergo only small sliding during the entire solution; arbitrary rotation is permitted.
- 2 --
Adaptive small sliding. The contact interface can undergo either small sliding or finite sliding within each substep based on the contact status at the beginning of the substep. If the contact status is closed, small sliding is used.
Table 177.1: CONTA177 Real Constants
No. | Name | Description | For more information, see this section in the Contact Technology Guide . . . |
---|---|---|---|
1 | R1 |
Target radius for cylinder, cone, or sphere, used for line-to-surface contact (rigid target) | |
Target radius for beam-to-beam contact | Real constants R1, R2 | ||
2 | R2 |
Target radius at second node of cone, used for line-to-surface contact (rigid target) | |
Contact radius for beam-to-beam contact | Real constants R1, R2 | ||
3 | FKN[1] | ||
4 | FTOLN |
Penetration tolerance factor [6] | |
5 | ICONT |
Initial contact closure | |
6 | PINB |
Pinball region | or |
7 | PZER | Pressure at zero penetration [2] [3] | Pressure-Penetration Relationship (KEYOPT(6) = 3) |
8 | CZER | Initial contact clearance | Pressure-Penetration Relationship (KEYOPT(6) = 3) |
9 | TAUMAX | ||
10 | CNOF | ||
11 | FKOP | ||
12 | FKT[1] | ||
13 | COHE |
Contact cohesion | |
21 | FACT |
Static/dynamic ratio | |
22 | DC |
Exponential decay coefficient | |
23 | SLTO |
Allowable elastic slip | |
24 | TNOP | ||
25 | TOLS |
Target edge extension factor | |
29 | COR |
Coefficient of restitution | |
30 | STRM |
Load step number for ramping penetration or Starting time for contact stiffness ramping | |
31 | FDMN | Normal stabilization damping factor [2] [3] | |
32 | FDMT | Tangential stabilization damping factor [2] [3] | |
35 | TBND | Critical bonding temperature [2] [3] | |
45 | BSRL | Original contact pair real constant ID (after contact splitting) | Real Constant Set IDs for Split Pairs |
46 | KSYM | Real constant ID of the associated companion pair for symmetric contact or self contact (after contact splitting) | Real Constant Set IDs for Split Pairs |
47 | TFOR | Pair-based force tolerance | |
48 | TEND | Ending time for ramping contact stiffness | Modeling Interference Fit |
For the contact force-based model (KEYOPT(3) = 0 or 4), the units of real constants FKN and FKT have a factor of AREA with respect to those used in the surface-to-surface contact elements. See Performing a 3D Line-to-Surface Contact Analysis for more information.
This real constant can be defined as a function of certain primary variables.
This real constant can be defined by the user subroutine USERCNPROP.F.
For the contact force-based model (KEYOPT(3) = 0 or 4), TNOP is the allowable tensile contact force. For the contact traction-based model (KEYOPT(3) = 1, 2, or 3), TNOP is the allowable tensile contact pressure.
When CONTA177 is used as part of a forced-distributed constraint and KEYOPT(7) = 2 on the TARGE170 element, FKN is used to define weighting factors in tabular format with node number as the primary variable.
When the relaxation option is enabled (KEYOPT(11) = 1 on the TARGE170 element), FKN and FKT are translational relaxation coefficient and rotational relaxation coefficient, respectively, and tabular input is not supported. In addition, FTOLN and TNOP are translational tolerance and rotational tolerance, respectively.
CONTA177 Output Data
The solution output associated with the element is in two forms:
Nodal displacements included in the overall nodal solution
Additional element output as shown in Table 177.2: CONTA177 Element Output Definitions.
A general description of solution output is given in Solution Output. See the Basic Analysis Guide for ways to view results.
The Element Output Definitions table uses the following notation:
A colon (:) in the Name column indicates that the item can be accessed by the Component Name method (ETABLE, ESOL). The O column indicates the availability of the items in the file jobname.out. The R column indicates the availability of the items in the results file.
In either the O or R columns, “Y” indicates that the item is always available, a letter or number refers to a table footnote that describes when the item is conditionally available, and “-” indicates that the item is not available.
Table 177.2: CONTA177 Element Output Definitions
Name | Definition | O | R |
---|---|---|---|
EL | Element Number | Y | Y |
NODES | Nodes I, J, K | Y | Y |
XC, YC, ZC | Location where results are reported (same as nodal location) | Y | Y |
TEMP | Temperature T(I) | Y | Y |
VOLU | Length | Y | Y |
NPI | Number of integration points | Y | - |
ITRGET | Target surface number (assigned by the program) | Y | - |
ISOLID [14] | Underlying beam or shell element number | Y | - |
CONT:STAT | Current contact statuses | 1 | 1 |
OLDST | Old contact statuses | 1 | 1 |
ISEG | Current contacting target element number | Y | Y |
OLDSEG | Underlying old target number | Y | - |
CONT:PENE | Current penetration (gap = 0; penetration = positive value) | Y | Y |
CONT:GAP | Current gap (gap = negative value; penetration = 0) | Y | Y |
NGAP | New or current gap at current converged substep (gap = negative value; penetration = positive value) | Y | - |
OGAP | Old gap at previously converged substep (gap = negative value; penetration = positive value) | Y | - |
IGAP | Initial gap at start of current substep (gap = negative value; penetration = positive value) | Y | Y |
GGAP | Geometric gap at current converged substep (gap = negative value; penetration = positive value) | - | Y |
CONT:PRES | Normal contact force/pressure | 2 | 2 |
TAUR/TAUS [7] | Tangential contact forces/stresses | 2 | 2 |
KN | Current normal contact stiffness (units: FORCE/LENGTH for contact force model, FORCE/LENGTH3 for contact traction model) | 5 | 5 |
KT | Current tangent contact stiffness (same units as KN) | 5 | 5 |
MU [8] | Friction coefficient | Y | Y |
TASS/TASR [7] | Total (algebraic sum) sliding in S and R directions | 3 | 3 |
AASS/AASR [7] | Total (absolute sum) sliding in S and R directions | 3 | 3 |
TOLN | Penetration tolerance | Y | Y |
CONT:SFRIC | Frictional force/stress, SQRT (TAUR**2+TAUS**2) | 2 | 2 |
CONT:STOTAL | Total force/stress, SQRT (PRES**2+TAUR**2+TAUS**2) | 2 | 2 |
CONT:SLIDE | Amplitude of total accumulated sliding, SQRT (TASS**2+TASR**2) | 3 | 3 |
FDDIS | Frictional energy dissipation rate | 6 | 6 |
ELSI | Total equivalent elastic slip distance | - | Y |
PLSI | Total (accumulated) equivalent plastic slip due to frictional sliding | - | Y |
GSLID | Amplitude of total accumulated sliding (including near-field) | - | 9 |
VREL | Equivalent sliding velocity (slip rate) | - | Y |
DBA | Penetration variation | Y | Y |
PINB | Pinball Region | - | Y |
CONT:CNOS | Total number of contact status changes during substep | Y | Y |
TNOP | Maximum allowable tensile contact force/pressure | 2 | 2 |
SLTO | Allowable elastic slip | Y | Y |
CAREA | Contacting area | - | Y |
R1 | Target radius for beam-to-beam contact | - | Y |
R2 | Contact radius for beam-to-beam contact | - | Y |
DTSTART | Load step time during debonding | Y | Y |
DPARAM | Debonding parameter | Y | Y |
DENERI [12] | Energy released due to separation in normal direction - mode I debonding | Y | Y |
DENERII [12] | Energy released due to separation in tangential direction - mode II debonding | Y | Y |
DENER [13] | Total energy released due to debonding | Y | Y |
CNFX [10] | Contact element force - X component | - | 4 |
CNFY [10] | Contact element force - Y component | - | 4 |
CNFZ [10] | Contact element force - Z component | - | 4 |
CNTX [11] | Contact element force due to tangential stresses - X component | - | 4 |
CNTY [11] | Contact element force due to tangential stresses - Y component | - | 4 |
CNTZ [11] | Contact element force due to tangential stresses - Z component | - | 4 |
SDAMP | Stabilization damping coefficient | - | Y |
The possible values of STAT and OLDST are:
0 = Open and not near contact 1 = Open but near contact 2 = Closed and sliding 3 = Closed and sticking For the force-based model (KEYOPT(3) = 0 or 4), the unit of this quantity is FORCE. For the traction-based model (KEYOPT(3) = 1, 2, 3), the unit is FORCE/AREA.
Only accumulates the sliding for sliding and closed contact (STAT = 2,3).
Contact element forces are defined in the global Cartesian system
For the force-based model (KEYOPT(3) = 0 or 4), the unit of stiffness is FORCE/LENGTH. For the traction-based model (KEYOPT(3) = 1, 2, 3), the unit is FORCE/LENGTH3.
FDDIS = (contact friction stress)*(sliding distance of substep)/(time increment of substep)
For the case of orthotropic friction in contact between beams (or shell edges) and a 3D surface, components are defined in the global Cartesian system.
For orthotropic friction, an equivalent coefficient of friction is output.
Accumulated sliding distance for near-field, sliding, and closed contact (STAT = 1,2,3).
The contact element force values (CNFX, CNFY, CNFZ) are calculated based on the individual contact element plus the surrounding contact elements. Therefore, the contact force values may not equal the contact element area times the contact pressure (CAREA * PRES).
CNTX, CNTY, and CNTZ report the total contact element forces due to tangential stresses. Since CNFX, CNFY, and CNFZ report the total contact element forces, the contact element forces due to normal pressure are (CNFX-CNTX), (CNFY-CNTY), and (CNFZ-CNTZ).
DENERI and DENERII are available only when one of the following material models is used: TB,CZM,,,,CBDD or TB,CZM,,,,CBDE.
DENER is available only when one of the following material models is used: TB,CZM,,,,BILI or TB,CZM,,,,EXPO.
The underlying element (ISOLID) can be obtained by a *GET command following a CNCHECK command. See Designating Contact and Target Surfaces for detail.
Contact results (including all element results) are generally not reported for elements that have a status of "open and not near contact" (far-field).
Table 177.3: CONTA177 (3D) Item and Sequence Numbers lists output available via the ETABLE command using the Sequence Number method. See Creating an Element Table in the Basic Analysis Guide and The Item and Sequence Number Table in this reference for more information. The following notation is used:
- Name
output quantity as defined in Table 177.2: CONTA177 Element Output Definitions
- Item
predetermined item label for ETABLE
- E
sequence number for single-valued or constant element data
- I, J, K
sequence number for data at nodes I, J, K (contact results for line-to-surface and parallel beams in contact)
- IP
sequence number for data of crossing beams in contact (contact results at intersection point)
Table 177.3: CONTA177 (3D) Item and Sequence Numbers
Output Quantity Name | ETABLE and ESOL Command Input | |||||
---|---|---|---|---|---|---|
Item | E | I | J | K | IP | |
PRES | SMISC | 13 | 1 | 2 | 3 | 4 |
TAUR | SMISC | - | 5 | 6 | 7 | 8 |
TAUS | SMISC | - | 9 | 10 | 11 | 12 |
FDDIS | SMISC | - | 18 | 19 | 20 | 21 |
STAT [1] | NMISC | 41 | 1 | 2 | 3 | 4 |
OLDST | NMISC | - | 5 | 6 | 7 | 8 |
PENE [2] | NMISC | - | 9 | 10 | 11 | 12 |
DBA | NMISC | - | 13 | 14 | 15 | 16 |
TASR | NMISC | - | 17 | 18 | 19 | 20 |
TASS | NMISC | - | 21 | 22 | 23 | 24 |
KN | NMISC | - | 25 | 26 | 27 | 28 |
KT | NMISC | - | 29 | 30 | 31 | 32 |
TOLN | NMISC | - | 33 | 34 | 35 | 36 |
IGAP | NMISC | - | 37 | 38 | 39 | 40 |
PINB | NMISC | 42 | - | - | - | |
CNFX | NMISC | 43 | - | - | - | |
CNFY | NMISC | 44 | - | - | - | |
CNFZ | NMISC | 45 | - | - | - | |
CNTX | NMISC | 186 | - | - | - | - |
CNTY | NMISC | 187 | - | - | - | - |
CNTZ | NMISC | 188 | - | - | - | - |
ISEG [3] | NMISC | - | 46 | 47 | 48 | 49 |
AASR | NMISC | - | 50 | 51 | 52 | 53 |
AASS | NMISC | - | 54 | 55 | 56 | 57 |
CAREA | NMISC | 58 | - | - | - | |
MU | NMISC | - | 62 | 63 | 64 | 65 |
DTSTART | NMISC | - | 66 | 67 | 68 | 69 |
DPARAM | NMISC | - | 70 | 71 | 72 | 73 |
CNOS | NMISC | - | 112 | 113 | 114 | 115 |
TNOP | NMISC | - | 116 | 117 | 118 | 119 |
SLTO | NMISC | - | 120 | 121 | 122 | 123 |
ELSI | NMISC | - | 136 | 137 | 138 | 139 |
DENERI or DENER | NMISC | - | 140 | 141 | 142 | 143 |
DENERII | NMISC | - | 144 | 145 | 146 | 147 |
GGAP | NMISC | - | 152 | 153 | 154 | 155 |
VREL | NMISC | - | 156 | 157 | 158 | 159 |
SDAMP | NMISC | - | 160 | 161 | 162 | 163 |
PLSI | NMISC | - | 164 | 165 | 166 | 167 |
GSLID | NMISC | - | 168 | 169 | 170 | 171 |
R1 | NMISC | - | 172 | 173 | 174 | 175 |
R2 | NMISC | - | 176 | 177 | 178 | 179 |
Element Status = highest value of status of integration points within the element
The floating point output format for large integers may lead to incorrect ISEG values. It is good practice to verify the NMISC values via the *GET command. For example, *GET,
Par
,ELEM,N
,NMISC,46 returns the ISEG value for node I of elementN
.
You can display or list contact results through several POST1 postprocessor commands. The contact-specific items for PLNSOL, PLESOL, PRNSOL, and PRESOL are:
STAT | Contact status |
PENE | Contact penetration |
PRES | Contact pressure for the traction-based model. Contact normal force for the force-based model. |
SFRIC | Contact friction stress for the traction-based model. Friction force for the force-based model. |
STOT | Contact total stress (pressure plus friction) for the traction-based model. Total contact force for the force-based model. |
SLIDE | Contact sliding distance |
GAP | Contact gap distance |
CNOS | Total number of contact status changes during substep |
Contact results (I, J, K columns in Table 177.3: CONTA177 (3D) Item and Sequence Numbers) are reported at contact nodes for 3D beam-to-surface and parallel beam-to-beam contact. Contact results from crossing beam-to-beam contact (IP column in Table 177.3: CONTA177 (3D) Item and Sequence Numbers) are reported at an intersection point of two crossing beams. When contact from crossing beams is detected, the associated contact pressure (PRES), the contact frictional stress (SFRIC), and the total stress (STOT) are superimposed on each nodes. Maximum values are reported for other contact results (STAT, PENE, SLIDE, GAP, CNOS).
CONTA177 Assumptions and Restrictions
For line-to-surface contact, the thickness effects of underlying beam elements on the contact side and underlying shell elements on the target side can be taken into account by setting KEYOPT(11) = 1.
This element is nonlinear and requires a full Newton iterative solution, regardless of whether large or small deflections are specified. An exception to this is when MPC bonded contact is specified (KEYOPT(2) = 2 and KEYOPT(12) = 5 or 6).
The normal contact stiffness factor (FKN) must not be so large as to cause numerical instability.
FTOLN, PINB, and FKOP can be changed between load steps or during restart stages.
You can use this element in nonlinear static or nonlinear full transient analyses.
In addition, you can use it in modal analyses, eigenvalue buckling analyses, and harmonic analyses. For these analysis types, the program assumes that the initial status of the element (that is, the status at the completion of the static prestress analysis, if any) does not change.
In a distributed-memory parallel (DMP) processing run, multiple target segments can interact with this element's contact detection points (KEYOPT(14) > 0) only when the pure mesh-based domain decomposition method or the frequency-based domain composition method is used. KEYOPT(14) > 0 is not supported when both the frequency-based and mesh-based domain decomposition methods are used together in a harmonic analysis (for example, DDOPTION,FREQ,2). See DDOPTION for more information about the domain decomposition methods used in DMP.
Only damping defined via KEYOPT(15) is supported. All other damping specifications (Rayleigh damping, DMPSTR, and so on) are not supported.
Certain contact features are not supported when this element is used in a general contact definition. For details, see General Contact in the Contact Technology Guide.
CONTA177 Product Restrictions
When used in the product(s) listed below, the stated product-specific restrictions apply to this element in addition to the general assumptions and restrictions given in the previous section.
Ansys Mechanical Pro —
Birth and death is not available.
Debonding is not available.
User-defined contact is not available.
User-defined friction is not available.
Ansys Mechanical Premium —
Debonding is not available.
User-defined contact is not available.
User-defined friction is not available.