任金波,張 翔,施火結(jié)
(福建農(nóng)林大學(xué) 機(jī)電工程學(xué)院,福建 福州 350002)
新能源汽車電池散熱風(fēng)扇軸向振動(dòng)分析與改進(jìn)
任金波,張 翔,施火結(jié)
(福建農(nóng)林大學(xué) 機(jī)電工程學(xué)院,福建 福州 350002)
針對(duì)福州某公司新能源汽車電池散熱風(fēng)扇軸向振動(dòng)較大的問題,研究了散熱風(fēng)扇系統(tǒng)的減振技術(shù)。分析了影響風(fēng)扇軸向振動(dòng)的因素,采用INV3160型智能信號(hào)采集儀和數(shù)據(jù)采集分析(DASP)軟件對(duì)風(fēng)扇樣機(jī)進(jìn)行了振動(dòng)測(cè)試試驗(yàn),得出影響風(fēng)扇軸向振動(dòng)的最大因素是直流無刷電機(jī)的換相脈動(dòng)和齒槽脈動(dòng)。對(duì)風(fēng)扇減振結(jié)構(gòu)進(jìn)行了改進(jìn)設(shè)計(jì),對(duì)外轉(zhuǎn)子電機(jī)進(jìn)行齒面開槽,并通過試驗(yàn)對(duì)改進(jìn)前后的結(jié)果進(jìn)行對(duì)比。試驗(yàn)結(jié)果表明:改進(jìn)減振結(jié)構(gòu)后,在低轉(zhuǎn)速時(shí),直流無刷電機(jī)的換相脈動(dòng)和齒槽脈動(dòng)總的功率譜幅值可降低50%以上;齒面開槽后,軸向振動(dòng)的諧波幅值下降幅度達(dá)45.8%,可有效減小直流無刷電機(jī)的高頻轉(zhuǎn)矩脈動(dòng)對(duì)風(fēng)扇振動(dòng)的影響。
新能源汽車;電池散熱風(fēng)扇;振動(dòng);減振
蓄電池是新能源汽車中油電混合動(dòng)力汽車和純電動(dòng)汽車的核心部件,在大電流工作或充電的過程中,會(huì)因電池內(nèi)阻而產(chǎn)生大量的熱能,從而引起蓄電池溫度快速升高和電解液的氣化,導(dǎo)致蓄電池爆炸而產(chǎn)生危險(xiǎn)。因而,需要采用風(fēng)扇進(jìn)行強(qiáng)制散熱以抑制蓄電池在新能源汽車行駛或充電過程中升溫。減振降噪是風(fēng)扇在工程應(yīng)用中急需解決的關(guān)鍵問題之一,國(guó)內(nèi)外學(xué)者在這方面進(jìn)行了較多研究。文獻(xiàn)[1-4]采用改變扇葉的形狀和角度等方法來降低汽車風(fēng)扇噪聲。文獻(xiàn)[5]通過改變消聲器的結(jié)構(gòu),來提高消聲性能,從而達(dá)到降低汽車空調(diào)連接管噪聲的目的。文獻(xiàn)[6]提出通過重構(gòu)氣隙磁密波形來降低噪聲的方法。文獻(xiàn)[7]研究了抑制電動(dòng)汽車用永磁同步電機(jī)轉(zhuǎn)矩脈動(dòng)的方法。
以往對(duì)風(fēng)扇減振降噪的研究,多基于一般汽車的發(fā)動(dòng)機(jī)冷卻系統(tǒng),對(duì)新能源汽車電池散熱方面的研究,目前主要集中在電池的熱效應(yīng),對(duì)風(fēng)扇具體工程案例的研究還較少。本文以福州某公司專為奇瑞汽車在研車型A21ISG油電混合汽車所設(shè)計(jì)的電池散熱風(fēng)扇為例,提高了現(xiàn)有風(fēng)扇樣機(jī)產(chǎn)品的可靠性、舒適性和節(jié)能環(huán)保性能,通過理論分析和振動(dòng)測(cè)試得到了影響風(fēng)扇軸向振動(dòng)的最大因素,對(duì)風(fēng)扇減振結(jié)構(gòu)和電機(jī)進(jìn)行改進(jìn),并通過試驗(yàn)對(duì)比驗(yàn)證改進(jìn)效果。
以福州某公司專為奇瑞汽車在研車型A21ISG油電混合汽車所設(shè)計(jì)的電池散熱風(fēng)扇為例,新能源汽車電池散熱風(fēng)扇樣機(jī)如圖1所示。風(fēng)扇機(jī)械部分主要由整體式外轉(zhuǎn)子葉輪、電機(jī)、減振結(jié)構(gòu)和法蘭盤底座等部件組成。其中,風(fēng)輪與轉(zhuǎn)子不是兩個(gè)獨(dú)立的零件,而是將電機(jī)外轉(zhuǎn)子的磁軛簡(jiǎn)化為一個(gè)簡(jiǎn)單的套筒,然后以緊配合的方式壓入風(fēng)輪,共同形成整體式的電機(jī)外轉(zhuǎn)子風(fēng)輪。電池散熱風(fēng)扇外轉(zhuǎn)子結(jié)構(gòu)和車用空調(diào)直流無刷風(fēng)機(jī)的外轉(zhuǎn)子結(jié)構(gòu)一致[8]。
新能源汽車電池散熱風(fēng)扇的葉輪固裝于無刷電機(jī)驅(qū)動(dòng)的轉(zhuǎn)軸上,電池散熱風(fēng)扇的工作原理如圖2所示。傳感器將采集到的軸向輸出信號(hào)傳送給控制器,經(jīng)脈沖寬度調(diào)制(pulse width modulation,PWM)后,驅(qū)動(dòng)控制電機(jī)運(yùn)轉(zhuǎn)從而帶動(dòng)葉輪不斷旋轉(zhuǎn),達(dá)到給電池組送風(fēng)散熱的目的。
圖1 新能源汽車電池散熱風(fēng)扇樣機(jī)圖2 電池散熱風(fēng)扇的工作原理
新能源汽車電池散熱風(fēng)扇振動(dòng)的原因復(fù)雜,軸向振動(dòng)、徑向振動(dòng)、扭擺振動(dòng)等均為常見的振動(dòng)形式。根據(jù)風(fēng)扇樣機(jī)的結(jié)構(gòu)特點(diǎn)和工作原理,影響樣機(jī)軸向振動(dòng)的主要因素有整體式外轉(zhuǎn)子質(zhì)量不平衡和直流無刷電機(jī)的各種轉(zhuǎn)矩脈動(dòng)。
2.1 整體式外轉(zhuǎn)子質(zhì)量不平衡
電池散熱風(fēng)扇整體式外轉(zhuǎn)子葉輪有41片葉片,制造誤差導(dǎo)致的葉片分布不均勻,葉輪孔和旋轉(zhuǎn)軸在裝配過程中存在的裝配誤差,葉輪內(nèi)圈上緊密貼合的磁軛加工誤差及分布不均等,都會(huì)造成外轉(zhuǎn)子質(zhì)量偏心,從而引起風(fēng)扇振動(dòng)。葉輪等外轉(zhuǎn)子質(zhì)量不平衡引起的振動(dòng)頻率分量[9]為:
(1)
其中:f為外轉(zhuǎn)子質(zhì)量不平衡時(shí)的振動(dòng)頻率,Hz;n為葉輪工作時(shí)轉(zhuǎn)速,r/min。
2.2 直流無刷電機(jī)的轉(zhuǎn)矩脈動(dòng)
2.2.1 換相脈動(dòng)
直流無刷電機(jī)換相時(shí)電磁轉(zhuǎn)矩的表達(dá)式[10]為:
(2)
其中:Te為直流無刷電機(jī)換相脈動(dòng)時(shí)的電磁轉(zhuǎn)矩,N·mm;E為電樞相繞組電勢(shì)的穩(wěn)態(tài)值,V;ωm為機(jī)械角速度,rad/s;I為各相電流的穩(wěn)態(tài)值,A;V為電機(jī)端電壓,V;L為定子繞組等效自感,mH;M為定子繞組等效互感,mH;t為換相時(shí)間,ms。
因此,換相時(shí)電流會(huì)發(fā)生脈動(dòng),設(shè)電機(jī)相數(shù)為k,對(duì)極數(shù)為p,電機(jī)換相時(shí)的電流脈動(dòng)頻率為:
(3)
樣機(jī)的直流無刷電機(jī)參數(shù)k=3,p=2,根據(jù)式(3)可知,風(fēng)扇樣機(jī)中直流無刷電機(jī)的換相脈動(dòng)頻率為轉(zhuǎn)子轉(zhuǎn)動(dòng)時(shí)工作頻率的12倍。
2.2.2 齒槽脈動(dòng)
齒槽轉(zhuǎn)矩是由于齒槽電樞鐵芯和永磁體磁場(chǎng)相互作用形成的。根據(jù)機(jī)電能量轉(zhuǎn)換中的能量方法,電樞繞組開路計(jì)算時(shí),齒槽轉(zhuǎn)矩可由電機(jī)磁場(chǎng)能量相對(duì)轉(zhuǎn)子轉(zhuǎn)角的導(dǎo)數(shù)[11-15]求得:
(4)
(5)
B(θ,θ0,l)=Λ(θ)F(θ,θ0,l);
(6)
(7)
其中:T為直流無刷電機(jī)的齒槽脈動(dòng)轉(zhuǎn)矩,N·mm;Wd為電機(jī)磁場(chǎng)能量,W;θ0為轉(zhuǎn)子轉(zhuǎn)角,rad;θ為電樞周向坐標(biāo),rad;l為軸向坐標(biāo),mm;B為等效氣隙磁密,T;g為氣隙長(zhǎng)度,mm;Dd為電樞直徑,mm;ha為磁鋼厚度,mm;μ0為真空磁導(dǎo)率;F(θ,θ0,l)為磁鋼磁勢(shì),A。
將式(7)代入式(5)可得:
(8)
式(8)經(jīng)過傅里葉變換后得:
(9)
(10)
其中:p為極對(duì)數(shù);Q為風(fēng)扇中電機(jī)定子槽數(shù);Λt為第t次磁導(dǎo)諧波幅值,V;fv為第v次磁勢(shì)諧波幅值,V;v為磁勢(shì)諧波的次數(shù)。將式(9)和式(10)代入式(5)得:
(11)
由式(11)可知:當(dāng)磁勢(shì)諧波、磁導(dǎo)諧波具有相同次數(shù)時(shí),會(huì)導(dǎo)致直流無刷電機(jī)的齒槽脈動(dòng),且產(chǎn)生力矩,此時(shí):
(12)
其中:n、β為Q和2p的公倍數(shù),是齒槽轉(zhuǎn)矩的次數(shù)。
新能源汽車電池散熱風(fēng)扇樣機(jī)中電機(jī)齒槽數(shù)Q為6,極數(shù)2p為4,其最小公倍數(shù)為12,可知風(fēng)扇樣機(jī)中直流無刷電機(jī)齒槽脈動(dòng)頻率為轉(zhuǎn)子轉(zhuǎn)動(dòng)時(shí)工作頻率的12倍。
將新能源電池散熱風(fēng)扇置于振動(dòng)試驗(yàn)臺(tái)上,啟動(dòng)風(fēng)扇,從1 200 r/min的正常工作轉(zhuǎn)速開始加速,轉(zhuǎn)速每增加100 r/min測(cè)量一次風(fēng)扇軸向振動(dòng)信號(hào),直至達(dá)到最高轉(zhuǎn)速3 000 r/min時(shí)結(jié)束。測(cè)試信號(hào)經(jīng)DLF-4型放大器放大、INV3160型智能信號(hào)采集儀低通濾波采樣后,用數(shù)據(jù)采集分析(data acquisition and signal processing,DASP)軟件對(duì)振動(dòng)信號(hào)進(jìn)行分析處理,得到風(fēng)扇在各穩(wěn)定轉(zhuǎn)速下的功率譜。風(fēng)扇以1 600 r/min和 2 000 r/min的轉(zhuǎn)速穩(wěn)定運(yùn)行時(shí),測(cè)試的軸向功率譜分別如圖3和圖4所示。
圖3 軸向功率譜(1600r/min)圖4 軸向功率譜(2000r/min)
由圖3和圖4可知:當(dāng)風(fēng)扇以1 600 r/min的轉(zhuǎn)速穩(wěn)定運(yùn)行時(shí),在320 Hz頻率分量處,軸向的功率譜有較大的幅值;當(dāng)風(fēng)扇以2 000 r/min的轉(zhuǎn)速穩(wěn)定運(yùn)行時(shí),在400 Hz頻率分量處出現(xiàn)較大幅值,表明在320 Hz和400 Hz頻率分量處,風(fēng)扇的振動(dòng)較大。經(jīng)過計(jì)算可知:該頻率分量均為軸向振動(dòng)測(cè)試時(shí)轉(zhuǎn)子轉(zhuǎn)動(dòng)時(shí)工作頻率的12倍。上文在分析影響風(fēng)機(jī)軸向振動(dòng)的因素時(shí),同樣得出電機(jī)中換相脈動(dòng)和齒槽脈動(dòng)的頻率為轉(zhuǎn)子轉(zhuǎn)動(dòng)時(shí)工作頻率的12倍。因此,影響風(fēng)扇軸向振動(dòng)的主要原因確定為電機(jī)的換相脈動(dòng)和齒槽脈動(dòng)。
圖5 風(fēng)扇軸向振動(dòng)時(shí)的時(shí)域波形
為進(jìn)一步驗(yàn)證測(cè)試試驗(yàn)和結(jié)論的正確性,用示波器測(cè)試風(fēng)扇軸向振動(dòng)時(shí)的時(shí)域波形,如圖5所示。圖5中:實(shí)測(cè)波形T1通道顯示外轉(zhuǎn)子周期;T2通道顯示外轉(zhuǎn)子的軸向齒槽脈動(dòng)波形。從圖5中可以看出:實(shí)測(cè)的諧波次數(shù)為12,與理論值轉(zhuǎn)子轉(zhuǎn)動(dòng)時(shí)工作頻率的12倍一致,并與功率譜圖的分析相符。
根據(jù)新能源汽車電池散熱風(fēng)扇的結(jié)構(gòu),可將其減振部分視為等效阻尼和彈簧系統(tǒng),將整體式電機(jī)外轉(zhuǎn)子風(fēng)輪等效為總體質(zhì)量。當(dāng)風(fēng)扇運(yùn)行時(shí),受外界激振力的作用,可能出現(xiàn)軸向振動(dòng)、徑向振動(dòng)和扭擺振動(dòng),本文主要研究其軸向振動(dòng)。在激振力作用下,軸向振動(dòng)可簡(jiǎn)化為單自由度受迫振動(dòng),其振幅和相位[16]分別為:
(13)
(14)
其中:A為單自由度振動(dòng)的振幅,mm;F0為激振力,N;k為剛度系數(shù),N·mm-1;δ為阻尼系數(shù),N·s·mm-1;ω0為干擾力圓頻率,rad/s;ω為風(fēng)扇自振圓頻率,rad/s;φ為相位角,rad。當(dāng)ω=ω0時(shí),振動(dòng)加大,系統(tǒng)發(fā)生共振;當(dāng)ω<ω0/2時(shí),系統(tǒng)振動(dòng)將會(huì)大大減小。
新能源汽車電池散熱風(fēng)扇的驅(qū)動(dòng)采用直流無刷電機(jī),但直流無刷電機(jī)普遍存在轉(zhuǎn)矩脈動(dòng)問題,如換相脈動(dòng)和齒槽脈動(dòng)。定子齒面開槽是降低直流無刷電機(jī)齒槽脈動(dòng)較簡(jiǎn)單有效的方法。開槽時(shí),槽面沿定子齒均勻分布。當(dāng)電機(jī)齒面開1個(gè)槽時(shí),其基本齒槽轉(zhuǎn)矩次數(shù)為12次;電機(jī)齒面開2個(gè)槽時(shí),基本齒槽轉(zhuǎn)矩次數(shù)變?yōu)?6次,脈動(dòng)次數(shù)增加,幅值隨之下降??紤]新能源汽車電池散熱風(fēng)扇的結(jié)構(gòu),本文主要通過改進(jìn)風(fēng)扇的減振結(jié)構(gòu),從而改變頻率和齒面開槽來減小風(fēng)機(jī)振動(dòng)。
圖6 散熱風(fēng)扇改進(jìn)前減振墊結(jié)構(gòu)
4.1 改進(jìn)減振結(jié)構(gòu)
新能源汽車電池散熱風(fēng)扇改進(jìn)前減振墊呈階梯圓柱狀結(jié)構(gòu),如圖6所示。采用3個(gè)減振墊通過內(nèi)中心孔與電機(jī)線架相連接,然后用螺釘通過減振蓋與法蘭盤底座固聯(lián)在一起。
通過多次試驗(yàn),將風(fēng)扇現(xiàn)有的3個(gè)階梯圓柱狀減振墊改為夾心式圓柱狀,減振材料不變。3個(gè)減振墊通過內(nèi)中心孔與電機(jī)線架相連接,上端外圓面與法蘭盤底座上的3個(gè)圓柱孔配合。改進(jìn)后的減振墊結(jié)構(gòu)及安裝位置如圖7所示。4.2 試驗(yàn)對(duì)比
改進(jìn)減振結(jié)構(gòu)后,再次將風(fēng)扇置于測(cè)試臺(tái),測(cè)試風(fēng)扇中直流無刷電機(jī)換相脈動(dòng)和齒槽脈動(dòng)在減振前后的功率譜,軸向振動(dòng)信號(hào)經(jīng)采樣處理分析后得到換相脈動(dòng)和齒槽脈動(dòng)總的功率譜,如圖8所示。從圖8中可以看出:當(dāng)風(fēng)扇改進(jìn)減振結(jié)構(gòu)后,在低轉(zhuǎn)速時(shí)直流無刷電機(jī)的換相脈動(dòng)和齒槽脈動(dòng)總的功率譜幅值有極大降低,降幅可達(dá)50%以上。在2 300 r/min以上高轉(zhuǎn)速狀態(tài)下,換相脈動(dòng)和齒槽脈動(dòng)總的功率譜基本呈無波動(dòng)近直線狀態(tài),且幅值和零點(diǎn)較接近,說明改進(jìn)減振結(jié)構(gòu)后風(fēng)扇的軸向振動(dòng)進(jìn)一步降低。
圖7 散熱風(fēng)扇改進(jìn)后的減振墊結(jié)構(gòu)及安裝位置圖8 改進(jìn)前后的換相脈動(dòng)和齒槽脈動(dòng)總的功率譜
4.3 齒面開槽
新能源汽車電池散熱風(fēng)扇中,電機(jī)為外轉(zhuǎn)子結(jié)構(gòu),對(duì)外轉(zhuǎn)子電機(jī)齒面開槽,如圖9所示。當(dāng)電機(jī)外齒面開2個(gè)矩形凹槽(槽寬2 mm,槽深1 mm)后,測(cè)試得到風(fēng)扇軸向振動(dòng)的時(shí)域波形,如圖10所示。由圖5和圖10可以看出:電機(jī)齒面開槽后,齒槽脈動(dòng)周期增加,諧波次數(shù)倍增,周期內(nèi)變?yōu)?6次。諧波幅值從 830 mV下降到450 mV,較齒面開槽前下降幅度達(dá)45.8%,表明齒面開槽對(duì)抑制直流無刷電機(jī)的齒槽脈動(dòng)有較好效果。
圖9 外轉(zhuǎn)子電機(jī)齒面開槽圖10 電機(jī)齒面開槽后的風(fēng)扇軸向振動(dòng)時(shí)域波形
(1)根據(jù)新能源汽車電池散熱風(fēng)扇的結(jié)構(gòu)特點(diǎn),通過理論分析和振動(dòng)測(cè)試得出,影響新能源汽車電池散熱風(fēng)扇軸向振動(dòng)的主要因素為風(fēng)扇中直流無刷電機(jī)的換相脈動(dòng)和齒槽脈動(dòng)。
(2)改進(jìn)風(fēng)扇的減振結(jié)構(gòu)后,低轉(zhuǎn)速時(shí)直流無刷電機(jī)的換相脈動(dòng)和齒槽脈動(dòng)總的功率譜幅值降低達(dá)50%以上,高轉(zhuǎn)速狀態(tài)下幅值可進(jìn)一步降低。
(3)外轉(zhuǎn)子電機(jī)齒面開槽后,軸向振動(dòng)的諧波幅值較開槽前下降幅度達(dá)45.8%,表明齒面開槽對(duì)抑制直流無刷電機(jī)的齒槽脈動(dòng)有較好效果。
[1] SINGH O P,KHILWANI R,SREENIVASULU T,et al.Parametric study of centrifugal fan performance:experiments and numerical simulation[J].International journal of advances in engineering & technology,2011,1(2):33-50.
[2] 王永生,楊愛玲,戴韌,等.低壓頭小流量軸流風(fēng)扇氣動(dòng)噪聲特性研究[J].工程熱物理學(xué)報(bào),2010,31(7):1131-1134.
[3] SCHEIT C,KARIC B,BECKER S.Effect of blade wrap angle on efficiency and noise of small radial fan impellers-a computational and experimental study[J].Journal of sound and vibration,2012,331(5):996-1010.
[4]SEUNGH,CHEOLUNGC,TAE-HOONK.DevelopmentoflownoisecentrifugalfansforarefregeratorusinginclinedS-shapedtrailingedge[J].Internationaljournalofrefrigeration,2011,34(8):2076-2091.
[5] 費(fèi)標(biāo)求.汽車空調(diào)連接管聲學(xué)及動(dòng)剛度特性的研究[D].合肥:合肥工業(yè)大學(xué),2015.
[6] 于慎波,姜菲菲,王輝,等.永磁同步電主軸分?jǐn)?shù)槽電機(jī)的徑向電磁力分析[J].組合機(jī)床與自動(dòng)化加工技術(shù),2014(6):15-18,24.
[7] 王碩,康勁松,鐘再敏,等.電動(dòng)汽車用永磁同步電機(jī)轉(zhuǎn)矩脈動(dòng)抑制方法綜述[J].電源學(xué)報(bào),2016,14(5):24-32.
[8] 張翔,胡冰樂,張祖峰.車用空調(diào)無刷風(fēng)機(jī)整體式外轉(zhuǎn)子風(fēng)輪:2006200099747[P].(2007-09-12)[2016-10-20].
[9] 陳虹微.旋轉(zhuǎn)機(jī)械振動(dòng)特征及診斷方法[J].噪聲與振動(dòng)控制,2009,29(1):134-136.
[10] 徐衡,周云山.電動(dòng)車無刷直流電機(jī)轉(zhuǎn)矩脈動(dòng)抑制[J].計(jì)算機(jī)仿真,2016,33(4):195-199.
[11] 于慎波,駱開軍.氣隙形狀對(duì)永磁同步電主軸齒槽轉(zhuǎn)矩的影響[J].機(jī)床與液壓,2016,44(13):10-13.
[12] 古海江,黃文美,王超,等.定子齒開槽對(duì)內(nèi)置式永磁電機(jī)齒槽轉(zhuǎn)矩的影響[J].電機(jī)與控制應(yīng)用,2016,43(8):40-45.
[13] 黃燕濤,郭新華,項(xiàng)雷軍.內(nèi)置式永磁電機(jī)齒槽轉(zhuǎn)矩的優(yōu)化設(shè)計(jì)[J].華僑大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,37(5):536-540.
[14] 王建設(shè),鄒海榮.利用轉(zhuǎn)子輔助孔削弱永磁電動(dòng)機(jī)齒槽轉(zhuǎn)矩的方法[J].上海電機(jī)學(xué)院學(xué)報(bào),2016,19(3):141-146.
[15] 尚海,竇滿峰,陳敏,等.永磁同步電機(jī)極槽組合優(yōu)化的研究[J].科學(xué)技術(shù)與工程,2014,14(11):27-30.
[16] 莫文輝.單自由度強(qiáng)迫振動(dòng)的可靠性[J].湖北工業(yè)職業(yè)技術(shù)學(xué)院學(xué)報(bào),2013,26(1):97-99.
福建省科技廳高校產(chǎn)學(xué)研重大專項(xiàng)基金項(xiàng)目(2011H6004);福建省教育廳專業(yè)改革基金項(xiàng)目(1114S1205);福建農(nóng)林大學(xué)機(jī)械工程學(xué)科整體學(xué)科水平提升計(jì)劃基金項(xiàng)目(612014049)
任金波(1982-),男,湖北隨州人,實(shí)驗(yàn)師,博士生,主要研究方向?yàn)闄C(jī)械振動(dòng)和優(yōu)化設(shè)計(jì);張翔(1957-),男,通信作者,福建福州人,教授,碩士生導(dǎo)師,主要研究方向?yàn)楝F(xiàn)代設(shè)計(jì)方法和優(yōu)化.
2016-11-21
1672-6871(2017)04-0025-05
10.15926/j.cnki.issn1672-6871.2017.04.006
TH113.1
A