形状係数
理論と物理
形状係数（Shape Factor）とは
先生、形状係数って何に使うんですか？
2次元や3次元の定常熱伝導問題を、1次元の熱抵抗に帰着させるための概念だ。形状係数 $S$ を使えば
と表現でき、$\Delta T$ から直接放熱量が分かる。熱抵抗は $R = 1/(kS)$ だ。
複雑な形状を1つの数値に集約できるんですね。
$S$ は形状と境界条件だけで決まる幾何学的な量で、単位は [m]（2D問題では [m/m] = 無次元/単位奥行）だ。
代表的な形状係数
| 形状 | $S$ | 適用条件 |
|---|---|---|
| 無限平板 | $A/L$ | 基本形 |
| 同心円筒 | $2\pi L / \ln(r_2/r_1)$ | 長さ $L$ |
| 同心球 | $4\pi r_1 r_2/(r_2 - r_1)$ | — |
| 埋設球（半無限体表面から深さ$z$） | $4\pi r / (1 - r/(2z))$ | $z > r$ |
| 埋設円筒（半無限体） | $2\pi L / \cosh^{-1}(z/r)$ | $z > r$, $L \gg r$ |
| 2つの平行円筒 | $2\pi L / \cosh^{-1}((d^2-r_1^2-r_2^2)/(2r_1 r_2))$ | 中心間距離 $d$ |
埋設球や埋設円筒の式が実用的そうですね。
地中埋設配管の放熱量計算、建物基礎からの地中伝熱、地熱ヒートポンプの設計などで使われる。手計算で概算が得られるのが形状係数の最大の利点だ。
形状係数の導出
ラプラス方程式 $\nabla^2 T = 0$ の解から導出する。同心円筒の場合
よって $S = 2\pi L / \ln(r_2/r_1)$ が得られる。
解析解から形状係数を逆算するんですね。
解析解が得られない形状ではFEMで数値的に $S = q/(k\Delta T)$ を算出する。
形状係数の定義と物理的意味
形状係数Sは複雑形状での熱流量をq=SkΔTで表す無次元係数だ。1950年代にCarslaw & Jaegerが「固体の熱伝導」で体系化し、埋設配管から球体・円柱まで60種以上の解析解が表として整備された。
各項の物理的意味
- 蓄熱項 $\rho c_p \partial T/\partial t$：単位体積あたりの熱エネルギー蓄積率。【日常の例】鉄のフライパンは熱しにくく冷めにくいが、アルミ鍋は熱しやすく冷めやすい——これは密度 $\rho$ と比熱 $c_p$ の積（熱容量）の違い。熱容量が大きい物体は温度変化が緩やかになる。水は比熱が非常に大きい（4,186 J/(kg·K)）ため、海沿いの気温は内陸より安定する。非定常解析ではこの項が温度の時間変化速度を決める。
- 熱伝導項 $\nabla \cdot (k \nabla T)$：フーリエの法則に基づく熱伝導。温度勾配に比例した熱流束。【日常の例】金属スプーンを熱い鍋に入れると持ち手まで熱くなる——金属は熱伝導率 $k$ が高いため、高温側から低温側へ素早く熱が伝わる。木製スプーンが熱くならないのは $k$ が小さいから。断熱材（グラスウール等）は $k$ が極めて小さく、温度勾配があっても熱が伝わりにくい。「温度差のあるところに熱が流れる」という自然の傾向を数式化したもの。
- 対流項 $\rho c_p \mathbf{u} \cdot \nabla T$：流体の運動に伴う熱輸送。【日常の例】扇風機に当たると涼しく感じるのは、風（流体の流れ）が体表面近くの暖かい空気を運び去り、新鮮な冷たい空気を供給するから——これが強制対流。暖房で部屋の天井付近が暖かくなるのは、暖められた空気が浮力で上昇する自然対流。PCのCPUクーラーのファンも強制対流で放熱している。対流は熱伝導よりも桁違いに効率的な熱輸送手段。
- 熱源項 $Q$：内部発熱（ジュール熱、化学反応熱、放射線吸収等）。単位: W/m³。【日常の例】電子レンジは食品内部のマイクロ波吸収（体積発熱）で加熱する。電気毛布のヒーター線はジュール発熱（$Q = I^2 R / V$）で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
数値解法と実装
数値的な形状係数の算出
複雑な形状の形状係数はどうやって求めるんですか？
FEMで定常熱伝導を解き、以下の手順で $S$ を算出する。
1. 高温面に $T_1$、低温面に $T_2$ のDirichlet条件を設定
2. $k = 1$ W/(m K) に設定（簡略化のため）
3. 解析を実行し、高温面（または低温面）の全熱流量 $q$ を取得
4. $S = q / (k \cdot \Delta T) = q / (T_1 - T_2)$ を算出
$k = 1$ にするのは計算を楽にするためですね。
そう。$S$ は $k$ に依存しない幾何量なので、$k$ の値は結果に影響しない。
メッシュ収束の確認
形状係数は積分量（全熱流量）なので局所的なメッシュ感度は小さいが、形状が複雑な場合は収束確認が必要だ。
| メッシュレベル | 要素数 | $S$ | 誤差 |
|---|---|---|---|
| 粗い | 1,000 | 15.2 m | — |
| 中間 | 10,000 | 15.8 m | 3.9% |
| 細かい | 100,000 | 15.9 m | 0.6% |
| 非常に細かい | 1,000,000 | 15.9 m | 0.0% |
積分量は比較的早く収束しますね。
温度の積分量は応力の局所値より収束が速い。1万要素でも実用精度が得られることが多い。
重ね合わせの原理
形状係数は線形問題なので重ね合わせが可能だ。複数の熱源がある場合、各熱源からの形状係数を足し合わせられる。また対称条件を使えば計算量を削減できる。
対称面で半分のモデルにして $S$ を2倍にすればいいんですね。
その通り。地中埋設円筒の場合、地表面を対称面（断熱面）として扱うと、半無限体から有限領域に問題を縮小できる。
形状係数の解析解一覧の使い方
地中埋設管（直径D、深さz、長さL）の形状係数はS=2πL/ln(4z/D)（z>>D時）。例えばD=100mm、z=1mの東京都水道管20mではS≈27m、k=1.5 W/m·Kの土壌でq≈243WのΔT毎の熱流量が得られる実用計算だ。
線形要素 vs 2次要素
熱伝導解析では線形要素でも十分な精度が得られることが多い。温度勾配が急な領域（熱衝撃等）では2次要素を推奨。
熱流束の評価
要素内の温度勾配から算出。節点応力と同様にスムージングが必要な場合がある。
対流-拡散問題
ペクレ数が高い（対流支配）場合、風上的安定化（SUPG等）が必要。純粋な熱伝導問題では不要。
非定常解析の時間刻み
熱拡散の特性時間 $\tau = L^2 / \alpha$（$\alpha$: 熱拡散率）に対して十分小さい刻みを設定。急激な温度変化には自動時間刻み制御が有効。
非線形収束
温度依存物性値による非線形性はマイルドな場合が多く、Picard反復（直接置換法）で十分なことが多い。放射の強非線形性ではニュートン法を推奨。
定常解析の判定
全節点の温度変化が閾値以下（$|\Delta T| / T_{max} < 10^{-5}$等）で収束と判定。
陽解法と陰解法のたとえ
陽解法は「今の情報だけで次を予測する天気予報」——計算は速いが大きな時間刻みでは不安定（嵐を見逃す）。陰解法は「未来の状態も考慮した予測」——大きな時間刻みでも安定するが、各ステップで方程式を解く手間がかかる。急激な温度変化がない問題では陰解法で大きな時間刻みを使う方が効率的。
実践ガイド
地中埋設配管の放熱計算
形状係数が最も活用される場面は何ですか？
地中埋設配管の放熱量計算だ。地中の温度場は半無限体として扱え、形状係数の公式がそのまま使える。
計算例
蒸気配管（外径114.3mm、断熱材外径214.3mm）を地表から1.5mの深さに埋設。地中温度15℃、断熱材外面80℃。地盤の $k_{\text{soil}} = 1.5$ W/(m K)。
1mあたり184Wの放熱ですか。
100mの配管なら18.4kWの熱損失だ。年間のエネルギーコストに換算して断熱厚の経済最適化を行う。
建物基礎の地中伝熱
スラブオングレード（土間コンクリート床）からの地中伝熱もISO 13370で形状係数ベースの計算が規定されている。
| パラメータ | 影響 |
|---|---|
| 基礎面積/周長比 | 大きいほど地中への放熱が少ない |
| 断熱材配置 | 基礎外周部の断熱が最も効果的 |
| 地盤の $k$ | 砂 1.5、粘土 1.0、泥炭 0.5 W/(m K) |
大きな建物ほど有利なんですね。
面積が $L^2$ に比例し周長が $L$ に比例するため、大きな建物は周長あたりの面積が大きく、地中伝熱の影響が相対的に小さくなる。
形状係数の参考文献
形状係数の網羅的な一覧はIncropera「Fundamentals of Heat and Mass Transfer」のTable 4.1、またはBejan「Heat Transfer」のAppendixに収録されている。特殊形状はHahne & Grigullの文献が充実している。
教科書に載っていない形状はFEMで求めるしかないですね。
その通り。一度FEMで求めた $S$ は設計式としてデータベース化しておけば、以降は手計算で使える。
地中埋設配管の熱損失評価
北海道の地域熱供給では地中80cmに埋設する100A鋼管から年間熱損失が約15 W/mに達する。形状係数を使った解析で断熱被覆の最適厚さを30mmと算出し、ポリウレタンフォーム（k=0.027 W/m·K）の採用で損失を5 W/m以下に低減した事例がある。
解析フローのたとえ
熱解析のフローは「お風呂の追い焚き設計」で考えてみましょう。浴槽の形（解析対象）を決め、お湯の初期温度（初期条件）と外気温（境界条件）を設定し、追い焚きの出力（熱源）を調整する。「2時間後にぬるくなっていないか？」を計算で予測する——これが非定常熱解析の本質です。
初心者が陥りやすい落とし穴
「放射を無視していいですか？」——室温付近なら大抵OK。でも数百度を超えたら話は別です。放射による熱伝達は温度の4乗に比例するため、高温では対流を圧倒します。晴れた日に日向と日陰で体感温度が全然違うのを経験したことがありますよね？ あれが放射の威力です。工業炉やエンジン周りの解析で放射を無視するのは、猛暑日に「日差しは関係ない」と言い張るようなものです。
境界条件の考え方
熱伝達係数 $h$ は「窓の断熱性能」だと思ってください。$h$ が大きい＝窓が薄い＝熱がどんどん逃げる。$h$ が小さい＝二重窓＝熱が逃げにくい。この数値1つで結果が大きく変わるため、文献値の引用や実験による同定が重要です。「とりあえず10 W/(m²·K)で…」と適当に入れていませんか？
ソフトウェア比較
商用ツールでの形状係数算出
形状係数を商用ツールで求めるにはどうしますか？
どのFEMツールでも基本手順は同じだ。
| ツール | 手順 | 熱流量の取得方法 |
|---|---|---|
| Ansys Mechanical | T固定面のReaction Heat Flowを取得 | FSUM, PRRSOL |
| Abaqus | Output > Node > RFL11（反力熱流量） | History Output |
| COMSOL | Surface Integration > Total Heat Flux | Derived Values |
APDL実装例
埋設円筒の形状係数を数値的に求めるAPDLコード。
```
/PREP7
ET,1,PLANE55,,,1 ! 軸対称ではなく平面
MP,KXX,1,1.0 ! k=1（形状係数算出用）
! 半無限体を有限領域で近似
RECTNG,0,5,0,-3 ! 5m幅 × 3m深さ
CYL4,0,-1.5,0.107,,,, ! 埋設円筒（深さ1.5m, r=107mm）
ASBA,1,2 ! 差分
ESIZE,0.02
AMESH,ALL
/SOL
D,円筒面NODE,,80 ! 管外面80℃
D,地表面NODE,,15 ! 地表15℃
! 側面・底面は断熱（デフォルト）
SOLVE
*GET,Q,FSUM,,HEAT ! 全熱流量
S_factor = Q / (1(80-15)) ! S = q/(kDT)
```
$k = 1$ にしてるから $S = Q/\Delta T$ ですね。
理論値 $S = 1.89$ m/m と比較して検証する。有限領域の大きさが十分か（側面の影響がないか）も確認が必要だ。
地熱設計専用ツール
地熱ヒートポンプの設計では専用ツールが使われる。
| ツール | 用途 |
|---|---|
| GLD (Ground Loop Design) | ボアホール熱交換器の設計 |
| EED (Earth Energy Designer) | 地中熱交換器の長期性能予測 |
| COMSOL Geothermal | 3D地中温度場のFEM解析 |
| FEFLOW | 地下水流を考慮した地中伝熱 |
地下水の流れも影響するんですね。
地下水流がある場合、対流による熱輸送で形状係数が大幅に変わる。純粋な伝導問題ではなくなるため、FEFLOWやCOMSOLでの数値解析が必要になる。
形状係数計算のCAEツール活用
Autodesk CFD 2025は地中埋設構造の形状係数を自動的にポスト処理で出力する機能を持つ。OpenFOAMのlaplacianFoamソルバーで定常熱伝導を解いた後、surfaceIntegrateコマンドで熱流束を積算してS値を算出するワークフローも確立されている。
選定で最も重要な3つの問い
- 「何を解くか」：形状係数に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
複合形状の形状係数
教科書に載っていない複雑な形状はどう扱うんですか？
複合形状を簡単な形状に分解して重ね合わせる方法がある。L字型の壁はストレート部分とコーナー部分に分解し、各部の形状係数を足し合わせる。
コーナーの形状係数はコーナー角度と壁厚で決まる補正値で、教科書に表がある。
パズルみたいに分解するんですね。
分解が難しい場合はFEMで直接算出する。一度求めた $S$ を設計式にフィッティングすれば、パラメトリックな設計計算に使える。
非定常問題への拡張
非定常問題では形状係数が時間の関数 $S(t)$ になる。地中埋設配管の起動時の過渡応答は
$\alpha$ は地盤の温度拡散率、$\gamma = 0.5772$ はオイラー定数。長時間後に定常の $S$ に漸近する。
地中の温度場が定常に達するまでどのくらいかかりますか？
配管サイズと地盤の熱拡散率による。典型的には数日〜数週間。地熱ヒートポンプでは数年かかることもある。
異方性媒体の形状係数
地盤が層状構造（水平方向と鉛直方向で $k$ が異なる）の場合、座標変換で等方性問題に帰着させる。
変換後の等方性媒体（$k_{\text{eff}} = \sqrt{k_x k_y}$）で形状係数を求め、元の座標系に戻す。
座標をスケーリングして等方性に変換するのは賢いですね。
地質の層状構造はまさにこのケースだ。水平方向の透水層は鉛直方向より $k$ が大きいことが多い。
数値形状係数の抽出手法
解析解のない複雑形状の形状係数はFEMで数値的に抽出できる。熱流量qをFEMで計算後、S=q/(kΔT)で逆算する手法をFENICS Projectで実装すると、任意の3D形状に対して1時間以内にSを取得できる。
トラブルシューティング
よくあるトラブルと対策
形状係数の計算で注意すべきことは何ですか？
典型的な問題を整理する。
1. 有限領域サイズの影響
問題: 半無限体をFEMで解くとき、有限領域で打ち切る。領域が小さすぎると $S$ が過小評価される。
対策: 領域の境界を関心領域の10倍以上離す。境界サイズを2倍にして $S$ が1%以内で変わらなければ十分だ。
2. 形状係数の適用条件の逸脱
教科書の公式をそのまま使って間違えるケースですね。
よくある誤用を挙げる。
| 公式 | 適用条件 | よくある誤用 |
|---|---|---|
| 埋設円筒 | $L \gg r$, $z > r$ | 短い管（$L/r < 10$）に適用 |
| 埋設球 | $z > r$ | 地表面に近い球（$z \approx r$）に適用 |
| 2平行円筒 | $d > r_1 + r_2$ | 円筒が接触する場合に適用 |
3. 地盤物性の不確実性
地中伝熱の最大の不確実性は地盤の熱物性だ。
| 地盤 | $k$ [W/(m K)] | 変動幅 |
|---|---|---|
| 乾燥砂 | 0.3 | $\pm$50% |
| 湿潤砂 | 2.0 | $\pm$30% |
| 粘土 | 1.0 | $\pm$40% |
| 岩盤 | 2.5 | $\pm$20% |
地盤の $k$ は含水率で数倍変わるんですね。
TRT（Thermal Response Test）で現地計測するのが最も確実だ。ボアホールに熱負荷を与えて温度応答から $k$ を逆算する。地熱ヒートポンプ設計では必須の工程だ。
4. 重ね合わせの限界
形状係数の重ね合わせは線形問題でのみ有効。放射を含む場合や温度依存 $k$ の場合は重ね合わせが使えない。
高温問題では要注意ですね。
その場合はFEMで直接解くしかない。形状係数は概算ツールであり、限界を理解して使うことが重要だ。
形状係数表の適用条件ミス
Carslaw & Jaegerの形状係数表には各解析解の有効条件（z/D比など）が明記されているが、範囲外適用によって誤差50%超になることがある。地中埋設管でz/D<2の浅い場合は表の式が使えず、Moody補正式か直接FEM解析が必要だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——形状係数の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告