Esope source files

   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 :