OS-HM-T:6020 E-N（ひずみ - 寿命）法を用いたFatigue Process Manager（FPM）

チュートリアルレベル：中級与えられた周期荷重の下で塑性ひずみが発生する場合、疲労寿命の予測にE-N（ひずみ - 寿命）法が選択されるべきです。S-N（応力 - 寿命）法は疲労挙動に塑性ひずみが中心的な役割を担う低サイクル疲労には適していません。

開始する前に、このチュートリアルで使用するファイルを作業ディレクトリにコピーします。
S-N解析で疲労寿命が 10,000 サイクル以下となった場合は、E-N法がより良い選択となり得る兆候です。E-N法は計算的にはS-N法よりも高価ですが、高サイクル疲労の場合においても妥当な推定を行うはずです。
図 1. S-N曲線の低サイクルと高サイクルの領域

rd2070_SN_curve
E-N理論は単軸ひずみを扱うため、それぞれの計算点、それぞれの時間ステップにおけるひずみ成分は一つの統合された値にする必要があり、そこで、等価な公称ひずみが用いられ、E-N曲線に適用されます（図 2）。
図 2. ひずみ-寿命曲線

rd2080a_strainlife_curve

OptiStructではさまざまなひずみの組み合わせタイプが利用可能で、デフォルトは “Absolute maximum principle strain （絶対値最大の主ひずみ）”です。一般的に脆性材料には "Absolute maximum principle stain"が推奨されますが、延性材料には "Signed von Mises stain"が推奨されます。符号付きパラメータの符号には最大絶対値主値の符号が用いられます。

本チュートリアルでは、E-N法による疲労寿命を評価することができます。

コントロールアームは、図 3に示すとおり、ブレーキ力と鉛直力を受けています。２つの荷重の時刻歴は 1 Hzで 2545 秒間図 4と図 5に示すように得られており、これが使われます。コントロールアームに用いられる材料はアルミニウムでそのE-N曲線を図 6に示します。亀裂は常に表面から始まるため、 表皮のシェルメッシュがソリッド要素を覆って置かれています。これにより、計算精度も良くなります。
図 3. 疲労解析に用いるコントロールアームモデル

rd2070a_control_arm
図 4. 鉛直力の荷重時刻歴

rd2070a_load_time
図 5. ブレーキ力の荷重時刻歴

rd2070a_vertical_force
図 6. アルミニウムのE-N曲線

rd2080a_en_curve

Altair HyperWorks / HyperMeshとProcess Managerの起動

  1. 起動メニューからAltair HyperWorksを起動します。
    New Sessionダイアログが開きます。
  2. HyperMeshラジオボタンを選択し、プロファイルをOptiStructに設定し、Create Sessionボタンをクリックします。
  3. Templatesリボンから、Analyzeメニューを選択し、Fatigue PMを選択します。
  4. New Session Nameに<my_session_name>と入力します。
  5. Working Folderに自身の作業フォルダーを選択します。
  6. Createをクリックします。
    これで、現在ロードされている疲労プロセステンプレートの内容を保存するための新しいファイルが生成されます。

Import the Model

  1. Make sure the task Import File is selected in the Fatigue Analysis tree.
  2. For the Model file type, select OptiStruct.
  3. Click the Open model file icon files_panel.
    A Select File browser window opens.
  4. Select the ctrlarm.fem file you saved to your working directory and click Open.
  5. Click Import.
    This loads the control arm model. It includes a whole definition of two static subcases, elements sets, and material static properties, etc.
  6. Click Apply.
    This guides you to the next task Fatigue Subcase of the Fatigue Analysis tree.
    図 7. Import a Finite Element Model file


モデルのセットアップ

Create a Fatigue Subcase

  1. Make sure the task Fatigue Subcase is selected in the Fatigue Analysis tree.
  2. In the Create new fatigue subcase field, enter fatsub_fpmtut.
  3. Click Create.
  4. For the Select existing fatigue subcase field, select the newly created fatigue subcase fatsub_fpmtut.
    fatsub_fpmtut is selected as the active fatigue subcase. Definitions in the following processes (analysis parameters, fatigue elements and properties, loading sequences, etc.) will be for this subcase.
  5. Optionally, you can choose to create all fatigue solver cards (such as FATPARM, FATDEF, FATEVNT, PFAT etc. that are created in subsequent steps) in a separate include file. For this, you should select the check box Create include file for fatigue solver cards (optional).
  6. Click Apply.
    This saves the current definitions and guides you to the next task Analysis Parameters of the Fatigue Analysis tree.
    図 8. Create and Select Active Fatigue Subcase to Process


Apply Fatigue Analysis Parameters

  1. Make sure the task Analysis Parameters is selected in the Fatigue Analysis tree.
  2. Select the following options:
    Analysis type
    E-N
    Stress combination method
    Signed von Mises
    FEA model unit
    MPA
    Mean stress correction
    SWT
    Rainflow type
    STRESS
    Plasticity correction
    NEUBER
  3. Enter the following values:
    Certainty of survival
    0.5
    Gate
    0.0
  4. Click Apply.
    This saves the current definitions and guides you to the next task Elements and Materials of the Fatigue Analysis tree. For details, consult the Altair HyperWorks 2025 help.
    図 9. Fatigue Analysis Parameters Definition


疲労要素と材料の追加

  1. Fatigue AnalysisツリーでタスクElements and Materialsが選択されていることを確認します。
  2. Add Materialをクリックします。
    Material Dataウィンドウが開きます。
  3. Material名にAluminumを選択します。
  4. Stress unitがMPAにセットされていることを確認します。
  5. Ultimate tensile strength (UTS)に、600と入力します。
  6. EN曲線を定義するInput methodにEstimate From UTSを選択します。
    図 10. 材料データ定義


  7. Show EN curve definitionアイコン05_sn_curveをクリックします。
    EN 材料パラメータの生成方法を紹介するEN method descriptionウィンドウが開きます。
    図 11. ENカーブ材料定義


  8. Closeをクリックします。
  9. Material TypeにAluminum and Titanium Alloysを選択し、Estimateをクリックします。
    EN 曲線定義の全てのデータが自動的に推定されます。
  10. ウィンドウの下部にあるPlot EN Curveをクリックし、EN曲線を表示させます。
  11. CloseでEN Curve plotウィンドウを閉じます。
  12. Saveをクリックし、選択された要素へのENデータ定義を保存します。
    図 12.


  13. Add Propertyをクリックします。
    Property Dataウィンドウが開きます。
  14. Property TypeにProperty - PSHELLを選択します。
  15. Property nameにshellを選択します。
    これはソリッドのコントロールアームを覆う表皮です。
  16. CreateをクリックしてPFATプロパティを作成します。
  17. LAYERにTOPを選択します。
  18. FINISHにNONEを選択します。
  19. TREATMENTにNONEを選択します。
  20. KFの後ろの欄（Fatigue strength reduction factor）は1.0のままにしておきます。
  21. Closeをクリックし、選択された要素へのENデータ定義を保存します。
  22. Applyをクリックします。
    これで現在の定義が保存され、Fatigue Analysisツリーの次のタスクLoad-Time Historyに導かれます。
    図 13. 特性データ定義

    rd2080b_mat_data
    図 14. 要素と材料の定義


荷重時刻歴の適用

  1. Fatigue AnalysisツリーでタスクLoad-Time Historyが選択されていることを確認します。
  2. Add by Fileをクリックします。
    Load Time Historyウィンドウが表示されます。
  3. Load-time history名にLTH1と入力します。
  4. Load-time historyタイプにCSVを選択します。
  5. Open load-time fileアイコンfiles_panelをクリックします。
    Open fileブラウザウィンドウが開きます。
  6. load1.csvを探して選択します。
  7. Open > Importをクリックします。
  8. Saveをクリックし、新しい荷重-時間履歴をHyperMeshデータベースに書き込みます。
  9. もう１つの荷重-時間履歴LTH2を load2.csvを読み込んで生成します。
    図 15. 荷重-時間履歴の読み込み


  10. Plot L-Tをクリックし、荷重-時間履歴を表示させます。
  11. Load Time Historyウィンドウを閉じます。
  12. Applyをクリックします。
    これで現在の定義が保存され、Fatigue Analysisツリーの次のタスクLoading Sequencesに導かれます。
    図 16. 荷重時刻歴の定義


荷重順序

このステップでは、2つの荷重時刻歴の1つのイベントが生成されます。言い換えれば、2つの荷重時刻歴で生じる応力の線形の重ね合わせが、解析の間に要求されたことになります。このイベントを用いて、1つの荷重順序が構成されます。
  1. Fatigue AnalysisツリーでタスクLoading Sequencesが選択されていることを確認します。
  2. Addをクリックします。
    Loading Mappingウィンドウが開きます。
  3. ChannelsにLTH1とLTH2を選択します。
  4. SubcaseにSUBCASE1とSUBCASE2を選択します。
  5. ラジオボタンAutoを有効にし、イベントの作成法はデフォルトのSingle Eventのままとします。
  6. +をクリックし、2つのサブケース、2つのチャンネルのシングルイベントを作成します。
  7. 下図のように、Scaleを 5.0に設定します。
    図 17. 荷重-時刻歴を静的サブケースと関連付けするための荷重マッピング


  8. Saveをクリックしてウィンドウを閉じ、選択されたサブケースとチャンネルを用いて疲労イベントを作成します。
    図 18. 荷重順序の定義


ジョブのサブミット

  1. Fatigue AnalysisツリーでタスクSubmit Analysisが選択されていることを確認します。
  2. Save .femファイルアイコンfiles_panelをクリックします。
    Save Asブラウザウィンドウが開きます。
  3. ファイルを保存するディレクトリをセットし、File nameにctrlarm_fpmtut.femと入力します。
  4. Saveをクリックし、ウィンドウを閉じます。
  5. Saveをクリックし、OptiStructモデルファイルを保存します。
  6. Run Optionにanalysisを選択します。
  7. Submitをクリックします。
    これでOptiStruct 2025が起動し、疲労解析が実行されます。ジョブが成功した場合、ctrlarm_fpmtut.femに選択されたディレクトリに新しい結果ファイルが現れます。
    そのディレクトリに書かれるデフォルトのファイルは：
    ctrlarm_fpmtut.0.3.fat 反復ステップでのそれぞれの疲労サブケースの疲労解析結果を含むASCIIフォーマットファイル。
    ctrlarm_fpmtut.h3d 静解析の結果と疲労解析の結果の両方の Hyper 3D バイナリ結果ファイル。
    ctrlarm_fpmtut.out ファイルのセットアップ、疲労問題のセットアップ、計算時間の情報などを含むOptiStruct出力ファイル。ワーニングやエラーもこのファイルで確認します。
    ctrlarm_fpmtut.stat 解析のプロセスの間のそれぞれのステップでのCPU情報を提供する、解析のプロセスの要約。
    注: filename.#.fatは、疲労最適化が実行される場合に限り、最初と最後の反復計算において各疲労サブケースについて作成されます。
    図 19. 疲労解析のサブミット


結果のポスト処理

  1. Fatigue AnalysisツリーでタスクPost-processingが選択されていることを確認します。
    前のサブミットの後、疲労解析が正常に終了すると、自動的にこのタスクに入ります。
  2. FatigueサブケースにSelect Subcaseが選択されていることを確実にします。
  3. Result TypeとData Componentに、コンター表示したいデータをドロップダウンメニューから選択します。
  4. Load H3D Results (HV)をクリックします。
    これでHyperViewが起動し、ctrlarm_fpmtut.h3d結果ファイルが読み込まれます。これは、選択された結果タイプとコンポーネントについて、結果のコンターを適用します。より詳細な結果については、HyperViewを使用することができます。
  5. Exitをクリックし、Fatigue Process Managerをアンロードします。
    図 20. ポスト処理


    図 21. HyperViewでのダメージコンター