Télécharger dogo88.eso

Retour à la liste

Numérotation des lignes :

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

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