$K_p=1.2\,\tau/(K\cdot L),\quad T_i=2L,\quad T_d=0.5L$
Kp・積分時間Ti・微分時間Tdをスライダーで操作し、閉ループステップ応答をリアルタイム表示。オーバーシュート・整定時間・IAEを即時計算して最適PIDゲインを直感的に学ぼう。
PID制御は、目標値と現在値の差(偏差 $e$)に対し、比例(P)・積分(I)・微分(D)の3つの操作を組み合わせて操作量 $u$ を決めます。
$u(t) = K_p\,e(t) + K_i\!\int_0^t\! e\,d\tau + K_d\,\dfrac{de}{dt}$
| 項 | 作用 | 効果と副作用 |
|---|---|---|
| 比例 P ($K_p$) | 偏差に比例 | 応答を速くするが定常偏差が残り、大きすぎると振動 |
| 積分 I ($K_i$) | 偏差の累積 | 定常偏差を除去するが、行き過ぎ(オーバーシュート)を招く |
| 微分 D ($K_d$) | 偏差の変化率 | 応答を安定化・先読みするが、ノイズに敏感 |
各ゲインの調整(チューニング)が制御性能を左右します。代表的なジーグラ・ニコルス法では、I・Dを切り、$K_p$ を上げて系が持続振動する限界ゲイン $K_u$ と振動周期 $T_u$ を求め、そこから経験式で $K_p, K_i, K_d$ を設定します。
一般に、立ち上がりを速く=$K_p$↑、定常偏差を消す=$K_i$↑、オーバーシュート/振動を抑える=$K_d$↑ が基本方針です。ただし過大なゲインは不安定化を招くため、応答(立ち上がり時間・行き過ぎ量・整定時間)を見ながら調整します。本シミュレーターで各ゲインを変え、ステップ応答の変化を確認できます。
化学・プラントプロセス制御:反応器の温度、蒸留塔の圧力、タンクの液面レベルなど、安定性が求められる連続プロセスで広く使用されます。PIDパラメータのチューニング不良は、製品品質のばらつきやエネルギー損失に直結します。
モーションコントロール:CNC工作機械や産業用ロボットの位置・速度制御に用いられます。ここでは微分動作(D)が重要で、位置決めのオーバーシュートを抑えつつ、高速で正確に目標位置に到達させるために調整されます。
自動車制御:クルーズコントロール(定速走行装置)は代表例です。設定速度と現在速度の偏差をPID制御で計算し、スロットル開度を調整します。坂道でも速度を一定に保つためには積分動作(I)が不可欠です。
家電・ビル設備:エアコンの室温制御、給湯器の湯温制御、ビルの空調システムなど身近なところにも応用されています。省エネと快適性の両立のために、外気温変化に対する応答性がチューニングのポイントになります。
まず、「積分時間(Ti)は短いほど良い」という誤解があります。確かにTiを短く(積分ゲインを強く)すると定常偏差は速く解消されますが、過度に短くするとシステムが不安定になり、持続的な振動(積分ワインドアップのリスクも含む)を引き起こします。例えば、Tiを0.1秒など極端に小さく設定すると、グラフは目標値を行ったり来たりして落ち着かなくなるはずです。実務では、まず比例ゲイン(Kp)で基本的な応答を整え、その後、振動を起こさない範囲でTiを調整するのがコツです。
次に、微分動作(D)は「万能の先読み機能」ではないという点。微分は確かに変化率を予測しますが、ノイズに対して極めて敏感です。実際のセンサ信号には必ずノイズが含まれており、これを微分するとノイズが巨大な制御信号に増幅されてしまいます。そのため、実機では微分ゲインを控えめに設定したり、微分動作にローパスフィルタを併用することがほとんどです。シミュレーターでは理想的な信号を使っているので、この危険性を実感しにくい点に注意が必要です。
最後に、Ziegler-Nichols法は「最適解」ではなく「出発点」です。この方法で得られたパラメータは、やや攻撃的(オーバーシュートが大きめ)で、ロバスト性(外乱への強さ)が必ずしも高くありません。例えば、ZN法で設定したパラメータでシミュレーションを実行し、その後、オーバーシュートを抑えたい場合はKpを少し下げ、Tiを少し上げるといった微調整が不可欠です。「一度設定したら終わり」ではなく、実際の応答を見ながら、トレードオフ(応答速度 vs 安定性)を考慮して手を加えるのが実践の現場です。
1次遅れ系プロセス(時定数τ=2秒、ゲイン=1)を想定します。目標値SP=1.0、Kp=2.0、Ti=3.0秒、Td=0.5秒を設定した場合、立上り時間tr≈4.1秒、整定時間ts≈8.2秒、オーバーシュート0%、IAE≈1.5となります。Kpを4.0に増加させると立上り時間がtr≈2.1秒、整定時間がts≈3.7秒に短縮され、IAEも約0.75まで低下します。むだ時間のない1次遅れ系では、Kpを上げるほど応答が速くなることが確認できます。
準拠/参考: 標準形(非干渉形)PID:u(t) = Kp[e + (1/Ti)∫e dτ + Td·de/dt](ISA / 制御工学の標準形)。Ziegler–Nichols ステップ応答法(FOPDT プラント K·e^(−Ls)/(τs+1)):Kp = 1.2τ/(K·L)、Ti = 2L、Td = 0.5L。
モデルの前提: 単一ループ SISO。プラントは一次遅れ+むだ時間/二次系/積分系を前進オイラー法で数値積分。微分項はノイズ抑制のためフィルタ付き(τ_f = Td/8)、積分項はアンチワインドアップでクランプ。検証では PI 制御で定常偏差ゼロ(目標 1.0 に対し 0.9999)を確認。
適用範囲・限界: 教育用の連続時間近似で、刻み dt=0.02 s 固定。アクチュエータ飽和は ±10 に固定し、計測ノイズ・量子化・離散化整定(Tustin 等)は扱わない。ZN 整定は出発点であり、実機では追い込みが前提。