Télécharger la plaquette Cast3M


fr en

Télécharger dohtif.eso

Retour à la liste

Numérotation des lignes : 

dohtif

  1. C DOHTIF    SOURCE    CHAT      05/01/12    22:54:10     5004
  2.       SUBROUTINE DOHTIF(VELA,MATE,IFOU,LHOOK,DDHOOK,IRET)
  3. C
  4. C=======================================================================
  5. C
  6. C  MATRICE DE HOOK DES POUTRES DE TIMOSHENKO
  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=======================================================================
  19. C     Pierre Pegon (ISPRA) Juillet/Aout 1993
  20. C=======================================================================
  21.       IMPLICIT INTEGER(I-N)
  22.       IMPLICIT REAL*8(A-H,O-Z)
  23.       CHARACTER*8 MATE
  24. C
  25.       DIMENSION VELA(*),DDHOOK(LHOOK,*)
  26. C
  27.       IRET=1
  28.       CALL ZERO(DDHOOK,LHOOK,LHOOK)
  29. C
  30. C  MATERIAU ISOTROPE
  31. C
  32.       IF(MATE.EQ.'SECTION') THEN
  33. C
  34.         IF(IFOU.EQ.2) THEN
  35. C
  36. C  LA MATRICE EST PLEINE
  37. C
  38.           DDHOOK(1,1)=VELA(1)
  39.           DDHOOK(1,5)=VELA(3)
  40.            DDHOOK(5,1)=VELA(3)
  41.           DDHOOK(1,6)=-VELA(2)
  42.             DDHOOK(6,1)=-VELA(2)
  43. C
  44.           DDHOOK(5,5)=VELA(6)
  45.           DDHOOK(5,6)=-VELA(5)
  46.             DDHOOK(6,5)=-VELA(5)
  47. C
  48.           DDHOOK(6,6)=VELA(4)
  49. C
  50.           DDHOOK(2,2)=VELA(8)
  51.           DDHOOK(2,4)=-VELA(10)
  52.             DDHOOK(4,2)=-VELA(10)
  53. C
  54.           DDHOOK(3,3)=VELA(7)
  55.           DDHOOK(3,4)=VELA(9)
  56.           DDHOOK(4,3)=VELA(9)
  57. C
  58.           DDHOOK(4,4)=VELA(11)+VELA(12)
  59. C
  60.         ELSEIF(IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.-3) THEN
  61. C
  62. C  LA MATRICE EST PLEINE
  63. C
  64.           DDHOOK(1,1)=VELA(1)
  65.           DDHOOK(1,3)=-VELA(2)
  66.           DDHOOK(3,1)=-VELA(2)
  67. C
  68.           DDHOOK(3,3)=VELA(3)
  69. C
  70.           DDHOOK(2,2)=VELA(4)
  71. C
  72.         ELSE
  73.           IRET=0
  74.         ENDIF
  75.       ELSE
  76.         IRET=0
  77.       ENDIF
  78. C
  79.       RETURN
  80.       END
  81.  
  82.  
  83.  

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