Télécharger dohmur.eso

Retour à la liste

Numérotation des lignes :

  1. C DOHMUR SOURCE CHAT 05/01/12 22:53:48 5004
  2. SUBROUTINE DOHMUR(VELA,MATE,IFOU,WORK,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 WORK = tabbleau contenant les caracteristiques du mur global
  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 PP 27/04/04
  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. SECT=WORK(1)
  33. XINZ=WORK(2)
  34. C
  35. C MATERIAU ISOTROPE
  36. C
  37. IF(MATE.EQ.'ISOTROPE') THEN
  38. C
  39. C LA MATRICE EST DIAGONALE
  40. C
  41. IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN
  42. YOU=VELA(1)
  43. XNU=VELA(2)
  44. CISA=YOU*0.5D0/(1.D0+XNU)
  45. C
  46. DDHOOK(1,1)=YOU *SECT
  47. DDHOOK(2,2)=YOU *XINZ
  48. DDHOOK(3,3)=CISA*SECT
  49. DDHOOK(4,4)=YOU *SECT
  50. DDHOOK(5,5)=YOU *XINZ
  51. C
  52. ELSE
  53. IRET=0
  54. ENDIF
  55. ELSE
  56. IRET=0
  57. ENDIF
  58. C
  59. RETURN
  60. END
  61.  
  62.  
  63.  

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