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
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,
where
G is the stiffness of the bond-slip relationship.
ii. After the failure of the spring element,
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):
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)
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:
Taking the limit of Equation (6) as
N approaches infinity, one obtains
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:
When
N approaches to infinity, the damage variable is able to be obtained as follows,
By substituting Equation (9) to Equation (7), a damage model for the bond stress-slip relationship is derived such that
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,
τ =
βG ∫
01 Δ(
x)d
x.
Hence, the mean value of the damage variable
d can be expressed as follows:
The square of the mean value of the damage variable is derived such that,
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:
where
P denotes the probability.
Therefore, the mean function of bond stress can be expressed as
2.5. Variance Value of Stochastic Mechanical Responses
The variance function of stochastic mechanical responses of the parallel system is expressed as follows,
By combining Equations (9) and (10), it is able to obtain
E(
τ2) such that,
By substituting Equations (11) and (12) to Equation (18),
E(
τ2) is able to be derived as follows
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,
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
Therefore, the distributions can be obtained respectively as follows,
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:
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.
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.