-- 有限要素法(3) (2009/10/09 00:48) - 2122PV

近似式が得られたところで、次は重み関数を考えます。

いきなりですが、今回使う重み関数を紹介します。



上が前回作った近似解。
下が今回設定した重み関数。

重み関数の条件としては、
・対象となる節点で1であること。
・それ以外の節点で0であること。
・全ての重み関数を足し合わせると、任意の点で1になることです。

最初の二つの条件については、グラフを見れば明らかです。
最後の条件については、グラフを見ても分からないです。

今回は節点が3つあります。
0とπ/4とπ/2です。
両端については境界条件で答えが出てるので考える必要ありません。

でも、両端にも重み関数があって、本当はそれを考慮する必要があります。
その重み関数も足し合わせると任意の点で1になります。

重み関数の理論とか、細かいことは置いておいて。
とりあえず、実際に計算してみることにします。



何を計算したら良いかと言えば。
解きたい偏微分方程式に、重み関数をかけたものでした。

このまま計算したら良いものですが、部分積分により少し展開します。



この必要性は不明。
そもそも不連続な関数を微分して良いのかが不明。
でも、これをしないと上手くいかない。



第1項はゼロになります。
なぜなら、重み関数は両端でゼロだからです。



結局のところ、この式を解けば良いということになります。
aの姿が見えませんが、これはaについての方程式です。

というわけで、とりあえず近似式u(x)と重み関数φ(x)の微分を求める必要がありそうです。
どちらも線形な関数なんで、xの係数をとってくるだけですが。



あとは、ぽいぽいぽいっと積分式に投げ込んで積分するだけです。



被積分関数は、要素ごとに定義されています。
そういうわけで、積分は要素ごとに分けて計算します。



そして、計算した結果がこちら。
計算途中をだらだらと書いても仕方ないので省略。
多項式の積分なので、高校生でもできるはず!

そしてもう、aの一次方程式になっているので、簡単に解けます。



これを実際に数値にしてみると、1.3882・・・ってところです。
厳密解は2の平方根で、1.4142・・・です。
まあまあ。



最後に厳密解と近似解をグラフで比較してみます。
節点が3つだけで、しかも線形という荒い近似でも、それらしい結果になりました。

精度を上げる方法は大きく分けて二通り。
・分割数を上げて、細かく計算する。
・直線でなく高次式で近似する。

ちなみに、今回は計算したい節点が1つでした。
だから、方程式は1本でした。
細かく分割すると、その数だけ方程式が出てきて、連立方程式になります。

実際に実用的な計算をしようと思うと、かなり巨大な連立方程式になります。
そして、その計算に、有限要素法の演算時間の大半をとられます。

ここまでで既に大変だけど、実際にはこの先がまた大変ということみたいです。

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

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

名前
コメント
更新情報
コメント「Pharmc378さん」(2018/10/16 01:33)
コメント「Pharmg834さん」(2018/10/16 01:31)
コメント「Pharmd230さん」(2018/10/16 01:31)
コメント「Pharme155さん」(2018/10/14 19:25)
コメント「Pharmd788さん」(2018/10/14 19:25)
コメント「Pharmk8さん」(2018/10/13 13:25)
コメント「Pharmg628さん」(2018/10/13 13:25)
コメント「Pharmf709さん」(2018/10/13 13:25)
コメント「Pharmb611さん」(2018/10/13 13:25)
コメント「Pharmb287さん」(2018/10/12 08:02)
コメント「Pharmf594さん」(2018/10/11 01:20)
コメント「Pharmb140さん」(2018/10/11 01:20)
コメント「Pharmg975さん」(2018/10/11 01:20)
コメント「Pharmg566さん」(2018/10/11 01:20)
コメント「Pharmk56さん」(2018/10/09 19:18)
コメント「Pharmg808さん」(2018/10/09 19:18)
コメント「Pharmk885さん」(2018/10/09 19:18)
コメント「Pharmb377さん」(2018/10/09 19:18)
コメント「Pharma706さん」(2018/10/08 13:15)
コメント「Pharmd718さん」(2018/10/08 13:15)
コメント「Pharme503さん」(2018/10/08 13:15)
コメント「Pharme612さん」(2018/10/08 13:15)
コメント「Johna911さん」(2018/10/07 20:26)
コメント「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)[1660]」 (2011/01/23 00:17)
電子工学「TwitCrew[2060]」 (2011/01/02 06:22)
コメント「かってぃさん」(2010/12/07 15:49)
コメント「通りすがりさん」(2010/12/02 23:37)
電子工学「ソーシャル目覚まし[2938]」 (2010/11/14 07:41)
その他「自己紹介[7675]」 (2010/10/16 01:36)
コメント「かってぃさん」(2010/10/15 02:03)
コメント「えんぺらあさん」(2010/10/15 02:00)
電子工学「性格の木[1816]」 (2010/10/15 01:43)
PR

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

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

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