On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is present...On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.展开更多
The complex variable reproducing kernel particle method （CVRKPM） of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape...The complex variable reproducing kernel particle method （CVRKPM） of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method （RKPM） and the element- free Galerkin （EFG） method.展开更多
As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle d...As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle displacements are decoupled in nature,thus making this method suitable for parallelization.The FPM also requires an acceleration strategy to overcome the heavy computational burden of its explicit framework for time-dependent dynamic analysis.To this end,a GPU-accelerated parallel strategy for the FPM is proposed in this paper.By taking advantage of the independence of each step of the FPM workflow,a generic parallelized computational framework for multiple types of analysis is established.Using the Compute Unified Device Architecture(CUDA),the GPU implementations of the main tasks of the FPM,such as evaluating and assembling the element equivalent forces and solving the kinematic equations for particles,are elaborated through careful thread management and memory optimization.Performance tests show that speedup ratios of 8,25 and 48 are achieved for beams,hexahedral solids and triangular shells,respectively.For examples consisting of explicit dynamic analyses of shells and solids,comparisons with Abaqus using 1 to 8 CPU cores validate the accuracy of the results and demonstrate a maximum speed improvement of a factor of 11.2.展开更多
An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes（SNN） and local-Delaunay-tri...An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes（SNN） and local-Delaunay-triangulation（LDT） techniques, which were suitable and effective for h-adaptivity analysis on 2-D problems with the regular or irregular distribution of the nodes. The results of multiresolution and h- adaptivity analyses on 2-D linear elastostatics and bending plate problems demonstrate that the improper high-gradient indicator will reduce the convergence property of the h- adaptivity analysis, and that the efficiency of the LDT node refinement strategy is better than SNN, and that the presented h-adaptivity analysis scheme is provided with the validity, stability and good convergence property.展开更多
A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow t...A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.展开更多
The cohesive solids in liquid flows are featured by the dynamic growth and breakage of agglomerates, and the difficulties in the development, design and optimization of these systems are related to this significant fe...The cohesive solids in liquid flows are featured by the dynamic growth and breakage of agglomerates, and the difficulties in the development, design and optimization of these systems are related to this significant feature.In this paper, discrete particle method is used to simulate a solid–liquid flow system including millions of cohesive particles, the growth rate and breakage rate of agglomerates are then systematically investigated. It was found that the most probable size of the agglomerates is determined by the balance of growth and breakage of the agglomerates the cross point of the lines of growth rate and breakage rate as a function of the particle numbers in an agglomerate, marks the most stable agglomerate size. The finding here provides a feasible way to quantify the dynamic behaviors of growth and breakage of agglomerates, and therefore offers the possibility of quantifying the effects of agglomerates on the hydrodynamics of fluid flows with cohesive particles.展开更多
Large deformation contact problems generally involve highly nonlinear behaviors,which are very time-consuming and may lead to convergence issues.The finite particle method(FPM)effectively separates pure deformation fr...Large deformation contact problems generally involve highly nonlinear behaviors,which are very time-consuming and may lead to convergence issues.The finite particle method(FPM)effectively separates pure deformation from total motion in large deformation problems.In addition,the decoupled procedures of the FPM make it suitable for parallel computing,which may provide an approach to solve time-consuming issues.In this study,a graphics processing unit(GPU)-based parallel algorithm is proposed for two-dimensional large deformation contact problems.The fundamentals of the FPM for planar solids are first briefly introduced,including the equations of motion of particles and the internal forces of quadrilateral elements.Subsequently,a linked-list data structure suitable for parallel processing is built,and parallel global and local search algorithms are presented for contact detection.The contact forces are then derived and directly exerted on particles.The proposed method is implemented with main solution procedures executed in parallel on a GPU.Two verification problems comprising large deformation frictional contacts are presented,and the accuracy of the proposed algorithm is validated.Furthermore,the algorithm’s performance is investigated via a large-scale contact problem,and the maximum speedups of total computational time and contact calculation reach 28.5 and 77.4,respectively,relative to commercial finite element software Abaqus/Explicit running on a single-core central processing unit(CPU).The contact calculation time percentage of the total calculation time is only 18%with the FPM,much smaller than that(50%)with Abaqus/Explicit,demonstrating the efficiency of the proposed method.展开更多
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.展开更多
In this paper, the complex variable reproducing kernel particle （CVRKP） method and the finite element （FE） method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE metho...In this paper, the complex variable reproducing kernel particle （CVRKP） method and the finite element （FE） method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.展开更多
The Reproducing Kernel Particle Method (RKPM) is one of several new meshless numerical methods de- veloped internationally in recent years. The ideal elasto-plastic constitutive model of material under a Taylor impact...The Reproducing Kernel Particle Method (RKPM) is one of several new meshless numerical methods de- veloped internationally in recent years. The ideal elasto-plastic constitutive model of material under a Taylor impact is characterized by the Jaumann stress- and strain-rates. An updated Lagrangian format is used for the calculation in a nu- merical analysis. With the RKPM, this paper deals with the calculation model for the Taylor impact and deduces the control equation for the impact process. A program was developed to simulate numerically the Taylor impact of projec- tiles composed of several kinds of material. The simulation result is in good accordance with both the test results and the Taylor analysis outcome. Since the meshless method is not limited by meshes, it is believed to be widely applicable to such complicated processes as the Taylor impact, including large deformation and strain and to the study of the dy- namic qualities of materials.展开更多
Most natural resources are processed as particle-fluid multiphase systems in chemical,mineral and material indus-tries,therefore,discrete particles methods(DPM)are reasonable choices of simulation method for engineeri...Most natural resources are processed as particle-fluid multiphase systems in chemical,mineral and material indus-tries,therefore,discrete particles methods(DPM)are reasonable choices of simulation method for engineering the relevant processes and equipments.However,direct application of these methods is challenged by the complex multiscale behavior of such systems,which leads to enormous computational cost or otherwise qualitatively inac-curate description of the mesoscale structures.The coarse-grained DPM based on the energy-minimization multi-scale(EMMS)model,or EMMS-DPM,was proposed to reduce the computational cost by several orders while main-taining an accurate description of the mesoscale structures,which paves the way for its engineering applications.Further empowered by the high-efficiency multi-scale DEM software DEMms and the corresponding customized heterogeneous supercomputing facilities with graphics processing units(GPUs),it may even approach realtime simulation of industrial reactors.This short review will introduce the principle of DPM,in particular,EMMS-DPM,and the recent developments in modeling,numerical implementation and application of large-scale DPM which aims to reach industrial scale on one hand and resolves mesoscale structures critical to reaction-transport coupling on the other hand.This review finally prospects on the future developments of DPM in this direction.展开更多
During splitting rolling simulation, re-meshing is necessary to prevent the effect of severe mesh distortion when the conventional finite element method is used. However, extreme deformation cannot be solved by the fi...During splitting rolling simulation, re-meshing is necessary to prevent the effect of severe mesh distortion when the conventional finite element method is used. However, extreme deformation cannot be solved by the finite element method in splitting rolling. The reproducing kernel particle method can solve this problem because the continuum body is discretized by a set of nodes, and a finite element mesh is unnecessary, and there is no explicit limitation of mesh when the metal is split. To ensure stability in the large deformation elastoplastic analysis, the Lagrange material shape function was introduced. The transformation method was utilized to impose the essential boundary conditions. The splitting rolling method was simulated and the simulation results were in accordance with the experimental ones in the literature.展开更多
The moving particle semi-implicit（MPS）method has demonstrated its usefulness in practical engineering applications.Although it has wide applicability,it is still hard to predict the pressure precisely using the MPS ...The moving particle semi-implicit（MPS）method has demonstrated its usefulness in practical engineering applications.Although it has wide applicability,it is still hard to predict the pressure precisely using the MPS method.A pressure-convection particle method based on the MPS method is proposed to overcome this problem.The improved performance of this new method is validated with computational and measured results.The approach is also applied to compute the problem of sloshing associated with floating body motion in waves.The pressure-convection MPS method demonstrated its capability to improve the prediction of pressure.展开更多
In this paper, the normal derivative of the radial basis function （RBF） is introduced into the reproducing kernel particle method （RKPM）, and the improved reproducing kernel particle method （IRKPM） is proposed. ...In this paper, the normal derivative of the radial basis function （RBF） is introduced into the reproducing kernel particle method （RKPM）, and the improved reproducing kernel particle method （IRKPM） is proposed. The method can decrease the errors on the boundary and improve the accuracy and stability of the algorithm. The proposed method is applied to the numerical simulation of piezoelectric materials and the corresponding governing equations are derived. The numerical results show that the IRKPM is more stable and accurate than the RKPM.展开更多
To study the non-linear fracture, a non-linear constitutive model for piezoelectric ceramics was proposed, in which the polarization switching and saturation were taken into account. Based on the model, the non-linear...To study the non-linear fracture, a non-linear constitutive model for piezoelectric ceramics was proposed, in which the polarization switching and saturation were taken into account. Based on the model, the non-linear fracture analysis was implemented using reproducing kernel particle method （RKPM）. Using local J-integral as a fracture criterion, a relation curve of fracture loads against electric fields was obtained. Qualitatively, the curve is in agreement with the experimental observations reported in literature. The reproducing equation, the shape function of RKPM, and the transformation method to impose essential boundary conditions for meshless methods were also introduced. The computation was implemented using object-oriented programming method.展开更多
In recent years,the network continues to enter people’s lives,followed by network security issues that continue to appear,causing substantial economic losses to the world.As an effective method to tackle the network ...In recent years,the network continues to enter people’s lives,followed by network security issues that continue to appear,causing substantial economic losses to the world.As an effective method to tackle the network security issues,intrusion detection system has been widely used and studied.In this paper,the NSL-KDD data set is used to reduce the dimension of data features,remove the features of low correlation and high interference,and improve the computational efficiency.To improve the detection rate and accuracy of intrusion detection,this paper introduces the particle method for the first time that we call it intrusion detection with particle(IDP).To illustrate the effectiveness of this method,experiments are carried out on three kinds of data-before dimension reduction,after dimension reduction and importing particle method based on dimension reduction.By comparing the results of DT,NN,SVM,K-NN,and NB,it is proved that the particle method can effectively improve the intrusion detection rate.展开更多
Stable and controllable solid flow is essential in circulating fluidized bed (CFB) systems. The L-valve is a typical non-mechanical valve that can provide flexible solid feeding. The investigation of the solid circula...Stable and controllable solid flow is essential in circulating fluidized bed (CFB) systems. The L-valve is a typical non-mechanical valve that can provide flexible solid feeding. The investigation of the solid circulation rate and the hydrodynamic characteristics of the L-valve is crucial to its design and operation. The gas-solid flow in the L-valve of a full-loop CFB is studied with the coarse-grained discrete particle method (EMMS-DPM). Good agreements on the solid circulation rate and the pressure drop through the L-valve are achieved between the simulated and experimental data. The solid circulation rate increases linearly with the aeration velocity until the stable particle circulation of the CFB is destroyed. The flow patterns in the horizontal section of L-valve are gas-solid slug flow above the stationary solid layer and the moving solid layer, respectively. The effects of L-valve geometric parameters on the solid flow characteristics are also investigated. The results indicate that reducing the diameter and length of the horizontal section of L-valve can improve the solid transport efficiency, especially at low aeration velocity. Besides, the solid conveying capacity and flow stability are improved when the sharp bend of L-valve is modified to be a gradual bend.展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
This study presents a structural analysis algorithm called the finite particle method （FPM） for kinematically indeterminate bar assemblies. Different from the traditional analysis method, FPM is based on the combina...This study presents a structural analysis algorithm called the finite particle method （FPM） for kinematically indeterminate bar assemblies. Different from the traditional analysis method, FPM is based on the combination of the vector mechanics and numerical calculations. It models the analyzed domain composed of finite particles. Newton＇s second law is adopted to describe the motions of all particles. A convected material flame and explicit time integration for the solution procedure is also adopted in this method. By using the FPM, there is no need to solve any nonlinear equations, to calculate the stiffness matrix or equilibrium matrix, which is very helpful in the analysis of kinematically indeterminate structures. The basic formulations for the space bar are derived, following its solution procedures for bar assemblies. Three numerical examples are analyzed using the FPM. Results obtained from both the straight pretension cable and the suspension cable assembly show that the FPM can produce a more accurate analysis result. The motion simulation of the four-bar space assembly demonstrates the capability of this method in the analysis ofkinematically indeterminate structures.展开更多
A mathematical model has been formulated based on the combined continuous and discrete particle method for investigating the sedimentation behaviour of microparticles in aqueous suspensions, by treating the fluid phas...A mathematical model has been formulated based on the combined continuous and discrete particle method for investigating the sedimentation behaviour of microparticles in aqueous suspensions, by treating the fluid phase as continuous and the particles phase as discrete, thus allowing the behaviour of individual particles to be followed and the evolution of the structure of the particle phase to be investigated as a function of time. The model takes into account most of the prevailing forces acting on individual particles including van der Waals attractive, electrostatic repulsive, gravitational, Brownian, depletion, steric, contact and drag forces. A code has also been developed based on the model. This paper reports some preliminary modelling results of mono-dispersed microparticles settling in aqueous suspensions under various conditions. The results show the short time dynamics of the fluid phase, which has a similar order of magnitude to the particle phase. Such short time dynamics could bear significance to processes such as particle aggregation when their size becomes very small. Preliminary analyses of the results have also been carried out on the evolution of particle settling based on a newly proposed parameter, local normalised volume fraction （LNVF）.展开更多
基金supported by the National Natural Science Foundation of China (Grant No.10871124)the Innovation Program of Shanghai Municipal Education Commission,China (Grant No.09ZZ99)
文摘On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.
基金supported by the National Natural Science Foundation of China (Grant No. 11171208)the Leading Academic Discipline Project of Shanghai City,China (Grant No. S30106)
文摘The complex variable reproducing kernel particle method （CVRKPM） of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method （RKPM） and the element- free Galerkin （EFG） method.
基金the financial support provided by the National Key Research and Development Program of China(Grant No.2016YFC0800200)the National Natural Science Foundation of China(Grant Nos.51578494 and 51778568)the Fundamental Research Funds for the Central Universities(Grant No.2019QNA4043).
文摘As a novel kind of particle method for explicit dynamics,the finite particle method(FPM)does not require the formation or solution of global matrices,and the evaluations of the element equivalent forces and particle displacements are decoupled in nature,thus making this method suitable for parallelization.The FPM also requires an acceleration strategy to overcome the heavy computational burden of its explicit framework for time-dependent dynamic analysis.To this end,a GPU-accelerated parallel strategy for the FPM is proposed in this paper.By taking advantage of the independence of each step of the FPM workflow,a generic parallelized computational framework for multiple types of analysis is established.Using the Compute Unified Device Architecture(CUDA),the GPU implementations of the main tasks of the FPM,such as evaluating and assembling the element equivalent forces and solving the kinematic equations for particles,are elaborated through careful thread management and memory optimization.Performance tests show that speedup ratios of 8,25 and 48 are achieved for beams,hexahedral solids and triangular shells,respectively.For examples consisting of explicit dynamic analyses of shells and solids,comparisons with Abaqus using 1 to 8 CPU cores validate the accuracy of the results and demonstrate a maximum speed improvement of a factor of 11.2.
文摘An h-adaptivity analysis scheme based on multiple scale reproducing kernel particle method was proposed, and two node refinement strategies were constructed using searching-neighbor-nodes（SNN） and local-Delaunay-triangulation（LDT） techniques, which were suitable and effective for h-adaptivity analysis on 2-D problems with the regular or irregular distribution of the nodes. The results of multiresolution and h- adaptivity analyses on 2-D linear elastostatics and bending plate problems demonstrate that the improper high-gradient indicator will reduce the convergence property of the h- adaptivity analysis, and that the efficiency of the LDT node refinement strategy is better than SNN, and that the presented h-adaptivity analysis scheme is provided with the validity, stability and good convergence property.
基金This work was supported by the National Natural Science Foundation of China （No. 50275094）.
文摘A meshless approach, called the rigid-plastic reproducing kernel particle method （RKPM）, is presented for three-dimensional （3D） bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM （finite element method） software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.
基金Supported by TOTAL(DS-2885)the National Natural Science Foundation of China(91434201,21422608)the “Strategic Priority Research Program” of the Chinese Academy of Sciences(XDA07080000)
文摘The cohesive solids in liquid flows are featured by the dynamic growth and breakage of agglomerates, and the difficulties in the development, design and optimization of these systems are related to this significant feature.In this paper, discrete particle method is used to simulate a solid–liquid flow system including millions of cohesive particles, the growth rate and breakage rate of agglomerates are then systematically investigated. It was found that the most probable size of the agglomerates is determined by the balance of growth and breakage of the agglomerates the cross point of the lines of growth rate and breakage rate as a function of the particle numbers in an agglomerate, marks the most stable agglomerate size. The finding here provides a feasible way to quantify the dynamic behaviors of growth and breakage of agglomerates, and therefore offers the possibility of quantifying the effects of agglomerates on the hydrodynamics of fluid flows with cohesive particles.
基金This work was supported by the National Key Research and Development Program of China[Grant No.2016YFC0800200]the National Natural Science Foundation of China[Grant Nos.51778568,51908492,and 52008366]+1 种基金Zhejiang Provincial Natural Science Foundation of China[Grant Nos.LQ21E080019 and LY21E080022]This work was also sup-ported by the Key Laboratory of Space Structures of Zhejiang Province(Zhejiang University)and the Center for Balance Architecture of Zhejiang University.
文摘Large deformation contact problems generally involve highly nonlinear behaviors,which are very time-consuming and may lead to convergence issues.The finite particle method(FPM)effectively separates pure deformation from total motion in large deformation problems.In addition,the decoupled procedures of the FPM make it suitable for parallel computing,which may provide an approach to solve time-consuming issues.In this study,a graphics processing unit(GPU)-based parallel algorithm is proposed for two-dimensional large deformation contact problems.The fundamentals of the FPM for planar solids are first briefly introduced,including the equations of motion of particles and the internal forces of quadrilateral elements.Subsequently,a linked-list data structure suitable for parallel processing is built,and parallel global and local search algorithms are presented for contact detection.The contact forces are then derived and directly exerted on particles.The proposed method is implemented with main solution procedures executed in parallel on a GPU.Two verification problems comprising large deformation frictional contacts are presented,and the accuracy of the proposed algorithm is validated.Furthermore,the algorithm’s performance is investigated via a large-scale contact problem,and the maximum speedups of total computational time and contact calculation reach 28.5 and 77.4,respectively,relative to commercial finite element software Abaqus/Explicit running on a single-core central processing unit(CPU).The contact calculation time percentage of the total calculation time is only 18%with the FPM,much smaller than that(50%)with Abaqus/Explicit,demonstrating the efficiency of the proposed method.
基金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.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Special Fund for Basic Scientific Research of Central Colleges of Chang’an University, China (Grant No. CHD2011JC080)
文摘In this paper, the complex variable reproducing kernel particle （CVRKP） method and the finite element （FE） method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.
基金Project /s50674002 supported by the National Natural Science Foundation of China
文摘The Reproducing Kernel Particle Method (RKPM) is one of several new meshless numerical methods de- veloped internationally in recent years. The ideal elasto-plastic constitutive model of material under a Taylor impact is characterized by the Jaumann stress- and strain-rates. An updated Lagrangian format is used for the calculation in a nu- merical analysis. With the RKPM, this paper deals with the calculation model for the Taylor impact and deduces the control equation for the impact process. A program was developed to simulate numerically the Taylor impact of projec- tiles composed of several kinds of material. The simulation result is in good accordance with both the test results and the Taylor analysis outcome. Since the meshless method is not limited by meshes, it is believed to be widely applicable to such complicated processes as the Taylor impact, including large deformation and strain and to the study of the dy- namic qualities of materials.
基金supported by the National Natural Sci-ence Foundation of China(Grant Nos.21978295,22078330,92034302 and 91834303)Innovation Academy for Green Manufacture,Chinese Academy of Sciences(Grant Nos.IAGM-2019-A03 and IAGM-2019-A13)+2 种基金Key Research Program of Frontier Sciences,Chinese Academy of Sciences(Grant No.QYZDJ-SSWJSC029)“Transformational Technologies for Clean Energy and Demonstration”Strategic Prior-ity Research Program of the Chinese Academy of Sciences(Grant No.XDA21030700)the Youth Innovation Promotion Association,Chinese Academy of Sciences(Grant No.2019050).
文摘Most natural resources are processed as particle-fluid multiphase systems in chemical,mineral and material indus-tries,therefore,discrete particles methods(DPM)are reasonable choices of simulation method for engineering the relevant processes and equipments.However,direct application of these methods is challenged by the complex multiscale behavior of such systems,which leads to enormous computational cost or otherwise qualitatively inac-curate description of the mesoscale structures.The coarse-grained DPM based on the energy-minimization multi-scale(EMMS)model,or EMMS-DPM,was proposed to reduce the computational cost by several orders while main-taining an accurate description of the mesoscale structures,which paves the way for its engineering applications.Further empowered by the high-efficiency multi-scale DEM software DEMms and the corresponding customized heterogeneous supercomputing facilities with graphics processing units(GPUs),it may even approach realtime simulation of industrial reactors.This short review will introduce the principle of DPM,in particular,EMMS-DPM,and the recent developments in modeling,numerical implementation and application of large-scale DPM which aims to reach industrial scale on one hand and resolves mesoscale structures critical to reaction-transport coupling on the other hand.This review finally prospects on the future developments of DPM in this direction.
基金Item Sponsored by National Natural Science Foundation of China(50474016)
文摘During splitting rolling simulation, re-meshing is necessary to prevent the effect of severe mesh distortion when the conventional finite element method is used. However, extreme deformation cannot be solved by the finite element method in splitting rolling. The reproducing kernel particle method can solve this problem because the continuum body is discretized by a set of nodes, and a finite element mesh is unnecessary, and there is no explicit limitation of mesh when the metal is split. To ensure stability in the large deformation elastoplastic analysis, the Lagrange material shape function was introduced. The transformation method was utilized to impose the essential boundary conditions. The splitting rolling method was simulated and the simulation results were in accordance with the experimental ones in the literature.
基金supported by the Science Council under grant No.NSC94-2611-E-002-016
文摘The moving particle semi-implicit（MPS）method has demonstrated its usefulness in practical engineering applications.Although it has wide applicability,it is still hard to predict the pressure precisely using the MPS method.A pressure-convection particle method based on the MPS method is proposed to overcome this problem.The improved performance of this new method is validated with computational and measured results.The approach is also applied to compute the problem of sloshing associated with floating body motion in waves.The pressure-convection MPS method demonstrated its capability to improve the prediction of pressure.
基金Project supported by the National Natural Science Foundation of China(Grant No.11271234)the Shandong Provincial Science Foundation,China(Grant No.ZR2017MA028)
文摘In this paper, the normal derivative of the radial basis function （RBF） is introduced into the reproducing kernel particle method （RKPM）, and the improved reproducing kernel particle method （IRKPM） is proposed. The method can decrease the errors on the boundary and improve the accuracy and stability of the algorithm. The proposed method is applied to the numerical simulation of piezoelectric materials and the corresponding governing equations are derived. The numerical results show that the IRKPM is more stable and accurate than the RKPM.
基金Foundation of Southwest Jiaotong Univer-sity (No.2005B25)
文摘To study the non-linear fracture, a non-linear constitutive model for piezoelectric ceramics was proposed, in which the polarization switching and saturation were taken into account. Based on the model, the non-linear fracture analysis was implemented using reproducing kernel particle method （RKPM）. Using local J-integral as a fracture criterion, a relation curve of fracture loads against electric fields was obtained. Qualitatively, the curve is in agreement with the experimental observations reported in literature. The reproducing equation, the shape function of RKPM, and the transformation method to impose essential boundary conditions for meshless methods were also introduced. The computation was implemented using object-oriented programming method.
文摘In recent years,the network continues to enter people’s lives,followed by network security issues that continue to appear,causing substantial economic losses to the world.As an effective method to tackle the network security issues,intrusion detection system has been widely used and studied.In this paper,the NSL-KDD data set is used to reduce the dimension of data features,remove the features of low correlation and high interference,and improve the computational efficiency.To improve the detection rate and accuracy of intrusion detection,this paper introduces the particle method for the first time that we call it intrusion detection with particle(IDP).To illustrate the effectiveness of this method,experiments are carried out on three kinds of data-before dimension reduction,after dimension reduction and importing particle method based on dimension reduction.By comparing the results of DT,NN,SVM,K-NN,and NB,it is proved that the particle method can effectively improve the intrusion detection rate.
基金the National Natural Science Foundation of China(grant No.22278404),and the Science Fund for Creative Research Groups of the National Natural Science Foundation of China(grant No.21921005)the National Key Research and Development Program of China(grant No.2023YFC2908002)the State Key Laboratory of Multiphase Complex Systems(grant No.MESO-23-A03).
文摘Stable and controllable solid flow is essential in circulating fluidized bed (CFB) systems. The L-valve is a typical non-mechanical valve that can provide flexible solid feeding. The investigation of the solid circulation rate and the hydrodynamic characteristics of the L-valve is crucial to its design and operation. The gas-solid flow in the L-valve of a full-loop CFB is studied with the coarse-grained discrete particle method (EMMS-DPM). Good agreements on the solid circulation rate and the pressure drop through the L-valve are achieved between the simulated and experimental data. The solid circulation rate increases linearly with the aeration velocity until the stable particle circulation of the CFB is destroyed. The flow patterns in the horizontal section of L-valve are gas-solid slug flow above the stationary solid layer and the moving solid layer, respectively. The effects of L-valve geometric parameters on the solid flow characteristics are also investigated. The results indicate that reducing the diameter and length of the horizontal section of L-valve can improve the solid transport efficiency, especially at low aeration velocity. Besides, the solid conveying capacity and flow stability are improved when the sharp bend of L-valve is modified to be a gradual bend.
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金supported by the National Natural Science Foundation of China (No. 50638050)the National High-Tech R&D (863) Program (No. 2007AA04Z441), China
文摘This study presents a structural analysis algorithm called the finite particle method （FPM） for kinematically indeterminate bar assemblies. Different from the traditional analysis method, FPM is based on the combination of the vector mechanics and numerical calculations. It models the analyzed domain composed of finite particles. Newton＇s second law is adopted to describe the motions of all particles. A convected material flame and explicit time integration for the solution procedure is also adopted in this method. By using the FPM, there is no need to solve any nonlinear equations, to calculate the stiffness matrix or equilibrium matrix, which is very helpful in the analysis of kinematically indeterminate structures. The basic formulations for the space bar are derived, following its solution procedures for bar assemblies. Three numerical examples are analyzed using the FPM. Results obtained from both the straight pretension cable and the suspension cable assembly show that the FPM can produce a more accurate analysis result. The motion simulation of the four-bar space assembly demonstrates the capability of this method in the analysis ofkinematically indeterminate structures.
文摘A mathematical model has been formulated based on the combined continuous and discrete particle method for investigating the sedimentation behaviour of microparticles in aqueous suspensions, by treating the fluid phase as continuous and the particles phase as discrete, thus allowing the behaviour of individual particles to be followed and the evolution of the structure of the particle phase to be investigated as a function of time. The model takes into account most of the prevailing forces acting on individual particles including van der Waals attractive, electrostatic repulsive, gravitational, Brownian, depletion, steric, contact and drag forces. A code has also been developed based on the model. This paper reports some preliminary modelling results of mono-dispersed microparticles settling in aqueous suspensions under various conditions. The results show the short time dynamics of the fluid phase, which has a similar order of magnitude to the particle phase. Such short time dynamics could bear significance to processes such as particle aggregation when their size becomes very small. Preliminary analyses of the results have also been carried out on the evolution of particle settling based on a newly proposed parameter, local normalised volume fraction （LNVF）.
