Stress Space Research Articles

Assessment of interaction the new construction object at full load and the surrounding building in the dense urban development conditions

The work is devoted to research the impact of a new building at full load after completion of its construction has been investigated through an assessment of the interaction with the building envelope, the ground mass of the base and the foundation of the existing five-story building, in comparison with the assessment of the impact of a deep excavation in the process of engineering preparation for this construction, based on the developed methodology for studying the interaction of new construction objects and nearby buildings in complex engineering and geological conditions and dense development conditions.
 For numerical studies of the stability of the combined space, the proposed methodology has been used, which is built on the initial relations of the variational equations of equilibrium and the equation of the load surface in the six-dimensional stress space, which is a development of the theory of the limit stress state of the soil half-space based on the introduction of an extended yield criterion for a plane problem of the nonlinear theory of elasticity and plasticity, which provides for the determination of the value of the second critical load, at which solid sections of the ultimate stress state appear in the soil half-space.
 From the analysis of the movements of the nodes located on the vertical face of the «wall in the soil» structure (retaining wall of the pit), it can be concluded that the transverse displacements of the retaining wall compared to the calculations according to the first option have changed to the opposite, which is associated with the formation of a sliding prism under the passive pressure on the retaining wall from the side of the pile foundation of the new building. This indicates a change in the stress-deformation state of the soil base of the existing building, from the influence of the deep excavation during the construction process, which was evaluated according to the first version of calculations, to the impact of the new building after the end of construction, which was analyzed according to the second version of calculations. In this way, the change in the stress-deformation state of the main urban building is determined, which is a criterion for determining the impact of adjacent construction in the urban area on the condition of the ground bases and foundations of adjacent buildings, their preservation and ensuring safe operation.

Double image stress point bounding surface model for monotonic and cyclic loading on anisotropic clays

In this paper, an anisotropic non-associative constitutive model, based on the well-established critical state soil mechanics framework, is developed for the natural and reconstituted clays subjected to monotonic and cyclic loadings. A robust non-elliptical yield surface is implemented in the model which enables it to capture a wide variety of yield stress points. In addition to a realistic soil stiffness simulation, the employed flexible plastic potential surface equips the proposed model to predict a more accurate coefficient of earth pressure at rest. The model uses a combination of conventional volumetric and comprehensive rotational hardening rules to control the evolution of the yield surface caused by plastic strain increments. The adopted rotational hardening rule uses the concept of governing plastic strain increment which not only ensures the existence of the equilibrium state of anisotropy but also takes the effect of plastic strains at different constant stress ratios into account more sensibly. The proposed model is enhanced with a novel double-image stress point bounding surface plasticity type theory to capture nonlinear behaviour inside the yield surface, and also to simulate the cyclic responses. The detailed model formulation is discussed in both triaxial stress space and the e-mathrm{ln}, p plane, and the important features of the model are elaborately explored. The capabilities of the model are also demonstrated by an extensive sensitivity analysis. In the end, the model is used to carry out simulations of monotonic and cyclic element tests on different clays and the results are compared with the available experimental data.

Experimental study of dynamic biaxial compressive properties of full grade aggregate concrete after freeze thaw cycles

Most hydraulic concrete structures are in complex stress state, this study attempts to take some test of full grade aggregate concrete under biaxial compression after 0, 25, 50, 75, and 100 cycles of freeze-thaw. Five principle compression stress ratios of 1:0, 1:0.25, 1:0.5, 1:0.75 and 1:1 and four different strain rates of 10−5/s, 10−4/s, 10−3/s and 10−2/s were the main variables. The mass loss of concrete after freeze-thaw cycles has been studied. With the increase of freeze-thaw cycles, the quality loss is little. The failure pattern, ultimate strength, peak strain and stress-strain curve under different freeze-thaw cycles, stress ratios and strain rates were obtained. The dependence of ultimate strength on freeze-thaw cycles and ultimate strength, and the relationship between freeze-thaw cycles, stress ratio and strain rate were analyzed according to the testing data. The unified failure criterion of full grade aggregate concrete after freeze-thaw cycles under dynamic biaxial compressive stress state was established in principal stress space and octahedral stress space, respectively. The comparative analysis of experimental results of different types of concrete(large aggregate concrete, ordinary concrete, lightweight aggregate concrete, wet screened concrete) under biaxial compression after freeze-thaw cycles are made. With the increase of freeze-thaw cycles, the ultimate compressive strength of common concrete, wet screened concrete, lightweight aggregate concrete and large aggregate concrete gradually decreases. The results may serve as an experimental and theoretical reference for the design and maintenance of hydraulic structures.

Failure of substrate-supported hexagonal boron nitride under complicated loading conditions

Buckling and fracture are two crucial failure modes of h-BN sheet and establishing the criterion for the failure modes is the prerequisite to guarantee the proper applications of h-BN sheet in practice. For the failure of h-BN sheet, previous studies mainly focused on the freestanding h-BN sheet under uniaxial or equi-biaxial loadings. However, for the more general case in practice, the critical failure of substrate-supported h-BN sheet under complicated loadings remains to be studied. In this paper, the substrate-supported h-BN sheet under complicated loadings is considered to reveal the failure mechanism and to establish the criteria for both the buckling and fracture. In order to predict the critical failure stresses efficiently and accurately for complicated loadings, a modified mode-independent energy-based analysis method (mMIEM) is proposed. Then, by analyzing the stress state and observing the failure features, the rational criteria for the buckling and fracture are established and the space of the allowable stress for the SiO2 substrate supported h- BN sheet is provided. The results indicate that the buckling of h-BN sheet is dominated by the second principal stress, namely the maximum compressive stress, while the fracture depends on the stresses along and vertical to the armchair direction rather than the principal stresses. Moreover, the substrate has a significant effect on the critical buckling of the h-BN sheet but has little effect on the critical fracture. The failure criterion and allowable stress space can be used in the design and evaluation of the h-BN based electronic devices.

Genesis of Bedding Fractures in Ordovician to Silurian Marine Shale in Sichuan Basin

The effective utilization of shale bedding fractures is of great significance to improve shale gas recovery efficiency. Taking the Wufeng–Longmaxi Formation shale in Sichuan Basin as the research object, the formation process and mechanism of bedding fractures in marine shale are discussed, based on field observation and description, high-resolution electron microscope scanning, fluid inclusion detection, and structural subsidence history analysis. The results show that the formation of bedding fractures is jointly controlled by sedimentary characteristics, hydrocarbon generation, and tectonic movement: the development degree of bedding (fractures) is controlled by the content of shale organic matter and brittle minerals, and bedding fractures formed in the layers with high organic matter; tectonic movement created stress environment and space for bedding fractures and promoted the opening of bedding fractures; the time for calcite vein to capture fluid is consistent with the time of oil-gas secondary pyrolysis stage. The formation of the calcite vein is accompanied by the opening of fractures. The acid and oil-gas generated in the hydrocarbon generation process occupied the opening space and maintained the bedding fractures open. The study of the formation process of bedding fractures is helpful to select a suitable method to identify bedding fractures, and then effectively use it to form complex fracture networks in the fracturing process to improve shale oil and gas recovery.

A stress-based shear fracture criterion considering the effect of stress triaxiality and Lode parameter

Failure in lightweight metal forming is a big challenge currently and prevents their widely application in weight reduction of automobile and aerospace structures. A stress-based shear ductile fracture criterion is introduced in this study to improve the prediction accuracy of failure for lightweight metals. The criterion (sDF2016) is developed based on the DF2016 criterion thereby endowing the criterion with the precise fracture predictability under wide stress states of shear, uniaxial tension, plane strain tension and equibiaxial tension. Besides, the criterion takes into account the stress state effect on fracture in a form of the stress triaxiality, the maximum shear stress and the Lode parameter. The sDF2016 criterion is also expected to be less sensitive to strain path changing effect because sDF2016 describe the onset of fracture in stress space. For the verification of the proposed sDF2016 criterion, the experiments are carried out for AA5182-O sheet with central hole, notched, in-plane shear and bulging specimens. Plastic deformation is accurately modeled by the Swift-Voce hardening law and the Drucker yield function. The sDF2016 criterion is then calibrated both by the direct and inverse engineering approaches, which is applied to predict fracture initiation under distinct stress states. The modeled result indicates that the calibrated sDF2016 criterion based on the inverse engineering method predicts the fracture stroke with a higher accuracy than the direct approach. The sDF2016 criterion is also used to depict fracture limits of AA2024-T351 alloy from compression to tension. The application shows that the proposed ductile fracture criterion is capable of modeling the fracture limits for sheet metals under various loading conditions from shear to plane strain tension.

Coupled hydro-mechanical modelling of unsaturated soils; numerical implementation and application to large deformation problems

This paper presents coupled hydro-mechanical modelling of unsaturated soil problems by incorporating some advanced numerical and constitutive models in a general-purpose commercial software package, Abaqus. Two different strategies for the interpretation of unsaturated soil behaviour with respect to the soil volume change are considered and the relevant constitutive models are implemented through user-defined subroutines. The first approach that is commonly used in geomechanics considers suction as an additional variable and uses a constitutive model in the space of effective stress–suction assuming that soil compressibility is a function of suction. The second approach treats the constitutive model in the space of effective stress and degree of saturation assuming that the soil compressibility is a function of the degree of saturation. The hydro-mechanical behaviour regarding the change in the degree of saturation caused by both suction and net stress changes is also considered, together with the effect of hydraulic hysteresis. Validation of the implemented algorithms is presented through several benchmark problems. Finally, the application and utility of the implemented procedures are illustrated by simulations of two challenging problems of unsaturated geomechanics, including a slope failure due to seepage and rainfall infiltration as well as cone penetration tests in unsaturated soil. A suitable mesh optimisation scheme is also incorporated to handle finite deformations.

Hodge duality between stress space and strain space in anisotropic media

AbstractDuality in a linear elasticity theory is studied based on a theory of differential form. A generalized expression of a Hodge star operator with an index is introduced. The index in the Hodge star operator means a superposition of ordinary Hodge star operators. By using the superposed Hodge star operator, a linear constitutive relation in anisotropic media can be expressed as a duality between a stress and strain spaces. Then, a basic elastodynamic equation is derived from the conservation law and the linear constitutive relation in the stress and strain spaces. Moreover, a geometric expression of equation of stress function and displacement function is derived by using the differential form and the dual relation of linear constitutive relation. These geometric results imply that the approach of differential forms is applicable to an analysis of deformation in anisotropic media. The superposed Hodge star operator discussed in this study has similar properties of a discrete Hodge operator, which represents the summation of differential forms on the decomposed regions. This means that the superposed Hodge star operator and the discrete Hodge operator express the constitutive relation in which physical quantities act on each other.

The effect of Lode parameter on the yield surface of Schoen's IWP triply periodic minimal surface lattice

Owing to the advancements in additive manufacturing and increased applications of additively manufactured structures, it is essential to fully understand both the elastic and plastic behavior of cellular materials, which include the mathematically-driven sheet lattices based on triply periodic minimal surface (TPMS) that have received significant attention recently. The compressive elastic and plastic behaviors have been well established for many TPMS latticed structures, but not under multiaxial loading. Furthermore, TPMS lattices are computationally expensive to model explicitly when used in latticing various structures for enhanced multi-functionality, and hence the need to develop accurate yield surfaces in order to capture their plastic behavior in a homogenized approach. The majority of previous yield surfaces developed for cellular materials originate from cellular foams, and limited attempts has been made to develop yield surfaces for TPMS lattices. In this study, a numerical modeling framework is proposed for developing the initial yield surface for cellular materials and is used to develop the initial yield surface for Schoen's IWP sheet-based TPMS cellular lattices. The effect of different loading conditions on the effective yield strength of the IWP sheet-based (IWP-s) TPMS lattice is numerically investigated, based on a single unit cell of IWP-s under fully periodic boundary conditions, assuming an elastic-perfectly plastic behavior of the base material, for relative densities (ρ‾) ranging from 7% to 28%. In order to account for different loading conditions, the stress state is characterized in a generalized fashion through the Lode parameter (L). The effect of L is studied over a range of mean stress values (σm) to understand the effect of both L and σm on the effective yield strength. An initial yield surface is developed incorporating the effect of L, σm and ρ‾, and is validated numerically showing rather good agreement. In the 3D principal stress space, the shape of the yield surface for the IWP-s lattice resembles the shape of a cocoa pod.

Some Issues with Statistical Crystal Plasticity Models: Description of the Effects Triggered in FCC Crystals by Loading with Strain-Path Changes.

The justification of the applicability of constitutive models to exploring technological processes requires a detailed analysis of their performance when they are used to describe loadings including the complex loading mode that is characteristic of these processes. This paper considers the effect of equivalent stress overshooting after the strain-path changes known to occur in metals and alloys. The macrophenomenological and multilevel models, which are based on crystal plasticity, account for this effect by applying anisotropic yield criteria at the macro- and mesolevels, respectively. We introduce a two-level constitutive statistical inelastic deformation model (identified for aluminum) that incorporates the popular simple phenomenological anisotropic hardening law for describing the behavior of FCC polycrystals. The results of the numerical simulation are in satisfactory agreement with existing experimental data. Statistical analysis of the motion of a mesostress in the stress space on the crystallite yield surface is performed. The obtained data are compared with the results found using the isotropic hardening law. The results clarify the simulation details of statistical crystal plasticity models under loading with strain-path changes in materials and demonstrate their suitability for describing the processes under consideration.

Geometric prior of multi-resolution yielding manifolds and the local closest point projection for nearly non-smooth plasticity

Elastoplasticity models often introduce a scalar-valued yield function to implicitly represent the boundary between elastic and plastic material states. This paper introduces a new alternative where the yield envelope is represented by a manifold of which the topology and the geometry are learned from a set of data points in a parametric space (e.g. principal stress space, π-plane). Here, deep geometric learning enables us to reconstruct a highly complex yield envelope by breaking it down into multiple coordinate charts. The global atlas that consists of these coordinate charts in return allows us to represent the yield surface via multiple overlapping patches, each with a specific local parametrization. This setup provides several advantages over the classical implicit function representation approach. For instance, the availability of coordinate charts enables us to introduce an alternative stress integration algorithm where the trial stress may project directly on a local patch and hence circumvent the issues related to non-smoothness and the lack of convexity of yield surfaces. Meanwhile, the local parametric approach also enables us to predict hardening/softening locally in the parametric space, even without complete knowledge of the yield surface. Comparisons between the classical yield function approach on the non-smooth plasticity and anisotropic cam-clay plasticity model are provided to demonstrate the capacity of the models for highly precise yield surface and the feasibility of the implementation of the learned model in the local stress integration algorithm.

Macromechanical Failure Criteria: Elasticity, Plasticity and Numerical Applications for the Non-Linear Masonry Modelling

Sometimes it is difficult to choose the most appropriate failure criterion for the problem analyzed. For brittle materials, attention must be paid to the availability of experimental data and the calibration of the representative parameters, within the chosen failure criterion. The work herein presented, starting with an overview on machromechanical failure criteria, analysed in the Haigh-Westergaard Stress Space, investigates the suitability of Mohr-Coulomb, Drucker-Prager and Concrete Damaged Plasticity failure criteria of masonry structures, underlining their specific characteristics and implementation in FEM simulations. The Pavia Door Wall experimental campaign under pseudo-static cyclic test is considered as benchmark study. The results of the experimental tests are compared with a FE model developed with ABAQUS computer program considering several failure criteria and equivalent frame approach. Among the investigated failure criteria Concrete Damaged Plasticity is able to capture the actual behaviour of the masonry walls under monotonic excitation. In particular, thanks to the adaptability of the Guo’s model in the definition and calibration of the uniaxial behavior, the model suitability in catching the variation of the cohesion and the evolution of the damage is better in comparison with the other addressed failure criteria.

Interpretation of mechanical and biodegradation behaviour of municipal solid waste

An advanced constitutive model for MSW (municipal solid waste) is proposed in this study to consider biodegradation and overconsolidation jointly based on sub-loading surface plasticity. A reference line, namely the BNCL (biodegradation-independent normally compression line), is introduced to quantify the total biodegradation time and overconsolidation state of MSW. Under unchanged stress level, the void ratio change due to biodegradation is interpreted as the change of overconsolidation. An isotropic stress–strain-biodegradation constitutive model is then proposed by considering biodegradation-dependent deformation through overconsolidation parameter. The current yield surface and the reference yield surface are introduced for assessing the effect of overconsolidation (biodegradation has been included). The isotropic model is then extended to a triaxial stress state in the space of mean effective stress, deviator stress and biodegradation variable. The performance of the constitutive model is discussed through various scenarios and validated by using experimental data from the literature.

Experimental analysis of the multiaxial failure stress locus of commercially pure titanium at low and high rates of strain

The mechanical response and failure mechanism of commercially pure titanium subjected to combined tension-torsion loading are studied experimentally at strain rates ranging from 10−3s−1 to 103s−1. A novel tension-torsion split Hopkinson bar (TTHB) equipped with a high speed camera was utilised during high-rate experiments, while quasi-static tests were conducted using a universal screw-driven machine. The multiaxial dynamic experiments demonstrate the ability of the developed TTHB apparatus to achieve synchronisation of longitudinal and torsional waves upon loading the specimen, to satisfy the dynamic equilibrium of the specimen and to attain constant strain rate loading. The failure envelope of commercially pure titanium was analysed over a wide range of stress states including pure torsion, shear-dominated combined tension-shear, tension-dominated combined tension-shear, and plain tension. The analyses of the loading paths show that these were nearly proportional in terms of strain. The multiaxial failure stress locus was constructed in the normal versus shear stress space from experiments conducted at low and high rates of strain.The Drucker-Prager criterion was employed to approximate the failure envelope and to assess its rate sensitivity. The failure stress locus of commercially pure titanium and its rate dependence are reported for the first time. The TTHB apparatus developed allows the definition of the failure stress locus of aerospace materials directly from experiments and, therefore, the evaluation of the existing failure/yielding criteria.

Modeling cyclic behavior of clay subjected to principal stress rotation

Experimental evidence indicates that principal stress rotation yields complex cyclic behavior in clay, which is significant in many cases but not adequately described by constitutive modeling. To address this issue, a constitutive model was developed to describe the cyclic behavior of clay within the framework of bounding surface elastoplasticity. To describe the plastic mechanism induced during reversal loading, the mapping rules for the relocatable projection center were formulated in the deviatoric stress space, which are also applicable to pure principal stress rotation. The effects of small strain stiffness and anisotropic elasticity were incorporated into the model to describe cyclic degradation and plastic accumulation behavior under principal stress rotation. To describe the non-coaxiality, a three-dimensional non-coaxial flow rule was incorporated into the model. It was assumed to be colinear with a non-coaxial stress rate orthogonal to a reference stress tensor and dependent on the stress ratio. The developed model was validated against the undrained hollow cylindrical torsional shear tests subjected to cyclic pure principal stress rotation for Shanghai and Hangzhou clay with various intermediate principal stress coefficients. The comparison between simulations and experimental results demonstrates that the proposed model can model the cyclic behavior of clay subjected to principal stress rotation.

Energy-based universal failure criterion and strength-Poisson's ratio relationship for isotropic materials

The prediction of material failures has been one of the most important topics in engineering design since the onset of human civilization. To date, a plethora of failure criteria have been developed, particularly for homogeneous isotropic materials. However, many of them have limited usage for specific groups of materials, are obscure in their origin, or different from the actual phenomena at the microscopic and macroscopic scales. Therefore, based on the most fundamental concept of energy, this study aims to establish a universal failure criterion for homogeneous full-density isotropic materials. The developed failure criterion requires three independent load-bearing capabilities of the material, that is, tensile (T), compressive (C), and shear (S) strengths, as criterion calibrators for applicability to a wide range of materials, including those with extremely ductile to exceedingly brittle behaviors. Energy and polynomial invariants expansion approaches were used to formulate the universal failure criterion, and both yielded a perfect dual relationship. The criterion in the energy form revealed a new physical insight into the zero-energy bound or fracture criterion, which corresponds to an intrinsic cut-off plane in a failure envelope. Comparisons between the experimental data across a broad spectrum of material types and the proposed failure criterion were in excellent agreement in both the stress and energy spaces. The strength-Poisson's ratio relationship was also derived for the first time, demonstrating that only two strength properties are independent and that three important bounding conditions for isotropic materials exist: (1) S/C = 0 when T/C = 0; (2) when T/C = 1 and S/C = 1/3, Poisson's ratio v = 0.5; and (3) v = −1 when T/S and C/S approach zero.

Progressive damage and failure analysis of three-dimensional braided composites under multiaxial loadings

A meso-mechanical computational model based on the finite element method (FEM) is developed to predict the progressive elastic-plastic damage and failure behavior of three-dimensional four directional (3D4D) braided composites under multiaxial loadings. Firstly, the constitutive laws for mesoscale constituents (i.e. yarn, matrix and interface) of braided representative volume element (RVE) are established and implemented by the user defined material subroutine UMAT in ABAQUS. Then, how to build the braided RVE and conduct the simulation are introduced. After that, the uniaxial tensile data is used to verify the proposed model. Finally, the damage and failure behavior of the RVE under multiaxial loadings is investigated, considering the effects of loading paths and interface. The results show that loading path and interface have little effect on the strengths. The predicted failure envelopes in different stress spaces for RVE are compared with the Tsai-Wu and Hoffman theories, which are much closer to the Tsai-Wu theory.

Elastoplastic Integration Method of Mohr-Coulomb Criterion

A new method for implicit integration of the Mohr-Coulomb non-smooth multisurface plasticity models is presented, and Koiter’s requirements are incorporated exactly within the proposed algorithm. Algorithmic and numerical complexities are identified and introduced by the nonsmooth intersections of the Mohr-Coulomb surfaces; then, a projection contraction algorithm is applied to solve the classical Kuhn–Tucker complementary equations which provide the only characterization of possible active yield surfaces as a special class of variational inequalities, and the actual active yield surface is further determined by iteration. The basic idea is to calculate derivatives of the yield and potential functions with the expressions in the principal stresses and perform the return manipulations in the general stress space. Based on the principal stress characteristic equation, partial derivatives of principal stresses are calculated. The proposed algorithm eliminates the error caused by smoothing the corner of Mohr-Coulomb surfaces, avoids the numerical singularity at the intersections in the general stress space, and does not require the stress transformation needed in the principal stress space method. Lastly, several numerical examples are given to verify the validity of the proposed method.

An elastoplastic damage constitutive model for hybrid steel-polypropylene fiber reinforced concrete

An appropriate constitutive model to characterize the mechanical behavior of fiber reinforced concrete (FRC) plays a critical role in accurate prediction of structural performance. In this study, an elastoplastic damage constitutive model is developed for hybrid steel-polypropylene fiber reinforced concrete (HFRC), which can predict the versatile mechanical responses of HFRC as a result of various hybrid fiber inclusions. For plasticity growth, the Willam-Warnke five-parameter failure envelope is modified in the effective stress space to govern the yield condition of HFRC, while a nonlinear Drucker-Prager type plastic potential function is proposed to control the non-associated plastic flow. Regarding the damage evolution, the damage criterion and evolution laws are established in the Cauchy stress space, in which the plastic hardening and unilateral behaviors are taken into consideration. Upon two parallel integration algorithms to describe the development of plasticity and damage, as well as the identification of fiber-dependent parameters, the proposed model is numerically implemented and then validated by the experimental results from available research. The comparisons between the numerical predictions and the test results at both material scale and structural scale solidly demonstrate the capacity of the model in capturing the main features of HFRC when subjected to different loading paths.

Experimental Study and Theoretical Analysis on the Compression-Shear Multiaxial Mechanical Properties of Recycled Concrete.

Recycled concrete, which is formed by replacing coarse aggregates in ordinary concrete with recycled aggregates (RA), is of great significance for the secondary utilization of waste building resources. In civil engineering, concrete structures are sometimes subjected to a compression–shear multiaxial stress state. Therefore, research on the compression–shear multiaxial mechanical properties of recycled concrete plays an important role in engineering practice. To explore the effect of RA replacement rate on the compression–shear properties of recycled concrete, an experimental study was carried out using a compression–shear testing machine and considering five RA replacement rates and five axial compression ratios. Consequently, the failure modes and mechanical property parameters under different working conditions were obtained and were used to analyze the effects of RA replacement rate and axial compression ratio on the shear stress of recycled concrete. Eventually, the following conclusions were reached: With the growth of axial compression ratio, the shear cracks exhibit a developing trend along the oblique direction, and the friction traces on the shear surface are gradually deepened. As the replacement rate increases, the number of shear cracks is gradually increased, accompanied by increasing broken fragments falling off from the shear interface. Since the action of the axial compression ratio can effectively improve the mechanical bite force and friction on the shear interface of recycled concrete, as the axial compression ratio increases, the shear stress is gradually increased. On the other hand, due to the initial damage of RA and its weak adhesion with cement mortar, the shear stress is gradually reduced with the increase of RA replacement rate. Meanwhile, the increase in shear stress shows a gradually decreasing trend with the growth of axial compression ratio. Specifically, for the RA replacement rates of 0% and 100%, the shear stress increased by 4.06 times and 3.21 times, respectively, under the influence of the axial compression ratio. Under different axial compression ratios, the shear stress was reduced by 43~46%, due to the increase of RA replacement rate. In addition, based on the octahedral stress space and the principal stress space, a compression–shear multiaxial failure criterion and shear stress calculation model for recycled concrete were proposed, by considering the effect of the RA replacement rate. The outcomes of this research are of great significance for engineering applications and the development of recycled concrete.