ACU-T：5200 チェックバルブの剛体力学

このチュートリアルでは、圧力チェックバルブの開口部のシミュレーションの設定、解析、および結果表示のための手順を説明します。このシミュレーションでは、AcuSolveを使用して、経時変化する流入の流れ場によりバルブにかかる力を計算し、これらの流体力によって生じるバルブの動きを計算します。このチュートリアルの目的は、剛体力学のシミュレーションを実行するために必要ないくつかのモデリングの概念を紹介することです。

前提条件

このチュートリアルを実行する前に、HyperWorksの入門チュートリアルであるACU-T：1000 HyperWorksユーザーインターフェースをすでに完了し、HyperWorks CFD、AcuSolve、およびHyperViewの基本を理解しているものとします。この解析を実行するには、ライセンス供与済みバージョンのHyperWorks CFDとAcuSolveにアクセスできる必要があります。

このチュートリアルを実行する前に、HyperWorksCFD_tutorial_inputs.zipを<Altair_installation_directory>\hwcfdsolvers\acusolve\win64\model_files\tutorials\AcuSolveから作業ディレクトリにコピーします。 ACU-T5200_pressureCheckValve.x_t をHyperWorksCFD_tutorial_inputs.zipから取り出します。

問題の説明

このチュートリアルで扱う問題は、図 1に図示しています。この問題は、仮想スプリング（形状には含まれていない）にシャッターが取り付けられたチェックバルブを通過して流れる水が含まれた円筒管で構成されています。入口圧力は時間とともに変化し、シャッターの動きは、スプリングの反力に対する流体力のバランスに応じて決まります。この問題は、長手方向の軸を中心とした30度ごとの回転周期性を持ち、生じる流れも回転周期性を持つことが想定されるため、くさび形の部分を使用したモデリングが可能になります。このチュートリアルでは、図に示すように30度分の形状をモデル化します。回転周期性を持つ形状の一部をモデル化することで、正確な解を維持しつつ、計算時間を短縮できます。



図 1. スプリング式シャッターがあるチェックバルブの図

管の入口直径は0.08mで、長さは0.3 mです。チェックバルブのアセンブリは、入口の下流0.085mにあります。これは、中心に配置された直径0.044mの開口部と開口部から0.005m離れた初期位置にあるシャッターを備えた0.005m厚の薄板で構成されており、ほぼ閉じられた状態をシミュレートします。シャッター薄板の直径は0.05mで、厚みは0.005mです。シャッター薄板は、長さが0.03mで直径が0.01mのステムに接合されています。シャッターとステムの質量は0.2kgであり、その動きは剛性が2162N/mである仮想スプリングの影響を受けます。バルブシャッターの動きは、シャッターの下流にある穴あき薄板に取り付けられた止め具によって制限されます。

AcuSolveの内部剛体力学ソルバーは、接触をシミュレートできないことに留意してください。したがって、この問題はバルブと止め具との間の接触を避けるために定式化されます。

30度の部分として形状をモデル化するには、流体モデルが剛体モデルと一致するように設定される必要があります。剛体の1/12のみがモデル化されるため、AcuSolveによって計算されるバルブシャッターに作用する力は、装置にかかる実際の力の1/12に相当します。したがって、シミュレーションでもこの点を考慮に入れる必要があります。そのための方法は2つあります。
  1. 1つ目は、AcuSolveによって計算される流体力を12倍にすることです。これにより、ボディの変位が計算される際にデバイスにかかる全荷重が表されます。

    この方法を使用すると、バルブスプリングの全剛性が剛体解で使用され、バルブの全質量が使用されます。

  2. 2つ目は、バルブの質量とスプリングの剛性を12分の1に縮小して、モデル化されるバルブ形状の比率と一致させることです。

    この方法を使用した場合は、剛体ソルバーに渡される荷重はスケーリングされません。

このチュートリアルではこの２つ目の方法が使用されるため、縮小された質量0.0167kgと剛性180.1667N/mが使用されます。

この問題で扱う流体は水であり、密度 (ρ)は1000 kg/m3、分子粘性 (μ)は1.781X 10-3kg/m-secです。

シミュレーションの開始時には、流れ場は静止しています。流れは、図 2に示す区分線形関数として経時変化する入口圧力によって引き起こされます。入口圧力が上昇するにつれて、バルブが開くため流れは加速されます。乱流粘性率は10と想定されています。

初期の入口圧力は0Paです。0.002秒で圧力が増加し始め、0.05秒で26,500Paに達します。圧力は26,500Paで維持され、0.2秒から減少し始めて、0.25秒で初期圧力に達し、シミュレーションの残りについてはそのままの圧力が維持されます。


図 2. 入口の過渡圧力

前に行ったこの形状のシミュレーションでは、入口での平均速度が最大の0.9 m/sに達することが示されています。この速度では、流れのレイノルズ数は72,000となります。レイノルズ数が4,000より大きい場合は、一般に流れを乱流としてモデル化することができます。

流れの初期状態は実際には層流ですが、バルブシャッター周囲の流量と流速の増大によって、乱流状態への急速な変化が生じると予想されます。したがって、このシミュレーションは、過渡的な乱流をモデル化するように設定されます。非定常解析を実行する際は、定義済みのスタッガー基準に基づいて時間ステップごとに収束に達します。メッシュモーションは、任意のメッシュ移動（ALEメッシュモーション）を使用してモデル化されます。

このケースでは、対象となるこの過渡挙動は、過渡圧力プロファイルによって与えられている、圧力が上昇してから再び低下するのにかかる時間内に生じます。スプリングが元の状態に戻るための時間を見越して、追加の時間がシミュレートされます。このチュートリアルでは、圧力が降下して初期状態に戻った後に0.1秒が追加されるため、合計時間は0.35秒となります。

非定常解析において極めて重要なもう１つの決定事項は、時間増分の選択です。時間増分は、シミュレーションの指定された時間ステップにおける時間変化です。対象となる流動特性の変化を捕捉するのに十分な短さで、不要な計算時間がかからない時間増分を選択することが重要です。

初期状態から最大値への入口圧力の変化は、0.048秒を超える時間にわたって発生します。時間増分が0.002秒の場合は、過大な計算時間を必要とすることなく、この過渡変化の優れた解析が可能になります。

HyperWorks CFDの起動とHyperMeshモデルデータベースの作成

  1. Windows のスタートメニューからスタート > Altair <バージョン> > HyperWorks CFDをクリックして HyperWorks CFDを起動します。
    HyperWorks CFDが読み込まれると、Geometryリボンが表示されます（デフォルト）。
  2. 以下の方法のいずれかで新規.hmデータベースを作成します。
    • メニューバーで、File > Saveをクリックします。
    • HomeツールのFilesツールグループからSave Asツールをクリックします。


      図 3.
  3. Save File Asダイアログで、データベースを保存したいディレクトリを指定します。
  4. データベース名として CheckValve_RBDを入力して、Saveをクリックします。
    このディレクトリが解析用ディレクトリになり、シミュレーションに関連するすべてのファイルがこの場所に保存されます。

形状のインポート検証

形状のインポート

  1. メニューバーからFile > Import > Geometry Modelをクリックします。
  2. Import Fileダイアログで、作業ディレクトリに移動し、 ACU-T5200_pressureCheckValve.x_tを選択してOpenをクリックします。
  3. Geometry Import Optionsダイアログで、すべてのオプションをデフォルト設定にしたままImportをクリックします。


    図 4.


    図 5.

形状の検証

  1. GeometryリボンからValidateツールをクリックします。


    図 6.
    Validateツールでは、モデル全体をスキャンしてサーフェスおよびソリッドのチェックを実行し、フリーエッジ、閉じたシェル、交差、重複、細長いサーフェスなど、形状内の欠陥にフラグを付けます。

    現在のモデルには、上記のような問題は存在しません。問題が見つかった場合は、ツール名の横の括弧内の数で示されます。

    Validateアイコンの左上に青色のチェックマークが表示されているのがわかります。これは、このツールでは形状モデルの問題は検出されなかったことを示しています。


    図 7.
  2. Escまたは、モデリングウィンドウで右クリックして、緑色のチェックマークの上でカーソルを右から左にスワイプしてください。
  3. データベースを保存します。

流れのセットアップ

シミュレーションパラメータとソルバー設定のセットアップ

  1. FlowリボンからPhysicsツールをクリックします。


    図 8.
    Setupダイアログが開きます。
  2. Physics modelsの設定で
    1. Time marchingを Transientに設定します。
    2. Time step sizeを0.002に設定し、Final timeを0.35に設定します。
    3. Turbulence モデルとしてSpalart-Allmarasを選択します。


    図 9.
  3. Solver controls設定をクリックし、Maximum stagger iterationsを3に設定します。
    図 10.
  4. ダイアログを閉じてモデルを保存します。

入口圧力のMultiplier Functionの作成

  1. Flowリボンから、Setupツールセットの横にある矢印をクリックし、Multipliersを選.択します。


    図 11.
  2. Multiplier Libraryダイアログで をクリックします。
  3. multiplier作成のダイアログで、左上隅をクリックしてMultiplier Functionの名前をInlet Pressureに変更します。
  4. TypeをPiecewise Linearに設定します。
  5. VariableがTimeに、EvaluationがTime Stepに設定されていることを確認します。
  6. を4回クリックし、テーブルの最後に4つの行を追加します。
  7. 下の画像に従ってテーブルの値を入力します。


    図 12.

材料プロパティの割り当て

  1. FlowリボンからMaterialsツールをクリックします。


    図 13.
  2. Pipeジオメトリの任意の場所をクリックします。
    ジオメトリ全体がハイライト表示されます。


    図 14.
  3. マイクロダイアログで、MaterialドロップダウンメニューからWaterを選択します。
  4. ガイドバーでをクリックすると、コマンドが実行されてツールが終了します。
  5. モデルを保存します。

流れ境界条件の定義

この手順では、入口、出口、および対称フェイスの境界条件を割り当てます。
  1. FlowリボンのPressureツールグループからStagnation Pressureツールをクリックします。


    図 15.
  2. モデリングウィンドウで、図 16 に示すように入口サーフェスをクリックします。
  3. マイクロダイアログで、Stagnation pressureに1N/m2という値を入力します。
  4. Multiplier functionドロップダウンをクリックし、リストから関数Inlet Pressureを選択します。


    図 16.
  5. マイクロダイアログ内のTurbulenceアイコンをクリックし、Turbulence input typeをViscosity Ratioに設定し、Turbulence viscosity ratioに10という値を入力します。


    図 17.
  6. ガイドバーでをクリックすると、コマンドが実行されますが、ツールは終了しません。
  7. Boundaries凡例でStagnation pressureを右クリックし、これをInletという名前に変更して、Enterキーを押します。
  8. Outletツールをクリックします。


    図 18.
  9. 下の図に示されている出口サーフェスをクリックします。


    図 19.
  10. ガイドバーでをクリックすると、コマンドが実行されてツールが終了します。
  11. Symmetryツールをクリックします。


    図 20.
  12. 下の図に示されている2つのサーフェスを選択します。


    図 21.
  13. ガイドバーでをクリックすると、コマンドが実行されますが、ツールは終了しません。
  14. Boundaries凡例で、SymmetryをFront_symmetryという名前に変更します。
  15. モデルを回転して、あと2つの対称フェイスを選択します。


    図 22.
  16. ガイドバーでをクリックすると、コマンドが実行されますが、ツールは終了しません。
  17. Boundaries凡例でSymmetryをBack_symmetryという名前に変更します。
  18. モデルを保存します。

モーションのセットアップ

この手順では、メッシュモーションをアクティブにし、バルブ壁の剛体モーションを定義します。その後、対称サーフェスのメッシュ変位境界条件を定義します。

メッシュモーションタイプの定義

  1. MotionリボンからSettingsツールをクリックします。


    図 23.
  2. Setupダイアログで、Mesh motionをArbitraryに変更します。


    図 24.

剛体モーションの定義

この手順では、バルブ壁の剛体メッシュモーションを定義します。
  1. MotionリボンからRigid Bodyツールをクリックします。


    図 25.
  2. ガイドバーで、Surfacesセレクターをクリックしてアクティブにします。
  3. 作業しやすいように、モデリングウィンドウでFront_symmetryとBack_symmetryのフェイスを非表示にします。このためには、これら4つのサーフェスを選択してキーボードのHキーを押すか、右クリックしてコンテキストメニューでHideを選択します。


    図 26.
  4. ウィンドウ選択（左マウスボタンを押したままドラッグします）により、下の図のようにバルブ壁サーフェスを選択します。
    選択されたサーフェスの数は8つのはずです。これはガイドバー上に表示されます。


    図 27.
  5. ガイドバーのDOF上でクリックし、モデリングウィンドウ内のバルブ壁サーフェス上でクリックします。
  6. 表示されるマイクロダイアログで、次の値を入力します：
    • Mass：0.0167kg
    • Center：(0, 0, 0)


    図 28.
  7. モデリングウィンドウで、X-displacement矢印をクリックします。
  8. マイクロダイアログで、ばねの剛性Kxxに180.1667kg/sec2の値を入力します。


    図 29.
  9. X-rotation矢印をクリックし、これを無効にします。


    図 30.
  10. 同様に、残りの4つの矢印（Y-displacement、Y-rotation、Z-displacement、およびZ-rotation）をクリックし、これらを無効にします。
    これにより、バルブ本体がx方向にのみ移動するようになります。
  11. ガイドバーでをクリックすると、コマンドが実行されてツールが終了します。
  12. すべてのサーフェスの表示をオンにして、モデルを保存します。

メッシュ変位境界条件の定義

  1. MotionリボンからPlanar Slipツールをクリックします。


    図 31.
  2. 図 32 に示されている2つの前方対称サーフェスを選択します。
  3. Mesh Motion凡例で、Planar Slipという名前をダブルクリックしてFront_symmetryに変更します。


    図 32.
  4. ガイドバーでをクリックすると、コマンドが実行されますが、ツールは終了しません。
  5. モデルを回転して、2つの後方対称フェイスを選択します。
  6. Mesh Motion凡例でPlanar SlipをBack_symmetryという名前に変更します。
  7. ガイドバーでをクリックすると、コマンドが実行されてツールが終了します。
  8. モデルを保存します。

メッシュの生成

このステップでは、メッシュコントロールを指定し、メッシュを作成します。

ゾーンメッシュコントロールの定義

ゾーンメッシュコントロールにより、円筒、球、立方体のゾーンなど、標準形状のゾーンを使用して、特定の領域のローカルボリュームメッシュサイズを定義できます。
  1. MeshリボンのZonesツールグループからCylinderツールをクリックします。


    図 33.
  2. モデリングウィンドウで、下の図に示されているポイントの周辺にマウスを置きます。円筒プレビューゾーンが管の軸と平行な場合、下に示されているポイントの近くのモデルをクリックします。
    このポイントは、円筒の前方フェイスの中心となります。次のいくつかの手順では、このポイントの座標を手動で編集します。


    図 34.
  3. マウスカーソルを中心から離し、再度クリックします。
    ゾーンのプレビューがマニピュレーターとともに画面に表示されます。マニピュレーターを使用して、ゾーンの位置および向きを変更することができます。


    図 35.
  4. マニピュレーターの中心をクリックし、その中心として（-0.06, 0.02, 0）の座標を入力します。


    図 36.
  5. ゾーンの円筒サーフェスをダブルクリックし、円筒の半径として0.022mと入力します。


    図 37.
  6. 円筒の底面をダブルクリックし、円筒の高さとして0.1mと入力します。


    図 38.
  7. モデリングウィンドウ内の空白部分をクリックします。
  8. マイクロダイアログのmesh sizeに0.0015mの値を入力します。
  9. Refinement zones凡例で、ゾーンの名前をValve region mesh refinementに変更します。


    図 39.
  10. モデルを保存します。

境界層コントロールの定義

  1. MeshリボンからBoundary Layerツールをクリックします。


    図 40.
  2. モデリングウィンドウで右クリックしてSelect > Advanced Select > By Boundaries > Default Wallの順に移動します。
    すべての壁サーフェスをハイライト表示する必要があり、BL指定用のダイアログが表示されます。
  3. ダイアログに次の値を入力します：
    1. First layer thickness：0.00035
    2. Total number of layers：3
    3. Growth method：Constant
    4. Growth rate：1.2
    5. Termination policy：Truncate
    6. Enable surface mesh modificationをアクティブにします。


    図 41.
  4. BL Controls凡例でBL controlをPipe walls BLという名前に変更します。
  5. ガイドバーでをクリックすると、コマンドが実行されますが、ツールは終了しません。
  6. 次のいくつかの手順では、バルブ壁の境界層メッシュコントロールを定義します。作業しやすいように、前方および後方の対称フェイスは非表示にします。
  7. ウィンドウ選択の手法を使用して、バルブ壁サーフェスを選択します。
  8. 表示されたダイアログで、BL仕様として次の値を入力します：
    1. First layer thickness：0.00015
    2. Total number of layers：3
    3. Growth method：Constant
    4. Growth rate：1.2
    5. Termination policy：Truncate
    6. Enable surface mesh modificationをアクティブにします。


    図 42.
  9. BL Controls凡例でBL controlをValve walls BLという名前に変更します。
  10. ガイドバーでをクリックすると、コマンドが実行されてツールが終了します。
  11. Mesh Controlsツールセットの横のドロップダウンをクリックし、ツールセットを選択し、Advancedを選択します。


    図 43.
  12. Advanced mesh settingsダイアログで、Boundary layer element typeをTetraに変更し、ダイアログを閉じます。


    図 44.
  13. すべてのサーフェスの表示をオンにして、モデルを保存します。

メッシュの生成

  1. MeshリボンからBatchツールをクリックします。


    図 45.
  2. Meshing Operationsダイアログで、Mesh sizeオプションをMaximum sizeに、Maximum element sizeを0.005に設定します。
  3. Curvature based surface refinementを無効にします。
  4. Meshをクリックし、メッシュを生成します。


    図 46.

    メッシュプロセスが開始されるとRun Statusダイアログが表示され、アプリケーションがSolutionリボンに移動します。

  5. メッシュが完了したら、Run Statusダイアログを閉じます。

解の計算

節点出力の定義

  1. SolutionリボンからFieldツールをクリックします。


    図 47.
  2. Field Outputsダイアログで、Solution variablesのTime step intervalを3に設定します。


    図 48.

AcuSolveの起動

  1. SolutionリボンからRunツールをクリックします。


    図 49.
  2. Launch AcuSolveダイアログで、Additional argumentsボックスに次のテキストを入力します：-tlog -lprobe
    これにより、実行が開始されたときに自動的にAcuProbeおよびAcuTailユーティリティが開始されます。
  3. Parallel processingオプションをIntel MPIに設定します。
  4. オプション: プロセッサーの数を、利用環境に合わせ、4または8に設定します。
  5. Automatically define pressure referenceオプションを無効にします。
  6. 他のオプションはデフォルト設定のままとして、Runをクリックして解析プロセスを開始します。


    図 50.
    Run Statusダイアログが開きます。解析が終了すると、ステータスが更新され、ダイアログを閉じることができます。

AcuProbeによる解析のモニター

AcuSolveの実行中は、AcuProbeを使用して入口圧力、バルブの変位、およびバルブの速度などの流動特性を確認できます。

  1. AcuProbeウィンドウに移動します。Data Treeで、Surface Output > Rigid Body – surface output > Geometryの順に拡張表示します。
  2. mesh_x_displacementを右クリックし、Plotを選択します。
  3. 同様にしてSurface OutputでInlet - Output > Pressureの順に拡張表示します。
  4. total_pressureを右クリックし、Plotを選択します。
    注: プロットを正しく表示するために、ツールバーで をクリックする必要がある場合があります。


    図 51.

HyperViewによる結果のポスト処理

この手順では、水がバルブを流れる際のバルブの動きのアニメーションを作成します。

HyperViewのオープンとモデルおよび結果の読み込み

  1. WindowsのスタートメニューからStart > All Programs > Altair <version> > HyperViewをクリックして、HyperViewを起動します。
    HyperViewウィンドウを読み込むと、デフォルトでLoad model and resultsパネルが開きます。このパネルが表示されない場合は、File > Open > Modelの順にクリックします。
  2. Load model and resultsパネルで、Load modelの隣にある をクリックします。
  3. Load Model Fileダイアログで、作業ディレクトリに移動して、ポスト処理する解析実行のAcuSolve .logファイルを選択します。この例で選択するファイルは、CheckValve_RBD.1.Logです。
  4. Openをクリックします。
  5. パネル領域でApplyをクリックしてモデルと結果を読み込みます。
    読み込むと、モデルが形状で色分けされます。

速度のアニメーションの作成

  1. ResultsブラウザでComponentsのリストを拡張表示し、Isolate Shownアイコン をクリックします。
  2. Front_symmetry - Outputコンポーネントをクリックして、前方対称サーフェスを除くすべてのコンポーネントの表示をオフにします。


    図 52.
  3. Standard Viewsツールバーの をクリックすることで、xy平面を正面から見た表示にします。
  4. Resultsツールバーで をクリックしてContourパネルを開きます。
  5. パネル領域でResult typeをVelocity (v)に設定し、下のドロップダウンがMag（大きさ）に設定されていることを確認します。
  6. Componentsエンティティセレクターをクリックします。


    図 53.
  7. Extended Entity SelectionダイアログでDisplayedを選択します。
  8. パネル領域でApplyをクリックします。
  9. Legendタブに移動し、つづいてEdit Legendをクリックします。
  10. Edit Legendダイアログで、TypeをDynamic scale、Numeric formatをFixedにそれぞれ変更し、OKをクリックします。
  11. Animationツールバーで をクリックし、前方対称面上の速度のアニメーションを再生します。


    図 54.
  12. ImageCaptureツールバーで、Capture Graphics Area Videoアイコン をクリックします。
  13. Save Graphics Area Video Asダイアログでアニメーションを保存するディレクトリを参照し、このビデオにvelocity magなどの名前を付けてSaveをクリックします。

断面上の圧力コンターおよび速度ベクトルの表示

次の手順では、中央Z平面に断面を作成し、その断面上の圧力コンターと速度ベクトルを表示します。

  1. Resultsブラウザで、Fluidコンポーネント（Fluid_Main_SolidBody_2_1）以外のすべてのコンポーネントの表示をオフにします。


    図 55.
  2. Resultsツールバーで をクリックしてContourパネルを開きます。
  3. パネル領域でResult typeをPressure (s)に設定します。
  4. Components エンティティセレクターをクリックします。Extended Entity Selection ダイアログでDisplayedを選択します。
  5. Applyをクリックして圧力のコンタープロットを作成します。
  6. パネル領域のResultタブの下で、Overlay result displayチェックボックスをアクティブにします（設定されていない場合）。


    図 56.
  7. HV-DisplayツールバーのSection cutアイコン をクリックします。
  8. パネル領域でAddをクリックし、Section 1という名前の新しい断面を作成します。
  9. Define planeセクションで、軸をZ Axisに設定し、Applyをクリックします。
  10. Z Base coordinateを0に設定し、Enterキーを押します。
  11. DisplayオプションをClipping planeからCross sectionに変更します。
  12. Gridlineをクリックします。Gridline OptionsダイアログのGrid lineの下のShowチェックボックスを無効にし、OKをクリックします。
  13. ResultsツールバーのVectorアイコン をクリックして、Vectorパネルを開きます。
  14. パネル領域でResult typeをVelocity (v)に設定します。
  15. Selectionドロップダウンをクリックし、オプションリストからSectionsを選択します。


    図 57.
  16. Sectionsエンティティセレクターをクリックし、Allを選択します。
  17. パネル領域で、Overlay result displayチェックボックスをアクティブにします（設定されていない場合）。
  18. Applyをクリックします。
  19. Plotタブの下で、X+Y+Z Resultantオプションのみが選択されていることを確認します。
  20. Displayタブに移動し、Size scalingオプションをUniformに設定して、サイズ欄に0.002の値を入力します。
  21. Color byオプションをDirectionに設定し、X+Y+Zの色を白に設定します。


    図 58.
  22. Sectionタブに移動し、Projectedチェックボックスをアクティブにして、Applyをクリックします。
    ベクトルプロットは下の図のようになるはずです。0.156秒での結果が示されています。


    図 59.

要約

このチュートリアルでは、剛体メッシュモーションを伴うCFDシミュレーションを設定して解析する方法を知ることができました。まず、形状をインポートし、流れの設定を定義しました。次に、メッシュモーション設定を定義し、メッシュを生成しました。解が計算されたら、AcuProbeとHyperViewを使用して結果をポスト処理しました。