/INTER/TYPE25

ソリッド要素は、接触ギャップの板厚がゼロです。接触入力は、単一サーフェス、サーフェス対サーフェス、または節点対サーフェスを指定できます。

このインターフェースは、インターフェースTYPE3、TYPE5、TYPE7、TYPE19またはTYPE24に代えて使用できます。

このオプションは陰解法解析では使用できません。

フォーマット

定義

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

定義

コメント

1. 接触メイン / セカンダリペアは以下の3つの方法で定義できます：
   • 1つの自己接触サーフェスのみ： surf_ID₁ > 0かつsurf_ID₂ = 0
   • 対称的なサーフェス対サーフェス： surf_ID₁ > 0かつsurf_ID₂ > 0
   • 節点対サーフェス： grnd_ID_s > 0、surf_ID₁ = 0かつsurf_ID₂ > 0

   grnd_ID_s > 0は、節点対サーフェス接触タイプを定義するために使用されますが、他の接触タイプでも使用できます。その場合、節点グループは、単純に補足的セカンダリ節点として追加されます。これは、ユーザーがスプリング要素節点、剛体のメイン節点などを接触に（セカンダリ節点として）追加する場合に役立ちます。

   サーフェスがシェルを使用して定義されている場合は、反対の法線方向を持つ2つの接触セグメント（半分の板厚（t）でシフトされたもの）が生成されます：

   
   
   図 1.

   SPMDの場合、surf_IDi （i=1, 2）によって定義される各メインセグメントを1つの要素に関連付ける必要があります（ボイド要素の場合もあります）。

   2次要素が使用されている場合は、2次要素の中間節点が接触処理で使用されるため、/SURF/PART/EXTを使用してサーフェスを定義することが推奨されます。

   サーフェス定義/SURF/PART/ALLは、TYPE25では使用できません。

2. 接触剛性 $K$ は以下のように計算されます：(1)
   $$K=\max \left[S{t}_{\min },\min \left(S{t}_{\max },K_n\right)\right]$$
   ここで、 $K_n$ はI_stfに依存：

   • I_stf = 1000、 $K_n=K_m$
   • I_stf = 2、 $K_n=\frac{K_m+K_s}{2}$
   • I_stf = 3、 $K_n=\max \left(K_m,K_s\right)$
   • I_stf = 4、 $K_n=\min \left(K_m,K_s\right)$
   • I_stf = 5、 $K_n=\frac{K_m\cdot K_s}{K_m+K_s}$
   $K_m$ ： メインセグメントの剛性で、次のように計算されます：

   • $K_m=\mathit{Stfac}\cdot 0.5\cdot E\cdot t$ 、メインセグメントがシェル上にある場合。
   • $K_m=\mathit{Stfac}\cdot B\cdot \frac{S^2}{V}$ メインセグメントがソリッド上に存在する場合。
   • $K_m=\max \left(\mathit{Stfac}\cdot 0.5\cdot E\cdot t,\mathit{Stfac}\cdot B\cdot \frac{S^2}{V}\right)$ 、メインセグメントがシェルとソリッドによって共有されている場合。
   $K_s$ ：セカンダリ節点剛性は、インターフェースTYPE25として考慮される相当節点剛性で、次のように計算されます：

   • $K_m=\mathit{Stfac}\cdot 0.5\cdot E\cdot t$ 、節点がシェル要素に結合されている場合
   • $K_s=Stfac\cdot B\cdot \sqrt[3]{V}$ 、節点がソリッド要素に結合されている場合
   ここで、

   $S$

   セグメント面積

   $V$

   ソリッドの体積

   $B$

   体積弾性率

   $t$

   シェルの板厚

   Stfacの値は1.0より大きい値であることが可能です。剛性係数の値に制限はありません（値が1.0より大きいと、最初の時間ステップが短くなる場合があります）。

   /PROP/VOIDと/MAT/VOIDを使用する際、ボイド材料の材料プロパティと板厚が入力されなければなりません。そうしないと、ボイド要素の接触合成はゼロになります。シェル要素の剛性が接触計算に使用されるため、これは、ボイドシェル要素がソリッド要素と要素を共有する場合に特に重要です。

3. ギャップは衝撃のそれぞれに対して次のように自動的に計算されます：
   • I_gap = 1の場合、可変ギャップは次のように計算されます：(2)
     $$\min (g_s,Gap\_\max \_s)+\min (g_m,Gap\_\max \_m)$$
   • I_gap = 2の場合、可変ギャップは次のように計算されます：(3)
     $$\min (g_s,Gap\_\max \_s)+\min (g_m,Gap\_\max \_m)$$
     ここで、要素サイズがギャップ値よりも小さい場合、セカンダリ節点の非アクティブ化

     
     
     図 2.

     自己接触では、Curvilinear Distance（メインセグメントの節点からセカンダリ節点まで）が $\sqrt{2}\cdot Gap$ よりも小さい場合（初期構成において）、このセカンダリ節点はこのメインセグメントで考慮されず、他のメインセグメントの接触からは削除されません。

   • I_gap = 3の場合、可変ギャップは次のように計算されます：(4)
     $$\min \left[\min (g_s,Gap\_\max \_s)+\min (g_m,Gap\_\max \_m),\%mesh\_size\cdot \left(g_{s\_l}+g_{m\_l}\right)\right]$$
     ここで、

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

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

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

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

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

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

       $g_s=Gap\_scale*\frac{\sqrt{S}}{2}$ 、セカンダリ節点がトラス要素またはビーム要素に結合されている場合。ここで、 $S$ は1次元要素の断面

       フリーエッジ上のセカンダリシェル節点のギャップI_gap0が1の場合： $g_s$ は、セカンダリー節点がセカンダリサーフェスのフリーエッジ上にある場合、ゼロにリセットされます。自由エッジ上のセカンダリシェル節点のギャップ修正フラグは、セカンダリ節点がオプションの節点グループ（grnod_ID）で定義されている場合は何の影響も及ぼしません。

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

   • $gm\_l$ ： メインセグメントの最も短いエッジの長さ
   • $gs\_l$ ： セカンダリ節点がメインサーフェスに属している場合、 $gs\_l$ はそのセカンダリ節点に結合しているメインセグメントの最も短いエッジの長さ、そうでない場合は $gs\_l$ =1E+30。

     いかなる場合でも、 $g_m$ と $g_s$ は、ギャップの計算前にGap_max_mとGap_max_sによって別々に制限されます。

     セカンダリ節点がメインサーフェスに属していない場合、ギャップは下記のままとなります (5)

     $$\min (g_s,Gap\_\max \_s)+\min (g_m,Gap\_\max \_m)$$
4. 節点とサーフェスの接触では、ギャップがサーフェスの外部境界からセカンダリ節点のギャップ以上に伸びることはありません。I_shape は、このギャップの形状が正方形または円形で接触力（法線）方向であるかどうかを決定します。I_shape は、エッジ対エッジ接触の場合、ギャップとその形状には影響しません。
   I_shapeに応じて、メインサーフェスの外部境界での接触に使用されるギャップとその結果としての力の方向。

   
   
   図 3. I_shape=1（正方形ギャップ）

   
   
   図 4. I_shape=2（円形ギャップ）

   I_shape =1I_gap =3とは使用できず、I_shape =2にリセットされます。

5. シェル要素のエッジ対エッジ接触では、ギャップは丸くなります。フリーエッジのメイン側の接触ギャップは、エッジがシェルセグメントからはみ出さないようにフリーエッジをシフトさせます。

   
   
   図 5. エッジ接触メイン側

   フリーエッジにおけるセカンダリ接触ギャップは、図に示すように、I_gap0の値に依存します。

   
   
   図 6. エッジ接触セカンダリ側

6. I_edge=11、I_edge=13のときのソリッド要素では、セカンダリ側はEdge_angleよりも小さい角度の鋭角なエッジのみで構成されます。I_edge=22では、ソリッド要素からのすべてのエッジがセカンダリ側にあると考えられます。メイン側では、3つのI_edgeのケースすべてについて、ソリッド要素からのすべてのエッジが含まれています。

   
   
   図 7. I_edge=11とI_edge=13のセカンダリ側エッジ

7. fric_IDが定義されている場合、接触摩擦は/FRICTIONで定義され、この入力カード内の摩擦入力（I_fric、C₁など）使用されません。
   摩擦力は次のとおりです：(6)

   $$F_t^{new}=\min \left(\mu F_n,F_{adh}\right)$$

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

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

   ここで、 $\mu$ はクーロン摩擦で、以下のように定義されます：

   • フラグI_fricの場合（デフォルト）：

     $\mu =Fric$ ここで $F_T\le \mu \cdot F_N$ （クーロン摩擦）

   • フラグI_fric > 1の場合、新しい摩擦モデルが導入されます。この場合、摩擦係数は次の関数によって設定されます：

     $\mu =\text{μ}(\rho ,V)$

     ここで、

     $\rho$

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

     $V$

     メインセグメントに相対するセカンダリ節点の接線速度

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

   以下の定式化を使用できます：

   • I_fric = 1（汎用の粘性摩擦則）：(7)
     $$\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則）：(8)
     $$\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$ 、静的摩擦係数は ${\mu }_{\min }<{\mu }_s<{\mu }_{\max }$
     • $C_2={\mu }_d$ 、動的摩擦係数は ${\mu }_{\min }<{\mu }_d<{\mu }_{\max }$
     • $C_3={\mu }_{\max }$ 、最大摩擦係数
     • $C_4={\mu }_{\min }$ 、最小摩擦係数
     • $C_5=V_{cr1}\ne 0$ 、第1臨界速度は、> 0
     • $C_6=V_{cr2}$ 、第2臨界速度は $>V_{cr1}$

     • 第1臨界速度 $V_{cr1}=C_5$ は、第2臨界速度 $V_{cr2}=C_6\left(C_5<C_6\right)$ より小さくする必要があります。
     • 静止摩擦係数 $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$ ）。

       表 1. 摩擦定式化の単位

       Ifric	Fric	C1	C2	C3	C4	C5	C6
       1		$\left[\frac{1}{\text{P}\text{a}}\right]$	$\left[\frac{\text{s}}{\text{m}}\right]$	$\left[\frac{\text{s}}{\text{Pa}\cdot \text{m}}\right]$	$\left[\frac{1}{{\text{Pa}}^2}\right]$	$\left[\frac{{\text{s}}^2}{{\text{m}}^2}\right]$
       2		$\left[\frac{1}{{\text{Pa}}^2}\right]$	$\left[\frac{\text{s}}{\text{m}}\right]$	$\left[\frac{1}{\text{P}\text{a}}\right]$	$\left[\frac{\text{s}}{\text{m}}\right]$		$\left[\frac{\text{s}}{\text{m}}\right]$
       3						$\left[\frac{\text{m}}{\text{s}}\right]$	$\left[\frac{\text{m}}{\text{s}}\right]$

8. 摩擦フィルタリング
   I_filtr = 1、2または3の場合は、接線力がフィルタを使用して以下のようにスムージングされます：(9)

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

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

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

9. InactiおよびIpen_max、初期貫通の取り扱い：
   • Inacti = 1000： 初期貫通が無視されます。接触力は加えられませんが、節点は接触から非アクティブ化されません。節点が接触しなくなり、後で接触し直した場合は、そのときに接触力が加えられます。

     
     
     図 8.

   • Inacti = -1： 初期荷重はすべての貫通節点上に加えられます。高い初期貫通は回避する必要があります。これにより、大きな接触力が生じ、計算開始時に高エネルギーエラーとなる可能性があるからです。
   • Inacti = 5：= 5： メインセグメントは初期貫通値（ $P_0$ ）でシフトされるため、時間0においては初期荷重は付与されません。

   メインセグメントの位置が復元されるのは、 $P_0$ より大きい反発の場合のみです。

   逆に、セカンダリ節点が引き続き貫通する場合は、貫通は以下のように計算されます：(10)

   $$P\text{'}=P-P_0$$

   
   
   図 9.

   • 交差および大きい初期貫通（Inacti= -1および5）：

     シェル： 初期交差は回避する必要があります。これにより、接触力方向に誤りが生じ、セカンダリ節点の固定につながる可能性があるからです。

10. インターフェースのアクティブ化 / 非アクティブ化のためにsens_IDが定義されている場合、T_startおよびT_stopは考慮されません。
11. 出力される力について：

    接触タイプが非対称的なサーフェス対サーフェスの場合、これら2つのサーフェスが適切に分離されていれば、時刻歴内の出力法線接触力が正しく計算されます。

12. IVIS2=-1は、法線方向に接着を、接線方向に粘性抵抗力を加えるために用いられます。これは、熱塑性複合材成形のモデル化に使用できます。
    使用される際、接触ギャップの半分は接着領域、残りの半分は物理的な接触領域とみなされます。したがって、同じ物理的接触ギャップを保つには、接触板厚はGap_scaleを用いて2倍にされる必要があります。

    
    
    図 10.

    接着力は、セカンダリ節点が物理的接触領域に入り、接着領域に戻った後でのみ適用されます。接着力は、節点が接着領域からはみ出してしまうことを防ぐよう機能し、法線方向に適用されます。(11)

    $$F_N=\frac{SigMaxAdh\cdot Area}{\frac{1}{2}Gap}(\frac{1}{2}Gap-P_{adh})$$

    ここで、

    $Area$

    セカンダリサーフェスの面積

    $P_{adh}$

    接着領域への貫通

    $Gap$

    コメントにあるとおり計算された接触ギャップ 3

    接着スプリングは節点が接着領域から出ると破断し、節点が再度接着領域に入ると再生されます。

    粘性抵抗力は、セカンダリ節点が接着領域に入った際に、接線方向に適用されます。反接線方向の粘着力は摩擦力の代わりに適用され、次のように計算されます：(12)

    $$F_T=-(ViscAdhFact)\frac{ViscFluid\cdot Area}{\frac{1}{2}Gap}{V}_{rel}$$

    ここで、

    $Area$

    セカンダリサーフェスの面積

    $V_{rel}$

    接着領域への貫通

    $Gap$

    接触ギャップ

    
    
    図 11.