Télécharger pfluage.procedur

Retour à la liste

Numérotation des lignes :

  1. * PFLUAGE PROCEDUR CB215821 16/06/17 21:15:02 8979
  2. DEBP PFLUAGE NOMREG*MOT MOBET*MMODEL MABET*MCHAML
  3. SIGEPS*MCHAML RAYSECH*MCHAML AGECHAB*FLOTTANT
  4. TEMPCOUR*FLOTTANT RHOS*MCHAML TCOULB*FLOTTANT
  5. FLOT1/FLOTTANT FC28/FLOTTANT ;
  6. *
  7. * RAYSECH : Rayon de sechage en cm
  8. * AGECHAB : instant de la mise en chargement(meme base que les autres)
  9. * TEMPCOUR : Instant courant en jours
  10. * RHOS : Taux d'armatures
  11. * TCOULA : Instants de coulage du beton en jours
  12. * (Defaut : 0.)
  13. * FLOT1 : Humidite du milieu environnant (formulation BPEL)
  14. * FC28 : Resistance caracteristique en traction a 28 jours en Pa
  15. * (Defaut : 40E6)
  16. *
  17. CHAMUNIT = MANU CHML MOBET 'SCAL' 1. TYPE 'CONTRAINTES' STRESSES ;
  18. RAYSECH = EXCO EPAI RAYSECH SCAL ;
  19. AGECHAR = (AGECHAB - tcoulb) *CHAMUNIT ;
  20. TCOULA = TCOULB*CHAMUNIT ;
  21. RHOS = RHOS 'NOMC' 'SCAL' ;
  22. *
  23. * MODIFICATIONS B. CIREE DU 25/02/03
  24. * AJOUT DES FORMULATIONS BPEL MISE A JOUR 99
  25. * CORRECTION D'UNE CONFUSION ENTRE L'AGE DU BETON ET
  26. * L'INSTANT COURANT POUR LES FORMULATIONS BPEL
  27. *
  28. SI (EGA NOMREG 'BPEL99') ;
  29. *
  30. MESS 'Calcul du Fluage avec la formulation BPEL99' ;
  31. MESS '-------------------------------------------' ;
  32. *
  33. SI (NON (EXIS FLOT1)) ;
  34. MESS '--- Pas de donnees pour l humidite du milieu environnant' ;
  35. MESS ' pour la formulation BPEL99 - On quitte la procedure' ;
  36. QUIT PFLUAGE ;
  37. FINSI ;
  38. SI (NON (EXIS FC28)) ;
  39. MESS '--- Pas de Resistance caracteristique du beton a 28 jours' ;
  40. MESS ' On prend 40.E6 Pa par defaut ' ;
  41. FC28 = 40.E6 ;
  42. FINSI ;
  43. *
  44. MESS 'Instant de calcul du fluage :' TEMPCOUR 'jours' ;
  45. MESS 'Instant de coulage du beton :'
  46. (MINI TCOULA) (MAXI TCOULA) 'jours' ;
  47.  
  48. MESS 'Age du beton au chargement :'
  49. (MINI AGECHAR) (MAXI AGECHAR) 'jours' ;
  50. *
  51. * --- FLUAGE PROPRE
  52. *
  53. HUMREL = FLOT1 ;
  54. EI28 = 11000.E6 * ((FC28 /1.E6) ** (1. / 3.)) ;
  55. *
  56. AGEBETON = (TEMPCOUR * CHAMUNIT) - TCOULA ;
  57. *
  58. A1 = AGEBETON * (((4.76 * CHAMUNIT) + (0.83 * AGEBETON)) ** -1) ;
  59. FCJ = ((A1 MASQ 'EGSUPE' 1.1) * 1.1 * FC28)
  60. + ((A1 MASQ 'INFERIEUR' 1.1) * A1 * FC28) ;
  61. *
  62. TD2 = (AGEBETON - AGECHAR) ** 0.5 ;
  63. A1 = 0.40 * (EXP (3.1 * (FCJ / FC28))) ;
  64. EPSFP = 1.4 * (EI28 ** -1)
  65. * TD2 * ((TD2 + (A1 * CHAMUNIT)) ** -1) ;
  66. *
  67. *
  68. * --- FLUAGE DE DESSICCATION
  69. *
  70. KFC28 = 18. ;
  71. A1 = KFC28 * ((72. * (EXP (-0.046E-6 * FC28))) + (75. - HUMREL))
  72. * 1.E-6 ;
  73.  
  74. * CB215821 : Formulation pour eviter la division par zero...
  75. * A2 = 1. + ((8.4 * (RAYSECH ** 2)) * (AGEBETON ** -1) ) ;
  76. * EPSRD1 = A1 * (A2 ** -1.) ;
  77. A2 = AGEBETON + (8.4 * (RAYSECH ** 2)) ;
  78. EPSRD1 = A1 * (AGEBETON * (A2 ** -1)) ;
  79.  
  80.  
  81.  
  82. * CB215821 : Formulation pour eviter la division par zero...
  83. * A2 = 1. + ((8.4 * (RAYSECH ** 2)) * (AGECHAR ** -1.) ) ;
  84. * EPSRD0 = A1 * (A2 ** -1.) ;
  85. A2 = AGECHAR + (8.4 * (RAYSECH ** 2)) ;
  86. EPSRD0 = A1 * (AGECHAR * (A2 ** -1)) ;
  87. *
  88. EPSFD = 3.2E3 * (EI28 ** -1) * (EPSRD1 - EPSRD0) ;
  89. *
  90. *
  91. * --- FLUAGE TOTAL
  92. *
  93. MCHAMC = (EPSFP + EPSFD) * ((1. + (15. * RHOS)) ** -1) ;
  94. *
  95. MESS 'Module elastique du beton a 28 jours en Pa :' EI28 ;
  96. MESS 'Resistance caracteristique a 28 jours en Pa :' FC28 ;
  97. MESS 'Resistance caracteristique a' TEMPCOUR ' jours :'
  98. (MINI FCJ) (MAXI FCJ) ;
  99. MESS 'EPSRD1 ' (MINI EPSRD1) (MAXI EPSRD1) ;
  100. MESS 'EPSRD0 ' (MINI EPSRD0) (MAXI EPSRD0) ;
  101. MESS 'Deformation de Fluage de dessiccation unitaire : '
  102. (MINI EPSFD) (MAXI EPSFD) ;
  103. MESS 'Deformation de Fluage propre unitaire : '
  104. (MINI EPSFP) (MAXI EPSFP) ;
  105. MESS 'Deformation de Fluage total unitaire : '
  106. (MINI MCHAMC) (MAXI MCHAMC) ;
  107. FINSI ;
  108. *
  109. *
  110. * La procedure d'utilisation du BPEL pour le fluage a ete modifiee
  111. * lorsqu'il y a une division de deux MCHAML. Il faut faire la meme
  112. * chose pour Eurocode 2 et LCPC-GRANGER. Chose non faite par
  113. * A. KANJI NANJI de SOCOTEC INDUSTRIE.
  114. *
  115. SI (EGA NOMREG 'BPEL91') ;
  116. *
  117. MESS 'Calcul du Fluage avec la formulation BPEL91' ;
  118. MESS '-------------------------------------------' ;
  119. *
  120. SI (NON (EXIS TCOULA)) ;
  121. TCOULA = ZERO MOBET 'STRESSES' ;
  122. MESS '--- Pas de donnees d Instants de coulage du beton ' ;
  123. MESS ' On prend ZERO par defaut ' ;
  124. FINSI ;
  125. SI (NON (EXIS FLOT1)) ;
  126. MESS '--- Pas de donnees pour l humidite du milieu environnant' ;
  127. MESS ' pour la formulation BPEL91 - On quitte la procedure' ;
  128. QUIT PFLUAGE ;
  129. FINSI ;
  130. *
  131. MESS 'Instant de calcul du fluage :' TEMPCOUR 'jours' ;
  132. MESS 'Instant de coulage du beton :'
  133. (MINI TCOULA) (MAXI TCOULA) 'jours' ;
  134. MESS 'Age du beton au chargement :'
  135. (MINI AGECHAR) (MAXI AGECHAR) 'jours' ;
  136. *
  137. HUMREL = FLOT1 ;
  138. EI28 = 11000.E6 * ((FC28 /1.E6) ** (1. / 3.)) ;
  139. *
  140. MESS 'Module elastique du beton a 28 jours en Pa :' EI28 ;
  141. *
  142. * RHOS ETANT UN CHAMPS PAR ELEMENT ATTENTION AU CALCUL DE KS
  143. *
  144. *
  145. KS = (1. + (20.*RHOS)) * -1 ;
  146. * MESS 'FLUAGE MINI DE KS' (MINI KS) 'MAXI KS' (MAXI KS) ;
  147. * VBID = RHOS**-1. ;
  148. * VBID1 = RHOS1 * VBID ;
  149. * KS0 = VBID1 + (20. * RHOS) ;
  150. * KS = KS0**-1. ;
  151. *
  152. * ------- Gtv
  153. *
  154. KE = 0.4*CHAMUNIT ;
  155. * MESS 'FLUAGE MINI DE KE' (MINI KE) 'MAXI KE' (MAXI KE) ;
  156. *
  157. * ------- KC = Gech(rm)
  158. *
  159. KC = ((120. - HUMREL)/30.) +
  160. (((2./3.)*(100. - HUMREL))*((20. + RAYSECH) ** -1)) ;
  161. * MESS 'FLUAGE MINI DE KC' (MINI KC) 'MAXI KC' (MAXI KC) ;
  162. *
  163. * ------- Gac(w)
  164. *
  165. K_T1 = 100.*((100. + AGECHAR)**-1) ;
  166. * MESS 'FLUAGE MINI DE KT' (MINI K_T1) 'MAXI KT' (MAXI K_T1) ;
  167. *
  168. * ------- Ffluage(t-taux) avec taux = AGECHAR
  169. *
  170. AGEBETON = (TEMPCOUR * CHAMUNIT) - TCOULA ;
  171. NUMER = (AGEBETON - AGECHAR)**0.5 ;
  172. CONST = RAYSECH**0.5 ;
  173. DENOM = NUMER + (5.*CONST) ; DENOM = DENOM**(-1) ;
  174. F_DT = NUMER * DENOM ;
  175. * MESS 'FLUAGE MINI DE ft' (MINI F_DT) 'MAXI FT' (MAXI F_DT) ;
  176. *
  177. * -- KS*(KE + (KC*K_T1)) = KS*(Gtv + Gech(rm)*Gac(w) = Kfluage(w)
  178. * -- E_BETON = E_BETON28
  179. *
  180. * -- MCHAMC = (1./E_BETON28)* Kfluage(taux)*Ffluage(t-taux)
  181. *
  182. MCHAMC = (EI28 ** -1)*KS*(KE + (KC * K_T1)) * F_DT ;
  183. MESS 'Deformation de Fluage total unitaire : '
  184. (MINI MCHAMC) (MAXI MCHAMC) ;
  185. FINSI ;
  186. *
  187. SI (EGA NOMREG 'EC2') ;
  188. E_BETON = CHAN STRESSES MOBET
  189. (CHAN TYPE (EXCO YOUN MABET SCAL) 'CONTRAINTES') ;
  190. HUMREL = FLOT1 ;
  191. FCM = FLOT2 ;
  192. *
  193. PHIHR = CHAMUNIT + ((1. - (HUMREL/100.)) /
  194. (0.46*(RAYSECH/100.)**(1./3.))) ;
  195. BETAFCM = 5.3/((FCM/10.)**0.5) ;
  196. BETATAU = 1./((0.1*CHAMUNIT) + (AGECHAR**0.2)) ;
  197. BETAH = 1500.*CHAMUNIT ;
  198. BETACIN = (((TEMPCOUR*CHAMUNIT) - AGECHAR) /
  199. (((TEMPCOUR*CHAMUNIT) - AGECHAR) + BETAH))**0.3 ;
  200. *
  201. MCHAMC = (E_BETON** -1 )*PHIHR*BETAFCM*BETATAU*BETACIN ;
  202. FINSI ;
  203. *
  204. SI (EGA NOMREG 'LG') ;
  205. *
  206. SI (NON (EGA (TYPE FLOT1) 'FLOTTANT')) ;
  207. KSA = 1./(1. + (16.*1.071E-2)) ;
  208. KSIGP = 1./(1. + (20.*3.790E-3)) ;
  209. SINON ;
  210. KSA = FLOT1 ;
  211. KSIGP = FLOT2 ;
  212. FINSI ;
  213. *
  214. K2 = (2.*160./(3.*RAYSECH)) ** 2 ;
  215. *
  216. RAPVTAU = (0.1 + (28.**0.2))/((0.1*CHAMUNIT) + (AGECHAR**0.2)) ;
  217. *
  218. EEPREUVE = 31500.E6 ;
  219. EPSFP = 65.63E-6 ;
  220. AFP = 0.4969242 ;
  221. BFPJ = 21.97750113 ;
  222. EPSFT = 262.6E-6 ;
  223. AFT = 0.30071 ;
  224. BFTJ = 17.39335116 ;
  225. *
  226. FPTMTO = EPSFP*EEPREUVE*
  227. ((TEMPCOUR - AGECHAR)**AFP)/
  228. (((TEMPCOUR - AGECHAR)**AFP) + BFPJ) ;
  229. *
  230. FTK2TMTO = EPSFT*EEPREUVE*
  231. ((K2*(TEMPCOUR - AGECHAR))**AFT)/
  232. (((K2*(TEMPCOUR - AGECHAR))**AFT) + BFTJ) ;
  233. *
  234. FPK2TMTO = EPSFP*EEPREUVE*
  235. ((K2*(TEMPCOUR - AGECHAR))**AFP)/
  236. (((K2*(TEMPCOUR - AGECHAR))**AFP) + BFPJ) ;
  237. *
  238. MCHAMC = 1.E-6 * (E_BETON** -1 )*KSA*KSIGP*RAPVTAU*
  239. (FPTMTO + FTK2TMTO - FPK2TMTO) ;
  240. FINSI ;
  241. *
  242. *
  243. * --- Construction du tenseur de deformation differee de fluage
  244. *
  245. * | A D E |
  246. * EPSFLU = | D B F |
  247. * | E F C |
  248. *
  249. EPSFLUXX = EXCO SCAL ((EXCO SMXX SIGEPS 'SCAL')*MCHAMC) EPXX ;
  250. EPSFLUYY = EXCO SCAL ((EXCO SMYY SIGEPS 'SCAL')*MCHAMC) EPYY ;
  251. EPSFLUZZ = EXCO SCAL ((EXCO SMZZ SIGEPS 'SCAL')*MCHAMC) EPZZ ;
  252. EPSFLUXY = EXCO SCAL (2.*(EXCO SMXY SIGEPS 'SCAL')*MCHAMC) GAXY ;
  253. EPSFLUXZ = EXCO SCAL (2.*(EXCO SMXZ SIGEPS 'SCAL')*MCHAMC) GAXZ ;
  254. EPSFLUYZ = EXCO SCAL (2.*(EXCO SMYZ SIGEPS 'SCAL')*MCHAMC) GAYZ ;
  255. *
  256. EPSFLU = EPSFLUXX ET EPSFLUYY ET EPSFLUZZ ET
  257. EPSFLUXY ET EPSFLUXZ ET EPSFLUYZ ;
  258. EPSFLU = CHAN 'TYPE' EPSFLU 'DEFORMATIONS' ;
  259. *
  260. MESS 'Deformation de Fluage (min/max) : '
  261. (MINI (EPSFLUXX ET EPSFLUYY ET EPSFLUZZ ET
  262. ((EPSFLUXY ET EPSFLUXZ ET EPSFLUYZ) / 2.)))
  263. (MAXI (EPSFLUXX ET EPSFLUYY ET EPSFLUZZ ET
  264. ((EPSFLUXY ET EPSFLUXZ ET EPSFLUYZ) / 2.))) ;
  265. *
  266. FINP EPSFLU ;
  267.  
  268.  
  269.  
  270.  

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