/INTER/TYPE7

ブロックフォーマットキーワード インターフェースTYPE7は、メインサーフェスとセカンダリ節点グループ間の接触をモデル化する多用途衝撃インターフェースです。熱伝導および熱摩擦を考慮することもできます。

内容

インターフェースTYPE3、TYPE4、およびTYPE5で課せられていたすべての制限が、このインターフェースでは解決されています。
  • 1つの節点を同時にセカンダリ節点とメイン節点にすることができます。
  • 各セカンダリ節点は、現在のセグメントに結合しているメインセグメントを除き、各メインセグメントに衝撃を与えることができます。
  • 1つの節点が複数のセグメントに衝撃を与えることができます。
  • また、1つの節点が、各セグメントの2つの側面、エッジ、およびコーナーに衝撃を与えることができます。
  • 制約なしの高速検索アルゴリズムです。
このインターフェースの主な制限を以下に示します。
  • 高衝撃速度または小規模ギャップを伴う接触の場合は、時間ステップが短縮されます。
  • 高衝撃速度での剛体または小規模ギャップを伴う剛体で使用すると、正しく機能しません。
  • エッジ間接触は解決されません(これを解決するには、/INTER/TYPE11をTYPE7と一緒に使用する必要があります)。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/INTER/TYPE7/inter_ID/unit_ID
inter_title
grnd_IDs surf_IDm Istf Ithe Igap   Ibag Idel Icurv Iadm
Fscalegap Gapmax Fpenmax     ITIED  
Stmin Stmax %mesh_size dtmin Irem_gap Irem_i2
Icurv = 1または2の場合挿入
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
node_ID1 node_ID2                
必要なフィールド
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Stfac Fric Gapmin Tstart Tstop
IBC     Inacti VISs VISF Bumult
Ifric Ifiltr Xfreq Iform sens_ID fct_IDF AscaleF fric_ID
Ifric > 0の場合
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
C1 C2 C3 C4 C5
Ifric > 1の場合
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
C6                
Iadm = 2の場合
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
NRadm Padm Angladm          
Ithe = 1の場合
(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
Kthe fct_IDK   Tint Ithe_form AscaleK  
Frad Drad Fheats Fheatm    

定義

フィールド 内容 SI単位の例
inter_ID インターフェースの識別子

(整数、最大10桁)

 
unit_ID unit_ID

(整数、最大10桁)

 
inter_title インターフェースのタイトル

(文字、最大100文字)

 
grnd_IDs セカンダリ節点グループ識別子

(整数)

 
surf_IDm メインサーフェスの識別子

(整数)

 
Istf インターフェース剛性定義フラグ 3
SPHでは、Istf=01および1000のみ可能。
= 0
/DEFAULT/INTER/TYPE7で定義した値に設定されます。
= 1
インターフェース剛性は、Stfacとして入力されます。
= 2
インターフェース剛性は、メインおよびセカンダリ剛性の平均です。
= 3
インターフェース剛性は、メインおよびセカンダリ剛性の最大値です。
= 4
インターフェース剛性は、メインおよびセカンダリ剛性の最小値です。
= 5
インターフェース剛性は、直列のメインおよびセカンダリ剛性です。
= 1000 /DEFAULT/INTER/TYPE7が定義されていない場合のデフォルト
インターフェース剛性はメイン側の剛性のみに基づきます。

(整数)

 
Ithe 熱接触フラグ。
=0
熱伝導も熱摩擦もなし
=1
熱伝導または熱摩擦がアクティブ

(整数)

 
Igap ギャップ/要素オプションフラグ 11
= 0
/DEFAULT/INTER/TYPE7で定義した値に設定されます。
= 1
可変ギャップは被衝撃メインサーフェスおよび衝撃セカンダリ節点の特性に応じて変化します。
= 2
計算するギャップの可変ギャップ + ギャップスケール補正
= 3
可変ギャップ+算出されたギャップのギャップスケール修正+自己接触中の初期貫通回避のために考慮されるメッシュサイズ
= 1000 /DEFAULT/INTER/TYPE7が定義されていない場合のデフォルト
一定のギャップ; 最小ギャップGapminと一致

(整数)

 
Ibag 接触時のエアバッグベントホール閉鎖フラグ
このフラグは、対応するソリッド要素が接触している場合にSol2SPH粒子をアクティブにするためにも使用されます。
= 0
/DEFAULT/INTER/TYPE7で定義した値に設定されます。
= 1
閉鎖
= 2 /DEFAULT/INTER/TYPE7が定義されていない場合のデフォルト
閉鎖なし

(整数)

 
Idel 節点およびセグメントの削除フラグ。
= 0
/DEFAULT/INTER/TYPE7で定義した値に設定されます。
= 1
1つのセグメントに関連付けられたすべての要素4節点シェル、3節点シェル、ソリッド)が削除されると、メイン側の面からそのセグメントを削除します。このセグメントはEngineファイルでRadioss Engineキーワード/DELを使用した明示的な削除の場合にも、削除されます。
さらに、結合されていない節点をインターフェースのセカンダリ側から削除します。
= 2
4節点シェル、3節点シェル、またはソリッド要素が削除されると、メイン側の面から、対応するセグメントを削除します。このセグメントはEngineファイルでRadioss Engineキーワード/DELを使用した明示的な削除の場合にも、削除されます。
さらに、結合されていない節点をインターフェースのセカンダリ側から削除します。
= 1000 /DEFAULT/INTER/TYPE7が定義されていない場合のデフォルト
削除しません。
= -1
結合されていない節点をインターフェースのセカンダリ側から削除する以外は、=1と同様です。
= -2
結合されていない節点をインターフェースのセカンダリ側から削除する以外は、=2と同様です。

(整数)

注: Idel=1および-1は、Idel=2および-2と比較してCPUコストが高くなります。
 
Icurv 曲率を伴うギャップエンベロープ 5
= 0(デフォルト)
曲率なし。
=1
球面。
=2
円筒面。
=3
自動双3次曲面。

(整数)

 
Iadm アダプティブメッシュでの局所曲率計算フラグ。 6 7 8
= 0(デフォルト)
非アクティブ。
=1
メッシュサイズに応じたインターフェース更新。
=2
メッシュサイズ、貫通、および角度に応じたインターフェース更新。

(整数)

 
Fscalegap ギャップスケールファクター(Igap = 2および3の場合にのみ使用)

デフォルト = 1.0(実数)

 
Gapmax 最大ギャップ(Igap = 2および3の場合にのみ使用)

(実数)

[ m ]
Fpenmax 初期貫通の最大比率。 13

(実数)

 
ITIED Tie接触節点フラグ。
= 0(デフォルト)
時間ステップのため節点同士は結合されていません。
=1
接触時間ステップが、接触で定義された/DT/INTER/DELまたはdtmin未満の場合、セカンダリ節点は結合され、リバウンド可能です。
=2
接触時間ステップが、接触で定義された/DT/INTER/DELまたはdtmin未満の場合、セカンダリ節点は結合され、リバウンドはしません。

(整数)

 
Stmin 最小剛性(Istf = 234または5の場合にのみ使用されます)

(実数)

[ N m ]
Stmax 最大剛性(Istf = 234または5の場合にのみ使用されます)

デフォルト = 1030(実数)

[ N m ]
%mesh_size メッシュサイズのパーセンテージ(Igap = 3の場合にのみ使用)

デフォルト = 0.4(実数)

 
dtmin 最小インターフェース時間ステップ 23

(実数)

[ s ]
Irem_gap 自己接触において要素サイズがギャップ値より小さい場合にセカンダリ節点を非アクティブ化するためのフラグ。 14
= 0
/DEFAULT/INTER/TYPE7で定義した値に設定されます。
= 1 /DEFAULT/INTER/TYPE7が定義されていない場合のデフォルト
セカンダリ節点の非アクティブ化なし
= 2
セカンダリ節点の非アクティブ化あり

(整数)

 
Irem_i2 インターフェースTYPE2で同じ接触ペア(節点)が定義されている場合に、セカンダリ節点を非アクティブにするためのフラグ。
= 0
/DEFAULT/INTER/TYPE7で定義した値に設定されます。
= 1
/INTER/TYPE2タイド接触内のセカンダリ節点は、この接触から削除されます。
= 3 /DEFAULT/INTER/TYPE7が定義されていない場合のデフォルト
セカンダリ節点に変更なし。
 
node_ID1 第1節点識別子

(整数)

 
node_ID2 第2節点識別子(Icurv = 1の場合は無視されます)

(整数)

 
Stfac インターフェースのメイン側に適用される剛性スケールファクター(Istf 1の場合)

デフォルト = 1.0(実数)

 
インターフェース剛性(Istf = 1の場合)

デフォルト = 1.0(実数)

[ N m ]
Fric Coulomb摩擦(fct_IDF = 0の場合)

デフォルト = 0.0(実数)

 
クーロン摩擦のスケールファクター(fct_IDF 0の場合) 20

デフォルト = 1.0(実数)

Gapmin 衝撃アクティブ化の最小ギャップ

デフォルト値については、10を参照のこと(実数)

[ m ]
Tstart 開始時間。

(実数)

[ s ]
Tstop 一時的な非アクティブ化の時間。

(実数)

[ s ]
IBC 接触時の境界条件の非アクティブ化フラグ

(ブーリアン)

 
Inacti 初期貫通時の剛性の非アクティブ化フラグ 13
= 0
/DEFAULT/INTER/TYPE7で定義した値に設定されます。
= 1
節点の剛性を非アクティブ化します。
= 2
要素の剛性を非アクティブ化します。
= 3
初期貫通が発生しないように、節点の座標を変更します。
= 5
ギャップは時間とともに変化し、初期ギャップは次のように調整されます:
gap 0 = Gap P 0 ここで P 0 は初期貫通
= 6
ギャップは時間とともに変化しますが、初期ギャップは次のように調整されます(節点の貫通が少し解除されます):
gap 0 = Gap P 0 5 % ( Gap P 0 )
= 1000 /DEFAULT/INTER/TYPE7が定義されていない場合のデフォルト
操作は行われません。

(整数)

 
VISs インターフェース剛性の臨界減衰係数

デフォルトは0.05に設定されます(実数)

 
VISF 接触インターフェース摩擦の臨界減衰係数。 22

デフォルトは1.0に設定されます(実数)

 
Bumult ソート係数はソートアルゴリズムの高速化に使用されます。また、マシンによって異なります。

デフォルトは0.20に設定されます(実数)

 
Ifric 摩擦定式化フラグ。 1617
fric_IDが定義されていない場合のみ使用されます。
= 0(デフォルト)
静的なCoulomb摩擦則。
= 1
汎用の粘性摩擦則。
= 2
(修正)Darmstad摩擦則。
= 3
Renard摩擦則。
= 4
指数関数的減衰摩擦則。

(整数)

 
Ifiltr 摩擦フィルタリングフラグ。 18
= 0(デフォルト)
フィルターを使用しません。
= 1
単純な数値フィルター。
= 2
フィルタリング時間による標準の-3dBフィルター。
= 3
カットオフ周波数による標準の-3dBフィルター。

(整数)

 
Xfreq フィルタリング係数。

デフォルト値はIfiltrに依存します(実数)

 
Iform 摩擦ペナルティ定式化のタイプ。 19
= 0
/DEFAULT/INTER/TYPE7で定義した値に設定されます。
= 1 /DEFAULT/INTER/TYPE7が定義されていない場合のデフォルト
粘性(合計)の定式化。
= 2
剛性(増分)の定式化。

(整数)

 
sens_ID インターフェースをアクティブ化 / 非アクティブ化するためのセンサーの識別子 24

識別子センサーが定義されている場合は、インターフェースのアクティブ化 / 非アクティブ化はセンサーに基づき、TstartまたはTstopには基づきません。

(整数)

 
fct_IDF 温度による摩擦係数の関数識別子

デフォルト = 0(整数)

 
AscaleF fct_IDFの横軸のスケールファクター

デフォルト = 1.0(実数)

 
fric_ID 選択されたパートのペアについての摩擦定義の摩擦識別子
= 0(デフォルト)
このインターフェース内で定義されている摩擦パラメータを使用します。
0
/FRICTION/fric_IDを使用します。

(整数)

 
C1 - C6 摩擦則係数。

(実数)

参照: 表 1
Kthe 熱交換係数(fct_IDK = 0の場合)

デフォルト = 0.0

[ W m 2 K ]
熱交換係数(fct_IDK 0の場合) 20

デフォルト = 1.0(実数)

fct_IDK 接触圧力に対する熱交換定義の関数識別子

デフォルト = 0(整数)

 
AscaleK fct_IDKの横軸のスケールファクター

デフォルト = 1.0(実数)

[ Pa ]
Tint インターフェース温度 20

(実数)

[ K ]
Ithe_form 熱接触定式化フラグ
= 0(デフォルト)
インターフェース(一定温度)とシェル(セカンダリ側)の間のみで交換
= 1
すべての接触片間の熱交換

(整数)

 
NRadm 90度の半径を経由する要素の数

(整数)

 
Padm 貫通の割合に関する基準

デフォルト = 1.0(実数)

 
Angladm 角度基準

(実数)

[ deg ]
Frad 放射係数

(実数)

[ W m 2 K 4 ]
Drad 放射計算の最大距離

(実数)

[ m ]
Fheats セカンダリの摩擦加熱ファクター 21

(実数)

 
Fheatm セカンダリの摩擦加熱ファクター 21

(実数)

 

境界条件の非アクティブ化フラグ:IBC

(1)-1 (1)-2 (1)-3 (1)-4 (1)-5 (1)-6 (1)-7 (1)-8
          IBCX IBCY IBCZ

定義

フィールド 内容 SI単位の例
IBCX
=1
衝撃時のX境界条件の非アクティブ化フラグ

(ブーリアン)

 
IBCY
=1
衝撃時のY境界条件の非アクティブ化フラグ

(ブーリアン)

 
IBCZ
=1
衝撃時のZ境界条件の非アクティブ化フラグ

(ブーリアン)

 

コメント

  1. SPMDの場合、surf_IDmによって定義される各メインセグメントを1つの要素に関連付ける必要があります(ボイド要素の場合もあります)。
  2. フラグIbagについては、モニター体積オプション(モニター体積(エアバッグ))をご参照ください。
  3. 接触剛性Kは以下のように計算されます:
    Istf =1の場合:(1)
    K = S t f a c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbGaey ypa0Jaam4uaiaadshacaWGMbGaamyyaiaadogaaaa@3CBE@
    Istf = 234または5の場合:(2)
    K = max [ S t min , min ( S t max , K n ) ]
    Istf =1000の場合:(3)
    K = K m

    ここで、 K n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaaaa@384D@ は、メインセグメントの剛性 K m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaaaa@384D@ とセカンダリ節点の剛性 K s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaaaa@384D@ の両方から計算されます:

    Istf = 2 K n = K m + K s 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaOGaeyypa0ZaaSaaaeaacaWGlbWaaSbaaSqa aiaad2gaaeqaaOGaey4kaSIaam4samaaBaaaleaacaWGZbaabeaaaO qaaiaaikdaaaaaaa@3F01@

    Istf = 3 K n = max ( K m , K s ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaOGaeyypa0JaciyBaiaacggacaGG4bWaaeWa aeaacaWGlbWaaSbaaSqaaiaad2gaaeqaaOGaaiilaiaadUeadaWgaa WcbaGaam4CaaqabaaakiaawIcacaGLPaaaaaa@4260@

    Istf = 4 K n = min ( K m , K s ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaOGaeyypa0JaciyBaiaacMgacaGGUbWaaeWa aeaacaWGlbWaaSbaaSqaaiaad2gaaeqaaOGaaiilaiaadUeadaWgaa WcbaGaam4CaaqabaaakiaawIcacaGLPaaaaaa@425E@

    Istf = 5 K n = K m K s K m + K s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaOGaeyypa0ZaaSaaaeaacaWGlbWaaSbaaSqa aiaad2gaaeqaaOGaeyyXICTaai4samaaBaaaleaacaGGZbaabeaaaO qaaiaadUeadaWgaaWcbaGaamyBaaqabaGccqGHRaWkcaWGlbWaaSba aSqaaiaadohaaeqaaaaaaaa@4479@

    K m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaaaa@384D@ はメインセグメントの剛性で、次のように計算されます:

    メインセグメントがシェル上に存在するか、シェルとソリッドによって共有されている場合:(4)
    K m = Stfac 0.5 E t
    メインセグメントがソリッド上に存在する場合:(5)
    K m = Stfac B S 2 V
    ここで、
    S MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGtbaaaa@39AF@
    セグメント面積
    V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGtbaaaa@39AF@
    ソリッドの体積
    B MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGtbaaaa@39AF@
    体積弾性率

    K s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbWaaS baaSqaaiaad6gaaeqaaaaa@384D@ はインターフェースTYPE7として考慮される相当節点剛性で、次のように計算されます:

    節点がシェル要素に結合されている場合:(6)
    K s = 1 2 E t
    節点がソリッド要素に結合されている場合:(7)
    K s = B V 3

    剛性スケールファクターの値に対する制限はありません(ただし、値が1.0より大きいと、初期時間ステップが短くなる場合があります)。

    /PROP/VOID/MAT/VOIDを使用する際、ボイド材料の材料プロパティと板厚が入力されなければなりません。そうしないと、ボイド要素の接触合成はゼロになります。シェル要素の剛性が接触計算に使用されるため、これは、ボイドシェル要素がソリッド要素と要素を共有する場合に特に重要です。

  4. Istf = 234または5は、SPH定式化と適合性がありません。
  5. Icurv = 1の場合は、node_ID1(球の中心)を使用してギャップの球形曲率が定義されます。

    Icurv = 2の場合は(円筒形の軸上にある)node_ID1node_ID2を使用してギャップの円筒曲率が定義されます。

    Icurv = 3の場合は、メインサーフェスの形状が、座標の連続性とセグメント間の法線に配慮した双三次補間を使用して取得されます。曲率の変化が速くて大きい場合は、この定式化が不安定になる可能性があります(将来のバージョンで改善される予定です)。

    inter_type7_Icurv
    図 1.
  6. アダプティブメッシングとIadm = 1の場合:
    曲率半径が要素サイズ(セカンダリ側)未満のゾーン(メイン側)で接触が発生した場合は、セカンダリ側の要素が分割されます(最大レベルに到達していない場合)。

    inter_type7_Iadm
    図 2.
  7. アダプティブメッシングとIadm = 2の場合:

    曲率半径が要素サイズのNRadm倍(セカンダリ側)未満のゾーン(メイン側)で接触が発生した場合は、セカンダリ側の要素が分割されます(最大レベルに到達していない場合)。

    法線間の角度がAngladmを超えており、貫通の割合がPadmを超えているゾーン(メイン側)で接触が発生した場合は、セカンダリ側の要素が分割されます(最大レベルに到達していない場合)。


    図 3.
  8. 係数NRadmPadmおよびAngladmは、アダプティブメッシングとIadm=2の場合にのみ使用されます。
  9. Gapmax=0の場合、ギャップの最大値が存在しません。
  10. Gapmin=0または空白の場合、デフォルト値は以下のように計算されます:

    メインセグメントがシェルおよびソリッド要素である場合、Gapmin = min ( t m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiDamaaBe aaleaacaWGTbaabeaaaaa@380E@ , l min 2 )。

    ここで、
    t m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiDamaaBe aaleaacaWGTbaabeaaaaa@380E@
    メインシェル要素の平均板厚(Igap=0の場合)
    t m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiDamaaBe aaleaacaWGTbaabeaaaaa@380E@
    メインシェル要素の最小板厚(Igap=12または3の場合)
    l min MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGSbWaaS baaSqaaiGac2gacaGGPbGaaiOBaaqabaaaaa@3A4D@
    すべてのメインセグメント(シェルまたは3次元ソリッド)の中の最小辺長

    メインセグメントがすべてソリッド要素である場合、Gapmin = l min 10

    ここで、 l min MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGSbWaaS baaSqaaiGac2gacaGGPbGaaiOBaaqabaaaaa@3A4D@ は、すべてのメイン3次元ソリッドセグメントの中で最小の辺です。

  11. 可変ギャップ
    Igap = 1の場合、可変ギャップは次のように計算されます:(8)
    max [ Gap min , ( g s + g m ) ]
    Igap = 2の場合、可変ギャップは次のように計算されます:(9)
    max { Gap min , min [ Fscal e g a p ( g s + g m ) , Gap max ] }
    Igap = 3の場合、可変ギャップは自己接触に対して次のように計算されます:(10)
    max { Gap min , min [ Fscal e gap ( g s + g m ) , % m e s h _ s i z e ( g s _ l + g m _ l ) , Gap max ] }
    ここで、
    • g m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaad2gaaeqaaaaa@3868@ : メイン要素のギャップ

      g m = t 2 : ここで t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG0baaaa@3758@ は、シェル要素に対するメイン要素の板厚

      g m = 0 : 3次元ソリッド要素の場合

    • g s MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaadohaaeqaaaaa@386E@ : セカンダリ節点のギャップ

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

      g s = t 2 : ここで t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG0baaaa@3758@ は、セカンダリ節点に結合されているシェル要素の最大板厚です。

      g s = 1 2 S : トラスおよびビーム要素の場合。ここで、 S MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG0baaaa@3758@ はトラスおよびビーム要素の断面

    • g m _ l MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaad2gacaGGFbGaamiBaaqabaaaaa@3A3C@ : 要素の短いエッジの長さ
    • g s _ l MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGNbWaaS baaSqaaiaad2gacaGGFbGaamiBaaqabaaaaa@3A3C@ : セカンダリ要素に結合された要素の短いエッジ長さ

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

    可変ギャップは常にGapmin以上です。

  12. 境界条件の非アクティブ化は、セカンダリ節点グループ(grnd_IDs)に適用されます。
  13. 節点がスプリング要素に属している場合にInacti = 3であると、初期エネルギーが発生することがあります。
    インターフェースに対する高周波数の影響を回避するためには、Inacti =5ではなく、Inacti = 6が推奨されます。


    図 4.

    Fpenmaxが0でない場合、次の条件が満たされると、節点剛性は非アクティブになります:

    P e n e t r a t i o n F p e n m a x G a p これはInactiの値に左右されません。

  14. Irem_gap = 2を使用すると、要素サイズをギャップ値よりも小さくできます:

    inter_type7_master_seg
    図 5. 節点対サーフェス接触から消去されたセカンダリ節点

    自己接触では、Curvilinear Distance(メインセグメントの節点からセカンダリ節点まで)が 2 G a p よりも小さい場合(初期構成において)、このスセカンダリ節点はこのメインセグメントで考慮されず、他のメインセグメントの接触からは削除されません。

  15. 1つの節点が2つのサーフェスに同時に属することができます。
  16. fric_IDが定義されている場合、接触摩擦は/FRICTIONで定義され、この入力カード内の摩擦入力(IfricC1など)使用されません。

    摩擦定式化の場合:

    摩擦フラグIfric=0(デフォルト)の場合は、従来の静的摩擦係数の定式化が使用されます:

    F t μ F n ここで、 μ (は、クーロン摩擦係数)
    • fct_IDF = 0の場合:
      IfricはCoulomb摩擦。(11)
      μ = Fric
    • fct_IDF0の場合:
      Fricは、温度に依存するクーロン摩擦係数のスケールファクターになります。(12)
      μ = Fric f F ( Ascale F , T interface )
    ここで、 T interface は、セカンダリとメインの平均温度として取得される温度:(13)
    T interface = T secondary + T main 2

    フラグIfric > 0の場合、新しい摩擦モデルが導入されます。この場合、摩擦係数は関数によって設定されます。 μ = μ ( ρ , V )

    ここで、
    ρ
    メインセグメントの垂直抗力の圧力
    V MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGwbaaaa@373A@
    メインセグメントに相対するセカンダリ節点の接線速度
  17. 現在は、係数C1C6を使用して、新しい摩擦定式化の可変摩擦係数 μ を定義しています。
    以下の定式化を使用できます:
    • Ifric = 1(汎用の粘性摩擦則):(14)
      μ = Fric + C 1 . p + C 2 V + C 3 . p V + C 4 p 2 + C 5 V 2
    • Ifric = 2(修正Darmstad則):(15)
      μ = Fric + C 1 . e ( C 2 V ) . p 2 + C 3 . e ( C 4 V ) . p + C 5 . e ( C 6 V )
    • Ifric = 3(Renard則):(16)
      μ = C 1 + ( C 3 C 1 ) V C 5 ( 2 V C 5 ) if V [ 0 , C 5 ]
      (17)
      μ = C 3 ( ( C 3 C 4 ) ( V C 5 C 6 C 5 ) 2 ( 3 2 V C 5 C 6 C 5 ) ) if V [ C 5 C 6 ]
      (18)
      μ = C 2 1 1 C 2 C 4 + ( V C 6 ) 2 if V C 6
      ここで、Ifric=3の場合:
      • C 1 = μ s 、静摩擦係数、右記である必要があります; μ min < μ s < μ max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH8oqBda WgaaWcbaGaciyBaiaacMgacaGGUbaabeaakiabgYda8iabeY7aTnaa BaaaleaacaWGZbaabeaakiabgYda8iabeY7aTnaaBaaaleaaciGGTb GaaiyyaiaacIhaaeqaaaaa@44BF@
      • C 2 = μ d 、動摩擦係数、右記である必要があります; μ min < μ d < μ max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH8oqBda WgaaWcbaGaciyBaiaacMgacaGGUbaabeaakiabgYda8iabeY7aTnaa BaaaleaacaWGKbaabeaakiabgYda8iabeY7aTnaaBaaaleaaciGGTb GaaiyyaiaacIhaaeqaaaaa@44B0@
      • C 3 = μ max 、最大摩擦係数
      • C 4 = μ min 、最小摩擦係数
      • C 5 = V cr 1 、第1臨界速度、> 0である必要があります
      • C 6 = V c r 2 、第2臨界速度、右記である必要があります; > V c r 1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqGH+aGpca WGwbWaaSbaaSqaaiaadogacaWGYbGaaGymaaqabaaaaa@3B08@
    • Ifric = 4(指数関数的減衰摩擦則):
      摩擦係数は、接触面の相対速度 V MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOvaaaa@36D1@ に依存すると仮定し、次式に従う:(19)
      μ = C 1 + F r i c C 1 e C 2 V MathType@MTEF@5@5@+= feaahqart1ev3aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH8oqBcq GH9aqpcaWGdbWaaSbaaSqaaiaaigdaaeqaaOGaey4kaSYaaeWaaeaa caWGgbGaamOCaiaadMgacaWGJbGaeyOeI0Iaam4qamaaBaaaleaaca aIXaaabeaaaOGaayjkaiaawMcaaiabgwSixlaadwgadaahaaWcbeqa amaabmaabaGaeyOeI0Iaam4qamaaBaaameaacaaIYaaabeaalmaaem aabaGaamOvaaGaay5bSlaawIa7aaGaayjkaiaawMcaaaaaaaa@4F0A@
    表 1. 摩擦定式化の単位
    Ifric Fric C1 C2 C3 C4 C5 C6
    1   [ 1 P a ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaaGymaaqaaiaaccfacaGGHbaaaaGaay5waiaaw2faaaaa @3AD5@ [ s m ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaae4Caaqaaiaab2gaaaaacaGLBbGaayzxaaaaaa@3A46@ [ s Pa m ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaae4CaaqaaiaabcfacaqGHbGaeyyXICTaaeyBaaaaaiaa wUfacaGLDbaaaaa@3E47@ [ 1 Pa 2 ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaaGymaaqaaiaabcfacaqGHbWaaWbaaSqabeaacaaIYaaa aaaaaOGaay5waiaaw2faaaaa@3BC6@ [ s 2 m 2 ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaae4CamaaCaaaleqabaGaaGOmaaaaaOqaaiaab2gadaah aaWcbeqaaiaaikdaaaaaaaGccaGLBbGaayzxaaaaaa@3C2C@  
    2   [ 1 Pa 2 ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaaGymaaqaaiaabcfacaqGHbWaaWbaaSqabeaacaaIYaaa aaaaaOGaay5waiaaw2faaaaa@3BC6@ [ s m ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaae4Caaqaaiaab2gaaaaacaGLBbGaayzxaaaaaa@3A46@ [ 1 P a ] MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaaGymaaqaaiaaccfacaGGHbaaaaGaay5waiaaw2faaaaa @3AD5@ [ s m ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaae4Caaqaaiaab2gaaaaacaGLBbGaayzxaaaaaa@3A46@   [ s m ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaae4Caaqaaiaab2gaaaaacaGLBbGaayzxaaaaaa@3A46@
    3           [ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@ [ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
    4     [ s m ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWadaqaam aalaaabaGaae4Caaqaaiaab2gaaaaacaGLBbGaayzxaaaaaa@3A46@        
  18. 摩擦フィルタリング:
    Ifiltrフラグ 0の場合は、接線力がフィルタを使用して以下のようにスムージングされます:(20)
    F T f = α F T ( t ) + ( 1 α ) F T f ( t d t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGubGaamOzaaqabaGccqGH9aqpcqaHXoqycaWHgbWaaSba aSqaaiaadsfaaeqaaOGaaiikaiaadshacaGGPaGaey4kaSYaaeWaae aacaaIXaGaeyOeI0IaeqySdegacaGLOaGaayzkaaGaaCOramaaBaaa leaacaWGubGaamOzaaqabaGccaGGOaGaamiDaiabgkHiTiaadsgaca WG0bGaaiykaaaa@4D30@
    ここで、
    F T f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGubGaamOzaaqabaaaaa@38B5@
    フィルタリングされた接線力。
    F T ( t ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGubaabeaakiaacIcacaWG0bGaaiykaaaa@3A26@
    フィルター前の時間tにおける計算された接線力。
    F Tf (tdt) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGubGaamOzaaqabaGccaGGOaGaamiDaiabgkHiTiaadsga caWG0bGaaiykaaaa@3DE0@
    前の時間ステップでフィルタリングされた接線力
    t MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizaiaads haaaa@37D8@
    現在のシミュレーション時間
    dt MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizaiaads haaaa@37D8@
    現在のシミュレーション時間ステップ
    α MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3796@
    フィルタリング係数
    ここで、 α MathType@MTEF@5@5@+= feaahqart1ev3aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@3792@ 係数は、以下のように計算されます:
    • Ifiltr = 1の場合: α = X f r e q 、単純な数値フィルターで、値の範囲は0~1です。
    • Ifiltr = 2の場合: α = 2 π X f r e q 、標準の-3dBフィルターで、フィルター処理する時間ステップ数は X freq = d t T として定義されます。 T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGubaaaa@39B0@ はフィルター処理周期です。
    • Ifiltr = 3の場合: α = 2 π X freq dt 標準の-3dBフィルターで、Xfreq = カット周波数。
  19. 摩擦ペナルティ定式化Iform
    • Iform = 1(デフォルト)の粘性定式化の場合、摩擦力は次のとおりです:(21)
      F t = min ( μ F n , F adh )
      このとき、粘着力は以下のように計算されます:(22)
      F adh = C V t with C = VIS F 2 Km
    • Iform = 2の粘性定式化の場合、摩擦力は次のとおりです:(23)
      F t new = min ( μ F n , F adh )

      このとき、粘着は以下のように計算されます:

      F a d h = F t o l d + Δ F t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGgbWaaS baaSqaaiaadggacaWGKbGaamiAaaqabaGccqGH9aqpcaWGgbWaa0ba aSqaaiaadshaaeaacaWGVbGaamiBaiaadsgaaaGccqGHRaWkcaqGuo GaamOramaaBaaaleaacaWG0baabeaaaaa@43D7@ ここで、 Δ F t = K V t d t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaqGuoGaam OramaaBaaaleaacaWG0baabeaakiabg2da9iaadUeacqGHflY1caWG wbWaaSbaaSqaaiaadshaaeqaaOGaeyyXICTaamizaiaadshaaaa@43C9@

      V t は、メインセグメントに相対するセカンダリ節点の接線速度

      陰解析と陽解析には、Iform = 2が推奨されています。

  20. 熱交換:
    Ithe=1(熱伝導をアクティブ化)にすることで、接触の熱交換と熱摩擦を考慮します。
    • Ithe_form = 0の場合、熱交換はシェルと一定温度接触Tintの間で行われます。
    • Ithe_form = 1の場合、熱交換はすべての接触片間で行われます。

      TintIthe_form=0の場合にのみ使用されます。この場合、メイン側の温度は一定(Tintに等しい)と想定されます。Ithe_form=1の場合、Tintは考慮されません。このため、メイン側の節点温度が考慮されます。

    熱交換係数:
    • fct_IDK = 0の場合、Ktheは熱交換係数で、熱交換は熱交換サーフェスのみに依存します。
    • fct_IDK0の場合、Ktheはスケールファクターで、熱交換は接触圧力に依存します:(24)
      K = K the f K ( Ascale K , P )
    • ここで、 f K fct_IDKの関数です。
  21. 熱摩擦:
    • インターフェースType 7では、Ithe > 0である場合、摩擦エネルギーは熱に変換されます。
    • FheatsFheatmは摩擦エネルギーの比率で、エネルギーがこの比率でそれぞれセカンダリ側とメイン側に分配されます。したがって一般的に:(25)
      Fheat s + Fheat m 1.0

      FheatsFheatmがいずれも0であると、摩擦スライディングエネルギーから熱への変換はアクティブになりません。

    • 摩擦熱QFricは、以下のように定義されます:
      • Iform= 2の場合(剛性定式化):
        セカンダリ側: (26)
        Q Fric = Fheat s ( F adh F t ) K F t
        メイン側: (27)
        Q Fric = Fheat m ( F adh F t ) K F t
        Ithe_form=1
      • Iform= 1の場合(ペナルティ定式化):
        セカンダリ側: (28)
        Q Fric = Fheat s C V t 2 d t
        メイン側: (29)
        Q Fric = Fheat m C V t 2 d t
        Ithe_form=1
  22. 放射:
    放射が接触で考慮されるのは、 F r a d 0 で、セカンダリ節点のメインセグメントまでの距離、 d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGebWaaS baaSqaaiaadkhacaWGHbGaamizaaqabaaaaa@3A1A@ が次の場合です:(30)
    Gap < d < D rad
    ここで、 D r a d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGebWaaS baaSqaaiaadkhacaWGHbGaamizaaqabaaaaa@3A1A@ は放射計算の最大距離です。 D r a d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGebWaaS baaSqaaiaadkhacaWGHbGaamizaaqabaaaaa@3A1A@ のデフォルト値は、以下の最大値として計算されます:
    • 全節点中のGapの上限値(時間0における)
    • セカンダリ要素の最小辺長

    D r a d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9 Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGebWaaS baaSqaaiaadkhacaWGHbGaamizaaqabaaaaa@3A1A@ の値は、Radioss Engineのパフォーマンスの減少を引き起こすかもしれないため、高すぎる値に設定しないことを推奨します。

    熱伝導の放射は以下のように計算されます:(31)
    h rad = F rad ( T m 2 + T s 2 ) ( T m + T s )
    ここで、(32)
    F rad = σ 1 ε 1 + 1 ε 2 1
    ここで、
    σ = 5.669 × 10 8 [ W m 2 K 4 ]
    シュテファンボルツマン定数
    ε 1
    セカンダリサーフェスの輻射率
    ε 2
    メインサーフェスの輻射率
  23. この接触内のセカンダリ節点の時間ステップがdtminより小さくなると、セカンダリ節点は接触から消去され、警告のメッセージが出力ファイルに書き出されます。このdtmin値は、/DT/INTER/DEL内に入力されるあらゆるモデルインターフェースの最小タイムステップより優先されます。
  24. インターフェースのアクティブ化 / 非アクティブ化のためにsens_IDが定義されている場合、TstartおよびTstopは考慮されません。
  25. SOL2SPH定式化(/PROP/TYPE14(SOLID))を使用するパートでは、オプションIdel= -1または-2をアクティブにする必要があります。