User Tools

Site Tools


elements:faif2

FAIF2D

Interface 2D fault

Description

Type: 210

The element consists of 2, 3, 4, or 5 nodes (line, parabola - DoF 1, 2, 3, 4, 5) relating to one side of the fault (= the structure), and 2, 3, 4, or 5 additional nodes to discretize the flow of water, gas, or temperature inside the fault (DoF 3, 4, 5). It also relates to a predefined foundation.
This element is dedicated to soil mechanics and in particular to structure-structure contact (fault). It includes pressure, friction, flow axial to the fault (water, gas, temperature), flow (water, gas, temperature) transversal to the fault (along the thickness).

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 mechanical law INTME2 or a classic mechanical law and flow law INTFL2.

The element is compatible with plane stress state, plane strain state, and generalized plane state, as well as axisymmetric state. It can be used for thermo-mechanical analysis with 5 DoFs.

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

Implemented by: J.P. Radu, 2007 (taken from FAIN2)

Files

Prepro: FAIF2A.F
Lagamine: FAIF2B.F

Input file

Title (A5)
TITLE“FAIF2” in the first 5 columns
Control data (3I5)
NELEMNumber of elements
INSIG= 0 no initial stresses
= 1 or 2 → initial stresses
IDOMANumber of the ES group to which the elements are connected (necessary in case of remeshing)
Initial stresses - only if INSIG = 1
If law = INTME2/INTME3 (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 ≠ INTME2/INTME3 (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 SIGnn)
DTAU
RIGMPenalty
= 1.0 by default
= 1/e inverse of the element thickness
Definition of the elements (6I5/10I5)
NINTENumber of integration points (1 to 10, 5 max can be drawn by DESFIN - recommanded value: 3)
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.
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.
NODES(10)List of nodes: 1, 2, (3), (4), (5), 1', 2', (3'), (4'), (5')
Nodes 3 and 3' can be 0; in that case the element is a straight line
Nodes 4 and 4' can be 0; in that case the element is a parabola
Nodes 5 and 5' can be 0; in that case the element is a cubic

Results

If the law is INTME2/INTME3:
Stresses (20):

  • Pressure
  • Tangent stress
  • 0
  • 0
  • Longitudinal water flow,
  • Stored water flow,
  • “FI” (between Foundation and Inside) transversal water flow
  • “IS” (between Inside and Structure) transversal water flow
  • Longitudinal gas flow,
  • Stored water flow,
  • “FI” (between Foundation and Inside) transversal gas flow
  • “IS” (between Inside and Structure) transversal gas flow
  • Longitudinal thermal flow,
  • Stored thermal flow,
  • “FI” (between Foundation and Inside) transversal thermal flow
  • “IS” (between Inside and Structure) transversal thermal flow
  • Longitudinal vapor flow,
  • Stored vapor flow,
  • “FI” (between Foundation and Inside) transversal vapor flow
  • “IS” (between Inside and Structure) transversal vapor flow

If the law is not INTME2/INTME3:
Stresses:

  • $\sigma_{11}$, $\sigma_{22}$, $\sigma_{12}$, $\sigma_{33}$ in local axes
  • Longitudinal water flow, Stored water flow, “FI” (between Foundation and Inside) transversal water flow, “IS” (between Inside and Structure) transversal water flow
  • Longitudinal gas flow, Stored water flow, “FI” (between Foundation and Inside) transversal gas flow, “IS” (between Inside and Structure) transversal gas flow
  • Longitudinal thermal flow, Stored thermal flow, “FI” (between Foundation and Inside) transversal thermal flow, “IS” (between Inside and Structure) transversal thermal flow
  • Longitudinal vapor flow, Stored vapor flow, “FI” (between Foundation and Inside) transversal vapor flow, “IS” (between Inside and Structure) transversal vapor flow

Internal variables:
The first values are the ones relating to the law
For INTME2/INTME3, this will be: plasticity indicator, dissipation, fault closing $V$, fault opening $D$).
The 6 next 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
  • The relative sliding

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

  • The relative permeability to water
  • The relative permeability to gas
  • The porosity
  • The saturation to water
  • The suction
  • The specific mass of water
  • The specific mass of gas
  • Peclet number
  • The water content
  • The specific mass of vapor
  • The vapor pressure
  • The relative humidity
  • The mass of liquid water per unit volume
  • The mass of dry air per unit volume of soil
  • The mass of vapor per unit volume of soil
  • The intrinsic longitudinal permeability
  • The gas saturation
  • The partial pressure of hydrogen (not used)
  • The water pressure inside the fault
  • The gas pressure inside the fault
  • The temperature inside the fault
  • The transversal transmissivity of water
  • The transversal transmissivity of gas
  • The thermal transversal transmissivity
  • The thermal conductivity
elements/faif2.txt · Last modified: 2020/08/25 15:46 (external edit)