Development Of Law Research Articles

Mechanism of uneven densification in PBX compression molding

The uneven density distribution during compression molding of polymer bonded explosive (PBX ) seriously affects the mechanical properties of high-energy material components and the precise output ability of shock waves. Investigating the multi-scale densification evolution laws and density non-uniformity mechanism of PBX composite powders is crucial for understanding the density non-uniformity phenomenon of energetic material components and evaluating the differences in macroscopic mechanical response. To this end, we developed a new dual-scale 3D discrete element method (DEM) model for PBX compression molding, which can increase the computational quantity of macroscopic powders by about an order of magnitude while considering the true morphology of mesoscopic explosive crystals. First, the densification behavior at different scales were analyzed using the developed 3D DEM models, including the rearrangement of macroscopic powder during the low strain stage, the deformation of macroscopic powders and mesoscopic bonding failure inside the powders during the high strain stage. Furtherly, the types and percentages of microcracks at the crystal-binder interface and inside the binder at the mesoscale during the macro powder deformation in the compression densification process were quantified. The characteristics of pore distribution and evolution laws of porosity and density uniformity during compression densification were revealed by simulations and X-ray μCT scanning. Finally, the mechanism of density non-uniformity was elucidated from multiple perspectives including displacement, force chain, contact fabric, and stress distribution. Subsequently, different degrees of energy dissipation behavior and the same gradient attenuation pattern were captured at different scales. This indicates that axial stress attenuation and density inhomogeneity are co-dominant results of the gradient attenuation of different scales of energy during layer-by-layer transfer along the loading direction. The potential induced mechanism of uneven density and stress attenuation were explained for the first time from the perspective of energy dissipation by this finding.

Spatio-temporal evolution law of anisotropic shear damage on rock mass joint surface affected by joint morphology

The spatial–temporal evolution law of shear damage on the joint surface is closely related to its morphological characteristics. Understanding the heterogeneous failure behavior of the joint surface is essential for revealing the controlling role of the joint morphology. The diorite from the Guanshan Tunnel is investigated through anisotropic direct shear tests in this study. Based on shear damage quantification and acoustic emission (AE) localization techniques, the heterogeneous damage evolution law of joint surfaces under different joint compression strength (JCS), normal stress (σn), and morphology parameter (M) is studied. The results indicate that the joint surface morphology, shear damage, and AE parameters exhibit similar anisotropic characteristics in different shear directions. An increase in the morphology parameter leads to an increase in shear failure, with shear fracture energy and the number of fracture events showing an increasing trend. Meanwhile, the occurrence of a large number of fracture points shifts from the pre-peak nonlinear period to the post-peak period. However, the evolution process of fracture events exhibits certain similarities, though the morphology parameters are different in two opposite shear directions. Additionally, shear damage volume and fracture energy demonstrate consistency in spatial distribution and quantitative characteristics, while the linear scale factor between cumulative shear fracture energy and damage volume decreases with increasing joint surface strength. These findings assist in a better understanding of the anti-sliding property of joint surfaces and provide crucial clues for further exploration of joint surface failure mechanisms.

Analysis of the Stability of a High Fill Slope under Different Gradients and Precipitation Conditions

The stability problem of high fill slopes has always been a research hotspot. Its failure mechanism is complex and prominent, featuring strong concealment, a short occurrence time and great harmfulness. In this paper, the stability of a high fill slope under rainfall conditions will be studied by using indoor tests, numerical simulations, etc. The study is based on a high fill slope in Yichang City. The evolution law of high fill slope stability under the maximum rainfall condition is revealed. The results show the following: The influence of moisture content on stress–strain curves is reflected in both the curve’s shape and the peak value of deviatoric stress. Under the constraint of confining pressure, the curve decreases and the peak value of deviatoric stress decreases with the increase of moisture content at the same confining pressure. The safety factor obtained by a rigid body limit equilibrium analysis and numerical calculation indicates that the safety factor for a 30° slope meets the requirements for slope stability evaluation and remains in a fundamentally stable state. An on-site investigation suggests that surface failure and shallow failure may be primary failure modes for this slope; therefore, it is recommended to implement slope protection measures. This study provides valuable references for similar high fill slopes.

Dynamic mechanical properties of sandstone with two circular inclusions under impact loading

Inclusions are often present in rock defects, which can compensate for some of the problems of low rock strength and poor stability caused by hole defects. Specimens containing inclusions exhibit different mechanical properties under dynamic loading, and the strength change, energy evolution law and fracture mechanism deserve to be clearly discussed. In this work, a dynamic impact test of sandstone with circular inclusions was carried out. The failure process of each sample was recorded by a high-speed camera, and the influences of the inclusion strength and strain rate on the dynamic impact performance of the sandstone were analysed. Increasing the strain rate and inclusion strength can significantly improve the dynamic compressive strength of a sample. The dissipated energy density and impact toughness index are positively correlated with the strain rate. The dissipated energy density increased the most when the sample type was FA, at 10.6%. When the sample type was FD, the impact toughness index increased the most, which was 49.1%. In the early stage of loading, the localized strain is mainly concentrated inside the filler and gradually evolves into deformation of the filler and extension of shear cracks in the late stage of loading. The failure modes of a sample are not the result of a single factor but are influenced by both the strain rate and the filler. The types of cracks in the samples are mainly shear cracks and shear–tension composite cracks.

Unveiling the nexus: exploring the influence of terrorism on energy trade in China and the Belt and Road countries

An examination of the scale of energy trade between China and the countries along the Belt and Road, as well as the terrorist activities in these countries from 2008 to 2019, has been conducted in this study. The spatial and temporal evolution law is revealed through the use of exploratory spatial data analysis. The spatial Durbin model is employed to investigate the impacts and spatial spillover effects of terrorist activities along the Belt and Road on China’s energy trade with these countries. The findings of the study are as follows: Spatial clustering is exhibited by the energy scale of China’s imports and exports with countries along the routes and terrorist activities along the routes, with energy imports concentrated in West Asia and North Africa, energy exports concentrated in Southeast Asia, and terrorist activities concentrated in West Asia and North Africa, South Asia, and Central Asia. A negative effect on the scale of China’s energy imports from the countries along the route is exerted by the terrorism index, and a positive spatial spillover effect is present. The number of terrorist activities and deaths have a negative effect on the scale of China’s energy imports from the countries along the route, but there is no spatial spillover effect. The terrorism index does not affect the scale of energy exports from China to countries along the route, but there is a negative spatial spillover effect. There are differences in the spatial spillover effects of the number of terrorist activities and deaths on the scale of China’s energy exports. The level of economic development, liner transportation convenience, and energy exports dependence of the countries along the route also significantly influence China’s energy trade with these countries.

Effects of vetiver root on cracking of expansive soils and its mechanistic analysis

The study investigated the reinforcing effect of vetiver root on soil by conducting outdoor planting tests and indoor root tests. The cracking indexes of soil specimens with varying root contents were analyzed, and a statistical model was established to determine the relationship between the cracking indexes, the number of dry and wet cycles, and the root content. The study revealed the crack evolution law of vetiver-reinforced expansive soil. The study explored the mechanism of the vegetation root in inhibiting the cracking of expansive soil and determined the optimal planting density of vetiver grass through outdoor planting tests. The results indicate that: The surface crack rate (CR), total crack length (CL), and crack number (CN) in the root-soil specimen exhibited exponential growth with an increase in the number of wet and dry cycles. This growth was more pronounced during the first and second cycles. The vetiver root could effectively reduce soil crack formation, and the specimen's cracking resistance is positively correlated with the root content. With the root content increased, the CR, CN, and CL decreased. The logistic model is suited to the CL of added root soil. The logistic model is more suitable for the growth model of the CR of the expansive soil with low root content, while the Boltzmann model is more suitable for the growth model of the CR of the expansive soil with high root content. Width of crack (CW) is better suited to the DoseResp growth model. The Boltzmann model is more applicable to the CN in expansive soils with low reinforcement, while the logistic growth model is more suitable for the development of CN above 0.21% root content. The development of the crack network was influenced by two key factors: the root content and the number of wet and dry cycles. Under the condition of planting roots, the development of crack networks in expansive soil differs from that of expansive soil with added roots, and there is no clear pattern to follow. The inhibitory effect of the vetiver root on cracking of expansive soil is related to the planting density of vetiver.

Fault diagnosis of early internal short circuit for power battery systems based on the evolution of the cell charging voltage slope in variable voltage window

Internal short circuit (ISC) is one of the main causes of thermal runaway (TR) accident in power battery systems, to effectively avoid the development of early stage ISC towards TR, this paper innovatively proposes an ISC fault diagnosis method based on the evolution of the cell charging voltage slope (CCVS) in variable voltage window (VVW). Firstly, the ISC characterisation parameter, i.e. CCVS, is extracted through battery cell aging and equivalent ISC experimental studies. This parameter has two advantages: on the one hand, the charging voltage slopes of different cycles of one cell in different voltage windows keep the same evolution law, so that the mechanism of VVW can be established, and the generalisation ability of the algorithm is improved. On the other hand, the charging voltage data of one cell as the research object effectively avoids the interference of inconsistent factors between battery cells. Then, the battery data are preprocessed using wavelet denoising, and a sliding window strategy is introduced to calculate the positional ranking of the CCVS in the historical data, which represents the capacity ranking level of the cell in this charging segment, defined as the capacity rating factor (CRF). Further, a double-layer diagnostic strategy based on the 3 σ criterion and distance factor assessment is used to localise the fault, which can effectively avoid the occurrence of false alarms. Finally, experimental and real-world vehicle data validate the effectiveness of the method, different types of battery data verify the applicability of the method, and the comparison results of the same type of model show that the proposed method is significantly superior in terms of robustness, accuracy and computational efficiency.

Study on the macro-mechanical behavior and micro-structure evolution law of broken rock mass under triaxial compression

Study on the macro-mechanical behavior and micro-structure evolution law of broken rock mass under triaxial compression

Research on improving the flexural performance of alkali-activated geopolymer cemented coal gangue through layered addition of fibers

The brittle characteristics of geopolymers during failure seriously hinder their practical engineering applications. To improve the toughness of low-energy alkali-activated geopolymer cementing coal gangue (AGCCG), the flexural performance of alkali-activated geopolymer cementing coal gangue (AGCCG) was enhanced by layering polypropylene fibers. Based on DIC digital speckle technology, the four-point flexure test of AGCCG was carried out, and the strength formation mechanism of the material was analyzed in depth by combining SEM-EDS and XRD. 3D visualization and analysis of the inside of the specimen with the help of CT scanning technology. The internal crack evolution law of the sample was illustrated by PFC numerical simulation. The results showed that the crack formation threshold, peak flexural strength and residual flexural strength of the specimens reached a local maximum at 3 % NaOH doping. With the increase of NaOH dosage, the cementitious product becomes more and more dense and homogeneous, and the area of void defects in the hydration product within the field of view becomes smaller and smaller. The peak bending strength of specimen B4-#4 is 79.43 % of that of specimen B5-#1–4, in which the cost of specimen B5-#1–4 is four times of that of specimen B4-#4, and the way of fiber assignment in specimen B4-#4 is more reasonable when considered together. The layered addition of fibers has a directional strengthening effect on the mechanical properties of geopolymer, and the overall addition of fibers can improve the overall strength and toughness of geopolymer. The numerical simulation results basically coincide with the macroscopic test results, and are consistent with the progressive damage law of the specimens obtained from the numerical scattering test. The results of the study can provide a theoretical basis for improving the bearing capacity of geopolymer and optimizing its structural performance. As well as reducing the impact of solid waste on the environment and improving the efficiency of resource use.

Study on signal characteristics of burst tendency coal under different loading rates

In order to study the mechanics, acoustic emission (AE) and electromagnetic emission (EME) response law of bursting liability coal at different loading rates, uniaxial compression tests were carried out on coal mass from Konggu Coal Mine. The corresponding relations among mechanical properties, AE and EME signals in the process of coal failure under loading were analyzed, and the energy evolution law of coal failure with bursting liability under loading rate was discussed. The results show that within a certain range of loading rate, the higher the loading rate, the higher the compressive strength and peak load of bursting liability coal, and the shorter the time for coal to reach the peak load. Under different loading rates, the mechanics, AE and EME signals of coal samples can be well corresponded. When the loading rate is low, the number of blocks destroyed of coal sample is large and the block size is relatively small, and the blocks are mainly scattered around the test platform. When the loading rate is high, the number of damaged blocks is relatively small and the block size is relatively large, and the blocks are far away from the test bench. When loading at a low rate, the internal cracks in coal can be fully developed and connected, and the energy release rate is relatively uniform in the process of loading and failure of coal sample. In the case of high loading rate, the energy release rate of coal sample in the loading process is much smaller than that in the moment of failure. Combining the above test results with the actual situation of the working face, it can be concluded that the total energy stored in the coal of fast mining increases and the threshold of impact decreases compared with that of slow mining. Therefore, under the disturbance of external dynamic load, rapid mining is more likely to induce rock burst.

Study on seismic failure mechanism and damage model of unconstrained adobe block walls

Through the quasi-static cyclic in-plane loading test of three kinds of unconstrained adobe block walls, namely plain adobe wall, modified mud adobe wall and modified adobe wall, the failure mode and seismic performance of each specimen were clarified. Using residual deformation and cumulative energy dissipation as failure parameters, a two-parameter seismic damage model of unconstrained adobe block walls was established. and a calculation method for the adjustment coefficients of energy consumption items was provided. The results showed that the failure mode of the typical unconstrained adobe block wall under earthquake action was not affected by the modification of adobe materials, and the final failure mode was shear failure along the mud crack and the inclined crack extending through the block. Furthermore, the proposed seismic damage model can well reflect the evolution law of earthquake damage of adobe block walls, and the recommended damage index boundary values for five performance levels were given, providing a reference for performance-based seismic design and damage assessment of adobe block walls.

Development of coupled fluid-flow/geomechanics model considering storage and transport mechanism in shale gas reservoirs with complex fracture morphology

Field observations frequently demonstrate stress fluctuations resulting from the reservoir depletion. The development of reservoirs, particularly the completion of infill wells and refracturing, can be significantly impacted by stress changes in and around drainage areas. Previous studies mainly focus on plane fractures and few studies consider the influence of complex transport and storage mechanism and irregular fracture geometry on stress evolution in shale gas reservoirs. Based on the embedded discrete fracture model (EDFM) and finite-volume method (FVM), a coupled geomechanics/fluid model has been successfully developed considering the adsorption, desorption, diffusion and slippage of shale gas. This model achieves coupling simulation of natural fractures, hydraulic fractures with complex geometry, storage and transport mechanism, reservoir stress, and pore-elastic effect. The open-source software OpenFOAM is used as the main solver for this model. The stress calculation and productivity simulation of the model are verified by the classical poroelasticity problem and the simulation results of published research and commercial simulator with EDFM respectively. The simulation results indicate that σxx, σyy, σxy and Δσ changes with time and space due to the time effect and anisotropy of formation pressure depletion; Due to the influence of different mechanisms on shale gas storage and transport, the reservoir pressure and stress distribution under different mechanisms are different; Among them, the stress with full mechanisms differs the most compared to the stress without any mechanism. The reservoir with stronger stress sensitivity (smaller Biot coefficient) is less sensitive to formation pressure depletion, and the stress variation range is smaller. For reservoirs with weak stress sensitivity, formation pressure depletion is more likely to lead to stress reversal. Under the influence of fracture geometry, the pressure depletion regions caused by the three types of fracture geometry are approximately rectangular, parallelogram and square, respectively. The corresponding σxx, σyy and Δσ also have great differences in spatial distribution and values. Therefore, the time effect, shale gas storage and transport mechanism and the influence of complex fracture geometry should be considered when predicting pressure depletion induced stress under the condition of simultaneous production. This study is of great significance for understanding the evolution law of stress induced by pressure consumption, as well as the design of infill wells and repeated fracturing.

Deformation behavior of plastic zone at crack tip of 304 stainless steel

AbstractIn this paper, the deformation behavior of plastic zone at crack tip during crack propagation of 304 stainless steel was studied. Firstly, the fatigue crack propagation tests of 304 stainless steel material were carried out. Combined with digital image correlation (DIC) technology, the strain field data of plastic zone at crack tip, the size of plastic zone, and the strain evolution law of crack tip region at different positions were obtained. Then, combined with the theoretical model and the numerical model, the theoretical, numerical, and experimental results were compared to verify the accuracy of the model. Finally, the effects of cyclic loading, specimen thickness, and weld position on plastic zone at crack tip were analyzed. The research results provide a reference for accurately predicting the fatigue crack growth of 304 stainless steel.

Mechanisms of Thick-Hard Roof and Thin Aquifer Zone Floor Destruction and the Evolution Law of Water Inrush

Mechanisms of Thick-Hard Roof and Thin Aquifer Zone Floor Destruction and the Evolution Law of Water Inrush

Fatigue failure mechanism of high-speed train bearing steel after long-term service

The remaining useful life prediction of the bearings after long-term service is a crucial issue. Exploring the fatigue failure mechanism of bearing is helpful to grasp the dynamic evolution law during its whole period. In this paper, the tensile and fatigue property degradation law of bearings for high-speed train under different service mileages were explored. The results display that there is a downward trend for the tensile/fatigue strength and elongation of the bearing as prolonging its service mileage. The area proportion of dimples declines while that of cleavage plane increases in bearing steel after long-term service. The variations of secondary crack length and equivalent diameters of cleavage plane indicate that more cracks initiated in damaged bearing under cyclic loading. Based on the energy principle, the fatigue failure mechanism of bearings after long-term service was elaborated with the plastic deformation, crack initiation and growth resistance. The results would provide theoretical reference for the fatigue property degradation law of high-speed train bearings.

Study on the carbon evolution laws and carbon accounting methods of coal-fired boilers under deep peaking

The power grid needs to rapidly modify its power output to counterbalance the variable nature of new energy generation. This necessitates enhanced deep-peaking capabilities in coal-fired boilers. This research investigates the CO2 emission evolution in a 330 MW tangentially fired pulverized coal boiler under varying loads, secondary air distributions, and coal blending to elucidate the CO2 emission evolution dynamics of coal-fired boilers. A comprehensive three-dimensional model of the gas-solid two-phase flow in the furnace, coupled with the coal-firing chemical reactions, is established. Finally, a novel carbon accounting method is proposed, which enhances the accuracy of carbon emission results under various conditions. Results indicate that as the load decreases, combustion becomes progressively less complete. At a 30 % load reduction, carbon conversion into CO2 is minimized, and the amounts of CO and unburned carbon in the fly ash peak. Using equal air distribution and layering coal blending—inferior in layers A, B, and C; general in D; and superior in E—across five burners enhances CO2 production under low-load conditions. Consequently, the levels of CO and unburned carbon in the fly ash are reduced. These research findings serve as a reference for precise accounting of boiler carbon emissions.

Research on the evolution law of aggregate micro-texture during long-term wearing of asphalt pavement

The micro-texture of aggregates plays a very important role in anti-skid performance of asphalt pavement. Three kinds of aggregates, limestone, basalt and sandstone, were used to prepare asphalt mixture to study the evolution law of micro-texture of aggregates under long-term wear. One Third Scaled Model Mobile Load Simulator (MMLS) was used to conduct wear tests on asphalt mixtures under two working conditions: drying and water immersion. Laser confocal scanning microscopy (LCSM) was used to test and analyze the micro-texture evolution of the original aggregate, mixed aggregate, and accelerated loaded aggregate surfaces. The mechanism of texture evolution was analyzed based on the mineral composition, Vickers hardness, and micro-structure of the aggregate. The results show that wear can lead to varying degrees of polishing of the aggregate micro-texture, with limestone>basalt>sandstone, and the drying condition is more severe than the water immersion condition. The wear degree of the texture is closely related to the mineral composition and hardness of the aggregate. The more complex the mineral composition, the greater the hardness and hardness dispersion, the less likely the texture will be polished. The micro-texture indexes Ra and Rq have the highest correlation with anti-skid performance, and the relationship between them is obtained by fitting. Scanning electron microscopy (SEM) shows that the fine particles in sandstone are easy to fall off, which is conducive to the regeneration of texture, and is also the internal reason for its best anti-skid durability.

Enhanced strength, toughness and reliability in crab exoskeleton–inspired 3D-printed porous thermoplastics

PurposeFused deposition modeling enables multiscale structure control. However, most of this structural space is unexplored. Specifically, the impact of biomimetic porous structures on the mechanical behavior and reliability of common thermoplastics are unclear. In this work, porous structures inspired by the multifunctional crab exoskeleton were 3D-printed with different raster orientations, including fully rotating rasters similar to Bouligand structures found in biological materials. Tensile tests and simulations were performed to observe the stochastic behavior of fracture properties and to reveal the underlying origins of mechanical reliability in biomimetic porous systems.Design/methodology/approachTensile tests were performed on 3D-printed porous structures with four different rasters. These rasters were biomimetic Bouligand, semi-Bouligand, 00 raster and 45°/−45° raster. In addition, two different sets were manufactured to observe the impact of contours on the mechanical behavior. A total of 137 tensile tests were performed. A total of 88 finite element simulations were executed using Abaqus built-in Hashin damage initiation criterion and energy-based damage evolution law. Weibull analyses were performed to quantify the stochastic properties.FindingsBiomimetic Bouligand structure is effective in increasing fracture strength. Average fracture strength of the Bouligand structure was 33% higher compared to the default 45°/−45° and 10% higher compared to 00 rasters. Variations in strength were lower in Bouligand structure compared to the default 45°/−45° raster. However, 00 raster had the highest Weibull modulus m = 54 compared to Bouligand m = 25 and 45°/−45° m = 17. Simulations showed that Bouligand structure is effective in increasing the mechanical reliability through local damage accumulation around the holes. The simulated Weibull modulus of the Bouligand structure was 40 compared to the moduli of other rasters that ranged from 18 to 25.Practical implicationsThe mechanical reliability of porous Bouligand structures is higher compared to other rasters, which makes the biomimetic structure a better choice for industrial applications. Contours decrease the strength and strain at failure for 3D-printed porous structures. Bouligand structures with rotating raster orientations increased strength and strain at failure when contours are present in the porous structure.Originality/valueTo the best of the authors’ knowledge, this is the first study showing the effects of biomimetic raster orientations on the mechanical behavior and the effects of contours on the tensile fracture properties of 3D-printed porous acrylonitrile butadiene styrene using tensile tests and fracture simulations. This is the first study applying composite fracture model to anisotropic porous 3D-printed polymers.

Dynamic Evolution Law of Production Stress Field in Fractured Tight Sandstone Horizontal Wells Considering Stress Sensitivity of Multiple Media

Inter-well frac-hit has become an important challenge in the development of unconventional oil and gas resources such as fractured tight sandstone. Due to the presence of hydraulic fracturing fractures, secondary induced fractures, natural fractures, and other seepage media in real formations, the acquisition of stress fields requires the coupling effect of seepage and stress. In this process, there is also stress sensitivity, which leads to unclear dynamic evolution laws of stress fields and increases the difficulty of the staged multi-cluster fracturing of horizontal wells. The use of a multi-stage stress-sensitive horizontal well production stress field prediction model is an effective means of analyzing the influence of natural fracture parameters, main fracture parameters, and multi-stage stress sensitivity coefficients on the stress field. This article considers multi-stage stress sensitivity and, based on fractured sandstone reservoir parameters, establishes a numerical model for the dynamic evolution of the production stress field in horizontal wells with matrix self-supporting fracture-supported fracture–seepage–stress coupling. The influence of various factors on the production stress field is analyzed. The results show that under constant pressure production, for low-permeability reservoirs, multi-stage stress sensitivity has a relatively low impact on reservoir stress, and the amplitude of principal stress change in the entire fracture length direction is only within the range of 0.27%, with no significant change in stress distribution; The parameters of the main fracture have a significant impact on the stress field, with a variation amplitude of within 2.85%. The ability of stress to diffuse from the fracture tip to the surrounding areas is stronger, and the stress concentration area spreads from an elliptical distribution to a semi-circular distribution. The random natural fracture parameters have a significant impact on pore pressure. As the density and angle of the fractures increase, the pore pressure changes within the range of 3.32%, and the diffusion area of pore pressure significantly increases, making it easy to communicate with the reservoir on both sides of the fractures.

Dynamic analysis of a class of fractional‐order dry friction oscillators

This article investigates a class of Duffing nonlinear dynamic systems with fractional‐order dry friction and conducts in‐depth research on the stability, chaotic characteristics, and erosion of the safety basin of this system; the results are verified through numerical simulation. First, the average method is used to approximate the amplitude–frequency relationship of the system, and the accuracy of the analytical results is verified through numerical experiments. Second, the Melnikov method is used to obtain the conditions for the system to enter chaos in the Smale horseshoe sense, and the Melnikov curve is drawn for further verification. Then, bifurcation diagrams are drawn for the changes in various parameters in the system, with a focus on analyzing the influence of friction factors on chaotic bifurcation. By applying the definition and calculation principle of the maximum Lyapunov exponent, and drawing and utilizing the maximum Lyapunov exponent graph, the chaotic state that the system enters under different parameters is more clearly defined. Finally, the evolution law of the safety basin under various parameter changes, especially dry friction changes, is analyzed, and the erosion and bifurcation mechanism of the safety basin is studied. Comparing with the bifurcation diagram, it reveals that chaos primarily contributes to the erosion of the safety basin.