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