Télécharger la plaquette Cast3M


fr en

Télécharger dohti2.eso

Retour à la liste

Numérotation des lignes : 

dohti2

  1. C DOHTI2    SOURCE    CHAT      05/01/12    22:54:07     5004
  2. C DOHTI2    SOURCE    AM        93/01/12    21:18:38     781
  3. C DOHPTR    SOURCE    PLA2      90/07/11    21:19:35
  4.       SUBROUTINE DOHTI2(VELA,MATE,IFOU,WORK,LHOOK,DDHOOK,IRET)
  5. C
  6. C=======================================================================
  7. C
  8. C  MATRICE DE HOOK DES POUTRES DE TIMOSHENKO
  9. C
  10. C  ENTREES
  11. C     VELA() = materiau dans un tableau de travail
  12. C     MATE   = Nom du materiau
  13. C     IFOU   = num{ro d'harmonique de  fourier: IFOUR de CCOPTIO
  14. C     WORK   = tabbleau contenant les caracteristiques de poutres
  15. C     LHOOK  = taille de la matrice de hooke
  16. C
  17. C  SORTIES
  18. C     DDHOOK(LHOOK,LHOOK) = matrice de hooke
  19. C     IRET = 1 si option existante 0 SINON
  20. C
  21. C     PP 4/9/92
  22. C=======================================================================
  23.       IMPLICIT INTEGER(I-N)
  24.       IMPLICIT REAL*8(A-H,O-Z)
  25.       CHARACTER*8 MATE
  26. C
  27.       DIMENSION VELA(*),DDHOOK(LHOOK,*),WORK(*)
  28. C
  29.       IRET=1
  30.       CALL ZERO(DDHOOK,LHOOK,LHOOK)
  31. C
  32. C  ON RECUPERE LES CARACTERISTIQUES
  33. C
  34.       XINZ=WORK(2)
  35.       SD=WORK(1)
  36.       SREDY=WORK(3)
  37. C
  38. C  SECTION "REDUITE"
  39. C
  40.       IF ( SREDY. NE.0. ) THEN
  41.            SRY=SREDY
  42.       ELSE
  43.            SRY=SD
  44.       ENDIF
  45. C
  46. C  MATERIAU ISOTROPE
  47. C
  48.       IF(MATE.EQ.'ISOTROPE') THEN
  49. C
  50. C  LA MATRICE EST DIAGONALE
  51. C
  52.       IF(IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.-3) THEN
  53.       YOU=VELA(1)
  54.       XNU=VELA(2)
  55.       CISA=YOU*0.5D0/(1.D0+XNU)
  56. C
  57.       DDHOOK(1,1)=YOU*SD
  58.       DDHOOK(2,2)=CISA*SRY
  59.       DDHOOK(3,3)=YOU*XINZ
  60. C
  61.           ELSE
  62.           IRET=0
  63.           ENDIF
  64.       ELSE
  65.       IRET=0
  66.       ENDIF
  67.       RETURN
  68.       END
  69.  
  70.  
  71.  

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