李百華1，郭燦彬1,鐘其水2，涂濤2

　?。?.廣東機電職業技術學院，廣東 廣州 510515；2.電子科技大學 航空航天學院，四川 成都 611731）

       摘要：首先分析了鋰離子電池的特性和充放電原理，介紹了鋰離子電池的戴維南（Thevinin）等效電路模型，并采用脈沖法和遞推最小二乘法相結合對戴維南等效電路模型參數進行辨識，實現了戴維南等效電路模型參數的實時在線辨識，為電動汽車鋰電池的等效電路模型的分析研究提供了一種可行方式。

　　關鍵詞：電動汽車；等效電路；脈沖法；最小二乘法

　　中圖分類號：U469.72+2文獻標識碼：ADOI： 10.19358/j.issn.1674-7720.2017.01.025

　　引用格式：李百華，郭燦彬，鐘其水，等. 電動汽車鋰電池戴維南等效電路模型參數辨識研究［J］.微型機與應用，2017,36（1）：83-85，88.

0引言

　　隨著汽車需求不斷上升，能源短缺與環境污染日趨嚴重，電動汽車具有無污染、電能來源多樣的優點，其發展得到廣泛重視。電動汽車鋰電池存在過充或過放后將產生不可恢復的破壞性影響，不正常工作狀態下將會導致自燃等安全問題［1］。通過建立有效理論模型對其工作狀態、特性和性能作有效評價是必要的。

　　目前，常用的電池模型有多種，戴維南等效電路模型由于模型參數簡單、計算量適中得到廣泛采用。本文通過分析戴維南等效電路模型，研究等效電路模型參數辨識，得到該等效電路模型參數的確定方法［2］。

1鋰離子電池充放電原理

　　鋰離子電池正極一般采用鈷酸鋰、錳酸鋰、磷酸鐵鋰、鎳酸鋰等鋰離子化合物，負極一般采用石墨等可以嵌入鋰離子的化合物［3］。磷酸鐵鋰離子電池正極充放電變化如圖1所示，鋰離子通過電解液和電池隔膜在電池的正負極間來回轉移。

2戴維南（Thevinin）等效電路模型

　　電池組等效為電壓源、內阻、一階RC串聯的電路，如圖2所示，其中R1是電池的歐姆內阻，一階RC表征電池充放電過程中極板的極化電容和電阻［4］。這樣電池組充放電過程中電壓的突變特性可由內阻R1表征，電壓的漸變特性可由Rp和Cp表征。

　　戴維南模型中，歐姆內阻R1由電池電解質、正極鋁箔、負極銅箔等組成，用一階的RC串聯電路反應電池充放電時電壓的漸變特性，一階RC電路由電池充放電過程中正極和負極的極化效應產生。

3戴維南等效電路模型參數辨識

　　3.1遞推最小二乘法

　　對于多元線性方程，在未知數個數小于方程個數，方程無解的情況下［5］，可以通過最小二乘法計算方程的最小二乘解。M抽頭濾波器的權向量應滿足的線性方程如式(1)所示。

　　Aw=b(1)

　　式中A是列滿秩矩陣，且A的行數大于列數，即式(1)中，獨立方程數大于未知數的數目，因此方程無解。但是，應用最小二乘法可以計算方程的一個近似解，通過取得極小值并求導后可以得到式(2)：

　　LS=(AHA)－1AHb(2)

　　式(2)即為方程的最小二乘解。將最小二乘解推廣到遞推最小二乘(RLS)解需要假定系數矩陣A如式(3)所示。

　　其中M為濾波器系數的個數或者線性方程的未知數個數，N－M+1為觀測數據的個數，或者線性方程的個數。為了充分利用觀測數據，可以將矩陣A擴展為式(4)所示的結果。

　　式(1)中相關矩陣如式(5)所示:

　　3.2復頻域分析

　　鋰電池的戴維南模型等效電路如圖2所示，將該模型的電路數學方程轉換到復頻域可得電路輸出方程如下所示:

　　式(9)中，U(s)為電池端電壓復頻域值，I(s)為電流復頻域值，R1為電池內阻值，Rp和Cp為電池極化電阻和極化電容。同時得到電池阻抗如下所示:

　　式(10)中， E(s)=Uoc(s)－U(s)，根據雙線性變換原理可得:

　　將式(11)帶入式(10)可得z域下鋰電池戴維南模型參數表達式如下所示:

　　從而可得：

　　E(k)=k1I(k)+k2I(k－1)+k3E(k－1)(14)

　　3.3脈沖放電實驗和遞推最小二乘法結合辨識戴維南模型參數

　　由式(13)可知，參數k3只與采樣周期和時間常數有關，時間常數可由脈沖放電實驗獲得，脈沖放電實驗的系統如圖3所示。

　　圖4是磷酸鐵鋰電池以10 A電流放電的端電壓變化曲線。

　　圖4中電池端電壓的漸變過程是由電池的極化電阻和極化電容引起的，通過MATLAB進行曲線擬合可得電壓漸變段電池端電壓與時間的方程關系，如下所示：

　　U=0.051 9exp(-0.092 5t)+3.269 0(15)

　　由式(15)可知，電壓漸變曲線的時間常數如下式所示:

　　式(19)中，b(n)=E(n)－0.687 76E(n－1)，φ(n)=［I(n)I(n－1)］，θ(n)=［k1k2］T。b(k)的值與E(k)相關，而E(k)的值與電池的開路電壓值Uoc(k)相關，電池的開路電壓值Uoc(k)可由電池剩余電量SOC(k)求得，因此，這里需要獲得電池SOC值和開路電壓值Uoc間的擬合曲線。

　　為了獲取鋰電池SOC值和開路電壓Uoc間的擬合曲線關系，可以對鋰電池進行雙脈沖放電實驗，獲取指定SOC值下的開路電壓值，從而進行擬合［6］。鋰電池雙脈沖放電的端電壓曲線如圖5所示，實驗中，充放電脈沖電流大小都是10 A，脈沖寬度為40%。

　　圖5中，a點是電池脈沖放電結束后，靜置了30 s的端電壓值，在靜置過程中，電池極化電壓將逐漸降低，電池端電壓值逐漸趨近電池開路電壓值。b點是電池充電結束后，靜置了30 s的端電壓值，在靜置過程中，電池充電極化電壓逐漸降低，電池端電壓值逐漸趨近電池開路電壓值。由于電池脈沖充放電的時間很短，因此，電池的SOC可以認為不變,a點和b點分別是電池放電后靜置一段時間和充電后靜置一段時間的端電壓值，所以，電池的開路電壓值一定介于a點和b點電壓值之間，可以取a、b點電壓的平均值作為電池的開路電壓值。通過雙脈沖法測得不同SOC值的開路電壓值如表1所示。

　　根據表1，通過分段擬合可得磷酸鐵鋰電池SOC和開路電壓間的擬合曲線如圖6所示。

　　圖中a~b段采用5次擬合，b~c段采用3次擬合，c~d段采用二次擬合，d~e段采用5次擬合。通過分段擬合的方式獲取鋰離子電池SOC值和開路電壓值間的曲線關系，可以大大提高通過電池SOC值獲取開路電壓值的精度。

　　電流I(k)可以通過實驗測得，令λ=1,u(n)=I(n)，帶入式(8),可迭代計算k1、k2的值，將k1、k2的值代入式(17)即可求得R1、Rp的迭代值。通過實驗得到R1、Rp的估算值如圖7所示。

4結論

　　本文首先分析了鋰離子電池的特性和充放電原理，介紹了鋰離子電池戴維南等效電路模型，并采用脈沖法和遞推最小二乘法相結合對戴維南等效電路模型參數進行辨識。通過雙脈沖法測得了不同SOC值下磷酸鐵鋰電池的開路電壓值，并對磷酸鐵鋰電池的SOC和開路電壓值進行了分段擬合。通過脈沖放電法和遞推最小二乘法實現了戴維南等效電路模型參數的實時在線辨識。該方法為電動汽車鋰電池的等效電路模型的分析研究提供了一種可行方式。

參考文獻

　?。?］ 張禹軒. 電動汽車動力電池模型參數在線辨識及 SOC 估計［D］.長春：吉林大學，2014.

　?。?］ 周美蘭, 趙強, 周永勤. 改進的PSOBP神經網絡估算磷酸鐵鋰電池SOC［J］. 哈爾濱理工大學學報, 2015, 20(4):88-92.

　?。?］ GREGORY L P. LiPB dynamic cell models for Kalmanfilter SOC estimation［A］. Proceedings of the 19th International Electric Vehicle Symposium［C］. Busan:EVS,2002：1860-1871.

　　［4］ 張東華,馬燕，陳思琪，等.鋰電池模型參數估計與荷電狀態估算研究［J］.武漢理工大學學報,2015,37(2):179-182.

　?。?］ 范興明, 曾求勇,張鑫. 基于改進安時積分法的電動汽車電池SOC估計與仿真研究［J］.電氣應用, 2015(8):111-115.

　?。?］ Shang Yunlong, Zhang Chenghui, Cui Naxin, et al. A celltocell battery equalizer with zerocurrentswitching and zerovoltage gap based on quasiresonant LC converter and boost converter［C］. IEEE Transactions on Power Electronics, 2015, 30(7):3731-3747.