3.1. Validation of Models
Figure 4 shows the load vs. vertical displacement of the reinforced beam span obtained through numerical simulation and the values obtained experimentally (3 reinforced beams, entitled R1, R2, and R3) by Luca and Marano [
14]. The curve obtained via numerical simulation was plotted to the maximum point where the convergence of the model occurred with a tolerance of 10
−5. The values of the applied force for displacement in the service limit state condition (L/200—L is the useful distance between the end of the static bending test, 1800 mm—
Figure 1) established by the Brazilian standard NBR 7190 [
19] and for the rupture values of the reinforced beams are shown in
Table 6, and it should be noted that T1, T2, and T3 are the beams with no reinforcement (reference), also tested by Luca and Marano [
14]. The ductility, which is the mechanical behavior expected of timber beams by incorporating steel bars, was more noticeable in beam R3. Thus, the results of the load vs. displacement curve of beam R3 were used for comparison with the results obtained from the simulation (
Table 6).
Considerable agreement can be seen in
Figure 4 and
Table 6 between the numerical and experimental results, not only in the linear phase but also in the nonlinear phase when the plasticization mechanism of the materials starts. The maximum experimental and numerical load difference was not more than 8.28%. The behavior of the load vs. displacement curve obtained by the numerical model presented a behavior similar to the curves of the experimental values. In addition, the post-rupture ductile behavior was similar to the one Luca and Marano [
12] observed, as the cross-section of the GLT beam showed residual resistance to the loads. The differences observed between the numerical and experimental curves may be due to the non-consideration of imperfections, such as knots and distortions in the fibers. In addition, the timber properties were defined from empirical relationships due to the lack of experimental values.
In the experimental models, the non-reinforced beams failed at the end of the tensile region, causing a sudden loss of resistant section, and the reinforced beams presented a mixed failure of tensile and compression, with cracks in the tensile area and crushing in the compressed one.
Figure 5a,b show the stresses obtained by numerical analysis of the models with no reinforcement and reinforcement, respectively. Considerable plastification of the beam is observed, in which the central region of the span presented stresses equal to or greater than the compressive strength parallel to the fibers (
Table 1).
Reinforced beams initially fail in the tensile zone and then crush in the compression zone [
14]. Several authors observed such behaviors [
2,
4,
6,
13,
16,
39,
40,
41,
42]. Thus, the stresses presented by the numerical model are consistent with the stresses obtained from the experimental model.
Figure 6a,b show the stresses in the reinforcement and the maximum vertical displacements of the reinforced beam for the maximum load.
The results show that the proposed numerical models were able to simulate the mechanical behavior of reinforced and unreinforced beams with steel bars in the linear and nonlinear phases. The maximum numerical displacements in the service phase (L (span)/200) and the maximum displacements presented values close to the experimental ones and the respective loads. In addition, the models showed stress concentrations in the experimental rupture regions, demonstrating the model’s agreement.
3.2. Parametric Study Results
Figure 7 shows the Tukey test results of the factors (variables considered independent)
L,
%s,
h,
b,
d, and
δs over the response variables (considered dependent)
Fs,
Fu, and
δu, respectively. It must be noted that the variables b (base of the beam cross-section) and
d (diameter of the bars) were properly categorized into 5 different classes with an amplitude of 50 mm for variable b and 10 mm for variable
d.
The reinforcement showed an increase of approximately 32% and 49% in the service force values for reinforcement rates of 2% and 4%, respectively, which are lower than the values obtained by Soriano et al. [
4], with an increase in the load service of 53.1% and 79% for reinforcement rates of 2% and 4%, respectively. The values of the modulus of elasticity of the wood and the steel reinforcement used by Soriano et al. [
4] were different from those defined in this work; however, both studies showed significant increases in the load capacity of the beam. It can be observed that the different reinforcement rates showed statistically equivalent mean service load values, both values being higher than those with no reinforcement, indicating the efficiency of the use of steel bars. Thus, reinforced elements allow the solution of certain load and height limits of the cross-section related to design restrictions. In addition, it was noted that the beam height had a significant influence on the service load, with the heights of 300 and 350 mm showing statistically equivalent mean values. However, the average supported by the 350 mm high beams was substantially higher than the 300 mm beams, resulting in a value of approximately 44%.
The reinforcement values of 2% and 4% presented average ultimate load values of 42.9% and 66.9%, respectively, when compared with beams with no reinforcement. The values found in this research were higher than those observed by Luca and Marano [
14], in which the increase in the ultimate load was 48.1% for a reinforcement rate of 0.82%. The results were also higher than the values obtained by Yang et al. [
6] and Chaudhari and Chakrabarti [
20], in which the application of fiberglass bars in the stretched and compressed regions increased the final load by approximately 30%.
The increase observed in the present study was also greater than that observed by Raftery and Kelly [
13] in beams reinforced with basaltic fibers (an increase of 23% in maximum capacity). In addition, the increase in the ultimate load obtained in this research was similar to that presented by Yang et al. [
6], in which the reinforcement with a carbon fiber plate located in the tensile zone showed an increase of approximately 56% in the ultimate load (reinforcement rate of 1.85). Finally, the reinforcement with fiberglass bars showed better performance, increasing the final capacity of the beam by 68% with 1.4% of reinforcement [
16]. The fact that the different reinforcement rates have statistically equivalent averages and service loads shows the efficiency of the use of reinforcement, even with a reinforcement rate of 2%.
The results illustrated indicate that the heights of 300 and 350 mm did not show statistical differences from the ultimate load, and did not significantly increase the value of the ultimate load of the structure. Thus, for the models proposed in this research, the height of 300 mm was the most efficient because it presented a load approximately 4% lower than the 350 mm beam. Base (b) and diameter (d) parameters showed similar behaviors, with categories 1 and 2 showing equivalent average values. As a result, category 1 has more efficient values because it requires less material.
The results obtained show the efficiency of reinforcement with steel bars in the tensile and compression areas, with a significant increase in the ultimate load. In addition, steel bars are easily found, unlike other reinforcements such as synthetic and natural fibers, and they have similar behavior in compression and tensile. All these factors make steel bar reinforcements an efficient alternative to increase the load capacity of GLT timber elements.
Analyzing the last displacements as a function of height, the consideration of 100 mm, 150 mm, or 200 mm in the height value did not promote significant changes in the mean values of δu because they were statistically equivalent to each other. The 250 mm and 300 mm measurements were statistically equivalent to each other but less than 100 mm, 150 mm, and 200 mm. Finally, the height of 350 mm showed the lowest average value.
As expected, beams with a lower height and, consequently, a lower moment of inertia, had the highest values of ultimate displacements. It was also noted that the use of a reinforcement rate of 2% did not significantly modify the average values of δu, being statistically equivalent to each other. In contrast, beams with a reinforcement rate of 4% showed lower mean values of δu. The beams with reinforcement rates showed average displacements lower than 4.5% and 16.6% for reinforcement rates of 2% and 4% when compared to beams with no reinforcement. In addition, reinforcements of categories 4 and 5 in diameter were used in beams with higher heights (300 mm and 350 mm), which justifies the considerably lower mean values of δu.
For a better understanding of the values considered in this research for the measures of the base (
b) and the diameters (
d) of the steel bars, which are of fundamental importance for the generation of the regression models, in
Figure 8, the respective frequency histograms are presented, with M being the mean, SD the standard deviation and N the number of observations.
Figure 8a,b show the frequency of distribution of variables
b and
h considered in the simulations (synthesis of the results generated from the simulations), showing that the two variables (
b and
d) showed normality in their respective distributions.
For the use of the regression models presented below as a way of estimating
Fs, Fu, and
δu, the limits of the variables considered dependent must be respected, being these: 2000 mm ≤ mm
L ≤ 10000 mm, 33.33 mm ≤ mm
b ≤ 233.33 mm; 100 mm ≤ mm
h ≤ 350 mm, 0% ≤
%s ≤ 4%, and 0 mm ≤ mm
d ≤ 45.60 mm. Equations (13)–(16) express the models for estimating
δs, Fs, Fu, and
δu, respectively, together with the adjusted determination coefficients (R
2 adj), and in
Table 7 are presented, for each of the four adjustments, the significance (
p-value < 0.05) or not (
p-value ≥ 0.05) of the models and their component terms, and it should be noted that the closer to 0, the more significant the terms and the model are.
From
Table 7, none of the evaluated models were considered significant by ANOVA. The low value associated with
R2 adj (3.48%) of the regression model in Equation (13) is because the limit displacements (
δs), regardless of other factors, are always equal to 10 mm, 25 mm, or 50 mm (
L/200). Estimates of better precision were achieved with the models in Equations (6) and (7), and of lower precision in estimating the ultimate displacement (Equation (11)). Even though
R2 adj is close to 70% in the models in Equations (14) and (15), these were still considered insignificant by ANOVA, which implies that variations in the independent variables cannot significantly explain the variations suffered by the dependent variables.
Knowing that the displacements in the service condition are based on the useful span of the beams [
19], in an attempt to improve the models previously obtained, the
δs was also considered as an independent variable, which resulted in new adjustments for
Fs, Fu, and
δu, expressed by Equations (17)–(19), respectively, and the results of ANOVA in
Table 8.
Based on the terms considered significant by ANOVA presented in
Table 8 and the hierarchy of the isolated terms, new regression models (hierarchical models) were generated, expressed in Equations (20)–(22).
The large reduction in terms from 21 to 6 in the Equation (17) model compared to the Equation (20) model had little impact on the adjusted determination coefficient, which decreased from 97.57% to 90.24%, showing the robustness of ANOVA as an analysis tool for the sensitivity of the terms of the models, and the same occurred with the models of Equations (13) and (14), which were changed to Equations (13) and (14), respectively. Regarding the final displacement, the only independent variable considered significant by ANOVA was the limit displacement of the service condition.
Some authors have developed analytical and theoretical models to assess the increase in stiffness and the total load of beams reinforced with different materials. This research presented equations to estimate the service and ultimate loads and displacement of beams reinforced with steel bars positioned in the tensile and compression regions. The expressions presented coefficients of determination adjusted above 90% to estimate the last loads and above 70% to determine the ultimate displacement. With simplified expressions (Equations (20)–(22)), it was possible to estimate with good precision the service load and the ultimate load using parameters of the cross-section and the service limit displacement, showing the high applicability of the formulations. The results indicate that the expressions can be used to estimate load values of reinforced beams with parameters within the intervals analyzed in this work.
The results presented contribute to deepening the knowledge about beams reinforced with steel bars, analyzing the influence of different parameters on the load and displacement, and proposing expressions to estimate the loading of reinforced beams. However, more experimental and numerical studies with different reinforcement rates, different strengths of timber and reinforcement bars, and different beam dimensions are needed, aiming at more information to collaborate in the design of the structural project.