Télécharger la plaquette Cast3M


fr en

Télécharger devolj.eso

Retour à la liste

Numérotation des lignes : 

devolj

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

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