/INTER/TYPE20（廃止）

エッジ-エッジの接触も使用できます。ペナルティ剛性は一定であるため、時間ステップの影響は受けません（標準ペナルティ剛性の場合）。このインターフェースは、インターフェースTYPE3、TYPE5、TYPE7、TYPE11またはTYPE19に代えて使用できます。インターフェースは、基本的に1つまたは2つのサーフェスによって定義されます。1つのサーフェスのみが使用される場合、このサーフェスは自己接触を行います。2つのサーフェスが定義される場合、サーフェス2の節点がサーフェス1に衝撃を与えます。対称処理をアクティブにすることができます。サーフェス1および2のエッジは、接触において考慮されます。サーフェスに節点を追加できます。

フォーマット

定義

境界条件の非アクティブ化フラグ：IBC

定義

コメント

1. 最も単純な入力は、1つの自己接触サーフェスsurf_ID₁のみを入力する方法です。対称性処理のフラグI_symは節点-サーフェス接触に使用されますが、エッジ-エッジ接触は常に対称です。
2. インターフェースTYPE7またはTYPE11の入力を模擬するには、以下のように入力します：

   TYPE20でTYPE7の入力を模擬	TYPE20でTYPE11の入力を模擬
   (TYPE20) surf_ID1 = surf_IDm (TYPE7) (TYPE20) grnd_ID = grnd_IDs (TYPE7)	(TYPE20)line_ID1 = line_IDs (TYPE11) (TYPE20)line_ID2 = line_IDm (TYPE11)
   (TYPE20)line_ID1 = 0 (TYPE20)line_ID2 = 0	(TYPE20) surf_ID1 = 0 (TYPE20) grnd_ID = 0
   (TYPE20) surf_ID2 = 0 (TYPE20)Isym = 2 (TYPE20)Iedge = 0	(TYPE20) surf_ID2 = 0 (TYPE20)Isym = 0 (TYPE20)Iedge = 0

3. SPMDの場合、surf_ID_mによって定義される各メインセグメントを1つの要素に関連付ける必要があります（ボイド要素の場合もあります）。
4. フラグI_bagについては、モニター体積オプション（モニター体積（エアバッグ））をご参照ください。
5. フラグ I_del = 1のCPUコストはI_del = 2よりも高くなります。
6. I_gap = 0の場合、Gap₀に使用されるデフォルト値。次のように計算されます：(1)
   $$Ga{p}_0=min\left(t,\frac{l}{10},\frac{l_{min}}{2}\right)$$
   ここで、

   t

   メインシェル要素の平均板厚

   l

   メインソリッド要素の平均辺長

   l_min

   すべてのメインセグメント（シェルまたは3次元ソリッド）の中の最小辺長

7. I_gap = 1の場合、ギャップは衝撃のそれぞれに対して下記の式で計算されます：
   (2)

   $$g_s+g_m$$
   ここで、

   • g_m： メイン要素のギャップ：

     $g_m=\frac{t}{2}$ 、ここでtは、シェル要素に対するメイン要素の板厚

     g_m = 0（3次元ソリッド要素の場合）

   • g_s： セカンダリ節点のギャップ：

     g_s = 0（セカンダリ節点がどの要素にも結合されていないか、3次元ソリッド要素またはスプリング要素にのみ結合されている場合）

     $g_s=\frac{t}{2}$ ここでtは、セカンダリ節点に結合されているシェル要素の最大板厚です。

     $g_s=\frac{1}{2}\sqrt{S}$ （トラスおよびビーム要素の場合）。ここで、Sは要素の断面とします。

     セカンダリ節点が複数のシェルおよび/またはビームまたはトラスに結合されている場合は、計算された中で最も大きいセカンダリギャップが使用されます。

     シェル要素のフリーエッジが接触している場合、I_gapは、以下のようにフリーエッジの境界シェルのギャップをシフトできます。

     
     
     図 1.

8. 境界条件の非アクティブ化は、セカンダリ節点グループ（surf_ID_s）に適用されます。
9. 節点がスプリング要素に属している場合にInacti = 3であると、初期エネルギーが発生することがあります。
   Inacti = 5：

   
   
   図 2.

10. 最大貫通値は、実際のギャップ（可変ギャップを含む）の割合として設定されます。

    $Penetration\ge Fpenmax\cdot Gap$

    セカンダリ節点の初期貫通が計算された最大値（Fpenmax）より大きい場合、その節点はインターフェースで非アクティブとなります（節点剛性の非アクティブ化）。

11. 1つの節点が2つのサーフェスに同時に属することができます。
12. 剛性係数の値に制限はありません（ただし、値が1.0より大きいと、最初の時間ステップが短くなる場合があります）。
13. 摩擦定式化の場合：
    • 摩擦フラグ = 0（デフォルト）の場合は、従来の静的摩擦係数の定式化が使用されます：

      $F_t\le \mu \cdot F_n$ ここで、 $\mu =\mathit{Fric}$ （ $\mu$ は、クーロン摩擦係数）

    • フラグI_fric > 0の場合、新しい摩擦モデルが導入されます。この場合、摩擦係数は関数 $\mu =\text{μ}(\rho ,V)$ によって設定されます。
      ここで、

      $\rho$

      メインセグメントの垂直抗力の圧力

      $V$

      セカンダリ節点の接線速度

14. 現在は、係数C₁～C₆を使用して、新しい摩擦定式化の可変摩擦係数 $\mu$ を定義しています。
    以下の定式化を使用できます：

    • I_fric = 1（汎用の粘性摩擦則）：(3)
      $$\mu =\mathit{Fric}+C_1\cdot p+C_2\cdot V+C_3\cdot p\cdot V+C_4\cdot p^2+C_5\cdot V^2$$
    • I_fric = 2（修正Darmstad則）：(4)
      $$\mu =Fric+C_1\cdot e^{\left(C_{2}V\right)}\cdot p^2+C_3\cdot e^{\left(C_{4}V\right)}\cdot p+C_5\cdot e^{\left(C_{6}V\right)}$$
    • I_fric = 3（Renard則）：

      $\mu =C_1+\left(C_3-C_1\right)\cdot \frac{V}{C_5}\cdot \left(2-\frac{V}{C_5}\right)$ の、 $V\in \left[0,C_5\right]$

      $\mu =C_3-\left(\left(C_3-C_4\right)\cdot {\left(\frac{V-C_5}{C_6-C_5}\right)}^2\cdot \left(3-2\cdot \frac{V-C_5}{C_6-C_5}\right)\right)$ 、右記の場合； $V\in \left[C_5,C_6\right]$

      $\mu =C_2-\frac{1}{\frac{1}{C_2-C_4}+{\left(V-C_6\right)}^2}$ 、右記の場合； $V\ge C_6$

      ここで、

      $C_1={\mu }_s$ $C_2={\mu }_d$

      $C_3={\mu }_{\max }$ $C_4={\mu }_{\min }$

      $C_5=V_{\mathit{cr}1}$ $C_6=V_{cr2}$

    • 第1臨界速度 $V_{cr1}=C_5$ は、0以外にする必要があります（ $C_5\ne 0$ ）。
    • 第1臨界速度 $V_{cr1}=C_5$ は、第2臨界速度 $V_{cr2}=C_6$ より小さくする必要があります（ $C_5<C_6$ ）。
    • 静止摩擦係数 $C_1$ と動摩擦係数 $C_2$ は、最大摩擦 $C_3$ より小さくする必要があります（ $C_1\le C_3$ かつ $C_2\le C_3$ ）。
    • 最小摩擦係数 $C_4$ は、静止摩擦係数 $C_1$ および動摩擦係数 $C_2$ より小さくする必要があります（ $C_4\le C_1$ かつ $C_4\le C_2$ ）。
15. 摩擦フィルタリング：
    I_filtr ≠ 0の場合は、接線力がフィルタを使用して以下のようにスムージングされます：(5)

    $$F_t=\alpha \cdot {F^{\prime }}_t+\left(1-\alpha \right)\cdot {F^{\prime }}_t{}^{-1}$$
    ここで、α係数は、以下のように計算されます：

    • I_filtr= 1の場合、 $\alpha =X_{\mathit{freq}}$ 、単純な数値フィルター
    • I_filtr = 2 ➤ $\alpha =\frac{2\cdot \pi }{X_{\mathit{freq}}}$ の場合、標準の-3dBフィルター。ここで、 $X_{\mathit{freq}}=\frac{dt}{T}$ で、Tはフィルタリング期間
    • I_filtr = 3 ➤ $\alpha =2\cdot \pi \cdot X_{freq}\cdot dt$ の場合、カット周波数X_freqを使用する標準の-3dBフィルター

      フィルタリング係数X_freqは、0～1の値にする必要があります。

16. 摩擦ペナルティ定式化I_form
    • I_form = 1（デフォルト）の粘性定式化の場合、摩擦力は次のとおりです：(6)
      $$F_t=\min \left(\mu F_n,F_{adh}\right)$$

      このとき、粘着力は以下のように計算されます：

      $F_{adh}=C\cdot V_t$ ここで、 $C=VI{S}_F\cdot \sqrt{2Km}$

    • I_form = 2の粘性定式化の場合、摩擦力は次のとおりです：(7)
      $$F_t^{new}=\min \left(\mu F_n,F_{adh}\right)$$

      このとき、粘着力は以下のように計算されます：

      $F_{adh}=F_t^{old}+\text{Δ}{F}_t$ ここで、 $\text{Δ}{F}_t=K\cdot V_t\cdot dt$

      ここで、V_tは、接触接線速度です。