Télécharger endame.eso

Retour à la liste

Numérotation des lignes :

  1. C ENDAME SOURCE CHAT 05/01/12 23:38:48 5004
  2. SUBROUTINE ENDAME(IC,BETINSA)
  3. C
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. SEGMENT BETINSA
  7. REAL*8 RT,RC,YOUN,XNU,GFT,GFC,HC
  8. REAL*8 DKT,DKC,SEQT,SEQC,ENDT,ENDC
  9. INTEGER IFIS,IPLA,IBB,IGAU
  10. ENDSEGMENT
  11. C
  12. *--------------------------------------------------------------------*
  13. * *
  14. * CALCUL DE L'ENDOMMAGEMENT *
  15. * _________________________ *
  16. * *
  17. * ENTREES: *
  18. * *
  19. * DKT PARAMETRE D'ECROUISSAGE EN TRACTION *
  20. * DKC PARAMETRE D'ECROUISSAGE EN COMPRESSION *
  21. * IC =1 : COURBE de TRACTION *
  22. * =2 : COURBE de COMPRESSION *
  23. * *
  24. * SORTIES: *
  25. * *
  26. * ENDC ENDOMMAGEMENT DE COMPRESSION *
  27. * ENDT ENDOMMAGEMENT DE TRACTION *
  28. * *
  29. *--------------------------------------------------------------------*
  30. C
  31. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  32. C Initialisation
  33. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  34. C
  35. BB=3.D0
  36. DT = 0.5D0
  37. DC = 0.25D0
  38. ENDL = 0.99D0
  39. C
  40. C
  41. C
  42. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  43. C Direction de traction
  44. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  45. C
  46. IF (IC.EQ.1) THEN
  47. FT=RT
  48. AT= -0.5D0
  49. BT=FT*(HC/GFT)*(1.D0+AT/2.D0)
  50. CT=BT*LOG(1.D0-DT)/(LOG(((1+AT)-SQRT(1+AT*AT))/(2.D0*AT)))
  51. C
  52. ENDT=1- EXP(-CT*DKT)
  53. C
  54. IF (ENDT.GE.ENDL) THEN
  55. IFIS=2
  56. ENDIF
  57. IF (ENDT.GT.0.D0) THEN
  58. IFIS=1
  59. ENDIF
  60. ENDIF
  61. C
  62. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  63. C Direction de compression
  64. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  65. C
  66. IF (IC.EQ.2) THEN
  67. FC=RC
  68. AC=(2.D0*BB-1.D0)+2.D0* SQRT(BB**2-BB)
  69. BC=FC*(HC/GFC)*(1.D0+AC/2.D0)
  70. CC=BC*LOG(1.D0-DC)/(LOG((1+AC)/(2.D0*AC)))
  71. C
  72. ENDC=1- EXP(-CC*DKC)
  73. C
  74. IF (ENDC.GE.ENDL) THEN
  75. IPLA=3
  76. ENDIF
  77. IF (ENDC.GE.DC) THEN
  78. IPLA=2
  79. ENDIF
  80. IF (ENDC.GT.0.D0) THEN
  81. IPLA=1
  82. ENDIF
  83. ENDIF
  84. C
  85. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  86. C CAS OU ENDO=ENDO_LIMITE
  87. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  88. C
  89. IF(ENDT.GE.ENDL) ENDT= ENDL
  90. IF(ENDC.GE.ENDL) ENDC= ENDL
  91. C
  92. IF ((ENDT.LT.0.D0).OR.(ENDC.LT.0.D0)) THEN
  93. WRITE(*,*)'PROBLEME DANS DAMAGE'
  94. WRITE(*,*)'ENDT=',ENDT,'ENDC=',ENDC
  95. C STOP
  96. ENDIF
  97. C
  98. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  99. C
  100. RETURN
  101. END
  102.  
  103.  
  104.  

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