Esope source files

   1 : C CAMFLO    SOURCE    AM        13/12/16    21:15:23     7832           
   2 :       SUBROUTINE CAMFLO(IFOUR,SIG,XINV,PI,XMAT,DFDS,DFDPI,H)
   3 : *----------------------------------------------------------------
   4 : *                      MODELE CAM-CLAY
   5 : *            QUANTITES NECESSAIRES A L'ECOULEMENT
   6 : *----------------------------------------------------------------
   7 : *   ENTREE :
   8 : *          IFOUR     CF CCOPTIO
   9 : *          SIG       CONTRAINTES
  10 : *          PI        VARIABLE INTERNE
  11 : *          XMAT      PROPRIETES MATERIELLES
  12 : *
  13 : *   SORTIE :
  14 : *          DFDS      DF / DSIG
  15 : *          DFDPI     DF / DPI
  16 : *          H         LOI D'ECROUISSAGE
  17 : *----------------------------------------------------------------
  18 : *
  19 :       IMPLICIT INTEGER(I-N)
  20 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  21 :       DIMENSION SIG(*),XINV(*),XMAT(*),DFDS(*)
  22 : *
  23 :       E0   = XMAT(5)
  24 :       COHE = XMAT(7)
  25 :       GM   = XMAT(8)
  26 :       XLAMD= XMAT(9)
  27 :       XKAPA= XMAT(10)
  28 : *
  29 :       TETA=(1.D0+E0)/(XLAMD-XKAPA)
  30 :       UNGM2= 1.D0/(GM*GM)
  31 :       AUX=(COHE-PI)*0.5D0 -XINV(1)
  32 :       Y =AUX**2 + (XINV(2)/GM)**2
  33 :       Y = SQRT(Y)
  34 :       UNSY=1.D0/Y
  35 :       FAC =AUX/3.D0
  36 : *
  37 :       DFDS(1)=UNSY*(UNGM2*(SIG(1)-0.5D0*(SIG(2)+SIG(3))) - FAC)
  38 :       DFDS(2)=UNSY*(UNGM2*(SIG(2)-0.5D0*(SIG(1)+SIG(3))) - FAC)
  39 :       DFDS(3)=UNSY*(UNGM2*(SIG(3)-0.5D0*(SIG(1)+SIG(2))) - FAC)
  40 :       DFDS(4)=UNSY*UNGM2*3.D0*SIG(4)
  41 :       IF(IFOUR.GE.1) THEN
  42 :          DFDS(5)=UNSY*UNGM2*3.D0*SIG(5)
  43 :          DFDS(6)=UNSY*UNGM2*3.D0*SIG(6)
  44 :       ENDIF
  45 : *
  46 :       DFDPI=-0.5D0*(1.D0 + AUX*UNSY)
  47 : *
  48 :       H=TETA*(PI+COHE)*AUX*UNSY
  49 : *
  50 :       RETURN                                                                    
  51 :       END                                                                       
  52 :  
  53 :