1. Introduction
The splitting test is primarily utilized for measuring the tensile strength of brittle materials. Compared to traditional direct tension tests, it offers the advantages of simplicity in operation and minimal deviation coefficients. The types of specimens for splitting tests can be selected as circular or rectangular sections; the circular section is known as the Brazilian splitting test. To ensure reliable results in the Brazilian splitting test, it is crucial to maintain uniform tensile stress on the splitting plane and avoid eccentric loading [
1,
2,
3]. This test demands high standards in the collection and fabrication of test blocks, as well as in the application of loads. The rectangular section splitting test is more convenient in terms of specimen preparation and loading compared to the circular section splitting test. Additionally, the shear-strength parameters such as cohesion (c) and the angle of internal friction (φ) of concrete can be assessed based on geotechnical strength criteria to determine the safety of stress–strain conditions in geotechnical engineering. These parameters are typically required with specific values in numerical modeling. The existing determination of concrete shear-strength parameters is generally conducted through triaxial tests, where the Mohr–Coulomb law is applied to fit and calculate the data [
4,
5]. This process is more complex than the splitting test and requires specific rock-direct shear equipment. If it were possible to obtain the cohesion (c) and the angle of internal friction (φ) of concrete materials through common testing equipment like the rectangular section splitting test, it would significantly simplify the process.
Currently, to determine the shear-strength parameters of geotechnical materials, numerous scholars have proposed various methods for measuring the shear-strength parameters of such materials. Labuz, J.F. [
6], based on the Mohr–Coulomb failure criterion, derived formulas for calculating the cohesion (c) and the angle of internal friction (φ) by solving the stress envelope of triaxial tests and verified these through experiments; Seung-Jun Kwon [
7] proposed an integrated model for the shear strength of concrete interfaces from the perspective of the plastic limit theory of concrete. This model takes into account the influence of additional axial stress and lateral reinforcement on the shear friction at the interface shear crack, and its reliability was verified by comparing it with previous empirical formulas and 103 different push-off test samples from various literature sources; Lelovic, S. [
8] introduced a novel method for determining cohesion at preset shear deformation angles, which involves first measuring the indirect tensile strength of the material, then calculating the Mohr circle values and determining the cohesion as the shear stress value on the Mohr circle where the normal stress is zero. In the numerical computation section, the shear friction angle from the Mohr–Coulomb (MC) model was used, and the ultimate load theory was applied to the Brazilian test model, matching experimental results with numerical calculations to verify the applicability of this method in numerical analysis; Komurlu, E. [
9] used the double shear jaws (DSJ) test to determine the cohesion values of different types of rock materials through experimental and numerical analysis and investigated the effects of different test conditions such as jaw-plate size, three-block gap, and jaw-plate contact conditions on the cohesion test results and validity; Lelovic, S [
10] proposed a straightforward method for determining two material parameters of concrete, cohesion (c) and internal friction angle (φ), based on indirect tension tests and a theoretical framework. Cohesion (c) and internal friction angle (φ) were derived from experimental results of concrete deformation within the realms of both elastic and plastic theories. Numerical analysis was conducted using the finite element method, and the boundary value problems were solved using a step-by-step integration approach to validate the experimental outcomes. Gong, F. [
11] conducted preset angle shear tests, direct shear tests, and triaxial compression tests on two representative rocks, marble and red sandstone. The cohesion (c) and internal friction angle (φ) under different test conditions were determined, and the stress–strain relationships, failure modes, and shear-strength parameters of the marble and red sandstone samples were analyzed and compared. Based on the test results, a rational method for estimating shear-strength parameters was proposed. Ebrahim A.S, A [
12] introduced a micro-triple-axis test (MTT) device and used it as an index test to apply triaxial loads to small cubic rock samples to estimate their shear-strength parameters. The repeatability of the results was determined by conducting 27 MTTs on 6 × 6 × 6 mm polymethyl methacrylate (PMMA) cubic samples. The results showed that the values measured by the device were repeatable, with a similarity exceeding 97%. The internal friction angles (φ) obtained from conventional triaxial tests (CTTs) on cylindrical samples and a MTT on 6 × 6 × 6 mm samples for mortar samples were similar, thereby proving the feasibility of using small-scale triaxial tests to estimate the shear-strength parameters of intact limestone.
In addition to determining material shear-strength parameters through improved testing techniques, some scholars have also used machine learning to establish predictive models for determination. Shahani, N.M [
13] developed four advanced machine learning-based intelligent predictive models for estimating the cohesion (c) and internal friction angle (φ) of rocks, using P-wave velocity, density, uniaxial compressive strength, and tensile strength as input parameters. An iterative five-fold cross-validation method was employed, with the coefficient of determination, mean absolute error, mean squared error, and root mean squared error serving as performance metrics to evaluate the optimal predictive model. The results indicated that uniaxial compressive strength and tensile strength are the parameters with the greatest influence on predicting cohesion (c) and internal friction angle (φ). The findings of this study contribute to the assessment of long-term stability and deformation of rock masses in surface and underground excavations. Fathipour-Azar, H. [
14] employed computational data analysis methods to propose a novel technical approach for measuring the strength parameters of geomaterials. The study primarily utilized Random Forest, Alternating Model Trees, and Support Vector Machine (SVM) technologies to predict the Mohr–Coulomb shear-strength parameters of geomaterials. The models were trained and tested using a published database of geomaterials, including uniaxial compressive strength, uniaxial tensile strength, and various stress conditions. The average error of the established cohesion strength prediction model was less than 1.1 during both training and testing phases, and the average root mean square error was less than 1.5 MPa, demonstrating the effectiveness of this method in predicting the shear-strength parameters of sandstone. This research outcome aids in the long-term stability and deformation assessment of rock masses in surface and underground excavations. Danial, J. A. [
15] proposed an indirect method for measuring the shear-strength parameters of shale through stratigraphic index testing. Triaxial compression tests were conducted to obtain the shear-strength parameters of 230 samples, along with relevant rock parameters. A Particle Swarm Optimization Artificial Neural Network (PSO-ANN) integrated model was established, using the results of rock index tests as inputs and shear-strength parameters as the model’s output. The determination coefficients of 0.966 and 0.944 obtained for the training and testing datasets indicate that the model can be used to predict the shear-strength parameters of shale with high precision.
However, these research contents all pose certain difficulties in terms of experimental equipment and implementation processes, and a series of challenging issues arise during the specific operation. Based on this, the present study aims to combine common splitting tests with numerical simulations to develop a method for determining the shear-strength parameters of geotechnical materials such as concrete and to provide relevant formulas. This approach can further ascertain their shear strength and simplify the process of determining the shear-strength parameters of concrete materials, while also obtaining their tensile strength. For geotechnical materials, which exhibit good compressive properties but poor tensile performance, the main causes of failure are the materials’ tensile strengths and linear strains reaching their ultimate limits. The failure theories are primarily the first strength theory, which is the maximum tensile stress theory, and the second strength theory, which is the maximum tensile strain theory [
16,
17]. The rectangular section splitting method can obtain the tensile stress considered in the strength theory and indirectly measure its shear-strength parameters. This is of reference value for evaluating material characteristics and their applicable scope, providing a basis for the verification of structural strength, stiffness, and stability.
The study on determining the shear-strength parameters of concrete using the rectangular section splitting method introduces a novel testing approach for assessing the shear-strength parameters of concrete and other analogous materials. This research is of significant theoretical importance and practical value for the further refinement of mechanical performance testing for concrete materials. Consequently, it serves as a reference for researchers both domestically and internationally who are involved in the investigation of splitting tensile strength.
2. The Basic Principle of the Rectangular Section Splitting Method
The shear strength of certain geological materials, such as soil and rock, is typically determined by the Mohr–Coulomb failure theory, which encompasses two parameters: the material’s cohesion (c) and the angle of internal friction (φ). These parameters are generally established through indoor triaxial tests or direct shear tests on rock. The Mohr–Coulomb theory is extensively applied to determine the shear strength of soils and is also suitable for cement-stabilized soils and granular base courses. However, triaxial tests require very expensive equipment and a significant amount of time for sample preparation and testing, making them unsuitable for routine experiments or field tests. A novel approach is proposed to determine the shear-strength parameters of concrete materials using the rectangular section splitting method, which combines the rectangular section splitting test with the unconfined compressive strength test. By employing the Mohr–Coulomb theory, this method aims to determine the cohesion and the angle of internal friction of concrete materials.
The Mohr–Coulomb theory is represented by Mohr’s circles, where a series of circles depicting the stress state at the initial stage of failure under increasing confining pressure are plotted. The common tangent to these Mohr’s circles is then drawn, and the cohesion (c) and the angle of internal friction (φ) of the material are determined from this tangent. By combining the splitting test and the unconfined compressive strength test, Mohr’s circles are constructed as shown in
Figure 1.
where
σc is the unconfined compressive strength of concrete,
σx is the tensile stress at the center of the concrete specimen, σ
y is the compressive stress at the center of the concrete specimen, c is the cohesion stress of the concrete material, and φ is the internal friction angle of the concrete.
Huang, Y.G. [
18] proposed that during the splitting process of rock-like materials, there is a corresponding relationship between the tensile stress and compressive stress at the center of the specimen, as well as between the splitting strength of the specimen and the horizontal tensile stress at the center. Assuming the relationship between the splitting strength (
σIDT) of the specimen and the horizontal tensile stress (
σx) at the center, with a modification factor of
k*, the equation is as follows:
The relationship between the horizontal tensile stress (
σx) and the compressive stress (
σy) at the center of the specimen is:
where
σIDT is the splitting strength of concrete,
is the influence coefficient between the splitting strength (
σIDT) of the specimen and the horizontal tensile stress (
σx) at the center, and
is the influence coefficient between the horizontal tensile stress (
σx) and the compressive stress (
σy) at the center of the specimen.
Equation (6) is substituted into Equation (3) to solve, and then:
Utilizing the aforementioned theory to determine the cohesion (c) and the angle of internal friction (φ) of the material, it is only necessary to measure the unconfined compressive strength and the splitting strength of the specimen to obtain the shear-strength parameters. This method is simple and efficient, and it does not require complex triaxial testing equipment or other specialized devices.
3. Stress Relationship in Splitting Test Specimen
From the theoretical derivations mentioned above, it can be deduced that to obtain the shear-strength parameters of concrete materials through the rectangular section splitting method, one must first establish the relationship between the tensile stress and compressive stress at the center of the specimen and the specimen’s splitting strength. That is, it is necessary to determine the correction factor (
) for the relationship between the splitting strength (
σIDT) and the horizontal tensile stress at the center (
σx), as well as the influence factor (
) for the relationship between the horizontal tensile stress at the center (
σx) and the compressive stress (
σy). The determination of
and
will now be discussed. SłowiK, M. [
19] conducted a study on the influence of cylindrical specimen length on the tensile strength of concrete determined by the Brazilian splitting test. The experiment involved two different sizes of cylindrical specimens with a diameter of 150 mm and lengths of 150 mm and 300 mm, respectively, to test the tensile strength of concrete. Large-scale statistical analysis was performed using the Statistica(version 13, StatSoft Polska Sp. z o.o, Krakow, Poland) software. The results indicated that the shape and size of the specimens have a certain impact on the splitting tensile strength. Yong Yu [
20] utilized discrete element simulations to investigate the effects of different water–cement ratios on the compressive strength, elastic modulus, peak strain, and failure modes of Recycled Aggregate Concrete (RAC). The conclusions drawn were that RAC is influenced by size effects, albeit not as significantly as Normal Aggregate Concrete (NAC) with the same water–cement ratio, but more so than mortar. Deng Huafeng [
21] carried out splitting tensile strength tests on over 100 sandstone disc specimens with varying thickness-to-diameter ratios. The test results demonstrated that as the thickness-to-diameter ratio decreased, the tensile strength gradually increased. When the thickness-to-diameter ratio was approximately less than 0.3, the tensile strength increased and approached a relatively stable value. Three-dimensional finite element models with different thickness-to-diameter ratios were also established to study the distribution patterns of equivalent stress along the central axis of the discs in detail. The results showed that the maximum tensile stress along the central axis of the specimens occurred at the center of the specimen’s end faces, and the larger the thickness-to-diameter ratio of the specimens, the greater the tensile stress at the midpoint of the disc’s end face. Ince, R. [
22] conducted Brazilian splitting tests on cubic and diagonal cubic specimens with similar shapes but different sizes, considering various load distribution methods. The test results were analyzed using a modified size effect law, a predictive formula was proposed, and historical test data were compared to ascertain that the size effect of concrete and the loading method have a significant influence on the test outcomes.
Equation (5) assumes a coefficient relationship between compressive stress and tensile stress, denoted as
. Based on existing research outcomes [
23,
24,
25,
26], it can be derived that for disk-shaped specimens, the compressive stress at the center is three times the tensile stress. For rectangular cross-section specimens, this ratio varies and is primarily influenced by the specimen’s loading width (i.e., ratio of the bearing strip), the aspect ratio of height to length, and the thickness. To obtain an expression for the influence coefficient
, the effects of loading width, aspect ratio, and thickness on the correction coefficient were considered. Finite element software(Abaqus/CAE 2022) was utilized for simulation, and fitting curves were further plotted to provide the corresponding formulas for each factor’s impact on
.
3.1. Numerical Modeling
The numerical simulation process is conducted using the Abaqus/CAE 2022 finite element software for splitting tests. The material selected for analysis is C35 concrete. The applied load is in the form of a distributed load, with a design value of 20 kN, directed vertically downward. The specific load application is as shown in
Figure 2a. The finite element model uses a cubic specimen with dimensions of 150 mm × 150 mm × 150 mm, with a pad width of 30 mm; thus, the pad ratio is set to 0.2. The pad is modeled as a rigid body, and the interaction between the pad and concrete is achieved through contact. The tangential behavior of the contact is set to rough contact, while the normal behavior is set to hard contact. The Poisson’s ratio for concrete is taken as 0.2, and the elastic modulus is set to 3.0 × 10
3 MPa. The Mohr–Coulomb elastoplastic model is adopted for the material. The material density is 2200 kg/m
3. The cohesion of the material is set to 4.81 MPa, and the internal friction angle is set to 55°. The element type used is C3D20, and the finite element model is as shown in
Figure 2b.
3.2. Effect of Bearing Strip Ratio on ks
Considering the influence of the distributed load width on the splitting tensile strength (
), i.e., the effect of the pad ratio on
, common splitting test loading conditions are simulated numerically. In the numerical simulation process, the pad ratio is set to 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, and 0.5, corresponding to pad widths of 15 mm, 22.5 mm, 30 mm, 37.5 mm, 45 mm, 60 mm, and 75 mm, respectively. The numerical simulation results extract the numerical relationship between the compressive stress and tensile stress ratio at the center of the specimen, as shown in
Table 1.
The data in
Table 1 were drawn and fit, and the fitting results are shown in
Figure 3.
From the fitting of
Figure 3, it can be observed that the stress ratio at the center of the specimen has a quadratic relationship with the strip ratio
a/l. The stress correction coefficient
due to the loading width is hypothesized to have the following expression:
3.3. Effect of Specimen Height-to-Length Ratio on ks
Considering the influence of the aspect ratio on the stress distribution within the specimen, based on common splitting test loading conditions, the aspect ratio is set to 0.25, 0.5, 0.75, 1.0, 1.5, and 2.0 in the numerical simulation process. This corresponds to specimen heights of 37.5 mm, 75 mm, 112.5 mm, 150 mm, 225 mm, and 300 mm, respectively. The relationship between the compressive stress and tensile stress ratio at the center of the specimen is extracted and presented as shown in
Table 2.
The data in
Table 2 were drawn and fit, and the fitting results are shown in
Figure 4.
From
Figure 4, it is evident that the stress ratio at the center of the specimen has a complex relationship with the height-to-length ratio h/l, generally following an exponential curve. When the height-to-length ratio approaches 1.0, the stress ratio at the center of the specimen is at its minimum. The stress correction coefficient
, due to different height-to-length ratios of the specimens, is hypothesized to have the following expression:
3.4. Effect of Specimen Thickness on ks
Considering the impact of thickness on the stress distribution within the specimen, based on common splitting test loading conditions, the thickness-to-length ratio is set to 0.25, 0.5, 0.75, 1.0, 1.5, 2.0, 2.5, and 3 in the numerical simulation process. This corresponds to specimen thicknesses of 37.5 mm, 75 mm, 112.5 mm, 150 mm, 225 mm, 300 mm, 375 mm, and 450 mm, respectively. The relationship between the compressive stress and tensile stress ratio at the center of the specimen is extracted and presented as shown in
Table 3.
Table 3 is drawn, and the data in it are fitted, and the fitting results are shown in
Figure 5.
From
Figure 5, it can be observed that the relationship between the stress ratio at the center of the specimen and the thickness is quite complex. The stress ratio is maximum when the thickness-to-height ratio is 1. The rate of decrease in the stress ratio varies with different thicknesses. The stress correction coefficient
due to the specimen thickness is hypothesized to have the following expression:
Taking into account all the influencing factors, the correction coefficient
for the relationship between the compressive stress and tensile stress at the center of the splitting test specimen can be expressed as:
3.5. Determination of Correction Coefficient k*
When conducting splitting tests on specimens, the tensile stress generated at the center of the specimen under the action of external loads exhibits a certain relationship in materials of the same nature and under the action of external loads. The value of this tensile stress (
σx) can be expressed as:
where
σIDT is the splitting tensile strength of old and new concrete (MPa), P is tensile failure load (N), and A is the area of the specimen’s splitting surface (mm
2).
The relationship between the splitting strength (
σIDT) of the specimen and the horizontal tensile stress at the center, considering various influencing factors, is assumed to have a correction factor
. Numerical simulations are conducted by altering the magnitude of different applied loads, and the calculated results are presented in
Table 4.
From
Table 4, it can be observed that under different loading conditions, the relationship between the tensile stress at the center of the specimen and the splitting strength remains approximately constant.
can be approximated as 0.885.