過渡現象

複雑な回路の過渡現象をラプラス変換で解く!

皆さまお疲れさまです。ケンタ(@den1_tanaoroshi)です。

 

少し前に、ラプラス変換で過渡現象を解くと計算が捗るという記事を書きました。

電験2種 一次試験 理論 過渡現象 ラプラス変換
ラプラス変換で解く過渡現象皆さまお疲れさまです。ケンタ(@den1_tanaoroshi)です。 Twitterや2chでは合格・不合格の阿鼻叫喚となってい...

先日、ラプラス変換で解いてみると面白そうな過去問をtwitterで見つけましたので、今回はその解答例を紹介してみたいと思います。

過渡現象→微分方程式という固定観念

解説に入る前に、、、過渡現象は微分方程式でということが参考書の主流になっているので、その解説ももちろん微分方程式となることが殆どです。

今回の穴埋め問題は、誘導の流れに乗って微分方程式を使った方が答えやすいのかもしれませんが、問題によっては誘導文に微分方程式が出てこないことがあるので、微分方程式でもラプラス変換でもどちらを使っても構わないような時があります。それでも頑なに微分方程式で解説するのは、何故なんでしょう?

微分方程式で解く場合は立式も面倒ですし、特性方程式や指数関数を使うのも若干面倒です。一方でラプラス変換は、s回路図を作ってしまえばオームの法則っぽく立式できますし、最後のラプラス変換公式に式変形していくのが若干面倒ではありますすが、微積分を使わない分計算が少し楽に感じます。

ケンタ
ケンタ
ラプラス変換は自動制御だけという固定観念も有るような気がします。

まぁ、文句はこのくらいにして、解説を作っていきましょう!

解説の全体像

問題の誘導を気持ちいいくらいに完全に無視して、先ずはラプラス変換で一気にVC(t),i1(t),i2(t)を計算してから、穴埋め問題毎に形を整えて加工していきます。

理論問題を解くときは特にそうですが、私は誘導を完全に無視して、計算した結果から適宜解答を切り出すようなことを偶にします笑 去年の電験1種二次試験でもそうしていました。誘導に乗ると窮屈に感じるときがありますので。

稀に勝手に計算しすぎて、誘導に戻ってこれなくなる時がありますが、その時はその時です。

VC(t)の計算

IC(s)が分かると、VC(t),i1(t),i2(t)が芋づる式に分かっていきます。

 

先ずは、回路図をs回路図に変換してみましょう。

次に、条件式を立てます。

E1s=R1I1(s)+IC(s)Cs …式(1)

E2s=R2I2(s)+IC(s)Cs …式(2)

I1(s)+I2(s)=IC(s) …式(3)

 

立式する際に気にしておいてほしいこと

式を立てる上で気にしておいてほしいことに、未知数の数と式の数というのがあります。この2つの数が同じか式の数の方が多くないと、基本的に連立方程式は解けません。今回は未知数I1(s),I2(s),IC(s)の3つに対して、式が3つあるので計算が大変かもしれませんが、解けるということは漠然と分かります。

 

また、闇雲に式変形をするのではなく、未知数を消去していくという目算を立ててから式変形をしていって下さい。今回で言うと

 式(1) → I1(s)=(IC(s)の式) …式(1)’

 式(2) → I2(s)=(IC(s)の式) …式(2)’

としてから、式(3)に代入して

 式(3) → IC(s)だけの式に変形できる

    → VC(t),i1(t),i2(t)が分かる

というようにです。

問題に戻ります

式(1)より  I1(s)=1R1(E1sIC(s)Cs)

式(2)より  I2(s)=1R2(E2sIC(s)Cs)

を式(3)に代入して

1R1(E1sIC(s)Cs)+1R2(E2sIC(s)Cs)=IC(s)

1C(1R1+1R2)は表記が面倒ですので、時定数τを使って、1/τに置き換えます。

1s(E1R1+E2R2)=IC(s)(1+1τs)IC(s)=11+1τs1s(E1R1+E2R2)VC(s)=IC(s)Cs=1s(s+1τ)1C(E1R1+E2R2)=τ(1s1s+1τ)1C(E1R1+E2R2)L1τ(1etτ)1C(E1R1+E2R2)=τ1C(E1R1+E2R2)(1etτ)

 

式変形に絶対の自信があれば、τ1C(E1R1+E2R2)を定常値VCと置き替えてしまって構いませんが、念のため確認してみます。

式からは

 τ1C(E1R1+E2R2)=CR1R2R1+R21C(E1R1+E2R2)=R1R1+R2E1+R2R1+R2E2

 

一方で回路図からは、定常状態ではコンデンサに電流が流れていないので、

I=E1E2R1+R2

より、

VC=E1R1I=R1R1+R2E1+R2R1+R2E2

となるので、式変形に間違いが無かったことが確認されました。

 

よって、

VC(t)=VC(1etτ)

より、(1)は

τ=1C(1R1+1R2)

となります。

i1(t),i2(t)の計算

VC(s)=VC(1s1s+1τ)

を式(1)’を少し変形したところに代入して、

I1(s)=1R1(E1s+VC(s))=1R1(E1sVCs+VCs+1τ)=1R1(E1VCs+VCs+1τ)L11R1[(E1VC)+VCetτ]

となります。

 

I2(s)についても同様に、式(2)’から

I2(s)L11R2[(E2VC)+VCetτ]

となります。

 

以上から(2)は

i1(0)+i2(0)=1R1[(E1VC)+VC]+1R1[(E1VC)+VC]=E1R1+E2R2

となります。

 

また、(3)は

I1=limti1(t)=E1VCR1=1R1(R1R1+R2E1+R2R1+R2E2)=E1E2R1+R2

となります。

t=t0について

一般的にIC(t)=i1(t)+i2(t)となるので、IC(t0)=i1(t0)のときは、

IC(t0)=i1(t0)+i2(t0)i1(t0)=i1(t0)+i2(t0)i2(t0)=0

となります。

 

ここから必要な条件を抽出していきます。

i2(t0)=01R2[(E2VC)+VCet0τ]=0VCet0τ=VCE2

 

これをIC(t0)に代入して、(4)は

IC(t0)=i1(t0)+i2(t0)=1R1[(E1VC)+VCet0τ]+1R2[(E2VC)+VCet0τ]=1R1[(E1VC)+(VCE2)]+1R2[(E2VC)+(VCE2)]=E1E2R1

となります。

コンデンサへの充電について

最後です!

 

先ずは、コンデンサへ流れ込むIC(t)を求めます。

IC(t)=i1(t)+i2(t)=1R1[(E1VC)+VCetτ]+1R2[(E2VC)+VCetτ]=(1R1+1R2)VCetτ

よって(5)は

Δq=t0IC(t)dt=(1R1+1R2)VC[τetτ]t0=(1R1+1R2)VCet0ττ=(1R1+1R2)(VCE2)C11R1+1R2=(VCE2)C=CR2(E1E2)R1+R2

となります。

まとめ

微分方程式を使うかラプラス変換を使うかは、好みによるところが大きいです。

 

ただ、理想としてはどちらも使えるようにしておいたほうが良いです。今回のように微分方程式を前提とした誘導文だと、微分方程式を使って解いていったほうが省エネっぽそうですし。

手札は数が多い事に越したことはありませんからね。

 

それでは次回!

【2022.9.16更新】「電験2種 過渡現象をラプラス変換で解く27年間 2021度版」は好評発売中です!
「電験2種 過渡現象をラプラス変換で解く30年間 2025年版」のご紹介電験2種一次試験の過渡現象を煩雑な微分方程式ではなく、簡易なラプラス変換で解説しています。本書を理解することで、電験2種二次試験の自動制御だけではなく、電験1種一次試験の過渡現象にも対応できるようになります。...
こちらの記事もおすすめ!

COMMENT

メールアドレスが公開されることはありません。 が付いている欄は必須項目です

CAPTCHA