1 : * RETRAIT PROCEDUR CHAT 05/10/11 21:16:03 5209 2 : DEBP RETRAIT NOMREG*MOT MOBET*MMODEL RAYSECH*MCHAML 3 : TCOULB*FLOTTANT TEMPCOUR*FLOTTANT 4 : FLOT1*FLOTTANT RHOS*MCHAML FC28/FLOTTANT ; 5 : * 6 : * RAYSECH : Rayon de sechage en cm 7 : * TCOULA : Instants de coulage du beton en jours 8 : * TEMPCOUR : Instant courant en jours 9 : * FLOT1 : Humidite du milieu environnant (formulation BPEL) 10 : * RHOS : Taux d'armatures 11 : * FC28 : Resistance caracteristique a 28 jours en Pa (Defaut : 40E6) 12 : * 13 : CHAMUNIT = MANU CHML MOBET SCAL 1. TYPE 'CONTRAINTES' STRESSES ; 14 : TCOULA=TCOULB*CHAMUNIT; 15 : FLOT2=RHOS; 16 : * 17 : RAYSECH = EXCO EPAI RAYSECH SCAL ; 18 : * 19 : * MODIFICATIONS B. CIREE DU 25/02/03 20 : * AJOUT DES FORMULATIONS BPEL MISE A JOUR 99 21 : * 22 : SI (EGA NOMREG 'BPEL99') ; 23 : * 24 : MESS 'Calcul du Retrait avec la formulation BPEL99' ; 25 : MESS '--------------------------------------------' ; 26 : * 27 : SI (NON (EXIS FLOT1)) ; 28 : MESS '--- Pas de donnees pour l humidite du milieu environnant' ; 29 : MESS ' pour la formulation BPEL99 - On quitte la procedure' ; 30 : QUIT RETRAIT ; 31 : FINSI ; 32 : SI (NON (EXIS FC28)) ; 33 : MESS '--- Pas de Resistance caracteristique du beton a 28 jours' ; 34 : MESS ' On prend 40.E6 Pa par defaut ' ; 35 : FC28 = 40.E6 ; 36 : FINSI ; 37 : * 38 : MESS 'Instant de calcul du retrait :' TEMPCOUR 'jours' ; 39 : MESS 'Instant de calcul de coulage du beton :' 40 : (MINI TCOULA) (MAXI TCOULA) 'jours' ; 41 : HUMREL = FLOT1 ; 42 : KFC28 = 18. ; 43 : * 44 : * --- RETRAIT ENDOGENE 45 : * 46 : EPSRE = (FC28 - 20.E6) * 1.E-12 47 : * ((2.8 * CHAMUNIT) 48 : - (1.1 * (EXP (((TEMPCOUR*CHAMUNIT) - TCOULA) / -96.)))); 49 : * 50 : * 51 : * --- RETRAIT DE DESSICCATION 52 : * 53 : A1 = KFC28 * ((72. * (EXP (-0.046E-6 * FC28))) + (75. - HUMREL)) 54 : * 1.E-6 ; 55 : A2 = CHAMUNIT + ((8.4 * RAYSECH * RAYSECH) * 56 : (((TEMPCOUR*CHAMUNIT) - TCOULA) ** -1.) ) ; 57 : EPSRD = A1 * (A2 ** -1.) ; 58 : * 59 : * 60 : * --- RETRAIT TOTAL 61 : * 62 : MCHAMC = (-1.) * ((EPSRE * CHAMUNIT) + EPSRD) 63 : * ((CHAMUNIT + (15. * RHOS)) ** -1.) ; 64 : * 65 : MESS 'retrait endogene : ' (MINI EPSRE) (MAXI EPSRE) ; 66 : MESS 'retrait de dessiccation :' (MINI EPSRD) (MAXI EPSRD) ; 67 : MESS 'retrait total : ' (MINI MCHAMC) (MAXI MCHAMC) ; 68 : * 69 : FINSI ; 70 : * 71 : SI (EGA NOMREG 'BPEL91') ; 72 : * 73 : MESS 'Calcul du Retrait avec la formulation BPEL91' ; 74 : MESS '--------------------------------------------' ; 75 : * 76 : SI (NON (EXIS FLOT1)) ; 77 : MESS '--- Pas de donnees pour l humidite du milieu environnant' ; 78 : MESS ' pour la formulation BPEL91 - On quitte la procedure' ; 79 : QUIT RETRAIT ; 80 : FINSI ; 81 : * 82 : HUMREL = FLOT1 ; 83 : KS = 1.*((CHAMUNIT + (20. * RHOS))** -1. ) ; 84 : * 85 : * ____ MODIFICATION DU 30 MARS 2001 F. PICARD 86 : * 87 : EPSI01 = (100. - HUMREL) ; 88 : A1 = (6.*CHAMUNIT) ; 89 : B1 = (10.*CHAMUNIT) ; 90 : C1 = (3.*RAYSECH) ; 91 : B1 = B1 + C1 ; 92 : B1 = B1**-1. ; 93 : B1 = 80. * B1 ; 94 : EPSI02 = A1 + B1 ; 95 : EPSI0 = EPSI01 * EPSI02 * 1.E-6 ; 96 : * 97 : NUMER = ((TEMPCOUR*CHAMUNIT) - TCOULA) ; 98 : DENOM = ((TEMPCOUR*CHAMUNIT) - TCOULA) + (9.*RAYSECH) ; 99 : DENOM = DENOM**(-1) ; 100 : R_T = NUMER * DENOM ; 101 : 102 : * 103 : * --- DEFORMATION DIFFEREE DE RETRAIT : SCALAIRE 104 : * 105 : MCHAMC = (-1.) * KS*EPSI0*R_T ; 106 : MESS 'retrait total : ' (MINI MCHAMC) (MAXI MCHAMC) ; 107 : * 108 : FINSI ; 109 : * 110 : SI (EGA NOMREG 'EC2') ; 111 : * 112 : HUMREL = FLOT1 ; 113 : FCM = FLOT2 ; 114 : EPSIFCM = (160. + (10.*5.*(9. - (FCM/10.))))*1.E-6 ; 115 : SI (HUMREL < 99.) ; 116 : BETAHR = -1.55*(1. - ((HUMREL/100.)**3.)) ; 117 : SINON ; 118 : BETAHR = 0.25 ; 119 : FINSI ; 120 : * 121 : ----------- MODIFICATION DU 29 MARS 2001 F. PICARD 122 : * 123 : XX1 = ((TEMPCOUR*CHAMUNIT) - TCOULA) ; 124 : XX2 = ((350.*((RAYSECH/100.)**2.)) + 125 : (TEMPCOUR*CHAMUNIT) - TCOULA)**0.5 ; 126 : XX2 = XX2** -1. ; 127 : BETACIN = XX1 * XX2 ; 128 : 129 : MCHAMC = EPSIFCM*BETAHR*BETACIN ; 130 : * 131 : FINSI ; 132 : * 133 : SI (EGA NOMREG 'LG') ; 134 : * 135 : SI (NON (EGA (TYPE FLOT1) 'FLOTTANT')) ; 136 : KSA = 1./(1. + (16.*1.071E-2)) ; 137 : KRETRP = 1./(1. + (16.*3.790E-3)) ; 138 : SINON ; 139 : KSA = FLOT1 ; 140 : KRETRP = FLOT2 ; 141 : FINSI ; 142 : * 143 : K2 = (2.*160.*((3.*RAYSECH)**-1.))**2. ; 144 : * 145 : EPSRE = -81.355E-6 ; 146 : ARE = 0.88023 ; 147 : BREJ = 182.9674339*CHAMUNIT ; 148 : EPSRD = -581.7077E-6 ; 149 : ARD = 0.9688749 ; 150 : BRDJ = 36.51898484*CHAMUNIT ; 151 : * 152 : RET = EPSRE * (((TEMPCOUR*CHAMUNIT) - TCOULA)**ARE) * 153 : (((((TEMPCOUR*CHAMUNIT) - TCOULA)**ARE) + BREJ)**-1.) ; 154 : * 155 : RDK2T = EPSRD * ((K2*((TEMPCOUR*CHAMUNIT) - TCOULA))**ARD) 156 : ((((K2*((TEMPCOUR*CHAMUNIT) - TCOULA))**ARD) + BRDJ)**-1.) ; 157 : * 158 : MCHAMC = KSA*KRETRP*(RET + RDK2T) ; 159 : * 160 : FINSI ; 161 : * 162 : * --- CONSTRUCTION DU TENSEUR DE DEFORMATION DIFFEREE DIAGONALE 163 : * DE RETRAIT 164 : * 165 : * | A 0 0 | 166 : * EPSRET = | 0 A 0 | AVEC A = MCHAMC EN SCALAIRE C'EST A DIRE : 167 : * | 0 0 A | 168 : * 169 : * EPSRET(I,J) = MCHAMC*D(I,J), D(I,J) = 1 SI I = J ET 0 SINON 170 : * 171 : EPSRETXX = EXCO SCAL MCHAMC EPXX ; 172 : EPSRETYY = EXCO SCAL MCHAMC EPYY ; 173 : EPSRETZZ = EXCO SCAL MCHAMC EPZZ ; 174 : EPSRETXY = EXCO SCAL (0.*MCHAMC) GAXY ; 175 : EPSRETXZ = EXCO SCAL (0.*MCHAMC) GAXZ ; 176 : EPSRETYZ = EXCO SCAL (0.*MCHAMC) GAYZ ; 177 : * 178 : EPSRET = EPSRETXX ET EPSRETYY ET EPSRETZZ ET 179 : EPSRETXY ET EPSRETXZ ET EPSRETYZ ; 180 : EPSRET = CHAN 'TYPE' EPSRET 'DEFORMATIONS' ; 181 : * 182 : FINP EPSRET ; 183 :
© Cast3M 2003 - All rights reserved.
Disclaimer