PCNTX24

フォーマット

例

定義

コメント

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

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

   陰解法解析では、スケーリングファクターが縮小される（<varname>STFAC</varname> = 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. ギャップは、衝撃ごとにgs + gmで自動的に計算されます（PCNTX7でのIGAP = VARに似ています）。
   ここで：

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

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

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

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

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

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

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

   gmおよびgsは、ギャップ計算の前にGAPMAXmおよびGAPMAXsによって別個に制限されます。

7. 係数 C1 conref="./../bank/solvers_symbols_equations_r_b.dita#reference_xh2_vkk_pw/sym_c_6" id="pcntx24_bulk_r_ph_rtr_tsg_x1b"/> は可変摩擦係数 $\mu$ を決定するのに使用されます。
8. IFRICは摩擦モデルを定義します。

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

   IFRIC > 0の場合、摩擦係数は関数（ $\mu =\mu (p,V)$ ）によって設定されます。

   ここで、pはメインセグメント上の垂直抗力の圧力、Vはセカンダリ節点の接線速度です。

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

   • IFRIC = GEN - 汎用の粘性摩擦則(1)
     $$\mu =\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 = DARM - Darmstad則(2)
     $$\mu ={\text{C}}_1*\text{e}({\text{C}}_2\text{V})*{\text{p}}^2+{\text{C}}_3*\text{e}({\text{C}}_4\text{V})*\text{p}+{\text{C}}_5*\text{e}({\text{C}}_6\text{V})$$
   • IFRIC = REN - 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_{\mathit{cr}1},C_6=C6=V_{\mathit{cr}2}\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）。
9. IFILTRは、摩擦フィルタリング係数を算出する方法を定義します。IFILTR ≠ NOの場合、接線摩擦力はフィルターを使用して下記のようにスムージングされます。

   F_T = α * F'_T + (1 conref="../../bank/solvers_symbols_equations_r_b.dita#reference_xh2_vkk_pw/o2_sym_alpha"/>) * F'_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はカットオフ周波数

10. IFORMは、2種類の接触摩擦ペナルティ定式化を選択します。

    粘性（合計）の定式化（IFORM = VISC）では、接着力は下記のように算出されます。

    F_adh = VISF * ✓(2KM) * V_T

    F_T = min ( $\mu$ F_N, F_adh)

    剛性（増分）の定式化（IFORM = STIFF）では、接着力は下記のように算出されます。

    F_adh = F_Told + $\text{Δ}$ F_T

    $\text{Δ}$ F_T = K * V_T * dt

    F_Tnew = min ( $\mu$ F_N, F_adh)

11. インターフェースのアクティブ化 / 非アクティブ化のためにSENSIDが定義されている場合は、TSTARTおよびTSTOPは考慮されません。
12. 接触タイプが対称サーフェス間であり、２つのサーフェスが適切に分離されている場合は、THファイルの出力垂直接触力は正しく計算されます。
13. 陰解法テストの場合： インターフェースTYPE24は、現在はSMPでのみ使用可能です。ISTFのデフォルトは、4に設定されます。デフォルトのINACTIは、-1に設定されます。
14. このカードは、HyperMesh内のPCONTプロパティの拡張として表現されます。