8.11. Assigning Contact Properties via Real Constants

The program uses real constants to control contact parameters and properties for both general contact and pair-based contact. However, the assignment of these properties is handled differently for the two types of contact.

In a pair-based definition, the real constant ID is assigned to the contact and target elements as one of the element attributes. In a general contact definition, the real constant ID (REALID) is assigned via the GCDEF command which controls the contact properties at a specific general contact interface identified by section IDs (SECT1 and SECT2). Unlike pair-based contact, the real constant ID (REALID) used on the GCDEF command is not a contact element attribute.

If the REALID is not specified or is set to zero on GCDEF, all default real constant values are used for the specified interface throughout the entire analysis.

Most real constants used in a pair-based contact definition (see Real Constants) are also available in a general contact definition, with the following exceptions:

  • R1, R2, PPCN, FPAT, COR, FDMD, FDMS, BSRL, KSYM, and TFOR are not supported in a general contact definition.

  • General contact provides default values for TCC, ECC, MCC, PCC, and DCC based on the highest conductivity and element sizes of underlying elements from both sides of the contact interface. The default values of these real constants are not provided in the pair-based definition.

  • The program-provided defaults for certain real constants may not be identical for pair-based contact and general contact. For example, the resulting contact stiffness (FKN) and pinball radius (PINB) vary between the two contact definitions. The resulting pinball (PINB) in the general contact definition is no longer important for standard contact. However, it is critical in some situations such as for bonded contact (TB,INTER,,,ABOND) and for resolving initial interference fit (TB,INTER with TBDATA,,2). Specifying an absolute PINB value (negative input value) is highly recommended in these situations.

The table below compares default real constant values between pair-based contact and general contact and also indicates which real constants are not used in a general contact definition.

Table 8.2: Real Constant Defaults for Pair-Based vs. General Contact

Real ConstantsDescriptionPair-Based ContactGeneral Contact
No.Name
1R1Radius associated with target geometry

0

Not supported [5]

Target radius (CONTA177)

Calculated by program

2R2

Radius associated with target geometry

0

Not supported [5]

Superelement thickness

1

Contact radius (CONTA177)

Calculated by program

3FKNNormal penalty stiffness factor11 [1]
4FTOLNPenetration tolerance factor0.10.1 [1]
5ICONTInitial contact closure00
6PINBPinball region[2][2] [1]
7PZERPressure at zero penetration[6][6]
8CZERInitial contact clearance factor0.010.01
9TAUMAXMaximum friction stress1.00E+201.00E+20
10CNOFContact surface offset00
11FKOPContact opening stiffness11 [1]
12FKTTangent penalty stiffness factor11 [1]
13COHEContact cohesion00
14TCCThermal contact conductance0Auto [3]
15FHTGFrictional heating factor11
16SBCTStefan-Boltzmann constant00
17RDVFRadiation view factor11
18FWGTHeat distribution weighing factor0.50.5
19ECCElectric contact conductance0Auto [3]
20FHEGJoule dissipation weighting factor11
21FACTStatic/dynamic ratio11
22DCExponential decay coefficient00
23SLTOAllowable elastic slip1%1% [1]
24TNOPMaximum allowable tensile contact pressure[4][4]
25TOLSTarget edge extension factor22
26MCCMagnetic contact permeance0Auto [3]
27PPCNPressure-penetration criterion0Not supported
28FPATFluid penetration acting time0.01Not supported
29CORCoefficient of restitution1Not supported
30STRMLoad step number for ramping penetration1 [7]1 [7]
31FDMNNormal stabilization damping factor11 [1]
32FDMTTangential stabilization damping factor0.0010.001 [1]
33FDMDDestabilizing squeal damping factor1Not supported
34FDMSStabilizing squeal damping factor0Not supported
35TBNDCritical bonding temperatureNo defaultNo default
36WBIDInternal contact pair ID (used only by the Results Tracker in Workbench Mechanical)n/an/a
37PCCPore fluid contact permeability coefficient0Auto [3]
38PSEEPore fluid seepage coefficient00
39ABPPAmbient pore pressure 00
40FPFTGap pore fluid flow participation factor 00
41FPWTGap pore fluid flow distribution weighting factor0.50.5
42DCCContact diffusivity coefficient0Auto [3]
43DCONDiffusive convection coefficient00
44ABDCAmbient concentration00
45BSRLReal ID of base contact pair (after contact pair-splitting)0Not supported
46KSYMReal ID of the companion pair for symmetric contact or self contact (after contact pair-splitting)0Not supported
47TFORPair-based force convergence toleranceNo defaultNot supported
48TENDEnding time for ramping contact stiffness00

  1. The internal default may be different for pair-based contact and general contact.

  2. The program-provided defaults for pinball region varies between the two contact definitions.

  3. The default value is calculated as a function of highest conductivity and mesh size at the interface.

  4. TNOP defaults to the force convergence tolerance divided by contact area at contact nodes.

  5. The R1 and R2 real constants are not used in a general contact definition. However, you can input radius values required by contact and target elements via a section definition (SECTYPE and SECDATA commands). This includes the radius of a general contact surface that represents 3D beam-to-beam or 3D edge-to-edge contact, and radius values associated with rigid target segments. See Applying Surface Geometry Correction for details.

  6. Program-calculated value.

  7. If TEND is specified as a positive value, STRM defaults to zero.

Examples

When you assign a non-zero REALID on the GCDEF command, you will typically input the same number in the real constant set ID field (NSET) on the R and RMODIF commands.

The following example defines a contact stiffness factor (FKN) of 10 at a specific contact interface. REALID = 14 is used to identify the non-default property:

GCDEF,,SECT1,SECT2,,14
R,14,,,10

If the REALID has been specified by a previous GCDEF command, you can use the RMODIF command to add any non-default real constant values regardless of whether or not an associated R command was issued.

Use NSET = GCN on the RMODIF command to define or modify real constants for all real constant sets that were previously assigned by GCDEF. For example:

RMODIF,GCN,STLOC,V1,V2,V3,V4,V5,V6

Note that the GCN label only acts on real constant sets predefined by GCDEF commands. The default real constant set (REALID = 0 or undefined) is not affected.

For example, the following command defines a contact stiffness factor (FKN) of 10 for all real constant sets previously specified by GCDEF commands:

RMODIF,GCN,3,10