過去問一覧第二種電気主任技術者2012年度 機械・制御内部正帰還補償器を含む制御系の応答・外乱抑圧・安定判別
一覧に戻る

内部正帰還補償器を含む制御系の応答・外乱抑圧・安定判別

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

図のフィードバック制御系について,次の問に答えよ。ただし,R(s) R\left( s\right) ,E(s) E\left( s\right) ,U(s) U\left( s\right) ,D(s) D\left( s\right) 及びY(s) Y\left( s\right) は,目標値r(t) r\left( t\right) ,偏差e(t) e\left( t\right) ,制御入力u(t) u\left( t\right) ,外乱d(t) d\left( t\right) 及び出力y(t) y\left( t\right) をそれぞれラプラス変換したものであり,C(s) C\left( s\right) は点線で囲んだ補償器内の補償要素の伝達関数を表す。

問題図

(1) D(s)=0 D\left( s\right) =0 の場合,制御対象だけを取り出したとき,u(t) u\left( t\right) として単位ステップ応答を加えた時の出力応答y(t) y\left( t\right) を求めよ。

(2) 点線で囲んだ補償器だけを取り出したとき,E(s) E\left( s\right) からU(s) U\left( s\right) までの伝達関数を求めよ。

(3) 図のフィードバック制御系において,R(s)=0 R\left( s\right) =0 のとき,D(s) D\left( s\right) からE(s) E\left( s\right) までの伝達関数を求めよ。

(4) R(s)=0 R\left( s\right) =0 の場合,C(s) C\left( s\right) として,C(s)=sTs+1 \displaystyle C\left( s\right) =\frac {s}{Ts+1} を選んだとき,外乱d(t) d\left( t\right) がランプ関数d(t)=t (t≧0) d\left( t\right) =t \ \left( t≧0\right) で与えられるときの定常速度偏差を求めよ。

(5) (4)のC(s) C\left( s\right) を選んだとき,外乱d(t) d\left( t\right) の影響が偏差e(t) e\left( t\right) に現れないようにするには,C(s) C\left( s\right) の時定数T T をどのように選べばよいかを説明せよ。

(6) 点線で囲んだ補償器をK1+K2s \displaystyle K_{1}+\frac {K_{2}}{s} に置き換えたときのフィードバック制御系が安定となる条件を求めよ。

解説

本問は、内部正帰還ループを持つ補償器を含むフィードバック制御系について、時間応答解析・伝達関数の導出・定常外乱応答・安定判別までを総合的に問う良問です。

【重要公式・定理の確認】

  • 部分分数分解とラプラス逆変換:重根を持つ伝達関数 1s2(s+1)=As+Bs2+Cs+1\dfrac{1}{s^2(s+1)} = \dfrac{A}{s} + \dfrac{B}{s^2} + \dfrac{C}{s+1} の展開
  • ブロック線図の基本関係:正帰還フィードバック要素の等価伝達関数 G1−GH\dfrac{G}{1 - GH}
  • 最終値の定理:系が安定なとき、定常偏差は lim⁡t→∞e(t)=lim⁡s→0sE(s)\displaystyle \lim_{t \to \infty} e(t) = \lim_{s \to 0} s E(s)
  • ラウスの安定判別法:3次特性多項式 s3+a2s2+a1s+a0=0s^3 + a_2 s^2 + a_1 s + a_0 = 0 が安定となる条件は、全係数が正かつ a2a1−a0>0a_2 a_1 - a_0 > 0

(1) 制御対象単体の単位ステップ応答 y(t)y(t)

外乱 D(s)=0D(s)=0 とし、制御対象 1s(s+1)\dfrac{1}{s(s+1)} の入力として単位ステップ関数 u(t)=1u(t)=1(ラプラス変換 U(s)=1sU(s) = \dfrac{1}{s})を与えたときの出力 Y(s)Y(s) は、次のようになります。

Y(s)=1s(s+1)U(s)=1s2(s+1)Y(s) = \frac{1}{s(s+1)} U(s) = \frac{1}{s^2(s+1)}

これを部分分数分解します。

1s2(s+1)=As+Bs2+Cs+1\frac{1}{s^2(s+1)} = \frac{A}{s} + \frac{B}{s^2} + \frac{C}{s+1}

両辺に s2(s+1)s^2(s+1) を掛けると、

1=As(s+1)+B(s+1)+Cs2=(A+C)s2+(A+B)s+B1 = A s(s+1) + B(s+1) + C s^2 = (A + C)s^2 + (A + B)s + B

各次数の係数を比較すると、

  • 定数項:B=1B = 1
  • s1s^1 の係数:A+B=0  ⟹  A=−1A + B = 0 \implies A = -1
  • s2s^2 の係数:A+C=0  ⟹  C=1A + C = 0 \implies C = 1

したがって、Y(s)Y(s) は次のように展開されます。

Y(s)=−1s+1s2+1s+1Y(s) = -\frac{1}{s} + \frac{1}{s^2} + \frac{1}{s+1}

各項をラプラス逆変換すると、出力応答 y(t)y(t) が求まります。

y(t)=t−1+e−t(t≧0)y(t) = t - 1 + \mathrm{e}^{-t} \quad (t \geqq 0)

(2) 補償器単体の伝達関数 U(s)E(s)\dfrac{U(s)}{E(s)}

問題図

問題図
の点線で囲まれた補償器に注目します。
内部の加算点では、外部からの偏差 E(s)E(s) に、出力 U(s)U(s) を積分要素 1s\dfrac{1}{s} に通した信号が正帰還(加算)されています。

したがって、関係式は以下の通りです。

U(s)=C(s){E(s)+1sU(s)}U(s) = C(s) \left\{ E(s) + \frac{1}{s} U(s) \right\}

U(s)U(s) について整理します。

U(s)−C(s)sU(s)=C(s)E(s)U(s) - \frac{C(s)}{s} U(s) = C(s) E(s)(1−C(s)s)U(s)=C(s)E(s)\left( 1 - \frac{C(s)}{s} \right) U(s) = C(s) E(s)s−C(s)sU(s)=C(s)E(s)\frac{s - C(s)}{s} U(s) = C(s) E(s)

両辺を E(s)E(s) で割り、伝達関数を求めます。

U(s)E(s)=sC(s)s−C(s)\frac{U(s)}{E(s)} = \frac{s C(s)}{s - C(s)}

(3) 外乱 D(s)D(s) から偏差 E(s)E(s) までの伝達関数 E(s)D(s)\dfrac{E(s)}{D(s)}

目標値 R(s)=0R(s) = 0 のとき、偏差は E(s)=−Y(s)E(s) = -Y(s)(すなわち Y(s)=−E(s)Y(s) = -E(s))となります。
補償器全体の伝達関数を Gc(s)=sC(s)s−C(s)G_c(s) = \dfrac{s C(s)}{s - C(s)}、制御対象を P(s)=1s(s+1)P(s) = \dfrac{1}{s(s+1)} と置きます。

出力 Y(s)Y(s) は、制御入力からの応答と外乱の和であるため、

Y(s)=P(s)U(s)+D(s)=P(s)Gc(s)E(s)+D(s)Y(s) = P(s) U(s) + D(s) = P(s) G_c(s) E(s) + D(s)

Y(s)=−E(s)Y(s) = -E(s) を代入すると、

−E(s)=P(s)Gc(s)E(s)+D(s)-E(s) = P(s) G_c(s) E(s) + D(s){1+P(s)Gc(s)}E(s)=−D(s)\{ 1 + P(s) G_c(s) \} E(s) = -D(s)E(s)D(s)=−11+P(s)Gc(s)\frac{E(s)}{D(s)} = -\frac{1}{1 + P(s) G_c(s)}

ここで、一巡伝達関数 P(s)Gc(s)P(s) G_c(s) を計算します。

P(s)Gc(s)=1s(s+1)⋅sC(s)s−C(s)=C(s)(s+1){s−C(s)}P(s) G_c(s) = \frac{1}{s(s+1)} \cdot \frac{s C(s)}{s - C(s)} = \frac{C(s)}{(s+1)\{s - C(s)\} }

分母を共通通分して整理します。

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

分母を展開・整理すると、

(s+1){s−C(s)}+C(s)=s(s+1)−(s+1)C(s)+C(s)=s(s+1)−sC(s)=s{s+1−C(s)}(s+1)\{s - C(s)\} + C(s) = s(s+1) - (s+1)C(s) + C(s) = s(s+1) - s C(s) = s\{s + 1 - C(s)\}

したがって、求める伝達関数は次のようになります。

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

(4) ランプ外乱入力に対する定常速度偏差

C(s)=sTs+1C(s) = \dfrac{s}{Ts+1} を代入します。
ここで補償器全体の伝達関数 Gc(s)G_c(s) を計算すると、

Gc(s)=sC(s)s−C(s)=s⋅sTs+1s−sTs+1=s2s(Ts+1)−s=s2Ts2=1TG_c(s) = \frac{s C(s)}{s - C(s)} = \frac{s \cdot \dfrac{s}{Ts+1}}{s - \dfrac{s}{Ts+1}} = \frac{s^2}{s(Ts+1) - s} = \frac{s^2}{Ts^2} = \frac{1}{T}

となり、補償器全体は純粋な比例ゲイン 1T\dfrac{1}{T} として機能することが分かります。
これより、外乱から偏差への伝達関数は次のようにシンプルに表せます。

E(s)D(s)=−11+P(s)Gc(s)=−11+1Ts(s+1)=−Ts(s+1)Ts(s+1)+1=−Ts(s+1)Ts2+Ts+1\frac{E(s)}{D(s)} = -\frac{1}{1 + P(s) G_c(s)} = -\frac{1}{1 + \dfrac{1}{T s(s+1)}} = -\frac{Ts(s+1)}{Ts(s+1) + 1} = -\frac{Ts(s+1)}{Ts^2 + Ts + 1}

外乱がランプ関数 d(t)=t (t≧0)d(t) = t \ (t \geqq 0) のとき、D(s)=1s2D(s) = \dfrac{1}{s^2} であるため、

E(s)=−Ts(s+1)Ts2+Ts+1⋅1s2=−T(s+1)s(Ts2+Ts+1)E(s) = -\frac{Ts(s+1)}{Ts^2 + Ts + 1} \cdot \frac{1}{s^2} = -\frac{T(s+1)}{s(Ts^2 + Ts + 1)}

特性多項式 Ts2+Ts+1=0Ts^2 + Ts + 1 = 0 は T>0T > 0 において安定(極がすべて左半平面)であるため、最終値の定理が適用できます。

ess=lim⁡t→∞e(t)=lim⁡s→0sE(s)=lim⁡s→0(−T(s+1)Ts2+Ts+1)=−Te_{\mathrm{ss}} = \lim_{t \to \infty} e(t) = \lim_{s \to 0} s E(s) = \lim_{s \to 0} \left( -\frac{T(s+1)}{Ts^2 + Ts + 1} \right) = -T−T(偏差の絶対値としては T)-T \quad (\text{偏差の絶対値としては } T)

(5) 外乱の影響を偏差に現れないようにするための時定数 TT の選定

(4)で求めた外乱から偏差への伝達関数は、

E(s)D(s)=−Ts(s+1)Ts2+Ts+1\frac{E(s)}{D(s)} = -\frac{Ts(s+1)}{Ts^2 + Ts + 1}

であり、定常速度偏差は −T-T です。
したがって、外乱 d(t)d(t) の影響が偏差 e(t)e(t) に現れないようにするには、伝達関数の分子を 00 に近づけるため、次のように選定します。

時定数 TT をできるだけ小さく選ぶ(理想的には T→0T \to 0 とする)。

※ T→0T \to 0 とすることで E(s)D(s)→0\dfrac{E(s)}{D(s)} \to 0 となり、外乱の影響を過渡的にも定常的にも完全に抑圧できます(物理的には補償器の比例ゲイン 1T\dfrac{1}{T} を極めて大きくすることに相当します)。

(6) 補償器を K1+K2sK_1 + \dfrac{K_2}{s} に置き換えたときの安定条件

点線で囲まれた補償器全体を Gc(s)=K1+K2s=K1s+K2sG_c(s) = K_1 + \dfrac{K_2}{s} = \dfrac{K_1 s + K_2}{s} に置き換えます。
開ループ(一巡)伝達関数 L(s)L(s) は、

L(s)=Gc(s)P(s)=K1s+K2s⋅1s(s+1)=K1s+K2s2(s+1)L(s) = G_c(s) P(s) = \frac{K_1 s + K_2}{s} \cdot \frac{1}{s(s+1)} = \frac{K_1 s + K_2}{s^2(s+1)}

フィードバック制御系の閉ループ特性方程式 1+L(s)=01 + L(s) = 0 は、

s2(s+1)+(K1s+K2)=0  ⟹  s3+s2+K1s+K2=0s^2(s+1) + (K_1 s + K_2) = 0 \implies s^3 + s^2 + K_1 s + K_2 = 0

ラウスの安定判別法を適用するためにラウス数表を作成します。

s31K1s21K2s11⋅K1−1⋅K21=K1−K20s0K2\begin{array}{c|cc} s^3 & 1 & K_1 \\[8pt] s^2 & 1 & K_2 \\[8pt] s^1 & \dfrac{1 \cdot K_1 - 1 \cdot K_2}{1} = K_1 - K_2 & 0 \\[8pt] s^0 & K_2 & \end{array}

制御系が安定であるための必要十分条件は、ラウス数表の第1列のすべての要素が同符号(正)であることです。

  1. s3s^3 の係数:1>01 > 0
  2. s2s^2 の係数:1>01 > 0
  3. s1s^1 の要素:K1−K2>0  ⟹  K1>K2K_1 - K_2 > 0 \implies K_1 > K_2
  4. s0s^0 の要素:K2>0K_2 > 0

これらを同時に満たす範囲をまとめると、求める条件が得られます。

K1>K2>0K_1 > K_2 > 0

【最終解答まとめ】

  • (1) y(t)=t−1+e−t(t≧0)y(t) = t - 1 + \mathrm{e}^{-t} \quad (t \geqq 0)
  • (2) U(s)E(s)=sC(s)s−C(s)\dfrac{U(s)}{E(s)} = \dfrac{s C(s)}{s - C(s)}
  • (3) E(s)D(s)=−(s+1){s−C(s)}s{s+1−C(s)}\dfrac{E(s)}{D(s)} = -\dfrac{(s+1)\{s - C(s)\}}{s\{s + 1 - C(s)\} }
  • (4) −T-T
  • (5) 時定数 TT をできるだけ小さく選ぶ(理想的には T→0T \to 0 とする)
  • (6) K1>K2>0K_1 > K_2 > 0