www.kxcad.net Home > CAE Index > ANSYS Index > Release 11.0 Documentation for ANSYS
CONTA176 is used to represent contact and sliding between 3-D line segments (TARGE170) and a deformable line segment, defined by this element. The element is applicable to 3-D beam-beam structural contact analyses. This element is located on the surfaces of 3-D beam or pipe elements with or without midside nodes (BEAM4, BEAM24, BEAM188, BEAM189, PIPE16, PIPE20). Contact occurs when the element surface penetrates one of the 3-D straight line or parabolic line segment elements (TARGE170) on a specified target surface. Coulomb and shear stress friction is allowed. This element also allows separation of bonded contact to simulate interface delamination. See CONTA176 in the Theory Reference for ANSYS and ANSYS Workbench for more details about this element. To model beam-to-surface contact, use the line-to-surface contact element, CONTA177.
The geometry and node locations are shown in Figure 176.1: "CONTA176 Geometry". The element is defined by two nodes (if the underlying beam element does not have a midside node) or three nodes (if the underlying beam 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. See Generating Contact Elements in the Contact Technology Guide for more information on generating elements automatically using the ESURF command.
Three different scenarios can be modeled by CONTA176:
Internal contact where one beam (or pipe) slides inside another hollow beam (or pipe) (see Figure 176.2: "Beam Sliding Inside a Hollow Beam")
External contact between two beams that lie next to each other and are roughly parallel (see Figure 176.3: "Parallel Beams in Contact")
External contact between two beams that cross (see Figure 176.4: "Crossing Beams in Contact")
Use KEYOPT(3) = 0 for the first two scenarios (internal contact and parallel beams). In both cases, the contact condition is only checked at contact nodes.
Use KEYOPT(3) = 1 for the third scenario (beams that cross). In this case, the contact condition is checked along the entire length of the beams. The beams with circular cross sections are assumed to come in contact in a point-wise manner. Each contact element can potentially contact no more than one target element.
The 3-D line-to-line contact elements are associated with the target line segment elements (LINE or PARA segment types for TARGE170) via a shared real constant set. The contact/target surface is assumed to be the surface of a cylinder. For a general beam cross section, use an equivalent circular beam (see Figure 176.5: "Equivalent Circular Cross Section"). Use the first real constant, R1, to define the radius on the target side (target radius rt). Use the second real constant, R2, to define the radius on the contact side (contact radius rc). Follow these guidelines to define the equivalent circular cross section:
Determine the smallest cross section along the beam axis.
Determine the largest circle embedded in that cross section.
The target radius can be entered as either a negative or positive value. Use a negative value when modeling internal contact (a beam sliding inside a hollow beam, or pipe sliding inside another pipe), with the input value equal to the inner radius of the outer beam (see Figure 176.2: "Beam Sliding Inside a Hollow Beam"). Use a positive value when modeling contact between the exterior surfaces of two cylindrical beams.
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 176.4: "Crossing Beams in Contact"). Contact occurs for negative values of g.
ANSYS looks for contact only between contact and target surfaces with the same real constant set. For either rigid-flexible or flexible-flexible contact, one of the deformable surfaces must be represented by a contact surface. See Designating Contact and Target Surfaces in the Contact Technology Guide for more information. If more than one target surface will make contact with the same boundary of beam elements, you must define several contact elements that share the same geometry but relate to separate targets (targets which have different real constant numbers), or you must combine the two target surfaces into one (targets that share the same real constant numbers).
CONTA176 supports isotropic and orthotropic Coulomb friction. For isotropic friction, specify a single coefficient of friction, MU, using either TB command input (recommended) or the MP command. For orthotropic friction, specify two coefficients of friction, MU1 and MU2, in two principal directions using TB command input. (See Contact Friction for more information.) The local element coordinates based on the nodal connectivity are used as principal directions. Local element coordinates defined using the ESYS command are ignored.
To model separation of bonded contact with KEYOPT(12) = 4, 5, or 6, use the TB command with the CZM label. 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. " 3-D Beam-to-Beam Contact" discusses CONTA176 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.
I, J, (K)
| UX, UY, UZ |
| R1, R2, FKN, FTOLN, ICONT, PINB, |
| PMAX, PMIN, TAUMAX, CNOF, FKOP, FKT, |
| COHE, (Blank), (Blank), (Blank), (Blank), (Blank), |
| (Blank), (Blank), FACT, DC, SLTO, TNOP, |
| TOLS, (Blank) |
| See Table 176.1: "CONTA176 Real Constants" for descriptions of the real constants. |
| DAMP, MU (MP command) |
| FRIC (TB command; see Contact Friction) |
| CZM (TB command; see Cohesive Zone Materials Used for Debonding in the Contact Technology Guide) |
| Nonlinear |
| Large deflection |
| Isotropic or orthotropic friction |
| Debonding |
| Birth and death |
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.
Selects degrees of freedom. Currently, the default (UX, UY, UZ) is the only valid option:
UX, UY, UZ
Contact algorithm:
Augmented Lagrangian (default)
Penalty function
Multipoint constraint (MPC); see "Multipoint Constraints and Assemblies" in the Contact Technology Guide for more information
Lagrange multiplier on contact normal and penalty on tangent
Pure Lagrange multiplier on contact normal and tangent
Beam contact type:
Parallel beams or beam inside beam
Crossing beams
Type of surface-based constraint (see Surface-based Constraints for more information):
Rigid surface constraint
Force-distributed constraint
CNOF/ICONT Automated adjustment:
No automated adjustment
Close gap with auto CNOF
Reduce penetration with auto CNOF
Close gap/reduce penetration with auto CNOF
Auto ICONT
Contact stiffness variation (used to enhance stiffness updating when KEYOPT(10) > 0):
Use default range for stiffness updating
Make a nominal refinement to the allowable stiffness range
Make an aggressive refinement to the allowable stiffness range
Element level time incrementation control:
No control
Automatic bisection of increment
Change in contact predictions are made to maintain a reasonable time/load increment
Change in contact predictions made to achieve the minimum time/load increment whenever a change in contact status occurs
For KEYOPT(7) = 2 or 3, includes automatic bisection of increment. It is activated only if SOLCONTROL,ON,ON is issued at the procedure level.
Asymmetric contact selection:
No action
ANSYS internally selects which asymmetric contact pair is used at the solution stage (used only when symmetry contact is defined).
Effect of initial penetration or gap:
Include both initial geometrical penetration or gap and offset
Exclude both initial geometrical penetration or gap and offset
Include both initial geometrical penetration or gap and offset, but with ramped effects
Include offset only (exclude initial geometrical penetration or gap)
Include offset only (exclude initial geometrical penetration or gap), but with ramped effects
For KEYOPT(9) = 1, 3, or 4, the indicated initial gap effect is considered only if KEYOPT(12) = 4 or 5.
Contact Stiffness Update:
Each load step if FKN is redefined during load step (pair based).
Each substep based on mean stress of underlying elements from the previous substep (pair based).
Each iteration based on current mean stress of underlying elements (pair based).
Each load step if FKN is redefined during load step (individual element based).
Each substep based on mean stress of underlying elements from the previous substep (individual element based).
Each iteration based on current mean stress of underlying elements (individual element based).
KEYOPT(10) = 0, 1, and 2 are pair based, meaning that the stiffness and settings for ICONT, FTOLN, PINB, PMAX, and PMIN are averaged across all the contact elements in a contact pair. For KEYOPT(10) = 3, 4, and 5, the stiffness and settings are based on each individual contact element (geometry and material behaviors).
Behavior of contact surface:
Standard
Rough
No separation (sliding permitted)
Bonded
No separation (always)
Bonded (always)
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 impact on the local coordinate system of the contact element nodes. See Specifying a Local Coordinate System in the Contact Technology Guide for more information.
Table 176.1 CONTA176 Real Constants
| No. | Name | Description | For more information, see this section in the Contact Technology Guide . . . |
|---|---|---|---|
| 1 | R1 | Target radius | |
| 2 | R2 | Contact radius | |
| 3 | FKN[1] | Normal penalty stiffness factor | |
| 4 | FTOLN | Penetration tolerance factor | |
| 5 | ICONT | Initial contact closure | |
| 6 | PINB | Pinball region | or |
| 7 | PMAX | Upper limit of initial allowable penetration | |
| 8 | PMIN | Lower limit of initial allowable penetration | |
| 9 | TAUMAX | Maximum friction stress | |
| 10 | CNOF | Contact surface offset | |
| 11 | FKOP | Contact opening stiffness or contact damping | |
| 12 | FKT[1] | Tangent penalty stiffness factor | |
| 13 | COHE | Contact cohesion | |
| 21 | FACT | Static/dynamic ratio | |
| 22 | DC | Exponential decay coefficient | |
| 23 | SLTO | Allowable elastic slip | |
| 24 | TNOP | Maximum allowable tensile contact force | |
| 25 | TOLS | Target edge extension factor |
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 3-D Beam-to-Beam Contact Analysis for more information.
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 176.2: "CONTA176 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 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 number refers to a table footnote that describes when the item is conditionally available, and a - indicates that the item is not available.
Table 176.2 CONTA176 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 ANSYS) | Y | - |
| ISOLID | Underlying beam element number | Y | - |
| CONT:STAT | Current contact statuses | 1 | 1 |
| OLDST | Old contact statuses | 1 | 1 |
| ISEG | Underlying current target 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 |
| CONT:PRES | Normal contact force | 2 | 2 |
| TAUR/TAUS[7] | Tangential contact stresses | 2 | 2 |
| KN | Current normal contact stiffness (units: Force/Length) | 5 | 5 |
| KT | Current tangent contact stiffness (same units as KN) | 5 | 5 |
| MU[8] | Friction coefficient | 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 stress SQRT (TAUR**2+TAUS**2) | 2 | 2 |
| CONT:STOTAL | Total stress SQRT (PRES**2+TAUR**2+TAUS**2) | 2 | 2 |
| CONT:SLIDE | Total sliding SQRT (TASS**2+TASR**2) | Y | Y |
| DBA | Penetration variation | Y | Y |
| PINB | Pinball Region | - | Y |
| CNFX[4] | Contact element force-X component | - | Y |
| CNFY | Contact element force-Y component | - | Y |
| CNFZ | Contact element force-Z component | - | Y |
| CAREA | Contacting area | - | Y |
| FDDIS | Frictional energy dissipation | 6 | 6 |
| CNOS | Total number of contact status changes during substep | Y | Y |
| TNOP | Maximum allowable tensile contact force | 2 | 2 |
| SLTO | Allowable elastic slip | Y | Y |
| ELSI | Elastic slip distance for sticking contact within a substep | - | Y |
| DTSTART | Load step time during debonding | Y | Y |
| DPARAM | Debonding parameter | Y | Y |
| DENERI | Energy released due to separation in normal direction - mode I debonding | Y | Y |
| DENERII | Energy released due to separation in tangential direction - mode II debonding | Y | 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 |
Contact element forces are defined in the global Cartesian system
FDDIS = (contact friction stress)*(sliding distance of substep)/(time increment of substep)
For the case of orthotropic friction in contact between beams, components are defined in the global Cartesian system.
For orthotropic friction, an equivalent coefficient of friction is output.
The following table lists output available through 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 manual for more information.
output quantity as defined in Table 176.2: "CONTA176 Element Output Definitions"
predetermined item label for ETABLE command
sequence number for single-valued or constant element data
sequence number for data at nodes I, J, K
Table 176.3 CONTA176 (3-D) Item and Sequence Numbers
| Output Quantity Name | ETABLE and ESOL Command Input | ||||
|---|---|---|---|---|---|
| Item | E | I | J | K | |
| PRES | SMISC | 13 | 1 | 2 | 3 |
| TAUR | SMISC | - | 5 | 6 | 7 |
| TAUS | SMISC | - | 9 | 10 | 11 |
| FDDIS | SMISC | - | 18 | 19 | 20 |
| STAT[1] | NMISC | 41 | 1 | 2 | 3 |
| OLDST | NMISC | - | 5 | 6 | 7 |
| PENE[2] | NMISC | - | 9 | 10 | 11 |
| DBA | NMISC | - | 13 | 14 | 15 |
| TASR | NMISC | - | 17 | 18 | 19 |
| TASS | NMISC | - | 21 | 22 | 23 |
| KN | NMISC | - | 25 | 26 | 27 |
| KT | NMISC | - | 29 | 30 | 31 |
| TOLN | NMISC | - | 33 | 34 | 35 |
| IGAP | NMISC | - | 37 | 38 | 39 |
| PINB | NMISC | 42 | - | - | - |
| CNFX | NMISC | 43 | - | - | - |
| CNFY | NMISC | 44 | - | - | - |
| CNFZ | NMISC | 45 | - | - | - |
| ISEG | NMISC | - | 46 | 47 | 48 |
| AASR | NMISC | - | 50 | 51 | 52 |
| AASS | NMISC | - | 54 | 55 | 56 |
| CAREA | NMISC | - | 58 | 59 | 60 |
| MU | NMISC | - | 62 | 63 | 64 |
| DTSTART | NMISC | - | 66 | 67 | 68 |
| DPARAM | NMISC | - | 70 | 71 | 72 |
| CNOS | NMISC | - | 112 | 113 | 114 |
| TNOP | NMISC | - | 116 | 117 | 118 |
| SLTO | NMISC | - | 120 | 121 | 122 |
| ELSI | NMISC | - | 136 | 137 | 138 |
| DENERI | NMISC | - | 140 | 141 | 142 |
| DENERII | NMISC | - | 144 | 145 | 146 |
The main restriction is the assumption of constant circular beam cross section. The contact radius is assumed to be the same for all elements in the contact pair.
For KEYOPT(3) = 1 (crossing beams), contact between the beams is pointwise, and each contact element contacts no more than one target element.
This element is nonlinear and requires a full Newton iterative solution, regardless of whether large or small deflections are specified.
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.
The value of FKN can be smaller when combined with the Augmented Lagrangian method, for which TOLN must be used.
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 (i.e., the status at the completion of the static prestress analysis, if any) does not change.
This element allows birth and death and will follow the birth and death status of the underlying beam, pipe, or target elements.