This work develops a Hermitian C^(2) 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 C^(2) 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 C^(2) 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 ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) 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.展开更多
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.展开更多
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 ex...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.展开更多
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.展开更多
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.展开更多
An interpolating reproducing kernel particle method for two-dimensional (2D) scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in th...An interpolating reproducing kernel particle method for two-dimensional (2D) scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.展开更多
In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE metho...In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.展开更多
It is well known that the problem on the stability of the solutions for Fredholm integral equation of the first kind is an ill-posed problem in C[a, b] or L2 [a, b]. In this paper, the representation of the solution f...It is well known that the problem on the stability of the solutions for Fredholm integral equation of the first kind is an ill-posed problem in C[a, b] or L2 [a, b]. In this paper, the representation of the solution for Fredholm integral equation of the first kind is given if it has a unique solution. The stability of the solution is proved in the reproducing kernel space, namely, the measurement errors of the experimental data cannot result in unbounded errors of the true solution. The computation of approximate solution is also stable with respect to ||· ||c or ||L2· A numerical experiment shows that the method given in this paper is stable in the reproducing kernel space.展开更多
The reproducing kernel particle method (RKPM) has been efficiently applied to problems with large deformations, high gradients and high modal density. In this paper, it is extended to solve a nonlocal problem modele...The reproducing kernel particle method (RKPM) has been efficiently applied to problems with large deformations, high gradients and high modal density. In this paper, it is extended to solve a nonlocal problem modeled by a fractional advectiondiffusion equation (FADE), which exhibits a boundary layer with low regularity. We formulate this method on a moving least-square approach. Via the enrichment of fractional-order power functions to the traditional integer-order basis for RKPM, leading terms of the solution to the FADE can be exactly reproduced, which guarantees a good approximation to the boundary layer. Numerical tests are performed to verify the proposed approach.展开更多
In this paper,an efficient multi-step scheme is presented based on reproducing kernel Hilbert space(RKHS)theory for solving ordinary stiff differential systems.The solution methodology depends on reproducing kernel fu...In this paper,an efficient multi-step scheme is presented based on reproducing kernel Hilbert space(RKHS)theory for solving ordinary stiff differential systems.The solution methodology depends on reproducing kernel functions to obtain analytic solutions in a uniform formfor a rapidly convergent series in the posed Sobolev space.Using the Gram-Schmidt orthogonality process,complete orthogonal essential functions are obtained in a compact field to encompass Fourier series expansion with the help of kernel properties reproduction.Consequently,by applying the standard RKHS method to each subinterval,approximate solutions that converge uniformly to the exact solutions are obtained.For this purpose,several numerical examples are tested to show proposed algorithm’s superiority,simplicity,and efficiency.The gained results indicate that themulti-step RKHSmethod is suitable for solving linear and nonlinear stiffness systems over an extensive duration and giving highly accurate outcomes.展开更多
The main objective of this work is to decompose orthogonally the reproducing kernels Hilbert space using any conditionally positive definite kernels into smaller ones by introducing the theory of power kernels, and to...The main objective of this work is to decompose orthogonally the reproducing kernels Hilbert space using any conditionally positive definite kernels into smaller ones by introducing the theory of power kernels, and to show how to do this decomposition recur- sively. It may be used to split large interpolation problems into smaller ones with different kernels which are related to the original kernels. To reach this objective, we will reconstruct the reproducing kernels Hilbert space for the normalized and the extended kernels and give the recursive algorithm of this decomposition.展开更多
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.展开更多
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.展开更多
Consider the design problem for estimation and extrapolation in approximately linear regression models with possible misspecification. The design space is a discrete set consisting of finitely many points, and the mod...Consider the design problem for estimation and extrapolation in approximately linear regression models with possible misspecification. The design space is a discrete set consisting of finitely many points, and the model bias comes from a reproducing kernel Hilbert space. Two different design criteria are proposed by applying the minimax approach for estimating the parameters of the regression response and extrapolating the regression response to points outside of the design space. A simulated annealing algorithm is applied to construct the minimax designs. These minimax designs are compared with the classical D-optimal designs and all-bias extrapolation designs. Numerical results indicate that the simulated annealing algorithm is feasible and the minimax designs are robust against bias caused by model misspecification.展开更多
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.展开更多
During splitting rolling simulation, re-meshing is necessary to prevent the effect of severe mesh distortion when the conventional finite element method is used. However, extreme deformation cannot be solved by the fi...During splitting rolling simulation, re-meshing is necessary to prevent the effect of severe mesh distortion when the conventional finite element method is used. However, extreme deformation cannot be solved by the finite element method in splitting rolling. The reproducing kernel particle method can solve this problem because the continuum body is discretized by a set of nodes, and a finite element mesh is unnecessary, and there is no explicit limitation of mesh when the metal is split. To ensure stability in the large deformation elastoplastic analysis, the Lagrange material shape function was introduced. The transformation method was utilized to impose the essential boundary conditions. The splitting rolling method was simulated and the simulation results were in accordance with the experimental ones in the literature.展开更多
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.展开更多
In this paper, the normal derivative of the radial basis function (RBF) is introduced into the reproducing kernel particle method (RKPM), and the improved reproducing kernel particle method (IRKPM) is proposed. ...In this paper, the normal derivative of the radial basis function (RBF) is introduced into the reproducing kernel particle method (RKPM), and the improved reproducing kernel particle method (IRKPM) is proposed. The method can decrease the errors on the boundary and improve the accuracy and stability of the algorithm. The proposed method is applied to the numerical simulation of piezoelectric materials and the corresponding governing equations are derived. The numerical results show that the IRKPM is more stable and accurate than the RKPM.展开更多
基金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 C^(2) 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 C^(2) 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 ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) 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 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 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.
基金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.
基金This work was supported by the National Natural Science Foundation of China (No. 50275094).
文摘A meshless approach, called the rigid-plastic reproducing kernel particle method (RKPM), is presented for three-dimensional (3D) bulk metal forming simulation. The approach is a combination of RKPM with the flow theory of 3D rigid-plastic mechanics. For the treatments of essential boundary conditions and incompressibility constraint, the boundary singular kernel method and the modified penalty method are utilized, respectively. The arc-tangential friction model is employed to treat the contact conditions. The compression of rectangular blocks, a typical 3D upsetting operation, is analyzed for different friction conditions and the numerical results are compared with those obtained using commercial rigid-plastic FEM (finite element method) software Deform^3D. As results show, when handling 3D plastic deformations, the proposed approach eliminates the need of expensive meshing and remeshing procedures which are unavoidable in conventional FEM and can provide results that are in good agreement with finite element predictions.
基金supported by the National Natural Science Foundation of China(Grant No.11171208)the Natural Science Foundation of Shanxi Province,China(Grant No.2013011022-6)
文摘An interpolating reproducing kernel particle method for two-dimensional (2D) scatter points is introduced. It elim- inates the dependency of gridding in numerical calculations. The interpolating shape function in the interpolating repro- ducing kernel particle method satisfies the property of the Kronecker delta function. This method offers a mathematics basis for recognition technology and simulation analysis, which can be expressed as simultaneous differential equations in science or project problems. Mathematical examples are given to show the validity of the interpolating reproducing kernel particle method.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Special Fund for Basic Scientific Research of Central Colleges of Chang’an University, China (Grant No. CHD2011JC080)
文摘In this paper, the complex variable reproducing kernel particle (CVRKP) method and the finite element (FE) method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.
文摘It is well known that the problem on the stability of the solutions for Fredholm integral equation of the first kind is an ill-posed problem in C[a, b] or L2 [a, b]. In this paper, the representation of the solution for Fredholm integral equation of the first kind is given if it has a unique solution. The stability of the solution is proved in the reproducing kernel space, namely, the measurement errors of the experimental data cannot result in unbounded errors of the true solution. The computation of approximate solution is also stable with respect to ||· ||c or ||L2· A numerical experiment shows that the method given in this paper is stable in the reproducing kernel space.
基金supported partly by the National Natural Science Foundation of China (Grant 11521202)support from the Chinese Scholarship Councilpartially support by an Army Research Office (Grant W911NF-15-10569)
文摘The reproducing kernel particle method (RKPM) has been efficiently applied to problems with large deformations, high gradients and high modal density. In this paper, it is extended to solve a nonlocal problem modeled by a fractional advectiondiffusion equation (FADE), which exhibits a boundary layer with low regularity. We formulate this method on a moving least-square approach. Via the enrichment of fractional-order power functions to the traditional integer-order basis for RKPM, leading terms of the solution to the FADE can be exactly reproduced, which guarantees a good approximation to the boundary layer. Numerical tests are performed to verify the proposed approach.
文摘In this paper,an efficient multi-step scheme is presented based on reproducing kernel Hilbert space(RKHS)theory for solving ordinary stiff differential systems.The solution methodology depends on reproducing kernel functions to obtain analytic solutions in a uniform formfor a rapidly convergent series in the posed Sobolev space.Using the Gram-Schmidt orthogonality process,complete orthogonal essential functions are obtained in a compact field to encompass Fourier series expansion with the help of kernel properties reproduction.Consequently,by applying the standard RKHS method to each subinterval,approximate solutions that converge uniformly to the exact solutions are obtained.For this purpose,several numerical examples are tested to show proposed algorithm’s superiority,simplicity,and efficiency.The gained results indicate that themulti-step RKHSmethod is suitable for solving linear and nonlinear stiffness systems over an extensive duration and giving highly accurate outcomes.
文摘The main objective of this work is to decompose orthogonally the reproducing kernels Hilbert space using any conditionally positive definite kernels into smaller ones by introducing the theory of power kernels, and to show how to do this decomposition recur- sively. It may be used to split large interpolation problems into smaller ones with different kernels which are related to the original kernels. To reach this objective, we will reconstruct the reproducing kernels Hilbert space for the normalized and the extended kernels and give the recursive algorithm of this decomposition.
基金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.
文摘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.
基金Supported by National Natural Science Foundation of China(11471216,11301332)E-Institutes of Shanghai Municipal Education Commission(E03004)+1 种基金Central Finance Project(YC-XK-13105)Shanghai Municipal Science and Technology Research Project(14DZ1201902)
文摘Consider the design problem for estimation and extrapolation in approximately linear regression models with possible misspecification. The design space is a discrete set consisting of finitely many points, and the model bias comes from a reproducing kernel Hilbert space. Two different design criteria are proposed by applying the minimax approach for estimating the parameters of the regression response and extrapolating the regression response to points outside of the design space. A simulated annealing algorithm is applied to construct the minimax designs. These minimax designs are compared with the classical D-optimal designs and all-bias extrapolation designs. Numerical results indicate that the simulated annealing algorithm is feasible and the minimax designs are robust against bias caused by model misspecification.
基金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.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant No.11271234)the Shandong Provincial Science Foundation,China(Grant No.ZR2017MA028)
文摘In this paper, the normal derivative of the radial basis function (RBF) is introduced into the reproducing kernel particle method (RKPM), and the improved reproducing kernel particle method (IRKPM) is proposed. The method can decrease the errors on the boundary and improve the accuracy and stability of the algorithm. The proposed method is applied to the numerical simulation of piezoelectric materials and the corresponding governing equations are derived. The numerical results show that the IRKPM is more stable and accurate than the RKPM.