(57)【特許請求の範囲】
【請求項１】
複数の振動要素を有する振動系の振動解析装置であって、
各振動要素の加速度、振動系の加振力、及び各振動要素を伝わる伝達力の変数のうち、いずれか複数を観測変数として取得する観測変数取得部と、
複数の振動要素間の接続関係を表して振動系をモデル化した場合における各振動要素の運動方程式を用いて、前記変数のうち観測変数以外の未知変数を、観測変数取得部で取得された観測変数に基づいて推定する未知変数推定部と、
を備え、
未知変数に関わる列ベクトルをｘunknown、観測変数に関わる列ベクトルをｘmeasured
とし、前記運動方程式を行列Ｈ1，Ｈ2を用いて、
【数１】
0006260477.tif 000020
で表した場合に、観測変数取得部は、行列Ｈ1が列フルランクになる条件を満たすように観測変数を取得する、振動解析装置。
【請求項２】
請求項１に記載の振動解析装置であって、
振動要素を伝わる伝達力に基づいて振動要素間の弾性係数及び減衰係数を推定するパラメータ推定部を備える、振動解析装置。
【請求項３】
請求項１又は２に記載の振動解析装置であって、
観測変数取得部は、各振動要素の加速度を観測変数として取得し、
未知変数推定部は、観測変数取得部で取得された各振動要素の加速度に基づいて、振動系の加振力及び各振動要素を伝わる伝達力を未知変数として推定する、振動解析装置。
【請求項４】
複数の振動要素を有する振動系の振動解析方法であって、
各振動要素の加速度、振動系の加振力、及び各振動要素を伝わる伝達力の変数のうち、いずれか複数を観測変数として取得する観測変数取得処理と、
複数の振動要素間の接続関係を表して振動系をモデル化した場合における各振動要素の運動方程式を用いて、前記変数のうち観測変数以外の未知変数を、観測変数取得処理で取得された観測変数に基づいて推定する未知変数推定処理と、
を含み、
未知変数に関わる列ベクトルをｘunknown、観測変数に関わる列ベクトルをｘmeasured
とし、前記運動方程式を行列Ｈ1，Ｈ2を用いて、
【数２】
0006260477.tif 000021
で表した場合に、観測変数取得部は、行列Ｈ1が列フルランクになる条件を満たすように観測変数を取得する、振動解析方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数の振動要素を有する振動系の振動解析装置及び振動解析方法に関する。
【背景技術】
【０００２】
下記特許文献１では、入出力間の振動伝達特性を求めるために、実稼働状態にて入出力ともに加速度を計測し、入出力間の伝達率（出力加速度／入力加速度）を計算している。下記特許文献２では、入出力間の振動伝達特性を求めるために、入力として加速度と力の両方を計測するとともに出力として加速度を計測し、加速度入力に対しては伝達率を計算し、力入力に対しては伝達関数（加速度／力）を計算している。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００７−５７４６０号公報
【特許文献２】特開２０１１−８１５４３号公報
【特許文献３】特開２００９−９７９７３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
特許文献１，２において計算される伝達率や伝達関数は、振動系（システム）の入出力間の周波数毎のノンパラメトリックモデルであるため、伝達率や伝達関数を計算するだけでは、振動系（システム）内部の個々の振動要素間の弾性係数や減衰係数のパラメータを求めることは困難であり、パラメトリックモデルを構築することは困難である。振動系の弾性係数や減衰係数のパラメータを求めるためには、振動系の振動要素を伝わる伝達力を計測する必要があるが、力は加速度に比べて直接的な計測が困難である。
【０００５】
本発明は、振動系の振動要素を伝わる伝達力の推定を可能にすることを目的とする。
【課題を解決するための手段】
【０００６】
本発明に係る振動解析装置及び振動解析方法は、上述した目的を達成するために以下の手段を採った。
【０００７】
本発明に係る振動解析装置は、複数の振動要素を有する振動系の振動解析装置であって、各振動要素の加速度、振動系の加振力、及び各振動要素を伝わる伝達力の変数のうち、いずれか複数を観測変数として取得する観測変数取得部と、複数の振動要素間の接続関係を表して振動系をモデル化した場合における各振動要素の運動方程式を用いて、前記変数のうち観測変数以外の未知変数を、観測変数取得部で取得された観測変数に基づいて推定する未知変数推定部と、を備えることを要旨とする。上記の加速度については、並進加速度であってもよいし、角加速度であってもよい。そして、上記の加振力及び伝達力についても、並進方向の力であってもよいし、回転方向の力（トルク）であってもよい。
【０００８】
本発明の一態様では、振動要素を伝わる伝達力に基づいて振動要素間の弾性係数及び減衰係数を推定するパラメータ推定部を備えることが好適である。
【０００９】
本発明の一態様では、未知変数に関わる列ベクトルをｘunknown、観測変数に関わる列ベクトルをｘmeasuredとし、前記運動方程式を行列Ｈ1，Ｈ2を用いて、
【数１】
0006260477.tif 000002
で表した場合に、観測変数取得部は、行列Ｈ1が列フルランクになる条件を満たすように観測変数を取得することが好適である。
【００１０】
本発明の一態様では、観測変数取得部は、各振動要素の加速度を観測変数として取得し、未知変数推定部は、観測変数取得部で取得された各振動要素の加速度に基づいて、振動系の加振力及び各振動要素を伝わる伝達力を未知変数として推定することが好適である。
【００１１】
また、本発明に係る振動解析方法は、複数の振動要素を有する振動系の振動解析方法であって、各振動要素の加速度、振動系の加振力、及び各振動要素を伝わる伝達力の変数のうち、いずれか複数を観測変数として取得する観測変数取得処理と、複数の振動要素間の接続関係を表して振動系をモデル化した場合における各振動要素の運動方程式を用いて、前記変数のうち観測変数以外の未知変数を、観測変数取得処理で取得された観測変数に基づいて推定する未知変数推定処理と、を含むことを要旨とする。
【発明の効果】
【００１２】
本発明によれば、複数の振動要素間の接続関係を表して振動系をモデル化した場合における各振動要素の運動方程式を用いて、観測変数に基づいて未知変数を推定することで、振動要素を伝わる伝達力を推定することができる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施形態に係る振動解析装置の適用対象となる車両の概略構成を示す図である。
【図２】本発明の実施形態に係る振動解析装置の概略構成を示す図である。
【図３】原理モデルの一例を示す図である。
【図４】変速機と駆動輪間の伝達トルクを推定した結果の一例を示す図である。
【図５】原理モデルの他の例を示す図である。
【図６】車体とサブフレーム間の伝達力、車体とパワートレインブロック間の伝達力、及び車体慣性力を算出した結果の一例を示す図である。
【図７】駆動輪と路面間の駆動力、サブフレームと駆動輪間の伝達力、パワートレインブロックとサブフレーム間の伝達力、及びサブフレーム慣性力を算出した結果の一例を示す図である。
【発明を実施するための形態】
【００１４】
以下、本発明を実施するための形態（以下実施形態という）を図面に従って説明する。
【００１５】
図１は本発明の実施形態に係る振動解析装置の適用対象となる車両１０の概略構成を示す図であり、図２は本発明の実施形態に係る振動解析装置の概略構成を示す図である。振動系としての車両１０は、複数の振動要素として車体１１とパワートレインブロック１２（エンジン１３及び変速機１４を含む）と駆動輪１５を含む。エンジン１３は変速機１４に連結され、変速機１４は駆動輪１５に連結され、エンジン１３及び変速機１４はパワートレインブロック１２に含まれている。パワートレインブロック１２はマウント１６を介して車体１１に支持され、駆動輪１５はサスペンションブッシュ１７を介して車体１１に支持されている。車両１０の実稼働状態において、エンジン１３はトルク（回転力）を加振力として発生し、エンジン１３と変速機１４間でトルクが伝達され、変速機１４と駆動輪１５間でトルクが伝達され、エンジン１３及び変速機１４とパワートレインブロック１２間でトルクが伝達される。その際に、パワートレインブロック１２とエンジン１３と変速機１４と駆動輪１５が回転方向に互いに相対的に振動する。さらに、マウント１６を介してパワートレインブロック１２と車体１１間で車両前後方向の力が伝達され、サスペンションブッシュ１７を介して駆動輪１５と車体１１間で車両前後方向の力が伝達される。その際に、車体１１とパワートレインブロック１２と駆動輪１５が車両前後方向に互いに相対的に振動する。
【００１６】
回転方向に振動する振動要素の角加速度を検出するために、パワートレインブロック１２の角加速度ｄ2θb／ｄｔ2を検出するための角加速度センサ２２がパワートレインブロック１２に付設され、エンジン１３の角加速度ｄ2θe／ｄｔ2を検出するための角加速度センサ２３がエンジン１３に付設され、変速機１４の角加速度ｄ2θm／ｄｔ2を検出するための角加速度センサ２４が変速機１４に付設され、駆動輪１５の角加速度ｄ2θw／ｄｔ2を検出するための角加速度センサ２５が駆動輪１５に付設されている。また、車両前後方向に振動する振動要素の車両前後方向並進加速度を検出するために、車体１１の車両前後方向並進加速度ｄ2ｘv／ｄｔ2を検出するための並進加速度センサ３１が車体１１に付設され、パワートレインブロック１２の車両前後方向並進加速度ｄ2ｘb／ｄｔ2を検出するための並進加速度センサ３２がパワートレインブロック１２に付設され、駆動輪１５の車両前後方向並進加速度ｄ2ｘw／ｄｔ2を検出するための並進加速度センサ３５が駆動輪１５に付設されている。車両１０の実稼働状態における、角加速度センサ２２によるパワートレインブロック１２の角加速度ｄ2θb／ｄｔ2、角加速度センサ２３によるエンジン１３の角加速度ｄ2θe／ｄｔ2、角加速度センサ２４による変速機１４の角加速度ｄ2θm／ｄｔ2、角加速度センサ２５による駆動輪１５の角加速度ｄ2θw／ｄｔ2、並進加速度センサ３１による車体１１の車両前後方向並進加速度ｄ2ｘv／ｄｔ2、並進加速度センサ３２によるパワートレインブロック１２の車両前後方向並進加速度ｄ2ｘb／ｄｔ2、及び並進加速度センサ３５による駆動輪１５の車両前後方向並進加速度ｄ2ｘw／ｄｔ2は、振動解析装置７０に入力される。
【００１７】
本実施形態では、振動解析装置７０で車両１０の振動解析を行うために、例えば図３に示すように、複数の振動要素（車体１１とパワートレインブロック１２とエンジン１３と変速機１４と駆動輪１５）間の接続関係を表して振動系（車両１０）をモデル化する。以下、複数の振動要素間の接続関係を表した振動系のモデルを原理モデルとする。原理モデルにおいて、エンジン１３が発生するトルク（振動系の加振力）をτe、エンジン１３と変速機１４間で伝達されるトルクをτm、変速機１４と駆動輪１５間で伝達されるトルクをτds、パワートレインブロック１２に伝達されるトルクをτb、車体１１とパワートレインブロック１２間で伝達される車両前後方向の力をＦb、車体１１と駆動輪１５間で伝達される車両前後方向の力をＦs、駆動輪１５と路面間に作用する車両前後方向の力をＦdとすると、パワートレインブロック１２の回転運動方程式は以下の（１）式で表され、エンジン１３の回転運動方程式は以下の（２）式で表され、変速機１４の回転運動方程式は以下の（３）式で表され、駆動輪１５の回転運動方程式は以下の（４）式で表される。そして、車両前後方向に関して、駆動輪１５の並進運動方程式は以下の（５）式で表され、車体１１の並進運動方程式は以下の（６）式で表され、パワートレインブロック１２の並進運動方程式は以下の（７）式で表される。
【００１８】
【数２】
0006260477.tif 000003
【００１９】
（１）〜（７）式において、Ｊbはパワートレインブロック１２の慣性モーメント、Ｊeはエンジン１３の慣性モーメント、Ｊmは変速機１４の慣性モーメント、Ｊwは駆動輪１５の慣性モーメント、ｍwは駆動輪１５の質量、ｍvは車体１１の質量、ｍbはパワートレインブロック１２の質量、ｄmは変速機１４の変速比、ｒwは駆動輪１５の半径である。ただし、（１）〜（７）式では、時間微分ｄ／ｄｔを・（ドット）で表し（以下の他式も同様）、パワートレインブロック１２の角加速度ｄ2θb／ｄｔ2をθbの上に・・を付して表し、エンジン１３の角加速度ｄ2θe／ｄｔ2をθeの上に・・を付して表し、変速機１４の角加速度ｄ2θm／ｄｔ2をθmの上に・・を付して表し、駆動輪１５の角加速度ｄ2θw／ｄｔ2をθwの上に・・を付して表し、駆動輪１５の車両前後方向並進加速度ｄ2ｘw／ｄｔ2をｘwの上に・・を付して表し、車体１１の車両前後方向並進加速度ｄ2ｘv／ｄｔ2をｘvの上に・・を付して表し、パワートレインブロック１２の車両前後方向並進加速度ｄ2ｘb／ｄｔ2をｘbの上に・・を付して表している。
【００２０】
（１）〜（７）式については、以下の（８）式のように書き換えることができる。（８）式において、Ｍａは各振動要素の慣性力に関わる７行の列ベクトルで以下の（９）式により表され、ｆは各振動要素に作用する力（加振力及び伝達力）に関わる７行の列ベクトルで以下の（１０）式により表され、Ｃは係数に関わる７行７列の行列で以下の（１１）式により表される。
【００２１】
【数３】
0006260477.tif 000004
【００２２】
振動解析装置７０の機能ブロック図の一例を図２に示す。振動解析装置７０は、ＣＰＵを中心としたコンピュータとして構成可能であり、コンピュータを以下に説明する観測変数取得部７１、設定定数記憶装置７２、未知変数推定部７３、パラメータ推定部７４、及び伝達関数演算部７５として機能させる。
【００２３】
振動解析装置７０において、設定定数記憶装置７２には、車体１１の質量ｍv、パワートレインブロック１２の質量ｍb、駆動輪１５の質量ｍw、パワートレインブロック１２の慣性モーメントＪb、エンジン１３の慣性モーメントＪe、変速機１４の慣性モーメントＪm、駆動輪１５の慣性モーメントＪw、変速機１４の変速比ｄm、及び駆動輪１５の半径ｒwが記憶されている。設定定数記憶装置７２に記憶するｍv，ｍb，ｍw，Ｊb，Ｊe，Ｊm，Ｊw，ｄm，ｒwの値については、予め計測しておくか設計値を用いる。（１）〜（７）式におけるｍv，ｍb，ｍw，Ｊb，Ｊe，Ｊm，Ｊw，ｄm，ｒwは、設定定数記憶装置７２に記憶された値から既知となる。
【００２４】
観測変数取得部７１は、（１）〜（７）式における各振動要素の角加速度ｄ2θb／ｄｔ2，ｄ2θe／ｄｔ2，ｄ2θm／ｄｔ2，ｄ2θw／ｄｔ2及び並進加速度ｄ2ｘw／ｄｔ2，ｄ2ｘv／ｄｔ2，ｄ2ｘb／ｄｔ2、振動系の加振力（エンジン１３のトルク）τe、及び加振力τeにより各振動要素を伝わる伝達力τm，τds，τb，Ｆb，Ｆs，Ｆdの変数のうち、いずれか複数を観測変数として取得する。ここでは、角加速度センサ２２，２３，２４，２５及び並進加速度センサ３１，３２，３５による角加速度ｄ2θb／ｄｔ2，ｄ2θe／ｄｔ2，ｄ2θm／ｄｔ2，ｄ2θw／ｄｔ2及び並進加速度ｄ2ｘv／ｄｔ2，ｄ2ｘb／ｄｔ2，ｄ2ｘw／ｄｔ2を観測変数として取得する。
【００２５】
未知変数推定部７３は、（１）〜（７）式の各振動要素の運動方程式を用いて、（１）〜（７）式の変数のうち観測変数以外の未知変数を、観測変数取得部７１で取得された観測変数に基づいて推定する。ここでは、観測変数として各振動要素の角加速度ｄ2θb／ｄｔ2，ｄ2θe／ｄｔ2，ｄ2θm／ｄｔ2，ｄ2θw／ｄｔ2及び並進加速度ｄ2ｘw／ｄｔ2，ｄ2ｘv／ｄｔ2，ｄ2ｘb／ｄｔ2に基づいて、未知変数として振動系の加振力τe及び加振力τeにより各振動要素を伝わる伝達力τm，τds，τb，Ｆb，Ｆs，Ｆdを推定する。
【００２６】
ここで、（８）式を以下の（１２）式のように書き換える。（１２）式において、Ｈは以下の（１３）式で表されるｎ行ｍ列の行列であり、ｘは以下の（１４）式で表されるｍ行の列ベクトルである。
【００２７】
【数４】
0006260477.tif 000005
【００２８】
（１１）式を（１３）式に代入すると、Ｈは以下の（１５）式で表される７行１４列の行列となり、（９）、（１０）式を（１４）式に代入すると、ｘは以下の（１６）式で表される１４行の列ベクトルとなる（ｎ＝７，ｍ＝１４）。（１５）式において、列ベクトルｈ1はτbに対応する係数を表し、列ベクトルｈ2はτeに対応する係数を表し、列ベクトルｈ3はτmに対応する係数を表し、列ベクトルｈ4はτdsに対応する係数を表し、列ベクトルｈ5はＦdに対応する係数を表し、列ベクトルｈ6はＦsに対応する係数を表し、列ベクトルｈ7はＦbに対応する係数を表す。そして、列ベクトルｈ8はＪb×ｄ2θb／ｄｔ2に対応する係数を表し、列ベクトルｈ9はＪe×ｄ2θe／ｄｔ2に対応する係数を表し、列ベクトルｈ10はＪm×ｄ2θm／ｄｔ2に対応する係数を表し、列ベクトルｈ11はＪw×ｄ2θw／ｄｔ2に対応する係数を表し、列ベクトルｈ12はｍw×ｄ2ｘw／ｄｔ2に対応する係数を表し、列ベクトルｈ13はｍv×ｄ2ｘv／ｄｔ2に対応する係数を表し、列ベクトルｈ14はｍb×ｄ2ｘb／ｄｔ2に対応する係数を表す。また、Ｈの１行目は（１）式に対応する係数を表し、Ｈの２行目は（２）式に対応する係数を表し、Ｈの３行目は（３）式に対応する係数を表し、Ｈの４行目は（４）式に対応する係数を表し、Ｈの５行目は（５）式に対応する係数を表し、Ｈの６行目は（６）式に対応する係数を表し、Ｈの７行目は（７）式に対応する係数を表す。
【００２９】
【数５】
0006260477.tif 000006
【００３０】
（１４）式のベクトルｘの成分を未知変数と観測変数に分け、（１２）式を以下の（１７）式のように書き換える。（１７）式において、ｘunknownは未知変数に関わるｍ1行の列ベクトルであり、ｘmeasuredは観測変数に関わるｍ2行の列ベクトルであり、Ｈ1はｎ行ｍ1列の行列で未知変数に対応する係数を表し、Ｈ2はｎ行ｍ2列の行列で観測変数に対応する係数を表し、ｍ＝ｍ1＋ｍ2，ｍ1≦ｎ，ｍ2≦ｎである。
【００３１】
【数６】
0006260477.tif 000007
【００３２】
（１）〜（７）式において、未知変数はτb，τe，τm，τds，Ｆd，Ｆs，Ｆbであり、観測変数はｄ2θb／ｄｔ2，ｄ2θe／ｄｔ2，ｄ2θm／ｄｔ2，ｄ2θw／ｄｔ2，ｄ2ｘw／ｄｔ2，ｄ2ｘv／ｄｔ2，ｄ2ｘb／ｄｔ2であり、ｍ1＝７，ｍ2＝７となる。したがって、ｘunknownは以下の（１８）式で表され、ｘmeasuredは以下の（１９）式で表される。その場合、Ｈ1は以下の（２０）式で表され、Ｈ2は以下の（２１）式で表される。
【００３３】
【数７】
0006260477.tif 000008
【００３４】
（２０）式で表されるＨ1のランクは７でＨ1の列の数ｍ2＝７に等しく、Ｈ1が列フルランクとなる。その場合、ｘunknownは、ｘmeasured，Ｈ1，Ｈ2を用いて以下の（２２）式により算出される。未知変数推定部７３は、観測変数取得部７１で各サンプリング時刻ｋ毎に取得された観測変数ｄ2θb（ｋ）／ｄｔ2，ｄ2θe（ｋ）／ｄｔ2，ｄ2θm（ｋ）／ｄｔ2，ｄ2θw（ｋ）／ｄｔ2，ｄ2ｘw（ｋ）／ｄｔ2，ｄ2ｘv（ｋ）／ｄｔ2，ｄ2ｘb（ｋ）／ｄｔ2に基づいて（２２）式で算出されたｘunknownによって、未知変数τb（ｋ），τe（ｋ），τm（ｋ），τds（ｋ），Ｆd（ｋ），Ｆs（ｋ），Ｆb（ｋ）を各サンプリング時刻ｋ毎に推定することができる。ただし、Ｈ1が正方行列である場合は、（２２）式は、当然、ｘunknown＝−Ｈ1-1Ｈ2ｘmeasuredでよい。
【００３５】
【数８】
0006260477.tif 000009
【００３６】
未知変数推定部７３で変速機１４と駆動輪１５間の伝達トルクτds（ｋ）を推定した結果の一例を図４に示す。図４では、（２２）式で算出された推定トルクτds（ｋ）の時系列波形を、トルクセンサで計測された変速機１４と駆動輪１５間の計測トルクの時系列波形と比較している。図４に示すように、未知変数推定部７３で推定された推定トルクτds（ｋ）がトルクセンサでの計測トルクにほぼ一致していることがわかる。つまり、未知変数推定部７３により、変速機１４と駆動輪１５間の伝達トルクτds（ｋ）を精度よく推定できていることがわかる。
【００３７】
伝達関数演算部７５は、観測変数取得部７１で取得された観測変数及び未知変数推定部７３で推定された未知変数のいずれかに基づいて、振動系の伝達関数を演算する。例えばτe（ｋ）を入力、ｄ2ｘv（ｋ）／ｄｔ2を出力とする場合における入出力の関係を表す伝達関数ｄ2ｘv／ｄｔ2／τeを算出する。また、τe（ｋ）を入力、τds（ｋ）を出力とする場合における入出力の関係を表す伝達関数τds／τeを算出することも可能である。
【００３８】
パラメータ推定部７４は、未知変数推定部７３で推定された振動要素を伝わる伝達力に基づいて、振動要素間の弾性係数及び減衰係数を推定する。例えば変速機１４と駆動輪１５間で伝達されるトルクτds（ｋ）については、変速機１４と駆動輪１５間の回転方向の弾性係数（ばね定数）ｋds及び減衰係数（粘性係数）ｃdsを用いて以下の（２３）式で表すことが可能である。（２３）式において、θm（ｋ）は変速機１４の回転角、θw（ｋ）は駆動輪１５の回転角である。ただし、（２３）式でも、変速機１４の角速度ｄθm（ｋ）／ｄｔをθm（ｋ）の上に・を付して表し、駆動輪１５の角速度ｄθw（ｋ）／ｄｔをθw（ｋ）の上に・を付して表している。（２３）式において、変速機１４と駆動輪１５間の伝達トルクτds（ｋ）については、未知変数推定部７３で推定することができる。変速機１４の角速度ｄθm（ｋ）／ｄｔについては、角加速度センサ２４による変速機１４の角加速度ｄ2θm（ｋ）／ｄｔ2を積分演算して取得することができ、変速機１４の回転角θm（ｋ）については、変速機１４の角速度ｄθm（ｋ）／ｄｔを積分演算して取得することができる。駆動輪１５の角速度ｄθw（ｋ）／ｄｔについては、角加速度センサ２５による駆動輪１５の角加速度ｄ2θw（ｋ）／ｄｔ2を積分演算して取得することができ、駆動輪１５の回転角θw（ｋ）については、駆動輪１５の角速度ｄθw（ｋ）／ｄｔを積分演算して取得することができる。
【００３９】
【数９】
0006260477.tif 000010
【００４０】
（２３）式を各サンプリング時刻ｋ＝１，２，〜，Ｎ毎に行列化すると、以下の（２４）式が得られる。パラメータ推定部７４は、各サンプリング時刻ｋ毎に得られる、τds（ｋ），ｄθm（ｋ）／ｄｔ，ｄθw（ｋ）／ｄｔ，θm（ｋ），θw（ｋ）を用いて、（２４）式における未知パラメータｃds，ｋdsを最小二乗法により同定する。
【００４１】
【数１０】
0006260477.tif 000011
【００４２】
（２４）式を以下の（２５）式のように書き換える。ただし、（２５）式のＡ，ｂは以下の（２６）、（２７）式でそれぞれ表される。
【００４３】
【数１１】
0006260477.tif 000012
【００４４】
（２６）、（２７）式はＮサンプリング点からなる既知の縦ベクトルであり、未知パラメータｃds，ｋdsは以下の（２８）式により算出される。パラメータ推定部７４は、τds（ｋ），ｄθm（ｋ）／ｄｔ，ｄθw（ｋ）／ｄｔ，θm（ｋ），θw（ｋ）に基づく（２８）式により、変速機１４と駆動輪１５間の弾性係数ｋds及び減衰係数ｃdsを推定することができる。
【００４５】
【数１２】
0006260477.tif 000013
【００４６】
また、（２４）式を時間微分すると、以下の（２９）式が得られる。（２９）式を（２５）式のように書き換えると、未知パラメータｃds，ｋdsは以下の（３０）式により算出される。ただし、（３０）式のＡ，ｂは以下の（３１）、（３２）式でそれぞれ表される。パラメータ推定部７４は、ｄτds（ｋ）／ｄｔ，ｄ2θm（ｋ）／ｄｔ2，ｄ2θw（ｋ）／ｄｔ2，ｄθm（ｋ）／ｄｔ，ｄθw（ｋ）／ｄｔに基づく（３０）式によっても、変速機１４と駆動輪１５間の弾性係数ｋds及び減衰係数ｃdsを推定することができる。
【００４７】
【数１３】
0006260477.tif 000014
【００４８】
同様に、エンジン１３と変速機１４間の回転方向の弾性係数ｋm及び減衰係数ｃmについても、（２８）式のｃds，ｋdsをｃm，ｋmにそれぞれ置き換えた式に以下の（３３）、（３４）式を代入するか、（３０）式のｃds，ｋdsをｃm，ｋmにそれぞれ置き換えた式に以下の（３５）、（３６）式を代入することで、推定することができる。（３３）式において、θe（ｋ）はエンジン１３の回転角である。
【００４９】
【数１４】
0006260477.tif 000015
【００５０】
同様に、車体１１とパワートレインブロック１２間の車両前後方向の弾性係数ｋb及び減衰係数ｃbについても、（２８）式のｃds，ｋdsをｃb，ｋbにそれぞれ置き換えた式に以下の（３７）、（３８）式を代入するか、（３０）式のｃds，ｋdsをｃb，ｋbにそれぞれ置き換えた式に以下の（３９）、（４０）式を代入することで、推定することができる。（３７）式において、ｘv（ｋ）は車体１１の車両前後方向変位、ｘb（ｋ）はパワートレインブロック１２の車両前後方向変位である。
【００５１】
【数１５】
0006260477.tif 000016
【００５２】
同様に、車体１１と駆動輪１５間の車両前後方向の弾性係数ｋs及び減衰係数ｃsについても、（２８）式のｃds，ｋdsをｃs，ｋsにそれぞれ置き換えた式に以下の（４１）、（４２）式を代入するか、（３０）式のｃds，ｋdsをｃs，ｋsにそれぞれ置き換えた式に以下の（４３）、（４４）式を代入することで、推定することができる。（４１）式において、ｘw（ｋ）は駆動輪１５の車両前後方向変位である。
【００５３】
【数１６】
0006260477.tif 000017
【００５４】
このようにして弾性係数及び減衰係数をパラメータ推定部７４で推定できるため、（１）〜（７）式については、以下の（４５）式のように書き換えることができる。（４５）式において、Ｍは質量ｍv，ｍb，ｍw及び慣性モーメントＪb，Ｊe，Ｊm，Ｊwに関わる行列であり、Ｃは減衰係数に関わる行列であり、Ｋは弾性係数に関わる行列であり、ｘ，ｕは以下の（４６）式で表される。
【００５５】
【数１７】
0006260477.tif 000018
【００５６】
以上説明した本実施形態では、複数の振動要素間の接続関係だけを表して振動系をモデル化する。そして、各振動要素の加速度、振動系の加振力、及び各振動要素を伝わる伝達力の変数のうち、いずれか複数（ここでは各振動要素の加速度ｄ2θb／ｄｔ2，ｄ2θe／ｄｔ2，ｄ2θm／ｄｔ2，ｄ2θw／ｄｔ2，ｄ2ｘv／ｄｔ2，ｄ2ｘb／ｄｔ2，ｄ2ｘw／ｄｔ2）を観測変数として取得する処理と、（１）〜（７）式の運動方程式を用いて観測変数ｄ2θb／ｄｔ2，ｄ2θe／ｄｔ2，ｄ2θm／ｄｔ2，ｄ2θw／ｄｔ2，ｄ2ｘv／ｄｔ2，ｄ2ｘb／ｄｔ2，ｄ2ｘw／ｄｔ2に基づいて未知変数（ここでは振動系の加振力τe及び各振動要素を伝わる伝達力τm，τds，τb，Ｆb，Ｆs，Ｆd）を推定する処理とをコンピュータに実行させる。その場合、Ｈ1が列フルランクになる条件を満たすため、未知変数（振動系の加振力τe及び各振動要素を伝わる伝達力τm，τds，τb，Ｆb，Ｆs，Ｆd）を推定することができ、振動伝達経路解析を実稼働データから実施することができる。その際には、加速度ｄ2θb／ｄｔ2，ｄ2θe／ｄｔ2，ｄ2θm／ｄｔ2，ｄ2θw／ｄｔ2，ｄ2ｘv／ｄｔ2，ｄ2ｘb／ｄｔ2，ｄ2ｘw／ｄｔ2を観測変数とすることで、加速度に比べて直接的な計測の難しい加振力τe及び伝達力τm，τds，τb，Ｆb，Ｆs，Ｆdの計測が不要となる。また、車体１１を含むパッシブ系とエンジン１３及び変速機１４を含むアクティブ系の全系のシステムを扱うため、アクティブ系を考慮する必要がある低周波数領域でも利用することが可能である。また、特許文献１，２の伝達率や伝達関数のような線形モデルに限定されずに、非線形モデルを扱うことも可能であり、また、定常現象のみならず、車両１０の加速時のような過渡現象も扱うことが可能である。
【００５７】
さらに、本実施形態では、振動要素を伝わる伝達力に基づいて振動要素間の弾性係数及び減衰係数を推定する処理をコンピュータに実行させる。これによって、振動系の弾性係数及び減衰係数のパラメータを求めてパラメトリックモデルを構築することができ、振動系の設計に利用することができる。さらに、本実施形態では、観測変数及び未知変数に基づいて振動系の伝達関数を演算する処理をコンピュータに実行させる。伝達関数は伝達率とは異なりシステム極（共振点）の情報が含まれており、共振点をずらす（例えば上げる）ための弾性係数及び減衰係数の設計に利用することができる。
【００５８】
上記の説明では、観測変数ｄ2θb（ｋ）／ｄｔ2，ｄ2θe（ｋ）／ｄｔ2，ｄ2θm（ｋ）／ｄｔ2，ｄ2θw（ｋ）／ｄｔ2，ｄ2ｘw（ｋ）／ｄｔ2，ｄ2ｘv（ｋ）／ｄｔ2，ｄ2ｘb（ｋ）／ｄｔ2に基づいて、未知変数τb（ｋ），τe（ｋ），τm（ｋ），τds（ｋ），Ｆd（ｋ），Ｆs（ｋ），Ｆb（ｋ）を推定する例について説明した。ただし、観測変数及び未知変数については、Ｈ1が列フルランクとなる条件を満たせば、上記に説明したもの以外であってもよい。例えば観測変数として、ｄ2θb／ｄｔ2，ｄ2θe／ｄｔ2，τe，ｄ2θw／ｄｔ2，ｄ2ｘw／ｄｔ2，ｄ2ｘv／ｄｔ2，ｄ2ｘb／ｄｔ2を観測変数取得部７１で取得することも可能である。その場合、ｘunknownは以下の（４７）式で表され、ｘmeasuredは以下の（４８）式で表され、Ｈ1は以下の（４９）式で表され、Ｈ2は以下の（５０）式で表される。（４９）式で表されるＨ1のランクも７でＨ1の列の数ｍ2＝７に等しく、Ｈ1が列フルランクとなるため、（４８）式のｘmeasured、（４９）式のＨ1、及び（５０）式のＨ2を用いて（２２）式により未知変数τb（ｋ），ｄ2θm／ｄｔ2，τm（ｋ），τds（ｋ），Ｆd（ｋ），Ｆs（ｋ），Ｆb（ｋ）を推定することができる。
【００５９】
【数１８】
0006260477.tif 000019
【００６０】
また、原理モデルの他の例を図５に示す。図５の原理モデルでは、図３の原理モデルと比較して、車体１１とパワートレインブロック１２が連結されているとともに、サブフレーム１８が車体１１、パワートレインブロック１２、及び駆動輪１５と連結されている。図５の原理モデルを用いて、車体１１とサブフレーム１８間の車両前後方向伝達力Ｆsf（ｋ）、車体１１とパワートレインブロック１２間の車両前後方向伝達力Ｆrod（ｋ）、及び車体１１の車両前後方向慣性力を算出した結果を図６に示し、駆動輪１５と路面間の車両前後方向駆動力Ｆd（ｋ）、サブフレーム１８と駆動輪１５間の車両前後方向伝達力Ｆs（ｋ）、パワートレインブロック１２とサブフレーム１８間の車両前後方向伝達力Ｆb（ｋ）、及びサブフレーム１８の車両前後方向慣性力を算出した結果を図７に示す。
【００６１】
以上の実施形態では、振動系が車両１０である例について説明したが、本実施形態は車両以外の振動系に対しても適用可能である。
【００６２】
以上、本発明を実施するための形態について説明したが、本発明はこうした実施形態に何等限定されるものではなく、本発明の要旨を逸脱しない範囲内において、種々なる形態で実施し得ることは勿論である。
【符号の説明】
【００６３】
１０  車両、１１  車体、１２  パワートレインブロック、１３  エンジン、１４  変速機、１５  駆動輪、２２，２３，２４，２５  角加速度センサ、３１，３２，３５  並進加速度センサ、７０  振動解析装置、７１  観測変数取得部、７２  設定定数記憶装置、７３  未知変数推定部、７４  パラメータ推定部、７５  伝達関数演算部。
【図１】
0006260477.tif 000022
【図２】
0006260477.tif 000023
【図３】
0006260477.tif 000024
【図４】
0006260477.tif 000025
【図５】
0006260477.tif 000026
【図６】
0006260477.tif 000027
【図７】
0006260477.tif 000028