The Vortex Particle Method(VPM)is a meshless Lagrangian vortex method.Its low numerical dissipation is exceptionally suitable for wake simulation.Nevertheless,the inadequate numerical stability of VPM prevents its wid...The Vortex Particle Method(VPM)is a meshless Lagrangian vortex method.Its low numerical dissipation is exceptionally suitable for wake simulation.Nevertheless,the inadequate numerical stability of VPM prevents its widespread application in high Reynolds number flow and shear turbulence.To better simulate these flows,this paper proposes the stability-enhanced VPM based on a Reformulated VPM(RVPM)constrained by conservation of angular momentum,integrating a relaxation scheme to suppress the divergence of the vorticity field,and further coupling the Sub-Grid Scale(SGS)model to account for the turbulence dissipation caused by vortex advection and vortex stretching.The validity of the RVPM is confirmed by simulating an isolated vortex ring's evolution.The results also demonstrate that the relaxation scheme of vorticity enhances the numerical stability of the VPM by mitigating the divergence of the vorticity field.The leapfrogging vortex rings simulation demonstrates that the RVPM with the present SGS model can more precisely feature the leapfrog and fusion of vortex rings and has improved numerical stability in high Reynolds number flows.The round turbulent jet simulation confirms that the stability-enhanced VPM can stably simulate shear turbulence and accurately resolve fluctuating components and Reynolds stresses in the turbulence.展开更多
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-trian...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) tech-niques, 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 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.展开更多
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.展开更多
The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte ...The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.展开更多
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.展开更多
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.展开更多
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.展开更多
The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at t...The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at the mesoscopic scale, such as loading frequency, stress amplitude, mean stress, confining pressure and loading sequence, are also investigated with sinusoidal waveform in detail. The related numerical results have demonstrated that： 1） the loading waveform has a certain effect on rock failure processes. The square waveform has the most damage within these waveforms, while the triangle waveform has less damage than sinusoidal waveform. In each cycle, the number of microscopic cracks increases in the loading stage, while it keeps nearly constant in the unloading stage. 2） The loading frequency, stress amplitude, mean stress, confining pressure and loading sequence have considerable effects on rock damage subjected to cyclic loading. The higher the loading frequency, stress amplitude and mean stress, the greater the damage the rock accumulated; in contrast, the lower the confining pressure, the greater the damage the rock has accumulated. 3） There is a threshold value of mean stress and stress amplitude, below which no further damage accumulated after the first few cycle loadings. 4） The high-to-low loading sequence has more damage than the low-to-high loading sequence, suggesting that the rock damage is loading-path dependent.展开更多
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.展开更多
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.展开更多
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.展开更多
Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian g...Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian gridless particle method, and investigated its performance and stability to simulate large deformation of geomaterials. A calculation method was developed using geomaterials modeled as Bingham fluids to improve the original MPS method and enhance its stability. Two numerical tests showed that results from the improved MPS method was in good agreement with the theoretical value.Furthermore, numerical simulations were calibrated by laboratory experiments. It showed that the simulation results matched well with the experimentally observed free-surface configurations for flowing sand. In addition, the model could generally predict the time-history of the impact force. The MPS method could be a useful tool to evaluate large deformation of geomaterials.展开更多
A robust airfoil optimization platform is constructed based on the modified particle swarm optimization method （i.e., the second-order oscillating particle swarm method）, which consists of an efficient optimization ...A robust airfoil optimization platform is constructed based on the modified particle swarm optimization method （i.e., the second-order oscillating particle swarm method）, which consists of an efficient optimization algorithm, a precise aerodynamic analysis program, a high accuracy surrogate model, and a classical airfoil parametric method. There are two improvements for the modified particle swarm method compared with the standard particle swarm method. First, the particle velocity is represented by the combination of the particle position and the variation of position, which makes the particle swarm algorithm a second-order precision method with respect to the particle po- sition. Second, for the sake of adding diversity to the swarm and enlarging the parameter searching domain to improve the global convergence performance of the algorithm, an oscillating term is introduced to the update formula of the particle velocity. At last, tak- ing two airfoils as examples, the aerodynamic shapes are optimized on this optimization platform. It is shown from the optimization results that the aerodynamic characteristic of the airfoils is greatly improved in a broad design range.展开更多
Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and ther...Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and thereby are only applicable only to simple,single,or multiple degree-of-freedom structures.The current approaches to optimization procedures take a specific damper with its properties and observe the effect of applying time history data to the building;however,there are many different dampers and isolators that can be used.Furthermore,there is a lack of studies regarding the optimum location for various viscous and wall dampers.The main aim of this study is hybridization of the particle swarm optimization(PSO) and gravitational search algorithm(GSA) to optimize the performance of earthquake energy dissipation systems(i.e.,damper devices) simultaneously with optimizing the characteristics of the structure.Four types of structural dampers device are considered in this study:(ⅰ) variable stiffness bracing(VSB) system,(ⅱ) rubber wall damper(RWD),(ⅲ) nonlinear conical spring bracing(NCSB) device,(iv) and multi-action stiffener(MAS) device.Since many parameters may affect the design of seismic resistant structures,this study proposes a hybrid of PSO and GSA to develop a hybrid,multi-objective optimization method to resolve the aforementioned problems.The characteristics of the above-mentioned damper devices as well as the section size for structural beams and columns are considered as variables for development of the PSO-GSA optimization algorithm to minimize structural seismic response in terms of nodal displacement(in three directions) as well as plastic hinge formation in structural members simultaneously with the weight of the structure.After that,the optimization algorithm is implemented to identify the best position of the damper device in the structural frame to have the maximum effect and minimize the seismic structure response.To examine the performance of the proposed PSO-GSA optimization method,it has been applied to a three-story reinforced structure equipped with a seismic damper device.The results revealed that the method successfully optimized the earthquake energy dissipation systems and reduced the effects of earthquakes on structures,which significantly increase the building’s stability and safety during seismic excitation.The analysis results showed a reduction in the seismic response of the structure regarding the formation of plastic hinges in structural members as well as the displacement of each story to approximately 99.63%,60.5%,79.13% and 57.42% for the VSB device,RWD,NCSB device,and MAS device,respectively.This shows that using the PSO-GSA optimization algorithm and optimized damper devices in the structure resulted in no structural damage due to earthquake vibration.展开更多
The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze th...The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze the combustion and expansion processes of the fresh particles, and the thermodynamic cycle process of CRDE. In a 3D CRDE flow field, as the radius of the annulus increases, the no-injection area proportion increases, the non-detonation proportion decreases, and the detonation height decreases. The flow field parameters on the 3D mid annulus are different from in the 2D flow field under the same chamber size. The non-detonation proportion in the 3D flow field is less than in the 2D flow field. In the 2D and 3D CRDE, the paths of the flow particles have only a small fluctuation in the circumferential direction. The numerical thermodynamic cycle processes are qualitatively consistent with the three ideal cycle models, and they are right in between the ideal F–J cycle and ideal ZND cycle. The net mechanical work and thermal efficiency are slightly smaller in the 2D simulation than in the 3D simulation. In the 3D CRDE, as the radius of the annulus increases, the net mechanical work is almost constant, and the thermal efficiency increases. The numerical thermal efficiencies are larger than F–J cycle, and much smaller than ZND cycle.展开更多
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.展开更多
In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avo...In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.展开更多
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.展开更多
基金co-supported by the National Natural Science Foundation of China(No.12402272)the Natural Science Basic Research Program of Shaanxi Province,China(No.2024JC-YBQN-0024)the Fundamental Research Funds for the Central Universities,China(No.D5000240030)。
文摘The Vortex Particle Method(VPM)is a meshless Lagrangian vortex method.Its low numerical dissipation is exceptionally suitable for wake simulation.Nevertheless,the inadequate numerical stability of VPM prevents its widespread application in high Reynolds number flow and shear turbulence.To better simulate these flows,this paper proposes the stability-enhanced VPM based on a Reformulated VPM(RVPM)constrained by conservation of angular momentum,integrating a relaxation scheme to suppress the divergence of the vorticity field,and further coupling the Sub-Grid Scale(SGS)model to account for the turbulence dissipation caused by vortex advection and vortex stretching.The validity of the RVPM is confirmed by simulating an isolated vortex ring's evolution.The results also demonstrate that the relaxation scheme of vorticity enhances the numerical stability of the VPM by mitigating the divergence of the vorticity field.The leapfrogging vortex rings simulation demonstrates that the RVPM with the present SGS model can more precisely feature the leapfrog and fusion of vortex rings and has improved numerical stability in high Reynolds number flows.The round turbulent jet simulation confirms that the stability-enhanced VPM can stably simulate shear turbulence and accurately resolve fluctuating components and Reynolds stresses in the turbulence.
基金Project supported by the National Natural Science Foundation of China (No.10202018)
文摘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) tech-niques, 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.
基金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.
基金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.
文摘The low diffusion （LD） particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo （DSMC） method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.
基金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.
基金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.
基金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.
基金Projects(11702235,51641905,41472269) supported by the National Natural Science Foundation of ChinaProject(2017JJ3290) supported by the Natural Science Foundation of Hunan Province,China+1 种基金Project(17C1540) supported by the Scientific Research Foundation of Education Department of Hunan Province,ChinaProject(16GES07) supported by the Open Research Fund of Hunan Key Laboratory of Geomechanics and Engineering Safety,China
文摘The particle simulation method is used to study the effects of loading waveforms （i.e. square, sinusoidal and triangle waveforms） on rock damage at mesoscopic scale. Then some influencing factors on rock damage at the mesoscopic scale, such as loading frequency, stress amplitude, mean stress, confining pressure and loading sequence, are also investigated with sinusoidal waveform in detail. The related numerical results have demonstrated that： 1） the loading waveform has a certain effect on rock failure processes. The square waveform has the most damage within these waveforms, while the triangle waveform has less damage than sinusoidal waveform. In each cycle, the number of microscopic cracks increases in the loading stage, while it keeps nearly constant in the unloading stage. 2） The loading frequency, stress amplitude, mean stress, confining pressure and loading sequence have considerable effects on rock damage subjected to cyclic loading. The higher the loading frequency, stress amplitude and mean stress, the greater the damage the rock accumulated; in contrast, the lower the confining pressure, the greater the damage the rock has accumulated. 3） There is a threshold value of mean stress and stress amplitude, below which no further damage accumulated after the first few cycle loadings. 4） The high-to-low loading sequence has more damage than the low-to-high loading sequence, suggesting that the rock damage is loading-path dependent.
基金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.
基金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.
基金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.
文摘Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian gridless particle method, and investigated its performance and stability to simulate large deformation of geomaterials. A calculation method was developed using geomaterials modeled as Bingham fluids to improve the original MPS method and enhance its stability. Two numerical tests showed that results from the improved MPS method was in good agreement with the theoretical value.Furthermore, numerical simulations were calibrated by laboratory experiments. It showed that the simulation results matched well with the experimentally observed free-surface configurations for flowing sand. In addition, the model could generally predict the time-history of the impact force. The MPS method could be a useful tool to evaluate large deformation of geomaterials.
文摘A robust airfoil optimization platform is constructed based on the modified particle swarm optimization method （i.e., the second-order oscillating particle swarm method）, which consists of an efficient optimization algorithm, a precise aerodynamic analysis program, a high accuracy surrogate model, and a classical airfoil parametric method. There are two improvements for the modified particle swarm method compared with the standard particle swarm method. First, the particle velocity is represented by the combination of the particle position and the variation of position, which makes the particle swarm algorithm a second-order precision method with respect to the particle po- sition. Second, for the sake of adding diversity to the swarm and enlarging the parameter searching domain to improve the global convergence performance of the algorithm, an oscillating term is introduced to the update formula of the particle velocity. At last, tak- ing two airfoils as examples, the aerodynamic shapes are optimized on this optimization platform. It is shown from the optimization results that the aerodynamic characteristic of the airfoils is greatly improved in a broad design range.
基金University Putra Malaysia under Putra Grant No.9531200。
文摘Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and thereby are only applicable only to simple,single,or multiple degree-of-freedom structures.The current approaches to optimization procedures take a specific damper with its properties and observe the effect of applying time history data to the building;however,there are many different dampers and isolators that can be used.Furthermore,there is a lack of studies regarding the optimum location for various viscous and wall dampers.The main aim of this study is hybridization of the particle swarm optimization(PSO) and gravitational search algorithm(GSA) to optimize the performance of earthquake energy dissipation systems(i.e.,damper devices) simultaneously with optimizing the characteristics of the structure.Four types of structural dampers device are considered in this study:(ⅰ) variable stiffness bracing(VSB) system,(ⅱ) rubber wall damper(RWD),(ⅲ) nonlinear conical spring bracing(NCSB) device,(iv) and multi-action stiffener(MAS) device.Since many parameters may affect the design of seismic resistant structures,this study proposes a hybrid of PSO and GSA to develop a hybrid,multi-objective optimization method to resolve the aforementioned problems.The characteristics of the above-mentioned damper devices as well as the section size for structural beams and columns are considered as variables for development of the PSO-GSA optimization algorithm to minimize structural seismic response in terms of nodal displacement(in three directions) as well as plastic hinge formation in structural members simultaneously with the weight of the structure.After that,the optimization algorithm is implemented to identify the best position of the damper device in the structural frame to have the maximum effect and minimize the seismic structure response.To examine the performance of the proposed PSO-GSA optimization method,it has been applied to a three-story reinforced structure equipped with a seismic damper device.The results revealed that the method successfully optimized the earthquake energy dissipation systems and reduced the effects of earthquakes on structures,which significantly increase the building’s stability and safety during seismic excitation.The analysis results showed a reduction in the seismic response of the structure regarding the formation of plastic hinges in structural members as well as the displacement of each story to approximately 99.63%,60.5%,79.13% and 57.42% for the VSB device,RWD,NCSB device,and MAS device,respectively.This shows that using the PSO-GSA optimization algorithm and optimized damper devices in the structure resulted in no structural damage due to earthquake vibration.
文摘The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze the combustion and expansion processes of the fresh particles, and the thermodynamic cycle process of CRDE. In a 3D CRDE flow field, as the radius of the annulus increases, the no-injection area proportion increases, the non-detonation proportion decreases, and the detonation height decreases. The flow field parameters on the 3D mid annulus are different from in the 2D flow field under the same chamber size. The non-detonation proportion in the 3D flow field is less than in the 2D flow field. In the 2D and 3D CRDE, the paths of the flow particles have only a small fluctuation in the circumferential direction. The numerical thermodynamic cycle processes are qualitatively consistent with the three ideal cycle models, and they are right in between the ideal F–J cycle and ideal ZND cycle. The net mechanical work and thermal efficiency are slightly smaller in the 2D simulation than in the 3D simulation. In the 3D CRDE, as the radius of the annulus increases, the net mechanical work is almost constant, and the thermal efficiency increases. The numerical thermal efficiencies are larger than F–J cycle, and much smaller than ZND cycle.
基金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.
基金the National Natural Science Foundation of China(Grant No.41807223)the Fundamental Research Funds for the Central Universities(Grant No.B210202096)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDA 23090202).
文摘In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.
基金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.
