The current study examines the important class of Chebyshev’s differential equations via the application of the efficient Adomian Decomposition Method (ADM) and its modifications. We have proved the effectiveness of ...The current study examines the important class of Chebyshev’s differential equations via the application of the efficient Adomian Decomposition Method (ADM) and its modifications. We have proved the effectiveness of the employed methods by acquiring exact analytical solutions for the governing equations in most cases;while minimal noisy error terms have been observed in a particular method modification. Above all, the presented approaches have rightly affirmed the exactitude of the available literature. More to the point, the application of this methodology could be extended to examine various forms of high-order differential equations, as approximate exact solutions are rapidly attained with less computation stress.展开更多
A series solution of displacement response of the ground surface in the presence of underground twin tunnels subjected to excitation of incident plane SV waves is derived by using Fourier-Bessel series expansion metho...A series solution of displacement response of the ground surface in the presence of underground twin tunnels subjected to excitation of incident plane SV waves is derived by using Fourier-Bessel series expansion method.The numerical parametric study shows that underground twin tunnels significantly amplify the nearby surface ground motion.It is suggested that the effect of subways on ground motion should be considered when the subways are planned and designed.展开更多
A series solution for surface motion amplification due to underground group cavities for incident plane P waves is derived by Fourier-Bessel series expansion method. It is shown that underground group cavities signifi...A series solution for surface motion amplification due to underground group cavities for incident plane P waves is derived by Fourier-Bessel series expansion method. It is shown that underground group cavities significantly am-plify the surface ground motion nearby. It is suggested that the effect of subways on ground motion should be con-sidered when the subways are planned and designed.展开更多
From the point of view of approximate symmetry, the modified Korteweg-de Vries-Burgers (mKdV-Burgers) equation with weak dissipation is investigated. The symmetry of a system of the corresponding partial differentia...From the point of view of approximate symmetry, the modified Korteweg-de Vries-Burgers (mKdV-Burgers) equation with weak dissipation is investigated. The symmetry of a system of the corresponding partial differential equations which approximate the perturbed mKdV-Burgers equation is constructed and the corresponding general approximate symmetry reduction is derived; thereby infinite series solutions and general formulae can be obtained. The obtained result shows that the zero-order similarity solution to the mKdV-Burgers equation satisfies the Painleve II equation. Also, at the level of travelling wave reduction, the general solution formulae are given for any travelling wave solution of an unperturbed mKdV equation. As an illustrative example, when the zero-order tanh profile solution is chosen as an initial approximate solution, physically approximate similarity solutions are obtained recursively under the appropriate choice of parameters occurring during computation.展开更多
Starting from Lie symmetry theory and combining with the approximate symmetry method, and using the package LieSYMGRP proposed by us, we restudy the perturbed Kuramoto-Sivashinsky (KS) equation. The approximate symm...Starting from Lie symmetry theory and combining with the approximate symmetry method, and using the package LieSYMGRP proposed by us, we restudy the perturbed Kuramoto-Sivashinsky (KS) equation. The approximate symmetry reduction and the infinite series symmetry reduction solutions of the perturbed KS equation are constructed. Specially, if selecting the tanh-type travelling wave solution as initial approximate, we not only obtain the general formula of the physical approximate similarity solutions, but also obtain several new explicit solutions of the given equation, which are first reported here.展开更多
The multipole moment method not only conduces to the understanding of the deformation of the space-time, but also serves as an effective tool to approximately solve the Einstein field equation with. However, the usual...The multipole moment method not only conduces to the understanding of the deformation of the space-time, but also serves as an effective tool to approximately solve the Einstein field equation with. However, the usual multipole moments are recursively determined by a sequence of symmetric and trace-free tensors, which is inconvenient for practical resolution. In this paper, we develop a simplified procedure to generate the series solutions to the metric of the stationary vacuum with axisymmetry, and show its validity. In order to understand the free parameters in the solution, we propose to take the Schwarzschild metric as a standard ruler, and some well- known examples are analysed and compared with the series solutions in detail.展开更多
Nonlinear dynamic equation is a common engineering model.There is not precise analytical solution for most of nonlinear differential equations.These nonlinear differential equations should be solved by using approxima...Nonlinear dynamic equation is a common engineering model.There is not precise analytical solution for most of nonlinear differential equations.These nonlinear differential equations should be solved by using approximate methods.Classical perturbation methods such as LP method,KBM method,multi-scale method and the averaging method on weakly nonlinear vibration system is effective,while the strongly nonlinear system is difficult to apply.Approximate solutions of primary resonance for forced Duffing equation is investigated by means of homotopy analysis method (HAM).Different from other approximate computational method,the HAM is totally independent of small physical parameters,and thus is suitable for most nonlinear problems.The HAM provides a great freedom to choose base functions of solution series,so that a nonlinear problem may be approximated more effectively.The HAM provides us a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter and the auxiliary function.Therefore,HAM not only may solve the weakly non-linear problems but also may be suitable for the strong non-linear problem.Through the approximate solution of forced Duffing equation with cubic non-linearity,the HAM and fourth order Runge-Kutta method of numerical solution were compared,the results show that the HAM not only can solve the steady state solution,but also can calculate the unsteady state solution,and has the good computational accuracy.展开更多
The present manuscript examines the circular restricted gravitational three-body problem (CRGTBP) by the introduction of a new approach through the power series method. In addition, certain computational algorithms wi...The present manuscript examines the circular restricted gravitational three-body problem (CRGTBP) by the introduction of a new approach through the power series method. In addition, certain computational algorithms with the aid of Mathematica software are specifically designed for the problem. The algorithms or rather mathematical modules are established to determine the velocity and position of the third body’s motion. In fact, the modules led to accurate results and thus proved the new approach to be efficient.展开更多
The Schroedinger equation involving the phenomenon of the localization and entanglement for an exciton in a quantum dot molecule by an ac electric field is analytically investigated. New exact series solutions for the...The Schroedinger equation involving the phenomenon of the localization and entanglement for an exciton in a quantum dot molecule by an ac electric field is analytically investigated. New exact series solutions for the Schroedinger equation have been obtained for the first time. The analytical expressions can further describe the dynamical behaviors of an interacting electron-hole pair in a double coupled quantum dot molecule under an ac electric field accurately.展开更多
In this paper, the(2+1)-dimensional perturbed Boussinesq equation is transformed into a series of two-dimensional(2 D) similarity reduction equations by using the approximate symmetry method. A step-by-step proce...In this paper, the(2+1)-dimensional perturbed Boussinesq equation is transformed into a series of two-dimensional(2 D) similarity reduction equations by using the approximate symmetry method. A step-by-step procedure is used to acquire Jacobi elliptic function solutions to these similarity equations, which generate the truncated series solutions to the original perturbed Boussinesq equation. Aside from some singular area, the series solutions are convergent when the perturbation parameter is diminished.展开更多
In this article, we study on the existence of solution for a singularities of a system of nonlinear fractional differential equations (FDE). We construct a formal power series solution for our considering FDE and pr...In this article, we study on the existence of solution for a singularities of a system of nonlinear fractional differential equations (FDE). We construct a formal power series solution for our considering FDE and prove convergence of formal so- lutions under conditions. -We use the Caputo fractional differential operator and the nonlinearity depends on the fractional derivative of an unknown function.展开更多
Precision, Productivity and Performance def ine the GERBERcutter? Z7Tolland, Conn., USA – Gerber Technology, a business unit of Gerber Scientific, Inc. (NYSE: GRB), and the world leader in providing innovative integr...Precision, Productivity and Performance def ine the GERBERcutter? Z7Tolland, Conn., USA – Gerber Technology, a business unit of Gerber Scientific, Inc. (NYSE: GRB), and the world leader in providing innovative integrated software and hardware automation systems to展开更多
An analytical technique, namely the homotopy analysis method (HAM), is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small phys...An analytical technique, namely the homotopy analysis method (HAM), is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.展开更多
Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix d...Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.展开更多
The dynamic buckling of an elastic-plastic column subjected to axial impact by a rigid body has been discussed in this paper. The whole traveling process of elastic-plastic waves under impact action is analyzed with t...The dynamic buckling of an elastic-plastic column subjected to axial impact by a rigid body has been discussed in this paper. The whole traveling process of elastic-plastic waves under impact action is analyzed with the characteristics method. The regularity of stress changes in both column ends and the first separating time of a rigid body and column are obtained. By using the energy principle and taking into account the propagation and reflection of stress waves the lateral disturbance equation is derived and the power series solution is given. In addition, the critical buckling condition can be obtained from the stability analysis of the solution. By numerical computation and analysis, the relationship among critical velocity and impact mass, hardening modulus, and buckling time is given.展开更多
In this study, we proposed an analytical solution for eddy currents as well as electromagnetic forces of a conductive circular plate in a time varying magnetic field. Specifically, an analytical series solution for ed...In this study, we proposed an analytical solution for eddy currents as well as electromagnetic forces of a conductive circular plate in a time varying magnetic field. Specifically, an analytical series solution for eddy currents in a circular plate subjected to an axisymmetrie time varying magnetic field has been proposed based on the T-method that has been widely used in the eddy current analysis of conductive and superconductive structures. Accordingly, the dynamic response, the dynamic instability and the magnetic damping of a circular plate in a transverse transient magnetic field as well as a stationary in-plane magnetic field have also been obtained. The analytical series solution proposed in this work as well as the subsequent numerical analysis not only confirmed the emergence of dynamic instability of a circular plate in a strong transverse magnetic field, but also demonstrated the existence of magneto-damping of a circular conductive plate in an in-plane magnetic field. The method developed in this paper provides a potential new possible way by which the analysis of the electromagnetic coupling problems of conductive structures can be simplified.展开更多
Melting heat transfer in the boundary layer flow of a couple stress fluid over a stretching surface is investigated. The developed differential equations are solved for homotopic solutions. It is observed that the vel...Melting heat transfer in the boundary layer flow of a couple stress fluid over a stretching surface is investigated. The developed differential equations are solved for homotopic solutions. It is observed that the velocity and the boundary layer thickness are decreasing functions of the couple stress fluid parameter. However, the temperature and surface heat transfer increase when the values of the couple stress fluid parameter increase. The velocity and temperature fields increase with an increase in the melting process of the stretching sheet.展开更多
A homotopy analysis method(HAM)is presented for the primary resonance of multiple degree-of-freedom systems with strong non-linearity excited by harmonic forces.The validity of the HAM is independent of the existenc...A homotopy analysis method(HAM)is presented for the primary resonance of multiple degree-of-freedom systems with strong non-linearity excited by harmonic forces.The validity of the HAM is independent of the existence of small parameters in the considered equation.The HAM provides a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter.Two examples are presented to show that the HAM solutions agree well with the results of the modified Linstedt-Poincar'e method and the incremental harmonic balance method.展开更多
This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions an...This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions and similarity reduction equations of different orders, showing that the approximate homotopy direct method yields more general approximate similarity reductions than the approximate homotopy symmetry method. The homotopy series solutions to the generalized Kawahara equation are consequently derived.展开更多
Free and steady state forced transverse vibrations of non-uniform beams are investigated with a proposed method, leading to a series solution. The obtained series is verified to be convergent and linearly independent ...Free and steady state forced transverse vibrations of non-uniform beams are investigated with a proposed method, leading to a series solution. The obtained series is verified to be convergent and linearly independent in a convergence test and by the non-zero value of the corresponding Wronski determinant, respectively. The obtained solution is rigorous, which can be reduced to a classical solution for uniform beams. The proposed method can deal with arbitrary non-uniform Euler-Bernoulli beams in principle, but the methods in terms of special functions or elementary functions can only work in some special cases.展开更多
文摘The current study examines the important class of Chebyshev’s differential equations via the application of the efficient Adomian Decomposition Method (ADM) and its modifications. We have proved the effectiveness of the employed methods by acquiring exact analytical solutions for the governing equations in most cases;while minimal noisy error terms have been observed in a particular method modification. Above all, the presented approaches have rightly affirmed the exactitude of the available literature. More to the point, the application of this methodology could be extended to examine various forms of high-order differential equations, as approximate exact solutions are rapidly attained with less computation stress.
基金National Natural Science Foundation of China(50378063)EYTP of MOESRF for ROCS,MOE
文摘A series solution of displacement response of the ground surface in the presence of underground twin tunnels subjected to excitation of incident plane SV waves is derived by using Fourier-Bessel series expansion method.The numerical parametric study shows that underground twin tunnels significantly amplify the nearby surface ground motion.It is suggested that the effect of subways on ground motion should be considered when the subways are planned and designed.
基金Supported by National Natural Science Foundation of China (50378063), Excellent Young Teachers Program of MOE and SRF for ROCS, MOE.
文摘A series solution for surface motion amplification due to underground group cavities for incident plane P waves is derived by Fourier-Bessel series expansion method. It is shown that underground group cavities significantly am-plify the surface ground motion nearby. It is suggested that the effect of subways on ground motion should be con-sidered when the subways are planned and designed.
基金supported by the National Natural Science Foundations of China(Grant Nos 10735030,10475055,and 90503006)the National Basic Research Program of China(Grant No 2007CB814800)+1 种基金the Science Foundation for Post Doctorate Research from the Ministry of Science and Technology of China(Grant No 20070410727)the Natural Science Basic Research Plan in Shaanxi Province,China(Grant No SJ08A09)
文摘From the point of view of approximate symmetry, the modified Korteweg-de Vries-Burgers (mKdV-Burgers) equation with weak dissipation is investigated. The symmetry of a system of the corresponding partial differential equations which approximate the perturbed mKdV-Burgers equation is constructed and the corresponding general approximate symmetry reduction is derived; thereby infinite series solutions and general formulae can be obtained. The obtained result shows that the zero-order similarity solution to the mKdV-Burgers equation satisfies the Painleve II equation. Also, at the level of travelling wave reduction, the general solution formulae are given for any travelling wave solution of an unperturbed mKdV equation. As an illustrative example, when the zero-order tanh profile solution is chosen as an initial approximate solution, physically approximate similarity solutions are obtained recursively under the appropriate choice of parameters occurring during computation.
基金The project supported by National Natural Science Foundations of China under Grant Nos. 10735030, 10475055, and 90503006; the Natural Science Research Plan in Shaanxi Province under Grant No. SJ08A09; the Research Fund of Postdoctoral of China under Grant No. 20070410727;the Research Found of Shaanxi Normal University
文摘Starting from Lie symmetry theory and combining with the approximate symmetry method, and using the package LieSYMGRP proposed by us, we restudy the perturbed Kuramoto-Sivashinsky (KS) equation. The approximate symmetry reduction and the infinite series symmetry reduction solutions of the perturbed KS equation are constructed. Specially, if selecting the tanh-type travelling wave solution as initial approximate, we not only obtain the general formula of the physical approximate similarity solutions, but also obtain several new explicit solutions of the given equation, which are first reported here.
文摘The multipole moment method not only conduces to the understanding of the deformation of the space-time, but also serves as an effective tool to approximately solve the Einstein field equation with. However, the usual multipole moments are recursively determined by a sequence of symmetric and trace-free tensors, which is inconvenient for practical resolution. In this paper, we develop a simplified procedure to generate the series solutions to the metric of the stationary vacuum with axisymmetry, and show its validity. In order to understand the free parameters in the solution, we propose to take the Schwarzschild metric as a standard ruler, and some well- known examples are analysed and compared with the series solutions in detail.
基金supported by Fundamental Research Funds for the Central Universities of China (Grant No. N090405009)
文摘Nonlinear dynamic equation is a common engineering model.There is not precise analytical solution for most of nonlinear differential equations.These nonlinear differential equations should be solved by using approximate methods.Classical perturbation methods such as LP method,KBM method,multi-scale method and the averaging method on weakly nonlinear vibration system is effective,while the strongly nonlinear system is difficult to apply.Approximate solutions of primary resonance for forced Duffing equation is investigated by means of homotopy analysis method (HAM).Different from other approximate computational method,the HAM is totally independent of small physical parameters,and thus is suitable for most nonlinear problems.The HAM provides a great freedom to choose base functions of solution series,so that a nonlinear problem may be approximated more effectively.The HAM provides us a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter and the auxiliary function.Therefore,HAM not only may solve the weakly non-linear problems but also may be suitable for the strong non-linear problem.Through the approximate solution of forced Duffing equation with cubic non-linearity,the HAM and fourth order Runge-Kutta method of numerical solution were compared,the results show that the HAM not only can solve the steady state solution,but also can calculate the unsteady state solution,and has the good computational accuracy.
文摘The present manuscript examines the circular restricted gravitational three-body problem (CRGTBP) by the introduction of a new approach through the power series method. In addition, certain computational algorithms with the aid of Mathematica software are specifically designed for the problem. The algorithms or rather mathematical modules are established to determine the velocity and position of the third body’s motion. In fact, the modules led to accurate results and thus proved the new approach to be efficient.
基金The project partially supported by National Natural Science Foundation of China under Grant No.10247008 and the Science Foundation of Northwest Normal University of China under Grant No. NWNU-KJCXGC-02-04
文摘The Schroedinger equation involving the phenomenon of the localization and entanglement for an exciton in a quantum dot molecule by an ac electric field is analytically investigated. New exact series solutions for the Schroedinger equation have been obtained for the first time. The analytical expressions can further describe the dynamical behaviors of an interacting electron-hole pair in a double coupled quantum dot molecule under an ac electric field accurately.
基金Project supported by the National Natural Science Foundation of China(Grant No.11505094)the Natural Science Foundation of Jiangsu Province,China(Grant No.BK20150984)
文摘In this paper, the(2+1)-dimensional perturbed Boussinesq equation is transformed into a series of two-dimensional(2 D) similarity reduction equations by using the approximate symmetry method. A step-by-step procedure is used to acquire Jacobi elliptic function solutions to these similarity equations, which generate the truncated series solutions to the original perturbed Boussinesq equation. Aside from some singular area, the series solutions are convergent when the perturbation parameter is diminished.
文摘In this article, we study on the existence of solution for a singularities of a system of nonlinear fractional differential equations (FDE). We construct a formal power series solution for our considering FDE and prove convergence of formal so- lutions under conditions. -We use the Caputo fractional differential operator and the nonlinearity depends on the fractional derivative of an unknown function.
文摘Precision, Productivity and Performance def ine the GERBERcutter? Z7Tolland, Conn., USA – Gerber Technology, a business unit of Gerber Scientific, Inc. (NYSE: GRB), and the world leader in providing innovative integrated software and hardware automation systems to
文摘An analytical technique, namely the homotopy analysis method (HAM), is used to solve problems of nonlinear oscillations with parametric excitation. Unlike perturbation methods, HAM is not dependent on any small physical parameters at all, and thus valid for both weakly and strongly nonlinear problems. In addition, HAM is different from all other analytic techniques in providing a simple way to adjust and control convergence region of the series solution by means of an auxiliary parameter h. In the present paper, a periodic analytic approximations for nonlinear oscillations with parametric excitation are obtained by using HAM, and the results are validated by numerical simulations.
基金Project supported by the National Natural Science Foundation of China(Nos.60273048and60174023).
文摘Nonlinear dynamic equations can be solved accurately using a precise integration method. Some algorithms exist, but the inversion of a matrix must be calculated for these al- gorithms. If the inversion of the matrix doesn’t exist or isn’t stable, the precision and stability of the algorithms will be afected. An explicit series solution of the state equation has been pre- sented. The solution avoids calculating the inversion of a matrix and its precision can be easily controlled. In this paper, an implicit series solution of nonlinear dynamic equations is presented. The algorithm is more precise and stable than the explicit series solution and isn’t sensitive to the time-step. Finally, a numerical example is presented to demonstrate the efectiveness of the algorithm.
基金Project supported by the National Natural Science Foundation of China (No. 10472076).
文摘The dynamic buckling of an elastic-plastic column subjected to axial impact by a rigid body has been discussed in this paper. The whole traveling process of elastic-plastic waves under impact action is analyzed with the characteristics method. The regularity of stress changes in both column ends and the first separating time of a rigid body and column are obtained. By using the energy principle and taking into account the propagation and reflection of stress waves the lateral disturbance equation is derived and the power series solution is given. In addition, the critical buckling condition can be obtained from the stability analysis of the solution. By numerical computation and analysis, the relationship among critical velocity and impact mass, hardening modulus, and buckling time is given.
基金supported by the National Natural Science Foundation of China (No.10672070)the Program for New Century Excellent Talents in University (No. NCET-06-0896)
文摘In this study, we proposed an analytical solution for eddy currents as well as electromagnetic forces of a conductive circular plate in a time varying magnetic field. Specifically, an analytical series solution for eddy currents in a circular plate subjected to an axisymmetrie time varying magnetic field has been proposed based on the T-method that has been widely used in the eddy current analysis of conductive and superconductive structures. Accordingly, the dynamic response, the dynamic instability and the magnetic damping of a circular plate in a transverse transient magnetic field as well as a stationary in-plane magnetic field have also been obtained. The analytical series solution proposed in this work as well as the subsequent numerical analysis not only confirmed the emergence of dynamic instability of a circular plate in a strong transverse magnetic field, but also demonstrated the existence of magneto-damping of a circular conductive plate in an in-plane magnetic field. The method developed in this paper provides a potential new possible way by which the analysis of the electromagnetic coupling problems of conductive structures can be simplified.
基金supported by the Deanship of Scientific Research(DSR),King Abdulaziz University,Jeddah,Saudi Arabia
文摘Melting heat transfer in the boundary layer flow of a couple stress fluid over a stretching surface is investigated. The developed differential equations are solved for homotopic solutions. It is observed that the velocity and the boundary layer thickness are decreasing functions of the couple stress fluid parameter. However, the temperature and surface heat transfer increase when the values of the couple stress fluid parameter increase. The velocity and temperature fields increase with an increase in the melting process of the stretching sheet.
基金supported by the Fundamental Research Funds for the Central Universities(No.N090405009)
文摘A homotopy analysis method(HAM)is presented for the primary resonance of multiple degree-of-freedom systems with strong non-linearity excited by harmonic forces.The validity of the HAM is independent of the existence of small parameters in the considered equation.The HAM provides a simple way to adjust and control the convergence region of the series solution by means of an auxiliary parameter.Two examples are presented to show that the HAM solutions agree well with the results of the modified Linstedt-Poincar'e method and the incremental harmonic balance method.
基金Project supported by the National Natural Science Foundations of China(Grant Nos.10735030,10475055,10675065 and 90503006)the National Basic Research Program of China(Grant No.2007CB814800)
文摘This paper studies the generalized Kawahara equation in terms of the approximate homotopy symmetry method and the approximate homotopy direct method. Using both methods it obtains the similarity reduction solutions and similarity reduction equations of different orders, showing that the approximate homotopy direct method yields more general approximate similarity reductions than the approximate homotopy symmetry method. The homotopy series solutions to the generalized Kawahara equation are consequently derived.
基金Project supported by the National Natural Science Foundation of China(Nos.11072157,11272219,11227201,and 10932006)the National Basic Research Program of China(No.2012CB723301)the Training Program for Leading Talent in University Innovative Research Team in Hebei Province(No.LJRC006)
文摘Free and steady state forced transverse vibrations of non-uniform beams are investigated with a proposed method, leading to a series solution. The obtained series is verified to be convergent and linearly independent in a convergence test and by the non-zero value of the corresponding Wronski determinant, respectively. The obtained solution is rigorous, which can be reduced to a classical solution for uniform beams. The proposed method can deal with arbitrary non-uniform Euler-Bernoulli beams in principle, but the methods in terms of special functions or elementary functions can only work in some special cases.