基于流固耦合的大型礦用灑水車礦區(qū)行駛防浪涌分析

金 斌，姜永正，戴德志，李梅龍

1湘電重型裝備有限公司 湖南湘潭 411100

2湖南科技大學(xué) 湖南湘潭 411201

大 型露天礦山作業(yè)環(huán)境灰塵多，嚴(yán)重影響到礦山作業(yè)設(shè)備的使用以及工作人員的身體健康，因此，礦山需要配備一種裝載量大、噴灑效果好、工作效率高且行駛穩(wěn)定性好的大型礦山灑水車，以滿足礦山作業(yè)環(huán)境的除塵要求。

大型礦山灑水車主要由水箱總成、水路系統(tǒng)、氣路系統(tǒng)、液壓系統(tǒng)和電氣系統(tǒng)等部分組成，適用于礦山礦場道路灑水噴霧，降塵、消防滅火、輸送工業(yè)用水。灑水車工作時的安全穩(wěn)定性主要取決于水箱的設(shè)計，因此筆者將對灑水車水箱的結(jié)構(gòu)運行安全性以及穩(wěn)定性進(jìn)行分析，研究何種工況條件下對水箱結(jié)構(gòu)的強度及穩(wěn)定性影響最大。

水箱內(nèi)液體沖擊的主要研究方法有準(zhǔn)靜態(tài)方法[1]、流體動力學(xué)法[2-3]、等效機械模型法試驗[4-6]和仿真方法[7-8]。液體浪涌與灑水車在轉(zhuǎn)向、緊急制動等工況下的穩(wěn)定性一方面與車輛整體結(jié)構(gòu)布局、懸架、輪胎、質(zhì)心位置等自身設(shè)計參數(shù)有關(guān)外，另一方面也與罐內(nèi)液體的晃動情況有密切關(guān)系[9-10]。目前，防止浪涌現(xiàn)象的方法一般有添加隔板[11]或者防波板[12]，也有研究使用一種彈性膜[13]對液面晃動進(jìn)行抑制。根據(jù)對實際工作環(huán)境的考慮，以及各種方法的可行性分析，判定采取隔板的方法對浪涌現(xiàn)象進(jìn)行抑制，并開展相關(guān)仿真工作。以一款容積為 120 m3電動輪灑水車水箱為研究對象，采用 ANSYS 軟件對該水箱工作時的各種不同工況進(jìn)行仿真分析。考慮到其裝載的工作介質(zhì)是水，在行駛過程中，任何改變車輛運行狀態(tài)的因素都會導(dǎo)致箱體中的水產(chǎn)生浪涌，從而在水箱中產(chǎn)生沖擊，使水箱及其隔板受到?jīng)_擊載荷的作用。針對勻速、轉(zhuǎn)彎、啟停和過坑條件下的 8 種工況，進(jìn)行水箱結(jié)構(gòu)強度分析，研究不同工況對水箱結(jié)構(gòu)強度的影響。

1 水箱的結(jié)構(gòu)及主要參數(shù)

1.1 水箱設(shè)計

選取 SF32601 型灑水車作為搭載車輛，其基礎(chǔ)參數(shù)如表 1 所列。

表1 SF32601 型灑水車基礎(chǔ)參數(shù)Tab.1 Basic parameters of SF32601 sprinkling vehicle

1.2 水箱結(jié)構(gòu)

水箱采用隔板設(shè)計，內(nèi)部設(shè)有縱橫隔板，用以減輕浪涌對水箱壁的沖擊，隔板底部有通孔，便于水在水箱中的流動。水箱結(jié)構(gòu)模型如圖 1 所示。

圖1 水箱結(jié)構(gòu)模型Fig.1 Model of tank structure

1.3 材料參數(shù)

水箱板材材料為 Q355，其力學(xué)性能如表 2 所列。

表2 Q355 力學(xué)性能Tab.2 Mechanical property parameters of Q355

2 水箱流固耦合有限元分析

2.1 水箱典型工況

SF32601 型灑水車主要用于澳洲礦區(qū)，結(jié)合澳州重卡行駛規(guī)范手冊《LOADRESTRAINT 2018》，以澳洲礦區(qū)的實際路況進(jìn)行有限元分析，充分考慮了行駛過程中的常規(guī)勻速、轉(zhuǎn)彎、剎車等工況，各工況如表 3 所列。

表3 工況計算列表Tab.3 Calculation in various operation modes

在考慮側(cè)向離心力作用時，取轉(zhuǎn)彎行駛速度為30 km/h，轉(zhuǎn)彎半徑為 12.5 m，由此計算側(cè)向加速度為 5.51 m/s2。重力加速度設(shè)為 9.8 m/s2。

在工況 8 情況下，同時考慮了轉(zhuǎn)彎以及制動工況。考慮側(cè)向離心力作用，取轉(zhuǎn)彎行駛速度為 15km/h，轉(zhuǎn)彎半徑為 12.5 m，由此計算側(cè)向加速度為1.39 m/s2，行駛方向加速度取 -2.1 m/s2。

2.2 水箱流固耦合有限元建模

2.2.1 水箱建模

結(jié)合水箱整體尺寸，由于水箱板材主體尺寸遠(yuǎn)大于厚度方向尺寸，故在有限元軟件 ANSYS/LS-DYNA中決定采用殼單元 Shell163 對水箱板結(jié)構(gòu)進(jìn)行離散分析。水箱模型如圖 2 所示。

圖2 水箱模型Fig.2 Tank model

2.2.2 流體建模

水箱主要承受 120 m3水重力載荷，水結(jié)構(gòu)模型設(shè)定為實體 ALE 單元，考慮到水箱內(nèi)流體涉及到水和空氣 2 種流體介質(zhì)，采用多物質(zhì) ALE 算法進(jìn)行求解，單元算法采用單點 Euler/ALE 多物質(zhì)單元算法。多物質(zhì)單元是指這種單元劃分的網(wǎng)格中，允許多種物質(zhì)的流動，在同一個網(wǎng)格中，可以包含多種材料的物質(zhì)，因此，求解物質(zhì)在網(wǎng)格仿真模型中的流體材料包括水和空氣 2 種介質(zhì)，建模時將二者間處理為共面，劃分網(wǎng)格時 2 種材料間即為共節(jié)點。水和空氣均采用Null 材料模型，采用 Gruneisen 狀態(tài)方程。流體網(wǎng)格劃分如圖 3 所示。

圖3 流體網(wǎng)格劃分Fig.3 Mesh generation of fluid

2.3 各工況下邊界條件

2.3.1 行駛工況約束

在前 6 種工況下水箱約束方式相同，水箱行駛工況約束如圖 4 所示。水箱底部縱梁直接與橡膠墊之間定義接觸，舉升支座處由于安裝有液壓缸支撐，故約束其垂向位移，斗銷孔處與車架之間定義轉(zhuǎn)動副進(jìn)行約束，4 個懸架支點進(jìn)行零位移約束。

圖4 水箱行駛工況約束Fig.4 Constraints of tank in various operation modes

2.3.2 過坑工況約束

第 7 種工況即過坑工況，約束方式與前 6 種工況的區(qū)別在于：釋放前左位置的懸架約束，同時施加垂向 200 mm 位移。過坑位移曲線如圖 5 所示。

圖5 過坑位移曲線Fig.5 Variation curve of displacement in pit-crossing operation mode

3 結(jié)果與分析

8 種工況下水箱的計算結(jié)果如表 4 所列。由表 4可知，最大應(yīng)力工況為工況 1，位置在隔板開孔處，工況 1 的安全系數(shù)最小，對水箱沖擊影響最大。

表4 8 種工況下的最大應(yīng)力及分布Tab.4 Maximum stress and distribution in eight operation modes

不同工況下，水箱內(nèi)水的最大流速如圖 6 所示。在工況 7 水箱處于過坑工況，水箱內(nèi)水流速度最大為9 897 mm/s。工況 7 狀態(tài)下的水流形態(tài)及流速分布如圖 7 所示。

圖6 水箱內(nèi)水的最大流速Fig.6 Maximum flow velocity of water in tank

圖7 水流形態(tài)及流速分布Fig.7 Water flow morphology and velocity distribution

8 種工況下，水箱的整體側(cè)向位移最大出現(xiàn)在工況 7，如圖 8 所示。在工況 7 的工作條件下，水箱的側(cè)向位移最大為 93.2 mm。工況 7 的水箱側(cè)向位移分布如圖 9 所示。

圖8 水箱最大側(cè)向位移Fig.8 Maximum lateral displacement of tank

圖9 水箱側(cè)向位移分布Fig.9 Lateral displacement distribution of tank

由圖 6 可知，過坑工況水的流速最大，達(dá) 9.9 m/s，但是其流速是由于水箱整體向下做剛體位移導(dǎo)致的，因此對結(jié)構(gòu)產(chǎn)生的應(yīng)力并非最大。此外，在轉(zhuǎn)彎時水的流速也較大，此時水向側(cè)向移動，對縱隔板及側(cè)板產(chǎn)生沖擊，從而導(dǎo)致縱隔板與底板的焊縫處出現(xiàn)最大應(yīng)力。在制動及啟動工況中，水對橫隔板及前后板會產(chǎn)生沖擊，但是由于制動及啟動的加速度小于轉(zhuǎn)彎加速度，所以制動及啟動工況的應(yīng)力小于轉(zhuǎn)彎工況。

8 種工況下，水箱隔板、蓋板以及底板的最大應(yīng)力如圖 10 所示。隔板上的應(yīng)力最大出現(xiàn)在工況 1，為275 MPa，也是水箱整體的最大應(yīng)力；蓋板最大應(yīng)力出現(xiàn)在工況 7，為 264 MPa；底板的最大應(yīng)力也出現(xiàn)在工況 7，為 164 MPa。

圖10 水箱最大應(yīng)力Fig.10 Maximum stress of tank

通過對比各工況對水箱側(cè)向位移、流速、應(yīng)力的影響，發(fā)現(xiàn)轉(zhuǎn)彎工況和過坑工況對各項參數(shù)影響顯著；對比工況 1 與工況 8，轉(zhuǎn)彎時制動減速可明顯減小水箱內(nèi)隔板的應(yīng)力集中；底板在各工況下，應(yīng)力相對較小。

4 結(jié)語

在 8 種典型行駛工況下，通過灑水車水箱內(nèi)部浪涌對結(jié)構(gòu)影響的對比分析可知：

（1) 轉(zhuǎn)彎工況下，水箱內(nèi)水對箱體的影響尤為明顯，水箱隔板的應(yīng)力很大，且隨著充液比的提高，水箱內(nèi)部隔板的應(yīng)力也隨之增大。

（2) 過坑工況下，水箱的最大側(cè)向位移較其他工況下會有明顯增大，且該工況下水箱的最大應(yīng)力為265 MPa。

（3) 水箱在灑水車做轉(zhuǎn)向運動時，其最大應(yīng)力分布在水箱內(nèi)隔板處，制動、啟動、過坑工況時，其最大應(yīng)力主要分布在水箱蓋板上。

（4) 為保障灑水車的正常工作，水箱及灑水車設(shè)計時，應(yīng)充分考慮轉(zhuǎn)向與過坑工況對水箱結(jié)構(gòu)以及車輛穩(wěn)定性的影響，并且車輛轉(zhuǎn)向時，應(yīng)適當(dāng)降低車速。