No.3ベストアンサー
- 回答日時:
もうだいぶ時間が経っているので解決したと思いますが,一応・・・
非定常熱伝導の問題で時間積分にクランクニコルソン法を
使うのは常套手段で非常に良い選択だと思いますが,
まずは,陽的オイラー法で行われてはいかがでしょうか.
時間刻みを充分小さくすれば全く問題無く解けます.
また,質量行列も集中化して陽的に解きます.
こうすれば,境界条件によって節点の温度ベクトル
に直接値を入れることで計算できます.
拡散係数行列をいじらなくとも計算できます.
一方,陰解法の場合は,係数行列の方を工夫します.
係数行列は温度の既知量を代入したところに対応する
係数行列の対角成分を「1」にして計算を行います.
このとき,非対角成分の行方向だけ0にする方法と
非対角成分の行と列の両方を0にするやり方があり
ます.
前者は非対称行列,後者は対称行列となります.
後者の方を使う場合には,温度境界ベクトルをさらに
工夫する必要があります.とりあえずやるならば,
前者の方をお勧めします.
回答ありがとうございます.
実はもう解決していたのですが締め切るのを忘れていました;)
僕は後者の,行と列を0にするやり方でやっています.
興味深いので陽的オイラー法でも試してみようと思います.
No.2
- 回答日時:
No1 ency です。
> すでにそのように計算しているのですがうまくいきません。
そうでしたか。。。
あとは、時間項の近似法とか、時間ステップをどのくらいの大きさでとっているのか、誤差が大きくなる理由としては、そのくらいしか思いつきません。。。
時間項の近似法には、陽解法、陰解法…他にもいろいろあったような気がします。
そのあたりも疑ってみてはいかがですか?
この回答への補足
近似法にはクランク・ニコルソン法を用いています。
今、許容できる最大刻み時間を計算してみたのですが、
その1/10ぐらいの刻み時間でやっているので問題はないと思います。
やはり原因は他にあるのかもしれません。
でも一応この質問は締め切らずにしばらく他の回答を待ってみようと思います。
No.1
- 回答日時:
非定常問題の場合、こんな感じで計算を進めるんじゃありませんでしたっけ?
(ステップ1)
境界条件として固定温度を、境界の節点温度として入力して計算する。
→結果として、各節点温度が出力される。
(ステップ2)
ステップ1 の結果を入力として計算する。
# 固定端は同じ条件を入力する。
→結果として、各節点温度が出力される。
……(以下繰返し)
有限要素法をやっていたのは、かれこれ十年近く前の話ですし、私自身は定常問題しかやったことはなかったので、見当違いなことを書いてしまっているかもしれません。
# 周りで非定常問題や、動的解析をやっている人がいたので、そのとき聞いた
# 話を思い出しつつ…なのもので。。。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 物理学 物理 7 2023/08/05 11:51
- 工学 電磁気学の質問です。 電流による電磁束密度ベクトルポテンシャルの計算ですが、ベクトルポテンシャルのポ 1 2022/04/19 17:23
- 数学 条件付き極値問題といわれる問題です。ラグランジュの乗数法 について、質問したいことがあります。 条件 3 2023/05/15 21:38
- 化学 化学の気体の法則について教えてください! 3 2023/06/02 22:26
- その他(ソフトウェア) F-BASICで計算中の実行が中途で勝手に止まり、大変困っています。 2 2023/03/02 16:15
- 弁護士・行政書士・司法書士・社会保険労務士 ★行政書士試験の行政法についての質問になります。 地方自治法についての質問になります。 問1 都道府 1 2023/06/18 09:00
- その他(職業・資格) 弁理士試験の勉強方法について 1 2022/09/11 07:32
- 発達障害・ダウン症・自閉症 中学の時にIQ82の境界知能と診断されました。 今の私も、やはり境界知能でしょうか? そしてこれは、 3 2023/02/19 00:37
- 数学 『4色問題③』 2 2022/11/14 00:31
- 工学 PID制御 1 2022/07/27 11:50
おすすめ情報
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
おすすめ情報