Next Article in Journal
Laser Micro-Texturing of Sintered Tool Materials Surface
Previous Article in Journal
3D Multi-Branched SnO2 Semiconductor Nanostructures as Optical Waveguides
Previous Article in Special Issue
Mesh-Based and Meshfree Reduced Order Phase-Field Models for Brittle Fracture: One Dimensional Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Stochastic Damage Model for Bond Stress-Slip Relationship of Rebar-Concrete Interface under Monotonic Loading

1
School of Civil Engineering, Central South University, Changsha 410075, China
2
National Engineering Laboratory for High Speed Railway Construction, Changsha 410075, China
3
Engineering Technology Research Center for Prefabricated Construction Industrialization of Hunan Province, Changsha 410075, China
*
Author to whom correspondence should be addressed.
Materials 2019, 12(19), 3151; https://doi.org/10.3390/ma12193151
Submission received: 4 July 2019 / Revised: 12 September 2019 / Accepted: 24 September 2019 / Published: 26 September 2019
(This article belongs to the Special Issue Randomness and Uncertainty)

Abstract

:
The stochastic bond stress-slip behavior is an essential topic for the rebar-concrete interface. However, few theoretical models incorporating stochastic behavior in current literature can be traced. In this paper, a stochastic damage model based on micro-mechanical approach for bond stress-slip relationship of the interface under monotonic loading was proposed. In order to describe the mechanical behaviors of the rebar-concrete interface, a microscopic damage model was proposed. By introducing a micro-element consists of parallel spring element, friction element and a switch element, the model is formulated. In order to reflect the randomness of the bond stress-slip behavior contributed by the micro-fracture in the interface, a series of paralleled micro-elements are adopted with the failure threshold of individual spring element is set as a random variable. The expression of both mean and variance for the bond stress-slip relationship was derived based on statistical damage mechanics. Furthermore, by utilizing a search heuristic global optimization algorithm (i.e., a genetic algorithm), parameters of the proposed model are able to be identified from experimental results, which a lognormal distribution has adopted. The prediction was verified against experimental results, and it reveals that the proposed model is capable of capturing the random nature of the micro-structure and characterizing the stochastic behavior.

1. Introduction

Although concrete is the most widely material in construction field around world, the bond stress-slip behavior of the rebar-concrete interface has not been fully understood. Among the several focused topics towards the bond stress-slip behavior, the damage of the interface is an essential issue, which has been found in conclusion to be the main reason for degradation of concrete structures [1,2]. Thus, misestimating could be occurred regarding safety and economic design/analysis of relevant structures when consideration of such behavior is absent.
The complexity of bond stress-slip behavior for rebar-concrete interface commonly originated from two essential characteristics; that is, the nonlinearity and randomness. In recent decades, the nonlinearity for bond stress-slip behavior is treated as a fundamental problem, which the constitutive relationship has drawn intensified study and several celebrated works have been conducted by researchers. However, referring to the randomness, research on stochastic properties of relevant material is comparatively blank. In detail, the causes for such behavior can be mainly attributed to the complex inhomogeneous microstructure of the interface, due to the irregular shape, random size and distribution of aggregates, the sedimentation of cement particles, the accumulation of pore water, the curing conditions-induced undetermined hydration, the random corrosion of the rebar, the random size of the ribs, the different location of the concrete and rebar in the structures, and the varied skill of the workers.
Presently, unlike the concrete materials, research on the mechanical behaviors of bond stress-slip to the rebar-concrete interface is relatively rare. Besides, few literature evolved with the theoretical model regarding stochastic behavior can be traced. From a summary standpoint, currently, the researches refer to the theoretical model of the mechanical behaviors of bond stress-slip are able to be classified into three main categories: empirical models, theoretical analysis models and macro-mechanical models.
In detail, the empirical models are usually developed based on the observation of the experimental results, and could also be sub-divided into two types: the segmental function model and continuous function model. Precisely, the most representative contributions of the segmental function models listed in this work were conducted by Alsiwat, et al. [3], Haraji, et al. [4], Eligehausen, et al. [1], Xu [5] and Wang, et al. [6]. For example, Xu [5] divided the bond stress-slip curves into five sections: micro-slip section, slip section, splitting section, descending section and residual section, and thus, the segmental function model of the corresponding sections was proposed. In addition, the representative works of continuous function models were developed by Lutz and Gergely [7], Nilson [8], Mirza and Houde [9], Kankam [10], Di [11], Jin, et al. [12] and Teng, et al. [13]. Although most of the empirical models are able to describe the bond stress-slip behavior, the comprehensive micro-damage mechanism in a physical sense was absence to some extent.
Although aforementioned empirical models may have the advantage of simplifying relative design/analysis process towards certain problems, such formulated equations can be only used for representing specific observed phenomenon. Therefore, a number of theoretical attempts have also been devoted to the modelling of the bond stress-slip behaviors. For instance, considering the force balance of the interface under loading, the stress-strain relationship of the steel bars and concrete and the relationship between slip and strains, Somayaji and Shah [14] proposed the differential equations, corrected the solution of these equations, and obtained the relation model between the local slip and the buried rebar length. Yankelevsky [15] developed a model for interface between un-cracked concrete and deformed steel bars, derived the second-order differential equations dependent on the tensile force of the steel bar, combined the equations with the boundary conditions, further the relationship between bond stress and distance in the longitudinal direction of steel bars was predicted. Zhao and Xiao [16] developed a bond stress-slip model before peak bond stress, formulated based on a wedge-shaped model and mechanism for the interface. Song and Zhao [17] proposed a model based on the stress balance, deformation coordination, and physical conditions by considering the influence of material characteristics, concrete layer thickness, crack spacing, and distance from the crack section. The theoretical analysis models [14,15,16,17,18] considered the micro mechanisms of the interface behaviors, however, they are too complex to apply to engineering including finite element analysis. Nevertheless, there are a number of certain conditional assumptions for the model development. Unfortunately, these assumptions generally did not consider the random nature of the interface micro-structures.
Besides the above-mentioned two categories of theoretical models, the establishment of macroscopic mechanical models were also attracting extensive attentions from researchers. For instance, Alfano, et al. [19,20] established a thermodynamic-based cohesive zone model considering the damage-friction evolution with unilateral contact, which can be used for analysis of rebar-concrete interfaces and cracks in concrete-rock foundation interfaces of concrete dams, etc. Based on continuum damage mechanics and considering the accumulation of interfacial damage of reinforced concrete, Soh, et al. [21] developed a bond stress-slip constitutive model by adopting Weibull distribution function. In view of above introduced models, macroscopic mechanical models, especially the continuum damage mechanical models have been gradually focused by relevant researchers in recent. These models have considered different mechanisms for interface deterioration of reinforced concrete (e.g., damage, friction and interlocking mechanism) based on thermodynamics. However, even though the nonlinearity was effectively addressing and modelling by such models, they are still in deficiency of describing the stochastic properties exhibited in the bond stress-slip behavior.
Therefore, this work aims at developing a stochastic damage constitutive model for bond stress-slip relationship of the rebar-concrete interface in order to capture the random nature of the interface and characterize the stochastic constitutive behaviors. The outline is listed as follows, the characteristics of the microscopic element model of rebar-concrete interface is first introduced, and then a stochastic damage model is proposed based on such microscopic model. Afterward, the expressions of mean and variance of bond stress-slip relationship is derived in Section 2. In addition, random variable parameter identification is conducted by adopting a search heuristic global optimization algorithm (i.e., a genetic algorithm), which the detailed illustration is listed in Section 3. In Section 4, the proposed model is verified against experimental results. Our conclusions are finally given in Section 5.

2. Stochastic Damage Model of Bond Stress-Slip Relationship

2.1. Microscopic Element Model

For the purpose of mimicking the mechanical properties of the interface between the rebars and concrete evolved with the interclocking, damage (cracking) and friction, a new bond stress-slip model was proposed based on microscopic elements (see Figure 1). Concretely, the single microscopic element consists of a spring element and a friction element. Therefore, the elastic deformation, debonding and damage (cracking) is able to be modelled by the spring element and the mutual friction and slipping behavior can be characterized by the friction element.
It is assumed that when external load is applied on the interface, the spring element initially undergoes an elastic deformation process, and then breaks (debonding) when the deformation reaches its fracture threshold. As a result, the friction element then experiences a slipping process, which represents the friction between the rebar and concrete. A switch element was also introduced for controlling the work of the spring element and friction element. During the elastic deformation process, the switch element is in opening state, which reveals that the friction element was isolated. After the failure of the spring element, the switch element was closed and caused the friction element began to slide. Therefore, the entire response of the interface are able to be characterized when subject to external loading. Physically, it is worthy of note that such a switch element is not a mechanical element and the aim for introducing it is only for purpose of demonstrating the transformation process from debonding to the frictional sliding in microscopic scale.
Specifically, in this work, in order to describe the randomness of micro-fracture behaviors in the interface between the rebar and concrete, the random variable theory is applied. By presuming the fracture threshold Δ of spring elements obey certain probability densities, the randomness of the failure events are able to be depicted. Therefore, the state of the switch element (open or close) can be determined by comparison between the value of relative displacement between the rebars and concrete s and the fracture threshold Δ. That is, when s < Δ, the switch element is in an open state, otherwise it is closed.
Hence, the control function of the switch element can be expressed as
S W I T C H = H ( s Δ ) = { 1 s Δ , switch on 0 s < Δ , switch off
where H (x) is the Heaviside function, s is the slip.

2.2. Mechanical Behavior of Individual Microscopic Element

Thus, the bond stress-slip relationship for individual microscopic element is able to be divided into two phases:
i. Before the failure of the spring element,
τ e = G s , s < Δ
where G is the stiffness of the bond-slip relationship.
ii. After the failure of the spring element,
τ r = β τ u = β G Δ , s Δ
where β is the coefficient of friction and τu is the failure stress of the spring τu = GΔ.
Therefore, the entire bond stress-slip relationship of individual microscopic element can be expressed as follows (see Figure 2):
τ = H ( Δ s ) G s + β G Δ H ( s Δ )

2.3. Mechanical Behavior of Parallel System of Microscopic Elements

In order to mimic the corresponded relationship of the rebar-concrete interface, a representative interface element (RIE) in mesoscale was firstly introduced in this work, which consists of a parallel system of microscopic elements (see Figure 3). In detail, based on Equation (4), the bond stress of the k-th microscopic element was defined such that (k = 1, 2, 3, ···, N)
τ k = H [ Δ ( x k ) s ] G k s + H [ s Δ ( x k ) ] β k G k Δ
where xk denotes the location of the k-th microscopic element in the coordinate.
By assuming that the stiffness G and the coefficient of friction β of individual element is equal, the average bond stress for parallel system is obtained as follows:
τ = 1 N k = 1 N H [ Δ ( x k ) s ] G s + 1 N k = 1 N H [ s Δ ( x k ) ] β G Δ ( x k )
Taking the limit of Equation (6) as N approaches infinity, one obtains
τ = 0 1 H [ Δ ( x ) s ] G s d x + 0 1 β G Δ H [ Δ ( x ) s ] d x
In addition, the fracture threshold Δ is set as a random variable, with density function f(Δ) and distribution function F(Δ) = ʃ0+hf(Δ)dΔ, as illustrated in Figure 4.
Since the failures of spring element represent the damage of the interface, the damage variable can be defined as the ratio of damage area to total area of the parallel system [22,23] as follows:
d = A d A = 1 N k = 1 N H [ s Δ ( x k ) ]
When N approaches to infinity, the damage variable is able to be obtained as follows,
d ( s ) = 0 1 H [ s Δ ( x ) ] d x
By substituting Equation (9) to Equation (7), a damage model for the bond stress-slip relationship is derived such that
τ = ( 1 d ) G s + β G 0 1 Δ H ( s Δ ) d x
especially if the damage variable d = 0 implies the interface is in a non-damage state and τ = Gs. When the damage variable d = 1, it represents that the interface is completely damaged and only residual frictional stress is existed, that is, τ =βG01 Δ(x)dx.
Hence, the mean value of the damage variable d can be expressed as follows:
E [ d ( s ) ] = 0 0 1 H [ s Δ ( x ) ] f ( Δ ) d x d Δ = 0 1 0 s f ( Δ ) d x d Δ = 0 1 F ( s ) d x = F Δ ( s )
The square of the mean value of the damage variable is derived such that,
E [ d 2 ( s ) ] = 0 0 [ 0 1 0 1 H ( s Δ 1 ) H ( s Δ 2 ) d x 1 d x 2 ] f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 = 0 1 0 1 [ 0 s 0 s f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 ] d x 1 d x 2 = 0 s 0 s f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2
where fΔ1, Δ2) is the two-dimensional joint probability density function of the fracture threshold for spring elements.

2.4. Mean Value of Stochastic Mechanical Responses

Considering the properties between the expected operator and the integral operator, the following results can be obtained:
μ τ ( s ) = μ [ 0 1 H ( Δ s ) G s d x + 0 1 β G Δ H ( s Δ ) d x ] = μ [ 0 1 H ( Δ s ) G s d x ] + μ [ 0 1 β G Δ H ( s Δ ) d x ] = 0 1 μ [ H ( Δ s ) G s ] d x + 0 1 μ [ β G Δ H ( s Δ ) ]   d x
where
μ { H [ s Δ ( x ) ] } = 1 × P [ H ( s Δ ( x ) ) = 1 ] + 0 × P [ H ( s Δ ( x ) ) = 0 ] = P [ H ( s Δ ( x ) ) = 1 ] = P [ s Δ ( x ) 0 ] = P [ Δ ( x ) s ] = 0 s f ( Δ ) d Δ = F ( s )
μ { H [ Δ ( x ) s ] } = 1 0 s f ( Δ ) d Δ = 1 F ( s )
P denotes the probability.
Therefore, the mean function of bond stress can be expressed as
μ τ ( s ) = 0 1 μ [ H ( Δ s ) G s ] d x + 0 1 μ [ β G Δ H ( s Δ ) ] d x = G s 0 1 ( 1 0 s f ( Δ ) d Δ ) d x + β G 0 1 0 s Δ f ( Δ ) d Δ d x

2.5. Variance Value of Stochastic Mechanical Responses

The variance function of stochastic mechanical responses of the parallel system is expressed as follows,
V 2 ( τ ) = E ( τ 2 ) E 2 ( τ )
By combining Equations (9) and (10), it is able to obtain E(τ2) such that,
E ( τ 2 ) = E { ( 1 d ) 2 G 2 s 2 + β 2 G 2 [ 0 1 Δ H ( s Δ ) d x ] 2 + 2 β G 2 s ( 1 d ) 0 1 Δ H ( s Δ ) d x } = G 2 s 2 [ 1 + E ( d 2 ) 2 E ( d ) ] + β 2 G 2 0 0 [ 0 1 Δ 1 H ( s Δ 1 ) d x 1 0 1 Δ 2 H ( s Δ 2 ) d x 2 ] f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 = 2 β G 2 s { 0 0 1 Δ H ( s Δ ) f ( Δ ) d x d Δ 0 0 [ 0 1 Δ H 1 ( s Δ 1 ) d x 1 0 1 H ( s Δ 2 ) d x 2 ] f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 }
By substituting Equations (11) and (12) to Equation (18), E(τ2) is able to be derived as follows
E ( τ 2 ) = G 2 s 2 [ 1 + 0 s 0 s f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 2 0 s f ( Δ ) d Δ ] + β 2 G 2 0 s 0 s Δ 1 Δ 2 f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 + 2 β G 2 s 0 s Δ f ( Δ ) d Δ 2 β G 2 s 0 s 0 s Δ 1 f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2
Therefore, by substituting Equations (16) and (19) to Equation (17), the variance value of the stochastic mechanical responses is able to be obtained as follows,
V 2 ( τ ) = G 2 s 2 [ 1 + 0 s 0 s f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 2 0 s f ( Δ ) d Δ ] + β 2 G 2 0 s 0 s Δ 1 Δ 2 f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 + 2 β G 2 s 0 s Δ f ( Δ ) d Δ 2 β G 2 s 0 s 0 s Δ 1 f Δ ( Δ 1 , Δ 2 ) d Δ 1 d Δ 2 [ G s 0 1 ( 1 0 s f ( Δ ) d Δ ) d x + β G 0 1 0 s Δ f ( Δ ) d Δ d x ] 2
By presuming the fracture threshold Δ of obeys the lognormal distribution [22,24] and defining the mean and the variance of the fracture threshold as μΔ and σΔ, let Z(x) be a homogeneous normal function with the mean value and standard deviation (λ, ζ2), it is derived that
Z ( x ) = ln Δ ( x )
λ = E [ ln Δ ( x ) ] = ln ( μ 1 + σ Δ 2 / μ Δ 2 )
ζ 2 = var [ ln Δ ( x ) ] = ln ( 1 + σ Δ 2 / μ Δ 2 )
Therefore, the distributions can be obtained respectively as follows,
f ( Δ ) = 1 2 π Δ ζ exp [ 1 / 2 ( ln Δ λ ζ ) 2 ]
f Δ ( Δ 1 , Δ 2 ) = 1 2 π ζ 1 ζ 2 1 ρ 2 Δ 1 Δ 2 exp { 1 2 ( 1 ρ 2 ) [ ( ln Δ 1 λ 1 ) 2 ζ 1 2 2 ρ ( ln Δ 1 λ 1 ) ( ln Δ 2 λ 2 ) ζ 1 ζ 2 + ( ln Δ 2 λ 2 ) 2 ζ 2 2 ] }
Consequently, the properties of the fracture threshold can be described by the parameters λ, ζ, λ1, λ2, ζ1, ζ2 and ρ.

3. Random Variable Parameter Identification Based on Genetic Algorithm

To verify the foregoing analytical model proposed in this work, a random variable parameter identification process is developed based on genetic algorithm. Genetic algorithms (GAs) is a search method relative to the theory of evolution, including the processes of reproduction, crossover, mutation, and selection, create populations of solutions for optimizing an objective function. The applicability of GAs was widely accepted by relevant researchers for dealing with the problems of optimizing analysis.
In detail, by observing the τ- s relationship from a group of pullout tests, the mean and variance of τ- s curves are able to be determined. Secondly, by discretizing the mean and variance of τ- s curves into ith intervals, corresponded points (si,μ(τio)) and (si,v(τio)) are selected. At last, by calculating the corresponded points (si,μ(τi)) and (si,v(τi)) from the proposed model, an objective function R is able to be constructed, which is illustrated as follows:
R μ = i = 1 N [ μ ( τ i ) μ ( τ i o ) ] 2 a p p r o a c h m i n
R V 2 = i = 1 N [ V 2 ( τ i ) V 2 ( τ i o ) ] 2 a p p r o a c h m i n
Specifically, an optimization criterion is introduced for determining the applicability of the objective function R, that is, the value of R is approaching minimum. The identification flow chart is listed in Figure 5.

4. Model Validation and Discussion

4.1. Bond Stress-Slip Behaviors in the Sense of Mean Value

To validate the effectiveness of the proposed model for the mean τ- s curve in this work, the experimental results of bond stress-slip relationship under monotonic loading [25,26] were used in this section. The parameters of the proposed model were determined by using the above-mentioned random identification process (i.e., a genetic algorithm), and the results were shown in Table 1 and Figure 6, Table 2 and Figure 7, respectively. It is observed from Figure 6 and Figure 7 that the predicted results agrees well with the experimental results.

4.2. Stochastic Bond Stress-Slip Behaviors

In order to verify the effectiveness of the proposed model for characterizing the stochastic behaviors, the experimental results of a series of pull-out tests [23,27,28] were examined in this section (see Figure 8).
By using the proposed identification process, the parameters and predicted results were shown in Table 3 and Table 4 and Figure 9 and Figure 10.
Specifically, Figure 10 illustrates the effectiveness of the proposed model for characterizing the stochastic mechanical behaviors of the interface. It shows that the proposed model enable us to capture the random nature of the interface by the derived expression of the variance, based on the micro mechanics and random variable theory [22,24]. In detail, the proposed model is able to describe both the mean and variance value of the bond stress-slip behaviors (see Figure 6 and Figure 7, Figure 9 and Figure 10) and promotes the development and application of statistical mechanics and random variable theory by deriving the expression of variance in Section 2, Section 3, Section 4 and Section 5. Moreover, the literature [24] has considered that this expression is very difficult to be obtained and further study is needed in the future, due to the complexity of related theoretical analysis.

4.3. The Comparison between the Proposed Model with a Non-Stochastic Model

A comparison between the predicted mean τ–ѕ curve by using the proposed model and a non-stochastic model obtained by literature [27] is listed in Figure 11. The observation shows that the predicted τ–s curve generated by proposed model was closer to the experimental results, while an underestimation of the stress in most load region is found from the prediction obtained by the non-stochastic model. Additionally, by considering a lack of the capability of characterizing the stochastic behaviors (variance), the applicability of the proposed model is more widely to some extent compared with the non-stochastic model mentioned in this work.

5. Conclusions

In this paper, a stochastic damage model for bond stress-slip relationship of rebar-concrete interface was proposed based on the micro mechanical method and random variable theory. The conclusions can be drawn as follows:
A microscopic damage model was proposed to mimic the micro-damage mechanical behaviors of rebar-concrete interface under monotonic loading. The model consisted of a spring element and a friction element in parallel, and the failure/friction behavior of the elements was proposed to be controlled by a switching function. When the slip of the spring element reached a certain threshold, failures occurred and resulted in the friction element slipping instantly. In spite of the conventional theoretical model, by setting the fracture threshold of the spring element as a random variable, the stochastic properties were able to be depicted. More precisely, based on the introduced methodology in this work, the derivation of the expressions for the mean and variance of the bond stress-slip relationship were obtained through the help of statistical mechanics.
A search heuristic global optimization algorithm (i.e., genetic algorithm) was adopted to identify the random variable parameters evolved in the proposed model based on the experimental results, and the lognormal distribution was utilized.
The validation of the proposed model was performed by comparison between the predictions and the experimental results. It has been revealed that the proposed model can effectively describe the stochastic constitutive relationship of the bond stress-slip of rebar-concrete interface under monotonic loading.
This work may be of significance for studying and modelling the stochastic constitutive relationship of bond stress-slip relationship of rebar-concrete interface under uniaxial monotonic loading, as well as providing a better understanding and awareness of its uncertain effects on relevant engineering applications.

Author Contributions

Conceptualization, X.L., Z.Y., Z.S. and J.Y.; methodology, X.L., Z.Y., Z.S. and J.Y.; software, X.L., Z.Y., Z.S. and J.Y.; validation, X.L., Z.Y., Z.S. and J.Y.; formal analysis, X.L., Z.Y., Z.S. and J.Y.; investigation, X.L., Z.Y., Z.S. and J.Y.; visualization, X.L., Z.Y., Z.S. and J.Y.; resources, X.L., Z.Y., Z.S. and J.Y.; data curation, X.L., Z.Y., Z.S. and J.Y.; writing—original draft preparation, X.L., Z.Y., Z.S. and J.Y.; writing—review and editing, X.L., Z.Y., Z.S. and J.Y.; supervision, Z.S. and Z.Y.; project administration, Z.S. and Z.Y.; funding acquisition, Z.S. and Z.Y.

Funding

This research was supported by the National Key R&D Program of China, grant number 2018YFD1100401; the National Natural Science Foundation of China, grant numbers 51808558, 51820105014, U1434204, 51378506, and 51478478; and the Natural Science Foundation of Hunan Province, China, grant number 2019JJ50800.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eligehausen, R.; Popov, E.P.; Bertero, V.V. Local Bond Stress-Slip Relationships of Deformed Bars Under Generalized Excitations; no.UCB/EERC-83/23; University of California, Berkeley. Earthquake Engineering Research Center: Berkeley, CA, USA, 1983. [Google Scholar]
  2. Ragueneau, F.; Dominguez, N.; Ibrahimbegović, A. Thermodynamic-based interface model for cohesive brittle materials: Application to bond slip in RC structures. Comput. Methods Appl. Mech. Eng. 2006, 195, 7249–7263. [Google Scholar]
  3. Alsiwat, J.M.; Saatcioglu, M. Reinforcement Anchorage Slip under Monotonic Loading. J. Struct. Eng. 1992, 118, 2421–2438. [Google Scholar] [CrossRef]
  4. Harajli, M.H.; Salloukh, K.A. Effect of fibers on development/splice strength of reinforcing bars in tension. ACI Mater. J. 1997, 94, 317–324. [Google Scholar]
  5. Xu, Y.L. Simplified model of bond-slip constitutive relationship of reinforced concrete. Eng. Mech. 1997, 2, 34–38. [Google Scholar]
  6. Wang, B.; Bai, G.L.; Dai, H.J. Experimental and mechanical analysis of bond-slip performance between recycled concrete and rebar. Eng. Mech. 2013, 30, 54–64. [Google Scholar]
  7. Lutz, L.A.; Gergely, P. Mechanics of bond and slip of deformed bars in concrete. Am. Concr. INST J. Proc. 1967, 64, 711–721. [Google Scholar]
  8. Nilson, A.H. Nonlinear analysis of reinforced concrete by the finite element method. ACI J. Proc. 1968, 65, 757–766. [Google Scholar]
  9. Mirza, S.M.; Honde, J. Study of bond stress-Slip Relationship in Reinfored Concrete. ACI J. Proc. 1979, 76, 19–46. [Google Scholar]
  10. Kankam, C.K. Relationship of Bond Stress, Steel Stress, and Slip in Reinforced Concrete. J. Struct. Eng. 1997, 123, 79–85. [Google Scholar]
  11. Di, S.L. Experimental study on the relationship between grip force–slip of reinforced concrete. Ph.D. Thesis, Southeast University, Jiangsu, China, 1981. [Google Scholar]
  12. Jin, Z.S.; Zhu, W.F.; Pang, T.H. Study on bonding properties of rebars and concrete. J. Nanjing Inst. Technol. 1985, 2, 73–85. [Google Scholar]
  13. Teng, Z.M.; Zhang, H.G. Experimental study on splitting bond failure and extension length of rebars in reinforced concrete beams. J. Civ. Eng. 1989, 2, 33–43. [Google Scholar]
  14. Somayaji, S.; Shah, S.P. Bond stress versus slip relationship and cracking response of tension members. J. Am. Concr. Inst. 1981, 78, 217–225. [Google Scholar]
  15. Yankelevsky, D.Z. Bond action between concrete and a deformed bar–A new model. J. Am. Concr. Inst. 1985, 82, 154–161. [Google Scholar]
  16. Zhao, W.P.; Xiao, J.Z. On bond-slip constitutive model between ribbed steel bars and concrete. Eng. Mech. 2011, 8, 164–171. [Google Scholar]
  17. Song, Y.P.; Zhao, G.F. Study on bond-slip performance between rebars and concrete. J. Dalian Inst. Technol. 1987, 26, 93–100. [Google Scholar]
  18. Ruiz, M.; Muttoni, A.; Gambarova, P.G. Analytical modeling of the pre- and post-yield behaviour of bond in reinforced concrete. J. Struct. Eng. 2007, 133, 1364–1372. [Google Scholar] [CrossRef]
  19. Alfano, G.; Marfia, S.; Sacco, E. A cohesive damage–friction interface model accounting for water pressure on crack propagation. Comput. Methods Appl. Mech. Eng. 2006, 196, 192–209. [Google Scholar] [CrossRef]
  20. Alfano, G.; Sacco, E. Combining interface damage and friction in a cohesive-zone model. Int. J. Numer. Methods Eng. 2006, 68, 542–582. [Google Scholar] [CrossRef]
  21. Soh, C.K.; Liu, Y.; Dong, Y.X.; Lu, X.Z. Damage Model Based Reinforced-Concrete Element. J. Mater. Civ. Eng. 2003, 15, 371–380. [Google Scholar] [CrossRef]
  22. Li, J.; Wu, J.Y.; Chen, J.B. Random Damage Mechanics of Concrete; Science Press: Beijing, China, 2015. [Google Scholar]
  23. Shan, Z. Stochastic damage model for concrete and its application. Ph.D. Thesis, Central South University, Changsha, China, 2017. [Google Scholar]
  24. Feng, D.; Ren, X.; Li, J. Stochastic damage hysteretic model for concrete based on micromechanical approach. Int. J. Non-linear Mech. 2016, 83, 15–25. [Google Scholar] [CrossRef]
  25. Campione, G.; Cucchiara, C.; La Mendola, L.; Papia, M. Steel–concrete bond in lightweight fiber reinforced concrete under monotonic and cyclic actions. Eng. Struct. 2005, 27, 881–890. [Google Scholar] [CrossRef]
  26. Oh, B.H.; Kim, S.H. Realistic Models for Local Bond Stress-Slip of Reinforced Concrete under Repeated Loading. J. Struct. Eng. 2007, 133, 216–224. [Google Scholar] [CrossRef]
  27. Sun, M. Experimental study and finite element analysis on bond-slip constitutive behavior of reinforced concrete. Master’s Thesis, Zhejiang University, Hangzhou, China, 2015. [Google Scholar]
  28. Wu, L.Y. Research on seismic behavior of precast concrete columns considering micro bond-slip damage. Master’s Thesis, Central South University, Changsha, China, 2017. [Google Scholar]
Figure 1. Microscopic element model.
Figure 1. Microscopic element model.
Materials 12 03151 g001
Figure 2. Sketch of two statuses of individual microscopic element under loading.
Figure 2. Sketch of two statuses of individual microscopic element under loading.
Materials 12 03151 g002
Figure 3. Sketch of parallel elements system.
Figure 3. Sketch of parallel elements system.
Materials 12 03151 g003
Figure 4. Sketch of the fracture threshold.
Figure 4. Sketch of the fracture threshold.
Materials 12 03151 g004
Figure 5. Flow chart of algorithm for parameter identification of random variable. (a) Main program module of genetic algorithm; (b) Main program module of the objective function
Figure 5. Flow chart of algorithm for parameter identification of random variable. (a) Main program module of genetic algorithm; (b) Main program module of the objective function
Materials 12 03151 g005
Figure 6. Bond stress-slip behavior in the sense of mean value. Experimental results were obtained from literature [25]. (a) Comparison between experimental and predicted results; (b) The probability density function; (c) Damage evolution; (d) The graph of objective function R changing with optimization algebra.
Figure 6. Bond stress-slip behavior in the sense of mean value. Experimental results were obtained from literature [25]. (a) Comparison between experimental and predicted results; (b) The probability density function; (c) Damage evolution; (d) The graph of objective function R changing with optimization algebra.
Materials 12 03151 g006
Figure 7. Bond stress-slip behavior in the sense of mean value. Experimental results were obtained from literature [26]. (a) Comparison between experimental and predicted results; (b) The probability density function; (c) Damage evolution; (d) The graph of objective function R changing with optimization algebra.
Figure 7. Bond stress-slip behavior in the sense of mean value. Experimental results were obtained from literature [26]. (a) Comparison between experimental and predicted results; (b) The probability density function; (c) Damage evolution; (d) The graph of objective function R changing with optimization algebra.
Materials 12 03151 g007aMaterials 12 03151 g007b
Figure 8. Experimental results of pullout tests [23,24,25].
Figure 8. Experimental results of pullout tests [23,24,25].
Materials 12 03151 g008
Figure 9. Bond stress-slip behavior in the sense of mean value. Experimental results were obtained from literature [23,27,28]. (a) Comparison between experimental and predicted results; (b) The probability density function; (c) Damage evolution; (d) The graph of objective function R changing with optimization algebra.
Figure 9. Bond stress-slip behavior in the sense of mean value. Experimental results were obtained from literature [23,27,28]. (a) Comparison between experimental and predicted results; (b) The probability density function; (c) Damage evolution; (d) The graph of objective function R changing with optimization algebra.
Materials 12 03151 g009
Figure 10. Stochastic bond stress-slip behaviors. (a) Comparison between predicted and experimental results for standard deviation value; (b) Curves of predicted stochastic behaviors; (c) Comparison between predicted and experimental results; (d) Cloud map of predicted stochastic behaviors.
Figure 10. Stochastic bond stress-slip behaviors. (a) Comparison between predicted and experimental results for standard deviation value; (b) Curves of predicted stochastic behaviors; (c) Comparison between predicted and experimental results; (d) Cloud map of predicted stochastic behaviors.
Materials 12 03151 g010aMaterials 12 03151 g010b
Figure 11. The comparison between the proposed model and a non-stochastic model.
Figure 11. The comparison between the proposed model and a non-stochastic model.
Materials 12 03151 g011
Table 1. Identification results of random variable parameters. Experimental results were obtained from literature [26].
Table 1. Identification results of random variable parameters. Experimental results were obtained from literature [26].
λζβG(MPa)R
6.111.080.070.0912.86
Table 2. Identification results of random variable parameters. Experimental results were obtained from literature [26].
Table 2. Identification results of random variable parameters. Experimental results were obtained from literature [26].
λζβG(MPa)R
0.8720.50.05615.2528.82
Table 3. Identification results of random variable parameters of the mean value.
Table 3. Identification results of random variable parameters of the mean value.
λζβG(MPa)R
5.0761.3280.1020.0992.03
Table 4. Identification results of random variable parameters of the variance value.
Table 4. Identification results of random variable parameters of the variance value.
λ1.λ2ζ1ζ2ρR
5.0494.961.2921.3580.0223.3

Share and Cite

MDPI and ACS Style

Lv, X.; Yu, Z.; Shan, Z.; Yuan, J. A Stochastic Damage Model for Bond Stress-Slip Relationship of Rebar-Concrete Interface under Monotonic Loading. Materials 2019, 12, 3151. https://doi.org/10.3390/ma12193151

AMA Style

Lv X, Yu Z, Shan Z, Yuan J. A Stochastic Damage Model for Bond Stress-Slip Relationship of Rebar-Concrete Interface under Monotonic Loading. Materials. 2019; 12(19):3151. https://doi.org/10.3390/ma12193151

Chicago/Turabian Style

Lv, Xiaoyong, Zhiwu Yu, Zhi Shan, and Ju Yuan. 2019. "A Stochastic Damage Model for Bond Stress-Slip Relationship of Rebar-Concrete Interface under Monotonic Loading" Materials 12, no. 19: 3151. https://doi.org/10.3390/ma12193151

APA Style

Lv, X., Yu, Z., Shan, Z., & Yuan, J. (2019). A Stochastic Damage Model for Bond Stress-Slip Relationship of Rebar-Concrete Interface under Monotonic Loading. Materials, 12(19), 3151. https://doi.org/10.3390/ma12193151

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