Télécharger la plaquette Cast3M


fr en

Télécharger dohptr.eso

Retour à la liste

Numérotation des lignes : 

dohptr

  1. C DOHPTR    SOURCE    CHAT      05/01/12    22:54:04     5004
  2.       SUBROUTINE DOHPTR(VELA,MATE,IFOU,WORK,LHOOK,DDHOOK,IRET)
  3. C
  4. C=======================================================================
  5. C
  6. C  MATRICE DE HOOK DES POUTRES
  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     WORK   = tabbleau contenant les caracteristiques de poutres
  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     Passage en nouveau CHAMELEM par I.Monnier le 4 05 90
  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,*),WORK(*)
  26. C
  27.       IRET=1
  28.       CALL ZERO(DDHOOK,LHOOK,LHOOK)
  29. C
  30. C  ON RECUPERE LES CARACTERISTIQUES
  31. C
  32.       XINX=WORK(1)
  33.       XINY=WORK(2)
  34.       XINZ=WORK(3)
  35.       SD=WORK(4)
  36.       SREDY=WORK(5)
  37.       SREDZ=WORK(6)
  38. C
  39. C  MATERIAU ISOTROPE
  40. C
  41.       IF(MATE.EQ.'ISOTROPE') THEN
  42. C
  43. C  LA MATRICE EST DIAGONALE
  44. C
  45.       IF(IFOU.EQ.2) THEN
  46.       YOU=VELA(1)
  47.       XNU=VELA(2)
  48.       CISA=YOU*0.5D0/(1.D0+XNU)
  49. C
  50.       DDHOOK(1,1)=YOU*SD
  51.       DDHOOK(2,2)=CISA*SREDY
  52.       DDHOOK(3,3)=CISA*SREDZ
  53.       DDHOOK(4,4)=CISA*XINX
  54.       DDHOOK(5,5)=YOU*XINY
  55.       DDHOOK(6,6)=YOU*XINZ
  56. C
  57.           ELSE
  58.           IRET=0
  59.           ENDIF
  60.       ELSE
  61.       IRET=0
  62.       ENDIF
  63.       RETURN
  64.       END
  65.  
  66.  

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