孫躍 楊芳勛 戴欣
(重慶大學(xué)自動(dòng)化學(xué)院,重慶400030)
無線電能傳輸本質(zhì)上是一種借助于空間無形軟介質(zhì)(如磁場、電場、激光、微波等)將電能由源極傳遞至受電極的全新電能供給模式[1-4],其改變了傳統(tǒng)的只能依靠電導(dǎo)體直接傳輸電能的供電方式,是電能傳輸?shù)囊环N革命性進(jìn)步.感應(yīng)耦合電能傳輸(ICPT)作為無線電能傳輸?shù)囊环N,具有安全、便捷、綠色、靈活的特征,已引起了眾多國內(nèi)外學(xué)者的關(guān)注[5-8],并在很多領(lǐng)域得到了廣泛的應(yīng)用,如軌道交通[9]、生物醫(yī)療[10]、移動(dòng)電子設(shè)備[11-12]等.2007 年,美國麻省理工學(xué)院的研究團(tuán)隊(duì)通過磁共振原理實(shí)現(xiàn)了2 m范圍內(nèi)60 W的無線電能傳輸,效率約為40%[13].
縱觀目前的ICPT技術(shù),普遍是采用“廣播式”的供電模式,即一個(gè)源極發(fā)射電能,直接為分布于源極機(jī)構(gòu)一定尺度范圍內(nèi)的單個(gè)或多個(gè)受電設(shè)備供電.然而,ICPT的電能傳輸效率與傳輸距離之間的矛盾一直是阻礙該項(xiàng)技術(shù)推廣應(yīng)用的瓶頸,為確保距離電源供給點(diǎn)較遠(yuǎn)的拾取機(jī)構(gòu)能獲取電能,通常采用提高傳輸頻率的辦法來加以解決.通過提升發(fā)射磁場頻率固然可以適當(dāng)緩解這一矛盾,但隨著頻率的增大(大尺度傳輸?shù)念l率在GHz級,中尺度的在MHz級)[14],能量在外界空間的耗散也將急劇增加,且電磁干擾(EMI)將隨之增強(qiáng).同時(shí),隨著傳輸距離的增加,傳輸效率將急劇下降,甚至可導(dǎo)致拾取電能失效.因此,以“廣播式”供電模式實(shí)現(xiàn)中、大尺度ICPT在實(shí)際應(yīng)用中存在著巨大的挑戰(zhàn).另外,對于一定空間范圍內(nèi)任意分布的多用電設(shè)備的供電問題,由于用電設(shè)備位置的隨機(jī)性與時(shí)變性、電能容量以及耗能水平的差異性,采用“廣播式”供電顯然存在空間限制、實(shí)時(shí)性差等局限.但是一定距離的無線電能傳輸對于分布于空間范圍內(nèi)的多移動(dòng)設(shè)備實(shí)時(shí)供電或間歇充電是極其重要的,例如機(jī)器人足球賽中位于場地中的多個(gè)移動(dòng)機(jī)器人的電能實(shí)時(shí)動(dòng)態(tài)補(bǔ)給、無線傳感器網(wǎng)絡(luò)的節(jié)點(diǎn)能量保障、分布式衛(wèi)星群能量保障等.
基于以上分析,針對空間上任意分布的多負(fù)載因ICPT傳輸空間尺度限制無法實(shí)時(shí)動(dòng)態(tài)供電和電能傳輸效率低下的問題,文中提出一種無線電能傳輸網(wǎng)架構(gòu).基于此架構(gòu),分析了無線電能傳輸網(wǎng)的電能傳輸效率特性和電能傳輸路徑優(yōu)化模型,并提出一種基于改進(jìn)型蟻群算法的動(dòng)態(tài)無線電能傳輸組網(wǎng)方法,以解決電能傳輸中的動(dòng)態(tài)路由問題,并通過仿真驗(yàn)證該算法的有效性.
在文中提出的無線電能傳輸網(wǎng)中,用電設(shè)備作為網(wǎng)絡(luò)節(jié)點(diǎn),可分為具有3種功能的節(jié)點(diǎn),即源節(jié)點(diǎn)、受電節(jié)點(diǎn)、中繼節(jié)點(diǎn).源節(jié)點(diǎn)作為電源發(fā)射電能,受電節(jié)點(diǎn)作為負(fù)載接收電能,中繼節(jié)點(diǎn)完成能量傳輸?shù)闹欣^任務(wù).通過中繼節(jié)點(diǎn)的“接力”作用構(gòu)成一條源節(jié)點(diǎn)與受電節(jié)點(diǎn)之間的由一系列無線移動(dòng)節(jié)點(diǎn)組成的電能傳輸通道,實(shí)現(xiàn)多跳鏈路電能傳輸.若干條傳輸通道集合便形成一種網(wǎng)內(nèi)無中心、節(jié)點(diǎn)平等的無線電能傳輸網(wǎng),如圖1所示.所有節(jié)點(diǎn)均具有雙向無線電能傳輸功能[8],既可接收電能也可傳輸電能.
圖1 無線電能傳輸網(wǎng)的物理拓?fù)銯ig.1 Physical topology of wireless power transfer network
圖2 無線電能傳輸過程拓?fù)銯ig.2 Topology of wireless power transfer
如圖2所示,受電節(jié)點(diǎn)根據(jù)自身的能量需求向網(wǎng)絡(luò)中的節(jié)點(diǎn)提出能量供給請求,該請求經(jīng)Zigbee通信方式向網(wǎng)絡(luò)中鄰近節(jié)點(diǎn)發(fā)出(過程①).鄰近節(jié)點(diǎn)將自身狀態(tài)信息(包括能量狀況、與受電節(jié)點(diǎn)的相對位置)向受電節(jié)點(diǎn)回饋(過程②).受電節(jié)點(diǎn)根據(jù)鄰近節(jié)點(diǎn)的狀態(tài)信息,動(dòng)態(tài)決定當(dāng)前電能傳輸鏈路,并向相關(guān)節(jié)點(diǎn)下達(dá)傳輸令牌(過程③),完成無線電能傳輸過程.
把無線電能傳輸網(wǎng)抽象為有權(quán)圖C(V,E),如圖3所示.V(C)為圖C的節(jié)點(diǎn)集,E(C)為圖C的邊集,元素v∈V為圖 C的一個(gè)節(jié)點(diǎn),元素 eij(eij∈E,eij=vivj)為圖C的一條從vi到vj的邊.令元素e∈E具有一組有序數(shù)列(wp,we)作為e的屬性,稱作弧的度量,wp為兩端的剩余電能,we為端到端的電能傳輸效率.為了便于分析,在研究中通常使用對稱網(wǎng)絡(luò)模型以減少弧的數(shù)量,即有eij=eji.在圖C中,如果 vi,vi+1∈E,i=1,2,…,d,則稱 P=(v1,vd)為圖 C的一條從v1到vd的電能傳輸路徑.
圖3 有權(quán)圖Fig.3 Weighted chart
在本模型中wp為可加性度量,其值越大,整條傳輸路徑抗毀性越強(qiáng),we為可乘性度量.對于路徑P=(v1,vd),有
路徑P=(v1,vd)上的多目標(biāo)優(yōu)化函數(shù)為
約束條件為
式中,γp、γe分別為剩余電能和傳輸效率的權(quán)重因子,m為每次迭代過程中搜索到目標(biāo)節(jié)點(diǎn)的螞蟻個(gè)數(shù),Pmax為某條路徑上的最大剩余電能.若該函數(shù)值越小,則電能傳輸路徑越優(yōu)越.因此,無線電能傳輸網(wǎng)多目標(biāo)電能傳輸路徑優(yōu)化問題可以描述為:在圖C中,搜尋一條電能傳輸路徑P=(v1,vd),使得該路徑上的目標(biāo)函數(shù)值最小.
文中采用蟻群優(yōu)化思想實(shí)現(xiàn)無線電能傳輸網(wǎng)混合多目標(biāo)電能傳輸路徑尋優(yōu).
蟻群算法為對自然界蟻群覓食尋徑方式進(jìn)行數(shù)學(xué)模擬而得出的一種仿生算法.在t時(shí)刻,某只螞蟻k從節(jié)點(diǎn)vi轉(zhuǎn)移到vj的轉(zhuǎn)移概率函數(shù)為
式中:τij(t)為t時(shí)刻在路徑(vi,vj)上殘留的信息素;ηij(t)為 t時(shí)刻在路徑(vi,vj)上的啟發(fā)因子;τiv(t)為t時(shí)刻螞蟻k在路徑(vi,v)上殘留的信息素;ηiv(t)為t時(shí)刻螞蟻k在路徑(vi,v)上的啟發(fā)因子;α、β為路由選擇中路徑上信息素殘留和啟發(fā)因子的權(quán)重;因?yàn)閳D搜索蟻群算法具有記憶功能,tabuk(k=1,2,…,m)稱為禁忌表,用以保存螞蟻 k已走過的節(jié)點(diǎn)集合.螞蟻在經(jīng)歷Δt時(shí)間后,從源節(jié)點(diǎn)v1到達(dá)目標(biāo)節(jié)點(diǎn)vd,所經(jīng)過的路徑采用全局狀態(tài)信息素更新策略:
式中,ρ為更新信息素?fù)]發(fā)系數(shù),改變其大小可調(diào)整路徑信息素增長速度.
在Δt時(shí)間內(nèi),設(shè)每次迭代過程中共有m只螞蟻搜索到了較優(yōu)路徑,則路徑上信息素增量可通過式(7)-(9)求得:
式中,αk為第k較優(yōu)螞蟻在路徑P上信息素增加的權(quán)重,K為常數(shù)為螞蟻k剩余電能.定義局部信息素更新規(guī)則和局部信息素啟發(fā)因子:
為抑制蟻群算法早熟收斂,使得算法在初始階段能夠?qū)尚杏蜻M(jìn)行全面搜索,文中采用動(dòng)態(tài)更新機(jī)制對信息素更新[15]
圖4 c、m、γ 的關(guān)系Fig.4 Relationship among c,m and γ
為提升無線電能傳輸網(wǎng)節(jié)點(diǎn)的電磁能空間傳輸能力,文中擬采用基于磁共振耦合模式的電能傳輸方式,系統(tǒng)的電能傳輸原理框圖如圖5所示.
圖5 磁共振耦合電能傳輸原理框圖Fig.5 Basic structure of magnetic resonant-coupling power transmission
磁共振耦合電能傳輸系統(tǒng)由4個(gè)線圈組成,即激發(fā)線圈、特征形狀及特性參數(shù)均一致的兩個(gè)共振線圈,共振線圈具有完全相同的共振頻率、接收線圈.激發(fā)線圈與共振線圈之間通過傳統(tǒng)的電磁感應(yīng)耦合方式將電能傳輸?shù)焦舱窬€圈;根據(jù)共振原理,如果兩個(gè)振蕩電路具有相同頻率,則在波長范圍內(nèi),通過近場瞬逝波耦合.感應(yīng)器產(chǎn)生的駐波在遠(yuǎn)遠(yuǎn)小于損耗時(shí)間的周期內(nèi),允許能量高效地從一個(gè)物體傳到另一物體.兩個(gè)共振線圈以相同的共振頻率振蕩,磁場會(huì)通過類似于“隧道”的傳能形式將電能傳輸?shù)骄哂邢嗤C振頻率的共振線圈.接收線圈與共振線圈之間以傳統(tǒng)的感應(yīng)耦合方式完成電能傳輸,接收線圈獲得的電能經(jīng)過功率調(diào)節(jié)后提供給負(fù)載.為實(shí)現(xiàn)雙向無線電能傳輸,電能發(fā)射端與接收端在拓?fù)渖贤耆珜Φ?,功能互補(bǔ).
文中以原、副邊均為串聯(lián)補(bǔ)償方式分析磁場共振耦合電能傳輸系統(tǒng),其電路拓?fù)浣Y(jié)構(gòu)如圖6所示.以電磁感應(yīng)耦合理論分析系統(tǒng)的數(shù)學(xué)模型,不失一般性,磁場共振耦合電能傳輸系統(tǒng)可以集中電路元素(L、C、R)來表示,4個(gè)線圈之間的交叉耦合作用很弱,所以,在以下的分析過程忽略交叉耦合.
圖6 磁共振耦合拓?fù)浣Y(jié)構(gòu)Fig.6 Topology of magnetic resonant coupling
基于Kirchhoff電壓定理可得4個(gè)回路的電壓:
假定電壓源逆變電路工作在零電流軟開關(guān)(ZCS)模式下,Vp-v可用式(14)計(jì)算[4].
式中:κij為線圈之間的耦合系數(shù).通過方程(13)-(15)可求出負(fù)載上的電壓表達(dá)式.激發(fā)線圈與共振線圈通過耦合系數(shù)κ12連接,共振線圈與接收線圈通過κ34連接.在系統(tǒng)中,一般地,κ12與 κ34是固定不變的.κ23是兩個(gè)共振線圈之間距離的函數(shù),隨著距離的改變而變化.為簡化分析,忽略互感系數(shù)κ13、κ24、κ14的影響.ω為系統(tǒng)工作角頻率.Mij為線圈Li與Lj之間的互感,Li為線圈i的自感.系統(tǒng)要實(shí)現(xiàn)磁場共振耦合,共振頻率ω0可由式(16)確定:
用Zij表示線圈i反射到線圈j的反射等效阻抗.諧振時(shí),根據(jù)文獻(xiàn)[16]可得Zij的計(jì)算公式如下.
從接收線圈和共振線圈L2、L3反射到激發(fā)線圈的等效反射阻抗Zr為
在共振的條件下,由于共振線圈的共振頻率相等,滿足
因此,在共振頻率下,假設(shè)共振線圈品質(zhì)因素Q很高,忽略共振線圈內(nèi)阻,反射阻抗Zr為
由式(20)可知,反射阻抗Zr類似于傳統(tǒng)ICPT的反射阻抗,只是耦合系數(shù)有所變化而已.共振耦合系統(tǒng)可等效為如圖7所示的傳統(tǒng)感應(yīng)耦合系統(tǒng).
圖7 ICPT網(wǎng)絡(luò)的拓?fù)浣Y(jié)構(gòu)Fig.7 Topology of ICPT system
設(shè)定激發(fā)線圈與接收線圈的自感相等,即
則Zr的計(jì)算公式為
分析式(22)可知,理論上通過調(diào)節(jié)系統(tǒng)4個(gè)線圈之間的距離,系統(tǒng)能獲得很高的耦合效率.
忽略接收線圈內(nèi)阻和共振線圈內(nèi)阻時(shí),在共振頻率下,磁共振耦合系統(tǒng)的傳輸功率與傳輸效率分別為
由上式可知,磁共振耦合電能傳輸系統(tǒng)的傳輸效率隨著耦合系數(shù)的減小呈指數(shù)下降趨勢.
仿真時(shí)所用的算法參數(shù)設(shè)置如表1所示,為驗(yàn)證所提無線電能傳輸網(wǎng)組網(wǎng)算法的動(dòng)態(tài)性與有效性,在兩種不同的仿真試驗(yàn)條件下進(jìn)行比較分析:一種情況為最大無線電能傳輸距離為1m,另外一種情況為最大無線電能傳輸距離為2m.
文中擬以機(jī)器人足球賽中電能耗盡的機(jī)器人向場內(nèi)其它機(jī)器人請求電能供給為例來論證所提出的無線電能傳輸組網(wǎng)算法.假定網(wǎng)絡(luò)拓?fù)浣Y(jié)構(gòu)如圖1所示,通過節(jié)點(diǎn)的無線通信模塊可得個(gè)體剩余電能和個(gè)體與其鄰居節(jié)點(diǎn)的距離,由第3節(jié)分析可知節(jié)點(diǎn)間的能量傳輸效率.由于場內(nèi)個(gè)體所剩電能不同且動(dòng)態(tài)變化,同時(shí)無線電能傳輸效率隨著距離變大而降低,因此以某一時(shí)刻下節(jié)點(diǎn)v1請求無線電能供給,與節(jié)點(diǎn)v1相距較遠(yuǎn)的節(jié)點(diǎn)v11作為無線電能補(bǔ)給點(diǎn),建立無線能量多跳路由.如無中繼節(jié)點(diǎn)的接力作用,則節(jié)點(diǎn)v1和節(jié)點(diǎn)v11之間無線傳能失效.基于以上仿真環(huán)境運(yùn)行組網(wǎng)算法,結(jié)果分別如表2、圖8(a)、8(b)所示.
表1 仿真參數(shù)1)Table 1 Simulation parameters
表2 仿真結(jié)果比較Table 2 Comparison of simulation results
圖8 不同最大傳輸距離時(shí)的迭代過程Fig.8 Process of the iterative optimization at different maximum transmission distances
由表2可知,不同最大傳輸距離下所搜索到的電能傳輸鏈路有所差異,改進(jìn)蟻群算法平均搜索用時(shí)要低于蟻群算法,更易找到最優(yōu)解.
縱坐標(biāo)的目標(biāo)函數(shù)值為在每次迭代過程中,所有找到目標(biāo)節(jié)點(diǎn)的螞蟻中最優(yōu)螞蟻的目標(biāo)函數(shù)值.由圖可知,蟻群算法在迭代15次時(shí)收斂,改進(jìn)蟻群算法在迭代次數(shù)約為13次時(shí)收斂,因此改進(jìn)蟻群算法搜尋到最優(yōu)電能傳輸鏈路的計(jì)算效率要高于蟻群算法,減少了尋優(yōu)時(shí)間.
由于傳統(tǒng)的ICPT系統(tǒng)不能有效地解決電能傳輸距離和系統(tǒng)整體電能傳輸效率之間的突出矛盾,文中提出無線電能傳輸網(wǎng)的概念,用于實(shí)現(xiàn)中尺度范圍內(nèi)電能無線傳輸.文章分析并建立了磁共振耦合電能傳輸?shù)南到y(tǒng)模型,研究了系統(tǒng)的傳輸效率特性.采用動(dòng)態(tài)更新機(jī)制對信息素進(jìn)行更新,以增強(qiáng)蟻群算法的全局搜索能力.基于改進(jìn)蟻群算法的無線電能傳輸網(wǎng)組網(wǎng)方法可以實(shí)現(xiàn)效率最大化且鏈路剩余電能最多的最優(yōu)電能傳輸路徑,以達(dá)到延伸電能傳輸距離和保證電能傳輸路徑穩(wěn)定的目的.仿真結(jié)果驗(yàn)證了改進(jìn)蟻群算法在無線電能傳輸網(wǎng)組網(wǎng)方面的優(yōu)越性.
[1]Boys J T,Green A W.Inductively coupled power transmission concept-design and application [J].IPENZ Transactions,1995,22(1):1-9.
[2]Green A W,Boys J T.10 kHz inductively coupled power transfer-concept and control[J].Power Electronics and Variable Speed Drives,1994,399:694-699.
[3]Tang Chun-sen,Sun Yue,Su Yu-gang,et al.Determining multiple steady-state ZCS operating points of a switchmode contactless power transfer system [J].IEEE Transactions on Power Electronics,2009,24(2):416-425.
[4]Li H L,Hu A P,Covic G A,et al.Optimal coupling condition of IPT system for achieving maximum power transfer[J].Electronics Letters,2009,45(1):76-77.
[5]Wu H H,Hu A P,Si P,et al.A push-pull resonant converter with dual coils for transcutaneous energy transfer systems[C]∥Industrial Electronics and Applications.Xi'an:IEEE,2009.
[6]Si P,Hu A P,Malpas S,et al.A frequency control method for regulating wireless power to implantable devices[J].IEEE Transactions on Biomedical Circuits and Systems,2008,2(1):22-29.
[7]陳敏,周鄧燕,徐德鴻.注入高次諧波電流的磁懸浮列車非接觸供電方法[J].中國電機(jī)工程學(xué)報(bào),2005,25(6):104-108.Chen Min,Zhou Deng-yan,Xu De-hong.Contactless power supply of maglev using harmonic injection method[J].Proceedings of the CSEE,2005,25(6):104-108.
[8]戴欣,孫躍,蘇玉剛,等.非接觸電能雙向推送模式研究 [J].中國電機(jī)工程學(xué)報(bào),2010,30(18):55-61.Dai Xin,Sun Yue,Su Yu-gang,et al.Study on contactless power bi-directional push mode[J].Proceedings of the CSEE,2010,30(18):55-61.
[9]王璐,陳敏,徐德鴻.磁懸浮列車非接觸緊急供電系統(tǒng)的工程化設(shè)計(jì)[J].中國電機(jī)工程學(xué)報(bào),2007,27(18):67-70.Wang Lu,Chen Min,Xu De-hong.The engineering design of contactless emergency power supply in maglev[J].Proceedings of the CSEE,2007,27(18):67-70.
[10]Kissin M L G,Covic G A,Boys J T.Estimating the output power of flat pickups in complex IPT systems[C]∥Proceedings of IEEE PESC.Island of Rhodes:IEEE,2008.
[11]Zhou Wen-qi,Ma Hao.Winding structure and circuit design of contactless power transfer platform[C]∥Proceedings of IEEE IECON.Orlando:IEEE,2008.
[12]Liu X,Hui S Y.Optimal design of a hybrid winding structure for planar contactless battery charging platform[J].IEEE Transactions on Power Electronics,2008,23(1):455-463.
[13]Kurs A,Karatis A,Moffatt R,et al.Wireless power transfer via strongly coupled magnetic resonances[J].Science,2007,317(7):83-85.
[14]Karalis A,Joannopoulos J D,Soljacic M.Efficient wireless non-radiative mid-range energy transfer[J].Annals of Physics,2008,323(1):34-48.
[15]Song Y H,Irving M R.An overview of heuristic optimization techniques for power system expansion planning and design[J].Power Engineering Society General Meeting,2004,1(6):933.
[16]Wang Chwei-Sen,Stielau Oskar H,Grant A.Covic design considerations for a contactless electric vehicle battery charger[J].IEEE Transactions on Industrial Electronics,2005,52(5):1308-1314.