Table of Contents

FAIF3

3D fault interface

Description

Type: 227

The element consists of 4 or 8 nodes (DoF 1, 2, 3, 4, 5, and 6) relating to one side of the fault (= the structure), and 4 or 8 additional nodes to discretize the flow of water, gas, or temperature inside the fault (DoF 4, 5, and 6). It also relates to a predefined foundation.
The normal (red side in DESFIN) must go out of the structure.

The foundation, defined by a series of triangular or quadrangular sides, is defined in the “FOUND” section of the *.lag file. If the foundation is defined by a series of sides, the normal is given by the product of vectors $\vec{12}$ and $\vec{13}$ built on the nodes of the foundation $\vec{n}=\vec{12} \wedge\vec{13}$ of each side. This normal must be directed inwards the foundation.

This element is dedicated to soil mechanics and in particular to structure-structure contact (fault). It includes pressure, friction, flows axial to the fault (water, gas, thermal), and two (x3) flows (water, gas, thermal) transversal to the fault (along the thickness).

The element can work in coupled mechanical-water-air-temperature analysis, with 6 DoF.

This element can use mechanical law INTME3 or a classic law and flow law INTFL3.

Warning: When using renumbering of the DoFs, only directional renumbering works.

Implemented by: J.P. Radu, 2007

Files

Prepro: FAIF3A.F
Lagamine: FAIF3B.F

Input file

Title (A5)
TITLE“FAIF3” in the first 5 columns
Control data (3I5)
NELEMNumber of elements
INSIG= 0 no initial stresses
= 1 → initial stresses
Initial stresses - only if INSIG = 1 (6G10.0)
The pressure varies as: $PRESS = PRES0 + (Z*DPRES)$
Tau varies as: $TAUn = TAUn0+(Z*DTAUn)$
PRES0Pressure of the contact at the axis origin
DPRESCoefficient of variation of the pressure along Y (= 0 → constant pressure)
TAU10TAU along direction 1 at the axis origin
DTAU1Coefficient of variation of TAU1 along Z (= 0 → constant TAU1)
TAU20TAU along direction 2 at the axis origin
DTAU2Coefficient of variation of TAU2 along Z (= 0 → constant TAU2)
Definition of the elements (5I5/8I5/8I5)
NINTENumber of integration points (1, 4, or 9)
LMATMMechanical law number
LMATFFlow law
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.
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.
NODES(1 to 8)List of nodes on the structure side. Nodes 5 to 8 can be 0; in that case the element is of linear degree
NODES(9 to 16)List of nodes inside the fault. Nodes 5' to 8' can be 0; in that case the element is of linear degree

Results

Stresses (24):

Internal variables (4+2+25 = 31):
The first values are the ones relating to the law
For INTME2/INTME3, this will be:

The 2 next values correspond to the contact geometry. They are:

The last values are those relating to the flow law; for INTFL3, they are: