周智勇 劉愛華 楊占錄
(海軍潛艇學(xué)院機(jī)電研究所,青島266042)
潛艇蓄電池的電量有限,儲(chǔ)能降低后潛艇必須在適當(dāng)?shù)臅r(shí)機(jī)浮起到水面或通氣管狀態(tài)進(jìn)行充電。合理管理蓄電池能量,降低潛艇的暴露率,首先就需要建立了較為準(zhǔn)確的蓄電池?cái)?shù)學(xué)模型。
蓄電池的充放電伴隨著復(fù)雜的化學(xué)變化和物理變化過程,使得影響蓄電池剩余容量和荷電狀態(tài)(SOC)的因素眾多。蓄電池端電壓、電解液密度和溫度、攪拌系統(tǒng)、使用周期、蓄電池內(nèi)阻等都對(duì)剩余容量產(chǎn)生影響,并且都是非線性關(guān)系[1]。因此,建立其數(shù)學(xué)模型的難度較大。本文以實(shí)驗(yàn)數(shù)據(jù)為基礎(chǔ),結(jié)合神經(jīng)網(wǎng)絡(luò)技術(shù),建立了比較準(zhǔn)確的蓄電池模型。
為準(zhǔn)確、全面地研究潛艇蓄電池的電特性,本文設(shè)計(jì)并建立了功能完善的蓄電池充放電實(shí)驗(yàn)室。實(shí)驗(yàn)室主要包括蓄電池、蓄電池附屬系統(tǒng)(包括低壓氣攪拌系統(tǒng)、蒸餾水冷卻系統(tǒng))、充放電裝置、蓄電池與充放電裝置間的可變換式連接系統(tǒng)、通風(fēng)系統(tǒng)、氫氣監(jiān)測(cè)裝置等。
充放電模塊是充放電裝置中的執(zhí)行機(jī)構(gòu)。三相交流電經(jīng)變壓器降壓后進(jìn)行整流。整流后的直流電壓(5~9 V)加在晶體管群上。主控模塊根據(jù)設(shè)定的充放電電流的大小,通過控制晶體管群的基極電壓,控制晶體管群上的集電極電流。這個(gè)電流經(jīng)充放電轉(zhuǎn)換模塊控制后,加在被測(cè)蓄電池上。同時(shí),數(shù)據(jù)采集模塊采集被測(cè)蓄電池上的電壓、電流、溫度參數(shù),經(jīng)過主控模塊將數(shù)據(jù)傳到上位機(jī)。充放電模塊的工作原理如圖1所示。
圖1 充放電模塊基本原理圖
為了與實(shí)艇蓄電池的充放電條件相近并保證蓄電池組可靠工作,實(shí)驗(yàn)室設(shè)置了蓄電池電解液攪拌系統(tǒng)、蒸餾水冷卻系統(tǒng)、通風(fēng)系統(tǒng)和氫氣監(jiān)測(cè)裝置。
蓄電池溫度過高時(shí)(溫度值可以設(shè)定)用蒸餾水冷卻系統(tǒng)進(jìn)行冷卻,共四個(gè)冷卻回路,每組蓄電池一個(gè),回路內(nèi)冷卻水的流向可以控制。
電解液攪拌系統(tǒng)用低壓氣攪拌電解液,使其上下兩部分的密度均勻。每塊蓄電池一個(gè)攪拌支路,共12個(gè)支路。每個(gè)支路供給每只蓄電池壓縮空氣,為其內(nèi)循環(huán)裝置提供動(dòng)力。
密度的測(cè)量是本試驗(yàn)的關(guān)鍵環(huán)節(jié),為了保證測(cè)量精度,實(shí)驗(yàn)采用了奧地利 Anton Paar公司DMA35型便攜式密度計(jì),密度測(cè)量范圍是0.8~2 g/cm3,精度為 0.001 g/cm3,分辨率為 0.0001 g/cm3,該密度計(jì)同時(shí)可測(cè)量已抽取液體的溫度,測(cè)量范圍是為0~40 ℃,精度是0.2 ℃,分辨率為0.1 ℃。
實(shí)驗(yàn)室的通風(fēng)系統(tǒng)、氫氣監(jiān)測(cè)裝置、蓄電池與充放電裝置間的連接系統(tǒng)等不再介紹。
BP神經(jīng)網(wǎng)絡(luò)是一種單向傳播的多層前向網(wǎng)絡(luò),它除了有一個(gè)輸入層和一個(gè)輸出層外,還有一層或多層的隱層,同層節(jié)點(diǎn)中沒有任何耦合。輸入信號(hào)從輸入層節(jié)點(diǎn)依次傳過各隱層節(jié)點(diǎn),然后傳到輸出節(jié)點(diǎn),每一層節(jié)點(diǎn)的輸入只影響下一層節(jié)點(diǎn)的輸出。BP神經(jīng)網(wǎng)絡(luò)具有的非線性映射能力保證其能夠成功逼近各種非線性函數(shù),其并行結(jié)構(gòu)加速了運(yùn)算速度。因此在目前的人工神經(jīng)網(wǎng)絡(luò)的實(shí)際應(yīng)用中,大部分的神經(jīng)網(wǎng)絡(luò)模型是采用BP神經(jīng)網(wǎng)絡(luò)及它的變化形式[2]。
如果模型的輸入輸出均已確定,則 BP神經(jīng)網(wǎng)絡(luò)的輸入層、輸出層神經(jīng)元的個(gè)數(shù)是確定的,關(guān)鍵是確定隱含層的數(shù)量及每個(gè)隱含層神經(jīng)元的個(gè)數(shù)。經(jīng)過多次實(shí)驗(yàn),證明隱含層確定為3層能獲得較好的預(yù)測(cè)精度。每個(gè)隱含層的數(shù)量由人工確定則費(fèi)時(shí)費(fèi)力,可以編寫合適的軟件,由其自動(dòng)判斷[3]。實(shí)踐表明,總神經(jīng)元的數(shù)量過少時(shí)不能達(dá)到設(shè)定的訓(xùn)練目標(biāo),神經(jīng)元數(shù)量過多則會(huì)造成網(wǎng)絡(luò)訓(xùn)練時(shí)的精度很高但用于預(yù)測(cè)時(shí)的精度反而下降,因此選擇適當(dāng)數(shù)量的神經(jīng)元至關(guān)重要[4]。
編寫程序建立了神經(jīng)網(wǎng)絡(luò)。該程序設(shè)置每個(gè)隱含層神經(jīng)元的個(gè)數(shù)為3~9個(gè),并實(shí)現(xiàn)了隱含層神經(jīng)元數(shù)量的自動(dòng)選擇。如果訓(xùn)練完成的網(wǎng)絡(luò)在測(cè)試時(shí)得到的測(cè)試誤差絕對(duì)值的最大值<0.005則停止訓(xùn)練,并記錄下當(dāng)前工作空間作為最終的網(wǎng)絡(luò),該網(wǎng)絡(luò)完全滿足模型計(jì)算的精度要求。同時(shí)每次訓(xùn)練后的測(cè)試都計(jì)算仿真測(cè)試誤差的平均絕對(duì)誤差mae和均方根誤差mse。并分別記錄下全部訓(xùn)練過程中 mae最小的網(wǎng)絡(luò)Fine_BP_FD_DL_net_Mae和 mse最小的網(wǎng)絡(luò)Fine_BP_FD_DL_net_Mse。如果最終網(wǎng)絡(luò)的計(jì)算精度也不能達(dá)到要求,則改變隱含層神經(jīng)元的數(shù)量限制重新計(jì)算。以下為神經(jīng)網(wǎng)絡(luò)建立與測(cè)試的源代碼[5]。
通過該程序獲得的 mae、mse最小的神經(jīng)網(wǎng)絡(luò)有三個(gè)隱含層,每個(gè)隱含層神經(jīng)元數(shù)量分別為8、6、9,將該網(wǎng)絡(luò)命名為Fine_BP_FD_DL_net_Mae。
該網(wǎng)絡(luò)的總體結(jié)構(gòu)如圖2所示。Layer1層的結(jié)構(gòu)如圖 3所示。Layer2、Layer3、Layer4層的結(jié)構(gòu)與Layer1層的結(jié)構(gòu)相似。
圖2 神經(jīng)網(wǎng)絡(luò)的總體結(jié)構(gòu)
用1286組測(cè)試數(shù)據(jù)進(jìn)行測(cè)試,得到的誤差統(tǒng)計(jì)如表1所示。其中絕對(duì)誤差小于0.005g/cm的數(shù)據(jù)占總測(cè)試數(shù)據(jù)的 85.77%;絕對(duì)誤差小于0.01g/cm3的數(shù)據(jù)占總測(cè)試數(shù)據(jù)的 97.05%。與恒流放電相似,預(yù)測(cè)時(shí)間較長(zhǎng)時(shí)密度的計(jì)算誤差相對(duì)較大,預(yù)測(cè)時(shí)間較短時(shí)計(jì)算誤差較小。
圖3 Layer1層的總體結(jié)構(gòu)
表1 恒功率放電密度預(yù)測(cè)誤差分布統(tǒng)計(jì)
本文所建立的鉛酸蓄電池充放電模型可以完整地反映該型蓄電池的充放電過程,以此模型為基礎(chǔ)研制的蓄電池電量預(yù)測(cè)裝置經(jīng)實(shí)艇航行驗(yàn)證,精度達(dá)到了設(shè)計(jì)要求。
[1]朱松然.蓄電池手冊(cè)[M].天津: 天津大學(xué)出版社,1998: 1~47.
[2]薛建軍.人工神經(jīng)網(wǎng)路法預(yù)測(cè) MH-Ni蓄電池容量[J].電源技術(shù),2003, (03): 305~307.
[3]雷肖.電動(dòng)車電池SOC估計(jì)的徑向基函數(shù)神經(jīng)網(wǎng)絡(luò)方法[J].電工技術(shù)學(xué)報(bào),2008, (05): 81~87.
[4]李立偉.基于改進(jìn)灰色模型的蓄電池剩余容量預(yù)測(cè)[J].電源技術(shù),2006, (12): 1006~1008.
[5]葛哲學(xué),孫志強(qiáng).神經(jīng)網(wǎng)絡(luò)理論與 MATLABR2007實(shí)現(xiàn)[M].北京: 電子工業(yè)出版社,2007: 108~116.