古典積層理論（CLT）
理論と物理
CLTとは
先生、古典積層理論（CLT: Classical Lamination Theory）って何ですか？
CLTは繊維強化複合材（FRP）の積層板の力学を記述する基礎理論だ。キルヒホッフ板理論を異方性積層板に拡張したもので、各層の材料特性と繊維角から積層板全体の剛性を計算する。
等方性の板理論とどう違うんですか？
等方性板では曲げ剛性 $D = Et^3/12(1-\nu^2)$ の1つのパラメータで全て決まる。CLTではABD行列と呼ばれる6×6のマトリクスが必要。膜剛性、曲げ剛性、膜-曲げ連成を全て含む。
ABD行列
CLTの中心概念であるABD行列：
ここで：
- $\{N\}$ — 膜力合力（$N_x, N_y, N_{xy}$）
- $\{M\}$ — 曲げモーメント合力（$M_x, M_y, M_{xy}$）
- $\{\varepsilon^0\}$ — 中立面のひずみ
- $\{\kappa\}$ — 曲率
$[A]$, $[B]$, $[D]$ はそれぞれ何を表しますか？
$[B]$ が面白いですね。等方性板では $[B] = 0$ ですよね？
その通り。等方性板や対称積層（例: $[0/90]_s$）では $[B] = 0$。非対称積層では $[B] \neq 0$ で、引張をかけると板が反るという不思議な挙動が起きる。
ABD行列の計算
各マトリクスの計算式：
ここで $\bar{Q}_{ij}^{(k)}$ は第 $k$ 層の変換剛性マトリクス（繊維角を考慮）、$z_k$ は層の位置。
各層の位置と繊維角から、積層板全体のABD行列が計算できるんですね。
そう。CLTは「各層の材料特性＋繊維角＋積層順序」から「積層板の巨視的剛性」を導出する理論だ。FEMのシェル要素で複合材を扱う場合、内部ではCLTが使われている。
積層表記
「$[0/90/\pm 45]_s$」のような表記はどう読むんですか？
繊維角を下から順に記述する：
- $[0/90/\pm 45]_s$ = $[0/90/+45/-45/-45/+45/90/0]$
- $s$ は対称（symmetric）を意味。中立面で折り返す
- 対称積層なら $[B] = 0$（膜-曲げ連成なし）
対称積層にするのが設計の基本ですか？
基本的にはそう。非対称積層は成形後に反り（curing warpage）が発生するため、製造上も問題がある。航空機の構造は原則として対称積層だ。
まとめ
CLTの理論を整理します。
要点：
- ABD行列が積層板の剛性を記述 — $[A]$: 膜、$[B]$: 連成、$[D]$: 曲げ
- 各層の $\bar{Q}$ と位置 $z_k$ からABD行列を計算 — 層の材料特性と繊維角が入力
- $[B] = 0$ なら膜-曲げ連成なし — 対称積層で達成
- 対称積層が設計の基本 — 反りを防止
- FEMの複合材シェル要素は内部でCLTを使用 — CLTはFEMの基盤理論
CLTを理解していないと、FEMの複合材解析の設定も結果の解釈もできないんですね。
まさにそう。CLTは複合材設計の「アルファベット」だ。これなしに複合材のFEM解析を語ることはできない。
CLTの発祥とReissner・Mindlinの貢献
古典積層理論（CLT）は1960〜70年代にWhitney・Leissaらが確立した。CLTは各積層を平面応力状態と仮定し、Kirchhoff板理論（面外変形は法線を剛直に保つ）をベースに積分して全体剛性行列[ABD]を求める。Reissner（1945年）のせん断変形理論とMindlinの改良（1951年）はCLTの「せん断剛性なし」の限界を補う一次せん断変形理論（FSDT）に発展した。
各項の物理的意味
- 慣性項（質量項）：$\rho \ddot{u}$、つまり「質量×加速度」。急ブレーキで体が前に投げ出された経験はありませんか？ あの「持っていかれる感じ」がまさに慣性力です。重い物体ほど動き出しにくく、動き出したら止まりにくい。地震で建物が揺れるのも、地面が急に動いたのに建物の質量が「置いていかれる」から。静解析ではこの項をゼロにしますが、それは「ゆっくり力をかけるから加速度は無視できる」という仮定です。衝撃荷重や振動問題では絶対に省略できません。
- 剛性項（弾性復元力）：$Ku$ や $\nabla \cdot \sigma$。ばねを引っ張ると「戻ろうとする力」を感じますよね？ あれがフックの法則 $F=kx$ であり、剛性項の本質です。では質問——鉄の棒とゴム紐、同じ力で引っ張るとどちらが伸びるでしょうか？ 当然ゴムです。この「伸びにくさ」がヤング率 $E$ であり、剛性を決めます。よくある勘違い：「剛性が高い＝強い」ではありません。剛性は「変形しにくさ」、強度は「壊れにくさ」で、別の概念です。
- 外力項（荷重項）：体積力 $f_b$（重力など）と表面力 $f_s$（圧力、接触力など）。こう考えてみてください——橋の上のトラックの重さは「中身全体にかかる力」（体積力）、タイヤが路面を押す力は「表面だけにかかる力」（表面力）。風圧、水圧、ボルトの締付力…すべて外力です。ここでありがちな失敗：荷重の方向を間違える。「引張」のつもりが「圧縮」になっていた——笑い話に聞こえますが、3D空間で座標系が回転していると実際に起こります。
- 減衰項：レイリー減衰 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。ギターの弦を弾いてみてください。音は鳴り続けますか？ いいえ、徐々に小さくなりますよね。振動エネルギーが空気抵抗や弦の内部摩擦で熱に変わるからです。車のショックアブソーバーも同じ原理——わざと振動エネルギーを吸収して乗り心地を良くしています。もし減衰がゼロだったら？ 建物は地震の後いつまでも揺れ続けることになります。実際にはそうならないので、適切な減衰の設定が重要です。
仮定条件と適用限界
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 変位 $u$ | m（メートル） | mm入力時は荷重・弾性率もMPa/N系に統一すること |
| 応力 $\sigma$ | Pa（パスカル）= N/m² | MPa = 10⁶ Pa。降伏応力との比較時に単位系の不一致に注意 |
| 歪み $\varepsilon$ | 無次元（m/m） | 工学歪みと対数歪みの区別に注意（大変形時） |
| 弾性率 $E$ | Pa | 鋼: 約210 GPa、アルミ: 約70 GPa。温度依存性に注意 |
| 密度 $\rho$ | kg/m³ | mm系ではtonne/mm³（= 10⁻⁹ tonne/mm³ for 鋼） |
| 力 $F$ | N（ニュートン） | mm系ではN、m系ではNで統一 |
数値解法と実装
FEMでのCLT実装
FEMのシェル要素でCLTはどう実装されていますか？
複合材のシェル要素では、各層の $\bar{Q}$ マトリクスと積層情報（繊維角、板厚、位置）からABD行列を構成し、要素の剛性マトリクスに反映する。板厚方向の積分点は各層に配置される。
Nastran
```
PCOMP, 1, , , , SYM
, 1, 0.125, 0., YES,
, 1, 0.125, 90., YES,
, 1, 0.125, 45., YES,
, 1, 0.125, -45., YES
```
PCOMPカードで各層の材料ID、板厚、繊維角を指定。SYMで対称積層。
Abaqus
```
*SHELL SECTION, COMPOSITE, ELSET=panel
0.125, 3, CFRP, 0.
0.125, 3, CFRP, 90.
0.125, 3, CFRP, 45.
0.125, 3, CFRP, -45.
```
各行が1層を定義。板厚、積分点数、材料名、繊維角。
Ansys
Workbenchでは「ACP（Ansys Composite PrePost）」で積層を定義。GUI操作で直感的にドレーピングや繊維角を設定できる。
NastranのPCOMPが航空宇宙で広く使われている理由は？
PCOMPは1980年代からの歴史があり、航空宇宙の認証（型式証明）で膨大な検証実績がある。各層のひずみ・応力を直接出力でき、破壊判定（Tsai-Wu, Hashin等）との連携も充実している。
板厚方向の積分点
各層に何個の積分点が必要ですか？
各層の板厚方向に最低3点（Simpson積分）が推奨。弾塑性を扱う場合は5点以上。
全層で $n$ 層 × 3点 = $3n$ 点。20層の積層板なら60点。板厚方向の積分点が多いと計算コストが増えるが、精度には重要。
材料座標系
繊維角はどの座標系に対して定義しますか？
材料座標系は各層ごとに定義される。通常はシェル要素の面内方向（1軸方向）に対する角度で繊維角を指定。要素の向きが変わると繊維角の基準も変わるため、材料方向の確認が不可欠。
ドレーピング（曲面上での繊維角の変化）を正確に表現するには、各要素で材料座標系を個別に設定する必要がある。ACP（Ansys）やFibersim（Siemens）がこの作業を自動化する。
まとめ
CLTの数値手法、整理します。
要点：
- PCOMP（Nastran）, *SHELL SECTION COMPOSITE（Abaqus）, ACP（Ansys） — 積層定義
- 各層に板厚方向3点以上の積分点 — 精度のため
- 材料座標系の確認が不可欠 — 繊維角の基準方向
- ドレーピングの自動化 — ACP, FibersimなどのCAEツール
- NastranのPCOMPが航空宇宙で圧倒的実績 — 認証での検証実績
[ABD]行列の計算手順
CLTの[ABD]行列はAij（面内剛性）・Bij（カップリング）・Dij（曲げ剛性）の3部分からなる。Aij=Σ(Q̄ij)k×hkで、Q̄は座標変換された直交異方性剛性行列で計算する。Bij≠0は「非対称積層」でのカップリング（引張で面外変形）を表し、設計上は対称積層（Bij=0）にすることが多い。炭素繊維強化プラスチック（CFRP）[0/90/±45]s積層ならB=0となりsymmetric balanced laminateとなる。
線形要素（1次要素）
節点間を線形補間。計算コストは低いが、応力の精度が低い。せん断ロッキングに注意（低減積分やB-bar法で緩和）。
2次要素（中間節点付き）
曲線的な変形を表現可能。応力精度が大幅に向上するが、自由度は約2〜3倍に増加。推奨：応力評価が重要な場合。
完全積分 vs 低減積分
完全積分：過剰拘束（ロッキング）のリスク。低減積分：アワーグラスモード（零エネルギーモード）のリスク。適材適所で選択。
アダプティブメッシュ
誤差指標（ZZ推定量等）に基づく自動細分化。応力集中部の精度を効率的に向上。h法（要素分割）とp法（次数増加）がある。
ニュートン・ラフソン法
非線形解析の標準的手法。接線剛性マトリクスを毎反復更新。収束半径内で2次収束するが、計算コストが高い。
修正ニュートン・ラフソン法
接線剛性マトリクスを初期値または数反復毎に更新。各反復のコストは低いが、収束速度は線形的。
収束判定基準
力の残差ノルム: $||R|| / ||F_{ext}|| < \epsilon$（一般に $\epsilon = 10^{-3}$〜$10^{-6}$）。変位増分ノルム: $||\Delta u|| / ||u|| < \epsilon$。エネルギーノルム: $\Delta u \cdot R < \epsilon$
荷重増分法
全荷重を一度に負荷せず、小刻みに増加させる。弧長法（Riks法）は荷重-変位関係の極値点を越えて追跡可能。
直接法 vs 反復法のたとえ
直接法は「連立方程式を筆算で正確に解く」方法——確実だが大規模問題では時間がかかりすぎる。反復法は「当て推量を繰り返して正解に近づく」方法——最初は大雑把な答えだが、反復するたびに精度が上がる。辞書で言葉を探すとき、最初のページから順番に探す（直接法）より、見当をつけて開き、前後に調整する（反復法）方が効率的なのと同じ原理。
メッシュの次数と精度の関係
1次要素は「定規で曲線を近似する」——直線の折れ線で表現するため精度に限界がある。2次要素は「フレキシブルカーブ」——曲線的な変化を表現でき、同じメッシュ密度でも格段に精度が向上する。ただし、1要素あたりの計算コストは増えるため、トータルのコスト対効果で判断する。
実践ガイド
CLTの実務適用
CLTは実務でどう使われていますか？
複合材構造の設計は全てCLTが基盤だ。航空機、自動車、風車、スポーツ用品…繊維強化複合材を使う全ての分野で必須の理論。
積層設計のルール
積層設計にはルールがありますか？
航空宇宙の設計ルール（典型的）：
| ルール | 理由 |
|---|---|
| 対称積層 | 反りの防止（$[B] = 0$） |
| バランス積層 | $+\theta$ と $-\theta$ を同数で。面内せん断-膜連成の防止 |
| 10%ルール | 各方向（0°, 90°, ±45°）に最低10%の層を確保 |
| 隣接層の角度差 | 隣接する層の繊維角の差を45°以下に（層間剥離防止） |
| 外層を±45°に | 耐損傷性の向上 |
10%ルールって何ですか？
どの方向にも最低10%の繊維を入れるルール。例えば100層なら0°, 90°, +45°, -45°がそれぞれ最低10層。特定方向にゼロの層があると、その方向の強度が極端に低くなり、予想外の破壊が起きる。
等価材料特性
積層板の等価ヤング率はどう計算しますか？
$[A]$ マトリクスから等価の面内特性が計算できる：
ここで $h$ は全板厚。
等方性材料のように1つの $E$ ではなく、$x$ 方向と $y$ 方向で異なるんですね。
$[0/90]_s$ の等価積層板は $E_x \neq E_y$（直交異方性）。$[0/±60]_s$ のような擬似等方性積層なら $E_x \approx E_y$ で等方性に近づく。
実務チェックリスト
CLTのチェックリストをお願いします。
積層設計のルールが多いですね。金属構造にはないチェック項目ばかり。
複合材は自由度が高い分、設計ルールも多い。積層順序1つで剛性も強度も大きく変わる。この自由度が複合材の魅力であり難しさでもある。
CFRP翼構造の積層設計最適化
Airbus A350のCFRP主翼は[45/0/-45/90]4sのような[A]平均積層を基本に、荷重方向に合わせた最適積層比率（0°: 57%、±45%: 33%、90%: 10%など）が採用されている。CLTによる解析で翼下面引張側は0°繊維を増やし、端部せん断が大きい桁ウェブは±45°を多くする設計最適化を行い、アルミニウム翼比較で30%の重量削減を達成した。
解析フローのたとえ
解析の流れは、実は料理とそっくりです。まず材料を買い出し（CADモデルの準備）、下ごしらえをして（メッシュ生成）、火にかけて（ソルバー実行）、最後に盛り付ける（後処理で可視化）。ここで大事な問いかけ——料理で一番失敗しやすい工程はどこでしょう？ 実は「下ごしらえ」なんです。メッシュの品質が悪いと、どんなに優秀なソルバーを使っても結果はめちゃくちゃになります。
初心者が陥りやすい落とし穴
あなたはメッシュ収束性を確認していますか？ 「計算が回った＝結果が正しい」と思っていませんか？ これ、実はCAE初心者が最も陥りやすい罠です。ソルバーは与えられたメッシュで「それなりの答え」を必ず返します。でもメッシュが粗すぎれば、その答えは現実から大きくずれている。最低3段階のメッシュ密度で結果が安定することを確認する——これを怠ると「コンピュータが出した答えだから正しいはず」という危険な思い込みに陥ります。
境界条件の考え方
境界条件の設定は、試験の「問題文を書く」のと同じです。問題文が間違っていたら？ どんなに正確に計算しても答えは間違いますよね。「この面は本当に完全固定なのか」「この荷重は本当に一様分布なのか」——現実の拘束条件を正しくモデル化することが、実は解析全体で最も重要なステップだったりします。
ソフトウェア比較
CLTのツール
CLTの計算や複合材解析に使えるツールを教えてください。
CLT計算ツール
| ツール | 特徴 |
|---|---|
| eLamX² | ドレスデン工科大学開発。無料。ABD行列、ひずみ、破壊判定 |
| CLT Calculator（各種） | 無料Webツール。ABD行列の計算 |
| MATLAB/Python | CLTの式を自作。研究用 |
| ESAComp | ESA公認。複合材の設計・解析・最適化 |
eLamX²が無料で使えるのはいいですね。
eLamX²はJavaベースで、ABD行列の計算、各層のひずみ・応力、破壊判定（Tsai-Wu, Hashin, Puck等）まで全て無料で計算できる。FEMを使う前の「積層設計のスクリーニング」に最適だ。
FEMソルバーの比較
| 機能 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 積層定義 | PCOMP/PCOMPG | *SHELL SECTION COMPOSITE | ACP |
| 破壊判定 | 内蔵（Tsai-Wu等） | USDFLD/VUMAT | ACP Post |
| 層間剥離 | CVISC/CZM | CZM（Cohesive） | CZM |
| ドレーピング | 外部ツール連携 | Abaqus/CAE | ACP（統合） |
| 積層最適化 | SOL 200 | Pythonスクリプト | Workbench + OptiSlang |
AnsysのACPが使いやすそうですね。
ACPは複合材の前後処理に特化したモジュールで、GUI操作で積層定義、ドレーピング、破壊判定まで一貫して行える。特にWorkbenchとの統合がスムーズだ。
NastranのPCOMPは航空宇宙で圧倒的実績。Abaqusは非線形（層間剥離、プログレッシブ損傷）に強い。Ansys ACPは使いやすさで優位。
選定ガイド
複合材はツールの選択肢が多いですね。
複合材は異方性・積層・破壊判定・ドレーピングと独自の要素が多いから、専用ツールの価値が大きい。CLTの基礎を理解した上で適切なツールを選ぶこと。
Ansys CompositePrepPost（ACP）の活用
Ansys ACP（旧ESAComp）はCLT計算と連続FEMの橋渡しをするCFRP設計の専門ツールだ。積層角度・厚さ・材料の3次元的なばらつきを入力にFEMのシェル要素に積層情報を自動割り当てし、CLT計算結果のプロット・破損指標の可視化を一括実行する。Saab空軍機のCFRP部品設計に採用され、積層エラーの検出率を50%向上させた。
選定で最も重要な3つの問い
- 「何を解くか」：古典積層理論（CLT）に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
CLTの先端研究
CLTに関する最新の研究はありますか？
CLTは1960年代に確立された理論だが、新材料と新製造法で進化がある。
可変剛性複合材（Variable Stiffness）
AFP（Automated Fiber Placement）で繊維の配置角を場所ごとに変える可変剛性複合材。従来のCLTは「各層の繊維角が一様」を仮定するが、可変剛性ではCLTを局所的に適用する。各要素で異なるABD行列を持つ。
座屈荷重を30〜50%向上できるとの研究がありましたよね。
そう。応力の流れに沿って繊維を配置することで、材料の利用効率を飛躍的に高められる。ただし設計変数が膨大になるため、最適化手法（遺伝的アルゴリズム、勾配法）の発展が不可欠だ。
高次積層理論
CLTはキルヒホッフの仮定（せん断変形無視）に基づく。厚い複合材板では高次せん断変形理論（HSDT: Higher-Order Shear Deformation Theory）が必要。Reddy（1984）の第3次理論やCarreraのCUFが代表的。
複合材は層間のせん断が重要ですよね。
金属板では板厚方向のせん断は副次的だが、複合材では繊維方向と直交方向のせん断剛性が桁違いに異なるため、層間せん断が破壊（層間剥離）の主因になる。HSDTやCUFでこの層間せん断を正確に評価する。
デジタルツインと複合材
航空機の複合材構造のデジタルツインでは、製造時の繊維角のばらつき、ボイド率、板厚変動をFEMモデルに反映する。CLTのABD行列を「実測の積層情報」で更新し、個体ごとの構造特性を予測する。
まとめ
CLTの先端研究、まとめます。
CLTは「枯れた理論」ではなく、複合材構造の全ての解析・設計の基盤であり続けている。
非線形CLT：幾何学的非線形と材料非線形の組合せ
厚みの大きな積層板（L/t<10）ではCLTは精度不足でHigher-Order Shear Deformation Theory（HSDT）が必要だ。さらにCFRP積層板の初期破損後の「ノックダウン」を評価するためにCLT＋Progressive Damage解析の連成が必要で、Linde・Ladeveze・Hashinの各破損基準が使われる。NASA製複合材シェルの設計ではこの連成解析が2005年から型式証明の要件となっている。
トラブルシューティング
CLTのトラブル
複合材解析でよくあるトラブルを教えてください。
複合材特有のトラブルは「材料座標系」と「積層順序」に集中している。
材料座標系の間違い
繊維角を正しく設定したはずなのに、結果がおかしいです。
材料座標系の向きが間違っている可能性が高い。繊維角は要素の材料座標系の1軸方向に対する角度。要素の向きが意図と異なっていると、繊維角の基準がずれる。
確認方法：
- プリプロセッサで材料方向ベクトルを表示 — 矢印がどこを向いているか目視確認
- 0°層のみのモデルで引張試験をシミュレーション — 繊維方向に引っ張って $E_1$ が出るか確認
- 曲面では要素ごとに材料方向が変わる — ドレーピングツールの使用を検討
曲面で材料方向がずれるのは厄介ですね。
曲面上では要素の面内方向（1軸）が連続的に変化するため、繊維角の基準も変わる。材料方向の可視化なしに複合材解析をやるのは危険だ。
積層順序の間違い
積層の順番を間違えるとどうなりますか？
$[A]$ マトリクスは積層順序に依存しない（面内剛性は層の合計で決まる）。しかし $[D]$ マトリクスは積層順序に強く依存する。外層に0°を配置するか内層に配置するかで曲げ剛性が大きく変わる。
同じ層構成でも順番で曲げ剛性が変わる！
$[D]$ は $z^3$ の重み付きだから、外層（$z$ が大きい）の層が曲げ剛性に大きく寄与する。外層に0°繊維を置くと$x$方向の曲げ剛性が最大になる。
結果の解釈ミス
結果の解釈で気をつけることは？
複合材ではグローバル座標系の応力を見てはいけないんですか？
破壊判定（Tsai-Wu, Hashin等）は材料座標系（1: 繊維方向、2: 直交方向、3: 板厚方向）で評価する。グローバル座標系のvon Mises応力は複合材には使えない。
まとめ
CLTのトラブル対処、整理します。
「材料座標系の可視化」が複合材解析の最も重要な確認項目ですね。
間違いない。繊維方向を正しくモデル化できなければ、全ての結果が無意味だ。複合材解析では材料方向の確認に時間の半分を使うくらいの心構えが必要だ。
CLTと実測剛性の15%以上の乖離
CLTの予測と実測のねじれ剛性・曲げ剛性が15%以上ずれる場合、積層の「巻き取り角度誤差」または「繊維体積分率のばらつき」が主因だ。CFRP製造では±2°の角度誤差が生じ、CLT計算で剛性が5〜10%変化する。デジタルカメラ断面観察で実際の積層角度を測定し、CLT入力に反映する「実製造品評価」が精度向上の近道だ。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——古典積層理論（CLT）の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告