Télécharger la plaquette Cast3M


fr en

Télécharger redu3d.eso

Retour à la liste

Numérotation des lignes : 

redu3d

  1. C REDU3D    SOURCE    PV090527  23/01/27    21:16:02     11574          
  2.         subroutine redu3d(NBRINC,NINC,NDIMG,DEPSE1,DEPSM1,DEPSK1,
  3.      #  DEPST1,DEPSC1,DEPSL1,LogReducGEL,reduc,dbgvd,Mbg,dbgpg,
  4.      #  precision3d,ierr1,affiche)
  5.  
  6. c       reduction de l increment de la solution dans l espace generalise
  7.  
  8.         implicit real*8 (a-h,o-z)
  9.         implicit integer (i-n) 
  10.  
  11.         integer NBRINC,NINC,NDIMG
  12. c       vecteur deformations solutions       
  13.         real*8 DEPSE1(NDIMG),DEPSM1(NDIMG),DEPSK1(NDIMG)
  14.         real*8 Mbg(0:NBRINC,0:1)
  15. c       poro meca gels
  16.         real*8 dbgPg(0:NBRINC)
  17. c       vecteur solution def anelastiques      
  18.         real*8 DEPST1(NDIMG),DEPSC1(NDIMG)
  19.         real*8 DEPSL1(NDIMG)
  20.         real*8 reduc,precision3d
  21.         integer ierr1
  22.         logical affiche        
  23.         integer idir
  24.         logical affiche_local
  25.         logical LogReducGEL
  26.         real*8 dbgvd(0:NBRINC)       
  27.  
  28.         affiche_local=affiche    
  29.  
  30.         if ((reduc.gt.0.d0).and.
  31.      #      (reduc.lt.1.d0))then          
  32. c           application de la reduction     
  33.             do idir=1,NDIMG
  34. c               elastiques                     
  35.                 DEPSE1(idir)=DEPSE1(idir)*reduc
  36. c               maxwell                
  37.                 DEPSM1(idir)=DEPSM1(idir)*reduc
  38. c               kelvin                
  39.                 DEPSK1(idir)=DEPSK1(idir)*reduc
  40. c               plasticite traction                
  41.                 DEPST1(idir)=DEPST1(idir)*reduc 
  42. c               cisiallement                
  43.                 DEPSC1(idir)=DEPSC1(idir)*reduc 
  44. c               localisee                
  45.                 DEPSL1(idir)=DEPSL1(idir)*reduc                   
  46.             end do
  47.             if(LogReducGEL) then
  48.                 do iphase=0,ninc
  49.                     dbgvd(iphase)=reduc*dbgvd(iphase)
  50.                     dbgpg(iphase)=-MBG(iphase,1)*dbgvd(iphase)
  51.                 end do
  52.             end if
  53.             if(affiche_local) then
  54.                 print*,'Dans redu3d reduc=',reduc
  55.             end if
  56.         else
  57.             if( (reduc.gt.1.d0).or.
  58.      #          (reduc.lt.0.d0)) then
  59.                 print*,'Probleme dans reducinc3d',reduc
  60.                 ierr1=1
  61.                 return
  62.             end if
  63.         end if
  64.  
  65.         return
  66.         end
  67.  
  68.  
  69.  
  70.  

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