Télécharger bbgrdk.eso

Retour à la liste

Numérotation des lignes :

  1. C BBGRDK SOURCE CHAT 05/01/12 21:36:32 5004
  2. SUBROUTINE BBGRDK(QSI,ETA,GEOM,BB)
  3. C |====================================================================|
  4. C | |
  5. C | EVALUATION DE LA MATRICE BB(2,9) AU POINT QSI ETA |
  6. C | POUR LES TERMS DE GRADIAN DE L'ELEMENT DKT |
  7. C | DUXDZ DUYDZ DUZDX DUZDY MAIS |
  8. C | DUXDZ=-DUZDX=+BETA(X)=BB(1,I) (I=1,2...9) |
  9. C | DUYDZ=-DUZDY=+BETA(Y)=BB(2,I) (I=1,2...9) |
  10. C | |
  11. C | CODE SUO X.Z. 23/01/1987 |
  12. C |====================================================================|
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8 (A-H,O-Z)
  15. DIMENSION GEOM(*),BB(2,*)
  16. C
  17. DATA TRDE/1.5D0/,DEUX/2.D0/,QUATRE/4.D0/,UN/1.D0/
  18. C
  19. A4=GEOM(1)
  20. A5=GEOM(2)
  21. A6=GEOM(3)
  22. B4=GEOM(4)
  23. B5=GEOM(5)
  24. B6=GEOM(6)
  25. C4= GEOM(7)
  26. C5= GEOM(8)
  27. C6= GEOM(9)
  28. D4=GEOM(10)
  29. D5=GEOM(11)
  30. D6=GEOM(12)
  31. E4=GEOM(13)
  32. E5=GEOM(14)
  33. E6=GEOM(15)
  34. AUX=UN-QSI-ETA
  35. C
  36. C FONCTIONS DE FORME DE TRIANGLE A 6 NOEUDS
  37. C
  38. SN1=-AUX*(1-DEUX*AUX)
  39. SN2=-QSI*(1-DEUX*QSI)
  40. SN3=-ETA*(1-DEUX*ETA)
  41. SN4=QUATRE*QSI*ETA
  42. SN5=QUATRE*ETA*AUX
  43. SN6=QUATRE*QSI*AUX
  44. C
  45. C FONCTIONS DE FORME DE TRIANGLE DKT
  46. C
  47. BB(1,1)= TRDE*(A6*SN6-A5*SN5)
  48. BB(1,2)= B5*SN5+B6*SN6
  49. BB(1,3)= SN1-C5*SN5-C6*SN6
  50. BB(1,4)= TRDE*(A4*SN4-A6*SN6)
  51. BB(1,5)= B6*SN6+B4*SN4
  52. BB(1,6)= SN2-C6*SN6-C4*SN4
  53. BB(1,7)= TRDE*(A5*SN5-A4*SN4)
  54. BB(1,8)= B4*SN4+B5*SN5
  55. BB(1,9)= SN3-C4*SN4-C5*SN5
  56. BB(2,1)= TRDE*(D6*SN6-D5*SN5)
  57. BB(2,2)= -SN1+E5*SN5+E6*SN6
  58. BB(2,3)= -B5*SN5-B6*SN6
  59. BB(2,4)= TRDE*(D4*SN4-D6*SN6)
  60. BB(2,5)= -SN2+E6*SN6+E4*SN4
  61. BB(2,6)= -B6*SN6-B4*SN4
  62. BB(2,7)= TRDE*(D5*SN5-D4*SN4)
  63. BB(2,8)= -SN3+E4*SN4+E5*SN5
  64. BB(2,9)= -B4*SN4-B5*SN5
  65. C
  66. C
  67. RETURN
  68. END
  69.  
  70.  

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