In order to study the effect of rotation on the combustion in the underwater vehicle,a two-phase turbulent combustion process is described with Reynolds stress turbulence model,eddy-dissipation turbulent combustion mo...In order to study the effect of rotation on the combustion in the underwater vehicle,a two-phase turbulent combustion process is described with Reynolds stress turbulence model,eddy-dissipation turbulent combustion model,P-1 radiation model and particle tracking model of liquid. The flow in the rotating combustion chamber is simulated at two different working speeds,0?r/min and 1?000?r/min by Fluent software. The temperature,gas velocity,static pressure of wall and fuel concentration are computed and compared. The results show that the combustion in rotating combustor is faster and more effective.展开更多
In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fi...In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.展开更多
The 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 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.展开更多
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.展开更多
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.展开更多
We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficienc...We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.展开更多
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.展开更多
The present study is devoted to researching the thermal security problems of large-scale solid rocket motor with Ammonium perchlorate/Hydroxyl-terminated polybutadiene(AP/HTPB). A two-dimensional axisymmetric model fo...The present study is devoted to researching the thermal security problems of large-scale solid rocket motor with Ammonium perchlorate/Hydroxyl-terminated polybutadiene(AP/HTPB). A two-dimensional axisymmetric model for the cook-off of solid rocket motor is established. The reaction kinetics for the cook-off process of AP/HTPB is described by the two-step global chemical mechanism. Numerical predictions of the cook-off behavior for the propellant are conducted at fast heating rate of 1.45-2.45 K/s,and slow heating rate of 0.001-0.003 K/s, respectively. The results show that in the fast cook-off condition. the initial ignition position of AP/HTPB occurs in the annular region of the outer wall of propellant without exception, and the center point in the region is(889.1,149.5). For the region, the axial width is1.8 mm and radial thickness is 0.8 mm. However, in the slow cook-off condition, the ignition center position is shifted along the axial direction toward the right end face of the propellant with the increase of heating rate. Therefore, the influence of heating rate on ignition temperature and ignition delay time is nonnegligible within a certain range.展开更多
Based on the engineering background of a soft rock roadway in Qinan Coal Mine 82 Area,Huaibei Mining Group,three creep models with different support patterns in soft rock roadway were established by using geotechnical...Based on the engineering background of a soft rock roadway in Qinan Coal Mine 82 Area,Huaibei Mining Group,three creep models with different support patterns in soft rock roadway were established by using geotechnical software of FLAC2D.According to the calculation results of different models,the change law of mechanical properties with the time of bolt-grouting support structure was obtained.Furthermore,for the test bolt-grouting support roadway,the deformation law of surrounding rock got by underground industrial experiment and field observation accords with the creep law got by numerical calculation.The results of numerical calculation and field observation show that,compared with other supports,the creep of bolt-grouting support roadway enters the steady-state creep stage from tertiary creep stage ahead,the deformations of roof,floor and two sides are decreased greatly,the plastically deforming area in surrounding rock is reduced obviously,and the distribution ranges of maximum and minimum principal stress are shrank obviously.All those fully show that the bolt-grouting support has its remarkable advantages in controlling surrounding rock creep and improving the whole strength of surrounding rock and self-bearing capacity.展开更多
Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and opt...Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and optimized remeshing function, especially for the large deformation. In order to use this function sufficiently, simulation of cold rolling involutes spline can be implemented indirectly. The relationship between die and workpiece, forming force and characteristic of deformation in the forming process of cold rolling involutes spline are analyzed and researched. Meanwhile, reliable proofs for the design of dies and deforming equipment are provided.展开更多
Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model a...Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.展开更多
On the basis of the two dimensional finite element analysis model, the pile foundations' mechanical effect of the rigid pile composite foundation under the dynamic load was researched. Through the research, the de...On the basis of the two dimensional finite element analysis model, the pile foundations' mechanical effect of the rigid pile composite foundation under the dynamic load was researched. Through the research, the development law and deformation property of axial force of pile body, shaft resistance of pile, and cumulative settlement of pile head under vertical cyclic dynamic loads were concluded. Through the comparison and analysis of the test results of dynamic models, the test results of Poulos(1989) and cumulative settlement model of the single pile under cyclic loads were confirmed. Based on the above research, Fortran language was adopted to introduce the soil attenuation factor, the secondary development of relevant modules of ABAQUS was carried out, and the effect of soil attenuation factor on dynamic property of pile-soil was discussed further.展开更多
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.展开更多
High temperature air combustion (HTAC) is an attractive technology of saving energy and controlling environment. The mathematical models of turbulent jet flame under the highly preheated air combustion condition are c...High temperature air combustion (HTAC) is an attractive technology of saving energy and controlling environment. The mathematical models of turbulent jet flame under the highly preheated air combustion condition are conducted in the paper. The mixture fraction/probability density function model is employed. The results show that the maximum flame temperature is decreased, the temperature in the HTAC furnace is more uniform than that in the conventional furnace, and the NO x emission is low. The numerical results are partially validated by some experimental measurements.展开更多
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.展开更多
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.展开更多
The design of solid armature of railgun should take full account of its operating conditions and material properties because the armature is subjected to dynamic loading conditions and experiences a complicated electr...The design of solid armature of railgun should take full account of its operating conditions and material properties because the armature is subjected to dynamic loading conditions and experiences a complicated electrical,thermal and mechanical process in the interior ballistic cycle.In this paper present,we first introduced a multi-physical field model of railgun,followed by several examples to investigate the launching process.Especially,we used the explicit finite element method,in which material nonlinearity and geometric nonlinearity were accounted,to investigate the deform behaviors of solid armature.The results show that the dynamic mechanical process of armature is dependent on the armature geometry,material and exciting electric current.By the numerical simulation,the understanding of the fracture mechanism of solid armature was deepened.展开更多
The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is establishe...The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program.Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25-0.5 cm and 1-3 cm.The Euler--fluxcorrected transport solver is found to be suitable for the far-field calculation after mapping.A numerical simulation is conducted to study the formation,propagation,and interaction of the shock wave of the composite charge for different initiation modes.It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge.Additionally,it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation,and the peak pressure of the shock wave exponentially decays,fitting the calculation formula of the peak overpressure attenuation under different initiation modes,which is obtained and verified by experiment.The difference between numerical and experimental results is less than 10%,and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.展开更多
文摘In order to study the effect of rotation on the combustion in the underwater vehicle,a two-phase turbulent combustion process is described with Reynolds stress turbulence model,eddy-dissipation turbulent combustion model,P-1 radiation model and particle tracking model of liquid. The flow in the rotating combustion chamber is simulated at two different working speeds,0?r/min and 1?000?r/min by Fluent software. The temperature,gas velocity,static pressure of wall and fuel concentration are computed and compared. The results show that the combustion in rotating combustor is faster and more effective.
基金Project(51038004) supported by the National Natural Science Foundation of China
文摘In order to verify the validity of finite element numerical simulation method for asphalt mixture, which consists of aggregates, mastic （where mastic is a kind of fine mixture composed of asphalt binder mixed with fines and fine aggregates） and air voids, based on three-dimensional （3D） heterogeneous specimen, X-ray computerized tomography （X-ray CT） was used to scan the asphalt specimens to obtain the real internal microstrnctures of asphalt mixture. CT images were reconstructed to build up 3D digital specimen, and the viscoelastic properties of mastic were described with Burgers model The uniaxial creep numerical simulations of three different levels of aggregate gradation were conducted. The simulation results agree well with the experimental results. The numerical simulation of asphalt mixture incorporated with real 3D microstructure based on finite element method is a promising application to conduct research of asphalt concrete. Additionally, this method can increase the mechanistic understanding of global viscoelastic properties of asphalt mixtures by linking the real 3D microstructure.
文摘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 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 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.
基金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.
基金Open Fund of State Key Laboratory of Geological Processes and Mineral Resources, China University of Geo-sciences (GPMR0750)National Natural Science Foundation of China (40437018)
文摘We apply the forward modeling algorithm constituted by the convolutional Forsyte polynomial differentiator pro-posed by former worker to seismic wave simulation of complex heterogeneous media and compare the efficiency and accuracy between this method and other seismic simulation methods such as finite difference and pseudospec-tral method. Numerical experiments demonstrate that the algorithm constituted by convolutional Forsyte polyno-mial differentiator has high efficiency and accuracy and needs less computational resources, so it is a numerical modeling method with much potential.
文摘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.
文摘The present study is devoted to researching the thermal security problems of large-scale solid rocket motor with Ammonium perchlorate/Hydroxyl-terminated polybutadiene(AP/HTPB). A two-dimensional axisymmetric model for the cook-off of solid rocket motor is established. The reaction kinetics for the cook-off process of AP/HTPB is described by the two-step global chemical mechanism. Numerical predictions of the cook-off behavior for the propellant are conducted at fast heating rate of 1.45-2.45 K/s,and slow heating rate of 0.001-0.003 K/s, respectively. The results show that in the fast cook-off condition. the initial ignition position of AP/HTPB occurs in the annular region of the outer wall of propellant without exception, and the center point in the region is(889.1,149.5). For the region, the axial width is1.8 mm and radial thickness is 0.8 mm. However, in the slow cook-off condition, the ignition center position is shifted along the axial direction toward the right end face of the propellant with the increase of heating rate. Therefore, the influence of heating rate on ignition temperature and ignition delay time is nonnegligible within a certain range.
基金Project(106084) supported by the Scientific and Technological Research of the Ministry of EducationProject(BK2007701) supported by the Natural Science Foundation of Jiangsu Province,ChinaProject(2006CB2022010) supported by the National Basic Research Program of China and the Qing-lan Project of Jiangsu Province,China
文摘Based on the engineering background of a soft rock roadway in Qinan Coal Mine 82 Area,Huaibei Mining Group,three creep models with different support patterns in soft rock roadway were established by using geotechnical software of FLAC2D.According to the calculation results of different models,the change law of mechanical properties with the time of bolt-grouting support structure was obtained.Furthermore,for the test bolt-grouting support roadway,the deformation law of surrounding rock got by underground industrial experiment and field observation accords with the creep law got by numerical calculation.The results of numerical calculation and field observation show that,compared with other supports,the creep of bolt-grouting support roadway enters the steady-state creep stage from tertiary creep stage ahead,the deformations of roof,floor and two sides are decreased greatly,the plastically deforming area in surrounding rock is reduced obviously,and the distribution ranges of maximum and minimum principal stress are shrank obviously.All those fully show that the bolt-grouting support has its remarkable advantages in controlling surrounding rock creep and improving the whole strength of surrounding rock and self-bearing capacity.
文摘Design of forming dies and whole process of simulation of cold rolling involutes spline can be realized by using of CAD software of PRO-E and CAE software of DEFORM-3D. Software DEFORM-3D provides an automatic and optimized remeshing function, especially for the large deformation. In order to use this function sufficiently, simulation of cold rolling involutes spline can be implemented indirectly. The relationship between die and workpiece, forming force and characteristic of deformation in the forming process of cold rolling involutes spline are analyzed and researched. Meanwhile, reliable proofs for the design of dies and deforming equipment are provided.
基金Project supported by the Scientific Research Foundation of Hebei University of Technology of China
文摘Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.
基金Projects(51478178,51508181) supported by the National Natural Science Foundation of China
文摘On the basis of the two dimensional finite element analysis model, the pile foundations' mechanical effect of the rigid pile composite foundation under the dynamic load was researched. Through the research, the development law and deformation property of axial force of pile body, shaft resistance of pile, and cumulative settlement of pile head under vertical cyclic dynamic loads were concluded. Through the comparison and analysis of the test results of dynamic models, the test results of Poulos(1989) and cumulative settlement model of the single pile under cyclic loads were confirmed. Based on the above research, Fortran language was adopted to introduce the soil attenuation factor, the secondary development of relevant modules of ABAQUS was carried out, and the effect of soil attenuation factor on dynamic property of pile-soil was discussed further.
基金Project(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.
文摘High temperature air combustion (HTAC) is an attractive technology of saving energy and controlling environment. The mathematical models of turbulent jet flame under the highly preheated air combustion condition are conducted in the paper. The mixture fraction/probability density function model is employed. The results show that the maximum flame temperature is decreased, the temperature in the HTAC furnace is more uniform than that in the conventional furnace, and the NO x emission is low. The numerical results are partially validated by some experimental measurements.
文摘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.
基金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.
文摘The design of solid armature of railgun should take full account of its operating conditions and material properties because the armature is subjected to dynamic loading conditions and experiences a complicated electrical,thermal and mechanical process in the interior ballistic cycle.In this paper present,we first introduced a multi-physical field model of railgun,followed by several examples to investigate the launching process.Especially,we used the explicit finite element method,in which material nonlinearity and geometric nonlinearity were accounted,to investigate the deform behaviors of solid armature.The results show that the dynamic mechanical process of armature is dependent on the armature geometry,material and exciting electric current.By the numerical simulation,the understanding of the fracture mechanism of solid armature was deepened.
基金funded by the National Natural Science Foundation of China under NO.11202103Qing-lan Project of Jiangsu Province。
文摘The numerical simulation of a blast wave of a multilayer composite charge is investigated.A calculation model of the near-field explosion and far-field propagation of the shock wave of a composite charge is established using the AUTODYN finite element program.Results of the near-field and far-field calculations of the shock wave respectively converge at cell sizes of 0.25-0.5 cm and 1-3 cm.The Euler--fluxcorrected transport solver is found to be suitable for the far-field calculation after mapping.A numerical simulation is conducted to study the formation,propagation,and interaction of the shock wave of the composite charge for different initiation modes.It is found that the initiation mode obviously affects the shock-wave waveform and pressure distribution of the composite charge.Additionally,it is found that the area of the overpressure distribution is greatest for internal and external simultaneous initiation,and the peak pressure of the shock wave exponentially decays,fitting the calculation formula of the peak overpressure attenuation under different initiation modes,which is obtained and verified by experiment.The difference between numerical and experimental results is less than 10%,and the peak overpressure of both internal and external initiation is 56.12% higher than that of central single-point initiation.
