らい・ぶらり資料請求番号：TS36 TS41 TS91
一般化した複雑な電気回路や機械力学の微分方程式をラプラス変換で解く
こちら↓の記事でラプラス変換の基礎的な説明と簡単な演習を行いました。
今回は、ラプラス変換による微分方程式の求解に慣れてきた人のために、より一般化した、すなわち複雑な微分方程式の例題を取り扱っていきたいと思います。
RLC直列回路のステップ応答
こちらの記事ではRL直列回路を取り扱いましたが、今回はコンデンサCも含む、RLC直列回路を取り扱います。回路方程式は以下の通りです。
さらに、今回はステップ関数Θを使用し、ある時刻t=0において、回路のスイッチを入れたとき、その応答はどのようになるか？を計算します。
このようなステップや、後ほど取り扱うインパルスに対する過渡応答の解析と言うのはラプラス変換の得意とするところです。
上の式を微分します。ステップ関数は微分するとインパルス関数になります。
さらに、両辺をLで割ったこの式をラプラス変換していきます。
まず、①をラプラス変換すると以下のようになります。（※ここで？になった方はこちらをチェック！）
求めるべき関数i(t)のラプラス変換をI(s)とおき、
回路のスイッチを入れる瞬間t=0では電流が流れていないことから
よって
同様に他の項もラプラス変換すると、以下のような代数方程式が得られます。
代数方程式を解きます。
このとき、
と置きました。
RLC直列回路の過渡応答は
と置いたとき、b2の値によってその解の形が変化しますので、
b2>0, b2=0, b2<0に場合分けして
ラプラス逆変換を行い、解を求めます。
b2＞0の場合
をラプラス逆変換します。変換表
出典：斎藤制海・徐 粒 共著，計測と制御シリーズ 制御工学 －フィートバック制御の考え方ー，第１版第１刷，2003年，pp.19
より、
であるから、記号を置き換えて
を得ます。
従って、求める電流の時間変化i(t)は
となります。
b2＝0の場合
より、
代数方程式の解は以下のように変形できます。
ラプラス変換表から
となるので、
ラプラス逆変換を行い、以下の解を得ます。
b2＜0の場合
代数方程式の解を以下のように変形します。
さらに、
と置いているので、
これをラプラス逆変換します。
ラプラス変換表には記載がありませんが、
このような公式がありますので、今回のケースに当てはめて考えると
これによりラプラス逆変換を行い、解を得ます。
機械力学（バネーダッシュポット系のインパルス応答）
次は、以下の方程式をラプラス変換によって解くことを考えます。
これは、質量mの物体にバネを取り付け、その物体に一瞬だけF/mω0なる刺激を与えたら物体はどのように運動するか？という問題です。一瞬の刺激ではなく、持続的に力を加えた場合の応答は
で表現でき、この微分方程式を↓こちら↓
のページでルンゲクッタ法により解きました。
δ(t)をインパルス関数と言い、ある時間t=t0の時だけ値が∞になって、あとはゼロになる関数です。この関数はラプラス変換をとると1になります。
従って、
として、初期条件
のもと、それぞれの項をラプラス変換すると
となり、代数方程式
を得ます。これを解くと
となりますので、これをラプラス逆変換します。
変換表
出典：斎藤制海・徐 粒 共著，計測と制御シリーズ 制御工学 －フィートバック制御の考え方ー，第１版第１刷，2003年，pp.19
より、
であるので、これを今回得られた代数方程式に合うよう変形していきます。
このとき、分母のγ2がが余計になるので、ω0を√ω02–γ2と置き換えて、γ2を消去します。
これにより、
の関係を得られるので代数方程式
の逆ラプラス変換が可能になります。
これにより、バネーダッシュポット系の要素に対し、インパルスを与えたときの応答を求めることができました。
コメントを残す