2D contact element
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)
Prepro: CNTCPA.F
Lagamine: CNTCPB.F
Title (A5) | |
---|---|
TITLE | “CNTCP” in the first 5 columns |
Control data (3I5) | |
NELEM | Number of elements |
INSIG | = 0 no initial stresses |
= 1 initial stresses (see below) | |
IDOMA | Number 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)$ |
|
PRES0 | Pressure of the contact at the axis origin |
DPRES | Coefficient of variation of the pressure along Y (= 0 → constant pressure) |
TAU0 | TAU at the axis origin |
DTAU | Coefficient of variation of TAU along Y (= 0 → constant TAU) |
If law ≠ COU2DC (7G10.0) | |
All stresses vary as: $SIGnn=SI0nn+(Y*DSIG)$ | |
SI011 | Stresses at the axis origin (y=0) expressed in local axes |
SI022 | |
SI012 | |
SI033 | |
DSIG | Coefficient of stress variation along Y ( = 0 → constant stress) |
DTAU | |
RIGM | Penalty = 1 by default = 1/e inverse of the element thickness |
Definition of the elements (5I5/3I5) | |
NINTE | Number of integration points (1 to 10, 5 max can be drawn by DESFIN - recommanded value: 2) |
LMATE | Law number |
IFOND | Number 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. |
INTYP | Type of numerical integration (recommanded value: 0) 0 → Gauss 1 → Lobatto 2 → Newton-Cote |
IRIGF | Type 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) |
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:
If the law is not COU2DC:
Stresses: $\sigma_{11}$, $\sigma_{22}$, $\sigma_{12}$, $\sigma_{33}$ in local axes