Télécharger douo88.eso

Retour à la liste

Numérotation des lignes :

  1. C DOUO88 SOURCE CHAT 05/01/12 22:54:30 5004
  2. SUBROUTINE DOUO88(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET)
  3. C
  4. C=======================================================================
  5. C
  6. C MATRICE DE HOOKE DES ELEMENTS JOINTS JOI4 ET JOT3 (3D)
  7. C
  8. C ENTREES
  9. C VELA() = materiau dans un tableau de travail
  10. C MATE = Nom du materiau
  11. C IFOU = num{ro d'harmonique de fourier: IFOUR de CCOPTIO
  12. C LHOOK = taille de la matrice de hooke
  13. C
  14. C SORTIES
  15. C DDHOOK(LHOOK,LHOOK) = matrice de hooke
  16. C IRET = 1 si option existante 0 SINON
  17. C
  18. C REMARQUE : ATTENTION : DANS LES CAS CONTRAINTES PLANES,
  19. C DEFORMATIONS PLANES ET AXISYMETRIQUE,
  20. C LA MATRICE DE HOOK N'A PLUS LA DIMENSION
  21. C (3X3), MAIS (2X2). C'EST DONC LES DIMENSIONS DU CAS
  22. C BIDIMENSIONNEL.
  23. C
  24. C CECI EST UN DOUBLE DE DOHO88, MAIS AVEC LA THEORIE DE PARVIS
  25. C POUR LE CALCUL DE LA MATRICE DE TRANSFORMATION. LA MATRICE DE
  26. C ROTATION ROTHOO TRANSFORME UN VECTEUR (EN L'OCCURENCE LE VECTEUR
  27. C DEPLACEMENT) DU REPERE D'ORTHOTROPIE AU REPERE GLOBAL. IL FAUT
  28. C BIEN VOIR QUE CECI EST UNE TRANSFORMATION D'UN VECTEUR EN UN AUTRE
  29. C VECTEUR, ET NON PAS LA TRANSFORMATION D'UN TENSEUR (EX: TENSEUR
  30. C DES CONTRAINTES) EN UN AUTRE TENSEUR. C'EST DONC CETTE SUBROUTINE
  31. C QUI EST LA BONNE, CAR JE RAISONNE EN TERMES DE DEPLACEMENTS, ET
  32. C NON EN TERMES DE CONTRAINTES, DANS LE CAS DES JOINTS (cf. LA
  33. C RELATION FORCES/DEPLACEMENTS ET NON PAS CONTRAINTES/DEFORMATIONS)
  34. C
  35. C=======================================================================
  36. C
  37. IMPLICIT INTEGER(I-N)
  38. IMPLICIT REAL*8(A-H,O-Z)
  39. CHARACTER*8 MATE
  40. C
  41. DIMENSION VELA(*),DDHOOK(LHOOK,*)
  42. REAL*8 D3HOO1(3,3),D3HOO2(3,3),ROTHOO(3,3)
  43. DATA DEUX/2.0D0/
  44. C
  45. IRET=1
  46. CALL ZERO(DDHOOK,LHOOK,LHOOK)
  47. CALL ZERO(D3HOO1,LHOOK,LHOOK)
  48. C
  49. C MATERIAU ISOTROPE
  50. C
  51. IF(MATE.EQ.'ISOTROPE') THEN
  52. C
  53. C CAS TRIDIMENSIONNEL
  54. C
  55. IF(IFOU.EQ.2) THEN
  56. ZK1 =VELA(1)
  57. ZK2 =VELA(1)
  58. ZK3 =VELA(2)
  59. C
  60. DDHOOK(1,1)=ZK1
  61. DDHOOK(2,2)=ZK2
  62. DDHOOK(3,3)=ZK3
  63. C
  64. C CAS CONTRAINTES PLANES DANS LE PLAN (S1,N)
  65. C DEFORMATIONS PLANES DANS LE PLAN (S1,N)
  66. C AXISYMETRIE DANS LE PLAN (S1,N)
  67. C
  68. ELSE IF(IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  69. ZK1 =VELA(1)
  70. ZK2 =VELA(2)
  71. C
  72. DDHOOK(1,1)=ZK1
  73. DDHOOK(2,2)=ZK2
  74. C
  75. C OPTION NON DEFINIE
  76. C
  77. ELSE
  78. IRET=0
  79. ENDIF
  80. C
  81. C MATERIAU ORTHOTROPE
  82. C
  83. ELSE IF(MATE.EQ.'ORTHOTRO') THEN
  84. C
  85. C CAS TRIDIMENSIONNEL
  86. C
  87. IF(IFOU.EQ.2) THEN
  88. ZK1 =VELA(1)
  89. ZK2 =VELA(2)
  90. ZK3 =VELA(3)
  91. COSA=VELA(4)
  92. SINA=VELA(5)
  93. C
  94. D3HOO1(1,1)=ZK1
  95. D3HOO1(2,2)=ZK2
  96. D3HOO1(3,3)=ZK3
  97. C
  98. ROTHOO(1,1) = COSA
  99. ROTHOO(1,2) = -SINA
  100. ROTHOO(1,3) = 0.0D0
  101. ROTHOO(2,1) = SINA
  102. ROTHOO(2,2) = COSA
  103. ROTHOO(2,3) = 0.0D0
  104. ROTHOO(3,1) = 0.0D0
  105. ROTHOO(3,2) = 0.0D0
  106. ROTHOO(3,3) = 1.0D0
  107. C
  108. C PASSAGE DANS LE REPERE DE L'ELEMENT:
  109. C
  110. CALL PRODT (D3HOO2,D3HOO1,ROTHOO,3,3)
  111. DO 400 J=1,3
  112. DO 400 I=1,3
  113. DDHOOK(I,J) = D3HOO2(I,J)
  114. 400 CONTINUE
  115. C
  116. C OPTION NON DEFINIE
  117. C
  118. ELSE
  119. IRET=0
  120. ENDIF
  121. C
  122. C OPTION NON DEFINIE
  123. C
  124. ELSE
  125. IRET=0
  126. ENDIF
  127. C
  128. RETURN
  129. END
  130.  
  131.  
  132.  
  133.  

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