(57)【要約】【課題】劣化したときの蓄電デバイスの特性を精度よく推定することができるシミュレーション方法及びシミュレーション装置を提供する。【解決手段】このシミュレーション方法は、蓄電デバイスの等価回路モデルを用いてシミュレーションを行う方法であって、等価回路モデルを流れる電流Ｉ（ｔ）に基づいて等価回路モデルの端子電圧Ｖ（ｔ）を計算するステップを含む。等価回路モデルは、複数の特性パラメータを含む。少なくとも一つの特性パラメータは、蓄電デバイスの劣化の影響を表す時間関数を含む。時間関数は、蓄電デバイスの暗電流放電時の充電率及び休止時の充電率の一方または双方の時間積分と、時間積分に乗算された劣化速度を表す係数とを含む項を有する。【選択図】図５
2018009940.tif 000002
【特許請求の範囲】
【請求項１】
蓄電デバイスの等価回路モデルを用いてシミュレーションを行う方法であって、
前記等価回路モデルを流れる電流に基づいて前記等価回路モデルの端子電圧を計算するステップを含み、
前記等価回路モデルが複数の特性パラメータを含んでおり、
少なくとも一つの前記特性パラメータが、前記蓄電デバイスの劣化の影響を表す時間関数を含んでおり、
前記時間関数が、前記蓄電デバイスの暗電流放電時の充電率及び休止時の充電率の一方または双方の時間積分と、前記時間積分に乗算された劣化速度を表す係数とを含む項を有する、シミュレーション方法。
【請求項２】
前記係数が前記蓄電デバイスの温度に応じて変化する、請求項１に記載のシミュレーション方法。
【請求項３】
前記蓄電デバイスが鉛蓄電池である、請求項１または２に記載のシミュレーション方法。
【請求項４】
蓄電デバイスの等価回路モデルを用いてシミュレーションを行う装置であって、
前記等価回路モデルを流れる電流に基づいて前記等価回路モデルの端子電圧を計算する電圧計算部を含み、
前記等価回路モデルが複数の特性パラメータを含んでおり、
少なくとも一つの前記特性パラメータが、前記蓄電デバイスの劣化の影響を表す時間関数を含んでおり、
前記時間関数が、前記蓄電デバイスの暗電流放電時の充電率及び休止時の充電率の一方または双方の時間積分と、前記時間積分に乗算された劣化速度を表す係数とを含む項を有する、シミュレーション装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、シミュレーション方法及びシミュレーション装置に関する。
【背景技術】
【０００２】
特許文献１には、電池モデル同定方法が記載されている。この方法では、電池に入力される電流波形を電流センサを用いて測定し、電池の端子電圧の電圧波形を電圧センサを用いて測定する。そして、システム同定演算部が、これらの電流波形及び電圧波形に基づいて、電池モデルのシステム同定を行う。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１６−３９６３号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
現在、車載用蓄電デバイスとして、鉛蓄電池、リチウムイオン電池、リチウムイオンキャパシタといった様々な蓄電デバイスが用いられている。そして、これらの電池を組み合わせて、様々な種類のハイブリッド方式が実用化されている。例えば、メイン蓄電デバイスとは別に、減速の際のエネルギー回生及び補機への給電のためのサブ蓄電デバイスを設ける、いわゆるμＨＥＶ（Hybrid Electric Vehicle）方式が近年特に有用とされている。
【０００５】
一方、例えば車両の燃費シュミレーションにおいては、エンジン及び蓄電デバイスといった様々な動力源並びに負荷をモデル化し、規定の走行パターンを該モデルに入力して燃費を算出することが行われている。このような燃費シュミレーション等に含まれる蓄電デバイスのシミュレーションにおいて、ハイブリッド方式における複雑化した蓄電デバイス構成を正確にモデル化することは、シミュレーションを精度よく行うために極めて重要である。しかしながら、従来のシュミレーションにおいては、蓄電デバイスの使用による劣化を考慮せずに蓄電デバイスをモデル化しているので、劣化したときの蓄電デバイスの特性を精度良く推定するためには、モデルの特性パラメータを同定する際に、実際に劣化した蓄電デバイスを用意する必要があるという問題がある。
【０００６】
本発明は、実際に劣化した蓄電デバイスを用いなくても、劣化したときの蓄電デバイスの特性を精度良く推定することができるシミュレーション方法及びシミュレーション装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
本発明の一実施形態によるシミュレーション方法は、蓄電デバイスの等価回路モデルを用いてシミュレーションを行う方法であって、等価回路モデルを流れる電流に基づいて等価回路モデルの端子電圧を計算するステップを含み、等価回路モデルが複数の特性パラメータを含んでおり、少なくとも一つの特性パラメータが、蓄電デバイスの劣化の影響を表す時間関数を含んでおり、時間関数が、蓄電デバイスの暗電流放電時のＳＯＣ及び休止時のＳＯＣの一方または双方の時間積分と、該時間積分に乗算された劣化速度を表す係数とを含む項を有する。
【０００８】
また、本発明の一実施形態によるシミュレーション装置は、蓄電デバイスの等価回路モデルを用いてシミュレーションを行う装置であって、等価回路モデルを流れる電流に基づいて等価回路モデルの端子電圧を計算する電圧計算部を含み、等価回路モデルが複数の特性パラメータを含んでおり、少なくとも一つの特性パラメータが、蓄電デバイスの劣化の影響を表す時間関数を含んでおり、時間関数が、蓄電デバイスの暗電流放電時のＳＯＣ及び休止時のＳＯＣの一方または双方の時間積分と、該時間積分に乗算された劣化速度を表す係数とを含む項を有する。
【０００９】
これらのシミュレーション方法及びシミュレーション装置では、少なくとも一つの特性パラメータが、蓄電デバイスの劣化の影響を表す時間関数を含む。このような時間関数に蓄電デバイスの使用期間（トータルサイクル時間）内の適切な時間（例えば動作状態を表す数値に関わる動作時間）を入力することにより、該使用時間経過時における蓄電デバイスの劣化度合いを特性パラメータに反映させることができる。そして、本発明者の知見によれば、暗電流放電時または休止時のＳＯＣによって蓄電デバイスの劣化の度合いは大きく変化する。上記の方法及び装置では、時間関数が、劣化速度を表す係数と、蓄電デバイスの暗電流放電時のＳＯＣ及び休止時のＳＯＣの一方または双方の時間積分とを乗算した項を含む。これにより、蓄電デバイスの暗電流放電時または休止時のＳＯＣに基づく劣化度合いを精度良く表すことができる。従って、上記の方法及び装置によれば、実際に劣化した蓄電デバイスを用いなくても、劣化したときの蓄電デバイスの入出力特性を精度良く推定することができ、劣化した蓄電デバイスを用いた燃費シミュレーションなどを精度良く行うことができる。なお、この蓄電デバイスは、例えば鉛蓄電池である。
【００１０】
上記のシミュレーション方法において、係数は蓄電デバイスの温度に応じて変化してもよい。これにより、蓄電デバイスの温度に応じて変化する劣化度合いを精度良く表すことができる。
【発明の効果】
【００１１】
本発明によるシミュレーション方法及びシミュレーション装置によれば、実際に劣化した蓄電デバイスを用いなくても、劣化したときの蓄電デバイスの特性を精度良く推定することができる。
【図面の簡単な説明】
【００１２】
【図１】図１は、燃費シミュレーションを行う装置の概略構成を示す図である。
【図２】図２は、等価回路モデルの例を示す図である。
【図３】図３は、一実施形態に係るシミュレーション装置の概略構成を示す図である。
【図４】図４は、図３のシミュレーション装置のハードウェア構成の例を示す図である。
【図５】図５は、シミュレーション装置が実行する端子電圧の計算処理を説明する図である。
【図６】図６（ａ）〜図６（ｃ）は、蓄電デバイスを流れる電流の波形の例を概念的に示すグラフである。
【図７】図７は、リチウムイオン電池を例として図６（ａ）〜図６（ｃ）の電流波形を入力したときの直流抵抗の変化を示すグラフである。
【図８】図８は、図６（ａ）〜図６（ｃ）の電流波形に対応する、通電劣化項の時間積分の計算値をプロットしたグラフである。
【図９】図９は、図６（ａ）〜図６（ｃ）の電流波形に対応する、自己発熱劣化項の時間積分の計算値をプロットしたグラフである。
【図１０】図１０は、図６（ａ）〜図６（ｃ）の電流波形に対応する、休止劣化項の時間積分の計算値をプロットしたグラフである。
【図１１】図１１（ａ）〜図１１（ｃ）は、カーブフィッティングの様子を概念的に示すグラフである。
【図１２】図１２は、各特性パラメータの係数の数値例を示す図表である。
【図１３】図１３は、車両の燃費シミュレーションにおいて、蓄電デバイスの等価回路モデルに一実施形態の特性パラメータを使用した場合と、劣化による影響を考慮しない特性パラメータを使用した場合とにおける燃費誤差の最大値を比較した結果を示すグラフである。
【図１４】図１４は、各特性パラメータの同定に用いた電流波形を示す図である。
【図１５】図１５は、一実施形態による燃費シミュレーション結果の一例を示すグラフである。
【発明を実施するための形態】
【００１３】
以下、添付図面を参照しながら本発明によるシミュレーション方法及びシミュレーション装置の実施の形態を詳細に説明する。なお、図面の説明において同一の要素には同一の符号を付し、重複する説明を省略する。以下の説明において、蓄電デバイスを流れる電流とは、蓄電デバイスに入力される充電電流および蓄電デバイスから出力される放電電流の双方を指し、電流の符号が正である場合は充電を表し、電流の符号が負である場合は放電を表す。
【００１４】
一実施形態に係るシミュレーション方法及びシミュレーション装置は、たとえば、蓄電デバイスが搭載された車両の燃費シミュレーションにおいて用いられる。この蓄電デバイスは、例えばμＨＥＶ方式を採用した車両に搭載されるメイン蓄電デバイス、若しくはメイン蓄電デバイスとは別に設けられたサブ蓄電デバイスである。蓄電デバイスがサブ蓄電デバイスである場合、蓄電デバイスは、車両に搭載された１２Ｖ系の補機の消費電流を賄うために用いられる。
【００１５】
［燃費シミュレーション装置の概要］
図１は、燃費シミュレーションを行う装置の概略構成を示す図である。図１に示されるように、燃費シミュレーション装置９０は、その機能ブロックとして、入力部９１と、制御部９２と、出力部９３とを含む。
【００１６】
入力部９１は、燃費シミュレーションに必要なデータを入力する。入力データの例は、車両の走行パターンである。それ以外にも、車両に搭載されるエンジンなどの各種デバイスの特性を定めるパラメータ、蓄電デバイスの充放電の制御方法の種類、蓄電デバイスの構成、車両に搭載される補機の消費電力、および車両の重量などのデータが入力され得る。
【００１７】
制御部９２は、入力部９１によって入力されたデータを用いて、燃費シミュレーションを行う。燃費シミュレーションの具体的な手法は特に限定されないが、たとえば、次のような手順で行われる。
【００１８】
まず、制御部９２は、入力部９１によって入力された走行パターンなどから、たとえば区間ごとに、車両が走行するために要求されるパワー（以下、単に「要求パワー」という）および補機の消費電流を算出する。区間としては、停止区間、加速区間、定速走行区間、および減速区間などがある。要求パワーは、加速区間では比較的大きく、定速走行区間では比較的小さい。要求パワーは、停止区間および減速区間では０であってもよい。補機の消費電流は、補機の種類によって異なる。たとえばオーディオ機器など連続的に使用される補機の消費電流の大きさは、区間によらずほぼ一定である。これに対し、エンジンの点火装置など一時的に使用される補機の消費電流の大きさは、使用時のみ大きくなる。
【００１９】
次に、制御部９２は、区間ごとのエンジンの出力を算出する。エンジンの出力は、たとえば、停止区間ではエンジンが停止して０となり、それ以外の区間では所定の出力とされる。エンジンの出力のうち、要求パワーを上回る分の出力が、オルタネータによって電力に変換され、オルタネータから補機および蓄電デバイスに向かって供給される。オルタネータから供給される電力が補機の消費電力を上回ると、オルタネータから蓄電デバイスに電流が流れ、蓄電デバイスが充電される。オルタネータから供給される電力が補機の消費電力を下回ると、蓄電デバイスから補機に電流が流れ、蓄電デバイスが放電する。ここで、蓄電デバイスの端子電圧は、蓄電デバイスのＳＯＣおよび充放電電流の大きさなどに依存する。この蓄電デバイスの端子電圧が、たとえば、蓄電デバイスの端子電圧を計算するための等価回路モデルを用いて推定される。端子電圧の推定の詳細については後述する。蓄電デバイスの充放電電流および蓄電デバイスの端子電圧から、制御部９２は、区間ごとの蓄電デバイスの充放電電力も算出する。
【００２０】
その後、制御部９２は、全区間におけるエンジンの出力および蓄電デバイスの充放電電力の積算値を算出する。全区間におけるエンジンの出力の積算値は、入力部９１によって入力された走行パターンで車両が走行した場合に、エンジンが消費するであろうエネルギー量を示す。全区間における蓄電デバイスの充放電電力の積算値は、入力部９１によって入力された走行パターンで車両が走行した場合に、蓄電デバイスにおいて増減するであろうエネルギー量の大きさを示す。エンジンが消費するであろうエネルギー量と、蓄電デバイスにおいて減少するであろうエネルギー量との合計のエネルギー量は、入力部９１によって入力された走行パターンの車両の走行に要するエネルギー量となる。走行パターンから車両の走行距離も分かるので、当該走行距離とそれに要するエネルギー量とに基づいて、制御部９２は、所定エネルギー量当たりに走行可能な距離を燃費として算出する。
【００２１】
出力部９３は、制御部９２によって算出された燃費を出力する。これにより、入力部９１によって入力された走行パターンなどに基づく燃費シミュレーションの結果が得られる。
【００２２】
上述のように、燃費シミュレーションにおいては、蓄電デバイスの端子電圧が推定される。蓄電デバイスの端子電圧の推定精度を向上させることによって燃費シミュレーションの精度も向上するので、たとえば燃費の計算精度を向上させることを目的として、実施形態に係るシミュレーション装置（蓄電デバイスシミュレータ）が用いられてもよい。なお、以下の説明において、蓄電デバイスとしては、単一の鉛蓄電池が用いられる。蓄電デバイスは、鉛蓄電池に限られず、他の蓄電デバイスであってもよく、複数の蓄電デバイスを組み合わせた複合型の蓄電デバイスであってもよい。
【００２３】
本実施形態では、シミュレーション装置は、蓄電デバイスの端子電圧を計算するための蓄電デバイスモデルを用いて、蓄電デバイスの端子電圧を推定する。本実施形態では、蓄電デバイスモデルとして、蓄電デバイスの等価回路モデルを用いることとする。まず、等価回路モデルの例について、図２を参照して説明する。
【００２４】
［蓄電デバイスの等価回路モデル］
図２に示される例では、等価回路モデル４０は、互いに逆極性のノードＮ１およびノードＮ２の間に直列に接続された、回路１０と、回路２０と、定電圧源３０とを含む。
【００２５】
ノードＮ１およびノードＮ２は、蓄電デバイスの外部の要素と電気的に接続される部分であり、等価回路モデル４０に発生する電圧を与える。等価回路モデル４０に発生する電圧は、蓄電デバイスの端子電圧Ｖ（ｔ）である。ノードＮ１はアノードであり、蓄電デバイスを流れる電流Ｉ（ｔ）を与える。なお、電圧および電流などの時間変化する物理量を示す符号に（ｔ）などを付す場合があるが、このように示された物理量は、時刻ｔにおける当該物理量の値を意味するものとする。また、時刻ｔは、０以上の整数であり、端子電圧Ｖ（ｔ）の推定の開始時刻からの経過時間を示す。時刻ｔ＝０は、端子電圧Ｖ（ｔ）の推定の開始時刻である。
【００２６】
回路１０は、蓄電デバイスの直流インピーダンス（直流抵抗成分）を模擬する直流抵抗部である。回路１０は、抵抗器１１を含む。抵抗器１１は、蓄電デバイスの線形直流抵抗成分を模擬している。線形直流抵抗成分としては、電極の抵抗が挙げられる。抵抗器１１の抵抗値は定数である。回路１０の抵抗器１１の抵抗値によって、回路１０のインピーダンスが定まる。回路１０のインピーダンスが定まれば、等価回路モデル４０に電流Ｉ（ｔ）が流れたときに、その電流Ｉ（ｔ）が回路１０にも流れるので、電流Ｉ（ｔ）と回路１０のインピーダンスとから、回路１０に発生する電圧が計算できる。回路１０に発生する電圧を、直流抵抗電圧Ｖｄｃ（ｔ）と称し図示する。
【００２７】
回路２０は、蓄電デバイスの分極インピーダンス成分を模擬する分極モデル部である。回路２０は、並列接続された抵抗器およびコンデンサ（ＲＣ並列回路）を含む。図２に示される例では、２つのＲＣ並列回路が直列に接続されている。具体的に、並列接続された抵抗器２１およびコンデンサ２２（第１のＲＣ並列回路）と、並列接続された抵抗器２３およびコンデンサ２４（第２のＲＣ並列回路）とが、直列に接続されている。第１のＲＣ並列回路を構成する抵抗器２１の抵抗値およびコンデンサ２２の容量値は定数である。抵抗器２１は、蓄電デバイスの第１の分極抵抗成分を模擬し、コンデンサ２２は、蓄電デバイスの第１の分極容量成分を模擬している。第２のＲＣ並列回路を構成する抵抗器２３の抵抗値およびコンデンサ２４の容量値は定数である。抵抗器２３は蓄電デバイスの第２の分極抵抗成分を模擬し、コンデンサ２４は蓄電デバイスの第２の分極容量成分を模擬している。
【００２８】
なお、図２に示される例では回路２０は、第１及び第２の２つのＲＣ並列回路を含むが、回路２０は、少なくとも第１のＲＣ並列回路（抵抗器２１およびコンデンサ２２）を含んでいればよい。また、回路２０は、３つ以上のＲＣ並列回路を含んでいてもよい。
【００２９】
回路２０の各抵抗器の抵抗値および各コンデンサの容量値によって、回路２０のインピーダンスが定まる。回路２０のインピーダンスが定まれば、等価回路モデル４０に電流Ｉ（ｔ）が流れたときに、その電流Ｉ（ｔ）が回路２０にも流れるので、電流Ｉ（ｔ）と回路２０のインピーダンスとから、回路２０に発生する電圧が計算できる。回路２０に発生する電圧を、分極電圧Ｖｐｏｌ（ｔ）と称し図示する。
【００３０】
分極電圧Ｖｐｏｌは、抵抗器２１およびコンデンサ２２に発生する電圧と、抵抗器２３およびコンデンサ２４に発生する電圧との合計電圧である。抵抗器２１およびコンデンサ２２に発生する電圧を、第１分極電圧Ｖｐ１（ｔ）と称し図示する。抵抗器２３およびコンデンサ２４に発生する電圧を、第２分極電圧Ｖｐ２（ｔ）と称し図示する。すなわち、回路２０において、以下の関係式（１）が成立する。
【数１】
2018009940.tif 000003
【００３１】
ここで、抵抗器２１およびコンデンサ２２から構成される第１のＲＣ並列回路の時定数を時定数τ１とすると、時定数τ１は、抵抗器２１の抵抗値とコンデンサ２２の容量値とを乗じた値として定められる。時定数τ１は、抵抗器２１およびコンデンサ２２に発生する第１分極電圧Ｖｐ１（ｔ）の時間変化に反映される。たとえば、時定数τ１が大きいほど、第１分極電圧Ｖｐ１（ｔ）の時間変化は遅くなる。同様に、抵抗器２３およびコンデンサ２４から構成される第２のＲＣ並列回路の時定数を時定数τ２とすると、時定数τ２は、抵抗器２３の抵抗値とコンデンサ２４の容量値とを乗じた値として定められる。時定数τ２は、抵抗器２３およびコンデンサ２４に発生する第２分極電圧Ｖｐ２（ｔ）の時間変化に反映される。時定数τ１及びτ２は互いに異なる値に設定されてよい。回路２０が複数の異なる時定数を有するＲＣ並列回路を含むことで、分極電圧Ｖｐｏｌ（ｔ）の電圧の時間変化をより正確に表すことができる。各時定数は、たとえば、時定数τ１＜時定数τ２となるように設定されてよい。
【００３２】
定電圧源３０は、一定の直流（ＤＣ）電圧を有する。定電圧源３０の有する電圧は、蓄電デバイスの開放電圧（ＯＣＶ：Open Circuit Voltage）である。定電圧源３０のインピーダンスは０である。蓄電デバイスの開放電圧を、開放電圧Ｖｏｃｖ（ｔ）と称し図示する。開放電圧Ｖｏｃｖ（ｔ）は、たとえば、蓄電デバイスのＳＯＣから求められる。その場合、開放電圧Ｖｏｃｖ（ｔ）は、ＳＯＣを引数とする関数となる。蓄電デバイスの温度なども、引数に含まれてもよい。
【００３３】
以上説明した回路１０に発生する直流抵抗電圧Ｖｄｃ（ｔ）、回路２０に発生する分極電圧Ｖｐｏｌ（ｔ）および定電圧源３０が有する開放電圧Ｖｏｃｖ（ｔ）と、端子電圧Ｖ（ｔ）との間には、以下の関係式（２）が成立する。
【数２】
2018009940.tif 000004
【００３４】
以上説明した蓄電デバイスの等価回路モデル４０を用いて、実施形態に係るシミュレーション装置は、蓄電デバイスの端子電圧Ｖ（ｔ）を推定する。
【００３５】
図３は、一実施形態に係るシミュレーション装置の概略構成を示す図である。シミュレーション装置１は、その機能ブロックとして、入力部２と、ＳＯＣ計算部３と、パラメータ設定部４と、直流抵抗計算部５と、分極計算部６と、ＯＣＶ計算部７と、端子電圧計算部８とを含む。
【００３６】
図４は、図３のシミュレーション装置１のハードウェア構成の例を示す図である。図４に示されるように、シミュレーション装置１は、物理的には、一または複数のＣＰＵ（Central Processing Unit）１０１と、主記憶装置であるＲＡＭ（Random Access Memory）１０２およびＲＯＭ（Read Only Memory）１０３と、データ送受信デバイスである通信モジュール１０４と、ハードディスクおよびフラッシュメモリなどの補助記憶装置１０５と、キーボードなどのユーザの入力を受け付ける入力装置１０６と、ディスプレイなどの出力装置１０７と、を備えるコンピュータとして構成されている。図３に示されるシミュレーション装置１の各機能は、ＣＰＵ１０１およびＲＡＭ１０２などのハードウェア上に一または複数の所定のコンピュータソフトウェアを読み込ませることにより、ＣＰＵ１０１の制御のもとで通信モジュール１０４、入力装置１０６、および出力装置１０７を動作させるとともに、ＲＡＭ１０２および補助記憶装置１０５におけるデータの読み出しおよび書き込みを行うことで実現される。なお、上記の説明はシミュレーション装置１のハードウェア構成として説明したが、燃費シミュレーション装置９０がＣＰＵ１０１、ＲＡＭ１０２およびＲＯＭ１０３などの主記憶装置、通信モジュール１０４、補助記憶装置１０５、入力装置１０６、および出力装置１０７などを含む通常のコンピュータシステムとして構成されてもよい。
【００３７】
再び図３を参照して、シミュレーション装置１の各機能の詳細を説明する。入力部２は、蓄電デバイスへの指定値（ｂａｔ＿ｄｅｍａｎｄ）を入力する部分である。指定値は、たとえば上述の燃費シミュレーション装置９０による燃費計算において蓄電デバイスに要求される、充放電電流の大きさ、および充放電電力の大きさなどを含む。入力部２は、入力した指定値を直流抵抗計算部５に出力する。
【００３８】
ＳＯＣ計算部３は、蓄電デバイスのＳＯＣを計算する部分である。たとえば、蓄電デバイスの初期のＳＯＣ（０）と、その後の蓄電デバイスの充放電電気量とから、蓄電デバイスのＳＯＣ（ｔ）が計算される。蓄電デバイスの初期のＳＯＣ（０）の値は特に限定されず、適宜設定されてよい。蓄電デバイスの充放電電気量は、蓄電デバイスの充放電電流を充放電時間で積算することによって求められる。蓄電デバイスのＳＯＣ（ｔ）は、時刻ｔにおける蓄電デバイスの充放電電気量と蓄電デバイスの満充電容量とに基づいて求められる。時刻ｔのＳＯＣ（ｔ）の計算において、蓄電デバイスに流れる電流として、等価回路モデル４０を時刻０から時刻ｔ−１までに流れた電流Ｉが用いられ得る。この場合、ＳＯＣ計算部３は、たとえば以下の式（３）によってＳＯＣ（ｔ）を計算する。ＳＯＣ計算部３は、計算したＳＯＣ（ｔ）をパラメータ設定部４、分極計算部６、およびＯＣＶ計算部７にそれぞれ出力する。
【数３】
2018009940.tif 000005
【００３９】
パラメータ設定部４は、蓄電デバイスの端子電圧の推定に必要な種々の特性パラメータの値を設定する部分である。特性パラメータは、例えば、抵抗器１１の抵抗値（直流抵抗）、抵抗器２１の抵抗値（第１の分極抵抗）、時定数τ１（第１の分極時定数）、抵抗器２３の抵抗値（第２の分極抵抗）、及び時定数τ２（第２の分極時定数）である。なお、各特性パラメータの値は、蓄電デバイスのＳＯＣに応じて変更されてもよい。
【００４０】
パラメータ設定部４は、たとえば、各パラメータの値を記述するルックアップテーブルを参照することによって、各パラメータの値を設定する。ルックアップテーブルは、パラメータごとに設けられる。ルックアップテーブルは、たとえばＳＯＣと各パラメータの値とが対応付けられたテーブルである。この場合、パラメータ設定部４は、各ルックアップテーブルを参照することによって、ＳＯＣ計算部３から受け取ったＳＯＣ（ｔ）に対応付けられた各パラメータの値を取得し、取得した値を各パラメータの値に設定する。なお、各ルックアップテーブルは、蓄電デバイスの温度ごとに準備されていてもよい。その場合には、さらに、蓄電デバイスの温度も考慮して、各パラメータの値が設定される。また、各パラメータの値は予め定められていてもよい。パラメータ設定部４は、設定した各パラメータの値を直流抵抗計算部５および分極計算部６に出力する。
【００４１】
直流抵抗計算部５は、等価回路モデル４０中の回路１０に発生する直流抵抗電圧Ｖｄｃ（ｔ）を計算する部分である。また、直流抵抗計算部５は、入力部２によって入力された指定値（ｂａｔ＿ｄｅｍａｎｄ）から、等価回路モデル４０に流れる電流Ｉ（ｔ）を計算する部分でもある。分極計算部６は、等価回路モデル４０中の回路２０に発生する分極電圧Ｖｐｏｌ（ｔ）を計算する部分である。ＯＣＶ計算部７は、蓄電デバイスの開放電圧Ｖｏｃｖ（ｔ）を計算する部分である。先に説明したように、開放電圧Ｖｏｃｖ（ｔ）は、蓄電デバイスのＳＯＣから求められる。たとえば、各ＳＯＣの値と開放電圧Ｖｏｃｖの値とを対応付けたテーブルが予め準備されている。ＯＣＶ計算部７は、当該テーブルを参照することによって、ＳＯＣ計算部３から受け取ったＳＯＣ（ｔ）から開放電圧Ｖｏｃｖ（ｔ）を計算する。なお、上述のテーブルが、温度ごとに準備されていてもよく、その場合には、さらに、蓄電デバイスの温度も考慮して、開放電圧Ｖｏｃｖ（ｔ）が計算される。
【００４２】
端子電圧計算部８は、蓄電デバイスの端子電圧Ｖ（ｔ）を計算する部分である。先に説明したように、直流抵抗計算部５によって計算された直流抵抗電圧Ｖｄｃ（ｔ）、分極計算部６によって計算された分極電圧Ｖｐｏｌ（ｔ）、およびＯＣＶ計算部７によって計算された開放電圧Ｖｏｃｖ（ｔ）が端子電圧計算部８に送られる。端子電圧計算部８は、直流抵抗電圧Ｖｄｃ（ｔ）、分極電圧Ｖｐｏｌ（ｔ）、および開放電圧Ｖｏｃｖ（ｔ）に基づいて、端子電圧Ｖ（ｔ）を計算する。具体的には、端子電圧計算部８は、上記式（２）に示されるように、直流抵抗電圧Ｖｄｃ（ｔ）、分極電圧Ｖｐｏｌ（ｔ）、および開放電圧Ｖｏｃｖ（ｔ）を加算し、その合計電圧を端子電圧Ｖ（ｔ）として計算する。端子電圧計算部８は、計算した端子電圧Ｖ（ｔ）をシミュレーション装置１の外部および直流抵抗計算部５に出力する。
【００４３】
次に、図５を参照して、シミュレーション装置１が実行する端子電圧Ｖ（ｔ）の計算処理（シミュレーション方法）を説明する。図５は、シミュレーション装置１が実行する端子電圧Ｖ（ｔ）の計算処理の例を示すフローチャートである。図５に示されるフローチャートの処理は、たとえば燃費シミュレーション装置９０の燃費計算において、ある時刻ｔにおける蓄電デバイスの端子電圧を推定する際に実行される。
【００４４】
まず、入力部２が指定値（ｂａｔ＿ｄｅｍａｎｄ）を入力する（ステップＳ０１）。たとえば、入力部２は、シミュレーション装置１の外部装置から指定値を受け取ることにより、その指定値を入力する。そして、入力部２は、入力した指定値を直流抵抗計算部５に出力する。そして、ＳＯＣ計算部３は、蓄電デバイスのＳＯＣを計算する（ステップＳ０２）。ＳＯＣ計算部３は、たとえば、上述された式（３）を用いてＳＯＣ（ｔ）を計算する。そして、ＳＯＣ計算部３は、計算したＳＯＣ（ｔ）をパラメータ設定部４、分極計算部６、およびＯＣＶ計算部７に出力する。
【００４５】
続いて、パラメータ設定部４は、等価回路モデル４０の各特性パラメータを設定する（ステップＳ０３）。ステップＳ０３において設定される特性パラメータは、たとえば、抵抗器１１の抵抗値、抵抗器２１の抵抗値、時定数τ１、抵抗器２３の抵抗値、及び時定数τ２である。パラメータ設定部４は、たとえば、各特性パラメータの値を記述するルックアップテーブルを参照することによって、ＳＯＣ計算部３から受け取ったＳＯＣ（ｔ）に対応付けられた各パラメータの値を取得し、取得した値を各パラメータの値に設定する。そして、パラメータ設定部４は、設定したパラメータを直流抵抗計算部５および分極計算部６に出力する。
【００４６】
続いて、直流抵抗計算部５は、パラメータ設定部４から提供された抵抗器１１の抵抗値を用いて、電流Ｉ（ｔ）および直流抵抗電圧Ｖｄｃ（ｔ）を計算する（ステップＳ０４）。直流抵抗計算部５は、充放電モードが定電流放電モード（端子電圧Ｖ（ｔ）によらず、一定の電流を流すモード）である場合には、入力部２によって入力された指定値に含まれる指定電流を電流Ｉ（ｔ）に設定する。そして、直流抵抗計算部５は、この電流Ｉ（ｔ）に基づいて直流抵抗電圧Ｖｄｃ（ｔ）を計算する。また、直流抵抗計算部５は、充放電モードが定電圧充電モード（蓄電デバイスを充電するための電圧源（たとえばオルタネータ）の出力電圧を一定にした状態で蓄電デバイスを充電するモード）である場合には、まず直流抵抗電圧Ｖｄｃ（ｔ）を計算する。そして、この直流抵抗電圧Ｖｄｃ（ｔ）に基づいて、等価回路モデル４０に流れる電流Ｉ（ｔ）を計算する。
【００４７】
続いて、分極計算部６は、分極電圧Ｖｐｏｌ（ｔ）を計算する（ステップＳ０５）。具体的には、分極計算部６は、パラメータ設定部４から提供された抵抗器２１の抵抗値、時定数τ１、抵抗器２３の抵抗値、及び時定数τ２を用いて、第１分極電圧Ｖｐ１（ｔ）および第２分極電圧Ｖｐ２（ｔ）を計算する。そして、分極計算部６は、それら第１分極電圧Ｖｐ１（ｔ）および第２分極電圧Ｖｐ２（ｔ）の合計値を、分極電圧Ｖｐｏｌ（ｔ）として計算する。
【００４８】
続いて、ＯＣＶ計算部７は、開放電圧Ｖｏｃｖ（ｔ）を計算する（ステップＳ０６）。たとえば、ＯＣＶ計算部７は、各ＳＯＣの値と開放電圧Ｖｏｃｖの値とを対応付けたテーブルを参照することによって、ＳＯＣ計算部３から受け取ったＳＯＣ（ｔ）から開放電圧Ｖｏｃｖ（ｔ）を計算する。そして、ＯＣＶ計算部７は、計算した開放電圧Ｖｏｃｖ（ｔ）を端子電圧計算部８に出力する。
【００４９】
続いて、端子電圧計算部８は、端子電圧Ｖ（ｔ）を計算する（ステップＳ０７）。具体的には、端子電圧計算部８は、直流抵抗計算部５によって計算された直流抵抗電圧Ｖｄｃ（ｔ）、分極計算部６によって計算された分極電圧Ｖｐｏｌ（ｔ）、およびＯＣＶ計算部７によって計算された開放電圧Ｖｏｃｖ（ｔ）に基づいて、端子電圧Ｖ（ｔ）を計算する。より具体的には、端子電圧計算部８は、上記式（２）に示されるように、直流抵抗電圧Ｖｄｃ（ｔ）、分極電圧Ｖｐｏｌ（ｔ）、および開放電圧Ｖｏｃｖ（ｔ）を加算し、その合計電圧を端子電圧Ｖ（ｔ）として計算する。そして、端子電圧計算部８は、計算した端子電圧Ｖ（ｔ）をシミュレーション装置１の外部、および直流抵抗計算部５に出力する。以上のようにして、時刻ｔにおける端子電圧Ｖ（ｔ）の計算処理が終了する。
【００５０】
なお、ステップＳ０５の処理とステップＳ０６の処理とは、並行して行われてもよく、実施される順番が逆になってもよい。
【００５１】
ここで、等価回路モデル４０を構成する特性パラメータについて詳細に説明する。前述したように、等価回路モデル４０は、例えば抵抗器１１の抵抗値（直流抵抗）、抵抗器２１の抵抗値（第１の分極抵抗）、時定数τ１（第１の分極時定数）、抵抗器２３の抵抗値（第２の分極抵抗）、及び時定数τ２（第２の分極時定数）といった複数の特性パラメータを有する。通常、これらの特性パラメータは、蓄電デバイスの使用による劣化を考慮せずに算出される。その場合、劣化したときの蓄電デバイスの入出力特性を精度よく推定することができないという問題がある。
【００５２】
そこで、本実施形態では、複数の特性パラメータのうち少なくとも一つの特性パラメータを以下の数式（４）のように設定する。すなわち、任意の特性パラメータＡを
【数４】
2018009940.tif 000006
としてモデル化する。右辺の第１項Ａ0は蓄電デバイス未使用時に対応する特性パラメータの初期値であり、右辺の第２項Ａ1は、蓄電デバイスを或る時間使用した後における特性パラメータの変化分である。この変化分は、蓄電デバイスの劣化の影響を表す時間関数であり、使用期間内の適切な時間を入力することによって、当該使用時間経過後に対応する劣化した特性パラメータが得られる。なお、この時間は、例えば数百時間ないし数千時間といった長さである。
【００５３】
特性パラメータＡを数式（４）のように表現することにより、例えば未使用状態の蓄電デバイスを用いて初期特性パラメータＡ0を同定するだけで、任意の期間経過後の特性パラメータＡを計算によって求めることができる。
【００５４】
上記の第２項Ａ1は、蓄電デバイスの種類によって異なってもよい。蓄電デバイスが鉛蓄電池である場合、第２項Ａ1は例えば次の数式（５）のように定義される。なお、係数ａ1、ａ2、及びａ3は定数であり、ＳＯＣ1（ｔ）は暗電流放電時のＳＯＣであり、ＳＯＣ2（ｔ）は休止時のＳＯＣである。これらのＳＯＣは時間ｔの関数となっている。
【数５】
2018009940.tif 000007
【００５５】
数式（５）に示されるＡ1は、３つの項を含んでいる。第１項は、
【数６】
2018009940.tif 000008
である。この項は、電流Ｉ（ｔ）の絶対値の時間積分と、該時間積分に乗算された係数ａ1とを含む。この場合、電流Ｉ（ｔ）は蓄電デバイスの動作状態を表す数値であって、電流Ｉ（ｔ）の絶対値の時間積分は即ち、使用期間中に蓄電デバイスを流れる総電流量を表し、係数ａ1は、総電流量に対する蓄電デバイスの劣化の速度を表す。従って、数式（６）で表される項は、蓄電デバイスを流れる総電流量に基づく劣化（以下、通電劣化という）を表す。なお、ｔ１は充放電時間である。また、蓄電デバイスの特性パラメータは複数あり、劣化の進行に伴ってこれらの特性パラメータの値が変化するが、通電劣化による変化分は、特性パラメータによって増加傾向の場合と減少傾向の場合とがある。従って、特性パラメータごとに係数ａ1の符号が決定される。
【００５６】
第２項は、
【数７】
2018009940.tif 000009
である。この項は、電流Ｉ（ｔ）の絶対値の時間積分と、該時間積分に乗算された係数ａ2及び放電深度（ＤＯＤ）とを含む。ＤＯＤは、例えば定数である。或いは、ＤＯＤは時間ｔの関数であってもよく、その場合、Ｉ（ｔ）の絶対値とＤＯＤ（ｔ）との積が時間積分される。また、係数ａ2は、ＤＯＤに対する蓄電デバイスの劣化の速度を表す。従って、数式（７）で表される項は、蓄電デバイスのＤＯＤに基づく劣化（ＤＯＤ劣化）を表す。なお、ＤＯＤ劣化による特性パラメータの変化分は、特性パラメータによって増加傾向の場合と減少傾向の場合とがある。従って、特性パラメータごとに係数ａ２の符号が決定される。
【００５７】
第３項及び第４項は、
【数８】
2018009940.tif 000010
である。第３項は、ＳＯＣ1（ｔ）の時間積分と、この時間積分に乗算された係数ａ3とを含む。また、第４項は、ＳＯＣ2（ｔ）の時間積分と、この時間積分に乗算された、第３項と共通の係数ａ3とを含む。ＳＯＣ1（ｔ）及びＳＯＣ2（ｔ）は蓄電デバイスの動作状態を表す数値である。係数ａ3は、暗電流放電時及び休止時の各ＳＯＣに対する蓄電デバイスの劣化の速度を表す。数式（８）で表される第３項及び第４項は、それぞれ蓄電デバイスの暗電流放電時及び休止時の各ＳＯＣに基づく劣化（暗電流放電劣化および休止劣化）を表す。暗電流放電劣化および休止劣化による特性パラメータの変化分は、特性パラメータによって増加傾向の場合と減少傾向の場合とがある。従って、特性パラメータごとに係数ａ３の符号が決定される。なお、暗電流放電とは、車両のエンジンが停止している状態（すなわちオルタネータが発電していない状態。但し、アイドリングストップ時を除く）において、カーナビゲーションシステム及び時計などに供給される微弱な電流をいい、暗電流放電時とはこのような微弱が流れている期間をいう。また、休止とは、暗電流が全く流れていない状態をいい、休止時とは休止状態である期間をいう。なお、ｔ２は暗電流放電時間であり、ｔ３は休止時間である。数式（８）においては、必要に応じて、ＳＯＣ1（ｔ）及びＳＯＣ2（ｔ）の各時間積分項の一方（すなわち第３項及び第４項の一方）を省略してもよい。
【００５８】
なお、通常の充放電時においてはＳＯＣ1（ｔ）＝０、ＳＯＣ2（ｔ）＝０である。暗電流放電時においてはＩ（ｔ）＝０、ＳＯＣ2（ｔ）＝０である。休止時においてはＩ（ｔ）＝０、ＳＯＣ1（ｔ）＝０である。
【００５９】
また、蓄電デバイスがリチウムイオン電池である場合、数式（４）の第２項Ａ1は例えば次の数式（９）のように定義される。
【数９】
2018009940.tif 000011
【００６０】
数式（９）に示されるＡ1は、３つの項を含んでいる。第１項は通電劣化項である。第２項は、
【数１０】
2018009940.tif 000012
である。この項は、電流Ｉ（ｔ）の絶対値の時間積分と、該時間積分に乗算された係数ａ2と、該時間積分に乗算された自己発熱とを含む。すなわち、次の数式（１１）によって表される時間τの関数Ｔは、蓄電デバイスの自己発熱を表す。
【数１１】
2018009940.tif 000013
自己発熱Ｔは、数式（１１）に示されるように、電流Ｉ（τ）の二乗の時間積分によって求められる。また、係数ａ2は、自己発熱Ｔに対する蓄電デバイスの劣化の速度を表す。従って、数式（１０）で表される項は、蓄電デバイスの自己発熱に基づく劣化（以下、自己発熱劣化という）を表す。
【００６１】
第３項は、
【数１２】
2018009940.tif 000014
である。この項は、休止時のＳＯＣ2（ｔ）の時間積分と、該時間積分に乗算された係数ａ3とを含む。これらの意味付けは、前述した鉛蓄電池の場合の第４項と同様である。すなわち、この項は休止劣化項である。
【００６２】
上記の例では、数式（８）において暗電流放電劣化項が存在するが、数式（１２）においては同様の暗電流放電劣化項が存在しない。それは次の理由による。例えば蓄電デバイスが単一で用いられる場合には、車両のエンジンが停止しても暗電流が流れる期間が必ず存在する。このような単一で用いられる蓄電デバイスとしては、鉛蓄電池が挙げられる。これに対し、例えばμＨＥＶ方式のようにメイン蓄電デバイス及びサブ蓄電デバイスが用いられる場合には、メイン蓄電デバイスのみから暗電流が供給され、サブ蓄電デバイスからは暗電流が供給されない状況が考えられる。そのような状況では、サブ蓄電デバイスにおいて暗電流状態は生じない。リチウムイオン電池は、このようなサブ蓄電デバイスとして用いられることが多い。従って、数式（１２）においては暗電流放電劣化項が省略されている。
【００６３】
また、蓄電デバイスがニッケル亜鉛電池である場合、数式（４）の第２項Ａ1は例えば次の数式（１３）のように定義される。
【数１３】
2018009940.tif 000015
数式（１３）に示されるＡ1は、３つの項を含んでいる。第１項は通電劣化項である。第２項はＤＯＤ劣化項である。第３項は休止劣化項である。ニッケル亜鉛電池もまた、上述したリチウムイオン電池と同様に、サブ蓄電デバイスとして用いられることが多い。従って、数式（１３）においては暗電流放電劣化項が省略されている。
【００６４】
ここで、上述した係数ａ1，ａ2，ａ3の算出方法について説明する。３つの係数ａ1，ａ2，ａ3を算出するためには、３つの異なるモデル式が必要となる。そこで、互いに時間波形が異なる３つの電流Ｉ（ｔ）を特性パラメータの式に入力して３つのモデル式を立て、それらに基づいて係数ａ1，ａ2，ａ3を求める。図６（ａ）〜図６（ｃ）は、そのような電流Ｉ（ｔ）の波形の例を概念的に示すグラフである。これらの図において、縦軸は電流Ｉ（ｔ）を表し、横軸は時間を表す。これらの電流波形は、定電圧充電期間Ｔａ、及び定電流放電期間Ｔｂ、及び休止期間Ｔｃを含む。図６（ａ）は、定電圧充電期間Ｔａにおける電流が比較的大きい場合を示し、充電電流が大きい分だけ定電流放電期間Ｔｂが長くなっている。図６（ｂ）は、定電圧充電期間Ｔａにおける電流が比較的小さい場合を示し、充電電流が小さい分だけ定電流放電期間Ｔｂが短くなっている。図６（ｃ）は、図６（ａ）及び図６（ｂ）と較べて休止期間Ｔｃが長い場合を示している。
【００６５】
図７は、リチウムイオン電池を例として図６（ａ）〜図６（ｃ）の電流波形を入力したときの、複数の特性パラメータのうち直流抵抗（図２の抵抗器１１の抵抗値）の変化を示すグラフである。縦軸は直流抵抗（単位：ｍΩ）を表し、横軸はトータルサイクル時間（使用時間、単位：時間）を示す。また、同図において、菱形のプロットＰ１は図６（ａ）に示された電流波形を入力した場合を示し、正方形のプロットＰ２は図６（ｂ）に示された電流波形を入力した場合を示し、三角形のプロットＰ３は図６（ｃ）に示された電流波形を入力した場合を示す。図７に示されるように、入力される電流Ｉ（ｔ）の波形が異なると、それに伴い直流抵抗の劣化分（数式（４）のＡ1）が変化することがわかる。
【００６６】
図８は、図６（ａ）〜図６（ｃ）の電流波形に対応する、通電劣化項の時間積分の計算値をプロットしたグラフである。図９は、図６（ａ）〜図６（ｃ）の電流波形に対応する、自己発熱劣化項の時間積分の計算値をプロットしたグラフである。図１０は、図６（ａ）〜図６（ｃ）の電流波形に対応する、休止劣化項の時間積分の計算値をプロットしたグラフである。なお、図８〜図１０において、縦軸は時間積分値を表し、横軸はトータルサイクル時間（単位：時間）を表す。また、これらの図において、菱形のプロットＰ４は図６（ａ）に示された電流波形を入力したときの数値を示し、正方形のプロットＰ５は図６（ｂ）に示された電流波形を入力したときの数値を示し、三角形のプロットＰ６は図６（ｃ）に示された電流波形を入力したときの数値を示す。
【００６７】
図８〜図１０に示されるように、図６（ａ）〜図６（ｃ）の電流波形に基づいて、通電劣化項、自己発熱劣化項、及び休止劣化項それぞれにおける時間積分の値が計算される。従って、係数ａ1，ａ2，ａ3を変数として含む互いに独立した３つの関数を作成でき、実験値とのカーブフィッティングによる最適化を行うことで、係数ａ1，ａ2，ａ3を求めることができる。図１１（ａ）〜図１１（ｃ）は、カーブフィッティングの様子を概念的に示すグラフである。図１１（ａ）〜図１１（ｃ）は、それぞれ図６（ａ）〜図６（ｃ）の電流波形により得られたａ1，ａ2，ａ3の関数と実験値とのフィッティングの様子を示している。図中のプロットＰ７〜Ｐ９が実験値であり、曲線Ｒ１〜Ｒ３が関数からの推定値である。
【００６８】
図１２は、上述した方法によって得られた、各特性パラメータの係数ａ1，ａ2，ａ3の数値例を示す図表である。図１２に示されるように、上述した方法によって係数ａ1，ａ2，ａ3が好適に求められる。なお、図１２には、フィッティング誤差（％）が併せて示されている。第１分極抵抗及び第２分極抵抗の誤差が比較的大きくなっているが、これは、試験期間が短く、劣化の進行度合いが小さい段階であることが原因と考えられる。
【００６９】
以上に説明した本実施形態によるシミュレーション方法およびシミュレーション装置１によって得られる効果について説明する。本実施形態のシミュレーション方法及びシミュレーション装置１では、数式（４）に示したように、少なくとも一つの特性パラメータＡが、蓄電デバイスの劣化の影響を表す時間関数Ａ1を含む。このような時間関数Ａ1に蓄電デバイスの使用時間（トータルサイクル時間）を入力することにより、該使用時間経過時における蓄電デバイスの劣化度合いを特性パラメータＡに反映させることができる。
【００７０】
そして、本発明者の知見によれば、暗電流放電時または休止時のＳＯＣによって蓄電デバイスの劣化の度合いは大きく変化する。本実施形態では、時間関数Ａ1が、劣化速度を表す係数ａ3と、蓄電デバイスの暗電流放電時のＳＯＣ（ＳＯＣ1（ｔ））及び休止時のＳＯＣ（ＳＯＣ2（ｔ））の一方又は双方の時間積分とを乗算した項（例えば数式（８）、数式（１２）など）を含む。これにより、蓄電デバイスの暗電流放電時または休止時の劣化度合いを精度良く表すことができる。従って、本実施形態によれば、実際に劣化した蓄電デバイスを用いなくても、劣化したときの蓄電デバイスの入出力特性を精度よく推定することができ、劣化した蓄電デバイスを用いた燃費シミュレーションなどを精度良く行うことができる。
【００７１】
図１３は、車両の燃費シミュレーションにおいて、リチウムイオン電池の等価回路モデルに本実施形態の特性パラメータを使用した場合（グラフＧ１１）と、劣化による影響を考慮しない特性パラメータ（すなわち数式（４）の右辺第２項Ａ1がないもの）を使用した場合（グラフＧ１２）とにおける燃費誤差の最大値を比較した結果を示すグラフである。図１３において、縦軸は燃費誤差の最大値（単位：％）を表し、横軸は試験期間（単位：日）を表す。なお、特性パラメータの初期値（すなわち数式（４）の右辺第１項Ａ0）の同定に用いる電流Ｉ（ｔ）としては、図１４に示される電流波形を用いた。この電流波形は、定電圧充電期間Ｔ１、定電圧充電期間Ｔ１後の定電流放電期間Ｔ２、及び定電流放電期間Ｔ２後のクランキング期間Ｔ３を含む第１〜第３の期間ＴＡ〜ＴＣを繰り返し含んでいる。なお、これら第１〜第３の期間ＴＡ〜ＴＣにおける定電圧充電期間Ｔ１の電圧値は１４（Ｖ）で一定であり、定電流放電期間Ｔ２及びクランキング期間Ｔ３の時間はそれぞれ５９秒、１秒で一定である。また、第１〜第３の期間ＴＡ〜ＴＣにおける電流値は次の通りである。
＜第１の期間ＴＡ＞
定電圧充電期間Ｔ１：１００（Ａ）
定電流放電期間Ｔ２：−２０（Ａ）
クランキング期間Ｔ３：−３００（Ａ）
＜第２の期間ＴＢ＞
定電圧充電期間Ｔ１：２００（Ａ）
定電流放電期間Ｔ２：−４５（Ａ）
クランキング期間Ｔ３：−３００（Ａ）
＜第３の期間ＴＣ＞
定電圧充電期間Ｔ１：５０（Ａ）
定電流放電期間Ｔ２：−１０（Ａ）
クランキング期間Ｔ３：−３００（Ａ）
【００７２】
また、燃費誤差の最大値とは、特性パラメータを含む等価回路モデルを用いて燃費シミュレーションを行った結果と、実際に測定された燃費との差の最大値である。燃費誤差は、電池の特性パラメータ抽出時の電圧誤差（蓄電デバイスの端子電圧の実測値とモデルによる端子電圧の推定値との差）に依存する。特性パラメータは複数存在し、それらの値の組合せに応じて電圧誤差は変化する。従って、同じ電圧誤差の値となる特性パラメータの値の組合せは無数に存在するので、同じ電圧誤差の値であっても特性パラメータの値の組合せによって燃費計算結果は異なる。上述した「最大値」とは、様々な特性パラメータの値の組合せに対する電圧誤差の最大値に対応する、燃費誤差の最大値を意味する。
【００７３】
図１３に示されるように、試験期間が長くなるほど燃費誤差は大きくなるが、劣化による影響を考慮しない場合には、試験期間が６０日を過ぎると燃費誤差の最大値が０．３％を超えている。これに対し、休止時等の劣化による影響を考慮した本実施形態では、試験期間が６０日を過ぎても燃費誤差の最大値が０．０５％以下に収まっている。このように、本実施形態の方法および装置によれば、試験期間が長くなるほど、蓄電デバイスの劣化の状態を精度よく燃費シミュレーション結果に反映させることができる。
【００７４】
本実施形態による燃費シミュレーション結果は、例えば、車両に採用される蓄電デバイス容量の選択に応用することができる。一般に、車両の燃費は搭載する蓄電デバイスの容量が大きいほど良くなる。一方、使用開始からの年数を経るほど、蓄電デバイスの性能が劣化し、車両の燃費は低下する。従来は、シミュレーションにおいて蓄電デバイスの劣化により性能がどれほど低下するかが不明であったため、燃費シミュレーション結果から選択される蓄電デバイス容量よりも十分に余裕のある蓄電デバイス容量が選択されていた。このような選択方法では、蓄電デバイス容量が必要以上に大きくなり易く、車両コスト低減の妨げとなるおそれがある。
【００７５】
そのような問題に対し、本実施形態のシミュレーション方法およびシミュレーション装置１では、蓄電デバイスの劣化度合いに応じたシミュレーションを行うことができるので、使用開始からの年数を考慮した燃費シミュレーションを精度良く行うことができる。従って、所定の年数を経た後の推定燃費に基づいて、所定の燃費条件を満足できる蓄電デバイス容量を的確に選択することができる。
【００７６】
図１５は、本実施形態による燃費シミュレーション結果の一例を示すグラフである。図１５において、縦軸は燃費（単位：ｋｍ／ｌ）を表し、横軸は使用年数（単位：年）を表す。また、図中の菱形のプロットＰ１０、正方形のプロットＰ１１、及び三角形のプロットＰ１２は、蓄電デバイス初期容量がそれぞれ３Ａｈ、５Ａｈ、及び７Ａｈである場合を示す。図１５に示されるように、車両の燃費は蓄電デバイス容量が大きいほど良いが、使用年数が長くなるほど車両の燃費は低下する。そこで、例えば使用年数が５年経過した時点での燃費を３０（ｋｍ／ｌ）以上としたい場合、このグラフによれば、車両に搭載する蓄電デバイス初期容量を５Ａｈとすれば良いことがわかる。このように、本実施形態によれば、所定の年数を経た後の推定燃費に基づいて、所定の燃費条件を満足できる蓄電デバイス容量を的確に選択することができる。
【００７７】
なお、蓄電デバイス容量の選択は、推定燃費を基準として行う場合に限られない。本実施形態によれば、所定の年数を経た後の蓄電デバイスの推定特性に基づいて、所定の条件を満足できる蓄電デバイス容量を的確に選択することができる。
【００７８】
（変形例）
上記実施形態では、係数ａ1，ａ2，ａ3と蓄電デバイスの温度との関係については述べていない。すなわち、係数ａ1，ａ2，ａ3は、蓄電デバイスの温度によらず一定であってもよい。しかしながら、多くの場合において、好適な係数ａ1，ａ2，ａ3の値は蓄電デバイスの温度に依存する。従って、係数ａ1，ａ2，ａ3は、蓄電デバイスの温度に応じて変化してもよい。これにより、蓄電デバイスの温度に応じて変化する劣化度合いを精度良く表すことができる。
【００７９】
具体的には、係数ａ1，ａ2，ａ3が温度ＴＨの関数ａ1（ＴＨ），ａ2（ＴＨ），ａ3（ＴＨ）であってもよく、或いは、複数の温度毎に異なる係数ａ1，ａ2，ａ3が設定されてもよい。そのために、実験値とのカーブフィッティングによる最適化を行う際に、蓄電デバイスの温度を変えながら実験値を取得するとよい。
【００８０】
本発明によるシミュレーション方法及びシミュレーション装置は、上述した実施形態及び変形例に限られるものではなく、他に様々な変形が可能である。例えば、上述した実施形態では、蓄電デバイスの動作状態を表し時間積分される数値として、蓄電デバイスを流れる電流Ｉ（ｔ）、暗電流放電時のＳＯＣ1（ｔ）、及び休止時のＳＯＣ2（ｔ）を例示したが、本発明における当該数値としては、蓄電デバイスの動作状態を表すものであればこれら以外にも様々な数値を採用し得る。例えば、ＳＯＣ1（ｔ）、ＳＯＣ2（ｔ）に代わる数値として、蓄電デバイスの端子電圧を用いてもよい。
【符号の説明】
【００８１】
１…シミュレーション装置、２…入力部、３…ＳＯＣ計算部、４…パラメータ設定部、５…直流抵抗計算部、６…分極計算部、７…ＯＣＶ計算部、８…端子電圧計算部、１０，２０…回路、１１，２１，２３…抵抗器、２２，２４…コンデンサ、３０…定電圧源、４０…等価回路モデル、９０…燃費シミュレーション装置、９１…入力部、９２…制御部、９３…出力部、Ｎ１，Ｎ２…ノード。
【図１】
2018009940.tif 000016
【図２】
2018009940.tif 000017
【図３】
2018009940.tif 000018
【図４】
2018009940.tif 000019
【図５】
2018009940.tif 000020
【図６】
2018009940.tif 000021
【図７】
2018009940.tif 000022
【図８】
2018009940.tif 000023
【図９】
2018009940.tif 000024
【図１０】
2018009940.tif 000025
【図１１】
2018009940.tif 000026
【図１２】
2018009940.tif 000027
【図１３】
2018009940.tif 000028
【図１４】
2018009940.tif 000029
【図１５】
2018009940.tif 000030