Télécharger la plaquette Cast3M


fr en

Télécharger hyxm3d.eso

Retour à la liste

Numérotation des lignes : 

hyxm3d

  1. C HYXM3D    SOURCE    PV090527  23/01/27    21:15:39     11574          
  2. c***********************************************************************
  3.       subroutine hyxm3d(xmat0,xmat00,xmat1,hydra1,hydras,hydrar,
  4.      # non_lin,erreur)
  5. c     effet de l hydratation sur les caracteristiques materiau
  6. c     declaration externe
  7.       real*8 xmat0,xmat1,hydra1,hydras,hydrar,xmat00,non_lin
  8.       integer erreur
  9. c     declaration locale      
  10.       real*8 yy,zz,zzmin,hydra
  11. c     avant le seuil d hydratation on a 1e-3 des caracteristiques       
  12.       parameter (zzmin=0.d0)
  13. c     initialisation controle d erreur      
  14. c      erreur=0
  15.       if((hydras.lt.0.d0).or.(hydras.gt.1.d0).
  16.      #       or.(hydrar.le.hydras)) then
  17.              print*,'------->Donnees incoherentes dans hyxm3d'
  18.              print*,'0<HYDR<1  0<HYDS<1 HREF>HYDS cf. hydra3d.eso'
  19.              print*,xmat0,xmat00,xmat1,hydra1,hydras,hydrar,non_lin
  20. c             read*
  21.              xmat1=xmat0
  23.              return
  24.       end if      
  25.       hydra=min(hydra1,1.d0)
  26.       hydra=max(hydra,0.d0)
  27.       if(hydra.ne.hydrar)then
  28.          if (hydra.gt.hydras) then
  29. c            intrpoation de DeSchutter 
  30.              if(hydrar-hydras.le.0.) then
  31.                  print*,'Pb dans hydrxamt3d de fluendo3'
  32.                  print*,'HREF=HYDS impossible'
  34.                  return
  35.              end if                 
  36.              yy=(hydra-hydras)/(hydrar-hydras)
  37.              xmat1=(xmat00-xmat0)*(yy**non_lin)+xmat0
  38.          else
  39. c            valeur constante tant que hydra1<hydras         
  40.              xmat1=xmat0
  41.          end if
  42.       else
  43. c        on adopte directement la velaur de reference      
  44.          xmat1=xmat00
  45.       end if
  46.  
  47. c      print*,'ds hydrxmat'
  48. c      print*,xmat0,xmat1,hydra1,hydras,non_lin,erreur
  49.  
  50.       return
  51.       end
  52.  
  53.  
  54.  

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