Microsoft Word - numero_56_art_17_2991 K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 203 Experimental study and modeling of the mechanical behavior of recycled aggregates-based high-strength concrete Kheira Camellia Nehar Laboratory of Mechanics and Materials Development at the University of Djelfa, POBox 3117, 17000 Djelfa - Algeria camellia90@hotmail.fr , c.nehar@univ-djelfa.dz , https://orcid.org/0000-0002-9188-2496 Dalila Benamara Laboratory of Civil Engineering at the University of Djelfa, POBox 3117, 17000 Djelfa - Algeria benamaradalila2018@gmail.com , d.benamra@univ-djelfa.dz ABSTRACT. The concrete constitutes a material of prime importance which is widely used in different works and infrastructures. Among the raw materials component the concrete are the aggregates which take the first place. It is widely acknowledged that the consumption of natural aggregates increases with the growth in the amount of concrete needed. It has recently emerged that concrete waste can be recycled and reused in civil engineering works after a series of treatments. Moreover, in order to protect the environment and based on the principles of sustainable development, it was considered urgent to produce a High Strength Concrete incorporating recycled materials in addition to silica fume and a high-efficiency superplasticizer. This would certainly help to establish a harmonious sustainable development that guaranties the ecological balance and environmental protection, and prevents the depletion of natural resources. This study is part of a larger research program that that seeks to recover, recycle and valorize construction and demolition wastes. The main objective sought in this article is firstly to use aggregates from demolition concrete in the manufacture of a new concrete with high mechanical and rheological performance, and secondly, to model the behavior of this type of concrete using the Finite Element Method. This modeling aims to evaluate the maximum compressive strengths and compare them with those obtained experimentally. KEYWORDS. High strength concrete; Recycled aggregates; Environment; Modelization; Finite element method. Citation: Nehar, K.C.., Benamara D., Experimental study and modeling of the mechanical behavior of recycled aggregates- based high-strength concrete, Frattura ed Integrità Strutturale, 56 (2021) 203-216. Received: 19.01.2021 Accepted: 18.03.2021 Published: 01.04.2021 Copyright: © 2021 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. https://youtu.be/bsKlrmifwEA K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 204 INTRODUCTION he refurbishment and reconstruction of civil engineering works produce enormous amounts of waste every year. Algeria, like many other countries around the world, produces between 20 and 30 million tons of demolition debris per year [1]. Considering the accelerated depletion of building material resources, in addition to the new economic, environmental, and political requirements, it has been necessary to envisage recycling waste building materials with a view to using them as substitute materials to manufacture a new type of concrete. In general, many substitute aggregates have been developed through recycling throughout the world in order to preserve the natural aggregates which are constantly decreasing. The first quality benchmark in this field dates from the Second World War (Nixon, 1977). The two most important and interactive points in these studies are the fields of application of recycled concrete waste and their classification with respect to their quality. Recycled materials, sometimes called secondary materials, are numerous and this is why it has become urgent to classify them according to their characteristics [2]. Various types of recyclable materials are currently used in civil engineering applications. These include, recycled concrete, fly and bottom ash, blast- furnace slag, steel slag, cement kiln dust, silica fume, crushed glass, reclaimed asphalt pavement, and rice husk ash. Reutilization of these recyclable materials is especially beneficial in civil engineering applications that require large volumes of materials. When these waste products are used in place of other conventional materials, natural resources and energy are preserved and expensive and/or potentially harmful waste disposal is avoided. It has been established that it is no longer required to use high quality elements to manufacture strong and efficient materials. Nowadays, more and more attention is paid to the choice of materials to be used in the design and construction of these structures. Furthermore, it is no longer necessary to use high quality materials to make a strong and resistant structure, because safety and economic factors do not have the same importance for the different types of structures. Indeed, economic factors (i.e. Supply, demand and market saturation) beside the concept of recycling, which consists of avoiding the exhaustion of supply sources, are now inciting us to revalorize recyclable materials of similar sectors [3]. The previously mentioned factors tend to increasingly encourage the recovery of construction and demolition waste to be recycled in the form of substitute aggregates for the development of new concretes. Therefore, using recycled products is perfectly justified. In this context, the right approach consists in using the right material for the right purpose under the right conditions [4]. In construction, the demand annually amounts to more than 30 million tons of cubic meters of aggregates. Furthermore, much of the waste from construction sites is increasingly transformed into aggregates to be utilized in the formulation of concrete. It can be predicted that in the future about 4 ‰ of buildings will be destroyed each year, which is expected to generate roughly 1 million of cubic meters of concrete waste [5]. In Algeria, the entire production of concrete in construction and public works is made up only of natural aggregates. Over the past few years, constraints related to the growth in demand for aggregates have surfaced, which incited researchers and practitioners to design and develop new aggregates as alternatives to replace natural aggregates. Moreover, it is worth indicating that demolition or construction materials are being increasingly produced in very large quantities. These wastes are generally used as backfill or discharged to the environment which engenders a very bad impact on the ecosystem [6]. As these new aggregates can be recycled on site, all the expenses and costs associated with the transport of natural aggregates on the one hand, and the disposal of demolition waste, on the other, can be saved, which represents a considerable economic advantage. Consequently, this procedure would contribute to reducing the depletion of natural resources and enhancing the concept of sustainable development. The valorization and the use of aggregates issues from concrete for the manufacturing of a resistant and an efficient hydraulic concrete would help to extend their current field. To the best of our knowledge, few studies have been carried out to date on the possibility of recycling construction demolition waste in order to formulate a high-strength concrete (HSC). Moreover, not many researchers have tried to numerically model this type of concrete. [7]. For these reasons, we have endeavored to explore the above-mentioned subject in order to provide a deeper insight into the topic of construction waste recovery. According to some authors, Zhang et al. [8], and Abid et al. [9], the effect of bentonite on mechanical and durability behaviour of both natural aggregates concrete (NAC) and recycled aggregates concrete (RAC) was investigated and compared. Results of testing revealed that with the incorporation of bentonite, RAC showed significant improvement in the durability and strength. To understand the degradation of mechanical properties of recycled aggregate concrete (RAC) after acid rain exposure, the effect of recycled coarse aggregate (RCA) replacement ratios, the acidity (pH value) of the simulated acid rain solution and mineral admixtures on compressive strength and elastic modulus of RAC are investigated in the study of Lu et al. [10]. T K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 205 However, several researchers have conducted studies on in the numerical modeling of concrete sections in recent years. Deshpande et al. [11] presented the findings of a study which carried out for modeled compressive strength of concrete using the techniques of Artificial Neural Network (ANN), Model Tree (MT) and Non-linear Regression. The prediction of compressive strength of concrete by data-driven models was used by Khademi et al. [12]. Nitka et Tejchman [13] studied the modelling of the behaviour of plain concrete during uniaxial compression and tension using the discrete element method DEM. An advanced finite element model using the general-purpose FE software ANSYS to simulate the behaviour of the initially damaged concrete beams was studied by Berezoug et al. [14]. Peng et al. [15] have studied the numerical simulation of recycled concrete using convex aggregate model and Base Force Element Method (BFEM), they show this combination can be used for simulating the relationship between microstructure and mechanical properties concrete. The numerical modeling of a circular concrete columns has presented by Havlásek [16] he used the nonlinear finite element simulations which allow for arbitrary loading configuration and the topology into the design process. All these research works have conducted the studies on a numerical modeling of concrete behaviour using various approaches based on software simulations. However, this modeling can be done by using the programming of the finite element method (FEM) theory [17, 18]. This method is now widely used in the engineering field and which is adopted in this work to model the recycled aggregate concrete [19]. The objective of the present work is, on the one hand, to supply an economical resistant concrete, based on recycled aggregates, having very high mechanical, physical and elastic properties, by adding silica fume and a superplasticizer and, on the other hand, to evaluate and compare the maximum compressive strengths using the finite element method. To achieve these goals, it was decided to divide the work into four stages: • The first step consisted in optimizing the high-strength concrete (HSC) with 100% natural aggregates; this was used as reference concrete, • The second step aimed to optimize the recycled HSC while varying the contents of recycled and natural aggregates, • The third step was about determining an optimal formulation that could give the maximum strength, • The last step consisted in comparing the strength obtained with that provided by the finite element method; it also sought to identify the distribution of stresses in the various samples of concrete used. USED MATERIALS Aggregates wo types of aggregates were used in this study:  Natural aggregates (NA) Class 3/8 and 8/15 crushed limestone gravels were used. These gravels have a specific gravity equal to 2.65, with 3.70% impurity and 0.30% humidity with, and a Los Angeles coefficient of 19.20%. As for sand, it has a specific weight of 2.70 and a fineness modulus of 3.20.  Recycled Aggregates The recycled aggregates (RA) were prepared in the Civil Engineering Laboratory at the University of Djelfa, in Algeria. The following table shows the chemical characteristics of the aggregates. Chemical characteristics P.A.F. % SiO2 % AI2O3 % Fe2O3 % CaO % MgO % SO3 % K2O % Na2O % NA 1.47 11.35 4.59 3.52 78.89 1.10 1.22 0.41 0.13 RA - 47.17 -0.63 0.85 49 1.02 1.30 0.50 - Table 1: Chemical characteristics of aggregates. Tab. 1 above summarizes the chemical composition of the different aggregates. It was found that the recycled sand included silica and limestone (47% of silica and 49% of carbonate) with a small amount of Alumina (Al2O3) which was due to the presence of cement paste around the sand grains. However, natural aggregates were calcareous in nature, with a fairly high percentage of carbonates (around 78%). Therefore, it can be said that in general recycled aggregates are of a lime-silicate nature. However, this is not always true because it can vary according to the type of the recycled waste. For example, some researchers [20, 21] found out that their recycled aggregates had a rather siliceous nature. T K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 206  Physical and mechanical characteristics: The water absorption test was carried out according to Standard NF EN 18-554, and the results obtained are summarized in the table below. Type of Aggregate Los Angeles (%) Water absorption (%) Sand equivalent (%) Absolute density Apparent density NA RA NA RA NA RA NA RA NA RA 8/15 19.20 31 1.35 1.1 2.48 2.45 - - 0.78 5.20 3/8 21 33 1.38 1.21 2.57 2.60 - - 0.92 6.30 0/5 - - 1.20 0.9 2.89 2.80 80 76 - - Table 2: Physical and mechanical characteristics of aggregates. The findings revealed that the densities decreased as the granular fraction augmented. In addition, natural aggregates and recycled aggregates exhibited nearly the same absolute densities. On the other hand, the apparent densities of recycled aggregates were found lower than those of natural aggregates due to the presence of pores in the recycled aggregates because the concrete waste had a high porosity. Furthermore, the water absorption test helped to see the difference between the two types of aggregates from the density viewpoint, since the recycled aggregates showed a greater water absorption coefficient (approximately 8%) than that of natural aggregates (about 2%). This difference is primarily attributed to the fact that recycled aggregates had a higher porosity than that of natural aggregates [4, 6]. Cement The cement (C) used is a Portland cement with high resistance to sulphates (CEM I 42.5 HRS). Its mineralogical composition is such that C3S = 60%, C2S = 25.55%, C3A = 2% and C4AF = 12.45%. The density of anhydrous cement is 3.23 g / cm3 and its specific surface is 3880 (cm2 / g). Adjuvant The adjuvant used is a new generation high-water-reducing superplasticizer of the MEDAFLOW 145 type. It is a solution of polycarboxylates, 31% dry extract, light pitch color, and pH between 5 and 6. Silica fume The actual density of silica fume is 2194 (Kg/m3) and its Blaine specific surface is greater than 15000 (cm²/g). The chemical compositions of cement and silica fume (SF) used in this study are given in Tab. 3. The samples P.A.F. % SiO2 % AI2O3 % Fe2O3 % CaO % MgO % SO3 % K2O % Na2O % Cement (C) 1.47 22.31 4.39 5.52 63.89 0.80 1.11 0.41 0.09 Silica fume (SF) - 93.17 0.60 1.25 1.40 1.02 2.30 1.00 - Table 3: Chemical compositions of the constituents. The Blaine specific surface area of cement is SSB = 3304 (cm² / gr). The Blaine Specific Surface Area of silica fume is SSB > 20000 (cm² / gr). The abbreviations used: NC: Normal Concrete RC: Recycled Concrete C (N + R): Normal Concrete + Recycled Concrete HSC: High Strength Concrete A: Aggregates NA: Natural Aggregates RA: Recycled Aggregates SF: Silica Fume K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 207 Adj: Adjuvant C: Cement (W/C): (Water/Cement) ratio (S/G): (Sand/Gravel) ratio (W/B): (Water/Binder) ratio CS 28D: Compressive strength at 28 Days. FEM: Finite Element Method. EXPERIMENTAL STUDY OF THE COMPOSITION OF CONCRETE Formulation of reference concrete enerally, the formulation methods do not all provide the same results, because each of them must simplify the hypotheses and rely only on a few of them. The assumptions are not always the same. Note, however, that none of these hypotheses leads with certainty to an accurate composition. This is certainly due to the fact that the numerous parameters which have an impact on concrete are difficult to evaluate. With regard to the formulation of the reference concrete, it was decided to use the Baron-Lesage method, but the aggregate content should be known in advance in order to have an optimal mineral skeleton. According to Baron and Lesage [22], the following two assumptions must be respected: • For a fixed dosage of cement and water at the beginning (fixed ratio W/C), the most workable and compact concrete is the one that is the most compact with an optimal skeleton. • The optimal dosage of aggregates does not depend on the contents of water and cement. It is worth noting that in order to achieve the optimal composition of concrete, several parameters, such as the cement content, water content and aggregate content (S/G = Sand/Gravel) should be considered. Therefore, in order to simplify our study and satisfy the hypotheses of the method, only the (S/G) ratio was varied. The cement dosage and the ratio (W/C) were kept unchanged at the beginning. According to Gorisse [23], the optimal ratio (S/G) is between 0.50 and 0.83. In the present work, seven values of this ratio were selected and tested (step: 0.05) for the purpose of obtaining seven different compositions. The Abrams cone sag, corresponding to the values of the ratio (S/G) = 0.50, 0.55, 0.60, 0.65, 0.70, 0.75, and 0.80, was measured for each composition. Types of concrete C (kg/m3) A (kg/m3) W (L/m3) W/B 8/15 3/8 0/5 NC (100%) 350 905 200 670 192.5 0.55 RC (100%) 350 905 200 670 245 0.70 (70%NA+30%RA) 350 905 200 670 210 0.60 (50%NA+50%RA) 350 905 200 670 227.5 0.65 Table 4: Optimal composition of reference concrete. The compressive strength at 28 days with these types of concrete are presented in the Fig. 1. Formulation of the concretes under study Once the reference concrete was prepared, a formulation method was established in order to find the optimal composition of the high performance concretes (HPC). It should be noted that the formulation method adopted for the mixture of constituents entering into the composition of high performance concrete is the one that was developed by Aïtcin [2] at the University of Sherbrooke. This is a simple and practical method. It consists in combining the empirical results and those obtained using the absolute volume method. The amount of water contained in the superplasticizer is considered to be part of the amount of mixing water. The data used for the optimal formulation which gave the maximum strength with acceptable workability for one cubic meter of concrete are presented in the Tab. 5. G K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 208 Figure 1: The compressive strength at 28 days of reference concrete. Types of concrete C (kg/m3) A (kg/m3) SF (%) Adj. (%) W/B Sagging (cm) Fresh density8/15 3/8 0/5 NC+SF 450.70 950 200 670 60.2 10.5 0.30 5.0 2567 RC+SF 522.46 950 200 670 70.4 14.2 0.37 5.7 2389 (70%NA+30%RA+SF) 472.13 950 200 670 65.6 11.7 0.32 4 2465 (50%NA+50%RA+SF) 484.60 950 200 670 72.3 13.4 0.35 3 2332 Table 5: Composition of concrete mixes. Figure 2: The compressive strength at 28 days of concrete mixes. Fig. 2 shows the compressive strength at 28 days obtained with this optimal formulation. N C (1 00 % ) RC (1 00 % ) (7 0% N A+ 30 % RA ) (5 0% N A+ 50 % RA ) 0 10 20 30 40 C o m p re ss iv e st re n gt h a t 28 d ay s (M P a) N C+ SF RC +S F (7 0% N A+ 30 % RA +S F) (5 0% N A+ 50 % RA +S F) 0 10 20 30 40 50 60 70 80 90 C o m p re ss iv e st re n gt h a t 28 d ay s (M P a) K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 209 Figure 3: Flowchart of the numerical simulation model. Discussion and interpretation The table above shows that: - It is possible to manufacture high strength concretes with a compressive strength exceeding 70 MPa (at 28 days) and adequate plasticity, provided that the recycled aggregates are used up to 30% of the total amount of aggregates in ordinary high-strength concrete. - Partial substitution of natural aggregates with recycled aggregates for the manufacture of the new concrete offers a new aggregate resource and allows saving natural materials. - The findings were encouraging since the strength drop of concrete incorporating recycled gravel only does not exceed 17% for a target characteristic strength of 70 MPa. These results confirm the trend evoked by the antecedent studies of Revilla-Cuesta et al. [25], and Masood et al. [26], their researches is presented into recycled concrete aggregate properties and the mix-design of the self-compacting concretes that contain them, as well as relevant results on the fresh state (workability, rheology), the hardened state (compressive strength, splitting tensile and flexural strength, modulus of elasticity, density, and porosity), durability (resistance to aggressive agents), long-term properties of concrete (shrinkage, creep), and structural elements manufactured with self- - Mesh construction. - Introduction of applied loading. - Form the force vector {F} - Computation of displacement {U} = [K] /{F} - Computation of stresses - Computation of the elementary stiffness matrix  [ ] [ ] [ ][ ] T e V K B C B dV - Assembly: Construction of the global stiffness matrix [K]. - Introduction of boundary conditions: Embedding, single or double support depending on the study case. BEGIN Data input: Geometric: - Geometry dimension b, h, ep. - Number of elements in X and Y direction Mechanics: - Young's modulus E. - Poisson's ratio ν. TREATMENT END -Initialization of the elementary stiffness matrix [Ke]. - Computation of interpolation functions [N]. - Computation of the derivatives of interpolation functions [dN]. - Computation of the matrix [B]. - Computation of the elasticity matrix [C]. K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 210 compacting concrete containing recycled concrete aggregate. The results under review reaffirm that the incorporation of recycled concrete aggregate can produce a suitable self-compacting recycled concrete, on the basis of careful designs that are essential for successful performance. NUMERICAL STUDY OF CONCRETE SAMPLES Modeling by the finite element method he finite element approximation was programmed under MATLAB [27] for the finite element modeling of concrete specimens. It is represented by the organization chart in Fig. 3. Numerical application The samples investigated in the present work were modeled and evaluated using our finite element calculation code. These are concrete elements of width b = 10 cm, height H = 10 cm and thickness ep = 10 cm, as shown in Fig. 4 (a), and plane strain conditions are adopted. In the design procedure, it was decided to apply a compression loading similar to that used in the experimental part with a quadrangle mesh density of 100 x 100 elements for all types of concrete (Fig. 4 (b)). Figure 4: (a) Geometric configuration of the concrete samples studied, (b) Modeled sample. Several types of concrete were prepared and investigated in this study. Their mechanical characteristics are summarized in the following table: Types of concrete Young’s Modulus E(MPa) Poisson's ratio ν NC (100%) 29000 0.20 RC (100%) 24000 0.19 (70%NA+30%RA) 27000 0.20 (50%NA+50%RA) 25000 0.20 (NA+ SF) 46000 0.25 (RA +SF) 38000 0.23 (70%NA+30%RA+SF) 42000 0.24 (50%NA+50%RA+SF) 39000 0.22 Table 6: Mechanical properties of concretes used in this study. Results and discussion The finite element method (FEM) was programmed and applied to the previously prepared specimens. The results obtained were compared with those provided by the experimental study for the purpose of highlighting the performance of the developed program. T (b) (a) K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 211 Fig. 5 displays the results of the comparison between the experimental compressive strengths and those predicted numerically by the finite element method (FEM), at 28 days, for all types of concrete studied. Figure 5: Experimental and numerical compressive strength values at 28 days. It is quite clear that the compressive strength values obtained at 28 days using our finite element calculation code are in good agreement with those resulting from the experimental tests. Consequently, the data obtained for the compressive strength are quite satisfactory and encouraging. In order to better understand the distribution of compressive stresses in our concrete samples, the evolution of these stresses was presented for each type of concrete as follows: (NC (100%), RC (100%), (70% NA + 30% RA), (50% NA + 50% RA), (NA + SF), (RA + SF), (70% NA + 30% RA + SF), (50% NA + 50% RA + SF)), and which are shown in Figs. 6 (a), (b), (c), (d), (e), (f), (g), (h), respectively. Figs. 6 (a), (b), (c), (d), (e), (f), (g) and (h) present the mapping of the equivalent stress field resulting from numerical simulations using the finite element method (FEM). The distribution of stress is symmetrical about the loading application point. The stress field is maximal at the center of the test specimen. However, the mapping shows the presence of stress concentrations (circled in black in Fig. 6) in the loading application regions which are located at the edge of the specimens. These concentrations were in fact expected due to the presence of the compressive loading at the edge of the specimens. Furthermore, this mapping also showed that the addition of 30% natural aggregates along with silica fume in the concrete makes it possible to achieve a better distribution of stresses compared to normal concrete. CONCLUSIONS he experimental results obtained indicate that: - It is quite possible to manufacture a highly resistant concrete, by the partial substitution of natural aggregates with recycled aggregates (70% NA + 30% RA), with a compressive strength of about 70.80 MPa at 28 days. This value is not too far from that of concrete incorporating natural aggregates only (100% NA); its compressive strength was found approximately equal to 83 MPa for the same age. These findings are attributed to the fact that the super plasticizer was very effective in reducing the water-to-binder (W/B) ratio below 0.30. In addition, the use of additions, such as the finely ground silica fume, considerably helped to make the cement paste denser and to further reduce the amount of water, especially that the cement used contains a very small amount of Tricalcium aluminate (CA3). - Recycled aggregates could be used up to 30% of the total amount of aggregates for most environments. If recycled sand is used, then the amount used should be less than 30%. T NC (10 0% ) RC (10 0% ) (70 %N A+ 30% RA ) (50 %N A+ 50% RA ) (N A+ SF ) (RA +S F) (70 %N A+ 30% RA +S F) (50 %N A+ 50% RA +S F) 0 10 20 30 40 50 60 70 80 90 C o m p re ss iv e st re n gt h a t 28 d ay s (M P a) Experimental results Numerical results K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 212 Figure 6: Distribution of compressive stresses in the samples using MATLAB [27]: (a) NC (100%), (b) RC (100%), (c) (70% NA+ 30% RA), (d) (50% NA + 50 % RA), (e) (NA+ SF), (f) (RA + SF), (g) (70% NA+ 30% RA + SF), and (h) (50% NA+ 50 % RA + SF). (a) (b) (c) (d) (e) (f) (g) (h) K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 213 - Moreover, good workability can easily be achieved despite the low ratio (W/B) = 0.35 thanks to the good quality superplasticizer (Media Flow 145) used. It is useful to remember that this superplasticizer has fairly high rheological and mechanical performance. - The use of recycled aggregates is of paramount importance from the environmental point of view. This allows getting rid of the materials resulting from the demolition of old buildings and the destruction of structures following natural disasters. - Moreover, this study allowed us to model the concrete samples prepared, using the finite element method. - A simple comparison helped to confirm the good correlation between the results obtained by our calculation code and those obtained by the experimental study, for all types of concrete. It can therefore be concluded that the addition of 30% of recycled aggregates and silica fume into concrete allows achieving a better distribution of stresses in comparison with ordinary concrete. REFERENCES [1] Ministère de L’aménagement du Territoire et de L’Environnement (MATE), Algeria. [2] Tayeh, B. A., Al Saffar, D. M. and Alyousef, R. (2020). The utilization of recycled aggregate in high performance concrete: a review. J. of Mat. Res. and Tech., 9(4), pp. 8469–8481. DOI: 10.1016/j.jmrt.2020.05.126. [3] Ali, A. A. M., Zidan, R. S. and Ahmed, T. W. (2020). Evaluation of high-strength concrete made with recycled aggregate under effect of well water. Case Stu. in Const. Mat., 12, e00338. DOI:10.1016/j.cscm. 2020.e00338. [4] Jin, R., Li, B., Elamin, A., Wang, S., Tsioulou, O. and Wanatowski, D. (2018). Experimental investigation of properties of concrete containing recycled construction wastes. Int. J. of Civ. Eng., 16(11), pp. 1621–1633. DOI: 10.1007/s40999-018-0301-4. [5] Guo, H., Shi, C., Guan, X., Zhu, J., Ding, Y., Ling, T. C., Zhang, H. and Wang, Y. (2018). Durability of recycled aggregate concrete–a review. Cem. and conc. comp., 89, pp. 251–259. DOI: 10.1016/j.cemconcomp.2018.03.008. [6] Verian, K. P., Ashraf, W. and Cao, Y. (2018). Properties of recycled concrete aggregate and their influence in new concrete production. Res. Conserv. and Recyc., 133, pp. 30–49. DOI:10.1016/j.resconrec.2018.02.005. [7] Siempu, R. and Pancharthi, R. K. (2017). Bond characteristics of concrete made of recycled aggregates from building demolition waste. Mag. of Conc. Res., 69(13), pp. 665–682. DOI:10.1680/jmacr.16.00400. [8] Zhang, Y., Luo, W., Wang, J., Wang, Y., Xu, Y., & Xiao, J. (2019). A review of life cycle assessment of recycled aggregate concrete. Const. and Build. Mat., 209, pp. 115-125. DOI: 10.1016/j.conbuildmat.2019.03.078. [9] Abid, S. R., Nahhab, A. H., Al-aayedi, H. K. and Nuhair, A. M. (2018). Expansion and strength properties of concrete containing contaminated recycled concrete aggregate. Case Stud. in Constr. Mat., 9, pp. e00201. DOI: 10.1016/j.cscm.2018.e00201. [10] Lu, C., Wang, W., Zhou, Q., Wei, S. and Wang, C. (2020). Mechanical behavior degradation of recycled aggregate concrete after simulated acid rain spraying. Jour. of Clea. Produc., 262, pp. 121237. DOI: 10.1016/j.jclepro.2020.121237. [11] Deshpande, N., Londhe, S., and Kulkarni, S. (2014). Modeling compressive strength of recycled aggregate concrete by Artificial Neural Network, Model Tree and Non-linear Regression. Int. J. of Sust. Bui. Env., 3(2), pp. 187-198. DOI: 10.1016/j.ijsbe.2014.12.002. [12] Khademi, F., Akbari, M., and Jamal, S. M. M. (2015). Prediction of compressive strength of concrete by data-driven models. I-Manager’s J Civ Eng, 5, pp. 16. [13] Nitka, M. and Tejchman, J. (2015). Modelling of concrete behaviour in uniaxial compression and tension with DEM. Granu. Mat., 17(1), pp. 145–164. DOI: 10.1007/s10035-015-0546-4. [14] Berrezoug Zidani, M., Belakhdar, K. and Tounsi, A., Adda Bedia, E. (2015). Finite element analysis of initially damaged beams repaired with FRP plates. Compo. struc., 134, pp. 429–439. DOI: 10.1016/j.compstruct.2015.07.124. [15] Peng, Y., Chu, H., and Pu, J. (2016). Numerical simulation of recycled concrete using convex aggregate model and base force element method. Adv. in Mat. Scien. and Eng., 2016. DOI: 10.1155/2016/5075109. [16] Havlásek, P. (2020). Numerical modeling of axially compressed circular concrete columns. Eng. Struc., 227, pp. 111445. DOI: 10.1016/j.engstruct.2020.111445. [17] Pradhan, K. K. and Chakraverty, S. (2019). Chapter Four—Finite element method. In Computational structural mechanics, Academic Press, pp. 25–28. K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 214 [18] Baaser, H. (2010). Development and Application of the Finite Element Method based on MatLab. Springer Science & Business Media. [19] Tekkaya, A.E. and Soyarslan, C. (2018). Finite Element Method. In: Int. Academy for Production Engineering, Chatti S., Tolio T. (eds) CIRP Encyc. of Prod. Eng. Springer, Berlin, Heidelberg. DOI:10.1007/978-3-642-35950-7_16699-3 [20] Xiao, J. (2018). Recycled aggregate concrete structures, Springer-Verlag Berlin Heidelberg, pp. 65–98. DOI: 10.1007/978-3-662-53987-3. [21] Babu, V. S., Mullick, A. K., Jain, K. K. and Singh, P. K. (2015). Strength and durability characteristics of high-strength concrete with recycled aggregate-influence of processing. J. of Sust. Cem.-Bas. Mat., 4(1), pp. 54-71. DOI: 10.1080/21650373.2014.976777. [22] Baron, J. (1997). Les additions normalisées pour béton, Ecole française du béton, édition Eyrolles. [23] Gorisse, D., Festa, J. (1998). Nouveau guide du béton et de ses constituants, Huitième édition, Edition Eyrolles, France. [24] Aïtcin, P.C. (2001). Bétons haute performance, Edition Eyrolles, Paris. [25] Revilla-Cuesta, V., Skaf, M., Faleschini, F., Manso, J. M. and Ortega-López, V. (2020). Self-compacting concrete manufactured with recycled concrete aggregate: An overview. Jour. of Clea. Produc., 262, pp. 121362. DOI: 10.1016/j.jclepro.2020.121362. [26] Masood, B., Elahi, A., Barbhuiya, S. and Ali, B. (2020). Mechanical and durability performance of recycled aggregate concrete incorporating low calcium bentonite. Constru. and Build. Mat., 237, pp. 117760. DOI: 10.1016/j.conbuildmat.2019.117760. [27] MATLAB (1984). Script language with a development environmental for numerical calcul, Version 16.0, MathWorks, California, United States. ANNEX FEM Modeling In order to model the concrete samples, the finite element approximation is used [12]. We consider a quadrilateral element (Fig. A.1) the element is defined by four nodes in natural coordinates (ξ, η). The coordinates are interpolated and given as follows [15]:      4 4 1 1 i i i i i i u N u v N v (1) Where u, v represents the displacements fields at any point of the element and ui, vi; the nodal displacements. Figure A.1. Quadrilateral Q4 element in natural coordinates (ξ, η). Ni represents the standard interpolation functions which are given by:         1 1 1 1 ( , ) ( ) ( ) (1 )(1 ) 4 N l l (2)         2 2 1 1 ( , ) ( ) ( ) (1 )(1 ) 4 N l l (3) K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 215         3 2 2 1 ( , ) ( ) ( ) (1 )(1 ) 4 N l l (4)         4 1 2 1 ( , ) ( ) ( ) (1 )(1 ) 4 N l l (5) In matrix form:                          (1 )(1 ) (1 )(1 )1 (1 )(1 )4 (1 )(1 ) T iN (6) The derivatives of the interpolation functions Ni are:                      (1 ) (1 )1 (1 )4 (1 ) T iN ,                      (1 ) (1 )1 (1 )4 (1 ) T iN (7) These derivatives are transformed by a function F(x, y) with respect to ξ and η:                                                                                               [ ] yF yF F x x F F yx x x J F FF F x y yF x y yx y (8) where J is the Jacobian operator, connecting the natural and global coordinates. The derivatives with respect to the global coordinates can be found as follows:                              1[ ] FF x J F F y (9) The stiffness matrix is then obtained by:  [ ] [ ] [ ][ ] T V K B C B dV (10) where [B] is the strain matrix derived from the interpolation functions:                                                   1 0 [ ] 0 , [ ] i ii i i i i i N NNx N x B J N Ny yN N x y (11) K.C. Nehar et alii, Frattura ed Integrità Strutturale, 56 (2021) 203-216; DOI: 10.3221/IGF-ESIS.56.17 216 [C] is the elasticity matrix:              2 1 0 [ ] 1 0 1 0 0 2(1 ) E C , in plane strains (12)                        1 0 [ ] 1 0 (1 ) 1 2 (1 2 )0 0 2 E C , in plane stresses (13) The volume of the element dV is given by:   .det . .dV h J d d (14) where detJ is the determinant of the Jacobian matrix and h the thickness of the element. The surface dξ.dη at the point (ξ, η) of the reference element is transformed into the surface dA at the point (x ((ξ, η)), y ((ξ, η))) of the natural element. Were:   det . .dA J d d (15) The integral in the stiffness matrix on the real element is:  [ ] [ ]. V K F dV (16) Therefore, becomes on the reference element:  [ ] .[ ]. A K h F dA (17) Taking: [ ] [ ] [ ][ ]TF B C B (18) The integral in the stiffness matrix is then calculated numerically by Gaussian quadrature in two dimensions. In the four- node element, we can use a 2 × 2 numerical integration for exact integration.                1 1 1 11 1 , . . , p q i j i j i j F d d w w F (19) where p, q are the number of integration points in the directions ξ, η, respectively, and wi, wj are the corresponding weights. The stiffness matrix can be calculated using 2 × 2 Gauss points as:                                   1 1 1 1 2 2 1 1 . . . det . . . . det . T T e A T i j i j K h B C B dA h B C B J d d h B C B J w w (20)