1 : * PFLUAGE PROCEDUR CHAT 08/04/29 21:15:00 6098 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 : A2 = CHAMUNIT + ((8.4 * RAYSECH * RAYSECH) * (AGEBETON ** -1.) ) ; 74 : EPSRD1 = A1 * (A2 ** -1.) ; 75 : * 76 : A2 = CHAMUNIT + ((8.4 * RAYSECH * RAYSECH) * (AGECHAR ** -1.) ) ; 77 : EPSRD0 = A1 * (A2 ** -1.) ; 78 : * 79 : EPSFD = 3.2E3 * (EI28 ** -1.) * (EPSRD1 - EPSRD0) ; 80 : * 81 : * 82 : * --- FLUAGE TOTAL 83 : * 84 : MCHAMC = (EPSFP + EPSFD) * ((CHAMUNIT + (15. * RHOS)) ** -1.) ; 85 : * 86 : MESS 'Module elastique du beton a 28 jours en Pa :' EI28 ; 87 : MESS 'Resistance caracteristique a 28 jours en Pa :' FC28 ; 88 : MESS 'Resistance caracteristique a' TEMPCOUR ' jours :' 89 : (MINI FCJ) (MAXI FCJ) ; 90 : MESS 'EPSRD1 ' (MINI EPSRD1) (MAXI EPSRD1) ; 91 : MESS 'EPSRD0 ' (MINI EPSRD0) (MAXI EPSRD0) ; 92 : MESS 'Deformation de Fluage de dessiccation unitaire : ' 93 : (MINI EPSFD) (MAXI EPSFD) ; 94 : MESS 'Deformation de Fluage propre unitaire : ' 95 : (MINI EPSFP) (MAXI EPSFP) ; 96 : MESS 'Deformation de Fluage total unitaire : ' 97 : (MINI MCHAMC) (MAXI MCHAMC) ; 98 : FINSI ; 99 : * 100 : * 101 : * La procedure d'utilisation du BPEL pour le fluage a ete modifiee 102 : * lorsqu'il y a une division de deux MCHAML. Il faut faire la meme 103 : * chose pour Eurocode 2 et LCPC-GRANGER. Chose non faite par 104 : * A. KANJI NANJI de SOCOTEC INDUSTRIE. 105 : * 106 : SI (EGA NOMREG 'BPEL91') ; 107 : * 108 : MESS 'Calcul du Fluage avec la formulation BPEL91' ; 109 : MESS '-------------------------------------------' ; 110 : * 111 : SI (NON (EXIS TCOULA)) ; 112 : TCOULA = ZERO MOBET 'STRESSES' ; 113 : MESS '--- Pas de donnees d Instants de coulage du beton ' ; 114 : MESS ' On prend ZERO par defaut ' ; 115 : FINSI ; 116 : SI (NON (EXIS FLOT1)) ; 117 : MESS '--- Pas de donnees pour l humidite du milieu environnant' ; 118 : MESS ' pour la formulation BPEL91 - On quitte la procedure' ; 119 : QUIT PFLUAGE ; 120 : FINSI ; 121 : * 122 : MESS 'Instant de calcul du fluage :' TEMPCOUR 'jours' ; 123 : MESS 'Instant de coulage du beton :' 124 : (MINI TCOULA) (MAXI TCOULA) 'jours' ; 125 : MESS 'Age du beton au chargement :' 126 : (MINI AGECHAR) (MAXI AGECHAR) 'jours' ; 127 : * 128 : HUMREL = FLOT1 ; 129 : EI28 = 11000.E6 * ((FC28 /1.E6) ** (1. / 3.)) ; 130 : * 131 : MESS 'Module elastique du beton a 28 jours en Pa :' EI28 ; 132 : * 133 : * RHOS ETANT UN CHAMPS PAR ELEMENT ATTENTION AU CALCUL DE KS 134 : * 135 : * 136 : KS = 1.*((CHAMUNIT + (20.*RHOS))**-1.) ; 137 : * MESS 'FLUAGE MINI DE KS' (MINI KS) 'MAXI KS' (MAXI KS) ; 138 : * VBID = RHOS**-1. ; 139 : * VBID1 = RHOS1 * VBID ; 140 : * KS0 = VBID1 + (20. * RHOS) ; 141 : * KS = KS0**-1. ; 142 : * 143 : * ------- Gtv 144 : * 145 : KE = 0.4*CHAMUNIT ; 146 : * MESS 'FLUAGE MINI DE KE' (MINI KE) 'MAXI KE' (MAXI KE) ; 147 : * 148 : * ------- KC = Gech(rm) 149 : * 150 : KC = (((120. - HUMREL)/30.) * CHAMUNIT) + 151 : (((2./3.)*(100. - HUMREL))*(((20.*CHAMUNIT) + RAYSECH)**-1.)) ; 152 : * MESS 'FLUAGE MINI DE KC' (MINI KC) 'MAXI KC' (MAXI KC) ; 153 : * 154 : * ------- Gac(w) 155 : * 156 : K_T1 = 100.*(((100.*CHAMUNIT) + AGECHAR)**-1.) ; 157 : * MESS 'FLUAGE MINI DE KT' (MINI K_T1) 'MAXI KT' (MAXI K_T1) ; 158 : * 159 : * ------- Ffluage(t-taux) avec taux = AGECHAR 160 : * 161 : AGEBETON = (TEMPCOUR * CHAMUNIT) - TCOULA ; 162 : NUMER = (AGEBETON - AGECHAR)**0.5 ; 163 : CONST = RAYSECH**0.5 ; 164 : DENOM = NUMER + (5.*CONST) ; DENOM = DENOM**(-1) ; 165 : F_DT = NUMER * DENOM ; 166 : * MESS 'FLUAGE MINI DE ft' (MINI F_DT) 'MAXI FT' (MAXI F_DT) ; 167 : * 168 : * -- KS*(KE + (KC*K_T1)) = KS*(Gtv + Gech(rm)*Gac(w) = Kfluage(w) 169 : * -- E_BETON = E_BETON28 170 : * 171 : * -- MCHAMC = (1./E_BETON28)* Kfluage(taux)*Ffluage(t-taux) 172 : * 173 : MCHAMC = (EI28 ** (-1.))*KS*(KE + (KC * K_T1)) * F_DT ; 174 : MESS 'Deformation de Fluage total unitaire : ' 175 : (MINI MCHAMC) (MAXI MCHAMC) ; 176 : FINSI ; 177 : * 178 : SI (EGA NOMREG 'EC2') ; 179 : E_BETON = CHAN STRESSES MOBET 180 : (CHAN TYPE (EXCO YOUN MABET SCAL) 'CONTRAINTES') ; 181 : HUMREL = FLOT1 ; 182 : FCM = FLOT2 ; 183 : * 184 : PHIHR = CHAMUNIT + ((1. - (HUMREL/100.)) / 185 : (0.46*(RAYSECH/100.)**(1./3.))) ; 186 : BETAFCM = 5.3/((FCM/10.)**0.5) ; 187 : BETATAU = 1./((0.1*CHAMUNIT) + (AGECHAR**0.2)) ; 188 : BETAH = 1500.*CHAMUNIT ; 189 : BETACIN = (((TEMPCOUR*CHAMUNIT) - AGECHAR) / 190 : (((TEMPCOUR*CHAMUNIT) - AGECHAR) + BETAH))**0.3 ; 191 : * 192 : MCHAMC = (E_BETON**(-1.))*PHIHR*BETAFCM*BETATAU*BETACIN ; 193 : FINSI ; 194 : * 195 : SI (EGA NOMREG 'LG') ; 196 : * 197 : SI (NON (EGA (TYPE FLOT1) 'FLOTTANT')) ; 198 : KSA = 1./(1. + (16.*1.071E-2)) ; 199 : KSIGP = 1./(1. + (20.*3.790E-3)) ; 200 : SINON ; 201 : KSA = FLOT1 ; 202 : KSIGP = FLOT2 ; 203 : FINSI ; 204 : * 205 : K2 = (2.*160./(3.*RAYSECH))**2. ; 206 : * 207 : RAPVTAU = (0.1 + (28.**0.2))/((0.1*CHAMUNIT) + (AGECHAR**0.2)) ; 208 : * 209 : EEPREUVE = 31500.E6 ; 210 : EPSFP = 65.63E-6 ; 211 : AFP = 0.4969242 ; 212 : BFPJ = 21.97750113*CHAMUNIT ; 213 : EPSFT = 262.6E-6 ; 214 : AFT = 0.30071 ; 215 : BFTJ = 17.39335116*CHAMUNIT ; 216 : * 217 : FPTMTO = EPSFP*EEPREUVE* 218 : (((TEMPCOUR*CHAMUNIT) - AGECHAR)**AFP)/ 219 : ((((TEMPCOUR*CHAMUNIT) - AGECHAR)**AFP) + BFPJ) ; 220 : * 221 : FTK2TMTO = EPSFT*EEPREUVE* 222 : ((K2*((TEMPCOUR*CHAMUNIT) - AGECHAR))**AFT)/ 223 : (((K2*((TEMPCOUR*CHAMUNIT) - AGECHAR))**AFT) + BFTJ) ; 224 : * 225 : FPK2TMTO = EPSFP*EEPREUVE* 226 : ((K2*((TEMPCOUR*CHAMUNIT) - AGECHAR))**AFP)/ 227 : (((K2*((TEMPCOUR*CHAMUNIT) - AGECHAR))**AFP) + BFPJ) ; 228 : * 229 : MCHAMC = 1.E-6 * (E_BETON**(-1.))*KSA*KSIGP*RAPVTAU* 230 : (FPTMTO + FTK2TMTO - FPK2TMTO) ; 231 : FINSI ; 232 : * 233 : * 234 : * --- Construction du tenseur de deformation differee de fluage 235 : * 236 : * | A D E | 237 : * EPSFLU = | D B F | 238 : * | E F C | 239 : * 240 : EPSFLUXX = EXCO SCAL ((EXCO SMXX SIGEPS 'SCAL')*MCHAMC) EPXX ; 241 : EPSFLUYY = EXCO SCAL ((EXCO SMYY SIGEPS 'SCAL')*MCHAMC) EPYY ; 242 : EPSFLUZZ = EXCO SCAL ((EXCO SMZZ SIGEPS 'SCAL')*MCHAMC) EPZZ ; 243 : EPSFLUXY = EXCO SCAL (2.*(EXCO SMXY SIGEPS 'SCAL')*MCHAMC) GAXY ; 244 : EPSFLUXZ = EXCO SCAL (2.*(EXCO SMXZ SIGEPS 'SCAL')*MCHAMC) GAXZ ; 245 : EPSFLUYZ = EXCO SCAL (2.*(EXCO SMYZ SIGEPS 'SCAL')*MCHAMC) GAYZ ; 246 : * 247 : EPSFLU = EPSFLUXX ET EPSFLUYY ET EPSFLUZZ ET 248 : EPSFLUXY ET EPSFLUXZ ET EPSFLUYZ ; 249 : EPSFLU = CHAN 'TYPE' EPSFLU 'DEFORMATIONS' ; 250 : * 251 : MESS 'Deformation de Fluage (min/max) : ' 252 : (MINI (EPSFLUXX ET EPSFLUYY ET EPSFLUZZ ET 253 : ((EPSFLUXY ET EPSFLUXZ ET EPSFLUYZ) / 2.))) 254 : (MAXI (EPSFLUXX ET EPSFLUYY ET EPSFLUZZ ET 255 : ((EPSFLUXY ET EPSFLUXZ ET EPSFLUYZ) / 2.))) ; 256 : * 257 : FINP EPSFLU ; 258 : 259 :
© Cast3M 2003 - All rights reserved.
Disclaimer