針對前混合式磨料射流磨料加速過程運(yùn)動復(fù)雜、實驗研究困難及有限元處理超大變形存在網(wǎng)格畸變等問 題，基于光滑粒子（ＳＰＨ）耦合有限元（ＦＥＭ）方法模擬前混合式磨料射流噴嘴不同階段磨料粒子加速特征及磨料射流破 碎靶體全過程。其中水介質(zhì)用 ＳＰＨ建模，磨料粒子、噴嘴、靶體等用 ＦＥＭ建模。揭示磨料粒子群在噴嘴中的運(yùn)動軌跡及 噴嘴結(jié)構(gòu)對磨料粒子加速影響規(guī)律。研究表明，磨料粒子進(jìn)入噴嘴收斂段之前已基本達(dá)到與水介質(zhì)相同速度，進(jìn)入收斂 段后因與水介質(zhì)存在速度差使其獲得加速，但與水介質(zhì)速度差逐漸增大；進(jìn)入直線段后水介質(zhì)與磨
磨料水射流作為新型特種加工技術(shù)，將磨料與高速流水或高壓水互相混合形成液固兩相介質(zhì)射流［１－３］。據(jù)射流混合方式不同，磨料水射流可分為前、后混合式磨料水射流［４］。前者因能使磨料與高速流水充分紊動混合，較后者切割、破碎效率更高。磨料粒子的加速效率為決定磨料射流破碎結(jié)構(gòu)體效果的主要原因，因此研究前混合式磨料射流中磨料粒子加速過程，對提高水介質(zhì)對磨料的加速效率，優(yōu)化噴嘴結(jié)構(gòu)具有重要意義。 前混合式磨料水射流磨料在噴嘴內(nèi)運(yùn)動復(fù)雜，實驗研究需在尺寸較小的噴嘴內(nèi)測量超高速射流運(yùn)動、磨料粒子運(yùn)動軌跡及兩者間混合時相互作用現(xiàn)象極其復(fù)雜、困難，但隨計算機(jī)技術(shù)、理論發(fā)展，可用數(shù)值模擬計算手段對該問題分析研究。Ｋｕｍａｒ等［５－８］用非線性有限元 方 法 模 擬 磨 料 粒 子 破 碎 結(jié) 構(gòu) 體 過 程；Ｗａｎｇ等［９－１１］用 ＳＰＨ算法模擬水射流破巖過程；司鵠等［１２］用ＡＬＥ流 －固耦合罰函數(shù)算法模擬磨料水射流破煤巖過程。以上研究僅局限于噴嘴外射流模擬，側(cè)重點是靶體受作用后的力學(xué)分析，且對磨料、水均賦予相同初速度，兩者間無直接能量傳遞。Ｎgｕｙｅｎ等［１３－１７］對兩相流在噴嘴內(nèi)、外流場數(shù)值模擬，獲得較好模擬結(jié)果，且在一定程度上可獲得磨料粒子的運(yùn)動情況。王建明等［1８］基于 ＳＰＨ與 ＦＥＭ耦合法模擬磨料水射流中單磨粒加速過程，但僅針對后混合磨料射流中單顆粒磨料在水刀直線加速段的加速過程進(jìn)行模擬。 本文采用 ＳＰＨ與 ＦＥＭ耦合算法，將前混合磨料水射流的形成過程考慮成柱塞以一定運(yùn)動速度作用于水，使水以一定流量向前運(yùn)動，帶動事先隨機(jī)分布于水中的磨料，通過收縮噴嘴形成高速磨料水射流。重點研究磨料粒子在噴嘴內(nèi)部及距出口一定范圍內(nèi)的加速過程，分析噴嘴結(jié)構(gòu)對磨料粒子加速效果影響規(guī)律。１ ＳＰＨ方法相關(guān)理論 磨料水射流形成過程中高速流動的水涉及大變形問題，采用傳統(tǒng)拉格朗日法模擬會出現(xiàn)網(wǎng)格畸變導(dǎo)致計算終止。ＳＰＨ法不用單元而用固定質(zhì)量的可動點，無網(wǎng)格畸變。磨料粒子用 ＦＥＭ模擬，利用有限元計算精度高的特點能準(zhǔn)確獲得加速過程中運(yùn)動特性。因此，ＳＰＨ與ＦＥＭ耦合算法既可較好模擬磨料水射流的形成過程，又能準(zhǔn)確獲得磨料粒子在噴嘴不同階段的加速特性。１.１ ＳＰＨ基本理論 ＳＰＨ方法基礎(chǔ)為插值理論［１９－２０］。在 ＳＰＨ中任一宏觀變量均能方便借助一組無序點值表示成積分插值計算獲得。核近似函數(shù)為
１．２ ＳＰＨ-ＦＥＭ 耦合原理
ＳＰＨ粒子在 ＬＳ-ＤＹＮＡ中被視為特殊節(jié)點單元，控制參數(shù)為節(jié)點編號、質(zhì)量及空間位置。ＳＰＨ粒子與ＦＥＭ耦合通過罰函數(shù)方式將質(zhì)點力作用于有限元單元表面，因此 ＳＰＨＦＥＭ耦合采用 ＮＯＤＥ＿ＴＯ＿ＳＵＲＦＡＣＥ點面接觸方式。本文 ＳＰＨ與 ＦＥＭ耦合處理中，將 ＳＰＨ粒子定義為從節(jié)點，與 ＳＰＨ粒子接觸界面的有限元單元表面定義為主面。接觸耦合算法見圖 １［２１］。
２ 計算模型描述
為對問題進(jìn)行一定簡化，本文設(shè)：① 磨料水射流形成過程僅涉及水、磨料、噴嘴、推動水前進(jìn)的柱塞及磨料水射流沖蝕靶體 ５種物質(zhì)；② 磨料粒子為等直徑球體并事先隨機(jī)分布于水中；③ 高壓水由一定速度前進(jìn)的柱塞推動水形成；④靶體為連續(xù)介質(zhì)體。
２．１ 水、磨料模型
水采用 ＳＰＨ算法建模，在 ＬＳ-ＤＹＮＡ中用 ＮＵＬＬ材料模型，并對水材料模型賦予 Ｍｉｅ-Ｇｒｕｅｉｓｅｎ狀態(tài)方程；磨料粒子為陶粒，采用 ｌａｇｒａｎｇｅ有限元建模，并用剛體材料模型，直徑均為 ０.９ｍｍ（對應(yīng)２０目）。建模過程中為盡可能接近實際工況，建立密集的 ＳＰＨ水粒子束，將磨料粒子分布于水粒子中；為防止水粒子對磨料粒子初始穿透，將每顆磨料粒子位置的ＳＰＨ水粒子逐一刪除，以“挖”出磨料粒子的存在空間。水與磨料模型參數(shù)分別見表 １、表 ２，其中 ρ０為材料密度；Ｅ為彈性模量；υ為泊松比；?為直徑。Ｍｉｅ-Ｇｒｕｅｉｓｅｎ狀態(tài)方程為
式中：Ｅ為單位體積內(nèi)能；Ｃ為 ｖｓ－ｖｐ曲線截距；ｓ１，ｓ２，ｓ３為 ｖｓ－ｖｐ曲線斜率系數(shù)；γ０ 為 Ｇｒｕｅｉｓｅｎ常數(shù)；ａ為一介體積修正量。
２．２ 噴嘴、柱塞模型
將前混合磨料水射流形成過程考慮為柱塞以一定速度作用于水，使水以一定流量向前運(yùn)動，帶動事先隨機(jī)分布于水中的磨料，通過收縮噴嘴形成高速磨料水射流。噴嘴、柱塞均采用 ｌａｇｒａｎｇｅ有限元建模，采用剛體材料模型，其結(jié)構(gòu)示意圖見圖 ２。
２．３ 靶體模型
靶體為鋼材料采用 ｌａｇｒａｎｇｅ有限元建模，選用ＰＬＡＳＴＩＣ＿ＫＩＮＥＭＡＴＩＣ材料模型，該材料模型按照單元的應(yīng)變值是否超過材料失效應(yīng)變值確定是否刪除單元，可在一定程度上模擬靶體受磨料水射流作用后沖蝕孔的發(fā)展過程。材料參數(shù)見表 ３，其中 ρ０ 為材料密度；Ｅ為彈性模量；υ為泊松比；σ為失效應(yīng)力；ε為失效應(yīng)變。
２．４ 耦合處理
在 ＳＰＨ-ＦＥＭ耦合算法中涉及 ＳＰＨ粒子及 ｌａｇｒａｎｇｅ有限元，兩種算法間相互作用需進(jìn)行耦合處理。此模擬中水采用ＳＰＨ算法建模，而磨料、噴嘴、柱塞及靶
體則采用ｌａｇｒａｎｇｅ有限元法建模。為實現(xiàn)在柱塞推動下水?dāng)y帶磨料粒子通過噴嘴形成高速磨料水射流作用于靶體上目的，水與柱塞、水與磨料粒子、水與噴嘴、水與靶體
間均通過 ＮＯＤＥ＿ＴＯ＿ＳＵＲＦＡＣＥ關(guān)鍵字實現(xiàn)耦合；該模擬過程涉及到磨料粒子與噴嘴及磨料粒子與靶體的碰撞，考慮磨料粒子單元相對其余兩者尺寸較小，也通過ＮＯＤＥ＿ＴＯ＿ＳＵＲＦＡＣＥ關(guān)鍵字實現(xiàn)相互作用；磨料粒子之間碰撞通過設(shè)定 ＳＩＮＧＬＥ＿ＳＵＲＦＡＣＥ關(guān)鍵字實現(xiàn)。
２．５ 幾何模型及邊界條件
據(jù)以上建模方法建立幾何模型見圖 ３。圖 ３（ｂ）為（ａ）中管路水柱段內(nèi)部水粒子與磨料粒子分布。限制柱塞自由度，使其只在 Ｚ方向移動，柱塞直徑 １４ｍｍ，厚度 １ｍｍ；限制噴嘴在各方向移動，其中管路水柱段直徑 １４ｍｍ，長度 １００ｍｍ，噴嘴收斂段長度 ２０ｍｍ，噴嘴直線段直徑 ３ｍｍ，長度 １０ｍｍ；靶體為 ２０ｍｍ×２０ｍｍ×１５ｍｍ長方體，對底面、外圍 ４面采用 ＮＯＮ＿ＲＥ-ＦＬＥＣＴＩＮＧ非反射約束模擬無限邊界，同時約束地面節(jié)點在各自由度的位移。水模型含 １３７００個 ＳＰＨ粒子，７２個磨料粒子，共計 １８４３２個六面體單元；噴嘴模型含 ３３８０個六面體單元；柱塞模型含 ３８４個六面體單元；靶體模型含 ６０００個六面體單元。
３ 計算結(jié)果分析
據(jù)建模方法進(jìn)行兩次模擬計算，即側(cè)重于磨料水射流形成后對靶體作用及將模型中靶體刪除并著重分析研究磨料粒子在噴嘴不同階段加速過程。
流量 ７４Ｌ／ｍｉｎ的高壓水?dāng)y帶磨料粒子群經(jīng)噴嘴形成高速磨料水射流作用于靶體后效果圖及沖蝕坑截面形狀演化圖見圖 ４、圖 ５，該流量高壓水?dāng)y帶磨料粒子群在噴嘴收斂段、直線段及出口段運(yùn)動全過程見圖 ６。
３．１ 靶體沖孔形態(tài)分析
由圖 ５看出，在磨料水射流沖蝕作用中心區(qū)域的靶體單元迅速破壞失效，形成初始孔徑，隨沖蝕過程延續(xù)一部分粒子繼續(xù)加深沖蝕深度，同時與孔底撞擊后向周圍發(fā)散飛濺的另部分粒子對初始孔壁沖蝕并擴(kuò)大形成“Ｖ”形剖面，并使沖蝕坑不斷加深。此與文獻(xiàn)［１２
－２２］結(jié)果吻合較好，表明本文建模、分析方法可行。
３．２ 磨料粒子加速過程分析
由于數(shù)值模擬中建立的磨料粒子非單獨(dú)一顆，在磨料粒子的加速過程中涉及到磨料粒子間碰撞，為盡可能獲得完整加速過程，在磨料粒子群即將進(jìn)入噴嘴收斂段時提取其前端靠近噴嘴軸心線的一顆磨料粒子，選水介質(zhì)部分中與提取的磨料粒子處于 Ｚ軸方向同一斷面上一顆 ＳＰＨ水粒子，繪制噴嘴軸線上速度 －時間曲線，分別按磨料粒子、水粒子運(yùn)動歷程將噴嘴不同階段標(biāo)注于速度 －時間曲線上，見圖 ７。由圖 ７看出，磨料粒子進(jìn)入噴嘴收斂段前已基本達(dá)到與水相同的運(yùn)動速度，進(jìn)入收斂段后在流量一定情況下，由于過流斷面積變小，水流速度急劇增大使磨料粒子與水之間存在速度滑移，該速度差的存在使磨料粒子獲得加速，但由于磨料粒子質(zhì)量大，慣性大，故磨料粒子提速相對緩慢，導(dǎo)致在收斂段中其與水速度差逐漸增大；進(jìn)入直線段后，由于收斂段積聚的壓能部分釋放，水流速度繼續(xù)增大，但增量有限并在直線段末端趨于穩(wěn)定，磨料粒子由于在收斂段導(dǎo)致速度滯后較大，使其在進(jìn)入直線段后與水存在巨大速度差，直線段過流斷面較小，使磨料粒子在直線段劇烈加速，并隨時間延續(xù)逐漸趨向水流的速度；在離開噴嘴后的短距離范圍內(nèi)，不考慮噴嘴外空氣阻力情況，水流有小幅度速度增加并迅速穩(wěn)定，主要因水介質(zhì)在噴嘴中由于瞬間高壓存在微量的體積壓縮，離開噴嘴后能量釋放導(dǎo)致速度微量增加，而磨料粒子在噴嘴外核心段射流的繼續(xù)作用下仍有速度增加，但最終因磨料粒子運(yùn)動逐漸脫離核心段射流部分導(dǎo)致速度趨于穩(wěn)定。數(shù)值模擬結(jié)果與文獻(xiàn)［１５，２３－２７］研究成果吻合較好。
３．３ 磨料粒子群運(yùn)動軌跡分析
為獲得磨料粒子在噴嘴中的運(yùn)動形態(tài)，選５顆粒子追蹤其在噴嘴中的運(yùn)動軌跡，見圖８。由圖８看出，磨料粒子運(yùn)動軌跡表明其在收斂段的運(yùn)動較劇烈，因過流斷面變窄，其與噴嘴壁面及磨料粒子之間均存在碰撞；進(jìn)入直線段后因滯留時間較短，磨料粒子彼此之間碰撞機(jī)會變小，運(yùn) 動 相 對 收 斂 段平緩。
３．４ 噴嘴結(jié)構(gòu)對磨料
影響磨料粒子加速因素眾多，現(xiàn)有研究結(jié)果表明，噴嘴收斂段長度、直線段長度對磨料加速性能有重要影響。本文通過控制變量法，分別在只改變收斂段長度或只改變直線段長度下研究噴嘴結(jié)構(gòu)對磨料粒子加速影響規(guī)律。
按本文模擬方法，在流量一定即 ７４Ｌ／ｍｉｎ，噴嘴進(jìn)出口直徑一定即１４ｍｍ、３ｍｍ，噴嘴直線段長度１０ｍｍ下改變收斂段長度，選４顆磨料粒子獲得噴嘴出口（直線段末端）平均速度，見圖９。由圖 ９看出，噴嘴收斂段加長下由于流量一定，噴嘴進(jìn)出口速度一定，水流經(jīng)噴嘴收斂段的平均速度變小，由此可以推斷水流經(jīng)收斂段加速度變小，從而對磨料的加速較弱，但收斂段加長會延長磨料粒子的加速歷程，使其有相對充裕的加速時間。隨收斂段加長，磨料粒子速度呈增大趨勢，收斂段長度由１０ｍｍ增加到１５ｍｍ，速度增加１ｍ／ｓ；收斂段長度由１５ｍｍ到３０ｍｍ每增加５ｍｍ，速度增加０．５ｍ／ｓ。由此，收斂段長度變化２０ｍｍ 速度增加 ２．５ｍ／ｓ，增加效果有限。
流量一定即７４Ｌ／ｍｉｎ，噴嘴進(jìn)出口直徑一定即１４ｍｍ、３ｍｍ，噴嘴收斂段長度 ２０ｍｍ情況下，改變直線段長度，選４顆磨料粒子獲得噴嘴出口（直線段末端）平均速度，見圖１０。由于收斂段中磨料粒子已與水存在速度差，且直線段中水仍有加速，故直線段加長會致磨料粒子加速空間更長。由圖１０看出，隨直線段加長磨料粒子在噴嘴出口速度增大，但隨直線段長度逐漸增大，磨料粒子速度增加幅度逐漸變小，直線段長度由５ｍｍ到１５ｍｍ，速度增加較大，平均每增加５ｍｍ速度增大２．7５ｍ／ｓ；直線段長度由１５ｍｍ到２０ｍｍ，速度增大１．５ｍ／ｓ；直線段長度由２０ｍｍ到２５ｍｍ，速度增大０．５ｍ／ｓ。直線段長度變化 ２０ｍｍ速度增加７．５ｍ／ｓ，增加幅度較收斂段有所增大，但當(dāng)磨料粒子速度增加到與水流速度接近后，其速度增長幅度隨直線段延長變小，且在實際應(yīng)用中直線段過長會使水與磨料粒子的沿程阻力增大導(dǎo)致能耗增加，故在滿足工況設(shè)計前提下，可適當(dāng)延長噴嘴直線段長度，但不能過多延長。
４ 結(jié) 論
（１）本文基于 ＳＰＨＦＥＭ耦合算法模擬磨料粒子在噴嘴不同階段的加速全過程，既接近實際工況又使高流速水因大變形導(dǎo)致模擬計算困難問題得以解決。
（２）通過分析前混合式磨料射流磨料粒子不同階段的加速特征知，磨料粒子進(jìn)入噴嘴收斂段前已基本達(dá)到與水介質(zhì)相同的運(yùn)動速度，進(jìn)入后因與水介質(zhì)存在速度差獲得加速，因磨料粒子質(zhì)量大，與水介質(zhì)速度差逐漸加大；進(jìn)入直線段水介質(zhì)與磨料粒子一直加速，且水介質(zhì)速度在末端趨于穩(wěn)定；在噴嘴短距離內(nèi)水介質(zhì)在噴嘴中的壓縮獲得釋放，速度微量增加且穩(wěn)定較快，磨料粒子在噴嘴外核心段射流作用下繼續(xù)加速，最終趨于穩(wěn)定。
（３）獲得噴嘴結(jié)構(gòu)對加速效果影響規(guī)律，即流量一定，磨料粒子速度隨噴嘴收斂段延長而增加。收斂段長度變化 ２０ｍｍ速度增加 ２．５ｍ／ｓ，增加效果有限；隨直線段的延長而增加，直線段長度變化 ２０ｍｍ速度增加 ７．５ｍ／ｓ，速度增加相對顯著。
（４）本文所提數(shù)值模擬方法可用于前混合磨料水射流中磨料粒子加速特性及水射流撞擊結(jié)構(gòu)體的研究。對不同特征的磨料粒子（如密度、尺寸、形狀等）、不同結(jié)構(gòu)噴嘴尚待進(jìn)一步研究。
首頁
產(chǎn)品中心
電話咨詢
售后服務(wù)
售前服務(wù)
關(guān)注公眾號