Télécharger dohti2.eso

Retour à la liste

Numérotation des lignes :

  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