Télécharger newbe2.eso

Retour à la liste

Numérotation des lignes :

  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