Esope source files

   1 : C DEVOLP    SOURCE    CHAT      05/01/12    22:46:54     5004           
   2 :       SUBROUTINE DEVOLP(XEL,SHP,MFR,NBBB,NBNO,IFOU,NIFOU,
   3 :      +     DIM3,RR,DJAC)
   4 : C=======================================================================
   5 : C     
   6 : C    CALCULE L'ELEMENT DE VOLUME A PARTIR DU JACOBIEN POUR
   7 : C    MILIEUX POREUX
   8 : C     
   9 : C=======================================================================
  10 : C  ENTREES
  11 : C     XEL(3,NBNO)=COORDONNEES LOCALES DE L ELEMENT
  12 : C     SHP(6,NBNO)=DERIVEES PAR RAPPORT A LA GEOMETRIE DE REFERENCE
  13 : C     MFR=FORMULATION DE L'ELEMENT FINI
  14 : C     NBNO=NOMBRE DE FONCTIONS DE FORME
  15 : C     NBBB=NOMBRE DE NOEUDS
  16 : C     IFOU=IFOUR DE CCOPTIO
  17 : C     NIFOU=NIFOUR DE CCOPTIO
  18 : C     DIM3=EPAISSEUR DE L'ELEMENT (CAS MASSIF CONTRAINTES PLANES)
  19 : C     DJAC =JACOBIEN
  20 : C     
  21 : C  SORTIE
  22 : C     SHP(6,NBNO)=DERIVEES PAR RAPPORT A LA GEOMTRIE DE L ELEMENT
  23 : C     RR=RAYON AU POINT DE GAUSS (POUR L AXISYMMETRIQUE)
  24 : C     DJAC =ELEMENT DE VOLUME
  25 : C     
  26 : C=======================================================================
  27 :       IMPLICIT INTEGER(I-N)
  28 :       IMPLICIT REAL*8(A-H,O-Z)
  29 : C    Include contenant quelques constantes dont XPI :
  30 : -INC CCREEL
  31 : *     
  32 :       DIMENSION XEL(3,*),SHP(6,*)
  33 : *     
  34 : *     
  35 :       IF (IFOU.EQ.2) THEN
  35 :          CALL JACOBP(XEL,SHP,3,NBBB,NBNO,DJAC)
  36 :       ELSE
  36 :          CALL JACOBP(XEL,SHP,2,NBBB,NBNO,DJAC)
  37 :       ENDIF
  38 : *     
  39 :       IF(IFOU.EQ.-2) DJAC=DJAC*DIM3
  40 : *     
  41 :       IF(IFOU.EQ.0.OR.IFOU.EQ.1) THEN
  41 :          CALL DISTRR(XEL,SHP,NBBB,RR)
  42 :          IF(IFOU.EQ.0.OR.(IFOU.EQ.1.AND.NIFOU.EQ.0)) THEN
  43 :             DJAC=DJAC*RR*2*XPI
  44 :          ELSE
  45 :             DJAC=DJAC*RR*XPI
  46 :          ENDIF
  47 :       ENDIF
  48 : *     
  49 :       RETURN
  50 :       END
  51 : 
  52 : 
  53 : 
  54 : 
  55 :  
  56 :  
  57 :