Li-ion cell Hysteresis Parameter Estimation

Hello,

I am currently using 3-RC ECM state equation to Estimate SoC of a li-ion cell. I intend to incorporate a hysteresis component into this model, described by the following differential equation:

I am familiar with Parameter estimation for a 3-RC ECM viz: R0, R1, R2, R3, Tau1, Tau2 and Tau3. I would like to know how to estimate gamma, M (hysteresis components) using the table based battery.

Could you provide guidance or references on setting up a parameter estimation workflow in MATLAB/Simulink that includes the hysteresis component.

I appreciate your assistance in this matter and look forward to your guidance

Regards

Answer 1

[Javier Gazzarri]

Hello Chinmaya,

In 2025a, Model Based Calibration (MBC) includes a template for ECMs with hysteresis.

Calibrate Powertrain Blockset and Simscape Battery Blocks - MATLAB & Simulink

To do this, you must use the Battery Equivalent Circuit block or, alternatively, MBC's BatteryEstimSOCT_ECB.slx template as shown in the documentation linked above.

It may be a good idea to design an experiment dedicated exclusively to this end, such as a low C-rate back and forth charge and discharge, as explained in Prof. Plett's paper: (ECE5710-Notes02.pdf, Section 2.4)

Please let me know if you have any questions.

Best regards,

Javier

[Chinmaya]

Hello Sir,

Thank you for your detailed response.

I had a couple of follow-up questions regarding the experiment. Would performing the charge-discharge cycle at C/2 be acceptable for capturing hysteresis characteristics, or would you recommend using a lower C-rate such as C/3 or C/5 for better accuracy?

Additionally, I’ve attached a sample charge-discharge profile for your reference. Could you kindly confirm if this approach aligns with the recommended test procedure?

Looking forward to your guidance.

Regards,

Chinmaya

Answer 2

[Javier Gazzarri]

Hello Chinmaya,

I'd recommend C/20 or lower. C/5 is still too high to avoid ohmic losses.

Your profile goes back to 100% in the charging direction. I'd recommend the following SOC breakpoints:

100 - 0 - 90 - 10 - 80 - 20 - 70 - 30 - 60 - 40 - 50

In this way you will capture the hysteresis behavior in both directions at all SOCs.

Best regards,

Javier