Variational Integrators for Electrical Circuits
|
|
- Eugenia Hunt
- 6 years ago
- Views:
Transcription
1 Variational Integrators for Electrical Circuits Sina Ober-Blöbaum California Institute of Technology Joint work with Jerrold E. Marsden, Houman Owhadi, Molei Tao, and Mulin Cheng Structured Integrators Workshop May 8, 29 Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.1
2 Motivation Mechanical systems: variational integrators well established based on discrete variational principle in mechanics Electrical systems: L 1 C 1 R 1 U V C 2 variational principle required Lagrangian formulation with degenerate symplectic form constraints discrete variational scheme Future goal: powerful unified variational scheme for the simulation of electromechanical systems Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.2
3 Outline basics on circuit modeling variational modeling, constraints, and degeneracy approaches to handle constraints and degeneracy construction of integrators standard approach: Modified Nodal Analysis (MNA) and Linear Multistep Methods numerical examples future work Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.3
4 Circuit modeling C 1 L 1 U V C 2 R 1 device linear nonlinear resistor i R = Gu R i R = g(u R, t) capacitor i C = C d dt u C i C = d dt q C (u C, t) inductor u L = L d dt i L u L = d dt φ L(i L, t) device independent controlled voltage source u V = v(t) u V = v(u ctrl, i ctrl, t) current source i I = i(t) i I = i(u ctrl, i ctrl, t) characteristic equations for basic elements + Kirchoff constraint laws Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.4
5 Kirchoff constraint laws Network topology graph consisting of n branches, m + 1 nodes incidence matrix K R n,m describes brach-node connection 1 branch i connected inward to node j K ij = +1 branch i connected outward to node j otherwise K = [K T L K T C K T R K T V K T I ] T, K J R n J,m, J {L, C, R, V, I } Constraints current law (KCL) K T I (t) = I (t): branch currents voltage law (KVL) K u(t) = U(t) u(t): node voltages U(t): branch voltages Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.5
6 Variational formulation mechanical system electrical circuit (linear relations) configuration charge q Q R n velocity current v T q Q R n momentum flux p T q Q R n kinetic energy magnetic energy 1 2 v T Lv R potential energy electrical energy 1 2 qt Cq R friction dissipation R v R n external force external sources E R n with L = diag(l 1,..., L n ), C = diag( 1 C 1,..., 1 C n ), R = diag(r 1,..., R n ), E = (ɛ 1,..., ɛ n ) Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.6
7 Variational formulation Lagrangian L(q, v) = 1 2 v T Lv 1 2 qt Cq on TQ forces f = R v + E T q Q distribution Q (q) = {v T q Q K T v = } = N (K T ) annihilator Q (q) = {w T q Q w, v = v Q (q)} = R(K) Legendre transformation FL(q, v) = (q, L/ v) if FL is not invertible degenerate system Lagrange-d Alembert-Pontryagin Principle on TQ T Q gives implicit Euler-Lagrange equations δ T (L(q, v) + p, q v ) dt + T f δqdt =, δq Q (q) q = v charge conservation (n eq.) ṗ = L q + f + K λ KVL with node voltages λ Rm (n eq.) = L v p flux conservation (n eq.) = K T v KCL (m eq.) Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.7
8 Literature overview Lagrangian / Hamiltonian formulation of circuits MacFarlane (1967): tree structures (capacitor charges or inductor fluxes as generalized coordinates) Chua, McPherson (1973): mixed set of coordinates Chua, McPherson and Milic, Novak: enlarge applicability of Lagrangian methods to electrical networks Degenerate Lagrangian / Hamiltonian Dirac, Bergmann (1956): construction of submanifolds Gotay, Nester (1979): Lagrangian viewpoint van der Schaft (1995): implicit Hamiltonian systems Yoshimura, Marsden (26): implicit Lagrangian systems Leok, Ohsawa (29): discrete version Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.8
9 Reduction integrable distribution Q (q) = {v T q Q K T v = } = N (K T ) configuration submanifold C = {q Q K T (q q ) = } K 2 R n,n m with R(K) R(K 2 ): q = K 2 q + α q Q R n m α = for consistent q v = K 2 ṽ ṽ T q Q R n m constrained Lagrangian L c ( q, ṽ) := L(K 2 q, K 2 ṽ) constrained Legendre transformation FL c ( q, ṽ) = ( q, L c / ṽ) forces f = K T 2 f, momenta p = K T 2 p T q Q R n m reduced L-d A-Pontryagin Principle on T Q T Q δ T ( L c ( q, ṽ) + p, q ṽ ) T dt + f δ qdt = q = ṽ p = Lc q = Lc ṽ + f p Is Lc ṽ invertible? Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.9
10 Cancellation of degeneracy by constraints linear RCLV circuit FL c ( q, ṽ) = ( q, K T 2 LK 2ṽ) constrained Lagrangian L c is non-degenerate if R(K 2 ) N (L) = {} it holds N (L) = {v T q Q v L = } with R(K 2 ) R(K) and R(K) N (K T ) we have R(K 2 ) = N (K T ) = {v T q Q K T L v L +K T C v C +K T R v R +K T V v V = } and thus R(K 2 ) N (L) = {v T q Q [K T C K T R K T V ] [v C v R v V ] T =, v L = } it follows that for N ([K T C K T R K T V ]) = {} the constrained Lagrangian L c is non-degenerate Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.1
11 Examples L 1 L 1 C 1 C2 L 2 C1 C2 C 3 K T L = 1 1 «, K T C = «K T L = 1 «, K T C = N (KC T ) = {} N (K C T ) {} non-degenerate constrained Lagrangian degenerate constrained Lagrangian non-degeneracy iff number of capacitors = number of independent constraints involving capacitors equivalent statements, when resistors and voltage sources are included «Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.11
12 Projection for degenerate systems identify Dirac constraints: denote by FL c M : T Q M the restriction of FL c : T Q T Q to its image primary submanifold M = {( q, p) T Q φ( q, p) = R l } l: degree of degeneracy consistency with dynamics leads to more constraints (motion constrained to lie in M) Dirac-Bergmann: M M 2 M k T Q Gotay-Nester: N N 2 N k T Q knowledge about existence of solutions if constrained algorithm converges consistent equations of motion on N k useful for circuit modeling? Can we find a regular Lagrangian on the final constraint manifold N k? Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.12
13 Example: linear RLCV circuit q = ṽ K T 2 LK 2 ṽ = K T 2 C q K T 2 Rṽ + Ẽ = case 1: no resistors and no external sources φ( q, ṽ, Ẽ) =, linear case: P T q + U T ṽ + V T Ẽ = holonomic constraint defines constraint manifolds C 2 = { q Q P T q = }, T q C 2 = {ṽ T q Q PT ṽ = } reparametrizations q = P 2ˆq, ˆq ˆQ R n m l and ṽ = P 2ˆv, ˆv Tˆq ˆQ R n m l constrained Lagrangian L c2 (ˆq, ˆv) := L(K 2 P 2ˆq, K 2 P 2ˆv) case 2: with resistors and external sources affine non-holonomic constraint defines constraint manifold P 2 = {ṽ T q Q P T q + U T ṽ + V T Ẽ = } reparametrization ṽ = U 2ˆv + γ( q, Ẽ), ˆv R n m l, q R n m constrained Lagrangian L c2 ( q, ˆv) := L(K 2 q, K 2 U 2ˆv + K 2 γ( q, Ẽ)) Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.13
14 The constrained Variational Integrator same reduction steps can be performed on the discrete level to construct a discrete constrained Lagrangian discrete variational principle for discrete constrained Lagrangian case 1: N 1 ( δ L c2 (ˆq k, ˆv k ) + ˆp k, ˆq ) N 1 k ˆq k 1 ˆv k + ˆf k δˆq k dt = h k= with variations δˆq k, δˆv k, δˆp k case 2: N 1 δ k= ( L c2 ( q k, ˆv k ) + with variations δ q k P 2, δˆv k, δ p k k= p k, q ) k q k 1 (U 2ˆv k + γ( q k, Ẽ k ) h gives well-defined scheme without any constraints N 1 + f k δ q k dt = k= Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.14
15 The constrained Variational Integrator straight forward for linear systems loss of physical meaning of generalized coordinates for nonlinear systems the additional constraints are hard to identify no global reparametrization are there other ways to deal with the degeneracy? Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.15
16 Keep degeneracy Variational Principle results into DAE system q = v ṗ = L q + f + K λ = L v p = K T v I I q v p λ = F (q, v, p, λ, f ) for a linear circuit (no external sources) we have I I } {{ } E q v p λ } {{ } ẋ = I C R K L I K T } {{ } A q v p λ } {{ } x Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.16
17 Existence and uniqueness for linear DAE Definition: The pair (E, A) of square matrices is regular, if µ R s.t. det(µe A). Theorem: Let the pair (E, A) of square matrices be regular. Then it holds 1. The DAE system E ẋ(t) = A x(t) is solvable. 2. The set of consistent initial conditions x is nonempty. 3. Every initial value problem with consistent initial conditions is uniquely solvable. (Kunkel, Mehrmann: Differential-Algebraic equations) Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.17
18 Existence and uniqueness for linear DAE Lagrangian system: for µ = 1 we have µe A = E A = I I C R I K L I K T det(e A) N (K T ) N (L + C + R) = {} a linear RLC circuit has a unique solution (it holds N (L + C + R) = {}) a linear RLCV circuit has a unique solution if N (K T V ) = {} Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.18
19 Variational Integrator I discrete Lagrange-d Alembert-Pontryagin Principle N 1 ( δ L(q k, v k ) + p k, q ) N 1 k q k 1 v k + f k δq k dt = h k= k= δq k Q (q k ) q k+1 = q k + hv k+1 p k+1 = p k hcq k hrv k + K λ k = Lv k+1 p k+1 I hi I K L I K T } {{ } E d q k+1 v k+1 p k+1 λ k = K T v k+1 } {{ } x k+1 = I hc hr I } {{ } A d unique discrete solution exists if N (K T ) N (L) = {} integrator applicable if KCL cancels degeneracy q k v k p k λ k 1 } {{ } x k Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.19
20 Variational Integrator II discrete Lagrange-d Alembert-Pontryagin Principle N 1 ( δ L(q k, v k ) + p k, q ) N 1 k+1 q k v k + f k δq k dt = h k= k= δq k Q (q) q k+1 = q k + hv k p k+1 = p k hcq k+1 hrv k+1 + K λ k+1 = Lv k+1 p k+1 I hc hr I K L I K T } {{ } E d = K T v k+1 q k+1 v k+1 p k+1 = λ k+1 } {{ } x k+1 I hi I } {{ } A d unique discrete solution exists if N (K T ) N (L + hr) = {} q k v k p k λ k } {{ } x k Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.2
21 Variational Integrator III higher order discrete Lagrange-d Alembert-Pontryagin Principle of second order resulting into implicit midpoint rule 1 q k+1 = q k + hv k+ 1 2 p k+1 = p k hc q k +q k+1 2 hr v k +v k K λ k+ 1 2 = Lv k+ 1 p k +p k = K T v k+ 1 2 I hi q k+1 h 2 C hr I K I L 2 v k+ 1 2 p k+1 = K T λ k+ 1 2 }{{}}{{} E d x k+1 I h 2 C I I 2 } {{ } A d unique solution exists if N (K T ) N (L + hc + hr) = {} integrator applicable if continuous system has unique solution 1 cf. Nawaf s thesis: Variational Partitioned RK method q k v k 1 2 p k λ k 1 2 }{{} x k Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.21
22 Nonlinear case e.g. Variational Integrator II q k+1 + q k + hv k p = k+1 + p k + hd 1 L(q k+1, v k+1 ) + f k+1 + K λ k+1 D 2 L(q k+1, v k+1 ) p k+1 K T v k+1 = F (x k, x k+1 ) find x k+1 such that F (x k, x k+1 ) = for given x k Jacobian of F w.r.t. x k+1 has to be regular (e.g. Newton scheme) applicability dependent on update rule trade off: the more implicit the better but the more expensive what is state of the art (e.g. simulation software SPICE)? Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.22
23 Modified Nodal Analysis (MNA) (charge-flux) apply KCL to every node except ground insert representation for the branch current of resistors, capacitors and current sources add representation for inductors and voltage sources explicitely to the system KC T C q(k C u, t) + KR T g(k Ru, t) + KL T v L + KV T v V + KI T v I (Ku, q, v L, v V, t) = φ(v L, t) K L u = linear relations v V (Ku, q C, v L, v V, t) K V u = K T C CK C u + K T R GK Ru + K T L v L + K T V v V + K T I v I (Ku, q, v L, v V, t) = L v L K L u = v V (Ku, q C, v L, v V, t) K V u = DAE for node voltages u, inductor currents v L Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.23
24 Numerical integration schemes DAE system A[d(x(t), t)] + b(x(t), t) = conventional approach: Implicit linear multi-step formulas k α i d(x n+1 i, t n+1 i ) = h i= k β i d n+1 i BDF method: β 1 = = β k =, α = 1 k β α 1 α 2 scheme d(x n+1 ) d(x n ) = hd n+1 (implicit Euler) d(x n+1 ) 4 3 d(x n) d(x n 1) = h 2 3 d n+1 advantage: low computational cost (1 function evaluation per step) compared to RK methods disadvantage: bad stability properties (Voigt: General linear methods for integrated circuit design, PhD thesis, 26) i= Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.24
25 Example I C 1 C2 L 1 L = L 1 L 2 L C = B 1 C 1 1 C 2 1 C A, K T L = C A, K T C = ««KCL: [KL T KC T ] [v L v C ] T = N (KC T ) = {} non-degenerate constrained Lagrangian N (KC T ) = {} all variational integrators are applicable Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.25
26 Example I results q C1 q C2 f L1 f L q [C] f [A] t [sec] t [sec] charges on capacitors currents on inductors u 1 u 2 1 x KCL charge conservation flux conservation 1 6 u [V] conservation properties t [sec] t [sec] voltages on nodes conserved quantities Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.26
27 E E Example I results VI/FD VI/BD VI/MP MNA/BDF VI/FD VI/BD VI/MP MNA/BDF t [sec] t [sec] Energy plots: Comparison with BDF1 (left) and BDF2 (right) (h =.1) Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.27
28 E E Example I results VI/FD VI/BD VI/MP MNA/BDF VI/FD VI/BD VI/MP MNA/BDF t [sec] t [sec] Energy plots: Comparison with BDF1 (left) and BDF2 (right) (h =.1) Hairer, Lubich: multi-step methods are in general not symplectic (Tang (1993): The underlying one-step method of a linear multistep method cannot be symplectic.) same discrete solutions for full, reduced and projected system, BUT condition numbers full system reduced system projected system O(h 2 ) O(h 1 ) a Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.28
29 Example II L 1 C 1 R 1 L L 1 1 «A 1, KL T = 1 «C 1 A 1 C 1, KC T = 1 R 1 «A, KR T = 1 R 1 N ([K T C K T R )] = {} non-degenerate constrained Lagrangian constrained Lagrangian L = 1 2 L 1v C 1, dissipation R 1 v harmonic oscillator q 2 Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.29
30 E Example II results 1.8 q C v L q [C] v [A] t [sec] t [sec] charge on capacitor current on inductor u 1 u VI/FD VI/BD VI/MP MNA/BDF1 MNA/BDF u [V] t [sec] t [sec] voltages on nodes energy decay (h=.2) Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.3
31 Example III L 1 C 1 R 1 U V C 2 voltage source U V (t) = sin t resistor R 1 =.1 KL T 1 1 A, KC T A KR T 1 A, KV T 1 1 A 1 1 N ([KC T K R T K V T )] {} degenerate constrained Lagrangian and VI I not applicable N ([K T C K T V )] {} VI II not applicable N (KV T ) = {} VI III applicable Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.31
32 Example III results 1 1 q C v L E q [C] v [A] t [sec] t [sec] charges on capacitors currents on inductor 1 u u 2 u u [V] t [sec] 1 5 VI/MP and MNA/BDF (h=.1) VI/MP (h=.1) MNA/BDF2 (h=.1) MNA/BDF1 (h=.1) t [sec] voltages on nodes energy Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.32
33 Conclusions discrete variational approach provides symplectic integrator with good energy behavior different approaches to treat degeneracy: reduction to final constraint manifold + minimal number of variables (reduces cpu time for large circuits) + condition number of scheme independent on step size - difficult for nonlinear systems - no real physical meaning of generalized coordinates keep degeneracy + applicable to nonlinear systems + easily scalable + full information about charges, currents, node voltages - huge number of variables - condition becomes worse for small step sizes Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.33
34 Future work Open questions: determine degree of degeneracy dependent on circuit topology relation between degree of degeneracy and index of DAE systems Construction of integrators: can we construct multistep variational integrators that are symplectic to be computational more efficient? better way to find generalized coordinates, e.g. exploit graph structure (Todd) identify generate and degenerate sub-circuits and apply different integrators parallelize simulation of sub-circuits Future applications: application to noisy circuits (work in progress with Houman and Molei) application to electromechanical systems Sina Ober-Blöbaum Variational Integrators for Electrical Circuits p.34
Variational integrators for electric and nonsmooth systems
Variational integrators for electric and nonsmooth systems Sina Ober-Blöbaum Control Group Department of Engineering Science University of Oxford Fellow of Harris Manchester College Summerschool Applied
More informationDiscrete Dirac Mechanics and Discrete Dirac Geometry
Discrete Dirac Mechanics and Discrete Dirac Geometry Melvin Leok Mathematics, University of California, San Diego Joint work with Anthony Bloch and Tomoki Ohsawa Geometric Numerical Integration Workshop,
More informationarxiv: v1 [math.na] 8 Mar 2011
Variational integrators for electric circuits Sina Ober-Blöbaum a, Molei Tao b, Mulin Ceng c, Houman Owadi b,c, Jerrold E. Marsden b,c, a Computational Dynamics and Optimal Control, University of Paderborn,
More informationDirac Structures and the Legendre Transformation for Implicit Lagrangian and Hamiltonian Systems
Dirac Structures and the Legendre Transformation for Implicit Lagrangian and Hamiltonian Systems Hiroaki Yoshimura Mechanical Engineering, Waseda University Tokyo, Japan Joint Work with Jerrold E. Marsden
More informationA geometric Birkhoffian formalism for nonlinear RLC networks
Journal of Geometry and Physics 56 (2006) 2545 2572 www.elsevier.com/locate/jgp A geometric Birkhoffian formalism for nonlinear RLC networks Delia Ionescu Institute of Mathematics, Romanian Academy of
More informationarxiv:math/ v1 [math.ds] 5 Sep 2006
arxiv:math/0609153v1 math.ds 5 Sep 2006 A geometric Birkhoffian formalism for nonlinear RLC networks Delia Ionescu, Institute of Mathematics of the Romanian Academy P.O. Box 1-764, RO-014700, Bucharest,
More informationAutomatic Formulation of Circuit Equations
ECE 570 Session 3 IC 752-E Computer Aided Engineering for Integrated Circuits Automatic Formulation of Circuit Equations Objective: Basics of computer aided analysis/simulation Outline:. Discussion of
More informationStochastic Variational Integrators
Stochastic Variational Integrators Jeremy Schmitt Variational integrators are numerical geometric integrator s derived from discretizing Hamilton s principle. They are symplectic integrators that exhibit
More informationDirac Structures in Lagrangian Mechanics
Dirac Structures in Lagrangian Mechanics Part II: Variational Structures Hiroaki Yoshimura Jerrold E. Marsden Department of Mechanical Engineering Control and Dynamical Systems Waseda University California
More informationVariational principles and Hamiltonian Mechanics
A Primer on Geometric Mechanics Variational principles and Hamiltonian Mechanics Alex L. Castro, PUC Rio de Janeiro Henry O. Jacobs, CMS, Caltech Christian Lessig, CMS, Caltech Alex L. Castro (PUC-Rio)
More informationarxiv: v1 [math-ph] 12 Apr 2017
Dirac structures in nonequilibrium thermodynamics arxiv:1704.03935v1 [math-ph] 12 Apr 2017 François Gay-Balmaz Hiroaki Yoshimura CNRS, LMD, IPSL School of Science and Engineering Ecole Normale Supérieure
More informationGeometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics
Geometric Mechanics and Global Nonlinear Control for Multi-Body Dynamics Harris McClamroch Aerospace Engineering, University of Michigan Joint work with Taeyoung Lee (George Washington University) Melvin
More informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
More informationDiscrete Dirac Structures and Implicit Discrete Lagrangian and Hamiltonian Systems
Discrete Dirac Structures and Implicit Discrete Lagrangian and Hamiltonian Systems Melvin Leok and Tomoki Ohsawa Department of Mathematics, University of California, San Diego, La Jolla, CA 92093-0112,
More informationDiscrete Mechanics and Optimal Control for Constrained Systems (DMOCC)
Discrete Mechanics and Optimal Control for Constrained Systems (DMOCC) Sigrid Leyendecker Aeronautics and Control and Dynamical Systems California Institute of Technology 5-45 Pasadena, CA 9115 email:
More informationThe Geometry Underlying Port-Hamiltonian Systems
The Geometry Underlying Port-Hamiltonian Systems Pre-LHMNC School, UTFSM Valparaiso, April 30 - May 1, 2018 Arjan van der Schaft Jan C. Willems Center for Systems and Control Johann Bernoulli Institute
More informationHamilton Jacobi Theory for Degenerate Lagrangian Systems with Holonomic and Nonholonomic Constraints. and Diana Sosa
Hamilton Jacobi Theory for Degenerate Lagrangian Systems with Holonomic and Nonholonomic Constraints Melvin Leok, 1, a) Tomoki Ohsawa, 1, b) 2, c) and Diana Sosa 1) Department of Mathematics, University
More informationModule 24: Outline. Expt. 8: Part 2:Undriven RLC Circuits
Module 24: Undriven RLC Circuits 1 Module 24: Outline Undriven RLC Circuits Expt. 8: Part 2:Undriven RLC Circuits 2 Circuits that Oscillate (LRC) 3 Mass on a Spring: Simple Harmonic Motion (Demonstration)
More informationNetwork Topology-2 & Dual and Duality Choice of independent branch currents and voltages: The solution of a network involves solving of all branch currents and voltages. We know that the branch current
More informationModel order reduction of electrical circuits with nonlinear elements
Model order reduction of electrical circuits with nonlinear elements Andreas Steinbrecher and Tatjana Stykel 1 Introduction The efficient and robust numerical simulation of electrical circuits plays a
More informationQuantize electrical circuits
Quantize electrical circuits a lecture in Quantum Informatics the 4th and 7th of September 017 Thilo Bauch and Göran Johansson In this lecture we will discuss how to derive the quantum mechanical hamiltonian
More informationCHAPTER 10: Numerical Methods for DAEs
CHAPTER 10: Numerical Methods for DAEs Numerical approaches for the solution of DAEs divide roughly into two classes: 1. direct discretization 2. reformulation (index reduction) plus discretization Direct
More informationPort-Hamiltonian systems: network modeling and control of nonlinear physical systems
Port-Hamiltonian systems: network modeling and control of nonlinear physical systems A.J. van der Schaft February 3, 2004 Abstract It is shown how port-based modeling of lumped-parameter complex physical
More informationElectric Circuit Theory
Electric Circuit Theory Nam Ki Min nkmin@korea.ac.kr 010-9419-2320 Chapter 8 Natural and Step Responses of RLC Circuits Nam Ki Min nkmin@korea.ac.kr 010-9419-2320 8.1 Introduction to the Natural Response
More informationReducing round-off errors in symmetric multistep methods
Reducing round-off errors in symmetric multistep methods Paola Console a, Ernst Hairer a a Section de Mathématiques, Université de Genève, 2-4 rue du Lièvre, CH-1211 Genève 4, Switzerland. (Paola.Console@unige.ch,
More informationGeneralized Coordinates, Lagrangians
Generalized Coordinates, Lagrangians Sourendu Gupta TIFR, Mumbai, India Classical Mechanics 2012 August 10, 2012 Generalized coordinates Consider again the motion of a simple pendulum. Since it is one
More informationModeling of Electromechanical Systems
Page 1 of 54 Modeling of Electromechanical Systems Werner Haas, Kurt Schlacher and Reinhard Gahleitner Johannes Kepler University Linz, Department of Automatic Control, Altenbergerstr.69, A 4040 Linz,
More informationVariational Collision Integrators and Optimal Control
Variational Collision Integrators and Optimal Control David Pekarek, and Jerrold Marsden 1 Abstract This paper presents a methodology for generating locally optimal control policies for mechanical systems
More informationThe Harmonic Balance Method
For Nonlinear Microwave Circuits Hans-Dieter Lang, Xingqi Zhang Thursday, April 25, 2013 ECE 1254 Modeling of Multiphysics Systems Course Project Presentation University of Toronto Contents Balancing the
More informationPort-Hamiltonian systems: a theory for modeling, simulation and control of complex physical systems
Port-Hamiltonian systems: a theory for modeling, simulation and control of complex physical systems A.J. van der Schaft B.M. Maschke July 2, 2003 Abstract It is shown how port-based modeling of lumped-parameter
More informationPort-Hamiltonian Systems: from Geometric Network Modeling to Control
Port-Hamiltonian Systems: from Geometric Network Modeling to Control, EECI, April, 2009 1 Port-Hamiltonian Systems: from Geometric Network Modeling to Control, EECI, April, 2009 2 Port-Hamiltonian Systems:
More informationModified Equations for Variational Integrators
Modified Equations for Variational Integrators Mats Vermeeren Technische Universität Berlin Groningen December 18, 2018 Mats Vermeeren (TU Berlin) Modified equations for variational integrators December
More informationSuperconductors: Quantum circuits
Superconductors: Quantum circuits J. J. García-Ripoll IFF, CSIC Madrid (20-4-2009) Mesoscopic QIPC Small systems So far we have only seen small systems to store and process QI Individual atoms As trapped
More informationSolution: Based on the slope of q(t): 20 A for 0 t 1 s dt = 0 for 3 t 4 s. 20 A for 4 t 5 s 0 for t 5 s 20 C. t (s) 20 C. i (A) Fig. P1.
Problem 1.24 The plot in Fig. P1.24 displays the cumulative charge q(t) that has entered a certain device up to time t. Sketch a plot of the corresponding current i(t). q 20 C 0 1 2 3 4 5 t (s) 20 C Figure
More informationLagrange-d Alembert integrators for constrained systems in mechanics
Lagrange-d Alembert integrators for constrained systems in mechanics Conference on Scientific Computing in honor of Ernst Hairer s 6th birthday, Geneva, Switzerland Laurent O. Jay Dept. of Mathematics,
More informationElectrical Circuits I
Electrical Circuits I This lecture discusses the mathematical modeling of simple electrical linear circuits. When modeling a circuit, one ends up with a set of implicitly formulated algebraic and differential
More informationBasics of Network Theory (Part-I)
Basics of Network Theory (PartI). A square waveform as shown in figure is applied across mh ideal inductor. The current through the inductor is a. wave of peak amplitude. V 0 0.5 t (m sec) [Gate 987: Marks]
More informationQUESTION BANK SUBJECT: NETWORK ANALYSIS (10ES34)
QUESTION BANK SUBJECT: NETWORK ANALYSIS (10ES34) NOTE: FOR NUMERICAL PROBLEMS FOR ALL UNITS EXCEPT UNIT 5 REFER THE E-BOOK ENGINEERING CIRCUIT ANALYSIS, 7 th EDITION HAYT AND KIMMERLY. PAGE NUMBERS OF
More informationTWO APPROACHES UTILIZING INVARIANT MANIFOLDS TO DESIGN TRAJECTORIES FOR DMOC OPTIMIZATION
TWO APPROACHES UTILIZING INVARIANT MANIFOLDS TO DESIGN TRAJECTORIES FOR DMOC OPTIMIZATION INTRODUCTION Ashley Moore Invariant manifolds of the planar circular restricted 3-body problem are used to design
More informationLecture I: Constrained Hamiltonian systems
Lecture I: Constrained Hamiltonian systems (Courses in canonical gravity) Yaser Tavakoli December 15, 2014 1 Introduction In canonical formulation of general relativity, geometry of space-time is given
More informationUNIVERSITY OF CALIFORNIA, SAN DIEGO
UNIVERSITY OF CALIFORNIA, SAN DIEGO Structured approaches to large-scale systems: Variational integrators for interconnected Lagrange-Dirac systems and structured model reduction on Lie groups A dissertation
More informationarxiv: v2 [math.na] 15 Apr 2016
Zeitschrift für Angewandte Mathematik und Mechanik, 8/4/206 C -continuous space-time discretization based on Hamilton s law of varying action Janine C. Mergel, Roger A. Sauer,, and Sina Ober-Blöbaum 2
More informationModeling of electromechanical systems
Modeling of electromechanical systems Carles Batlle II EURON/GEOPLEX Summer School on Modeling and Control of Complex Dynamical Systems Bertinoro, Italy, July 18-22 2005 Contents of this lecture Basic
More informationENGR 2405 Chapter 8. Second Order Circuits
ENGR 2405 Chapter 8 Second Order Circuits Overview The previous chapter introduced the concept of first order circuits. This chapter will expand on that with second order circuits: those that need a second
More informationNetwork Graphs and Tellegen s Theorem
Networ Graphs and Tellegen s Theorem The concepts of a graph Cut sets and Kirchhoff s current laws Loops and Kirchhoff s voltage laws Tellegen s Theorem The concepts of a graph The analysis of a complex
More informationarxiv: v1 [math.oc] 1 Feb 2015
HIGH ORDER VARIATIONAL INTEGRATORS IN THE OPTIMAL CONTROL OF MECHANICAL SYSTEMS CÉDRIC M. CAMPOS, SINA OBER-BLÖBAUM, AND EMMANUEL TRÉLAT arxiv:152.325v1 [math.oc] 1 Feb 215 Abstract. In recent years, much
More informationRobotics. Dynamics. Marc Toussaint U Stuttgart
Robotics Dynamics 1D point mass, damping & oscillation, PID, dynamics of mechanical systems, Euler-Lagrange equation, Newton-Euler recursion, general robot dynamics, joint space control, reference trajectory
More informationREVIEW. Hamilton s principle. based on FW-18. Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws!
Hamilton s principle Variational statement of mechanics: (for conservative forces) action Equivalent to Newton s laws! based on FW-18 REVIEW the particle takes the path that minimizes the integrated difference
More informationP441 Analytical Mechanics - I. RLC Circuits. c Alex R. Dzierba. In this note we discuss electrical oscillating circuits: undamped, damped and driven.
Lecture 10 Monday - September 19, 005 Written or last updated: September 19, 005 P441 Analytical Mechanics - I RLC Circuits c Alex R. Dzierba Introduction In this note we discuss electrical oscillating
More informationDIRAC COTANGENT BUNDLE REDUCTION HIROAKI YOSHIMURA JERROLD E. MARSDEN. (Communicated by Juan-Pablo Ortega)
JOURNAL OF GEOMETRIC MECHANICS doi:10.3934/jgm.2009.1.87 c American Institute of Mathematical Sciences Volume 1, Number 1, March 2009 pp. 87 158 DIRAC COTANGENT BUNDLE REDUCTION HIROAKI YOSHIMURA Applied
More informationDiscretizations of Lagrangian Mechanics
Discretizations of Lagrangian Mechanics George W. Patrick Applied Mathematics and Mathematical Physics Mathematics and Statistics, University of Saskatchewan, Canada August 2007 With Charles Cuell Discretizations
More informationAPPLICATION TO TRANSIENT ANALYSIS OF ELECTRICAL CIRCUITS
EECE 552 Numerical Circuit Analysis Chapter Nine APPLICATION TO TRANSIENT ANALYSIS OF ELECTRICAL CIRCUITS I. Hajj Application to Electrical Circuits Method 1: Construct state equations = f(x, t) Method
More information9. Introduction and Chapter Objectives
Real Analog - Circuits 1 Chapter 9: Introduction to State Variable Models 9. Introduction and Chapter Objectives In our analysis approach of dynamic systems so far, we have defined variables which describe
More informationWhen Gradient Systems and Hamiltonian Systems Meet
When Gradient Systems and Hamiltonian Systems Meet Arjan van der Schaft Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, the Netherlands December 11, 2011 on the
More informationEE292: Fundamentals of ECE
EE292: Fundamentals of ECE Fall 2012 TTh 10:00-11:15 SEB 1242 Lecture 14 121011 http://www.ee.unlv.edu/~b1morris/ee292/ 2 Outline Review Steady-State Analysis RC Circuits RL Circuits 3 DC Steady-State
More informationSource-Free RC Circuit
First Order Circuits Source-Free RC Circuit Initial charge on capacitor q = Cv(0) so that voltage at time 0 is v(0). What is v(t)? Prof Carruthers (ECE @ BU) EK307 Notes Summer 2018 150 / 264 First Order
More informationDISCRETE MECHANICS AND OPTIMAL CONTROL
DISCRETE MECHANICS AND OPTIMAL CONTROL Oliver Junge, Jerrold E. Marsden, Sina Ober-Blöbaum Institute for Mathematics, University of Paderborn, Warburger Str. 1, 3398 Paderborn, Germany Control and Dynamical
More informationHandout 10: Inductance. Self-Inductance and inductors
1 Handout 10: Inductance Self-Inductance and inductors In Fig. 1, electric current is present in an isolate circuit, setting up magnetic field that causes a magnetic flux through the circuit itself. This
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationLecture 3 BRANCHES AND NODES
Lecture 3 Definitions: Circuits, Nodes, Branches Kirchoff s Voltage Law (KVL) Kirchoff s Current Law (KCL) Examples and generalizations RC Circuit Solution 1 Branch: BRANCHES AND NODES elements connected
More informationSeries & Parallel Resistors 3/17/2015 1
Series & Parallel Resistors 3/17/2015 1 Series Resistors & Voltage Division Consider the single-loop circuit as shown in figure. The two resistors are in series, since the same current i flows in both
More informationHAMILTON S PRINCIPLE
HAMILTON S PRINCIPLE In our previous derivation of Lagrange s equations we started from the Newtonian vector equations of motion and via D Alembert s Principle changed coordinates to generalised coordinates
More informationCanonical transformations (Lecture 4)
Canonical transformations (Lecture 4) January 26, 2016 61/441 Lecture outline We will introduce and discuss canonical transformations that conserve the Hamiltonian structure of equations of motion. Poisson
More informationSPECTRAL VARIATIONAL INTEGRATORS
SPECTRAL VARIATIONAL INTEGRATORS JAMES HALL AND MELVIN LEOK Abstract. In this paper, we present a new variational integrator for problems in Lagrangian mechanics. Using techniques from Galerkin variational
More informationChapter 1W Basic Electromagnetic Concepts
Chapter 1W Basic Electromagnetic Concepts 1W Basic Electromagnetic Concepts 1W.1 Examples and Problems on Electric Circuits 1W.2 Examples on Magnetic Concepts This chapter includes additional examples
More informationProblem 1: Lagrangians and Conserved Quantities. Consider the following action for a particle of mass m moving in one dimension
105A Practice Final Solutions March 13, 01 William Kelly Problem 1: Lagrangians and Conserved Quantities Consider the following action for a particle of mass m moving in one dimension S = dtl = mc dt 1
More informationMath 692A: Geometric Numerical Integration
Math 692A: Geometric Numerical Integration Melvin Leok Mathematics, Purdue University. mleok@math.purdue.edu http://www.math.purdue.edu/ mleok/ Mathematics, Purdue University NSF DMS-54747, 714223 and
More informationRunge-Kutta type methods for differentialalgebraic equations in mechanics
University of Iowa Iowa Research Online Theses and Dissertations Spring 211 Runge-Kutta type methods for differentialalgebraic equations in mechanics Scott Joseph Small University of Iowa Copyright 211
More informationPreamble. Circuit Analysis II. Mesh Analysis. When circuits get really complex methods learned so far will still work,
Preamble Circuit Analysis II Physics, 8 th Edition Custom Edition Cutnell & Johnson When circuits get really complex methods learned so far will still work, but they can take a long time to do. A particularly
More informationModel reduction of nonlinear circuit equations
Model reduction of nonlinear circuit equations Tatjana Stykel Technische Universität Berlin Joint work with T. Reis and A. Steinbrecher BIRS Workshop, Banff, Canada, October 25-29, 2010 T. Stykel. Model
More informationDefect-based a-posteriori error estimation for implicit ODEs and DAEs
1 / 24 Defect-based a-posteriori error estimation for implicit ODEs and DAEs W. Auzinger Institute for Analysis and Scientific Computing Vienna University of Technology Workshop on Innovative Integrators
More informationDynamics. Basilio Bona. Semester 1, DAUIN Politecnico di Torino. B. Bona (DAUIN) Dynamics Semester 1, / 18
Dynamics Basilio Bona DAUIN Politecnico di Torino Semester 1, 2016-17 B. Bona (DAUIN) Dynamics Semester 1, 2016-17 1 / 18 Dynamics Dynamics studies the relations between the 3D space generalized forces
More informationClassical Mechanics and Electrodynamics
Classical Mechanics and Electrodynamics Lecture notes FYS 3120 Jon Magne Leinaas Department of Physics, University of Oslo December 2009 2 Preface These notes are prepared for the physics course FYS 3120,
More informationVariational Integrators for Interconnected Lagrange Dirac Systems
DOI 10.1007/s00332-017-9364-7 Variational Integrators for Interconnected Lagrange Dirac Systems Helen Parks 1 Melvin Leok 1 Received: 4 March 2016 / Accepted: 28 January 2017 Springer Science+Business
More informationNoether Symmetries and Conserved Momenta of Dirac Equation in Presymplectic Dynamics
International Mathematical Forum, 2, 2007, no. 45, 2207-2220 Noether Symmetries and Conserved Momenta of Dirac Equation in Presymplectic Dynamics Renato Grassini Department of Mathematics and Applications
More informationFactorization of Indefinite Systems Associated with RLC Circuits. Patricio Rosen Wil Schilders, Joost Rommes
Factorization of Indefinite Systems Associated with RLC Circuits Patricio Rosen Wil Schilders, Joost Rommes Outline Circuit Equations Solution Methods (Schilders Factorization Incidence Matrix Decomposition
More informationDelhi Noida Bhopal Hyderabad Jaipur Lucknow Indore Pune Bhubaneswar Kolkata Patna Web: Ph:
Serial : CH_EE_B_Network Theory_098 Delhi Noida Bhopal Hyderabad Jaipur Lucknow ndore Pune Bhubaneswar Kolkata Patna Web: E-mail: info@madeeasy.in Ph: 0-56 CLASS TEST 08-9 ELECTCAL ENGNEENG Subject : Network
More informationElectric Circuits I. Nodal Analysis. Dr. Firas Obeidat
Electric Circuits I Nodal Analysis Dr. Firas Obeidat 1 Nodal Analysis Without Voltage Source Nodal analysis, which is based on a systematic application of Kirchhoff s current law (KCL). A node is defined
More information15. Hamiltonian Mechanics
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 15. Hamiltonian Mechanics Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationElectromagnetic Oscillations and Alternating Current. 1. Electromagnetic oscillations and LC circuit 2. Alternating Current 3.
Electromagnetic Oscillations and Alternating Current 1. Electromagnetic oscillations and LC circuit 2. Alternating Current 3. RLC circuit in AC 1 RL and RC circuits RL RC Charging Discharging I = emf R
More informationKirchhoff's Laws and Circuit Analysis (EC 2)
Kirchhoff's Laws and Circuit Analysis (EC ) Circuit analysis: solving for I and V at each element Linear circuits: involve resistors, capacitors, inductors Initial analysis uses only resistors Power sources,
More information2.004 Dynamics and Control II Spring 2008
MT OpenCourseWare http://ocwmitedu 200 Dynamics and Control Spring 200 For information about citing these materials or our Terms of Use, visit: http://ocwmitedu/terms Massachusetts nstitute of Technology
More informationSinusoidal Steady State Analysis (AC Analysis) Part I
Sinusoidal Steady State Analysis (AC Analysis) Part I Amin Electronics and Electrical Communications Engineering Department (EECE) Cairo University elc.n102.eng@gmail.com http://scholar.cu.edu.eg/refky/
More informationNumerische Mathematik
Numer. Math. DOI 1.17/s211-14-679- Numerische Mathematik Spectral variational integrators James Hall Melvin Leok Received: 19 November 212 / Revised: 16 October 214 Springer-Verlag Berlin Heidelberg 214
More informationJ10M.1 - Rod on a Rail (M93M.2)
Part I - Mechanics J10M.1 - Rod on a Rail (M93M.2) J10M.1 - Rod on a Rail (M93M.2) s α l θ g z x A uniform rod of length l and mass m moves in the x-z plane. One end of the rod is suspended from a straight
More informationNotes for course EE1.1 Circuit Analysis TOPIC 10 2-PORT CIRCUITS
Objectives: Introduction Notes for course EE1.1 Circuit Analysis 4-5 Re-examination of 1-port sub-circuits Admittance parameters for -port circuits TOPIC 1 -PORT CIRCUITS Gain and port impedance from -port
More informationIntroduction to Hamiltonian systems
Introduction to Hamiltonian systems Marlis Hochbruck Heinrich-Heine Universität Düsseldorf Oberwolfach Seminar, November 2008 Examples Mathematical biology: Lotka-Volterra model First numerical methods
More informationElectric and Magnetic Forces in Lagrangian and Hamiltonian Formalism
Electric and Magnetic Forces in Lagrangian and Hamiltonian Formalism Benjamin Hornberger 1/26/1 Phy 55, Classical Electrodynamics, Prof. Goldhaber Lecture notes from Oct. 26, 21 Lecture held by Prof. Weisberger
More informationOptimal control problems with PDE constraints
Optimal control problems with PDE constraints Maya Neytcheva CIM, October 2017 General framework Unconstrained optimization problems min f (q) q x R n (real vector) and f : R n R is a smooth function.
More informationEnergy-conserving formulation of RLC-circuits with linear resistors
Energy-conserving formulation of RLC-circuits with linear resistors D. Eberard, B.M. Maschke and A.J. van der Schaft Abstract In this paper firstly, the dynamics of LC-circuits without excess elements
More informationIntroduction to the. Geometry of Classical Dynamics
Introduction to the Geometry of Classical Dynamics Renato Grassini Dipartimento di Matematica e Applicazioni Università di Napoli Federico II HIKARI LT D HIKARI LTD Hikari Ltd is a publisher of international
More informationA Very General Hamilton Jacobi Theorem
University of Salerno 9th AIMS ICDSDEA Orlando (Florida) USA, July 1 5, 2012 Introduction: Generalized Hamilton-Jacobi Theorem Hamilton-Jacobi (HJ) Theorem: Consider a Hamiltonian system with Hamiltonian
More informationAnalysis and Control of Multi-Robot Systems. Elements of Port-Hamiltonian Modeling
Elective in Robotics 2014/2015 Analysis and Control of Multi-Robot Systems Elements of Port-Hamiltonian Modeling Dr. Paolo Robuffo Giordano CNRS, Irisa/Inria! Rennes, France Introduction to Port-Hamiltonian
More informationAlgebraic flux correction and its application to convection-dominated flow. Matthias Möller
Algebraic flux correction and its application to convection-dominated flow Matthias Möller matthias.moeller@math.uni-dortmund.de Institute of Applied Mathematics (LS III) University of Dortmund, Germany
More informationLecture IV: Time Discretization
Lecture IV: Time Discretization Motivation Kinematics: continuous motion in continuous time Computer simulation: Discrete time steps t Discrete Space (mesh particles) Updating Position Force induces acceleration.
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationModeling of Dynamic Systems: Notes on Bond Graphs Version 1.0 Copyright Diane L. Peters, Ph.D., P.E.
Modeling of Dynamic Systems: Notes on Bond Graphs Version 1.0 Copyright 2015 Diane L. Peters, Ph.D., P.E. Spring 2015 2 Contents 1 Overview of Dynamic Modeling 5 2 Bond Graph Basics 7 2.1 Causality.............................
More informationAppendix A Installing QUCS
Appendix A Installing QUCS In this appendix, we will discuss how to install QUCS [1]. Note that QUCS has a lot of components, many of which we will not use. Nevertheless, we will install all components
More informationAn Algorithmic Framework of Large-Scale Circuit Simulation Using Exponential Integrators
An Algorithmic Framework of Large-Scale Circuit Simulation Using Exponential Integrators Hao Zhuang 1, Wenjian Yu 2, Ilgweon Kang 1, Xinan Wang 1, and Chung-Kuan Cheng 1 1. University of California, San
More informationGauge Fixing and Constrained Dynamics in Numerical Relativity
Gauge Fixing and Constrained Dynamics in Numerical Relativity Jon Allen The Dirac formalism for dealing with constraints in a canonical Hamiltonian formulation is reviewed. Gauge freedom is discussed and
More information