OS-T:9000 OptiStructとVABSの統合

構造FEAソルバーであるOptiStructとVABSはAltair Simulationプラットフォームで統合されており、細長い構造物を解析するために、任意の断面形状と材料のビーム断面特性の完全なセットをアドホックな運動学的仮定なしに計算できる後者の機能を利用しています。

HyperMeshがVABSの入力を生成するための主要な前処理ツールであることに変わりはありません。

VABS(Variational Asymptotic Beam Sectional Analysis)は、1次元のビーム特性を計算し、細長い複合構造物(および等方性材料)の3次元応力 / ひずみを復元するための断面解析ツールです。

OptiStructおよびAltair Simulationフレームワークのこの機能強化により、ユーザーはOptiStructの単一のシームレスな実行で複合材ビームの解析を実行することができ、ここでVABSが内部的に呼び出されます。VABSライブラリはAltair Connectからダウンロードできます。
注: Altair Connectには、ダウンロード可能なVABSのバージョンが2つあり、 OptiStruct-VABS Integrationでは、選択するVABSのバージョンは“VABS-# Packaged with OptiStruct”という名前で、 # はVABSソフトウェアのバージョン番号です。
ダウンロードしたVABSライブラリは、以下の場所にあるOptiStructのインストールディレクトリに置いてください。
  • Windowsの場合:

    <install_directory>\hwsolvers\optistruct\lib\win64\VABS

  • Linuxの場合:

    <install_directory>/hwsolvers/optistruct/lib/linux64/VABS

HyperMeshインターフェースには、構造特性や慣性特性を含む断面特性を計算するための入力として、形状や材料のすべての詳細情報を含む断面の有限要素メッシュを生成するための特別なユーティリティがあります。VABS互換の入力ファイルは、所定の作業ディレクトリに保存されます。統合されたワークフローは、さらに、同じHyperMeshセッション内で残差モデル用のOptiStruct入力デックを生成することを可能にします。

ASSIGNコマンドを使用して、OptiStructはVABSの入力を特定し、VABS実行ファイルを起動して等価剛性マトリックスを解いて生成します。OptiStructはさらにVABSの出力を読み込み、完全なソルバーランを実行します。


図 1. ハイレベルのワークフロー
本チュートリアルでは、均一な板厚の複合材パイプを検討します。プライ順序と配向は:


図 2. 複合材のパイプ
以下のプライ配向(単位は度)を有する直交異方性材料:
  • E11 = 141.693
  • E22 = E33 = 9.79056 (GPa)
  • G12 = G23 = G31 = 5.99844e9 (GPa)
  • v 12 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamODamaaBa aaleaacaaIXaGaaGOmaaqabaaaaa@3894@ = v 23 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamODamaaBa aaleaacaaIXaGaaGOmaaqabaaaaa@3894@ = v 31 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamODamaaBa aaleaacaaIXaGaaGOmaaqabaaaaa@3894@ = 0.42
1D要素(CBEAM)は、ビーム軸がパイプのX軸に一致するようにパイプの中心に沿って理想化されています。より良い収束のために、合計25個のCBEAM要素が選択されています。


図 3. 1D(CBEAM)要素を複合材パイプのX軸中心に合わせて配置

HyperMeshの起動とプロファイルの設定

本演習では、図 3に示されるモデルが使用されます。2D(SHELL)メッシュおよび1D(CBEAM)要素が既に理想化されています。

  1. HyperMesh Desktopを起動します。
    User Profilesダイアログが現れます。
  2. ApplicationにEngineering Solutionsを選択します。
  3. Aerospace-OptiStructに切り替え、OKをクリックします。

モデルのオープン

  1. File > Open > Modelをクリックします。
  2. optistruct.zipファイルから自身の作業ディレクトリに保存したosVABS_example.hmファイルを開きます。モデルファイルへのアクセスをご参照ください。
  3. Openをクリックします。
  4. ModelブラウザでPropertiesフォルダーのReviewを右クリックし、プロパティの割り当てとプライレイアウトをします。

    2つのプロパティ、wall_top_bottomおよびwall_sidesを確認します。

  5. いずれかのプロパティをクリックします。
  6. エンティティエディターで、Number of Pliesを展開します。
  7. テーブルエディタをクリックして、プライ順序と配向を確認します。

モデルのセットアップ

VABSの断面プロパティの更新と*.datファイルの作成

与えられた断面の断面プロパティを作成し、VABS入力ファイルを作業ディレクトリにエクスポートするための特別なユーティリティが開かれます。

  1. Aerospace > Beams > Beams from linesをクリックします。
    Beam Toolダイアログが開きます。
  2. ActionでValueの下にUpdateをドロップダウンメニューから選択します。
  3. Entitiesに0 Elementsをクリックし、黄色いElementsボタンをクリックします。
    これにより、パネルから要素が選択できるようになります。
  4. パネルで、黄色いelemsボタンをクリックします。
  5. by idを選択し、グラフィックスウィンドウから任意の1D(CBEAM)要素を1つ選択します(演習問題では、要素ID: 5761を選択します)。
  6. Proceedをクリックします。
  7. Property/Sectionsで、Update Section/Propertyの横にあるチェックボックスをクリックします。
  8. SectionでValue欄をクリックし、ドロップダウンメニューからElementsを選択します。
    これにより、パネルから要素を選択できるようになるはずです (上記の手順と同様)。
  9. パネルからElements > By Collectorをクリックし、シェル要素を含む4つのコンポーネント(pipe_top、pipe_bottom、pipe_right、pipe_left)を選択します。
  10. Proceedをクリックします。
  11. Export VABSの下で、VABS入力ファイル(*.dat)がエクスポートされる場所を指定します。


    図 4. Beam Toolダイアログで設定を定義
  12. Applyをクリックします。
    これは、可視化のみを目的として、断面の平面にシェルメッシュを生成します。


    図 5. 断面におけるシェルメッシュの可視化
  13. Closeをクリックし、Beam Toolダイアログを閉じます。
    断面のグラフィックス(断面メッシュ)をクリアするためのユーザー警告ウィンドウが表示されます。
  14. Yesをクリックします。

コンポーネントの更新

  1. Modelブラウザで、Components フォルダーを展開します。
  2. beamコンポーネントを右クリックし、Assignを選択します。
    Assign to Componentsダイアログが開きます。
  3. 新たに作成されたプロパティBeam_ID_Vabs_5761_Propを選択します。
  4. ApplyOKをクリックします。


    図 6. ビームコンポーネントへのプロパティの割り当て
    注:
    • この演習では板厚が均一なので、1つのCBEAM要素(エレメントID:5761)だけに断面プロパティを作成し、新たに作成したビームプロパティを残りの24個のCBEAM要素に割り当てました(Step 2)。
    • パイプの断面に変更がある場合は、複数の断面を作成する必要があります。
    • Beam Toolは、Step 3で複数のCBEAM要素が選択された場合、選択されたCBEAM要素ごとに個別の断面積を生成し、そのため個別のVABS入力ファイル(*.dat)が生成されます。
  5. Modelブラウザで、Componentsフォルダーを展開します。
  6. シェル要素を含む4つのコンポーネント((pipe_toppipe_bottompipe_rightpipe_left)を選択します。
  7. 右クリックし、Deleteを選択します。

材料の定義

MAT8エントリでは、E2NU12、およびHyperMeshはデフォルトで定義されており、現時点でE3NU13およびNU23を0として*.datファイルにエクスポートします。したがって、*.datファイルに以下のデータを手動で設定します。

  1. テキストエディタで、Step 4で保存した *.dat ファイルを開きます。
  2. .datファイルの最後までスクロールして次の箇所を探します:
    1.41963000e+11  9.79056000e+09  0.00000000e+00
    5.99844000e+09  5.99844000e+09  5.99844000e+09
    4.20000000e-01  0.00000000e+00  0.00000000e+00
  3. 次のように読めるよう、手動で編集します:
    1.41963000e+11  9.79056000e+09  9.79056000e+09
    5.99844000e+09  5.99844000e+09  5.99844000e+09
    4.20000000e-01  4.20000000e-01  4.20000000e-01
    注: 上記のようなE3、NU13、NU23のデフォルト設定は、HyperMeshの将来のバージョンで自動化される予定です .

断面プロパティとASSIGNカードの確認

  1. Modelブラウザで、Propertiesフォルダーを展開します。
  2. Beam_ID_Vabs_5761_Propを右クリックします。
  3. Card Editをクリックし、PBEAMLのGroupとTypeを確認します。
    1. GroupがVABSにセットされていることを確認します。
    2. TypeはVabs5761(断面積生成に選択されたCBEAMの要素IDの接尾部)に設定されます。


      図 7. PBEAML VABSプロパティの確認
  4. Modelブラウザで、Cardsフォルダーを展開します。
  5. ASSIGNを右クリックします。
  6. エンティティエディターで、以下のことを確認すます:
    1. TypeをVABSに設定します。
    2. Section_NameはVabs5761に設定します。
    3. Section_Pathは、VABS入力*.datがエクスポートされるディレクトリに設定されます。


      図 8. ASSIGNエントリの確認

荷重コレクターの作成

ここでは、モーダル解析の境界条件が定義されます。

  1. Modelブラウザ内で右クリックし、Create > Load Collectorを選択します。
  2. NameにSPCと入力します。
  3. Card Imageには、NONEを選択します。
  4. ツールバーからBCs > Create > Constraintsを選択します。
  5. 節点ID: 1を選択し、1から6までのすべての自由度にチェックを入れます。
  6. createをクリックします。
    CBEAM要素の節点1にSPCが作成されます。
  7. をクリックし、制約条件作成のパネルを終了します。
  8. Modelブラウザ内で右クリックし、Create > Load Collectorを選択します。
  9. NameにEigenと入力します。
  10. Card Imageには、EIGRLを選択します。
  11. NDに10と入力します。

荷重ステップの作成

このステップでは、モーダル荷重ケースを定義します。

  1. Modelブラウザ内で右クリックし、Create > Load Stepを選択します。
  2. NameにModeShapesと入力します。
  3. Analysis typeには、Normal modesを選択します。
  4. SPCに、荷重コレクターのリストからSPCを選択します。
  5. Method(STRUCT)に、Eigenを選択します。
    モーダル解析の荷重ケースが作成されます。

ジョブのサブミット

  1. AnalysisページからOptiStructをクリックします。
  2. input file欄に続くsave asをクリックします。
    Save Asダイアログが開きます。
  3. File nameにosVabs_example.femと入力します。
  4. input file:欄は、osVabs_example.femの場所に設定されます。
  5. Saveをクリックします。
  6. OptiStructをクリックし、解析をサブミットします。

結果の確認

OptiStructは自動的にVABS実行ファイルを起動し、*.datファイルを実行して等価剛性マトリックスを生成します。OptiStructはVABSの出力を読み込み、完全なソルバーランを実行します。

  1. 解析プロセスが終了したら、HyperViewをクリックして、結果を読み込みます。
  2. Resultsタブで、サブケースフィールドからSubcase 1 (ModeShapes) を選択します。
  3. HyperViewContour panelに進みます。
  4. Result typeに、Eigen Mode(v)を選択します。
  5. Applyをクリックし、コンター図を表示します。
    Subcaseの下のモード番号を変更することで、異なるモードをコンター化することができます。


    図 9. HyperViewでのモード形状結果

    等方性のビームとは異なり、複合材のビームは様々な変形の間に強い連成を示し、このケースでは、伸長-ねじれ連成と曲げ-せん断連成がモード形状から明らかです。