Télécharger la plaquette Cast3M


fr en

Télécharger chgrep1.eso

Retour à la liste

Numérotation des lignes : 

chgrep1

  1. C CHGREP1   SOURCE    CHAT      05/01/12    21:56:44     5004
  2.       SUBROUTINE CHGREP1(COSA,SINA,DHOOK,LHOOK)
  3. *---------------------------------------------------------------------*
  4. *                                                                     *
  5. *  PASSAGE DU REPERE LOCAL DES ELEMENTS AU REPERE D'ORTHOTROPIE       *
  6. *  POUR LA MATRICE DE HOOKE DANS LE CAS DE L'ELEMENT DST              *
  7. *                                                                     *
  8. *---------------------------------------------------------------------*
  9. *                                                                     *
  10. *   ENTREES :                                                         *
  11. *   ________                                                          *
  12. *                                                                     *
  13. *        COSA, SINA = A angle de rotation entre le repère local et le *
  14. *                     repère d'orthotropie                            *
  15. *        DHOOK      = matrice de HOOKE                                *
  16. *        LHOOK      = taille de la matrice de HOOKE                   *
  17. *                                                                     *
  18. *   SORTIES :                                                         *
  19. *   ________                                                          *
  20. *                                                                     *
  21. *        DHOOK      = matrice de HOOKE après le changement de repère  *
  22. *                                                                     *
  23. *---------------------------------------------------------------------*
  24.       IMPLICIT INTEGER(I-N)
  25.       IMPLICIT REAL*8(A-H,O-Z)
  26. *
  27.       DIMENSION DHOOK(LHOOK,LHOOK),RO1HOO(3,3),D3HOO1(3,3)
  28.       DIMENSION D3HOO2(3,3),D6HOO1(3,3)
  29.       DIMENSION D6HOO2(3,3),RO4HOO(3,3),RO3HOO(3,3)
  30.       DIMENSION RO2HOO(2,2),D2HOO1(2,2),D2HOO2(2,2)
  31.       PARAMETER (DEUX=2.D0)
  32. *
  33.       CALL ZERO(D3HOO1,3,3)
  34.       CALL ZERO(D3HOO2,3,3)
  35.       CALL ZERO(D6HOO1,3,3)
  36.       CALL ZERO(D6HOO2,3,3)
  37. *
  38.             DO 1000 I=1,3
  39.              DO 1000 J=1,3
  40.                D3HOO1(I,J) = DHOOK(I,J)
  41.                D6HOO1(I,J) = DHOOK(I+3,J+3)
  42.  1000       CONTINUE
  43. *
  44.  
  45.             COS2 = COSA**2
  46.             SIN2 = SINA**2
  47.             SINCOS = SINA * COSA
  48.             RO1HOO(1,1) = COS2
  49.             RO1HOO(1,2) = SIN2
  50.             RO1HOO(1,3) = - SINCOS
  51.             RO1HOO(2,1) = SIN2
  52.             RO1HOO(2,2) = COS2
  53.             RO1HOO(2,3) =  SINCOS
  54.             RO1HOO(3,1) = DEUX*SINCOS
  55.             RO1HOO(3,2) = - DEUX*SINCOS
  56.             RO1HOO(3,3) = COS2 - SIN2
  57. C
  58. C        PASSAGE DANS LE REPERE DE L'ELEMENT:
  59.             CALL PRODT (D3HOO2,D3HOO1,RO1HOO,3,3)
  60.             CALL PRODT (D6HOO2,D6HOO1,RO1HOO,3,3)
  61.             DO 1150 I=1,3
  62.              DO 1150 J=1,3
  63.                DHOOK(I,J) = D3HOO2(I,J)
  64.                DHOOK(I+3,J+3) = D6HOO2(I,J)
  65.  1150       CONTINUE
  66. *
  67. C
  68.             DO 1100 I=1,2
  69.              DO 1100 J=1,2
  70.                D2HOO1(I,J) = DHOOK(I+6,J+6)
  71.  1100       CONTINUE
  72. *
  73.             RO2HOO(1,1) = COSA
  74.             RO2HOO(1,2) = -SINA
  75.             RO2HOO(2,1) = SINA
  76.             RO2HOO(2,2) = COSA
  77.             CALL PRODT (D2HOO2,D2HOO1,RO2HOO,2,2)
  78.             DO 1200 I=1,2
  79.              DO 1200 J=1,2
  80.                 DHOOK(6+I,6+J) = D2HOO2(I,J)
  81.  1200       CONTINUE
  82. *
  83.       RETURN
  84.       END
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  

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