A DSP-based AC Electronic Load for Unintentional Islanding Tests by Wei Feng B.Sc., Tsinghua University, Beijing, China, 2004 M.Sc., Institute of Electrical Engineering, CAS, China, 2007 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF MASTER OF APPLIED SCIENCE in The Faculty of Applied Science (Electrical and Computer Engineering) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) August 2009 c Wei Feng 2009 Abstract The distributed generation (DG) system for the renewable energy resources is being used more and more widely around the world in the past decades. In order to protect the whole power system from the islanding situation, the interface inverters of the DG system have to pass the unintentional islanding test. In this study, a single phase AC electronic load is introduced to simulate the local loads under unintentional islanding conditions. This thesis proposes a proper schematic of the electronic loads based on the H-bridge dc-dc converter and analyzes the performance of this system with PowerSIM. Then, a set of specifications for different ranges of the electronic loads are calculated through theoretical formulas and the PSIM simulations. Meanwhile, the transfer function of the system is also derived to analyze the stability of the PI control system. According to the theoretical analysis and simulations, a control program is implemented based on the Texas Instruments (TI) Digital Signal Processor (DSP) TMS320F2407A for the different kinds of electronic loads. A test circuit is then built to validate the performance of the system. Some experiments are performed for the resitive, inductive and capacitive loads respectively. ii Table of Contents Abstract . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . ii Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . iii List of Tables vi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . List of Figures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . vii Acknowledgements . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 1 Introduction . . . . . . . . . . . . . . . . . . . . 1.1 Project Background . . . . . . . . . . . . . . 1.2 System Overview . . . . . . . . . . . . . . . 1.2.1 Unintentional Islanding Test Scheme 1.2.2 AC Electronic Loads for Unintentional 1.3 Thesis Objective . . . . . . . . . . . . . . . . 2 2 3 3 4 4 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Islanding Test . . . . . . . . . 2 Principles of the Simulated Loads . . . . . . . . . . . 2.1 Schematic of Electronic loads . . . . . . . . . . . . . . 2.1.1 H-Bridge Bidirectional Current Source . . . . 2.1.2 The Range of Inductor Current . . . . . . . . 2.2 AC Small-Signal Modeling and PI Control Method . . 2.2.1 AC Small-Signal Model . . . . . . . . . . . . . 2.2.2 The Transfer Function of Control System . . . 2.2.3 Control System Stability Analysis . . . . . . . 2.3 System Performance Analysis . . . . . . . . . . . . . . 2.3.1 Current Ripple . . . . . . . . . . . . . . . . . . 2.3.2 Relationship between Response Time and Zsim 2.3.3 Conclusions of System Performance . . . . . . 2.4 Calculation of Circuit Parameters . . . . . . . . . . . 2.4.1 Introduction of PSIM Simulation Diagram . . 2.4.2 DC Bus Voltage: Vdc . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 6 7 8 9 10 11 12 13 13 17 17 18 18 19 iii Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 20 20 22 25 27 28 28 32 35 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 36 36 37 38 41 44 45 46 47 48 4 Software Implementation . . . . . . . . . . . . 4.1 Program Overview . . . . . . . . . . . . . . . 4.2 Variables Normalization . . . . . . . . . . . . 4.2.1 Q Format . . . . . . . . . . . . . . . . 4.2.2 Normalization of Voltage and Current 4.3 System Initialization . . . . . . . . . . . . . . 4.3.1 Event Manager Modules Initialization 4.3.2 ADC Module Initialization . . . . . . 4.3.3 I/O Ports Module Initialization . . . 4.4 Interrupt Service Routine (ISR) . . . . . . . 4.5 Generation of Sine Wave . . . . . . . . . . . 4.6 Digital PI Control . . . . . . . . . . . . . . . 4.7 RLC Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49 49 51 51 51 52 52 53 54 54 55 57 58 . . . . . . . . . . . . . . . . . . . . . . . . . . . at L = 2.6mH . . . . . . . . . . . . . . . . 61 61 62 62 2.5 2.4.3 Series Resistor: R1 . . . . . . . 2.4.4 Switching Frequency: fpwm . . . 2.4.5 Inductor: L . . . . . . . . . . . 2.4.6 Calculation of the Range of Zsim 2.4.7 Duty Ratio D . . . . . . . . . . 2.4.8 Specifications of Electronic Loads Simulations of Electronic Loads . . . . 2.5.1 Resistive Load Simulation . . . 2.5.2 Inductive Load Simulation . . . 2.5.3 Capacitive Load Simulation . . 3 Hardware Implementation 3.1 Introduction of System . 3.2 IGBT Inverter Module . 3.3 IGBT Driver Circuit . . . 3.4 Voltage Measuring Circuit 3.5 Current Measuring Circuit 3.6 Temperature Sensor . . . 3.7 Protection Circuits . . . 3.8 DSP Control Board . . . 3.9 Power Supply . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5 Experimental Results . . . . . . . . . . 5.1 Experimental Conditions . . . . . . . 5.2 Resistive Load Results . . . . . . . . 5.2.1 The Range of Resistance Loads . . . . . . . . . . . . . . . iv Table of Contents . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64 65 67 69 . . . . . . . . . . . . . . . . . . . . . . . . . . . . 71 71 72 73 Bibliography . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 75 5.3 5.4 5.2.2 Switching Frequency and Current Ripple 5.2.3 Series Inductor L and Current Ripple . . Inductive Load Results . . . . . . . . . . . . . . Capacitive Load Results . . . . . . . . . . . . . 6 Conclusions and Future Work 6.1 Synopsis . . . . . . . . . . . 6.2 Conclusions . . . . . . . . . . 6.3 Suggestions for Future Work . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Appendices A Schematics . . . . . . . . . . . . A.1 Schematics of Control Board A.2 Schematics of Power Board . A.3 PCB Layouts . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76 76 81 83 B Schematics for 3-Phase Islanding Test . . . . . . . . . . . . . 89 C Unintentional Islanding Test Conditions 91 . . . . . . . . . . . v List of Tables 2.1 2.2 2.4 2.5 2.6 2.7 2.8 2.9 Inductor Values for the Resistive Simulated Loads . . . . . . Inductor Values for the Inductive or Capacitive Simulated Loads . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Inductor Values for the Inductive or Capacitive Simulated Loads (Rounding) . . . . . . . . . . . . . . . . . . . . . . . . Minimum Rsim for the Given Series Inductor . . . . . . . . . Minimum Lsim (or Csim ) for the Given Series Inductor . . . . The Resolution of Rsim . . . . . . . . . . . . . . . . . . . . . The Resolution of Lsim (or Csim ) . . . . . . . . . . . . . . . . Specifications for the Resistive Simulated Loads . . . . . . . . Specifications for the Inductive or Capacitive Simulated Loads 21 25 25 26 27 28 28 3.1 Values of the dead times of RC networks . . . . . . . . . . . . 40 5.1 Comparisons of Current Ripple between Simulation Results and Experimental Results . . . . . . . . . . . . . . . . . . . . Comparisons of Current Ripple at 20kHz . . . . . . . . . . . Comparisons of Current Ripple at 20kHz . . . . . . . . . . . 64 65 66 2.3 5.2 5.3 20 21 vi List of Figures 1.1 1.2 Generic System for Unintentional Islanding Test . . . . . . . Electronic Loads for Unintentional Islanding Test . . . . . . . 3 4 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 2.10 2.11 2.12 2.13 2.14 2.15 2.16 2.17 2.18 2.19 2.20 2.21 2.22 2.23 2.24 2.25 2.26 2.27 2.28 Diagram of Electronic Loads . . . . . . . . . . . . . . . . . . Controllable Current Source with Regulated DC Bus Voltage H-Bridge Inverter With Regulated DC Bus: Switching Status Schematic Diagram of PWM Modulator . . . . . . . . . . . . Block Diagram of the System Transfer Function . . . . . . . . Bode Diagram of the Open-Loop Transfer Function . . . . . . Resistance Load Simulation . . . . . . . . . . . . . . . . . . . Unit Step Response at Rsim = 2p.u. . . . . . . . . . . . . . . Voltage across Rsim = 2p.u. . . . . . . . . . . . . . . . . . . . Unit Step Response at Rsim = 1p.u. . . . . . . . . . . . . . . Voltage across Rsim = 1p.u. . . . . . . . . . . . . . . . . . . . Unit Step Response at Rsim = 20p.u. . . . . . . . . . . . . . . Relationship between Rsim and D . . . . . . . . . . . . . . . Resolution at Rsim = 2000pu . . . . . . . . . . . . . . . . . . Rsim = 0.5p.u. at fpwm = 20kHz . . . . . . . . . . . . . . . . Current Ripple at Rsim = 0.5p.u. . . . . . . . . . . . . . . . . Rsim = 1p.u. at fpwm = 20kHz . . . . . . . . . . . . . . . . . Current Ripple at Rsim = 1p.u. . . . . . . . . . . . . . . . . . Rsim = 2p.u. at fpwm = 20kHz . . . . . . . . . . . . . . . . . Current Ripple at Rsim = 2p.u. . . . . . . . . . . . . . . . . . Rsim = 2p.u. at fpwm = 20kHz . . . . . . . . . . . . . . . . . Current Ripple at Rsim = 2p.u. . . . . . . . . . . . . . . . . . Rsim = 20p.u. at fpwm = 20kHz . . . . . . . . . . . . . . . . Current Ripple at Rsim = 20p.u. . . . . . . . . . . . . . . . . fpwm = 10kHz, Rsim = 2p.u. and L = 2.6mH . . . . . . . . Current Ripple at Rsim = 2p.u. and fpwm = 10kHz . . . . . . Control Loop for Inductance Loads Simulations . . . . . . . . Simulated Inductor Voltages at Lsim = 0.5p.u. and 1p.u. . . 6 7 8 11 12 13 19 22 23 23 24 24 26 27 29 29 29 29 30 30 30 30 31 31 31 31 32 33 vii List of Figures 2.29 2.30 2.31 2.32 2.33 2.34 Simulated Inductor Voltages at Lsim = 1p.u. and 10p.u. . . Simulated Inductor Voltages at Lsim = 10p.u. and 100p.u. . Current Ripple at Lsim = 1p.u. . . . . . . . . . . . . . . . . . Current Ripple at Lsim = 10p.u. . . . . . . . . . . . . . . . . Simulated Capacitor Voltages at Csim = 1p.u. and 10p.u. . Simulated Capacitor Voltages at Csim = 10p.u. and 100p.u. 34 34 34 34 35 35 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 3.9 3.10 3.11 3.12 3.13 3.14 Hardware System of Electronic Loads . . . . . . . . . . . . . Structure of IGBT Module . . . . . . . . . . . . . . . . . . . Structure of IGBT Driver Module: 6SD106EI . . . . . . . . . HalfBridge Mode Selection . . . . . . . . . . . . . . . . . . . . Protection Circuit of Driver Module 6SD106EI . . . . . . . . Voltage Sensor Connection . . . . . . . . . . . . . . . . . . . . Voltage Measuring Circuit . . . . . . . . . . . . . . . . . . . . Voltage Scaling from Primary current to DSP Analog Voltage Structure of DSP On-chip ADC . . . . . . . . . . . . . . . . . Current Measuring Circuit . . . . . . . . . . . . . . . . . . . . Current Scaling from Primary Current to DSP Analog Voltage Temperature Sensor Circuit . . . . . . . . . . . . . . . . . . . Temperature Hysteresis Loop . . . . . . . . . . . . . . . . . . System Protection Circuit . . . . . . . . . . . . . . . . . . . . 36 37 38 39 40 41 42 43 43 44 45 45 46 47 4.1 4.2 4.3 4.4 4.5 4.6 4.7 Program Flow Chart . . . . . . . . . . . . . . . . ADC Conversion Time . . . . . . . . . . . . . . . Timer 1 and 2 ISR Sequence . . . . . . . . . . . Flowchart of P DP IN T A ISR . . . . . . . . . . . Sine Look-up Table . . . . . . . . . . . . . . . . . Phase Detection Diagram . . . . . . . . . . . . . The Flow Chart of the Current Calculation Block . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50 53 54 55 56 59 60 5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8 5.9 5.10 Control Circuit Board . . . . . . Power Circuit Board . . . . . . . The Whole Circuits . . . . . . . Experimental System . . . . . . . Rsim = 0.5p.u. at fpwm = 20kHz Current Ripple at Rsim = 0.5p.u. Rsim = 1p.u. at fpwm = 20kHz . Current Ripple at Rsim = 1p.u. . Rsim = 2p.u. at fpwm = 20kHz . Current Ripple at Rsim = 2p.u. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 61 61 61 61 63 63 63 63 63 63 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . viii List of Figures 5.11 5.12 5.13 5.14 5.15 5.16 5.17 5.18 5.19 5.20 5.21 5.22 5.23 5.24 5.25 5.26 5.27 5.28 5.29 5.30 Rsim = 2p.u. at fpwm = 10kHz and L = 2.6mH . . Current Ripple at Rsim = 2p.u. and fpwm = 10kHz Rsim = 2p.u. at fpwm = 20kHz and L = 2.6mH . . Current Ripple at Rsim = 2p.u. and fpwm = 20kHz L = 2.6mH, Rsim = 2p.u. and fpwm = 20kHz . . . Current Ripple at L = 2.6mH . . . . . . . . . . . . L = 26mH, Rsim = 2p.u. and fpwm = 20kHz . . . Current Ripple at L = 26mH . . . . . . . . . . . . Lsim = 0.5p.u. at fpwm = 20kHz . . . . . . . . . . Lsim = 1p.u. at fpwm = 20kHz . . . . . . . . . . . L = 2.6mH ,Lsim = 1p.u. and fpwm = 20kHz . . . L = 26mH ,Lsim = 1p.u. and fpwm = 20kHz . . . Current Ripple at L = 2.6mH . . . . . . . . . . . . Current Ripple at L = 26mH . . . . . . . . . . . . Csim = 0.5p.u. at fpwm = 20kHz . . . . . . . . . . Csim = 1p.u. at fpwm = 20kHz . . . . . . . . . . . L = 2.6mH ,Csim = 1p.u. and fpwm = 20kHz . . . L = 26mH ,Csim = 1p.u. and fpwm = 20kHz . . . Current Ripple at L = 2.6mH . . . . . . . . . . . . Current Ripple at L = 26mH . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64 64 65 65 66 66 66 66 67 67 68 68 68 68 69 69 70 70 70 70 6.1 6.2 Simplified Electronic Load Diagram . . . . . . . . . . . . . . Combination Electronic Load Diagram . . . . . . . . . . . . . 72 73 A.1 PowerSupply AD.SCHDOC . . . . A.2 VoltageSensor.SCHDOC . . . . . . A.3 CurrentSensor.SCHDOC . . . . . . A.4 TMPSensor.SCHDOC . . . . . . . A.5 Driver.SCHDOC . . . . . . . . . . A.6 MUBW20.SCHDOC . . . . . . . . A.7 Control Board PCB . . . . . . . . A.8 Control Board PCB Top Layer . . A.9 Control Board PCB Bottom Layer A.10 Power Board PCB . . . . . . . . . A.11 Power Board PCB Top Layer . . . A.12 Power Board PCB Bottom Layer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77 78 79 80 81 82 83 84 85 86 87 88 B.1 Schematics for 3-phase Islanding Test . . . . . . . . . . . . . 90 ix Acknowledgements I wish to express my deepest gratitude to my supervisor, Dr. William G. Dunford. During the thesis research, Dr. William G. Dunford provided enormous help and encouragement. This work would not be possible without his broad knowledge and patient instruction. I would like to thank Dr. Juri Jatskevich for his insights and valuable suggestions toward the completeness of this thesis. I would like to express my thanks to Michel AlSharidah who worked with me together in this research and shared with me his experience and optimism in the islanding tests and simulations. Thanks are extended to my colleges Tom De Rybel, Dean Chen, Yong Zhang and Mehmet Sucu for their great suggestions. Finally, I am expressing my sincerest gratitude to my parents and my girlfriend for their love and support during my studies. 1 Chapter 1 Introduction 1.1 Project Background Distributed Generation (DG) which is gaining more and more interest around the world generally refers to small-scale (typically 1kW – 50M W ) electric power generators that produce electricity at a site close to customers or that are tied to an electric distribution system. Among DG technologies, some are tried and tested but a number are new and developing technologies which use renewable energy sources, such as wind energy, biofuels and photovoltaic (PV) energy. Compared with the traditional larger power stations, DG promises significant reductions in the power transportation cost and reduces the capital and operating costs which will encourage their further growth. The integration of Distributed Generation into the main electricity networks doesn’t only bring in numerous benefits but some security issues. The anti-islanding protection[1] is the most frequent issue raised by DG interconnection to the utility networks because the island is an unregulated power system. Its behaviour is unpredictable due to the power mismatch between the load and generation and the lack of voltage and frequency control. Therefore, it is creating the possibility of damage to customer equipment in a situation over which the utility has no control. In wind and PV energy applications, the inverter is the most frequently used interface[2] which interconnects DG to the grids. According to IEEE standard 1547, the inverters used as interface in the distributed generation systems have to pass a series of tests and evaluation before interconnection to the utility. The unintentional islanding test is to verify that the anti-islanding protection of the DG interconnection component is capable of detecting the islanding conditions and ceasing to energize the electric power system. In the unintentional islanding test, a resistance, inductance and capacitance (RLC) load is employed as simulated local loads and shall be tuned at the nominal frequency of electricity networks. With increasing of the power rating of DG, the size of RLC loads used in the unintentional islanding test become bigger and bigger which makes 2 1.2. System Overview the tuning process of RLC circuits very time-consuming and inconvenient. Therefore, a compact and programmable electronic load is needed which can simulate different RLC networks. 1.2 1.2.1 System Overview Unintentional Islanding Test Scheme As described in the previous section, an interface inverter has to pass the unintentional islanding test before put into practical use. IEEE standard 1547 defines a generic system for unintentional islanding test which is illustrated in Figure 1.1. ∆P, ∆Q P, Q V, f P+∆P, Q+∆Q Figure 1.1: Generic System for Unintentional Islanding Test In the normal operation mode, switches S1, S2 and S3 are closed and the electricity network, local loads and DG interface inverter are interconnected together. The real and reactive power consumed in local loads are provided by both utility network and distributed generation system. Once the switch S1 is open, the electricity network is unable to provide power for customers and local loads are connected solely to DG system which is called ”islanding situation”. If the power generated by the DG system matches the RLC load power, i.e., in Figure 1.1, ∆P = ∆Q = 0. Under this condition, even if switch S1 is open, the voltage and frequency at the Point of Common Coupling (PCC) 3 1.3. Thesis Objective will not be changed and an unintentional islanding situation occurs which is considered the worst case for islanding detection. 1.2.2 AC Electronic Loads for Unintentional Islanding Test In this thesis, a single phase AC electronic load is designed to replace the RLC load shown in Figure 1.2 so that the tuning of RLC circuits can be achieved more easily by adjusting the values of electronic loads in the software. ∆P, ∆Q P, Q V, f P+∆P, Q+∆Q Figure 1.2: Electronic Loads for Unintentional Islanding Test The AC electronic loads consist of three parts: • A bidirectional H-Bridge based current source. • Voltage and current sensor circuits. • A TMS320LF2407A DSP control board. 1.3 Thesis Objective The objective of this thesis is to clarify the practical approaches needed to set up a programmable AC electronic load used in unintentional islanding test. The specific objectives include: • To find a suitable structure of the bidirectional current source 4 1.3. Thesis Objective • To find the design criteria of H-bridge based current source for the different ranges of the resistive, inductive and capacitive electronic loads with the 60Hz AC source. • To build a DSP-controlled H-bridge based current source. • To develop a control method for the programmable RLC electronic loads • To generate a DSP program for the TMS320LF2407A DSP. • To simulate different ranges of resistive, inductive and capacitive loads with a 60Hz sinusoidal voltage excitation and verify the theoretical analysis using the prototype. The thesis is organized into six chapters. Chapter 1 gives a brief introduction of the system and outlines the objectives of this thesis. In Chapter 2, the principles of the AC electronic loads are proposed and the PSIM simulations for different electronic loads are performed to give a set of specifications of the system. Chapter 3 is focused on hardware setup for the electronic loads. Chapter 4 will be dealing with software implementation. Some design illustration concerning software is presented. Selective experimental results are included in Chapter 5. The last chapter concludes the design and the implementation and proposes some work needed to be done in the future. 5 Chapter 2 Principles of the Simulated Loads In this chapter the principles of the simulated loads will be introduced step by step. First of all, the schematic of the simulated loads will be proposed carefully. Secondly, the ac small-signal model of the circuit will be derived as well as the stability of the system. Then, the relationship between the parameters of the circuit and the range of the simulated loads will be analyzed. At the end, a series of PSIM simulations for different kinds of loads will be performed and several groups of parameters will be given for different ranges of the simulated loads. 2.1 Schematic of Electronic loads Figure 2.1: Diagram of Electronic Loads Figure 2.1 shows the diagram of the electronic loads. It can be seen that the different load will draw the different current I from the input voltage source and thus the load characteristics can be changed if the output current of the source is under control. For an AC input source, the direction of the output current could be either positive or negtive, which means a bidirectional controllable current source is needed to control the output current of the source so that it can 6 2.1. Schematic of Electronic loads simulate different load characteristics. In this project, an H-bridge inverterbased bidirectional current source will be designed as an electronic load. 2.1.1 H-Bridge Bidirectional Current Source L g dc 2 R L 1 Figure 2.2: Controllable Current Source with Regulated DC Bus Voltage The schematic of the H-bridge inverter-based bidirectional current source is illustrated in Figure 2.2. It is very similar to the traditional H-bridge dcdc converter but the difference is that the port “AB” is used as the input port instead of the output port. The parameters of the circuit is explained as follows: • DC bus voltage. In order to make this circuit work properly, the DC bus voltage is regulated by a DC power supply which must be positive so that there are no short-circuits caused by the flywheel diodes which may lead to the serious damage of IGBTs. • Resistor R1 . In the real system, a series resistor R1 always exists because of the equivalent series resistance of the inductor L and the resistance of the conductors. • Resistor R2 . The resistor R2 is a dumping resistor and used to provide a closed loop for the inductor current IL . • Inductor L. the series inductor L is capable of smoothing the current IL so that this circuit can work as a current source. 7 2.1. Schematic of Electronic loads 2.1.2 The Range of Inductor Current In the following parts, the range of the current IL will be derived in order to analyze the electronic load characteristics. In simple terms of, a DC voltage source Vg will be applied in the circuit as an input source. However, the results are also available for low-frequency AC voltage sources, such as 60Hz. As depicted in Figure 2.2, the duty ratio D is the fraction of time that the switches S1 and S4 spend in the conducting position, and is a number between zero and one. The complement of the duty ratio, D , is defined as (1 − D), which means the fraction of time that the switches S2 and S3 spend in the conducting status. The small ripple approximation and the principles of inductor voltsecond balance[3] will be applied here to find the steady state of the inductor current IL . Figure 2.3 illustrates the operation topologies of the H-bridge inverter with regulated DC bus when the switches commutate. R L L DC 2 1 R g L L DC 2 1 g Figure 2.3: H-Bridge Inverter With Regulated DC Bus: (a) the first subinterval: D, (b) the second subinterval: D As illustrated in the Figure 2.3(a), during the first subinterval, DT , the switches S1 and S4 are on and the other switches S2 and S3 are off. The inductor voltage for this subinterval is given by vL = Vg − Vdc − iL × R1 (2.1) Use of the small ripple approximation, iL ≈ IL , leads to vL = Vg − Vdc − IL × R1 (2.2) 8 2.2. AC Small-Signal Modeling and PI Control Method Similarly, during the second subinterval, D T , the status of the switches are illustrated in the Figure 2.3(b). The inductor voltage for this subinterval is vL = Vg + Vdc − iL × R1 (2.3) Use of the small ripple approximation again leads to vL = Vg + Vdc − IL × R1 (2.4) According to the principles of inductor volt-second balance, the total volt-seconds applied to the inductor over one switching period is < vL >= 0 (2.5) Substitution of Eq.(2.2) and Eq.(2.4) into Eq.(2.5), and the result is D (Vg − Vdc − IL R1 ) + D (Vg + Vdc − IL R1 ) = 0 (2.6) Upon collecting terms, one obtains Vg = (2D − 1) Vdc + IL R1 (2.7) From Eq.(2.7), the average inductor current IL can be obtained IL = Vg + (1 − 2D) Vdc R1 (2.8) Therefore, the controllable range of the current flowing through the inductor L can be derived easily according to Eq.(2.8). When the switches S2 and S3 are always on, the inductor current IL reaches to the maximum value; when the switches S1 and S4 are always off instead, the minimum value can be obtained. The results are given by IL max = (Vg + Vdc )/R1 IL min = (Vg − Vdc )/R1 2.2 (2.9) AC Small-Signal Modeling and PI Control Method In this section, the ac small-signal modeling technique is employed for the derivation of the system transfer function and a conventional PI controller is designed accordingly for the PWM current controller. 9 2.2. AC Small-Signal Modeling and PI Control Method 2.2.1 AC Small-Signal Model The ac small-signal model is always used in the analysis of the dynamic performance of power electronic circuits. With this model, the transfer function of the system can be derived which is very useful for the stability analysis of the control system. However, this model only works at the lowfrequency and the disturbance of the switching frequency will be neglected. According to Figure 2.3, the voltage equations during these two intervals can be derived respectively. d : vg = vdc + iL (t) Ts · R1 + L d iL (t) dt Ts (2.10) d : vg = −vdc + iL (t) T s · R1 + L d iL (t) dt Ts Where iL (t) T s denotes the average of iL (t) over an interval of length Ts which represents one PWM switching period. Using Eq.(2.10), the average of voltage vg over one switching period equals: d + d vg = d − d vdc + d + d iL (t) Ts · R1 + L d iL (t) dt Ts (2.11) Upon collecting items, it can be rewritten as: vg = (2d − 1) vdc + iL (t) Ts · R1 + L d iL (t) dt Ts (2.12) Perturbing the inductor current, duty-cycle, input voltage and DC bus ∧ ∧ ∧ ∧ voltage with small signals of iL , d, vg and vdc respectively; also, assuming the quiescent value of these variables are IL , D, Vg and Vdc then: ∧ iL (t) Ts = IL + i L ∧ d=D+d vg = vdc = ∧ V g + vg ∧ Vdc + vdc (2.13) (2.14) (2.15) (2.16) Insert Eq.(2.13) to (2.16) into Eq.(2.12), one obtains: 10 2.2. AC Small-Signal Modeling and PI Control Method ∧ ∧ ∧ V g + vg = 2 D + d − 1 Vdc + vdc ∧ d IL + iL ∧ + IL + iL R1 + L dt (2.17) Linearizing Eq.(2.17), the ac small signal equation can be derived: ∧ vg = (2D − ∧ ∧ ∧ d iL ∧ 1) vdc +2Vdc d +R1 iL +L (2.18) dt Taking the Laplace transform for Eq.(2.18), the inductor current to duty ratio transfer function can be obtained: ∧ iL Gid (s) = d 2.2.2 = ∧ ∧ ∧ −2Vdc R1 + sL (2.19) vdc ,vg =0 The Transfer Function of Control System In this system, a traditional PI compensator is chose for the current feedback loop due to the uncompensated system containing a single pole. The transfer function of the compensator can be expressed as: ki (2.20) s The function of the pulse-width modulator is to produce the duty ratio d that is proportional to the output voltage of the compensator. Figure 2.4 shows the schematic diagram of the pulse-width modulator. Gc (s) = kp + saw c M s s s Figure 2.4: Schematic Diagram of PWM Modulator 11 2.2. AC Small-Signal Modeling and PI Control Method In terms of the waveform above, the duty cycle d will be a linear function of Vc . Hence, we can write vc (t) (2.21) VM Perturbation and linearization of Eq.(2.21), the transfer function of the pulse-width modulator is: d (t) = ∧ Gpwm (s) = d (2.22) ∧ vc Based on the analysis above, the system transfer function can be derived as follows (Figure 2.5). Note that the pwm modulator outputs d in our design, and thus a negative transfer function of the circuit is applied. In addition, a saturation block is also added to the PI compensator to eliminate the effect of integral saturation. err c p i M id ref sim Figure 2.5: Block Diagram of the System Transfer Function 2.2.3 Control System Stability Analysis In this project, a set of PI parameters of the control system is obtained using PSIM simulations and practical experiments: • PI parameters: kp = 1 and ki = 5 × 105 . • Saturation block: -20 to 20. • PWM modulator: VM = 40. According to the analysis above, the stability of these PI parameters can be investigated using Simulink7.0 and LTIViewer in MATLAB2008a. The bode diagram of the open-loop transfer function is depicted in Figure 2.6. It can be seen that the phase margin is always larger than zero so the system is stable. 12 2.3. System Performance Analysis From: In1 (pt. 1) To: Gain1 (pt. 1) 100 Magnitude (dB) 50 0 Phase (deg) -50 -90 -135 -180 10 2 10 3 10 4 10 5 10 6 10 7 Frequency (rad/sec) Figure 2.6: Bode Diagram of the Open-Loop Transfer Function 2.3 System Performance Analysis The electronic load designed in this project works in the switching mode, so the triangle waveforms are used to track the given sinusoidal waveforms. In order to meet the accuracy requirements, the effects of circuit parameters on the switching ripples and the system response time have to be investigated carefully. Therefore, this section will explore the relationship between the circuit parameters and the system performance. 2.3.1 Current Ripple Figure 2.2 shows that the inductor L which prevents the current from being changed suddenly plays the most important role in the current control. In this section, the effect of circuit parameters on the current ripple will be investigated carefully and the formula of the current ripple will be derived for the calculation of circuit parameters according to the given accuracy demands . The current ripple used here is defined as ∆iL /2IL . Where ∆iL is the peak-to-peak current change and IL is the average current over one switching period. The following parts will introduce the relationship between the current ripple and the circuit parameters for the simulated resistive, inductive 13 2.3. System Performance Analysis and capacitive loads separately. Current Ripple for Resistive Loads In simple terms, a DC voltage source Vg is employed as the input voltage source as shown in Figure 2.2 and the circuit is controlled to simulate the resistance characteristic. To the extent of an low-frequency AC voltage source, such as 60Hz, the analysis still works fine because the AC source can be discretized into a series of DC voltage source at low frequency. As illustrated in Figure 2.2, the current change ∆iL is related to the inductor L and the current ripple ∆iL /2IL should be as small as possible so that the high frequency noise of the electronic loads can be decreased greatly. Only consider the intervals D T that S2 and S3 are on. According to the inductance characteristic, one obtains L diL = vL dt (2.23) Where vL = Vdc + Vg − vR1 When R1 (2.24) Rsim (Rsim : simulated resistor), vR1 can be ignored vL ≈ Vdc + Vg (2.25) Insert Eq.(2.25) into Eq.(2.23), one obtains L diL ≈ Vdc + Vg = (1 + k) Vg dt (2.26) Where Vdc = Vg , then the current change ∆iL is equal to ∆iL = (1 + k) Vg × ∆t L (2.27) Using Eq.(2.8), the steady state of inductor current IL is given by IL = Vg +(1−2D)Vdc R1 = Vg +(1−2D)kVg R1 = 1+(1−2D)k Vg R1 (2.28) 14 2.3. System Performance Analysis Because the circuit is working as resistance loads, the average current of inductor IL is also equal to IL = Vg Rsim (2.29) The duty ratio D can be derived through combination of Eq.(2.28) and Eq.(2.29) Vg 1+(1−2D)k × Vg = Rsim R1 (2.30) R1 1 ⇒ D = 2k k + 1 − Rsim The time interval ∆t can be obtained using Eq.(2.30) ∆t = (1 − D) × Tpwm (2.31) = (k−1)Rsim +R1 2·k·Rsim × Tpwm Combining Eq.(2.27) and Eq.(2.29), the ratio of the current change ∆iL to inductor current IL is equal to (k + 1) Rsim ∆iL = × ∆t IL L (2.32) Insert Eq.(2.31) into Eq.(2.32), the current ripple is obtained ∆iL 2IL = 1 2 = 1 2 × (k+1)Rsim L × Tpwm L × (k−1)Rsim +R1 2·k·Rsim × Tpwm (2.33) × (k2 −1)Rsim +(k+1)R1 2k The equation above indicates that the current ripple ∆iL /2IL is a function of the simulated resistive load Rsim , the PWM period Tpwm and the value of the inductor L when k and series resistor R1 are constant. According to this function, the following relationship can be concluded: • The maximum impedance of the simulated loads which is achievable on this system can be derived by Eq.(2.33) when all the circuit parameters and the current ripple requirement are given. • The larger inductor L, the smaller the current ripple. Hence, the minimum value of the inductor L can be calculated with Eq.(2.33) when the simulated resistive load Rsim and the current ripple are given. • The higher switching frequency fpwm , the smoother the current ripple. 15 2.3. System Performance Analysis Current Ripple for Inductive and Capacitive Loads The expression of current ripple for the inductive and capacitive loads are exactly same, so here we only take the inductive loads for example to derive the current ripple expression. From the schematic of the current source, it can be concluded that the largest current ripple will occur when the AC input source reaches zero. Therefore, following the similar process shown in the previous section, one obtains L diL ≈ VL dt (2.34) Where vL = Vdc − vR1 (2.35) When R1 Zsim (Zsim : the impedance of the simulated inductive or capacitive loads), vR1 can be ignored vL ≈ Vdc (2.36) Insert Eq.(2.36) into Eq.(2.34), one obtains diL ≈ Vdc = kVg dt Then the current change ∆iL is equal to L (2.37) kVg × ∆t (2.38) L According to Eq.(2.8), when Vg = 0, the average current IL over one switching cycle is given by ∆iL = IL = 1 − 2D kVg R1 (2.39) When Vg = 0, the average current IL also equals IL = Vg Zsim (2.40) Combining Eq.(2.39) and (2.40), the duty cycle D is equal to D= 1 2 1− R1 k Zsim (2.41) 16 2.3. System Performance Analysis When R1 Zsim , the duty cycle D is approximately 0.5 and the time interval ∆t = 0.5Tpwm . Therefore, combining Eq.(2.38) and (2.40), the current ripple for inductive or capacitive loads is Tpwm ∆iL 1 k Zsim = × × 2IL 2 L 2 (2.42) The expression above shows the similar function as the resistive loads. However, compared with the function of the resistive loads, the current ripple of the simulated inductive or capacitive loads will be larger with the same circuit parameters. 2.3.2 Relationship between Response Time and Zsim To simulate electronic loads under low-frequency voltage source, the response time of the system should be small enough so that the system is capable of tracking the required load characteristic curves. The following derivations are based on the resistive simulated loads and the conclusion is also available for the inductive or capacitive simulated loads. According to Eq.(2.27) and (2.31), the current change function can be derived (1 + k) Vg (k − 1) Rsim + R1 ∆iL = × × Tpwm (2.43) L 2k × Rsim It can be derived that if the circuit parameters keep constant, the current change ∆iL will increase with the decreasing of Rsim . However, the series inductor L will prevent the current from changing suddenly and thus the response time will increase greatly when the simulated load Rsim decreases. Hence, for a given series inductor L and other given parameters, the minimum achievable impedance of the simulated loads can be found considering the proper response time. In the following part, the unit step response of the system will be investigated to find the proper minimum values for the given parameters previously. 2.3.3 Conclusions of System Performance The previous sections explain the relationship between the circuit parameters and the system performance. Therefore, a general design procedure can be used according to the requirements of the system performance: • First of all, determine the proper values of k and fpwm and the maximum impedance of the simulated loads. 17 2.4. Calculation of Circuit Parameters • Secondly, according to the current ripple requirement, calculate the minimum inductor L with Eq.(2.33) or Eq.(2.42). • At the end, find the minimum impedance of the simulated loads according to the given circuit parameters so that the response time of the system is small enough for the 60Hz AC source. 2.4 Calculation of Circuit Parameters The schematic and performance of the system have been explained previously. This section will introduce how to calculate the circuit parameters according to the given accuracy demands and a set of values of circuit parameters will be listed for the different ranges of the simulated loads. The analysis proposed here is all based on the PSIM simulations and the theoretical derivation. 2.4.1 Introduction of PSIM Simulation Diagram Powersim (PSIM) is a powerful circuit simulation software which is specially suited to power electronic circuit simulations. All kinds of semiconductor parts and the control blocks have been integrated in the software package, which makes it very easy to use. In order to analyze the circuit parameters, a PSIM simulation program for the simulated resistive loads will be used in this section and the conclusions obtained from this system is also suitable for the inductive and capacitive loads with some minor changes. Figure 2.7 shows the schematic of the resistive load simulation consisting of the H-Bridge current source and the control program. The control program senses the inductor current iL and the input voltage source Vg to calculate the error between the actual voltage and the demanded voltage using equation: error = Vsource − Isense × Rsim (2.44) Where Rsim represents the resistance of the simulated resistive loads. With a set of proper PI parameters (given in the previous section), the error Verr can be obtained which is used to compared with a saw-tooth waveform to generate PWM signals. In the PSIM simulation, the per-unit system is employed and the voltage and current base values are equal to 30V and 0.78125A separately which are also used in the DSP program and thus the impedance base value is 38.4Ω. The choice of the base values will be explained in the chapter four. 18 2.4. Calculation of Circuit Parameters Figure 2.7: Resistance Load Simulation 2.4.2 DC Bus Voltage: Vdc As discussed in the previous sections, the DC bus voltage Vdc must be positive so that there are no short-circuits and according to Eq.(2.9), Vdc should be greater than or equal to Vg to increase the range of the inductor current, in other words, to increase the range of the electronic loads. A parameter of k is defined as the ratio of Vdc to Vg and should be larger than or equal to 1. However, if k is 1, it’s going to take a very long period to decrease the current due to the zero reverse voltage across the series inductor; if k is much larger than 1, the current ripple will be increased significantly. Through the practical experiments, the parameter of k is kept to be 1.3 to improve the dynamic performance of the device as well as the current ripple. 2.4.3 Series Resistor: R1 From Eq.(2.9), IL max = (Vg + Vdc )/R1 , it can be seen that the minimum impedance of electronic loads is decided by the series resistor R1 . When 19 2.4. Calculation of Circuit Parameters Vdc = 1.3Vg , the minimum simulated impedance |Zsim |min = R1 /2.3. In our system, the series resistor R1 is 17Ω and the minimum limitation of the simulated impedance is about 7.4Ω. 2.4.4 Switching Frequency: fpwm As discussed previously, the current ripple depends upon the switching frequency fpwm and the higher switching frequency, the smaller current ripple. In this project, the switching frequency fpwm is set to 20kHz which is the most typical frequency of IGBTs. 2.4.5 Inductor: L It has been discussed in the previous sections that the value of inductor L is determined by the requirement of the current ripple (Eq.(2.33) and Eq.(2.42)). On the other hand, once the inductor value and other parameters are chosen, the range of the impedance of the simulated loads can be calculated with the consideration of both the current ripple and the response time. The following parts will show how to calculate the inductor value and the range of the simulated loads. Calculation of Inductor L For resistive loads, three different ranges of the simulated loads are given in Table 2.1. From Eq.(2.33), a set of minimum values of the series inductor L can be calculated for the resistive loads according to the given current ripple and the maximum impedance of the simulated loads. Assume the current ripple ∆iL /2IL should be less than 0.3, the minimum inductors for the simulated resistive loads are shown as follows. Range of Rsim 0 − 2p.u. 0 − 20p.u. 0 − 200p.u. Lmin 2.6mH 26mH 260mH Maximum Current Ripple 0.34 (actual value: 0.3) 0.21 0.19 Table 2.1: Series Inductors for Different Ranges of Resistance Loads: switching frequency fpwm = 20kHz, voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω 20 2.4. Calculation of Circuit Parameters It can be found that the current ripple is larger than 0.3 at Rsim = 2p.u. but the real current ripple is equal to 0.3 because the equation we used here is derived with the approximation of R1 Rsim but Rsim = 2p.u. is comparable with R1 . Therefore, the actual current ripple at small simulated resistance is less than the calculated value (the real current ripple can be found using the PSIM simulations in Figure 2.20). The similar results for the inductive or capacitive simulated loads can be obtained with Eq.(2.42). Considering the hardware consistency, the same set of values of the series inductor L as shown in Table 2.1 will be used in the calculations, so the maximum impedance of the simulated loads and the associated current ripple are calculated as follows. Maximum Zsim 1.25p.u. 12.5p.u. 125p.u. Series Inductor L 2.6mH 26mH 260mH Maximum Current Ripple 0.30 0.30 0.30 Table 2.2: Series Inductors for Different Ranges of Inductive or Capacitive Loads: switching frequency fpwm = 20kHz, voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω Practically, with the consideration of simplicity, another set of values of Zsim is given in Table 2.3 after the rounding adjustment. The parameters listed above are all based on the assumption that the current ripple is less than 0.3. If a smaller current ripple is required for the specific application, the simplest approach is just to increase the switching frequency fpwm proportionally and retain other parameters. Maximum Zsim 1p.u. 10p.u. 100p.u. Series Inductor L 2.6mH 26mH 260mH Maximum Current Ripple 0.24 0.24 0.24 Table 2.3: Series Inductors for Different Ranges of Inductive or Capacitive Loads: switching frequency fpwm = 20kHz, voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω 21 2.4. Calculation of Circuit Parameters 2.4.6 Calculation of the Range of Zsim As discussed above, once the current ripple requirement and the series inductor are given, the maximum achievable impedance of the simulated loads is also set up. However, not all impedances less than the maximum value are able to be simulated with this system. In this section, the response time of the system will be analyzed and the minimum achievable impedance of the simulated loads can be determined by the response time. Simulation Results of the Response Time Here we takes L = 26mH for example to find the minimum limitation of the impedance. In the PSIM simulations shown in Figure 2.7, the parameters are: fpwm = 20kHz, Rsim = 2p.u.(76.8Ω), R1 = 17Ω and L = 26mH. To explore the unit step response, a unit step input voltage source Vg is applied. Figure 2.8 depicts the unit step response of the system at Rsim = 2p.u.. Replacing the unit step source with a 60Hz AC source, the voltage across Rsim (2p.u.) is plotted in Figure 2.9. It can seen that the rise time is equal to 0.46ms and the settling time is 3.50ms which is fast enough for the application of 60Hz. V_Rsim (2pu) 2.00 Vg rise 1.50 1.00 0.50 settle 0.0 200.00 202.00 204.00 206.00 208.00 210.00 212.00 Time (ms) Figure 2.8: Unit Step Response at Rsim = 2p.u. However, if the resistive simulated load Rsim keeps going down to 1p.u., 22 2.4. Calculation of Circuit Parameters V_Rsim Vg 1.00 1.0 0.50 0.5 0.0 0.0 -0.50 -0.5 -1.00 -1.0 -1.50 -1.5 1.50 1.5 1.00 1.0 0.50 0.5 0.0 0.0 -0.50 -0.5 -1.00 -1.0 -1.50 -1.5 V_Rsim (2pu) V_Rsim Vg Vg 1.0 1.00 0.5 0.50 0.0 0.0 -0.50 -0.5 -1.00 -1.0 -1.50 -1.5 120.00 130.00 140.00 150.00 160.00 Time (ms) Figure 2.9: Voltage across Rsim = 2p.u. the unit step response time will increase consistently as shown in Figure 2.10. Replacing the unit step input source with a 60Hz AC source, the voltage waveform across Rsim is depicted in Figure 2.11. It can be seen that the average value of VRsim is no longer equal to the input voltage due to the long response time. V_Rsim (1pu) Vg rise 1.50 10% error 1.00 0.50 0.0 settle 200.00 202.50 205.00 207.50 210.00 212.50 Time (ms) Figure 2.10: Unit Step Response at Rsim = 1p.u. 23 2.4. Calculation of Circuit Parameters V_Rsim (1pu) Vg V_Rsim (1pu) Vg 1020.00m 1.00 1010.00m Zoom In 1000.00m 0.50 990.00m 170.75 170.80 170.85 Time (ms) 170.90 0.0 -0.50 -1.00 150.00 160.00 170.00 180.00 Time (ms) 190.00 Figure 2.11: Voltage across Rsim = 1p.u. On the other hand, if Rsim increases, the response time will keep decreasing and the response time for Rsim = 20p.u. is illustrated in Figure 2.12. Compared with the response time at Rsim = 2p.u., the response time at Rsim = 20p.u. is much faster but more ripples. V_Rsim Vg 3.00 rise 48% error 2.00 1.00 0.0 settle -1.00 100.00 102.00 104.00 106.00 Time (ms) Figure 2.12: Unit Step Response at Rsim = 20p.u. 24 2.4. Calculation of Circuit Parameters The Range of Zsim Based upon the discussion above it can be concluded that the minimum achievable impedance of the simulated loads is determined by the response time if the value of the series inductor is given. With the help of PSIM simulations, a set of values of the minimum Rsim is found in Table 2.4. Series Inductor L 2.6mH 26mH 260mH Minimum Rsim ≥ 0.2p.u. ≥ 2pu ≥ 20pu Table 2.4: Minimum Rsim for the Given Series Inductor: switching frequency fpwm = 20kHz, voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω Similarly, the minimum impedance of the inductive or capacitive simulated loads is also able to be found in the table below with the PSIM simulation programs. Series Inductor L 2.6mH 26mH 260mH Minimum Lsim (or Csim ) ≥ 0.1p.u. ≥ 1pu ≥ 10pu Table 2.5: Minimum Lsim (or Csim ) for the Given Series Inductor: switching frequency fpwm = 20kHz, voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω 2.4.7 Duty Ratio D This section will introduce the effect of the duty cycle D on the range and resolution of the simulated loads. According to Eq.(2.28) and Eq.(2.29), a resistive simulated load can be represented as: Rsim = Vg R1 = IL 1 + (1 − 2D) · k (2.45) In a continuous system, the value of Rsim could be infinite when the denominator is approaching zero. However, the real system we used in this 25 2.4. Calculation of Circuit Parameters project is a digital system with a 16-bit DSP so the step size of the duty ratio D is 1 216 . With this step size, the range of the resistive load is illustrated in Figure 2.13 for the given parameters above. 6 3 x 10 Relationship between D and R sim 2.5 2 Rsim 1.5 1 0.5 0 -0.5 -1 0 0.2 0.4 0.6 Duty Ratio D 0.8 1 Figure 2.13: Relationship between Rsim and D It can be found that the slope of the function Rsim (D) will increase significantly when the duty cycle D approaches the value which makes the denominator equal zero (Here it is 0.8846). This phenomenon indicates that if the step size of D is fixed, the resolution of the simulated resistive load is going down when D approaches that sensitive value. In order to achieve a higher resolution of the simulated resistive load, the simulated impedance should fall into a reasonable range depending on the required resolution. Using Eq.(2.45), the resolution can be calculated for the simulated resistive loads shown in the previous section (See Table 2.6). Range of Rsim < 2p.u. < 20p.u. < 200p.u. Resolution (%) ≤ 0.018 ≤ 0.18 ≤ 1.8 Table 2.6: The Resolution of Rsim : switching frequency fpwm = 20kHz, voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω 26 2.4. Calculation of Circuit Parameters For the impedance larger than 200p.u., the resolution becomes too small to be applied in the practical system. Figure 2.14 depicts that the resolution of Rsim at 2000p.u. is going up to 18.7% and this error is too big for a practical application. x 10 5 Resolution of Rsim (%) 4 Rsim (Ohms) 3 2 X: 0.8845 Y: 8.004e+004 1 X: 0.8845 Y: 6.744e+004 0 Resolution 18.7% at R = 2000pu sim -1 -2 0.8841 0.8842 0.8843 0.8844 Duty Cycle D 0.8845 Figure 2.14: Resolution at Rsim = 2000pu According to Eq.(2.41), the similar results can be derived for the inductive or capacitive simulated loads (as shown in Table 2.7) Range of Lsim (or Csim ) < 1p.u. < 10p.u. < 100p.u. Resolution (%) ≤ 0.0089 ≤ 0.090 ≤ 0.90 Table 2.7: The Resolution of Lsim (or Csim ): switching frequency fpwm = 20kHz, voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω 2.4.8 Specifications of Electronic Loads From the previous discussion, a series values of the circuit parameters can be derived for the resistive, inductive and capacitive simulated loads shown in the Table 2.8 and 2.9 respectively. 27 2.5. Simulations of Electronic Loads Rsim 0.5 − 2p.u. 2 − 20p.u. 20 − 200p.u. Inductor L 2.6mH 26mH 260mH Resolution (%) ≤ 0.018 ≤ 0.18 ≤ 1.8 Current Ripple ≤ 0.3 < 0.3 < 0.3 Table 2.8: Specifications for the Resistive Simulated Loads: switching frequency fpwm = 20kHz, , voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω Lsim (Csim ) 0.5 − 1p.u. 1 − 10p.u. 10 − 100p.u. Inductor L 2.6mH 26mH 260mH Resolution (%) ≤ 0.00089 ≤ 0.090 ≤ 0.90 Current Ripple ≤ 0.24 ≤ 0.24 ≤ 0.24 Table 2.9: Specifications for the Inductive or Capacitive Simulated Loads: switching frequency fpwm = 20kHz, , voltage ratio k = 1.3, series resistor R1 = 17Ω and the impedance base value Z = 38.4Ω These parameters are all obtained based on the 0.3 current ripple assumption and increasing the switching frequency is able to achieve the smaller current ripple. 2.5 Simulations of Electronic Loads In this section, RLC electronic loads will be simulated using PSIM8.0 separately to verify the principles derived in the previous sections and the calculations of the specifications suggested above. 2.5.1 Resistive Load Simulation The PSIM simulation program for the resistive simulated loads has been depicted in Figure 2.7. Based upon the previous introduction and the parameters, the simulations can be performed and the waveforms are shown below for the current ripple analysis. Current Ripple and Inductor L According to Eq.(2.33), the current ripple can be calculated when R1 Rsim . But when the resistance of the simulated load is comparable to R1 , the 28 2.5. Simulations of Electronic Loads actual ripples are smaller than the calculated values obtained with Eq.(2.33). For example, assume the switching frequency fpwm = 20kHz, voltage ratio k = 1.3, series resistor R1 = 17Ω and series inductor L = 2.6mH. The calculated ripple is equal to 0.19 at Rsim = 0.5p.u., 0.24 at Rsim = 1p.u. and 0.34 at Rsim = 2p.u.. However, in the PSIM simulations, the actual ripples can be obtained (See Figure 2.15 to 2.20) and equal 0.13, 0.19 and 0.3 separately which are smaller than the calculated values. V_Rsim (0.5pu) V_Rsim (0.5pu) Vg Vg 1.10 V_Rsim Zoom In 1.00 ΔiL/2IL=0.13 1.05 Vg 0.50 1.00 0.0 -0.50 0.95 -1.00 0.90 154.10154.12154.14154.16154.18154.20154.22154.24 Time (ms) 130.00 140.00 150.00 160.00 170.00 Time (ms) Figure 2.15: Rsim = 0.5p.u. at fpwm = 20kHz Figure 2.16: Current Ripple at Rsim = 0.5p.u. V_Rsim (1pu) V_Rsim (1pu) V_Rsim 1.00 Vg Vg 1.20 Zoom In 1.10 Vg ΔiL/2IL=0.19 0.50 1.00 0.0 -0.50 0.90 -1.00 0.80 140.00 160.00 Time (ms) 180.00 Figure 2.17: Rsim = 1p.u. at fpwm = 20kHz 154.02154.04154.06154.08154.10154.12154.14 Time (ms) Figure 2.18: Current Ripple at Rsim = 1p.u. 29 2.5. Simulations of Electronic Loads V_Rsim (2pu) V_Rsim (2pu) 1.50 Vg V_Rsim 1.00 Vg 1.30 1.20 Zoom In Vg ΔiL/2IL=0.30 1.10 0.50 0.0 1.00 -0.50 0.90 -1.00 0.80 -1.50 0.70 140.00 160.00 Time (ms) 154.016 154.048 180.00 Figure 2.19: Rsim = 2p.u. at fpwm = 20kHz 154.08 154.112 154.144 Time (ms) Figure 2.20: Current Ripple at Rsim = 2p.u. To investigate the current ripples for 2p.u. ≤ Rsim ≤ 20p.u., the 2.6mH inductor is replaced with a 26mH inductor to simulate Rsim = 2p.u. and 20p.u. respectively. Using Eq.(2.33) again, the current ripple equals 0.21 at Rsim = 20p.u.. The waveforms of voltage across Rsim and the associated current ripples are depicted in Figure 2.21 to 2.24. Figure 2.24 shows that the current ripple at Rsim = 20p.u. equals 0.20 which is very close to the calculated value 0.21 because Rsim = 20p.u. = 768Ω is much bigger than R1 = 17Ω. In Figure 2.22,the current ripple is 0.03 at Rsim = 2p.u. and L = 26mH. Compared with the value when L is 2.6mH, the current ripple goes down to 10% which is inversely proportional to the series inductor L same as the relationship shown in Eq.(2.33). V_Rsim (2pu) 1.00 V_Rsim Vg 0.50 V_Rsim (2pu) Vg Vg Zoom In 1.02 1020.00m ΔiL/2IL=0.03 0.0 1.00 1000.00m -0.50 980.00m 0.98 -1.00 130.00 140.00 150.00 160.00 170.00 Time (ms) Figure 2.21: Rsim = 2p.u. at fpwm = 20kHz 154.105 Figure 2.22: Rsim = 2p.u. 154.142 154.179 Time (ms) Current Ripple at 30 2.5. Simulations of Electronic Loads V_Rsim (20pu) Vg V_Rsim (20pu) 1.00 Zoom In V_Rsim Vg 1.20 Vg 0.50 1.10 0.0 1.00 -0.50 ΔiL/2IL=0.20 0.90 -1.00 0.80 -1.50 130.00 140.00 150.00 160.00 170.00 Time (ms) Figure 2.23: Rsim = 20p.u. at fpwm = 20kHz 154.123 154.154 154.185 154.216 154.246 Time (ms) Figure 2.24: Current Ripple at Rsim = 20p.u. Current Ripple and Switching Frequency fpwm Using Eq.(2.33), the current ripple is inversely proportional to the switching frequency fpwm and the PSIM simulation results are shown in following figures. It can be seen in Figure 2.26 that the current ripple ∆iL /2IL equals to 0.61 and is 2 times of the current ripple at 20kHz which proves the inversely proportional relationship between the current ripple and the switching frequency. V_Rsim (2pu) at 10kHz Vg 1.60 Zoom In V_Rsim 1.00 V_Rsim (2pu) at 10kHz Vg Vg ΔiL/2IL=0.61 1.20 0.0 0.80 -1.00 0.40 130.00 140.00 150.00 Time (ms) 160.00 170.00 Figure 2.25: fpwm = 10kHz, Rsim = 2p.u. and L = 2.6mH 153.92 154.00 154.08 Time (ms) 154.16 Figure 2.26: Current Ripple at Rsim = 2p.u. and fpwm = 10kHz 31 2.5. Simulations of Electronic Loads 2.5.2 Inductive Load Simulation Voltage Magnitude Phase Angle Figure 2.27: Control Loop for Inductance Loads Simulations For the inductance load simulations, no changes have to be done in the circuit but the control algorithm which is shown in Figure 2.27. The control program senses the inductor current iL and input voltage Vg at first. Then, the magnitude and phase angle of Vg are detected and with 90 degrees phase delay, a new voltage magnitude is calculated which can be compared with the inductor current iL directly. Waveforms for Different Ranges of the Inductive Loads With the control program mentioned above, the inductive loads simulations can be achieved using the parameters given previously. The parameters used in the simulations are: Vg = 10V , fpwm = 20kHz, Lsim = 0.5p.u. and 1p.u., R1 = 17Ω and L = 2.6mH. The waveforms of simulated inductor voltages are depicted in Figure 2.28. It can be found in Figure 2.28 that the AC input voltage Vg is leading the current flowing through the series inductor and resistor by 90 degrees which 32 2.5. Simulations of Electronic Loads V_LsimVg(1pu) Vg V_Lsim (1pu) 15.0 15.00 10.00 10.0 5.00 5.0 0.0 0.0 -5.00 -5.0 -10.00 -10.0 -15.00 -15.0 V_cmd V_sensedV_cmd (1pu) V_sensed (1pu) 0.40 0.4 0.20 0.2 0.00.0 -0.20 -0.2 -0.40 -0.4 100.00 100.00 110.00 110.00 V_Lsim Vg(0.5pu) V_Lsim (0.5pu) 120.00 120.00 Time (ms) 130.00 130.00 140.00 140.00 150.00 150.00 Vg 10.00 10.0 5.00 5.0 0.0 0.0 -5.00 -5.0 -10.00 -10.0 V_sensedV_cmd (0.5pu) V_sensed (0.5pu) V_cmd 0.40 0.4 0.20 0.2 0.00.0 -0.20 -0.2 -0.40 -0.4 100.00 100.00 110.00 110.00 120.00 130.00 120.00 130.00 Time (ms) Time (ms) 140.00 140.00 150.00 150.00 Figure 2.28: Simulated Inductor Voltages: Vsensed is the simulated inductor voltages in per-unit system and Vcmd is the AC input voltage source in perunit system with 90 degrees phase delay represents the inductive load characteristics. The waveforms for other ranges of the inductive simulated loads are illustrated in Figure 2.29 and 2.30. Current Ripple Waveforms From the previous discussion, when Zsim ≥ R1 , the current ripple can be calculated using Eq.(2.42); when Zsim is comparable with the series resistor R1 , the actual current ripple should be smaller than the calculated value. Figure 2.31 and 2.32 show the current ripples of the inductive simulated loads when L = 26mH and fpwm = 20kHz. As shown in these figures, the current ripple at Lsim = 10p.u. is 0.23 33 2.5. Simulations of Electronic Loads V_Lsim (10pu) V_Lsim (1pu) 10.00 10.0 Vg 10.00 10.0 Vg V_Lsim Vg Vg V_Lsim 5.00 5.0 5.00 5.0 0.0 0.0 0.0 0.0 -5.00 -5.0 -5.00 -5.0 -10.00 -10.0 -10.00 -10.0 100.00 100.0 110.00 110.0 120.00 120.0 Time (ms) 130.00 130.0 140.00 140.0 150.00 150.0 100.00 100.0 V_Lsim (100pu) V_Lsim (10pu) 10.00 10.0 110.00 110.0 120.00 120.0 130.00 130.0 Time (ms) 140.00 140.0 150.00 150.0 160.00 160.0 Vg Vg V_Lsim Vg 10.00 10.0 V_Lsim Vg 5.00 5.0 5.00 5.0 0.0 0.0 0.00.0 -5.00 -5.0 -5.00 -5.0 -10.00 -10.0 -10.00 -10.0 100.00 100.0 110.00 110.0 120.00 120.0 130.0 130.00 Time (ms) 140.00 140.0 150.00 150.0 100.00 100.0 160.00 160.0 110.00 110.0 Time (ms) Figure 2.29: Simulated Inductor Voltages at Lsim = 1p.u. and 10p.u. 335.00m 140.00 140.0 150.00 150.0 160.00 160.0 Figure 2.30: Simulated Inductor Voltages at Lsim = 10p.u. and 100p.u. V Vg cmd VVg cmd 0.40 340.00m 337.50m 130.00 130.0 Time (ms) Time (ms) V_Lsim (10pu) V_Lsim (1pu) 120.00 120.0 ΔiL/2IL=0.21 ΔiL/2IL=0.23 0.35 332.50m 0.30 330.00m 327.50m 0.25 325.00m 191.62191.64191.66191.68191.70191.72191.74 Time (ms) Figure 2.31: Current Ripple at Lsim = 1p.u. 208.32 208.36 208.40 Time (ms) 208.44 Figure 2.32: Current Ripple at Lsim = 10p.u. which is very close to the calculated value 0.24; the current ripple at Lsim = 1p.u. is 0.21 which is smaller than the calculated value 0.24 because 1p.u. is comparable with the series resistor R1 . 34 2.5. Simulations of Electronic Loads 2.5.3 Capacitive Load Simulation The control program of capacitance loads simulations is almost the same as the program used in the inductance loads simulation except that the phase angle of series inductor current iL leads that of input voltage Vg by 90 degrees instead of 90 degrees delay. Using the parameters given previously, Figures 2.33 and 2.34 show the waveforms of the different ranges of capacitive loads. V_Lsim V_Csim(1pu) (1pu) V_Csim (10pu) Vg Vg V_Csim Vg 10.00 10.00 Vg V_Csim Vg 5.00 5.00 0.0 0.0 -5.00 -5.00 -10.00 -10.00 110.00 120.00 130.00 140.00 Time (ms) V_Lsim V_Csim(10pu) (10pu) 150.00 130.00 140.00 150.00 160.00 170.00 Time (ms) 160.00 V_Csim (100pu) VgVg V_Csim V_Csim Vg 10.00 10.00 Vg 5.00 5.00 0.0 0.0 -5.00 -5.00 -10.00 -10.00 110.00 Vg 120.00 130.00 140.00 Time (ms) 150.00 160.00 Figure 2.33: Simulated Capacitor Voltages at Csim = 1p.u. and 10p.u. 130.00 140.00 150.00 160.00 170.00 Time (ms) Figure 2.34: Simulated Capacitor Voltages at Csim = 10p.u. and 100p.u. 35 Chapter 3 Hardware Implementation Based on the theory discussed in chapter 2, a practical AC electronic load system will be built for experimental results. In this chapter, emphasis will be given on how to choose suitable components and integrate them to form a prototype of AC electronic load. 3.1 Introduction of System The whole system is depicted in Figure 3.1. It consists of a regulated DC bus, an H-bridge inverter, a current filter, an IGBT driver module, sensor circuits, a DSP control board and a power supply circuit. 2 Voltage, Current and Temperature Signals Isolated IGBT Driver Sensor Circuits Power Supply DSP Control Board Figure 3.1: Hardware System of Electronic Loads The sensor circuits sense voltage, current and temperature signals. Voltage and current signals will be used to calculate current reference in the DSP 36 3.2. IGBT Inverter Module software, and temperature signal will be used for overheat protection. An on-board power supply with some accessorial components provides ±15V , +15V and +5V for the whole hardware system. In the following parts, each subsystem will be introduced in detail. 3.2 IGBT Inverter Module In this design, an IGBT module, MUBW 20-06 A7[4], produced by IXYS Corporation is employed to set up the H-bridge inverter. As depicted in Figure 3.2[4], the IGBT module includes a three phase full-bridge rectifier, MUBW 35-12 A7 a three phase DC-AC inverter and a NTC temperature sensor. For the application of single phase AC electronic loads, only the 3-phase inverter and NTC temperature sensor will be used. Converter - Brake - Inverter Module (CBI2) 21 D11 D13 22 D7 D15 7 1 2 3 D12 D14 D16 T1 D1 16 15 6 T7 14 23 T2 11 10 T3 18 17 T4 D2 D3 T5 D5 20 19 5 D4 12 4 T6 D6 13 24 8 NTC 9 Three Phase Figure Brake Three 3.2:Chopper Structure of Phase IGBT Rectifier Inverter Module VRRM = 1600V VCES = 1200 V VCES = 1200 V IDAVM = 44 A IC25 = 35 A IC25 = 50 A According to the supposed conditions of unintentional islanding test IFSM = 400 A VCE(sat) = 2.3 V VCE(sat) = 2.6 V (Appendix C), a voltage rating above 120VRM S and a current rating above 8.33ARM S are required. The chosen IGBT module, MUBW 20-06 A7, which Application: AC motor drives with Input Bridge D1135A, - D16 meets these requirements. can stand upRectifier to 600V and Input from single or three phase grid The Symbol integratedConditions temperature sensor consists of an NTC thermistor whose Maximum Ratings Three phase synchronous or asynchronous motor resistance gets smaller with the increasing of temperature. This sensor will V 1600 V electric braking operation be used Iin the design to prevent T = 80°C; sine 180° IGBTs damage from overheating. 30 A ● ● RRM FAV ● C IDAVM IFSM TC = 80°C; rectangular; d = 1/3 TVJ = 25°C; t = 10 ms; sine 50 Hz 29 400 A A Ptot TC = 25°C 120 W Features ● ● ● Symbol Conditions Characteristic Values (TVJ = 25°C, unless otherwise specified) min. typ. max. VF IF = 35 A; TVJ = 25°C TVJ = 125°C 1.4 1.4 IR VR = VRRM; TVJ = 25°C TVJ = 125°C 2.0 ● 1.7 V V ● 0.2 mA mA ● High level of integration - only one power semiconductor module required for the whole drive Fast rectifier diodes for enhanced EMC behaviour NPT IGBT technology with low saturation voltage, low switching losses, high RBSOA and short circuit ruggedness Epitaxial free wheeling diodes with Hiperfast and soft reverse recovery Industry standard package with insulated copper base plate and soldering pins for PCB mounting Temperature sense included 37 3.3. IGBT Driver Circuit 3.3 IGBT Driver Circuit Main purpose of the IGBT Driver is to provide enough power to the gate to make the IGBT switch properly. In this system, a compact IGBT driver module (6SD106EI) by Concept Technologie[5] was used at hand which is able to output high gate current of ±6A. This driver module can operate in two different modes: direct mode and half-bridge mode and is capable &$/( ULYHU driver signals or three pairs of driver of outputting either six independent signals with dead time. Safety is another important issue in IGBT driver Data Sheet 6SD106E selection and the driver 6SD106EI integrate a short circuit and over-current Block Diagram protection circuit inside. 6 ' Rth IGD VDD Rg Viso1 LDI Rth GND IGD Rg Viso2 VDC Viso1 PWM oscillator Viso2 GND Interface on Electronic Level Electrical Isolation Driver on Power Level Power Semiconductor (external) SCALE Driver Module Fig. 1 Block diagram shows 2 channels (i.e. one third) of the 6SD106E block Structure diagram shows two (i.e. Driver one third) ofModule: the 6SD106E 6SD106EI six-pack FigureThe3.3: ofchannels IGBT driver. There is only one PWM oscillator, whereas all other components are present in triplicate. 2 I n t e r n e t : of w wawtwo-channel . C T - C O N C E P T .driver com Figure3.3[5]Page shows the structure on the chip set. For a three-phase version, there is only one PWM oscillator, all other components are present in triplicate. For each channel, the driver module contains the electrical isolation between the control and power sides which prevents damage of control sides from high voltage and current of power sides. The internal block IGD contains an over-current and short-circuit protection circuit for the power transistors, a feed monitoring circuit as well as a status acknowledgement circuit. 38 6&$/( 'ULYHU Description and Application Manual +15V PWM Input 0V 3.3. IGBT Driver Circuit +15V Gate G2 -15V An electrically separated power supply can provide ±15V for the drive elec+15V tronics via an integrated DC/DC converter. Gate G1 -15V Fig. 8 Signals curves of the circuit as per Fig. 7 Mode Selection As mentioned above,mode the driver module Half-bridge with dead time6SD106EI can operate in two different modes. Two sets of complementary drive signals are needed for the In half-bridge mode, two channels are always operated as a half bridge. In this H-Bridgemode, inverter. Hence, in this specific system, the driver module will the SCALE driver can generate the required dead times directly in a range work in the which is able three pairs from half-bridge about 100ns upmode[5] to several microseconds. Only to twooutput external RC networks are of comrequired (see page 26 signals. for dimensioning). All power semiconductors can be turned off plementary IGBT drive by switching the release input (InB) to low. +VCC + DC Link Channel 2 4k7 C2 VCC (LDI) 4V7 Rth MOD IGD GND Reset VL/Reset Enable (TTL) Rth2 Rg G2 InB PWM Input (TTL) E2 InA Inductive load +VCC LDI 15k 10k 10k SO Channel 1 C1 SO2 Rth SO1 RC1 IGD RC2 GND 100p GND 100p GND GND SCALE Driver Module Rth1 Rg G1 E1 - DC Link Fig. 9 Application example for half-bridge mode with dead time generation Figure 3.4: HalfBridge Mode Selection Page to 20 select the Ihalf-bridge n t e r n e t : mode, w w w .the C T -pin CON C E P Tshould . c o m be conIn order MOD nected to GND and inputs pins of RC1 and RC2 must be connected to RC networks. The dimensioning will be discussed further in the next section. Dead Time Because any real power electronic devices do not turn on or off instantaneously, it is necessary to include a protection time, called the dead time, to avoid cross conduction of two switching devices in the same leg of the H-bridge inverter. The dead time can be set either in DSP software or in the driver module 6SD106EI by selecting the half-bridge operation mode. In this design, the dead time is set in the driver module. In the half-bridge mode, RC networks must be connected to the pins RC1 and RC2. According to Table.3.1, 22kΩ and 150pF are selected and the dead time is set to 2.1µs. 39 3.3. IGBT Driver Circuit R C typ. dead time 10k 47pF ≈ 200ns 10k 100pF ≈ 500ns 15k 120pF ≈ 1.1µs Description and Application Manual 22k 150pF ≈ 2.1µs 33k 220pF ≈ 4.6µs Pin Ex (emitter terminal) 6&$/( 'ULYHU The “x” in “Ex”3.1: standsValues for the number of the drive times channels of in multi-channel drivers. Table of the dead RC networks This terminal should be connected to the emitter or source terminal of the power transistor. The connection must be as short as possible and be run directly to the emitter or source terminal of the power element. This terminal should be used in Short-Circuit and Over-Current Protection modules with auxiliary emitters or an auxiliary source. This terminal is also used as the low end of the reference resistor Rthx. Where possible, this should be connected directly toofthethe terminal Ex ofmodule the driver. is equipped with a Vce monitoring Every channel driver circuit as illustrated in Figure 3.5[5]. A resistor (Rthx) is connected to the pin Rthx Pin Cx (collector as a reference. It definessense) the maximum voltage drop across the turned-on power transistor atstands which protection function of the drivers. drive circuit is The “x” in “Cx” for thethe number of the drive channels in multi-channel activatedThis and thus the power transistor is turned off. terminal is used to measure the voltage drop across the turned-on power transistor in order to ensure protection from short circuit and overload. It should be noted that it V+ 1,4mA V+ 150uA 4 OVERCURRENT Rm Dm (2 x 1N4007) Cx Ca 5 5WK[ MEASURING RGx Gx Rthx IGD 001 Ex SCALE Driver Module Fig. 13 Principle of the collector sense circuit Figure 3.5: Protection Circuit of Driver Module 6SD106EI Page 28 Internet: www.CT-CONCEPT.com If the voltage at Cx exceeds the voltage at Rthx, a Vce or under-voltage error occurs and the protection function is activated. The driver blocks the power semiconductor and accepts no drive signals. Meanwhile, the status outputs SOx (See Figure3.4, The “x” in “SOx” stands for the number of the drive channel in the driver module.) for the corresponding channels 40 3.4. Voltage Measuring Circuit are pulled to the logic low level which will be used in the DSP protection routine. The pin Cx must be never connected directly to the collector of the power transistor and a high-blocking diode is needed to protect the measuring terminal Cx. For 600V IGBT modules, one diode of type 1N 4007 is connected to pin Cx. It is recommended that the peak repetitive reverse voltage of these diodes be over-dimensioned by at least 40%. According to the application manual of 6SD106EI, a threshold voltage Dimensions LV 25-P (in mm. 1 mm = 0.0394 inch) Vth = 5.85V is recommended to use in IGBT drives. The reference resistor Rth can be calculated Bottom view Right viewas follows: Top view Rth = 3.4 Vth 5.85V = = 39kΩ 150µA 150µA Voltage Measuring Circuit (3.1) swiss made Voltage measuring circuit consists of a voltage sensor, voltage level shifter circuits and a voltage buffer. In the following section, these blocks will be introduced one by one in detail. Standard 00 or N° SP.. Year Week Voltage Sensor Secondary terminals Voltage is measured with a voltage transducer LV 25-P (LEM Inc.). Its range Terminal : supply voltage 12 .. 15 V output of this sensor goes up to 500V which is suited+ for 120V the RM S +and Terminal M : measure is a current with the ratio of 2.5 : 1 of the primary current. This module is Terminal - : supply voltage - 12 .. 15 V supplied by ±15V power and the connection is depicted in Figure 3.6. Connection Back view Figure 3.6: Voltage Sensor Connection Remarks Mechanical characteristics • General tolerance • Fastening & connection of primary ± 0.2 mm • IS is positive when VP is applied on terminal +HT. 2 pins • This is a standard model. For different versions (supply For xvoltage current proportional to the measured volt0.635 0.635 mmmeasurements, voltages, turnsa ratios, unidirectional measurements...), • Fastening & connection of secondary pins ∅ be 1 mm please contact age3must passed through an us. external resistor R1 which is selected by the • Recommended PCB hole 1.2 mm user and installed in series with the primary circuit of the transducer. In our case, R1 = 6kΩ was selected to measure 60V voltage which corresponds Instructions for use of the voltage transducer model LV 25-P 41 Primary resistor R 1 : the transducer’s optimum accuracy is obtained at the nominal primary current. As far as possible, R 1 should be calculated so that the nominal voltage to be measured corresponds to a primary current of 10 mA . Example: Voltage to be measured V PN = 250 V a) R 1 = 25 kΩ / 2.5 W, IP = 10 m A b) R 1 = 50 kΩ / 1.25 W, IP = 5 m A Accuracy = ± 0.8 % of V PN (@ TA = + 25°C) Accuracy = ± 1.6 % of V PN (@ TA = + 25°C) Operating range (recommended) : taking into account the resistance of the primary windings (which must remain low compared to R 1, in order to keep thermal deviation as low as possible) and the isolation, this transducer is suitable for measuring nominal voltages from 10 to 500 V. LEM reserves the right to carry out modifications on its transducers, in order to improve them, without previous notice. 3.4. Voltage Measuring Circuit to 10mA primary current in the testing experiments. For the real islanding test,a 17kΩ resistor will be used to measure 120VRM S voltage. Voltage Level Shifter A voltage level shifting circuit is required because the input voltage range of ADC module on the ezDSPlf2407A board used in this case is 0V − 3.3V which is different from that of the voltage sensor output. This function can be implemented with a two-level operational amplifier (OPA) circuit as shown in Figure 3.7. -15i VT LV 25-P 3 2 1 R1 240 R3 R2 30K C3 30K 150pF +15i C1 1 1 11 2 LM324AM U4A C2 R4 12K 9 6 5 7 2 LM324AM U4B 10 3 8 VM1 LM324AM U4C -15i -15I 100nF -15i 1 3 NC U5 LM4040A30IDBZR-3.0 +15i +15I 4 2 3 3.3K 100nF 4 R6 4 5.1K R5 11 +15i -15i 11 CN7 Figure 3.7: Voltage Measuring Circuit The current signal from LV 25-P is transformed into a voltage signal with R1, the range of which is −6V to 6V . The voltage across R1 is then offset with the first level OPA circuit and the output voltage of U4A is given by: R3 R3 Vout1 = VREF 1 + − VIN × (3.2) R2 R2 The voltage VREF is voltage across the precision voltage reference LM4040A and equals to 3V . Through the second level OPA circuit, the output voltage of U4B is scaled by: R5 Vout = − × Vout1 (3.3) R4 42 3.4. Voltage Measuring Circuit Plugging in values for the components in Figure 3.7, the scaling of signals from the primary current to DSP analog voltage is illustrated in Figure 3.8. Current Sensor Offset IN Coefficient out1 out Figure 3.8: Voltage Scaling from Primary current to DSP Analog Voltage Voltage Buffer The ADC module in DSP is an unbuffered multiplexed ratiometric ADC which contains a multiplexer, a sample/hold circuit and an ADC comparator. This type of ADC has no internal buffer amplifiers to introduece input offset and gain errors. The structure of ADC is shown in Figure 3.9. A disadvantage of this kind of on-chip ADC is that the sample capacitor within the ADC is directly charged by the external signal, and the charge left on the sample capacitor by the previous conversion of a channel can affect the accuracy of the channel currently being converted if inadequate settling time is allowed for a given source impedance. This phenomenon is referred to as channel-to-channel crosstalk[6]. sample MUX Figure 3.9: Structure of DSP On-chip ADC 43 3.5. Current Measuring Circuit In order to avoid the crosstalk, the source impedance of input signals should be much less than the input impedance of ADC. Therefore, a voltage buffer U4C in Figure 3.7 is added to decrease the output impedance of the scaling circuit so that ADC module can work properly. 3.5 Current Measuring Circuit The current measuring circuit is mostly the same as the voltage measuring circuit except a few changes of parameters. In this case, the current is measured with a LEM current transducer LT 100-S, the range of which goes up to ±100A. The output of this module is a current with a ratio of 1 : 1000 of the primary current. With the consideration of accuracy, the current should be close to 100A. The current measuring circuit is illustrated in Figure 3.10. -15i CT LT 100-S 3 2 1 R1 62 R3 C3 30K +15i 150pF C1 R6 2 1 +15i +15I 11 C2 R4 12K 9 6 5 7 2 LM324AM U1B 10 3 8 CM1 LM324AM U1C -15i -15I 100nF -15i 1 3 NC U5 LM4040A30IDBZR-3.0 1 LM324AM U1A 2 3 3.3K 100nF 4 4 5.1K R5 4 32K 11 +15i -15i R2 11 CN1 Figure 3.10: Current Measuring Circuit Substitution of parameters shown in the figure above into Eq.(3.2) and (3.3), a solution for scaling of current signals is depicted in Figure 3.11. In the testing experiments, the maximum current is below 3A and 32 turns of coil is used so that the range of measured current is close to 100A, the primary nominal current of the current sensor. 44 3.6. Temperature Sensor Current Sensor Offset Coefficient IN out1 out Figure 3.11: Current Scaling from Primary Current to DSP Analog Voltage 3.6 Temperature Sensor A temperature sensor circuit is provided on the PCB to achieve IGBT overheat protection. The temperature sensor used in this project is an NTC thermistor integrated in the IGBT module. According to the datasheet of this module, the IGBT operating temperature has a range from −40 ◦ C to 125 ◦ C and the NTC thermistor resistance equals to 300Ω at 125 ◦ C. Because the NTC thermistor resistance varies with temperature significantly, a circuit detecting the variation of thermistor resistance can be applied for temperature sensing. In the practical application, a failure recovery from an over-heat fault will take a period of time. Therefore, a hysteresis comparator circuit shown in Figure 3.12 is used as temperature sensor circuit to aviod frequent alarms. +5 5K R12 3K R13 R16 12 1.5K R14 Vin 5 U3A 2 VREF Vout 4 3 LM339AM C31 2.5v D4 300 R15 +5 100nF Figure 3.12: Temperature Sensor Circuit 45 3.7. Protection Circuits As depicted in Figure 3.12, in normal conditions, the temperature keeps below 125 ◦ C and the thermistor resistance R16 > 300Ω. The output of comparator U3A is zero. When the temperature rises above the maximum value, R16 is less than 300Ω and the output of U3A jumps up to high voltage level (VH = 5V ). Due to the hysteresis loop R12 , the lower limit of Vin is solved by: R16 R15 Vin min = VREF − × VH (3.4) R16 R15 + R12 Substitution of parameters shown in Figure 3.4 into Eq.(3.4) yields Vin min = 2.35V , which means Vout will not drop to zero unless Vin is less than 2.35V instead of 2.5V . Convert voltage to resistance and a range of resistance from 300Ω to 338Ω can be obtained. Looking up the characteristic of thermistor resistance versus temperature leads to a temperature hysteresis loop in Figure 3.13. ℃ ℃ ℃ Figure 3.13: Temperature Hysteresis Loop 3.7 Protection Circuits The pin P DP IN T A of ezDSPlf2407A[7] is used to prevent the system from overheat and over-current faults. If this pin is driven low, a fault occurs in DSP program and all PWM output pins will be put in the high-impedance state immediately. A logic circuit collecting all fault signals of the system is applied as an interface shown in Figure 3.14. 46 3.8. DSP Control Board U2A 6 U1A 5 4 6 PDPINT 2 1 SN74ALS21AD 1 Temperature Sensor SN74LVC2GU04 SO_W SO_V SO_U Figure 3.14: System Protection Circuit The pin “Temperature Sensor” is the output of temperature sensor circuit which is in high voltage level when an overheat fault occurs. In order to fit the logic level of P DP IN T A, an inverter SN74LVC2GU04 is used to invert the logic level. Pins of “SO x” (x = U, V orW ) are the status pins of driver module. When an over-current fault occurs, the logic levels of “SO x” are pulled down to zero and a P DP IN T A interrupt is activated in DSP program. 3.8 DSP Control Board An ezDSPlf2407A board manufactured by TI company is employed in this case as a control unit. It contains a TMS320LF2407A digital signal processor, an on-board JTAG connector which provides interface to emulators for program debugging, and other peripherals. The major features of the TMS320F2407A include: • A high-speed CPU with 40MHz clock rate and capable of processing 40 million instructions per second (MIPS). • 64K words of on-chip data/program RAM, 32K words of on-chip program ROM or Flash EEPROM, 64K words of program, 64Kwords of data and 64K words of I/O space of addressing space. • Sixteen multiplexed analog inputs 10-bit ADC core with built-in Sample and Hold (S/H) circuit and fast conversion time (S/H + Conversion): 375ns. It supports up to four trigger sources for start-ofconversion (SOC) sequence and autosequencing function. • Two Event Manager (EV) modules provide a broad range of functions and features that are particularly useful in motion control and motor control applications. 47 3.9. Power Supply • 40 multiplexed general-purpose I/O pins. • Watchdog Timer which monitors software and hardware operations, and implements system reset functions upon CPU disruption. • Controller Area Network (CAN). • Serial communications interface (SCI). The EV modules in this DSP are very versatile and the major features are: • Two general-purpose (GP) timers. • Three general-purpose up and up/down timers, each with a 16-bit compare unit capable of generating one independent PWM output. • Pulse-width modulation (PWM) circuits that include space vector PWM circuits, dead-band generation units, and output logic. • Three 16-bit simple compare units capable of generating 4 independent PWM outputs. • Three capture units. • Quadrature encoder pulse (QEP) circuit. 3.9 Power Supply The printed circuit boards are supplied by a commercial switching power supply VOF-65-15. Two other commercial dc dc converters were used to provide proper supply voltage for the circuits. A 5V dc dc converter (CC101205SF-E) supplies the ezDSP board, the temperature sensor circuit and the protection circuit; a ±15V supply (CC10-1212DF-E) provides power for the op amp circuits; the switching power supply (VOF-65-15) feeds the driver circuits. All these converters were equipped with the suggested input and output filters as stated in the corresponding data sheets and application notes. The input filter consists in all three cases of a Π-type filter, consisting of two capacitors and one inductor. On the output side, there are only smoothing capacitors and a zener diode, to protect devices from an unexpected overvoltage. The complete power supply circuit is shown in the Figure A.1 in the Appendix A.1. 48 Chapter 4 Software Implementation The hardware system of AC electronic loads has been described in the previous chapter. In this chapter, a software applied to implement the control algorithm will be introduced in detail. All implementations of software are based on TI’s (Texas Instruments) DSP TMS320LF2407A and the associated assembly language. 4.1 Program Overview As the previous introduction, the DSP TMS320LF2407A has many peripherals and is enough for this application. The peripherals used in the control program contains: event managers (EV), general purpose timers 1 and 2, the PWM generation unit, the 16-channel ADC unit and the digital I/O ports module[8],[9],[10]. The AC electronic loads control program has several major functions: accurate voltage and current acquisition, current command calculation and PWM duty cycle update. With the consideration of software portability, the program was divided into several sub-function blocks. The flowchart of the control program is shown in Figure 4.1. The control program consists of four routines: • System initialization and Main Loop. All initializations are finished in this routine and an endless loop is waiting for different ISRs to achieve all major features of AC electronic loads. • Timer1 Interrupt service routine (ISR). This routine calculates the required current command and update PWM duty ratios. • Timer2 Interrupt service routine (ISR). Average voltage and current sampling are implemented in this routine. 49 4.1. Program Overview • P DP IN T A Interrupt service routine (ISR). Over-heat and over-current protection is achieved here and makes the whole system safer and more reliable. Start T1 ISR Start Initialize Variables and Constants Start Timer2 and ADC System Configuration First Sampling? Yes No Setup Watchdog and I/O Pins Read ADC Data Calculate Average Values Setup EV, GP Timers and PWM Pins Voltage/Current Normalization Configure ADC Current Command Calculation PI Control PWM Duty Update Interrupts Initialization Back to Main Loop Clear INT Flags Enable Interrupts T2 ISR Start Main Loop Reset T2 Interrupt Timer1 ISR PDPINT ISR Timer2 ISR Read Voltage/Current Data Save Accumulation Back to Main Loop Figure 4.1: Program Flow Chart The following sections will explain the details of these routines. 50 4.2. Variables Normalization 4.2 Variables Normalization TMS320LF2407A is a 16-bit fixed-point DSP and unable to represent floatingpoint numbers directly. Therefore, the normalization of variables is required in the program so that the real values of voltage and current signals can be represented correctly. 4.2.1 Q Format In fixed-point DSPs, all the numbers must be represented as a collection of bits. Each bit represents either “0” or “1”, hence the number system naturally used in microprocessors is the binary system. In order to represent fractional numbers, fixed-point system requires the programmer to create a virtual decimal place in between two bit locations for a given number. Qformat is such kind of fractional fixed-point representation suitable for DSPs algorithms[11]. The labeling convention of Q-format is as follows: Q [QI] . [QF ] (4.1) Where QI is the number of integer bits and QF is the number of fractional bits. In a 16-bit DSP used in our system, the sum of QI and QF equals to 15 and the MSB is the sign bit. In addition and subtraction of two Q-format numbers, the fixed-point decimal places must be aligned and an appropriate dynamic range has to be chosen to handle the overflow of the addition. In multiplication, the number of integer and fractional bits in the product is the sum of the corresponding multiplier and multiplicand Q-format numbers as described the following equations: QIproduct = QImultiplicand + QImultiplier QFproduct = QFmultiplicand + QFmultiplier 4.2.2 (4.2) Normalization of Voltage and Current As mentioned above, there is a trade-off between dynamic range and precision for Q-format numbers. Higher precision leads to a narrower dynamic range, hence the real voltage and current values need to be normalized and thus the Per Unit (PU) system will be applied to achieve the high precision as well as the wide dynamic range. Because the measured current in the experiments is less than 3A, 32 turns of the current will be measured to achieve higher accuracy. 51 4.3. System Initialization In this system, the base values of variables are as follows: Vbase = 30V Ibase = 100A/(4 · 32turns) = 0.78125A |Z|base = Vbase /Ibase = 38.4Ω With the consideration of both precision and dynamic range, the PU values will be represented in Q[3].[12] format (or Q12). The maximum values of current and voltage in the experiments are imax = 3.125A and Vmax = 60V . When the current reaches the maximum value, the ADC result ibinary is 512. The corresponding per-unit value of the Q-format number is given by iP U Q12 = (imax /ibase ) × 212 = 214 , and a conversion ratio Kcurrent = iP U Q12 ibinary = 214 512 = 32. All measured values of current can be converted into the PU values in Q12 format: iP U Q12 = Kcurrent × ibinary (4.3) Following the same steps, the voltage conversion ratio Kvoltage equals to 16 and the conversion formulas is vP U Q12 = Kvoltage × vbinary 4.3 (4.4) System Initialization To implement the specific functions, the initialization module performs the following tasks. 4.3.1 Event Manager Modules Initialization As mentioned above, event manager modules are the most significant function blocks which provide a flexible method for controlling both dedicated I/O and shared pin functions. The following blocks are initialized in this module: PWM generation block To drive the H-bridge inverter, four PWM pins 1-4 are employed and the polarity of PWM pins are set in Compare Action Control Register A (ACTRA) as follows: • PWM 1 and 4 are active high. 52 4.3. System Initialization • PWM 2 and 3 are active low. PWM 1-4 corresponds to switches 1-4 and “active high” means the output of PWM pins are set to one when a compare match happens. Additionally, asymmetric PWM waveforms are used in this program and associated timers are set to continuous-up count mode. General-Purpose (GP) Timers 1&2 GP timer 1 is employed in the software to control the switching frequency of PWM signals. Usually, PWM frequencies are in the range of 10kHz. In this project, a PWM frequency of 20kHz has been chosen. GP timer 2 is used to setup the voltage and current sampling frequency. In order to achieve higher precision, the sampling frequency is set to 160kHz. 4.3.2 ADC Module Initialization ADC Clock Prescaler ADC Clock Prescaler ADC conversion time consists of two time segments: Sample/Hold (S/H) window and conversion window, as shown in Figure 4.2. The S/H window 7.3 ADC Clock Prescaler can be tailored to accommodate the variation in source impedances by the The S/H block in the 240xA ADC can be tailored to accomodate the variation ACQ PS3-ACQ PS0 inbits and the CPS in the source impedances. This isbit achieved by theADCTR1 ACQ PS3−ACQregister. PS0 bits and the ACQ CPS bit inPS3-ACQ the ADCTR1 register. analog-to-digital conversion process In this project, the PS0The bits are set to “0111b” and the can be divided into two time segments, as shown in Figure 7−6. CPS bit is set to “0b” to accommodate a 2290Ω source impedance. Figure 7−6. ADC Conversion Time •••••• S/H window (2 * PS) Conversion (11 * ACLK) 1 complete ADC conversion PS = a prescaled CPU clock Figure 4.2: ADC Conversion Time PS will be the same as the CPU clock if the prescaler = 1 (i.e., ACQ PS3−ACQ PS0 bits are all zero) and if CPS = 0. For any other value of the prescaler, the magnitude of PS will be magnified (effectively increasing the S/H window time) as described by the “Acquisition Time Window” column in the bit description for ACQ PS3−ACQ PS0. If the CPS bit is made 1, the S/H window is doubled. This doubling of the S/H window is in addition to the “stretching” provided by the prescaler. Figure 7−7 shows the role played by the various prescaler bits in the ADC module. Note that PS and ACLK will be equal to CPU clock if CPS = 0. 53 Analog-to-Digital Converter (ADC) 7-17 4.4. Interrupt Service Routine (ISR) ADC Conversion Sequencer The ADC module in TMS320LF2407A contains a 16-state sequencer which can perform a series of conversions without external intervention and the conversion sequence is stored in the ADC input channel select sequencing control registers (CHSELSEQn). In this project, a timer2 underflow interrupt is used as a trigger to start an auto-conversion sequencer. The number of auto-conversions is two and a single sampling of voltage and current signals is performed for each trigger. 4.3.3 I/O Ports Module Initialization The commercial driver module contains an enable pin InB (shown in Figure 3.4) which is able to block all PWM channels if pulled to zero. Therefore, the P DP IN T A interrupt service routine employs the digital I/O port IOPB4 as an enable pin for the commercial driver module to achieve the protection function. The pin IOPB4 is configured as an output port with the I/O mux control register A (MCRA) and set to zero when a fault occurs. 4.4 Interrupt Service Routine (ISR) As shown in Figure 4.1, three ISRs are applied to perform the current control and protection functions. Timer 1&2 ISRs The major features of Timer 1&2 ISRs have been described in Figure 4.1. Here the discussion will be focused on the time sequence between two ISRs. In these two routines, the timer underflow events are configured to generate interrupts which means an interrupt event will occur once the counter of timer 1 or 2 reaches “0000h”. PWM Figure 4.3: Timer 1 and 2 ISR Sequence 54 4.5. Generation of Sine Wave Figure 4.3 depicts the time sequence of two ISRs (TxU ISR represents the associated interrupt routine, x = 1or2). It can be seen that each underflow interrupt event of timer2 will perform two conversions: one is voltage sampling and the other is current sampling. Due to the ADC autoconversion sequencer, the pointer of the sequencer will move to the address of next two channels in the sequence. Because the frequency of timer2 is eight times of the timer1’s, 16 conversions are implemented in one PWM period and the pointer of the ADC sequence is set to the original address automatically. Additionally, timer 1 and 2 must be synchronized to keep the correct time sequence between two timers. Therefore, the timer2 is reset and started at the beginning of each T1U ISR. P DP IN T A ISR With the protection circuit shown in Figure 3.14, the system is capable of monitoring the overheat and overcurrent faults. Figure 4.4 shows that the pin IOPB4 is set to low level which will block all PWM outputs of the driver module so that the whole system is able to survive a severe accident. Start PDPINTA ISR Save Interrupt Context Set IOPB4 to Zero Restore Interrupt Context End ISR Figure 4.4: Flowchart of P DP IN T A ISR 4.5 Generation of Sine Wave To simulate an inductive or capacitive load with different input sources, a general approach based on LdiL /dt = vL or Cdvc /dt = ic might be possible. In each PWM switching period, ∆iL or ∆vc can be calculated using these 55 4.5. Generation of Sine Wave two equations and after comparison with the actual ∆iL or ∆vc , a new duty cycle D can be calculated and used to control the current to track the required load characteristics. However, this prototype only focuses on the 60Hz sinusoidal voltage source and thus a sine look-up table will be implemented to achieve the inductive or capacitive loads. Compared with the general approach, the sine look-up table needs less operation time. The table contains 512 words to represent sine values of phase angles in the range of [0◦ , 360◦ ]. As a result, the resolution on phase angle of electronic loads is limited to 360/512 = 0.703125◦ . In the binary system of DSP, the phase angle θ varies from 0 to 4095. As only 512 words are available to represent this range, θ is divided by 8 and stored into the variable index that will be used to address the lookup table. The content of the table row pointed by the index is fetched in indirect addressing mode via AR5 auxiliary register (Figure 4.5). This content coded in Q12 format is stored in the variable sincos that will be used in the phase angle calculation of electronic loads. Sine Table Address 0 50 101 0 ... θ >>3bits Index 4095 4096 4096 4096 4095 π/2 ... 101 50 0 65486 65435 π ... 61441 61440 61440 61440 61441 3π/2 ... 65435 65486 2π Figure 4.5: Sine Look-up Table Note that 90◦ is added to θ to get the cosine value of the phase angle. This operation corresponds to addition of 128 (512/4) to the value of index. 56 4.6. Digital PI Control 4.6 Digital PI Control The PI controller is a generic control loop feedback mechanism widely used in industrial control systems. It is an effective approach in nearly all types of the feedback system, especially for systems containing a single pole. In our control program, a digital PI control will be implemented. An analog PI controller can be expressed by: u (t) = KP e (t) + 1 TI t e (t) dt + u0 (4.5) 0 Where KP is the proportional gain, TI is the time constant of the controller and thus the integral gain KI = KP /TI ; e (t) and u (t) are the input and output of the controller respectively. The proportional term means that the system responds in proportion to the deviation from the set point and determine the sensitivity of the controller to the error. A high proportional gain results in a large change in the output for a given change in the error. If the proportional gain is too high, the system can become unstable. In contrast, a small gain results in a small output response to a large input error, and a less sensitive controller. The integral term means that the system will respond in proportion to the integral of the error over time. A small time constant accelerates the movement of the process towards setpoint and eliminates the residual steady-state error that occurs with a proportional only controller. However, since the integral term is responding to accumulated errors from the past, it can cause the present value to overshoot the setpoint value. In order to implement the PI controller in a DSP, it has to be converted into digital form. A common way of doing this is to discretize the controller to approximate the continuous time derivatives as follows: k T pwm uk = KP ek + (4.6) ej + u0 TI j=0 Where k represents the kth sampling, Tpwm is the switching period and u0 is the initial value of the output. Considering the kth and (k − 1)th samplings, the change of output at the kth sampling can be derived: ∆uk = uk − uk−1 = KP (ek − ek−1 ) + Tpwm KI ek (4.7) Upon collecting items, one obtains uk = uk−1 + (KP + Tpwm KI ) ek − KP ek−1 (4.8) 57 4.7. RLC Algorithm With Eq.(4.8), a digital PI controller can be implemented easily in the DSP program. Additionally, a saturation block is also implemented in the program to protect the control system from the integration saturation. 4.7 RLC Algorithm The device we designed outputs different kinds of load characteristics by controlling the inductor current. Hence, a current command calculation block shown in Figure 4.1 is required for different electronic loads and the algorithms of different loads will be described here. For the resistive loads, the current command calculation is quite straightforward and easy to implement in the DSP program. As a resistor, the current flowing through and the voltage across it are in phase, and thus the inductor current iL is proportional to the source voltage Vg at each sampling point. Therefore, in the software, the current command can be obtained by simply dividing Vg by the simulated resistor Rsim . When it comes to the inductive or capacitive loads, however, such a simple algorithm couldn’t be used any more due to the phase difference between iL and Vg . Because the inductive and the capacitive loads have similar V − I characteristics, only the algorithm of the inductive loads will be discussed here in detail. Phase Detection Block In order to simulate the inductance V − I characteristic with respect to the reference voltage Vg , it is necessary to detect the phase angle of Vg so that the control software can output the correct current with 90-degree phase delay (90-degree lead for the capacitance loads). In our program, the first cycle of Vg is used for phase detection without any other operations. As shown in Figure 4.6, checking the zero points of Vg is an effective approach of the phase detection. In the program, two variables are defined to save the previous and current voltage values: V1 is the old value and V2 is the current value. When V1 > 0 and V2 ≤ 0, Z1 can be found and the phase angle is 90◦ (vg = Vg × cos (ωt)); when V1 < 0 and V2 ≥ 0, Z2 can be found and the phase angle is 270◦ . Once the zero-crossing and the phase angle are detected, the correct phase shift φ can be generated using ωt + φ (φ = 90◦ for the inductive loads). 58 4.7. RLC Algorithm g L 2 1 1 2 1 2 Figure 4.6: Phase Detection Diagram Frequency Detection Block In the first cycle of Vg , the frequency is also detected to calculate the phase angle using θ = ωt. Digitalizing this equation, one obtains θ = ω · ∆t = 2πf · N · Tpwm (4.9) Where f is the frequency of Vg , Tpwm is the switching period and N represents the Nth switching cycle. The frequency f also equals f= 1 Nperiod × Tpwm (4.10) Where Nperiod represents the period of Vg which can be got by Nperiod × Tpwm . Insert Eq.(4.10) into Eq.(4.9) and θ can be derived θ = 2π × N Nperiod (4.11) In the program, once the first zero point of Vg is found, a counter F RQ CN T will be started immediately until the second zero point occurs and the period Nperiod is equal to 2 × F RQ CN T . Then, the phase angle can be derived using Eq.(4.11) and the program will start to generate the correct current command from the second zero point. Magnitude Detection Block The current command icmd equals Icmd · cos (ωt + φ), where Icmd is the magnitude of the current and can be obtained by Icmd = Vg Z (4.12) 59 4.7. RLC Algorithm Z represents the impedance of the inductive load (or the capacitive load) and the voltage magnitude can be detected in the first cycle of Vg by checking the maximum value of Vg . Figure 4.7 shows the flowchart of the current calculation block. The variable P hD F LAG indicates whether the phase detection is finished or not and the program will jump over the current calculation block unless P hD F LAG = 0. PhaseDetection Start PhD_FLAG=0? Yes No 1st. Zero Point Found? No Yes FRQ_CNT ++ Save |Vg|max If larger than the previous value No 2nd. Zero Point Found? Yes Calculate Icmd Using Sine Look-up Table END Figure 4.7: The Flow Chart of the Current Calculation Block 60 Chapter 5 Experimental Results In the previous discussions, the principles, hardware and software implementations of this AC electronic load system have been introduced in detail. The following chapter will show a series of experimental results to prove the principles and real hardware system. 5.1 Experimental Conditions On the basis of theory study and software simulation a prototype is researched and developed shown in Figure 5.4. Figure 5.1: Control Circuit Board Figure 5.2: Power Circuit Board Figure 5.3: Circuits Figure 5.4: Experimental System The Whole 61 5.2. Resistive Load Results The experimental system consists of: • A system control board including: voltage and current sensor circuits, protection circuits, an ezDSP board and power supply circuits. • A IGBT-based power board containing: an IGBT driver circuit and an H-bridge IGBT inverter. • A 3-phase AC source and a 3-phase transformer used to adjust the voltage level of AC source. • A DC power supply HPD 30-10 by XANTREX which can provide 60V and 3A DC source • Two series inductors, 2.6mH and 26mH, are used for ripple studies. • A 315Ω power resistor is used as the parallel resistor. In the experiments, a 4-channel 500MHz Oscilloscope ( Lecroy 6050A) is used to measure and save the results. 5.2 Resistive Load Results With the consideration of simplicity and typicality, the electronic load experiments will start with resistance loads. The peak voltage of AC source Vgpeak = 10V , the voltage ratio k = 1.3 and thus DC bus voltage Vdc = 13V . 5.2.1 The Range of Resistance Loads at L = 2.6mH In the PSIM simulations, it is found that the value of the series inductor L has to be changed with the variation of electronic load impedances. According to Table 2.8, an inductor of 2.6mH is selected to achieve resistance load experiments, the range of which is from 0.5p.u. to 2p.u.. The experimental results are depicted in Figure 5.5 to 5.10. The variable VRsim shown in following figures represents the voltage across Rsim which is calculated by VRsim = iL × Rsim and supposed to equal to the AC source Vg . It can be seen that VRsim is mostly same as Vg , both in magnitude and in phase. On the other hand, the current ripple increases with increasing of Rsim and introduces more high frequency interference. The current ripples for different Rsim s can be obtained in the following Figures 5.6, 5.8 and 5.10. 62 5.2. Resistive Load Results V Rsim and V (Rsim = 0.5pu) g VRsim and Vg (Rsim = 0.5pu) 10 VRsim (V) V Rsim (V) -4 5 0 -5 -10 -0.1 -0.08 -0.06 -0.04 -0.02 -5 -6 -7 -8 -9 -0.0634 0 10 ΔiL/2IL =0.22 -0.0633 -0.0632 -0.0631 -0.063 -0.0629 -0.0628 -4 Vg (V) g V (V) 5 0 -5 -10 -0.1 -0.08 -0.06 -0.04 -0.02 -6 ΔVg/Vg =0.26 -8 -10 -0.0634 0 -0.0633 -0.0632 Figure 5.5: Rsim fpwm = 20kHz V Rsim = 0.5p.u. at g VRsim (V) (V) Rsim V 0 -5 -10 -0.04 -0.02 0 ΔiL/2IL=0.33 -0.0303 -0.0302 -0.0301 -0.03 -0.0299 -0.0298 -0.0303 -0.0302 -0.0301 -0.03 -0.0299 -0.0298 -5 -6 5 Vg (V) g -8 -0.0304 10 V (V) -6 -10 -0.06 15 0 -5 -7 -8 -9 -10 -0.08 -0.06 -0.04 -0.02 -10 -0.0304 0 Time (seconds) Time (seconds) Figure 5.7: Rsim = 1p.u. at fpwm = 20kHz Figure 5.8: Rsim = 1p.u. VRsim and Vg (Rsim = 2pu) Current Ripple at -5 5 VRsim (V) Rsim (V) 10 V ΔVg/Vg =0.18 VRsim and Vg (Rsim = 2pu) 15 0 -5 -10 -15 -0.1 -0.0628 VRsim and Vg (Rsim = 1pu) -4 5 -15 -0.1 -0.0629 and V (Rsim = 1pu) 10 -0.08 -0.063 Figure 5.6: Current Ripple at Rsim = 0.5p.u. 15 -15 -0.1 -0.0631 Time (seconds) Time (seconds) -0.08 -0.06 -0.04 -0.02 -10 -15 0 15 ΔiL/2IL=0.42 -0.0359-0.0358-0.0357-0.0356-0.0355-0.0354-0.0353-0.0352-0.0351-0.035 -7 -8 5 Vg (V) g V (V) 10 0 -5 -10 -15 -0.1 -9 ΔVg/Vg =0.22 -10 -0.08 -0.06 -0.04 -0.02 0 -11 -0.0359-0.0358-0.0357-0.0356-0.0355-0.0354-0.0353-0.0352-0.0351-0.035 Time (seconds) Time (seconds) Figure 5.9: Rsim = 2p.u. at fpwm = 20kHz Figure 5.10: Rsim = 2p.u. Current Ripple at 63 5.2. Resistive Load Results Table 5.1 shows the comparisons of current ripple between the results of PSIM simulations and the experimental results. It can be seen that the experimental results have 30% to 40% error with respect to the simulation results. The main reason is because of the noise introduced by the AC source Vg and the noise ∆Vg /Vg shown in Figure 5.6, 5.8 and 5.10 is around 20%. Rsim (p.u.) 0.5 1.0 2.0 Simulation Results 0.13 0.19 0.30 Experimental Results 0.22 0.33 0.42 Error(%) 40.9 42.4 28.6 Table 5.1: Comparisons of Current Ripple between Simulation Results and Experimental Results 5.2.2 Switching Frequency and Current Ripple According to Eq.(2.33), the PWM switching frequency is inversely proportional to the current ripple. In this section, the effect of fpwm will be studied at 10kHz and 20kHz separately and Rsim is set to 2p.u.. Figures 5.11 to 5.14 show the current ripples at different frequencies. VRsim and Vg (Rsim = 2pu & f pwm = 10kHz) 20 0 -10 -20 -0.1 -0.08 -0.06 -0.04 -0.02 -10 -15 ΔiL/2IL=0.95 -0.0318 -0.0316 -0.0314 -0.0312 -0.031 -0.0318 -0.0316 -0.0314 -0.0312 -0.031 -6 Vg (V) Vg (V) -5 -20 -0.032 0 10 0 -8 -10 -10 -0.1 VRsim and Vg (Rsim = 2pu & f pwm = 10kHz) 0 VRsim (V) VRsim (V) 5 10 -0.08 -0.06 -0.04 -0.02 0 Time (seconds) Figure 5.11: Rsim = 2p.u. at fpwm = 10kHz and L = 2.6mH -0.032 Time (seconds) ΔVg/Vg =0.25 Figure 5.12: Current Ripple at Rsim = 2p.u. and fpwm = 10kHz Table 5.2 shows the comparisons of experimental results at different switching frequencies. It can be derived that the current ripple at 10kHz is about 2.6 times of the ripple at 20kHz.Considering the measurement errors 64 5.2. Resistive Load Results VRsim and Vg (Rsim = 2pu & f pwm = 20kHz) VRsim and Vg (Rsim = 2pu & f pwm = 20kHz) -5 VRsim (V) VRsim (V) 10 0 -10 -0.1 -0.08 -0.06 -0.04 -0.02 -10 -15 0 ΔiL/2IL=0.42 -0.0354-0.0353-0.0352-0.0351-0.035-0.0349-0.0348-0.0347-0.0346-0.0345 -6 Vg (V) Vg (V) 10 0 -8 -10 ΔVg/Vg =0.22 -10 -0.1 -0.08 -0.06 -0.04 -0.02 0 Time (seconds) Figure 5.13: Rsim = 2p.u. at fpwm = 20kHz and L = 2.6mH -12 -0.0354-0.0353-0.0352-0.0351-0.035-0.0349-0.0348-0.0347-0.0346-0.0345 Time (seconds) Figure 5.14: Current Ripple at Rsim = 2p.u. and fpwm = 20kHz and the noise introduced by Vg , the experimental result agrees well to the theoretical value “2 times”. fP W M (kHz) 10 20 Simulation Results 0.6 0.3 Experimental Results 0.95 0.42 Error(%) 36.8 28.6 Table 5.2: Comparisons of Current Ripple at 20kHz 5.2.3 Series Inductor L and Current Ripple The inductance of the series inductor L is also supposed to be inversely proportional to the current ripple and the associated experimental results are illustrated in Figure 5.15 to 5.18 to prove this relationship. For the experiments, two inductors of 2.6mH and 26mH are selected and the switching frequency equals to 20kHz and Rsim is still set to 2p.u.. It can be seen that the current ripple with the 2.6mH inductor is 7.7 times of the current ripple with the 26mH inductor which has about 20% error compared with the theoretical value “10 times”. With the consideration of measurement errors and the noise of the input source, the series inductance is inversely proportional to the current ripple(Table 5.3). 65 5.2. Resistive Load Results L(mH) 26 2.6 Simulation Results 0.03 0.30 Experimental Results 0.06 0.46 Error(%) 50.0 34.8 Table 5.3: Comparisons of Current Ripple at 20kHz VRsim and Vg (Rsim = 2pu & L = 2.6mH) VRsim and Vg (Rsim = 2pu & L = 2.6mH) -5 VRsim (V) VRsim (V) 10 0 -10 -0.1 -0.08 -0.06 -0.04 -0.02 -10 ΔiL/2IL=0.46 -15 0 -0.0354 -0.0352 -0.035 -0.0348 -0.0346 -0.0344 -0.0352 -0.035 -0.0348 -0.0346 -0.0344 -6 Vg (V) Vg (V) 10 0 -8 -10 ΔVg/Vg =0.22 -10 -0.1 -0.08 -0.06 -0.04 -0.02 -12 -0.0354 0 Time (seconds) Figure 5.15: L = 2.6mH, Rsim = 2p.u. and fpwm = 20kHz Figure 5.16: Current Ripple at L = 2.6mH VRsim and Vg (Rsim = 2pu & L = 26mH) 10 VRsim and Vg (Rsim = 2pu & L = 26mH) VRsim (V) 5 0 -5 -10 -0.1 -0.08 -0.06 -0.04 -0.02 -7 -8 ΔiL/2IL=0.06 -9 -0.048 0 10 -5 5 -6 VRsim (V) VRsim (V) -6 VRsim (V) Time (seconds) 0 -5 -10 -0.1 -0.08 -0.06 -0.04 -0.02 0 Time (seconds) Figure 5.17: L = 26mH, Rsim = 2p.u. and fpwm = 20kHz -0.0478 -0.0476 -0.0474 -0.0472 -0.047 -7 -8 ΔVg/Vg =0.18 -9 -10 -0.048 -0.0478 -0.0476 -0.0474 Time (seconds) -0.0472 -0.047 Figure 5.18: Current Ripple at L = 26mH 66 5.3. Inductive Load Results 5.3 Inductive Load Results The principles and PSIM simulations of inductance have been introduced previously and the experimental results will be shown in this section. From Table 2.9, the parameters for different inductance ranges can be found. Here, two inductors of 2.6mH and 26mH are selected, the switching frequency is 20kHz and the DC bus voltage Vdc = 1.5 × Vg . The range of Inductance loads at L = 2.6mH First of all, the range of inductance loads at L = 2.6mH will be depicted in Figures 5.19 and 5.20. VLsim and Vg (0.5pu) 6 V Lsim V Voltage (V) 4 g 2 0 -2 -4 -6 -0.1 -0.08 -0.06 -0.04 -0.02 0 Time (seconds) Figure 5.19: Lsim = 0.5p.u. at fpwm = 20kHz Voltage (V) VLsim and Vg (1pu) 6 V 4 V Lsim g 2 0 -2 -4 -6 -0.1 -0.08 -0.06 -0.04 -0.02 Time (seconds) Figure 5.20: Lsim = 1p.u. at fpwm = 20kHz 67 5.3. Inductive Load Results The variable VLsim in the figures represents the voltage across Lsim and can be calculated by VLsim = iL × Lsim . It can be found that the magnitude of VLsim is mostly equal to Vg ’s and has 90 degrees phase delay to Vg . Similarly with the resistance loads, the current ripple goes up with the increasing of the impedance Lsim . Series Inductor L and Current Ripple In order to obtain better performance, the current ripple should be reduced to a certain range and the most efficient method is increasing the series inductance. In the experiments, 2.6mH and 26mH inductors are used separately to show the improvement on current ripples (Figure 5.21 to 5.24) Figure 5.21: L = 2.6mH ,Lsim = 1p.u. and fpwm = 20kHz Figure 5.22: L = 26mH ,Lsim = 1p.u. and fpwm = 20kHz VLsim Vg 1pu) VLsim and Vgand (26mH VLsim Vg (1pu) VLsim and and Vg (2.6mH 1pu) 4.7 Voltage (V) Voltage (V) 6 ΔiL/2IL=0.5 5 4 3 ΔiL/2IL=0.05 4.5 Vg 4.4 4.3 4.2 2 -0.0885 VLsim 4.6 -0.0885 -0.0884 -0.0884 Time (seconds) -0.0883 Figure 5.23: Current Ripple at L = 2.6mH -0.0954 -0.0953 -0.0953 Time (seconds) -0.0952 Figure 5.24: Current Ripple at L = 26mH It can be seen in the figures above that the current ripple decreases significantly with the increasing of series inductance. The current ripple at 68 5.4. Capacitive Load Results L = 2.6mH is 10 times of the ripple at L = 26mH, which means the current ripple is inversely proportional to the series inductance. 5.4 Capacitive Load Results The capacitance loads have the same features as the inductance loads and the only difference is that the phase angle of capacitance loads leads the AC source Vg by 90 degrees instead of 90 degrees delay. The following figures depict the range of capacitance loads at L = 2.6mH. VCsim and Vg (0.5pu) 6 V Csim Voltage (V) 4 V g 2 0 -2 -4 -6 -0.09 -0.08 -0.07 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0 Time (seconds) Figure 5.25: Csim = 0.5p.u. at fpwm = 20kHz Voltage (V) VCsim and Vg (1pu) 6 V 4 V Csim g 2 0 -2 -4 -6 -0.06 -0.05 -0.04 -0.03 -0.02 -0.01 0 0.01 0.02 0.03 Time (seconds) Figure 5.26: Csim = 1p.u. at fpwm = 20kHz Figures 5.27 to 5.30 show the the relationship between current ripple and series inductance. From the experimental results, it can be derived that the 69 5.4. Capacitive Load Results current ripple is inversely proportional to the series inductance as same as the conclusion of inductive loads. Figure 5.27: L = 2.6mH ,Csim = 1p.u. and fpwm = 20kHz Figure 5.28: L = 26mH ,Csim = 1p.u. and fpwm = 20kHz VCsim Vg (1pu)1pu) VLsim andand Vg (2.6mH 4.8 Vg ΔiL/2IL=0.95 3 Voltage (V) Voltage (V) andVV (26mH1pu) 1pu) VVLsim and (2.6mH gg Csim VCsim 4 2 1 VCsim 4.6 Vg ΔiL/2IL=0.1 4.4 4.2 4 -0.0397 -0.0397 Time (seconds) -0.0396 Figure 5.29: Current Ripple at L = 2.6mH -0.0774 -0.0774 -0.0773 -0.0772 Time (seconds) -0.0772 Figure 5.30: Current Ripple at L = 26mH 70 Chapter 6 Conclusions and Future Work In this thesis, a DSP-controlled programmable AC electronic loads used in the unintentional islanding test has been implemented. A number of issues regarding the principles and the specifications of electronic loads have been investigated. In this final chapter a summary of the contributions contained in this thesis is made and some conclusions from this work are also presented. In the final section ideas for extending the results of this thesis are presented. 6.1 Synopsis The chapter 2 introduced the basic principles of the AC electronic loads designed in this project. A schematic of AC electronic loads based on the H-bridge dc-dc converter with a regulated DC bus voltage was presented first. Then, some theoretical analysis of the limitations of this circuit was developed and the relationship between the parameters of the circuit and the range of the electronic loads was introduced to determine the specifications of this circuit for different electronic loads. At the end, some simulations using PSIM were performed to verify the specifications for different ranges of electronic loads obtained from the principles. Additionally, the ac smallsignal model of this system was derived and a set of PI parameters was obtained using MATLAB and PSIM. In the following chapter the hardware implementation was presented in detail. The features of the main parts, such as the IGBT module, the driver module and the ezDSP board, were introduced as well as the configurations of these parts. Then, the voltage and current measuring circuits, the temperature sensor and the protection circuits were presented. The chapter 4 described the realization of the control program on TI’s DSP TMS320LF2407A. At the beginning, a Q-format floating point representation for the fixed-point DSP was introduced and the normalization of voltage and current in the per-unit system was developed. A digital PI 71 6.2. Conclusions control program was also presented here. At the end, the current control algorithm for different types of loads was described. The chapter 5 showed the experimental results for different loads and conditions. The current ripples of the experimental results were analyzed Application Note and matched well with the theoretical analysis in the chapter 2. PIC-004 6.2 Conclusions Test Using Electronic Load This thesis hasStart-Up explored the design of an AC electronic loads used for the unintentional islanding test and hardware systemInc. has been developed Zach Zhang, Alphaa &real Omega Semiconductor, and tested with a 60Hz AC voltage source. According to the experimental 1. Introduction Load this system is able to simulate different types of results, it to is Electronic proved that loads and in this project it was programmed to simulate resistive, inductive Electronic loads, such as Chroma 6310 series, can be configured as constant current (CC) mode, constant voltage (CV) mode constant resistance (CR) mode. In CC mode, the load will sink a current in accordance the programmed andand capacitive loads. In order to simulate a complex circuit,withsome minor value regardless of input voltage. In CV mode, the load will sink current to control the voltage source in programmed value. modifications of the current command algorithm areprogrammed required. In CR mode, the load sinks a current linearly proportional to thecalculation input voltage in accordance with the resistance. Internalthis feedback control could circuit controls loadacurrent based on its mode set up and input voltage.loads Therefore, device workthe as programmable AC electronic the real unintentional islanding Afor simplified electronics load diagram is shown in Figuretest. 1. Use CC mode as example, Iset is the load current set value. Isense is the actual load current sense signal. The feedback consists errorH-bridge amplifier and inverter feedback network f. The However, because this system circuit is based onof an andZworks amplifier compares the Iset and Isense, the output signal is used to control power transistor Q1 base voltage. As a result, the in the switching mode, current ripple is much higher that of the power transistor Q1 only pulls current the Iload equal to Iset, into the electronic load. When Iload < Ithan set, feedback circuit will turn on Q1 harder to pull more current until I = I . When I > I , feedback circuit will turn on Q1 less to pull less current until load set load set traditional electronic load which works in the linear mode. Iload = Iset. Figure 1: Simplified electronics load diagram Figure 6.1: Simplified Electronic Load Diagram Besides operation modes, there are other parameters that need to be programmed before each use. They are Turn-on Voltage (Von), Von Latch switch and Slew Rate, etc. The electronic load starts to sink current when the input voltage reaches Von voltage. Enabled Von latch means the load will sink current continuously when the input voltage reaches Von voltage. Disabled 6.1 Von latch means the the load will stop sinking current whenof its input is below Von voltage level. The Figure shows simplified structure thevoltage traditional electronic slew rate is defined as current change over time. Regardless the input voltage, the load will sink current at programmed slew load. Iset is the load current set value. Isense is the actual load current rate. August 08 Tel: 408.830.9742 • Fax: 408.830.9749 • www.aosmd.com 72 1 6.3. Suggestions for Future Work sense signal. The feedback circuit consists of error amplifier and feedback network Zf . The amplifier compares the Iset and Isense , the output signal is used to control power transistor Q1 base voltage. As a result, the power transistor Q1 only pulls current Iload equal to Iset , into the electronic load. When Iload < Iset , feedback circuit will turn on Q1 harder to pull more current until Iload = Iset . When Iload > Iset , feedback circuit will turn on Q1 less to pull less current until Iload = Iset . Therefore, a combination of the switching mode and the linear mode might be a better solution for the higher accuracy requirement (shown in Figure 6.2). The switching mode block could track the required load characteristics and the linear mode block can compensate the current ripple induced by the switching operation to achieve the smoother load characteristics. Iload Switching Mode AC Linear Mode Figure 6.2: Combination Electronic Load Diagram 6.3 Suggestions for Future Work It has been proved that the system designed in this thesis is able to work as an electronic load. But in order to finish the real unintentional islanding test, there is still a lot of work to do. • A DC power supply which can output up to 220V is required for the real islanding test. As the previous analysis, the system need a 73 6.3. Suggestions for Future Work regulated DC bus voltage. Therefore, the islanding test needs a high voltage DC power supply to provide a high enough DC bus voltage. • A three-phase four-leg inverter is required for the three-phase unintentional islanding test. In our design, the system is an H-bridge inverter which is only able to do the single phase islanding test. For the threephase test, a four-leg inverter should be used and the schematic is shown in the Appendix B. • The output filter has to be designed for the real test. Due to the switching parts are used in the circuit, a lot of high frequency interference is introduced. But in the real islanding test, too much noise will lead to the testing failure. Therefore, an output filter is needed to reduce the high frequency noise. (Shown in the Appendix B) • A phase-locked loop (PLL) circuit is required for the frequency detection. In our design, the frequency detection was implemented in the DSP program. However, it is not very accurate and unable to perform a real-time frequency detection. The PLL circuit can detect the real-time frequency of the applied voltage source and provide a high precision. • Modify the current calculation algorithm for the islanding test. The islanding test requires a variable LC circuit with respect to the different output reactive power generate by the device under test. So a new current calculation algorithm has to be implemented. 74 Bibliography [1] Niklas Strath. Islanding Detection in Power Systems. Lund University, 2005. [2] Konrad Mauch Wilsun Xu and Sylvain Martel. An Assessment of Distributed Generation Islanding Detection Methods and Issues for Canada. CETC-Varennes, 2004. [3] Robert W. Erickson and Dragan Maksimovic. Fundamentals of Power Electronics. Springer, 2000. [4] MUBW20-06A7 Datasheet. IXYS Corporation, 2004. [5] Six-pack SCALE Driver 6SD106E for IGBTs and Power MOSFETs. CONCEPT: www.IGBT-Driver.com, 2000. [6] Alex Tessarolo. SPRA989A: F2810, F2811, and F2812 ADC Calibration. Texas Instruments, 2004. [7] TMS320LF/LC240xA DSP Controllers Reference Guide: System and Peripherals. Texas Instruments, 2006. [8] TMS320F/C24x DSP Controllers Reference Guide: CPU and Instruction Set. Texas Instruments, 1999. [9] TMS320C1x/C2x/C2xx/C5x Assembly Language Tools: User’s Guide. Texas Instruments, 1995. [10] David M. Alter. Getting Started in C and Assembly Code With the TMS320LF240x DSP. Texas Instruments, 2002. [11] Erick L. Oberstar. Fixed-Point Representation and Fractional Math. Oberstar Consulting, 2007. 75 Appendix A Schematics A.1 Schematics of Control Board 76 C24 10uF C10 10uF C23 10uF C9 10uF C22 10uF C8 10uF C21 10uF C7 10uF C20 10uF +15 C6 10uF +15 +15 1 2 3 4 CN1 L1 3.5uH L2 3.5uH C26 10uF C25 10uF C12 10uF C11 10uF C19 10uF C5 10uF C18 10uF C4 10uF C17 10uF C3 10uF 3 2 1 3 2 1 CC10-1205SF-E -Vin RC +Vin CC2 CC10-1212DF-E -Vin RC +Vin CC1 NC -Vout Trim +Vout -Vout Com Trim +Vout 4 5 6 7 4 5 6 7 C2 C1 C15 0.1uF C16 C13 0.1uF 0.1uF 10uF C14 10uF 10uF 5V 16V 16V D3 D2 D1 +5 -15i +15i A.1. Schematics of Control Board Figure A.1: PowerSupply AD.SCHDOC 77 A.1. Schematics of Control Board -15i 3 2 1 R16 Res1 240 R17 R18 Res1 30K -15i C32 30K +15i 150pF C33 R19 100nF 3.3K 2 1 1 LM324AM U4A 11 R21 12K 5 C34 100nF 2 9 6 7 2 10 LM324AM U4B 11 3 +15i +15I 4 4 5.1K R20 4 +15i 8 3 VM1 LM324AM U3C 11 CN7 -15i -15I U5 -15i LM4040A30IDBZR-3.0 -15i R22 Res1 240 CN8 +15i R23 R24 C35 30K 30K 150pF R25 8 3 R26 12K 13 13 12 -15I 14 4 12 LM324AM U4D 11 11 LM324AM U4C 14 4 VM2 LM324AM U3D 11 9 10 +15i +15I 4 4 3.3K 4 +15I -15i -15I -15i 3 2 1 R27 Res1 240 CN9 +15i R28 R29 C36 30K 150pF +15i +15i 4 +15I 1 1 LM324AM U6A C38 R31 12K 9 6 5 7 2 LM324AM U6B 10 8 3 VM3 LM324AM U6C 11 2 3 R30 3.3K 100nF 4 4 C37 11 30K 11 1 3 NC 3 2 1 -15i -15I 100nF -15i Figure A.2: VoltageSensor.SCHDOC 78 A.1. Schematics of Control Board -15i 3 2 1 R1 Res1 62 CN2 +15i R2 R3 Res1 32K C27 30K 150pF +15i C28 4 4 LM324AM U1A C29 11 2 1 1 R6 12K 2 6 5 U2 LM4040A30IDBZR-3.0 7 2 11 2 3 +15i 3.48K +15I 100nF 4 R5 5.1K R4 3 LM324AM U1B -15I 1 1 11 -15i CM1 LM324AM U3A -15i NC 100nF 1 3 -15i -15i 3 2 1 R7 Res1 62 CN3 +15i 30K C30 150pF R10 3.48K +15I 8 R11 LM324AM U1C 12K 11 3 -15I 6 13 12 14 4 LM324AM U1D 11 9 10 +15i 4 4 +15I 4 R8 Res1 32K 5 7 2 CM2 LM324AM U3B 11 R9 -15i -15I Figure A.3: CurrentSensor.SCHDOC 79 2.5v D4 R14 1.5K +5 300 R15 TMP9 4 5 5K R12 12 +5 3 2 100nF C31 LM339AM U3A R13 3K TMP 4 3 2 1 2 1 4 3 2 1 2 4 6 8 10 12 14 16 18 20 Header 4 CN6 JP1 1 3 5 7 9 11 13 15 17 19 CN5 CN4 VM3 VM2 VM1 CM2 CM1 SW FRQCN FRQCN +5 SW TMP9 TMP +15 +5 A.1. Schematics of Control Board Figure A.4: TMPSensor.SCHDOC 80 4 3 2 1 40 38 36 34 32 30 28 26 24 22 20 18 16 14 12 10 8 6 4 2 JP1 39 37 35 33 31 29 27 25 23 21 19 17 15 13 11 9 7 5 3 1 6 7 2 1 5 4 14 0.1uF C2 SN74ALS21AD U1A 10uF C1 6 +5 1 TMP 4V7 4V7 4V7 SN74LVC2GU04 U2A D1 16V +15 +5 D4 D3 D2 2 1 +5 +15 +15 +15 +5 +15 +15 +15 +5 +15 +15 Header 2 P2 4.7K 4.7K R12 22K R11 22K R10 4.7K R9 C8150pF C7150pF C6 150pF 22K R722K R84.7K C5 150pF 150pF 150pF C4 C3 R6 R1 4.7K R5 4.7K 22K R4 22K R3 R2 TMP9 TMP +15 34 33 32 31 30 29 28 27 26 25 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 6SD106E GND W_SO1 W_VL W_RC1 W_InA W_InB W_RC2 W_Mod W_SO2 VDD VDD GND GND V_SO1 V_VL V_RC1 V_InA V_InB V_RC2 V_Mod V_SO2 VDC VDC GND GND U_SO1 U_VL U_RC1 U_InA U_InB U_RC2 U_Mod U_SO2 GND Concept1 U_C2 U_Rth2 U_E2 U_G2 Free Free U_C1 U_Rth1 U_E1 U_G1 Free Free V_C2 V_Rth2 V_E2 V_G2 Free Free V_C1 V_Rth1 V_E1 V_G1 Free Free W_C2 W_Rth2 W_E2 W_G2 Free Free W_C1 W_Rth1 W_E1 W_G1 68 67 66 65 64 63 62 61 60 59 58 57 56 55 54 53 52 51 50 49 48 47 46 45 44 43 42 41 40 39 38 37 36 35 D6 180 R24 R23 R22 R21 R20 R19 R18 R17 180 39K D10 180 39K D9 180 39K D8 180 39K D7 39K R16 180 R15 R14 D5 R13 39K WG1 WE1 WC1 WG2 WE2 WC2 VG1 VE1 VC1 VG2 VE2 VC2 UG1 UE1 UC1 UG2 UE2 UC2 15V 15V D22 15V D21 15V D20 15V D19 15V D18 15V D17 15V D16 15V D15 15V D14 15V D13 15V D12 D11 A.2 Header 4 P1 A.2. Schematics of Power Board Schematics of Power Board Figure A.5: Driver.SCHDOC 81 A.2. Schematics of Power Board C9 * 2700uF C10 R25 * Thermistor 5 1.5uF C11 * 1.5uF UC2 UE2 UG2 UC1 UE1 UG1 VC2 VE2 VG2 VC1 VE1 VG1 Inv1 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 3 2 1 24 23 22 21 20 19 18 17 16 15 14 13 12 11 10 3 2 1 4 5 6 7 8 9 4 5 6 7 8 9 1 2 3 JP3 +5 TMP9 MUBW20-06A7 WC2 WE2 WG2 WC1 JP2 1 2 3 WE1 WG1 Figure A.6: MUBW20.SCHDOC 82 A.3. PCB Layouts A.3 PCB Layouts Figure A.7: Control Board PCB 83 A.3. PCB Layouts Figure A.8: Control Board PCB Top Layer 84 A.3. PCB Layouts Figure A.9: Control Board PCB Bottom Layer 85 A.3. PCB Layouts Figure A.10: Power Board PCB 86 A.3. PCB Layouts Figure A.11: Power Board PCB Top Layer 87 A.3. PCB Layouts Figure A.12: Power Board PCB Bottom Layer 88 Appendix B Schematics for 3-Phase Islanding Test 89 Appendix B. Schematics for 3-Phase Islanding Test Figure B.1: Schematics for 3-phase Islanding Test 90 Appendix C Unintentional Islanding Test Conditions Unintentional Islanding Test Test conditions: Single Phase 120V /60H Z /1kW 1. PF=[ 1, 0.37, 0.707 ] ⇔ Q f =[ 0, 2.5, 1 ] 2. Values Table: Criteria Value Unit Pload 1000 W Qload 0 VAR Matched LC VEPS 120 V RMS fEPS 60 Hz VIUT 120 V PIUT 1000 W IIUT 8.33 A PFIUT 0.95 Qf 1 Rload 14.4 OHM LLoad 38.197 mH CLoad 184.207 F Qf 2.5 Rload 14.4 OHM LLoad 15.28 mH R / 2 f o Q f , iL(0)=11.785A CLoad 460.52 F Q f / 2 f o R Q f =R Notes Inverter Under Test RMS =18.195o RLC below designed for this Qf V2/P R / 2 f o Q f , iL(0)=11.785A Q f / 2 f o R RLC below designed for this Qf V2/P 1 C for a parallel RLC load, and Q F = −1 . L PF 2 P IUT =P LoadP EPS & QIUT =Q Load QEPS 91
- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- A DSP based AC electronic load for unintentional islanding...
Open Collections
UBC Theses and Dissertations
Featured Collection
UBC Theses and Dissertations
A DSP based AC electronic load for unintentional islanding tests Feng, Wei 2009
pdf
Page Metadata
Item Metadata
Title | A DSP based AC electronic load for unintentional islanding tests |
Creator |
Feng, Wei |
Publisher | University of British Columbia |
Date Issued | 2009 |
Description | The distributed generation (DG) system for the renewable energy resources is being used more and more widely around the world in the past decades. In order to protect the whole power system from the islanding situation, the interface inverters of the DG system have to pass the unintentional islanding test. In this study, a single phase AC electronic load is introduced to simulate the local loads under unintentional islanding conditions. This thesis proposes a proper schematic of the electronic loads based on the H-bridge dc-dc converter and analyzes the performance of this system with PowerSIM. Then, a set of specifications for different ranges of the electronic loads are calculated through theoretical formulas and the PSIM simulations. Meanwhile, the transfer function of the system is also derived to analyze the stability of the PI control system. According to the theoretical analysis and simulations, a control program is implemented based on the Texas Instruments (TI) Digital Signal Processor (DSP) TMS320F2407A for the different kinds of electronic loads. A test circuit is then built to validate the performance of the system. Some experiments are performed for the resitive, inductive and capacitive loads respectively. |
Extent | 5941123 bytes |
Genre |
Thesis/Dissertation |
Type |
Text |
File Format | application/pdf |
Language | eng |
Date Available | 2009-08-11 |
Provider | Vancouver : University of British Columbia Library |
Rights | Attribution-NonCommercial-NoDerivatives 4.0 International |
DOI | 10.14288/1.0067479 |
URI | http://hdl.handle.net/2429/11991 |
Degree |
Master of Applied Science - MASc |
Program |
Electrical and Computer Engineering |
Affiliation |
Applied Science, Faculty of Electrical and Computer Engineering, Department of |
Degree Grantor | University of British Columbia |
Graduation Date | 2009-11 |
Campus |
UBCV |
Scholarly Level | Graduate |
Rights URI | http://creativecommons.org/licenses/by-nc-nd/4.0/ |
Aggregated Source Repository | DSpace |
Download
- Media
- 24-ubc_2009_fall_feng_wei.pdf [ 5.67MB ]
- Metadata
- JSON: 24-1.0067479.json
- JSON-LD: 24-1.0067479-ld.json
- RDF/XML (Pretty): 24-1.0067479-rdf.xml
- RDF/JSON: 24-1.0067479-rdf.json
- Turtle: 24-1.0067479-turtle.txt
- N-Triples: 24-1.0067479-rdf-ntriples.txt
- Original Record: 24-1.0067479-source.json
- Full Text
- 24-1.0067479-fulltext.txt
- Citation
- 24-1.0067479.ris
Full Text
Cite
Citation Scheme:
Usage Statistics
Share
Embed
Customize your widget with the following options, then copy and paste the code below into the HTML
of your page to embed this item in your website.
<div id="ubcOpenCollectionsWidgetDisplay">
<script id="ubcOpenCollectionsWidget"
src="{[{embed.src}]}"
data-item="{[{embed.item}]}"
data-collection="{[{embed.collection}]}"
data-metadata="{[{embed.showMetadata}]}"
data-width="{[{embed.width}]}"
async >
</script>
</div>
Our image viewer uses the IIIF 2.0 standard.
To load this item in other compatible viewers, use this url:
http://iiif.library.ubc.ca/presentation/dsp.24.1-0067479/manifest