Acta Polytechnica CTU Proceedings doi:10.14311/APP.2018.15.0012 Acta Polytechnica CTU Proceedings 15:12–19, 2018 © Czech Technical University in Prague, 2018 available online at http://ojs.cvut.cz/ojs/index.php/app COMPARISON OF DRYING SHRINKAGE AND DRYING CREEP KINETICS IN CONCRETE Lenka Dohnalová∗, Petr Havlásek Department of Mechanics, Faculty of Civil Engineering, Czech Technical University in Prague, Thákurova 7, 166 29 Prague 6, Czech Republic ∗ corresponding author: lenka.dohnalova@fsv.cvut.cz Abstract. The aim of this paper is to show and compare the time evolution of drying shrinkage and drying creep in concrete from three different perspectives. The first one analyzes the basic equations defined in the most common design codes and prediction models for the description of the long-term behavior of concrete (ACI 209, EC2, Model Code 2010, B3, B4). Next, the evolution of drying creep and shrinkage is examined by processing suitable experimental data available from the database developed at the Northwestern University. Finally, the last point of view investigates the results obtained from the finite element simulations employing the material point approach, in particular, the material model based on the Microprestress-Solidification theory. Keywords: Drying creep, shrinkage, concrete, design code, prediction model. 1. Introduction Concrete drying is accompanied by its gradual con- tractive volume changes referred to as shrinkage and additional compliance called drying creep or the Pick- ett effect. On contrary to basic creep (i.e. creep of hygrally sealed specimens at room temperature), both drying shrinkage and drying creep are believed to be bounded. The complex interplay between shrinkage and creep determines the evolution of stresses and so its thor- ough analysis can help to identify potential cracking and thus can contribute to estimation of the lifetime of the concrete structure. This is the case for members restrained either externally (statically indeterminate structures) or internally (nonuniform stress field gen- erated by the steep gradient of the profile of relative humidity) but very often by the combination of both. This paper summarizes and extends the final thesis of the first author [1]. The first part of this paper compares the time evo- lution of drying creep and drying shrinkage according to the most common design codes and prediction mod- els. These models operate on the cross-sectional level and thus they provide the estimated average values of strain and compliance which are assumed to be uniform across the cross-section. Nowadays, most of the prediction models use as the input parameters the recipe of the concrete mixture, the 28-day compressive strength fc, and the duration of curing (or the onset of drying). In order to estimate the rate of drying (and so the rate of shrinkage and drying creep), it is necessary to provide the model with the size (usually by the quantity called “equivalent thickness” com- puted as the ratio of volume to drying surface) and the shape of the cross-section. The second part of the paper examines and pro- cesses the available experimental data gathered in the database developed and maintained at the North- western University (NU) [2]. The next part compares the drying shrinkage and drying creep kinetics ob- tained from the results of the finite element simula- tions. These simulations operate on the material point level which implicates the non-uniform stress field. If such model is correctly calibrated, it is applicable not only for the analysis of complex and creep-sensitive structures, but also for the assessment of potential surface or internal cracking stemming from the severe drying conditions or excessive temperature generated by cement hydration in massive structures. The results are compared and discussed in the last part of the paper. 2. Design Codes and Prediction Models Two design codes (ACI 209.2R–08 [3] and Eurocode 2 [4]) and three prediction models for the long–term behavior of concrete (B3 model [5], B4 model [6] and fib Model Code 2010 [7]) were selected for a detailed analysis of time evolution of drying creep and drying shrinkage. The investigated prediction models split concrete creep additively into basic creep and drying creep components allowing to evaluate drying creep directly. On the other hand, the analyzed design codes (ACI 209, Eurocode 2) provide only the formula for the total compliance; therefore, the drying creep compliance was obtained by subtracting the compliance at sealed conditions from the total compliance corresponding to the drying conditions. The total shrinkage is usually considered as the sum of autogenous and drying shrinkage. The B3 model neglects the autogenous shrinkage, for this reason the drying shrinkage is set equal to the total shrinkage. 12 http://dx.doi.org/10.14311/APP.2018.15.0012 http://ojs.cvut.cz/ojs/index.php/app vol. 15/2018 Comparison of drying shrinkage and drying creep kinetics in concrete The ACI 209 does not neglect the autogenous shrink- age but provides only the formula for total shrinkage which is nil at 100 % relative humidity. Therefore, similarly to the B3 model, the total shrinkage is taken as the drying shrinkage. Even though the individual models use different functions to define the time evolution of drying creep and shrinkage, the shape of the resulting curves is very similar. Bigger differences are observed in the final values but these are not in the scope of the present paper. Next, we present a brief description of the design and prediction models used in this study. Throughout the different models εsh,d stands for the drying shrinkage and henv for the ambient relative humidity (expressed as decimal). The expressions for the drying creep compliance Jd are for simplicity presented with t′ = t0, i.e. with the age at loading t′ coinciding with the onset of drying t0. The duration of loading/drying is denoted as t̂ = t− t′ = t− t0. In some models the evolution of drying is derived from V/S which is the ratio of the volume of the concrete member to its surface exposed to the ambient relative humidity. B3 model was developed by Bažant and co–workers at the Northwestern University and uses, similarly to the other prediction models, empirical formulae to de- termine the time development of creep and shrinkage. The parameters in the model can be estimated based on the composition of the concrete mixture, the mean compressive strength and in the case of drying on the ambient relative humidity and on the size and shape of the concrete member. The autogenous shrinkage is neglected in the model and so the range of appli- cability is limited to the water-to-cement ratio above 0.35. The drying creep compliance and drying shrinkage are defined by the following two equations Jd(t̂) = q5  e8 ( (1−henv)·tanh √ t̂ τsh −1 ) − 1  0.5 (1) εsh(t̂) = −ε∞sh ·kh · tanh √ t̂ τsh (2) Here, τsh is the shrinkage half–time depending on the size and shape of the member, ε∞sh is the absolute value of the ultimate drying shrinkage at henv = 0, kh is a humidity-dependent factor, and q5 is an adjustable parameter. B4 model builds on the previous B3 model and uses more input parameters of the concrete mixture (e.g. aggregate type, additives and admixtures) as well as more elaborate formulae for the prediction of the model parameters. Nevertheless, the structure of basic equations of the B4 model is very similar to its predecessor. The autogenous shrinkage and drying shrinkage are treated in the B4 model separately. fib Model Code 2010 The next two equations present a simplified form of the drying creep compli- ance and shrinkage functions Jd(t̂, t′) = 1 E28 β∗dc [ t̂ βD + t̂ ] √t′ 3.5+ √ 5.29 t′ (3) εsh,d(t̂) = εsh,d0 βRH [ t̂ 0.035 D2 + t̂ ]0.5 (4) where E28 is the modulus of elasticity at the age of 28 days, D is the effective cross section thickness (twice the volume-surface ratio), βD is a size dependence factor, βRH is a humidity dependence factor, and β∗dc and εsh,d0 are model parameters depending on compressive strength, relative humidity and the age at loading. Eurocode 2 does not specify an explicit formula for the evaluation of the compliance function J from the creep coefficient ϕ, the conversion was done using a modified expression from Model Code 1990 [8]. J(t̂, t′) = 1 Ecm(t′) + ϕ 1.05Ecm(28) [ t̂ βH + t̂ ]0.3 (5) εsh,d(t̂) = kD ·εsh0 · t̂ t̂ + 0.04 ·D1.5 (6) where Ecm(t′) is the secant modulus of elasticity of concrete at the age of loading, βH is a parameter depending on relative humidity and member size, kD ·εsh0 expresses the ultimate value of drying shrink- age. ACI 209.2R–08 On contrary to EC2, the Amer- ican standard determines the compliance function from the creep coefficient and the secant modulus of elasticity at the age at loading Ecm(t′) J(t̂, t′) = 1 Ecm(t′) · [ 1 + t̂ψ d + t̂ψ ·ϕu ] (7) εsh(t̂) = εshu · t̂α f + t̂α (8) where d, ψ, f, and α are parameters, and ϕu and εshu are the ultimate creep coefficient and the ultimate shrinkage, respectively, both influenced by the ambient relative humidity. There are two distinct options how to select the model parameters defining the time evolution of drying creep and shrinkage. The first one, as simplified in [3] from the former version of the code [9], is to use the average and fixed values of the parameters, then ψ = 0.6, α = 1.0 and in general d 6= f. The second possibility, which reflects the size and shape of the structural member, is to set the exponents α = ψ = 1.0, and d = f which are determined from V/S of the drying member. 13 Lenka Dohnalová, Petr Havlásek Acta Polytechnica CTU Proceedings Basic parameters fc Ecm28 t0 t′ c w/c a/c Cem. henv 2V/S MPa GPa day day kg m−3 - - type - mm Example from ACI [3] 33.3 28.2 7 7 409 0.50 4.23 I 0.7 200 Experiment Troxell [10] 16.5 NA 28 28 320 0.59 5.67 I 0.5/0.7 45 Experiment Bryant [11] 50.1 29.8 8 8 390 0.47 5.09 I 0.6 150/75 ∗) Experiment Keeton [12] 45 25.9 8 8 452 0.46 3.74 III 0.2/0.5 65 ∗) 150 mm slab, 75 mm prism Table 1. Details of the experimental setup in the creep and shrinkage experiments. 10-1 100 101 102 103 104 0 10 20 30 40 - - 10-1 100 101 102 103 104 0 100 200 300 400 500 600 - - Figure 1. Drying creep compliance curves (top) and drying shrinkage curves (bottom) given by the ana- lyzed design codes and prediction models. The time evolution of drying creep and shrinkage shown in Fig. 1 was obtained with the parameters computed for the concrete mixture and the specimen properties defined in Appendix C of ACI 209 [3]. The parameters are listed in Tab. 1. The specimen is con- sidered as a slab with 200 mm thickness and moist cured for 7 days. The aggregate type used in the B4 model is the quartzite. The curves in Fig. 1 are evalu- ated for the duration of 100 years which is the typical design period of most of the bridges and engineering structures. Both figures show the horizontal time axis in the logarithmic scale which is more suitable for processes with initially very high rate and which gradually decreases. Fig. 2 presents the time evolution of normalized dry- 100 102 104 0 0.5 1 100 102 104 0 0.5 1 100 102 104 0 0.5 1 100 102 104 0 0.5 1 Figure 2. Time evolution of normalized drying creep and shrinkage according to the selected design codes and prediction models. The functions are normalized with respect to the ultimate value. Decreasing drying creep from EC2 model is physically unjustified and stems from the model formulation. ing creep, (J̄d, solid red lines) and drying shrinkage (ε̄sh,d, dashed blue lines) for the selected design codes and prediction models. The evolution is normalized with respect to the ultimate value. The presented response of the models is based on the input param- eters specified in the first line of Tab. 1. Naturally, the curves will be different for a modified composi- tion, dimensions or other conditions, but the dominant trends will prevail provided that the ambient relative humidity is kept constant (see the first part in Section 5). A different way of comparison is shown in Fig. 3 where the normalized drying creep is plotted against normalized drying shrinkage. However, here, the real kinetics of both phenomena is not reflected. It must be noted that in the present example (shown in Figs. 1–3), the ACI model uses the recommended fixed values of parameters (i.e. the model is size independent). 14 vol. 15/2018 Comparison of drying shrinkage and drying creep kinetics in concrete 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 - - Figure 3. Comparison of normalized drying shrink- age and normalized drying creep evolution according to the selected design codes and prediction models. Normalization is done with respect to the ultimate value. 3. Experimental Data The freely available experimental database developed at the Northwestern University [2] was selected as the source of suitable experimental evidence of dry- ing shrinkage and drying creep. Some experimental data comprised in this database were used for the calibration of the models presented in the previous section, and so it can be expected that the processed experimental data will very likely resemble the curves in Figs. 2 and 3. The database contains 61,930 values measured on 3,308 specimens (1,439 for creep and 1,869 for shrink- age). These values come from 362 experimental sur- veys (172 creep, 190 shrinkage). Some experimental surveys analyzed more than one concrete composition; in total it was 378 in the case of creep and 652 dif- ferent compositions in the case of shrinkage. Only 68 experimental surveys are relevant to this study be- cause only those used the same recipe of the concrete mixture for both creep and shrinkage specimens. For the present purpose it was necessary to further narrow the selection only to those experiments which tested simultaneously drying shrinkage on unloaded specimens, basic creep, and total creep matter. Ad- ditionally, the loading in the creep experiments had to be applied simultaneously with the beginning of drying of both loaded and unloaded specimens. The database contains only very few surveys which provide the information on autogenous shrinkage and sometimes it is difficult to distinguish between experi- ments carried out at sealed conditions and experiments where the specimens are immersed in water or ex- posed to very high relative humidity and thus exhibit swelling. In order to eliminate this deficiency an addi- tional criteria were set on the value of water-to-cement ratio (minimum value w/c = 0.4) and the compres- sive strength (maximum value fc = 60 MPa), then the autogenous shrinkage can be treated as negligible compared to the magnitude of the drying shrinkage. Only 97 sets from 9 experimental surveys met the aforementioned criteria. Unfortunately, it was nec- essary to further reduce the selection because of in- correct or missing pieces of information, or too short duration of the experiment. The duration was judged as sufficient if the rate of drying creep and shrinkage was minimal at the end of the experiment. This was satisfied only in 6 experimental surveys (with 71 spec- imens) out of which 3 representative are partially presented in this paper. The value of drying creep was determined as the difference between total and basic creep; when needed (in the case of different times of measurement), the value was computed using linear interpolation between two adjacent measurements. The drying creep and shrinkage was afterwards normalized with respect to the last value attained. The suitable experiments are summarized in the following paragraphs. Experiments of Troxell [10] adopted in this study are part of an extensive experimental survey focused on the influence of different aspects on the develop- ment of creep and shrinkage. The experiment started in 1928 and its duration exceeded 20 years. The data used in this study come from the experimental series #3 which examined the influence of relative humidity; the specimens were cylinders with 102 mm diameter and 356 mm height. The basic properties and the composition of the concrete mixture are specified in Tab. 1. The NU database does not contain data for t̂ < 12 days and thus completely misses the initial evolution of creep and shrinkage. Two triplets of specimens were suitable for the comparison of drying shrinkage and drying creep kinetics, one for 50 %, the second one for 70 % relative humidity. The results are shown in Fig. 4. Experiments of Bryant [11] was initially focused on study of creep and shrinkage in bridge structures. Only specimens with dimensions 150 × 150 × 600 mm were suitable for the present study because only this size was loaded at the age when the drying began which was not the case of the remaining sizes. The specimens should have realistically represented struc- tural members - slabs and pillars, to achieve this, one set of specimens had both ends and two opposite lat- eral sides hygrally sealed (specimens referred to as “slab”) which resulted in unidirectional drying. The other set was sealed only at the ends and thus the specimens were drying evenly from all lateral sides, this set is denoted as “prism”. 15 Lenka Dohnalová, Petr Havlásek Acta Polytechnica CTU Proceedings 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 Figure 4. Comparison of normalized drying shrink- age and normalized drying creep evaluated from the experiments of Troxell and Bryant. Only the points represent the experimental data. At the age of 8 days all specimens were transferred to the environment with 60 % relative humidity. Un- fortunately, during the experiment the aluminum foil used for sealing has corroded and especially the en- tirely sealed specimens started unintentionally dry- ing resulting into an increase in creep or shrinkage rate. For this reason the autogenous shrinkage was neglected and the total shrinkage was taken as the drying shrinkage. At the later stage of the experiment, the measurements of the basic creep (spoiled by the sealing leak) partially influenced the determination of drying creep. The processed results are for both groups (slabs and prisms) shown in Fig. 4. Experiments of Keeton [12] The suitable experi- mental data come from the part of research studying the influence of ambient relative humidity, specimen size, and the load level on creep and shrinkage. The research started in 1957 and lasted for 897 days. The data are not from the NU database but instead from the original paper (and are different). In several cases the increments of basic creep are larger than the incre- ments of total creep. This was the problem especially of smaller specimens with lower level of loading and obviously specimens exposed to higher of relative hu- midity. Additionally, both loaded and unloaded com- panion specimens were not sealed but were subjected to 100 % relative humidity. Owing to the aforemen- tioned reasons the smaller specimens, the specimens with low loading level and specimens drying at 75 % relative humidity are omitted from the present analy- sis. Fig. 5 shows processed results obtained from exper- iments on cylinders 152 × 457 mm drying at 20 % and 50 % relative humidity and loaded at 30 % and 40 % of 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 Figure 5. Comparison of normalized drying shrink- age and normalized drying creep evaluated from the experiments of Keeton at 20 % and 50 % relative hu- midity and two levels of loading. f′c is the compressive strength at the age of loading. Points are the experi- mental data. the compressive strength at the age of loading, f′c. In this figure, the higher initial level of shrinkage stems from the first measurement in which the total creep was smaller than the basic creep. 4. Data from Finite Elements Simulations Except for the prismatic specimens 150×150×600 mm shown in the previous section of this study, the Bryant’s experimental survey comprises the informa- tion on shrinkage and drying creep measured on geo- metrically specimens ranging from 100×100×400 mm to 400 × 400 × 1600 mm prepared from a single con- crete batch. For this reason this experiment was in [13] and [14] selected as an ideal opportunity for verifi- cation of the material model based on Microprestress- Solidification theory (MPS) [15]. It turned out that the original MPS gives the oppo- site size effect on drying creep than the experiment and for this reason the MPS model was modified [13, 14]. Here we present the comparison of drying creep and shrinkage for both approaches. Hereafter, the mod- ified version of MPS model is referred to as “MPS improved”. The presented data use the processed results obtained from the finite element simulations carried out in the open source FE package [16]. The problem was analyzed as one-way coupled with the structural problem depending on the moisture trans- port according to the Bažand and Najjar model [17] for diffusion of water vapor. In order to obtain the shrinkage and drying creep kinetics, it was necessary perform, similarly to the real experiment, 3 analy- 16 vol. 15/2018 Comparison of drying shrinkage and drying creep kinetics in concrete 100 102 104 0 0.5 1 100 102 104 0 0.5 1 100 102 104 0 0.5 1 100 102 104 0 0.5 1 Figure 6. The processed results of FE simulations of Bryant’s experiment showing the time evolution of normalized drying creep and shrinkage. The top two figures are for 150 mm slabs bottom two figures are for 150 × 150 mm prism. The curves are normalized with respect to the ultimate value. ses: basic creep, drying creep and shrinkage. The values of the material parameters in the models were selected in order to obtain a reasonably accurate fit of all specimens (“prisms” and “slabs” of different sizes), therefore, the agreement of the simulation of 150 mm specimen with the experimental data cannot be expected to be satisfactory. The normalized drying creep and shrinkage is for the original and improved MPS model and the two geometries shown in Fig. 6. 5. Discussion In general, according to the prediction models and design codes, the time evolution of drying creep and drying shrinkage depends dominantly on the ambient relative humidity and on the size of the concrete mem- ber, the influence of other factors is not so significant. When evaluating the differences in kinetics of drying creep and drying shrinkage, each model is sensitive only to certain input parameters. In the case of the B3 model, the only factor which affects the relationship between drying creep and dry- ing shrinkage is the ambient relative humidity. At low relative humidity, the drying creep is delayed be- hind shrinkage, but on contrary to this, at high levels (henv > 0.8), the delay is in shrinkage. However, the difference is in both cases relatively small, see Fig. 7. On the other hand, the fib MC2010 is completely insensitive to henv. The kinetics is influenced domi- nantly by the size of the member (its V/S), see Fig. 8. The larger the size the more is the curve in Fig. 2 bent downwards and the bigger is the delay of shrinkage behind drying creep. This is caused by the fact that with decreasing V/S the drying shrinkage becomes 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 - - Figure 7. The influence of different levels of ambi- ent relative humidity on the drying creep vs. drying shrinkage kinetics according to the B3 model. The re- maining parameters are from the first line in Tab. 1. 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 - - Figure 8. The influence of different V /S ratios on the drying creep vs. drying shrinkage kinetics according to the fib MC 2010 model. The remaining parameters are from the first line in Tab. 1 more accelerated than the drying creep. To certain extent the kinetics is also influenced by the onset of drying, but for reasonable values the differences are not substantial. The Eurocode 2 predicts, except for extreme com- binations of parameters, delay in shrinkage behind drying creep. The model is sensitive both to V/S a and henv. For smaller specimens and higher rela- tive humidity the drying creep curve becomes more non-monotonic. 17 Lenka Dohnalová, Petr Havlásek Acta Polytechnica CTU Proceedings 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 - - - - - Figure 9. The comparison of the experimental data, prediction models and FE simulations of the Bryant’s experiment, 2V /S = 150 mm (slab). ACI 209 is the exception to the other presented models. For first approach with the average recom- mended values (according to [9]) the kinetics of drying creep and shrinkage remains constant, with an initial small delay of shrinkage behind drying creep. Ac- cording to the second approach, which sets equal the exponents α and ψ as well as the parameters d and f, the kinetics of drying creep and drying shrinkage be- come identical (i.e. the function would be a diagonal in Fig. 3). As presented in Figs. 4 and 5, the selected suitable experimental data give the opposite trend compared the prediction models (except for B3/4 at lower rela- tive humidity and ACI at the later stage of drying). This means that in the experiments the drying shrink- age is a little bit ahead of the drying creep. The points in these figures correspond to the experimental data which are interpolated by the curves of the same color. Two issues have to be pointed out. First, the result- ing shape of the curve in the normalized plot shrinkage vs. drying creep depends on the duration of the ex- periment. Insufficiently long experiments can lead to misleading conclusions. Second, as demonstrated in Fig. 5, the readings should be documented more frequently in the initial phase of drying and the speci- mens should be instrumented as soon as possible. In the case of Keeton’s experiment, the first displayed data point corresponds to the drying duration 7 days and the shrinkage strain reaches almost 40 % of the last measured value (after approx. 3 years). (The first data points for t̂ = 1 day which would have yielded negative Jd is omitted here.) According to the computations with the original ver- sion of the material model based on the Microprestress- Solidification theory, the drying creep evolves in an 0 0.2 0.4 0.6 0.8 1 0 0.2 0.4 0.6 0.8 1 - - - - - Figure 10. The comparison of the experimental data, prediction models and FE simulations of the Bryant’s experiment, 2V /S = 75 mm (prism). unexpected fashion. First, as shown in the two left graphs in Fig. 6 and in Figs. 9 and 10, the drying creep evolves initially faster than drying shrinkage, then, suddenly, the creep stops and after certain time it continues with a considerable delay after shrinkage. This effect is more pronounced in the case of prismatic shape which dries more rapidly. The response of the improved MPS model is more promising, see the two right graphs in Fig.6. Similarly to the cross-sectional approaches, it gives a small ini- tial delay of shrinkage after drying creep, afterwards, the relationship is almost linear. However, the initial delay is reversed compared to the analyzed experimen- tal data. Very probably, in the case of the models working on the cross-sectional approach, the differences of the shrinkage and drying creep kinetics origin from their optimization and calibration procedure. Very likely the models were calibrated against the entire database with evenly distributed weights and without emphasizing that some specimens are companions belonging to the same experimental survey. 6. Conclusions The presented study analyzed the time evolution of drying shrinkage εsh,d and drying creep Jd from three different perspectives: the most common design and prediction models, experimental data measured on laboratory concrete specimens, and the finite element simulations. The following conclusions can be drawn: • Even though the design codes and prediction mod- els are calibrated on a very similar experimental database, Jd and εsh,d are described in each model 18 vol. 15/2018 Comparison of drying shrinkage and drying creep kinetics in concrete differently and their kinetics is different. Addition- ally, in each model, the relationship between the Jd and εsh,d depends on different factors. Three of the analyzed models, B3, B4, and ACI 209, give for the standard conditions comparable kinetics of these processes. However, in the case of fib MC2010 and EC2 the delay of εsh,d behind Jd is substantial. According to EC2, Jd is a non-monotonic function. • A closer inspection of the experimental database which contains a vast amount of data revealed that only very few experiments fulfill the needs of the present research. Unfortunately, most of the ex- periments does not last sufficiently long or certain essential pieces of information are missing. The an- alyzed data indicate that the Jd is initially slightly delayed behind εsh,d, which means that the trend is exactly opposite compared to the prediction models. • The improved MPS model gives more realistic kinet- ics of Jd and εsh,d compared to its predecessor which unexpectedly gives a significant decrease in drying creep rate in the middle of drying. Both models predict, similarly to the cross-sectional models, an initial small delay of εsh,d behind Jd. • Sufficiently long experiments on normal-strength concrete are vital for the development and proper calibration of the material models for creep and concrete, however, they are very scarce. Acknowledgements The authors gratefully acknowledge financial support from the Grant Agency of the Czech Technical University in Prague, project number SGS18/037/OHK1/1T/11. The numerical analyses have been performed with OOFEM, an open-source object-oriented finite element program [16, 18]. References [1] L. Dohnalová. Comparison of shrinkage and drying creep kinetics of concrete, Bc. Thesis. Czech Technical University in Prague, 2018. [2] M. Hubler, R. Wendner, Z. Bažant. Comprehensive database for concrete creep and shrinkage: analysis and recommendations for testing and recording. ACI 112:547–558, 2015. [3] ACI Committee 209. Guide for Modeling and Calculating Shrinkage and Creep in Hardened Concrete. ACI, 2008. [4] Český normalizační institut. Eurokód 2: Navrhování betonových konstrukcí - Část 1-1: Obecná pravidla a pravidla pro pozemní stavby. Český normalizační institut, 1996. [5] Z. Bažant, S. Baweja. Creep and shrinkage prediction model for analysis and design of concrete structures: Model B3. Adam Neville Symposium: Creep and Shrinkage - Structural Design Effects 2000. [6] RILEM Technical Committee TC-242-MDC (Z.P. Bažant, chair). Model B4 for creep, drying shrinkage and autogenous shrinkage of normal and high-strength concretes with multi-decade applicability (RILEM Technical Committee TC-242-MDC multi-decade creep and shrinkage of concrete: material model and structural analysis). Materials and Structures 48:753–770, 2015. [7] Fédération Internationale du Béton. Model Code 2010. No. vol. 65 in fib Bulletin. International Federation for Structural Concrete (fib), 2012. [8] Comité Euro–International du Béton (CEB). CEB–FIP Model Code 1990: Design Code. T. Telford, London, 1993. [9] A. C. 209. Prediction of Creep, Shrinkage, and Temperature Effects in Concrete Structures. American Concrete Institute, 1997. [10] G. Troxel, J. Raphael, R. Davis. Long-time creep and shrinkage tests of plain and reinforced concrete. Proc ASTM 58 pp. 1101–1120, 1958. [11] A. H. Bryant, C. Vadhanavikkit. Creep, shrinkage-size, and age at loading effects. ACI Materials Journal 84:117–123, 1987. [12] Keeton. Study of Creep in Concrete, Phase 1 (I-beam). U. S. Naval Civil Engineering Laboratory, Port Hueneme, California, 1965. [13] P. Havlásek. Creep and Shrinkage of Concrete Subjected to Variable Environmental Conditions, PhD. Thesis. Czech Technical University in Prague, 2014. [14] Z. Bažant, P. Havlásek, M. Jirásek. Microprestress-solidification theory: Modeling of size effect on drying creep. In N. Bicanic, H. Mang, G. Meschke, R. de Borst (eds.), Computational Modelling of Concrete Structures, pp. 749–758. CRC Press/Balkema, EH Leiden, The Netherlands, 2014. [15] Z. P. Bažant, A. P. Hauggaard, S. Baweja, F. J. Ulm. Microprestress solidification theory for concrete creep. I: Aging and drying effects. Journal of Engineering Mechanics 123:1188–1194, 1997. [16] B. Patzák. OOFEM home page, 2000. Http://www.oofem.org. [17] Z. P. Bažant, L. J. Najjar. Nonlinear water diffusion in nonsaturated concrete. Materials and Structures 5:3–20, 1972. doi:10.1007/BF02479073. [18] B. Patzák. OOFEM - an object-oriented simulation tool for advanced modeling of materials and structures. Acta Polytechnica 52(6):59–66, 2012. 19 http://dx.doi.org/10.1007/BF02479073 Acta Polytechnica CTU Proceedings 15:12–19, 2018 1 Introduction 2 Design Codes and Prediction Models 3 Experimental Data 4 Data from Finite Elements Simulations 5 Discussion 6 Conclusions Acknowledgements References