1 : C NEWBE2 SOURCE CHAT 12/04/06 21:15:21 7348 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 : C 54 : C**************----------------------------------------------- 55 : C SLOPI -> ELASTIC YOUNG MODULUS 56 : C ET___ -> STRAIN AT THE PICK 57 : C**************----------------------------------------------- 58 : SLOPI = 2.D0*FC___/EZERO 59 : C 60 : EPTRA = ET___/EZERO 61 : STRET = FLT__/FC___ 62 : STREL = FTPL_/FC___ 63 : C 64 : SCALE = EZERO/FC___ 65 : C 66 : IF (EPMAX .GE. 0.D0 .AND. EPMIN .LE. 0.D0) THEN 67 : EPMIN = EPTRA 68 : STMIN = STRET 69 : ENDIF 70 : C**************----------------------------------------------- 71 : C SLOPE -> DISCHARGE YOUNG MODULUS 72 : C STMAX -> LAST COMPRESSION STRESS ON THE ENVELOPE 73 : C EPLA1 -> PLASTIC STRAIN FOR THE FIRST DISCHARGE 74 : C**************----------------------------------------------- 75 : SLOPE = SLOPI*(1.D0-EPMAX**2/(1.D0-EPMAX+EPMAX**2))*SCALE 76 : C 77 : IF (EPMAX .GE. -1.D0) THEN 78 : STMAX = EPMAX*(2.D0+EPMAX) 79 : STRX2 = 0.92D0*STMAX 80 : ELSE IF (EPMAX .LE. EPS85) THEN 81 : STMAX = -1.D0-ZETA_*EZERO*(EPS85+1.D0) 82 : STRX2 = STMAX+0.08D0 83 : ELSE 84 : STMAX = -1.D0-ZETA_*EZERO*(EPMAX+1.D0) 85 : STRX2 = STMAX+0.08D0 86 : ENDIF 87 : C 88 : IF (FACX2 .GT. 0.D0) STRX2 = STMAX/FACX2 89 : EPLA1 = EPMAX-STMAX/SLOPE 90 : C 91 : IF (STRX2 .GT. 0.D0) STRX2 = 0.D0 92 : EPSX2 = EPMAX-(STMAX-STRX2)/SLOPE 93 : C 94 : TANGD = SLOPE 95 : C 96 : IF (EPSON .LE. EPMAX) THEN 97 : GOTO 100 98 : ELSEIF (EPSON .GE. EPMIN) THEN 99 : IF (STMIN .GT. 1.D0) GOTO 300 100 : C 101 : GOTO 200 102 : ELSE 103 : GOTO 300 104 : ENDIF 105 : C 106 : C ####################################################### 107 : C 108 : C ...... LOADING IN COMPRESSION ALONG THE ENVELOPE ...... 109 : C 110 : C ####################################################### 111 : C 112 : 100 EPMAX = EPSON 113 : C 114 : IF (EPMAX .GE. -1.D0) THEN 115 : C**************----------------------------------------------- 116 : C ASCENDING CURVE 117 : C**************----------------------------------------------- 118 : STRES = EPMAX*(2.D0+EPMAX) 119 : TANGE = 2.D0*(FC___/EZERO)*(1.D0+EPMAX) 120 : C 121 : STRX3 = 0.92D0*STRES 122 : ELSE IF (EPMAX .LE. EPS85) THEN 123 : C**************----------------------------------------------- 124 : C PLATEAU IN COMPRESSION 125 : C**************----------------------------------------------- 126 : STRES = -1.D0-ZETA_*EZERO*(EPS85+1.D0) 127 : TANGE = 0.D0 128 : C 129 : STRX3 = STRES+0.08D0 130 : ELSE 131 : C**************----------------------------------------------- 132 : C DESCENDING CURVE 133 : C**************----------------------------------------------- 134 : STRES = -1.D0-ZETA_*EZERO*(EPMAX+1.D0) 135 : TANGE = -ZETA_*FC___ 136 : C 137 : STRX3 = STRES+0.08D0 138 : ENDIF 139 : C 140 : STMAX = STRES 141 : IF (FACX2 .GT. 0.D0) STRX3 = STMAX/FACX2 142 : C 143 : SLOPE = SLOPI*(1.D0-EPMAX**2/(1.D0-EPMAX+EPMAX**2))*SCALE 144 : EPLA2 = EPMAX-STMAX/SLOPE 145 : C 146 : IF (STMIN .GT. 0.D0) THEN 147 : C**************----------------------------------------------- 148 : C DAMAGE IN TENSION 149 : C**************----------------------------------------------- 150 : IF (EPLA2 .GT. -1.D0) THEN 151 : STMI1 = (1.D0+EPLA2)/(1.D0+EPLA1)*STMIN 152 : ELSE 153 : STMI1 = 0.D0 154 : ENDIF 155 : C 156 : IF (STRX3 .GT. 0.D0) STRX2 = 0.D0 157 : EPSX3 = EPMAX-(STMAX-STRX3)/SLOPE 158 : C 159 : EPMI1 = STMI1/STMIN*(EPMIN-EPLA1)+EPLA2 160 : EPMIN = (STMI1-STRX3)/(STMIN-STRX2)*(EPMIN-EPSX2)+EPSX3 161 : C 162 : IF (EPMI1 .LT. EPMIN) EPMIN = EPMI1 163 : C 164 : STMIN = STMI1 165 : C 166 : SLOPX = STMIN/(EPMIN-EPLA2) 167 : IF (SLOPX .GT. SLOPE) EPMIN = STMIN/SLOPE+EPLA2 168 : ELSE 169 : EPMIN = EPLA2 170 : ENDIF 171 : C 172 : GO TO 2000 173 : C 174 : C ####################################################### 175 : C 176 : C ........ LOADING IN TENSION ALONG THE ENVELOPE ........ 177 : C 178 : C ####################################################### 179 : C 180 : 200 TANGE = SLOPI/(1.D0-TRAFA) 181 : STRES = (TANGE*SCALE)*(EPSON-EPMIN)+STMIN 182 : C 183 : IF (EPLA1 .LE. -1.D0) THEN 184 : C**************----------------------------------------------- 185 : C NO MORE TENSILE FORCES ALLOWED 186 : C**************----------------------------------------------- 187 : STRES = 0.D0 188 : TANGE = 0.D0 189 : ELSEIF (STRES .LT. STREL) THEN 190 : C**************----------------------------------------------- 191 : C PLATEAU FOR TENSILE FORCES 192 : C**************----------------------------------------------- 193 : IF (STMIN .LT. STREL) THEN 194 : STRES = STMIN 195 : ELSE 196 : STRES = STREL 197 : ENDIF 198 : C 199 : TANGE = 0.D0 200 : ENDIF 201 : C 202 : STMIN = STRES 203 : EPMIN = EPSON 204 : C 205 : GOTO 2000 206 : C 207 : C ####################################################### 208 : C 209 : C ....... LOADING OR UNLOADING UNDER THE ENVELOPE ....... 210 : C 211 : C ####################################################### 212 : C 213 : 300 RFACT = 0.1D0+0.9D0/(1.D0+(EPMIN-EPLA1)) 214 : IF (FACLO .GT. 0.D0) RFACT = FACLO 215 : C-JP- 216 : IF (STMIN .GT. 1.D0) THEN 217 : STRES = SLOPE*(EPSON-EPSO1)+STRE1 218 : TANGE = SLOPE/SCALE 219 : C-JP- 220 : IF ((STRET .LT. 1.D0) .AND. (STRES .GT. 0.D0)) THEN 221 : C-JP- EPMIN = EPTRA 222 : C-JP- STMIN = SLOPE*(EPMIN-EPSO1)+STRE1 223 : STMIN = STRET 224 : EPMIN = (STMIN-STRE1)/SLOPE+EPSO1 225 : C 226 : IF (EPSON .LT. EPMIN) THEN 227 : EPMIN = EPSON 228 : STMIN = SLOPE*(EPMIN-EPSO1)+STRE1 229 : ENDIF 230 : C 231 : GOTO 200 232 : ENDIF 233 : C-JP- 234 : GOTO 2000 235 : ENDIF 236 : C-JP- 237 : IF (EPSON .LE. EPSX2) THEN 238 : STRES = SLOPE*(EPSON-EPSX2)+STRX2 239 : TANGE = SLOPE/SCALE 240 : C 241 : GOTO 2000 242 : ELSEIF (EPSON .GT. EPSO1) THEN 243 : RFACT = RFACT/2.D0 244 : ENDIF 245 : C 246 : IF (STMIN .LE. 0.D0) THEN 247 : EPMI1 = EPLA1+RFACT*(EPMIN-EPLA1) 248 : ELSE 249 : EPMI1 = EPMIN 250 : ENDIF 251 : C 252 : STRES = TANGD*(EPSON-EPSO1)+STRE1 253 : TANGE = TANGD/SCALE 254 : C 255 : IF (EPSON .GT. EPMI1) THEN 256 : IF (EPSON .GT. EPSO1) THEN 257 : IF (STRES .GT. STMIN) THEN 258 : IF (STMIN .GT. 0.D0) THEN 259 : EPMIN = (STMIN-STRE1)/TANGD+EPSO1 260 : C 261 : GOTO 200 262 : ELSE 263 : STRES = 0.D0 264 : TANGE = 0.D0 265 : ENDIF 266 : ENDIF 267 : ELSE 268 : STRES = 0.D0 269 : TANGE = 0.D0 270 : ENDIF 271 : C 272 : GOTO 2000 273 : ENDIF 274 : C**************----------------------------------------------- 275 : C TANGENT MODULUS AT "EPSON=EPMI1" EQUAL TO "SLOP2" 276 : C WHEN "FACT1 -> +00", "SLOP2 -> 0" 277 : C WHEN "FACT1 -> 1+", THE CURVED LINE BECOMES 278 : C A STRAIGHT LINE 279 : C TANGENT MODULUS AT "EPSON=EPSX2" EQUAL TO "SLOP3" 280 : C WHEN "FACT2 -> +00", "SLOP2 -> SLOPE" 281 : C WHEN "FACT2 -> 1+", THE CURVED LINE BECOMES 282 : C A STRAIGHT LINE 283 : C NOTE: "FACT1" AND "FACT2" MUST BE ALWAYS ">1" 284 : C**************----------------------------------------------- 285 : C TANGENT MODULUS FOR THE CHARGE CURVE: 286 : C SLOP2 -> AT THE POINT EPSON = EPMI1 287 : C SLOP3 -> AT THE POINT EPSON = EPSX2 288 : C**************----------------------------------------------- 289 : SLOP1 = (STRX2-STMIN)/(EPSX2-EPMI1) 290 : C 291 : IF (EPMI1 .GT. EPLA1) THEN 292 : SLOP2 = STMIN/(0.90D0*(EPMI1-EPLA1)) 293 : ELSE 294 : SLOP2 = SLOPE 295 : ENDIF 296 : C 297 : SLOP3 = -1.D0*STRX2/(0.10D0*(EPMI1-EPLA1)+(EPLA1-EPSX2)) 298 : C 299 : SLOP4 = STRX2/FACT1/(EPSX2-EPMI1) 300 : SLOP5 = STRX2*FACT2/(EPSX2-EPMI1) 301 : C 302 : IF (SLOP2 .LT. SLOP4) SLOP2 = SLOP4 303 : IF (SLOP3 .GT. SLOP5) SLOP3 = SLOP5 304 : C 305 : IF (SLOP2 .GE. SLOP1) THEN 306 : STRES = SLOP1*(EPSON-EPSO1)+STRE1 307 : TANGE = SLOP1/SCALE 308 : C 309 : GOTO 2000 310 : ENDIF 311 : C 312 : STRX4 = (STRX2-STMIN)-SLOP2*(EPSX2-EPMI1) 313 : C 314 : FACTT = (SLOP3-SLOP2)*(EPSX2-EPMI1)/STRX4 315 : STREE = STRX4*(((EPSON-EPMI1)/(EPSX2-EPMI1))**FACTT)+ 316 : . SLOP2*(EPSON-EPMI1)+STMIN 317 : TANGG = (SLOP3-SLOP2)*(((EPSON-EPMI1)/ 318 : . (EPSX2-EPMI1))**(FACTT-1.D0))+SLOP2 319 : C 320 : C**************----------------------------------------------- 321 : C CHECK THE ACTIVE CURVE 322 : C**************----------------------------------------------- 323 : IF (EPSON .GT. EPSO1) THEN 324 : IF (STREE .LT. STRES) THEN 325 : C**************----------------------------------------------- 326 : C PARABOLIC CURVE 327 : C**************----------------------------------------------- 328 : STRES = STREE 329 : TANGE = TANGG/SCALE 330 : ENDIF 331 : ELSE 332 : IF (STREE .GT. STRES) THEN 333 : C**************----------------------------------------------- 334 : C PARABOLIC CURVE 335 : C**************----------------------------------------------- 336 : STRES = STREE 337 : TANGE = TANGG/SCALE 338 : ENDIF 339 : ENDIF 340 : C 341 : 2000 RETURN 342 : END 343 : 344 : 345 :
© Cast3M 2003 - All rights reserved.
Disclaimer