時間ステップ

陽解法は、小さな時間増分または時間ステップで結果を計算することによって解かれます。時間ステップのサイズは多くの要因に依存しますが、Radiossによって自動的に計算されます。

結果は各時間ステップ毎、もしくはシミュレーション内のサイクルについて計算されます。したがって、時間ステップが小さいほどシミュレーションが解を得るまでの時間は長くかかります。これは、より多くのサイクルと計算が行われるためです。Radioss Theory ManualDynamic Analysisで説明しているとおり、運動方程式を解くために、直接積分法が使用されます。Radiossで使用される直接積分法は、Newmark時間積分スキームから得られます。この手法は、数値的に安定した時間ステップ Δ t MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0baaaa@3B37@ を使ってステップ毎の手順で運動方程式を解きます。 Numerical Stability of Undamped SystemsRadioss Theory Manual内)は、 Δ t 2 ω max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacqqHuoarcaWG0bGaeyizIm6aaSaaaeaacaaIYaaabaGaeqyYdC3a aSbaaSqaaiGac2gacaGGHbGaaiiEaaqabaaaaaaa@4285@ である場合、減衰なしの系が安定していることを示します。ここで、 ω max は系の最大角振動数です。離散有限要素シミュレーションでは、メッシュを通過する衝撃波が1つの時間ステップ中に1つ以上の要素を通過しない場合、解は安定を保ちます。この方法では、衝撃波はメッシュを通過する際にいかなる節点も失わず、したがって、有限要素メッシュ内のすべての振動数を励起します。材料 c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbwvMCKf MBHbqefqvATv2CG4uz3bIuV1wyUbqedmvETj2BSbqefm0B1jxALjhi ov2DaebbnrfifHhDYfgasaacH8srps0lbbf9q8WrFfeuY=Hhbbf9v8 qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq=He9 q8qqQ8frFve9Fve9Ff0dmeaacaGacmGadaWaaiqacaabaiaafaaake aacaWGJbaaaa@39BF@ の音速および有限要素の特性要素長 l c を使うことで、1つの要素長を通過する波動の時間は次のようになります:(1)
Δ t = l c c MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iDaiabg2da9maalaaabaGaamiBamaaBaaaleaacaWGJbaabeaaaOqa aiaadogaaaaaaa@3C63@
離散解を安定に保つには、時間ステップを、波動が1つの要素を通過するために必要な時間以下にする必要があります:(2)
Δt l c c MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaam iDaiabgsMiJoaalaaabaGaamiBamaaBaaaleaacaWGJbaabeaaaOqa aiaadogaaaaaaa@3D11@

この安定性の基準は、1928年のCourantによって初めて行われた研究に因んでクーラン条件と呼ばれています。 1