任云鵬 李安帥 李旭東 宋 芳

(1.沈陽建筑大學機械工程學院,遼寧 沈陽 110168;2.山東黃金礦業(yè)(萊州)有限公司三山島金礦,山東 萊州 261442)

高壓輥磨機是一種新型破磨粉碎設備,具有處理能力強、單位破碎能耗低、設備工作效率高等特征[1]。該設備最先應用于水泥領域,隨著技術的不斷發(fā)展,已經廣泛應用于金屬礦山領域[2-3]。高壓輥磨機的相關研究涉及機械運動學、動力學、離散元等諸多方面,當前研發(fā)設計中仍存在現(xiàn)有模型不準確、破碎過程中重要參數無法定量表示等問題。

虛擬仿真技術是改進粉碎設備設計和操作的一項重要方法,它能夠比較不同參數設置而不需要物理原型或生產損失的相關成本,其仿真結果的可用性取決于模型的保真度。BILGILI 等[4]用嚴格的數學模型研究了多個接觸顆粒的粉碎,得出了粉碎過程矩陣,但無法考慮到設備參數、顆粒材料性質對粉碎效果的影響。TORRES 等[5]以礦石性質、高壓輥磨機工作參數、粉碎工藝為基礎,研究了一系列高壓輥磨機粉碎的數學模型,包括處理能力模型、能耗模型,但在預測準確度上有所不足,無法應用到工業(yè)領域。DANIEL 等[6]通過大量實驗開展了關于高壓輥磨機數學模型及應用的專項研究,驗證和擴展了高壓輥磨機的產品粒度分布、處理量以及能量利用率的數學模型。KUMAR[7]使用EDEM 模擬了實驗室規(guī)格的高壓輥磨機,研究了顆粒破碎過程,以及工藝參數和操作參數對處理量和能耗的影響。祖大磊等[8]和何劍偉等[9]根據高壓輥磨機破碎過程中的工藝參數、礦石單軸抗壓強度、給料粒度及粉碎產品粒度分布對比能耗的影響,提出了利用高壓輥磨機比能耗進行計算評估的DUCS 數學模型。鮑諾[10]利用EDEM 軟件研究了高壓輥磨機工作過程中顆粒運動狀態(tài)與工作參數之間的關系,得到了不同輥速和粒徑下顆粒所受壓力隨時間變化的規(guī)律。張樂[11]利用EDEM 軟件模擬了高壓輥磨機破碎礦石過程,通過正交試驗方法和單因素試驗方法研究壓輥轉速、壓輥間隙以及壓輥與礦石間的靜摩擦因數對礦石破碎效果的影響,但單一的DEM 仿真模型忽視了設備的動態(tài)響應,無法真實模擬實際工作狀態(tài)。在高壓輥磨機中,壓輥間隙不是一個直接控制的變量,是所施加的壓力和顆粒床之間相互作用的結果。以上基于理論的數學模型和基于離散元法的仿真模型由于忽略了顆粒對高壓輥磨機運行的影響,無法描述輥隙的動態(tài)變化對仿真結果的影響,導致在預測處理量、單位能耗和破碎比時均存在較大誤差。

針對高壓輥磨機現(xiàn)有模型存在誤差和無法描述設備的動態(tài)輥隙隨顆粒載荷變化規(guī)律的問題,為了減小誤差,降低企業(yè)研發(fā)設計成本,在離散元法基礎上,提出采用離散元(DEM)和多體動力學(MBD)聯(lián)合仿真技術對高壓輥磨機進行仿真分析,通過對比數學模型、單一離散元仿真模型以及實際工作數據,證實聯(lián)合仿真模型有著較為理想的結果。

1 高壓輥磨機數學模型

1.1 動態(tài)輥隙模型

動態(tài)輥隙模型是在動輥動力學模型基礎上推導而來的,高壓輥磨機主要由兩個壓輥和液壓系統(tǒng)組成,液壓系統(tǒng)部分蓄能器中的壓縮氮氣充當液壓系統(tǒng)的彈簧。在破碎過程中,動輥和定輥的旋轉速度是恒定的,但動輥會因受到物料的擠壓力和液壓系統(tǒng)提供的力往復平移。動輥的橫向移動可以由彈簧阻尼系統(tǒng)來描述,如圖1所示。

圖1 高壓輥磨機動態(tài)輥隙模型Fig.1 Dynamic roll gap model of high-pressure grinding rolls

高壓輥磨機的輥隙動態(tài)變化可由下式表示:

1.2 處理量模型

高壓輥磨機的處理量數學模型可以基于設備的設計參數和物料的物理性質進行推導。最常見的處理量數學模型是基于高壓輥磨機的輥隙寬度、壓輥直徑、壓輥旋轉速度、壓輥面長度、料餅密度等參數來進行建模的。圖2 是高壓輥磨機的粉碎示意圖。

圖2 高壓輥磨機粉碎示意Fig.2 High-pressure grinding rolls crushing diagram

根據質量守恒定律,物料在擠壓區(qū)域內,質量不發(fā)生改變,并且假設輥面與物料之間不存在相對滑動,則高壓輥磨機的單位處理量G可以表示為:

式中,G(θ)為任意圓心角θ下高壓輥磨機的單位處理量,t/s;ρ(θ)為任意圓心角θ下物料密度,t/m3;S(θ)為任意圓心角θ下的壓輥間隙,m;L為輥寬,m;v為壓輥線速度,m/s。

在實際應用中,高壓輥磨機的處理量一般指出料口位置處,即θ=0 時,高壓輥磨機的理論處理量G為[12]:

式中,δ為料餅密度,一般為物料密度的0.85[13],t/m3;s為出料口間隙,m。

1.3 單位能耗模型

圖3 是高壓輥磨機壓輥受力示意圖。

圖3 高壓輥磨機壓輥受力示意Fig.3 Force diagram of pressure roller of high-pressure grinding rolls

如圖所示,陰影區(qū)域是工作壓力的影響區(qū)域,在該區(qū)域內壓輥對物料的擠壓力F可以通過壓輥的工作壓力求得。由于高壓輥磨機的工作特點,物料只在壓縮區(qū)受到壓輥的擠壓力,故相應的受力區(qū)域為壓輥的上半部分,因而壓輥的擠壓力F為:

式中,F為壓輥對物料的擠壓力,kN;P為壓輥的工作壓力,MPa;D為輥徑,m。

則垂直方向上的擠壓力F引起的轉矩T為:

式中,T為轉矩,kN;β為發(fā)生層壓粉碎角度的一半。

高壓輥磨機具有2 個壓輥,那么總功率是壓輥垂直方向上擠壓力F和壓輥線速度v乘積的2 倍,故功率P為:

單位能耗為總功率P與處理量G的比值,故高壓輥磨機的理論單位能耗W為[14]:

式中,W為理論單位能耗,kW·h/t;G為處理量,t/h。

1.4 破碎比模型

平均破碎比是破碎機的一個重要性能參數,用來衡量高壓輥磨機對物料的破碎效果。平均破碎比i通常定義為物料經過高壓輥磨機破碎后,進料的平均粒徑與出料的平均粒徑之比,該種方法求得的破碎比能較為真實地反映破碎程度。將顆粒劃分為體積或質量相等的兩部分的粒徑d50亦被稱為中位徑或中值粒徑,常用來表征顆粒樣品的平均粒徑[15]。高壓輥磨機出料的平均粒徑可以用d50表示,d50可通過產品粒徑分布函數羅遜—拉姆勒(Rosin-Rammler)求得。該函數是一種經驗公式,可以用來描述顆粒大小的分布情況。該函數的形式如下[16]:

式中,?(d)為累計粒度產率分布函數;d為顆粒粒徑,mm;c,n,dmax是根據經驗數據確定的參數,其中dmax為破碎產品中最大粒徑,mm。

則高壓輥磨機的理論平均破碎比i為:

式中,D50和d50分別表示破碎前后累計粒度分布百分數達到50%時所對應的粒徑。

2 DEM-MBD 聯(lián)合仿真模型建立

考慮到高壓輥磨機是一個復雜的多體動力學系統(tǒng),其運行狀態(tài)會受到多個因素的影響,當物料與壓輥之間的擠壓作用力導致物料材料性質發(fā)生變化時,如粉碎過程中的物料磨損和破碎,則必須考慮設備的動態(tài)響應,如輥隙的動態(tài)變化,因為它會影響使用DEM 進行模擬仿真的效果。對于單一DEM 仿真模型,只能考慮壓輥固定間隙,無法反映輥隙的動態(tài)變化對系統(tǒng)的影響。而輥隙的變化又是影響高壓輥磨機處理量、單位能耗、破碎比的關鍵因素,所以現(xiàn)有的DEM 仿真模型在預測高壓輥磨機的性能上具有較大誤差。高壓輥磨機的MBD 仿真需要對壓輥施加瞬態(tài)載荷,以便描述輥隙隨載荷變化而變化的特性。DEM-MBD 聯(lián)合仿真利用兩軟件之間的數據交互功能,將物料破碎過程中引起的設備實時動態(tài)響應又反饋給破碎過程。DEM-MBD 聯(lián)合仿真模型引入了輥隙的動態(tài)變化,可以更真實、更準確地模擬高壓輥磨機的運行狀態(tài)和性能,從而提高仿真結果的可靠性和有效性,故本文提出采用DEM-MBD 聯(lián)合仿真方法預測高壓輥磨機性能指標。

2.1 聯(lián)合仿真方法

高壓輥磨機聯(lián)合仿真的流程圖見圖4。當仿真開始后,MBD 部分將定義完材料參數和運動參數的動力學模型傳遞給DEM 部分,DEM 部分顆粒工廠產生顆粒,然后由EDEM 軟件計算確定每個顆粒的位置和速度,以及顆粒之間和顆粒與幾何體之間的作用力。DEM 部分結束后,EDEM 將幾何體上的等效粒子力和力矩再反饋給MBD 部分,MBD 部分由Motionview 軟件實現(xiàn)求解MBD 的控制方程,并計算幾何體的位置和速度。

圖4 DEM-MBD 聯(lián)合仿真流程圖Fig.4 DEM-MBD co-simulation flow chart

2.2 多體動力學模型

高壓輥磨機幾何機構的精確構建是聯(lián)合仿真的關鍵,以KHD375?1400-1000 型高壓輥磨機為研究對象,利用SolidWorks 軟件構建幾何機構模型,包括定輥、動輥、機架、液壓缸、氮氣蓄能器、喂料裝置六大部分,如圖5所示。

圖5 高壓輥磨機幾何模型Fig.5 Geometry model of high-pressure grinding rolls

2.2.1 約束關系構建

將構建好的KHD375?1400-1000 型高壓輥磨機幾何模型導入Motionview,此時模型之間只存在空間位置關系而不存在運動學上的約束關系。因此為保證虛擬高壓輥磨機的準確運動,首先應根據各零部件的實際運動關系構建模型之間的約束關系。其約束關系如表1所示。

表1 約束關系Table 1 Binding relationships

高壓輥磨機零部件的幾何關系和約束關系構建完后,需要進一步施加作用在壓輥上的旋轉運動,從而使高壓輥磨機能夠實現(xiàn)破碎作業(yè)。

2.2.2 液壓系統(tǒng)參數理論計算

高壓輥磨機中的液壓系統(tǒng)由3 個主要部件組成:液壓缸、氮氣蓄能器和液壓油。蓄能器中的壓縮氮氣充當液壓系統(tǒng)的彈簧,其剛度可通過初始氮氣壓力進行調整,初始油壓和氮氣壓力的組合有效地決定了彈簧的非線性剛度響應。圖6 為液壓系統(tǒng)力學模型圖,當力從f1增加到f2時,它會引起壓輥位移Δs,該位移被傳遞到橫截面積為ap的活塞上,使氮氣蓄能器中的壓力從p1增加到p2。

圖6 高壓輥磨機液壓系統(tǒng)力學模型Fig.6 Mechanical model of hydraulic system of high-pressure grinding rolls

液壓系統(tǒng)部件參數與彈簧剛度k和阻尼系數c有關,可通過以下公式估算[17]:

式中,n是多變指數;p0是初始氮氣壓力,MPa;V0是初始氮氣體積,m3;p1是初始液壓,MPa;ap是液壓缸的橫截面積,m2;np是氣缸數,個;g 是重力加速度,m/s2。該模型考慮了液壓系統(tǒng)的主要部件,但未考慮其他部件,如閥門、油管、油黏度等。

根據實際高壓輥磨機運行參數,依據上述公式計算出液壓系統(tǒng)的彈簧剛度k和阻尼系數c,通過在Motionview 中添加彈簧系統(tǒng)部件并設置彈簧剛度k和阻尼系數c,實現(xiàn)液壓系統(tǒng)的代替。

2.3 離散元模型

2.3.1 顆粒破碎模型的選用

目前EDEM 中有2 種用于模擬顆粒破碎的模型:Bonding 鍵模型和Tavares 破碎能模型。Bonding鍵模型是將小顆粒通過Bonding 鍵粘結成大顆粒,以鍵的斷裂來模擬顆粒的破碎。Bonding 鍵模型在顆粒形狀和大小相對均勻的情況下表現(xiàn)較好,但在顆粒形狀和大小差異較大時,可能會出現(xiàn)較大誤差,并且計算量比較大,破碎后無法對顆粒粒徑進行統(tǒng)計。Tavares 破碎能模型是一種基于經驗的破碎模型,當沖擊能大于破碎能時,顆粒發(fā)生破碎,沒有破碎的顆粒也會受到損傷,產生更小的破碎能。該模型基于經典的Kick、Rittinger 和Bond 破碎理論,考慮了影響顆粒內部的應力分布和破碎的各種因素,如顆粒的強度、形狀、大小、載荷等[18]。Tavares 破碎模型適用于各種形狀和大小的顆粒,且在顆粒破碎前后的形狀不會發(fā)生太大變化的情況下表現(xiàn)較好,并且計算速度快,可對破碎后顆粒進行粒徑分布、獲取破碎比等后處理,能直觀展示破碎效果。

高壓輥磨機仿真破碎過程顆粒數量巨大,且需要分析顆粒破碎后的破碎比、處理量和單位能耗等。相比Bonding 鍵破碎模型,Tavares 破碎模型能夠方便對顆粒粒徑進行統(tǒng)計,并且能同時滿足其他的仿真需求。Tavares 破碎模型原理如圖7所示。

圖7 Tavares 破碎原理Fig.7 Tavares crushing principle

在Tavares 破碎模型中,每個顆粒都有特定的破碎能量,根據其大小、平均值和標準偏差進行分配。這個能量根據公式[19-20]所描述的分布而變化。

式中,E為破碎能量分布,J/kg;Emax為能量分布的上限值,J/kg;E50為能量分布的中值,J/kg;σE為標準差。

破碎能量的中值由以下公式給出:

式中,E∞為極限破碎能,J/kg;φ是根據試驗數據調整的參數;KP和Kst分別為顆粒和鋼的硬度,GPa;dp為顆粒的大小,mm;d0為礦石中特征顆粒的大小,mm。

當顆粒受到擠壓而未發(fā)生破碎時,它將受到損傷,之后顆粒會產生新的破碎能E′f,新的破碎能會低于其原破碎能Ef:

式中,Ef為顆粒的破碎能,J;eEk為有效碰撞能量,J;D為損傷系數;γ為損傷累積系數。

Tavares 破碎模型主要是基于參數t10,t10表示產品中小于母粒徑1/10 的碎片比例,由公式描述為:

式中,A和b是根據試驗數據調整的參數。t10值越大,新產生的顆粒就越細。

根據Tavares 團隊研究成果,參考RODRIGUEZ等[21]校核的Tavares 破碎模型參數來定義物料的破碎參數,具體破碎參數如表2所示。

表2 Tavares 破碎模型參數Table 2 Tavares crushing model parameters

2.3.2 仿真時間步長理論計算

在EDEM 仿真中首先需要合理地設定時間步長,時間步長的設定需要考慮多個因素,如模擬的物理現(xiàn)象、計算機的性能、仿真時間等等。時間步長的大小對于仿真結果的準確性和計算速度都有很大的影響,如果設置時間步長過大,則會令仿真結果不收斂,并且對于硬球顆粒,過大的時間步長可能會使某些顆粒穿過設備,使得仿真結果不準確。如果設置的時間步長太小,則會要求計算機較高的配置并且耗費時間較長。因此,針對不同的仿真工況需要選取適當的時間步長。

顆粒在發(fā)生碰撞時,百分之七十的能量是通過瑞利波消耗的,因此臨界的仿真時間步長需要根據瑞利波的傳播速度決定。瑞利波是指當兩個單元發(fā)生碰撞時,因單元表面受到變應力而產生的沿單元表面?zhèn)鞑サ钠癫āＨ鹄▊鞑r,單元表面的能量最強,并且沿著球心方向上顯著減弱。各向同性材料中,瑞利波呈指數形式衰減;各向異性材料中,則隨著深度呈振蕩衰減,振蕩幅度的包絡線呈指數形式。通常情況下,如果傳播深度超過兩倍波長時,振幅就已非常小[22]。

彈性固體顆粒表層的瑞利波速度計算公式為:

式中,ν是顆粒的泊松比,則由上式求出β1的近似解為:

瑞利波速度進一步寫成:

由于兩接觸顆粒的接觸作用力僅限于發(fā)生在碰撞的兩顆粒上,而不應該通過瑞利波傳遞到其他顆粒上,因此時間步長必須小于瑞利波傳遞半球面的時間,故可求得時間步長的理論計算公式為:

當由不同粒徑的顆粒組成的系統(tǒng)時,時間步長應為:

上述公式是在兩顆粒為靜態(tài)或者相對運動速度較小的情況下得出的結論,由于高壓輥磨機的顆粒運行并不劇烈,故通過上式取得步長可使仿真模型穩(wěn)定。但是考慮到高壓輥磨機仿真破碎顆粒數量巨大,為使計算穩(wěn)定可取上式的1%～10%。

2.4 仿真參數設定

三山島金礦KHD375?1400-1000 高壓輥磨機設備參數和操作參數如表3所示。

表3 高壓輥磨機設備參數和操作參數Table 3 Equipment parameters and operating parameters of high-pressure grinding rolls

表4 為該仿真中使用的材料參數和接觸參數,這些參數參考了BARRIOS 等[23]和RODRIGUEZ 等[21]的研究,并根據EDEM 軟件中材料數據庫(通用材料模型數據庫-GEMM)的安息角結果進行驗證。

表4 離散元仿真參數Table 4 Discrete element simulation parameters

3 聯(lián)合仿真模型驗證與分析

為了驗證DEM 和MBD 聯(lián)合仿真在高壓輥磨機預測處理量、破碎比和單位能耗上的準確性,以及證明該聯(lián)合仿真方法的優(yōu)勢,以三山島金礦KHD375?1400-1000 型高壓輥磨機實際工作數據作為對比目標,對輥隙的動態(tài)特性進行驗證與分析,以及對前文所述的數學模型、單一離散元仿真模型和聯(lián)合仿真模型進行對比分析。試驗設備、控制系統(tǒng)和整體仿真模型如圖8所示,在EDEM 軟件中出料口位置設置GeometryBin(幾何選擇區(qū)域)和Mass Flow sensor(質量流傳感器),用于監(jiān)測破碎出料顆粒粒徑大小和排料流量。

高壓輥磨機進料的一個重要的特性是,需要一定高度的料柱給位于壓輥之間的顆粒床施加壓力[24]。為了避免進料的堵塞,根據CLEARY 等[25]的相關結論,通過DEM 模擬證明,料斗中的物料高度至少應等于其寬度,超過這個高度,物料額外的壓力就會逐步傳遞到料斗壁上,而不影響高壓輥磨機的性能。因此,仿真中采用的喂料高度為壓輥的寬度,只有當DEM 顆粒達到此高度時,仿真的結果才視為有效結果。并且在仿真過程中,保持相同的顆粒產生速度,以保持破碎過程的穩(wěn)定。

3.1 輥隙的動態(tài)變化

通過在Motionview 中添加動輥和定輥之間的位移輸出,從而得到輥隙在破碎過程中的動態(tài)變化,仿真結果具體如圖9所示。

圖9 高壓輥磨機輥隙變化Fig.9 Roll gap change of high-pressure grinding rolls

從圖9 可以看出,物料在0.5 s 時接觸壓輥,1.5 s 后喂料高度到達理論要求,在穩(wěn)定后輥隙保持在35.5 mm 左右。根據礦山提供的數據,破碎過程中KHD375?1400-1000 型高壓輥磨機的輥隙在30～40 mm之間,該仿真結果符合工程實際。通過DEMMBD 聯(lián)合仿真可以模擬高壓輥磨機輥隙的動態(tài)響應,研究不同工藝參數下輥隙的變化規(guī)律,能夠幫助設計人員更好地了解系統(tǒng)的性能和特性,優(yōu)化設計方案。

3.2 處理量

通過EDEM 導出的質量流量數據和根據理論計算數據以及高壓輥磨機工作數據,繪制各模型預測處理量,結果如圖10所示。

圖10 各模型處理量結果曲線Fig.10 Results curve of treatment volume of each model

比較了DEM-MBD 聯(lián)合仿真、DEM 仿真以及處理量數學模型各自的預測值與礦山實際處理量數據。由圖10 可以明顯看出,通過數學模型預測的處理量數據和實際工作數據偏差較大,可能是數學模型無法考慮到工作壓力、輥隙變化、喂料粒徑大小、物料物理性質等其他參數對設備的影響。DEM-MBD 仿真結果與DEM 仿真結果和實際工作數據比較接近,但DEM 仿真結果比實際數據偏低,原因可能是設備的處理量和輥隙有直接關系,在破碎過程中,輥隙因受到物料的作用力而發(fā)生改變,DEM 模型無法仿真輥隙的動態(tài)變化,造成誤差。

為減小誤差和避免偶然性,每個模型仿真2 次,各模型預測處理量數據誤差如表5所示。

表5 各模型處理量結果對比Table 5 Comparison of treatment volume results by each model

分析表5 可知,DEM-MBD 聯(lián)合仿真模型預測處理量準確度最高,其次是DEM 模型,數學模型準確度最低。DEM-MBD 聯(lián)合仿真模型預測值與實際測定值相比,雖然誤差更小,但仍然偏低,可能是因為在實際生產中高壓輥磨機輥面是柱釘面,而仿真中采取的是光滑面,降低了輥面對礦石的一部作用。在前期仿真試驗過程中發(fā)現(xiàn),數量大且尺寸較小的柱釘會極大增加仿真中細小單元格數量,影響仿真進度,考慮到計算機性能,在不對仿真結果造成較大誤差的情況下,將柱釘輥面簡化成光滑輥面。

3.3 單位能耗

在EDEM 軟件中導出壓輥的轉矩,再利用公式(8)和(9)求得設備的單位能耗,將此結果與單位能耗數學模型求得的單位能耗,以及工業(yè)生產中的實際單位能耗數據進行對比,結果如圖11所示。

圖11 各模型單位能耗曲線Fig.11 Unit energy consumption curve for each model

由圖11 可以看出,DEM-MBD 聯(lián)合仿真模型預測的單位能耗和實際數據擬合度非常高,DEM 仿真模型和數學模型預測值比較接近,但都比實際數據大?？赡苁且驗樵谄渌麠l件一致的情況下,單位能耗和輥隙大小成反比,輥隙越大,單位能耗越低。DEM模型和數學模型無法模擬輥隙的變化,所以造成預測值偏大。

仿真2 次,各模型預測單位能耗數據誤差見表6。

表6 各模型單位能耗結果對比Table 6 Comparison of energy consumption results per unit for each model

分析表6 可知,DEM-MBD 聯(lián)合仿真模型在預測單位能耗時,準確度最高,最小誤差僅有3.51%,具有比較滿意的結果,為后續(xù)預測單位能耗提供了一種更為準確的方式。

3.4 破碎比

表7 為各模型預測平均破碎比與實際設備測定平均破碎比的相對誤差。比較表中數據,發(fā)現(xiàn)DEMMBD 聯(lián)合仿真模型預測的平均破碎比與實際數據擬合度最好,其次是DEM 模型,但DEM 模型預測數據比實際數據要高,這是因為在高壓輥磨機中,壓輥通過高壓將物料擠壓成片,使其在輥隙中受到高壓力的擠壓和剪切作用。輥隙越小,物料在擠壓和剪切過程中受到的力度越大,破碎效果越好,破碎產品的細度也就越高。DEM 模型沒有考慮到輥隙的動態(tài)變化,一直以初始輥隙運行,所以造成破碎比偏大。

表7 各模型平均破碎比對比Table 7 Comparison of the average crushing ratio of each model

表8 分析總結了高壓輥磨機各性能指標的不同模型仿真結果與實際數據的平均誤差值,在初步建模階段,應該認為仿真結果誤差在一定范圍內并顯示正確的趨勢,則DEM-MBD 聯(lián)合仿真模型是成功的。

表8 各模型的性能指標相對誤差Table 8 Relative error of performance indicators for each model%

4 結論

(1)以KHD375?1400-1000 型高壓輥磨機為研究對象,為提高高壓輥磨機處理量、單位能耗、破碎比預測的準確度,本文首次將DEM-MBD 聯(lián)合仿真技術應用在高壓輥磨機虛擬仿真上,并詳細介紹聯(lián)合仿真方法和相關參數設定,研究表明DEM-MBD 仿真能夠更真實地模擬設備的運行狀態(tài)。

(2)為評價該聯(lián)合仿真模型對高壓輥磨機性能指標的預測準確度,將DEM-MBD 模型分別與DEM模型以及前文所述的處理量、單位能耗、破碎比數學模型和實際工作數據進行比較。通過對比處理量、單位能耗、平均破碎比三大性能指標,綜合考慮各模型的精確度和穩(wěn)定性,結果表明DEM-MBD 模型預測準確度最為理想。DEM-MBD 聯(lián)合仿真模型比現(xiàn)有的單一DEM 仿真模型在預測處理量、單位能耗、破碎比上的誤差值分別減小了13.41、38.10、6.74 個百分點。

(3)高壓輥磨機DEM-MBD 聯(lián)合仿真相比于DEM 仿真能夠更真實、更準確地模擬設備的運行狀態(tài)和性能,從而使高壓輥磨機仿真模型的準確度得到顯著改善,為此類仿真提供合理參考,降低了高壓輥磨機的研發(fā)成本,也為高壓輥磨機的設計和優(yōu)化提供了更科學的依據。然而,對高壓輥磨機產品粒徑分布的仿真還有未盡之處,后續(xù)將針對此開展進一步研究。