1 : * @MATETHM PROCEDUR FANDEUR 10/12/17 21:18:00 6427 2 : * PROCEDURE @MATETHM 3 : *---------------------------------------------------------------------* 4 : * NOM : @MATETHM * 5 : * * 6 : * DESCRIPTION : * 7 : * * 8 : * SYNTAXE : * 9 : * * 10 : * * 11 : * LANGAGE : GIBIANE-CAST3M * 12 : * AUTEUR : S. PASCAL (CEA/DEN/DM2S/SEMT/LM2S) * 13 : * COURRIEL : serge.pascal@cea.fr * 14 : *---------------------------------------------------------------------* 15 : * VERSION : v1, xx/xx/xx08, version initiale * 16 : * HISTORIQUE : v1, xx/xx/xx08, creation * 17 : * HISTORIQUE : * 18 : * HISTORIQUE : * 19 : *---------------------------------------------------------------------* 20 : * Priere de PRENDRE LE TEMPS de completer les commentaires * 21 : * en cas de modification de ce sous-programme afin de faciliter * 22 : * la maintenance ! * 23 : *---------------------------------------------------------------------* 24 : 'DEBP' @MATETHM ; 25 : * * 26 : 'ARGU' MOD1*'MMODEL' CH0*'CHPOINT' ; 27 : * * 28 : CH0 = 'REDU' ('EXTR' MOD1 'MAIL') CH0 ; 29 : * * 30 : CHPG0 = 'EXCO' CH0 'PG' ; 31 : CHPC0 = 'EXCO' CH0 'PC' ; 32 : CHTK0 = 'EXCO' CH0 'T ' ; 33 : GRCH0 = 'GRAD' MOD1 CH0 ; 34 : * * 35 : * LISTMOTS des composantes du gradient : * 36 : VDIM1 = 'VALE' 'DIME' ; 37 : 'SI' ('EGA' VDIM1 3) ; 38 : LMGPG1 = 'MOTS' 'PG,X' 'PG,Y' 'PG,Z' ; 39 : LMGPC1 = 'MOTS' 'PC,X' 'PC,Y' 'PC,Z' ; 40 : LMGTK1 = 'MOTS' 'T,X' 'T,Y' 'T,Z' ; 41 : LMSCA1 = 'MOTS' 'SCAL' 'SCAL' 'SCAL' ; 42 : 'SINO' ; 43 : VMOD1 = 'VALE' 'MODE' ; 44 : 'SI' ('EGA' VMOD1 'AXIS') ; 45 : LMGPG1 = 'MOTS' 'PG,R' 'PG,Z' ; 46 : LMGPC1 = 'MOTS' 'PC,R' 'PC,Z' ; 47 : LMGTK1 = 'MOTS' 'T,R' 'T,Z' ; 48 : 'SINO' ; 49 : LMGPG1 = 'MOTS' 'PG,X' 'PG,Y' ; 50 : LMGPC1 = 'MOTS' 'PC,X' 'PC,Y' ; 51 : LMGTK1 = 'MOTS' 'T,X' 'T,Y' ; 52 : 'FINS' ; 53 : LMSCA1 = 'MOTS' 'SCAL' 'SCAL' ; 54 : 'FINS' ; 55 : * * 56 : * D'ou les gradients de chaque composante : * 57 : GRPG0 = 'EXCO' GRCH0 LMGPG1 ; 58 : GRPC0 = 'EXCO' GRCH0 LMGPC1 ; 59 : GRTK0 = 'EXCO' GRCH0 LMGTK1 ; 60 : * 61 : TC0 = 273.15 ; 62 : CHTC0 = CHTK0 - TC0 ; 63 : CHPGSC0 = CHPG0 'NOMC' 'SCAL' ; 64 : CHPCSC0 = CHPC0 'NOMC' 'SCAL' ; 65 : CHTKSC0 = CHTK0 'NOMC' 'SCAL' ; 66 : CHTCSC0 = CHTC0 'NOMC' 'SCAL' ; 67 : * * 68 : *---------------------------------------------------------------------* 69 : * * 70 : * CARACTERISTIQUES MATERIELLES * 71 : * * 72 : *---------------------------------------------------------------------* 73 : * * 74 : *-------------------- CARACTERISTIQUES CONSTANTES --------------------* 75 : * * 76 : * RGP1 : Constante des gaz parfaits (J/mol/K) * 77 : * MMGA1 : Masse molaire de l'air (kg/m3) * 78 : * MMGW1 : Masse molaire de l'eau (kg/m3) * 79 : * RHOS1 : Masse volumique du solide (kg/m3) * 80 : * CPW1 : Capacite calorifique a P constante de l'eau (J/kg/K) * 81 : * CPGW1 : Capacite calorifique a P constante de la vapeur (J/kg/K) * 82 : * CPGA1 : Capacite calorifique a P constante de l'air (J/kg/K) * 83 : * DHHD1 : Enthalpie de deshydratation (J/kg) * 84 : RGP1 = 8.3145 ; 85 : MMGA1 = 28.9645E-3 ; 86 : MMGW1 = 18.0153E-3 ; 87 : RHOS1 = 2625.8 ; 88 : CPW1 = 4181.0 ; 89 : CPGW1 = 1805.0 ; 90 : CPGA1 = 1005.7 ; 91 : DHHD1 = 0.5E6 / 0.24 ; 92 : * * 93 : * Variables de travail : 94 : CHTKSCM1 = CHTKSC0 ** -1 ; 95 : FGPGA0 = MMGA1 / RGP1 * CHTKSCM1 ; 96 : FGPGW0 = MMGW1 / RGP1 * CHTKSCM1 ; 97 : CHPGSCM1 = CHPGSC0 ** -1 ; 98 : CHPGSCM2 = CHPGSC0 ** -2 ; 99 : * * 100 : *------------------ CAPACITE CALORIFIQUE DU SOLIDE -------------------* 101 : * Ok ! * 102 : * ENTREE * 103 : * - Temperature (K) * 104 : * SORTIE * 105 : * - Capacite calorifique a P constante du solide (J/kg/K) : CPS0 * 106 : * PARAMETRE * 107 : * CPS1 : Capacite calorifique a P constante du solide a 0°C * 108 : * ACHI1 : Coef. de variation lineaire en temperature * 109 : CPS1 = 940. ; 110 : ACHI1 = 0.35 ; 111 : CPS0 = CPS1 + (ACHI1 * CHTCSC0) ; 112 : * * 113 : *---------------------------- SATURATION -----------------------------* 114 : * * 115 : * ENTREES * 116 : * - Temperature (K) * 117 : * - Pression capillaire (Pa) * 118 : * SORTIES * 119 : * - Satiration en eau (-) : SW0 * 120 : * - Derivee -/- Pc : DSWDPC0 * 121 : * - Derivee -/- T : DSWDT0 * 122 : * - Satiration en gaz () : SG0 * 123 : * PARAMETRES * 124 : * @SATURAT : Procedure * 125 : SW0 = @SATURAT CHTKSC0 CHPCSC0 ; 126 : *'MESS' 'SW0 =' ('MINI' SW0) ; 127 : * * 128 : * Derivee de la staturation par rapport a Pc : Ok ! 129 : DSWDPC0 = (@SATURAT CHTKSC0 (CHPCSC0 + 50.)) 130 : - (@SATURAT CHTKSC0 (CHPCSC0 - 50.)) ; 131 : DSWDPC0 = 0.01 * DSWDPC0 ; 132 : * * 133 : * Derivee de la staturation par rapport a T : 134 : DSWDT0 = (@SATURAT (CHTKSC0 + 0.05) CHPCSC0) 135 : - (@SATURAT (CHTKSC0 - 0.05) CHPCSC0) ; 136 : DSWDT0 = 10. * DSWDT0 ; 137 : * * 138 : * Saturation en phase gazeuse : SG0 * 139 : SG0 = 1. - SW0 ; 140 : * * 141 : *----------------------------- POROSITE ------------------------------* 142 : * * 143 : * ENTREE * 144 : * - Temperature (Celsius) * 145 : * SORTIES * 146 : * - Porosite (-) : N0 * 147 : * PARAMETRES * 148 : * N1 : porosite initiale (a TREF1) * 149 : * APN : coef. de variation avec la temperature * 150 : N1 = 0.1368 ; 151 : APN1 = 7.8E-5 ; 152 : * * 153 : * Porosite : 154 : N0 = N1 + (APN1 * CHTCSC0) ; 155 : * * 156 : * Derivee -/- T : 157 : DNDT1 = APN1 ; 158 : * * 159 : *----------------- MASSE VOLUMIQUE DE L'EAU LIQUIDE ------------------* 160 : * * 161 : * ENTREES * 162 : * - Temperature (Celsius) * 163 : * SORTIES * 164 : * - Masse volumique de l'eau (kg/m3) : RHOW0 * 165 : * - Derivee -/- T (kg/m3/K) : DRHOWDT0 * 166 : * - Derivee -/- Pg ou Pc (kg/m3/Pa) : DRHOWDP1 * 167 : * PARAMETRES * 168 : * PAi ,PBi : facteurs polynomiaux : * 169 : * * 170 : PA0 = 4.8863E-7 ; 171 : PA1 = -1.6528E-9 ; 172 : PA2 = 1.8621E-12 ; 173 : PA3 = 2.4266E-13 ; 174 : PA4 = -1.5996E-15 ; 175 : PA5 = 3.3703E-18 ; 176 : PB0 = 1.0213E+3 ; 177 : PB1 = -7.7377E-1 ; 178 : PB2 = 8.7696E-3 ; 179 : PB3 = -9.2118E-5 ; 180 : PB4 = 3.3534E-7 ; 181 : PB5 = -4.4034E-10 ; 182 : PAX = (PA5 * CHTCSC0) + PA4 ; 183 : PAX = (PAX * CHTCSC0) + PA3 ; 184 : PAX = (PAX * CHTCSC0) + PA2 ; 185 : PAX = (PAX * CHTCSC0) + PA1 ; 186 : PAX = (PAX * CHTCSC0) + PA0 ; 187 : PBX = (PB5 * CHTCSC0) + PB4 ; 188 : PBX = (PBX * CHTCSC0) + PB3 ; 189 : PBX = (PBX * CHTCSC0) + PB2 ; 190 : PBX = (PBX * CHTCSC0) + PB1 ; 191 : PBX = (PBX * CHTCSC0) + PB0 ; 192 : * * 193 : * La masse volumique de l'eau liquide : Ok ! * 194 : RHOW0 = (PAX * -1.E7) + PBX ; 195 : * * 196 : * Sa derivee -/- T : Ok ! 197 : DAX = (4. * PA4) + (5. * PA5 * CHTCSC0) ; 198 : DAX = (3. * PA3) + (DAX * CHTCSC0) ; 199 : DAX = (2. * PA2) + (DAX * CHTCSC0) ; 200 : DAX = PA1 + (DAX * CHTCSC0) ; 201 : DBX = (4. * PB4) + (5. * PB5 * CHTCSC0) ; 202 : DBX = (3. * PB3) + (DBX * CHTCSC0) ; 203 : DBX = (2. * PB2) + (DBX * CHTCSC0) ; 204 : DBX = PB1 + (DBX * CHTCSC0) ; 205 : DRHOWDT0 = (DAX * -1.E7) + DBX ; 206 : * * 207 : * Sa derivee -/- Pg et Pc : * 208 : DRHOWDP1 = 0. ; 209 : * * 210 : *------------------- PRESSION DE VAPEUR SATURANTE --------------------* 211 : * * 212 : * ENTREE * 213 : * - Temperature (K) * 214 : * SORTIES * 215 : * - Pression de vapeur saturante (Pa) : CHPVSAT0 * 216 : * - Derivee -/- T (Pa/K) : DPSATDT0 * 217 : * PARAMETRES * 218 : * GCi : facteurs polynomiaux * 219 : GC1 = -5.8002206E+03 ; 220 : GC2 = 1.3914993E+00 ; 221 : GC3 = -4.8640239E-02 ; 222 : GC4 = 4.1764768E-05 ; 223 : GC5 = -1.4452093E-08 ; 224 : GC6 = 6.5459673E+00 ; 225 : * * 226 : * Pression de vapeur saturante : * 227 : CHPVSAT0 = (GC1 * CHTKSCM1) + GC2 + (GC3 * CHTKSC0) 228 : + (GC4 * (CHTKSC0 ** 2)) + (GC5 * (CHTKSC0 ** 3)) 229 : + (GC6 * ('LOG' CHTKSC0)) ; 230 : CHPVSAT0 = 'EXP' CHPVSAT0 ; 231 : * * 232 : * Derivee -/- T : * 233 : DPSATDT0 = (-1. * GC1 * (CHTKSCM1 ** 2)) + GC3 + (2. * GC4 * CHTKSC0) 234 : + (3. * GC5 * (CHTKSC0 ** 2)) + (GC6 * CHTKSCM1) ; 235 : DPSATDT0 = DPSATDT0 * CHPVSAT0 ; 236 : * * 237 : *---------------------- TAUX DE DESHYDRATATION -----------------------* 238 : * * 239 : * ENTREE * 240 : * - Temperature (Celcius) * 241 : * SORTIES * 242 : * - Taux de deshydratation (kg/s) : MDH0 * 243 : * - Derivee -/- T (kg/s/K) : DMDHDT0 * 244 : * PARAMETRES * 245 : * FS1 : facteur stoechiometrique * 246 : * PM1 : parametre prenanyt en compte l'age du beton * 247 : * GC1 : une constante * 248 : FS1 = 0.4 ; 249 : PM1 = 0.24 ; 250 : GC1 = 200. ; 251 : FTC1 = 'EXP' (-0.004 * (CHTCSC0 - 105.)) ; 252 : DFTC1 = -0.004 * FTC1 ; 253 : FTC1 = 0.5 * PI * (1. - (2. * FTC1)) ; 254 : DFTC1 = -1. * PI * DFTC1 ; 255 : DFTC1 = 0.5 * ('COS' (180. / PI * FTC1)) * DFTC1 ; 256 : FTC1 = 0.5 * (1. + ('SIN' (180. / PI * FTC1))) ; 257 : * * 258 : * Taux de deshydratation : Ok ! * 259 : MAQST1 = CHTCSC0 'MASQ' 'SUPE' 105. ; 260 : MDH0 = FS1 * PM1 * GC1 * FTC1 * MAQST1 ; 261 : * * 262 : * Derivee -/- T : Ok ! * 263 : DMDHDT0 = FS1 * PM1 * GC1 * DFTC1 * MAQST1 ; 264 : * * 265 : *--------------------- ENTHALPIE DE VAPORISATION ---------------------* 266 : * * 267 : * ENTREE * 268 : * - Temperature (K) * 269 : * SORTIES * 270 : * - Enthalpie de vaporisation (J) : HVAP0 * 271 : * - Derivee -/- T (J/K) : DHVAPDT1 * 272 : * PARAMETRES * 273 : * TCRI1 : Temperature critique (idem que saturation) * 274 : * PN1 : exposant * 275 : * KHVAPT1 : coeff. de proportionalite en temperature * 276 : TCRI1 = 647.15 ; 277 : PN1 = 0.38 ; 278 : KHVAPT1 = 2.672E5 ; 279 : DCHTK1 = TCRI1 - CHTKSC0 ; 280 : HVAP0 = KHVAPT1 * (DCHTK1 ** PN1) ; 281 : DHVAPDT1 = -1. * PN1 * HVAP0 * (DCHTK1 ** -1) ; 282 : * * 283 : *--------------------- PERMEABILITE INTRINSEQUE ----------------------* 284 : * * 285 : * ENTREE * 286 : * - Temperature (K) * 287 : * - Pression de gaz (Pa) * 288 : * SORTIES * 289 : * - Permeabilite intrinseque (---) : KINT0 * 290 : * PARAMETRES * 291 : * T0 : Temperature de reference (25 degres Celcius) * 292 : * P0 : Pression de reference (1 atm.) * 293 : * KK0 : Permeabilite a (T0,P0) * 294 : * AT0 : coef. sur temperature * 295 : * AP0 : coef. sur pression * 296 : T0 = 298.15 ; 297 : P0 = 101325. ; 298 : KK0 = 3.2E-18 ; 299 : AT0 = 0.005 ; 300 : AP0 = ('LOG' (5. / 3.)) / ('LOG' 4.) ; 301 : KINT0 = KK0 * ('EXP' (('LOG' 10.) * AT0 * (CHTKSC0 - T0))) 302 : * ((CHPGSC0 / P0) ** AP0) ; 303 : * * 304 : *------------------- PERMEABILITE RELATIVE AU GAZ --------------------* 305 : * * 306 : * ENTREE * 307 : * - Saturation (-) : SW0 * 308 : * SORTIE * 309 : * - Permeabilite relative au gaz (---) : KRG0 * 310 : * PARAMETRES * 311 : * SCR1 : Saturation critique * 312 : * AG1 : exposant * 313 : SCR1 = 1. ; 314 : AG1 = 1. ; 315 : KRG0 = 1. - ((SW0 / SCR1) ** AG1) ; 316 : * * 317 : *------------------- PERMEABILITE RELATIVE A L'EAU -------------------* 318 : * * 319 : * ENTREE * 320 : * - Saturation (-) : SW0 * 321 : * SORTIE * 322 : * - Permeabilite relative a l'eau (---) : KRW0 * 323 : * PARAMETRES * 324 : * AW1 : exposant * 325 : AW1 = 6 ; 326 : KRW0 = SW0 ** AW1 ; 327 : * * 328 : *------------------------- VISCOSITE DE L'EAU ------------------------* 329 : * * 330 : * ENTREE * 331 : * - Temperature (K) * 332 : * SORTIE * 333 : * - Viscosite de l'eau (---) : MUW0 * 334 : * PARAMETRES * 335 : * T0 : Temperature de reference (K) * 336 : * ALPHW1 : coef. sur temperature * 337 : * EXPOW1 : exposant 338 : T0 = 229. ; 339 : ALPHGW1 = 0.6612 ; 340 : EXPOW1 = -1.562 ; 341 : MUW0 = ALPHGW1 * ((CHTKSC0 - T0) ** EXPOW1) ; 342 : * * 343 : *------------------- VISCOSITE DE LA VAPEUR D'EAU --------------------* 344 : * * 345 : * ENTREE * 346 : * - Temperature (C) * 347 : * SORTIE * 348 : * - Viscosite de la vapeur d'eau (---) : MUGW0 * 349 : * PARAMETRES * 350 : * MUGW1 : Viscosite a 0 degre Celcius * 351 : * ALPHGW1 : coef. sur temperature * 352 : MUGW1 = 8.85E-6 ; 353 : ALPHGW1 = 3.633E-8 ; 354 : MUGW0 = MUGW1 + (ALPHGW1 * CHTCSC0) ; 355 : * * 356 : *---------------------- VISCOSITE DE L'AIR SEC -----------------------* 357 : * * 358 : * ENTREE * 359 : * - Temperature (C) * 360 : * SORTIE * 361 : * - Viscosite de l'air sec (---) : MUGA0 * 362 : * PARAMETRES * 363 : * MUGA1 : Viscosite a 0 degre Celcius * 364 : * ALPHGA1 : coef. sur temperature * 365 : * ALPHGA2 : coef. sur temperature au carre * 366 : MUGA1 = 17.17E-6 ; 367 : ALPHGA1 = 4.733E-8 ; 368 : ALPHGA2 = -2.222E-11 ; 369 : MUGA0 = MUGA1 + (ALPHGA1 * CHTCSC0) + (ALPHGA2 * (CHTCSC0 ** 2)) ; 370 : * * 371 : *-------------- DIFFUSITE DE LA VAPEUR D'EAU DANS l'AIR --------------* 372 : * * 373 : * ENTREES * 374 : * - Temperature (K) * 375 : * - Pression gazeuse (Pa) * 376 : * SORTIE * 377 : * - Diffusivite de la vapeur d'eau dans l'air (---) : DGWGA0 * 378 : * PARAMETRES * 379 : * T0 : Temperature de reference (0 degre C) * 380 : * P0 : Pression e reference (1 atm.) * 381 : * DGWGA1 : Diffusite a (T0,P0) * 382 : * ET1 : Exposant sur la temperature * 383 : T0 = 273.15 ; 384 : P0 = 101325. ; 385 : DGWGA1 = 2.17E-5 ; 386 : ET1 = 1.88 ; 387 : DGWGA0 = DGWGA1 * (CHPGSCM1 * P0) * ((CHTKSC0 / T0) ** ET1) ; 388 : 389 : *------------------------ DIFFUSITE EFFECTIVE ------------------------* 390 : * * 391 : * ENTREES * 392 : * - Saturation (-) : SW0 * 393 : * - Porosite (-) : N0 * 394 : * - Diffusivite de la vapeur d'eau dans l'air (---) : DGWGA0 * 395 : * SORTIE * 396 : * - Diffusivite effective (---) : DEFF0 * 397 : UNTIERS1 = 1. / 3. ; 398 : TORT0 = (N0 ** UNTIERS1) * (SG0 ** (7. * UNTIERS1)) ; 399 : DEFF0 = SG0 * N0 * TORT0 * DGWGA0 ; 400 : * * 401 : *----------------- CONDUCTIVITE THERMIQUE DU SOLIDE ------------------* 402 : * * 403 : * ENTREES * 404 : * - Temperature (K) * 405 : * SORTIE * 406 : * - Conductivite thermique du solide (W/m/K) : KS0 * 407 : * PARAMETRES * 408 : * T0 : Temperature de reference (25 degres C) * 409 : * KS1 : Conductivite thermique a T0 * 410 : * AKS1 : Coeff. de proportionalite en temperature * 411 : T0 = 298. ; 412 : KS1 = 2.1 ; 413 : AKS1 = -1.017E-3 ; 414 : KS0 = KS1 + (AKS1 * (CHTKSC0 - T0)) ; 415 : * * 416 : *---------------------------------------------------------------------* 417 : * * 418 : * CALCUL DES TERMES DE LA MATRICE DU SYSTEME * 419 : * * 420 : *---------------------------------------------------------------------* 421 : * * 422 : *---------------------- Quantites intermediaires ---------------------* 423 : * * 424 : * CHPGWSC0 : Pression de vapeur d'eau (Pa, issue formule de Kelvin) : * 425 : CHXX = -1. * FGPGW0 * (RHOW0 ** -1) ; 426 : CHYY = 'EXP' (CHXX * CHPCSC0) ; 427 : CHPGWSC0 = CHPVSAT0 * CHYY ; 428 : * DPGWDPC0 : Derivee -/- Pc (-) : * 429 : DPGWDPC0 = CHXX * CHPGWSC0 ; 430 : * DPGWDT0 : Derivee -/- T (Pa/K) : * 431 : DPGWDT0 = (DPSATDT0 * CHYY) + (CHPVSAT0 * CHYY * CHPCSC0 432 : * ( (FGPGW0 * (RHOW0 ** -1) * CHTKSCM1) 433 : + (FGPGW0 * (RHOW0 ** -2) * DRHOWDT0))) ; 434 : * * 435 : * MUG0 : Viscosite de la phase gazeuse : * 436 : MASQ1 = (CHPGSC0 - CHPGWSC0) 'MASQ' 'SUPE' 0. ; 437 : CHXX = (1. - (CHPGWSC0 * CHPGSCM1)) * MASQ1 ; 438 : MUG0 = MUGW0 + ((MUGA0 - MUGW0) * (CHXX ** 0.6083)) ; 439 : * * 440 : * RHOGA0 : Masse volumique de l'air sec (kg/m3) : Ok ! * 441 : RHOGA0 = (CHPGSC0 - CHPGWSC0) * FGPGA0 ; 442 : * DRHOGADT : Derivee -/- T (kg/m3/K) : Ok ! * 443 : DRHOGADT = -1. * ((RHOGA0 * CHTKSCM1) + (FGPGA0 * DPGWDT0)) ; 444 : * * 445 : * RHOGW0 : Masse volumique de la vapeur (kg/m3) : * 446 : RHOGW0 = FGPGW0 * CHPGWSC0 ; 447 : * DRHOGWDT : Derivee -/- T (kg/m3/K) : Ok ! * 448 : DRHOGWDT = (FGPGW0 * DPGWDT0) - (RHOGW0 * CHTKSCM1) ; 449 : * * 450 : * RHOG0 : Masse volumique de la phase gazeuse (kg/m3) : Ok ! * 451 : RHOG0 = (FGPGA0 * CHPGSC0) + ((FGPGW0 - FGPGA0) * CHPGWSC0) ; 452 : * * 453 : * CPG0 : Capacite calorifique a P constante du gaz (J/kg/K) Ok ! * 454 : CPG0 = (RHOG0 * CPGA1) + (RHOGW0 * (CPGW1 - CPGA1)) ; 455 : * * 456 : * Masse Molaire du gaz (kg/mol) : * 457 : MMG0 = MMGA1 + ((MMGW1 - MMGA1) * (CHPGWSC0 * CHPGSCM1)) ; 458 : * * 459 : *------------------------ Matrice de Capacite ------------------------* 460 : * * 461 : * Pression de Gaz (PG) : * 462 : CGG0 = N0 * SG0 * FGPGA0 ; 463 : * * 464 : CGC0 = (RHOGA0 * DSWDPC0) + (SG0 * FGPGA0 * DPGWDPC0) ; 465 : CGC0 = -1. * N0 * CGC0 ; 466 : * * 467 : 468 : CGT0 = (N0 * ((SG0 * DRHOGADT) - (RHOGA0 * DSWDT0))) 469 : + (DMDHDT0 * RHOGA0 / RHOS1 * SG0) ; 470 : * * 471 : CGU0 = SG0 * RHOGA0 ; 472 : * * 473 : * Pression Capillaire (PC) : * 474 : CCG0 = 0. ; 475 : * * 476 : DRHOW0 = RHOW0 - RHOGW0 ; 477 : CCC0 = N0 * ((SG0 * FGPGW0 * DPGWDPC0) + (DRHOW0 * DSWDPC0)) ; 478 : * * 479 : CCU0 = (RHOGW0 * SG0) + (RHOW0 * SW0) ; 480 : * * 481 : * Remarque : signe termes avec DMDHDT0 opposes dans eq. 37 du rapport * 482 : * SdP et premiere integrale de l'eq. ne se retrouve pas ci-dessous. * 483 : CCT0 = (SG0 * DRHOGWDT) + (SW0 * DRHOWDT0) + (DRHOW0 * DSWDT0) ; 484 : CCT0 = (N0 * CCT0) - DMDHDT0 + (DMDHDT0 * CCU0 / RHOS1) ; 485 : * * 486 : * Temperature : * 487 : CTG0 = 0. ; 488 : * * 489 : CTC0 = -1. * HVAP0 * N0 * RHOW0 * DSWDPC0 ; 490 : * * 491 : * Remarque : corrections dans rapport SdP. * 492 : CTT0 = ((1. - N0) * RHOS1 * CPS0) 493 : + (N0 * ((SW0 * RHOW0 * CPW1) + (SG0 * CPG0))) 494 : + (DMDHDT0 * DHHD1) ; 495 : CTT0 = CTT0 - (HVAP0 * (N0 * ((RHOW0 * DSWDT0) + (SW0 * DRHOWDT0)) 496 : + (RHOW0 / RHOS1 * SW0 * DMDHDT0) - DMDHDT0)) ; 497 : * * 498 : CTU0 = -1. * HVAP0 * RHOW0 * SW0 ; 499 : * * 500 : *'MESS' 'CTC0 =' ('MINI' CTC0) ; 501 : *'MESS' 'CTT0 =' ('MINI' CTT0) ; 502 : *'MESS' 'CTU0 =' ('MINI' CTU0) ; 503 : *'MESS' 'CGG0 =' ('MINI' CGG0) ; 504 : *'MESS' 'CGC0 =' ('MINI' CGC0) ; 505 : *'MESS' 'CGT0 =' ('MINI' CGT0) ; 506 : *'MESS' 'CGU0 =' ('MINI' CGU0) ; 507 : *'MESS' 'CCC0 =' ('MINI' CCC0) ; 508 : *'MESS' 'CCT0 =' ('MINI' CCT0) ; 509 : *'MESS' 'CCU0 =' ('MINI' CCU0) ; 510 : * * 511 : CGG0 = 'CHAN' 'CHAM' CGG0 MOD1 ; 512 : CGC0 = 'CHAN' 'CHAM' CGC0 MOD1 ; 513 : CGT0 = 'CHAN' 'CHAM' CGT0 MOD1 ; 514 : CCC0 = 'CHAN' 'CHAM' CCC0 MOD1 ; 515 : CCT0 = 'CHAN' 'CHAM' CCT0 MOD1 ; 516 : CTC0 = 'CHAN' 'CHAM' CTC0 MOD1 ; 517 : CTT0 = 'CHAN' 'CHAM' CTT0 MOD1 ; 518 : * * 519 : *---------------------- Matrice de conductivite ----------------------* 520 : * * 521 : * Pression de Gaz (PG) : * 522 : RPERMGA0 = RHOGA0 * KINT0 * KRG0 * (MUG0 ** -1) ; 523 : CHFDIFF0 = RHOG0 * MMGA1 * MMGW1 * (MMG0 ** -2) * DEFF0 ; 524 : KGG0 = RPERMGA0 + (CHFDIFF0 * CHPGWSC0 * CHPGSCM2) ; 525 : * * 526 : SGN1 = 1. ; 527 : * Remarque : signe KGC0 oppose dans eq. 31 rapport SdP. * 528 : KGC0 = SGN1 * CHFDIFF0 * DPGWDPC0 * CHPGSCM1 ; 529 : * * 530 : * Remarque : signe KGT0 oppose dans eq. 32 rapport SdP. * 531 : KGT0 = SGN1 * CHFDIFF0 * DPGWDT0 * CHPGSCM1 ; 532 : * * 533 : KGU0 = 0. ; 534 : * * 535 : * Pression Capillaire (PC) : * 536 : RPERMG0 = KINT0 * KRG0 * (MUG0 ** -1) ; 537 : RPERMGW0 = RHOGW0 * RPERMG0 ; 538 : RPERMG0 = RHOG0 * RPERMG0 ; 539 : RPERMW0 = RHOW0 * KINT0 * KRW0 * (MUW0 ** -1) ; 540 : * * 541 : * Remarque : 1er terme de signe oppose dans eq. 39 rapport SdP, ainsi * 542 : * que terme RPERMGW0. 543 : KCG0 = (CHFDIFF0 * CHPGWSC0 * CHPGSCM2) - RPERMGW0 + RPERMW0 ; 544 : * * 545 : * Remarque : Dans rapport SdP, KCC0 de signe oppose. * 546 : KCC0 = (-1. * SGN1 * KGC0) - RPERMW0 ; 547 : * * 548 : * Remarque : idem KCC0 pour KCT0. * 549 : KCT0 = -1. * SGN1 * KGT0 ; 550 : * * 551 : KCU0 = 0. ; 552 : * * 553 : * Temperature : * 554 : * ATTENTION : termes en GRAD PG/PC/T ! 555 : KTG0 = -1. * HVAP0 * RPERMW0 ; 556 : KTGGT0 = -1. * ((CPW1 * RPERMW0) + (CPG0 * RPERMG0)) ; 557 : *KG21 = KTGGT0 * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'T,X' 'SCAL')) ; 558 : *KG22 = KTGGT0 * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'T,Y' 'SCAL')) ; 559 : * * 560 : KTC0 = -1. * KTG0 ; 561 : KTCGT0 = CPW1 * RPERMW0 ; 562 : *KC21 = KTCGT0 563 : * * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'T,X' 'SCAL')) ; 564 : *KC22 = KTCGT0 565 : * * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'T,Y' 'SCAL')) ; 566 : 567 : * * 568 : * Remarque : gravite negligee. * 569 : KTT0 = KS0 570 : * ((4. * N0 * SW0 * RHOW0 * ((1. - N0) ** -1) / RHOS1) + 1.) ; 571 : KTTGPG0 = KTGGT0 ; 572 : KTTGPC0 = KTCGT0 ; 573 : *KT21 = (KTTGPG0 * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'PG,X' 'SCAL'))) 574 : * + (KTTGPC0 * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'PC,X' 'SCAL'))) ; 575 : *KT22 = (KTTGPG0 * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'PG,Y' 'SCAL'))) 576 : * + (KTTGPC0 * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'PC,Y' 'SCAL'))) ; 577 : * * 578 : KTU0 = 0. ; 579 : * * 580 : *'MESS' 'KGG0 =' ('MINI' KGG0) ; 581 : *'MESS' 'KGC0 =' ('MINI' KGC0) ; 582 : *'MESS' 'KGT0 =' ('MINI' KGT0) ; 583 : *'MESS' 'KCG0 =' ('MINI' KCG0) ; 584 : *'MESS' 'KCC0 =' ('MINI' KCC0) ; 585 : *'MESS' 'KCT0 =' ('MINI' KCT0) ; 586 : *'MESS' 'KTG0 =' ('MINI' KTG0) ; 587 : *'MESS' 'KTC0 =' ('MINI' KTC0) ; 588 : *'MESS' 'KTT0 =' ('MINI' KTT0) ; 589 : * ATTENTION : valeurs des termes gradient differentes de ZH car 590 : * ces dernieres sont fausses (capacite cal. pas en USI) * 591 : *'MESS' 'KTGGT0 =' ('MINI' KTGGT0) ; 592 : *'MESS' 'KTCGT0 =' ('MINI' KTCGT0) ; 593 : *'MESS' 'KTTGPG0 =' ('MINI' KTTGPG0) ; 594 : *'MESS' 'KTTGPC0 =' ('MINI' KTTGPC0) ; 595 : * * 596 : * Sauvegarde CHPO pour 2nd membre : * 597 : KGC1 = KGC0 ; 598 : KGT1 = KGT0 ; 599 : KCG1 = KCG0 ; 600 : KCT1 = KCT0 ; 601 : KTC1 = KTC0 ; 602 : * * 603 : KGG0 = 'CHAN' 'CHAM' KGG0 MOD1 ; 604 : KGC0 = 'CHAN' 'CHAM' KGC0 MOD1 ; 605 : KGT0 = 'CHAN' 'CHAM' KGT0 MOD1 ; 606 : KCG0 = 'CHAN' 'CHAM' KCG0 MOD1 ; 607 : KCC0 = 'CHAN' 'CHAM' KCC0 MOD1 ; 608 : KCT0 = 'CHAN' 'CHAM' KCT0 MOD1 ; 609 : KTG0 = 'CHAN' 'CHAM' KTG0 MOD1 ; 610 : KTC0 = 'CHAN' 'CHAM' KTC0 MOD1 ; 611 : KTT0 = 'CHAN' 'CHAM' KTT0 MOD1 ; 612 : *KG21 = 'CHAN' 'CHAM' KG21 MOD1 ; 613 : *KG22 = 'CHAN' 'CHAM' KG22 MOD1 ; 614 : *KC21 = 'CHAN' 'CHAM' KC21 MOD1 ; 615 : *KC22 = 'CHAN' 'CHAM' KC22 MOD1 ; 616 : *KT21 = 'CHAN' 'CHAM' KT21 MOD1 ; 617 : *KT22 = 'CHAN' 'CHAM' KT22 MOD1 ; 618 : KTGGT0 = 'CHAN' 'CHAM' KTGGT0 MOD1 ; 619 : KTCGT0 = 'CHAN' 'CHAM' KTCGT0 MOD1 ; 620 : *KTTGPG0 = 'CHAN' 'CHAM' KTTGPG0 MOD1 ; 621 : *KTTGPC0 = 'CHAN' 'CHAM' KTTGPC0 MOD1 ; 622 : * * 623 : *----- Construction du MCHAML de Caracteristiques : appel a MATE -----* 624 : * * 625 : MAT1 = 'MATE' MOD1 626 : 'CGG' CGG0 'CGC' CGC0 'CGT' CGT0 627 : 'CCG' CCG0 'CCC' CCC0 'CCT' CCT0 628 : 'CTG' CTG0 'CTC' CTC0 'CTT' CTT0 629 : 'KGG' KGG0 'KGC' KGC0 'KGT' KGT0 630 : 'KCG' KCG0 'KCC' KCC0 'KCT' KCT0 631 : 'KTG' KTG0 'KTC' KTC0 'KTT' KTT0 'KTGG' KTGGT0 'KTCG' KTCGT0 ; 632 : * * 633 : 'RESP' MAT1 ; 634 : * * 635 : *---------------------------------------------------------------------* 636 : * * 637 : * CALCUL DES TERMES DU SECOND MEMBRE * 638 : * * 639 : *---------------------------------------------------------------------* 640 : * * 641 : * Acceleration de la gravite : * 642 : AGRAV1 = 9.81 ; 643 : * * 644 : * Composante "gravite" a preciser car depend du produit scalaire avec * 645 : * gradient de la fonction d'interpolation support inconnue Pg. * 646 : *FG0 = RHOG0 * RPERMGA0 * AGRAV1 ; 647 : *FG0 = FG0 'NOMC' 'Q' ; 648 : * * 649 : * Composante "gravite" a preciser car depend du produit scalaire avec * 650 : * gradient de la fonction d'interpolation support inconnue Pc. * 651 : *FC0 = ((RPERMGW0 * RHOG0) + (RPERMW0 * RHOW0)) * AGRAV1 ; 652 : *FC0 = FC0 'NOMC' 'Q' ; 653 : * * 654 : * Composante "gravite" a preciser car depend du produit scalaire avec * 655 : * gradient de la fonction d'interpolation support inconnue T. * 656 : *FT0 = ((CPW1 * RPERMW0 * RHOW0) + (CPG0 * RPERMG0 * RHOG0)) 657 : * * ('CHAN' 'CHPO' MOD1 ('EXCO' GRCH0 'T,Y' 'SCAL')) * AGRAV1 ; 658 : *FT0 = (HVAP0 * RPERMW0 * RHOW0 * AGRAV1) + FT0 ; 659 : *FT0 = FT0 'NOMC' 'Q' ; 660 : * * 661 : * Laplaciens de Pg, Pc, T : * 662 : GRCH1 = 'CHAN' 'CHPO' MOD1 GRCH0 ; 663 : GRXPG0 = 'EXCO' GRCH1 ('EXTR' LMGPG1 1) 'PG' ; 664 : GRXPC0 = 'EXCO' GRCH1 ('EXTR' LMGPC1 1) 'PC' ; 665 : GRXTK0 = 'EXCO' GRCH1 ('EXTR' LMGTK1 1) 'T' ; 666 : GRYPG0 = 'EXCO' GRCH1 ('EXTR' LMGPG1 2) 'PG' ; 667 : GRYPC0 = 'EXCO' GRCH1 ('EXTR' LMGPC1 2) 'PC' ; 668 : GRYTK0 = 'EXCO' GRCH1 ('EXTR' LMGTK1 2) 'T' ; 669 : IDIM3 = 'EGA' VDIM1 3 ; 670 : 'SI' IDIM3 ; 671 : GRZPG0 = 'EXCO' GRCH1 ('EXTR' LMGPG1 3) 'PG' ; 672 : GRZPC0 = 'EXCO' GRCH1 ('EXTR' LMGPC1 3) 'PC' ; 673 : GRZTK0 = 'EXCO' GRCH1 ('EXTR' LMGTK1 3) 'T' ; 674 : 'FINS' ; 675 : * * 676 : FG0X = GRXPG0 + (KGC1 * GRXPC0) + (KGT1 * GRXTK0) ; 677 : FG0Y = GRYPG0 + (KGC1 * GRYPC0) + (KGT1 * GRYTK0) ; 678 : FG0X = 'GRAD' MOD1 FG0X ; 679 : FG0Y = 'GRAD' MOD1 FG0Y ; 680 : CPGX1 = 'EXTR' LMGPG1 1 ; 681 : CPGY1 = 'EXTR' LMGPG1 2 ; 682 : CPCX1 = 'EXTR' LMGPC1 1 ; 683 : CPCY1 = 'EXTR' LMGPC1 2 ; 684 : CTKX1 = 'EXTR' LMGTK1 1 ; 685 : CTKY1 = 'EXTR' LMGTK1 2 ; 686 : FG0G = ('EXCO' FG0X CPCX1 'SCAL') + ('EXCO' FG0X CTKX1 'SCAL') 687 : + ('EXCO' FG0Y CPCY1 'SCAL') + ('EXCO' FG0Y CTKY1 'SCAL') ; 688 : 'SI' IDIM3 ; 689 : FG0Z = GRZPG0 + (KGC1 * GRZPC0) + (KGT1 * GRZTK0) ; 690 : FG0Z = 'GRAD' MOD1 FG0Z ; 691 : FG0G = ('EXCO' FG0Z 'PC,Z' 'SCAL') + ('EXCO' FG0Z 'T,Z' 'SCAL') 692 : + FG0G ; 693 : 'FINS' ; 694 : FG0 = 'SOUR' MOD1 FG0G 'QG' ; 695 : * * 696 : FC0X = (KCG1 * GRXPG0) + GRXPC0 + (KCT1 * GRXTK0) ; 697 : FC0Y = (KCG1 * GRYPG0) + GRYPC0 + (KCT1 * GRYTK0) ; 698 : FC0X = 'GRAD' MOD1 FC0X ; 699 : FC0Y = 'GRAD' MOD1 FC0Y ; 700 : FC0G = ('EXCO' FC0X CPGX1 'SCAL') + ('EXCO' FC0X CTKX1 'SCAL') 701 : + ('EXCO' FC0Y CPGY1 'SCAL') + ('EXCO' FC0Y CTKY1 'SCAL') ; 702 : 'SI' IDIM3 ; 703 : FC0Z = (KCG1 * GRZPG0) + GRZPC0 + (KCT1 * GRZTK0) ; 704 : FC0Z = 'GRAD' MOD1 FC0Z ; 705 : FC0G = ('EXCO' FC0Z 'PG,Z' 'SCAL') + ('EXCO' FC0Z 'T,Z' 'SCAL') 706 : + FC0G ; 707 : 'FINS' ; 708 : FC0 = 'SOUR' MOD1 FC0G 'QC' ; 709 : * * 710 : FT0X = (KTC1 * (GRXPG0 + GRXPC0)) + GRXTK0 ; 711 : FT0Y = (KTC1 * (GRYPG0 + GRYPC0)) + GRYTK0 ; 712 : FT0X = 'GRAD' MOD1 FT0X ; 713 : FT0Y = 'GRAD' MOD1 FT0Y ; 714 : FT0G = ('EXCO' FT0X CPGX1 'SCAL') + ('EXCO' FT0X CPCX1 'SCAL') 715 : + ('EXCO' FT0Y CPGY1 'SCAL') + ('EXCO' FT0Y CPCY1 'SCAL') ; 716 : 'SI' IDIM3 ; 717 : FT0Z = (KTC1 * (GRZPG0 + GRZPC0)) + GRZTK0 ; 718 : FT0Z = 'GRAD' MOD1 FT0Z ; 719 : FT0G = ('EXCO' FT0Z 'PG,Z' 'SCAL') + ('EXCO' FT0Z 'PC,Z' 'SCAL') 720 : + FT0G ; 721 : 'FINS' ; 722 : FT0 = 'SOUR' MOD1 FT0G 'Q' ; 723 : * * 724 : F0 = FG0 + FC0 + FT0 ; 725 : * * 726 : *'MESS' 'FG0 =' ('MINI' FG0) ; 727 : *'MESS' 'FC0 =' ('MINI' FC0) ; 728 : *'MESS' 'FT0 =' ('MINI' FT0) ; 729 : * * 730 : 'RESP' F0 ; 731 : * * 732 : *---------------------------------------------------------------------* 733 : * * 734 : 'FINP' ; 735 : *---------------------------------------------------------------------* 736 : * FIN PROCEDURE @MATETHM 737 : 738 :
© Cast3M 2003 - All rights reserved.
Disclaimer