このページの翻訳は最新ではありません。ここをクリックして、英語の最新版を参照してください。
三相整流器の高調波解析
この例では、Simscape™ のログ データを解析して高調波の振幅を取得する関数、全高調波歪みの割合を計算する関数、高調波の振幅をプロットする関数の使用方法を説明します。この解析対象のモデルは、三相整流器です。説明する関数は次のとおりです。
ee_getHarmonics
ee_calculateThdPercent
ee_plotHarmonics
モデルを開く
モデルを開きます。
open_system( 'ee_harmonics_rectifier' );
パラメーターの指定
パラメーターは次のとおりです。
定格 AC 電圧:
定格 AC 周波数:
AC 有効電力負荷:
DC 有効電力負荷:
合計皮相電力:
テスト回路では、AC 負荷は 5 MW を消費するように設定され、DC 負荷は約 5 MW を消費するように設定されます。
V_Rated = 4160; F_Rated = 60; P_AC = 5e6; P_DC = 5e6; S_Rated = P_AC + P_DC;
電源インピーダンスの計算
パラメーターは次のとおりです。
X/R 比:
pu インピーダンス:
pu ベース インピーダンス:
pu ベース インダクタンス:
電源直列抵抗:
電源直列インダクタンス:
XR = 15; Z_pu = 0.01; Z_base = ((V_Rated/sqrt(3))^2)/(S_Rated/3); L_base = Z_base/(2*pi*F_Rated); R_series = cos(atan(XR))*Z_pu*Z_base; L_series = sin(atan(XR))*Z_pu*L_base;
DC 抵抗の計算
パラメーターは次のとおりです。
損失を無視して計算した DC 平均電圧:
DC 抵抗:
(式の導出は、付録を参照)
V_DC = 3*sqrt(2)*V_Rated/pi; R_DC = V_DC^2/P_DC; disp( [ 'DC resistance required to draw ', num2str( P_DC ), ' W on DC side = ', num2str( R_DC ), ' Ohm' ] );
DC resistance required to draw 5000000 W on DC side = 6.3123 Ohm
Simscape モデルでのパラメーターの使用
これで、計算したパラメーターを Simscape モデル ee_harmonics_rectifier で使用できます。シミュレーションを実行すると、Simscape のログ変数 simlog_ee_harmonics_rectifier を作成するようにモデルが設定されます。
sim( 'ee_harmonics_rectifier' );
Voltage_Source_Currents = simlog_ee_harmonics_rectifier.Voltage_Source.I;
高調波データの取得
関数 ee_getHarmonics を使用して、Simscape のログ変数から高調波の次数、高調波の振幅、基本波の周波数の詳細を取得できます。
[ harmonicOrder, harmonicMagnitude, fundamentalFrequency ] = ee_getHarmonics( Voltage_Source_Currents );
基本波のピーク値の計算
基本波のピーク値を抽出できます。
fundamentalPeak = harmonicMagnitude( harmonicOrder==1 ); disp( [ 'Peak value of fundamental = ', num2str( fundamentalPeak ), ' A' ] );
Peak value of fundamental = 1945.806 A
小さい高調波の削除
基本波の 1000 分の 1 より大きい高調波を検出して保持します。
threshold = fundamentalPeak ./ 1e3; aboveThresold = harmonicMagnitude > threshold; harmonicOrder = harmonicOrder( aboveThresold )'; harmonicMagnitude = harmonicMagnitude( aboveThresold )';
表形式の高調波データの表示
高調波データは、MATLAB® の table に格納することができます。
harmonicRms = harmonicMagnitude./sqrt(2); harmonicPercentage = 100.*harmonicMagnitude./harmonicMagnitude( harmonicOrder==1 ); harmonicTable = table( harmonicOrder,... harmonicMagnitude,... harmonicRms,... harmonicPercentage,... 'VariableNames', {'Order','Magnitude','RMS','Percentage'}); display( harmonicTable );
harmonicTable = 10x4 table Order Magnitude RMS Percentage _____ _________ ______ __________ 1 1945.8 1375.9 100 5 218.86 154.75 11.248 7 105.83 74.835 5.439 11 85.135 60.2 4.3753 13 57.599 40.729 2.9602 17 50.417 35.65 2.5911 19 37.612 26.596 1.933 23 33.859 23.942 1.7401 25 26.507 18.743 1.3622 29 23.979 16.955 1.2323
全高調波歪み
関数 ee_calculate_ThdPercent を使用して、高調波データから THD (全高調波歪み) の割合を計算します。
thdPercent = ee_calculateThdPercent( harmonicOrder, harmonicMagnitude ); disp( [ 'Total Harmonic Distortion percentage = ' num2str( thdPercent ), ' %' ] );
Total Harmonic Distortion percentage = 14.1721 %
高調波のプロット
MATLAB の関数 bar を使用して、高調波データをプロットできます。便宜上、関数 ee_plotHarmonics は Simscape のログ変数から直接、棒グラフをプロットします。
ee_plotHarmonics( Voltage_Source_Currents ); h1_ee_rectifier_thd = gcf;
まとめ
この例では、MATLAB コマンド ラインから呼び出し可能な 3 つの関数を使用する方法を説明しています。これらの関数は Simscape のログ データを解析して、高調波の振幅の取得、全高調波歪みの割合の計算、高調波の振幅のプロットを実行します。
付録 - DC 抵抗値の計算式
三相整流器の AC ピーク入力 と DC 平均出力
の関係は、損失を無視して次のように計算できます。
定格電圧 (線間 RMS) とピーク相電圧
の関係は次の式で表されます。
定格 AC 電圧と平均 DC 電圧の関係は次の式で表されます。
特定の電力 を取り出すために必要な抵抗
は、次の式で表されます。