Computer Simulation Research Articles

Analysis of Liquefaction in Tailings Deposits by Fem Modeling of Undrained Cyclic Triaxial

In this article, a numerical calibration procedure for undrained cyclic triaxial tests is presented to evaluate the liquefaction potential in sand and silt samples from mining tailings in northern Chile. The numerical modeling of an axisymmetric specimen involved two stages: isotropic consolidation using the Hardening Soil Small (HSS) model and a cycling phase employing the UBC3D-PLM model to simulate the onset of liquefaction using the criterion that the excess pore pressure ratio Ru should exceed 0.8. The results demonstrate that the UBC3D-PLM modeling calibrated with experimental data from cyclic triaxial tests effectively represents the excess pore pressure in both sandy and silty soils from mining tailings. The accuracy of the modeling decreases when a single set of parameters is applied to the same soil at different cyclic stress ratios (CSR), highlighting the need for specific calibrations for each loading.

Analysis of Flood Resistance of Masonry Arch Bridges Combining Numerical Simulation and Probabilistic Analysis

ABSTRACT In the context of global climate change, extreme weather events such as heavy rainfall and floods are becoming more frequent. The safety of masonry arch bridges is facing unprecedented challenges. This study aims to explore the collapse mechanism of bridges in flood environments by combining numerical simulation and stochastic analysis. Taking the collapse of the Original Zhenhai Bridge caused by heavy rainfall in July 2020 as a specific case, a three-dimensional numerical model is established to analyze the mechanical response of the bridge under the combined action of flood scour and water flow loads. Additionally, predictions of scour depth are made using both theoretical formulas and stochastic methods, comparing the impact of stochastic analysis on the mechanical performance of arch bridges under the same flood conditions. The study reveals that compared to pure scour, the structural characteristics of the bridge undergo significant changes under the coupled effects of scour and water flow loads. In extreme flood events, stochastic analysis has a significant impact on the structural safety under high flow conditions. The research results provide case references and scientific methods for the safety analysis of masonry arch bridges in flood environments.

Stochastic Multiscale Modeling of Electrical Conductivity of Carbon Nanotube Polymer Nanocomposites: An Interpretable Machine Learning Approach

This study introduces an interpretable machine learning (ML) framework for efficiently predicting the electrical conductivity of carbon nanotube (CNT)/polymer nanocomposites. A stochastic multiscale numerical model based on representative volume element (RVE) is employed to generate a representative dataset. This dataset is used to train three ML models, including random forest, XGBoost, and artificial neural networks (ANN). The dataset includes six input features: CNT length, aspect ratio, intrinsic CNT conductivity, number of CNT conduction channels, energy barrier height, and volume fraction, with the electrical conductivity of the nanocomposites as the output feature. The findings highlight the exceptional accuracy of the ANN model in predicting electrical conductivity at significantly lower computational costs. Furthermore, the use of Shapley additive explanations (SHAP) enhances the interpretability of these ML models, identifying the volume fraction, energy barrier height, and intrinsic CNT conductivity as the most influential factors affecting conductivity. This approach sets the stage for rapid and efficient modeling of CNT/polymer nanocomposites facilitating the design of materials with tailored electrical properties for diverse applications.

Numerical Investigation of Heterogeneous Calcite Distributions in MICP Processes

Microbially induced calcite precipitation (MICP) is a sustainable and environmentally friendly technology with applications in soil stabilization, concrete crack repair, and wastewater treatment. This study presents an improved Darcy-scale numerical model to simulate the MICP processes in heterogeneous porous media. It focuses on the effects of porosity heterogeneity, characterized by average porosity and correlation length, as well as injection strategies. Both average porosity and correlation length are critical factors influencing mass transport and calcite distribution during MICP treatment. An increase in average porosity leads to significant reductions in transport distance and total calcite mass. Notably, in the case of low averaged porosity, a larger correlation length results in more heterogeneous calcite distributions. However, there exists an upper threshold value of the initial averaged porosity (ϕ0=0.45) above which the heterogeneity of the calcite does not present clear dependence on the correlation length. Additionally, injection strategies significantly impact the consolidation effects. Compared to continuous injection, using the phased injection strategy can greatly improve the precipitated calcite area and mass due to its high utility and the efficiency of reactants.

Analysis of the Effectiveness of Using a High-Accuracy Numerical Algorithm in the Simulation of Two-Phase Flow Problems

The continuous technical development of weapons and ammunition, as well as the growing requirements placed on new equipment, require the use of modern tools for their designing and testing. The analysis of the operation of powder propellant systems, which is the focus of interior ballistics, requires particularly specialised tools due to the complexity of the phenomena occurring during a shot. The development of numerical methods and the increase in the computing power of modern computers have enabled the modelling and simulation of highly complex physical problems and the development of numerical schemes with a theoretically arbitrary high order of accuracy. In this study, an attempt was made to develop a one-dimensional model of two-phase flows with a high order of solution accuracy in time, based on the finite volume method. The use of a two-phase flow model in the context of interior ballistics allows, among other things, the analysis of wave phenomena occurring during firing, which is crucial for safety reasons, especially in tank guns. The developed algorithm for solving two-phase systems was checked using recognised test problems for this type of model. The results of the work will be used to develop a numerical model to solve the main problem of interior ballistics.

Experimental characterization of acoustic damping materials

AbstractDue to their ease of use and low cost, passive damping methods are a preferred mean for the reduction of noise in many engineering applications. This applies in particular to foam materials, which exhibit good acoustic and mechanical damping properties. The selection of suitable materials is usually carried out experimentally and can be very labor‐intensive and time‐consuming. For this reason, it is helpful to develop qualified numerical methodsthat can be used for material design and selection. Hence, foam materials must be characterized experimentally in order to enable a later comparison with vibroacoustic simulations. This contribution, therefore, aims at providing suitable parameters for the use in numerical analyses and their validation. Firstly, the microstructure of a foam specimen is captured by means of a CT scan. In addition to providing the geometry for the multi‐physics simulations, this measurement is also used to determine characteristic foam features, for example, strut thickness and pore size distribution. In the second step, the frequency‐dependent stiffness and damping properties of the material are determined by a special experimental setup utilizing an electrodynamic shaker. Here, the dynamic system is approximated as a single‐mass oscillator, which is sufficiently accurate for low frequencies. These properties will later be used in the numerical model to evaluate different parameter identification approaches. In the third and last step of the experimental campaign, measurements with an impedance tube are conducted to obtain the coefficient of absorption. This material parameter is particularly suitable for comparing experiments and simulations. Finally, the correlation between the experimental results is examined to provide a deeper understanding of the foam materials.

Geological hazards when justifying the regulations of urban development

Relevance. Territorial planning of cities determines the reduction or prevention of possible or existing losses of the population, economic facilities and the environment from technoprirodnye hazards of different genesis. The effectiveness of planning regulations for urban development depends on the degree of validity of identification and forecasts of the occurrence and evolution of hazardous natural and technogenic processes in time and space. This, in its turn, is determined by the reliability of our ideas about the patterns of their formation and regional distribution. Aim. To identify geological hazards and propose the main engineering and geological constraints in urban planning on the example of Gomel. Objects. Natural and technical system "geological environment – technogenic impacts – hazardous natural and technoprirodnye processes". Methods. Systematic approach to the analysis of patterns of formation of hazardous natural and technoprirodic processes in the city; numerical modeling of geofiltration and geomigration processes based on the geofiltration model "GOMEL". Results. The most characteristic dangerous natural and technoprirodnye processes have been established on the example of the city of Gomel. Priority areas for the introduction of planning restrictions: tectonic situation; man-made flooding, were identified. Soil conditions are considered separately as one of the factors effecting the precipitation of natural and artificial bases under the influence of loads from civil and industrial buildings and structures. The authors have revealed the presence of disjunctive, as well as plicative in the form of flexures, age-varying and multi-scale dislocations forming the block structure of the upper part of the earth's crust of the city. It is shown that from the point of view of the stability of engineering structures, the absolute values of the velocities of long-period unidirectional block displacements are dangerous during long-term operation of structures located in interblock active zones. At the same time, short-period multidirectional movements in the active geodynamic zones of the articulation of blocks cause a change in the slope and bending of the foundations of structures. The paper introduces the engineering-geological features of the most common surface deposits: the terrigenous gray-colored formation of the Paleogene; glacial, periglacial and extraglacial formations of anthropogenic origin and their facies-genetic complexes. The authors identified the potential natural and technological hazards associated with the marked deposits. The paper demonstrates the causes and mechanisms of the formation of technogenic flooding of the city, as well as the genesis of the decrease in the deformation properties of moraine sandy loam with a technogenic increase in humidity. It is concluded that when establishing land use regulations for a city, as a rule, general measures cannot be taken, both for its entire territory and for certain stages of construction activity. The paper introduces the corresponding restrictions on the use of subsoil.

EXPERIMENTAL AND OPTIMIZATION-BASED ANALYSIS TO MAXIMIZE MECHANICAL-RELATED ATTRIBUTES OF FRICTION STIR WELDED CDA101 CU ALLOY JOINTS

This paper deals with the identification of optimized process parameters for FSW of 6[Formula: see text]mm thick CDA 101 Cu alloy plates. Three distinctive parameters, the tool’s traverse speed, rotational speed, axial load, have been considered as self-reliant variables for formulating a numerical model based on nonlinear regression-type approach, with the objective of anticipating mechanical attributes of joints. From the surface- and contour-type plots, it was revealed that the axial load in combination with rotational speed has played a dominant role in influencing the tensile and yield strength of the CDA101 Cu alloy joints. At the same time, the traverse speed of the tool has played a dominant role in enhancing the ductility of the CDA101 joints through the generation of ample volumes of frictional heat, which in turn have facilitated both coarsening of grain structures and softening of materials in an ideal manner. Employing tool rotational speed, traverse speed and axial load in the ranges of 1750–1900[Formula: see text]rpm, 23–27[Formula: see text]mm/min and 5.75–6.5[Formula: see text]kN, respectively, was found to enhance the tensile strength of the CDA101 Cu alloy joints. FSW process parameters attained through the numerical meta-model of RSM have led to the attainment of the highest value of tensile strength of 209.5[Formula: see text]MPa together with a 141.4[Formula: see text]MPa yield strength with very minimal fitting errors of 0.28% and 0.52%, respectively. A collection of additional experimental data has been employed to verify the formulated empirical-based formula and additionally, an entirely different optimization technique namely PSO (i.e. particle swarm optimization) was implemented for comparison and validation purposes. The formulated meta-model of RSM generated the largest value of % of elongation of 14.23%, though the relative percent error was around 1.988%, which was slightly larger than that of the value (i.e. 14.11%) generated by PSO, with a relative error of 0.58%.

Numerical investigation for the influence of suspended sediment on salinity distribution in the Qiantang Estuary, China

Numerous studies have demonstrated that high suspended sediment concentration (SSC) can change density distribution and affect water mixing, but few scholars have investigated this impact on numerical simulations of estuarian salinity distribution. The Qiantang Estuary is a macro-tidal estuary with high SSC, which has a more significant influence on water density than that of salinity. Therefore, this paper established a three-dimensional (3D) numerical model coupling flow, salinity and SSC based on Delft3D and first analyzes the impact of SSC on salinity distribution under different runoff and tidal conditions in the Qiantang Estuary. The results indicated that simulated salinity generally decreases when considering the impact of SSC, suggesting a weakening effect on saltwater intrusion. The distribution of salinity difference (ΔS) and SSC show a strong spatial and temporal correlation, and ΔS peak increases and shifts upstream as the tidal range increases or runoff discharge decreases. The mechanism of SSC influencing saltwater intrusion can be summarized as follows: On the one hand, SSC increases the water density, which weakens the driving force for saltwater to move upstream, causing a decrease in flood current velocity and water level, and thereby diminishing the advective transport of salinity. On the other hand, SSC enhances density stratification, which weakens vertical turbulence and reduces the dispersive transport of salinity. These combined effects reduce both the advective and diffusive salinity fluxes during the flood tide, ultimately leading to a decrease in upstream salinity. Therefore, neglecting this effect in estuaries with high SSC can cause significant deviations in salinity simulation results, especially under low-flow and high-tide conditions.

Partitioning the drivers of Antarctic glacier mass balance (2003–2020) using satellite observations and a regional climate model

We investigate the mass changes of Antarctic glaciers from 2003 to 2020, partitioning them into the contributions of surface mass balance (SMB) and ice discharge, using high-resolution ice mass change estimates derived from the combination of two different types of satellite observations (gravimetry and altimetry) and outputs from a regional climate model. Our analysis indicates that changes in ice discharge have played a dominant role in ongoing ice mass trends and their accelerations, especially in glaciers near the Amundsen and Bellingshausen Seas in West Antarctica. In particular, mass losses of the Thwaites and Pine Island Glaciers have been mostly (>90%) controlled by ice discharge, while the contribution of SMB has been relatively minor. In East Antarctica, SMB accounts for significant portions (>50%) of ice mass imbalances of glaciers in e.g., Dronning Maud Land and Wilkes Land. Ice discharge has also played a notable role in overall mass gain in the region. While our ice discharge estimates agree well with previous estimates from satellite imagery in West Antarctica, notable differences are found in glaciers of East Antarctica and the Antarctic Peninsula. This highlights the need for more observations and improved numerical models to refine these estimates.

Nd Isotopic Equilibration During Channelized Melt Transport Through the Lithosphere: A Feasibility Study Using Idealized Numerical Models

AbstractThis study is motivated by the observed variability in trace element isotopic and chemical compositions of primitive (Si52 wt %) basalts in southwest North America (SWNA) during the Cenozoic transition from subduction to extension. Specifically, we focus on processes that may explain the enigmatic observation that in some localities, basalts with low Ta/Th, consistent with parental melts in a subduction setting, have signatures consistent with continental lithospheric mantle (CLM). In locations with the oldest CLM (Proterozoic and Archean), Cenozoic basalts with low Ta/Th have well below zero. We model channelized magma transport through the CLM using simple 1D transport models to explore the extent to which diffusive and reactive mass exchange can modify Nd isotopic compositions via open system melt‐wallrock interactions. For geologically reasonable channel spacings and volume fractions, we quantify the reactive assimilation rates required for incoming melt with a different than the wall‐rock to undergo a substantial isotopic shift during 10 km channelized melt transport. In the presence of grain boundaries, enhanced diffusion between melt‐rich channels and melt‐poor surrounding rock contributes to isotopic equilibration, however this effect is not enough to explain observations; our models suggest a significant contribution from reactive assimilation of wall‐rock. Additionally our models support the idea that the observed covariability in Ta/Th and in Cenozoic basalts cannot be attributed to transport alone and must also reflect the transition from subduction‐related to extension‐related parental melts in SWNA.

A review of multiscale numerical modeling of rock mechanics and rock engineering

AbstractRock is geometrically and mechanically multiscale in nature, and the traditional phenomenological laws at the macroscale cannot render a quantitative relationship between microscopic damage of rocks and overall rock structural degradation. This may lead to problems in the evaluation of rock structure stability and safe life. Multiscale numerical modeling is regarded as an effective way to gain insight into factors affecting rock properties from a cross‐scale view. This study compiles the history of theoretical developments and numerical techniques related to rock multiscale issues according to different modeling architectures, that is, the homogenization theory, the hierarchical approach, and the concurrent approach. For these approaches, their benefits, drawbacks, and application scope are underlined. Despite the considerable attempts that have been made, some key issues still result in multiple challenges. Therefore, this study points out the perspectives of rock multiscale issues so as to provide a research direction for the future. The review results show that, in addition to numerical techniques, for example, high‐performance computing, more attention should be paid to the development of an advanced constitutive model with consideration of fine geometrical descriptions of rock to facilitate solutions to multiscale problems in rock mechanics and rock engineering.

ANALYSIS OF THE CONSTRUCTION OF THE CAR TRAILER FRAME IN TERMS OF CHANGING THE ASSEMBLY TECHNOLOGY

Car trailers are quite a popular means of transport and are offered in many versions, from single axle light trailers with a maximum permissible weight of 750 kg, through two- or more-axle specialized trailers. The issues of research car trailers focus on two directions: testing the driving properties and analysing the strength of the supporting frame system. Issues related to the construction of light trailers are often common to trailers used in agriculture or general transport. In this article, based on a mass-produced car trailer, an analysis was carried out regarding the choice in the technology of making the supporting structure consisting of a lower frame and an upper frame. The term upper frame should be understood as the structure on which the lifted load box rests. The costs of materials, assembly and technological possibilities of small-scale production were taken into account. In addition, strength analyses of the numerical models were carried out for critical areas of the frame during operation. After considering the unit costs for each of the analysed assembly technologies, it was shown that riveting would be the cheapest. However, the most suitable method of assembly is welding, as it allows the use of standard profiles.

Application of artificial intelligence technologies at substantiation of effective control algorithms for the electrical engineering complex of urban electric transport traction electrical equipment

Background. Studies on the current state and development trends of urban ground public electric transportation in Russia highlight the urgent need for innovative technologies. These innovations should focus on designing and operating new types of rolling stock, traction electrical equipment, and promising types of electric traction. At the same time, special attention should be paid to building autonomous control systems for electric transportation using artificial intelligence technologies. Aim. This study aims to explore the application of neural networks to develop algorithms for effectively controlling the electrical engineering complex of traction electrical equipment in urban ground rail electric transport. Materials and Methods. The research utilized data from studies on the traction electrical equipment modes of rolling stock. These studies were conducted through both computer simulations and field experiments under real operating conditions of urban electric transport using different control algorithms. By applying probability theory and mathematical statistics, the analysis identified the correlations between operational and energy parameters of rolling stock movement and the operational mode of its traction electrical equipment. Results. The identified correlation dependencies informed the design of an effective network architecture, including its size and complexity, as well as the composition of their training samples. This led to the development of an original, simplified algorithm for determining effective control parameters for the electrical complex of traction electrical equipment during the movement of a vehicle on a given section of track. Conclusion. The research concluded that using “simple” neural networks for calculating the parameters of effective control of traction electrical equipment operation modes in urban electric transport provides higher speed and sufficient accuracy compared to complex neural network models. These results are valuable for developers of intelligent control systems for streetcar transportation.

Influence of the quality of electrical energy on the operation of the main elements of electric arc plasma systems for thermochemical fuel preparation

The work is devoted to the study of the influence of the quality of electrical energy on the service life of electrodes of plasma thermochemical fuel treatment systems used in thermal power plants running on coal fuel. The issues of moving the near-electrode plasma section using magnetic wave scanning are discussed in detail. The ANSYS Maxwell software was used as the main research tool. The computer simulation was performed with the following main experimental parameters: the coefficient of asymmetry in the zero sequence (K0U) is 2% and 4%; the coefficient of asymmetry in the negative sequence (K2U) is 2% and 4%; the voltage deviation dU(+) is 5% and 10%; the materials used for manufacturing electrodes: copper (Cu), a pseudo-alloy of tungsten, nickel and copper – (W + Ni + Cu), a pseudo-alloy of molybdenum, tungsten and copper – (Mo + W+ Cu). Based on the obtained simulation results, graphical representations of changes in the trajectories of the near-electrode plasma section under various distorting factors, dependences of changes in the values of specific erosion of various electrode materials, graphical dependences of changes in the service life of electrodes of plasma systems are constructed. The studies performed using computer modeling based on the ANSYS Maxwell software product made it possible to quantify the effect of voltage distortions on the trajectory of the near-electrode plasma section and, consequently, the service life of the plasmatron electrodes. In particular, in the course of the study, a detailed assessment and analysis of the degree of influence of the following indicators of electrical energy quality was performed: the coefficients of asymmetry in the negative and zero sequence, voltage deviations on the operation and technical condition of the electrodes of plasma systems. The results of the study were discussed, and recommendations were formulated for the use of plasma thermochemical fuel preparation systems used at thermal power plants to ignite the fuel of pulverized coal boilers.

Robustness and Scalability of Incomplete Virtual Pheromone Maps for Stigmergic Collective Exploration

The Swarm Guiding and Communication System (SGCS) is a decision-making and information-sharing framework for robot swarms that only needs close-range peer-to-peer communication and no centralized control. Each robot makes decisions based on an incomplete virtual pheromone map that is updated on each interaction with another robot, imitating ant colonial behavior. Similar systems rely on continuous communication with no range limitations, environment modification, or centralized control. A computer simulation is developed to assess the effectiveness and robustness of the framework in covering an area. Consistency and the time needed for 99% coverage are compared with an unbiased random walk. The pheromone approach is shown to outperfom the unbiased one regardless of number of agents. Innate resilience to individual failures is also demonstrated.

Multi-Scale Integrated Corrosion-Adjusted Seismic Fragility Framework for Critical Infrastructure Resilience

This study presents a novel framework for integrating corrosion effects into critical infrastructure seismic risk assessment, focusing on reinforced concrete (RC) structures. Unlike traditional seismic fragility curves, which often overlook time-dependent degradation such as corrosion, this methodology introduces an approach incorporating corrosion-induced degradation into seismic fragility curves. This framework combines time-dependent corrosion simulation with numerical modeling, using the finite–discrete element method (FDEM) to assess the reduction in structural capacity. These results are used to adjust the seismic fragility curves, capturing the increased vulnerability due to corrosion. A key novelty of this work is the development of a comprehensive risk assessment that merges the corrosion-adjusted fragility curves with seismic hazard data to estimate long-term seismic risk, introducing a cumulative risk ratio to quantify the total risk over the structure’s lifecycle. This framework is demonstrated through a case study of a one-story RC moment frame building, evaluating its seismic risk under various corrosion scenarios and locations. The simulation results showed a good fit, with a 3% to 14% difference between the case study and simulations up to 75 years. This fitness highlights the model’s accuracy in predicting structural degradation due to corrosion. Furthermore, the findings reveal a significant increase in seismic risk, particularly in moderate and intensive corrosion environments, by 59% and 100%, respectively. These insights emphasize the critical importance of incorporating corrosion effects into seismic risk assessments, offering a more accurate and effective strategy to enhance infrastructure resilience throughout its lifecycle.

Selecting Erosion- and Deposition-Dominated Zones in the Jezero Delta Using a Water Flow Model for Targeting Future In Situ Mars Surface Missions

Identifying surface sites with significant astrobiological potential on Mars requires a comprehensive understanding of past geological processes and conditions there, including the shallow subsurface region. Numerical modelling could distinguish between regions dominated by erosion and those characterized by sediment accumulation in ancient wet environments. The target area of Jezero Crater is relatively well explored and thus is an ideal site to evaluate model calculations; however, important works are still missing on expectations related to its shallow subsurface . In this work, the best available approaches were followed, and only surface morphology was considered (supposedly formed by the last fluvial episode). The shallow subsurface became an important target recently, and this model could provide new inputs in this area. Erosion–accumulation models are suitable for terrestrial surface features, but few have been applied to Mars. This work addresses this challenge using the SIMWE (SIMulated Water Erosion) model on the Jezero Crater delta, the landing site of the Perseverance rover. For calculations, the average grain size according to the THEMIS TI data was applied to the target area. The flow depth varied between 1.89 and 34.74 m (average of 12.66 m). The water-filled channel width ranged from 35.3 to 341.42 m. A flow velocity of 0.008–11.6 m/s, a maximum erosion rate of 5.98 g/m2/h, and a deposition 4.07 g/m2/h were estimated. These calculated values are close to the range of estimations from other authors assuming precipitation of 1–20 mm/h and discharges of 60–400 m3/s. The model was able to distinguish between erosion- and accumulation-dominated areas about 1 m above Jezero Crater’s delta that are not visible from above. This model helps to identify the accumulation-dominated areas with the finest grain size with good preservation capability for the shallow but invisible subsurface.

Quantification of ocean microplastic fragmentation processes in the Sea of Japan using a combination of field observations and numerical particle tracking model experiments

This study estimated the fragmentation rate of microplastics (MiPs) in the Sea of Japan by analyzing MiP size over time following their generation from macroplastics (MaPs). A 5-year particle-tracking model was used to simulate the MaP and MiP motions driven by ocean currents, Stokes drift, the windage of MaPs, beaching, re-drifting, the conversion process from MaPs to MiPs, and the removal of MiPs from the upper ocean. MiP sizes decreased downstream in the Tsushima Current flowing northeastward. The highest correlation between MiP size and elapsed time occurred in the simulation where MiP fragmentation also occurred in the ocean, at 20 % of the rate on beaches. The apparent fragmentation rate in nature was estimated to approximately 1.0 mm/100 days. This study demonstrated that incorporating spatiotemporal information from the simulation into observed size results could further our understanding of fragmentation of MiPs degraded in marine environments.

Design of Roadway Support for Soft Coal Seam Back Mining in Wangxingzhuang Coal Mine

Wangxingzhuang coal mine has soft coal seams, difficult mining, and complex mining disturbed stress field, resulting in serious deformation and damage and disrepair of the perimeter rock of the back-mining roadway. On the basis of geomechanical investigation, combined with relevant engineering cases, a feasible perimeter rock control scheme and support parameters are proposed for the deformation and damage characteristics of the back-mining roadway in Wangxingzhuang coal mine. A numerical simulation model is established to verify the support parameters. The results prove that the support scheme proposed in this paper can meet the needs of perimeter rock control in the back-mining roadway of Wangxingzhuang coal mine, and has certain theoretical significance and engineering practice value.