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.展开更多
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 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 order to forecast projectile impact points quickly and accurately,aprojectile impact point prediction method based on generalized regression neural network（GRNN）is presented.Firstly,the model of GRNN forecasting ...In order to forecast projectile impact points quickly and accurately,aprojectile impact point prediction method based on generalized regression neural network（GRNN）is presented.Firstly,the model of GRNN forecasting impact point is established;secondly,the particle swarm algorithm（PSD）is used to optimize the smooth factor in the prediction model and then the optimal GRNN impact point prediction model is obtained.Finally,the numerical simulation of this prediction model is carried out.Simulation results show that the maximum range error is no more than 40 m,and the lateral deviation error is less than0.2m.The average time of impact point prediction is 6.645 ms,which is 1 300.623 ms less than that of numerical integration method.Therefore,it is feasible and effective for the proposed method to forecast projectile impact points,and thus it can provide a theoretical reference for practical engineering applications.展开更多
The flash points of organic compounds were estimated using a hybrid method that includes a simple group contribution method (GCM) implemented in an artificial neural network (ANN) with particle swarm optimization (PSO...The flash points of organic compounds were estimated using a hybrid method that includes a simple group contribution method (GCM) implemented in an artificial neural network (ANN) with particle swarm optimization (PSO). Different topologies of a multilayer neural network were studied and the optimum architecture was determined. Property data of 350 compounds were used for training the network. To discriminate different substances the molecular structures defined by the concept of the classical group contribution method were given as input variables. The capabilities of the network were tested with 155 substances not considered in the training step. The study shows that the proposed GCM+ANN+PSO method represent an excellent alternative for the estimation of flash points of organic compounds with acceptable accuracy (AARD = 1.8%; AAE = 6.2 K).展开更多
An interpolating reproducing kernel particle method for two-dimensional （2D） scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in th...An interpolating reproducing kernel particle method for two-dimensional （2D） scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.展开更多
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.展开更多
In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the met...In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.展开更多
In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of th...In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of the mappings which solves some variational inequality. The results presented extend the corresponding results of Shimizu and Takahashi IT. Shimizu, W. Takahashi, Strong convergence to common fixed point of families of nonexpansive mappings, J. Math. Anal. Appl. 211 （1997）, 71-83], and Yao and Chen [Y. Yao, R. Chert, Convergence to common fixed points of average mappings without commutativity assumption in Hilbert spaces, Nonlinear Analysis 67（2007）, 1758-1763].展开更多
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.展开更多
A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(F...A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.展开更多
The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their ...The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report.展开更多
In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method （CHIIP） for convex nonlinear programming problems. For any convex n...In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method （CHIIP） for convex nonlinear programming problems. For any convex nonlinear programming, without strict convexity for the logarithmic barrier function, we get different solutions of the convex programming in different cases by CHIIP method.展开更多
Objective:To objectively evaluate the clinical effect of resisting method of acupuncture at movementpainful points on treating the refractory cases of frozen shoulder during adhesion period,and to compare with that of...Objective:To objectively evaluate the clinical effect of resisting method of acupuncture at movementpainful points on treating the refractory cases of frozen shoulder during adhesion period,and to compare with that of the conventional acupuncture method.Methods:According to the random grouping principle,90 patients with frozen shoulder during adhesion period were divided into"resisting method acupuncture group,general acupuncture group,and physiotherapy group",30 patients in each.In the resisting method acupuncture group,the resisting method of needling was applied at the"movement-painful points",and during acupuncture the patient was asked to do shoulder movements.In the general acupuncture group,the conventional needling method of acupuncture was applied at the conventional three shoulder-acupoints.In the physiotherapy group,the patients were treated with shock wave at the affected shoulder.All groups were treated once a day,5 times a week,a total of 2 weeks as a course of treatment.After one course of treatment,the clinical effect of the three groups was evaluated.The shoulder pain and shoulder function scale(Constant-Murley),Rating Scale of the American Shoulder and Elbow Surgeons(ASES),and infrared thermal imaging temperature change score were observed.Results:After treatment,the shoulder pain was relieved and movement improved in all three groups,whilst the cured and remarkably effective rate of the resisting method acupuncture group was 96.7%(29/30),better than those of the general acupuncture group 86.2%(25/29)and the physiotherapy group 73.3%(22/30),with a statistically significant difference(P<0.05);the VAS score of the resisting method acupuncture group was significantly lower than those of the other two groups(P<0.05),the ConstantMurley score and the ASES score were significantly lower than those of the other two groups(all P<0.05),and the infrared thermal imaging temperature change score was significantly improved in comparing with those of the other two groups,with a statistically significant difference(both P<0.05).Conclusion:The resisting method of needling has an obvious analgesic effect and is able to restore the joint movement function in the treatment of refractory frozen shoulder during adhesion period.展开更多
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.展开更多
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.展开更多
Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of...Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.展开更多
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.展开更多
基金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.
基金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.
基金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.
基金Project Funded by Chongqing Changjiang Electrical Appliances Industries Group Co.,Ltd
文摘In order to forecast projectile impact points quickly and accurately,aprojectile impact point prediction method based on generalized regression neural network（GRNN）is presented.Firstly,the model of GRNN forecasting impact point is established;secondly,the particle swarm algorithm（PSD）is used to optimize the smooth factor in the prediction model and then the optimal GRNN impact point prediction model is obtained.Finally,the numerical simulation of this prediction model is carried out.Simulation results show that the maximum range error is no more than 40 m,and the lateral deviation error is less than0.2m.The average time of impact point prediction is 6.645 ms,which is 1 300.623 ms less than that of numerical integration method.Therefore,it is feasible and effective for the proposed method to forecast projectile impact points,and thus it can provide a theoretical reference for practical engineering applications.
文摘The flash points of organic compounds were estimated using a hybrid method that includes a simple group contribution method (GCM) implemented in an artificial neural network (ANN) with particle swarm optimization (PSO). Different topologies of a multilayer neural network were studied and the optimum architecture was determined. Property data of 350 compounds were used for training the network. To discriminate different substances the molecular structures defined by the concept of the classical group contribution method were given as input variables. The capabilities of the network were tested with 155 substances not considered in the training step. The study shows that the proposed GCM+ANN+PSO method represent an excellent alternative for the estimation of flash points of organic compounds with acceptable accuracy (AARD = 1.8%; AAE = 6.2 K).
基金supported by the National Natural Science Foundation of China(Grant No.11171208)the Natural Science Foundation of Shanxi Province,China(Grant No.2013011022-6)
文摘An interpolating reproducing kernel particle method for two-dimensional （2D） scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.
基金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.
基金supported by the National Natural Science Foundation of China (11132004 and 51078145)the Natural Science Foundation of Guangdong Province (9251064101000016)
文摘In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.
基金the Thailand Research Fund for financial support under Grant BRG5280016
文摘In this article, we introduce a new viscosity iterative method for two nonexpansive mappings in Hilbert spaces. We also prove, without commutativity assumption, that the iterates converge to a common fixed point of the mappings which solves some variational inequality. The results presented extend the corresponding results of Shimizu and Takahashi IT. Shimizu, W. Takahashi, Strong convergence to common fixed point of families of nonexpansive mappings, J. Math. Anal. Appl. 211 （1997）, 71-83], and Yao and Chen [Y. Yao, R. Chert, Convergence to common fixed points of average mappings without commutativity assumption in Hilbert spaces, Nonlinear Analysis 67（2007）, 1758-1763].
基金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.
基金co-supported by the National Key R&D Program of China(Nos.2018YFF01012401-05)the National Natural Science Foundation of China(No.51975243)+2 种基金Jilin Provincial Department of Education(No.JJKH20180084KJ),Chinathe Fundamental Research Funds for the Central Universities and Jilin Provincial Department of Science&Technology Fund Project,China(Nos.20170101043JC and 20180520072JH)Graduate Innovation Fund of Jilin University,China(No.101832018C184).
文摘A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.
文摘The exact calculation of point kinetic parameters is very important in nuclear reactor safety assessment, and most sophisticated safety codes such as RELAP5, PARCS,DYN3D, and PARET are using these parameters in their dynamic models. These parameters include effective delayed neutron fractions as well as mean generation time.These parameters are adjoint-weighted, and adjoint flux is employed as a weighting function in their evaluation.Adjoint flux calculation is an easy task for most of deterministic codes, but its evaluation is cumbersome for Monte Carlo codes. However, in recent years, some sophisticated techniques have been proposed for Monte Carlo-based point kinetic parameters calculation without any need of adjoint flux. The most straightforward scheme is known as the ‘‘prompt method'' and has been used widely in literature. The main objective of this article is dedicated to point kinetic parameters calculation in Tehran research reactor(TRR) using deterministic as well as probabilistic techniques. WIMS-D5B and CITATION codes have been used in deterministic calculation of forward and adjoint fluxes in the TRR core. On the other hand, the MCNP Monte Carlo code has been employed in the ‘‘prompt method''scheme for effective delayed neutron fraction evaluation.Deterministic results have been cross-checked with probabilistic ones and validated with SAR and experimental data. In comparison with experimental results, the relativedifferences of deterministic as well as probabilistic methods are 7.6 and 3.2%, respectively. These quantities are10.7 and 6.4%, respectively, in comparison with SAR report.
文摘In this paper, on the basis of the logarithmic barrier function and KKT conditions, we propose a combined homotopy infeasible interior-point method （CHIIP） for convex nonlinear programming problems. For any convex nonlinear programming, without strict convexity for the logarithmic barrier function, we get different solutions of the convex programming in different cases by CHIIP method.
基金Doctoral Scientific Research Start-up Funding Project of Fujian Health College:2016-1-1。
文摘Objective:To objectively evaluate the clinical effect of resisting method of acupuncture at movementpainful points on treating the refractory cases of frozen shoulder during adhesion period,and to compare with that of the conventional acupuncture method.Methods:According to the random grouping principle,90 patients with frozen shoulder during adhesion period were divided into"resisting method acupuncture group,general acupuncture group,and physiotherapy group",30 patients in each.In the resisting method acupuncture group,the resisting method of needling was applied at the"movement-painful points",and during acupuncture the patient was asked to do shoulder movements.In the general acupuncture group,the conventional needling method of acupuncture was applied at the conventional three shoulder-acupoints.In the physiotherapy group,the patients were treated with shock wave at the affected shoulder.All groups were treated once a day,5 times a week,a total of 2 weeks as a course of treatment.After one course of treatment,the clinical effect of the three groups was evaluated.The shoulder pain and shoulder function scale(Constant-Murley),Rating Scale of the American Shoulder and Elbow Surgeons(ASES),and infrared thermal imaging temperature change score were observed.Results:After treatment,the shoulder pain was relieved and movement improved in all three groups,whilst the cured and remarkably effective rate of the resisting method acupuncture group was 96.7%(29/30),better than those of the general acupuncture group 86.2%(25/29)and the physiotherapy group 73.3%(22/30),with a statistically significant difference(P<0.05);the VAS score of the resisting method acupuncture group was significantly lower than those of the other two groups(P<0.05),the ConstantMurley score and the ASES score were significantly lower than those of the other two groups(all P<0.05),and the infrared thermal imaging temperature change score was significantly improved in comparing with those of the other two groups,with a statistically significant difference(both P<0.05).Conclusion:The resisting method of needling has an obvious analgesic effect and is able to restore the joint movement function in the treatment of refractory frozen shoulder during adhesion period.
基金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.
基金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.
文摘Seawater intrusion caused by groundwater over-exploitation from coastal aquifers poses a severe problem in many regions. Formulation of proper pumping strategy using a simulation model can assure sustainable supply of fresh water from the coastal aquifers. The focus of the present study is on the development of a numerical model based on Meshfree (MFree) method to study the seawater intrusion problem. For the simulation of seawater intrusion problem, widely used models are based on Finite Difference (FDM) and Finite Element (FEM) Methods, which demand well defined grids/meshes and considerable pre-processing efforts. Here, MFree Point Collocation Method (PCM) based on the Radial Basis Function (RBF) is proposed for the simulation. Diffusive interface approach with density-dependent dispersion and solution of flow and solute transport is adopted. These equations are solved using PCM with appropriate boundary conditions. The developed model has been verified with Henry’s problem, and found to be satisfactory. Further the model has been applied to another established problem and an attempt is made to examine the influence of important system parameters including pumping and recharge on the seawater intrusion. The PCM based MFree model is found computationally efficient as preprocessing is avoided when compared to other numerical methods.
基金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.
