複合材料の熱伝導
理論と物理
複合材料の熱伝導特性
CFRPやGFRPのような複合材料の熱伝導って、均質材料とどう違うんですか？
繊維強化複合材料は強い異方性を持つ。繊維方向の熱伝導率は繊維の寄与が大きく、直交方向はマトリクス（樹脂）が支配する。結果として方向によってkが10倍以上異なることがある。
有効熱伝導率の理論
繊維方向（並列モデル）と直交方向（直列モデル）の有効熱伝導率はこうなる。
ここで $V_f$ は繊維体積率、$k_f$ は繊維の熱伝導率、$k_m$ はマトリクスの熱伝導率だ。
具体的な数値でいうとどのくらいの差ですか？
PAN系炭素繊維（$k_f=10$ W/(mK)）、エポキシ樹脂（$k_m=0.2$ W/(mK)）、$V_f=0.6$ の場合:
- $k_{\parallel} = 0.6 \times 10 + 0.4 \times 0.2 = 6.08$ W/(mK)
- $k_{\perp} \approx 0.45$ W/(mK)
13倍以上の差がある。これを無視して等方性で解析すると温度分布が全く変わる。
Hashin-Shtrikmanの限界
より精密な評価にはHashin-Shtrikmanの上下界を使う。
実測値はこの上下界の間に入る。Halpin-Tsaiモデルや有限要素均質化（RVE解析）でさらに精密な予測が可能だ。
繊維配向がランダムな短繊維材料だとどうなりますか？
ランダム配向の場合は等方的に近づくが、射出成形品では流動方向に繊維が配向するため部分的な異方性が生じる。Moldflow等の射出成形シミュレーションから繊維配向テンソルをエクスポートして、熱伝導率テンソルにマッピングする手法が実用化されている。
複合材料の法則、1850年代から
Maxwell（1873年）が球状粒子を分散させた複合材の等価熱伝導率を初めて理論化した。彼の式は今日でもポリマー基板に銅粒子を混ぜたTIM（熱界面材料）の設計に使われ、λeff ≈ λm(λp+2λm+2φ(λp−λm))/(λp+2λm−φ(λp−λm))という形で教科書に残る。
各項の物理的意味
- 蓄熱項 $\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$）で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
数値解法と実装
RVE解析による均質化
繊維レベルまでモデル化して解析するんですか？
代表体積要素（RVE）を用いたマルチスケール解析が標準手法だ。繊維直径7umの炭素繊維を六角配列で並べたRVEを作成し、各方向に温度差を与えて有効熱伝導率テンソルを求める。
RVEのサイズはどう決めるんですか？
繊維直径の10〜20倍が目安だ。RVEサイズを大きくして結果が変わらなければ十分。COMSOLやDigimatでは自動RVE生成とパラメトリック均質化が可能だ。
積層板のモデリング
積層板（例: [0/90/45/-45]s）は各プライの熱伝導率テンソルを積層方向に回転させて重ね合わせる。AbaqusではORIENTATIONとSHELL SECTIONで各層の配向角を定義する。
| 積層構成 | 面内k [W/(mK)] | 面直k [W/(mK)] |
|---|---|---|
| UD [0]8 | 6.0 / 0.45 | 0.45 |
| クロス [0/90]2s | 3.2 / 3.2 | 0.45 |
| 準等方 [0/45/90/-45]s | 3.2 / 3.2 | 0.45 |
クロスと準等方は面内が均一になるんですね。
そうだ。ただし面直方向はどの構成でもマトリクス支配で低いままだ。面直方向の熱経路確保がCFRP構造の熱設計で最大の課題になる。Zピンやカーボンナノチューブの面直導入で改善する研究が進んでいる。
並列・直列則の精度の違い
繊維強化樹脂（CFRP）の面内熱伝導率は並列則（混合則）で±5%以内に収まるが、厚さ方向は直列則でも誤差20%超になることが多い。1990年代にNASAランレーが炭素繊維/エポキシで実測し、厚さ方向にはHashin-Shtrikmanの上下限モデルを推奨した報告（NASA TM-4756）がある。
線形要素 vs 2次要素
熱伝導解析では線形要素でも十分な精度が得られることが多い。温度勾配が急な領域（熱衝撃等）では2次要素を推奨。
熱流束の評価
要素内の温度勾配から算出。節点応力と同様にスムージングが必要な場合がある。
対流-拡散問題
ペクレ数が高い（対流支配）場合、風上的安定化（SUPG等）が必要。純粋な熱伝導問題では不要。
非定常解析の時間刻み
熱拡散の特性時間 $\tau = L^2 / \alpha$（$\alpha$: 熱拡散率）に対して十分小さい刻みを設定。急激な温度変化には自動時間刻み制御が有効。
非線形収束
温度依存物性値による非線形性はマイルドな場合が多く、Picard反復（直接置換法）で十分なことが多い。放射の強非線形性ではニュートン法を推奨。
定常解析の判定
全節点の温度変化が閾値以下（$|\Delta T| / T_{max} < 10^{-5}$等）で収束と判定。
陽解法と陰解法のたとえ
陽解法は「今の情報だけで次を予測する天気予報」——計算は速いが大きな時間刻みでは不安定（嵐を見逃す）。陰解法は「未来の状態も考慮した予測」——大きな時間刻みでも安定するが、各ステップで方程式を解く手間がかかる。急激な温度変化がない問題では陰解法で大きな時間刻みを使う方が効率的。
実践ガイド
実務での注意点
複合材料の熱解析で特に気をつけることは何ですか？
最も重要なのは材料座標系の正しい定義だ。プライごとに繊維方向が異なるので、各要素の座標系を正確に設定する必要がある。
Ansys ACPとの連携
Ansys Composite PrepPost (ACP) を使えば、複合材積層の定義から熱解析モデルへの転写が自動化できる。
1. ACPで積層構成（プライ順序、配向角、厚さ）を定義
2. 材料座標系が自動生成される
3. Steady-State Thermalにデータ転送
4. 各プライの異方性kが自動適用
手動で座標系を設定しなくていいのは大きいですね。
AbaqusならAbaqus/CAEのComposite Layup機能で同等のことができる。COMSOLではComposite Materials Moduleが対応する。
実測値との比較
複合材料の熱伝導率は測定方法によってばらつく。
| 測定法 | 適用 | 精度 |
|---|---|---|
| レーザーフラッシュ法 (LFA) | 面直方向 | ±5% |
| 定常法（保護熱板法） | 面直方向 | ±3% |
| Angstrom法 | 面内方向 | ±10% |
面内と面直で測定方法が違うんですね。
レーザーフラッシュ法は面直方向の熱拡散率を測定し、比熱と密度から $k = \alpha \rho c_p$ で換算する。面内方向の測定はサンプル加工が難しく精度も劣る。解析と実測の比較ではこの不確かさを考慮する必要がある。
スマホ放熱材の進化
2010年代初頭のスマートフォンは銅箔グラファイトシート（λ≈400 W/m·K）を貼るだけだったが、2019年以降のハイエンド機（Samsung Galaxy S10以降）では気相堆積グラファイト（VGCFs）複合材が採用され、面内熱伝導率は1500 W/m·Kを超えた。
解析フローのたとえ
熱解析のフローは「お風呂の追い焚き設計」で考えてみましょう。浴槽の形（解析対象）を決め、お湯の初期温度（初期条件）と外気温（境界条件）を設定し、追い焚きの出力（熱源）を調整する。「2時間後にぬるくなっていないか？」を計算で予測する——これが非定常熱解析の本質です。
初心者が陥りやすい落とし穴
「放射を無視していいですか？」——室温付近なら大抵OK。でも数百度を超えたら話は別です。放射による熱伝達は温度の4乗に比例するため、高温では対流を圧倒します。晴れた日に日向と日陰で体感温度が全然違うのを経験したことがありますよね？ あれが放射の威力です。工業炉やエンジン周りの解析で放射を無視するのは、猛暑日に「日差しは関係ない」と言い張るようなものです。
境界条件の考え方
熱伝達係数 $h$ は「窓の断熱性能」だと思ってください。$h$ が大きい＝窓が薄い＝熱がどんどん逃げる。$h$ が小さい＝二重窓＝熱が逃げにくい。この数値1つで結果が大きく変わるため、文献値の引用や実験による同定が重要です。「とりあえず10 W/(m²·K)で…」と適当に入れていませんか？
ソフトウェア比較
専用ツールの比較
複合材料の熱伝導に特化したツールはありますか？
均質化（ミクロ）とマクロ解析の両方に対応するツールを比較する。
| ツール | 均質化 | マクロ解析 | 特徴 |
|---|---|---|---|
| Digimat (e-Xstream/MSC) | RVE自動生成 | FEA連携 | 射出成形材の繊維配向連携 |
| COMSOL Composite | 内蔵RVE | FEM | マルチフィジックス対応 |
| Ansys ACP + Mechanical | Ansys Materialsで均質化 | FEM | Workbenchワークフロー |
| Abaqus + Micromechanics Plugin | RVE解析 | FEM | Python自動化が強力 |
Digimatが射出成形との連携に強いんですね。
DigimatはMoldflowやMoldex3Dから繊維配向テンソルを読み込み、各要素の異方性kを自動マッピングできる。短繊維強化PA66やPPS等の樹脂成形品の熱解析には事実上の標準ツールだ。
データベースの活用
複合材料の熱物性値はばらつきが大きいので、信頼性の高いデータベースを使うことが重要だ。
- AGATE/NCAMP: 航空宇宙用プリプレグの認定データ
- Granta (Ansys): 材料データベース、複合材料プロパティ含む
- CAMPUS: 熱可塑性樹脂の物性DB
自分で測定しなくてもデータベースから引ける場合があるんですね。
航空宇宙分野では認定プリプレグのデータが充実しているが、新規材料は必ず実測が必要だ。データベースの値はあくまで参考値として、設計の初期段階で使う。
Ansys Composite PrepPostの登場
Ansys（旧ANSYS）は2010年にACP（Ansys Composite PrepPost）をリリースし、積層角度ごとの熱伝導テンソルを自動変換する機能を実装した。それ以前はユーザーがPythonスクリプトで座標変換行列を手書きするのが一般的で、入力ミスによる解析事故が後を絶たなかった。
選定で最も重要な3つの問い
- 「何を解くか」：複合材料の熱伝導に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
CNT/グラフェン強化複合材料
カーボンナノチューブやグラフェンで熱伝導率を上げる研究って進んでいますか？
個々のCNTの熱伝導率は3000 W/(mK)以上だが、複合材料中での界面熱抵抗（Kapitza抵抗）により、実効的な改善は限定的だ。現状では面直kが0.5→1.5 W/(mK)程度に改善される。
$G_{interface}$ はCNT-マトリクス界面のコンダクタンスで、10〜100 MW/(m2K)程度だ。
多機能複合材料
構造と熱機能を兼ねる多機能複合材料の設計にマルチスケール熱解析が活用されている。例えばEVバッテリーケースのCFRP化では、構造強度と放熱性を同時に最適化する。
構造と熱の両方を満たす設計って難しそうですね。
多目的最適化（パレートフロント）で構造重量と最大温度のトレードオフを可視化し、積層構成と厚さを決定する。Ansys optiSLangやmodeFRONTIERが使われる。
3Dプリント複合材料
連続繊維3Dプリント（Markforged等）では、印刷パスが繊維方向を決定する。印刷パスデータから直接熱伝導率テンソルをマッピングするワークフローが開発されている。
製造プロセスと解析が直結するんですね。
Design-for-AM（積層造形設計）と熱解析の統合は今後の大きなトレンドだ。印刷方向による異方性を積極的に活用して、放熱経路を最適化する設計が可能になる。
CNT複合材の熱伝導限界
カーボンナノチューブ（CNT）単体の熱伝導率は3000 W/m·K超だが、CNT/エポキシ複合材では界面熱抵抗（カピタンス抵抗）が支配的になり、実測値は5〜10 W/m·Kにとどまる。2004年にPhilip Kimらがナノスケール熱輸送の界面抵抗を定量化した論文（Nano Lett.）が転換点となった。
トラブルシューティング
よくある問題
複合材料の熱解析で特に失敗しやすい点は何ですか？
1. 材料座標系の不整合
問題: 繊維方向をグローバルX方向と仮定してkを入力したが、実際の部品形状では繊維がドレープに沿って湾曲している。
対策: 曲面部品ではAnsys ACPやAbaqus Composite Layupで面に沿った座標系を自動生成する。平面投影で座標系を定義すると角部でずれる。
2. 等方性仮定の誤り
「異方性は面倒だから等方性の平均値で入れよう」はダメですか？
完全にダメだ。面内k=3 W/(mK)、面直k=0.45 W/(mK)を平均して1.7で入力すると、面直方向の温度降下を3.8倍も過小評価する。面直の熱抵抗が支配的な電子基板のような問題では致命的な誤差になる。
3. プライドロップの熱影響
厚さ遷移部（プライドロップ）では断面が変化し、局所的な温度勾配集中が生じる。メッシュがプライドロップ形状に追従していないと温度場が不正確になる。
対策: プライドロップ領域のメッシュを手動で細分化する。ACPのドレープ機能でプライドロップ位置を正確に反映させる。
4. 繊維体積率のばらつき
実部品では $V_f$ が設計値から±5%程度ばらつく。$V_f$ が5%変わると $k_{\parallel}$ は約5%、$k_{\perp}$ は約15%変化する。感度が高い面直方向に注意して、ばらつきを考慮した設計マージンを設ける。
面直方向がいつもボトルネックなんですね。
そうだ。複合材料の熱設計では「面直方向の熱抵抗をいかに下げるか」が永遠のテーマだ。
異方性を無視した設計ミス
2008年、某衛星メーカーがCFRP構体の熱解析を等方性（λ=5 W/m·K）で実施し、実際の厚さ方向λ≈0.6 W/m·Kを見落としてパネル温度を30℃低く見積もった。打ち上げ後に温度超過が発覚し、軌道上でヒータ出力を再プログラムする事態になった実例が報告されている。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——複合材料の熱伝導の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告