Télécharger la plaquette Cast3M


fr en

Télécharger epsw3d.eso

Retour à la liste

Numérotation des lignes : 

epsw3d

  1. C EPSW3D    SOURCE    PV090527  23/01/27    21:15:33     11574          
  2.       subroutine epsw3d(PORO,SRW,MW,M,K,G,EPSW,SIGW,ERR1)
  3.  
  4. c      sous programme de calcul de la deformation due a la capillarite
  5.  
  6. c   ********************************************************************
  7.        implicit real*8 (a-h,o-z)
  8.        implicit integer (i-n)
  9. c   ********************************************************************
  10.  
  11.       real*8 PORO,SRW,MW,M,PCW,BIOT,SIGW,EPSW
  12.  
  13. c     coeffs de compressibilite ; cisaillement, dilatation thermique
  14.       real*8 G,K
  15.       integer err1
  16.  
  17. c     variable locale      
  18.       real*8 precision3d 
  19.       parameter (precision3d=1.d-6)
  20.  
  21. c   **pression capillaire suivant Van Genuchten ************************
  22.       if((Srw.gt.0.).and.(Srw.le.1.)) then
  23.          t2 = Srw ** (0.1D1 / m)
  24.          t6 = (0.1D1 / t2 - 0.1D1) ** (0.1D1 - m)
  25.          Pcw = t6 * Mw
  26.          call bitm3d(K,G,poro,biot,err1,precision3d)
  27. c        contribution de la contrainte hydrique -biot*sw*pw si pg=0           
  28.          sigw=biot*Srw*Pcw
  29. c        deformation induite dans la matrice solide       
  30.          epsw=-sigw/(3.d0*K)
  31. C          print*,'dans epsw3d:'
  32. C          print*,'sigw=',sigw
  33. C          print*,'epsw=',epsw
  34.       else
  35.          if((Srw.le.0.).or.(Srw.gt.1.)) then
  36.             print*,'Pb dans epsw3d, valeur Srw incoherente:',Srw
  37.             err1=1
  38.             return
  39.          else
  40.             epsw=0.d0
  41.             sigw=0.d0
  42.          end if
  43.       end if
  44.  
  45.       return
  46.       end
  47.  
  48.  
  49.  
  50.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales