初期温度0度の縦10mm、横20mmの長方形の左右の辺に、いきなり300度の温度を与えた与えた続けた場合の解析です。ansysではステップ加熱を選択しています。このansysでの解析計算のグラフとエクセルでの計算結果のグラフの比較したいのですがエクセルで使う計算式が分かりません。グラフは縦5mmの場所から横方向に0.5mm幅で温度を取りたいです。この時の計算式を教えて欲しいです。今回は非定常なので1秒ごとに解析を行い、1秒から10秒まで解析を行いました。材料特性は熱伝導率λ=50(W/m・K)、密度ρ=7300(kg/m^3)、比熱Cp=520(J/kg・K)です。
A 回答 (7件)
- 最新から表示
- 回答順に表示
No.7
- 回答日時:
面倒くさいので、スプレッドシートを作って公開しました。
計算結果だけならば、h=0.001m、1mm間隔で、時間は0.025秒の場合のものを、以下のリンク先で閲覧可能です。
https://docs.google.com/spreadsheets/d/1eyGcSD-- …
以下のアドレスをクリックすると、エクセル形式のファイルがダウンロードされます。
https://docs.google.com/spreadsheets/d/1eyGcSD-- …
わたしはエクセルとは違う表計算ソフトを使っているので互換性が完全ではないようで、フィルを変換した時に、計算結果のグラフの表示が変わってしまっているようですが、エクセルでも計算は可能だと思います。
対称性を考慮して計算させているので、計算領域は加熱面から中心部までになっています。
マクロは一切使っていないので、ウィルス感染の心配は無用かと思いますが、ファイルをダウンロードするときには覚悟してください。
それっぽい値が出ていますから、計算は間違っていないと思いますが・・・。
No.6
- 回答日時:
この問題は、両端が300℃で対称なので、この対称性を考慮すると、計算量が半分になります。
で、この対称性を考慮し、先のワークシートを少し変更し、1mm間隔、つまり、h=0.001mで計算した結果は次のようになるようです。
陽解法ですと、時間間隔Δtを小さくしないと、解が振動して大変なことになりますね(^^ゞ
参考にしてみてください。
それっぽい値が出ていますから、間違ってはいないと思います。
表計算ソフトの基本的な機能を使うだけでも、結構、いろいろなことができますね。
驚いております。
No.4
- 回答日時:
例えば、表計算ソフトで次のような表を作る。
ちなみに、B1からL1は計算点のx座標。
B2からL2までは時刻t=0の計算点xにおける温度。
B3以降は、時刻、時刻に応じた計算点の温度。
計算するときには、初期値として、B2、L2は300、C2〜K2は0を入力しておかないといけない。
また、B列のB2〜Bホニャララ、L列のL2〜Lホニャララまでには境界条件として300を事前にセットしておかないといけない。
計算式は、添付した写真に出ているけれど、
=(1-2*$O$7)*C2 +$O$7*(B2+D2)
これをC3からK3までにコピーぺして、
その後、B3からK3までを、B4〜K4までをBホニャララ〜Kホニャララまで計算できるはずです。
それらしい値が出ていますから、たぶん、間違ってはいないと思いますが、計算式は確かめてください。
添付した計算結果は、h=0.002mだから、2mm間隔ということになります。
h=0.0005だとセルの数が膨大になるので、表計算ソフトでの計算は厳しいかと思います。
この計算法ですと、θ=κΔt/h²<1/2にしないと、一般的に数値解が振動するので、θ<1/2におさめるようにΔtを設定し直す必要があります。
プログラミング言語を何かご存知でしたら、それを使ったほうが賢明かと思います。
先に紹介したものは陽解法と呼ばれるもので、連立方程式を解く陰解法を使うと、θ<1/2という制約条件はなくなり、無条件安定で計算できますが、このためには連立方程式を解くプログラムを作らないといけない。
No.3
- 回答日時:
ちょっと質問しますが、加熱していない両側は断熱ですよね。
ここで、例えば、対流や輻射などによって外部と熱のやりとりがあるわけじゃないですよね。
だとしたら、話がまったく違ってきます。
いま、ちょこっと表計算ソフトを計算してみたところ、計算は可能ですね。
簡単にできるので、チャレンジしてください。
わたしは、「教えてgoo」で回答をしないことにしているのですけれど、乗りかかった船ですから、最後まで面倒をみます。
わからないところがありましたら、遠慮せずに聞いてください。
No.2
- 回答日時:
この長方形の片側だけを加熱し、しかも、この長方形が半無限の場合は解析解が存在するけれど、この質問のケースについては解析的な解はないと思う。
あったとしても、その解は無限三角級数(フーリエ級数)の形で与えられるので、エクセルで関数の形で与えて解くことはできないと思う。温度が300度にどんと上がってから時間が短い場合、半無限の解析解を使って計算することはできるけれど、1秒とか10秒という長い時間の場合、この式は使えない。
ですから、エクセルを使ってこの問題を解こうとする場合、エクセルのプログラム記述用の言語VBAを使って、この偏微分方程式を数値的に解くことができるプログラムを自力で作り、その結果を使うということになる。
このためには、例えば、
「Excelで◯◯数値解析」
といった本を購入し、勉強をする必要がある。
こうした本には、1次非定常熱伝導方程式を解くプログラムが公開されていると思う。
ですが、
エクセルの表計算の機能を使って、陽解法で、この問題を解くことは可能。
以下のサイトなどが参考になると思います。
http://chemeng.in.coocan.jp/ExcelCe/exhc2.html
ワークシートまで公開されているので、数値を変更するだけで計算できるかもしれない。
ここで使っている解法は陽解法なので、ひょっとしたら、数値解が振動し、とんでもない結果になるかもしれないですね。
この場合、時間間隔を小さくして、数値的に解いた値が振動しないようにしないといけない。
この程度のワークシートならば、自力で簡単にできると思いますので、チャレンジしてみてください。
実は、もっと簡単に解く方法があるんだけどね。
お探しのQ&Aが見つからない時は、教えて!gooで質問しましょう!
似たような質問が見つかりました
- 中途・キャリア さて、社会人になってから13年目。 今まで、機械設計、機械製品の品質管理、部門の経営企画をやってきま 4 2023/07/17 08:59
- 化学 次の問題について聞きたいことがあります。 様々な温度(t)で五酸化二炭素(N2O5)の分解反応の反応 3 2023/06/16 15:35
- 物理学 風力発電での音 1 2023/04/16 08:55
- 数学 情報処理詳しい人!! A4縦のレポート文書に4:3の大きさの横向きの写真画像を貼り付けることにした。 2 2022/12/18 02:30
- 統計学 一変量分析(度数分布表)は、結果をもとに特に検定せずに断定してもよいですか? 7 2022/11/24 23:20
- 統計学 統計の課題なのですが、エクセルで1990年から2022年の8月の平均気温を分析したいのですがどんな表 2 2022/11/06 16:27
- Excel(エクセル) エクセルの大きなシートでグラフを見つける 4 2022/07/28 10:07
- 数学 数2Bの数列の問題です。 自分は、 まず数列 an=ar^(n-1)と置き こちらの問題の、y= の 1 2022/07/07 16:26
- 宇宙科学・天文学・天気 AIが答えた方程式 1 2023/02/20 00:12
- 統計学 統計学、エクセルがわかりません!解答と詳しい解説をお願いします! (1)それぞれの地域別に記述統計量 9 2022/08/21 16:30
このQ&Aを見た人はこんなQ&Aも見ています
おすすめ情報
このQ&Aを見た人がよく見るQ&A
デイリーランキングこのカテゴリの人気デイリーQ&Aランキング
-
勉強のことをかんがえる
-
マイクロブラックホール
-
状態ベクトルはなぜ線形空間の...
-
原子内の電子について
-
虚数時間と無限宇宙
-
遠心力は?
-
数学でこのMを横にしたのは、な...
-
高校物理解説お願いします!
-
S・ホーキングは世界的に有名で...
-
現時点での応力エネルギーテン...
-
LC回路の問題について 添付写真...
-
大学物理の問題について
-
電界中の電子運動の問題につい...
-
1次元とは?2次元とは?3次元...
-
エネルギー管理士の空気線図が...
-
タービンの内部効率って? エネ...
-
車の電気(電流)の流れ
-
空間以外のベクトル
-
波動関数の時間反転操作でなぜ...
-
化学で使う、水流ポンプ(アスピ...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
温かい水が下に行かない理由
-
刃物って当てただけでは切れま...
-
物質中の電子の速度の上限
-
1次元とは?2次元とは?3次元...
-
インダクターがエネルギーを一...
-
v^2-v0^2=2ax 今日この式を習っ...
-
I2Cのクロック信号を見たところ...
-
難しい理論を解りやすく説明す...
-
量子力学は物理学の最終形でし...
-
赤いブラックホール
-
ブラックホールを解明した
-
微分とは何か
-
電気は質量が無いので電力に周...
-
地球反ニュートリノ観測
-
アインシュタイン方程式
-
波動関数の時間反転操作でなぜ...
-
宇宙の中ではバグは出ないので...
-
コンデンサーは電圧降下 起電力??
-
エントロピーの減少
-
液体と気体の相違点は?
おすすめ情報
加熱していない上辺と下辺は断熱です。表計算ソフトでチャレンジしてみたんですが計算式するための式が分からず、とまっています。良かったら表計算ソフトで計算に使った式を教えて頂けたらうれしいです。
加熱していない上辺と下辺は断熱です。表計算ソフトでチャレンジしてみたんですが計算式するための式が分からず、とまっています。良かったら表計算ソフトで計算に使った式を教えて頂けたらうれしいです。