Télécharger la plaquette Cast3M


fr en

Télécharger douco2.eso

Retour à la liste

Numérotation des lignes : 

douco2

  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