Télécharger devolj.eso

Retour à la liste

Numérotation des lignes :

  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