Télécharger indica.eso

Retour à la liste

Numérotation des lignes :

indica
  1. C INDICA SOURCE CHAT 05/01/13 00:35:55 5004
  2. SUBROUTINE INDICA(DKT,DKC,IFIS,IPLA,IC,BETJEF)
  3. C
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8 (A-H,O-Z)
  6. C
  7. *
  8. SEGMENT BETJEF
  9. REAL*8 AA,BETA,RB,ALPHA,YOUN,XNU,GFC,GFT,HC,ETA,TDEF,
  10. & TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  11. INTEGER ICC1,ICC2,IMOD,IVIS,ITER,
  12. & ISIM,IBB,IGAU,IZON
  13. ENDSEGMENT
  14. *
  15. * COMMON /DBETJEF/AA,BETAT,RB,ALPHA,YOUN,XNU,GFC,GFT,HC,ETA,TDEF,
  16. * & TCON,DPSTF1,DPSTF2,TETA,PDT,ICC1,ICC2,IMOD,IVIS,ITER,
  17. * & ISIM,IBB,IGAU,IZON
  18. C
  19. C *************************************************
  20. C ** MISE A JOUR DES INDICATEURS D'ENDOMMAGEMENT **
  21. C *************************************************
  22. C-------------------------------------------------------------------
  23. C IC=1 : COURBE de TRACTION
  24. C IC=2 : COURBE de COMPRESSION
  25. C IC=3 : COURBES TRACTION COMPRESSION
  26. IRET=0
  27. C
  28. C ******************* direction de traction *******************
  29. C
  30. IF (IC.EQ.1.OR.IC.EQ.3) THEN
  31. Ft=ALPHA*RB
  32. C
  33. C--------------- Courbe bilineaire --------------------------
  34. C
  35. IF (ICC1.EQ.1) THEN
  36. DKUT=(2.D0*GFT)/(HC*Ft)
  37. DKUT1=0.95*DKUT
  38. IF (DKT.GE.DKUT1) THEN
  39. IFIS=2
  40. ENDIF
  41. ENDIF
  42. C
  43. C--------------- Courbe exponentielle -----------------------
  44. C
  45. IF (ICC1.EQ.2) THEN
  46. DKUT=(3.D0*GFT)/(HC*Ft*(1.D0-EXP(-3.D0)))
  47. IF (DKT.GE.DKUT) THEN
  48. IFIS=2
  49. ENDIF
  50. ENDIF
  51. C
  52. C------------------------------------------------------------
  53. C
  54. C *************** direction de compression ******************
  55. C
  56. ENDIF
  57. IF (IC.EQ.2.OR.IC.EQ.3) THEN
  58. C
  59. Rb=RB
  60. C
  61. C------------------- Courbe Feenstra -----------------------
  62. C
  63. IF (ICC2.EQ.2) THEN
  64. DKE=(4.D0*Rb)/(3.D0*YOUN)
  65. DKU=1.5*(GFC/(HC*Rb))-(11.D0/48.D0)*DKE
  66. IF (DKC.GE.DKE) THEN
  67. IPLA=2
  68. ENDIF
  69. IF (DKC.GE.DKU) THEN
  70. IPLA=3
  71. ENDIF
  72. ENDIF
  73. C
  74. C--------------- Courbe post pic lineaire -----------------
  75. C
  76. IF (ICC2.EQ.1) THEN
  77. DKE=(4.D0*Rb)/(3.D0*YOUN)
  78. DKU=((2.D0*GFC)/(HC*Rb))+DKE
  79. IF (DKC.GE.DKE) THEN
  80. IPLA=2
  81. ENDIF
  82. IF (DKC.GE.DKU) THEN
  83. IPLA=3
  84. ENDIF
  85. ENDIF
  86. C--------------- Courbe bilineaire --------------------------
  87. C
  88. IF (ICC2.EQ.0) THEN
  89. DKU=(2.D0*GFC)/(HC*Rb)
  90. DKLIM=Rb/YOUN
  91. IF (DKC.GE.DKU) THEN
  92. IPLA=3
  93. ENDIF
  94. ENDIF
  95. ENDIF
  96. C
  97. 10 CONTINUE
  98. RETURN
  99. END
  100.  
  101.  
  102.  
  103.  
  104.  

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