Télécharger chgrep1.eso

Retour à la liste

Numérotation des lignes :

  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