車用鋰電池組熱流場(chǎng)特性數(shù)值模擬與優(yōu)化設(shè)計(jì)

韓鋒鋼， 潘懷杰， 彭倩， 林有淮， 劉勝

(1. 廈門理工學(xué)院 機(jī)械與汽車工程學(xué)院， 福建 廈門 361024；2. 廈門威迪思汽車設(shè)計(jì)服務(wù)有限公司， 福建 廈門 361024)

車用鋰電池組熱流場(chǎng)特性數(shù)值模擬與優(yōu)化設(shè)計(jì)

韓鋒鋼1，2， 潘懷杰1， 彭倩1，2， 林有淮2， 劉勝1

(1. 廈門理工學(xué)院 機(jī)械與汽車工程學(xué)院， 福建 廈門 361024；2. 廈門威迪思汽車設(shè)計(jì)服務(wù)有限公司， 福建 廈門 361024)

針對(duì)車用電池溫升過高、電池組溫差大的問題，開展電池包熱流場(chǎng)分析與優(yōu)化設(shè)計(jì).根據(jù)Bernardi的生熱速率方程式，建立由電池電解液、正負(fù)極柱和隔膜四部分組成的單體電池?zé)狁詈夏Ｐ图俺山M電池傳熱模型；利用Fluent軟件分析鋰電池單體在自然對(duì)流環(huán)境下的溫升特性，研究成組電池在強(qiáng)制對(duì)流條件下的熱流場(chǎng)特性；通過增加導(dǎo)流板優(yōu)化電池箱內(nèi)流場(chǎng)結(jié)構(gòu)，并評(píng)估導(dǎo)流板對(duì)電池組散熱效率的作用.結(jié)果表明：?jiǎn)误w鋰電池在自然對(duì)流下溫升明顯，電池內(nèi)核溫度遠(yuǎn)高于正負(fù)極柱溫度；電池箱進(jìn)出風(fēng)口位置及結(jié)構(gòu)決定箱內(nèi)空氣的流向和成組電池的散熱效果；通過對(duì)進(jìn)、出風(fēng)口位置的設(shè)計(jì)及增加導(dǎo)流板，可有效改進(jìn)電池組熱流場(chǎng)的均勻性，從而提高散熱效果.

鋰離子電池組； 新能源汽車； 熱模型； 風(fēng)冷； 導(dǎo)流板； 溫度場(chǎng)

Abstract： In order to solve the problem of temperature sharp rising and significant thermal difference among the battery modules of electric vehicles， a study of thermal flow analysis and optimization on battery pack was carried out. A thermal coupling model with four parts, namely a single cell electrolyte， a positive plate collector, a negative plate collector and a battery plate was established. A heating transferring model of individual battery was also established based on Bernardi heat generation rate model; then using Fluent software the heat generation characteristics of lithium battery monomer in natural convection environment was investigated and the effect of flow structure on the heat radiation was assessed with an optimized battery pack model. The re-sults indicate that the individual lithium battery temperature rises sharply under natural convection， while temperature increment of positive pole core is far higher than that of negative pole core； the inlet structure greatly influences the velocity field of cooling air and temperature distribution of the battery； and the uniformity of the battery thermal flow can be significantly improved by reasonably designing the air inlet and outlet positions and adding a guide plate, which can further enhance the cooling effect.

Keywords： lithium-ion battery pack； new energy vehicles； thermal model； air cooling； guide plate； temperature field

由于與新能源電池性能有關(guān)系的相關(guān)參數(shù)，如工作電壓、放電容量、循環(huán)壽命等都與溫度息息相關(guān)[1-2],電池組充放電過程中電流大、箱內(nèi)布置密集緊湊且工作環(huán)境存在較大不確定性，電池組的工作效率、穩(wěn)定性和使用壽命會(huì)由于箱內(nèi)熱量耗散的問題受到嚴(yán)重制約.當(dāng)電池溫度不在40 ℃到20 ℃的范圍內(nèi)時(shí)，磷酸鐵鋰離子動(dòng)力電池組(LiFePO4)的放電效率會(huì)有一個(gè)快速下降的表現(xiàn)，若要使其放電效率達(dá)到80%以上，電池組理想的工作溫度最好是處于18～43 ℃之間[3].國(guó)內(nèi)眾多的學(xué)者針對(duì)單體電池生熱模型及電池箱體散熱效果等問題開展了研究.李騰等[4]介紹了在鋰離子新能源電池方面所開展的電-熱及電化學(xué)-熱耦合模型等的研究；陳燕虹等[5]采取數(shù)值模擬的方式研究了成組電池生熱和散熱特性；羅玉濤等[6]則提出了以變壓器油作為冷卻介質(zhì)直接接觸液體的另一種冷卻方式.但目前大部分的研究都集中在散熱方式、相變材料的研究上，而在風(fēng)冷方式的研究中，還缺乏對(duì)電池箱內(nèi)流場(chǎng)起導(dǎo)向作用的進(jìn)、出風(fēng)口的設(shè)計(jì).本文對(duì)車用鋰電池組熱流場(chǎng)特性進(jìn)行了數(shù)值模擬與改進(jìn)設(shè)計(jì)，并通過進(jìn)、出風(fēng)口位置設(shè)計(jì)及增加導(dǎo)流板優(yōu)化電池箱體內(nèi)流場(chǎng)走向，以改善成組電池?zé)崃鲌?chǎng)的散熱性能.

1 數(shù)學(xué)模型

由于電池箱內(nèi)緊密布置結(jié)構(gòu)件及電池組，箱內(nèi)氣流為湍流流動(dòng)，其換熱形式主要為強(qiáng)迫對(duì)流換熱及固體間熱傳導(dǎo).

文中在建立流體仿真計(jì)算的數(shù)學(xué)模型時(shí)采用k-ε湍流模型(k為湍動(dòng)能，ε為耗散率)[7-8]，該模型包括5個(gè)主要方程.

1) 質(zhì)量守恒方程，其表達(dá)式為

式(1)中:vr，vz，vφ分別為圓柱坐標(biāo)系下徑向、軸向和周向的相對(duì)流體速度；ρ為氣體密度；t為時(shí)間.

2) 動(dòng)量方程(忽略氣體重力作用)為

式(2)中：V為流體速度向量；μeff為等效粘度，μeff=μ+μt，其中,μ為流體粘度，μt為湍流粘度，μt=Cμρk2/ε，Cμ為常數(shù)，取為0.09；pt為修正壓力，pt=p+2ρk/3，p為壓力.

3) 能量方程為

式(3)中：ht為氣體總焓；λ為導(dǎo)熱系數(shù)；T為氣體溫度.

4) 湍動(dòng)能方程(k方程)為

5) 湍動(dòng)耗散率方程(ε方程)為

式(5)中：Cε1，Cε2，σk，σε均為常數(shù).根據(jù)Launder等[9]的研究結(jié)果及后來相關(guān)的試驗(yàn)，在標(biāo)準(zhǔn)k-ε模型中，Cε1，Cε2，σk，σε分別取Cε1=1.44，Cε2=1.92，σk=1.0，σε=1.3.

2 鋰離子電池生熱機(jī)理

充放電過程中，電池內(nèi)數(shù)目相等的鋰離子和電子會(huì)有嵌入和脫嵌的現(xiàn)象，該動(dòng)作產(chǎn)生的熱量稱之為反應(yīng)熱.充放電中的該過程為可逆反應(yīng)，大小相等符號(hào)相反[10]，記為Qr；根據(jù)熵增原理，實(shí)際上該動(dòng)作中，還包括極化反應(yīng)熱Qp、過充及過放副反應(yīng)熱Qs、電解質(zhì)分解和由于自放電生成的熱量Qj，以及電池內(nèi)阻焦耳熱Qt.因此，上述關(guān)系簡(jiǎn)化為Qt=Qr+Qp+Qs+Qj.

因?yàn)锽MS對(duì)每個(gè)電池的電流、電壓和溫度數(shù)據(jù)都會(huì)進(jìn)行檢測(cè)，防止過充過放電，所以Qs作用很小，可以忽略.此外，可以利用等效極化內(nèi)阻產(chǎn)生的熱量來代替極化反應(yīng)熱[11]，因此，該過程可簡(jiǎn)化為

式(6)中：R=Rz+Rp，Rz是電池自身的歐姆內(nèi)阻，Rp是電池的等效極化內(nèi)阻.

3 數(shù)值模型建立

3.1鋰離子電池生熱模型

在電池包的熱行為仿真計(jì)算中，將電池單體視為一個(gè)單獨(dú)的封閉子系統(tǒng)，和外界僅有熱量方面的交換.假設(shè)單電池內(nèi)部的材料是相同的材料，比熱容和X，Y，Z三個(gè)方向的熱導(dǎo)率等于電池的導(dǎo)熱系數(shù)，充電狀態(tài)和實(shí)時(shí)溫度不會(huì)對(duì)其產(chǎn)生影響，則成組電池生熱率模型為

運(yùn)行中的電動(dòng)汽車電池組的發(fā)熱遵守?zé)崃W(xué)平衡方程，即

通常，電池組中單體質(zhì)量和材料都是完全相同的，因此，電池組吸收熱量計(jì)算式可簡(jiǎn)化為

式(9)中：cp為電池的平均比熱容；m為電池平均質(zhì)量；ΔTi為單體電池的溫差.

當(dāng)單體電池以不同倍率放電時(shí)，由于材料不同導(dǎo)致電池的內(nèi)核與極柱生熱的速率是有差異的.當(dāng)生熱速率比散熱速率大時(shí)，會(huì)導(dǎo)致電池的溫度持續(xù)上升；溫度過高將影響電池的循環(huán)使用壽命和安全性.因此，計(jì)算散熱率和散熱速率對(duì)各種工況下的仿真計(jì)算都有著重要的意義[12].為了便于數(shù)值模擬，文中只針對(duì)0.5C的放電倍率恒流放電的情況進(jìn)行仿真，如表1所示.表1中：C為放電倍率；I為單體電池電流；Q為生熱率.

表1 不同放電倍率下單體電池各部分生熱率Tab.1 Heat generation rate of each part of single cell at different discharge rate

3.2鋰離子電池?zé)崽匦詤?shù)

采用的鋰離子單體電池為方體，尺寸為140 mm×64 mm×211 mm，鋰單體電池劃分為電解液區(qū)域、正負(fù)極柱和隔膜四部分，每個(gè)區(qū)域的等效參數(shù)，如表2所示.表2中：h為材料厚度；ρ為材料密度；正極、負(fù)極及隔膜的參數(shù)為浸滿電解液工況下的參考值，由實(shí)驗(yàn)獲得;正極、負(fù)極、隔膜及電解液的比熱容(c)和熱導(dǎo)率(λ)值可由文獻(xiàn)[2]得到.

表2 鋰離子電池材料特征參數(shù)Tab.2 Characteristic parameters of lithium ion battery materials

圖1 鋰單體電池有限元模型Fig.1 Finite element model of lithium battery

4 車用電池箱熱流場(chǎng)分析

4.1單體電池?zé)崃鲌?chǎng)分析

基于電池的幾何參數(shù)，利用CREO建模軟件對(duì)單體電池進(jìn)行幾何建模并導(dǎo)入到Hypermesh軟件中進(jìn)行流體網(wǎng)格的劃分，如圖1所示.采用流體模擬軟件Fluent對(duì)單體電池的熱耦合模型進(jìn)行求解，得到單體電池的溫度分布.根據(jù)相關(guān)文獻(xiàn)，對(duì)于空氣中自然散熱的電池表面，在計(jì)算時(shí)取其極耳傳熱系數(shù)為3.9 W·(m·K)-1，極板表面?zhèn)鳠嵯禂?shù)為0.427 W·(m·K)-1.電池單體在自然冷卻，0.5C和1.0C放電倍率下的恒流放電溫度云圖，如圖2，3所示.由圖2，3可知：0.5C放電倍率時(shí),單體電池正、負(fù)極柱最高溫度分別為36.4，35.7 ℃；而1.0C放電倍率時(shí)則為37.1，36.3 ℃.說明正極柱溫差高于負(fù)極柱溫差，單體電池溫度由質(zhì)心向四周殼體逐級(jí)遞減.由于極柱傳熱系數(shù)大于電池殼體的傳熱系數(shù)，因此,最大溫度出現(xiàn)在電池中心偏上的位置，在自然冷卻環(huán)境中，單體電池最大溫差達(dá)到13.3 ℃，說明電池的熱積累是電動(dòng)汽車使用的關(guān)鍵問題.

(a) Z-X對(duì)稱面 (b) 單體電池整體 (a) Z-X對(duì)稱面 (b) 單體電池整體 圖2 單體電池0.5C放電倍率的溫度云圖 圖3 單體電池1.0C放電倍率的溫度云圖 Fig.2 Temperature nephogram of single battery at 0.5C charge Fig.3 Temperature nephogram of single battery at 1.0C charge

圖4 不同放電倍率下電池組的溫度曲線Fig.4 Temperature curve of battery pack at different discharge rate

分別監(jiān)測(cè)單體電池Z-X對(duì)稱面自上向下的5個(gè)等距離監(jiān)測(cè)點(diǎn)的溫度，結(jié)果如圖4所示.

4.2電池組熱流場(chǎng)分析

當(dāng)汽車以較高的速度運(yùn)行，如大角度爬坡時(shí)，電池組的放電倍率是最大的，也是電池生熱量最大的工況.若僅靠箱體自然冷卻會(huì)造成熱量積累，特別是電池組中心聚集的熱量無法很快地耗散，過熱的溫度不僅使電池性能下降，而且造成的損傷是無法修復(fù)的，使電池容量和性能下降.因此,必須采取強(qiáng)制冷卻的方式，使電池箱內(nèi)的溫度保持在能使鋰離子電池正常工作的范圍內(nèi).以某車用鋰電池組為例，該電池箱內(nèi)置電池由4列5排共20個(gè)電池單體單層排列構(gòu)成，大小為953.5 mm×556.5 mm×211.0 mm.圖5為冷卻風(fēng)扇及電池組布置圖,圖6為其有限元模型.用于送風(fēng)的風(fēng)扇有兩個(gè)，置于距離電池組前端50 mm，出風(fēng)口有一個(gè)，置于距離電池組25 mm處.

數(shù)值模擬邊界條件采用速度入口和壓力出口，箱體其余外部邊界采用默認(rèn)壁面邊界，空氣入口速度為2.2 m·s-1.電池箱體與外部空氣的能量交換主要途徑是對(duì)流和輻射換熱，由于熱管理系統(tǒng)的存在，可忽略其對(duì)成組電池?zé)崃鲌?chǎng)特性的影響.參考文獻(xiàn)[13]可得對(duì)流換熱系數(shù)hx=15 W·(m2·K)-1，輻射換熱系數(shù)hr=0.2.在外部環(huán)境攝氏溫度為25 ℃情況下，以0.5C放電倍率下,恒流放電完成后箱體內(nèi)冷卻空氣流線，如圖7所示.

圖5 冷卻風(fēng)扇及電池組布置圖 圖6 電池組網(wǎng)格模型 Fig.5 Diagram of cooling fan and battery pack Fig.6 Meshing schematic diagram of battery pack

由圖7(a)可知:進(jìn)風(fēng)端空氣流速(v)較出風(fēng)端低，冷卻空氣流速由進(jìn)口端沿出口端均勻平順，空氣流動(dòng)能量損失較小.由圖7(b)可知：在進(jìn)口端與電池組之間存在多個(gè)漩渦，造成冷卻空氣能量的損失，當(dāng)冷卻空氣到達(dá)電池箱上端面后，有一部分空氣經(jīng)單體電池之間的縫隙流向箱體底部，隨后往電池箱出風(fēng)口流動(dòng)；各排之間空氣流速大小及均勻性相差不大，而各列之間空氣流速大小及均勻性相差較大.

優(yōu)化前的電池組溫度分布，如圖8所示.由圖8可知：第3，4排電池組溫度較高，結(jié)合圖7可知此處的空氣是流速最低的地方；溫度較低處出現(xiàn)在進(jìn)風(fēng)口最前排頂部和出風(fēng)口第5排中間區(qū)域；最高溫度44 ℃出現(xiàn)在第3排中下部位置，最低溫度30 ℃出現(xiàn)在第1排進(jìn)風(fēng)口處，成組電池的溫差最大為14 ℃，可見電池組熱量分布均勻性不理想.

(a) 軸視圖 (b) 流場(chǎng)云圖圖7 電池箱冷卻空氣流場(chǎng)云圖Fig.7 Velocity streamlines of cooling air in battery box

(a) 電池組溫度云圖 (b) 電池組溫度曲線圖8 優(yōu)化前的電池組溫度分布Fig.8 Temperature distribution of battery pack before optimization

4.3電池箱結(jié)構(gòu)優(yōu)化分析

針對(duì)上述仿真結(jié)果存在的不足之處對(duì)電池箱結(jié)構(gòu)進(jìn)行改進(jìn)設(shè)計(jì)，在進(jìn)風(fēng)口與電池組之間增加導(dǎo)流板，其結(jié)構(gòu)如圖9，10所示.為了加大對(duì)電池組底部區(qū)域的冷卻效果，將出風(fēng)口向電池箱底部下移30 mm，其他布局與原結(jié)構(gòu)一致.優(yōu)化后，以0.5C放電倍率時(shí)的冷卻空氣速度流線云圖，如圖11所示；而以0.5C和1.0C放電倍率時(shí)的電池組溫度分布，如圖12,13所示.

由圖12，13可知：該電池組最高和最低溫度分別為35.2，26.8 ℃，均在單體電池理想的工作溫度范圍內(nèi)，并且將電池組的溫差縮小至8.4 ℃.為滿足使用壽命及安全性，要求留出了足夠的余量.

圖9 導(dǎo)流板結(jié)構(gòu)示意圖 圖10 導(dǎo)流板位置Fig.9 Schematic diagram of guide plate Fig.10 Guide plate position

(a) 俯視圖 (b) 流線云圖圖11 優(yōu)化后以0.5C放電倍率時(shí)的冷卻空氣速度流線云圖Fig.11 Velocity streamline of cooling air at 0.5C discharge rate after optimization

(a) 0.5C放電倍率 (b) 1.0C放電倍率圖12 優(yōu)化后不同放電倍率下的電池組溫度云圖Fig.12 Optimized battery temperature at different discharge magnifications

(a) 0.5C放電倍率 (b) 1.0C放電倍率圖13 優(yōu)化后不同放電倍率下的電池組溫度曲線Fig.13 Optimized battery temperature curves at different discharge magnifications

4.4不同放電倍率下的電池組溫度

以改進(jìn)后電池箱結(jié)構(gòu)為基準(zhǔn)，模擬計(jì)算了不同放電倍率下電池組的溫度場(chǎng)，結(jié)果如表3所示.由表3可知： 隨著運(yùn)行工況不同導(dǎo)致放電倍率逐漸增大的工況下， 電池組的溫升和最大溫差都在逐漸增加，成組電池組在1.0C恒流放電時(shí)的溫差最高可達(dá)到38.56 ℃，接近鋰離子電池最舒適的工作環(huán)境溫度(40 ℃).由此可見，當(dāng)電動(dòng)汽車長(zhǎng)時(shí)間在復(fù)雜環(huán)境中運(yùn)行時(shí)，會(huì)使電池組長(zhǎng)期處于高倍率放電狀態(tài)，應(yīng)考慮加強(qiáng)對(duì)電池組的冷卻措施.

表3 不同放電倍率下電池組的溫度特性Tab.3 Temperature characteristics of battery at different discharge rates

5 結(jié)論

通過對(duì)單體鋰離子電池的熱特性和電池包熱流場(chǎng)的分析與研究，得到以下3個(gè)主要結(jié)論.

1) 成組電池單體電池間的最大溫差與放電倍率成正比.考慮到實(shí)際汽車運(yùn)行工況需要，尤其在爬坡、加速等大功率使用的工況下，應(yīng)設(shè)計(jì)電池?zé)峁芾硐到y(tǒng)以保證電池組的熱穩(wěn)定性，減小功率的損失和對(duì)電池造成的損害.

2) 成組電池的溫度場(chǎng)與單體電池所組成的溫度場(chǎng)沒有直接的疊加關(guān)系.在0.5C放電時(shí)，電池單體的最大溫差為13.30 ℃，而電池組則達(dá)到7.40 ℃.因此，電池組溫度場(chǎng)的均勻性對(duì)于電池組壽命和電池性能尤其重要.

3) 通過CFD的數(shù)值模擬仿真方法，可詳細(xì)了解電池箱箱體內(nèi)部的流場(chǎng)流向和溫度的總體分布情況.根據(jù)所得出的結(jié)果，對(duì)電池箱體的通風(fēng)冷卻結(jié)構(gòu)作出部分優(yōu)化，提出了在箱體進(jìn)口端增加導(dǎo)流板的改進(jìn)方案，有效地提高了溫度場(chǎng)分布的均勻性.結(jié)果表明，在0.5C放電倍率下的各電池單體最大溫差為7.40 ℃；而在1.0C放電倍率下的各電池單體最大溫差為8.10 ℃.

[1] SIEVERS M，SIEVERS U，MAO S S.Thermal modeling of new Li-ion cell design modifications[J].Forsch Ingenieurwes，2010，74(4):215-231.DOI:10.1007/s10010-010-0127-y.

[2] CHEN S C，WAN C C，WANG Y Y.Thermal analysis of lithium-ion batteries[J].Journal of Power Sources，2005，140(1):111-124.DOI:10.1016/j.jpowsour.2004.05.064.

[3] 辛乃龍.純電動(dòng)汽車鋰離子動(dòng)力電池組熱特性分析及仿真研究[D].長(zhǎng)春:吉林大學(xué)，2012:29-30.

[4] 李騰，林成濤，陳全世.鋰離子電池?zé)崮Ｐ脱芯窟M(jìn)展[J].電源技術(shù)，2009，33(10):927-932.

[5] 陳燕虹，吳偉靜，劉宏偉，等.純電動(dòng)汽車電池箱的熱特性[J].吉林大學(xué)學(xué)報(bào)(工學(xué)版)，2014，44(4):925-953.

[6] 羅玉濤，羅卜爾思，郎春艷.鋰離子動(dòng)力電池組的直接接觸液體冷卻方法研究[J].汽車工程，2016，38(7):909-914.

[7] 王福軍.計(jì)算流體動(dòng)力學(xué)分析CFD軟件原理與應(yīng)用[M].北京:清華大學(xué)出版社，2004：78-79.

[8] AMANO R S.Development of a turbulence near-wall model and its application to separated and reattached flows[J].Numerical Heat Transfer，2004，7(1):59-75.DOI:10.1080/10407798408546956.

[9] LAUNDER BE，SAMARAWEER D S A.Application of a second-moment turbulence closure to heat and mass transport in thin shear flows (Ⅰ): Two dimensional transport[J].International Journal of Heat and Mass Transfer，1999，22(12):1631-1643.DOI:10.1016/0017-9310(79)90080-2.

[10] 彭強(qiáng).電動(dòng)汽車用鋰離子動(dòng)力電池?zé)嵝?yīng)研究[D].長(zhǎng)春:吉林大學(xué)，2012：9-10.

[11] 車杜蘭.電動(dòng)汽車鋰離子電池包熱特性研究與優(yōu)化[D].武漢:武漢理工大學(xué)，2009：21-22.

[12] 張?jiān)圃?，白潔，張?guó)慶.大倍率放電時(shí)電動(dòng)汽車用鋰離子電池的熱性能[J].汽車安全與節(jié)能學(xué)報(bào)，2015，6(1):97-101.

[13] 馮竟陽，戴作強(qiáng)，張紀(jì)鵬，等.基于Ansys Workbench 12.0的磷酸鐵鋰動(dòng)力電池溫度場(chǎng)特性研究[J].青島大學(xué)學(xué)報(bào)(工程技術(shù)版)，2013，28(1):51-55.

(責(zé)任編輯： 錢筠英文審校： 崔長(zhǎng)彩)

NumericalSimulationandOptimizationDesignofThermalFlowFieldCharacteristicsofLithium-IonBatteryPackforVehicles

HAN Fenggang1，2， PAN huaijie1， PENG Qian1，2， LIN Youhuai2， LIU Sheng1

(1. School of Mechanical and Automotive Engineering， Xiamen University of Technology， Xiamen 361024， China；2. Xiamen Vehicle Design & Services Company Limited， Xiamen 361024， China)

10.11830/ISSN.1000-5013.201704076

2017-04-24

韓鋒鋼(1965-)，男，副教授，博士，主要從事新能源汽車研發(fā)、客車整車及零部件CAE分析及輕量化的研究.E-mail:hanfg@xmut.edu.cn.

國(guó)家自然科學(xué)基金資助項(xiàng)目(51505403)； 國(guó)家高端外國(guó)專家項(xiàng)目(GDT20153600065)； 福建省教育廳客車及特種車輛研發(fā)協(xié)同創(chuàng)新中心2011計(jì)劃項(xiàng)目(2016AYF004)； 福建省高校杰出青年科研人才培育計(jì)劃項(xiàng)目(2016年度)； 廈門理工學(xué)院研究生科技創(chuàng)新計(jì)劃項(xiàng)目(YKJCX2016001)

U 469.72

A

1000-5013(2017)05-0625-07