Predicting Concrete Compressive Strength and Modulus of Rupture Using Different NDT Techniques

1Facultad de Ingenieŕıa Civil, Laboratorio de Materiales, Ciudad Universitaria, Universidad Michoacana de San Nicolás de Hidalgo, Avenida Francisco J. Múgica s/n, Colonia Felicitas del Rı́o, 58040 Morelia, MICH, Mexico 2Facultad de Ingenieŕıa, Ciudad Universitaria, Universidad Autónoma de Querétaro, Cerro de las Campanas s/n, 76010 Santiago de Querétaro, QRO, Mexico 3Instituto Mexicano del Transporte, 76703 Sanfandila, QRO, Mexico 4INFOTEC, Centro de Investigación e Innovación en Tecnologı́as de la Información y Comunicación, Cátedras-CONACyT, Circuito Tecnopolo Sur, No. 112 Colonia Fracc. Tecnopolo Pocitos, 20313 Aguascalientes, AGS, Mexico 5Facultad de Ciencias, Universidad Autónoma del Estado de México, Campus Universitario “El Cerrillo”, El Cerrillo Piedras Blancas, 50200 Toluca, MEX, Mexico


Introduction
In the building industry, the measurement of certain physical and mechanical properties of concrete is particularly important, such as compressive strength (  ) and modulus of rupture (MOR).Test methods to obtain    consist of applying a compressive axial load to molded cylinders or cores.The    of the specimen is calculated by dividing the maximum load reached during the test by the cross-sectional area of the specimen [1]; the method to determine MOR includes the determination of the flexural strength of concrete from simply supported unreinforced concrete beams, where a central load is applied to the beam until fracture (also called a three-point load test: the load applied at the center of the beam and the other two are the supports).The result of such test is reported as the MOR, and might be used to determine compliance with the materials' specifications or as a basis for proportioning, mixing, and placement operations for slabs and pavements [2].Traditionally, these estimates have been obtained following destructive methods, which are the standard ones to determine the mechanical resistance (i.e., compressive, tension, and bending strength) [3].However, these tests are often costly and time consuming.For example, to determine the compressive strength of concrete, it is necessary to process a large amount of testing samples (at least fifteen) [4]; then a compressometer is required in order 2 Advances in Materials Science and Engineering to determine the elastic modulus; these samples are deformed and destroyed in order to determine the compressive strength and the static modulus of elasticity [5].As it can be seen, this procedure modifies the samples, and, consequently, it is impossible to repeat the test in the same sample, in case it would be necessary [6].
On the other hand, it is not easy to determine the mechanical characteristics of the materials using only NDT procedures.This problem is more evident since the models are obtained with laboratory settings and under strict controlled conditions.As a consequence, it is necessary to evaluate the applicability of these models in uncontrolled environments and using different conditions from those used in the formulation of an empirical mechanistic model [7].
One particular characteristic of most of these NDTs is that they can be conducted at the work place and can be used to determine certain properties of already built structures; this process is known as in-place testing [8].
In accordance with ACI 228-1R-03, the objective of inplace testing is to estimate the compressive strength of concrete by the use of NDTs.In order to validate the results, it is necessary to have a known relationship between these results and the mechanical resistance obtained by DTs.This relationship is usually established by performing DTs of drilled cores taken from adjacent locations where the in-place test is conducted.The accuracy of the estimated compressive strength depends on the correlation between the strength of the concrete and the quantity measured by the in-place test.As a general index of concrete strength, the compressive strength of concrete    is important in the performance assessment of existing reinforced concrete structures [9].On the other hand, the MOR is very useful, especially in relation to road slab design and airfield runways, because the flexure tension is a critical factor to determine concrete's performance in such applications [10].
When NDT first appeared for concrete application, it was common to model the prediction of the concrete strength using only one of the above tests.However, as found by many researchers, the combination of different NDTs can lead to more accurate and reliable models [9,[11][12][13].In order to approximate the appropriate estimates of indicators (response variables), for example, Ploix et al. [14] performed a combination or fusion of NDT data as radar test, electrical resistivity and capacity, infrared thermography, impact echo, and ultrasound to predict the simultaneous estimation of the water saturation and porosity rate in concrete for undamaged concretes measured in laboratory conditions.They measured different NDT parameters, such as frequencies, velocities, and attenuations, and after that they conducted an statistical analysis with the purpose of identifying empirical correlations, linking each experimental parameter to both searched indicators.There were identified correlations with bilinear regressions, as a first approximation, but since some different NDT measurements are very sensitive to material heterogeneity, variability, and experimental noise, relative disagreement or conflict between sources of information happened (i.e., measurements); thus the prediction models rapidly became ill-conditioned.Therefore, it was decided to combine several NDT techniques (at least three) providing more suitable information.Kaftandjian et al. [15] used a data fusion approach to improve weld inspection by NDT, combining the use of evidence theory and fuzzy logic, in their framework.X-ray and ultrasonic inspection were modeled in order to generate 3D images from both NDTs for evaluating defect detection and defect sizing of welded joints for industrial purposes, showing that the fusion of data significantly enhances the reliability of the defect recognition.Maierhofer et al. [16] used radar, ultrasonic, impact echo, and thermography data fusion for predicting deterioration of concrete vehicular bridges, noting that the data fusion technique increases reliability, and improved the accuracy of the quantitative results, exploiting the different, and sometimes complementary, physical effects provided for different NDTs.Shah and Ribakov [17] described a testing procedure based on the dynamic response of the cracked or flawed concrete structures subjected to impact loading and on analyzing the ability of the resistivity measurement to detect and to locate cracks and spalling with in situ measurements on a damaged concrete slab.Huang et al. [9] developed a SonReb (the combination of rebound hammer and ultrasonic pulse velocity tests), for compressive strength prediction of concrete samples, using a probabilistic multivariate regression model, taking into account the data of the proportions of the concrete mix (curing conditions and the age of the concrete), assessing the model from Bayesian updating rule.Lui et al. [18] provide a very wide state-of-the-art report regarding NDT data fusion techniques.They established the use of multiple NDT methods to increase the reliability and reduce the uncertainty of testing and evaluation.Nonetheless, there are cases in which the use of all available tests creates a model that only works in the samples used to empirically correlate the obtained data, and thus the model does not generalize.In order to deal with this problem, it is necessary for an automatic method of variable selection (here variables are the NDT) to establish a simple model that achieves the prediction/explanation of the response variable, together with an efficient estimation of the coefficients.
Given that there are only 6 different variables to create the model, in this contribution, a decision was made to test all the linear combinations of these variables, that is, create 64 models.The model proposed was selected from these 64 models by using a cross-validation technique.Furthermore, to illustrate the generality of the model, few samples were removed from the model formulation and selection process, and the results presented on this contribution are computed using the removed samples to validate the predicted model.Clearly, there are situations where it is not possible to compare all possible models.To provide a complete picture of this scenario, it was decided to use a cross-validation technique to select the best from the different 64 models.
In addition to this, in order to compare the proposed model against the ones described in the literature, a comparison was made with more than ten different models against the ones proposed here.The results show that the model proposed here is competitive with the ones proposed in previous investigations.In addition, the independent variables used in this investigation were used also with the other models found in the literature.With this extended set of independent variables, all the possible linear modes were initiated and used the selected mechanism to identify the best model.The result observed was that this novel model is competitive when compared against other models previously reported.However, the proposed model has a better correlation factor and was simpler in the sense that it has fewer degrees of freedom.
In this study, different frameworks to predict the    and the MOR of hydraulic concrete were constructed, based on data from electrical resistivity test (ERT), ultrasonic pulse (UPV), resonance frequency test (RFT), and rebound number test (RNT), obtained from cylindrical and prismatic specimens.
The objective of this research was to predict the    and MOR of concrete, since these are index parameters that determine the quality of the concrete, but with the use of models based on NDT.By employing more than one NDT, it has been possible to generate highly reliable and accurate models, as demonstrated herein.Furthermore, the obtained models can help engineers to understand how each one of the different variables, in the determination of the compressive strength and modulus of rupture, is important, either in laboratory samples or in samples from the place (field samples).This research was conducted in two stages: the first was controlling workability of the concrete, and the second one was controlling water/cement (w/c) ratio.

Electrical Resistivity Test (ERT).
The application of DC current to quantify resistivity was performed by Conard Schlumberger in 1912.It was reported as one of the most successful experimental approach of electrical resistivity survey.In USA, the idea was developed by Frank Wenner in 1915.After that, the method has undergone significant improvement in last three decades [19].
Later Robert W. Fox in 1930 observed electric currents flowing in the copper mines in Cornwall, as a result of chemical reactions within the veins of the deposits [20].In terms of materials and equipment, the ERT can be performed in laboratory specimens, in specimens taken from the analyzed structure, or directly in the structure at the site field [21].ERT measurements are made by applying a current to a porous body, through two current electrodes, and measuring the resulting difference voltage at two potential electrodes.The test must be performed in saturated samples, under the mechanism of electrical conductivity through saturated connected pores, which is based on the conductivity of the water through the pores of the material [22].Liu et al. [23] conducted a study in soil-cement, used to improve the behavior of soft soils in road construction engineering; they built a model to verify the quality control of soil-cement in the field.The model predicted the electrical resistivity (ER) of soil-cement, under certain conditions of curing time and the w/c ratio.It was found that the ER had a good relationship with the    of the material in study.Cao and Chung [24] presented a research in dynamic response of cement mortars, subjected to repeated cycles of    by measuring ER.Lataste et al. [25] used the ERT technique to characterize damage in field concrete structures, associated with the aperture and depth of surface cracks.On the other hand, for routine quality control measures, more rapid but also reliable methods are needed.Recent studies confirm that electrical resistivity measurement is a simple, reliable, and rapid test method for quality control of concrete [26,27].Concrete has a resistivity varying substantially depending on a number of factors.The electrical current is carried by the dissolved charged ions flowing through the pore solution in the concrete.Therefore, all the factors such as water/cement ratio, cement type, pozzolanic admixtures, and degree of hydration that is affecting the pore structure of the concrete are also affecting the electrical resistivity of the concrete.Environmental factors such as temperature and moisture conditions can also have a large impact on the electrical resistivity of the concrete [26].

Ultrasonic Pulse Velocity (UPV).
The UPV has been successfully used to evaluate the quality of hydraulic concrete for over 60 years.The    of concrete specimens and of concrete in situ can also be estimated through the test UPV [28,29].Ultrasonic waves are mechanical waves that are generated by exciting a piezoelectric crystal with a high voltage pulse.The wave is then transmitted through the tested material, which is in contact with the transducer containing the crystal.When this wave impinges upon a receiving transducer, it produces an output voltage [30].There are three possible ways in which the transducers may be arranged: direct transmission, semidirect transmission, and indirect transmission.In the direct method the transducers are placed on opposite faces of the tested element.In the semidirect method the transducers are placed at a 90 ∘ angle whereas in the indirect method both transducers are placed on one face of the tested element.Since the transmission path length and direction are not well defined in the indirect method, the results produced are less satisfactory.Thus, the indirect method does not give any information on deeper concrete; it only collects information on the surface [31].
In order to use the UPV, a pulse of ultrasonic wave through the concrete is created at a point in the test object surface, and the time it takes to travel from that point to another is measured.Knowing the distance between these two points, the pulse wave velocity can be determined [28].A lot of research has been conducted in order to use this technique to predict the mechanical properties of concrete and rock; for example, Wu and Lin [32] and Selleck et al. [33] investigated the behavior of UPV associated with the damage of concrete.When a solid medium is altered by a dynamic or vibratory load, three kinds of mechanical waves are generated; these waves can be classified by the way they propagate, (also called stress waves): (1) compression waves (also called longitudinal or P waves), (2) shear waves (also called S waves), and (3) surface waves (also known as Rayleigh waves) [34].The velocity of such waves depends on the elastic properties of the medium, so that, knowing the speed of sound and the solid mass, it can estimate the elastic properties of the medium, which can be associated with the quality parameters of the material [35].The relationship between the    and UPV is not unique, and it is affected by many factors such as the size, type, and content of aggregate, w/c ratio, and content of humidity.The effect of these factors has been studied by many researchers, who have clearly indicated that there should be estimated    of the concrete values from the UPV, unless these have similar correlations previously established for the type of concrete studied [36][37][38][39][40].

The Resonance Frequencies Test (RFT).
The method for determining the dynamic modulus of elasticity of solid bodies, from their resonance frequencies, has been used in the last 55 years.However, in recent years, methods of the RFT have been used almost exclusively in laboratory studies [28].
In the late 1930s and early 1940s, other researchers improved the method by using electronic equipment for measuring the resonant frequency [41][42][43]; currently, this type of equipment is required by ASTM C 215-85 Standard Test Method [44].The apparatus consists mainly of two sections: one generates mechanical vibrations and the other detects these vibrations [45].The mechanical vibrations are detected by a piezoelectric transducer and it converts the mechanical vibrations into electric voltage alternating current of the same frequencies [28].
The great advantage of the vibrational methods of testing, apart from their nondestructive nature, is that the dynamic values are obtained from transient loads far below the elastic limit, so that the results are free from time-dependent inelastic strains and hysteresis and are directly related to the internal structure of the material [46].
Vibration-based techniques can be used to monitor the concrete structures.Bagchi et al. applied cost-effective and easy to implement vibration-based damage identification (VBDI) techniques for structural health monitoring of a bridge, based on changes in the dynamic characteristics of a structure to determine the location and extent of damage in the structure.Hsieh et al. described the use of vibrational monitoring in the field of structural analysis for detecting and locating structural damage for the purpose of structural health monitoring [47,48].
On the other hand, the measurement of the dynamic modulus of elasticity with resonance frequency method predictions can be made about the damage degree and deleterious progress, respectively.The dynamic moduli and the damping capacity are the two important properties of engineering materials obtained from vibrational methods of testing.To ensure stability and long life it is necessary to avoid the large amplitudes of vibration associated with resonant frequencies.However, where the nature of the structure and its environment are beyond the control of the designer, the amplitude of vibration at resonance depends on the damping capacity of the structure and its mounting.Knowledge of both the dynamic moduli and the damping characteristics of cementitious materials is thus an important aid to design [46].

Rebound Number Test (RNT).
The RNT method is performed through the device called sclerometer (Schmidt Hammer Rebound).It is based on the principle that the rebound of an elastic mass depends on the hardness of the surface against which the hammer hits; there is a relationship of proportionality between the resistance of concrete and the rebound height of the hammer mass, which is measured by an arbitrary scale of integer numbers.Within certain limits, empirical correlations have been established between the properties of strength of concrete and the RNT [45].The RNT provides information about the strength of the material near the surface where the test is performed; rebound number (RN) is an indicator of the mechanical properties of the material.Although a large rebound number represent a concrete with a higher compressive strength than a concrete with a low rebound, the test is only useful if a relationship can be established between the rebound number and the analyzed resistance of the concrete, assuming that concrete samples have been made with the same type of stone aggregates [10].
For estimation of compressive strength of concrete in structures that are using different testing procedures, compressive strength of concrete can be found as follows: destructive by taking core bore samples from a structure; nondestructive by Schmidt rebound hammer or ultrasonic pulse velocity method; semidestructive test methods or a combination of destructive and nondestructive testing.Standard procedures have been established and are described in detail in European and international standards as European standards: CSN EN 12504-2 [49] and Czech technical standard: CSN 731373 [50] and standards as ASTM C 805-97 [51], JGJ/T 23-2001 [52].
The combination of NDT such as UVP and RNT can be conducted to assess the resistance of the concrete in situ.This combination has an advantage; the Sclerometer reveals information regarding the concrete surface, while ultrasound provides information from the inside of the material [9].In technical literature, there are models proposed by several authors [9,11,[53][54][55][56][57][58][59][60][61][62].Table 1 summarizes the proposed equations to have a fair comparison.
Carbonation of concrete occurs when the carbon dioxide, in the atmosphere in the presence of moisture, reacts with hydrated cement minerals to produce carbonates, for example, calcium carbonate.In this case, the research was carried out with hydraulic concrete with ages from 3 to 28 days.Period during which did not exist the necessary conditions for carbonation should arise in concrete.However we can say that, in older concrete, the carbonation depth can be several millimeters thick and, in extreme cases, up to 20 mm thick.In such cases the rebound numbers can be up to 50% higher than those obtained on an uncarbonated concrete surface.In such cases it is necessary to make an adjustment to the obtained rebounds [10].Moisture condition type of aggregate and carbonation depth of the material can affect the strength estimate [63] and some standards (e.g., Chinese standards JGJ/T 23-2001) offer recommendations for compensating these effects [52].

Significance of Current Research
Nowadays it has been recognized that NDT plays an important role in the condition monitoring of civil infrastructures.To properly maintain the civil infrastructures integrity, engineers required new methods of inspection.Better inspection techniques are needed for damaged infrastructure [64].Being aware of this problem, in this study, four nondestructive    =  0 +  1  2 +  2 UPV 3 +  0.614  1 [53,71]    =  0 +  1  +  2 UPV 0.621  2 [54]    =  0 +  1  3 +  2 UPV 0.616  3 [11]    =  0 +  1  +  2 UPV 4 0.612  4 [55]    =  0 +  1  +  2 UPV +  3 UPV 2 0.596  5 [56,58,61] l n    =  0 +  1  +  2 UPV 0.613  6 [59,60] l n    =  0 +  1 ln  +  2 ln UPV 0.617 techniques were used to predict two important properties to consider in a structure needed for proper maintenance or to know the service conditions of the structure of hydraulic concrete (   and MOR).By combining several nondestructive techniques, robust models that improve the accuracy of the model were generated.In addition, in situ performance conditions are different from controlled laboratory conditions and, so, it is difficult to keep strict quality control.In this paper, the conditions of both field and laboratory for concrete mixtures under study were obtained as reported in Section 4; the forecast mathematical frameworks were generated taking into account these conditions and finally they were fed with data from nondestructive and destructive tests, so these models will have better attachment to the reality.Furthermore a methodology that could be applied to any situation to generate and apply a particular model in the prediction of the properties of the structure under study could be chosen by the particular characteristic of the research material.
And inasmuch as that we can demonstrate the reliability of nondestructive techniques for the analysis and inspection of concrete structures, we will be providing an important tool for engineers that can be exploited to save and optimize resources.

Experimental
The experimental study was divided in two stages.The first one developed three types of specimens: type I: 70 cylinders with dimensions of 100 mm (3.94 in) diameter and 200 mm (7.87 in) height; type II: 70 cylinders with dimensions of 150 mm (5.91 in) diameter and 300 mm (1.81 in) height; type III: 20 beams with dimensions of 150 mm × 150 mm × 600 mm (5.91 in × 5.91 in × 23.62 in).All samples were prepared per ASTM C-143 requirements [65] as reported in Table 2.The w/c ratio varied from 0.55 to 0.76, simulating field concrete specimens.In the second stage, three other types of specimens were developed with the same dimensions to that in the first stage, but with other w/c ratios: type I: 70 cylinders; type II: 70 cylinders; and type III: 40 beams, half of each type specimens were made up with w/c = 0.65 and the rest with w/c = 0.50.These specimens were labeled as concrete made in the laboratory.Table 3 summarizes the measurements of the NDT performed on the specimens, that is, ERT, RFT, UPV, and RNT; the DTs values are also included in the table,    and MOR.It should be noted that all the test specimens were performed under saturated conditions; for structures in situ special equipment is needed for locating the presence of the reinforcing steel and performing the readings of NDT on places away from reinforcing steel to avoid possible errors due to the presence of it [66,67].When NDTs were completed some of the specimens were selected and scheduled for destructive testing.In the first stage of the investigation, 5 cylinders type I, 5 cylinders type II, and 2 beams from each different mixture (see Table 2) were tested.In the second stage only 3 cylinders type I, 3 cylinders type II, and 2 beams were tested due to the assumption that these samples were made in laboratory conditions; therefore, it can be assumed that there is less dispersion of the data from "laboratory" with respect to the information from the "field".The tests were performed at 3, 7, 14, 21, and 28 days.The    and MOR values of hydraulic concrete were calculated in accordance with ASTM [1,2]; the tests were performed

Modeling the Compressive Strength (𝑓 󸀠 𝑐) and Modulus of Rupture (MOR)
As mentioned above, the contribution of this paper is to predict the    and MOR, by proposing a linear equation ( 1), which combines data obtained through four different NDTs: ERT, UPV, RFT, and RNT.The following model was proposed: where   are the coefficients to be calculated and  is either the compressive strength    or the modulus of rupture MOR in MPa.
Even though it was decided to produce the model using a linear equation, it is the possibility that one or some of these NDTs were not correlated with the characteristic being modeled (in this case    or MOR).If this were the case, their inclusion in the model would not be beneficial.For example, let us suppose that all the available variables are used; however, one of these is not correlated with the characteristic being predicted.In this case, it would see that the accuracy of the model on the specimens used to identify the coefficient is the poorest.Unfortunately, given that one test is not correlated with response, then the accuracy of this model on different specimens is not as good as the performance of a model that does not include this noncorrelated test; that is, a model with poor generality can be calculated.
In order to overcome this problem, was decided to create all the linear models that can be built using the available independent variables.Then, our proposed model was identified by following the next procedure.That is, in order to identify the final model, a one-leave-out technique [70] is used in a subset of the total number of specimens.This is performed as follows: let  be the number of specimens, in this subset, used to identify the coefficients of a particular model; then  − 1 specimens are used to identify the coefficients of that model using ordinary least squares, and, then, this model is used to predict the characteristic of the missing specimen.This process is repeated for all the specimens; at the end of this process, the predicted values are compared against the actual values, and, as consequence, the generality of the model is assessed.This procedure is repeated for all the available models and at the end the model selected is the one that presents the best generality.

Results and Discussion
Comparative graphs of the mechanical properties measurements against nondestructive tests carried out were plotted; they were constructed with the point values of the readings recorded for each mixture at different ages (3, 7, 14, 21, and 28 days).Figure 1(a) shows the scatter plot of the ERT measurements versus the    values of cylinders types I and II for the concrete considered as made in situ.This figure is representative of the results obtained in this stage of research.In these samples the ERT increases with the age; in spite of the fact that this is a general trend, where a substantial dispersion of the plotted data is shown.Figure 1(b) shows the representative results of the cylinders types I and II of the second stage (concrete made in the laboratory); unlike the results shown in Figure 1(a), there is a definite upward trend, where the mechanical properties of concrete increase with age and therefore also increase the ERT.The mixtures of the second stage with w/c = 0.65 have a lower ERT than mixtures with w/c = 0.50 due to the content of interconnected saturated pores through which the electric current flows; this is being reflected in the results shown.The behavior described by the ERT is very similar to the one observed for the UPV (Figure 2).The results obtained with the RFT showed no trend.Figure 3 shows representative results obtained on cylinders type I and II at both stages of the investigation, this graph corresponds to cylinders type II with w/c ratio = 0.50 and ages of 3 to 28 days; any trend is not observed in such a graph.
On the other hand, the graphs that correspond to the beams of the first stage have a FRT in longitudinal mode that describes an upward trend with respect to the increase  of the age (see Figure 4(a)), as in the case of ERT versus    (Figure 1(a)) and UPV versus    (Figure 2(a)); where a substantial dispersion of the plotted data is shown that also can be attributed to not controlling the w/c ratio as described above (simulating concrete made in situ).
The FRT and longitudinal and transverse mode observed on beams of the second stage showed an increase of both the MOR and the FRT values as the age of the samples increased, the values for the longitudinal FRT of the beams with a w/c ratio = 0.50 are higher than the frequencies measured on the beams with a w/c ratio = 0.65.Figure 4(b) shows representative results obtained in beams of the second stage; there is an upward trend over the time.
The results obtained from    versus RNT in cylinders type I are shown in Figure 5; it is possible to try to define a trend indicating a rebound increase associated with the age of the concrete; however it can be observed a significant dispersion of data due to the fact that the w/c ratio was not controlled to simulate concrete made in situ; therefore, there are almost the same values of    at ages of 14, 21, and 28 days.
In the common situation where an assessment of material strength is required, it is unfortunate that the complexity of correlation tends to be greatest for the test methods which cause the least damage.Although rebound hammer test causes little damage, is cheap and quick, and is ideal for comparative and uniformity assessments, its correlation to absolute strength prediction poses many problems.Core tests provide the most reliable in situ strength assessment but also cause the most damage and are slow and expensive.However one of the advantages of this paper is the combination of several nondestructive tests to get models that provide better predictions.
Using the data recorded from first stage samples by the NDT, different models found in technical literature [9,11,[53][54][55][56][57][58][59][60][61][62]71] were tested in order to compare other frameworks against the models developed herein.The coefficients of these models are identified using ordinary least squares and the predictions are computed using the Leave-one-out crossvalidation (LOOCV) technique.Their respective correlations factors were obtained using the information of the destructive tests (see Table 1).In order to facilitate the comparison, in this table, our proposed models   and   were also included.It is important to note that these proposed models were generated with the information of the first stage since at this stage we have the RNT.The models proposed in Table 5 (models  0 ,  1 , and  2 ) were generated with the information of the two stages, but eliminating the RNT because there were cylinders with RNT (first stage) and cylinders without RNT (second stage).
The model   corresponds to a proposed model that takes into consideration only the four NDTs; it is a first degree equation with four variables and four constants, with a correlation factor of 0.636.The model   was generated taking the model  0 [9] as a basis, but the ERT and RFT were added to the model; it is a third degree equation with six variables and seven constants, with a correlation factor of 0.642; the variables considered for this model are the NDT and the age and the w/c ratio of the concrete.As can be observed, correlation factors for models   and   are very near to the values of the other frameworks; therefore, these two models can be considered as accurate enough for the predictions of the mechanical properties of concrete.    of the total of 117); due to the fact that the RNT was not performed on the second stage, it was necessary to remove such information in order to generate the models.Figures 6(a) and 6(b) present the correlation coefficient of these models in the training phase.In order to measure the generality of the models selected, a validation was conducted with information of NDT and DTs performed in the 27 remaining cylinders.The best models are summarized in Table 5.
Leave-one-out cross-validation (LOOCV) implies isolation of the data for each iteration, in such a manner that a single sample for the test data and the rest forming the training data are obtained.The evaluation is given by the error; in this type of cross-validation error is very low, but, in contrast, it is computationally very time consuming because the method has to make a large number of iterations, as many as n samples and for each test data (test data and training data) [72].
From the cross-validation analysis, 64 models were obtained as a result of all possible combinations of the variables considered, 32 of type   and 32 of type   ; the two best to predict the    of the hydraulic concrete ( 0 and  1 ) are shown in Table 5.The model  0 is that in which the variables (NDT data) are multiplied by constants; the correlation factor is 0.886 which is greater than the correlation factor obtained by the model  1 , with a correlation factor of 0.858.In case of the beams, models were generated to predict the MOR; in this case all the samples were used to generate and validate the models due to the small number of specimens.32 models were generated obtaining their correlation coefficients (see Figure 6(c)).The best model is  2 with a correlation factor of 0.913 as shown in Table 5.The predictions of the    obtained with the models  0 and  1 are shown in Figures 7(a Figure 8 presents the boxplot of the errors produced by the different models presented on this contribution, that is,  0 ,  1 , and   .This figure is complemented with the boxplot of the error presented by model  0 (see Figure 9).Comparing the boxplot on Figures 8(a), 8(b) and 9-which correspond to models of   -it is observed that these boxplots are very similar indicating that these models are equivalent with respect to their error.Model  1 has a higher variability indicated by the lines extending from the box; however, seeing the size of the box one can consider them similar.

Validation of the Models Using Test of Hypotheses.
According to Walpole et al. [73], "A statistical hypothesis is an assertion or conjecture concerning one or more populations, " in our case, this assertion or conjecture is relative to defining whether or not the information predicted by the different calculated models is equivalent to the real values estimated by destructive testing.In order to determine the validity of the prediction models we assumed a paired observations test of hypotheses; computation of the confidence interval for mean1−mean2 in the situation with paired observations is based on the random variable, presented in (2) [73]: where  and   are random variables representing the sample mean and standard deviation of the differences of the observations between paired observations,  0 is the expected differences between means of the two samples, and  is the number of pairs.Critical regions are constructed using the t-distribution with  − 1 degrees of freedom.
The null hypothesis for this model is that the means of both samples are equal; therefore,  0 becomes zero (alternative hypothesis is the means are different).For a confidence interval of 95% the critical values from t-distribution are −2.056< t > 2.056, and for a confidence interval of 99% critical values are −2.779< t > 2.779.Table 6 shows the computed test statistic (t) for each model.
It can be observed from data in Table 6 that the test statistic lies between critical values for both confidence intervals; therefore the null hypothesis cannot be rejected and it can be assumed that the mean of the predicted values when compared with real data is equivalent; therefore, the models are reliable for these two confidence intervals.

Conclusions
The measurements of the UPV showed the most uniform data when compared against other NDTs; the samples of the second stage indicate velocities of pulse in different ranges; for samples of w/c ratio = 0.50 the velocity is within a range from 3.6 to 3.9 km/s, while for samples of w/c ratio = 0.65 the velocity is within the range of 3.3 to 3.6 km/s.The resonant frequency is the NDT more affected by the morphology or geometry of the analyzed samples.Therefore, the resonance frequency should be used only if prismatic samples that comply with the dimensions indicated in the literature are tested [31].Otherwise it will be difficult to obtain accurate frequency values; furthermore, this test cannot be applied in the field.
The rebound method should be used to determine the homogeneity of the concrete but not its    unless a correlation between the rebound number and    has been made.
The prediction models obtained are  0 and  1 with a correlation factor of 0.895 and 0.889, respectively, to predict the    and the model  2 with a correlation factor of 0.913 to predict the MOR.These models are linear and provide better results than models formed by higher degree polynomials and also provide improved estimations than previous models proposed in the literature.
Nondestructive tests have a fairly wide scope, both in engineering and in other areas of knowledge; therefore it is recommended to continue with investigations using nondestructive techniques for better documentation and reliability thereof.
With the proposed models, it is possible to take readings of nondestructive testing and predicting    or MOR of the structure in the field and samples are not necessary for that purpose.

Figure 1 :Figure 2 :
Figure 1: Compressive strength versus electrical resistivity of cylinder type I with w/c ratio = (a) 0.55, 0.57, 0.61, and 0.76 (concrete made in situ) and (b) 0.50 and 0.65 (concrete made in laboratory) for ages of 3 days to 28 days.

Figure 3 :
Figure 3: Compressive strength versus resonance frequency in longitudinal mode obtained in cylinders type II with w/c ratio = 0.50 and 0.65 for ages of 3 days to 28 days (concrete made in laboratory).

Figure 4 :
Figure 4: Modulus of rupture versus fundamental resonance frequency in longitudinal mode in beams with dimensions of 150 mm × 150 mm × 600 mm, with w/c ratios = (a) 0.55, 0.57, 0.61, and 0.76 (concrete made in situ) and (b) 0.50 and 0.65 (concrete made in laboratory) for ages of 7 day to 28 days.

Figure 5 :
Figure 5: Compressive strength versus hammer rebound obtained in cylinders type I with w/c ratio = 0.57 and 0.76 for ages of 7 day to 28 days (concrete made in situ).

Figure 6 :
Figure 6: Correlation factor of the generated models: (a) of the type   and (b) of the type   with combined cylinders and (c) on beams.

Figure 7 :
Figure 7: Predictions obtained with the models (a)  0 and (b)  1 for cylinders and (c) with the model  2 for beams.

Figure 9 :
Figure 9: Tolerance bounds of the  0 model, model proposed by other authors.

Table 2 :
Slumps of mixtures made.

Table 3 :
Nondestructive and destructive testing applied to the specimens.

Table 4 :
Materials used in concrete mixtures.The mixtures were designed for    = 24.52MPa and    = 16.67 MPa.

Table 5
summarizes our models for prediction of the    and MOR; the correlation factor values are indicated.The model was constructed with NDT data from both stages of the experimental study (90 cylinders randomly selected 55, 0.57, 0.61, and 0.76 (concrete made in situ) and (b) 0.50 and 0.65 (concrete made in laboratory) for ages of 7 day to 28 days.

Table 6 :
Computed test statistic () for each of the prediction models.