In this paper, we consider the accuracy of integration algorithms such as the implicit Euler and the trapezoidal ones, which are largely employed in the time domain circuit analysis. These algorithms require to make h...
详细信息
In this paper, we consider the accuracy of integration algorithms such as the implicit Euler and the trapezoidal ones, which are largely employed in the time domain circuit analysis. These algorithms require to make hypotheses on the intersample shape and on the "energy content" of the sampled waveforms. For example, the implicit Euler algorithm supposes functions to be piecewise constant. When these hypotheses are violated, some errors are introduced by the integration process into the solution waveform. We consider the energy of the sampled functions, and through energy balance equations, estimate the accuracy of the integration algorithm. Furthermore, we propose an implicit algorithm to determine an adequate integration time step during numerical time domain analysis, This algorithm is based on a global energy balance equation and not on the conventional estimation of the local truncation error. It avoids the "cut and try" mechanism used in SPICE to determine the time step that satisfies the desired error tolerance.
This study investigates the use of explicit/implicit, step/linear space extrapolation as numerical outflow boundary conditions for turbulent flow calculations. The emphasis is on the interplay between the use of curvi...
详细信息
This paper describes a numerical scheme for computing optimal solutions to a class of nonlinear optimal control problems in which parameter uncertainty may be a feature of the state dynamics or objective function. Con...
详细信息
This research relates to a numerical integrator with post-stabilization of several constraints for an autonomous dynamical system. A generally analytical approach shows that the total energy correction is not valid in...
详细信息
This research relates to a numerical integrator with post-stabilization of several constraints for an autonomous dynamical system. A generally analytical approach shows that the total energy correction is not valid in most cases, while post-stabilization of each independent energy is. As a typical test example, we consider a non-integrable Hamiltonian system of three degrees of freedom, which can be split into two independent pieces, one 1D harmonic oscillator and another 2D non-integrable system, by using a transformation of variables. Phase portraits on Poincaxe sections about the 2D system manifest that our analysis is reasonable. In addition, a problem how to compute Lyapunov exponents in constrained systems is proposed. As a suggestion, it is best to stabilize all constraints involving each energy integral and its corresponding variation in order to avoid spurious Lyapunov exponents. Because an appropriately larger time step is acceptable in this sense, it is not expensive to use the fast Lyapunov indicators to describe the global dynamics of phase space for the 3D system, where regions of chaos and order are clearly identified.
The Residual Stress Decomposition Method for Shakedown (RSDM-S) is a new iterative direct method to estimate the shakedown load in a 2-dimensional (2D) loading domain. It may be implemented to any existing finite elem...
详细信息
The Residual Stress Decomposition Method for Shakedown (RSDM-S) is a new iterative direct method to estimate the shakedown load in a 2-dimensional (2D) loading domain. It may be implemented to any existing finite element code, without the need to use a mathematical programming algorithm. An improved and enhanced RSDM-S is proposed herein. A new convergence criterion is presented that makes the procedure almost double as fast. At the same time, the procedure is formulated in a 3 -dimensional (3D) polyhedral loading domain, consisting of independently varying mechanical and thermal loads. Using a cyclic loading program that follows the outline of this domain, it is shown that there is hardly any increase in the computational time when passing from a 2D to a 3D domain. Finally, keeping the efficiency, using an alternative cyclic loading program, an automation of the approach to any n-dimensional loading domain is presented. Examples of application are included. (C) 2017 Elsevier Ltd. All rights reserved.
作者:
Kushner, HJBrown Univ
Dept Appl Math Lefschetz Ctr Dynam Syst Providence RI 02912 USA
Consider the problem of value iteration for solving Markov stochastic games. One simply iterates backward, via a Jacobi-like procedure. The convergence of the Gauss-Seidel form of this procedure is shown for both the ...
详细信息
Consider the problem of value iteration for solving Markov stochastic games. One simply iterates backward, via a Jacobi-like procedure. The convergence of the Gauss-Seidel form of this procedure is shown for both the discounted and ergodic cost problems, under appropriate conditions, with extensions to problems where one stops when a boundary is hit or if any one of the players chooses to stop, with associated costs. Generally, the Gauss-Seidel procedure accelerates convergence.
We develop a finite-difference scheme for approximation of a system of nonlinear PDEs describing the Q-switching process. We construct it by using staggered grids. The transport equations are approximated along charac...
详细信息
We develop a finite-difference scheme for approximation of a system of nonlinear PDEs describing the Q-switching process. We construct it by using staggered grids. The transport equations are approximated along characteristics, and quadratic nonlinear functions are linearized using a special selection of staggered grids. The stability analysis proves that a connection between time and space steps arises only due to approximation requirements in order to follow exactly the directions of characteristics. The convergence analysis of this scheme is done in two steps. First, some estimates of the uniform boundedness of the discrete solution are proved. This part of the analysis is done locally, in some neighborhood of the exact solution. Second, on the basis of the obtained estimates, the main stability inequality is proved. The second-order convergence rate with respect to the space and time coordinates follows from this stability estimate. Using the obtained convergence result, we prove that the local stability analysis in the selected neighborhood of the exact solution is sufficient.
This paper focuses on the problem of computing Lyapunov matrices for integral delay systems. It is shown that these Lyapunov matrices cannot be computed by means of the existing methods for Lyapunov matrices of differ...
详细信息
This letter introduces the NOnSmooth numerical Optimal Control (NOSNOC) open-source software package. It is a modular MATLAB tool based on CasADi and IPOPT for numerically solving Optimal Control Problems (OCP) with p...
详细信息
This letter introduces the NOnSmooth numerical Optimal Control (NOSNOC) open-source software package. It is a modular MATLAB tool based on CasADi and IPOPT for numerically solving Optimal Control Problems (OCP) with piecewise smooth systems (PSS). The tool supports: 1) automatic reformulation of systems with state jumps into PSS (via the time-freezing reformulation [1]) and of PSS into computationally more convenient forms, 2) automatic discretization of the OCP via, e.g., the recently introduced Finite Elements with Switch Detection [2] which enables high accuracy optimal control and simulation of PSS, 3) solution methods for the resulting discrete-time OCP. The nonsmooth discrete-time OCP are solved with techniques of continuous optimization in a homotopy procedure, without the use of integer variables. This enables the treatment of a broad class of nonsmooth systems in a unified way. Two tutorial examples are given. A benchmark shows that NOSNOC provides both faster and more accurate solutions than conventional approaches, including mixed-integer formulations.
We present a novel reformulation of nonsmooth differential equations with state jumps enabling their easier simulation and use in optimal control problems without the need for integer variables. The main idea is to in...
详细信息
We present a novel reformulation of nonsmooth differential equations with state jumps enabling their easier simulation and use in optimal control problems without the need for integer variables. The main idea is to introduce an auxiliary differential equation to mimic the state jump map. Thereby, a clock state is introduced which does not evolve during the runtime of the auxiliary system. The pieces of the trajectory that correspond to the parts when the clock state was evolving recover the solution of the original system with jumps. Our reformulation results in nonsmooth ordinary differential equations where the discontinuity is in the first time derivative of the trajectory, rather than in the trajectory itself. This class of systems is easier to handle both theoretically and numerically. The reformulation is suitable for partially elastic mechanical impact problems. We provide numerical examples demonstrating the ease of use of this reformulation in both simulation and optimal control. In the optimal control example, we solve a sequence of nonlinear programming problems (NLPs) in a homotopy penalization approach and recover a time-optimal trajectory with state jumps.
暂无评论