せん断弾性係数 — CAE用語解説

良い質問だ。せん断弾性係数Gは、材料が「ずれ変形」に対してどれだけ抵抗するかを表す定数だ。立方体の上面を横方向に引っ張り、下面を固定したときの、ずれ角γとせん断応力τの比例定数として定義される。式で書くと、

その通り、等方性の線形弾性体では、三つの材料定数は独立ではなく、次の関係で結ばれている。

鋭い指摘だ。等方性材料ではせん断弾性係数は1つだが、異方性材料では複数必要になる。最も一般的な直交異方性材料では、3つの主せん断弾性係数

その核心は、せん断ひずみの計算が変位の「微分」に依存する点にある。一次要素（線形要素）では変位場が線形なので、せん断ひずみは要素内で一定値になってしまう。これは実際の応力分布を過小評価する。特に「せん断ロッキング」という問題が発生し、曲げ変形時に不自然に硬い結果が出る。これを防ぐために、二次要素（2次アイソパラメトリック要素）や、一次要素でも低減積分や非適合モードを導入した要素が開発された。

数学的に言うと、せん断ひずみエネルギー項の数値積分を「過剰評価」することが原因だ。完全積分（2x2）の一次要素では、この過剰評価が顕著で、要素が曲げ変形を起こす際に莫大なせん断ひずみエネルギーが発生したと計算され、硬くなる。低減積分（1点）では、この積分評価が要素中心の1点だけになるため、過剰評価が起こりにくくなる。ただし、代わりに「アワーモグラスモード」というゼロエネルギー変形が発生する可能性があり、安定化処理が必要になる。

実務では以下の手順を推奨する。まず、要素選択：AbaqusならC3D8R（低減積分一次要素）より、C3D10（二次四面体）やC3D20R（二次六面体低減積分）を優先する。AnsysならSOLID186（二次六面体）が信頼性が高い。後処理では、節点平均化されたせん断応力をそのまま信用するのではなく、要素内の分布を確認すること。特に応力集中部では、メッシュを細かくして収束性を確認する「メッシュコンバージェンススタディ」が必須だ。また、せん断応力は座標系に依存するので、主せん断応力やトレースを見ることも有効だ。

まず、材料規格やメーカーのデータシートを確認すること。JIS G 4051（機械構造用炭素鋼）などの規格では、EとνからGを計算して使うことが暗黙の了解だ。ソフトウェア（Ansys Mechanical, Abaqus/CAEなど）の材料プロパティ入力画面では、Eとνを入力する欄があり、Gの欄は空白でも自動計算されるか、入力しなくても内部で

もちろんある。代表的なのは「ねじり試験」だ。丸棒試験片の一端を固定し、他端にトルクを加えてねじれ角を測定する。せん断応力τとせん断ひずみγの関係からGを直接求められる。ASTM E143という標準試験法がある。もう一つの方法は、片持ち梁の先端に荷重をかけてたわみを測る「三点曲げ試験」から間接的に求める方法だ。この場合、たわみの理論式にEとGの両方が含まれるため、Eを別の引張試験で求めておけば、Gを逆算できる。

重要な視点だ。チェックポイントは三つ。第一に、出典が明確か（JIS、ISO、ASTM規格値か、メーカー公表値か）。第二に、Eとνの値も同時に記載されているか。先の式で計算したGとデータベースのGが大きく異なる場合は、どちらかが間違っている可能性が高い。第三に、温度依存性のデータがあるか。特にプラスチックでは、ガラス転移温度付近でGが1～2桁急激に低下する。信頼性の高いデータベースとしては、CES EduPack (Granta Design) や、各CAEベンダー（Ansys Granta MIなど）が提供するデータベースを参照することを勧める。

基本的な考え方は同じだが、インターフェースとオプションに違いがある。Abaqus/CAEでは、材料プロパティマネージャーの「Elastic」タイプで、Typeを「Isotropic」とし、「Young's Modulus」と「Poisson's Ratio」を入力する。Gを直接入力するフィールドはデフォルトでは表示されない（隠しオプションのようなものだ）。一方、Ansys Mechanicalでは、Engineering Dataで「Isotropic Elasticity」を追加し、「Young's Modulus」と「Poisson's Ratio」を入力する。こちらもGの入力欄はない。両者とも、内部でEとνからGを計算している。

ここで違いが顕著になる。Abaqusでは、「Engineering Constants」を用いた直交異方性定義が強力だ。E1, E2, E3, ν12, ν13, ν23, G12, G13, G23の9定数を直接入力する。超弾性（ゴム）では、せん断弾性は初期せん断弾性率μとして、Mooney-Rivlin定数などから定義される。Ansysも同様に直交異方性を定義できるが、入力インターフェースがマトリクス形式だったりする。COMSOL Multiphysicsは「線形弾性材料」ノードで、等方性、直交異方性、完全異方性と柔軟に選択でき、各Gijを直接入力するUIが非常に明確だ。用途によって使い分けが発生する部分である。

大きな違いがある。Abaqus/Viewerでは、「S12」「S13」「S23」がせん断応力成分として直接出力され、カスタム場出力で「Max. In-Plane Shear」や「Tresca Stress」も簡単にプロットできる。Ansys Mechanicalの「Equivalent Stress」はデフォルトでフォンミーゼス応力だが、プローブ機能で任意のせん断応力成分を抽出できる。また、Ansys APDL（コマンド）では、「RSYS」コマンドで結果座標系を変更できるので、ボルトの軸に沿ったせん断応力を簡単に見られる。一方、COMSOLは「ストレス」ノード下に「せん断応力」コンポーネントが最初から用意されており、直感的だ。いずれも、せん断応力は符号（方向）を持つので、絶対値や大きさを確認するのがコツだ。

まず疑うべきは「せん断ロッキング」だ。特に、板厚方向に1層しか要素がないソリッド要素（一次要素）で薄板をモデル化すると高確率で発生する。理論的には薄板曲げではせん断変形は無視できる（Euler-Bernoulli梁理論）が、一次要素ではせん断変形が支配的になり、異常に硬く、かつせん断応力が過大評価される。対策は、板厚方向に最低3層以上の要素を使う、またはシェル要素（S4Rなど）を使用することだ。シェル要素は板の曲げ理論に基づいており、せん断変形を適切に扱える。

そのような挙動は「コヒーシブゾーンモデル（CZM）」または「接触にせん断剛性を付与する」方法でモデル化する。Abaqusでは「Cohesive Elements」または「Surface-based Cohesive Behavior」を使用し、法線方向とせん断方向の剛性（および強度）を別々に定義する。せん断剛性Ksは、接着剤のせん断弾性係数Gを厚みtで割った

これは重大なエラーで、入力した材料定数の組み合わせが物理的に不可能であることを意味する。直交異方性材料では、材料定数が以下の不等式を満たす必要がある（熱力学的安定性の条件）。