In this paper, we made a new breakthrough, which proposes a new recursion–transform(RT) method with potential parameters to evaluate the nodal potential in arbitrary resistor networks. For the first time, we found ...In this paper, we made a new breakthrough, which proposes a new recursion–transform(RT) method with potential parameters to evaluate the nodal potential in arbitrary resistor networks. For the first time, we found the exact potential formulae of arbitrary m × n cobweb and fan networks by the RT method, and the potential formulae of infinite and semi-infinite networks are derived. As applications, a series of interesting corollaries of potential formulae are given by using the general formula, the equivalent resistance formula is deduced by using the potential formula, and we find a new trigonometric identity by comparing two equivalence results with different forms.展开更多
We gave the localized solutions,the interaction solutions and the mixed solutions to a reduced(3+1)-dimensional nonlinear evolution equation.These solutions were characterized by superposition formulas of positive qua...We gave the localized solutions,the interaction solutions and the mixed solutions to a reduced(3+1)-dimensional nonlinear evolution equation.These solutions were characterized by superposition formulas of positive quadratic functions,the exponential and hyperbolic functions.According to the known lump solution in the outset,we obtained the superposition formulas of positive quadratic functions by plausible reasoning.Next,we constructed the interaction solutions between the localized solutions and the exponential function solutions with the similar theory.These two kinds of solutions contained superposition formulas of positive quadratic functions,which were turned into general ternary quadratic functions,the coefficients of which were all rational operation of vector inner product.Then we obtained linear superposition formulas of exponential and hyperbolic function solutions.Finally,for aforementioned various solutions,their dynamic properties were showed by choosing specific values for parameters.From concrete plots,we observed wave characteristics of three kinds of solutions.Especially,we could observe distinct generation and separation situations when the localized wave and the stripe wave interacted at different time points.展开更多
Based on the nonlocal theory and Mindlin plate theory,the governing equations(i.e.,a system of partial differential equations(PDEs)for bending problem)of magnetoelectroelastic(MEE)nanoplates resting on the Pasternak e...Based on the nonlocal theory and Mindlin plate theory,the governing equations(i.e.,a system of partial differential equations(PDEs)for bending problem)of magnetoelectroelastic(MEE)nanoplates resting on the Pasternak elastic foundation are first derived by the variational principle.The polynomial particular solutions corresponding to the established model are then obtained and further employed as basis functions with the method of particular solutions(MPS)to solve the governing equations numerically.It is confirmed that for the present bending model,the new solution strategy possesses more general applicability and superior flexibility in the selection of collocation points.The effects of different boundary conditions,applied loads,and geometrical shapes on the bending properties of MEE nanoplates are evaluated by using the developed method.Some important conclusions are drawn,which should be helpful for the design and applications of electromagnetic nanoplate structures.展开更多
The boundary element method is used for he modal analysis of freevibration of 2-D composite structures in this paper. Since theparticular solution method is used to treat the terms of body forces(inertial forces) in t...The boundary element method is used for he modal analysis of freevibration of 2-D composite structures in this paper. Since theparticular solution method is used to treat the terms of body forces(inertial forces) in the equation of motion, only static fundamentalsolutions are needed in solving the problem. For an isotropiccantilever beam, the numerical results obtained by using the BEMpresented in this paper are in good agreement, with those of usingFEM or other BEM, but this BEM can also be used to analyze problemsfor anisotropic materials.展开更多
The aim of this paper is to propose a fast meshless numerical scheme for the simulation of non-linear Schrodinger equations.In the proposed scheme,the implicit-Euler scheme is used for the temporal discretization and ...The aim of this paper is to propose a fast meshless numerical scheme for the simulation of non-linear Schrodinger equations.In the proposed scheme,the implicit-Euler scheme is used for the temporal discretization and the localized method of approximate particular solution(LMAPS)is utilized for the spatial discretization.The multiple-scale technique is introduced to obtain the shape parameters of the multiquadric radial basis function for 2D problems and the Gaussian radial basis function for 3D problems.Six numerical examples are carried out to verify the accuracy and efficiency of the proposed scheme.Compared with well-known techniques,numerical results illustrate that the proposed scheme is of merits being easy-to-program,high accuracy,and rapid convergence even for long-term problems.These results also indicate that the proposed scheme has great potential in large scale problems and real-world applications.展开更多
In this paper we develop an efficient meshless method for solving inhomogeneous elliptic partial differential equations.We first approximate the source function by Chebyshev polynomials.We then focus on how to find a ...In this paper we develop an efficient meshless method for solving inhomogeneous elliptic partial differential equations.We first approximate the source function by Chebyshev polynomials.We then focus on how to find a polynomial particular solution when the source function is a polynomial.Through the principle of the method of undetermined coefficients and a proper arrangement of the terms for the polynomial particular solution to be determined,the coefficients of the particular solution satisfy a triangular system of linear algebraic equations.Explicit recursive formulas for the coefficients of the particular solutions are derived for different types of elliptic PDEs.The method is further incorporated into the method of fundamental solutions for solving inhomogeneous elliptic PDEs.Numerical results show that our approach is efficient and accurate.展开更多
The boundary element method for the modal analysis of freevibration for 3-D anisotropic structures using particular solutionshas been developed. The complete polynomials of order two are used toconstruct the particula...The boundary element method for the modal analysis of freevibration for 3-D anisotropic structures using particular solutionshas been developed. The complete polynomials of order two are used toconstruct the particular solutions for general anisotropic materials.The numerical results for 3-D free vibra- tion analysis of anisotropic cantilever beam by the method presented is in goodagreement with the results us- ing the Ritz technique. Foranisotropic materials, the numerical results calculated form theproposed method are in good agreement with the results from MSC.NASTRAN.展开更多
On the basis of von Karnmnequations, the axisymmetric buckling and post-buckling of annular plates on anelastic foundation is so wematically discussed byusing shooting
The stress rate integral equations of elastoplasticity are deduced based on Ref. [1] by consistent methods. The point at which the stresses and/or displacements are calculated can be in the body or on the boundary, an...The stress rate integral equations of elastoplasticity are deduced based on Ref. [1] by consistent methods. The point at which the stresses and/or displacements are calculated can be in the body or on the boundary, and in the plastic region or elastic one. The existence of the principal value integral in the plastic region is demonstrated strictly, and the theoretical basis is presented for the paticular solution method by unit initial stress fields. In the present method, programming is easy and general, and the numerical results are excellent.展开更多
In this study,the polynomial expansion method(PEM)and the polynomial method of particular solutions(PMPS)are applied to solve a class of linear elliptic partial differential equations(PDEs)in two dimensions with const...In this study,the polynomial expansion method(PEM)and the polynomial method of particular solutions(PMPS)are applied to solve a class of linear elliptic partial differential equations(PDEs)in two dimensions with constant coefficients.In the solution procedure,the sought solution is approximated by the Pascal polynomials and their particular solutions for the PEM and PMPS,respectively.The multiple-scale technique is applied to improve the conditioning of the resulted linear equations and the accuracy of numerical results for both of the PEM and PMPS.Some mathematical statements are provided to demonstrate the equivalence of the PEM and PMPS bases as they are both bases of a certain polynomial vector space.Then,some numerical experiments were conducted to validate the implementation of the PEM and PMPS.Numerical results demonstrated that the PEM is more accurate and well-conditioned than the PMPS and the multiple-scale technique is essential in these polynomial methods.展开更多
A meshless method based on the method of fundamental solutions(MFS)is proposed to solve the time-dependent partial differential equations with variable coefficients.The proposed method combines the time discretization...A meshless method based on the method of fundamental solutions(MFS)is proposed to solve the time-dependent partial differential equations with variable coefficients.The proposed method combines the time discretization and the onestage MFS for spatial discretization.In contrast to the traditional two-stage process,the one-stage MFS approach is capable of solving a broad spectrum of partial differential equations.The numerical implementation is simple since both closed-form approximate particular solution and fundamental solution are easy to find than the traditional approach.The numerical results show that the one-stage approach is robust and stable.展开更多
In this paper,we applied the polyharmonic splines as the basis functions to derive particular solutions for the differential operator ∆^(2) ± λ^(2).Similar to the derivation of fundamental solutions,it is non-tr...In this paper,we applied the polyharmonic splines as the basis functions to derive particular solutions for the differential operator ∆^(2) ± λ^(2).Similar to the derivation of fundamental solutions,it is non-trivial to derive particular solutions for higher order differential operators.In this paper,we provide a simple algebraic factorization approach to derive particular solutions for these types of differential operators in 2D and 3D.The main focus of this paper is its simplicity in the sense that minimal mathematical background is required for numerically solving higher order partial differential equations such as thin plate vibration.Three numerical examples in both 2D and 3D are given to validate particular solutions we derived.展开更多
Choosing particular solution source and its position have great influence on accu- racy of sound field prediction in distributed source boundary point method. An optimization method for determining the position of par...Choosing particular solution source and its position have great influence on accu- racy of sound field prediction in distributed source boundary point method. An optimization method for determining the position of particular solution sources is proposed to get high accu- racy prediction result. In this method, tripole is chosen as the particular solution. The upper limit frequency of calculation is predicted by setting 1% volume velocity relative error limit using vibration velocity of structure surface. Then, the optimal position of particular solution sources, in which the relative error of volume velocity is minimum, is determined within the range of upper limit frequency by searching algorithm using volume velocity matching. The transfer matrix between pressure and surface volume velocity is constructed in the optimal position. After that, the sound radiation of structure is calculated by the matrix. The results of numerical simulation show that the calculation error is significantly reduced by the proposed method. When there are vibration velocity measurement errors, the calculation errors can be controlled within 5% by the method.展开更多
Using the complex variable method,an elastic analytical solution of the ground displacement caused by a shallow circular tunneling is derived.Non-symmetric deformation relative to the horizontal center line of the tun...Using the complex variable method,an elastic analytical solution of the ground displacement caused by a shallow circular tunneling is derived.Non-symmetric deformation relative to the horizontal center line of the tunnel cross-section is used as a boundary condition.A comparison between the proposed analytical method and the Finite Element Method is carried out to validate the rationality of the obtained analytical solution.Two parameters in the Peck formula,namely the maximum settlement of the ground surface center and the width coefficient of settlement curve,are fitted and determined.We propose a modified Peck formula by considering three input parameters,namely the tunnel depth,tunnel radius,and the tunnel gap.The influence of these three parameters on the modified Peck formula is analyzed.The applicability of the modified Peck formula is further investigated by reference to six engineering projects.The ground surface displacement obtained by the explicit Peck formula is in good agreement with the field data,and the maximum error is only 1.3 cm.The proposed formula can quickly and reasonably predict the ground surface settlement caused by tunnelling.展开更多
In this paper,we propose efficient algorithms for approximating particular solutions of second and fourth order elliptic equations.The approximation of the particular solution by a truncated series of Chebyshev polyno...In this paper,we propose efficient algorithms for approximating particular solutions of second and fourth order elliptic equations.The approximation of the particular solution by a truncated series of Chebyshev polynomials and the satisfaction of the differential equation lead to upper triangular block systems,each block being an upper triangular system.These systems can be solved efficiently by standard techniques.Several numerical examples are presented for each case.展开更多
基金Project supported by the Natural Science Foundation of Jiangsu Province,China(Grant No.BK20161278)
文摘In this paper, we made a new breakthrough, which proposes a new recursion–transform(RT) method with potential parameters to evaluate the nodal potential in arbitrary resistor networks. For the first time, we found the exact potential formulae of arbitrary m × n cobweb and fan networks by the RT method, and the potential formulae of infinite and semi-infinite networks are derived. As applications, a series of interesting corollaries of potential formulae are given by using the general formula, the equivalent resistance formula is deduced by using the potential formula, and we find a new trigonometric identity by comparing two equivalence results with different forms.
基金the National Natural Science Foundation of China(Grant No.12061054)Program for Young Talents of Science and Technology in Universities of Inner Mongolia Autonomous Region of China(Grant No.NJYT-20A06)。
文摘We gave the localized solutions,the interaction solutions and the mixed solutions to a reduced(3+1)-dimensional nonlinear evolution equation.These solutions were characterized by superposition formulas of positive quadratic functions,the exponential and hyperbolic functions.According to the known lump solution in the outset,we obtained the superposition formulas of positive quadratic functions by plausible reasoning.Next,we constructed the interaction solutions between the localized solutions and the exponential function solutions with the similar theory.These two kinds of solutions contained superposition formulas of positive quadratic functions,which were turned into general ternary quadratic functions,the coefficients of which were all rational operation of vector inner product.Then we obtained linear superposition formulas of exponential and hyperbolic function solutions.Finally,for aforementioned various solutions,their dynamic properties were showed by choosing specific values for parameters.From concrete plots,we observed wave characteristics of three kinds of solutions.Especially,we could observe distinct generation and separation situations when the localized wave and the stripe wave interacted at different time points.
基金Project supported by the National Natural Science Foundation of China(Nos.11872257 and 11572358)the German Research Foundation(No.ZH 15/14-1)。
文摘Based on the nonlocal theory and Mindlin plate theory,the governing equations(i.e.,a system of partial differential equations(PDEs)for bending problem)of magnetoelectroelastic(MEE)nanoplates resting on the Pasternak elastic foundation are first derived by the variational principle.The polynomial particular solutions corresponding to the established model are then obtained and further employed as basis functions with the method of particular solutions(MPS)to solve the governing equations numerically.It is confirmed that for the present bending model,the new solution strategy possesses more general applicability and superior flexibility in the selection of collocation points.The effects of different boundary conditions,applied loads,and geometrical shapes on the bending properties of MEE nanoplates are evaluated by using the developed method.Some important conclusions are drawn,which should be helpful for the design and applications of electromagnetic nanoplate structures.
文摘The boundary element method is used for he modal analysis of freevibration of 2-D composite structures in this paper. Since theparticular solution method is used to treat the terms of body forces(inertial forces) in the equation of motion, only static fundamentalsolutions are needed in solving the problem. For an isotropiccantilever beam, the numerical results obtained by using the BEMpresented in this paper are in good agreement, with those of usingFEM or other BEM, but this BEM can also be used to analyze problemsfor anisotropic materials.
基金The authors thank the editor and anonymous reviewers for their constructive comments on the manuscript.The research of the authors was supported by the Natural Science Foundation of Jiangsu Province(No.BK20150795)the Fundamental Research Funds for the Central Universities(No.2018B16714)+3 种基金the National Natural Science Foundation of China(Nos.11702083,51679150,51579153,51739008,51527811)the State Key Laboratory of Mechanics and Control of Mechanical Structures(Nanjing University of Aeronautics and Astronautics)(No.MCMS-0218G01)the National Key R&D Program of China(No.2016YFC0401902)the Fund Project of NHRI(Nos.Y417002,Y417015).
文摘The aim of this paper is to propose a fast meshless numerical scheme for the simulation of non-linear Schrodinger equations.In the proposed scheme,the implicit-Euler scheme is used for the temporal discretization and the localized method of approximate particular solution(LMAPS)is utilized for the spatial discretization.The multiple-scale technique is introduced to obtain the shape parameters of the multiquadric radial basis function for 2D problems and the Gaussian radial basis function for 3D problems.Six numerical examples are carried out to verify the accuracy and efficiency of the proposed scheme.Compared with well-known techniques,numerical results illustrate that the proposed scheme is of merits being easy-to-program,high accuracy,and rapid convergence even for long-term problems.These results also indicate that the proposed scheme has great potential in large scale problems and real-world applications.
文摘In this paper we develop an efficient meshless method for solving inhomogeneous elliptic partial differential equations.We first approximate the source function by Chebyshev polynomials.We then focus on how to find a polynomial particular solution when the source function is a polynomial.Through the principle of the method of undetermined coefficients and a proper arrangement of the terms for the polynomial particular solution to be determined,the coefficients of the particular solution satisfy a triangular system of linear algebraic equations.Explicit recursive formulas for the coefficients of the particular solutions are derived for different types of elliptic PDEs.The method is further incorporated into the method of fundamental solutions for solving inhomogeneous elliptic PDEs.Numerical results show that our approach is efficient and accurate.
基金the National Natural Science Foundation of China (No.19872011)
文摘The boundary element method for the modal analysis of freevibration for 3-D anisotropic structures using particular solutionshas been developed. The complete polynomials of order two are used toconstruct the particular solutions for general anisotropic materials.The numerical results for 3-D free vibra- tion analysis of anisotropic cantilever beam by the method presented is in goodagreement with the results us- ing the Ritz technique. Foranisotropic materials, the numerical results calculated form theproposed method are in good agreement with the results from MSC.NASTRAN.
基金Project support by the State Education Commission of the People’s Republic of China
文摘On the basis of von Karnmnequations, the axisymmetric buckling and post-buckling of annular plates on anelastic foundation is so wematically discussed byusing shooting
基金The project supported by the National Natural Science Foundation of China
文摘The stress rate integral equations of elastoplasticity are deduced based on Ref. [1] by consistent methods. The point at which the stresses and/or displacements are calculated can be in the body or on the boundary, and in the plastic region or elastic one. The existence of the principal value integral in the plastic region is demonstrated strictly, and the theoretical basis is presented for the paticular solution method by unit initial stress fields. In the present method, programming is easy and general, and the numerical results are excellent.
基金The Ministry of Science and Technology of Taiwan is gratefully acknowledged for providing financial support to carry out the present work under the Grant No.MOST 109-2221-E-992-046-MY3.
文摘In this study,the polynomial expansion method(PEM)and the polynomial method of particular solutions(PMPS)are applied to solve a class of linear elliptic partial differential equations(PDEs)in two dimensions with constant coefficients.In the solution procedure,the sought solution is approximated by the Pascal polynomials and their particular solutions for the PEM and PMPS,respectively.The multiple-scale technique is applied to improve the conditioning of the resulted linear equations and the accuracy of numerical results for both of the PEM and PMPS.Some mathematical statements are provided to demonstrate the equivalence of the PEM and PMPS bases as they are both bases of a certain polynomial vector space.Then,some numerical experiments were conducted to validate the implementation of the PEM and PMPS.Numerical results demonstrated that the PEM is more accurate and well-conditioned than the PMPS and the multiple-scale technique is essential in these polynomial methods.
基金the support of NATO Collaborative Linkage Grant under reference ESP.CLG.982891。
文摘A meshless method based on the method of fundamental solutions(MFS)is proposed to solve the time-dependent partial differential equations with variable coefficients.The proposed method combines the time discretization and the onestage MFS for spatial discretization.In contrast to the traditional two-stage process,the one-stage MFS approach is capable of solving a broad spectrum of partial differential equations.The numerical implementation is simple since both closed-form approximate particular solution and fundamental solution are easy to find than the traditional approach.The numerical results show that the one-stage approach is robust and stable.
文摘In this paper,we applied the polyharmonic splines as the basis functions to derive particular solutions for the differential operator ∆^(2) ± λ^(2).Similar to the derivation of fundamental solutions,it is non-trivial to derive particular solutions for higher order differential operators.In this paper,we provide a simple algebraic factorization approach to derive particular solutions for these types of differential operators in 2D and 3D.The main focus of this paper is its simplicity in the sense that minimal mathematical background is required for numerically solving higher order partial differential equations such as thin plate vibration.Three numerical examples in both 2D and 3D are given to validate particular solutions we derived.
文摘Choosing particular solution source and its position have great influence on accu- racy of sound field prediction in distributed source boundary point method. An optimization method for determining the position of particular solution sources is proposed to get high accu- racy prediction result. In this method, tripole is chosen as the particular solution. The upper limit frequency of calculation is predicted by setting 1% volume velocity relative error limit using vibration velocity of structure surface. Then, the optimal position of particular solution sources, in which the relative error of volume velocity is minimum, is determined within the range of upper limit frequency by searching algorithm using volume velocity matching. The transfer matrix between pressure and surface volume velocity is constructed in the optimal position. After that, the sound radiation of structure is calculated by the matrix. The results of numerical simulation show that the calculation error is significantly reduced by the proposed method. When there are vibration velocity measurement errors, the calculation errors can be controlled within 5% by the method.
基金National Natural Science Foundation of China(Grant Nos.52025084 and 52208396)China Postdoctoral Science Foundation(No.2021M701934).
文摘Using the complex variable method,an elastic analytical solution of the ground displacement caused by a shallow circular tunneling is derived.Non-symmetric deformation relative to the horizontal center line of the tunnel cross-section is used as a boundary condition.A comparison between the proposed analytical method and the Finite Element Method is carried out to validate the rationality of the obtained analytical solution.Two parameters in the Peck formula,namely the maximum settlement of the ground surface center and the width coefficient of settlement curve,are fitted and determined.We propose a modified Peck formula by considering three input parameters,namely the tunnel depth,tunnel radius,and the tunnel gap.The influence of these three parameters on the modified Peck formula is analyzed.The applicability of the modified Peck formula is further investigated by reference to six engineering projects.The ground surface displacement obtained by the explicit Peck formula is in good agreement with the field data,and the maximum error is only 1.3 cm.The proposed formula can quickly and reasonably predict the ground surface settlement caused by tunnelling.
文摘In this paper,we propose efficient algorithms for approximating particular solutions of second and fourth order elliptic equations.The approximation of the particular solution by a truncated series of Chebyshev polynomials and the satisfaction of the differential equation lead to upper triangular block systems,each block being an upper triangular system.These systems can be solved efficiently by standard techniques.Several numerical examples are presented for each case.