User Tools

Site Tools


elements:ehmic

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
elements:ehmic [2023/11/23 11:01]
gilles
elements:ehmic [2023/12/12 16:02] (current)
gilles [Description]
Line 15: Line 15:
   * [[laws:​hmic|HMIC]]:​ Hydraulic micro law with WA coupling and implicit mechanics   * [[laws:​hmic|HMIC]]:​ Hydraulic micro law with WA coupling and implicit mechanics
 \\  \\ 
-Implmented ​by: F. Bertrand & G. CORMAN, ​2019+Implemented ​by: F. Bertrand & G. Corman (2019
 +\\ \\ 
 +The framework definition of this element can be found in Corman (2024)((Corman,​ G. (2024). Hydro-mechanical modelling of gas transport processes in clay host rocks in the context of a nuclear waste repository. PhD thesis, University of Liège. https://​hdl.handle.net/​2268/​307996)).
 ==== Files ==== ==== Files ====
 Prepro: EHMICA.F ​ \\ Prepro: EHMICA.F ​ \\
Line 25: Line 27:
 ^Control data (4I5)^^ ^Control data (4I5)^^
 |NELEM|Number of elements| |NELEM|Number of elements|
-|ISPSMAS|= 0 → Nothing \\ = 1 → Take into account the specific mass if and only if NTANA<0| 
 |INSIG|= 0 → No initial stress \\ = 1 or 2 → Initial stresses| |INSIG|= 0 → No initial stress \\ = 1 or 2 → Initial stresses|
-|INBIO|= 0 → No Biot coefficient \\ = 1 → Isotropic Biot coefficient | 
-^Specific mass in dynamic analysis - Only if ISPMAS = 1 (1G10.0)^^ 
-|SPEMAS|Specific mass| 
 ^Initial stresses - Only if INSIG > 0 (4G10.0)^^ ^Initial stresses - Only if INSIG > 0 (4G10.0)^^
 |If INSIG=1: $\sigma_y=\sigma_{y0}+yd\sigma_{y}$ \\ If INSIG=2: $\sigma_y=min(\sigma_{y0}+yd\sigma_y,​0)$|| |If INSIG=1: $\sigma_y=\sigma_{y0}+yd\sigma_{y}$ \\ If INSIG=2: $\sigma_y=min(\sigma_{y0}+yd\sigma_y,​0)$||
Line 37: Line 35:
 |AK0Z|$k_0$ ratio $\sigma_z/​\sigma_y$ (if AK0Z=0, AK0Z=AK0X)| |AK0Z|$k_0$ ratio $\sigma_z/​\sigma_y$ (if AK0Z=0, AK0Z=AK0X)|
 |The computation of SIGY0 and DSIGY must take into account the apparent specific mass, defined as \[\rho_a'​=[(1-n)\rho_s+nS_w\rho_w]-\rho_w\] where: \\ $\rho_s$ is the solid specific mass - this represents the specific mass of a fictive sample where ther is no porosity, i.e. where the grains occupy the whole volume of the sample \\  $\rho_w$ is the fluid specific mass \\ $n$ is the porosity defined in the flow law related to the element \\ $S_w$ fluid saturation, ∈ [0,1]|| |The computation of SIGY0 and DSIGY must take into account the apparent specific mass, defined as \[\rho_a'​=[(1-n)\rho_s+nS_w\rho_w]-\rho_w\] where: \\ $\rho_s$ is the solid specific mass - this represents the specific mass of a fictive sample where ther is no porosity, i.e. where the grains occupy the whole volume of the sample \\  $\rho_w$ is the fluid specific mass \\ $n$ is the porosity defined in the flow law related to the element \\ $S_w$ fluid saturation, ∈ [0,1]||
-^Biot coefficient - Only if INBIO = 1 (1G10.0)^^ 
-|CBIOT|Biot coefficient| 
 ^Definition of the elements (6I5/​16I5(/​9I5)) ^^ ^Definition of the elements (6I5/​16I5(/​9I5)) ^^
-|NNODM|Number of nodes for the mechancial description: ​3, 4, 6, 8, 15, or 25+|NNODM|Number of nodes for the mechancial description:​ 8| 
-|NINTM|Number of integration point (1, 3, 4, 7, 9, 12, or 16) for the mechanical description| +|NINTM|Number of integration point for the mechanical description: 4
-|LMATM|Mechanical ​material+|LMATM|Mechanical ​law
-|NNODP|Number of nodes for the flow description: ​3, 4, 6, 8, 15, or 25+|NNODP|Number of nodes for the flow description:​ 8| 
-|NINTP|Number of integration points ​(1, 3, 4, 7, 9, 12, or 16) for the flow description \\ Must be equal to NINTM| +|NINTP|Number of integration points for the flow description: 4 \\ Must be equal to NINTM| 
-|LMATF|Flow ​material|+|LMATF|Flow ​law|
 |NODES(NNODEM)|List of nodes| |NODES(NNODEM)|List of nodes|
 +^Definition of the hydraulic micro-elements in the microstructure (1I5/4I5) ^^ 
 +|NUMEL2|Number of hydraulic micro-elements| 
 +|IELEM2|No. of the hydraulic micro-element| 
 +|ILAW|Type of element: 1=bedding plane, 2=bundle of tubes, 3=bridging plane| 
 +|NDUN|No. of the micro-node at one side of the micro-element| 
 +|NDDEUX|No. of the micro-node at the other side of the micro-element| 
 +^Definition of the hydraulic micro-nodes in the microstructure (1I5/10I5) ^^ 
 +|NUMNDH|Number of hydraulic micro-nodes| 
 +|INOD2|No. of the hydraulic micro-node| 
 +|IEDGE|No. of the microstructure boundary to which belong the micro-node: 0=none, 1=left or bottom, 2=right or top| 
 +|IELUN|No. of the 1st micro-element connected to the micro-node| 
 +|IELDEUX|No. of the 2nd micro-element connected to the micro-node| 
 +|IELTROIS|No. of the 3rd micro-element connected to the micro-node| 
 +|IELQUTR|No. of the 4th micro-element connected to the micro-node| 
 +|IELCINQ|No. of the 5th micro-element connected to the micro-node| 
 +|IELSIX|No. of the 6th micro-element connected to the micro-node| 
 +|IELSEPT|No. of the 7th micro-element connected to the micro-node| 
 +|IELHUIT|No. of the 8th micro-element connected to the micro-node|
 ===== Results ===== ===== Results =====
   * Stresses (in global axes)   * Stresses (in global axes)
     * Mechanical stresses $\sigma_x$, $\sigma_y$, $\sigma_{xy}$,​ $\sigma_z$ ​     * Mechanical stresses $\sigma_x$, $\sigma_y$, $\sigma_{xy}$,​ $\sigma_z$ ​
-    * Flow in water $f_{wx}$, $f_{wy}$, $f_{w,​stored}$, 0 +    * Flow in water $f_{wx}$, $f_{wy}$, $f_{w,​stored}$ 
-    * Flow in air $f_{ax}$, $f_{ay}$, $f_{a,​stored}$, 0 +    * Flow in air $f_{ax}$, $f_{ay}$, $f_{a,​stored}$ 
-    * Thermal flow $f_{tx}$, $f_{ty}$$f_{t,stored}$, 0+    * Advection dissolved gas flux $f_{adx}$, $f_{ady}$ 
 +    * Diffusion dissolved gas flux $f_{addx}$$f_{addy}$ 
 +    * Advection gaseous gas flux $f_{agx}$, $f_{agy}$
   * Internal variables:   * Internal variables:
     * Internal variables of the mechanical law     * Internal variables of the mechanical law
     * Internal variables of the flow law     * Internal variables of the flow law
  
elements/ehmic.1700733705.txt.gz · Last modified: 2023/11/23 11:01 by gilles