Télécharger la plaquette Cast3M


fr en

Télécharger newbe2.eso

Retour à la liste

Numérotation des lignes : 

newbe2

  1. C NEWBE2    SOURCE    PV        14/03/18    21:15:01     8001
  2. C NEWBE2    SOURCE    AM        94/02/28    21:33:04
  3.       SUBROUTINE NEWBE2
  4.      .       ( FC___, EZERO, EPS85, ZETA_, FLT__, FTPL_,
  5.      .         TRAFA, EPSON, EPSO1, EPMAX, EPMIN, STRES,
  6.      .         STRE1, STMIN, TANGE, TANGD, FACT1, FACT2,
  7.      .         FACX2, FACLO, EPLA1, ET___ )
  8. *
  9. *---------------------------------------------------------------------
  10. *
  11. *   CONCRETE BEHAVIOUR MODEL - THE BEST !!!
  12. *                              ------------
  13. *   (WITH TENSION ALLOWED AND COMPRESSION PLATEAU AFTER SOFTENING)
  14. ***** CLEANED *****
  15. *   (CRACK CLOSING MODEL CONSIDERED) 22/08/94
  16. *   (LINEAR BEHAVIOUR CURVE FOR SOFTENING IN TENSION)
  17. *
  18. *---------------------------------------------------------------------
  19. *
  20. *   INPUT :
  21. *   -----
  22. *               FC___ : MAXIMUM COMPRESSION STRESS
  23. *               EZERO : STRAIN OF THE MAXIMUM COMPRESSION STRESS
  24. *               FLT__ : MAXIMUM TENSILE STRESS
  25. *               ZETA_ : SLOPE OF THE DESCENDING PART OF THE CURVE
  26. *               EPSON : ACTUAL STRAIN VALUE
  27. *               EPSO1 : LAST STRAIN VALUE
  28. *               EPMAX : GRATEST COMPRES. STRAIN VALUE EVER REACHED
  29. *-JP-
  30. *               EPS85 : PLATEAU STRAIN AFTER SOFTENING
  31. *-JP-
  32. *-JP2-
  33. *               STRE1 : LAST STRESS VALUE
  34. *               TRAFA : DESCENDING PART OF THE TRACTION CURVE FACTOR
  35. *-JP2-
  36. *
  37. *   OUTPUT :
  38. *   ------
  39. *               EPMIN : GRATEST TENSILE STRAIN VALUE EVER REACHED
  40. *               TANGE : ACTUAL TANGENT MODULUS
  41. *               STRES : ACTUAL STRESS VALUE
  42. *
  43. *       NOTE : POSITIVE STRESSES = TENSILE STRESSES
  44. *              NEGATIVE STRESSES = COMPRESSION STRESSES
  45. *
  46. *              POSITIVE STRAIN = TENSILE STRAIN
  47. *              NEGATIVE STRAIN = COMPRESSION STRAIN
  48. *
  49. *---------------------------------------------------------------------
  50. C
  51.       IMPLICIT INTEGER(I-N)
  52.         IMPLICIT REAL*8 (A-H,O-Z)
  53. -INC CCREEL
  54. C
  55. C**************-----------------------------------------------
  56. C       SLOPI -> ELASTIC YOUNG MODULUS
  57. C       ET___ -> STRAIN AT THE PICK
  58. C**************-----------------------------------------------
  59.         SLOPI = 2.D0*FC___/EZERO
  60. C
  61.         EPTRA = ET___/EZERO
  62.         STRET = FLT__/FC___
  63.         STREL = FTPL_/FC___
  64. C
  65.         SCALE = EZERO/FC___
  66. C
  67.         IF (EPMAX .GE. 0.D0 .AND. EPMIN .LE. 0.D0) THEN
  68.            EPMIN = EPTRA
  69.            STMIN = STRET
  70.         ENDIF
  71. C**************-----------------------------------------------
  72. C       SLOPE -> DISCHARGE YOUNG MODULUS
  73. C       STMAX -> LAST COMPRESSION STRESS ON THE ENVELOPE
  74. C       EPLA1 -> PLASTIC STRAIN FOR THE FIRST DISCHARGE
  75. C**************-----------------------------------------------
  76.         SLOPE = SLOPI*(1.D0-EPMAX**2/(1.D0-EPMAX+EPMAX**2))*SCALE
  77. C
  78.         IF (EPMAX .GE. -1.D0) THEN
  79.            STMAX = EPMAX*(2.D0+EPMAX)
  80.            STRX2 = 0.92D0*STMAX
  81.         ELSE IF (EPMAX .LE. EPS85) THEN
  82.            STMAX = -1.D0-ZETA_*EZERO*(EPS85+1.D0)
  83.            STRX2 = STMAX+0.08D0
  84.         ELSE
  85.            STMAX = -1.D0-ZETA_*EZERO*(EPMAX+1.D0)
  86.            STRX2 = STMAX+0.08D0
  87.         ENDIF
  88. C
  89.         IF (FACX2 .GT. 0.D0) STRX2 = STMAX/FACX2
  90.         EPLA1 = EPMAX-STMAX/SLOPE
  91. C
  92.         IF (STRX2 .GT. 0.D0) STRX2 = 0.D0
  93.         EPSX2 = EPMAX-(STMAX-STRX2)/SLOPE
  94. C
  95.         TANGD = SLOPE
  96. C
  97.         IF (EPSON .LE. EPMAX) THEN
  98.            GOTO 100
  99.         ELSEIF (EPSON .GE. EPMIN) THEN
  100.            IF (STMIN .GT. 1.D0) GOTO 300
  101. C
  102.            GOTO 200
  103.         ELSE
  104.            GOTO 300
  105.         ENDIF
  106. C
  107. C       #######################################################
  108. C
  109. C       ...... LOADING IN COMPRESSION ALONG THE ENVELOPE ......
  110. C
  111. C       #######################################################
  112. C
  113.   100   EPMAX = EPSON
  114. C
  115.         IF (EPMAX .GE. -1.D0) THEN
  116. C**************-----------------------------------------------
  117. C          ASCENDING CURVE
  118. C**************-----------------------------------------------
  119.            STRES = EPMAX*(2.D0+EPMAX)
  120.            TANGE = 2.D0*(FC___/EZERO)*(1.D0+EPMAX)
  121. C
  122.            STRX3 = 0.92D0*STRES
  123.         ELSE IF (EPMAX .LE. EPS85) THEN
  124. C**************-----------------------------------------------
  125. C          PLATEAU IN COMPRESSION
  126. C**************-----------------------------------------------
  127.            STRES = -1.D0-ZETA_*EZERO*(EPS85+1.D0)
  128.            TANGE = 0.D0
  129. C
  130.            STRX3 = STRES+0.08D0
  131.         ELSE
  132. C**************-----------------------------------------------
  133. C          DESCENDING CURVE
  134. C**************-----------------------------------------------
  135.            STRES = -1.D0-ZETA_*EZERO*(EPMAX+1.D0)
  136.            TANGE = -ZETA_*FC___
  137. C
  138.            STRX3 = STRES+0.08D0
  139.         ENDIF
  140. C
  141.         STMAX = STRES
  142.         IF (FACX2 .GT. 0.D0) STRX3 = STMAX/FACX2
  143. C
  144.         SLOPE = SLOPI*(1.D0-EPMAX**2/(1.D0-EPMAX+EPMAX**2))*SCALE
  145.         EPLA2 = EPMAX-STMAX/SLOPE
  146. C
  147.         IF (STMIN .GT. 0.D0) THEN
  148. C**************-----------------------------------------------
  149. C          DAMAGE IN TENSION
  150. C**************-----------------------------------------------
  151.            IF (EPLA2 .GT. -1.D0) THEN
  152.               STMI1 = (1.D0+EPLA2)/(1.D0+EPLA1)*STMIN
  153.            ELSE
  154.               STMI1 = 0.D0
  155.            ENDIF
  156. C
  157.            IF (STRX3 .GT. 0.D0) STRX2 = 0.D0
  158.            EPSX3 = EPMAX-(STMAX-STRX3)/SLOPE
  159. C
  160.            EPMI1 = STMI1/STMIN*(EPMIN-EPLA1)+EPLA2
  161.            EPMIN = (STMI1-STRX3)/(STMIN-STRX2)*(EPMIN-EPSX2)+EPSX3
  162. C
  163.            IF (EPMI1 .LT. EPMIN) EPMIN = EPMI1
  164. C
  165.            STMIN = STMI1
  166. C
  167.            SLOPX = STMIN/(EPMIN-EPLA2)
  168.            IF (SLOPX .GT. SLOPE) EPMIN = STMIN/SLOPE+EPLA2
  169.         ELSE
  170.            EPMIN = EPLA2
  171.         ENDIF
  172. C
  173.         GO TO 2000
  174. C
  175. C       #######################################################
  176. C
  177. C       ........ LOADING IN TENSION ALONG THE ENVELOPE ........
  178. C
  179. C       #######################################################
  180. C
  181.   200   TANGE = SLOPI/(1.D0-TRAFA)
  182.         STRES = (TANGE*SCALE)*(EPSON-EPMIN)+STMIN
  183. C
  184.         IF (EPLA1 .LE. -1.D0) THEN
  185. C**************-----------------------------------------------
  186. C          NO MORE TENSILE FORCES ALLOWED
  187. C**************-----------------------------------------------
  188.            STRES = 0.D0
  189.            TANGE = 0.D0
  190.         ELSEIF (STRES .LT. STREL) THEN
  191. C**************-----------------------------------------------
  192. C          PLATEAU FOR TENSILE FORCES
  193. C**************-----------------------------------------------
  194.            IF (STMIN .LT. STREL) THEN
  195.               STRES = STMIN
  196.            ELSE
  197.               STRES = STREL
  198.            ENDIF
  199. C
  200.            TANGE = 0.D0
  201.         ENDIF
  202. C
  203.         STMIN = STRES
  204.         EPMIN = EPSON
  205. C
  206.         GOTO 2000
  207. C
  208. C       #######################################################
  209. C
  210. C       ....... LOADING OR UNLOADING UNDER THE ENVELOPE .......
  211. C
  212. C       #######################################################
  213. C
  214.   300   RFACT = 0.1D0+0.9D0/(1.D0+(EPMIN-EPLA1))
  215.         IF (FACLO .GT. 0.D0) RFACT = FACLO
  216. C-JP-
  217.         IF (STMIN .GT. 1.D0) THEN
  218.            STRES = SLOPE*(EPSON-EPSO1)+STRE1
  219.            TANGE = SLOPE/SCALE
  220. C-JP-
  221.            IF ((STRET .LT. 1.D0) .AND. (STRES .GT. 0.D0)) THEN
  222. C-JP-              EPMIN = EPTRA
  223. C-JP-              STMIN = SLOPE*(EPMIN-EPSO1)+STRE1
  224.               STMIN = STRET
  225.               EPMIN = (STMIN-STRE1)/SLOPE+EPSO1
  226. C
  227.               IF (EPSON .LT. EPMIN) THEN
  228.                  EPMIN = EPSON
  229.                  STMIN = SLOPE*(EPMIN-EPSO1)+STRE1
  230.               ENDIF
  231. C
  232.               GOTO 200
  233.            ENDIF
  234. C-JP-
  235.            GOTO 2000
  236.         ENDIF
  237. C-JP-
  238.         IF (EPSON .LE. EPSX2) THEN
  239.            STRES = SLOPE*(EPSON-EPSX2)+STRX2
  240.            TANGE = SLOPE/SCALE
  241. C
  242.            GOTO 2000
  243.         ELSEIF (EPSON .GT. EPSO1) THEN
  244.            RFACT = RFACT/2.D0
  245.         ENDIF
  246. C
  247.         IF (STMIN .LE. 0.D0) THEN
  248.            EPMI1 = EPLA1+RFACT*(EPMIN-EPLA1)
  249.         ELSE
  250.            EPMI1 = EPMIN
  251.         ENDIF
  252. C
  253.         STRES = TANGD*(EPSON-EPSO1)+STRE1
  254.         TANGE = TANGD/SCALE
  255. C
  256.         IF (EPSON .GT. EPMI1) THEN
  257.            IF (EPSON .GT. EPSO1) THEN
  258.               IF (STRES .GT. STMIN) THEN
  259.                  IF (STMIN .GT. 0.D0) THEN
  260.                     EPMIN = (STMIN-STRE1)/TANGD+EPSO1
  261. C
  262.                     GOTO 200
  263.                  ELSE
  264.                     STRES = 0.D0
  265.                     TANGE = 0.D0
  266.                  ENDIF
  267.               ENDIF
  268.            ELSE
  269.               STRES = 0.D0
  270.               TANGE = 0.D0
  271.            ENDIF
  272. C
  273.            GOTO 2000
  274.         ENDIF
  275. C**************-----------------------------------------------
  276. C       TANGENT MODULUS AT "EPSON=EPMI1" EQUAL TO "SLOP2"
  277. C                 WHEN "FACT1 -> +00", "SLOP2 -> 0"
  278. C                 WHEN "FACT1 ->  1+", THE CURVED LINE BECOMES
  279. C                                             A STRAIGHT LINE
  280. C       TANGENT MODULUS AT "EPSON=EPSX2" EQUAL TO "SLOP3"
  281. C                 WHEN "FACT2 -> +00", "SLOP2 -> SLOPE"
  282. C                 WHEN "FACT2 ->  1+", THE CURVED LINE BECOMES
  283. C                                             A STRAIGHT LINE
  284. C       NOTE: "FACT1" AND "FACT2" MUST BE ALWAYS ">1"
  285. C**************-----------------------------------------------
  286. C       TANGENT MODULUS FOR THE CHARGE CURVE:
  287. C          SLOP2 -> AT THE POINT EPSON = EPMI1
  288. C          SLOP3 -> AT THE POINT EPSON = EPSX2
  289. C**************-----------------------------------------------
  290.         SLOP1 = (STRX2-STMIN)/(EPSX2-EPMI1)
  291. C
  292.         IF (EPMI1 .GT. EPLA1) THEN
  293.            SLOP2 = STMIN/(0.90D0*(EPMI1-EPLA1))
  294.         ELSE
  295.            SLOP2 = SLOPE
  296.         ENDIF
  297. C
  298.         SLOP3 = -1.D0*STRX2/(0.10D0*(EPMI1-EPLA1)+(EPLA1-EPSX2))
  299. C
  300.         SLOP4 = STRX2/FACT1/(EPSX2-EPMI1)
  301.         SLOP5 = STRX2*FACT2/(EPSX2-EPMI1)
  302. C
  303.         IF (SLOP2 .LT. SLOP4) SLOP2 = SLOP4
  304.         IF (SLOP3 .GT. SLOP5) SLOP3 = SLOP5
  305. C
  306.         IF (SLOP2 .GE. SLOP1*(1.d0 - xzprec)) THEN
  307.            STRES = SLOP1*(EPSON-EPSO1)+STRE1
  308.            TANGE = SLOP1/SCALE
  309. C
  310.            GOTO 2000
  311.         ENDIF
  312. C
  313.         STRX4 = (STRX2-STMIN)-SLOP2*(EPSX2-EPMI1)
  314. C
  315.         FACTT = (SLOP3-SLOP2)*(EPSX2-EPMI1)/STRX4
  316.         STREE = STRX4*(((EPSON-EPMI1)/(EPSX2-EPMI1))**FACTT)+
  317.      .          SLOP2*(EPSON-EPMI1)+STMIN
  318.         TANGG = (SLOP3-SLOP2)*(((EPSON-EPMI1)/
  319.      .          (EPSX2-EPMI1))**(FACTT-1.D0))+SLOP2
  320. C
  321. C**************-----------------------------------------------
  322. C       CHECK THE ACTIVE CURVE
  323. C**************-----------------------------------------------
  324.         IF (EPSON .GT. EPSO1) THEN
  325.            IF (STREE .LT. STRES) THEN
  326. C**************-----------------------------------------------
  327. C             PARABOLIC CURVE
  328. C**************-----------------------------------------------
  329.               STRES = STREE
  330.               TANGE = TANGG/SCALE
  331.            ENDIF
  332.         ELSE
  333.            IF (STREE .GT. STRES) THEN
  334. C**************-----------------------------------------------
  335. C             PARABOLIC CURVE
  336. C**************-----------------------------------------------
  337.               STRES = STREE
  338.               TANGE = TANGG/SCALE
  339.            ENDIF
  340.         ENDIF
  341. C
  342.  2000   RETURN
  343.         END
  344.  
  345.  
  346.  
  347.  
  348.  

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