Télécharger devolp.eso

Retour à la liste

Numérotation des lignes :

  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