-- 有限要素法 (2009/10/08 00:01) - 2242PV

有限要素法でヘルムホルツ方程式を解いてみます。



復習ですが、これがヘルムホルツ方程式。
変数自体と、その2階微分を含む微分方程式です。



簡単にするために、1次元に落としました。
3次元の有限要素法を説明しようと思ったら、それこそ本になります。

有限要素法に使うのは、「近似解」と「重み関数」です。
近似解が欲しくて有限要素法を使おうとしているのに、近似解を要求されるなんて。
それは何故か?

まず、微分方程式じゃなくて、ただの方程式の数値解析法を思い出してみます。
分からないなら、エクセルのゴールシークを想像します。
それも分からないなら、電卓だけで方程式を解く方法を想像してみます。



まず、解はこれくらいかなっていう「適当な数値」を方程式に代入してみます。
答えがぴったりゼロになれば、それは厳密に方程式の解です。
そうでなければ、違う数値を入れてみて計算を繰り返します。

ぴったりでなくとも、0.0001くらいになれば、十分に良い近似解と言えます。



同様に、微分方程式を解くときも、まず「適当な近似解」を用意します。
それを微分方程式に代入してみて、0に近い値になれば解と認められます。
そういうわけで、まずは「とりあえずの近似解」を準備するわけです。

ただ、色々な近似解を代入して試してみるわけではありません。
定数を含む近似解を作って、その定数を定めることで、解を見つけ出します。
「色々な」の部分を、色々な数値に変化できる「定数」に肩代わりさせる様なイメージ。

そして、その定数の定め方がこれ。



近似解に重み関数をかけて、領域積分し、ゼロと置く。
これを計算すると求めたい定数を含む方程式が得られます。
その方程式を解いて定数を定めれば解が得られる・・・というわけ。

さて、話は少し戻って、「重み付け残差法」の「残差」とは何か?
「とりあえずの近似解」を微分方程式に代入してみたらどうなるかといえば。



この式に限って、解u(x)が近似式であることを強調するために、ダッシュをつけています。

「とりあえずの近似解」が厳密解にぴったり一致していない限り、残差が残りますよね。
だから、残差です。

繰り返しになりますが、この残差を小さくするように、近似解中の定数を定めます。
その定め方は、重み関数をかけて、領域積分して、ゼロと置くことです。

これが重み付け残差法!

-- 参考文献
よくわかる有限要素法

-- 関連項目
1階線形常微分方程式
1階線形常微分方程式(2)
1階線形常微分方程式(3)
1階線形常微分方程式(4)
抵抗器のモデリング
よくわかる有限要素法
微分方程式の数値解析
ヘルムホルツ方程式
有限要素法(2)
有限要素法(3)
静電場の支配方程式
有限要素法の精度向上法
有限要素法の精度向上法(2)
MATLABクローン

名前
コメント
更新情報
コメント「Johnf318さん」(2015/11/30 06:26)
コメント「もちもちさん」(2015/10/15 01:25)
コメント「しもさんさん」(2014/07/11 12:57)
コメント「ljilchさん」(2014/06/28 02:40)
コメント「ledogoさん」(2013/06/25 23:11)
コメント「tcrwefcxtjさん」(2013/06/25 21:00)
コメント「gihjfldrwさん」(2013/06/23 08:07)
コメント「yfaewyqviさん」(2013/06/23 03:12)
コメント「Papan3世さん」(2013/05/09 11:48)
コメント「Papan3世さん」(2013/05/09 11:46)
コメント「Papan3世さん」(2013/05/09 11:45)
コメント「Papan3世さん」(2013/05/09 11:40)
コメント「Papan3世さん」(2013/05/09 11:37)
コメント「さん」(2013/04/21 01:57)
コメント「上田朋哉さん」(2013/04/08 02:43)
コメント「syazaqxさん」(2013/02/14 11:36)
コメント「seqdspさん」(2012/11/08 02:34)
コメント「kzghepさん」(2012/11/03 21:45)
コメント「さん」(2012/10/16 08:05)
コメント「yflpfhftbvpさん」(2012/07/21 14:51)
コメント「szodydさん」(2012/06/22 15:38)
コメント「nkhkajuttzさん」(2012/06/20 03:23)
コメント「jjomlllxさん」(2012/02/13 12:34)
コメント「さん」(2011/10/28 14:11)
コメント「onmqdqobrpdさん」(2011/10/15 03:57)
コメント「さん」(2011/08/20 14:05)
コメント「さん」(2011/06/26 20:59)
コメント「bfixkuohhzoさん」(2011/06/10 00:04)
コメント「さん」(2011/05/24 13:47)
コメント「dyecrnさん」(2011/05/12 20:06)
コメント「atkyzmさん」(2011/03/02 07:59)
日記「新年を迎えまして(4)[1632]」 (2011/01/23 00:17)
電子工学「TwitCrew[2038]」 (2011/01/02 06:22)
コメント「かってぃさん」(2010/12/07 15:49)
コメント「通りすがりさん」(2010/12/02 23:37)
電子工学「ソーシャル目覚まし[2899]」 (2010/11/14 07:41)
その他「自己紹介[7633]」 (2010/10/16 01:36)
コメント「かってぃさん」(2010/10/15 02:03)
コメント「えんぺらあさん」(2010/10/15 02:00)
電子工学「性格の木[1794]」 (2010/10/15 01:43)
PR

プロフィール
名前: かってぃ
職業: 学生、性別: 男性、年齢: 22歳、住所: 横浜市
理系用語
コンテンツ
トップページ
自己紹介
電子工学

料理
勉強
日記
ノート
自転車
その他
本棚
一言メモ
電子回路コレクション
ヘッドホンアンプのまとめ
リンク
三度の飯とエレクトロン
choilog
TwitCrew
Questo
Ropi.me
studylog
つい、投票。
tvtw - ツイッターでテレビ実況
理系用語
理系用語bot on twitter

ぷよぷよ!!クエスト 攻略
気ままでいたいがあんまり気ままではない感じのブログ
魚料理と簡単レシピ
素晴らしき哉,電気工学的人生!
アンチ・バベルの塔
大学生ブログランキング