複合材料の熱伝導
複合材料の熱伝導の理論基礎
複合材料の熱伝導特性
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))という形で教科書に残る。
複合材料の熱伝導の数値計算手法
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）がある。
複合材料の熱伝導の実務適用
実務での注意点
複合材料の熱解析で特に気をつけることは何ですか？
最も重要なのは材料座標系の正しい定義だ。プライごとに繊維方向が異なるので、各要素の座標系を正確に設定する必要がある。
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を超えた。
複合材料の熱伝導のソフトウェア比較
専用ツールの比較
複合材料の熱伝導に特化したツールはありますか？
均質化（ミクロ）とマクロ解析の両方に対応するツールを比較する。
| ツール | 均質化 | マクロ解析 | 特徴 |
|---|---|---|---|
| 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スクリプトで座標変換行列を手書きするのが一般的で、入力ミスによる解析事故が後を絶たなかった。
複合材料の熱伝導の先端研究
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℃低く見積もった。打ち上げ後に温度超過が発覚し、軌道上でヒータ出力を再プログラムする事態になった実例が報告されている。
なった
詳しく
報告