# Description of Concrete Creep under Time-Varying Stress Using Parallel Creep Curve

Authors: Yeong-Seong Park; Yong-Hak Lee; Youngwhan Lee
Journal: Advances in Materials Science and Engineering
 (2016)

---

## Abstract

An incremental format of creep model was presented to take account of the development of concrete creep due to loading at different ages. The formulation was attained by introducing a horizontal parallel assumption of creep curves and combining it with the vertical parallel creep curve of the rate of creep method to remedy the disadvantage of the rate of creep method that significantly underestimates the amount of creep strain, regardless of its simple format. Two creep curves were combined by introducing an ageing parameter whose value was obtained from two sets of time-dependent laboratory experiments on cylindrical specimens. The presented creep description takes the advantage that a single creep curve due to the initial loading describes the entire development of creep under the persistent change of creep-causing stress. Further, the creep formulation takes advantage of being consistent with the incremental format of age-dependent constitutive formulation. The performance of the presented creep equation was investigated with time-dependent laboratory experiments on cylindrical specimens and compared with the performances of four existing creep models.


---

## Body

## 1. Introduction

The restraint of creep as well as shrinkage strains causes the mechanical strain and becomes a source of persistent change in the creep-causing mechanism. This type of creep development forms a circulating loop with a mechanical significance that is the same as the creep strain developed under a time-varying stress history. Mathematical modeling of the creep mechanism is formulated in terms of the age and concrete properties at loading. Aside from the effect of time-varying stress history on the creep model, another aspect concerning the creep model is the formulation format where the creep model is combined with the age-dependent constitutive model that relates the change of stress with the mechanical strain developed due to the restraint of creep and shrinkage. When the creep model is combined with the constitutive model, a consistent formulation framework is required between the two models, such as the total time-based or incremental time-based formulation frameworks. The consistency requirement in the formulation framework is further extended to the global equilibrium equations to compute the nodal displacements whose formulation strongly depends on the type of constitutive descriptions. A number of studies regarding the formulation of global equilibrium equations have been presented in the total time-based format of the finite element analysis scheme and applied to the time-dependent behaviors of concrete structures [1–3]. In this paper, a creep model is formulated within the incremental framework of the formulation to be consistent with the incremental format of the constitutive model.Most studies regarding creep models have been performed to model the creep phenomena under various conditions of mix proportions, curing environments, ages of loading, and geometrical shape and dimensions [4–11]. However, focusing on the circulating loop phenomenon occurring in the time-dependent analysis of concrete structures, the creep model concentrates on depicting the creep strain developed under the time history loadings applied at different ages. A simple and robust method to model creep strain under time-varying stress history is a step-by-step method based on the principle of superposition. The formulation subdivides the entire stress history with stepwise variations of stress, which are called incremental stresses, applied at the end of a small time interval [11]. However, it is a difficult task to define the number of creep curves corresponding to the number of time intervals with stress increments. A great deal of effort has been made to model creep curves. The related studies are well introduced in the books of the ACI Committee 209 [5], Bažant and Wittmann [12], CEB [9], Ghali et al. [13], Gilbert [14], Gilbert and Ranzi [11], Jirasek and Bažant [15], and Rusch et al. [16]. Among the number of studies presented in the literature, four creep models draw attention to depict the circulating loop phenomenon: the rate of creep method (RCM) [17], step-by-step method [11], effective modulus method [18], and ageing coefficient method [19]. The performance of the models had been extensively studied by Gilbert and Ranzi [11] through a comparison of the mathematical formulations and numerical applications. Each model shows a different performance depending on the assumption made regarding the ageing effects of concrete loaded at different ages. Among these four models presented in the literature, special attention is given to RCM [17] because of its simple single-curve representation of creep behavior under a time-varying stress history. However, this method significantly underestimates the creep strain when it is applied to the creep problem subjected to persistent change in creep-causing stress. This paper presents a model to improve the disadvantage observed in RCM by introducing a horizontal parallel creep curve.RCM is often referred to as the parallel creep curve method because the creep curves due to the loadings at different ages are assumed to be parallel. The term “parallel” in this method means that the tangents of the creep curves are identical along a vertical line. However, the method suffers from a significant underestimation of creep strain because of the parallel assumption when the sustained loads are applied at different ages. To remedy this inherent limitation, a horizontal parallel creep curve assumption is introduced and combined with the vertical parallel creep assumption of RCM. This type of formulation provides a creep strain bounded between upper and lower limits, where the lower limit is defined by RCM and the upper limit is defined by the horizontal parallel creep curve assumption. As a result, the presented creep formulation has the advantage of a single-curve representation for depicting creep strains under a time-varying stress history.Two sets of laboratory experiments were sequentially conducted on cylindrical specimens to obtain the model parameters and to investigate the performance of the presented, two-way parallel creep curve formulation. Both sets of experiments included two cases of axial loads, including both constant and stepwise loads, where the stepwise loads were designed to depict the time-varying stress history. The performance of the presented creep model was compared with the effective modulus method, ageing coefficient method, step-by-step method, and RCM.


## 2. Formulation of Creep

The RCM, originally proposed by Glanville [17], assumes that the change in the rate of creep with time is independent of the age of loading. This means that the creep curves for concrete loaded at different times are assumed to be parallel. The meaning of “parallel” is illustrated in Figure 1(a), which shows creep curves corresponding to loads applied at different times to,t1,t2,t3,…,tn-1. RCM underestimates creep strain because of the inherent limitation imposed by the parallel assumption, as observed from Figure 1(a). Therefore, a parallel creep curve along the horizontal direction is introduced, as shown in Figure 1(b), and combined with the vertical parallel creep curve of RCM. The relationship between the creep function J′(t,to), creep coefficient φ(t,to), and creep compliance function J(t,to) under a constant stress σ(to) applied at time to follows the convention of (1)εt=1+φt,toEctoσcto=1Ecto+J′t,toσcto=Jt,toσcto.Figure 1
Parallel creep curves: (a) vertical parallel creep curve and (b) horizontal parallel creep curve.
(a)
(b)Denoting the vertical parallel creep function asJv′(t,tn-1) and the horizontal parallel creep function as Jh′(t,tn-1), the actual creep function Jact′(t,tn-1) due to load applied at time tn-1 ranges between Jv′(t,tn-1) and Jh′(t,tn-1). Figure 2 shows the relation among the three creep functions due to the load applied at time t=tn-1. The relation among the three creep functions Jv′(t,tn-1)<Jact′(t,tn-1)<Jh′(t,tn-1) was linearly interpolated in the subsequent model formulation by introducing the ageing factor α(t) as follows:(2)Jact′t,tn-1=1-αtn-1Jv′t,tn-1+αtn-1Jh′t,tn-1,where the ageing factor α(tn-1) is used to adjust the decrease in creep over time and ranges between 0 and 1. The remaining part of this paper is focused on formulating the two creep functions Jv′(t,tn-1) and Jh′(t,tn-1) in the form of the initial creep function J′(t,to), defining the creep strain due to the first load applied at time to, which leads to the single creep curve under a time-varying stress history.Figure 2
Concept of the two-way creep model.


### 2.1. One-Way Rate of Creep Formulation

Based on the RCM shown in Figure1(a), the creep strain εcri-1(tn) at time tn due to the loading at time ti-1 is calculated with the vertical parallel creep function Jv′(tn,ti-1) as (3)εcri-1tn=Jv′tn,ti-1Δσti-1,where Δσ(ti-1) is the stress increment at time ti-1. Denoting the creep function due to the first loading at the initial time to as J′(t,to), the vertical parallel creep function Jv′(tn,ti-1) due to the loading at time ti-1 is identical to J′(tn,to)-J′(ti-1,to). Thereby, the creep strain εcr(tn) at time tn due to the loadings at times to,t1,t2,t3,…,tn-1 can be computed by summing the creep strains calculated from (3) as(4)εcrtn=∑i=1nJ′tn,to-J′ti-1,toΔσti-1.The increment of creep strainΔεcr(tn) during the time change from tn-1 to tn is computed by expressing (4) in incremental form:(5)Δεcrtn=ΔJ′t=tn-1σtn-1,where ΔJ′t=tn-1=J′(tn,tn-1)=J′(tn,to)-J′(tn-1,to).The integral form of (4) can be derived by expressing the incremental form of (5) as the rate form: (6)ε˙crt=J˙′tσt;J˙′t=dJ′tdt.Integrating (6) from time to to time t leads to the integral form of the creep strain εcr(t):(7)εcrt=∫totσtdJ′t.In this paper, the creep equation of (5) is termed the one-way RCM because the tangents of the creep curves are identical along a vertical line.


### 2.2. Two-Way Rate of Creep Formulation

A creep curve is introduced under the self-similarity assumption of Figure1(b), where the tangents of the creep curves are identical along the horizontal direction. According to this assumption, the creep strain εcri-1(tn) at time tn due to the loading at time ti-1 is calculated with the horizontal parallel creep function Jh′(tn,ti-1) as(8)εcri-1tn=Jh′tn,ti-1Δσti-1.Representing the horizontal parallel creep function Jh′(tn,ti-1) in the form of J′(t,to), the creep strain εcr(tn) at time tn due to the loadings at times to,t1,t2,t3,…,tn-1 can be computed by summing the creep strains calculated from (8) as(9)εcrtn=∑i=1nJ′tn-ti-1+to,toΔσti-1.The relation of Jh′(tn,ti-1)=J′tn-ti-1+to,to in (8) and (9) is obtained from the geometry of Figure 1(b), where the creep curve due to the loading at time ti-1 is horizontally projected to the initial creep curve J′(t,to).The incremental expressionΔεcr(tn) of (9) can be obtained similarly to the RCM:(10)Δεcrtn=∑i=1nΔJi-1′tnΔσti-1,where ΔJi-1′tn=J′tn-ti-1+to,to-J′tn-1-ti-1+to,to represents the increment of the creep function due to the loading at time ti-1, from time tn-1 to time tn. The full expression of (10) is derived by expressing Δεcr(tn) in (10) as the rate form as follows:(11)ε˙cri-1t=σ˙ti-1J˙i-1′t,where ε˙cri-1(t) denotes the rate of the creep strain at time t due to the load applied at time ti-1. Integrating the rate of creep strain of (11) from time to to time t as well as summing all the creep curves corresponding to the loadings at times to,t1,t2,t3,…,tn-1 leads to the integral form:(12)εcrt=∑i=1nΔσti-1∫tot-ti-1+todJ′t.Since the actual creep curve shown in Figure2 is bounded between the two limits of (4) and (9) (or (7) and (12), if expressed in integral form), the two limit cases of the creep strain are combined by the ageing factor α(t) defined in (2). The resulting incremental expression of creep strain can be obtained by combining (5) and (10) as follows:(13)Δεcrtn=1-αnΔJ′t=tn-1σtn-1+αn∑i=1nΔJi-1′tnΔσti-1,where αn=α(tn). The creep equation of (13) combines the two equality assumptions for the tangents of the creep curves along the vertical and horizontal directions. The equality assumption along the vertical direction corresponds to the case of α=0 and provides the lower limit of the creep strain, whereas that along the horizontal direction corresponds to the case of α=1 and provides the upper limit. The presented creep expression of (13) models the particular creep behavior under multiple time-varying loadings by combining the two ultimate creep formulations. The value of the ageing factor α(t) in (13) can be determined from a set of experiments, including creep tests under constant and stepwise loads. It is noted that the two-way rate of creep formulation of (13) only requires a single creep curve that depicts the creep strain under a constant load applied at time to to calculate the creep strain under a time-varying stress history.


## 2.1. One-Way Rate of Creep Formulation

Based on the RCM shown in Figure1(a), the creep strain εcri-1(tn) at time tn due to the loading at time ti-1 is calculated with the vertical parallel creep function Jv′(tn,ti-1) as (3)εcri-1tn=Jv′tn,ti-1Δσti-1,where Δσ(ti-1) is the stress increment at time ti-1. Denoting the creep function due to the first loading at the initial time to as J′(t,to), the vertical parallel creep function Jv′(tn,ti-1) due to the loading at time ti-1 is identical to J′(tn,to)-J′(ti-1,to). Thereby, the creep strain εcr(tn) at time tn due to the loadings at times to,t1,t2,t3,…,tn-1 can be computed by summing the creep strains calculated from (3) as(4)εcrtn=∑i=1nJ′tn,to-J′ti-1,toΔσti-1.The increment of creep strainΔεcr(tn) during the time change from tn-1 to tn is computed by expressing (4) in incremental form:(5)Δεcrtn=ΔJ′t=tn-1σtn-1,where ΔJ′t=tn-1=J′(tn,tn-1)=J′(tn,to)-J′(tn-1,to).The integral form of (4) can be derived by expressing the incremental form of (5) as the rate form: (6)ε˙crt=J˙′tσt;J˙′t=dJ′tdt.Integrating (6) from time to to time t leads to the integral form of the creep strain εcr(t):(7)εcrt=∫totσtdJ′t.In this paper, the creep equation of (5) is termed the one-way RCM because the tangents of the creep curves are identical along a vertical line.


## 2.2. Two-Way Rate of Creep Formulation

A creep curve is introduced under the self-similarity assumption of Figure1(b), where the tangents of the creep curves are identical along the horizontal direction. According to this assumption, the creep strain εcri-1(tn) at time tn due to the loading at time ti-1 is calculated with the horizontal parallel creep function Jh′(tn,ti-1) as(8)εcri-1tn=Jh′tn,ti-1Δσti-1.Representing the horizontal parallel creep function Jh′(tn,ti-1) in the form of J′(t,to), the creep strain εcr(tn) at time tn due to the loadings at times to,t1,t2,t3,…,tn-1 can be computed by summing the creep strains calculated from (8) as(9)εcrtn=∑i=1nJ′tn-ti-1+to,toΔσti-1.The relation of Jh′(tn,ti-1)=J′tn-ti-1+to,to in (8) and (9) is obtained from the geometry of Figure 1(b), where the creep curve due to the loading at time ti-1 is horizontally projected to the initial creep curve J′(t,to).The incremental expressionΔεcr(tn) of (9) can be obtained similarly to the RCM:(10)Δεcrtn=∑i=1nΔJi-1′tnΔσti-1,where ΔJi-1′tn=J′tn-ti-1+to,to-J′tn-1-ti-1+to,to represents the increment of the creep function due to the loading at time ti-1, from time tn-1 to time tn. The full expression of (10) is derived by expressing Δεcr(tn) in (10) as the rate form as follows:(11)ε˙cri-1t=σ˙ti-1J˙i-1′t,where ε˙cri-1(t) denotes the rate of the creep strain at time t due to the load applied at time ti-1. Integrating the rate of creep strain of (11) from time to to time t as well as summing all the creep curves corresponding to the loadings at times to,t1,t2,t3,…,tn-1 leads to the integral form:(12)εcrt=∑i=1nΔσti-1∫tot-ti-1+todJ′t.Since the actual creep curve shown in Figure2 is bounded between the two limits of (4) and (9) (or (7) and (12), if expressed in integral form), the two limit cases of the creep strain are combined by the ageing factor α(t) defined in (2). The resulting incremental expression of creep strain can be obtained by combining (5) and (10) as follows:(13)Δεcrtn=1-αnΔJ′t=tn-1σtn-1+αn∑i=1nΔJi-1′tnΔσti-1,where αn=α(tn). The creep equation of (13) combines the two equality assumptions for the tangents of the creep curves along the vertical and horizontal directions. The equality assumption along the vertical direction corresponds to the case of α=0 and provides the lower limit of the creep strain, whereas that along the horizontal direction corresponds to the case of α=1 and provides the upper limit. The presented creep expression of (13) models the particular creep behavior under multiple time-varying loadings by combining the two ultimate creep formulations. The value of the ageing factor α(t) in (13) can be determined from a set of experiments, including creep tests under constant and stepwise loads. It is noted that the two-way rate of creep formulation of (13) only requires a single creep curve that depicts the creep strain under a constant load applied at time to to calculate the creep strain under a time-varying stress history.


## 3. Existing Creep Models

### 3.1. Effective Modulus Method

The effective modulus method (EMM), presented by Faber [18], is the oldest and simplest method to calculate the creep strain of concrete subjected to a time-varying stress history. This method assumes that the magnitude of the creep function due to the loading at time ti-1 is the same as the magnitude of the initial creep function J′(tn,to) due to the loading at time to. Therefore, the creep strain εcr(tn) and the corresponding total strain ε(tn) due to the total stress σc(tn) at time tn can be calculated as(14)εcrtn=φtn,toEctoσctn=Jtn,toσctn,(15)εtn=1+φtn,toEctoσctn+εshtn=1Ecetn,toσctn+εshtn,where Ece(to) is the effective modulus. The total stress σc(tn)=∑i=1nΔσ(ti-1) is the sum of the loads applied at different times to,t1,t2,t3,…,tn-1. This method computes the largest value for the creep strain among the five creep models considered in this paper in the case of a time-varying stress history.


### 3.2. Ageing Coefficient Method

An ageing coefficientχ(t,to) was introduced to the effective modulus method of Faber by Bažant [19] to account for the age-dependency of creep at loading of which the earlier the concrete is loaded, the greater the creep strain is (referred to as the age-adjusted effective modulus method (AEMM)). The creep strain εcr(tn) and total strain ε(tn) at time tn can be presented with the age-adjusted effective modulus E-ce(to) as follows:(16)εcrtn=σctoEctoφtn,to+ΔσctnEctoχtn,toφtn,to,(17)εtn=1+φtn,toEctoσcto+1+χtn,toφtn,toEctoΔσctn+εsh=σctoEcetn,to+ΔσctnE-cetn,to+εsh,where E-ce(t,to) is the age-adjusted effective modulus and Δσ(tn) is the stress variation during the time change from to to tn. In (16), the creep model of Faber is recovered when χ(t,to)=1. Since the ageing coefficient χ(t,to) is dependent on the stress level Δσ(tn), its value is obtained from the creep test under the time-varying stress history. In the special case of pure relaxation where only an initial strain is applied, the value of χ(t,to) is computed from [11]. The most typical value of χ(t,to)=0.8 was suggested from the general range of 0.5 to 1 [15].


### 3.3. Step-by-Step Method

The step-by-step method (SSM) superimposes the creep strains due to the loadings at each time increment. In SSM, the continuously varying stress is divided by the specified time intervals, and the creep strain due to a stress increment at a time interval is calculated by the known creep function. Assuming the creep functionJi-1′(tn,ti-1) due to the stress increment Δσ(ti-1) at time interval ti-1, the creep strain εcri-1(tn) at time tn is calculated from (18)εcri-1tn=Ji-1′tn,ti-1Δσti-1.Equation (18) is similar to (3) and (8). The difference between these equations, however, is found in the definition of the creep function, where, in (3) and (8), the initial creep curve due to the first loading at time to defines the subsequent creep curves due to the loads at different times, while (18) requires the same number of creep functions as there are loads. Upon determination of the creep functions, the total creep strain εcr(tn) at time tn can be obtained by superimposing the creep contributions from each time interval as(19)εcrtn=∑i=1nJi-1′tn,ti-1Δσti-1.


## 3.1. Effective Modulus Method

The effective modulus method (EMM), presented by Faber [18], is the oldest and simplest method to calculate the creep strain of concrete subjected to a time-varying stress history. This method assumes that the magnitude of the creep function due to the loading at time ti-1 is the same as the magnitude of the initial creep function J′(tn,to) due to the loading at time to. Therefore, the creep strain εcr(tn) and the corresponding total strain ε(tn) due to the total stress σc(tn) at time tn can be calculated as(14)εcrtn=φtn,toEctoσctn=Jtn,toσctn,(15)εtn=1+φtn,toEctoσctn+εshtn=1Ecetn,toσctn+εshtn,where Ece(to) is the effective modulus. The total stress σc(tn)=∑i=1nΔσ(ti-1) is the sum of the loads applied at different times to,t1,t2,t3,…,tn-1. This method computes the largest value for the creep strain among the five creep models considered in this paper in the case of a time-varying stress history.


## 3.2. Ageing Coefficient Method

An ageing coefficientχ(t,to) was introduced to the effective modulus method of Faber by Bažant [19] to account for the age-dependency of creep at loading of which the earlier the concrete is loaded, the greater the creep strain is (referred to as the age-adjusted effective modulus method (AEMM)). The creep strain εcr(tn) and total strain ε(tn) at time tn can be presented with the age-adjusted effective modulus E-ce(to) as follows:(16)εcrtn=σctoEctoφtn,to+ΔσctnEctoχtn,toφtn,to,(17)εtn=1+φtn,toEctoσcto+1+χtn,toφtn,toEctoΔσctn+εsh=σctoEcetn,to+ΔσctnE-cetn,to+εsh,where E-ce(t,to) is the age-adjusted effective modulus and Δσ(tn) is the stress variation during the time change from to to tn. In (16), the creep model of Faber is recovered when χ(t,to)=1. Since the ageing coefficient χ(t,to) is dependent on the stress level Δσ(tn), its value is obtained from the creep test under the time-varying stress history. In the special case of pure relaxation where only an initial strain is applied, the value of χ(t,to) is computed from [11]. The most typical value of χ(t,to)=0.8 was suggested from the general range of 0.5 to 1 [15].


## 3.3. Step-by-Step Method

The step-by-step method (SSM) superimposes the creep strains due to the loadings at each time increment. In SSM, the continuously varying stress is divided by the specified time intervals, and the creep strain due to a stress increment at a time interval is calculated by the known creep function. Assuming the creep functionJi-1′(tn,ti-1) due to the stress increment Δσ(ti-1) at time interval ti-1, the creep strain εcri-1(tn) at time tn is calculated from (18)εcri-1tn=Ji-1′tn,ti-1Δσti-1.Equation (18) is similar to (3) and (8). The difference between these equations, however, is found in the definition of the creep function, where, in (3) and (8), the initial creep curve due to the first loading at time to defines the subsequent creep curves due to the loads at different times, while (18) requires the same number of creep functions as there are loads. Upon determination of the creep functions, the total creep strain εcr(tn) at time tn can be obtained by superimposing the creep contributions from each time interval as(19)εcrtn=∑i=1nJi-1′tn,ti-1Δσti-1.


## 4. Experiments and Numerical Predictions

Two sets of time-dependent laboratory experiments of types A and B were sequentially conducted on cylindrical concrete specimens to determine the value of the ageing factorα(t) defined in (2) and (13). Two load cases, including constant and stepwise loads, were taken into account to find the value of the ageing factor in type A test. The same load cases were employed to type B test to validate the value of the ageing factor obtained.


### 4.1. Cylindrical Concrete Specimens

All cylindrical concrete specimens were cast with 150 mm diameter and 300 mm height. For each type of test, thirty-five cylindrical concrete specimens were cast, of which twelve specimens were used for the creep tests in uniaxial compression, three specimens were used to measure the shrinkage, and twenty specimens were used to measure the development of the elastic modulus. Table1 lists the concrete mixture proportions for the cylindrical specimens for the two test types where different slump values were taken to provide more experimental evidence for the performance of the presented creep formulation. All the specimens were covered with nonwoven fabric so that they were not exposed to the air and were fully moisturized by sprinkling water until the formworks were removed. The specimens were removed from the molds 24 hours after casting and then cured in water. At the age of 10 days (7 days for type B test), the initial value of the elastic modulus was measured from a direct compression test. Test data, including creep, shrinkage, and elastic modulus, were obtained from strain gages attached at the midheight of the specimen at equal angles of 90°. Six creep-testing apparatuses were used for the creep tests where two specimens were set in series in each apparatus. Two of the testing apparatuses were used to measure the creep strain under a constant load. The other four were used to measure creep strains under the two cases of stepwise loads. Figures 3(a) and 3(b) show the load histories of the two stepwise load cases applied to type A and type B tests, respectively.Table 1
Mixture proportions for the two test types of concrete.
Test type
Max. size of aggr. (mm)
W/C ratio (%)
Slump (mm)
Unit weight,  (N/m3)
Water
Cement
Fine aggr.
Coarse aggr.
Type A
20
54
110
1,880
3,480
7,360
10,200
Type B
20
57
200
2,120
3,720
9,450
7,610Figure 3
Stepwise loads applied to (a) type A test and (b) type B test.
(a)
(b)Hydraulic pressure was applied by a hydraulic jack. The pressure level was monitored by a pressure meter attached to the jack and a load cell mounted at the top of the specimen throughout the testing period (see Figure4). Since the leak of hydraulic pressure causes undesirable specimen deformation, the locking bolts in Figure 4 were locked immediately after the hydraulic pressure was applied. However, the locking of the bolts did not interfere with the release of the spring that was induced, owing to the progressive deformation of the specimen. The loss of axial pressure due to the release of the spring was supplemented every two days to maintain the pressure level, as designated in Figures 3(a) and 3(b).Figure 4
Schematic drawing of creep test setup.


### 4.2. Basic Age-Dependent Properties

All experiments were conducted in a controlled room with a constant temperature of22±1°C and relative humidity of 60±2% RH for 120 days. Gage readings from the specimens were automatically stored every minute in a data acquisition system with 81 data channels. All stored data were extracted from the acquisition system every 7 days to avoid an accidental loss of data. For the creep test under constant load, an initial pressure of 5 MPa was applied at the age of 10 days for type A test (7 days for type B test). For the creep test under stepwise loads, initial pressure values of 5 and 6 MPa were applied, respectively, for the two stepwise load cases of type A test, and the same initial pressure of 5 MPa was applied for the two load cases of type B test. The final axial pressure was 11 MPa at the end of the stepwise loadings for both type A and type B tests (Figures 3(a) and 3(b)).When concrete is loaded with a time-varying stress history, the present single-curve formulation requires a creep curve that describes the creep strain under a constant load. The creep curve was calibrated from the constant load case for both type A and type B tests. The resulting three equations of creep strain to defineJ′(t,to), shrinkage strain, and development of elastic modulus are listed in Table 2. The values of the elastic modulus were measured at the ages of 7 (10 in type A test), 14, 28, and 56 days. The parameter εshu=7.8×10-4 in Table 2 defines the shrinkage strain at t=∞ according to [5, 9], and fc′ is the uniaxial compressive strength at the age of 28 days. Figures 5(a) and 5(b) show the total, creep, and shrinkage strains, as well as development of the elastic modulus for type A test, respectively, where the creep strain was obtained by subtracting the immediate elastic strain and shrinkage strain from the total strain measurements.Table 2
Three empirical equations for the two types of tests.
Test type
Parameters
a
1
b
1
a
2
b
2
a
3
b
3
t
o (days)
E
c
t
=
t
o (MPa)
E
c
t
=
28 (MPa)
f
c
′ (MPa)
A
9.5
1.7
51
1.45
4
0.85
10
29,400
32,800
30
B
11.4
5.1
48.4
0.73
2.6
0.9
7
22,700
25,300
28
J
′
t
,
t
o
=
t
-
t
o
0.6
a
1
+
t
-
t
o
0.6
b
1
E
c
t
=
t
o,  εsht,to=t-toa2+b2t-toεshu,  Ect=ta3+b3tEct=28daysFigure 5
Basic age-dependent properties of type A test: (a) total, creep, and shrinkage strains, (b) development of elastic modulus.
(a)
(b)


### 4.3. Determination of the Ageing Factor Value

The value of the ageing factorα(t) defining the creep formulation in (2) and (13) was determined from the creep tests under the stepwise loads depicted in Figure 3(a). Figures 6(a) and 6(b) compare the measured and predicted total strains where the measured total strain for each load case was obtained by averaging the total strains measured from four specimens. The numerical predictions were conducted by increasing the value of α from 0 to 1 in 0.25 increments, where α=α(t) was assumed to be independent of time. The creep strains were extracted by subtracting the immediate elastic strain and shrinkage strain from the total strains and compared in Figures 7(a) and 7(b) for the two load cases 1 and 2, respectively. It is observed from Figures 6 and 7 that the creep strain is significantly underestimated in the vertical parallel creep curve method corresponding to α=0, whereas it is overestimated in the horizontal parallel creep curve assumption corresponding to α=1. The value of the ageing factor α=0.75, obtained from the comparisons, indicates that, in the current creep formulation, the creep strain is predominantly predicted by the horizontal parallel creep curve assumption, rather than the vertical creep curve assumption for concrete loaded at a relatively early age.Figure 6
Predicted and measured total strains for the two stepwise load cases of type A test: (a) load case 1 and (b) load case 2.
(a)
(b)Figure 7
Predicted and measured creep strains for the two stepwise load cases of type A test: (a) load case 1 and (b) load case 2.
(a)
(b)Strain measurements of the two stepwise load cases for type B test in Figure3(b) were predicted by the presented creep model to verify whether the value of the ageing factor α=0.75, obtained from type A test, is acceptable. The value of the ageing factor was increased in the prediction from α=0 to α=1 with 0.25 increments. Three basic equations used in the analysis are listed in Table 2 where the creep function was obtained from the creep test under a constant pressure of 5 MPa. Total strains measured for the two stepwise load cases 1 and 2 are compared to the predicted strains in Figures 8(a) and 8(b), respectively. The creep strains are compared in Figures 9(a) and 9(b) for the two load cases 1 and 2, respectively. It is noted from Figures 8 and 9 that the measured strains match well with the predicted strains for the two stepwise load cases when α=0.75, and the development of creep is similar to type A test.Figure 8
Predicted and measured total strains for the two stepwise load cases of type B test: (a) load case 1 and (b) load case 2.
(a)
(b)Figure 9
Predicted and measured creep strains for the two stepwise load cases of type B test: (a) load case 1 and (b) load case 2.
(a)
(b)


### 4.4. Prediction by Existing Creep Models

Strain measurements for the two stepwise load cases of type A and type B tests were predicted and compared using five creep models, including EMM, AEMM, RCM, SSM, and the method presented in this study. A constant value of the ageing coefficientχ(t,to)=0.8 was used in the AEMM case. The values of the elastic moduli Ec(to=10days) and Ec(to=7days) for type A and type B tests in Table 2, respectively, defined the elastic modulus in (15) and (17).In the SSM case, the continuous creep functionJ′(t,to) for the two stepwise load cases in Figure 3 were obtained based on ACI 209R-92 [2]. According to ACI 209R-92, the creep function J′(t,ti-1) due to the loading at time to is computed from(20)J′t,ti-1=t-ti-1ψa+t-ti-1ψϕuti-1Ecti-1,where ϕu(ti-1) is the ultimate creep coefficient defined with six correction factors in the case of nonstandard conditions. The coefficient ϕu(ti-1) accounts for the correction factor for the time at loading under moist curing conditions as follows:(21)ϕuti-1=bγc,ti-1,where γc,ti-1=1.25ti-1-0.118 is the correction factor due to the time at loading. Substituting (21) into (20), the creep function J′(t,ti-1) at loading time ti-1 is obtained as follows:(22)J′t,ti-1=t-ti-1ψa+t-ti-1ψ1.25bti-1-0.118Eti-1.The parameter values a, ψ, and b in (22) were determined through regression analyses for the creep tests of types A and B at loading times ti-1=10 days and ti-1=7 days, respectively (Table 2). The final forms of the creep functions, according to ACI 209R-92, were obtained as follows:(23)J′t,ti-1=2.23ti-1-0.118t-ti-10.6Eti-19.5+t-ti-10.6type  A,J′t,ti-1=6.42ti-1-0.118t-ti-10.6Eti-111.4+t-ti-10.6type  B.The types of variables and the model format were compared among the five creep models to identify the model performance depending on those modeling factors. Three modeling factors were verified in the comparisons: whether the creep function accounts for the time-varying stress history, whether the development of the elastic modulus is taken into account in the model formulation, and whether the formulation format is based on incremental time or total time. Regarding the time-varying stress history, three models, AEMM, SSM, and the presented model, account for the time-varying stress history in their creep functions. Regarding the development of elastic modulus, two models, SSM and the presented model, take into account the change of elastic modulus at the age of loadings while the other three models, AEMM, EMM, and RCM, use the elastic modulus evaluated at the age of initial loading. Regarding the formulation format, two models, EMM and AEMM, are based on the total format and three models, RCM, SSM, and the presented model, are based on the incremental format.Figures10 and 11 compare the total strains predicted by the five creep models for the two stepwise load histories of type A and type B tests, respectively. Figures 10 and 11 show the largest strain for the EMM case among the five creep models. This is because the model uses the values of the modulus of elasticity E(to) and creep function J′(t,to) evaluated at the initial times to=10 days and 7 days for type A and type B tests, respectively, through the entire time period. RCM computes the lowest value of strain, as observed from Figures 10 and 11. The strains predicted by the other three creep models are bounded between EMM and RCM. The strain predicted by AEMM is lower than that predicted by EMM but higher than those by the other three methods. However, this is not an issue in the current numerical prediction because the ageing coefficient value χ(t,to) was set at a constant value of 0.8 through the entire time period. This means that, in this particular creep test program, the value of χ(t,t0) may be less than 0.8. In the cases of EMM, AEMM, and RCM, the value of the elastic modulus was taken as the value evaluated at the initial times of to=10 days and 7 days for type A and type B tests, respectively, while, in the cases of SSM and the presented creep model, a variable value of elastic modulus was used in the analyses. The difference between the two cases is identified in Figures 10 and 11, where a gradual increase of time-dependent strain is observed in the latter two methods, while an acute increase of time-dependent strain is observed in the former three methods when the stepwise load was applied.Figure 10
Comparison of the total strain predicted by the five creep models (type A test): (a) load case 1 and (b) load case 2.
(a)
(b)Figure 11
Comparison of the total strains predicted by the five creep models (type B test): (a) load case 1 and (b) load case 2.
(a)
(b)To investigate the long-term performances of the five models, an axial load of 2 MPa was numerically applied to both of the specimens in type A and type B tests at the ages of 90 days and 103 days, respectively, when the experiments were completed. Strains were computed for 500 days for the five creep models and compared in Figures12 and 13 for type A and type B tests, respectively. In the case of type A test in Figure 12, the strain predicted by the presented model is approximately the same as that predicted by AEMM, and the strain predicted by SSM is lower than both of those. However, the case of type B test in Figure 13 shows that the presented model predicts a strain larger than that of AEMM and SSM. It is noted from the long-term analyses that the AEMM, SSM, and the presented creep models predict creep strains between those of the EMM and RCM models and agree well with each other in the long-term analysis, as well as the short-term analysis.Figure 12
Extended predictions and comparison of the total strain by the five creep models (type A test): (a) load case 1 and (b) load case 2.
(a)
(b)Figure 13
Extended predictions and comparison of the total strain by the five creep models (type B test): (a) load case 1 and (b) load case 2.
(a)
(b)


## 4.1. Cylindrical Concrete Specimens

All cylindrical concrete specimens were cast with 150 mm diameter and 300 mm height. For each type of test, thirty-five cylindrical concrete specimens were cast, of which twelve specimens were used for the creep tests in uniaxial compression, three specimens were used to measure the shrinkage, and twenty specimens were used to measure the development of the elastic modulus. Table1 lists the concrete mixture proportions for the cylindrical specimens for the two test types where different slump values were taken to provide more experimental evidence for the performance of the presented creep formulation. All the specimens were covered with nonwoven fabric so that they were not exposed to the air and were fully moisturized by sprinkling water until the formworks were removed. The specimens were removed from the molds 24 hours after casting and then cured in water. At the age of 10 days (7 days for type B test), the initial value of the elastic modulus was measured from a direct compression test. Test data, including creep, shrinkage, and elastic modulus, were obtained from strain gages attached at the midheight of the specimen at equal angles of 90°. Six creep-testing apparatuses were used for the creep tests where two specimens were set in series in each apparatus. Two of the testing apparatuses were used to measure the creep strain under a constant load. The other four were used to measure creep strains under the two cases of stepwise loads. Figures 3(a) and 3(b) show the load histories of the two stepwise load cases applied to type A and type B tests, respectively.Table 1
Mixture proportions for the two test types of concrete.
Test type
Max. size of aggr. (mm)
W/C ratio (%)
Slump (mm)
Unit weight,  (N/m3)
Water
Cement
Fine aggr.
Coarse aggr.
Type A
20
54
110
1,880
3,480
7,360
10,200
Type B
20
57
200
2,120
3,720
9,450
7,610Figure 3
Stepwise loads applied to (a) type A test and (b) type B test.
(a)
(b)Hydraulic pressure was applied by a hydraulic jack. The pressure level was monitored by a pressure meter attached to the jack and a load cell mounted at the top of the specimen throughout the testing period (see Figure4). Since the leak of hydraulic pressure causes undesirable specimen deformation, the locking bolts in Figure 4 were locked immediately after the hydraulic pressure was applied. However, the locking of the bolts did not interfere with the release of the spring that was induced, owing to the progressive deformation of the specimen. The loss of axial pressure due to the release of the spring was supplemented every two days to maintain the pressure level, as designated in Figures 3(a) and 3(b).Figure 4
Schematic drawing of creep test setup.


## 4.2. Basic Age-Dependent Properties

All experiments were conducted in a controlled room with a constant temperature of22±1°C and relative humidity of 60±2% RH for 120 days. Gage readings from the specimens were automatically stored every minute in a data acquisition system with 81 data channels. All stored data were extracted from the acquisition system every 7 days to avoid an accidental loss of data. For the creep test under constant load, an initial pressure of 5 MPa was applied at the age of 10 days for type A test (7 days for type B test). For the creep test under stepwise loads, initial pressure values of 5 and 6 MPa were applied, respectively, for the two stepwise load cases of type A test, and the same initial pressure of 5 MPa was applied for the two load cases of type B test. The final axial pressure was 11 MPa at the end of the stepwise loadings for both type A and type B tests (Figures 3(a) and 3(b)).When concrete is loaded with a time-varying stress history, the present single-curve formulation requires a creep curve that describes the creep strain under a constant load. The creep curve was calibrated from the constant load case for both type A and type B tests. The resulting three equations of creep strain to defineJ′(t,to), shrinkage strain, and development of elastic modulus are listed in Table 2. The values of the elastic modulus were measured at the ages of 7 (10 in type A test), 14, 28, and 56 days. The parameter εshu=7.8×10-4 in Table 2 defines the shrinkage strain at t=∞ according to [5, 9], and fc′ is the uniaxial compressive strength at the age of 28 days. Figures 5(a) and 5(b) show the total, creep, and shrinkage strains, as well as development of the elastic modulus for type A test, respectively, where the creep strain was obtained by subtracting the immediate elastic strain and shrinkage strain from the total strain measurements.Table 2
Three empirical equations for the two types of tests.
Test type
Parameters
a
1
b
1
a
2
b
2
a
3
b
3
t
o (days)
E
c
t
=
t
o (MPa)
E
c
t
=
28 (MPa)
f
c
′ (MPa)
A
9.5
1.7
51
1.45
4
0.85
10
29,400
32,800
30
B
11.4
5.1
48.4
0.73
2.6
0.9
7
22,700
25,300
28
J
′
t
,
t
o
=
t
-
t
o
0.6
a
1
+
t
-
t
o
0.6
b
1
E
c
t
=
t
o,  εsht,to=t-toa2+b2t-toεshu,  Ect=ta3+b3tEct=28daysFigure 5
Basic age-dependent properties of type A test: (a) total, creep, and shrinkage strains, (b) development of elastic modulus.
(a)
(b)


## 4.3. Determination of the Ageing Factor Value

The value of the ageing factorα(t) defining the creep formulation in (2) and (13) was determined from the creep tests under the stepwise loads depicted in Figure 3(a). Figures 6(a) and 6(b) compare the measured and predicted total strains where the measured total strain for each load case was obtained by averaging the total strains measured from four specimens. The numerical predictions were conducted by increasing the value of α from 0 to 1 in 0.25 increments, where α=α(t) was assumed to be independent of time. The creep strains were extracted by subtracting the immediate elastic strain and shrinkage strain from the total strains and compared in Figures 7(a) and 7(b) for the two load cases 1 and 2, respectively. It is observed from Figures 6 and 7 that the creep strain is significantly underestimated in the vertical parallel creep curve method corresponding to α=0, whereas it is overestimated in the horizontal parallel creep curve assumption corresponding to α=1. The value of the ageing factor α=0.75, obtained from the comparisons, indicates that, in the current creep formulation, the creep strain is predominantly predicted by the horizontal parallel creep curve assumption, rather than the vertical creep curve assumption for concrete loaded at a relatively early age.Figure 6
Predicted and measured total strains for the two stepwise load cases of type A test: (a) load case 1 and (b) load case 2.
(a)
(b)Figure 7
Predicted and measured creep strains for the two stepwise load cases of type A test: (a) load case 1 and (b) load case 2.
(a)
(b)Strain measurements of the two stepwise load cases for type B test in Figure3(b) were predicted by the presented creep model to verify whether the value of the ageing factor α=0.75, obtained from type A test, is acceptable. The value of the ageing factor was increased in the prediction from α=0 to α=1 with 0.25 increments. Three basic equations used in the analysis are listed in Table 2 where the creep function was obtained from the creep test under a constant pressure of 5 MPa. Total strains measured for the two stepwise load cases 1 and 2 are compared to the predicted strains in Figures 8(a) and 8(b), respectively. The creep strains are compared in Figures 9(a) and 9(b) for the two load cases 1 and 2, respectively. It is noted from Figures 8 and 9 that the measured strains match well with the predicted strains for the two stepwise load cases when α=0.75, and the development of creep is similar to type A test.Figure 8
Predicted and measured total strains for the two stepwise load cases of type B test: (a) load case 1 and (b) load case 2.
(a)
(b)Figure 9
Predicted and measured creep strains for the two stepwise load cases of type B test: (a) load case 1 and (b) load case 2.
(a)
(b)


## 4.4. Prediction by Existing Creep Models

Strain measurements for the two stepwise load cases of type A and type B tests were predicted and compared using five creep models, including EMM, AEMM, RCM, SSM, and the method presented in this study. A constant value of the ageing coefficientχ(t,to)=0.8 was used in the AEMM case. The values of the elastic moduli Ec(to=10days) and Ec(to=7days) for type A and type B tests in Table 2, respectively, defined the elastic modulus in (15) and (17).In the SSM case, the continuous creep functionJ′(t,to) for the two stepwise load cases in Figure 3 were obtained based on ACI 209R-92 [2]. According to ACI 209R-92, the creep function J′(t,ti-1) due to the loading at time to is computed from(20)J′t,ti-1=t-ti-1ψa+t-ti-1ψϕuti-1Ecti-1,where ϕu(ti-1) is the ultimate creep coefficient defined with six correction factors in the case of nonstandard conditions. The coefficient ϕu(ti-1) accounts for the correction factor for the time at loading under moist curing conditions as follows:(21)ϕuti-1=bγc,ti-1,where γc,ti-1=1.25ti-1-0.118 is the correction factor due to the time at loading. Substituting (21) into (20), the creep function J′(t,ti-1) at loading time ti-1 is obtained as follows:(22)J′t,ti-1=t-ti-1ψa+t-ti-1ψ1.25bti-1-0.118Eti-1.The parameter values a, ψ, and b in (22) were determined through regression analyses for the creep tests of types A and B at loading times ti-1=10 days and ti-1=7 days, respectively (Table 2). The final forms of the creep functions, according to ACI 209R-92, were obtained as follows:(23)J′t,ti-1=2.23ti-1-0.118t-ti-10.6Eti-19.5+t-ti-10.6type  A,J′t,ti-1=6.42ti-1-0.118t-ti-10.6Eti-111.4+t-ti-10.6type  B.The types of variables and the model format were compared among the five creep models to identify the model performance depending on those modeling factors. Three modeling factors were verified in the comparisons: whether the creep function accounts for the time-varying stress history, whether the development of the elastic modulus is taken into account in the model formulation, and whether the formulation format is based on incremental time or total time. Regarding the time-varying stress history, three models, AEMM, SSM, and the presented model, account for the time-varying stress history in their creep functions. Regarding the development of elastic modulus, two models, SSM and the presented model, take into account the change of elastic modulus at the age of loadings while the other three models, AEMM, EMM, and RCM, use the elastic modulus evaluated at the age of initial loading. Regarding the formulation format, two models, EMM and AEMM, are based on the total format and three models, RCM, SSM, and the presented model, are based on the incremental format.Figures10 and 11 compare the total strains predicted by the five creep models for the two stepwise load histories of type A and type B tests, respectively. Figures 10 and 11 show the largest strain for the EMM case among the five creep models. This is because the model uses the values of the modulus of elasticity E(to) and creep function J′(t,to) evaluated at the initial times to=10 days and 7 days for type A and type B tests, respectively, through the entire time period. RCM computes the lowest value of strain, as observed from Figures 10 and 11. The strains predicted by the other three creep models are bounded between EMM and RCM. The strain predicted by AEMM is lower than that predicted by EMM but higher than those by the other three methods. However, this is not an issue in the current numerical prediction because the ageing coefficient value χ(t,to) was set at a constant value of 0.8 through the entire time period. This means that, in this particular creep test program, the value of χ(t,t0) may be less than 0.8. In the cases of EMM, AEMM, and RCM, the value of the elastic modulus was taken as the value evaluated at the initial times of to=10 days and 7 days for type A and type B tests, respectively, while, in the cases of SSM and the presented creep model, a variable value of elastic modulus was used in the analyses. The difference between the two cases is identified in Figures 10 and 11, where a gradual increase of time-dependent strain is observed in the latter two methods, while an acute increase of time-dependent strain is observed in the former three methods when the stepwise load was applied.Figure 10
Comparison of the total strain predicted by the five creep models (type A test): (a) load case 1 and (b) load case 2.
(a)
(b)Figure 11
Comparison of the total strains predicted by the five creep models (type B test): (a) load case 1 and (b) load case 2.
(a)
(b)To investigate the long-term performances of the five models, an axial load of 2 MPa was numerically applied to both of the specimens in type A and type B tests at the ages of 90 days and 103 days, respectively, when the experiments were completed. Strains were computed for 500 days for the five creep models and compared in Figures12 and 13 for type A and type B tests, respectively. In the case of type A test in Figure 12, the strain predicted by the presented model is approximately the same as that predicted by AEMM, and the strain predicted by SSM is lower than both of those. However, the case of type B test in Figure 13 shows that the presented model predicts a strain larger than that of AEMM and SSM. It is noted from the long-term analyses that the AEMM, SSM, and the presented creep models predict creep strains between those of the EMM and RCM models and agree well with each other in the long-term analysis, as well as the short-term analysis.Figure 12
Extended predictions and comparison of the total strain by the five creep models (type A test): (a) load case 1 and (b) load case 2.
(a)
(b)Figure 13
Extended predictions and comparison of the total strain by the five creep models (type B test): (a) load case 1 and (b) load case 2.
(a)
(b)


## 5. Conclusions

A single-curve creep model was presented to depict the creep behavior of concrete subjected to time-varying stress history. The formulation was accomplished by introducing a horizontal parallel creep curve assumption and combining with the vertical parallel creep curve concept of RCM. An ageing factor, defining the effect of concrete age on loading, was introduced to combine the two parallel creep curve concepts. Two sets of laboratory experiments on cylindrical specimens cast with different mixture proportions were conducted to determine the value of the ageing factor. The performance of the presented creep formulation was investigated by comparing the measured strains to the predicted strains and then comparing to the existing four creep models. The following conclusions were drawn:(1)
The presented creep model uses a single creep curve to calculate creep strain under a time-varying stress history. This representation of creep simplifies the age-dependent formulation of the long-term behavior of concrete structures and requires no memory of stress history.(2)
The presented formulation combines the two creep curve concepts that define the upper and lower bounds of creep strain, thereby providing a reliable level of accuracy for the predicted creep strain.(3)
The presented creep model was formulated in an incremental format that allowed accounting for the time-dependent development of the elastic modulus. Furthermore, the incremental form of the presented creep model paves the way for being associated with the incremental form of the time-dependent finite element formulation.(4)
A single value of 0.75 was obtained for the ageing factor through the creep tests under stepwise loads. The value of the ageing factor indicated that the creep strain is predominantly predicted by the horizontal parallel creep curve assumption, rather than the vertical creep curve assumption, for concrete loaded at a relatively early age.(5)
The predicted strains in the presented creep model were compared to the measured strains and those of the existing four creep models. This indicated that the predicted strains of ageing coefficient method, step-by-step method, and the presented model agreed well with each other and ranged between those of EMM and RCM.



---

*Source: 9370514-2016-02-14.xml*