董 星， 楊 威， 李光煜

(黑龍江科技大學(xué) 機(jī)械工程學(xué)院，哈爾濱150022)

噴丸強(qiáng)化作為一種表面強(qiáng)化技術(shù)廣泛應(yīng)用于眾多工程領(lǐng)域，以提高碳鋼、合金鋼、鋁合金和鈦合金等材料的疲勞抗力［1－5］。后混合水射流噴丸強(qiáng)化作為一種新型的噴丸強(qiáng)化方法而受到人們的關(guān)注。而噴頭作為射流噴丸中實(shí)現(xiàn)能量轉(zhuǎn)換的關(guān)鍵元件，其內(nèi)流場(chǎng)多相流混合機(jī)制的好壞是決定噴丸質(zhì)量的關(guān)鍵因素之一［6－10］。因此，筆者采用FLUENT 軟件對(duì)后混合水射流噴丸噴頭內(nèi)流場(chǎng)進(jìn)行數(shù)值模擬，探究其內(nèi)流場(chǎng)分布規(guī)律，為提高噴丸質(zhì)量奠定技術(shù)基礎(chǔ)。

1 模型與方法

1.1 幾何模型

采用Pro/E 軟件建立后混合水射流噴丸噴頭三維實(shí)體結(jié)構(gòu)幾何模型。模型的總長(zhǎng)為93 mm，其中，水噴嘴入口圓錐收斂段長(zhǎng)為4 mm，入口直徑為1 mm，水噴嘴出口圓柱段長(zhǎng)為2 mm，出口直徑為0.3 mm;彈丸供丸管直徑為6 mm，彈丸供丸管與混合室軸線(xiàn)成45°角;混合室長(zhǎng)為23 mm，直徑為6.4 mm;彈丸噴嘴圓錐收斂段長(zhǎng)9 mm，彈丸噴嘴出口段長(zhǎng)為64 mm，圓柱段直徑為0.8 mm。

1.2 有限元模型

將幾何模型導(dǎo)入FLUENT 軟件中的GAMBIT 模塊進(jìn)行前處理，采用Hex 和Tet 網(wǎng)格對(duì)其進(jìn)行網(wǎng)格劃分，網(wǎng)格劃分后的單元數(shù)為191 292 個(gè)，節(jié)點(diǎn)數(shù)為101 572 個(gè)，有限元模型如圖1 所示。

圖1 有限元模型Fig．1 Finite element model

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

由于所研究后混合水射流噴丸過(guò)程中彈丸粒子的體積分?jǐn)?shù)大于10%，因此，采用歐拉模型進(jìn)行求解，并補(bǔ)充湍流模型使方程封閉。

1.3.1 控制方程

連續(xù)性方程。假設(shè)第一相為液相水，第二相為固相彈丸，其連續(xù)性方程為

式中:αl、αs——水和彈丸的體積分?jǐn)?shù);

ρl、ρs——水和彈丸的密度;

vl、vs——水和彈丸的速度向量。

動(dòng)量方程。歐拉模型中液相水和固相彈丸的動(dòng)量方程分別為

學(xué)生作為教學(xué)的主體，只有充分重視學(xué)生的主體地位，發(fā)揮其能動(dòng)性，才能更好調(diào)動(dòng)其學(xué)習(xí)熱情，使學(xué)生積極參與到教學(xué)中來(lái)，進(jìn)而實(shí)現(xiàn)高效教學(xué)。抽象思維能力是學(xué)生應(yīng)具備的一項(xiàng)素質(zhì)，同時(shí)也是物理教學(xué)的重要目標(biāo)，而實(shí)現(xiàn)這一目標(biāo)，關(guān)鍵在于重視學(xué)生的主體作用，圍繞學(xué)生來(lái)開(kāi)展教學(xué)，讓學(xué)生參與到課堂教學(xué)中來(lái) 。如，在教習(xí)“平拋運(yùn)動(dòng)”這一知識(shí)時(shí)，教師可以讓學(xué)生去做平拋運(yùn)動(dòng)的實(shí)驗(yàn)，讓學(xué)生記錄實(shí)驗(yàn)過(guò)程，并提出自己在實(shí)驗(yàn)過(guò)程中遇到的問(wèn)題，然后就學(xué)生遇到的問(wèn)題在課堂上展開(kāi)互動(dòng)討論。實(shí)驗(yàn)既鍛煉了學(xué)生的實(shí)際操作能力，同時(shí)也有效地發(fā)展了學(xué)生的思維，有利于學(xué)生思維能力的提升。

式中:p——水和彈丸所有相共享的壓力;

τl、τs——水和彈丸的壓力應(yīng)變量;

g——重力加速度;

Kls、Ksl——水和彈丸之間的動(dòng)量交換系數(shù)，Kls=Ksl;

Fl、Fs——水和彈丸的外部體積力;

Flift，l、Flift，s——水和彈丸的升力;

ps——彈丸的壓力。

標(biāo)準(zhǔn)k －ε 兩方程湍流模型的湍動(dòng)能k 和耗散率ε 方程分別為

式中:k——單位質(zhì)量流體湍動(dòng)能;

μt——湍動(dòng)能黏性系數(shù)，μt=Cμρk2/ε;

Cμ——系數(shù);

ε——單位質(zhì)量流體湍動(dòng)能耗散率;

C1、C2——系數(shù);

σk、σε——k 和ε 的湍流Prandtl 數(shù)。

模擬時(shí)，取Cμ=0.09，C1=1.44，C2=1.92，σk=1.0，σε=1.3。

1.4 計(jì)算方法

選擇三維單精度求解器，采用一階精度迎風(fēng)差分格式，基于液固兩相體積分?jǐn)?shù)的湍流控制方程，并通過(guò)SIMPLE 算法求解壓力速度耦合。計(jì)算過(guò)程中壓力、動(dòng)量、k 和ε 等的亞松弛因子均取默認(rèn)值，收斂判據(jù)為流場(chǎng)迭代的殘差值R≤10－3。

1.5 邊界條件

將固相彈丸視為擬流體。

(1)入口邊界條件:對(duì)水噴嘴入口取為壓力入口邊界條件，給定入口全壓，湍動(dòng)能k 按來(lái)流平均動(dòng)能的0.5%給出，湍動(dòng)能耗散率由式ε=k3/2/(0.2d)計(jì)算出，d 為水噴嘴入口直徑。固相彈丸材料為玻璃，彈丸直徑為0.25 mm，彈丸入口取為速度入口邊界條件，速度大小為5 m/s，固相體積分?jǐn)?shù)為13%。

(2)出口邊界條件:取壓力出口邊界條件，給定靜壓力為當(dāng)?shù)卮髿鈮?01 325 Pa。

(3)壁面邊界條件:對(duì)連續(xù)相，在固體壁面上速度滿(mǎn)足無(wú)滑移條件，近壁面區(qū)域采用標(biāo)準(zhǔn)壁面函數(shù)法。

2 模擬結(jié)果及分析

2.1 噴丸壓力對(duì)內(nèi)流場(chǎng)的影響

圖2 不同噴丸壓力下內(nèi)部水流場(chǎng)軸向速度分布云圖Fig．2 Inside water flow field under different peening pressure axial velocity distribution nephogram

圖3 不同噴丸壓力下內(nèi)部彈丸流場(chǎng)的軸向速度分布云圖Fig．3 Inside pill flow field under different peening pressure axial velocity distribution nephogram

圖2 和圖3 分別給出了在彈丸噴嘴長(zhǎng)l 為55 mm，噴丸壓力p 分別為6、8、10、12 和14 MPa 時(shí)內(nèi)部水流場(chǎng)和彈丸流場(chǎng)的軸向速度va分布云圖。

由圖2 可知，在不同的噴丸壓力下，噴頭內(nèi)部水流場(chǎng)的軸向速度變化規(guī)律基本相同，水噴嘴入口處軸向速度隨噴丸壓力的增加而增大，在水噴嘴入口段的圓錐收斂段軸向速度有明顯的變化，在混合室內(nèi)部軸向速度基本沒(méi)有變化，在混合室的圓錐收斂段軸向速度開(kāi)始增加，在彈丸噴嘴出口段水流場(chǎng)速度繼續(xù)增加，在近壁面附近由于噴嘴壁面的黏滯力作用使得近壁面速度遠(yuǎn)小于軸線(xiàn)速度，從壁面沿徑向方向射流速度依次增大，且軸向速度分布呈現(xiàn)對(duì)稱(chēng)性。隨著噴丸壓力的增加，水流場(chǎng)在彈丸噴嘴出口處的軸向速度逐漸增大，在噴丸壓力分別為6、8、10、12 和14 MPa 的情況下，水流場(chǎng)在彈丸噴嘴出口處的最大軸向速度分別為55.0、60.0、63.3、65.2 和68.4 m/s。

由圖3 可知，彈丸流場(chǎng)的軸向速度變化規(guī)律與水流場(chǎng)的基本相同，彈丸流場(chǎng)在彈丸噴嘴出口處的最大軸向速度分別為49.9、52.2、56.9、62.0 和65.5 m/s。通過(guò)比較可知，在相同的噴丸壓力下彈丸噴嘴出口處水流場(chǎng)的軸向速度大于彈丸流場(chǎng)的軸向速度。

噴丸過(guò)程中彈丸噴嘴出口處的軸向動(dòng)壓強(qiáng)也隨著噴丸壓力的增加而增大，在噴丸壓力分別為6、8、10、12 和14 MPa 的情況下，水流場(chǎng)在彈丸噴嘴出口處的最大軸向動(dòng)壓強(qiáng)分別為1.56、1.79、1.97、2.12和2.34 MPa;彈丸流場(chǎng)在彈丸噴嘴出口處的最大軸向動(dòng)壓強(qiáng)分別為2.81、3.40、3.94、4.13 和4.63 MPa。通過(guò)比較可知，在相同的噴丸壓力下水流場(chǎng)的軸向動(dòng)壓強(qiáng)小于彈丸流場(chǎng)的軸向動(dòng)壓強(qiáng)。

2.2 彈丸噴嘴長(zhǎng)度對(duì)內(nèi)流場(chǎng)的影響

圖4 和圖5 分別給出了在噴丸壓力為10 MPa，彈丸噴嘴長(zhǎng)度分別為25、35、45、55 和65 mm 時(shí)，內(nèi)部水流場(chǎng)和彈丸流場(chǎng)的軸向速度分布云圖。

由圖4 可知，其水流場(chǎng)的變化規(guī)律與圖2 基本相同，在彈丸噴嘴中，其軸向速度沿軸線(xiàn)方向呈對(duì)稱(chēng)分布，且軸線(xiàn)上速度最大。另外，由于水在彈丸噴嘴內(nèi)運(yùn)動(dòng)過(guò)程中存在沿程阻力損失，消耗了一定的能量，因此，隨著彈丸噴嘴長(zhǎng)度的增加水流場(chǎng)的軸向速度逐漸減小;在彈丸噴嘴長(zhǎng)分別為25、35、45、55 和65 mm 的情況下，水流場(chǎng)在彈丸噴嘴出口處的最大軸向速度分別為69.3、68.2、67.0、63.3 和62.8 m/s。

由圖5 可知，在彈丸噴嘴長(zhǎng)分別為25、35、45、55 和65 mm 的情況下，彈丸流場(chǎng)在彈丸噴嘴出口處的最大軸向速度分別為49.3、51.8 、54.2 、56.9 和57.4 m/s;通過(guò)比較可知，在相同的彈丸噴嘴長(zhǎng)度下水流場(chǎng)的軸向速度大于彈丸流場(chǎng)的軸向速度，但隨著彈丸噴嘴長(zhǎng)度的增加水流場(chǎng)的軸向速度逐漸減小，彈丸流場(chǎng)的軸向速度逐漸增大。

噴丸產(chǎn)生的軸向動(dòng)壓強(qiáng)隨著彈丸噴嘴長(zhǎng)度的增加而逐漸減小，在彈丸噴嘴長(zhǎng)分別為25、35、45、55 和65 mm 的情況下，水流場(chǎng)在彈丸噴嘴出口處的最大軸向動(dòng)壓強(qiáng)分別為2.31、2.11、2.04、1.97和1.89 MPa;在彈丸噴嘴長(zhǎng)度分別為25、35、45、55 和65 mm 的情況下，彈丸流場(chǎng)在彈丸噴嘴出口處的最大軸向動(dòng)壓強(qiáng)分別為4.37、4.20、4.02、3.94 和3.76 MPa。通過(guò)比較可知，在相同的彈丸噴嘴長(zhǎng)度下水流場(chǎng)的軸向動(dòng)壓強(qiáng)小于彈丸流場(chǎng)的軸向動(dòng)壓強(qiáng)。

圖4 不同彈丸噴嘴長(zhǎng)度下內(nèi)部水流場(chǎng)的軸向速度分布云圖Fig．4 Inside water flow field under different pill nozzle length axial velocity distribution nephogram

圖5 不同彈丸噴嘴長(zhǎng)度下內(nèi)部彈丸流場(chǎng)的軸向速度分布云圖Fig．5 Inside pill flow field under different pill nozzle length axial velocity distribution nephogram

3 結(jié) 論

(1)在相同的噴丸壓力下彈丸噴嘴出口處水流場(chǎng)的軸向速度大于彈丸流場(chǎng)的軸向速度，水流場(chǎng)的軸向動(dòng)壓強(qiáng)小于彈丸流場(chǎng)的軸向動(dòng)壓強(qiáng);隨著噴丸壓力的增加，水流場(chǎng)的軸向速度和軸向動(dòng)壓強(qiáng)均是逐漸增大。

(2)在相同的彈丸噴嘴長(zhǎng)度下彈丸噴嘴出口處水流場(chǎng)的軸向速度大于彈丸流場(chǎng)的軸向速度，但隨著噴嘴長(zhǎng)度的增加水流場(chǎng)的軸向速度逐漸減小，彈丸流場(chǎng)的軸向速度逐漸增大，彈丸流場(chǎng)的軸向速度逐漸趨近于水流場(chǎng)的軸向速度;水流場(chǎng)的軸向動(dòng)壓強(qiáng)小于彈丸流場(chǎng)的軸向動(dòng)壓強(qiáng)。

［1］董 星，段 雄．高壓水射流噴丸強(qiáng)化技術(shù)［J］．表面技術(shù)，2005，34(1):48 －49．

［2］呂國(guó)慶，張新民．前混合水射流強(qiáng)化技術(shù)及發(fā)展［J］．機(jī)械工程與自動(dòng)化，2010，21(3):215 －216．

［3］CHILLMAN，RAMULU，HASHISH．Water jet peening and surface preparation at 600MPa:A preliminary experimental study［J］．American Society of Mechanical Engineers，2007，129(3):485 －490．

［3］董 星，陳長(zhǎng)霞，吳斌斌，等．柱面零件噴丸強(qiáng)化殘余應(yīng)力場(chǎng)數(shù)值模擬研究［J］．黑龍江科技學(xué)院學(xué)報(bào)，2012，22(3):287 －292．

［5］AZMIR AZHARI，CHRISTIAN SCHINDLER，BO LI．Effect of waterjet peeing on aluminum alloy 5005［J］．Int J Adv Manuf Technol，2012，19(3):259 －263

［6］鈴木貞次．ウォータージェットを用いた歯車(chē)の高強(qiáng)度化法［J］．熱処理，1998，38(5):273 －276．

［7］AROLA D，ALADE A E，WEBER W．Improving fatigue strength of metals using abrasive waterjet peening［J］．Machining Science and Technology，2006，10(3):197 －218．

［8］SADASIVAM B，HIZAL A，AROLA D．Abrasive waterjet peening with elastic prestress:Effect of boundary conditions［C］//Proceedings of the ASME international mechanical engineering congress and exposition．Boston:Massachusetts，2008:672 －675．

［9］董 星，任樹(shù)棟，孫劍峰，等．后混合水射流噴丸噴頭外流連續(xù)相數(shù)值模擬［J］．黑龍江科技學(xué)院學(xué)報(bào)，2011，21(3):206 －209．

［10］BALAJI SADASIVAM，DWAYNE AROLA．An examination of abrasive waterjet peening with elastic pre-stress and the effects of boundary conditions［J］．Machining Science and Technology．2012，16(1):71－95．