# Temps de depart: scalaire T_0 0 # Temps final scalaire T_F 600000000 # Pas de temps initial scalaire PasDeTempsInit 0.001 # Pas de temps minimal scalaire PasDeTempsMin 0.001 # Pas de temps maximal (pas le limite = 0) scalaire PasDeTempsMax 0 # Nombre maximal de pas de temps scalaire NbPasDeTemps 9999 # frequence du calcul en elasticite scalaire freqCalcul 1 # frequence d'impression des champs scalaire freqImpression 1 # tolerance pour le calcul des pas de temps scalaire TolerancePasVariable 5e-3 # Tolerance pour la convergence des solveurs scalaire ToleranceNlin 1e-9 scalaire ToleranceLin 5e-9 scalaire Divergence 1e+10 # nombre maximum d'iterations pour les solveur scalaire MaxItNlin 100 scalaire MaxItLin 5000 # temps impose pour le calcul booleen TempsImpose v # Propriete des matrices booleen MDependTemps f booleen MDependSolution f booleen KDependTemps f booleen KDependSolution v booleen KDependCL v ## Solveur PETSc ## solveurlinpetsc typesolveur "SolveurLin" gmres solveurlinpetsc typeprecond "SolveurLin" ilu solveurlinpetsc suiviconvergence "SolveurLin" aucun solveurlinpetsc tolerance "SolveurLin" [ToleranceLin,ToleranceLin,Divergence,MaxItLin] solveurinstnlinpetsc typesolveur "SolveurInst" solveurinstnlinpetsc solveurlin "SolveurInst" SolveurLin solveurinstnlinpetsc pasdetempsvariable "SolveurInst" [PasDeTempsMin,PasDeTempsMax,T_F,TolerancePasVariable] solveurinstnlinpetsc parametresnlin "SolveurInst" [ToleranceNlin,MaxItNlin] solveurinstnlinpetsc parametresinst "SolveurInst" [PasDeTempsInit, NbPasDeTemps, T_0, eulerimplicit] solveurlinpetsc typesolveur "SolveurMec" direct #solveurlinpetsc typeprecond "SolveurMec" ilu solveurlinpetsc suiviconvergence "SolveurMec" defaut solveurlinpetsc tolerance "SolveurMec" [ToleranceLin,ToleranceLin,Divergence,MaxItLin] ## Definition des champs du probleme champgeolin "ChampGeo" scalaire "zero" 0. scallin "M" 6.5 scallin "MPrec" 6.5 scallin "M_t" 0. v3dquad "U" [0, 0, 0] booleen "MtQuad" false booleen "MtLin" true ## ## Les valeurs exprimes en fct de la densite comme la table 1 article no4 ## ## Les valeurs de densité possibles (/100) scalaire "densite1" 5.40 scalaire "densite2" 6.50 scalaire "densite3" 8.00 scalaire "densite4" 9.00 ## Niveau d'humidite de référence pour les valeurs scalaire MREF0 6.6 scalaire MREF1 9.0 ## Module d'Young (lineaire alors on donne E au point de reference M = MREF et la pente de la droite i.e. (E_1(M=9)-E110)/(M=9-MREF)) ## E1 scalaire E110 7.4e+8 scalaire E111 6.3e+8 scalaire DE11 -4.58333e+7 scalaire E120 1.3e+9 scalaire E121 1.1e+9 scalaire DE12 -8.33333e+7 scalaire E130 2.2e+9 scalaire E131 1.9e+9 scalaire DE13 -1.25e+8 scalaire E140 2.8e+9 scalaire E141 2.5e+9 scalaire DE14 -1.25e+8 ## E3 scalaire E310 1.7e+7 scalaire E311 1.2e+7 scalaire DE31 -2.08333e+6 scalaire E320 3.5e+7 scalaire E321 2.6e+7 scalaire DE32 -3.75e+6 scalaire E330 6.1e+7 scalaire E331 4.4e+7 scalaire DE33 -7.08333e+6 scalaire E340 7.8e+7 scalaire E341 5.6e+7 scalaire DE34 -9.16667e+6 ## Poisson scalaire "Nu12" 0.3 scalaire "Nu13" 0.2 ## Cisaillament (lineaire alors on donne G au point de reference et la pente de la droite) ## G1 scalaire G110 2.8e+8 scalaire G111 2.4e+8 scalaire DG11 -1.66667e+7 scalaire G120 4.9e+8 scalaire G121 4.5e+8 scalaire DG12 -1.66667e+7 scalaire G130 8.5e+8 scalaire G131 7.5e+8 scalaire DG13 -4.16667e+7 scalaire G140 1.06e+9 scalaire G141 9.4e+8 scalaire DG14 -6.66667e+7 ## G3 scalaire G310 8.4e+7 scalaire G311 7.12e+7 scalaire DG31 -5.41667e+6 scalaire G320 1.17e+8 scalaire G321 9.1e+7 scalaire DG32 -1.25e+7 scalaire G330 1.34e+8 scalaire G331 1.02e+8 scalaire DG33 -1.25e+7 scalaire G340 1.92e+8 scalaire G341 1.4e+8 scalaire DG34 -2.08333e+7 ## Diffusion scalaire D10 1.0e-9 scalaire D11 9.9e-11 scalaire DD1 -3.7541e-10 scalaire D20 9.3e-10 scalaire D21 8.9e-11 scalaire DD2 -3.50416e-10 scalaire D30 8.0e-10 scalaire D31 8.0e-11 scalaire DD3 -3.0e-10 scalaire D40 7.2e-10 scalaire D41 7.6e-11 scalaire DD4 -2.68333e-10 ## Les valeurs de dilatation ## Absorption scalaire "Beta11Abs" 1.7e-4 scalaire "Beta12Abs" 2.3e-4 scalaire "Beta13Abs" 3.6e-4 scalaire "Beta14Abs" 3.9e-4 scalaire "Beta31Abs" 7.4e-3 scalaire "Beta32Abs" 7.7e-3 scalaire "Beta33Abs" 8.4e-3 scalaire "Beta34Abs" 8.7e-3 ## Dessorption scalaire "Beta11Des" 3.4e-4 scalaire "Beta12Des" 4.8e-4 scalaire "Beta13Des" 7.7e-4 scalaire "Beta14Des" 10.0e-4 scalaire "Beta31Des" 5.4e-3 scalaire "Beta32Des" 6.2e-3 scalaire "Beta33Des" 9.5e-3 scalaire "Beta34Des" 1.1e-2 ## ## Fin du tableau on construit les fonctions ## ## Pour les CL ## scalaire "h" -3.2e-5 #scalaire "M_infini" 12.5 scalaire "M_infini" f()=12.5 + ((t>540000)*(t<626400)+(t>2160000)*(t<2332800))*1.5 ### A Priori le reste ne doit pas etre modifie SAUF la densiste ( on suppose 12 couches et 4 densites) {"d_b1" scalaire densite f(densite1)=densite1 } {"d_b2" scalaire densite f(densite2)=densite2 } {"d_b3" scalaire densite f(densite3)=densite3 } {"d_b4" scalaire densite f(densite4)=densite4 } ## la densite est defini par couches: #CAS 1: #scalsurgeom "densite" \ #couche1=d_b1 couche2=d_b1 couche3=d_b1 couche4=d_b1 couche5=d_b1 couche6=d_b1 \ #couche7=d_b1 couche8=d_b1 couche9=d_b1 couche10=d_b1 couche11=d_b1 couche12=d_b1 defaut=d_b1 #CAS 2 #scalsurgeom "densite" \ #couche1=d_b2 couche2=d_b2 couche3=d_b2 couche4=d_b2 couche5=d_b2 couche6=d_b2 \ #couche7=d_b2 couche8=d_b2 couche9=d_b2 couche10=d_b2 couche11=d_b2 couche12=d_b2 defaut=d_b2 #CAS 3 #scalsurgeom "densite" \ #couche1=d_b3 couche2=d_b3 couche3=d_b3 couche4=d_b3 couche5=d_b3 couche6=d_b3 \ #couche7=d_b3 couche8=d_b3 couche9=d_b3 couche10=d_b3 couche11=d_b3 couche12=d_b3 defaut=d_b3 #CAS 4 scalsurgeom "densite" \ couche1=d_b4 couche2=d_b4 couche3=d_b3 couche4=d_b3 couche5=d_b2 couche6=d_b1 \ couche7=d_b1 couche8=d_b2 couche9=d_b3 couche10=d_b3 couche11=d_b4 couche12=d_b4 defaut=d_b1 #CAS 5 #scalsurgeom "densite" \ #couche1=d_b2 couche2=d_b2 couche3=d_b2 couche4=d_b2 couche5=d_b2 couche6=d_b2 \ #couche7=d_b2 couche8=d_b2 couche9=d_b2 couche10=d_b2 couche11=d_b3 couche12=d_b3 defaut=d_b2 ## caracteristique associee a chaque densite scalaire "Khi1" f(densite,densite1)=(densite<=densite1)*(densite>=densite1) scalaire "Khi2" f(densite,densite2)=(densite<=densite2)*(densite>=densite2) scalaire "Khi3" f(densite,densite3)=(densite<=densite3)*(densite>=densite3) scalaire "Khi4" f(densite,densite4)=(densite<=densite4)*(densite>=densite4) scalaire "beta1Abs" f(Khi1,Khi2,Khi3,Khi4,Beta11Abs,Beta12Abs,Beta13Abs,Beta14Abs)=Khi1*Beta11Abs + Khi2*Beta12Abs + Khi3*Beta13Abs + Khi4*Beta14Abs scalaire "beta3Abs" f(Khi1,Khi2,Khi3,Khi4,Beta31Abs,Beta32Abs,Beta33Abs,Beta34Abs)=Khi1*Beta31Abs + Khi2*Beta32Abs + Khi3*Beta33Abs + Khi4*Beta34Abs scalaire "beta1Des" f(Khi1,Khi2,Khi3,Khi4,Beta11Des,Beta12Des,Beta13Des,Beta14Des)=Khi1*Beta11Des + Khi2*Beta12Des + Khi3*Beta13Des + Khi4*Beta14Des scalaire "beta3Des" f(Khi1,Khi2,Khi3,Khi4,Beta31Des,Beta32Des,Beta33Des,Beta34Des)=Khi1*Beta31Des + Khi2*Beta32Des + Khi3*Beta33Des + Khi4*Beta34Des scalaire D1 f(M,D10,D11,DD1,MREF0,MREF1)= (D10 + DD1*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire D2 f(M,D20,D21,DD2,MREF0,MREF1)= (D20 + DD2*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire D3 f(M,D30,D31,DD3,MREF0,MREF1)= (D30 + DD3*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire D4 f(M,D40,D41,DD4,MREF0,MREF1)= (D40 + DD4*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire E11 f(M,E110,E111,DE11,MREF0,MREF1)= (E110 + DE11*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire E12 f(M,E120,E121,DE12,MREF0,MREF1)= (E120 + DE12*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire E13 f(M,E130,E131,DE13,MREF0,MREF1)= (E130 + DE13*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire E14 f(M,E140,E141,DE14,MREF0,MREF1)= (E140 + DE14*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire E31 f(M,E310,E311,DE31,MREF0,MREF1)= (E310 + DE31*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire E32 f(M,E320,E321,DE32,MREF0,MREF1)= (E320 + DE32*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire E33 f(M,E330,E331,DE33,MREF0,MREF1)= (E330 + DE33*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire E34 f(M,E340,E341,DE34,MREF0,MREF1)= (E340 + DE34*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire G11 f(M,G110,G111,DG11,MREF0,MREF1)= (G110 + DG11*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire G12 f(M,G120,G121,DG12,MREF0,MREF1)= (G120 + DG12*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire G13 f(M,G130,G131,DG13,MREF0,MREF1)= (G130 + DG13*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire G14 f(M,G140,G141,DG14,MREF0,MREF1)= (G140 + DG14*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire G31 f(M,G310,G311,DG31,MREF0,MREF1)= (G310 + DG31*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire G32 f(M,G320,G321,DG32,MREF0,MREF1)= (G320 + DG32*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire G33 f(M,G330,G331,DG33,MREF0,MREF1)= (G330 + DG33*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire G34 f(M,G340,G341,DG34,MREF0,MREF1)= (G340 + DG34*(M-MREF0)*(M>MREF0))*(M=MREF1) scalaire "E1" f(E11,E12,E13,E14,Khi1,Khi2,Khi3,Khi4)=Khi1*E11+Khi2*E12+Khi3*E13+Khi4*E14 scalaire "E3" f(E31,E32,E33,E34,Khi1,Khi2,Khi3,Khi4)=Khi1*E31+Khi2*E32+Khi3*E33+Khi4*E34 scalaire "G1" f(G11,G12,G13,G14,Khi1,Khi2,Khi3,Khi4)=Khi1*G11+Khi2*G12+Khi3*G13+Khi4*G14 scalaire "G3" f(G31,G32,G33,G34,Khi1,Khi2,Khi3,Khi4)=Khi1*G31+Khi2*G32+Khi3*G33+Khi4*G34 scalaire "K" f(D1,D2,D3,D4,Khi1,Khi2,Khi3,Khi4)=(Khi1*D1+Khi2*D2+Khi3*D3+Khi4*D4) scalaire "beta1" f(M_t,beta1Abs,beta1Des)=(M_t>=0.)*beta1Abs + (M_t<0.)*beta1Des scalaire "beta3" f(M_t,beta3Abs,beta3Des)=(M_t>=0.)*beta3Abs + (M_t<0.)*beta3Des ## Hydrique ## empilementto2 "K_M" [zero,zero,K,zero,zero,zero] ## Mecanique ## to4elasticite "Eijkl" [E1,E1,E3,Nu12,Nu13,Nu13,G1,G3,G3] empilementto2 "Beta" [beta1,beta1,beta3,zero,zero,zero]