/PROP/TYPE45 (KJOINT2)

フォーマット

定義

ジョイントタイプリスト

例（回転）

#RADIOSS STARTER
/UNIT/2
unit for prop
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/TYPE45/2/2
Revolute 
### - Define angle limit < 0.52rad 
### - And define friction moment 100GPa to block angle (if it reached the limit)
#     Type                  KN                 SCF                  CR  SENSORID
         2                   0                   0                   0         0
#                KR1Func_ID_Kr                SA1-                SA1+  Icomb_r1
                   0         0                   0                 .52         0
#                CR1Func_ID_Cr
                   0         0
#               KFR1                 FM1   FCT_FM1
                   0                 100         0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA

例（並進）

#RADIOSS STARTER
/UNIT/2
unit for prop
                  kg                  mm                  ms
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/TYPE45/12/2
Translational 
### - define displacement limit: -100mm ~ 100 mm
#     Type                  KN                 SCF                  CR  SENSORID
         6                   0                   0                  .2         0
#                KX1Func_ID_Kx                SD1-                SD1+  Icomb_t1
                   0         0                -100                 100         0
#                CX1Func_ID_Cx
                   0         0
#               KFX1                 FF1   FCT_FF1
                1000                   0         0
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA

コメント

1. KJOINT2スプリングは通常、剛体を接続するために使用されますが、変形可能な要素の節点を接続するためにも使用できます。変形可能な要素の節点を接続する際、要素の節点に回転自由度（DOF）がない場合は、モーメントは伝達されません。これには以下のケースが含まれます：ソリッド要素とトラス要素の節点に対するすべての回転自由度、シェルプロパティI_dril=2（デフォルト）を持つシェル要素のドリリング自由度、およびビーム要素とスプリング要素の節点に対する何らかの自由度。
2. スプリング要素：
   • その識別子は、それぞれの要素ファミリー内で一意である必要がありますが、グローバルモデルにおいて、各要素タイプが固有の要素識別子を持つことが推奨されます。
   • パートの定義には複数のスプリングブロックを使用できます。
3. スプリングの自由度：
   • ジョイントプロパティはジョイント要素の局所座標系で定義されます。
   • 局所座標系のフレームで計算されるジョイントの自由度数は６です：

     ${\delta }_{{\mathbf{X}}^{\prime }},{\delta }_{{\mathbf{Y}}^{\prime }},{\delta }_{{\mathbf{Z}}^{\prime }},{\theta }_{{\mathbf{X}}^{\prime }},{\theta }_{{\mathbf{Y}}^{\prime }},{\theta }_{{\mathbf{Z}}^{\prime }}$

   • 固定自由度とフリー自由度はそれぞれのジョイントタイプで区別されます。
   • 固定自由度には剛性が一定しているという特徴があります。デフォルトでは、剛性値が入力されていない場合、剛性の値が時間ステップを保つように自動計算されます。
   • 並進および回転自由度は次のように定義されます：(1)
     $$\delta =d{x}_2-d{x}_1$$
     ここで、 $d{x}_1$ と $d{x}_2$ は局所座標系における2つのジョイント節点の全変位です。(2)

     $$\theta ={\theta }_2-{\theta }_1$$

     ここで、 ${\theta }_1$ と ${\theta }_2$ は、ジョイントの局所座標系を基準とした、2つの結合された物体軸の全相対回転です。

   • sens_IDが定義された場合、ジョイントはセンサーがアクティブになった時（全ての自由度が）完全に固定されます。
4. 力とモーメントの計算：
   • 方向の力 $\mathbf{\delta }$ は次のように計算されます：
     線形スプリング：(3)

     $$\mathbf{F}={\mathbf{K}}_{ti}\mathbf{\delta }+{\mathbf{C}}_{ti}\dot{\mathbf{\delta }}$$

     $K_{ti}$ ： 並進剛性 $\left(K_{tx},K_{ty},K_{tz}\right)$

     $C_{ti}$ ： 並進粘性 $\left(C_{tx},C_{ty},C_{tz}\right)$

     非線形スプリング：(4)

     $$F=K_{ti}\mathrm{f}\left(\delta \right)+C_{ti}\mathrm{g}\left(\dot{\delta }\right)$$
   • $\mathbf{\theta }$ 方向のモーメントは次のように計算されます：
     線形スプリング：(5)

     $$\mathbf{M}={\mathbf{K}}_{\mathit{ri}}\mathbf{\theta }+{\mathbf{C}}_{\mathit{ri}}\dot{\mathbf{\theta }}$$

     $K_{ri}$ ： 回転剛性（ $\left(K_{rx},K_{ry},K_{rz}\right)$ ）

     $C_{ri}$ ： 回転粘性（ $\left(C_{rx},C_{ry},C_{rz}\right)$ ）

     非線形スプリング：(6)

     $$\mathbf{M}={\mathbf{K}}_{\mathit{ri}}\mathrm{f}\left(\mathbf{\theta }\right)+{\mathbf{C}}_{\mathit{ri}}\mathrm{g}\left(\dot{\mathbf{\theta }}\right)$$
   • ジョイントの長さは0にすることもできますが、常に0に等しくなければならないということはありません。球面ジョイントやユニバーサルジョイントの定義には長さ0のスプリングを用いることを推奨します。
   • 一般的なケースでモーメントの全体的なバランスを保つため、回転自由度の補正項は次のように計算されます：(7)
     $${\mathbf{M}}_{\theta x}={\mathbf{M}}_{\theta x}+{\mathbf{L}}_y\times {\mathbf{F}}_z-{\mathbf{L}}_z\times {\mathbf{F}}_y$$
     (8)
     $${\mathbf{M}}_{\theta y}={\mathbf{M}}_{\theta y}+{\mathbf{L}}_z\times {\mathbf{F}}_x-{\mathbf{L}}_x\times {\mathbf{F}}_z$$
     (9)
     $${\mathbf{M}}_{\theta z}={\mathbf{M}}_{\theta z}+{\mathbf{L}}_x\times {\mathbf{F}}_y-{\mathbf{L}}_y\times {\mathbf{F}}_x$$

     ジョイントはユーザー定義の質量や慣性をもつことができません、このため節点時間ステップが常に用いられます。節点の時間ステップは、ジョイントが接続するエンティティに基づいています。剛体のメイン節点に取り付けられている場合、剛体の時間ステップが使用されます。変形可能要素の節点に取り付けられている場合、節点時間は変形可能な要素に基づいて計算されます。

5. スプリングの剛性
   • ユーザー定義の関数がない場合、係数 $K_{ti}$ と $K_{ri}$ は一定の剛性として用いられます。関数が定義されている場合、対応する剛性係数は関数のスケールファクターになります。
   • K_n = 0の場合、固定剛性は時間ステップを保つように併進と回転自由度に対して計算開始時に自動的に計算されます。これらの値はモデルの物理現象に応じて選択することもできますが、その剛性は隣接要素よりも高い必要があります。
   • K_n = 0の場合、Scfは並進と回転固定剛性の両方に作用するスケールファクターです。K_n > 0の場合、Scfは固定回転剛性にのみ作用します。このパラメータは固定回転剛性をマニュアルで調整するために用いられます。
6. スプリングの粘性
   • ユーザー定義の関数がない場合、係数 $C_{ti}$ と $C_{ri}$ は線形粘性係数として用いられます。関数が定義されている場合、対応する係数は関数のスケールファクターになります。
   • 粘性減衰を組み込む方法は以下の2つです:
     • 臨界減衰の定義（固定自由度の場合のみ）：

       臨界減衰係数の項で粘性減衰を定義します。臨界減衰係数は、要素の固定剛性値と換算質量（接続された節点の質量の調和平均の半分）を使用して計算されます。すべての固定自由度には同じ臨界減衰が適用されます。

     • ユーザー定義の一定減衰または非線形減衰：

       フリー自由度毎に独立した減衰パラメータを定義できます。

     剛性および減衰の詳細については、ユーザーズガイドのスプリングジョイントTYPE45（/PROP/KJOINT2）をご参照ください。

7. 摩擦
   • K_ftiまたはK_friが定義されていない場合、摩擦はアクティブになりません。
   • FFiとFMiはユーザー定義の関数がない場合に、一定の摩擦力とモーメントとして用いられます。摩擦関数番号が0でない場合、FM_iとFF_iは関数のスケールファクターとなります（デフォルト値 = 1.0）。
8. スプリングの限界
   • SD_i-またはSD_i+に0以外の値が指定された場合、変位が正の変位のためのSD_i+を超えることや、負の変位のためのSD_i-を超えることがないよう、さらにペナルティ力がかけられます。このペナルティ力はK_ftiを用いて計算されます。
   • SA_i-またはSA_i+に0以外の値が指定された場合、回転が正の回転のためのSA_i+を超えることや、負の回転のためのSA_i-を超えることがないよう、さらにペナルティモーメントがかけられます。このペナルティ力はKを用いて計算されます_fri
9. スプリングのスキュー
   • Skew_ID₁が定義されている場合、初期局所座標系はSkew_ID₁によって定義されます。Skew_ID₁ = 0である場合、局所座標系はスプリングの追加節点に従って計算されます。詳細については、/SPRINGをご参照ください。
   • Skew_ID₂が0ではない場合、ジョイントの回転角は、Skew_ID₁とSkew_ID₂の間の回転に従って初期化されます。初期角度の値は、Starter出力ファイルで確認できます。回転ジョイント（TYPE2）、円筒ジョイント（TYPE3）、および平面ジョイント（TYPE4）では ${\theta }_X$ のみを初期化できます。Skew_ID₁とSkew_ID₂の最初の軸は平行であることが必要です。

     
     
     図 1. 回転ジョイント（タイプ2）の初期回転角の計算

   • Skew_ID₂が0ではなく、Skew_ID₁ =0である場合、Skew_ID₁は全体座標系です。
   • 局所座標系の向きは、node_ID₁（/SPRING）の回転に伴って更新されます。
10. 停止変位 / 角度の結合のフラグ：
    • Icomb_t_i = 0:の場合： 停止変位は独立（デフォルト）です（図 2に示すとおり）。停止変位は、各方向で別々にチェックされます。
    • Icomb_t_i = 1の場合： すべての自由な並進自由度の停止変位が結合されます。停止基準は変位に適用されなくなりますが、結合される変位のノルムについては、図 3をご参照ください。
    • 2つまたは3つの停止変位が結合されます。
    • 停止変位が結合されると、同じSDi-とSDi+の値が各自由度に使用される必要があります。
    • 上記の停止変位の記述は、Icomb_r_iを用いて定義された停止角度にも適用されます。

      
      
      図 2. 独立停止変位を伴う平面ジョイント、Icomb_t_i=0

      
      
      図 3. 2つの結合された停止変位を伴う平面ジョイント、Icomb_t_i=1