Sources Esope | Cast3M

jepg3d
C JEPG3D    SOURCE    FD218221  24/02/07    21:15:17     11834          
      subroutine jepg3d(NBRINC,NINC,NDIMG,NDIMA,NDIMPL,FRAC,GCH,Mbg,bg,
     # JEA,cg,AFFICHE,ERR1)
 
c     calcul du jacobien  d eps e / d epsp generalise (NDIMG*NDIMPL) 
c     stocke dans cg() (cas avec une seule inclusion)    
 
      implicit real*8 (a-h,o-z)
      implicit integer (i-n)
 
c     nbre d inc max et numero de la contrainte du critere actif de traction       
      integer NBRINC,ERR1,NDIMG,NINC,NDIMPL
c     Jacobienne dVe/dVA
      real*8 JEA(NDIMG,NDIMA) 
c     fraction des phases
      real*8 FRAC(0:NBRINC)
c     indicateur de diffusion produits chimiques dans fissure
      real*8 GCH(0:NBRINC)
c     coeff de biot et module de biot des gels
      real*8 bg(0:NBRINC,0:1),Mbg(0:NBRINC,0:1)   
c     jacobienne Je1epg=cg (cf maple)
      real*8 cg(NDIMG,NDIMPL)  
c     options d affichage     
      logical affiche,affiche_local
c     fraction inclusion 
      real*8 f           
      if(NBRINC.ne.1) then
        print*,'Nbr inc inattendu ds jbgp3d'
        err1=1
        return
      end if
      affiche_local=affiche
      err1=0 
 
      f=frac(1)
c     pour la suite cf fichier maple