Télécharger la plaquette Cast3M


fr en

Télécharger plaend.eso

Retour à la liste

Numérotation des lignes : 

plaend

  1. C PLAEND    SOURCE    CB215821  16/04/21    21:18:01     8920
  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=ABS(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'
  54.          CALL DECHAT(EPSRX,STRNX,STRNRX,SIGMRX,TANG,S1X,EDT,EDC,
  55.      1   EPST,DEFR,IFISU,IPLA,EQSTR,RTM,EPSC,IREFE,EX,RB,ALPHA,EPSU,
  56.      2   EPO1,ICAL,IBB1,IGAU1,IDIR,PASDT)
  57.          IF(IFISU.EQ.0.OR.IREFE.EQ.1) RETURN
  58. C
  59. C        TENSION-SOFTENING
  60. C
  61.          CALL FISSU(S1X,STRNX,TANG,RTM,PENT)
  62.          RETURN
  63.       ENDIF
  64. C======================================================================
  65.       IF(IFISU.EQ.0.AND.JFISU.EQ.1) THEN
  66.  
  67. C       COMPORTEMENT DU BETON FISSURE (FISSURE FERMEE)
  68. C
  69.          CALL BETFIS(EPSRX,STRNX,STRNRX,SIGMRX,SIGRX,TANG,S1X,EDT,EDC,
  70.      1   EPST,EPSC,DEFR,RTM,IFISU,JFISU,IPLA,EQSTR,PENT,EX,RB,ALPHA,
  71.      2   EMAX,EPO1,ICAL,IBB1,IGAU1,IDIR)
  72.         RETURN
  73.       ENDIF
  74. C==========================================================================
  75. C
  76. C       BETON INTEGRE (NON FISSURE)
  77. C
  78. C==========================================================================
  79.       IF(ABS(STRNRX).LT.EPSEQ) THEN
  80. C
  81. C       ************ DOMAINE DES DECHARGES ************
  82. C                    (DECHARGE ENDOMMAGEE)
  83. C   ZONES DE COMPRESSION OU ZONES DE TRACTION APRES DECHARGE EN COMPRESSION
  84. C
  85.         PASDC= 'PLAEND22'
  86.         CALL DECHAC(S1X,EDC,EPST,EDT,RTM,EPSC,DEFR,SIGRX,STRNX,
  87.      1  STRNRX,EPSRX,IFISU,JFISU,SIGMRX,IPLA,TANG,PENT,EX,RB,ALPHA,
  88.      2  EPO1,IBB1,IGAU1,IDIR,PASDC)
  89.       RETURN
  90.       ENDIF
  91. C==========================================================================
  92.  
  93.       IF(ABS(STRNRX).GE.EPSEQ .AND. SIGMRX .LT. 0.D0) THEN
  94. C
  95. C       ************ ECROUISSAGE EN COMPRESSION ************
  96. C
  97.       CALL CDCY(EPEQ,SEQ,IPLA,TANG,EX,RB,EMAX,EPO1,ICAL)
  98.       EQSTR=SEQ
  99.       EPSEQ=EPEQ
  100.       S1X=-SEQ
  101.       RETURN
  102.       ENDIF
  103. C==========================================================================
  104.       IF(ABS(STRNRX).GE.EPSEQ.AND.SIGMRX.GE.0.D0.AND.EPSRX.GE.0.D0)
  105.      + THEN
  106. C
  107. C       ***********  DOMAINE DES TRACTIONS  ***********
  108. C       ********** POINT INITIALLEMENT TENDU **********
  109. C
  110.        IF(SIGMRX.GE.RTM) THEN
  111. C
  112. C              *************************************
  113. C              *    CE POINT VIENT DE FISSURER     *
  114. C              *************************************
  115. C
  116.           CALL FISSU2(S1X,SIGRX,STRNX,RTM,PENT,TANG)
  117.           IFISU=1
  118.           JFISU=1
  119. C
  120. C       Remarque : … l it‚ration suivante chemin ---> BETFIS
  121. C
  122.        ELSE
  123. C              ****************************************
  124. C              *  POINT ELASTIQUE INTEGRE (TRACTION)  *
  125. C              ****************************************
  126. C
  127.           S1X=SIGMRX
  128.        ENDIF
  129.       RETURN
  130.       ENDIF
  131. C==========================================================================
  132.       IF(ABS(STRNRX).GE.EPSEQ.AND.SIGMRX.GE.0.D0.AND.EPSRX.LT.0.D0)
  133.      + THEN
  134. C
  135. C       *********** DOMAINE DES TRACTIONS  ***********
  136. C       ******** POINT INITIALLEMENT COMPRIME **********
  137. C
  138.         PASDC= 'PLAEND23'
  139.         CALL DECHAC(S1X,EDC,EPST,EDT,RTM,EPSC,DEFR,SIGRX,STRNX,
  140.      1  STRNRX,EPSRX,IFISU,JFISU,SIGMRX,IPLA,TANG,PENT,EX,RB,ALPHA,
  141.      2  EPO1,IBB1,IGAU1,IDIR,PASDC)
  142.       RETURN
  143.       ENDIF
  144. C--------------------------------------------------------------------------
  145.       END
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  

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