The impact of demographic dynamics on food consumption and its environmental outcomes:Evidence from China

With increasing population and changing demographics,food consumption has experienced a significant transition in quantity and quality.However,a dearth of knowledge remains regarding its environmental impacts and how it responds to demographic dynamics,particularly in emerging economies like China.Using the two-stage Quadratic Almost Demand System(QUAIDS)model,this study empirically examines the impact of demographic dynamics on food consumption and its environmental outcomes based on the provincial data from 2000 to 2020 in China.Under various scenarios,according to changes in demographics,we extend our analysis to project the long-term trend of food consumption and its environmental impacts,including greenhouse gas(GHG)emissions,water footprint(WF),and land appropriation(LA).The results reveal that an increase in the proportion of senior people significantly decreases the consumption of grain and livestock meat and increases the consumption of poultry,egg,and aquatic products,particularly for urban residents.Moreover,an increase in the proportion of males in the population leads to higher consumption of poultry and aquatic products.Correspondingly,in the current scenario of an increased aging population and sex ratio,it is anticipated that GHG emissions,WF,and LA are likely to decrease by 1.37,2.52,and 3.56%,respectively.More importantly,in the scenario adhering to the standards of nutritional intake according to the Dietary Guidelines for Chinese Residents in 2022,GHG emissions,WF,and LA in urban areas would increase by 12.78,20.94,and 18.32%,respectively.Our findings suggest that changing demographics should be considered when designing policies to mitigate the diet-environment-health trilemma and achieve sustainable food consumption.

Investigation of Projectile Impact Behaviors of Graphene Aerogel Using Molecular Dynamics Simulations

Graphene aerogel(GA),as a novel solid material,has shown great potential in engineering applications due to its unique mechanical properties.In this study,the mechanical performance of GA under high-velocity projectile impacts is thoroughly investigated using full-atomic molecular dynamics(MD)simulations.The study results show that the porous structure and density are key factors determining the mechanical response of GA under impact loading.Specifically,the impact-induced penetration of the projectile leads to the collapse of the pore structure,causing stretching and subsequent rupture of covalent bonds in graphene sheets.Moreover,the effects of temperature on the mechanical performance of GA have been proven to be minimal,thereby highlighting the mechanical stability of GA over a wide range of temperatures.Finally,the energy absorption density(EAD)and energy absorption efficiency(EAE)metrics are adopted to assess the energy absorption capacity of GA during projectile penetration.The research findings of this work demonstrate the significant potential of GA for energy absorption applications.

Study on dynamic response of multi-degree-of-freedom explosion vessel system under impact load

In order to study the dynamic response and calculate the axial dynamic coefficient of the monolayer cylindrical explosion vessel,the wall of vessel is simplified as a multi-degree-of-freedom(MDoF) undamped elastic foundation beam.Decoupling the coupled motion equation and using Duhamel's integrals,the solutions in generalized coordinates of the equations under exponentially decaying loads,square wave loads and triangular wave loads are calculated.These solutions are consistent in form with the solutions of single-degree-of-freedom(SDoF) undamped forced vibration simplified model.Based on the model,equivalent MDoF design method(also called MDoF dynamic coefficient method) of cylindrical explosion vessel is proposed.The traditional method can only predict the dynamic coefficient of torus portion around the explosion center,but this method can predict that of the vessel wall at any axial n dividing point position.It is verified that the prediction accuracy of this model is greatly improved compared with the SDoF model by comparing the results of this model with SDoF model and numerical simulation in different working conditions.However,the prediction accuracy decreases as the scaled distance decreases and approaches the end of the vessel,which is related to the accuracy of the empirical formula of the implosion load,the simplification of the explosion load direction,the boundary conditions,and the loading time difference.

Numerical simulation of dynamic response of operating metro tunnel induced by ground explosion

To evaluate the effects of possible ground explosion on a shallow-buried metro tunnel, this paper attempts to analyze the dynamic responses of the operating metro tunnel in soft soil, using a widely applied explicit dynamic nonlinear finite element software ANSYS/LS-DYNA. The blast induced wave propagation in the soil and the tunnel, and the von Mises effective stress and acceleration of the tunnel lining were presented, and the safety of the tunnel lining was evaluated based on the failure criterion. Besides, the parametric study of the soil was also carried out. The numerical results indicate that the upper part of the tunnel lining cross-section with directions ranging from 0° to 22.5° and horizontal distances 0 to 7 m away from the explosive center are the vulnerable areas, and the metro tunnel might be safe when tunnel depth is more than 7 m and TNT charge on the ground is no more than 500 kg, and the selection of soil parameters should be paid more attentions to conduct a more precise analysis.

Partition method for impact dynamics of flexible multibody systems based on contact constraint

The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.

THE APPLICATION OF DISCRETE ELEMENT METHOD IN SOLVING THREE-DIMENTIONAL IMPACT DYNAMICS PROBLEMS

A three-dimensional discrete element model of the connective type is presented. Moreover,a three-dimensional numerical analysis code,which can carry out the transitional pro- cess from connective model(for continuum)to contact model(for non-continuum),is developed for simulating the mechanical process from continuum to non-continuum.The wave propagation process in a concrete block(as continuum)made of cement grout under impact loading is numer- ically simulated with this code.By comparing its numerical results with those by LS-DYNA,the calculation accuracy of the model and algorithm is proved.Furthermore,the failure process of the concrete block under quasi-static loading is demonstrated,showing the basic dynamic tran- sitional process from continuum to non-continuum.The results of calculation can be displayed by animation.The damage modes are similar to the experimental results.The two numerical examples above prove that our model and its code are powerful and efficient in simulating the dynamic failure problems accompanying the transition from continuum to non-continuum.It also shows that the discrete element method(DEM)will have broad prospects for development and application.

Dynamic responses of reinforced concrete beams under double-endinitiated close-in explosion

The reinforced concrete(RC) structural component might suffer a great damage under close-in explosion.Different from distant explosions, blast loads generated by the close-in explosion are non-uniformly distributed on the structural component and may cause both local and structural failure. In this study,an experimental study was conducted to investigate the dynamic responses of RC beams under doubleend-initiated close-in explosions. The experimental results show that the distribution of blast loads generated by the double-end-initiated explosion is much more non-uniform than those generated by single-point detonation, which is caused by the self-Mach-reflection effects. A 3 D finite element model was developed and validated in LS-DYNA by employing the modified K&C model. Intensive numerical calculations were conducted to study the influences of the initiation way, scaled distance and longitudinal reinforcement ratio on the dynamic responses and failure modes of RC beams. Numerical results show that the RC beam suffers greater damage as the cylindrical explosive is detonated at its double ends than the scenario in which the cylindrical explosive is detonated at its central point. RC beams mainly suffer flexural failure and flexure-shear failure under the double-end close-in explosion, and the failure modes of RC beams change from the flexural damage to flexure-shear damage as the scaled distance or the longitudinal reinforcement ratio decreases. The direct shear failure mode is not usually observed in the double-end-initiated explosion, since the intense blast loads is basically concentrated in the midspan of RC beam, which is due to self-Mach-reflection enhancement.

Impact dynamics analysis of free-floating space manipulator capturing satellite on orbit and robust adaptive compound control algorithm design for suppressing motion

The impact dynamics, impact effect, and post-impact unstable motion sup- pression of free-floating space manipulator capturing a satellite on orbit are analyzed. Firstly, the dynamics equation of free-floating space manipulator is derived using the sec- ond Lagrangian equation. Combining the momentum conservation principle, the impact dynamics and effect between the space manipulator end-effector and satellite of the cap- ture process are analyzed with the momentum impulse method. Focusing on the unstable motion of space manipulator due to the above impact effect, a robust adaptive compound control algorithm is designed to suppress the above unstable motion. There is no need to control the free-floating base position to save the jet fuel. Finally, the simulation is proposed to show the impact effect and verify the validity of the control algorithm.

Partition method and experimental validation for impact dynamics of flexible multibody system

The impact problem of a flexible multibody system is a non-smooth, high-transient, and strong-nonlinear dynamic process with variable boundary. How to model the contact/impact process accurately and efficiently is one of the main difficulties in many engineering applications. The numerical approaches being used widely in impact analysis are mainly from two fields： multibody system dynamics （MBS） and computational solid mechanics （CSM）. Approaches based on MBS provide a more efficient yet less accurate analysis of the contact/impact problems, while approaches based on CSM are well suited for particularly high accuracy needs, yet require very high computational effort. To bridge the gap between accuracy and efficiency in the dynamic simulation of a flexible multibody system with contacts/impacts, a partition method is presented considering that the contact body is divided into two parts, an impact region and a non-impact region. The impact region is modeled using the finite element method to guarantee the local accuracy, while the non-impact region is modeled using the modal reduction approach to raise the global efficiency. A three-dimensional rod-plate impact experiment is designed and performed to validate the numerical results. The principle for how to partition the contact bodies is proposed： the maximum radius of the impact region can be estimated by an analytical method, and the modal truncation orders of the non-impact region can be estimated by the highest frequency of the signal measured. The simulation results using the presented method are in good agreement with the experimental results. It shows that this method is an effec-rive formulation considering both accuracy and efficiency. Moreover, a more complicated multibody impact problem of a crank slider mechanism is investigated to strengthen this conclusion.

A novel approach to the dynamic response analysis of Euler-Bernoulli beams resting on a Winkler soil model and subjected to impact loads

This work presents a novel approach to the dynamic response analysis of a Euler-Bernoulli beam resting on a Winkler soil model and subjected to an impact loading.The approach considers that damping has much less importance in controlling the maximum response to impulsive loadings because the maximum response is reached in a very short time,before the damping forces can dissipate a significant portion of the energy input into the system.The development of two sine series solutions,relating to different types of impulsive loadings,one involving a single concentrated force and the other a distributed line load,are presented.This study revealed that when a simply supported Euler-Bernoulli beam,resting on a Winkler soil model,is subject to an impact load,the resulting vertical displacements,bending moments and shear forces produced along the span of the beam are considerably affected.In particular,the quantification of this effect is best observed,relative to the corresponding static solution,via an amplification factor.The computed impact amplification factors,for the sub-grade moduli used in this study,were in magnitude greater than 2,thus confirming the multiple-degree-of-freedom nature of the problem.

Numerical and experimental studies on impact dynamics of a planar flexible multibody system

In this paper a computational methodology on impact dynamics of the flexible multibody system is presented. First, the floating frame of reference approach and nodal coordinates on the basis of finite element formulation are used to describe the kinematics of planar deformable bodies. According to the kinematic description of contact conditions, the contact constraint equations of planar flexible bodies are derived. Based on the varying topology technique the impact dynamic equations for a planar multibody system are established. Then the initial conditions of the equations in each contact stage are determined according to the discontinuity theory in continuum mechanics. The experiments between the aluminum rods are performed to check the correctness of the proposed method. Through the comparison between the numerical and experimental results the proposed method is validated. Experimental results also show that the impulse momentum method cannot accurately predict the complex impact dynamic phenomena and the continuous model may lead to a serious error when used to simulate the impact problems with significant wave propagation effects.

Quantitative prediction and ranking of the shock sensitivity ofexplosives via reactive molecular dynamics simulations

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.

Missile guidance law design based on free-time convergent error dynamics

To solve the finite-time error-tracking problem in mis-sile guidance,this paper presents a unified design approach through error dynamics and free-time convergence theory.The proposed approach is initiated by establishing a desired model for free-time convergent error dynamics,characterized by its independence from initial conditions and guidance parameters,and adjustable convergence time.This foundation facilitates the derivation of specific guidance laws that integrate constraints such as leading angle,impact angle,and impact time.The theoretical framework of this study elucidates the nuances and synergies between the proposed guidance laws and existing methodologies.Empirical evaluations through simulation comparisons underscore the enhanced accuracy and adaptability of the proposed laws.

Micro defects formation and dynamic response analysis of steel plate of quasi-cracking area subjected to explosive load

As the protective component,steel plate had attracted extensive attention because of frequently threats of explosive loads.In this paper,the evolution of microstructure and the mechanism of damage in the quasi-cracking area of steel plate subjected to explosive load were discussed and the relationships between micro defects and dynamic mechanical response were revealed.After the explosion experiment,five observation points were selected equidistant from the quasi-cracking area of the section of the steel plate along the thickness direction,and the characteristics of micro defects at the observation points were analyzed by optical microscope(OM),scanning electron microscope(SEM) and electron backscattered diffraction(EBSD).The observation result shows that many slip bands(SBs) appeared,and the grain orientation changed obviously in the steel plate,the two were the main damage types of micro defects.In addition,cracks,peeling pits,grooves and other lager micro defects were appeared in the lower area of the plate.The stress parameters of the observation points were obtained through an effective numerical model.The mechanism of damage generation and crack propagation in the quasicracking area were clarified by comparing the specific impulse of each observation point with the corresponding micro defects.The result shows that the generation and expansion of micro defects are related to the stress area(i.e.the upper compression area,the neutral plane area,and the lower tension area).The micro defects gather and expand at the grain boundary,and will become macroscopic damage under the continuous action of tensile stress.Besides,the micro defects at the midpoint of the section of the steel plate in the direction away from the explosion center(i.e.the horizontal direction) were also studied.It was found that the specific impulse at these positions were much smaller than that in the thickness direction,the micro defects were only SBs and a few micro cracks,and the those decreased with the increase of the distance from the explosion center.

Evaluation of tidal stream energy and its impacts on surrounding dynamics in the Eastern Region of Pingtan Island, China

Using an improved FVCOM numerical model, combined with the momentum-sinking scheme based on the structural characteristics of specific turbines, this study analyzed the temporal and spatial distributions of tidal energy resources before and after the deployment of tidal turbines near Pingtan Island, China. Considering factors such as the distribution of tidal stream energy, bathymetry, topography, and the design parameters of the turbines, an appropriate location for a demonstration tidal turbine was selected and the corresponding energy resource was evaluated. Several sites with strong tidal streams were considered： south of the northern cape, east of the southem cape, and the southern end of Haitan Bay. The former was thought most suitable for the deployment of a tidal energy turbine, with projected power generation for approximately 470 h per month. The average power of this demonstration was about 2.4 kW, and the annual electricity output was approximately 17.47 MWh. The intervention of the turbine device had little influence on the near-field tidal stream or water level. The tidal stream was reduced slightly in the area south of the northern cape, although the effect weakened further from the turbine. Conversely, the velocity increased slightly on both sides of the demonstration site. The difference in current speed with and without the turbine was greater at slack tide than still tide. The influence of turbine operation on water level was minor. The method adopted in this study can be considered a reference for the selection of sites for the demonstration of tidal stream energy. However, the method is unable describe the dynamic characteristics of the turbulent flow surrounding the deployed turbines, which has an important role regarding the optimal designs of the turbine blade and pile foundations. Therefore, we will continue to work to improve this model in future research.

Dynamic impact properties of deep sandstone under thermal-hydraulicmechanical coupling loads

The deep rock mass within coal mines situated in a challenging environment are characterized by high ground stress,high geotemperature,high osmotic water pressure,and dynamic disturbances from mechanical excavation.To investigate the impact of this complex mechanical environment on the dynamic characteristics of roof sandstone in self-formed roadways without coal pillars,standard specimens of deep sandstone from the 2611 upper tunnel working face of the Yongmei Company within the Henan Coal Chemical Industry Group in Henan,China were prepared,and an orthogonal test was designed.Using a self-developed geotechnical dynamic impact mechanics test system,triaxial dynamic impact tests under thermal-hydraulicmechanical coupling conditions were conducted on deep sandstone.The results indicate that under high confining pressure,deep sandstone exhibits pronounced brittle failure at low temperatures,with peak strength gradually decreasing as temperature and osmotic water pressure increase.Conversely,under low confining pressure and low temperature,the brittleness of deep sandstone weakens gradually,while ductility increases.Moreover,sandstone demonstrates higher peak strength at low temperatures under high axial pressure conditions,lower peak strength at high temperatures,and greater strain under low axial pressure and high osmotic water pressure.Increases in impact air pressure and osmotic water pressure have proportionally greater effects on peak stress and peak strain.Approximately 50%of the input strain energy is utilized as effective energy driving the sandstone fracture process.Polar analysis identifies the optimal combination of factors affecting the peak stress and peak strain of sandstone.Under the coupling effect,intergranular and transgranular fractures occur within the sandstone.SEM images illustrate that the damage forms range from minor damage with multiple fissures to extensive fractures and severe fragmentation.This study elucidates the varied dynamic impact mechanical properties of deep sandstones under thermal-hydraulic-mechanical coupling,along with multifactor analysis methods and their optimal factor combinations.

Experimental investigation on dynamic response and damage models of circular RC columns subjected to underwater explosions

Reinforced concrete(RC) columns are widely used as supporting structures for high-piled wharfs.The study of damage model of a RC column due to underwater explosion is a critical issue to assess the wharfs antiknock security.In this study,the dynamic response and damage model of circular RC columns subjected to underwater explosions were investigated by means of scaled-down experiment models.Experiments were carried out in a 10.0 m diameter tank with the water depth of 2.25 m,under different explosive quantities(0.025 kg-1.6 kg),stand-off distances(0.0 m-7.0 m),and detonation depths(0.25 m-2.0 m).The shock wave load and dynamic response of experiment models were measured by configuring sensors of pressure,acceleration,strain,and displacement.Then,the load distribution characteristics,time history of test data,and damage models related to present conditions were obtained and discussed.Three damage models,including bending failure,bending-shear failure and punching failure,were identified.In addition,the experie nce model of shock wave loads on the surface of a RC column was proposed for engineering application.

An Analytical Solution for Dynamic Response of the Plate with Different Impedances Subjected to Underwater Explosion

In this paper, the transmitted part of the incident wave is considered to revise Taylor＇s solution, which is used to extend its application for analytical models to predict the response of the plate with different material properties. The influence of the material properties and the boundary condition of the plate on fluid and structural dynamics is systematically investigated. The analytical results are compared with those of detailed dynamic FE simulations and the two are in good agreement. The results indicate that the analytical method is valid and suitable for the plates with different material properties subjected to underwater explosion. It is found that Taylor＇s results of the plate with small impedance are invalid, which indicates a potential application field for the analytical method.

Microscopic defects formation and dynamic mechanical response analysis of Q345 steel plate subjected to explosive load

As the basic protective element, steel plate had attracted world-wide attention because of frequent threats of explosive loads. This paper reports the relationships between microscopic defects of Q345 steel plate under the explosive load and its macroscopic dynamics simulation. Firstly, the defect characteristics of the steel plate were investigated by stereoscopic microscope(SM) and scanning electron microscope(SEM). At the macroscopic level, the defect was the formation of cave which was concentrated in the range of 0-3.0 cm from the explosion center, while at the microscopic level, the cavity and void formation were the typical damage characteristics. It also explains that the difference in defect morphology at different positions was the combining results of high temperature and high pressure. Secondly, the variation rules of mechanical properties of steel plate under explosive load were studied. The Arbitrary Lagrange-Euler(ALE) algorithm and multi-material fluid-structure coupling method were used to simulate the explosion process of steel plate. The accuracy of the method was verified by comparing the deformation of the simulation results with the experimental results, the pressure and stress at different positions on the surface of the steel plate were obtained. The simulation results indicated that the critical pressure causing the plate defects may be approximately 2.01 GPa. On this basis, it was found that the variation rules of surface pressure and microscopic defect area of the Q345 steel plate were strikingly similar, and the corresponding mathematical relationship between them was established. Compared with Monomolecular growth fitting models(MGFM) and Logistic fitting models(LFM), the relationship can be better expressed by cubic polynomial fitting model(CPFM). This paper illustrated that the explosive defect characteristics of metal plate at the microscopic level can be explored by analyzing its macroscopic dynamic mechanical response.

Molecular Dynamics Simulation of Shock Response of CL-20 Co-crystals Containing Void Defects

To investigate the effect of void defects on the shock response of hexanitrohexaazaisowurtzitane(CL-20)co-crystals,shock responses of CL-20 co-crystals with energetic materials ligands trinitrotoluene(TNT),1,3-dinitrobenzene(DNB),solvents ligands dimethyl carbonate(DMC) and gamma-butyrolactone(GBL)with void were simulated,using molecular dynamics method and reactive force field.It is found that the CL-20 co-crystals with void defects will form hot spots when impacted,significantly affecting the decomposition of molecules around the void.The degree of molecular fragmentation is relatively low under the reflection velocity of 2 km/s,and the main reactions are the formation of dimer and the shedding of nitro groups.The existence of voids reduces the safety of CL-20 co-crystals,which induced the sensitivity of energetic co-crystals CL-20/TNT and CL-20/DNB to increase more significantly.Detonation has occurred under the reflection velocity of 4 km/s,energetic co-crystals are easier to polymerize than solvent co-crystals,and are not obviously affected by voids.The results show that the energy of the wave decreases after sweeping over the void,which reduces the chemical reaction frequency downstream of the void and affects the detonation performance,especially the solvent co-crystals.