We present and explore a new shock-capturing particle hydrodynamics approach.Our starting point is a commonly used discretization of smoothed particle hydrodynamics.We enhance this discretization with Roe’s approx-im...We present and explore a new shock-capturing particle hydrodynamics approach.Our starting point is a commonly used discretization of smoothed particle hydrodynamics.We enhance this discretization with Roe’s approx-imate Riemann solver,we identify its dissipative terms,and in these terms,we use slope-limited linear reconstruction.All gradients needed for our method are calculated with linearly reproducing kernels that are constructed to enforce the two lowest-order consistency relations.We scrutinize our reproducing kernel implementation carefully on a“glass-like”particle distribution,and we find that constant and linear functions are recovered to machine precision.We probe our method in a series of challenging 3D benchmark problems ranging from shocks over instabilities to Schulz-Rinne-type vorticity-creating shocks.All of our simulations show excellent agreement with analyticeference solutions.展开更多
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.展开更多
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.展开更多
Deep learning algorithms based on neural networks make remarkable achievements in machine fault diagnosis,while the noise mixed in measured signals harms the prediction accuracy of networks.Existing denoising methods ...Deep learning algorithms based on neural networks make remarkable achievements in machine fault diagnosis,while the noise mixed in measured signals harms the prediction accuracy of networks.Existing denoising methods in neural networks,such as using complex network architectures and introducing sparse techniques,always suffer from the difficulty of estimating hyperparameters and the lack of physical interpretability.To address this issue,this paper proposes a novel interpretable denoising layer based on reproducing kernel Hilbert space(RKHS)as the first layer for standard neural networks,with the aim to combine the advantages of both traditional signal processing technology with physical interpretation and network modeling strategy with parameter adaption.By investigating the influencing mechanism of parameters on the regularization procedure in RKHS,the key parameter that dynamically controls the signal smoothness with low computational cost is selected as the only trainable parameter of the proposed layer.Besides,the forward and backward propagation algorithms of the designed layer are formulated to ensure that the selected parameter can be automatically updated together with other parameters in the neural network.Moreover,exponential and piecewise functions are introduced in the weight updating process to keep the trainable weight within a reasonable range and avoid the ill-conditioned problem.Experiment studies verify the effectiveness and compatibility of the proposed layer design method in intelligent fault diagnosis of machinery in noisy environments.展开更多
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.展开更多
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.展开更多
This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)metho...This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.展开更多
This work develops a Hermitian C2 differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredepend...This work develops a Hermitian C2 differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C2 differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC2DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C2 DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.展开更多
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.展开更多
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.展开更多
Tomodel amultibody systemcomposed of shell components,a geometrically exact Kirchho-Love triangular shell element is proposed.The middle surface of the shell element is described by using the DMS-splines,which can exa...Tomodel amultibody systemcomposed of shell components,a geometrically exact Kirchho-Love triangular shell element is proposed.The middle surface of the shell element is described by using the DMS-splines,which can exactly represent arbitrary topology piecewise polynomial triangular surfaces.The proposed shell element employs only nodal displacement and can automatically maintain C1 continuity properties at the element boundaries.A reproducing DMS-spline kernel skill is also introduced to improve computation stability and accuracy.The proposed triangular shell element based on reproducing kernel DMS-splines can achieve an almost optimal convergent rate.Finally,the proposed shell element is validated via three static problems of shells and the dynamic simulation of aexible multibody system undergoing both overall motions and large deformations.展开更多
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 finit...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.展开更多
A reproducing kernel collocation method based on strong formulation is introduced for transient dynamics. To study the stability property of this method, an algorithm based on the von Neumann hypothesis is proposed to...A reproducing kernel collocation method based on strong formulation is introduced for transient dynamics. To study the stability property of this method, an algorithm based on the von Neumann hypothesis is proposed to predict the critical time step. A numerical test is conducted to validate the algorithm. The numerical critical time step and the predicted critical time step are in good agreement. The results are compared with those obtained based on the radial basis collocation method, and they axe in good agreement. Several important conclusions for choosing a proper support size of the reproducing kernel shape function are given to improve the stability condition.展开更多
By combining the wavelet decomposition with kernel method, a practical approach of universal multiscale wavelet kernels constructed in reproducing kernel Hilbert space (RKHS) is discussed, and an identification sche...By combining the wavelet decomposition with kernel method, a practical approach of universal multiscale wavelet kernels constructed in reproducing kernel Hilbert space (RKHS) is discussed, and an identification scheme using wavelet support vector machines (WSVM) estimator is proposed for nordinear dynamic systems. The good approximating properties of wavelet kernel function enhance the generalization ability of the proposed method, and the comparison of some numerical experimental results between the novel approach and some existing methods is encouraging.展开更多
How to solve the partial differential equation has been attached importance to by all kinds of fields. The exact solution to a class of partial differential equation with variable-coefficient is obtained in reproducin...How to solve the partial differential equation has been attached importance to by all kinds of fields. The exact solution to a class of partial differential equation with variable-coefficient is obtained in reproducing kernel space. For getting the approximate solution, give an iterative method, convergence of the iterative method is proved. The numerical example shows that our method is effective and good practicability.展开更多
Many mechanical problems can be induced from differential equations with boundary conditions; there exist analytic and numerical methods for solving the differential equations. Usually it is not so easy to obtain anal...Many mechanical problems can be induced from differential equations with boundary conditions; there exist analytic and numerical methods for solving the differential equations. Usually it is not so easy to obtain analytic solutions. So it is necessary to give numerical solutions. The reproducing kernel particle (RKP) method is based on the Carlerkin Meshless method. According to the Sobolev space and Fourier transform, the RKP shape function is mathematically proved in this paper.展开更多
The spherical approximation between two nested reproducing kernels Hilbert spaces generated from different smooth kernels is investigated. It is shown that the functions of a space can be approximated by that of the s...The spherical approximation between two nested reproducing kernels Hilbert spaces generated from different smooth kernels is investigated. It is shown that the functions of a space can be approximated by that of the subspace with better smoothness. Furthermore, the upper bound of approximation error is given.展开更多
It is well-known that using the traditional reproducing kernel method(TRKM) for solving the fractional partial differential equation(FPDE) is very intractable. In this study, the adaptive single piecewise interpolatio...It is well-known that using the traditional reproducing kernel method(TRKM) for solving the fractional partial differential equation(FPDE) is very intractable. In this study, the adaptive single piecewise interpolation reproducing kernel method(ASPIRKM) is determined to solve the FPDE. This improved method not only improves the calculation accuracy, but also reduces the waste of time. Two numerical examples show that the ASPIRKM is a more time-saving numerical method than the TRKM.展开更多
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 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.展开更多
基金supported by the Swedish Research Council(VR)under grant number 2020-05044by the research environment grant"Gravitational Radiation and Electromagnetic Astrophysical Transients"(GREAT)funded by the Swedish Research Council(VR)under Dnr 2016-06012+2 种基金by the Knut and Alice Wallenberg Foundation under grant Dnr.KAW 2019.0112by the Deutsche Forschungsgemeinschaft(DFG,German Research Foundation)under Germany's Excellence Strategy-EXC 2121"Quantum Universe"-390833306by the European Research Council(ERC)Advanced Grant INSPIRATION under the European Union's Horizon 2020 Research and Innovation Programme(Grant agreement No.101053985).
摘要We present and explore a new shock-capturing particle hydrodynamics approach.Our starting point is a commonly used discretization of smoothed particle hydrodynamics.We enhance this discretization with Roe’s approx-imate Riemann solver,we identify its dissipative terms,and in these terms,we use slope-limited linear reconstruction.All gradients needed for our method are calculated with linearly reproducing kernels that are constructed to enforce the two lowest-order consistency relations.We scrutinize our reproducing kernel implementation carefully on a“glass-like”particle distribution,and we find that constant and linear functions are recovered to machine precision.We probe our method in a series of challenging 3D benchmark problems ranging from shocks over instabilities to Schulz-Rinne-type vorticity-creating shocks.All of our simulations show excellent agreement with analyticeference solutions.
基金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.
基金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 National Natural Science Foundation of China(Grant Nos.12072188,11632011,11702171,11572189,51121063)Shanghai Municipal Natural Science Foundation of China(Grant No.20ZR1425200).
摘要Deep learning algorithms based on neural networks make remarkable achievements in machine fault diagnosis,while the noise mixed in measured signals harms the prediction accuracy of networks.Existing denoising methods in neural networks,such as using complex network architectures and introducing sparse techniques,always suffer from the difficulty of estimating hyperparameters and the lack of physical interpretability.To address this issue,this paper proposes a novel interpretable denoising layer based on reproducing kernel Hilbert space(RKHS)as the first layer for standard neural networks,with the aim to combine the advantages of both traditional signal processing technology with physical interpretation and network modeling strategy with parameter adaption.By investigating the influencing mechanism of parameters on the regularization procedure in RKHS,the key parameter that dynamically controls the signal smoothness with low computational cost is selected as the only trainable parameter of the proposed layer.Besides,the forward and backward propagation algorithms of the designed layer are formulated to ensure that the selected parameter can be automatically updated together with other parameters in the neural network.Moreover,exponential and piecewise functions are introduced in the weight updating process to keep the trainable weight within a reasonable range and avoid the ill-conditioned problem.Experiment studies verify the effectiveness and compatibility of the proposed layer design method in intelligent fault diagnosis of machinery in noisy environments.
摘要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.
基金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.
基金the National Natural Science Foundation of China(Grant Nos.71961022,11902163,12265020,and 12262024)the Natural Science Foundation of Inner Mongolia Autonomous Region of China(Grant Nos.2019BS01011 and 2022MS01003)+5 种基金2022 Inner Mongolia Autonomous Region Grassland Talents Project-Young Innovative and Entrepreneurial Talents(Mingjing Du)2022 Talent Development Foundation of Inner Mongolia Autonomous Region of China(Ming-Jing Du)the Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region Program(Grant No.NJYT-20-B18)the Key Project of High-quality Economic Development Research Base of Yellow River Basin in 2022(Grant No.21HZD03)2022 Inner Mongolia Autonomous Region International Science and Technology Cooperation High-end Foreign Experts Introduction Project(Ge Kai)MOE(Ministry of Education in China)Humanities and Social Sciences Foundation(Grants No.20YJC860005).
摘要This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.
基金supported by a grant from the National Science and Technology Council of the Republic of China(Grant Number:MOST 112-2221-E-006-048-MY2).
摘要This work develops a Hermitian C2 differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C2 differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC2DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C2 DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.
基金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.
基金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 in part by the National Natural Science Foundations of China under Grants 11290151,11672034 and 11902363。
摘要Tomodel amultibody systemcomposed of shell components,a geometrically exact Kirchho-Love triangular shell element is proposed.The middle surface of the shell element is described by using the DMS-splines,which can exactly represent arbitrary topology piecewise polynomial triangular surfaces.The proposed shell element employs only nodal displacement and can automatically maintain C1 continuity properties at the element boundaries.A reproducing DMS-spline kernel skill is also introduced to improve computation stability and accuracy.The proposed triangular shell element based on reproducing kernel DMS-splines can achieve an almost optimal convergent rate.Finally,the proposed shell element is validated via three static problems of shells and the dynamic simulation of aexible multibody system undergoing both overall motions and large deformations.
基金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.
基金Project supported by the Western Transport Technical Project of Ministry of Transport of China(No. 2009318000046)
摘要A reproducing kernel collocation method based on strong formulation is introduced for transient dynamics. To study the stability property of this method, an algorithm based on the von Neumann hypothesis is proposed to predict the critical time step. A numerical test is conducted to validate the algorithm. The numerical critical time step and the predicted critical time step are in good agreement. The results are compared with those obtained based on the radial basis collocation method, and they axe in good agreement. Several important conclusions for choosing a proper support size of the reproducing kernel shape function are given to improve the stability condition.
基金the National 973 Key Fundamental Research Project of China (Grant No.2002CB312200)
摘要By combining the wavelet decomposition with kernel method, a practical approach of universal multiscale wavelet kernels constructed in reproducing kernel Hilbert space (RKHS) is discussed, and an identification scheme using wavelet support vector machines (WSVM) estimator is proposed for nordinear dynamic systems. The good approximating properties of wavelet kernel function enhance the generalization ability of the proposed method, and the comparison of some numerical experimental results between the novel approach and some existing methods is encouraging.
基金Project supported by the National Natural Science Foundation of China(No.10461005)
摘要How to solve the partial differential equation has been attached importance to by all kinds of fields. The exact solution to a class of partial differential equation with variable-coefficient is obtained in reproducing kernel space. For getting the approximate solution, give an iterative method, convergence of the iterative method is proved. The numerical example shows that our method is effective and good practicability.
摘要Many mechanical problems can be induced from differential equations with boundary conditions; there exist analytic and numerical methods for solving the differential equations. Usually it is not so easy to obtain analytic solutions. So it is necessary to give numerical solutions. The reproducing kernel particle (RKP) method is based on the Carlerkin Meshless method. According to the Sobolev space and Fourier transform, the RKP shape function is mathematically proved in this paper.
基金the NSFC(60473034)the Science Foundation of Zhejiang Province(Y604003).
摘要The spherical approximation between two nested reproducing kernels Hilbert spaces generated from different smooth kernels is investigated. It is shown that the functions of a space can be approximated by that of the subspace with better smoothness. Furthermore, the upper bound of approximation error is given.
基金Natural Science Foundation of Inner Mongolia Autonomous Region of China (No.2019BS01011)Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region,China (No.NJYT-20-B18)2022 Talent Development Foundation of Inner Mongolia Autonomous Region,China。
摘要It is well-known that using the traditional reproducing kernel method(TRKM) for solving the fractional partial differential equation(FPDE) is very intractable. In this study, the adaptive single piecewise interpolation reproducing kernel method(ASPIRKM) is determined to solve the FPDE. This improved method not only improves the calculation accuracy, but also reduces the waste of time. Two numerical examples show that the ASPIRKM is a more time-saving numerical method than the TRKM.
基金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.
基金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.