Télécharger la plaquette Cast3M


fr en

Télécharger dohtim.eso

Retour à la liste

Numérotation des lignes : 

dohtim

  1. C DOHTIM    SOURCE    CHAT      05/01/12    22:54:14     5004
  2. C DOHPTR    SOURCE    PLA2      90/07/11    21:19:35
  3.       SUBROUTINE DOHTIM(VELA,MATE,IFOU,WORK,LHOOK,DDHOOK,IRET)
  4. C
  5. C=======================================================================
  6. C
  7. C  MATRICE DE HOOK DES POUTRES DE TIMOSHENKO
  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     WORK   = tabbleau contenant les caracteristiques de poutres
  14. C     LHOOK  = taille de la matrice de hooke
  15. C
  16. C  SORTIES
  17. C     DDHOOK(LHOOK,LHOOK) = matrice de hooke
  18. C     IRET = 1 si option existante 0 SINON
  19. C
  20. C     PP 4/9/92
  21. C=======================================================================
  22.       IMPLICIT INTEGER(I-N)
  23.       IMPLICIT REAL*8(A-H,O-Z)
  24.       CHARACTER*8 MATE
  25. C
  26.       DIMENSION VELA(*),DDHOOK(LHOOK,*),WORK(*)
  27. C
  28.       IRET=1
  29.       CALL ZERO(DDHOOK,LHOOK,LHOOK)
  30. C
  31. C  ON RECUPERE LES CARACTERISTIQUES
  32. C
  33.       XINX=WORK(1)
  34.       XINY=WORK(2)
  35.       XINZ=WORK(3)
  36.       SD=WORK(4)
  37.       SREDY=WORK(5)
  38.       SREDZ=WORK(6)
  39. C
  40. C  SECTION "REDUITE"
  41. C
  42.       IF ( SREDZ. NE.0. ) THEN
  43.            SRZ=SREDZ
  44.       ELSE
  45.            SRZ=SD
  46.       ENDIF
  47.       IF ( SREDY. NE.0. ) THEN
  48.            SRY=SREDY
  49.       ELSE
  50.            SRY=SD
  51.       ENDIF
  52. C
  53. C  MATERIAU ISOTROPE
  54. C
  55.       IF(MATE.EQ.'ISOTROPE') THEN
  56. C
  57. C  LA MATRICE EST DIAGONALE
  58. C
  59.       IF(IFOU.EQ.2) THEN
  60.       YOU=VELA(1)
  61.       XNU=VELA(2)
  62.       CISA=YOU*0.5D0/(1.D0+XNU)
  63. C
  64.       DDHOOK(1,1)=YOU*SD
  65.       DDHOOK(2,2)=CISA*SRY
  66.       DDHOOK(3,3)=CISA*SRZ
  67.       DDHOOK(4,4)=CISA*XINX
  68.       DDHOOK(5,5)=YOU*XINY
  69.       DDHOOK(6,6)=YOU*XINZ
  70. C
  71.           ELSE
  72.           IRET=0
  73.           ENDIF
  74.       ELSE
  75.       IRET=0
  76.       ENDIF
  77.       RETURN
  78.       END
  79.  
  80.  

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