Creep and shrinkage of concrete
Creep and shrinkage of concrete are two physical properties of concrete. The creep of concrete, which originates from the calcium silicate hydrates (C-S-H) in the hardened Portland cement paste (which is the binder of mineral aggregates), is fundamentally different from the creep of metals and polymers. Unlike the creep of metals, it occurs at all stress levels and, within the service stress range, is linearly dependent on the stress if the pore water content is constant. Unlike the creep of polymers and metals, it exhibits multi-months aging, caused by chemical hardening due  to hydration which stiffens the microstructure, and multi-year aging, caused by long-term relaxation of self-equilibrated micro-stresses in the nano-porous microstructure of the C-S-H. If concrete is fully dried, it does not creep, but it is next to impossible to dry concrete fully without severe cracking. Changes of pore water content due to drying or wetting processes cause significant volume changes of concrete in load-free specimens. They are called the shrinkage (typically causing strains between 0.0002 and 0.0005, and in low strength concretes even 0.0012) or swelling (< 0.00005 in normal concretes, <  0.00020 in high strength concretes). To separate shrinkage from creep, the compliance function formula_1, defined as the stress-produced strain formula_2 (i.e., the total strain minus shrinkage) caused at time t by a unit sustained uniaxial stress formula_3 applied at age formula_4, is measured as the strain difference between the loaded and load-free specimens. The multi-year creep evolves logarithmically in time (with no final asymptotic value), and over the typical structural lifetimes it may attain values 3 to 6 times larger than the initial elastic strain. When a deformation is suddenly imposed and held constant, creep causes relaxation  of critically produced elastic stress. After unloading, creep recovery takes place, but it is partial, because of aging. In practice, creep during drying is inseparable from shrinkage. The rate of creep increases with the rate of change of pore humidity (i.e., relative vapor pressure in the pores). For small specimen thickness, the creep during drying greatly exceeds the sum of the drying shrinkage at no load and the creep of a loaded sealed specimen (Fig. 1 bottom). The difference, called the drying creep or Pickett effect (or stress-induced shrinkage), represents a hygro-mechanical coupling between strain and pore humidity  changes. Drying shrinkage at high humidities (Fig. 1 top and middle) is caused mainly by compressive stresses in the solid microstructure which balance the increase in capillary tension and surface tension on the pore walls. At low pore humidities (<75%), shrinkage is caused by a decrease of the disjoining pressure across nano-pores less than about 3 nm thick, filled by adsorbed water. The chemical processes of Portland cement hydration lead to another type of shrinkage, called the autogeneous shrinkage, which is observed in sealed specimens, i.e., at no moisture loss. It is caused partly by chemical volume changes , but mainly by self-desiccation due to loss of water consumed by the hydration reaction. It amounts to only about 5% of the drying shrinkage in normal concretes, which self-desiccate to about 97% pore humidity. But it can equal the drying shrinkage in modern high-strength concretes with very low water-cement ratios, which may self-desiccate to as low as 75% humidity. The creep originates in the calcium silicate hydrates (C-S-H) of hardened Portland cement paste. It is caused by slips due to bond ruptures, with bond restorations at adjacent sites. The C-S-H is  strongly hydrophilic, and has a colloidal microstructure disordered from a few nanometers up. The paste has a porosity of about 0.4 to 0.55 and an enormous specific surface area, roughly 500 m/cm. Its main component is the tri-calcium silicate hydrate gel (3 CaO · 2 SiO · 3 HO, in short CSH). The gel forms particles of colloidal dimensions, weakly bound by van der Waals forces. The physical mechanism and modeling are still being debated. The constitutive material model in the equations that follow is not the only one available but has at present the strongest theoretical foundation and fits  best the full range of available test data. Stress–strain relation at constant environment. In service, the stresses in structures are < 50% of concrete strength, in which case the stress–strain relation is linear, except for corrections due to microcracking when the pore humidity changes. The creep may thus be characterized by the compliance function formula_1 (Fig. 2). As formula_4 increases, the creep value for fixed formula_7 diminishes. This phenomenon, called aging, causes that formula_8 depends not only on the time lag formula_7 but on both formula_10 and formula_4 separately. At variable stress formula_12, each stress increment formula_13 applied  at time formula_4 produces strain history formula_15. The linearity implies the principle of superposition (introduced by Boltzmann and for the case of aging, by Volterra). This leads to the (uniaxial) stress–strain relation of linear aging viscoelasticity: Here formula_16 denotes shrinkage strain formula_17 augmented by thermal expansion, if any. The integral is the Stieltjes integral, which admits histories formula_18 with jumps; for time intervals with no jumps, one may set formula_19 to obtain the standard (Riemann) integral. When history formula_20 is prescribed, then Eq.(1) represents a Volterra integral equation for formula_18. This equation is not analytically integrable for realistic forms  of formula_1, although numerical integration is easy. The solution formula_18 for strain formula_24 imposed at any age formula_25 (and for formula_26) is called the relaxation function formula_27. To generalize Eq. (1) to a triaxial stress–strain relation, one may assume the material to be isotropic, with an approximately constant creep Poisson ratio, formula_28. This yields volumetric and deviatoric stress–strain relations similar to Eq. (1) in which formula_1 is replaced by the bulk and shear compliance functions: At high stress, the creep law appears to be nonlinear (Fig. 2) but Eq. (1) remains applicable if the inelastic strain due to  cracking with its time-dependent growth is included in formula_30. A viscoplastic strain needs to be added to formula_30 only in the case that all the principal stresses are compressive and the smallest in magnitude is much larger in magnitude than the uniaxial compressive strength formula_32. In measurements, Young's elastic modulus formula_33 depends not only on concrete age formula_4 but also on the test duration because the curve of compliance formula_1 versus load duration formula_7 has a significant slope for all durations beginning with 0.001 s or less. Consequently, the conventional Young's elastic modulus should be obtained as  formula_37, where formula_38 is the test duration. The values formula_39 day and formula_40 days give good agreement with the standardized test of formula_33, including the growth of formula_33 as a function of formula_4, and with the widely used empirical estimate formula_44 formula_45 formula_46. The zero-time extrapolation formula_47 happens to be approximately age-independent, which makes formula_48 a convenient parameter for defining formula_1. For creep at constant total water content, called the basic creep, a realistic rate form of the uniaxial compliance function (the thick curves in Fig. 1 bottom) was derived from the solidification theory: where formula_50; formula_51 =  flow viscosity, which dominates multi-decade creep; formula_52 = load duration; formula_53 = 1 day, formula_54, formula_55; formula_56 = volume of gel per unit volume of concrete, growing due to hydration; and formula_57 = empirical constants (of dimension formula_58). Function formula_59 gives age-independent delayed elasticity of the cement gel (hardened cement paste without its capillary pores) and, by integration, formula_60. Integration of formula_61 gives formula_1 as a non-integrable binomial integral, and so, if the values of formula_1 are sought, they must be obtained by numerical integration or by an approximation formula (a good formula exists). However, for computer  structural analysis in time steps, formula_1 is not needed; only the rate formula_61 is needed as the input. Equations (3) and (4) are the simplest formulae satisfying three requirements: 1) Asymptotically for both short and long times formula_52, formula_61, should be a power function of time; and 2) so should the aging rate, given by formula_68) (power functions are indicated by self-similarity conditions); and 3) formula_69 (this condition is required to prevent the principle of superposition from giving non-monotonic recovery curves after unloading which are physically objectionable). Creep at variable environment. At variable mass formula_70 of evaporable (i.e. , not chemically bound) water per unit volume of concrete, a physically realistic constitutive relation may be based on the idea of microprestress formula_71, considered to be a dimensionless measure of the stress peaks at the creep sites in the microstructure. The microprestress is produced as a reaction to chemical volume changes and to changes in the disjoining pressures acting across the hindered adsorbed water layers in nanopores (which are < 1 nm thick on the average and at most up to about ten water molecules, or 2.7 nm, in thickness), confined between the C-S-H sheets. The disjoining pressures  develop first due to unequal volume changes of hydration products. Later, they relax due to creep in the C-S-H so as to maintain thermodynamic equilibrium (i.e., equality of chemical potentials of water) with water vapor in the capillary pores, and build up due to any changes of temperature or humidity in these pores. The rate of bond breakages may be assumed to be a quadratic function of the level of microprestress, which requires Eq. (4) to be generalized as A crucial property is that the microprestress is not appreciably affected by the applied load (since pore water is  much more compressible than the solid skeleton and behaves like a soft spring coupled in parallel with a stiff framework). The microprestress relaxes in time and its evolution at each point of a concrete structure may be solved from the differential equation where formula_72 = positive constants (the absolute value ensures that formula_71 could never become negative). The microprestress can model the fact that drying and cooling, as well as wetting and heating, accelerate creep. The fact that changes of formula_70 or formula_75 produce new microprestress peaks and thus activate new creep sites explains the drying creep effect. A part  of this effect, however, is caused by the fact that microcracking in a companion load-free specimen renders its overall shrinkage smaller than the shrinkage in an uncracked (compressed) specimen, thus increasing the difference between the two (which is what defines creep). The concept of microprestress is also needed to explain the stiffening due to aging. One physical cause of aging is that the hydration products gradually fill the pores of hardened cement paste, as reflected in function formula_76 in Eq. (3). But hydration ceases after about one year, yet the effect of the age at loading formula_4 is strong  even after many years. The explanation is that the microstress peaks relax with age, which reduces the number of creep sites and thus the rate of bond breakages. At variable environment, time formula_10 in Eq. (3) must be replaced by equivalent hydration time formula_79 where formula_80 = decreasing function of formula_75 (0 if formula_82 about 0.8) and formula_83 formula_84. In Eq. (4), formula_85 must be replaced by formula_86 where formula_87 = reduced time (or maturity), capturing the effect of formula_75 and formula_89 on creep viscosity; formula_90 = function of formula_75 decreasing from 1 at formula_92 to 0 at formula_93; formula_94 , formula_95 5000 K. The evolution of humidity profiles formula_96 (formula_97 = coordinate vector) may be approximately considered as uncoupled from the stress and deformation problem and may be solved numerically from the diffusion equation formula_98 div[formula_99grad formula_100} where formula_101 = self-desiccation caused by hydration (which reaches about 0.97 in normal concretes and about 0.80 in high strength concretes), formula_102 = diffusivity, which decreases about 20 times as formula_75 drops from 1.0 to 0.6. The free (unrestrained) shrinkage strain rate is, approximately, where formula_104 = shrinkage coefficient. Since the formula_105-values at various points are incompatible, the calculation of the  overall shrinkage of structures as well as test specimens is a stress analysis problem, in which creep and cracking must be taken into account. For finite element structural analysis in time steps, it is advantageous to convert the constitutive law to a rate-type form. This may be achieved by approximating formula_59 with a Kelvin chain model (or the associated relaxation function with a Maxwell chain model). The history integrals such as Eq. 1 then disappear from the constitutive law, the history being characterized by the current values of the internal state variables (the partial strains or stresses of the  Kelvin or Maxwell chain). Conversion to a rate-type form is also necessary for introducing the effect of variable temperature, which affects (according to the Arrhenius law) both the Kelvin chain viscosities and the rate of hydration, as captured by formula_107. The former accelerates creep if the temperature is increased, and the latter decelerates creep. Three-dimensional tensorial generalization of Eqs. (3)-(7) is required for finite element analysis of structures. Approximate cross-section response at drying. Although multidimensional finite element calculations of creep and moisture diffusion are nowadays feasible, simplified one-dimensional analysis of concrete beams or girders based on  the assumption of planar cross sections remaining planar still reigns in practice. Although (in box girder bridges) it involves deflection errors of the order of 30%. In that approach, one needs as input the average cross-sectional compliance function formula_108 (Fig. 1 bottom, light curves) and average shrinkage function formula_109 of the cross section (Fig. 1 left and middle) (formula_110 = age at start of drying). Compared to the point-wise constitutive equation, the algebraic expressions for such average characteristics are considerably more complicated and their accuracy is lower, especially if the cross section is not under centric compression. The  following approximations have been derived and their coefficients optimized by fitting a large laboratory database for environmental humidities formula_111 below 98%: where formula_112 = effective thickness, formula_113 = volume-to-surface ratio, formula_114 = 1 for normal (type I) cement; formula_115 = shape factor (e.g., 1.0 for a slab, 1.15 for a cylinder); and formula_116, formula_117 = constant; formula_118 (all times are in days). Eqs. (3) and (4) apply except that formula_119 must be replaced by where formula_120 and formula_121. The form of the expression for shrinkage halftime formula_122 is based on the diffusion theory. Function 'tanh' in Eq. 8  is the simplest function satisfying two asymptotic conditions ensuing from the diffusion theory: 1) for short times formula_123, and 2) the final shrinkage must be approached exponentially. Generalizations for the temperature effect exist, too. Empirical formulae have been developed for predicting the parameter values in the foregoing equations on the basis of concrete strength and some parameters of the concrete mix. However, they are very crude, leading to prediction errors with the coefficients of variation of about 23% for creep and 34% for drying shrinkage. These high uncertainties can be drastically reduced by updating certain coefficients of the formulae according  to short-time creep and shrinkage tests of the given concrete. For shrinkage, however, the weight loss of the drying test specimens must also be measured (or else the problem of updating formula_124 is ill-conditioned). A fully rational prediction of concrete creep and shrinkage properties from its composition is a formidable problem, far from resolved satisfactorily. Engineering applications. The foregoing form of functions formula_125 and formula_126 has been used in the design of structures of high creep sensitivity. Other forms have been introduced into the design codes and standard recommendations of engineering societies. They are simpler though less realistic , especially for multi-decade creep. Creep and shrinkage can cause a major loss of prestress. Underestimation of multi-decade creep has caused excessive deflections, often with cracking, in many of large-span prestressed segmentally erected box girder bridges (over 60 cases documented). Creep may cause excessive stress and cracking in cable-stayed or arch bridges, and roof shells. Non-uniformity of creep and shrinkage, caused by differences in the histories of pore humidity and temperature, age and concrete type in various parts of a structures may lead to cracking. So may interactions with masonry or with steel parts, as in  cable-stayed bridges and composite steel-concrete girders. Differences in column shortenings are of particular concern for very tall buildings. In slender structures, creep may cause collapse due to long-time instability. The creep effects are particularly important for prestressed concrete structures (because of their slenderness and high flexibility), and are paramount in safety analysis of nuclear reactor containments and vessels. At high temperature exposure, as in fire or postulated nuclear reactor accidents, creep is very large and plays a major role. In preliminary design of structures, simplified calculations may conveniently use the dimensionless creep coefficient formula_127 = formula_128. The  change of structure state from time formula_129 of initial loading to time formula_10 can simply, though crudely, be estimated by quasi-elastic analysis in which Young's modulus formula_33 is replaced by the so-called age-adjusted effective modulus formula_132. The best approach to computer creep analysis of sensitive structures is to convert the creep law to an incremental elastic stress–strain relation with an eigenstrain. Eq. (1) can be used but in that form the variations of humidity and temperature with time cannot be introduced and the need to store the entire stress history for each finite element  is cumbersome. It is better to convert Eq. (1) to a set of differential equations based on the Kelvin chain rheologic model. To this end, the creep properties in each sufficiently small time step may be considered as non-aging, in which case a continuous spectrum of retardation moduli of Kelvin chain may be obtained from formula_125 by Widder's explicit formula for approximate Laplace transform inversion. The moduli formula_134 (formula_135) of the Kelvin units then follow by discretizing this spectrum. They are different for each integration point of each finite element in each time step. This way the creep  analysis problem gets converted to a series of elastic structural analyses, each of which can be run on a commercial finite element program. For an example see the last reference below.