最近の寒波で電力供給の危機が叫ばれています。原子力発電所が震災以来(ほぼ?)停止していて、CO2の対策として石炭火力や石油での火力もだいぶ止まっているような話を聞きます。そこへ来て、大雪によってソーラーパネルに雪が積もって発電できないのですから当たり前といえば当たり前です。
電力供給が厳しくなると発電所から見た負荷が重くなるので、電源周波数の変動として見えてくるはずです。そこで、電源周波数の変動を見てみよう、ということになります。
電源周波数の取り出し
ちゃんと見ようと思うと、どうしても商用電源へのアクセスが必要で感電防止とかを考えると大がかりになりがちですが、今回は100円ショップで売っているニッケル水素電池の充電器を使うことにしました。100円ショップで入手したものは、以前分解したときに(電池を入れない状態では)電池のマイナス端子に内部のトランスの端子がそのまま出ていることに気づきました。そこで、ここから商用電源の信号を取ることにします。(注意: 感電の恐れが大きい部分に触れることなく取り出せますが、この充電器の本来の使い方とは異なりますので、まねする場合には十分注意して行うとともに、自己責任で行ってください。似た形状で回路構成が違うものもあるかもしれませんし。)
このようにオシロスコープのプローブを2つの電池のマイナス端子に接続して波形観測すると、Peak-to-Peakで12Vくらいあることがわかります。ここにフォトカプラを入れて信号を取り出すことを考えます。
手持ちのフォトカプラPC817のデータシートを見ると、以下の特性を持っています。
フォトカプラの内蔵LEDにはIFは5mAも流してやれば良さそうです。VF=1.2V程度なので、ピーク時の6Vの際に10mA程度をターゲットに470Ωをつけてやることにします。抵抗には10mA程度流れるときに5Vの電圧がかかるので、50mWの消費電力・・・と思ったけど直流ではないので平均はもっと少ないですが、1/6Wのものを使うことにしました。あとはフォトカプラの内蔵LEDと逆向きに適当なLEDをつけておきます。これは何でもよく、AC電圧が逆相の際にフォトカプラの内蔵LEDに逆電圧がかかるのを防止するためです。
フォトカプラのフォトトランジスタ側の負荷抵抗は50Hzや60Hzでは周波数特性もへったくれもないので、IC=0.5mAも流れていれば十分ですから4.7kΩとしました。(よく考えたら、この抵抗は無しでIO8の内蔵プルアップを有効にするだけで十分かもしれない)
回路はこんな感じ。左のAC信号源の部分がニッケル水素充電器の2つのマイナス端子。真ん中のダイオードはLEDでフォトカプラの保護用。
ブレッドボード上に組んで波形を観測してみました。
右側のボードは秋月のAE-ATmega328ボードがボードのみで売られていた時代に組んだものです。ATmega328を載せてブートローダを書き込んだのだと思います。クリスタルは16MHzのものが載せてあって、Arduino Pro互換として動作していると思います。フォトカプラの出力はIO8へ接続しています。ATmega328のTCC1のカウント値をキャプチャする機能の入力信号ICP1として使用するためです。
上が充電器の2つの端子の電圧で、下がArduino互換ボードへの信号です。狙ったとおりの信号が得られています。
周期の計測
これをATmega328で周波数を計測します。TCC1の16ビットタイマをCPUのクロックを8分周したもの(つまり2MHz=0.5uS)でフリーランさせて、インプットキャプチャ機能でICP1信号(IO8)の下がりエッジでその瞬間の値をICR1レジスタに取得します。取得したら割り込みでただちにICR1レジスタの値を読み出して保存します。実際には1秒周期で平均周波数を求めてホストに送信するので、割り込み内ではlongの変数に積算するとともに、積算回数をカウントします。また、よりデータを取る場合にはその1秒間での最短周期や最長周期も取得しておきます。
メインのループ側では1秒毎に積算値から平均の周期もしくは周波数を求めてシリアルでホストに送信します。それをArduino IDEのシリアルプロッタの機能でグラフ化するとこんな感じになります。
縦軸は周波数で、単位は1/1000Hzです。50000が50Hzぴったりということになります。このグラフの例では縦軸は49.92Hzから50.16Hzまでが表示されていて、電源の周波数は49.94Hzくらいから50.1Hzくらいまで変動しているということがわかります。右側はサンプル数ですが、データ送信は1秒間に1回なので、過去500秒間の電源周波数変動が見えます。
これをブレッドボードではなくて、ユニバーサル基板で組んでみました。
Arduino用のユニバーサル基板もあったのですが、回路規模が小さいのでより小さなユニバーサル基板で組んでみました。この範囲には5Vがないので、IO12をH出力してフォトカプラの内蔵フォトトランジスタの負荷抵抗の電源端子として使用することにしました。(注:実際には、IO8の内蔵プルアップを有効にするだけでも十分そうです。頑張れば更に小さいユニバーサル基板でも入ったかも。)
さらにソフトを修正して、1秒間の間の最短周期と最大周期も送信するようにしたところ、以下のような感じになりました。
理由はきちんと確認していませんが、1秒間の中での周波数変動が大きい時間と小さい時間があるようです。おそらくはインバータ機器などが動作している/していない期間なのでしょう。
ソフトウェア
平均のみをホストに送信する場合
1 | #define TCCCLK 2000000 // CPUクロックを8分周した値。 |
3 | #define ICP1 8 // ICP1の端子を定義 |
4 | #define PWR 12 // フォトカプラの負荷抵抗に電源を供給するための端子 |
10 | volatile uint16_t PreCap,NowCap; |
17 | unsigned long current_time; |
20 | ISR(TIMER1_CAPT_vect) { |
23 | CapV = NowCap - PreCap; |
32 | pinMode(LED_BUILTIN, OUTPUT); |
46 | digitalWrite(PWR, HIGH); |
49 | TIMSK1 |= (1 << ICIE1); |
50 | current_time = millis(); |
53 | #define FREQ // 計測を周波数で行う場合 |
57 | static uint8_t skp = 3; |
59 | if (millis() - current_time >= 1000) |
63 | tSumV = SumV; tNumV = NumV; |
66 | tAveV = ( float )tSumV / ( float )tNumV; |
68 | digitalWrite(LED_BUILTIN, !digitalRead(LED_BUILTIN)); |
69 | current_time = millis(); |
74 | Serial.println(( double )TCCCLK/( double )tAveV*1000.); |
76 | Serial.println(tAveV/2); |
1秒間の中の最大/最小も送信する場合
1 | #define TCCCLK 2000000 // CPUクロックを8分周した値。 |
3 | #define ICP1 8 // ICP1の端子を定義 |
4 | #define PWR 12 // フォトカプラの負荷抵抗に電源を供給するための端子 |
7 | volatile uint16_t MaxV,MinV,NumV; |
10 | volatile uint16_t PreCap,NowCap; |
13 | uint16_t tMaxV,tMinV,tNumV; |
17 | unsigned long current_time; |
20 | ISR(TIMER1_CAPT_vect) { |
23 | CapV = NowCap - PreCap; |
25 | if (CapV > MaxV) MaxV = CapV; |
26 | if (CapV < MinV) MinV = CapV; |
34 | pinMode(LED_BUILTIN, OUTPUT); |
48 | digitalWrite(PWR, HIGH); |
51 | TIMSK1 |= (1 << ICIE1); |
52 | current_time = millis(); |
55 | #define FREQ // 計測を周波数で行う場合 |
59 | static uint8_t skp = 3; |
61 | if (millis() - current_time >= 1000) |
65 | tMaxV = MaxV; tMinV = MinV; tSumV = SumV; tNumV = NumV; |
66 | MaxV = 0; MinV = 65535; SumV = 0; NumV = 0; |
68 | tAveV = ( float )tSumV / ( float )tNumV; |
70 | digitalWrite(LED_BUILTIN, !digitalRead(LED_BUILTIN)); |
71 | current_time = millis(); |
76 | Serial.print(( double )TCCCLK/( double )tMaxV*1000.); |
78 | Serial.print(( double )TCCCLK/( double )tMinV*1000.); |
80 | Serial.println(( double )TCCCLK/( double )tAveV*1000.); |
82 | Serial.print(tMaxV/2); |
84 | Serial.print(tMinV/2); |
86 | Serial.println(tAveV/2); |
おまけ
実は、AC電源から信号を取ってこなくてもAC電源の周波数変動は大雑把には計測できます。(バッテリで動作しているノートPCの場合は無理かも)
ICP1端子に適当なケーブルをつないで、ICP1端子をプルアップなしの入力端子にしてやると、周りとの静電誘導にて電源周波数のノイズが回り込みます。これで一応電源周波数のモニターができます。
ただし、外部からのノイズに非常に弱いと思われますし、静電気などで破壊されるリスクも増えると思います。
1 | #define TCCCLK 2000000 // CPUクロックを8分周した値。 |
3 | #define ICP1 8 // ICP1の端子を定義 |
9 | volatile uint16_t PreCap,NowCap; |
16 | unsigned long current_time; |
19 | ISR(TIMER1_CAPT_vect) { |
22 | CapV = NowCap - PreCap; |
31 | pinMode(LED_BUILTIN, OUTPUT); |
41 | TIMSK1 |= (1 << ICIE1); |
42 | current_time = millis(); |
45 | #define FREQ // 計測を周波数で行う場合 |
49 | static uint8_t skp = 3; |
51 | if (millis() - current_time >= 1000) |
55 | tSumV = SumV; tNumV = NumV; |
58 | tAveV = ( float )tSumV / ( float )tNumV; |
60 | digitalWrite(LED_BUILTIN, !digitalRead(LED_BUILTIN)); |
61 | current_time = millis(); |
66 | Serial.println(( double )TCCCLK/( double )tAveV*1000.); |
68 | Serial.println(tAveV/2); |
周波数観測サイトのWEB側のプログラムの解説記事を作っていただけないでしょうか??
とりあえず作ってみました。
https://tomono.tokyo/2021/03/08/9373/
chart.js を使って表示するHTMLを1秒毎に吐いているだけで、実は大したことはしてないのですけどね。