OS-T:1300 直接法による平板の周波数応答解析

本チュートリアルでは、既存の平板のFEモデルを読み込み、境界条件を与え、有限要素解析を実行する方法について実際に示します。

開始する前に、このチュートリアルで使用するファイルを作業ディレクトリにコピーします。

平板は、直接法での周波数が変化する単位荷重の加振を受けます。ポスト処理では、HyperViewHyperGraphでモードシェイプの応答、および周波数-位相の出力特性の可視化が行われます。

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

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

モデルの読み込み

  1. File > Import > Solver Deckをクリックします。
    Importタブがタブメニューに追加されます。
  2. File typeにOptiStructを選択します。
  3. Filesアイコンfiles_panelを選択します。
    Select OptiStruct Fileブラウザが開きます。
  4. 自身の作業ディレクトリに保存したdirect_response_flat_plate_input.femファイルを選択します。
  5. Openをクリックします。
  6. Import、続いてCloseをクリックし、Importタブを閉じます。

モデルのセットアップ

荷重と境界条件の適用

続くステップで、モデルの1辺が拘束されます。鉛直の単位荷重が、板の自由端の角にz-方向の正の向き(上向き)に作用します。
  1. Modelタブをクリックします。
  2. Model Browser内で右クリックし、Create > Load Collectorを選択します。
  3. Nameにspcsと入力します。
  4. Colorをクリックし、カラーパレットから色を選択します。
  5. Card ImageをNoneに設定します。
    新しい荷重コレクターspcsが生成されます。
  6. Model Browser内で右クリックし、Create > Load Collectorを選択します。
  7. Nameにunit-loadと入力します。
  8. Colorをクリックし、カラーパレットから別の色を選択します。
    新しい荷重コレクターunit-loadが生成されます。

拘束の作成

  1. Model Browserで、Load Collectorを展開してspcs > Make Currentをクリックします。
    1.

    rd2000_spc_load_collector
  2. Display NumbersアイコンinfoNumbers-24をクリックします。
  3. nodes > displayedをクリックします。
  4. on(緑のボタン)を選択します。
    平板上の全ての節点番号が今度は表示されるはずです。
  5. returnをクリックし、メインメニューに戻ります。
  6. BCs > Create > Constraintsをクリックし、Constraintsメニューを開きます。
  7. エンティティ選択スイッチをクリックし、ポップアップメニューからnodesを選択します。
  8. nodesをクリックし、節点 5、29、30、31、32 を選択します(図 2)。
    2. 単点拘束の付与に選択される節点

    os1100_pic1
  9. dof6のみチェックを外し、dof1dof2dof3dof4、およびdof5を拘束します。
    • チェックマークのついている自由度は拘束され、ついていないものはフリーとなります。
    • Dofs 1、2および3は、x、y、z方向の並進自由度を表します。
    • Dof 4、5、および6は、x、y、z方向の回転自由度を表します。
  10. createをクリックします。
    選択された節点は dof6 がチェックされていないことから、z-軸周りの回転のみ自由となります。
  11. returnをクリックし、メインメニューに戻ります。

Create a Unit Load at a Point on the Flat Plate

  1. In the Model Browser, right-click on the load collector unit-load and select Make Current.
  2. From the Analysis page, click load types.
  3. Select constraint = and select DAREA from the extended entity selection menu.
  4. Click return to exit the Load Types panel.
  5. Click BCs > Create > Constraints to open the Constraints menu.
  6. Select node number 19 on the plate by clicking on it (図 3).
    3. Node Selected for Creating Unit Vertical Load

    os1100_pic2
  7. Uncheck all the dof's except dof3 and click the = to the right of dof3 and enter a value of 20.
  8. Click load types= and verify that DAREA is selected from the extended entity selection menu.
  9. Click create, and then click return.
    The unit load is applied to the selected node.

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) = 1000.0, y(2) = 1.0.
  4. Close the Curve Editor window.
  5. From Curves, select tabled1.
  6. For Type, select TABLED1 from the drop-down menu.
    This provides a frequency range of 0.0 to 1000.0 with a constant 1.0 over this range.

Create a Frequency Dependent Dynamic Load

  1. In the Model Browser, right-click and select Create > Load Step Inputs.
  2. For Name, enter rload2.
  3. For Config type, select Dynamic Load – Frequency Dependent from the drop-down list.
  4. For Type, and select RLOAD2 from the drop-down list.
  5. For Excited, click Unspecified > Loadcol.
  6. In the Select Loadcol dialog, select unit-load from the list of load collectors and click OK to complete the selection.
  7. For TB, select the tabled1 curve.
    The type of excitation can be an applied load (force or moment), an enforced displacement, velocity or acceleration. The field Type in the RLOAD2 load step input defines the type of load. The type is set to applied load by default.

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 20.0.
    2. For DF, enter 20.0.
    3. For NDF, enter 49.
  7. Click Close.
    This provides a set of frequencies beginning with 20.0, incremented by 20.0 and 49 frequencies increments.

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 (direct) from the drop-down menu.
  4. For SPC, select Unspecified > Loadcol.
  5. From the Select Loadcol dialog, select SPCS.
  6. For DLOAD, select rload2 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 in the load collector spc and the unit load in the load collector step input rload2 with a set of frequencies defined in load collector freq1

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 None.
  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 with IDs 15, 17 and 19.
  7. Click proceed.

出力セットとマスファクターの生成

  1. Setup > Create > Control Cardsをクリックし、Control Cardsパネルを開きます。
  2. GLOBAL_OUTPUT_REQUESTを選択し、DISPLACEMENTの横のボックスにチェックを入れます。
  3. FORM(1)の下で、ポップアップメニューからPHASEを選択します。
  4. OPTION(1)の下で、ポップアップメニューからSIDを選択します。
    黄色い新しい欄が現われます。
  5. SID(1)ボックスをダブルクリックし、SETAを選択します。
    値1がSID入力ボックスの下に現れます。これは、set 1内の節点に限った出力を設定しています。
    4.

    rd2000_displacement
  6. returnをクリックし、GLOBAL_OUTPUT_REQUESTSメニューを終了します。
  7. Control CardsパネルからFORMATを選択します。
    作業領域のスクリーンに新しいウィンドウが現われます。
  8. number_of_formats =をクリックし、値2を入力します。
  9. 作業領域の拡張メニューで、1つ目のFORMAT_V1欄をクリックし、ポップアップメニューからOPTIを選択します。
    OPTIの使用で、 .disp.strsなどの OptiStruct ASCII結果ファイルが計算終了時に生成されます。これらのファイルはポスト処理で用いられます。
  10. 2つ目のフィールドボックスがH3Dにセットされていることを確認してください。
  11. returnをクリックしてFormatメニューから抜け、Control Cardsメニューに戻ります。
  12. nextをクリックしてPARAMサブパネルを選択します。
  13. 左隅の矢印を用いてリストをスクロールダウンし、COUPMASSの隣のボックスにチェックを入れます。
    新しいPARAMカードが作業領域のスクリーンに現われます。
  14. COUPM_V1の下のNOをクリックし、ポップアップメニュー選択からYESを選択します。
    YESの選択で整合質量マトリックスアプローチが用いられます。
  15. 左隅の矢印を用いてリストをスクロールダウンし、Gの横のボックスにチェックを入れます。
    新しいPARAMカードが作業領域のスクリーンに現われます。
  16. G_V1の下をクリックし、入力ボックスに値0.06を入力します。
    この値は一様な構造減衰の係数で、臨界減衰[C/C 0]比に2.0を掛けることにより得られます。
  17. 左隅の矢印を用いてリストをスクロールダウンし、WTMASSの横のボックスにチェックを入れます。
    作業領域のスクリーンに新しいウィンドウが現われます。
  18. WTM_V1の下をクリックし、入力ボックスに値0.00259を入力します。
    作業領域にポップアップメニューが現われ、3つのPARAMステートメントが表示されます。この係数は全ての質量入力が重量単位の場合に入力します。このPARAMで、全ての質量マトリックスの項にこの係数が掛けられます。
    5.

    rd2000_param
  19. returnをクリックし、PARAMメニューを終了します。
  20. OUTPUTサブパネルを選択します。
  21. KEYWORDがHGFREQにセットされていることを確認します。
    HGFREQの使用により、周波数がHyperGraph形式で出力されます。
  22. FREQの下のボックスをクリックし、ポップアップからALLを選択して、すべての周波数についてすべての出力リクエストを選びます。
  23. number_of_outputsは1にセットしたままにします。
  24. returnをクリックし、OUTPUTから抜けます。
  25. returnをクリックし、Control Cardsパネルを終了します。

ジョブのサブミット

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

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

結果の表示

このステップでは変位結果(.mvw ファイル)をHyperGraphで表示する方法を示し、この実行での変位出力(.dispファイル)についても説明します。

HyperView結果(.h3d)は、節点セット出力で指定された3節点の変位結果のみを含んでいます。
  1. OptiStructパネルから、HyperViewをクリックします。
    HyperViewが起動され、結果が読み込まれます。HyperViewにモデルと結果が正しく読み込まれたことを示すメッセージウィンドウが現われます。
  2. 表示されたら、Closeをクリックし、メッセージウィンドウを閉じます。
  3. HyperViewウィンドウで、File > Open > Sessionをクリックします。
    Open Session Fileウィンドウが表示されます。
  4. ジョブが実行されたディレクトリを選び、ファイルflat_plate_direct_response_freq.mvwを選択します。
  5. Openをクリックします。
    既にある内容を廃棄するかどうかを問うワーニングが表示されます。
  6. Yesをクリックします。
    ページ毎に2つのグラフで、合計3ページが表示されます。ページ 1 のグラフタイトルは サブケース 1 の節点 15 の変位であることを示しています。

    このページには2セットの結果があります。上のグラフは、位相角vs周波数(log)、下のグラフは、grid 15における変位の絶対値vs周波数(log)です(図 7参照)。

    7. 節点15の周波数応答

    rd2000_results
  7. Next PageアイコンpageNext-24をクリックします。
    これで、Subcase 1 (subcase1) - Displacement of grid 17を示す、ページ 2が表示されます(図 8)。
    8. 節点17の周波数応答

    rd2000_results2
  8. 再度Next PageアイコンpageNext-24を選択し、Subcase 1 (subcase1) - Displacement of grid 19を含むページ 3を表示させます(図 9)。
    9. 節点19の周波数応答

    rd2000_results3
    これでHyperGraphでの結果のポスト処理は終わりです。
  9. 変位ファイル(.disp)をテキストエディタを用いて開きます。
    第2行の最初のフィールドは反復番号を示し、第2のフィールドはデータの点数を、第3のフィールドはその反復の周波数を示しています。

    第3行の最初のフィールドは節点番号、次に x、y、z 方向変位の大きさ、x、y、z 方向の回転変位の大きさを示しています。

    第4行の最初のフィールドは節点番号、次に x、y、z 方向変位の位相角、x、y、z 方向回転角を示しています。