TWI760246B - 一種可縮短充電時間之多階段定電流充電方法 - Google Patents
一種可縮短充電時間之多階段定電流充電方法 Download PDFInfo
- Publication number
- TWI760246B TWI760246B TW110120920A TW110120920A TWI760246B TW I760246 B TWI760246 B TW I760246B TW 110120920 A TW110120920 A TW 110120920A TW 110120920 A TW110120920 A TW 110120920A TW I760246 B TWI760246 B TW I760246B
- Authority
- TW
- Taiwan
- Prior art keywords
- charging
- current
- battery
- remaining power
- stage
- Prior art date
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
- Y02E60/10—Energy storage using batteries
Landscapes
- Charge And Discharge Circuits For Batteries Or The Like (AREA)
- Secondary Cells (AREA)
Abstract
一種可縮短充電時間之多階段定電流充電方法,其包括以下步驟: 在一電池之初始剩餘電量和100%電量之間設置四個目標剩餘電量以定義出六個剩餘電量值,並依所述六個剩餘電量值定義出五個階段,其中各所述階段之期初剩餘電量和期末剩餘電量係由兩個相鄰的所述剩餘電量值決定;在所述五個階段中對應設定五個充電電流;依所述六個剩餘電量值及所述五個充電電流執行五次充電時間計算以產生五個充電時間,其中各所述充電時間計算表為:(該電池之額定容量)*( 兩個相鄰的所述剩餘電量值之絕對差值)/(一對應之所述充電電流);以及依所述五個充電電流在對應之所述五個充電時間內對該電池充電。
Description
本發明係關於二次電池充電之演算法,特別是一種基於多階段定電流充電法之二次電池充電之演算法。
隨著科技進步,3C產品在積體電路快速發展之影響下,逐漸走向小型化且具有更強大之工作性能,卻也造成電量消耗困擾與電力來源供給問題。如今多數可攜式產品,如手機、筆記型電腦、行動電源等,其電力來源均為可重複充電的二次電池,近年來由於各國禁用燃油車之宣示,電動車產業的發展,使得二次電池的製造材料技術、充電平衡技術及電能轉換與管理技術發展備受矚目。
在二次電池當中,鋰離子電池因為具有高能量密度、循環壽命長、體積輕巧、充放電速度快、無記憶效應、平均工作電壓高及自放電率低等特點,所以能成為二次電池中的主流,許多大電力儲能系統也開始投入大容量鋰離子電池的研究與發展。
相對於上述優點,鋰離子電池仍具有不耐過充、過放電及易受溫度影響等缺點,故需針對鋰離子電池之潛在問題進行考量,如何達到縮短充電時間、降低充電溫升、提高充電效率及延長電池的使用時間等目標,是鋰離子電池充電技術的討論重點。
二次電池為求輕、薄致使電池容量受到限制,因此快速充電技術變得十分重要,然而為了追求快速充電必須將充電電流加大,這將使二次電池的溫度提升、電池壽命減短。溫升對電池來說非常重要,充電時溫升過高除了有發生爆炸之安全顧慮以外,也會造成電池加速老化、電池組之間電量不平衡、容量衰退、和循環壽命減短等問題。上述老化之現象係因為在較高的操作電流下,電化學反應過度激烈,引起電池內阻逐漸變大所致。
因此,充電技術對二次電池來說十分重要。充電技術關係到二次
電池之充電速度、充電效率、電池溫升值,電池循環壽命等因素,習知技術之充電法最普遍被使用係定電流定電壓充電法(Constant Current Constant Voltage,CC-CV),該充電法分為兩階段,先以定電流對電池充電直到電池端電壓達到額定截止端電壓後,再使用截止端電壓對電池進行充電,此時充電電流會因端電壓與內電壓之差逐漸下降,當電流下降至額定截止電流時即可視為充飽,其優點為簡單容易實現,缺點為充電時間較久。
因此許多文獻提出改良之充電方法。有文獻利用雙迴路控制,使用正、負迴授控制電池電壓,可以得到與定電流-定電壓相似的充電曲線,且不需使用電流感測器,故較定電流定電壓充電法簡單且成本更低;另有文獻提出一開始使用高於電池上限額定電壓(如4.3V)進行充電之升壓式定電流-定電壓充電法,在定電壓週期過後切換至定電流定電壓充電法,此法能在短時間將電池充至額定容量的30%,所以充電時間較短,但缺點為充電前電池必須完全放電;亦有文獻提出主動式充電狀態結合模糊控制充電法及灰預測控制電池充電系統,可在相同樣時間內充入較多的電量。
再有文獻提出以智慧演算法決定多階段定電流充電法各階段之充電電流;或有文獻提出以電池模型為基礎之最佳化多階段充電法,其優點為用簡單數學運算決定充電時間為最短的充電電流,缺點為無法優化其他因素如充電損失、充電效率等;又有文獻透過改變電流大小和脈衝寬度及脈衝之間的休息時間使充電的方法有多種變化;尚有文獻提出模糊五階段鋰離子電池充電法,並使用田口方法來決定模糊控制之歸屬函數,藉此優化充電方法之表現。
然而上述文獻皆需進行多次實際實驗以找出最佳化電流設定,因此本領域亟需一新穎的可縮短充電時間之多階段定電流充電方法。
本發明之一目的在於揭露一種可縮短充電時間之多階段定電流充電方法,其係透過交流阻抗分析以建立鋰離子電池之等效模型,藉由庫倫積分與等效電路模型求得充電時間與充電損失,並利用粒子群演算法尋找基於剩餘容量之多階段定電流最佳化充電設定,藉此實現縮短充電時間及降低充電損
失之最佳化多階段充電法。
本發明之另一目的在於揭露一種可縮短充電時間之多階段定電流充電法,係利用粒子群演算法尋找最佳化充電電流設定值,僅需使用電池等效模型,不需制訂田口直交表之每階段初始電流值之設置及多次耗時費工之實際實驗即能得到最佳化電流設定。
本發明之再一目的在於揭露一種可縮短充電時間之多階段定電流充電法,相較習知技術之CC-CV充電法可將電池容量充至CC-CV所得容量的97%以上,其充電時間能縮短14.12%且平均溫升低0.69℃。
為達前述目的,一種可縮短充電時間之多階段定電流充電方法乃被提出,其包括以下步驟:在一電池之初始剩餘電量和100%電量之間設置四個目標剩餘電量以定義出六個剩餘電量值,並依所述六個剩餘電量值定義出五個階段,其中各所述階段之期初剩餘電量和期末剩餘電量係由兩個相鄰的所述剩餘電量值決定;在所述五個階段中對應設定五個充電電流;依所述六個剩餘電量值及所述五個充電電流執行五次充電時間計算以產生五個充電時間,其中各所述充電時間計算表為:(該電池之額定容量)*(兩個相鄰的所述剩餘電量值之絕對差值)/(一對應之所述充電電流);以及依所述五個充電電流在對應之所述五個充電時間內對該電池充電。
在一實施例中,所述四個目標剩餘電量及五個充電電流係依一粒子群演算法決定。
在一實施例中,該粒子群演算法包括:以所述四個目標剩餘電量及五個充電電流代表一粒子空間;在一記憶體中儲存該粒子空間中之n個不同粒子的目前位置、目前速度、目前適應值及個別最佳位置,以及所述n個不同粒子之一群體最佳位置及一群體最佳適應值,n為大於1之整數;依所述目前位置、目前速度、個別最佳位置及群體最佳位置對各所述粒子進行一位置更新程序以產生一新位置以更新所述目前位置,依該新位置之內含資訊進行一適應值計算程序以獲得一適應值,及在該適應值小於該目前適應值時以該適應值更新該目前適應值;以及在n個所述目前適應值中之最小者小於該群體目前最佳
適應值時,以該最小者更新該群體最佳適應值,並以與該最小者對應之所述目前位置更新該群體最佳位置。
在一實施例中,該位置更新程序包括以下公式:
x ij (n+1)=x ij (n)+v ij (n+1),其中,w為權重值(介於0~1之間),C 1和C 2分別為個體和群體學習因子,均介於1~4之間,x ij 代表粒子的位置,i為第幾個粒子,j為其粒子維度,v ij 表第i粒子j維度的速度,pBest ij 表個體粒子最佳的位置,gBest ij 表群體最佳的位置,r 1和r 2為介於[0,1]之間均勻分布亂數,n為疊代次數。
其中,s為階段數,I charge,s 為第s階段之充電電流,SOC charge,s 為第s階段之充入電量,SOC 0 為起始充電之剩餘容量,Q rate 為電池之額定容量,R s 為串聯等效電阻,R p 為並聯等效電阻,Objective為該適應值,α為權重值,(T now ,L now )為一目前總充電時間及一目前充電總損失,(T min ,L min )為一理想總充電時間及一理想充電總損失,(T max ,L max )為一上限總充電時間及一上限充電總損失。
在一實施例中,該權重值α為0.5。
其中,R eq 為等效阻抗,SOC為電池剩餘容量。
在一實施例中,該預設之額定電壓為4.2V。
為使 貴審查委員能進一步瞭解本發明之結構、特徵及其目的,茲附以圖式及較佳具體實施例之詳細說明如後。
步驟a:在一電池之初始剩餘電量和100%電量之間設置四個目標剩餘電量以定義出六個剩餘電量值,並依所述六個剩餘電量值定義出五個階段,其中各所述階段之期初剩餘電量和期末剩餘電量係由兩個相鄰的所述剩餘電量值決定
步驟b:在所述五個階段中對應設定五個充電電流
步驟c:依所述六個剩餘電量值及所述五個充電電流執行五次充電時間計算以產生五個充電時間,其中各所述充電時間計算表為:(該電池之額定容量)*(兩個相鄰的所述剩餘電量值之絕對差值)/(一對應之所述充電電流)
步驟d:依所述五個充電電流在對應之所述五個充電時間內對該電池充電
圖1繪示本發明之可縮短充電時間之多階段定電流充電方法之流程圖。
圖2繪示習知技術之多階段定電流充電法之充電示意圖。
圖3a繪示戴維寧電池等效模之示意圖。
圖3b繪示圖3a的改良電池等效電路模型。
圖4繪示電池之交流阻抗特性示意之奈氏圖。
圖5a繪示交流阻抗之分析流程圖。
圖5b繪示交流阻抗測試之實驗流程圖。
圖6a繪示EC-Lab交流阻抗分析之操作介面示意圖。
圖6b繪示交流阻抗R s 、R p 與等效阻抗R eq 對剩餘容量百分比之關係曲線圖(環境溫度5℃)。
圖6c繪示交流阻抗R s 、R p 與等效阻抗R eq 對剩餘容量百分比之關係曲線圖(環境溫度25℃)。
圖6d繪示交流阻抗R s 、R p 與等效阻抗R eq 對剩餘容量百分比之關係曲線圖(環境溫度45℃)。
圖7繪示粒子群演算法之粒子在解空間移動之示意圖。
圖8a繪示本發明之電池剩餘容量與等效阻抗之擬合曲線圖。
圖8b繪示量測阻抗與擬合曲線之誤差量。
圖8c繪示0.5C CC-CV充電法之模擬結果。
圖8d繪示2C CC-CV充電法之模擬結果。
圖9a繪示本發明採用之充電機系統架構示意圖。
圖9b繪示五階段充電法之操作流程圖。
圖10a繪示本發明進行充電實測之實驗平台示意圖。
圖10b繪示圖10a之五階段充電法之操作流程圖。
圖10c繪示習知技術之CC-CV實測0.5C(1.3 A)之充電電壓、電流波形圖。
圖10d繪示習知技術之CC-CV實測1C(2.6 A)之充電電壓、電流波形圖。
圖10e繪示習知技術之CC-CV實測2C(5.2 A)之充電電壓、電流波形圖。
圖11a繪示Case1(本發明充電法)之實測波形圖
圖11b繪示Case2(習知技術一)之實測波形圖。
圖11c繪示Case3(習知技術二)之實測波形圖。
圖11d繪示Case4(習知技術三)之實測波形圖。
圖11e繪示Case5(微調剩餘容量轉態點)之實測波形圖。
圖11f繪示Case6(微調各階段之電流值)之實測波形圖。
圖11g繪示0.5C(1.3 A)、1C(2.6 A)、2C(5.2 A)、Case1至Case6之充電溫升曲線圖。
圖11h繪示Case1至Case6之充電溫升曲線圖。
請參照圖1,其繪示本發明之可縮短充電時間之多階段定電流充電方法之流程圖。
如圖1所示,本發明之可縮短充電時間之多階段定電流充電方法包括以下步驟:在一電池之初始剩餘電量和100%電量之間設置四個目標剩餘電量以定義出六個剩餘電量值,並依所述六個剩餘電量值定義出五個階段,其中各所述階段之期初剩餘電量和期末剩餘電量係由兩個相鄰的所述剩餘電量值決定(步驟a);在所述五個階段中對應設定五個充電電流(步驟b);依所述六個剩餘電量值及所述五個充電電流執行五次充電時間計算以產生五個充電時
間,其中各所述充電時間計算表為:(該電池之額定容量)*(兩個相鄰的所述剩餘電量值之絕對差值)/(一對應之所述充電電流)(步驟c);以及依所述五個充電電流在對應之所述五個充電時間內對該電池充電(步驟d)。
在步驟a中,所述四個目標剩餘電量及五個充電電流係依一粒子群演算法決定。該粒子群演算法包括:以所述四個目標剩餘電量及五個充電電流代表一粒子空間;在一記憶體中儲存該粒子空間中之n個不同粒子的目前位置、目前速度、目前適應值及個別最佳位置,以及所述n個不同粒子之一群體最佳位置及一群體最佳適應值,n為大於1之整數;依所述目前位置、目前速度、個別最佳位置及群體最佳位置對各所述粒子進行一位置更新程序以產生一新位置以更新所述目前位置,依該新位置之內含資訊進行一適應值計算程序以獲得一適應值,及在該適應值小於該目前適應值時以該適應值更新該目前適應值;以及在n個所述目前適應值中之最小者小於該群體最佳適應值時,以該最小者更新該群體最佳適應值,並以與該最小者對應之所述目前位置更新該群體最佳位置。
另外,該位置更新程序包括以下公式:
x ij (n+1)=x ij (n)+v ij (n+1),其中,w為權重值(介於0~1之間),C 1 和C 2 分別為個體和群體學習因子,均介於1~4之間,x ij 代表粒子的位置,i為第幾個粒子,j為其粒子維度,v ij 表第i粒子j維度的速度,pBest ij 表個體粒子最佳的位置,gBest ij 表群體最佳的位置,r 1和r 2為介於[0,1]之間均勻分布亂數,n為疊代次數。
其中,s為階段數,I charge,s 為第s階段之充電電流,SOC charge,s 為第s階段之充入電量,SOC 0 為起始充電之剩餘容量,Q rate 為電池之額定容量,R s 為串聯等效電阻,R p 為並聯等效電阻,Objective為該適應值,α為權重值,(T now ,L now )為一目前總充電時間及一目前充電總損失,(T min ,L min )為一理想總充電時間及一理想充電總損失,(T max ,L max )為一上限總充電時間及一上限充電總損失。
另外,該權重值α例如但不限於為0.5。
其中,R eq 為等效阻抗,SOC為電池剩餘容量。
以下將針對本發明之原理進行說明:
請參照圖2,其繪示習知技術之多階段定電流(Multi-Stage Constant Current,MSCC)充電法之充電示意圖。
如圖所示,多階段定電流充電法進行充電時,當電池電壓達到上限額定電壓4.2V時電流會切換到下一階段,然後使用小於上一階段的充電電流繼續充電。若多階段定電流充電法之充電時間較短且電池溫升低,則可增加電池的循環壽命。但由於電池內部電化學特性相當複雜,因此為得到最佳化充電電流,習知技術利用模糊控制法,以充電時的電池溫升作為模糊控制的輸入,經過模糊控制計算後,輸出充電電流;或使用田口方法(Taguchi Methods,TM)來尋找最佳的充電電流;或使用蟻群系統(Ant colony system,ACS)求取最佳化充
電電流。然而由於上述方法均需複雜運算,必須藉由電腦來實現。
鋰離子電池等效模型:
由於電池內部電化學反應的影響會使得其輸出電壓隨使用時間增加而下降,因此建立正確電池等效模型有利於實驗之分析。最簡單之鋰電池等效電路模型可視為一恆電壓源,即輸出電壓維持不變,而電池在實際充放電時,由於電池內部化學反應之影響,使其輸出電壓會隨使用時間上升或下降。所以一個精確的電池等效電路模型必須考量自放電現象、電池內阻及暫態響應等因素。
本測試選用三星電子(SAMSUNG)推出的ICR18650-26H鋰離子電池,其參數規格如表1所示。
請一併參照圖3a至3b,其中圖3a其繪示戴維寧電池等效模之示意圖,圖3b其繪示圖3a的改良電池等效電路模型。
如圖3a所示,該模型與線性電池模型不同處在於增加一組R p -C p 並聯阻抗,因為能藉由C p 與R p 模擬電池的動態響應,而在模擬電池之充放電時較為精確,本發明將以此模型來推導多階段定電流充電法之最佳化充電模式。
如圖3b所示,本電池等效電路模型之組成包括電池等效電容(C eq )、串聯等效電阻(R s )及並聯等效電阻(R p )與電容(C p ),其中R p 與C p 係用來描述電池之暫態響應,跨於兩端之電壓為V transient ,而C eq 上之電壓V oc 為開路電壓(Open-circuit voltage,OCV),V T 為電池端電壓(Terminal voltage),其中電池等效電路模型之相關參數值與剩餘容量和溫度有關。
根據庫倫積分法,電池的剩餘容量SOC如方程式(1)所示。
其中SOC0為起始充電之剩餘容量,Q rate 為電池之額定容量,I charge 為當時的充電電流。由方程式(1)可得到充電時間△t如方程式(2)所示。
根據克希荷夫電壓定律,電池端電壓V T 如方程式(3)所示。
V T =V oc +I charge ×R s +V transient (3)
其中V transient =V Rp =V Cp 。
由於多階段定電流充電法之充電時間係由每一階段之充電電流透過庫倫積分所推得,因此完整之充電時間T s 如方程式(4)或方程式(5)所示。
其中n代表多階段定電流充電法之階段數,I charge,s 和SOC charge,s 分別為第s階段之充電電流和所充入的電量。
又根據庫倫積分法,每階段充電後之安培小時(△Ah s )容量和總安培小時容量Ah T 分別如方程式(6)及(7)所示。
△Ah s =I charge,s ×△t s (6)
而在多階段定電流充電法中,關於充電電流階數(n)的選擇,雖能藉由增加充電階數可增加充電效率及改善充電時間,但根據先前文獻研究指出階段數大於5之後的改善效益不高,且會增加控制電路複雜度及成本,故本發明充電階段數選定為五段,但不以此為限。
而各階段之充電損失△L s 和總充電損失L s 分別如方程式(8)及(9)
所示。
鋰電池之交流阻抗分析:
在二次電池之充放電技術中,交流阻抗(Alternating Current Internal Resistance,ACIR)分析係一種電化學分析之方式,主要為分析電池在不同狀態下之化學反應,藉此得到電池在不同狀態下之等效阻抗,而能得到鋰離子電池等效電路模型之參數R s 、R p 和C p 。而交流阻抗分析係利用小振幅交流弦波電壓或電流對電池端電極進行擾動,從中獲取交流阻抗數據,並以此數據建立等效電路。此外,藉由改變交流弦波的頻率,可以得到實、虛阻抗對頻率之響應變化曲線,稱為電化學阻抗頻譜(Electrochemical Impedance Spectrum,EIS)或奈氏圖(Nyquist-plot)。
請參照圖4,其繪示電池之交流阻抗特性示意之奈氏圖。
如圖所示,電池之交流阻抗特性之奈氏圖可劃分為低頻、中頻及高頻三個部份加以探討,分別為低頻的質傳作用(Mass Transport Effects),中頻的荷傳及電雙層作用(Double-Layer Effects)以及高頻時電磁效應(Electric and Magnetic Effects)。
量測交流阻抗時,需要在電池端電極加上交流電壓或電流進行擾動,以下說明如何進行交流阻抗分析之恆電位偵測(Potentio Electrochemistry Impedance Spectroscopy,PEIS),
PEIS為由交流阻抗分析儀產生一組變頻弦波電壓主動訊號對電池進行擾動,直到所有設定擾動頻率範圍處理完畢。恆電位模式下的擾動電壓需要注意擾動量不能過大,以避免干擾電池的平衡狀態,造成量測失真。
請一併參照圖5a至5b,其中圖5a其繪示交流阻抗之分析流程圖,圖5b其繪示交流阻抗測試之實驗流程圖。
在進行恆電位擾動偵測時,電池因擾動而產生相對應之電流變化,因此可偵測出電池電流的振幅及相位角度,接著將該電流參數進行訊號調
節及轉換,再利用擾動電壓與相對產生的電流即可計算出阻抗大小及相位角,如方程式(10)所示。
當確認電池規格符合實驗需求後及完成所有頻率響應量測後,則可執行交流阻抗參數分析,分析流程如圖5a所示。
本發明之交流阻抗分析儀選用Bio-Logic公司型號VSP-300之多功能模組化恆電位儀,並搭配EC-Lab軟體介面進行多種電池分析及實驗,但不以此為限。
由於電池之剩餘容量和電池交流阻抗大小有相對應之關係,實驗進行時,電池的剩餘容量精確度與實驗所需花費之時間需做權衡,因為每次測量前均需要靜置一小時等電化學反應穩定後才能進行,若以1%之剩餘容量作為量測間隔,則整個過程至少需100個小時之休息時間。為了權衡電池等效電路模型參數之精確度與所需花費之量測時間,將以每1%的剩餘容量作為每筆資料的間隔。
如圖5b所示,實驗開始前需將待測電池完全充飽,充飽後以定電流放電至所需之剩餘容量,並使電池休息1小時再進行下一次的交流阻抗分析。
電池之剩餘容量又與電池交流阻抗大小有著相對應的關係,而在實驗過程中,電池剩餘容量精確度與實驗所需花費時間需要做權衡,因為每一次的測量電池都需要靜置約一小時後才能進行,所以以每1%的剩餘容量作為間隔的所有資料,至少需花費100個小時的休息時間進行靜置。
請一併參照圖6a至6d,其中圖6a其繪示EC-Lab交流阻抗分析之操作介面示意圖,圖6b其繪示交流阻抗R s 、R p 與等效阻抗R eq 對剩餘容量百分比之關係曲線圖(環境溫度5℃),圖6c其繪示交流阻抗R s 、R p 與等效阻抗R eq 對剩餘容量百分比之關係曲線圖(環境溫度25℃),圖6d其繪示交流阻抗R s 、R p 與等效阻抗R eq 對剩餘容量百分比之關係曲線圖(環境溫度45℃)。
電池之交流阻抗測量完成,可得不同剩餘容量之奈奎斯特阻抗
圖,如圖6a所示,EC-Lab中的功能Z Fit提供了多組的等效電路模型,可提供使用者分析奈奎斯特阻抗圖。
本發明採用戴維寧電池等效電路模型,其對應的參數鑑別操作步驟如下:
Step 1:選擇量測的數據呈現方式,奈奎斯特圖需選擇Nyquist Impedance,並選擇Z cycle顯示不同剩餘容量下的奈奎斯特圖。
Step 2:Z Fit功能可於使用者操作介面之Analysis清單中的Electrochemical Impedance Spectroscopy選擇,而本發明採用戴維寧電池等效電路模型,故選用兩個電阻R 1、R 2及一個電容C 2組成之等效電路模型進行參數擬合。
Step 3:交流阻抗為每1%之剩餘容量測量一次,故先選擇欲分析剩餘容量下之奈氏圖後,再圈選曲線擬合範圍。
Step 4:執行曲線擬合,先將曲線擬合最小化(Minimize)後,再計算(Calculate),得到曲線的參數。
Step 5:擬合後的參數將顯示於此欄,可將其紀錄下來。
如圖6b至6d,將操作於環境溫度分別設在5℃、25℃及45℃,以1C(2.6 A)定電流充電至端電壓4.2V,再以4.2V電壓浮充至充電電流小於0.02C(52mA)後視為電池充飽。電池充飽後,以每1%之剩餘容量放電後,再休息1小時取其開路電壓(V oc )及在不同溫度下由參數鑑別操作步驟量測等效電路內之交流阻抗(串聯阻抗R s 、並聯阻抗R p 、並聯電容C p )參數值,其中等效阻抗(R eq )=串聯阻抗(R s )+並聯阻抗(R p )。
本發明採用粒子群演算法(Particle Swarm Optimization,PSO)尋找最佳化充電電流設定值:
粒子群演算法之基本精神來自於群體動物間相互溝通及交換資訊等行為,並將其用於搜尋最佳解之相關問題。PSO藉由互相傳遞訊息來達到個體和群體找尋食物、聚集等目的,群體中的個體稱為粒子(Particle),在解空間中所有粒子都有本身所對應之適應值,並且每個粒子都知道自己目前為止之最佳適應值及其最佳位置,稱之為個體粒子的最佳值(Particle best value,pBest),
這項資訊為每個粒子自己所擁有的經驗,同時每個粒子亦知道全部粒子之最佳解以及其位置,稱之為群體粒子的最佳值(Globe best value,gBest)。經過每次疊代更新,粒子會以個體經驗及群體經驗作為參考,並且更新個體之速度與位置。
請參照圖7,其繪示粒子群演算法之粒子在解空間移動之示意圖。
如圖所示,一開始將所有粒子隨機散佈在解空間中,若有粒子能接近區域內最佳解附近,該區域的粒子將往此區域之最佳解進行搜尋,但此區域有可能只是區域最佳解,因此必須由各粒子群的搜索結果修正群體最佳的適應值及位置,進而接近全域最佳解。
粒子之移動速度如方程式(11)所示。
其中,其中,w為權重值(介於0~1之間),C 1 和C 2 分別為個體和群體學習因子,均介於1~4之間,x ij 代表粒子的位置,i為第幾個粒子,j為其粒子維度,v ij 表第i粒子j維度的速度,pBest ij 表個體粒子最佳的位置,gBest ij 表群體最佳的位置,r 1介於[0,1]之間均勻分布亂數,r 2介於[0,1]之間均勻分布亂數,n為疊代次數。
PSO一開始使用亂數來初始化整個粒子群使其散佈於所要搜尋之解空間,而粒子點會往個體最佳適應值及群體最佳適應值之方向移動,經過每次的疊代更新,粒子會以個體粒子的經驗以及群體的經驗作為參考,並且更新粒子個體的速度與位置,粒子經過疊代後更新之位置如方程式(12)所示。
x ij (t+1)=x ij (t)+v ij (t+1) (12)
要將PSO應用於搜尋最佳化SOC轉態值與充電電流值,須先定義粒子與解空間及要解決之最佳化問題。粒子位置可以二維表示為各階段SOC轉態設定值及充電電流之C數,搜尋空間為每1%之SOC與介於0~1C間之充電電流組合,而本發明以縮短充電時間及降低充電損失作為最佳化目標,故需
找出影響最佳化目標的因素。
請一併參照圖8a至8d,其中圖8a其繪示本發明之電池剩餘容量與等效阻抗之擬合曲線圖,圖8b其繪示量測阻抗與擬合曲線之誤差量,圖8c其繪示0.5C CC-CV充電法之模擬結果,圖8d其繪示2C CC-CV充電法之模擬結果。
影響充電時間及充電損失之因素包含各階段充電電流(I1~I5)及不同SOC下之交流阻抗大小,本發明將依此尋找基於SOC之多階段充電法之最佳化設定。
因為充電損失與電池內阻有關,因此將所量測之交流阻抗數據匯入MATLAB後,並使用曲線擬合(cftool)功能建立剩餘容量與阻抗之擬合曲線如方程式(13)所示。
其中SOC為電池剩餘容量,R eq 為等效阻抗,而各擬合係數如表2所示。
如圖8a至8b所示,實際量測阻抗與擬合曲線之誤差量,誤差範圍均在0.5%以內,因此可證應用曲線擬合能得知不同剩餘容量下的交流阻抗大小。
PSO搜尋解空間中所有粒子都是藉由價值函數來決定各粒子之
適應值(fitness value),因此適應值之評估方式極為重要,將影響到最佳化之結果,為同時考量充電時間及充電損失兩種不同量之參數,本發明將此兩參數正規化並加上權重(weight)以強調此二參數之重要性。
本發明之評估方法為利用兩點間直線之距離做為適應值的分數,當A點(T s_now ,L s_now )與理想最佳解(T s_min ,L s_min )之距離d越小,其適應值評估越佳,兩點間距離d如方程式(14)所示。
其中T s_now 和L s_now 代表目前充電時間和損失,T s_min 和L s_min 代表最小的充電時間和損失。
為了調整此二參數之重要性,加入權重係數α,改寫如方程式(15)所示。
其中,首先以MATLAB模擬CC-CV充電法,而模擬條件為選用電池廠商建議的標準充電電流範圍作為最大充電電流,充電截止電流為0.02 C所得之結果如表3所示。
可以得知,方程式(15)的T s_max 、T s_min 、L s_max 及L s_min 分別為0.5 C與2 C條件下所模擬出之最大及最小的充電時間與充電損失,作為最佳化演算法模擬的重要參數,即正規化採用之參考資料點(T s_min 、T s_max 、L s_min 、L s_max )為(3241、8790、792.16、2872.43),0.5 C CC-CV與2 C CC-CV充電法模擬結果分別如圖8c與圖8d所示。
本發明之基於SOC轉態值之最佳化多階段定電流充電法,一開始先初始化所有粒子之速度和位置(SOC1-4,I1-5),接著計算充電時間、充電損失和評估適應值,並更新pBest和gBest,接著更新每個粒子之速度和位置,若到達終止條件或疊代次數,則停止解空間搜尋,並輸出最佳解。
請一併參照圖9a至9b,其中圖9a其繪示本發明實作採用之充電機系統架構示意圖,圖9b其繪示五階段充電法之操作流程圖。
如圖9a所示,本發明實作採用之充電機系統架構採用全橋相移轉換器,並使用Texas Instruments公司推出之TMS320F280049數位訊號處理器作為電池充電機之數位控制核心,利用取樣電路取樣輸出電壓及電流,將訊號透過類比/數位轉換後再送入數位濾波器進行濾波,接著將濾波結果透過PID補償器計算功率開關之相移量,最後由PWM模組將PWM訊號輸出至功率開關控制全橋相移轉換器,進而對電池進行充電。另外也使用LabVIEW人機介面透過USB-6009資料擷取卡讀取電池充電過程的電壓電流,方便使用者即時監控電池充電過程,充電機電路設計規格如表4所示。
本發明係採用TMS320F280049微控制器來實現數位有限脈衝響應(FIR)濾波器及增量型比例-積分-微分(Proportional-Integral-Derivative,PID)控制器,程式可分為主程式、ADC中斷副程式及TIMER1中斷副程式三個部分,一開始時首先會針對程式所需的變數進行宣告,並將系統CLK初始化,接著依序對GPIO模組、類比/數位轉換器ADC模組、TIMER1模組,以及脈波寬度調變模組進行模組設定與致能,之後進入無窮迴圈等待中斷向量旗標發生;ADC中斷副程式部分又可以細分為取樣濾波、五階段定電流充電及增量型PID三部分,TIMER1中斷副程式設定1秒產生中斷一次,用來計數每階段的充電時間,並等待計數秒數到達,重置下一階段TIMFR1計數。
如圖9b所示,本發明藉由PSO最佳化演算法求得五個階段的定電流值與轉態之剩餘容量設定值,接著將結果寫入數位信號處理器,並以全橋相移轉換器作為充電機對電池充電來實現,當電池接上充電機時,開始以定電流進行充電,其設定第一階段充電電流I 1,且TIMER1進入中斷,開始計數至第一階段的充電時間,接著重置TIMER1計數,進入下一階段電流I 2,並TIMER1再次進入中斷,開始計數至第二階段的充電時間,接著重置TIMER1計數,以此類推,進而完成五階段電池充電。
以下將依序從演算法之模擬結果、實測環境及實測結果進行說明與分析,驗證本技術所提最佳化充電法之優點,並提供充電過程之實際量測波形。
模擬結果:
本發明模擬時選用電池之規格與適應值評估之參數如表5所示。
表中定義SOC於100%與0%時之開路電壓為環境溫度在25℃下所量測而得,而適應值評估權重係數α設為0.5,代表充電時間與充電損失具有相同的重要性。
影響充電時間與損失之因素包含各階段充電電流(I 1~I 5)及每1% SOC下所量測出之交流阻抗大小,另外,實測分析得知交流阻抗大小也會受測試環境溫度影響。因此,基於充電安全考量,多階段定電流充電法之各階段電流介於0C至1C間,本發明之粒子群範圍設定如表6所示,其中粒子的維度為9,包含四個轉態之剩餘容量(SOC 1~SOC 4)與五個階段之電流值(I 1~I 5)。
結果比較:
由於影響充電時間及充電損失之因素包含了各階段充電電流(I 1~I 5)及每1%之剩餘容量下的交流阻抗大小,Case1(本發明充電法)將各階段之SOC轉態點亦列入PSO加以考量。此外,為搜尋最佳化充電結果,本發明設定各階段的充電電流值無須小於上一階段電流,使各階段之充電電流值依演算法收斂至最佳位置,本發明另附加限制條件使電池內電壓不得超過4.25V,否則不更新粒子之最佳位置,因此本發明亦考量到電池實際充電之安全性,確保電池不會有過充的問題,其模擬輸出結果(包含充電時間、充電損失、充電後之剩餘容量和評估適應值)如表7a所示。
Case2(習知技術一)提出以固定之剩餘容量作為轉態點,並以田口法來搜尋各階段之電流值,在制定田口直交表時,其限制下一階段之充電電流值須小於上一階段,以此為基礎,並以固定20%間隔之剩餘容量作為轉態點,應用PSO求取最佳化電流,其模擬輸出結果(包含充電時間、充電損失、充電後之剩餘容量和評估適應值)如表7b所示。
Case3(習知技術二)係只考量充電時間最佳化之案例,其各階段電流之關係與Case2相同,其模擬輸出結果(包含充電時間、充電損失、充電後之剩餘容量和評估適應值)如表7c所示。
Case4(習知技術三)將電流保持不變,但SOC轉態點改為20%、40%、60%以及80%,其模擬輸出結果(包含充電時間、充電損失、充電後之剩餘容量和評估適應值)如表7d所示。
Case5與Case6則用來驗證本發明(Case1)是否為最佳解,其分別將剩餘容量轉態點與各階段之電流值進行微調,其模擬輸出結果(包含充電時間、充電損失、充電後之剩餘容量和評估適應值)分別如表7e及7f所示。
六種充電法之模擬輸出結果比較如表7g所示,可以得知本發明之適應值(0.2287)確實最佳。
實驗環境設置:
請一併參照圖10a至10e,其中圖10a其繪示本發明進行充電實測之實驗平台示意圖,圖10b其繪示圖10a之五階段充電法之操作流程圖,圖10c其繪示習知技術之CC-CV實測0.5C(1.3 A)之充電電壓、電流波形圖,圖10d其繪示習知技術之CC-CV實測1C(2.6 A)之充電電壓、電流波形圖,圖10e其繪示習知技術之CC-CV實測2C(5.2 A)之充電電壓、電流波形圖。
以下對電池進行充電實測以驗證本發明所得設定值為最佳,實測平台之電池充放電機係採用Bio-Logic公司之可程控充放電機BCS-815及監控介面,其能設定多種充電法則,且能即時監控電池的電壓、電流及溫度等,並可將充放電資料回傳至電腦紀錄,而電池測試環境皆放置於RI-45恆溫箱內並控制在環境溫度為25℃。
如圖10b所示,實測一開始使用I 1充電,當預設之各階段充電時間到達即使用下個階段電流充電,直到五個階段完成。完成一次充電後,電池休息約1小時以緩和電池內部之電化學反應,再以定電流0.2C(0.52A)進行放電,直到端電壓小於等於3V完成實驗流程。
因為充電損失不易量測,一般而言,電池的充電損失與充電溫升有直接關聯,且充電損失越高時,電池溫升也會升高,因此將以電池之表面溫升作為實測評比的標準。在實測結果中,比較項目包含:充電時間、充電溫升、充放電容量及充電效率,並將本發明與習知技術之CC-CV進行比較。如圖10c至10e所示,其充電截止條件皆為電流低於0.02C(52mA)視為充飽,而0.5C(1.3 A)、1C(2.6 A)及2C(5.2 A)之完整充電時間分別為9889秒、6621秒及5208秒。
請一併參照圖11a至11f,其中圖11a其繪示Case1(本發明充電法)之實測波形圖,圖11b其繪示Case2(習知技術一)之實測波形圖,圖11c其繪示Case3(習知技術二)之實測波形圖,圖11d其繪示Case4(習知技術三)之實測波形圖,圖11e其繪示Case5(微調剩餘容量轉態點)之實測波形圖,圖11f其繪示Case6(微調各階段之電流值)之實測波形圖,圖11g其繪示0.5C(1.3 A)、1C(2.6 A)、2C(5.2 A)、Case1至Case6之充電溫升曲線圖,圖11h其繪示Case1至Case6之充電溫升曲線圖。
各方法之實測充電時間比較如表8a所示。
上述實驗結果之數據比較如表8b所示,充電效率定義為(電池放出容量)/(電池充入容量)。可以得知本發明能將電池容量充至CC-CV所得容量的97%以上。此外,由於電池表面之最高溫升無法表示完整充電過程溫度的變化,故本實測採用平均溫升進行。
上述實驗結果之平均溫升與適應值比較如表8c所示。
可以得知,本發明(Case1)之適應值(0.3493)為最佳的,與Case2(習知技術一)相比充電時間快3.83%,且平均溫升低0.29℃;與Case3(習知技術二)相比充電時間快14.12%,且平均溫升低0.69℃,本發明能在不過充電的情況下,優化充電時間與充電損失;而由Case4、Case5、Case6實測結果可驗證,各階段SOC轉態設定值與充電電流皆會影響最佳化目標,綜上所述,本發明為最佳解。又根據實測數據,表8c之適應值函數如方程式(16)所示。
藉由前述所揭露的設計,本發明乃具有以下的優點:
1.本發明揭露一種可縮短充電時間之多階段定電流充電方法,其係透過交流阻抗分析以建立鋰離子電池之等效模型,藉由庫倫積分與等效電路模型求得充電時間與充電損失,並利用粒子群演算法尋找基於剩餘容量之多階段定電流最佳化充電設定,藉此實現縮短充電時間及降低充電損失之最佳化多階段定電流充電法。
2.本發明揭露一種可縮短充電時間之多階段定電流充電法,係利用粒子群演算法尋找最佳化充電電流設定值,僅需使用電池等效模型,不需制
訂田口直交表之每階段初始電流值之設置及多次耗時費工之實際實驗即能得到最佳化電流設定。
3.本發明揭露一種可縮短充電時間之多階段定電流充電法,相較習知技術之CC-CV充電法可將電池容量充至CC-CV所得容量的97%以上,其充電時間能縮短14.12%且平均溫升低0.69℃。
本發明所揭示者,乃較佳實施例,舉凡局部之變更或修飾而源於本發明之技術思想而為熟習該項技藝之人所易於推知者,俱不脫本發明之專利權範疇。
綜上所陳,本發明無論就目的、手段與功效,在在顯示其迥異於習知之技術特徵,且其首先發明合於實用,亦在在符合發明之專利要件,懇請 貴審查委員明察,並祈早日賜予專利,俾嘉惠社會,實感德便。
步驟a:在一電池之初始剩餘電量和100%電量之間設置四個目標剩餘電量以定義出六個剩餘電量值,並依所述六個剩餘電量值定義出五個階段,其中各所述階段之期初剩餘電量和期末剩餘電量係由兩個相鄰的所述剩餘電量值決定
步驟b:在所述五個階段中對應設定五個充電電流
步驟c:依所述六個剩餘電量值及所述五個充電電流執行五次充電時間計算以產生五個充電時間,其中各所述充電時間計算表為:(該電池之額定容量)*(兩個相鄰的所述剩餘電量值之絕對差值)/(一對應之所述充電電流)
步驟d:依所述五個充電電流在對應之所述五個充電時間內對該電池充電
Claims (5)
- 一種可縮短充電時間之多階段定電流充電方法,其包括以下步驟:在一電池之初始剩餘電量和100%電量之間設置四個目標剩餘電量以定義出六個剩餘電量值,並依所述六個剩餘電量值定義出五個階段,其中各所述階段之期初剩餘電量和期末剩餘電量係由兩個相鄰的所述剩餘電量值決定;在所述五個階段中對應設定五個充電電流;依所述六個剩餘電量值及所述五個充電電流執行五次充電時間計算以產生五個充電時間,其中各所述充電時間計算表為:(該電池之額定容量)*(兩個相鄰的所述剩餘電量值之絕對差值)/(一對應之所述充電電流);以及依所述五個充電電流在對應之所述五個充電時間內對該電池充電;其中所述四個目標剩餘電量及五個充電電流係依一粒子群演算法決定,且該粒子群演算法包括:以所述四個目標剩餘電量及五個充電電流代表一粒子空間;在一記憶體中儲存該粒子空間中之n個不同粒子的目前位置、目前速度、目前適應值及個別最佳位置,以及所述n個不同粒子之一群體最佳位置及一群體最佳適應值,n為大於1之整數;依所述目前位置、目前速度、個別最佳位置及群體最佳位置對各所述粒子進行一位置更新程序以產生一新位置以更新所述目前位置,依該新位置之內含資訊進行一適應值計算程序以獲得一適應值,及在該適應值小於該目前適應值時以該適應值更新該目前適應值;以及在n個所述目前適應值中有小於該群體最佳適應值時,以該最小者更新該群體最佳適應值,並以與該最小者對應之所述目前位置更新該群體最佳位置。
- 如申請專利範圍第3項所述之可縮短充電時間之多階段定電流充電方法,其中該權重值α為0.5。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW110120920A TWI760246B (zh) | 2021-06-09 | 2021-06-09 | 一種可縮短充電時間之多階段定電流充電方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW110120920A TWI760246B (zh) | 2021-06-09 | 2021-06-09 | 一種可縮短充電時間之多階段定電流充電方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
TWI760246B true TWI760246B (zh) | 2022-04-01 |
TW202248873A TW202248873A (zh) | 2022-12-16 |
Family
ID=82198750
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW110120920A TWI760246B (zh) | 2021-06-09 | 2021-06-09 | 一種可縮短充電時間之多階段定電流充電方法 |
Country Status (1)
Country | Link |
---|---|
TW (1) | TWI760246B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI465939B (zh) * | 2013-01-07 | 2014-12-21 | Univ Lunghwa Sci & Technology | A Multi - stage Fast Charge Method for Optimizing Lithium Batteries |
TWI527337B (zh) * | 2013-12-18 | 2016-03-21 | Univ Lunghwa Sci & Technology | The method of determining the charging current of multi - stage lithium battery using Taguchi orthogonal table and fuzzy algorithm |
TWI686617B (zh) * | 2019-03-20 | 2020-03-01 | 龍華科技大學 | 一種最佳化輸出電流值之多階段定電流充電法 |
-
2021
- 2021-06-09 TW TW110120920A patent/TWI760246B/zh active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI465939B (zh) * | 2013-01-07 | 2014-12-21 | Univ Lunghwa Sci & Technology | A Multi - stage Fast Charge Method for Optimizing Lithium Batteries |
TWI527337B (zh) * | 2013-12-18 | 2016-03-21 | Univ Lunghwa Sci & Technology | The method of determining the charging current of multi - stage lithium battery using Taguchi orthogonal table and fuzzy algorithm |
TWI686617B (zh) * | 2019-03-20 | 2020-03-01 | 龍華科技大學 | 一種最佳化輸出電流值之多階段定電流充電法 |
Non-Patent Citations (2)
Title |
---|
期刊韕Thanh Tu Vo, et al., 韕 韕 " New charging strategy for lithium-ion batteries based on the integration of Taguchi method and state of charge estimation", 韕Vol. 273,韕Journal of Power Sources, 韕28 Sep., 2014,韕pp. 413-422.韕; * |
期刊韕Y. Liu, et al., 韕 韕"Search for the Optimal Charging Pattern of Multi-stage Constant Current Charging Method Using Particle Swarm Optimization", 韕2018 7th 韕International Congress on Advanced Applied Informatics (IIAI-AAI), 韕2018,韕pp. 536-541.韕 * |
Also Published As
Publication number | Publication date |
---|---|
TW202248873A (zh) | 2022-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106451600B (zh) | 为电池快速充电的设备和方法 | |
Shen et al. | Accurate state of charge estimation with model mismatch for li-ion batteries: A joint moving horizon estimation approach | |
Chen et al. | State-of-charge estimation of lithium-ion batteries based on improved H infinity filter algorithm and its novel equalization method | |
CN110632528B (zh) | 一种基于内阻检测的锂电池soh估算方法 | |
CN109991554B (zh) | 一种电池电量检测方法、装置及终端设备 | |
CN107690585A (zh) | 用于确定锂硫电池组的健康状况和充电状态的方法和装置 | |
KR101402802B1 (ko) | 배터리 셀의 전압 변화 거동을 이용한 셀 밸런싱 장치 및 방법 | |
US9645200B2 (en) | Battery power measuring method, measuring device and battery-powered equipment | |
CN111426956B (zh) | 考虑温度及迟滞效应的分数阶动力电池soc估算方法 | |
RU2006119641A (ru) | Устройство и способ для оценки уровня заряженности батареи с использованием нейронной сети | |
CN104956233A (zh) | 电池状态推断装置 | |
CN111812515A (zh) | 基于XGboost模型的锂离子电池荷电状态估算 | |
CN113359044A (zh) | 测量电池剩余容量的方法、装置及设备 | |
CN108829911A (zh) | 一种开路电压与soc函数关系优化方法 | |
CN114091282B (zh) | 基于分数阶模型的锂离子电池状态估计方法及系统 | |
CN110244236A (zh) | 一种锂离子电池包soc估计方法、装置及终端设备 | |
CN114928133A (zh) | 一种用于锂离子电池的充电控制方法、系统、设备和介质 | |
CN114755582A (zh) | 一种环境温度不同时锂离子电池健康状态预估方法 | |
TWI760246B (zh) | 一種可縮短充電時間之多階段定電流充電方法 | |
Huang et al. | State of charge estimation of lithium-ion battery based on second-order extended Kalman filter | |
TWI686617B (zh) | 一種最佳化輸出電流值之多階段定電流充電法 | |
Shen et al. | A state of charge estimation method based on APSO-PF for lithium-ion battery | |
Wang et al. | An Online Method for Power State Estimation of Lithium-Ion Batteries under the Constraints of the Fusion Model Considering Temperature Effect | |
CN113625183A (zh) | 一种电池包寿命预测方法及电池包模拟系统 | |
Wu et al. | Multistage fast charging optimization protocol for lithium-ion batteries based on the biogeography-based algorithm |