Télécharger la plaquette Cast3M


fr en

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