Tied接触の失敗

接着剤よりもはるかに強固で、数学的には「拘束条件」に近いです。主従面の対応する節点の自由度を等価にする、つまり剛性マトリクスの中で強制的に同じ変位・回転をさせます。例えば、主面の節点mと従面の節点sがあれば、その変位関係は

その通りです。現実のボルト締結やスポット溶接でも、局所的には高い剛性ですが有限です。Tied接触はこの局所的な柔らかさを無視するため、応力集中が不自然に高くなったり、接触面直下の要素が極端な変形を強いられて「歪みエネルギー発散」を起こす原因になります。特に、面の曲率が大きく異なる部分や、メッシュサイズが急激に変化する境界でこの問題が顕著です。

メッシュを細かくするほど、計算される歪みエネルギーが無限大に発散してしまう非物理的な現象です。理論的には、一点に無限大の応力が生じることを意味します。Tied接触では、主面と従面の節点が完全に同一の動きを強制されるため、接触エッジ部分で変形の不連続性が生じ、これが特異点となって発散を引き起こします。対策として、Abaqusの「Tie接触」には「位置の許容値」設定があり、一定距離内の面を自動的に拘束しますが、この許容値を大きくしすぎると意図しない部分がくっついてしまいます。

主に2つの方法があります。「ラグランジュ乗数法」と「ペナルティ法」です。Ansys Mechanicalではデフォルトで「MPC（多点拘束）アプローチ」が使われ、これはラグランジュ乗数法の一種で、拘束条件を厳密に満たします。一方、Abaqus/Explicitでは安定性のために「キネマティック（運動学的）カップリング」や「ペナルティ法」が多用されます。ペナルティ法では、拘束条件

全体の剛性マトリクスが数値的に「硬い」状態になり、連立一次方程式のソルバーが収束しにくくなります。反復解法では収束までの反復回数が異常に増え、直接法では丸め誤差の影響が大きくなり、結果として「ソルバーが停止する」「負のピボットエラーが出る」「応力結果がNaN（非数）になる」といった失敗に直結します。Ansysでは「ピボット警告」、Abaqusでは「ゼロピボット」や「過度の歪み」というメッセージでこれが報告されます。

これがTied接触失敗の最大の原因の一つです。多くのソルバーは「サーフェス・トゥ・サーフェス」のアルゴリズムを使い、従面の節点を主面の要素面上に投影し、その要素の形状関数を使って重み付け平均した主面節点の変位を従面節点に割り当てます。式で書くと、従面節点sの変位は、投影先の主面要素の節点m1, m2, m3...の変位を用いて

以下のチェックリストを必ず実行してください。1. **幾何クリアランス**: 拘束したい面同士の距離が、ソフトウェアのデフォルトの位置許容値（例：Abaqusでは「スレッショルド」、デフォルトはモデル全体サイズの0.05%程度）以内か。2. **メッシュ整合性**: 接触境界でのメッシュサイズ比を5倍以内に収める。3. **面の向き**: 面法線ベクトルが互いに向き合っているか。4. **エッジ/頂点の接触**: エッジが面に、または頂点が面に接触する設定は避け、面-面接触のみとする。

その場合はTied接触の使用を再考すべきです。代替案として、(1) **サブモデリング**を用い、全体モデルでは粗いメッシュで結合し、局部モデルで細かいメッシュに切り替える。(2) **ボルト/リベット/スポット溶接要素**（Ansysの`CONTA178`、Abaqusの`FASTENER`）で離散的な結合点を定義する。(3) どうしてもTied接触を使うなら、粗いメッシュ側を「主面」に設定する。一般的に、より剛性が高い、またはメッシュが粗い面を主面に選択することが原則です。

必ず「単体剛性テスト」を行ってください。具体的には、Tied接触を設定したモデルに、片方を完全固定し、もう片方の端面にごく小さな強制変位（例えば0.001 mm）を与える静解析を実行します。この時、以下のことを確認します：計算が正常に収束するか、接触面付近の応力が異常な高値（降伏応力を遥かに超える等）ではないか、反力がつり合っているか。このテストで失敗する場合は、本解析を実行する意味がありません。

コンセプトは同じですが、実装とデフォルト設定に重要な違いがあります。

これは非常に重要です。「位置の許容値」は、従面節点を探すために主面から広げる検索距離です。例えば0.1mmと設定すれば、主面から0.1mm以内の従面節点を探し出します。一方、「表面のスレッショルド」は、見つかった従面節点のうち、実際に拘束から除外する割合です。デフォルトは0.05（5%）で、これは「拘束すると最も歪みエネルギーが大きくなる上位5%の節点を、自動的に拘束対象から外す」機能です。位置許容値を大きくしすぎると遠くの面までくっつき、表面スレッショルドを大きくしすぎると肝心な部分の拘束が外れてしまいます。実務では、まず位置許容値を幾何クリアランスよりわずかに大きい値（例：0.05mm）に手動設定し、表面スレッショルドはデフォルトのまま試すのが基本です。

まず**Tied接触面の直下の要素**を疑ってください。対策は以下の順で試します：

**絶対に信用してはいけません**。それは先ほど話した数値的特異性による「見かけ上の応力」です。特に、メッシュを細かくするほどその値が増加するなら、ほぼ間違いありません。このような結果を設計判断に使うのは危険です。対策としては、(1) 接触面から少し離れた位置（例えば、要素1~2個分内側）の応力を評価する「パスプロット」を作成する。(2) ノッチの応力評価で使われる「応力グラディエント法」や「構造応力法」を参考に、特異点の影響を除外した方法で応力を平均化して評価する。実務では、Tied接触面の応力そのものではなく、そこから伝達される力やモーメントが全体としてつり合っているかを確認する方が重要です。

それは「重複拘束」または「過剰拘束」が原因です。一つの従面節点が、複数の主面要素（または別々のTied接触定義）から同時に拘束されようとすると、ソルバーは矛盾した指示を受けて計算不能に陥ります。例えば、3つのプレートを重ねて両端のプレートを中央にTied接触する場合、中央のプレートは上下両方から拘束されるため過剰拘束になります。対策は、**拘束の連鎖を作らない**ことです。この例なら、「プレート1 → プレート2」「プレート2 → プレート3」というように、隣接する部品同士のみでTied接触を定義し、連鎖的に結合させます。Ansysの「Contact Tool」やAbaqusの「Interaction Manager」で接触ペアを可視化し、ネットワークがループしていないかを確認する習慣をつけましょう。