Télécharger la plaquette Cast3M


fr en

Télécharger betfis.eso

Retour à la liste

Numérotation des lignes : 

betfis

  1. C BETFIS    SOURCE    CB215821  16/04/21    21:15:20     8920
  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(ABS(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=ABS(S1X)
  30.           IF(S1X .GE.0.D0) EQSTR=0.D0
  31.          ELSE
  32.           EPEQ=ABS(STRNRX)
  33.           CALL CDCY(EPEQ,SEQ,IPLA,TANG,EX,RB,EMAX,EPO1,ICAL)
  34.           EQSTR=SEQ
  35.           S1X=-SEQ
  36.          ENDIF
  37.         RETURN
  38.         ENDIF
  39. C==================================================================
  40.          IF(EPSRX.GT.EPST) THEN
  41. C
  42. C          PENTE DE DECHARGE ---> EDT
  43. C
  44.           S1X=SIGRX+EDT*STRNX
  45.           TANG=EDT
  46. C
  47. C       REOUVERTURE DE LA FISSURE
  48. C
  49.            IF(S1X.GT.RTM) THEN
  50.             CALL FISSU2(S1X,SIGRX,STRNX,RTM,PENT,TANG)
  51.             IFISU=1
  52.            ENDIF
  53.          RETURN
  54. C
  55.          ELSE
  56. C
  57. C          PENTE DE DECHARGE ---> EDC
  58. C
  59.          PASDC= 'BETFIS1'
  60.          CALL DECHAC(S1X,EDC,EPST,EDT,RTM,EPSC,DEFR,SIGRX,STRNX,
  61.      1   STRNRX,EPSRX,IFISU,JFISU,SIGMRX,IPLA,TANG,PENT,EX,RB,ALPHA,
  62.      2   EPO1,IBB1,IGAU1,IDIR,PASDC)
  63.          RETURN
  64.          ENDIF
  65. C        ENDIF
  66.       ENDIF
  67. C======================================================================
  68. C
  69. C       POINT AYANT SUBIT UNE DECHARGE EN COMPRESSION (EDC .NE. 0)
  70. C
  71.        IF(ABS(SIGMRX).GT.EQSTR.AND.SIGMRX.LE.0.D0) THEN
  72. C======================================================================
  73. C    ************************
  74. C    * DOMAINES DES CHARGES *
  75. C    ************************
  76. C
  77. C       ECROUISSAGE ?
  78. C
  79.         IF(STRNRX.LT.EPSC) THEN
  80. C
  81. C        ON EST SUR LA COURBE DE COMPRESSION (CHARGE)
  82. C
  83.           EPEQ=ABS(STRNRX)
  84.           CALL CDCY(EPEQ,SEQ,IPLA,TANG,EX,RB,EMAX,EPO1,ICAL)
  85.           EQSTR=SEQ
  86.           S1X=-SEQ
  87.         ELSE
  88. C
  89. C        ON EST SUR LA COURBE DE COMPRESSION ---> EDC
  90. C                                         ou ---> EDT
  91.             IF(STRNRX.LT.EPST) THEN
  92.              S1X=EDC*(STRNRX-DEFR)
  93.              TANG=EDC
  94.             ELSE
  95.              S1X=SIGRX+EDT*STRNX
  96.              TANG=EDT
  97.             ENDIF
  98.         ENDIF
  99. C
  100.       ELSE
  101. C======================================================================
  102. C    **************************
  103. C    * DOMAINES DES DECHARGES *
  104. C    **************************
  105.        IF(EPSRX.LT.EPSC) THEN
  106. C
  107. C       DECHARGE A PARTIR DE LA COURBE DE COMPRESSION
  108. C       LE CALCUL DES NOUVELLES VALEURS EPSC/EPST/EDC sera effectu‚
  109. C                                                     dans DECHARC
  110. C
  111.          PASDC= 'BETFIS2'
  112.          CALL DECHAC(S1X,EDC,EPST,EDT,RTM,EPSC,DEFR,SIGRX,STRNX,
  113.      1   STRNRX,EPSRX,IFISU,JFISU,SIGMRX,IPLA,TANG,PENT,EX,RB,ALPHA,
  114.      2   EPO1,IBB1,IGAU1,IDIR,PASDC)
  115.          RETURN
  116.        ELSE
  117.         IF(STRNRX.LT.EPST) THEN
  118. C
  119. C       ON EST SUR LA COURBE DE DECHARGE EDC
  120. C
  121.          IF(EPSRX.GT.EPST) THEN
  122. C
  123. C         PASSAGE DE LA COURBE EDT A LA COURBE EDC
  124. C
  125.          S1X=EDC*(STRNRX-DEFR)
  126.          TANG=EDC
  127.          ELSE
  128. C
  129. C         ON RESTE SUR LA COURBE EDC
  130. C
  131.          S1X=SIGRX+EDC*STRNX
  132.          TANG=EDC
  133.          ENDIF
  134.         ELSE
  135. C
  136. C         ON EST SUR LA COURBE DE DECHARGE EDT
  137. C
  138.          TANG=EDT
  139.          IF(EPSRX.GT.EPST) THEN
  140. C
  141. C         ON RESTE SUR LA COURBE EDT
  142. C
  143.           S1X=SIGRX+EDT*STRNX
  144.          ELSE
  145. C
  146. C         PASSAGE DE LA COURBE EDC A LA COURBE EDT
  147. C
  148.           S1X=-RBT+EDT*(STRNRX-EPST)
  149.          ENDIF
  150.         ENDIF
  151.        ENDIF
  152.       ENDIF
  153. C==================================================================
  154.       IF(S1X.GT.RTM) THEN
  155. C
  156. C       REOUVERTURE DE LA FISSURE
  157. C
  158.        CALL FISSU2(S1X,SIGRX,STRNX,RTM,PENT,TANG)
  159.        IFISU=1
  160.       RETURN
  161.       ENDIF
  162. C==================================================================
  163.       RETURN
  164.       END
  165.  
  166.  
  167.  
  168.  
  169.  

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