(57)【要約】      （修正有）【課題】電池に詳しくない技術者であっても蓄電デバイスの熱現象を容易にシミュレートできるシミュレーション方法、シミュレーション装置及びコンピュータプログラムを提供する。【解決手段】シミュレーションシステムにおいて方法は、クライアント装置が、蓄電デバイスに関するシミュレーション条件を受付け、受付けたシミュレーション条件に基づき、サーバ装置が、短絡電流を計算して蓄電デバイスから外部への熱現象をシミュレートする。蓄電デバイスに関する条件は、蓄電デバイスにおける内部短絡の発生場所及び短絡部の抵抗又は外部短絡における抵抗値に関連する情報又は外部から加熱する加熱箇所、熱量及び環境温度を含む。【選択図】図７
2020161115.tif 000002
【特許請求の範囲】
【請求項１】
蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートする
シミュレーション方法。
【請求項２】
前記シミュレーション条件は、前記蓄電デバイスにおける内部短絡の発生箇所、及び短絡部の抵抗値を含み、
前記内部短絡に伴う前記熱現象をシミュレートする
請求項１に記載のシミュレーション方法。
【請求項３】
前記蓄電デバイスは、巻回セルを含み、
前記巻回電極体を仮想的に展開した状態にて短絡電流を計算する
請求項２に記載のシミュレーション方法。
【請求項４】
前記シミュレーション条件は、前記蓄電デバイスの外部短絡における抵抗値に関連する情報を含み、
前記外部短絡に伴う前記熱現象をシミュレートする
請求項１に記載のシミュレーション方法。
【請求項５】
前記シミュレーション条件は、前記蓄電デバイスを外部から加熱する際の加熱箇所、熱量、及び環境温度を含み、
前記蓄電デバイスの加熱に伴う前記熱現象をシミュレートする
請求項１から請求項４の何れか１つに記載のシミュレーション方法。
【請求項６】
蓄電デバイスに関するシミュレーション条件を受付け、
前記シミュレーション条件は、前記蓄電デバイスを外部から加熱する際の加熱箇所を含み、
受付けたシミュレーション条件に基づき、前記蓄電デバイスの加熱に伴う前記蓄電デバイスから外部への熱現象をシミュレートする
シミュレーション方法。
【請求項７】
前記蓄電デバイスの示差熱分析によって得られる温度と発熱量との関係を定式化し、
定式化して得られる温度と発熱量との関係式に基づき、前記蓄電デバイスの材料分解反応における発熱速度を計算する
請求項１又は請求項６に記載のシミュレーション方法。
【請求項８】
蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、前記蓄電デバイスの材料分解反応に伴うガスの発生をシミュレートする
シミュレーション方法。
【請求項９】
前記材料分解反応の反応速度に基づき、前記ガスの発生速度、及び熱量の発生速度の少なくとも１つを計算する
請求項８に記載のシミュレーション方法。
【請求項１０】
前記ガスの発生速度は、前記材料分解反応の反応速度に比例するように計算される
請求項９に記載のシミュレーション方法。
【請求項１１】
第一の蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、短絡電流を計算して前記第一の蓄電デバイスから外部への熱現象をシミュレートし、
前記第一の蓄電デバイスから外部への熱現象による、第二の蓄電デバイスの加熱に伴う、前記第二の蓄電デバイスから外部への熱現象をシミュレートする
シミュレーション方法。
【請求項１２】
第一の蓄電デバイスに関するシミュレーション条件を受付け、
前記シミュレーション条件は、前記第一の蓄電デバイスを外部から加熱する際の加熱箇所を含み、
受付けたシミュレーション条件に基づき、前記第一の蓄電デバイスの加熱に伴う前記第一の蓄電デバイスから外部への熱現象をシミュレートし、
前記第一の蓄電デバイスから外部への熱現象による、第二の蓄電デバイスの加熱に伴う、前記第二の蓄電デバイスから外部への熱現象をシミュレートする
シミュレーション方法。
【請求項１３】
外部端末を用いたユーザ認証の後に前記外部端末から送信されるシミュレーション条件を受信し、
受信したシミュレーション条件に基づくシミュレーション結果または前記シミュレーション条件に基づくシミュレーションプログラムを前記外部端末へ送信する
請求項１から請求項１２の何れか１つに記載のシミュレーション方法。
【請求項１４】
蓄電デバイスに関するシミュレーション条件を受付ける受付部と、
受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートするシミュレーション実行部と、
前記シミュレーション実行部によるシミュレーション結果または前記シミュレーション条件に基づくシミュレーションプログラムを出力する出力部と
を備えるシミュレーション装置。
【請求項１５】
コンピュータに、
蓄電デバイスに関するシミュレーション条件を受付けるユーザインタフェースを提示させ、
受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートする
処理を実行させるためのコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、コンピュータ実装されるシミュレーション方法、シミュレーション装置、及びコンピュータプログラムに関する。
【背景技術】
【０００２】
近年、自動車業界を始めとした各業界で、ＭＢＤ（モデルベース開発）が盛んに導入されており、シミュレーションに基づいた製品開発が浸透している（例えば、特許文献１を参照）。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開平１１−１４５０７号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
モデルベース開発において、例えば、開発要素の１つである特定の蓄電デバイスについて熱的な安全性をシミュレートする場合が考えられる。蓄電デバイスの内部にて起こる物理現象や化学現象に応じて様々な条件設定をする必要がある。蓄電デバイスの安全性に関わる事象は、化学反応、伝熱、電流、電気化学、流体力学などの複数の物理現象が互いに関連しており、メカニズムや物性値が未知である場合も少なくない。このため、電池に詳しくない技術者が蓄電デバイスの安全性をシミュレートすることは困難である。しかし、昨今の電気自動車、再生可能エネルギ、スマートグリッドなどの目覚ましい開発進展を考え、高性能・高安全な蓄電デバイスに対する期待は大きく、シミュレーションを高度に活用した安全性設計の意義は大きい。
【０００５】
本発明は、斯かる事情に鑑みてなされたものであり、電池に詳しくない技術者であっても蓄電デバイスの熱現象を容易にシミュレートできるシミュレーション方法、シミュレーション装置、及びコンピュータプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００６】
シミュレーション方法は、蓄電デバイスに関するシミュレーション条件を受付け、受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートする。
【０００７】
シミュレーション装置は、蓄電デバイスに関するシミュレーション条件を受付ける受付部と、受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートするシミュレーション実行部と、前記シミュレーション実行部によるシミュレーション結果または前記シミュレーション条件に基づくシミュレーションプログラムを出力する出力部とを備える。
【０００８】
コンピュータプログラムは、コンピュータに、蓄電デバイスに関するシミュレーション条件を受付けるユーザインタフェースを提示させ、受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートする処理を実行させるためのコンピュータプログラムである。
【発明の効果】
【０００９】
上記構成によれば、電池に詳しくない技術者であっても蓄電デバイスの熱現象を容易にシミュレートできる。
【図面の簡単な説明】
【００１０】
【図１】実施の形態に係るシミュレーションシステムの全体構成を説明する模式図である。
【図２】サーバ装置の内部構成を説明するブロック図である。
【図３】電池テーブルの一例を示す概念図である。
【図４】クライアント装置の内部構成を説明するブロック図である。
【図５】シミュレーション条件を受付ける受付画面の一例を示す模式図である。
【図６】シミュレーション手法の概略を説明する説明図である。
【図７】サーバ装置及びクライアント装置が実行する処理の手順を説明するフローチャートである。
【図８】実施の形態３に係る蓄電デバイスの構成部品を説明する説明図である。
【図９】巻回セルの、電極タブを有するタイプの電極体を示す図である。
【図１０】巻回電極体における貫通短絡での短絡部の出現箇所を説明する説明図である。
【図１１】巻回電極体における部分短絡での短絡部の出現箇所を説明する説明図である。
【図１２】複数の蓄電デバイスを含むシステム（組電池）を示す模式図である。
【図１３】図１２のシステムにおいて複数の蓄電デバイスが連鎖的にガスを排出する様子を可視化したグラフィック表示（動画表示）である。
【発明を実施するための形態】
【００１１】
シミュレーション方法は、蓄電デバイスに関するシミュレーション条件を受付け、受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートする。外部とは、例えば蓄電デバイスまたは複数の蓄電デバイスからなる蓄電システムを収納したコンテナや、車両の蓄電池格納部、或いは蓄電デバイスを収納した建屋の外部の空間であってもよい。
この構成によれば、受付けたシミュレーション条件に基づき、蓄電デバイスから外部への熱現象をシミュレートできる。蓄電デバイスから外部への熱現象は、例えば、内部短絡や外部短絡に伴うジュール熱及び材料分解反応による発熱現象、蓄電デバイスを外部から加熱した場合の材料分解反応による発熱現象、及び蓄電デバイスの発熱に伴うガスの発生を含む。本シミュレーション方法は、蓄電デバイスの発熱速度、ガスの発生速度等をシミュレーション条件に基づき計算する。
【００１２】
シミュレーション方法は、外部端末を用いたユーザ認証の後に前記外部端末から送信されるシミュレーション条件を受信し、受信したシミュレーション条件に基づくシミュレーションの実行結果を前記外部端末へ送信してもよい。この構成によれば、ユーザが蓄電デバイスの挙動を表す理論に精通していない場合であっても、シミュレーション条件を受付けるだけで、蓄電デバイスの外部に現れる熱現象のシミュレーション結果またはシミュレーションプログラムをユーザに提供できる。
【００１３】
前記シミュレーション条件は、前記蓄電デバイスにおける内部短絡の発生箇所を含んでもよく、シミュレーション方法は、前記内部短絡に伴う前記熱現象をシミュレートしてもよい。前記シミュレーション条件は、短絡部の抵抗値に関連する情報を更に含んでもよい。短絡部の抵抗値に関連する情報は、内部短絡を生じさせた物質の名称（例えば、ニッケル、鉄）と形状（例えば、円形、四角形、それらのサイズ）、内部短絡のモード（例えば、クラッシュ、釘刺し、異物混入）などを含んでもよく、その情報から短絡部の抵抗値を計算できるものであることが好ましい。短絡部の抵抗値に関連する情報は、上述の間接的な情報に代えて、短絡部と蓄電デバイスを構成する部材（例えば正極集電箔など）との接触抵抗など、短絡部の抵抗値を直接的に示すものでもよい。この構成によれば、内部短絡の発生箇所を与えることによって、異なる発生箇所での内部短絡に伴う熱現象をシミュレートできる。
【００１４】
前記蓄電デバイスは、巻回電極体を含んでもよく、シミュレーション方法は、前記巻回電極体を仮想的に展開した状態にて短絡電流を計算してもよい。この構成によれば、巻回セルの構造を考慮して短絡電流を計算するので、貫通短絡又は部分短絡に応じた熱現象を的確にシミュレートできる。
【００１５】
前記シミュレーション条件は、前記蓄電デバイスの外部短絡における抵抗値に関連する情報を含んでもよく、シミュレーション方法は、前記外部短絡に伴う前記熱現象をシミュレートしてもよい。この構成によれば、外部短絡における抵抗値に関連する情報を与えることによって、外部短絡に伴う熱現象をシミュレートできる。
【００１６】
他の実施形態に係るシミュレーション方法は、蓄電デバイスに関するシミュレーション条件を受付け、前記シミュレーション条件は、前記蓄電デバイスを外部から加熱する際の加熱箇所を含み、受付けたシミュレーション条件に基づき、前記蓄電デバイスの加熱に伴う前記蓄電デバイスから外部への熱現象をシミュレートする。前記シミュレーション条件は、外部加熱の熱量を更に含んでもよい。前記シミュレーション条件は、環境温度を更に含んでもよい。この構成によれば、蓄電デバイスを外部から加熱する際の加熱箇所を与えることによって、加熱に伴う熱現象をシミュレートできる。
【００１７】
シミュレーション方法は、前記蓄電デバイスにおける電気化学反応と、前記蓄電デバイスの材料分解反応における発熱反応との連成解析により前記熱現象をシミュレートしてもよい。電池の発熱反応及びジュール発熱は独立した物理現象ではなく、伝熱といった物理現象を介して相互に関連しながら進行するので、連成解析を行うことにより、電池内部で起こっている現象を的確に反映させて電池の挙動をシミュレートすることができる。
【００１８】
シミュレーション方法は、前記蓄電デバイスの示差熱分析によって得られる温度と発熱量との関係を定式化し、定式化して得られる温度と発熱量との関係式に基づき、前記蓄電デバイスの材料分解反応における発熱速度を計算してもよい。この構成によれば、示差熱分析によって得られる温度と発熱量との関係を定式化するので、温度及び発熱量の関係を、時刻及び発熱量の関係に換算する手間が省ける。
【００１９】
他の実施形態に係るシミュレーション方法は、蓄電デバイスに関するシミュレーション条件を受付け、受付けたシミュレーション条件に基づき、前記蓄電デバイスの材料分解反応に伴うガスの発生をシミュレートする。蓄電デバイスにおいて内部短絡などの事象が発生し、安全機構がうまく機能しなかった場合、材料分解反応が進行して蓄電デバイスの内部から高温のガスが噴出する可能性がある。本シミュレーション方法では、材料分解反応に伴うガスの発生をシミュレートできる。
【００２０】
シミュレーション方法は、前記材料分解反応の反応速度に基づき、前記ガスの発生速度、及び熱量の発生速度の少なくとも１つを計算してもよい。前記ガスの発生速度は、前記材料分解反応の反応速度に比例するように計算されてもよい。この構成によれば、材料分解反応の反応速度に基づき、ガスの発生速度、及び熱量の発生速度の少なくとも１つが計算される。
【００２１】
シミュレーション装置は、蓄電デバイスに関するシミュレーション条件を受付ける受付部と、受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートするシミュレーション実行部と、前記シミュレーション実行部によるシミュレーション結果または前記シミュレーション条件に基づくシミュレーションプログラムを出力する出力部とを備える。
この構成によれば、受付けたシミュレーション条件に基づき、蓄電デバイスから外部への熱現象をシミュレートできる。蓄電デバイスから外部への熱現象は、例えば、内部短絡や外部短絡に伴う材料分解反応による発熱現象、蓄電デバイスを外部から加熱した場合の材料分解反応による発熱現象、及び蓄電デバイスの発熱に伴うガスの発生を含む。本シミュレーション装置は、蓄電デバイスの発熱速度、ガスの発生速度等をシミュレーション条件に基づき計算する。
【００２２】
コンピュータプログラムは、コンピュータに、蓄電デバイスに関するシミュレーション条件を受付けるユーザインタフェースを提示させ、受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートする処理を実行させる。
この構成によれば、受付けたシミュレーション条件に基づき、蓄電デバイスから外部への熱現象をシミュレートできる。蓄電デバイスから外部への熱現象は、例えば、内部短絡や外部短絡に伴う材料分解反応による発熱現象、蓄電デバイスを外部から加熱した場合の材料分解反応による発熱現象、及び蓄電デバイスの発熱に伴うガスの発生を含む。本コンピュータプログラムは、蓄電デバイスの発熱速度、ガスの発生速度等をシミュレーション条件に基づき計算する。
【００２３】
以下、本発明をその実施の形態を示す図面に基づいて具体的に説明する。
（実施の形態１）
図１は実施の形態に係るシミュレーションシステムの全体構成を説明する模式図である。本実施の形態に係るシミュレーションシステムは、通信網Ｎを介して互いに通信可能に接続されるサーバ装置１００とクライアント装置２００とを備える。サーバ装置１００は、クライアント装置２００からの要求に応じて、蓄電デバイスから外部に現れる熱現象をシミュレートし、シミュレーション結果をクライアント装置２００へ提供する。ここで、シミュレーション対象の蓄電デバイスは、鉛蓄電池、リチウムイオン電池等の二次電池、又はキャパシタ等の再充電可能な蓄電素子（セル）を含む。また、シミュレーション対象の蓄電デバイスは、複数のセルを直列に接続したモジュール、複数のモジュールを直列に接続したバンク（直列接続電池群）、複数のバンクを並列に接続したドメイン（並列接続電池群）等を含んでもよい。
【００２４】
クライアント装置２００は、ユーザによって利用されるパーソナルコンピュータ、スマートフォン、タブレット端末などの端末装置である。クライアント装置２００には、サーバ装置１００にアクセスするためのソフトウェア（アプリケーションプログラム）がインストールされている。サーバ装置１００は、クライアント装置２００からのアクセスを受付けた際に例えばユーザＩＤ及びパスワードに基づくユーザ認証を行い、ユーザ認証に成功した場合、クライアント装置２００に対して適宜のサービスを提供する。
【００２５】
本実施の形態に係るサーバ装置１００は、ユーザ認証の後、クライアント装置２００のユーザによる各種入力を受付けるためのインタフェース画面をクライアント装置２００へ送信する。このインタフェース画面には、例えば、シミュレーション条件を受付けるための受付画面が含まれる。サーバ装置１００は、受付けた条件に基づいて実行したシミュレーション結果をクライアント装置２００へ送信する。
【００２６】
サーバ装置１００がクライアント装置２００に対して送信するシミュレーション結果は、シミュレーションの実行結果として得られる数値データ、グラフ等のデータを含む。サーバ装置１００がクライアント装置２００に対して送信するシミュレーション結果は、シミュレーションの実行結果として得られる数理モデル、または、シミュレーション条件に基づくシミュレーションプログラムを含んでもよい。
【００２７】
本実施の形態では、クライアント装置２００においてシミュレーション条件を受付け、受付けたシミュレーション条件等をサーバ装置１００へ送信してシミュレーションを実行する。代替的に、サーバ装置１００において、シミュレーション条件を受付け、受付けたシミュレーション条件等に基づきシミュレーションを実行し、サーバ装置１００にてシミュレーション結果を表示してもよい。
【００２８】
図２はサーバ装置１００の内部構成を説明するブロック図である。サーバ装置１００は、制御部１０１、記憶部１０２、通信部１０３、操作部１０４及び表示部１０５を備える。
【００２９】
制御部１０１は、ＣＰＵ（Central Processing Unit）、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）などにより構成される。制御部１０１が備えるＣＰＵは、ＲＯＭ又は記憶部１０２に記憶されている各種コンピュータプログラムをＲＡＭ上に展開して実行することにより、装置全体を本願のシミュレーション装置として機能させる。サーバ装置１００は、シミュレーション装置の一実施形態に過ぎず、クライアント装置２００と通信可能に接続された任意の情報処理装置であればよい。
【００３０】
制御部１０１は、上記の構成に限定されるものではなく、複数のＣＰＵ、マルチコアＣＰＵ、ＧＰＵ（Graphics Processing Unit）、マイコン、揮発性又は不揮発性のメモリ等を備える任意の処理回路又は演算回路であってもよい。制御部１０１は、計測開始指示を与えてから計測終了指示を与えるまでの経過時間を計測するタイマ、数をカウントするカウンタ、日時情報を出力するクロック等の機能を備えていてもよい。
【００３１】
記憶部１０２は、ＨＤＤ（Hard Disk Drive）、ＳＳＤ（Solid State Drive）等を用いた記憶装置を備える。記憶部１０２には、制御部１０１によって実行される各種コンピュータプログラム、及びコンピュータプログラムの実行に必要なデータ等が記憶される。記憶部１０２に記憶されるコンピュータプログラムは、蓄電デバイスから外部に現れる熱現象をシミュレートするシミュレーションプログラムを含む。シミュレーションプログラムは、例えば実行バイナリである。シミュレーションプログラムの元となる理論式は、蓄電デバイスから外部に現れる熱現象を表す代数方程式又は微分方程式によって記述される。シミュレーションプログラムは単一のコンピュータプログラムであってもよく、複数のコンピュータプログラムから構成されるプログラム群であってもよい。シミュレーションプログラムは、ＭＡＴＬＡＢ（登録商標）、Ａｍｅｓｉｍ（登録商標）、Ｔｗｉｎ  Ｂｕｉｌｄｅｒ（登録商標）、ＭＡＴＬＡＢ＆Ｓｉｍｕｌｉｎｋ（登録商標）、Ｓｉｍｐｌｏｒｅｒ（登録商標）、ＡＮＳＹＳ（登録商標）、Ａｂａｑｕｓ（登録商標）、Ｍｏｄｅｌｉｃａ（登録商標）、ＶＨＤＬ−ＡＭＳ（登録商標）、Ｃ言語、Ｃ＋＋、Ｊａｖａ（登録商標）などの市販の数値解析ソフトウェア又はプログラミング言語によって記述されてもよい。数値解析ソフトウェアは、１Ｄ−ＣＡＥと称される回路シミュレータ、または、３Ｄ形状による有限要素法や有限体積法などのシミュレーションであってもよい。これらに基づいた縮退モデル（ＲＯＭ：Ｒｅｄｕｃｅｄ−Ｏｒｄｅｒ  Ｍｏｄｅｌ）を用いてもよい。
【００３２】
記憶部１０２に記憶されるプログラムは、当該プログラムを読み取り可能に記録した非一時的な記録媒体Ｍにより提供されてもよい。記録媒体Ｍは、例えば、ＣＤ−ＲＯＭ、ＵＳＢ（Universal Serial Bus）メモリ、ＳＤ（Secure Digital）カード、マイクロＳＤカード、コンパクトフラッシュ（登録商標）などの可搬型メモリである。この場合、制御部１０１は、不図示の読取装置を用いて記録媒体Ｍからプログラムを読み取り、読み取ったプログラムを記憶部１０２にインストールする。記憶部１０２に記憶されるプログラムは、通信部１０３を介した通信により提供されてもよい。この場合、制御部１０１は、通信部１０３を通じてプログラムを取得し、取得したプログラムを記憶部１０２にインストールする。
【００３３】
記憶部１０２には、シミュレーションの結果として得られる数理モデルが記憶されてもよい。数理モデルは、例えば、プログラミング言語又は数値解析ソフトウェアにより実行される実行コードである。数理モデルは、プログラミング言語又は数値解析ソフトウェアにより参照される、定義情報若しくはライブラリファイルであってもよい。
【００３４】
記憶部１０２は、蓄電デバイス（例えば、二次電池）の情報をユーザＩＤに関連付けて記憶する電池テーブルを有していてもよい。図３は電池テーブルの一例を示す概念図である。電池テーブルは、例えば、電池を識別する電池ＩＤ、ユーザを識別するユーザＩＤ、及び電池情報を関連付けて記憶する。電池テーブルに登録される電池情報は、例えば、正極及び負極の情報、電解液の情報、集電体の情報などを含む。正極及び負極の情報とは、正極及び負極の活物質名、厚み、幅、奥行き、開回路電位などの情報である。電解液及び集電体の情報とは、イオン種、輸率、拡散係数、導電率などの情報である。電池テーブルには、蓄電デバイスの物理的性質、動作状態、回路構成等の情報を参照するリンクが含まれてもよい。電池テーブルに記憶される情報は、サーバ装置１００の管理者によって登録されてもよく、クライアント装置２００を介してユーザによって登録されてもよい。電池テーブルに記憶されている情報は、蓄電デバイスの熱現象をシミュレートする際に、シミュレーション条件の一部として利用される。
【００３５】
通信部１０３は、通信網Ｎを通じてクライアント装置２００と通信を行うためのインタフェースを備える。通信部１０３は、クライアント装置２００へ送信すべき情報が制御部１０１から入力された場合、入力された情報をクライアント装置２００へ送信すると共に、通信網Ｎを通じて受信したクライアント装置２００からの情報を制御部１０１へ出力する。
【００３６】
操作部１０４は、キーボード、マウスなどの入力インタフェースを備えており、ユーザによる操作を受付ける。表示部１０５は、液晶ディスプレイ装置などを備えており、ユーザに対して報知すべき情報を表示する。本実施の形態では、サーバ装置１００が操作部１０４及び表示部１０５を備える構成としたが、操作部１０４及び表示部１０５は必須ではなく、サーバ装置１００の外部に接続されたコンピュータを通じて操作を受付け、通知すべき情報を外部のコンピュータへ出力する構成であってもよい。
【００３７】
図４はクライアント装置２００の内部構成を説明するブロック図である。クライアント装置２００は、パーソナルコンピュータ、スマートフォン、タブレット端末等であり、制御部２０１、記憶部２０２、通信部２０３、操作部２０４及び表示部２０５を備える。
【００３８】
制御部２０１は、ＣＰＵ、ＲＯＭ、ＲＡＭなどにより構成される。制御部２０１が備えるＣＰＵは、ＲＯＭ又は記憶部２０２に記憶されている各種コンピュータプログラムをＲＡＭ上に展開して実行することにより、装置全体の制御を実行する。
【００３９】
制御部２０１は、上記の構成に限定されるものではなく、複数のＣＰＵ、マルチコアＣＰＵ、マイコン等を含む任意の処理回路又は演算回路であってもよい。制御部２０１は、計測開始指示を与えてから計測終了指示を与えるまでの経過時間を計測するタイマ、数をカウントするカウンタ、日時情報を出力するクロック等の機能を備えていてもよい。
【００４０】
記憶部２０２は、ＥＥＰＲＯＭ（Electronically Erasable Programmable Read Only Memory）などの不揮発性メモリにより構成されており、各種コンピュータプログラム及びデータを記憶する。記憶部２０２に記憶されるコンピュータプログラムは、サーバ装置１００と情報の授受を行うために用いられる汎用又は専用のアプリケーションを含む。汎用のアプリケーションプログラムの一例は、ウェブブラウザである。ウェブブラウザを用いてサーバ装置１００にアクセスする場合、ユーザＩＤ及び認証コードを用いたユーザ認証を行うことが好ましく、ユーザ認証に成功した場合にのみ、サーバ装置１００とクライアント装置２００との間の通信を許可すればよい。
【００４１】
通信部２０３は、通信網Ｎを通じてサーバ装置１００と通信を行うためのインタフェースを備える。通信部２０３は、サーバ装置１００へ送信すべき情報が制御部２０１から入力された場合、入力された情報をサーバ装置１００へ送信すると共に、通信網Ｎを通じて受信したサーバ装置１００からの情報を制御部２０１へ出力する。
【００４２】
操作部２０４は、キーボード、マウス、タッチパネルなどの入力インタフェースを備えており、ユーザによる操作を受付ける。表示部２０５は、液晶ディスプレイ装置などを備えており、ユーザに対して報知すべき情報を表示する。本実施の形態では、クライアント装置２００が操作部２０４を備える構成としたが、クライアント装置２００にキーボード、マウス等の入力インタフェースが接続される構成であってもよい。
【００４３】
以下、サーバ装置１００において、蓄電デバイスの内部短絡に伴う熱現象をシミュレートする構成について説明する。
【００４４】
サーバ装置１００は、蓄電デバイスの内部短絡に伴う熱現象をシミュレートする場合、シミュレーション条件として、内部短絡の発生箇所、及び短絡部の抵抗値に関連する情報を受付ける。このとき、サーバ装置１００は、シミュレーション条件を受付けるための受付画面をクライアント装置２００の表示部２０５に表示させ、表示させた受付画面を通じてシミュレーション条件を受付けてもよい。
【００４５】
図５はシミュレーション条件を受付ける受付画面の一例を示す模式図である。図５に示す受付画面２１０は、シミュレーション対象の電池を選択する選択欄２１１、短絡部の抵抗値に関連する情報を入力する入力欄２１２、及び、内部短絡の発生箇所を指定する指定欄２１３を備える。選択欄２１１には、プルダウンメニュー２１１ａが配置され、プルダウンメニュー２１１ａによってシミュレーション対象の電池の種類を受付ける。入力欄２１２には、入力ボックス２１２ａが配置され、例えば操作部２０４を用いた数値入力により、抵抗値の入力を受付ける。指定欄２１３には、内部短絡の発生箇所を受付けるために、蓄電デバイス（セル）の立体的な模式図と、内部短絡の発生箇所を示すアイコン２１３ａとが表示される。代替的に、モジュール、バンク、又はドメインの立体的な模式図を表示し、表示した模式図において短絡箇所を受け付けてもよい。受付画面２１０は、操作部２０４を用いたアイコン２１３ａの移動操作（ドラッグ操作）により、内部短絡の発生箇所を受付ける。
【００４６】
受付画面２１０において、送信ボタン２１４が操作された場合、入力されたシミュレーション条件（電池に関する情報、短絡部の抵抗値に関連する情報、内部短絡の発生箇所）がクライアント装置２００からサーバ装置１００へ送信される。サーバ装置１００は、クライアント装置２００から送信されるシミュレーション条件に基づき、蓄電デバイスのシミュレーションを行う。
【００４７】
以下、蓄電デバイスのシミュレーション手法について説明する。
図６はシミュレーション手法の概略を説明する説明図である。本実施の形態に係るサーバ装置１００は、ジュール発熱と、材料分解による発熱反応とに基づき、蓄電デバイスの挙動をシミュレートする。発熱には、この他に電気化学反応に伴うエンタルピー反応熱（可逆反応熱）、電気化学反応熱（不可逆反応熱）を含めてもよい。
【００４８】
サーバ装置１００は、電気化学反応に関して、例えばＮｅｗｍａｎモデルを用いることができる。Ｎｅｗｍａｎモデルは、正極および負極の電極において、均質かつ単一径の球が近接して並んでいると想定する電気化学モデルである。Ｎｅｗｍａｎモデルは、以下において説明するＮｅｒｎｓｔ−Ｐｌａｎｃｋ式、電荷保存式、拡散方程式、Ｂｕｔｌｅｒ−Ｖｏｌｍｅｒ式、及びＮｅｒｎｓｔ式により記述される
【００４９】
Ｎｅｒｎｓｔ−Ｐｌａｎｃｋ式は、電解質や多孔電極におけるイオン泳動とイオン拡散とを解くための方程式であり、次式により表される。以降の数１から数４に表される各種パラメータは、構成材料の気孔率の関数として、バルクでの値と多孔体での値とが適宜換算されるようになっていてもよい。
【００５０】
【数１】
2020161115.tif 000003
【００５１】
ここで、ｉl は液相電流密度（Ａ／ｍ2 ）、σl,effは液相伝導率（Ｓ／ｍ）、φl は液相電位（Ｖ）、Ｒは気体定数（Ｊ／（Ｋ・ｍｏｌ））、Ｔは温度（Ｋ）、Ｆはファラデー定数（Ｃ／ｍｏｌ）、ｆは活量係数、Ｃl は電解質のイオン濃度（ｍｏｌ／ｍ3 ）、ｔ+ はカチオン輸率、ｉtot は体積あたり反応電流密度（Ａ／ｍ3 ）である。
【００５２】
電荷保存式は、活物質や集電箔での電子伝導を表す式であり、次式により表される。
【００５３】
【数２】
2020161115.tif 000004
【００５４】
ここで、φs は固相電位（Ｖ）、σs は固相伝導率（Ｓ／ｍ）、ｉs は固相電流密度（Ａ／ｍ2 ）、ｉtot は体積あたり反応電流密度（Ａ／ｍ3 ）である。
【００５５】
拡散方程式は、活粒子中での活物質の拡散を表す方程式であり、次式により表される。
【００５６】
【数３】
2020161115.tif 000005
【００５７】
ここで、Ｃs は固相中の活物質濃度（ｍｏｌ／ｍ3 ）、ｔは時刻（ｓ）、Ｄs は固相中の拡散係数（ｍ2 ／ｓ）である。
【００５８】
Ｂｕｔｌｅｒ−Ｖｏｌｍｅｒ式は、固相と液相との界面で起こる電荷移動反応での活性化過電圧を表す式、Ｎｅｒｎｓｔ式は、開回路電位の定義式であり、それぞれ次式により表される。
【００５９】
【数４】
2020161115.tif 000006
【００６０】
ここで、ｉioc は反応電流密度（Ａ／ｍ2 ）、ｉ0 は交換電流密度（Ａ／ｍ2 ）、αa ，αc は移行係数、ηは活性化過電圧（Ｖ）、Ｅeqは平衡電位（Ｖ）、Ｅ0 は標準平衡電位（Ｖ）、ａox は酸化剤濃度（ｍｏｌ／ｍ3 ）、ａred は還元剤濃度（ｍｏｌ／ｍ3 ）である。交換電流密度ｉ0 は、例えば電解液のイオン濃度や活物質の濃度の関数として定義してもよい。数４に記載した理論式ではなく、実験結果に基づいた数値を用いてもよい。例えばリチウムイオン二次電池の開回路電位においては、Ｎｅｒｎｓｔ式に代替して、充電状態（ＳＯＣ）と開回路電位（ＯＣＰまたはＯＣＶ）の実測データを用いてもよい。
【００６１】
数５に活粒子の表面における、固相中の活物質濃度と電荷移動反応に関わる活物質フラックスの関係式を示す。ｒ0 は活粒子の半径（ｍ）を表し、Ｊs は活物質のフラックス（ｍｏｌ／ｍ2 s）である。換言すれば、Ｊs は電荷移動反応によって消失生成する、単位面積単位時間当たりの活物質の量である。
【００６２】
【数５】
2020161115.tif 000007
【００６３】
数６は、活物質のフラックスＪs と反応電流密度ｉloc との関係を述べた式である。
【００６４】
【数６】
2020161115.tif 000008
【００６５】
数７は、反応電流密度ｉloc と体積あたり反応電流密度ｉtot との関係を述べた式である。
【００６６】
【数７】
2020161115.tif 000009
【００６７】
本実施の形態では、蓄電デバイスの電気化学現象を表すモデルの一例としてＮｅｗｍａｎモデルを示した。代替的には、電極を単一の活粒子で表現する単粒子モデルや、ＮＴＧＫモデルに代表されるような開回路電位及び内部抵抗を温度及び充電状態（ＳＯＣ）の関数として表した多項式モデルを用いてもよく、等価回路モデルを使用してもよい。単粒子モデルについては非特許文献『Ｃｙｃｌｅ Ｌｉｆｅ Ｍｏｄｅｌｉｎｇ ｏｆ Ｌｉｔｈｉｕｍ−Ｉｏｎ Ｂａｔｔｅｒｉｅｓ , Ｇａｎｇ Ｎｉｎｇ ａｎｄ Ｂｒａｎｋｏ Ｎ. Ｐｏｐｏｖ，Ｊｏｕｒｎａｌ ｏｆ Ｔｈｅ Ｅｌｅｃｔｒｏｃｈｅｍｉｃａｌ Ｓｏｃｉｅｔｙ, １５１  （１０） Ａ１５８４−Ａ１５９１ (２００４)』に詳しく記載されている。
【００６８】
続いて、材料分解により生じる発熱反応について説明する。蓄電デバイスを構成する物質は、昇温により材料分解などの反応が開始することが知られている。例えばリチウムイオン電池の正極や負極材料は、典型的には２００〜３００℃程度になると、材料分解が始まり、発熱と共にガスを生じる。このような温度依存性のある反応速度を表すため、当該反応を以下のようなアレニウス反応式で表現することができる。
【００６９】
【数８】
2020161115.tif 000010
【００７０】
ここで、ｒは反応速度（１／ｓ）、ｋ0 は反応速度定数（１／ｓ）、Ｅa は活性化エネルギ（Ｊ／ｍｏｌ）、Ｒは気体定数（Ｊ／（Ｋ・ｍｏｌ））、Ｔは温度（Ｋ）、ｘf は反応率、ｐ，ｑ，Ｃ0 は定数である。Ｑは発熱密度（Ｗ／ｍ3 ）、ρは密度（ｋｇ／ｍ3 ）、Ｈp は反応熱（Ｊ／ｋｇ）である。
【００７１】
蓄電デバイスにおいて、内部短絡が発生した場合、蓄電デバイス全体から短絡部に向かって電流が流れ込む。この内部短絡に伴う電流によりジュール熱が発生する。ジュール熱の発生に応じて、前記の材料分解などの発熱反応が進行する。この発熱反応によって昇温した場合、電気抵抗が変化するので、短絡部に流れ込む電流の大きさも変化する。以上のように、ジュール発熱と材料分解発熱反応とは独立した物理現象ではなく、伝熱といった物理現象を介して相互に関連しながら進行する。
【００７２】
そこで、本実施の形態に係るサーバ装置１００は、材料分解に伴う発熱反応と、ジュール発熱とを連成させて解析し、例えば昇温による電気抵抗の変化を考慮して蓄電デバイスから外部への熱現象をシミュレートする。
【００７３】
以下、サーバ装置１００及びクライアント装置２００の動作について説明する。
【００７４】
図７はサーバ装置１００及びクライアント装置２００が実行する処理の手順を説明するフローチャートである。クライアント装置２００の制御部２０１は、ユーザ認証の後にサーバ装置１００から送信される表示画面用のデータを受信し、シミュレーション条件を受付けるための受付画面２１０を表示部２０５に表示する（ステップＳ１０１）。制御部２０１は、表示部２０５に表示した受付画面２１０を通じて、シミュレーション条件を受付ける（ステップＳ１０２）。具体的には、制御部２０１は、電池の種別の選択、内部短絡における短絡部の抵抗値に関連する情報、および、内部短絡の発生箇所を受付ける。
【００７５】
次いで、制御部２０１は、シミュレーション条件の送信指示を受付けたか否かを判断する（ステップＳ１０３）。制御部２０１は、図５に示す受付画面２１０において、送信ボタン２１４が操作された場合、送信指示を受付けたと判断する。送信指示を受け付けていない場合（Ｓ１０３：ＮＯ）、制御部２０１は、送信指示を受付けるまで待機する。
【００７６】
送信指示を受付けたと判断した場合（Ｓ１０３：ＹＥＳ）、制御部２０１は、ステップＳ１０２で受付けたシミュレーション条件を通信部２０３よりサーバ装置１００へ送信する（ステップＳ１０４）。
【００７７】
サーバ装置１００は、クライアント装置２００から送信されるシミュレーション条件を通信部１０３にて受信する（ステップＳ１０５）。サーバ装置１００の制御部１０１は、通信部１０３を通じて受信した条件に基づき、シミュレーションを実行する（ステップＳ１０６）。このとき、制御部１０１は、シミュレーション対象の挙動に対応するシミュレーションプログラムを実行することによって、蓄電デバイスの熱現象をシミュレートする。ユーザによって入力されたシミュレーション条件は、シミュレーションプログラムを実行する際に適用される。
【００７８】
蓄電デバイスにおいて、内部短絡が発生した場合、蓄電デバイス全体から短絡部に向かって電流が流れ込む。この内部短絡に伴う電流によりジュール熱が発生する。ジュール熱の発生に応じて、材料分解などの発熱反応が進行する。制御部１０１は、ジュール熱反応と材料分解の発熱反応とを連成させ、蓄電デバイスにおける起電力及び内部抵抗などを例えばＮｅｗｍａｎモデルにより計算し、発熱反応における反応速度を数５に示すアレニウス型反応式により計算する。以上により、制御部１０１は、昇温による電気抵抗の変化を考慮して電流計算を実行することができ、昇温による通電の停止、発熱反応の進行により材料が分解反応により発熱する現象を関連させながら蓄電デバイスにおける熱現象をシミュレートすることができる。例えば、材料分解反応が一定程度進行した部位は導電性を失い、通電が停止すると仮定してもよい。
【００７９】
シミュレーションが完了した場合、制御部１０１は、通信部１０３を通じて、シミュレーション結果をクライアント装置２００へ送信する（ステップＳ１０７）。ステップＳ１０７で送信するシミュレーション結果は、数値データであってもよく、数値データから生成したグラフ、コンター図及び動画などであってもよい。ステップＳ１０７で送信するシミュレーション結果は、シミュレーションの結果として得られる数理モデルであってもよい。数理モデルは、単なる理論モデルではなく、選択された蓄電デバイスについてシミュレーションが実行され、各種のパラメータが調整された後のモデルを表す。数理モデルは、例えば、ＭＡＴＬＡＢ（登録商標）、Ａｍｅｓｉｍ（登録商標）、Ｔｗｉｎ  Ｂｕｉｌｄｅｒ（登録商標）、ＭＡＴＬＡＢ＆Ｓｉｍｕｌｉｎｋ（登録商標）、Ｓｉｍｐｌｏｒｅｒ（登録商標）、ＡＮＳＹＳ（登録商標）、Ａｂａｑｕｓ（登録商標）、Ｍｏｄｅｌｉｃａ（登録商標）、ＶＨＤＬ−ＡＭＳ（登録商標）、Ｃ言語、Ｃ＋＋、Ｊａｖａなどの市販の数値解析ソフトウェア又はプログラミング言語で用いられるライブラリ、モジュール等のフォーマットにより提供される。
【００８０】
クライアント装置２００は、サーバ装置１００から送信されるシミュレーション結果を通信部２０３にて受信する（ステップＳ１０８）。クライアント装置２００の制御部２０１は、受信したシミュレーション結果を表示部２０５に表示させる（ステップＳ１０９）。
【００８１】
以上のように、本実施の形態では、サーバ装置１００は、蓄電デバイスの材料分解による発熱反応と、ジュール発熱とを考慮して、蓄電デバイスから外部に現れる熱現象をシミュレートすることができる。ユーザは、蓄電デバイスにおける物理現象について詳しくなくても、クライアント装置２００から蓄電デバイスの種類、内部短絡の発生箇所、及び短絡部の抵抗値に関連する情報を入力することによって、複雑なシミュレーションの設定なしに、蓄電デバイスの熱現象に関するシミュレーション結果を取得することができる。
【００８２】
本実施の形態では、蓄電デバイスにおいて内部短絡が発生した場合における熱現象をシミュレートする構成について説明したが、外部短絡が発生した場合における熱現象もシミュレートすることが可能である。この場合、サーバ装置１００は、クライアント装置２００を通じて外部短絡における抵抗値に関連する情報（例えば、外部短絡における抵抗値を直接的に示す数値。代替的に、外部短絡における抵抗値を間接的に示す情報（具体的には、外部短絡を生じさせた原因（スパナ等の工具による正極端子と負極端子との外部短絡、配線の絶縁（被覆）破壊による電池外部での短絡、クラッシュによる電池外部での短絡、スイッチ故障による外部短絡など）））の入力を受付け、受付けた外部短絡の抵抗値に関連する情報に基づき、熱現象をシミュレートする。具体的には、サーバ装置１００の制御部１０１は、外部短絡により、電池の正極端子から負極端子に流れる電流を計算する。そして、制御部１０１は、電流に伴うジュール熱と、ジュール熱の発生に応じて材料分解などの発熱反応とを連成解析することによって、蓄電デバイスおよび蓄電デバイスから外部に現れる熱現象をシミュレートすることができる。
【００８３】
（実施の形態２）
実施の形態２では、蓄電デバイスを外部から加熱した場合の熱現象についてシミュレートする。蓄電デバイスは、周囲環境の異常に起因して、加熱される場合がある（蓄電デバイスを搭載した移動体のクラッシュ時、蓄電デバイス設備における冷却装置の故障時など）。そのような状況において蓄電デバイスがどのような挙動を示し、外部へどのような熱現象を発現するかを事前にシミュレートしたいというニーズがある。
【００８４】
蓄電デバイスが直方体形状（プリズマティックセル）である場合、蓄電デバイスの側面が加熱されるか、上面が加熱されるかに応じて、挙動が異なることがある。蓄電デバイスがパウチセルである場合や、円筒セルである場合も同様である。複数の蓄電デバイスを含むシステム（組電池）では、内方に位置する蓄電デバイスと外方に位置する蓄電デバイスとは、熱源からの距離の違いや熱ごもりの影響により、挙動が異なることがある。そこで、シミュレーションの際に、熱源による蓄電デバイスの加熱箇所を受付ける。
【００８５】
サーバ装置１００は、クライアント装置２００を通じて、蓄電デバイスを外部から加熱する際の加熱箇所、熱量、及び環境温度を含むシミュレーション条件を受付ける。加熱箇所は、図５に類似するグラフィック表示を含むユーザインタフェースを介して受け付けてもよい。サーバ装置１００は、受付けたシミュレーション条件に基づき、熱現象をシミュレートする。具体的には、サーバ装置１００の制御部１０１は、外部から与えられる熱に応じて進行する材料分解反応をアレニウス型反応の式を用いて解析することにより、蓄電デバイスから外部に現れる熱現象をシミュレートすることができる。
【００８６】
（実施の形態３）
実施の形態３では、巻回セル（巻回電極体を容器に収容したセル）の内部短絡についてシミュレーションを実行する構成について説明する。
【００８７】
図８は実施の形態３に係る蓄電デバイス１の構成部品を説明する説明図である。実施の形態３に係る蓄電デバイス１は、巻回電極体１０ｖ、正極端子１１、正極集電体１１ａ、負極端子１２、及び負極集電体１２ａを備え、例えばそれらが中空直方体状の容器２０に収容される。電極体１０ｖは、例えば、銅箔からなる負極集電箔に負極活物質を設けた帯状の負極と、負極の幅方向にずらして配置されたアルミニウム箔からなる正極集電箔に正極活物質を設けた帯状の正極との間に、多孔性の樹脂フィルムからなるセパレータを配置し、これらを巻回することにより巻回電極体１０ｖを構成する。
【００８８】
正極集電体１１ａは、正極端子１１と電極体１０ｖの正極（活物質非形成部、正極集電箔）とを電気的に接続する。負極集電体１２ａは、負極端子１２と電極体１０ｖの負極（活物質非形成部、負極集電箔）とを電気的に接続する。正極集電体１１ａは、電極体１０ｖの正極と少なくとも１点で電気的に接続されており、負極集電体１２ａは、電極体１０ｖの負極と少なくとも１点で電気的に接続されている。正極集電体１１ａ及び負極集電体１２ａの位置は、電極体１０ｖの正極及び負極との電気的接続が変わらない範囲で、様々な形態が考えられる。図８の例では、電極体１０ｖの巻回中心の向きは紙面においてＸ軸方向を向いている。代替的に、電極体１０ｖの巻回中心の向きは紙面のＹ軸方向またはＺ軸方向を向いていてもよい。
【００８９】
図９は、巻回セルの、電極タブを有するタイプの電極体を示す。負極は、帯状の負極集電箔の長手方向に延びる一辺に、間隔を隔てて負極タブ１２ｂ（活物質非形成部）が複数設けられている。正極は、帯状の正極集電箔の長手方向に延びる一辺に、間隔を隔てて正極タブ１１ｂ（活物質非形成部）が複数設けられている。負極と正極との間にセパレータを配置し、これらを巻回することにより巻回電極体を構成する。複数の負極タブ１２ｂ、複数の正極タブ１１ｂは、それぞれ束ねられて、図示しない負極集電体、正極集電体に電気的に接続される。
【００９０】
図示しないが、積層セル（シート状の複数の集電箔を積層した積層電極体を容器に収容したセル）では、短絡時に短絡層以外の電極からタブのみを介して短絡する。それに対して、巻回セルでは、短絡層からタブ（正極タブ１１ｂまたは負極タブ１２ｂ）を介して流れる電流と、隣接電極との中間にある湾曲部から流れ込む電流との２種類が存在する。このため、巻回セルでは、積層セルで用いられる計算手法を流用することはできない。
【００９１】
そこで、本実施の形態では、巻回電極体を仮想的に展開して広げた状態にして短絡電流の計算を行う。内部短絡のシミュレーションにおいて実施される典型例として、釘刺し試験に代表される貫通短絡と、ニッケル片混入試験に代表される部分短絡とがある。以降ではそれぞれについて詳細を説明する。
【００９２】
本実施の形態では、貫通短絡および部分短絡の双方における短絡電流の計算を、巻回電極体を展開した形状で行う場合について詳述する。
【００９３】
図１０は、図９に示した巻回電極体における貫通短絡での短絡部の出現箇所を説明する説明図である。図１０では説明のため、巻回電極体１０、正極タブ１１ｂ、及び負極タブ１２ｂ（破線で描画する）の形状及び大きさをデフォルメして記載している。図１０の例は、３巻回、すなわち全６層の巻回セルにおいて、釘などの導電体が巻回電極体を貫通した状態を示している。展開されたモデルは、紙面手前から奥に向けて、正極集電箔、正極活物質、セパレータ、負極活物質、及び負極集電箔が配置された３次元モデルである。正極タブ１１ｂは正極集電箔と同じ面内（奥行き方向で同じ位置）に存在し、展開された電極体１０の帯状の正極集電箔と６箇所で電気的に接続されている。同じように、負極タブ１２ｂは負極集電箔と同じ面内（奥行き方向で同じ位置）に存在し、展開された電極体１０の帯状の負極集電箔と６箇所で電気的に接続されている。正極タブ１１ｂと負極タブ１２ｂは直接には電気的に接続されていない。複数の正極タブ１１ｂ（複数の負極タブ１２ｂ）を、紙面横方向につなぐ部分は、図９に示したようなタブの結束による複数タブの電気的接続を表す仮想的な導電路である。図１０に示すように、貫通短絡の場合、巻回電極体を展開すると、セルの長手方向に複数の短絡部が現れる。本実施の形態では、タブ（正極タブ１１ｂまたは負極タブ１２ｂ）又は隣接電極から各短絡部に向かって流れ込む電流を計算することによって、巻回セルの内部で起こっている短絡現象を適切に表現することができ、巻回セルから外部に現れる熱現象を正確にシミュレートすることができる。巻回電極体には図１０中に示す湾曲部１０ａと呼ばれる電極が曲がった部分があるが、この箇所においては平坦部１０ｂとは異なる物性値（例えば電気伝導率、気孔率、液相伝導率など）を与えてもよい。タブの結束による電気的接続を表す仮想的な導電路についても、計算の収束し易さを考えて適当な物性値（例えば電子伝導率を１．０×１０10Ｓ／ｍなど）を与えてもよい。
【００９４】
図１０に記載した貫通短絡の例は巻回電極体の中央で短絡が生じた場合を示しているが、このシミュレーション手法であれば、短絡位置は巻回電極体の中央に限定されることはなく、例えばタブの近傍や湾曲部１０ａ近くであっても同じモデリング手法で計算が可能である。
【００９５】
図１１は巻回電極体における部分短絡での短絡部の出現箇所を説明する説明図である。図１１では説明のため、巻回電極体１０、正極タブ１１ｂ、及び負極タブ１２ｂ（破線で描画する）の形状及び大きさをデフォルメして記載している。図１１の例は、ニッケル片などの導電体が巻回電極体の一部を短絡した状態を示している。展開されたモデルは、紙面手前から奥に向けて、正極集電箔、正極活物質、セパレータ、負極活物質、及び負極集電箔が配置された３次元モデルである。この場合、巻回電極体を展開すると、セルの長手方向に１箇所の短絡部が現れる。本実施の形態では、タブ（正極タブ１１ｂまたは負極タブ１２ｂ）又は隣接電極から各短絡部に向かって流れ込む電流を計算することによって、巻回セルの内部で起こっている短絡現象を適切に表現することができ、巻回セルから外部に現れる熱現象を正確にシミュレートすることができる。巻回電極体には図１１中に示す湾曲部１０ａと呼ばれる電極が曲がった部分があるが、この箇所においては平坦部１０ｂとは異なる物性値（例えば電気伝導率、気孔率、液相伝導率など）を与えてもよい。タブの結束による電気的接続を表す仮想的な導電路についても、計算の収束し易さを考えて適当な物性値（例えば電子伝導率を１．０×１０10Ｓ／ｍなど。）を与えてもよい。
【００９６】
図１１に記載した部分短絡の例は巻回電極体の中央の最外周層で短絡が生じた場合を示しているが、このシミュレーション手法であれば、短絡位置は巻回電極体の中央に限定されることはなく、最外周層である必要もない。さらには、短絡箇所は複数あってもよい。例えば巻回層２層が短絡した場合、複数個所にニッケル片が混入した場合などであっても、対応する展開形状でシミュレーションを実施することにより、同じモデリング手法で計算が可能である。
【００９７】
以上は短絡電流の計算手法に関する説明であるが、伝熱や発熱反応については実物通りの形状で実施してよい。すなわち、前述のように仮想的に平面に展開する必要はない。
本実施の形態で述べたように、巻回セルの安全性シミュレーションは、巻回電極体を展開した形状での短絡電流のシミュレーションと、実物通りの形状で行う伝熱と発熱反応とのシミュレーションを連成して実施することによって、より実物を的確に表現した計算が実施可能になる。
【００９８】
（実施の形態４）
実施の形態４では、蓄電デバイスの材料分解反応における発熱速度を計算する。
【００９９】
実施の形態４に係るサーバ装置１００は、材料分解反応の発熱反応を記述する式として、Ｄａｈｎモデルとして知られる前述の数５のアレニウス型反応式を用いる。サーバ装置１００の制御部１０１は、数５のアレニウス型反応式を用いて反応率ｘf を計算する。
【０１００】
数５のアレニウス型反応式を利用して反応率ｘf 、発熱密度Ｑを計算する際、反応速度定数ｋ0 、活性化エネルギＥa 、定数ｐ，ｑ，Ｃ0 の値が必要となる。特開２００６−１０６４８号公報に開示された手法では、示差熱分析で得られた温度−発熱量の関係のグラフを、時刻−発熱量の関係に換算してからフィッテングを行い、上述の各パラメータを取得していた。これに対し、本実施の形態では、示差熱分析で得られた温度−発熱量の関係のグラフをそのまま利用して、各パラメータを取得する方法を開示する。
【０１０１】
制御部１０１は、まず、示差熱分析によって得られる温度−発熱量のデータをローレンツ関数又はガウス関数などによってフィッテングし、温度−発熱量の関係式を取得する。ローレンツ関数又はガウス関数によるフィッテングは適宜の最適化ツールを用いて実施される。
【０１０２】
次いで、制御部１０１は、次式で示す微分連鎖律の式を用いて、温度の関数として得られる発熱密度Ｑ（Ｔ）を、時間の関数として表される発熱密度Ｑ（ｔ）に変換する。以下の数６において、温度Ｔの時間偏微分で表される因数は計算の過程で算出される。
【０１０３】
【数９】
2020161115.tif 000011
【０１０４】
以上のように、本実施の形態では、温度−発熱量の関係のグラフをそのまま使用して、時間の関数として表される発熱密度Ｑ（ｔ）を計算するので、温度−発熱量の関係を、時間−発熱量の関係に換算する手間が省ける。
【０１０５】
（実施の形態５）
実施の形態５では、ガス発生をシミュレートする方法について説明する。
【０１０６】
蓄電デバイスにおいて内部短絡などの事象が発生し、安全機構がうまく機能しなかった場合、材料分解反応が進行して蓄電デバイスの内部からガスが噴出する可能性がある。例えば、液系リチウムイオン電池の場合、ガスは昇温により正極活物質中から脱離した酸素が電解液と反応することにより生じる。通常、蓄電デバイスの筐体には圧力で作動する破裂弁などの安全機構が設けられており、筐体内部のガス発生により内圧が上がることで、破裂弁が開き、ガスが外部へ噴出する。
【０１０７】
このガスは高温であるため、隣接セルへの類焼や構造部材の焼損を引き起こす原因となる。また、蓄電デバイスから噴出するガスには、一酸化炭素などの有害なガスが含まれることもあるので、温度、流速、ガス濃度を考慮した上で、周囲の温度を予測することは蓄電デバイスおよび蓄電デバイスを含む製品全体の安全性設計の上で重要である。
【０１０８】
そこで、実施の形態５に係るサーバ装置１００は、材料分解反応による発熱反応に基づき、ガス発生をシミュレートする。ガスの発生に関わる諸量を計算することにより、例えば、熱流体シミュレーションの境界条件として用いることで、電気自動車や発電プラントなどの熱設計および安全性設計に活用することができる。
【０１０９】
ガスの発生に関与する反応として、例えば正極活物質中から脱離した酸素と電解液とによるガス生成や、電極中に含まれる有機助剤の熱分解によるガス生成があることが知られている。これらの反応に、便宜的に１、２、…、ｉ、…と番号を付ける。実用上は、電池のガス発生に関連する素反応過程を詳細に分析考察することは難しい。そこで、示差熱分析（ＤＳＣ）チャートやガス量計測の結果などを照合して、反応を分離する方法を採用することもある。
【０１１０】
制御部１０１は、総ノルマルガス体積Ｖtot を、Ｖtot =Σｖi ｘfiにより計算する。ここでｖi は反応ｉが完全に進行した場合に発生するガスのノルマル体積であり、ｘfiは反応ｉの反応率である。
【０１１１】
制御部１０１は、蓄電デバイス筐体内部の圧力を計算する。筐体内部のガスは理想気体であると仮定する。制御部１０１は、セルの内圧Ｐinを、Ｐin＝Ｐ0 ×（Ｖtot ／Ｖgap ）×（Ｔ／Ｔ0 ）として計算する。ここで、Ｐ0 は初期の蓄電デバイス筐体内部圧力（Ｐａ）であり、通常は１（ａｔｍ）である。Ｖgap は蓄電デバイス筐体内部の気体存在領域の体積（ｍ3 ）、Ｔは気体温度（Ｋ）、Ｔ0 は基準温度（Ｋ）である。
【０１１２】
蓄電デバイスの破裂弁はＰth＜Ｐinの条件を満たす場合に開くように構成される。ここで、Ｐthは破裂弁が開く内圧の閾値である。破裂弁が開いた後は、反応により発生したガスは開いた破裂弁から外部に排出される。
【０１１３】
制御部１０１は、ノルマルガスの発生体積速度ｑnorm,totを、ｑnorm,tot＝Σｖi ｒi により計算する。ここで、ｑnorm,totはノルマルガスの発生体積速度（ｍ3 ／ｓ）、ｖi はｉ番目の反応により生じるノルマルガスの体積（ｍ3 ）、ｒi はｉ番目の反応の反応速度（１／ｓ）であり、ｒi ＝（ｄ／ｄｔ）ｘfi、すなわち反応ｉにおける反応率の時間微分である。
【０１１４】
破裂弁の開口部の面積をＳ（ｍ2 ）としたとき、破裂弁の開口部からの噴出ガス速度は、熱膨張も考慮して、ｖvent＝ｑnorm,tot／Ｓ×（Ｔ／Ｔ0 ）によって計算される。ここで、ｖventは破裂弁の開口部から噴出するガスの噴出ガス速度（ｍ／ｓ）である。この値を、そのまま破裂弁の開口部からのガス噴出の速度境界条件としてもよく、放物線状の速度分布を与えてもよい。流体の計算モデルとしてｋ−εモデルやｋ−ωモデルなどの乱流モデルを使用する場合は、乱流エネルギや乱流消失率などの乱流モデル特有の項に、収束性を考慮した値を与えてもよい。
【０１１５】
制御部１０１は、ガスの噴出速度とともに、発熱反応による発熱速度も計算する。この発熱速度の計算方法は、実施の形態１に記載した方法と同じであってよい。
【０１１６】
発熱反応により発生した熱量は、噴出ガスと蓄電デバイスに適切に按分される。この時の按分比率として、一例としては、ガスと蓄電デバイスが瞬時に熱平衡に至ることを仮定し、各々の熱容量の比で按分する方法を採用してもよい。代替的に、実験結果に整合するように噴出ガスと蓄電デバイスへの熱量の按分比率を与えてもよい。
【０１１７】
開示された実施形態は、全ての点において例示であって、制限的なものではない。本発明の範囲は、特許請求の範囲によって示され、特許請求の範囲と均等の意味及び範囲内での全ての変更が含まれる。
【０１１８】
例えば、実施形態１から５においては単一の蓄電デバイスを例として説明したが、複数の蓄電デバイスからなるシステム（組電池など）についてもシミュレーションを実行可能である。図１２は複数の蓄電デバイスを含むシステム（組電池）を示す模式図であり、図１３は図１２のシステムにおいて複数の蓄電デバイスが連鎖的にガスを排出する様子を可視化したグラフィック表示（動画表示）である。例えば、蓄電デバイスＡが内部短絡などの要因で昇温し、蓄電デバイスＡから排出されたガスが蓄電デバイスＢを昇温させて、蓄電デバイスＢが連鎖的にガスを排出する、というようなガス噴出の連鎖事象をシミュレートすることもできる（図１３参照）。
要約すると、このシミュレーション方法は、
１）第一の蓄電デバイスに関するシミュレーション条件を受付け、受付けたシミュレーション条件に基づき、短絡電流を計算して前記第一の蓄電デバイスから外部への熱現象をシミュレートし、
２）前記第一の蓄電デバイスから外部への熱現象による、第二の蓄電デバイスの加熱に伴う、前記第二の蓄電デバイスから外部への熱現象をシミュレートしている。
他の実施形態として、
１）第一の蓄電デバイスに関するシミュレーション条件を受付け、前記シミュレーション条件は、前記第一の蓄電デバイスを外部から加熱する際の加熱箇所を含み、受付けたシミュレーション条件に基づき、前記第一の蓄電デバイスの加熱に伴う前記第一の蓄電デバイスから外部への熱現象をシミュレートし、
２）前記第一の蓄電デバイスから外部への熱現象による、第二の蓄電デバイスの加熱に伴う、前記第二の蓄電デバイスから外部への熱現象をシミュレートする、シミュレーション方法も考えられる。
これらシミュレーション方法は、シミュレーション装置やコンピュータプログラムとして実現されてもよい。
これらシミュレーション方法により、複数の蓄電デバイスを含むシステムについて、類焼の様子を可視化できる。時系列にどの蓄電デバイスがどのような熱現象を発現するか、どの蓄電デバイスとどの蓄電デバイスが連動して熱現象を発現するかなどを把握できる。
【０１１９】
本明細書において、シミュレーションの対象は蓄電デバイスが外部に及ぼす熱影響に主眼をおいて記載しているが、蓄電デバイスの内部に着目してもよい。例えば、蓄電デバイスから噴出するガスは、その蓄電デバイスの昇温を抑える効果がある。そのため、全く同一の考え方で蓄電デバイス内部の状態を計算することも可能である。
【０１２０】
本明細書において記載した安全性シミュレーションには様々な活用方法が考えられる。例えば、蓄電デバイスを所望の熱条件にさらした場合に、材料分解によるガス発生と内圧上昇により、破裂弁が開くに至るか否かの判定に本願のシミュレーション方法を活用できる。仮に破裂弁が開きガスの噴出が起こったとして、周辺の蓄電デバイスへの類焼が起こるか否かの検討にも本願のシミュレーション方法を活用できる。さらには、例えば、類焼防止のための断熱材や耐火材などの効果を計算で確認することができ、モデルベース開発の理念に基づいた製品設計を強力に推進する手段となる。
【０１２１】
実施形態１ではサーバとクライアント間の通信によりシミュレーションを実施する形態を例示したが、サーバ管理者がシミュレーションプログラムをＤＶＤ−ＲＯＭなどの記憶媒体の手段でクライアント利用者に提供し、クライアント端末ローカルでシミュレーションが実施される形態であってもよい。提供手段として、通信を介したダウンロード形式であってもよい。
【符号の説明】
【０１２２】
１００  サーバ装置
１０１  制御部
１０２  記憶部
１０３  通信部
１０４  操作部
１０５  表示部
２００  クライアント装置
２０１  制御部
２０２  記憶部
２０３  通信部
２０４  操作部
２０５  表示部
Ｎ  通信網
【図１】
2020161115.tif 000012
【図２】
2020161115.tif 000013
【図３】
2020161115.tif 000014
【図４】
2020161115.tif 000015
【図５】
2020161115.tif 000016
【図６】
2020161115.tif 000017
【図７】
2020161115.tif 000018
【図８】
2020161115.tif 000019
【図９】
2020161115.tif 000020
【図１０】
2020161115.tif 000021
【図１１】
2020161115.tif 000022
【図１２】
2020161115.tif 000023
【図１３】
2020161115.tif 000024
【手続補正書】
【提出日】令和2年1月30日(2020.1.30)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートする
シミュレーション方法。
【請求項２】
前記シミュレーション条件は、前記蓄電デバイスにおける内部短絡の発生箇所を含み、
前記内部短絡に伴う前記熱現象をシミュレートする
請求項１に記載のシミュレーション方法。
【請求項３】
前記蓄電デバイスは、巻回電極体を含み、
前記巻回電極体を仮想的に展開した状態にて短絡電流を計算する
請求項１又は請求項２に記載のシミュレーション方法。
【請求項４】
前記シミュレーション条件は、前記蓄電デバイスの外部短絡における抵抗値に関連する情報を含み、
前記外部短絡に伴う前記熱現象をシミュレートする
請求項１に記載のシミュレーション方法。
【請求項５】
前記シミュレーション条件は、前記蓄電デバイスを外部から加熱する際の加熱箇所、熱量、及び環境温度を含み、
前記蓄電デバイスの加熱に伴う前記熱現象をシミュレートする
請求項１から請求項４の何れか１つに記載のシミュレーション方法。
【請求項６】
蓄電デバイスに関するシミュレーション条件を受付け、
前記シミュレーション条件は、前記蓄電デバイスを外部から加熱する際の加熱箇所を含み、
受付けたシミュレーション条件に基づき、前記蓄電デバイスの加熱に伴う前記蓄電デバイスから外部への熱現象をシミュレートする
シミュレーション方法。
【請求項７】
前記蓄電デバイスの示差熱分析によって得られる温度と発熱量との関係を定式化し、
定式化して得られる温度と発熱量との関係式に基づき、前記蓄電デバイスの材料分解反応における発熱速度を計算する
請求項１又は請求項６に記載のシミュレーション方法。
【請求項８】
蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、前記蓄電デバイスの材料分解反応に伴うガスの発生をシミュレートする
シミュレーション方法。
【請求項９】
前記材料分解反応の反応速度に基づき、前記ガスの発生速度、及び熱量の発生速度の少なくとも１つを計算する
請求項８に記載のシミュレーション方法。
【請求項１０】
前記ガスの発生速度は、前記材料分解反応の反応速度に比例するように計算される
請求項９に記載のシミュレーション方法。
【請求項１１】
第一の蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、短絡電流を計算して前記第一の蓄電デバイスから外部への熱現象をシミュレートし、
前記第一の蓄電デバイスから外部への熱現象による、第二の蓄電デバイスの加熱に伴う、前記第二の蓄電デバイスから外部への熱現象をシミュレートする
シミュレーション方法。
【請求項１２】
第一の蓄電デバイスに関するシミュレーション条件を受付け、
前記シミュレーション条件は、前記第一の蓄電デバイスを外部から加熱する際の加熱箇所を含み、
受付けたシミュレーション条件に基づき、前記第一の蓄電デバイスの加熱に伴う前記第一の蓄電デバイスから外部への熱現象をシミュレートし、
前記第一の蓄電デバイスから外部への熱現象による、第二の蓄電デバイスの加熱に伴う、前記第二の蓄電デバイスから外部への熱現象をシミュレートする
シミュレーション方法。
【請求項１３】
外部端末を用いたユーザ認証の後に前記外部端末から送信されるシミュレーション条件を受信し、
受信したシミュレーション条件に基づくシミュレーション結果または前記シミュレーション条件に基づくシミュレーションプログラムを前記外部端末へ送信する
請求項１から請求項１２の何れか１つに記載のシミュレーション方法。
【請求項１４】
蓄電デバイスに関するシミュレーション条件を受付ける受付部と、
受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートするシミュレーション実行部と、
前記シミュレーション実行部によるシミュレーション結果または前記シミュレーション条件に基づくシミュレーションプログラムを出力する出力部と
を備えるシミュレーション装置。
【請求項１５】
コンピュータに、
蓄電デバイスに関するシミュレーション条件を受付けるユーザインタフェースを提示させ、
受付けたシミュレーション条件に基づき、短絡電流を計算して前記蓄電デバイスから外部への熱現象をシミュレートする
処理を実行させるためのコンピュータプログラム。
【手続補正２】
【補正対象書類名】図面
【補正対象項目名】図９
【補正方法】変更
【補正の内容】
【図９】
2020161115.tif 000025
【手続補正書】
【提出日】令和2年3月18日(2020.3.18)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
複数のタブを束ねて電極タブが形成された電極体を備える蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、前記電極タブを介して短絡部に流れ込む短絡電流を計算し、内部短絡に伴う前記蓄電デバイスから外部への熱現象をシミュレートする
処理をコンピュータが実行するシミュレーション方法。
【請求項２】
前記シミュレーション条件は、前記蓄電デバイスにおける内部短絡の発生箇所を含み、
前記コンピュータは、前記内部短絡に伴う前記熱現象をシミュレートする
請求項１に記載のシミュレーション方法。
【請求項３】
前記蓄電デバイスは、巻回電極体を含み、
前記コンピュータは、前記巻回電極体を仮想的に展開した状態にて、前記電極タブを介して前記短絡部に流れ込む短絡電流と、前記巻回電極体の湾曲部から前記短絡部に流れ込む短絡電流とを計算する
請求項１又は請求項２に記載のシミュレーション方法。
【請求項４】
蓄電デバイスの外部短絡における抵抗値に関連する情報を含むシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、短絡電流を計算し、前記外部短絡に伴う前記蓄電デバイスから外部への熱現象をシミュレートする
処理をコンピュータが実行するシミュレーション方法。
【請求項５】
蓄電デバイスに関するシミュレーション条件を受付け、
前記シミュレーション条件は、前記蓄電デバイスを外部から加熱する際の加熱箇所を含み、
受付けたシミュレーション条件に基づき、前記蓄電デバイスの加熱に伴う前記蓄電デバイスから外部への熱現象をシミュレートする
処理をコンピュータが実行するシミュレーション方法。
【請求項６】
蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、前記蓄電デバイスの材料分解反応に伴うガスの発生をシミュレートする
処理をコンピュータが実行するシミュレーション方法。
【請求項７】
前記コンピュータは、前記材料分解反応の反応速度に基づき、前記ガスの発生速度、及び熱量の発生速度の少なくとも１つを計算する
請求項６に記載のシミュレーション方法。
【請求項８】
前記ガスの発生速度は、前記材料分解反応の反応速度に比例するように計算される
請求項７に記載のシミュレーション方法。
【請求項９】
複数のタブを束ねて電極タブが形成された電極体を備える第一の蓄電デバイスに関するシミュレーション条件を受付け、
受付けたシミュレーション条件に基づき、前記電極タブを介して短絡部に流れ込む短絡電流を計算し、内部短絡に伴う前記第一の蓄電デバイスから外部への熱現象をシミュレートし、
前記第一の蓄電デバイスから外部への熱現象による、第二の蓄電デバイスの加熱に伴う、前記第二の蓄電デバイスから外部への熱現象をシミュレートする
処理をコンピュータが実行するシミュレーション方法。
【請求項１０】
第一の蓄電デバイスに関するシミュレーション条件を受付け、
前記シミュレーション条件は、前記第一の蓄電デバイスを外部から加熱する際の加熱箇所を含み、
受付けたシミュレーション条件に基づき、前記第一の蓄電デバイスの加熱に伴う前記第一の蓄電デバイスから外部への熱現象をシミュレートし、
前記第一の蓄電デバイスから外部への熱現象による、第二の蓄電デバイスの加熱に伴う、前記第二の蓄電デバイスから外部への熱現象をシミュレートする
処理をコンピュータが実行するシミュレーション方法。
【請求項１１】
前記コンピュータは、
外部端末を用いたユーザ認証の後に前記外部端末から送信されるシミュレーション条件を受信し、
受信したシミュレーション条件に基づくシミュレーション結果または前記シミュレーション条件に基づくシミュレーションプログラムを前記外部端末へ送信する
請求項１から請求項１０の何れか１つに記載のシミュレーション方法。
【請求項１２】
複数のタブを束ねて電極タブが形成された電極体を備える蓄電デバイスに関するシミュレーション条件を受付ける受付部と、
受付けたシミュレーション条件に基づき、前記電極タブを介して短絡部に流れ込む短絡電流を計算し、内部短絡に伴う前記蓄電デバイスから外部への熱現象をシミュレートするシミュレーション実行部と、
前記シミュレーション実行部によるシミュレーション結果または前記シミュレーション条件に基づくシミュレーションプログラムを出力する出力部と
を備えるシミュレーション装置。
【請求項１３】
コンピュータに、
複数のタブを束ねて電極タブが形成された電極体を備える蓄電デバイスに関するシミュレーション条件を受付けるユーザインタフェースを提示させ、
受付けたシミュレーション条件に基づき、前記電極タブを介して短絡部に流れ込む短絡電流を計算し、内部短絡に伴う前記蓄電デバイスから外部への熱現象をシミュレートする
処理を実行させるためのコンピュータプログラム。
【請求項１４】
コンピュータに、
蓄電デバイスの外部短絡における抵抗値に関連する情報を含むシミュレーション条件を受付けるユーザインタフェースを提示させ、
受付けたシミュレーション条件に基づき、短絡電流を計算し、前記外部短絡に伴う前記蓄電デバイスから外部への熱現象をシミュレートする
処理を実行させるためのコンピュータプログラム。
【請求項１５】
コンピュータに、
蓄電デバイスに関するシミュレーション条件を受付けるユーザインタフェースを提示させ、
受付けたシミュレーション条件に基づき、前記蓄電デバイスの材料分解反応に伴うガスの発生をシミュレートする
処理を実行させるためのコンピュータプログラム。