OS-T：1330 1/2車室モデルの音場解析

本チュートリアルの目的は、流体-構造連成を受ける1/2車室モデルの振動特性を評価することにあります。参照される流体は、airです。特に、車室内での流体内の主な応答の場所である、ドライバーの耳位置近くでの騒音または音響レベルが評価されます。

開始する前に、このチュートリアルで使用するファイルを作業ディレクトリにコピーします。
図 1. １/２車室モデル

rd2060_half_car_model

1/2車室モデルが、Figure 1の赤い拘束記号（三角形）で示したように底部で加振されます。加振は、車室の高さ方向（Z-軸）に沿った単位荷重で与えられます。

HyperMeshの起動とOptiStructユーザープロファイルの設定

  1. HyperMeshを起動します。
    User Profilesダイアログが現れます。
  2. OptiStructを選択し、OKをクリックします。
    これで、ユーザープロファイルが読み込まれます。ユーザープロファイルには、適切なテンプレート、マクロメニュー、インポートリーダーが含まれており、OptiStructモデルの生成に関連したもののみにHyperMeshの機能を絞っています。

モデルのオープン

  1. File > Open > Modelをクリックします。
  2. 自身の作業ディレクトリに保存したHalf_Car.hmファイルを選択します。
  3. Openをクリックします。
    Half_Car.hmデータベースが現在のHyperMeshセッションに読み込まれます。

モデルのセットアップ

材料とプロパティの生成およびそれらの構造と流体要素への割り当て

  1. Model Browser内で右クリックし、Create > Materialを選択します。
  2. Curve editorウィンドウが表示されます。Newをクリックし、Name=にMAT1_shellsと入力してproceedをクリックします。
  3. Card Imageに、ドロップダウンメニューからMAT1を選択します。
  4. E、Nu、Rhoの欄にそれぞれ 2.1e04、0.33、 8.0e-10と入力します。
    図 2.

    os_1330_05
  5. Model Browser内で右クリックし、Create > Materialを選択します。
  6. Curve editorウィンドウが表示されます。Newをクリックし、Name=にMAT10_Solidsと入力してproceedをクリックします。
  7. Card Imageに、ドロップダウンリストからMAT10を選択します。
  8. RhoとCの欄にそれぞれ1.2e-13と3.4e5を入力します。
  9. Model Browser内で右クリックし、Create > Propertyを選択します。
  10. NameにShellsと入力します。
  11. Card Imageに、ドロップダウンメニューからPSHELLを選択します。
  12. Materialに、Unspecified > Materialをクリックします。
  13. Select Materialダイアログで、材料のリストからMAT1_shellsを選択し、OKをクリックして選択を完了します。
  14. シェルコンポーネントの板厚として、Tをクリックし、2.0と入力します。
  15. Model Browser内で右クリックし、Create > Propertyを選択します。
  16. NameにSolidsと入力します。
  17. Card Imageに、ドロップダウンメニューからPSOLIDを選択します。
  18. Materialに、Unspecified > Materialをクリックします。
  19. Select Materialダイアログで、MAT10_Solidsを選択します。
  20. FCTNにPFLUIDを選択します。
  21. fluidコンポーネントをクリックします。
    コンポーネントのエントリがEntity Editorに表示されます。
  22. Propertyに、Unspecified > Propertyをクリックします。
  23. Select Propertyダイアログで、プロパティsolidsを選択します。
  24. structureコンポーネントをクリックします。
    コンポーネントのエントリがEntity Editorに表示されます。
  25. Propertyに、Unspecified > Propertyをクリックします。
  26. Select Propertyダイアログで、プロパティshellsを選択します。

Create Load Collectors

In this step the model is unconstrained and a unit vertical load is applied acting upwards in the positive z-direction at a point on the base of the car (shown in page 1). The model can be unconstrained as the solver applies PARAM, INREL -2 by default to avoid the model from experiencing a rigid body motion.

  1. Model Browserで右クリックしcontext menuからCreate > Load Collectorを選択します。
    デフォルトの荷重コレクターがEntity Editorに表示されます。
  2. For Name, enter unit-load.
  3. Colorをクリックし、カラーパレットから色を選択します。
  4. Set Card Image to None and click Close.
  5. Verify unit_load is the current load collector. If it is not the current load collector, right-click on unit_load in the Model Browser and select Make Current from the context menu.
    図 3.


    ヒント: In the Model Browser, Load Collectors folder, the current load collector is bold.
  6. Click unit-load > return.

1つの点における単位荷重の生成

  1. Analysisページからconstraintsをクリックします。
  2. パネル左側のラジオボタンを用いてcreateサブパネルを選択します。
  3. nodes >> by idをクリックし、車両の節点19072を選択します。
  4. dof3を除いた全ての自由度からチェックマークを外します。
  5. dof3の右にある=をクリックし、値1を入力します。
  6. Load Types =に、拡張エンティティ選択メニューからDAREAを選択します。
  7. createをクリックします。
    これで選択された節点の単位荷重が付与されます。
    図 4.

    rd2060_unit_load_point
  8. returnをクリックします。

Create a Frequency Range Table

  1. In the Model Browser, right-click and select Create > Curve.
    A new window opens.
  2. For Name, enter tabled1.
  3. In the table, enter x(1) = 0.0, y(1) = 1.0, x(2) = 200, y(2) = 1.0.
  4. Close the Curve Editor window.
  5. From Curves, select tabled1.
  6. Click Color and select from the palette.
  7. For Card Image, select TABLED1 from the drop-down menu.
    This provides a frequency range of 0.0 to 200 with a constant 1.0 over this range.
    図 5.


周波数依存動的荷重の作成

  1. Model Browser内で右クリックし、Create > Load Step Inputsを選択します。
  2. Curve editorウィンドウが表示されます。Newをクリックし、Name=にrload1と入力してproceedをクリックします。
  3. Config typeに、ドロップダウンメニューからDynamic Load – Frequency Dependentを選択します。
  4. TypeにドロップダウンメニューからRLOAD1を選択します。
  5. EXCITEIDに、Loadcolをクリックします。
  6. Select Loadcolダイアログで、荷重コレクターのリストから unit-loadを選択し、OKをクリックして選択を完了します。
  7. TCフ欄で、カーブtabled1を選択します。

    加振のタイプには、荷重（力またはモーメント）、強制変位、速度、または加速度を与えることができます。RLOAD1荷重ステップ入力のTYPE欄で荷重のタイプを定義します。タイプはデフォルトでapplied loadにセットされます。

    典型的なRLOAD1カードが下のように現われます。
    図 6.

    OS_1330_02

  1. Model Browserで右クリックしてCreate > Load Collectorを選択します。
  2. Curve editorウィンドウが表示されます。Newをクリックし、Name=にfreq1と入力してproceedをクリックします。
  3. Colorをクリックし、カラーパレットから色を選択します。
  4. Card Imageに、ドロップダウンメニューからFREQiを選択します。
  5. FREQ1の隣のボックスにチェックを入れます。
  6. NUMBER_OF_FREQ1に値1を入力し、キーを押します。
  7. Data欄の横のtable_pencilをクリックし、F1= 0.0、DF= 1.0、NDF= 200と入力します。
    これで周波数0.0から1.0刻みで200の周波数増分を与えたことになり、GUIに下のようにカードが現われます。
    図 7.

    OS_1330_03

  1. Model Browserで右クリックしてCreate > Load Collectorを選択します。
  2. Curve editorウィンドウが表示されます。Newをクリックし、Name=にeigrl1と入力してproceedをクリックします。
  3. Colorをクリックし、カラーパレットから色を選択します。
  4. Card Imageに、ドロップダウンメニューからEIGRLを選択します。
  5. V2に値600.0を入力します。
  6. NDに値50を入力します。

    これで Lanczos 法を用いた固有値の取り出しに最初の周波数から 600 Hz の間の周波数範囲を指定したことになります。

  7. 同様に、手順1.から6.までに従って、 eigrl2という名称の荷重コレクターを生成します。
    図 8.

    OS_1330_04

  1. Model Browserで右クリックしてCreate > Load Stepを選択します。
  2. Curve editorウィンドウが表示されます。Newをクリックし、Name=にsubcase1と入力してproceedをクリックします。
  3. Colorをクリックし、カラーパレットから色を選択します。
  4. Analysis typeに、ドロップダウンメニューからFreq.resp (modal)を選択します。
  5. METHOD(STRUCT)に、荷重コレクターのリストからeigrl1を選択します。
  6. METHOD(FLUID)に、荷重コレクターのリストからeigrl2を選択します。
  7. DLOADに、荷重コレクターのリストからrload1を選択します。
  8. FREQに、荷重コレクターのリストからfreq1を選択します。

    拘束条件と単位荷重として荷重コレクターrload1を周波数定義の荷重コレクターfreq1とモーダル法が定義された荷重コレクターeigrlと共に参照したOptiStructサブケースが生成されました。

  1. Model Browserで右クリックしてCreate > Setを選択します。
  2. NemeにSETAと入力します。
  3. Card Imageに、ドロップダウンメニューからNoneを選択します。
  4. Set Typeスイッチはnon-orderedにセットされたままにしておきます。
  5. Entity IDsに、黄色いNodesパネルをクリックし、ID 18881の節点を選択します。
  6. proceedをクリックします。

Create a Set of Frequencies

  1. In the Model Browser, right-click and select Create > Load Collector.
  2. For Name, enter freq1.
  3. Colorをクリックし、カラーパレットから色を選択します。
  4. For Card Image, select FREQi from the drop-down menu.
  5. Check the FREQ1 option and enter 1 in the NUMBER_OF_FREQ1 field.
  6. Update the following fields in the pop-out window.
    1. For F1, enter 0.0.
    2. For DF, enter 1.0.
    3. For NDF, enter 200.
  7. Click Close.
    This provides a set of frequencies beginning with 0.0, incremented by 1.0 and 200 frequencies increments and the card appears as shown below on the GUI.

固有値解析のためのモーダル法の生成

  1. Model Browser内で右クリックし、Create > Load Step Inputsを選択します。
  2. Curve editorウィンドウが表示されます。Newをクリックし、Name=にeigrl1と入力してproceedをクリックします。
  3. Config typeにReal Eigen Value Extractionを選択します。
  4. TypeにドロップダウンメニューからEIGRLを選択します。
  5. V2に値600.0を入力します。
  6. NDに値50を入力します。

    これで Lanczos 法を用いた固有値の取り出しに最初の周波数から 600 Hz の間の周波数範囲を指定したことになります。

  7. 同様に、eigrl2という名称の荷重ステップ入力を作成します。
    図 9.

    OS_1330_04

Create a Load Step

  1. In the Model Browser, right-click and select Create > Load Step.
    A default load step template is now displayed in the Entity Editor below the Model Browser.
  2. For Name, enter subcase1.
  3. For Analysis type, select Freq.resp (modal) from the drop-down menu.
  4. For METHOD(STRUCT), select Unspecified > Load step inputseigrl1.
  5. For METHOD(FLUID), select eigrl2 from the list of load step inputs.
  6. For DLOAD, select rload1 from the Select Load Step Inputs pop-out window.
  7. For FREQ, click Unspecified > Loadcol
  8. From the Select Loadcol dialog, select freq1.
    An OptiStruct subcase has been created which references the constraints, the unit load in the load step input rload1 with a set of frequencies defined in load collector freq1 and modal method defined in the load step input eigrl.

Create a Set of Nodes

  1. In the Model Browser, right-click and select Create > Set.
  2. For Name, enter SETA.
  3. For Card Image, select SET_GRID.
  4. Leave the Set Type switch set to non-ordered type.
  5. For Entity IDs, select Nodes from the selection switch.
  6. Click Nodes and select nodes 18881.
  7. Click proceed.

出力セットの生成

  1. Analysisページからcontrol cardsをクリックします。
  2. ACMODLをクリックします。
    これでモデルの流体-構造相互作用のパラメータを定義します。
  3. [INTER] をクリックしDIFFを選択します。
  4. [INFOR]をクリックしALLを選択します。
  5. returnをクリックし、このメニューを終了します。
  6. GLOBAL_OUTPUT_REQUESTを選択します。続いて、DISPLACEMENTの左のボックスにチェックマークを入れます。
    作業領域のスクリーンに新しいウィンドウが現われます。
  7. 入力ボックスFORMをクリックし、ポップアップメニューからPHASEを選択します。
  8. 入力ボックスOPTIONをクリックし、ポップアップメニューからSIDを選択します。
    黄色い新しい欄が現われます。
  9. 黄色いSIDボックスをダブルクリックし、左下隅のポップアップ選択からSETAを選択します。
    値1がSID入力ボックスの下に現れます。これは、set 1内の節点に限った出力を設定しています。
  10. returnをクリックし、このメニューを終了します。
  11. GLOBAL_CASE_CONTROLを選択します。
  12. FREQの横のボックスにチェックを入れます。
  13. FREQをクリックし、荷重コレクターfreq1を選択します。
  14. returnをクリックしてこのメニューを終了し、nextをクリックします。
  15. OUTPUTサブパネルを選択します。
    作業領域に新しいウィンドウが現われます。
  16. number of outputs = 4と指定します。
  17. KEYWORDがHGFREQにセットされていることを確認します。
    HGFREQの使用により、周波数がHyperGraph形式で出力されます。
  18. FREQの下のボックスをクリックし、ポップアップ選択からALLを選択します。
    ALLを選択することにより、全ての周波数の結果が出力されます。
  19. KEYWORDがOPTIにセットされていることを確認します。
  20. FREQの真下のボックスをクリックし、ポップアップ選択からALLを選択します。
  21. 同様にKEYWORDの下でPUNCHとH3Dを選択します。
  22. returnをクリックし、このメニューを終了します。
  23. PARAMを選択します。
  24. AUTOSPCをクリックします。
  25. スクロールダウンしてGの横のボックスにチェックを入れます。
    作業領域のスクリーンに新しいウィンドウが現われます。
  26. G_V1入力ボックスに値0.06を入力します。
    この値は一様な構造減衰の係数を指定し、臨界減衰[ C / C 0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qaiaac+ cacaWGdbWaaSbaaSqaaiaaicdaaeqaaaaa@391F@ ]比に2.0を掛けることにより得られます。
  27. GFLの横のボックスにチェックを入れます。
  28. [VALUE]をクリックし、0.12と入力します。
  29. returnをクリックし、PARAMメニューを終了します。
  30. returnをクリックし、control cardsメニューを終了します。

ジョブのサブミット

  1. AnalysisページからOptiStructパネルをクリックします。
    図 10. OptiStructパネルへのアクセス

    OS_1000_13_17
  2. save asをクリックします。
  3. Save Asダイアログで、OptiStructモデルファイルを書き出す場所を指定し、ファイル名としてHalf_carと入力します。
    OptiStruct入力ファイルには、拡張子 .femが推奨されます。
  4. Saveをクリックします。
    入力ファイル欄には、Save Asダイアログで指定されたファイル名と場所が表示されます。
  5. export optionsのトグルをallにセットします。
  6. run optionsのトグルをanalysisにセットします。
  7. memory optionsのトグルはmemory defaultにセットします。
  8. OptiStructをクリックし、OptiStructジョブを開始します。
ジョブが成功した場合、Half_car.femが書き出されたディレクトリに新しい結果ファイルがあるはずです。何らかのエラーがある場合、Half_car.outファイルはデバッグを手助けするエラーメッセージを探すのに良い場所です。
そのディレクトリに書かれるデフォルトのファイルは：
Half_car.html
問題の定式と解析結果のサマリーに関する解析のHTMLレポート。
Half_car.out
ファイルの設定、最適化問題の設定、実行に必要なRAMおよびディスクスペースの推定量、各最適化反復計算の情報、解析時間等、特定の情報を含むOptiStructの出力ファイル。ワーニングおよびエラーに関しては、このファイルを確認すること。
Half_car.h3d
HyperViewバイナリ結果ファイル。
Half_car.res
HyperMeshバイナリ結果ファイル。
Half_car.stat
解析のプロセスの間のそれぞれのステップでのCPU情報を提供する、解析のプロセスの要約。

結果の確認

このステップでは変位結果（.mvwファイル）をHyperGraphで表示する方法を示します。HyperView結果ファイル（.h3d）は節点セット出力で指定された節点の変位結果のみを含んでいます。
  1. OptiStructパネルから、HyperViewをクリックします。
    HyperViewが起動され、結果が読み込まれます。HyperViewにモデルと結果が正しく読み込まれたことを示すメッセージウィンドウが現われます。
  2. 表示されたら、Closeをクリックし、メッセージウィンドウを閉じます。
  3. HyperViewウィンドウで、File > Open > Sessionをクリックします。
    Open Session Fileウィンドウが開きます。
  4. ジョブが実行されたディレクトリを選び、ファイルHalf_car_freq.mvwを選択します。
  5. Openをクリックします。
    旧データを廃棄するかどうかのワーニングが現われます。
  6. Yesをクリックします。
    ページ毎に２つのグラフで、合計１ページがHyperGraphに表示されます。グラフタイトルは、Subcase 1 (subcase 1) pressure at grid 18881 を示しています。
  7. AxisツールバーアイコンannotateAxes-24をクリックします。
  8. AxisがPrimaryとHorizontalにセットされていることを確認してください。
  9. Scale and Ticsタブをクリックします。
  10. トグルがLinearにセットされていることを確認してください。
    図 11.


  11. Axisで、HorizontalからVerticalに切り替えます。
  12. Scale and Tics (Magnitude)タブをクリックします。
  13. トグルがdB10にセットされていることを確認してください。
    図 12.


    このページには２セットの結果があります。上のグラフは、位相角vs周波数（log）、下のグラフは節点 18881の圧力のMagnitude verses Frequency (log)（下図参照）を示しています。
    図 13.


これでHyperGraphでの結果のポスト処理は終わりです。