Télécharger la plaquette Cast3M


fr en

Télécharger rag3d.eso

Retour à la liste

Numérotation des lignes : 

rag3d

  1. C RAG3D     SOURCE    FD218221  24/02/07    21:15:24     11834          
  2.       subroutine rag3d(taar,nrjg,tref0,aar0,sr1,srsrag,teta1,
  3.      # teta2,dt,aar1,vrag00,Kb,Kgel,bgel,xmgel,dphigel,bg0,bg1,
  4.      # vvgel,cgel,rt,alat,err1)
  5.  
  6. c      sous programme de calcul de l avancement chimique de rag
  7.  
  8.  
  9. c   ********************************************************************
  10.        implicit real*8 (a-h,o-z)
  11.        implicit integer (i-n)
  12. c   ********************************************************************
  13.  
  14.       real*8 taar,nrjg,tref0,aar0,sr1,srsrag,teta1,teta2     
  15.       real*8 alpharag,Ear,temp1,aar1,vvgel,cgel,rt
  16.       real*8 vrag00,Kb,Kgel,bgel,xmgel,dphigel,bg0,bg1
  17.       integer err1
  18.  
  19.       real*8 vgel0,vgel1,vgel
  20.  
  21. c      print*,"taar,nrjg,tref0,aar0,sr1,srsrag,teta1,dt"
  22. c      print*, "vrag00,aar1,vrag1"
  23. c      print*,taar,nrjg,tref0,aar0,sr1,srsrag,teta1,dt,
  24. c     # vrag00,aar1,vrag1 
  25.  
  26. c     **** avancement chimique ***************************************** 
  27.  
  28.       call avch3d(vrag00,taar,nrjg,tref0,sr1,srsrag,teta1,
  29.      # teta2,dt,aar0,aar1,err1)
  30.  
  31. c     **** coefficients poromecanique **********************************
  32.       if(vrag00.gt.0.) then
  33.          if(aar0.gt.alat) then
  34.              vgel1=(aar1-alat)*vrag00/(1.d0-alat)
  35.              vgel0=(aar0-alat)*vrag00/(1.d0-alat)                      
  36.          else
  37.              vgel0=0.d0       
  38.              if(aar1.gt.alat) then
  39.                  vgel1=(aar1-alat)*vrag00
  40. c                /(1.d0-alat)
  41.              else
  42.                  vgel1=0.d0
  43.              end if
  44.          end if
  45.          dphigel=vgel1-vgel0
  46. c        evaluation des coeffs de Biot pour le gel de rag
  47. c        vgel=0.5d0*(vgel0+vgel1)
  48.          vgel=vgel1
  49.          bg1=2.d0*vgel1/(1.d0+vgel1)
  50. c        bgel=0.5d0*(bg0+bg1)
  51.          bgel=bg1
  52.          if(bgel.ge.1.) then
  53.              print*,'bgel trop grand dans rag3d bgel=',bgel
  54.              err1=1
  55.              return
  56.          end if 
  57. c        module de Biot gel-beton        
  58.          call mbio3d(Kb,vgel,bgel,Kgel,xmgel,err1)
  59. c        prise en compte des vides connectee         
  60.          xmgel=xmgel/(1.0d0+xmgel*Vvgel*Cgel/Rt)       
  61.       else
  62. c       pas de rag      
  63.         dphigel=0.d0
  64.         bgel=0.d0
  65.         bg1=0.d0
  66.         xmgel=0.d0
  67.       end if
  68.  
  69.  
  70. c      print*,'rag3d',bgel,xmgel,kb,vgel,kgel,cgel      
  71.  
  72.       return
  73.  
  74.       end
  75. c***********************************************************************      
  76.  
  77.  
  78.  
  79.  

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