61次元熱伝導解析「Therm1D」-CQ出版社.pdfVIP

118 第 6 章 1次元熱伝導解析「Therm1D」 汎用的な1次元熱伝導解析ソフトウェアを添付しました．本ソフトウェアは，差分法により，1 次元の非定常熱伝導問題を解くものです． 1次元熱伝導解析は，1方向のみの熱伝導を取り扱うため，基礎的な熱伝導問題の検討に適して おり，入力データ作成の手間や計算時間も少なくて済みます． 6.1 差分法? 差分法は，対象を平行メッシュで有限個の差分要素に分割し，要素毎にその要素変数に関する 差分方程式を作成し，それを合わせて全体の方程式として解く解析方法です．以下では，1次元の 熱伝導に限定して説明します． （1）1次元熱伝導モデル 本ソフトウェアで取り扱うモデルを図6.1に示します．同図において，境界1と境界2の間に1次 元に熱伝導する複数の層があります．そして，その複数の各層はそれぞれ均等分割して差分方程 式を組み立てます．もちろん，複数層でなく単一層であっても構いません． また，同図に示すように，層を構成する材料の物性は，密度は固定で，比熱と熱伝導率は温度 区分により変えることができます．さらに，境界条件は，断熱，温度固定，熱伝達，熱流束，熱 放射の5種類があります． （2）差分法の基本計算式 上述のように，本ソフトウェアでは複数層の計算が可能ですが，差分法の基本計算式は複数層 図6.1●1次元熱伝導のモデル 複数層（各均等分割）? 境界1 境界2 ［材料物性］①密度（固定）? 　　　　　　②比熱（温度区分有り）? 　　　　　　③熱伝導率（温度区分有り）? ［境界条件］1）断熱　　　2）温度固定　　3）熱伝達? 　　　　　　4）熱流速　　5）熱放射? 6.1 差分法 119 図6.2●差分法の単一層モデル n 等分? 境界1 境界2 T0 T1 T2 T3 Tn 0 LX Δ X ［材料物性］? ρ：密度?c ：比熱?λ? ：熱伝導率? a ：温度伝導率（＝?λ?/ρ?c ）? でも単一層でも同様なので，図6.2の単一層のモデルで説明します． 同図において，境界1と境界2の間に1次元に熱伝導する単一層があり，同層はn等分されていま す．また，同図と重複しますが，次のように記号を定めます． Χ：X座標（熱伝導方向） ?Χ：メッシュ幅 Τ Τ i：各メッシュの中心温度（i =1, 2, ???, n） 0：境界の表面温度 ρ：密度 c：比熱 λ：熱伝導率　　a：温度伝導率（＝λ/ρc） さて，1次元の非定常熱伝導問題の支配方程式は，次式で表されます