直交異方性 — CAE用語解説
理論と物理
直交異方性の物理的意味
「直交異方性」という言葉をよく見かけますが、等方性と具体的に何が違うんですか?
根本的な違いは、材料の特性が方向によって異なるかどうかです。等方性材料は、どの方向から引っ張ってもヤング率やポアソン比が同じです。直交異方性材料は、互いに直交する3つの主方向(例えば、繊維方向、面内直交方向、厚さ方向)で特性が異なります。例えば、CFRP(炭素繊維強化プラスチック)では、繊維方向のヤング率は約150 GPaなのに対し、繊維に垂直な方向では10 GPa程度と、10倍以上の差が出ます。
3つの方向で違うのは分かりました。では、その特性を数式で表すと、等方性と比べてパラメータが増えるということですか?
その通りです。等方性の線形弾性は、ヤング率Eとポアソン比νの2つで記述できます。一方、直交異方性の完全な線形弾性構成則には、独立な材料定数が9つ必要です。3方向のヤング率
9つもあると、実験で全部求めるのは大変じゃないですか?実際の設計ではどうしているんですか?
良い指摘です。実際、厚さ方向(3方向)の特性
数値解法と実装
CAEソフトウェアでの入力と座標系
CAEソフトで直交異方性材料を定義する時、9つの定数をそのまま入力する欄があるんですか?
はい、主要なソフトウェアにはそのためのインターフェースがあります。例えば、Abaqusでは「Engineering Constants」として、
「材料座標系」とは、具体的にどう設定するんですか?部品の形状が複雑だと、繊維方向(1方向)を全体で定義するのは難しそうです。
それが直交異方性解析のキーポイントです。ソフトウェアごとに方法は異なりますが、代表的なのは「局所座標系の割り当て」です。Ansysでは「Coordinate Systems」を定義して材料に割り当て、Abaqusでは「Orientation」という機能を使います。例えば、航空機の複合材パネルでは、パネルの面内方向を1,2方向、法線方向を3方向とする座標系を、要素ごとに定義します。曲がったシェルでは、シェルの法線が自動的に3方向になるように設定します。
要素ごとに座標系が違う場合、ソルバーは内部でどう処理しているんですか?全体の剛性マトリクスを作る時に、毎回変換している?
その通りです。要素剛性マトリクス
実践ガイド
解析ワークフローと検証
実際に直交異方性材料で解析を始める時、最初に何を確認すべきですか?
まずは「材料定数の物理的妥当性」のチェックです。入力した9定数が熱力学的に許容されるか、つまり、どんな負荷をかけてもひずみエネルギーが正になるか確認します。具体的には、ヤング率とせん断弾性係数は全て正
材料定数は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キーワードを使用します。ただし、商用ソフトのようなグラフィカルで直感的な積層板定義機能や高度な破壊則は標準では備えていない場合が多く、ユーザーがスクリプトや手動で詳細を設定する必要があります。
トラブルシューティング
よくあるエラーとその原因
解析を実行したら「負の対角項」や「非正定値マトリクス」というエラーでソルバーが停止しました。材料定数は正の値しか入れていないのですが。
それは、個々の定数が正でも、組み合わせが熱力学的に不安定(ひずみエネルギーが負になりうる)な状態だからです。特にポアソン比の値が問題です。先述の条件
応力集中部の結果を見ると、せん断応力が異常に高くなっています。材料のせん断強度を超える値が出ているのですが、信頼できますか?
それは「直交異方性材料における応力特異点」の典型的な現象かもしれません。異方性材料、特に異方性比(
シェル要素を使った積層板解析で、面外せん断変形が無視できないと言われました。どう対処すれば?
それは「厚みせん断変形」の問題です。通常の薄板理論(Kirchhoff-Loveの仮定)は面外せん断変形を無視しますが、複合材は
関連トピック
なった
詳しく
報告