王昊利, 汪 兵
(1. 金陵科技學(xué)院 機(jī)電工程學(xué)院, 南京 211169; 2. 中國計(jì)量大學(xué), 杭州 310018)
微小通道沿程壓力測量方法研究
王昊利1,2,*, 汪 兵2
(1. 金陵科技學(xué)院 機(jī)電工程學(xué)院, 南京 211169; 2. 中國計(jì)量大學(xué), 杭州 310018)
設(shè)計(jì)并搭建了一套微小通道沿程壓力的測量系統(tǒng),包括PMMA通道和壓力方腔、微應(yīng)變傳感器及多通道應(yīng)變儀等。利用注射泵的推進(jìn)方法提供微通道靜壓,采用FCO510型高精度微差壓計(jì)的測量值作為標(biāo)準(zhǔn)壓力,通過多通道應(yīng)變儀測量微通道方腔中各個(gè)應(yīng)變片的應(yīng)變值,從而建立標(biāo)準(zhǔn)壓力和應(yīng)變之間的標(biāo)定函數(shù)。分別對3種微壓芯片在80、70、60及50mL/min等4種不同流量下的壓力分布進(jìn)行了測量,壓力分布具有良好的線性規(guī)律。不確定度分析表明壓力誤差的相對擴(kuò)展不確定度范圍為0.15%~6.82%,測量結(jié)果的有效性和可靠性較高。
微小通道;沿程壓力;應(yīng)變測量;不確定度
微流體(尺度約為10-4~10-3m量級)技術(shù)是在微觀尺度下操作、控制及檢測復(fù)雜流體的技術(shù),近年來已經(jīng)在生物、化學(xué)、醫(yī)學(xué)以及生命科學(xué)等領(lǐng)域內(nèi)產(chǎn)生了重要影響。流場與壓力分布是開展微小流體元件優(yōu)化設(shè)計(jì)的2個(gè)重要流體動(dòng)力學(xué)參數(shù)。到目前為止,流場結(jié)構(gòu)的研究已經(jīng)有了較為完善的測量技術(shù),但微流體壓力分布測量仍然存在諸多難點(diǎn)需要進(jìn)行解決。
微流體壓力的測量之所以受到關(guān)注,一個(gè)重要的原因在于微流體條件下的結(jié)構(gòu)受力與流動(dòng)壓力密不可分。特別地,對于生物、醫(yī)學(xué)、人體生理學(xué)、微流控芯片和動(dòng)力器件所涉及的流體輸運(yùn)都需要清楚了解流場中的壓力分布。例如紅細(xì)胞在血管中的運(yùn)動(dòng),氣泡、液滴的運(yùn)動(dòng)等典型微小尺度流動(dòng)的受力和復(fù)雜三維運(yùn)動(dòng)都必須首先清楚流場中的壓力分布狀況,才能通過計(jì)算分析獲得。此外,針對微流體動(dòng)力裝置的設(shè)計(jì)與動(dòng)力性能研究也與流體壓力的分布密切相關(guān)。因此,壓力分布的測量在微流體動(dòng)力學(xué)分析以及微小動(dòng)力機(jī)械優(yōu)化設(shè)計(jì)中具有重要應(yīng)用。
遺憾的是,目前有關(guān)微流體壓力分布的精確測量還存在諸多困難。近年來國內(nèi)外針對這一問題開展了研究以發(fā)展可靠性的壓力傳感器技術(shù)為主。國內(nèi)方面,楊梅和于煒等[1]設(shè)計(jì)研制了基于梁膜結(jié)構(gòu)的微壓傳感器,能夠?qū)崿F(xiàn)0~100Pa量級的壓力測量。?,摵婉R炳和等[2]研發(fā)了以柔性襯底基陣列的微型壓力傳感器,安裝于翼型外表面實(shí)現(xiàn)壓力分布測量。國外不少學(xué)者對微尺度下壓力傳感器陣列作了深入的研究和開發(fā)。1994年,Pong等[3]開發(fā)了微通道沿程壓力測量的芯片,針對氮?dú)夂秃膺M(jìn)行了測量,結(jié)果表明沿程壓力分布和距離呈非線性關(guān)系,受努森數(shù)(Kn)的影響顯著。Ko和Liu等[4]在復(fù)雜的微通道中集成了利用一種高分子聚合物材料制作的微壓傳感器陣列來研究流體的流動(dòng)特性。Wang等[5]將成熟的硅體微加工技術(shù)和低溫PDMS工藝方法結(jié)合起來制成柔性測壓元件。Li和Luo等[6]提出了利用導(dǎo)電性的PDMS制作微壓傳感器的微流體應(yīng)用研究。Foland和Liu等[7]設(shè)計(jì)和制作了一種基于導(dǎo)模共振的壓力傳感器,能夠檢測出PDMS微通道中的微壓力的變化。Jung和Yang[8]提出了一種基于液態(tài)金屬的薄膜壓力傳感器,能夠通過軟光刻技術(shù)集成到微流控系統(tǒng)中。Tsai和Nakamura等[9]設(shè)計(jì)了一種新的用于檢測微流體裝置內(nèi)局部壓力的方法,其原理是利用聚合物PDMS薄膜材料的變形,由于變形室中的變形使得感應(yīng)區(qū)顏色強(qiáng)度發(fā)生變化,通過亮度和壓力之間的對應(yīng)關(guān)系來測量微流體裝置中局部壓力。Yeo和Yu等[10]提出了一種基于液體的薄膜微流體壓力觸覺傳感器,并且具有很高的靈敏度、穩(wěn)定性和靈活性。Song和Gillies等[11]開發(fā)了一種基于電感耦合的微流體壓力傳感器。Wu和Liao等[12]研發(fā)了一種基于離子液體電讀出壓力傳感器,通過測量微流控裝置內(nèi)壓力引起的電路特性的變化來實(shí)現(xiàn)壓力測量。
本文對微小通道內(nèi)壓力分布測量方法開展研究。制作了以微應(yīng)變片為壓力傳感元件的微壓芯片,利用高精度二次測量儀表采集通道微應(yīng)變片陣列受到流體壓力產(chǎn)生的應(yīng)變信號,在開展壓力-應(yīng)變標(biāo)定實(shí)驗(yàn)的基礎(chǔ)上獲得微通道沿程壓力分布。
設(shè)計(jì)并搭建了一套基于變形傳感原理的微壓測量裝置系統(tǒng),如圖1所示。其組成部分包括,微量注射泵、微壓芯片、無線AP、多通道應(yīng)變儀及計(jì)算機(jī)等。采用的應(yīng)變儀型號為DH3818N-2,其最高分辨率可達(dá)0.5με(微應(yīng)變),主要技術(shù)指標(biāo)由表1給出。系統(tǒng)設(shè)計(jì)原理和工作流程如下:在被測微通道兩側(cè)均勻分布應(yīng)變式微壓力傳感元件(本文采用應(yīng)變片作為測壓元件,參數(shù)等詳見下文),其引線分別接入多通道應(yīng)變儀的對應(yīng)測量端口。流體介質(zhì)(氣體或液體)在微通道中流動(dòng)產(chǎn)生的壓力使應(yīng)變片產(chǎn)生變形,其微應(yīng)變值通過應(yīng)變儀測量獲得,利用壓力-應(yīng)變標(biāo)定函數(shù)實(shí)現(xiàn)壓力測量。工作時(shí),利用微量注射泵作為動(dòng)力,以精確、均勻、恒定的速度推動(dòng)其上固定的注射器向微通道注入流體介質(zhì)。改變微量注射泵流量,應(yīng)變片在不同流量下受到流體靜壓而產(chǎn)生變形,應(yīng)變值由多通道應(yīng)變儀采集并通過無線AP與計(jì)算機(jī)進(jìn)行數(shù)據(jù)傳輸。由于應(yīng)變值并非壓力測量的目標(biāo)數(shù)據(jù),因此,需要通過標(biāo)定實(shí)驗(yàn)建立應(yīng)變與壓力之間的函數(shù)關(guān)系。
參數(shù)技術(shù)指標(biāo)測量通道數(shù)20適用應(yīng)變片電阻值50~10000Ω應(yīng)變片靈敏度系數(shù)1.0~3.0自動(dòng)修正采樣速率2Hz/通道測量應(yīng)變范圍±30000με系統(tǒng)示值誤差不大于0.5%最高分辨率0.5με
采用聚甲基丙烯酸甲酯(PMMA)材料加工尺度為毫米到亞毫米量級的微細(xì)通道,利用熱壓成型工藝制成。在通道兩側(cè)具有一定間隔的位置處加工與通道聯(lián)通的方腔,其開口面覆蓋微應(yīng)變片,制成微壓芯片。在方腔處采集與其連接點(diǎn)的壓力信號,從而獲得通道沿程各點(diǎn)處的壓力。
設(shè)計(jì)了3組不同尺寸的微壓芯片,芯片結(jié)構(gòu)與制作完成的實(shí)物圖如圖2所示,尺寸參數(shù)由表2給出。芯片由上下2層PMMA板高溫?zé)釅烘I合而成,每組微壓芯片上層板加工長直通道,通道兩側(cè)分別加工若干形狀大小相同且均勻分布的方腔,相鄰2個(gè)腔體之間的間距為5mm。矩形腔體覆蓋微應(yīng)變片,所測壓力為流體壓力與大氣壓之間的壓差,即表壓。
(a) 結(jié)構(gòu)設(shè)計(jì) (b) 實(shí)物圖
圖2 微壓芯片
采用應(yīng)變片作為壓力傳感元件,有如下諸多優(yōu)點(diǎn),包括:測量應(yīng)變的靈敏度和精確度高,性能穩(wěn)定、可靠,誤差小于1%;應(yīng)變片尺寸小,質(zhì)量輕,結(jié)構(gòu)簡單,響應(yīng)速度快;測量范圍大,可測量彈性變形,也可測量塑性變形,變形范圍可從1%~20%;適應(yīng)性強(qiáng);可在高溫、超低溫、高壓、強(qiáng)磁場等惡劣環(huán)境下使用;便于多點(diǎn)測量、遠(yuǎn)距離測量和控制等。利用應(yīng)變片感受流體靜壓作用產(chǎn)生變形,輸出應(yīng)變信號,通過標(biāo)定后轉(zhuǎn)化為流體壓力值。本文選用日本生產(chǎn)的金屬箔式應(yīng)變片,根據(jù)被測對象的尺寸要求,定制了KFG-1和KFG-5這2種型號,主要參數(shù)如表3所示。
表3 2種型號應(yīng)變片的主要參數(shù)Table 3 The main parameters of two strain gauges
根據(jù)應(yīng)變測壓原理,需建立壓力與應(yīng)變之間的函數(shù)關(guān)系,因此在進(jìn)行測量之前,先要開展標(biāo)定實(shí)驗(yàn)。標(biāo)定系統(tǒng)由微注射泵、高精度微差壓計(jì)、多通道應(yīng)變儀以及計(jì)算機(jī)等組成,如圖3所示。其中圖3(a)為標(biāo)定系統(tǒng)示意圖,圖3(b)為該系統(tǒng)的實(shí)物圖。
為了獲得穩(wěn)定的標(biāo)定壓力,將注射泵與微壓芯片通道的進(jìn)口端連接,對通道出口端進(jìn)行密封。由于通道中的壓力為數(shù)十到百帕量級,相應(yīng)的標(biāo)準(zhǔn)壓力源很難找到。因此本文采用將注射泵施加推力產(chǎn)生的靜壓作為標(biāo)準(zhǔn)壓力,解決了標(biāo)準(zhǔn)壓力源問題。由于出口端密封,通道內(nèi)壓力與注射泵推力容易實(shí)現(xiàn)平衡,測量平衡穩(wěn)定狀態(tài)下的壓力值作為當(dāng)前標(biāo)準(zhǔn)壓力。與此同時(shí),利用多通道應(yīng)變儀測量各分布點(diǎn)應(yīng)變片在該壓力時(shí)的應(yīng)變值。標(biāo)定中,通過改變注射泵推進(jìn)參數(shù)來改變標(biāo)準(zhǔn)壓力值,則基于一組標(biāo)準(zhǔn)壓力值與應(yīng)變值之間的對應(yīng)關(guān)系即可以確定,最后通過數(shù)據(jù)擬合的方法獲得壓力-應(yīng)變函數(shù)關(guān)系。
為了獲得標(biāo)準(zhǔn)壓力值,高精度微壓差計(jì)的使用是關(guān)鍵。本文采用英國Furness Controls公司生產(chǎn)的FCO510型高精度微差壓計(jì),其壓力測量范圍為0~2000Pa,分辨率高達(dá)0.01Pa。
依次對3種通道的微壓芯片開展標(biāo)定實(shí)驗(yàn)。前2組芯片標(biāo)定的最大壓強(qiáng)是50Pa,每間隔5Pa作為1個(gè)標(biāo)定區(qū)間,得到了10組標(biāo)定結(jié)果。第3組微壓芯片標(biāo)定的最大壓強(qiáng)為130Pa,每組標(biāo)定之間的間隔為10Pa。最后得到全部35組應(yīng)變值對應(yīng)壓強(qiáng)值的標(biāo)定結(jié)果。標(biāo)定曲線如圖4所示,其中圖4(a)和(b)分別有11個(gè)標(biāo)定點(diǎn),圖4(c)共13個(gè)標(biāo)定點(diǎn)。
Fig.4Calibrationcurvesandfunctions. (a)~(c)forthreekindsofchannels,respectively
在完成標(biāo)定實(shí)驗(yàn)的基礎(chǔ)上,利用微壓測量系統(tǒng),開展了微壓芯片沿程壓力的測量。本文采用氣體作為流動(dòng)介質(zhì),相應(yīng)的實(shí)驗(yàn)方法同樣適用于液體流動(dòng)。實(shí)驗(yàn)時(shí),分別取80、70、60及50mL/min等4種流量,利用注射泵進(jìn)行輸入。在不同流量下,測量應(yīng)變片的應(yīng)變值,利用圖4的標(biāo)定函數(shù),將應(yīng)變值轉(zhuǎn)化為壓力值。對各點(diǎn)所獲壓力通過最小二乘擬合,獲得沿程壓力分布曲線,并和理論結(jié)果進(jìn)行比對。
由于輸入流量較小,屬于低速流,因此采用不可壓縮管流理論進(jìn)行推導(dǎo),得到沿程壓力分布為,
式中:p為距離通道出口位置為l(mm)的壓強(qiáng)(Pa);q為氣體流量(mL/min);η是空氣的動(dòng)力粘度(Pa·s);a和b分別是通道的寬度和深度(mm);ξ為局部阻力系數(shù);ρ為密度(kg/m3)。
式(1)給出了距離通道出口位置壓強(qiáng)的函數(shù)關(guān)系。由于通道沿程存在周期性測壓腔體的分布,與完全封閉的長直通道結(jié)構(gòu)有較大差別。為了簡化分析,本文將這些腔體結(jié)構(gòu)對流動(dòng)的影響計(jì)入局部阻力損失,加上通道進(jìn)出口的局部損失,成為通道阻力損失的重要組成部分。結(jié)合實(shí)測結(jié)果,對局部阻力系數(shù)進(jìn)行調(diào)整,以修正長直通道的理論解。
對3個(gè)通道進(jìn)行理論分析,分別取局部阻力損失占沿程阻力損失的60%、51%及32%。實(shí)驗(yàn)結(jié)果、擬合直線與理論直線由圖5給出??梢钥吹?,測量點(diǎn)的壓力呈現(xiàn)良好的線性分布,修正后的理論壓力分布與之相吻合,平均偏差在5%以內(nèi)。實(shí)驗(yàn)表明應(yīng)變片的安裝、應(yīng)變測量以及標(biāo)定等實(shí)驗(yàn)關(guān)鍵環(huán)節(jié)達(dá)到了預(yù)期結(jié)果。由此說明采用應(yīng)變片方法測量微小尺度流動(dòng)的壓力具有可行性。此外,利用實(shí)驗(yàn)結(jié)果修正理論壓力分布,也為實(shí)驗(yàn)確定局部損失提供了一種新思路。
Fig.5Measurementresultsofpressuresalongthechannels.(a)~(c)forthreekindsofchannels,respectively
為進(jìn)一步開展定量分析,對測量結(jié)果進(jìn)行了不確定度的評定。不確定度評定過程如圖6所示[13]。
采用壓力誤差Δp作為被測量建立測量模型,即:
式中:n為應(yīng)變信號采樣數(shù);pi為對應(yīng)的第i次采樣的壓力測量值。
單次實(shí)驗(yàn)標(biāo)準(zhǔn)偏差為:
式中:ε是該點(diǎn)的應(yīng)變測量值。
由u1和u2進(jìn)行合成獲得合成不確定度。首先對式(1)求各變量的偏導(dǎo)數(shù)值獲得靈敏系數(shù)C,
則合成標(biāo)準(zhǔn)不確定度按下式計(jì)算:
擴(kuò)展不確定度是求測量結(jié)果區(qū)間的量,由合成標(biāo)準(zhǔn)不確定度的倍數(shù)進(jìn)行表示,兩者的乘積稱為總不確定度,計(jì)算如下:
這里的倍數(shù)即為包含因子k,本文取為k=2。
由于實(shí)驗(yàn)共涉及3組不同微通道尺寸的壓力測量,而且每組微通中距離流道進(jìn)口不同距離的壓力測量結(jié)果也不相同,因此需要對每組測壓芯片的所有測量點(diǎn)進(jìn)行不確定度的評定。根據(jù)以上公式,計(jì)算了所有測量點(diǎn)壓力誤差的不確定度,結(jié)果列入表4~6中。從表中可以看到,第1組通道的相對擴(kuò)展不確定度在0.26%~6.53%之間,第2組和第3組的相對擴(kuò)展不確定度區(qū)間分別為0.29%~6.82%及0.15%~6.12%??梢姕y量誤差不確定度分布區(qū)間較為集中,均未高于7%,體現(xiàn)了各壓力測量點(diǎn)總體的一致性,測量結(jié)果具有較高的可信度。
表4 第1組所有測點(diǎn)壓力測量不確定度評定結(jié)果Table 4 Uncertainty evaluation results of all pressure measurement points of type 1 channel
表5 第2組不同測點(diǎn)壓力測量不確定度評定結(jié)果Table 5 Uncertainty evaluation results of all pressure measurement points of type 2 channel
表6 第3組不同測點(diǎn)壓力測量不確定度評定結(jié)果Table 6 Uncertainty evaluation results of all pressure measurement points of type 3 channel
本文針對微流體分布式壓力方法展開研究,自主設(shè)計(jì)并搭建了一套分布式微壓測量系統(tǒng)。利用PMMA材料設(shè)計(jì)并制作了3組不同尺寸微通道芯片。采用高精度微壓計(jì)對該微壓芯片測量系統(tǒng)進(jìn)行了標(biāo)定,獲得了壓力和應(yīng)變之間的標(biāo)定曲線與標(biāo)定函數(shù)。分別對3種微壓芯片開展了4種不同流量的測量實(shí)驗(yàn),將理論值與實(shí)驗(yàn)結(jié)果進(jìn)行比對分析,結(jié)果表明3種微壓芯片不同流量的實(shí)驗(yàn)點(diǎn)所獲壓力呈線性函數(shù)關(guān)系,與經(jīng)過局部阻力修正的理論直線基本吻合。不確定度評定結(jié)果表明,壓力誤差的相對擴(kuò)展不確定度在0.15%~6.82%之間,測量結(jié)果具有較高的有效性和可信度。本文的實(shí)驗(yàn)方法可用于毫米及亞毫米尺度通道的氣體和液體以及多相流動(dòng)的壓力分布測量。
[1]楊梅, 于煒, 張瑩, 等. 梁-膜結(jié)構(gòu)微壓傳感器研制[J]. 實(shí)驗(yàn)流體力學(xué), 2010, 24 (2): 74-76. Yang M, Yu W, Zhang Y, et al. Development of a beam-membrane structure micro-pressure sensor[J]. Journal of Experiments in Fluid Mechanics, 2010, 24 (2): 74-76.
[2]常瑩, 馬炳和, 鄧進(jìn)軍, 等. 基于微型壓力傳感器陣列的翼面壓力分布直接測量系統(tǒng)[J]. 實(shí)驗(yàn)流體力學(xué), 2008, 22(3): 89-93. Chang Y, Ma B H, Deng J J, et al. Direct measurement system of pressure distribution on airfoil surface using micro pressure sensor array[J]. Journal of Experiments in Fluid Mechanics, 2008, 22(3): 89-93.
[3]Pong K C, Ho C M, Liu J Q, et al. Non-linear pressure distribution in uniform microchannels[C]. American Society of Mechanical Engineers, Fluids Engineering Division (Publication) FED, 1994, 197 : 51-56.
[4]Ko H S, Liu C W, Gau C, et al. Flow characteristics in a microchannel system integrated with arrays of micro-pressure sensors using a polymer material[J]. Journal of Micromechanics and Microengineering, 2008, 18(7): 75016.
[5]Wang L, Ding T, Wang P. Thin flexible pressure sensor array based on carbon black/silicone rubber nanocomposite[J]. IEEE Sensors Journal, 2009, 9(9): 1130-1135.
[6]Li H, Luo C X, Ji H, et al. Micro-pressure sensor made of conductive PDMS for microfluidic applications[J]. Microelectronic Engineering, 2010, 87(5-8): 1266-1269.
[7]Foland S, Liu K, Macfarlane D, et al. High-sensitivity microfluidic pressure sensor using a membrane-embedded resonant optical grating[J]. Sensor, 2011: 101-104.
[8]Jung T, Yang S. Highly stable liquid metal-based pressure sensor integrated with a microfluidic channel[J]. Sensors, 2014, 15(5): 11823-11829.
[9]Tsai C H D, Nakamura T, Kaneko M. An on-chip, electricity-free and single-layer pressure sensor for microfluidic applications[C]. International Conference on Intelligent Robots and Systems, IEEE, 2015.
[10]Yeo J C, Yu J, Zhao M K, et al. Wearable pressure sensors based on flexible microfluidics[J]. Lab on a Chip, 2016, 16(17): 3244-3250.
[11]Song S H, Gillies G T, Begley M R, et al. Inductively coupled microfluidic pressure meter for monitoring of cerebrospinal fluid shunt function[J]. Journal of Medical Engineering and Technology, 2012, 36(3): 156-162.
[12]Wu C Y, Liao W H, Tung Y C. A seamlessly integrated microfluidic pressure sensor based on an ionic liquid electrofluidic circuit[C]. International Conference on Micro Electro Mechanical Systems. 2011: 1087-1090.
[13]張有康, 甘蓉. 壓力傳感器測量中不確定度的評定[J]. 中國測試技術(shù), 2005, 31(6): 25-26.Zhang Y K, Gan R. Evaluation of strain gauge measurement uncertainty of pressure transducer[J]. China Measurement Technology, 2005, 31(6): 25-26.
Researchofthemeasurementmethodforthepressuredistributionalongthemicro/mini-channel
Wang Haoli1,2,*, Wang Bing2
(1. College of Electrical Engineering, Jinling Institute of Technology, Nanjing 211169, China; 2. China Jiliang University, Hangzhou 310018, China)
In this study, a system was designed and built up for the pressure measurement along micro/mini-channel, which mainly includes micro/mini-channel and square pressure cavities on PMMA chip, micro-strain sensors and multi-channel strain instruments. The static pressure in micro/mini-channel was provided by the syringe pump, and the high precision pressure values measured by the micromanometer of FCO510 was employed as the standard pressure. The strain values from the strain sensors installed on the square pressure cavities were obtained by multi-channels strain gauge, and the calibration functions between the standard pressure and the strain were established. The pressure distributions of three kinds of micro pressure chips were measured under flow rates of 80, 70, 60 and 50mL/min, respectively. The pressures have good linear distribution.The uncertainty analysis indicates that the relative uncertainty of the pressure error is between 0.15% and 6.82%. The validity and reliability of pressure measurement are high.
micro/mini-channel; pressure along the channel; strain measurement; uncertainty
2017-04-25;
2017-06-21
國家自然科學(xué)基金(11472261, 11172287); 金陵科技學(xué)院高層次人才啟動(dòng)項(xiàng)目(2016)
*通信作者 E-mail: whl@cjlu.edu.cn
WangHL,WangB.Researchofthemeasurementmethodforthepressuredistributionalongthemicro/mini-channel.JournalofExperimentsinFluidMechanics, 2017, 31(6): 56-61. 王昊利, 汪 兵. 微小通道沿程壓力測量方法研究. 實(shí)驗(yàn)流體力學(xué), 2017, 31(6): 56-61.
1672-9897(2017)06-0056-07
10.11729/syltlx20170050
O352
A
王昊利(1972-),男,山西芮城人,博士,教授。研究方向:微納尺度流動(dòng)可視化實(shí)驗(yàn)技術(shù)。通信地址:江蘇南京江寧區(qū)弘景大道99號金陵科技學(xué)院機(jī)電工程學(xué)院(211169)。E-mail: whl@cjlu.edu.cn
(編輯:張巧蕓)