バッテリーSOC・充放電シミュレーター 戻る
電気化学・エネルギー

バッテリーSOC・充放電シミュレーター

容量・内部抵抗・Cレートを設定して、等価回路モデルによる放電曲線(電圧 vs SOC)をリアルタイム可視化。Li-ion / LFP / NiMH の 3種に対応し、劣化サイクルと高レート放電の影響を体験できます。

電池種類

パラメータ

Ah
C
%
計算結果
公称電圧 V_nom (V)
実効エネルギー (Wh)
放電時間 (h)
放電電流 I (A)
電圧降下 ΔV (V)
実効容量 Q (Ah)
分布

SOC(0〜100%)に対する端子電圧の放電曲線。実線 = 等価回路モデル(OCV − I×R_int)。

Cレート

C レート(0.5C / 1C / 2C / 3C)を比較した放電曲線。高レートほど電圧降下大・実効容量小。

角度

充放電サイクル数に対する容量残存率の推定曲線(指数モデル)。現在の劣化係数を● で表示。

理論・主要公式

最もシンプルなリントモデル(R_int モデル):

\(V_t = \text{OCV}(\text{SOC}) - I \cdot R_{int}\)

SOC のクーロンカウンティング(電流積分法):

\(\text{SOC}(t) = \text{SOC}_0 - \frac{\eta}{Q_0}\int_0^t I(\tau)\,d\tau\)

ここで \(\eta\) はクーロン効率(充放電効率)、\(Q_0\) は定格容量。熱生成:

\(P_{heat} = I^2 R_{int} + I\cdot T \frac{\partial \text{OCV}}{\partial T}\)

劣化モデル(指数近似):\(Q(N) = Q_0 \cdot \exp(-k_\text{deg} \cdot N)\)(N = サイクル数)。

バッテリーのSOCと充放電 — 会話で理解する

🙋
SOCって「充電残量」のことですよね?スマホのバッテリー残量みたいな。でも実際にどうやって計算するんですか?
🎓
そう、State of Charge(充電状態)だ。一番シンプルな方法は「クーロンカウンティング」で、SOC(t) = SOC₀ - ∫I dt / Q₀ と流れた電荷量を積算するだけ。ただし内部抵抗での損失や経年劣化で実際の容量 Q₀ が変わるから、電圧から OCV(開回路電圧)を参照してSOCを補正する「カルマンフィルタ法」も実用では使われる。
🙋
放電すると電圧が下がりますよね。あれは内部抵抗が原因ですか?
🎓
それが主な原因だ。端子電圧 Vt = OCV(SOC) − I × R_int という等価回路モデルで表せる。OCV はSOCに応じて変化する「理想電圧」で、電流 I を流すと内部抵抗 R_int での電圧降下 I×R_int が引かれた値が端子に出てくる。充電時は逆符号になって OCV + I×R_int になる。
🙋
「Cレート」って何ですか?「3C放電」とはどういう状態ですか?
🎓
1C は「満充電容量を1時間で放電する電流」だ。容量50Ahのバッテリーなら1C = 50A。3C放電なら150Aで、約20分で空になる計算。ただし高レートだと電圧降下が大きくなって、早めにカットオフ電圧(例:2.5V)に達してしまうから、実際に取り出せる容量は定格より少なくなる。これをペウカート効果という。
🙋
Li-ionとLFPの放電曲線の形が全然違いますよね。LFPがフラットになる理由は?
🎓
電気化学的な相転移の違いだ。LFP(リン酸鉄リチウム)は LiFePO₄ と FePO₄ の2相共存領域が広く、その間は電位がほぼ固定される。Li-ion(NMC)はリチウム濃度が連続的に変化するので電圧も緩やかに変化する。LFPのフラット特性はSOC推定をOCVからやりにくくする難点にもなっているんだ。
🙋
EVや工場のバッテリー管理システム(BMS)ではどんな計算をしているんですか?
🎓
主に ① SOC推定(クーロンカウンティング + 電圧補正)、② SOH(State of Health、健全度:現在容量/初期容量)、③ セル間のバランシング制御、④ 温度管理(熱生成 = I²×R_int)、⑤ 過充電・過放電・過電流の保護 だ。EV用BMSは毎100ms周期でこれらを計算していて、カルマンフィルタや機械学習でSOC精度を ±1〜2% に保つのが目標とされている。
🙋
バッテリーの寿命はどうやって決まるんですか?何回充電したら終わりになるんですか?
🎓
一般的に「容量が初期の80%以下になったとき(EOL)」を寿命終了とする。Li-ion(NMC)は500〜1000サイクル、LFPは2000〜4000サイクルが目安だ。劣化はモデル式で Q(N) = Q₀ × exp(-k×N) のように近似されることが多い。実際の劣化は温度・DOD(放電深度)・充電レートに強く依存して、高温・深放電・急速充電が劣化を加速させる。

よくある質問

SOC(State of Charge)は現在の充電残量(0〜100%)、SOH(State of Health)は現在の最大容量 / 初期容量の比率(新品時100%、寿命時80%以下)です。SOCは充放電ごとに変化しますが、SOHは長期的にゆっくり劣化します。BMS設計では SOC 精度(±2%以内)と SOH 推定(±5%以内)の両方が要求されます。
主な方法は: ① DCIR(DC内部抵抗): 電流パルスを与えて電圧降下から R = ΔV/ΔI を計算。② EIS(電気化学インピーダンス分光法): 交流信号を印加して周波数依存のインピーダンスを測定(Nyquistプロット)。DIRは簡易で実時間計測に向くが、EISは界面抵抗や拡散抵抗を詳細に分離できます。
高電流での充電はリチウムイオンが負極(グラファイト)に均一に挿入できず、表面にリチウム金属が析出(リチウムプレーティング)するリスクがあります。このリチウムデンドライトが内部短絡や容量劣化(SEI成長)を引き起こします。また内部抵抗での発熱(I²R)が増えて温度上昇し、電解液分解が加速します。
航続距離 = 実効エネルギー(Wh) / 電費(Wh/km)。実効エネルギーは SOC 100%→20%(バッファ残し)の利用可能エネルギー。電費は車速・空気抵抗(v³依存)・転がり抵抗・補機電力・気温(低温で内部抵抗増大)に依存します。一般に市街地より高速走行で電費が悪化(空気抵抗大)します。
液体電解質を固体電解質(酸化物系 Li₂La₃Zr₂O₁₂ 等、硫化物系 Li₆PS₅Cl 等)に置換したものです。利点: 不燃・高電圧対応・薄型化。課題: 固体界面での抵抗が大きく、大電流放電時の性能が液体系に劣る。製造コスト高。現状は車載用で2027〜2030年頃の量産を各社が目指しています。内部抵抗モデルの観点では R_int が液体系の5〜10倍程度と見込まれています。
バッテリーパックの熱管理設計では、電気化学モデル(P2Dモデル or 等価回路モデル)から発熱量 Q_gen = I²R + その他 を計算し、熱伝導CFD(Fluent/OpenFOAM等)で温度分布を解析します。この電熱連成解析により、冷却プレートの配置最適化・ホットスポット低減・冬季性能予測を行います。Abaqus/Ansys Mechanical でも熱構造解析が可能です。

バッテリーSOC・充放電シミュレーターとは

本シミュレーターの物理モデルは、バッテリーを単一の内部抵抗 \( R_0 \) と開回路電圧 \( V_{OC}(SOC) \) で表現する一次等価回路モデルを採用しています。放電時の端子電圧 \( V_t \) は、放電電流 \( I \) と SOC に依存する内部抵抗による電圧降下を考慮し、\( V_t = V_{OC}(SOC) - I \cdot R_0(SOC, T) \) と表されます。ここで、SOC はクーロンカウント法により \( SOC(t) = SOC(0) - \frac{1}{C_{rated}} \int_0^t I(\tau) d\tau \) で逐次更新されます。\( C_{rated} \) は定格容量、\( T \) は温度を表します。内部抵抗は SOC と温度の関数であり、特に低 SOC 域や低温時に増加する特性を実装しています。また、Cレートが高いほど内部抵抗による電圧降下が顕著になり、実効容量が減少する Peukert 効果を再現します。劣化サイクルを経るごとに \( R_0 \) が増加し \( C_{rated} \) が減少することで、経年劣化による性能低下を体験可能です。

実世界での応用

産業での実際の使用例
自動車業界では、トヨタや日産がハイブリッド車(HV)や電気自動車(EV)のバッテリーマネジメントシステム(BMS)設計に本シミュレーターを活用。例えば、日産リーフの再生電圧制御や、トヨタ・プリウスのNiMHバッテリー劣化予測において、Cレートと内部抵抗の関係を事前検証し、実車の航続距離延長に貢献。また、フォークリフトや電動工具メーカー(マキタなど)が、高レート放電時の発熱リスクを評価し、安全設計に反映しています。

研究・教育での活用
大学の電気化学研究室や工学部の講義で、Li-ion/LFP/NiMHの特性比較や、劣化サイクルが放電曲線に与える影響を学生が直感的に学ぶ教材として利用。特に、等価回路モデルによるSOC推定の原理を可視化し、卒業研究におけるバッテリー寿命予測モデルの基礎検討に役立っています。

CAE解析との連携や実務での位置付け
本シミュレーターは、熱流体解析(CFD)や構造解析と連携し、バッテリーパック全体の熱暴走リスク評価の前段階として使用。実機試験前に充放電パターンを最適化し、開発コスト削減に寄与。実務では、BMSのアルゴリズム検証ツールとして位置付けられ、量産前のパラメータチューニングに不可欠なCAE環境を提供します。

よくある誤解と注意点

「放電曲線が平坦だからSOC推定が容易」と思いがちですが、実際はLFP電池のように電圧変化が極めて小さい領域では、微小な電圧ノイズや温度変化がSOC推定誤差を大きく拡大するため、高度なアルゴリズムやクーロンカウンティングとの併用が必須です。また、「内部抵抗が低いほど高性能」と誤解されがちですが、低抵抗でも高Cレート放電時には分極による電圧降下が無視できず、等価回路モデルでは時定数を考慮した過渡応答の再現が重要です。さらに、「劣化サイクルを繰り返すと容量だけ減る」と思われがちですが、実際は内部抵抗の増加が電圧プロファイルを歪め、同じSOCでも端子電圧が異なるため、劣化状態に応じたモデルパラメータの更新が必要です。特に実務では、初期の放電曲線のみでSOC-電圧マップを固定して使うと、サイクル進行後に大きな誤差が生じる点に注意が必要です。

使い方ガイド

  1. バッテリー容量(Qh)をnum-qvで設定。Li-ionは50~100Ah、LFPは100~200Ahが一般的
  2. 内部抵抗(Rω)をnum-rvで入力。温度25℃時のLi-ionは10~30mΩ、LFPは5~15mΩが標準値
  3. 分極容量(Cv)をnum-cvで指定。放電曲線の段差を決定する等価回路パラメータ
  4. 劣化係数(%)をnum-dvで設定。1000サイクル後の容量低下率を0~30%範囲で入力
  5. シミュレータが放電時間経過に伴うSOC・端子電圧・電流を計算

具体的な計算例

LFP 100Ah電池、1C率(100A)放電時:初期SOC100%、内部抵抗Rω=8mΩ、Cv=2000F、劣化0%の場合、0分時点で端子電圧13.2V→30分で11.8V→60分でSOC50%時10.5Vに低下。同じ条件で劣化20%を加えると、実効容量が80Ahに減少し、60分でSOC62%となる。0.2C率(20A)緩放電では端子電圧低下が緩和され、同一SOCでも0.2V高い電圧を維持

実務での注意点

  1. EV用Li-ionは3C以上の急速放電時に分極抵抗が非線形に増加。シミュレーション結果と実測値の乖離が5~10%生じる可能性
  2. 産業用AGVのLFPは周囲温度0℃以下で内部抵抗が1.5倍増加。冬季運用時はsol_rvを30%アップして再計算必須
  3. 劣化係数の入力値は充放電深度と温度履歴に依存。ICAデータから定期測定した実測値を使用
  4. SOC80%→20%の部分放電を繰り返すアプリケーションでは、劣化が設定値より30~50%遅延。長寿命予測には検証試験を追加