Télécharger la plaquette Cast3M


fr en

Télécharger bbgrdk.eso

Retour à la liste

Numérotation des lignes : 

bbgrdk

  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