/INTER/TYPE18

ブロックフォーマットのキーワード このインターフェースは、FSIアプリケーションについてALEとLagrangeの間の非直接連成を管理します。ペナルティ法が使用され、剛性とギャップ値を必要とします。

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/INTER/TYPE18/inter_ID
inter_title
  surf_ID grbric_ID       Ipres Idel   Istf
Stfval Vref Gap Tstart Tstop

定義

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

(整数、最大10桁)

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

(文字、最大100文字)

 
surf_ID Lagrangeサーフェスの識別子

(整数)

 
grbric_ID Euler 3次元ソリッド要素またはALE 3次元ソリッド要素のグループの識別子。3D解析に使用。

(整数)

 
Ipres 圧力修正フラグ。
= 0
圧力は修正されません。
= 1
圧力の中央値分減圧されます。
= 2
平均サーフェス圧力分減圧されます。
= 3
平均予想体積圧力分減圧されます。

(整数)

 
Idel 節点およびセグメントの削除フラグ。
= 0(デフォルト)
削除しません。
= 1
1つのセグメントに関連付けられたすべての要素(4節点シェル、3節点シェル、ソリッド)を削除すると、そのセグメントがインターフェースのLagrange側から除去されます。
また、結合していない節点がインターフェースのEuler側またはALE側から除去されます。
= 2
4節点シェル要素、3節点シェル要素、またはソリッド要素を削除すると、インターフェースのLagrange側からその要素に対応するセグメントが除去されます。
また、結合していない節点がインターフェースのEuler側またはALE側から除去されます。

(整数)

 
Istf インターフェース剛性定義フラグ
= 1(デフォルト)
計算された剛性値とStfvalを入力する必要があります。
= 2
インターフェースの剛性値は、Vrelとスケールファクターから自動的に計算されます。

(整数)

 
Stfval インターフェースの剛性値。 2 3 4
Istf = 0または1の場合
インターフェースの剛性値。
デフォルト = 0.0
Istf=2の場合
インターフェース剛性スケール値。
デフォルト = 1.0

(実数)

[ N m ] またはなし
Vref ALE材料とLagrangeサーフェスの間の推定相対速度(Istf=2の場合にのみ使用)。 4

(実数)

[ m s ] MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaamWaaeaada Wcaaqaaiaab2gaaeaacaqGZbaaaaGaay5waiaaw2faaaaa@39DE@
Gap インターフェースのギャップ 5
= 0
一定のグローバルギャップが自動的に計算されます。

(実数)

[ m ]
Tstart 開始時間

(実数)

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

(実数)

[ s ]

コメント

  1. Vrefは基準速度であり、ALE材料とLagrangeサーフェスの間の推定相対速度と見なすことができます。シミュレーション時に予想される最大値を使用することが推奨されます。例えば、着水の適用例では、構造物の初期衝突速度を使用できます。IEDまたは鉱山発破の適用例では、構造に影響を与える衝撃波の最大材料速度を使用できます。
  2. Istf=1の場合、Stfvalの推奨値は次のとおりです:(1)
    Stfval= ρ V ref 2 S el Gap MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaiaads hacaWGMbGaamODaiaadggacaWGSbGaeyypa0ZaaSaaaeaacqaHbpGC cqGHflY1caWGwbWaaSbaaSqaaiaadkhacaWGLbGaamOzaaqabaGcda ahaaWcbeqaaiaaikdaaaGccqGHflY1caWGtbWaaSbaaSqaaiaadwga caWGSbaabeaaaOqaaiaadEeacaWGHbGaamiCaaaaaaa@4D4E@
    ここで、
    ρ
    (最も高い)流体密度。
    V r e f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOvamaaBa aaleaacaWGYbGaamyzaiaadAgaaeqaaaaa@39C9@
    相対速度。
    S e l
    Lagrange要素の表面積平均。
  3. Istf=2の場合、インターフェースの剛性値は、VrelとスケールファクターStfvalから次のように自動的に計算されます:(2)
    S t f = S t f v a l . ρ V r e f 2 S e l G a p MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4uaiaads hacaWGMbGaeyypa0Jaam4uaiaadshacaWGMbGaamODaiaadggacaWG SbGaaiOlamaalaaabaGaeqyWdiNaeyyXICTaamOvamaaBaaaleaaca WGYbGaamyzaiaadAgaaeqaaOWaaWbaaSqabeaacaaIYaaaaOGaeyyX ICTaam4uamaaBaaaleaacaWGLbGaamiBaaqabaaakeaacaWGhbGaam yyaiaadchaaaaaaa@50BC@
  4. インターフェース剛性のあるモデルの感度スタディを実施することが推奨されます。そのために、Istf=2を用いると、剛性スケールファクターを使用して、シミュレーションごとに異なる値に設定できます(例えば、値を1から2、4、10、および100に変更できます)。
  5. Gapは、インターフェースのギャップ。推奨値は、1.5 x 接触部分に対する法線方向に沿った流体要素のサイズ( 1.5 Δ x f l u i d MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaGymaiaac6 cacaaI1aGaaGPaVlabfs5aejaadIhadaWgaaWcbaGaamOzaiaadYga caWG1bGaamyAaiaadsgaaeqaaaaa@40E9@ )です。ギャップが定義されていない場合は、以下を使用して計算されます:(3)
    g a p = 3 2 d i a g max MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4zaiaadg gacaWGWbGaeyypa0ZaaSaaaeaadaGcaaqaaiaaiodaaSqabaaakeaa caaIYaaaaiaadsgacaWGPbGaamyyaiaadEgadaWgaaWcbaGaciyBai aacggacaGG4baabeaaaaa@421A@

    ここで、 d i a g max MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamizaiaadM gacaWGHbGaam4zamaaBaaaleaaciGGTbGaaiyyaiaacIhaaeqaaaaa @3C9F@ は流体要素の最大対角線であり、これは通常のメッシュの 1.5Δ x fluid MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaGymaiaac6 cacaaI1aGaaGPaVlabfs5aejaadIhadaWgaaWcbaGaamOzaiaadYga caWG1bGaamyAaiaadsgaaeqaaaaa@40E9@ とまったく同じです。

  6. 2020より前のバージョンでは、流体のEuler節点やALE節点は、1行目の第1フィールドでgrnd_IDsを使用して定義していました。この入力フォーマットは、レガシーモデルには今もサポートされています。この入力フォーマットは、LAW5、6または51のレガシーモデルでは引き続きサポートされます。
  7. ALE側のメッシュサイズを小さくすると、インターフェースの精度が向上します。したがって、インターフェースの精度を上げるために、接触ゾーンで、少なくともLagrangeサーフェスのグローバル法線方向でALEメッシュを細分化することが推奨されます。
  8. 精度を上げるために、ギャップと剛性の値が異なる複数のTYPE18インターフェースを、異なる接触ゾーンに使用できます(例えば、メッシュサイズや相対速度が大幅に変化する場合)。
  9. Lagrangeパートで材料破壊の場合の要素削除が有効になっている場合は、Idel=1または2を使用することを強く推奨します。
  10. SPMD並列化の場合は、surf_IDで定義しているLagrangeセグメントをそれぞれ1つの要素に関連付ける必要があります(この要素には/PROP/VOIDとなるプロパティが割り当てられます)。