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    CB215821  26/08/24    21:15:30     12622          
  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 1201 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.  1201       CONTINUE
  44. *
  45.  
  46.             COS2 = COSA**2
  47.             SIN2 = SINA**2
  48.             SINCOS = SINA * COSA
  49.             RO1HOO(1,1) = COS2
  50.             RO1HOO(1,2) = SIN2
  51.             RO1HOO(1,3) = - SINCOS
  52.             RO1HOO(2,1) = SIN2
  53.             RO1HOO(2,2) = COS2
  54.             RO1HOO(2,3) =  SINCOS
  55.             RO1HOO(3,1) = DEUX*SINCOS
  56.             RO1HOO(3,2) = - DEUX*SINCOS
  57.             RO1HOO(3,3) = COS2 - SIN2
  58. C
  59. C        PASSAGE DANS LE REPERE DE L'ELEMENT:
  60.             CALL PRODT (D3HOO2,D3HOO1,RO1HOO,3,3)
  61.             CALL PRODT (D6HOO2,D6HOO1,RO1HOO,3,3)
  62.             DO 1202 I=1,3
  63.              DO 1150 J=1,3
  64.                DHOOK(I,J) = D3HOO2(I,J)
  65.                DHOOK(I+3,J+3) = D6HOO2(I,J)
  66.  1150       CONTINUE
  67.  1202       CONTINUE
  68. *
  69. C
  70.             DO 1203 I=1,2
  71.              DO 1100 J=1,2
  72.                D2HOO1(I,J) = DHOOK(I+6,J+6)
  73.  1100       CONTINUE
  74.  1203       CONTINUE
  75. *
  76.             RO2HOO(1,1) = COSA
  77.             RO2HOO(1,2) = -SINA
  78.             RO2HOO(2,1) = SINA
  79.             RO2HOO(2,2) = COSA
  80.             CALL PRODT (D2HOO2,D2HOO1,RO2HOO,2,2)
  81.             DO 1204 I=1,2
  82.              DO 1200 J=1,2
  83.                 DHOOK(6+I,6+J) = D2HOO2(I,J)
  84.  1200       CONTINUE
  85.  1204       CONTINUE
  86. *
  87.       RETURN
  88.       END
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  

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