Télécharger dohtim.eso

Retour à la liste

Numérotation des lignes :

  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