Télécharger mfdkt.eso

Retour à la liste

Numérotation des lignes :

mfdkt
  1. C MFDKT SOURCE CHAT 05/01/13 01:44:32 5004
  2. SUBROUTINE MFDKT(QSI,ETA,XEL,VNM)
  3. C=======================================================================
  4. C
  5. C EVALUATION DE LA MATRICE DE MASSE DKT D.D.L. FLEXION
  6. C
  7. C ENTREES
  8. C QSI ETA COORDONNEES REDUITES DU POINT DE GAUSS
  9. C XEL(3,3)=COORDONNEES LOCALES DE L ELEMENT
  10. C SORTIES
  11. C VNM(9) =FONCTION DE FORME POUR LA MASSE
  12. C EBERSOLT AOUT 85
  13. C=======================================================================
  14. IMPLICIT INTEGER(I-N)
  15. IMPLICIT REAL*8(A-H,O-Z)
  16. DIMENSION XEL(3,*),VNM(*)
  17. DATA UN,DEUX,TROIS/1.D0,2.D0,3.D0/
  18. DATA XZER/0.D0/
  19. C
  20. C MISE A ZERO DE BB
  21. C
  22. DO 10 IA=1,9
  23. VNM(IA)=XZER
  24. 10 CONTINUE
  25. C
  26. C CALCUL DES FONCTIONS D INTERPOLATION
  27. C
  28. AL=UN-QSI-ETA
  29. A=QSI*ETA*AL
  30. C
  31. VNM(1)=(AL**2)*(TROIS-DEUX*AL)+DEUX*A
  32. VNM(2)=(AL**2)*QSI+(A/DEUX)
  33. VNM(3)=(AL**2)*ETA+(A/DEUX)
  34. VNM(4)=(TROIS-DEUX*QSI)*(QSI**2)+DEUX*A
  35. VNM(5)=(QSI-UN)*(QSI**2)-A
  36. VNM(6)=(QSI**2)*ETA+(A/DEUX)
  37. VNM(7)=(TROIS-DEUX*ETA)*(ETA**2)+DEUX*A
  38. VNM(8)=(ETA**2)*QSI+(A/DEUX)
  39. VNM(9)=(ETA-UN)*(ETA**2)-A
  40. C
  41. C VARIABLES DANS LE REPERE X,Y
  42. C
  43. X21=XEL(1,2)-XEL(1,1)
  44. X13=XEL(1,1)-XEL(1,3)
  45. Y13=XEL(2,1)-XEL(2,3)
  46. C
  47. C =VNM(2)*X21-VNM(3)*X13
  48. VNM(2)=-VNM(3)*Y13
  49. VNM(3)=-C
  50. C =VNM(5)*X21-VNM(6)*X13
  51. VNM(5)=-VNM(6)*Y13
  52. VNM(6)=-C
  53. C =VNM(8)*X21-VNM(9)*X13
  54. VNM(8)=-VNM(9)*Y13
  55. VNM(9)=-C
  56. C
  57. RETURN
  58. END
  59.  
  60.  

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