直交異方性 — CAE用語解説

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

理論と物理

直交異方性の物理的意味

🧑‍🎓

「直交異方性」という言葉をよく見かけますが、等方性と具体的に何が違うんですか?

🎓

根本的な違いは、材料の特性が方向によって異なるかどうかです。等方性材料は、どの方向から引っ張ってもヤング率やポアソン比が同じです。直交異方性材料は、互いに直交する3つの主方向(例えば、繊維方向、面内直交方向、厚さ方向)で特性が異なります。例えば、CFRP(炭素繊維強化プラスチック)では、繊維方向のヤング率は約150 GPaなのに対し、繊維に垂直な方向では10 GPa程度と、10倍以上の差が出ます。

🧑‍🎓

3つの方向で違うのは分かりました。では、その特性を数式で表すと、等方性と比べてパラメータが増えるということですか?

🎓

その通りです。等方性の線形弾性は、ヤング率Eとポアソン比νの2つで記述できます。一方、直交異方性の完全な線形弾性構成則には、独立な材料定数が9つ必要です。3方向のヤング率

$$ E_1, E_2, E_3 $$
、せん断弾性係数
$$ G_{23}, G_{31}, G_{12} $$
、そしてポアソン比
$$ \nu_{12}, \nu_{23}, \nu_{31} $$
です。ここで
$$ \nu_{ij} $$
は、j方向の応力によるi方向のひずみを表します。

🧑‍🎓

9つもあると、実験で全部求めるのは大変じゃないですか?実際の設計ではどうしているんですか?

🎓

良い指摘です。実際、厚さ方向(3方向)の特性

$$ E_3, G_{23}, G_{31} $$
は測定が非常に困難な場合が多い。そのため、薄板構造(シェル要素)の解析では、面内特性(1,2方向)と面内せん断
$$ G_{12} $$
に焦点を当て、厚さ方向のせん断は「横せん断補正係数」として簡略化して扱うことがあります。航空宇宙分野では、複合材料の試験方法を定めたASTM D3039(引張試験)やASTM D3518(面内せん断試験)などの規格に従ってデータ取得します。

数値解法と実装

CAEソフトウェアでの入力と座標系

🧑‍🎓

CAEソフトで直交異方性材料を定義する時、9つの定数をそのまま入力する欄があるんですか?

🎓

はい、主要なソフトウェアにはそのためのインターフェースがあります。例えば、Abaqusでは「Engineering Constants」として、

$$ E_1, E_2, E_3, \nu_{12}, \nu_{13}, \nu_{23}, G_{12}, G_{13}, G_{23} $$
の9つを直接入力します。Ansys Mechanicalでも同様に「Orthotropic Elasticity」の項目で定義可能です。ここで最も重要なのは、これらの定数が「材料座標系」に対して定義されていることです。

🧑‍🎓

「材料座標系」とは、具体的にどう設定するんですか?部品の形状が複雑だと、繊維方向(1方向)を全体で定義するのは難しそうです。

🎓

それが直交異方性解析のキーポイントです。ソフトウェアごとに方法は異なりますが、代表的なのは「局所座標系の割り当て」です。Ansysでは「Coordinate Systems」を定義して材料に割り当て、Abaqusでは「Orientation」という機能を使います。例えば、航空機の複合材パネルでは、パネルの面内方向を1,2方向、法線方向を3方向とする座標系を、要素ごとに定義します。曲がったシェルでは、シェルの法線が自動的に3方向になるように設定します。

🧑‍🎓

要素ごとに座標系が違う場合、ソルバーは内部でどう処理しているんですか?全体の剛性マトリクスを作る時に、毎回変換している?

🎓

その通りです。要素剛性マトリクス

$$ [k_e] $$
は、まず材料座標系(ローカル系)で計算されます。その後、要素の材料座標系から全体座標系(グローバル系)への方向余弦マトリクス
$$ [T] $$
を用いて変換
$$ [k_e]_{global} = [T]^T [k_e]_{local} [T] $$
が行われ、全体剛性マトリクスに組み込まれます。この変換により、繊維方向が部品の形状に沿って曲がっているような複雑なケースも解析可能になります。

実践ガイド

解析ワークフローと検証

🧑‍🎓

実際に直交異方性材料で解析を始める時、最初に何を確認すべきですか?

🎓

まずは「材料定数の物理的妥当性」のチェックです。入力した9定数が熱力学的に許容されるか、つまり、どんな負荷をかけてもひずみエネルギーが正になるか確認します。具体的には、ヤング率とせん断弾性係数は全て正

$$ E_i > 0, G_{ij} > 0 $$
であり、ポアソン比は
$$ |\nu_{ij}| < \sqrt{(E_i / E_j)} $$
を満たす必要があります。これを怠ると、解析が発散したり、非現実的な結果が出たりします。

🧑‍🎓

材料定数はOKでも、結果がおかしい場合、次に疑うべきポイントは?

🎓

ほぼ間違いなく「材料座標系の設定ミス」です。これが直交異方性解析で最も多いエラー原因です。検証方法として、単純な引張りや曲げをかけた時、剛性の高い方向(繊維方向)に応力が集中しているか、変形が小さいかを確認します。例えば、CFRPの板をモデル化し、繊維方向(1方向)に引張りをかけた時の変位は、垂直方向(2方向)に引張りをかけた時の1/10以下になるはずです。この差が出なければ座標系が間違っています。

🧑‍🎓

複合材積層板(ラミネート)のように、層ごとに繊維角度が違う場合はどう扱うんですか?

🎓

その場合は、各単層(ラミナ)を直交異方性材料として定義し、積層順序(Ply Stacking Sequence)と各層の繊維角度(例えば、0°, 45°, -45°, 90°)を指定する「複合材モデル」機能を使います。Abaqusの「Composite Layup」、Ansysの「ACP (ANSYS Composite PrepPost)」や「Shell Layup」が該当します。これらのツールは、各層の材料座標系を指定された角度だけ回転させ、積層全体の等価剛性を自動計算してくれます。

ソフトウェア比較

主要ソフトの機能と特徴

🧑‍🎓

Ansys、Abaqus、COMSOLで、直交異方性材料の扱いに違いはありますか?

🎓

基本的な材料定義の機能はどれも備えていますが、強みの分野が異なります。Ansysは、構造解析(Mechanical)に加え、専用の複合材前後処理モジュール「ACP」が強力です。特に自動車・航空宇宙産業で広く使われ、積層板の詳細な破壊解析までカバーします。Abaqusも「Composite Layup」機能が充実しており、非線形挙動(損傷進展など)の高度なモデリングに定評があります。

🧑‍🎓

COMSOLはどういう位置付けなんですか?

🎓

COMSOL Multiphysicsの最大の特徴は、直交異方性を「構造」以外の物理場に簡単に拡張できる点です。例えば、熱伝導解析で、繊維方向と垂直方向で熱伝導率が異なる「熱直交異方性」を定義できます。電気伝導や拡散現象でも同様です。また、材料座標系の定義が非常に柔軟で、任意の数式や関数で方向を指定できるため、複雑な織物構造などのモデリングに適しています。

🧑‍🎓

無料や低価格のCAEソフト(例えば、CalculiX、Code_Aster)では直交異方性は扱えるんですか?

🎓

はい、扱えます。例えば、CalculiX(フリーのソルバー)では入力ファイル(.inp)に*ELASTIC, TYPE=ORTHOキーワードで9定数を定義できます。Code_Aster(Salome-MECA環境)でも、材料定義コマンド(ELAS)でORTHキーワードを使用します。ただし、商用ソフトのようなグラフィカルで直感的な積層板定義機能や高度な破壊則は標準では備えていない場合が多く、ユーザーがスクリプトや手動で詳細を設定する必要があります。

トラブルシューティング

よくあるエラーとその原因

🧑‍🎓

解析を実行したら「負の対角項」や「非正定値マトリクス」というエラーでソルバーが停止しました。材料定数は正の値しか入れていないのですが。

🎓

それは、個々の定数が正でも、組み合わせが熱力学的に不安定(ひずみエネルギーが負になりうる)な状態だからです。特にポアソン比の値が問題です。先述の条件

$$ |\nu_{12}| < \sqrt{E_1/E_2} $$
をチェックしてください。例えば、
$$ E_1=150GPa, E_2=10GPa $$
の場合、
$$ \sqrt{150/10} \approx 3.87 $$
なので、
$$ \nu_{12} $$
は3.87未満でなければなりません。通常ポアソン比は0.5以下ですが、直交異方性ではこの上限がもっと緩和される点に注意が必要です。

🧑‍🎓

応力集中部の結果を見ると、せん断応力が異常に高くなっています。材料のせん断強度を超える値が出ているのですが、信頼できますか?

🎓

それは「直交異方性材料における応力特異点」の典型的な現象かもしれません。異方性材料、特に異方性比(

$$ E_1/E_2 $$
)が大きい材料では、切欠きや孔の縁で、等方性材料よりもはるかに高い応力勾配が生じることが理論的に知られています。メッシュを細かくするほどせん断応力が増加し、収束しない場合もあります。このような場合は、線形弾性解析の応力結果をそのまま強度評価に使うのは危険で、平均化応力を評価するか、初めから損傷則を組み込んだ非線形解析を行う必要があります。

🧑‍🎓

シェル要素を使った積層板解析で、面外せん断変形が無視できないと言われました。どう対処すれば?

🎓

それは「厚みせん断変形」の問題です。通常の薄板理論(Kirchhoff-Loveの仮定)は面外せん断変形を無視しますが、複合材は

$$ G_{13}, G_{23} $$
(面外せん断剛性)が比較的小さいため、この影響が無視できません。対処法は二つあります。第一に、厚みせん断変形を考慮したシェル要素(AbaqusのS4R, AnsysのSHELL181の「Shear Correction」ON)を使用すること。第二に、より高精度を求めるなら、連続体シェル要素(Continuum Shell)やソリッド要素を用いて厚さ方向も離散化することです。ただし計算コストは大幅に増加します。

この記事の評価
ご回答ありがとうございます!
参考に
なった
もっと
詳しく
誤りを
報告
参考になった
0
もっと詳しく
0
誤りを報告
0
Written by NovaSolver Contributors
Anonymous Engineers & AI — サイトマップ
プロフィールを見る