Télécharger la plaquette Cast3M


fr en

Télécharger douo88.eso

Retour à la liste

Numérotation des lignes : 

douo88

  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