1. Reference problem#
1.1. Geometry#
Point coordinates (in meters):
\(A\) |
|
|
|
|
\(x\) |
0.5 |
|||
\(y\) |
0.5 |
|||
\(z\) |
0.5 |
1.2. Material properties for model LETK#
PA = 0.1
NELAS = 0.
SIGMA_C = 12.
H0_ EXT = 1.10292
GAMMA_CJS = 0.8
XAMS = 0.1
ETA = 0.04
A_0 = 0.25
A_E = 0.60
A_ PIC = 0.4
S_0 = 0.0005
M_0 = 0.01
M_E = 2.
M_ PIC = 6.
M_ ULT = 0.61
XI_ULT = 0.365
XI_E = 0.028
XI_PIC = 0.015
MV_MAX = 3.
XIV_MAX = 0.0039
A = 1.5e-12
N = 4.5
SIGMA_P1 = 57.8
MU0_V = 0.1
XI0_V = 0.3
MU1 = 0.1
XI1 = 0.3
1.3. Material properties for model LKR#
PA = .1E6
GAMMA =.85
M_0 = .5
F_P = 0.136047510046
M_1 = 9.69880017363
SIGMA_C = 10.9985715832E6
A_2 = 0.580184800258
Q_I = 100.000648048E6
V_1 = 1.5
V_2 = 1.5
XI_1 = 1.e-2
XI_2 = 1.8e-2
XI_5 = 1.6e-2
A = 1.e-14
N = 3.5
RHO_1 = 1.
RHO_2 = 0.1
RHO_4 = 1.10668567265
R_Q = 1.e-6
R_M = 1.e-6
R_S = 1.e-6
R_X1 = 1.e-6
R_X2 = 1.e-6
R_X5 = 1.e-6
Z = 1.e-6
COUPLAGE_P_VP = 1
A_ SIGC = 0.155495602806
B_ SIGC = 4.69721443803
1.4. Material properties for model NLH_CSRM#
Material parameters are given in the International System of Units (SI).
YoungModulus=7.0E9
Fish ratio = 0.3
Isocomplaslim=50.0E6
Isotenselaslim=0.1E6
MCCSlopeCSL =0.5
NLHIndex =1.0
MbigocritCoef=10.0
abigocritCoef=0.75
IncompIndex=15.0
Tau=2.0e2
PerzynaExpo=2.0
NLHModulusP =7.0e9/2.5
NLHModulusV =0.01*7.0e9
1.5. Initial conditions, boundary conditions, and loading#
Phase 1:
The sample is brought to a homogeneous state: \({\sigma }_{\text{xx}}^{0}={\sigma }_{\text{yy}}^{0}={\sigma }_{\text{zz}}^{0}\), by imposing the corresponding confinement pressure on the front, right lateral and upper faces. The movements are blocked on the back (\({u}_{x}=0\)), left side (\({u}_{y}=0\)) and bottom (\({u}_{z}=0\)) faces.
Phase 2:
The movements are kept blocked on the rear (\({u}_{x}=0\)), left lateral (\({u}_{y}=0\)) and lower (\({u}_{z}=0\)) faces. On all sides, the water pressure is zero.
An imposed displacement is applied on the upper face in order to obtain a deformation \({\mathrm{\epsilon }}_{\text{zz}}=-6\text{\%}\) (counted from the start of phase 2). On the front and right lateral faces, a \(5\mathit{MPa}\) constraint is imposed.
For the E, F and G models, an increasing temperature is imposed via command AFFE_CHAR_MECA between 0°C and 50°C in phase 1 and 50°C and 100°C in phase 2.
For H modeling, the deformation at the end of the test is a hundred times greater than that obtained at the end of the confinement phase: –0.00028571428571428416 * 100.