Numerical Software Research Articles

A Study on the Impact of Different Delay Times on Rock Mass Throwing and Movement Characteristics Based on the FEM–SPH Method

Burst morphology is a crucial indicator for evaluating the effectiveness of blasting, as it directly reflects the actual state of the blasting results. The results of rock displacement following blasting partially reflect the effectiveness of throw blasting, while the rock ejection process serves as the macroscopic manifestation of the blasting method. To accurately assess the impact of different delay times on burst formation, this study addressed the issues of rock movement and ejection in underground blasting. Using three-dimensional modeling, we constructed a FEM–SPH model and utilized LS-DYNA numerical simulation software to investigate the movement patterns of rock in precise delayed blasting scenarios underground. This study explored the spatiotemporal evolution characteristics of rock movement post-blasting. Digital electronic detonators were used to set precise inter-row delay times of 25 ms, 50 ms, and 75 ms. The results revealed that the ejection distance of blasted rock in underground mining increased with longer inter-row delay times, while the slope angle of the blasted muck pile decreased as the delay time increased. Furthermore, at a micro level, the study found that a 75 ms delay created new free surfaces, providing effective compensation space for subsequent blasts, thereby improving blasting outcomes. Analysis of the 25 ms and 50 ms delay periods indicated a clamping effect on rock movement. Field comparisons of blasting results were conducted to validate the influence of precise delay times on the movement patterns and spatiotemporal evolution characteristics of blasted rock.

Simulations on Evacuation Strategy and Evacuation Process of the Subway Train Under the Fire

This study focuses on the safe evacuation strategy and evacuation process in the subway train under the fires. The subway station evacuation mode should be adopted if the power system of a subway train is normal on fire. While, the tunnel evacuation mode should be adopted if the power system of the train fails because of the effects of fire. Under the tunnel evacuation mode, the direction of tunnel smoke should be opposite to that of most passengers, and passengers should be evacuated toward the fresh wind. By using the numerical simulation software Pathfinder and PyroSim, the passenger evacuation time under different conditions is calculated, and the safety of the evacuation process is evaluated. The results show that the evacuation time of the station evacuation mode is obviously shorter than that of the tunnel evacuation mode. With the same conditions, the evacuation time of the tunnel evacuation mode is 2193 s, which is about four times as much as the evacuation time of the station evacuation mode (526 s). The total evacuation time increases with the total number of passengers and the proportion of older people and children. Under an oil pool fire, which is an extreme fire condition, the fire environment inside the train may reach a level threatening the passengers’ safety before the evacuation is complete, even before the door opens; therefore, special attention should be paid to the safety issues in stage from the fire begins to the evacuation complete.

An analytical I-V model of SiC double-gate junctionless MOSFETs

Silicon carbide(SiC) double gate junctionless metal oxide semiconductor field-effect transistors(DG JL MOSFETs) have attracted significant attention due to their ideal high temperature characteristics and radiation resistance. Therefore, it is meaningful to exploit an I-V model for SiC DG JL MOSFETs. In this article, we make a linear approximation to describe the relationship between the surface mobile charge density and the surface electron concentration of the device. Based on this approximation and using the one-dimensional Poisson’s equation, we solve for the potential distribution of a SiC DG JL MOSFET in the subthreshold region. From this solution, we derived a functional relationship between the surface mobile charge density in the channel and the channel quasi-Fermi potential. Then we successfully developed a unified I-V model for the SiC DG JL MOSFETs. Based on the drain to source current calculation formula, the calculation expressions for the device’s transconductance and output conductance are derived. By comparing our model with the results from the two-dimensional numerical simulation software Silvaco Atlas, our model’s calculations closely match the two-dimensional numerical simulation results from the subthreshold region to the accumulation region. This model has reference significance for SiC DG JL MOSFETs in the high temperature electronic circuit application field.

Acoustic Design of a Modular Studio Box for Use as Control and Mixing Room for Immersive Audio

This article presents a detailed analysis of the acoustic design process of a soundproofed box intended for control room and mixing for immersive audio. The main objectives of the project were defined based on the quality criteria required for critical listening environments, particularly for Dolby Atmos rooms. The adopted methodological approach relies on predictive evaluations of key acoustic parameters, supported by studies and calculations conducted using advanced numerical models and software. Key phases of the design process include the optimisation of the box dimensions, the modal analysis for low-frequency response, the design of the listening layout and early reflections, the specification of the internal acoustic treatment. The technical solutions adopted aim to ensure an accessible and cost-effective environment, achieved through optimization in material and space selection while still meeting requirements.

Research of Heat Transfer Characteristics of Supercritical CO2 in Spiral Tube Heat Exchangers

The spiral tube heat exchanger has the advantages of stable structure, high space utilization rate, and strong heat transfer ability. Supercritical carbon dioxide (s-CO2) has the advantages of high thermal conductivity, high specific heat capacity, and low viscosity. When the spiral tube heat exchanger uses supercritical carbon dioxide as the heat transfer medium, the two complement each other’s advantages and can further improve the heat transfer efficiency and optimize the space utilization rate. However, under the coupling influence of factors such as centrifugal force and buoyancy force, the heat transfer process of s-CO2 in the spiral tube is extremely complex. At present, there is no design criterion or empirical correlation with strong universality in the industrial field. Based on this, this paper mainly uses ANSYS FLUENT, a numerical simulation software suitable for numerical simulation, to study the influence of internal vortex field and secondary flow intensity on the thickness of thermal boundary layers through a combination of numerical simulation and theoretical analysis. At the same time, by changing various operating and structural parameters, the heat transfer characteristics of supercritical s-CO2 in the spiral tube heat exchanger are analyzed and summarized. In addition, we also established a new buoyancy factor, Fu, and tested that when Fu < 1.6 × 10−5, the effect of the buoyancy lift in the pipeline disappeared. At the same time, a correlation formula for predicting the heat transfer performance of spiral tube heat exchangers with a Nu number error less than 20% is established in this paper, which can provide solid theoretical guidance for industrial application.

MLTL Multi-type: A Typed Logic for Cyber-Physical Systems

Modern cyber-physical systems-of-systems (CPSoS) operate in complex systems-of-systems that must seamlessly work together to control safety- or mission-critical functions. Linear Temporal Logic (LTL) and Mission-time Linear Temporal logic (MLTL) intuitively express CPSoS requirements for automated system verification and validation. However, both LTL and MLTL presume that all signals populating the variables in a formula are sampled over the same rate and type (e.g., time or distance), and agree on a standard “time” step. Formal verification of cyber-physical systems-of-systems needs validate-able requirements expressed over (sub-)system signals of different types, such as signals sampled at different timescales, distances, or levels of abstraction, expressed in the same formula. Previous works developed more expressive logics to account for types (e.g., timescales) by sacrificing the intuitive simplicity of LTL. However, a legible direct one-to-one correspondence between a verbal and formal specification will ease validation, reduce bugs, increase productivity, and linearize the workflow from a project’s conception to actualization. Validation includes both transparency for human interpretation, and tractability for automated reasoning, as CPSoS often run on resource-limited embedded systems. To address these challenges, we introduced Mission-time Linear Temporal Logic Multi-type (Hariharan et al., Numerical Software Verification Workshop, 2022), a logic building on MLTL. MLTLM enables writing formal requirements over finite input signals (e.g., sensor signals and local computations) of different types, while maintaining the same simplicity as LTL and MLTL. Furthermore, MLTLM maintains a direct correspondence between a verbal requirement and its corresponding formal specification. Additionally, reasoning a formal specification in the intended type (e.g., hourly for an hourly rate, and per second for a seconds rate) will use significantly less memory in resource-constrained hardware. This article extends the previous work with (1) many illustrated examples on types (e.g., time and space) expressed in the same specification, (2) proofs omitted for space in the workshop version, (3) proofs of succinctness of MLTLM compared to MLTL, and (4) a minimal translation to MLTL of optimal length.

Existence of traveling wave solutions in continuous optimal velocity models

In traffic flow theory, hydrodynamic models, a subset of macroscopic models, can be derived from microscopic-level car-following models. Self-organized wave propagation, which characterizes congestion, has been replicated in these macroscopic models. However, the existence of wave propagation has only been validated using numerical technique or formal analyses and has not yet been rigorously proven. Therefore, analytical approaches are necessary to ensure their validity rigorously. This study investigates the properties of solutions corresponding to congestion with sparse and dense waves. Specifically, we demonstrate the existence of traveling back/front, traveling pulse, and periodic traveling wave solutions in macroscopic models. All theorems are proven using phase-plane analysis without local bifurcation theory. The key to the proofs is the monotonicity of solution trajectories concerning implicit parameters that naturally appear in the models. We also examine the global bifurcation structure for heteroclinic, homoclinic, and periodic orbits, which correspond to traveling back/front, traveling pulse, and periodic traveling wave solutions, via the numerical continuation package HomCont/AUTO.

Experimental study on the effect of dual fire sources elevation on flame merging characteristics in road tunnels with natural ventilation

A series of experiments was carried out using a 1:10 small-scale model tunnel to investigate the effect of different fire source elevations on the flame merging characteristics of dual fire sources under natural ventilation. CFD numerical simulation software was used to analyze the impact of fire source elevation on flame merging characteristics and air entrainment from the perspectives of lateral velocity, longitudinal velocity, vertical velocity and air entrainment rates. The effects of fire source spacing and heat release rate on the probability of flame merging and the height of the flame merging point under different fire source elevations were investigated. It is found that as the fire source spacing S increases, the flame tilt angle θ first increases and then stabilizes for the first time. With a further increase in S, θ increases again and stabilizes for the second time. The first stabilization is mainly controlled by the interaction between the fire sources, while the second stabilization is primarily influenced by natural ventilation from both ends of the tunnel. As the fire source elevation increases, the flame merging probability decreases with the addition of air entrainment between the fire sources. What is more, when the fire source elevation exceeds 8 cm, the effect of the fire source elevation on the flame merging probability is no longer obvious. The prediction model of flame merging probability under different fire source elevations is proposed. Meanwhile, the dimensionless fire source elevation h/D and the modified dimensionless heat release rate Q̇DS∗ are used to express the dimensionless height of flame merging point Zm/Hef. The piecewise function of Zm/Hef and Q̇DS∗ is obtained. WhenQ̇DS∗≤88,h/D≤8, the relationship between Zm/Hef and Q̇DS∗ could be described with exponential function and whenhD>0.8, the fire source elevation has little impact on Zm/Hef. WhenQ̇DS∗>8, Zm/Hef is approximately constant and increases with fire source elevation.

Design of the electrode structure for the dust removal efficiency optimization of electrostatic precipitator

Abstract Air dust pollution prevention and control is a global environmental concern. Electrostatic precipitators (ESPs) play a crucial role in dust pollution control, and its electrode structure significantly affects plasma parameters and the movement of charged dust particles, which in turn influences the dust removal efficiency. This study focuses on optimizing the commonly used wire-plate dust removal device, experimental research is conducted first, revealing that the hole-hole electrode structure of an ESP exhibits higher peak currents at the same voltage compared to both plate-hole and plate-plate electrode configurations. For the removal efficiency of particles with a radius greater than 1 μm, the ESP with a hole-hole electrode structure performs better than those with plate-hole and plate-plate electrode structures. Moreover, the particle collection efficiency is positively correlated with particle radius, larger particles correspond to higher dust removal efficiency. Based on experimental findings, simulations are conducted to provide a deeper analysis and understanding of the dust removal mechanisms of ESPs. This paper primarily utilizes COMSOL Multiphysics numerical simulation software to simulate traditional wire-plate ESPs and a new wire-hole ESP. It explores the multiphysical characteristics of ESPs under different electrode structure configurations, and investigates in depth the impact of electrode structure on dust removal efficiency. The simulation results indicate that compared to flat collecting electrodes, porous collecting electrodes exhibit a significant increase in electric field intensity at the openings, and they have a lower surface charge density than flat plates, thereby reducing reverse corona phenomena. Ion wind affects the internal airflow dynamics of ESPs, thereby influencing their efficiency in capturing fine particulate matter. This negative impact of ion wind can be mitigated by introducing perforations in the collection plates.

Numerical Simulation on the Influence of the Distribution Characteristics of Cracks and Solution Cavities on the Wellbore Stability in Carbonate Formation

The development of cracks and solution cavities in carbonate reservoirs can notably reduce the rock’s mechanical properties, leading to a severe wellbore collapse problem during drilling operations. To clarify the influence of the characteristics of cracks and solution cavities on the wellbore stability in the Dengying Formation carbonate reservoir in the Gaoshiti–Moxi area of Sichuan, the mechanical properties of carbonate rock were analyzed. Then, the influences of the attitude and width of cracks, the size and quantity of solution cavities, and their connectivity on wellbore stability were studied using FLAC3D 6.00 numerical simulation software. Our results show the following: (1) The cracks and solution cavities in the Dengying Formation carbonate rock cause significant differences in the rock’s mechanical properties. (2) The equivalent drilling fluid density of collapse pressure (ρc) considering the effects of cracks and solution cavities is 6.4% higher than without these effects, which is in good accordance with engineering practice. Additionally, cracks play a more significant role than solution cavities in affecting the wellbore stability. (3) When the orientation of a crack is closer to the direction of maximum horizontal stress, and the dip angle and width of the crack increase, the stress and deformation at the intersection of the crack and wellbore gradually increase, and correspondingly, ρc also increases. (4) The stress and displacement of various points around the solution cavities gradually increase with the increases in diameter and quantity of solution cavities, and ρc also increases. (5) Compared with the situation where cracks and solution cavities are not interconnected, the stress disturbance area around the wellbore is larger, and ρc is greater when cracks and solution cavities are interconnected.

Numerical simulation study on the force of overwintering foundation support structure of unsaturated seasonal permafrost under indoor experiments.

When analysing the effect of negative temperature on overwintering pit constructions of unsaturated soil, using the mechanical parameter of saturated soil at room temperature leads to an inaccuracy in the research findings. The strength parameters are obtained through indoor experiments. The foundation pit model is created using FLAC3D numerical simulation software based on the indoor experimental data. The influence of different parameters on the stress and deformation of the overwintering deep foundation pit supporting structure is analysed. The numerical simulation results obtained are compared with the actual monitoring data. According to research, the matric suction of the silty clay in its natural state in the Changchun area is 70 kPa. As the temperature decreases, the total cohesion of the unsaturated soil increases, and the internal friction angle tends to decrease. The numerical simulation results are consistent with the actual monitoring data changes. With the excavation, the horizontal displacement of the supporting structure increases first and then decreases, reaching the maximum displacement at two-thirds of the foundation pit. Compared with room temperature, the deformation of the supporting structure is larger under a negative temperature condition. The deformation of the supporting structure simulated by the actual temperature mechanical parameters is larger than that under the condition of normal temperature mechanical parameters. The frost-heaving force increases with the overall excavation, and a surge occurs at the bottom of the pit. The frost-heaving force changes most significantly under the condition of freezing at -20°C for 30 days.

Pavement performance analysis considering the influence of suction and temperature on subgrade soil modulus

Subgrade soils can experience varying suction (equivalent to moisture) and temperature, which could affect the modulus of subgrade soil and the response of flexible pavements. The influence of these two factors on pavement performance has not been explored on purpose despite their importance in a changing climate. In this study, a simple approach was proposed for analyzing pavement performance with consideration of suction and temperature effects on soil modulus. Modulus values at different conditions of suction and temperature are calculated using a semi-empirical equation. Then, they are utilized in the numerical software KENLAYER to evaluate the stress and strain distributions in pavements composed of an asphalt layer, base layer, and subgrade layer based on linear elastic theory. The computed results are used to predict the rutting and fatigue cracking using some field-calibrated semi-empirical equations. In parametric studies using this method, the suction and temperature-dependent modulus of a clayey soil determined from laboratory experiments was applied. The results reveal that an increase in subgrade soil modulus due to drying and cooling can reduce the total rut depth mainly due to a reduction of the vertical compressive strain in the subgrade layer. For instance, a change in suction from 0 to 150 kPa (equivalent to optimum moisture content of the clayey soil) can lead to a 34 % decrease in the total rut depth. However, the effects of an increase in soil subgrade modulus on fatigue cracking are influenced by the asphalt thickness. A reduced subgrade modulus because of wetting and heating decreases horizontal tensile strain (more susceptible to fatigue cracking) at the lowest section of relatively thinner asphalt but increases horizontal tensile strain for thicker asphalt.

Examining the Effectiveness of Continuous Integration and Continuous Deployment (CI/CD) Within SDLC Frameworks

Continuous Integration and Continuous Deployment frameworks have been a subject of significant discussion among numerous software development cycles however their impact on Software Development Life Cycles has not been investigated in detail in this research. With survey data collected from the subjects and interviews conducted on the pre-researched parameters of software quality, time-to-market of the product, and interteam collaboration, the study compares the firms that use CI/CD with those that use more conventional software development approaches. The study shows that CI/CD integration enhances software quality since defects are reduced from 1.5 to 0.5 per thousand lines of code. Furthermore, CI/CD teams have 10 releases per month against 3 of traditional and time-to-market that reduces from 6 weeks to 2 weeks. The improvement of collaboration is very visible, the CI/CD teams said that they had more meetings and their job satisfaction scores increased. As with any change, transitioning to CI/CD has potential drawbacks such as resistance to change, but these can be managed especially through training and what has been said in the article specific practices. Therefore, this research offers important insights to organisations aiming at improving their software development processes and preparedness in a relatively dynamic market environment

Bio-inspired Design and Performance Evaluation of Lawnmower Cutter Blade

This study presents a bio-inspired design approach for the engineering applications of agricultural implements, with practical implications for the field. The CAD modeling of bio-inspired serrated lawnmower cutter blades was followed by structural and fatigue analysis with a comparative evaluation of the performance of conventional plain cutter blades. This method involved the bionic transformation of a biological model of the serrated leg of the grasshopper, traced using MATLAB edge detection and image processing software. The solid modeling CAD software for the Bionic Blade lawnmower cutter blade converts this biological model into an engineering model. After designing the geometric shape of the bionic tooth profile's serrated cutting edge, the boundary conditions are evaluated for numerical simulation based on the machine specifications of the conventional lawnmower cutter. After analyzing the power of the motor, the rotational speed, and the torque, the static structural and fatigue analysis was carried out in numerical simulation software to synthesize the performance of the bionic lawnmower cutter blade for engineering analysis of deformation, shear stress, and fatigue life. The bionic design of the blade and conventional plain blades were manufactured using an Amada CO2 laser cutting machine and sheet metal bending processes for Stainless Steel (SS 304) material. The experimental validation of software simulation was carried out by field trials using the Bionic design of a lawnmower cutter blade to evaluate the performance improvement in cutting time and yield of grass cut in kg for 60 to 100 meters distance cutting. After comparing the performance of the plain blade and the design of the bionic blade, there was a reduction in abrasive wear, a 12.87% improvement in the cutting time, and a 14.8% increase in the yield of grass compared to the conventional blade. This performance improvement was attributed to the unique bionic design of serrated cutting-edge and ground clearance. The bio-inspired design approach has immense practical applications in designing agricultural implements, as confirmed by this performance evaluation of this case study of bio-inspired lawnmower cutter blades.

Introduction to the Investigation of reproduction of the real geometry of UBM panels

The article concisely describes UBM technology (Ultimate Building Machine), that can be used as a solution for buildings and roof structures. The structure of this technology is made of double corrugated thin-walled steel profiles manufactured on the construction site by a self-contained UBM manufacturing factory on wheels. These panels serve as both the building envelope and the structural system. The specificity of the construction poses many design problems, especially the determination of the strength parameters and stiffness of the double-corrugated panels from which the structure is made. The article presents the results of spatial scanning tests of double-corrugated steel sheets, which were carried out using commonly available 3D scanning devices: Leica 3D Disto and MagiScan app. Additionally, results of numerical analyses performed on scanned samples and a comparison of these results with preliminary laboratory tests have been presented in the article. The purpose of scanning was to obtain an accurate and real geometry of the UBM panels, to implement it into a numerical software, and then to perform numerical analyses. Commonly available 3D scanning devices were used because using advanced 3D scanners is not popular nowadays for economic reasons, and hand-built geometric models pose a lot of problems and are not accurate enough. Promising results were obtained, which form the basis for further research.

γ-Cascade V4: A semi-analytical code for modeling cosmological gamma-ray propagation

Since the universe is not transparent to gamma rays with energies above around one hundred GeV, it is necessary to account for the interaction of high-energy photons with intergalactic radiation fields in order to model gamma-ray propagation. Here, we present a public numerical software for the modeling of gamma-ray observables. This code computes the effects on gamma-ray spectra from the development of electromagnetic cascades and cosmological redshifting. The code introduced here is based on the original γ-Cascade, and builds on it by improving its performance at high redshifts, introducing new propagation modules, and adding many more extragalactic radiation field models, which enables the ability to estimate the uncertainties inherent to EBL modeling. We compare the results of this new code to existing Monte Carlo electromagnetic transport models, finding good agreement within EBL uncertainties.

Design of gas control lane of 9# coal seam in Wuhushan Mine based on layer layout optimization

In order to address the issue of gas over limit in the upper corner of the working face of the 9# coal seam in Wuhushan Mine, a series of theoretical and numerical simulation analyses were conducted to evaluate the optimal configuration for the gas control lane of the 9# coal seam. In accordance with the "O" circle theory and the lithology of the overlying rock strata of the 9# coal seam, the height range of the fallout zone and fissure zone in the working face mining area was determined by employing empirical formulas. The change rule and distribution characteristics of the porosity of the fissure zone and the fall zone in the mining area were analyzed based on the characteristics of rock movement and fall. The determination method was also provided. The numerical simulation software was employed to simulate and analyze the gas concentration field in the air-mining zone under conditions of no extraction and six distinct layer positions of the gas control lane. The optimal layer position of the gas control lane in the 9# coal seam was determined and subsequently implemented in the field. The results demonstrate that the overlying rock layer in the 9# coal seam exhibits a height range of 6.86 ~ 11.26 m, while the fissure zone displays a height range of 30.11 ~ 41.31 m. When the gas control road is situated in close proximity to the working face, the gas concentration field exhibits a markedly low concentration. When the distance between the gas control lane and the return airway of the working face is 20 m and the distance from the top of the coal seam is 20 m, the gas concentration in the upper corner and the return airway is 0.35% and 0.26%, respectively. These values are close to the lowest concentration observed in the layout scheme. Additionally, the gas extraction concentration and the pure volume of the gas control lane are 23.7% and 38.3 m3 min−1, respectively. These values represent the highest concentrations observed in the various layout schemes. The application of the gas management lane in the field, based on the numerical simulation results, demonstrated a successful extraction effect, which was consistent with the numerical simulation results. This effectively managed the issue of an over-limit of gas in the upper corner of the working face of the 9# coal seam.

Static output feedback control for positive 2D discrete systems in Roesser model

Abstract This paper investigates the output stabilization problem for multi-input multi-output positive 2D systems described by a linear discrete-time Roesser system. This kind of systems have the property that the horizontal and vertical states take non-negative values whenever the initial boundaries are non-negative. Conditions for the existence of desired static output feedback controllers guaranteeing the resultant closed-loop system to be positive and asymptotically stable are obtained. Also, the synthesis of static output feedback controllers, including the requirement of positivity of the controllers, is solved. These results are then extended to the case of uncertain plants. All of the proposed conditions are expressed in terms of Linear Programming, which can be easily verified by using standard numerical software. Finally, several numerical examples are included to illustrate the validity and the effectiveness of the developed results.

Influence of reservoir water saturation on coalbed methane development: Based on multi-field coupling numerical simulation

Water in coal reservoirs serves as both the driving force of reservoir fluid and the resistance to coalbed methane (CBM) migration. However, the impact and control mechanism of water saturation on CBM development in the field have not been revealed yet. The evolution of reservoir dynamics and gas production rate under various initial water saturation conditions is simulated based on the COMSOL multiphysics numerical simulation software. The simulation results show that the reduction of initial water saturation can significantly enhance the effective permeability of the reservoir gas phase, which can be increased by up to 4764%. This increases the rate of gas production, which is the main mechanism for controlling production by water saturation in the reservoir. For high water saturation reservoirs, the positive effect of inter-well interference should be reasonably utilized, and high well network density should be used to rapidly increase gas-effective permeability. For low water saturation reservoirs, the negative effect of inter-well interference should be avoided, and low well network density should be used to increase the control of single well reserves.

Failure characteristics and pressure relief effectiveness of non-persistent jointed rock mass with holes

In tunnel engineering, the rock mass contains a significant number of irregularly distributed joints, and typically exhibits high energy accumulation, thereby posing a risk of rockburst occurrence. Therefore, it is of paramount importance to investigate the fracture propagation behavior in jointed rock masses and assess the impact of borehole pressure relief on mitigating rockburst occurrences for effective prevention and control measures. This paper focuses on failure characteristics and pressure relief effectiveness of non-persistent jointed rock mass with holes through laboratory testing and numerical simulation. In laboratory experiments, rock samples are prepared to include a range of crack dip angles and circular holes. Then, the crack propagation law of crack inclination and circular hole is studied by AE and DIC technology. The experimental results show that with the increase of fracture dip angle, the peak strength and energy change of the sample decrease first and then increase. Due to the existence of holes, the crack propagation direction of the original crack is changed. After drilling, the strain energy of the sample is obviously reduced, which shows that the drilling pressure relief effect is obvious, which can effectively reduce the energy accumulated inside the rock mass and reduce the risk of rockburst. Finally, the PFC numerical simulation software is used to analyze the micro-failure process and energy change law of the sample from three aspects: the relative position of cracks and holes, the diameter of boreholes and the spacing of boreholes. Further understanding of the energy dissipation law and mechanical behavior characteristics of jointed rock mass provides a reference for exploring the pressure relief effect of rock mass and preventing rockburst.