The influences of machining and misalignment errors play a very critical role in the performance of the anti-backlash double-roller enveloping hourglass worm gear(ADEHWG).However,a corresponding efficient method for e...The influences of machining and misalignment errors play a very critical role in the performance of the anti-backlash double-roller enveloping hourglass worm gear(ADEHWG).However,a corresponding efficient method for eliminating or reducing these errors on the tooth profile of the ADEHWG is seldom reported.The gear engagement equation and tooth profile equation for considering six different errors that could arise from the machining and gear misalignment are derived from the theories of differential geometry and gear meshing.Also,the tooth contact analysis(TCA) is used to systematically investigate the influence of the machining and misalignment errors on the contact curves and the tooth profile by means of numerical analysis and three-dimensional solid modeling.The research results show that vertical angular misalignment of the worm wheel(Δβ) has the strongest influences while the tooth angle error(Δα) has the weakest influences on the contact curves and the tooth profile.A novel efficient approach is proposed and used to minimize the effect of the errors in manufacturing by changing the radius of the grinding wheel and the approaching point of contact.The results from the TCA and the experiment demonstrate that this tooth profile design modification method can indeed reduce the machining and misalignment errors.This modification design method is helpful in understanding the manufacturing technology of the ADEHWG.展开更多
To understand the engine main bearings' working condition is important in order to improve the performance of engine. However, thermal effects and thermal effect deformations of engine main bearings are rarely consid...To understand the engine main bearings' working condition is important in order to improve the performance of engine. However, thermal effects and thermal effect deformations of engine main bearings are rarely considered simultaneously in most studies. A typical finite element model is selected and the effect of thermoelastohydrodynamic(TEHD) reaction on engine main bearings is investigated. The calculated method of main bearing's thermal hydrodynamic reaction and journal misalignment effect is finite difference method, and its deformation reaction is calculated by using finite element method. The oil film pressure is solved numerically with Reynolds boundary conditions when various bearing characteristics are calculated. The whole model considers a temperature-pressure-viscosity relationship for the lubricant, surface roughness effect, and also an angular misalignment between the journal and the bearing. Numerical simulations of operation of a typical I6 diesel engine main bearing is conducted and importance of several contributing factors in mixed lubrication is discussed. The performance characteristics of journal misaligned main bearings under elastohydrodynamic(EHD) and TEHD loads of an I6 diesel engine are received, and then the journal center orbit movement, minimum oil film thickness and maximum oil film pressure of main bearings are estimated over a wide range of engine operation. The model is verified through the comparison with other present models. The TEHD performance of engine main bearings with various effects under the influences of journal misalignment is revealed, this is helpful to understand EHD and TEHD effect of misaligned engine main bearings.展开更多
In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference me...In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.展开更多
This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element...This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element method based on piecewise linear approximating functions and we also use the barycentric quadrature rule to compute the stiffness matrix and the L2-norm.展开更多
To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltai...To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltaic and wind resources,has been carried out.This has been done using a new meta-heuristic algorithm,improved artificial rabbits optimization(IARO).In this study,the uncertainty of solar and wind energy sources is modeled using Hang’s two-point estimating method(TPEM).The IARO algorithm is applied to calculate the best capacity of hub energy equipment,such as solar and wind renewable energy sources,combined heat and power(CHP)systems,steamboilers,energy storage,and electric cars in the day-aheadmarket.The standard ARO algorithmis developed to mimic the foraging behavior of rabbits,and in this work,the algorithm’s effectiveness in avoiding premature convergence is improved by using the dystudynamic inertia weight technique.The proposed IARO-based scheduling framework’s performance is evaluated against that of traditional ARO,particle swarm optimization(PSO),and salp swarm algorithm(SSA).The findings show that,in comparison to previous approaches,the suggested meta-heuristic scheduling framework based on the IARO has increased energy profit in day-ahead electricity,gas,and heating markets by satisfying the operational and energy hub limitations.Additionally,the results show that TPEM approach dependability consideration decreased hub energy’s profit by 8.995%as compared to deterministic planning.展开更多
To solve the receding horizon control (RHC) problem in an online manner, a novel numerical method called the indirect Radau pseudospectral method (IRPM) is proposed in this paper. Based on calculus of variations a...To solve the receding horizon control (RHC) problem in an online manner, a novel numerical method called the indirect Radau pseudospectral method (IRPM) is proposed in this paper. Based on calculus of variations and the first-order necessary optimality condition, the RHC problem for linear time-varying (LTV) system is transformed into the two-point boundary value problem (TPBVP). The Radau pseudospectral approximation is employed to discretize the TPBVP into well-posed linear algebraic equations. The resulting linear algebraic equations are solved via a matrix partitioning approach afterwards to obtain the optimal feedback control law. For the nonlinear system, the linearization method or the quasi linearization method is employed to approximate the RHC problem with successive linear approximations. Subsequently, each linear problem is solved via the similar method which is used to solve the RHC problem for LTV system. Simulation results of three examples show that the IRPM is of high accuracy and of high compu- tation efficiency to solve the RHC problem and the stability of closed-loop systems is guaranteed.展开更多
The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″(t) + q(t)f(t,u(t),u′(t)) = 0, 0 〈 t 〈 1, u(0) = u″(0) = u′(1) = 0 is studied....The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″(t) + q(t)f(t,u(t),u′(t)) = 0, 0 〈 t 〈 1, u(0) = u″(0) = u′(1) = 0 is studied. The iterative schemes for approximating the solutions are obtained by applying a monotone iterative method.展开更多
This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matr...This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.展开更多
Based on the structural and mechanics analysis of aero-engines rotor system, the dynamic model of the flexible rotor system with multi-supports are presented in order to solve the bearing misalignment problem of rotor...Based on the structural and mechanics analysis of aero-engines rotor system, the dynamic model of the flexible rotor system with multi-supports are presented in order to solve the bearing misalignment problem of rotor system in aero-engines. The motion equations are derived through Lagrange method. The relationship between structural and mechanics characteristics parameters are built up. Finally, the dynamic influence of bearing misalignment on rotor system are divided into three kinds: additional rotor bending rigidity, additional bearing misalignment excitation force and additional imbalance. The equations suggest that additional imbalance excitation force activates the nonlinearity on rotor system and an additional 2x excitation force might appear.展开更多
To overcome the inaccuracy problem of the traditional wedge evaluation of steel plates and strips caused by the ran-domness of the thicknesses of two local points and improve the reliability of the wedge index,the dou...To overcome the inaccuracy problem of the traditional wedge evaluation of steel plates and strips caused by the ran-domness of the thicknesses of two local points and improve the reliability of the wedge index,the double-centroid method for the wedge evaluation was proposed,and a model based on the centroid theory was established.Meanwhile,an integral model for the discrete thickness values of the cross-section profiles was derived.The discussion focused on the distinct characteristics of the two-point method,asymmetric method,and double-centroid method in evaluating the asymmetric distribution of cross-sections.The three methods were employed to evaluate the wedge values of both the theoretical and measured cross-sections of steel plates and strips,and the accuracies of three wedge evaluation models were analyzed and discussed.The results showed that the double-centroid method objectively reflects the degree and variation characteristics of the wedge values of the cross-sections of steel plates and strips,and this method is feasible,reliable,and outstanding.展开更多
In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, w...In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.展开更多
This paper presents an enhanced version of the standard shooting method that enables problems with two unknown parameters to be solved.A novel approach is applied to the analysis of the natural vibrations of Euler-Ber...This paper presents an enhanced version of the standard shooting method that enables problems with two unknown parameters to be solved.A novel approach is applied to the analysis of the natural vibrations of Euler-Bernoulli beams.The proposed algorithm,named as two-parameter multiple shooting method,is a new powerful numerical tool for calculating the natural frequencies and modes of multi-segment prismatic and non-prismatic beams with different boundary conditions.The impact of the axial force and additional point masses is also taken into account.Due to the fact that the method is based directly on the fourth-order ordinary differential equation,the structures do not have to be divided into many small elements to obtain an accurate enough solution,even though the geometry is very complex.To verify the proposed method,three different examples are considered,i.e.,a three-segment non-prismatic beam,a prismatic column subject to non-uniformly distributed compressive loads,and a two-segment beam with an additional point mass.Numerical analyses are carried out with the software MATHEMATICA.The results are compared with the solutions computed by the commercial finite element program SOFiSTiK.Good agreement is achieved,which confirms the correctness and high effectiveness of the formulated algorithm.展开更多
In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis sh...In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis shows that the postproeess improves the order of convergence. Consequently, we obtain asymptotically exact aposteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.展开更多
For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan's element energy projection (EEP) method has ...For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan's element energy projection (EEP) method has the accuracy O(h^min{2k,k+4}) The theoretical analysis coincides the reported numerical results.展开更多
We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the p...We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the problem into an equivalent system of weakly singular integral equations.Then,a Legendre-based spectral collocation method is developed for solving the transformed system.Therefore,we can make good use of the advantages of the Gauss quadrature rule.We present the construction and analysis of the collocation method.These results can be indirectly applied to solve fractional optimal control problems by considering the corresponding Euler–Lagrange equations.Two numerical examples are given to confirm the convergence analysis and robustness of the scheme.展开更多
Based on the precise integration method (PIM), a coupling technique of the high order multiplication perturbation method (HOMPM) and the reduction method is proposed to solve variable coefficient singularly pertur...Based on the precise integration method (PIM), a coupling technique of the high order multiplication perturbation method (HOMPM) and the reduction method is proposed to solve variable coefficient singularly perturbed two-point boundary value prob lems (TPBVPs) with one boundary layer. First, the inhomogeneous ordinary differential equations (ODEs) are transformed into the homogeneous ODEs by variable coefficient dimensional expansion. Then, the whole interval is divided evenly, and the transfer ma trix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efficient.展开更多
This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singula...This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations (ODEs) are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.展开更多
Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Br...Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Bratu’s equation, Troesch’s problems) occurs engineering and science, including the modeling of chemical reactions diffusion processes and heat transfer. An analytical expression pertaining to the concentration of substrate is obtained using Homotopy perturbation method for all values of parameters. These approximate analytical results were found to be in good agreement with the simulation results.展开更多
The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and...The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.展开更多
基金supported by National Natural Science Foundation of China(Grant Nos. 50775190No.51275425)+2 种基金Spring Sunshine Plan of Ministry of Education of China(Grant No. 10202258)Talent Introduction of Xihua UniversityChina(Grant No. Z1220217)
文摘The influences of machining and misalignment errors play a very critical role in the performance of the anti-backlash double-roller enveloping hourglass worm gear(ADEHWG).However,a corresponding efficient method for eliminating or reducing these errors on the tooth profile of the ADEHWG is seldom reported.The gear engagement equation and tooth profile equation for considering six different errors that could arise from the machining and gear misalignment are derived from the theories of differential geometry and gear meshing.Also,the tooth contact analysis(TCA) is used to systematically investigate the influence of the machining and misalignment errors on the contact curves and the tooth profile by means of numerical analysis and three-dimensional solid modeling.The research results show that vertical angular misalignment of the worm wheel(Δβ) has the strongest influences while the tooth angle error(Δα) has the weakest influences on the contact curves and the tooth profile.A novel efficient approach is proposed and used to minimize the effect of the errors in manufacturing by changing the radius of the grinding wheel and the approaching point of contact.The results from the TCA and the experiment demonstrate that this tooth profile design modification method can indeed reduce the machining and misalignment errors.This modification design method is helpful in understanding the manufacturing technology of the ADEHWG.
基金Supported by National Science and Technology Support Program of China:Vibration and Noise Reduction Technology Research and Application of Bulldozers and Other Earth Moving Machinery(Grant No.2015BAF07B04)
文摘To understand the engine main bearings' working condition is important in order to improve the performance of engine. However, thermal effects and thermal effect deformations of engine main bearings are rarely considered simultaneously in most studies. A typical finite element model is selected and the effect of thermoelastohydrodynamic(TEHD) reaction on engine main bearings is investigated. The calculated method of main bearing's thermal hydrodynamic reaction and journal misalignment effect is finite difference method, and its deformation reaction is calculated by using finite element method. The oil film pressure is solved numerically with Reynolds boundary conditions when various bearing characteristics are calculated. The whole model considers a temperature-pressure-viscosity relationship for the lubricant, surface roughness effect, and also an angular misalignment between the journal and the bearing. Numerical simulations of operation of a typical I6 diesel engine main bearing is conducted and importance of several contributing factors in mixed lubrication is discussed. The performance characteristics of journal misaligned main bearings under elastohydrodynamic(EHD) and TEHD loads of an I6 diesel engine are received, and then the journal center orbit movement, minimum oil film thickness and maximum oil film pressure of main bearings are estimated over a wide range of engine operation. The model is verified through the comparison with other present models. The TEHD performance of engine main bearings with various effects under the influences of journal misalignment is revealed, this is helpful to understand EHD and TEHD effect of misaligned engine main bearings.
基金heprojectissupportedbyNNSFofChina (No .1 9972 0 39) .
文摘In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.
文摘This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element method based on piecewise linear approximating functions and we also use the barycentric quadrature rule to compute the stiffness matrix and the L2-norm.
基金This research is supported by the Deputyship forResearch&Innovation,Ministry of Education in Saudi Arabia under Project Number(IFP-2022-35).
文摘To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltaic and wind resources,has been carried out.This has been done using a new meta-heuristic algorithm,improved artificial rabbits optimization(IARO).In this study,the uncertainty of solar and wind energy sources is modeled using Hang’s two-point estimating method(TPEM).The IARO algorithm is applied to calculate the best capacity of hub energy equipment,such as solar and wind renewable energy sources,combined heat and power(CHP)systems,steamboilers,energy storage,and electric cars in the day-aheadmarket.The standard ARO algorithmis developed to mimic the foraging behavior of rabbits,and in this work,the algorithm’s effectiveness in avoiding premature convergence is improved by using the dystudynamic inertia weight technique.The proposed IARO-based scheduling framework’s performance is evaluated against that of traditional ARO,particle swarm optimization(PSO),and salp swarm algorithm(SSA).The findings show that,in comparison to previous approaches,the suggested meta-heuristic scheduling framework based on the IARO has increased energy profit in day-ahead electricity,gas,and heating markets by satisfying the operational and energy hub limitations.Additionally,the results show that TPEM approach dependability consideration decreased hub energy’s profit by 8.995%as compared to deterministic planning.
基金supported by the National Natural Science Foundation of China(Nos.61174221 and 61402039)
文摘To solve the receding horizon control (RHC) problem in an online manner, a novel numerical method called the indirect Radau pseudospectral method (IRPM) is proposed in this paper. Based on calculus of variations and the first-order necessary optimality condition, the RHC problem for linear time-varying (LTV) system is transformed into the two-point boundary value problem (TPBVP). The Radau pseudospectral approximation is employed to discretize the TPBVP into well-posed linear algebraic equations. The resulting linear algebraic equations are solved via a matrix partitioning approach afterwards to obtain the optimal feedback control law. For the nonlinear system, the linearization method or the quasi linearization method is employed to approximate the RHC problem with successive linear approximations. Subsequently, each linear problem is solved via the similar method which is used to solve the RHC problem for LTV system. Simulation results of three examples show that the IRPM is of high accuracy and of high compu- tation efficiency to solve the RHC problem and the stability of closed-loop systems is guaranteed.
基金Supported by the Natural Science Foundation of Zhejiang Province (Y605144)the XNF of Zhejiang University of Media and Communications (XN080012008034)
文摘The existence of nondecreasing positive solutions for the nonlinear third-order twopoint boundary value problem u′″(t) + q(t)f(t,u(t),u′(t)) = 0, 0 〈 t 〈 1, u(0) = u″(0) = u′(1) = 0 is studied. The iterative schemes for approximating the solutions are obtained by applying a monotone iterative method.
基金Project supported by the National Natural Science Foundation of China(No.10672194)the China-Russia Cooperative Project(the National Natural Science Foundation of China and the Russian Foundation for Basic Research)(No.10811120012)
文摘This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.
文摘Based on the structural and mechanics analysis of aero-engines rotor system, the dynamic model of the flexible rotor system with multi-supports are presented in order to solve the bearing misalignment problem of rotor system in aero-engines. The motion equations are derived through Lagrange method. The relationship between structural and mechanics characteristics parameters are built up. Finally, the dynamic influence of bearing misalignment on rotor system are divided into three kinds: additional rotor bending rigidity, additional bearing misalignment excitation force and additional imbalance. The equations suggest that additional imbalance excitation force activates the nonlinearity on rotor system and an additional 2x excitation force might appear.
基金supported by the National Natural Science Foundation of China (Grant No.52305405)the Natural Science Foundation Research Program of Shanxi Province (Grant No.202203021222121)+4 种基金the Major Project of Science and Technology of Shanxi Province (Grant No.20181102016)the Chinese Postdoctoral Science Foundation (Grant No.2021M702544)the Central Government Guides the Special Fund Projects of Local Scientific and Technological Development (YDZX20191400002149)the Open Project of Research Institute of Hai'an-Taiyuan University of Technology (Grant No.2023HA-TYUTKFYF008)the School Fund of Taiyuan University of Technology (Grant No.2022QN007).
文摘To overcome the inaccuracy problem of the traditional wedge evaluation of steel plates and strips caused by the ran-domness of the thicknesses of two local points and improve the reliability of the wedge index,the double-centroid method for the wedge evaluation was proposed,and a model based on the centroid theory was established.Meanwhile,an integral model for the discrete thickness values of the cross-section profiles was derived.The discussion focused on the distinct characteristics of the two-point method,asymmetric method,and double-centroid method in evaluating the asymmetric distribution of cross-sections.The three methods were employed to evaluate the wedge values of both the theoretical and measured cross-sections of steel plates and strips,and the accuracies of three wedge evaluation models were analyzed and discussed.The results showed that the double-centroid method objectively reflects the degree and variation characteristics of the wedge values of the cross-sections of steel plates and strips,and this method is feasible,reliable,and outstanding.
基金The Project was supported by National Natural Science Foundation of China
文摘In this paper, we consider the following problem:The quadratic spline collocation, with uniform mesh and the mid-knot points are taken as the collocation points for this problem is considered. With some assumptions, we have proved that the solution of the quadratic spline collocation for the nonlinear problem can be written as a series expansions in integer powers of the mesh-size parameter. This gives us a construction method for using Richardson’s extrapolation. When we have a set of approximate solution with different mesh-size parameter a solution with high accuracy can he obtained by Richardson’s extrapolation.
文摘This paper presents an enhanced version of the standard shooting method that enables problems with two unknown parameters to be solved.A novel approach is applied to the analysis of the natural vibrations of Euler-Bernoulli beams.The proposed algorithm,named as two-parameter multiple shooting method,is a new powerful numerical tool for calculating the natural frequencies and modes of multi-segment prismatic and non-prismatic beams with different boundary conditions.The impact of the axial force and additional point masses is also taken into account.Due to the fact that the method is based directly on the fourth-order ordinary differential equation,the structures do not have to be divided into many small elements to obtain an accurate enough solution,even though the geometry is very complex.To verify the proposed method,three different examples are considered,i.e.,a three-segment non-prismatic beam,a prismatic column subject to non-uniformly distributed compressive loads,and a two-segment beam with an additional point mass.Numerical analyses are carried out with the software MATHEMATICA.The results are compared with the solutions computed by the commercial finite element program SOFiSTiK.Good agreement is achieved,which confirms the correctness and high effectiveness of the formulated algorithm.
基金supported partially by the innovation fund of Shanghai Normal Universitysupported partially by NSERC of Canada under Grant OGP0046726.
文摘In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis shows that the postproeess improves the order of convergence. Consequently, we obtain asymptotically exact aposteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.
基金Project supported by the National Natural Science Foundation of China (Nos. 10571046, 10371038)
文摘For a class of two-point boundary value problems, by virtue of onedimensional projection interpolation, it is proved that the nodal recovery derivative obtained by Yuan's element energy projection (EEP) method has the accuracy O(h^min{2k,k+4}) The theoretical analysis coincides the reported numerical results.
基金The Russian Foundation for Basic Research(RFBR)Grant No.19-01-00019.
文摘We target here to solve numerically a class of nonlinear fractional two-point boundary value problems involving left-and right-sided fractional derivatives.The main ingredient of the proposed method is to recast the problem into an equivalent system of weakly singular integral equations.Then,a Legendre-based spectral collocation method is developed for solving the transformed system.Therefore,we can make good use of the advantages of the Gauss quadrature rule.We present the construction and analysis of the collocation method.These results can be indirectly applied to solve fractional optimal control problems by considering the corresponding Euler–Lagrange equations.Two numerical examples are given to confirm the convergence analysis and robustness of the scheme.
基金Project supported by the National Natural Science Foundation of China(Key Program)(Nos.11132004 and 51078145)
文摘Based on the precise integration method (PIM), a coupling technique of the high order multiplication perturbation method (HOMPM) and the reduction method is proposed to solve variable coefficient singularly perturbed two-point boundary value prob lems (TPBVPs) with one boundary layer. First, the inhomogeneous ordinary differential equations (ODEs) are transformed into the homogeneous ODEs by variable coefficient dimensional expansion. Then, the whole interval is divided evenly, and the transfer ma trix in each sub-interval is worked out through the HOMPM. Finally, a group of algebraic equations are given based on the relationship between the neighboring sub-intervals, which are solved by the reduction method. Numerical results show that the present method is highly efficient.
基金supported by the National Natural Science Foundation of China(Key Program)(Nos.11132004 and 51078145)
文摘This paper presents a high order multiplication perturbation method for sin- gularly perturbed two-point boundary value problems with the boundary layer at one end. By the theory of singular perturbations, the singularly perturbed two-point boundary value problems are first transformed into the singularly perturbed initial value problems. With the variable coefficient dimensional expanding, the non-homogeneous ordinary dif- ferential equations (ODEs) are transformed into the homogeneous ODEs, which are then solved by the high order multiplication perturbation method. Some linear and nonlinear numerical examples show that the proposed method has high precision.
文摘Several problems arising in science and engineering are modeled by differential equations that involve conditions that are specified at more than one point. The non-linear two-point boundary value problem (TPBVP) (Bratu’s equation, Troesch’s problems) occurs engineering and science, including the modeling of chemical reactions diffusion processes and heat transfer. An analytical expression pertaining to the concentration of substrate is obtained using Homotopy perturbation method for all values of parameters. These approximate analytical results were found to be in good agreement with the simulation results.
文摘The present paper proposes a mathematical method to numerically treat a class of third-order linear Boundary Value Problems (BVPs). This method is based on the combination of the Adomian Decomposition Method (ADM) and, the modified shooting method. A complete derivation of the proposed method has been provided, in addition to its numerical implementation and, validation via the utilization of the Runge-Kutta method and, other existing methods. The method has been applied to diverse test problems and turned out to perform remarkably. Lastly, the simulated numerical results have been graphically illustrated and, also supported by some absolute error comparison tables.