This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatia...We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.展开更多
In this work,a consistent and physically accurate implementation of the general framework of unified second-order time accurate integrators via the well-known GSSSS framework in the Discrete Element Method is presente...In this work,a consistent and physically accurate implementation of the general framework of unified second-order time accurate integrators via the well-known GSSSS framework in the Discrete Element Method is presented.The improved tangential displacement evaluation in the present implementation of the discrete element method has been derived and implemented to preserve the consistency of the correct time level evaluation during the time integration process in calculating the algorithmic tangential displacement.Several numerical examples have been used to validate the proposed tangential displacement evaluation;this is in contrast to past practices which only seem to attain the first-order time accuracy due to inconsistent time level implementation with different algorithms for normal and tangential directions.The comparisons with the existing implementation and the superiority of the proposed implementation are given in terms of the convergence rate with improved numerical accuracy in time.Moreover,several schemes via the unified second-order time integrators within the framework of the GSSSS family have been carried out based on the proposed correct implementation.All the numerical results demonstrate that using the existing state-of-the-art implementation reduces the time accuracy to be first-order accurate in time,while the proposed implementation preserves the correct time accuracy to yield second-order.展开更多
There are two models in use today to analyze structural responses when subjected to earthquake ground motions, the Displacement Input Model (DIM) and the Acceleration Input Model (AIM). The time steps used in dire...There are two models in use today to analyze structural responses when subjected to earthquake ground motions, the Displacement Input Model (DIM) and the Acceleration Input Model (AIM). The time steps used in direct integration methods for these models are analyzed to examine the suitability of DIM. Numerical results are presented and show that the time-step for DIM is about the same as for AIM, and achieves the same accuracy. This is contrary to previous research that reported that there are several sources of numerical errors associated with the direct application of earthquake displacement loading, and a very small time step is required to define the displacement record and to integrate the dynamic equilibrium equation. It is shown in this paper that DIM is as accurate and suitable as, if not more than, AIM for analyzing the response of a structure to uniformly distributed and spatially varying ground motions.展开更多
The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficu...The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method (e.g., an algorithm with an arbitrary order of accuracy) is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.展开更多
In this paper,we discuss on the convergence and approximation of an α times integrated semigroups. The Trotter kato theorems for an α times integrated semigroups are obtained.
To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of...To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.展开更多
Time delay and integration (TDI) charge coupled device (CCD) noise sets a fundamental limit on image sensor performance, especially under low illumination in remote sensing applications. After introducing the comp...Time delay and integration (TDI) charge coupled device (CCD) noise sets a fundamental limit on image sensor performance, especially under low illumination in remote sensing applications. After introducing the complete sources of CCD noise, we study the effects of TDI operation mode on noise, and the relationship between different types of noise and number of the TDI stage. Then we propose a new technique to identify and measure sources of TDI CCD noise employing mathematical statistics theory, where theoretical analysis shows that noise estimated formulation converges well. Finally, we establish a testing platform to carry out experiments, and a standard TDI CCD is calibrated by using the proposed method. The experimental results show that the noise analysis and measurement methods presented in this paper are useful for modeling TDI CCDs.展开更多
Finite element(FE) is a powerful tool and has been applied by investigators to real-time hybrid simulations(RTHSs). This study focuses on the computational efficiency, including the computational time and accuracy...Finite element(FE) is a powerful tool and has been applied by investigators to real-time hybrid simulations(RTHSs). This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time(TET) decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method(CDM), the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ(λ = 4) method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.展开更多
This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditiona...This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai's approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.展开更多
The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integ...The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.展开更多
This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order...This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a "major compo- nent" and a "high order small quantity" and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly.展开更多
The computational uncertainty principle states that the numerical computation of nonlinear ordinary differential equations(ODEs) should use appropriately sized time steps to obtain reliable solutions.However,the int...The computational uncertainty principle states that the numerical computation of nonlinear ordinary differential equations(ODEs) should use appropriately sized time steps to obtain reliable solutions.However,the interval of effective step size(IES) has not been thoroughly explored theoretically.In this paper,by using a general estimation for the total error of the numerical solutions of ODEs,a method is proposed for determining an approximate IES by translating the functions for truncation and rounding errors.It also illustrates this process with an example.Moreover,the relationship between the IES and its approximation is found,and the relative error of the approximation with respect to the IES is given.In addition,variation in the IES with increasing integration time is studied,which can provide an explanation for the observed numerical results.The findings contribute to computational step-size choice for reliable numerical solutions.展开更多
Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in trac...Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in tracking long-term dynamics.For improving such a higher-order accurate algorithm,this paper proposes a two sub-step higher-order algorithm with unconditional stability and controllable dissipation.In the proposed algorithm,a time step interval[t_(k),t_(k)+h]where h stands for the size of a time step is divided into two sub-steps[t_(k),t_(k)+γh]and[t_(k)+γh,t_(k)+h].A non-dissipative fourth-order algorithm is used in the rst sub-step to ensure low-frequency accuracy and a dissipative third-order algorithm is employed in the second sub-step to lter out the contribution of high-frequency modes.Besides,two approaches are used to design the algorithm parameterγ.The rst approach determinesγby maximizing low-frequency accuracy and the other determinesγfor quickly damping out highfrequency modes.The present algorithm usesρ_(∞)to exactly control the degree of numerical dissipation,and it is third-order accurate when 0≤ρ_(∞)<1 and fourth-order accurate whenρ_(∞)=1.Furthermore,the proposed algorithm is self-starting and easy to implement.Some illustrative linear and nonlinear examples are solved to check the performances of the proposed two sub-step higher-order algorithm.展开更多
The optimal tracking performance for integrator and dead time plant in the case where plant uncertainty and control energy constraints are to be considered jointly is inrestigated. Firstly, an average cost function of...The optimal tracking performance for integrator and dead time plant in the case where plant uncertainty and control energy constraints are to be considered jointly is inrestigated. Firstly, an average cost function of the tracking error and the plant input energy over a class of stochastic model errors are defined. Then, we obtain an internal model controller design method that minimizes the average performance and further studies optimal tracking performance for integrator and dead time plant in the simultaneous presence of plant uncertainty and control energy constraint. The results can be used to evaluate optimal tracking performance and control energy in practical designs.展开更多
For the constrained nonlinear optimal control problem, by taking the first term of Taylor series, the dynamic equation is linearized. Thus by, introducing into the dual variable (Lagrange multiplier vector), the dynam...For the constrained nonlinear optimal control problem, by taking the first term of Taylor series, the dynamic equation is linearized. Thus by, introducing into the dual variable (Lagrange multiplier vector), the dynamic equation can be transformed into Hamilton system from Lagrange system on the basis of the original variable. Under the whole state, the problem discussed can be described from a new view, and the equation can be precisely solved by, the time precise integration method established in linear dynamic system. A numerical example shows the effectiveness of the method.展开更多
Due to the widespread application of the PID controller in industrial control systems, it is desirable to know the complete set of all the stabilizing PID controllers for a given plant before the controller design and...Due to the widespread application of the PID controller in industrial control systems, it is desirable to know the complete set of all the stabilizing PID controllers for a given plant before the controller design and tuning. In this paper, the stabilization problems of the classical proportionalintegral-derivative (PID) controller and the singleparameter PID controller (containing only one adjustable parameter) for integral processes with time delay are investigated, respectively. The complete set of stabilizing parameters of the classical PID controller is determined using a version of the Hermite-Biehler Theorem applicable to quasipolynomials. Since the stabilization problem of the singie-parameter PID controller cannot be treated by the Hermite-Biehler Theorem, a simple method called duallocus diagram is employed to derive the stabilizing range of the single-parameter PID controller. These results provide insight into the tuning of the PID controllers.展开更多
A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iter...A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iterative errors generated within one time step is constructed.With an estimate of temporal and spatial error using an embedded RungeKutta scheme and a higher order spatial discretization,an adaptive time-stepping strategy is proposed based on the idea that the time step should be the maximum without obviously infuencing the total error of the discretization.The designed adaptive time-stepping strategy is then tested in various types of problems including isentropic vortex convection,steady-state fow past a fat plate,Taylor-Green vortex and turbulent fow over a circular cylinder at Re=3900.The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efciency is obtained for unsteady and steady,well-resolved and under-resolved simulations.展开更多
We introduce a new transmit/receive dipole pair array to obtain a compact quasi\|monostatic antenna structure for ground penetrating radar systems. And we analyze this transmit/receive dipole ...We introduce a new transmit/receive dipole pair array to obtain a compact quasi\|monostatic antenna structure for ground penetrating radar systems. And we analyze this transmit/receive dipole pair array in time domain. The numerical results show that if the distance between the transmit antenna and receive antenna is appropriate the array configuration is adoptable.展开更多
A time domain electric al field integral equation (TDEFIE) is formulated for the problem of a thin wire antenna in the presence of conductor bodies, and this equation is solved by the me...A time domain electric al field integral equation (TDEFIE) is formulated for the problem of a thin wire antenna in the presence of conductor bodies, and this equation is solved by the method of time marching algorithm. The analysis is valid for any arbitrarily shaped, oriented and positioned wire antennas relative to arbitrarily shaped conductor bodies. Current at the excited point, input admittance and radiation pattern are given and agree with the results computed by the method in frequency domain.展开更多
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
文摘We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.
文摘In this work,a consistent and physically accurate implementation of the general framework of unified second-order time accurate integrators via the well-known GSSSS framework in the Discrete Element Method is presented.The improved tangential displacement evaluation in the present implementation of the discrete element method has been derived and implemented to preserve the consistency of the correct time level evaluation during the time integration process in calculating the algorithmic tangential displacement.Several numerical examples have been used to validate the proposed tangential displacement evaluation;this is in contrast to past practices which only seem to attain the first-order time accuracy due to inconsistent time level implementation with different algorithms for normal and tangential directions.The comparisons with the existing implementation and the superiority of the proposed implementation are given in terms of the convergence rate with improved numerical accuracy in time.Moreover,several schemes via the unified second-order time integrators within the framework of the GSSSS family have been carried out based on the proposed correct implementation.All the numerical results demonstrate that using the existing state-of-the-art implementation reduces the time accuracy to be first-order accurate in time,while the proposed implementation preserves the correct time accuracy to yield second-order.
文摘There are two models in use today to analyze structural responses when subjected to earthquake ground motions, the Displacement Input Model (DIM) and the Acceleration Input Model (AIM). The time steps used in direct integration methods for these models are analyzed to examine the suitability of DIM. Numerical results are presented and show that the time-step for DIM is about the same as for AIM, and achieves the same accuracy. This is contrary to previous research that reported that there are several sources of numerical errors associated with the direct application of earthquake displacement loading, and a very small time step is required to define the displacement record and to integrate the dynamic equilibrium equation. It is shown in this paper that DIM is as accurate and suitable as, if not more than, AIM for analyzing the response of a structure to uniformly distributed and spatially varying ground motions.
基金financial support from Hunan Provincial Natura1 Science Foundation of China,Grant Number:02JJY2085,for this study
文摘The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method (e.g., an algorithm with an arbitrary order of accuracy) is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.
文摘In this paper,we discuss on the convergence and approximation of an α times integrated semigroups. The Trotter kato theorems for an α times integrated semigroups are obtained.
基金The project supported by the National Key Program for Developing Basic Sciences (G1999043408 and G1998040901-1)the National Natural Sciences Foundation of China (40175024 and 40035010)
文摘To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme.
基金Project supported by the National High Technology Research and Development Program of China (Grant No. 2006AA06A208)
文摘Time delay and integration (TDI) charge coupled device (CCD) noise sets a fundamental limit on image sensor performance, especially under low illumination in remote sensing applications. After introducing the complete sources of CCD noise, we study the effects of TDI operation mode on noise, and the relationship between different types of noise and number of the TDI stage. Then we propose a new technique to identify and measure sources of TDI CCD noise employing mathematical statistics theory, where theoretical analysis shows that noise estimated formulation converges well. Finally, we establish a testing platform to carry out experiments, and a standard TDI CCD is calibrated by using the proposed method. The experimental results show that the noise analysis and measurement methods presented in this paper are useful for modeling TDI CCDs.
基金National Natural Science Foundation of China under Grant Nos.51639006 and 51725901
文摘Finite element(FE) is a powerful tool and has been applied by investigators to real-time hybrid simulations(RTHSs). This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time(TET) decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method(CDM), the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ(λ = 4) method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.
基金The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China (G2000048702, 2003CB716707)the National Science Fund for Distinguished Young Scholars (10025208)+1 种基金 the National Natural Science Foundation of China (Key Program) (10532040) the Research Fund for 0versea Chinese (10228028).
文摘This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai's approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.
文摘The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.
基金supported by the National Natural Science Foun-dation of China (11172334)
文摘This paper presents a high order symplectic con- servative perturbation method for linear time-varying Hamil- tonian system. Firstly, the dynamic equation of Hamilto- nian system is gradually changed into a high order pertur- bation equation, which is solved approximately by resolv- ing the Hamiltonian coefficient matrix into a "major compo- nent" and a "high order small quantity" and using perturba- tion transformation technique, then the solution to the orig- inal equation of Hamiltonian system is determined through a series of inverse transform. Because the transfer matrix determined by the method in this paper is the product of a series of exponential matrixes, the transfer matrix is a sym- plectic matrix; furthermore, the exponential matrices can be calculated accurately by the precise time integration method, so the method presented in this paper has fine accuracy, ef- ficiency and stability. The examples show that the proposed method can also give good results even though a large time step is selected, and with the increase of the perturbation or- der, the perturbation solutions tend to exact solutions rapidly.
基金supported by the National Natural Science Foundation of China[grant numbers 41375110,11471244]
文摘The computational uncertainty principle states that the numerical computation of nonlinear ordinary differential equations(ODEs) should use appropriately sized time steps to obtain reliable solutions.However,the interval of effective step size(IES) has not been thoroughly explored theoretically.In this paper,by using a general estimation for the total error of the numerical solutions of ODEs,a method is proposed for determining an approximate IES by translating the functions for truncation and rounding errors.It also illustrates this process with an example.Moreover,the relationship between the IES and its approximation is found,and the relative error of the approximation with respect to the IES is given.In addition,variation in the IES with increasing integration time is studied,which can provide an explanation for the observed numerical results.The findings contribute to computational step-size choice for reliable numerical solutions.
基金supported by the National Natural Science Foundation of China(Grant Numbers 11872090,11672019,11472035).
文摘Based on the weighted residual method,a single-step time integration algorithm with higher-order accuracy and unconditional stability has been proposed,which is superior to the second-order accurate algorithms in tracking long-term dynamics.For improving such a higher-order accurate algorithm,this paper proposes a two sub-step higher-order algorithm with unconditional stability and controllable dissipation.In the proposed algorithm,a time step interval[t_(k),t_(k)+h]where h stands for the size of a time step is divided into two sub-steps[t_(k),t_(k)+γh]and[t_(k)+γh,t_(k)+h].A non-dissipative fourth-order algorithm is used in the rst sub-step to ensure low-frequency accuracy and a dissipative third-order algorithm is employed in the second sub-step to lter out the contribution of high-frequency modes.Besides,two approaches are used to design the algorithm parameterγ.The rst approach determinesγby maximizing low-frequency accuracy and the other determinesγfor quickly damping out highfrequency modes.The present algorithm usesρ_(∞)to exactly control the degree of numerical dissipation,and it is third-order accurate when 0≤ρ_(∞)<1 and fourth-order accurate whenρ_(∞)=1.Furthermore,the proposed algorithm is self-starting and easy to implement.Some illustrative linear and nonlinear examples are solved to check the performances of the proposed two sub-step higher-order algorithm.
基金the High Technology Research and Development (863) Program (2003AA517020).
文摘The optimal tracking performance for integrator and dead time plant in the case where plant uncertainty and control energy constraints are to be considered jointly is inrestigated. Firstly, an average cost function of the tracking error and the plant input energy over a class of stochastic model errors are defined. Then, we obtain an internal model controller design method that minimizes the average performance and further studies optimal tracking performance for integrator and dead time plant in the simultaneous presence of plant uncertainty and control energy constraint. The results can be used to evaluate optimal tracking performance and control energy in practical designs.
文摘For the constrained nonlinear optimal control problem, by taking the first term of Taylor series, the dynamic equation is linearized. Thus by, introducing into the dual variable (Lagrange multiplier vector), the dynamic equation can be transformed into Hamilton system from Lagrange system on the basis of the original variable. Under the whole state, the problem discussed can be described from a new view, and the equation can be precisely solved by, the time precise integration method established in linear dynamic system. A numerical example shows the effectiveness of the method.
基金National Science Foundation of China (60274032) SRFDP (20030248040) SRSP (04QMH1405)
文摘Due to the widespread application of the PID controller in industrial control systems, it is desirable to know the complete set of all the stabilizing PID controllers for a given plant before the controller design and tuning. In this paper, the stabilization problems of the classical proportionalintegral-derivative (PID) controller and the singleparameter PID controller (containing only one adjustable parameter) for integral processes with time delay are investigated, respectively. The complete set of stabilizing parameters of the classical PID controller is determined using a version of the Hermite-Biehler Theorem applicable to quasipolynomials. Since the stabilization problem of the singie-parameter PID controller cannot be treated by the Hermite-Biehler Theorem, a simple method called duallocus diagram is employed to derive the stabilizing range of the single-parameter PID controller. These results provide insight into the tuning of the PID controllers.
基金Zhen-Guo Yan acknowledges supports from the National Natural Science Foundation of China(Grant no.11902344)National Numerical Windtunnel Project.The development of the implicit solver in Nektar++has been supported by EPSRC grant(EP/R029423/1)UK Turbulence Consortium grant(EP/R029326/1).
文摘A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iterative errors generated within one time step is constructed.With an estimate of temporal and spatial error using an embedded RungeKutta scheme and a higher order spatial discretization,an adaptive time-stepping strategy is proposed based on the idea that the time step should be the maximum without obviously infuencing the total error of the discretization.The designed adaptive time-stepping strategy is then tested in various types of problems including isentropic vortex convection,steady-state fow past a fat plate,Taylor-Green vortex and turbulent fow over a circular cylinder at Re=3900.The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efciency is obtained for unsteady and steady,well-resolved and under-resolved simulations.
文摘We introduce a new transmit/receive dipole pair array to obtain a compact quasi\|monostatic antenna structure for ground penetrating radar systems. And we analyze this transmit/receive dipole pair array in time domain. The numerical results show that if the distance between the transmit antenna and receive antenna is appropriate the array configuration is adoptable.
文摘A time domain electric al field integral equation (TDEFIE) is formulated for the problem of a thin wire antenna in the presence of conductor bodies, and this equation is solved by the method of time marching algorithm. The analysis is valid for any arbitrarily shaped, oriented and positioned wire antennas relative to arbitrarily shaped conductor bodies. Current at the excited point, input admittance and radiation pattern are given and agree with the results computed by the method in frequency domain.