Energy Evolution Behavior and Mesodamage Mechanism of Crumb Rubber Concrete

School of Civil Engineering, Shandong University, Jinan, Shandong 250061, China Institute of Marine Science and Technology, Shandong University, Qingdao, Shandong 266237, China Key Laboratory of Disaster Prevention and Reduction of Civil Engineering in Shandong Province, Qingdao, Shandong 266590, China School of Civil Engineering and Architecture, Shandong University of Science and Technology, Qingdao, Shandong 266590, China


Introduction
CRC (crumb rubber concrete), referred as elastic composite material concrete, has the excellent performance on light weight, elastic damping, ductility and toughness, and environmental protection.CRC has been broadly applied in the fields of road, highway, bridge deck, and underground engineering, while it often suffers damage and failure on a larger scale due to the generation and propagation of microcracks, which are often neglected.erefore, understanding and mastering the mesodamage mechanism and failure rules of CRC materials is very important to ensure the quality and safety of the materials applied in engineering.
Current research often focuses on improving the overall strength, deformation, or impact resistance of CRC materials by designing better ratio schemes through laboratory tests, but ignores the occurrence and expansion process of material failure that often occurs in engineering.Eldin and Senouci [1], Khatib and Bayomy [2], and Yu and Zhu [3] conducted laboratory experiments to study the law of the physical and mechanical properties of composite material mixed with different sizes and different amounts of rubber particles.Segre and Joekes [4] and Wang et al. [5] studied the bonding effect of the rubber modified by chemical solution and obtained the properties of mortar-rubber interface and the physical and mechanical properties law of the aggregate concrete.Topcu [6][7][8] carried out a series of researches on waste tire concrete, in which the physical and mechanical properties of the concrete at different ages were studied in detail by preparing samples with different size and dosage of rubber particles content.Hernandez-Olivares et al. [9] prepared aggregate concrete by adding the crushed tire rubber and polypropylene fiber of a certain volume fraction into the concrete and studied the cementation and compatibility of cement-rubber through SEM (scanning electron microscope), and the tensile and bending resistance tests and dynamic compression tests were performed at different curing times of samples.Youssf et al. [10] took the effect of rubber particles' preprocessing time, silicon powder content, and cement content on the amount of compression, compression strength, and tensile strength of CRC into account, and the best mix ratio that can improve the mechanical properties of CRC was obtained.
Zhu et al. [11][12][13] conducted extensive research on the mechanical properties of crumb rubber concrete and established a regression relationship between compressive strength, elastic modulus, and influencing factors, and the failure mechanism of crumb rubber concrete was discussed.Yang [14] made a theoretical and analytical exploration of CRC from the perspective of microstructure and application.Liu et al. [15][16][17] established an aggregate model of a random distribution by discretizing the CRC at a mesoscopic level and conducted a simulation study on the deformation of CRC based on the finite element method.Liu [18] studied the frost-resistance durability and microscopic frostresistance mechanism of rubber concrete mixtures in different erosion media.Yuan and Zheng [19], Xu et al. [20], and Guo et al. [21] carried out the dynamic performance test of rubber concrete, respectively, and different impact resistance between ordinary concrete and CRC was analyzed.
However, the internal mesomechanical response processes and laws in materials cannot be revealed through conventional laboratory experiments.In addition, the limitations of laboratory conditions and the discrete physicalmechanical properties of the samples make it impossible to systematically understand and master the energy evolution rules and mesodamage mechanism of materials.
In the present study, the damage mechanism and energy evolution of CRC materials in the process of deformation and failure were studied from a microscopic perspective by referring to the particle discrete element method introduced by Cundall and Strack [22] in 1979.First, the deformation and failure characteristics of CRC with different rubber content were briefly described and analyzed, and the correlation between CRC strength, deformation parameters, and rubber content of CRC with different rubber content was obtained by data fitting method.
en, the micromechanical parameters of CRC composites were determined by inverse analysis based on the laboratory research results, and the rationality and validity of the micromechanical parameters were verified by the experimental results of other groups of samples.Next, the mechanism of generation and propagation of microcracks during CRC failure was analyzed, and the correlation between mesoscopic damage variable and dissipated energy at different stages was studied, and the evolution of CRC damage threshold with rubber content was discussed.And then, the energy release ratio is defined after the composite materials reach the peak strength, which is used to analyze and evaluate the flexibility and ductility of the composite concrete materials, and the energy release ratio of CRC with different rubber content was calculated and analyzed by data fitting.Finally, taking the CRC group with the maximum damage threshold as the representative, the essential mechanism of different failure characteristics of CRC and NC was discussed from the perspective of mesoscopic energy.is process can not only effectively solve the limitations caused by the laboratory test and the dispersion of experimental results, but also can provide useful reference and supplement for the study on mechanics, progressive deformation, and failure characteristics of composite materials.

Laboratory Experiment
is test is mainly to study the effect of different blending content of rubber particles on the mechanical properties of CRC. e following materials were used in this work: P.O 42.5 cement and polyacrylate emulsion as binders; continuous grading ordinary sand with maximum size of 3 mm; rubber particles with a size of 3-5 mm; naphthalene superplasticizer as an admixture; and saturated sodium hydroxide solution (NaOH).
e rubber particles were surface-pretreated with saturated NaOH solution at room temperature for 30 minutes while stirring.It not only removed dust particles on the rubber surface, but effectively decomposed carboxyl and acid-based components on the rubber surface [7].e adhesion of rubber-mortar and the overall strength of the material were also enhanced.Rubber particles before and after surface-pretreated with NaOH saturated solution are shown in Figure 1.
According to the orthogonal test results of our previous work [23], the amount or proportion of materials used in samples are as follows: the cement dosage was 380 kg/m 3 , the water/cement ratio was 0.35, the cement/aggregate ratio was 0.2, the water reducer dosage was 1.5% of the cement dosage, and the binder dosage was 15% of the rubber particle dosage.
ere were eight groups with different rubber content in the experiment, they were 0%, 5%, 10%, 15%, 20%, 25%, 30%, and 40% of the sample mass, and the group of 0% rubber content, that is, normal concrete, was used as the control group.Based on the mixing principle of materials, the mixing ratio of each group was determined as shown in Table 1.en, the same preparation process was used to complete the samples preparation in the same laboratory environment, and the strength tests were conducted after curing in the standard curing room for 28 days.
e test specimen was the standard cylindrical shape with a specification of Φ50 mm × 100 mm.ree specimens were prepared for each group and the average value was taken as the representative value of the corresponding mix proportion scheme.
e electrohydraulic servo compression testing machine was used to test the strength and deformation characteristics of the CRC. Figure 2 shows the loading process of CRC specimens, the loading rate is set to 0.1 mm/min.

2
Advances in Materials Science and Engineering

Experimental Results and Analysis.
e typical failure form of uniaxial compression of each group is shown in Figure 3, and it can be seen that normal concrete shows shear failure of single sloping plane and no other obvious small cracks appear on the surface.Group A and B also show shear failure of single sloping plane with large perforative crack.With the increase of rubber particles content, the failure form of the samples from C to H was gradually changed from a form of single sloping and complete plane to a form of short and dispersed cracks.And the serial macrocrack began to disappear, and many fine discontinuous microcracks began to appear on the surface.e axial deformation of samples gradually increased, and the integrity of the sample was better after reaching its ultimate bearing capacity.
e correlation between the mechanical deformation parameters and the rubber particles content in every sample group is shown in Figure 4.It shows that the compressive strength of CRC is 69.9% lower than that of NC, and CRC has a slow decrease in compressive strength with the increase of rubber particles content.
e fitting function y 1 � (−9375 * x + 1.873e5)/(x 2 + 2.572e4 * x + 3.994e4) can be used to describe the relationship between the compressive strength and the rubber particles content, with a variance of 0.9917.With the increase of rubber particles content, the elastic modulus shows a similar attenuation trend as the compressive strength, and the elastic modulus of CRC has a decrease of 56.5% compared to that of NC.Compared with the trend of change in compressive strength, the decrease in elastic modulus is smaller and the downward trend is slower in the later period.e function y 2 � (−3.691e8 * x + 5.985e8)/(x 2 + 3.692e5 * x + 6.994e5) can be used to fit and analyze the relationship between the elastic modulus and the rubber particles content, with a variance of 0.9848.

Numerical Simulation Experiment
3.1.Particle Discrete Element Method.Particle discrete element method is a particular implementation of a broader class of methods known as discrete element methods, which are defined as methods that allow finite displacements and rotations of discrete bodies, including complete detachment, and recognize new contacts automatically as the calculation progresses.
e CRC material is simulated by a dense packing of non-uniform-sized spherical particles that are bonded together at their contact points, and its mechanical behavior is described by the movement of each particle and the force and moment acting at each contact based on Newton's second law and force-displacement law.e PBM (parallel-bond model) is implemented in the PFC3D using the DEM and can be used to simulate and analyze the conventional mechanics, deformation, and failure characteristics of rock and concrete [24].e contact model and failure criterion are shown in Figure 5.

Advances in Materials Science and Engineering
As shown in Figure 5, the PBM is composed of cemented substance and adjacent rigid particles, and its mechanical model consists of several mechanical components that have axial or tangential forces along the line joining the particle centers.In the axial direction, the normal linear contact force of particles is produced by linear springs with constant normal stiffness k n , the normal dashpot force is produced by dashpots with normal viscous rheological β n , the tensile strength σ c of cement is produced by bonding block, and the normal linear contact force of cement is produced by the other linear springs with constant normal stiffness k n .In the tangential direction, the tangential linear contact force of particles is produced by linear springs with the tangential stiffness k s , the tangential dashpot force is produced by dashpots with tangential viscous rheological β s , the shearing strength c and Φ of cement is produced by shear blocks, the tangential linear contact force of cement is produced by the other linear springs with constant tangential stiffness k s , and the friction between adjacent particles is produced by sliders with friction coefficient μ.
Under the action of load on both ends of the particles, the normal stress σ and tangential stress τ of the bond in the parallel-bond model can be expressed, respectively, as where the moment-contribution factor β ∈ [0, 1], R is the radius of bond, A is the cross-sectional area of bond, I is the inertia moment of bonding cross section, and J is the polar inertia moment of the cross section of bond.e force acting on the bond can be resolved into the normal force F n and the tangent force F s , and the moment of bond can be resolved into the torsional moment M t and the bending moment M b .

Numerical Simulation Scheme.
In order to better simulate the mechanics and contact behavior between particles in the CRC, the three-dimensional numerical models with different rubber particles content were established according to the size of the laboratory samples, and its specification was Φ50 mm × 100 mm.e yellow particles in models represent the sands with the minimum radius of 0.8 mm, and the black particles represent the rubber particles with the minimum radius of 1.2 mm. e particle diameter ratio is 1.66 [25].e eight numerical models established by particle discrete element method are shown in Figure 6, in which the rubber mass fraction is 0, 5%, 10%, 15%, 20%, 25%, 30%, and 40%, respectively.
In order to eliminate the "cyclo-hoop effect" at both ends of the sample, the friction coefficient between the loading plate and the specimen is set to 0 during the numerical experiment to keep consistent with the sample boundary conditions in laboratory tests.e compiled code using the FISH language was used to achieve the uniaxial compression numerical experiment of the samples.e microdistribution of the particles and bond morphology in CRC materials are shown in Figure 7, in which the black parts represent the rubber particles, the yellow parts represent the sand particles and the white parts represent the bond formed by cement and coagulant.

Reverse Analysis and Determination of Mesoparameters.
Based on the results of a large number of numerical experiments, combined with conclusions of other scholars' papers and monography [26][27][28], the mechanical mechanism of PBM and calibration basis of mesoscopic parameters can be summarized as follows: the macroelastic modulus is mainly related to the contact stiffness k n of particle surfaces and the stiffness k n of parallel bond, and the Poisson ratio is related to the normal/tangent stiffness ratio k n /k s of adjacent particles and the normal/tangent stiffness ratio k n /k s of bond.And the peak strength is related to the shearing strength c, Φ of bond, and the friction coefficient μ between adjacent particles.e inversion calibration of the CRC mesoparameters was conducted based on the results of laboratory experiments.e materials and processing technologies used in the preparation process of laboratory specimen are the same, and the only difference among them lies in the content of rubber particles; thus, the mesoparameters of particles and bonds in materials are the same in each sample.Firstly, the mesoparameters inversion calibration of sand particles based on the failure pattern and mechanical behavior of the NC was implemented, and then based on the established mesoparameters of the sand particles as well as the failure morphology and mechanical behavior of the sample with 40% rubber contents, the mesoparameters of rubber particles are established.e calibration curve of constitutive relation of materials is shown in Figure 8.
By combining the results of numerical tests and indoor tests and comparing the macromechanical parameters (Table 2) and failure modes (Figure 9) of the two materials, it can be concluded that the numerical model and experiment can ideally simulate and reveal the mechanical properties of the NC or the CRC, the established mesoparameters can reflect the mesomechanics and interaction of the internal particles within the composite material, and the mesoscopic physical and mechanical parameters of materials are shown in Table 3.
In order to verify the rationality and validity of the determined mesomechanical parameters, numerical simulations were conducted on other CRC samples based on the determined mesomechanical parameters, and the constitutive curves of each group were obtained.e comparison of constitutive curves obtained by laboratory experiments and those obtained by numerical simulations is shown in Figure 10.It can be found from the comparative analysis of the curve in the figure that the constitutive curves obtained from numerical experimentation are basically consistent with the laboratory test results, and the characteristics are as follows: with the increase of the rubber particles content, the compressive strength of the material shows a decreasing trend, the results are basically the same as that from laboratory tests, and the peak strain increases accordingly.It can be seen from the constitutive curve that the peak area of the stress-strain curve is gradually broaden with the increase of rubber particle content and the curve descending rate after the peak gradually become gentle.Since spherical particles are used in numerical simulation and the accidental factors in laboratory experiments cannot be excluded, the established model and microscopic parameters are considered reasonable and can reflect the microscopic mechanical response process of CRC composite materials with different rubber particles content, although there are some minor    Advances in Materials Science and Engineering differences between the numerical simulation and the laboratory test in the peak phase and after peak phase of some materials.

Energy Evolution Rules and Microscopic Damage Mechanism
e failure process of materials is essentially the process of energy input, transformation, and dissipation.e failure process of CRC in uniaxial compression test is also an evolution process driven by energy as described below.Firstly, the loading plates of the testing machine input boundary energy to the concrete.As the loading progresses, the input energy can be gradually transformed into the releasable elastic strain energy, the sliding friction energy needed for the microcrack propagation and the fracture surface formation, and other dissipative energy such as particle kinetic energy and damping dissipation energy during the loading of samples.

e Microcracks Propagation in Damage Process of CRC.
In the numerical model, the spherical particles are used to simulate aggregates in the CRC and the parallel bonds are used to simulate the cement and agglomerant between particles.When the shear stress between the particles is higher than its bond strength, the bond breaks and the particles slide, which can be considered as one microcrack generating inside the concrete.
e process and rules of microcracks propagation of the CRC are simulated and analyzed by tracing the formation and number of broken bonds between particles.
Figure 11 shows the curve of the initiation and evolution of microcracks in CRC with different rubbers content.It can be seen from the figure that the microcrack number grows as "S" shape with the increase of strain during the loading process.e generating rate of cracks is relatively smooth in the beginning, then increases rapidly, and gradually becomes slow and stable after reaching the peak intensity.e time of cracks generation gradually lags behind with the increase of rubber contents, and the speed of cracks generation gradually decreases, resulting in a corresponding reduction in the total number of generated microcracks.
It can be seen from Figure 12 that the axial deformation at the initial microcrack formation of different CRC samples presents a "concave-upward" relationship with the of rubber content.e strain when the initial crack generation in the NC sample (without rubber) is 0.424%, and the corresponding strains of five groups of samples with 5%, 10%, 15%, 20%, and 25% rubber content are basically stable at about 0.2%, while the corresponding strains of the two groups of samples with 30% and 40% rubber content are 0.28% and 0.48%, respectively, showing an upward trend.
e reasons for that are as follows: Since the higher cohesive force among sand particles in NC, the first microcrack generates when the stress between particles exceeds the ultimate strength of bond, resulting in the relatively large deformation of the NC. e addition of rubber reduces the bonding strength of the particles between sands and rubbers, which leads to the microcrack occurs when the strain is relatively small (e.g., five groups of samples with 5%, 10%, 15%, 20% and 25% rubber content).With the increase of rubber content, the deformation of rubber particles begins to be the main factor of high strain of samples, and the initial crack occurs when the samples have a relatively large deformation (e.g., two groups of samples with 30%, 40% rubber content).Advances in Materials Science and Engineering e overall relationship between the number of microcracks generated and the content of rubber particles presents the "concave-downward" shape when the axial strain of samples is 2%.Accordingly, when the strain of materials is the same 2%, the number of cracks in the concrete without rubber is 7250, and the number of cracks in samples with 10% rubber and 20% rubber are, respectively, 8550 and 8250.en the amount of cracks of CRC begins decrease with the increase of rubber contents, and the group with the largest rubber content has the least cracks which are 3230.e results show that the deformation and flexibility of rubber particles can effectively reduce the cracking of normal concrete and improve the overall stability of the material.

Rules of Energy Dissipation and Damage Mechanism.
By using particle discrete element method to track and record the mechanical body energies, mechanical contact energies, and the formation of microcracks in the process of material destruction, the correlation between the energy dissipation and the damage mechanism of materials is analyzed.In the particle discrete element theory, the input energy U can be expressed as where U pre is the total input energy at the end of the last time step, and F 1 , F 2 are the load of the up and down loading plates of the time step, respectively, and ΔU 1 , ΔU 2 are the displacements of the loading plates, respectively, corresponding to the load at the time step.As for the PBM, the elastic strain energy U e stored in the particle aggregate consists of strain energy E k of contacting particles and the  The number of microcracks (1) Axial strain (%) 5% rubber 10% rubber 15% rubber 20% rubber 25% rubber 30% rubber 40% rubber  8 Advances in Materials Science and Engineering strain energy E k of parallel bond, and their formulas are as follows respectively: where F l n is the normal force of contacting particles, F l s is the tangential force of contacting particles, k n is the normal stiffness of contacting particles, k s is the shear stiffness of contacting particles, F n is the normal force of bond, F s is the tangential force of bond, M t is the torsional moment, M b is the bending moment, A is the cross-sectional area of bond, I is the moment of inertia of cross section of bond, and J is the polar moment of inertia of cross section of bond.
According to the second law of thermodynamics, the testing system is assumed to be a closed system without any heat exchange, the total input energy U consists of the elastic strain energy U e stored in the system and the dissipated energy U d which is necessary for the sample's internal damage and plastic deformation, and the formula is as follows: ere are many forms to describe the macroscopic damage variable of rock and soil, and usually, it can be described as crack density, energy accumulation, fatigue life, mechanical parameter, and so on, and the method of degradation of dynamic elastic modulus [29] based on the effective stress and the strain equivalence principle is widely used, and the formula of the damage variable D 0 is where E 0 is the elastic modulus of materials at undamaged state, and E is the elastic modulus at the damage state.It can be seen from the expression that when E changes to E 0 , D 0 is 0, and it is the original and undamaged state.us, the lager D 0 is, the more severe the damage of material is.
In the analysis of the particle discrete element model, based on the microstructure and the fracture mechanism of bond, the ratio of the number of fractured bonds N to the total number of bonds between particles N 0 is regarded as the mesoscopic damage variable D m [30], and the expression is e damage variables and dissipation energy of eight groups of samples at different stages of the loading process are recorded and analyzed, as shown in Table 4, from which the relationship between them can be obtained as follows: When the CRC samples reach the peak strength, the variation trend of dissipation energy ratio is basically consistent with that of damage variable with the increase of rubber particles content, showing a trend of increase first and then decrease.e ratio of dissipation energy to total energy in the CRC group with 25% rubber content reaches the maximum of 41.4% and the corresponding damage variable is also the largest at 0.3, indicating that the overall structural stability of the sample is relatively the best and it needs to dissipate enough energy to reach its damage threshold.And other samples with a rubber content of 0%, 5%, or 10% have the relatively high integral strength and low resist deformation ability due to the less rubber content, and the brittle failure occurs when the damage variable is low.Similarly, samples with a rubber content of 30% and 40% also suffer destruction when the damage variable is relatively small due to the low bonding strength between the particles caused by the incorporation of more rubber particles, and the corresponding dissipation energy ratio is also relatively low.
e statistical analysis of damage variables and dissipation energy is carried out in the postpeak stage of material destruction, which shows that the damage variable and the dissipative energy ratio after the destruction of materials increase with the increase of rubber content.e proportion of the dissipated energy increases from 81.7% to 98.4% as the rubber content increases, showing that most of the total input energy of CRC is dissipated after the material reaches destructive stability and the rest energy is stored in the retained samples in the form of elastic strain energy among particles.
e microdamage variable at the peak strength is used as the material's damage threshold D cri , which can be seen at formula (7), and it is defined as the ratio between the number of broken bonds N cri at the peak strength and the total number of bonds N 0 .e damage threshold represents the circumstance of internal damage when the material reaches the ultimate strength, that is, the critical value of structural damage in the material under external load.And when the damage degree of material exceeds the value, the structure is completely destroyed and loses its bearing capacity.
Figure 13 shows the relationship between the damage threshold and the number of microcracks at the peak strength of CRC samples with different rubber contents.It can be seen from the figure that with the increase of rubber content, the total number of bonds in samples shows a decreasing trend with the numbers being 51510, 39896, 38364, 37480, 35034, 31712, 29564, and 22847, respectively, and the corresponding number of microcracks at the peak presents a "increase first and then decrease" trend with the numbers being 1510, 2950, 4215, 4420, 6310, 9610, 7490, and 1520, respectively.e damage threshold also presents the "increase first and then decrease" trend with the addition of the rubber contents, and the damage thresholds are 0.029, 0.074, 0.11, 0.118, 0.18, 0.3, 0.253, and 0.067, respectively.Wherein the group with rubber content of 25% has the largest damage threshold which is 0.3, showing that this sample reaches its ultimate strength when the number of cracks accounts for 30% of the number of total bonds and has the maximum overall utilization; that is, only when the overall damage degree of this sample reaches 30%, can it lose the bearing capacity.

Advances in Materials Science and Engineering
As for the group without rubber contents, it reaches the ultimate strength when the damage thresholds value is 0.029, showing that the material in this group has reached the ultimate bearing capacity when the overall damage degree comes to 2.9%.Similarly, the two groups with rubber content of 5% and 40% have been completely destroyed when their overall damage degrees are less than 10%.
e elastic strain energy release ratio E rt is defined from the perspective of energy to describe the energy release state of the material after the peak.E rt is defined as the ratio that the absolute value of difference value between the maximum elastic strain energy stored in the materials E p and the elastic strain energy E s at the time of the stability of destroy compared with the absolute value of the corresponding strain difference value, which represents the release value of the elastic strain energy in the unit strain after the sample reaches the peak strength.e bigger the E rt is, the bigger the energy release value in the unit strain is after material the reaches the ultimate bearing capacity, and the more brittle the material is, the more severe the degree of damage is.Oppositely, the smaller the value of E rt is, the more pliable and tough the material is and the more smooth and steady the degree of damage is after the material reaches the peak strength.
Figure 14 shows the relationship between the different rubber contents in CRC and the corresponding elastic strain energy release ratio.It can be seen that the elastic strain energy release ratio comes to 48.8 J/% after the concrete reaches the peak strength without adding rubbers.When rubber particles are added to the concrete, the energy release ratio decreases sharply, and presents negative exponential function type with the downtrend of first sharp and then slow, and other groups' energy release ratios are 8.83 J/ %, 7.54 J/%, 7.36 J/%, 6.17 J/%, 5.47 J/%, 3.29 J/%, and 1.15 J/ %, respectively, which can quantitatively describe the relation between the energy release rate and the rubber particle contents by using the fitting function y � 37.53 exp(−5.622x)+ 11.27 exp(−0.036x) in the figure .And the fitting effect is good with the fitting variance value of 0.9971.e analysis of the fitting curve shows that the presence or absence of the rubber particles has more sensitive effect to the energy release ratio than the content of the rubber particles, indicating that rubber materials can significantly reduce the brittleness of concrete and improve the overall flexibility and ductility of materials.

4.3.
e Correlation of Material Damage and Energy Evolution. Figure 15 shows the correlation rule between damage process and energy evolution in the group without rubber contents (0% of the rubber particle contents, (a)) and the group with maximum damage threshold (25% of the rubber particle contents, (b)).e damage mechanism and the energy evolution rule of the CRC composite materials and normal concrete materials are analyzed by taking the group with maximum damage threshold and the group without rubber contents as representatives.
e analysis of Figure 15 can be obtained: In the initial loading phase (OA), the boundary input energy of the two materials both shows a slowly increasing trend, and basically coincide with the elastic strain energy.
e sliding friction energy, kinetic energy, and other dissipated energy are basically zero and can be neglected.e boundary input energy of materials is completely transformed into the elastic strain energy stored in the materials, and no microcracks are generated in this phase.The number of microcracks (10 4   ) Thedamage threshold (1)

Mass fraction (%)
The damage threshold The total microcracks of material The microcracks when reaching the peak intensity  Advances in Materials Science and Engineering In the linear elastic phase (AB), the input energy and elastic strain energy begin to separate, and both of them show a slow development trend in the shape of "upward concave."e elastic strain energy accumulates and stores gradually, the microcracks in the materials begin to occur, and the sliding friction energy and other dissipated energy begin to increase slowly.e dissipation energy curve basically coincides with the frictional energy curve, showing the dissipated energy is mainly used for sliding friction consumption of materials at this phase.Macroscopically, the stress-strain curve shows the linear elasticity, which is the stable development stage of the materials.
In the yield phase (BC) of materials, the curves of the input energy and the elastic strain energy are completely separated; both of them show a rapid linear growth in general, and the growth rate of input energy is larger than that of the elastic strain energy.e curves of sliding friction energy and dissipated energy begin to separate and grow rapidly in a shape of "upward concave," and a large number of microcracks are rapidly generated inside the material.It indicates that a large amount of damage occurs in the microstructure of the material at this phase, and the particles are mutually slipping and moving, and the dissipated energy represented by the sliding friction energy between particles increases significantly, which shows the materials reach the ultimate bearing capacity at the point C.
In the postpeak phase (after the point C), the curve of material constitutive relation lags behind that of energy Advances in Materials Science and Engineering evolution as energy is the driving factor for damage and fracture of materials.e boundary input energy continues to increase in a linear trend and tends to be gentle in the later phase, while the elastic strain energy continues to grow in a shape of "upward convex," and begins to decline after reaching the peak at the point D. e sliding friction energy, kinetic energy, and other dissipated energy increase rapidly and gradually tend to be stable in the later phase, and the generation rate of microcracks gradually reduces to zero.

Discussion
e stress-strain relationship of CRC materials is basically consistent with the development trend of their elastic strain energy, but the stress-strain curve lags behind the development of elastic strain energy, which indicates that energy is an important driving factor for the damage and fracture of materials.e CRC with 25% of the rubber mass has a larger yield phase and smaller elastic strain energy release ratio than the NC, which can be explained as the addition of rubber particles makes more external input energy to be converted into dissipative energy required for microcracks propagation and sliding friction between particles and released step by step.
e purpose of CRC materials research is to improve the deformation capacity and toughness as far as possible on the basis of maintaining a certain compressive strength to prevent sudden failure of the project.
e size of rubber particles used in this study is 3-5 mm, and when the rubber mass fraction is 25%, the coordination between the strength and deformation of the material was the best.In the later research, experiments can be carried out on different CRC rubber particle sizes and different rubber treatment methods to obtain a more ideal CRC ratio scheme with better mechanical and deformation capacity.
e acoustic emission experiment can be used to study the micromechanical failure mechanism of CRC materials to verify the effectiveness and accuracy of conclusions obtained by this method.
Different pretreatment methods of rubber particles before CRC production may affect the overall compatibility, strength, and other properties of the material.We will try a variety of different treatment methods of rubber particles in the later research, including untreated rubber particles, acid treatment, and mixed treatment, etc., to study the effect of functionalization of rubber particle surface on the compatibility of CRC material components.us, a better treatment method of rubber particles may be obtained to improve the overall strength and deformation performance of CRC.

Conclusions
In this paper, the energy evolution law and mesodamage mechanism of CRC (crumb rubber concrete) are investigated both by laboratory experiments and the particle discrete element method, and the energy evolution behavior and damage characteristics of NC (normal concrete) are compared and analyzed.e results of this investigation are summarized as follows: (1) With the increase of rubber particles content, the axial deformation of CRC gradually increases, the integrity of the samples is better after reaching its ultimate bearing capacity, and the failure form is gradually changed from a form of single sloping and complete plane to a form of short and dispersed cracks.e proposed fitting functions can well fit the relationship among the compressive strength, elastic modulus, and the rubber particles content, respectively.(2) e mesoscopic physical and mechanical parameters of the CRC with rubber particles and sand particles as the main aggregate are analyzed and determined, which can better simulate the macroscopic mechanics and deformation characteristics of CRC materials.(3) It is considered that the mesoscopic damage threshold of CRC materials with different rubber contents has a strong correlation with the proportion of dissipated energy at the peak strength, and both of them show the tendency of increasing first then decreasing with the increase of rubber mass.(4) e damage threshold of the group with 25% of the rubber mass is the largest, and it is 30%, and the corresponding dissipated energy ratio at the peak strength is 41.4%, which indicates that the composite materials need to dissipate 41.4% of the input energy, and when its damage state reaches 30% of the overall integrity, it will reach the ultimate bearing capacity.(5) e elastic strain energy release ratio is proposed to describe the elastic strain energy released by the unit strain after the material reaches the ultimate bearing capacity.e relationship between elastic strain energy release ratio of CRC and rubber particle contents is obtained and fitted by the negative exponential function y � 37.53 exp(−5.622x)+ 11.27 exp(−0.036x).
It is considered that the rubber composition can significantly reduce the energy release ratio and obviously improve the flexibility and ductility of the composite concrete materials.

Figure 1 :
Figure 1: Rubber particles before and after surface-pretreated with NaOH saturated solution.

Figure 3 :Figure 4 :
Figure 3: e shape of the destruction (from left to right: O, A, B, C, D, E, F, H).
Dashpot force (F d ), not shown Linear force (F 1 ), linear elastic (no tension) and frictional Bone load (F and M)

Figure 5 :Figure 6 :Figure 7 :
Figure 5: Mesoscopic contact model and failure criterion of the PBM.

Figure 10 :
Figure 10: Comparison of constitutive curves obtained by laboratory experiments (a) and those obtained by numerical simulations (b).

Figure 11 :
Figure 11: e initiation and evolution of microcracks.

Figure 13 :
Figure 13: Relationship between damage threshold, number of microcracks, and content of rubber particles.

Figure 15 :
Figure 14: e relationship between energy release ratio of CRC and rubber content.

Table 1 :
Mix proportion of rubber concrete in experiment (unit: kg/m 3 ).

Table 2 :
Mechanical parameters and PFC results.

Table 3 :
Mesomechanical parameters of the CRC.

Table 4 :
Statistics for the damage variable and dissipation energy.