張 帥 , 曾銳利
(1.軍事交通學院 研究生管理大隊,天津 300161;2. 軍事交通學院 軍用車輛系,天津 300161)
基于振動信號的柴油發(fā)動機轉速測量
張 帥1, 曾銳利2
(1.軍事交通學院 研究生管理大隊,天津 300161;2. 軍事交通學院 軍用車輛系,天津 300161)
利用振動信號測量發(fā)動機轉速是一種不解體檢測的重要方法,但其對發(fā)動機振動頻率要求較高,因此轉速急劇變化時,對測量精度影響較大。提出了基于希爾波特包絡和自相關的新方法;利用小波變換消除基線漂移,再通過短時時間窗進行類穩(wěn)態(tài)分析,并采用希爾伯特包絡和三次自相關對信號進行時域分析,計算短時窗內的發(fā)動機轉速,通過移窗法得到連續(xù)的發(fā)動機轉速。通過在發(fā)動機臺上進行變工況試驗,同步采集振動信號和轉速傳感器信號進行對比分析,在急加速條件下轉速測量最大誤差小于1.4%,能夠有效的滿足柴油發(fā)動機轉速的實時性和精確性要求,在車輛不解體檢測和無負荷測功方面有很大應用前景。
柴油發(fā)動機轉速;振動信號;基線漂移;短時窗;希爾伯特變換;自相關
發(fā)動機轉速是發(fā)動機的一個重要參數(shù),是發(fā)動機運行狀況的的綜合體現(xiàn)[1]。非穩(wěn)態(tài)轉速的測量對于評價發(fā)動機性能有著極為重要的作用,現(xiàn)在的絕大多數(shù)測量手段(光電法和霍爾法)都需要在發(fā)動機機體上安裝傳感器,安裝復雜,費用高,不滿足柴油車輛的不解體檢測要求,而通過振動信號進行發(fā)動機非穩(wěn)態(tài)轉速的測量可以極大地降低車輛發(fā)動機性能的檢測要求門檻,能夠滿足便捷性和精度要求?,F(xiàn)有的一些基于振動信號測量轉速的技術手段多為汽油發(fā)動機設計,而柴油發(fā)動機較汽油發(fā)動機來說,相同運轉條件下,振動信號更為復雜,目前的方法不能很好的解決這些問題,另外,現(xiàn)有的關于柴油發(fā)動機轉速測量的方式也存在局限性,在非穩(wěn)態(tài)轉速條件下,尤其在發(fā)動機急加速條件下,頻率混疊比較嚴重,直接進行頻域分析得到的結果不能作為計算發(fā)動機振動基頻的依據[2-5];另外,引起發(fā)動機振動的因素很多,在對算法進行驗證時,將實車振動信號簡單等同于正弦信號與噪聲信號的疊加顯然是不夠準確的;文獻[6]進行誤差分析時,利用OBD(On Board Diagnostics)接口采集的轉速作為真實轉速進行驗證,這樣的轉速是經過發(fā)動機ECU(Electronic Control Unit)處理后發(fā)送出來的,傳輸上有一定的延時,對比驗證時存在誤差。
本文提出了一種基于振動信號的發(fā)動機轉速測量方法,但不同于已經發(fā)表的其他轉速測量的文章,本文采用的方法利用小波消除基線漂移,通過截取短時間的振動信號,利用希爾伯特包絡提取振動信號中缸體的振動規(guī)律,再利用三次自相關找到發(fā)動機每缸間的時間差,根據發(fā)動機的缸數(shù)計算極短時間內的轉速,最后通過移窗得到轉速-時間曲線。
這種轉速測量的新方法,振動傳感器可安裝在發(fā)動機的任意缸體上,簡單高效,適用于柴油發(fā)動機實時轉速的測量。
1.1 振動信號的希爾伯特包絡
在信號處理領域,利用信號的包絡對信號進行提取是一種常用的信號處理方法。常用的信號包絡提取方法有:極值法、檢波濾波法、Hilbert解調法和復解析小波變換等[7-9],它們各自的特點如表1所示。
表1 四種包絡求解方法比較
復解析小波變換需要針對信號的特點進行小波類型的選擇,并根據感興趣頻段進行時間平移和尺度伸縮因子的選擇,對不同的振動信號的適應性不高;而Hilbert解調法操作簡單,精度高,但保留了信號中的所有頻率,故本文采用先進行低通濾波再進行Hilbert包絡,這樣就能有效的去除信號中的高頻干擾。
Hilbert變換在本質上是一種全通濾波器,Hilbert變換巧妙地應用解析表達式中的實部與虛部的正弦和余弦關系,定義出任意時刻的瞬時頻率、瞬時相位及瞬時幅度,使得對于短信號和復雜信號的瞬時參數(shù)的提取成為可能,從而能更有效地、真實地獲取信號中所含的信息,有利于分析短時間的柴油發(fā)動機振動信號。利用Hilbert變換將實信號變換成解析信號,其模即為包絡信號。Hilbert變換提取實信號包絡的方法如下
實信號s(t)的復解析表達式如式(1)所示
(1)
(2)
1.2 三次自相關分析
經過Hilbert包絡計算,在短時間內包絡波形的頻率與原始振動信號中各缸工作振動波形的頻率一致,因此采用自相關方法分析包絡信號中的周期信號可以用來分析發(fā)動機振動信號中的周期。
對同一隨機振動樣本函數(shù)隨時間坐標移動進行相似程度計算稱為自相關分析,得到的函數(shù)稱為自相關函數(shù)。自相關函數(shù)就是信號x(t)和它自己的相關性,表達式如式(3)所示。
Rxx(t,τ)=E[x(t)x(t+τ)]=
(3)
將測得的振動信號表示為:x(t)=s(t)+n(t);s(t)為有用信號;n(t)為噪聲信號[10]。
則振動信號的自相關函數(shù)如式(4)所示。
Rxx(t,τ)=E[x(t)x(t+τ)]=
E[(s(t)+n(t))(s(t+τ)+n(t+τ))]=
Rss(t,τ)+Rsn(t,τ)+Rns(t,τ)+Rnn(t,τ)
(4)
由于噪聲信號與其他信號沒有相關性,故
Rsn(t,τ)=0,Rns(t,τ)=0,Rxx(t,τ)=Rss(t,τ)+Rnn(t,τ)。
再對振動信號的自相關函數(shù)求自相關,如式(5)所示。
(Rss(t,τ1+τ)+Rnn(t,τ1+τ)))=
Rss(Rss(t,τ1))+Rnn(Rnn(t,τ1))=
(5)
進而可得到式(6)
(6)
發(fā)動機在平穩(wěn)轉速下,其振動信號具有明顯的周期性。通過自相關分析,可以有效的濾除信號內噪聲信號,提高了振動信號的信噪比。構造頻率為10 Hz的信號,其中加入高斯噪聲,試驗表明,自相關的次數(shù)越多,其降噪效果越好,但同時計算量增加,當進行三次自相關運算之后,噪聲已經被極大程度的削弱如圖1所示。
圖1 不同次數(shù)自相關頻譜對比Fig.1 Spectral contrast of autocorrelation with different times
進行柴油發(fā)動機的轉速測量需要經過小波消除基線漂移、選取起始點和分析點數(shù)、濾波、Hilbert包絡、三次自相關分析等過程,其算法流程如圖2所示。
2.1 基線漂移的濾除
在發(fā)動機振動信號的實時處理中,往往有很多的噪聲干擾,而對信號影響最大的就是基線漂移噪聲。在發(fā)動機處于穩(wěn)態(tài)時,其振動圖像對稱性較好;當發(fā)動機處于非穩(wěn)態(tài),尤其在急加速狀態(tài)下,振動信號容易出現(xiàn)漂移,對后續(xù)的處理造成較大干擾。而小波變換是一種窗口大小不變但形狀可變,即時間窗和頻率窗都可以改變的時頻局部化分析方法,這一特點使得小波分析方法具有對信號的自適應性。
圖2 柴油發(fā)動機轉速測量算法流程Fig.2 Diesel engine speed measurement algorithms process
本文是利用coif3小波多尺度分解后所得到的低頻逼近信號充分逼近發(fā)動機振動信號中的基線漂移噪聲的特性,對其進行均一化處理。具體過程是:選擇適當?shù)男〔ê瘮?shù)和分解水平n,將原始的振動信號進行分解,得到逼近信號為cAn,使得cAn盡可能逼近振動信號中的基線漂移噪聲[11],將該逼近信號置零,以此來消除逼近信號所包含的基線漂移的信息,再將其與分解時所得的細節(jié)信號進行重構,則得到消除基線漂移干擾后的振動信號。
通過小波分解的逼近信號來消除振動信號中的基線漂移,效果如圖3所示。原有的基線漂移在很大程度上被消除了。
(a) 出現(xiàn)零點漂移的信號
(b) 經多尺度分解后的漂移噪聲
(c) 消除基線漂移后的信號
2.2 時間窗長度的選取
(7)
(8)
式中,n為柴油發(fā)動機的缸數(shù)。
由于振動信號要進行三次自相關分析,故要在保證滿足采樣振動信號的時間足夠短和時間窗的長度足夠進行三次自相關分析的要求,如式(9)所示
N0(N0(n,3 000),N0(n,600))
(9)
根據三次自相關后要保留至少三個主峰,故時間窗內的采樣點數(shù)如式(10)所示
N≥5×N0(n,600)
(10)
2.3 信號濾波及基于Hilbert變換的包絡分析
在發(fā)動機轉速測量中,我們只對信號中的低頻信號感興趣,為了抑制高頻噪聲對信號的干擾,采用低通濾波對時間窗內信號進行濾波。對經濾波處理后的信號進行基于Hilbert變換的包絡分析,如圖4所示。經過包絡分析之后,原信號的主要尖峰特征都很好的凸顯出來,極大的提高了“興趣”數(shù)據在信號中的比例。
(a) 原始振動信號 (b) 低通濾波后信號 (c) 包絡后信號圖4 信號經濾波、包絡處理后效果對比Fig.4 Comparison of signal filtering and envelope
2.4 三次自相關
對信號求取三次自相關,此時得到的曲線的周期性很明顯,如圖5所示。取前三尖峰所對應的偏移量分別為l1、l2、l3,則周期τ(k)和速度v(k)分別為
(11)
式中,k為當前求解的速度點的序號。
(a) 一次自相關 (b) 二次自相關 (c) 三次自相關 (d) 四次自相關圖5 不同次數(shù)自相關函數(shù)Fig.5 Different degree autocorrelation function
2.5 移窗方法和速度-時間曲線生成
進行完一個時間窗的分析后,需要對時間窗進行移動以測取不同時間點的轉速值,其移窗值如式(12)
(12)
速度-時間曲線由多個關鍵速度點vpk(tk,vk)構成,由于假定短時時間窗內轉速穩(wěn)定,所以以時間窗的中點作為速度點,取值如式(13)
(13)
式中,k為當前求解的速度點的序號。
而柴油發(fā)動機在非穩(wěn)態(tài)轉速下,振動信號的周期特征不明顯,可以截取極短時間內的振動信號,認為該極短時間內的轉速是穩(wěn)定的。在此基礎上,對振動信號進行三次自相關分析。
臺架試驗的發(fā)動機采用了濰柴WD615型六缸四沖程柴油發(fā)動機,振動傳感器采用的是IMI-603C01型加速度傳感器,轉速傳感器采用霍爾式轉速傳感器。振動傳感器吸附在發(fā)動機任意缸位置,轉速傳感器、振動傳感器及采集系統(tǒng)的安裝情況見圖6~圖8。利用NI數(shù)據采集系統(tǒng)采集振動傳感器和轉速傳感器的信號,采樣頻率為60 000 Hz,根據式(1),選取每次采樣點數(shù)為10 000點。
圖6 安裝轉速傳感器Fig.6 Install the speed sensor
圖7 安裝振動傳感器Fig.7 Install the vibration sensor
圖8 信號采集系統(tǒng)Fig.8 Signal acquisition system
臺架驗證時,均在空載條件下進行。首先對平穩(wěn)轉速信號進行驗證,如表2所示。
然后進行發(fā)動機非平穩(wěn)狀態(tài)下的驗證,試驗中采用的非平穩(wěn)狀態(tài)為手動調節(jié)油門實現(xiàn)的發(fā)動機加速狀態(tài)。
表2穩(wěn)態(tài)下的轉速測量誤差值
Tab.2Errorvalueofrotationalspeedmeasurementundersteadystate
真實轉速/(r·min-1)計算轉速值/(r·min-1)誤差/%8008040.5120012050.42140014070.5160016040.25
經處理后,柴油發(fā)動機轉速還原情況如圖9所示。
從圖9可知,本算法計算的轉速曲線與真實轉速曲線基本吻合,能夠滿足發(fā)動機非穩(wěn)態(tài)條件下的轉速測量,通過對多次測試的結果進行分析,得到本算法在非穩(wěn)態(tài)條件下的最大誤差為1.4%,符合轉速測量的精度要求。
(a)急加速條件下轉速恢復 (b)緩加速條件下轉速恢復 (c)急加速條件誤差曲線 (d)緩加速條件誤差曲線圖9 加速條件下轉速恢復及誤差分析Fig.9 Speed recovery and error analysis under accelerated conditions
利用小波消除了基線漂移,采用Hilbert包絡和三次自相關能夠很好的提取出振動信號中的轉速信息,測量精度高,時間分辨率高,能夠滿足穩(wěn)態(tài)和非穩(wěn)態(tài)條件下的轉速測量,適應性好,具體如下:
(1)本算法能夠滿足車輛不解體條件下的轉速測量,硬件安裝方便,可將振動傳感器吸附在發(fā)動機的任一缸體上,對車輛適應性好。
(2)本算法采用小波變換、希爾伯特包絡和三次自相關分析,算法結構簡單易行,執(zhí)行效率高,精度誤差小于1.4%,可滿足穩(wěn)態(tài)和非穩(wěn)態(tài)轉速條件下的實時轉速測量要求。
(3)本算法的時間分辨率高,發(fā)動機每個工作循環(huán)可進行n(發(fā)動機缸數(shù))次速度點的計算??梢跃珳诗@取發(fā)動機狀態(tài),在車輛不解體檢測和無負荷測功中有很大的應用前景。
[ 1 ] 李華英,果霖,付宏財. 基于單片機的發(fā)動機轉速測量系統(tǒng)設計[J]. 云南農業(yè)大學學報(自然科學版),2015, 30(2): 294-297.
LI Huaying, GUO Lin,FU Hongcai. Engine rotational speed measuring system design based on single chip computer[J]. Journal of Yunnan Agricultural University (Natural Science),2015,30(2):294-297.
[ 2 ] 丁柏,宋翔. 基于振動基頻的一種發(fā)動機轉速測量新方法[J]. 長江大學學報(自然科學版),2013, 10(10): 76-79.
DING Bai, SONG Xiang. A new engine speed measurement method based on vibration fundamental frequency[J]. Journal of Yangtze University(Natural Science),2013, 10(10): 76-79.
[ 3 ] LIN Hinbin, DING Kang. A new method for measuring engine rotational speed based on the vibration and discrete spectrum correction technique[J]. Measurement, 2013,46(7):2056-2064.
[ 4 ] SUN Yiquan, ZHANG Yingtang, LI Zhining, et al. Design on instantaneous rotational speed measuring device of engine[C]∥ International Conference on Instrumentation. Beijing:IEEE, 2011:126-129.
[ 5 ] 蔣強,孫利生. 基于振動信號和小波變換的轉速測量方法研究[J]. 汽車科技,2011,38(1):39-41.
JIANG Qiang,SUN Lisheng. Research on the speed measurement based on the vibration signal and wavelet Trousform[J]. Auto Mobile Science amp; Technology, 2011,38(1):39-41.
[ 6 ] 林慧斌,丁康. 基于振動的發(fā)動機轉速高精度動態(tài)測量新方法[J]. 振動工程學報,2012,25(5):595-601.
LIN Huibin, DING Kang. A high accuracy dynamical measurement method for engine rotation speed based on vibration[J]. Journal of Vibration Engineering, 2012,25(5):595-601.
[ 7 ] 張永生.基于Hilbert包絡信號提取技術的LCR波超聲應力檢測系統(tǒng)研究[D].成都:電子科技大學,2013.
[ 8 ] 張盈盈,潘宏俠,鄭茂遠.基于小波包和Hilbert包絡分析的滾動軸承故障診斷方法[J].電子測試, 2010,17(6): 20-23.
ZHANG Yingying, PAN Hongxia, ZHENG Maoyuan. Rolling bearings fault diagnosis based on wavelet packet and Hilbert envelope analysis[J]. Electronic Test, 2010,17(6): 20-23.
[ 9 ] ZHOU Jiadong, OUYANG Guangyao, GAO Hongbin. Demodulation of instantaneous rotational speed of diesel engine based on hilbert transform[C]∥ Proceedings of 2009 9th International Conference on Electronic Measurement amp; Instruments.Beijing:IEEE, 2010.
[10] 衰文文. 基于振動分析的汽車發(fā)動機在線轉速測量研究[D].天津:河北工業(yè)大學,2011.
[11] 張文瓊,劉肖琳,吳濤. 一種利用小波變換逼近信號濾除心電圖基線漂移的方法[J]. 計算機工程與應用, 2005,41(20):222-224.
ZHANG Wenqiong, LIU Xiaolin, WU Tao. A method using wavelet approximation to remove the electrocardiogram baseline wander[J]. Computer Engineering and Applications, 2005,41(20):222-224.
《振動與沖擊》征稿簡則
《振動與沖擊》是中國振動工程學會、上海交通大學、上海市振動工程學會聯(lián)合主辦的反映和交流我國力學、航空、航天、造船、車輛、機械、土木建筑、水利、電力、礦山、化工、冶金、儀器儀表、海洋工程、輕工、環(huán)境保護等領域的振動、沖擊和噪聲方面的成果及經驗的綜合性學術刊物。本刊為國家科委中國科技論文統(tǒng)計源期刊、《中國核心期刊要目總覽》核心期刊,《中國知識資源總庫中國科技期刊精品數(shù)據庫》、《中國科學引文數(shù)據庫》、《中國物理學文獻數(shù)據庫》、《中國學術期刊綜合評價數(shù)據庫》、《中國學術期刊文摘》源期刊,美國工程信息公司Ei-Compendex數(shù)據庫收錄期刊,美國《劍橋科學文摘》等的來源期刊。
1 刊登內容
本刊刊登上述領域的結構動力分析、模態(tài)分析、參數(shù)分析、參數(shù)識別、隨機振動、非線性振動、振動控制、轉子動力學、結構動力穩(wěn)定性、流固耦合振動、減振、隔振、抗沖擊、噪聲防治、環(huán)境試驗、模擬技術、測試技術、信號處理、計算機軟件工程、消聲器材等方面的論文和短文。
本刊將優(yōu)先發(fā)表有創(chuàng)新性工作、工程應用背景的論文;刊登在本學科中能解決工程問題而具有較大經濟效益的應用性文章,同時也將刊登一些介紹國內外新技術、新設備的文章和國內外學術動態(tài)。
2 投稿須知
2.1 稿件內容必須具有創(chuàng)新性、學術性、科學性和準確性、規(guī)范性和可讀性。
2.2 請勿一稿兩投或多投。文稿自收到之日起,3個月內本刊將發(fā)出是否錄用通知;逾期請及時向編輯部查詢。稿件不退,請自留底稿。
2.3 編輯部將按照規(guī)范的程序,聘請有關同行專家評審(二審制,部分稿件須三審)和主編終審,根據評審意見公平、公正地決定稿件的取舍。
2.4 稿件文責自負。編輯部對來稿有權作技術性和文字性的修改,實質性內容修改須征得作者同意。
2.5 稿件的電子文檔通過網址:http:∥jvs.sjtu.edu.cn或E-mail: jvs@sjtu.edu.cn發(fā)到編輯部。
2.6 稿件內容不得涉及國家機密,如涉及的請?zhí)峁┍C茏C明材料。
3 稿件要求和注意事項
3.1 來稿要求論點明確、數(shù)據可靠、邏輯嚴密、文字簡煉。來稿必須包括題名、作者姓名、單位及所在城市、郵編、中英文摘要和關鍵詞、中圖分類號[請查閱中國圖書分類號(第4版)]、第一和通訊作者簡介(姓名、性別、學位、職稱、出生年月、聯(lián)系電話、E-mail),英文摘要放在中圖分類號和文獻標識碼下面。請在稿件首頁地腳處注明論文屬何基金項目資助及項目編號。
3.2 論文摘要的內容獨立于正文而存在,它能準確、具體、完整地概括原文的創(chuàng)新之處。摘要內容在200字~300字左右,英文摘要須與中文摘要相對應,中英文摘要一律采用第三人稱表述,不使用“本文”、“作者”等作為主語。
3.3 文中圖、表應有自明性,且隨文出現(xiàn)。插圖須符合制圖規(guī)范。圖中文字、縱橫坐標中的標值、標值線(向內)必須寫清,標目應使用法定計量單位。圖名、表名用中英文表示。文中表格一律使用“三線表”,必要時可加輔助線,表的內容切忌與圖和文字內容重復。
3.4 文稿章節(jié)編號采用三級標題頂格排序。一級標題形如1,2,3…排序; 二級標題形如1.1,1.2…;2.1,2.2…排序;三級標題形如1.1.1,1.1.2…;2.1.1,2.1.2…排序;引言不排序。
3.5 參考文獻應是文中直接引用的公開出版物,以6篇以上為宜,其中80%應為期刊或會議論文,80%以上為近5年出版的文獻,40%以上為外文文獻。須按在文中出現(xiàn)的順序編號,列于文后。文獻作者3名以內全部列出, 3名以上則列前3名,后加“,等”(外文加“,et al”)。外文作者采用姓前名后格式,名用縮寫,不加縮寫點。參考文獻的中文期刊須使用中、英雙語著錄。各類文獻的著錄格式如下:①專著:[序號]著者.書名[M].版本(第1版免著).出版地:出版者,出版年.②期刊:[序號]作者.題名[J].刊名,出版年,卷(期):起止頁碼.③論文集、會議錄:[序號]作者.題名[C]//論文集名.出版地:出版者,出版年.④學位論文:[序號]作者.題名[D].保存地點:保存單位,年份.⑤科技報告:[序號]主要責任者.題名[R].出版地:出版單位,出版年.⑥專利:[序號]專利所有者.專利題名[P].專利國別:專利號,出版日期.⑦標準:[序號]標準編號,標準名稱[S].⑧電子文獻:[序號]主要責任者.題名[文獻類型標識/載體類型標識].電子文獻的出處或可獲得的地址,發(fā)表或更新日期/引用日期.其中:文獻類型標識為M、C、D、R、J、S、P、DB(數(shù)據庫)、CP (計算機程序)、EB (電子公告);載體類型標識為MT (磁帶)、DK (磁盤)、CD (光盤)、OL (聯(lián)機網絡)。
4 網上檢索
本刊網址:http:∥jvs.sjtu.edu.cn
中國期刊網:http:/www.chinajournal.net.cn 或http://www.cnki.net
萬方數(shù)據系統(tǒng)科技期刊群:http://202.120.13.32:83
《美國工程索引》(Ei):http://www.engineeringvillage2.org.cn
本刊已許可萬方數(shù)據、重慶維普、中國學術期刊(光盤版)電子雜志社在中國知網及其系列數(shù)據庫產品中以數(shù)字化方式復制、匯編、發(fā)行、信息網絡傳播本刊全文。作者向本刊提交文章發(fā)表的行為即視為同意我社上述聲明。
Speedmeasurementofdieselenginesbasedonvibrationsignals
ZHANG Shuai1, ZENG Ruili2
(1. Postgraduate Training Brigade,Military Transportation University,Tianjin 300161, China;2. Military Vehicle Department,Military Transportation University,Tianjin 300161, China)
It is an important method to measure the engine speed by using vibration signals. Engine vibration frequency has a great influence on the measurement accuracy when the speed changes dramatically. A new method based on Hilbert envelope and autocorrelation was proposed. Wavelet transform was used to eliminate the baseline drift firstly. Then the quasi-steady-state analysis was carried out by means of a short-time window, and the Hilbert envelope and three autocorrelations were used in time domain analysis to calculate the rotational speed of the engine in a short-time window. Finally, the continuous engine speed was obtained by the moving window method. Through a variable working condition experiment on the engine bench, the synchronous acquisition vibration signal and the rotating speed sensor signal were compared and analyzed, and the maximum error of speed measurement under the condition of rapid acceleration was less than 1.4%. It can effectively meet the requirements of real-time and accuracy of diesel engine's speed, and has great application potential in the field of non-disassembly tests and unloaded power detection.
diesel engine speed; vibration signal; baseline drift; short-time window; Hilbert transform; autocorrelation
國防預研項目(40407030401)
2016-11-07 修改稿收到日期:2017-01-12
張帥 男,碩士生,1992年生
曾銳利 男,博士,副教授,1976年生
TK427
A
10.13465/j.cnki.jvs.2017.22.039