過去問一覧第二種電気主任技術者2011年度 機械・制御2自由度制御系の伝達関数導出とPID補償器の安定条件判別
一覧に戻る

2自由度制御系の伝達関数導出とPID補償器の安定条件判別

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

図のフィードバック制御系について,次の問に答えよ。ただし,R(s) R\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) ,操作量u(t) u\left( t\right) ,外乱d(t) d\left( t\right) ,出力y(t) y\left( t\right) のラプラス変換を表す。また,G(s) G\left( s\right) は制御対象の伝達関数,F(s) F\left( s\right) 及びK(s) K\left( s\right) は補償器の伝達関数を表す。

(1) R(s)=0 R\left( s\right) =0 のとき,D(s) D\left( s\right) からY(s) Y\left( s\right) までの伝達関数を求めよ。

(2) D(s)=0 D\left( s\right) =0 のとき,R(s) R\left( s\right) からY(s) Y\left( s\right) までの伝達関数を求めよ。

(3) 図において,G(s)=1s2 \displaystyle G\left( s\right) =\frac {1}{s^{2}} ,F(s)=cs2+as+b \displaystyle F\left( s\right) =\frac {c}{s^{2}+as+b} ,K(s)=KP(1+1TIs+TDs) \displaystyle K\left( s\right) =K_{\mathrm {P}}\left( 1+\frac {1}{T_{\mathrm {I}}s}+T_{\mathrm {D}}s \right) とおく。D(s)=0 D\left( s\right) =0 のとき,R(s) R\left( s\right) からY(s) Y\left( s\right) までの応答特性として,単位ステップ関数の目標値r(t)=1 r\left( t\right) =1 に対して出力y(t) y\left( t\right) の定常値が1 1 となり,かつ,減衰定数が0.8 0.8 ,固有角周波数が10 [rad/s] 10 \ \mathrm {[rad / s]} を満たす2 2 次系の補償器F(s) F\left( s\right) の係数a a ,b b ,c c を求めよ。

(4) 上記(3)の補償器K(s) K\left( s\right) の名称を答えよ。また,各係数KP K_{\mathrm {P}} ,TI T_{\mathrm {I}} ,TD T_{\mathrm {D}} の名称についても答えよ。

(5) 上記(3)において,F(s) F\left( s\right) は安定な補償器であり,図の制御系全体の安定性はF(s) F\left( s\right) にはよらない。制御系全体が安定となるために補償器K(s) K\left( s\right) の係数KP K_{\mathrm {P}} ,TI T_{\mathrm {I}} ,TD T_{\mathrm {D}} が満たさなければならない条件を求めよ。ただし,KP>0 K_{\mathrm {P}}\gt 0 ,TI>0 T_{\mathrm {I}}\gt 0 ,TD>0 T_{\mathrm {D}}\gt 0 とする。

問題図

解説

本問は、目標値に対する追従特性と外乱に対する抑制特性を独立に設計できる代表的な制御方式である「2自由度制御系」を題材とした問題です。

ブロック線図の信号関係の丁寧な整理、2次遅れ系の標準形、最終値の定理、および3次特性方程式に対するラウスの安定判別法の基礎がしっかりと身についていれば、確実に完答できる良問です。

【ブロック線図における基本関係式の整理】

添付図

問題図
より、目標値 R(s)R(s)、外乱 D(s)D(s)、操作量 U(s)U(s)、出力 Y(s)Y(s) の関係を数式化します。

下側の偏差信号(加算点の出力)は F(s)R(s)−Y(s)F(s)R(s) - Y(s) です。したがって、操作量 U(s)U(s) を作る加算点の入力は F(s)G(s)R(s)\displaystyle \frac{F(s)}{G(s)}R(s) と K(s){F(s)R(s)−Y(s)}K(s)\{F(s)R(s) - Y(s)\} となり、関係式は以下のようになります。

U(s)=F(s)G(s)R(s)+K(s){F(s)R(s)−Y(s)}⋯①Y(s)=G(s)U(s)+D(s)⋯②\begin{aligned} U(s) &= \frac{F(s)}{G(s)}R(s) + K(s)\{F(s)R(s) - Y(s)\} & \quad \cdots ① \\[5pt] Y(s) &= G(s)U(s) + D(s) & \quad \cdots ② \end{aligned}

(1) R(s)=0R(s)=0 のときの外乱伝達関数 Y(s)D(s)\frac{Y(s)}{D(s)}

R(s)=0R(s)=0 を式①に代入すると、制御系の入力側からの信号はすべて零となるため、

U(s)=−K(s)Y(s)U(s) = -K(s)Y(s)

これを式②に代入して整理します。

Y(s)=G(s){−K(s)Y(s)}+D(s)Y(s)=−G(s)K(s)Y(s)+D(s){1+G(s)K(s)}Y(s)=D(s)\begin{aligned} Y(s) &= G(s)\{-K(s)Y(s)\} + D(s) \\[5pt] Y(s) &= -G(s)K(s)Y(s) + D(s) \\[5pt] \{1 + G(s)K(s)\}Y(s) &= D(s) \end{aligned}

したがって、D(s)D(s) から Y(s)Y(s) までの伝達関数は次のようになります。

Y(s)D(s)=11+G(s)K(s)\frac{Y(s)}{D(s)} = \frac{1}{1 + G(s)K(s)}

(2) D(s)=0D(s)=0 のときの目標値伝達関数 Y(s)R(s)\frac{Y(s)}{R(s)}

D(s)=0D(s)=0 を式②に代入すると Y(s)=G(s)U(s)Y(s) = G(s)U(s) となります。式①の両辺に G(s)G(s) を乗じると、

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

両辺の {1+G(s)K(s)}\{1 + G(s)K(s)\} で除すると、閉ループのフィードバック補償器 K(s)K(s) が相殺され、目標値伝達関数は F(s)F(s) そのものになります。

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

※これにより、目標値応答の特性を F(s)F(s) のみで独立に設計できるという2自由度制御系の利点が確認できます。

(3) 補償器 F(s)F(s) の係数 a,b,ca, b, c の導出

(2)より、D(s)=0D(s)=0 のときの伝達関数は F(s)=cs2+as+b\displaystyle F(s) = \frac{c}{s^2 + as + b} です。

減衰定数(減衰比)ζ\zeta、固有角周波数 ωn\omega_{\mathrm{n}} をもつ標準的な2次系の分母多項式は s2+2ζωns+ωn2s^2 + 2\zeta \omega_{\mathrm{n}}s + \omega_{\mathrm{n}}^2 で表されます。

問題の条件 ζ=0.8\zeta = 0.8、ωn=10 [rad/s]\omega_{\mathrm{n}} = 10\ \mathrm{[rad/s]} を代入すると、

a=2ζωn=2×0.8×10=16b=ωn2=102=100\begin{aligned} a &= 2\zeta \omega_{\mathrm{n}} = 2 \times 0.8 \times 10 = 16 \\[5pt] b &= \omega_{\mathrm{n}}^2 = 10^2 = 100 \end{aligned}

次に、単位ステップ入力 r(t)=1  ⟺  R(s)=1sr(t)=1 \iff R(s) = \frac{1}{s} に対する出力 y(t)y(t) の定常値を最終値の定理を用いて求めます。

lim⁡t→∞y(t)=lim⁡s→0sY(s)=lim⁡s→0s(F(s)1s)=F(0)=cb\lim_{t\to\infty} y(t) = \lim_{s\to 0} s Y(s) = \lim_{s\to 0} s \left( F(s) \frac{1}{s} \right) = F(0) = \frac{c}{b}

定常値が 11 であることから、cb=1  ⟹  c=b=100\displaystyle \frac{c}{b} = 1 \implies c = b = 100 と求まります。

a=16,b=100,c=100a = 16, \quad b = 100, \quad c = 100

(4) 補償器 K(s)K(s) および各係数の名称

補償器 K(s)=KP(1+1TIs+TDs)K(s) = K_{\mathrm{P}}\left( 1 + \frac{1}{T_{\mathrm{I}}s} + T_{\mathrm{D}}s \right) は、比例(Proportional)、積分(Integral)、微分(Derivative)の3動作を行う制御器です。

K(s)K(s) の名称:PID補償器(または PID調節器、PID制御器)
KPK_{\mathrm{P}} の名称:比例ゲイン
TIT_{\mathrm{I}} の名称:積分時間
TDT_{\mathrm{D}} の名称:微分時間

(5) 制御系全体が安定となるための条件

(1)で求めた伝達関数より、フィードバック制御ループの特性方程式は 1+G(s)K(s)=01 + G(s)K(s) = 0 です。

与えられた G(s)=1s2G(s) = \frac{1}{s^2} および K(s)=KP(1+1TIs+TDs)K(s) = K_{\mathrm{P}}\left( 1 + \frac{1}{T_{\mathrm{I}}s} + T_{\mathrm{D}}s \right) を代入します。

1+1s2⋅KP(1+1TIs+TDs)=01 + \frac{1}{s^2} \cdot K_{\mathrm{P}}\left( 1 + \frac{1}{T_{\mathrm{I}}s} + T_{\mathrm{D}}s \right) = 0

両辺に s2s^2 を乗じて整理すると、

s2+KP+KPTIs+KPTDs=0s^2 + K_{\mathrm{P}} + \frac{K_{\mathrm{P}}}{T_{\mathrm{I}}s} + K_{\mathrm{P}}T_{\mathrm{D}}s = 0

さらに両辺に ss を乗じて次数の高い順に並べ替えます。

s3+KPTDs2+KPs+KPTI=0s^3 + K_{\mathrm{P}}T_{\mathrm{D}}s^2 + K_{\mathrm{P}}s + \frac{K_{\mathrm{P}}}{T_{\mathrm{I}}} = 0

この3次多項式に対してラウスの数表を作成します。

s31KPs2KPTDKPTIs1KP2TD−KPTIKPTD0s0KPTI\begin{array}{c|c} s^3 & 1 & K_{\mathrm{P}} \\[8pt] s^2 & K_{\mathrm{P}}T_{\mathrm{D}} & \displaystyle\frac{K_{\mathrm{P}}}{T_{\mathrm{I}}} \\[8pt] s^1 & \displaystyle\frac{K_{\mathrm{P}}^2 T_{\mathrm{D}} - \displaystyle\frac{K_{\mathrm{P}}}{T_{\mathrm{I}}}}{K_{\mathrm{P}}T_{\mathrm{D}}} & 0 \\[8pt] s^0 & \displaystyle\frac{K_{\mathrm{P}}}{T_{\mathrm{I}}} & \end{array}

系が安定であるための必要十分条件は、ラウス表の第1列の要素がすべて正(同符号)であることです。KP>0,TI>0,TD>0K_{\mathrm{P}}>0, T_{\mathrm{I}}>0, T_{\mathrm{D}}>0 の条件より、第1行、第2行、第4行の要素は正であるため、第3行の要素が正となる条件を求めます。

KP2TD−KPTIKPTD>0\frac{K_{\mathrm{P}}^2 T_{\mathrm{D}} - \displaystyle\frac{K_{\mathrm{P}}}{T_{\mathrm{I}}}}{K_{\mathrm{P}}T_{\mathrm{D}}} > 0

分母 KPTD>0K_{\mathrm{P}}T_{\mathrm{D}} > 0 であるため、分子が正であれば成立します。

KP2TD−KPTI>0KP(KPTD−1TI)>0\begin{aligned} K_{\mathrm{P}}^2 T_{\mathrm{D}} - \frac{K_{\mathrm{P}}}{T_{\mathrm{I}}} &> 0 \\[5pt] K_{\mathrm{P}} \left( K_{\mathrm{P}}T_{\mathrm{D}} - \frac{1}{T_{\mathrm{I}}} \right) &> 0 \end{aligned}

KP>0K_{\mathrm{P}} > 0 および TI>0T_{\mathrm{I}} > 0 より、

KPTITD>1K_{\mathrm{P}}T_{\mathrm{I}}T_{\mathrm{D}} > 1

(※ KPTD>1TIK_{\mathrm{P}}T_{\mathrm{D}} > \displaystyle\frac{1}{T_{\mathrm{I}}} や TITD>1KPT_{\mathrm{I}}T_{\mathrm{D}} > \displaystyle\frac{1}{K_{\mathrm{P}}} などの同値な表現でも正解です。)

【解答のまとめ】

  • (1) Y(s)D(s)=11+G(s)K(s)\displaystyle \frac{Y(s)}{D(s)} = \frac{1}{1 + G(s)K(s)}
  • (2) Y(s)R(s)=F(s)\displaystyle \frac{Y(s)}{R(s)} = F(s)
  • (3) a=16,b=100,c=100a = 16, \quad b = 100, \quad c = 100
  • (4) K(s)K(s):PID補償器、KPK_{\mathrm{P}}:比例ゲイン、TIT_{\mathrm{I}}:積分時間、TDT_{\mathrm{D}}:微分時間
  • (5) KPTITD>1K_{\mathrm{P}}T_{\mathrm{I}}T_{\mathrm{D}} > 1