Differences

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

--- laws:licha [2019/06/25 15:18]
helene [Case 1: 2-D state (LICHA or SUCHA)]
+++ laws:licha [2021/07/28 16:44] (current)
helene [The model]
@@ Line 5: / Line 5: @@
 ==== The model ====
-Definition of a uniformly distributed load (whether normal or tangent, whether in global axis) on a line (LICHA element) or on a surface (SUCHA element).
+Definition of a uniformly distributed load (whether normal or tangent, whether in global axis) on a line ([[elements:licha|LICHA]] element) or on a surface ([[elements:sucha|SUCHA]] element).
 ==== Files ====
-Write here the names of the main subroutines of the law (those called by loi2 for Lagamine)
+Write here the names of the main subroutines of the law (those called by loi2 for Lagamine) \\
 Prepro: LLICHA.F \\
 Lagamine: LICHAB.F/SUCHAB.F (element subroutine)
@@ Line 140: / Line 140: @@
 {{ :laws:licha_4.png?500 |}}
 During the non-linear analysis, PRESSF and TAUF will be multiplied by the load factor ALAMBF (= FMULT cumulated) while PRESSD and TAUD will be multiplied by the load factor ALAMBD (= DMULT cumulated).
+=== 1.5 Y-dependent functions – Type 2 ===
+^Line 1 (I5)^^
+|IVAL |= ± 10 in local axis|
+|:::|= ± 11 in global axis|
+^Line 2 (4G10.0)^^
+|PRESSF|  Normal pressure = (PRESSF + PRESSA + PRESSB*y + PRESSC*y²)*FMULT + PRESSD*DMULT   |
+|TAUF|:::|
+|PRESSD |:::|
+|TAUD| :::|
+^Line 3 (6G10.0)^^
+|PRESSA|  Tangent load = (TAUF + TAUA + TAUB*y + TAUC*y2)*FMULT + TAUD*DMULT \\ y is the ordinate of the node where these stresses are applied  |
+|TAUA|:::|
+|PRESSB |:::|
+|TAUB| :::|
+|PRESSC |:::|
+|TAUC| :::|
+If $F_n(y).F_n(y=0) < 0$ and IVAL = + 10 => $F_n= 0$ and $F_r= 0$ \\
+If $F_n(y).F_n(y=0) > 0$ and IVAL = - 10 => $F_n= 0$ and $F_r= 0$ \\
+During the non-linear analysis, PRESSF and TAUF will be multiplied by the load factor ALAMBF (= FMULT cumulated) while PRESSD and TAUD will be multiplied by the load factor ALAMBD (= DMULT cumulated).
+=== 1.6 Imposition gravity by SUCHA ===
+^Line 1 (4I5)^^
+|IVAL |= 9|
+|IMLIC| = 0 normal load for are-length method|
+|:::|= 1 dead load|
+|ISRW| = Formulation index for Sr,w|
+|IMULT| = FMULT/DMULT multiplicator index|
+|:::|= 0 the gravity will be multiplied by FMULT|
+|:::|= 1 the gravity will be multiplied by DMULT (Used only if ISRW ≠ 0)|
+**If ISRW = 0**
+^Line 2 (2G10.0) ^^
+|PRESSF | Normal pressure|
+|PRESSD | Normal pressure|
+**If ISRW ≠ 0**
+^Line 2 (3G10.0)^^
+|GAMMAS|= $\rho_s.g$ [N/m<sup>3</sup>]|
+|GAMMAW| = $\rho_w.g$ [N/m<sup>3</sup>]|
+|POROS| = Soil porosity n|
+^Line 3 (7G10.0)^^
+|CSR1| = 1st coefficient of the function S<sub>r,w</sub>|
+|CSR2| = 2nd coefficient of the function S<sub>r,w</sub>|
+|CSR3| = 3rd coefficient of the function S<sub>r,w</sub>|
+|CSR4 |= 4th coefficient of the function S<sub>r,w</sub>|
+|SRES |= Residual saturation degree ( = S<sub>res</sub> )|
+|SRFIELD |= Field saturation degree ( = S<sub>r,field</sub> )|
+|AIREV |= Air entry value [Pa]|
+__Remarks__: \\
+With IVAL = 9 and ISRW = 0, the gravity ($\rho g$) is imposed by mean of "PRESSF" multiplied by ALAMBF or "PRESSD" multiplied by ALAMBD. \\
+With IVAL = 9 and ISRW ≠ 0, the gravity ($\rho g= (1-n) \rho_s*g+S_r*n*\rho_w *g$) is imposed and multiplied by ALAMBF if IMULT = 0, or by ALAMBD if IMULT = 1. \\
+The obtained value (positive) will be used to impose gravity forces in the –Y direction. The SUCHA element is used so of course in a 2D analysis.
 ==== Case 2: 3-D state (SUCHA) ====
+The ξ and η directions correspond to the intrinsic coordinates of the SUCHA element. \\
+Reminder: for SUCHA element, the positive axis of the tangential stresses is in the opposite direction of the axes ξ (positively oriented from IP (integration point) 1 to IP 2) and η (positively oriented from IP 1 to IP 3). The positive axis of the pressure is also in the opposite direction of the positive normal to the element (given by the corkscrew rule).
+{{ :elements:sucha_1.png?400 |}}
+^Line 1 (I5)^^
+|IVAL| = 0, 10 index telling if the function is constant (IVAL=0) or Z-dependent function, in local axis (IVAL=10)|
+|:::|= 11, 12 or 13 index telling if the function is constant, in global axis (be careful, at least one component or the normal unit of the SUCHA element must be equal to zero): \\ 11 = SUCHA with normal unit is in YZ plane. \\ 12 = SUCHA with normal unit is in ZX plane. \\ 13 = SUCHA with normal unit is in XY plane.|
+|:::|= 14: idem 11, 12, or 13, but no restriction on the orientation of the unit normal vector (more general approach).|
+=== 2.1 If IVAL = 0 (constant function, load in local axis) ===
+^Line 2 (9G10.0)^^
+|PRESSF | normal pressure|
+|TKSIF | tangent load in the ξ direction|
+|TETAF | tangent load in the η direction|
+|PRESSD | normal pressure|
+|TKSID | tangent load in the ξ direction|
+|TETAD | tangent load in the η direction|
+|PRESSC | constant normal pressure|
+|TIME1 | time of decreasing beginning|
+|TIME2 | time of decreasing end|
+These data are reference values for the distributed loads. \\
+During the non-linear analysis PRESSF, TKSIF, TETAF will be multiplied by the load factor ALAMBF while PRESSD, TKSID, TETAD will be multiplied by the load factor ALAMBD. The constant normal pressure PRESSC will be also added. \\
+  * If time < TIME1:
+    * PRESS = PRESS
+    * TKSI = TKSI
+    * TETA = TETA
+  * If TIME1 < time < TIME2:
+    * PRESS = PRESS*(TIME2-time)/(TIME2-TIME1)
+    * TKSI = TKSI *(TIME2-time)/(TIME2-TIME1)
+    * TETA = TETA *(TIME2-time)/(TIME2-TIME1)
+  * If time > TIME2:
+    * PRESS = 0
+    * TKSI = 0
+    * TETA = 0
+=== 2.2 If IVAL = 10 ===
+^1 Blank line^^
+^Line 3 (2G10.0)^^
+|PRESSA | PRESS = (PRESSA + PRESSB*ZINI)* ALAMBF  |
+|PRESSB | :::|
+=== 2.3 If IVAL = 11, 12, 13 (constant function, load in global axis) ===
+^Line 2 (6G10.0)^^
+|SIGXF | (*Fmult) Total stress in X axis|
+|SIGYF | (*Fmult) Total stress in Y axis|
+|SIGZF | (*Fmult) Total stress in Z axis|
+|SIGXD | (*Dmult) Total stress in X axis|
+|SIGYD | (*Dmult) Total stress in Y axis|
+|SIGZD | (*Dmult) Total stress in Z axis|
+=== 2.4 If IVAL = 14 (constant function, load in global axis) ===
+^Line 2 (6G10.0)^^
+|SIGXF | (*Fmult) Total stress in X axis|
+|SIGYF | (*Fmult) Total stress in Y axis|
+|SIGZF | (*Fmult) Total stress in Z axis|
+|SIGXD | (*Dmult) Total stress in X axis|
+|SIGYD | (*Dmult) Total stress in Y axis|
+|SIGZD | (*Dmult) Total stress in Z axis|
+^Line 3 (5G10.0)^^
+|SIGXC | Constant total stress in X axis|
+|SIGYC | Constant total stress in Y axis|
+|SIGZC | Constant Total stress in Z axis|
+|TIME1 | Time of decreasing beginning|
+|TIME2 | Time of decreasing end|
+  * If time < TIME1:
+    * SIGX = SIGX
+    * SIGY = SIGY
+    * SIGZ = SIGZ
+  * If TIME1 < time < TIME2:
+    * SIGX = SIGX *(TIME2-time)/(TIME2-TIME1)
+    * SIGY = SIGY *(TIME2-time)/(TIME2-TIME1)
+    * SIGZ = SIGZ *(TIME2-time)/(TIME2-TIME1)
+  * If time > TIME2:
+    * SIGX = 0
+    * SIGY = 0
+    * SIGZ = 0
 ===== Stresses =====
 ==== Number of stresses ====

Lagamine

User Tools

Site Tools

Differences

Page Tools