A new method for approximating the inerse Laplace transform is presented. We first change our Laplace transform equation into a convolution type integral equation, where Tikhonov regularization techniques and the Four...A new method for approximating the inerse Laplace transform is presented. We first change our Laplace transform equation into a convolution type integral equation, where Tikhonov regularization techniques and the Fourier transformation are easily applied. We finally obtain a regularized approximation to the inverse Laplace transform as finite sum展开更多
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.展开更多
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%.展开更多
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.展开更多
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.展开更多
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.展开更多
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).展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes ...The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes fan be assumed equal to arbitrary points,where the integrand function f is known; iii) the number of the requested evaluations of f at the nodes is low,iv) a satisfactory convergence theory can be proved.展开更多
The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, w...The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares (MLS) approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.展开更多
In this work, by choosing an orthonormal basis for the Hilbert space L^2[0, 1], an approximation method for finding approximate solutions of the equation (I + K)x = y is proposed, called Haar wavelet approximation ...In this work, by choosing an orthonormal basis for the Hilbert space L^2[0, 1], an approximation method for finding approximate solutions of the equation (I + K)x = y is proposed, called Haar wavelet approximation method (HWAM). To prove the applicabifity of the HWAM, a more general applicability theorem on an approximation method (AM) for an operator equation Ax = y is proved first. As an application, applicability of the HWAM is obtained. Fhrthermore, four steps to use the HWAM are listed and three numerical examples are given in order to illustrate the effectiveness of the method.展开更多
An industrial building is a non-classically damped system due to the different damping properties of the primary structure and equipment.The objective of this paper is to quantify the range of applicability of the rea...An industrial building is a non-classically damped system due to the different damping properties of the primary structure and equipment.The objective of this paper is to quantify the range of applicability of the real model superposition approximation method to the seismic response calculation of industrial buildings.The analysis using lumped mass-and-shear spring models indicates that for the equipment-to-structure frequency ratiosγf>1.1 orγf<0.9,the non-classical damping effect is limited,and the real mode superposition approximation method provides accurate estimates.For 0.9<γf<1.1,the system may have a pair of closely spaced frequency modes,and the non-zero off-diagonal damping terms have a non-negligible effect on the damping ratios and mode shape vectors of these modes.For 0.9<γf<1.1 and the equipment-to-structure mass ratiosγm<0.07,the real mode superposition approximation method results in large errors,while the approximation method can provide an accurate estimation for 0.9<γf<1.1 andγm>0.07.Furthermore,extensive parametric analyses are conducted,where both steel structures and reinforced concrete structures with equipment with various damping ratios are considered.Finally,the finite element analysis of a five-story industrial building is adopted to validate the proposed range of applicability.展开更多
In the reliability analysis of complex structures,response surface method(RSM)has been suggested as an efficient technique to estimate the actual but implicit limit state function.A set of sample points are needed to ...In the reliability analysis of complex structures,response surface method(RSM)has been suggested as an efficient technique to estimate the actual but implicit limit state function.A set of sample points are needed to fit to the implicit function.It has been noted that the accuracy of RSM depends highly on the so-called sample points.However,the technique for point selection has had little attention.In the present study,an improved response surface method(IRSM)based on two sample point selection techniques,named the direction cosines projected strategy(DCS)and the limit step length iteration strategy(LSS),is investigated.Since it uses the sampling points selected to be located in the region close to the original failure surface,and since it needs only one response surface,the IRSM should be accurate and simple in practical structural problems.Applications to several typical examples have helped to elucidate the successful working of the IRSM.展开更多
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, 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.展开更多
文摘A new method for approximating the inerse Laplace transform is presented. We first change our Laplace transform equation into a convolution type integral equation, where Tikhonov regularization techniques and the Fourier transformation are easily applied. We finally obtain a regularized approximation to the inverse Laplace transform as finite sum
基金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.
文摘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%.
基金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.
文摘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 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.
文摘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).
文摘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.
基金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.
基金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.
基金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.
基金Work sponsored by"Ministero dell' University"CNR of Italy
文摘The purpose of this paper is to propose and study local spline approximation methods for singular product integration,for which;i)the precision degree is the highest possible using splint approximation; ii) the nodes fan be assumed equal to arbitrary points,where the integrand function f is known; iii) the number of the requested evaluations of f at the nodes is low,iv) a satisfactory convergence theory can be proved.
基金supported by the National Natural Science Foundation of China (11171208)Shanghai Leading Academic Discipline Project (S30106)
文摘The paper presents the improved element-free Galerkin (IEFG) method for three-dimensional wave propa- gation. The improved moving least-squares (IMLS) approx- imation is employed to construct the shape function, which uses an orthogonal function system with a weight function as the basis function. Compared with the conventional moving least-squares (MLS) approximation, the algebraic equation system in the IMLS approximation is not ill-conditioned, and can be solved directly without deriving the inverse matrix. Because there are fewer coefficients in the IMLS than in the MLS approximation, fewer nodes are selected in the IEFG method than in the element-free Galerkin method. Thus, the IEFG method has a higher computing speed. In the IEFG method, the Galerkin weak form is employed to obtain a dis- cretized system equation, and the penalty method is applied to impose the essential boundary condition. The traditional difference method for two-point boundary value problems is selected for the time discretization. As the wave equations and the boundary-initial conditions depend on time, the scal- ing parameter, number of nodes and the time step length are considered for the convergence study.
基金support by the NSFC(11371012,11401359,11471200)the FRF for the Central Universities(GK201301007)the NSRP of Shaanxi Province(2014JQ1010)
文摘In this work, by choosing an orthonormal basis for the Hilbert space L^2[0, 1], an approximation method for finding approximate solutions of the equation (I + K)x = y is proposed, called Haar wavelet approximation method (HWAM). To prove the applicabifity of the HWAM, a more general applicability theorem on an approximation method (AM) for an operator equation Ax = y is proved first. As an application, applicability of the HWAM is obtained. Fhrthermore, four steps to use the HWAM are listed and three numerical examples are given in order to illustrate the effectiveness of the method.
基金Fund of China National Industrial Building Diagnosis and Reconstruction Engineering Technology Research Center under Grant No.YZA2017Ky03the Beijing Natural Science Foundation under Grant No.JQ18029the National Natural Science Foundation of China under Grant No.52078277。
文摘An industrial building is a non-classically damped system due to the different damping properties of the primary structure and equipment.The objective of this paper is to quantify the range of applicability of the real model superposition approximation method to the seismic response calculation of industrial buildings.The analysis using lumped mass-and-shear spring models indicates that for the equipment-to-structure frequency ratiosγf>1.1 orγf<0.9,the non-classical damping effect is limited,and the real mode superposition approximation method provides accurate estimates.For 0.9<γf<1.1,the system may have a pair of closely spaced frequency modes,and the non-zero off-diagonal damping terms have a non-negligible effect on the damping ratios and mode shape vectors of these modes.For 0.9<γf<1.1 and the equipment-to-structure mass ratiosγm<0.07,the real mode superposition approximation method results in large errors,while the approximation method can provide an accurate estimation for 0.9<γf<1.1 andγm>0.07.Furthermore,extensive parametric analyses are conducted,where both steel structures and reinforced concrete structures with equipment with various damping ratios are considered.Finally,the finite element analysis of a five-story industrial building is adopted to validate the proposed range of applicability.
文摘In the reliability analysis of complex structures,response surface method(RSM)has been suggested as an efficient technique to estimate the actual but implicit limit state function.A set of sample points are needed to fit to the implicit function.It has been noted that the accuracy of RSM depends highly on the so-called sample points.However,the technique for point selection has had little attention.In the present study,an improved response surface method(IRSM)based on two sample point selection techniques,named the direction cosines projected strategy(DCS)and the limit step length iteration strategy(LSS),is investigated.Since it uses the sampling points selected to be located in the region close to the original failure surface,and since it needs only one response surface,the IRSM should be accurate and simple in practical structural problems.Applications to several typical examples have helped to elucidate the successful working of the IRSM.
文摘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.
基金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.