Flicker en harmonischen
Inhoud
1. Flicker: meten volgens IEC 61000-3-3 en 61000-3-11 2. Flickermeter volgens IEC 61000-4-15 3. Implementatie van een flickermeter in Simulink 4. Harmonischen: meten volgens IEC 61000-3-2 en 61000-3-12 5. Harmonischen: praktische metingen
1
1. Flicker: meten volgens IEC 61000-3-3 en 61000-3-11 1.1. Proefopstelling
R
X
RT
XT
R
X
RT
XT
R
X
RT
XT
~ M 3~
~ ~
Fasespanning
PC
Flickermeter
Oscilloscoop
Lijnstroom
PC
1.1. Proefopstelling
2
1.2. Resultaten flickermeter
Pst number Pst value 1 2 3 4 5 6 7 8 9 10 11 12
Berekening PLT door flickermeter 3,291 3,24 5,248 3,196 0,309 5,228 0,39 0,389 0,345 0,354 0,203
3.291 3.24 5.248 3.196 0.309 5.228 0.39 0.389 0.345 0.354 0.203 0.239
35,64377 34,01222 144,5378 32,64527 0,029504 142,8916 0,059319 0,058864 0,041064 0,044362 0,008365
0,239 0,013652 Plt
3,191213
12
PLT =
Plt = 3.191
3
∑ i =1
PSTi 12
3
1.2. Resultaten oscilloscoop Spanning
Stroom
3
1.2. Resultaten oscilloscoop
Niet correcte eerste waarde, doordat de start niet op een nuldoorgang ligt. U(t) : 268.0792 241.4294 242.6910 225.0022 190.8998 187.2540 199.5139 218.7647 227.0236 235.6336 243.4847 239.6999 237.7963 238.4744 242.0083 240.2081 239.2281 238.6521 241.1892 240.0378
Tf
Tt
10 ms
1.2. Resultaten oscilloscoop
Tf = 20 ms Tt = 50 ms F = 0.94 F = 0.94
4
1.2. Resultaten oscilloscoop
2 ms, 36°
1.2. Resultaten oscilloscoop 1 2 3 4 5 6 7
Tf 20 10 10 10 10 20 10
Tt 50 60 60 60 60 40 60
F 0,94 0,94 0,94 0,94 0,94 0,94 0,94
∆ U (max) ∆ U (max) top 77,78 54,99876544 76,79 54,29872973 79,33 56,09478095 86,26 60,99503095 85,418 60,39964704 85 60,1040764 87,39 61,79406161
Un 242 242 241,5 241,5 242,5 244 244,3
d (max) 22,726763 22,437492 23,227653 25,256742 24,907071 24,632818 25,294335
tf 41367,30302 39705,87213 44356,16557 57988,83641 55458,67283 53528,13257 58265,49576 Plt form 1 Plt form 2 flick
(
t f = 2 . 3 F .d max P ST = 2 . 3
∑ T
t p
f
)
3 .2
1 3 .2
PLT
Tp = 600s
PST PST tot Gemeten Pst 3,754335773 3,29 3,706549807 3,24 3,837080035 4,172274346 4,982801 5,25 4,114510684 3,2 4,069205644 4,178484651 5,122218 5,23 3,367941944 3,300913188
∑ tf = 2 .3 T p 12
PLT =
3
∑ i =1
3,191
1
3 .2 Tp = 7200s
PSTi 12
3
5
2. Flickermeter volgens IEC 61000-4-15 2.1. Architectuur Architectuur kan in twee delen verdeeld worden: -simulatie van de “lamp-eye-brain chain respons” -on-line statistische analyse van het flickersignaal en de output van de resultaten Implementatie bestaat uit 5 blokken: -blok 1: input and calibration -blok 2, 3 and 4: chain response -blok 5: statistical analysis
2.1. Architectuur - let op het feit dat de meter als een analoge meter beschreven wordt, met een digitalisering in blok 5
6
2.1. Architectuur Blok 1: •signaalgenerator voor calibratie •input voltage adaptor: flicker wordt relatief uitgedrukt (verhouding) schalen van de input heeft geen invloed
Blok 2: •kwadrateren: simulatie van het gedrag van een lamp
2.1. Architectuur Blok 3: •Weighting filters: eerste filter: eliminatie van DC en netharmonischen tweede filter: weighting filter voor simulatie van de “lamp-eye chain” •Range selector: zie verder: van belang voor de analoge meter Transfer functie voor de weighting filters:
7
2.1. Architectuur Blok 4: •Squaring multiplier: Eye-brain chain •Eerste orde sliding mean filter Opslageffect in de hersenen Tijdsconstante van 300 ms Blok 5: •On-line statistische analyse: Sampling of the flicker level Minimum sample rate is 50 s-1 Verdeel flicker level in klassen (bins) Tel hoeveel keer het flicker niveau voorkomt op een zeker niveau (frequency distribution function) Bereken de cumulative probability function in the time-at-level statistical method
2.1. Architectuur Blok 5: voorbeeld met 10 klassen
8
2.2. Evaluatie Short-term flicker evaluation: * Observatie tijd Tst = 10 min * Formule:
de percentielen Pi zijn de flicker levels die overschreden worden gedurende i% van de tijd. Pis zijn de uitgemiddelde waarden, bijv.: Long-term flicker evaluation: * combinatie van N short-term periodes:
2.3. Performance testing Betekenis van ∆V/V: 50 Hz, 1V gemiddelde spanning, ∆V/V van 40%, 8.8 Hz rechthoekige modulatie:
9
2.3. Performance testing De volgende waarden zouden een Pst = 1.00 ± 0.05 moeten geven:
waarde x 3 → Pst = 3
3. Implementatie van een flickermeter in Simulink 3.1. Matlab & Simulink Matlab: • MATrix LABoratory • high-performance language for technical computing • integrates: •mathematics • visualisation • programmation → M-files
10
3.1. Matlab & Simulink Simulink:
• software package for modelling, simulating and analysing • building models as block diagrams • linked to Matlab by the workspace (variables, postprocessing, …)
3.1. Matlab & Simulink Simulink building blocks: • basic building blocks: sources, sinks, functions, subsystems
• user functions: ~ M-files ~ S-function M-files ~ C-MEX S-functions (compiled C-coded S-functions)
• toolboxes and third-party libraries: ~ dSpace RTI1102
11
3.2. Model
Blok testsignaal:
•simulatie: input via generator •testsignaal volgens norm (Pst=1) •werkelijke meting: toevoer via datafile
12
Blok 1: voltage adapter en sliding filter
•voltage adapter: verwachte input van het model is 1Vrms •berekenen van RMS-waarde (M-function) •sliding filter voor het wegfilteren van trage spanningsvariaties •nadeel: opstart transiënt (Simulink ziet dit de start als een aangelegde stap) •de opstart transiënt wordt als flicker beschouwd oplossing: * transiënt verminderen door andere startwaarde * eerste x aantal seconden weglaten
Blok 1: voltage adapter en sliding filter •welke startwaarde: * RMS-waarde is altijd ongeveer 1V (andere waarden zijn aanpasbaar via de voltage adapter) * test gebeurde met 7.07V en 9.9V als initialisatiewaarde voor sftijd10s Trunc van eerste x seconden
Init = 0
Init = 10/sqrt(2)
Init = 9.99
Evalflicker1(0)
0.4350
2.1904
2.0834
Evalflicker1(10)
3.3324
2.3797
1.5251
Evalflicker1(20)
1.0863
1.6114
1.3131
Evalflicker1(30)
2.6046
1.3801
1.2666
Evalflicker1(40)
1.0000
1.3529
1.2666
Evalflicker1(50)
2.8038
1.3306
1.2666
Evalflicker1(60)
2.1172
1.3216
1.2666
Evalflicker1(70)
2.5160
1.3002
1.2666
Evalflicker1(80)
2.1341
1.2817
1.2666
Evalflicker1(90)
1.1808
1.2666
1.2666
Te oscillerend St.st vanaf 90s St.st. vanaf 30 s
•gebruikte methode: 9.9 neem de eerste 40s als ongeldige meetpunten, dus 10 min + 40s meten
13
Blok 2: squaring operator
•RMS-waarde en vermogen van de lamp zijn gerelateerd: U²/R
Blok 3: Weighting filters
14
Blok 4: Sliding mean filter en versterker voor calibratie
•averaging filter: simulatie van het “opslag” effect in het menselijk geheugen •versterker: ~ calibratie ~ testsignalen moeten Pst = 1 geven
Blok 5: Statistische analyse
•sampler: 100 Hz sampling •sfhistogram: S-function voor het berekenen van een histogram •memory en feedback: cumulatie per bin
15
Blok 6: Data: •blok voor het oplossen van Matlab-problemen •probleem: Simulink onthoudt alles (400 bins, gesampled op 1kHz, 10 min) Beperken van geheugenverbruik: - eerste versie:
Matlab loopt tegen 640 MB geheugengebruik Oplossing: alles schrijven naar file i.p.v. variabelen
- tweede versie:
Bestand is 1 GB en kan niet door Matlab geopend worden Oplossing: schrijven enabled pas in de laatste seconde Probleem, we hebben de 40ste seconde ook nodig. Daarom volgend systeem om ook de 40ste seconde op te slaan
Naverwerking van de data: Evalueren met evalflicker.m Resultaat is de CPF en de Pst-waarde
16
Handleiding: -laat het model lopen -laat evalflicker lopen -kijk naar fig. 1 -als de waarde op het maximaal aantal bins (bins is bijvoorbeeld 6400) niet 0 is, dan moet je max verhogen (maximum hiervoor is 2200 voor 20% variatie) - let op:bij max*2 moet aantal bins ook *2 voor dezelfde calibratie!
Invloed van samplefrequentie (Simulink) op het resultaat:
Basis simulatiestap voor simulaties: - testsignaal is 50 Hz en 13.5 Hz modulatie - gemeenschappelijk veelvoud is 337.5 Hz of 3 ms Test met volgende frequenties: - 200 Hz geeft Pst = 1.0368 - 1 kHz geeft Pst = 0.9785 - 2 kHz geeft Pst = 0.9791 - 337.5 Hz geeft Pst = 0.9797 - 650 Hz geeft Pst = 1.0237 Voor echte metingen: - samplefrequentie is 2 x gewenste maximale flickerfrequentie - norm spreekt over 33.3 Hz, dit geeft 66.6 Hz - Pst = 0.9998 - intern zit echter een sampler op 100 Hz, hiervoor nemen we een simulatiestap van 5e-3 of 200 Hz
17
3.3. Performance testing Testwaarden (norm):
- juiste waarde is Pst = 1
(versterkingswaarde 3e5) - maximale fout moet < 5% (gemeten < 3%) - simulatiestap 5e-3 Variaties/min
DV/V %
Pst juist
Pst gemeten
bins
maxhistinput
1 square
2.724
1
1.0282
1200
30
2 square
2.211
1
1.0066
800
20
7 square
1.459
1
1.0119
400
10
39 square
0.906
1
1.0189
400
10
110 square
0.725
1
1.0156
400
10
1620 square
0.402
1
0.9951
400
10
1620 square
1.206
3
2.9960
800
20
8.8 sine
20
(max. vereiste meetwaarde)
57
6400
6555
3.3. Performance testing Probleem: * niet altijd real time: * afhankelijk van de computersnelheid * afhankelijk van sampletijd (bijv. op 3.2 GHz) 200 Hz: rekentijd 64% van real time 1 kHz: rekentijd 153% van real time * voorlopige oplossing: 10 min. loggen, daarna verwerken * volgende stap: m-files (interpreter) omzetten naar gecompileerde C-files
18
4. Harmonischen: meten volgens IEC 61000-3-2 en 61000-3-12 4.1. IEC 61000-3-2 en IEC 61000-3-12 •
harmonic currents
•
scope: 61000-3-2: equipment having an input current ≤ 16A per phase (low voltage public network; no limits for network <220 Vphase) 61000-3-12: equipment having an input current > 16A and ≤ 75A per phase (max. 240V single phase, 690V three phase)
•
important definitions: Total harmonic current: total rms value of the harmonic components of orders 2 to 40: 40
∑ I n2 n=2
4.2. IEC 61000-3-2 •
Classification: Class A: - balanced three phase equipment - household appliances - tools (non-portable) - dimmers - audio Class B: - portable tools - non-prof arc welding equipment Class C: - lighting equipment Class D: - PC, monitor, TV
19
4.3. Limits 61000-3-2: Class A: table 1 Class B: table 1 x 1.5 Class C > 25 W: table 2 Class C < 25 W: table 3 column 2 OR 3de harm < 86% and 5th harm < 61% Class D: table 3
4.4. Meten Voedingsspanning: stabiel beperkte harmonische vervorming:
Testcircuit:
Enkelfasige en driefasig Zs en Zm zijn niet gedefinieerd, maar moeten voldoende laag zijn
20
4.4. Meten Meetprincipe: - voor elke harmonische: 1.5 s smoothed rms current bepalen - observatietijd: voldoende voor een herhaalbaarheid van 5% - meetkundig gemiddelde van de gemeten waarden Toepassen van de limieten: - individuele harmonische, meetkundig gemiddelde ≤ limiet - voor elke harmonische orde, elke 1,5 s smoothed value ≤ 150% limiet (200% onder bepaalde voorwaarden) Smoothed value: - beschrijving in IEC 61000-4-7 - digitale LPF met tijdsconstante 1.5 s
4.5. IEC 61000-3-12 Belangrijke definities: Three phase short circuit power Ssc = U2nom / Z Rated apparent power:
Short circuit ratio:
21
4.6. Limieten:
4.7. Meten Meetprincipe: - voor elke harmonische: 1,5 s smoothed rms current bepalen - observatietijd: voldoende voor een herhaalbaarheid van 5% (10% boven 7de orde, 1% voor grondgolf) - meetkundig gemiddelde van de gemeten waarden Toepassen van de limieten: - individuele harmonische, meetkundig gemiddelde ≤ limiet - voor elke harmonische orde, elke 1,5 s smoothed value ≤ 150% limiet Smoothed value: - beschrijving in IEC 61000-4-7: * per harmonische, de RMS sturen door een digitale LPF met tijdsconstante 1.5 s
22
5. Harmonischen: praktische metingen 5.1. Oscilloscoop met fft-functie
Principe: neem een geheel aantal periodes en voor een fft uit let op:
- oscilloscoop geeft amplitude en geen RMS - Micro VIP geeft RMS
5. Harmonischen: praktische metingen 5.1. Oscilloscoop met fft-functie Tijdsbasis: * voor onze oscilloscoop: - 40 000 punten per reading - mogelijk aantal sec/div: 1s, 2s, 5s, 10s - aantal seconden per reading: 10s, 20s, 50s, 100s - samplefrequenties (40 000 samples per aantal seconden): 4 kHz, 2 kHz, 800 Hz, 400 Hz - 40ste orde harmonische is 2 kHz, dus volgens Shannon Fs = 4 kHz - besluit: max. 1sec/div
23
5. Harmonischen: praktische metingen 5.1. Oscilloscoop met fft-functie Probleem: - wat met de 1.5s smoothing? - voert een fft een smoothing uit?
Signaal met subharmonische. Subharmonische bevat geen geheel aantal perioden binnen de 1.5s.
-volgens de norm: •per periode de rms waarde van elke harmonische bepalen •smoothing
5. Harmonischen: praktische metingen 5.1. Oscilloscoop met fft-functie Vgl. zonder subharmonische en met subharmonische:
Zonder
Met
Grondgolf
230 V
232.2 V
3de
23 V
23.22 V
9de
11.5 V
11.60 V
Is dit hetzelfde als smoothing? Indien niet, is de fout toelaatbaar?
24
5. Harmonischen: praktische metingen 5.2. MicroVIP Driefasige energie analyzer Extentie: harmonic box 4 kHz sampling 1 meting / seconde Extentie: memory pack Loggen: 1 record/min
25
26
27
28
