Material point method simulation of hydro-mechanical behaviour in twophase porous geomaterials: A state-of-the-art review

The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.

A hybrid contact approach for modeling soil-structure interaction using the material point method

The grid-based multi-velocity field technique has become increasingly popular for simulating the Material Point Method(MPM)in contact problems.However,this traditional technique has some shortcomings,such as(1)early contact and contact penetration can occur when the contact conditions are unsuitable,and(2)the method is not available for contact problems involving rigid-nonrigid materials,which can cause numerical instability.This study presents a new hybrid contact approach for the MPM to address these limitations to simulate the soil and structure interactions.The approach combines the advantages of point-point and point-segment contacts to implement contact detection,satisfying the impenetrability condition and smoothing the corner contact problem.The proposed approach is first validated through a disk test on an inclined slope.Then,several typical cases,such as granular collapse,bearing capacity,and deformation of a flexible retaining wall,are simulated to demonstrate the robustness of the proposed approach compared with FEM or analytical solutions.Finally,the proposed method is used to simulate the impact of sand flow on a deformable structure.The results show that the proposed contact approach can well describe the phenomenon of soil-structure interaction problems.

Mobility and dynamic erosion process of granular flow:insights from numerical investigation using material point method

In order to understand the dynamics of granular flow on an erodible base soil,in this paper,a series of material point method-based granular column collapse tests were conducted to investigate numerically the mobility and dynamic erosion process of granular flow subjected to the complex settings,i.e.,the aspect ratio,granular mass,friction and dilatancy resistance,gravity and presence of water.A set of power scaling laws were proposed to describe the final deposit characteristics of granular flow by the relations of the normalized run-out distance and the normalized final height of granular flow against the aspect ratio,being greatly affected by the complex geological settings,e.g.,granular mass,the friction and dilatancy resistance of granular soil,and presence of water in granular flow.An index of the coefficient of friction of granular soil was defined as a ratio of the target coefficient of friction over the initial coefficient of friction to quantify the scaling extent of friction change(i.e.,friction strengthening or weakening).There is a characteristic aspect ratio of granular column corresponding to the maximum mobility of granular flow with the minimum index of the apparent coefficient of friction.The index of the repose coefficient of friction of granular flow decreased gradually with the increase in aspect ratio because higher potential energy of granular column at a larger aspect ratio causes a larger kinetic energy of granular soil to weaken the friction of granular soil as a kind of velocity-related friction weakening.An increase in granular mass reduces gradually the indexes of the apparent and repose coefficients of friction of granular soil to enhance the mobility of granular flow.The mobility of granular flow increases gradually with the decrease in friction angle or increase in dilatancy angle of granular soil.However,the increase of gravity accelerates granular flow but showing the same final deposit profile without any dependence on gravity.The mobility of granular flow increases gradually by lowering the indexes of the apparent and repose coefficients of friction of granular flow while changing the surroundings,in turn,the dry soil,submerged soil and saturated soil,implying a gradually increased excessive mobility of granular flow with the friction weakening of granular soil.Presence of water in granular flow may be a potential catalyzer to yield a long run-out granular flow,as revealed in comparison of water-absent and water-present granular flows.In addition,the dynamic erosion and entrainment of based soil induced by granular flow subjected to the complex geological settings,i.e.,the aspect ratio,granular mass,gravity,friction and dilatancy resistance,and presence of water,were comprehensively investigated as well.

Transfer matrix method for free and forced vibrations of multi-level functionally graded material stepped beams with different boundary conditions

Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.

Multi-layer perceptron-based data-driven multiscale modelling of granular materials with a novel Frobenius norm-based internal variable

One objective of developing machine learning(ML)-based material models is to integrate them with well-established numerical methods to solve boundary value problems(BVPs).In the family of ML models,recurrent neural networks(RNNs)have been extensively applied to capture history-dependent constitutive responses of granular materials,but these multiple-step-based neural networks are neither sufficiently efficient nor aligned with the standard finite element method(FEM).Single-step-based neural networks like the multi-layer perceptron(MLP)are an alternative to bypass the above issues but have to introduce some internal variables to encode complex loading histories.In this work,one novel Frobenius norm-based internal variable,together with the Fourier layer and residual architectureenhanced MLP model,is crafted to replicate the history-dependent constitutive features of representative volume element(RVE)for granular materials.The obtained ML models are then seamlessly embedded into the FEM to solve the BVP of a biaxial compression case and a rigid strip footing case.The obtained solutions are comparable to results from the FEM-DEM multiscale modelling but achieve significantly improved efficiency.The results demonstrate the applicability of the proposed internal variable in enabling MLP to capture highly nonlinear constitutive responses of granular materials.

Existence and numerical approximation of a solution to frictional contact problem for electro-elastic materials

In this paper,a frictional contact problem between an electro-elastic body and an electrically conductive foundation is studied.The contact is modeled by normal compliance with finite penetration and a version of Coulomb’s law of dry friction in which the coefficient of friction depends on the slip.In addition,the effects of the electrical conductivity of the foundation are taken into account.This model leads to a coupled system of the quasi-variational inequality of the elliptic type for the displacement and the nonlinear variational equation for the electric potential.The existence of a weak solution is proved by using an abstract result for elliptic variational inequalities and a fixed point argument.Then,a finite element approximation of the problem is presented.Under some regularity conditions,an optimal order error estimate of the approximate solution is derived.Finally,a successive iteration technique is used to solve the problem numerically and a convergence result is established.

Friction Characteristics Between Marine Clay and Construction Materials

Structure-soil interface friction characteristics is of importance to investigate the interaction between engineering structures and soils,especially for offshore structures.The interface friction behavior between marine clay and structural materials with different roughness was studied in this paper by using 3D optical scanning tests,a modified direct shear device and numerical simulation.Relationships between the surface roughness of structures,water content and interface friction angle were presented by model tests.The increase of water contents decreased the interface friction angles.For interfaces with different roughness,the interface friction angles will be smaller than that of the soil when the water content exceeds a certain value.The roughness of the interface and the water content of the soil are mutually coupled to influence the coefficient of friction(COF).This paper proposed a Finite Element Method(FEM)to simulate the interface direct shear tests of structures with different roughness.The surface models with different roughness are established based on the structure data obtained by 3D scanning.The Coupled Eulerian-Lagrangian(CEL)approach was employed to analyse soils sheared by irregular surfaces.The interface behavior for interfaces with different roughness under cyclic shear stresses was analyzed by FEM.

Sol-gel synthesis of nanometer silicon/silicon suboxide/carbon anode material

A stacked Si/SiO_(x)/C composite anode material with carbon-coated structure was prepared by sol-gel method combined with carbothermal reduction using organic silicon.The results of X-ray diffractometry, scanning electron microscopy, and elemental analysis show that the Si/SiO_(x)/C material is a secondary particle with a porous micronanostructure, and the presence of nanometer silicon does not affect the carbothermal reduction and carbon coating.Electrochemical test results indicate that the specific capacity and first coulombic efficiency of SiO_(x)/C composite with nanometer silicon can be increased to 1 946.05 mAh/g and 76.49%,respectively.The reversible specific capacity of Si/SiO_(x)/C material blended with graphite is 749.69 mAh/g after 100 cycles at a current density of 0.1 C,and the capacity retention rate is up to 89.03%.Therefore, the composite has excellent electrochemical cycle stability.

Design,preparation,application of advanced array structured materials and their action mechanism analyses for high performance lithium-sulfur batteries

Lithium-sulfur battery(LSB)has brought much attention and concern because of high theoretical specific capacity and energy density as one of main competitors for next-generation energy storage systems.The widely commercial application and development of LSB is mainly hindered by serious“shuttle effect”of lithium polysulfides(Li PSs),slow reaction kinetics,notorious lithium dendrites,etc.In various structures of LSB materials,array structured materials,possessing the composition of ordered micro units with the same or similar characteristics of each unit,present excellent application potential for various secondary cells due to some merits such as immobilization of active substances,high specific surface area,appropriate pore sizes,easy modification of functional material surface,accommodated huge volume change,enough facilitated transportation for electrons/lithium ions,and special functional groups strongly adsorbing Li PSs.Thus many novel array structured materials are applied to battery for tackling thorny problems mentioned above.In this review,recent progresses and developments on array structured materials applied in LSBs including preparation ways,collaborative structural designs based on array structures,and action mechanism analyses in improving electrochemical performance and safety are summarized.Meanwhile,we also have detailed discussion for array structured materials in LSBs and constructed the structure-function relationships between array structured materials and battery performances.Lastly,some directions and prospects about preparation ways,functional modifications,and practical applications of array structured materials in LSBs are generalized.We hope the review can attract more researchers'attention and bring more studying on array structured materials for other secondary batteries including LSB.

Modeling footing resting on anisotropic sand using material point method

Sand typically exhibits anisotropic internal structure which may significantly influence its mechanical behavior. The material point method (MPM) can eliminate mesh distortion and thus is suitable for investigating geotechnical problems with large deformation. In this study, an advanced anisotropic critical state theory (ACST)-based soil model is implemented in MPM to study the response of strip footing resting on anisotropic sand. The capability of the model is verified by simulating several element tests and strip footing tests with different soil densities and fabric bedding plane orientations. For the footing problem with a vertical load, as the fabric bedding plane orientation increases, the bearing capacity decreases and its corresponding settlement increases. The failure pattern becomes asymmetrical when the bedding plane orientation or the loading direction is inclined. A comparison between the simulation results predicted by the anisotropic and isotropic models is made, which demonstrates that neglecting the fabric anisotropy may lead to the overestimation of the bearing capacity.

Modeling pipe-soil interaction under vertical downward relative offset using B-spline material point method

To analyze the pipeline response under permanent ground deformation,the evolution of resistance acting on the pipe during the vertical downward offset is an essential ingredient.However,the efficient simulation of pipe penetration into soil is challenging for the conventional finite element(FE)method due to the large deformation of the surrounding soils.In this study,the B-spline material point method(MPM)is employed to investigate the pipe-soil interaction during the downward movement of rigid pipes buried in medium and dense sand.To describe the density-and stress-dependent behaviors of sand,the J2-deformation type model with state-dependent dilatancy is adopted.The effectiveness of the model is demonstrated by element tests and biaxial compression tests.Afterwards,the pipe penetration process is simulated,and the numerical outcomes are compared with the physical model tests.The effects of pipe size and burial depth are investigated with an emphasis on the mobilization of the soil resistance and the failure mechanisms.The simulation results indicate that the bearing capacity formulas given in the guidelines can provide essentially reasonable estimates for the ultimate force acting on buried pipes,and the recommended value of yield displacement may be underestimated to a certain extent.

Constitutive modelling of idealised granular materials using machine learning method

Predicting the constitutive response of granular soils is a fundamental goal in geomechanics.This paper presents a machine learning(ML)framework for the prediction of the stress-strain behaviour and shearinduced contact fabric evolution of an idealised granular material subject to triaxial shearing.The MLbased framework is comprised of a set of mini-triaxial tests which provide a benchmark for the setup and validation of the discrete element method(DEM)model of the granular materials,a parametric DEM simulation programme of virtual triaxial tests which provides datasets of micro-and macro-mechanical information,as well as a multi-layer perceptron(MLP)neural network which is trained and tested using the DEM-based datasets.The ML model only requires the initial void ratio of the granular sample as the input for predicting its constitutive response.The excellent agreement between the ML model prediction and experimental test and DEM simulation results indicates that the MLebased modelling approach is capable of capturing accurately the effects of initial void ratio on the constitutive behaviour of idealised granular materials,bypassing the need to incorporate the complex micromechanics underlying the macroscopic mechanical behaviour of granular materials.Lastly,a detailed comparison between the used MLP model and long short-term memory(LSTM)model was made from the perspective of technical algorithm,prediction accuracy,and computational efficiency.

Research Progress of High Entropy Ceramic Materials

High-entropy materials(HEMs)have better mechanical,thermal,and electrical properties than traditional materials due to their special"high entropy effect".They can also adjust the performance of high entropy ceramics by adjusting the proportion of raw materials,and have broad application prospects in many fields.This article provides a review of the high entropy effect,preparation methods,and main applications of high entropy ceramic materials,especially exploring relevant research on high entropy perovskite ceramics.It is expected to provide reference for the promotion of scientific research and the development of further large-scale applications of high-entropy ceramic materials.

Topology Optimization of Sound-Absorbing Materials for Two-Dimensional Acoustic Problems Using Isogeometric Boundary Element Method

In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.

Displacement Field Variable Modeling Method for Heterogeneous Materials in Wind Power Blade Core Plates

In order to study the mechanical properties of the heterogeneous core plate of the wind turbine blade,a modeling method of the core plate based on displacement field variables is proposed.Firstly,the wind turbine blade core plate was modeled according to the theory of modeling heterogeneous material characteristics.Secondly,the three-point bending finite element model of the wind turbine blade core plate was solved by the display dynamic equation to obtain the deformation pattern and force-deformation relationship of the core plate.Finally,the three-point bending static test was conducted to compare with the finite element analysis.The test results show that:the damage form of the wind turbine blade core plate includes elasticity,yield,and failure stages.The main failure modes are plastic deformation,core material collapse,and panel-core delamination.The failure load measured by the test is 1.59 kN,which is basically consistent with the load-displacement result obtained by the simulation,with a difference of only 1.9%,which verifies the validity and reliability of the model.It provides data references for wind turbine blade structure design.

Coseismic site response and slope instability using periodic boundary conditions in the material point method

This paper proposed the explicit generalized-a time scheme and periodic boundary conditions in the material point method(MPM)for the simulation of coseismic site response.The proposed boundary condition uses an intuitive particle-relocation algorithm ensuring material points always remain within the computational mesh.The explicit generalized-a time scheme was implemented in MPM to enable the damping of spurious high frequency oscillations.Firstly,the MPM was verified against finite element method(FEM).Secondly,ability of the MPM in capturing the analytical transfer function was investigated.Thirdly,a symmetric embankment was adopted to investigate the effects of ground motion arias intensity(I_(a)),geometry dimensions,and constitutive models.The results show that the larger the model size,the higher the crest runout and settlement for the same ground motion.When using a Mohr-Coulomb model,the crest runout increases with increasing I_(a).However,if the strain-softening law is activated,the results are less influenced by the ground motion.Finally,the MPM results were compared with the Newmark sliding block solution.The simplified analysis herein highlights the capabilities of MPM to capture the full deformation process for earthquake engineering applications,the importance of geometry characterization,and the selection of appropriate constitutive models when simulating coseismic site response and subsequent large deformations.

Advances in Mn‑Based Electrode Materials for Aqueous Sodium‑Ion Batteries

Aqueous sodium-ion batteries have attracted extensive attention for large-scale energy storage applications,due to abundant sodium resources,low cost,intrinsic safety of aqueous electrolytes and eco-friendliness.The electrochemical performance of aqueous sodium-ion batteries is affected by the properties of electrode materials and electrolytes.Among various electrode materials,Mn-based electrode materials have attracted tremendous attention because of the abundance of Mn,low cost,nontoxicity,eco-friendliness and interesting electrochemical performance.Aqueous electrolytes having narrow electrochemical window also affect the electrochemical performance of Mn-based electrode materials.In this review,we introduce systematically Mn-based electrode materials for aqueous sodium-ion batteries from cathode and anode materials and offer a comprehensive overview about their recent development.These Mn-based materials include oxides,Prussian blue analogues and polyanion compounds.We summarize and discuss the composition,crystal structure,morphology and electrochemical properties of Mn-based electrode materials.The improvement methods based on electrolyte optimization,element doping or substitution,optimization of morphology and carbon modification are highlighted.The perspectives of Mn-based electrode materials for future studies are also provided.We believe this review is important and helpful to explore and apply Mn-based electrode materials in aqueous sodium-ion batteries.

Accuracy and Precision of Energy Dispersive X-Ray Fluorescence (EDXRF) Analysis of Trace and Major Elements in Rock Standard Reference Materials Using Fine Powder

In this work the performance of a screening analytical method for Energy Dispersive X-Ray Fluorescence (EDXRF) analysis in terms of accuracy and precision was evaluated through analysis of rock standard reference materials. The method allowed the division of elements into four groups taking into account the excitation energies and measurement conditions of the sample. Two standard reference materials were used and 15 sample replicates were prepared and analyzed, then statistics were applied to assess the precision and accuracy of analytical results. The obtained results show that major compounds or elements (SiO2, P2O5, K2O, CaO, Fe2O3, Ti) can be determined in fine powder sample with a deviation lower than 15%, and a relative standard deviation in the range (1 - 10)%. The deviation was found to be lower than 5% for major compounds such as K2O, and CaO, which suggest that the EDXRF is accurate in evaluating major elemental concentrations in rock samples. It was also found that the method seems to be more accurate and precise for major elements than for trace element investigation. This screening analytical method can be used for routine analysis with acceptable results, even though the method should be optimized to increase its precision and accuracy.

Inverse design for material anisotropy and its application for a compact X-cut TFLN on-chip wavelength demultiplexer

Inverse design focuses on identifying photonic structures to optimize the performance of photonic devices.Conventional scalar-based inverse design approaches are insufficient to design photonic devices of anisotropic materials such as lithium niobate(LN).To the best of our knowledge,this work proposes for the first time the inverse design method for anisotropic materials to optimize the structure of anisotropic-material based photonics devices.Specifically,the orientation dependent properties of anisotropic materials are included in the adjoint method,which provides a more precise prediction of light propagation within such materials.The proposed method is used to design ultra-compact wavelength division demultiplexers in the X-cut thin-film lithium niobate(TFLN)platform.By benchmarking the device performances of our method with those of classical scalar-based inverse design,we demonstrate that this method properly addresses the critical issue of material anisotropy in the X-cut TFLN platform.This proposed method fills the gap of inverse design of anisotropic materials based photonic devices,which finds prominent applications in TFLN platforms and other anisotropicmaterial based photonic integration platforms.

A two-scale method for identifying mechanical parameters of composite materials with periodic configuration

In this paper, a two-scale method （TSM） is presented for identifying the mechanics parameters such as stiffness and strength of composite materials with small periodic configuration. Firstly, a formulation is briefly given for two-scale analysis （TSA） of the composite materials. And then a two-scale computation formulation of strains and stresses is developed by displacement solution with orthotropic material coefficients for three kinds of such composites structures, i.e., the tension column with a square cross section, the bending cantilever with a rectangular cross section and the torsion column with a circle cross section. The strength formulas for the three kinds of structures are derived and the TSM procedure is discussed. Finally the numerical results of stiffness and strength are presented and compared with experimental data. It shows that the TSM method in this paper is feasible and valid for predicting both the stiffness and the strength of the composite materials with periodic configuration.