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.展开更多
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.展开更多
During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configura...During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configuration of the optical path within the internal channel necessitates complex and time-consuming efforts to assess the impact of thermal blooming effect on the optical path.To meet the engineering need for rapid evaluation of thermal blooming effect in optical paths,this study proposed a rapid simulation method for the thermal blooming effect in internal optical paths based on the finite element method.This method discretized the fluid region into infinitesimal elements and employed finite element method for flow field analysis.A simplified analytical model of the flow field region in complex internal channels was established,and regions with similar thermal blooming effect were divided within this model.Based on the calculated optical path differences within these regions,numerical simulations of phase distortion caused by thermal blooming were conducted.The calculated result were compared with those obtained using the existing methods.The findings reveal that for complex optical paths,the discrepancy between the two approaches is less than 3.6%,with similar phase distortion patterns observed.For L-type units,this method and the existing methods identify the same primary factors influencing aberrations and exhibit consistent trends in their variation.This method was used to analyze the impact of thermal blooming effect in a straight channel under different gravity directions.The results show that phase distortion varies with changes in the direction of gravity,and the magnitude of the phase difference is strongly correlated with the component of gravity perpendicular to the optical axis.Compared to the existing methods,this approach offers greater flexibility,obviates the need for complex custom analysis programming.The analytical results of this method enable a rapid assessment of the thermal blooming effect in optical paths within the internal channel.This is especially useful during the engineering design.These results also provide crucial references for developing strategies to suppress thermal blooming effect.展开更多
Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical cha...Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical charges on the tunnel blast wave loads and to develop a quantitative calculation method, this study carried out experimental and numerical research. Initially, external explosion experiments were conducted using both 35 kg spherical charges and cylindrical charges with aspect ratio of 4.8 at two different distances from the tunnel entrance. Comparative analysis of the blast wave parameters in the tunnel revealed that the explosive equivalent of the cylindrical charges was significantly higher than that of the spherical charges. To address this, an equivalent coefficient κ based on the spherical charges was proposed for the cylindrical charges. Subsequently, numerical simulations were conducted for the experimental conditions, and the numerical simulation results match the experiments well. Through numerical calculations, the reliability of the equivalent coefficient κ under the experimental conditions was verified, and comparison analysis indicated that the explosion energy of cylindrical charges spreads more radially, resulting in more explosion energy entering the tunnel, which is the fundamental reason for the increase in tunnel blast wave loads. Additionally, analyzing the explosion energy ratio entering the tunnel is an effective method for calculating the equivalent coefficient κ. Finally, through more than one hundred sets of numerical calculation results, the impact of the proportional distance λ and the ratio of charge mass to the tunnel cross-section dimension φ on the equivalence coefficients κ was investigated. An empirical formula for the equivalence coefficient κ was derived through fitting, and the accuracy of the formula was validated through literature experimental results. The research findings of this paper will provide valuable guidance for the calculation of blast wave loads in tunnel.展开更多
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.展开更多
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.展开更多
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展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
High-speed train running in the sand environment is different from the general environment. In the former situation, there will be sand load applied on high-speed train(SLAHT) caused by sand particles hitting train su...High-speed train running in the sand environment is different from the general environment. In the former situation, there will be sand load applied on high-speed train(SLAHT) caused by sand particles hitting train surface. This will have a great impact on the train stability, running drag and surface corrosion. Numerical simulation method of SLAHT in sand environment is studied. The velocity and mass flow rate models of saltation and suspension sand particles and the calculation model of SLAHT caused by sand particles hitting train surface are established. The discrete phase method is adopted for numerical simulating the process of saltation and suspension sand particles moving to train surface and generating sand load. By comparison with the field tests, the numerical simulation reliability is analysed. The theoretical formula of SLAHT changing with cross-wind and train speed is proposed. SLAHT changing law is analyzed. Research results indicate that SLAHT changing with cross-wind and train speed is a quadratic relationship. When train speed is constant, SLAHT increases quadratically with cross-wind speed improvement. When cross-wind speed is constant, SLAHT increases quadratically with train speed improvement.展开更多
基金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.
基金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.
文摘During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configuration of the optical path within the internal channel necessitates complex and time-consuming efforts to assess the impact of thermal blooming effect on the optical path.To meet the engineering need for rapid evaluation of thermal blooming effect in optical paths,this study proposed a rapid simulation method for the thermal blooming effect in internal optical paths based on the finite element method.This method discretized the fluid region into infinitesimal elements and employed finite element method for flow field analysis.A simplified analytical model of the flow field region in complex internal channels was established,and regions with similar thermal blooming effect were divided within this model.Based on the calculated optical path differences within these regions,numerical simulations of phase distortion caused by thermal blooming were conducted.The calculated result were compared with those obtained using the existing methods.The findings reveal that for complex optical paths,the discrepancy between the two approaches is less than 3.6%,with similar phase distortion patterns observed.For L-type units,this method and the existing methods identify the same primary factors influencing aberrations and exhibit consistent trends in their variation.This method was used to analyze the impact of thermal blooming effect in a straight channel under different gravity directions.The results show that phase distortion varies with changes in the direction of gravity,and the magnitude of the phase difference is strongly correlated with the component of gravity perpendicular to the optical axis.Compared to the existing methods,this approach offers greater flexibility,obviates the need for complex custom analysis programming.The analytical results of this method enable a rapid assessment of the thermal blooming effect in optical paths within the internal channel.This is especially useful during the engineering design.These results also provide crucial references for developing strategies to suppress thermal blooming effect.
文摘Most of the existing studies on tunnel blast wave are based on spherical or grouped charges, however,conventional weapons are mostly cylindrical rather than spherical. In order to analyze the impact of cylindrical charges on the tunnel blast wave loads and to develop a quantitative calculation method, this study carried out experimental and numerical research. Initially, external explosion experiments were conducted using both 35 kg spherical charges and cylindrical charges with aspect ratio of 4.8 at two different distances from the tunnel entrance. Comparative analysis of the blast wave parameters in the tunnel revealed that the explosive equivalent of the cylindrical charges was significantly higher than that of the spherical charges. To address this, an equivalent coefficient κ based on the spherical charges was proposed for the cylindrical charges. Subsequently, numerical simulations were conducted for the experimental conditions, and the numerical simulation results match the experiments well. Through numerical calculations, the reliability of the equivalent coefficient κ under the experimental conditions was verified, and comparison analysis indicated that the explosion energy of cylindrical charges spreads more radially, resulting in more explosion energy entering the tunnel, which is the fundamental reason for the increase in tunnel blast wave loads. Additionally, analyzing the explosion energy ratio entering the tunnel is an effective method for calculating the equivalent coefficient κ. Finally, through more than one hundred sets of numerical calculation results, the impact of the proportional distance λ and the ratio of charge mass to the tunnel cross-section dimension φ on the equivalence coefficients κ was investigated. An empirical formula for the equivalence coefficient κ was derived through fitting, and the accuracy of the formula was validated through literature experimental results. The research findings of this paper will provide valuable guidance for the calculation of blast wave loads in tunnel.
基金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(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(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 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(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.
基金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(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.
文摘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 (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(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.
基金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.
基金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.
文摘High-speed train running in the sand environment is different from the general environment. In the former situation, there will be sand load applied on high-speed train(SLAHT) caused by sand particles hitting train surface. This will have a great impact on the train stability, running drag and surface corrosion. Numerical simulation method of SLAHT in sand environment is studied. The velocity and mass flow rate models of saltation and suspension sand particles and the calculation model of SLAHT caused by sand particles hitting train surface are established. The discrete phase method is adopted for numerical simulating the process of saltation and suspension sand particles moving to train surface and generating sand load. By comparison with the field tests, the numerical simulation reliability is analysed. The theoretical formula of SLAHT changing with cross-wind and train speed is proposed. SLAHT changing law is analyzed. Research results indicate that SLAHT changing with cross-wind and train speed is a quadratic relationship. When train speed is constant, SLAHT increases quadratically with cross-wind speed improvement. When cross-wind speed is constant, SLAHT increases quadratically with train speed improvement.
