![](http://oshiete.xgoo.jp/images/v2/pc/qa/question_title.png?8acaa2e)
初期温度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件)
- 最新から表示
- 回答順に表示
![](http://oshiete.xgoo.jp/images/v2/common/profile/M/noimageicon_setting_07.png?8acaa2e)
No.7
- 回答日時:
面倒くさいので、スプレッドシートを作って公開しました。
計算結果だけならば、h=0.001m、1mm間隔で、時間は0.025秒の場合のものを、以下のリンク先で閲覧可能です。
https://docs.google.com/spreadsheets/d/1eyGcSD-- …
以下のアドレスをクリックすると、エクセル形式のファイルがダウンロードされます。
https://docs.google.com/spreadsheets/d/1eyGcSD-- …
わたしはエクセルとは違う表計算ソフトを使っているので互換性が完全ではないようで、フィルを変換した時に、計算結果のグラフの表示が変わってしまっているようですが、エクセルでも計算は可能だと思います。
対称性を考慮して計算させているので、計算領域は加熱面から中心部までになっています。
マクロは一切使っていないので、ウィルス感染の心配は無用かと思いますが、ファイルをダウンロードするときには覚悟してください。
それっぽい値が出ていますから、計算は間違っていないと思いますが・・・。
![](http://oshiete.xgoo.jp/images/v2/common/profile/M/noimageicon_setting_07.png?8acaa2e)
No.6
- 回答日時:
この問題は、両端が300℃で対称なので、この対称性を考慮すると、計算量が半分になります。
で、この対称性を考慮し、先のワークシートを少し変更し、1mm間隔、つまり、h=0.001mで計算した結果は次のようになるようです。
陽解法ですと、時間間隔Δtを小さくしないと、解が振動して大変なことになりますね(^^ゞ
参考にしてみてください。
それっぽい値が出ていますから、間違ってはいないと思います。
表計算ソフトの基本的な機能を使うだけでも、結構、いろいろなことができますね。
驚いております。
![「非定常熱伝導解析のエクセルよるグラフの作」の回答画像6](http://oshiete.xgoo.jp/_/bucket/oshietegoo/images/media/2/1058352_58315d99b52f0/M.png)
![](http://oshiete.xgoo.jp/images/v2/common/profile/M/noimageicon_setting_07.png?8acaa2e)
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という制約条件はなくなり、無条件安定で計算できますが、このためには連立方程式を解くプログラムを作らないといけない。
![「非定常熱伝導解析のエクセルよるグラフの作」の回答画像4](http://oshiete.xgoo.jp/_/bucket/oshietegoo/images/media/f/1058352_58314a3bc492c/M.png)
![](http://oshiete.xgoo.jp/images/v2/common/profile/M/noimageicon_setting_07.png?8acaa2e)
No.3
- 回答日時:
ちょっと質問しますが、加熱していない両側は断熱ですよね。
ここで、例えば、対流や輻射などによって外部と熱のやりとりがあるわけじゃないですよね。
だとしたら、話がまったく違ってきます。
いま、ちょこっと表計算ソフトを計算してみたところ、計算は可能ですね。
簡単にできるので、チャレンジしてください。
わたしは、「教えてgoo」で回答をしないことにしているのですけれど、乗りかかった船ですから、最後まで面倒をみます。
わからないところがありましたら、遠慮せずに聞いてください。
![](http://oshiete.xgoo.jp/images/v2/common/profile/M/noimageicon_setting_07.png?8acaa2e)
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ランキング
-
ラプラス方程式
-
単位体積あたりについて
-
相対論的熱力学
-
大学理数科目の式変形、ゼロ除...
-
特殊相対性理論について
-
光において、波長がかなり違う...
-
量子化には二つの方法があると...
-
トンネルの中の質点
-
飛び降りた人はうつぶせになっ...
-
場の量子論での繰り込み処方に...
-
強い力は、核の外まで働いてい...
-
電磁気学の流束の問題
-
【高校物理】斜方投射の問題
-
ソレノイドコイルのコイル1巻あ...
-
熱電対を振ると温度が下がる理...
-
永久磁石の着磁後はどのように...
-
自然単位系の説明でc=ħ=1という...
-
重力が作用するものは、1つに繋...
-
極座標の運動方程式の計算の間...
-
摂動論での摂動項は厳密に解け...
マンスリーランキングこのカテゴリの人気マンスリーQ&Aランキング
-
熱湯500mlって、どうやって用意...
-
ローレンツ変換の説明で添字の...
-
この世で一番危険な放射性物質
-
物理 模試の問題
-
https://hooktail.sub.jp/mecha...
-
電子レンジでコップに入ったサ...
-
提灯や行灯って、蝋燭などの炎...
-
半球表面の問題
-
波動の光を視認する方法
-
トンネルの中の質点
-
Caの最外殻電子はなぜ2なのでし...
-
質量Mの物質を糸で吊るし、質量...
-
宇宙の広さってどこまであるん...
-
てこの計算問題
-
光において、波長がかなり違う...
-
電子レンジや冷蔵庫ってファラ...
-
倒れ始める瞬間では左下の角の...
-
余剰次元を減らすことは可能か?
-
二次元流れにおいて、流体のx、...
-
コンセントにあるアースターミ...
おすすめ情報
加熱していない上辺と下辺は断熱です。表計算ソフトでチャレンジしてみたんですが計算式するための式が分からず、とまっています。良かったら表計算ソフトで計算に使った式を教えて頂けたらうれしいです。
加熱していない上辺と下辺は断熱です。表計算ソフトでチャレンジしてみたんですが計算式するための式が分からず、とまっています。良かったら表計算ソフトで計算に使った式を教えて頂けたらうれしいです。