Because of my carelessness,Eq.（1）in the paper "An approximate method for calculating the fluid force and response of a circular cylinder at lock-in"(China Ocean Engineering,22（3）,2008,pp.373)should be f...Because of my carelessness,Eq.（1）in the paper "An approximate method for calculating the fluid force and response of a circular cylinder at lock-in"(China Ocean Engineering,22（3）,2008,pp.373)should be f’-1.0/U’-5.0=f’;-1.0/5.75f’;-5.0,not f’=U’/5.75. My apology is hereby given.展开更多
In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem....In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem. the evaluation of the objective function is very difficult, so that only their approximate values can be obtained. This algorithm is obtained by combining penalty function method and approximation in bilevel programming. The presented algorithm is completely different from existing methods. That convergence for this algorithm is proved.展开更多
This paper 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.展开更多
The performance of structures with active variable stiffness (AVS) systems exhibits strong nonlinearity due to the variety with time of the stiffness of each storey unit,in which the AVS system is installed.Hence,the ...The performance of structures with active variable stiffness (AVS) systems exhibits strong nonlinearity due to the variety with time of the stiffness of each storey unit,in which the AVS system is installed.Hence,the classical dynamic analysis method for linear structures,such as the mode-superposition method,is not applicable to structures with AVS systems.In this paper,an approximate analysis method is proposed for displacement responses of structures with AVS systems.Firstly,an equivalent relationship between single-degree-of-freedom (SDOF) structures equipped with AVS systems and so-called fictitious linear structures is established.Then,an approximate mode-superposition (AMS) method is presented for multi-degree-of-freedom (MDOF) structures equipped with AVS systems.The accuracy of this method is investigated through extensive parametrical study using different types of earthquake excitations,and some modification is made to the method. Numerical calculation results indicate that the modified AMS method is effective for estimating the maximum displacements relative to the ground and the maximum interstorey drifts of MDOF structures equipped with AVS systems.展开更多
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.展开更多
The Kawahara equation is studied through the approximate homotopy symmetry method. Under this method we get the similarity reduction solutions of the Kawahara equation, leading to the corresponding homotopy series sol...The Kawahara equation is studied through the approximate homotopy symmetry method. Under this method we get the similarity reduction solutions of the Kawahara equation, leading to the corresponding homotopy series solutions. Furthermore, the similarity solutions of the corresponding reduced linear ordinary differential equations are also considered.展开更多
The approximate direct reduction method is applied to the perturbed mKdV equation with weak fourth order dispersion and weak dissipation. The similarity reduction solutions of different orders conform to formal cohere...The approximate direct reduction method is applied to the perturbed mKdV equation with weak fourth order dispersion and weak dissipation. The similarity reduction solutions of different orders conform to formal coherence, accounting for infinite series reduction solutions to the original equation and general formulas of similarity reduction equations. Painleve Ⅱ type equations, hyperbolic secant and Jacobi elliptic function solutions are obtained for zeroorder similarity reduction equations. Higher order similarity reduction equations are linear variable coefficient ordinary differential equations.展开更多
The perturbed Kaup-Kupershmidt equation is investigated in terms of the approximate symmetry perturbationmethod and the approximate direct method.The similarity reduction solutions of different orders are obtainedfor ...The perturbed Kaup-Kupershmidt equation is investigated in terms of the approximate symmetry perturbationmethod and the approximate direct method.The similarity reduction solutions of different orders are obtainedfor both methods, series reduction solutions are consequently derived.Higher order similarity reduction equations arelinear variable coefficients ordinary differential equations.By comparison, it is find that the results generated from theapproximate direct method are more general than the results generated from the approximate symmetry perturbationmethod.展开更多
The electron impact excitation(EIE) cross sections of an atom/ion in the whole energy region are needed in many research fields, such as astrophysics studies, inertial confinement fusion researches and so on. In the p...The electron impact excitation(EIE) cross sections of an atom/ion in the whole energy region are needed in many research fields, such as astrophysics studies, inertial confinement fusion researches and so on. In the present work, an effective method to calculate the EIE cross sections of an atom/ion in the whole energy region is presented. We use the EIE cross sections of helium as an illustration example. The optical forbidden 1^(1)S–n^(1)S(n = 2–4) and optical allowed 1^(1)S–n^(1)P(n = 2–4) excitation cross sections are calculated in the whole energy region using the scheme that combines the partial wave R-matrix method and the first Born approximation. The calculated cross sections are in good agreement with the available experimental measurements. Based on these accurate cross sections of our calculation, we find that the ratios between the accurate cross sections and Born cross sections are nearly the same for different excitation final states in the same channel. According to this interesting property, a universal correction function is proposed and given to calculate the accurate EIE cross sections with the same computational efforts of the widely used Born cross sections,which should be very useful in the related application fields. The datasets presented in this paper are openly available at https://www.doi.org/10.57760/sciencedb.j00113.00142.展开更多
Based on the boundary layer corrective method, a class of generalized nonlinear perturbed model in the critical case is studied. The asymptotic solution for the original equation is constructed. And the method is of s...Based on the boundary layer corrective method, a class of generalized nonlinear perturbed model in the critical case is studied. The asymptotic solution for the original equation is constructed. And the method is of significance to seek approximate solutions to other nonlinear models.展开更多
The behavior of beams with variable stiffness subjected to the action of variable loadings (impulse or harmonic) is analyzed in this paper using the successive approximation method. This successive approximation metho...The behavior of beams with variable stiffness subjected to the action of variable loadings (impulse or harmonic) is analyzed in this paper using the successive approximation method. This successive approximation method is a technique for numerical integration of partial differential equations involving both the space and time, with well-known initial conditions on time and boundary conditions on the space. This technique, although having been applied to beams with constant stiffness, is new for the case of beams with variable stiffness, and it aims to use a quadratic parabola (in time) to approximate the solutions of the differential equations of dynamics. The spatial part is studied using the successive approximation method of the partial differential equations obtained, in order to transform them into a system of time-dependent ordinary differential equations. Thus, the integration algorithm using this technique is established and applied to examples of beams with variable stiffness, under variable loading, and with the different cases of supports chosen in the literature. We have thus calculated the cases of beams with constant or variable rigidity with articulated or embedded supports, subjected to the action of an instantaneous impulse and harmonic loads distributed over its entire length. In order to justify the robustness of the successive approximation method considered in this work, an example of an articulated beam with constant stiffness subjected to a distributed harmonic load was calculated analytically, and the results obtained compared to those found numerically for various steps (spatial h and temporal τ ¯ ) of calculus, and the difference between the values obtained by the two methods was small. For example for ( h=1/8 , τ ¯ =1/ 64 ), the difference between these values is 17%.展开更多
This paper investigates an improved SIR model for COVID-19 based on the Caputo fractional derivative. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system...This paper investigates an improved SIR model for COVID-19 based on the Caputo fractional derivative. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system. Secondly, the stability of the system is discussed, among other things. Then, the GMMP method is introduced to obtain numerical solutions for the COVID-19 system. Numerical simulations were conducted using MATLAB, and the results indicate that our model is valuable for studying virus transmission.展开更多
In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseu...In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseudo-contractive mappings, and the set of solutions of the variational inclusion problem with multi-valued maximal monotone mappings and inverse-strongly monotone mappings in Hilbert space. Under suitable conditions, some strong convergence theorems are proved. Our results extends the recent results in G.L.Acedo and H.K.Xu [2], Zhang, Lee and Chan [8], Wakahashi and Toyoda [9], Takahashi and Takahashi [I0] and S. S. Chang, H. W. Joseph Lee and C. K. Chan [II], S.Takahashi and W.Takahashi [12]. Moreover, the method of proof adopted in this article is different from those of [4] and [12].展开更多
The aim of this paper is to employ fractional order proportional integral derivative(FO-PID)controller and integer order PID controller to control the position of the levitated object in a magnetic levitation system(M...The aim of this paper is to employ fractional order proportional integral derivative(FO-PID)controller and integer order PID controller to control the position of the levitated object in a magnetic levitation system(MLS),which is inherently nonlinear and unstable system.The proposal is to deploy discrete optimal pole-zero approximation method for realization of digital fractional order controller.An approach of phase shaping by slope cancellation of asymptotic phase plots for zeros and poles within given bandwidth is explored.The controller parameters are tuned using dynamic particle swarm optimization(d PSO)technique.Effectiveness of the proposed control scheme is verified by simulation and experimental results.The performance of realized digital FO-PID controller has been compared with that of the integer order PID controllers.It is observed that effort required in fractional order control is smaller as compared with its integer counterpart for obtaining the same system performance.展开更多
Let K be a closed convex subset of a real reflexive Banach space E, T：K→K be a nonexpansive mapping, and f：K→K be a fixed weakly contractive （may not be contractive） mapping. Then for any t∈（0, 1）, let x1∈K ...Let K be a closed convex subset of a real reflexive Banach space E, T：K→K be a nonexpansive mapping, and f：K→K be a fixed weakly contractive （may not be contractive） mapping. Then for any t∈（0, 1）, let x1∈K be the unique fixed point of the weak contraction x1→tf（x）＋（1-t）Tx. If T has a fixed point and E admits a weakly sequentially continuous duality mapping from E to E^＊, then it is shown that {xt} converges to a fixed point of T as t→0. The results presented here improve and generalize the corresponding results in （Xu, 2004）.展开更多
An intrinsic extension of Pad′e approximation method, called the generalized Pad′e approximation method, is proposed based on the classic Pad′e approximation theorem. According to the proposed method, the numerator...An intrinsic extension of Pad′e approximation method, called the generalized Pad′e approximation method, is proposed based on the classic Pad′e approximation theorem. According to the proposed method, the numerator and denominator of Pad′e approximant are extended from polynomial functions to a series composed of any kind of function, which means that the generalized Pad′e approximant is not limited to some forms, but can be constructed in different forms in solving different problems. Thus, many existing modifications of Pad′e approximation method can be considered to be the special cases of the proposed method. For solving homoclinic and heteroclinic orbits of strongly nonlinear autonomous oscillators, two novel kinds of generalized Pad′e approximants are constructed. Then, some examples are given to show the validity of the present method. To show the accuracy of the method, all solutions obtained in this paper are compared with those of the Runge–Kutta method.展开更多
Systemreliability sensitivity analysis becomes difficult due to involving the issues of the correlation between failure modes whether using analytic method or numerical simulation methods.A fast conditional reduction ...Systemreliability sensitivity analysis becomes difficult due to involving the issues of the correlation between failure modes whether using analytic method or numerical simulation methods.A fast conditional reduction method based on conditional probability theory is proposed to solve the sensitivity analysis based on the approximate analytic method.The relevant concepts are introduced to characterize the correlation between failure modes by the reliability index and correlation coefficient,and conditional normal fractile the for the multi-dimensional conditional failure analysis is proposed based on the two-dimensional normal distribution function.Thus the calculation of system failure probability can be represented as a summation of conditional probability terms,which is convenient to be computed by iterative solving sequentially.Further the system sensitivity solution is transformed into the derivation process of the failure probability correlation coefficient of each failure mode.Numerical examples results show that it is feasible to apply the idea of failure mode relevancy to failure probability sensitivity analysis,and it can avoid multi-dimension integral calculation and reduce complexity and difficulty.Compared with the product of conditional marginalmethod,a wider value range of correlation coefficient for reliability analysis is confirmed and an acceptable accuracy can be obtained with less computational cost.展开更多
In this paper, we firstly define a decreasing sequence {Pn(S)} by the generation of the Sierpinski gasket where each Pn(S) can be obtained in finite steps. Then we prove that the Hausdorff measure Hs(S) of the Sierpin...In this paper, we firstly define a decreasing sequence {Pn(S)} by the generation of the Sierpinski gasket where each Pn(S) can be obtained in finite steps. Then we prove that the Hausdorff measure Hs(S) of the Sierpinski gasket S can be approximated by {Pn(S)} with Pn(S)/(l + l/2n-3)s≤Hs(S)≤ Pn(S). An algorithm is presented to get Pn(S) for n ≤5. As an application, we obtain the best lower bound of Hs(S) till now: Hs(S)≥0.5631.展开更多
A formulation of a differential equation as projection and fixed point pi-Mem alloivs approximations using general piecnvise functions. We prone existence and uniqueness of the up proximate solution* convergence in th...A formulation of a differential equation as projection and fixed point pi-Mem alloivs approximations using general piecnvise functions. We prone existence and uniqueness of the up proximate solution* convergence in the L2 norm and nodal supercnnvergence. These results generalize those obtained earlier by Hulme for continuous piecevjise polynomials and by Delfour-Dubeau for discontinuous pieceuiise polynomials. A duality relationship for the two types of approximations is also given.展开更多
文摘Because of my carelessness,Eq.（1）in the paper "An approximate method for calculating the fluid force and response of a circular cylinder at lock-in"(China Ocean Engineering,22（3）,2008,pp.373)should be f’-1.0/U’-5.0=f’;-1.0/5.75f’;-5.0,not f’=U’/5.75. My apology is hereby given.
文摘In this paper, a new algorithm-approximate penalty function method is designed, which can be used to solve a bilevel optimization problem with linear constrained function. In this kind of bilevel optimization problem. the evaluation of the objective function is very difficult, so that only their approximate values can be obtained. This algorithm is obtained by combining penalty function method and approximation in bilevel programming. The presented algorithm is completely different from existing methods. That convergence for this algorithm is proved.
基金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.
基金National Natural Science foundation of China,Grant number 59895410
文摘The performance of structures with active variable stiffness (AVS) systems exhibits strong nonlinearity due to the variety with time of the stiffness of each storey unit,in which the AVS system is installed.Hence,the classical dynamic analysis method for linear structures,such as the mode-superposition method,is not applicable to structures with AVS systems.In this paper,an approximate analysis method is proposed for displacement responses of structures with AVS systems.Firstly,an equivalent relationship between single-degree-of-freedom (SDOF) structures equipped with AVS systems and so-called fictitious linear structures is established.Then,an approximate mode-superposition (AMS) method is presented for multi-degree-of-freedom (MDOF) structures equipped with AVS systems.The accuracy of this method is investigated through extensive parametrical study using different types of earthquake excitations,and some modification is made to the method. Numerical calculation results indicate that the modified AMS method is effective for estimating the maximum displacements relative to the ground and the maximum interstorey drifts of MDOF structures equipped with AVS systems.
基金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.
基金Supported by the National Natural Science Foundations of China under Grant Nos.10735030,10475055,10675065,and 90503006National Basic Research Program of China (973 Program 2007CB814800)
文摘The Kawahara equation is studied through the approximate homotopy symmetry method. Under this method we get the similarity reduction solutions of the Kawahara equation, leading to the corresponding homotopy series solutions. Furthermore, the similarity solutions of the corresponding reduced linear ordinary differential equations are also considered.
基金supported by the National Natural Science Foundations of China (Grant Nos 10735030,10475055,10675065 and 90503006)National Basic Research Program of China (Grant No 2007CB814800)+2 种基金PCSIRT (Grant No IRT0734)the Research Fund of Postdoctoral of China (Grant No 20070410727)Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No 20070248120)
文摘The approximate direct reduction method is applied to the perturbed mKdV equation with weak fourth order dispersion and weak dissipation. The similarity reduction solutions of different orders conform to formal coherence, accounting for infinite series reduction solutions to the original equation and general formulas of similarity reduction equations. Painleve Ⅱ type equations, hyperbolic secant and Jacobi elliptic function solutions are obtained for zeroorder similarity reduction equations. Higher order similarity reduction equations are linear variable coefficient ordinary differential equations.
基金Supported by the National Natural Science Foundation of China under Grant Nos.10735030,10475055,10675065,and 90503006National Basic Research Program of China (973 Program 2007CB814800)
文摘The perturbed Kaup-Kupershmidt equation is investigated in terms of the approximate symmetry perturbationmethod and the approximate direct method.The similarity reduction solutions of different orders are obtainedfor both methods, series reduction solutions are consequently derived.Higher order similarity reduction equations arelinear variable coefficients ordinary differential equations.By comparison, it is find that the results generated from theapproximate direct method are more general than the results generated from the approximate symmetry perturbationmethod.
基金Project supported by the National Natural Science Foundation of China (Grant No. 12241410)。
文摘The electron impact excitation(EIE) cross sections of an atom/ion in the whole energy region are needed in many research fields, such as astrophysics studies, inertial confinement fusion researches and so on. In the present work, an effective method to calculate the EIE cross sections of an atom/ion in the whole energy region is presented. We use the EIE cross sections of helium as an illustration example. The optical forbidden 1^(1)S–n^(1)S(n = 2–4) and optical allowed 1^(1)S–n^(1)P(n = 2–4) excitation cross sections are calculated in the whole energy region using the scheme that combines the partial wave R-matrix method and the first Born approximation. The calculated cross sections are in good agreement with the available experimental measurements. Based on these accurate cross sections of our calculation, we find that the ratios between the accurate cross sections and Born cross sections are nearly the same for different excitation final states in the same channel. According to this interesting property, a universal correction function is proposed and given to calculate the accurate EIE cross sections with the same computational efforts of the widely used Born cross sections,which should be very useful in the related application fields. The datasets presented in this paper are openly available at https://www.doi.org/10.57760/sciencedb.j00113.00142.
基金Supported by the National Science Foundation of China（11071075）
文摘Based on the boundary layer corrective method, a class of generalized nonlinear perturbed model in the critical case is studied. The asymptotic solution for the original equation is constructed. And the method is of significance to seek approximate solutions to other nonlinear models.
文摘The behavior of beams with variable stiffness subjected to the action of variable loadings (impulse or harmonic) is analyzed in this paper using the successive approximation method. This successive approximation method is a technique for numerical integration of partial differential equations involving both the space and time, with well-known initial conditions on time and boundary conditions on the space. This technique, although having been applied to beams with constant stiffness, is new for the case of beams with variable stiffness, and it aims to use a quadratic parabola (in time) to approximate the solutions of the differential equations of dynamics. The spatial part is studied using the successive approximation method of the partial differential equations obtained, in order to transform them into a system of time-dependent ordinary differential equations. Thus, the integration algorithm using this technique is established and applied to examples of beams with variable stiffness, under variable loading, and with the different cases of supports chosen in the literature. We have thus calculated the cases of beams with constant or variable rigidity with articulated or embedded supports, subjected to the action of an instantaneous impulse and harmonic loads distributed over its entire length. In order to justify the robustness of the successive approximation method considered in this work, an example of an articulated beam with constant stiffness subjected to a distributed harmonic load was calculated analytically, and the results obtained compared to those found numerically for various steps (spatial h and temporal τ ¯ ) of calculus, and the difference between the values obtained by the two methods was small. For example for ( h=1/8 , τ ¯ =1/ 64 ), the difference between these values is 17%.
文摘This paper investigates an improved SIR model for COVID-19 based on the Caputo fractional derivative. Firstly, the properties of the model are studied, including the feasible domain and bounded solutions of the system. Secondly, the stability of the system is discussed, among other things. Then, the GMMP method is introduced to obtain numerical solutions for the COVID-19 system. Numerical simulations were conducted using MATLAB, and the results indicate that our model is valuable for studying virus transmission.
基金supported by Scientific Research Fund of Sichuan Provincial Education Department (09ZB102)Scientific Research Fund of Science and Technology Deportment of Sichuan Provincial (2011JYZ011)
文摘In this article, we introduce a hybrid iterative scheme for finding a common element of the set of solutions for a generalized equilibrium problems, the set of common fixed point for a family of infinite k-strict pseudo-contractive mappings, and the set of solutions of the variational inclusion problem with multi-valued maximal monotone mappings and inverse-strongly monotone mappings in Hilbert space. Under suitable conditions, some strong convergence theorems are proved. Our results extends the recent results in G.L.Acedo and H.K.Xu [2], Zhang, Lee and Chan [8], Wakahashi and Toyoda [9], Takahashi and Takahashi [I0] and S. S. Chang, H. W. Joseph Lee and C. K. Chan [II], S.Takahashi and W.Takahashi [12]. Moreover, the method of proof adopted in this article is different from those of [4] and [12].
基金supported by the Board of Research in Nuclear Sciences of the Department of Atomic Energy,India(2012/36/69-BRNS/2012)
文摘The aim of this paper is to employ fractional order proportional integral derivative(FO-PID)controller and integer order PID controller to control the position of the levitated object in a magnetic levitation system(MLS),which is inherently nonlinear and unstable system.The proposal is to deploy discrete optimal pole-zero approximation method for realization of digital fractional order controller.An approach of phase shaping by slope cancellation of asymptotic phase plots for zeros and poles within given bandwidth is explored.The controller parameters are tuned using dynamic particle swarm optimization(d PSO)technique.Effectiveness of the proposed control scheme is verified by simulation and experimental results.The performance of realized digital FO-PID controller has been compared with that of the integer order PID controllers.It is observed that effort required in fractional order control is smaller as compared with its integer counterpart for obtaining the same system performance.
文摘Let K be a closed convex subset of a real reflexive Banach space E, T：K→K be a nonexpansive mapping, and f：K→K be a fixed weakly contractive （may not be contractive） mapping. Then for any t∈（0, 1）, let x1∈K be the unique fixed point of the weak contraction x1→tf（x）＋（1-t）Tx. If T has a fixed point and E admits a weakly sequentially continuous duality mapping from E to E^＊, then it is shown that {xt} converges to a fixed point of T as t→0. The results presented here improve and generalize the corresponding results in （Xu, 2004）.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11172093 and 11372102)the Hunan Provincial Innovation Foundation for Postgraduate,China(Grant No.CX2012B159)
文摘An intrinsic extension of Pad′e approximation method, called the generalized Pad′e approximation method, is proposed based on the classic Pad′e approximation theorem. According to the proposed method, the numerator and denominator of Pad′e approximant are extended from polynomial functions to a series composed of any kind of function, which means that the generalized Pad′e approximant is not limited to some forms, but can be constructed in different forms in solving different problems. Thus, many existing modifications of Pad′e approximation method can be considered to be the special cases of the proposed method. For solving homoclinic and heteroclinic orbits of strongly nonlinear autonomous oscillators, two novel kinds of generalized Pad′e approximants are constructed. Then, some examples are given to show the validity of the present method. To show the accuracy of the method, all solutions obtained in this paper are compared with those of the Runge–Kutta method.
基金This research is supported by National Key Research and Development Project(Grant Number 2019YFD0901002)Also Natural Science Foundation of Liaoning Province(Grant Number 20170540105)Liaoning Province Education Foundation(Grant Number JL201913)are gratefully acknowledged.
文摘Systemreliability sensitivity analysis becomes difficult due to involving the issues of the correlation between failure modes whether using analytic method or numerical simulation methods.A fast conditional reduction method based on conditional probability theory is proposed to solve the sensitivity analysis based on the approximate analytic method.The relevant concepts are introduced to characterize the correlation between failure modes by the reliability index and correlation coefficient,and conditional normal fractile the for the multi-dimensional conditional failure analysis is proposed based on the two-dimensional normal distribution function.Thus the calculation of system failure probability can be represented as a summation of conditional probability terms,which is convenient to be computed by iterative solving sequentially.Further the system sensitivity solution is transformed into the derivation process of the failure probability correlation coefficient of each failure mode.Numerical examples results show that it is feasible to apply the idea of failure mode relevancy to failure probability sensitivity analysis,and it can avoid multi-dimension integral calculation and reduce complexity and difficulty.Compared with the product of conditional marginalmethod,a wider value range of correlation coefficient for reliability analysis is confirmed and an acceptable accuracy can be obtained with less computational cost.
文摘In this paper, we firstly define a decreasing sequence {Pn(S)} by the generation of the Sierpinski gasket where each Pn(S) can be obtained in finite steps. Then we prove that the Hausdorff measure Hs(S) of the Sierpinski gasket S can be approximated by {Pn(S)} with Pn(S)/(l + l/2n-3)s≤Hs(S)≤ Pn(S). An algorithm is presented to get Pn(S) for n ≤5. As an application, we obtain the best lower bound of Hs(S) till now: Hs(S)≥0.5631.
基金This research has been supported in part by the Natural Sciences and Engineering Research Council of Canada(Grant OGPIN-336)and by the"Ministere de l'Education du Quebec"(FCAR Grant-ER-0725)
文摘A formulation of a differential equation as projection and fixed point pi-Mem alloivs approximations using general piecnvise functions. We prone existence and uniqueness of the up proximate solution* convergence in the L2 norm and nodal supercnnvergence. These results generalize those obtained earlier by Hulme for continuous piecevjise polynomials and by Delfour-Dubeau for discontinuous pieceuiise polynomials. A duality relationship for the two types of approximations is also given.
