www.kxcad.net Home > CAE Index > ANSYS Index > Release 11.0 Documentation for ANSYS


Your Ad Here

CONTA173

3-D 4-Node Surface-to-Surface Contact

MP ME ST PR PRN DS <> <> <> <> <> PP <>

Product Restrictions

CONTA173 Element Description

CONTA173 is used to represent contact and sliding between 3-D "target" surfaces (TARGE170) and a deformable surface, defined by this element. The element is applicable to 3-D structural and coupled field contact analyses. This element is located on the surfaces of 3-D solid or shell elements without midside nodes (SOLID5, SOLID45, SOLID46, SOLID65, SOLID69, SOLID70, SOLID96, SOLID185, SOLSH190, VISCO107, SHELL28, SHELL41, SHELL43, SHELL57, SHELL63, SHELL131, SHELL157, SHELL181, and MATRIX50). It has the same geometric characteristics as the solid or shell element face with which it is connected (see Figure 173.1: "CONTA173 Geometry"). Contact occurs when the element surface penetrates one of the target 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 CONTA173 in the Theory Reference for ANSYS and ANSYS Workbench for more details about this element. Other surface-to-surface contact elements (CONTA171, CONTA172, CONTA174) are also available.

Figure 173.1  CONTA173 Geometry

R = Element x-axis for isotropic friction

xo = Element axis for orthotropic friction if ESYS is not supplied (parallel to global X-axis)

x = Element axis for orthotropic friction if ESYS is supplied

CONTA173 Input Data

The geometry and node locations are shown in Figure 173.1: "CONTA173 Geometry". The element is defined by four nodes (the underlying solid or shell element has no midside nodes). If the underlying solid or shell elements do have midside nodes, use CONTA174. The node ordering is consistent with the node ordering for the underlying solid or shell element. The positive normal is given by the right-hand rule going around the nodes of the element and is identical to the external normal direction of the underlying solid or shell element surface. For shell elements, the same nodal ordering between shell and contact elements defines upper surface contact; otherwise, it represents bottom surface contact. Remember the target surfaces must always be on its outward normal direction. See Generating Contact Elements in the Contact Technology Guide for more information on generating elements automatically using the ESURF command.

CONTA173 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.)

For isotropic friction, the applicable coordinate system is the default element coordinate system (noted by the R and S axes in the above figure).

For orthotropic friction, the principal directions are determined as follows. The global coordinate system is used by default, or you may define a local element coordinate system with the ESYS command. (These are depicted by the xoand x axes in the above figure.) The first principal direction is defined by projecting the first direction of the chosen coordinate system onto the contact surface. 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 the first direction of that system is within 45° of the tangent to the contact surface.

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.

The 3-D contact surface elements are associated with the 3-D target segment elements (TARGE170) via a shared real constant set. ANSYS looks for contact only between 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 solid 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).

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.

A summary of the element input is given in "CONTA173 Input Summary". A general description of element input is given in Element Input.

CONTA173 Input Summary

Nodes

I, J, K, L

Degrees of Freedom
UX, UY, UZ (if KEYOPT(1) = 0)
UX, UY, UZ, TEMP (if KEYOPT(1) = 1)
TEMP (if KEYOPT(1) = 2)
UX, UY, UZ, TEMP, VOLT (if KEYOPT(1) = 3)
TEMP, VOLT (if KEYOPT(1) = 4)
UX, UY, UZ, VOLT (if KEYOPT(1) = 5)
VOLT (if KEYOPT(1) = 6)
MAG (if KEYOPT(1) = 7)
Real Constants
R1, R2, FKN, FTOLN, ICONT, PINB,
PMAX, PMIN, TAUMAX, CNOF, FKOP, FKT,
COHE, TCC, FHTG, SBCT, RDVF, FWGT,
ECC, FHEG, FACT, DC, SLTO, TNOP,
TOLS, MCC
See Table 173.1: "CONTA173 Real Constants" for descriptions of the real constants.
Material Properties
DAMP, MU, EMIS (MP command)
FRIC (TB command; see Contact Friction)
CZM (TB command; see Cohesive Zone Materials Used for Debonding in the Contact Technology Guide)
Surface Loads
Convection, Face 1 (I-J-K-L)
Heat Flux, Face 1 (I-J-K-L)
Special Features
Nonlinear
Large deflection
Isotropic or orthotropic friction
Debonding
Birth and death
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:

0 -- 

UX, UY, UZ

1 -- 

UX, UY, UZ, TEMP

2 -- 

TEMP

3 -- 

UX, UY, UZ, TEMP, VOLT

4 -- 

TEMP, VOLT

5 -- 

UX, UY, UZ, VOLT

6 -- 

VOLT

7 -- 

MAG

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

KEYOPT(4)

Location of contact detection point:

0 -- 

On Gauss point (for general cases)

1 -- 

On nodal point - normal from contact surface

2 -- 

On nodal point - normal to target surface

Note

When using the multipoint constraint (MPC) approach to define surface-based constraints, use KEYOPT(4) in the following way: set KEYOPT(4) = 1 for a force-distributed constraint, set KEYOPT(4) = 2 for a rigid surface constraint. See Surface-based Constraints for more information.

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)

Contact stiffness variation (used to enhance stiffness updating when KEYOPT(10) > 0):

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

KEYOPT(7)

Element level time incrementation control:

0 -- 

No control

1 -- 

Automatic bisection of increment

2 -- 

Change in contact predictions 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

Note

For KEYOPT(7) = 2 or 3, includes automatic bisection of increment. Activated only if SOLCONTROL,ON,ON at the procedure level.

KEYOPT(8)

Asymmetric contact selection:

0 -- 

No action

2 -- 

ANSYS internally selects which asymmetric contact pair is used at the solution stage (used only when symmetry contact is defined).

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

Note

For KEYOPT(9) = 1, 3, or 4, the indicated initial gap effect is considered only if KEYOPT(12) = 4 or 5.

KEYOPT(10)

Contact stiffness update:

0 -- 

Each load step if FKN is redefined during load step (pair based).

1 -- 

Each substep based on mean stress of underlying elements from the previous substep (pair based).

2 -- 

Each iteration based on current mean stress of underlying elements (pair based).

3 -- 

Each load step if FKN is redefined during load step (individual element based).

4 -- 

Each substep based on mean stress of underlying elements from the previous substep (individual element based).

5 -- 

Each iteration based on current mean stress of underlying elements (individual element based).

Note

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).

KEYOPT(11)

Shell thickness effect:

0 -- 

Exclude

1 -- 

Include

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)

Note

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 173.1  CONTA173 Real Constants

No.NameDescriptionFor more information, see this section in the Contact Technology Guide . . .
1R1

Target radius for cylinder, cone, or sphere

Defining the Target Surface

2R2

Target radius at second node of cone

Defining the Target Surface

3FKN

Normal penalty stiffness factor

Determining Contact Stiffness and Penetration

4FTOLN

Penetration tolerance factor

Determining Contact Stiffness and Penetration

5ICONT

Initial contact closure

Adjusting Initial Contact Conditions

6PINB

Pinball region

Determining Contact Status and the Pinball Region

or

Defining Influence Range (PINB)

7PMAX

Upper limit of initial allowable penetration

Adjusting Initial Contact Conditions

8PMIN

Lower limit of initial allowable penetration

Adjusting Initial Contact Conditions

9TAUMAX

Maximum friction stress

Choosing a Friction Model

10CNOF

Contact surface offset

Adjusting Initial Contact Conditions

11FKOP

Contact opening stiffness or contact damping

Selecting Surface Interaction Models

12FKT

Tangent penalty stiffness factor

Determining Contact Stiffness

13COHE

Contact cohesion

Choosing a Friction Model

14TCC

Thermal contact conductance

Modeling Conduction

15FHTG

Frictional heating factor

Modeling Heat Generation Due to Friction

16SBCT

Stefan-Boltzmann constant

Modeling Radiation

17RDVF

Radiation view factor

Modeling Radiation

18FWGT

Heat distribution weighing factor

Modeling Heat Generation Due to Friction (thermal)

or

Heat Generation Due to Electric Current (electric)

19ECC

Electric contact conductance

Modeling Surface Interaction

20FHEG

Joule dissipation weight factor

Heat Generation Due to Electric Current

21FACT

Static/dynamic ratio

Static and Dynamic Friction Coefficients

22DC

Exponential decay coefficient

Static and Dynamic Friction Coefficients

23SLTO

Allowable elastic slip

Using FKT and SLTO

24TNOPMaximum allowable tensile contact pressure

Chattering Control Parameters

25TOLS

Target edge extension factor

Selecting Location of Contact Detection

26MCC

Magnetic contact permeance

Modeling Magnetic Contact

CONTA173 Output Data

The solution output associated with the element is in two forms:

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 173.2: "CONTA173 Element Output Definitions" gives element output at the element level. In the results file, the nodal results are obtained from its closest integration point.

Table 173.2  CONTA173 Element Output Definitions

NameDefinitionOR
ELElement NumberYY
NODESNodes I, J, K, LYY
XC, YC, ZCLocation where results are reportedY5
TEMPTemperatures T(I), T(J), T(K), T(L)YY
VOLUAREAYY
NPINumber of integration pointsY-
ITRGETTarget surface number (assigned by ANSYS)Y-
ISOLIDUnderlying solid or shell element numberY-
CONT:STATCurrent contact statuses11
OLDSTOld contact statuses11
ISEGUnderlying current target numberYY
OLDSEGUnderlying old target numberY-
CONT:PENECurrent penetration (gap = 0; penetration = positive value)YY
CONT:GAPCurrent gap (gap = negative value; penetration = 0)YY
NGAPNew or current gap at current converged substep (gap = negative value; penetration = positive value)Y-
OGAPOld gap at previously converged substep (gap = negative value; penetration = positive value)Y-
IGAPInitial gap at start of current substep (gap = negative value; penetration = positive value)YY
CONT:PRESNormal contact pressureYY
TAUR/TAUS[7]Tangential contact stressesYY
KNCurrent normal contact stiffness (Force/Length3)YY
KTCurrent tangent contact stiffness (Force/Length3)YY
MU[8]Friction coefficientY-
TASS/TASR[7]Total (algebraic sum) sliding in S and R directions33
AASS/AASR[7]Total (absolute sum) sliding in S and R directions33
TOLNPenetration toleranceYY
CONT:SFRICFrictional stress SQRT (TAUR**2+TAUS**2)YY
CONT:STOTALTotal stress SQRT (PRES**2+TAUR**2+TAUS**2)YY
CONT:SLIDETotal sliding SQRT (TASS**2+TASR**2)YY
DBAPenetration variationYY
PINBPinball Region-Y
CNFX[4]Contact element force-X component-Y
CNFYContact element force-Y component-Y
CNFZContact element force-Z component-Y
CONVConvection coefficientYY
RACRadiation coefficientYY
TCCConductance coefficientYY
TEMPSTemperature at contact pointYY
TEMPTTemperature at target surfaceYY
FXCVHeat flux due to convectionYY
FXRDHeat flux due to radiationYY
FXCDHeat flux due to conductanceYY
FDDISFrictional energy dissipation66
FLUXTotal heat flux at contact surfaceYY
FXNPFlux input-Y
CNFHContact element heat flow-Y
CAREAContacting area-Y
JCONTContact current density (Current/Unit Area)YY
CCONTContact charge density (Charge/Unit Area)YY
HJOUContact power/areaYY
ECURTCurrent per contact element-Y
ECHARCharge per contact element-Y
ECCElectric contact conductance (for electric current DOF), or electric contact capacitance per unit area (for piezoelectric or electrostatic DOFs)YY
VOLTSVoltage on contact nodesYY
VOLTTVoltage on associated targetYY
CNOSTotal number of contact status changes during substepYY
TNOPMaximum allowable tensile contact pressureYY
SLTOAllowable elastic slipYY
MCCMagnetic contact permeanceYY
MFLUXMagnetic flux densityYY
MAGSMagnetic potential on contact nodeYY
MAGTMagnetic potential on associated targetYY
ELSIElastic slip distance for sticking contact within a substep-Y
DTSTARTLoad step time during debondingYY
DPARAMDebonding parameterYY
DENERIEnergy released due to separation in normal direction - mode I debondingYY
DENERIIEnergy released due to separation in tangential direction - mode II debondingYY
  1. 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
  2. ANSYS will evaluate model to detect initial conditions.

  3. Only accumulates the sliding when contact occurs.

  4. Contact element forces are defined in the global Cartesian system

  5. Available only at centroid as a *GET item.

  6. FDDIS = (contact friction stress)*(sliding distance of substep)/(time increment of substep)

  7. For the case of orthotropic friction, components are defined in the global Cartesian system (default) or in the local element coordinate system specified by ESYS.

  8. For orthotropic friction, an equivalent coefficient of friction is output.

Note

If ETABLE is used for the CONT items, the reported data is averaged across the element.

Table 173.3: "CONTA173 Item and Sequence Numbers" 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. The following notation is used in Table 173.3: "CONTA173 Item and Sequence Numbers":

Name

output quantity as defined in the Table 173.2: "CONTA173 Element Output Definitions"

Item

predetermined Item label for ETABLE command

E

sequence number for single-valued or constant element data

I,J,K,L

sequence number for data at nodes I,J,K,L

Table 173.3  CONTA173 Item and Sequence Numbers

Output Quantity NameETABLE and ESOL Command Input
ItemEIJKL
PRESSMISC131234
TAURSMISC-5678
TAUSSMISC-9101112
FLUXSMISC-14151617
FDDISSMISC-18192021
FXCVSMISC 22232425
FXRDSMISC-26272829
FXCDSMISC-30313233
FXNPSMISC-34353637
JCONTSMISC-38394041
CCONTSMISC-38394041
HJOUSMISC-42434445
MFLUXSMISC-46474849
STAT[1]NMISC411234
OLDSTNMISC-5678
PENE[2]NMISC-9101112
DBANMISC-13141516
TASRNMISC-17181920
TASSNMISC-21222324
KNNMISC-25262728
KTNMISC-29303132
TOLNNMISC-33343536
IGAPNMISC-37383940
PINBNMISC42----
CNFXNMISC43----
CNFYNMISC44----
CNFZNMISC45----
ISEGNMISC-46474849
AASRNMISC-50515253
AASSNMISC-54555657
CAREANMISC-58596061
MUNMISC-62636465
DTSTARTNMISC-66676869
DPARAMNMISC-70717273
TEMPSNMISC-78798081
TEMPTNMISC-82838485
CONVNMISC-86878889
RACNMISC-90919293
TCCNMISC-94959697
CNFHNMISC98----
ECURTNMISC99----
ECHARNMISC99----
ECCNMISC-100101102103
VOLTSNMISC-104105106107
VOLTTNMISC-108109110111
CNOSNMISC-112113114115
TNOPNMISC-116117118119
SLTONMISC-120121122123
MCCNMISC-124125126127
MAGSNMISC-128129130131
MAGTNMISC-132133134135
ELSINMISC-136137138139
DENERINMISC-140141142143
DENERIINMISC-144145146147
  1. Element Status = highest value of status of integration points within the element

  2. Penetration = positive value, gap = negative value

You can display or list contact results through several POST1 postprocessor commands. The contact specific items for the PLNSOL, PLESOL, PRNSOL, and PRESOL commands are listed below:

STATContact status
PENEContact penetration
PRESContact pressure
SFRICContact friction stress
STOTContact total stress (pressure plus friction)
SLIDEContact sliding distance
GAPContact gap distance
FLUXTotal heat flux at contact surface
CNOSTotal number of contact status changes during substep

CONTA173 Assumptions and Restrictions

  • The 3-D contact element must coincide with the external surface of the underlying solid or shell element or with the original elements comprising the superelement.

  • 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 Lagrangian multiplier 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 solid, shell, beam, or target elements.

CONTA173 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 Professional. 

  • The MU material property is not allowed

  • The birth and death special feature is not allowed.

  • The DAMP material property is not allowed.

ANSYS Structural. 

  • The VOLT DOF (KEYOPT(1) = 3 through 6) is not allowed.

  • The MAG DOF (KEYOPT(1) = 7) is not allowed.

ANSYS Mechanical. 

  • The MAG DOF (KEYOPT(1) = 7) is not allowed.

Your Ad Here