(57)【要約】【課題】複雑な構造物でも、正確かつ短時間で振動モードを判定することができる解析モデルの振動モード判定装置及び方法を提供する。【解決手段】構造物の振動試験結果及び解析結果を求め（Ｓ１〜Ｓ２）、振動試験結果内の１つの振動モードの固有振動数の±１０％内にある固有振動数を有する振動モードを解析結果から第１の振動モードとして抽出し（Ｓ３）、振動試験結果内の１つの振動モードの変位が大きい部位と、第１の振動モードの変位が大きい部位とを比較し、同じ部位の変位が大きい振動モードを第１の振動モードから第２の振動モードとして選択し（Ｓ４）、振動試験結果内の１つの振動モードのモードベクトルと第２の振動モードのモードベクトルとからＭＡＣを求め、当該ＭＡＣが最大となる第２の振動モードを、振動試験結果内の１つの振動モードに対応する振動モードと判定する（Ｓ５〜Ｓ６）。【選択図】図２
2016177530.tif 000002
【特許請求の範囲】
【請求項１】
解析対象の構造物の解析モデルを数値解析して、少なくとも振動モード、モードベクトル及び固有振動数を含む解析結果を求める解析部と、
前記構造物に実際の振動試験を行って求められた少なくとも振動モード、モードベクトル及び固有振動数を含む振動試験結果と前記解析結果とを用いて、前記振動試験結果内の振動モードと前記解析結果内の振動モードとの対応関係を判定するモード判定部とを有し、
前記モード判定部は、
前記振動試験結果内の任意の１つの振動モードについて、当該振動モードの固有振動数の所定範囲内にある固有振動数を有する振動モードを前記解析結果から第１の振動モードとして抽出し、
前記任意の１つの振動モードの変位が大きい部位と、前記第１の振動モードの変位が大きい部位とを比較して、同じ部位の変位が大きい振動モードを前記第１の振動モードから第２の振動モードとして選択し、
前記任意の１つの振動モードのモードベクトルと、前記第２の振動モードのモードベクトルとから、ＭＡＣ（Modal Assurance Criterion）を求め、当該ＭＡＣが最大となる前記第２の振動モードを、前記任意の１つの振動モードに対応する振動モードと判定する
ことを特徴とする解析モデルの振動モード判定装置。
【請求項２】
請求項１に記載の解析モデルの振動モード判定装置において、
更に、前記解析モデルの精度を評価する評価部を有し、
前記評価部は、
対応すると判定された振動モード同士の全てについて、当該振動モード同士のモードベクトルから前記ＭＡＣを求め、求めた全ての前記ＭＡＣが所定値以上である場合、前記解析モデルの精度が高いと評価する
ことを特徴とする解析モデルの振動モード判定装置。
【請求項３】
請求項２に記載の解析モデルの振動モード判定装置において、
更に、前記解析モデルをチューニングするチューニング部を有し、
前記チューニング部は、
前記評価部で評価した前記ＭＡＣの少なくとも１つが所定値未満である場合、対応すると判定された振動モード同士の対応関係に基づいて、前記解析モデルの剛性値を最適化する
ことを特徴とする解析モデルの振動モード判定装置。
【請求項４】
解析対象の構造物に実際の振動試験を行って、少なくとも振動モード、モードベクトル及び固有振動数を含む振動試験結果を求め、
前記構造物の解析モデルを数値解析して、少なくとも振動モード、モードベクトル及び固有振動数を含む解析結果を求め、
前記振動試験結果内の任意の１つの振動モードについて、当該振動モードの固有振動数の所定範囲内にある固有振動数を有する振動モードを前記解析結果から第１の振動モードとして抽出し、
前記任意の１つの振動モードの変位が大きい部位と、前記第１の振動モードの変位が大きい部位とを比較して、同じ部位の変位が大きい振動モードを前記第１の振動モードから第２の振動モードとして選択し、
前記任意の１つの振動モードのモードベクトルと、前記第２の振動モードのモードベクトルとから、ＭＡＣ（Modal Assurance Criterion）を求め、当該ＭＡＣが最大となる前記第２の振動モードを、前記任意の１つの振動モードに対応する振動モードと判定する
ことを特徴とする解析モデルの振動モード判定方法。
【請求項５】
請求項４に記載の解析モデルの振動モード判定方法において、
対応すると判定された振動モード同士の全てについて、当該振動モード同士のモードベクトルから前記ＭＡＣを求め、求めた全ての前記ＭＡＣが所定値以上である場合、前記解析モデルの精度が高いと評価する
ことを特徴とする解析モデルの振動モード判定方法。
【請求項６】
請求項５に記載の解析モデルの振動モード判定方法において、
評価した前記ＭＡＣの少なくとも１つが所定値未満である場合、対応すると判定された振動モード同士の対応関係に基づいて、前記解析モデルの剛性値を最適化する
ことを特徴とする解析モデルの振動モード判定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、解析モデルのチューニングの際の振動モード判定装置及び方法に関する。
【背景技術】
【０００２】
航空機においては、例えば、そのフラッタ速度を解析するために、ＦＥＭ（Finite Element Method；有限要素法）の解析モデルの振動特性と実機の振動特性を合わせる必要がある。そのために解析モデルをチューニングする作業を実施している。
【０００３】
航空機の実機の振動特性はＧＶＴ（Ground Vibration Test；全機地上振動試験）にて計測される。このＧＶＴとは、航空機の振動特性（振動数、振動モード形、構造減衰率等）を計測する試験である。解析モデルチューニングでは、解析モデルの固有振動数と振動モード形がＧＶＴ結果に合うように、解析モデルの剛性値を変更する。この際、ＧＶＴ結果と解析結果で各振動モードの固有振動数、振動モード形が全体的に合うようにするため、最適化ツールを使用している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第４６２６３５１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
解析モデルチューニングにおいて、固有振動数は最適化ツールの評価関数として使用するための表現が容易であるが、振動モード形は定量化が困難である。振動モード形を定量化する指標として、以下に示すＭＡＣ（Modal Assurance Criterion）があり、この値を評価関数として使用しているが、近傍の振動数帯に類似の振動モードがあると、振動モード判定がうまくいかないという問題がある。
【０００６】
ここで、ＭＡＣとは、振動モード形（モードベクトル）について、解析値と試験値の相関を判断する指標の１つであり、次式で定義される。
【０００７】
【数１】
2016177530.tif 000003
【０００８】
上述した問題を、図３を参照して説明する。なお、図３は、試験及び解析における航空機の振動モードを説明する図であるが、ここでは、図を簡単にするため、航空機の胴体と片方の主翼のみ図示している。
【０００９】
例えば、航空機の試験結果のモード１が主翼の動きが顕著な主翼系の振動モードであり、解析結果のモード１Ａが胴体の動きが顕著な胴体系の振動モードであり、解析結果のモード１Ｂが主翼の動きが顕著な主翼系の振動モードであるとする。このような場合において、モード１に対するモード１Ａ及びモード１ＢのＭＡＣを各々計算したとき、両者とも同等の値となることがある。しかしながら、モード１Ａは胴体の動きが顕著な振動モードであり、試験結果のモード１とは異なる振動モードである。このように、ＭＡＣだけで機械的にモード判定すると、異なる振動モードを同じ振動モードと判定してしまうおそれがある。
【００１０】
上記特許文献１では、振動モード判定の自動化にＭＡＣを用いているが、この特許文献１は、単純な軸対称形状が対象であり、航空機のような複雑形状の構造体を対象としておらず、そのままでの適用は不可であり、もし適用したとしても、上述したように、異なる振動モードを同じ振動モードと判定してしまうおそれがある。
【００１１】
本発明は上記課題に鑑みなされたもので、複雑な構造物でも、正確かつ短時間で振動モードを判定することができる解析モデルの振動モード判定装置及び方法を提供することを目的とする。
【課題を解決するための手段】
【００１２】
上記課題を解決する第１の発明に係る解析モデルの振動モード判定装置は、
解析対象の構造物の解析モデルを数値解析して、少なくとも振動モード、モードベクトル及び固有振動数を含む解析結果を求める解析部と、
前記構造物に実際の振動試験を行って求められた少なくとも振動モード、モードベクトル及び固有振動数を含む振動試験結果と前記解析結果とを用いて、前記振動試験結果内の振動モードと前記解析結果内の振動モードとの対応関係を判定するモード判定部とを有し、
前記モード判定部は、
前記振動試験結果内の任意の１つの振動モードについて、当該振動モードの固有振動数の所定範囲内にある固有振動数を有する振動モードを前記解析結果から第１の振動モードとして抽出し、
前記任意の１つの振動モードの変位が大きい部位と、前記第１の振動モードの変位が大きい部位とを比較して、同じ部位の変位が大きい振動モードを前記第１の振動モードから第２の振動モードとして選択し、
前記任意の１つの振動モードのモードベクトルと、前記第２の振動モードのモードベクトルとから、ＭＡＣ（Modal Assurance Criterion）を求め、当該ＭＡＣが最大となる前記第２の振動モードを、前記任意の１つの振動モードに対応する振動モードと判定する
ことを特徴とする。
【００１３】
上記課題を解決する第２の発明に係る解析モデルの振動モード判定装置は、
上記第１の発明に記載の解析モデルの振動モード判定装置において、
更に、前記解析モデルの精度を評価する評価部を有し、
前記評価部は、
対応すると判定された振動モード同士の全てについて、当該振動モード同士のモードベクトルから前記ＭＡＣを求め、求めた全ての前記ＭＡＣが所定値以上である場合、前記解析モデルの精度が高いと評価する
ことを特徴とする。
【００１４】
上記課題を解決する第３の発明に係る解析モデルの振動モード判定装置は、
上記第２の発明に記載の解析モデルの振動モード判定装置において、
更に、前記解析モデルをチューニングするチューニング部を有し、
前記チューニング部は、
前記評価部で評価した前記ＭＡＣの少なくとも１つが所定値未満である場合、対応すると判定された振動モード同士の対応関係に基づいて、前記解析モデルの剛性値を最適化する
ことを特徴とする。
【００１５】
上記課題を解決する第４の発明に係る解析モデルの振動モード判定方法は、
解析対象の構造物に実際の振動試験を行って、少なくとも振動モード、モードベクトル及び固有振動数を含む振動試験結果を求め、
前記構造物の解析モデルを数値解析して、少なくとも振動モード、モードベクトル及び固有振動数を含む解析結果を求め、
前記振動試験結果内の任意の１つの振動モードについて、当該振動モードの固有振動数の所定範囲内にある固有振動数を有する振動モードを前記解析結果から第１の振動モードとして抽出し、
前記任意の１つの振動モードの変位が大きい部位と、前記第１の振動モードの変位が大きい部位とを比較して、同じ部位の変位が大きい振動モードを前記第１の振動モードから第２の振動モードとして選択し、
前記任意の１つの振動モードのモードベクトルと、前記第２の振動モードのモードベクトルとから、ＭＡＣ（Modal Assurance Criterion）を求め、当該ＭＡＣが最大となる前記第２の振動モードを、前記任意の１つの振動モードに対応する振動モードと判定する
ことを特徴とする。
【００１６】
上記課題を解決する第５の発明に係る解析モデルの振動モード判定方法は、
上記第４の発明に記載の解析モデルの振動モード判定方法において、
対応すると判定された振動モード同士の全てについて、当該振動モード同士のモードベクトルから前記ＭＡＣを求め、求めた全ての前記ＭＡＣが所定値以上である場合、前記解析モデルの精度が高いと評価する
ことを特徴とする。
【００１７】
上記課題を解決する第６の発明に係る解析モデルの振動モード判定方法は、
上記第５の発明に記載の解析モデルの振動モード判定方法において、
評価した前記ＭＡＣの少なくとも１つが所定値未満である場合、対応すると判定された振動モード同士の対応関係に基づいて、前記解析モデルの剛性値を最適化する
ことを特徴とする。
【発明の効果】
【００１８】
本発明によれば、複雑な構造物でも、正確かつ短時間で振動モードを判定することができ、その結果、振動解析のモデルチューニング作業が効率的になり、その作業時間を短縮することができる。
【図面の簡単な説明】
【００１９】
【図１】（ａ）は、本発明に係る解析モデルの振動モード判定方法を実施する処理装置を説明する構成図であり、（ｂ）は、当該処理装置における機能を説明するブロック図である。
【図２】本発明に係る解析モデルの振動モード判定方法を説明するフローチャートである。
【図３】ＭＡＣを用いた振動モード判定の問題点を説明する図である。
【発明を実施するための形態】
【００２０】
以下、図１及び図２を参照して、本発明に係る解析モデルの振動モード判定装置及び方法の実施形態を説明する。
【００２１】
（実施例１）
図１（ａ）は、本実施例の解析モデルの振動モード判定方法を実施する処理装置を説明する構成図であり、図１（ｂ）は、当該処理装置における機能を説明するブロック図である。
【００２２】
本実施例の処理装置１０は、バス又はネットワークを介して相互に接続されるＣＰＵ（Central Processing Unit）１１と、ＲＡＭ（Random Access Memory）１２と、入力装置１３と、出力装置１４と、記憶装置１５等を有する。つまり、コンピュータの構成を有している。
【００２３】
入力装置１３は、例えば、キーボードやマウス等の入力用ユーザインタフェースであり、ユーザに操作されることにより、各種データをＣＰＵ１１や記憶装置１５に入力するものである。この入力装置１３は、後述する振動試験結果２１や解析モデルの解析結果２４が蓄積された他の記憶装置から当該結果を授受するインタフェースを含んでも良い。出力装置１４は、例えば、モニタやプリンタなどであり、ＣＰＵ１１から出力される情報をユーザが視認可能に出力するものである。記憶装置１５は、例えば、ハードディスク、メモリ、ディスク等の記憶装置である。この記憶装置１５は、ネットワークを介して外部に接続されたストレージシステムでも良い。
【００２４】
そして、図１（ｂ）に示すように、処理装置１０は、その機能として、解析対象の構造物の解析モデル２２を用いて解析を行う解析部２３と、振動試験結果２１と解析部２３の解析結果２４とを用いて振動モードの判定を行うモード判定部２５と、モード判定部２５で判定した振動モードを用いて解析モデル２２の評価を行う評価部２６と、評価部２６での評価に基づいて、解析モデル２２のチューニングを行うチューニング部２７とを有している。
【００２５】
解析モデル２２、振動試験結果２１及び解析結果２４はデータであり、ＲＡＭ１２や記憶装置１５に記憶されており、必要に応じて、ＣＰＵ１１がデータを読み込んで使用している。解析モデル２２には、解析対象の構造物の形状、材料、剛性値などのデータが含まれており、振動試験結果２１及び解析結果２４には、振動モード、振動モード形、固有振動数などのデータが含まれている。振動試験結果２１は、実際の解析対象の構造物に対する加振などによる振動試験の結果であり、解析結果２４は、解析モデル２２に対する数値解析による結果である。
【００２６】
又、解析部２３、モード判定部２５、評価部２６及びチューニング部２７はプログラムであり、ＲＡＭ１２や記憶装置１５に記憶されており、必要に応じて、ＣＰＵ１１がプログラムを読み込んで実行している。解析部２３では、解析モデル２２の解析を行うため、例えば、ＦＥＭなどの数値解析法のプログラムが用いられている。チューニング部２７では、解析モデル２２のチューニングを行うため、例えば、遺伝的アルゴリズムなどの最適化ツールのプログラムが用いられている。なお、モード判定部２５及び評価部２６は、後述する動作を行うための特有のプログラムとなっている。
【００２７】
次に、図２に示すフローチャートも参照して、解析部２３、モード判定部２５、評価部２６及びチューニング部２７の動作を詳細に説明する。なお、ここでは、一例として、解析対象の構造物を航空機とし、チューニング対象とする振動モードを主翼曲げ１次モードとして以降の説明を行う。又、航空機は、主な構成として、主翼、水平尾翼、垂直尾翼及び胴体などの各部位から構成されるものとする。
【００２８】
（ステップＳ１）
実際の航空機に対する振動試験を行い、振動試験結果２１を計算する。具体的には、実際の航空機に対して、上述したＧＶＴ（全機地上振動試験）を行い、実際の航空機の各部の加速度などのデータを計測し、計測したデータを用いて、振動モードの固有振動数、振動モード形（モードベクトル）を計算する。振動試験結果２１としては、例えば、下記表１に示すようなデータが計算される。
【００２９】
【表１】
2016177530.tif 000004
【００３０】
（ステップＳ２）
航空機の解析モデル２２を数値解析して、解析結果２４を計算する（解析部２３）。具体的には、ＦＥＭなどの数値解析法を用いて、振動モードの固有振動数、振動モード形（モードベクトル）を計算する。最初は、解析モデル２２の剛性値としては設計値（初期値）を用いる。解析結果２４としては、例えば、下記表２に示すようなデータが計算される。
【００３１】
【表２】
2016177530.tif 000005
【００３２】
振動試験結果２１と解析結果２４が計算された時点では、振動試験結果２１と解析結果２４の振動モードについて、どれとどれが対応しているかは分からない。例えば、振動モードａが振動モードＡ（主翼曲げ１次モード）に対応しているとは限らない。そこで、本実施例では、モード判定部２５が、振動試験結果２１と解析結果２４に基づき、以下のステップＳ３〜Ｓ７の手順を用いて、振動モード判定を行っている。
【００３３】
（ステップＳ３）
ここでは、上述したように、主翼曲げ１次モードをチューニング対象にしているので、振動試験結果２１の主翼曲げ１次モードである振動モードＡの固有振動数の±１０％以内の固有振動数を有する振動モードを解析結果２４から第１の振動モードとして抽出する。解析結果２４から抽出されなかった振動モードは、主翼曲げ１次モードの判定対象から除外し、後述するステップＳ７で処理を行う。なお、ここでは、一例として、抽出範囲を±１０％以内としているが、これは標準的な閾値であり、状況に応じて、適宜に変更可能である。
【００３４】
（ステップＳ４）
振動試験結果２１の振動モードＡは、主翼の構造グリッドの変位が最大であり、主翼の動きが支配的な振動モードである。そして、解析結果２４から抽出した第１の振動モードについても、主翼、水平尾翼、垂直尾翼及び胴体を構成する構造グリッドの最大変位を評価し、主翼の構造グリッドが最大変位の７０％以上であれば、主翼の動きが支配的な振動モードであると判定し、第１の振動モードから第２の振動モードとして選択する。つまり、振動試験結果２１の振動モードＡの支配的振動部位の変位が、解析結果２４から抽出した第１の振動モードでも大きいか確認している。主翼の動きが支配的でないと判定した第１の振動モードについては、主翼曲げ１次モードの判定対象から除外し、後述するステップＳ７で処理を行う。なお、ここでは、一例として、支配的振動部位の変位を最大変位の７０％以上としているが、これは標準的な閾値であり、状況に応じて、適宜に変更可能である。
【００３５】
（ステップＳ５）
ステップＳ４で主翼の動きが支配的と判定した解析結果２４の第２の振動モードのモードベクトルを用いて、上述したＭＡＣを計算し、計算したＭＡＣが最大値のものを選択する。例えば、振動試験結果２１の振動モードＡの固有振動数の±１０％以内の固有振動数であり、かつ、主翼の動きが支配的と判定した解析結果２４の第２の振動モードを、振動モードａと振動モードｂとすると、振動モードａのモードベクトルａと振動モードｂのモードベクトルｂの各々に対して、振動試験結果２１の振動モードＡのモードベクトルＡとのＭＡＣを計算し、ＭＡＣが大きい第２の振動モードを選択することになる。選択されなかった第２の振動モードは、主翼曲げ１次モードの判定対象から除外し、後述するステップＳ７で処理を行う。
【００３６】
（ステップＳ６）
ステップＳ５で計算したＭＡＣが最大となる解析結果２４の第２の振動モードを、主翼曲げ１次モードと判定し、振動試験結果２１の振動モードＡに対応すると判定する。このとき、ＭＡＣの値により、対応すると判定された解析結果２４の第２の振動モードが振動試験結果２１の振動モードＡをどの程度正しく表現しているかを評価している。ＭＡＣは０〜１の値を持ち、ＭＡＣが１に近いほど、正しく表現していることになる。例えば、計算したＭＡＣが０．９以上であれば、正しく表現していると考えて良い。なお、ここでは、一例として、ＭＡＣを０．９以上としているが、これは標準的な閾値であり、状況に応じて、適宜に変更可能である。
【００３７】
このように、上述したステップＳ３〜Ｓ６では、振動試験結果２１の任意の１つの振動モードと解析結果２４の振動モードとのＭＡＣを用いて、振動試験結果２１の任意の１つの振動モードに対応する解析結果２４の振動モードを正確に選択すると共に、選択した解析結果２４の振動モードが振動試験結果２１の任意の１つの振動モードをどの程度正しく表現されているか評価している。これにより、振動モード判定の自動化を行うと共に、振動モード形（モードベクトル）の評価を実施している。
【００３８】
（ステップＳ７）
ステップＳ３で解析結果２４から抽出されなかった振動モード、ステップＳ４で主翼の動きが支配的でないと判定された振動モード（第１の振動モード）、ステップＳ５で選択されなかった振動モード（第２の振動モード）については、振動試験結果２１のチューニング対象の振動モードを主翼曲げ１次モードとは異なる振動モードに変更して、上述したステップＳ３〜Ｓ７を実施すれば、振動試験結果２１の異なる振動モードに対応する解析結果２４の振動モードを判定することができる。例えば、振動試験結果２１のチューニング対象の振動モードを振動モードＢ（主翼曲げ２次モード）に変更すれば、振動モードＢ（主翼曲げ２次モード）に対応する振動モードを解析結果２４の中から正確に選択すると共に評価することができる。
【００３９】
同様に、振動試験結果２１のチューニング対象の振動モードを振動モードＣ、Ｄ、Ｅ・・・に順次変更して、上述したステップＳ３〜Ｓ７を実施すれば、振動試験結果２１のチューニング対象の振動モードに対応する解析結果２４の振動モードを正確に選択すると共に評価することができる。
【００４０】
（ステップＳ８）
振動試験結果２１の全ての振動モードＡ、Ｂ、Ｃ、Ｄ、Ｅ・・・について、対応する解析結果２４の振動モードの選択及び評価を行うと、この際に計算した複数のＭＡＣに基づいて、解析モデル２２の精度を評価することができる（評価部２６）。例えば、計算した複数のＭＡＣが全て０．９以上であれば、解析モデル２２の精度が高いと評価して良い。なお、ここでも、一例として、ＭＡＣを０．９以上としているが、これは標準的な閾値であり、状況に応じて、適宜に変更可能である。
【００４１】
（ステップＳ９〜Ｓ１０）
計算した複数のＭＡＣが所定値（例えば、０．９）以上であれば、所望の精度が出ていると判断し、一連の処理手順を終了するが、所望の精度が出ていないときは、ステップＳ１０へ進む。ステップＳ１０では、遺伝的アルゴリズム等を利用した最適化ツールにより、解析モデル２２の剛性値を変更して（チューニング部２７）、ステップＳ２へ戻る。このとき、最適化ツールには、上述した手順により得られた振動モード判定の結果（振動試験結果２１と解析結果２４の振動モード同士の対応関係）が初期値として与えられる。
【００４２】
そして、剛性値が変更された新たな解析モデル２２を用いて、解析結果２４を計算し直して、上述したステップＳ３〜Ｓ８を再び実施することになる。これにより、解析モデル２２をより精度が高い解析モデルへチューニング（合わせ込み）をすることが可能となる。
【００４３】
従来はＭＡＣだけで機械的にモード判定していたので、異なる振動モードを同じ振動モードと判定してしまうおそれがあったが、本実施例では、以上説明した手順により、振動試験結果２１と解析結果２４の振動モード同士の対応関係を判定しているので、振動モードの対応関係を正確かつ短時間で判定することができる。これにより、最適化ツールのルーチン内の振動モード判定の確からしさが向上し、解析モデルのチューニング作業を自動化できるため、作業時間短縮に寄与することになる。
【産業上の利用可能性】
【００４４】
本発明は、複雑な構造物、例えば、航空機（戦闘機、民間機）などの解析モデルのチューニングの際に、振動モードを判定するときに好適なものである。
【符号の説明】
【００４５】
１０  処理装置
１１  ＣＰＵ
１２  ＲＡＭ
１３  入力装置
１４  出力装置
１５  記憶装置
２１  振動試験結果
２２  解析モデル
２３  解析部
２４  解析結果
２５  モード判定部
２６  評価部
２７  チューニング部
【図１】
2016177530.tif 000006
【図２】
2016177530.tif 000007
【図３】
2016177530.tif 000008