The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the bas...The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.展开更多
This paper presents a new approach to the structural topology optimization of continuum structures. Material-point independent variables are presented to illustrate the existence condition,or inexistence of the materi...This paper presents a new approach to the structural topology optimization of continuum structures. Material-point independent variables are presented to illustrate the existence condition,or inexistence of the material points and their vicinity instead of elements or nodes in popular topology optimization methods. Topological variables field is constructed by moving least square approximation which is used as a shape function in the meshless method. Combined with finite element analyses,not only checkerboard patterns and mesh-dependence phenomena are overcome by this continuous and smooth topological variables field,but also the locations and numbers of topological variables can be arbitrary. Parameters including the number of quadrature points,scaling parameter,weight function and so on upon optimum topological configurations are discussed. Two classic topology optimization problems are solved successfully by the proposed method. The method is found robust and no numerical instabilities are found with proper parameters.展开更多
Based on the moving least square （MLS） approximations and the boundary integral equations （BIEs）, a meshless algorithm is presented in this paper for elliptic Signorini problems. In the algorithm, a projection ope...Based on the moving least square （MLS） approximations and the boundary integral equations （BIEs）, a meshless algorithm is presented in this paper for elliptic Signorini problems. In the algorithm, a projection operator is used to tackle the nonlinear boundary inequality conditions. The Signorini problem is then reformulated as BIEs and the unknown boundary variables are approximated by the MLS approximations. Accordingly, only a nodal data structure on the boundary of a domain is required. The convergence of the algorithm is proven. Numerical examples are given to show the high convergence rate and high computational efficiency of the presented algorithm.展开更多
By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive...By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive linear operators. As an example, Hermite Fejér interpolation polynomial operators are analysed and studied, and a general conclusion is obtained.展开更多
In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem....In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem. the evaluation of the objective function is very difficult, so that only their approximate values can be obtained. This algorithm is obtained by combining penalty function method and approximation in bilevel programming. The presented algorithm is completely different from existing methods. That convergence for this algorithm is proved.展开更多
This paper is devoted to the study of approximation of the solution for the differential equation whose coefficients are almost period functions. To this end the authors establish the estimation of the solution of gen...This paper is devoted to the study of approximation of the solution for the differential equation whose coefficients are almost period functions. To this end the authors establish the estimation of the solution of general linear differential equation for infinite interval case. For finite interval case, this equation was investigated by G. Tamarkin([1]) applying the Picard method of successive approximation.展开更多
The independent continuous mapping（ICM） method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of ...The independent continuous mapping（ICM） method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of the essential boundary condition and derivative of various sensitivities,a singular weight function in element free Galerkin method is introduced.Material point variable is defined to illustrate the condition of material point and its vicinity instead of element or node.The topological variables field is constructed by moving least square approximation which inherits the continuity and smoothness of the weight function.Due to reciprocal relationships between the topological variables and design variables,various structural responses sensitivities are derived according to the method for calculating the partial derivatives of compound functions.Numerical examples indicate that checkerboard pattern and mesh-dependence phenomena are overcome without additional restriction methods.展开更多
This paper presents an interpolation-based method(IBM)for approximating some trigonometric functions or their integrals as well.It provides two-sided bounds for each function,which also achieves much better approximat...This paper presents an interpolation-based method(IBM)for approximating some trigonometric functions or their integrals as well.It provides two-sided bounds for each function,which also achieves much better approximation effects than those of prevailing methods.In principle,the IBM can be applied for bounding more bounded smooth functions and their integrals as well,and its applications include approximating the integral of sin(x)/x function and improving the famous square root inequalities.展开更多
Asymptotic eigenvalues and eigenfunctions for the Orr-Sommerfeld equation in two-dimensional and three-dimensional incompressible flows on an infinite domain and on a semi-infinite domain are obtained. Two configurati...Asymptotic eigenvalues and eigenfunctions for the Orr-Sommerfeld equation in two-dimensional and three-dimensional incompressible flows on an infinite domain and on a semi-infinite domain are obtained. Two configurations are considered, one in which a short-wave limit approximation is used, and another in which a long-wave limit approximation is used. In the short-wave limit, Wentzel-Kramers-Brillouin (WKB) methods are utilized to estimate the eigenvalues, and the eigenfunctions are approximated in terms of Green’s functions. The procedure consists of transforming the Orr-Sommerfeld equation into a system of two second order ordinary differential equations for which the eigenvalues and the eigenfunctions can be approximated. In the long-wave limit approximation, solutions are expressed in terms of generalized hypergeometric functions. Our procedure works regardless of the values of the Reynolds number.展开更多
This paper generalizes the basic principle of multiplier-enlargement approach to approximating any nonbounded continuous functions with positive linear operators, and as an example, Bernstein polynomial operators are ...This paper generalizes the basic principle of multiplier-enlargement approach to approximating any nonbounded continuous functions with positive linear operators, and as an example, Bernstein polynomial operators are analysed and studied. This paper gives a certain theorem as a general rule to approximate any nonbounded continuous functions.展开更多
A least-squares mixed finite element (LSMFE) method for the numerical solution of fourth order parabolic problems analyzed and developed in this paper. The Ciarlet-Raviart mixed finite element space is used to approxi...A least-squares mixed finite element (LSMFE) method for the numerical solution of fourth order parabolic problems analyzed and developed in this paper. The Ciarlet-Raviart mixed finite element space is used to approximate. The a posteriori error estimator which is needed in the adaptive refinement algorithm is proposed. The local evaluation of the least-squares functional serves as a posteriori error estimator. The posteriori errors are effectively estimated. The convergence of the adaptive least-squares mixed finite element method is proved.展开更多
In order to deal with the issue of huge computational cost very well in direct numerical simulation, the traditional response surface method (RSM) as a classical regression algorithm is used to approximate a functiona...In order to deal with the issue of huge computational cost very well in direct numerical simulation, the traditional response surface method (RSM) as a classical regression algorithm is used to approximate a functional relationship between the state variable and basic variables in reliability design. The algorithm has treated successfully some problems of implicit performance function in reliability analysis. However, its theoretical basis of empirical risk minimization narrows its range of applications for...展开更多
General neural network inverse adaptive controller has two flaws: the first is the slow convergence speed; the second is the invalidation to the non-minimum phase system. These defects limit the scope in which the neu...General neural network inverse adaptive controller has two flaws: the first is the slow convergence speed; the second is the invalidation to the non-minimum phase system. These defects limit the scope in which the neural network inverse adaptive controller is used. We employ Davidon least squares in training the multi-layer feedforward neural network used in approximating the inverse model of plant to expedite the convergence, and then through constructing the pseudo-plant, a neural network inverse adaptive controller is put forward which is still effective to the nonlinear non-minimum phase system. The simulation results show the validity of this scheme.展开更多
An M-metric least square method for polynomial analogy is presented. The relative normal eqUation is of diagonal form, such that the concise solution formula is explicit, and it is suitable to Parallel computation. On...An M-metric least square method for polynomial analogy is presented. The relative normal eqUation is of diagonal form, such that the concise solution formula is explicit, and it is suitable to Parallel computation. On the other hand, by error analysis of a typical example, we can see that the presented method is reliable.展开更多
Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid...Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid.Moreover,the scaling between these two grid sizes is super-linear.Approximation,stability and convergence aspects of a fully discrete scheme are analyzed.At last a numrical example is given whose results show that the algorithm proposed in this paper is effcient.展开更多
In this paper, the normal approximation rate and the random weighting approximation rate of error distribution of the kernel estimator of conditional density function f(y|x) are studied. The results may be used to...In this paper, the normal approximation rate and the random weighting approximation rate of error distribution of the kernel estimator of conditional density function f(y|x) are studied. The results may be used to construct the confidence interval of f(y|x) .展开更多
In this paper, the method proposed recently by the author for the solution of probability density function (PDF) of nonlinear stochastic systems is presented in detail and extended for more general problems of stochas...In this paper, the method proposed recently by the author for the solution of probability density function (PDF) of nonlinear stochastic systems is presented in detail and extended for more general problems of stochastic differential equations (SDE), therefore the Fokker Planck Kolmogorov (FPK) equation is expressed in general form with no limitation on the degree of nonlinearity of the SDE, the type of δ correlated excitations, the existence of multiplicative excitations, and the dimension of SDE or FPK equation. Examples are given and numerical results are provided for comparing with known exact solution to show the effectiveness of the method.展开更多
The multilevel characteristic basis function method(MLCBFM)with the adaptive cross approximation(ACA)algorithm for accelerated solution of electrically large scattering problems is studied in this paper.In the convent...The multilevel characteristic basis function method(MLCBFM)with the adaptive cross approximation(ACA)algorithm for accelerated solution of electrically large scattering problems is studied in this paper.In the conventional MLCBFM based on Foldy-Lax multiple scattering equations,the improvement is only made in the generation of characteristic basis functions(CBFs).However,it does not provide a change in impedance matrix filling and reducing matrix calculation procedure,which is time-consuming.In reality,all the impedance and reduced matrix of each level of the MLCBFM have low-rank property and can be calculated efficiently.Therefore,ACA is used for the efficient generation of two-level CBFs and the fast calculation of reduced matrix in this study.Numerical results are given to demonstrate the accuracy and efficiency of the method.展开更多
In this paper, we apply the two-time Green's function method, and provide a simple way to study themagnetic properties of one-dimensional spin-(S, s) Heisenberg ferromagnets.The magnetic susceptibility and correla...In this paper, we apply the two-time Green's function method, and provide a simple way to study themagnetic properties of one-dimensional spin-(S, s) Heisenberg ferromagnets.The magnetic susceptibility and correlationfunctions are obtained by using the Tyablikov decoupling approximation.Our results show that the magnetic susceptibilityand correlation length are a monotonically decreasing function of temperature regardless of the mixed spins.It isfound that in the case of S = s, our results of one-dimensional mixed-spin model is reduced to be those of the isotropicferromagnetic Heisenberg chain in the whole temperature region.Our results for the susceptibility are in agreement withthose obtained by other theoretical approaches.展开更多
In this note, we experimentally demonstrate, on a variety of analytic and nonanalytic functions, the novel observation that if the least squares polynomial approximation is repeated as weight in a second, now weighted...In this note, we experimentally demonstrate, on a variety of analytic and nonanalytic functions, the novel observation that if the least squares polynomial approximation is repeated as weight in a second, now weighted, least squares approximation, then this new, second, approximation is nearly perfect in the uniform sense, barely needing any further, say, Remez correction.展开更多
文摘The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.
文摘This paper presents a new approach to the structural topology optimization of continuum structures. Material-point independent variables are presented to illustrate the existence condition,or inexistence of the material points and their vicinity instead of elements or nodes in popular topology optimization methods. Topological variables field is constructed by moving least square approximation which is used as a shape function in the meshless method. Combined with finite element analyses,not only checkerboard patterns and mesh-dependence phenomena are overcome by this continuous and smooth topological variables field,but also the locations and numbers of topological variables can be arbitrary. Parameters including the number of quadrature points,scaling parameter,weight function and so on upon optimum topological configurations are discussed. Two classic topology optimization problems are solved successfully by the proposed method. The method is found robust and no numerical instabilities are found with proper parameters.
基金supported by the National Natural Science Foundation of China(Grant No.11101454)the Natural Science Foundation of Chongqing CSTC,China(Grant No.cstc2014jcyjA00005)the Program of Innovation Team Project in University of Chongqing City,China(Grant No.KJTD201308)
文摘Based on the moving least square （MLS） approximations and the boundary integral equations （BIEs）, a meshless algorithm is presented in this paper for elliptic Signorini problems. In the algorithm, a projection operator is used to tackle the nonlinear boundary inequality conditions. The Signorini problem is then reformulated as BIEs and the unknown boundary variables are approximated by the MLS approximations. Accordingly, only a nodal data structure on the boundary of a domain is required. The convergence of the algorithm is proven. Numerical examples are given to show the high convergence rate and high computational efficiency of the presented algorithm.
文摘By combining the classical appropriate functions “1, x, x 2” with the method of multiplier enlargement, this paper establishes a theorem to approximate any unbounded continuous functions with modified positive linear operators. As an example, Hermite Fejér interpolation polynomial operators are analysed and studied, and a general conclusion is obtained.
文摘In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem. the evaluation of the objective function is very difficult, so that only their approximate values can be obtained. This algorithm is obtained by combining penalty function method and approximation in bilevel programming. The presented algorithm is completely different from existing methods. That convergence for this algorithm is proved.
文摘This paper is devoted to the study of approximation of the solution for the differential equation whose coefficients are almost period functions. To this end the authors establish the estimation of the solution of general linear differential equation for infinite interval case. For finite interval case, this equation was investigated by G. Tamarkin([1]) applying the Picard method of successive approximation.
基金Sponsored by the Ministerial Level Advanced Research Foundation (010896367)
文摘The independent continuous mapping（ICM） method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of the essential boundary condition and derivative of various sensitivities,a singular weight function in element free Galerkin method is introduced.Material point variable is defined to illustrate the condition of material point and its vicinity instead of element or node.The topological variables field is constructed by moving least square approximation which inherits the continuity and smoothness of the weight function.Due to reciprocal relationships between the topological variables and design variables,various structural responses sensitivities are derived according to the method for calculating the partial derivatives of compound functions.Numerical examples indicate that checkerboard pattern and mesh-dependence phenomena are overcome without additional restriction methods.
基金Supported by the National Natural Science Foundation of China(61672009,61502130).
文摘This paper presents an interpolation-based method(IBM)for approximating some trigonometric functions or their integrals as well.It provides two-sided bounds for each function,which also achieves much better approximation effects than those of prevailing methods.In principle,the IBM can be applied for bounding more bounded smooth functions and their integrals as well,and its applications include approximating the integral of sin(x)/x function and improving the famous square root inequalities.
文摘Asymptotic eigenvalues and eigenfunctions for the Orr-Sommerfeld equation in two-dimensional and three-dimensional incompressible flows on an infinite domain and on a semi-infinite domain are obtained. Two configurations are considered, one in which a short-wave limit approximation is used, and another in which a long-wave limit approximation is used. In the short-wave limit, Wentzel-Kramers-Brillouin (WKB) methods are utilized to estimate the eigenvalues, and the eigenfunctions are approximated in terms of Green’s functions. The procedure consists of transforming the Orr-Sommerfeld equation into a system of two second order ordinary differential equations for which the eigenvalues and the eigenfunctions can be approximated. In the long-wave limit approximation, solutions are expressed in terms of generalized hypergeometric functions. Our procedure works regardless of the values of the Reynolds number.
文摘This paper generalizes the basic principle of multiplier-enlargement approach to approximating any nonbounded continuous functions with positive linear operators, and as an example, Bernstein polynomial operators are analysed and studied. This paper gives a certain theorem as a general rule to approximate any nonbounded continuous functions.
文摘A least-squares mixed finite element (LSMFE) method for the numerical solution of fourth order parabolic problems analyzed and developed in this paper. The Ciarlet-Raviart mixed finite element space is used to approximate. The a posteriori error estimator which is needed in the adaptive refinement algorithm is proposed. The local evaluation of the least-squares functional serves as a posteriori error estimator. The posteriori errors are effectively estimated. The convergence of the adaptive least-squares mixed finite element method is proved.
基金National High-tech Research and Development Pro-gram (2006AA04Z405)
文摘In order to deal with the issue of huge computational cost very well in direct numerical simulation, the traditional response surface method (RSM) as a classical regression algorithm is used to approximate a functional relationship between the state variable and basic variables in reliability design. The algorithm has treated successfully some problems of implicit performance function in reliability analysis. However, its theoretical basis of empirical risk minimization narrows its range of applications for...
基金Tianjin Natural Science Foundation !983602011National 863/CIMS Research Foundation !863-511-945-010
文摘General neural network inverse adaptive controller has two flaws: the first is the slow convergence speed; the second is the invalidation to the non-minimum phase system. These defects limit the scope in which the neural network inverse adaptive controller is used. We employ Davidon least squares in training the multi-layer feedforward neural network used in approximating the inverse model of plant to expedite the convergence, and then through constructing the pseudo-plant, a neural network inverse adaptive controller is put forward which is still effective to the nonlinear non-minimum phase system. The simulation results show the validity of this scheme.
文摘An M-metric least square method for polynomial analogy is presented. The relative normal eqUation is of diagonal form, such that the concise solution formula is explicit, and it is suitable to Parallel computation. On the other hand, by error analysis of a typical example, we can see that the presented method is reliable.
文摘Two-level finite element approximation to stream function form of unsteady Navier-Stokes equations is studied.This algorithm involves solving one nonlinear system on a coarse grid and one linear problem on a fine grid.Moreover,the scaling between these two grid sizes is super-linear.Approximation,stability and convergence aspects of a fully discrete scheme are analyzed.At last a numrical example is given whose results show that the algorithm proposed in this paper is effcient.
基金Supported by Natural Science Foundation of Beijing City and National Natural Science Foundation ofChina(2 2 30 4 1 0 0 1 30 1
文摘In this paper, the normal approximation rate and the random weighting approximation rate of error distribution of the kernel estimator of conditional density function f(y|x) are studied. The results may be used to construct the confidence interval of f(y|x) .
文摘In this paper, the method proposed recently by the author for the solution of probability density function (PDF) of nonlinear stochastic systems is presented in detail and extended for more general problems of stochastic differential equations (SDE), therefore the Fokker Planck Kolmogorov (FPK) equation is expressed in general form with no limitation on the degree of nonlinearity of the SDE, the type of δ correlated excitations, the existence of multiplicative excitations, and the dimension of SDE or FPK equation. Examples are given and numerical results are provided for comparing with known exact solution to show the effectiveness of the method.
基金supported by the National Natural Science Foundation of China (No.61401003)the Specialized Research Fund for the Doctoral Program of Higher Education of China (No.20123401110006)the Natural Science Research Project of Anhui Education ( No. KJ2015A436)
文摘The multilevel characteristic basis function method(MLCBFM)with the adaptive cross approximation(ACA)algorithm for accelerated solution of electrically large scattering problems is studied in this paper.In the conventional MLCBFM based on Foldy-Lax multiple scattering equations,the improvement is only made in the generation of characteristic basis functions(CBFs).However,it does not provide a change in impedance matrix filling and reducing matrix calculation procedure,which is time-consuming.In reality,all the impedance and reduced matrix of each level of the MLCBFM have low-rank property and can be calculated efficiently.Therefore,ACA is used for the efficient generation of two-level CBFs and the fast calculation of reduced matrix in this study.Numerical results are given to demonstrate the accuracy and efficiency of the method.
基金Supported by the Natural Science Foundation of Guangdong Province under Grant No.8151009001000055
文摘In this paper, we apply the two-time Green's function method, and provide a simple way to study themagnetic properties of one-dimensional spin-(S, s) Heisenberg ferromagnets.The magnetic susceptibility and correlationfunctions are obtained by using the Tyablikov decoupling approximation.Our results show that the magnetic susceptibilityand correlation length are a monotonically decreasing function of temperature regardless of the mixed spins.It isfound that in the case of S = s, our results of one-dimensional mixed-spin model is reduced to be those of the isotropicferromagnetic Heisenberg chain in the whole temperature region.Our results for the susceptibility are in agreement withthose obtained by other theoretical approaches.
文摘In this note, we experimentally demonstrate, on a variety of analytic and nonanalytic functions, the novel observation that if the least squares polynomial approximation is repeated as weight in a second, now weighted, least squares approximation, then this new, second, approximation is nearly perfect in the uniform sense, barely needing any further, say, Remez correction.
