The numerical calculation of the limit cycle of oscillators with resonators exhibiting a high-quality factor Q such as quartz crystals is a difficult task in the time domain. Time domain integration formulas, when not...
详细信息
The numerical calculation of the limit cycle of oscillators with resonators exhibiting a high-quality factor Q such as quartz crystals is a difficult task in the time domain. Time domain integration formulas, when not carefully selected, introduce numerical damping that leads to erroneous limit cycles or spurious oscillations. A novel class of adaptive multistep integration formulas based on finite difference (FD) schemes is derived, which circumvent the aforementioned problems. The results are compared with the well-known harmonic balance (HB) technique. Moreover, the range of absolute stability is derived for these methods. The resulting discretized system by FD methods is sparser than that of HB and, therefore, easier to solve and easier to implement.
Fractional Ginzburg-Landau equations as generalizations of the classical one have been used to describe various physical phenomena. In this paper, we propose a numericalintegration method for space fractional Ginzbur...
详细信息
Fractional Ginzburg-Landau equations as generalizations of the classical one have been used to describe various physical phenomena. In this paper, we propose a numericalintegration method for space fractional Ginzburg-Landau equations based on a dynamical low-rank approximation. We first approximate the space fractional derivatives by using a fractional centered difference method. Then, the resulting matrix differential equation is split into a stiff linear part and a nonstiff (nonlinear) one. For solving these two subproblems, a dynamical low-rank approach is employed. The convergence of our method is proved rigorously. numerical examples are reported which show that the proposed method is robust and accurate. (C) 2021 Elsevier Inc. All rights reserved.
Coordinate-free Lie-group integration method of arbitrary (and possibly higher) order of accuracy for constrained multibody systems (MBS) is proposed in the paper. Mathematical model of MBS dynamics is shaped as a DAE...
详细信息
Coordinate-free Lie-group integration method of arbitrary (and possibly higher) order of accuracy for constrained multibody systems (MBS) is proposed in the paper. Mathematical model of MBS dynamics is shaped as a DAE system of equations of index 1, whereas dynamics is evolving on the system state space modeled as a Lie-group. Since the formulated integration algorithm operates directly on the system manifold via MBS elements' angular velocities and rotational matrices, no local rotational coordinates are necessary, and kinematical differential equations (that are prone to singularities in the case of three-parameter-based local description of the rotational kinematics) are completely avoided. Basis of the integration procedure is the Munthe-Kaas algorithm for ODE integration on Lie-groups, which is reformulated and expanded to be applicable for the integration of constrained MBS in the DAE-index-1 form. In order to eliminate numerical constraint violation for generalized positions and velocities during the integration procedure, a constraint stabilization projection method based on constrained least-square minimization algorithm is introduced. Two numerical examples, heavy top dynamics and satellite with mounted 5-DOF manipulator, are presented. The proposed Lie-group DAE-index-1 integration scheme is easy-to-use for an MBS with kinematical constraints of general type, and it is especially suitable for dynamics of mechanical systems with large 3D rotations where standard (vector space) formulations might be inefficient due to kinematical singularities (three-parameter-based rotational coordinates) or additional kinematical constraints (redundant quaternion formulations).
This article presents a new numerical method for integration of ordinary differential equations based on state variable quantization. Using the idea of implicit integration, the new method called BQSS (Backward Quanti...
详细信息
This article presents a new numerical method for integration of ordinary differential equations based on state variable quantization. Using the idea of implicit integration, the new method called BQSS (Backward Quantized State Systems) allows integrating stiff systems in an efficient way. Being the first quantization-based method for stiff systems, BQSS is itself an explicit method, so the contribution is important in the general context of numericalintegration. Besides introducing the method, the article studies its main theoretical properties, discusses some practical issues related to the algorithm implementation and presents simulation results. Copyright (C) 2007 CEA-IFAC.
It is known that Gauss-Raclau quadrature rule integral(1)(-1) f(x)dx similar or equal to Sigma(n)(i=1) a(i)f(b(i)) + pf(-1) (or qf(1)), is exact for polynomials of degree at most 2n. In this paper we intend to find a ...
详细信息
It is known that Gauss-Raclau quadrature rule integral(1)(-1) f(x)dx similar or equal to Sigma(n)(i=1) a(i)f(b(i)) + pf(-1) (or qf(1)), is exact for polynomials of degree at most 2n. In this paper we intend to find a formula which is nearly exact for monomial functions x(j), j = 0, 1,..., 2n + 2, instead of being analytically exact for the basis space x(j), j = 0, 1,..., 2n. In this way, several examples are also given to show the numerical superiority of the presented rules with respect to usual Gauss-Radau quadrature rules. (c) 2004 Elsevier Inc. All rights reserved.
This article presents a family of variational integrators from a continuous time point of view. A general procedure for deriving symplectic integration schemes preserving an energy-like quantity is shown, which is bas...
详细信息
This article presents a family of variational integrators from a continuous time point of view. A general procedure for deriving symplectic integration schemes preserving an energy-like quantity is shown, which is based on the principle of virtual work. The framework is extended to incorporate holonomic constraints without using additional regularization. In addition, it is related to well-known partitioned Runge-Kutta methods and to other variational integration schemes. As an example, a concrete integration scheme is derived for the planar pendulum using both polar and Cartesian coordinates. Copyright (C) 2016 John Wiley & Sons, Ltd.
One of the quadrature rules is the "Equal coefficients quadrature rules" represented by integral(h)(a) w(x)f(x)dx similar or equal to C-''Sigma(n)(i=1)f(x(i)) a where C-n is a constant number and w(x...
详细信息
One of the quadrature rules is the "Equal coefficients quadrature rules" represented by integral(h)(a) w(x)f(x)dx similar or equal to C-''Sigma(n)(i=1)f(x(i)) a where C-n is a constant number and w(x) is a weight function on [a, b]. In this work, we show that the precisian degree of above formula can be increased by taking the upper and lower bounds of the integration formula as unknowns. This causes to numerically be extended the monomial space {1, x,..., x(n)} to {1, x,..., x(n+2).} We use a matrix proof to show that the resulting nonlinear system for the basis f(x) = x(j), j = 0,...,n + 2 has no analytic solution. Thus, we solve this system numerically to find unknowns x(1),x(2),x(n), C-n, a and b. Finally, some examples will be given to show the numerical superiority of the new developed method. (c) 2005 Elsevier Inc. All rights reserved.
This paper discusses on numerical improvement of the Newton-Cotes integration rules, which are in forms of: integral(b=a+nh)(a) f(x) dx similar or equal to Sigma(n)(k=0) B-k((n)) integral(a + kh). It is known that the...
详细信息
This paper discusses on numerical improvement of the Newton-Cotes integration rules, which are in forms of: integral(b=a+nh)(a) f(x) dx similar or equal to Sigma(n)(k=0) B-k((n)) integral(a + kh). It is known that the precision degree of above formula is n + 1 for even n's and is n for odd n's. However, if its bounds are considered as two additional variables (i.e. a and h in fact) we reach a nonlinear system that numerically improves the precision degree of above integration formula up to degree n + 2. In this way, some numerical examples are given to show the numerical superiority of our approach with respect to usual Newton-Cotes integration formulas. (c) 2004 Elsevier Inc. All rights reserved.
The effects of considering the soil-structure interaction on the structural response can be determined using the cone model. In the general approach of analysis using the cone model, the analysis is performed in the f...
详细信息
The effects of considering the soil-structure interaction on the structural response can be determined using the cone model. In the general approach of analysis using the cone model, the analysis is performed in the frequency domain;then, the results are reverted to the time domain through a Fourier transformation. This indirect procedure is only applicable for linear systems;whereas, for considering the nonlinear behavior of the system, the analysis should be performed in the time domain. Thus, in this study, for considering the effects of interaction, cone model is used directly in the time domain. Formulations of the motion equations in the cone model are provided directly in the time domain for two types of superstructures, including rigid block and one degree-of-freedom superstructure. In order to demonstrate the effects of soil layering on the formulation of cone models, the sub-structure half-space is considered as a layer on the rigid and flexible rocks. In order to determine the discrepancy between the results of the time and frequency domains, the results of analyses in two domains are compared. Based on the outcomes of this research, the concurrency of the system responses, both in the direct method of analysis in the time domain and the indirect method of analysis, is ascertained. Thus, the direct method of analysis in the time domain is applicable for manifold soil layering.
In many memristor models the internal variable is related to the width of a conductive layer and is therefore limited by the physical dimensions of the device. The necessity of keeping the internal variable within all...
详细信息
In many memristor models the internal variable is related to the width of a conductive layer and is therefore limited by the physical dimensions of the device. The necessity of keeping the internal variable within allowable limits may introduce discontinuities in right hand sides of differential equations describing dynamics of circuits with memristors. Such discontinuities are difficult to handle using standard numerical integration methods. An approach based on the concept of a Poincare map is proposed to solve these difficulties. Two examples are discussed to show the usefulness of the proposed technique.
暂无评论