This paper presents a simple solution of the dynamic buckling of stiffened plates under in-plane impact loading. Based on large deflection theory, a discretely stiffened plate model has been used. The tangential stres...This paper presents a simple solution of the dynamic buckling of stiffened plates under in-plane impact loading. Based on large deflection theory, a discretely stiffened plate model has been used. The tangential stresses of stiffeners and in-plane displacement are neglected. Appling the Hamilton's principle, the motion equations of stiffened plates are obtained. The deflection of the plate is taken as Fourier series, and using Galerkin method the discrete equations can be deduced, which can be solved easily by Runge-Kutta method. The dynamic buckling loads of the stiffened plates are obtained form Budiansky-Roth criterion.展开更多
In order to correct the test error caused by the dynamic characteristics of pressure sensor and avoid the influence of the error of sensor's dynamic model on compensation results,a dynamic compensation method of the ...In order to correct the test error caused by the dynamic characteristics of pressure sensor and avoid the influence of the error of sensor's dynamic model on compensation results,a dynamic compensation method of the pressure sensor is presented,which is based on quantum-behaved particle swarm optimization(QPSO)algorithm and the mean square error(MSE).By using this method,the inverse model of the sensor is built and optimized and then the coefficients of the optimal compensator are got.This method is verified by the dynamic calibration with shock tube and the dynamic characteristics of the sensor before and after compensation are analyzed in time domain and frequency domain.The results show that the working bandwidth of the sensor is extended effectively.This method can reduce dynamic measuring error and improve test accuracy in actual measurement experiments.展开更多
In this paper, we study the vanishing viscosity limit for non-isentropic gas dy- namics with interacting shocks. Given any entropy solution of non-isentropic gas dynamics which consists of two different families of sh...In this paper, we study the vanishing viscosity limit for non-isentropic gas dy- namics with interacting shocks. Given any entropy solution of non-isentropic gas dynamics which consists of two different families of shocks interacting at some positive time, we show that such solution is the vanishing viscosity limit of a family of smooth global solutions for a viscous system of conservation law. We remark that, after the interacting time, not only shocks but also contact discontinuity are generated.展开更多
BACKGROUND:Hypothermia is associated with poor outcome in trauma patients;however,hemorrhagic shock(HS)model with anesthetized swine was different from that of clinical reality.To identify the effects of environmental...BACKGROUND:Hypothermia is associated with poor outcome in trauma patients;however,hemorrhagic shock(HS)model with anesthetized swine was different from that of clinical reality.To identify the effects of environmental hypothermia on HS,we investigated hemodynamics and oxygen dynamics in an unanesthetized swine model of HS under simulating hypothermia environment.METHODS:Totally 16 Bama pigs were randomly divided into ambient temperature group(group A)and low temperature group(group B),8 pigs in each group.Venous blood(30 mL/kg)was continuously withdrawn for more than 15 minutes in conscious swine to establish a hemorrhagic shock model.Pulmonary arterial temperature(Tp),heart rate(HR),mean arterial pressure(MAP),pulmonary arterial pressure(PAP),pulmonary arterial wedge pressure(PAWP),central venous pressure(CVP),cardiac output(CO),hemoglobin(Hb),saturation of mixed venous blood(SvO_2)and blood gas analysis were recorded at the baseline and different hemorrhagic shock time(HST).The whole body oxygen delivery indices,DO_2l and VO_2l,and the O_2 extraction ratio(O_2ER)were calculated.RESULTS:Core body temperature in group A decreased slightly after the hemorrhagic shock model was established,and environmental hypothermia decreased in core body temperature.The mortality rate was significantly higher in group B(50%)than in group A(0%).DO_2l and VO_2l decreased significantly after hemorrhage.No difference was found in hemodynamics,DO_2l and VO_2l between group A and group B,but the difference in pH,lactic acid and O_2ER was significant between the two groups.CONCLUSION:Environmental hypothermia aggravated the disorder of oxygen metabolism after hemorrhagic shock,which was associated with poor prognosis.展开更多
The power field of shock waves and fragments is analyzed and set up, and the damage modes of stiffened plates are put forward. According to the structural characters of the stiffened plates investigated and the proper...The power field of shock waves and fragments is analyzed and set up, and the damage modes of stiffened plates are put forward. According to the structural characters of the stiffened plates investigated and the properties of the shock waves and fragments, the experiments on the shock waves acting on the stiffened plates (penetrated and non-penetrated by fragments) are mainly conducted. The dynamic response rules of stiffened plates with holes under shock waves and fragments loading are obtained. The results show that the penetration of fragments into stiffened plates hardly affects their deformation produced by shock waves..展开更多
An analytical method is developed to determine the transient response of dynamic thermostress in a two-layered cylinder with initial interface pressure. At first, the initial interface pressure in a two-layered cylind...An analytical method is developed to determine the transient response of dynamic thermostress in a two-layered cylinder with initial interface pressure. At first, the initial interface pressure in a two-layered cylinder caused by a heat-assembling method is considered as the initial condition of a thermal elastodynamic equilibrium equation. Thus, a thermal elastodynamic solution for a separate hollow cylinder with the initial stress field is found out by means of a series of simply mathematical transform. By making use of the boundary conditions and continuity conditions of a layered cylinders, the solution for the thermal shock exerting an influence on the initial interface pressure in a two-layered cylinder is also discussed.展开更多
A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simu...A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simulation models for high-speed piston impacts on explosive supercells were established.Simulations were also performed to investigate shock-induced reactions of various high-energy explosives.The fraction of reacted explosive molecules in an initial supercell changed linearly with the propagation distance of the shock-wave front.The corresponding slope could be used as a reaction rate for a specific shock-loading velocity.Reaction rates that varied with the shock-loading pressure exhibited two-stage linearities with different slopes.The two inflection points corresponded to the initial and accelerated reactions,which respectively correlated to the thresholds of shock-induced ignition and detonation.Therefore,the ignition and detonation critical pressures could be determined.The sensitivity could then be a quantitative prediction of the critical pressure.The accuracies of the quantitative shock sensitivity predictions were verified by comparing the impact and shock sensitivities of common explosives and the characteristics of anisotropic shock-induced reactions.Molecular dynamics simulations quantitatively predict and rank shock sensitivities by using only crystal structures of the explosives.Overall,this method will enable the design and safe use of explosives.展开更多
A finite-difference algorithm is proposed for numerical modeling of hydrodynamic flows with rarefaction shocks, in which the fluid undergoes a jump-like liquid-gas phase transition. This new type of flow discontinuity...A finite-difference algorithm is proposed for numerical modeling of hydrodynamic flows with rarefaction shocks, in which the fluid undergoes a jump-like liquid-gas phase transition. This new type of flow discontinuity, unexplored so far in computational fluid dynamics, arises in the approximation of phase-flip(PF) hydrodynamics, where a highly dynamic fluid is allowed to reach the innermost limit of metastability at the spinodal, upon which an instantaneous relaxation to the full phase equilibrium(EQ) is assumed. A new element in the proposed method is artificial kinetics of the phase transition, represented by an artificial relaxation term in the energy equation for a "hidden"component of the internal energy, temporarily withdrawn from the fluid at the moment of the PF transition. When combined with an appropriate variant of artificial viscosity in the Lagrangian framework, the latter ensures convergence to exact discontinuous solutions, which is demonstrated with several test cases.展开更多
Using molecular dynamics (MD) simulation, we study the thermal shock behavior of tungsten (W), which has been used for the plasma facing material (PFM) of tokamaks. The thermo-elastic stress wave, corresponding ...Using molecular dynamics (MD) simulation, we study the thermal shock behavior of tungsten (W), which has been used for the plasma facing material (PFM) of tokamaks. The thermo-elastic stress wave, corresponding to the collective displacement of atoms, is analyzed with the Lagrangian atomic stress method, of which the reliability is also analyzed. The stress wave velocity corresponds to the speed of sound in the material, which is not dependent on the thermal shock energy. The peak pressure of a normal stress wave increases with the increase of thermal shock energy. We analyze the temperature evolution of the thermal shock region according to the Fourier transformation. It can be seen that the “obvious” velocity of heat propagation is less than the velocity of the stress wave; further, that the thermo-elastic stress wave may contribute little to the transport of kinetic energy. The heat propagation can be described properly by the heat conduction equation. These results may be useful for understanding the process of the thermal shock of tungsten.展开更多
The mechanical response of a single crystal titanium sample against(0001)α surface impact was investigated using molecular dynamics simulation.Remarkably,non-uniform plastic deformation was observed in the sample.At ...The mechanical response of a single crystal titanium sample against(0001)α surface impact was investigated using molecular dynamics simulation.Remarkably,non-uniform plastic deformation was observed in the sample.At high strain rates,amorphization occurred near the edge of the contact region where severe shear strain induced a large number of stacking faults(SFs)and dislocations.In contrast,the central part of the contact region underwent less deformation with significantly fewer dislocations.Moreover,instead of amorphization by consuming SFs and dislocations,there was a gradual increase in the density of dislocations and SFs during the process of amorphization.These local amorphous regions eventually grew into shear bands.展开更多
Chemical inclusions significantly alter shock responses of crystalline explosives in macroscale gap experiments but their microscale dynamics origin remains unclear.Herein shock-induced energy localization,overall phy...Chemical inclusions significantly alter shock responses of crystalline explosives in macroscale gap experiments but their microscale dynamics origin remains unclear.Herein shock-induced energy localization,overall physical responses,and reactions in a-1,3,5-trinitro-1,3,5-triazinane(a-RDX)crystal entrained various chemical inclusions were investigated by the multi-scale shock technique implemented in the reactive molecular dynamics method.Results indicated that energy localization and shock reaction were affected by the intrinsic factors within chemical inclusions,i.e.,phase states,chemical compositions,and concentrations.The atomic origin of chemical-inclusions effects on energy localization is dependent on the dynamics mechanism of interfacial molecules with free space volume,which includes homogeneous intermolecular compression,interfacial impact and shear,and void collapse and jet.As introducing various chemical inclusions,the initiation of those dynamics mechanisms triggers diverse decay rates of bulk RDX molecules and hereby impacts on growth speeds of final reactions.Adding chemical inclusions can reduce the effectiveness of the void during the shock impacting.Under the shockwave velocity of 9 km/s,the parent RDX decay rate in RDX entrained amorphous carbon decreases the most and is about one fourth of that in RDX with a vacuum void,and solid HMX and TATB inclusions are more reactive than amorphous carbon but less reactive than dry air or acetone inclusions.The lessdense shocking system denotes the greater increases in local temperature and stress,the faster energy liberation,and the earlier final reaction into equilibrium,revealing more pronounced responses to the present intense shockwave.The quantitative models associated with the relative system density(RD_(sys))were proposed for indicating energy-localization mechanisms and evaluating initiation safety in the shocked crystalline explosive.RD_(sys)is defined by the density ratio of defective RDX to perfect crystal after dynamics relaxation and reveals the global density characteristic in shocked systems filled with chemical inclusions.When RD_(sys)is below 0.9,local hydrodynamic jet initiated by void collapse dominates upon energy localization instead of interfacial impact.This study sheds light on novel insights for understanding the shock chemistry and physical-based atomic origin in crystalline explosives considering chemical-inclusions effects.展开更多
The gun-track launch system is a new special launch device that connects the track outside the muzzle.Because it is constrained by the track,the characteristics of development of the muzzle jet differ from those of th...The gun-track launch system is a new special launch device that connects the track outside the muzzle.Because it is constrained by the track,the characteristics of development of the muzzle jet differ from those of the traditional muzzle jet.Specifically,it changes from freely developing to doing so in a constrained manner,where this results in an asymmetric direction of flow as well as spatio-temporal coupling-induced interference between various shock waves and the formation of vortices.In this background,the authors of this article formulate and consider the development and characteristics of evolution of the muzzle jet as it impacts a constrained moving body.We designed simulations to test the gun-track launch system,and established a numerical model based on the dynamic grid method to explore the development and characteristics of propagation of disturbances when the muzzle jet impacted a constrained moving body.We also considered models without a constrained track for the sake of comparison.The results showed that the muzzle jet assumed a circumferential asymmetric shape,and tended to develop in the area above the muzzle.Because the test platform was close to the ground,the muzzle jet was subjected to reflections from it that enhanced the development and evolution of various forms of shock waves and vortices in the muzzle jet to exacerbate its rate of distortion and asymmetric characteristics.This in turn led to significant differences in the changes in pressure at symmetric points that would otherwise have been identical.The results of a comparative analysis showed that the constrained track could hinder the influence of reflections from the ground on the muzzle jet to some extent,and could reduce the velocity of the shock waves inducing the motion of the muzzle as well as the Mach number of the moving body.The work here provides a theoretical basis and the requisite technical support for applications of the gun-track launch system.It also sheds light on the technical bottlenecks that need to be considered to recover high-value warheads.展开更多
A critical challenge of any blast simulation facility is in producing the widest possible pressure-impulse range for matching against equivalent high-explosive events.Shock tubes and blast simulators are often constra...A critical challenge of any blast simulation facility is in producing the widest possible pressure-impulse range for matching against equivalent high-explosive events.Shock tubes and blast simulators are often constrained with the lack of effective ways to control blast wave profiles and as a result have a limited performance range.Some wave shaping techniques employed in some facilities are reviewed but often necessitate extensive geometric modifications,inadvertently cause flow anomalies,and/or are only applicable under very specific configurations.This paper investigates controlled venting as an expedient way for waveforms to be tuned without requiring extensive modifications to the driver or existing geometry and could be widely applied by existing and future blast simulation and shock tube facilities.The use of controlled venting is demonstrated experimentally using the Advanced Blast Simulator(shock tube)at the Australian National Facility of Physical Blast Simulation and via numerical flow simulations with Computational Fluid Dynamics.Controlled venting is determined as an effective method for mitigating the impact of re-reflected waves within the blast simulator.This control method also allows for the adjustment of parameters such as tuning the peak overpressure,the positive phase duration,and modifying the magnitude of the negative phase and the secondary shock of the blast waves.This paper is concluded with an illustration of the potential expanded performance range of the Australian blast simulation facility when controlled venting for blast waveform tailoring as presented in this paper is applied.展开更多
Products are often subject to dynamic environmental conditions in field use.When stress transition occurs,products may be exposed to instantaneous shocks that result in shock damages to the products,causing a permanen...Products are often subject to dynamic environmental conditions in field use.When stress transition occurs,products may be exposed to instantaneous shocks that result in shock damages to the products,causing a permanent change of the degradation signals.Meanwhile,under some conditions,instantaneous shocks also lead to stress drift,causing a temporary change of the degradation signals.In this paper,a degradation model is proposed to assess the reliability and predict the residual lifetime of products operating in a dynamic environment considering shock damage and stress drift.The model is established based on a Wiener process which combines a stress-dependent degradation rate function,a shock damage function and a stress drift function in response to the dynamic environment.The shock damage function is established as a linear function of the stress transition start level and the stress level increment.The stress drift function is established as the difference value of a specified function at the stress transition start and end levels.A simulation study is presented to demonstrate the application of the model,and a case study for miniature light bulbs is used to validate the effectiveness of the proposed model.展开更多
Dynamic strength behavior of Zr51Ti5NiloCu25A19 bulk metallic glass (BMG) up to 66 GPa was investigated in a series of plate impact shock-release and shock-reload experiments. Particle velocity profiles measured at ...Dynamic strength behavior of Zr51Ti5NiloCu25A19 bulk metallic glass (BMG) up to 66 GPa was investigated in a series of plate impact shock-release and shock-reload experiments. Particle velocity profiles measured at the sample/LiF window interface were used to estimate the shear stress, shear modulus, and yield stress in shocked BMG. Beyond confirm- ing the previously reported strain-softening of shear stress during the shock loading process for BMGs, it is also shown that the softened Zr-BMG still has a high shear modulus and can support large yield stress when released or reloaded from the shocked state, and both the shear modulus and the yield stress appear as strain-hardening behaviors. The work provides a much clearer picture of the strength behavior of BMGs under shock loading, which is useful to comprehensively understand the plastic deformation mechanisms of BMGs.展开更多
In the paper, the problems of dynamie response of multi-rigidbody systems with external impulsive forces are diseussed, and a set of equations with form of Lagrange method is presented. These equations are easy to be ...In the paper, the problems of dynamie response of multi-rigidbody systems with external impulsive forces are diseussed, and a set of equations with form of Lagrange method is presented. These equations are easy to be caleulated by computer.展开更多
The spall behavior of uranium is investigated using direct laser ablation loading experiments. The uranium targets are cut and ground to 0.05 mm, 0.1 mm, and 0.15 mm in thickness. Laser energies are varied to yield a ...The spall behavior of uranium is investigated using direct laser ablation loading experiments. The uranium targets are cut and ground to 0.05 mm, 0.1 mm, and 0.15 mm in thickness. Laser energies are varied to yield a constant peak pressure. This results in different strain rates and varying degrees of damage to the uranium targets. The spall strength is calculated and analyzed from the free surface velocity histories recorded using a line velocity interferometer for any reflections system. The spall strength increases from 4.3 GPa to 9.4 GPa with strain rates ranging from 4.0 × 10~6 s^(-1) to 1.7 × 10~7 s^(-1). Post-mortem analysis is performed on the recovered samples, revealing the twin-matrix interfaces together with the inclusions to be the primary factor governing the spall fracture of uranium.展开更多
With the rapid development of computer science and artificial intelligence technology, the complexity and intelligence of the neural network models constructed by people have been greatly improved. When the complex ne...With the rapid development of computer science and artificial intelligence technology, the complexity and intelligence of the neural network models constructed by people have been greatly improved. When the complex neuron system is subjected to the impact of "catastrophic", its original characteristics may be changed, and the consequences are difficult to predict. Catastrophe dynamics mainly studies the source of the sudden violent change of nature and human society and its evolution. The impact of the system can be divided into endogenous and exogenous shocks. In this article, catastrophe theory is used to study the neuron system. Based on the mean field model of Hurst and Sornette, introducing the weight parameters, mathematical models are constructed to study the response characteristics of the neuron system in face of exogenous shocks, endogenous shocks, and integrated shocks. The time characteristics of the shock response of the neuron system are discussed too, such as the instantaneous and long-term response of the system in face of shocks, the different response forms according to the weight or linear superposition, and the influence of adjusting parameters on the neuron system. The research result shows that the authoritarian coefficient and weight coefficient have a very important influence on the response of neuron system; By adjusting the two coefficients, the purpose of disaster prevention, self-healing protection and response reducing can be well achieved.展开更多
An intersecting cavern is a common structural form used in underground engineering,and its safety and stability performance directly control the service performance of the whole project.The dynamic responses of the th...An intersecting cavern is a common structural form used in underground engineering,and its safety and stability performance directly control the service performance of the whole project.The dynamic responses of the three kinds of crossing type(+-shaped,T-shaped,L-shaped)caverns subjected to ground shock were studied by numerical simulation.The velocity plus force mode boundary setting method was proposed in the coupled static and dynamic analysis of a deep underground cavern.The results show that,among the three types of crossing caverns,the+-shaped cavern is the most significantly affected by the dynamic action,followed by T-shaped,and then L-shaped caverns.The vault settlement,straight wall deformation,vault peak particle velocity,effective plastic strain of surrounding rock,and maximum principal stress and strain at the bottom of the lining of the straight wall increase with the increase of cavern span.The vault settlement,straight wall deformation,effective plastic strain of surrounding rock,and the maximum principal stress and strain at the bottom of lining to the straight wall decrease with the increase of lateral pressure coefficient,and the peak particle velocity at the vault increases.The variation is small compared with the change of cavern span.The influence range of the underground cavern intersection is two cavern diameters from the intersection centre.The bottom of the straight wall at the intersection is the weak part.It is suggested to thicken the support locally to improve the stability of the cavern.展开更多
文摘This paper presents a simple solution of the dynamic buckling of stiffened plates under in-plane impact loading. Based on large deflection theory, a discretely stiffened plate model has been used. The tangential stresses of stiffeners and in-plane displacement are neglected. Appling the Hamilton's principle, the motion equations of stiffened plates are obtained. The deflection of the plate is taken as Fourier series, and using Galerkin method the discrete equations can be deduced, which can be solved easily by Runge-Kutta method. The dynamic buckling loads of the stiffened plates are obtained form Budiansky-Roth criterion.
基金The 11th Postgraduate Technology Innovation Project of North University of China(No.20141147)
文摘In order to correct the test error caused by the dynamic characteristics of pressure sensor and avoid the influence of the error of sensor's dynamic model on compensation results,a dynamic compensation method of the pressure sensor is presented,which is based on quantum-behaved particle swarm optimization(QPSO)algorithm and the mean square error(MSE).By using this method,the inverse model of the sensor is built and optimized and then the coefficients of the optimal compensator are got.This method is verified by the dynamic calibration with shock tube and the dynamic characteristics of the sensor before and after compensation are analyzed in time domain and frequency domain.The results show that the working bandwidth of the sensor is extended effectively.This method can reduce dynamic measuring error and improve test accuracy in actual measurement experiments.
基金Xiaoding Shi was supported by National Natural Sciences Foundation of China(11471321)Yan Yong was supported by National Natural Sciences Foundation of China(11201301)
文摘In this paper, we study the vanishing viscosity limit for non-isentropic gas dy- namics with interacting shocks. Given any entropy solution of non-isentropic gas dynamics which consists of two different families of shocks interacting at some positive time, we show that such solution is the vanishing viscosity limit of a family of smooth global solutions for a viscous system of conservation law. We remark that, after the interacting time, not only shocks but also contact discontinuity are generated.
基金supported by a grant from 11th Five-Year Plan Medical Science Scientific Research Fund of the Chinese People's Liberation Army(08G002)
文摘BACKGROUND:Hypothermia is associated with poor outcome in trauma patients;however,hemorrhagic shock(HS)model with anesthetized swine was different from that of clinical reality.To identify the effects of environmental hypothermia on HS,we investigated hemodynamics and oxygen dynamics in an unanesthetized swine model of HS under simulating hypothermia environment.METHODS:Totally 16 Bama pigs were randomly divided into ambient temperature group(group A)and low temperature group(group B),8 pigs in each group.Venous blood(30 mL/kg)was continuously withdrawn for more than 15 minutes in conscious swine to establish a hemorrhagic shock model.Pulmonary arterial temperature(Tp),heart rate(HR),mean arterial pressure(MAP),pulmonary arterial pressure(PAP),pulmonary arterial wedge pressure(PAWP),central venous pressure(CVP),cardiac output(CO),hemoglobin(Hb),saturation of mixed venous blood(SvO_2)and blood gas analysis were recorded at the baseline and different hemorrhagic shock time(HST).The whole body oxygen delivery indices,DO_2l and VO_2l,and the O_2 extraction ratio(O_2ER)were calculated.RESULTS:Core body temperature in group A decreased slightly after the hemorrhagic shock model was established,and environmental hypothermia decreased in core body temperature.The mortality rate was significantly higher in group B(50%)than in group A(0%).DO_2l and VO_2l decreased significantly after hemorrhage.No difference was found in hemodynamics,DO_2l and VO_2l between group A and group B,but the difference in pH,lactic acid and O_2ER was significant between the two groups.CONCLUSION:Environmental hypothermia aggravated the disorder of oxygen metabolism after hemorrhagic shock,which was associated with poor prognosis.
文摘The power field of shock waves and fragments is analyzed and set up, and the damage modes of stiffened plates are put forward. According to the structural characters of the stiffened plates investigated and the properties of the shock waves and fragments, the experiments on the shock waves acting on the stiffened plates (penetrated and non-penetrated by fragments) are mainly conducted. The dynamic response rules of stiffened plates with holes under shock waves and fragments loading are obtained. The results show that the penetration of fragments into stiffened plates hardly affects their deformation produced by shock waves..
文摘An analytical method is developed to determine the transient response of dynamic thermostress in a two-layered cylinder with initial interface pressure. At first, the initial interface pressure in a two-layered cylinder caused by a heat-assembling method is considered as the initial condition of a thermal elastodynamic equilibrium equation. Thus, a thermal elastodynamic solution for a separate hollow cylinder with the initial stress field is found out by means of a series of simply mathematical transform. By making use of the boundary conditions and continuity conditions of a layered cylinders, the solution for the thermal shock exerting an influence on the initial interface pressure in a two-layered cylinder is also discussed.
基金supported by the National Natural Science Foundation of China(Grant No.11832006).
文摘A deep understanding of explosive sensitivities and their factors is important for safe and reliable applications.However,quantitative prediction of the sensitivities is difficult.Here,reactive molecular dynamics simulation models for high-speed piston impacts on explosive supercells were established.Simulations were also performed to investigate shock-induced reactions of various high-energy explosives.The fraction of reacted explosive molecules in an initial supercell changed linearly with the propagation distance of the shock-wave front.The corresponding slope could be used as a reaction rate for a specific shock-loading velocity.Reaction rates that varied with the shock-loading pressure exhibited two-stage linearities with different slopes.The two inflection points corresponded to the initial and accelerated reactions,which respectively correlated to the thresholds of shock-induced ignition and detonation.Therefore,the ignition and detonation critical pressures could be determined.The sensitivity could then be a quantitative prediction of the critical pressure.The accuracies of the quantitative shock sensitivity predictions were verified by comparing the impact and shock sensitivities of common explosives and the characteristics of anisotropic shock-induced reactions.Molecular dynamics simulations quantitatively predict and rank shock sensitivities by using only crystal structures of the explosives.Overall,this method will enable the design and safe use of explosives.
文摘A finite-difference algorithm is proposed for numerical modeling of hydrodynamic flows with rarefaction shocks, in which the fluid undergoes a jump-like liquid-gas phase transition. This new type of flow discontinuity, unexplored so far in computational fluid dynamics, arises in the approximation of phase-flip(PF) hydrodynamics, where a highly dynamic fluid is allowed to reach the innermost limit of metastability at the spinodal, upon which an instantaneous relaxation to the full phase equilibrium(EQ) is assumed. A new element in the proposed method is artificial kinetics of the phase transition, represented by an artificial relaxation term in the energy equation for a "hidden"component of the internal energy, temporarily withdrawn from the fluid at the moment of the PF transition. When combined with an appropriate variant of artificial viscosity in the Lagrangian framework, the latter ensures convergence to exact discontinuous solutions, which is demonstrated with several test cases.
基金Project supported by the National Magnetic Confinement Fusion Science Program of China(Grant No.2013GB109004)the National Natural Science Foundation of China(Grant Nos.51071095 and 50971077)
文摘Using molecular dynamics (MD) simulation, we study the thermal shock behavior of tungsten (W), which has been used for the plasma facing material (PFM) of tokamaks. The thermo-elastic stress wave, corresponding to the collective displacement of atoms, is analyzed with the Lagrangian atomic stress method, of which the reliability is also analyzed. The stress wave velocity corresponds to the speed of sound in the material, which is not dependent on the thermal shock energy. The peak pressure of a normal stress wave increases with the increase of thermal shock energy. We analyze the temperature evolution of the thermal shock region according to the Fourier transformation. It can be seen that the “obvious” velocity of heat propagation is less than the velocity of the stress wave; further, that the thermo-elastic stress wave may contribute little to the transport of kinetic energy. The heat propagation can be described properly by the heat conduction equation. These results may be useful for understanding the process of the thermal shock of tungsten.
基金the National Natural Science Foundation of China(U2241245,91960202 and 52271012)the National Key Laboratory Foundation of Science and Technology on Materials under Shock and Impact(6142902220301)+2 种基金the Aeronautical Science Foundation of China(2022Z053092001)the Shanghai Engineering Research Center of High-Performance Medical Device Materials(20DZ2255500)the Opening Project of National Key Laboratory of Shock Wave and Detonation Physics(2022JCJQLB05702).
文摘The mechanical response of a single crystal titanium sample against(0001)α surface impact was investigated using molecular dynamics simulation.Remarkably,non-uniform plastic deformation was observed in the sample.At high strain rates,amorphization occurred near the edge of the contact region where severe shear strain induced a large number of stacking faults(SFs)and dislocations.In contrast,the central part of the contact region underwent less deformation with significantly fewer dislocations.Moreover,instead of amorphization by consuming SFs and dislocations,there was a gradual increase in the density of dislocations and SFs during the process of amorphization.These local amorphous regions eventually grew into shear bands.
基金the financial support from National Natural Science Foundation of China(Grant Nos.11872119,12172051,and 11972329)Natural Science Foundation of Hubei Province(Grant No.2021CFB120)。
文摘Chemical inclusions significantly alter shock responses of crystalline explosives in macroscale gap experiments but their microscale dynamics origin remains unclear.Herein shock-induced energy localization,overall physical responses,and reactions in a-1,3,5-trinitro-1,3,5-triazinane(a-RDX)crystal entrained various chemical inclusions were investigated by the multi-scale shock technique implemented in the reactive molecular dynamics method.Results indicated that energy localization and shock reaction were affected by the intrinsic factors within chemical inclusions,i.e.,phase states,chemical compositions,and concentrations.The atomic origin of chemical-inclusions effects on energy localization is dependent on the dynamics mechanism of interfacial molecules with free space volume,which includes homogeneous intermolecular compression,interfacial impact and shear,and void collapse and jet.As introducing various chemical inclusions,the initiation of those dynamics mechanisms triggers diverse decay rates of bulk RDX molecules and hereby impacts on growth speeds of final reactions.Adding chemical inclusions can reduce the effectiveness of the void during the shock impacting.Under the shockwave velocity of 9 km/s,the parent RDX decay rate in RDX entrained amorphous carbon decreases the most and is about one fourth of that in RDX with a vacuum void,and solid HMX and TATB inclusions are more reactive than amorphous carbon but less reactive than dry air or acetone inclusions.The lessdense shocking system denotes the greater increases in local temperature and stress,the faster energy liberation,and the earlier final reaction into equilibrium,revealing more pronounced responses to the present intense shockwave.The quantitative models associated with the relative system density(RD_(sys))were proposed for indicating energy-localization mechanisms and evaluating initiation safety in the shocked crystalline explosive.RD_(sys)is defined by the density ratio of defective RDX to perfect crystal after dynamics relaxation and reveals the global density characteristic in shocked systems filled with chemical inclusions.When RD_(sys)is below 0.9,local hydrodynamic jet initiated by void collapse dominates upon energy localization instead of interfacial impact.This study sheds light on novel insights for understanding the shock chemistry and physical-based atomic origin in crystalline explosives considering chemical-inclusions effects.
文摘The gun-track launch system is a new special launch device that connects the track outside the muzzle.Because it is constrained by the track,the characteristics of development of the muzzle jet differ from those of the traditional muzzle jet.Specifically,it changes from freely developing to doing so in a constrained manner,where this results in an asymmetric direction of flow as well as spatio-temporal coupling-induced interference between various shock waves and the formation of vortices.In this background,the authors of this article formulate and consider the development and characteristics of evolution of the muzzle jet as it impacts a constrained moving body.We designed simulations to test the gun-track launch system,and established a numerical model based on the dynamic grid method to explore the development and characteristics of propagation of disturbances when the muzzle jet impacted a constrained moving body.We also considered models without a constrained track for the sake of comparison.The results showed that the muzzle jet assumed a circumferential asymmetric shape,and tended to develop in the area above the muzzle.Because the test platform was close to the ground,the muzzle jet was subjected to reflections from it that enhanced the development and evolution of various forms of shock waves and vortices in the muzzle jet to exacerbate its rate of distortion and asymmetric characteristics.This in turn led to significant differences in the changes in pressure at symmetric points that would otherwise have been identical.The results of a comparative analysis showed that the constrained track could hinder the influence of reflections from the ground on the muzzle jet to some extent,and could reduce the velocity of the shock waves inducing the motion of the muzzle as well as the Mach number of the moving body.The work here provides a theoretical basis and the requisite technical support for applications of the gun-track launch system.It also sheds light on the technical bottlenecks that need to be considered to recover high-value warheads.
基金funded partially by the Australian Government through the Australian Research Council’s Linkage Infrastructure,Equipment and Facilities (LIEF)funding scheme (LE130100133)。
文摘A critical challenge of any blast simulation facility is in producing the widest possible pressure-impulse range for matching against equivalent high-explosive events.Shock tubes and blast simulators are often constrained with the lack of effective ways to control blast wave profiles and as a result have a limited performance range.Some wave shaping techniques employed in some facilities are reviewed but often necessitate extensive geometric modifications,inadvertently cause flow anomalies,and/or are only applicable under very specific configurations.This paper investigates controlled venting as an expedient way for waveforms to be tuned without requiring extensive modifications to the driver or existing geometry and could be widely applied by existing and future blast simulation and shock tube facilities.The use of controlled venting is demonstrated experimentally using the Advanced Blast Simulator(shock tube)at the Australian National Facility of Physical Blast Simulation and via numerical flow simulations with Computational Fluid Dynamics.Controlled venting is determined as an effective method for mitigating the impact of re-reflected waves within the blast simulator.This control method also allows for the adjustment of parameters such as tuning the peak overpressure,the positive phase duration,and modifying the magnitude of the negative phase and the secondary shock of the blast waves.This paper is concluded with an illustration of the potential expanded performance range of the Australian blast simulation facility when controlled venting for blast waveform tailoring as presented in this paper is applied.
基金supported by the National Natural Science Foundation of China(NSFC71601009)the Technical Foundation Program from the Ministry of Industry and Information Technology of China(JSZL2015601B010)
文摘Products are often subject to dynamic environmental conditions in field use.When stress transition occurs,products may be exposed to instantaneous shocks that result in shock damages to the products,causing a permanent change of the degradation signals.Meanwhile,under some conditions,instantaneous shocks also lead to stress drift,causing a temporary change of the degradation signals.In this paper,a degradation model is proposed to assess the reliability and predict the residual lifetime of products operating in a dynamic environment considering shock damage and stress drift.The model is established based on a Wiener process which combines a stress-dependent degradation rate function,a shock damage function and a stress drift function in response to the dynamic environment.The shock damage function is established as a linear function of the stress transition start level and the stress level increment.The stress drift function is established as the difference value of a specified function at the stress transition start and end levels.A simulation study is presented to demonstrate the application of the model,and a case study for miniature light bulbs is used to validate the effectiveness of the proposed model.
基金Project supported by the National Natural Science Foundation of China(Grant No.11172281)
文摘Dynamic strength behavior of Zr51Ti5NiloCu25A19 bulk metallic glass (BMG) up to 66 GPa was investigated in a series of plate impact shock-release and shock-reload experiments. Particle velocity profiles measured at the sample/LiF window interface were used to estimate the shear stress, shear modulus, and yield stress in shocked BMG. Beyond confirm- ing the previously reported strain-softening of shear stress during the shock loading process for BMGs, it is also shown that the softened Zr-BMG still has a high shear modulus and can support large yield stress when released or reloaded from the shocked state, and both the shear modulus and the yield stress appear as strain-hardening behaviors. The work provides a much clearer picture of the strength behavior of BMGs under shock loading, which is useful to comprehensively understand the plastic deformation mechanisms of BMGs.
文摘In the paper, the problems of dynamie response of multi-rigidbody systems with external impulsive forces are diseussed, and a set of equations with form of Lagrange method is presented. These equations are easy to be caleulated by computer.
基金Supported by the Science Foundation of China Academy of Engineering Physics under Grant No A090504
文摘The spall behavior of uranium is investigated using direct laser ablation loading experiments. The uranium targets are cut and ground to 0.05 mm, 0.1 mm, and 0.15 mm in thickness. Laser energies are varied to yield a constant peak pressure. This results in different strain rates and varying degrees of damage to the uranium targets. The spall strength is calculated and analyzed from the free surface velocity histories recorded using a line velocity interferometer for any reflections system. The spall strength increases from 4.3 GPa to 9.4 GPa with strain rates ranging from 4.0 × 10~6 s^(-1) to 1.7 × 10~7 s^(-1). Post-mortem analysis is performed on the recovered samples, revealing the twin-matrix interfaces together with the inclusions to be the primary factor governing the spall fracture of uranium.
基金Project(CX2016B142)supported by the Hunan Provincial Innovation Foundation for Postgraduate,China
文摘With the rapid development of computer science and artificial intelligence technology, the complexity and intelligence of the neural network models constructed by people have been greatly improved. When the complex neuron system is subjected to the impact of "catastrophic", its original characteristics may be changed, and the consequences are difficult to predict. Catastrophe dynamics mainly studies the source of the sudden violent change of nature and human society and its evolution. The impact of the system can be divided into endogenous and exogenous shocks. In this article, catastrophe theory is used to study the neuron system. Based on the mean field model of Hurst and Sornette, introducing the weight parameters, mathematical models are constructed to study the response characteristics of the neuron system in face of exogenous shocks, endogenous shocks, and integrated shocks. The time characteristics of the shock response of the neuron system are discussed too, such as the instantaneous and long-term response of the system in face of shocks, the different response forms according to the weight or linear superposition, and the influence of adjusting parameters on the neuron system. The research result shows that the authoritarian coefficient and weight coefficient have a very important influence on the response of neuron system; By adjusting the two coefficients, the purpose of disaster prevention, self-healing protection and response reducing can be well achieved.
基金This research is funded by the National Natural Science Foundation of China(No.51478469).
文摘An intersecting cavern is a common structural form used in underground engineering,and its safety and stability performance directly control the service performance of the whole project.The dynamic responses of the three kinds of crossing type(+-shaped,T-shaped,L-shaped)caverns subjected to ground shock were studied by numerical simulation.The velocity plus force mode boundary setting method was proposed in the coupled static and dynamic analysis of a deep underground cavern.The results show that,among the three types of crossing caverns,the+-shaped cavern is the most significantly affected by the dynamic action,followed by T-shaped,and then L-shaped caverns.The vault settlement,straight wall deformation,vault peak particle velocity,effective plastic strain of surrounding rock,and maximum principal stress and strain at the bottom of the lining of the straight wall increase with the increase of cavern span.The vault settlement,straight wall deformation,effective plastic strain of surrounding rock,and the maximum principal stress and strain at the bottom of lining to the straight wall decrease with the increase of lateral pressure coefficient,and the peak particle velocity at the vault increases.The variation is small compared with the change of cavern span.The influence range of the underground cavern intersection is two cavern diameters from the intersection centre.The bottom of the straight wall at the intersection is the weak part.It is suggested to thicken the support locally to improve the stability of the cavern.