Accueil | Cast3M

   1 : $$$$ K_PRO    NOTICE  CHAT      11/09/12    21:16:47     7124           
   2 :                                              DATE     11/09/12
   3 : 
   4 :   Procedure K_PRO                         Voir aussi : DARCYSAT 
   5 :      ---------------                                      HT_PRO  
   6 :         K1 = K_PRO TAB1 TAB2 ;                   
   7 : 
   8 : 
   9 :      Objet :           
  10 :      _______       
  11 :            
  12 :      Cette procedure permet de calculer la permeabilite a l'eau d'un 
  13 :      milieu poreux non sature
  14 :                                        
  15 : 
  16 :     
  17 : 
  18 :          K1 :  CHPOINT contenant la permeabilite
  19 : 
  20 :          TAB1 : table contenant les variables du problemes
  21 : 
  22 :          TAB2 : table contenant  les parametres de la loi de permeabilite
  23 :  
  24 :     Commentaires :                                                             
  25 :      ____________                                                             
  26 : 
  27 :    Les indices de la table TAB1 sont les suivants
  28 : 
  29 :                 .'PNS_PROV' : CHPOINT des pressions negatives  
  30 :                 .'TH2O_PROV' : CHPOINT des teneurs en eau 
  31 :                 .'SATURATION_PROV' : CHPOINT des taux de saturation 
  32 :                                                 
  33 :    Dans le cadre d'un appel par DARCYSAT, ces CHPOINTs sont fournies.
  34 :    L'utilisateur peut ainsi creer une procedure personnelle decrivant
  35 :    une loi K(pression), K(teneur en eau) ou  K(taux de saturation).
  36 : 
  37 :    Dans le cadre d'une utilisation hors contexte DARCYSAT, l'utilisateur 
  38 :    doit seulement fournir la saturation pour la version standard.
  39 : 
  40 :    La loi de permeabilite standard est de la forme : 
  41 :                                        K = K_sat * (S **  alpha)         
  42 :    avec K_sat, permeabilite a saturation et S, taux de saturation
  43 : 
  44 :    Les indices de la table TAB2 sont les suivants :
  45 : 
  46 :            .'ALPHA' (type FLOTTANT OR CHPOINT) exposant de la loi
  47 :            .'PERMSAT' (type FLOTTANT OR CHPOINT) permeabilite a saturation
  48 :  
  49 :     Remarques :       
  50 :      ____________   
  51 : 
  52 :   Dans le cas de l'utilisation d'une procedure K_PRO personnelle,
  53 :   TAB2 doit etre de soustype PERSONNELLE. L'utilisateur a alors toute 
  54 :   liberte de choix sur les indices de la table TAB2.  
  55 : 
  56 :   Les parametres sont de type FLOTTANT ou CHPOINT issu d'une operations
  57 :   KCHT ou KOPS. Pour un CHPOINT et dans le cadre d'une utilisation
  58 :   via la procedure DARCYSAT, le support geometrique des CHPOINT doit
  59 :   etre choisi en fonction de l'option d'homogeneisation (indice
  60 :   'HOMOGENEISATION' de la table transmise a DARCYSAT) : maillage
  61 :   'FACE' pour l'option decentree et maillage 'CENTRE' pour l'option 
  62 :    centree.
  63 : 
  64 :      Exemple :                                                             
  65 :      _______                                                             
  66 : 
  67 :    Procedure personnelle K_PRO rassemblant 2 lois de forme analytique 
  68 :    differentes s'appliquant sur deux zones distinctes : 
  69 : 
  70 : ********************************************
  71 : 'DEBPROC' KR_PRO TAB1*'TABLE' PRECED*'TABLE' ;
  72 :   si (EGA PRECED.'NOMZONE' 'SITE') ;  
  73 :     K1 = 'KOPS' TAB1.'SATURATION_PROV' '**' PRECED.'ALPHA' ;
  74 :     K1 = 'KOPS' K1 '*' PRECED.'PERMSAT' ;
  75 :   finsi ;
  76 :   si (EGA PRECED.'NOMZONE' 'BO') ;  
  77 :     K1 = (kops TAB1.'PNS_PROV' '*' PRECED.'ALPHA') EXP  ; 
  78 :     K1 = kops K1 '*' PRECED.'PERMSAT' ;
  79 :   finsi ;
  80 : 'FINPROC' K1 ;
  81 : ********************************************
  82 : 
  83 :   Cette procedure est exploitable par DARCYSAT si les parametres
  84 :   des deux lois sont donnes a la table SATUR  argument de la procedure 
  85 :   DARCYSAT de la maniere suivante (cf. notice DARCYSAT) :
  86 : 
  87 : *---- definition de la loi de permeabilite
  88 : SATUR.'LOI_PERMEABILITE' = TABLE 'MULTIZONE' ;
  89 : *- pour la zone site
  90 : SATUR.'LOI_PERMEABILITE'. 'SITE' = TABLE 'STANDARD' ;
  91 : SATUR.'LOI_PERMEABILITE'. 'SITE'. 'ALPHA' = ... ;
  92 : SATUR.'LOI_PERMEABILITE'. 'SITE'. 'PERMSAT'= ... ;
  93 : SATUR.'LOI_PERMEABILITE'. 'SITE'. 'MODELE' = ... ;
  94 : *- pour la zone BO
  95 : SATUR.'LOI_PERMEABILITE'. 'BO' = TABLE 'PERSONNELLE' ;
  96 : SATUR.'LOI_PERMEABILITE'. 'BO'. 'ALPHA' = ... ;
  97 : SATUR.'LOI_PERMEABILITE'. 'BO'. 'PERMSAT'= ... ;
  98 : SATUR.'LOI_PERMEABILITE'. 'BO'. 'MODELE' = ... ;
  99 : ... 
 100 : ...
 101 : *----calcul
 102 : DARCYSAT SATUR
 103 :