Comparative Analyses of Selected Neural Networks for Prediction of Sustainable Cementitious Composite Subsurface Tensile Strength

: The article assesses comparative analyses of some selected machine-learning algorithms for the estimation of the subsurface tensile strength of cementitious composites containing waste granite powder. Any addition of material to cementitious composites causes their properties to differ; therefore, there is always a need to prepare a precise model for estimating these properties’ values. In this research, such a model of prediction of the subsurface tensile strength has been carried out by using a hybrid approach of using a nondestructive method and neural networks. Moreover, various topologies of neural networks have been evaluated with different learning algorithms and number of hidden layers. It has been proven by the very satisfactory results of the performance parameters that such an approach might be used in practice. The errors values (MAPE, NRMSE, and MAE) of this model range from 10 to 12%, which, in the case of civil engineering practice, proves that this model is sufﬁcient for being used. This novel approach can be a reasonable alternative for evaluating the properties of spacious cementitious composite elements where there is a need to analyse not only the compressive strength but also its subsurface tensile strength.


Introduction
A composite is defined as a material made up of at least two components with different properties such that it has superior and/or new properties over those components taken separately.However, the given definition is not generally accepted and precise, but is intended to generalize the most common descriptive definitions.On its basis, it can be assumed that a cement composite is a composite of whom one of the basic components is cement.The most popular composites in which cement acts as a binder are: grout, mortar, and concrete [1].By the term grout, we mean the material that results from the combination of water and cement.The leaven itself is rarely used in construction.Cement mortar is used much more often, e.g., for bricklaying, plastering, or making floor underlays.Standard cement mortar is a mixture of the appropriate proportions of cement, water, and sand.The composition of the mortar can be supplemented with additives and admixtures added to give the appropriate properties to the mixture or hardened mortar.Concrete is another material.It is the most widespread cement composite in the world.In addition to water, sand, and cement, concrete contains coarse aggregate.In concrete, as in the case of mortar, we use admixtures and additives, which, in standard concretes, allow us, for example, to obtain better working parameters of the mixture with the use of less water, or to improve the physical parameters of the hardened composite [2].
An additive to concrete can be called a fine-grained component added to the concrete mix, whose task is to improve its properties or obtain the special parameters of hardened concrete.The amount of additive is generally from 5% to 15% of the weight of the cement.The content of the additive should be taken into account in subsequent iterations of the Appl.Sci.2023, 13, 4817 2 of 22 calculations of the designed mixtures.There are two main types of additives.The first type is nearly inert additives, which are in turn divided into fillers and pigments.Pigments are used to tint the concrete to the desired color.The role of fillers is to fill the pores of the concrete mixture [3].An additive typical for this group is limestone powder.Additives with latent hydraulic or pozzolanic properties are the second type of concrete additives.Their use leads to an improvement in the properties of concrete or the acquisition of special properties.Representatives of this group are silica fume or fly ash.The most popular methods of adding fillers to a cement composite are in the form of (i) cement replacement, (ii) aggregate replacement, and (iii) grout replacement.These methods are successfully applied to waste materials in the form of quartz dust, marble powder, powdered glass, or fly ash, added to concrete and mortar as fillers or additives that modify the properties of the mixture or the hardened cement composite [4].
The mere production of cement as a basic component of concretes and mortars also has an undesirable impact on the environment.During the production of 1 ton of this material, about 900 kg of carbon dioxide are emitted; therefore, the partial replacement of cement with another component may contribute to reducing CO 2 emissions to the atmosphere [5].
Granite powder is a waste material that results from the crushing or cutting of granite rocks for industrial purposes.In recent years, an increase in the amount of granite powder waste has been observed as a result of the development of the granite industry.The powder fraction of this material alone is currently not applicable, which exacerbates the problem of the proper management of this product.The improper storage of dry material in the form of heaps, which can easily spread in the environment, causes water pollution, changes in the pH of land intended for agricultural purposes, and air pollution in the vicinity of landfills.Insisting that the effects of the environmental degradation by granite powder can be felt directly by humans, it is necessary to take the appropriate measures to reduce the amount of this production residue in the environment.One way to reduce its harmful effects is to use this by-product for the sustainable production of cement composites as the most common building material and use it as a supplementary cementitious material (SCM) [6].
However, when another material is added to cementitious composites, it is necessary to evaluate how it affects their properties.Most of the tests used for evaluating the materials' properties in civil engineering are destructive; therefore, they are costly and time-consuming, and they produce wastes.Therefore, sometimes for this purpose, the nondestructive tests are reasonable substitutes [7].It is also important while performing such tests in buildings which are operating or which are characterized as historical heritage sites and cannot be tested destructively [8].It was previously shown that an NDT method such as the Schmidt hammer method was successful in evaluating the compressive strength of concrete [9,10].However, there is an issue with evaluating the other properties of cementitious composites such as subsurface tensile strength.For this purpose, very often, a method such as the pull-off test is used, but, after the test, the sample needs to be repaired, which makes this method not fully nondestructive [11].However, it is still possible to try to correlate the compressive strength with the subsurface tensile strength.Therefore, it is possible to try to evaluate the opposite, which means that it is possible to evaluate the subsurface tensile strength using a Schmidt hammer [12].Evaluating the subsurface tensile strength in cementitious composites such as mortars dedicated for floor substrate is very important because of the fact that, very often, these parameters are responsible for the durability of the element [13].
Recently, very often, soft computing techniques such as artificial neural networks are used for this purpose, while they are able to strengthen the results of NDT methods by improving the performance of these methods, which might be necessary in this case [14].Such approaches have previously been successful when using the combination of a Schmidt hammer and neural networks to predict the unconfined compressive strength of rocks [15], to predict the damage analyses of concrete subjected to high temperature [16], or to assess the load-carrying capacity of concrete anchor bolts [17].
Taking the aforementioned thoughts into account, it might be reasonable to use such a combination for evaluating the subsurface tensile strength of cementitious composites using the combination of a Schmidt hammer and neural networks.Therefore, for this purpose, the authors propose the hybrid method of using the nondestructive Schmidt hammer method and neural networks for evaluating the properties of cementitious composites containing granite powder that can be used in the civil-engineering practice.This solution might be beneficial for evaluating the properties of existing floors made of mortars which are located in operating buildings.

Materials
In order to create the database used for ANN modeling, the research part was necessary.Four different cementitious composite mixes were prepared with varying compositions in terms of the amount of granite powder used, as shown in Table 1.The first mix, which is the reference mix designated REF, contained only Portland cement as a binder, and is typical mortar mixture used in floors.In the other mixes, Portland cement was partially replaced by granite powder, obtained from a local stone plant in Strzegom, in proportions of 10%, 20%, and 30% (designated GP10, GP20, and GP30, respectively).The amount of water was assumed to be constant for all mixes and was determined for a reference composition for a water-cement ratio of 0.5.CEM I 42.5 cement, which consisted of 95-100% Portland clinker and 0-5% secondary components, was used for the mortars.Its normal compressive strength was ≥42.5 MPa and ≤62.5 MPa.The sand used in the study had a specific density of about 1600 kg/m 3 and a water absorption of W (A,24) = 0.6%.The samples were made into 500 mm × 500 mm × 40 mm slabs, which were then cured for the required time (56 days or 90 days).All samples were compacted and manually trowelled.The samples were cured in air (AIR) and in a moist environment (WET), created by applying moist sponges to the surface of the samples, at a temperature of 18 ± 3 • C. The condition of the environment in which the samples were cured was maintained for 28 days, which, in the case of the wet environment, consisted of checking the water saturation of the sponges and re-wetting them if necessary.It was recommended that the sponges not be saturated, so, after wetting, they were left for 15-30 min to drain excess water.In order to better evaluate the surface layer of the cement mortars made, a Schmidt hammer test was performed to correlate the results with the compressive strength of the surface layer of the specimens.The test is carried out by placing the hammer at right angles to the test surface and gradually increasing the pressure until the hammer strikes.At the moment of impact, the number of rebounds is recorded, representing a test of the hardness of the material.Adjacent test points shall not be closer than 25 mm to either each other or to the edge.If the test results in crushing or damage to the mortar at the measurement point, the result shall be rejected.The hammer was checked before the test on a calibration steel anvil-6 control measurements were taken, each indicating a reflection number of 80 ± 2. This allowed the artificial neural networks to increase the database on which they worked.The test was conducted according to the method of [18].Twelve impact tests were performed, 3 at each pull-off strength test location.
The approximate strength f c was taken as equal to the value of f L .For the obtained reflection number (L), the compressive strength value was calculated from Equations ( 1) and ( 2) based on the base curve (Figure 1a): Appl.Sci.2023, 13, x FOR PEER REVIEW 4 of 27 measurement point, the result shall be rejected.The hammer was checked before the test on a calibration steel anvil-6 control measurements were taken, each indicating a reflection number of 80 ± 2. This allowed the artificial neural networks to increase the database on which they worked.The test was conducted according to the method of [18].Twelve impact tests were performed, 3 at each pull-off strength test location.The approximate strength fc was taken as equal to the value of fL.For the obtained reflection number (L), the compressive strength value was calculated from Equations ( 1) and ( 2) based on the base curve (Figure 1a):

Pull-Off Tests
Four boreholes were drilled on the prepared panels after each curing time as shown in Figure 2a.The holes were drilled with a 50 mm inside diameter, using diamond saw to a depth of 8 mm.This will clarify the tests directly on the near-surface layer, which was defined as a range from 0 to 4 mm depth (Figure 2b).According to the standard [19], to achieve meaningful results, it was necessary to place the edge of the disc a minimum of 50 mm from the edge of the specimen.After drilling, the surface of the discs was vacuumed and cleaned with acetone without compromising the integrity of the coating.Then, in order to connect the Proceq Dy-216 test device to the specimen, the steel discs (cleaned with acetone) were glued to the wells with Poxipol two-component adhesive.The glue was applied to the surface of the puck and the surface of the sample, then positioning the puck by carefully pressing it to remove air bubbles.The adhesive had no effect on the test specimen.Once the adhesive was fully cured, the test was proceeded.The connection between the puck and the test instrument was screwed into the puck.The machine was levelled and programmed to complete the test in a maximum of 100 s at a loading rate of 0.05 MPa/s.The result of the test is recorded on the instrument along with a graph of load increase over time.After entering the relevant data into the apparatus, the pull-off strength was automatically calculated from Equation (3).Regarding the work of [20] and [21] where the pull-off strength between the substrate layer and the overlay layer is predicted using ANN, the network performed on the basis of the obtained results is supposed to allow for us to estimate the pull-off strength of the near-surface layer as the most exposed layer.

Pull-Off Tests
Four boreholes were drilled on the prepared panels after each curing time as shown in Figure 2a.The holes were drilled with a 50 mm inside diameter, using diamond saw to a depth of 8 mm.This will clarify the tests directly on the near-surface layer, which was defined as a range from 0 to 4 mm depth (Figure 2b).According to the standard [19], to achieve meaningful results, it was necessary to place the edge of the disc a minimum of 50 mm from the edge of the specimen.After drilling, the surface of the discs was vacuumed and cleaned with acetone without compromising the integrity of the coating.Then, in order to connect the Proceq Dy-216 test device to the specimen, the steel discs (cleaned with acetone) were glued to the wells with Poxipol two-component adhesive.The glue was applied to the surface of the puck and the surface of the sample, then positioning the puck by carefully pressing it to remove air bubbles.The adhesive had no effect on the test specimen.Once the adhesive was fully cured, the test was proceeded.The connection between the puck and the test instrument was screwed into the puck.The machine was levelled and programmed to complete the test in a maximum of 100 s at a loading rate of 0.05 MPa/s.The result of the test is recorded on the instrument along with a graph of load increase over time.After entering the relevant data into the apparatus, the pull-off strength was automatically calculated from Equation (3).Regarding the work of [20] and [21] where the pull-off strength between the substrate layer and the overlay layer is predicted using ANN, the network performed on the basis of the obtained results is supposed to allow for us to estimate the pull-off strength of the near-surface layer as the most exposed layer.

Neural Networks
As we work to reduce the use of cement in the construction of subfloors, replacing it with waste granite powder, it is necessary to determine the optimal values of the components while maintaining the appropriate parameters of the final product.In order to optimize the process of designing the composition of cement mortars with cement replacements, artificial neural networks (ANN) in the form of multilayer perceptron (MLP) were used.MLP is one of the most recognized ANNs used when solving design problems in the construction industry.Karthiyaini et al. [22] predicted the mechanical strength of concrete, doped with fibers, using a multiple regression model (MRA) and an artificial neural network (ANN).The purpose of their work was to compare the two models with each other.The work of Mehdi et al. [23] describes the use of evolutionary neural networks (EANN), as a combination of artificial neural network (ANN) and evolutionary search procedures, in predicting the compressive strength of concrete.Panagiotis G. Asteris et al. [24] trained a database, in conventional machine learning models, of an artificial neural network (ANN) to predict the compressive strength of concrete.Czarnecki et al. [25] presented the application of artificial neural networks (ANNs) to nondestructively evaluate the adhesion of a repair overlay to a concrete substrate.The subject of the paper by Czarnecki et al. [26] was the prediction of the strength of a cementitious composite with the addition of ground granulated blast-furnace slag.There, they compared an artificial neural network (ANN) with a self-organizing feature map (SOFM).This tool is eagerly used by researchers to work with cementitious composites because it can handle inconsistent or "noisy" information which is a common occurrence when studying these materials.According to Figure 3, an artificial neural network in the form of an MLP consists of three layers: input, output, and one or two hidden.Each layer can have one or more neurons.Before interpreting the results, such networks should be subjected to training.

Neural Networks
As we work to reduce the use of cement in the construction of subfloors, replacing it with waste granite powder, it is necessary to determine the optimal values of the components while maintaining the appropriate parameters of the final product.In order to optimize the process of designing the composition of cement mortars with cement replacements, artificial neural networks (ANN) in the form of multilayer perceptron (MLP) were used.MLP is one of the most recognized ANNs used when solving design problems in the construction industry.Karthiyaini et al. [22] predicted the mechanical strength of concrete, doped with fibers, using a multiple regression model (MRA) and an artificial neural network (ANN).The purpose of their work was to compare the two models with each other.The work of Mehdi et al. [23] describes the use of evolutionary neural networks (EANN), as a combination of artificial neural network (ANN) and evolutionary search procedures, in predicting the compressive strength of concrete.Panagiotis G. Asteris et al. [24] trained a database, in conventional machine learning models, of an artificial neural network (ANN) to predict the compressive strength of concrete.Czarnecki et al. [25] presented the application of artificial neural networks (ANNs) to nondestructively evaluate the adhesion of a repair overlay to a concrete substrate.The subject of the paper by Czarnecki et al. [26] was the prediction of the strength of a cementitious composite with the addition of ground granulated blast-furnace slag.There, they compared an artificial neural network (ANN) with a self-organizing feature map (SOFM).This tool is eagerly used by researchers to work with cementitious composites because it can handle inconsistent or "noisy" information which is a common occurrence when studying these materials.According to Figure 3, an artificial neural network in the form of an MLP consists of three layers: input, output, and one or two hidden.Each layer can have one or more neurons.Before interpreting the results, such networks should be subjected to training.

Neural Networks
As we work to reduce the use of cement in the construction of subfloors, replacing it with waste granite powder, it is necessary to determine the optimal values of the components while maintaining the appropriate parameters of the final product.In order to optimize the process of designing the composition of cement mortars with cement replacements, artificial neural networks (ANN) in the form of multilayer perceptron (MLP) were used.MLP is one of the most recognized ANNs used when solving design problems in the construction industry.Karthiyaini et al. [22] predicted the mechanical strength of concrete, doped with fibers, using a multiple regression model (MRA) and an artificial neural network (ANN).The purpose of their work was to compare the two models with each other.The work of Mehdi et al. [23] describes the use of evolutionary neural networks (EANN), as a combination of artificial neural network (ANN) and evolutionary search procedures, in predicting the compressive strength of concrete.Panagiotis G. Asteris et al. [24] trained a database, in conventional machine learning models, of an artificial neural network (ANN) to predict the compressive strength of concrete.Czarnecki et al. [25] presented the application of artificial neural networks (ANNs) to nondestructively evaluate the adhesion of a repair overlay to a concrete substrate.The subject of the paper by Czarnecki et al. [26] was the prediction of the strength of a cementitious composite with the addition of ground granulated blast-furnace slag.There, they compared an artificial neural network (ANN) with a self-organizing feature map (SOFM).This tool is eagerly used by researchers to work with cementitious composites because it can handle inconsistent or "noisy" information which is a common occurrence when studying these materials.According to Figure 3, an artificial neural network in the form of an MLP consists of three layers: input, output, and one or two hidden.Each layer can have one or more neurons.Before interpreting the results, such networks should be subjected to training.For the purpose of this study, the MATLAB program was used for designing the neural networks' structures and learning algorithms as: gradient descent algorithm, conjugate gra-dient algorithm, and the Broyden-Fletcher-Goldfarb-Shanno and Levenberg-Marquardt algorithms.It also allows the use of 5 types of hidden layer activation functions and output layer activation functions.During ANN modelling, calculations were made for all possible combinations of the above-mentioned network elements.In addition, the number of hidden layer neurons was changed in the range of 1 to 15.We changed the number of layers of hidden neurons to 2 and again edited the number of neurons in layer 1 and layer 2 in the range of 1 to 15.Each of these functions processes the data using the corresponding Formulaes ( 4)-( 8): Based on the database created, each input variable has 7 parameters, including: (1) the amount of cement (C), ( 2) the amount of sand (S), (3) the amount of water (W), ( 4) the amount of granite powder (GP), ( 5) the curing method (dry or wet), ( 6) the curing time, and ( 7) the compressive strength tested with a Schmidt hammer (fc).The output parameter consists of the tensile strength values of the floor mortar substrate measured using pulloff method.

Schmidt Hammer Test
The samples were compared with each other based on the amount of granite powder added to partially replace the cement, as well as based on the maturation conditions of the samples.Figure 4 shows the average compressive strengths obtained when testing each of the sleepers.The basic requirement for cementitious composite substrates in terms of compressive strength is to achieve a minimum value of 20 MPa.It can be seen that all of the samples made reached the minimum value; therefore, the results obtained were used to make a database.From the graph, it can be seen that, as the amount of granite powder in the mix increases, the compressive strength of the substrate decreases, but even for the samples with the highest filler content, the substrate achieved the required strength.

Pull-Off Tests
The samples were compared to each other due to the amount of granite powder added to partially replace the cement, as well as the conditions and maturation time of the samples.Figure 5 shows the average pull-off strengths of each mixture.First of all, it can be seen that there is a large difference in the results due to the different maturation conditions of the samples.Variation in testing time also affected the results.Likewise, the

Pull-Off Tests
The samples were compared to each other due to the amount of granite powder added to partially replace the cement, as well as the conditions and maturation time of the samples.Figure 5 shows the average pull-off strengths of each mixture.First of all, it can be seen that there is a large difference in the results due to the different maturation conditions of the samples.Variation in testing time also affected the results.Likewise, the amount of granite powder replacing the cement made a difference.It can be noted that, with the amount of 10% filler, the obtained results for the air-matured samples increased, which suggests the validity of using this modern waste approach to mortar design.

Pull-Off Tests
The samples were compared to each other due to the amount of granite powder added to partially replace the cement, as well as the conditions and maturation time of the samples.Figure 5 shows the average pull-off strengths of each mixture.First of all, it can be seen that there is a large difference in the results due to the different maturation conditions of the samples.Variation in testing time also affected the results.Likewise, the amount of granite powder replacing the cement made a difference.It can be noted that, with the amount of 10% filler, the obtained results for the air-matured samples increased, which suggests the validity of using this modern waste approach to mortar design.

Statistical Analyses of the Collected Data
The whole dataset is presented in Appendix A. The minimum, maximum, mean, standard deviation, and range of the parameters included in the database are shown in Table 2.The correlation coefficients between all possible variables are shown in Table 3. High positive or negative correlation coefficient values between input variables can result in low efficiency and difficulty in determining the impact of these variables on the response.There are no significant correlations between the independent input variables, which was to be expected.As expected, it is clear that there is a strong correlation between the mortar f h and input parameters, such as the amount of cement and granite powder, water-to-cement ratio (W/C), and curing method.In addition, Figure 6 shows the frequency histogram of the six input parameters.These figures are extremely useful because they can help identify the range of parameter values where the data are insufficient and where more data are required.The whole dataset is presented in Appendix A. The minimum, maximum, mean, standard deviation, and range of the parameters included in the database are shown in Table 2.The correlation coefficients between all possible variables are shown in Table 3. High positive or negative correlation coefficient values between input variables can result in low efficiency and difficulty in determining the impact of these variables on the response.There are no significant correlations between the independent input variables, which was to be expected.As expected, it is clear that there is a strong correlation between the mortar fh and input parameters, such as the amount of cement and granite powder, water-to-cement ratio (W/C), and curing method.In addition, Figure 6 shows the frequency histogram of the six input parameters.These figures are extremely useful because they can help identify the range of parameter values where the data are insufficient and where more data are required.

Neural Network Analyses
Using numerical methods of design and prediction, it is necessary to create a database.Based on the results of the study, a database was created containing information on the material composition of the four previously mentioned mortars, how and when they were matured, and the results of the tests presented in Sections 2 and 3. To meet the requirements of the PN-EN 12504-2 standard of the Schmidt hammer tests, the authors decided to performed the four pull-off tests for each composition and maturing conditions, and, in each of these areas, the three Schmidt hammer tests were performed.Because these three tests inside one pull-off area is less than that stated in the standard, the authors have not calculated the average value out of them.With this, a database of 192 data sets was collected, which were divided into a set for the learning process (136 data sets-70%), the testing process (28 data sets-15%), and the validation process (28 data sets-15%), in a machine-learning model with artificial neural networks (ANN).
Using the ranking method, the most effective artificial neural network was determined, whose calculated values of R 2 (9), NRMSE (10), MAE (11), and MAPE (12) were, at the same time, closest to the ideal values of the given indicators.
In the given formulae, x i is the i-th actual value obtained in the experimental part, while xi is the i-th value predicted by the network.The total number of data is described by r, and x mean is the average of the real values, i.e., the experimental value f h of the cement composite.
Figure 7 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the gradient descent learning algorithm with one layer of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with two hidden-layer neurons.It obtained the following network fit evaluation indices: R 2 = 0.7807, NRMSE = 0.1407, MAE = 0.1460, and MAPE = 11.4751%. Figure 6 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match, while the dashed lines represent the 20% error area.About 16% of the results are outside this area.
Figure 8 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the Broyden-Fletcher-Goldfarb-Shanno learning algorithm with one layer of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with nine hidden-layer neurons.It obtained the following network fit evaluation indices: R 2 = 0.8399, NRMSE = 0.1028, MAE = 0.1290, and MAPE = 9.8658%.Figure 7 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match while the dashed lines represent the 20% error area.About 6% of the results are outside this area.Figure 8 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the Broyden-Fletcher-Goldfarb-Shanno learning algorithm with one layer of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with nine hidden-layer neurons.It obtained the following network fit evaluation indices: = 0.8399 , = 0.1028 , = 0.1290, and = 9.8658%.Figure 7 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function = represents a perfect match while the dashed lines represent the 20% error area.About 6% of the results are outside this area.Figure 9 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the Levenberg-Marquardt learning algorithm with one layer of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with seven hidden-layer neurons.It obtained the following network fit evaluation indices: = 0.7934 , = 0.1261 , = 0.1410 , and = 10.4077%.Figure 8 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function = represents a perfect match while the dashed lines represent the 20% error area.About 9% of the results are outside this area.Figure 9 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the Levenberg-Marquardt learning algorithm with one layer of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with seven hidden-layer neurons.It obtained the following network fit evaluation indices: R 2 = 0.7934, NRMSE = 0.1261, MAE = 0.1410, and MAPE = 10.4077%.Figure 8 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match while the dashed lines represent the 20% error area.About 9% of the results are outside this area.Figure 10 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the conjugate gradient learning algorithm with one layer of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with three hidden-layer neurons.It obtained the following network fit evaluation indices: = 0.8134 , = 0.1109 , = 0.1391 , and = 10.3214%. Figure 9 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function = represents a perfect match while the dashed lines represent the 20% error area.About 10% of the results are outside this area.Figure 10 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the conjugate gradient learning algorithm with one layer of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with three hidden-layer neurons.It obtained the following network fit evaluation indices: R 2 = 0.8134, NRMSE = 0.1109, MAE = 0.1391, and MAPE = 10.3214%. Figure 9 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match while the dashed lines represent the 20% error area.About 10% of the results are outside this area.Figure 11 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the gradient descent learning algorithm with two layers of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with 13 neurons in the first hidden layer and 15 neu-  Figure 10 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match while the dashed lines represent the 20% error area.About 18% of the results are outside this area.Figure 12 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the Broyden-Fletcher-Goldfarb-Shanno learning algorithm with two layers of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with four neurons in the first hidden layer and three neurons in the second hidden layer.It obtained the following network fit evaluation indices: = 0.7657 , = 0.1382 , = 0.1507 , and = 11.3525%.Figure 11 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function = represents a perfect match while the dashed lines represent the 20% error area.About 11.5% of the results are outside this area.Figure 12 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the Broyden-Fletcher-Goldfarb-Shanno learning algorithm with two layers of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with four neurons in the first hidden layer and three neurons in the second hidden layer.It obtained the following network fit evaluation indices: R 2 = 0.7657, NRMSE = 0.1382, MAE = 0.1507, and MAPE = 11.3525%.Figure 11 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match while the dashed lines represent the 20% error area.About 11.5% of the results are outside this area.
Figure 13 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the Levenberg-Marquardt learning algorithm with two layers of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with four neurons in the first hidden layer and 10 neurons in the second hidden layer.It obtained the following network fit evaluation indices: R 2 = 0.8046, NRMSE = 0.1112, MAE = 0.1481, and MAPE = 11.0211%.Figure 12 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match while the dashed lines represent the 20% error area.About 8% of the results are outside this area.Figure 13 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the Levenberg-Marquardt learning algorithm with two layers of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with four neurons in the first hidden layer and 10 neurons in the second hidden layer.It obtained the following network fit evaluation indices: = 0.8046 , = 0.1112 , = 0.1481 , and = 11.0211%.Figure 12 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function = represents a perfect match while the dashed lines represent the 20% error area.About 8% of the results are outside this area.Figure 14 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the conjugate gradient learning algorithm with two layers of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with 11 neurons in the first hidden layer and eight neurons in the second hidden layer.It obtained the following network fit evaluation indices: = 0.7734, = 0.1332, = 0.1488, and = 10.8728%. Figure 13 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function = represents a perfect match while the dashed lines represent the 20% error area.About 11% of the results are outside this area.Figure 14 illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the conjugate gradient learning algorithm with two layers of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with 11 neurons in the first hidden layer and eight neurons in the second hidden layer.It obtained the following network fit evaluation indices: R 2 = 0.7734, NRMSE = 0.1332, MAE = 0.1488, and MAPE = 10.8728%. Figure 13 shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match while the dashed lines represent the 20% error area.About 11% of the results are outside this area.

Comparison of the Models
The network allows us to estimate the strength based on the material composition of the floor mortar and the compressive strength value determined by the Schmidt hammer.This allows the use of a nondestructive method to evaluate the peel strength, thereby not exposing the tested surface to the repairs necessary after using the pull-off method.Subsequently, using the Schmidt hammer, it is possible to obtain both the compressive strength and pull-off strength.These are the two most important characteristics of floor substrates.
The machine-learning algorithms used in the article are accurate in predicting the subsurface tensile strength of green cementitious composites.This can be concluded from the values of the parameters determining their accuracy, which were R 2 , NRMSE, MAE, and MAPE.Figures 15-19 and Table 4 compare the above algorithms.It can be seen that all the models tested predict the tensile strength of green cement composites well.The BFGS algorithm with one layer of hidden neurons performed best.However, since none of the models were perfectly accurate (R 2 would be equal to 1; and NRMSE, MAE, MAPE, and relative error would be equal to 0), there is still room for improvement in the algorithms by building other databases or using other algorithms.
The network based on the BFGS algorithm with nine hidden-layer neurons performed best in predicting the pull-off strength of cement floor underlayments.It obtained the following network fit evaluation indices: = 0.8399 , = 0.1637 , = 0.1290, and = 9.8658 %.According to the article [10], the values of determination coefficients achieved by the selected network are within the range of a very good match with the results obtained experimentally because > 0.81.The MAPE absolute percentage error values of up to 10% are defined as a good fit.As can be seen in Figure 6, the presented network correctly predicts the pull-off strength of the surface layer.The strength ranging from 0.8 MPa to 2.6 MPa is predicted with similar accuracy.The result obtained for the experimental strength above 3 MPa significantly deviates from the rest of the results.This may be due to insufficient batch data from the 2.8 MPa to 3.2 MPa range.

Comparison of the Models
The network allows us to estimate the strength based on the material composition of the floor mortar and the compressive strength value determined by the Schmidt hammer.This allows the use of a nondestructive method to evaluate the peel strength, thereby not exposing the tested surface to the repairs necessary after using the pull-off method.Subsequently, using the Schmidt hammer, it is possible to obtain both the compressive strength and pull-off strength.These are the two most important characteristics of floor substrates.
The machine-learning algorithms used in the article are accurate in predicting the subsurface tensile strength of green cementitious composites.This can be concluded from the values of the parameters determining their accuracy, which were R 2 , NRMSE, MAE, and MAPE.Figures 15-19 and Table 4 compare the above algorithms.It can be seen that all the models tested predict the tensile strength of green cement composites well.The BFGS algorithm with one layer of hidden neurons performed best.However, since none of the models were perfectly accurate (R 2 would be equal to 1; and NRMSE, MAE, MAPE, and relative error would be equal to 0), there is still room for improvement in the algorithms by building other databases or using other algorithms.It should be noted that the prediction accuracy in the learning, testing, and validation processes of the given network was similar.It should be noted that the prediction accuracy in the learning, testing, and validation processes of the given network was similar.The network based on the BFGS algorithm with nine hidden-layer neurons performed best in predicting the pull-off strength of cement floor underlayments.It obtained the following network fit evaluation indices: R 2 = 0.8399, RMSE = 0.1637, MAE = 0.1290, and MAPE = 9.8658 %.According to the article [10], the values of determination coefficients R 2 achieved by the selected network are within the range of a very good match with the results obtained experimentally because R 2 > 0.81.The MAPE absolute percentage error values of up to 10% are defined as a good fit.As can be seen in Figure 6, the presented network correctly predicts the pull-off strength of the surface layer.The strength ranging from 0.8 MPa to 2.6 MPa is predicted with similar accuracy.The result obtained for the experimental strength above 3 MPa significantly deviates from the rest of the results.This may be due to insufficient batch data from the 2.8 MPa to 3.2 MPa range.It should be noted that the prediction accuracy in the learning, testing, and validation processes of the given network was similar.

Conclusions
In this work, the authors proposed a hybrid neural-nondestructive model of subsurface tensile strength prediction in cementitious composites containing granite powder.For this purpose, the experimental tests using a Schmidt hammer were performed, combined with numerical analyses using neural networks.It has been proven that it is possible to design such accurate models using a neural network with the Broyden-Fletcher-Goldfarb-Shanno learning algorithm.The performance was very good, which was shown by the very high value of the coefficient of determination R 2 , equal to 0.8399, and very low values of the errors presented as NRMSE, MAE, and MAPE, equal to 0.1028, 0.1290, and 9.8658%, respectively.According to the results of the correlation between the input variables and output parameter, it can be seen that the amount of cement content and curing conditions contributed the most to the output parameter, whether the samples were kept in a dry environment or wet.
One of the advantages of this method is the fact that it can be used in existing structures, without a need to obtain the samples and destroy them.Moreover, it can still be used after almost 100 days of concreting.Because different curing conditions were taken into account, this method is more universal in comparison to standardized equations.
This novel approach can be a reasonable alternative for evaluating the properties of spacious mortar elements in which there is a need to analyse not only the compressive strength but also its subsurface tensile strength.Moreover, it is worth emphasizing that using a Schmidt hammer during the tests allows us to evaluate the compressive strength, which is beneficial in order to shorten the time of analysing the cementitious composite element.
However, there are some limitations which should be considered before using it in practice.Firstly, it has not been tested yet in the field; therefore, there is need for experimental verification.Second, it has not been used to evaluate the subsurface tensile strength of composites containing admixtures other than granite powder, which might be interesting from a cognitive point of view.Moreover, even though the mortars have been tested, there is still an open question if this model is suitable for concrete containing core aggregates.
Taking the above into account, the authors are still investigating the possibility of using these methods for different cementitious composite mixtures containing such admixtures as fly ash and blast-furnace slag.Moreover, the authors are trying different algorithms such as SVM, linear regression, or random forest, which might be more accurate than those presented in the article.Because this is still an open problem, the authors submit the data to the article for other researchers to use.

Figure 3 .
Figure 3. Block diagram of an artificial neural network.

Figure 4 .
Figure 4. Average compressive strength of the near-surface layer fc.

Figure 4 .
Figure 4. Average compressive strength of the near-surface layer f c .

Figure 4 .
Figure 4. Average compressive strength of the near-surface layer fc.

Figure 5 .
Figure 5. Average pull-off strength fh of the near-surface layer fc.

Figure 5 .
Figure 5. Average pull-off strength f h of the near-surface layer f c .

Figure 6 .
Figure 6.The relation between input variables and tensile strength.

Figure 7 .
Figure 7. Relationship between the predicted value and the experimental value of the pull-off strength fh for ANN test, train, and validation process for the GD algorithm with one layer of hidden neurons.

Figure 7 .Figure 8 .
Figure 7. Relationship between the predicted value and the experimental value of the pull-off strength f h for ANN test, train, and validation process for the GD algorithm with one layer of hidden neurons.Appl.Sci.2023, 13, x FOR PEER REVIEW 12 of 27

Figure 8 .
Figure 8. Relationship between the predicted value and the experimental value of the pull-off strength f h for ANN test, train, and validation process for the BFGS algorithm with one layer of hidden neurons.

Figure 9 .
Figure 9. Relationship between the predicted value and the experimental value of the pull-off strength fh for ANN test, train, and validation process for the LM algorithm with one layer of hidden neurons.

Figure 9 .
Figure 9. Relationship between the predicted value and the experimental value of the pull-off strength f h for ANN test, train, and validation process for the LM algorithm with one layer of hidden neurons.

Figure 10 .
Figure 10.Relationship between the predicted value and the experimental value of the pull-off strength fh for ANN test, train, and validation process for the CG algorithm with one layer of hidden neurons.

Figure 10 .
Figure 10.Relationship between the predicted value and the experimental value of the pull-off strength f h for ANN test, train, and validation process for the CG algorithm with one layer of hidden neurons.

Figure 11
Figure 11  illustrates the results of the predicted values of the pull-off strength in relation to the results obtained by the experimental method using the gradient descent learning algorithm with two layers of hidden neurons.With this algorithm, the best ability to predict the subsurface tensile strength of cementitious composite containing granite powder was obtained for a network with 13 neurons in the first hidden layer and 15 neurons in the second hidden layer.It obtained the following network fit evaluation indices: R 2 = 0.7386, NRMSE = 0.1262, MAE = 0.1582, and MAPE = 11.8732%.Figure10shows the results for the learning process (blue rhombuses), the testing process (pink triangles), and the validation process (green circles).The red solid line with the function y = x represents a perfect match while the dashed lines represent the 20% error area.About 18% of the results are outside this area.

27 Figure 11 .
Figure 11.Relationship between the predicted value and the experimental value of the pull-off strength fh for ANN test, train, and validation process for the GD algorithm with two layers of hidden neurons.

Figure 11 .
Figure 11.Relationship between the predicted value and the experimental value of the pull-off strength f h for ANN test, train, and validation process for the GD algorithm with two layers of hidden neurons.

27 Figure 12 .
Figure 12.Relationship between the predicted value and the experimental value of the pull-off strength fh for ANN test, train, and validation process for the BFGS algorithm with two layers of hidden neurons.

Figure 12 .
Figure 12.Relationship between the predicted value and the experimental value of the pull-off strength f h for ANN test, train, and validation process for the BFGS algorithm with two layers of hidden neurons.

Figure 13 .
Figure 13.Relationship between the predicted value and the experimental value of the pull-off strength fh for ANN test, train, and validation process for the LM algorithm with two layers of hidden neurons.

Figure 13 .
Figure 13.Relationship between the predicted value and the experimental value of the pull-off strength f h for ANN test, train, and validation process for the LM algorithm with two layers of hidden neurons.

27 Figure 14 .
Figure 14.Relationship between the predicted value and the experimental value of the pull-off strength fh for ANN test, train, and validation process for the CG algorithm with two layers of hidden neurons.

Figure 14 .
Figure 14.Relationship between the predicted value and the experimental value of the pull-off strength f h for ANN test, train, and validation process for the CG algorithm with two layers of hidden neurons.

Figure 15 .
Figure 15.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of R 2 .

Figure 16 .
Figure 16.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of NRMSE.

Figure 15 .
Figure 15.Comparison of models for predicting the pull-off strength of the surface layer f h due to the value of R 2 .

Figure 15 .
Figure 15.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of R 2 .

Figure 16 .
Figure 16.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of NRMSE.

Figure 16 .Figure 17 .
Figure 16.Comparison of models for predicting the pull-off strength of the surface layer f h due to the value of NRMSE.Appl.Sci.2023, 13, x FOR PEER REVIEW 20 of 27

Figure 17 .Figure 17 .
Figure 17.Comparison of models for predicting the pull-off strength of the surface layer f h due to the value of MAE.

Figure 18 .
Figure 18.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of MAPE.

Figure 19 .
Figure 19.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of average relative error.

Figure 18 .
Figure 18.Comparison of models for predicting the pull-off strength of the surface layer f h due to the value of MAPE.

Figure 17 .
Figure 17.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of MAE.

Figure 18 .
Figure 18.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of MAPE.

Figure 19 .
Figure 19.Comparison of models for predicting the pull-off strength of the surface layer fh due to the value of average relative error.

Figure 19 .
Figure 19.Comparison of models for predicting the pull-off strength of the surface layer f h due to the value of average relative error.

Table 1 .
Compositions of individual floor mortar mixes.

Table 2 .
Range of input and output variables.

Table 3 .
Correlation matrix of the input and output variables.

Table 2 .
Range of input and output variables.

Table 3 .
Correlation matrix of the input and output variables.

Table 4 .
Comparison of accuracy of compared neural networks.

Informed Consent Statement: Not applicable. Data Availability Statement:
Data available upon request.The authors declare no conflict of interest.

Table A1 .
Data collected during studies.