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