周波数応答解析
定常状態の振動荷重に対する構造の応答を計算します。
一般的な用途は、自動車、回転機械、トランスミッション、パワートレインシステムから発生する騒音振動(NVH)の解析です。
周波数応答解析を使用して、静的な周波数領域における構造の応答を計算します(実際には時刻歴解析となります)。荷重は正弦波荷重を使用します。単純なケースは、特定の周波数における所定の振幅の荷重の場合です。応答は同一の周波数で計算されますが、減衰によって位相のずれが発生します(図 1)。
荷重としては、力、変位、速度および加速度が使用できます。これらは、加振周波数 に依存します。
周波数応答解析の結果は変位、速度、加速度、力、応力、およびひずみが計算されます。通常、応答は複素数であり、絶対値および位相角として、または実数および虚数として算出されます。
図 1. 周波数応答解析の荷重と応答
直接法による周波数応答解析
直接法による周波数応答解析を使用すると、一連の複素行列の方程式を解くことで、離散的な励振周波数 に対する構造の応答を直接計算できます。
- は質量マトリックス
- は減衰マトリックス
- は剛性マトリックス
- は変位ベクトル
- は荷重ベクトル
- 一定の構造減衰係数を使用する方法。
- 材料特性で定義した減衰係数、ブッシュおよびスプリング要素の特性定義のを使用し、構造要素減衰が定義されます。これらはマトリックスを形成します。
- ダンパー要素によって粘性減衰が生成されます。これらはマトリックスを形成します。
運動方程式は、複素代数を使用して直接解析されます。
直接法による周波数応答解析の実行
荷重および境界条件は、入力デックのバルクデータエントリセクションで定義します。これらをサブケース情報セクション内で参照する必要があり、これにはSUBCASEでSPCステートメントおよびDLOADステートメントを使用します。
OptiStructは、直接法による周波数応答解析の慣性リリーフをサポートしていません。これを試みると、エラーが出されます。
周波数セットは、FREQUENCYステートメントを使用して参照する必要があります。
さまざまな減衰要素と材料減衰のほか、一定の構造減衰を適用することもできます。これにはPARAM,Gを使用します。
モーダル法による周波数応答解析
モーダル法では、まずノーマルモード解析を実行し、システムの固有値 およびこれに対応するシステムの固有ベクトルを取得します。
モード質量マトリックスとモード剛性マトリックスは、対角マトリックスです。固有ベクトルが質量マトリックスを基準に正規化されている場合、モーダル質量マトリックスは単位マトリックス、モーダル剛性マトリックスはシステムの固有値を保持する対角マトリックスです。これによってシステム方程式は、の各コンポーネントについての一連の非連成の方程式に単純化されるので、容易に解析できます。
したがって、マトリックス と は、一般には対角マトリックスとはなりません。この結果、連成問題は、直接法によって解析されるシステムと似たものになりますが、自由度ははるかに小さくなります。これは、直接法を使用して解析できます。
運動方程式の評価は、方程式の非連成が維持されている場合は、より高速になります。非連成は、各モードに個別に減衰を適用する場合に維持できます。これには、減衰値 と固有振動数の関係を示す減衰テーブルTABDMP1を使用します。このアプローチが使用される場合、構造要素と粘性減衰は定義できません。
- モード減衰比
- モード固有値
- 構造減衰
- 臨界減衰
- 品質係数
モーダル法では、ノーマルモード解析を制御するためにMETHODステートメントが必要です。そのMETHODステートメントはEIGRLまたはEIGRAバルクデータエントリのどちらかを参照します。
剰余ベクトル生成(精度の向上)
モーダル法の精度は、固有ベクトルのマトリックスへの動的荷重に基づいて静解析の変位ベクトルを追加することによって著しく向上します。これらのベクトルは剰余ベクトルと呼ばれ、この方法はモーダル加速度と呼ばれます。
- unit load methodは、動的荷重の自由度における単位ベクトルである静的荷重に基づき、剰余ベクトルを生成します。すなわち、剰余ベクトル生成のための静的荷重は、動的荷重が適用される自由度における単位ベクトルです。剰余ベクトルの数は、適用された自由度の数と等しくなります。一般的に、この方法がより正確であることから、これがデフォルトの方法となっています。
- applied load methodは、周波数ゼロにおける動的荷重を使用して、最大2つの剰余ベクトルを生成します。動的荷重の実数部分と虚数部分が同じである場合、またはそれらのいずれかがゼロである場合、どちらか一方のみが使用されます。
図 2. モーダル法による周波数応答解析(FRA)の結果精度に剰余ベクトルの使用が及ぼす影響の. 正確な直接法の結果との比較
強制運動を使用したモーダル法による周波数応答解析
- 相対的な方法:
この場合、解析が2段階で進められます。まず、強制運動を使用した静解析が行われ、静的変位が求められます。次に、先に計算された静的変位と固有ベクトルを使用して動的解析が行われます。この方法は、あまり効率的ではありませんが、より正確な解が得られるデフォルトの方法です。
- 全体的 / 絶対的な方法:
この場合、解析が1段階で進められ、静的変位の計算が不要です。モード動的荷重の寄与が、SPCDの自由度で適用された変位 / 速度 / 加速度から直接もたらされます。この方法は、静的変位ベクトルの計算が必要ないため、計算効率が高くなります。
これらの方法を使用して計算を制御するには、PARAM, ENFMETHをご参照ください。
モーダル法による周波数応答解析の実行
入力
荷重および境界条件は、入力デックのバルクデータエントリセクションで定義します。これらをサブケース情報セクション内で参照する必要があり、これにはSUBCASEでSPCステートメントおよびDLOADステートメントを使用します。
周波数セットは、FREQUENCYステートメントを使用して参照する必要があります。モーダル法では、ノーマルモード解析を制御するためにMETHODステートメントが必要です。計算結果を保存するために、それより以前に保存された固有ベクトルを、EIGVRETRIEVEサブケースステートメントを使って読み込むことができます。
さまざまな減衰要素と材料減衰のほか、一定の構造減衰を適用することもできます。これにはPARAM, Gを使用します。
モード減衰を適用するには、減衰テーブルTABDMP1のSDAMPING参照を使用します。減衰テーブルを適用する方法を定義するには、PARAM,KDAMPを使用します。
直接法およびモーダル法周波数応答解析では、MATiエントリ上の対応するフィールドのTABLEDiエントリを通して、周波数依存材料(MATFiバルクデータエントリ)を使用できます。MATF1、MATF2、MATF3、MATF8、MATF9、MATF10バルクデータエントリを使用して、現在使用可能な周波数依存材料を定義できます。
PBUSHTエントリ上の対応するフィールドのTABLEDiエントリを通して、周波数応答解析で周波数依存特性(PBUSHTバルクデータエントリ)を使用することもできます。
モーダルFRF解析用のデフォルトソルバー
デフォルトでは、モーダル法による周波数応答解析には、標準の内部ソルバーが使用されます。ただし、一部のモデルのクラスについては、高速なモーダル解法(FASTFR)または高速周波数応答ソルバー(FastFRS)を使用することにより、パフォーマンスの向上が得られる場合があります。
剰余ベクトル
剰余ベクトルはモーダルFRF / 音場 / 過渡解析に関係します。これらは解析の精度を上げるため、デフォルトで計算されます。RESVEC計算はケースコントロール文で制御できます。
RESVEC(APPLOD/UNITLOD,DAMPLOD/NODAMP)=Value
ここで、ValueはYesまたはNoのいずれかです。ValueがNoの場合、カッコ内のキーワードは無視され、全てのRESVEC計算がオフになります。キーワードAPPLODではモーダルFRF / 音場 / 過渡解析の動的荷重を元にRESVECが生成されます。キーワードUNITLOD動的荷重の自由度への単位荷重を元にRESVECが生成されます。キーワードDAMPLODでは、粘性減衰自由度への単位荷重を元に粘性減衰のRESVECが生成されます。キーワードNODAMPでは粘性減衰のRESVEC生成がOFFになり、ない場合はデフォルトで生成されます。DAMPLODとNODAMPがケースコントロール内でオプションであったとしても、それらは、モデル内のすべてのモーダルFRF / 音場 / 過渡サブケースに適用されるグローバルスイッチです。
基本的な固有値がLanczos法で計算されると 、デフォルトのRESVECが作用荷重と粘性減衰の自由度に基づいて生成されます。基本的な固有値解析がAMSESまたはAMLSを用いてなされた場合、デフォルトの RESVECは荷重自由度への単位荷重と粘性減衰自由度に基づいて生成されます。強制変位、速度または加速度が定義されている場合、剰余ベクトルは常に生成されます。加えて、USET U6データがある場合、AMSESまたはAMLS固有値ソルバー使用時には剰余ベクトルが計算されます。USET U6 剰余ベクトルは、Lanczos 固有値ソルバーが用いられている場合には計算されません。
剰余ベクトルが含まれる場合、拘束されていないモデルにデフォルトで慣性リリーフが適用されます。慣性リリーフがRESVECにとって好ましくない場合、PARAM, INREL, 0を用いてOFFにできます。
剰余ベクトルが含まれる場合、FRF / 過渡サブケースの基本的な固有値解析による固有モードが慣性リリーフで使用されます。限界値(FZERO)より下の全ての固有値のモードは慣性リリーフ解析での剛体モードとして用いられます。FZEROより下に固有モードがない場合、 モデル形状に基づいて6までの全体の剛体モードが内部的に生成され慣性リリーフに用いられます。PARAM,FZERO, Valueを用いてFZEROを設定することが可能です。FZEROのデフォルトの値は0.1です。
出力
周波数応答解析の結果は、変位、速度、加速度、力、応力、およびひずみが計算されます。STRESS、STRAIN、DISPLACEMENTなどの通常の出力エントリを使用して、対応する出力値を要求できます。
DISP(MODAL)を使用すると、モーダルFRF解析 / 過渡応答で固有ベクトルのみを出力することができます。DISP(MODAL,NODAL)を使用すると、固有ベクトルとそれに対応するFRF解析 / 過渡応答の結果を両方出力することができます。
PARAM, ENFMOTN, RELを使用して、指定された強制運動に対する変位、速度、加速度の出力を生成できます。このような場合、応力や力などの後続で計算される出力も、指定された強制運動に対して生成されます。PARAM, ENFMOTN, TOTAL/ABSを使用して、指定された強制モーションを含む合計出力値を生成できます(デフォルトはTOTAL/ABS)。