FEMで複合材料を解析するとき、繊維の向きってどう入力するんですか？

シェル要素ベースが一般的だ。各プライに繊維方向角（例：0°、90°、±45°）と厚さを与えて積層定義する。AbaqusのComposite Shell SectionやAnsys ACPを使えばGUIで積層シーケンスを視覚的に定義できる。材料定数はE1（繊維方向）、E2（横方向）、G12、ν12の4定数が最低限必要で、これを各プライの局所座標系に紐づける。マルチプライ解析では各プライの応力とTsai-Wu破壊指数を確認するポスト処理が重要だ。

繊維一本一本をモデル化するんですか？それとも全体を均質な材料として扱うんですか？

通常は均質化された等価材料定数を使う。ミクロスケール（繊維径10μm程度）のRVE（代表体積要素）に周期境界条件を課して均質化計算をすると、等価弾性定数が導出できる。これをマクロスケール（部品レベル）のFEMに入力するのがマルチスケール法だ。局所的な繊維破断・マトリクスクラックが知りたい場合はRVEに戻って詳細解析する——いわゆる「ズームイン」アプローチだ。Digimat（MSC）やMicro2Macroがこの種のツールとして知られている。

複合材料の疲労ってどう評価するんですか？金属疲労とは違うんですか？

かなり違う。金属疲労は単一き裂の進展（Paris則）で評価するけど、複合材料疲労は複数損傷モードの複合体——マトリクスクラック、繊維切断、層間はく離が絡み合って進展する。解析的には残留剛性低下モデルや損傷力学アプローチが使われる。実務ではSN曲線（各積層方向・応力比R別）を試験で取得してPaldyなどの疲労損傷モデルに入力する手法が多い。試験コストが非常に高いため、少ない試験点から最大情報を引き出す試験計画法（DoE）との組み合わせが重要だ。

複合材料 — CAE用語解説

Q: 先生、複合材料ってCFRPだけじゃないんですよね？ CAE解析的にはどう分類するんですか？

複合材料は二種以上の材料を組み合わせて単体では得られない特性を実現したものの総称だ。CAEで扱う代表例は①繊維強化複合材（CFRP・GFRP・アラミド繊維）——異方性が強いので積層板理論（CLT）や3D直交異方性弾性が必要。②粒子強化複合材（セラミック粒子入り金属マトリクス）——均質化法で等価弾性定数を推定。③サンドイッチ構造（CFRP面材+アルミハニカムコア）——曲げ剛性と軽量を両立した構造だ。

カテゴリ: 用語集 | 2026-01-15

CAE visualization for composite material - technical simulation diagram

複合材料

🧑‍🎓

先生、複合材料ってCFRPだけじゃないんですよね？ CAE解析的にはどう分類するんですか？

理論と物理

複合材料の基本概念と支配方程式

🧑‍🎓

複合材料の「異方性」とは、具体的にどのような性質の違いを指すのですか？金属のような等方性材料と何が根本的に違うのでしょうか。

🎓

根本的な違いは、材料定数の数です。等方性材料はヤング率とポアソン比の2つで弾性を記述できますが、一般的な直交異方性を持つ複合材ラミナでは、独立な弾性定数は9つ必要です。例えば、炭素繊維強化プラスチック（CFRP）の一方向材では、繊維方向のヤング率

$$ E_1 $$

は約140GPa、垂直方向

$$ E_2 $$

は約10GPaと、10倍以上の差が出ます。この方向依存性が設計の自由度であり、解析の難しさの源です。

🧑‍🎓

9つの弾性定数とは具体的に何ですか？また、それらはどのような実験で決まるのでしょう。

🎓

直交異方性の構成則は、

$$ \{\sigma\} = [C]\{\epsilon\} $$

で表され、剛性マトリックス[C]の独立成分は、

$$ E_1, E_2, E_3, \nu_{12}, \nu_{13}, \nu_{23}, G_{12}, G_{13}, G_{23} $$

の9つです。これらはASTM D3039（引張試験）やASTM D3518（面内せん断試験）などの規格に基づく試験で決定します。実務では、メーカーカタログ（例えばTorayca T800Sのデータシート）にこれらの値が掲載されています。

🧑‍🎓

ラミナとラミネートの違いは何ですか？CAEではどちらを扱うことが多いのですか。

🎓

ラミナは単一の繊維配向を持つ層（プレプレグ1枚分）を指し、ラミネートはそれらを積層した全体を指します。CAEでは、まず各ラミナの性質を定義し、積層順序（例えば[0/90/±45]s）と各層の厚みを入力して、ラミネート全体の等価剛性を古典積層理論（CLT）で計算します。実機設計では、航空機主翼スキンなど数十層からなるラミネートを扱うのが一般的です。

🧑‍🎓

古典積層理論（CLT）では、どのようにしてラミネートの剛性を計算するのですか？単純な平均ではないですよね。

🎓

その通り、厚み方向の積分になります。各ラミナの剛性マトリックス

$$ [Q] $$

をその配向角

$$ \theta $$

で座標変換した

$$ [\bar{Q}] $$

を用い、ラミネートの面内剛性

$$ [A] $$

、カップリング剛性

$$ [B] $$

、曲げ剛性

$$ [D] $$

を計算します。式で書くと、

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

のようになります。ここで

$$ z_k $$

は層の位置です。非対称積層では

$$ [B] $$

がゼロでなく、引張と曲げが連成するのが特徴です。

数値解法と実装

FEMにおける複合材モデリング

🧑‍🎓

FEMで複合材をモデル化する時、ソリッド要素とシェル要素、どちらを使うべきかの判断基準は何ですか？

🎓

基本的な判断基準は、スパン（長さ）に対する厚みの比です。航空機や風力発電ブレードの様な薄肉構造では、厚みがスパンの1/10以下であることが多く、シェル要素（AbaqusのS4R, AnsysのShell181）が標準です。一方、接着層や厚み方向の応力集中を評価する必要がある場合は、連続体シェル要素やソリッド要素（各層を別要素でモデル化）を用いますが、要素数が爆発的に増えるため、解析対象を局所化するのがコツです。

🧑‍🎓

シェル要素で積層を定義する「積層シェル」機能について教えてください。入力するのはどのようなデータですか？

🎓

積層シェルでは、各層（ラミナ）の材料、厚み、配向角、積層順序を定義します。Ansys Composite PrepPost (ACP) や Abaqus/CAEのComposite Layup機能では、テーブル形式で入力します。具体的には、Layer 1: Material=「T800S/Epoxy」, Thickness=0.2mm, Orientation=0度、のように定義していきます。ソルバーは内部でCLTを用いて

$$ [A], [B], [D] $$

マトリックスを計算し、要素剛性に組み込みます。

🧑‍🎓

要素座標系と材料座標系の違いがよくわかりません。配向角0度とは、どの座標系に対する角度ですか？

🎓

これは非常に重要です。要素座標系は要素の形状（例えばシェル要素の面内方向）で決まるローカルな座標系です。材料座標系（または材料主方向）は、繊維が配向している方向を定義する座標系です。配向角（例えば0度）は、「要素座標系の第1軸から材料座標系の第1軸（繊維方向）への回転角」として定義されます。Ansysでは「SHELL181の要素X軸」、Abaqusでは「シェル平面の1方向」が基準になります。この設定を間違えると、剛性が全く違う方向に働いてしまいます。

🧑‍🎓

離散化の際、メッシュサイズは複合材のどのような特徴サイズを基準に決めるべきですか？

🎓

複合材では、層の厚み（通常0.1〜0.3mm）よりも、面内の特徴長が重要です。特に、応力集中が発生するボルト穴周りや急激な形状変化部では、メッシュを細かくする必要があります。一つの目安は、せん断遅れ長さ（shear lag length）です。これは、

$$ \lambda^{-1} = \sqrt{\frac{G_{13}}{E_1} \cdot \frac{t}{h}} $$

のような形で表され、荷重が有効に伝達される長さを表します。ここで

$$ t $$

は層厚、

$$ h $$

は全体厚さです。この長さの1/2〜1/3程度のサイズでメッシュを切ると、層間応力の評価にある程度の精度が得られます。

実践ガイド

解析ワークフローと検証

🧑‍🎓

複合材構造の線形静解析を実行する際の、具体的なワークフローのステップを教えてください。

🎓

実務的なワークフローは以下の通りです。

1. **幾何モデルの単純化**: リブやスティフナーのみを面として抽出し、中間面モデルを作成。 2. **材料プロパティ定義**: カタログ値（例：Hexcel IM7/8552の

$$ E_1=171 GPa, E_2=9.08 GPa, G_{12}=5.29 GPa, \nu_{12}=0.32 $$

）を入力。強度値（

$$ X_t, X_c, Y_t, Y_c, S $$

）も破壊判定用に準備。 3. **積層定義**: 設計図面の積層表（[45/0/-45/90]sなど）を、ソフトウェアの積層定義ツールに転記。 4. **メッシュ生成**: 曲率変化の大きい領域はメッシュを細かく（要素サイズ5mm以下）、平坦部は粗く（20mm程度）。 5. **境界条件・荷重**: 実際の支持条件を単純化しすぎない（例えば、ピン支持ではなく面拘束）。 6. **求解**: 線形ソルバーで実行。 7. **後処理**: 変位、層ごとの応力、破壊指標（後述）を評価。

🧑‍🎓

解析結果を信頼するための「検証」は、具体的に何をすればいいですか？単に変形形状を見るだけでは不十分ですよね。

🎓

その通りです。最低限実施すべき検証は3つあります。

1. **質量・体積のチェック**: ソフトウェアが計算したラミネートの総質量が、設計値（層数×層厚×密度×面積）と数%以内で一致するか。 2. **反力のチェック**: 全ての反力の合力が、印加した荷重とつり合っているか（通常、誤差1e-3以下）。 3. **単純ケースでの検証**: 例えば、一方向材の単純梁を片持ち支持し、先端に荷重をかけた時、先端たわみが梁理論

$$ \delta = PL^3/(3E_1I) $$

で計算した値と一致するか。ここで

$$ I $$

は断面二次モーメントです。この検証で材料定数と配向角の設定が正しいか確認できます。

🧑‍🎓

複合材の破壊はどうやって判定するのですか？金属のように単一の等価応力では評価できないと思います。

🎓

複合材には普遍的な破壊基準はなく、いくつかの破壊則を使い分けます。最も広く使われるのはTsai-Wu則です。これは、

$$ F_i \sigma_i + F_{ij} \sigma_i \sigma_j = 1 $$

というテンソル形式の式で、引張・圧縮強度が異なることも考慮できます。他には、モードごとの破壊を判定するPuck則やHashin則があります。実務では、Ansys ACPの「Failure Criteria」やAbaqusの「Hashin Damage」など、ソフトウェアに実装された機能を用い、各層ごとに破壊指標（0から1の値、1で破壊）を出力し、最も危険な層とその破壊モードを特定します。

🧑‍🎓

層間はがれ（デラミネーション）の解析は、特別な手法が必要ですか？

🎓

はい、通常の線形静解析では評価できず、破壊力学に基づく手法か、コヒーシブゾーンモデル（CZM）が必要です。実用的なのは、層間に非常に薄いコヒーシブ要素を挿入し、その界面の強度と破壊エネルギー（例：モードIの

$$ G_{Ic} = 0.2 N/mm $$

）を定義する方法です。AbaqusではCohesive Element、AnsysではINTER202要素などを使います。これは非線形解析となり、計算コストが高くなります。初期検討では、層間せん断応力（

$$ \tau_{13}, \tau_{23} $$

）が許容値（例：40MPa）を超えないかを簡易チェックします。

ソフトウェア比較

主要CAEソフトウェアの機能比較

🧑‍🎓

Ansys、Abaqus、MSC Nastranで複合材解析を行う場合、それぞれの強みや特徴的な機能は何ですか？

🎓

良い質問です。ソフトウェア選定はプロジェクトのフェーズによります。

- **Ansys (ACP + Mechanical)**: 前処理ツール「ACP」が強力です。CADから直接積層領域を定義し、繊維配向を可視化しながら積層設計できます。製造性評価（ドレープシミュレーション）も可能で、設計初期段階でよく使われます。 - **Abaqus/CAE**: 非線形解析（接触、大変形、破壊）との連成が得意です。VUMAT（ユーザー定義材料サブルーチン）を使って、独自の複合材破壊則を実装する研究開発現場で重宝されます。 - **MSC Nastran**: 航空宇宙業界の歴史的な標準です。SOL 200を用いた複合材積層の最適化（層厚や配向角を設計変数）が強み。大規模な線形解析（固有値、周波数応答）が高速で安定しています。

🧑‍🎓

「ドレープシミュレーション」とは具体的に何を計算するのですか？解析結果の使い道は？

🎓

プレプレグシートを曲面金型に貼り付ける（ドレープする）際に、生地が伸びたり、繊維配向が変わったり、シワが発生したりする現象をシミュレーションします。Ansys ACPやPAM-FORMなどの専用ソフトが使われます。出力は、設計した配向角からのずれ（スキュー角）や、繊維の伸び率です。これらは製造上の許容値（スキュー角は通常±5度以内）と比較され、設計が製造可能か、または金型形状を修正するかの判断材料になります。設計と製造の橋渡しをする重要な工程です。

🧑‍🎓

専用の複合材プリポストツール（ACPなど）を使わず、汎用プリポストで複合材解析を行う場合の限界は何ですか？

🎓

主に「効率性」と「表現力」に限界があります。汎用ツールでは、数十層に及ぶ積層定義を一つ一つの要素プロパティとして手入力するのは現実的ではありません。また、後処理で「特定の層の繊維方向応力」だけを抽出してコンター表示する、といった操作が煩雑です。さらに、積層順序を変えた複数の設計案を簡単に比較できません。専用ツールは、積層表を一括管理し、層ごとの結果を自動的にマッピングし、レポート生成までを効率化します。小規模な検討以外では、専用ツールの導入が必須と言えます。

🧑‍🎓

オープンソースのFEMソフト（CalculiX, Code_Aster）で複合材解析は可能ですか？商用ソフトとのギャップは？

🎓

コアとなる求解機能（異方性材料の線形静解析）は可能です。CalculiXやCode_Asterはシェル要素での積層定義をサポートしています。しかし、ギャップは大きく3点あります。

1. **GUIと前処理**: 商用ソフトのような直感的な積層定義インターフェースがなく、入力ファイル（.inp, .comm）を手書きまたはスクリプトで生成する必要がある。 2. **検証済みの破壊則**: Tsai-Wuなどの標準的な破壊則がプリセットで実装されていない場合が多く、ユーザーが自分で実装（ユーザーサブルーチン）する必要がある。 3. **技術サポートと保証**: 航空機などの安全係数に関わる分野では、解析結果の保証とサポートが求められるため、商用ソフトが選ばれます。研究や教育、概念検証レベルではオープンソースも有力な選択肢です。

トラブルシューティング

よくあるエラーと対策

🧑‍🎓

複合材の解析で「剛性マトリックスが非正定値です」というエラーが出ました。原因と対策を教えてください。

🎓

これは最もよくあるエラーの一つで、材料定数の入力ミスまたは物理的に不可能な値が原因です。具体的には：

1. **ヤング率やせん断弾性係数に負の値が入っている**：単純な入力ミス。 2. **ポアソン比の値が異方性の安定条件を満たしていない**：直交異方性材料では、

$$ \nu_{12}^2 < E_1/E_2, \quad \nu_{13}^2 < E_1/E_3, \quad \nu_{23}^2 < E_2/E_3 $$

などの条件を満たす必要があります。例えば、

$$ E_1=140GPa, E_2=10GPa $$

の場合、

$$ \nu