/ALE/GRID/STANDARD

ブロックフォーマットキーワード ALE グリッド速度計算の標準定式化を記述します。

これはエッジスプリングと反せん断スプリングに基づいた改良された/ALE/GRID/SPRING定式化です。 1

ale_standard
図 1.

フォーマット

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10)
/ALE/GRID/STANDARD
α γ η l c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGJbaabeaaaaa@37FC@  
空白のフォーマット

定義

フィールド 内容 SI単位の例
α 最大剛性に対するスケールファクター 2

デフォルト = 0.9(実数)

 
γ エッジスプリング剛性の非線形倍率 3

デフォルト = 1e-2  (実数)

 
η 減衰係数。 4

デフォルト = 1e-2  (実数)

 
l c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGJbaabeaaaaa@37FC@ 特性長さ

(実数)

[ m ]

コメント

  1. ソリッド要素に仮想のスプリングが導入されグリッド速度を制御します。
    このスプリングは非線形の粘弾性です。安定性を確実にするため、その剛性は時間ステップから計算されます。2種類のスプリングがエッジと反せん断スプリングです。
    1. エッジスプリング
      エッジスプリングの力は時間での長さ変化の関数です。(1)
      Δ F e d g e = k ( h ) ( w 2 w 1 ) d t

      ここで、

      w 1 , w 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaC4DamaaBa aaleaacaaIXaaabeaakiaacYcacaWH3bWaaSbaaSqaaiaaikdaaeqa aaaa@3A80@ はそれぞれ、節点N1N2のグリッド速度。

      h は反対の面からのN1 の距離です。

      d t は時間ステップ。

      k ( h ) はスプリング剛性 k ( h ) = k c r i t i c a l

      h が特性長さ l c よりも小さく、N1 が反対の面に達する場合は、次のようになります:(2)
      k ( h ) = 1 λ 2 [ γ + ( γ 1 ) ( h l c l c ) 3 ] k c r i t i c a l
      1 λ 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSGaaeaaca aIXaaabaGaeq4UdW2aaWbaaSqabeaacaaIYaaaaaaaaaa@3960@ は、減衰係数 β 、スケールファクター α 、および右記の時間ステップを考慮する安定化係数です; d t 4

      ale_standard_towards
      図 2.
      それ以外の場合は、 k ( h ) = γ λ 2 k c r i t i c a l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaci4Aamaabm aabaGaamiAaaGaayjkaiaawMcaaiabg2da9maalaaabaGaeq4SdCga baGaeq4UdW2aaWbaaSqabeaacaaIYaaaaaaakiaadUgadaWgaaWcba Gaam4yaiaadkhacaWGPbGaamiDaiaadMgacaWGJbGaamyyaiaadYga aeqaaaaa@4751@

      ale_standard_eq
      図 3.
    2. 反せん断スプリング

      反せん断力 F s h e a r MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGZbGaamiAaiaadwgacaWGHbGaamOCaaqabaaaaa@3B9E@ は節点の貫通から計算されます。gapは反対の面からの l c s MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSGaaeaaca WGSbWaaSbaaSqaaiaadogaaeqaaaGcbaGaam4Caaaaaaa@3910@ です。

      F s h e a r MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGZbGaamiAaiaadwgacaWGHbGaamOCaaqabaaaaa@3B9E@ の値は:

      (3)
      F s h e a r = g a p k ( h ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWGZbGaamiAaiaadwgacaWGHbGaamOCaaqabaGccqGH9aqp caWGNbGaamyyaiaadchacqGHflY1ciGGRbWaaeWaaeaacaWGObaaca GLOaGaayzkaaaaaa@4526@
      および (4)
      k ( h ) = γ λ 2 [ γ + ( γ 1 ) ( h l c l c ) 3 ] k c r i t i c a l MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaci4Aamaabm aabaGaamiAaaGaayjkaiaawMcaaiabg2da9maalaaabaGaeq4SdCga baGaeq4UdW2aaWbaaSqabeaacaaIYaaaaaaakmaadmaabaGaeq4SdC Maey4kaSYaaeWaaeaacqaHZoWzcqGHsislcaaIXaaacaGLOaGaayzk aaWaaeWaaeaadaWcaaqaaiaadIgacqGHsislcaWGSbWaaSbaaSqaai aadogaaeqaaaGcbaGaamiBamaaBaaaleaacaWGJbaabeaaaaaakiaa wIcacaGLPaaadaahaaWcbeqaaiaaiodaaaaakiaawUfacaGLDbaaca WGRbWaaSbaaSqaaiaadogacaWGYbGaamyAaiaadshacaWGPbGaam4y aiaadggacaWGSbaabeaaaaa@5929@

      ale_standard_critical
      図 4.
    3. 粘性減衰
      粘性力は剛性の上限に対応する臨界減衰から計算されます: 1 λ 2 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaSGaaeaaca aIXaaabaGaeq4UdW2aaWbaaSqabeaacaaIYaaaaaaaaaa@3961@ (5)
      F v i s c o u s = β α ( 1 + β 2 β ) ( w 2 w 1 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaCOramaaBa aaleaacaWG2bGaamyAaiaadohacaWGJbGaam4BaiaadwhacaWGZbaa beaakiabg2da9iabek7aIjabeg7aHnaabmaabaWaaOaaaeaacaaIXa Gaey4kaSIaeqOSdi2aaWbaaSqabeaacaaIYaaaaaqabaGccqGHsisl cqaHYoGyaiaawIcacaGLPaaadaqadaqaaiaahEhadaWgaaWcbaGaaG OmaaqabaGccqGHsislcaWH3bWaaSbaaSqaaiaaigdaaeqaaaGccaGL OaGaayzkaaaaaa@50A2@
    4. グリッド速度
      グリッド速度は次に以下の式に従って更新されます:(6)
      w n + 1 = w n + ( Δ F e d g e + F s h e a r + F v i s c o u s ) d t m

      ここで、 m はスプリングからの仮想の節点質量です(Starterの実行中に自動的に計算されます)。

  2. α = 1 が増加すると、最大剛性が大きくなります。スケールファクター α は、長さ0で指定されたスプリングの最大剛性を決定します。スケールファクターにより、限界安定値を超えないようにすることができます(時間ステップの減少を避けるため)。
  3. このフラグは剛性形状に作用します。剛性は γ = 0で線形です。また、 γ が増加すると、エッジスプリングの下限剛性は大きくなります。スプリングは長さ0の時の限界剛性を持ちます(これは係数1倍に相当します)。長さが特性長さ以上の場合、スプリング剛性は限界剛性値に γ をかけたものになります。
  4. β には小さな値を用いることを推奨します。そうでない場合、減衰が限界を超えるかもしれません。安定化係数は次のようになります:(7)
    λ = d t α 0 ( 1 + β 2 β ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4UdWMaey ypa0ZaaSaaaeaacaWGKbGaamiDaaqaaiabeg7aHnaaBaaaleaacaaI WaaabeaakmaabmaabaWaaOaaaeaacaaIXaGaey4kaSIaeqOSdi2aaW baaSqabeaacaaIYaaaaaqabaGccqGHsislcqaHYoGyaiaawIcacaGL Paaaaaaaaa@458A@
  5. l c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGJbaabeaaaaa@37FC@ は以下で長さを定義します:
    • エッジスプリング剛性を増加 h < l c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiAaiabgY da8iaadYgadaWgaaWcbaGaam4yaaqabaaaaa@39ED@
    • 反せん断スプリングをアクティブ化 h < l c 5 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiAaiabgY da8maaliaabaGaamiBamaaBaaaleaacaWGJbaabeaaaOqaaiaaiwda aaaaaa@3AC8@
  6. 全てのこれらのパラメータはEngineリスタートの際に修正できます(/ALE/GRID/STANDARD)。
  7. メッシュ自動修正。反せん断力の重みを次のどちらかでより大きくすることができます:
    • l c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiBamaaBa aaleaacaWGJbaabeaaaaa@37FC@ パラメータをメッシュサイズに近く設定
    • α パラメータに負の値を設定(現在の実行の最初のサイクルでエッジの弾性力が0に設定されます)
  8. この方法は節点周りでの均一なスプリングの再分配を仮定します。このれは2つのメッシュのトポロジーが異なる場合はこのケースではありません。