Study on Finite Element Model of the Prefabricated Reinforced Concrete Beam-Column Joints with Grouted Sleeves

Based on the monotonic tensile test of grouted sleeve specimens conducted, this paper uses multifactor regression analysis to construct the equivalent constitutive relationship of grouted sleeve specimens under uniaxial tension. )e study based on this constitutive relationship of grouted sleeves and the effect of bond and slip between steel and concrete were considered. )e prefabricated reinforced concreted beam-column joints with grouted sleeves were presented with finite element software ABAQUS. )e seismic behavior of prefabricated reinforced concreted beam-column joints with grouted sleeves under lowfrequency cyclic loading was then investigated. In addition, parametric studies via finite element analysis were performed to examine the influence of various parameters on the strength and energy dissipation capacity of the specimens. )e simulation results show that plastic deformation was mainly observed near the beam-column interface; the hysteretic curve of this joint was plump. )e test results showed that good energy dissipation and displacement ductility capacities can be achieved. )e error of yield load between the numerical simulated and experimental result was 7.11%, the error of peak load was 6.88%, the error of ultimate load was 3.76%, and the error of displacement ductility was 7.84%. Results showed that the calculated results obtained in the paper agree well with test results from the references. )e finite element model adopted in this paper can reflect the seismic behavior of the prefabricated reinforced concrete beam-column joints with grouted sleeves by using equivalent constitutive relation.


Introduction
Construction with the precast concrete (PC) method has gathered much attention as a substitution for convention cast-in-situ concrete structures because of its good quality, cost-effectiveness, accelerated construction speed, and lower energy consumption [1][2][3][4]. However, the relatively low reliability, high cost of the connections, and the low seismic performance of PC structures have limited their further application [3]. e seismic performance of the PC structure mainly depends on the connection type, the connection reliability between different precast members, and the arrangement of the precast members and cast-in-place concrete [5,6]. To date, many studies have been conducted on the connection between precast concrete members. e mainly used connections include traditional lap connections, improved lap connections, and grouted sleeve connection. Traditional lap connections, such as mechanical connections, welding connections, and lap connections, have been widely used in cast-in-situ concrete structures. However, even with longer lap lengths than other connection types, the seismic behavior of the traditional lap joints is still relatively low because of the short bond strength [7,8].
Besides, the complexity of connecting steel bars between precast components also hinders the wide application of precast concrete structures. In the last decades, more work has been done to solve the problems of the connections in the precast concrete frame by different researchers [9]. Among all the connections, the grouted sleeve has increased its popularity in the PC industry because of its simplicity and convenience for construction. e grouted sleeve effectively solves the technical problem of joint connection with precast concrete components [10]. e fully grouted sleeve and semigrouted sleeve are the two main types of grouted sleeves, while the semigrouted one is more commonly used in practice [11]. e semigrouted sleeve has two ends: one end with threads and the threaded steel bar from the precast part installed in it, and the other end is plain with a deformed steel bar from the postcast part in the sleeve. After the two parts are placed in the sleeve, high-strength grout will be filled in the sleeve and the precast and postcast parts are connected. e seismic behavior of the precast concrete frame is significantly different from the traditional cast-in-situ reinforced concrete because the precast and postcast concrete parts are connected by the grouted sleeves instead of steel reinforcement. Extensive experimental and numerical work has been done to study the seismic behavior of the grouted sleeve connections.
ere are two widely used methods to simulate the grouted sleeves in the structure: the first one uses the solid element to model all the connecting components [12][13][14], which is accurate but with significantly high computational efforts and the other method assumes the sleeve connection is always reliable and models the sleeve as steel reinforcement [15]; this method is computationally efficient and easy to operate, but it fails to capture the actual force in the sleeve and the influence of the sleeve joint on the whole structure. us, the simulation results are usually not accurate enough. Besides, in the current simulations, the effects of the interface at the fresh concrete and existing concrete is usually ignored, which may also reduce the simulation accuracy.
In this paper, the tensile constitutive relationship of the grouted sleeve connection is firstly constructed using multifactor linear statistical analysis on the pull-out test results. en, this relation is applied to the finite element model of a prefabricated beam-column joint with grouted sleeves, and the influences of the fresh and hardened concrete interface on the joint are also considered. It shows that the simulation fits the experiment results well, so the refined finite element model based on the equivalent grouted sleeve constitutive relation can accurately capture the seismic performance of the beam-to-column joint.

Uniaxial Tensile Test and Equivalent
Constitutive Relation of Grouted Sleeve Reinforcement

Uniaxial Tensile Test Result and Analysis.
e strainstress relationship of the grouted sleeve is firstly received based on the pull-out test results by Liu et al. [16]. 26 pull-out tests have been done to get the tensile constitutive relationship of the grouted sleeves. e sleeves were made by either cast iron or steel and the pull-out tests were carried out by a 100 kN universal testing machine with a constant loading rate at 1 kN/s; a few typical stress-strain curves of the grouted sleeves were selected as shown in Figure 1. It can be found that the tensile behavior of the grouted sleeve joint specimens was similar to steel, with firstly a linear increase stage, followed by a yield plateau, and then reached the ultimate strength and failure. e details of the grouted sleeve joints in the test are listed in Table 1. According to the observations of the test results, the parameters that may affect the constitutive relation of grouted sleeve joints include the diameter of reinforcement d, length of sleeve L 1 , length of grouted segment L 2 , the external diameter of the sleeve D 1 , the inner diameter of the sleeve D 2 , the yield strength of reinforcement, and ultimate strength of reinforcement; details of grouted sleeve are shown in Figure 2. Based on the experimental results, an equivalent constitutive relation of grouted sleeve joint similar to reinforcement is proposed.

e Equivalent Constitutive Relation of Reinforcement
Joint with Grouted Sleeve. As discussed in Section 2.1, the grouted sleeve joint specimens have similar characteristics as the reinforcement; therefore, a trilinear constitutive model is selected to simulate the reinforcement joint with grouted sleeve under uniaxial tension, as shown in Figure 3. In Figure 3, point A and Point C are the equivalent yield point and ultimate load point of the grouted sleeve, respectively; segment AB is the yield plateau of the grouted sleeve.

e Value of Characteristic Point Parameter.
e characteristic points A∼C in the trilinear constitutive model for the grouted sleeve joint were decided based on the statistical analysis on the strain-stress data using a statistical software SPSS; the results are as listed in Tables 2 and 3. e parameters considered in the analysis include anchorage length, yield strain, the diameter of reinforcement, the total length of sleeve, the outer diameter of the sleeve, the inner diameter of the sleeve, and the diameter of reinforcement. Value t in Tables 2 and 3 represents the test statistic regression coefficient, sig. is the associated probability value, the significance of the stress at the characteristic point is the yield strength of the reinforcement f y , and the stress at the characteristic point is an increasing function of the yield strength of reinforcement. e yield strength of reinforcement was selected as an independent variable for multiple regression; the results are shown in Table 4. e goodnessof-fitting square (R 2 ) is 0.668, the statistic F � 19.338, the associated probability value sig. � 0.00, the regression is relatively significant, and the fitting result is obtained as shown in (1) Table 5 shows the fitting results of strain at each characteristic point; similarly, the results of other characteristic points can be obtained as shown in the following equations:  17 20 0.0084 0.0538 * In the joint number: the first number "1" represents steel and "2" represents cast iron, the second number represents the diameter of the steel, the third number represents the anchorage length, and the last number "1∼3" represents the three repetitive specimens; if there is an "s" before the last number, it represents the steel is placed eccentrically in the sleeve. For example, "1-16-100-1" represents the first steel sleeve specimen with a diameter of 16 mm and an anchorage length of 100 mm.  Figure 3: e trilinear constitutive model of the grouted sleeve.

Details of the Test.
A fine finite element model based on the grouted sleeve joint pull-out tests was established to simulate the low-frequency cyclic test on a beam-to-column connection, as shown in Figures 4 and 5. In the tests, the precast beam and column ends were jointed by the grouted sleeve and then connected as a whole using postcast concrete. e size and reinforcement details of the specimen are shown in Figure 4. e designed axial load ratio was 0.3, the designed strength grade of precast concrete member was C30 (compressive strength � 30 MPa), the designed strength grade of postcast concrete was C40 (compressive strength � 40 MPa), the longitudinal reinforcements of the beam and column were HRB335, and the stirrups were HPB300. e mechanical parameters of the materials in the specimen are listed in Table 6. To ensure the bonding and anchoring strength of the steel reinforcement, the diameter of the sleeve was 40 mm, the length of the sleeve was 159 mm, the length of the grouted segment was 110 mm, the concrete cover thickness of the grouted sleeves was 20 mm, and the concrete cover thickness of the steel reinforcement was 34 mm. e loading condition of this test [17] was as shown in Figure 5, the bottom of the specimen was fixed by a directional hinge, and the horizontal movement of the specimen top was constrained by a rod. Vertical loads were applied by the servo machine at the column ends. Force control was firstly used, following the scheme as shown in Figure 5 until the specimen yield [18] (when the inflection point first emerged in skeleton curve); then, displacement control was adapted, with the yield displacement as the grade-distinction, cyclic loading each grade-distinction twice. e test ended when the load decreased to be below 85% of the peak load. In Figure 5, P y is the yield load, Δ is the loading displacement at the column ends, and Δ y is the yield displacement. e failure mode, deformation capacity, energy dissipation capacity, ductility, and other parameters of beam-to-column joints of postcast concrete with the grouted sleeve were analyzed.

Fine Finite Element Modeling
e concrete damage plastic model in ABAQUS was used for the concrete; and the reinforcement behavior was modeled as bilinear, with the slope of the elastic stage taken as the elastic modulus of steel and the slop of the plastic stage as 0.01, the yield strength and ultimate strength were obtained by the material test. e grouted sleeve constitutive model adopts the equivalent constitutive model proposed in equations (1)∼(5).

Bond-Slip of Reinforced Concrete.
Spring A was used to simulate the bond and slip between the concrete and the longitudinal reinforcement, and the axial spring stiffness was determined by the Eligehausen model [19]. e spring element in the radial direction of reinforcement mainly simulates the dowel action and the vertical displacement of the reinforcement nodes and the concrete nodes inside the stirrup, which has a low influence on the bond-slip between the concrete and longitudinal reinforcement. e maximum elastic modulus of steel and concrete was used as the spring stiffness along the radial direction of the reinforcement; in this paper, the space between the springs along the longitudinal reinforcement was set to be 50 mm.

Simulation of Precast and Postcast Concrete Interface.
Liu and Zhao [20] and Zhang [21] found that the ultimate tensile strength of the concrete with precast and postcast interface was only 80% of the tensile strength of normal concrete cast once as a whole, and the load dropped rapidly after reaching the ultimate load. An ideal elastic-plastic model was assumed as the compressive behavior at the concrete to reinforcement interface, and the ultimate compressive strength is the compressive strength of concrete. us, the normal force-displacement relationship of the joint spring is shown in Figure 6.
A large number of experimental studies have been conducted on the interface shear strength of unconfined precast and postcast concrete; the shear strength at the interface is mainly composed of three parts: internal cohesive force in the concrete, friction force, and shear resistance of reinforcement. Fan [22] established the formula for the interfacial shear strength at the precast and postcast concrete interface after comprehensively considering several factors affecting the shear strength of the interface, as shown in equation (6), and then determined the calculation parameters by regression analysis. e tangential force-displacement relationship of the spring elements can be obtained by substituting the relevant parameters from the test to the equation:

Element Selection.
In this paper, eight-node hexahedral (brick) elements (C3D8R) were used for concrete, twonode linear truss elements (T3D2) were used to model reinforcements [23], the grouted sleeve joints were simulated by nonlinear Spring A, and the tangential and normal connection at the precasting and postcasting concrete interface were simulated by nonlinear spring elements. e fine finite element model is shown in Figure 7. e nonlinear relationship of the springs is reflected by the respective force-displacement (F-D) curves. e F-D curve is calculated by the normal constitutive and tangential constitutive of the precasting and postcasting concrete interface. beam-column joint conforms to the Strut-and-Tie model as shown in Figure 8(a). Concrete crushing firstly occurs at the beam end near the core area, as shown in Figure 8(b). When the load reaches the crack load, due to the uneven distribution of beam ends stiffness, vertical cracks appear at the interface between the end of the grouted sleeve and the beam-to-column joint. After the concrete at the beam end crashed, the force in the beam decreases sharply, and the force in the structure is mainly taken by the steel frame and concrete column; the stress in the core area of the beamcolumn joint no longer follows the Strut-and-Tie model. As the loading increases, the cracking gradually expands to the middle part of the beam. When the load reaches 30 kN, the crack developed through the beam height, and inclined cracks appear in the core area. By the end of the tests, the cracks are mainly distributed at the beam-column interface, the two ends of the grouted sleeve, and the region within 500 mm from the column end. After the concrete at the beam end yield, the stress at the core stirrups increases sharply and reaches the yield stress, as shown in Figure 8(d).   When loading at the peak displacement, most of the core stirrups are yield, and the stress in the beam is relatively small; the structure enters the failure stage, and crushing failure happens in the concrete at the beam ends. e finite element simulation successfully captures the cracking development in the structure, and the failure modes in the model are identical to those in the test, shown in Figure 9.

Bearing Capacity and Ductility.
e skeleton curves of the finite element simulation and the tested beam-to-column joints of the prefabricated concrete structure with grouted sleeve joints are as shown in Figure 10. e numerical simulation result of the ascending segment of the skeleton curve fits well with the experimental result, and the initial stiffness and peak load of the numerical simulation are close to the test results. However, as the plastic damage of the grouted sleeve is not considered in the model at the descending segment of the skeleton curve, the loading capacity from the numerical simulation is smaller than the test results. e yield load, yield displacement, peak load, peak displacement, ultimate load, and ultimate displacement from the numerical model are shown in Table 7. In Table 7, the displacement ductility coefficient of the specimen [24] is  Advances in Civil Engineering calculated as the ratio of the ultimate displacement to the yield displacement, that is, μ � (Δ u /Δ y ). e difference between the finite element simulation result and the test result of the specimen is within 8%; the bearing capacity of the structure can be effectively simulated; the simulation errors of yield load, peak load, and ultimate load are 7.11%, 3.76%, and 4.72%, respectively, and the error of displacement ductility coefficient is 7.84%. So, the simulation can capture the experiments' results quite well.

Load-Displacement Hysteretic Curve.
e hysteresis curves of the finite element model and the tests on the postcast concrete beam-to-column joints with grouted sleeve joints are shown in Figure 11, and the results are shown in Table 8. Comparison of the analytical and experimental results showed that the hysteresis curves obtained from ABAQUS analyses were quite similar. e shape of the hysteresis curves is relatively plump, which shows that the prefabricated beam-column joints with grouted sleeve joints have a high energy consumption capacity.
Due to the fact that the spring element and steel truss element of the simulated sleeve cannot define the constitutive relation under cyclic load in ABAQUS, when applying forward load, the finite element simulation results are 3.76% lower than the test results, the area enclosed by the hysteresis curve is relatively small when the load is reversed, and the load capacity from the finite element simulation is 13    Advances in Civil Engineering hysteresis loop is similar to the test specimen and tends to be arched, it can be seen from the hysteresis loop area ratio in Table 8 that the energy dissipation capacity error of the component is relatively small when approaching the failure.

Strain at the Joint of the Grouted Sleeve.
e strain at the grouted sleeve was measured as shown in Figure 12. Figure 13 shows the strain value at the maximum cyclic load when the grouted sleeve joint is in a tensile state. In Figure 13, δ is the loading displacement, and δ y is the yield displacement. At the initial loading stage, the axial strain of the grouted sleeve and the reinforcement strain at both ends of the grouted sleeve are basically the same. Due to the large stiffness of the grouted sleeve section, the load is transferred from the beam end to the reinforcement through the grouted sleeve section, the steel bars at both ends of the grouted sleeve joints yield, with the increase in loading, the strain at both ends of the grouted sleeve (measuring points C and D, as shown in Figure 12) near the core of the node increases obviously. As a result, the strain increased at the section with a large bending moment of the beam end; however, the grouted sleeve has been kept at the elastic stage, and the stress of the grouted sleeve is also similar to the experimental result.

Parametric Study
Based on the demonstrated accuracy of the developed finite element model, this section presents a parametric investigation [25] as shown in Table 9 to study the effect of different design variables (the diameter of the sleeve; the length of the sleeve; and the length of the grouted segment) on the connection performance. JD-1 model was chosen as the control specimen for the parametric analysis. Each model is briefly described in terms of variables that were adjusted, and the results are discussed. e key feature of comparison focused on the skeleton curves.    10 Advances in Civil Engineering

Effect of the Sleeve External Diameter.
e external diameter of the sleeve of JD-1, JD-2, and JD-3 is 40 mm, 44 mm, and 48 mm; other parameters are the same. Figure 14 shows the skeleton curves of the beam-column joints with a different external diameter of the sleeve. e comparisons of the curves show that the external diameter of the sleeve has a remarkable effect on the skeleton curves. Table 10 indicates that when the external diameter of the sleeve increases from 40 mm to 44 mm, the ultimate bearing capacity is found to increase by 2.39%, in contrast, the ductility factor of displacement decreases by 6.68%; when the external diameter of the sleeve increases from 44 mm to 48 mm, the ultimate bearing capacity is found to decrease by 8.87%, and the ductility factor of displacement decreases by 23.87%. As the external diameter of the sleeve increased, the peak bearing capacity and the stiffness have obvious increase, the ductility factor of displacement is in decline. Within a certain range, the ultimate bearing capacity increases with the increase of the external diameter of the sleeve; when the diameter exceeds a certain value, the ultimate bearing capacity of the joint decreases.

Effect of the Grouted Segment Length.
e lengths of the grouted segment of JD-1, JD-4, and JD-5 are 110 mm, 100 mm, and 130 mm; other parameters are the same. Figure 15 shows the skeleton curves of the beam-column joints with different lengths of the grouted segment. Figure 15 indicates that the length of the grouted segment has no significant effect on the skeleton curves. e increase in the length of grouted segment leads to a slight increase in the peak bearing capacity, and the ultimate bearing capacity basically has little change. Table 10 indicates that when the length of the grouted segment increases from 100 mm to 110 mm or 130 mm, the ultimate bearing capacity is found to increase by 0.03% and 0.81%, respectively. Additionally, when the length of the grouted segment increases from 100 mm to 110 mm or 130 mm, the ductility of the specimens did not have significant change.

Effect of the Sleeve Length.
e lengths of the sleeve of JD-1, JD-6, and JD-7 are 159 mm, 150 mm, and 143 mm; other parameters are the same. Figure 16 shows the skeleton curves of the beam-column joints with different lengths of the sleeve. Figure 16 indicates that the length of the sleeve has no significant effect on the peak bearing capacity. As the sleeve length increased, the initial stiffness has an obvious decrease. Table 10 indicates that when the length of the sleeve decreases from 159 mm to 150 mm or 143 mm, the ultimate bearing capacity is found to decrease by 0.78% and 6.31%, respectively. Additionally, when the length of sleeve length decreases from 159 mm to 150 mm or 143 mm, Table 10 shows that the length of sleeve length has little influence on the ductility.

Conclusions
In this paper, an equivalent constitutive relation is constructed for the sleeve grouted joints, and the finite element model of prefabricated beam-column joints with a grouted sleeve joint is established using a nonlinear spring element. In addition, the bond and slip between reinforced concrete and the influence of the precast and postcast concrete interface to joint performance are also considered, and the experimental results were compared with the existing research results. e main conclusions are as follows: (1) e equivalent constitutive model under unidirectional tensile load is constructed based on the classification and analysis of the data obtained from the monotonic tensile test of the grouted sleeve joint specimen. e influencing factors of the constitutive model were obtained by statistical analysis software, the stress of the grouted sleeve is significantly correlated with the yield strength of the reinforcement, and the strain is significantly correlated with the yield strain of the reinforcement and the anchorage length of the reinforcement. Besides, the parameters of the relevant influencing factors are chosen. Advances in Civil Engineering 13 (2) e fracture morphology, development process, failure mode of the specimen, and stress of grouted sleeve in the loading process of the finite element model are roughly the same as the experiments, which verifies the equivalent constitutive structure of the grouted sleeve joint and the simulation method. (3) e differences in the yield load, peak load, and ultimate load between the finite element model and test are 7.11%, 3.76%, and 4.72%, respectively. When the structure is close to failure, the skeleton curve begins to fall, which is also similar to the experiment. (4) e hysteretic curve of the finite element model is also similar to that of the test specimens which tends to be arched. e error of the envelope area of hysteretic loop reflecting energy dissipation capacity is also below 13%, in the late loading stage, due to the fact that the plastic damage of the sleeve cannot be defined. (5) In this paper, the equivalent constitutive method based on the nonlinear spring element definition is proposed. Compared with the solid element simulation, this method not only greatly reduces the calculation time of finite element software but also effectively avoids the complicated sleeve modeling and convergence problems, which is meaningful for the grouted sleeve simulation. (6) e finite element parameter analysis results show that the external diameter of the sleeve has a remarkable effect on the seismic performance. As the external diameter of the sleeve increases, the bearing capacity and stiffness increase significantly, and the ductility decreases. Grouted segment length has little effect on the seismic performance. Increasing the length of the sleeve length can improve the bearing capacity and rigidity of the specimen, but the ductility is reduced, and bearing capacity decreases quickly.

Data Availability
e data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.