//
// Creation du namespace erable_Siau
//
{ Erable_Siau
//
// Materiau anisotrope. On suppose que 
// L == axe des x
// T == axe des y
// R == axe des z
//

// Temperature dans la piece
scalaire TInfini 21

// Suivant labo
scalaire specific_gravity .634

// Suivant Siau 
scalaire porosity f(specific_gravity,M)=1-specific_gravity*(0.653+M/100)
scalaire D_BT f(M,TInfini)=0.000007*exp( (70*M-9200)/(2.*(TInfini+273.15)) )
scalaire D_T f(D_BT,porosity)=D_BT/((1-porosity)*(1-sqrt(porosity)))
scalaire D_L f(D_T)=2.5*D_T

//
// On fabrique le module de Young, le cisaillement la difussion et les 
// coefficients de dilatation hydrique à partir du tableau Blanchet
//
// Valeurs du tableau:

scalaire        rho             f(specific_gravity)=specific_gravity*1000

scalaire        "DL"            f(D_L)=D_L
scalaire        "DR"            f(D_T)=D_T
scalaire        "DT"            f(D_T)=D_T

scalaire        "EL"            13.810e9
scalaire        "ER"            1.311e9
scalaire        "ET"            0.678e9

scalaire        "GLR"           1.013e9
scalaire        "GLT"           0.753e9
scalaire        "GRT"           0.255e9

scalaire        "VTR"           0.42
scalaire        "VLR"           0.46
scalaire        "VLT"           0.50

scalaire        "BdeL"            1.8e-4
scalaire        "BdeR"            1.9e-3
scalaire        "BdeT"            3.39e-3

scalaire        "BadL"            1.5e-4
scalaire        "BadR"            2.1e-3
scalaire        "BadT"            3.63e-3

// Pour la dilatation qui est BetaM
scalaire "BL" f(M,MP,BdeL,BadL)=(M <= MP)*BdeL+(M > MP)*BadL
scalaire "BT" f(M,MP,BdeT,BadT)=(M <= MP)*BdeT+(M > MP)*BadT
scalaire "BR" f(M,MP,BdeR,BadR)=(M <= MP)*BdeR+(M > MP)*BadR

scalaire        C_M             0.01

//
// Les données pour le programme:
//

scalaire      C_MRho f(rho,C_M)=rho*C_M

scalaire        "KL"            f(DL,C_MRho)=DL*C_MRho
scalaire        "KR"            f(DR,C_MRho)=DR*C_MRho
scalaire        "KT"            f(DT,C_MRho)=DT*C_MRho

empilementto2 K_M    [KL,KT,KR,Zero,Zero,Zero]

to4elasticite Eijkl  [EL,ET,ER,VLT,VTR,VLR,GLT,GRT,GLR]

empilementto2 BetaM  [BL,BT,BR,Zero,Zero,Zero]

//
// Pour les C.L.
//

scalaire  H_M     -3.2e-4
scalaire  MAmbiant 5.8

}