PSIM User Manual
PSIM User Manual
PSIM User Manual
Powersim Inc.
Chapter : -7
PSIM Users Guide
Version 10.0
Release 3
April 2015
Disclaimer
Powersim Inc. (Powersim) makes no representation or warranty with respect to the adequacy or accuracy of this
documentation or the software which it describes. In no event will Powersim or its direct or indirect suppliers be liable for any
damages whatsoever including, but not limited to, direct, indirect, incidental, or consequential damages of any character
including, without limitation, loss of business profits, data, business information, or any and all other commercial damages or
losses, or for any damages in excess of the list price for the licence to the software and documentation.
Powersim Inc.
email: info@powersimtech.com
http://www.powersimtech.com
-6 Chapter :
Contents
1 General Information
1.1 Introduction 1
1.2 Circuit Structure 3
1.3 Software/Hardware Requirement 3
1.4 Installing the Program 3
1.5 Simulating a Circuit 4
1.6 Simulation Control 4
1.7 Component Parameter Specification and Format 5
Chapter : -5i
3.7 View Menu 37
3.8 Option Menu 38
3.9 Label Menu 38
3.10 Exporting Data 39
ii
-4 Chapter :
4.6.1 Reference Direction of Mechanical Systems 77
4.6.2 DC Machine 84
4.6.3 Induction Machines 80
4.6.4 Induction Machine with Saturation 82
4.6.5 DC Machine 84
4.6.6 Synchronous Machine with External Excitation 89
4.6.7 Permanent Magnet Synchronous Machine 92
4.6.8 Permanent Magnet Synchronous Machine with Saturation 94
4.6.9 Switched Reluctance Machine 97
4.6.10 Nonlinear Switched Reluctance Machine 99
4.6.11 Motor Control Blocks 100
4.6.11.1 Maximum-Torque-Per-Ampere Control 100
4.6.11.2 Field Weakening Control 101
4.7 MagCoupler Module 104
4.7.1 MagCoupler-DL Block 104
4.7.2 MagCoupler Block 105
4.8 MagCoupler-RT Module 108
4.9 Mechanical Elements and Sensors 112
4.9.1 Mechanical Elements and Sensors 112
4.9.1.1 Constant-Torque Load 112
4.9.1.2 Constant-Power Load 112
4.9.1.3 Constant-Speed Load 113
4.9.1.4 General-Type Load 113
4.9.1.5 Externally-Controlled Load 114
4.9.2 Gear Box 114
4.9.3 Mechanical Coupling Block 115
4.9.4 Mechanical-Electrical Interface Block 115
4.9.5 Speed/Torque Sensors 116
4.9.6 Position Sensors 118
4.9.6.1 Absolute Encoder 118
4.9.6.2 Incremental Encoder 119
4.9.6.3 Resolver 119
4.9.6.4 Hall Effect Sensor 120
4.10 Thermal Module 120
4.10.1 Diode Thermal Model 121
4.10.1.1 Diode Device in Database 121
4.10.1.2 Diode Loss Calculation 122
4.10.2 IGBT Thermal Model 123
4.10.2.1 IGBT Device in Database 123
4.10.2.2 IGBT Loss Calculation 125
4.10.3 MOSFET Thermal Model 125
4.10.3.1 MOSFET Device in Database 127
4.10.3.2 MOSFET Loss Calculation 129
4.10.4 Inductor Thermal Model 131
4.10.4.1 Inductor in Database 131
4.10.4.2 Inductor Loss Calculation 131
4.10.5 Device Database Editor 132
4.10.6 Adding a Switching Device to Database 134
4.11 Renewable Energy Module 138
4.11.1 Solar Modules 138
4.11.1.1 Solar Module - Functional Model 138
4.11.1.2 Solar Module - Physical Model 139
4.11.1.3 Solar Module - cSi and Thin-File Models 141
4.11.2 Wind Turbine 143
4.11.3 Lithium-Ion Battery Model 145
Chapter : iii
-3
4.11.4 Ultracapacitor Model 146
iv
-2 Chapter :
5.5.6 Convolution Block 180
5.5.7 Memory Read Block 180
5.5.8 Data Array 181
5.5.9 Stack 181
5.5.10 Multi-Rate Sampling System 182
5.6 SimCoupler Module 182
5.6.1 Set-up in PSIM and Simulink 182
5.6.2 Solver Type and Time Step Selection in Simulink 184
5.7 CosiMate Links 186
5.8 Design Suite Blocks 186
5.8.1 Torque Control (PMSM) 186
5.8.2 Dynamic Torque Limit Control (PMSM) 187
5.8.3 Dynamic Torque Limit Control (Nonlinear PMSM) 188
5.8.4 Voltage Control (PMSM) 189
5.8.5 DC-DC Charging Control 190
5.8.6 DC-DC Discharging Control 191
5.8.7 DC-DC Regeneration Control 192
6 Other Components
6.1 Switch Controllers 193
6.1.1 On-Off Switch Controller 193
6.1.2 Alpha Controller 193
6.1.3 PWM Lookup Table Controller 194
6.2 Sensors 196
6.3 Probes, Meters, and Scopes
6.3.1 Probes and Meters 196
6.3.2 Voltage/Current Scopes 198
6.4 Function Blocks 200
6.4.1 Control-Power Interface Block 200
6.4.2 Initial Values 213
6.4.2.1 abc-dqo Transformation 201
6.4.2.2 abc-/ Transformation 202
6.4.2.3 /-to-dq Transformation 203
6.4.2.4 Cartesian-Polar Transformation 204
6.4.3 Math Function Blocks 204
6.4.4 Lookup Tables 205
6.4.5 C Block 207
6.4.6 Simplified C Block 208
6.4.7 External DLL Blocks 209
6.4.8 Embedded Software Block 210
6.5 IC Models 211
6.5.1 PWM IC 211
6.5.2 Driver IC 212
6.5.3 555 Timer 213
6.6 Initial Values 213
6.7 AC Analysis 215
6.7.1 AC Sweep 215
6.7.2 AC Sweep Multi-Sine 218
6.8 Parameter Sweep 219
7 Sources
Chapter : -1v
7.1 Constant 221
7.2 Time 221
7.3 Ground 221
7.4 Voltage and Current Sources 222
7.4.1 DC Source 222
7.4.2 Sinusoidal Source 222
7.4.3 Square-Wave Source 223
7.4.4 Triangular/Sawtooth Sources 224
7.4.5 Step Sources 225
7.4.6 Piecewise Linear Source 225
7.4.7 Random Source 226
7.4.8 Math Function Source 227
7.4.9 Voltage/Current-Controlled Sources 227
7.4.10 Nonlinear Voltage-Controlled Sources 229
Index 235
vi
0 Chapter :
1
General Information
1.1 Introduction
PSIM1 is a simulation software specifically designed for power electronics, motor drives, and power conversion
systems. With fast simulation speed and friendly user interface, PSIM provides a powerful simulation
environment to meed your simulation and development needs.
PSIM includes the basic package, as well as the following add-on options:
Motor Drive Module: Built-in electric machine models and mechanical load models for motor drive
system studies.
Digital Control Module: Discrete library elements such as zero-order hold, z-domain transfer function
blocks, quantization blocks, digital filters, for digital control system analysis.
SimCoupler Module: Interface between PSIM and Matlab/Simulink2 for co-simulation.
Thermal Module: Library elements and functions calculate semiconductor device losses and
inductor losses.
Renewable Energy Module: Library elements, such as solar module, wind turbine, battery, and
ultracapacitor models, for renewable energy applications.
SimCoder3 Module: Function for automatic code generation capability.
F2833x Target: Library elements for automatic code generation for TI F2833x series DSP.
F2803x Target: Library elements for automatic code generation for TI F2803x series DSP.
MagCoupler Module: Interface between PSIM and the electromagnetic field analysis software
JMAG4 for co-simulation.
MagCoupler-RT Module: Link between PSIM and JMAG-RT4 data files.
ModCoupler5 Modules: Interface between PSIM and ModelSim6 for co-simulation. There are two
versions of the interface: ModCoupler-VHDL that supports VHDL code, and
ModCoupler-Verilog that supports Verilog code.
Motor Control Design Suite: Pre-built templates for induction motor and linear/nonlinear PMSM drives.
HEV Design Suite: Pre-built templates for hybrid electric vehicle (HEV) powertrain system
design.
In addition, a link is provided between PSIM and CosiMate7, a co-simulation framework that supports various
software such as Matlab/Simulink, ModelSim, Saber (from Synopsys), Easy5 and Adams (from MSC
Software), Inventor (from Autodesk), AMESim (from LMS), GT-Power (from Gamma Technologies), etc.
Through CosiMate, PSIM can perform co-simulation with these software. For more information on CosiMate,
please visit www.chiastek.com.
Also, PSIM links with the software SmartCtrl5 for control loop design. SmartCtrl is designed specifically for
power converter applications. For more information on SmartCtrl, please refer to SmartCtrl Users Guide.
1. PSIM and SIMVIEW are registered trademarks of, and are copyright by, Powersim Inc., 2001-2015
2. Matlab and Simulink are registered trademarks of the MathWorks, Inc.
3. SimCoder is a trademark of Powersim Inc., and is copyright by Powersim Inc., 2008-2015
4. JMAG and JMAG-RT are copyright by JSOL Corporation, 1997-2015
5. ModCoupler and SmartCtrl are copyright by Universidad Carlos III de Madrid, 2011-2015
6. ModelSim is a registered trademark of Mentor Graphics Co.
7. CosiMate is copyright by ChiasTek Inc., 2010-2015
Hardware
Design Simulation Implementation
F2833x &
F2803x
PSIM Targets
SmartCtrl Thermal Auto Code DSP
Analysis Power Generation Hardware
Thermal Electronics SimCoder
Motor Control
Design Suite
Solar/Wind ModCoupler
Power ModelSim
HEV Design Suite Control Renewable Energy Motor Drives
Digital Control Motor Drive
FPGA
SimCoupler MagCoupler/ DLL
MagCoupler-RT
Matlab/ 3rd-party
Simulink JMAG CosiMate
Software
Link to Saber, AMESim, etc.
The PSIM simulation environment consists of the circuit schematic program PSIM, the simulator engine, and
the waveform processing program SIMVIEW1. The simulation process is illustrated as follows.
This manual covers both PSIM and all add-on Modules except SimCoder Module, ModCoupler Modules, HEV
Design Suite, and Motor Control Design Suite. The document SimCoder User Manual describes the use of the
SimCoder Module and associated hardware targets. The document ModCoupler User Manual describes the use
of the ModCoupler Module. The document "Tutorial - HEV Design Suite.pdf" describes the use of the HEV
Design Suite. The document "Tutorial - Motor Control Design Suite.pdf" describes the use of the Motor Control
Design Suite.
The organization of this manual is as follows:
Power Circuit
Switch
Sensors
Controllers
Control Circuit
The power circuit consists of switching devices, RLC branches, transformers, and coupled inductors. The
control circuit is represented in block diagram. Components in s-domain and z-domain, logic components (such
as logic gates and flip flops), and nonlinear components (such as multipliers and dividers) are used in the
control circuit. Sensors are used to measure power circuit quantities and pass them to the control circuit. Gating
signals are then generated from the control circuit and sent back to the power circuit through switch controllers
to control switches.
The parameters in the Parameters tab are used in the simulation. The information in the Other Info tab, on the
other hand, is not used in the simulation. It is for reporting purposes only and will appear in the parts list in
View >> Element List in PSIM. Information such as device rating, manufacturer, and part number can be
stored under the Other Info tab.
PSIMs schematic program provides interactive and user-friendly interface for circuit schematic entry and
editing. The PSIM user interface consists of an integrated set of windows, tools, menus, toolbars, and other
elements that allow you to create, simulate, and refine your circuits in one places.
By default, the menu bar and the standard toolbar appear on top of the window, while the frequently used
element bar appears at the bottom, and the Project View is on the left hand side.
The Project View window provides an organized tree view of users projects and their related files, as well as
the simulation result graphs. Each PSIM circuit is treated as a project. The following content are displayed in
stacking tiers in the project view:
Project Name: Usually this is the same as the top level PSIM circuit file name.
Documents: Any files related to the PSIM circuit: descriptive, parameter, lookup table, etc.
Schematic: The PSIM schematic files, top level and subcircuits.
Graphs: All probes in the PSIM schematic are included in the graph list.
On the right hand side is the Design window. This is a graphic editor where users can build and editor their
simulation circuit schematics. User may arrange schematics in the Design window in tiles (as shown above), or
in tabs (as shown below).
Place Once an element is selected from the menu, the image of the element will appear on the
screen and move with the mouse. Click the left button of the mouse to place the element at
desired location on schematic.
Select Element(s) To select an existing element on a schematic, click on the element. A rectangle will appear
around the element. To select a section of a circuit, keep the left button of a mouse
pressed and drag the mouse until the rectangle covers the selected area.
Rotate Before the element is placed, right click to rotate the element. After an element is
selected, select Edit >> Rotate to rotate the element.
Wire To connect a wire between two nodes, select Edit >> Wire. The image of a pen will
appear on the screen. To draw a wire, keep the left button of the mouse pressed and drag
the mouse. A wire always starts from and end at a grid intersection.
For easy inspection, a floating node is displayed as a circle, and a junction node is
displayed as a solid dot.
Label If two or more nodes are connected to the same label, they are connected. It is equivalent
as though they were connected by wire. Using labels will reduce the cross-wiring and
improve the schematic layout.
The text of a label can be moved. To select the text, left click on the label, then press the
Tab key.
Assign To assign the parameters of an element, double click on the element. A dialog box will
appear. Specify the values and hit the <Return> key or click on OK.
Move To move an element or a circuit block, select the element/circuit block and drag the mouse
while keeping the left button pressed.
Pan Schematic To scroll schematic, right click and drag the mouse.
The diamonds on the four sides represent the connection nodes and the positions of the subcircuit. They
correspond to the connection nodes of the subcircuit block on the right. There are no diamonds at
the four corners since connections to the corners are not permitted.
When a diamond is selected, it is colored red. By default, the left diamond at the top is selected and
marked with red color. Click on the desired diamond to select and to specify the port name.
In this example, in the main circuit chop.sch, there are four linking nodes, two on the left side and
The creation of the subcircuit is now complete. Save the subcircuit, and go back to the main circuit.
File: sub.sch
File: main.sch
In the window, the diamonds marked red are the connection nodes of the subcircuit block, in exactly
the same positions as appearing in the main circuit.
- Use the drawing tool to create/edit the image for the subcircuit block. If the drawing tool is not
already displayed, go to the View menu and check Drawing Tools. Click on Zoom In and Zoom
Out icons on the toolbar to adjust the size of the image working area.
After the image is created, the pop-out window will appear as follows.
Gain of the
PI controller
- Place the cursor inside the input field of the dialog window for the gain, and change the gain either
by pressing on the upper/down arrow keys on the keyboard, or by entering a new value and then
clicking on Apply. Watch how the waveforms change as the gain is changed.
Other parameters, such as current reference, dc input voltage, inductance, capacitance, and load
resistance, can be changed in the similar way.
- Branch currents can also be displayed in the free-run mode. To display the inductor current, for
example, right click on top of the inductor, and a menu will appear. Choose Current Scopes and
the branch current name.
- An image of the current scope (similar to the voltage scope image, but without connection
terminals) will appear. Double click on the scope to expand and view the inductor waveform.
Below is how the window would look like with both the voltage scope and the current scope.
Other branch currents, such as capacitor current, load current, diode current, or MOSFET switch
current, can be displayed in the similar way.
Voltage scope
Toolbar
icon area
Predefined
icon images
Icon
editing area
Choose View >> Custom Toolbars. The Custom Toolbars dialog window will appear. Choose New Tool-
bar, and the following window will appear.
Specify the Toolbar Name as "new".
Draw the AND gate icon in the icon editing area. Or if the icon is already available in the predefined icon
images, select the icon and copy it to the icon editing area.
Choose View >> Custom Keyboard. The Custom Toolbars dialog window will appear. Choose New Tool-
bar, and the following window will appear.
In the section Add Shortcut Key, select the option Elements. Then find and highlight the element "Resis-
tor".
Move the cursor into the input field of Press new shortcut key, and press the key "r" on the keyboard. Then
click on Assign.
The key "r" will be assigned to the resistor, and the definition will appear in the Current Shortcut Key list.
Netlist name
Standard image
Secondary image
In this case, the inductance and capacitance values will be defined through the interface, and need to appear in
the property window of the new LC filter element. Therefore, the parameter value for the inductance needs to be
defined as a variable, in this case, L, and the value for the capacitance as C.
Then from Subcircuit >> Edit Default Variable List, add the variables L and C as the default variables. This
step is necessary as the new element obtains the parameter information from the default variable list. The
default variable list window should appears as follows.
Name: Name of the new inductor element as it appears in the PSIM library
Description: Description of the new inductor element
File Path: The location of the DLL file "inductor_model.dll" that models the new inductor element.
The DLL file must be placed in the "lib" sub-folder in the PSIM directory.
Input Nodes: Number of input terminals of the new element.
Output Nodes: Number of output terminals of the new element.
Hide (menu): Leave this box unchecked. If this box is checked, this element will not appear in the
library.
Help File: On-line help file associated with this element. This file must be placed in the "help" sub-
folder in the PSIM directory. When the Help button is clicked in the property dialog window, this
file will be displayed. This file can be a text file (which can be opened by a text editor such as
NotePad) or a HTML file.
In the next dialog window, set the new element size as: Width = 5, and Height = 2. Then create an image for
this element, or accept the default image.
Click on the buttons Save Image Library and Update Menu. The new element will appear in the library
and will be ready to use.
The information regarding the number of parameters and the parameter description for the new inductor
element is obtained from the DLL file automatically. In this case, the new element will have one parameter as
"Inductance".
Click on the Pins tab. In this tab, the connection pins will be defined. Click on the Add button, and define
Location as "1 - Left", Name as "A", and Number as "1". Leave the boxes for Dot, clk, and Overline
unchecked.
If the boxes are checked, for Dot, there will be a circle between the connection pin and the block; for
clk, a symbol > will be added inside the block to indicate that it is a clock; and for Overline, an
overline will be added to the name.
Click on the Add button, and define Location as "2 - Left", Name as "B", and Number as "2". Repeat the
same process for the rest of the pins. The dialog windows below show the definitions of all the connection
pins.
SIMVIEW is PSIMs waveform display and post-processing program. The following shows simulation
waveforms in the SIMVIEW environment.
SIMVIEW reads data in either ASCII text format or SIMVIEW binary format. The following shows a sample
text data file:
Time Isa Isc Isb Tem_IM
5.000000000E-006 0.000000000E+000 0.000000000E+000 0.000000000E+000 7.145888260E-048
1.000000000E-005 0.000000000E+000 0.000000000E+000 0.000000000E+000 1.082981714E-046
1.500000000E-005 0.000000000E+000 0.000000000E+000 0.000000000E+000 5.408644357E-046
2.000000000E-005 1.139566166E-001 -2.279132474E-001 1.139566166E-001 1.613605209E-017
2.500000000E-005 5.072914178E-001 -1.014582858E+000 5.072914178E-001 3.598226665E-015
... ... ... ...
If the Auto-Grid box is checked, the number of axis divisions will be automatically determined. The data range
and grid division, however, can be manually set.
The following figure shows a sine waveform is chosen as the X-axis versus a cosine waveform in the Y axis.
Edit Box
All the data variables available for display are in the Variables Available box, and the variables currently being
displayed are in the Variables for Display box. After a variable is highlighted in the Variables Available box, it
can be added to the Variables for Display box by clicking on Add ->. Similarly, a variable can be removed
from display by highlighting the variable and clicking on <- Remove.
In the Edit Box, mathematical expressions can be specified. The mathematical expressions may contain
brackets and is not case sensitive. The following math functions are allowed:
+ addition
- subtraction
* multiplication
/ division
^ to the power of [Example: 2^3 = 2*2*2]
SQRT square-root function
SIN sine function
COS cosine function
TAN tangent function
ATAN inverse tangent function
EXP exponential (base e) [Example: EXP(x) = ex]
LOG logarithmic function (base e) [Example: LOG(x) = ln (x)]
LOG10 logarithmic function (base 10)
ABS absolute function
SIGN sign function [Example: SIGN(1.2) = 1; SIGN(-1.2) = -1]
AVG moving average function, calculates the average of the curve y up to the measured
point.
AVGX periodic average function AVGX(y,Tp) where y is the curve name and Tp is the time
interval where the average is calculated. Calculates the average of the curve y in each
time segment with interval of the time segment as Tp. For example, AVGX(V1,
0.016667) will calculate the average of the curve V1 at the 60-Hz interval.
INT integration function
Type an expression in the Edit Box and click the "Add ->" button to add the calculated curve to the screen. .
Highlight the expression on the right and click the "<- Remove" button, the expression will be moved into the
Edit Box for further editing.
Once Measure is selected, an individual curve can be selected by clicking on the pull-down menu
on the Measure toolbar. The functions, Max, Min, Next Max, Next Min, Avg, and rms, can
be used to evaluate the curve. Note that these functions are only enabled in the Measure mode.
Memory
Expression
Result
One key feature of the calculator is that it provides 9 memory spaces. By double clicking on a number in the
The following chapters provide the detailed information of the PSIM library elements. This chapter focuses on
the Power circuit components.
Parameters Description
Resistance Resistance, in Ohm
Voltage Rating Voltage rating of the component, in V, for Level 2 model only.
Current Flag Flag for branch current waveform display.
If this flag is set 1, the current flowing through this component will be saved to
display in SIMVIEW and Runtime Graphs.
The current is positive when it flows into the dotted terminal of the branch.
Attributes of Level-2 Model:
Parameters Description
Resistance Resistance, in Ohm
Voltage Rating Voltage rating of the component, in V.
Power Rating Power rating of the component, in W.
Inductance ESL Equivalent series inductance (ESL), in H.
Parallel Capacitance Equivalent parallel parasitic capacitance, in F.
Current Flag Flag for branch current waveform display.
If this flag is set 1, the current flowing through this component will be saved to
display in SIMVIEW and Runtime Graphs.
The current is positive when it flows into the dotted terminal of the branch.
Parameters Description
Inductance Inductance, in H.
Initial Current Initial inductor current, in A
Current Flag Flag for branch current waveform display
Attributes of Level-2 Model
Parameters Description
Inductance Inductance, in H.
Peak Current Rating Peak current rating of the inductor, in A.
Series Resistance Equivalent series resistance (ESR), in Ohm.
Parallel Resistance Leakage resistance in parallel with the inductor, in Ohm.
Parallel Capacitance Parallel parasitic capacitance, in F.
Initial Current Initial inductor current, in A
Current Flag Flag for branch current waveform display.
Equivalent Circuit:
Level 2 Model
Level 1 Model
Rp
Parallel
Resistance
Cp
Parallel
Capacitance
L ESR
Parameters Description
Capacitance Capacitance, in F.
Initial Voltage Initial capacitor voltage, in V
Current Flag Flag for branch current display
Attributes of Level-2 Model:
Parameters Description
Capacitance Capacitance, in F.
Voltage Rating Voltage rating of the component, in V.
RMS Current Rating RMS current rating, in A.
Resistance ESR Equivalent series resistance (ESR), in Ohm.
Inductance ESL Equivalent series inductance (ESL), in H.
Leakage Resistance Leakage resistance, in Ohm.
Initial Voltage Initial capacitor voltage, in V
Current Flag Flag for branch current waveform display
Equivalent Circuit:
Level 1 Model Level 2 Model
Parallel
Rp
Resistance
Capacitor
ESR C ESL
Parallel
Rp Resistance
Capacitor (Electrolytic)
ESR C ESL
Parameters Description
Resistance Resistance, in Ohm
Inductance Inductance, in H
Capacitance Capacitance, in F
Initial Current Initial inductor current, in A
Initial Cap. Voltage Initial capacitor voltage, in V
Current Flag Flag for branch current waveform display.
Note: For combined R-L-C brand, at least one of the parameters R, L, or C has to be a non-zero value.
Parameters Description
Resistance Resistance, in Ohm
Inductance Inductance, in H
Capacitance Capacitance, in F
Current Flag_A; Current Current flags for Phase A, B, and C of three-phase branches, respectively.
Flag_B; Current Flag_C If this flag is set 1, the current flowing through this branch will be saved to the file
for display in SIMVIEW. It will also be available for display in the runtime graphs
(under Simulate >> Runtime Graphs).
The dot in the image indicates phase A. The current is positive when it flows into
the dotted terminal of the branch.
A a
B b
C c
N n
Attributes:
Parameters Description
Cable Length Length of the cable, in m
Operating Frequency Operating frequency of the cable, in Hz
Resistance Rd Positive sequence resistance Rd of the cable, in Ohm
Reactance Xd Positive sequence reactance Xd of the cable, in Ohm
Resistance Ro Zero sequence resistance Ro of the cable, in Ohm
Reactance Xo Zero sequence reactance Xo of the cable, in Ohm
Capacitance Co Zero sequence capacitance Co of the cable, in F
In the image, the bottom node N or n on each side of the cable is the cable screen. It is generally connected to
the ground on the side of the cable.
The cable parameters should be available from manufacturers datasheet. When they are not available,
assuming that each phase of the cable has a resistance of R, self inductance of L, and mutual inductance of M,
the positive sequence and zero sequence parameters can be calculated as follows:
4.1.7 Rheostat
A rheostat is a resistor with a tap.
Image:
t
k m
Attributes:
Parameters Description
Total Resistance Total resistance of the rheostat R (between Node k and m), in Ohm
Tap Position (0 to 1) The tap position Tap. The resistance between Node k and t is: R*Tap.
Current Flag Flag for the current that flows into Node k.
Attributes:
Parameters Description
Current vs. Inductance Characteristics of the current versus the inductance (i1, L1), (i2, L2), etc.
Current Flag Flag for the current display
The nonlinear B-H curve is represented by piecewise linear approximation. Since the flux density B is
proportional to the flux linkage and the magnetizing force H is proportional to the current i, the B-H curve can
be represented by the -i curve instead, as shown below.
(B)
3
2 Inductance L = / i
1
i1 i2 i3 i (H)
The inductance is defined as: L = / i, the ratio of vs. i at each point. The saturation characteristics are defined
by a series of data points as: (i1, L1), (i2, L2), (i3, L3), etc.
In the images, the small circle indicates the input terminal of inductor 1 while the small square, triangle, + sign,
x sign, and * sign indicate the input terminal of inductor 2, 3, 4, 5, and 6.
Attributes:
Parameters Description
Lii (self) Self inductance of the inductor i, in H
Lij (mutual) Mutual inductance between Inductor i and j, in H
Initial Current ii Initial current in Inductor i
Current Flag_i Flag for the current printout in Inductor i
The following shows a coupled inductor with two branches.
+ v1 -
i1
i2
+ v2 -
Let L11 and L22 be the self-inductances of Branch 1 and 2, and L12 and L21 the mutual inductances, the branch
voltages and currents have the following relationship:
v1 L 11 L 12 d i1
= ------
v2 L 21 L 22 dt i
2
The mutual inductances between two windings are assumed to be always equal, i.e., L12 = L21.
In the SPICE software, coupled inductors are defined by the self inductances and the coupling factors. For the
coupled inductor with two branches, the coupling factor K is defined as:
L 12 = K L 11 L 22
The coupling factor has a value between 0 and 1. When the coupling factor is 1, it is perfect coupling. Note that
in PSIM, perfect coupling is not allowed. That is, the coupling factor cannot be equal to 1.
Also, when the leakage inductances are ignored (i.e. the coupling factor is close to 1), the self inductances are
proportional to the turns ratio squared. That is, if Branch 1 has N1 turns, and Branch 2 has N2 turns:
2
L1 N
----- = ----12-
L2 N2
Example:
Two mutually coupled inductors have the self inductances and mutual inductance as: L11 = 1 mH, L22 = 1.1 mH,
and L12 = L21 = 0.9 mH. The coupling factor is: K = 0.86. The specification of this element will be:
L11 (self) 1m
L12 (mutual) 0.9m
L22 (self) 1.1m
Attributes:
For resistance-type elements:
Parameters Description
Expression f(i) or f(i,x) Expression of v in terms of i and x [v = f(i) or v = f(i,x)]
Expression df/di The derivative of the voltage v versus current i, i.e. df(i)/di
Initial Value io The initial value of the current i
Lower Limit of i The lower limit of the current i
Upper Limit of i The upper limit of the current i
Parameters Description
Expression f(v) or f(v,x) Expression of i in terms of v and x [i = f(v) or i = f(v,x)]
Expression df/dv Derivative of the current i versus voltage v, i.e. df(v)/dv
Initial Value vo The initial value of the voltage v
Lower Limit of v The lower limit of the voltage v
Upper Limit of v The upper limit of the voltage v
A good initial value and lower/upper limits will help the convergence of the solution.
Example: Nonlinear Diode
The nonlinear element (NONI) in the circuit above models a nonlinear diode. The diode current is expressed as
a function of the voltage as: i = 10-14 * (e 40*v-1). In PSIM, the specifications of the nonlinear element will be:
4.2.1 Diode
The conduction of a diode or LED is determined by circuit operating conditions. PSIMs diode has 2 model
levels:
Level -1 Model: A diode is turned on when the positive bias voltage is greater than the threshold voltage,
and is turned off when the current drops to zero.
Level-2 Model: Includes leads parasitic inductance and reverse recovery parameters.
Image:
+ v -
Parameters Description
Forward Voltage Forward voltage Vd_th, in V. A diode starts to conduct when the positive bias
voltage is greater than Vd_th.
Resistance On-resistance Rd, in Ohm, after it starts to conduct.
Initial Position Flag for the initial diode position. If it is 0, the diode is off. If it is 1, the diode is on.
Current Flag Current flag for diode current display.
Attributes of the Level-2 Model:
Parameters Description
Forward Voltage Forward voltage Vd_th, in V. A diode starts to conduct when the positive bias
voltage is greater than Vd_th.
Resistance On-resistance Rd, in Ohm, after it starts to conduct.
Parasitic Inductance Series parasitic inductance, in H
Parallel Capacitance Parallel capacitance between the diode anode and cathode, in F. If the capacitance
value is 0, the capacitor is ignored and is removed from simulation.
Forward Current Forward current Ifwd under test conditions of diode reverse recovery characteristic,
in A.
Peak Reverse Current Peak reverse current Irm under test conditions, in A
Current Slope Current slope dI/dt under test conditions, in A/sec.
Reverse Recovery Time Reverse recovery time trr under test conditions, in sec.
Initial Position Flag for the initial diode position. If it is 0, the diode is off. If it is 1, the diode is on.
Current Flag Current flag for diode current waveform display.
Slope = 1/Rd
0 v
Vd_th
The reverse recovery time trr consists of two parts: ta and tb. During the period ta, the current increases from 0
to -Irm. Based on JEDEC (Joint Electron Device Engineering Council), tb is defined as the interval from the
time of -Irm to the time when the straight line from -Irm through -0.25*Irm (in red in the diagram) intersects
with the time axis.
For further information on how to use the Level-2 model, please refer to the tutorial "Tutorial - Diode model
with reverse recovery.pdf".
4.2.2 LED
A light-emitting diode (LED) emits light when it conducts. Its v-i characteristic is the same model as diodes
Level-1 model.
Image:
+ v -
Attributes:
Parameters Description
Forward Voltage LED forward threshold voltage Vd_th, in V. The LED starts to conduct and emit
light when the positive bias voltage is greater than Vd_th.
Resistance LED on-resistance Rd, in Ohm, after it starts to conduct.
Initial Position Flag for the initial LED position. If the flag is 0, the LED is off. If it is 1, the LED is
on.
Current Flag Current flag for LED current display.
VB
Attributes:
Parameters Description
Breakdown Voltage Breakdown voltage VB of the zener diode, in V
Forward Threshold Threshold voltage of the forward conduction (from anode to cathode), in V
Voltage
Forward Resistance On resistance of the forward conduction, in Ohm
Current Flag Flag for zener current output (from anode to cathode)
When the zener diode is positively biased, it behaviors as a regular diode. When it is reverse biased, it will block
the conduction as long as the cathode-anode voltage VKA is less than the breakdown voltage VB. When VKA
exceeds VB, the voltage VKA will be clamped to VB.
Note that when the zener is clamped, since the diode is modelled with an on-resistance of 10, the cathode-
anode voltage will in fact be equal to: VKA = VB + 10 * IKA. Therefore, depending on the value of IKA, VKA
will be slightly higher than VB. If IKA is very large, VKA can be substantially higher than VB.
DIAC:
A DIAC is a bi-directional diode. A DIAC does not conduct until the breakover voltage is reached. After that,
the DIAC goes into avalanche conduction, and the conduction voltage drop is the breakback voltage.
Image:
Attributes:
Parameters Description
Breakover Voltage Voltage at which breakover occurs and the DIAC begins to conduct, in V
Breakback Voltage Conduction voltage drop, in V
Current Flag Current flag
A K
Gate
Gate
Attributes:
Parameters Description
Voltage Drop Thyristor conduction voltage drop, in V
Holding Current Minimum conduction current below which the device stops conducting and returns to
the OFF state (for thyristor only)
Latching Current Minimum ON state current required to keep the device in the ON state after the
triggering pulse is removed (for thyristor only)
Initial Position Flag for the initial switch position (for thyristor only)
Current Flag Flag for switch current output
Note that for the TRIAC device, the holding current and latching current are set to zero.
There are two ways to control a thyristor or TRIAC. One is to use a gating block, and the other is to use a switch
controller. The gate node of a thyristor or TRIAC must be connected to either a gating block or a switch
controller.
The following examples illustrate the control of a thyristor switch.
Examples: Control of a Thyristor Switch
Gating Block
Alpha
Controller
This circuit on the left uses a switching gating block. The switching gating pattern and the frequency are pre-
defined, and remain unchanged throughout the simulation. The circuit on the right uses an alpha switch
controller. The delay angle alpha, in deg., is specified through the dc source in the circuit.
Parameters Description
Saturation Voltage The saturation voltage Vce_sat for npn, or Vec_sat for pnp, in V
Initial Position Initial switch position flag (0: off; 1: on)
Current Flag Switch current flag (0: no display; 1: display)
Example: Control of a npn Bipolar Junction Transistor
The circuit on the left uses a gating block, and the one on the right uses an on-off switch controller.
If the switch is controlled by a control circuit signal, as shown in the circuit on the right side, an on-off switch
controller must be used at the gate node.
The following shows another example of controlling the BJT switch. The circuit on the left shows how a BJT
switch is controlled in the real life. In this case, the gating voltage VB is applied to the transistor base drive
circuit through a transformer, and the base current determines the conduction state of the transistor.
This circuit can be modelled and implemented in PSIM as shown on the right. A diode, Dbe, with a conduction
voltage drop of 0.7V, is used to model the pn junction between the base and the emitter. When the base current
exceeds 0 (or a certain threshold value, in which case the base current will be compared to a dc source), the
comparator output will be 1, applying the turn-on pulse to the transistor through the on-off switch controller.
R DS ( on ) = R DS ( on )_ b ( 1 + K T ( T j T j _ b ) )
A switch can be controlled by either a gating block or a switch controller. They must be connected to the gate
(base) node of the switch.
Images:
MOSFET p-MOSFET MOSFET p-MOSFET
(RDS(on)) (RDS(on))
Tj Tj
Parameters Description
On Resistance The on resistance Rds_on of the MOSFET, in Ohm
Diode Forward Anti-parallel diode forward threshold voltage, in V
Voltage
Diode Resistance Anti-parallel diode on resistance, in Ohm
Initial Position Initial switch position flag for the transistor only (0: off; 1: on)
Current Flag Switch current flag for the whole module (the transistor plus the diode) (0: no
display; 1: display)
Parameters Description
Tj at Test Junction temperature Tj_b at test, in oC
RDS(on) at Test On-resistance RDS(on)_b at test, in Ohm
Temperature Temperature coefficient KT of the on-resistance, in 1/oC
Coefficient
Diode Forward Anti-parallel diode forward threshold voltage, in V
Voltage
Diode Resistance Anti-parallel diode on resistance, in Ohm
Initial Position Initial switch position flag for the transistor only (0: off; 1: on)
Current Flag Switch current flag for the whole module (the transistor plus the diode)
The following examples illustrate the control of a MOSFET switch.
On-off Controller
The circuit on the left uses a gating block, and the one on the right uses an on-off switch controller. The gating
signal is determined by the comparator output.
4.2.7 IGBT
An IGBT switch consist of an active switch with an anti-parallel diode. It is turned on when the gating signal is
high (when a voltage of 1V or higher is applied to the gate node) and the switch is positively biased (collector-
emitter voltage is positive). It is turned off whenever the gating signal is low or the current drops to zero.
Images:
Attributes:
Parameters Description
Saturation Voltage Saturation voltage Vce_sat of the IGBT, in V
Transistor Resistance Transistor on resistance, in Ohm
Diode Forward Voltage Anti-parallel diode forward threshold voltage, in V
4.2.8 GTO
A GTO switch is a symmetrical device with both forward-blocking and reverse-blocking capabilities.
Images:
Attributes:
Parameters Description
Voltage Drop GTO conduction voltage drop, in V
Initial Position Initial switch position flag (0: off; 1: on)
Current Flag Switch current flag (0: no display; 1: display)
Attributes:
Parameters Description
Initial Position Initial switch position flag (for single switch only).
Initial Position (for Phase A/B/ Initial switch position for Phase A, or B, or C.
C)
Switch Position Switch position can be either "On" or "Off" (for push button switch only)
Current Flag Switch current flag (for single switch only)
Current Flag (for Phase A/B/C) Switch current flag for Phase A, or B, or C.
In the 3-phase switch image, the phase with a dot is Phase A.
For single and 3-phase bi-directional switches, they are on when the gating signal is high, and are off when the
gating signal is low, regardless of the voltage bias conditions. For the push button switch, the switch position is
set directly from the parameter input.
Parameters Description
Current Gain beta Transistor current gain , defined as: =Ic/Ib
Bias Voltage Vr Forward bias voltage, in V, between base and emitter for the npn transistor, or
between emitter and base for the pnp transistor.
Vce,sat [or Vec,sat for pnp] Saturation voltage, in V, between collector and emitter for the npn transistor, and
between emitter and collector for the pnp transistor.
Attributes for MOSFET (3-state) and p-MOSFET (3-state) Level-1 Model:
Parameters Description
On Resistance On resistance RDS(on) of the MOSFET, in Ohm
Threshold Voltage Vgs(th) Threshold voltage of the gate-to-source voltage, in V, beyond which the MOSFET
starts to conduct.
Transconductance gm Transconductance gm of the MOSFET.
Diode Forward Voltage Forward voltage of the anti-parallel diode, in V
Diode Resistance On resistance of the anti-parallel diode, in Ohm
Attributes for MOSFET (3-state) and p-MOSFET (3-state) Level-2 Model:
Parameters Description
Vbreakdown (drain- Maximum drain-source voltage without causing avalanche breakdown, in V.
source)
On Resistance On resistance RDS(on) of the MOSFET, in Ohm
Gate Threshold Voltage Threshold voltage of the gate-to-source voltage, in V, above which, the MOSFET
Vgs(th) starts to conduct.
Internal Gate Resistance Internal gate resistance, in Ohm.
Transconductance gm Transconductance gm of the MOSFET.
Capacitance Cgs Gate-to-source intrinsic capacitance, in F.
Capacitance Cgd Gate-to-drain intrinsic capacitance, in F.
NPN_1
NPN_1
Attributes:
Parameters Description
Frequency Operating frequency of the switch or switch module connected to the gating block, in
Hz
No. of Points Number of switching points (for the Gating Block element only)
Switching Points Switching points, in deg. If the frequency is zero, the switching points is in second. (for
the Gating Block element only)
File for Gating Name of the file that stores the gating table (for the Gating Block (file) element only)
Table
The number of switching points is defined as the total number of switching actions in one period. Each turn-on
or turn-off action is counted as one switching point. For example, if a switch is turned on and off once in one
cycle, the number of switching points will be 2.
For the Gating Block (file) element, the file for the gating table must be in the same directory as the schematic
file. The gating table file has the following format:
n
G1
G2
... ...
Gn
where G1, G2, ..., Gn are the switching points.
Example:
Assume that a switch operates at 2000 Hz and has the following gating pattern in one period:
The specification of the Gating Block element for this switch will be:
Frequency 2000.
No. of Points 6
Switching Points 35. 92. 175. 187. 345. 357.
The gating pattern has 6 switching points (3 pulses). The corresponding switching angles are 35o, 92o, 175o,
187o, 345o, and 357o, respectively.
Frequency 2000.
File for Gating Table test.tbl
The file test.tbl will contain the following:
6
35.
92.
175.
187.
345.
357.
A-
A-
4 2 A- 4 2
A- DC- DC- DC-
DC-
Ct
Attributes:
Parameters Description
Diode Threshold Voltage Threshold voltage of the diode, or forward voltage drop of the thyristor, in V
or Voltage Drop
Diode Resistance On resistance of the diode, in Ohm (for diode bridges only)
Init. Position_i Initial position for Switch i
Current Flag_i Current flag for Switch i
Node Ct at the bottom of the thyristor module is the gating control node for Switch 1. For the thyristor module,
only the gating signal for Switch 1 needs to be specified. The gating signals for other switches will be derived
internally in the program.
Similar to the single thyristor switch, a thyristor bridge can also be controlled by either a gating block or an
alpha controller, as shown in the following examples.
Examples: Control of a Thyristor Bridge
The gating signal for the circuit on the left is specified through a gating block, and the gating signal for the
circuit on the right is provided through an alpha controller. A major advantage of the alpha controller is that the
delay angle alpha of the thyristor bridge, in deg., can be directly controlled.
C 3
C
Ct A6 6
Ct
DC-
CSI3 CSI3
DC+
DC+ A 1 3 5
Ct
B A
B
C
DC- C
4 6 2
Ct DC-
Parameters Description
Diode Threshold Voltage Threshold voltage drop of the diode, in V
Diode Resistance On resistance of the diode, in Ohm
Init. Position_i Initial position for Switch i
Current Flag_i Current flag for Switch i
Attributes for thyristor bridges:
Parameters Description
Voltage Drop Conduction voltage drop of the thyristor, in V
Init. Position_i Initial position for Switch i
Current Flag_i Current flag for Switch i
Attributes for VSI3 bridge (with MOSFET switches):
Parameters Description
On Resistance On resistance of the MOSFET transistor, in Ohm
Diode Threshold Voltage Threshold voltage drop of the diode, in V
Diode Resistance On resistance of the diode, in Ohm
Init. Position_i Initial position for Switch i
Current Flag_i Current flag for Switch i
Attributes for VSI3 bridge (with IGBT switches):
Parameters Description
Saturation Voltage Saturation voltage Vce_sat of the IGBT transistor, in Ohm
Transistor Resistance On resistance of the IGBT transistor, in Ohm
Diode Threshold Voltage Threshold voltage drop of the anti-parallel diode, in V
Diode Resistance On resistance of the anti-parallel diode, in Ohm
Init. Position_i Initial position for Switch i
Current Flag_i Current flag for Switch i
Attributes for the CSI3 bridge:
Parameters Description
Voltage Drop Forward conduction voltage drop of the switch, in V
Resistance Forward on resistance of the switch, in Ohm
Init. Position_i Initial position for Switch i
Current Flag_i Current flag for Switch i
Similar to single-phase modules, only the gating signal for Switch 1 need to be specified for three-phase
modules. Gating signals for other switches will be automatically derived. For the 3-phase half-wave thyristor
bridge, the phase shift between two consecutive switches is 120o. For all other bridges, the phase shift is 60o.
Thyristor bridges can be controlled by an alpha controller. Similarly, voltage/current source inverters can be
The thyristor circuit on the left uses an alpha controller. For a three-phase circuit, the zero-crossing of the
voltage Vac corresponds to the moment when the delay angle alpha is equal to zero. This signal is used to
provide synchronization to the controller.
The circuit on the right uses a PWM lookup table controller. The PWM patterns are stored in a lookup table in a
text file. The gating pattern is selected based on the modulation index. Other inputs of the PWM lookup table
controller include the delay angle, the synchronization, and the enable/disable signal. A detailed description of
the PWM lookup table controller is given in the Switch Controllers section.
4.3 Transformers
Np Ns Np Ns
The winding with the larger dot is the primary, and the other winding is the secondary.
Attributes:
Parameters Description
Np (primary) No. of turns of the primary winding
Ns (secondary) No. of turns of the secondary winding
Since the turns ratio is equal to the ratio of the rated voltages, the number of turns can be replaced by the rated
voltage at each side.
In the images, p refers to primary, s refers to secondary, and t refers to tertiary. The winding with the largest dot
is the primary winding or first primary winding. For the multiple winding transformers, the sequence of the
windings is from the top to the bottom. For the transformers with 2 or 3 windings, the attributes are as follows.
Attributes:
Parameters Description
Rp (primary); Resistance of the primary/secondary/tertiary winding, in Ohm
Rs (secondary);
Rt (tertiary)
Lp (pri. leakage); Leakage inductance of the primary/secondary/tertiary winding, in H (seen from the
Ls (sec. leakage); primary)
Lt (ter. leakage)
Lm (magnetizing) Magnetizing inductance, in H
Np (primary); No. of turns of the primary/secondary/tertiary winding
Ns (secondary);
Nt (tertiary)
All the resistances and inductances are referred to the primary winding side. If there are multiple primary
windings, they are referred to the first primary winding side.
For a transformer with more than 1 primary winding or more than 3 secondary windings, the attributes are as
follows.
Attributes:
Parameters Description
Rp_i (primary i); Resistance of the ith primary/secondary/tertiary winding, in Ohm
Rs_i (secondary i)
Primary Lm Secondary
Ideal
In the circuit, Rp and Rs are the primary and secondary winding resistances; Lp and Ls are the primary and
secondary winding leakage inductances; and Lm is the magnetizing inductance. All the values are referred to the
primary side.
Example:
A single-phase two-winding transformer has a winding resistance of 0.002 Ohm and leakage inductance of 1
mH at both the primary and the secondary side (all the values are referred to the primary). The magnetizing
inductance is 100 mH, and the turns ratio is Np:Ns = 220:440. The transformer will be specified as:
Rp (primary) 2m
Rs (secondary) 2m
Lp (primary) 1m
Ls (secondary) 1m
Lm (magnetizing) 100m
Np (primary) 220
Ns (secondary) 440
b a+
B b B b B a-
b+ Winding 2
C c C c C c b-
c+
N n N c-
Y/Y/D Y/D/D 2-winding (unconnected) Winding 1
n Winding 3
a+ (primary)
a a A+
b A- a- A+
A b A b+ A-
B+
c c B+
B B B- b- B- Winding 4
aa aa C+ c+ C+
C C bb C- c- C-
bb
cc cc
Winding 5
N N
B b B b B b B b
C c C c C c C c
N N
Connection diagrams:
a a a a
b b b b
c c c c
Np Ns1 Ns2 Np Ns1 Ns2 Np Ns1 Ns2 Np Ns1 Ns2
Parameters Description
Rp (primary); Rs (secondary); Resistance of the primary/secondary/tertiary winding, in Ohm
Rt (tertiary)
Lp (pri. leakage); Ls (sec. leakage); Leakage inductance of the primary/secondary/tertiary winding, in H
Lt (ter. leakage)
Lm (magnetizing) Magnetizing inductance, in H (seen from the primary side)
Np (primary); Ns (secondary); No. of turns of the primary/secondary/tertiary winding
Nt (tertiary)
Parameters Description
Rp (primary); Rs1 (secondary 1); Resistance of the primary/1st secondary/2nd secondary winding, in
Rs2 (secondary 2) Ohm
In the images, P refers to primary, S refers to secondary, and T refers to tertiary. All resistances and inductances
are referred to the primary or the first primary winding side.
Three-phase transformers are modeled in the same way as single-phase transformers.
For the 3-phase phase-shifting transformers, let the angle be the phase difference between the secondary line
voltage Vab and the primary line voltage VAB. The relationship between the angle and the winding turns are:
For the Y/Z1 transformer:
Ns2 o
sin ( 30 ) o o
--------------------- = -------------------------------
-, 0 30
Ns1 + Ns2 sin ( 30 + )
o
Np 1 V AB
--------------------- = ----------------------------------- ------
-
Ns 1 + Ns 2 2 sin ( 30 + ) V ab
o
Np 1 V AB
--------------------- = -------------------------------------- ------
-
Ns 1 + Ns 2 2 sin ( 30 + ) V ab
o
Np 3 V AB
--------------------- = -------------------------------------
- ------
-
Ns 1 + Ns 2 2 sin ( 60 ) V ab
o
Np 3 V AB
--------------------- = ------------------------ -------
Ns1 + Ns2 2 sin ( ) V ab
A+ a+
A- a-
B+ b+
B- b-
C+ c+
C- c-
en
Attributes:
Parameters Description
Operating Voltage Transformer primary side rms voltage per phase under the circuit operating
(phase) condition, in V
Operating Frequency Operating frequency of the transformer, in Hz
Rp (primary) Resistance of the primary winding, in Ohm
Lp (pri. leakage) Leakage inductance of the primary winding, in H
Rs (secondary) Resistance of the secondary winding, in Ohm (referred to the primary side)
Ls (sec. leakage) Leakage inductance of the secondary winding, in H (referred to the primary side)
Rm (core losses) Resistance that represents transformer core losses, in Ohm
Im vs. Lm Values of the magnetizing current Im (in A) vs. the magnetizing inductance Lm (in
H) that defines the magnetizing characteristics seen from the primary winding
Phase A Residual Flux Phase A residual flux, in per unit value
Phase B Residual Flux Phase B residual flux, in per unit value
Phase C Residual Flux Phase C residual flux, in per unit value
Np (primary) Number of turns of the primary winding
Ns (secondary) Number of turns of the secondary winding
4.4.1 Winding
A winding element provides the interface between the electric circuit and the magnetic equivalent circuit.
Image:
M1
E1
E2
M2
Attributes:
Parameters Description
Number of Turns No. of turns of the winding
Winding Resistance Winding resistance, in Ohm
This element represents a winding on a magnetic core. The two electric nodes (E1 and E2) are connected to an
electric circuit, while the two magnetic nodes (M1 and M2) are connected to other magnetic elements (such as
leakage flux path, air gap, and magnetic core).
M1 M2
Attributes:
Parameters Description
Inductance Factor AL Inductance factor AL, defined as the inductance per turn squared
Resistance for Losses Resistance R, in Ohm, that represents the losses due to the leakage flux.
Current Flag Display flag of the current that flows through the resistor R
The resistance R represents the losses due to the leakage flux.
Assuming that the mmf (magnetomotive force) applied across the leakage flux path is F, the electric equivalent
circuit of the leakage flux path is as follows:
+ i
AL
F
R
-
The mmf, in the form of a voltage source, applies across the capacitor (the capacitance is AL) and the resistor R.
Let the current flowing through this branch be i, and the rms value be Irms, the relationship between the losses
due to the leakage flux and the resistance R is:
2
P loss = I rms R
M1 M2
The input parameters of the air gap can be defined in two ways. One is to define the air gap length and the cross
section area, and the other is to define the inductance factor AL. They are as follows.
Attributes:
For the element Air Gap:
Parameters Description
Air Gap Length The length of the air gap, lg, in m
Parameters Description
Inductance Factor AL Inductance factor AL, defined as the inductance per turn squared
Resistance for Losses Resistance R, in ohm, that represents the losses due to the air gap fringing effect
Current Flag Display flag of the current that flows through the resistor R
The resistance R represents the losses due to the air gap fringing effect. Assuming that the mmf (magnetomotive
force) applied across the air gap is F, the electric equivalent circuit of the air gap is as follows:
+ i
AL
F
R
-
The mmf, in the form of a voltage source, applies across the capacitor (the capacitance has the value of the
inductance factor AL) and the resistor R. For the element Air Gap, the inductance factor can be calculated from
the air gap length and the cross section area as:
o Ac
A L = --------------
-
lg
where o= 4107.
The losses on the resistor represents the losses due to the fringing effect, which can be expressed as:
2
P loss = I rms R
where Irms is the rms value of the current i flowing through the resistor.
M1 M2
Attributes:
Parameters Description
Inductance Factor AL Inductance factor AL of the core, defined as the inductance per turn squared
If the length of the core is Llength and the cross section area is Ac, the inductance factor AL is expressed as:
o r Ac
A L = -------------------------
L length
Attributes:
Parameters Description
Inductance Factor AL Inductance factor AL of the core, defined as the inductance per turn squared
Resistance for Losses Resistance R, in Ohm, that represents the core losses
Coefficient phi_sat Coefficient sat for the core B-H curve, in Weber
Coefficient K1 Coefficient K1 for the core B-H curve
Coefficient Kexp1 Coefficient Kexp1 for the core B-H curve
Coefficient K2 Coefficient K2 for the core B-H curve
Coefficient Kexp2 Coefficient Kexp2 for the core B-H curve
Initial Flux phi_o Initial flux o of the core, in Weber
Current Flag Display flag of the electric current that flows through the resistor R. If the rms value
of the current is Irms, the core losses can be calculated as: Pcore_loss = Irms2 * R.
In the element image, the nodes M1 and M2 are the two nodes that connect the core to other magnetic elements
(such as winding, flux leakage path, air gap, etc.). The node marked with a dot is Node M2. Node C1 is a control
output node, which shows the flux (in Weber) flowing through the core, from Node M2 to M1.
The initial flux of the core, o, is the flux flowing from Node M2 to M1 at the beginning, corresponding to the
field intensity H = 0. The initial flux density Bo can be calculated as: Bo = o/Ac (where Ac is the core cross
section area), and the initial mmf Fo is: Fo = o/AL.
A differential voltage probe connected between Node M2 to M1 will measure the mmf (in ampere*turn) applied
to the core.
The coefficients sat, K1, Kexp1, K2, and Kexp2 are used to fit the B-H curve of an actual magnetic material. A
good initial guess of sat is the maximum flux of the B-H curve in deep saturation. To calculate this flux,
multiply the corresponding flux density B by the cross section area of the core. Coefficient K1 usually varies
between 0.7 and 1, depending on the core material. Coefficient Kexp1 mainly affects the rate of the core
saturation, and is in the range between 10 and 200 (10 for low permeability ferrite, and 200 for metglas).
The coefficients K2 and Kexp2 are used in very rare occasions, such as for ferroresonant regulators. They are
normally set as follows to keep them from affecting the B-H curve:
K2 > 2
Kexp2 > 20
A tutorial on how to define the core coefficients is provided in the "doc" sub-folder in the PSIM directory. Also,
a program is provided to plot the B-H curve of the core given a set of coefficients. To launch this program, in
PSIM, go to Utilities >> B-H Curve.
V+ Ro
gnd
where
V+; V- - noninverting and inverting input voltages
Vo - output voltage
Ao - op. amp. dc gain (Ao is set to 100,000.)
Ro - output resistance (Ro is set to 80 Ohms)
Attributes:
Parameters Description
Voltage Vs+ Upper voltage source level of the op. amp.
Voltage Vs- Lower voltage source levels of the op. amp.
The difference between the element Op. Amp. and Op. Amp. (ground) or Op. Amp. (ground/inverted) is that, for
the Op. Amp. element, the reference ground of the op. amp. model is connected to the power ground, whereas
for Op. Amp. (ground) or Op. Amp. (ground/inverted), the reference ground node of the model is accessible and
can be floating.
Note that the image of an op. amp. is similar to that of a comparator. For the op. amp., the inverting input is at
the upper left and the noninverting input is at the lower left. For the comparator, it is the opposite.
Example: A Boost Power Factor Correction Circuit
The figure below shows a boost power factor correction circuit. The PI regulators of both the inner current loop
and the outer voltage loop are implemented using op. amp.
Attributes:
Parameters Description
Input Resistance Rin Input resistance of the op. amp., in Ohm
DC Gain Ao DC gain of the op. amp.
Unit Gain Frequency Frequency at which the op. amp. gain is 1, in Hz.
Output Resistance Ro Output resistance of the op. amp., in Ohm
Maximum Output Maximum current that the op. amp. output can provide, in A
Current
Voltage Vs+ Upper voltage source level of the op. amp.
Voltage Vs- Lower voltage source levels of the op. amp.
The figure below shows the definition of the dc gain Ao and the unit gain frequency.
0
f (Hz)
Cathode
Reference
Anode
The TL431 regulator maintains the reference node voltage to be around 2.495V. In this model, the limited
bandwidth of the TL431 is modelled. That is, the ac frequency response of the TL431 model resembles the
following figure from the datasheet.
Test Circuit
4.5.3 Opto-Coupler
The image and parameters of a opto-coupler are shown below.
Image:
id ic
Parameters Description
Current Transfer Ratio The current transfer ratio CTR between the transistor current ic and the diode
current id, that is, CTR = ic/id.
Diode Resistance Diode resistance Rd, in Ohm. Note that the resistance must be greater than 0.
Diode Threshold The forward diode threshold voltage Vd_th, in V
Voltage
Transistor Vce_sat The transistor saturation voltage Vce_sat, in V
Transistor-side Capacitance Cp across the collector and emitter of the transistor, in F
Capacitance
These parameters can either be read directly from the manufacturer datasheet, or can be calculated from the
datasheet information. Take the Motorola optocoupler 4N25 as an example. From the datasheet, we can read the
current transfer ratio and the transistor saturation voltage as:
CTR = 70% = 0.7; Vce_sat = 0.15 V
From the "LED Forward Voltage versus Forward Current" characteristics on the datasheet, one can calculate the
diode resistance Rd and the threshold voltage Vd_th. We assume that the diode forward voltage is approximated
by a straight line. That is,
vd = Vd_th + Rd * id
With TA = 25oC, we can obtain two readings from the curve: id = 10 mA and vd = 1.15 V; id = 40 mA and vd =
1.3 V. Based on these two points, we can calculate Rd and Vd_th as: Rd = 5 Ohm; Vd_th = 1.1 V.
The optocoupler model does not take into account the delay, but it does model the turn-on/turn-off transient
through the capacitor Cp across the collector and emitter of the transistor. The capacitance value can be obtained
from the fall time tfall and the switching time test condition as:
t fall
C p = ------------------
2.2 R L
where RL is the load resistor in the switching time test circuit. For 4N25, for example, from the datasheet, we
can obtain: tfall = 1.3 s, and RL = 100 Ohm. The capacitance is calculated as: Cp = 6 nF.
The output of the dv/dt block is equal to the derivative of the input voltage versus time. It is calculated as:
V in ( t ) V in ( t t )
V o = --------------------------------------------------
t
where Vin(t) and Vin(t-t) are the input values at the current and previous time step, and t is the simulation time
step.
Attributes:
Parameters Description
Rated Coil Voltage Rated voltage of the relay coil, in V
Coil Resistance Resistance of the coil, in Ohm
Operate Voltage Voltage at which the relay will operate, in V
Release Voltage Voltage at which the relay will return to the default position, in V
Operate Time Time from the moment the operate voltage is reached to the moment the switches
operate, in sec.
Release Time Time from the moment the release voltage is reached to the moment the switches
return to the default positions, in sec.
The relays have one normally open switch and one normally closed switch. When a dc voltage is applied to the
relay coil in the polarity as shown in the image, if the voltage reaches the Operate Voltage, after a time delay
defined by the Operate Time, the NO switch will be closed, and the NC switch will be opened.
When the coil voltage reduces to the Release Voltage, after a time delay defined by the Release Time, the two
switches will return to their default positions.
The system consists of two induction machines, IM1 and IM2, connected back-to-back. One operates as a
motor, and the other as a generator. From the point of view of the first machine IM1, the mechanical equation
can be written as:
The reference direction of each element and the reference direction of the overall mechanical system determine
how the element interacts with the mechanical system.
For example, if the reference direction of a machine is along the same direction as the reference direction of the
mechanical system, the developed torque of the machine will contribute to the shaft rotation in the positive
direction. However, if the reference direction of the machine is opposite to that of the mechanical system, the
developed torque will contribute to the shaft rotation in the negative direction.
In the two-machine example above, using the notation of the "reference direction", if we define the machine
IM1 as the master unit, the reference direction of the overall mechanical system will be from left to right, as
shown below. Based on this direction, the machine IM1 will be along the reference direction, and the machine
IM2 will be opposite to the reference direction. This leads to the equivalent circuit of the mechanical system as
Equivalent Circuit
Reference direction
Master Unit (J1+J2)*dWm/dt = Tem1 - Tem2
Similarly, if we define the machine IM2 as the master unit, the reference direction of the overall mechanical
system will be from right to left, as shown below. Based on this direction, the machine IM1 will be opposite to
the reference direction, and the machine IM2 will be along the reference direction. This leads to the equivalent
circuit of the mechanical system as shown on the right.
Master Unit
Equivalent Circuit
The following shows another mechanical system with sensors and loads connected in different ways.
In this mechanical system, the machine on the left is the master unit. The reference direction of the mechanical
system is from left to the right along the mechanical shaft. Comparing this direction with the reference direction
of each element, Load 1, Speed Sensor 1, and Torque Sensor 1, will be along the reference direction, and Load
2, Speed Sensor 2, and Torque Sensor 2 will be opposite to the reference direction of the mechanical system.
Therefore, if the speed of the machine is positive, Speed Sensor 1 reading will be positive, and Speed Sensor 2
reading will be negative.
Similarly, the two constant-torque mechanical loads, with the amplitudes of TL1 and TL2, interact with the
machine in different ways. Load 1 is along the reference direction, and the loading torque of Load 1 to the
master machine will be TL1. On the other hand, Load 2 is opposite to the reference direction, and the loading
torque of Load 2 to the machine will be -TL2.
Attributes:
Parameters Description
Rs (stator) Stator winding resistance, in Ohm
Ls (stator) Stator winding leakage inductance, in H
Rr (rotor) Rotor winding resistance referred to the stator side, in Ohm
Lr (rotor) Rotor winding leakage inductance referred to the stator side, in H
Lm (magnetizing) Magnetizing inductance, in H
Ns/Nr Turns Ratio Stator and rotor winding turns ratio (for wound-rotor machine only)
No. of Poles Number of poles P of the machine (an even integer)
Moment of Inertia Moment of inertia J of the machine, in kg*m2
Torque Flag Output flag for internal torque Tem
Master/Slave Flag Master/slave flag of the machine (1: master; 0: slave)
All the parameters are referred to the stator side.
For more details on the definition and use of the master/slave flag, refer to Section 4.6.1.
The models of the squirrel-cage induction machine with and without the neutral are the same internally.
v a, s v a, r i a, s i a, r
v abc, s = v b, s v abc, r = v b, r i abc, s = i b, s i abc, r = i b, r
v c, s v c, r i c, s i c, r
For squirrel-cage machines, va,r = vb,r = vc,r= 0. The parameter matrices are defined as:
Rs 0 0 Rr 0 0
Rs = 0 Rs 0 Rr = 0 Rr 0
0 0 Rs 0 0 Rr
M sr M sr M sr M sr
L s + M sr ------- ------- L r + M sr ------- -------
2 2 2 2
M M M M
Ls = sr
------
- L s + M sr sr
------
- Lr = sr
------
- L r + M sr sr
------
-
2 2 2 2
M sr M sr M sr M sr
------
- ------
- L s + M sr ------
- ------
- L r + M sr
2 2 2 2
2 2
cos cos + ------ cos ------
3 3
M sr = M sr cos 2
------ cos
2
cos + ------
3 3
2 2
cos + ------ cos ------ cos
3 3
where Msr is the mutual inductance between the stator and rotor windings, and is the mechanical angle. The
mutual inductance is related to the magnetizing inductance as:
3
L m = --- M sr
2
The mechanical equation is expressed as:
d
J --------m- = T em T L
dt
where the developed torque Tem is defined as:
P T d
T em = --- i abc, s ------ M sr i abc, r
2 d
For a symmetrical squirrel-cage induction machine, the steady state equivalent circuit is shown below. In the
figure, s is the slip.
Lm Rr(1-s)/s
Speed Torque
Sensor Sensor
Speed
Tem
Tload
SPWM
3-phase currents
Parameters Description
Rs (stator) Stator winding resistance, in Ohm
Ls (stator) Stator winding leakage inductance, in H
Rr (rotor) Rotor winding resistance referred to the stator side, in Ohm
Lr (rotor) Rotor winding leakage inductance referred to the stator side, in H
Ns/Nr Turns Ratio Stator and rotor winding turns ratio (for wound-rotor machine only)
No. of Poles Number of poles P of the machine (an even integer)
Moment of Inertia Moment of inertia J of the machine, in kg*m2
Torque Flag Output flag for internal torque Tem
Master/Slave Flag Master/slave flag of the machine (1: master; 0: slave)
Im vs. Lm (Im1,Lm1) ... Characteristics of the magnetizing current Im versus the magnetizing inductance
[(Im1,Lm1) (Im2,Lm2) ...] where Im is in A and Lm is in H.
2 2
1 ---
1 1
--- cos cos + ------ cos ------
2 2 3 3
1 1 2 2
abc, s = M sr --- 1 --- i abc, s + M sr cos ------ cos cos + ------ i abc, r
2 2 3 3
1 1 2 2
--- ---
2 2
1 cos + ------ cos ------ cos
3 3
2 2
cos cos ------ cos + ------ 1 ---
1 1
---
3 3 2 2
2 2 + M sr 1 1
abc, s = M sr cos + ------ cos cos ------ i abc, s --- 1 --- i abc, r
3 3 2 2
2 2 1 1
cos ------ cos + ------ cos --- ---
2 2
1
3 3
In this case, the inductance Msr is no longer constant, but a function of the magnetizing current Im. The
magnetizing current Im and inductance Lm are defined by a series of data points in pair, for example, "(0, 0.041)
(5, 0.035) (20, 0.03)". Between two data points, the inductance Lm is interpolated linearly. When the
magnetizing current is less than the minimum or greater than the maximum, the inductance will be clamped to
the value at the first point or the last point.
+
Armature
Shaft Node
Winding
-
+
Field
Winding
-
Attributes:
Parameters Description
Ra (armature) Armature winding resistance, in Ohm
La (armature) Armature winding inductance, in H
Rf (field) Field winding resistance, in Ohm
Lf (field) Field winding inductance, in H
Moment of Inertia Moment of inertia of the machine, in kg*m2
Vt (rated) Rated armature terminal voltage, in V
Ia (rated) Rated armature current, in A
n (rated) Rated mechanical speed, in rpm
If (rated) Rated field current, in A
Torque Flag Output flag for internal torque Tem
Master/Slave Flag The master/slave flag of the machine (1: master; 0: slave)
When the torque flag is set to 1, the internal torque generated by the machine will be saved to the output file for
display.
For more details on the definition and use of the master/slave flag, refer to Section 4.6.1.
The operation of a dc machine is described by the following equations:
di
v t = E a + i a R a + L a -------a-
dt
di
v f = i f R f + L f -------f-
dt
Ea = k m
T em = k i a
d
J ---------m = T em T L
dt
where vt, vf, ia, and if are the armature and field winding voltage and current, respectively; Ea is the back emf,
E a = L af i f m
T em = L af i f i a
where Laf is the mutual inductance between the armature and the field windings. It can be calculated based on
the rated operating conditions as:
( Vt Ia Ra )
L af = ------------------------------
-
If m
Note that the dc machine model assumes magnetic linearity. Saturation is not considered.
Example: A DC Motor with a Constant-Torque Load
The circuit below shows a shunt-excited dc motor with a constant-torque load TL. Since the load is along the
reference direction of the mechanical system, the loading torque to the machine is TL. Also, the speed sensor is
along the reference direction. It will give a positive output for a positive speed.
The simulation waveforms of the armature current and the speed are shown on the right.
Speed
Sensor
Armature current
Constant-
Torque
Load Speed (in rpm)
b Shaft Node
n sa sb sc
6-pulse Hall Effect Position Sensor
Parameters Description
R (stator resistance) Stator phase resistance R, in Ohm
L (stator self ind.) Stator phase self inductance L, in H
Parameters Description
Resistance (phase-phase) Phase-to-phase (or line-to-line) resistance, in Ohm
Inductance (phase-phase) Phase-to-phase (or line-to-line) inductance, in H
Speed Constant Speed constant Kv, defined as the ratio between the speed and the applied voltage,
in rpm/V
Torque Constant Torque constant Kt, defined as the ratio between the generated torque and the
applied current, in N*m/A
No. of Poles P Number of poles P
Moment of Inertia Moment of inertia J of the machine, in kg*m2
No Load Speed The motor speed at no load with the nominal voltage applied, in rpm
Ec = ke _ c m
The coefficients ke_a, ke_b, and ke_c are dependent on the rotor angle r. In this model, an ideal trapezoidal
waveform profile is assumed, as shown below for Phase A. Also shown is the Phase A hall effect sensor signal
Sa.
ke_a
Sa
Kpk
o
180
o r
360
-Kpk
V pk krpm 1
where Kpk is the peak trapezoidal value, in V/(rad./sec.), which is defined as: K pk = ------------------------ --------------------------------- .
2 1000 2 60
Given the values of Vpk/krpm and Vrms/krpm, the angle is determined automatically in the program.
The developed torque of the machine is:
T em = ( E a i a + E b i b + E c i c ) m
The shaft time constant shaft, therefore, reflects the effect of the friction and windage of the machine. Note that
when the shaft time constant is set to 0, the friction term is B*m is ignored.
To better understand the definition of the shaft time constant, we can represent the mechanical equation with the
following equivalent circuit:
m
1/B Tload
Tem J (or shaft/J)
This circuit shows that the shaft time constant shaft is equal to the RC time constant of the resistor 1/B and the
capacitor J. Therefore, the shaft time constant can be measured by the following test:
- Connect the machine to an external mechanical source. With the stator side in open circuit, drive the
machine to a certain speed.
- Remove the mechanical source. The shaft time constant will be equal to the time that it takes the
machine to decelerate to 36.8% of its initial speed.
More Explanation on the Hall Effect Sensor:
A hall effect position sensor consists of a set of hall switches and a set of trigger magnets.
The hall switch is a semiconductor switch (e.g. MOSFET or BJT) that opens or closes when the magnetic field
is higher or lower than a certain threshold value. It is based on the hall effect, which generates an emf
proportional to the flux-density when the switch is carrying a current supplied by an external source. It is
common to detect the emf using a signal conditioning circuit integrated with the hall switch or mounted very
closely to it. This provides a TTL-compatible pulse with sharp edges and high noise immunity for connection to
the controller via a screened cable. For a three-phase brushless dc motor, three hall switches are spaced 120
electrical deg. apart and are mounted on the stator frame.
The set of trigger magnets can be a separate set of magnets, or it can use the rotor magnets of the brushless
motor. If the trigger magnets are separate, they should have the matched pole spacing (with respect to the rotor
magnets), and should be mounted on the shaft in close proximity to the hall switches. If the trigger magnets use
the rotor magnets of the machine, the hall switches must be mounted close enough to the rotor magnets, where
they can be energized by the leakage flux at the appropriate rotor positions.
Example: Defining Brushless DC Motor Parameters from Manufacturer Datasheet
This example illustrates how to define brushless dc motor parameters from manufacturer datasheet. Below is
the information provided on the datasheet of the brushless dc motor Maxon EC-22-16730 (32V, 50W) from
Maxon Motor.
Values at nominal voltage
Nominal voltage (V) 32
No load speed (rpm) 38700
No load current (mA) 327
Brushless DC Motor
Speed
Tem
3-phase currents
b Shaft Node
c
n
field+ field-
Attributes:
Parameters Description
Rs (stator) Stator winding resistance, in Ohm
Ls (stator) Stator leakage inductance, in H
Ldm (d-axis mag. ind.) d-axis magnetizing inductance, in H
Lqm (q-axis mag. ind.) q-axis magnetizing inductance, in H.
Rf (field) Field winding resistance referred to the stator side, in Ohm
Lfl (field leakage ind.) Field winding leakage inductance referred to the stator side, in H
Rdr (damping cage) Rotor damping cage d-axis resistance referred to the stator side, in Ohm
Ldrl (damping cage) Rotor damping cage d-axis leakage inductance referred to the stator side, in H
Rqr (damping cage) Rotor damping cage q-axis resistance referred to the stator side, in Ohm
Lqrl (damping cage) Rotor damping cage q-axis leakage inductance referred to the stator side, in H
Ns/Nf (effective) Stator-field winding effective turns ratio
Number of Poles P Number of Poles P
Moment of Inertia Moment of inertia J of the machine, in kg*m2
Torque Flag Output flag for internal developed torque Tem
Master/slave Flag Master/slave flag of the machine (1: master; 0: slave)
All the parameters are referred to the stator side.
For more details on the definition and use of the master/slave flag, refer to Section 4.6.1.
The equations of the synchronous machine can be expressed as follows:
d-
V = R I + -----
dt
where
T
V = va vb vc vf 0 0 I = i a i b i c i f i dr i qr
L 11 L 12
L = T
L 12 L 22
and
L 2 L 2
L s + L o + L 2 cos ( 2 r ) ----o- + L 2 cos 2 r ------ ----o- + L 2 cos 2 r + ------
2 3 2 3
L L
L 11 = ----o- + L 2 cos 2 r 2 2
------ L s + L o + L 2 cos 2 r + ------ ----o- + L 2 cos ( 2 r )
2 3 3 2
Lo 2 Lo 2
----- + L 2 cos 2 r + ------ ----- + L 2 cos ( 2 r ) L s + L o + L 2 cos 2 r ------
2 3 2 3
L f L fdr 0
L 22 = L fdr L dr 0
0 0 L qr
b Shaft Node
Attributes:
Parameters Description
Rs (stator resistance) Stator winding resistance, in Ohm
Ld (d-axis ind.) Stator d-axis inductance, in H
Lq (q-axis ind.) Stator q-axis inductance, in H.
The d-q coordinate is defined such that the d-axis passes through the center of the
magnet, and the q-axis is in the middle between two magnets. The q-axis is leading
the d-axis.
Vpk / krpm Peak line-to-line back emf constant, in V/krpm (mechanical speed).
The value of Vpk/krpm should be available from the machine data sheet. If this data
is not available, it can be obtained through an experiment by operating the machine
as a generator at 1000 rpm and measuring the peak line-to-line voltage.
No. of Poles P Number of poles P
Moment of Inertia Moment of inertia J of the machine, in kg*m2
Shaft Time Constant Shaft time constant shaft
Torque Flag Output flag for internal developed torque Tem
Master/slave Flag Master/slave flag of the machine (1: master; 0: slave)
The node assignments of the image are: Nodes a, b, and c are the stator winding terminals for Phase a, b, and c,
respectively. The stator windings are Y connected, and Node n is the neutral point. The shaft node is the
connecting terminal for the mechanical shaft. They are all power nodes and should be connected to the power
circuit.
For more details on the definition and use of the master/slave flag, refer to Section 4.6.1.
The equations of the permanent-magnet synchronous machine are:
where va, vb, vc, and ia, ib, and ic, and a, b, c are the stator phase voltages, currents, and flux linkages,
respectively, and Rs is the stator phase resistance. The flux linkages are further defined as:
cos ( r )
a L aa L ab L ac ia 2
cos r ------
b = L ba L bb L bc i b + pm 3
c L ca L cb L cc ic 2
cos r + ------
3
where r is the rotor electrical angle, and pm is the peak stator phase flux linkage defined as:
60 V pk krpm
pm = --------------------------------------
-
3 P 1000
where P is the number of poles.
The stator self and mutual inductances are rotor position dependent, and are defined as:
L aa = L s + L o + L 2 cos ( 2 r )
2
L bb = L s + L o + L 2 cos 2 r + ------
3
2
L cc = L s + L o + L 2 cos 2 r ------
3
L 2
L ab = L ba = ----o- + L 2 cos 2 r ------
2 3
L 2
L ac = L ca = ----o- + L 2 cos 2 r + ------
2 3
L
L bc = L cb = ----o- + L 2 cos ( 2 r )
2
where Ls is the stator leakage inductance. The variables in the abc frame can be converted to the dq0 frame
using the following transformation:
2 2
cos ( r ) cos r ------ cos r + ------
3 3
xd xa
2
xq = --- sin ( ) sin 2
------ sin r + 2
------ x b
3 r r 3 3
x0 xc
1 1 1
--- --- ---
2 2 2
2 2
sin ( 2 r ) sin 2 r ------ sin 2 r + ------ sin ( r )
3 3
ia 2
P
= --- L 2 i a i b i c sin 2 r 2 2 P sin r ------
------ sin 2 r + ------ sin ( 2 r ) i b --2- pm i a i b i c 3
2 3 3
ic 2
2 2 sin r + ------
sin 2 r + ------ sin ( 2 r ) sin 2 r ------ 3
3 3
where B is a coefficient, Tload is the load torque, and P is the number of poles. The coefficient B is calculated
from the moment of inertia J and the shaft time constant shaft as below:
J
B = --------------
shaft
b Shaft Node
c
n
Attributes:
Parameters Description
Rs (stator resistance) Stator winding resistance, in Ohm
Ls (stator leakage ind.) Stator d-axis inductance, in H
L d = L s + L dm
L q = L s + L qm
where Ls is the stator leakage inductance. Since Ls is normally very small, Ld can be considered equivalent to
Ldm, and Lq can be considered equivalent to Lqm.
The Transformation Flag defines the transformation convention between the abc frame and the dq frame. When
the Transformation Flag is 0:
2 2
cos ( r ) cos r ------ cos r + ------ ia
Id 2 3 3
= --- ib
Iq 3 2 2
sin ( r ) sin r ------ sin r + ------ ic
3 3
2 2
Im = Id + Iq
m = atan 2 ( I q, I d )
2 2
cos ( r ) cos r ------ cos r + ------ ia
Id 2 3 3
= --- ib
Iq 3 2 2
sin ( r ) sin r ------ sin r + ------ ic
3 3
m = atan 2 ( I d, I q )
b+
b- Shaft Node
c+
c-
Attributes:
Parameters Description
Resistance Stator phase resistance R, in Ohm
Inductance Lmin Minimum phase inductance Lmin, in H
Inductance Lmax Maximum phase inductance Lmax, in H
Theta_min (deg.) Duration of the interval min where the inductance is at the minimum, in deg.
Theta_max (deg.) Duration of the interval max where the inductance is at the maximum, in deg.
Stator Pole Number Number of stator poles (teeth)
Rotor Pole Number Number of rotor poles (teeth)
Moment of Inertia Moment of inertia J of the machine, in kg*m2
Torque Flag Output flag for internal torque Tem
Master/Slave Flag Master/slave flag of the machine (1: master; 0: slave)
For more details on the definition and use of the master/slave flag, refer to Section 4.6.1.
The node assignments are: Nodes a+, a-, b+, b-, c+, c-, etc. are the stator winding terminals for Phase a, b, c, d,
and e, respectively. The shaft node is the connecting terminal for the mechanical shaft. They are all power nodes
and should be connected to the power circuit.
The equation of the switched reluctance machine for one phase is:
d(L i)
v = i R + ------------------
dt
where v is the phase voltage, i is the phase current, R is the phase resistance, and L is the phase inductance. The
phase inductance L is a function of the rotor angle (in electrical deg.), as shown in the following figure. The
inductance profile is half-wave symmetrical. That is, it repeats itself after 180o.
Lmin
rf min
The rotor angle is defined such that, when the stator and the rotor teeth are completely aligned, = 0. The value
of the inductance can be in either rising stage, flat-top stage, falling stage, or flat-bottom stage.
If we define the constant k as the rising slope of the inductance from Lmin to Lmax (which is the same as the
falling slope from Lmax to Lmin), and the interval of the rising/falling slope is rf, we can express the inductance
L as a function of the rotor angle (from 0 to 180o) as follows:
max
L = L max for 0 ---------
2
-
L = L max k ---------
max max max
2
-
for ---------
2
- ---------
2
- + rf
max
max
L = L min for ---------
2
- + rf ---------
2
- + rf + min
L = L min + k ---------
- rf min for ---------
max max max
- + rf + min ---------
- + 2 rf + min
2 2 2
max
L = L max for ---------
2
- + 2 rf + min 180 o
b+ Shaft Node
b-
c+
c-
Attributes:
Parameters Description
Resistance Stator phase resistance R, in Ohm
Stator Pole Number Number of stator poles (teeth)
Rotor Pole Number Number of rotor poles (teeth)
Inductance Table File The file that stores the 2-dimensional table of the inductance versus the rotor
position and the current.
Moment of Inertia Moment of inertia J of the machine, in kg*m2
Torque Flag Output flag for internal torque Tem
Master/Slave Flag Master/slave flag of the machine (1: master; 0: slave)
The inductance table file has the following format:
m, n
1, 2, ..., m
I1, I2, ..., In
L1,1, L1,2, ..., L1,n
L2,1, L2,2, ..., L2,n
... ... ...
Lm,1, Lm,2, ..., Lm,n
where m is the number of rows and n is the number of columns; 1 to m is the row vector for the rotor position
(in deg.), and I1 and In is the column vector for the phase current (in A); and Li,j is the inductance value (in H) at
Row i and Column j. For example, when the rotor position is 2 and the phase current is I1, the inductance is
L2,1.
Because of the half-wave symmetry, only inductance values from 0 to 180o need to be provided. Note that the
row vector and column vector must be monotonically increasing.
The machine equations are the same as these in the linear switched reluctance machine model.
Parameters Description
Ld (d-axis inductance) d-axis inductance of the PMSM machine, in H
Lq (q-axis inductance) q-axis inductance of the PMSM machine, in H
Vpk / krpm Peak line-to-line back emf constant of the machine, in V/krpm (mechanical speed)
Number of Poles Number of poles of the machine
Maximum Inverter Maximum inverter output current amplitude (peak), in A
Current
Parameters Description
Maximum Inverter Maximum inverter output current amplitude (peak), in A
Current
Parameters Description
Rs (stator resistance) Stator winding resistance of the machine, in Ohm
Ld (d-axis inductance) d-axis inductance of the PMSM machine, in H
Lq (q-axis inductance) q-axis inductance of the PMSM machine, in H
Vpk / krpm Peak line-to-line back emf constant of the machine, in V/krpm (mechanical speed)
Number of Poles Number of poles of the machine
Maximum Inverter Maximum inverter output current amplitude (peak), in A
Current
Maximum Inverter Maximum inverter output voltage amplitude (phase peak), in V
Voltage
Base Voltage Value System base voltage value, in V
Base Current Value System base current value, in A
Base Mechanical Speed System base mechanical speed, in rad/sec
The field weakening (IPM) control block accepts the following input signals, all in per-unit value.
Is: Inverter current amplitude reference
Vdc: DC bus voltage feedback
Wm: Motor mechanical speed, in rad/sec.
It has the following output signals, all in per-unit value:
Id: d-axis current reference
Iq: q-axis current reference
Parameters Description
Number of Poles Number of poles of the machine
Maximum Inverter Maximum inverter output current amplitude (peak), in A
Current
Maximum Inverter Maximum inverter output voltage amplitude (phase peak), in V
Voltage
Base Voltage Value System base voltage value, in V
Base Current Value System base current value, in A
Base Mechanical Speed System base mechanical speed, in rad/sec
The block has the following inputs and outputs, all in per unit value except Ld, Lq, and Lambda. If base voltage
value Vb, Ib, and Wmb are all set to 1, all input and output quantities are in true values.
Input Signals:
Is: Inverter current amplitude reference.
Ld: d-axis inductance at the rated operating conditions, in real value in H
Lq: q-axis inductance at the rated operating conditions, in real value in H
Lambda: Peak stator phase flux linkage at the rated operating conditions, in real value in Weber
Vdc: DC bus voltage feedback
Wm: Motor mechanical speed in
Output Signals:
Id: d-axis current reference.
Iq: q-axis current reference.
This block is for the control of nonlinear interior PMSM (IPM) machines only. It is not for surface mounted
PMSM (SPM) machines.
When the machine speed is higher than a certain value, the machine will not be able to generate the maximum
torque. Instead, it will be limited by the machine power rating. Given the current amplitude reference, dc bus
voltage, and the motor speed, the field weakening control block will calculate the d-axis and q-axis current
references Id and Iq to operate in the constant power region.
Attributes for Field Weakening (SPM):
Parameters Description
Ld (d-axis inductanc) The inductance on d-axis, in H
Vpk/krpm The ratio of peak voltage versus speed. in V/krpm
Number of Poles Number of poles of the machine
Maximum Inverter Maximum inverter output current amplitude (peak), in A
Current
Maximum Inverter Maximum inverter output voltage amplitude (phase peak), in V
Voltage
Parameters Description
Rs (stator) Stator winding resistance of the machine, in Ohm
Lls (stator leakage) Stator leakage inductance of the PMSM machine, in H
Rr (rotor) Rotor winding resistance of the machine, in Ohm
Llr (rotor leakage) Rotor leakage inductance of the machine, in H
Lm (magnetizing) Magnetizing inductance, in H
Number of Poles Number of poles of the machine
Rated Speed (rpm) Rated speed of the machine, in rpm
AC Source Frequency AC source frequency, in Hz
AC Voltage (line-line rms) AC source voltage, line-to-line rms, in V
Base Voltage Value System base voltage value, in V
Base Current Value System base current value, in A
The block has the following inputs and outputs, all in per unit value. If base values Vb, Ib, and Wmb are all set
at 1.0, all input and output quantities are in real values.
Input Signals:
Vdc: DC bus voltage
We: Stator electrical speed (in per unit, with the base value of Wmb)
Output Signal:
Id: d-axis current reference
This block is for the control of linear induction machines only.
When the machine speed is higher than a certain value, the machine will not be able to generate the maximum
torque. Instead, it will be limited by the machine power rating. Given the dc bus voltage and the stator speed,
the field weakening control block will calculate the d-axis current reference Id to operate in the constant power
region.
Attributes:
Parameter Description
Netlist XML File The file defines the interface between PSIM and JMAG. The file extension is .xml.
JMAG Input File The JMAG input data file. It has the .jcf extension. Note that the .xml and .jcf files
must be in the same directory.
JMAG Case Text Comments for the JMAG circuit
Terminal Names Terminal names of the block
The power circuit nodes and mechanical shaft nodes, as well as the rest of the interface between PSIM and the
JMAG data files, are defined in the Netlist XML File. This file is in XML format, and is generated by JMAG. To
specify this file, click on the browse button at the right of the edit field.
Parameter Description
dt_of_JMAG Time step used in JMAG.
Ir Rated current. It is used in JMAG
dl_threshold Current threshold. It is used in JMAG
dr_threshold Displacement threshold. It is used in JMAG
Moment of Inertia Moment of inertia of the machine, in kg*m2.
Shaft Time Constant Shaft time constant of the machine, in sec.
Back emf Flag Display flag for all FEM coil currents (1: display; 0: no display)
Rotor Angle Flag Display flag for all shaft angles
Speed Flag Display flag for all shaft speeds
Torque Flag Display flag for all developed torques
shaft1 Master Flag Master/slave flag of the machine (1: master; 0: slave)
Please consult JMAG manuals and documents on how to set up the co-simulation on the JMAG side.
Attributes:
Parameter Description
Link Table File The XML file that defines the interface between PSIM and JMAG. It has the .xml
extension.
JMAG Input File The JCF input data file for JMAG. It has the .jcf extension. Note that the .xml file and
the .jcf file must be in the same directory.
The setup process of calling JMAG in PSIM through the MagCoupler block is as follows.
In JMAG:
- In the JMAG circuit, connect a voltage function to the right of the FEM coil. Under Electrical
Potential in the property window, choose Constant Value, and set Constant Value [V] to 0.
- Connect a current probe to the left of the FEM coil.
- Connect another voltage function to the left of the current probe (the circuit will look like what is
shown above). In the property window, choose Cooperates with an external circuit simulator.
- Highlight the inductor structure window. Go to the menu Conditions >> Create Conditions. From
the Conditions List, highlight Coupled External Circuit Simulator, and click Modify.
- On the Coupled External Circuit Simulator dialog window, there are two lists. The list on the right,
called JMAG, contains all the functions that can be used to interface with PSIM. The list on the
left, called External Circuit Simulator, contains the functions that are selected to interface with
PSIM. In this case, there are two items in the JMAG list, one is the Voltage Function, and the other
is the Current Probe.
- Highlight the Voltage Function, and click on the <- button to move the item from the list on the right
to the list on the left. Repeat the same step to the Current Probe. Now both items should appear in
the list on the left.
- Highlight the Voltage Function, and change the terminal name to VL. Also, change the Current Probe
terminal name to iL. Close the dialog window.
- Go to the menu File >> Export and select JCF.... With the JCF file name defined as "inductor", the
JCF file "inductor.jcf" and the link table file "inductor_csl.xml" will be generated.
- Copy the JCF file "inductor.jcf" and the link table file "inductor_csl.xml" to the folder containing the
PSIM schematic file "inductor_jmag.sch". Rename the link table file to "inductor_jmag.xml". Note
that the XML file does not have to be renamed, and both the JCF and XML files do not have to be
moved to the folder of the schematic file. They are done here for the simplicity of file management
In PSIM:
- After the rest of the power circuit is created, go to Elements >> Power >> MagCoupler Module,
and select MagCoupler Block. Place the block on the schematic.
M- M+ M- M+
M- M+ M- M+
Parameter Description
Netlist XML File The XML file that defines the interface between PSIM and JMAG-RT. It has the .xml
extension.
JMAG-RT Input File The JMAG-RT data file. It has the .rtt extension. Note that the .xml file and the .rtt
file must be in the same directory.
JMAG Case Text Comments for the JMAG-RT circuit
Terminal Names Terminal names of the block
In the MagCoupler-RT block images, the power circuit nodes (such as A, B, C, A+, A-, B+, B-, C+, C-, D+, and
D- as shown above) are placed at the top of the block, arranged from the left to the right. The mechanical shaft
nodes are placed on the left and right of the block, with the first shaft node (such as M+ as shown above) on the
right, and the second shaft node (such as M-) on the left.
The electric nodes and rotor shaft nodes, as well as the rest of the interface between PSIM and the JMAG-RT
data files (.rtt file), are defined in the Netlist XML File. This file is in XML format, and is generated by the
JMAG-RT Manager. To specify this file, click on the browse button at the right of the edit field.
Note that for each type of MagCoupler-RT block, the "NetlistElement" must be the same as the predefined
name. The predefined names are:
- "PM Synchronous Machine" for 3-phase PMSM
- "Step Machine" for 2-phase step machine
- "Linear Synchronous Machine" for linear synchronous machine
- "Linear Solenoid" for linear solenoid
For example, the first several lines of a netlist XML file for a PMSM may look like the following:
<?xml version="1.0" encoding="UTF-8"?>
<CircuitElement file="imp.rtt" name="MagCoupler-RT">
<NetlistElement type="JMAGRT" element="PM Synchronous Machine" pole_number="4">
... ... ...
The JMAG-RT Input File is the JMAG-RT data file for the device modeled. The file has the .rtt extension and is
defined in the netlist XML file. Note that the .rtt file and the .xml file must be in the same directory.
The JMAG Case Text is a text identifying the specific JMAG-RT study case. It can be any text.
The Terminal Names are the names of the interface nodes.
The parameters of a 3-phase PMSM is shown below:
Parameter Description
shaft1_MomentofInertia Moment of inertia of the machine, in kg*m2
shaft1_ShaftTimeConstant Shaft time constant of the machine, in sec.
RU_resistance Resistance of Phase U, in Ohm
RV_resistance Resistance of Phase V, in Ohm
RW_resistance Resistance of Phase W, in Ohm
OffsetAngle Initial rotor angle, in mechanical deg.
coef_inductance Inductance coefficient (used in JMAG-RT)
coef_flux Flux coefficient (used in JMAG-RT)
coef_torque Torque coefficient (used in JMAG-RT)
coef_magnet Magnet coefficient (used in JMAG-RT)
Parameter Description
RA+ Resistance of Phase A+, in Ohm
RA- Resistance of Phase A-, in Ohm
RB+ Resistance of Phase B+, in Ohm
RB- Resistance of Phase B-, in Ohm
MomentofInertia Moment of inertia of the machine, in kg*m2
ShaftTimeConstant Shaft time constant of the machine, in sec.
OffsetAngle Initial rotor angle, in mechanical deg.
turns_coil2 Coil 2 turns (used in JMAG-RT)
coef_inertia1 Inertia 1 coefficient (used in JMAG-RT)
coef_inertia2 Inertia 2 coefficient (used in JMAG-RT)
Current Flag Display flag for 3-phase FEM coil currents (1: display; 0: no display)
Back emf Flag Display flag for 3-phase FEM coil back emf
Rotor Angle Flag Display flag for the rotor angle in mechanical deg.
Speed Flag Display flag for the shaft mechanical speed in rad./sec.
Torque Flag Display flag for the developed torque in N*m
shaft1 Master Flag Master/slave flag of the machine (1: master; 0: slave)
The parameters of a linear solenoid is shown below:
Parameter Description
RA Resistance of the solenoid, in Ohm
Mass Mass of the solenoid, in kg
ShaftTimeConstant Shaft time constant of the solenoid, in sec.
SpringConstant Spring constant of the solenoid (used in JMAG-RT)
DispLimitMax Maximum limit of the displacement of the solenoid, in m
DispLimitMin Minimum limit of the displacement of the solenoid, in m
OffsetDisp Initial displacement, in m
Parameter Description
RU Resistance of Phase U, in Ohm
RV Resistance of Phase V, in Ohm
RW Resistance of Phase W, in Ohm
Mass Mass of the machine, in kg
ShaftTimeConstant Shaft time constant of the solenoid, in sec.
OffsetDisp Initial displacement, in m
coef_inductance Inductance coefficient (used in JMAG-RT)
coef_flux Flux coefficient (used in JMAG-RT)
coef_force Torque coefficient (used in JMAG-RT)
coef_magnet Magnet coefficient (used in JMAG-RT)
coef_material Material coefficient (used in JMAG-RT)
turns_coil1 Coil 1 turns (used in JMAG-RT)
turns_coil2 Coil 2 turns (used in JMAG-RT)
coef_mass1 Mass 1 coefficient (used in JMAG-RT)
coef_mass2 Mass 2 coefficient (used in JMAG-RT)
Current Flag Display flag for 3-phase FEM coil currents (1: display; 0: no display)
Back emf Flag Display flag for 3-phase FEM coil back emf
Position Flag Display flag for the solenoid position in m
Velocity Flag Display flag for the solenoid velocity in m/sec.
Force Flag Display flag for the developed force in N
mass1 Master Flag Master/slave flag of the solenoid (1: master; 0: slave)
In all the parameter definitions above, except the flags (current flag, back emf flag, rotor angle/position flag,
speed/velocity flag, torque/force flag, and master flag), all other parameters are read from the .rtt file defined in
the XML file. These values can be changed, and for these parameters noted with "used in JMAG-RT", the
changed values will be sent back to JMAG-RT for calculation.
Several MagCoupler-RT examples are provided under the "examples\MagCoupler-RT" folder.
Attributes:
Parameters Description
Constant Torque Torque constant Tconst, in N*m
Moment of Inertia Moment of inertia of the load, in kg*m2
If the reference direction of a mechanical system enters the dotted terminal, the load is along the reference
direction, and the loading torque to the master machine is Tconst. Otherwise the loading torque will be -Tconst.
See Section 2.6.1 for more detailed explanation on the reference direction.
A constant-torque load is expressed as:
T L = T const
The torque does not depend on the speed direction.
Attributes:
Parameters Description
Maximum Torque Maximum torque Tmax of the load, in N*m
Base Speed Base speed nbase of the load, in rpm
Moment of Inertia Moment of inertia of the load, in kg*m2
The torque-speed curve of a constant-power load is shown below:
Torque
(N*m)
When the mechanical speed is less than the base speed nbase, the load torque is:
T L = T max
When the mechanical speed is above the base speed, the load torque is:
P
T L = ---------
m
Attributes:
Parameters Description
Constant Speed (rpm) Speed constant, in rpm
Moment of Inertia Moment of inertia of the load, in kg*m2
A constant-speed mechanical load defines the speed of a mechanical system, and the speed will remain
constant, as defined by the speed constant.
Attributes:
Parameters Description
Tc Constant torque term
k1 (coefficient) Coefficient for the linear term
Attributes:
Parameters Description
Speed Flag Flag for speed dependency
Flag = 0: The load is frictional and is always against the rotational direction;
Flag = 1: The load is independent of the rotational direction.
Moment of Inertia Moment of inertia of the load, in kg*m2
The value of the mechanical load is defined by the voltage value at the control node (1V corresponds to 1 N*m).
This node is a control circuit node.
Shaft 1
Shaft 2
Attribute:
Parameter Description
Gear Ratio The gear ratio a
Shaft 1 Master/Slave Flag Master/slave flag for Shaft 1
Shaft 2 Master/Slave Flag Master/slave flag for Shaft 2
The shaft with the bigger dot is Shaft 1.
If the numbers of teeth of the first gear and the second gear are n1 and n2, respectively, the gear ratio a is defined
as: a = n1 / n2. Let the radius, torque, and speed of these two gears be: r1, r2, T1, T2, 1, and 2, we have: T1 /
This block is used in situations where both mechanical systems have a device in the master mode, and they must
be connected together to form one system. Based on the connection convention in PSIM, a mechanical system
can have only one master device. In this case, the mechanical coupling block can be inserted in between, and the
mechanical system on each side of the coupling block can have its own device in the master mode.
Attribute:
Parameter Description
Master/Slave Flag Flag for the master/slave mode (1: master, 0: slave)
Similar to electric machines, the mechanical-electrical interface block can be used to define the reference
direction of a mechanical system through the master/slave flag. When the interface block is set to the master
mode, the reference direction is along the mechanical shaft, away from the mechanical node, and towards the
rest of the mechanical elements.
Lets assume that a drive system consists of a motor (with a developed torque of Tem and a moment of inertia of
J1) and a mechanical load (with a load torque of Tload and a moment of inertia of J2). The equation that
describes the mechanical system is:
d
( J 1 + J 2 ) ---------m = T em T load
dt
where m is the shaft mechanical speed. In PSIM, this equation is modelled by an equivalent circuit as shown
below.
m speed node
Tem J1 J2 Tload
In this circuit, the two current sources have the values of Tem and Tload, and the capacitors have the values of J1
and J2. The node-to-ground voltage (speed node voltage) represents the mechanical speed m. This is
analogous to C*dV/dt = i for a capacitor where C = J1+J2, V = m, and i = Tem-Tload.
Wm
Mechanical
speed
Attribute:
Parameter Description
Gain Gain of the sensor
The output of the speed sensor is in rpm.
The torque sensor measures the torque difference between the dotted side of the sensor and the undotted side.
To understand the physical meaning of the torque sensor measurement, we use the diagram below as an
illustration.
The figure on the left shows a torque sensor connected with a 10-N*m mechanical load, and the reference
direction of the mechanical system is from left to right. Based on the reference direction, if we use the right-
hand method, by pointing the thumb in the reference direction and rotating the right hand, the direction of the
fingers will show the direction of the positive speed and torque. The physical interpretation of the system is
shown on the right.
10 Wm Load
10
Torque sensor
In this case, the direction of the positive speed and torque is in the clockwise direction. The dotted side of the
sensor is on the left, and the load is in such a way that it tries to slow down the shaft (the load torque is in the
counter-clockwise direction).
The physical meaning of the torque sensor is that, if the dotted side of the sensor is fixed, the sensor will
measure the torque tension on the undotted side of the sensor, and a positive sensor output would mean that the
torque is opposite to the direction of the speed reference. Therefore, for the example above, the positive speed
reference is in the clockwise direction, and the load torque is in the counter-clockwise direction. This will give a
torque sensor reading of 10 N*m.
Similarly, if the undotted side of the sensor is fixed, the sensor will measure the torque tension on the dotted
side of the sensor, in the positive direction of the speed reference. For example, in the system below, the torque
sensor is flipped with the dotted side on the right. If the undotted side is fixed, the load torque is applied to the
dotted side of the sensor, in the opposite direction of the speed reference. The torque sensor output will be -10
N*m instead.
10 Wm Load
10
Torque sensor
To understand how the torque sensor is modeled in the equivalent circuit of the mechanical system, we use the
Load 1 Load 2
Sensor 1 Sensor 2
The system consists of one machine, 2 torque sensors, and 2 mechanical loads. The torques and moment of
inertia for the machine and the loads are as labelled in the diagram. The reference direction of this mechanical
system is from left to right. The equation for this system can be written as:
d
( J + J L 1 + J L 2 ) --------m- = T em T L 1 T L 2
dt
The equivalent electrical circuit of the equation is shown below:
Sensor 1 Sensor 2
m
The node voltage in the circuit represents the mechanical speed m. The current probe on the left represents the
reading of the Torque Sensor 1. Similarly, the current probe on the right represents the reading of the Torque
Sensor 2. Note that the second current probe is from right to left since Sensor 2 is opposite to the reference
direction of the mechanical system.
The equivalent circuit also illustrates how mechanical power is transferred. The multiplication of the current to
the voltage, which is the same as the torque times the mechanical speed, represents the mechanical power. If the
power is positive, it is transferred in the direction of the speed m.
Parameter Description
Initial Position (deg.) Initial shaft position, in deg.
No. of Bits of Resolution Number of Bits of resolution N
The encoder output resolution is determined by the number of bits N. The encoder has two outputs: one is the
number of counts (the range is from 0 to 2N-1), and the other is the position, in mechanical deg. (the range is
from 0 to 360o).
An example of a PMSM drive system using the absolute encoder is given in the sample file "Absolute Encoder
PMSM Drive.sch".
Attribute:
Parameter Description
Initial Position (deg.) Initial shaft position, in deg.
No. of Lines Number of lines that are in the code pattern of the code disk of the encoder
The two quadrature outputs are A and A (A is the inverse of A), and B and B. They are offset by 90o. In
addition, the encoder provides separate index signal output Z and Z that provide one count per revolution.
An example of an induction motor drive system using the incremental encoder is given in the sample file
"Incremental Encoder INDM Drive.sch".
4.9.6.3 Resolver
A resolver is essentially a rotary transformer with one rotor winding and two stator windings. These two stator
windings, referring to as the COS winding and SIN winding, are located 90o apart.
As the shaft rotates, the output voltages of the COS and SIN windings vary as the cosine and sine functions of
the shaft angle.
Image:
cos+
sin+
cos-
sin-
Parameter Description
Initial Position (deg.) Initial shaft position, in deg.
No. of Poles Number of poles of the resolver
The resolver has four outputs: cos+, cos- (the inverse of cos+), sin+, and sin- (the inverse of sin+). The peak
amplitude of all the outputs is 1.
An example of a PMSM drive system using the resolver is given in the sample file "Resolver PMSM
Drive.sch".
Attribute:
Parameter Description
Initial Position (deg.) Initial shaft position, in deg.
No. of Poles Number of poles of the sensor
The hall effect sensor provides three logic signal outputs A, B, and C, which are spaced 120 electrical deg.
apart.
The hall effect sensor is the same as the built-in hall effect sensor in the brushless dc machine.
Examples of BDCM motor drive systems using the hall effect sensor are given in the "examples" sub-folder of
the PSIM directory.
Psw
Pcond
In the images, beside the diode anode and cathode terminals, there are two extra nodes. The node with a dot is
for the diode conduction losses Pcond, and the node without a dot is for the diode switching losses Psw.
The style of the package (for example, TO-220, TO-247, etc.) can be specified in the Style field.
Absolute Maximum Ratings:
Vrrm,max (V): Peak reverse blocking voltage
IF,max (A): Maximum dc current
Tj,max (oC): Maximum junction temperature
Electrical Characteristics:
Vd vs. IF: Forward conduction voltage drop Vd vs. forward current IF
trr vs. IF: Reverse recovery time trr vs. current IF
Irr vs. IF: Peak reverse recovery current Irr vs. current IF
Qrr vs. IF: Reverse recovery charge Qrr vs. current IF
Err vs. IF: Reverse recovery energy losses Err vs. current IF
Thermal Characteristics:
Rth(j-c): Junction-to-case thermal resistance, in oC/W
Rth(c-s): Case-to-sink thermal resistance, in oC/W
Dimensions and Weight:
Length (mm): Length of the device, in mm
Height (mm): Height of the device, in mm
Weight (g): Weight of the device, in g
Parameters Description
Device The specific device selected from the device database
Frequency Frequency, in Hz, under which the losses are calculated
Pcond Calibration Factor The calibration factor Kcond of the conduction losses Pcond
Psw Calibration Factor The calibration factor Ksw of the switching losses Psw
Number of Parallel Number of identical diode devices in parallel
Devices
The parameter Frequency refers to the frequency under which the losses are calculated. For example, if the
device operates at the switching frequency of 10 kHz, and the parameter Frequency is also set to 10 kHz, the
losses will be the values for one switching period. However, if the parameter Frequency is set to 60 Hz, then the
losses will be the value for a period of 60 Hz.
The parameter Pcond Calibration Factor is the correction factor for the conduction losses. For the example, if
the calculated conduction losses before the correction is Pcond_cal, then
Pcond = Kcond * Pcond_cal
Similarly, the parameter Psw Calibration Factor is the correction factor for the switching losses. For the
example, if the calculated switching losses before the correction is Psw_cal, then
Psw = Ksw * Psw_cal
When several identical diodes are in parallel, one should have just one device in the schematic, and set the
correct number of devices in the parameter input. This is because when several identical devices are in parallel
in the schematic, the device currents may not be exactly equal due to small differences in the simulation. When
the number of parallel devices is greater than one, the total current through the group of the parallel devices will
be equally divided among the devices. The total losses are then obtained by multiplying the losses in each
device by the number of parallel devices. The currents flowing out of Nodes Pcond and Psw are the total losses of
all the parallel devices combined.
The voltage at the conduction losses node Pcond or the switching losses node Psw represents the calculated
junction temperature Tj of the diode, and this junction temperature is used when the database curves are used for
loss calculation. If the calculated Tj is between the junction temperatures of two curves, interpolation will be
used. If the calculated Tj is lower than the lowest Tj or higher than the highest Tj in the curves, the curve
corresponding to the lowest or highest Tj will be used. If there is only one curve, that curve is used regardless of
the calculated Tj.
Conduction Losses:
The diode conduction losses are calculated as:
Conduction Losses = Vd * IF
Q1
Pcond_Q
Psw_Q
Q4
Pcond_D
Psw_D
Q1 Q4
In the images, beside the IGBT and diode terminal nodes, there are four extra nodes from the top to the bottom
(or from the left to the right on the top for the 6-pack package). They are the nodes for transistor conductor
losses Pcond_Q (the node with a circle), for transistor switching losses Psw_Q, for diode conductor losses Pcond_D
(the node with a square), and for diode switching losses Psw_D, respectively.
The style of the package (for example, TO-247, TO-268, etc.) can be specified in the Style field.
Absolute Maximum Ratings:
Vce,max (V): Maximum collector-emitter voltage
Ic,max (A): Maximum collector current
Tj,max (oC): Maximum junction temperature
Electrical Characteristics - Transistor:
Vce(sat) vs. Ic: Collector-emitter saturation voltage Vce(sat) vs. collector current Ic
Eon vs. Ic: Turn-on energy losses Eon vs. collector current Ic
Eoff vs. Ic: Turn-off energy losses Eoff vs. collector current Ic
Electrical Characteristics - Diode (or Anti-Parallel Diode):
Vd vs. IF: Forward conduction voltage drop Vd vs. forward current IF
trr vs. IF: Reverse recovery time trr vs. current IF
Irr vs. IF: Peak reverse recovery current Irr vs. current IF
Qrr vs. IF: Reverse recovery charge Qrr vs. current IF
Err vs. IF: Reverse recovery energy losses Err vs. current IF
Electrical Characteristics - Free-Wheeling Diode (for IGBT-Diode device only):
Vd vs. IF: Forward conduction voltage drop vs. forward current IF
trr vs. IF: Reverse recovery time trr vs. current IF
Irr vs. IF: Peak reverse recovery current Irr vs. current IF
Qrr vs. IF: Reverse recovery charge Qrr vs. current IF
Err vs. IF: Reverse recovery charge Err vs. current IF
Thermal Characteristics:
Rth(j-c) (transistor): Transistor junction-to-case thermal resistance, in oC/W
Parameters Description
Device The specific device selected from the device database
Frequency Frequency, in Hz, under which the losses are calculated
Pcond_Q Calibration The calibration factor Kcond_Q of the transistor conduction losses Pcond_Q
Factor
Psw_Q Calibration The calibration factor Ksw_Q of the transistor switching losses Psw_Q
Factor
Pcond_D Calibration The calibration factor Kcond_D of the diode conduction losses Pcond_D
Factor
Psw_D Calibration The calibration factor Ksw_D of the diode switching losses Psw_D
Factor
Number of Parallel Number of identical IGBT devices in parallel
Devices
The parameter Frequency refers to the frequency under which the losses are calculated. For example, if the
device operates at the switching frequency of 10 kHz, and the parameter Frequency is also set to 10 kHz, the
losses will be the values for one switching period. However, if the parameter Frequency is set to 60 Hz, then the
losses will be the value for a period of 60 Hz.
The parameter Pcond_Q Calibration Factor is the correction factor for the transistor conduction losses. For the
example, if the calculated conduction losses before the correction is Pcond_Q_cal, then
Pcond_Q = Kcond_Q * Pcond_Q_cal
Similarly, the parameter Psw_Q Calibration Factor is the correction factor for the transistor switching losses.
For the example, if the calculated switching losses before the correction is Psw_Q_cal, then
Psw_Q = Ksw_Q * Psw_Q_cal
Parameters Pcond_D Calibration Factor and Psw_D Calibration Factor work in the same way, except that they
Pcond_Q
Psw_Q
Pcond_D
Psw_D
The calculated junction temperatures Tj_Q and Tj_D are used when the database curves are used for loss
calculation. If the calculated junction temperature is between the junction temperatures of two curves,
interpolation will be used. If the calculated junction temperature is lower than the lowest junction temperature
or higher than the highest junction temperature in the curves, the curve corresponding to the lowest or highest
junction temperature will be used. If there is only one curve, that curve is used regardless of the calculated
junction temperature.
Conduction Losses:
The transistor conduction losses is calculated as:
Transistor Conduction Losses = Vce(sat) * Ic
where Vce(sat) is the transistor collector-emitter saturation voltage, and Ic is the collector current. When the
transistor is conducting periodically with an on duty cycle of D, the conduction losses are calculated as:
Transistor Conduction Losses = Vce(sat) * Ic * D
Switching Losses:
The transistor turn-on losses is calculated as:
Transistor Turn-on Losses = Eon * f * Vcc / Vcc_datasheet
where Eon is the transistor turn-on energy losses, f is the frequency as defined in the input parameter Frequency,
Vcc is the actual dc bus voltage, and Vcc_datasheet is the dc bus voltage in the Eon and Eoff characteristics of the
datasheet, defined as "DC bus voltage (V)" in the test conditions.
The transistor turn-off losses is calculated as:
Transistor Turn-off Losses = Eoff * f * Vcc / Vcc_datasheet
where Eoff is the transistor turn-off energy losses.
The loss calculation for the anti-parallel diode or free-wheeling diode is the same as described in the section for
the diode device.
The losses Pcond_Q, Psw_Q, Pcond_D, and Psw_D, in watts, are represented in the form of currents which flow out
of these nodes. Therefore, to measure and display the losses, an ammeter should be connected between the
nodes and the ground. When they are not used, these nodes cannot be floating and must be connected to ground.
(n-channel) (p-channel)
Pcond_Q
Q1
Psw_Q
Pcond_D Q4
Psw_D
Q1 Q4
The base values are normally obtained at the test conditions of 25oC. Using the normalized value of RDS(on)
based on the base value, we have:
R ( DS ) on _normalized 1
K T = -----------------------------------------------
-
T j 25
The on-resistance RDS(on) is calculated at every time step and is used in the simulation.
Also, the forward conduction voltage drop Vd of the diode changes depending on the current. The new
conduction voltage value is used in the simulation.
Parameters Description
Device The specific device selected from the device database
Frequency Frequency, in Hz, under which the losses are calculated
VGG+ (upper level) Upper level of the gate source voltage, in V
VGG- (lower level) Lower level of the gate source voltage, in V
Rg_on (turn-on) Gate resistance during turn-on
Rg_off (turn-off) Gate resistance during turn-off. In most cases, the turn-on gate resistance Rg_on
and the turn-off gate resistance Rg_off are identical.
RDS(on) Calibration Factor The calibration factor of the on-state resistance RDS(on)
gfs Calibration Factor The calibration factor of the forward transconductance gfs
Pcond_Q Calibration Factor The calibration factor Kcond_Q of the transistor conduction losses Pcond_Q
Psw_Q Calibration Factor The calibration factor Ksw_Q of the transistor switching losses Psw_Q
Pcond_D Calibration Factor The calibration factor Kcond_D of the diode conduction losses Pcond_D
Psw_D Calibration Factor The calibration factor Ksw_D of the diode switching losses Psw_D
Number of Parallel Devices Number of identical devices in parallel
The parameter Frequency refers to the frequency under which the losses are calculated. For example, if the
device operates at the switching frequency of 10 kHz, and the parameter Frequency is also set to 10 kHz, the
losses will be the values for one switching period. However, if the parameter Frequency is set to 60 Hz, then the
losses will be the value for a period of 60 Hz.
The parameter Pcond_Q Calibration Factor is the correction factor for the transistor conduction losses. For the
example, if the calculated conduction losses before the correction is Pcond_Q_cal, then
Pcond_Q = Kcond_Q * Pcond_Q_cal
Electrical Characteristics:
L (uH): Inductance, in uH
Maximum Ratings:
Imax,rms (A): Maximum rms current rating, in A.
Core Information:
Core Type and Size: Specify core shape and size
Core Material: Specify core material and magnetic characteristics
Winding Information:
Winding Type and Size: Specify conductor shape, size, and insulation type
Winding Distribution: Specify winding number of turns, layers, distance between layers, distance
between wires, and distance between core and winding
Distribution of Parallel Wires: Specify parallel wire distribution in radial and height directions
Gap Information:
Gap Data: Specify the number of gaps and the length of each gap.
A new inductor can be added to the database by selecting Device >> New Inductor. Databases of core
materials, cores, and winding can be added and edited by selecting Device >> Edit Inductor Core Database,
Edit Winding Database, and Edit Core Material Database.
Parameters Description
Device The specific inductor selected from the device database.
Frequency Fundamental frequency of the inductor, in Hz.
Temperature Flag The flag that defines how core temperature is defined. It can be either
Calculated or Fixed. If it is Calculated, the core temperature will be
calculated based on circuit operating conditions and this temperature will
be used for loss calculation. If it is Fixed, the core temperature will be
specified by users. This will give users a clear idea the losses at a specific
temperature.
Ambient Temperature The ambient temperature of the inductor, in oC.
Convection Type of cooling. It can be Forced for forced convection or Natural for
natural convection.
Air Speed The air speed, in m/sec., if forced convection is selected.
Estimated Core Temp. Estimated core temperature, in oC, if the Temperature Flag is set to Fixed.
Estimated Winding Temp. Estimated winding temperature, in oC, if the Temperature Flag is set to
Fixed.
Loss Calculation Flag If the flag is 0, the loss calculation will be performed from the beginning of
the simulation. If the flag is 1, the loss calculation will be performed only in
the last fundamental cycle of the simulation. Since loss calculation will
slow down simulation, calculating the losses only in the last cycle will
speed up the simulation.
Current Flag Display flag of the inductor current.
The losses Pcore and Pwinding, in watts, are represented in the form of currents which flow out of the nodes. To
measure and display the losses, an ammeter should be connected between the nodes and the ground. When they
are not used, these nodes cannot be floating and must be connected to ground.
In the parameters, Frequency refers to the main frequency under which the inductor is excited. For example, if
an inductor conducts a 60-Hz current in an ac circuit which also has 20kHz switching harmonic components,
the frequency should be set to 60 Hz. If an inductor conducts dc current in a dc circuit which has 20kHz
switching harmonic components, the frequency should be set to 20kHz.
For more information on how to use the Thermal Module inductor, please refer to the tutorial "Tutorial -
Inductor loss calculation in Thermal Module.pdf".
Device
information
Device
list
On the left are the device database files that are loaded into the database editor, and the list of the devices. The
devices can be displayed based on either Device Type or Manufacturer. Also, the device list can be sorted by
Part Number, Voltage rating, or Current rating, by clicking on the title bars of the list.
To create a new device file, choose File >> New Device File.
To load a device files into the editor, choose File >> Open Device File.
To unload a device file from the editor, choose File >> Close Device File.
On the right is the information of each device. In general, the following information is required to define a
switching device for thermal simulation:
- Manufacturer and Part Number
- Package type and style
- Absolute maximum ratings
- Electrical characteristics
- Thermal characteristics
- Dimension and weight
The following information are required to define a database inductor for thermal simulation.
- Manufacturer and Part Number
- Package type and style
- Maximum ratings
- Electrical characteristics
- Core, winding, and gap information
Five types of elements can be added to a database device file: inductor and four types of switching devices:
diode, IGBT, dual IGBT-diode modules, and MOSFET.
To create a new device, go to the Device menu, and choose either New Diode, New IGBT, New IGBT-Diode,
New MOSFET, or New Inductor.
To make a copy of an existing device in the same database file, highlight the device in the list, and choose
Device >> Save Device As.
Graph
wizard
icons
Help area
X and Y
axis settings
Data area
X/Y axis
multiplying
factor
Graph area
Click on the forward wizard icon to paste the screen image into the graph window in the database editor.
Position the graph image properly in the graph window by dragging the left mouse so that the complete
graph is displayed within the window.
If the graph image is either too large or too small, go back to the previous step by clicking on the backward
wizard icon. Then resize the image of the graph in the Adobe Acrobat, and copy the screen image to the
clipboard again.
Opposite
end of the
origin
Origin
of the
graph
In this step, the border of the graph area is defined by first left clicking at the origin of the graph (usually the
lower left corner), then left clicking again at the opposite corner of the origin (usually the upper right cor-
ner). Note that the graph origin does not have to be the lower left corner, and it can be any one of the four
corners.
To locate the origin of the corner more accurately, right mouse click to zoom in, and press the Esc (escape)
key to exit the zoom. After this, a blue rectangle will appear around the border of the graph, and the dialog
window will appear as follows.
Data
Points
To see the x and y axis values of a particular data point on the graph, place the cursor inside the graph area.
The cursor image will change to a cross image, and the x and y coordinates of the cursor will be displayed
at the upper right corner of the dialog window. Place the cursor on top of the curve to read the x and y axis
readings.
If there are multiple curves for different junction temperatures, repeat the same process and enter the junc-
tion temperature value for each curve accordingly. Use the drop-down menu at the upper right corner of the
dialog window to show one particular curve, or to show all the curves at the same time.
With the same process, define the reverse recovery characteristics trr vs. IF, Irr vs. IF, and Qrr vs. IF.
Enter the Thermal Characteristics as:
Rth(j-c): 0.6
Rth(c-s): 0.4
Enter the Dimension and Weight as:
Length (mm):53
Width (mm):36
Height (mm):29
Choose Device >> Save Device to save the device information. This completes the process of adding the
diode into the database.
In the image, The nodes marked with the "+" and "-" signs are the positive and negative terminals. The node on
the top is theoretical maximum power (in W) given the operating conditions.
While the positive and negative terminal nodes are power circuit nodes, the other nodes are all control circuit
nodes.
Attributes:
Parameter Description
Open Circuit Voltage Voc Voltage measured when solar cell terminals are open circuit, in V.
Short Circuit Current Isc Current measured when the solar cell terminals are short circuited, in A.
Maximum Power Voltage Vm Solar cell terminal voltage when the output power is at the maximum, in V
Maximum Power Current Im Solar cell terminal current when the output power is at the maximum, in A
Using the four input parameters, the functional model will create the i-v curve of the solar cell. The figure
below shows the i-v and p-v curves of a typical solar cell.
Voltage
The curves show that the solar cell output power reaches the maximum at a specific voltage level. Many control
In the image, The nodes with the "+" and "-" signs are the positive and negative terminals. The node with the
letter "S" refers to the light intensity input (in W/m2), and The node with the letter "T" refers to the ambient
temperature input (in oC). The node on the top is theoretical power (in W) given the operating conditions.
While the positive and negative terminal nodes are power circuit nodes, the other nodes are all control circuit
nodes.
Attributes
Parameter Description
Number of Cells Ns Number of cells Ns of the solar module. A solar module consists of Ns solar
cells in series.
Standard Light Intensity S0 Light intensity S0 under the standard test conditions, in W/m2. The value is
normally 1000 W/m2 in manufacturer datasheet.
Ref. Temperature Tref Temperature Tref under the standard test conditions, in oC
Series Resistance Rs Series resistance Rs of each solar cell, in Ohm
Shunt Resistance Rsh Shunt resistance Rsh of each solar cell, in Ohm
Short Circuit Current Isc0 Short circuit current Isc0 of each solar cell at the reference temperature Tref,
in A
Saturation Current Is0 Diode saturation current Is0 of each solar cell at the reference temperature
Tref, in A
Band Energy Eg Band energy of each solar cell, in eV. It is around 1.12 for crystalline silicon,
and around 1.75 for amorphous silicon.
Ideality Factor A Ideality factor A of each solar cell, also called emission coefficient. It is
around 2 for crystalline silicon, and is less than 2 for amorphous silicon.
+
id ir +
iph vd v/Ns
Rsh
--
--
Some of the parameters of the physical model can be obtained from manufacturer datasheet, and the rest of the
parameters can be obtained by trial-and-error. A utility tool Solar Module (physical model) under the Utilities
T T
In the image, the nodes marked with the "+" and "-" signs are the positive and negative terminals. The node with
the letter "S" refers to the light intensity input (in W/m2), and The node with the letter "T" refers to the ambient
temperature input (in oC).
While the positive and negative terminal nodes are power circuit nodes, the other nodes are all control circuit
nodes.
Attributes:
Parameter Description
Maximum Power Solar cell maximum output power, in W.
Maximum Power Voltage Solar cell terminal voltage when the output power is at the maximum, in V
Test Condition Temperature Test condition temperature, in degree C.
The current-voltage characteristic of cSi and thin-film models equations are:
( Vpv + IpvRs ) ( mVr )
I pv = I ph I 0 ( e 1 ) ( V pv + I pv R s ) R p
Where,
Vpv
---------
3 Vr
I 0 = C 0 T mod e
and
kT mod
V T = ------------
-
e0
V mpp 2
G = 200 ( W m ) 0.95 0.98 +/- 1%
---------------------------------------
-
V mpp 2
G = 1000 ( W m )
V mppSTC
-----------------
V ocSTC 0.8 0.72 < 1%
I mppSTC
----------------- 0.9 0.8 < 1%
I scSTC
Where the temperature of the PV generator should follow the ambient conditions as follows:
k
T PV = T amb + T 0 + ------------------- G
1+s
where
TPV Computed PV generator temperature;
Voltage ratio from VMPP at an irradiance of 200W/m2 to VMPP at an irradiance of 1000W/m2 is given as:
V MPP 2
G = 200 ( W m )
V L 2 H = ---------------------------------------
-
V MPP 2
G = 1000 ( W m )
The parameters of the PV generator model must be set as listed in the table below:
cSi Thin-Film Tolerance
FFV 0.8 0.72 < 1%
FFI 0.9 0.8 < 1%
CV 8.593E-02 8.419E-02 -
0.04 0.02
-0.4 -0.2
Wind Turbine
Wind Speed
Pitch Angle
Parameter Description
Nominal Output Power The maximum output power of the wind turbine, in W, at 0o pitch angle. This
power is considered as the maximum power operating point of the turbine, and
it should not exceed the rated power of the generator.
Base Wind Speed The base wind speed that would produce the nominal output power, in m/s
Base Rotational Speed The base rotational speed of the turbine that would produce the nominal output
power, in rpm
Initial Rotational Speed The initial rotational speed of the turbine, in rpm
Moment of Inertia Moment of inertia of the wind turbine blade, in kg*m2
Torque Flag Flag to display the internal torque of the wind turbine (0: no display; 1:
display)
Master/Slave Flag Master/slave flag for the connected mechanical system (0: slave; 1: master)
In the image, the node with the letter "w" is for the wind speed input in m/s, and the node with the letter "p" is
for the blade pitch angle input in deg. Both nodes are control circuit nodes.
The power generated by a wind turbine can be expressed as:
1 3
P = --- A v wind C p
2
2), vwind is the wind speed (in m/sec.), is the air density (it is
where A is the area of the rotor blade (in m
3),
approximately 1.225 kg/m and Cp is the power coefficient. The power coefficient Cp is a function of the tip
speed ratio and the blade pitch angle . It can be expressed as [Heier1]:
x c6
Cp = c1 ( c2 c3 c4 c5 ) e + c7
where m is the rotor rotational speed (in rad/sec.) and Rblade is the radius of the rotor blade (in m).
The relationship between the power coefficient Cp and the tip speed ratio and the blade pitch angle can be
plotted in the figure below.
1. S. Heier and R. Waddington, Grid Integration of Wind Energy Conversion Systems, Wiley, 2006
Cp 0.20 =10
0.10
0.0
The figure shows that the power coefficient Cp reaches the maximum of 0.49 when the tip speed ratio is 8.18.
We choose the values of Cp and at the maximum as the nominal values, i.e., Cp_nom = 0.49, nom = 8.18. One
way to control the wind turbine, therefore, is to maintain the tip speed ratio at or close to the nominal value so
that the power that the wind turbine generates is at the maximum.
Several elements related to renewable energy applications are included here.
SOC
Attributes:
Parameter Description
No. of Cells in Series Number of cells Ns in series of a battery pack
No. of Cells in Parallel Number of cells Np in parallel of a battery pack
Voltage Derating Factor Voltage derating factor Ks, from 0 (100% derating) to 1 (no derating)
Capacity Derating Factor Capacity derating factor Kp, from 0 (100% derating) to 1 (no derating)
Rated Voltage Rated voltage Erated of the battery cell, in V
Discharge Cut-off Voltage Battery voltage Ecut corresponding to the maximum capacity, in V
Rated Capacity Rated capacity Qrated of the battery cell, in Ah
Internal Resistance Internal resistance Rbattery of the battery cell, in Ohm
Full Voltage Full (or maximum) voltage Efull of the battery cell, in V
Exponential Point Voltage Battery voltage Etop at the end of the exponential region in the discharge curve,
in V
Attributes:
Parameter Description
Number of Cells in Series Number of cells in series
Number of Cells in Parallel Number of cells in parallel
Capacitance per Cell Nominal capacitance per cell, in F.
Coefficient Kv Voltage coefficient
Resistance R1 Resistance coefficient R1, in Ohm
Capacitance C1 Capacitance coefficient C1, in F
This chapter provides descriptions for the components in PSIMs element library.
Two types of transfer function blocks are provided: one with zero initial values (the element is called s-domain
Transfer Function in the PSIM library) and the other with initial values as input parameters (called s-domain
Transfer Function (initial value) in the PSIM library).
Image:
Attributes:
Parameters Description
Order n Order n of the transfer function
Gain Gain k of the transfer function
Coeff. Bn...Bo Coefficients of the numerator (from Bn to Bo)
Coeff. An...Ao Coefficients of the denominator (from An to Ao)
Initial Values xn..x1 Initial values of the state variables xn to x1 (for the element s-domain Transfer
Function (initial value) only)
Let Y(s) = G(s)*U(s) where Y(s) is the output and U(s) is the input, we can convert the s-domain expression into
the differential equation form as follows:
x1 0 0 0 ... 0 A0 An x1 B0 A0 Bn An
x2 1 0 0 ... 0 A1 An x2 B1 A1 Bn An
d
------ x k
= 0 1 0 ... 0 A2 An x3 + ----- B2 A2 Bn An u
dt 3 An
... ... ... ... ... ... ... ... ...
xn 0 0 0 ... 1 A n 1 A n xn Bn 1 An 1 Bn An
Order n 2
Gain 1.5
Coeff. Bn...Bo 0. 0. 400.e3
Coeff. An...Ao 1. 1200. 400.e3
Attribute:
Parameter Description
Gain Gain k of the transfer function
5.1.2 Integrator
There are three types of integrators: regular integrator, external resettable integrator, and internal resettable
integrator.
Images:
Regular Integrator External Resettable Integrator Internal Resettable Integrator
Attributes:
For Regular Integrator:
Parameters Description
Time Constant Time constant T of the integrator, in sec.
Initial Output Value Initial value of the output
For External Resettable Integrator:
Parameters Description
Time Constant Time constant T of the integrator, in sec.
Initial Output Value Initial value of the output
Reset Flag Reset flag (0: edge reset; 1: level reset)
Parameters Description
Time Constant Time constant T of the integrator, in sec.
Initial Output Value Initial value of the output
Lower Output Limit Lower limit of the output
Upper Output Limit Upper limit of the output
The transfer function of an integrator is:
1
G ( s ) = ------
sT
The Bode plot of the amplitude |G| and the phase angle of the integrator is shown as below:
|G|
1/T
0 (rad/sec)
-20dB/dec
-90o
The output of the external resettable integrator can be reset by an external control signal (at the bottom of the
block). For the edge reset (reset flag = 0), the integrator output is reset to zero at the rising edge of the control
signal. For the level reset (reset flag = 1), the integrator output is reset to zero as long as the control signal is
high (1).
The output of the internal resettable integrator is reset to 0 when the output reaches either the lower limit or the
upper limit. It works in the same way as the external resettable integrator with the edge reset, except that in this
case users do not need to set up the external reset circuit.
When a limiter is connected to the integrator output, anti-windup is implemented automatically.
Example:
The following circuit illustrates the use of the resettable integrator. The input of the integrator is a dc quantity.
The control input of the integrator is a pulse waveform which resets the integrator output at the end of each
cycle. The reset flag is set to 0.
Vd
vo
vctrl
Attribute:
Parameter Description
Time Constant Time constant T of the differentiator, in sec.
Since sudden changes of the input will generate spikes at the output, it is recommended that a low-pass filter be
placed at the input of the differentiator.
Attributes:
Parameters Description
Gain Gain k of the PI controller
Time Constant Time constant T of the PI controller
The transfer function of a PI controller is defined as:
1 + sT
G ( s ) = k ---------------
sT
The Bode plot of the amplitude |G| and the phase angle of the PI controller is shown as below:
|G| -20dB/dec
0 (rad/sec)
1/T
-90o
Attributes:
Parameters Description
Gain Gain k of the controller
Pole Frequency Frequency of the pole fc, in Hz
where c = 2 f c .
The Bode plot of the amplitude |G| and the phase angle of the single-pole controller is shown as below:
|G|
-20dB/dec
0 c (rad/sec)
-90o
Attributes:
Parameters Description
Gain Gain k of the PI controller
Time Constant Time constant T of the PI portion of the controller
Pole Frequency Frequency of the pole fc, in Hz
0 (rad/sec)
1/T 1/Tc -20dB/dec
-90o
Attributes:
Parameters Description
Gain Gain k of the PI controller
Frequency fz1 Frequency fz1 of the first zero, in Hz
Frequency fz2 Frequency fz2 of the second zero, in Hz
Frequency fp1 Frequency fp1 of the first pole, in Hz
Frequency fp2 Frequency fp2 of the second pole, in Hz
1 1 1 1
where T z 1 = -------------- , T z 2 = -------------- , T p 1 = -------------- , and T p 2 = -------------- .
2 f z 1 2 f z 2 2 f p 1 2 f p 2
The Bode plot of the amplitude |G| and the phase angle of the Type-3 controller is shown as below:
0 (rad/sec)
1/Tz1 1/Tz2 1/Tp1 1/Tp2
90o
0
-90o
2nd-order low-pass 2nd-order high-pass 2nd-order band-pass 2nd-order band-stop 1st-order low-pass
Parameters Description
Gain Gain k
Damping Ratio Damping ratio
Cut-off Frequency Cut-off frequency fc ( f c = -----c- ) for low-pass and high-pass filters, in Hz
2
Center Frequency Center frequency fo ( f o = -----o- ) for band-pass and band-stop filter, in Hz
2
Passing Band; B
Frequency width fb ( f b = ------ ) of the passing/stopping band for band-pass/band-
Stopping Band 2
stop filters, in Hz
Attributes for first-order low-pass filter:
Parameters Description
Gain Gain k
Cut-off Frequency Cut-off frequency fc ( f c = -----c- ) of the low-pass filter, in Hz
2
The transfer function of these filters are listed below.
For second-order low-pass filter:
2
c
G ( s ) = k ----------------------------------------
2
-
2
s + 2 c s + c
For second-order high-pass filter:
2
s
G ( s ) = k ----------------------------------------
2
-
2
s + 2 c s + c
For second-order band-pass filter:
Bs
G ( s ) = k ----------------------------------
2
-
2
s + B s + o
For second-order band-stop filter:
2 2
s + o
G ( s ) = k ----------------------------------
2
-
2
s + B s + o
For first-order low-pass filter:
c
G ( s ) = k ---------------
s + c
5.2.1 Summer
The input of a one-input summer or two-input summer can be either a scalar or a vector. The input of a three-
input summer can only be a scalar.
Images:
1-input 2-input 2-input 3-input
Input 1
Input 1 Input 1
Input 2
Attribute:
Parameter Description
Gain_i Gain ki for the ith input
For the three-input summer, the input with a dot is the first input.
If the inputs are scalar, the output of a summer with n inputs is defined as:
V o = k 1 V 1 + k 2 V 2 + ... + k n V n
If the input is a vector, the output of a two-input summer will also be a vector, which is defined as:
V1 = [a1 a2 ... an]
V2 = [b1 b2 ... bn]
Vo = V1 + V2 = [a1+b1 a2+b2 ... an+bn]
The output of a one-input summer, however, will still be a scalar which is equal to the summation of the input
vector elements, that is, Vo = a1 + a2 + ... an.
Multiplier Divider
Numerator
Denominator
For the divider, the dotted node is for the numerator input.
The input of a multiplier can be either a vector or a scalar. If the two inputs are vectors, their dimensions must
be equal. Let the two inputs be:
V1 = [a1 a2 ... an]
V2 = [b1 b2 ... bn]
The output, which is a scalar, will be:
Vo = V1 * V2T = a1*b1 + a2*b2 + an*bn
Parameters Description
Coefficient k1 Coefficient k1
Coefficient k2 Coefficient k2
For example, if k1 = 1, k2 = 2.718281828, and Vin = 2.5, then Vo = e2.5 where e is the base of the natural
logarithm.
The output of a power function block is defined as:
k2
V o = k 1 V in
The function block LOG gives the natural logarithm (base e) of the input, and the block LOG10 gives the
common logarithm (base 10) of the input.
where T = 1/fb. The output is only updated at the beginning of each period.
Image:
Attribute:
Parameter Description
Base frequency Base frequency fb, in Hz
In the arctangent and arctangent 2 blocks, the input node marked with x is the real part, and the other input node
marked with y is the imaginary part.
The output of the arctangent and arctangent 2 blocks is the inverse tangent of the ratio between the imaginary
1 y
part y and the real part x, i.e. = tg ---. The output of the arctangent block is in degree, and the range is from
x
o. On the other hand, the output of the arctangent 2 block is in radian, and the range is from - to +.
0 to 360
The arctangent 2 block behaves in the same way as the function atan2(y,x) in the C language.
v in ( n ) v in n + --- e
N
v o = ---
N 2
n=0
Image:
Amplitude
Phase Angle
Parameters Description
No. of Sampling Points No. of sampling points N
Fundamental Frequency Fundamental frequency fb, in Hz.
The dotted node of the block refers to the output of the amplitude. Note that the phase angle output has been
internally adjusted such that a sine function Vm*sin(t) will give a phase angle output of 0.
Note that the FFT block only works if the input signal contains dc, the fundamental component, and harmonics
of the fundamental components (with the frequencies at multiple integer times of the fundamental frequency). If
the input signal contains another component that is not multiple integer of the fundamental frequency, the FFT
block will not give the correct result.
Also, the FFT block only measures the fundamental component, not any other harmonics. For example, if the
input signal contains 60 Hz (fundamental) and 180 Hz, if the parameter "Fundamental Frequency" is set to 180,
the FFT block output will be incorrect.
Example:
In the circuit below, the voltage vin contains a fundamental component v1 (100 V at 60 Hz), a 5th harmonic
voltage v5 (25 V at 300 Hz), and a 7th harmonic v7 (25 V at 420 Hz). After one cycle, the FFT block output
reaches the steady state with the amplitude of 100 V and the phase angle of 0o.
v1
vin vamp v1
vin
v5 Angle
v7 vamp
Angle
Attributes:
Parameters Description
Number of Inputs Number of inputs
Block Function Type Function type of the block. It can be either Maximum or Minimum.
The number of inputs can be changed through the drop-down menu.
5.3.1 Comparator
The output of a comparator is high when the positive input is higher than the negative input. When the positive
input is lower, the output is zero. If the two input are equal, the output is undefined and it will keep the previous
value.
Image:
Note that the comparator image is similar to that of the op. amp. For the comparator, the noninverting input is at
the upper left and the inverting input is at the lower left. For the op. amp., however, it is the opposite.
5.3.2 Limiters
Four types of limiters are provided to limit the amplitude of a signal: lower/upper limiter, lower limiter, upper
limiter, and range limiter.
Images:
Lower/upper Limiter Lower Limiter Upper Limiter Range Limiter
Attributes:
Parameters Description
Lower Limit Lower limit Vlower of the limiter
Upper Limit Upper limit Vupper of the limiter
For a lower/upper limiter, or lower limiter, or upper limiter, the output of the limiter is clamped to the upper or
lower limit whenever the input exceeds the limiter range. If the input is within the limit, the output is equal to
the input.
A range limiter limits the output vo between the lower limit and the upper limit, that is:
V lower v o < V upper
Let the range of the upper limit and the lower limit be Vrange. When the output is equal to or greater than the
upper limit, the output is subtracted by the range Vrange until it is within the range. When the output is below the
lower limit, it is added by range Vrange until it is within the range. When the input is within the limit, the output
is equal to the input.
Parameter Description
dv/dt Limit Limit of the rate of change (dv/dt) of the input
Parameters Description
Rising Angle theta Rising angle , in deg.
Peak Value Peak value Vpk of the waveform
Parameter Description
Pulse Width (deg.) Pulse width in half cycle, in deg.
The waveforms of these two blocks are shown below. Note that the input vin is in deg., and can be in the range
of -360o to 360o. Both waveforms are half-wave and quarter-wave symmetrical.
Vpk 1
o vin vin
180
0 0 o o
o
360 180 360
-Vpk -1
Input
Control
The difference between this block and the zero-order hold block (ZOH) in Digital Control Module is that this
block is treated as a continuous element and sampling moments can be controlled externally; whereas the zero-
order hold block is a discrete element and the sampling moments are fixed and of equal distance.
vctrl
Attributes:
Parameters Description
No. of Digits No. of digits N after the decimal point
Truncation Flag Truncation flag (1: truncation; 0: round-off)
Let the input of the round-off block be Vin. The input is first scaled based on the following expression:
N
V in, new = V in 10
If the truncation flag is 1, the output will be equal to Vin,new truncated, and then divided by 10N. Otherwise, the
output will be equal to Vin,new rounded off to the nearest integer, and then divided by 10N.
Examples:
If Vin = 34.5678, N = 0, truncation flag = 0, then we have the output Vout = 35. Similarly, if Vin = 34.5678, N = 0,
truncation flag = 1, then Vout = 34.
If Vin = 34.5678, N = 1, truncation flag = 1, then Vout = 34.5. If Vin = 34.5678, N = -1, truncation flag = 1, then
Vout = 30.
Parameter Description
Delay Time Delay time, in sec. (for the Time Delay block only)
A time delay block delays the input signal by a specified amount of time interval. It, for example. can be used to
model the propagation delay of a logic element.
A unit time delay block delays the input signal by just one simulation time step. If the simulation time step is
changed, the delay time will be changed accordingly.
Note that the difference between the unit time delay block and the unit delay block in Digital Control Module is
that this block is a continuous element and the delay time is one simulation time step; whereas the unit delay
block is a discrete element and the delay time is equal to the sampling period.
For a discrete system, the unit delay block should be used.
Example:
In this circuit, the first time delay block has a delay time of 1 ms, and the second block has a delay time of 4 ms.
This example illustrates that the input of the time delay block can be either an analog or a digital signal.
1 ms
vin1
vo1
vin2
vo2
4 ms
vin2
vo2
5.3.8 Multiplexer
The output of a multiplexer is equal to a selected input depending on the control signal. Two types of
multiplexers are provided: One with binary control inputs, and the other with individual control inputs.
Images:
Binary Control Inputs Individual Control Inputs
2-input 4-input 8-input 3-input 4-input
d0 d0
d0 d0 d1
d0 d1 d1 MUX MUX
MUX Y MUX Y d2
d1 d2 MUX Y d2 d3
d3
s0 d7 s1 s0 s2 s1 s0
s1 s0
s2 s1 s0
In the images, d0...d7 are the data inputs; and s0..s2 are the control inputs. The truth tables of the multiplexers
are as follows.
V V rms
2 V 12
THD = ----h- = -------------------------
-
V1 V1
where V1 is the fundamental component (rms), Vh is the harmonic rms value, and Vrms is the overall rms value
of the waveform. The THD block is modelled as shown below.
Vrms Vh
THD vin(t) THD
vin(t)
v1(t) V1
v1(t)
A second-order band-pass filter is used to extract the fundamental component. The center frequency and the
passing band of the band-pass filter need to be specified.
Attributes:
Parameters Description
Fundamental Frequency Fundamental frequency of the input, in Hz
Passing Band Passing band of the band-pass filter, in Hz
Example:
In the single-phase thyristor circuit below, a THD block is used to measure the THD of the input current. The
delay angle of the thyristor bridge is chosen as 30o. For the THD block, the fundamental frequency is set at 60
Hz and the passing band of the filter is set at 20 Hz. The simulation results are shown on the right.
vs
alpha=30 deg.
is THD
is1
One of the THD block output is the input current fundamental component is1. By comparing the phase
difference between the input voltage vs and the current is1, one can calculate the input displacement power
factor. This, together with the THD value, can be used to calculate the input power factor.
Attribute:
Parameter Description
Trigger Flag Trigger flag (0: edge-triggered; 1: level-triggered)
An edge-triggered flip-flop only changes the states at the rising edge of the set/reset input. The truth table of an
edge-triggered flip-flop is:
S R Q Q
0 0 no change
0 0 1
0 1 0
not used
A level-triggered flip-flop, on the other hand, changes the states based on the input level. The truth table of a
level-triggered set-reset flip-flop is:
S R Q Q
0 0 no change
0 1 0 1
1 0 1 0
1 1 not used
A J-K flip-flop is positive edge triggered, and the truth table is:
S R J K Clock Q Q
0 1 x x x 1 0
1 0 x x x 0 1
0 0 x x x 0 0
1 1 0 0 no change
1 1 0 1 0 1
1 1 1 0 1 0
1 1 1 1 Toggle
x: Do not care
5.4.4 D Flip-Flops
Two types of D flip-flops are provided: One without the set/reset inputs, and the other with the set input (S) and
reset input (R). For the element without the set and reset inputs, it is assumed that both set/reset inputs are high
(1).
Images:
Attribute:
Parameter Description
Pulse Width On-time pulse width, in sec.
The input node at the bottom of the controlled monostable block is for the pulse width input.
Preset Enable
Preset Value
Clock Output
Up/Down
Reset
Attribute:
Parameter Description
No. of Bits Number of bits N
When the Up/Down input is 0, the counter decrements, and when the Up/Down input is 1, the counter
increments.
The Reset input resets the counter to 0 when it is high (1). The Preset Enable input sets the counter to the preset
value when it is high.
Vin Vo Vin Vo
Clock
Let N be the number of bits. The output of the A/D converter is calculated as:
N
2
V o = ---------- V in
V ref
Attribute:
Parameter Description
Sampling Frequency Sampling frequency of the zero-order hold, in Hz
Like all other discrete elements, the zero-order hold has a free-running timer which determines the moment of
sampling. The sampling moment is synchronized with the origin of the simulation time. For example, if the
zero-order hold has a sampling frequency of 1000 Hz, the input will be sampled at 0, 1 msec., 2 msec., 3 msec.,
and so on.
Example:
In the following circuit, the zero-order hold sampling frequency is 1000 Hz. The input and output waveforms
are shown on the left.
Note that in above circuit, a continuous-domain integrator is also connected to the input sine source. This makes
it a mixed continuous-discrete circuit, and a simulation time step selected for the continuous circuit will be
used. With this time step, the familiar staircase-like waveform can be observed at the zero-order hold output.
Without the integrator, the circuit becomes a discrete circuit. Since only the calculation at the discrete sampling
points is needed, the simulation time step will be equal to the sampling period, and only the results at the
sampling points are available. The waveforms, as shown below, appear continuous. In fact the waveforms are
discrete, and the connection between two sampling points makes it look like continuous.
If a0 = 1, the expression Y(z) = H(z) * U(z) can be expressed in difference equation as:
y ( n ) = b 0 u ( n ) + b 1 u ( n 1 ) + ... + b N u ( n N )
[ a 1 y ( n 1 ) + a 2 y ( n 2 ) + ... + a N y ( n N ) ]
Image:
Attributes:
Parameters Description
Order N Order N of the transfer function
Coeff. b0...bN Coefficients of the numerator (from b0 to bN)
Coeff. a0...aN Coefficients of the denominator (from a0 to aN)
Sampling Frequency Sampling frequency, in Hz
Example:
The following is a second-order transfer function:
3
400. e
H ( z ) = ------------------------------------------------------
2
-
3
z + 1200 z + 400. e
Assuming a sampling frequency of 3 kHz, the specification will be:
Order N 2
Coeff. b0...bN 0. 0. 400.e3
Coeff. a0...aN 1. 1200. 400.e3
Sampling Frequency 3000.
Attribute:
Parameters Description
Algorithm Flag Flag for integration algorithm
0: trapezoidal rule
1: backward Euler
2: forward Euler
Initial Output Value Initial output value
Reset Flag Reset flag (0: edge reset; 1: level reset) (for external resettable integrator only)
Lower Output Limit Lower limit of the output (for internal resettable integrator only)
Upper Output Limit Upper limit of the output (for internal resettable integrator only)
Sampling Frequency Sampling frequency, in Hz
The output of an external resettable integrator can be reset by an external control signal (at the bottom of the
block). With the edge reset (reset flag = 0), the integrator output is reset to zero at the rising edge of the control
signal. With the level reset (reset flag = 1), the integrator output is reset to zero as long as the control signal is
high (1).
The output of an internal resettable integrator is reset to 0 whenever the output reaches either the lower limit or
the upper limit. The integrator works in the same way as the external resettable integrator with the edge reset,
except that users do not need to set up the external reset circuit in this case.
If we define u(t) as the input, y(t) as the output, Ts as the sampling period, and H(z) as the discrete transfer
function, the input-output relationship of an integrator can be expressed under different integration algorithms
as follows.
With trapezoidal rule:
T z+1
H ( z ) = ----s- ------------
2 z1
T
y ( n ) = y ( n 1 ) + ----s- ( u ( n ) + u ( n 1 ) )
2
With backward Euler:
z
H ( z ) = T s ------------
z1
y ( n ) = y ( n 1 ) + Ts u ( n )
Attribute:
Parameter Description
Sampling Frequency Sampling frequency, in Hz
Attribute:
Parameters Description
Gain Gain k of the PI controller
Time Constant Time constant T of the PI controller
Lower Output Limit Lower limit V_lower of the output
Upper Output Limit Upper limit V_upper of the output
Sampling Frequency Sampling frequency fs, in Hz
Attributes:
For 1st-order low-pass filter and 2nd-order low-pass filter:
Parameters Description
Gain Gain k of the filter
Cut-off Frequency Cut-off frequency fc ( f c = -----c- ), in Hz
2
Damping Ratio Damping ratio of the filter (for 2nd-order low-pass filter only)
Sampling Frequency Sampling frequency fs, in Hz
For General Digital Filter or FIR Filter that reads the coefficients directly:
Parameters Description
Order N Order N of the transfer function
Coeff. b0...bN Coefficients of the numerator (from b0 to bN)
Coeff. a0...aN Coefficients of the denominator (from a0 to aN)
Sampling Frequency Sampling frequency fs, in Hz
For General Digital Filter or FIR Filter that reads the coefficients from a text file:
Parameters Description
File for Coefficients Name of the file storing the filter coefficients
Sampling Frequency Sampling frequency fs, in Hz
y ( n ) = b 0 u ( n ) + b 1 u ( n 1 ) + ... + b N u ( n N )
[ a 1 y ( n 1 ) + a 2 y ( n 2 ) + ... + a N y ( n N ) ]
If the denominator coefficients a0..aN are not zero, this type of filter is called infinite impulse response (IIR)
filter.
The transfer function of a FIR filter is expressed in polynomial form as:
1 ( N 1 ) N
H ( z ) = b 0 + b 1 z + ... + b N 1 z + bN z
If a0 = 1, the output y and input u can be expressed in difference equation form as:
y ( n ) = b 0 u ( n ) + b 1 u ( n 1 ) + ... + b N u ( n N )
The coefficient file for block Digital Filter (file) and FIR Filter (file) has the following format:
For FIR Filter (file):
N
b0
b1
... ... ...
bN
For Digital Filter (file), the format can be either one of the following:
N or N
b0 b0, a0
b1 b1, a1
... ... ... ... ... ...
bN bN, aN
a0
a1
... ... ...
aN
Example:
To design a 2nd-order low-pass Butterworth digital filter with the cut-off frequency fc = 1 kHz, assuming the
sampling frequency fs = 10 kHz, using MATLAB, we have:
Nyquist frequency fn = fs / 2 = 5 kHz
Normalized cut-off frequency fc* = fc/fn = 1/5 = 0.2
[B,A] = butter (2, fc*)
which will give:
B = [0.0201 0.0402 0.0201 ] = [b0 b1 b2]
A = [ 1 -1.561 0.6414 ] = [a0 a1 a2]
The transfer function is:
Order N 2
Coeff. b0...bN 0.0201 0.0402 0.0201
Coeff. a0...aN 1. -1.561 0.6414
Sampling Frequency 10000.
If the coefficients are stored in a file, the file content will be:
2
0.0201
0.0402
0.0201
1.
-1.561
0.6414
Or the file can also have the content as follows:
2
0.0201, 1
0.0402, -1.561
0.0201, 0.6414
Attribute:
Parameter Description
Initial Output Value Initial value of the output of the block
Sampling Frequency Sampling frequency, in Hz
The difference between the unit delay block and the time delay block is that the unit delay block is a discrete
element and it delays the sampled points by one sampling period, whereas the time delay block is a continuous
element and it delays the whole waveform by the delay time specified.
Attributes:
Parameters Description
No. of Bits Number of bits N
Vin_min Lower limit of the input value Vin,min
Vin_max Upper limit of the input value Vin,max
Vo_min Lower limit of the output value Vo,min
Vo_max Output value Vo,max that corresponds to the input value Vin,max
Sampling Frequency Sampling frequency, in Hz
To understand how the quantization blocks work, we will look at a special case where N = 3, Vin,min = 0, and
Vin,max = 1. The input will be divided into 23 (or 8) stages. Depending on the input level, the output will be one
of the 3-bit discrete numbers 000, 001, 010, 011, 100, 101, 110, and 111. The output will have 23 (or 8) levels,
and 23-1 (or 7) steps. If the binary number 000 corresponds to the decimal number 0, and the binary number 111
corresponds to the decimal number 0.875, and Vo,min = 0, and Vo,max = 1, we can plot the waveforms of the
input, output, and the quantization error Vin-Vo of the two quantization blocks as below.
The waveforms on the left are from the quantization block without the offset, and the waveforms on the right
are from the quantization block with the offset.
(110) (110)
(101) (101)
(100) (100)
(011) (011)
(010) (010)
(001) (001)
(000) (000)
As shown from the waveforms on the left, the quantization error of the quantization block is from 0 to 1 LSB
(least significant bit), or 1/23 (or 0.125).
For the quantization block with offset, a value of 0.5*LSB is added to the input before performing the
quantization. This reduces the quantization error to -0.5 LSB to +0.5 LSB (except when the input is close to the
maximum value), as shown from the waveforms on the right.
The figure below shows the input-output relationship of a general quantization block with offset.
Vo,max N
2 Binary numbers
N
Actual output limit 2 -1
Vo
Vo
Vo,min
0
Vin,min Vin,max
Vin Vin
The number of bits determines the quantization resolution. The input range, Vin,max - Vin,min, is divided into 2N
stages, with the width of each stage as:
V o = V o, min + ( k 1 ) V o
Attributes:
Parameters Description
Buffer Length The length of the buffer
Sampling Frequency Sampling frequency, in Hz
A circular buffer stores data in a buffer. When the pointer reaches the end of the buffer, it will start again from
the beginning.
The circular buffer with single output is a type of first-in-first-out memory storage device, and the output is
equal to the value being pushed out.
The output of the circular buffer with vector output is a vector array with the length equal to the buffer length.
To access each memory location, use the memory read block.
Example:
If a circular buffer has a buffer length of 4 and a sampling frequency of 10 Hz, we have the buffer storage at
different time as follows:
where
ci = [ ak+1 * bj-k], k=0, ..., m+n-1; j=0, ..., m+n-1; i=1, ..., m+n-1
Example:
If A = [1 2 3] and B = [4 5], we have m = 3; n = 2; and the convolution of A and B is: C = [4 13 22 15].
Attribute:
Parameter Description
Memory Index Offset Offset from the starting memory location
A memory read block allows one to access the memory location of elements such as convolution block, vector
array, and circular buffer. The index offset defines the offset from the starting memory location.
Example:
Let a vector be A = [2 4 6 8]. If index offset is 0, the memory read block output will be 2. If the index offset is 2,
the output will be 6.
Attributes:
Parameters Description
Array Length The length of the data array N (for the element Array only)
Values Values of the array (for the element Array only)
File for Coefficients Name of the file storing the array (for the element Array (file) only)
If the array is read from a file, the file will have the following format:
N
a1
... ... ...
aN
where N is the length of the array, and a1..aN are the array values.
Example:
To define an array A = [2 4 6 8], we will have: Array Length = 4; Values = 2 4 6 8. If the array is to be read from
a file, the file will be:
4
2.
4.
6.
8.
5.5.9 Stack
A stack is a first-in-last-out register.
Image:
Vin
push Vo
pop
Attribute:
Parameter Description
Stack Depth The stack depth
The rising edge triggers the push or pop action. When a pop action is performed and the stack is empty, the
output remains unchanged. When a push action is performed and the stack is already full, the data at the bottom
of the stack will be pushed out and will be lost.
Vo
In PSIM, the In Link nodes receive values from Simulink, and the Out Link nodes send the values to Simulink.
They are all control elements and can be used in the control circuit only. In Simulink, the SimCoupler model
block is connected to the rest of the system through input/output ports.
File: pmsm_psim.sch
Control
in SimuLink
File: pmsm_simulink.mdl
The following are the steps to set up SimCoupler for PSIM-Matlab/Simulink co-simulation for the example
above.
Adding the SimCoupler Block to the Simulink Library:
Run the program "SetSimPath.exe" to add the SimCoupler block to the Simulink
library and set up the SimCoupler Module for co-simulation of PSIM and Matlab/Simulink. After
the execution, the SimCoupler block will appear as "S-function SimCoupler" in the Simulink
Library Browser.
Note that this step is necessary, otherwise Simulink will not be able to find PSIM.
With this, it is also not necessary to manually add the PSIM folder to the Matlab path.
Also, this program needs to be run only once. It needs to be run again only if the
PSIM folder or Matlab folder has changed.
In PSIM:
- After the rest of the power circuit is created, connect three Out Link nodes to the low-pass filters of
Phase A, B, and C currents, and rename them as Ia, Ib, and Ic; and connect one Out Link
node to the speed sensor output and rename it as Wrpm.
- Connect three In Link nodes to the positive inputs of the comparators, and rename them as Va,
Vb, and Vc.
- Go to the Simulate menu, and select Arrange SLINK Nodes. A dialog window will appear. Arrange
the order of the In Link nodes (SLINK In) and Out Link nodes (SLINK Out) to be the same as how
the input/output ports would appear in the SimCoupler model block in Simulink (the order of the
ports is from the top to the bottom). In this example, the order will be Va, Vb, and Vc for the
In Link nodes, and Ia, Ib, Ic, and Wrpm for the Out Link nodes.
- Save the schematic file. In this example, we assume that the file is saved to
C:\PSIM\pmsm_psim.sch.
In Simulink:
- Start Matlab.
- Launch Simulink. Open an existing file or create a new file. After the rest of the system is created, go
to the menu "S-function SimCoupler" in the Simulink Library Browser, select the SimCoupler
block, and place it on the schematic.
- In the PMSM example file, double click on the SimCoupler block, and click on the Browser button to
There are different ways of setting up Simulink to perform co-simulation. The recommend approach is to set the
Solve Type to Fixed-step and define the Fixed step size to be the same or close to PSIMs time step. The figure
below shows this option.
It is recommended that Simulink use the same time step as PSIM, although we have found that, even if the
Simulink time step is slightly larger than PSIM time step, satisfactory results are obtained. In this case, for
When the Simulink Solver type is set to Variable-step, in order to obtain correct results, a zero-order-hold must
be placed at the input of the SimCoupler model block. Moreover, the zero-order-hold sample time must be the
same or close to PSIM time step. The figure below shows the configuration.
Control in Simulink Solver Type: Variable-step
ZOH Sample Time: 2 us
Therefore, Simulink must be set up to have the Solver Type as Fixed-step with the time step the same or close
to the PSIM time step, or if the Solver Type is Variable-step, a zero-order-hold must be used with the sample
time the same or close to PSIM time step
In Simulink, double click on the SimCoupler block to open the property dialog, and click on Add Variable.
Click on the newly added row in the list and type the variable name and value, as shown below.
Torque Control
Id Is
Iq
Te Tes
Attributes:
Parameters Description
Ld (d-axis inductance) d-axis inductance of the machine, in H
Lq (q-axis inductance) q-axis inductance of the machine, in H
Vpk / krpm Peak line-to-line back emf constant of the machine, in V/krpm (mechanical
speed)
Number of Poles Number of poles of the machine
Integrator Gain Gain of the torque loop integrator controller
Maximum Motor Torque Maximum motor torque, in Nm
Base Current Value Base current value of the system, in A
Base Torque Value Base torque value, in Nm
Sampling Frequency Sampling frequency of the torque loop integrator controller, in Hz
The torque control block has the following inputs and outputs:
Id: d-axis current feedback
Iq: q-axis current feedback
Te: Torque reference
Is: Current amplitude reference
Tes: Estimated torque developed by the motor
This block estimates the motor developed torque from the current feedback and motor parameters. A control
loop based on a discrete integrator is used to regulate the motor torque and generate the motor current reference.
The block requires the parameters of the PMSM controlled and the inverter current rating.
Parameters Description
Ld (d-axis inductance) d-axis inductance of the PMSM machine, in H
Lq (q-axis inductance) q-axis inductance of the PMSM machine, in H
Vpk / krpm Peak line-to-line back emf constant of the machine, in V/krpm (mechanical
speed)
Number of Poles Number of poles of the machine
Maximum Motor Torque Maximum motor torque, in N*m
Maximum Motor Speed (rpm) Maximum motor speed, in rpm
Maximum Motor Power Maximum motor power, in W
Base Voltage Value Base voltage value of the system, in V
Base Current Value Base current value of the system, in A
Base Mechanical Speed Base mechanical speed, in rad/sec
Base Torque Value Maximum inverter output voltage amplitude (phase peak), in V
The torque control block has the following inputs and outputs:
Id: d-axis current feedback
Iq: q-axis current feedback
Vdc: DC bus voltage feedback
Wm: Motor mechanical speed, in rad/sec.
Tcmd: Torque command
Te: Torque reference
nmb: Calculated speed limit of the maximum torque region, in rpm
FW: Flag of field weakening (1: in field weakening; 0: not in field weakening)
This block calculates the speed limit of the maximum torque region. When the motor speed is less than this
speed limit, the motor operates in the maximum torque region. Otherwise, it operates in the maximum power
region with field weakening control.
The block requires the parameters of the PMSM controlled and the inverter voltage rating.
Parameters Description
Number of Poles Number of poles of the machine
Filter Cut-Off Frequency Cut-off frequency fc of the internal low-pass filter, in Hz
Maximum Motor Torque Maximum motor torque, in N*m
Maximum Motor Speed Maximum motor speed, in rpm
Maximum Motor Power Maximum motor power, in W
Base Voltage Value Base voltage value of the system, in V
Base Current Value Base current value of the system, in A
Base Mechanical Speed Base mechanical speed, in rad/sec
Base Torque Value Maximum inverter output voltage amplitude (phase peak), in V
Sampling Frequency Sampling frequency of the internal low-pass filter, in Hz
The dynamic torque limit control block is for the control of nonlinear PMSM machines only. It has the
following inputs and outputs, all in per unit value except Ld, Lq, and Lambda. When base values Vb, Ib, and
Wmb are all 1, all input and output quantities are in real values.
Input Signals:
Id: d-axis current feedback
Iq: q-axis current feedback
Ld: d-axis inductance Ld, in H
Lq: q-axis inductance Lq, in H
Lambda: Peak stator phase flux linkage, in Weber
Vdc: DC bus voltage feedback
Wm: Motor mechanical speed
Tcmd: Torque reference input
Output Signals:
Te: Torque reference output
Wm_th: Calculated threshold speed of the constant torque region
FW: Flag of field weakening (1: in field weakening; 0: not in field weakening)
This block calculates the threshold speed of the constant torque region. When the motor speed is less than this
speed, the motor operates in the constant torque region. Otherwise, it operates in the constant power region with
field weakening control.
A second-order low-pass digital filter, with a cut-off frequency of fc and a damping ratio of 0.7, is used to
smooth out the change of the calculated threshold speed.
Parameters Description
Vpk / krpm Peak line-to-line back emf constant of the machine, in V/krpm (mechanical
speed)
Number of Poles Number of poles of the machine
Maximum Motor Torque Maximum motor torque, in N*m
Number of Cells in Parallel Number of cells Np in parallel in the battery pack
PI Gain Gain of the voltage loop PI controller
PI Time Constant Time constant of the voltage loop PI controller, in sec.
Sampling Frequency Sampling frequency of the voltage PI controller, in Hz
The voltage control block has the following inputs and outputs:
Vdc*: DC bus voltage reference
Vdc: DC bus voltage feedback
Idc: DC bus current feedback
Wm: Motor mechanical speed in rad/sec.
Is: Current amplitude reference
This block uses a digital PI controller to regulate the dc bus voltage. Together with the dc bus current and the
machine speed, it generates the machine current reference.
The block requires the parameters of the PMSM machine controlled.
Charging Control
Vbatt
Vm
Ibatt
Attributes:
Parameters Description
Converter Rated Power Rated power of the dc-dc converter, in W
Battery-side Rated Voltage Converter rated voltage on the battery side, in V
Number of Cells in Series Number of cells Ns in series in the battery pack
Number of Cells in Parallel Number of cells Np in parallel in the battery pack
Voltage Derating Factor Voltage derating factor Ks, from 0 (100% derating) to 1 (no derating)
Full Battery Voltage Full (or maximum) voltage of the battery cell, in V
Battery Resistance Internal resistance of the battery cell, in Ohm
Current PI Gain Gain of the current loop PI controller
Discharging Control
Vdc
Vm
Ibatt
Attributes:
Parameters Description
Current/Voltage Mode Operation mode of the discharge control (0: current mode; 1: voltage mode)
Converter Rated Power Rated power of the dc-dc converter, in W
Battery-side Rated Voltage Converter rated voltage on the battery side, in V
DC Bus Voltage Reference DC bus voltage reference, in V
Current PI Gain Gain of the current loop PI controller
Current PI Time Constant Time constant of the current loop PI controller, in sec.
Voltage PI Gain Gain of the voltage loop PI controller
Voltage PI Time Constant Time constant of the voltage loop PI controller, in sec.
Control Block Output Limit The upper limit of the control block output Vm. The lower limit is 0.
Sampling Frequency Sampling frequency of the voltage and current PI controller, in Hz
The dc-dc discharging control block has the following inputs and outputs:
Vdc: DC bus voltage feedback
Regen Control
Vdc Rgn
Tes
Wm
Attributes:
Parameters Description
Regen Power Threshold Threshold of the motor power level. Beyond the threshold, regeneration is
enabled.
Regen Enable Time Time from which regeneration is enabled. Before this time, regeneration is
disabled.
Sampling Frequency Sampling frequency of the low-pass filter for speed feedback, in Hz
The dc-dc regeneration control block has the following inputs and outputs:
Vdc: DC bus voltage feedback
Tes: Estimated developed torque of the traction motor
Wm: Mechanical speed of the traction motor in rad/sec.
Rgn: Regeneration flag (0: no regeneration; 1: regeneration)
This block generates a regeneration flag based on the dc bus voltage and the traction motor power. When the
motor power is negative (in generation mode) and the power amplitude exceeds the threshold, the regeneration
flag is set to 1 to enable regeneration.
To avoid false triggering due to start-up transient, a regeneration enable time is defined. Regeneration is
possible only after this time has passed.
This chapter provides descriptions for the compoments in PSIM element librarys Other section.
Example:
The circuit below implements the step change of a load. In the circuit, the on-off switch controller is used to
control the bi-directional switch. The step voltage source, which is connected to the controller input, changes
from 0 to 1 at the time of 12 ms. The closure of the switch results in the short-circuit of the resistor across the
switch and the increase of the current.
On-off
Controller
Enable/Disable
Attributes:
Parameters Description
Frequency Operating frequency of the controlled switch/switch module, in Hz
Pulse Width On-time pulse width of the switch gating, in deg.
The input for the delay angle alpha is in deg.
Example:
The figure below shows a thyristor circuit using delay angle control. In the circuit, the zero-crossing of vs,
which corresponds to the moment that the thyristor would start conducting naturally, is used to provide the
synchronization. The delay angle is set at 30o. The gating signal is delayed from the rising edge of the
synchronization signal by 30o.
vs iRL1
vsync
Enable/Disable
Parameters Description
Frequency Switching frequency, in Hz
Update Angle Update angle, in deg., based on which the gatings are internally updated. If the
angle is 360o, the gatings are updated at every cycle. If it is 60o, the gatings are
updated at every 60o.
File Name Name of the file storing the PWM gating pattern
A lookup table, which is stored in a file, contains the gating patterns. It has the following format:
n, m1, m2, ..., mn
k1
G1,1, G1,2, ..., G1,k1
... ... ...
kn
Gn,1, Gn,2, ..., Gn,kn
where n is the number of gating patterns; mi is the modulation index correspondent to Pattern i; and ki is the
number of switching points in Pattern i. The modulation index array m1 to mn should be monotonically
increasing. The output will select the ith pattern if the input is smaller than or equal to mi. If the input exceeds
mn, the last pattern will be selected.
The following table shows an example of a PWM pattern file with five modulation index levels and 14
switching points.
5, 0.901, 0.910253, 0.920214, 1.199442, 1.21
14
7.736627 72.10303 80.79825 99.20176 107.8970 172.2634 180.
187.7366 252.1030 260.7982 279.2018 287.8970 352.2634 360.
14
7.821098 72.27710 80.72750 99.27251 107.7229 172.1789 180.
187.8211 252.2771 260.7275 279.2725 287.7229 352.1789 360.
14
7.902047 72.44823 80.66083 99.33917 107.5518 172.0979 180.
187.9021 252.4482 260.6608 279.3392 287.5518 352.0980 360.
14
10.186691 87.24225 88.75861 91.24139 92.75775 169.8133 180.
190.1867 267.2422 268.7586 271.2414 272.7578 349.8133 360.
14
10.189426 87.47009 88.97936 91.02065 92.52991 169.8106 180.
190.1894 267.4701 268.9793 271.0207 272.5299 349.8106 360.
In this example, if the modulation index input is 0.8, the controller will select the first gating pattern. If the
modulation index is 0.915, the controller will select the third pattern.
Example:
This example shows a three-phase voltage source inverter (file: vsi3pwm.sch). The PWM for the converter
uses the selected harmonic elimination. The gating patterns are described above and are pre-stored in File
vsi3pwm.tbl. The gating pattern is selected based on the modulation index. The waveforms of the line-to-line
voltage and the three-phase load currents are shown below.
Parameter Description
Gain Gain of the sensor
a a
b b
c c
Attributes:
Parameters Description
Operating Frequency Operating frequency or fundamental frequency of the ac meter, in Hz
Cut-off Frequency Cut-off frequency of the low-pass/high-pass filter, in Hz
kWh Meter Start Time Time from which the kilowatt hour meter starts, in sec. (for W/kWh meters only)
kWh Meter Stop Time Time at which the kilowatt hour meter stops, in sec. (for W/kWh meters only)
VA Display Flag Display flag for apparent power (0: no display; 1: display) (for VA-Power Factor
meters only)
PF Display Flag Display flag for power factor (0: no display; 1: display) (for VA-Power Factor
meters only)
DPF Display Flag Display flag for displacement power factor (0: no display; 1: display) (for VA-
Power Factor meters only)
In the images for the single-phase/3-phase wattmeter / kWh meter, the node with the letter "W" is for the real
power output, in watt, and the node with the letter "kWh" is for the kilowatt hour output, in kilowatt-hour.
A low-pass filter is used in the dc meter and wattmeter models to filter out high-frequency components, whereas
a high-pass filter is used in the ac meter models to filter out the dc component. The cut-off frequency determines
the transient response of the filter.
Except the voltage and current probes, the readings of all the meters are meaningful only when the readings
reach the steady state.
The single-phase and 3-phase VAR meter measure the reactive power (VAR) at the fundamental frequency; and
the VA-Power Factor meters measure the apparent power (VA), the total power factor (PF), and the
displacement power factor (DPF). The kilowatt hour meter reading gives the integration of the real power over
a time interval defined by the start time and the stop time.
For a single-phase circuit, the real power P, reactive power Q, apparent power S, total power factor PF, and the
displacement power factor DPF are defined as follows.
Assume both the voltage and current contains harmonics, i.e.
v(t) = 2 V 1 sin ( 1 t + 1 ) + 2 V 2 sin ( 2 t + 2 ) + ...
where 1 is the fundamental frequency and all others are harmonic frequencies. We have the rms values of the
Note that the reactive power is for the fundamental component only.
The apparent power S is defined as:
S = V rms I rms
The total power factor PF and the displacement power factor DPF are then defined as follow:
P
PF = ---
S
DPF = cos ( 1 1 )
A B A B C D
Interface for current scope, and
Interface for 2-channel scope 1-channel and 4-channel voltage scopes
The 1-channel and 4-channel voltage scope and the current scope have the same interface. For the 4-channel
scope, the channel can be selected through the drop-down channel selection, and the settings apply that the
channel selected.
The scope is designed to operate in a similar way as the actual oscilloscope in the lab. It has 3 main sections:
Timebase section, Channel section, and Trigger section.
In the Timebase section, the scale of time (x axis) is defined.
In the Channel section, the scale of the Y axis, as well as the offset and the color of the waveform, are defined.
The channel display mode can be either DC, AC, or Gnd (for ground). When the display mode is in DC, the
entire waveform is displayed. If the display mode is in AC, only the ac portion of the waveform is displayed.
When the display mode is Gnd, the waveform will show 0.
In the Trigger section, the trigger conditions are defined. The trigger can be set to either ON or OFF. When the
trigger is off, the waveform is free-running, and the display of the waveform in the scope may vary from one
frame to another. If the trigger is on, the waveform display will only start when the trigger conditions are met.
This will lead to a steady waveform display.
There are three trigger modes: rising-edge triggering, falling-edge triggering, and one-shot triggering (if the
once checkbox is checked, the one-shot triggering is selected). The one-shot triggering will trigger only once,
and it is useful, for example, in situations where a transient needs to be captured.
The trigger level sets the level at which the triggering occurs. For example, if Channel A is selected with the
rising-edge triggering and the trigger level of 0V, whenever the Channel A input crosses over 0 from negative to
positive, triggering will occur and the waveform display will start from that instant.
On the scope, if the Auto scale checkbox is checked, the scale of all the channels will be automatically adjusted
so that the waveforms are within the display of the scope.
Note that voltage scopes have connecting terminals which can be connected to either power circuit nodes or
control circuit nodes. The scopes will display the node-to-ground voltages at these nodes.
The current scope, on the other hand, has no connecting terminals. It is associated with any element that has the
parameter of the current flag, and it is enabled by right clicking on top of the element, and selecting the branch
current under the Current Scopes, as shown below. After the branch current is selected, a check mark will
appear in front of the branch current name.
C/P
P
k=1
Example:
The following circuit illustrates how a control circuit signal can be passed to the power circuit. As seen from the
power circuit, the CTOP block behaviors as a grounded voltage source.
Attributes:
Parameters Description
Transformation Flag Transformation flag (0: The q-axis leads the d-axis; 1: The q-axis lags the d-axis)
The angle at the bottom of the blocks is in rad.
It should be noted that, in power circuit, currents must first be converted into voltage quantities (using current-
controlled voltage sources) before they can be transformed. Also, if an input terminal is not used (such as in the
dqo-to-abc transformation block where only Phase d and q are not used, and Phase o is not used), it must be
connected to ground.
Depending on the transformation flag, the transformation equations are as follows.
From abc to dqo, with the transformation flag set to 0:
2 2
cos cos ------ cos + ------
3 3
vd va
2
vq = --- sin sin 2 2
------ sin + ------ v b
3 3 3
vo vc
1
--- 1
--- 1
---
2 2 2
Note that the coefficient 2/3 in front of the transformation matrix means that the transformation is amplitude
invariant, not power invariant. That is, the dq vector amplitude is equal to the peak amplitude of the three-phase
sinusoidal waveforms. But the power in the dq frame must be multiplied by 3/2 to have the same power, i.e.
Power = va*ia + vb*ib + vc*ic = 3/2 *(vd*id+vq*iq). To have the power invariant transformation, replace 2/3 with
2 , and 1/2 with 1.
--- ---
3 2
From dqo to abc, with the transformation flag set to 0:
cos sin 1
va 2 2 vd
cos ------ sin ------ 1
vb = 3 3 vq
vc 2 2 vo
cos + ------ sin + ------ 1
3 3
cos sin 1
va 2 2 vd
cos ------ sin ------ 1
vb = 3 3 vq
vc 2 2 vo
cos + ------ sin + ------ 1
3 3
Example:
In this example, three symmetrical ac waveforms are transformed into dqo quantities. The angle is defined as
= t where = 2*60. Since the angle changes linearly with time, a piecewise linear voltage which has a
ramp waveform is used to represent . The simulation waveforms show the three-phase ac (top), the angle
(middle), and the dqo output. In this example, the q component is constant, and both the d and the o
components are zero.
In the images, the letter "al" refers to , and the letter "be" refers to . The transformation equations are:
From abc to :
1 1 va
1 --- ---
v 2 2 2
= --- vb
v 3 3 3
0 ------- ------- vc
2 2
From ab to :
3
--- 0
v 2 v
= a
v vb
------2- 2
2
From ac to :
3
--- 0
v 2 va
=
v 2 vc
------- 2
2
From to abc:
1 0
va
2 1 - ------3- v
vb = --- --2 2
3 v
vc 1 3
--- -------
2 2
In the images, the letter "al" refers to , and the letter "be" refers to .
The transformation equations are:
From dq to :
v cos sin v d
=
v sin cos vq
In the images, the letter "r" refers to the amplitude, and "a" refers to the phase angle , in rad.
The transformation equations are:
From Cartesian to polar:
2 2
r = x +y
y
= atan ---
x
From polar to Cartesian:
x = r cos
y = r sin
Parameters Description
Expression f(x1,x2,...,xn) Expression of the output versus inputs where n is the number of inputs
Expression df/dxi Expression of the derivative of the function f versus the ith input
Parameter Description
Graph File Name of the Simview file to be used as the lookup table data
Input Column Input column of the lookup table, selected from a drop-down menu after a graph file is
loaded.
Output Column Output column of the lookup table, selected from a drop-down menu after a graph file is
loaded.
For 2-dimensional lookup tables, the node on the left is for the row index input, and the node on the top is for
the column index input.
A one-dimensional lookup table has one input and one output. Two data arrays, corresponding to the input and
output arrays, are stored in a lookup table.
A 2-dimensional lookup table with integer input has two inputs. The output data is stored in a 2-dimensional
matrix. The two inputs correspond to the row and column indices of the matrix. For example, if the row index is
3 and the column index is 4, the output will be A(3,4) where A is the data matrix.
A 2-dimensional lookup table with floating-point input is similar to the lookup table with integer input. The
difference is that inputs are floating-point values, and interpolation is used to calculate the output.
Except the one-dimensional lookup table (Simview graph) where data is defined in a Simview file, there are two
ways to define lookup table data. One is to enter the data directly in the dialog window in the Lookup Table
tab. Another is to prepare the lookup table externally using a text editor, and define the file in the dialog
window.
To enter the data directly in the dialog window, define the number of rows (and columns for 2-dimensional
lookup tables), and click on Set. Then enter the values in the data cells provided. For 2-dimensional lookup
tables with floating-point input, also enter row input array in the left most column and column input array in the
top most row.
To prepare the data externally in a text file, define the data in the format as described below. Then click on
1.1 -2 4
Row 2 2.091 3.826 4.818
2.2 3 5
6.4.5 C Block
A C block allows users to enter C code directly without compiling the code, unlike in external DLL blocks
where the code needs to be compiled into a DLL using an external compiler. The code of the C block will be
interpreted and executed at runtime by a built-in C interpreter in PSIM.
The interface of the C block dialog window is shown below.
Input/output
ports
Area for
custom code
2-input / 3-output
in[0] out[0]
in[1] out[2]
Input/output
ports
Area for
custom code
1 1 1 1
2 2
3 3
6 6
input output
Attribute:
Parameter Description
File Name Name of the DLL file
The node with a dot is for the first input (in[0]). The input/output node sequence is from the top to the bottom.
The images and parameters of a general DLL block are shown below.
Image (for a block with 2 inputs and 3 outputs):
1
1
input 2 output
2 3
Parameter Description
DLL File Name of the DLL file
Input Data File Name of the input data file that DLL reads (optional)
Number of Input Nodes Number of input nodes (optional)
Number of Output Nodes Number of output nodes (optional)
IN Nodes List of input nodes (optional)
OUT Nodes List of output nodes (optional)
Parameter 1 Parameter to be passed from PSIM into DLL (optional)
Parameter 2 Parameter to be passed from PSIM into DLL (optional)
Edit Image (button) Click this button to edit and customize the image of the DLL block.
Display File (button) Click this button to display the content of the Input Data File (optional).
Read File (button) If the Input Data File is modified, click this button to reload the data file
(optional).
The node with a dot is the first input (in[0]). The input/output node sequence is from the top to the bottom.
By default, users define the number of inputs and outputs. But the number of inputs and outputs, the node
names, as well as the number of parameters and the parameter names can all be defined inside the DLL routine.
For more details on defining and programming for the general DLL block, please refer to the help file "Help
General DLL Block.pdf" and related examples.
The name of the DLL file can be arbitrary. The DLL file can be placed in one of the two places, in the order of
precedence: in the PSIM directory, and in the same directory as the schematic file that uses the DLL file.
Note: When a DLL file is used in multiple DLL blocks in a schematic, if global or static variables are declared
and used in the DLL code, these global/static variables will be the same and be shared among all the DLL
blocks. If this is not what users intended, the DLL calculation can be incorrect, and users should avoid using
global/static variables in the code in this case.
Several examples using the simple DLL blocks and the general DLL blocks are provided in the
"examples\custom DLL" sub-folder in the PSIM directory.
Parameter Description
DLL File Name of the DLL file that defines the functionality and the interface of the block
Number of Nodes Total number of input and output nodes
An Embedded Software Block has similar functionality as a general external DLL block. However, unlike the
general DLL block whose connection nodes are predefined as either inputs or outputs, the Embedded Software
Block allows the node types to be programmed as needed. Also, additional information, such as the exact
instant at which the state of a variable changes, can be calculated and passed to and from PSIM.
The Embedded Software Block is a control circuit element, and can be used in the control circuit only.
For more information on how to use the Embedded Software Block, please refer to the document "Help
Embedded Software Block.pdf".
6.5.1 PWM IC
The following PWM IC models are provided:
Images:
UC3823A/B UC3825A/B UC3842/3843 UC3844/3845 UC3846 UC3854
1 16 1 16 1 8 1 8 1 16 1 16
INV VREF INV VREF COMP VREF COMP VREF ILIM STDN GND GDRV
NI VCC NI VCC VFB VCC VFB VCC VREF VIN PKLM VCC
EA OUTB EA OUTB IS OUT IS OUT CS- BOUT CAO CT
CLK VC CLK VC GND IS
4 Rt/Ct 5 4
Rt/Ct GND 5 CS+ VC SS
RT PGND RT PGND NI GND MOUT RSET
CT OUTA CT OUTA INV AOUT IAC VS
RAMP GND RAMP GND COMP SYNC VAO ENA
SS ILIM SS ILIM CT RT VRMS VREF
8 9 8 9 8 9 8 9
Attribute:
Parameter Description
Model Level All these IC models have two levels:
- Level 1: The PWM gating outputs are control signals with logic level of 1
(high) and 0 (low)
- Level 2: the PWM gating outputs are power signals which can directly
drive three-state switches.
PWM IC Descriptions:
UC3823A/B, High speed PWM controller for high frequency switch-mode power
UC3825A/B supplies. The UC3823A and UC3823B and the UC3825A and UC3825B
family of PWM controllers are improved versions of the standard UC3823
and UC3825 family.
UC3842/3843, UC3844/3845 Current mode fixed frequency PWM controller for off-line or dc-dc
converters. The UC3842/3/4/5 family of control devices provides the
necessary functions to implement off-line or dc-to-dc fixed frequency
current mode control schemes.
UC3846 Current mode fixed frequency PWM controller for off-line or dc-dc
converters. The UC3846 family of control ICs provides the necessary
functions to implement fixed frequency, current mode control schemes.
UC3854, UC3854A/B High power factor preregulator providing active power factor correction.
UC3854A/B products are pin compatible enhanced versions of the
UC3854. It provides all the functions necessary for active power factor
corrected preregulators.
6.5.2 Driver IC
PSIM has implemented the following power switching devices driver models:
Images:
IR21834 IRS21867 TC4423A TC4425A TC4427A
1 14 1 8 1 8 1 8 1 8
HIN Vcc VB NC NC NC NC NC NC
LIN VB HIN HO INA OUTA INA OUTA INA OUTA
VSS HO LIN Vs GND VDD GND VDD GND VDD
DT VS COM LO INB OUTB INB OUTB INB OUTB
4 5 4 5 4 5 4 5
COM
LO TC4424A TC4426A TC4428A
1 8 1 8 1 8
Vcc NC NC NC NC NC NC
7 8
INA OUTA INA OUTA INA OUTA
GND VDD GND VDD GND VDD
INB OUTB INB OUTB INB OUTB
4 5 4 5 4 5
Driver IC Description:
IR21834 Half-bridge driver for high voltage, high speed power MOSFET and IGBT
with dependent high and low side referenced output channels.
IRS21867 High and low side driver for voltage, high speed power MOSFET and
IGBT driver with independent high and low side referenced output
channels.
TC4423A, TC4424A, TC4425A 3A dual output high-speed power MOSFET drivers.
TC4423A: dual inverting.
TC4424A: dual non-inverting.
TC4425A: one inverting and one non-inverting.
TC4426A, TC4427A, TC4428A 1.5A dual output high-speed power MOSFET drivers.
TC4426A: dual inverting.
TC4427A: dual non-inverting.
TC4428A: one inverting and one non-inverting.
Attribute:
Parameter Description
Initial Value Initial value of the node voltage
A parameter file is a text file created by a text editor. It supports the following format:
// Text from the character // to the end of the line is treated as comment
k1 = 12.3 // Define the value of the variable k1
k1 = a + b*c // Define k1 in a math expression
(global) k1 = 12.3// Define k1 as a global variable (for use in SimCoder)
L1 = 3m // power-of-ten suffix is supported. L1 = 3e-3
C1 = 100uF // C1 = 100e-6. The unit "F" is ignored.
Note that, as compared to Version 9.1, the following format is no longer supported:
LIMIT var1 Vlow Vupper // This should be replaced by the "if...else if..." statement
% as comment. Now comments must start with double slash //. % is now used as modulo operation.
var1 value (for example "R1 2.5"). Now an equal sign must be used to define a variable.
6.8 AC Analysis
6.8.1 AC Sweep
The frequency response of a circuit or a control loop can be obtained with the ac analysis. A key feature of the
ac analysis in PSIM is that, a circuit can be in its original switchmode form, and no average model is required.
Nevertheless, with the average model, it takes much shorter time to perform the ac analysis.
The following are the steps to set up the ac analysis:
- Identify a sinusoidal voltage source as the excitation source for the ac sweep.
- Place ac sweep probes at the desired output location. To measure the loop response of a closed control
loop, use the node-to-node probe.
- Place the AC Sweep block on the schematic, and define the parameters of the ac sweep.
- Run the simulation.
Attributes:
Parameters Description
Start Frequency Start frequency of the ac sweep, in Hz
End Frequency End frequency of the ac sweep, in Hz
No. of Points Number of data points
Flag for Points Flag to define how the data points is generated.
Flag = 0: Points are distributed linearly in LOG10 scale
Flag = 1: Points are distributed linearly in linear scale
Source Name Name of the excitation source
Start Amplitude Excitation source amplitude at the start frequency
End Amplitude Excitation source amplitude at the end frequency
Freq. for extra Points Frequencies of additional data points. If the frequency-domain characteristics
change rapidly at a certain frequency range, one can add extra points in this region
to obtain better data resolution.
The principle of the ac analysis is that a small ac excitation signal is injected into the system as the perturbation,
and the signal at the same frequency is extracted at the output. To obtain accurate ac analysis results, the
excitation source amplitude must be set properly. The amplitude must be small enough so that the perturbation
stays in the linear region. On the other hand, the excitation source amplitude must be large enough so that the
output signal is not affected by numerical errors.
In general, a physical system has low attenuation in the low frequency range and high attenuation in the high
frequency range. A good selection of the excitation source amplitude would be to have a relatively small
amplitude at the low frequency, and a relatively large amplitude at the high frequency.
Sometimes, after ac analysis is complete, a warning message is displayed as follows:
Warning: The program did not reach the steady state after 60 cycles. See File message.txt for more details.
This message occurs when the software fails to detect the steady state at the ac sweep output after 60 cycles. To
address this problem, one may increase damping in the circuit (by including parasitic resistances, for example),
or adjust the excitation source amplitude, or reduce simulation time step. The file message.txt gives the
information on the frequency at which this occurs and the relative error. The relative error will indicate how far
the data point is from reaching the steady state.
Example: Open-Loop Response of a Buck Converter
The circuit on the left is an one-quadrant buck converter. An excitation source is injected to the modulation
signal, and the output voltage is measured. The result of the ac analysis, on the right, shows the open-loop
response of the output voltage versus the modulation signal.
AC
Sweep
m-sine
Attributes:
Parameters Description
Start Frequency Start frequency of the ac sweep, in Hz
End Frequency End frequency of the ac sweep, in Hz
No. of Points Number of data points
Flag for Points Flag to define how the data points is generated.
Flag = 0: Points are distributed linearly in LOG10 scale
Flag = 1: Points are distributed linearly in linear scale
Source Name Name of the excitation source
Source Amplitude Excitation source peak amplitude at the start frequency.
Number of Cycles The number of cycles of the excitation source amplitude to be simulated
Steady State Time The time predicted when the circuit has completed transient and reached steady
state, in second. User may make approximated estimation of this value or run a
single time simulation to determine it.
Attributes:
Parameters Description
Name Name of the parameter to be swept
Start Starting value of the parameter
End End value of the parameter
Increment Increment step value of the parameter
Enabled When checked, the parameter line will be enabled.
For example, let the resistance of a resistor be Ro. To sweep the resistance from 2 Ohm to 10 Ohm, with an
increment step value of 2 Ohm, check the Enabled checkbox, and define the following:
Name Ro
Start 2
End 10
Increment 2
Parameter sweep will produce two plots: One is the outputs versus time, and the other is the outputs at the last
simulation point versus the swept parameter.
For example, a circuit has two outputs V1 and V2, and the resistance Ro is swept. The total simulation time is
0.1 sec. After the simulation, there will be two plots in Simview. One is V1 and V2 versus time. The other is V1
and V2 versus Ro. The V1 and V2 values used in the second plot are the values at the last simulation point at 0.1
sec.
Several types of independent voltage/current sources are available in PSIM. The notation of a current source
direction is: the current flows out of the higher-potential node, through the external circuit, and back into the
lower-potential node of the source.
Note that current sources can be used in the power circuit only.
7.1 Constant
The constant element allows users to define a constant. It behaviors as a grounded voltage source.
Image:
Attribute:
Parameter Description
Value Value of the constant
7.2 Time
The Time element is a special case of the piecewise linear voltage source. It is treated as a grounded voltage
source, and the value is equal to the simulation time, in sec.
Image:
7.3 Ground
There are three different ground elements in PSIM library. Although they are in different shape, they are
electrically connected together. The different shape of these grounds provide the convenience for user to
separate the grounding in different functional sections of the circuits.
Images:
7.4.1 DC Source
A dc source has a constant amplitude. The reference of the grounded dc voltage sources is the ground.
Images:
DC DC (battery) Grounded DC (circle) Grounded DC (T) Current
Attribute:
Parameter Description
Amplitude Amplitude of the source
Series Resistance Source series resistance, in Ohm (for DC and DC (battery) voltage sources only)
Series Inductance Source series Inductance, in H (for DC and DC (battery) voltage sources only)
The series resistance and inductance represent the source impedance.
t
1/f
Images:
Voltage Current
Attributes:
Parameters Description
Peak Amplitude Peak amplitude Vm
Frequency Frequency f, in Hz
Phase Angle Initial phase angle , in deg.
DC Offset DC offset Voffset
Attributes:
Parameters Description
V (line-line rms) Line-to-line rms voltage amplitude
Frequency Frequency f, in Hz
Initial Angle (phase A) Initial angle for Phase A, in deg.
Series Resistance Source series resistance, in Ohm
Series Inductance Source series inductance, in H
Attributes:
Parameters Description
Vpeak-peak Peak-to-peak amplitude Vpp
Frequency Frequency, in Hz
Duty Cycle Duty cycle D of the high-potential interval
DC Offset DC offset Voffset
Phase Delay Phase delay of the waveform, in deg.
D*T Voffset
0 t
T
(T=1/f)
When the phase delay is positive, the waveform is shifted to the right along the time axis.
Attributes:
Parameters Description
Vpeak-peak Peak-to-peak amplitude Vpp
Frequency Frequency, in Hz
Duty Cycle Duty cycle D of the rising slope interval
DC Offset DC offset Voffset
Phase Delay Phase delay of the waveform, in deg.
D*T Vpp
0 t
Voffset
T
(T=1/f)
When the phase delay is positive, the waveform is shifted to the right along the time axis.
A sawtooth-wave voltage source is a special case of the triangular source with the duty cycle of 1 and dc offset
and phase delay of 0, and with one node grounded. It is defined as below:
Image:
Sawtooth-wave
Parameters Description
Vpeak Peak amplitude
Frequency Frequency, in Hz
Voltage Current
Attributes:
For the Step type source:
Parameters Description
Vstep Value Vstep after the step change
Tstep Time Tstep at which the step change occurs
Parameters Description
Vstep1 Value Vstep1 before the step change
Vstep2 Value Vstep2 after the step change
Tstep Time Tstep at which the step change occurs
T_transition Transition time Ttransition from Vstep1 to Vstep2
Vstep1
Ttransition
0 Tstep t 0 Tstep t
Attributes:
For the sources that define the values and times separately:
Parameters Description
Frequency Frequency of the waveform, in Hz
No. of Points n No. of points
Values V1...Vn Values at each point
Time T1...Tn Time at each point, in sec.
For the sources that define the values and times in pair:
Parameters Description
Frequency Frequency of the waveform, in Hz
Times, Values (t1,v1) ... Time and value at each point
The time and value pair must be enclosed by left and right brackets. The time and value can be separated by
either a comma, such as (1.2m,5.5), or a space, such as (1.2m 5.5), or both, such as (1.2m, 5.5).
Example:
The following is a non-periodic piecewise linear source. It has 3 segments which can be defined by four points
(marked in the figure).
3
2
1
Time (sec.)
0 0.1 0.2 0.3
The specification for the piecewise linear voltage source will be:
Frequency 0.
No. of Points n 4
Values V1...Vn 1. 1. 3. 3.
Times T1...Tn 0. 0.1 0.2 0.3
The specification for the piecewise linear (in pair) voltage source will be:
Frequency 0.
Times, Values (t1,v1)... (0., 1) (0.1, 1) (0.2, 3) (0.3, 3)
where Vm is the peak-to-peak amplitude of the source, n is a random number in the range of 0 to 1, and Voffset is
the dc offset.
Images:
Voltage Current
Attributes:
Parameters Description
Peak-Peak Amplitude Peak-to-peak amplitude of the source
DC Offset DC offset
Attributes:
Parameters Description
Expression The mathematical expression of the source
Tstart Start time of the source
In the expression, T or t represents time. For example, to implement a sinusoidal source, the expression will
be: sin(2*3.14159*60*t+2.09).
vin1 vin2
Current-controlled Variable-gain
Voltage-controlled Current-controlled voltage-controlled
(flowing through)
vin1 vin2
Attribute:
Parameter Description
Gain Gain of the source
For voltage-controlled sources, the controlling voltage is from the positive node (+) to the negative node (-).
For current-controlled sources (with an arrow pointing from one node to another), the control nodes are
connected across a RLC branch, and the direction of the controlling current is indicated by the arrow.
For current-controlled sources (with a wire connecting the two nodes), the controlling current flows from one
control node to another, as indicated by the arrow. A 10-uOhm resistor is used to sense the controlling current.
The output of a controlled source, except variable-gain controlled sources, is equal to the gain multiplied by the
controlling voltage or current. For the variable-gain controlled sources, the output is equal to the following:
v o = ( k v in 2 ) v in 1
i o = ( k v in 2 ) v in 1
Input 1 is on the side with the multiplication sign, and Input 2 is on the side with the letter k.
The difference between a variable-gain controlled source and a nonlinear source with multiplication is that, for
the nonlinear source with multiplication, values of both vin1 and vin2 at the current time step are used to calculate
the output and are updated in each iteration. But for the variable-gain controlled source, it is assumed that the
change of vin2 is small from one time step to the next, and the value of vin2 at the previous time step is used at
the current time step. This assumption is valid as long as vin2 changes at a much slower rate as compared to vin1
and the time step is small as compared to the change of vin2. Variable-gain controlled sources can be used in
circuits which may otherwise have convergence problem with nonlinear sources with multiplication.
Example:
The circuits below illustrates the use of current controlled voltage sources. In the circuit on the left, the voltage
source is controlled by the inductor current is. With a gain of 1, the waveform of the voltage vis is equal to that
of is. In this way, a current quantity can be converted to a voltage quantity.
Vis Vis
is is
In the nonlinear power source, the term sign(vin) is 1 if vin is positive, and it is -1 if vin is negative.
Note that these nonlinear sources can be used in the power circuit only.
Images:
vin1 vin2
vin1 vin2
Attributes:
For all the sources except the nonlinear power source:
Parameter Description
Gain Gain k of the source
For the nonlinear power source:
Parameters Description
Gain Gain k of the source
Coefficient k1 Coefficient k1
Coefficient k2 Coefficient k2
Q0 Q0 Q1
Q1
clock
clock
1V
1V
In the circuit on the left, the initial values of both Q0 and Q1 are assumed to be zero. At the clock rising edge, Q0
will change to 1. Without delay, the position of Q1, which should remain at 0, will toggle to 1 at the same time.
To prevent this, a time delay element with the delay period of one time step needs to be inserted between Q0 and
the input (J) of the second flip-flop.
It should be noted that, in PSIM, the power circuit and the control circuit are solved separately. There is one
time step delay between the power and the control circuit solutions.
W-2 Warning!!! The program did not reach the steady state after 60 cycles when performing the ac
sweep.
This warning occurs when the program fails to reach the steady state after 60 cycles when performing
the ac sweep. The cause of the problem could be that the system is poorly damped at that particular
frequency or the signal amplitude is too small.
You may try the following to isolate and solve the problem:
- Run the time-domain simulation with the excitation source at that frequency and see if time-
domain waveforms are oscillatory.
- Increase the excitation voltage amplitude for larger signal level, or
- Reduce the time step for better accuracy and resolution.
2.1 Debugging
Some of the approaches in debugging a circuit is discussed in the following.
Symptom:
Simulation results show sudden changes (discontinuity) of inductor currents and capacitor voltages.
Solution:
This may be caused by the interruption of inductor current path and short-circuit of capacitor (or
capacitor-voltage source) loops. Check the switch gating signals. If necessary, include overlap or
dead time pulses to avoid open-circuit or shooting-through.
If an initial current is assigned to an inductor, initial switch positions should be set such that a path is
provided for the current flow. Otherwise, the inductor current will be forced to start from zero.
Symptom:
Simulation waveforms look incorrect or inaccurate, or the waveform resolution is poor.
Solution:
This may be caused by two reasons. One is the time step. Since PSIM uses the fixed time step during
the entire simulation, one should make sure that the time step is sufficiently small. As a rule of
thumb, the time step should be several tens times smaller than the switching period.
Another reason is the problem of waveform display. One should make sure that the print step is not too
big. To display all the data points, set the print step to 1.
Index 235
parameter 6, 16, 213, 215 s-domain
filter external resettable 150
band-pass 155, 165 internal resettable 150
band-stop 155 z-domain
digital 174 external resettable 172
high-pass 155, 197 internal resettable 172
low-pass 152, 155, 174, 183, 197 interface
second-order 155 between power and control circuits 200
flag interface block 116
load 4 control-power 200, 231
master/slave 80, 83, 84, 86, 87, 90, 92, 95, 97, 99, 114, mechanical-electrical 78, 112, 115
115
save 4 J
flip-flop JMAG 1, 104
D 167 JMAG-RT 1, 108, 109
J-K 167
S-R 166 L
format label 9, 11, 28, 38
ASCII text 33 leakage flux 63, 69, 70, 88
binary 21, 33 LED 50
free run 4, 19 library
functional model 138, 141 adding new DLL element 28
adding new subcircuit element 27
G edit 11, 25, 26, 28, 29
gate library browser 8, 11, 183
AND 166 light intensity 139
logic 3, 166, 231 limiter
NAND 166 gradient (dv/dt) 160
NOR 166 lower 160
NOT 166 lower/upper 160
OR 166 range 160
XOR 166 upper 160
gating block 52, 54, 58, 59, 60 LOG 6, 35, 157
gear box 78, 112, 114, 115 LOG10 6, 35, 157, 216, 218
logic gate 166
H lookup table
hardware target 2-D 205
PE-Pro/F28335 5 loss calculation 120, 121
TI F28335 5 diode 123
HEV 1, 186 IGBT 125, 127, 131
HEV Design Suite 186 MOSFET 129
HEV generator 12
HEV traction motor 12 M
hybrid electric vehicle 1 machine
brushless dc 85, 87, 120
I dc 84, 85
image induction 77, 80, 82, 116
secondary 26 nonlinear induction 80
Image Editor 11, 29 nonlinear permanent magnet synchronous 94
inductor 28, 41, 106, 233 nonlinear switched reluctance 99
coupled 46, 65 permanent magnet synchronous 92, 94
saturable 45 squirrel-cage induction 80, 81, 82
initial value 47, 48, 149, 150, 151, 213 switched reluctance 97
integrator synchronous (external excitation) 89
236 Index
wound-rotor induction 80, 82, 83 enter 21
MagCoupler 105 PE-Pro/F28335 5
MagCoupler block 104, 105, 106, 107 PHEV 12
MagCoupler Module 104, 107, 108, 112 physical model 140
MagCoupler-DL block 104 pitch angle 144
MagCoupler-RT Module 108, 112 port
magnetic elements 24, 69, 72 bi-directional 13
math function block 204, 205 input signal 13, 158
Matlab/Simulink 1, 18, 24, 182, 183, 184 output signal 13
maximum power 138 subcircuit 14, 15
maximum/minimum function block 159 power coefficient 144, 145
Maximum-Torque-Per-Ampere Control 100 power function block 157
mechanical coupling block 112, 115 print step 4, 232, 233
mechanical load probe 196
constant-power 112 ac sweep 217
constant-speed 113 ac sweep (loop) 216
constant-torque 112 current 12, 18, 106, 107, 118, 196
externally-controlled 114 voltage 18, 196
general-type 113 propagation delay 163, 231
mechanical load model 1, 77, 112, 116 proportional-integral controller 152
memory read block 179, 180 Pulse 168
message PWM
error 22, 208, 232 space vector 165
warning 2, 18, 22, 216, 232
meter Q
3-phase VAR 198 quantization block 1, 177
VA-power factor 196
VAR 196, 198 R
monostable 168 reference direction
Motor Drive Module 1, 77, 112 of mechanical system 78, 79, 85, 112, 115, 117, 118
multiplexer 163 Regeneration Control 192
multiplier 3, 156 relay 77
multi-rate sampling system 182 renewable energy 145
resistor 6, 13, 24, 25, 41, 45, 70, 71, 72, 196, 200, 213,
N 219, 228, 233
node resolver 118, 119, 120
link 182, 183 rheostat 45
SLINK 18, 183 root-mean-square function block 157
nonlinear element 47, 48 round-off block 162
conductance-type 47, 48 rubber band 21
conductance-type with additional input 47 runtime graph 18, 21, 41, 44
resistance-type 47
resistance-type with additional input 47 S
sampling/hold block 161
O scope
operational amplifier current 4, 11, 18, 19, 198, 199, 200
ideal 73 voltage 4, 18, 198
non-ideal 74 sensor
optocoupler 75 current 196
hall effect 85, 86, 87, 88, 89, 120
P position 77, 86, 87, 88, 89, 104, 108, 112, 118, 119,
Parameter Tool 24, 215 120
password 10 speed 116, 117, 118, 183
disable 21 torque 116, 117, 118
Index 237
voltage 196, 200, 231 T
shunt regulator 75 THD block 164, 165
sign function block 6, 35, 158 Thermal Module 120, 122, 127
SimCoder 2, 13, 18, 22, 214 TI F28335 5
SimCoupler Module 3, 18, 24, 182, 183, 184 time 209, 221
simulation control 4, 17, 19 idle 22
solar cell 139 print 4
source 221 total 4, 17, 20
constant 221 time delay block 162, 163, 176, 231
current-controlled current 227 time step 4, 5, 17, 18, 20, 76, 104, 105, 184, 209, 231
current-controlled voltage 106, 227, 228 tip speed ratio 144, 145
dc 52, 53, 222 TL431 75
math function 227 toolbar 8, 11, 21, 23, 37
nonlinear voltage-controlled 229 Torque Control 186
piecewise linear 221 transconductance 57, 128, 129, 131
random 226 transfer function block
sawtooth 224 s-domain 149
sinusoidal 222, 223 z-domain 171
square wave 162 transformation 100, 186, 201
square-wave 223 abc-alpha/beta 202
step 225 abc-dqo 100, 101, 201
triangular 224 alpha/beta-dq 203
voltage-controlled current 106, 227 Cartesian-Polar 204
voltage-controlled voltage 106, 227 transformer 3, 63
Space 165 ideal 63
square waveform block 161 saturation 68
square-root function block 157 single-phase 64
stack 181 three-phase 65
subcircuit 10, 12, 13, 15, 17 transistor
connecting 15 BJT (npn) 48, 49, 53
creating BJT (pnp) 48, 49
in the main circuit 14 GTO 48, 56
in the subcircuit 14 IGBT 48
image 16 MOSFET (n-channel) 48, 49
summer 156 MOSFET (p-channel) 48, 49, 54
sweep trapezoidal waveform block 161
ac 215, 216, 218, 233 trigonometric function 158
parameter 219
switch U
bi-directional 48, 56, 193 unit delay block 163, 176
DIAC 48, 51
linear 49, 57, 58 V
self-commutated 48, 53, 54, 55
variable
thyristor 48, 52, 193
default 13, 16, 27
transistor 48, 124, 125, 128, 129, 130
passing 15
TRIAC 48, 52
runtime 11, 19
switch controller 3, 52, 58, 193
voltmeter
alpha 52, 60, 62, 63, 193
ac 196
on-off 53, 55, 193
dc 196
PWM lookup table 63, 194
switch module W
single-phase 60
wattmeter 196, 197, 198
three-phase 61
3-phase 198
wind turbine 143, 144
238 Index
wire 9, 11, 14, 22, 232
Z
zero-order hold 1, 161, 162, 170, 182, 209
Index 239