Esope source files

   1 : C BETFIS    SOURCE    CHAT      05/01/12    21:38:48     5004           
   2 : C BETDJE     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE BETFIS(EPSRX,STRNX,STRNRX,SIGMRX,SIGRX,TANG,S1X,EDT,
   4 :      1 EDC,EPST,EPSC,DEFR,RTM,IFISU,JFISU,IPLA,EQSTR,PENT,EX,RB,ALPHA,
   5 :      2 EMAX,EPO1,ICAL,IBB1,IGAU1,IDIR)
   6 : C
   7 :       IMPLICIT INTEGER(I-N)
   8 :       IMPLICIT REAL*8(A-H,O-Z)
   9 :       CHARACTER*8 PASDC
  10 : C
  11 : C       CE SOUS-PROGRAMME GERE LE COMPORTEMENT DU BETON EN UN POINT
  12 : C       OU LA FISSURE EST REFERMEE (IFISU=0 ET JFISU=1)
  13 : C
  14 :       RBT=ALPHA*RB
  15 :       IF(DABS(EDC).LT.0.1D-05) EDC=0.D0
  16 : C==================================================================
  17 :       IF(EDC.EQ.0.D0.OR.(EQSTR.LE.RBT.AND.IPLA.EQ.0)) THEN
  18 : C-------------------------------------------------------------
  19 : C       POINT INITIALLEMENT TENDU OU PAS COMPRIME PLUS LOIN
  20 : C       QUE -RBT
  21 : C
  22 :         IF(STRNRX.LT.EPSRX) THEN
  23 : C
  24 : C        COMPRESSION DE LA FISSURE (CHARGE)
  25 : C
  26 :          IF(STRNRX.GT.EPST) THEN
  27 :           S1X=SIGRX+EDT*STRNX
  28 :           TANG=EDT
  29 :           EQSTR=DABS(S1X)
  30 :           IF(S1X .GE.0.D0) EQSTR=0.D0
  31 :          ELSE
  32 :           EPEQ=DABS(STRNRX)
  32 :           CALL CDCY(EPEQ,SEQ,IPLA,TANG,EX,RB,EMAX,EPO1,ICAL)
  33 :           EQSTR=SEQ
  34 :           S1X=-SEQ
  35 :          ENDIF
  36 :         RETURN
  37 :         ENDIF
  38 : C==================================================================
  39 :          IF(EPSRX.GT.EPST) THEN
  40 : C
  41 : C          PENTE DE DECHARGE ---> EDT
  42 : C
  43 :           S1X=SIGRX+EDT*STRNX
  44 :           TANG=EDT
  45 : C
  46 : C       REOUVERTURE DE LA FISSURE
  47 : C
  48 :            IF(S1X.GT.RTM) THEN
  48 :             CALL FISSU2(S1X,SIGRX,STRNX,RTM,PENT,TANG)
  49 :             IFISU=1
  50 :            ENDIF
  51 :          RETURN
  52 : C
  53 :          ELSE
  54 : C
  55 : C          PENTE DE DECHARGE ---> EDC
  56 : C
  57 :          PASDC= 'BETFIS1'
  57 :          CALL DECHAC(S1X,EDC,EPST,EDT,RTM,EPSC,DEFR,SIGRX,STRNX,
  58 :      1   STRNRX,EPSRX,IFISU,JFISU,SIGMRX,IPLA,TANG,PENT,EX,RB,ALPHA,
  59 :      2   EPO1,IBB1,IGAU1,IDIR,PASDC)
  60 :          RETURN
  61 :          ENDIF
  62 : C        ENDIF
  63 :       ENDIF
  64 : C======================================================================
  65 : C
  66 : C       POINT AYANT SUBIT UNE DECHARGE EN COMPRESSION (EDC .NE. 0)
  67 : C
  68 :        IF(DABS(SIGMRX).GT.EQSTR.AND.SIGMRX.LE.0.D0) THEN
  69 : C======================================================================
  70 : C    ************************
  71 : C    * DOMAINES DES CHARGES *
  72 : C    ************************
  73 : C
  74 : C       ECROUISSAGE ?
  75 : C
  76 :         IF(STRNRX.LT.EPSC) THEN
  77 : C
  78 : C        ON EST SUR LA COURBE DE COMPRESSION (CHARGE)
  79 : C
  80 :           EPEQ=DABS(STRNRX)
  80 :           CALL CDCY(EPEQ,SEQ,IPLA,TANG,EX,RB,EMAX,EPO1,ICAL)
  81 :           EQSTR=SEQ
  82 :           S1X=-SEQ
  83 :         ELSE
  84 : C
  85 : C        ON EST SUR LA COURBE DE COMPRESSION ---> EDC
  86 : C                                         ou ---> EDT
  87 :             IF(STRNRX.LT.EPST) THEN
  88 :              S1X=EDC*(STRNRX-DEFR)
  89 :              TANG=EDC
  90 :             ELSE
  91 :              S1X=SIGRX+EDT*STRNX
  92 :              TANG=EDT
  93 :             ENDIF
  94 :         ENDIF
  95 : C
  96 :       ELSE
  97 : C======================================================================
  98 : C    **************************
  99 : C    * DOMAINES DES DECHARGES *
 100 : C    **************************
 101 :        IF(EPSRX.LT.EPSC) THEN
 102 : C
 103 : C       DECHARGE A PARTIR DE LA COURBE DE COMPRESSION
 104 : C       LE CALCUL DES NOUVELLES VALEURS EPSC/EPST/EDC sera effectu‚
 105 : C                                                     dans DECHARC
 106 : C
 107 :          PASDC= 'BETFIS2'
 107 :          CALL DECHAC(S1X,EDC,EPST,EDT,RTM,EPSC,DEFR,SIGRX,STRNX,
 108 :      1   STRNRX,EPSRX,IFISU,JFISU,SIGMRX,IPLA,TANG,PENT,EX,RB,ALPHA,
 109 :      2   EPO1,IBB1,IGAU1,IDIR,PASDC)
 110 :          RETURN
 111 :        ELSE
 112 :         IF(STRNRX.LT.EPST) THEN
 113 : C
 114 : C       ON EST SUR LA COURBE DE DECHARGE EDC
 115 : C
 116 :          IF(EPSRX.GT.EPST) THEN
 117 : C
 118 : C         PASSAGE DE LA COURBE EDT A LA COURBE EDC
 119 : C
 120 :          S1X=EDC*(STRNRX-DEFR)
 121 :          TANG=EDC
 122 :          ELSE
 123 : C
 124 : C         ON RESTE SUR LA COURBE EDC
 125 : C
 126 :          S1X=SIGRX+EDC*STRNX
 127 :          TANG=EDC
 128 :          ENDIF
 129 :         ELSE
 130 : C
 131 : C         ON EST SUR LA COURBE DE DECHARGE EDT
 132 : C
 133 :          TANG=EDT
 134 :          IF(EPSRX.GT.EPST) THEN
 135 : C
 136 : C         ON RESTE SUR LA COURBE EDT
 137 : C
 138 :           S1X=SIGRX+EDT*STRNX
 139 :          ELSE
 140 : C
 141 : C         PASSAGE DE LA COURBE EDC A LA COURBE EDT
 142 : C
 143 :           S1X=-RBT+EDT*(STRNRX-EPST)
 144 :          ENDIF
 145 :         ENDIF
 146 :        ENDIF
 147 :       ENDIF
 148 : C==================================================================
 149 :       IF(S1X.GT.RTM) THEN
 150 : C
 151 : C       REOUVERTURE DE LA FISSURE
 152 : C
 152 :        CALL FISSU2(S1X,SIGRX,STRNX,RTM,PENT,TANG)
 153 :        IFISU=1
 154 :       RETURN
 155 :       ENDIF
 156 : C==================================================================
 157 :       RETURN
 158 :       END
 159 :  
 160 :