(57)【要約】【課題】アクチュエータの内部温度を高精度に推定可能な能動型防振装置を提供する。【解決手段】ＡＣＭ−ＥＣＵ２０６は、エンジン１２の振動情報からアクチュエータ２１０の動作指令値Ｃを算出し、アクチュエータ２１０の内部温度Ｔａｃｍに応じた補正値Ｉｃｒｔ１を算出し、補正値Ｉｃｒｔ１を用いて動作指令値Ｃを補正し、算出対象区間Ｚｔａｒの電圧デューティ比ＤＵＴの平均としての平均デューティ比ＤＵＴａｖｅに基づいて内部温度Ｔａｃｍを推定する。また、ＡＣＭ−ＥＣＵ２０６は、指令値補正部２５４における補正値Ｉｃｒｔ１を調整する。【選択図】図３
2019100382.tif 000002
【特許請求の範囲】
【請求項１】
駆動源と車体との間に配置されたアクチュエータと、
前記アクチュエータに能動的振動を生成させることで前記駆動源から前記車体への振動伝達を抑制するコンピュータと、
を備える能動型防振装置であって、
前記コンピュータは、
前記駆動源の振動情報から前記アクチュエータの動作指令値を算出する指令値算出部と、
前記アクチュエータの内部温度に応じた補正値を算出し、該補正値を用いて前記動作指令値を補正する指令値補正部と、
補正した前記動作指令値に基づく電圧デューティ比を用いて前記アクチュエータに駆動電圧を印加するアクチュエータ駆動制御部と、
所定区間の前記電圧デューティ比の平均としての平均デューティ比に基づいて、前記内部温度を推定する内部温度推定部と、
前記指令値補正部における前記補正値を調整する補正値調整手段と、
を備えることを特徴とする能動型防振装置。
【請求項２】
請求項１に記載の能動型防振装置において、
前記指令値補正部は、補正の実行状態を切り替え可能に構成され、
前記補正値調整手段は、前記指令値補正部により補正の実行状態が切り替わるタイミングで前記補正値を調整する
ことを特徴とする能動型防振装置。
【請求項３】
請求項２に記載の能動型防振装置において、
前記補正値調整手段は、推定した前記内部温度を調整することで、前記補正値の計算に供される計算用温度を算出する内部温度調整部であり、
前記指令値補正部は、前記内部温度推定部からの前記内部温度に代えて、前記内部温度調整部の調整により得た前記計算用温度に応じて前記動作指令値を補正する
ことを特徴とする能動型防振装置。
【請求項４】
請求項３に記載の能動型防振装置において、
前記内部温度調整部は、補正の実行状態が切り替わるタイミングで推定された前記内部温度を、段階的に目標温度に近づける調整を行うことを特徴とする能動型防振装置。
【請求項５】
請求項４に記載の能動型防振装置において、
前記内部温度調整部は、補正の実行状態が切り替わるタイミングで推定された前記内部温度を、時間の経過につれて増加又は減少させながら前記目標温度に近づける調整を行うことを特徴とする能動型防振装置。
【請求項６】
請求項５に記載の能動型防振装置において、
前記内部温度推定部は、演算周期おきに前記内部温度を推定し、
前記内部温度調整部は、補正の実行状態が切り替わるタイミングで推定された前記内部温度を、該内部温度にかかわらず一定の変化量だけ、１つ又は複数の演算周期にわたって増加又は減少させながら前記目標温度に近づける調整を行う
ことを特徴とする能動型防振装置。
【請求項７】
請求項４〜６のいずれか１項に記載の能動型防振装置において、
前記内部温度調整部は、補正のオン状態からオフ状態に切り替わるタイミングで推定された前記内部温度を、前記補正値がゼロと算出される前記目標温度に近づける調整を行うことを特徴とする能動型防振装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、駆動源から車体への振動の伝達を抑制する能動型防振装置に関する。
【背景技術】
【０００２】
特許文献１では、エンジンマウントのボイスコイル３１に流れている制御電流値に基づいてマウント本体３０の温度Ｔ（又はボイスコイル３１の温度）を算出する装置が提案されている（要約、［００１４］）。この際、マウント本体３０の温度が上昇してボイスコイル３１の温度が上昇すると、ボイスコイル３１の抵抗値が上昇するので、制御電流値が減少することを利用する旨の記載がある（［００１４］）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平０８−１７７９６５号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ところが、特許文献１のように、ボイスコイル３１の温度変化に応じた抵抗値の変化を利用する場合、温度に対する抵抗値の関係性のばらつきにより、温度の推定精度が低下するおそれがある。
【０００５】
本発明は上記した問題を解決するためになされたものであり、アクチュエータの内部温度を高精度に推定可能な能動型防振装置を提供することを目的とする。
【課題を解決するための手段】
【０００６】
本発明に係る能動型防振装置は、駆動源と車体との間に配置されたアクチュエータと、前記アクチュエータに能動的振動を生成させることで前記駆動源から前記車体への振動伝達を抑制するコンピュータと、を備える能動型防振装置であって、前記コンピュータは、前記駆動源の振動情報から前記アクチュエータの動作指令値を算出する指令値算出部と、前記アクチュエータの内部温度に応じた補正値を算出し、該補正値を用いて前記動作指令値を補正する指令値補正部と、補正した前記動作指令値に基づく電圧デューティ比を用いて前記アクチュエータに駆動電圧を印加するアクチュエータ駆動制御部と、所定区間の前記電圧デューティ比の平均としての平均デューティ比に基づいて、前記内部温度を推定する内部温度推定部と、前記指令値補正部における前記補正値を調整する補正値調整手段と、を備える。
【０００７】
アクチュエータの内部温度に対する平均デューティ比の特性曲線には、高い相関性及び高い再現性が部分的にみられる。この特性曲線の傾向を考慮して、所定区間の平均デューティ比に基づいて動作指令値を補正することで、アクチュエータの内部温度を高精度に推定することができる。しかも、動作指令値の補正に供される補正値を調整することで、状況に応じて発生する補正の不具合を予防することができる。
【０００８】
また、前記指令値補正部は、補正の実行状態を切り替え可能に構成され、前記補正値調整手段は、前記指令値補正部により補正の実行状態が切り替わるタイミングで前記補正値を調整してもよい。一般的に言えば、補正の実行状態が切り替わるタイミングにて補正の不具合が発生する可能性が高いので、本発明による予防効果が顕著に現われる。
【０００９】
また、前記補正値調整手段は、推定した前記内部温度を調整することで、前記補正値の計算に供される計算用温度を算出する内部温度調整部であり、前記指令値補正部は、前記内部温度推定部からの前記内部温度に代えて、前記内部温度調整部の調整により得た前記計算用温度に応じて前記動作指令値を補正してもよい。これにより、指令値補正部が実行する演算処理を変更せずに計算用温度のみを変更することで、間接的に補正値を調整することができる。
【００１０】
また、前記内部温度調整部は、補正の実行状態が切り替わるタイミングで推定された前記内部温度を、段階的に目標温度に近づける調整を行ってもよい。これにより、補正の実行状態が切り替わるタイミングで発生する制御の不連続性を抑制可能となり、例えば過補正に起因する防振装置の性能低下を防止することができる。
【００１１】
また、前記内部温度調整部は、補正の実行状態が切り替わるタイミングで推定された前記内部温度を、時間の経過につれて増加又は減少させながら前記目標温度に近づける調整を行ってもよい。これにより、円滑かつ緩やかに補正値を調整することができる。
【００１２】
また、前記内部温度推定部は、演算周期おきに前記内部温度を推定し、前記内部温度調整部は、補正の実行状態が切り替わるタイミングで推定された前記内部温度を、該内部温度にかかわらず一定の変化量だけ、１つ又は複数の演算周期にわたって増加又は減少させながら前記目標温度に近づける調整を行ってもよい。これにより、内部温度と目標温度の乖離にかかわらず、略同じ変化度合いで補正値を調整することができる。
【００１３】
また、前記内部温度調整部は、補正のオン状態からオフ状態に切り替わるタイミングで推定された前記内部温度を、前記補正値がゼロと算出される前記目標温度に近づける調整を行ってもよい。これにより、内部温度の調整を通じて、補正をオン状態からオフ状態に円滑に移行させることができる。
【発明の効果】
【００１４】
本発明に係る能動型防振装置によれば、アクチュエータの内部温度を高精度に推定することができる。
【図面の簡単な説明】
【００１５】
【図１】本発明の一実施形態における能動型防振装置を搭載した車両の概略構成図である。
【図２】図１のエンジンマウントの内部構成を示す図である。
【図３】図１及び図２のＡＣＭ−ＥＣＵの演算部の機能を示すブロック図である。
【図４】図３のエンジン状態変数算出部の機能を示す詳細なブロック図である。
【図５】図３の指令値算出部の機能を示すブロック図である。
【図６】図５の振動周期ポイント数算出部が用いるポイント数マップを示す図である。
【図７】図３の指令値補正部の機能を示す詳細なブロック図である。
【図８】図３の目標電流波形算出部の機能を示す詳細なブロック図である。
【図９】図８の１次目標電流波形算出部が用いる目標電流波形マップを示す図である。
【図１０】図１０Ａ〜図１０Ｃは、図８の２次目標電流波形算出部により２次目標電流波形を算出する第１、第２、第３状態を示す図である。
【図１１】図１１Ａは、１次目標電流波形算出部が算出した１次目標電流波形を示す図である。図１１Ｂは、２次目標電流波形算出部が算出した２次目標電流波形を示す図である。図１１Ｃは、波形合成部が１次目標電流波形と２次目標電流波形を合成して算出した目標電流波形を示す図である。
【図１２】図３のアクチュエータ駆動制御部の機能を示す詳細なブロック図である。
【図１３】アクチュエータ駆動制御部が出力する駆動電圧の一例を示す図である。
【図１４】図１２の駆動回路の一部を示す回路図である。
【図１５】図１〜図３のアクチュエータの目標電流波形及び駆動電流の波形の例を示す図である。
【図１６】低温時及び高温時における目標電流波形と駆動電流との関係の一例を示す図である。
【図１７】本実施形態における温度推定の原理を説明するための図である。
【図１８】図１８Ａは、目標電流配列の一例を示す。図１８Ｂは、図１８Ａの目標電流配列を用いた場合の電圧デューティ比を示す図である。図１８Ｃは、図１８Ｂの一部を拡大して示す図である。
【図１９】図３の内部温度推定部の機能を示す詳細なブロック図である。
【図２０】図１９の基準平均デューティ比算出部が用いる基準平均デューティ比マップを示す図である。
【図２１】ポイント数と平均デューティ比との関係の一例を示す図である。
【図２２】図１９の平均デューティ比変化率算出部が用いる平均デューティ比変化率マップを示す図である。
【図２３】比較例における温度補正の実行状態の遷移を示すタイムチャートである。
【図２４】図３に示す内部温度調整部の動作説明に供されるフローチャートである。
【図２５】図２４のフローチャートに基づく内部温度の調整結果を模式的に示す図である。
【図２６】実施例における温度補正の実行状態の遷移を示すタイムチャートである。
【発明を実施するための形態】
【００１６】
以下、本発明における能動型防振装置について好適な実施形態を挙げ、添付の図面を参照しながら説明する。
【００１７】
［車両１０の概略構成］
図１は、本発明の一実施形態における能動型防振装置２００を搭載した車両１０の概略構成図である。本図に示す車両１０は、駆動源（原動機）としてのエンジン１２を有するエンジン車両である。この車両１０は、エンジン１２に加えて、走行モータを有するハイブリッド車両であってもよい。
【００１８】
エンジン１２は、その回転軸が車幅方向とされた状態において、エンジンマウント２０２ｆ、２０２ｒを介して車体１４に支持されている。エンジンマウント２０２ｆ、２０２ｒは、アクチュエータ２１０を駆動することによりエンジン１２からの振動（以下、エンジン振動）を能動的に抑制する能動型防振装置２００の一部を構成する。
【００１９】
車両１０は、能動型防振装置２００に加え、エンジン１２の制御に関連するエンジン制御系１００と、バッテリ１６とを有する。なお、車両１０の基本的な構成要素については、例えば、上記した特許文献１（特開平０８−１７７９６５号公報）、特開２０１１−２５２５５３号公報又は特開２０１４−１３７００３号公報と同様のものを用いることができる。
【００２０】
エンジン制御系１００は、エンジン１２に関連する構成要素として、クランクセンサ１０２と、上死点センサ（以下、ＴＤＣセンサ１０４）と、スタータモータ１０６と、外気温センサ１０８と、燃料噴射電子制御装置（以下、ＦＩ−ＥＣＵ１１０）と、を含んで構成される。
【００２１】
クランクセンサ１０２は、図示しないクランクシャフトの回転位置（以下、クランク回転位置θｃｒｋ）を検出し、クランク回転位置θｃｒｋを示す信号（以下、ＣＲＫパルス信号Ｓｃｒｋ）をＦＩ−ＥＣＵ１１０に出力する。ＴＤＣセンサ１０４は、図示しないエンジンピストンが上死点に来たタイミング（以下、上死点タイミング）を検出し、上死点タイミングを示す信号（以下、ＴＤＣパルス信号Ｓｔｄｃ）をＦＩ−ＥＣＵ１１０に出力する。
【００２２】
スタータモータ１０６は、エンジン１２のモータリングに用いられ、バッテリ１６からの電力に基づいてエンジン１２に対して駆動力を伝達する。外気温センサ１０８は、外気温Ｔｅｘ（単位：℃）を検出してＦＩ−ＥＣＵ１１０及び後述するＡＣＭ−ＥＣＵ２０６にそれぞれ出力する。
【００２３】
ＦＩ−ＥＣＵ１１０は、ＣＲＫパルス信号Ｓｃｒｋ、ＴＤＣパルス信号Ｓｔｄｃ等の各種入力信号に基づいてエンジン１２を制御する。例えば、ＦＩ−ＥＣＵ１１０は、ＣＲＫパルス信号Ｓｃｒｋに基づいて、単位時間当たりのエンジン１２の回転数（以下、エンジン回転速度Ｎｅ；単位はｒｐｍ）を算出して用いる。後述する能動型防振装置２００のＡＣＭ電子制御装置（以下、ＡＣＭ−ＥＣＵ２０６）と同様、ＦＩ−ＥＣＵ１１０は、図示しない入出力部、演算部及び記憶部を有する。
【００２４】
ＦＩ−ＥＣＵ１１０は、エンジン１２の気筒休止モード（気筒作動モード）を制御可能に構成される。この場合、ＦＩ−ＥＣＵ１１０は、ＣＲＫパルス信号Ｓｃｒｋと、ＴＤＣパルス信号Ｓｔｄｃと、気筒の作動状態（全気筒作動モード又は気筒休止モード）を示す気筒信号ＳｃｙとをＡＣＭ−ＥＣＵ２０６に送信する。
【００２５】
［能動型防振装置２００の構成］
能動型防振装置２００は、エンジン１２から車体１４への振動伝達を抑制するための装置である。この能動型防振装置２００は、上述のエンジンマウント２０２ｆ、２０２ｒ及びＡＣＭ−ＥＣＵ２０６に加え、電流センサ２０４ｆ、２０４ｒを有する。
【００２６】
電流センサ２０４ｆ、２０４ｒは、ＡＣＭ−ＥＣＵ２０６を介してバッテリ１６からエンジンマウント２０２ｆ、２０２ｒに供給される電流（以下、駆動電流Ｉｄｆ、Ｉｄｒ、或いは総称して駆動電流Ｉｄ）を検出して、ＡＣＭ−ＥＣＵ２０６に出力する。
【００２７】
図２は、図１のエンジンマウント２０２ｆ、２０２ｒの内部構成を示す図である。エンジンマウント２０２ｆ、２０２ｒは、エンジン振動を能動的に抑制するいわゆるＡＣＭ（Active Control Mount）として機能する。具体的には、エンジンマウント２０２ｆ、２０２ｒはそれぞれ、アクチュエータ２１０、加振板２１２及びゴム板２１４を有する。
【００２８】
アクチュエータ２１０は、駆動軸２１６及びコイル２１８を有し、エンジン振動を相殺する相殺振動を生成する。駆動軸２１６は、コイル２１８の通電に伴う電磁力に応じて進退する。加振板２１２は、駆動軸２１６の進退に応じて進退して、エンジンマウント２０２ｆ、２０２ｒ内に封入された液体を付勢する。ゴム板２１４は、固定されている加振板２１２の動きに合わせて一体的に変位する。
【００２９】
ＡＣＭ−ＥＣＵ２０６は、アクチュエータ２１０を駆動制御する装置であり、入出力部２２０、演算部２２２及び記憶部２２４（図１）を有する。ＡＣＭ−ＥＣＵ２０６がアクチュエータ２１０を駆動させることにより、車体１４へのエンジン振動の伝達を抑制するための振動抑制制御を行う。
【００３０】
入出力部２２０は、ＡＣＭ−ＥＣＵ２０６と他の各部との間の信号の入出力を行う。演算部２２２は、記憶部２２４に記憶されているプログラムを実行することにより、エンジンマウント２０２ｆ、２０２ｒを制御する装置であり、例えば、中央処理装置（ＣＰＵ）を含む。演算部２２２の詳細については、図３を参照して後述する。
【００３１】
記憶部２２４は、演算部２２２が利用するプログラム及びデータを記憶する。記憶部２２４は、揮発性メモリ又は不揮発性メモリからなるＲＡＭ（Random Access Memory）を備える。また、記憶部２２４は、ＲＡＭに加え、ＲＯＭ（Read Only Memory）を有してもよい。
【００３２】
なお、本実施形態では、演算部２２２が用いるプログラム及びデータは、車両１０の記憶部２２４に記憶されていることを想定している。しかしながら、例えば、入出力部２２０に含まれる無線装置を介して外部サーバ（いずれも不図示）からプログラム及びデータの一部を取得してもよい。
【００３３】
［演算部２２２の機能］
＜全体概要＞
図３は、図１及び図２のＡＣＭ−ＥＣＵ２０６の演算部２２２の機能を示すブロック図である。図３に示す演算部２２２の各部は、記憶部２２４に記憶されたプログラムの実行により実現される機能部である。
【００３４】
この演算部２２２は、エンジン状態変数算出部２５０と、指令値算出部２５２と、指令値補正部２５４と、目標電流波形算出部２５６と、アクチュエータ駆動制御部２５８と、内部温度推定部２６０と、内部温度調整部２６２と、補正実行判断部２６４を有する。
【００３５】
エンジン状態変数算出部２５０は、エンジン状態変数として、エンジン１２が振動する振幅（以下、振動振幅Ａｖ）及びエンジン１２が回転・振動する周期（以下、振動周期Ｐｖ）を算出する。なお、エンジン状態変数算出部２５０は、振動周期Ｐｖに代えて周波数を算出してもよい。
【００３６】
指令値算出部２５２は、エンジン状態変数算出部２５０からの振動振幅Ａｖ及び振動周期Ｐｖに基づいて、アクチュエータ２１０の動作指令値Ｃを算出する。ここで、指令値算出部２５２は、補正実行判断部２６４からの補正モードフラグＦｌｇｍｏｄに応じて、後述する２種類の算出モードを切り替えながら動作指令値Ｃを算出することができる。
【００３７】
指令値補正部２５４は、指令値算出部２５２からの動作指令値Ｃと、外気温センサ１０８からの外気温Ｔｅｘと、内部温度調整部２６２からの計算用温度Ｔｃａｌとに基づいて動作指令値Ｃを補正する。ここで、指令値補正部２５４は、補正実行判断部２６４からの補正モードフラグＦｌｇｍｏｄに応じて、補正の実行状態（ＯＮ／ＯＦＦ）を切り替えることができる。
【００３８】
目標電流波形算出部２５６は、指令値補正部２５４が補正した動作指令値Ｃに基づいて目標電流波形Ｗｉを算出する。アクチュエータ駆動制御部２５８は、目標電流波形Ｗｉに基づいて駆動電圧Ｖｄを印加することで、アクチュエータ２１０を駆動する。
【００３９】
内部温度推定部２６０は、アクチュエータ駆動制御部２５８からの電圧デューティ比ＤＵＴに基づいてアクチュエータ２１０の内部温度Ｔａｃｍを推定する。内部温度調整部２６２は、内部温度推定部２６０が推定した内部温度Ｔａｃｍを必要に応じて調整することで、内部温度Ｔａｃｍと同一の又は異なる温度（以下、計算用温度Ｔｃａｌ）を算出する。
【００４０】
補正実行判断部２６４は、ＡＣＭ−ＥＣＵ２０６の外部情報又は内部情報に基づいて、指令値補正部２５４による動作指令値Ｃの補正の実行要否を判断する。外部情報は、ＡＣＭ−ＥＣＵ２０６に入力される各種情報であり、例えば、外気温Ｔｅｘ、駆動電流Ｉｄを含む。内部情報は、ＡＣＭ−ＥＣＵ２０６が生成する各種情報であり、例えば、内部温度Ｔａｃｍ、振動周期Ｐｖ、ポイント数Ｎｐ、補正モードフラグＦｌｇｍｏｄを含む。
【００４１】
＜エンジン状態変数算出部２５０の詳細＞
図４は、図３のエンジン状態変数算出部２５０の機能を示す詳細なブロック図である。このエンジン状態変数算出部２５０は、エンジン振動振幅算出部２７０と、エンジン振動周期算出部２７２と、を含んで構成される。
【００４２】
エンジン振動振幅算出部２７０は、ＣＲＫパルス信号Ｓｃｒｋ及びＴＤＣパルス信号Ｓｔｄｃに基づいて振動振幅Ａｖを算出する。具体的には、エンジン振動振幅算出部２７０は、ＣＲＫパルス信号Ｓｃｒｋの１周期における特定時点（例えば、立ち上がり時点又は立ち下がり時点）の間隔Ｄｃｒｋの累積値∫Ｄｃｒｋを算出する。次いで、エンジン振動振幅算出部２７０は、累積値∫Ｄｃｒｋを上死点タイミングの間隔Ｄｔｄｃで割って平均値Ｄｃｒｋａｖｅ（＝（∫Ｄｃｒｋ）／Ｄｔｄｃ）を算出する。
【００４３】
そして、エンジン振動振幅算出部２７０は、各特定時点の間隔Ｄｃｒｋと平均値Ｄｃｒｋａｖｅの偏差ΔＤｃｒｋを算出し、偏差ΔＤｃｒｋの最大値と最小値の差を振動振幅Ａｖとする。振動振幅Ａｖは、実際のエンジン振動振幅と高い相関性を持つため、エンジン振動振幅の代わりとして用いられる。
【００４４】
なお、エンジン振動振幅算出部２７０中の括弧内には、演算処理の入出力特性が模式的に示されている。具体的には、エンジン振動振幅算出部２７０がＣＲＫパルス信号Ｓｃｒｋ及びＴＤＣパルス信号Ｓｔｄｃを主たる入力値として振動振幅Ａｖを算出することを示している。「Ｇ」は関数を意味している。他の構成要素についても同様である。
【００４５】
エンジン振動周期算出部２７２は、ＣＲＫパルス信号Ｓｃｒｋ及び気筒信号Ｓｃｙに基づいて振動周期Ｐｖを算出する。具体的には、エンジン振動周期算出部２７２は、ＣＲＫパルス信号Ｓｃｒｋの１周期における特定時点の間隔を振動周期Ｐｖとして算出する。なお、気筒信号Ｓｃｙは、気筒の作動状態に応じて振動周期Ｐｖの変化を判定するために用いられる。
【００４６】
＜指令値算出部２５２の詳細＞
図５は、図３の指令値算出部２５２の機能を示す詳細なブロック図である。この指令値算出部２５２は、１次指令電流算出部２８０と、１次位相算出部２８２と、２次指令電流算出部２８４と、２次位相算出部２８６と、振動周期ポイント数算出部２８８と、２次電流乗数設定部２９０と、乗算器２９２と、を含んで構成される。
【００４７】
１次指令電流算出部２８０は、振動振幅Ａｖと振動周期Ｐｖの組合せ毎に１次指令電流Ｉｃ１（単位：Ａ）を算出する。ここで、１次指令電流Ｉｃ１は、アクチュエータ２１０に供給する電流の指令値の１次成分である。より具体的には、１次指令電流算出部２８０は、振動振幅Ａｖと振動周期Ｐｖの組合せと１次指令電流Ｉｃ１の関係を規定したマップ（１次指令電流マップ）を有する。そして、１次指令電流算出部２８０は、この１次指令電流マップを参照することで、振動振幅Ａｖと振動周期Ｐｖの組合せに対応する１次指令電流Ｉｃ１を読み出して指令値補正部２５４に出力する。
【００４８】
１次位相算出部２８２は、振動振幅Ａｖと振動周期Ｐｖの組合せ毎に１次位相Ｐｃ１（単位；ｒａｄ）を算出する。ここで、１次位相Ｐｃ１は、１次指令電流Ｉｃ１の位相である。より具体的には、１次位相算出部２８２は、振動振幅Ａｖと振動周期Ｐｖの組合せと１次位相Ｐｃ１の関係を規定したマップ（１次位相マップ）を有する。そして、１次位相算出部２８２は、この１次位相マップを参照することで、振動振幅Ａｖと振動周期Ｐｖの組合せに対応する１次位相Ｐｃ１を読み出して指令値補正部２５４に出力する。
【００４９】
２次指令電流算出部２８４は、振動振幅Ａｖと振動周期Ｐｖの組合せ毎に２次指令電流Ｉｃ２（単位：Ａ）を算出する。ここで、２次指令電流Ｉｃ２は、アクチュエータ２１０に供給する電流の指令値の２次成分である。より具体的には、２次指令電流算出部２８４は、振動振幅Ａｖと振動周期Ｐｖの組合せと２次指令電流Ｉｃ２の関係を規定したマップ（２次指令電流マップ）を有する。そして、２次指令電流算出部２８４は、この２次指令電流マップを参照することで、振動振幅Ａｖと振動周期Ｐｖの組合せに対応する２次指令電流Ｉｃ２を読み出して指令値補正部２５４に出力する。
【００５０】
２次位相算出部２８６は、振動振幅Ａｖと振動周期Ｐｖの組合せ毎に２次位相Ｐｃ２（単位：ｒａｄ）を算出する。ここで、２次位相Ｐｃ２は、２次指令電流Ｉｃ２の位相である。より具体的には、２次位相算出部２８６は、振動振幅Ａｖと振動周期Ｐｖの組合せと２次位相Ｐｃ２の関係を規定したマップ（２次位相マップ）を有する。そして、２次位相算出部２８６は、この２次位相マップを参照することで、振動振幅Ａｖと振動周期Ｐｖの組合せに対応する２次位相Ｐｃ２を読み出して指令値補正部２５４に出力する。
【００５１】
振動周期ポイント数算出部２８８は、振動周期Ｐｖに応じて目標電流Ｉｔａｒのポイント数Ｎｐ（単位：個）を算出する。ここで、ポイント数Ｎｐは、１波形当たりの目標電流Ｉｔａｒのプロット数であり、振動周期Ｐｖに比例して増加する。より詳しくは、振動周期ポイント数算出部２８８は、振動周期Ｐｖとポイント数Ｎｐの関係を規定したマップ（以下、ポイント数マップ３００）を有する。そして、振動周期ポイント数算出部２８８は、図６に示すポイント数マップ３００を参照することで、振動周期Ｐｖに対応するポイント数Ｎｐを読み出して指令値補正部２５４に出力する。
【００５２】
２次電流乗数設定部２９０は、補正モードフラグＦｌｇｍｏｄに応じて、２次指令電流算出部２８４の出力側に設けられた乗算器２９２の乗数を設定する。例えば、２次電流乗数設定部２９０は、振動周期Ｐｖに対応する１次成分からなる動作指令値Ｃを算出する「第１算出モード」の場合には乗数を０に設定する。一方、振動周期Ｐｖに対応する１次及び複数次成分（ここでは、２次）からなる動作指令値Ｃを算出する「第２算出モード」の場合には乗数を正値（例えば、１）に設定する。
【００５３】
＜指令値補正部２５４の詳細＞
図７は、図３の指令値補正部２５４の機能を示す詳細なブロック図である。ここでは、指令値補正部２５４は、内部温度Ｔａｃｍに応じて１次指令電流Ｉｃ１を補正する。なぜならば、後述する目標電流波形Ｗｉが同じであっても、内部温度Ｔａｃｍの変化に応じて、アクチュエータ２１０による発生力Ｆが異なるためである。
【００５４】
指令値補正部２５４は、フラグ判定部３１０と、基準発生力算出部３１２と、発生力低下率算出部３１４と、発生力低下量算出部３１６と、補正値算出部３１８と、１次指令電流補正部３２０とを含んで構成される。
【００５５】
フラグ判定部３１０は、補正モードフラグＦｌｇｍｏｄの値を参照し、温度補正の実行状態を判定する。フラグ判定部３１０は、温度補正が「ＯＮ」である場合に、１次指令電流Ｉｃ１を基準発生力算出部３１２に出力する。一方、フラグ判定部３１０は、温度補正が「ＯＦＦ」である場合に、１次指令電流Ｉｃ１をそのまま目標電流波形算出部２５６に出力する。
【００５６】
基準発生力算出部３１２は、１次指令電流Ｉｃ１に基づいてアクチュエータ２１０の発生力Ｆの基準値（以下、基準発生力Ｆｒｅｆ）を算出する。
【００５７】
発生力低下率算出部３１４は、内部温度調整部２６２からの計算用温度Ｔｃａｌと、基準温度Ｔｒｅｆとの関係から、アクチュエータ２１０の発生力Ｆの低下度合いを示す発生力低下率Ｒｆｄｗｎを算出する。
【００５８】
発生力低下量算出部３１６は、基準発生力算出部３１２からの基準発生力Ｆｒｅｆと、発生力低下率算出部３１４からの発生力低下率Ｒｆｄｗｎとの積を、発生力低下量Ｆｄｗｎとして算出する（Ｆｄｗｎ＝Ｆｒｅｆ・Ｒｆｄｗｎ）。
【００５９】
補正値算出部３１８は、発生力低下量算出部３１６からの発生力低下量Ｆｄｗｎを電流に換算した後、１次指令電流Ｉｃ１に対する補正値Ｉｃｒｔ１を算出する。
【００６０】
１次指令電流補正部３２０は、温度補正が「ＯＮ」である場合、１次指令電流Ｉｃ１に対して、補正値Ｉｃｒｔ１を加えることで１次指令電流Ｉｃ１を補正する（Ｉｃ１（補正後）＝Ｉｃ１（補正前）＋Ｉｃｒｔ１）。
【００６１】
＜目標電流波形算出部２５６の詳細＞
図８は、図３の目標電流波形算出部２５６の機能を示す詳細なブロック図である。この目標電流波形算出部２５６は、１次目標電流波形算出部３３０と、２次目標電流波形算出部３３２と、波形合成部３３４と、ゼロ点補正部３３６と、を含んで構成される。
【００６２】
１次目標電流波形算出部３３０は、振動周期Ｐｖ、１次指令電流Ｉｃ１及び１次位相Ｐｃ１に基づいて１次目標電流波形Ｗｉ１（単位：Ａ）を算出する。より具体的には、１次目標電流波形算出部３３０は、図９の目標電流波形マップ３５０を参照することで、振動周期Ｐｖと１次指令電流Ｉｃ１の組合せに対応する目標電流配列Ｍｉを読み出す。そして、１次目標電流波形算出部３３０は、読み出した目標電流配列Ｍｉに１次位相Ｐｃ１を反映させて１次目標電流波形Ｗｉ１を算出する。
【００６３】
図９は、図８の１次目標電流波形算出部３３０が用いる目標電流波形マップ３５０を示す図である。目標電流波形マップ３５０は、振動周期Ｐｖと１次指令電流Ｉｃ１の組合せ毎に設けられた目標電流配列Ｍｉの集合体に相当する。各々の目標電流配列Ｍｉは、目標電流Ｉｔａｒの時系列を示すデータセットであり、目標電流波形Ｗｉを構成する。
【００６４】
２次目標電流波形算出部３３２は、振動周期Ｐｖ、２次指令電流Ｉｃ２及び２次位相Ｐｃ２に基づいて２次目標電流波形Ｗｉ２（単位：Ａ）を算出する。より具体的には、２次目標電流波形算出部３３２は、図９のマップと同一の又は異なる目標電流波形マップ３５０を参照することで、振動周期Ｐｖと２次指令電流Ｉｃ２の組合せに対応する目標電流配列Ｍｉを読み出す。なお、同一の目標電流波形マップ３５０を用いる場合、振動周期Ｐｖの半分に相当する周期と、２次指令電流Ｉｃ２の組合せに対応する目標電流配列Ｍｉを読み出す点に留意する。
【００６５】
図１０Ａ〜図１０Ｃは、図８の２次目標電流波形算出部３３２により２次目標電流波形Ｗｉ２を算出する第１、第２、第３状態を示す図である。
【００６６】
図１０Ａには、振動周期Ｐｖ及び２次指令電流Ｉｃ２に基づいて算出された２次目標電流波形Ｗｉ２が示されている。この２次目標電流波形Ｗｉ２は、１回の波形演算周期Ｐｗｃ内において、読み出された目標電流配列Ｍｉ（振動周期Ｐｖの半周期）を２回分だけ順次配置されてなる点に留意する。
【００６７】
図１０Ｂには、図１０Ａの２次目標電流波形Ｗｉ２に対して２次位相Ｐｃ２を反映した状態が示されている。この２次目標電流波形Ｗｉ２は、図１０Ａに示す波形に対して、２次位相Ｐｃ２だけ遅れるように調整されている。これにより、図１０Ｂの２次目標電流波形Ｗｉ２を構成する目標電流Ｉｔａｒのうち白丸で示した部分（右側）は、今回の波形演算周期Ｐｗｃからはみ出す。
【００６８】
図１０Ｃには、図１０Ｂではみ出した目標電流Ｉｔａｒ（白丸）を、同じ波形演算周期Ｐｗｃの開始位置に平行移動させた状態が示されている。これにより、目標電流Ｉｔａｒの周期性が保証される。
【００６９】
波形合成部３３４は、１次目標電流波形Ｗｉ１と２次目標電流波形Ｗｉ２を合成することで目標電流波形Ｗｉ（合成目標電流波形）を算出する。以下、具体的な合成結果について、図１１Ａ〜図１１Ｃを参照しながら説明する。
【００７０】
図１１Ａは、１次目標電流波形算出部３３０が算出した１次目標電流波形Ｗｉ１を示す。図１１Ｂは、２次目標電流波形算出部３３２が算出した２次目標電流波形Ｗｉ２を示す。図１１Ｃは、波形合成部３３４が１次目標電流波形Ｗｉ１と２次目標電流波形Ｗｉ２を合成して算出した目標電流波形Ｗｉを示す。
【００７１】
図１１Ｃの目標電流波形Ｗｉに関して、１つの波形演算周期Ｐｗｃは、電流のボトム値からピーク値までの区間（以下、電流上がり区間Ｚｕｐ）と、電流のピーク値からボトム値までの区間（以下、電流下り区間Ｚｄｗｎ）と、に区分される。
【００７２】
ゼロ点補正部３３６は、目標電流波形算出部２５６が算出した目標電流波形Ｗｉに対してゼロ点補正を行う。図１１Ｃに示すように、波形演算周期Ｐｗｃの全区間において、目標電流波形Ｗｉがゼロよりも大きくなる場合がある。ゼロ点補正は、目標電流波形Ｗｉのボトム値（最小値）がゼロになるように目標電流波形Ｗｉ全体を下げる補正に相当する。ゼロ点補正を行った目標電流波形Ｗｉは、後述する図１８Ａに示されている。
【００７３】
＜アクチュエータ駆動制御部２５８の概要＞
図１２は、図３のアクチュエータ駆動制御部２５８の詳細を示すブロック図である。このアクチュエータ駆動制御部２５８は、次回目標電流値算出部３６０と、電圧デューティ比算出部３６２と、駆動回路３６４と、を含んで構成される。本図では、１組の電流センサ２０４及びアクチュエータ２１０を表記しているが、実際には、図１に示すように複数組設けられる。
【００７４】
次回目標電流値算出部３６０は、目標電流波形算出部２５６で算出された目標電流波形Ｗｉ（図１１Ｃ）から順次選択した目標電流Ｉｔａｒを、次回目標電流Ｉｔａｒｎｘｔとして算出する。
【００７５】
電圧デューティ比算出部３６２は、次回目標電流値算出部３６０からの次回目標電流値Ｉｔａｒｎｘｔと、電流センサ２０４からの駆動電流Ｉｄに基づいて、電圧デューティ比ＤＵＴを算出する。電圧デューティ比ＤＵＴは、スイッチング周期Ｐｓｗに対する駆動電圧Ｖｄの印加時間Ｔａの割合であり、下記の式（１）により定義される。
ＤＵＴ＝Ｔａ／Ｐｓｗ      （１）
【００７６】
図１３は、アクチュエータ駆動制御部２５８が出力する駆動電圧Ｖｄの一例を示す図である。スイッチング周期Ｐｓｗは、波形演算周期Ｐｗｃの整数倍に相当する。また、駆動電圧Ｖｄは、時間によらず変化しない固定電圧である。
【００７７】
電圧デューティ比算出部３６２は、例えばＰＩＤ制御を用いる場合、下記の式（２）により電圧デューティ比ＤＵＴを算出する。
ＤＵＴ（ｎ）＝Ｐ（ｎ）＋Ｉ（ｎ）＋Ｄ（ｎ）      （２）
【００７８】
上記式（２）において、ＰはＰ項（比例項）であり、ＩはＩ項（積分項）であり、ＤはＤ項（微分項）である。Ｐ、Ｉ、Ｄは、それぞれ下記の式（３）〜（５）で定義される。
Ｐ（ｎ）＝Ｋｐ｛Ｉｔａｒ（ｎ＋１）−Ｉｄ（ｎ）｝      （３）
Ｉ（ｎ）＝Ｋｉ｛Ｉｔａｒ（ｎ）−Ｉｄ（ｎ）｝      （４）
Ｄ（ｎ）＝Ｄ（ｎ−１）＋Ｋｄ｛Ｉｔａｒ（ｎ）−Ｉｄ（ｎ）｝      （５）
【００７９】
上記式（２）〜（５）において、「ｎ」は今回のスイッチング周期Ｐｓｗにおける値を示し、「ｎ−１」は前回のスイッチング周期Ｐｓｗにおける値を示し、「ｎ＋１」は次回のスイッチング周期Ｐｓｗにおける値を示す。ＫｐはＰ項ゲインであり、ＫｉはＩ項ゲインであり、ＫｄはＤ項ゲインである。
【００８０】
駆動回路３６４は、電圧デューティ比算出部３６２からの電圧デューティ比ＤＵＴを用いて、アクチュエータ２１０に駆動電圧Ｖｄを印加する。ここでは、パルス幅変調（ＰＷＭ：pulse width modulation）を用いるため、駆動電圧Ｖｄは固定電圧である。なお、駆動回路３６４は、図１２に示すようにアクチュエータ駆動制御部２５８に含まれてもよいし、エンジンマウント２０２ｆ、２０２ｒの一部として構成してもよい。
【００８１】
図１４は、図１２の駆動回路３６４の一部を示す回路図である。この駆動回路３６４は、３つのスイッチング素子３８０ａ、３８０ｂ、３８０ｃと、４つのダイオード３８２ａ、３８２ｂ、３８２ｃ、３８２ｄと、を含んで構成される。
【００８２】
Ｖｉｎはバッテリ１６（図１）からの入力電圧を示すと共に、Ｖｏｕｔはバッテリ１６側への出力電圧を示している。入力電圧Ｖｉｎは、バッテリ１６の出力電圧（バッテリ電圧Ｖｂａｔ）そのものではなく、図示しない昇圧回路により昇圧された電圧であってもよい。
【００８３】
スイッチング素子３８０ａ及びダイオード３８２ａは、バッテリ１６とコイル２１８の間で並列配置される。スイッチング素子３８０ｂ及びダイオード３８２ｂは、コイル２１８とグラウンドの間で並列配置される。スイッチング素子３８０ｃ及びダイオード３８２ｄは、コイル２１８とグラウンドの間に並列配置される。ダイオード３８２ｃは、バッテリ１６とコイル２１８の間に配置される。
【００８４】
この駆動回路３６４は、スイッチング素子３８０ａ〜３８０ｃのオン・オフ動作に伴って、３つの通電経路（以下、ループ１〜３）を形成する。破線で示すループ１は、スイッチング素子３８０ａ、３８０ｃがオンとなり、スイッチング素子３８０ｂがオフとなった際のループである。実線で示すループ２は、スイッチング素子３８０ｂ、３８０ｃがオンとなり、スイッチング素子３８０ａがオフとなった際のループである。一点鎖線で示すループ３は、スイッチング素子３８０ａ〜３８０ｃがいずれもオフとなった際のループである。
【００８５】
＜各々の演算項Ｋｐ、Ｋｉ、Ｋｄの切替え＞
図１５は、図１〜図３のアクチュエータ２１０の目標電流波形Ｗｉ及び駆動電流Ｉｄの波形の例を示す図である。目標電流波形Ｗｉは、アクチュエータ２１０が１回進退する動作周期（つまり、波形演算周期Ｐｗｃ）を複数のポイントに分割して設定される。目標電流波形Ｗｉ（又は目標電流配列Ｍｉ）は、複数の目標電流Ｉｔａｒにより形成される。
【００８６】
目標電流Ｉｔａｒの区切りは、スイッチング素子３８０ａ〜３８０ｃのスイッチング周期Ｐｓｗ（単位：ｓｅｃ）に対応する。エンジンマウント２０２ｆ、２０２ｒに対して、同一の又は異なる振幅で位相をずらすように目標電流Ｉｔａｒが設定される。
【００８７】
なお、上記した通り、アクチュエータ駆動制御部２５８は、駆動回路３６４によるループ１〜３を用いてアクチュエータ２１０を制御する。電流上がり区間Ｚｕｐ（図１１Ｃ）にて駆動電流Ｉｄを増加させる際は、ループ１とループ２を交互に入れ替える。一方、電流下り区間Ｚｄｗｎ（図１１Ｃ）にて駆動電流Ｉｄを減少させる際は、ループ３とループ２を交互に入れ替える。
【００８８】
つまり、ＡＣＭ−ＥＣＵ２０６は、ループ１〜３を以下のように切り替えると共に、ループ１〜３の種類に応じて演算項Ｋｐ、Ｋｉ、Ｋｄを切り替える。これにより、図１５に示す電流上がり区間Ｚｕｐ及び電流下り区間Ｚｄｗｎを実現することができる。
Ｉｔａｒ−Ｉｄ≧０の場合              →ループ１
Ｉｔａｒ−Ｉｄ＜０且つＤＵＴ＞０の場合→ループ２
Ｉｔａｒ−Ｉｄ＜０且つＤＵＴ≦０の場合→ループ３
【００８９】
＜内部温度Ｔａｃｍの推定方法＞
続いて、本実施形態における内部温度Ｔａｃｍの推定方法について説明する。ここでは、立ち下がり時の駆動電流Ｉｄを用いて内部温度Ｔａｃｍを推定する。
【００９０】
（内部温度Ｔａｃｍの推定モデル）
図１６は、低温時及び高温時における目標電流波形Ｗｉと駆動電流Ｉｄとの関係の一例を示す図である。本図では、目標電流波形Ｗｉ及び駆動電流Ｉｄが減少している際（換言すると、電流下り区間Ｚｄｗｎ）の様子が示されている。
【００９１】
低温時の駆動電流Ｉｄと高温時の駆動電流Ｉｄとでは、傾き（すなわち、駆動電流Ｉｄの時間微分値又は電流傾き偏差ΔＩ（単位：Ａ／ｓｅｃ））が異なる。特に、ループ２により駆動電流Ｉｄが減少する際の傾きが大きく異なる。そこで、内部温度推定部２６０は、ループ２における駆動電流Ｉｄの傾きの変化（ここでは、電圧デューティ比ＤＵＴの変化度合い）を利用して内部温度Ｔａｃｍを推定する。
【００９２】
図１７は、本実施形態における温度推定の原理を説明するための図である。内部温度推定部２６０は、平均デューティ比ＤＵＴａｖｅ及びポイント数Ｎｐを用いて内部温度Ｔａｃｍを推定する。ここで、平均デューティ比ＤＵＴａｖｅは、所定区間における電圧デューティ比ＤＵＴの平均である。
【００９３】
ここで、所定区間は、電流下り区間Ｚｄｗｎの全部又は一部に設定される。電流下り区間Ｚｄｗｎの一部とは、例えば、電圧デューティ比ＤＵＴが第１デューティ比閾値ＴＨｄｕｔ１から第２デューティ比閾値ＴＨｄｕｔ２（＜ＴＨｄｕｔ１）に下がるまでの区間に相当する。
【００９４】
内部温度Ｔａｃｍが基準温度Ｔｒｅｆ以上である場合、電流下り区間Ｚｄｗｎの平均デューティ比ＤＵＴａｖｅは、内部温度Ｔａｃｍに対して概ね線形的に変化する傾向がある。そこで、ＡＣＭ−ＥＣＵ２０６は、この傾向を考慮して、この平均デューティ比ＤＵＴａｖｅから内部温度Ｔａｃｍを推定する。
【００９５】
具体的には、ポイント数Ｎｐ（ここでは、振動周期Ｐｖ１〜Ｐｖ３）に対応付けられた内部温度−平均デューティ比特性（以下、特性曲線５００）を特定する（例えば、複数の特性曲線５００ａ〜５００ｃのうち１つの特性曲線５００ｂ）と共に、電圧デューティ比ＤＵＴに基づいて平均デューティ比ＤＵＴａｖｅ（例えば、ＤＵＴａｖｅ１）を特定する。
【００９６】
そして、特性曲線５００において、平均デューティ比ＤＵＴａｖｅに対応する内部温度Ｔａｃｍ（例えば、Ｔａｃｍ１）を、現在の内部温度Ｔａｃｍとする。なお、平均デューティ比変化率Ｒｄｕｔａｖｅは、平均デューティ比ＤＵＴａｖｅが１増加したときの内部温度Ｔａｃｍの変化量を示す。
【００９７】
特性曲線５００から理解されるように、駆動電流Ｉｄの立ち下がり時（換言すると、ループ２とループ３を交互に実行する際）には、平均デューティ比ＤＵＴａｖｅが内部温度Ｔａｃｍに応じて略線形で変化する。なぜならば、内部温度Ｔａｃｍが上昇すると、アクチュエータ２１０内のエアが膨張して駆動軸２１６とコイル２１８の間のエアギャップが広がり、コイル２１８のインダクタンスが減少して応答性が変化するためである。
【００９８】
一方、内部温度Ｔａｃｍが基準温度Ｔｒｅｆ未満である場合、特性曲線５００を１次直線に近似することが困難である。なぜならば、内部温度Ｔａｃｍが低温のときには、ゴム板２１４（図２）の硬化度合いが大きくなり、内部温度Ｔａｃｍと発生力Ｆａｃｍが比例しないためである。
【００９９】
なお、駆動電流Ｉｄの立ち上がり時（換言すると、ループ１とループ３を交互に実行する際）には、上記した応答性の変化は相対的に小さい。よって、平均デューティ比ＤＵＴａｖｅは、内部温度Ｔａｃｍにかかわらず略一定である。
【０１００】
（平均デューティ比ＤＵＴａｖｅの算出方法）
図１８Ａは、目標電流波形Ｗｉの一例を示す。図１８Ａの目標電流波形Ｗｉ（目標電流配列Ｍｉ）は、図１１Ｃの目標電流波形Ｗｉにゼロ点補正を行ったものである。図１８Ｂは、図１８Ａの目標電流配列Ｍｉを用いた場合の電圧デューティ比ＤＵＴを示す図である。このため、図１８Ａ及び図１８Ｂの横軸は同じである。図１８Ｃは、図１８Ｂの一部を拡大して示す図である。
【０１０１】
アクチュエータ駆動制御部２５８は、算出対象区間Ｚｔａｒにわたる電圧デューティ比ＤＵＴの平均値（以下、平均デューティ比ＤＵＴａｖｅ）を算出する。ここで、算出対象区間Ｚｔａｒは、電圧デューティ比ＤＵＴが、ＴＨｄｕｔ１≦ＤＵＴ≦ＴＨｄｕｔ２の関係を満たす区間に相当する。
【０１０２】
ここでは、アクチュエータ駆動制御部２５８は、算出対象区間Ｚｔａｒにおける電圧デューティ比ＤＵＴから回帰直線５１０を算出した後、回帰直線５１０の最大値ＤＵＴｍａｘと最小値ＤＵＴｍｉｎの平均値（（ＤＵＴｍａｘ＋ＤＵＴｍｉｎ）／２）を、出力平均デューティ比ＤＵＴａｖｅｏｕｔとして算出する。
【０１０３】
＜内部温度推定部２６０の詳細＞
図１９は、図３の内部温度推定部２６０の機能を示す詳細なブロック図である。この内部温度推定部２６０は、出力波形特定部５３０と、出力平均デューティ比算出部５３２と、基準平均デューティ比算出部５３４と、平均デューティ比変化率算出部５３６と、内部温度算出部５３８と、を含んで構成される。
【０１０４】
出力波形特定部５３０は、アクチュエータ駆動制御部２５８の電圧デューティ比算出部３６２（図１２）から個々の電圧デューティ比ＤＵＴを取得する。そして、出力波形特定部５３０は、個々の電圧デューティ比ＤＵＴを波形演算周期Ｐｗｃ毎にまとめて、出力デューティ比配列Ｍｉｏｕｔとして出力する。
【０１０５】
出力平均デューティ比算出部５３２は、出力波形特定部５３０からの出力デューティ比配列Ｍｉｏｕｔに基づいて、出力平均デューティ比ＤＵＴａｖｅｏｕｔを算出する。図１８Ａ〜図１８Ｃを参照して説明した算出方法を用いることができる。
【０１０６】
基準平均デューティ比算出部５３４は、１次指令電流Ｉｃ１（図５）とポイント数Ｎｐ（図５及び図６）に基づいて、基準平均デューティ比ＤＵＴａｖｅｒｅｆ（図１７）を算出する。より具体的には、基準平均デューティ比算出部５３４は、図２０の基準平均デューティ比マップ５５０を参照することで、１次指令電流Ｉｃ１とポイント数Ｎｐの組合せに対応する基準平均デューティ比ＤＵＴａｖｅｒｅｆを読み出す。
【０１０７】
図２０は、図１９の基準平均デューティ比算出部５３４が用いる基準平均デューティ比マップ５５０を示す。基準平均デューティ比マップ５５０は、１次指令電流Ｉｃ１とポイント数Ｎｐの組合せと基準平均デューティ比ＤＵＴａｖｅｒｅｆの関係を規定する。
【０１０８】
図２１は、ポイント数Ｎｐと平均デューティ比ＤＵＴａｖｅとの関係の一例を示す図である。ここでは、ポイント数Ｎｐが奇数である場合と、偶数である場合の両方を示している。ポイント数Ｎｐは、エンジン回転速度Ｎｅの逆数（１／Ｎｅ）に概ね比例する。
【０１０９】
本図に示すように、ポイント数Ｎｐが奇数である場合と偶数である場合とで、ポイント数Ｎｐと平均デューティ比ＤＵＴａｖｅとの関係が異なっている。これは、ポイント数Ｎｐが奇数の場合と偶数の場合とで、ポイントの位置がずれるためと考えられる。この場合、ポイント数Ｎｐの偶奇性に応じた基準平均デューティ比マップ５５０を準備してもよい。
【０１１０】
平均デューティ比変化率算出部５３６は、１次指令電流Ｉｃ１（図５）とポイント数Ｎｐ（図５及び図６）に基づいて平均デューティ比変化率Ｒｄｕｔａｖｅ（図１７）を算出する。より具体的には、平均デューティ比変化率算出部５３６は、図２２の平均デューティ比変化率マップ５６０を参照することで、１次指令電流Ｉｃ１とポイント数Ｎｐの組合せに対応する平均デューティ比変化率Ｒｄｕｔａｖｅを読み出す。
【０１１１】
図２２は、図１９の平均デューティ比変化率算出部５３６が用いる平均デューティ比変化率マップ５６０を示す。平均デューティ比変化率マップ５６０は、１次指令電流Ｉｃ１とポイント数Ｎｐの組合せと平均デューティ比変化率Ｒｄｕｔａｖｅの関係を規定する。なお、基準平均デューティ比マップ５５０の場合と同様に、ポイント数Ｎｐの偶奇性に応じた平均デューティ比変化率マップ５６０を準備してもよい。
【０１１２】
内部温度算出部５３８は、出力平均デューティ比ＤＵＴａｖｅｏｕｔと、基準平均デューティ比ＤＵＴａｖｅｒｅｆと、平均デューティ比変化率Ｒｄｕｔａｖｅと、基準温度Ｔｒｅｆとに基づいて、アクチュエータ２１０の内部温度Ｔａｃｍを算出する。具体的には、内部温度算出部５３８は、下記の式（６）に基づいて内部温度Ｔａｃｍを算出する。なお、この式（６）は、特性曲線５００（図１７）には、内部温度Ｔａｃｍが基準温度Ｔｒｅｆ以上において高いリニアリティがみられることを前提としている。
Ｔａｃｍ＝（ＤＵＴａｖｅｏｕｔ−ＤＵＴａｖｅｒｅｆ）・Ｒｄｕｔａｖｅ＋Ｔｒｅｆ      （６）
【０１１３】
＜内部温度調整部２６２の詳細＞
（１．従来方法の問題点）
図２３は、比較例における温度補正の実行状態の遷移を示すタイムチャートである。このタイムチャートは、上側から順に、［１］補正判断（補正モードフラグＦｌｇｍｏｄ）、［２］補正電流（補正値Ｉｃｒｔ１）、［３］出力電流（目標電流波形Ｗｉ）、及び［４］車体１４の振動量の時系列を示す。
【０１１４】
本図では、時間ｔ＝Ｔｓｔｒにて、補正の実行状態が「ＯＮ」から「ＯＦＦ」に切り替えられた場合を想定する。つまり、補正の実行状態は、ｔ＜Ｔｓｔｒの時間範囲では「ＯＮ」に相当し、ｔ≧Ｔｓｔｒの時間範囲では「ＯＦＦ」に相当する。
【０１１５】
ｔ＜Ｔｓｔｒの時間範囲では、内部温度Ｔａｃｍの変化に応じて、補正値Ｉｃｒｔ１及び目標電流波形Ｗｉがそれぞれ緩やかに変化する。その結果、アクチュエータ２１０の適切な駆動制御により、エンジン１２から車体１４への振動伝達（車体１４の振動量）が抑制された状態が維持される。
【０１１６】
ところが、時間ｔ＝Ｔｓｔｒにおいて、補正「ＯＦＦ」に伴って補正値Ｉｃｒｔ１の算出を停止（Ｉｃｒｔ１＝０に設定）することで、補正値Ｉｃｒｔ１及び目標電流波形Ｗｉが急激に低下する。その結果、アクチュエータ２１０の不連続な駆動制御（つまり、過補正）により、車体１４の振動量が一時的に強まる。
【０１１７】
このように、例えば、補正の実行状態が切り替わるタイミングにおいて、過補正に起因する能動型防振装置２００の一時的な性能低下が起こり易くなる。そこで、状況に応じて発生する補正の不具合を予防すべく、指令値補正部２５４における補正値Ｉｃｒｔ１を調整する手法を提案する。
【０１１８】
（２．動作の概略）
内部温度調整部２６２は、補正値Ｉｃｒｔ１の計算に供される計算用温度Ｔｃａｌを算出することで、補正値Ｉｃｒｔ１を調整する補正値調整手段として機能する。この場合、指令値補正部２５４は、内部温度推定部２６０からの内部温度Ｔａｃｍに代えて、内部温度調整部２６２の調整により得た計算用温度Ｔｃａｌに応じて動作指令値Ｃを補正する。これにより、指令値補正部２５４が実行する演算処理を変更せずに計算用温度Ｔｃａｌのみを変更することで、間接的に補正値Ｉｃｒｔ１を調整することができる。
【０１１９】
また、内部温度調整部２６２は、指令値補正部２５４により補正の実行状態が切り替わるタイミング（ＯＮ→ＯＦＦ又はＯＦＦ→ＯＮ）で補正値Ｉｃｒｔ１を調整する。一般的に言えば、補正の実行状態が切り替わるタイミングにて補正の不具合が発生する可能性が高いので、不具合の予防効果が顕著に現われる。
【０１２０】
（３．具体的動作）
続いて、内部温度調整部２６２の具体的動作（内部温度Ｔａｃｍの調整方法）について、図２４のフローチャートを参照しながら詳細に説明する。なお、このフローチャートは、波形演算周期Ｐｗｃおきに定期的に実行されるとする。
【０１２１】
図２４に示すように、内部温度調整部２６２は、内部温度Ｔａｃｍ及び補正モードフラグＦｌｇｍｏｄを取得した後（ステップＳ１）、補正モードフラグＦｌｇｍｏｄを参照することで実行状態の変化を判別する（ステップＳ２）。
【０１２２】
少なくとも今回の温度補正が「ＯＮ」である場合（ステップＳ２：ＯＮ→ＯＮ、ＯＦＦ→ＯＮ）、内部温度調整部２６２は、ステップＳ１で取得した内部温度Ｔａｃｍを計算用温度Ｔｃａｌに代入し、この計算用温度Ｔｃａｌを指令値補正部２５４に向けて出力する（ステップＳ３）。
【０１２３】
ステップＳ２に戻って、前回及び今回の温度補正が両方とも「ＯＦＦ」である場合（ステップＳ２：ＯＦＦ→ＯＦＦ）、内部温度調整部２６２は、後述する調整パラメータが算出されているか否かを判定する（ステップＳ４）。
【０１２４】
例えば、温度補正が継続して「ＯＦＦ」であり、調整パラメータが算出されていない場合（ステップＳ４：ＮＯ）、内部温度調整部２６２は、基準温度Ｔｒｅｆを計算用温度Ｔｃａｌに代入し、この計算用温度Ｔｃａｌを指令値補正部２５４に向けて出力する（ステップＳ５）。なお、基準温度Ｔｒｅｆは、補正値Ｉｃｒｔ１＝０（補正ＯＦＦ）と算出される温度（以下、目標温度）に相当する。
【０１２５】
ステップＳ２に戻って、前回の温度補正が「ＯＮ」であり、かつ今回の温度補正が「ＯＦＦ」である場合（ステップＳ２：ＯＮ→ＯＦＦ）、内部温度調整部２６２は、次のステップＳ６、Ｓ７に進み、内部温度Ｔａｃｍの調整パラメータ（ここでは、温度偏差ΔＴ、温度分割数Ｎｄｉｖ）を算出する。
【０１２６】
内部温度調整部２６２は、下記の式（７）に基づいて温度偏差ΔＴを、下記の式（８）に基づいて温度分割数Ｎｄｉｖをそれぞれ算出する（ステップＳ６、Ｓ７）。ここで、温度調整量ΔＴｏは、波形演算周期Ｐｗｃ当たりの計算用温度Ｔｃａｌの変化量に相当し、例えば、ΔＴｏ＝１０［℃］である。また、ＩＮＴ（・）は、括弧内の変数における整数部分を示す関数である。
ΔＴ＝Ｔａｃｍ−Ｔｒｅｆ      （７）
Ｎｄｉｖ＝１＋ＩＮＴ（ΔＴ／ΔＴｏ）      （８）
【０１２７】
その後、内部温度調整部２６２は、ステップＳ１で取得した内部温度Ｔａｃｍを計算用温度Ｔｃａｌに代入し、この計算用温度Ｔｃａｌを指令値補正部２５４に向けて出力する（ステップＳ３）。このフローチャートを次に実行する際、前回及び今回の温度補正が両方とも「ＯＦＦ」であり（ステップＳ２：ＯＦＦ→ＯＦＦ）、かつ、調整パラメータが算出されているので（ステップＳ４：ＹＥＳ）、ステップＳ８に進む。
【０１２８】
内部温度調整部２６２は、既に算出した調整パラメータを用いて、内部温度Ｔａｃｍを段階的に調整する（ステップＳ８）。具体的には、内部温度調整部２６２は、下記の式（９）に基づいて計算用温度Ｔｃａｌを更新する。
Ｔｃａｌ（今回）＝Ｔｃａｌ（前回）−ΔＴ／Ｎｄｉｖ      （９）
【０１２９】
内部温度調整部２６２は、内部温度Ｔａｃｍの調整が終了したか否かを判定する（ステップＳ９）。具体的には、内部温度調整部２６２は、（９）式に基づく計算用温度Ｔｃａｌの更新回数が温度分割数Ｎｄｉｖを下回る場合には、調整がまだ終了していないと判定する（ステップＳ９：ＮＯ）。
【０１３０】
一方、内部温度調整部２６２は、上記した更新回数が温度分割数Ｎｄｉｖに等しい場合に調整が終了したと判定し（ステップＳ９：ＹＥＳ）、調整パラメータを初期化する（ステップＳ１０）。
【０１３１】
図２５は、図２４のフローチャートに基づく内部温度Ｔａｃｍの調整結果を模式的に示す図である。グラフの横軸は調整時間（単位：ｓｅｃ）を示すと共に、グラフの縦軸は計算用温度Ｔｃａｌ（単位：℃）を示す。なお、調整時間は、調整の開始時点（図２３及び図２６のｔ＝Ｔｓｔｒに相当）からの経過時間に相当する。
【０１３２】
本図では、調整の開始時点における内部温度ＴａｃｍがＴａｃｍ１〜Ｔａｃｍ３（なお、Ｔｒｅｆ＜Ｔａｃｍ１＜Ｔａｃｍ２＜Ｔａｃｍ３）である、３種類の調整パターンを表記している。本図から理解されるように、内部温度調整部２６２は、内部温度Ｔａｃｍ＝Ｔａｃｍ１〜Ｔａｃｍ３を、一定の変化量（ΔＴｏ／Ｐｗｃ）で減少させながら複数の波形演算周期Ｐｗｃにわたって、基準温度Ｔｒｅｆに近づける調整を行う。
【０１３３】
このように、内部温度調整部２６２は、補正の実行状態が切り替わるタイミングで推定された内部温度Ｔａｃｍを、段階的に基準温度Ｔｒｅｆ（目標温度）に近づける調整を行ってもよい。これにより、補正の実行状態が切り替わるタイミングで発生する制御の不連続性を抑制可能となり、例えば過補正に起因する能動型防振装置２００の性能低下を防止することができる。
【０１３４】
また、内部温度調整部２６２は、補正の実行状態が切り替わるタイミングで推定された内部温度Ｔａｃｍを、時間の経過につれて増加又は減少させながら基準温度Ｔｒｅｆに近づける調整を行ってもよい。これにより、円滑かつ緩やかに補正値Ｉｃｒｔ１を調整することができる。この調整は、単調な増加又は減少であってもよいし、線形的な増加又は減少であってもよい。
【０１３５】
また、内部温度推定部２６０が演算周期（例えば、波形演算周期Ｐｗｃ）おきに内部温度Ｔａｃｍを推定する場合、内部温度調整部２６２は、補正の実行状態が切り替わるタイミングで推定された内部温度Ｔａｃｍを、内部温度Ｔａｃｍにかかわらず一定の変化量（温度調整量ΔＴｏ）だけ、１つ又は複数（ここでは、温度分割数Ｎｄｉｖ）の演算周期にわたって増加又は減少させながら基準温度Ｔｒｅｆに近づける調整を行ってもよい。これにより、内部温度Ｔａｃｍと基準温度Ｔｒｅｆの乖離にかかわらず、略同じ変化度合いで補正値Ｉｃｒｔ１を調整することができる。
【０１３６】
また、内部温度調整部２６２は、補正のオン状態からオフ状態に切り替わるタイミングで推定された内部温度Ｔａｃｍを、補正値Ｉｃｒｔ１＝０と算出される基準温度Ｔｒｅｆに近づける調整を行ってもよい。これにより、内部温度Ｔａｃｍの調整を通じて、補正をオン状態からオフ状態に円滑に移行させることができる。
【０１３７】
（４．補正値Ｉｃｒｔ１の調整による効果）
図２６は、実施例における温度補正の実行状態の遷移を示すタイムチャートである。本図は、図２３の場合と同様に定義されたタイムチャートを示し、時間ｔ＝Ｔｓｔｒにて補正の実行状態が「ＯＮ」から「ＯＦＦ」に切り替えられた場合を想定する。
【０１３８】
ｔ＜Ｔｓｔｒの時間範囲では、内部温度Ｔａｃｍの変化に応じて、補正値Ｉｃｒｔ１及び目標電流波形Ｗｉがそれぞれ緩やかに変化する。その結果、アクチュエータ２１０の適切な駆動制御により、エンジン１２から車体１４への振動伝達（車体１４の振動量）が抑制された状態が維持される。
【０１３９】
そして、時間ｔ＝Ｔｓｔｒにおいて、温度補正が「ＯＦＦ」に移行するに伴って、推定された内部温度Ｔａｃｍを段階的に基準温度Ｔｒｅｆに近づける調整を開始する。そうすると、Ｔｓｔｒ≦ｔ≦Ｔｅｎｄ（調整の終了時間）の時間範囲では、計算用温度Ｔｃａｌの変化に応じて、補正値Ｉｃｒｔ１及び目標電流波形Ｗｉがそれぞれ緩やかに変化する。その結果、アクチュエータ２１０の適切な駆動制御により、エンジン１２から車体１４への振動伝達（車体１４の振動量）が抑制された状態が維持される。
【０１４０】
［能動型防振装置２００による効果］
以上のように、この能動型防振装置２００は、エンジン１２（駆動源）と車体１４との間に配置されたアクチュエータ２１０と、アクチュエータ２１０に能動的振動を生成させることでエンジン１２から車体１４への振動伝達を抑制するＡＣＭ−ＥＣＵ２０６（コンピュータ）と、を備える装置である。
【０１４１】
ＡＣＭ−ＥＣＵ２０６は、［１］エンジン１２の振動情報（振動振幅Ａｖ、振動周期Ｐｖ）からアクチュエータ２１０の動作指令値Ｃを算出する指令値算出部２５２と、［２］アクチュエータ２１０の内部温度Ｔａｃｍに応じた補正値Ｉｃｒｔ１を算出し、補正値Ｉｃｒｔ１を用いて動作指令値Ｃを補正する指令値補正部２５４と、［３］補正した動作指令値Ｃに基づく電圧デューティ比ＤＵＴを用いてアクチュエータ２１０に駆動電圧Ｖｄを印加するアクチュエータ駆動制御部２５８と、［４］所定区間の電圧デューティ比ＤＵＴの平均としての平均デューティ比ＤＵＴａｖｅに基づいて、内部温度Ｔａｃｍを推定する内部温度推定部２６０と、［５］指令値補正部２５４における補正値Ｉｃｒｔ１を調整する内部温度調整部２６２（補正値調整手段）と、を備える。
【０１４２】
アクチュエータ２１０の内部温度Ｔａｃｍに対する平均デューティ比ＤＵＴａｖｅの関係を示す特性曲線５００には、高い相関性及び高い再現性が部分的にみられる。そこで、この傾向を踏まえて動作指令値Ｃを補正することで、アクチュエータ２１０の内部温度Ｔａｃｍを高精度に推定することができる。しかも、動作指令値Ｃの補正に供される補正値Ｉｃｒｔ１を調整することで、状況に応じて発生する補正の不具合を解消することができる。
【０１４３】
［補足］
なお、この発明は、上述した実施形態に限定されるものではなく、この発明の主旨を逸脱しない範囲で自由に変更できることは勿論である。或いは、技術的に矛盾が生じない範囲で各々の構成を任意に組み合わせてもよい。
【０１４４】
上記した実施形態では、エンジン車両である車両１０に能動型防振装置２００を搭載しているが、適用範囲はこれに限らない。例えば、ハイブリッド車両を含む車両全般であってもよいし、駆動源を備える移動体（船舶や航空機）、製造装置、ロボット又は家電製品に適用してもよい。
【０１４５】
上記した実施形態では、内部温度推定部２６０は、式（６）に基づいて内部温度Ｔａｃｍを推定しているが、推定モデルはこれに限らない。例えば、高い相関性及び高い再現性があることを条件に、特性曲線５００を２次以上の多項式関数や指数関数を含む非線形関数で近似した推定モデルを用いてもよい。
【０１４６】
上記した実施形態では、内部温度調整部２６２は、図２４のフローチャートに従って内部温度Ｔａｃｍを調整しているが、調整方法はこれに限られない。例えば、補正の実行状態が「ＯＦＦ」から「ＯＮ」に切り替わった場合であっても、上記と同様に内部温度Ｔａｃｍを調整してもよい。
【０１４７】
上記した実施形態では、補正値調整手段としての内部温度調整部２６２は、指令値補正部２５４とは別機能として設けられているが、機能構成はこれに限らない。例えば、補正値調整手段は、指令値補正部２５４に組み込まれており、直接的に補正値Ｉｃｒｔ１を調整可能に構成されてもよい。
【符号の説明】
【０１４８】
１０…車両                              １２…エンジン（駆動源）
１４…車体                              ２００…能動型防振装置
２０２ｆ、２０２ｒ…エンジンマウント    ２０６…ＡＣＭ−ＥＣＵ（コンピュータ）
２１０…アクチュエータ                  ２５０…エンジン状態変数算出部
２５２…指令値算出部                    ２５４…指令値補正部
２５６…目標電流波形算出部              ２５８…アクチュエータ駆動制御部
２６０…内部温度推定部                  ２６２…内部温度調整部（補正値調整手段）
２６４…補正実行判断部                  ３００…ポイント数マップ
３５０…目標電流波形マップ              ５００…特性曲線
５５０…基準平均デューティ比マップ      ５６０…平均デューティ比変化率マップ
Ａｖ…振動振幅（振動情報）              Ｃ…動作指令値
ＤＵＴ…電圧デューティ比                ＤＵＴａｖｅ…平均デューティ比
Ｆｌｇｍｏｄ…補正モードフラグ          Ｉｃ１…１次指令電流
Ｉｃ２…２次指令電流                    Ｉｃｒｔ１…補正値
Ｐｃ１…１次位相                        Ｐｃ２…２次位相
Ｐｓｗ…スイッチング周期                Ｐｖ…振動周期（振動情報）
Ｐｗｃ…波形演算周期（演算周期）        Ｔａｃｍ…内部温度
Ｖｄ…駆動電圧                          Ｗｉ…目標電流波形
Ｚｔａｒ…算出対象区間（所定区間）
【図１】
2019100382.tif 000003
【図２】
2019100382.tif 000004
【図３】
2019100382.tif 000005
【図４】
2019100382.tif 000006
【図５】
2019100382.tif 000007
【図６】
2019100382.tif 000008
【図７】
2019100382.tif 000009
【図８】
2019100382.tif 000010
【図９】
2019100382.tif 000011
【図１０】
2019100382.tif 000012
【図１１】
2019100382.tif 000013
【図１２】
2019100382.tif 000014
【図１３】
2019100382.tif 000015
【図１４】
2019100382.tif 000016
【図１５】
2019100382.tif 000017
【図１６】
2019100382.tif 000018
【図１７】
2019100382.tif 000019
【図１８】
2019100382.tif 000020
【図１９】
2019100382.tif 000021
【図２０】
2019100382.tif 000022
【図２１】
2019100382.tif 000023
【図２２】
2019100382.tif 000024
【図２３】
2019100382.tif 000025
【図２４】
2019100382.tif 000026
【図２５】
2019100382.tif 000027
【図２６】
2019100382.tif 000028