複合材料 積層板解析(CLT) 戻る
複合材料解析ツール

複合材料 積層板解析(CLT・ABDマトリクス)

古典積層理論(Classical Laminate Theory)によりABDマトリクス・等価弾性係数・繊維角依存剛性・ply別失敗指数をリアルタイム計算します。

材料・積層パラメータ
材料プリセット
積層構成
例: [0/90/±45]s → 対称 / [0/45/-45/90]T → 合計
積層構成の解析に失敗しました
Plyカウント: 8
適用荷重 Nx
N/mm
単位幅あたり膜力(引張方向)
適用荷重 Ny
N/mm
計算結果
Ex
GPa
Ey
GPa
Gxy
GPa
νxy
総厚さ h
mm
最大 FI(Ply#)
最大失敗指数
Ex(θ) — 剛性異方性
Ply別 失敗指数 FI
理論・主要公式

$$\mathbf{N} = \mathbf{A}\,\boldsymbol{\varepsilon}^0 + \mathbf{B}\,\boldsymbol{\kappa}, \quad \mathbf{M} = \mathbf{B}\,\boldsymbol{\varepsilon}^0 + \mathbf{D}\,\boldsymbol{\kappa}$$

CLT の基本方程式。\(\mathbf{A}\):面内剛性行列、\(\mathbf{B}\):曲げ-面内カップリング行列、\(\mathbf{D}\):曲げ剛性行列 [N·mm]。

$$A_{ij} = \sum_{k=1}^{N} \bar{Q}_{ij}^{(k)}\,(z_k - z_{k-1})$$

各プライの変換済み剛性 \(\bar{Q}_{ij}\) を厚み方向に積分して積層板剛性を算出。積層角度 \(\theta\) によって \(\bar{Q}\) が変化する。

$$\text{Tsai-Wu FI:} F_1\sigma_1 + F_2\sigma_2 + F_{11}\sigma_1^2 + F_{22}\sigma_2^2 + F_{66}\tau_{12}^2 + 2F_{12}\sigma_1\sigma_2 = 1$$

FI ≥ 1 で破壊。各係数は繊維方向・横方向の引張・圧縮強度と面内せん断強度から決まる。

複合材料積層板解析(CLT・ABDマトリクス)とは

🙋
複合材料の積層板って、なぜ「ABDマトリクス」なんて複雑なものが必要なんですか?一枚の板のヤング率みたいに、単純な値で表せないんですか?
🎓
大まかに言うと、積層板は「方向によって強さが違う層」が何枚も重なったサンドイッチだからだよ。一枚の等方性板と違って、引っ張ると曲がったり(連成)、曲げるとねじれたりする。ABDマトリクスは、その複雑な剛性の関係をすべてまとめて表す「設計のレシピ」みたいなものなんだ。このシミュレーターで、上の「積層構成」を[0/90]と[0/45]に変えてみて。ABDマトリクスの数値、特にBマトリクスがどう変わるか確認してみて。
🙋
確かに、積層を対称じゃない[0/90]にすると、Bマトリクスにゼロじゃない数字が出てきました!これが「引っ張ると曲がる」ということですか?でも、実際の設計ではそんな変な挙動は困りませんか?
🎓
その通り!Bマトリクスがゼロでないと、面内の力(Nx, Ny)がかかっただけで板が「わん」と反ってしまう。実務では、航空機の翼やF1のモノコックのように、軽量化が命の構造物では、この連成効果を積極的に利用して最適化することもあるけど、多くの一般的な構造物では「対称積層」にしてB=0に抑えることが多いね。シミュレーターの「適用荷重Nx」のスライダーを動かして、非対称積層板の「曲率」がどう変わるか確認してみよう。
🙋
「等価弾性係数Ex」って、ABDマトリクスから出てくる値ですよね?これはどう使うんですか?積層板全体を、普通の鉄板みたいに「概略として一つの材料」として扱えるということですか?
🎓
いいところに気が付いたね。その通りで、ExやEyは積層板全体を等方性板と「見なした」ときの代表値だ。例えば、自動車のドアパネルに複合材料を使うとき、全体の剛性を既存の鋼板と比較したい場合などに使う。でも注意点があって、この等価値は「面内」の引張り圧縮に対してだけ有効で、曲げ剛性は別物だ。シミュレーターで「材料プリセット」を炭素繊維とガラス繊維に切り替えてみて。層の構成が同じでも、等価Exが大きく変わるのがわかるよ。

よくある質問

ABDマトリクスから、積層板に任意の面内力や曲げモーメントを加えたときのひずみ・曲率を計算できます。また、等価弾性係数からは積層板を均質材料とみなした剛性が得られるため、FEAの材料定数として入力したり、引張・曲げ剛性の設計目標との比較に使えます。
連成剛性Bと曲げ剛性Dは、各層の面外方向の位置(中立面からの距離)に依存するため、積層順序が変わるとこれらの値が変化します。一方、面内剛性Aは層の順序に依存せず、各層の厚さと剛性の総和で決まります。
失敗指数が1を超えると、そのplyが破壊する可能性があります。対策として、①該当plyの繊維角度を変更する、②plyの厚みを増やす、③高剛性・高強度の材料に変更する、④積層構成を対称にして連成変形を抑える、などが有効です。
いいえ、非対称積層でも計算可能です。ただし、非対称積層では連成剛性Bがゼロにならず、面内力と曲げ変形が連成するため、実際の製造時に反りやねじれが生じやすくなります。設計上は対称積層が推奨されますが、解析自体は任意の積層構成に対応します。

実世界での応用

航空宇宙機体:旅客機の主翼や尾翼、ロケットのフェアリングには、軽量かつ高剛性が要求されるため炭素繊維複合材料の積層板が多用されます。CLTを用いて、荷重条件に応じた最適な繊維配向角(例: 0°, ±45°, 90°の組み合わせ)と積層順序を設計し、ABDマトリクスから部材の変形や座屈強度を予測します。

自動車部品:F1やスポーツカーのモノコック(車体骨格)やサスペンションアームに複合材料が適用されます。ここでは、衝突安全性と軽量化の両立が課題となり、CLT解析で部材ごとに異なる積層設計を行います。特に衝撃吸収領域では、特定の層を意図的に破壊させる設計も行われます。

スポーツ用具:テニスラケット、ゴルフシャフト、自転車フレームなどでは、剛性としなり(ダンピング)のバランスが性能を左右します。CLTシミュレーションで、打球感や振動特性に影響する曲げ剛性(Dマトリクス)やねじり剛性を最適化する積層設計が行われています。

風力発電ブレード:巨大なブレードは複合材料の積層構造で作られ、重力、遠心力、風圧による複雑な荷重に耐える必要があります。CLTはブレード全体の構造解析の基礎として用いられ、部位ごとに異なる積層構成(スキンとスパー)を設計し、軽量化と耐久性の両立を図ります。

よくある誤解と注意点

CLT解析を始めたばかりの頃に陥りがちな落とし穴をいくつか挙げておくよ。まず「等価弾性係数Exで曲げも評価できる」という誤解。Exはあくまで面内引張り圧縮の代表値だ。例えばExが同じでも、[0/90]s(対称)と[0/90](非対称)では曲げ剛性Dマトリクスが全く異なり、曲げた時のたわみは大きく変わる。曲げ性能を見たいなら、Dマトリクスや曲げ剛性$D_c = 12D_{11}/h^3$を直接確認しよう。

次に繊維角の入力順序と積層表記のズレ。ツールで[0/45/90]と入力した時、それが板の上面から下面への順序なのか、下面から上面への順序なのかを意識しておこう。古典積層理論では通常、下面(z=-h/2)から層を積み上げていくが、ソフトによっては逆の定義もある。順序を間違えると、特に非対称積層でBマトリクスの符号が逆転してしまうことがあるんだ。

最後に「ツールの出力が正しい」という盲信。入力パラメータ、特に単層の基本剛性$Q_{11}, Q_{12}, Q_{22}, Q_{66}$は材料データシートから正確に引いてこないと、全てが台無しになる。例えば、よくある間違いが横弾性係数$G_{12}$の値だ。経験則で$E_2$の何分の一、と適当に入れると、せん断変形の予測が大きく外れる。まずは単一繊維角(例えば[0]一様板)でツールの出力が単層の理論値と一致するか、必ず検算するクセをつけよう。