Télécharger la plaquette Cast3M


fr en

Télécharger doco88.eso

Retour à la liste

Numérotation des lignes : 

doco88

  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