(57)【要約】【課題】  電解液中の塩濃度の偏りに伴う劣化を抑制する。【解決手段】  二次電池（１０）の抵抗上昇量（Ｒｈ）が閾値（Ｒ＿ｔｈ）以上であるとき、抵抗上昇量および閾値の差を示す抵抗上昇抑制量（ΔＲｈ）を算出する。抵抗上昇量と、正極板および負極板における塩濃度の差と、発電要素内を移動する電解液の流速との対応関係を用いて、抵抗上昇量を抵抗上昇抑制量だけ低下させるときの流速である流速制限値（ｕ＿ｌｉｍ）を算出する。二次電池の温度から特定される、二次電池の電流値および流速の対応関係を用いて、流速制限値に対応した電流値である電流制限値（Ｉｂ＿ｍａｘ，Ｉｂ＿ｍｉｎ）を算出し、電流制限値に応じた上限電力値（Ｗｏｕｔ，Ｗｉｎ）を設定する。電流制限値に応じた上限電力値を設定して、電解液の流速を流速制限値に制限することにより、抵抗上昇抑制量の分だけ、抵抗上昇量を低下させて、塩濃度の偏りに伴う劣化を抑制できる。【選択図】  図１０
2016039010.tif 000002
【特許請求の範囲】
【請求項１】
充放電を行う発電要素および電解液が電池ケースに収容された二次電池と、
前記二次電池の温度を検出する温度センサと、
前記二次電池の充放電を制御するときの上限電力値を設定するコントローラと、を有し、
前記発電要素は、正極板と、負極板と、前記正極板および前記負極板の間に配置されたセパレータとを有しており、
前記コントローラは、
前記発電要素に含浸された前記電解液中の塩濃度の偏りに伴って発生する前記二次電池の抵抗上昇量が閾値以上であるとき、前記抵抗上昇量および前記閾値の差を示す抵抗上昇抑制量を算出し、
前記抵抗上昇量と、前記正極板および前記負極板における前記塩濃度の差と、前記発電要素内を移動する前記電解液の流速との対応関係を用いて、前記抵抗上昇量を前記抵抗上昇抑制量だけ低下させるときの前記流速である流速制限値を算出し、
前記二次電池の温度から特定される、前記二次電池の電流値および前記流速の対応関係を用いて、前記流速制限値に対応した前記電流値である電流制限値を算出し、
前記電流制限値に応じた前記上限電力値を設定する、
ことを特徴とする電池システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、電解液中の塩濃度の偏りに伴う劣化を抑制する電池システムに関する。
【背景技術】
【０００２】
特許文献１では、電解液中の塩濃度の偏りに伴う劣化を評価するための評価値を設定し、組電池の電流値に基づいて評価値を算出している。ここで、評価値が目標値よりも大きいとき、評価値および目標値の差分を積算し、この積算値が閾値よりも大きいとき、組電池の出力制限値を最大値よりも小さい値に設定している。これにより、塩濃度の偏りに伴う劣化が進行することを抑制するようにしている。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１３−２２５３９７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
電解液中の塩濃度の偏りに伴う劣化は、塩濃度の偏りだけでなく、正極、負極およびセパレータに含浸された電解液の流れ（流速）にも依存することが分かった。このため、塩濃度の偏りに伴う劣化を抑制するためには、電解液の流れを考慮する必要がある。
【課題を解決するための手段】
【０００５】
本発明の電池システムは、二次電池と、二次電池の温度を検出する温度センサと、二次電池の充放電を制御するときの上限電力値を設定するコントローラと、を有する。二次電池では、充放電を行う発電要素および電解液が電池ケースに収容されており、発電要素は、正極板と、負極板と、正極板および負極板の間に配置されたセパレータとを有する。
【０００６】
コントローラは、発電要素に含浸された電解液中の塩濃度の偏りに伴って発生する二次電池の抵抗上昇量が閾値以上であるとき、抵抗上昇量および閾値の差を示す抵抗上昇抑制量を算出する。コントローラは、抵抗上昇量と、正極板および負極板における塩濃度の差と、発電要素内を移動する電解液の流速との対応関係を用いて、抵抗上昇量を抵抗上昇抑制量だけ低下させるときの流速である流速制限値を算出する。コントローラは、二次電池の温度から特定される、二次電池の電流値および流速の対応関係を用いて、流速制限値に対応した電流値である電流制限値を算出し、電流制限値に応じた上限電力値を設定する。
【０００７】
本発明によれば、電流制限値に応じた上限電力値を設定して、電解液の流速を流速制限値に制限することにより、抵抗上昇抑制量の分だけ、抵抗上昇量を低下させて、塩濃度の偏りに伴う劣化を抑制できる。
【０００８】
ここで、抵抗上昇量は、正極板および負極板における塩濃度の差と、発電要素内を移動する電解液の流速とに依存するため、所定の塩濃度差が発生しているとき、抵抗上昇量を抵抗上昇抑制量の分だけ低下させるための流速（流速制限値）を算出できる。また、電解液の流速は、二次電池の電流値と、電解液の膨張／収縮とに依存し、電解液の膨張／収縮は、二次電池の温度に依存する。すなわち、電解液の流速は、二次電池の電流値および温度に依存する。したがって、二次電池の温度を検出すれば、流速制限値に対応した電流制限値を算出できる。
【図面の簡単な説明】
【０００９】
【図１】電池システムの構成を示す図である。
【図２】二次電池の構成を示す図である。
【図３】発電要素の展開図である。
【図４】発電要素の外観図である。
【図５】負極板および正極板が対向する方向における塩濃度の偏りを説明する図である。
【図６】負極板および正極板のそれぞれの表面上における塩濃度の偏りを説明する図である。
【図７】電解液の膨張および収縮によって電解液が流れる方向を説明する図である。
【図８】抵抗上昇量および塩濃度差の関係を示す図である。
【図９】塩濃度分布を算出する処理を説明するフローチャートである。
【図１０】二次電池の充放電を制御する処理を説明するフローチャートである。
【図１１】正負極間の塩濃度差が所定値であるとき、電解液の流速および抵抗上昇量の関係を示す図である。
【図１２】電解液の流速および電流値の関係を示す図である。
【図１３】電解液の流速および電流値の関係を示す図である。
【図１４】バネモデルを用いて二次電池の厚さの変化を説明する図である。
【発明を実施するための形態】
【００１０】
以下、本発明の実施例について説明する。
【実施例１】
【００１１】
本発明の電池システムについて、図１を用いて説明する。二次電池１０は、正極ラインＰＬおよび負極ラインＮＬを介して、負荷２０と接続されている。負荷２０は、二次電池１０から出力された電力を受けて動作する。また、負荷２０は、発電を行うこともでき、負荷２０によって生成された電力は、二次電池１０に供給される。これにより、二次電池１０が充電される。
【００１２】
図１に示す電池システムは、例えば、車両に搭載できる。この場合には、複数の二次電池１０を直列に接続した組電池を車両に搭載できる。また、負荷２０としては、モータ・ジェネレータを用いることができる。モータ・ジェネレータは、二次電池１０から出力された電力を受けて、車両を走行させるための動力を生成する。モータ・ジェネレータが生成した動力は、車輪に伝達される。モータ・ジェネレータは、車両の制動時に発生する運動エネルギを電力に変換し、この電力を二次電池１０に供給できる。
【００１３】
電圧センサ３１は、二次電池１０の電圧値Ｖｂを検出し、検出結果をコントローラ４０に出力する。電流センサ３２は、二次電池１０の電流値Ｉｂを検出し、検出結果をコントローラ４０に出力する。本実施例において、二次電池１０を放電しているときの電流値Ｉｂを正の値とし、二次電池１０を充電しているときの電流値Ｉｂを負の値としている。
【００１４】
コントローラ４０は、電圧値Ｖｂおよび電流値Ｉｂに基づいて、二次電池１０の充放電を制御できる。ここで、コントローラ４０は、電圧値Ｖｂや電流値Ｉｂに基づいて、二次電池１０のＳＯＣ（State of Charge）を算出できる。ＳＯＣとは、満充電容量に対する現在の充電容量の割合である。ＳＯＣの算出方法としては、公知の方法を適宜採用できるため、ＳＯＣの算出方法に関する詳細な説明は省略する。
【００１５】
温度センサ３３は、二次電池１０の温度（電池温度）Ｔｂを検出し、検出結果をコントローラ４０に出力する。コントローラ４０は、メモリ４１を有する。メモリ４１は、コントローラ４０が所定の処理（特に、本実施例で説明する処理）を行うときに用いられる情報を記憶する。なお、メモリ４１は、コントローラ４０の外部に設けることもできる。
【００１６】
次に、二次電池１０の構造について、図２を用いて説明する。図２において、Ｘ軸およびＺ軸は、互いに直交する軸である。本実施例では、鉛直方向に相当する軸をＺ軸としている。なお、Ｘ軸およびＺ軸と直交する軸をＹ軸とする。
【００１７】
二次電池１０は、電池ケース１１０および発電要素１２０を有する。電池ケース１１０は、発電要素１２０を収容している。電池ケース１１０は密閉状態となっており、電池ケース１１０の内部には電解液が注入されている。電池ケース１１０には負極端子１１１および正極端子１１２が固定されている。負極端子１１１および正極端子１１２は、発電要素１２０と電気的に接続されている。
【００１８】
発電要素１２０は、充放電を行う要素であり、図３に示すように、負極板（本発明の電極板に相当する）１２１と、正極板（本発明の電極板に相当する）１２２と、セパレータ１２３とを有する。図３は、発電要素１２０の一部を展開した図である。負極板１２１は、集電箔１２１ａと、集電箔１２１ａの表面に形成された負極活物質層１２１ｂとを有する。負極活物質層１２１ｂは、負極活物質、導電剤、バインダーなどを含んでいる。負極活物質層１２１ｂは、集電箔１２１ａの一部の領域に形成されており、集電箔１２１ａの残りの領域には、負極活物質層１２１ｂが形成されていない。
【００１９】
正極板１２２は、集電箔１２２ａと、集電箔１２２ａの表面に形成された正極活物質層１２２ｂとを有する。正極活物質層１２２ｂは、正極活物質、導電剤、バインダーなどを含んでいる。正極活物質層１２２ｂは、集電箔１２２ａの一部の領域に形成されており、集電箔１２２ａの残りの領域には、正極活物質層１２２ｂが形成されていない。
【００２０】
負極活物質層１２１ｂ、正極活物質層１２２ｂおよびセパレータ１２３には、電解液が含浸している。この電解液は、発電要素１２０の内部に存在する。一方、発電要素１２０の外部、言い換えれば、発電要素１２０および電池ケース１１０の間に形成されたスペースにも、余剰液としての電解液が存在している。
【００２１】
図３に示す順番で、負極板１２１、正極板１２２およびセパレータ１２３を積層し、この積層体をＸ軸の周りで図４に示す矢印Ｄの方向に巻くことにより、発電要素１２０が構成される。ここで、負極板１２１および正極板１２２の間には、セパレータ１２３が配置される。
【００２２】
Ｘ軸が延びる方向（Ｘ方向という）における発電要素１２０の一端では、負極板１２１の集電箔１２１ａだけが巻かれている。集電箔１２１ａだけが巻かれた部分は、図２に示す負極端子１１１と電気的に接続される。また、Ｘ方向における発電要素１２０の他端では、正極板１２２の集電箔１２２ａだけが巻かれている。集電箔１２２ａだけが巻かれた部分は、図２に示す正極端子１１２と電気的に接続される。
【００２３】
本実施例では、上述したように、積層体を巻くことにより、発電要素１２０を構成しているが、これに限るものではない。具体的には、積層体を巻かずに、負極板１２１、正極板１２２およびセパレータ１２３を積層しただけで、発電要素１２０を構成することもできる。
【００２４】
図４に示す領域Ａは、負極活物質層１２１ｂおよび正極活物質層１２２ｂがセパレータ１２３を挟んで互いに向かい合う領域である。領域Ａにおいて、二次電池１０（発電要素１２０）の充放電に応じた化学反応が行われる。
【００２５】
二次電池１０では、電解液中の塩濃度に偏りが発生することにより、二次電池１０の内部抵抗値が上昇してしまう。このような内部抵抗値の上昇量を抵抗上昇量Ｒｈとする。抵抗上昇量Ｒｈは、二次電池１０の劣化に伴う内部抵抗値の上昇量とは異なる。劣化に伴う内部抵抗値の上昇量は増加するだけであり、減少することはない。一方、抵抗上昇量Ｒｈは、塩濃度の偏りに依存するため、塩濃度が偏るほど、抵抗上昇量Ｒｈが増加し、塩濃度の偏りが緩和されるほど、抵抗上昇量Ｒｈが減少する。
【００２６】
塩濃度の偏りの状態としては、図５および図６に示す状態がある。図５は、負極板１２１および正極板１２２が対向する方向（Ｙ方向）において、塩濃度の偏りが発生する状態を示す。図５では、負極板１２１、正極板１２２およびセパレータ１２３の位置関係を示す概略図（図５の上側の図）と、塩濃度の分布（一例）を示す図（図５の下側の図）とを表している。
【００２７】
塩濃度分布を示す図において、縦軸は塩濃度であり、横軸はＹ方向における位置である。図５（上側の図）では、負極板１２１および正極板１２２がセパレータ１２３から離れているが、実際には、負極板１２１および正極板１２２がセパレータ１２３に接触している。図５（下側の図）に示すように、二次電池１０の充電時では、塩濃度の偏りとして、実線で示す塩濃度分布が発生することがある。また、二次電池１０の放電時では、塩濃度の偏りとして、一点鎖線で示す塩濃度分布が発生することがある。
【００２８】
なお、図５では、Ｙ方向における塩濃度の偏りを示しているが、これに限るものではない。上述したように、本実施例の発電要素１２０では、負極板１２１および正極板１２２がＸ軸の周りで巻かれているため、負極板１２１（負極活物質層１２１ｂ）および正極板１２２（正極活物質層１２２ｂ）が対向する方向において、図５と同様の塩濃度の偏りが発生する。
【００２９】
二次電池１０の充放電を行うときには、負極板１２１および正極板１２２の間において、負極板１２１および正極板１２２が対向する方向（例えば、図５に示すＹ方向）に塩が移動する。二次電池１０がリチウムイオン二次電池であるとき、この塩はリチウム塩となる。負極板１２１および正極板１２２が対向する方向に塩が移動することによって、負極板１２１および正極板１２２が対向する方向において、塩濃度の偏りが発生する。
【００３０】
図６は、負極板１２１および正極板１２２のそれぞれの表面上（領域Ａ内）において、塩濃度の偏りが発生する状態を示す。図５に示す塩濃度の偏りが発生することに応じて、図６に示す塩濃度の偏りが発生する。図６では、領域Ａを含む負極板１２１の一部と、領域Ａを含む正極板１２２の一部とを上下に分けて示している。ここで、負極板１２１の一部と、正極板１２２の一部とは、セパレータ１２３を挟んで対向している。
【００３１】
図６の矢印で示すように、塩濃度の偏りは、領域Ａ内において、Ｘ方向に発生しやすい。図６では、負極板１２１および正極板１２２のそれぞれにおいて、領域Ａ内における塩濃度の分布（一例）も示している。塩濃度分布を示す図において、縦軸は塩濃度であり、横軸はＸ方向における位置である。ここで、二次電池１０の充電時では、塩濃度の偏りとして、実線で示す塩濃度分布が発生することがある。また、二次電池１０の放電時では、塩濃度の偏りとして、一点鎖線で示す塩濃度分布が発生することがある。
【００３２】
上述したように、Ｘ方向における発電要素１２０の両端部では、負極板１２１（集電箔１２１ａ）や正極板１２２（集電箔１２２ａ）がＸ軸の周りで巻かれているだけである。このため、Ｘ方向における発電要素１２０の両端部では、電解液が通過しやすい。言い換えれば、発電要素１２０の内部から発電要素１２０の外部に向かって電解液が移動したり、発電要素１２０の外部から発電要素１２０の内部に向かって電解液が移動したりしやすい。
【００３３】
これにより、図６に示すように、領域Ａ内のＸ方向において、塩濃度の偏りが発生しやすくなる。上述したように、負極板１２１、正極板１２２およびセパレータ１２３を積層しただけの構成であっても、発電要素１２０の内部から発電要素１２０の外部に向かって電解液が移動したり、発電要素１２０の外部から発電要素１２０の内部に向かって電解液が移動したりしやすい。
【００３４】
図６に示す塩濃度の偏りは、電解液の流れによって発生することが分かった。また、この電解液の流れは、電解液の体積変化（膨張および収縮）によって発生することが分かった。具体的には、電解液が膨張すると、Ｘ方向において、発電要素１２０の内部から発電要素１２０の外部に向かう電解液の流れが発生する。すなわち、図７において、Ｘ方向における領域Ａの中心Ｃを基準として、矢印Ｘ１で示す方向に電解液が流れる。一方、電解液が収縮すると、Ｘ方向において、発電要素１２０の外部から発電要素１２０の内部に向かう電解液の流れが発生する。すなわち、図７において、領域Ａの中心Ｃに向かう方向（矢印Ｘ２で示す方向）に電解液が流れる。矢印Ｘ１，Ｘ２で示す方向に電解液が流れることにより、図６に示す塩濃度の偏りが発生する。
【００３５】
そこで、本実施例では、電解液の膨張および収縮によって発生する電解液の流れ（流速）を算出することにより、図６に示す塩濃度の偏り（塩濃度分布）を把握するようにしている。図６に示す塩濃度の偏りを把握すれば、例えば、この塩濃度の偏りによって発生する抵抗上昇量Ｒｈを把握できる。本発明における塩濃度の分布とは、図６に示す塩濃度分布である。
【００３６】
以下、電解液の流れ（流速）を算出する方法について説明する。
【００３７】
下記式（１）は電解液の流れを規定する方程式であり、Brinkman-Navier-Stokes方程式として知られている。下記式（２）は電解液の流れに関する連続式であり、質量保存則から導かれる式である。下記式（１），（２）は、本発明の液流れ方程式に相当する。
【００３８】
【数１】
2016039010.tif 000003
【００３９】
上記式（１），（２）において、ｕｊは電解液の流速、ρは電解液の密度、εｅ，ｊは電解液の体積分率、ｔは時刻である。また、上記式（１）において、μは電解液の粘度、Ｋｊは透過係数、ｐは電解液の圧力である。
【００４０】
ここで、添字ｊは、負極板１２１、正極板１２２およびセパレータ１２３を区別するために用いられ、添字ｊには「ｎ」、「ｐ」および「ｓ」が含まれる。添字ｊが「ｎ」であるときには、負極板１２１に関する値を示し、添字ｊが「ｐ」であるときには、正極板１２２に関する値を示し、添字ｊが「ｓ」であるときには、セパレータ１２３に関する値を示す。
【００４１】
上述したように、電解液は、負極板１２１（負極活物質層１２１ｂ）、正極板１２２（正極活物質層１２２ｂ）およびセパレータ１２３のそれぞれに含浸している。このため、電解液に関するパラメータ（流速ｕｊ、体積分率εｅ，ｊ、透過係数Ｋｊ）としては、負極板１２１、正極板１２２およびセパレータ１２３のそれぞれにおいて規定される。なお、本明細書では、上記式（１），（２）に示すパラメータ以外にも添字ｊを用いることがある。
【００４２】
粘度μとしては予め定めた固定値を用いたり、電解液の温度に応じて粘度μを変更したりすることができる。電解液の温度としては、温度センサ３３によって検出される電池温度Ｔｂが用いられる。粘度μおよび電池温度Ｔｂの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度Ｔｂを検出することにより粘度μを特定できる。粘度μおよび電池温度Ｔｂの対応関係を示す情報はメモリ４１に記憶される。
【００４３】
密度ρは、電解液の膨張および収縮を規定するパラメータであり、電解液の膨張および収縮に応じた値を示す。電解液の膨張および収縮は電解液の温度（すなわち、電池温度Ｔｂ）に依存し、密度ρも電解液の温度（電池温度Ｔｂ）に依存する。したがって、密度ρおよび電池温度Ｔｂの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度Ｔｂを検出することにより密度ρを特定できる。密度ρを特定することにより、電解液の膨張および収縮を把握できる。密度ρおよび電池温度Ｔｂの対応関係を示す情報はメモリ４１に記憶される。体積分率εｅ，ｊや透過係数Ｋｊは予め定めた固定値を用いることができる。
【００４４】
上記式（１）では、流速ｕｊおよび圧力ｐが未知数となるため、上記式（２）に示す連続式を規定して上記式（１），（２）の連立方程式を解くことにより、流速ｕｊおよび圧力ｐを算出できる。上記式（１），（２）には、電解液の膨張および収縮を規定する密度ρが含まれているため、上記式（１），（２）の連立方程式を解くことにより、電解液の膨張および収縮に応じた流速ｕｊを算出できる。この流速ｕｊには、図７の矢印Ｘ１で示す方向の流速と、図７の矢印Ｘ２で示す方向の流速とが含まれる。流速ｕｊを算出するときには、例えば、上記式（１），（２）を用いて収束計算を行うことができる。また、上記式（１），（２）を用いた演算は所定の周期で行われるが、前回の演算周期で算出された値を今回の演算周期で用いることにより、流速ｕｊを算出できる。
【００４５】
負極板１２１では、負極活物質層１２１ｂの内部において電解液が移動する。このため、流速ｕｊ（すなわち、流速ｕｎ）は、負極活物質層１２１ｂの内部における位置毎に算出される。図６に示す塩濃度の偏りを把握するときには、Ｘ方向で互いに異なる位置において、流速ｕｎが算出される。
【００４６】
正極板１２２では、正極活物質層１２２ｂの内部において電解液が移動する。このため、流速ｕｊ（すなわち、流速ｕｐ）は、正極活物質層１２２ｂの内部における位置毎に算出される。図６に示す塩濃度の偏りを把握するときには、Ｘ方向で互いに異なる位置において、流速ｕｐが算出される。セパレータ１２３では、この内部において電解液が移動する。このため、流速ｕｊ（すなわち、流速ｕｓ）は、セパレータ１２３の内部における位置毎に算出される。図６に示す塩濃度の偏りを把握するときには、Ｘ方向で互いに異なる位置において、流速ｕｓが算出される。
【００４７】
一方、上記式（１），（２）において、様々な条件を仮定した上で、上記式（１）を簡易化することもできる。以下、上記式（１）を簡易化するときの手法（一例）について説明する。
【００４８】
電解液の密度ρが流速ｕｊを算出する位置（Ｘ方向の位置を含む）に関わらず一定であると仮定すると、上記式（２）は下記式（３）で表される。
【００４９】
【数２】
2016039010.tif 000004
【００５０】
上記式（３）から下記式（４）が導き出せる。下記式（４）に示すｖｐ，ｖｎ，ｖｓは、正極板１２２、負極板１２１およびセパレータ１２３のそれぞれにおける電解液の動粘性係数（ｖ＝μ／ρ）である。下記式（４）に示すｘ，ｙは、Ｘ方向およびＹ方向における位置をそれぞれ示す。ここで、Ｙ方向とは、負極板１２１および正極板１２２がセパレータ１２３を挟んで対向する方向（図７の上下方向）である。
【００５１】
【数３】
2016039010.tif 000005
【００５２】
電解液の連続性を考慮すると、上記式（４）は下記式（５）で表される。
【００５３】
【数４】
2016039010.tif 000006
【００５４】
一方、上記式（１）において、下記式（６）に示す仮定を行うと、下記式（７）が得られる。
【００５５】
【数５】
2016039010.tif 000007
【００５６】
Ｘ方向における電解液の圧力分布に関して、負極板１２１、正極板１２２およびセパレータ１２３における圧力分布が互いに等しいと仮定すると、下記式（８）が得られる。
【００５７】
【数６】
2016039010.tif 000008
【００５８】
上記式（５），（８）によれば、下記式（９）が導き出せる。下記式（９）は、負極板１２１において、Ｘ方向の位置に応じた流速ｕｎを示している。下記式（９）に示すｘはＸ方向の位置を示している。
【００５９】
【数７】
2016039010.tif 000009
【００６０】
上記式（９）をｘで積分し、ｘが０であるときの流速ｕｎを０と仮定すると、下記式（１０）が得られる。ここで、ｘが０であるときの位置は、Ｘ方向における領域Ａの一端を示している。また、Ｘ方向における領域Ａの長さをＬとすると、Ｘ方向における領域Ａの他端については、ｘがＬになる。
【００６１】
【数８】
2016039010.tif 000010
【００６２】
流速ｕｎと同様に、流速ｕｐ，ｕｓは下記式（１１），（１２）で表される。
【００６３】
【数９】
2016039010.tif 000011
【００６４】
上記式（１０）〜（１２）は上記式（１）を簡易化した式となり、負極板１２１、正極板１２２およびセパレータ１２３のそれぞれにおいて、Ｘ方向の位置に応じた流速ｕｎ，ｕｐ，ｕｓを算出できる。上記式（１０）〜（１２）は、本発明における液流れ方程式に相当する。なお、電解液の流速ｕｊを算出するための方程式は、上記（１），（２），（１０）〜（１２）に限るものではない。電解液の流れは電解液の体積変化によって発生するため、電解液の体積変化を規定する電解液の密度をパラメータ（変数）として含み、電解液の流れを規定できる方程式であれば、本発明を適用できる。
【００６５】
流速ｕｊを算出すれば、下記式（１３）に基づいて、電解液中の塩濃度ｃｅ，ｊを算出できる。ここで、塩濃度ｃｅ，ｊを算出するときには、流速ｕｎ，ｕｐ，ｕｓのすべてを考慮してもよいし、流速ｕｎ，ｕｐ，ｕｓの一部（例えば、流速ｕｎ）だけを考慮してもよい。
【００６６】
【数１０】
2016039010.tif 000012
【００６７】
上記式（１３）において、Ｄｅ，ｊｅｆｆは、電解液の実効拡散係数であり、ｔ＋０は電解液中の塩の輸率である。Ｆはファラデー定数であり、ｊｊは、単位体積および単位時間において、電解液中の塩の生成量である。
【００６８】
上記式（１３）の左辺第１項は、所定時間Δｔにおける塩濃度の変化を規定している。上記式（１３）の左辺第２項は、電解液の流れ（流速ｕｊ）に依存する塩濃度の変化を規定している。上記式（１３）の右辺第１項は、電解液中の塩の拡散状態を規定している。上記式（１３）の右辺第２項は、塩の生成量を規定している。ここで、二次電池１０の放電時では、負極板１２１の表面（負極活物質層１２１ｂ）において塩が生成され、二次電池１０の充電時では、正極板１２２の表面（正極活物質層１２２ｂ）において塩が生成される。
【００６９】
上記式（１３）によれば、図５および図６に示す塩濃度の偏りを把握することができる。上記式（１３）を解くことにより塩濃度ｃｅ，ｊを算出できる。ここで、流速ｕｊとしては、Ｘ方向の位置に応じた流速ｕｊが用いられるため、上記式（１３）を解くことにより、Ｘ方向の位置に応じた塩濃度ｃｅ，ｊを算出できる。これにより、Ｘ方向における塩濃度ｃｅ，ｊの分布（図６参照）を算出できる。
【００７０】
なお、塩濃度ｃｅ，ｊの算出方法は、上記式（１３）を用いた方法に限るものではない。Ｘ方向において塩濃度ｃｅ，ｊのバラツキが発生している状態であれば、流速ｕｊに基づいて、Ｘ方向における塩濃度ｃｅ，ｊの分布を算出できる。上記式（１３）では、塩濃度ｃｅ，ｊのバラツキを把握するために、塩の拡散状態および塩の生成量を規定している。例えば、塩濃度ｃｅ，ｊのバラツキを予め設定すれば、塩の拡散状態や生成量を考慮せずに、流速ｕｊに基づいて塩濃度ｃｅ，ｊの分布を算出できる。
【００７１】
また、上記式（１３）によれば、電解液中の塩の拡散状態を規定しているが、塩の拡散状態を考慮しなくてもよい。塩の拡散に関する時定数は、塩の生成に関する時定数よりも大きいため、塩の拡散が発生していないこともある。そこで、塩濃度ｃｅ，ｊの分布を算出するときに、塩の拡散状態を考慮しないようにしてもよい。
【００７２】
上述したように塩濃度ｃｅ，ｊの分布を算出すれば、塩濃度ｃｅ，ｊの最大差（塩濃度差）Δｃｅ，ｊ＿ｍａｘを算出できる。塩濃度差Δｃｅ，ｊ＿ｍａｘは、塩濃度（最大値）ｃｅおよび塩濃度（最小値）ｃｅの差である。
【００７３】
図８に示すように、抵抗上昇量Ｒｈおよび塩濃度差Δｃｅ＿ｍａｘの対応関係を実験などによって予め求めておけば、塩濃度差Δｃｅ＿ｍａｘを算出することにより、この塩濃度差Δｃｅ＿ｍａｘに対応した抵抗上昇量Ｒｈを算出できる。図８に示すように、塩濃度差Δｃｅ＿ｍａｘが大きくなるほど、抵抗上昇量Ｒｈが大きくなる。言い換えれば、塩濃度差Δｃｅ＿ｍａｘが小さくなるほど、抵抗上昇量Ｒｈが小さくなる。上記式（１３）によれば、塩濃度ｃｅは、負極板１２１および正極板１２２のそれぞれで算出される。ここで、図８に示す塩濃度差Δｃｅ＿ｍａｘを算出するときには、まず、負極板１２１および正極板１２２における塩濃度ｃｅの分布を合算する。具体的には、負極板１２１および正極板１２２において、互いに対向する位置における塩濃度ｃｅを合算する。そして、合算した塩濃度ｃｅの分布において、塩濃度（最大値）ｃｅおよび塩濃度（最小値）ｃｅの差を塩濃度差Δｃｅ＿ｍａｘとして算出する。
【００７４】
抵抗上昇量Ｒｈおよび塩濃度差Δｃｅ＿ｍａｘの対応関係は、マップ又は演算式として表すことができる。そして、この対応関係を特定する情報はメモリ４１に記憶しておくことができる。上記式（１３）によれば、所定時間（演算周期）Δｔが経過するたびに塩濃度ｃｅ，ｊが算出されて塩濃度ｃｅ，ｊの分布を把握できるため、所定時間Δｔが経過するたびに抵抗上昇量Ｒｈが算出される。これに伴い、抵抗上昇量Ｒｈの変化を把握できる。
【００７５】
一方、抵抗上昇量Ｒｈは、平均塩濃度ｃｅ＿ａｖｅに依存することがある。そこで、平均塩濃度ｃｅ＿ａｖｅに基づいて、抵抗上昇量Ｒｈを算出することもできる。平均塩濃度ｃｅ＿ａｖｅとは、負極板１２１および正極板１２２における塩濃度（上述した合算値）ｃｅの分布を平均化した塩濃度（平均値）である。塩濃度差Δｃｅ＿ｍａｘと同様に、抵抗上昇量Ｒｈおよび平均塩濃度ｃｅ＿ａｖｅの対応関係（マップ又は演算式）を実験などによって予め求めておけば、平均塩濃度ｃｅ＿ａｖｅを算出することにより、この平均塩濃度ｃｅ＿ａｖｅに対応した抵抗上昇量Ｒｈを算出できる。
【００７６】
また、塩濃度差Δｃｅ＿ｍａｘおよび平均塩濃度ｃｅ＿ａｖｅに基づいて、抵抗上昇量Ｒｈを算出することもできる。この場合には、塩濃度差Δｃｅ＿ｍａｘ、平均塩濃度ｃｅ＿ａｖｅおよび抵抗上昇量Ｒｈの対応関係（マップ又は演算式）を予め求めておけばよい。塩濃度差Δｃｅ＿ｍａｘおよび平均塩濃度ｃｅ＿ａｖｅを算出すれば、算出した塩濃度差Δｃｅ＿ｍａｘおよび平均塩濃度ｃｅ＿ａｖｅに対応する抵抗上昇量Ｒｈを算出できる。
【００７７】
図９は、塩濃度分布を算出する処理を示すフローチャートである。図９に示す処理は、コントローラ４０によって実行される。
【００７８】
ステップＳ１０１において、コントローラ４０は、温度センサ３３を用いて電池温度Ｔｂを検出する。ステップＳ１０２において、コントローラ４０は、電解液の流速ｕｊを算出する。具体的には、上述したように、コントローラ４０は、Ｘ方向の位置に応じた流速ｕｊを算出する。流速ｕｊを算出するときには、ステップＳ１０１の処理で検出された電池温度Ｔｂが用いられる。
【００７９】
ステップＳ１０３において、コントローラ４０は、ステップＳ１０２の処理で算出された流速ｕｊに基づいて塩濃度ｃｅ，ｊの分布を算出する。塩濃度ｃｅ，ｊの分布を算出する方法は上述した通りである。
【００８０】
本実施例では、ステップＳ１０１の処理において、温度センサ３３を用いて電池温度Ｔｂを検出しているが、これに限るものではない。具体的には、ステップＳ１０１の処理において、二次電池１０を充放電したときの発熱量（電池温度Ｔｂの上昇量）と、二次電池１０の放熱量（電池温度Ｔｂの低下量）とを考慮することにより、電池温度Ｔｂを推定できる。
【００８１】
二次電池１０の発熱量は、二次電池１０の電流値Ｉｂおよび内部抵抗値から算出できる。二次電池１０の内部抵抗値は、電流値Ｉｂおよび電圧値Ｖｂから算出できる。二次電池１０の放熱量は、電池温度Ｔｂと、二次電池１０の周囲に存在する大気中の温度（環境温度）とに基づいて算出できる。ここで、電池温度Ｔｂが環境温度よりも高いときにおいて、電池温度Ｔｂおよび環境温度の差が広がるほど、放熱量が大きくなりやすい。
【００８２】
二次電池１０を充放電していなければ、電池温度Ｔｂは環境温度と等しくなる。そこで、温度センサ３３とは異なる温度センサを用いて環境温度を検出しておき、二次電池１０を充放電するたびに算出される発熱量と、二次電池１０の放熱量とに基づいて、現在の電池温度Ｔｂを算出（推定）できる。発熱量および放熱量に基づいて電池温度Ｔｂを推定する方法は、公知の方法を適宜採用できる。二次電池１０を車両に搭載したとき、車両には、車両の外部における温度を検出する温度センサが設けられている。この温度センサを用いて、環境温度を検出できる。
【００８３】
図９に示す処理によって塩濃度分布を算出したときには、この塩濃度分布に基づいて抵抗上昇量Ｒｈを算出し、この抵抗上昇量Ｒｈに基づいて二次電池１０の充放電を制御できる。ここで、二次電池１０の充放電を制御するときの処理（一例）について、図１０に示すフローチャートを用いて説明する。図１０に示す処理は、コントローラ４０によって実行される。
【００８４】
ステップＳ２０１において、コントローラ４０は、図９に示す処理で算出した塩濃度ｃｅ，ｊの分布に基づいて抵抗上昇量Ｒｈを算出する。抵抗上昇量Ｒｈを算出する方法は上述した通りである。ステップＳ２０２において、コントローラ４０は、抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるか否かを判別する。閾値Ｒｈ＿ｔｈは、抵抗上昇量Ｒｈの上限値であり、二次電池１０の劣化を抑制する観点に基づいて適宜設定できる。閾値Ｒｈ＿ｔｈを特定する情報は、メモリ４１に記憶しておくことができる。
【００８５】
抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるとき、コントローラ４０は、ステップＳ２０３において、抵抗上昇量Ｒｈおよび閾値Ｒｈ＿ｔｈに基づいて、抵抗上昇抑制量ΔＲｈを算出する。ここで、抵抗上昇抑制量ΔＲｈは、ステップＳ２０１の処理で算出された抵抗上昇量Ｒｈから閾値Ｒｈ＿ｔｈを減算した値、すなわち、抵抗上昇量Ｒｈおよび閾値Ｒｈ＿ｔｈの差である。
【００８６】
ステップＳ２０４において、コントローラ４０は、ステップＳ２０３の処理で算出した抵抗上昇抑制量ΔＲｈと、正負極間の塩濃度差Δｃｅ，ｎｐとに基づいて、流速制限値ｕ＿ｌｉｍを算出する。以下、流速制限値ｕ＿ｌｉｍを算出する方法について説明する。
【００８７】
正負極間の塩濃度差Δｃｅ，ｎｐとは、図５に示す塩濃度分布において、正極板１２２における塩濃度ｃｅと、負極板１２１における塩濃度ｃｅとの差である。上述したように、塩濃度差Δｃｅ＿ｍａｘを算出するときに、正極板１２２および負極板１２１における塩濃度ｃｅがそれぞれ算出されているため、これらの塩濃度ｃｅの差（すなわち、正負極間の塩濃度差Δｃｅ，ｎｐ）を算出することができる。
【００８８】
抵抗上昇量Ｒｈは、塩濃度差Δｃｅ，ｎｐおよび流速ｕに依存する。具体的には、塩濃度差Δｃｅ，ｎｐが大きいほど、また、流速ｕが高いほど、抵抗上昇量Ｒｈが大きくなる。言い換えれば、塩濃度差Δｃｅ，ｎｐが小さいほど、また、流速ｕが低いほど、抵抗上昇量Ｒｈが小さくなる。
【００８９】
図１１に示すように、正負極間の塩濃度差Δｃｅ，ｎｐが所定値であるとき、流速ｕを低下させれば、抵抗上昇量Ｒｈを低下させることができる。図１１において、縦軸は抵抗上昇量Ｒｈを示し、横軸は流速ｕを示す。抵抗上昇量Ｒｈ＿ｃは、ステップＳ２０１の処理で算出された抵抗上昇量Ｒｈを示し、流速ｕ＿ｃは、抵抗上昇量Ｒｈ＿ｃが発生するときの流速ｕを示す。すなわち、流速ｕ＿ｃは、図９に示すステップＳ１０２の処理で算出された流速ｕｊとなる。
【００９０】
図１１に示す抵抗上昇量Ｒｈ＿ｌｉｍは、流速制限値ｕ＿ｌｉｍに対応した抵抗上昇量Ｒｈである。すなわち、抵抗上昇量Ｒｈ＿ｌｉｍは、塩濃度差Δｃｅ，ｎｐが所定値であり、流速ｕが流速制限値ｕ＿ｌｉｍであるときの抵抗上昇量Ｒｈとなる。ここで、抵抗上昇量Ｒｈ＿ｃおよび抵抗上昇量Ｒｈ＿ｌｉｍの差が、抵抗上昇抑制量ΔＲｈとなる。
【００９１】
上述したように、流速ｕ＿ｃおよび抵抗上昇量Ｒｈ＿ｃは予め算出されているため、ステップＳ２０３の処理で算出された抵抗上昇抑制量ΔＲｈが得られるときの流速制限値ｕ＿ｌｉｍを算出できる。すなわち、抵抗上昇量Ｒｈ＿ｃから抵抗上昇抑制量ΔＲｈを減算すれば、抵抗上昇量Ｒｈ＿ｌｉｍを算出できる。塩濃度差Δｃｅ，ｎｐが所定値であるとき、抵抗上昇量Ｒｈが抵抗上昇量Ｒｈ＿ｌｉｍとなるときの流速ｕを算出できる。このときの流速ｕが流速制限値ｕ＿ｌｉｍとなる。
【００９２】
ステップＳ２０５において、コントローラ４０は、ステップＳ２０４の処理で算出した流速制限値ｕ＿ｌｉｍに基づいて、電流制限値Ｉｂ＿ｍａｘ，Ｉｂ＿ｍｉｎを算出する。具体的には、図１２に示すように、流速ｕ（Ｔｂ）および電流値Ｉｂの対応関係を予め求めておく。この対応関係は、マップ又は演算式として表され、この対応関係の情報はメモリ４１に記憶しておくことができる。
【００９３】
図１２において、縦軸は流速ｕを示し、横軸は電流値Ｉｂを示す。ここで、流速ｕが０よりも高いときには、図７に示す矢印Ｘ１の方向において電解液が流れているときの流速を示す。また、流速ｕが０よりも低いときには、図７に示す矢印Ｘ２の方向において電解液が流れているときの流速を示す。流速ｕが０であれば、図７に示す矢印Ｘ１，Ｘ２のいずれの方向においても電解液が流れていない。
【００９４】
図１２に示す流速ｕ（Ｔｂ）は、各電流値Ｉｂにおいて、電解液の膨張や収縮に応じた流速ｕである。電流値（絶対値）Ｉｂが大きいほど、発電要素１２０の発熱量が増加し、電解液が膨張しやすくなる。この電解液の膨張により、図７に示す矢印Ｘ１の方向に電解液が流れやすくなるため、電流値（絶対値）Ｉｂが大きいほど、流速ｕが高くなる。言い換えれば、電流値（絶対値）Ｉｂが小さいほど、流速ｕが低くなる。したがって、流速ｕ（Ｔｂ）および電流値Ｉｂの関係は、図１２に示す関係となる。電解液の膨張や収縮は、電池温度Ｔｂにも依存するため、図１２に示す関係は、電池温度Ｔｂ毎に用意される。
【００９５】
コントローラ４０は、電池温度Ｔｂを検出すれば、この電池温度Ｔｂに対応した図１２に示す対応関係を特定できる。このように特定された図１２に示す対応関係において、ステップＳ２０４の処理で算出された流速制限値ｕ＿ｌｉｍに対応する電流値Ｉｂが電流制限値Ｉｂ＿ｍａｘ，Ｉｂ＿ｍｉｎとなる。電流制限値（正の値）Ｉｂ＿ｍａｘは、二次電池１０を放電するときの電流制限値となり、電流制限値（負の値）Ｉｂ＿ｍｉｎは、二次電池１０を充電するときの電流制限値となる。図１２から分かるように、流速制限値ｕ＿ｌｉｍが低いほど、電流制限値Ｉｂ＿ｍａｘ，Ｉｂ＿ｍｉｎが０［Ａ］に近づく。
【００９６】
ステップＳ２０６において、コントローラ４０は、充電電力許容値Ｗｉｎや放電電力許容値Ｗｏｕｔを設定する。充電電力許容値Ｗｉｎは、二次電池１０の充電を許容する上限の電力値であり、放電電力許容値Ｗｏｕｔは、二次電池１０の放電を許容する上限の電力値である。
【００９７】
上述したように、二次電池１０を充電したときの電流値Ｉｂは負の値となるため、充電電力値も負の値となる。一方、二次電池１０を放電したときの電流値Ｉｂは正の値となるため、放電電力値は正の値となる。二次電池１０を充電するときには、充電電力値（絶対値）が充電電力許容値（絶対値）Ｗｉｎよりも高くならないように充電が制御される。また、二次電池１０を放電するときには、放電電力値が放電電力許容値Ｗｏｕｔよりも高くならないように放電が制御される。
【００９８】
コントローラ４０は、ステップＳ２０５の処理によって電流制限値Ｉｂ＿ｍａｘを算出すれば、この電流制限値Ｉｂ＿ｍａｘに対応した放電電力許容値Ｗｏｕｔを算出できる。ステップＳ２０６の処理では、このように算出した放電電力許容値Ｗｏｕｔが設定される。電流制限値Ｉｂ＿ｍａｘが０［Ａ］に近づくほど、放電電力許容値Ｗｏｕｔが低下し、二次電池１０の放電が制限されやすくなる。
【００９９】
また、コントローラ４０は、ステップＳ２０５の処理によって電流制限値Ｉｂ＿ｍｉｎを算出すれば、この電流制限値Ｉｂ＿ｍｉｎに対応した充電電力許容値Ｗｉｎを算出できる。ステップＳ２０６の処理では、このように算出した充電電力許容値Ｗｉｎが設定される。電流制限値Ｉｂ＿ｍｉｎが０［Ａ］に近づくほど、充電電力許容値（絶対値）Ｗｉｎが低下し、二次電池１０の充電が制限されやすくなる。
【０１００】
抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈよりも小さいとき、コントローラ４０は、図１０に示す処理を終了する。このとき、充電電力許容値Ｗｉｎとしては、基準値としての充電電力許容値Ｗｉｎ＿ｒｅｆが設定され、放電電力許容値Ｗｏｕｔとしては、基準値としての放電電力許容値Ｗｏｕｔ＿ｒｅｆが設定される。
【０１０１】
充電電力許容値Ｗｉｎ＿ｒｅｆや放電電力許容値Ｗｏｕｔ＿ｒｅｆは、電池温度Ｔｂや二次電池１０のＳＯＣに基づいて設定される。ステップＳ２０６の処理で設定される放電電力許容値Ｗｏｕｔは、放電電力許容値Ｗｏｕｔ＿ｒｅｆよりも低くなる。また、ステップＳ２０６の処理で設定される充電電力許容値（絶対値）Ｗｉｎ＿ｒｅｆは、充電電力許容値（絶対値）Ｗｉｎ＿ｒｅｆよりも低くなる。
【０１０２】
図１０に示す処理では、抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるとき、充電電力許容値（絶対値）Ｗｉｎおよび放電電力許容値Ｗｏｕｔを低下させているが、これに限るものではない。具体的には、充電電力許容値（絶対値）Ｗｉｎおよび放電電力許容値Ｗｏｕｔの一方だけを低下させることができる。
【０１０３】
抵抗上昇量Ｒｈとしては、二次電池１０の充電に起因した抵抗上昇量Ｒｈと、二次電池１０の放電に起因した抵抗上昇量Ｒｈとがある。図６に示すように、充電時および放電時において、塩濃度ｃｅ，ｊの分布が異なる。したがって、塩濃度ｃｅ，ｊの分布を把握すれば、二次電池１０の充電に起因した抵抗上昇量Ｒｈと、二次電池１０の放電に起因した抵抗上昇量Ｒｈとを区別できる。
【０１０４】
ここで、二次電池１０の充電によって塩濃度差Δｃｅ＿ｍａｘが発生し、この塩濃度差Δｃｅ＿ｍａｘから算出された抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるときには、充電電力許容値（絶対値）Ｗｉｎだけを低下させることができる。これにより、充電に起因した抵抗上昇量Ｒｈに関して、抵抗上昇量Ｒｈの増加を抑制したり、抵抗上昇量Ｒｈを減少させたりすることができる。
【０１０５】
一方、二次電池１０の放電によって塩濃度差Δｃｅ＿ｍａｘが発生し、この塩濃度差Δｃｅ＿ｍａｘから算出された抵抗上昇量Ｒｈが閾値Ｒｈ＿ｔｈ以上であるときには、放電電力許容値Ｗｏｕｔだけを低下させることができる。これにより、放電に起因した抵抗上昇量Ｒｈに関して、抵抗上昇量Ｒｈの増加を抑制したり、抵抗上昇量Ｒｈを減少させたりすることができる。
【実施例２】
【０１０６】
本発明の実施例２について説明する。本実施例において、実施例１で説明した構成要素と同じ構成要素については同一の符号を用い、詳細な説明は省略する。以下、実施例１と異なる点について、主に説明する。
【０１０７】
実施例１では、電解液の膨張および収縮によって発生する電解液の流れ（流速ｕｊ）を考慮して塩濃度分布を算出している。ここで、電解液の流れは、電解液の膨張および収縮に加えて、活物質（正極活物質や負極活物質）の膨張および収縮によっても発生することがある。本実施例では、電解液の膨張および収縮と、活物質の膨張および収縮とによって発生する電解液の流れ（流速ｕｊ）を考慮して塩濃度分布を算出している。
【０１０８】
活物質の膨張時には活物質の体積が増加し、活物質の収縮時には活物質の体積が減少する。活物質の体積は、電池温度Ｔｂや二次電池１０のＳＯＣに依存する。すなわち、電池温度ＴｂやＳＯＣが変化すれば、活物質の体積が変化する。
【０１０９】
負極活物質層１２１ｂや正極活物質層１２２ｂの体積が変化しないと仮定すると、負極活物質や正極活物質の体積が変化することにより、各活物質層１２１ｂ，１２２ｂに含浸された電解液の体積分率εｅ，ｊが変化する。したがって、電解液の流れ（流速ｕｊ）を算出するときには、体積分率εｅ，ｊの変化を考慮することが好ましい。体積分率εｅ，ｊの変化量Δεｅ，ｊは、下記式（１４）に基づいて算出できる。
【０１１０】
【数１１】
2016039010.tif 000013
【０１１１】
上記式（１４）において、βｊ（具体的にはβｎ又はβｐ）は活物質（負極活物質又は正極活物質）の体積膨張率であり、εｅ，ｊは前回の演算周期で算出された体積分率εｅ，ｊである。体積膨張率βｊは、活物質の膨張および収縮を規定する値であるため、上記式（１４）によれば、体積膨張率βｊに基づいて変化量Δεｅ，ｊを算出できる。今回の演算周期において変化量Δεｅ，ｊを算出したときには、前回の体積分率εｅ，ｊに対して今回の変化量Δεｅ，ｊを加算することにより、今回の演算周期における体積分率εｅ，ｊを算出できる。
【０１１２】
ここで、変化量Δεｅ，ｊを最初に算出するとき、上記式（１４）に示す体積分率εｅ，ｊとしては、電池温度Ｔｂや二次電池１０のＳＯＣに応じた値を用いることができる。電池温度ＴｂおよびＳＯＣの少なくとも一方と、体積分率εｅ，ｊとの対応関係を示す情報を実験などによって予め求めておけば、電池温度ＴｂやＳＯＣを特定することにより、体積分率εｅ，ｊを特定できる。
【０１１３】
体積膨張率βｊは線膨張率αｊから算出できる。すなわち、線膨張率αｊの３倍の値が体積膨張率βｊとなる。ここで、線膨張率αｊは、電池温度Ｔｂや二次電池１０のＳＯＣに依存するため、電池温度ＴｂおよびＳＯＣの少なくとも一方から線膨張率αｊを特定できる。例えば、線膨張率αｊおよび電池温度Ｔｂの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度Ｔｂの検出又は推定によって線膨張率αｊを特定できる。また、線膨張率αｊおよびＳＯＣの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、ＳＯＣの推定によって線膨張率αｊを特定できる。上述した対応関係を示す情報はメモリ４１に記憶しておくことができる。なお、線膨張率αｊを算出せずに、電池温度ＴｂおよびＳＯＣの少なくとも一方から体積膨張率βｊを直接算出することもできる。
【０１１４】
なお、各活物質層１２１ｂ，１２２ｂに関する変化量Δεｅ，ｎ，Δεｅ，ｐに加えて、セパレータ１２３に関する変化量Δεｅ，ｓを算出してもよい。変化量Δεｅ，ｓを算出するとき、上記式（１４）に示すβｊ（すなわち、βｓ）は、セパレータ１２３を構成する材料の体積膨張率となる。体積膨張率βｓは電池温度Ｔｂに依存する。このため、体積膨張率βｓおよび電池温度Ｔｂの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度Ｔｂの検出又は推定によって体積膨張率βｓを特定できる。セパレータ１２３に関して、変化量Δεｅ，ｓおよび体積分率εｅ，ｓを算出する方法は、上述した算出方法と同様である。
【０１１５】
上記式（１４）の変化量Δεｅ，ｊから算出された体積分率εｅ，ｊは、上記式（１），（２），（１０）〜（１３）に示す体積分率εｅ，ｊとして用いられる。そして、実施例１と同様に、電解液の流速ｕｊを算出して塩濃度分布を算出することにより、二次電池１０の抵抗上昇量Ｒｈを算出できる。本実施例において、流速ｕｊを算出するときには、体積膨張率βｊを算出するために、図９に示すステップＳ１０１の処理において、電池温度ＴｂおよびＳＯＣの少なくとも一方が特定される。
【０１１６】
本実施例によれば、活物質の膨張および収縮によって電解液の流れが発生するときには、体積分率εｅ，ｊの変化を考慮して流速ｕｊを算出することにより、流速ｕｊの算出精度を向上させることができる。これに伴い、塩濃度分布の算出精度を向上させたり、抵抗上昇量Ｒｈの算出精度を向上させたりすることができる。
【０１１７】
本実施例において、二次電池１０の充放電を制御するときには、図１０に示す処理と同様の処理を行うことができる。ただし、ステップＳ２０５の処理において、電流制限値Ｉｂ＿ｍａｘ，Ｉｂ＿ｍｉｎを算出するときには、電解液の膨張や収縮だけでなく、活物質の膨張や収縮も考慮している。この点について、図１３を用いて説明する。
【０１１８】
図１３は、図１２に対応した図である。本実施例において、電流制限値Ｉｂ＿ｍａｘ，Ｉｂ＿ｍｉｎを算出するときには、図１３の実線で示すように、流速ｕ（Ｔｂ，ＳＯＣ）および電流値Ｉｂの対応関係が用いられる。図１３の実線で示す対応関係は、図１３の一点鎖線で示す対応関係と、図１３の二点差線で示す対応関係とから特定される。
【０１１９】
図１３の一点鎖線で示す対応関係は、図１２に示す対応関係となる。図１３の二点差線で示す対応関係は、流速ｕ（ＳＯＣ）および電流値Ｉｂの対応関係である。この対応関係は、マップ又は演算式として表され、この対応関係の情報はメモリ４１に記憶しておくことができる。
【０１２０】
図１３に示す流速ｕ（ＳＯＣ）は、各電流値Ｉｂにおいて、活物質（負極活物質）の膨張や収縮に応じた流速ｕである。充電時の電流値（絶対値）Ｉｂが大きいほど、負極活物質が膨張して、図７に示す矢印Ｘ１の方向に電解液が流れやすくなる。このため、充電時の電流値（絶対値）Ｉｂが大きいほど、矢印Ｘ１の方向における流速ｕが高くなる。言い換えれば、充電時の電流値（絶対値）Ｉｂが小さいほど、矢印Ｘ１の方向における流速ｕが低くなる。
【０１２１】
一方、放電時の電流値Ｉｂが大きいほど、負極活物質が収縮して、図７に示す矢印Ｘ２の方向に電解液が流れやすくなる。このため、放電時の電流値Ｉｂが大きいほど、矢印Ｘ２の方向における流速ｕが高くなる。言い換えれば、放電時の電流値Ｉｂが小さいほど、矢印Ｘ２の方向における流速ｕが低くなる。したがって、流速ｕ（ＳＯＣ）および電流値Ｉｂの関係は、図１３の二点差線で示す関係となる。負極活物質の膨張や収縮は、二次電池１０のＳＯＣにも依存するため、図１３の二点差線で示す対応関係は、二次電池１０のＳＯＣ毎に用意される。
【０１２２】
コントローラ４０は、二次電池１０のＳＯＣを算出すれば、このＳＯＣに対応した、図１３の二点差線で示す対応関係を特定できる。図１３において、一点鎖線および二点差線で示す対応関係を特定すれば、実線で示す対応関係を算出することができる。実線で示す流速ｕ（Ｔｂ，ＳＯＣ）は、各電流値Ｉｂにおいて、一点鎖線で示す流速ｕ（Ｔｂ）に二点差線で示す流速ｕ（ＳＯＣ）を加算した値である。ここで、流速ｕ（ＳＯＣ）が負の値であるとき、流速ｕ（Ｔｂ）から流速ｕ（ＳＯＣ）を減算した値が流速ｕ（Ｔｂ，ＳＯＣ）となる。
【０１２３】
上述したように、コントローラ４０は、電池温度Ｔｂを検出し、二次電池１０のＳＯＣを算出すれば、図１３の実線で示す対応関係を特定できる。図１３の実線で示す対応関係において、図１０のステップＳ２０４の処理で算出された流速制限値ｕ＿ｌｉｍに対応する電流値Ｉｂが電流制限値Ｉｂ＿ｍａｘ，Ｉｂ＿ｍｉｎとなる。
【実施例３】
【０１２４】
本発明の実施例３について説明する。本実施例において、実施例１，２で説明した構成要素と同じ構成要素については同一の符号を用い、詳細な説明は省略する。以下、実施例１，２と異なる点について、主に説明する。
【０１２５】
実施例２では、負極活物質層１２１ｂや正極活物質層１２２ｂの体積が変化しないと仮定した上で、体積分率εｅ，ｊが変化するものとしている。一方、活物質の膨張や収縮によって各活物質層１２１ｂ，１２２ｂの体積が変化することがある。この場合には、各活物質層１２１ｂ，１２２ｂの体積変化を考慮する必要がある。
【０１２６】
そこで、本実施例では、各活物質層１２１ｂ，１２２ｂの体積変化を考慮して流速ｕｊを算出している。各活物質層１２１ｂ，１２２ｂの体積が変化すると、各活物質層１２１ｂ，１２２ｂにおいて、電解液が存在する空間の体積が変化することに応じて体積分率εｅ，ｊが変化する。この点を考慮して、体積分率εｅ，ｊを算出できる。以下、体積分率εｅ，ｊの算出方法について、具体的に説明する。
【０１２７】
本実施例では、各活物質層１２１ｂ，１２２ｂの体積変化を各活物質層１２１ｂ，１２２ｂの厚さの変化とみなしている。各活物質層１２１ｂ，１２２ｂの厚さとは、負極板１２１および正極板１２２が対向する方向（図７の上下方向）における各活物質層１２１ｂ，１２２ｂの長さである。
【０１２８】
各活物質層１２１ｂ，１２２ｂの厚さが変化した分だけ、発電要素１２０の厚さが変化するときには、発電要素１２０の厚さの変化量が下記式（１５）によって表される。発電要素１２０の厚さとは、図４において、Ｙ方向における発電要素１２０の長さである。
【０１２９】
【数１２】
2016039010.tif 000014
【０１３０】
上記式（１５）において、Δｙ０は発電要素１２０の厚さの変化量である。変化量Δｙ０は、正極活物質層１２２ｂの厚さの変化量Δｙｐと、負極活物質層１２１ｂの厚さの変化量Δｙｎと、セパレータ１２３の厚さの変化量Δｙｓとの総和になる。上記式（１５）では、セパレータ１２３の厚さの変化量Δｙｓも考慮している。セパレータ１２３の厚さとは、負極板１２１および正極板１２２が対向する方向（図７の上下方向）におけるセパレータ１２３の長さである。なお、変化量Δｙｓを無視することもできる。
【０１３１】
変化量Δｙｐは、正極活物質層１２２ｂの厚さ（変化前の厚さ）ｙｐに正極活物質の体積膨張率βｐを乗算することによって求められる。変化量Δｙｎは、負極活物質層１２１ｂの厚さ（変化前の厚さ）ｙｎに負極活物質の体積膨張率βｎを乗算することによって求められる。変化量Δｙｓは、セパレータ１２３の厚さ（変化前の厚さ）ｙｓにセパレータ１２３を構成する材料（すなわち、セパレータ１２３自体）の体積膨張率βｓを乗算することによって求められる。
【０１３２】
二次電池１０には拘束力が与えられることがある。この拘束力は、Ｙ方向において二次電池１０を挟む力である。二次電池１０を固定したり、発電要素１２０の入出力性能を確保したりするために、二次電池１０に拘束力を与えることがある。二次電池１０（発電要素１２０）に拘束力が与えられており、発電要素１２０の厚さの変化が制限されるとき、発電要素１２０の厚さの変化量Δｙは変化量Δｙ０よりも小さくなる。すなわち、変化量Δｙは、発電要素１２０の厚さが変化するときの力と、発電要素１２０に与えられた拘束力とが釣り合っているときの変化量となる。
【０１３３】
ここで、二次電池１０の厚さの変化について、図１４に示すバネモデルを用いて説明する。二次電池１０は、電解液が含浸されていない部分（非含浸部分という）と、電解液が含浸された部分（含浸部分という）とで構成されている。非含浸部分は、電池ケース１１０および集電箔１２１ａ，１２２ａである。含浸部分は、負極活物質層１２１ｂ、正極活物質層１２２ｂおよびセパレータ１２３である。図１４において、バネ定数ｋ１は、非含浸部分のバネ定数（固定値）であり、バネ定数ｋ２は、含浸部分のバネ定数（固定値）である。
【０１３４】
図１４に示す初期状態とは、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３が膨張する前の状態である。ここで、非含浸部分の厚さは初期値を示す。二次電池１０に拘束力が与えられていない状態（非拘束状態、すなわち自然状態）では、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３の膨張によって、発電要素１２０の厚さが初期状態の厚さよりも変化量Δｙ０だけ増加する。ここで、非含浸部分の厚さは初期値のままである。
【０１３５】
一方、二次電池１０に拘束力が与えられている状態（拘束状態）では、変化量Δｙ０に相当する力と拘束力とが釣り合う。拘束力によって発電要素１２０が縮められ、拘束状態での発電要素１２０の厚さは、非拘束状態での発電要素１２０の厚さよりも小さくなる。拘束状態において、含浸部分の厚さは、初期状態の厚さよりも変化量Δｙの分だけ大きくなる。ここで、変化量Δｙは変化量Δｙ０よりも小さい。
【０１３６】
図１４に示す拘束状態では、下記式（１６）に示す関係（力の釣り合い関係）が成り立つ。
【０１３７】
【数１３】
2016039010.tif 000015
【０１３８】
上記式（１６）において、Ｈは発電要素１２０の高さ（固定値）であり、ｐは拘束状態における電解液の圧力であり、ｘはＸ方向における発電要素１２０の位置である。バネ定数ｋ１，ｋ２は実験などによって予め求めておくことができる。高さＨとは、図２や図４において、Ｚ方向における発電要素１２０の長さであり、予め求めておくことができる。
【０１３９】
変化量Δｙ０は、電池温度Ｔｂや二次電池１０のＳＯＣに依存する。このため、電池温度ＴｂおよびＳＯＣの少なくとも一方と、変化量Δｙ０との対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、電池温度ＴｂやＳＯＣを特定することにより変化量Δｙ０を特定できる。変化量Δｙ０を特定すれば、変化量Δｙは、圧力ｐを変数として含む関数で表すことができる。
【０１４０】
変化量Δｙは、発電要素１２０の全体における厚さの変化量である。ここで、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３のそれぞれにおける厚さの変化量Δｙｐ，Δｙｎ，Δｙｓは、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３におけるバネ定数の比と、変化量Δｙとから算出できる。具体的には、変化量Δｙｐ，Δｙｎ，Δｙｓは、下記式（１７）によって表される。
【０１４１】
【数１４】
2016039010.tif 000016
【０１４２】
上記式（１７）において、ｋｎは負極活物質層１２１ｂのバネ定数（固定値）、ｋｐは正極活物質層１２２ｂのバネ定数（固定値）、ｋｓはセパレータ１２３のバネ定数（固定値）である。各バネ定数ｋｎ，ｋｐ，ｋｓは実験などによって予め求めておくことができる。
【０１４３】
変化量Δｙｐ，Δｙｎ，Δｙｓを用いれば、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３のそれぞれにおける厚さｙｐ，ｙｎ，ｙｓを把握できる。すなわち、前回の演算周期で算出された厚さｙｊに対して、上記式（１７）で算出された変化量Δｙｊ（Δｙｐ、Δｙｎ、Δｙｓ）を加算することにより、今回の厚さｙｊを算出できる。
【０１４４】
また、変化量Δｙｐ，Δｙｎ，Δｙｓを用いれば、正極活物質層１２２ｂ、負極活物質層１２１ｂおよびセパレータ１２３のそれぞれにおいて、電解液が存在する空間の体積の変化量ΔＶｅ，ｊを算出できる。ここで、Ｖｅ，ｊは、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３のそれぞれにおいて、電解液が存在する空間の体積である。変化量ΔＶｅ，ｊは下記式（１８）によって表される。
【０１４５】
【数１５】
2016039010.tif 000017
【０１４６】
上記式（１８）において、Ｖｓ，ｊは、各活物質層１２１ｂ，１２２ｂでは、各活物質層１２１ｂ，１２２ｂに含まれる活物質の体積（Ｖｓ，ｎ、Ｖｓ，ｐ）を示し、セパレータ１２３では、セパレータ１２３自体の体積（Ｖｓ，ｓ）を示す。そして、体積Ｖｓ，ｊの変化量がΔＶｓ，ｊとなる。Ｖａｌｌ，ｊは体積Ｖｓ，ｊおよび体積Ｖｅ，ｊの総和である。そして、体積Ｖａｌｌ，ｊの変化量がΔＶａｌｌ，ｊとなる。
【０１４７】
変化量ΔＶａｌｌ，ｊおよび変化量ΔＶｓ，ｊは、下記式（１９），（２０）によって表される。
【０１４８】
【数１６】
2016039010.tif 000018
【０１４９】
上記式（１９）に基づいて変化量ΔＶａｌｌ，ｊを算出するとき、体積Ｖａｌｌ，ｊおよび厚さｙｊとしては、前回の演算周期で算出された値が用いられる。また、変化量Δｙｊとしては、上記式（１７）が用いられる。上記式（２０）に示す体積膨張率βｊは、実施例２で説明したように、電池温度Ｔｂや二次電池１０のＳＯＣから算出できる。また、上記式（２０）に示す体積分率εｅ，ｊおよび体積Ｖａｌｌ，ｊとしては、前回の演算周期で算出された値が用いられる。
【０１５０】
上記式（１９），（２０）に基づいて変化量ΔＶａｌｌ，ｊ，ΔＶｓ，ｊを算出すれば、上記式（１８）に基づいて変化量ΔＶｅ，ｊを算出できる。ここで、上記式（２０）には体積膨張率βｊが含まれているため、体積膨張率βｊに基づいて変化量ΔＶｅ，ｊを算出できる。
【０１５１】
電解液の流速ｕｊを算出する式（上記式（１），（１０）〜（１２））では、変化量ΔＶｅ，ｊがパラメータとして含まれていない。ただし、電解液が存在する空間の体積が変化するときには、電解液の体積分率εｅ，ｊが変化するため、変化量ΔＶｅ，ｊを変化量Δεｅ，ｊに変換することができる。具体的には、体積Ｖａｌｌ，ｊが予め定めた一定の体積Ｖａｌｌ＿ｍｏｄｅｌ，ｊであると仮定したとき、下記式（２１）に示すように、体積Ｖａｌｌ＿ｍｏｄｅｌ，ｊおよび変化量ΔＶｅ，ｊから変化量Δεｅ，ｊを算出できる。
【０１５２】
【数１７】
2016039010.tif 000019
【０１５３】
変化量Δεｅ，ｊを算出すれば、前回の演算周期で算出された体積分率εｅ，ｊに変化量Δεｅ，ｊを加算することにより、今回の体積分率εｅ，ｊを算出できる。今回算出した体積分率εｅ，ｊを用いれば、上記式（１），（２）に基づいて、又は、上記式（１０）〜（１２）に基づいて、流速ｕｊを算出できる。また、上記式（１３）に基づいて塩濃度ｃｅ，ｊを算出するときにも、今回算出した体積分率εｅ，ｊを用いることができる。
【０１５４】
本実施例によれば、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３のそれぞれにおける体積（具体的には厚さ）が変化したときを考慮して、体積分率εｅ，ｊを算出している。これにより、上述した体積が変化するときにおいて、流速ｕｊの算出精度を向上させることができる。これに伴い、塩濃度分布の算出精度を向上させたり、抵抗上昇量Ｒｈの算出精度を向上させたりすることができる。
【０１５５】
透過係数Ｋｊは、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３のそれぞれにおける厚さｙｊ（すなわち、厚さｙｐ，ｙｎ，ｙｓ）に依存する。例えば、正極活物質層１２２ｂにおいて、正極活物質の体積が一定であるとき、正極活物質層１２２ｂの厚さｙｐが大きくなるほど、正極活物質層１２２ｂ中で電解液が移動しやすくなるため、透過係数Ｋｐが大きくなる。
【０１５６】
また、正極活物質の体積が一定であるとき、正極活物質層１２２ｂの厚さｙｐが小さくなるほど、正極活物質層１２２ｂ中で電解液が移動しにくくなるため、透過係数Ｋｐが小さくなる。同様に、セパレータ１２３自体の体積が一定であるとき、セパレータ１２３の厚さｙｓが大きくなるほど、セパレータ１２３中で電解液が移動しやすくなるため、セパレータ１２３の透過係数Ｋｓが大きくなる。また、セパレータ１２３自体の体積が一定であるとき、セパレータ１２３の厚さｙｓが小さくなるほど、セパレータ１２３中で電解液が移動しにくくなるため、透過係数Ｋｓが小さくなる。
【０１５７】
したがって、透過係数Ｋｊおよび厚さｙｊの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、厚さｙｊに応じた透過係数Ｋｊを特定できる。このように特定した透過係数Ｋｊは、上記式（１）又は上記式（１０）〜（１２）において用いることができる。
【０１５８】
一方、透過係数Ｋｊは体積膨張率βｊに依存する。例えば、正極活物質層１２２ｂにおいて、厚さｙｐが一定であるとき、正極活物質が膨張するほど、正極活物質層１２２ｂ中で電解液が移動しにくくなるため、透過係数Ｋｐが小さくなる。また、厚さｙｐが一定であるとき、正極活物質が収縮するほど、正極活物質層１２２ｂ中で電解液が移動しやすくなるため、透過係数Ｋｐが大きくなる。同様に、セパレータ１２３の厚さｙｓが一定であるとき、セパレータ１２３自体が膨張するほど、セパレータ１２３中で電解液が移動しにくくなるため、透過係数Ｋｓが小さくなる。また、厚さｙｓが一定であるとき、セパレータ１２３自体が収縮するほど、セパレータ１２３中で電解液が移動しやすくなるため、透過係数Ｋｓが大きくなる。
【０１５９】
したがって、透過係数Ｋｊおよび体積膨張率βｊの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、体積膨張率βｊに応じた透過係数Ｋｊを特定できる。このように特定した透過係数Ｋｊは、上記式（１）又は上記式（１０）〜（１２）において用いることができる。なお、体積膨張率βｊは、実施例２で説明したように、電池温度Ｔｂや二次電池１０のＳＯＣに基づいて算出できる。
【０１６０】
透過係数Ｋｊを特定するときには、厚さｙｊおよび体積膨張率βｊの少なくとも一方を考慮できる。厚さｙｊおよび体積膨張率βｊの両方を考慮するときには、厚さｙｊおよび体積膨張率βｊと、透過係数Ｋｊとの対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけばよい。
【０１６１】
一方、上述した説明では、上記式（１６）に示すバネ定数ｋ１，ｋ２を固定値としているが、バネ定数ｋ１，ｋ２を変更することもできる。バネ定数ｋ１，ｋ２を規定する部分（含浸部分および非含浸部分）が図１４に示す非拘束状態から縮められると、バネ定数ｋ１，ｋ２が大きくなる。含浸部分が膨張するほど、非含浸部分が縮められるため、変化量Δｙに基づいて、非含浸部分のバネ定数ｋ１を算出できる。具体的には、変化量Δｙおよびバネ定数ｋ１の対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、変化量Δｙを算出することにより、バネ定数ｋ１を算出できる。バネ定数ｋ１を算出するときには、前回の演算周期で算出された変化量Δｙを用いることができる。
【０１６２】
図１４を用いて説明したように、含浸部分は、変化量Δｙ０および変化量Δｙの差（Δｙ０−Δｙ）に相当する分だけ、縮められることになる。したがって、この差（Δｙ０−Δｙ）に基づいてバネ定数ｋ２を算出できる。具体的には、差（Δｙ０−Δｙ）およびバネ定数ｋ２の対応関係を示す情報（マップ又は演算式）を実験などによって予め用意しておけば、差（Δｙ０−Δｙ）を算出することにより、バネ定数ｋ２を算出できる。バネ定数ｋ２を算出するときには、前回の演算周期で算出された差（Δｙ０−Δｙ）を用いることができる。
【０１６３】
含浸部分は、負極活物質層１２１ｂ、正極活物質層１２２ｂおよびセパレータ１２３となるため、バネ定数ｋ２にはバネ定数ｋｎ，ｋｐ，ｋｓが含まれる。各バネ定数ｋｎ，ｋｐ，ｋｓは、バネ定数ｋ２を算出する場合と同様に、差（Δｙ０，ｊ−Δｙｊ）に基づいて算出できる。ここで、変化量Δｙ０，ｊは、変化量Δｙ０のうち、各活物質層１２１ｂ，１２２ｂおよびセパレータ１２３のそれぞれが占める変化量である。各バネ定数ｋｎ，ｋｐ，ｋｓを算出するときには、前回の演算周期で算出された差（Δｙ０，ｊ−Δｙｊ）を用いることができる。このように算出された各バネ定数ｋｎ，ｋｐ，ｋｓは上記式（１７）で用いられる。
【符号の説明】
【０１６４】
１０：二次電池、２０：負荷、３１：電圧センサ、３２：電流センサ、
３３：温度センサ、４０：コントローラ、４１：メモリ、１１０：電池ケース、
１１１：負極端子、１１２：正極端子、１２０：発電要素、１２１：負極板、
１２１ａ：集電箔、１２１ｂ：負極活物質層、１２２：正極板、１２２ａ：集電箔、
１２２ｂ：正極活物質層、１２３：セパレータ
【図１】
2016039010.tif 000020
【図２】
2016039010.tif 000021
【図３】
2016039010.tif 000022
【図４】
2016039010.tif 000023
【図５】
2016039010.tif 000024
【図６】
2016039010.tif 000025
【図７】
2016039010.tif 000026
【図８】
2016039010.tif 000027
【図９】
2016039010.tif 000028
【図１０】
2016039010.tif 000029
【図１１】
2016039010.tif 000030
【図１２】
2016039010.tif 000031
【図１３】
2016039010.tif 000032
【図１４】
2016039010.tif 000033