Télécharger foncre.eso

Retour à la liste

Numérotation des lignes :

  1. C FONCRE SOURCE CB215821 16/04/21 21:16:56 8920
  2. C SOURCE BERTHOLLET
  3. C
  4. SUBROUTINE FONCRE(FLU,FI,TFLU,DUR,ELAS0,BETJEF,BETFLU)
  5. C
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. -INC CCREEL
  9. C
  10. C*****************************************************************
  11. C
  12. SEGMENT BETJEF
  13. REAL*8 AA,BETA,FC,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
  14. & TCON,DPSTF1,DPSTF2,TETA,PDT,TP00
  15. INTEGER ICT,ICC,IMOD,IVIS,ITER,
  16. & ISIM,IBB,IGAU,IZON
  17. ENDSEGMENT
  18. SEGMENT BETFLU
  19. REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2,
  20. & TP0,TZER
  21. INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR
  22. ENDSEGMENT
  23. C
  24. C MODELE B3 ET BP-KX
  25. C*****************************************************************
  26. C
  27. CIMENT = 266.
  28. EAU = 161.
  29. AGREG = 1615.
  30. SABLE = 1133.
  31. RH = 50.
  32. C POUR LE CAS D UNE EPROUVETTE 16*32
  33. RAYON = 0.08
  34. HAUT = 0.32
  35. C U = PERIMETRE EN CONTACT AVEC L EXTERIEUR EN M
  36. U = 2*XPI*RAYON
  37. C
  38. C AC = SURFACE EN CONTACT AVEC L ATMOSPHERE EN M2
  39. AC = XPI*(RAYON**2)
  40. C
  41. C VOL = VOLUME DE L ELEMENT CONSIDERE EN M3
  42. VOL = XPI*(RAYON**2)*HAUT
  43. C
  44. SILICE = 40.3
  45. FINES = 57.0
  46. ADJUV = 9.98
  47. C
  48. C*****************************************************************
  49. C CHOIX DU MODELE DE FLUAGE
  50. C*****************************************************************
  51. C
  52. IF (ITYPE.LT.1.OR.ITYPE.GT.8) THEN
  53. WRITE(*,*)'MAUVAIS CHOIX DE FONCTION DE FLUAGE'
  54. GO TO 100
  55. ENDIF
  56. C
  57. C 1 / MODELE ACI209 :FLUAGE TOTAL
  58. IF (ITYPE.EQ.1) GO TO 10
  59. C
  60. C 2 / MODELE ESPION ET HALLEUX : FLUAGE TOTAL
  61. IF (ITYPE.EQ.2) GO TO 20
  62. C
  63. C 3 / MODELE B3 :FLUAGE PROPRE ET DE DESSICCATION
  64. IF (ITYPE.EQ.3) GO TO 30
  65. C
  66. C 4 / MODELE CEB FIP : BETON GENERAL
  67. IF (ITYPE.EQ.4) GO TO 40
  68. C
  69. C 5 / MODELE EUROCODE 2 : BETON GENERAL
  70. IF (ITYPE.EQ.5) GO TO 40
  71. C
  72. C 6 / MODELE BEPL91 : BHP
  73. IF (ITYPE.EQ.6) GO TO 40
  74. C
  75. C 7 / MODELE AFREM :BHP : FLUAGE PROPRE ET DE DESSICCATION
  76. IF (ITYPE.EQ.7) GO TO 50
  77. C
  78. C 8 / MODELE BP-KX :BETON ORDINAIRE
  79. IF (ITYPE.EQ.8) GO TO 60
  80. C
  81. C*****************************************************************
  82. C 1 : MODELE ACI209
  83. C*****************************************************************
  84. C*****************************************************************
  85. C
  86. 10 E28=4734*SQRT(FC)
  87. FLU = 1/E28*SQRT((4+0.85*TFLU)/TFLU)*
  88. *(1+2.35*(TFLU**(-0.118))*(DUR**(0.6))/(10+
  89. *(DUR**(0.6))))
  90. ELAS0 = E28*SQRT(TFLU/(4+0.85*TFLU))
  91. TFLU = 0.D0
  92. DUR = 0.D0
  93. C
  94. FI = FLU*ELAS0-1
  95. GO TO 100
  96. C
  97. C*****************************************************************
  98. C 2 : MODELE ESPION ET HALLEUX
  99. C*****************************************************************
  100. C*****************************************************************
  101. C
  102. 20 FLU = 1/26000*(1+2.867*(TFLU**(-0.118))*
  103. *(DUR**(0.6))/(15+(DUR**(0.6))))
  104. C
  105. E28=4734*SQRT(FC)
  106. ELAS0 = E28*SQRT(TFLU/(4+0.85*TFLU))
  107. FI = FLU*ELAS0-1
  108. GO TO 100
  109. C
  110. C*****************************************************************
  111. C 3 : MODELE B3
  112. C*****************************************************************
  113. C*****************************************************************
  114. C
  115. C DATCUR = TEMPS DE CURE MINIMAL DU BETON
  116. C DATSEC = AGE DE DEBUT DE SECHAGE DU BETON POUR LE TERME DE DESSICCATION
  117. C
  118. C VOL EN M3, AC EN M2, RH EN % (EX : RH = 50)
  119. C
  120. C
  121. C
  122. C VALIDITE DU MODELE
  123. C*****************************************************************
  124. C
  125. C LA COMPOSITION DE BETON EST REDUITE A L UTILISATION DE CIMENT
  126. C PORTLAND AVEC :
  127. C
  128. C 17MPA<FC<70MPA
  129. C 0.3<W/C<0.85
  130. C 160kg/m3<CIMENT<720kg/m3
  131. C 2.5<A/C<13.5
  132. C
  133. C DATCUR = JOURS DE CURE OBLIGATOIRE
  134. 30 IF (DATCUR.LT.1.D0) THEN
  135. WRITE(*,*)'LA CURE MINI DU BETON EST DE 1 JOUR'
  136. ENDIF
  137. C
  138. C
  139. C COMMENTAIRES
  140. C*****************************************************************
  141. C
  142. C FLU= FONCTION DE FLUAGE OU DE COMPLAISANCE POUR TFLU ET DUR
  143. C FI = COEFFICIENT DE FLUAGE : FI = ET*FLU-1
  144. C ET = MODULE D ELASTICITE A TFLU
  145. C TFLU = DATE D'APPLICATION DE LA CHARGE
  146. C DUR = DUREE D APPLICATION DE LA CHARGE
  147. C
  148. C FLU0= TERME DE FLUAGE PROPRE POUR TFLU ET DUR
  149. C FLUD= TERME DE FLUAGE DE DESSICCATION POUR TFLU ET DUR
  150. C
  151. C E28 = 4734*SQRT(FC) = MODULE D ELASTICITE
  152. C A 28 JOURS DU BETON (MPA)
  153. C
  154. C EAU : EAU DE COMPOSITION DE BETON : VALEUR EN KG/M3
  155. C CIMENT : VALEUR EN KG/M3
  156. C AGREGA : VALEUR EN KG/M3
  157. C FC = RESISTANCE CARACTERISTIQUE A 28 JOURS DU BETON (MPA)
  158. C WSURC = W/C = RAPPORT EAU SUR CIMENT DE LA COMPO. DE BETON
  159. C
  160. C
  161. C
  162. C COMPOSITION DE CIMENT DE LA CENTRALE DE CIVAUX / BHP
  163. C*****************************************************************
  164. C VALEURS EN kg/m3
  165. C EAU = 161.
  166. C CIMENT = 266.
  167. C SILICE = 40.3
  168. C AGREG = 1133.0 + 482.0
  169. C SABLE = 1133.0
  170. C FINES = 57.0
  171. C FC = 30.
  172. C ADJUV = 9.98
  173. C
  174. C
  175. C
  176. C CIMENT PORTLAND
  177. C*****************************************************************
  178. C LA VALEUR DE CIMENT S EXPRIME EN kg/m3
  179. C = CIMENT
  180. IF (C.LT.160.D0.OR.C.GT.720.D0) THEN
  181. WRITE(*,*)'QUANTITE DE CIMENT NON RESPECTEE'
  182. ENDIF
  183. C
  184. C
  185. C RAPPORT AGREGAT/CIMENT
  186. C*****************************************************************
  187. C LA VALEUR D AGREGAT S EXPRIME EN kg/m3
  188. A = AGREG
  189. ASURC = A/C
  190. IF (ASURC.GT.2.50.OR.ASURC.LT.13.50) THEN
  191. WRITE(*,*)'RAPPORT A/C NON RESPECTE'
  192. ENDIF
  193. C
  194. C
  195. C RAPPORT EAU/CIMENT
  196. C*****************************************************************
  197. C LA VALEUR D EAU S EXPRIME EN kg/m3
  198. W = EAU
  199. WSURC = W/C
  200. C
  201. IF (WSURC.EQ.0.D0.OR.WSURC.GT.1.D0) THEN
  202. WSURC = 1/((FC/22.8)+0.535)
  203. ELSE
  204. IF (WSURC.LT.0.35.OR.WSURC.GT.0.85) THEN
  205. WRITE(*,*)'RAPPORT W/C NON RESPECTE'
  206. ENDIF
  207. ENDIF
  208. C
  209. C
  210. C RESISTANCE CARACTERISTIQUE DU BETON
  211. C*****************************************************************
  212. C LA VALEUR DE RESISTANCE S EXPRIME EN MPa
  213. C
  214. IF (FC.LT.17.D0.OR.FC.GT.70.D0) THEN
  215. WRITE(*,*)'VALEUR DE RESISTANCE NON RESPECTEE'
  216. ENDIF
  217. C
  218. C
  219. C*****************************************************************
  220. C TERME DE FLUAGE PROPRE FLU0
  221. C*****************************************************************
  222. C
  223. E28=4734*SQRT(FC)
  224. Q1 = 0.6*1E6/E28
  225. Q2 = 185.4*SQRT(C)*FC**(-0.9)
  226. Q3 = 0.29*((WSURC)**4.)*Q2
  227. Q4 = 20.3*(ASURC)**(-0.7)
  228. C
  229. C*****************************************************************
  230. C DANS LE MODELE B3
  231. C BETON NORMAUX : C1=m=0.5 ET C2=n=0.1
  232. C1 = 0.5
  233. C2 = 0.1
  234. C*****************************************************************
  235. C
  236. R = 1.7*(TFLU**(0.12))+8
  237. QF = 1/(0.086*(TFLU**(2/9))+1.21*(TFLU**(4/9)))
  238. C
  239. Z = (TFLU**(-C1))*LOG(1+(DUR)**(C2))
  240. Q = QF*(1+QF/(Z**R))**(-1/R)
  241. C
  242. FLU0 = Q2*Q+Q3*LOG(1+DUR**(C2))
  243. *+Q4*LOG((DUR+TFLU)/TFLU)
  244. C
  245. C
  246. C*****************************************************************
  247. C TERME DE FLUAGE DE DESSICCATION : FLUD
  248. C*****************************************************************
  249. C ALPHA1 = COEFFICIENT DU TYPE DE CIMENT
  250. C ALPHA1 = 1.00 = CIMENT TYPE 1
  251. C ALPHA1 = 0.85 = CIMENT TYPE 2
  252. C ALPHA1 = 1.10 = CIMENT TYPE 3
  253. C
  254. ALPHA1 = 1.00
  255. C
  256. C ALPHA2 = TYPE DE CURE
  257. C ALPHA2 = 0.75 = EVAPORATION
  258. C ALPHA2 = 1.00 = CURE DANS L EAU OU A 100% HR
  259. C ALPHA2 = 1.20 = CURE NORMALE OU CONFINE AVEC PROTECTION
  260. C CONTRE LE SECHAGE
  261. C
  262. ALPHA2 = 1.00
  263. C
  264. C
  265. C TERME DE RETRAIT INFINI : ESINF
  266. C*****************************************************************
  267. ESINF = (-ALPHA1*ALPHA2*(1.9*1E-2*(W**(2.1))*
  268. *(FC**(-0.28))+270))*1E-6
  269. C
  270. COEFT = 8.5*(DATSEC**(-0.08))*(FC**(-1/4))
  271. C
  272. C
  273. C TERME DE FORME : COEFS
  274. C*****************************************************************
  275. C COEFS = 1.00 : PLAQUE INFINIE
  276. C COEFS = 1.15 : CYLINDRE INFINI
  277. C COEFS = 1.25 : PRISME
  278. C COEFS = 1.30 : SPHERE
  279. C COEFS = 1.35 : CUBE
  280. C
  281. COEFS = 1.00
  282. C
  283. C
  284. C TERME DE TEMPS DU A LA TAILLE : TAUSH
  285. C*****************************************************************
  286. C LES DIMENSIONS SONT DONNEES EN M3 & M2 : TRANSFORMATION EN MM
  287. C
  288. DIM = 2*(VOL*1E9)/(AC*1E6)
  289. TAUSH = COEFT*((COEFS*DIM)**2)
  290. C
  291. C
  292. C TERME DEPENDANT DU RETRAIT ULTIME : ESHINF
  293. C*****************************************************************
  294. E607 = E28*SQRT(607/(4+0.85*607))
  295. C
  296. ESEC = E28*SQRT((DATSEC+TAUSH)/
  297. *(4+0.85*(DATSEC+TAUSH)))
  298. C
  299. ESHINF = ESINF*E607/ESEC
  300. C
  301. C
  302. C TERME DEPENDANT DE LA DUREE DE SECHAGE : ST
  303. C*****************************************************************
  304. Q5 = 7.57*1E5/FC*((ABS(ESINF))**(-0.6))
  305. ST = TANH(SQRT((TFLU+DUR-DATSEC)/TAUSH))
  306. C
  307. IF (TFLU.GE.DATSEC) THEN
  308. STTFLU = TANH(SQRT((TFLU-DATSEC)/TAUSH))
  309. ELSE
  310. STTFLU = TANH(SQRT((DATSEC-DATSEC)/TAUSH))
  311. ENDIF
  312. C
  313. HTDUR = 1-(1-RH/100)*ST
  314. HTTFLU = 1-(1-RH/100)*STTFLU
  315. FLUD=Q5*SQRT(EXP(-8*HTDUR)-EXP(-8*HTTFLU))
  316. C
  317. C
  318. C COEFFICIENT DE FLUAGE FI
  319. C*****************************************************************
  320. ELAS0 = E28*SQRT(TFLU/(4+0.85*TFLU))
  321. FI = ELAS0*FLU-1
  322. C
  323. C
  324. C TYPE DE FLUAGE CONSIDERE
  325. C*****************************************************************
  326. C ICAS = 1 : CALCUL DU FLUAGE PROPRE SEULEMENT
  327. C ICAS = 2 : CALCUL DU FLUAGE PROPRE + FLUAGE DE DESSICCATION
  328. C
  329. ICAS = 1
  330. IF (ICAS.EQ.1) THEN
  331. FLU = FLU0
  332. ELSE IF (ICAS.EQ.2) THEN
  333. FLU = FLU0 + FLUD
  334. ENDIF
  335. C
  336. C*****************************************************************
  337. C 4 : MODELE CEB FIP
  338. C*****************************************************************
  339. C*****************************************************************
  340. C
  341. C ICODEA = TYPE D AGGREGAT
  342. C ICODEA = 1 = AGGREGAT QUARTZITIQUE
  343. C ICODEA = 2 = AUTRE AGGREGAT
  344. C
  345. 40 ICODEA = 1
  346. C
  347. C ANALYSE ELASTIQUE DES GRANULATS
  348. C IANA = 1 = OUI (SANS FLUAGE)
  349. C IANA = 0 = NON (VISCOELASTIQUE)
  350. C
  351. IANA = 0
  352. C
  353. C MODULE D ELASTICITE
  354. C*****************************************************************
  355. C CH0 = HAUTEUR CARACTERISTIQUE
  356. CH0 = 2*AC*1E6/(U*1E3)
  357. FCM = FC+8
  358. FCM0 = 10.D0
  359. E28 = 21500*((FCM/FCM0)**(1/3))
  360. C
  361. IF (ICODEA.EQ.2) THEN
  362. C
  363. IF (IANA.EQ.1) THEN
  364. E28 = E28*0.85
  365. ELSE
  366. CONTINUE
  367. ENDIF
  368. C
  369. E28MAX = E28*1.20
  370. E28MIN = E28*0.7
  371. E28MOY = (E28MAX+E28MIN)/2
  372. E28 = E28MOY
  373. C
  374. ELSE
  375. IF (IANA.EQ.1) THEN
  376. E28 = E28*0.85
  377. ELSE
  378. CONTINUE
  379. ENDIF
  380. CONTINUE
  381. ENDIF
  382. C
  383. C
  384. C RESISTANCE ET MODULE
  385. C*****************************************************************
  386. C CEMENT = COEFFICIENT FONCTION DU TYPE DE CIMENT
  387. C CIMALPHA = EFFET DU TYPE DE CIMENT SUR LE FLUAGE
  388. C
  389. C CIMENT A DURCISSEMENT LENT
  390. C*****************************************************************
  391. C CEMENT = 0.38
  392. C CIMALPHA = -1.D0
  393. C
  394. C CIMENT NORMAL ET DURCISSEMENT RAPIDE
  395. C*****************************************************************
  396. CEMENT = 0.25
  397. CIMALPHA = 0.D0
  398. C
  399. C CIMENT HAUTE RESISTANCE ET DURCISSEMENT RAPIDE
  400. C*****************************************************************
  401. C CEMENT = 0.20
  402. C CIMALPHA = 1.D0
  403. C
  404. C
  405. IF (ITYPE.EQ.4.OR.ITYPE.EQ.5) THEN
  406. TFLU = TFLU*(9/(2+TFLU**(1.2))+1)**(CIMALPHA)
  407. ENDIF
  408. C
  409. C
  410. C MODULE AU TEMPS T = TFLU+DUR ET RESISTANCE
  411. C*****************************************************************
  412. BCC = EXP(CEMENT*(1-SQRT(28/(TFLU+DUR))))
  413. BE = SQRT(BCC)
  414. ELAS = BE*E28
  415. FCMT = FCM*BCC
  416. C
  417. C
  418. C MODULE AU TEMPS TFLU ET RESISTANCE
  419. C*****************************************************************
  420. BCC0 = EXP(CEMENT*(1-SQRT(28/TFLU)))
  421. BE0 = SQRT(BCC0)
  422. ELAS0 = BE0*E28
  423. FCMT0 = FCM*BCC0
  424. C
  425. C
  426. C COEFFICIENT DE FLUAGE : FI
  427. C*****************************************************************
  428. BFCM = 5.3/(SQRT(FCMT0/FCM0))
  429. BT0 = 1/(0.1+TFLU**0.2)
  430. C
  431. REFH = 100.D0
  432. BH = 150*CH0/REFH*(1+(0.012*RH)**(18))+250
  433. C
  434. IF (BH.GT.1500) THEN
  435. BH = 1500.D0
  436. ENDIF
  437. C
  438. BCTT0 = (DUR/(DUR+BH))**(0.3)
  439. PHIRH = 1+(1-(RH/100))/(0.46*((CH0/REFH)**(1/3)))
  440. PHI0 = PHIRH*BFCM*BT0
  441. C
  442. IF (ITYPE.EQ.4) THEN
  443. FI = PHI0*BCTT0*BE0
  444. FICEB = FI/BE0
  445.  
  446. FLU = 1/ELAS0*(1+FI)
  447. GO TO 100
  448. C
  449. C*****************************************************************
  450. C 5 : MODELE EUROCODE2
  451. C*****************************************************************
  452. C*****************************************************************
  453. C
  454. C
  455. C VALIDITE DU MODELE
  456. C*****************************************************************
  457. C
  458. C LA COMPOSITION DE BETON EST REDUITE A L UTILISATION DE CIMENT
  459. C PORTLAND AVEC :
  460. C
  461. C 12 MPA < FC < 60 MPA
  462. C 40% < RH < 100%
  463. C
  464. ELSE IF (ITYPE.EQ.5) THEN
  465. C
  466. IF (FC.LT.12.OR.FC.GT.60) THEN
  467. WRITE(*,*)'MODELE NON VALIDE'
  468. ELSE IF (RH.LT.12.OR.RH.GT.60) THEN
  469. WRITE(*,*)'MODELE NON VALIDE'
  470. ENDIF
  471. C
  472. E28 = 1.05*9500*(FCM)**(1/3)
  473. FCMT0 = FCM*EXP(CEMENT*(1-SQRT(28/TFLU)))
  474. ELAS0 = 1.05*9500*(FCMT0)**(1/3)
  475. PHIRH = 1+(10*(1-RH/100))/(CH0**(1/3))
  476. BFCM = 16.8/SQRT(FCM)
  477. PHITT0 = PHIRH*BFCM*BT0*BCTT0
  478. FI = PHITT0*ELAS0/E28
  479. FLU = 1/ELAS0*(1+FI)
  480. GO TO 100
  481. C
  482. C*****************************************************************
  483. C 6 : MODELE BPEL91
  484. C*****************************************************************
  485. C*****************************************************************
  486. C
  487. ELSE IF (ITYPE.EQ.6) THEN
  488. C
  489. IF (FC.LE.40) THEN
  490. FCJ = TFLU/(4.76+0.83*TFLU)*FC
  491. ELSE
  492. FCJ = TFLU/(1.40+0.95*TFLU)*FC
  493. ENDIF
  494. C
  495. FCLIMI = 1.1*FC
  496. IF (FCJ.GT.FCLIMI) THEN
  497. FCJ = FCLIMI
  498. ENDIF
  499. C
  500. ELAS0 = 11000*(FCJ)**(1/3)
  501. E28 = 11000*(FC)**(1/3)
  502. C
  503. RM = AC*1E6/(U*1E3)
  504. U1 = 5.
  505. FTT0 = SQRT(DUR)/(SQRT(DUR)+U1*SQRT(RM))
  506. CASC = (120.D0-RH)/30+2*(100.D0-RH)/(3*(20+RM))
  507. CAST0 = 100.D0/(100.D0+TFLU)
  508. CASFL = 0.4D0+CASC*CAST0
  509. FI = CASFL*FTT0
  510. FLU = 1/ELAS0+FI/E28
  511. C
  512. GO TO 100
  513. C
  514. ENDIF
  515. C
  516. C*****************************************************************
  517. C 7 : MODELE AFREM
  518. C*****************************************************************
  519. C
  520. C FLUAGE PROPRE
  521. C*****************************************************************
  522. C ICOMP = 1 SI BETON AVEC FUMEE DE SILICE
  523. C ICOMP = 0 SI BETON SANS FUMEE DE SILICE
  524. C
  525. 50 ICOMP = 1
  526. C
  527. IF (FC.LE.40) THEN
  528. FCJ = TFLU/(4.76+0.83*TFLU)*FC
  529. ELSE
  530. FCJ = TFLU/(1.40+0.95*TFLU)*FC
  531. ENDIF
  532. C
  533. FCLIMI = 1.1*FC
  534. IF (FCJ.GT.FCLIMI) THEN
  535. FCJ = FCLIMI
  536. ENDIF
  537. C
  538. ELAS0 = 11000*(FCJ)**(1/3)
  539. E28 = 11000*(FC)**(1/3)
  540.  
  541. IF (ICOMP.EQ.1) THEN
  542. FLU0 = 3.60/E28*SQRT(DUR)/(SQRT(DUR)
  543. *+0.37*EXP(2.8*FCJ/FC))
  544. C
  545. ELSE IF (ICOMP.EQ.0) THEN
  546. FLU0 = 1.40/E28*SQRT(DUR)/(SQRT(DUR)
  547. *+0.40*EXP(3.1*FCJ/FC))
  548. ENDIF
  549. C
  550. C RETRAIT
  551. C*****************************************************************
  552. C
  553. IF (FC.GE.57) THEN
  554. CASFC28 = 30-0.21*FC
  555. ELSE
  556. CASFC28 = 18.
  557. ENDIF
  558. C
  559. IF (ICOMP.EQ.1) THEN
  560. ERET1 = CASFC28*(72*EXP(-0.046*FC)+75-RH/100)
  561. **DATSEC/(DATSEC+2.8*(RM**2))
  562. ERET = CASFC28*(72*EXP(-0.046*FC)+75-RH/100)
  563. **DUR/(DUR+2.8*(RM**2))
  564. GO TO 51
  565. C
  566. ELSE IF (ICOMP.EQ.0) THEN
  567. C
  568. ERET1 = CASFC28*(72*EXP(-0.046*FC)+75-RH/100)
  569. **DATSEC/(DATSEC+8.4*(RM**2))
  570. ERET = CASFC28*(72*EXP(-0.046*FC)+75-RH/100)
  571. **DUR/(DUR+8.4*(RM**2))
  572. GO TO 52
  573. C
  574. ENDIF
  575. C
  576. C FLUAGE DE DESSICCATION
  577. C*****************************************************************
  578. 51 IF (ICOMP.EQ.1) THEN
  579. FLUD = 3.2*1E3/E28*(ERET-ERET1)
  580. ENDIF
  581. C
  582. 52 IF (ICOMP.EQ.0) THEN
  583. FLUD = 1E3/E28*(ERET-ERET1)
  584. ENDIF
  585. C
  586. C
  587. C TYPE DE FLUAGE CONSIDERE
  588. C*****************************************************************
  589. C ICAS = 1 : CALCUL DU FLUAGE PROPRE SEULEMENT
  590. C ICAS = 2 : CALCUL DU FLUAGE PROPRE + FLUAGE DE DESSICCATION
  591. C
  592. ICAS = 1
  593. IF (ICAS.EQ.1) THEN
  594. FLU = FLU0
  595. ELSE IF (ICAS.EQ.2) THEN
  596. FLU = FLU0 + FLUD
  597. ENDIF
  598. C
  599. GO TO 100
  600. C
  601. C*****************************************************************
  602. C 8 : MODELE BP-KX
  603. C*****************************************************************
  604. C*****************************************************************
  605. C
  606. 60 E28=4734*SQRT(FC)
  607. Q1 = 4597/E28
  608. WSURC = EAU/CIMENT
  609. Q2 = 0.00045*WSURC**(0.8)*C**(1.5)
  610. **(1-A/2500)**(-0.9)*FC**(-0.5)-0.39
  611. Q3 = 0.025*Q2
  612. Q4 = 0.017*WSURC**(2.3)*C**(0.2)
  613. **(1-A/2500)**(0.39)*FC**(0.46)
  614. Q5 = 3.321*FC**(-0.5)
  615.  
  616. R = 1.7*(TFLU**(0.12))+8
  617. ZTT0 = TFLU**(-0.5)*LOG(1+DUR**(0.1))
  618. QFT0 = 1/(0.086*(TFLU**(2/9))+1.21
  619. **(TFLU**(4/9)))
  620. QTT0 = QFT0*((1+((QFT0/ZTT0)**(R)))**(-1/R))
  621. FLU0 = Q1+Q2*QTT0+Q3*LOG(1+DUR)
  622. *+Q4*LOG((TFLU+DUR)/TFLU)
  623. C
  624. C
  625. C TERME D HUMIDITE : COEFH
  626. C*****************************************************************
  627. HUMI = RH/100
  628. IF(HUMI.EQ.1.D0) THEN
  629. COEFH = -0.2
  630. ELSE IF(HUMI.LE.0.98) THEN
  631. COEFH = 1-(HUMI)**3
  632. ELSE
  633. COEFH = (1-0.98**3)-((1-0.98**3)
  634. *+0.2)/0.02*(HUMI-0.98)
  635. ENDIF
  636. C
  637. RAP = AGREG/SABLE
  638. IF (RAP.GT.2.6) THEN
  639. ALPHA4 = 0.7+0.3*(RAP)**(-3)
  640. ELSE
  641. ALPHA4 = 1.
  642. ENDIF
  643.  
  644. ALPHA1 = 0.0175*WSURC**(1.5)*C**(1.1)
  645. **FC**(-0.2)*(1-A/2500)*ALPHA4
  646. C
  647. C ALPHA2 = COEFFICIENT DU TYPE DE CIMENT
  648. C ALPHA2 = 1.00 = CIMENT TYPE 1
  649. C ALPHA2 = 0.85 = CIMENT TYPE 2
  650. C ALPHA2 = 1.15 = CIMENT TYPE 3
  651. C
  652. ALPHA2 = 1.00
  653. C
  654. C ALPHA3 = TYPE DE CURE
  655. C ALPHA3 = 0.74 = EVAPORATION OU CURE SECHE
  656. C ALPHA3 = 1.00 = CURE DANS L EAU OU A 100% HR
  657. C ALPHA3 = 1.40 = CURE CONFINE AVEC PROTECTION
  658. C CONTRE LE SECHAGE
  659. C
  660. ALPHA3 = 1.00
  661. C
  662. ESINF = (1.15*ALPHA1+0.16)*ALPHA2*ALPHA3
  663. C
  664. DIM = 2*(VOL*1E9)/(AC*1E6)
  665. C
  666. C
  667. C TERME DE FORME : COEFS
  668. C*****************************************************************
  669. C COEFS = 1.00 : PLAQUE INFINIE
  670. C COEFS = 1.15 : CYLINDRE INFINI
  671. C COEFS = 1.25 : PRISME
  672. C COEFS = 1.30 : SPHERE
  673. C COEFS = 1.35 : CUBE
  674. C
  675. COEFS = 1.00
  676. C
  677. IF (ALPHA3.EQ.0.74) THEN
  678. TCUR = 25+TFLU
  679. ELSE
  680. TCUR = TFLU
  681. ENDIF
  682. C
  683. C1TCUR = 10*(0.6+SQRT(4.5/TCUR))
  684. IF (C1TCUR.GT.18.D0) THEN
  685. C1TCUR = 18.D0
  686. ENDIF
  687. C
  688. TAUSH = 0.32*((COEFS*DIM)**2)/C1TCUR
  689. X1 = 17+TAUSH/40
  690. X2 = 12+TCUR
  691. GX1 = SQRT(X1/(4+0.9*X1))
  692. GX2 = SQRT(X2/(4+0.9*X2))
  693. ESHINF = ESINF/(GX1*GX2)
  694. C
  695. FLUD = Q5*COEFH*ESHINF*(TANH(SQRT((TFLU+DUR-DATSEC)
  696. **(2*TAUSH)))-TANH(SQRT((TFLU-DATSEC)*(2*TAUSH))))**(0.5)
  697. C
  698. C
  699. C TYPE DE FLUAGE CONSIDERE
  700. C*****************************************************************
  701. C ICAS = 1 : CALCUL DU FLUAGE PROPRE SEULEMENT
  702. C ICAS = 2 : CALCUL DU FLUAGE PROPRE + FLUAGE DE DESSICCATION
  703. C
  704. ICAS = 1
  705. IF (ICAS.EQ.1) THEN
  706. FLU = 145*FLU0*1E-6
  707. ELSE IF (ICAS.EQ.2) THEN
  708. FLU = 145*(FLU0 + FLUD)*1E-6
  709. ENDIF
  710. C
  711. GO TO 100
  712. C
  713. C*****************************************************************
  714. C
  715. 100 RETURN
  716. END
  717.  
  718.  
  719.  
  720.  
  721.  
  722.  
  723.  

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