By the way,
I was thinking that the problem could be associated to the MSP. I have
assigned the Elasticity properties as follows.
GINA - Material Solid Properties
#SOLID_PROPERTIES
$NAME
MSPFAULT
$DENSITY
1 2.85e+03
$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
$THERMAL
EXPANSION
1.000e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000
$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
$GRAVITY_CONSTANT
0
$SOLID_BULK_MODULUS
1.e-9
#SOLID_PROPERTIES
$NAME
MSPFRACT2
$DENSITY
1 2.85e+003
$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
$THERMAL
EXPANSION
1.00e-005
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000
; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
$GRAVITY_CONSTANT
0
$SOLID_BULK_MODULUS
1.e-9
#SOLID_PROPERTIES
$NAME
MSPFRACT3
$DENSITY
1 2.85e+03
$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
$THERMAL
EXPANSION
1.000e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000
; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
$GRAVITY_CONSTANT
0
$SOLID_BULK_MODULUS
1.e-9
#SOLID_PROPERTIES
$NAME
MSPFRACT4
$DENSITY
1 2.85e+03
$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+09
$THERMAL
EXPANSION
1.00e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000
; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
$GRAVITY_CONSTANT
0
$SOLID_BULK_MODULUS
1.e-9
#SOLID_PROPERTIES
$NAME
MSPROCKMATRIX
$DENSITY
1 2.850e+003
$ELASTICITY
POISSION 0.25
YOUNGS_MODULUS
1 6.00e+010
$THERMAL
EXPANSION
1.00e-05
CAPACITY
1 1.017e+003
CONDUCTIVITY
1 3.00e+000
; BISHOP-BIOT Model:
; Sigma_tot = Sigma_eff - Biot*bishop(S)*p_w
$BIOT_CONSTANT
1.0
BISHOP_COEFFICIENT
3 0.0 ; model_3: bishop=0.0 or 1.0
$GRAVITY_CONSTANT
0
$SOLID_BULK_MODULUS
1.e-9
#STOP
$PLASTICITY
MOHR-COULOMB
4903500
45 ;friction angle
0
2451750
0
----------------------------------------------------------------------------------------------------------
I was wondering what characteristic I am missing. Are the friction
angles in radiant? I was also thinking that the problem might be in the
numerical NUM. Do I have to consider the overall coupling so that the
thermal effects cause strain in my material?
Thanks!
On Wednesday, 13 July 2016 09:05:58 UTC+1, Wenqing Wang wrote:
Then you have to input the initial stress data as
#INITIAL_CONDITION
$PCS_TYPE
DEFORMATION
$PRIMARY_VARIABLE
STRESS_XX
$GEO_TYPE
SUB_DOMAIN
5
0 [expression of stress]
1 [expression of stress]
2 [expression of stress]
3 [expression of stress]
4 [expression of stress]
On 07/12/2016 07:35 PM, Jhonny Encalada wrote:
> Dear Wenqing Wang
>
> Thanks for your help,
>
> Regarding to your explanation about the sub-domain. I have 5
> material groups actually; however I decided to assign the initial
> conditions to the WHOLE DOMAIN.
> My material groups are numbered from 0,1,.. to 4 and four
> materials are discrete fractures. The material No4 is the rock
> matrix. Let me know if I am mistaken please.
> Thanks again!
>
> Jhonny
>
> On Tuesday, 12 July 2016 17:09:59 UTC+1, Wenqing Wang wrote:
>
> * The sign of stress or strain takes the conventional solid
> mechanics' one, i.e. + for tension.
> * z takes the value of z coordinate of an node defined in the
> mesh data.
>
> The initial stresses can be input as (in a benchmark directory
> of THM/init)
> GeoSys-IC: Initial Conditions
> ------------------------------------------------
> #INITIAL_CONDITION
> $PCS_TYPE
> DEFORMATION
> $PRIMARY_VARIABLE
> STRESS_XX
> $GEO_TYPE
> SUB_DOMAIN
> 3
> 0 -32.1e6+0.55e5*y
> 1 -32.1e6+0.55e5*y
> 2 -32.1e6+0.55e5*y
> #INITIAL_CONDITION
> $PCS_TYPE
> DEFORMATION
> $PRIMARY_VARIABLE
> STRESS_YY
> $GEO_TYPE
> SUB_DOMAIN
> 3
> 0 -13243500+26487*y
> 1 -13243500+26487*y
> 2 -13243500+26487*y
> #INITIAL_CONDITION
> $PCS_TYPE
> DEFORMATION
> $PRIMARY_VARIABLE
> STRESS_ZZ
> $GEO_TYPE
> SUB_DOMAIN
> 3
> 0 -10.6e6+0.2e5*y
> 1 -10.6e6+0.2e5*y
> 2 -10.6e6+0.2e5*y
> #STOP
>
> Since one one domain is shown in your data, you input can be,
> for example,
> #INITIAL_CONDITION
> $PCS_TYPE
> DEFORMATION
> $PRIMARY_VARIABLE
> STRESS_YY
> $GEO_TYPE
> SUB_DOMAIN
> 1
> 0 -1.176e06+22.491e03*z
>
> On 07/12/2016 05:06 PM, Jhonny Encalada wrote:
>> Dear Users and Developers,
>>
>> I am running a 3 model for HTM processes. I have completed
>> the HT models and seem to be working well; however, when I
>> implemented the third effect (M); I am not getting any
>> solution at all. I have assigned the initial conditions as
>> they are showed in the benchmarks, for example:
>>
>> GINA - Initial Condition
>> #INITIAL_CONDITION
>> $PCS_TYPE
>> HEAT_TRANSPORT
>> $PRIMARY_VARIABLE
>> TEMPERATURE1
>> $GEO_TYPE
>> DOMAIN
>> $DIS_TYPE
>> CONSTANT 453.15
>> #INITIAL_CONDITION
>> $PCS_TYPE
>> LIQUID_FLOW
>> $PRIMARY_VARIABLE
>> PRESSURE1
>> $GEO_TYPE
>> DOMAIN
>> $DIS_TYPE
>> GRADIENT 9810.00
>> #INITIAL_CONDITION
>> $PCS_TYPE
>> DEFORMATION
>> $PRIMARY_VARIABLE
>> STRESS_XX
>> $GEO_TYPE
>> DOMAIN
>> $DIS_TYPE
>> FUNCTION -1.854e06+14.513e03*z
>> #INITIAL_CONDITION
>> $PCS_TYPE
>> DEFORMATION
>> $PRIMARY_VARIABLE
>> STRESS_YY
>> $GEO_TYPE
>> DOMAIN
>> $DIS_TYPE
>> FUNCTION -1.176e06+22.491e03*z
>> #INITIAL_CONDITION
>> $PCS_TYPE
>> DEFORMATION
>> $PRIMARY_VARIABLE
>> STRESS_ZZ
>> $GEO_TYPE
>> DOMAIN
>> $DIS_TYPE
>> FUNCTION 0.3e06+25.60e03*z
>> #STOP
>>
>> Since my calculations are at 4km depth, I was wondering what
>> is the correct sign (+/-) for stresses that OGS recognizes.
>> Is it (-) for compression and (+) for tension liek in
>> convetional mechanics of materials. Or is it (+) for
>> compression and negative for tension. Furthermore, what are
>> the implications when using equations; is the depth(z) taken
>> with the sign I modeled the geometry e.g. -4000m? I assume
>> that all the units need to be transformed into Pascals like
>> the pressures for the Hydraulic simulation.
>> Please let me know your opinions and advice. I will be very glad.
>>
>> Jhonny Encalada
>> --
>> You received this message because you are subscribed to the
>> Google Groups "ogs-users" group.
>> To unsubscribe from this group and stop receiving emails from
>> it, send an email to ogs-users+...@googlegroups.com.
>> For more options, visit https://groups.google.com/d/optout
>> <https://groups.google.com/d/optout>.
>
> --
> You received this message because you are subscribed to the Google
> Groups "ogs-users" group.
> To unsubscribe from this group and stop receiving emails from it,
> send an email to ogs-users+...@googlegroups.com <javascript:>.
> For more options, visit https://groups.google.com/d/optout
> <https://groups.google.com/d/optout>.
--
You received this message because you are subscribed to the Google
Groups "ogs-users" group.
To unsubscribe from this group and stop receiving emails from it, send
an email to ogs-users+unsubscribe@googlegroups.com
<mailto:ogs-users+unsubscribe@googlegroups.com>.
For more options, visit https://groups.google.com/d/optout.