#### Abstract

This paper aims to theoretically and numerically assess the effect of setting artificial-induced joints (AIJs) during construction period of a mass concrete structure to release the early-stage thermal stress. With respect to the coupling influences of various factors such as size and boundary of AIJs, an analytical model for its cracking strength on the setting section of mass concrete is proposed based on double-parameter fracture theory. A kind of hyper-finite element analysis (FEA) for many array AIJs in simplified plane pate is also presented by using bilinear cohesive force distribution. The results from the present model and numerical simulation were compared to those of experimental data to prove the efficiency and accuracy of the analytical model and FEA. The model presented in this study for the cracking strength of AIJs provides a simple useful tool to accurately evaluate how many early stress AIJs reduced. The theoretical solution and FEA results could also be significantly contributed to find the “just” and “perfect” release of the temperature stress and to improve the design level of AIJs in mass concrete structure.

#### 1. Introduction

In the past decades, mass concrete structures such as super long underground sidewall, large basement slab, and roller compacted concrete (RCC) have been widely used in the practical construction like subway station [1], high-rise building [2], and RCC arch dam [3]. The early irregular cracks caused by the thermal stress of mass concrete during construction period are a worldwide problem to the structure [4, 5], and how to effectively control the early-age temperature cracks in the mass concrete structure is significantly concerned by researchers [6–8]. Up to now, many methods have been introduced and applied to eliminate these early random cracks, such as adopting deformation joints [9], setting post-poured strips [2], and using high performance materials [10–13]. Although the above methods to some extent can inhibit the growth of random early-age cracks, there are still problems during the realistic applications [2, 10, 14].

Recently, therefore, more attention has been paid to a new method of setting artificial-induced joints (AIJs) in construction period to reduce the early-age thermal stress of mass concrete [1, 15, 16]. Figure 1 illustrates the joint treatment by this new method for different types of mass concrete (i.e., super long underground sidewall, large basement slab, and RCC arch dam). As can be seen from Figure 1(c), the cross section of mass concrete consists of many embedded gaps phase-to-phase arrayed if the setting of AIJs in RCC is taken as an example in the present study. The size of each embedded gap is ; the vertical and horizontal distance between two adjacent embedded gap centres are and , respectively. Cross section embedded AIJs are expected to be cracked prior to avoiding random cracks in other parts/sections of the structures when the early temperature rises in mass concrete due to the heat of cement hydration (Figure 1(c)). Then backfill grouting is applied in the cracking section when the chemical reaction of hydration heat has been done during the construction period of mass concrete structures.

**(a) AIJs in super long underground sidewall**

**(b) AIJs in large basement slab**

**(c) AIJs in RCC**

The cracking strength of AIJs which is related to the size, the shape of each embedded gap, the distance between two cross section embedded AIJs, and so forth is the most key and important parameter in the study of AIJs in mass concrete. This is because only the right value of cracking strength can be expected to avoid random cracks in other parts/sections and thereby to “just” and “perfectly” release temperature stress in mass concrete. If the value of cracking strength is too small or large, the cross section embedded AIJs cannot be cracked firstly or cannot completely reduce the early stress of other parts/sections. In this situation, the setting of AIJs to inhibit the growth of random early-age cracks will to some extent be limited. Therefore it is of great importance to accurately predict the cracking strength of AIJs in evaluating the effect of controlling the early stress in mass concrete structures. Up to now, the accurately analytical solution of cracking strength to “just” and “perfectly” release temperature stress in mass concrete is still not yet available although some experimental investigation, theoretical model, and numerical simulation have been carried out [15–17]. This is because there are so many complex factors such as boundary of cross section, different size, shape of AIJs, concrete strength, and construction schedule which affect the cracking strength [16, 17]. Therefore, it is necessary and of great importance to illustrate the highly precise model and accurate FEA for the prediction of cracking strength of AIJs so as to be able to evaluate its effect on releasing the early stress [18–20].

In this paper, an analytical model of the cracking strength of AIJs in RCC is studied and improved after considering the various factors such as size and boundary of AIJs. Then a kind of hyper-finite element method along with the experimental data is present. The results from the present model and numerical simulation are discussed and further verified. The model presented in this study for the cracking strength is a simple and useful tool to accurately evaluate the reduction of early-age stress. The theoretical solution and FEA results could also be contributed to find the “just” and “perfect” release of the temperature stress and to improve the design level of AIJs in mass concrete structure.

#### 2. Analytical Model of Cracking Strength for AIJs

The analytical solution of cracking strength for AIJs in Figure 1(c) can be simplified as a problem of penetrated crack in an infinite plate due to the fact that the length of arch direction is often hundreds of meters (Figure 2). The vertical distance between two adjacent embedded gap centres in Figure 2 is . Similar to the development of other cracks in concretes, the AIJs also include cracking initiation, stable growth, and unstable propagation process for mass concrete. As assumed, is the tensile strength along the arch direction, and and represent initial and unstable cracking strength, respectively. The propagation of AIJs occurred when thermal stress away from the AIJs is up to , thereby weakening the cross sections. Similarly, AIJs start to be unstably fractured as thermal stress comes to unstable cracking strength . Apparently, one can find that there is a stable propagation stage of the cracking of AIJs between and .

As mentioned above, the values of and are influenced by many factors [16, 17]. In the present study, the coupling effect of these factors such as boundary of cross section, different size, and concrete strength on the cracking strength of AIJs can be expressed as the following equations [16, 17]: where is the ratio of initiation crack length to depth and is the fracture energy of concrete materials. and represent the tensile strength and elastic modulus of RCC, respectively. and are the functions of and can be varied from different types of structure.

Although the fracture energy in (3) can be generally used to describe the average amount of energy consumed during the breaking down process of AIJs in mass concrete, it cannot effectively express the respective consumption of energy during two important stages of cracking expansion, namely, the stable growth and unstable propagation process [21]. To characterize the energy release ratio in the above different crack expanding periods, is divided into two fracture parameters: initiation fracture energy and unstable fracture energy . is associated with cracking initiation stress and initial crack length , while is related to maximum stress and critical effective crack length . More information about the parameters is shown in Table 1.

The solution of can be determined by experimental test of the three-point bending beams [21]:where and are the specimen thickness and loading span, respectively. is the ratio of the critical effective crack length to beam depth , is the maximum load, and is a coefficient related to :

If the contribution of energy due to cohesive forces is defined as , the initiation fracture energy can also be given aswhere can be obtained by integration method of cohesive forces on the crack surface.

According to fracture theory, the softening models of concrete [21] can be simplified into the bilinear cohesive force model to describe the cracking behaviors of mass concrete, as illustrated in Figure 3. The expression is In (7), is the effective crack propagation length,* w* is the crack opening displacement, and is the critical value of . can be regarded as crevice genesis zone, which can be evaluated by the following expression:

The stress distribution of the cross section where AIJs embedded will be more affected by the plasticity of concrete materials with the lager values of the fracture energy of RCC. As a result, the initiation cracking strength and the unstable cracking strength of AIJs in mass concrete can be determined by the above values of and in (4) and (6). The initiation cracking strength can be given as

After the stable growth and unstable propagation stage of AIJs, the critical value of unstable cracking strength can also be written as

#### 3. The Hyper-Finite Element Analysis for AIJs

Among the fracture models of concrete materials proposed by researchers, the fictitious crack model is famous and widely used. It is especially suitable for FEM analysis of AIJs in mass concrete. However, there still exist many difficulties to obtain accurate stress distribution on the vicinity of the crack tip even though a refined mesh is used [21]. An effective way to solve the problem in this study is to apply a hyperelement on the vicinity of AIJs tip based on the Hamiltonian theory of elasticity. Thereafter, the cracking strength of FEM analysis for AIJs in mass concrete is described when the hyperelement is combined with general finite elements.

For bilinear cohesive force distribution in fictitious crack model shown in Figure 4, the propagation process of AIJs in mass concrete can be considered as propagation of elastic cracks, which is a function of a cohesive distribution close force along the propagation length . The region of cross section embedded AIJs is dispersed with a ring hyperelement surrounding the crack tip and also ordinary finite elements around the hyperelement, as is illustrated in Figure 4. The equation of variation in the form of Hamiltonian for the crack-ring-singular element loaded by the linear cohesive force can be expressed aswhere and are the radial and tangential displacement, respectively. is its displacement vector, , , and is the dual vector.

The parameters , , are obtained bywhere and .

Assuming that there are nodes on the hyperelement and the unknown variables of each node are displacements of except for the node at which has only one variable , so the hyperelement has degrees of freedom. If is Eigen value and is Eigen-function vector, the complete state function vector expansion is expressed aswhere are unknown generalized constants and is a special solution:

If and are the generalized constant vector and the nodal displacement vector, respectively, the relationship between and iswhere = , = , = ,, and is the value vector of special solution at the nodes.

The components of the inverse of matrix arewherewhere are tangential coordinates of node .

The solution can be carried out by iterating over (13) and (11). and load vector can be obtained by the integration of partial stiffness matrix of hyperelement:where and are

#### 4. Efficiency and Accuracy of the Present Method

To obtain parameters of fracture energy ( and ) for AIJs of mass concrete, two types of beams with sizes of mm and mm were prepared for three-point bending measurement. The numerical results and laboratory measurement were then compared to verify the analytical model developed in this study, as well as investigate the accuracy of FEM analysis for prediction. According to the experimental data of RCC and (4)-(5), the values for and are shown in Figure 5. It can be seen that the evaluated values of and are independent of both size and boundary. The average values of and for both series are 15.32 N/m and 64.82 N/m, respectively. And these values are applied to calculate the cracking strength for AIJs in the following discussions.

**(a)**Values for (N/m)

**(b)**Values for (N/m)##### 4.1. FEM Simulation of Three-Point Bending Beams

For the FEM simulation of three-point bending beams, the tangential stress distribution near the crack tip is illustrated in Figure 6. Here, the critical crack length for stable propagation and unstable propagation is 30 mm and 50 mm, respectively. The numerical results of unstable fracture energy for three-point bending beams, as well as its experimental results, are present in Table 2. It can be seen that the prediction given by the present FEM is in a good agreement with the experimental data, which demonstrates that the FEM analysis is efficient and accurate for prediction.

**(a)**

**(b)**

##### 4.2. Results of Cracking Strength for AIJs in Mass Concrete

Whether AIJs in the early period of mass concrete will induce crack or not is important because there is no longer effect if setting of AIJs on the cross section is carried out after thermal stress is released. As a result, the “just” and “perfectly” cracking strength of AIJs in depositing concrete construction period is a key parameter along with the increase of altitude of mass concrete like RCC. Taking into account representative altitude of RCC, the cross section embedded AIJs of mass concrete can be considered as penetrated cracks in infinite plate for the simplification of hyper-finite element analysis (Figure 2). The typical local mesh for artificial-induced joints is shown in Figure 7.

The cracking strength for AIJs in mass concrete predicted by the present model is listed in Table 3. It can be found from Table 3 that the results from FEA are in consistence with those from analytical models, indicating that the proposed model can be used to accurately assess the effects of AIJs on releasing the early stress in mass concrete. One can also see when the relative cracking strength of AIJs (i.e., and for m, m to 0.9 m) reaches about 30%–36% and 55%–58% separately, the AIJs will be initiated to propagate and to be in unstable propagation process, respectively. While for m, m to 1.5 m, the above two values of relative cracking strength slightly decrease to about 23%–28% and 43%–55%. These results can be considered as the basic findings to assess the reduction of early stress of AIJs in mass concrete structure.

#### 5. Conclusions

The effect of setting artificial-induced joints (AIJs) during the construction period of mass concrete structures to release the early-stage thermal stress is theoretically and numerically assessed in this study. An analytical model which is related to the cracking strength of AIJs on the setting section of mass concrete, along with its hyper-finite element analysis in simplified plane pate, is also proposed based on the consideration on coupling influence of various factors such as size and boundary effect. The numerical results, together with the experimental data, prove that the analytical model suggested in this study is efficient and accurate in assessing the cracking strength of AIJs. The basic findings can also be used to describe the effects of AIJs on “just” and “perfect” releasing the early-age stress as well as further improve the design level of AIJs in the mass concrete.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This research was supported by the National Natural Science Foundation of China (51278304), the National Basic Research Program (973 Program) of China (2011-CB013604), and the Technology Research and Development Program (basic research project) of Shenzhen, China (JCYJ20120613174456685 and JCYJ20130329143859418). The authors greatly appreciate these supports.