Size Effect of Shear Strength of Recycled Concrete Beam without Web Reinforcement: Testing and Explicit Finite Element Simulation

Recycled concrete is a form of low-carbon concrete with great importance. The explicit finite element method is an economical and feasible method for analyzing static concrete structures, such as those made of recycled concrete. The shear strength of regular concrete beams has size effects. In this study, a group of physical tests on the size effect of the shear strength of recycled concrete beams without web reinforcement was carried out under the condition of a constant shear span ratio. The research results show that the shear strength of the test beam generally decreases with the increase in beam section height, and a regression formula of the shear strength was obtained, which can formulate this effect. The rationale and feasibility of the explicit finite element method solving the ultimate load of concrete structures (which can derive the shear strength) were briefly demonstrated, and an explicit finite element simulation of test beams was carried out. Results showed an obvious and phenomenologically regular size effect of the shear strength of recycled concrete beams without web reinforcement, which can be simulated by the explicit finite element method. This research aims to promote the study of low-carbon recycled concrete structures to a certain extent and encourage the application of economic explicit finite element methods for the static analysis of concrete structures.


Introduction
Recycled concrete is an important form of low-carbon-emission concrete; similarly, numerical simulation methods are generally regarded as low-carbon analysis methods compared with physical tests. The mass formed by waste concrete processing is commonly utilized as coarse aggregate, and the material formed by mixing it with inorganic binding substances, such as fine aggregate, cement, water, etc., is called recycled concrete (ReC) [1]. It is generally considered that ReC is a green, low-carbon-emission material, which has a certain positive role in the protection of resources and the environment [2]. It is worth noting that, from a broad perspective of processes in research, production and living, among multiple behaviors that can achieve the same goal, those producing lower greenhouse gas emissions are collectively called low-carbon behaviors [3]. Physical tests and numerical simulation tests are two common methods of studying building structures, and the finite element method is a regular means to carry out structural numerical simulation tests [4]. In most cases, the numerical simulation method is obviously a relatively low-carbon analysis method.
in the incremental form [40]. The implicit algorithm is often recommended as a general method to solve nonlinear problems in commercial finite element software [41]. One of its limitations is that, when solving structural static problems with local instability phenomena (such as cracking when concrete structures are forced), it is generally difficult to obtain a convergent solution [42]. For concrete beams, when the external load has approached the peak load of SRRC [23], inclined cracks have already appeared in the shear area of the concrete beam. From the foregoing of two points, it is easy to infer that it is difficult for the implicit algorithm to research the SRRC size effect. The essence of the explicit algorithm is to convert the static problem of the original structure into the corresponding structural dynamic problem, perform a pseudo-static analysis on this problem, and approximate the result of the pseudo-static analysis as the solution to the original static problem [43]. The static problem that the explicit algorithm can approximately solve does not depend on the material, geometry, and continuity characteristics of the problem at all [44]. To date, there have been many reports on the use of explicit algorithms to carry out approximate static analysis of concrete structures [45][46][47]. For example, Yao et al. [46] used an explicit finite element method to analyze the four-point bending beam of concrete in order to verify the reliability of the explicit algorithm in solving the quasi-static response of concrete beams. There is good agreement with the test results. Yu et al. [47] used the arch effect to explain the shear mechanism of variable cross-section beams, conducted finite element analysis on four cantilever beams without web reinforcement, proposed the influence coefficient of the compression inclination on the arch effect and the method to determine the position of the check section, and established the relationship between the standard formula for shear resistance and the formula for calculating the shear resistance of beams with variable cross-sections. However, no study exists on the use of explicit algorithms to investigate the size effect of SRRC. Therefore, the present study aims to use the explicit finite element method to carry out simulation experiments on the SRRC size effect.
In summary, the main contribution of this work is the testing and explicit finite element simulation of the shear strength of a group of recycled concrete beams without web reinforcement under the condition of a constant shear span ratio.
The purpose of this paper is to report that SRRC (the concrete itemized shear resistance of recycled concrete beams) has a size effect, and this effect can be simulated by the explicit finite element method but is difficult to simulate using the implicit finite element method.
The paper is structured as follows. An overview of the physical tests is provided in Section 2, including the specimen design, the loading equipment, and the loading system. In Section 3, an overview of the simulation tests is given, including the physical discrete and contact settings, material constitutive model selection, loading system, and solution algorithm. In Section 4, some typical test results and related discussions are provided first. Next, a regression formula is presented, which can reflect the size effect of shear strength. Section 5 concludes this paper.

Specimen Design
A total of four recycled concrete beams numbered B120, B180, B240, and B300 were designed and poured with the same batch of recycled concrete. Before pouring the 4 test specimens, the test materials and test mix ratios need to be determined. They are introduced as follows [22].
Firstly, the test materials include natural fine aggregate, recycled coarse aggregate and cement.
In this study, the performance indicators for natural fine aggregate are presented in Table 1. When determining the recycled coarse aggregate, certain factors are also mainly taken into account, such as density, water absorption, moisture and size. For details of the performance indicators of the coarse recycled aggregate selected in this study, see Table 2. When determining the cement, some factors are mainly considered compressive, such as density. The cement used in this test is ordinary Portland cement with a strength level of 32.5. For some performance indicators of the cement selected in this study, see Table 3. Secondly, the mix ratio is designed. The ratio used to produce ReC, the strength level of which is C30, the slump is 70 mm, and which was used to cast the 4 beams tested in this study, is given by Table 4. After testing, the average uniaxial compressive strength of this recycled concrete was and the initial elastic modulus was Based on the purpose of the test, the test beam was not equipped with stirrups. However, so that insufficient bending strength would not affect the results of the shear strength test, two longitudinal bars were arranged at the bottom of the beam. The position, yield strength ( f y ), elastic modulus (E s ) and area (A s ) information of the longitudinal bars in each beam, and the size information of each beam, are shown in Figure 1 and Table 5. Only two supporting points of the test beam were located at the bottom of the beam, and the form was simply supported; the only concentrated force loading point was at the top of the beam. The specific positions of the supporting point and the loading point are shown in Figure 1 and Table 5. It is easily deduced that the beam's shear span ratio is λ = l 2 /h 0 = 2.45, where l 2 is the horizontal distance between the supporting point and the loading point, and h 0 is the effective height. After the estimation and according to the above-mentioned geometric dimensions, steel bar configuration, material strength and loading design, the failure mode of recycled concrete beams will be guaranteed to be shear failure, and thus the research purpose can be achieved. /ℎ = 2.45, where is the horizontal distance between the supporting point and the loading point, and ℎ is the effective height. After the estimation and according to the above-mentioned geometric dimensions, steel bar configuration, material strength and loading design, the failure mode of recycled concrete beams will be guaranteed to be shear failure, and thus the research purpose can be achieved.

Loading Equipment and Loading System
The loading details of the test are depicted in Figure 2, where Figure 2a shows the reaction frame, and Figure 2b includes the details of the loading configuration. During loading, the preset loading schedule was implemented by simultaneously 'controlling the oil output rate of the oil pump to slowly move the jack lifting sleeve down, and monitoring the force value of the force indicator connected to the force sensor'.

Loading Equipment and Loading System
The loading details of the test are depicted in Figure 2, where Figure 2a shows the reaction frame, and Figure 2b includes the details of the loading configuration. During loading, the preset loading schedule was implemented by simultaneously 'controlling the oil output rate of the oil pump to slowly move the jack lifting sleeve down, and monitoring the force value of the force indicator connected to the force sensor'. /ℎ = 2.45, where is the horizontal distance between the supporting point and the loading point, and ℎ is the effective height. After the estimation and according to the above-mentioned geometric dimensions, steel bar configuration, material strength and loading design, the failure mode of recycled concrete beams will be guaranteed to be shear failure, and thus the research purpose can be achieved.

Loading Equipment and Loading System
The loading details of the test are depicted in Figure 2, where Figure 2a shows the reaction frame, and Figure 2b includes the details of the loading configuration. During loading, the preset loading schedule was implemented by simultaneously 'controlling the oil output rate of the oil pump to slowly move the jack lifting sleeve down, and monitoring the force value of the force indicator connected to the force sensor'.    The whole loading history was divided into two stages (pre-loading and formal loading) with different loading schedules. Before formal loading, preload was applied at about 50% of the estimated cracking load, F cr,c = M cr,c /l 2 , where the cracking moment for the flexural load capacity of a normal section of beams without web reinforcement, M cr,c , was calculated according to the formula given in the literature [24]: In the formula, σ tp is the peak tensile stress, and based on the formula derived from literature [24] and literature [48], After calculating, σ tp equals to 2.46 MPa. The formal loading phase was divided into three subphases. Firstly, the increment of ∆F 1 ≈ 0.1F cr,c was used to increase the concentrated load F from 0 to the actual cracking load, F cr,c . Secondly, ∆F 2 ≈ 0.1 F q,c − F cr,t was applied to increase F = F cr,t approximately to F q,c . The latter, F q,c , is the quantile value of the estimated ultimate load (F u,c ) and in this study is equal to 0.7F u,c . This paper used the formula V u,c = 1.75σ tp bh 0 /(λ + 1) to estimate F u,c = 2V u,c , given by the literature [25] for the shear bearing capacity of the oblique section of beams without web reinforcement. Thirdly, ∆F 3 ≈ 0.1 F u,c − F q,c was used to increase F = F q,c to the measured limit load, F u,t . M cr,c , F cr,c , F u,c , F q,c , ∆F 1 , ∆F 2 , and ∆F 3 were calculated according to the above principles and F cr,t , as seen in Table 6, thus determining the load control points during the entire test loading history.

Physical Discrete Setting and Contact Conditions
In this paper, Abaqus software was used to carry out the geometrical discrete and contact settings of the finite element model of the specimen. The operation mainly involved six aspects. Firstly, the element type was determined. Both the beam entities and the steel cushion blocks were discrete, with an eight-node reduced integral format threedimensional entity element. The number for this element in the Abaqus element library is C3D8R. Reinforcing bars were discretized by three-dimensional two-node line elements, with the number T3D2. Secondly, the element size was determined. For the rectangular solid elements of discrete steel cushion blocks and beam entities, the three orthogonal edges on the surface were parallel to the three orthogonal edges on the beam surface (see Figure 3), and the dimension of x l along the beam length was l e , the dimension of the direction x h along the beam height was h e , and the dimension of x b along the beam width was b e . For the concrete element in the column of the cushion blocks, however, the dimension in the x l direction was l e,in ; for the concrete element from the beam end to the supporting cushion block area, the dimension in the x l direction was l e,out . The dimensions of all the above physical elements are shown in Table 7. For the linear element of the discrete steel bars, the dimension l es along the beam length direction x l is also shown in Table 7, and its position in the plane x b -x l is shown in Figure 1. Thirdly, it was determined how steel and concrete interact. This article assumed that there was no bond-slip behavior between the steel bar and the concrete, and thus the 'Embedded' constraint was selected to embed the steel bar in the concrete entity. Fourthly, the contact mode between the steel cushion block and the concrete entity was determined. Herein, the mechanical behaviors of the upper and lower neighborhoods of the contact surface between the cushion blocks and the beam are set as follows: the contact surface between the left and right supporting cushion block and the bottom of the beam, and the contact surface between the loading cushion block and the top of the beam, adopted the 'Tie' mode to combine into a continuous medium area. Fifthly, boundary conditions were identified. In order to better simulate actual working conditions, the lower surface of the left and right supporting cushion block only had rotational degrees of freedom around x b and x h . Sixthly, the constitutive response law of the neighboring concrete at the concentrated force was assumed. Extremely high stresses occur in concrete elements near the supporting/loading cushion blocks. In order to prevent these from exceeding the description range of the concrete damage constitutive model and subsequently causing calculation failures, the constitutive model of the steel cushion block was adopted for the concrete elements that are in contact with the left and right supporting cushion blocks and for the concrete element in contact with the loading cushion block (as shown in the red area in Figure 3). The mathematical expression of this model is shown in Section 3.2. According to the above six operational steps, a threedimensional view of the finite element model of the beam, numbered B-120, is shown in Figure 3.
behaviors of the upper and lower neighborhoods of the contact surfac ion blocks and the beam are set as follows: the contact surface betwe supporting cushion block and the bottom of the beam, and the cont the loading cushion block and the top of the beam, adopted the 'Tie into a continuous medium area. Fifthly, boundary conditions were id better simulate actual working conditions, the lower surface of the lef ing cushion block only had rotational degrees of freedom around constitutive response law of the neighboring concrete at the concen sumed. Extremely high stresses occur in concrete elements near the cushion blocks. In order to prevent these from exceeding the descripti crete damage constitutive model and subsequently causing calculatio stitutive model of the steel cushion block was adopted for the concre in contact with the left and right supporting cushion blocks and for t in contact with the loading cushion block (as shown in the red area in F ematical expression of this model is shown in Section 3.2. According erational steps, a three-dimensional view of the finite element mode bered B-120, is shown in Figure 3. is the concr the constitutive model of the steel cushion block).

Material Constitutive Model
Before static analysis, constitutive models of various componen ment discrete model need to be selected.
Firstly, the constitutive model of the steel cushion blocks was sel the local influence principle [49], the selection of the constitutive mod ion block has a limited effect on the simulation results of the static beh beam. Therefore, in static loading, it can be assumed that the steel pad elastic state. Therefore, this study selected a linear elastic constitutive cushion block, which can be expressed as where is the Cauchy stress, is the Cauchy strain tensor, and ness tensor of the material, which is determined by the elastic modul

Material Constitutive Model
Before static analysis, constitutive models of various components in the finite element discrete model need to be selected.
Firstly, the constitutive model of the steel cushion blocks was selected. According to the local influence principle [49], the selection of the constitutive model of the steel cushion block has a limited effect on the simulation results of the static behavior of the concrete beam. Therefore, in static loading, it can be assumed that the steel pad is always in a linear elastic state. Therefore, this study selected a linear elastic constitutive model for the steel cushion block, which can be expressed as where σ is the Cauchy stress, ε is the Cauchy strain tensor, and D e is the elastic stiffness tensor of the material, which is determined by the elastic modulus E b and the Poisson's ratio ν b of the steel cushion block [49]. In the simulation of this study, E b = 190 GPa, and ν b = 0.3.
Secondly, the constitutive model of steel bars was selected. The simulated beam was a simply supported beam without web reinforcement. When it failed due to the insufficient shear strength of the oblique section, the longitudinal reinforcement had not yet generally entered the yield stage. Therefore, this study employed an ideal (onedimensional) elastoplastic constitutive model for steel bars, which can be expressed as where σ s and ε s indicate stress scalar and strain scalar, respectively; and ε sy = f y /E s is the initial yield strain, where f y stands for the yield strength. In the simulation, the E s and f y of the longitudinal reinforcement in each simulated beam adopt the actual measured values of the longitudinal reinforcement in the corresponding test beam; see Table 5 for details. Thirdly, the recycled concrete constitutive model was chosen. 1 In the static test, there is a cracking phenomenon in the concrete, strain softening in the cracked area, and elastic-plastic strain unloading in the adjacent cracked area. Accordingly, this paper used the concrete damage plasticity (hereinafter referred to as CDP) constitutive model proposed by J. Lee et al. [50] for recycled concrete, which can describe the above behaviors and can be expressed as where d is the damage variable (scalar) and ε p is the plastic strain tensor. According to this constitutive model theory [51,52], in the analysis, the stress-strain curve under the uniaxial (tension or compression) condition, the equivalent compression plastic strain ε p c , and the equivalent tension plastic strain ε p t should be given. With the aid of several conversion rules, Formula (4) can be utilized to calculate the constitutive state of the integration point in the element. 2 The literature [53] has pointed out that the stress-strain curve shapes of recycled concrete and original concrete are highly similar. Therefore, this study applied the uniaxial compressive stress-strain curve of concrete, as proposed by Guo [54], and the uniaxial tensile stress-strain curve to construct the basic data required for numerical implementation of the above damage constitutive model. In the above two formulas, σ c and ε c are the compressive stress and compressive strain, respectively; σ cp is the peak compressive stress; ε cp is the strain corresponding to σ cp ; σ t and ε t are the tensile stress and tensile strain, respectively; σ tp is the peak tensile stress; ε tp is the strain corresponding to σ tp ; and α c , β c , α t , β t and ϕ t are all coefficients calibrated by the tests. Based on the tested σ cp (the value is shown in Formula (1)), first calculate σ tp according to Formula (2), and then calculate ε tp according to the formula ε tp = 65σ tp 0.54 × 10 −6 [54]. The values of these three data are listed in Table 8. At the same time, according to [53,54], the selected or calculated values of ε cp , α c , β c , α t , β t and ϕ t are also listed in Table 8. According to Formulas (5) and (6) and the parameters in Table 8, the stress-strain curve is presented in Figure 4. 3 As suggested by the literature [52], the formula is chosen to describe ε p c and ε p t . In the formula, E c is the initial elastic modulus of concrete (its value is shown in Formula (2)), and γ c and γ t are the coefficients calibrated by the experiment, empirically taking the values of γ c = 0.7 and γ t = 0.1, based on the literature [52]. ing to the formula = 65 . × 10 [54]. The values of these three data are listed Table 8. At the same time, according to [53,54], the selected or calculated values of , , and are also listed in Table 8. According to Formulas (5) and (6) and parameters in Table 8, the stress-strain curve is presented in Figure 4. ③ As suggested the literature [52], the formula is chosen to describe ̃ and ̃ . In the formula, is the initial elastic modulus of c crete (its value is shown in Formula (2)), and and are the coefficients calibrated the experiment, empirically taking the values of = 0.7 and = 0.1, based on the li ature [52]. In order to increase the readability of the paper, follow the data input requireme of the CDP model in Abaqus. Four curves on the relationships of yield-stress (under co pression) to inelastic-strain, compression-damage to inelastic-strain, yield-stress (un tension) to cracking-strain, and tension-damage to cracking-strain are shown in Figur   In order to increase the readability of the paper, follow the data input requirements of the CDP model in Abaqus. Four curves on the relationships of yield-stress (under compression) to inelastic-strain, compression-damage to inelastic-strain, yield-stress (under tension) to cracking-strain, and tension-damage to cracking-strain are shown in Figure 5. is yield-stress under compression, is inelastic-strain, is compression-damage, is yieldstress under tension, is cracking-strain, and is tension-damage).

Loading System and Solution Algorithm Configuration
Before performing the finite element simulation, it is also necessary to define the loading schedule and solving algorithm. In a similar way to the physical test, the displacement control loading mode is adopted in the simulation test, with a constant loading rate. The total loading time of each numerical beam model and the initial displacement/final displacement of the active loading point are shown in Table 9.
Configure the parameters of the implicit/explicit solver. For the implicit solver, the time period was set to 1, the loading point displacement changed linearly with time, the "Nlgeom" option was checked, the initial increment size was set to 0.001, the minimum increment size was set to 1 × 10 , the maximum increment size was set to 0.01, and the iterative solution of the nonlinear equation system was Newton's method. For the explicit solver, the time period and the characteristics of loading point displacement in the time domain are shown in Table 9; check the "Automatic" option to determine the loading point increment, check the "Element-by-element" option to determine the stable increment step size, the time scaling factor was set to 1, the linear bulk viscosity parameter was set to 0.06, the quadratic bulk viscosity parameter was set to 1.2, and other parameters used the default values.

Loading System and Solution Algorithm Configuration
Before performing the finite element simulation, it is also necessary to define the loading schedule and solving algorithm. In a similar way to the physical test, the displacement control loading mode is adopted in the simulation test, with a constant loading rate. The total loading time of each numerical beam model and the initial displacement/final displacement of the active loading point are shown in Table 9. Configure the parameters of the implicit/explicit solver. For the implicit solver, the time period was set to 1, the loading point displacement changed linearly with time, the "Nlgeom" option was checked, the initial increment size was set to 0.001, the minimum increment size was set to 1 × 10 −5 , the maximum increment size was set to 0.01, and the iterative solution of the nonlinear equation system was Newton's method. For the explicit solver, the time period and the characteristics of loading point displacement in the time domain are shown in Table 9; check the "Automatic" option to determine the loading point increment, check the "Element-by-element" option to determine the stable increment step size, the time scaling factor was set to 1, the linear bulk viscosity parameter was set to 0.06, the quadratic bulk viscosity parameter was set to 1.2, and other parameters used the default values.

Failure Phenomenon and Characteristic Stress Contour Plots
The progress of crack development in the four beams was consistent with common sense, and the damage mode was typically uniform. When the load reached a certain level, normal transverse bending cracks first appeared in the area of maximum bending moment (mid-span). As the load increased, the area with normal bending cracks appearing enlarged, and the already existing normal bending cracks elongated. When the load continued to rise, shear oblique cracks appeared near the neutral axis of the beam, the oblique cracks extended to both ends, and part of the oblique cracks connected with the normal cracks. As the load reached an extreme value, a main diagonal crack suddenly appeared to tear the beam in a similar manner to the main control crack (the red curve) in beam B-240, as shown in Figure 6a, and the bearing capacity of the beam was lost. The progress diagram of the cracks (the curve in the figure) in beams B-120, B-180, and B-300 are shown in Figure 6b-d, respectively, in which the red curve refers to the main control crack. To determine the plane position of the cracks, an equidistant grid was drawn along the horizontal and vertical directions of the beam, and the grid number was marked with numbers and letters. The value near several positions of the crack represents the value of the applied load (unit: kN) corresponding to the crack developed to reach this point. The crack progress diagram of beam B-240 was consistent with the characteristics of the first three beams.

Failure Phenomenon and Characteristic Stress Contour Plots
The progress of crack development in the four beams was consistent with common sense, and the damage mode was typically uniform. When the load reached a certain level, normal transverse bending cracks first appeared in the area of maximum bending moment (mid-span). As the load increased, the area with normal bending cracks appearing enlarged, and the already existing normal bending cracks elongated. When the load continued to rise, shear oblique cracks appeared near the neutral axis of the beam, the oblique cracks extended to both ends, and part of the oblique cracks connected with the normal cracks. As the load reached an extreme value, a main diagonal crack suddenly appeared to tear the beam in a similar manner to the main control crack (the red curve) in beam B-240, as shown in Figure 6a, and the bearing capacity of the beam was lost. The progress diagram of the cracks (the curve in the figure) in beams B-120, B-180, and B-300 are shown in Figure 6b-d, respectively, in which the red curve refers to the main control crack. To determine the plane position of the cracks, an equidistant grid was drawn along the horizontal and vertical directions of the beam, and the grid number was marked with numbers and letters. The value near several positions of the crack represents the value of the applied load (unit: kN) corresponding to the crack developed to reach this point. The crack progress diagram of beam B-240 was consistent with the characteristics of the first three beams. The typical stress contour plots provided by the finite element simulation can verify the aforementioned experimental phenomenon. The law of the main tension stress contour plots given by the finite element simulation was similar for all four beams, while this paper only selects the beam B-240 for analysis. In basic terms, Figure 7a shows the main tensile stress contour plot of the beam after removing the cushion blocks and the element immediately adjacent to the cushion blocks. From this figure, it is clear that there is a high principal tensile stress in the high-stress area near the cushion blocks and the high-bending effect area at the bottom of the mid-span beam, which is qualitatively consistent with the common sense of elasticity and material mechanics; the principal tensile stress in the vicinity of the neutral axis of the beam is substantial. This phenomenon can also be qualitatively described by the theory of material mechanics-the shear stress here is large, but the normal stress is zero. Subsequently, the distribution of the aforementioned shear stress peak area can also be described by the contour map of Tresca stress , as shown in Figure 7b ( = ( − ) 2 ⁄ ), where and are the maximum and minimum principal stresses, respectively. In particular, in the bending neutral axis area with high shear stress but low normal stress, the area with a large value will produce inclined cracks at an angle of 45 to the neutral axis, and this simulation results confirms the test results given in Figure 6.
In summary, crack development in the test is close to the result obtained by the simulation, and both belong to the shear failure mode.
(a) Principal tensile stress contour plot of beam B-240 The typical stress contour plots provided by the finite element simulation can verify the aforementioned experimental phenomenon. The law of the main tension stress contour plots given by the finite element simulation was similar for all four beams, while this paper only selects the beam B-240 for analysis. In basic terms, Figure 7a shows the main tensile stress contour plot of the beam after removing the cushion blocks and the element immediately adjacent to the cushion blocks. From this figure, it is clear that there is a high principal tensile stress in the high-stress area near the cushion blocks and the highbending effect area at the bottom of the mid-span beam, which is qualitatively consistent with the common sense of elasticity and material mechanics; the principal tensile stress in the vicinity of the neutral axis of the beam is substantial. This phenomenon can also be qualitatively described by the theory of material mechanics-the shear stress here is large, but the normal stress is zero. Subsequently, the distribution of the aforementioned shear stress peak area can also be described by the contour map of Tresca stress τ max , as shown in Figure 7b (τ max = (σ 1 − σ 3 )/2), where σ 1 and σ 3 are the maximum and minimum principal stresses, respectively. In particular, in the bending neutral axis area with high shear stress but low normal stress, the area with a large τ max value will produce inclined cracks at an angle of 45 to the neutral axis, and this simulation results confirms the test results given in Figure 6.
In summary, crack development in the test is close to the result obtained by the simulation, and both belong to the shear failure mode. , where and are the maximum and minimum principal stresses, respectively. In particular, in the bending neutral axis area with high shear stress but low normal stress, the area with a large value will produce inclined cracks at an angle of 45 to the neutral axis, and this simulation results confirms the test results given in Figure 6.
In summary, crack development in the test is close to the result obtained by the simulation, and both belong to the shear failure mode.

Ultimate Load in Test and Simulation
The ultimate loads given by the simulation and experiment were very similar. Explicit and implicit algorithms were used to carry out the simulation analysis. The ultimate load , was calculated by the explicit algorithm, and the ultimate load , was calculated by the implicit algorithm, as shown in Table 10. In the table, the subscript represents the test value, the subscript indicates the simulated value, the superscript represents the simulated value given by the explicit algorithm, and the superscript denotes the simulated value given by the implicit algorithm. In combination with the ultimate load , given by the test, the relative error is , , / , ( e, i).
Obviously, for each beam, is not only small but also much smaller than . This shows that the constitutive model used in this paper and the approximate solution algorithm of the selected nonlinear equations-the explicit algorithm-are reasonable. However, from the four sample points alone, there is clearly no qualitative relationship between and the height of the beam ℎ (See Table 5).

Law of Shear Strength and Its Value
Shear strength can be easily calculated based on the test and simulation results provided in the previous section. In concrete structure engineering, the shear strength of concrete beams is generally determined by the formula / ℎ [55], where is the beam width and ℎ is the effective height of the cross-section of the beam-the two values for each specimen are shown in Table 5-and represents the ultimate shear bearing capacity of the cross-section. According to Figure 1, /2, where is the ultimate load given by test or simulation (see Table 10 for specific values). According to the above, the calculated shear strength values are shown in Table 10.

Ultimate Load in Test and Simulation
The ultimate loads given by the simulation and experiment were very similar. Explicit and implicit algorithms were used to carry out the simulation analysis. The ultimate load F e u,s was calculated by the explicit algorithm, and the ultimate load F i u,s was calculated by the implicit algorithm, as shown in Table 10. In the table, the subscript t represents the test value, the subscript s indicates the simulated value, the superscript e represents the simulated value given by the explicit algorithm, and the superscript i denotes the simulated value given by the implicit algorithm. In combination with the ultimate load F u,t given by the test, the relative error is ε i err = F i u,s − F u,t /F u,t (i = e, i). Obviously, for each beam, ε e err is not only small but also much smaller than ε i err . This shows that the constitutive model used in this paper and the approximate solution algorithm of the selected nonlinear equations-the explicit algorithm-are reasonable. However, from the four sample points alone, there is clearly no qualitative relationship between ε err and the height of the beam h (See Table 5).

Law of Shear Strength and Its Value
Shear strength can be easily calculated based on the test and simulation results provided in the previous section. In concrete structure engineering, the shear strength of concrete beams is generally determined by the formula τ = V u /bh 0 [55], where b is the beam width and h 0 is the effective height of the cross-section of the beam-the two values for each specimen are shown in Table 5-and V u represents the ultimate shear bearing capacity of the cross-section. According to Figure 1, V u = F u /2, where F u is the ultimate load given by test or simulation (see Table 10 for specific values). According to the above, the calculated shear strength τ values are shown in Table 10.
The values of shear strength τ t given by the test show an obvious size effect phenomenon. This phenomenon not only has regression statistical characteristics but can also be confirmed by the explicit simulation results, τ s,e . The τ t (dependent variable), with h 0 as the independent variable under the constant shear span ratio (λ = 2.45) given by the experiment, is shown in Figure 8. Based on the set of discrete data points (h 0 , τ t ), the regression formula for describing the size effect is  The continuous curve of ℎ and obtained by the regression formula is shown in Figure 8; it is clearly seen that the fitting effect is qualitatively better. Either from the discrete data point set or from the fitting curve, it can be established that and ℎ are in a monotonous and nonlinear inverse proportional relationship: as ℎ increases, the deceleration rate of declines and shows a trend toward zero. Based on Table 10, the explicit simulation results are better than the implicit simulation results. The discrete data point set (ℎ , , ) given by the explicit simulation is shown in Figure 8. From the figure, it is easily derived that the law of explicit simulation data is the same as that of experimental discrete data, especially when ℎ is large. Herein, it is considered that this phenomenon has a certain scientific nature. As it is widely known, the accuracy of the approximate solution to the mathematical and physical problems (in this paper, the solid static problem) given by the structural finite element mainly depends on two factors: the ability to accurately describe the boundary conditions and the ability to accurately describe the general law of object behavior. Therefore, the explicit solution method provides a satisfactory solution for the finite element calculation of the shear strength of concrete beams described by the damage constitutive model.

Conclusions
In this paper, physical and simulation tests on the size effect of the shear strength of a group of recycled concrete beams without web reinforcement under the condition of a constant shear span ratio are carried out. Accordingly, the conclusions of this article are as follows: 1. From the test results, the shear strength of recycled concrete beams without web reinforcement has a size effect. In general, the shear strength is inversely proportional to the effective height of the section (ℎ )-the smaller the ℎ , the more obvious the size effect. In addition, the above-mentioned relationship can be better quantitatively described by the regression fitting formula. 2. Compared with the implicit finite element simulation results, the consistency between the explicit finite element simulation results and the experimental results is much higher. This situation not only shows that the damage constitutive model used in finite element modeling is reasonable, but it also shows that the explicit algorithm for solving this model is also reasonable. The continuous curve of h 0 and τ t obtained by the regression formula is shown in Figure 8; it is clearly seen that the fitting effect is qualitatively better. Either from the discrete data point set or from the fitting curve, it can be established that τ t and h 0 are in a monotonous and nonlinear inverse proportional relationship: as h 0 increases, the deceleration rate of τ t declines and shows a trend toward zero. Based on Table 10, the explicit simulation results are better than the implicit simulation results. The discrete data point set (h 0 , τ s,e ) given by the explicit simulation is shown in Figure 8. From the figure, it is easily derived that the law of explicit simulation data is the same as that of experimental discrete data, especially when h 0 is large. Herein, it is considered that this phenomenon has a certain scientific nature. As it is widely known, the accuracy of the approximate solution to the mathematical and physical problems (in this paper, the solid static problem) given by the structural finite element mainly depends on two factors: the ability to accurately describe the boundary conditions and the ability to accurately describe the general law of object behavior. Therefore, the explicit solution method provides a satisfactory solution for the finite element calculation of the shear strength of concrete beams described by the damage constitutive model.

Conclusions
In this paper, physical and simulation tests on the size effect of the shear strength of a group of recycled concrete beams without web reinforcement under the condition of a constant shear span ratio are carried out. Accordingly, the conclusions of this article are as follows:

1.
From the test results, the shear strength of recycled concrete beams without web reinforcement has a size effect. In general, the shear strength is inversely proportional to the effective height of the section (h 0 )-the smaller the h 0 , the more obvious the size effect. In addition, the above-mentioned relationship can be better quantitatively described by the regression fitting formula.

2.
Compared with the implicit finite element simulation results, the consistency between the explicit finite element simulation results and the experimental results is much higher. This situation not only shows that the damage constitutive model used in finite element modeling is reasonable, but it also shows that the explicit algorithm for solving this model is also reasonable.
This work will further promote the research of low-carbon-emission recycled concrete structures and the application of economic explicit finite element methods in static analysis of concrete structures.