In order to study the effect of rotation on the combustion in the underwater vehicle,a two-phase turbulent combustion process is described with Reynolds stress turbulence model,eddy-dissipation turbulent combustion mo...In order to study the effect of rotation on the combustion in the underwater vehicle,a two-phase turbulent combustion process is described with Reynolds stress turbulence model,eddy-dissipation turbulent combustion model,P-1 radiation model and particle tracking model of liquid. The flow in the rotating combustion chamber is simulated at two different working speeds,0?r/min and 1?000?r/min by Fluent software. The temperature,gas velocity,static pressure of wall and fuel concentration are computed and compared. The results show that the combustion in rotating combustor is faster and more effective.展开更多
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.展开更多
With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving larg...With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.展开更多
In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fi...In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.展开更多
The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to ...The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to achieving controllable stress-strain rate loading.In this study,we have,for the first time,combined one-dimensional fluid computational software with machine learning methods.We first elucidated the mechanisms by which GDI structures control stress and strain rates.Subsequently,we constructed a machine learning model to create a structure-property response surface.The results show that altering the loading velocity and interlayer thickness has a pronounced regulatory effect on stress and strain rates.In contrast,the impedance distribution index and target thickness have less significant effects on stress regulation,although there is a matching relationship between target thickness and interlayer thickness.Compared with traditional design methods,the machine learning approach offers a10^(4)—10^(5)times increase in efficiency and the potential to achieve a global optimum,holding promise for guiding the design of GDI.展开更多
To elucidate the influence of confining pressure on microcrack evolution and macroscopic failure mechanisms in granite,a multi-perspective approach was adopted.This approach combined triaxial compression tests,acousti...To elucidate the influence of confining pressure on microcrack evolution and macroscopic failure mechanisms in granite,a multi-perspective approach was adopted.This approach combined triaxial compression tests,acoustic emission(AE)monitoring,and PFC simulations.The results show that:1)Confining pressure exhibits a pronounced linear correlation with both yield strength and compressive strength.The enhancement of confining pressure significantly improves the deformability of granite and promotes a progressive shift in failure mechanism from brittle rupture to ductile deformation;2)Increasing confining pressure elevates the stress threshold for microcrack initiation and suppresses crack propagation.As a result,the proportion of shear cracks increases(based on AE analysis)from 18.71%to 61.2%,marking a transition in the dominant failure mode from tensile to shear;3)Confining pressure facilitates the development of grain boundary shear cracks(GBSCs),establishing the primary damage pathways.In addition,local stress concentrations under high confinement conditions trigger intragranular cracking.This highlights the regulatory effect of confining pressure on microcrack evolution.展开更多
Proper room and pillar sizes are both critical factors for safe mining and high ore recovery rate in shrinkage stoping mining of underground metal mines.The rock masses of Tangdan copper mine of China are fractured,wh...Proper room and pillar sizes are both critical factors for safe mining and high ore recovery rate in shrinkage stoping mining of underground metal mines.The rock masses of Tangdan copper mine of China are fractured,which needs much reinforcement and support prior to mining.Cement-sodium silicate grout technology was selected,then its related parameters such as grout pressure,diffusion radius and time were calculated and proposed.In order to test the effect of the pressured grout in the fractured No.4 ore block,field experiments were conducted.To optimize stoping configuration,three-dimensional numerical simulation with ANSYS and FLAC 3 D softwares was proposed.The results show that the drilling porosity and mechanical properties of the rock masses are increased obviously.After grout,ore recovery rate is increased by 10.2%employing the newly designed stoping configuration compared with the previous.Last,analyzed from the surface movements,roof subsidence and the maximum principal stress of the pillars,the mining safety is probable of being ensured.展开更多
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.展开更多
We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficienc...We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.展开更多
In recent twenty years, much numerical simulation work has been done on the evolution of Qinghai-Xizang (Tibetan) plateau. In this paper some principal numerical models and results are reviewed and analyzed. The earli...In recent twenty years, much numerical simulation work has been done on the evolution of Qinghai-Xizang (Tibetan) plateau. In this paper some principal numerical models and results are reviewed and analyzed. The earlier plane stress or plane strain model has much discrepancy with the actual deformation of Qinghai-Xizang plateau, such as the thickening of Tibetan crust and the lateral extrusion of Tibet along strike-slip faults. The thin viscous sheet model and the thin-plate model may simulate the change of the crustal thickness and the deformation pro-duced by gravitational force. It is suitable for studying the large-scale and long-time deformation. The influence of faults on the deformation of Tibetan plateau should be further studied.展开更多
Distillation temperature,as a pivotal thermodynamic parameter in vacuum purification of metals,governs impurity migration through volatility-stratified mechanisms.This study establishes theoretical distribution models...Distillation temperature,as a pivotal thermodynamic parameter in vacuum purification of metals,governs impurity migration through volatility-stratified mechanisms.This study establishes theoretical distribution models for high-volatility(Na and Se),medium-volatility(Fe and Cu),and low-volatility(Ni and Cr)impurities,revealing dual threshold effects on impurity removal:low-to-medium temperatures(£550℃)effectively suppress volatilization,while elevated temperatures promote co-evaporation.At the optimal 550℃,tellurium purity reaches 5 N 8 with>90%yield.Spatial fractionation analysis demonstrates high-volatility impurities enriching in the upper condensation zone(X/L<0.25),whereas medium/low-volatility impurities accumulate in the lower zone and residues.Remarkably,residual impurities show significant enrichment versus raw materials-Na 4.04 times,Fe 13.3 times,Cu 53.6 times,Ni 7.17 times,and Cr 15.12 times,through formation of non-volatile compounds/solid solutions.Temperature-space coupling effects drive distinct deviation patterns:fluctuations in high-volatility impurities vs.temperature-progressive deviations in medium/low-volatility species.The impurity concentration of model-experiment deviations(mean±SD)are quantified as:Se 0.265±0.12,Na 0.224±0.02,Cu 0.146±0.06,Fe 0.133±0.13,Cr 0.101±0.07,and Ni 0.122±0.08.展开更多
The present study is devoted to researching the thermal security problems of large-scale solid rocket motor with Ammonium perchlorate/Hydroxyl-terminated polybutadiene(AP/HTPB). A two-dimensional axisymmetric model fo...The present study is devoted to researching the thermal security problems of large-scale solid rocket motor with Ammonium perchlorate/Hydroxyl-terminated polybutadiene(AP/HTPB). A two-dimensional axisymmetric model for the cook-off of solid rocket motor is established. The reaction kinetics for the cook-off process of AP/HTPB is described by the two-step global chemical mechanism. Numerical predictions of the cook-off behavior for the propellant are conducted at fast heating rate of 1.45-2.45 K/s,and slow heating rate of 0.001-0.003 K/s, respectively. The results show that in the fast cook-off condition. the initial ignition position of AP/HTPB occurs in the annular region of the outer wall of propellant without exception, and the center point in the region is(889.1,149.5). For the region, the axial width is1.8 mm and radial thickness is 0.8 mm. However, in the slow cook-off condition, the ignition center position is shifted along the axial direction toward the right end face of the propellant with the increase of heating rate. Therefore, the influence of heating rate on ignition temperature and ignition delay time is nonnegligible within a certain range.展开更多
Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and opt...Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and optimized remeshing function, especially for the large deformation. In order to use this function sufficiently, simulation of cold rolling involutes spline can be implemented indirectly. The relationship between die and workpiece, forming force and characteristic of deformation in the forming process of cold rolling involutes spline are analyzed and researched. Meanwhile, reliable proofs for the design of dies and deforming equipment are provided.展开更多
On the basis of the two dimensional finite element analysis model, the pile foundations' mechanical effect of the rigid pile composite foundation under the dynamic load was researched. Through the research, the de...On the basis of the two dimensional finite element analysis model, the pile foundations' mechanical effect of the rigid pile composite foundation under the dynamic load was researched. Through the research, the development law and deformation property of axial force of pile body, shaft resistance of pile, and cumulative settlement of pile head under vertical cyclic dynamic loads were concluded. Through the comparison and analysis of the test results of dynamic models, the test results of Poulos(1989) and cumulative settlement model of the single pile under cyclic loads were confirmed. Based on the above research, Fortran language was adopted to introduce the soil attenuation factor, the secondary development of relevant modules of ABAQUS was carried out, and the effect of soil attenuation factor on dynamic property of pile-soil was discussed further.展开更多
Based on the engineering background of a soft rock roadway in Qinan Coal Mine 82 Area,Huaibei Mining Group,three creep models with different support patterns in soft rock roadway were established by using geotechnical...Based on the engineering background of a soft rock roadway in Qinan Coal Mine 82 Area,Huaibei Mining Group,three creep models with different support patterns in soft rock roadway were established by using geotechnical software of FLAC2D.According to the calculation results of different models,the change law of mechanical properties with the time of bolt-grouting support structure was obtained.Furthermore,for the test bolt-grouting support roadway,the deformation law of surrounding rock got by underground industrial experiment and field observation accords with the creep law got by numerical calculation.The results of numerical calculation and field observation show that,compared with other supports,the creep of bolt-grouting support roadway enters the steady-state creep stage from tertiary creep stage ahead,the deformations of roof,floor and two sides are decreased greatly,the plastically deforming area in surrounding rock is reduced obviously,and the distribution ranges of maximum and minimum principal stress are shrank obviously.All those fully show that the bolt-grouting support has its remarkable advantages in controlling surrounding rock creep and improving the whole strength of surrounding rock and self-bearing capacity.展开更多
In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step ...In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.展开更多
The calculation model was established by k-ε turbulence stress which reflects the change of indoor gas leak's volume fraction,and the CFD software was used to numerically simulate the volume fraction of gas after...The calculation model was established by k-ε turbulence stress which reflects the change of indoor gas leak's volume fraction,and the CFD software was used to numerically simulate the volume fraction of gas after the gas of continuity leakage,at the same time the changes of gas leak were studied. The results show that the process of gas leakage is different with the change of conditions of indoor and outdoor. Because of the different influencing factors,when the gas leak reaches a certain stable value,the volume fraction,velocity and the explosion of regional are different in the same state indoor. In some regions the gas will explode which meets the fire even if the mean volume fraction of the gas cannot achieve the explosion limit.展开更多
Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model a...Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.展开更多
Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was ...Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was usually observed near the excavation boundaries. The formation mechanism of the “11·28” rockburst, which was a typical rockburst and occurred in a drainage tunnel under a deep burial depth, high in-situ stress state and complex geological conditions, has been difficult to explain. Realistic failure process analysis(RFPA3D) software was adopted to numerically simulate the whole failure process of the surrounding rock mass around the tunnel subjected to excavation. The spatial distribution of acoustic emission derived from numerical simulation contributed to explaining the mechanical responses of the process. Analyses of the stress, safety reserve coefficient and damage degree were performed to reveal the effect of faults on the formation of rockbursts in the deep tunnel. The existence of faults results in the formation of stress anomaly areas between the tunnel and the fault. The surrounding rock mass failure propagates toward the fault from the initial failure, to different degrees. The relative positions and angles of faults play significant roles in the extent and development of surrounding rock mass failure, respectively. The increase in the lateral stress coefficient leads to the aggravation of the surrounding rock mass damage, especially in the roof and floor of the tunnel. Moreover, as the rock strength-stress ratio increases, the failure mode of the near-fault tunnel gradually changes from the stress-controlled type to the compound-controlled type. These findings were consistent with the microseismic monitoring results and field observations, which was helpful to understand the mechanical behavior of tunnel excavation affected by faults. The achievements of this study can provide some references for analysis of the failure mechanisms of similar deep tunnels.展开更多
High temperature air combustion (HTAC) is an attractive technology of saving energy and controlling environment. The mathematical models of turbulent jet flame under the highly preheated air combustion condition are c...High temperature air combustion (HTAC) is an attractive technology of saving energy and controlling environment. The mathematical models of turbulent jet flame under the highly preheated air combustion condition are conducted in the paper. The mixture fraction/probability density function model is employed. The results show that the maximum flame temperature is decreased, the temperature in the HTAC furnace is more uniform than that in the conventional furnace, and the NO x emission is low. The numerical results are partially validated by some experimental measurements.展开更多
文摘In order to study the effect of rotation on the combustion in the underwater vehicle,a two-phase turbulent combustion process is described with Reynolds stress turbulence model,eddy-dissipation turbulent combustion model,P-1 radiation model and particle tracking model of liquid. The flow in the rotating combustion chamber is simulated at two different working speeds,0?r/min and 1?000?r/min by Fluent software. The temperature,gas velocity,static pressure of wall and fuel concentration are computed and compared. The results show that the combustion in rotating combustor is faster and more effective.
文摘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.
基金Project(2025ZD1009704)supported by the National Science and Technology Major Project of ChinaProjects(2023JJ30659,2022JJ30706)supported by Hunan Provincial Natural Science Foundation,China。
文摘With the evolution of geophysical surveys from traditional two-dimensional(2 D)to three-dimensional(3 D)models,the resulting large data volumes pose significant challenges to inversion,particularly when resolving large-scale 3 D structures.A direct solver for solving an ill-conditioned linear system resulting from the finite-difference approximation of a boundary value problem requires more memory and time than iterative solvers.To overcome this limitation,an efficient iterative solver for 3 D finite-difference approach is introduced to calculate the 3 D gravitational potential and the associated gravitational field.Firstly,the boundary value problem associated with 3 D gravitational potential is discretized using central finite-difference technique based on right rectangular prismatic grids.The resulting large unsymmetric sparse systems are then solved using the generalized minimal residual algorithm(GMRES)iterative solver in combination with incomplete LU factorization.Secondly,to obtain high-accuracy partial derivatives of gravitational potential,a high-degree Lagrange interpolation scheme is employed.Finally,three density models are applied to test the accuracy,reliability,and flexibility of our 3 D finite-difference algorithm.All computational results demonstrate that our method provides an accurate approximation of the gravitational field and is applicable to 3 D forward modeling.
基金Project(51038004) supported by the National Natural Science Foundation of China
文摘In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.
基金supported by the Guangdong Major Project of Basic and Applied Basic Research(Grant No.2021B0301030001)the National Key Research and Development Program of China(Grant No.2021YFB3802300)the Foundation of National Key Laboratory of Shock Wave and Detonation Physics(Grant No.JCKYS2022212004)。
文摘The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to achieving controllable stress-strain rate loading.In this study,we have,for the first time,combined one-dimensional fluid computational software with machine learning methods.We first elucidated the mechanisms by which GDI structures control stress and strain rates.Subsequently,we constructed a machine learning model to create a structure-property response surface.The results show that altering the loading velocity and interlayer thickness has a pronounced regulatory effect on stress and strain rates.In contrast,the impedance distribution index and target thickness have less significant effects on stress regulation,although there is a matching relationship between target thickness and interlayer thickness.Compared with traditional design methods,the machine learning approach offers a10^(4)—10^(5)times increase in efficiency and the potential to achieve a global optimum,holding promise for guiding the design of GDI.
基金Projects(U23A2060,42177143)supported by the National Natural Science Foundation of China。
文摘To elucidate the influence of confining pressure on microcrack evolution and macroscopic failure mechanisms in granite,a multi-perspective approach was adopted.This approach combined triaxial compression tests,acoustic emission(AE)monitoring,and PFC simulations.The results show that:1)Confining pressure exhibits a pronounced linear correlation with both yield strength and compressive strength.The enhancement of confining pressure significantly improves the deformability of granite and promotes a progressive shift in failure mechanism from brittle rupture to ductile deformation;2)Increasing confining pressure elevates the stress threshold for microcrack initiation and suppresses crack propagation.As a result,the proportion of shear cracks increases(based on AE analysis)from 18.71%to 61.2%,marking a transition in the dominant failure mode from tensile to shear;3)Confining pressure facilitates the development of grain boundary shear cracks(GBSCs),establishing the primary damage pathways.In addition,local stress concentrations under high confinement conditions trigger intragranular cracking.This highlights the regulatory effect of confining pressure on microcrack evolution.
基金Projects(51374034,51674012)supported by the National Natural Science Foundation of ChinaProject(2013BAB02B05)supported by the China National Science and Technology Support Program during the 12th Five-Year Plan Period
文摘Proper room and pillar sizes are both critical factors for safe mining and high ore recovery rate in shrinkage stoping mining of underground metal mines.The rock masses of Tangdan copper mine of China are fractured,which needs much reinforcement and support prior to mining.Cement-sodium silicate grout technology was selected,then its related parameters such as grout pressure,diffusion radius and time were calculated and proposed.In order to test the effect of the pressured grout in the fractured No.4 ore block,field experiments were conducted.To optimize stoping configuration,three-dimensional numerical simulation with ANSYS and FLAC 3 D softwares was proposed.The results show that the drilling porosity and mechanical properties of the rock masses are increased obviously.After grout,ore recovery rate is increased by 10.2%employing the newly designed stoping configuration compared with the previous.Last,analyzed from the surface movements,roof subsidence and the maximum principal stress of the pillars,the mining safety is probable of being ensured.
基金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.
基金Open Fund of State Key Laboratory of Geological Processes and Mineral Resources, China University of Geo-sciences (GPMR0750)National Natural Science Foundation of China (40437018)
文摘We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.
文摘In recent twenty years, much numerical simulation work has been done on the evolution of Qinghai-Xizang (Tibetan) plateau. In this paper some principal numerical models and results are reviewed and analyzed. The earlier plane stress or plane strain model has much discrepancy with the actual deformation of Qinghai-Xizang plateau, such as the thickening of Tibetan crust and the lateral extrusion of Tibet along strike-slip faults. The thin viscous sheet model and the thin-plate model may simulate the change of the crustal thickness and the deformation pro-duced by gravitational force. It is suitable for studying the large-scale and long-time deformation. The influence of faults on the deformation of Tibetan plateau should be further studied.
基金Project(2023YFC2907904)supported by the National Key Research and Development Program of ChinaProject(2024RC3032)supported by the Science and Technology Innovation Program of Hunan Province,China+2 种基金Project(2024JJ4057)supported by the Natural Science Foundation of Hunan Province,ChinaProject(YDZX2024091)supported by the Shandong Provincial Technology Innovation Guidance Plan,ChinaProject(YPML-20240502090)supported by the Key R&D Program of Yunnan Precious Metals Laboratory,China。
文摘Distillation temperature,as a pivotal thermodynamic parameter in vacuum purification of metals,governs impurity migration through volatility-stratified mechanisms.This study establishes theoretical distribution models for high-volatility(Na and Se),medium-volatility(Fe and Cu),and low-volatility(Ni and Cr)impurities,revealing dual threshold effects on impurity removal:low-to-medium temperatures(£550℃)effectively suppress volatilization,while elevated temperatures promote co-evaporation.At the optimal 550℃,tellurium purity reaches 5 N 8 with>90%yield.Spatial fractionation analysis demonstrates high-volatility impurities enriching in the upper condensation zone(X/L<0.25),whereas medium/low-volatility impurities accumulate in the lower zone and residues.Remarkably,residual impurities show significant enrichment versus raw materials-Na 4.04 times,Fe 13.3 times,Cu 53.6 times,Ni 7.17 times,and Cr 15.12 times,through formation of non-volatile compounds/solid solutions.Temperature-space coupling effects drive distinct deviation patterns:fluctuations in high-volatility impurities vs.temperature-progressive deviations in medium/low-volatility species.The impurity concentration of model-experiment deviations(mean±SD)are quantified as:Se 0.265±0.12,Na 0.224±0.02,Cu 0.146±0.06,Fe 0.133±0.13,Cr 0.101±0.07,and Ni 0.122±0.08.
文摘The present study is devoted to researching the thermal security problems of large-scale solid rocket motor with Ammonium perchlorate/Hydroxyl-terminated polybutadiene(AP/HTPB). A two-dimensional axisymmetric model for the cook-off of solid rocket motor is established. The reaction kinetics for the cook-off process of AP/HTPB is described by the two-step global chemical mechanism. Numerical predictions of the cook-off behavior for the propellant are conducted at fast heating rate of 1.45-2.45 K/s,and slow heating rate of 0.001-0.003 K/s, respectively. The results show that in the fast cook-off condition. the initial ignition position of AP/HTPB occurs in the annular region of the outer wall of propellant without exception, and the center point in the region is(889.1,149.5). For the region, the axial width is1.8 mm and radial thickness is 0.8 mm. However, in the slow cook-off condition, the ignition center position is shifted along the axial direction toward the right end face of the propellant with the increase of heating rate. Therefore, the influence of heating rate on ignition temperature and ignition delay time is nonnegligible within a certain range.
文摘Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and optimized remeshing function, especially for the large deformation. In order to use this function sufficiently, simulation of cold rolling involutes spline can be implemented indirectly. The relationship between die and workpiece, forming force and characteristic of deformation in the forming process of cold rolling involutes spline are analyzed and researched. Meanwhile, reliable proofs for the design of dies and deforming equipment are provided.
基金Projects(51478178,51508181) supported by the National Natural Science Foundation of China
文摘On the basis of the two dimensional finite element analysis model, the pile foundations' mechanical effect of the rigid pile composite foundation under the dynamic load was researched. Through the research, the development law and deformation property of axial force of pile body, shaft resistance of pile, and cumulative settlement of pile head under vertical cyclic dynamic loads were concluded. Through the comparison and analysis of the test results of dynamic models, the test results of Poulos(1989) and cumulative settlement model of the single pile under cyclic loads were confirmed. Based on the above research, Fortran language was adopted to introduce the soil attenuation factor, the secondary development of relevant modules of ABAQUS was carried out, and the effect of soil attenuation factor on dynamic property of pile-soil was discussed further.
基金Project(106084) supported by the Scientific and Technological Research of the Ministry of EducationProject(BK2007701) supported by the Natural Science Foundation of Jiangsu Province,ChinaProject(2006CB2022010) supported by the National Basic Research Program of China and the Qing-lan Project of Jiangsu Province,China
文摘Based on the engineering background of a soft rock roadway in Qinan Coal Mine 82 Area,Huaibei Mining Group,three creep models with different support patterns in soft rock roadway were established by using geotechnical software of FLAC2D.According to the calculation results of different models,the change law of mechanical properties with the time of bolt-grouting support structure was obtained.Furthermore,for the test bolt-grouting support roadway,the deformation law of surrounding rock got by underground industrial experiment and field observation accords with the creep law got by numerical calculation.The results of numerical calculation and field observation show that,compared with other supports,the creep of bolt-grouting support roadway enters the steady-state creep stage from tertiary creep stage ahead,the deformations of roof,floor and two sides are decreased greatly,the plastically deforming area in surrounding rock is reduced obviously,and the distribution ranges of maximum and minimum principal stress are shrank obviously.All those fully show that the bolt-grouting support has its remarkable advantages in controlling surrounding rock creep and improving the whole strength of surrounding rock and self-bearing capacity.
基金Project(41630642)supported by the Key Project of National Natural Science Foundation of ChinaProject(51974360)supported by the National Natural Science Foundation of ChinaProject(2018JJ3656)supported by the Natural Science Foundation of Hunan Province,China。
文摘In the context of deep rock engineering,the in-situ stress state is of major importance as it plays an important role in rock dynamic response behavior.Thus,stress initialization becomes crucial and is the first step for the dynamic response simulation of rock mass in a high in-situ stress field.In this paper,stress initialization methods,including their principles and operating procedures for reproducing steady in-situ stress state in LS-DYNA,are first introduced.Then the most popular four methods,i.e.,explicit dynamic relaxation(DR)method,implicit-explicit sequence method,Dynain file method and quasi-static method,are exemplified through a case analysis by using the RHT and plastic hardening rock material models to simulate rock blasting under in-situ stress condition.Based on the simulations,it is concluded that the stress initialization results obtained by implicit-explicit sequence method and dynain file method are closely related to the rock material model,and the explicit DR method has an obvious advantage in solution time when compared to other methods.Besides that,it is recommended to adopt two separate analyses for the whole numerical simulation of rock mass under the combined action of in-situ stress and dynamic disturbance.
基金Projects(50704041, 50838009, 50678179) supported by the National Natural Science Foundation of ChinaProjects(2006BAJ02A09, 2006BAJ02A13-4) supported by the National Key Technologies R & D Program of China+1 种基金Project(CSTC2008BB6347) supported by the Natural Science Foundation Project of ChongqingProject(20080430752) supported by the Postdoctoral Science Foundation of China
文摘The calculation model was established by k-ε turbulence stress which reflects the change of indoor gas leak's volume fraction,and the CFD software was used to numerically simulate the volume fraction of gas after the gas of continuity leakage,at the same time the changes of gas leak were studied. The results show that the process of gas leakage is different with the change of conditions of indoor and outdoor. Because of the different influencing factors,when the gas leak reaches a certain stable value,the volume fraction,velocity and the explosion of regional are different in the same state indoor. In some regions the gas will explode which meets the fire even if the mean volume fraction of the gas cannot achieve the explosion limit.
基金Project supported by the Scientific Research Foundation of Hebei University of Technology of China
文摘Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.
基金Project(42177143) supported by the National Natural Science Foundation of ChinaProject(2020JDJQ0011) supported by the Science Foundation for Distinguished Young Scholars of Sichuan Province,China。
文摘Rockbursts were frequently encountered in the construction of deeply buried tunnels at the Jinping-II hydropower station, Southwest China. In those cases, the existence of large structural planes, such as faults, was usually observed near the excavation boundaries. The formation mechanism of the “11·28” rockburst, which was a typical rockburst and occurred in a drainage tunnel under a deep burial depth, high in-situ stress state and complex geological conditions, has been difficult to explain. Realistic failure process analysis(RFPA3D) software was adopted to numerically simulate the whole failure process of the surrounding rock mass around the tunnel subjected to excavation. The spatial distribution of acoustic emission derived from numerical simulation contributed to explaining the mechanical responses of the process. Analyses of the stress, safety reserve coefficient and damage degree were performed to reveal the effect of faults on the formation of rockbursts in the deep tunnel. The existence of faults results in the formation of stress anomaly areas between the tunnel and the fault. The surrounding rock mass failure propagates toward the fault from the initial failure, to different degrees. The relative positions and angles of faults play significant roles in the extent and development of surrounding rock mass failure, respectively. The increase in the lateral stress coefficient leads to the aggravation of the surrounding rock mass damage, especially in the roof and floor of the tunnel. Moreover, as the rock strength-stress ratio increases, the failure mode of the near-fault tunnel gradually changes from the stress-controlled type to the compound-controlled type. These findings were consistent with the microseismic monitoring results and field observations, which was helpful to understand the mechanical behavior of tunnel excavation affected by faults. The achievements of this study can provide some references for analysis of the failure mechanisms of similar deep tunnels.
文摘High temperature air combustion (HTAC) is an attractive technology of saving energy and controlling environment. The mathematical models of turbulent jet flame under the highly preheated air combustion condition are conducted in the paper. The mixture fraction/probability density function model is employed. The results show that the maximum flame temperature is decreased, the temperature in the HTAC furnace is more uniform than that in the conventional furnace, and the NO x emission is low. The numerical results are partially validated by some experimental measurements.
