Geo-interfaces refer to the contact surfaces between multiple media within geological strata,as well as the transition zones that regulate the migration of three-phase matter,changes in physical states,and the deforma...Geo-interfaces refer to the contact surfaces between multiple media within geological strata,as well as the transition zones that regulate the migration of three-phase matter,changes in physical states,and the deformation and stability of rock and soil masses.Owing to the combined effects of natural factors and human activities,geo-interfaces play crucial roles in the emergence,propagation,and triggering of geological disasters.Over the past three decades,the material point method(MPM)has emerged as a preferred approach for addressing large deformation problems and simulating soil-water-structure interactions,making it an ideal tool for analyzing geo-interface behaviors.In this review,we offer a systematic summary of the basic concepts,classifications,and main characteristics of the geo-interface,and provide a comprehensive overview of recent advances and developments in simulating geo-interface using the MPM.We further present a brief description of various MPMs for modeling different types of geo-interfaces in geotechnical engineering applications and highlight the existing limitations and future research directions.This study aims to facilitate innovative applications of the MPM in modeling complex geo-interface problems,providing a reference for geotechnical practitioners and researchers.展开更多
This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-...This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-oped using the material point method.To reduce the computational cost of Monte Carlo simulations,response surface models are created as surrogate models for the material point system to approximate its dynamic behavior.An adaptive randomized greedy algorithm is employed to construct a sparse polynomial chaos expansion model with a fixed order,effectively balancing the accuracy and computational efficiency of the surrogate model.Based on the sparse polynomial chaos expansion,sensitivity analysis is conducted using the global finite difference and Sobol methods.Several examples of structural dynamics are provided to demonstrate the effectiveness of the proposed method in addressing structural dynamics problems.展开更多
In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this sche...In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this scheme seamlessly combines the material point method(MPM)for accurately tracking solid-phase deformation and heat transport,and the Eulerian finite element method(FEM)for effectively capturing fluid flow and heat advection-diffusion behavior.The proposed approach circumvents the substantial challenges posed by large nonlinear equation systems with the monolithic solution scheme.The staggered solution process strategically separates each physical field through explicit or implicit integration.The characteristic-based method is used to stabilize advection-dominated heat flows for efficient numerical implementation.Furthermore,a fractional step approach is employed to decompose fluid velocity and pressure,thereby suppressing pore pressure oscillation on the linear background grid.The fracturing initiation and propagation are simulated by a rate-dependent phase field model.Through a series of quasi-static and transient simulations,the exceptional performance and promising potential of the proposed model in addressing THM fracturing problems in poro-elastic media is demonstrated.展开更多
This study proposed a deep learning-based nanoindentation simulation method to address the challenge of obtaining the mechanical parameters of rock-forming minerals and the complexity of regression analysis.This appro...This study proposed a deep learning-based nanoindentation simulation method to address the challenge of obtaining the mechanical parameters of rock-forming minerals and the complexity of regression analysis.This approach enables the accurate assessment of rock-forming minerals'mechanical parameters.A material database of nanoindentation load-depth(P-h)curves was generated using the material point method(MPM)to characterize the mechanical behavior of major rock-forming minerals(quartz,albite,and muscovite)in sandstone.We used Bayesian hyperparameter optimization to determine the optimal hyperparameters for training a deep neural network(DNN).The trained DNN model accurately predicted the material parameters of rock-forming minerals using experimental nanoindentation P-h data.Numerical simulations of the uniaxial compression of heterogeneous sandstones were conducted using the predicted parameters to assess the sandstones’macro-mechanical characteristics.The research findings provide new insights into the fundamental mechanical behavior of heterogeneous rock materials.展开更多
The Three Gorges Reservoir area hosts numerous high-risk columnar rock masses,but predicting their full-process of collapse dynamics in narrow-deep terrains remains challenging.This study investigated the Changdongzi ...The Three Gorges Reservoir area hosts numerous high-risk columnar rock masses,but predicting their full-process of collapse dynamics in narrow-deep terrains remains challenging.This study investigated the Changdongzi columnar dangerous rock mass(CCDRM)through field surveys and simulated its collapse dynamics with Material Point Method(MPM),and analyzed the collapse-diffusionsurge generation mechanisms.Field investigations revealed a 0.8 m-thick lotus-like extrusion structure at the base,severely degraded by cyclic water fluctuations,predisposing the rock mass to base-crushing failure.Numerical simulations indicate collapse velocities up to 42 m/s,amplified spreading with a 318 m deposition range and 11 m maximum thickness,and two major collisions causing 50%-60%velocity reduction.Gully topography buffers kinetic energy,retaining 5%debris on slopes and channeling 7%into the river.High-speed river impact generates a 20 m-diameter cavity,producing an 11 m surge,identifying the Guandukou Town waterfront 1.2 km downstream as a high-risk zone.These results clarify the collapse mechanisms of columnar dangerous rock mass and inform its hazard mitigation strategies.展开更多
Large deformation analysis of slope failure is important for hazard and risk assessment of infrastructure.Recent studies have revealed that spatial variability of soil properties can significantly affect the probabili...Large deformation analysis of slope failure is important for hazard and risk assessment of infrastructure.Recent studies have revealed that spatial variability of soil properties can significantly affect the probability of slope failure.However,due to limitations of traditional numerical tools,the influence of spatial variability of soil properties on the post-failure behavior of slopes has not been fully understood.Therefore,in this study,we aimed to investigate the effects of the cross-correlation between cohesion and the friction angle on the probability of slope failure and post-failure behavior(e.g.run-out distance,influence distance,and influence zone)using a random material point method(RMPM).The study showed that mesh size,strength reduction shape factor parameter,and residual strength all play critical roles in the calculated post-failure behavior of a slope.Based on stochastic Monte Carlo simulation,the effects of cross-correlation between cohesion and the friction angle on the probability of slope failure,and its run-out distance,influence distance,influence zone,and sliding volume were studied.The study also showed that material point method(MPM)has great advantages compared with the finite element method(FEM)in handling large deformations.展开更多
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 stat...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.展开更多
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 pip...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.展开更多
Three-dimensional rock fracture induced by blasting is a highly complex problem and has received considerable attention in geotechnical engineering.The material point method is firstly applied to treat this challengin...Three-dimensional rock fracture induced by blasting is a highly complex problem and has received considerable attention in geotechnical engineering.The material point method is firstly applied to treat this challenging task.Some inherent weaknesses can be overcome by coupling the generalized interpolation material point(GIMP)and the convected particle domain interpolation technique(CPDI).For the media in the borehole,unchanged GIMP-type particles are used to guarantee a homogenous blast pressure.CPDITetrahedron type particles are employed to avoid the fake numerical fracture near the borehole for the rock material.A blasting experiment using three-dimensional single-borehole rock was simulated to examine the applicability of the coupled model under realistic loading and boundary conditions.A good agreement was achieved between the simulation and experimental results.Moreover,the mechanism of three-dimensional rock fracture was analyzed.It was concluded that rock particle size and material parameters play an important role in rock damage.The reflected tensile waves cause severe damage in the lower part of the model.Rayleigh waves occur on the top face of the rock model to induce a hoop failure band.展开更多
We applied the material balance principle of the denudation volume and sedimentary flux to study the denudation-accumulation system between the Longmen Mountains(Mts.)and the foreland basin.The amount of sediment in e...We applied the material balance principle of the denudation volume and sedimentary flux to study the denudation-accumulation system between the Longmen Mountains(Mts.)and the foreland basin.The amount of sediment in each sedimentation stage of the basin was estimated to obtain the denudation volume,erosion thickness and deposit thickness since the Late Triassic Epoch,to enable us to recover the paleoelevation of the provenance and the sedimentary area.The results show the following:(1)Since the Late Triassic Epoch,the elevation of the surface of the Longmen Mts.has uplifted from 0 m to 2751 m,and the crust of the Longmen Mts.has uplifted by 9.8 km.Approximately 72%of the materials introduced have been denuded from the mountains.(2)It is difficult to recover the paleoelevation of each stage of the Longmen Mts.foreland basin quantitatively by the present-day techniques and data.(3)The formation of the Longmen Mts.foreland basin consisted of three stages of thrust belt tectonic load and three stages of thrust belt erosional unload.During tectonic loading stages(Late Triassic Epoch,Late Jurassic-Early Cretaceous,Late Cretaceous-Miocene),the average elevation of Longmen Mts.was lower(approximately 700-1700 m).During erosional unloading stages(Early and Middle Jurassic,Middle Cretaceous and Jiaguan,Late Cenozoic),the average elevation of Longmen Mts.was high at approximately 2000-2800m.展开更多
This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit m...This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit material point method(iMPM).The global stiffness matrix is formed explicitly and the Newton-Raphson iterative method is used to solve the equilibrium equations.Where possible,the implementation procedure closely follows standard finite element method(FEM)approaches to allow easy conversion of other FEM codes.The generalized interpolation function is assigned to eliminate the inherent cell crossing noise within conventional MPM.For the first time,the B-bar approach is used to overcome volumetric locking in standard GIMP method for near-incompressible non-linear geomechanics.The proposed i GIMP was tested and compared with i MPM and analytical solutions via a 1 D column compression problem.Results highlighted the superiority of the i GIMP approach in reducing stress oscillations,thereby improving computational accuracy.Then,elasto-plastic slope stabilities and rigid footing problems were considered,further illustrating the ability of the proposed method to overcome volumetric locking due to incompressibility.Results showed that the proposed i GIMP with B-bar approach can be used to simulate geotechnical problems with large deformations.展开更多
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...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.展开更多
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 ge...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.展开更多
Advances in numerical simulation techniques play an important role in helpingmining engineers understand those parts of the rock mass that cannot be readily observed.The Material Point Method(MPM)is an example of such...Advances in numerical simulation techniques play an important role in helpingmining engineers understand those parts of the rock mass that cannot be readily observed.The Material Point Method(MPM)is an example of such a tool that is gaining popularity for studying geotechnical problems.In recent years,the original formulation of MPM has been extended to not only account for simulating the mechanical behaviour of rock under different loading conditions,but also to describe the coupled interaction of pore water and solid phases in materials.These methods assume that the permeability of mediums is homogeneous,and we show that these MPM techniques fail to accurately capture the correct behaviour of the fluid phase if the permeability of the material is heterogeneous.In this work,we propose a novel implementation of the coupled MPM to address this problem.We employ an approach commonly used in coupled Finite Volume Methods,known as the Two Point Flux Approximation(TPFA).Our new method is benchmarked against two well-known analytical expressions(a one-dimensional geostatic consolidation and the so-called Mandel-Cryer effect).Its performance is compared to existing coupled MPM approaches for homogeneous materials.In order to gauge the possible effectiveness of our technique in the field,we apply ourmethod to a case study relating to a mine known to experience severe problems with pore water.展开更多
ZnO nanostructure materials doped with different La contents were synthesized by sonochemical method. The products were characterized by X-ray diffraction（XRD）, scanning electron microscopy（SEM）,transmission elect...ZnO nanostructure materials doped with different La contents were synthesized by sonochemical method. The products were characterized by X-ray diffraction（XRD）, scanning electron microscopy（SEM）,transmission electron microscopy（TEM）, Raman spectroscopy, and Fourier transform infrared spectroscopy（FTIR）. In this research, XRD patterns of pure ZnO and La-doped ZnO are specified as hexagonal wurtzite ZnO structure with no detection of La2O3 phase. SEM and TEM characterization revealed the flower shape of pure ZnO built-up from petals of hexagonal prisms with hexagonal pyramid tips. Upon doping with La, the flower-shaped ZnO is broken into individual 1D prism-like nanorods. Photocatalytic activities of the as-synthesized products were determined by measuring the degradation of methylene blue（MB） under ultraviolet–visible（UV） light irradiation.Among them, the 2.0 mol% La-doped ZnO shows better photocatalytic properties than any other products.展开更多
Representative volume element （RVE） method and asymptotic homogenization （AH） method are two widely used methods in predicting effective properties of pe- riodic materials. This paper develops a novel implementa- ...Representative volume element （RVE） method and asymptotic homogenization （AH） method are two widely used methods in predicting effective properties of pe- riodic materials. This paper develops a novel implementa- tion of the AH method, which has rigorous mathematical foundation of the AH method, and also simplicity as the RVE method. This implementation can be easily realized using commercial software as a black box, and can use all kinds of elements available in commercial software to model unit cells with rather complicated microstructures, so the model may remain a fairly small scale. Several examples were car- fled out to demonstrate the simplicity and effectiveness of the new implementation.展开更多
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 intuit...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.展开更多
Outgoing waves arising from high-velocity impacts between soil and structure can be reflected by the conventional truncated boundaries.Absorbing boundary conditions(ABCs),to attenuate the energy of the outward waves,a...Outgoing waves arising from high-velocity impacts between soil and structure can be reflected by the conventional truncated boundaries.Absorbing boundary conditions(ABCs),to attenuate the energy of the outward waves,are necessary to ensure the proper representation of the kinematic field and the accurate quantification of impact forces.In this paper,damping layer and dashpot ABCs are implemented in the material point method(MPM)with slight adjustments.Benchmark scenarios of different dynamic problems are modelled with the ABCs configured.Feasibility of the ABCs is assessed through the velocity fluctuations at specific observation points and the impact force fluctuations on the structures.The impact forces predicted by the MPM with ABCs are verified by comparison with those estimated using a computational fluid dynamics approach.展开更多
In a bird strike, the bird undergoes large deformation like flows; while most part of the structure is in small deformation, the region near the impact point may experience large deformations, even fail. This paper de...In a bird strike, the bird undergoes large deformation like flows; while most part of the structure is in small deformation, the region near the impact point may experience large deformations, even fail. This paper develops a coupled shell-material point method （CSMPM） for bird strike simulation, in which the bird is modeled by the material point method （MPM） and the aircraft structure is modeled by the Belytschko-Lin-Tsay shell element. The interaction between the bird and the structure is handled by a particle-to-surface contact algorithm. The distorted and failed shell elements will be eroded if a certain criterion is reached. The proposed CSMPM takes full advantages of both the finite element method and the MPM for bird strike simulation and is validated by several numerical examples.展开更多
To better simulate multi-phase interactions involving failure evolution, the material point method （MPM） has evolved for almost twenty years. Recently, a particle-based multiscale simulation procedure is being devel...To better simulate multi-phase interactions involving failure evolution, the material point method （MPM） has evolved for almost twenty years. Recently, a particle-based multiscale simulation procedure is being developed, within the framework of the MPM, to describe the detonation process of energetic nano-composites from molecular to continuum level so that a multiscale equation of state could be formulated. In this letter, a multiscale MPM is proposed via both hierarchical and concurrent schemes to simulate the impact response between two microrods with different nanostructures. Preliminary results are presented to illustrate that a transition region is not required between different spatial scales with the proposed approach.展开更多
基金supported by the National Science Fund for Distinguished Young Scholars of China(Grant No.42225702)the National Natural Science Foundation of China(Grant Nos.42461160266 and 52379106).
文摘Geo-interfaces refer to the contact surfaces between multiple media within geological strata,as well as the transition zones that regulate the migration of three-phase matter,changes in physical states,and the deformation and stability of rock and soil masses.Owing to the combined effects of natural factors and human activities,geo-interfaces play crucial roles in the emergence,propagation,and triggering of geological disasters.Over the past three decades,the material point method(MPM)has emerged as a preferred approach for addressing large deformation problems and simulating soil-water-structure interactions,making it an ideal tool for analyzing geo-interface behaviors.In this review,we offer a systematic summary of the basic concepts,classifications,and main characteristics of the geo-interface,and provide a comprehensive overview of recent advances and developments in simulating geo-interface using the MPM.We further present a brief description of various MPMs for modeling different types of geo-interfaces in geotechnical engineering applications and highlight the existing limitations and future research directions.This study aims to facilitate innovative applications of the MPM in modeling complex geo-interface problems,providing a reference for geotechnical practitioners and researchers.
基金support from the National Natural Science Foundation of China(Grant Nos.52174123&52274222).
文摘This paper presents a framework for constructing surrogate models for sensitivity analysis of structural dynamics behavior.Physical models involving deformation,such as collisions,vibrations,and penetration,are devel-oped using the material point method.To reduce the computational cost of Monte Carlo simulations,response surface models are created as surrogate models for the material point system to approximate its dynamic behavior.An adaptive randomized greedy algorithm is employed to construct a sparse polynomial chaos expansion model with a fixed order,effectively balancing the accuracy and computational efficiency of the surrogate model.Based on the sparse polynomial chaos expansion,sensitivity analysis is conducted using the global finite difference and Sobol methods.Several examples of structural dynamics are provided to demonstrate the effectiveness of the proposed method in addressing structural dynamics problems.
基金supported by National Natural Science Foundation of China(Grant No.42377149)the Research Grants Council of Hong Kong(General Research Fund Project No.17202423).
文摘In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this scheme seamlessly combines the material point method(MPM)for accurately tracking solid-phase deformation and heat transport,and the Eulerian finite element method(FEM)for effectively capturing fluid flow and heat advection-diffusion behavior.The proposed approach circumvents the substantial challenges posed by large nonlinear equation systems with the monolithic solution scheme.The staggered solution process strategically separates each physical field through explicit or implicit integration.The characteristic-based method is used to stabilize advection-dominated heat flows for efficient numerical implementation.Furthermore,a fractional step approach is employed to decompose fluid velocity and pressure,thereby suppressing pore pressure oscillation on the linear background grid.The fracturing initiation and propagation are simulated by a rate-dependent phase field model.Through a series of quasi-static and transient simulations,the exceptional performance and promising potential of the proposed model in addressing THM fracturing problems in poro-elastic media is demonstrated.
基金supported by the National Key Research and Development Program of China(Grant no.2023YFC3009005)the Chongqing Technology Innovation and Application Development Special Key Project(Grant no.CSTB2022TIAD-KPX0135)the Fundamental Research Funds for the Central Universities(Grant no.2023CDJKYJH068).
文摘This study proposed a deep learning-based nanoindentation simulation method to address the challenge of obtaining the mechanical parameters of rock-forming minerals and the complexity of regression analysis.This approach enables the accurate assessment of rock-forming minerals'mechanical parameters.A material database of nanoindentation load-depth(P-h)curves was generated using the material point method(MPM)to characterize the mechanical behavior of major rock-forming minerals(quartz,albite,and muscovite)in sandstone.We used Bayesian hyperparameter optimization to determine the optimal hyperparameters for training a deep neural network(DNN).The trained DNN model accurately predicted the material parameters of rock-forming minerals using experimental nanoindentation P-h data.Numerical simulations of the uniaxial compression of heterogeneous sandstones were conducted using the predicted parameters to assess the sandstones’macro-mechanical characteristics.The research findings provide new insights into the fundamental mechanical behavior of heterogeneous rock materials.
基金supported by the Regional Innovation Development Joint Fund of the National Natural Science Foundation of China(U23A2045)the scientific research project of the Chongqing Municipal Bureau of Planning and Natural Resources(KJ-2023046)the scientific and technological project of the Hubei Provincial Bureau of Geological Survey(KJ2024-34)。
文摘The Three Gorges Reservoir area hosts numerous high-risk columnar rock masses,but predicting their full-process of collapse dynamics in narrow-deep terrains remains challenging.This study investigated the Changdongzi columnar dangerous rock mass(CCDRM)through field surveys and simulated its collapse dynamics with Material Point Method(MPM),and analyzed the collapse-diffusionsurge generation mechanisms.Field investigations revealed a 0.8 m-thick lotus-like extrusion structure at the base,severely degraded by cyclic water fluctuations,predisposing the rock mass to base-crushing failure.Numerical simulations indicate collapse velocities up to 42 m/s,amplified spreading with a 318 m deposition range and 11 m maximum thickness,and two major collisions causing 50%-60%velocity reduction.Gully topography buffers kinetic energy,retaining 5%debris on slopes and channeling 7%into the river.High-speed river impact generates a 20 m-diameter cavity,producing an 11 m surge,identifying the Guandukou Town waterfront 1.2 km downstream as a high-risk zone.These results clarify the collapse mechanisms of columnar dangerous rock mass and inform its hazard mitigation strategies.
基金the Fund of Hong Kong Research GrantsCouncil(RGC)(No.16214519),China。
文摘Large deformation analysis of slope failure is important for hazard and risk assessment of infrastructure.Recent studies have revealed that spatial variability of soil properties can significantly affect the probability of slope failure.However,due to limitations of traditional numerical tools,the influence of spatial variability of soil properties on the post-failure behavior of slopes has not been fully understood.Therefore,in this study,we aimed to investigate the effects of the cross-correlation between cohesion and the friction angle on the probability of slope failure and post-failure behavior(e.g.run-out distance,influence distance,and influence zone)using a random material point method(RMPM).The study showed that mesh size,strength reduction shape factor parameter,and residual strength all play critical roles in the calculated post-failure behavior of a slope.Based on stochastic Monte Carlo simulation,the effects of cross-correlation between cohesion and the friction angle on the probability of slope failure,and its run-out distance,influence distance,influence zone,and sliding volume were studied.The study also showed that material point method(MPM)has great advantages compared with the finite element method(FEM)in handling large deformations.
基金The financial supports from National Outstanding Youth Science Fund Project of National Natural Science Foundation of China(Grant No.52022112)the International Postdoctoral Exchange Fellowship Program(Talent-Introduction Program,Grant No.YJ20220219)。
文摘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.
基金supported by the National Natural Science Foundation of China(Grant Nos.42225702,42077235 and 41722209).
文摘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.
基金This research was funded by the Natural Science Foundation of Sichuan,China(No.2022NSFSC1915)the National Natural Science Foundation of China(No.U19A2098)+1 种基金State Key Laboratory of Precision Blasting and Hubei Key Laboratory of Blasting Engineering,Jianghan University(No.PBSKL2022B06)the Fundamental Research Funds for the Central Universities。
文摘Three-dimensional rock fracture induced by blasting is a highly complex problem and has received considerable attention in geotechnical engineering.The material point method is firstly applied to treat this challenging task.Some inherent weaknesses can be overcome by coupling the generalized interpolation material point(GIMP)and the convected particle domain interpolation technique(CPDI).For the media in the borehole,unchanged GIMP-type particles are used to guarantee a homogenous blast pressure.CPDITetrahedron type particles are employed to avoid the fake numerical fracture near the borehole for the rock material.A blasting experiment using three-dimensional single-borehole rock was simulated to examine the applicability of the coupled model under realistic loading and boundary conditions.A good agreement was achieved between the simulation and experimental results.Moreover,the mechanism of three-dimensional rock fracture was analyzed.It was concluded that rock particle size and material parameters play an important role in rock damage.The reflected tensile waves cause severe damage in the lower part of the model.Rayleigh waves occur on the top face of the rock model to induce a hoop failure band.
基金the Project of the National Natural Science Foudation of China(Grant No.41372114,41340005,41172162,40972083)
文摘We applied the material balance principle of the denudation volume and sedimentary flux to study the denudation-accumulation system between the Longmen Mountains(Mts.)and the foreland basin.The amount of sediment in each sedimentation stage of the basin was estimated to obtain the denudation volume,erosion thickness and deposit thickness since the Late Triassic Epoch,to enable us to recover the paleoelevation of the provenance and the sedimentary area.The results show the following:(1)Since the Late Triassic Epoch,the elevation of the surface of the Longmen Mts.has uplifted from 0 m to 2751 m,and the crust of the Longmen Mts.has uplifted by 9.8 km.Approximately 72%of the materials introduced have been denuded from the mountains.(2)It is difficult to recover the paleoelevation of each stage of the Longmen Mts.foreland basin quantitatively by the present-day techniques and data.(3)The formation of the Longmen Mts.foreland basin consisted of three stages of thrust belt tectonic load and three stages of thrust belt erosional unload.During tectonic loading stages(Late Triassic Epoch,Late Jurassic-Early Cretaceous,Late Cretaceous-Miocene),the average elevation of Longmen Mts.was lower(approximately 700-1700 m).During erosional unloading stages(Early and Middle Jurassic,Middle Cretaceous and Jiaguan,Late Cenozoic),the average elevation of Longmen Mts.was high at approximately 2000-2800m.
基金the National Natural Science Foundation of China(Nos.41807223 and 51908175)the Fundamental Research Funds for the Central Universities(No.B210202096)+1 种基金the Natural Science Foundation of Guangdong Province(No.2018A030310346)the Water Conservancy Science and Technology Innovation Project of Guangdong Province(No.2020-11),China。
文摘This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit material point method(iMPM).The global stiffness matrix is formed explicitly and the Newton-Raphson iterative method is used to solve the equilibrium equations.Where possible,the implementation procedure closely follows standard finite element method(FEM)approaches to allow easy conversion of other FEM codes.The generalized interpolation function is assigned to eliminate the inherent cell crossing noise within conventional MPM.For the first time,the B-bar approach is used to overcome volumetric locking in standard GIMP method for near-incompressible non-linear geomechanics.The proposed i GIMP was tested and compared with i MPM and analytical solutions via a 1 D column compression problem.Results highlighted the superiority of the i GIMP approach in reducing stress oscillations,thereby improving computational accuracy.Then,elasto-plastic slope stabilities and rigid footing problems were considered,further illustrating the ability of the proposed method to overcome volumetric locking due to incompressibility.Results showed that the proposed i GIMP with B-bar approach can be used to simulate geotechnical problems with large deformations.
基金This work was supported by the National Natural Science Foundation of China(Grant no.U22A20603)Sichuan Science and Technology Program-China(Grant No.2023ZYD0149)CAS"Light of West China"Program-China(Grant No.Fangwei Yu).In addition,a special acknowledgement should be expressed to a famous Chinese television drama:My Chief and My Regiment that accompanied me(Dr.Fangwei Yu)through the loneliness time of completing this study.
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.52108359).
文摘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.
文摘Advances in numerical simulation techniques play an important role in helpingmining engineers understand those parts of the rock mass that cannot be readily observed.The Material Point Method(MPM)is an example of such a tool that is gaining popularity for studying geotechnical problems.In recent years,the original formulation of MPM has been extended to not only account for simulating the mechanical behaviour of rock under different loading conditions,but also to describe the coupled interaction of pore water and solid phases in materials.These methods assume that the permeability of mediums is homogeneous,and we show that these MPM techniques fail to accurately capture the correct behaviour of the fluid phase if the permeability of the material is heterogeneous.In this work,we propose a novel implementation of the coupled MPM to address this problem.We employ an approach commonly used in coupled Finite Volume Methods,known as the Two Point Flux Approximation(TPFA).Our new method is benchmarked against two well-known analytical expressions(a one-dimensional geostatic consolidation and the so-called Mandel-Cryer effect).Its performance is compared to existing coupled MPM approaches for homogeneous materials.In order to gauge the possible effectiveness of our technique in the field,we apply ourmethod to a case study relating to a mine known to experience severe problems with pore water.
基金financially supported by the National Research University Project for Chiang Mai University (CMU) from the Thailand’s Office of the Higher Education Commission, Thailand
文摘ZnO nanostructure materials doped with different La contents were synthesized by sonochemical method. The products were characterized by X-ray diffraction（XRD）, scanning electron microscopy（SEM）,transmission electron microscopy（TEM）, Raman spectroscopy, and Fourier transform infrared spectroscopy（FTIR）. In this research, XRD patterns of pure ZnO and La-doped ZnO are specified as hexagonal wurtzite ZnO structure with no detection of La2O3 phase. SEM and TEM characterization revealed the flower shape of pure ZnO built-up from petals of hexagonal prisms with hexagonal pyramid tips. Upon doping with La, the flower-shaped ZnO is broken into individual 1D prism-like nanorods. Photocatalytic activities of the as-synthesized products were determined by measuring the degradation of methylene blue（MB） under ultraviolet–visible（UV） light irradiation.Among them, the 2.0 mol% La-doped ZnO shows better photocatalytic properties than any other products.
基金supported by the National Natural Science Foundation of China(91216201)
文摘Representative volume element （RVE） method and asymptotic homogenization （AH） method are two widely used methods in predicting effective properties of pe- riodic materials. This paper develops a novel implementa- tion of the AH method, which has rigorous mathematical foundation of the AH method, and also simplicity as the RVE method. This implementation can be easily realized using commercial software as a black box, and can use all kinds of elements available in commercial software to model unit cells with rather complicated microstructures, so the model may remain a fairly small scale. Several examples were car- fled out to demonstrate the simplicity and effectiveness of the new implementation.
基金funded by National Science Foundation(NSF)(Grant No.CMMI-2211002).
文摘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.
基金the Key Science and Technology Plan of Power China Huadong Engineering Corporation(No.KY2018-ZD-01)China and the National Natural Science Foundations of China(No.51909248)。
文摘Outgoing waves arising from high-velocity impacts between soil and structure can be reflected by the conventional truncated boundaries.Absorbing boundary conditions(ABCs),to attenuate the energy of the outward waves,are necessary to ensure the proper representation of the kinematic field and the accurate quantification of impact forces.In this paper,damping layer and dashpot ABCs are implemented in the material point method(MPM)with slight adjustments.Benchmark scenarios of different dynamic problems are modelled with the ABCs configured.Feasibility of the ABCs is assessed through the velocity fluctuations at specific observation points and the impact force fluctuations on the structures.The impact forces predicted by the MPM with ABCs are verified by comparison with those estimated using a computational fluid dynamics approach.
基金Supported by the National Natural Science Foundation of China(11390363)
文摘In a bird strike, the bird undergoes large deformation like flows; while most part of the structure is in small deformation, the region near the impact point may experience large deformations, even fail. This paper develops a coupled shell-material point method （CSMPM） for bird strike simulation, in which the bird is modeled by the material point method （MPM） and the aircraft structure is modeled by the Belytschko-Lin-Tsay shell element. The interaction between the bird and the structure is handled by a particle-to-surface contact algorithm. The distorted and failed shell elements will be eroded if a certain criterion is reached. The proposed CSMPM takes full advantages of both the finite element method and the MPM for bird strike simulation and is validated by several numerical examples.
基金supported in part by the U.S.Defense Threat Reduction Agency(HDTRA1-10-1-0022)the National Basic Research Program of China(2010CB832704)+2 种基金the National Natural Science Foundation of China(10721062)the 111 Joint Program by the Chinese Ministry of EducationState Administration of Foreign Experts Affairs(B08014)
文摘To better simulate multi-phase interactions involving failure evolution, the material point method （MPM） has evolved for almost twenty years. Recently, a particle-based multiscale simulation procedure is being developed, within the framework of the MPM, to describe the detonation process of energetic nano-composites from molecular to continuum level so that a multiscale equation of state could be formulated. In this letter, a multiscale MPM is proposed via both hierarchical and concurrent schemes to simulate the impact response between two microrods with different nanostructures. Preliminary results are presented to illustrate that a transition region is not required between different spatial scales with the proposed approach.
