結果チェック

エネルギーエラー計算。

出力ファイル（Runname_0001.out）のエネルギーエラーはどのように計算されていますか？(1)
% E r r o r = 100 ( E k + E k r + E i E k , 1 + E k , 1 r + E i , 1 + E w k − E w k , 1 − 1 )
ここで、
E k
現時点の並進運動エネルギー
E k r
現時点の回転運動エネルギー
E i
内部エネルギー
E w k
外力の仕事量（システムに取り込まれたエネルギー）
E , 1
実行の開始時（t=0時点のではなく）のエネルギー

このエネルギーバランスではアワグラスエネルギーは考慮されず、通常は負のエネルギーエラーが発生します（アワグラスが存在しない、Ishell=24（QEPH）またはIshell=12（BATOZ）シェルと完全積分ソリッドを使用している場合を除く）。

%Error MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiyjaiaadw eacaWGYbGaamOCaiaad+gacaWGYbaaaa@3B43@ は± 99%の範囲に制限されます。

エネルギーエラーはRESTARTの後に毎回リセットされます。

ジョブを停止し、/CHKPTを含むコントロールファイルを使用して再実行することも可能です。この場合、エラーおよびエネルギー値は、前回の値から再開されます。

エネルギーエラーに関して

エネルギーエラーはどのように理解すべきですか？またどの値が妥当と見なされますか？

Radiossにより計算されるエネルギーエラーはパーセンテージです。
  • エラーが負数の場合、一部のエネルギーが消えたことを意味します。

    低減積分要素（Belytschkoシェル、積分点が１つのソリッド）の場合、エネルギーバランスで考慮されないため、アワグラスエネルギーで負のエネルギーエラーも説明できます。アワグラスエネルギーの通常の量は-10～-15%です。

  • エラーが正数の場合、エネルギーは生成されません。

    QEPHシェル定式化（Ishell=24）または完全積分要素を使用している場合、アワグラスエネルギーは存在せず、計算はより正確であるため、エネルギーエラーはわずかに正になります。+1%または+2%のエラーは許容可能です。

    正のエネルギー誤差が2%より大きい場合、このエネルギーのソースを識別する必要があります。互換性のない運動条件によりこのような状況が発生することがあります。

  • エネルギー誤差が増加し、 99%に達する場合、これはシミュレーションは発散していることを示します。ただし、特定の状況では、高いエネルギー誤差を許容できる場合があります。
    • 系の初期エネルギーが低い場合、シミュレーションの初期には大きなエネルギー誤差がある可能性があります。この誤差は、系のエネルギーが増加すると小さくなります。これは、エネルギーにおける小さな数値での差分が、大きなパーセンテージのエネルギー誤差を生じさせるためです。
    • 一部分での発散によりシミュレーションの最後で大きなエネルギー誤差が生じる場合もありますが、その他のモデルでは適切な結果が得られます。
    • 総エネルギーに対して接触エネルギーが大きいと、負のエネルギー誤差が大きくなる可能性があります。これは、接触エネルギーがエネルギー誤差方程式に含まれていないためです。シミュレーションに摩擦および多くのスライディング接触が含まれる場合、大きな接触エネルギーおよび結果として生じるエネルギー誤差は許容可能と見なすことができます。

付加質量の確認

付加質量を考慮することは必要ですか？また、どのようにして結果が大きく変更されていないことを確認できますか？

付加質量は、インターフェースTYPE2、Spotflag =1が設定されていることが原因で発生することがあります。この場合、付加質量はすべて時間t=0に生成されています。また、一定タイムステップのオプション（/DT/NODA/CSTまたは/DT/INTER/CST）による場合もあります。

モデル内の付加質量については、モデルの総質量に対して重要ではないかどうかをチェックする必要があります（Radioss Engine出力ファイル（Runname_nnnn.out）の最後の列のDM/M値を参照）。

また、局所的に大きすぎないかどうかをチェックするために、この付加質量をポスト処理することが重要です。大きすぎる場合、結果が正しくないことを示している可能性があります（これをチェックするには、/ANIM/NODA/DMASで書き出されたアニメーションを可視化する必要があります）。

また、/DT/NODA/CSTにより付加慣性が発生することもあります。エンジニアリングの観点から言うと、これは評価がさらに難しく、Radioss V9.0より前のバージョンではポスト処理できません。

多くの場合、/DT/NODA/CSTによる付加質量が計算の発散を引き起こします。

インターフェースの定義が不適切な場合（小さすぎるギャップ、初期貫通など）、/DT/NODA/CSTはタイムステップを増やすためのソリューションとはなりません。モデルを改善する必要があります。/DT/NODA/CSTと/DT/INTER/CSTのどちらも使用しないことが推奨されます。

最初に/DT/NODA/CSTを使用して、付加質量と共に短く実行し（１サイクルで十分）、その後、時間t=0に書き出されたアニメーションをポスト処理することにより付加質量をチェックすることが可能です。その後、より適切なタイムステップを得ながら、問題を減らすために、/DT/NODA/CSTを使用せずに計算を続行できます。