季楚凌,李長(zhǎng)俊
(西南石油大學(xué) 石油與天然氣工程學(xué)院,四川 成都 610500)
沖蝕磨損是指材料在受到小而稀松的流動(dòng)粒子的沖擊時(shí)其表面發(fā)生破壞的一類磨損現(xiàn)象,其定義可以描述為固體表面同含有固體粒子的流體接觸做相對(duì)運(yùn)動(dòng)時(shí),其表面材料所發(fā)生的損耗[1]。對(duì)于特定的油氣集輸管道來(lái)說(shuō),由于在開(kāi)采及其它環(huán)節(jié)中會(huì)產(chǎn)生一定的雜質(zhì),管輸介質(zhì)中常含有一定數(shù)量的固體顆粒,這些含砂流體沖擊管道內(nèi)壁并對(duì)其造成一定程度的沖蝕磨損,最終影響正常的生產(chǎn)并有可能造成管道的泄漏。而之前眾多專家學(xué)者的研究結(jié)果表明,當(dāng)含砂流體流經(jīng)彎管、異徑接頭、三通、閥門等管道局部件時(shí),其自身流域的物理形狀會(huì)使得固粒在更大程度上偏離流線,從而對(duì)壁面造成更為嚴(yán)重的沖蝕磨損[2-5]。
目前,在油氣集輸管道上常用的耐磨方法是涂敷耐磨涂層,且主要使用有機(jī)涂層,但是有機(jī)涂層在實(shí)際生產(chǎn)中所起到的耐磨作用較為有限,并不能對(duì)集輸管道內(nèi)壁起到很好的保護(hù)作用,而近年來(lái),不同材料耐磨涂層的制備及應(yīng)用情況不斷被報(bào)道。本文較為系統(tǒng)地總結(jié)了各類耐磨涂層的耐磨機(jī)理及在管道上的實(shí)際使用情況,在此基礎(chǔ)上進(jìn)一步介紹了提高各涂層耐磨性能的技術(shù)改進(jìn)方法,并對(duì)其在油氣集輸管道上的應(yīng)用進(jìn)行了展望。
現(xiàn)今國(guó)內(nèi)外油氣管道(包括油氣集輸管道)上常用的有機(jī)耐磨涂層的材料主要有粉末環(huán)氧樹(shù)脂、聚氨酯、環(huán)氧酚醛樹(shù)脂等[6],而其中比較著名的產(chǎn)品包括英國(guó)伍德公司的COPON、德國(guó)杜邦公司的Permecor 及我國(guó)自主研制的AW-01 等[7]。有機(jī)耐磨涂層在固體顆粒沖擊下的沖蝕磨損機(jī)理與塑形、脆性材料涂層均不相同,相關(guān)的研究結(jié)果表明[8-9],當(dāng)固體顆粒以一定的速度沖擊涂敷有機(jī)涂層的靶材時(shí),涂層表面所受的剪切應(yīng)力會(huì)使得其中心區(qū)域發(fā)生環(huán)狀開(kāi)裂,并引起初步的涂層與基體分層現(xiàn)象;而在此后,無(wú)論是在較小角度還是較大角度下,固體顆粒的持續(xù)沖擊會(huì)使得基體與涂層之間的結(jié)合部受到較大的壓應(yīng)力作用,當(dāng)壓應(yīng)力大于塑形管材基體的屈服強(qiáng)度時(shí),塑形基體與有機(jī)涂層均會(huì)發(fā)生屈服現(xiàn)象,而兩者屈服程度的不同則會(huì)導(dǎo)致之前分層區(qū)域的不斷擴(kuò)大,并最終引起有機(jī)涂層的脫落。
對(duì)于有機(jī)涂層制備方法的研究始于20 世紀(jì)40年代,并且經(jīng)過(guò)多年的發(fā)展,現(xiàn)有的有機(jī)耐磨涂層主要由基料、填料及固化劑組成,其中基料主要選用環(huán)氧樹(shù)脂,而填料直接決定著涂層的性能,因此可以通過(guò)添加不同的填料來(lái)改善原有有機(jī)涂層的耐磨性[10]。Wetzel 等[11]在環(huán)氧樹(shù)脂基料中添加了體積分?jǐn)?shù)為1% ~2%的微米級(jí)CaSiO3和納米級(jí)Al2O3顆粒,并通過(guò)測(cè)試證明了在加入上述粉末增強(qiáng)相后有機(jī)涂層的耐磨性能可大幅增強(qiáng)。而在耐磨性能及機(jī)理的研究方面,除了需要進(jìn)行相應(yīng)的實(shí)驗(yàn)測(cè)試外,Wang 等[12]于2009 年、Gong 等[13]于2012 年均采用將有限元算法及無(wú)網(wǎng)格算法相結(jié)合的方法來(lái)建立沖蝕磨損的耦合模型,運(yùn)用SPH 方法的光滑質(zhì)點(diǎn)對(duì)涂層及金屬基體的受沖擊中心區(qū)域進(jìn)行建模,而其余部分均采用傳統(tǒng)的有限元方法建模,通過(guò)將數(shù)值模擬結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比,可以發(fā)現(xiàn)這種方法可以解決原先采用有限元法時(shí)有機(jī)涂層受沖擊區(qū)域變形量較大時(shí)網(wǎng)格扭曲的問(wèn)題,因而對(duì)于有機(jī)涂層沖蝕磨損動(dòng)態(tài)過(guò)程的模擬更為準(zhǔn)確。
總的來(lái)說(shuō),雖然現(xiàn)有的有機(jī)涂層具有涂裝工藝簡(jiǎn)單、技術(shù)成熟等優(yōu)點(diǎn),但是,較之陶瓷涂層及金屬合金涂層,有機(jī)涂層的耐磨性仍然較差。因此,有機(jī)涂層作為耐磨材料主要適用于固粒濃度較小的長(zhǎng)輸管道(如西氣東輸管道),而對(duì)于固粒濃度較大的油氣集輸管道,應(yīng)在使用前對(duì)有機(jī)涂層在特定條件下的耐磨性能進(jìn)行評(píng)估,否則極易造成其成片分層脫落。
陶瓷涂層屬于典型的脆性材料涂層,因此其沖蝕磨損機(jī)理包括應(yīng)力疲勞斷裂及脆性斷裂兩種模型,并且沖蝕磨損機(jī)制的轉(zhuǎn)化主要與顆粒的法向沖擊速度相關(guān)[14-15]。如圖1 固體顆粒沖擊陶瓷材料表面示意圖中所示,當(dāng)法向沖擊速度較小時(shí),固體顆粒只能對(duì)陶瓷材料造成彈塑性變形,而并不能造成其磨損;當(dāng)固粒的法向沖擊速度達(dá)到一定的值時(shí),陶瓷材料由于受到固體顆粒長(zhǎng)時(shí)間的反復(fù)沖擊而使得其應(yīng)力超過(guò)本身的強(qiáng)度,橫向及徑向裂紋得以產(chǎn)生并以一定的速度擴(kuò)展,最終造成陶瓷材料的沖蝕磨損,此即應(yīng)力疲勞斷裂模型;而當(dāng)固粒的法向沖擊速度較大時(shí),陶瓷材料在少數(shù)幾次沖擊下就會(huì)產(chǎn)生橫向及徑向裂紋,并且這些裂紋的擴(kuò)展速度較快,最終造成材料流失,此即脆性斷裂模型。
圖1 固體顆粒沖擊陶瓷材料表面示意圖Fig.1 Illustration of solid particle impacting ceramic material surface
現(xiàn)有的陶瓷涂層的制備方法有很多種,其中,Shin 等[16]利用氣相沉積法制備得到了陶瓷涂層,并對(duì)其微觀組織結(jié)構(gòu)和機(jī)械耐磨性能進(jìn)行了相應(yīng)的研究。朱昱等[17]綜述了激光表面改性技術(shù)在金屬基陶瓷涂層制備中的應(yīng)用,指出利用激光表面改性技術(shù)在金屬基體表面上能制備出耐磨、耐腐蝕和耐高溫的陶瓷涂層,并對(duì)其未來(lái)的發(fā)展方向進(jìn)行了展望。對(duì)于陶瓷涂層的制備,除了上述兩種方法外,還有熱噴涂法、離心法、包埋法等多種方法。而對(duì)于陶瓷涂層耐磨性能及機(jī)理的研究主要是運(yùn)用實(shí)驗(yàn)研究與數(shù)值模擬相結(jié)合的方法來(lái)進(jìn)行,例如,Griffin 等[18]建立了多顆粒沖擊帶有氧化鋁陶瓷涂層的金屬基體的數(shù)值模型,運(yùn)用拉伸失效準(zhǔn)則作為判斷涂層單元是否失效的基準(zhǔn),通過(guò)統(tǒng)計(jì)失效單元來(lái)計(jì)算沖蝕磨損速率,并由此評(píng)價(jià)涂層的耐磨性能,通過(guò)將計(jì)算結(jié)果與實(shí)驗(yàn)數(shù)據(jù)進(jìn)行對(duì)比可以證明數(shù)值計(jì)算模型的準(zhǔn)確性。
由于陶瓷涂層的制備工藝相對(duì)較為復(fù)雜,因此在油氣集輸管線上進(jìn)行全線的涂裝較為困難,但對(duì)于那些處在顆粒低角沖擊(此時(shí)法向沖擊速度較小)下且沖蝕磨損情況較為嚴(yán)重的局部件,例如彎管、異徑接頭、大開(kāi)度下的閥門等,可以在其沖蝕磨損集中區(qū)制備陶瓷涂層,預(yù)計(jì)將會(huì)取得較好的耐磨效果。
金屬合金涂層屬于塑形材料涂層,其沖蝕磨損機(jī)理可以通過(guò)微切削模型[19]來(lái)進(jìn)行解釋。如圖2固體顆粒沖擊金屬合金材料表面示意圖中所示,對(duì)于塑性材料來(lái)說(shuō),當(dāng)受到高角沖擊時(shí),其自身的彈性應(yīng)變能夠吸收沖擊能量,從而使其材料表面免受損失;反之,當(dāng)顆粒對(duì)材料表面進(jìn)行低角沖擊時(shí),沖擊和摩擦產(chǎn)生的表面應(yīng)力會(huì)引起裂紋并造成材料損失。若顆粒棱角比較銳利(圖2(a)),則其向前的滑動(dòng)會(huì)使金屬合金表面的材料就像被車刀車削一樣從磨粒前方被去除;而若顆粒棱角不夠銳利(圖2(b)),那么金屬合金材料表面會(huì)發(fā)生犁溝變形,固體顆粒一邊向前推擠材料,一邊將材料犁向溝槽兩側(cè),同樣造成材料損失。
圖2 固體顆粒沖擊金屬合金材料表面示意圖Fig.2 Illustration of solid particle impacting metal alloy material surface
現(xiàn)有金屬合金涂層的制備方法也有很多,其中,Zhou 等[20]利用超音速火焰噴涂(HVOF)技術(shù)制備得到了Fe48Cr15Mo14C15B6Y2合金涂層,通過(guò)研究發(fā)現(xiàn)該涂層呈現(xiàn)致密的層狀結(jié)構(gòu),并且這種結(jié)構(gòu)上的優(yōu)點(diǎn)使得其具有優(yōu)良的耐磨性能。黃標(biāo)等[21]利用激光熔覆技術(shù)在45 號(hào)鋼表面制備了FeCoCrxNiB 高熵合金涂層,并對(duì)其組織結(jié)構(gòu)和耐磨性進(jìn)行了一系列的觀察和測(cè)試,研究結(jié)果表明FeCoCr0.5NiB 涂層的耐磨性最高。除了上述介紹的制備方法外,還有熔融法、沉積法、電弧噴涂、等離子噴涂等多種方法。而對(duì)金屬合金涂層耐磨性能及機(jī)理的研究主要是運(yùn)用實(shí)驗(yàn)研究與有限元顯示動(dòng)力學(xué)模擬相結(jié)合的方法來(lái)進(jìn)行的,其中,Wang 等[22]運(yùn)用有限元方法中的顯示動(dòng)力學(xué)模塊分析了沖蝕速度、沖蝕角、沖蝕顆粒尺寸等不同影響因素對(duì)于塑形及脆性材料沖蝕磨損特性的影響,通過(guò)與前人實(shí)驗(yàn)數(shù)據(jù)的對(duì)比驗(yàn)證了數(shù)值模型的準(zhǔn)確性,并由此進(jìn)一步分析塑形及脆性材料各自的沖蝕磨損機(jī)理及相關(guān)的影響因素。
與陶瓷涂層類似,現(xiàn)有工藝的復(fù)雜性使得在油氣集輸管線上進(jìn)行金屬合金涂層的全線涂裝較為困難,此外,若在集輸管道上運(yùn)用金屬合金涂層時(shí),需要重點(diǎn)分析其在特定條件下的耐腐蝕性能,尤其是當(dāng)集輸管道中輸送CO2、H2S 等含量較高的強(qiáng)腐蝕性介質(zhì)時(shí),不宜使用金屬合金涂層。但對(duì)于那些處在較弱腐蝕性介質(zhì)的環(huán)境中、顆粒高角沖擊下且沖蝕磨損情況較為嚴(yán)重的局部件,例如T 型三通、小開(kāi)度下的閥門等,若在其沖蝕磨損集中區(qū)制備金屬合金涂層預(yù)計(jì)會(huì)取得較好的耐磨效果。
材料改性處理是從材料本身的耐磨性能方面入手,通過(guò)對(duì)原涂層材料進(jìn)行改性來(lái)提高涂層的耐磨性或解決某些涂層在管道內(nèi)壁上制備工藝復(fù)雜的難題。有機(jī)涂層的改性處理除了之前文獻(xiàn)[11]中所介紹的添加CaSiO3、Al2O3等現(xiàn)今常用的固體粉末增強(qiáng)相外,一些專家學(xué)者嘗試使用橡膠等彈性體材料對(duì)有機(jī)涂層進(jìn)行處理,利用橡膠顆粒脫膠或斷裂后所形成的孔洞的塑性體膨脹和顆粒或孔洞所誘發(fā)的剪切屈服形變使原有有機(jī)涂層的韌性大大增強(qiáng),于晶[23]將橡膠材料混入有機(jī)環(huán)氧樹(shù)脂材料中,開(kāi)發(fā)出高強(qiáng)耐磨蝕復(fù)合材料,實(shí)驗(yàn)結(jié)果表明這種復(fù)合材料涂層可以有效緩解過(guò)流部件的磨蝕問(wèn)題;對(duì)于陶瓷涂層來(lái)說(shuō),對(duì)其進(jìn)行的改性處理主要包括提高耐磨性及簡(jiǎn)化制備工藝兩方面,在提高耐磨性方面,Richard 等[24]運(yùn)用熱噴涂法在金屬管材基體表面制備得到了Al2O3-TiO2納米陶瓷結(jié)構(gòu)涂層,并且納米陶瓷涂層表面結(jié)構(gòu)的致密性使其耐磨性能較傳統(tǒng)陶瓷涂層提高很多,但是在實(shí)際應(yīng)用時(shí)還需要解決在噴涂過(guò)程中納米顆粒易沉積以及如何抑制納米晶在噴涂過(guò)程中長(zhǎng)大等問(wèn)題。在簡(jiǎn)化制備工藝方面,為使陶瓷涂層能夠方便高效地涂裝于集輸管線上,可以考慮將陶瓷粉末加入膠粘劑溶液中制成粘稠狀膠體后再進(jìn)行涂裝,吳瑞娟[25]將骨料MgO、TiO2、Al2O3處理后加入計(jì)量的膠粘劑磷酸二氫鋁溶液中用強(qiáng)力電動(dòng)攪拌器分散至粘稠狀,混合均勻后制得陶瓷涂料,并將其涂覆在尺寸為D108 ×5 的X80 管線鋼上,測(cè)試結(jié)果表明這種陶瓷涂料的附著力與抗沖擊性均能達(dá)到相關(guān)的要求;而對(duì)于金屬合金涂層,常在其中混入WC、TiC 等無(wú)氧碳化物陶瓷粉末以增強(qiáng)其耐磨性,Zhou 等[26]運(yùn)用激光技術(shù)制備得到了鐵基WC 復(fù)合涂層,并觀察了其微觀組織結(jié)構(gòu),測(cè)試結(jié)果表明,在相同條件下涂覆該涂層后AISI 1045 鋼材的磨損速率比原來(lái)減小了2 倍。
非光滑耐磨表面是通過(guò)一定的成型方法加工得到的具有一定耐磨減阻功能的仿生表面,它通過(guò)改變?cè)泄腆w顆粒的運(yùn)動(dòng)軌跡來(lái)達(dá)到提高涂層耐磨性的目的。仿生非光滑耐磨表面在各個(gè)領(lǐng)域發(fā)展迅速,且主要的表面形態(tài)有橫槽型、縱槽型、凹坑型、凸包型等[27-28],但是在油氣管道上的應(yīng)用仍然較少,主要原因是在管道內(nèi)壁上加工這類表面還比較困難,而近三年來(lái)哈爾濱工程大學(xué)趙剛教授的課題組對(duì)油氣輸送管道中涂層上的非光滑表面的成型機(jī)理進(jìn)行了研究,并成功研制了管道內(nèi)壁加工機(jī)器人,該機(jī)器人通過(guò)運(yùn)用熱輥壓的方法在天然氣輸送管道的有機(jī)涂層上加工出微米級(jí)的凹坑型耐磨表面[29-30],預(yù)計(jì)該方法可以提高原有機(jī)涂層的耐磨減阻性能,但其具體的耐磨減阻效果還需實(shí)驗(yàn)驗(yàn)證。
各類耐磨涂層有其各自不同的耐磨機(jī)理,因此日后在油氣集輸管線上進(jìn)行耐磨處理時(shí)應(yīng)根據(jù)各集輸管道的實(shí)際沖蝕磨損環(huán)境、各類型管道局部件的不同特點(diǎn)、各類耐磨涂層的特點(diǎn)來(lái)選用最為合適的耐磨涂層,而為了進(jìn)一步拓展不同種類耐磨涂層在油氣集輸管道中的應(yīng)用,應(yīng)積極開(kāi)展以下工作:
(1)對(duì)陶瓷涂層及金屬合金涂層的制備方法進(jìn)行改進(jìn),以使其能在油氣集輸管線上進(jìn)行大范圍乃至全線的涂裝。
(2)深入研究經(jīng)改性處理后的復(fù)合材料耐磨涂層的沖蝕磨損有限元模型,通過(guò)數(shù)值模擬與實(shí)驗(yàn)研究相結(jié)合的方法從機(jī)理上分析并解釋其耐磨性能提高的原因。
(3)進(jìn)一步研究橫槽型、縱槽型、凸包型等其它非光滑表面的成型方法,通過(guò)現(xiàn)場(chǎng)測(cè)試檢驗(yàn)其在油氣集輸管線上的耐磨減阻效果,并優(yōu)選出耐磨性能最好的非光滑表面形態(tài)。
[1] 鄧建新,丁澤良. 陶瓷噴嘴及其沖蝕磨損[M]. 北京:科學(xué)出版社,2009.
[2] Habib M A,Ben-Mansour R,Badr H M,et al.Erosion and penetration rates of a pipe protruded in a sudden contraction[J].Computers & Fluids,2008,37(2):146-160.
[3] Zhang H,Tan Y,Yang D,et al.Numerical investigation of the location of maximum erosive wear damage in elbow:Effect of slurry velocity,bend orientation and angle of elbow[J].Powder Technology,2012,217(2):467-476.
[4] Wong C Y,Solnordal C,Swallow A,et al. Experimental and computational modelling of solid particle erosion in a pipe annular cavity[J].Wear,2013,303(1/2):109-129.
[5] Zhu H,Pan Q,Zhang W,et al.CFD simulations of flow erosion and flow-induced deformation of needle valve:Effects of operation,structure and fluidparameters[J].Nuclear Engineering and Design,2014,273(7):396-411.
[6] 沈遠(yuǎn).輸氣管道內(nèi)涂層制備及其減阻性能研究[D].上海:華東理工大學(xué)材料科學(xué)與工程學(xué)院,2014.
[7] 林竹,秦延龍,黃驍卓,等. 天然氣管道減阻耐磨涂料的研究和應(yīng)用進(jìn)展[J]. 腐蝕與防護(hù),2003,24(5):206-209.
[8] Mabrouki T,Raissi K,Cornier A. Numerical simulation and experimental study of the interaction between a pure high-velocity waterjet and targets:contribution to investigate the decoating process[J]. Wear,2000,239(2):260-273.
[9] 王宇飛.耐磨耐蝕層狀復(fù)合管的制備及其沖蝕模擬分析[D].上海:復(fù)旦大學(xué)材料科學(xué)系,2008.
[10]周先輝,孫友松,魏良模,等. 新一代環(huán)氧潤(rùn)滑耐磨涂層研制與應(yīng)用中需解決的若干問(wèn)題[J]. 材料保護(hù),2007,40(3):56-60.
[11]Wetzel B,Haupert F,Zhang M Q. Epoxy nanocomposites with high mechanical and tribological performance[J].Composites Science & Technology,2003,63 (14):2055-2067.
[12]Wang Y,Yang Z. A coupled finite element and meshfree analysis of erosive wear[J]. Tribology International,2009,42(2):373-377.
[13] Gong Y,Yang Z,Wang Y. Impact simulation on ductile metal pipe with polymer coating by a coupled finite element and meshfree method[J].Journal of Failure Analysis & Prevention,2012,12(3):267-272.
[14]孫軍龍.新型B4C 基復(fù)合陶瓷噴嘴研究開(kāi)發(fā)及其沖蝕機(jī)理研究[D].濟(jì)南:山東大學(xué)機(jī)械工程學(xué)院,2007.
[15]Deng J,Wu F,Zhao J. Wear mechanisms of gradient ceramic nozzles in abrasive air-jetmachining[J]. International Journal of Machine Tools and Manufacture,2007,47(12):2031-2039.
[16] Shin J,Choi K,Wang T,et al. Microstructure evolution and mechanical properties of Ti-B-N coatings deposited by plasma-enhanced chemical vapor deposition[J].Transactions of Nonferrous Metals Society of China,2012,22(35):722-728.
[17]朱昱,張宇,王小美,等. 激光表面改性在金屬基陶瓷涂層制備中的應(yīng)用[J]. 現(xiàn)代化工,2014,34(10):64-67.
[18] Griffin D,Daadbin A,Datta S. The development of a three-dimensional finite element model for solid particle erosion on an alumina scale/MA956 substrate[J].Wear,2004,256(9/10):900-906.
[19]Finnie I. Erosion of surfaces by solid particles[J].Wear,1960,3(2):87-103.
[20]Zhou Z,Wang L,He D Y,et al. Microstructure and wear resistance of Fe-based amorphous metallic coatings prepared by HVOF thermal spraying[J].Journal of Thermal Spray Technology,2010,19(6):1287-1293.
[21]黃標(biāo),張沖,程虎,等. 激光熔覆FeCoCrxNiB 高熵合金涂層的組織結(jié)構(gòu)與耐磨性[J]. 中國(guó)表面工程,2014,27(6):82-88.
[22]Wang Y,Yang Z.Finite element model of erosive wear on ductile and brittlematerials[J]. Wear,2008,265(5/6):871-878.
[23]于晶.高強(qiáng)耐磨蝕復(fù)合材料的合成與應(yīng)用研究[D].濟(jì)南:山東大學(xué)化學(xué)與化工學(xué)院,2010.
[24]Richard C,Kowandy C,Landoulsi J,et al. Corrosion and wear behavior of thermally sprayed nano ceramic coatings on commercially pure Titanium and Ti-13Nb-13Zr substrates[J]. International Journal of Refractory Metals &Hard Materials,2010,28(1):115-123.
[25]吳瑞娟.一種天然氣管道內(nèi)防腐涂料的制備及性能研究[D]. 成都:西南石油大學(xué)材料科學(xué)與工程學(xué)院,2014.
[26] Zhou S,Dai X,Zheng H. Microstructure and wear resistance of Fe-based WC coating by multi-track overlapping laser induction hybrid rapidcladding[J]. Optics & Laser Technology,2012,44(1):190-197.
[27]Ge Chao,Jiang Jialian,Ren Luquan.Gas-solid erosion on bionic configuration surface[J].Journal of Wuhan University of Technology-Mater Sci Ed,2011,26(2):305-310.
[28]Zhang J,Han Z,Yin W,et al.Numerical experiment of the solid particle erosion of bionic configuration blade of centrifugalfan[J]. Acta Metallurgica Sinica (English Letters),2013,26(1):16-24.
[29]饒軍. 管道內(nèi)壁減阻結(jié)構(gòu)加工機(jī)器人關(guān)鍵技術(shù)研究[D].哈爾濱:哈爾濱工程大學(xué)機(jī)電工程學(xué)院,2012.
[30]于偉波. 油氣輸送管道仿生減阻表面成型機(jī)理研究[D].哈爾濱:哈爾濱工程大學(xué)機(jī)電工程學(xué)院,2013.