Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing addit...Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.展开更多
Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On ...Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On the other hand,the entry of ice crystal particles into the combustion chamber can cause a decrease in temperature or even flameout,leading to engine surge or shutdown.Therefore,it is necessary to conduct multiphase flow tests on ice crystals for aircraft components such as aircraft engines.Conducting ice crystal multiphase flow tests on aircraft is an effective research method,but it requires the construction of an ice crystal multiphase flow test platform that meets relevant technical requirements.The paper focuses on the relevant experimental requirements and combines wind tunnel test structures to conduct multiphase flow numerical simulations on various forms of jet pipelines,obtaining particle motion distribution results.After comparison,the optimal form of jet structure is obtained,providing the best selection scheme for the design of relevant wind tunnel structures.展开更多
Regenerative catalytic oxidizers(RCO)are widely used to remove volatile organic compounds(VOCs)due to their energy-saving and stability.In this study,a multi-component catalytic reaction model was constructed to numer...Regenerative catalytic oxidizers(RCO)are widely used to remove volatile organic compounds(VOCs)due to their energy-saving and stability.In this study,a multi-component catalytic reaction model was constructed to numerically investigate the reaction process of hydrocarbon-containing VOCs in RCO using computational fluid dynamics(CFD)simulation.To obtain the conversion characteristics of multi-component hydrocarbons,the effects of intake load,equivalence ratio,and the composition of multi-component hydrocarbons on the flow,heat transfer,and conversion rate of the reactor were analyzed.A feasibility study plan targeting the hard-to-convert components was also proposed.The results indicated that as the load increases,the conversion rates of the various components decrease,while the reaction rates increase.Moreover,increasing the flow velocity intensifies turbulence and enhances the collision frequency between the gas and the wall surfaces.This,in turn,amplifies the resistance effect of the porous medium.As the equivalence ratio of VOCs to oxygen increases,the oxygen-deficient condition leads to a decrease in the molecular weight of the hydrocarbons involved in the reaction.The reaction temperature also shows a downward trend.A comparative analysis of the catalytic combustion characteristics of multi-component VOCs and single-component gases reveals that adding ethane and propane can facilitate methane oxidation.展开更多
Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, w...Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, while neural networks struggle to represent underlying microscopic material properties. To overcome these challenges, a meso-micro scale numerical method using a virtual node approach is developed in this study. A Wbraid/Al/Epoxy functional structural material is fabricated, and a representative periodic unit cell is identified based on its architecture. The complex structure is then discretized into nodes, and mechanical interactions are governed by pre-defined computation rules. This virtual node method is systematically compared against both multiscale simulation and a neural network algorithm, with validation provided through mechanical experiments. The results demonstrate that the nodal operation strategy significantly reduces computational resource requirements. By quantifying microscopic bonding with coefficients, explicit interface treatment is avoided, granting the method strong adaptability to lattice materials. The method can simulate extremely complex structures using parameters from simple tests and is suited for large systems. Compared to three-point bending experiments, errors for multiscale, virtual node, and neural network methods were 12.4%, 6.9%, and 34.5%, respectively. Under dynamic compression, the errors were 2.7%, 9.3%, and 15.43%. The virtual node method demonstrated superior accuracy under static conditions, enabling efficient prediction and auxiliary development of complex structural materials.展开更多
Hydraulic fracturing serves as a critical technology for reservoir stimulation in deep coalbed methane(CBM)development,where the mechanical properties of gangue layers exert a significant control on fracture propagati...Hydraulic fracturing serves as a critical technology for reservoir stimulation in deep coalbed methane(CBM)development,where the mechanical properties of gangue layers exert a significant control on fracture propagation behavior.To address the unclear mechanisms governing fracture penetration across coal-gangue interfaces,this study employs the Continuum-Discontinuum Element Method(CDEM)to simulate and analyze the vertical propagation of hydraulic fractures initiating within coal seams,based on geomechanical parameters derived from the deep Benxi Formation coal seams in the southeastern Ordos Basin.The investigation systematically examines the influence of geological and operational parameters on cross-interfacial fracture growth.Results demonstrate that vertical stress difference,elastic modulus contrast between coal and gangue layers,interfacial stress differential,and interfacial cohesion at coal-gangue interfaces are critical factors governing hydraulic fracture penetration through these interfaces.High vertical stress differences(>3 MPa)inhibit interfacial dilation,promoting predominant crosslayer fracture propagation.Reduced interfacial stress contrasts and enhanced interfacial cohesion facilitate fracture penetration across interfaces.Furthermore,smaller elastic modulus contrasts between coal and gangue correlate with increased interfacial aperture.Finally,lower injection rates effectively suppress vertical fracture propagation in deep coal reservoirs.This study elucidates the characteristics and mechanisms governing cross-layer fracture propagation in coal–rock composites with interbedded partings,and delineates the dynamic evolution laws and dominant controlling factors involved.Thefindings provide critical theoretical insights for the optimization of fracture design and the efficient development of deep coalbed methane reservoirs.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
Based on the surrounding rock arching and hingeless arch structure theories,a theoretical formula for the minimum overburden thickness was derived.By substituting different mechanical parameters of multiple tunnels at...Based on the surrounding rock arching and hingeless arch structure theories,a theoretical formula for the minimum overburden thickness was derived.By substituting different mechanical parameters of multiple tunnels at home and abroad into this formula,minimum self-supporting arch formulas under different surrounding rock classes were obtained.Based on the actual engineering case of a dual-mode shield tunnel,a numerical model for the tunnel boring machine excavation mode was established to verify the theoretical formulas.Next,three surrounding rock classes,four soil layer thickness gradients,and twelve overburden thickness gradients were designed,resulting in 144 models formed by the combination of the three factors.Uniform tests were conducted,and the pressure arch heights under different surrounding rock classes were obtained.The results show that in the theoretical formulas,the tunnel radius has a linear positive correlation with the pressure arch height,while the tunnel depth has a linear positive correlation with the square of the pressure arch height.According to numerical simulation results,the pressure arch height increases with the increase of the overburden thickness and then tends toward a critical value of twice the tunnel diameter.Finally,the results of the numerical model are in good agreement with those calculated using the theoretical formulas,verifying the rationality of the established theoretical formulas.展开更多
Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers of...Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers offer advantages such as reduced material usage,lower refrigerant charge,and compact structure.However,they also face challenges,including increased refrigerant pressure drop and smaller heat transfer area inside the tubes.This paper combines the advantages and disadvantages of both small and large-diameter tubes and proposes a combined-diameter heat exchanger,consisting of large and small diameters,for use in the indoor units of split-type air conditioners.There are relatively few studies in this area.In this paper,A theoretical and numerical computation method is employed to establish a theoretical-numerical calculation model,and its reliability is verified through experiments.Using this model,the optimal combined diameters and flow path design for a combined-diameter heat exchanger using R32 as the working fluid are derived.The results show that the heat transfer performance of all combined diameter configurations improves by 2.79%to 8.26%compared to the baseline design,with the coefficient of performance(COP)increasing from 4.15 to 4.27~4.5.These designs can save copper material,but at the cost of an increase in pressure drop by 66.86%to 131.84%.The scheme IIIH,using R32,is the optimal combined-diameter and flow path configuration that balances both heat transfer performance and economic cost.展开更多
The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we pro...The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we propose a new approach by combining centrifugal analog modeling with numerical simulation to simulate the tectonic uplift history of the plateau based on the lower crustal flow model,and to investigate the material migration characteristics and the influence of crustal motion velocity and ductile layer viscosity on the plateau tectonic geomorphology.The models reproduce steep-sided flat-topped geomorphic features and clockwise rotation of the material at eastern Himalayan Syntaxis,verifying the rationality of the models.The results show that the greater the crustal motion velocity and the greater the ductile layer viscosity,the steeper the terrain change;and conversely,the smaller the crustal motion velocity and the smaller the ductile layer viscosity,the gentler the terrain change.This study further indicates that the weak lower crust plays an important role in the formation of geomorphic features and material migration characteristics of Qinghai-Tibet Plateau,and provides a new insight for the study of the uplift mechanism of the Tibetan Plateau.展开更多
Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions and has become a growing focus of research and deployment.This paper presents numerical simulations of CO_(2) injection and stora...Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions and has become a growing focus of research and deployment.This paper presents numerical simulations of CO_(2) injection and storage in a depleted gas reservoir within the B Depression and evaluates associated CO_(2) trapping mechanisms.In the base case,a constant injection rate of 3500 m^(3)/d over fifteen years resulted in a cumulative injection of 19.2×10^(6) m^(3).The CO_(2) plume expanded radially during injection and subsequently migrated up-dip under buoyancy forces.The final stored mass of CO_(2) in the reservoir was 10.6 million tonnes(Mt),representing less than 10% of its theoretical capacity.The plume was projected to reach the entrapment crest and the top of the reservoir within a century,indicating secure long-term containment.Structural,stratigraphic,and residual trapping dominate in Reservoir A(approximately 90%).Anticlinal closures with thick overlying mudstones in the Zhujiang Formation provide effective seals,further enhancing storage security.Reservoir properties and heterogeneity play a crucial role in controlling CO_(2) storage.However,reservoir heterogeneity exerts only a limited influence when intrinsic properties are favorable.Overall,the study and implementation of CO_(2) capture,utilization,and storage(CCUS)in China's offshore basins show promising prospects.展开更多
Focusing on the unclear mechanism of aerodynamic interference in overlapping rotors of heavy-load electric vertical take-off and landing(eVTOL)aircraft,this paper aims to reveal the aerodynamic interference characteri...Focusing on the unclear mechanism of aerodynamic interference in overlapping rotors of heavy-load electric vertical take-off and landing(eVTOL)aircraft,this paper aims to reveal the aerodynamic interference characteristics and flow field evolution laws of overlapping rotor configurations in hovering conditions through numerical simulation methods.The research method involves constructing a computational model for rotor flow fields and aerodynamic characteristics based on the Reynolds-averaged Navier-Stokes(RANS)equations and the Spalart-Allmaras(S-A)turbulence model.The dynamic simulation of rotor rotational motion was achieved by using the moving nested grid technology.The reliability of the computational method was ensured through the grid independence verification and the comparison with experimental data.The research results indicate that in overlapping rotor systems,rotorⅡexperiences a decrease in thrust,significant power fluctuations,and reduced hovering efficiency due to continuous interference from the adjacent rotor’s wake and blade-vortex interactions.Blade-tip vortices undergo breakage,fusion,and secondary rolling in the overlapping region,forming large-scale turbulent structures that lead to attenuation of the induced velocity field and aerodynamic efficiency losses.Additionally,the interaction between the rotor downwash and the fuselage triggers a“fountain effect”and a sudden increase in surface pressure on the fuselage,exacerbating flow field distortion.Based on the aforementioned mechanisms,the safe flight of overlapping rotor configurations can be achieved by optimizing the configuration strategy of the rotational speed phase difference between adjacent blades.This study provides a theoretical basis for the rotor layout design and the aerodynamic performance enhancement of heavy-load eVTOL aircraft.展开更多
The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combi...The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combines numerical simulation with machine learning techniques to explore this issue.It presents a summary of special-shaped tunnel geometries and introduces a shape coefficient.Through the finite element software,Plaxis3D,the study simulates six key parameters—shape coefficient,burial depth ratio,tunnel’s longest horizontal length,internal friction angle,cohesion,and soil submerged bulk density—that impact uplift resistance across different conditions.Employing XGBoost and ANN methods,the feature importance of each parameter was analyzed based on the numerical simulation results.The findings demonstrate that a tunnel shape more closely resembling a circle leads to reduced uplift resistance in the overlying soil,whereas other parameters exhibit the contrary effects.Furthermore,the study reveals a diminishing trend in the feature importance of buried depth ratio,internal friction angle,tunnel longest horizontal length,cohesion,soil submerged bulk density,and shape coefficient in influencing uplift resistance.展开更多
Traditional oilfields face increasing extraction challenges, primarily due to reservoir quality degradation and production decline, which are further exacerbated by volatile international crude oil prices—illustrated...Traditional oilfields face increasing extraction challenges, primarily due to reservoir quality degradation and production decline, which are further exacerbated by volatile international crude oil prices—illustrated by Brent Crude’s trajectory from pandemic-induced negative pricing to geopolitically driven surges exceeding USD 100 per barrel. This study addresses these complexities through an integrated methodological framework applied to medium-permeability sandstone reservoirs in the Xinjiang oilfield by combining advanced numerical simulations with multivariate regression analysis. The methodology employs Latin Hypercube Sampling (LHS) to stratify geological parameter distributions and constructs heterogeneous reservoir models using Petrel software, rigorously validated through historical production data matching. Production forecasting integrates numerical simulation and Decline Curve Analysis (DCA), while investment estimation utilizes Ordinary Least Squares (OLS) regression to correlate engineering parameters with drilling and completion costs. Economic evaluation incorporates Discounted Cash Flow (DCF) modeling and breakeven analysis, establishing techno-economic boundaries via oil price sensitivity analysis ranging from USD 40 to 90 per barrel. Visualization tools, including 3D heatmaps, delineate nonlinear interactions among engineering, geological, and investment datasets under economic constraints. Key findings demonstrate that for the target reservoirs, as oil prices increase from USD 40 to USD 90 per barrel, the minimum economic thickness threshold decreases from approximately 5.7 m to about 2.5 m, with model prediction errors consistently below 25% across validation datasets. This framework provides scientifically grounded decision support for optimizing capital allocation and offers actionable insights to enhance undeveloped hydrocarbon development planning amid market uncertainty. Ultimately, it supports national energy security through technically robust and economically viable resource exploitation strategies.展开更多
This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the pred...This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the prediction of the movement track and intensity of Typhoon Kompasu in 2021 is examined.Additionally,the possible reasons for their effects on tropical cyclone(TC)intensity prediction are analyzed.Statistical results show that both parameterization schemes improve the predictions of Typhoon Kompasu’s track and intensity.The influence on track prediction becomes evident after 60 h of model integration,while the significant positive impact on intensity prediction is observed after 66 h.Further analysis reveals that these two schemes affect the timing and magnitude of extreme TC intensity values by influencing the evolution of the TC’s warm-core structure.展开更多
In sediment transport on a mobile erodible bed,near-bed particles tend to preferentially locate under specific flow regions and form an uneven bedform.These moving and mobile particles can significantly modulate turbu...In sediment transport on a mobile erodible bed,near-bed particles tend to preferentially locate under specific flow regions and form an uneven bedform.These moving and mobile particles can significantly modulate turbulence at various scales,from inertial large-scale motions to small viscous motions.This study analyzes the particle-resolved direct numerical simulation data of particle-laden two-phase flow with multilayers of particles in turbulent flows over static and mobile beds.The double-average method is adopted for energy transfer analysis.The result shows that the alternative streaky bedform in the spanwise directions correlates with the streak structures in the near-wall turbulence in the mobile bed case.The energy redistribution and exchange,as well as the dissipation,are analyzed in detail,and an energy transfer diagram is given in the last to summarize the energy transfer processes.In both the static and mobile bed cases,flow energy is introduced into the system via the work performed by volume forces acting on the mean flow.The viscous dissipations in the double mean and form-induced fields are more pronounced in the static bed case,and the work done by the fluid-particle interfacial stress in the double mean and turbulent fields is more pronounced in the mobile bed case.The prominent energy contribution in the form-induced field is the production by the form-induced stress on the mean strain in the mobile bed case.In addition,sediment transport involving a limited number of mobile particles is insufficient to capture the energy transfer processes that occur over the troughs and may intertwine the energy transfer processes over the mobile particles and the fixed particle bed.展开更多
This study presents an investigation into shock-induced exothermic reactions within three distinct aluminum-based energetic mixtures:aluminum/sulfur(Al/S),aluminum/copper oxide(Al/CuO),and aluminum/polytetrafluoroethy...This study presents an investigation into shock-induced exothermic reactions within three distinct aluminum-based energetic mixtures:aluminum/sulfur(Al/S),aluminum/copper oxide(Al/CuO),and aluminum/polytetrafluoroethylene(Al/PTFE).A challenge in current modeling efforts is accurately capturing the complex physical and chemical coupling under extreme loading,especially the influence of rapidly forming gaseous products in Al/PTFE mixtures on material integrity.To address this,a wide-range numerical model based on the Smoothed Particle Hydrodynamics(SPH)method was developed.This mesh-free approach manages large deformations and incorporates elastic-plastic flow,heat transfer,component diffusion,and chemical kinetics simulated using both zero-and first-order reaction schemes,favoring the latter for surface-reaction mechanisms.The proposed model takes into account gaseous reaction products,specifically aluminum fluoride(AlF3)to assess their impact on ampoule fracture dynamics.Numerical simulations,validated against experimental data,demonstrated that reaction rate,local pressure,and temperature are the primary controlling factors governing energy release and structural response.Comparative analysis revealed that although Al/CuO initiates reaction more readily(lower critical pressure/temperature),the Al/S mixture exhibits superior overall reaction efficiency under shock-wave loading,highlighting the significance of post-initiation kinetic factors.Furthermore,simulations using the conical ampoule geometry confirmed its effectiveness in generating a continuous pressure gradient,enabling systematic characterization of pressure-dependent reaction kinetics.This validated SPH model provides a powerful and predictive tool for understanding the complex behavior of energetic materials under shock-wave loading and aids in optimizing material composition for desired performance characteristics.展开更多
Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical si...Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical simulation plays a significant role in quantitatively evaluating current processes and making targeted improvements,but its limitations lie in the inability to dynamically reflect the formation outcomes of castings under varying process conditions,making real-time adjustments to gating and riser designs challenging.In this study,an automated design model for gating and riser systems based on integrated parametric 3D modeling-simulation framework is proposed,which enhances the flexibility and usability of evaluating the casting process by simulation.Firstly,geometric feature extraction technology is employed to obtain the geometric information of the target casting.Based on this information,an automated design framework for gating and riser systems is established,incorporating multiple structural parameters for real-time process control.Subsequently,the simulation results for various structural parameters are analyzed,and the influence of these parameters on casting formation is thoroughly investigated.Finally,the optimal design scheme is generated and validated through experimental verification.Simulation analysis and experimental results show that using a larger gate neck(24 mm in side length) and external risers promotes a more uniform temperature distribution and a more stable flow state,effectively eliminating shrinkage cavities and enhancing process yield by 15%.展开更多
Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized...Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized coal co-injection in blast furnace tuyere was established through numerical simulation,and the effect of hydrogen injection and oxygen enrichment interaction on pulverized coal combustion and raceway smelting was investigated.The simulation results indicate that when the coal injection rate decreased from 36 to 30t/h and the hydrogen injection increased from 0 to 3600 m^(3)/h,the CO_(2)emissions decreased from 1860 to 1551 kg/t,which represents a16.6%reduction,and the pulverized coal burnout decreased from 70.1%to 63.7%.The heat released from hydrogen combustion can not only promote the volatilization of pulverized coal but also affect the combustion reaction between volatilization and oxygen,which resulted in a decrease in the temperature at the end of the raceway.Co-injection of hydrogen with PCI increased the wall temperature near the upper half part of the raceway and at the outlet of the tuyere,which required a high cooling efficiency to extend the service life of the blast furnace.The increase in oxygen level compensated for the decreased average temperature in the raceway due to hydrogen injection.The increase in the oxygen content by 3%while maintaining constant hydrogen and PCI injection rates increased the burnout and average raceway temperature by 4.2%and 43 K,respectively.The mole fraction of CO and H_(2) production increased by 0.04 and 0.02,respectively.Burnout can be improved through optimization of the particle size distribution of pulverized coal.展开更多
Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a s...Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a suitable approach to predict the rockburst tendency at the preliminary stage becomes very important.In this study,an integrated methodology combining 3D initial stress inversion and rockburst tendency prediction was developed and subsequently applied to a case study of the Sangzhuling Tunnel on the Sichuan–Tibet Railway.The numerical modelling involved inverting the initial stress field using a multiple linear regression method.The tunnel excavation was simulated separately by FDM and DEM,based on a stress boundary condition from the inverted stress field.The comparative analysis demonstrates that the rockburst ratio calculated using DEM(76.70%)exhibits a slight increase compared to FDM(75.38%),and the rockburst location is consistent with the actual situation.This suggests that DEM is more suitable for simulating the stress redistribution during excavation in a jointed rock mass.The numerical simulation combined with the deviatoric stress approach effectively predicts rockburst tendency,meeting the engineering requirements.Despite its limitations,numerical simulation remains a reliable method for predicting rock bursts.展开更多
基金National Key Research and Development Program of China(2022YFB4600902)Shandong Provincial Science Foundation for Outstanding Young Scholars(ZR2024YQ020)。
文摘Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.
文摘Ice crystal icing is an important cause of accidents in aircraft engines.Ice formation in aircraft engines can cause internal blades to freeze,affecting the quality of the air flow field and blocking the flow path.On the other hand,the entry of ice crystal particles into the combustion chamber can cause a decrease in temperature or even flameout,leading to engine surge or shutdown.Therefore,it is necessary to conduct multiphase flow tests on ice crystals for aircraft components such as aircraft engines.Conducting ice crystal multiphase flow tests on aircraft is an effective research method,but it requires the construction of an ice crystal multiphase flow test platform that meets relevant technical requirements.The paper focuses on the relevant experimental requirements and combines wind tunnel test structures to conduct multiphase flow numerical simulations on various forms of jet pipelines,obtaining particle motion distribution results.After comparison,the optimal form of jet structure is obtained,providing the best selection scheme for the design of relevant wind tunnel structures.
基金supported by National Key Research&Development Program of China(2022YFB4101500).
文摘Regenerative catalytic oxidizers(RCO)are widely used to remove volatile organic compounds(VOCs)due to their energy-saving and stability.In this study,a multi-component catalytic reaction model was constructed to numerically investigate the reaction process of hydrocarbon-containing VOCs in RCO using computational fluid dynamics(CFD)simulation.To obtain the conversion characteristics of multi-component hydrocarbons,the effects of intake load,equivalence ratio,and the composition of multi-component hydrocarbons on the flow,heat transfer,and conversion rate of the reactor were analyzed.A feasibility study plan targeting the hard-to-convert components was also proposed.The results indicated that as the load increases,the conversion rates of the various components decrease,while the reaction rates increase.Moreover,increasing the flow velocity intensifies turbulence and enhances the collision frequency between the gas and the wall surfaces.This,in turn,amplifies the resistance effect of the porous medium.As the equivalence ratio of VOCs to oxygen increases,the oxygen-deficient condition leads to a decrease in the molecular weight of the hydrocarbons involved in the reaction.The reaction temperature also shows a downward trend.A comparative analysis of the catalytic combustion characteristics of multi-component VOCs and single-component gases reveals that adding ethane and propane can facilitate methane oxidation.
文摘Existing numerical methods for complex composites, such as multiscale simulation and neural network algorithms, face significant limitations. Multiscale techniques are often prohibitively expensive for large models, while neural networks struggle to represent underlying microscopic material properties. To overcome these challenges, a meso-micro scale numerical method using a virtual node approach is developed in this study. A Wbraid/Al/Epoxy functional structural material is fabricated, and a representative periodic unit cell is identified based on its architecture. The complex structure is then discretized into nodes, and mechanical interactions are governed by pre-defined computation rules. This virtual node method is systematically compared against both multiscale simulation and a neural network algorithm, with validation provided through mechanical experiments. The results demonstrate that the nodal operation strategy significantly reduces computational resource requirements. By quantifying microscopic bonding with coefficients, explicit interface treatment is avoided, granting the method strong adaptability to lattice materials. The method can simulate extremely complex structures using parameters from simple tests and is suited for large systems. Compared to three-point bending experiments, errors for multiscale, virtual node, and neural network methods were 12.4%, 6.9%, and 34.5%, respectively. Under dynamic compression, the errors were 2.7%, 9.3%, and 15.43%. The virtual node method demonstrated superior accuracy under static conditions, enabling efficient prediction and auxiliary development of complex structural materials.
文摘Hydraulic fracturing serves as a critical technology for reservoir stimulation in deep coalbed methane(CBM)development,where the mechanical properties of gangue layers exert a significant control on fracture propagation behavior.To address the unclear mechanisms governing fracture penetration across coal-gangue interfaces,this study employs the Continuum-Discontinuum Element Method(CDEM)to simulate and analyze the vertical propagation of hydraulic fractures initiating within coal seams,based on geomechanical parameters derived from the deep Benxi Formation coal seams in the southeastern Ordos Basin.The investigation systematically examines the influence of geological and operational parameters on cross-interfacial fracture growth.Results demonstrate that vertical stress difference,elastic modulus contrast between coal and gangue layers,interfacial stress differential,and interfacial cohesion at coal-gangue interfaces are critical factors governing hydraulic fracture penetration through these interfaces.High vertical stress differences(>3 MPa)inhibit interfacial dilation,promoting predominant crosslayer fracture propagation.Reduced interfacial stress contrasts and enhanced interfacial cohesion facilitate fracture penetration across interfaces.Furthermore,smaller elastic modulus contrasts between coal and gangue correlate with increased interfacial aperture.Finally,lower injection rates effectively suppress vertical fracture propagation in deep coal reservoirs.This study elucidates the characteristics and mechanisms governing cross-layer fracture propagation in coal–rock composites with interbedded partings,and delineates the dynamic evolution laws and dominant controlling factors involved.Thefindings provide critical theoretical insights for the optimization of fracture design and the efficient development of deep coalbed methane reservoirs.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金The National Natural Science Foundation of China(No.52478426)the Natural Science Foundation of Hunan Province(No.2024JJ5428).
文摘Based on the surrounding rock arching and hingeless arch structure theories,a theoretical formula for the minimum overburden thickness was derived.By substituting different mechanical parameters of multiple tunnels at home and abroad into this formula,minimum self-supporting arch formulas under different surrounding rock classes were obtained.Based on the actual engineering case of a dual-mode shield tunnel,a numerical model for the tunnel boring machine excavation mode was established to verify the theoretical formulas.Next,three surrounding rock classes,four soil layer thickness gradients,and twelve overburden thickness gradients were designed,resulting in 144 models formed by the combination of the three factors.Uniform tests were conducted,and the pressure arch heights under different surrounding rock classes were obtained.The results show that in the theoretical formulas,the tunnel radius has a linear positive correlation with the pressure arch height,while the tunnel depth has a linear positive correlation with the square of the pressure arch height.According to numerical simulation results,the pressure arch height increases with the increase of the overburden thickness and then tends toward a critical value of twice the tunnel diameter.Finally,the results of the numerical model are in good agreement with those calculated using the theoretical formulas,verifying the rationality of the established theoretical formulas.
基金supported by Supported by the Scientific Research Foundation for High-Level Talents of Zhoukou Normal University(ZKNUC2024018).
文摘Energy shortage has become one of themost concerning issues in the world today,and improving energy utilization efficiency is a key area of research for experts and scholars worldwide.Small-diameter heat exchangers offer advantages such as reduced material usage,lower refrigerant charge,and compact structure.However,they also face challenges,including increased refrigerant pressure drop and smaller heat transfer area inside the tubes.This paper combines the advantages and disadvantages of both small and large-diameter tubes and proposes a combined-diameter heat exchanger,consisting of large and small diameters,for use in the indoor units of split-type air conditioners.There are relatively few studies in this area.In this paper,A theoretical and numerical computation method is employed to establish a theoretical-numerical calculation model,and its reliability is verified through experiments.Using this model,the optimal combined diameters and flow path design for a combined-diameter heat exchanger using R32 as the working fluid are derived.The results show that the heat transfer performance of all combined diameter configurations improves by 2.79%to 8.26%compared to the baseline design,with the coefficient of performance(COP)increasing from 4.15 to 4.27~4.5.These designs can save copper material,but at the cost of an increase in pressure drop by 66.86%to 131.84%.The scheme IIIH,using R32,is the optimal combined-diameter and flow path configuration that balances both heat transfer performance and economic cost.
基金supported by Excellent Research Group Project for Multiphase Evolution in Hyper-Gravity of the National Natural Science Foundation of China(No.52588202)。
文摘The India-Asia collision resulted in the formation of Qinghai-Tibet Plateau.Lower crustal flow model was proposed to explain the mechanism of Cenozoic tectonic deformation of Qinghai-Tibet Plateau.In this study,we propose a new approach by combining centrifugal analog modeling with numerical simulation to simulate the tectonic uplift history of the plateau based on the lower crustal flow model,and to investigate the material migration characteristics and the influence of crustal motion velocity and ductile layer viscosity on the plateau tectonic geomorphology.The models reproduce steep-sided flat-topped geomorphic features and clockwise rotation of the material at eastern Himalayan Syntaxis,verifying the rationality of the models.The results show that the greater the crustal motion velocity and the greater the ductile layer viscosity,the steeper the terrain change;and conversely,the smaller the crustal motion velocity and the smaller the ductile layer viscosity,the gentler the terrain change.This study further indicates that the weak lower crust plays an important role in the formation of geomorphic features and material migration characteristics of Qinghai-Tibet Plateau,and provides a new insight for the study of the uplift mechanism of the Tibetan Plateau.
基金funded by the National Natural Science Foundation of China(Grant No.42102169)the Shaanxi Province key Research and Development Project(Grant No.2023-ZDLSF-64)+1 种基金the Scientific Research and Technological Development Project of China National Logging Corporation(Grant No.25ZYCJSG013-2504)the Youth Science and Technology Special Fund of PetroChina(Grant No.2024DQ03172)。
文摘Geological CO_(2) storage is a promising strategy for reducing greenhouse gas emissions and has become a growing focus of research and deployment.This paper presents numerical simulations of CO_(2) injection and storage in a depleted gas reservoir within the B Depression and evaluates associated CO_(2) trapping mechanisms.In the base case,a constant injection rate of 3500 m^(3)/d over fifteen years resulted in a cumulative injection of 19.2×10^(6) m^(3).The CO_(2) plume expanded radially during injection and subsequently migrated up-dip under buoyancy forces.The final stored mass of CO_(2) in the reservoir was 10.6 million tonnes(Mt),representing less than 10% of its theoretical capacity.The plume was projected to reach the entrapment crest and the top of the reservoir within a century,indicating secure long-term containment.Structural,stratigraphic,and residual trapping dominate in Reservoir A(approximately 90%).Anticlinal closures with thick overlying mudstones in the Zhujiang Formation provide effective seals,further enhancing storage security.Reservoir properties and heterogeneity play a crucial role in controlling CO_(2) storage.However,reservoir heterogeneity exerts only a limited influence when intrinsic properties are favorable.Overall,the study and implementation of CO_(2) capture,utilization,and storage(CCUS)in China's offshore basins show promising prospects.
基金supported by the National Natural Science Foundation of China(No.11872211)。
文摘Focusing on the unclear mechanism of aerodynamic interference in overlapping rotors of heavy-load electric vertical take-off and landing(eVTOL)aircraft,this paper aims to reveal the aerodynamic interference characteristics and flow field evolution laws of overlapping rotor configurations in hovering conditions through numerical simulation methods.The research method involves constructing a computational model for rotor flow fields and aerodynamic characteristics based on the Reynolds-averaged Navier-Stokes(RANS)equations and the Spalart-Allmaras(S-A)turbulence model.The dynamic simulation of rotor rotational motion was achieved by using the moving nested grid technology.The reliability of the computational method was ensured through the grid independence verification and the comparison with experimental data.The research results indicate that in overlapping rotor systems,rotorⅡexperiences a decrease in thrust,significant power fluctuations,and reduced hovering efficiency due to continuous interference from the adjacent rotor’s wake and blade-vortex interactions.Blade-tip vortices undergo breakage,fusion,and secondary rolling in the overlapping region,forming large-scale turbulent structures that lead to attenuation of the induced velocity field and aerodynamic efficiency losses.Additionally,the interaction between the rotor downwash and the fuselage triggers a“fountain effect”and a sudden increase in surface pressure on the fuselage,exacerbating flow field distortion.Based on the aforementioned mechanisms,the safe flight of overlapping rotor configurations can be achieved by optimizing the configuration strategy of the rotational speed phase difference between adjacent blades.This study provides a theoretical basis for the rotor layout design and the aerodynamic performance enhancement of heavy-load eVTOL aircraft.
基金Guangzhou Metro Scientific Research Project(No.JT204-100111-23001)Chongqing Municipal Special Project for Technological Innovation and Application Development(No.CSTB2022TIAD-KPX0101)Science and Technology Research and Development Program of China State Railway Group Co.,Ltd.(No.N2023G045)。
文摘The uplift resistance of the soil overlying shield tunnels significantly impacts their anti-floating stability.However,research on uplift resistance concerning special-shaped shield tunnels is limited.This study combines numerical simulation with machine learning techniques to explore this issue.It presents a summary of special-shaped tunnel geometries and introduces a shape coefficient.Through the finite element software,Plaxis3D,the study simulates six key parameters—shape coefficient,burial depth ratio,tunnel’s longest horizontal length,internal friction angle,cohesion,and soil submerged bulk density—that impact uplift resistance across different conditions.Employing XGBoost and ANN methods,the feature importance of each parameter was analyzed based on the numerical simulation results.The findings demonstrate that a tunnel shape more closely resembling a circle leads to reduced uplift resistance in the overlying soil,whereas other parameters exhibit the contrary effects.Furthermore,the study reveals a diminishing trend in the feature importance of buried depth ratio,internal friction angle,tunnel longest horizontal length,cohesion,soil submerged bulk density,and shape coefficient in influencing uplift resistance.
文摘Traditional oilfields face increasing extraction challenges, primarily due to reservoir quality degradation and production decline, which are further exacerbated by volatile international crude oil prices—illustrated by Brent Crude’s trajectory from pandemic-induced negative pricing to geopolitically driven surges exceeding USD 100 per barrel. This study addresses these complexities through an integrated methodological framework applied to medium-permeability sandstone reservoirs in the Xinjiang oilfield by combining advanced numerical simulations with multivariate regression analysis. The methodology employs Latin Hypercube Sampling (LHS) to stratify geological parameter distributions and constructs heterogeneous reservoir models using Petrel software, rigorously validated through historical production data matching. Production forecasting integrates numerical simulation and Decline Curve Analysis (DCA), while investment estimation utilizes Ordinary Least Squares (OLS) regression to correlate engineering parameters with drilling and completion costs. Economic evaluation incorporates Discounted Cash Flow (DCF) modeling and breakeven analysis, establishing techno-economic boundaries via oil price sensitivity analysis ranging from USD 40 to 90 per barrel. Visualization tools, including 3D heatmaps, delineate nonlinear interactions among engineering, geological, and investment datasets under economic constraints. Key findings demonstrate that for the target reservoirs, as oil prices increase from USD 40 to USD 90 per barrel, the minimum economic thickness threshold decreases from approximately 5.7 m to about 2.5 m, with model prediction errors consistently below 25% across validation datasets. This framework provides scientifically grounded decision support for optimizing capital allocation and offers actionable insights to enhance undeveloped hydrocarbon development planning amid market uncertainty. Ultimately, it supports national energy security through technically robust and economically viable resource exploitation strategies.
基金supported by the National Key R&D Program of China[grant number 2023YFC3008004]。
文摘This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the prediction of the movement track and intensity of Typhoon Kompasu in 2021 is examined.Additionally,the possible reasons for their effects on tropical cyclone(TC)intensity prediction are analyzed.Statistical results show that both parameterization schemes improve the predictions of Typhoon Kompasu’s track and intensity.The influence on track prediction becomes evident after 60 h of model integration,while the significant positive impact on intensity prediction is observed after 66 h.Further analysis reveals that these two schemes affect the timing and magnitude of extreme TC intensity values by influencing the evolution of the TC’s warm-core structure.
基金supported by the National Natural Science Foundation of China(Grant Nos.12388101 and 12472221)the Fundamental Research Funds for the Central Universities(Grant No.lzujbky-2024-oy10).
文摘In sediment transport on a mobile erodible bed,near-bed particles tend to preferentially locate under specific flow regions and form an uneven bedform.These moving and mobile particles can significantly modulate turbulence at various scales,from inertial large-scale motions to small viscous motions.This study analyzes the particle-resolved direct numerical simulation data of particle-laden two-phase flow with multilayers of particles in turbulent flows over static and mobile beds.The double-average method is adopted for energy transfer analysis.The result shows that the alternative streaky bedform in the spanwise directions correlates with the streak structures in the near-wall turbulence in the mobile bed case.The energy redistribution and exchange,as well as the dissipation,are analyzed in detail,and an energy transfer diagram is given in the last to summarize the energy transfer processes.In both the static and mobile bed cases,flow energy is introduced into the system via the work performed by volume forces acting on the mean flow.The viscous dissipations in the double mean and form-induced fields are more pronounced in the static bed case,and the work done by the fluid-particle interfacial stress in the double mean and turbulent fields is more pronounced in the mobile bed case.The prominent energy contribution in the form-induced field is the production by the form-induced stress on the mean strain in the mobile bed case.In addition,sediment transport involving a limited number of mobile particles is insufficient to capture the energy transfer processes that occur over the troughs and may intertwine the energy transfer processes over the mobile particles and the fixed particle bed.
基金conducted under the public contract for the Tomsk Scientific Center of the Siberian Branch of the Russian Academy of Sciences(Project No.FWRF-2024-0009).
文摘This study presents an investigation into shock-induced exothermic reactions within three distinct aluminum-based energetic mixtures:aluminum/sulfur(Al/S),aluminum/copper oxide(Al/CuO),and aluminum/polytetrafluoroethylene(Al/PTFE).A challenge in current modeling efforts is accurately capturing the complex physical and chemical coupling under extreme loading,especially the influence of rapidly forming gaseous products in Al/PTFE mixtures on material integrity.To address this,a wide-range numerical model based on the Smoothed Particle Hydrodynamics(SPH)method was developed.This mesh-free approach manages large deformations and incorporates elastic-plastic flow,heat transfer,component diffusion,and chemical kinetics simulated using both zero-and first-order reaction schemes,favoring the latter for surface-reaction mechanisms.The proposed model takes into account gaseous reaction products,specifically aluminum fluoride(AlF3)to assess their impact on ampoule fracture dynamics.Numerical simulations,validated against experimental data,demonstrated that reaction rate,local pressure,and temperature are the primary controlling factors governing energy release and structural response.Comparative analysis revealed that although Al/CuO initiates reaction more readily(lower critical pressure/temperature),the Al/S mixture exhibits superior overall reaction efficiency under shock-wave loading,highlighting the significance of post-initiation kinetic factors.Furthermore,simulations using the conical ampoule geometry confirmed its effectiveness in generating a continuous pressure gradient,enabling systematic characterization of pressure-dependent reaction kinetics.This validated SPH model provides a powerful and predictive tool for understanding the complex behavior of energetic materials under shock-wave loading and aids in optimizing material composition for desired performance characteristics.
基金financially supported by the National Key Research and Development Program of China (2022YFB3706802)。
文摘Automation and intelligence have become the primary trends in the design of investment casting processes.However,the design of gating and riser systems still lacks precise quantitative evaluation criteria.Numerical simulation plays a significant role in quantitatively evaluating current processes and making targeted improvements,but its limitations lie in the inability to dynamically reflect the formation outcomes of castings under varying process conditions,making real-time adjustments to gating and riser designs challenging.In this study,an automated design model for gating and riser systems based on integrated parametric 3D modeling-simulation framework is proposed,which enhances the flexibility and usability of evaluating the casting process by simulation.Firstly,geometric feature extraction technology is employed to obtain the geometric information of the target casting.Based on this information,an automated design framework for gating and riser systems is established,incorporating multiple structural parameters for real-time process control.Subsequently,the simulation results for various structural parameters are analyzed,and the influence of these parameters on casting formation is thoroughly investigated.Finally,the optimal design scheme is generated and validated through experimental verification.Simulation analysis and experimental results show that using a larger gate neck(24 mm in side length) and external risers promotes a more uniform temperature distribution and a more stable flow state,effectively eliminating shrinkage cavities and enhancing process yield by 15%.
基金financially supported by the National Natural Science Foundation of China(No.51904026)the Fundamental Research Funds for the Central Universities(No.06500108)。
文摘Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized coal co-injection in blast furnace tuyere was established through numerical simulation,and the effect of hydrogen injection and oxygen enrichment interaction on pulverized coal combustion and raceway smelting was investigated.The simulation results indicate that when the coal injection rate decreased from 36 to 30t/h and the hydrogen injection increased from 0 to 3600 m^(3)/h,the CO_(2)emissions decreased from 1860 to 1551 kg/t,which represents a16.6%reduction,and the pulverized coal burnout decreased from 70.1%to 63.7%.The heat released from hydrogen combustion can not only promote the volatilization of pulverized coal but also affect the combustion reaction between volatilization and oxygen,which resulted in a decrease in the temperature at the end of the raceway.Co-injection of hydrogen with PCI increased the wall temperature near the upper half part of the raceway and at the outlet of the tuyere,which required a high cooling efficiency to extend the service life of the blast furnace.The increase in oxygen level compensated for the decreased average temperature in the raceway due to hydrogen injection.The increase in the oxygen content by 3%while maintaining constant hydrogen and PCI injection rates increased the burnout and average raceway temperature by 4.2%and 43 K,respectively.The mole fraction of CO and H_(2) production increased by 0.04 and 0.02,respectively.Burnout can be improved through optimization of the particle size distribution of pulverized coal.
基金financially supported by the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection(Chengdu University of Technology)(Grant No.SKLGP2020Z007)。
文摘Deeply buried mountain tunnels are often exposed to the risk of rock bursts,which always cause serious damage to the supporting structures and threaten the safety of the engineers.Due to the limited data available,a suitable approach to predict the rockburst tendency at the preliminary stage becomes very important.In this study,an integrated methodology combining 3D initial stress inversion and rockburst tendency prediction was developed and subsequently applied to a case study of the Sangzhuling Tunnel on the Sichuan–Tibet Railway.The numerical modelling involved inverting the initial stress field using a multiple linear regression method.The tunnel excavation was simulated separately by FDM and DEM,based on a stress boundary condition from the inverted stress field.The comparative analysis demonstrates that the rockburst ratio calculated using DEM(76.70%)exhibits a slight increase compared to FDM(75.38%),and the rockburst location is consistent with the actual situation.This suggests that DEM is more suitable for simulating the stress redistribution during excavation in a jointed rock mass.The numerical simulation combined with the deviatoric stress approach effectively predicts rockburst tendency,meeting the engineering requirements.Despite its limitations,numerical simulation remains a reliable method for predicting rock bursts.