Télécharger la plaquette Cast3M


fr en

Télécharger cubksi.eso

Retour à la liste

Numérotation des lignes : 

cubksi

  1. C CUBKSI    SOURCE    CHAT      07/10/22    21:15:34     5921
  2.       SUBROUTINE CUBKSI(NPINT,XXG,BKSI)
  3.       implicit real*8(A-H,O-Z)
  4.       implicit integer (I-N)
  5.  
  6. *      SUBROUTINE CUB_KSI(NPINT,XXG,BKSI)
  7. *
  8. *     ------------------------------------------------------------------
  9. *
  10. *         DERIVEE P/R VARIABLES CANONIQUE AUX POINTS D INTEGRATION
  11. *          POUR LES ELEMENTS CUBES A 8 NOEUDS
  12. *                                                    H. BUNG     03-93
  13. *     ------------------------------------------------------------------
  14. *
  15. *
  16. *   ENTREES
  17. *     NPINT   : NBRE DE  PTS D INTEGRATION
  18. *     XG      : COOR. CANONIQUES DES  PTS D INTEGRATION
  19. *
  20. *   SORTIES :
  21. *     BKSI(3,8,NPINT) : LES DERIVEES
  22. *
  23. *      IMPLICIT NONE
  24. *
  25. *---    VARIABLES GLOBALES
  26. *
  27. *      INTEGER NPINT
  28. *      REAL *8 XG(3,5),XXG(5),GS(8,4)
  29. *      REAL *8 BKSI(3,8,5)
  30.        dimension XG(3,5),XXG(*),GS(8,4),BKSI(3,8,*)
  31.  
  32. *
  33. *---    VARIABLES LOCALES
  34. *
  35. *      REAL *8 AUX
  36. *      REAL *8 UN,UNS8
  37. *      INTEGER IP,I
  38.        dimension XP(3),XM(3)
  39.  
  40.       UN=1.D0
  41.       UNS8=0.125D0
  42. *
  43. C MODIF ANTOINE:
  44. C ON REMPLIT LE TABLEAU XG
  45. C
  46.       DO IP=1,NPINT
  47.       XG(1,IP)=0.d0
  48.       XG(2,IP)=0.d0
  49.       XG(3,IP)=XXG(IP)
  50.       END DO
  51. *
  52. *
  53.       DO IP=1,NPINT
  54. *-
  55.         DO  I=1,3
  56.           XP(I) = UN + XG(I,IP)
  57.           XM(I) = UN - XG(I,IP)
  58.         END DO
  59. *-
  60.         AUX= -UNS8*XM(3)
  61.         BKSI(1,1,IP)= AUX*XM(2)
  62.         BKSI(1,2,IP)=-BKSI(1,1,IP)
  63.         BKSI(1,3,IP)=-AUX*XP(2)
  64.         BKSI(1,4,IP)=-BKSI(1,3,IP)
  65.         BKSI(2,1,IP)= AUX*XM(1)
  66.         BKSI(2,2,IP)= AUX*XP(1)
  67.         BKSI(2,3,IP)=-BKSI(2,2,IP)
  68.         BKSI(2,4,IP)=-BKSI(2,1,IP)
  69.         AUX= -UNS8*XP(3)
  70.         BKSI(1,5,IP)= AUX*XM(2)
  71.         BKSI(1,6,IP)=-BKSI(1,5,IP)
  72.         BKSI(1,7,IP)=-AUX*XP(2)
  73.         BKSI(1,8,IP)=-BKSI(1,7,IP)
  74.         BKSI(2,5,IP)= AUX*XM(1)
  75.         BKSI(2,6,IP)= AUX*XP(1)
  76.         BKSI(2,7,IP)=-BKSI(2,6,IP)
  77.         BKSI(2,8,IP)=-BKSI(2,5,IP)
  78.         AUX= -UNS8*XM(2)
  79.         BKSI(3,1,IP)= AUX*XM(1)
  80.         BKSI(3,2,IP)= AUX*XP(1)
  81.         BKSI(3,5,IP)=-BKSI(3,1,IP)
  82.         BKSI(3,6,IP)=-BKSI(3,2,IP)
  83.         AUX= -UNS8*XP(2)
  84.         BKSI(3,3,IP)= AUX*XP(1)
  85.         BKSI(3,4,IP)= AUX*XM(1)
  86.         BKSI(3,7,IP)=-BKSI(3,3,IP)
  87.         BKSI(3,8,IP)=-BKSI(3,4,IP)
  88. *
  89.       END DO
  90. *
  91.       END
  92.  
  93.  
  94.  
  95.  

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