過去問一覧第二種電気主任技術者2018年度 機械・制御2自由度制御系の伝達関数導出とPIパラメータ設計
一覧に戻る

2自由度制御系の伝達関数導出とPIパラメータ設計

PID
第二種電気主任技術者 2018年度 機械・制御 問4

図1のような2 2 自由度制御系がある。ここでG(s) G(s) は制御対象,C(s) C(s) 及びF(s) F(s) は補償器である。また,R(s) R(s) は目標値,D(s) D(s) は外乱,Y(s) Y(s) は制御量,E(s) E(s) は制御偏差である。この制御系について,次の問に答えよ。

(1) 図1に示すフィードバック補償器C(s) C(s) の係数KP K_{\mathrm {P}} 及びTI T_{\mathrm {I}} の名称を答えよ。

(2) KP=10 K_{\mathrm {P}}=10 ,TI=0.1 T_{\mathrm {I}}=0.1 のとき,C(s) C(s) の角周波数ω [rad/s] \omega \ \mathrm {[rad/s]} に対するゲイン特性の概形を答案用紙に印刷されている図2に折れ線近似で図示せよ。

(3) R(s)=0 R(s)=0 として,外乱D(s) D(s) から制御偏差E(s) E(s) までの閉ループ伝達関数を求めよ。

(4) 上記小問(3)で求めた閉ループ伝達関数において,固有角周波数が5 rad/s 5 \ \mathrm {rad/s} ,減衰係数が0.7 0.7 となるように,KP K_{\mathrm {P}} とTI T_{\mathrm {I}} の値を定めよ。

(5) D(s)=0 D(s)=0 として,目標値R(s) R(s) から制御量Y(s) Y(s) までの閉ループ伝達関数をF(s) F(s) ,C(s) C(s) ,G(s) G(s) を用いて求めよ。

(6) 上記小問(3)で求めた閉ループ伝達関数はF(s) F(s) に依存しない。また,上記小問(5)で求めた閉ループ伝達関数はC(s) C(s) に依存しない。これは制御系にどんな特長をもたらすか答えよ。この性質を利用しているのが2 2 自由度制御系である。

問題図

問題図

解説

本問は、目標値追従特性と外乱抑制特性をそれぞれ独立に調整できる2自由度制御系に関する総合問題です。

【使用する基本公式と概念】

  • PI補償器(比例・積分補償器)の伝達関数:C(s)=KP(1+1TIs)C(s) = K_{\mathrm{P}}\left(1 + \frac{1}{T_{\mathrm{I}}s}\right)(KPK_{\mathrm{P}}:比例ゲイン、TIT_{\mathrm{I}}:積分時間)
  • ゲイン[dB\mathrm{dB}]の定義:g=20log⁡10∣C(jω)∣g = 20\log_{10}|C(\mathrm{j}\omega)|
  • 標準2次系の特性方程式:s2+2ζωns+ωn2=0s^2 + 2\zeta\omega_{\mathrm{n}}s + \omega_{\mathrm{n}}^2 = 0(ωn\omega_{\mathrm{n}}:固有角周波数、ζ\zeta:減衰係数)
  • 2自由度制御系の本質:目標値から出力までの伝達関数と、外乱から出力(または偏差)までの伝達関数を、フィードフォワード補償器とフィードバック補償器によって独立に設計可能であること。

(1) フィードバック補償器 C(s)C(s) の係数 KPK_{\mathrm{P}} 及び TIT_{\mathrm{I}} の名称

伝達関数 C(s)=KP(1+1TIs)=KP+KPTIsC(s) = K_{\mathrm{P}}\left(1 + \dfrac{1}{T_{\mathrm{I}}s}\right) = K_{\mathrm{P}} + \dfrac{K_{\mathrm{P}}}{T_{\mathrm{I}}s} は、比例動作と積分動作を組み合わせたPI補償器の標準的な表現形式です。

  • 第1項にかかる係数 KPK_{\mathrm{P}} は、偏差の大きさに比例した制御出力を決める比例ゲイン(または比例感度)です。
  • 積分項の分母にある係数 TIT_{\mathrm{I}} は、積分動作の強さを時間の次元で表す積分時間です。

KPK_{\mathrm{P}}:比例ゲイン、TIT_{\mathrm{I}}:積分時間

(2) C(s)C(s) のボード線図におけるゲイン特性の折れ線近似

KP=10K_{\mathrm{P}}=10、TI=0.1 sT_{\mathrm{I}}=0.1\ \mathrm{s} を代入すると、

C(s)=10(1+10.1s)=10⋅0.1s+10.1sC(s) = 10\left(1 + \frac{1}{0.1s}\right) = 10\cdot\frac{0.1s+1}{0.1s}

周波数伝達関数 C(jω)C(\mathrm{j}\omega) の大きさ(絶対値)は、

∣C(jω)∣=10∣1+10jω∣=101+(10ω)2|C(\mathrm{j}\omega)| = 10\left|1 + \frac{10}{\mathrm{j}\omega}\right| = 10\sqrt{1 + \left(\frac{10}{\omega}\right)^2}

この伝達関数の折点角周波数は、虚数部の大きさと実数部の大きさが等しくなる点、すなわち 1TI=10.1=10 rad/s\dfrac{1}{T_{\mathrm{I}}} = \dfrac{1}{0.1} = 10\ \mathrm{rad/s} です。

① 低周波域(ω≤10 rad/s\omega \le 10\ \mathrm{rad/s})の漸近線:
ω≪10\omega \ll 10 では積分の効果が支配的となるため、∣C(jω)∣≈10⋅10ω=100ω|C(\mathrm{j}\omega)| \approx 10\cdot\frac{10}{\omega} = \frac{100}{\omega}ゲインは、20log⁡10∣C(jω)∣≈20log⁡10100ω=40−20log⁡10ω [dB]20\log_{10}|C(\mathrm{j}\omega)| \approx 20\log_{10}\frac{100}{\omega} = 40 - 20\log_{10}\omega\ \mathrm{[dB]}したがって、傾き −20 dB/dec-20\ \mathrm{dB/dec} の直線となります。
・ω=1 rad/s\omega = 1\ \mathrm{rad/s} のとき:40 dB40\ \mathrm{dB}
・ω=10 rad/s\omega = 10\ \mathrm{rad/s} のとき:20 dB20\ \mathrm{dB}

② 高周波域(ω≥10 rad/s\omega \ge 10\ \mathrm{rad/s})の漸近線:
ω≫10\omega \gg 10 では比例項が支配的となるため、∣C(jω)∣≈10|C(\mathrm{j}\omega)| \approx 10ゲインは、20log⁡1010=20 [dB]20\log_{10}10 = 20\ \mathrm{[dB]}したがって、傾き 0 dB/dec0\ \mathrm{dB/dec}(一定値 20 dB20\ \mathrm{dB} の水平線)となります。

図2において、点 (1 rad/s, 40 dB)(1\ \mathrm{rad/s},\ 40\ \mathrm{dB}) から点 (10 rad/s, 20 dB)(10\ \mathrm{rad/s},\ 20\ \mathrm{dB}) まで傾き −20 dB/dec-20\ \mathrm{dB/dec} の直線を描き、ω=10 rad/s\omega = 10\ \mathrm{rad/s} から ω=100 rad/s\omega = 100\ \mathrm{rad/s} 以降はゲイン 20 dB20\ \mathrm{dB} の水平な直線を描く。

(3) 外乱 D(s)D(s) から制御偏差 E(s)E(s) までの閉ループ伝達関数

R(s)=0R(s)=0 と置くと、図1より目標値からの入力経路はすべて 00 となります。各信号の関係式を立てます。

V(s)=0E(s)=V(s)−Y(s)=−Y(s)U(s)=C(s)E(s)=−C(s)Y(s)Y(s)=G(s){U(s)+D(s)}=G(s){−C(s)Y(s)+D(s)}\begin{aligned} V(s) &= 0 \\[5pt] E(s) &= V(s) - Y(s) = -Y(s) \\[5pt] U(s) &= C(s)E(s) = -C(s)Y(s) \\[5pt] Y(s) &= G(s)\{U(s) + D(s)\} = G(s)\{-C(s)Y(s) + D(s)\} \end{aligned}

これより Y(s)Y(s) をまとめると、

{1+C(s)G(s)}Y(s)=G(s)D(s)  ⟹  Y(s)=G(s)1+C(s)G(s)D(s)\{1 + C(s)G(s)\}Y(s) = G(s)D(s) \implies Y(s) = \frac{G(s)}{1 + C(s)G(s)}D(s)

E(s)=−Y(s)E(s) = -Y(s) であるから、

E(s)D(s)=−G(s)1+C(s)G(s)\frac{E(s)}{D(s)} = -\frac{G(s)}{1 + C(s)G(s)}

ここに G(s)=14s+1G(s) = \dfrac{1}{4s+1}、C(s)=KP(TIs+1TIs)C(s) = K_{\mathrm{P}}\left(\dfrac{T_{\mathrm{I}}s+1}{T_{\mathrm{I}}s}\right) を代入します。

E(s)D(s)=−14s+11+KPTIs+1TIs⋅14s+1=−TIs(4s+1)TIs+KP(TIs+1)=−TIs4TIs2+TI(1+KP)s+KP\begin{aligned} \frac{E(s)}{D(s)} &= -\frac{\dfrac{1}{4s+1}}{1 + K_{\mathrm{P}}\dfrac{T_{\mathrm{I}}s+1}{T_{\mathrm{I}}s}\cdot\dfrac{1}{4s+1}} \\[10pt] &= -\frac{T_{\mathrm{I}}s}{(4s+1)T_{\mathrm{I}}s + K_{\mathrm{P}}(T_{\mathrm{I}}s+1)} \\[10pt] &= -\frac{T_{\mathrm{I}}s}{4T_{\mathrm{I}}s^2 + T_{\mathrm{I}}(1+K_{\mathrm{P}})s + K_{\mathrm{P}}} \end{aligned}

E(s)D(s)=−TIs4TIs2+TI(1+KP)s+KP\frac{E(s)}{D(s)} = -\frac{T_{\mathrm{I}}s}{4T_{\mathrm{I}}s^2 + T_{\mathrm{I}}(1+K_{\mathrm{P}})s + K_{\mathrm{P}}}

(4) KPK_{\mathrm{P}} と TIT_{\mathrm{I}} の値の決定

(3)で求めた伝達関数の分母多項式において、最高次 s2s^2 の係数を 11 にするために分母分子を 4TI4T_{\mathrm{I}} で割ります。

E(s)D(s)=−s4s2+1+KP4s+KP4TI\frac{E(s)}{D(s)} = -\frac{\dfrac{s}{4}}{s^2 + \dfrac{1+K_{\mathrm{P}}}{4}s + \dfrac{K_{\mathrm{P}}}{4T_{\mathrm{I}}}}

標準的な2次系の特性方程式 s2+2ζωns+ωn2=0s^2 + 2\zeta\omega_{\mathrm{n}}s + \omega_{\mathrm{n}}^2 = 0 と係数を比較します。

与えられた条件は、固有角周波数 ωn=5 rad/s\omega_{\mathrm{n}} = 5\ \mathrm{rad/s}、減衰係数 ζ=0.7\zeta = 0.7 です。

① ss の1次の係数の比較:

1+KP4=2ζωn=2×0.7×5=7  ⟹  1+KP=28  ⟹  KP=27\frac{1+K_{\mathrm{P}}}{4} = 2\zeta\omega_{\mathrm{n}} = 2\times 0.7\times 5 = 7 \implies 1+K_{\mathrm{P}} = 28 \implies K_{\mathrm{P}} = 27

② 定数項の比較:

KP4TI=ωn2=52=25\frac{K_{\mathrm{P}}}{4T_{\mathrm{I}}} = \omega_{\mathrm{n}}^2 = 5^2 = 25

KP=27K_{\mathrm{P}}=27 を代入すると、

274TI=25  ⟹  100TI=27  ⟹  TI=0.27 s\frac{27}{4T_{\mathrm{I}}} = 25 \implies 100T_{\mathrm{I}} = 27 \implies T_{\mathrm{I}} = 0.27\ \mathrm{s}

KP=27,TI=0.27 sK_{\mathrm{P}} = 27,\quad T_{\mathrm{I}} = 0.27\ \mathrm{s}

(5) 目標値 R(s)R(s) から制御量 Y(s)Y(s) までの閉ループ伝達関数

D(s)=0D(s)=0 として、図1の信号の流れから方程式を立てます。

V(s)=F(s)R(s)E(s)=V(s)−Y(s)=F(s)R(s)−Y(s)U(s)=F(s)G(s)R(s)+C(s)E(s)=F(s)G(s)R(s)+C(s){F(s)R(s)−Y(s)}Y(s)=G(s)U(s)\begin{aligned} V(s) &= F(s)R(s) \\[5pt] E(s) &= V(s) - Y(s) = F(s)R(s) - Y(s) \\[5pt] U(s) &= \frac{F(s)}{G(s)}R(s) + C(s)E(s) = \frac{F(s)}{G(s)}R(s) + C(s)\{F(s)R(s) - Y(s)\} \\[5pt] Y(s) &= G(s)U(s) \end{aligned}

Y(s)Y(s) の式に U(s)U(s) を代入します。

Y(s)=G(s)[F(s)G(s)R(s)+C(s){F(s)R(s)−Y(s)}]Y(s)=F(s)R(s)+C(s)G(s)F(s)R(s)−C(s)G(s)Y(s){1+C(s)G(s)}Y(s)=F(s){1+C(s)G(s)}R(s)\begin{aligned} Y(s) &= G(s)\left[\frac{F(s)}{G(s)}R(s) + C(s)\{F(s)R(s) - Y(s)\}\right] \\[5pt] Y(s) &= F(s)R(s) + C(s)G(s)F(s)R(s) - C(s)G(s)Y(s) \\[5pt] \{1 + C(s)G(s)\}Y(s) &= F(s)\{1 + C(s)G(s)\}R(s) \end{aligned}

両辺を 1+C(s)G(s)1+C(s)G(s) で割ると、

Y(s)R(s)=F(s)\frac{Y(s)}{R(s)} = F(s)

(なお、図1より F(s)=12s+1F(s) = \dfrac{1}{2s+1} です)

Y(s)R(s)=F(s)(=12s+1)\frac{Y(s)}{R(s)} = F(s)\quad \left(=\frac{1}{2s+1}\right)

(6) 2自由度制御系がもたらす特長

小問(3)より、外乱に対する閉ループ特性は F(s)F(s) に依存せずフィードバック補償器 C(s)C(s) によって決まります。一方、小問(5)より、目標値に対する追従特性は C(s)C(s) に依存せずフィードフォワード補償器 F(s)F(s) のみによって決定されます。

通常の1自由度制御系では、同一の補償器で目標値追従性と外乱抑制性(安定余裕)の双方を満足させなければならず、トレードオフが生じます。これに対し、本方式では両者を完全に分離して取り扱うことができます。

目標値に対する応答特性(目標値追従特性)と外乱に対する抑制特性(外乱抑圧・安定化特性)を互いに干渉することなく、独立に設計・調整できる特長をもつ。

【総括・解答一覧】

  • (1) KPK_{\mathrm{P}}:比例ゲイン、TIT_{\mathrm{I}}:積分時間
  • (2) ω=10 rad/s\omega=10\ \mathrm{rad/s} まで傾き −20 dB/dec-20\ \mathrm{dB/dec}(ω=1\omega=1 で 40 dB40\ \mathrm{dB}、ω=10\omega=10 で 20 dB20\ \mathrm{dB})、ω≥10 rad/s\omega \ge 10\ \mathrm{rad/s} では 20 dB20\ \mathrm{dB} の水平線
  • (3) E(s)D(s)=−TIs4TIs2+TI(1+KP)s+KP\displaystyle \frac{E(s)}{D(s)} = -\frac{T_{\mathrm{I}}s}{4T_{\mathrm{I}}s^2 + T_{\mathrm{I}}(1+K_{\mathrm{P}})s + K_{\mathrm{P}}}
  • (4) KP=27K_{\mathrm{P}} = 27、TI=0.27 sT_{\mathrm{I}} = 0.27\ \mathrm{s}
  • (5) Y(s)R(s)=F(s)=12s+1\displaystyle \frac{Y(s)}{R(s)} = F(s) = \frac{1}{2s+1}
  • (6) 目標値追従特性と外乱抑制特性をそれぞれ独立に設計できる