Télécharger plaend.eso

Retour à la liste

Numérotation des lignes :

  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