摩擦

Radiossではいくつかの摩擦定式化が利用可能です。最も単純なものは、最も広く用いられているもので、Coulomb摩擦則です。この定式化は衝突解析で精度のよい結果を与え、1つのパラメータ(Coulomb摩擦係数、 μ )だけを必要とします。


図 1. 節点に作用する鉛直力と接線力


図 2. 粘着力の計算
μ のデフォルト値は0(サーフェス間の摩擦なし)です。摩擦力の計算では、デフォルトの摩擦ペナルティ定式化は接線速度に基づいた粘性定式化です。スライディングの間貫通した節点は位置C0(時刻tの接触位置)からC1(時刻 t + Δ t の接触位置)に移動します。接触は粘性型であるため、粘性係数Cが粘着力の計算に導入されます:(1)
F a d h = C V t
ここで、
C = V I S F 2 K M
K
瞬間的なインターフェース剛性
VISF
インターフェース摩擦の臨界減衰係数
M
メイン節点の質量
粘着力(Fh)が計算された際、それが μ F n より小さい場合、摩擦力Fhに等しいままで変更されず、固着が起きます。粘着力が μ F n より大きい場合、摩擦力が低下し、 μ F n と等しくなります。(2)
F t = min ( μ F n , F a d h )


図 3. 摩擦力の計算
スライディングが非常に低い速度で発生すると(例えば準静的シミュレーション)、摩擦力は接線速度に基づいて計算されるため、粘性定式化は機能しません。この制限を克服するため、接線変位に基づく新しい摩擦ペナルティ定式化(剛性増分定式化)が利用可能です。この方法では人工的な剛性Kが摩擦力の変化を計算するために導入されます:(3)
Δ F t = K V t δ t
ここで、
δ t
接線変位
そのため、前の定式化とは異なり、剛性定式化では、低速であっても正しい摩擦力を計算できます。 図 4 がこの点を示しています。強制変位がパート(3D 立方体)に低速で( 0.01 m/s)作用した場合、粘性定式化はうまく働きませんが、接線変位に基づく剛性定式化ではうまくいきます。


図 4. 粘性定式化対剛性定式化
他の摩擦定式化も利用可能で、その原理はCoulomb摩擦則と同様です。Radiossでは、まず粘着力が計算され、それが次に μ F n と比較されます。その違いは摩擦係数( μ )にあり、それらは最早一定ではなく、メインセグメントの法線方向力の圧力と、セカンダリ節点の接線速度の関数となります。フラグIfricにより、 3つの新しい定式化が利用可能です:
一般化粘性摩擦則:
(4)
μ = F r i c + C 1 . p + C 2 V + C 3 . p V + C 4 p 2 + C 5 V 2
修正Darmstad摩擦則:
(5)
μ = F r i c + C 1 e ( C 2 V ) p 2 + C 3 e ( C 4 V ) p + C 5 e ( C 6 V )
Renard摩擦則:
(6)
μ = C 1 + ( C 3 C 1 ) V C 5 ( 2 V C 5 )
の、 (7)
V [ C 5 , C 6 ]
(8)
μ = C 3 ( ( C 3 C 4 ) ( V C 5 C 6 C 5 ) 2 ( 3 2 V C 5 C 6 C 5 ) )
の、 (9)
V [ C 5 , C 6 ]
(10)
μ = C 2 1 1 C 2 C 4 + ( V C 6 ) 2
の、 (11)
V C 6
(12)
P = C 1 μ + C 4 ρ C υ T = C 1 μ + α υ T
注: 全ての摩擦定式化に摩擦のフィルタリングが利用可能で摩擦力をスムーズにすることができます。詳細については、RadiossStarter入力をご参照ください。


図 5. Renard摩擦モデルのグラフィカルな表現