動的空力弾性解析

動的空力弾性解析は、力と加速度が時間に依存する、空力荷重がかかる柔構造でのたわみを調査するものです。

フラッター解析

空力弾性フラッターは、空力荷重、弾性荷重、および慣性荷重の相互作用に関連付けられた構造の動的不安定性です。

空力弾性システムのフラッター解析では、システムの速度（およびマッハ数）と、システムがフラッターの状態に達する振動周波数の特定が必要です。この現象では、弾性体にかかる空力荷重が、その固有振動モードと連成し、振幅が増加する振動運動を生じさせます。

これが、壊滅的な構造破壊につながる可能性があります。したがって、空力荷重がかかる構造はフラッターを避けるように慎重に設計する必要があります。

有限要素解析では、フラッターの予測で一連の複素固有値解析が行われます。OptiStructは、選択された周波数範囲内の構造振動モードが自由度として使用されるモーダルアプローチを使用します。

フラッターは亜音速（マッハ数 < 1.0、亜音速領域の速度）および超音速（マッハ数 > 1.0、超音速領域の速度）の両方でサポートされています。現時点ではTransonic regime (マッハ数 = 1.0)はサポートされていません。

手法

OptiStructでサポートされるフラッター解析の手法には、K、KE、PK、PKNLの4つの種類があります。

後述するこれらの手法のそれぞれで、特定のケースについて複素固有値（ $p$ ）が抽出されます。これにより、空気速度と減衰のペア（ $V$ - $p$ ）と、空気速度と周波数のペア（ $V$ - $f$ ）を特定できます。

固有値の正確な形式は、手法によって異なります。

K法とKE法では、固有値が次の式で与えられます：

p^{2} = \frac{- V^{2}}{1 + i g}

PK法とPKNL法では、固有値が次の式で与えられます：

p = ω (γ \pm i)

ここで、

$γ$: 過渡減衰速度係数。
$ω$: 円振動数。 $2 π f$ で与えられます。

K法：
- K法の背景にある考え方は、人工的な減衰を伴う（構造減衰項の形で）システムを導入することで、システムをフラッターポイントに持っていくというものです。
- FLFACTエントリからの密度、マッハ数、および換算振動数のすべての組み合わせで複素固有値のセットが出力されます。
- そのため、K法から得られる結果は、解釈がやや難しくなります。
KE法：
- KE法は、K法の派生手法で、次のような違いがあります。
  - 粘性減衰項は無視されます。
  - 複素モードは出力されません。
  - フラッター出力は、モードに基づいて配置され、固有値外挿手法を使用してソートされます。
- このような特徴から、KE法は、K法の計算コストが低く、解釈しやすいバージョンだと言えます。
PK法：
- PK法では、反復プロセスを使用したダブレット格子法による空気力学（単振動を前提とする）を使用して一般的なフラッター解析を実行することができます。
- この手法では、剛性マトリックスに対する虚数の寄与が無視されます。これは、構造減衰項とPARAM, KDAMP, -1によるモード減衰が考慮されないことを意味します。
- 固有値抽出は、FLFACTエントリからの密度、マッハ数、および速度のすべての組み合わせに対して実行されます。
- 換算振動数の初期推測（ $k$ ）は、出力が更新された換算振動数を返す複素固有値問題の解析（収束まで繰り返されるプロセス）に使用されます。モードは、左と右の複素固有ベクトルを使用してさまざまな空気速度にわたって追跡されます。
PKNL法：
- PKNL法は、‘ループなし’で設計されたPK法のバリアントです。
- この手法では、密度比、マッハ数、および速度に関するFLFACTデータ内のエントリの数を同じにする必要があり、固有値抽出が密度比、マッハ数、および速度の線形選択ごとに実行されます。

以下のシナリオで、PK法とPKNL法の違いを簡単に説明します：

シナリオ1：

フラッター解析が以下のポイントで実行されます：(0.5, 0.3, 100), (0.5, 0.3, 200), (1.0, 0.3, 100), (1.0, 0.3, 200), (0.5, 0.4, 100), (0.5, 0.4, 200), (1.0, 0.4, 100) and (1.0, 0.4, 200).

$--1---><---2--><---3--><---4--><--5---><--6---><---7--><--8---><---9-->
FLUTTER      103      PK       1       2       3       L       4
FLFACT         1     0.5     1.0
FLFACT         2     0.3     0.4
FLFACT         3   100.0   200.0

シナリオ2：

フラッター解析が以下のポイントで実行されます：(0.5, 0.3, 100) and (1.0, 0.4, 200).

$--1---><---2--><---3--><---4--><--5---><--6---><---7--><--8---><---9-->
FLUTTER      103    PKNL       1       2       3       L       4
FLFACT         1     0.5     1.0
FLFACT         2     0.3     0.4
FLFACT         3   100.0   200.0

入力

図 1はOptiStructでの空力弾性フラッター解析のワークフローを要約します。

下の表に、空力弾性フラッター解析での関連する入力ファイルエントリを示します。

表 1. バルクデータエントリ
エントリ	説明
AERO	飛行条件を定義します。
MKAERO1/MKAERO2	空力マトリックスの陽的計算用のマッハ数と換算振動数のペアを指定します。
FLFACT	フラッター解析のフラッターパラメータの値（密度比、速度、および換算振動数）を指定します。
FLUTTER	フラッター解析の手法（K/KE/PK/PKNL）とパラメータを選択します。このエントリは、FLFACTの定義も参照します。
EIGC	K法用の複素固有値手法を選択します。
EIGRL/EIGRA	構造モードの計算方法とそれらの数を選択します。構造モードの数は、`LMODES`/`LFREQ`/`HFREQ`を使用して変更できます。 FLUTTERエントリ内の`NVALUE`は、.fltファイルに出力される固有値の数を制限するために使用できます。
PARAM, VREF	出力速度をスケーリングするために使用されます：Vout = V/Vref。
DMI	実マトリックスデータブロックを定義します。

表 2. 入出力オプションおよびサブケース情報エントリ
エントリ	説明
FMETHOD	FLUTTERエントリを参照します。
CMETHOD	複素固有値抽出では、EIGCエントリを参照します（K法のみ）。

問題の設定

次の入力ファイルの内容は、典型的な空力弾性フラッター解析セットアップの例を示しています。

$ ************************************************************
$ SNIPPET OF AN INPUT FILE FOR AEROELASTIC FLUTTER ANALYSIS
$ ************************************************************
DISP = ALL
SUBCASE 101
   SPC = 101
   METHOD = 102 
   FMETHOD = 103
BEGIN BULK
$--1---><--2---><--3---><--4---><--5---><--6---><--7---><--8---><--9---><--10-->
EIGRL   1       0.0     100.    4                               MASS
FLUTTER 103     PK      1       2       3  

FLFACT  1       0.4
FLFACT  2       0.5
FLFACT  3       150     -175      -200
$ If there are certain velocities for which flutter eigenvector output is desired,
$ in PK/PKNL methods, a negative sign can be added as shown above. In this example, 
$ output to the result file will be available only for the velocities: 175 and 200. 
$ See Output Section for more details.
AERO            102890. 2200.   .123E-110       0
MKAERO1 0.5                                                              
+       .001    0.002   0.005   0.01    0.05    0.1     0.2     0.3

MKAERO1 0.5                                                             
        0.5     0.6     0.8     1.0     1.5     2.0     3.0     4.0
$ other aeroelastic and structural entries

出力

フラッター解析サマリー

複素固有値解析の根が.fltファイルに出力されます。.fltファイル内の典型的なフラッター解析のサマリーを以下に示します。

K法およびKE法の場合：

固有値（ $p$ は、式 1を使用して空気速度（ $V$ ）と減衰（ $g$ を求めるために使用されます。
注: この方程式は、実数成分用と虚数成分用の2つの方程式に分割して、2つの未知数を求めることができます。
$p$ の値が、 $V$ が物理的である（実数）場合は、速度、減衰、および周波数の計算値が出力されます。
$p$ の値が、 $V$ が物理的でない（複素数）場合は、速度、減衰、および周波数の値の安定性情報のみが出力されます。
このようなケースでは、

g < 0

の場合にシステムが不安定になり、

g > 0

の場合にシステムが安定します。
$V$ の決定値とユーザー指定の換算振動数（ $k$ ）および基準翼弦長（AEROエントリからのREFC）を組み合わせて、周波数（ $f$ ）を計算するために使用できます。
$k = \frac{ω R E F C}{2 V}$
式 3から計算された $ω$ を使用して、
$f = \frac{ω}{2 π}$
KE法はフラッターサマリーをモード別に配置するため、減衰の符号が切り替わるタイミングを目視で確認することで、フラッターポイントを簡単に特定できます。

PK法およびPKNL法の場合：

周波数（ $f$ ）と減衰（ $g$ ）は、式 2を使用した固有値の形式から抽出されます。
注: この方程式は、実数成分用と虚数成分用の2つの方程式に分割して、2つの未知数を求めることができます。
報告される減衰は次のとおりです：
$g = 2 γ$
実根のみの場合は、減衰が次のようになります：
$g = \frac{2 p R E F C}{\ln (2) V}$

.fltファイルを、HyperMesh“Flutter Curves”機能に読み込むことで外部のポスト処理ツールを使用しなくても、（

V

g

）および（

V

f

）曲線を簡単に生成することができます。

H3D出力

現時点で、フラッター解析では次の出力要求のみがサポートされます。サポートされる出力フォーマットについては、空力弾性フラッター解析をご参照ください。

DISPLACEMENT
SDISPLACEMENT - K法、PK法、およびPKNL法の場合にのみ使用可能

次の表に、さまざまな手法の固有ベクトル出力に関する情報をまとめます。


方法	コメント
K法	固有ベクトル出力をポイントの形式で使用できます。各ポイントが1つの換算周波数に対応します。特定のポイントで、各モードの結果を使用できます。
KE法	固有ベクトル出力はサポートされません。 .h3d出力ファイルは使用できません。
PK法 / PKNL法	固有ベクトル出力をポイントの形式で使用できます。各ポイントが1つのモード番号に対応します。速度のFLFACTエントリで速度が正の場合は、出力を使用できません。特定の速度の出力が必要な場合は、FLFACTバルクデータエントリでそれらの速度の符号を負にできます。このページの問題の設定に例を示します。これは表記上の規則に過ぎず、負の速度は物理的意味を持ちません。FLFACTエントリ内の速度の符号に関係なく、計算では必ず正の値が使用されます。