With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ...With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.展开更多
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.展开更多
To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based sim...To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.展开更多
In this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbul...In this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbulence structure is described using the shear-stress transport k-ω(SST k-ω) model, and the volume of fluid(VOF) method is used to track the complex air-liquid interface. The motion of spheres during water entry is simulated using an independent overset grid. The numerical model is verified by comparing the cavity evolution results from simulations and experiments. Numerical results reveal that the time interval between the twin water entries evidently affects cavity expansion and contraction behaviors in the radial direction. However, this influence is significantly weakened by increasing the lateral distance between the two spheres. In synchronous water entries, pressure is reduced on the midline of two cavities during surface closure, which is directly related to the cavity volume. The evolution of vortexes inside the two cavities is analyzed using a velocity vector field, which is affected by the lateral distance and time interval of water entries.展开更多
The complex structure of the bottom of a high-speed train is an important source of train aerodynamic drag.Thus,improving the bottom structure is of great significance to reduce the aerodynamic drag of the train.In th...The complex structure of the bottom of a high-speed train is an important source of train aerodynamic drag.Thus,improving the bottom structure is of great significance to reduce the aerodynamic drag of the train.In this study,computational fluid dynamics(CFD)based on three-dimensional steady incompressible Reynolds-average Naiver-Stokes(RANS)equations and Realizable k-εturbulence model were utilized for numerical simulations.Inspired by the concept of streamlined design and the idea of bottom flow field control,this study iteratively designed the bogies in a streamlined shape and combined them with the bottom deflectors to investigate the joint drag reduction mechanism.Three models,i.e.,single-bogie model,simplified train model,and eight-car high-speed train model,were created and their aerodynamic characteristics were analyzed.The results show that the single-bogie model with streamlined design shows a noticeable drag reduction,whose power bogie and trailer bogie experience 13.92%and 7.63%drag reduction,respectively.The range of positive pressure area on the bogie is reduced.The aerodynamic drag can be further reduced to 15.01%by installing both the streamlined bogie and the deflector on the simplified train model.When the streamlined bogies and deflectors are used on the eight-car model together,the total drag reduction rate reaches 2.90%.Therefore,the proposed aerodynamic kit for the high-speed train bottom is capable to improve the flow structure around the bogie regions,reduce the bottom flow velocity,and narrow the scope of the train’s influence on the surrounding environment,achieving the appreciable reduction of aerodynamic drag.This paper can provide a new idea for the drag reduction of high-speed trains.展开更多
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.展开更多
This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare ...This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.展开更多
Aerocapture is one of the key technologies for low-cost transportation,with high demands of autonomy,accuracy,and robustness of guidance and control,due to its high reliability requirements for only one chance of tryi...Aerocapture is one of the key technologies for low-cost transportation,with high demands of autonomy,accuracy,and robustness of guidance and control,due to its high reliability requirements for only one chance of trying.A unified numerical predictor-corrector guidance method based on characteristic models for aerocapture is proposed.The numerical predictor-corrector guidance method is used to achieve autonomy and high accuracy,and the characteristic model control method is introduced to achieve robustness.At the same time,by transforming path constraints,characteristic model equations including apogee deviation and altitude differentiation are established.Based on the characteristic model equations,a unified guidance law which can satisfy path constraints and guidance objectives simultaneously is designed.In guidance problems,guidance deviation is not directly obtained from the output of the dynamics at present,but is calculated through integral and algebraic equations.Therefore,the method of directly discretizing differential equations cannot be used to establish characteristic models,which brings great difficulty to characteristic modeling.A method for characteristic modeling of guidance problems is proposed,and convergence analysis of the proposed guidance law is also provided.Finally,a joint numerical simulation of guidance and control considering navigation deviation and various uncertainties is conducted to verify the effectiveness of the proposed method.The proposed unified method can be extended to general aerodynamic entry guidance designs,providing theoretical and methodological support for them.展开更多
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 utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,...The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.展开更多
In this paper,an intelligent control method applying on numerical virtual flight is proposed.The proposed algorithm is verified and evaluated by combining with the case of the basic finner projectile model and shows a...In this paper,an intelligent control method applying on numerical virtual flight is proposed.The proposed algorithm is verified and evaluated by combining with the case of the basic finner projectile model and shows a good application prospect.Firstly,a numerical virtual flight simulation model based on overlapping dynamic mesh technology is constructed.In order to verify the accuracy of the dynamic grid technology and the calculation of unsteady flow,a numerical simulation of the basic finner projectile without control is carried out.The simulation results are in good agreement with the experiment data which shows that the algorithm used in this paper can also be used in the design and evaluation of the intelligent controller in the numerical virtual flight simulation.Secondly,combined with the real-time control requirements of aerodynamic,attitude and displacement parameters of the projectile during the flight process,the numerical simulations of the basic finner projectile’s pitch channel are carried out under the traditional PID(Proportional-Integral-Derivative)control strategy and the intelligent PID control strategy respectively.The intelligent PID controller based on BP(Back Propagation)neural network can realize online learning and self-optimization of control parameters according to the acquired real-time flight parameters.Compared with the traditional PID controller,the concerned control variable overshoot,rise time,transition time and steady state error and other performance indicators have been greatly improved,and the higher the learning efficiency or the inertia coefficient,the faster the system,the larger the overshoot,and the smaller the stability error.The intelligent control method applying on numerical virtual flight is capable of solving the complicated unsteady motion and flow with the intelligent PID control strategy and has a strong promotion to engineering application.展开更多
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.展开更多
Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensi...Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.展开更多
The NATO agreement STANAG 4569 defines the protection levels for the occupants of logistic and light armored vehicle.The Allied Engineering Publication,AEP-55,Volume 2 document outlines the test conditions for underbe...The NATO agreement STANAG 4569 defines the protection levels for the occupants of logistic and light armored vehicle.The Allied Engineering Publication,AEP-55,Volume 2 document outlines the test conditions for underbelly improvised explosive device(IEDs),which must be buried in water-saturated sandy gravel.The use of sandy gravel has some drawbacks,for instance reproducibility,time consumption,and cost.This paper focuses on the investigation of four alternatives to sandy gravel,which could produce similar specific and cumulative impulses:a concrete pot filled with water,a concrete pot filled with quartz sand,a steel pot without filling and a concrete pot filled with glass spheres(diameter 200μm—300μm)and different water contents.The impulses are measured with a ring technology developed at the Fraunhofer EMI.A numerical soil model based on the work of Marrs,2014 and Fi serov a,2006 and considering the soil moisture was used to simulate the experiments with glass spheres at different water contents,showing much better agreement with the experiments than the classical Laine&Sandvik model,even for high saturation levels.These results can be used to create new test conditions at original scale that are more cost-effective,more reproducible and simpler to manage in comparison to the current tests carried out with STANAG sandy gravel.展开更多
Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlin...Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlinear solid media by the numerical method. Distinctly different from the case of numerical modeling of linear elastic wave, there may be several difficulties in the numerical treatment to the nonlinear partial differential equation, such as the steep gradients, shocks and unphysical oscillations. All of them are the great obstacles to the stability and conver-gence of numerical calculation. Fortunately, the comparative study on the modeling of nonlinear wave by the two FD schemes presented in the paper can provide us with an easy method to keep the stability and convergence in the calculation field when the product of the absolute value of nonlinear coefficient and the value of u/x are small enough, namely, the value of bu/x is much smaller than 1. Several results are founded in the numerical study of nonlinear P-wave propagation, such as the waveform aberration, the generation and growth of harmonic wave and the energy redistribution among different frequency components. All of them will be more violent when the initial amplitude A0 is larger or the nonlinearity of medium is stronger. Correspondingly, we have found that the nonlinear P-wave propagation velocity will change with different initial frequency f of source wave or the wave velocity c (equal to the P-wave velocity in the same medium without considering nonlinearity).展开更多
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.展开更多
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.展开更多
Pile foundation is widely used in the offshore engineering. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. The potentials of liquefaction and damages of pile foundation due to ...Pile foundation is widely used in the offshore engineering. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. The potentials of liquefaction and damages of pile foundation due to the liquefaction can be reduced by the implementation of the drainage in the liquefiable foundation. A patented pile technology, named rigidity-drain pile, was introduced. The partial section of the pile body was filled by materials with higher penetrability which forms some effective drainage channels in the pile. The principles and construction methods were presented. 3D models for both rigidity-drain pile and ordinary pile were built in FLAC3D code. The dynamic loadings were applied on the bottom of the model. According to the numerical results, in the case of the rigidity-drain pile, the water in the relevant distance range around the pile flows toward the pile drainage, the contour of the pore pressure shows a funnel form. Contrast to the ordinary pile, the rigidity-drain pile can dissipate the accumulated excess pore water, maintain effective stress and obviously reduce the possibility of surrounding soil liquefaction.展开更多
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.展开更多
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.展开更多
文摘With different structural forms of ventilation pipes have various attenuation effects on incident shock waves while meeting ventilation requirements.The attenuation mechanism and the propagation law of shock waves in ventilation pipes of different structures are investigated by experiments and numerical simulations.Furthermore,for the same structure,the effects of peak pressure and positive pressure time on the attenuation rate are discussed.It is found that the attenuation rate increases with the incident shock wave pressure,and the shock wave attenuation rate tends to reach its limiting value k for the same structure and reasonably short positive pressure time.Under the same conditions,the attenuation rate is calculated using the pressure of the shock wave as follows:diffusion chamber pipe,branch pipe and selfconsumption pipe;the attenuation rate per unit volume is calculated as follows:self-consumption pipe,branch pipe and diffusion chamber pipe.In addition,an easy method is provided to calculate the attenuation rate of the shock wave in single and multi-stage ventilation pipes.Corresponding parameters are provided for various structures,and the margin of error between the formulae and experimental results is within 10%,which is significant for engineering applications.
文摘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.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271317 and 52071149)the Fundamental Research Funds for the Central Universities(HUST:2019kfy XJJS007)。
文摘To address the problems of low accuracy by the CONWEP model and poor efficiency by the Coupled Eulerian-Lagrangian(CEL)method in predicting close-range air blast loads of cylindrical charges,a neural network-based simulation(NNS)method with higher accuracy and better efficiency was proposed.The NNS method consisted of three main steps.First,the parameters of blast loads,including the peak pressures and impulses of cylindrical charges with different aspect ratios(L/D)at different stand-off distances and incident angles were obtained by two-dimensional numerical simulations.Subsequently,incident shape factors of cylindrical charges with arbitrary aspect ratios were predicted by a neural network.Finally,reflected shape factors were derived and implemented into the subroutine of the ABAQUS code to modify the CONWEP model,including modifications of impulse and overpressure.The reliability of the proposed NNS method was verified by related experimental results.Remarkable accuracy improvement was acquired by the proposed NNS method compared with the unmodified CONWEP model.Moreover,huge efficiency superiority was obtained by the proposed NNS method compared with the CEL method.The proposed NNS method showed good accuracy when the scaled distance was greater than 0.2 m/kg^(1/3).It should be noted that there is no need to generate a new dataset again since the blast loads satisfy the similarity law,and the proposed NNS method can be directly used to simulate the blast loads generated by different cylindrical charges.The proposed NNS method with high efficiency and accuracy can be used as an effective method to analyze the dynamic response of structures under blast loads,and it has significant application prospects in designing protective structures.
基金China Academy of Launch Vehicle Technology(Grant No.CALT-2022-03)Science and Technology on Underwater Information and Control Laboratory(Grant No.2021-JCJQ-LB-030-05).
文摘In this study, a three dimensional(3D) numerical model of six-degrees-of-freedom(6DOF) is applied to simulate the water entries of twin spheres side-by-side at different lateral distances and time intervals.The turbulence structure is described using the shear-stress transport k-ω(SST k-ω) model, and the volume of fluid(VOF) method is used to track the complex air-liquid interface. The motion of spheres during water entry is simulated using an independent overset grid. The numerical model is verified by comparing the cavity evolution results from simulations and experiments. Numerical results reveal that the time interval between the twin water entries evidently affects cavity expansion and contraction behaviors in the radial direction. However, this influence is significantly weakened by increasing the lateral distance between the two spheres. In synchronous water entries, pressure is reduced on the midline of two cavities during surface closure, which is directly related to the cavity volume. The evolution of vortexes inside the two cavities is analyzed using a velocity vector field, which is affected by the lateral distance and time interval of water entries.
基金Project(2020YFA0710901)supported by the National Key Research and Development Program of ChinaProject(2023JJ30643)supported by the Natural Science Foundation of Hunan Province,China+1 种基金Project(12372204)supported by the National Natural Science Foundation of ChinaProject(2022ZZTS0725)supported by the Self-exploration and Innovation Project for Postgraduates of Central South University,China。
文摘The complex structure of the bottom of a high-speed train is an important source of train aerodynamic drag.Thus,improving the bottom structure is of great significance to reduce the aerodynamic drag of the train.In this study,computational fluid dynamics(CFD)based on three-dimensional steady incompressible Reynolds-average Naiver-Stokes(RANS)equations and Realizable k-εturbulence model were utilized for numerical simulations.Inspired by the concept of streamlined design and the idea of bottom flow field control,this study iteratively designed the bogies in a streamlined shape and combined them with the bottom deflectors to investigate the joint drag reduction mechanism.Three models,i.e.,single-bogie model,simplified train model,and eight-car high-speed train model,were created and their aerodynamic characteristics were analyzed.The results show that the single-bogie model with streamlined design shows a noticeable drag reduction,whose power bogie and trailer bogie experience 13.92%and 7.63%drag reduction,respectively.The range of positive pressure area on the bogie is reduced.The aerodynamic drag can be further reduced to 15.01%by installing both the streamlined bogie and the deflector on the simplified train model.When the streamlined bogies and deflectors are used on the eight-car model together,the total drag reduction rate reaches 2.90%.Therefore,the proposed aerodynamic kit for the high-speed train bottom is capable to improve the flow structure around the bogie regions,reduce the bottom flow velocity,and narrow the scope of the train’s influence on the surrounding environment,achieving the appreciable reduction of aerodynamic drag.This paper can provide a new idea for the drag reduction of high-speed trains.
文摘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.
基金This research has been conducted under SEGTRANS project,funded by the Centre for Industrial Technological Development(CDTI,Government of Spain).
文摘This paper explores the performances of a finite element simulation including four concrete models applied to a full-scale reinforced concrete beam subjected to blast loading. Field test data has been used to compare model results for each case. The numerical modelling has been, carried out using the suitable code LS-DYNA. This code integrates blast load routine(CONWEP) for the explosive description and four different material models for the concrete including: Karagozian & Case Concrete, Winfrith, Continuous Surface Cap Model and Riedel-Hiermaier-Thoma models, with concrete meshing based on 10, 15, and 20 mm. Six full-scale beams were tested: four of them used for the initial calibration of the numerical model and two more tests at lower scaled distances. For calibration, field data obtained employing pressure and accelerometers transducers were compared with the results derived from the numerical simulation. Damage surfaces and the shape of rupture in the beams have been used as references for comparison. Influence of the meshing on accelerations has been put in evidence and for some models the shape and size of the damage in the beams produced maximum differences around 15%. In all cases, the variations between material and mesh models are shown and discussed.
基金The National Key R&D Program of China(2018YFA0703800)。
文摘Aerocapture is one of the key technologies for low-cost transportation,with high demands of autonomy,accuracy,and robustness of guidance and control,due to its high reliability requirements for only one chance of trying.A unified numerical predictor-corrector guidance method based on characteristic models for aerocapture is proposed.The numerical predictor-corrector guidance method is used to achieve autonomy and high accuracy,and the characteristic model control method is introduced to achieve robustness.At the same time,by transforming path constraints,characteristic model equations including apogee deviation and altitude differentiation are established.Based on the characteristic model equations,a unified guidance law which can satisfy path constraints and guidance objectives simultaneously is designed.In guidance problems,guidance deviation is not directly obtained from the output of the dynamics at present,but is calculated through integral and algebraic equations.Therefore,the method of directly discretizing differential equations cannot be used to establish characteristic models,which brings great difficulty to characteristic modeling.A method for characteristic modeling of guidance problems is proposed,and convergence analysis of the proposed guidance law is also provided.Finally,a joint numerical simulation of guidance and control considering navigation deviation and various uncertainties is conducted to verify the effectiveness of the proposed method.The proposed unified method can be extended to general aerodynamic entry guidance designs,providing theoretical and methodological support for them.
基金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.
基金Project(NB-2020-JG-07)supported by the Research and Engineering Application of Key Technologies for New Building Industrialization Project of China Northwest Architectural Design and Research Institute Co.,Ltd.Project(2023-CXTD-29)supported by the Key Scientific and Technological Innovation Team of Shaanxi Province,ChinaProject supported by the K.C.Wong Education Foundation。
文摘The utilization of prefabricated light modular radiant heating system has demonstrated significant increases in heat transfer efficiency and energy conservation capabilities.Within prefabricated building construction,this new heating method presents an opportunity for the development of comprehensive facilities.The parameters for evaluating the effectiveness of such a system are the upper surface layer’s heat flux and temperature.In this paper,thermal resistance analysis calculation based on a simplified model for this unique radiant heating system analysis is presented with the heat transfer mechanism’s evaluation.The results obtained from thermal resistance analysis calculation and numerical simulation indicate that the thermal resistance analysis method is highly accurate with temperature discrepancies ranging from 0.44℃ to−0.44℃ and a heat flux discrepancy of less than 7.54%,which can meet the requirements of practical engineering applications,suggesting a foundation for the prefabricated radiant heating system.
文摘In this paper,an intelligent control method applying on numerical virtual flight is proposed.The proposed algorithm is verified and evaluated by combining with the case of the basic finner projectile model and shows a good application prospect.Firstly,a numerical virtual flight simulation model based on overlapping dynamic mesh technology is constructed.In order to verify the accuracy of the dynamic grid technology and the calculation of unsteady flow,a numerical simulation of the basic finner projectile without control is carried out.The simulation results are in good agreement with the experiment data which shows that the algorithm used in this paper can also be used in the design and evaluation of the intelligent controller in the numerical virtual flight simulation.Secondly,combined with the real-time control requirements of aerodynamic,attitude and displacement parameters of the projectile during the flight process,the numerical simulations of the basic finner projectile’s pitch channel are carried out under the traditional PID(Proportional-Integral-Derivative)control strategy and the intelligent PID control strategy respectively.The intelligent PID controller based on BP(Back Propagation)neural network can realize online learning and self-optimization of control parameters according to the acquired real-time flight parameters.Compared with the traditional PID controller,the concerned control variable overshoot,rise time,transition time and steady state error and other performance indicators have been greatly improved,and the higher the learning efficiency or the inertia coefficient,the faster the system,the larger the overshoot,and the smaller the stability error.The intelligent control method applying on numerical virtual flight is capable of solving the complicated unsteady motion and flow with the intelligent PID control strategy and has a strong promotion to engineering application.
基金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.
文摘Understanding the steady mechanism of biomass smoldering plays a great role in the utilization of smoldering technology.In this study numerical analysis of steady smoldering of biomass rods was performed.A two-dimensional(2D)steady model taking into account both char oxidation and pyrolysis was developed on the basis of a calculated propagation velocity according to empirical correlation.The model was validated against the smoldering experiment of biomass rods under natural conditions,and the maximum error was smaller than 31%.Parameter sensitivity analysis found that propagation velocity decreases significantly while oxidation area and pyrolysis zone increase significantly with the increasing diameter of rod fuel.
基金Bundeswehr Technical Center for Weapons and Ammunition WTD-91 GF-440 in Meppen,Germany for funding this work。
文摘The NATO agreement STANAG 4569 defines the protection levels for the occupants of logistic and light armored vehicle.The Allied Engineering Publication,AEP-55,Volume 2 document outlines the test conditions for underbelly improvised explosive device(IEDs),which must be buried in water-saturated sandy gravel.The use of sandy gravel has some drawbacks,for instance reproducibility,time consumption,and cost.This paper focuses on the investigation of four alternatives to sandy gravel,which could produce similar specific and cumulative impulses:a concrete pot filled with water,a concrete pot filled with quartz sand,a steel pot without filling and a concrete pot filled with glass spheres(diameter 200μm—300μm)and different water contents.The impulses are measured with a ring technology developed at the Fraunhofer EMI.A numerical soil model based on the work of Marrs,2014 and Fi serov a,2006 and considering the soil moisture was used to simulate the experiments with glass spheres at different water contents,showing much better agreement with the experiments than the classical Laine&Sandvik model,even for high saturation levels.These results can be used to create new test conditions at original scale that are more cost-effective,more reproducible and simpler to manage in comparison to the current tests carried out with STANAG sandy gravel.
基金Project of Knowledge Innovation Program from Chinese Academy of Sciences (KZCX2-109).
文摘Two central schemes of finite difference (FD) up to different accuracy orders of space sampling step Dx (Fourth order and Sixth order respectively) were used to study the 1-D nonlinear P-wave propagation in the nonlinear solid media by the numerical method. Distinctly different from the case of numerical modeling of linear elastic wave, there may be several difficulties in the numerical treatment to the nonlinear partial differential equation, such as the steep gradients, shocks and unphysical oscillations. All of them are the great obstacles to the stability and conver-gence of numerical calculation. Fortunately, the comparative study on the modeling of nonlinear wave by the two FD schemes presented in the paper can provide us with an easy method to keep the stability and convergence in the calculation field when the product of the absolute value of nonlinear coefficient and the value of u/x are small enough, namely, the value of bu/x is much smaller than 1. Several results are founded in the numerical study of nonlinear P-wave propagation, such as the waveform aberration, the generation and growth of harmonic wave and the energy redistribution among different frequency components. All of them will be more violent when the initial amplitude A0 is larger or the nonlinearity of medium is stronger. Correspondingly, we have found that the nonlinear P-wave propagation velocity will change with different initial frequency f of source wave or the wave velocity c (equal to the P-wave velocity in the same medium without considering nonlinearity).
基金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.
文摘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 (50639010) supported by the National Natural Science Foundation of China
文摘Pile foundation is widely used in the offshore engineering. The pile can be seriously destroyed by the soil liquefaction during strong earthquakes. The potentials of liquefaction and damages of pile foundation due to the liquefaction can be reduced by the implementation of the drainage in the liquefiable foundation. A patented pile technology, named rigidity-drain pile, was introduced. The partial section of the pile body was filled by materials with higher penetrability which forms some effective drainage channels in the pile. The principles and construction methods were presented. 3D models for both rigidity-drain pile and ordinary pile were built in FLAC3D code. The dynamic loadings were applied on the bottom of the model. According to the numerical results, in the case of the rigidity-drain pile, the water in the relevant distance range around the pile flows toward the pile drainage, the contour of the pore pressure shows a funnel form. Contrast to the ordinary pile, the rigidity-drain pile can dissipate the accumulated excess pore water, maintain effective stress and obviously reduce the possibility of surrounding soil liquefaction.
基金Project supported by the Scientific Research Foundation of Hebei University of Technology of China
文摘Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.
基金Project(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.
