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 study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyur...To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.展开更多
In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommend...In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.展开更多
The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distanc...The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distance flight airframe power.In this paper,a computational model of the scramjet magnetohydrody⁃namic channel is developed and verified by using the commercial software Fluent.It is found that when the mag⁃netic induction intensity is 1,2,3,4 T,the power generation efficiency is 22.5%,22.3%,22.0%,21.5%,and decreases with the increase of the magnetic induction intensity,and the enthalpy extraction rate is 0.026%,0.1%,0.21%,0.34%,and increases with the increase of the magnetic induction intensity.The deceleration ef⁃fect of electromagnetic action on the airflow in the power channel increases with the increase of magnetic induc⁃tion intensity.The stronger the magnetic field intensity,the more obvious the decreasing effect of fluid Mach num⁃ber in the channel.The power generation efficiency decreases as the magnetic induction intensity increases and the enthalpy extraction rate is reversed.As the local currents gathering at inlet and outlet of the power generation area,total temperature and enthalpy along the flow direction do not vary linearly,and there are maximum and minimum values at inlet and outlet.Increasing the number of electrodes can effectively regulate the percentage of Joule heat dissipation,which can improve the power generation efficiency.展开更多
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.展开更多
Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of...Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.展开更多
In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading ef...In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading effect of roof vibration on roadway surrounding rocks as well as the impact on stability. The results show that, on one hand, hard-thick roof will result in high stress concentration on mining surrounding rocks; on the other hand, the breaking of hard-thick roof will lead to mining seismicity, causing dynamic loading effect on coal and rock mass. High stress concentration and dynamic loading combination reaches to the mechanical conditions for the occurrence of rock burst, which will induce rock burst. The mining induced seismic events occurring in the roof breaking act on the mining surrounding rocks in the form of stress wave. The stress wave then has a reflection on the free surface of roadway and the tensile stress will be generated around the free surface. Horizontal vibration of roadway surrounding particles will cause instant changes of horizontal stress of roadway surrounding rocks; the horizontal displacement is directly related to the horizontal stress but is not significantly correlated with the vertical stress; the increase of horizontal stress of roadway near surface surrounding rocks and the release of elastic deformation energy of deep surrounding coal and rock mass are immanent causes that lead to the impact instability of roadway surrounding rocks. The most significant measures for rock burst prevention are controlling of horizontal stress and vibration strength.Key words展开更多
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.展开更多
To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mat...To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.展开更多
In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by usi...In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by using FLUENT software,and the numerical simulation analysis of tunnel drilling ventilation and dust removal parameters was carried out.The results show that it is difficult to control the dust pollution of the face by conventional ventilation,and the drilling dust is distributed in the range of 10 m from the face;after the introduction of the long pressure and short suction ventilation scheme,when the ratio of compressed air volume to exhaust air volume is 0.72,the height of the pressure fan is 2.5 m,the distance between the pressure fan and the palm face is 20 m,and the exhaust fan is 12 m away from the palm,the dust concentration control efficiency of the working face is increased by about 60%.Therefore,in the similar long-distance single head tunnel construction,it is appropriate to adopt the dust removal method of long-distance short suction and exhaust fan to ensure the working environment.展开更多
Coal mine belt fire develops very rapidly and is difficult to control. If not suppressed quickly, a belt fire could easily lead to airflow disorder and undermine the ventilation system. However, belt fire can be preve...Coal mine belt fire develops very rapidly and is difficult to control. If not suppressed quickly, a belt fire could easily lead to airflow disorder and undermine the ventilation system. However, belt fire can be prevented effectively by establishing fire airflow control system. In this work, the 5th belt roadway of Kongzhuang coal mine was taken as the object of investigation, where geometrical models of this roadway were established firstly. Then, based on mathematical model of fire smoke flow, the CO volume fraction, smoke density distribution, air temperature and pollutant velocity vector in the roadway before and after taking airflow control measures were simulated by using Fluent software. It can be known from the simulation that with the normal ventilation status in 5th belt roadway, the countercurrent of smoke does not happen when a fire occurs; the roadway's section is almost filled with CO at 10 m downstream from the fire source, and with air velocity getting stable gradually, the CO concentration reaches about 15 %. After taking airflow control measures, the effect range of temperature field which are harmful to the miners decreases from 69 m to 30 m; and the distance of the roadway fully filled with CO is 5 m farther than that before taking measures. Finally, according to the numerical simulation results and the actual condition of the belt roadway, the warning and automatic remote airflow control system with short-circuit method for the 5th belt roadway was designed to guarantee the safety production.展开更多
During the compaction of a road subgrade, the mechanical parameters of the soil mass change in real time, but current research assumes that these parameters remain unchanged. In order to address this discrepancy, this...During the compaction of a road subgrade, the mechanical parameters of the soil mass change in real time, but current research assumes that these parameters remain unchanged. In order to address this discrepancy, this paper establishes a relationship between the degree of compaction K and strain ε. The relationship between the compaction degree K and the shear strength of soil(cohesion c and frictional angle φ) was clearly established through indoor experiments. The subroutine UMAT in ABAQUS finite element numerical software was developed to realize an accurate calculation of the subgrade soil compaction quality. This value was compared and analyzed against the assumed compaction value of the model, thereby verifying the accuracy of the intelligent compaction calculation results for subgrade soil. On this basis, orthogonal tests of the influential factors(frequency, amplitude, and quality) for the degree of compaction and sensitivity analysis were carried out. Finally, the ‘acceleration intelligent compaction value’, which is based on the acceleration signal, is proposed for a compaction meter value that indicates poor accuracy. The research results can provide guidance and basis for further research into the accurate control of compaction quality for roadbeds and pavements.展开更多
Frac-packing technology has been introduced to improve the development effect of weakly consolidated sandstone.It has double effects on increasing production and sand control.However,determining operation parameters o...Frac-packing technology has been introduced to improve the development effect of weakly consolidated sandstone.It has double effects on increasing production and sand control.However,determining operation parameters of frac-packing is the key factor due to the particularity of weakly consolidated sandstone.In order to study the mechanisms of hydraulic fracture propagation and reveal the effect of fracturing parameters on fracture morphology in weakly consolidated sandstone,finite element numerical model of fluid-solid coupling is established to carry out numerical simulation to analyze influences of mechanical characteristics,formation permeability,fracturing fluid injection rate and viscosity on fracture propagation.The result shows that lower elastic modulus is favorable for inducing short and wide fractures and controls the fracture length while Poisson ratio has almost no effect.Large injection rate and high viscosity of fracturing fluid are advantageous to fracture initiation and propagation.Suitable fractures are produced when the injection rate is approximate to3–4m3/min and fluid viscosity is over100mPa?s.The leak-off of fracturing fluid to formation is rising with the increase of formation permeability,which is adverse to fracture propagation.The work provides theoretical reference to determine the construction parameters for the frac-packing design in weakly consolidated reservoirs.展开更多
Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out...Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out area was carried out with FLAC3D. The whole analytical process was divided into two steps, including the static analysis and the dynamic analysis which were used to simulate the influence of excavation process and blasting vibration respectively. The results show that the shape of right upper boundary is extremely irregular after excavation, and stress concentration occurs at many places and higher tensile stress appears. The maximum tensile stress is higher than the tensile strength of rock mass, and surrounding rock of right roof will be damaged with tension fracture. The maximum displacement of surrounding rock is 4.75 mm after excavation. However, the maximum displacement increases to 5.47 mm after the blasting dynamic load is applied. And the covering area of plastic zones expands obviously, especially at the foot of right upper slope. The analytical results are in basic accordance with the observed results on the whole. Damage and disturbance on surrounding rock to some degree are caused by excavation, while blasting dynamic load increases the possibility of occurrence of dynamic instability and destruction further. So the effective supporting and vibration reducing measures should be taken during mining.展开更多
Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated...Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method （FEM）. The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/（m2. K）. When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the： decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.展开更多
Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an...Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an orthogonal design with two factors, three levels and nine runs was proposed, followed by three-dimensional numerical simulation using ANSYS and FLAC3~. Results show that surface settlement after excavation is concentrically ringed, and increases with the decrease of pillar width and distances to stope gobs. In the meantime, the ore-control fault at the ore-rock boundary and the fractured argillaceous dolomite with intercalated slate at the hanging wall deteriorate the roof settlement. Additionally, stope stability is challenged due to pillar rheological yield and stress concentration, and both are induced by redistribution of stress and plastic zones after mining. Following an objective function and a constraint function, room and pillar configuration with widths of 14 m and 16 m, respectively, is presented as the optimization for improving the ore recovery rate while maintaining a safe working environment.展开更多
Due to its high strength,high density,high hardness and good penetration capabilities,Depleted uranium alloys have already shined in armor-piercing projectiles.There should also be a lot of room for improvement in the...Due to its high strength,high density,high hardness and good penetration capabilities,Depleted uranium alloys have already shined in armor-piercing projectiles.There should also be a lot of room for improvement in the application of fragment killing elements.Therefore,regarding the performance of the depleted uranium alloy to penetrate the target plate,further investigation is needed to analyze its advantages and disadvantages compared to tungsten alloy.To study the difference in penetration performance between depleted uranium alloy and tungsten alloy fragments,firstly,a theoretical analysis of the adiabatic shear sensitivity of DU and tungsten alloys was given from the perspective of material constitutive model.Then,taking the cylindrical fragment penetration target as the research object,the penetration process and velocity characteristics of the steel target plates penetrated by DU alloy fragment and tungsten alloy fragment were compared and analyzed,by using finite element software ANSYS/LS-DYNA and Lagrange algorithm.Lastly,the influence of different postures when impacting target and different fragment shapes on the penetration results is carried out in the research.The results show that in the penetration process of the DU and tungsten alloy fragments,the self-sharpening properties of the DU alloy can make the fragment head sharper and the penetrating ability enhance.Under the same conditions,the penetration capability of cylindrical fragment impacting target in vertical posture is better than that in horizontal posture,and the penetration capability of the spherical fragment is slightly better than that of cylindrical fragment.展开更多
The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) appr...The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.展开更多
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.展开更多
When underground cavities are subjected to explosive stress waves,a uniquely damaged zone may appear due to the combined effect of dynamic loading and static pre-load stress.In this study,a rate-dependent two-dimensio...When underground cavities are subjected to explosive stress waves,a uniquely damaged zone may appear due to the combined effect of dynamic loading and static pre-load stress.In this study,a rate-dependent two-dimensional rock dynamic constitutive model was established to investigate the dynamic fractures of rocks under different static stress conditions.The effects of the loading rate and peak amplitude of the blasting wave under different confining pressures and the vertical compressive coefficient(K_(0))were considered.The numerical simulated results reproduced the initiation and further propagation of primary radial crack fractures,which were in agreement with the experimental results.The dynamic loading rate,peak amplitude,static vertical compressive coefficient(K_(0))and confining pressure affected the evolution of fractures around the borehole.The heterogeneity parameter(m)plays an important role in the evolution of fractures around the borehole.The crack propagation path became more discontinuous and rougher in a smallerheterogeneity parameter case.展开更多
基金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.
基金the Natural Science Foundation of Jiangsu Province(Grant No.BK20200494)China Postdoctoral Science Foundation(Grant No.2021M701725)+3 种基金Jiangsu Postdoctoral Research Funding Program(Grant No.2021K522C)Fundamental Research Funds for the Central Universities(Grant No.30919011246)National Natural Science Foundation of China(Grant No.52278188)Natural Science Foundation of Jiangsu Province(Grant No.BK20211196)。
文摘To study the anti-explosion protection effect of polyurea coating on reinforced concrete box girder,two segmental girder specimens were made at a scale of 1:3,numbered as G(without polyurea coating)and PCG(with polyurea coating).The failure characteristics and dynamic responses of the specimens were compared through conducting explosion tests.The reliability of the numerical simulation using LS-DYNA software was verified by the test results.The effects of different scaled distances,reinforcement ratios,concrete strengths,coating thicknesses and ranges of polyurea were studied.The results show that the polyurea coating can effectively enhance the anti-explosion performance of the girder.The top plate of middle chamber in specimen G forms an elliptical penetrating hole,while that in specimen PCG only shows a very slight local dent.The peak vertical displacement and residual displacement of PCG decrease by 74.8% and 73.7%,respectively,compared with those of specimen G.For the TNT explosion with small equivalent,the polyurea coating has a more significant protective effect on reducing the size of fracture.With the increase of TNT equivalent,the protective effect of polyurea on reducing girder displacement becomes more significant.The optimal reinforcement ratio,concrete strength,thickness and range of polyurea coating were also drawn.
基金the funding supports of the National Key Research and Development Plan,China(Grant No.2022YFC3801800)National Natural Science Foundation of China(Grant Nos.52038010 and 52078368)。
文摘In order to improve the overall resilience of the urban infrastructures, it is required to conduct blast resistant design for important building structures in the city. For complex terrain in the city, it is recommended to determine the blast load on the structures via numerical simulation. Since the mesh size of the numerical model highly depends on the explosion scenario, there is no generally applicable approach for the mesh size selection. An efficient method to determine the mesh size of the numerical model of near-ground detonation based on explosion scenarios is proposed in this study. The effect of mesh size on the propagation of blast wave under different explosive weights was studied, and the correlations between the mesh size effect and the charge weight or the scaled distance was described. Based on the principle of the finite element method and Hopkinson-Cranz scaling law, a mesh size measurement unit related to the explosive weight was proposed as the criterion for determining the mesh size in the numerical simulation. Finally, the applicability of the method proposed in this paper was verified by comparing the results from numerical simulation and the explosion tests and was verified in AUTODYN.
文摘The reverse magnetohydrodynamic(MHD)energy bypass technology is a promising energy redis⁃tribution technology in the scramjet system,in augmented with a power generation equipment to supply the neces⁃sary long-distance flight airframe power.In this paper,a computational model of the scramjet magnetohydrody⁃namic channel is developed and verified by using the commercial software Fluent.It is found that when the mag⁃netic induction intensity is 1,2,3,4 T,the power generation efficiency is 22.5%,22.3%,22.0%,21.5%,and decreases with the increase of the magnetic induction intensity,and the enthalpy extraction rate is 0.026%,0.1%,0.21%,0.34%,and increases with the increase of the magnetic induction intensity.The deceleration ef⁃fect of electromagnetic action on the airflow in the power channel increases with the increase of magnetic induc⁃tion intensity.The stronger the magnetic field intensity,the more obvious the decreasing effect of fluid Mach num⁃ber in the channel.The power generation efficiency decreases as the magnetic induction intensity increases and the enthalpy extraction rate is reversed.As the local currents gathering at inlet and outlet of the power generation area,total temperature and enthalpy along the flow direction do not vary linearly,and there are maximum and minimum values at inlet and outlet.Increasing the number of electrodes can effectively regulate the percentage of Joule heat dissipation,which can improve the power generation efficiency.
基金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.
基金Project(51375498)supported by the National Natural Science Foundation of China
文摘Based on RNG k-ε turbulence model and sliding grid technique, solid-liquid two-phase three-dimensional(3-D) unsteady turbulence of full passage in slurry pump was simulated by means of Fluent software. The effects of unsteady flow characteristics on solid-liquid two-phase flow and pump performance were researched under design condition. The results show that clocking effect has a significant influence on the flow in pump, and the fluctuation of flow velocity and pressure is obvious, particularly near the volute tongue, at the position of small sections of volute and within diffuser. Clocking effect has a more influence on liquid-phase than on solid-phase, and the wake-jet structure of relative velocity of solid-phase is less obvious than liquid-phase near the volute tongue and the impeller passage outlet. The fluctuation of relative velocity of solid-phase flow is 7.6% smaller than liquid-phase flow at the impeller outlet on circular path. Head and radial forces of the impeller are 8.1% and 85.7% of fluctuation, respectively. The results provide a theoretical basis for further research for turbulence, improving efficient, reducing the hydraulic losses and wear. Finally, field tests were carried out to verify the operation and wear of slurry pump.
基金Project(51404243)supported by the National Natural Science Foundation of ChinaProject(2014QNB26)supported by the Fundamental Research Funds for the Central Universities,China
文摘In order to reveal the dynamic process of hard-thick roof inducing rock burst, one of the most common and strongest dynamic disasters in coal mine, the numerical simulation is conducted to study the dynamic loading effect of roof vibration on roadway surrounding rocks as well as the impact on stability. The results show that, on one hand, hard-thick roof will result in high stress concentration on mining surrounding rocks; on the other hand, the breaking of hard-thick roof will lead to mining seismicity, causing dynamic loading effect on coal and rock mass. High stress concentration and dynamic loading combination reaches to the mechanical conditions for the occurrence of rock burst, which will induce rock burst. The mining induced seismic events occurring in the roof breaking act on the mining surrounding rocks in the form of stress wave. The stress wave then has a reflection on the free surface of roadway and the tensile stress will be generated around the free surface. Horizontal vibration of roadway surrounding particles will cause instant changes of horizontal stress of roadway surrounding rocks; the horizontal displacement is directly related to the horizontal stress but is not significantly correlated with the vertical stress; the increase of horizontal stress of roadway near surface surrounding rocks and the release of elastic deformation energy of deep surrounding coal and rock mass are immanent causes that lead to the impact instability of roadway surrounding rocks. The most significant measures for rock burst prevention are controlling of horizontal stress and vibration strength.Key words
基金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.
基金Project(2006AA03Z523) supported by the National High-Tech Research and Development Program of ChinaProject(08C26224302178) supported by the Innovation Foundation of Central South University,China
文摘To predict three-dimensional temperature distribution of molten aluminum and its influencing factors inside an industrial aluminum holding furnace,a fluid-solid coupled method was presented.The fluid-solid coupled mathematics models of aluminum holding furnace in the premixed combustion processing were established based on mass conservation,moment conservation,momentum conservation,energy conservation and chemistry species conservation.Computational results agree well with the test data of the typical condition.The maximum combustion temperature is 1 850 K.The average temperature of the molten aluminum is 1 158 K,and the maximum temperature difference is about 240 K.The average temperature increases 0.3 ℃ while the temperature of combustion air increases 1 ℃.The optimal excess air ratio is 1.25-1.30.
基金Project(51874016)supported by the National Natural Science Foundation of China。
文摘In order to control the dust pollution produced by air leg rock drill in the trolley area during the excavation of long-distance single ended tunnel,the full-scale physical model of working face was established by using FLUENT software,and the numerical simulation analysis of tunnel drilling ventilation and dust removal parameters was carried out.The results show that it is difficult to control the dust pollution of the face by conventional ventilation,and the drilling dust is distributed in the range of 10 m from the face;after the introduction of the long pressure and short suction ventilation scheme,when the ratio of compressed air volume to exhaust air volume is 0.72,the height of the pressure fan is 2.5 m,the distance between the pressure fan and the palm face is 20 m,and the exhaust fan is 12 m away from the palm,the dust concentration control efficiency of the working face is increased by about 60%.Therefore,in the similar long-distance single head tunnel construction,it is appropriate to adopt the dust removal method of long-distance short suction and exhaust fan to ensure the working environment.
基金Project supported by Joint Innovative Center for Safe and Effective Mining Technology and Equipment of Coal Resources of Shandong Province,ChinaProject supported by Taishan Scholar Program of Shandong Province,China+2 种基金Project(2014JQJH106)supported by Shandong University of Science and Technology Research Fund,ChinaProject(LAK2012-1)supported by Science and Technology Development Program of Safety Production of Shandong Province,ChinaProject(2012ZHTD06)supported by Science Research Innovative Group of College of Mining and Safety Engineering of Shandong University of Science and Technology,China
文摘Coal mine belt fire develops very rapidly and is difficult to control. If not suppressed quickly, a belt fire could easily lead to airflow disorder and undermine the ventilation system. However, belt fire can be prevented effectively by establishing fire airflow control system. In this work, the 5th belt roadway of Kongzhuang coal mine was taken as the object of investigation, where geometrical models of this roadway were established firstly. Then, based on mathematical model of fire smoke flow, the CO volume fraction, smoke density distribution, air temperature and pollutant velocity vector in the roadway before and after taking airflow control measures were simulated by using Fluent software. It can be known from the simulation that with the normal ventilation status in 5th belt roadway, the countercurrent of smoke does not happen when a fire occurs; the roadway's section is almost filled with CO at 10 m downstream from the fire source, and with air velocity getting stable gradually, the CO concentration reaches about 15 %. After taking airflow control measures, the effect range of temperature field which are harmful to the miners decreases from 69 m to 30 m; and the distance of the roadway fully filled with CO is 5 m farther than that before taking measures. Finally, according to the numerical simulation results and the actual condition of the belt roadway, the warning and automatic remote airflow control system with short-circuit method for the 5th belt roadway was designed to guarantee the safety production.
基金Project(51878164) supported by the National Natural Science Foundation of ChinaProjects(BK20161421, BK20140109) supported by the Natural Science Foundation of Jiangsu Province, China+4 种基金Project(141076) supported by the Huoyingdong Foundation of the Ministry of Education of ChinaProject(BZ2017011) supported by the Science and Technology Support Project of Jiangsu Province, ChinaProject(2242015R30027) supported by the Fundamental Research Funds for the Central Universities, ChinaProject(grant number KFJ170106) supported by the Changsha University of Science & Technology via Open Fund of National Engineering Laboratory of Highway Maintenance Technology, ChinaProject(2018B51) supported by the Science and Technology Support Project of Qilu Transportation Development Group, China。
文摘During the compaction of a road subgrade, the mechanical parameters of the soil mass change in real time, but current research assumes that these parameters remain unchanged. In order to address this discrepancy, this paper establishes a relationship between the degree of compaction K and strain ε. The relationship between the compaction degree K and the shear strength of soil(cohesion c and frictional angle φ) was clearly established through indoor experiments. The subroutine UMAT in ABAQUS finite element numerical software was developed to realize an accurate calculation of the subgrade soil compaction quality. This value was compared and analyzed against the assumed compaction value of the model, thereby verifying the accuracy of the intelligent compaction calculation results for subgrade soil. On this basis, orthogonal tests of the influential factors(frequency, amplitude, and quality) for the degree of compaction and sensitivity analysis were carried out. Finally, the ‘acceleration intelligent compaction value’, which is based on the acceleration signal, is proposed for a compaction meter value that indicates poor accuracy. The research results can provide guidance and basis for further research into the accurate control of compaction quality for roadbeds and pavements.
基金Project(2016ZX05058-002-006)supported by National Science and Technology Major Projects of ChinaProject(2018CXTD346)supported by Innovative Research Team Program of Natural Science Foundation of Hainan Province,China
文摘Frac-packing technology has been introduced to improve the development effect of weakly consolidated sandstone.It has double effects on increasing production and sand control.However,determining operation parameters of frac-packing is the key factor due to the particularity of weakly consolidated sandstone.In order to study the mechanisms of hydraulic fracture propagation and reveal the effect of fracturing parameters on fracture morphology in weakly consolidated sandstone,finite element numerical model of fluid-solid coupling is established to carry out numerical simulation to analyze influences of mechanical characteristics,formation permeability,fracturing fluid injection rate and viscosity on fracture propagation.The result shows that lower elastic modulus is favorable for inducing short and wide fractures and controls the fracture length while Poisson ratio has almost no effect.Large injection rate and high viscosity of fracturing fluid are advantageous to fracture initiation and propagation.Suitable fractures are produced when the injection rate is approximate to3–4m3/min and fluid viscosity is over100mPa?s.The leak-off of fracturing fluid to formation is rising with the increase of formation permeability,which is adverse to fracture propagation.The work provides theoretical reference to determine the construction parameters for the frac-packing design in weakly consolidated reservoirs.
基金Project (50490272) supported by the National Natural Science Foundation of China project(NCET-05-0687) supportedby Programfor New Century Excellent Talents project (040109) supported bythe Doctor Degree Paper Innovation Engineering of CentralSouth University
文摘Dynamic analysis steps and general flow of fast lagrangian analysis of continua in 3 dimensions(FLAC3D) were discussed. Numerical simulation for influence of excavation and blasting vibration on stability of mined-out area was carried out with FLAC3D. The whole analytical process was divided into two steps, including the static analysis and the dynamic analysis which were used to simulate the influence of excavation process and blasting vibration respectively. The results show that the shape of right upper boundary is extremely irregular after excavation, and stress concentration occurs at many places and higher tensile stress appears. The maximum tensile stress is higher than the tensile strength of rock mass, and surrounding rock of right roof will be damaged with tension fracture. The maximum displacement of surrounding rock is 4.75 mm after excavation. However, the maximum displacement increases to 5.47 mm after the blasting dynamic load is applied. And the covering area of plastic zones expands obviously, especially at the foot of right upper slope. The analytical results are in basic accordance with the observed results on the whole. Damage and disturbance on surrounding rock to some degree are caused by excavation, while blasting dynamic load increases the possibility of occurrence of dynamic instability and destruction further. So the effective supporting and vibration reducing measures should be taken during mining.
基金Project(20080431380) supported by China Postdoctoral Science Foundation
文摘Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method （FEM）. The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/（m2. K）. When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the： decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.
基金Projects(50934002,51074013,51104100)supported by the National Natural Science Foundation of ChinaProject(IRT0950)supported by the Program for Changjiang Scholars and Innovative Research Team in University of China
文摘Room and pillar sizes are key factors for safe mining and ore recovery in open-stope mining. To investigate the influence of room and pillar configurations on stope stability in highly fractured and weakened areas, an orthogonal design with two factors, three levels and nine runs was proposed, followed by three-dimensional numerical simulation using ANSYS and FLAC3~. Results show that surface settlement after excavation is concentrically ringed, and increases with the decrease of pillar width and distances to stope gobs. In the meantime, the ore-control fault at the ore-rock boundary and the fractured argillaceous dolomite with intercalated slate at the hanging wall deteriorate the roof settlement. Additionally, stope stability is challenged due to pillar rheological yield and stress concentration, and both are induced by redistribution of stress and plastic zones after mining. Following an objective function and a constraint function, room and pillar configuration with widths of 14 m and 16 m, respectively, is presented as the optimization for improving the ore recovery rate while maintaining a safe working environment.
文摘Due to its high strength,high density,high hardness and good penetration capabilities,Depleted uranium alloys have already shined in armor-piercing projectiles.There should also be a lot of room for improvement in the application of fragment killing elements.Therefore,regarding the performance of the depleted uranium alloy to penetrate the target plate,further investigation is needed to analyze its advantages and disadvantages compared to tungsten alloy.To study the difference in penetration performance between depleted uranium alloy and tungsten alloy fragments,firstly,a theoretical analysis of the adiabatic shear sensitivity of DU and tungsten alloys was given from the perspective of material constitutive model.Then,taking the cylindrical fragment penetration target as the research object,the penetration process and velocity characteristics of the steel target plates penetrated by DU alloy fragment and tungsten alloy fragment were compared and analyzed,by using finite element software ANSYS/LS-DYNA and Lagrange algorithm.Lastly,the influence of different postures when impacting target and different fragment shapes on the penetration results is carried out in the research.The results show that in the penetration process of the DU and tungsten alloy fragments,the self-sharpening properties of the DU alloy can make the fragment head sharper and the penetrating ability enhance.Under the same conditions,the penetration capability of cylindrical fragment impacting target in vertical posture is better than that in horizontal posture,and the penetration capability of the spherical fragment is slightly better than that of cylindrical fragment.
基金Project(41174061) supported by the National Natural Science Foundation of ChinaProject(2011QNZT011) supported by the Free Exploration Program of Central South University,China
文摘The absorbing boundary is the key in numerical simulation of borehole radar.Perfect match layer(PML) was chosen as the absorbing boundary in numerical simulation of GPR.But CPML(convolutional perfect match layer) approach that we have chosen has the advantage of being media independent.Beginning with the Maxwell equations in a two-dimensional structure,numerical formulas of finite-difference time-domain(FDTD) method with CPML boundary condition for transverse electric(TE) or transverse magnetic(TM) wave are presented in details.Also,there are three models for borehole-GPR simulation.By analyzing the simulation results,the features of targets in GPR are obtained,which can provide a better interpretation of real radar data.The results show that CPML is well suited for the simulation of borehole-GPR.
基金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.
基金Projects(51878190,51779031,51678170)supported by the National Natural Science Foundation of China。
文摘When underground cavities are subjected to explosive stress waves,a uniquely damaged zone may appear due to the combined effect of dynamic loading and static pre-load stress.In this study,a rate-dependent two-dimensional rock dynamic constitutive model was established to investigate the dynamic fractures of rocks under different static stress conditions.The effects of the loading rate and peak amplitude of the blasting wave under different confining pressures and the vertical compressive coefficient(K_(0))were considered.The numerical simulated results reproduced the initiation and further propagation of primary radial crack fractures,which were in agreement with the experimental results.The dynamic loading rate,peak amplitude,static vertical compressive coefficient(K_(0))and confining pressure affected the evolution of fractures around the borehole.The heterogeneity parameter(m)plays an important role in the evolution of fractures around the borehole.The crack propagation path became more discontinuous and rougher in a smallerheterogeneity parameter case.
