Télécharger la plaquette Cast3M


fr en

Télécharger tstc3d.eso

Retour à la liste

Numérotation des lignes : 

tstc3d

  1. C TSTC3D    SOURCE    FD218221  24/02/07    21:15:29     11834          
  2.       subroutine TSTC3D(poro0,EPLCCV,deplcc3,precision3d,
  3.      #log_reduc,reduc)
  4.  
  5. c     test refermetures excessives de la porosite lors de l'ecoulement 
  6. c     et calcul du coeff  de reduction de retour radial pour camclay
  7.  
  8.       implicit real*8 (a-h,o-z)
  9.       implicit integer (i-n)
  10.  
  11.       real*8 poro0,EPLCCV,deplcc3(3),precision3d,reduc
  12.       logical log_reduc
  13.  
  14.       real*8 dporo,poro1
  15.  
  16. c       test consolidation maxi non dépassée
  17.         dporo=deplcc3(1)+deplcc3(2)+deplcc3(3)
  18.         poro1=poro0+EPLCCV+dporo
  19. c        print*,'ds tstc3d',poro0,poro1,EPLCCV,dporo
  20. c        read*
  21.         if(poro1.lt.(poro0*precision3d)) then
  22. c           on arrive a la limite de la consolidation possible  
  23. c           on reduit le retour pour ne pas avoir de poro nulle
  24.             if(abs(dporo).gt.(precision3d*poro0)) then
  25.                reduc=-(poro0+EPLCCV)/dporo
  26.             else
  27.                reduc=-(poro0+EPLCCV)/(precision3d*poro0)
  28.             end if
  29.             log_reduc=.true.
  30.             print*,'Testcc3d reduction  increment plastique CC',reduc
  31. c            print*,poro0,EPLCCV,deplcc3,precision3d,
  32. c     #      log_reduc,reduc
  33. c            read*
  34.         else
  35.             reduc=1.d0        
  36.         end if
  37.         return
  38.         end
  39.  
  40.  
  41.  

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