Modeling time-varying storage components in PSpice

Similar documents
A Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

Application of the homotopy perturbation method to a magneto-elastico-viscous fluid along a semi-infinite plate

SYNCHRONOUS SEQUENTIAL CIRCUITS

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Polynomial Inclusion Functions

THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE

Chapter 9 Method of Weighted Residuals

An inductance lookup table application for analysis of reluctance stepper motor model

4. CONTROL OF ZERO-SEQUENCE CURRENT IN PARALLEL THREE-PHASE CURRENT-BIDIRECTIONAL CONVERTERS

Chapter 3. Modeling with First-Order Differential Equations

ELEC3114 Control Systems 1

STATISTICAL LIKELIHOOD REPRESENTATIONS OF PRIOR KNOWLEDGE IN MACHINE LEARNING

LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION

Conservation laws a simple application to the telegraph equation

How the potentials in different gauges yield the same retarded electric and magnetic fields

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Adjoint Transient Sensitivity Analysis in Circuit Simulation

A PAC-Bayesian Approach to Spectrally-Normalized Margin Bounds for Neural Networks

Chapter 31: RLC Circuits. PHY2049: Chapter 31 1

Linear First-Order Equations

Lecture 6: Control of Three-Phase Inverters

Spurious Significance of Treatment Effects in Overfitted Fixed Effect Models Albrecht Ritschl 1 LSE and CEPR. March 2009

3.0 DETECTION THEORY. Figure 3-1 IF Receiver/Signal Processor Representation

SIMULATION OF POROUS MEDIUM COMBUSTION IN ENGINES

Dynamic Modeling and Analysis of Large-scale Power Systems in the DQ0 Reference Frame

Improving Estimation Accuracy in Nonrandomized Response Questioning Methods by Multiple Answers

Three-Variable Bracket Polynomial for Two-Bridge Knots

Sturm-Liouville Theory

12.11 Laplace s Equation in Cylindrical and

Stable and compact finite difference schemes

Harmonic Modelling of Thyristor Bridges using a Simplified Time Domain Method

Short Intro to Coordinate Transformation

Table of Common Derivatives By David Abraham

TEST 2 (PHY 250) Figure Figure P26.21

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 13. Khyruddin Akbar Ansari, Ph.D., P.E.

Physics 4B. Chapter 31: Questions: 2, 8, 12 Exercises & Problems: 2, 23, 24, 32, 41, 44, 48, 60, 72, 83. n n f

Chapter 6: Energy-Momentum Tensors

. Using a multinomial model gives us the following equation for P d. , with respect to same length term sequences.

Predictive control of synchronous generator: a multiciterial optimization approach

Diagonalization of Matrices Dr. E. Jacobs

AN3400 Application note

Ordinary Differential Equations: Homework 1

Implicit Differentiation

Suppression Method of Rising DC Voltage for the Halt Sequence of an Inverter in the Motor Regeneration

EE 330 Lecture 15. Devices in Semiconductor Processes. Diodes Capacitors MOSFETs

AN INTRODUCTION TO NUMERICAL METHODS USING MATHCAD. Mathcad Release 14. Khyruddin Akbar Ansari, Ph.D., P.E.

28.1 Parametric Yield Estimation Considering Leakage Variability

APPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France

2 GUANGYU LI AND FABIO A. MILNER The coefficient a will be assume to be positive, boune, boune away from zero, an inepenent of t; c will be assume con

State Space Analysis of Power System Stability Enhancement with Used the STATCOM

Implicit Differentiation. Lecture 16.

State-Space Model for a Multi-Machine System

An M/G/1 Retrial Queue with Priority, Balking and Feedback Customers

Reactive Power Compensation in Mechanical Systems

ensembles When working with density operators, we can use this connection to define a generalized Bloch vector: v x Tr x, v y Tr y

The Efficiency Optimization of Permanent Magnet Synchronous Machine DTC for Electric Vehicles Applications Based on Loss Model

Semiclassical analysis of long-wavelength multiphoton processes: The Rydberg atom

Situation awareness of power system based on static voltage security region

Final Exam Study Guide and Practice Problems Solutions

Simple Electromagnetic Motor Model for Torsional Analysis of Variable Speed Drives with an Induction Motor

6 General properties of an autonomous system of two first order ODE

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

Comprehensive Procedure for Fast and Accurate Coupled Oscillator Network Simulation

An Unique SPICE Model of Photodiode with Slowly Changeable Carriers Velocities

Hyperbolic Moment Equations Using Quadrature-Based Projection Methods

arxiv:hep-th/ v1 3 Feb 1993

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

EE 330 Lecture 14. Devices in Semiconductor Processes. Diodes Capacitors MOSFETs

Tractability results for weighted Banach spaces of smooth functions

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

IPA Derivatives for Make-to-Stock Production-Inventory Systems With Backorders Under the (R,r) Policy

Chapter 6. Electromagnetic Oscillations and Alternating Current

inflow outflow Part I. Regular tasks for MAE598/494 Task 1

Math 342 Partial Differential Equations «Viktor Grigoryan

Recommendations: Part 7: Transient Creep for service and accident conditions

Why Bernstein Polynomials Are Better: Fuzzy-Inspired Justification

ANALYSIS OF A GENERAL FAMILY OF REGULARIZED NAVIER-STOKES AND MHD MODELS

Solution to the exam in TFY4230 STATISTICAL PHYSICS Wednesday december 1, 2010

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

Modelling of Three Phase Short Circuit and Measuring Parameters of a Turbo Generator for Improved Performance

Discrete Operators in Canonical Domains

Calculus of Variations

FIRST-ORDER FILTERS GENERALIZED TO THE FRACTIONAL DOMAIN

Finite element analysis of electromagnetic bulging of sheet metals

Efficient Macro-Micro Scale Coupled Modeling of Batteries

THE USE OF KIRCHOFF S CURRENT LAW AND CUT-SET EQUATIONS IN THE ANALYSIS OF BRIDGES AND TRUSSES

Calculus in the AP Physics C Course The Derivative

Statics, Quasistatics, and Transmission Lines

Arm Voltage Estimation Method for Compensated Modulation of Modular Multilevel Converters

Thermal conductivity of graded composites: Numerical simulations and an effective medium approximation

θ x = f ( x,t) could be written as

Lecture 12. Energy, Force, and Work in Electro- and Magneto-Quasistatics

Synchronization of Diffusively Coupled Oscillators: Theory and Experiment

PRACTICE 4. CHARGING AND DISCHARGING A CAPACITOR

EVALUATING HIGHER DERIVATIVE TENSORS BY FORWARD PROPAGATION OF UNIVARIATE TAYLOR SERIES

Chapter 2 Governing Equations

Some vector algebra and the generalized chain rule Ross Bannister Data Assimilation Research Centre, University of Reading, UK Last updated 10/06/10

Thermal runaway during blocking

Design and Application of Fault Current Limiter in Iran Power System Utility

Transcription:

Moeling time-varying storage components in PSpice Dalibor Biolek, Zenek Kolka, Viera Biolkova Dept. of EE, FMT, University of Defence Brno, Czech Republic Dept. of Microelectronics/Raioelectronics, FEEC, Brno University of Technology, Czech Republic fax 973442987- e-mail alibor.biolek@unob.cz- http://user.unob.cz/biolek Abstract The paper escribes a metho of moeling linear non-stationary capacitors an inuctors in PSpice. The capacitance or inuctance is generally varying in time accoring to a law which can be escribe either by an analytical formula or by a set of points, acquire by measurement. This metho is verifie by comparing the simulation results of a sample circuit in OrCa PSpice with the results from Micro-Cap, which enables irect moeling of the time-varying capacitances an inuctances. 1 Introuction The impossibility of a irect simulation of circuits with time-varying capacitors an inuctors belongs to well-known limitations of the PSpice program [1]. This program enables only the moeling of nonlinear polynomial capacitance/voltage an inuctance/current relations [2]. By means of tools of behavioral moeling, general nonlinearities can be moele aitionally on the assumption that the component parameters are time invariant [3]. From this point of view, there is a problem if one nees a transient analysis of parametric circuits in which the storage components vary their parameters accoring to a function of time that is perioical in most cases. This function can be available in the analytical form or as a set of measure points. It is well-known that the common equations v = i, ic = C vc (1) are not true when capacitances an inuctances are time-varying [4]. More complicate equations must be use for such cases: v = ( ( i ) = ( i + i ( (2) ic = ( C( vc ) = C( vc + vc C( However, there are two reasons why the above equations are not a proper starting point for PSpice moeling: 1. They contain the time erivatives of variables an C. While moeling the an C variables versus time via a look-up table, then the piece-wise approximation in PSpice leas to time-iscontinuous erivatives. It can be a source of potential problems uring the simulation incluing the convergence problems. 2. The a C variables as functions of time still appear in the equations. That is why these components cannot be moele by conventional PSpice storage components.

In this paper, a simple proceure is escribe how to overcome both problems, incluing the possibility of efining the initial conitions, i.e. the inuctor current an the capacitor voltage at time. The proceure is vali for linear time-varying an C evices. 2 ( an C( moeling The integration of ifferential equations (2) with subsequent simple arrangement will yiel the equations of time-varying inuctor an capacitor in the integral form: i () i () + t C C = ( v C() v () +, (, vc =, C(. (3a, b) C( It follows from these equations that the inuctor can be moele by a controlle current source an the capacitor by a controlle voltage source. The value of inuctor current at a concrete time instance is given by the initial value of the prouct of current an inuctance, by the integral of voltage, an by the instantaneous value of inuctance. Similarly, the value of capacitor voltage at a concrete time instance is given by the initial value of the prouct of voltage an capacitance, by the integral of current, an by the instantaneous value of capacitance. It is necessary to ensure that neither the capacitance nor the inuctance are zero at any time instance uring the simulation. However, it is a common limitation in the PSpice program. A emonstration of the Spice subcircuit for moeling the inuctor base on Eq. (3a) is in Table 1. The following function is use for moeling the time-varying inuctance: = (1 +.8sin(2π f ), = 1mH, f 5kHz. (4) ( =.subckt inuctor + - params: I=.func (time) {1m*(1+.8*sin(2*pi*5k*time))} gcurr + - value={(s(v(+,-))+i*())/(time)}.ens Table 1: Spice subcircuit for moeling ( accoring to (4). The inuctance is moele via the user-efine function. The inuctor current accoring to (3a) is generate by the G-type current source. The stanar SDT function of PSpice is use to integrate the inuctor voltage. The subcircuit can be calle with the I parameter, which inicates the initial current of the inuctor. A similar subcircuit for moeling the capacitor accoring to Eq. (3b) is in Table 2. For illustration, the time epenence of the capacity is now moele via a look-up table. At time, the capacitance is 5nF, after 1ns it is proportionally increase to 1nF, keeping this value till the time of 1us. Then it falls linearly to 5nF within a time slot of 1ns..subckt capacitor + - params: VC=.func C(time) {TABE(time, +, 5nF, 1ns, 1nF, 1us, 1nF, 1us+1ns, 5nF)} Ec + - value={(s(i(ec))+vc*c())/c(time)}.ens Table 2: Spice subcircuit for moeling the time-varying capacitance by piece-wise-linear function. t i

A rawback of the TABE function consists in the inability of making the time-omain function perioical. This can be overcome by a trick in which the time-varying capacitance or inuctance is moele in PSpice as a signal in the form of voltage or current by an inepenent V- or I- source, associate with the PW attribute. A emonstration is in Table 3..subckt capacitor + - params: VC=.param C 5nF Vc c PW + REPEAT FOREVER + (, {C}) + (1ns, 1nF) + (1us, 1nF) + (1us+1ns, 5nF) + (2us, 5nF) + ENDREPEAT Ec + - value={(s(i(ec))+vc*c)/v(c)}.ens Table 3: Moel of perioically varying capacitance; the look-up table is efine within one repeating perio. 3 Examples of computer simulation et us consier the benchmark RC circuit in Fig. 1. In the first step, the capacitor has a fixe capacitance of C = 7.5nF whereas the inuctance varies in time accoring to formula (4). The transient phenomenon cause by applying the battery is analyze. The initial conitions are as follows: v C () = V an i ()=, 1mA, an 2mA. R out V 1V 1k C Fig. 1: Benchmark circuit with time-varying storage components. The analysis in OrCa PSpice A/D program v. 15.7 using the inuctor subcircuit from Table 1 leas to the results in Fig. 2. The voltage an current waveforms converge to a perioical steay state with a repeating perio of 2µs. It correspons to the frequency of variation of the inuctance values. Since the analytical solution of this circuit is complicate, the simulation results have been compare to those from Micro-Cap program v. 9 [5], in which the moels of time-varying storage components are irectly implemente. In Micro-Cap, two circuits have been analyze simultaneously: the circuit with the utilization of Micro-Cap moels, an the original circuit from PSpice, whose moel bas been exporte into Micro-Cap as a Spice netlist. The results are virtually ientical with small ifferences of less than fractions of percentage points.

4.mA A SE>> -4.mA 1.V I(X:1) V -1.V s 2us 4us 6us 8us 1us V(ou Time Fig. 2: Results of transient analysis of circuit from Fig. 1 for C = 7.5nF an time-varying inuctance accoring to Table 1, for initial conitions i ()=, 1mA, an 2mA. In the next step, the capacitor with C = 7.5nF was replace by a time-varying capacitor accoring to Table 3, when the capacitance is toggle every microsecon between the values of 5nF an 1nF. The resulting waveforms generate by PSpice are in Fig. 3. The agreement with Micro-Cap is again excellent.

4.mA A -4.mA 1.V I(X:1) V SE>> -1.2V s 2us 4us 6us 8us 1us V(ou Time Fig. 3: Results of transient analysis of circuit in Fig. 1 for time-varying inuctance an capacitance accoring to Table 1 an Table 3 an for zero initial conitions. 4 Conclusions An effective metho of moeling time-varying storage components in PSpice is escribe in the paper. The time epenence can be moele by mathematical formulae or via look-up tables, acquire e.g. by measurements. In all the above cases, it is possible to perioize these an C time-omain functions. The metho escribe can take avantage of the transient analysis of a number of parametric circuits whose operation is base on parametric control of storage components.

References [1] PSpice A/D Reference Guie, inclues PSpice A/D, PSpice A/D Basics, an Pspice. Prouct Version 15.7, July 26. In OrCAD_15.7\oc\pspcref\ [2] Vlaimirescu, A. The Spice Book. John Wiley&Sons, Inc., 1994. [3] A Nonlinear Capacitor Moel for Use in Pspice. Application Note, Caence, 1999. http://www.caence.com/appnotes/anonlinearcapacitormoelforuseinpspice.pf [4] Desoer, Ch. A. an Kuh, E.S. Basic Circuit Theory, McGraw-Hill Book Company, 1969. [5] http://www.spectrum-soft.com Acknowlegment This work is supporte by the Grant Agency of the Czech Republic uner grants No. 12/5/771 an 12/5/277, an by the research programmes of BUT MSM216353, MSM2163513, an UD Brno MO FVT43.