Next Article in Journal
Comparison of Screw-In Forces during Movement of Endodontic Files with Different Geometries, Alloys, and Kinetics
Next Article in Special Issue
Characterization of Titanium Nanotube Reinforced Cementitious Composites: Mechanical Properties, Microstructure, and Hydration
Previous Article in Journal / Special Issue
Evaluation of Carbon Nanotube Incorporation in Cementitious Composite Materials
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Creep Damage Constitutive Model of Concrete Based on Fractional Calculus Theory

1
Key Laboratory of Ministry of Education for Geomechanics and Embankment Engineering, Hohai University, Nanjing 210098, China
2
Jiangsu Research Center for Geotechnical Engineering Technology, Hohai University, Nanjing 210098, China
3
Henan Province Key Laboratory of Rock and Soil Mechanics and Structural Engineering, North China University of Water Resources and Electric Power, Zhengzhou 450045, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(9), 1505; https://doi.org/10.3390/ma12091505
Submission received: 27 March 2019 / Revised: 3 May 2019 / Accepted: 6 May 2019 / Published: 8 May 2019
(This article belongs to the Special Issue Concrete and Construction Materials)

Abstract

:
Concrete creep has become one of the major problems that threatens concrete structural development and construction. However, a reasonable and accurate calculation model for numerical analysis is the key to control and solve the creep deformation of concrete. To better describe the concrete nonlinear creep damage evolution rule, the visco-elasticity Plasticity Rheological Theory, Riemann Liouville Theory and Combined Model Theory are quoted, and the Able dashpot is used to reconstruct fractional-order soft-body composite elements to propose the expression of the stress-strain relationship of the elastomer, visco-elasticity plasticity body, and Viscoplasticity body, considering the evolution of the concrete compression damage process. A nonlinear creep damage constitutive model of concrete, based on fractional calculus theory, is conducted, and the parameters of the specific calculation method of the model are given. The influence of stress level σ, fractional order n and material parameter α on the concrete creep process is determined by a sensitivity analysis of the model parameters. The creep process and deformation amount of concrete in practical engineering can be effectively controlled by the results of the proposed sensitivity analysis. The research results can be used to provide guidance and reference for the safe construction of concrete engineering in actual practice.

1. Introduction

Concrete creep refers to the phenomenon in which concrete deformation increases with time under long-term external load [1]. Creep development increases stress loss in the concrete structure. It also redistributes the internal force of static and statically indeterminate structures, resulting in excessive structural deformation, a significant reduction in overall strength, and even a loss of bearing capacity. With the current wide range, high quality, and general use of concrete materials, concrete creep has become one of the major problems that threatens concrete structural development and the safety of concrete construction and long-term stability [2,3].
Creep strain has a linear relationship with applied stress in the case of low stress levels (less than 30% compressive strength). However, creep presents an unstable phenomenon as time and stress increase in the case of high stress levels (between 30% and 80% compressive strength); that is, it presents a nonlinear relationship, which is found in concrete creep theories and experimental research [4,5,6,7]. At present, many linear creep models have been proposed to predict concrete creep deformation, such as the CEB-FIP model [8], the EC2 model [9], the ACI209 model [10], the GL2000 model [11], etc. The B3 model is constructed by using the Humidity Diffusion Theory and Consolidation Theory, and the various material parameters are determined through the specific concrete tests of Bazant [12,13]. The B4 model is proposed based on the B3 model, considering the effects of ambient temperature and self-shrinkage described by Mija H. Hubler [14,15]. This calculation analysis method and the inverse analysis method were introduced to conduct a concrete creep-fatigue constitutive model, based on the application of continuous damage mechanics that was proposed by Farnam Ghasemzadeh [16] and Vitaliy M. Kindrachuk et al. [17]. However, concrete nonlinear creep is very complex because it can be regarded as a composite material. The basic Maxwell model is referenced to propose the linearized mathematical modeling of the geometric nonlinearity theory proposed by Grenacher and Thurlimann [18]. But, the accuracy of concrete structure creep deformation cannot be accurately described by the above approximate processing model.
The nonlinear creep process of reinforced concrete columns was simulated considering the high stress conditions by Bažant and Tsubaki [19]. A simple numerical model was conducted to simplify the concrete nonlinear creep calculation process of Kwak and Kim [20]. Geometrical nonlinearity and concrete fracture development were described using a layered method, and slender reinforced concrete column creep deformation was predicted in the model presented by Kwak and Kim. And, a numerical model for analysis of reinforced and prestressed concrete plates and shells, including creep, shrinkage, and aging of concrete, already developed by Radnić and Matešan [21], can and the model can be applied for all levels of concrete stresses, while its use for ultimate stress levels was still not fully tested. For most multi-coefficient nonlinear creep models, the range of parameters can be further extended, and the time effectiveness and applicability of nonlinear creep have yet to be further investigated and verified for concrete linear and nonlinear creep.
Some strain in creep is related to the corresponding instantaneous stress state, and also to the previous stress state, which is illustrated by a large number of concrete creep phenomena. The fractional derivative, based on time parameters, is an arbitrary order time differential operator, which can fully reflect the course of the function change over time [22,23]. Fractional Calculus Theory is widely used in the field of materials [24], and the Fractional Derivative was introduced by Papoulia, K.D. [25] to establish a Viscoelastic Rheological Model that reveal the viscoelastic behavioral and material mechanical response. Viscoelastic behavior and response are consistent with the experimental results. The visco-elasticity plasticity model, based on the fractional derivative, was established to improve the viscoelastic-plastic response characteristics of materials by Suzuki, J.L. et al. [26]. A high temperature concrete creep nonlinear thermal viscoelastic creep model based on the fractional derivative was proposed to better describe nonlinear creep under conditions of high stress and high temperature by Bouras, Y. et al. [27]. The concrete fractional order genetic aging model was conducted by Beltempo, Angela et al. [28] in the variable-order fractional order calculus framework. The semi-analytic expression of material parameters was determined by using fractional order calculus forms.
To better describe the concrete nonlinear creep damage evolution rule, Fractional Calculus Theory and Rheological Model Theory have been quoted to propose a non-linear creep damage constitutive model of concrete using Able dashpot to replace Newton dashpot and the concrete damage evolution rule. The rationality and applicability of the conducted constitutive work was verified. Finally, the calculation method and sensitivity analysis of the conducted model parameters provide a reference for safe construction and the long-term stable operation of concrete engineering.

2. Construction of Concrete Nonlinear Creep Damage Constitutive Model

2.1. Riemann Liouville Fractional Order Software Components

There are many fractional order calculus definitions, among which the creep characteristics of concrete are defined and described by Riemann Liouville theory [29,30]. This is how the Riemann Liouville fractional calculus of function f t is defined: f ( · ) is set on the (0, + ) continuum and can be accumulated on any subinterval of [0, + ], for t > 0, Re(n) > 0.
D n f ( t ) = 0 t ( t τ ) n 1 Γ ( n ) f ( τ ) d τ
The Riemann Liouville fractional order differential form is defined as.
D n f ( t ) = d a d t a 0 t f ( t τ ) a n 1 Γ ( a n ) f ( τ ) d τ
where Γ ( n ) = 0 t n 1 e t d t is the gamma function, n (n > 0) is the fractional order, a ( a 1 n a ) is a positive integer larger than n, and d is a differential operator. The ideal solid stress-strain relationship follows Hooke’s law.
σ t = E ε t
The ideal fluid stress-strain relationship follows Newton’s viscous law.
σ t = η d ε ( t ) d t
where σ ( t ) is a function of stress versus time, ε ( t ) is a function of strain versus time, E is the elastic modulus, and t is time.
Ideal solids and ideal fluids are idealized models, which can only be used for approximate calculations in practical engineering. Concrete is a material between an ideal solid and an ideal fluid, and a software element [31] can be used to describe an intermediate material between a pure elastomer and a Newtonian fluid. Therefore, the concrete rheological properties can be described by combining a model with the Able dashpot instead of the Newton dashpot. The Newton dashpot and the Able dashpot are shown in Figure 1.
The constitutive equation based on the fractional derivative [32] is established as.
σ t = η n d n ε ( t ) d t n , 0 n 1
where the stress σ remains constant; that is, σ = c o n s t . This software element can describe the creep deformation. Fractional order integrals are performed on both sides of the Equation (5), by using the Riemann Liouville fractional order integral definition; that is.
ε ( t ) = σ η n t n Γ ( 1 + n ) , 0 n 1

2.2. Establishment and Solution of Fractional Order Creep Damage Model

In current research, the basic rheological mechanics models include elastomers, the viscoelastic body, the visco-elasticity plasticity body, the viscous body and the viscoplastic body. The rheological properties of materials can be described by the combination of strings and parallelism among them. The combined model has been widely used in the creep simulation research of asphalt, rock, polymer materials, concrete, and other materials because of its visualization and intuition capabilities [33,34,35].
The typical concrete compression and creep process consists of three stages: deceleration creep, steady creep, and accelerated creep. In the current study, the elastic-visco-elasticity plasticity-viscoplastic creep model (MSSB, Modified-Schofield-Scott- Blair) is used to describe the instantaneous and isovelocity creep of concrete. However, when concrete reaches the yield limit, the creep curve will exhibit accelerated unsteady characteristics. At this time, the MSSB model cannot describe the nonlinear creep characteristics of concrete, and it can only be studied by nonlinear theory. In this paper, the creep constitutive Equation based on the fractional derivative is constructed by using the Able dashpot, instead of the Newton dashpot in the MSSB model to accurately describe the three stages of deceleration creep, steady creep and accelerated creep. The constitutive model combination is shown in Figure 2.
To derive the constitutive equation, the main hypothesis conditions are given.
  • The small damage produced is neglected during the deceleration creep and steady creep stages, and only the damage caused by the accelerated creep stage is considered in the concrete indoor creep test.
  • The concrete damage evolution in the accelerated creep process is consistent with the change in the exponential function.
Based on the above constitutive models and assumptions, the stress-strain relationship between elastomers, the visco-elasticity plasticity body, and the viscoplastic body are derived as follows.
  • Stress-strain relationship of elastomer.
    ε e = σ E 0
    where E 0 is the elastomer elastic modulus.
  • Stress-strain relationship of visco-elasticity plasticity body.The visco-elasticity plasticity body is represented in parallel with the elastomer element, the Able dashpot and the friction element. The stress on the friction element is expressed by σ P 1 .
    σ P 1 = σ σ < σ S 1 σ S 1 σ σ S 1
    where σ is the total stress on the visco-elasticity plasticity body, and σ S 1 is the ultimate stress on the visco-elasticity plasticity body friction slider.
The total stress on the visco-elasticity plasticity body can be obtained from the combined model theory.
σ = σ d 1 + σ P 1
where σ d 1 is the stress on the Able dashpot.
When σ σ S 1
σ d 1 = 0 , ε v e = 0
when σ > σ S 1 .
σ = σ S 1 + E 1 ε v e + η 1 n d n ε v e d t n
where E 1 is the elastic modulus, and η 1 is the Able dashpot viscous coefficient.
Equation (11) can also be expressed as follows.
σ σ S 1 η 1 n = E 1 η 1 n ε v e + d n ε v e d t n
Laplace transformation of Equation (12).
E S = σ σ S 1 η 1 n s η 1 n s n η 1 n + E 1
Laplace inverse transformation of Equation (13).
ε v e = σ σ S 1 η 1 n k = 0 ( 1 ) k ( E 1 / η 1 ) k t n ( 1 + k ) n ( 1 + k ) Γ n ( 1 + k )
where n is the order of the fractional calculus and Γ ( · ) is the gamma function.
(3) Stress-strain relationship of viscoplastic body. The viscoplastic body is expressed in parallel with the friction element by the Able dashpot. The stress σ P 2 on the friction element is expressed by the Able dashpot, which satisfies the following requirements.
σ P 2 = σ σ < σ S 2 σ S 2 σ σ S 2
where σ is the total stress on the visco-elasticity plasticity body, and σ S 2 is the ultimate stress of the friction slider in the viscoplastic body. The total stress on the viscoplastic body can be obtained from the combination model theory.
σ = σ d 2 + σ P 2
where σ d 2 is the stress on the Able dashpot.
When σ σ S 2
ε v p = 0
when σ > σ S 2 .
σ = σ S 2 + η 2 n d n ε v p d t n
where σ S 2 is the yield stress, and η 2 is the viscosity coefficient of the Able dashpot.
According to the theory of fractional calculus [32,36], the following can be obtained:
ε v p = σ σ S 2 η 2 n t n Γ ( 1 + n )
During the creep test process, especially in the accelerated creep stage, the creep rate of the concrete specimens increases continuously until failure, which is caused by the expansion of internal damage and crack development in the concrete. Therefore, the viscosity coefficient is no longer a fixed value, but rather is a variable that decreases with the increase in the degree of damage. The damage variable is a negative exponential function [37] that considers the influence of damage growth on the creep process, which can be expressed as.
D = 1 exp ( α t )
where D is the damage variable, and its value range is 1], which corresponds to the concrete failure state without damage, and the complete damage, respectively. α is a material parameter, which controls the damage evolution process and is related to the applied load. At this time, the viscous coefficient changes with the damage variable, which can be expressed as:
η 2 n ( t , D ) = η 2 n ( 1 D ) = η 2 n e α t
then
ε v p = σ σ S 2 η 2 n e α t t n Γ ( 1 + n )
Combining the above three parts, the strain and nonlinear creep damage constitutive model of concrete, based on fractional calculus theory, is obtained as
ε = σ E 0 + σ σ S 1 η 1 n k = 0 ( 1 ) k ( E 1 / η 1 ) k t n ( 1 + k ) n ( 1 + k ) Γ n ( 1 + k ) σ S 1 σ σ S 2 σ E 0 + σ σ S 1 η 1 n k = 0 ( 1 ) k ( E 1 / η 1 ) k t n ( 1 + k ) n ( 1 + k ) Γ n ( 1 + k ) + σ σ S 2 η 2 n e α t t n Γ ( 1 + n ) σ > σ S 2

3. Determination of Concrete Nonlinear Creep Damage Constitutive Model Parameters

3.1. Determination of Parameter E 0

The parameter E 0 is the elastomer elastic modulus in the instantaneous state, which is calculated by the elastic strain section of the creep test at a certain stress level. The expression of parameter E 0 is as follows.
E 0 = σ 0 ε 0
where σ 0 is the initial stress, and ε 0 is the strain in the elastic stage.

3.2. Determination of Parameters E 1 , σ S 1 , and η 1

When the creep curve has an unloading curve, the three parameters E 1 , σ S 1 and η 1 of the creep model can be calculated from the strain when the creep tends to be stable under a certain level of horizontal stress and under residual strain after unloading [38]. The expressions of E 1 , σ S 1 , η 1 are as follows.
E 1 = σ 0 ε a ( ) + ε b ( ) , σ S 1 = E 1 ε b ( ) , η 1 = E 1 t ln ( 1 ε a ( t ) / ε a ( ) )
where ε a ( ) is the strain when the creep tends to be stable under a certain stress loading. ε b ( ) is a creep variable that cannot be recovered after unloading under a certain stress load. t is the creep time ε a ( t ) is the creep value corresponding to time t, where ( t , ε a ( t ) ) is any point in the creep viscoelastic stage. η 1 can select the average value of the calculation results from multiple data points.
A nonlinear regression analysis method [39] can be used to solve for parameters E 1 , σ S 1 , and η 1 when there is no unloading curve.

3.3. Determination of Parameters η 2 and σ S 2

The parameters η 2 and σ S 2 in the Viscoplasticity body can be obtained from the creep test curve rate in two groups when there is constant creep data in the experimental data. The expression is as follows.
ε ˙ c = σ 1 σ S 1 η 2 , ε ˙ d = σ 2 σ S 2 η 2
which are Calculated by Equation (26).
η 2 = σ 2 σ 1 ε ˙ d ε ˙ c , σ S 2 = ε ˙ d σ 2 ε ˙ c σ 1 ε ˙ d ε ˙ c
where σ 1 and σ 2 are two stress values of the creep test. ε ˙ c and ε ˙ d are creep rates at constant velocity under two stress levels. The concrete creep process will appear as third stage creep when the stress of the concrete is greater than its long-term strength. According to Equation (21), the relationship can be described by the exponential relationship. The viscosity coefficient decreases with the increase of time, and its variation law is related to the magnitude of the applied stress. At the same time, the greater the applied stress, the smaller the viscous coefficient. When the applied stress is large, the nonlinear regression analysis method can be used to solve for the parameters η 2 and σ S 2 due to the lack of two stress level isometric creep data.

4. Verification of Concrete Nonlinear Creep Damage Constitutive Model

4.1. Experimental Research on Concrete Creep under High Stress

Uniaxial compressive nonlinear non-linear creep and residual strain tests of normal concrete with deep storage nuclear waste were carried out by Narintsoa Ranaivomanana, et al. [40]. The specimen size was 100 mm × 100 mm × 500 mm. A series of uniaxial compression nonlinear creep tests under different stresses and maintaining 20 °C were carried out. The stress was controlled at 30% (about 20.9 MPa) and 50% (about 33.9 MPa). The cement class was CEM I 52.5R(CEM I is the abbreviation for Portland cement.), the modulus of elasticity was 41925 MPa, the modulus of elasticity was 4295 MPa, and the average density was 2410 kg/m3. The experimental curves are shown in Figure 3.
As shown in Figure 3, the concrete specimen nonlinear creep process can be divided into two stages under an uniaxial compression of stress levels of 20.9 MPa and 33.9 MPa: (1) the deceleration creep stage, with increasing strain and decreasing strain rate, and (2) the isovelocity creep stage with constant strain rate. The experimental curve characteristics can be described by using the elastomer and viscoelastic body in the model. Now, the parameters E 0 , E 1 , σ S 1 , η 1 , and n need to be solved. The numerical fitting analysis of Narintsoa Ranaivomanana, et al.’s [40] experimental results is carried out using the conducted model. The calculated model parameters are shown in Table 1, and the fitting results are shown in Figure 4. Where ‘-’ means that there is no such value. Accelerate creep has not occurred when the axial pressure is 20.9 MPa and 33.9 MPa, the concrete creep process consists of deceleration and isometric creep.
Narintsoa Ranaivomanana et al.’s [40] experimental curves show, a trend that is consistent with constructed concrete nonlinear creep damage evolution, and the fitting correlation coefficient is more than 0.96, as shown in Table 1 and Figure 4. The fitting results show that the concrete instantaneous elastic deformation after loading can be correctly described by the constitutive model. Meanwhile, the initial creep of stage 1 and the steady creep of stage 2 can also be correctly expressed.

4.2. Experimental Research on Concrete Creep Under Different Stresses

Compressive nonlinear creep and residual strain tests on C40 concrete specimens under different stresses were carried out by Can Tang et al. [41]. The concrete specimens were prisms with the size of 100 mm × l00 mm × 400 mm. The stress was controlled at 12.1 MPa, 16.9 MPa, 20.7 MPa, 25.7 MPa, 30.0 MPa, and 35.7 MPa. The cubic compressive strength and axial compressive strength was 50.4 MPa and 40.9 MPa, respectively, under a 28-day standard. The nonlinear creep test curve is shown in Figure 5.
Under stress of 16.9 MPa and 25.7 MPa, the instantaneous creep occurred first and then the isometric creep, and the creep was velocity first increased and then decreased, as shown in Figure 5a. The creep process includes elastic strain, deceleration creep and isovelocity creep. The experimental curve characteristics can be described using the elastomer and viscoelastic body in the model. At this time, the parameters E 0 , E 1 , σ S 1 and η 1 , and the fractional order n are selected according to the curve characteristics. Under stress levels of 30.0 MPa and 35.7 MPa, as shown in Figure 5b, the concrete specimen creep process consists of instantaneous creep, isovelocity creep, and accelerated creep, and the applied stress is greater than σ S 2 . The experimental curve characteristics can be described using the elastomers, viscoelastics, and viscous bodies in the model. Now, the parameters E 0 , E 1 , σ S 1 , η 1 , σ S 2 , and η 2 need to be solved, and the fractional order n can be selected according to the curve characteristics. The numerical fitting analysis of the test results are carried out using the conducted model. The calculation model parameters are shown in Table 2. The fitting results are shown in Figure 6.
Can Tang et al.’s [41] test curves show the same trend as with the concrete nonlinear creep damage evolution, and the fitting correlation coefficient is 0.96, which has a good correlation, as shown in Table 2 and Figure 6. The fractional order creep model represents the deceleration and isovelocity creep of the concrete, and it also express the accelerated creep stage by comparison. It has good consistency for the entire concrete creep process. In summary, the concrete creep damage evolution rule, under high and different levels of stress, can be well reflected by the constructed concrete creep damage constitutive model proposed in this paper. It is reasonable and feasible to describe the three stages of concrete creep deceleration creep, constant, and accelerated creep and the model has certain application and reference value. Meanwhile, the application of fractional calculus theory to improve traditional elements is also further illustrated, as this theory has become one of the important means to study the concrete creep nonlinearity. The research results are consistent with Zhilei He [42] and HW Zhou [36,43].

5. Parameter Sensitivity Analysis of Concrete Nonlinear Creep Damage Constitutive Model

Parameters σ, n, and α are the main factors affecting the concrete creep in the concrete creep damage constitutive model; these parameters are obtained by analyzing the experimental data curve and the theoretical curve. Therefore, the sensitivity analysis of the model parameters will be carried out by the control variable method in this section to study the influence of the above model parameters on concrete creep. The model parameters are as follows: σ = 30 · M P a , E 0 = 32.3 G P a , E 1 = 254 G P a , σ S 1 = 18.7 M P a , η 1 n = 246 G P a h , σ S 2 = 26.3 M P a , η 2 n = 854 G P a h , α = 0.001 and n = 0.38 . The concrete compression creep and residual strain experimental data of Can Tang are selected, and the stress level is 30.0 MPa.

5.1. Effect of Stress Level σ

Parameters E 0 = 32.3 G P a , E 1 = 254 G P a , σ S 1 = 18.7 M P a , η 1 n = 246 G P a h , σ S 2 = 26.3 M P a , η 2 n = 854 G P a h , α = 0.001 and n = 0.38 are selected and brought into Equation (23). The stress level varies: 28 MPa, 30 MPa, 32 MPa, and 34 MPa. The creep test results at different stress levels are shown in Figure 7.
The concrete creep deformation is greatly affected by the stress level, and the creep deformation will be increased as the stress level increases. However, stress has no obvious effect on the creep velocity. The higher the stress level, the greater the instantaneous deformation in the elastic stage, and the shorter the time required to enter steady state creep at a constant velocity. In addition, the form of the concrete creep curves is gradually changed to the exponential type as the stress level increased. When the stress level exceeds a certain level, a sudden change will occur. At this point, micro-cracks will occur, eventually causing macro-crack penetration failure. This is consistent with Narintsoa Ranaivomanana et al.’s test [40] of how high stress can accelerate the creep process.

5.2. The Influence of Fractional Order n

In order to analyze the sensitivity of the fractional order, the control variable method is used to ensure that other parameters remain unchanged. The fractional order is gradually increased from 0.1 to 0.6, and a series of creep curves with different orders are obtained, as shown in Figure 8.
From Figure 8, we can see that the duration of instantaneous creep is shortened with the increase of the fractional order and the duration of deceleration and isovelocity creep, and of the concrete entering the stage of accelerated creep, is correspondingly also shortened. The concrete creep deformation is increased gradually as the fractional order increases. In addition, the concrete creep velocity is affected by the fractional order number; that is, the higher the fractional order, the greater the creep rate at each stage.

5.3. Effect of Material Parameter α

Considering the change in the damage variable in a viscoplastic body, the control variable method is used to obtain the creep curves under different α, by changing the value of α, as shown in Figure 9. As material parameter α increases, the sooner the accelerated creep phase begins, and the greater the creep velocity will be. The power index is the main factor that influences creep velocity. The larger the power exponential value, the more prone the concrete will be to accelerated creep.
The classical rheological model and the Fractional order theory are introduced to inform the concrete nonlinear creep damage constitutive model after the results of the above research. The conducted constitutive model is a viscoelastic-plastic model.
It is assumed that the damage in the deceleration creep stage and in the stable creep stage is neglected, and only the damage in the accelerated creep stage is considered. Also, the accelerated creep rate conforms with changes in the exponential function charge in the indoor concrete creep tests. In practice, the above two assumptions may not be agreed with each other. Therefore, the concrete fractional nonlinear creep damage constitutive model derived from this assumption differs from the actual situation. However, the concrete specimen damage data contributes little to the constitutive equation in either the initial stage or the stable creep stage in the actual uniaxial compression creep tests. The constructed concrete creep damage model is more able to meet the needs of actual working conditions and tests that consider accelerated creep stage damage as the main research object. In addition, the constructed concrete creep damage model parameters have specific physical significance in this paper, and the entire concrete creep process in practical engineering can be reflected and described by the proposed model, which has certain engineering practical and reference values.

6. Conclusions

To better describe the concrete non-linear creep damage evolution rule, and better solve the concrete creep in concrete engineering, Fractional Calculus Theory is introduced to inform a concrete nonlinear creep damage constitutive model. The conducted constitutive model is verified and analyzed for sensitivity by means of the concrete creep test. The following conclusions may be drawn:
  • The Visco-elasticity Plasticity Rheological Theory, the Riemann Liouville Theory and the Combined Model Theory are quoted to propose the stress-strain relationship expressions of the elastomer body, the visco-elasticity plasticity body, and the viscoplastic body. The nonlinear creep damage evolution equation of the concrete fractional order is conducted with consideration of damage evolution, and the parameter specific calculation method in the constructed model is given.
  • The constructed concrete creep damage constitutive model in this paper is validated in high stress and different stress situations that involve uniaxial compression creep; data from the residual strain tests of Narintsoa Ranaivomanana et al. [40] and Can Tang [41] is considered. The results show that the conducted model can accurately describe the entire concrete creep process, and that it could be used as a reference for concrete non-linear creep.
  • To determine the influence degree of the model parameters, the sensitivity analysis of the main influence parameters in the concrete creep damage constitutive model, including σ, material parameter α, and fractional order n, are analyzed. The results show that σ and n are positively correlated with the concrete creep deformation and rate. The creep deformation and rate will increase as σ and n increase. α is related to the time at which the concrete enters the accelerated creep stage, and the larger the value of material parameter α, the earlier the concrete will enter the accelerated creep stage.

Author Contributions

Conceptualization, C.Z., Z.Z., S.Z. and J.L.; Methodology, C.Z., Z.H., D.Z. and J.L.; Software, Z.H., D.Z. and J.L.; Validation, C.Z., Z.H. and J.L.; Formal Analysis, C.Z., Z.Z., S.Z. and D.Z.; Investigation, Z.H., D.Z. and S.M.; Resources, C.Z., S.Z. and Z.Z.; Data Curation, C.Z., Z.H. and D.Z.; Writing-Original Draft Preparation, C.Z.; Writing-Review & Editing, C.Z., Z.Z., S.Z., J.L. and Z.H.; Visualization, J.L. and D.Z.; Supervision, C.Z., Z.Z., S.Z. and Z.H.; Project Administration, Z.Z.; Funding Acquisition, C.Z., Z.Z., S.Z. and D.Z.

Funding

This research was partially carried out with financial support from [the National Basic Research Program of China (973 Program)] grant number [2015CB057903] and supported by [the Fundamental Research Funds for the Central Universities and Postgraduate Research & Practice Innovation Program of Jiangsu Province] grant number [2019B74214,SJKY19_04533] and project support from [the National Natural Science Foundation of China] grant number [41831278, 51579081 and 51379065].

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

References

  1. Aitcin, P.C. Concrete structure, properties and materials. Preticehall Int. 2006, 13, 499. [Google Scholar] [CrossRef]
  2. Luu, C.H.; Mo, Y.L. Development of CSMM-based shell element for reinforced concrete structures. Eng. Struct. 2017, 132, 778–790. [Google Scholar] [CrossRef]
  3. Bažant, Z.P.; Hubler, M.H.; Yu, Q. Damage in prestressed concrete. In Handbook of Damage Mechanics; Springer: New York, NY, USA, 2016; pp. 515–560. [Google Scholar] [CrossRef]
  4. Balevicius, R.; Marciukaitis, G. Linear and non-linear creep models for a multi-layered concrete composite. Arch. Civ. Mech. Eng. 2013, 13, 472–490. [Google Scholar] [CrossRef]
  5. Hamed, E.; Lai, C. Geometrically and materially nonlinear creep behaviour of reinforced concrete columns. Structures 2016, 5, 1–12. [Google Scholar] [CrossRef]
  6. Majorana, C.E. Mathematical modeling of creep and shrinkage of concrete. Int. J. Numer. Methods Biomed. Eng. 1989, 5, 423. [Google Scholar]
  7. Sprince, A.; Korjakins, A.; Pakrastinsh, L. Time-dependent behavior of high performance fiber-reinforced concrete. Adv. Mater. Res. 2013, 705, 75–80. [Google Scholar] [CrossRef]
  8. CEB-FIP Model Code 1990. Design Code; Thomas Teloford Services Ltd: London, UK, 1991.
  9. Eurocode BS EN 1992-1-1. Design of Concrete Structures; Department for Communities and Local Government: London, UK, 2004.
  10. ACI Committee 209. Prediction of Creep, Shrinkage and Temperature Effects in Concrete Structures; American Concrete Institute: Farmington Hills, MI, USA, 1999. [Google Scholar]
  11. Gardner, N.J.; Lockman, M.J. Design provisions for drying shrinkage and creep of normal-strength concrete. ACI Mater. J. 2001, 98, 159–167. [Google Scholar]
  12. Bazant, Z.P.; Baweja, S. Creep and shrinkage prediction model for analysis and design of concrete structures—Model b 3. Mater. Struct. 1995, 28, 357–365. [Google Scholar]
  13. Bazant, Z.P.; Baweja, S. Justification and refinements of model B3 for concrete creep and shrinkage 2. Updating and theoretical basis. Mater. Struct. 1995, 28, 488–495. [Google Scholar] [CrossRef]
  14. Hubler, M.H.; Wendner, R.; Bažant, Z.K. Statistical justification of model b4 for drying and autogenous shrinkage of concrete and comparisons to other models. Mater. Struct. 2015, 48, 797–814. [Google Scholar] [CrossRef]
  15. Wendner, R.; Hubler, M.H.; Bazant, Z.P. Optimization method, choice of form and uncertainty quantification of model b4 using laboratory and multi-decade bridge databases. Mater. Struct. 2015, 48, 771–796. [Google Scholar] [CrossRef]
  16. Ghasemzadeh, F.; Manafpour, A.; Sajedi, S.; Shekarchi, M.; Hatami, M. Predicting long-term compressive creep of concrete using inverse analysis method. Constr. Build. Mater. 2016, 124, 496–507. [Google Scholar] [CrossRef]
  17. Kindrachuk, V.M.; Galanov, B.A.; Kartuzov, V.V.; Dub, S.N. Refined model of elastic nanoindentation of a half-space by the blunted Berkovich indenter accounting for tangential displacements on the contact surface. J. Mater. Sci. 2009, 44, 2599–2609. [Google Scholar] [CrossRef]
  18. Grenacher, M.; Thurlimann, B. Creep buckling of reinforced concrete columns. Int. Assoc. Bridg. Struct. Eng. 1970, 6, 275–282. (In German) [Google Scholar]
  19. Bažant, Z.P.; Tsubaki, T. Nonlinear creep buckling of reinforced concrete columns. J. Struct. Div. 1980, 106, 2235–2257. [Google Scholar]
  20. Kwak, H.G.; Kim, J.K. Nonlinear behaviour of slender RC columns: Numerical formulation. Constr. Build. Mater. 2006, 20, 527–537. [Google Scholar] [CrossRef]
  21. Radnić, J.; Matešan, D.; Smilović, M. Nonlinear creep model for concrete in analysis of Plates and Shells. Adv. Struct. Mater. 2012, 1, 163–173. [Google Scholar]
  22. Briffaut, M.; Benboudjema, F.; Laborderie, C.; Torrenti, J.-M. Creep consideration effect on meso-scale modeling of concrete hydration process and consequences on the mechanical behavior. J. Eng. Mech. 2013, 139, 1808–1817. [Google Scholar] [CrossRef]
  23. Abdulkadir, C.A.; Arslan, A.; Rüstem, G. Mesoscale simulation of cement-based materials’ time-dependent behavior. Comput. Mater. Sci. 2008, 41, 20–26. [Google Scholar]
  24. Sherief, H.H.; El-Sayed, A.M.A. Fractional order theory of thermoelasticity. Int. J. Solids Struct. 2010, 47, 269–275. [Google Scholar] [CrossRef]
  25. Papoulia, K.D.; Panoskaltsis, V.P.; Kurup, N.V.; Korovajchuk, I. Rheological representation of fractional order viscoelastic material models. Rheol. Acta 2010, 49, 381–400. [Google Scholar] [CrossRef]
  26. Suzuki, J.L.; Zayernouri, M.; Bittencourt, M.L.; Karniadakis, G.E. Fractional-order uniaxial visco-elasto-plastic models for structural analysis. Comput. Methods Appl. Mech. Eng. 2016, 308, 443–467. [Google Scholar] [CrossRef]
  27. Bouras, Y.; Zorica, D.; Atanacković, T.M.; Vrcelj, Z. A non-linear thermo-viscoelastic rheological model based on fractional derivatives for high temperature creep in concrete. Appl. Math. Model. 2018, 55, 551–568. [Google Scholar] [CrossRef]
  28. Beltempo, A.; Zingales, M.; Bursi, O.S.; Deseri, L. A fractional-order model for aging materials: An application to concrete. Int. J. Solids Struct. 2018, 138, 13–23. [Google Scholar] [CrossRef]
  29. Adolfsson, K.; Enelund, M.; Olsson, P. On the fractional order model of viscoelasticity. Mech. Time-Depend. Mater. 2005, 9, 15–34. [Google Scholar] [CrossRef]
  30. Koeller, R.C. Applications of fractional calculus to the theory of viscoelasticity. Trans. ASME J. Appl. Mech. 1984, 51, 299–307. [Google Scholar] [CrossRef]
  31. Yin, D.; Ren, J.; He, C.; Chen, Q. A new creep model element for geomaterials. Chin. J. Rock Mech. Eng. 2007, 26, 1899–1903. [Google Scholar]
  32. Zhou, H.; Wang, C.; Duan, Z.; Zhang, M.; Liu, J. Time-based fractional derivative approach to creep constitutive model of salt rock. Sci. Sin. Phys. Mech. Astronomica 2012, 42, 310–318. [Google Scholar] [CrossRef]
  33. Gama, D.A.; Júnior, J.M.R.; De Melo, T.J.A.D.; Rodrigues, J.K.G. Rheological studies of asphalt modified with elastomeric polymer. Constr. Build. Mater. 2016, 106, 290–295. [Google Scholar] [CrossRef]
  34. Jun, S. Rock rheological mechanics and its advance in engineering applications. Chin. J. Rock Mech. Eng. 2007, 26, 1081–1106. [Google Scholar]
  35. Berdichevsky, V.L.; Herman, J.N. On rheology of cross-linked polymers: Slippage of polymer chains and its macroscopic modeling. Int. J. Eng. Sci. 2016, 100, 183–198. [Google Scholar] [CrossRef]
  36. Zhou, H.W.; Wang, C.P.; Han, B.B.; Duan, Z.Q. A creep constitutive model for salt rock based on fractional derivatives. Int. J. Rock Mech. Min. Sci. 2011, 48, 116–121. [Google Scholar] [CrossRef]
  37. Han, B.; Xie, H.B.; Zhu, L.; Jiang, P. Nonlinear model for early age creep of concrete under compression strains. Constr. Build. Mater. 2017, 147, 203–211. [Google Scholar] [CrossRef]
  38. Xia, C.; Xu, C.; Wang, X.; Zhang, C. Method for parameters determination with unified rheological machanical model. J. Rock Mech. Eng. 2009, 28, 425–432. [Google Scholar]
  39. Wang, P.; Jiang, L.; Li, X.; Qin, Q.; Wang, E. Physical Simulation of Mining Effect Caused by A Fault Tectonic. Arab. J. Geosci. 2018, 11, 741. [Google Scholar] [CrossRef]
  40. Ranaivomanana, N.; Stéphane, M.; Turatsinze, A. Tensile, compressive and flexural basic creep of concrete at different stress levels. Cem. Concr. Res. 2013, 52, 1–10. [Google Scholar] [CrossRef]
  41. Can, T. Research on compressive creep model and stress-strain relationship of concrete. Ph.D. Thesis, Harbin Institute of Technology, Harbin, China, 2016. [Google Scholar]
  42. He, Z.; Zhu, Z.; Ni, X.; Li, Z. Shear creep tests and creep constitutive model of marble with structural plane. Eur. J. Environ. Civ. Eng. 2017, 1–19. [Google Scholar] [CrossRef]
  43. Zhou, H.W.; Wang, C.P.; Mishnaevsky, L., Jr.; Duan, Z.Q.; Ding, J.Y. A fractional derivative approach to full creep regions in salt rock. Mech. Time Depend. Mater. 2013, 17, 413–425. [Google Scholar] [CrossRef]
Figure 1. (a) Newton dashpot; (b) Able dashpot.
Figure 1. (a) Newton dashpot; (b) Able dashpot.
Materials 12 01505 g001
Figure 2. Creep Constitutive model.
Figure 2. Creep Constitutive model.
Materials 12 01505 g002
Figure 3. The curve of the concrete creep tests under different stress.
Figure 3. The curve of the concrete creep tests under different stress.
Materials 12 01505 g003
Figure 4. The fitting curve of the test curve and the model at (a) 20.9 MPa and (b) 33.9 MPa).
Figure 4. The fitting curve of the test curve and the model at (a) 20.9 MPa and (b) 33.9 MPa).
Materials 12 01505 g004
Figure 5. The curve of concrete creep test under different stress. (a) 16.9 MPa and 25.7 MPa; (b) 30 MPa and 35.7 MPa.
Figure 5. The curve of concrete creep test under different stress. (a) 16.9 MPa and 25.7 MPa; (b) 30 MPa and 35.7 MPa.
Materials 12 01505 g005
Figure 6. Fitting curves of test curves and built models. (a) σ = 16.9 MPa; (b) σ = 25.7 MPa; (c) σ = 30.0 MPa; (d) σ = 35.7 MPa.
Figure 6. Fitting curves of test curves and built models. (a) σ = 16.9 MPa; (b) σ = 25.7 MPa; (c) σ = 30.0 MPa; (d) σ = 35.7 MPa.
Materials 12 01505 g006
Figure 7. Sensitivity analysis under different stress σ .
Figure 7. Sensitivity analysis under different stress σ .
Materials 12 01505 g007
Figure 8. Sensitivity analysis under different orders of n.
Figure 8. Sensitivity analysis under different orders of n.
Materials 12 01505 g008
Figure 9. Sensitivity analysis under different material parameters α.
Figure 9. Sensitivity analysis under different material parameters α.
Materials 12 01505 g009
Table 1. Parameter table of concrete creep model.
Table 1. Parameter table of concrete creep model.
Stress/MPa E 0 / GPa σ S 1 / MPa η 1 n / GPa · h E 1 / GPa σ S 2 / MPa η 2 n / GPa · h α/h-in
20.945.213.52352630.31
33.940.124.01913320.36
Table 2. Parameter table of concrete creep model.
Table 2. Parameter table of concrete creep model.
Stress/MPa E 0 / GPa σ S 1 / MPa η 1 n / GPa · h E 1 / GPa σ S 2 / MPa η 2 n / GPa · h α/h-in
16.938.58.71712810.3
25.742.814.01053170.29
30.032.318.724625426.38540.0010.38
35.737.820.223938430.96880.00130.48

Share and Cite

MDPI and ACS Style

Zhang, C.; Zhu, Z.; Zhu, S.; He, Z.; Zhu, D.; Liu, J.; Meng, S. Nonlinear Creep Damage Constitutive Model of Concrete Based on Fractional Calculus Theory. Materials 2019, 12, 1505. https://doi.org/10.3390/ma12091505

AMA Style

Zhang C, Zhu Z, Zhu S, He Z, Zhu D, Liu J, Meng S. Nonlinear Creep Damage Constitutive Model of Concrete Based on Fractional Calculus Theory. Materials. 2019; 12(9):1505. https://doi.org/10.3390/ma12091505

Chicago/Turabian Style

Zhang, Cong, Zhende Zhu, Shu Zhu, Zhilei He, Duan Zhu, Jinzhong Liu, and Songsong Meng. 2019. "Nonlinear Creep Damage Constitutive Model of Concrete Based on Fractional Calculus Theory" Materials 12, no. 9: 1505. https://doi.org/10.3390/ma12091505

APA Style

Zhang, C., Zhu, Z., Zhu, S., He, Z., Zhu, D., Liu, J., & Meng, S. (2019). Nonlinear Creep Damage Constitutive Model of Concrete Based on Fractional Calculus Theory. Materials, 12(9), 1505. https://doi.org/10.3390/ma12091505

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop