User Tools

Site Tools


elements:cntcp

CNTCP

2D contact element

Description

The element consists of 2 or 3 nodes (linear or parabolic) and is related to a foundation or a tool previously defined.

Convention: when the element is considered in increasing order of nodes (1 → 2 → 3), the first body is to the left, and the tool or foundation to the right.

This element can be used with law COU2DC as it includes friction and thermal exchange by radiation-convection or with thermal resistance.
The element can also be used with classic laws; it is compatible with plane stress state, plane strain state, and generalized plane state, as well as axisymmetric state. It can be used for mechanical, thermal, or thermo-mechanical analysis.

Type: 211

Implemented by: R. Charlier (1989), Zhu (1991), A-M. Habraken (1991)

Files

Prepro: CNTCPA.F
Lagamine: CNTCPB.F

Input file

Title (A5)
TITLE“CNTCP” in the first 5 columns
Control data (3I5)
NELEMNumber of elements
INSIG= 0 no initial stresses
= 1 initial stresses (see below)
IDOMANumber of the ES group to which the elements are connected (necessary in case of remeshing)
Initial stresses - only if INSIG = 1
If law = COU2DC (4G10.0)
The pressure varies as: $PRESS = PRESS0 + (Y*DPRES) \text{ for } P\geq 0$
Tau varies as: $TAU = TAU0+(Y*DTAU)$
PRES0Pressure of the contact at the axis origin
DPRESCoefficient of variation of the pressure along Y (= 0 → constant pressure)
TAU0TAU at the axis origin
DTAUCoefficient of variation of TAU along Y (= 0 → constant TAU)
If law ≠ COU2DC (7G10.0)
All stresses vary as: $SIGnn=SI0nn+(Y*DSIG)$
SI011Stresses at the axis origin (y=0)
expressed in local axes
SI022
SI012
SI033
DSIGCoefficient of stress variation along Y ( = 0 → constant stress)
DTAU
RIGMPenalty
= 1 by default
= 1/e inverse of the element thickness
Definition of the elements (5I5/3I5)
NINTENumber of integration points (1 to 10, 5 max can be drawn by DESFIN - recommanded value: 2)
LMATELaw number
IFONDNumber of the foundation or tool
If the foundation number is equal to 0, the boundary thermal flow is calculated with the ambient, without mechanical contact with any foundation.
INTYPType of numerical integration (recommanded value: 0)
0 → Gauss
1 → Lobatto
2 → Newton-Cote
IRIGFType of contact
0 → rigid foundation or tool
1 → uncoupled solid/solid contact
One contact element on each structure, the interpenetration distance is divided by 2.
Suitable for solids with similar stiffnesses.
2 → coupled solid/solid contact
Only one contact element must be defined on a solid, the other being its foundation. The computation of MBAND and NHICO must be actualized (see ISTRA(4)).
Suitable for solids of different stiffnesses, with one (the foundation) can be more roughly approximated.
3 → coupled solid/solid contact
One contact element on each structure, the force is divided by 2. None of the structure is privileged. The computation of MBAND and NHICO must be actualized (see ISTRA(4)).
Suitable for solids of different stiffnesses, both must be properly represented.
4 → Coupled solid/foundation piloted in displacement contact (only if pilot node)
Only one contact element defined on the structure.
Corresponds to the case where IRIGF = 2, but with a reduction of the size of LM and AK and without any necessity to repeat the calculation of MBAND and NHICO (see ISTRA(4))
The computation method of AK with perturbation is adapted to the notion of piloted foundation
5 → Coupled solid/foundation piloted in rotation contact (only if pilot node)
Only one contact element defined on the structure.
Corresponds to the case where IRIGF = 2, but with a reduction of the size of LM and AK and without any necessity to repeat the calculation of MBAND and NHICO (see ISTRA(4))
The computation method of AK with perturbation is adapted to the notion of piloted foundation
NODES(3)List of nodes (2 or 3)

Results

If the law is COU2DC:
Stresses: pressure, tangent stress, thermal flux
Internal variables: The first values are the ones relating to the law (for Coulomb law, this will be: plasticity indicator, dissipation, information on mechanical dissipation in thermal flow). The 5 last values correspond to the contact geometry. They are:

  • The segment number of the foundation in contact
  • The interpenetration distance
  • The jacobian
  • NOCO the contact indicator given by CALFON
  • The relative tangent rate


If the law is not COU2DC:
Stresses: $\sigma_{11}$, $\sigma_{22}$, $\sigma_{12}$, $\sigma_{33}$ in local axes

elements/cntcp.txt · Last modified: 2020/08/25 15:46 (external edit)