Télécharger doco88.eso

Retour à la liste

Numérotation des lignes :

  1. C DOCO88 SOURCE CHAT 05/01/12 22:53:09 5004
  2. SUBROUTINE DOCO88(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET)
  3. C
  4. C=======================================================================
  5. C
  6. C MATRICE DE HOOKE DES ELEMENTS JOINTS JCT3 ET JCI4 (2D)
  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 P. PEGON 2001 d'apres DOUO88
  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. IRET=1
  47. CALL ZERO(DDHOOK,LHOOK,LHOOK)
  48. CALL ZERO(D3HOO1,LHOOK,LHOOK)
  49. C
  50. C MATERIAU ISOTROPE
  51. C
  52. IF(MATE.EQ.'ISOTROPE') THEN
  53. C
  54. C CAS CONTRAINTES PLANES DANS LE PLAN (S1,N)
  55. C DEFORMATIONS PLANES DANS LE PLAN (S1,N)
  56. C AXISYMETRIE DANS LE PLAN (S1,N)
  57. C
  58. IF(IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  59. C
  60. C et juste le cisaillement!
  61. C
  62. ZK1 =VELA(1)
  63. C
  64. DDHOOK(1,1)=ZK1
  65. DDHOOK(2,2)=ZK1
  66. C
  67. ELSE
  68. IRET=0
  69. ENDIF
  70. ELSE
  71. IRET=0
  72. ENDIF
  73. C
  74. RETURN
  75. END
  76.  
  77.  
  78.  
  79.  
  80.  

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