PCONTX

フォーマット

例

定義

コメント

1. プロパティ識別番号は、既存のPCONTバルクデータエントリの番号である必要があります。特定のPCONTに関連付けることができるのは、１つのPCONTXプロパティ拡張のみです。
2. PCONTXは、ANALYSIS = EXPDYNで定義される幾何学的非線形解析サブケースでのみ適用されます。他のすべてのサブケースでは無視されます。
3. 幾何学的非線形解析の/INTERエントリでは、PCONTX/PCNTX#エントリでFRICが明示的に定義されていない場合は、CONTACTまたはPCONTエントリのMU1値がFRICに対して使用されます。定義されている場合は、PCONTX/PCNTX#のFRICにより、CONTACT/PCONTのMU1値が上書きされます。
4. 陰解法解析では、セカンダリとメインのセットが重なったり（自己接触）、重ならなかったりする接触に対して異なる接触定式化が使用されます。

   自己接触の場合、ギャップのゼロはあり得ないため、一定のギャップが使用されます。小さな初期ギャップの場合、の方が初期ギャップよりも大きければ、収束はより安定し、収束までの時間も短縮されます。小さな初期ギャップの場合、GAPの方が初期ギャップよりも大きければ、収束はより安定し、収束までの時間も短縮されます。

   陰解法解析では、スケーリングファクターが縮小される（STFAC = 0.01など）剛性や衝撃を受ける厚さが縮小される（剛体の場合）剛性を使用すると、力の不均衡を低減し、収束が向上する場合があります。これは特に、実効剛性が膜剛性よりはるかに小さい構造や、曲げられているシェル構造などで顕著です。ただし、値が小さ過ぎても発散につながる可能性があるので注意してください。

5. ISTF ≠ 1の場合、インターフェース剛性Kはメインセグメント剛性Kmやセカンダリセグメント剛性Ksから算出されます。

   メイン剛性は、Km = STFAC * B * S * S/V（ソリッドの場合）およびKm = 0.5 * STFAC * E * t（シェルの場合）によって算出されます。

   セカンダリ剛性は、Ks = STFAC * B * V-3（ソリッドの場合）とKs = 0.5 * STFAC * E * t（シェルの場合）から算出される等価節点剛性です。

   これらの等式において、Bは体積弾性率、Sはセグメント面積、そしてVはソリッド体積です。剛性係数の値に対する制限はありません（ただし、値が1.0より大きいと、初期時間ステップが短くなる場合があります）。

   インターフェース剛性は、K = max (STMIN, min (STMAX, K1))と次の式で求められます。

   • ISTF = 0の場合、K1 = Km
   • ISTF = 2の場合、K1 = 0.5 * (Km + Ks)
   • ISTF = 3の場合、K1 = max (Km, Ks)
   • ISTF = 4の場合、K1 = min (Km, Ks)
   • ISTF = 5の場合、K1 = Km * Ks / (Km + Ks)
6. 一定のギャップ（IGAP = CONST）のデフォルトは、下記の中で最も小さい値になります。
   • t（メインシェル要素の平均板厚）
   • l/10（lはメインソリッド要素の辺の平均長）
   • lmin/2（lminはすべてのメインセグメント（シェルまたはソリッド）の最も短い辺の長さ）
7. 可変ギャップ（IGAP = VAR）は、gs + gmで計算します。
   ここで：

   • gm - 下記のようなメイン要素

     gm = t/2（t: メイン要素の厚さ）（シェル要素の場合）

     gm = 0（ソリッド要素の場合）

   • gs - 下記のようなセカンダリ節点のギャップ：

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

     gs = t/2（tは、セカンダリ節点に結合されているシェル要素の最大厚）

     gs = 1/2✓S（トラスおよびビーム要素の場合）。Sは要素の断面です。

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

8. INACTI= 3、4の使用は、初期貫通が小規模な場合にのみ推奨します。また、下記のような理由から、その使用には注意が必要です。
   • 変更した座標を元に戻すことはできません。
   • 接触面に複数のサーフェス層を定義すると、他の初期貫通が発生する可能性があります。
   • 節点がスプリング要素に属している場合、初期エネルギーが発生することがあります。
   INACTI = 5の場合は、下記のようになります。

   
   
   図 1.

9. IFRICは摩擦モデルを定義します。

   IFRIC = COUL - F_T ≤ FRIC * F_Nを伴うCoulomb摩擦。

   IFRIC > 0の場合、摩擦係数は関数（ $\mu$ = $\mu$ (p, V)）によって設定されます。ここで、pはメインセグメント上の垂直抗力の圧力、Vはセカンダリ節点の接線速度です。

   下記の式を使用できます。

   • IFRIC = 1 - 汎用の粘性摩擦則(1)
     $$\text{m}=\text{FRIC}+\text{C}1*\text{p}+\text{C}2*\text{V}+\text{C}3*\text{p}*\text{v}+\text{C}4*{\text{p}}^2+\text{C}5*{\text{v}}^2$$
   • IFRIC = 2 - Darmstad則(2)
     $$\text{m}={\text{C}}_1\cdot {\text{e}}^{(\text{C}{\hspace{0.17em}}_2\text{V})}\cdot {\text{p}}^2+{\text{C}}_3\cdot {\text{e}}^{(\text{C}{\hspace{0.17em}}_4\text{V})}\cdot \text{p}+{\text{C}}_5\cdot {\text{e}}^{(\text{C}{\hspace{0.17em}}_6\text{V})}$$
   • IFRIC = 3 - Renard則

     $\mu =C_1+\left(C_3-C_1\right)\cdot \frac{V}{C_5}\cdot \left(2-\frac{V}{C_5}\right)$	0 ≤ V ≤ C5
     $\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)$	C5 ≤ V ≤ C6
     $\mu =C_2-\frac{1}{\frac{1}{C_2-C_4}+{\left(V-C_6\right)}^2}$	C6 ≤ V

     ここで：(3)

     $$\begin{array}{l}{C}_1=C1={\mu }_s,C_2=C2={\mu }_d\hfill \\ C_3=C3={\mu }_{\text{max}},C_4=C4={\mu }_{\text{min}}\hfill \\ C_5=C5=V_{cr1},C_6=C6=V_{cr2}\hfill \end{array}$$
   • 最初の臨界速度のV_cr1は0以外にする必要があります（C5 ≠ 0）。また、2番目の臨界速度V_cr2未満にする必要もあります（C5 < C6）。
   • 静止摩擦係数のC1と動摩擦係数のC2は、最大摩擦のC3より小さくする必要があります（C1 ≤ C3かつC2 ≤ C3）。
   • 最小摩擦係数C4は静摩擦係数のC1と動摩擦係数のC2より小さくする必要があります（C4 ≤ C1かつC4 ≤ C2）。
10. IFORMは、2種類の接触摩擦ペナルティ定式化を選択します。
    粘性（合計）の定式化（IFORM = VISC）では、接着力は下記のように算出されます。(4)

    $$\begin{array}{l}{\text{F}}_{\text{adh}}=\text{VISF}*\surd (2\text{KM})*{\text{V}}_{\text{T}}\\ {\text{F}}_{\text{T}}=\text{min}({\mu \text{F}}_{\text{N}},{\text{F}}_{\text{adh}})\\ \end{array}$$
    剛性（増分）の定式化（IFORM = STIFF）では、接着力は下記のように算出されます。(5)

    $$\begin{array}{l}{\text{F}}_{\text{adh}}={\text{F}}_{\text{Told}}+\Delta {\text{F}}_{\text{T}}\hfill \\ \Delta {\text{F}}_{\text{T}}=\text{K}*{\text{V}}_{\text{T}}*\text{dt}\hfill \\ {\text{F}}_{\text{Tnew}}=\text{min}({\mu \text{F}}_{\text{N}},{\text{F}}_{\text{adh}})\hfill \end{array}$$
11. IFILTRは、摩擦フィルタリング係数を算出する方法を定義します。IFILTR ≠ NOの場合、接線摩擦力はフィルターを使用して下記のようにスムージングされます。(6)
    $${\text{F}}_{\text{T}}=\alpha *{\text{F}\prime }_{\text{T}}+(1-\alpha )*{\text{F}\prime }_{\text{T}-1}$$

    ここで、

    F_T は接線力

    F'_Tは時間tにおける接線力

    F'_T-1は時間t-1における接線力

    αはフィルタリング係数

    IFILTR = SIMP conref="../../bank/solvers_symbols_equations_r_b.dita#reference_xh2_vkk_pw/o2_sym_alpha" /> = FFAC

    IFILTR = PER conref="../../bank/solvers_symbols_equations_r_b.dita#reference_xh2_vkk_pw/o2_sym_alpha" /> = 2πdt/FFAC, ここで、dt/T = FFAC, Tは、フィルタリング時間

    IFILTR = CUTF conref="../../bank/solvers_symbols_equations_r_b.dita#reference_xh2_vkk_pw/o2_sym_alpha" /> = 2π * FFAC * dt、ここで、FFACはカットオフ周波数

12. このカードは、HyperMesh内のPCONTプロパティの拡張として表現されます。