|
Easily Add Sweep and Optimize Part
Parameters — Within Sweep and Optimize dialogs you can now select
from a part list and edit part parameters without leaving the dialog
box.
|
Better Optimizer Information — After
you enter a percent tolerance for the optimized parameter, you
are able to view the minimum, maximum and mean value for that
parameter.
|
|
Modified B Element — Added the capability to
copy “B Element” expressions to the part comments, so they will be displayed
on the drawing. An Enter dialog was added for ease of entering very
long expressions.
New ICL Scripts Added — Expandx, expandy,
and expandxAllPlots ICL scripts were added to let you easily modify
your waveform viewing area. Family script was added to create eye diagrams.
Superfamily script combines multiple vector tags into one to reduce
presentation clutter. Drawing scripts allow you to easily format waveform
graphs.
Easy Location of Subcircuit Vectors — The Add Waveform Dialog
in IntuScope now features an expandable tree structure for subcircuit
vectors.
Dynamic Help Menu For Scope5 —
HELP now includes updated Movie Tutorials, and web access to the
latest documentation. |
|
|
Smith Charts — As
a Test Designer, ICAP/4 Professional, or RF user, you can now plot
Smith Charts by reading in a Touchstone formatted file. Smith Charts
may be used to calculate impedance determination, impedance matching,
noise figure optimization, or stability. |
The forward converter Z-Transform model extends the approach
used by Yaakov [1] in which behavioral relationships were directly implemented.
Rydley [2] demonstrated the need for modeling transport delay. However,
in this case, a Z-Transform model implements this affect directly. Various
loss mechanisms can be added to the large signal model in order to predict
DC-DC converter efficiency using an operating point analysis. The following
model elements will be discussed and the behavioral equations developed.
Table 1.
Efficiency Measurement
|
Rload |
Pload |
From OP Analyses |
From TRAN Analyses |
.2 |
125 |
62.5% |
62.7% (Ron set to match) |
.7 |
37.5 |
81.8% |
81.8% |
1 |
25.0 |
84.6% |
84.4% |
2 |
12.5 |
88.1 |
87.8 |
4.7 |
5.3 |
89.5% |
89.5% |
10 |
2.5 |
88.6% |
88.8% |
20 |
1.25 |
85.6% |
87.2% |
Model
Elements
1. Sampling theory
2. CM/VM Control law
3. Fraction Order Hold
4. Z Transform Integrator
5. CCM/DCM mode switching
6. Input and Output current/voltage equations
7. Magnetizing inductance effects
8. Leakage inductance induced losses
9. Transistor switching loss
10. Switched capacitor loss model
11. Resistive Losses
12. Core Loss
Sampling
Theory
To achieve a large signal model, the cycle-by-cycle
switching difference equation must be derived. It is not suitable to
assume conservation of power through the converter at each instant because
the reactive elements store and release energy. Z-Transform theory predicts
the time domain performance at each sampling instant. For the class
of PWM circuits, it is convenient to set this instant at the beginning
of each pulse as shown in Figure 1a and 1b. This is not very satisfying
to the end user because it under-estimates the output current. To overcome
this we calculate an average output current based on our knowledge of
the sampled data state.
Z
Transform Integrator
The waveforms for DCM and CCM are shown in Figure 1a
and 1b, respectively. We begin each cycle with the previous sample current,
I1. The calculated result at the end of the cycle, Io, is fed into the
input of the transmission line. The sample delay, Z^-1, is e^sT in continuous
time and it can be simply modeled as a transmission line. The transmission
line works for both AC and Transient analysis, thereby preserving the
small signal characteristics needed for stability analysis and also
the large signal description needed for the transient case. Idischarge
will change with the inductor conduction mode. The representation in
Figure 1c is actually a bootstrap integrator. If the gain exceeds 1.0,
the output will grow to infinity. The gain was set to .995 to prevent
numerical round off from making this circuit unstable. Reducing Gmin
increases numerical precision. As Gmin is reduced, the integrator gain
can be made to approach 1.0; however, DC convergence will eventually
fail for sufficiently small Gmin.
Figure 1a. The
Z-Transform sampling instant is shown at the beginning of
each pulse of this class of PWM circuits. Figure 1a shows
the waveform for DCM, where lout = Ddis*IDischarge/2.
Figure 1b. The
Z-Transform sampling instant is shown at the beginning of
each pulse of this class of PWM circuits. Here the waveform
for CCM, where lout = Ddismax * (10 - Idischarge/2), is shown.
|
CCM/DCM
Mode Switching
The terms used in Figure 1a and 1b are defined as:
Don |
The fraction of the period the switch
is ON. |
DDIS |
The number of periods to fully discharge the inductor |
DDISMAX |
The time available for discharge in the current period.
If Ddismax > Ddis then the mode is CCM, else its DCM. 1-2
MegHz switching Voltage Regulator Module.
|
DDIS = 1 – Don
Icharge = (Vin-Vout)*Don/(L*F)
If (Mode == CCM)
Idischarge = (I1 + Icharge)* DDISMAX / DDIS
else
Idischarge = Vout * DDIS/(L*F) |
You’ll find these equations modified slightly in the
actual model to account for snubber delay and to prevent divide by zero
during numerical iterations. Figure 2 illustrates how that was done.
Notice that the ModeCCM switch has a small hysteresis in order to prevent
numerical oscillation. You can see that variables like Don are nodes
so that the node voltage, used in the equations, becomes V(Don).
Figure 1c.
A transmission line models Z^-1.
Figure 2. Calculating
average output current.
|
Output
Current/Voltage Equations
Figure 3 shows the method to calculate average output
current. These currents are sent to the output in 2 parts; charging
and discharging. The charging component (B24) is dependent on Don while
the discharging part (B4) is controlled by the mode. In both cases,
the average is taken so the output voltage and currents will be at their
average values rather than the sample data points. Diode, D2, prevents
negative output. It has a low emission coefficient, 3m, which has its
internal temperature set to 27 deg. C. independent of circuit temperature.
One of the general problems with this model is that it can find negative
solutions that have no interest. It is even possible for a positively
clamped state variable to have a negative result computed during the
matrix solution. These negative values can result in instability as
well as divide by zero overflows. That’s why you see many expressions
with positive clamped results (expression > 0 ? expression : 0).
Figure 3. Equations
used to model load current, series and parallel losses.
|
Figure 4. The input voltage
and current equations.
|
On the input side, illustrated by Figure 4, conduction
occurs only during Don. B1 computes the basic relationship. Notice the
magnetizing current, Imo, term. Its equations will be developed later
on. B23 and B18 introduce losses. B18 reduces the input voltage, simulating
resistive losses in the power switch and transformer. B23 accounts for
leakage inductor, core and switched capacitor losses. B18 and B23 are
used to model resistive and conductive losses using models that are
well behaved around zero. The following equation for B18 shows how the
parameter del prevents a singularity at zero current while squaring
the current in the denominator forces a the term to be positive, thereby
eliminating any chance for a divide by zero. Vin is referenced to simulator
ground, allowing the external inputs to be offset from ground.
V = I*R = I(VM2)*1/(1/R + del),
I^2R=P; R = P/I^2
V = I(VM2)*1/(I(VM2)^2/P + del) |
CM/VM
Control Law
The control signal, Vc, is converted from an external
ground on the input side using B10. The control signal is processed
through a fractional order hold, accounting for transport delay. This
FOH and the Z-Transform integrator realistically model the phase degradation
near one half the sampling (Nyquist) frequency.
The control equations are implemented as shown in Figure
5.
RB is the current sense resistor and the equations
assume that 1 volt will cause the maximum current flow. Setting the
parameters RB=0 and MC=1 will result in voltage mode control where the
duty cycle is set by v(Cntrl).
For CM controllers, the following
equations apply:
I1+Icharge+Imag = Vcontrol/Rb
Icharge = Don*(Vin-Vout)/(LF)
I1 + Don*Vin/(L*F) + Imag = Vcontrol/Rb
Don = (Vcontrol/Rb - I1-Imag)*(L*F)/Vin
Don = (Vcontrol - Rb(*I1-Imag))/(RB*Vin + MC)
To operate in VM set:
MC=1, RB=0
For CM MC = slope compensation as fraction of the slope |
Figure 5. Control law for
both CM and VM controllers.
|
What's
*exempt all about?
We have many requests from semiconductor vendors to
increase the size limit on our demo software so that they can show how
new ICs work. The IC models and circuits have become so complex that
to increase our complexity test so these circuit simulations would run
amounts to giving our software away. In order to accommodate the showcasing
of these products, Intusoft invented a method to exempt certain portions
of a drawing from the IsSpice complexity test. This allows a simple
circuit configuration to simulate but prevents it from working when
the user adds parts and functionality needed to make a production circuit.
Here’s how it works. Right after the .subckt line in
the Spice netlist, this line is inserted in the netlist:
*exempt arg1 arg2 arg3
This causes the demo version of IsSpice to scan the
remainder of the subcircuit netlist. If the subcircuit netlist passes
the scanning test, then the subcircuit is not processed by the complexity
testing procedure. The scanning method prevents topology changes in
the subcircuit from being accepted by the simulator. Each subcircuit
is scanned according to the rules presented in its *exempt line so that
*exempt lines can’t be applied to other subcircuits.
We do this to showcase our software so that the combination
of software and models is synergistic. The whole is more valuable than
the sum of its parts, or so we hope. You will see this feature in SpiceNet
drawings from ON Semiconductor and TI in the next few months.
Intusoft recently did a Bode plot for an internal loop
model using our new Z-Transform model for the switched inductor.
The Z-Transform model shows gain greater than 1 at
harmonics of the switching frequency. This behavior suggests instability
that would fold back below the power supply Nyquist frequency.
Note: The frequencies above the Nyquist frequency are
generally taken to mean that these frequencies fold back beginning at
½ the switching frequency, extending into negative frequencies
until reversing. This behavior is well known and used in heterodyne
detection systems.
Figure 6. Bode Plot with
Loop Cut at Internal Duty Ratio.
|
There is no evidence of instability in the outer voltage
control loop, as shown in Figure 7.
The switching models transient output voltage shows
ripple primarily at the switching frequency; however, there is a small
chaotic noise component. This is shown in Figure 8.
Figure 7. Bode Plot with
Loop Cut at Voltage Control.
Figure 8. Switch model's
transient output voltage shows ripple at the switching frequency.
|
Questions
for Future Study
1. Is the unexpected transient ripple really there
or is this just numerical noise?
2. What is to be made of the AC response above the switching frequency?
Test Designer
Key to Structural Integrity
|
Test Designer is software which designs test strategies
to assure structural integrity. There have long been arguments over
structural versus functional testing; however, structural test MUST
be performed to assure reliability. Here’s why:
Did
You Know? You've designed redundant or fault tolerant systems
in order to ride through a single point failure. But over time it is
increasingly likely that a failure has occurred and that if undetected
a catastrophic event is likely. The availability of redundant systems
is fundamental to avoiding catastrophe. Either built in self test, BIST
or auxiliary testing is needed to perform structural testing that assures
proper component operation. This kind of testing effectively resets
the accumulated failure rate to zero for each part that passes the test.
After a successful test; then when you begin a mission, you have full
advantage of the redundant systems.
How
is it Done? It is necessary to analyze the results of each
likely component failure. Then a test strategy is necessary to detect
each fault. Work done in the USA for commercial and military aircraft
systems has produced methods for sequencing these tests to produce a
decision strategy that isolates faults. Putting this together with simulation
technology, it is possible to complete this apparently monumental task.
Say
that Again Undetected part failures make a catastrophic event
more and more likely the longer a system operates without fault diagnosis.
Fault diagnostic test design is possible today using simulation technology
and established test design techniques. That combination has been put
together with Intusoft’s Test Designer. You can even model electro-mechanical
components and perform system level diagnostic testing.
Show
me an Example Let’s try an automotive steering system. In
the near future electrical systems will replace the mechanical-hydraulic
method currently in use. The steering wheel may become something like
a force feedback joystick. It sends a voltage to a controller and road
feedback is transmitted back to the driver. The controller requires
input position, output position and a control law. If analog, many resistors
and capacitors are used to obtain a control transfer function. If digital,
properly functioning A/D and D/A converters and other IC chips must
be working. Both require a switch mode power controller to drive an
electric motor that steers the car. Using affordable parts, a failure
rate of 1 to 100 failures per 100 million hours per part can be expected.
That produces about 10 failures per million hours for a system. While
that sounds pretty good; if you take 10 million cars for twenty years,
you get 17.5 million failures! Even lower production rate systems like
commercial and military aircraft are exposed to incredible risk levels.
Fault tolerance is no longer an option — it is a requirement. Now here’s
the rub; if you want a mission success of 99.99% (that’s not unusual
for taking a trip to the store and returning alive!), then you can’t
go for long without testing the redundant systems in order to purge
broken parts. It’s no longer just an oil change every 3000 miles!
Did
you know that even today, your low-tech automobile has a dashboard
light that checks your redundant breaking hydraulic fluid level?
|
New Spice Transistor
Models
|
DarNPN(8) |
2N6283, 2N6284, 2N6301, 2N6350, 2N6351, 2N6353, 2N6384, 2N6385
|
DarPNP(3) |
2N6649, 2N6650, 2N6650 |
NPN(27) |
2N2060, 2N2060L, 2N2219AL, 2N2222AU, 2N2432, 2N2919, 2N2919L,
2N2920L, 2N3019S, 2N3418, 2N3418S, 2N3419, 2N3419S, 2N3420, 2N3420S,
2N3421, 2N3421S, 2N3501L, 2N3735, 2N3735L, 2N4150, 2N4150S, 2N5237S,
2N5238S, 2N5662, 2N5663, 3N3439 |
PNP(24) |
2N2605, 2N3251A, 2N3468, 2N3468L, 2N3634L, 2N3635, 2N3635L,
2N3636L, 2N3637, 2N3637L, 2N3763, 2N3763L, 2N3764, 2N3765, 2N3810,
2N3810L, 2N3811, 2N3811L, 2N3867, 2N3867S, 2N3868, 2N3868S, 2N5416,
2N5416S |
PowNPN(36) |
2N2151, 2N2814, 2N2880, 2N3442, 2N3716, 2N3739, 2N3749, 2N3996,
2N3997, 2N5004, 2N5038, 2N5039, 2N5157, 2N5237, 2N5238, 2N5250,
2N5251, 2N5302, 2N5660, 2N5664, 2N5665, 2N5666, 2N5667, 2N5672,
2N5685, 2N5686, 2N6033, 2N6249, 2N6250, 2N6251, 2N6308, 2N6338,
2N6546, 2N6547, 2N6676, 2N6678 |
PowPNP(11) |
2N3740, 2N3741, 2N4399, 2N5005, 2N5683, 2N5684, 2N5745, 2N6211,
2N6212, 2N6213, 2N6438 |
[1] Ben-Yaakov, S. “Average Simulation of PWM Converters
by Direct Implementation of Behavioral Relationships,” IEEE Applied
Power Electronics Conference (APEC, 1993), pp. 510-516
[2] Ridley, R. B. “A New Small-Signal Model for Current-Mode
Control,” PhD. dissertation, Virginia Polytechnic Institute and State
University, 1990.
|