/INTER/TYPE2

フォーマット

Spot_flag=20、21または22の場合、この入力を読み込みます：

Spot_flag=25、27または28の場合、この入力を読み込みます：

定義

コメント

1. インターフェースタイプ2は運動条件です。セカンダリサーフェスのいずれの節点上にも、他の運動条件を設定することはできません（ Spot_flag =25、27または28の場合を除く）
2. d_searchは次のように計算されます（Radioss Theory ManualのTied Interface (TYPE2)をご参照ください）：
   (1)

   $$d_{\mathit{search}}=\frac{1}{n}\cdot {\displaystyle \sum _{i=1}^n{d}_i}$$
   ここで、

   $n$

   メインセグメントの数

   $d_i$

   すべてのメイン側セグメントの合計長

3. 1つ目のインターフェースの階層レベルが、2つ目のインターフェースの階層レベルよりも低い場合に限り、1つ目のインターフェースであるインターフェースTYPE2のメイン節点を、2つ目のインターフェースTYPE2のセカンダリ節点にすることができます。階層レベルは、Spot_flag =2の場合のみ使用できます。Spot_flag =0またはSpot_flag =1の場合には機能しません。

   可能な対応策として、Spot_flag=2を使用できます。これは、/DT/NODA/CSTと適合性がないことを除き、デフォルトの定式化（Spot_flag=0）に相当します。

4. Spot_flag =2は、節点時間ステップ/DT/NODA/CSTとの適合性がないことを除き、定式化0と同等です。
5. SPH粒子をサーフェスに結合する際は、Spot_flag =4が推奨されます（粒子法流体力学（SPH）をご参照ください）。
6. Spot_flag = 20、21または22の場合は破壊を含み、接着結合のモデル化に使用されます。節点時間ステップ/DT/NODA/CSTとの適合性はありません。応力は、節点周囲の”等価”なサーフェスに従って、セカンダリ節点ごとに計算されます。
   この場合、セカンダリ節点の力は、下記のように計算される縮退力係数Fac_N (Fac_T)によってスケーリングされます：(2)

   $$\mathit{Fac}\_N=\text{min}\left\{\sqrt{\frac{{\sigma }_{N\_\text{max}}{}^2}{\text{max}\left[{\left({\sigma }_N\left(t\right)\right)}^2,{10}^{-20}\right]}},1\right\}$$

   (3)

   $$\mathit{Fac}\_T=\text{min}\left\{\sqrt{\frac{{\sigma }_{T\_\text{max}}{}^2}{\text{max}\left[{\left({\sigma }_T\left(t\right)\right)}^2,{10}^{-20}\right]}},1\right\}$$

   縮退力が最大値と比較されます:

   ${\sigma }_N<{\sigma }_{N\_\max }$ の場合は、Fac_N =1となります。これは、力が減少しないことを意味します。

   ${\sigma }_N>{\sigma }_{N\_max}$ の場合は、 $Fac\_N=\sqrt{\frac{{\sigma }_{N\_\max }^2}{\max \left[{\left({\sigma }_N\left(t\right)\right)}^2,{10}^{-20}\right]}}$ となります。これは、力が減少することを意味します。

   ここで、最大値は下記を利用してユーザーが定義します：(4)

   $${\sigma }_{N_{\max }}=\mathrm{Fscale}\left(\dot{\sigma }\right)\cdot {\mathrm{f}}_{sn}\left(\frac{\text{Δ}{X}_N}{Fscal{e}_{dist}}\right)$$

   (5)

   $${\sigma }_{T_{\max }}=\mathrm{Fscale}\left(\dot{\sigma }\right)\cdot {\mathrm{f}}_{st}\left(\frac{\text{Δ}{X}_T}{Fscal{e}_{dist}}\right)$$

   (6)

   $$\mathrm{Fscale}\left(\dot{\sigma }\right)=Fscal{e}_{stress}\cdot {\mathrm{f}}_{sr}\left(\frac{\dot{\sigma }}{Fscal{e}_{str\_rate}}\right)$$

   ここで、 ${\mathrm{f}}_{sn}$ 、 ${\mathrm{f}}_{st}$ 、 ${\mathrm{f}}_{sr}$ は、fct_ID_sn、fct_ID_st、fct_ID_srの関数です。

   破裂基準（Ruptによって定義）に達すると、接触は削除されます。

   ここで、

   • ${\sigma }_{N\_\max }$ は、fct_IDによって定義される最大法線応力値_sn
   • ${\sigma }_N\left(t\right)$ は、法線応力
   • $\sigma$ _{T_max}は、fct_IDによって定義される最大接線応力値_st
   • ${\sigma }_T\left(t\right)$ は、接線応力
   • Fscale_stressは、応力係数の入力定数
   • fct_ID_sr は、可変係数の入力値
   • fct_ID_snとfct_ID_s は、入力される応力-変位関数
   • I_symは、対称破断または非対称破断（引張り / 圧縮）のどちらかを選択できるようにします。メインサーフェスからセカンダリ節点への初期方向は、正の側（引張り）を定めます。その距離がゼロの場合（セカンダリ節点がメインサーフェス上にある）は、I_sym =1でも、破断は対称になります。

   この破壊オプション（Spot_flag = 20、21または22）は、陰解法では使用できません。

7. Spot_flag =30： メイン節点に対するセカンダリの質量 / 慣性 / 剛性の分布はKirschoffモデルをベースとし、線形関数（標準の定式化）の代わりに双3次形式の関数が使用されます。この場合、荷重 / モーメントの伝達で要素の形状曲率が考慮されるため、より柔らかな接触の挙動が可能になります。
   警告: この定式化は、回転自由度が必要となるため、ソリッド要素とは適合性がありません。
8. フラグI_del2 =2の場合、4節点シェル、3節点シェル、またはソリッド要素が削除されると、それがインターフェースのメイン側からも削除されます（関連セカンダリ節点上の運動条件が除去されます）。
9. オプションI_del2 =1とI_del2 =2オプションは、Radioss Engineでの明示的な削除（Radioss Engine入力のキーワード/DEL（/DEL/SHELL、/DEL/BRICKなど）を使用）を使用してメイン要素が削除された場合に機能します。
10. I_filtrが1に設定されている場合、垂直応力と接線応力はalphaフィルタによって、下記のようにフィルタリングされます：
    (7)

    $${\sigma }_N\left(t\right)=Alpha\cdot {\sigma }_N\left(t\right)+\left(1-Alpha\right){\sigma }_N\left(t-dt\right)$$

    (8)

    $${\sigma }_T\left(t\right)=Alpha\cdot {\sigma }_T\left(t\right)+\left(1-Alpha\right)\cdot {\sigma }_T\left(t-dt\right)$$
11. Spot_flag =25（ペナルティの定式化）は、実行時全体にわたってペナルティの定式化を保持します。（この接触の）セカンダリ節点を、剛体等の別の運動学的オプションのセカンダリ節点にすることもできます。

    ペナルティ剛性は一定で、メインおよびセカンダリ側の平均節点剛性として計算されます。これは、必要に応じて剛性係数Stfacを使用して修正することができます。ペナルティ剛性には、Stfacを掛けます。

    臨界粘性減衰係数（Visc）でインターフェース剛性に減衰を作用させることができます。

12. Ignore = 1、2、または3の場合、Starter実行中にメインセグメントが見つからなかったセカンダリ節点をインターフェースから削除します。
13. Ignore ≠ 1000の場合、d_searchが使用されます。
    Ignore = 2または3でd_search = 0の場合、d_searchはそれぞれのセカンダリ節点に対して次のように計算されます：(9)

    $${\delta }_1=0.6\left(thickness\_secondary\_node+thickness\_main\_segment\right)$$

    (10)

    $${\delta }_2=0.05\left(main\_segment\_diagonal\right)$$

    $d_{search}=\max \left({\delta }_1,{\delta }_2\right)$

    シェルの場合：

    • thickness_secondary_node = セカンダリのシェルの板厚
    • thickness_main_segment = セカンダリのシェルの板厚
    ソリッドの場合：

    • thickness_secondary_node = 0
    Ignore = 2の場合：

    • thickness_main_segment = $\frac{Element\_volume}{Segment\_area}$
    Ignore = 3の場合：

    • thickness_main_segment = 0
    Ignore = 2または3の場合：

    • 板厚が維持される順序は、まず/PART定義から、そして/SHELLまたは/SH3Nの定義から、そして/PROPの定義からとなります。
14. 接触は、Spot_flag=0の場合のみ、2次元の平面および軸対称のシミュレーションと適合性があります。ソリッド要素をSpot_flag=0で結合する場合は、モーメントは伝達されません。
15. I_del2 =1の場合に、メインセグメントに属しているすべての4節点シェル、すべての3節点シェル、およびすべてのソリッド要素が削除されると、このセグメントはインターフェースのメイン側からも削除されます（関連セカンダリ節点上の運動条件が除去されます）。
16. Spot_flag = 25、27または28： インターフェースのペナルティ剛性は、I_stfフラグに応じて、メインセグメントの剛性K_mとセカンダリ節点の剛性K_sの両方から次のように計算されます：
    • I_stf = 1： $K_n=Stfac\cdot K_m$
    • I_stf = 2（デフォルト）： $K_n=Stfac\cdot \frac{K_m+K_s}{2}$
    • I_stf = 3： $K_n=Stfac\cdot \max \left(K_m,K_s\right)$
    • I_stf = 4： $K_n=Stfac\cdot \min \left(K_m,K_s\right)$
    • I_stf = 5： $K_n=Stfac\cdot \frac{K_m\cdot K_s}{K_m+K_s}$
17. I_the >1の場合、セカンダリ側およびメインの材料は、熱伝導に有限要素定式化を使用する熱材料とする必要があります（/HEAT/MAT）。

    熱伝導は、セカンダリ節点が接触状態になった場合に計算されます。

    熱交換は、メインからセカンダリ、およびセカンダリからメインの方向で計算されます：(11)

    $${\varphi }_{cond}=K_{the}\left(T_s^{}-T_m^{}\right)$$
18. デフォルトでは、セカンダリ節点がメイン要素の外側に投影される際に正しくない質量分布を持たないよう、I_proj =1が使用されます。質量と慣性は、このエッジ上のセカンダリ節点の投影に基づき、最も近いエッジ上に分布されます。Radiossバージョン14.0またはそれ以前と同じ結果を得るには、I_proj =2を使用します。

    
    
    図 1.

19. ペナルティの定式化Spot_flag=25を使用する際、モーメントはセカンダリ節点からメインセグメントに伝えられません。したがってこれは、セカンダリ節点が回転自由度を擁する接触に使用することが推奨されます。これには、シェル-シェル、スプリング-シェル、シェルがセカンダリでソリッドがメインであるシェル-ソリッドが含まれます。この制限および運動条件定式化と比べて低いロバスト性のため、運動条件定式化とペナルティ定式化の混合Spot_flag =27および28を使用することが推奨されます。
20. Spot_flag =27および28は、運動条件定式化とペナルティ定式化の混合タイド接触です。デフォルトでは、運動条件定式化が使用されます。非適合の運動条件があるセカンダリ節点は、自動的にペナルティ定式化に切り替えられます。剛体、強制変位、強制速度、強制化速度、その他のタイド接触セカンダリ節点と不適合な運動条件または境界条件は、ペナルティ定式化への切り替えの要因となります。セカンダリ節点がペナルティ定式化に切り替えられると、Starter出力ファイルに警告のメッセージが出力されます。

    ペナルティ定式化剛性は一定で、I_stfとStfacの使用により計算されます。臨界粘性減衰係数 (Visc) でインターフェース剛性に減衰を作用させることができます。ペナルティ定式化はモーメントをセカンダリ節点からメインセグメントに伝えることができます。

21. Spot_flag =1とは異なり、タイド接触のメインサーフェスがシェル要素である際、Spot_flag =28は時間=0における質量を追加しません。メインサーフェスがソリッド要素である場合、多少の質量追加があるかもしれません。Spot_flag =27が使用される場合、質量は追加されません。