Esope source files

   1 : C PLAEND    SOURCE    CHAT      05/01/13    02:14:51     5004           
   2 : C PLAEND     SOURCE    INSL       24/10/96    
   3 :       SUBROUTINE PLAEND(EPSRX,SIGRX,STRNX,S1X,IFISU,IPLA,EQSTR,EPSEQ,
   4 :      1 JFISU,TANG,EPST,EPSC,EDC,EDT,RTM,DEFR,SIGMRX,PENT,EX,RB,ALPHA,
   5 :      2 EMAX,EPSU,EPO1,ICAL,IBB1,IGAU1,IDIR)
   6 : C
   7 : C      MODELE 2 : BETON PLASTIQUE-ENDOMMAGEABLE (OMAR)
   8 : C          ******  VERSION 8 JUILLET 1991  ******
   9 : C          ******  DERNIERE ACTUALISATION  19 JUILLET 1993  ******
  10 : C          ******  DERNIERE ACTUALISATION  20 JUILLET 1995 ******
  11 : C
  12 : C ENTREES : EPSRX : DEFORMATION TOTALE DERNIER POINT CONVERGE
  13 : C           SIGRX : CONTRAINTE DERNIER POINT CONVERGE
  14 : C           STRNX = DEPS(1): INCREMENT DE DEFORMATION
  15 : C
  16 : C  E/S : IFISU,IPLA,EQSTR,EPSEQ,JFISU,TANG,EPST,EPSC,EDC,EDT,RTM,DEFR
  17 : C       CE SONT TOUTES DES VARIABLES INTERNES : VOIR LEUR SIGNIFICATION
  18 : C       DANS LA PROCEDURE APPELANTE : FIBINSA
  19 : C
  20 : C SORTIES :  S1X : LA CONTRAINTE REELLE OU FINALE
  21 : C
  22 : C
  23 :       IMPLICIT INTEGER(I-N)
  24 :       IMPLICIT REAL*8(A-H,O-Z)
  25 :       CHARACTER*8 PASDT,PASDC
  26 : C
  27 : C-----------------------------------------------------------------------
  28 : C           ############################################
  29 : C           *      B E T O N    U N I A X I A L        *
  30 : C           ############################################
  31 :       STRNRX=EPSRX+STRNX
  32 :       EPEQ=DABS(STRNRX)
  33 : C=================================================================
  34 : C      SI LE POINT A DEPASSE LA LIMITE EMAX EN COMPRESSION
  35 : C      ===> LE POINT NE REPREND PLUS DE CONTRAINTES ET LA CONTRAINTE
  36 : C           TOTALE EST NULLE QUELQUE SOIT L ETAT DE DEFORMATION
  37 : C
  38 : C
  39 :       IF(IPLA .EQ. 4) THEN
  40 :        EPSEQ=EPEQ
  41 :        S1X=0.D0
  42 :        TANG=0.D0
  43 :        RETURN
  44 :       ENDIF
  45 : C=================================================================
  46 :       IF(IFISU .EQ. 1) THEN
  47 : C=================================================================
  48 : C              ***************************
  49 : C              *   POINT DEJA FISSURE    *
  50 : C              ***************************
  51 : C            TEST DE FERMETURE DE LA FISSURE
  52 : C
  53 :          PASDT= 'PLAEND21'
  53 :          CALL DECHAT(EPSRX,STRNX,STRNRX,SIGMRX,TANG,S1X,EDT,EDC,
  54 :      1   EPST,DEFR,IFISU,IPLA,EQSTR,RTM,EPSC,IREFE,EX,RB,ALPHA,EPSU,
  55 :      2   EPO1,ICAL,IBB1,IGAU1,IDIR,PASDT)
  56 :          IF(IFISU.EQ.0.OR.IREFE.EQ.1) RETURN
  57 : C
  58 : C        TENSION-SOFTENING
  59 : C
  59 :          CALL FISSU(S1X,STRNX,TANG,RTM,PENT)
  60 :          RETURN
  61 :       ENDIF
  62 : C======================================================================
  63 :       IF(IFISU.EQ.0.AND.JFISU.EQ.1) THEN
  64 : 
  65 : C       COMPORTEMENT DU BETON FISSURE (FISSURE FERMEE)
  66 : C
  66 :          CALL BETFIS(EPSRX,STRNX,STRNRX,SIGMRX,SIGRX,TANG,S1X,EDT,EDC,
  67 :      1   EPST,EPSC,DEFR,RTM,IFISU,JFISU,IPLA,EQSTR,PENT,EX,RB,ALPHA,
  68 :      2   EMAX,EPO1,ICAL,IBB1,IGAU1,IDIR)
  69 :         RETURN
  70 :       ENDIF
  71 : C==========================================================================
  72 : C
  73 : C       BETON INTEGRE (NON FISSURE)
  74 : C
  75 : C==========================================================================
  76 :       IF(DABS(STRNRX).LT.EPSEQ) THEN
  77 : C
  78 : C       ************ DOMAINE DES DECHARGES ************
  79 : C                    (DECHARGE ENDOMMAGEE)
  80 : C   ZONES DE COMPRESSION OU ZONES DE TRACTION APRES DECHARGE EN COMPRESSION
  81 : C
  82 :         PASDC= 'PLAEND22'
  82 :         CALL DECHAC(S1X,EDC,EPST,EDT,RTM,EPSC,DEFR,SIGRX,STRNX,
  83 :      1  STRNRX,EPSRX,IFISU,JFISU,SIGMRX,IPLA,TANG,PENT,EX,RB,ALPHA,
  84 :      2  EPO1,IBB1,IGAU1,IDIR,PASDC)
  85 :       RETURN
  86 :       ENDIF
  87 : C==========================================================================
  88 : 
  89 :       IF(DABS(STRNRX).GE.EPSEQ .AND. SIGMRX .LT. 0.D0) THEN
  90 : C
  91 : C       ************ ECROUISSAGE EN COMPRESSION ************
  92 : C
  92 :       CALL CDCY(EPEQ,SEQ,IPLA,TANG,EX,RB,EMAX,EPO1,ICAL)
  93 :       EQSTR=SEQ
  94 :       EPSEQ=EPEQ
  95 :       S1X=-SEQ
  96 :       RETURN
  97 :       ENDIF
  98 : C==========================================================================
  99 :       IF(DABS(STRNRX).GE.EPSEQ.AND.SIGMRX.GE.0.D0.AND.EPSRX.GE.0.D0)
 100 :      + THEN
 101 : C
 102 : C       ***********  DOMAINE DES TRACTIONS  ***********
 103 : C       ********** POINT INITIALLEMENT TENDU **********
 104 : C
 105 :        IF(SIGMRX.GE.RTM) THEN
 106 : C
 107 : C              *************************************
 108 : C              *    CE POINT VIENT DE FISSURER     *
 109 : C              *************************************
 110 : C
 110 :           CALL FISSU2(S1X,SIGRX,STRNX,RTM,PENT,TANG)
 111 :           IFISU=1
 112 :           JFISU=1
 113 : C
 114 : C       Remarque : … l it‚ration suivante chemin ---> BETFIS
 115 : C
 116 :        ELSE
 117 : C              ****************************************
 118 : C              *  POINT ELASTIQUE INTEGRE (TRACTION)  *
 119 : C              ****************************************
 120 : C
 121 :           S1X=SIGMRX
 122 :        ENDIF
 123 :       RETURN
 124 :       ENDIF
 125 : C==========================================================================
 126 :       IF(DABS(STRNRX).GE.EPSEQ.AND.SIGMRX.GE.0.D0.AND.EPSRX.LT.0.D0)
 127 :      + THEN
 128 : C
 129 : C       *********** DOMAINE DES TRACTIONS  ***********
 130 : C       ******** POINT INITIALLEMENT COMPRIME **********
 131 : C
 132 :         PASDC= 'PLAEND23'
 132 :         CALL DECHAC(S1X,EDC,EPST,EDT,RTM,EPSC,DEFR,SIGRX,STRNX,
 133 :      1  STRNRX,EPSRX,IFISU,JFISU,SIGMRX,IPLA,TANG,PENT,EX,RB,ALPHA,
 134 :      2  EPO1,IBB1,IGAU1,IDIR,PASDC)
 135 :       RETURN
 136 :       ENDIF
 137 : C--------------------------------------------------------------------------
 138 :       END
 139 : 
 140 :  
 141 :