Nonlinear wave dispersion in monoatomic chains with lumped and distributed masses:discrete and continuum models

The main objective of this paper is to investigate the influence of inertia of nonlinear springs on the dispersion behavior of discrete monoatomic chains with lumped and distributed masses.The developed model can represent the wave propagation problem in a non-homogeneous material consisting of heavy inclusions embedded in a matrix.The inclusions are idealized by lumped masses,and the matrix between adjacent inclusions is modeled by a nonlinear spring with distributed masses.Additionally,the model is capable of depicting the wave propagation in bi-material bars,wherein the first material is represented by a rigid particle and the second one is represented by a nonlinear spring with distributed masses.The discrete model of the nonlinear monoatomic chain with lumped and distributed masses is first considered,and a closed-form expression of the dispersion relation is obtained by the second-order Lindstedt-Poincare method(LPM).Next,a continuum model for the nonlinear monoatomic chain is derived directly from its discrete lattice model by a suitable continualization technique.The subsequent use of the second-order method of multiple scales(MMS)facilitates the derivation of the corresponding nonlinear dispersion relation in a closed form.The novelties of the present study consist of(i)considering the inertia of nonlinear springs on the dispersion behavior of the discrete mass-spring chains;(ii)developing the second-order LPM for the wave propagation in the discrete chains;and(iii)deriving a continuum model for the nonlinear monoatomic chains with lumped and distributed masses.Finally,a parametric study is conducted to examine the effects of the design parameters and the distributed spring mass on the nonlinear dispersion relations and phase velocities obtained from both the discrete and continuum models.These parameters include the ratio of the spring mass to the lumped mass,the nonlinear stiffness coefficient of the spring,and the wave amplitude.

Discrete Choice Models and Artificial Intelligence Techniques for Predicting the Determinants of Transport Mode Choice--A Systematic Review

Forecasting travel demand requires a grasp of individual decision-making behavior.However,transport mode choice(TMC)is determined by personal and contextual factors that vary from person to person.Numerous characteristics have a substantial impact on travel behavior(TB),which makes it important to take into account while studying transport options.Traditional statistical techniques frequently presume linear correlations,but real-world data rarely follows these presumptions,which may make it harder to grasp the complex interactions.Thorough systematic review was conducted to examine how machine learning(ML)approaches might successfully capture nonlinear correlations that conventional methods may ignore to overcome such challenges.An in-depth analysis of discrete choice models(DCM)and several ML algorithms,datasets,model validation strategies,and tuning techniques employed in previous research is carried out in the present study.Besides,the current review also summarizes DCM and ML models to predict TMC and recognize the determinants of TB in an urban area for different transport modes.The two primary goals of our study are to establish the present conceptual frameworks for the factors influencing the TMC for daily activities and to pinpoint methodological issues and limitations in previous research.With a total of 39 studies,our findings shed important light on the significance of considering factors that influence the TMC.The adjusted kernel algorithms and hyperparameter-optimized ML algorithms outperform the typical ML algorithms.RF(random forest),SVM(support vector machine),ANN(artificial neural network),and interpretable ML algorithms are the most widely used ML algorithms for the prediction of TMC where RF achieved an R2 of 0.95 and SVM achieved an accuracy of 93.18%;however,the adjusted kernel enhanced the accuracy of SVM 99.81%which shows that the interpretable algorithms outperformed the typical algorithms.The sensitivity analysis indicates that the most significant parameters influencing TMC are the age,total trip time,and the number of drivers.

Qualitative Properties of Equilibrium Point in a Discrete Predator-Prey Model

In this paper, the dynamic properties of a discrete predator-prey model are discussed. The properties of non-hyperbolic fixed points and hyperbolic fixed points of the model are analyzed. First, by using the classic Shengjin formula, we find the existence conditions for fixed points of the model. Then, by using the qualitative theory of ordinary differential equations and matrix theory we indicate which points are hyperbolic and which are non-hyperbolic and the associated conditions.

Algorithmic approach to discrete fracture network flow modeling in consideration of realistic connections in large-scale fracture networks

Analyzing rock mass seepage using the discrete fracture network(DFN)flow model poses challenges when dealing with complex fracture networks.This paper presents a novel DFN flow model that incorporates the actual connections of large-scale fractures.Notably,this model efficiently manages over 20,000 fractures without necessitating adjustments to the DFN geometry.All geometric analyses,such as identifying connected fractures,dividing the two-dimensional domain into closed loops,triangulating arbitrary loops,and refining triangular elements,are fully automated.The analysis processes are comprehensively introduced,and core algorithms,along with their pseudo-codes,are outlined and explained to assist readers in their programming endeavors.The accuracy of geometric analyses is validated through topological graphs representing the connection relationships between fractures.In practical application,the proposed model is employed to assess the water-sealing effectiveness of an underground storage cavern project.The analysis results indicate that the existing design scheme can effectively prevent the stored oil from leaking in the presence of both dense and sparse fractures.Furthermore,following extensive modification and optimization,the scale and precision of model computation suggest that the proposed model and developed codes can meet the requirements of engineering applications.

Application of discrete choice model in trip mode structure forecast:a case study of Bengbu

In order to find the main factors that influence the urban traffic structure,a relational model between the travelers＇ characteristics and the trip mode choice is built.The data of urban residents＇ characteristics are obtained from statistical data,while the trip mode split data is collected through a trip survey in Bengbu.In addition,the discrete choice model is adopted to build the functional relationship between the mode choice and the travelers＇ personal characteristics,as well as family characteristics and trip characteristics.The model shows that the relationship between the mode split and the personal,as well as family and trip characteristics is stable and changes little as the time changes.Deduced by the discrete model,the mode split result is relatively accurate and can be feasibly used for trip mode structure forecasts.Furthermore,the proposed model can also contribute to find the key influencing factors on trip mode choice,and restructure or optimize the urban trip mode structure.

Discrete element modeling of migration and evolution rules of coarse aggregates in the static compaction process

To investigate migration and evolution rules of coarse aggregates in the static compaction process, an algorithm of generating digital coarse aggregates that can reflect real morphology（ such as shape, size and fracture surface） of aggregate particles, is represented by polyhedral particles based on the discrete element method（ DEM）. A digital specimen comprised of aggregates and air voids is developed. In addition,a static compaction model consisting of a digital specimen and three plates is constructed and a series of evaluation indices such as mean contact force σMCF, wall stress in direction of zcoordinate σWSZZ, porosity and coordination numbers are presented to investigate the motion rules of coarse aggregates at different compaction displacements of 7. 5, 15 and 30 mm. The three-dimensional static compaction model is also verified with laboratory measurements. The results indicate that the compaction displacements are positively related to σMCF and σWSZZ, which increase gradually with the increase in iterative steps. When the compaction proceeds, the digital specimen porosity decreases, but the coordination number increases. The variation ranges of these four indices are different at different compaction displacements. This study provides a method to analyze the compaction mechanism of particle materials such as asphalt mixture and graded broken stone.

Evaluating the mechanical properties of anisotropic shale containing bedding and natural fractures with discrete element modeling

Natural fracture data from one of the Carboniferous shale masses in the eastern Qaidam Basin were used to establish a stochastic model of a discrete fracture network and to perform discrete element simulation research on the size efect and mechanical parameters of shale.Analytical solutions of fctitious joints in transversely isotropic media were derived,which made it possible for the proposed numerical model to simulate the bedding and natural fractures in shale masses.The results indicate that there are two main factors infuencing the representative elementary volume(REV)size of a shale mass.The frst and most decisive factor is the presence of natural fractures in the block itself.The second is the anisotropy ratio:the greater the anisotropy is,the larger the REV.The bedding angle has little infuence on the REV size,whereas it has a certain infuence on the mechanical parameters of the rock mass.When the bedding angle approaches the average orientation of the natural fractures,the mechanical parameters of the shale blocks decrease greatly.The REV representing the mechanical properties of the Carboniferous shale masses in the eastern Qaidam Basin were comprehensively identifed by considering the infuence of bedding and natural fractures.When the numerical model size is larger than the REV,the fractured rock mass discontinuities can be transformed into equivalent continuities,which provides a method for simulating shale with natural fractures and bedding to analyze the stability of a borehole wall in shale.

IN-PLANE WAVE MOTION IN FINITE ELEMENT MODEL

The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV waves in the discrete model are first obtained by means of separating the characteristic equation of the motion equation, and then used to analyse the properties of P-and SV-homogeneous, inhomogeneous waves and other types of motion in the model. The dispersion characters, cut-off frequencies of P and SV waves, the polarization drift and appendent anisotropic property of wave motion caused by the discretization are finally discussed.

Interval grey number sequence prediction by using non-homogenous exponential discrete grey forecasting model

This paper aims to study a new grey prediction approach and its solution for forecasting the main system variable whose accurate value could not be collected while the potential value set could be defined. Based on the traditional nonhomogenous discrete grey forecasting model（NDGM）, the interval grey number and its algebra operations are redefined and combined with the NDGM model to construct a new interval grey number sequence prediction approach. The solving principle of the model is analyzed, the new accuracy evaluation indices, i.e. mean absolute percentage error of mean value sequence（MAPEM） and mean percent of interval sequence simulating value set covered（MPSVSC）, are defined and, the procedure of the interval grey number sequence based the NDGM（IG-NDGM） is given out. Finally, a numerical case is used to test the modelling accuracy of the proposed model. Results show that the proposed approach could solve the interval grey number sequence prediction problem and it is much better than the traditional DGM（1,1） model and GM（1,1） model.

A modified discrete element model for sea ice dynamics

Considering the discontinuous characteristics of sea ice on various scales,a modified discrete element model（DEM） for sea ice dynamics is developed based on the granular material rheology.In this modified DEM,a soft sea ice particle element is introduced as a self-adjustive particle size function.Each ice particle can be treated as an assembly of ice floes,with its concentration and thickness changing to variable sizes under the conservation of mass.In this model,the contact forces among ice particles are calculated using a viscous-elastic-plastic model,while the maximum shear forces are described with the Mohr-Coulomb friction law.With this modified DEM,the ice flow dynamics is simulated under the drags of wind and current in a channel of various widths.The thicknesses,concentrations and velocities of ice particles are obtained,and then reasonable dynamic process is analyzed.The sea ice dynamic process is also simulated in a vortex wind field.Taking the influence of thermodynamics into account,this modified DEM will be improved in the future work.

Research on cubic polynomial acceleration and deceleration control model for high speed NC machining

To satisfy the need of high speed NC (numerical control) machining, an acceleration and deceleration (acc/dec) control model is proposed, and the speed curve is also constructed by the cubic polynomial. The proposed control model provides continuity of acceleration, which avoids the intense vibration in high speed NC machining. Based on the discrete characteristic of the data sampling interpolation, the acc/dec control discrete mathematical model is also set up and the discrete expression of the theoretical deceleration length is obtained furthermore. Aiming at the question of hardly predetermining the deceleration point in acc/dec control before interpolation, the adaptive acc/dec control algorithm is deduced from the expressions of the theoretical deceleration length. The experimental result proves that the acc/dec control model has the characteristic of easy implementation, stable movement and low impact. The model has been applied in multi-axes high speed micro fabrication machining successfully.

Very Large Eddy Simulation of Cavitation from Inception to Sheet/Cloud Regimes by A Multiscale Model

The cavitating flow in different regimes has the intricate flow structure with multiple time and space scales.The present work develops a multiscale model by coupling the volume of fluid(VOF)method and a discrete bubble model(DBM),to simulate the cavitating flow in a convergent-divergent test section.The Schnerr-Sauer cavitation model is used to calculate the mass transfer rate to obtain the macroscale phase structure,and the simplified Rayleigh-Plesset equation is applied to simulate the growing and collapsing of discrete bubbles.An algorithm for bridging between the macroscale cavities and microscale bubbles is also developed to achieve the multiscale simulation.For the flow field,the very large eddy simulation(VLES)approach is applied.Conditions from inception to sheet/cloud cavitation regimes are taken into account and simulations are conducted.Compared with the experimental observations,it is shown that the cavitation inception,bubble clouds formation and glass cavity generation are all well represented,indicating that the proposed VOF-DBM model is a promising approach to accurately and comprehensively reveal the multiscale phase field induced by cavitation.

Mathematical modelling of translational motion of rail-guided cart with suspended payload

In this paper modelling of the translational motion of transportation rail-guided cart with rope suspended payload is considered. The linearly moving cart,driven by a travel mechanism,is modelled as a discrete six degrees of freedom (DOF) dynamic system. The hoisting mechanism for lowering and lifting the payload is considered and is included in the dynamic model as one DOF system. Differential equations of motion of the cart elements are derived using Lagrangian dynamics and are solved for a set of real-life constant parameters of the cart. A two-sided interaction was observed between the swinging payload and the travel mechanism. Results for kinematical and force parameters of the system are obtained. A verification of the proposed model was conducted.

Research on SINS Error Recursive Model of Ballistic Missile

The strap-down inertial navigation system (SINS) error of ballistic missile is generated by the mutual influence of gyroscope and accelerometer, and the recursive model is completely different from that of gimbaled INS. In the paper, a discrete error recursive model was obtained by studying the applied SINS error model of ballistic missile, and the discrete Kalman filtering simulation based on the model was carried out. The simulated results show that the model can depict the SINS error exactly and provide the advantages for research on integrated guidance and improved hit accuracy.

An Integrated Dynamic Model of Ocean Mining System and Fast Simulation of Its Longitudinal Reciprocating Motion

An integrated dynamic model of China＇s deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.

Dynamic Analysis of the Seafloor Pilot Miner Based on Single-Body Vehicle Model and Discretized Track-Terrain Interaction Model

In order to achieve the complex dynamic analysis of the self-propelled seafloor pilot miner moving on the seafloor of extremely cohesive soft soil and further to make it possible to integrate the miner system with some subsystems to form the complete integrated deep ocean mining pilot system and perform dynamic analysis, a new method for the dynamic modeling and analysis of the miner is proposed and developed in this paper, resulting in a simplified 3D single-body vehicle model with three translational and three rotational degrees of freedom, while the track-terrain interaction model is built by partitioning the track-terrain interface into discrete elements with parameterized force dements built on the theory of terramechanics acting on each discrete dement. To evaluate and verify the correctness and effectiveness of this new modeling and analysis method, typical comparative studies with regard to computational efficiency and solution accuracy are carried out between the traditional modeling method of building the tracked vehicle as a multi-body model and the new modeling method. In full consideration of the particMar structure design of the pilot miner, the special characteristics of the seafioor soil and the hydrodynamic force of near-seafloor currnt, the dynamic simulation analysis of the miner is performed and discussed, which can provide useful guidance and reference for the practical miner system in design and operation. This new method can not only realize the rapid dynamic simulation analysis of the miner but also make possible the integration and rapid dynamic analysis of the complete integrated deep ocean mining pilot system in further researches.

DPM simulation in an underground entry: Comparison between particle and species models

The diesel particulate matter(DPM) emission from diesel powered equipment in underground mines can cause health hazards including cancer to the miners. The understanding of the DPM propagation pattern under realistic mining condition is required for selecting proper DPM control strategies and to improve working practices in underground mines. In this paper, three dimensional simulations of DPM emission from the exhaust tail pipe of a load-haul-dump(LHD) vehicle and its subsequent distribution inside an isolated zone in the typical underground mine are carried out using two different solution models available in Ansys Fluent. The incoming fresh air into the isolated zone is treated as a continuous phase and DPM is treated either as a continuous phase(gas) or as a secondary discrete phase(particle). Species transport model is used when DPM is treated as gas and discrete phase model is used when DPM is assumed to behave like a particle. The distributions of DPM concentration inside the isolated zone obtained from each method are presented and compared. From the comparison results, an accurate and economical solution technique for DPM evaluation can be selected.

Discrete element modeling of acoustic emission in rock fracture

The acoustic emission （AE） features in rock fracture are simulated numerically with discrete element model （DEM）. The specimen is constructed by using spherical particles bonded via the parallel bond model. As a result of the heterogeneity in rock specimen, the failure criterion of bonded particle is coupled by the shear and tensile strengths, which follow a normal probability distribution. The Kaiser effect is simulated in the fracture process, for a cubic rock specimen under uniaxial compression with a constant rate. The AE number is estimated with breakages of bonded particles using a pair of parameters, in the temporal and spatial scale, respectively. It is found that the AE numbers and the elastic energy release curves coincide. The range for the Kaiser effect from the AE number and the elastic energy release are the same. Furthermore, the frequency-magnitude relation of the AE number shows that the value of B determined with DEM is consistent with the experimental data.

A review of discrete modeling techniques for fracturing processes in discontinuous rock masses

The goal of this review paper is to provide a summary of selected discrete element and hybrid finitediscrete element modeling techniques that have emerged in the field of rock mechanics as simulation tools for fracturing processes in rocks and rock masses. The fundamental principles of each computer code are illustrated with particular emphasis on the approach specifically adopted to simulate fracture nucleation and propagation and to account for the presence of rock mass discontinuities. This description is accompanied by a brief review of application studies focusing on laboratory-scale models of rock failure processes and on the simulation of damage development around underground excavations.

Isotropic compression behavior of granular assembly with non-spherical particles by X-ray micro-computed tomography and discrete element modeling

The particle morphological properties,such as sphericity,concavity and convexity,of a granular assembly significantly affect its macroscopic and microscopic compressive behaviors under isotropic loading condition.However,limited studies on investigating the microscopic behavior of the granular assembly with real particle shapes under isotropic compression were reported.In this study,X-ray computed tomography(mCT)and discrete element modeling(DEM)were utilized to investigate isotropic compression behavior of the granular assembly with regard to the particle morphological properties,such as particle sphericity,concavity and interparticle frictions.The mCT was first used to extract the particle morphological parameters and then the DEM was utilized to numerically investigate the influences of the particle morphological properties on the isotropic compression behavior.The image reconstruction from mCT images indicated that the presented particle quantification algorithm was robust,and the presented microscopic analysis via the DEM simulation demonstrated that the particle surface concavity significantly affected the isotropic compression behavior.The observations of the particle connectivity and local void ratio distribution also provided insights into the granular assembly under isotropic compression.Results found that the particle concavity and interparticle friction influenced the most of the isotropic compression behavior of the granular assemblies.