A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Dynamic failure process of expanded polystyrene particle lightweight soil under cyclic loading using discrete element method

Expanded polystyrene(EPS)particle-based lightweight soil,which is a type of lightweight filler,is mainly used in road engineering.The stability of subgrades under dynamic loading is attracting increased research attention.The traditional method for studying the dynamic strength characteristics of soils is dynamic triaxial testing,and the discrete element simulation of lightweight soils under cyclic load has rarely been considered.To study the meso-mechanisms of the dynamic failure processes of EPS particle lightweight soils,a discrete element numerical model is established using the particle flow code(PFC)software.The contact force,displacement field,and velocity field of lightweight soil under different cumulative compressive strains are studied.The results show that the hysteresis curves of lightweight soil present characteristics of strain accumulation,which reflect the cyclic effects of the dynamic load.When the confining pressure increases,the contact force of the particles also increases.The confining pressure can restrain the motion of the particle system and increase the dynamic strength of the sample.When the confining pressure is held constant,an increase in compressive strain causes minimal change in the contact force between soil particles.However,the contact force between the EPS particles decreases,and their displacement direction points vertically toward the center of the sample.Under an increase in compressive strain,the velocity direction of the particle system changes from a random distribution and points vertically toward the center of the sample.When the compressive strain is 5%,the number of particles deflected in the particle velocity direction increases significantly,and the cumulative rate of deformation in the lightweight soil accelerates.Therefore,it is feasible to use 5%compressive strain as the dynamic strength standard for lightweight soil.Discrete element methods provide a new approach toward the dynamic performance evaluation of lightweight soil subgrades.

Dynamic Location Method for Shallow Ocean Bottom Nodes Using the Levenberg-Marquart Algorithm

Ocean bottom node(OBN)data acquisition is the main development direction of marine seismic exploration;it is widely promoted,especially in shallow sea environments.However,the OBN receivers may move several times because they are easily affected by tides,currents,and other factors in the shallow sea environment during long-term acquisition.If uncorrected,then the imaging quality of subsequent processing will be affected.The conventional secondary positioning does not consider the case of multiple movements of the receivers,and the accuracy of secondary positioning is insufficient.The first arrival wave of OBN seismic data in shallow ocean mainly comprises refracted waves.In this study,a nonlinear model is established in accordance with the propagation mechanism of a refracted wave and its relationship with the time interval curve to realize the accurate location of multiple receiver movements.In addition,the Levenberg-Marquart algorithm is used to reduce the influence of the first arrival pickup error and to automatically detect the receiver movements,identifying the accurate dynamic relocation of the receivers.The simulation and field data show that the proposed method can realize the dynamic location of multiple receiver movements,thereby improving the accuracy of seismic imaging and achieving high practical value.

Comparative study of nudged elastic band and molecular dynamics methods for diffusion kinetics in solid-state electrolytes

Considerable efforts are being made to transition current lithium-ion and sodium-ion batteries towards the use of solid-state electrolytes.Computational methods,specifically nudged elastic band(NEB)and molecular dynamics(MD)methods,provide powerful tools for the design of solid-state electrolytes.The MD method is usually the choice for studying the materials involving complex multiple diffusion paths or having disordered structures.However,it relies on simulations at temperatures much higher than working temperature.This paper studies the reliability of the MD method using the system of Na diffusion in MgO as a benchmark.We carefully study the convergence behavior of the MD method and demonstrate that total effective simulation time of 12 ns can converge the calculated diffusion barrier to about 0.01 eV.The calculated diffusion barrier is 0.31 eV from both methods.The diffusion coefficients at room temperature are 4.3×10^(-9) cm^(2)⋅s^(−1) and 2.2×10^(-9) cm^(2)⋅s^(−1),respectively,from the NEB and MD methods.Our results justify the reliability of the MD method,even though high temperature simulations have to be employed to overcome the limitation on simulation time.

Dynamic Characteristics of Functionally Graded Timoshenko Beams by Improved Differential Quadrature Method

This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.

Stability and accuracy of central difference method for real-time dynamic substructure testing considering mass participation coefficient

For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.

A flexible multiscale algorithm based on an improved smoothed particle hydrodynamics method for complex viscoelastic flows

Viscoelastic flows play an important role in numerous engineering fields,and the multiscale algorithms for simulating viscoelastic flows have received significant attention in order to deepen our understanding of the nonlinear dynamic behaviors of viscoelastic fluids.However,traditional grid-based multiscale methods are confined to simple viscoelastic flows with short relaxation time,and there is a lack of uniform multiscale scheme available for coupling different solvers in the simulations of viscoelastic fluids.In this paper,a universal multiscale method coupling an improved smoothed particle hydrodynamics(SPH)and multiscale universal interface(MUI)library is presented for viscoelastic flows.The proposed multiscale method builds on an improved SPH method and leverages the MUI library to facilitate the exchange of information among different solvers in the overlapping domain.We test the capability and flexibility of the presented multiscale method to deal with complex viscoelastic flows by solving different multiscale problems of viscoelastic flows.In the first example,the simulation of a viscoelastic Poiseuille flow is carried out by two coupled improved SPH methods with different spatial resolutions.The effects of exchanging different physical quantities on the numerical results in both the upper and lower domains are also investigated as well as the absolute errors in the overlapping domain.In the second example,the complex Wannier flow with different Weissenberg numbers is further simulated by two improved SPH methods and coupling the improved SPH method and the dissipative particle dynamics(DPD)method.The numerical results show that the physical quantities for viscoelastic flows obtained by the presented multiscale method are in consistence with those obtained by a single solver in the overlapping domain.Moreover,transferring different physical quantities has an important effect on the numerical results.

Analysis of mesoscopic mechanical dynamic characteristics of ballast bed with under sleeper pads

The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scanning techniques,and the 3D dynamic model of the rail-sleeper-ballast bed was constructed using the coupled discrete element method-multiflexible-body dynamics(DEM-MFBD)approach.We analyse the meso-mechanical dynamics of the ballast bed with USPs under dynamic load on a train and verify the correctness of the model in laboratory tests.It is shown that the deformation of the USPs increases the contact area between the sleeper and the ballast particles,and subsequently the number of contacts between them.As the depth of the granular ballast bed increases,the contact area becomes larger,and the contact force between the ballast particles gradually decreases.Under the action of the elastic USPs,the contact forces between ballast particles are reduced and the overall vibration level of the ballast bed can be reduced.The settlement of the granular ballast bed occurs mainly at the shallow position of the sleeper bottom,and the installation of the elastic USPs can be effective in reducing the stress on the ballast particles and the settlement of the ballast bed.

Dynamic Modeling and Experimental Verification of an RPR Type Compliant Paralle Mechanism with Low Orders

Efficiency of calculating a dynamic response is an important point of the compliant mechanism for posture adjustment.Dynamic modeling with low orders of a 2R1T compliant parallel mechanism is studied in the paper.The mechanism with two out-of-plane rotational and one lifting degrees of freedom(DoFs)plays an important role in posture adjustment.Based on elastic beam theory,the stiffness matrix and mass matrix of the beam element are established where the moment of inertia is considered.To improve solving efficiency,a dynamic model with low orders of the mechanism is established based on a modified modal synthesis method.Firstly,each branch of the RPR type mechanism is divided into a substructure.Subsequently,a set of hypothetical modes of each substructure is obtained based on the C-B method.Finally,dynamic equation of the whole mechanism is established by the substructure assembly.A dynamic experiment is conducted to verify the dynamic characteristics of the compliant mechanism.

Dynamic Analysis of A Deepwater Drilling Riser with A New Hang-off System

The safety of risers in hang-off states is a vital challenge in offshore oil and gas engineering.A new hang-off system installed on top of risers is proposed for improving the security of risers.This approach leads to a challenging problem:coupling the dynamics of risers with a new hang-off system combined with multiple structures and complex constraints.To accurately analyze the dynamic responses of the coupled system,a coupled dynamic model is established based on the Euler-Bernoulli beam-column theory and penalty function method.A comprehensive analysis method is proposed for coupled dynamic analysis by combining the finite element method and the Newmarkβmethod.An analysis program is also developed in MATLAB for dynamic simulation.The simulation results show that the dynamic performances of the risers at the top part are significantly improved by the new hang-off system,especially the novel design,which includes the centralizer and articulation joint.The bending moment and lateral deformation of the risers at the top part decrease,while the hang-off joint experiences a great bending moment at the bottom of the lateral restraint area which requires particular attention in design and application.The platform navigation speed range under the safety limits of risers expands with the new hang-off system in use.

Structure and dynamical properties during solidification of liquid aluminum induced by cooling and compression

The structural transformation from a liquid into a crystalline solid is an important subject in condensed matter physics and materials science. In the present study, first-principles molecular dynamics calculations are performed to investigate the structure and properties of aluminum during the solidification which is induced by cooling and compression. In the cooling process and compression process, it is found that the icosahedral short-range order is initially enhanced and then begin to decay, the face-centered cubic short-range order eventually becomes dominant before it transforms into a crystalline solid.

Quasi-three-dimensional hydrodynamics of the corona region of laser irradiation of a slab

This paper introduces and establishes a quasi-three-dimensional physical model of the interaction between a laser and a slab target.In contrast to previous one-dimensional analytical models,this paper innovatively fits the real laser conditions based on an isothermal,homogeneous expansion similarity solution of the ideal hydrodynamic equations.Using this simple model,the evolution law and analytical formulae for key parameters(e.g.,temperature,density and scale length)in the corona region under certain conditions are given.The analytical solutions agree well with the relevant results of computational hydrodynamics simulation.For constant laser irradiation,the analytical solutions provide a meaningful power-law scaling relationship.The model provides a set of mathematical and physical tools that give theoretical support for adjusting parameters in experiments.

Experimental and numerical study on dynamic mechanical behaviors of shale under true triaxial compression at high strain rate

High-energy gas fracturing of shale is a novel,high efficacy and eco-friendly mining technique,which is a typical dynamic perturbing behavior.To effectively extract shale gas,it is important to understand the dynamic mechanical properties of shale.Dynamic experiments on shale subjected to true triaxial compression at different strain rates are first conducted in this research.The dynamic stress-strain curves,peak strain,peak stress and failure modes of shale are investigated.The results of the study indicate that the intermediate principal stress and the minor principal stress have the significant influence on the dynamic mechanical behaviors,although this effect decreases as the strain rate increases.The characteristics of compression-shear failure primarily occur in shale subjected to triaxial compression at high strain rates,which distinguishes it from the fragmentation characteristics observed in shale under dynamic uniaxial compression.Additionally,a numerical three-dimensional Split Hopkinson Pressure Bar(3D-SHPB),which is established by coupling PFC3D and FLAC3D methods,is validated to replicate the laboratory characteristics of shale.The dynamic mechanical characteristics of shale subjected to different confining stresses are systematically investigated by the coupling PFC3D and FLAC3D method.The numerical results are in good agreement with the experimental data.

An Arbitrarily High Order and Asymptotic Preserving Kinetic Scheme in Compressible Fluid Dynamic

We present a class of arbitrarily high order fully explicit kinetic numerical methods in compressible fluid dynamics,both in time and space,which include the relaxation schemes by Jin and Xin.These methods can use the CFL number larger or equal to unity on regular Cartesian meshes for the multi-dimensional case.These kinetic models depend on a small parameter that can be seen as a"Knudsen"number.The method is asymptotic preserving in this Knudsen number.Also,the computational costs of the method are of the same order of a fully explicit scheme.This work is the extension of Abgrall et al.(2022)[3]to multidimensional systems.We have assessed our method on several problems for two-dimensional scalar problems and Euler equations and the scheme has proven to be robust and to achieve the theoretically predicted high order of accuracy on smooth solutions.

Dynamic response of buildings under debris flow impact

This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology considers the variables of debris flow depth and velocity and introduces the intensity index IDV(IDV = DV) to evaluate three different levels of debris flow impact intensity. The primary focus of this study is to investigate the dynamic response and failure mechanism of RC-frame buildings under debris flow impact, including structural failure patterns, impact force and column displacement. The results show that under a highintensity impact, a gradual collapse process of the RCframe building can be observed, and the damage mode of the frame column reflects shear failure or plastic hinge failure mechanism. First, the longitudinal infill walls are damaged owing to their low out-of-plane flexural capacity;the critical failure intensity index IDV value is approximately 7.5 m2/s. The structure cannot withstand debris flows with an intensity index IDV greater than 16 m2/s, and it is recommended that the peak impact force should not exceed 2100 k N. The impact damage ability of debris flow on buildings mostly originates from the impact force of the frontal debris flow, with the impact force of the debris flow body being approximately 42% lower than that of the debris flow head. Finally, a five-level classification system for evaluating the damage status of buildings is proposed based on the numerical simulation and investigation results of the disaster site.

Probing the electric double layer structure at nitrogen-doped graphite electrodes by constant-potential molecular dynamics simulations

Electric double layer(EDL)is a critical topic in electrochemistry and largely determines the working performance of lithium batteries.However,atomic insights into the EDL structures on heteroatom-modified graphite anodes and EDL evolution with electrode potential are very lacking.Herein,a constant-potential molecular dynamics(CPMD)method is proposed to probe the EDL structure under working conditions,taking N-doped graphite electrodes and carbonate electrolytes as an example.An interface model was developed,incorporating the electrode potential and atom electronegativities.As a result,an insightful atomic scenario for the EDL structure under varied electrode potentials has been established,which unveils the important role of doping sites in regulating both the EDL structures and the following electrochemical reactions at the atomic level.Specifically,the negatively charged N atoms repel the anions and adsorb Li~+at high and low potentials,respectively.Such preferential adsorption suggests that Ndoped graphite can promote Li~+desolvation and regulate the location of Li~+deposition.This CPMD method not only unveils the mysterious function of N-doping from the viewpoint of EDL at the atomic level but also applies to probe the interfacial structure on other complicated electrodes.

A Dynamical System-Based Framework for Dimension Reduction

We propose a novel framework for learning a low-dimensional representation of data based on nonlinear dynamical systems,which we call the dynamical dimension reduction(DDR).In the DDR model,each point is evolved via a nonlinear flow towards a lower-dimensional subspace;the projection onto the subspace gives the low-dimensional embedding.Training the model involves identifying the nonlinear flow and the subspace.Following the equation discovery method,we represent the vector field that defines the flow using a linear combination of dictionary elements,where each element is a pre-specified linear/nonlinear candidate function.A regularization term for the average total kinetic energy is also introduced and motivated by the optimal transport theory.We prove that the resulting optimization problem is well-posed and establish several properties of the DDR method.We also show how the DDR method can be trained using a gradient-based optimization method,where the gradients are computed using the adjoint method from the optimal control theory.The DDR method is implemented and compared on synthetic and example data sets to other dimension reduction methods,including the PCA,t-SNE,and Umap.

Hierarchical Controller Synthesis Under Linear Temporal Logic Specifications Using Dynamic Quantization

Linear temporal logic(LTL)is an intuitive and expressive language to specify complex control tasks,and how to design an efficient control strategy for LTL specification is still a challenge.In this paper,we implement the dynamic quantization technique to propose a novel hierarchical control strategy for nonlinear control systems under LTL specifications.Based on the regions of interest involved in the LTL formula,an accepting path is derived first to provide a high-level solution for the controller synthesis problem.Second,we develop a dynamic quantization based approach to verify the realization of the accepting path.The realization verification results in the necessity of the controller design and a sequence of quantization regions for the controller design.Third,the techniques of dynamic quantization and abstraction-based control are combined together to establish the local-to-global control strategy.Both abstraction construction and controller design are local and dynamic,thereby resulting in the potential reduction of the computational complexity.Since each quantization region can be considered locally and individually,the proposed hierarchical mechanism is more efficient and can solve much larger problems than many existing methods.Finally,the proposed control strategy is illustrated via two examples from the path planning and tracking problems of mobile robots.

Dynamic Response of A Group of Cylindrical Storage Tanks with Baffles Considering the Effect of Soil Foundation

The sloshing in a group of rigid cylindrical tanks with baffles and on soil foundation under horizontal excitation is studied analytically.The solutions for the velocity potential are derived out by the liquid subdomain method.Equivalent models with mass-spring oscillators are established to replace continuous fluid.Combined with the least square technique,Chebyshev polynomials are employed to fit horizontal,rocking and horizontal-rocking coupling impedances of soil,respectively.A lumped parameter model for impedance is presented to describe the effects of soil on tank structures.A mechanical model for the soil-foundation-tank-liquid-baffle system with small amount of calculation and high accuracy is proposed using the substructure technique.The analytical solutions are in comparison with data from reported literature and numerical codes to validate the effectiveness and correctness of the model.Detailed dynamic properties and seismic responses of the soil-tank system are given for the baffle number,size and location as well as soil parameter.

Numerical Study on the Aerodynamic and Fluid−Structure Interaction of An NREL-5MW Wind Turbine

A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid dynamics(CFD) approach, in which the unsteady, noncompressible Reynolds Averaged Navier-Stokes(RANS) method is used. The main focus of the study is to analyze the tower shadow effect on the aerodynamic performance of the wind turbine under different inlet flow conditions. Subsequently, the finite element model is established by considering fluid/structure interactions to study the structural stress, displacement, strain distributions and flow field information of the structure under the uniform wind speed. Finally, the fluid-structure interaction model is established by considering turbulent wind and the tower shadow effect. The variation rules of the dynamic response of the one-way and two-way fluid-structure interaction(FSI) models under different wind speeds are analyzed, and the numerical calculation results are compared with those of the centralized mass model. The results show that the tower shadow effect and structural deformation are the main factors affecting the aerodynamic load fluctuation of the wind turbine, which in turn affects the aerodynamic performance and structural stability of the blades. The structural dynamic response of the coupled model shows significant similarity, while the structural displacement response of the former exhibits less fluctuation compared with the conventional centralized mass model. The one-way fluid-structure interaction(FSI)model shows a higher frequency of stress-strain and displacement oscillations on the blade compared with the two-way FSI model.