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.展开更多
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 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.展开更多
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 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 intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity mo...The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity model used for the conventional intersection method is based on homogeneous or laterally homogeneous media, which is too simple. In order to improve the accuracy, we have modified the existing intersection method. In the modified approach, the earthquake loci are not assumed to be circular or hyperbolic and calculation accuracy is improved using a minimum traveltime tree algorithm for tracing rays. The numerical model shows that the modified method can locate earthquakes in complex velocity models.展开更多
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.展开更多
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.展开更多
Volume parameter is the basic content of a spatial body object morphology analysis.However,the challenge lies in the volume calculation of irregular objects.The point cloud slicing method proposed in this study effect...Volume parameter is the basic content of a spatial body object morphology analysis.However,the challenge lies in the volume calculation of irregular objects.The point cloud slicing method proposed in this study effectively works in calculating the volume of the point cloud of the spatial object obtained through three-dimensional laser scanning(3DLS).In this method,a uniformly spaced sequent slicing process is first conducted in a specific direction on the point cloud of the spatial object obtained through 3DLS.A series of discrete point cloud slices corresponding to the point cloud bodies are then obtained.Subsequently,the outline boundary polygon of the point cloud slicing is searched one by one in accordance with the slicing sequence and areas of the polygon.The point cloud slice is also calculated.Finally,the individual point cloud section volume is calculated through the slicing areas and the adjacent slicing gap.Thus,the total volume of the scanned spatial object can be calculated by summing up the individual volumes.According to the results and analysis of the calculated examples,the slice-based volume-calculating method for the point cloud of irregular objects obtained through 3DLS is correct,concise in process,reliable in results,efficient in calculation methods,and controllable on accuracy.This method comes as a good solution to the volume calculation of irregular objects.展开更多
As a Lagrangian meshless method, the material point method （MPM） is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method （FEM...As a Lagrangian meshless method, the material point method （MPM） is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method （FEM） for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with that of both FEM and MPM展开更多
A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a...A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a background mesh. It is shown that, in a point collocation approach, the remesh problem because of the mesh distortion in FEM （finite element method） and the low efficiency in Galerkin-based meshfree method are avoided. The corrected kernel functions are introduced to the stabilization of free-surface boundary conditions. The solution of symmetric ring compression problem is compared with a conventional finite element solution, and reasonable results have been obtained.展开更多
Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and the...Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.展开更多
Many approaches inquiring into variational inequality problems have been put forward,among which subgradient extragradient method is of great significance.A novel algorithm is presented in this article for resolving q...Many approaches inquiring into variational inequality problems have been put forward,among which subgradient extragradient method is of great significance.A novel algorithm is presented in this article for resolving quasi-nonexpansive fixed point problem and pseudomonotone variational inequality problem in a real Hilbert interspace.In order to decrease the execution time and quicken the velocity of convergence,the proposed algorithm adopts an inertial technology.Moreover,the algorithm is by virtue of a non-monotonic step size rule to acquire strong convergence theorem without estimating the value of Lipschitz constant.Finally,numerical results on some problems authenticate that the algorithm has preferable efficiency than other algorithms.展开更多
The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and ...The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and operations research. Barker and Pang[1] have given an excellent survey of theories, methods and applications of VIPs.展开更多
Inspired by inertial methods and extragradient algorithms,two algorithms were proposed to investigate fixed point problem of quasinonexpansive mapping and pseudomonotone equilibrium problem in this study.In order to e...Inspired by inertial methods and extragradient algorithms,two algorithms were proposed to investigate fixed point problem of quasinonexpansive mapping and pseudomonotone equilibrium problem in this study.In order to enhance the speed of the convergence and reduce computational cost,the algorithms used a new step size and a cutting hyperplane.The first algorithm was proved to be weak convergence,while the second algorithm used a modified version of Halpern iteration to obtain strong convergence.Finally,numerical experiments on several specific problems and comparisons with other algorithms verified the superiority of the proposed algorithms.展开更多
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.展开更多
基金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 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 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 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.
基金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.
基金This work is supported by the National Natural Science Foundation of China(40674044)the Special Foundation for Basic Professional Scientific Research (DQJB06A02)
文摘The intersection method is one of the basic approaches for locating earthquakes and is not only robust but also efficient. However, its location accuracy is not high, especially for focal depth because the velocity model used for the conventional intersection method is based on homogeneous or laterally homogeneous media, which is too simple. In order to improve the accuracy, we have modified the existing intersection method. In the modified approach, the earthquake loci are not assumed to be circular or hyperbolic and calculation accuracy is improved using a minimum traveltime tree algorithm for tracing rays. The numerical model shows that the modified method can locate earthquakes in complex velocity models.
文摘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.
文摘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.
文摘Volume parameter is the basic content of a spatial body object morphology analysis.However,the challenge lies in the volume calculation of irregular objects.The point cloud slicing method proposed in this study effectively works in calculating the volume of the point cloud of the spatial object obtained through three-dimensional laser scanning(3DLS).In this method,a uniformly spaced sequent slicing process is first conducted in a specific direction on the point cloud of the spatial object obtained through 3DLS.A series of discrete point cloud slices corresponding to the point cloud bodies are then obtained.Subsequently,the outline boundary polygon of the point cloud slicing is searched one by one in accordance with the slicing sequence and areas of the polygon.The point cloud slice is also calculated.Finally,the individual point cloud section volume is calculated through the slicing areas and the adjacent slicing gap.Thus,the total volume of the scanned spatial object can be calculated by summing up the individual volumes.According to the results and analysis of the calculated examples,the slice-based volume-calculating method for the point cloud of irregular objects obtained through 3DLS is correct,concise in process,reliable in results,efficient in calculation methods,and controllable on accuracy.This method comes as a good solution to the volume calculation of irregular objects.
基金supported by the National Basic Research Program of China (2010CB832701)
文摘As a Lagrangian meshless method, the material point method （MPM） is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method （FEM） for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with that of both FEM and MPM
基金the National Natural Science Foundation of China （No. 50275059）.
文摘A meshfree method based on reproducing kernel approximation and point collocation is presented for analysis of metal ring compression. The point collocation method is a true meshfree method without the employment of a background mesh. It is shown that, in a point collocation approach, the remesh problem because of the mesh distortion in FEM （finite element method） and the low efficiency in Galerkin-based meshfree method are avoided. The corrected kernel functions are introduced to the stabilization of free-surface boundary conditions. The solution of symmetric ring compression problem is compared with a conventional finite element solution, and reasonable results have been obtained.
基金Project(2010CB732103)supported by the National Basic Research Program of ChinaProject(51179092)supported by the National Natural Science Foundation of ChinaProject(2012-KY-02)supported by the State Key Laboratory of Hydroscience and Engineering,China
文摘Recently,the radial point interpolation meshfree method has gained popularity owing to its advantages in large deformation and discontinuity problems,however,the accuracy of this method depends on many factors and their influences are not fully investigated yet.In this work,three main factors,i.e.,the shape parameters,the influence domain size,and the nodal distribution,on the accuracy of the radial point interpolation method(RPIM)are systematically studied and conclusive results are obtained.First,the effect of shape parameters(R,q)of the multi-quadric basis function on the accuracy of RPIM is examined via global search.A new interpolation error index,closely related to the accuracy of RPIM,is proposed.The distribution of various error indexes on the R q plane shows that shape parameters q[1.2,1.8]and R[0,1.5]can give good results for general 3-D analysis.This recommended range of shape parameters is examined by multiple benchmark examples in 3D solid mechanics.Second,through numerical experiments,an average of 30 40 nodes in the influence domain of a Gauss point is recommended for 3-D solid mechanics.Third,it is observed that the distribution of nodes has significant effect on the accuracy of RPIM although it has little effect on the accuracy of interpolation.Nodal distributions with better uniformity give better results.Furthermore,how the influence domain size and nodal distribution affect the selection of shape parameters and how the nodal distribution affects the choice of influence domain size are also discussed.
文摘Many approaches inquiring into variational inequality problems have been put forward,among which subgradient extragradient method is of great significance.A novel algorithm is presented in this article for resolving quasi-nonexpansive fixed point problem and pseudomonotone variational inequality problem in a real Hilbert interspace.In order to decrease the execution time and quicken the velocity of convergence,the proposed algorithm adopts an inertial technology.Moreover,the algorithm is by virtue of a non-monotonic step size rule to acquire strong convergence theorem without estimating the value of Lipschitz constant.Finally,numerical results on some problems authenticate that the algorithm has preferable efficiency than other algorithms.
基金The NNSF (10071031) of China and National 973 Project.
文摘The finite-dimensional variational inequality problem (VIP) has been studied extensively in the literature because of its successful applications in many fields such as economics, transportation, regional science and operations research. Barker and Pang[1] have given an excellent survey of theories, methods and applications of VIPs.
文摘Inspired by inertial methods and extragradient algorithms,two algorithms were proposed to investigate fixed point problem of quasinonexpansive mapping and pseudomonotone equilibrium problem in this study.In order to enhance the speed of the convergence and reduce computational cost,the algorithms used a new step size and a cutting hyperplane.The first algorithm was proved to be weak convergence,while the second algorithm used a modified version of Halpern iteration to obtain strong convergence.Finally,numerical experiments on several specific problems and comparisons with other algorithms verified the superiority of the proposed algorithms.
文摘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.
