Télécharger la plaquette Cast3M


fr en

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