Télécharger douco2.eso

Retour à la liste

Numérotation des lignes :

  1. C DOUCO2 SOURCE LJ1 14/11/13 21:15:21 8249
  2.  
  3. SUBROUTINE DOUCO2(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET,XL,DIAM)
  4. C
  5. C=======================================================================
  6. C
  7. C MATRICE DE HOOKE DES ELEMENTS COAXIAL COA2 (2D ou 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=======================================================================
  37. C
  38. IMPLICIT INTEGER(I-N)
  39. IMPLICIT REAL*8(A-H,O-Z)
  40. CHARACTER*8 MATE
  41. C
  42. DIMENSION VELA(*),DDHOOK(LHOOK,*)
  43. REAL*8 D3HOO1(3,3),D3HOO2(3,3),ROTHOO(3,3)
  44. DATA DEUX/2.0D0/
  45. C
  46. XPI = 3.1415926535897931D0
  47. C
  48. C
  49. IRET=1
  50. CALL ZERO(DDHOOK,LHOOK,LHOOK)
  51. CALL ZERO(D3HOO1,LHOOK,LHOOK)
  52. C
  53. C MATERIAU ISOTROPE
  54. C
  55. IF(MATE.EQ.'ISOTROPE') THEN
  56. C
  57. C CAS TRIDIMENSIONNEL
  58. C
  59. IF(IFOU.EQ.2) THEN
  60. ZK1 =VELA(1)*xpi*xl*diam
  61. ZK2 =VELA(2)*xl*diam
  62. ZK3 =VELA(2)*xl*diam
  63. C
  64. DDHOOK(1,1)=ZK1
  65. DDHOOK(2,2)=ZK2
  66. DDHOOK(3,3)=ZK3
  67. c write(6,*) 'zk', zk1,zk2,zk3
  68. C
  69. C CAS CONTRAINTES PLANES DANS LE PLAN (S1,N)
  70. C DEFORMATIONS PLANES DANS LE PLAN (S1,N)
  71. C AXISYMETRIE DANS LE PLAN (S1,N)
  72. C
  73. ELSE IF(IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  74. ZK1 =VELA(1)
  75. ZK2 =VELA(2)
  76. C
  77. DDHOOK(1,1)=ZK1
  78. DDHOOK(2,2)=ZK2
  79. C
  80. C OPTION NON DEFINIE
  81. C
  82. ELSE
  83. IRET=0
  84. ENDIF
  85. C
  86. C OPTION NON DEFINIE
  87. C
  88. ELSE
  89. IRET=0
  90. ENDIF
  91. C
  92. RETURN
  93. END
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  

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