Télécharger la plaquette Cast3M


fr en

Télécharger devolp.eso

Retour à la liste

Numérotation des lignes : 

devolp

  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
  36.          CALL JACOBP(XEL,SHP,3,NBBB,NBNO,DJAC)
  37.       ELSE
  38.          CALL JACOBP(XEL,SHP,2,NBBB,NBNO,DJAC)
  39.       ENDIF
  40. *
  41.       IF(IFOU.EQ.-2) DJAC=DJAC*DIM3
  42. *
  43.       IF(IFOU.EQ.0.OR.IFOU.EQ.1) THEN
  44.          CALL DISTRR(XEL,SHP,NBBB,RR)
  45.          IF(IFOU.EQ.0.OR.(IFOU.EQ.1.AND.NIFOU.EQ.0)) THEN
  46.             DJAC=DJAC*RR*2*XPI
  47.          ELSE
  48.             DJAC=DJAC*RR*XPI
  49.          ENDIF
  50.       ENDIF
  51. *
  52.       RETURN
  53.       END
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  

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