/INTER/LAGDT/TYPE7

内容

フォーマット

定義

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

定義

コメント

1. フラグI_bagについては、モニター体積オプション（モニター体積（エアバッグ））をご参照ください。
2. フラグ I_del = 1のCPUコストはI_del = 2よりも高くなります。
3. I_gap = 2の場合、可変ギャップは次のように計算されます：(1)
   $$\text{max}\left\{{\mathit{Gap}}_{\text{min}},\text{min}\left[{\mathit{Fscale}}_{\mathit{gap}}\cdot \left(g_s+g_m\right),{\mathit{Gap}}_{\text{max}}\right]\right\}$$

   I_gap ≠ 2の場合、行4で指定された値は無視されます。

4. 接触剛性は以下のように計算されます：
   • I_stf = 0の場合の剛性

     $K=K_m$

   • I_stf > 1の場合の剛性(2)
     $$K=\text{max}\left[S{t}_{\text{min}},\text{min}\left(S{t}_{\text{max}},K_n\right)\right]$$
     ここで、

     • $K_n$ は右記の両方から次のように計算されます； メインセグメントの剛性 $K_m$ とセカンダリ節点の剛性 $K_s$

       I_stf = 2、 $K_n=\frac{K_m+K_s}{2}$

       I_stf = 3、 $K_n=\text{max}\left(K_m,K_s\right)$

       I_stf = 4、 $K_n=\text{min}\left(K_m,K_s\right)$

       I_stf = 5、 $K_n=\frac{K_m\cdot K_s}{K_m+K_s}$

     • $K_m$ はメインセグメントの剛性で、次のように計算されます：

       メインセグメントがシェル上に存在するか、シェルとソリッドによって共有されている場合(3)

       $$K_m=\mathit{Stfac}\cdot 0.5\cdot E\cdot t$$
       メインセグメントがソリッド上に存在する場合。(4)

       $$K_m=B\cdot \frac{S^2}{V}$$
       ここで、

       $S$

       セグメント面積

       $V$

       ソリッドの体積

       $B$

       体積弾性率

       $K_s$

       インターフェースTYPE7として考慮される相当節点剛性で、次のように計算されます：

       • 節点がシェル要素に結合されている場合：(5)
         $$K_s=\mathit{Stfac}\cdot \frac{1}{2}\cdot E\cdot t$$
       • 節点がソリッド要素に結合されている場合：(6)
         $$K_s=\mathit{Stfac}\cdot B\cdot \sqrt[3]{V}$$

   剛性係数の値に対する制限はありません（ただし、値が1.0より大きいと、初期時間ステップが短くなる場合があります）。

5. I_stf < 1の場合、行5で指定された値は無視されます。
6. Gap_minのデフォルト値は、下記の最大値として計算されます：(7)
   $$Ga{p}_{\text{min}}=\text{min}\left(t,\frac{l}{10},\frac{l_{\text{min}}}{2}\right)$$
   ここで、

   $t$

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

   $l$

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

   $l_{\min }$

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

7. ギャップは衝撃のそれぞれに対して下記の式で計算されます：(8)
   $${\mathit{Fscale}}_{gap}\cdot \left(g_s+g_m\right)$$
   ここで、

   • $g_m$ ： メイン要素のギャップ：(9)
     $$g_m=\frac{t}{2}$$

     ここで、 $t$ ： シェル要素に対するメイン要素の板厚

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

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

     $$g_s=\frac{t}{2}$$

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

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

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

   可変ギャップは常にGap_min以上です。

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

   Inacti = 5は、エアバッグシミュレーション展開に推奨されています。

   インターフェースに対する高周波数の影響を回避するためには、Inacti =5ではなく、Inacti = 6が推奨されます。

   
   
   図 1.

10. ソートアルゴリズムを高速化するには、ソート係数Bumultを使用します。
11. Bumultのデフォルト値は、150万個を超える節点を持つモデルでは0.30に、250万個を超える節点を持つモデルでは0.40に自動的に増やされます。
12. 1つの節点が2つのサーフェスに同時に属することができます。
13. 剛性係数の値に制限はありません（ただし、値が1.0より大きいと、最初の時間ステップが短くなる場合があります）。
14. 摩擦定式化の場合
    • 摩擦フラグI_fric = 0（デフォルト）の場合は、従来の静的摩擦係数の定式化が使用されます：

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

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

      $\rho$

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

      $V$

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

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

    • I_fric = 1（汎用の粘性摩擦則）：(11)
      $$\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則）：(12)
      $$\mu =C_1\cdot e^{(C_{2}V)}\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則）：(13)
      $$\mu =C_1+\left(C_3-C_1\right)\cdot \frac{V}{C_5}\cdot \left(2-\frac{V}{C_5}\right)\hspace{0.17em}\text{if}\hspace{0.17em}V\in \left[0,C_5\right]$$
      (14)
      $$\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)\hspace{0.17em}\text{if}\hspace{0.17em}V\in \left[C_5,C_6\right]$$
      (15)
      $$\mu =C_2-\frac{1}{\frac{1}{C_2-C_4}+{\left(V-C_6\right)}^2}\hspace{0.17em}\text{if}\hspace{0.17em}V\ge C_6$$

      ここで、

      $C_1={\mu }_s$

      $C_2={\mu }_d$

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

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

      $C_5=V_{cr1}$

      $C_6=V_{cr2}$

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

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

    • I_filtr = 1 ➤ $\alpha =X_{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_{\mathit{freq}}\cdot \mathit{dt}$ の場合、カット周波数X_freqを使用する標準の-3dBフィルター

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

17. 摩擦ペナルティ定式化I_form
    • I_form = 1（デフォルト）の粘性定式化の場合、摩擦力は次のとおりです：(17)
      $$F_t=\min \left(\mu F_n,F_{\text{adh}}\right)$$
    • このとき、粘着力は以下のように計算されます：(18)
      $$F_{\mathit{adh}}=C\cdot V_t\text{with}C={\mathit{VIS}}_F\cdot \sqrt{2\mathit{Km}}$$
    • I_form = 2の粘性定式化の場合、摩擦力は次のとおりです：(19)
      $$F_t^{\mathit{new}}=\min \left(\mu F_n,F_{\mathit{adh}}\right)$$
    • このとき、粘着力は以下のように計算されます：(20)
      $$F_{\mathit{adh}}=F_t^{\mathit{old}}+\mathrm{\text{Δ}}{F}_t\text{with}\mathrm{\text{Δ}}{F}_t=K\cdot V_t\cdot {\delta }_t$$

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