王學(xué)文 ，麻豪洲 ，李 博 ，張沛林 ，夏 蕊

（1.太原理工大學(xué) 機械與運載工程學(xué)院, 山西 太原 030024；2.太原理工大學(xué) 煤礦綜采裝備山西省重點實驗室, 山西 太原 030024）

0 引 言

刮板輸送機作為綜采工作面的輸運設(shè)備，承擔(dān)著煤巖散料的輸運工作，為適應(yīng)綜采智能化發(fā)展，刮板輸送機不斷向長運距、大運量、高可靠性、智能化方向發(fā)展[1]。但井下工作環(huán)境惡劣，難以安全有效開展相關(guān)研究，實驗室試驗由于客觀條件等限制無法充分再現(xiàn)原位條件。而仿真試驗可以有效彌補研究成本高、耗時長、安全性差的缺陷，也可人為控制試驗條件，成為一種有效的研究手段。學(xué)者們利用多種數(shù)值方法構(gòu)建機械與散料的剛散耦合仿真模型[2-3]，對煤機裝備進行相關(guān)的研究分析[4-6]。

離散元法由于其在分析散料的獨特優(yōu)勢，已經(jīng)廣泛應(yīng)用于煤機裝備的相關(guān)研究。學(xué)者們基于離散元法構(gòu)建刮板輸送機剛散耦合模型[7]，分析煤料與刮板輸送機的相互作用效應(yīng)，即剛散耦合效應(yīng)。XIA等[8-9]基于EDEM 研究了中部槽的磨損機理和磨損規(guī)律，分析了不同工況條件和煤料特性下中部槽的磨損。LI 和LIU[10]基于EDEM 研究了刮板、鏈條與中部槽在不同鏈間距、不同刮板距離、不同傾角下的受力及沖擊磨損特性。李康等[11]基于EDEM 研究了不同結(jié)構(gòu)形式和煤料參數(shù)下的磨損特性。李博等[12]基于EDEM 研究了煤料在中部槽的輸運特征。刮板輸送機各零部件之間存在復(fù)雜的運動關(guān)系，僅通過離散元法無法模擬輸運過程中非線性、時變等特性。為準確模擬刮板輸送機的運動，學(xué)者們基于多體動力學(xué)[13]通過構(gòu)建離散元-動力學(xué)模型分析研究刮板輸送機輸運過程的剛散耦合效應(yīng)。REN 等[14]基于EDEM 和ADAMS 建立刮板輸送機模型，研究了鏈條傳動系統(tǒng)在不同載荷和偏心載荷作用下的動力學(xué)特性。WANG 等[15]基于EDEM、RecurDyn 和ANSYS 建立了刮板輸送機中部槽模型，研究了裂紋深度、初始角度和拉伸載荷對應(yīng)力強度因子的單因素和交互影響。MA 等[16]基于EDME 和RecurDyn構(gòu)建了刮板輸送機模型，研究了刮板輸送機中部槽的力學(xué)效應(yīng)及磨損特性。

針對目前對于刮板輸送機的仿真研究多數(shù)仍采用單一軟件進行，構(gòu)建的仿真模型精度較低，同時考慮到現(xiàn)有研究忽視了煤料載荷的時變性及其與相關(guān)零部件的相互作用，為有效分析顆粒系統(tǒng)在設(shè)備運行中的機理，以地質(zhì)工況條件和運量工況條件作為不同工況條件，分別構(gòu)建離散元-多體動力學(xué)仿真模型，研究不同工作傾角、不同走向傾角和不同運量工況條件下刮板輸送機的剛散耦合效應(yīng)。研究可為后續(xù)開展不同工況條件下刮板輸送機輸運狀態(tài)監(jiān)測，關(guān)鍵零部件結(jié)構(gòu)優(yōu)化和改進提供參考，為后續(xù)開展刮板輸送機磨損特性和故障分析等研究從細觀角度分析煤料系統(tǒng)在刮板輸送機輸運過程中形成煤料-中板-刮板（鏈）三體接觸的機理提供模型構(gòu)建依據(jù)和仿真參考。

1 刮板輸送機剛散耦合模型構(gòu)建

1.1 刮板輸送機剛體模型構(gòu)建

1.1.1 刮板輸送機幾何模型

研究選用SGZ880/800 型號的刮板輸送機作為研究對象。為提高仿真效率，縮短研發(fā)周期，對仿真模型進行必要簡化。研究主要針對刮板輸送機剛散耦合效應(yīng)進行分析，簡化模型只保留首尾鏈輪、三節(jié)中部槽、刮板和鏈條。同時添加底槽，使刮板和鏈條可在鏈輪驅(qū)動下在中部槽內(nèi)平穩(wěn)運行。采用UG 構(gòu)建簡化后的刮板輸送機幾何模型（圖1），以刮板輸送機鋪設(shè)方向為y向，推移方向為x向，豎直方向為z向構(gòu)建坐標系。

圖1 刮板輸送機幾何模型Fig.1 Geometric model of scraper conveyor

1.1.2 刮板輸送機動力學(xué)模型

為真實模擬刮板輸送機鏈傳動系統(tǒng)的運動，將刮板輸送機幾何模型導(dǎo)入RecurDyn 中，并通過定義零部件的約束及其相關(guān)接觸參數(shù)施加邊界約束。將中部槽、鏈輪軸和底槽與地面添加固定副，首尾鏈輪處添加轉(zhuǎn)動副，用實體接觸副定義鏈環(huán)之間、鏈環(huán)與刮板、鏈環(huán)和刮板與鏈輪等部件間的接觸關(guān)系，通過定義接觸系數(shù)（剛度系數(shù)、阻尼系數(shù)、摩擦因數(shù)和穿透深度）使得各接觸零部件可以進行力和運動的傳遞，同時不會發(fā)生穿模，接觸系數(shù)根據(jù)文獻[17]設(shè)定。刮板輸送機各零部件接觸約束如圖2 所示。同時為滿足井下工況需要，刮板輸送機需要一定的預(yù)緊力。在鏈輪處添加移動副，通過改變鏈輪的中心距調(diào)節(jié)鏈條的預(yù)緊力。鏈條張緊后導(dǎo)出刮板輸送機動力學(xué)模型。并在鏈輪轉(zhuǎn)動副添加驅(qū)動函數(shù)，通過鏈輪與鏈環(huán)間的接觸、鏈環(huán)間的接觸和鏈環(huán)與刮板間的接觸實現(xiàn)鏈輪帶動鏈條（刮板）的運動過程。

圖2 刮板輸送機各零部件接觸約束Fig.2 Contact restraint of various parts of scraper conveyor

1.2 煤料散體模型構(gòu)建

煤料的形狀不規(guī)則，通過球形顆粒難以模擬真實井下工況。通過對破碎煤樣顆粒形狀進行統(tǒng)計觀察，顆粒形狀可以大致分為扁平狀、類塊狀和類錐狀[18]。在EDEM 中根據(jù)典型煤顆粒形狀構(gòu)建煤料散體仿真模型。由于煤礦井下煤層地質(zhì)條件不同，造成截割后的煤巖粒度范圍較大，研究為研究底板地質(zhì)工況和運量工況對刮板輸送機輸運過程的影響，根據(jù)GB/T 17608—2022《煤炭產(chǎn)品品種和等級劃分》標準選取中塊煤（粒度25～50 mm）為研究對象。結(jié)合構(gòu)建的煤顆粒模型（圖3），3 種形狀的煤顆粒平均粒度為40 mm。

圖3 煤顆粒模型形狀Fig.3 Particle model shape of coal

離散元通過定義參數(shù)賦予材料特性，為保證參數(shù)設(shè)定的準確性，需要通過相關(guān)標定試驗進行標定。離散元參數(shù)主要包括本征參數(shù)和接觸參數(shù)，研究中使用的參數(shù)根據(jù)文獻[19]標定獲得，具體參數(shù)見表1和表2。

表1 本征參數(shù)Table 1 Intrinsic parameters

1.3 刮板輸送機剛散耦合模型構(gòu)建及驗證

1.3.1 剛散耦合模型構(gòu)建

在EDEM 中靠近尾輪的中部槽上方處設(shè)置顆粒工廠用于產(chǎn)生煤料，模擬采煤機截割落料過程。顆粒工廠為850 mm×1 000 mm×200 mm 的立方體區(qū)域，位于中部槽上方900 mm 處。根據(jù)所研究的刮板輸送機型號確定配套的采煤機型號，根據(jù)文獻[20]計算出煤料脫離采煤機滾筒時的速度，并根據(jù)實際截割高度進行預(yù)仿真試驗確定在中板上方900 mm 處煤料的三向速度，通過定義顆粒工廠中顆粒的初速度實現(xiàn)落料。為進一步提高計算效率，在EDEM 中設(shè)定計算域為三節(jié)中部槽及包括顆粒工廠的上方區(qū)域。生成顆粒形狀根據(jù)圖3 設(shè)定，每種形狀各占比1/3，顆粒生成速率根據(jù)研究內(nèi)容設(shè)定。為真實模擬井下工況，在顆粒工廠位置處添加一個新的顆粒工廠用于生產(chǎn)矸石顆粒，矸石顆粒形狀根據(jù)圖3 的類塊狀設(shè)定。顆粒間的接觸模型選擇Hertz-Mindlin(no slip)模型。顆粒與幾何體間的接觸模型選擇Archard Wear 模型。研究結(jié)合選用刮板輸送機型號以及文獻[12]，設(shè)定滿載情況下中部槽的輸送量為220 kg/s，鏈速為1.1 m/s[12]，確定EDEM 中煤巖散料的顆粒生成速率為220 kg/s，為較真實符合井下工況，設(shè)定煤料基本生成速率為200 kg/s，矸石生成速度為20 kg/s，RecurDyn 中設(shè)定鏈輪角速度為3.62 rad/s。重力加速度為9.8 m/s2。

1.3.2 離散元模型驗證

為提高仿真的準確性，需要對構(gòu)建的模型進行驗證。研究基于ML-100C 改進型磨粒磨損機對離散元模型進行驗證（圖4）。該磨粒磨損試驗機經(jīng)過改進，可以模擬刮板在中部槽的運動[21]。在EDEM中，構(gòu)建相應(yīng)的試驗?zāi)Ｐ?，根?jù)表1 參數(shù)進行設(shè)定，進行仿真試驗。

圖4 磨粒磨損試驗機示意Fig.4 Abrasive wear testing machine

選取試樣受力和煤料輪廓曲線作為驗證指標，對比結(jié)果如圖5 所示。相比試驗的試樣受力曲線，仿真曲線由于顆粒與試樣在瞬時產(chǎn)生擠壓，存在一些突變點，而在試驗中由于傳感器的采樣頻率較低，這些突變點沒有被捕捉到。根據(jù)采集到的試驗和仿真數(shù)據(jù)（圖5a），分別求取平均值，試驗平均值為0.94 N，仿真的平均值為0.82 N，誤差為12.77%。選取相同位置的煤料截面輪廓曲線（圖5b），分別對試驗和仿真數(shù)據(jù)點進行曲線擬合，并對試驗曲線和仿真曲線求取相關(guān)系數(shù)，相關(guān)系數(shù)為0.989 6。試驗和仿真的試樣受力以及煤料輪廓曲線基本接近，可認為離散元模型滿足研究要求[22]。

圖5 離散元模型驗證試驗結(jié)果Fig.5 Result of discrete model verification test

1.3.3 動力學(xué)模型驗證

根據(jù)前文構(gòu)建的刮板輸送機動力學(xué)模型進行空載試驗仿真，選取尾輪處刮板的位移和速度進行分析（圖6）。通過刮板的y向位移可以看出，刮板運行較為平穩(wěn)，且在6 s 仿真時間內(nèi)可滿足刮板完成一次運煤。在啟動階段，鏈輪角速度逐漸增加至3.62 rad/s后保持不變，刮板的速度也從0 增大到1 m/s 附近波動。當(dāng)刮板開始進入首輪時，速度曲線減小，并反向增大至-1 m/s 附近。符合預(yù)設(shè)工況下的真實變動，可認為刮板輸送機的動力學(xué)模型滿足研究要求[22]。

圖6 刮板的y 向位移和速度曲線Fig.6 The y displacement and velocity curve of scraper

1.4 不同工況條件剛散耦合仿真設(shè)定

刮板輸送機工作環(huán)境十分惡劣，運行過程中會受地理條件、煤質(zhì)因素、機械因素等，呈現(xiàn)出不同的剛散耦合效應(yīng)。研究針對不同工作傾角、不同走向傾角（圖7）和不同運量工況，研究分析刮板輸送機的剛散耦合效應(yīng)，圖7 中，α為工作傾角，β為走向傾角。

圖7 刮板輸送機鋪設(shè)示意[23]Fig.7 Diagram of scraper conveyor laying[23]

研究定義水平運輸工況（工作傾角0°）、無傾斜（走向傾角0°）、煤料生產(chǎn)速率200 kg/s 作為正常工況，通過調(diào)節(jié)重力加速度分量實現(xiàn)對鋪設(shè)條件（工作傾角和走向傾角）的改變，鋪設(shè)工況條件設(shè)定見表3，示意如圖8 所示。調(diào)節(jié)顆粒工廠的生產(chǎn)速率實現(xiàn)對運量的改變，運量分別設(shè)置為200、300 和400 kg/s。

表3 鋪設(shè)工況條件設(shè)定Table 3 Setting of working conditions for laying

圖8 鋪設(shè)工況條件示意Fig.8 Diagram of working conditions for laying

2 不同工況條件下運動學(xué)效應(yīng)

2.1 下運工況

井下煤層底板條件復(fù)雜，刮板輸送機往往難以水平鋪設(shè)，本節(jié)對不同工作傾角（下運工況）和不同走向傾角（傾斜工況）進行仿真，選取第三節(jié)中部槽作為研究對象，分析不同鋪設(shè)條件工況下的運動學(xué)效應(yīng)。

刮板帶動煤料實現(xiàn)刮板輸送機對煤料的輸送，而刮板輸送機沿工作傾角的傾斜時，會由于重力影響煤料的運輸阻力，進而對煤料的運動產(chǎn)生影響。正常工況下煤料分布如圖9a—圖9c 所示，可以看出，刮板前側(cè)煤料堆積較高，在后側(cè)則形成明顯的凹陷（圖9a），煤料在中部槽的中間堆積最高，達到煤料休止角后，煤料向兩側(cè)槽幫流動（圖9b），且刮板后側(cè)煤料也明顯少于前側(cè)（圖9c），兩節(jié)刮板之間煤料的運動具有相似的趨勢，煤料整體呈現(xiàn)堆狀間歇分布。以圖9a 所示，測量最高堆積點位置，L為最高點距后一節(jié)刮板的距離，H為最高點距中板的距離，測量結(jié)果如圖10 所示，可以看出不同下運工況條件下H變化基本相同；而L則隨著工作角度逐漸減小，這是由于當(dāng)堆積煤料達到休止角后，煤料最高點H基本不變，但由于重力分量的影響，下運工況會使得堆積煤料向前運動，進而會改變原來的堆積形貌形成一個新的堆積，使煤料堆積最高點逐漸后移。

圖9 正常工況煤料分布Fig.9 Coal distribution under normal conditions

圖10 下運工況最高點位置Fig.10 Position of highest point in downward conditions

2.2 傾斜工況

傾斜工況下煤料分布如圖11 所示，可以看出，正常工況條件下，煤料堆積最高處基本位于中部槽的中線位置（圖9b），而傾斜工況下煤料則向傾斜一側(cè)堆積。分別測量煤料最高點距中板距離H和偏移距離Q，如圖12 所示，可以看出，最高點距中板的距離與下運工況類似，基本保持不變，而偏移距離則由于重力分量的影響，煤料向傾斜側(cè)有運動趨勢，導(dǎo)致最高點向傾斜側(cè)偏移。選取鏈條兩側(cè)區(qū)域煤料進行分析，傾斜工況下煤料質(zhì)量曲線如圖13 所示，可以看出正常工況下兩區(qū)域煤料質(zhì)量基本相同，隨著傾角增大，A 區(qū)域內(nèi)的煤料少于B 區(qū)域，且呈現(xiàn)一定的滯后性，中部槽內(nèi)煤料的分布不均。

圖11 傾斜工況煤料分布Fig.11 Coal distribution under inclined condition

圖12 傾斜工況最高點位置Fig.12 Position of highest point under inclined condition

圖13 煤料質(zhì)量曲線Fig.13 Quality curve of coal

2.3 不同運量工況

取刮板輸送機中間3 個刮板間的煤料分布（圖14），可以看出運量200 kg/s 時煤料在中部槽內(nèi)呈現(xiàn)堆狀間歇分布，隨著運量提高，煤料分布由堆狀間歇分布轉(zhuǎn)為階梯狀分布。不同運量工況下通過刮板輸送機中間位置的質(zhì)量流率如圖15a 所示，可以看出由于刮板帶動煤料運動，質(zhì)量流率存在明顯的周期性變化，200 kg/s 的質(zhì)量流率最大值在200 kg/s附近與運量基本一致，而300 kg/s 和400 kg/s 的質(zhì)量流率則沒有達到設(shè)定運量。結(jié)合煤料y向速度云圖（圖15b）可以看出煤料中間紅色區(qū)域表示煤料運動速度與鏈速基本一致（1 m/s），而藍色區(qū)域煤料基本為靜止狀態(tài)，即中間紅色區(qū)域為煤料輸送的有效輸送區(qū)域。由于運量增多，導(dǎo)致煤料無法有效輸送而在機尾處和槽幫處形成堆積，使得煤料分布呈現(xiàn)階梯狀，質(zhì)量流率達不到設(shè)定運量。

圖14 運量工況煤料分布Fig.14 Coal distribution under traffic condition

圖15 煤料質(zhì)量流率和y 向速度Fig.15 Coal mass flow rate and y-direction velocity

3 不同工況條件下力學(xué)效應(yīng)

3.1 下運工況

3.1.1 煤料受力情況

正常工況下煤料間受力如圖16 所示，可以看出中部槽內(nèi)的煤料受力明顯大于堆積在槽幫上部的煤料，刮板前側(cè)煤料受力大于刮板后側(cè)，鏈條兩側(cè)煤料受力大于鏈條中間。中部槽內(nèi)煤料會受到刮板、中板和槽幫的相互作用力，在刮板帶動下，刮板前側(cè)的煤料向前運動，同時推動上一節(jié)刮板后側(cè)的煤料，而堆積在上部的煤料則是依靠底部煤料施加的摩擦力向前運動，因此小于中部槽內(nèi)煤料的受力。中板與槽幫會對煤料施加向后的摩擦力，同樣會導(dǎo)致刮板后側(cè)煤料受力小于刮板前側(cè)。鏈條兩側(cè)相比鏈條間的空間更大，有利于煤料的填充和流動，使得鏈條間的煤料的作用力較小。下運工況改變了重力對煤料的作用力方向，重力的分量推動煤料的輸運。下運工況煤料受力如圖17 所示，可以看出隨著工作傾角不斷增大，刮板前側(cè)紅色顆粒和刮板后側(cè)藍色顆粒數(shù)目逐漸減小，刮板間綠色顆粒占比逐漸增加，說明刮板前側(cè)煤料的受力逐漸減小，煤料向前一個刮板后方運動，中部槽內(nèi)煤料的載荷分布更為均勻。

圖16 正常工況煤料受力Fig.16 Stress of coal under normal conditions

圖17 下運工況煤料受力Fig.17 Stress of coal under downward conditions

3.1.2 中板受力情況

1）煤料與中板。下運工況下第三節(jié)中板受力時間和平均受力見表4，初始受力時間是指中板第一次檢測到受力的時間，此刻僅有零散煤料落在中板，中板受力在零值波動，穩(wěn)定初始時間是中板受力持續(xù)開始的時間，此刻煤料進入穩(wěn)定運行階段，即此時刻后中板受力不存在零值數(shù)據(jù)點。可以看出隨著工作傾角的增大，初始受力時間和穩(wěn)定初始時間逐漸減小，且兩者的差值逐漸增大，這是由于重力分量使得煤料更早的向前輸運。工作傾角較小時，中板平均受力呈現(xiàn)波動趨勢，隨著工作傾角的增大，中板平均受力呈現(xiàn)上升趨勢，這是由于傾角較小時，重力的分量對煤料影響較小，煤料運動的隨機性較大，當(dāng)傾角進一步提升時，重力的分量對煤料影響大于煤料運動的隨機性影響，中板平均受力與煤料運動呈現(xiàn)一致性。

表4 中板受力時間和平均受力值Table 4 Stress time and average force of deck plate

2）刮板與中板。選擇顆粒工廠正下方的刮板作為研究對象，刮板與中板的受力如圖18 所示?？蛰d工況下，鏈條帶動刮板在中板上運動，由于鏈傳動系統(tǒng)的多邊形效應(yīng)，使得刮板與中板的受力呈現(xiàn)波動狀態(tài)。運載工況下刮板與中板的接觸力也呈現(xiàn)波動狀態(tài)，但刮板與中板的接觸力并不是一直存在。下運10°工況同一刮板的受力和Z向位移如圖19 所示，結(jié)合圖9b 可以看出，在運載工況下由于鏈傳動的多邊形效應(yīng)會使得煤料逐漸填充到刮板下方，進而使得刮板z向受力出現(xiàn)一個瞬時峰值，刮板向上方（z向）移動，與中板脫離接觸。

圖18 刮板與中板受力Fig.18 Force of scraper and deck plate

圖19 刮板受力與z 向位移Fig.19 Force and z-directional displacement of scraper

3.1.3 刮板及鏈條受力情況

下運工況下某刮板受力情況如圖20 所示，可以看出刮板開始輸運煤料后，受力逐漸增大，隨著下運角度的增加，刮板受力呈現(xiàn)滯后性，這是由于部分煤料在重力分量的影響下更早向前輸運導(dǎo)致，與中板初始受力時間趨勢相反；在到達落料區(qū)時受力達到最大值，之后在某一水平波動。由圖20 可以看出4.5 s 以后該刮板卸載煤料，為直觀分析刮板受力情況，分別計算刮板在0～4.5 s 內(nèi)的平均受力值，隨著工作傾角的增大，刮板的平均受力分別為1 212.97、920.53、 907.11、 821.48 N。這是由于重力分量對煤料運動的影響，使得煤料有向前運動的趨勢，導(dǎo)致刮板平均受力逐漸減小，與中板平均受力呈現(xiàn)相反的趨勢。同時可以看出由于鏈傳動系統(tǒng)多邊形效應(yīng)的影響，刮板受力也呈現(xiàn)一定的周期性波動，之后到達卸料區(qū)，受力減小至0。

圖20 下運工況刮板受力Fig.20 Force on scraper under downward conditions

隨著下運角度的增加，鏈條受力會呈現(xiàn)減小的趨勢。下運工況下經(jīng)濾波平滑處理后的鏈條受力情況如圖21 所示，可以看出鏈條受鏈傳動的多邊形效應(yīng)影響，呈現(xiàn)明顯的周期性，鏈條受力并沒有明顯的變化趨勢。根據(jù)圖21 中下運15°時4 s 后鏈條受力曲線可以看出，在輸運過程中部分煤料會卡在刮板（鏈）和中板間會導(dǎo)致鏈條受力突然增加。

圖21 下運工況鏈條受力Fig.21 Force on chain under downward conditions

3.1.4 中板磨損情況

下運工況下中板磨損情況如圖22 所示，磨損深度為中板所有單元磨損的總和?？梢钥闯龀跗? 種工況磨損基本相同，這是由于此時中部槽內(nèi)煤料較少，煤料和中板只形成兩體磨損，磨損量較小。1.6 s 左右時，正常工況下中板磨損深度出現(xiàn)一個階梯上升，這是由于在運行過程中鏈條和刮板帶動煤料向前輸送，煤顆粒-中板-刮板（鏈）形成三體接觸進而造成三體磨損導(dǎo)致[16]。5 s 時刻正常工況與下運5°工況磨損深度基本相同，下運10°工況磨損深度略低，這是由于重力的分量使得刮板輸送機上堆積煤料有向前運動的趨勢，煤料在中部槽內(nèi)分布均勻，形成三體磨損的顆粒較少。下運15°在4 s 后的磨損深度出現(xiàn)明顯上升趨勢，由圖21 可知是由于部分煤料卡在刮板（鏈）和中板之間導(dǎo)致磨損深度急速增加。

圖22 下運工況磨損深度Fig.22 Wear depth under downward conditions

3.2 傾斜工況

3.2.1 煤料受力情況

傾斜10°工況下煤料間受力如圖23 所示，可以看出鏈條兩側(cè)煤料受力情況不一致，B 區(qū)域煤料受力大于A 區(qū)域。考慮到煤料生成的隨機性影響，選用A 區(qū)域和B 區(qū)域平均受力比值（FA/FB）進行對比，正常工況下FA/FB=96%，傾斜5°工況下FA/FB=86%，傾斜10°工況下FA/FB=70%，隨著走向傾角的不斷增大，煤料向傾斜側(cè)堆積，造成傾斜側(cè)煤料受力大于另一側(cè)。

圖23 傾斜10°工況煤料受力Fig.23 Stress of coal under the condition of inclination of 10 °

3.2.2 鏈條受力情況

煤料載荷在x方向分布不均導(dǎo)致兩條鏈條受力也呈現(xiàn)不均勻情況，傾斜工況下鏈條受力如圖24 所示，鏈2 為靠近傾斜側(cè)的鏈?？梢钥闯鲈谡９r下，鏈條受力基本同步，并呈現(xiàn)周期性變化。隨著走向傾角的增大，鏈條受力在初始階段基本同步；在穩(wěn)定運輸階段，鏈條受力出現(xiàn)波動，且兩條鏈受力呈現(xiàn)較大的差異，鏈2 受力明顯高于鏈1；在卸料階段，由于靠近鏈輪，鏈條受力的又開始呈現(xiàn)同步性。

圖24 傾斜工況鏈條受力Fig.24 Chain stress under inclined working condition

3.2.3 鏈條與鏈輪受力情況

選取同一對鏈環(huán)分別分析正常工況與傾斜10°工況條件下鏈條與鏈輪的接觸受力（圖25）。在0～0.9 s 內(nèi)鏈環(huán)與鏈輪嚙合，由于此時中部槽內(nèi)生成煤料較少，鏈條與鏈輪的接觸力基本呈現(xiàn)一致的波動性。在0.9～5 s 內(nèi)鏈條帶動刮板輸送煤料，此時鏈環(huán)與鏈輪脫離接觸，接觸力為0。5.35 s 后刮板卸料，鏈環(huán)與鏈輪嚙合，此時中部槽內(nèi)鋪滿煤料，鏈條與鏈輪的接觸力波動幅度相較于之前變大，正常工況下兩條鏈和傾斜10°工況下鏈1 波動較為接近，傾斜10°工況條件下鏈2 的波動明顯高于鏈1。分別計算正常工況和傾斜10°工況下兩條鏈與首尾鏈輪接觸力絕對差值的平均值，正常工況下鏈條與首尾鏈輪的平均差值為54.13 N，傾斜工況為18 015.64 N。結(jié)合圖24c 可知，由于傾斜工況使得煤料在中部槽內(nèi)分布不均，導(dǎo)致刮板卸料后兩條鏈仍存在一定的張力差，使得鏈條與鏈輪的接觸力也存在一定的差值。

圖25 鏈條與鏈輪受力Fig.25 Force of chain and sprocket

3.3 不同運量工況

3.3.1 中板受力情況

不同運量工況下中部槽受力時間和平均受力見表5?？梢钥闯龅谝还?jié)中板初始受力時間基本一致，由于運量的增加，300 kg/s 和400 kg/s 工況下煤料很快與中板穩(wěn)定接觸。400 kg/s 工況下第一節(jié)中板平均受力小于300 kg/s，是由于運量較大，一部分煤料在槽幫處形成了堆積，且根據(jù)第二節(jié)中板初始受力時間可以看出，一部分煤料直接落在了第二節(jié)中板。第二節(jié)中部槽內(nèi)中板和槽幫的受力都隨著運量的增大而增大。第三節(jié)中板在300 kg/s 工況下受力最大，結(jié)合2.3 節(jié)對不同運量工況煤料的運動分析可知，煤料輸運存在有效輸運區(qū)域，400 kg/s 由于運量過大，導(dǎo)致煤料無法有效向前輸運，在槽幫處形成堆積，使得第三節(jié)中板受力小于300 kg/s 工況。

表5 中部槽受力時間和平均受力值Table 5 Stress time and average force of line pan

3.3.2 刮板受力情況

運量400 kg/s 時刮板受力如圖26 所示，依據(jù)刮板推動煤料的順序，分別記為刮板1、刮板2、刮板3、刮板4。刮板開始推動煤料時間依次增大，刮板受力也呈現(xiàn)增大趨勢。分別提取刮板受到煤料阻力后1.5s 內(nèi)的阻力平均值和對應(yīng)的中部槽內(nèi)煤料的質(zhì)量流率數(shù)據(jù)（表6），可以看出刮板受力的增大是由于中部槽內(nèi)煤料質(zhì)量的增加。對表6 數(shù)據(jù)利用Origin 進行擬合，擬合結(jié)果為y=0.320 7x2-131.2x+15 999，R2=0.994，即隨著煤料的增加，刮板阻力的增速逐漸增大。

表6 刮板平均受力及質(zhì)量流率Table 6 Average force of scraper and mass flow rate

圖26 運量400 kg/s 刮板受力Fig.26 Force on scraper under 400 kg/s traffic condition

3.3.3 中板磨損情況

運量300 kg/s 和400 kg/s 時中板磨損和受力如圖27 所示。隨著運量的增加，磨損深度顯著增加，且存在明顯的階梯上升趨勢。結(jié)合中板受力對比，中板磨損急速增加與中板受力的波動相關(guān)，這是由于運量的增加，煤顆粒-中板-刮板（鏈）形成三體磨損的概率增加，加劇了中部槽的磨損，當(dāng)卡在刮板（鏈）和中板之間的煤料在輸運過程中姿態(tài)改變，使得煤顆粒-中板-刮板（鏈）之間形成兩體接觸，中板受力出現(xiàn)明顯降低，中板磨損平穩(wěn)增加。

圖27 不同運量下中板磨損和受力Fig.27 Wear and force on the deck plate underdifferent traffic volumes

4 結(jié) 論

1）基于離散元法和多體動力學(xué)構(gòu)建的刮板輸送機剛散耦合模型可以有效模擬刮板輸送機的輸運工況。正常工況條件下，中部槽內(nèi)煤料運動呈現(xiàn)堆狀間歇分布，且刮板前側(cè)煤料堆積多于刮板后側(cè)。中部槽內(nèi)的煤料受力明顯大于堆積在槽幫上部的煤料，刮板前側(cè)煤料受力大于刮板后側(cè)，鏈條兩側(cè)煤料受力大于鏈條中間。

2）下運工況條件下，煤料堆積最高點距中板距離基本不變，隨著工作傾角的增大，煤料最高點逐漸向著后側(cè)刮板移動，中板平均受力逐漸增大，后側(cè)刮板平均受力逐漸減小。

3）傾斜工況條件下，煤料堆積最高點距中板的距離基本不變，隨著走向傾角的增大，最高點逐漸向傾斜側(cè)偏移，造成中部槽內(nèi)煤料的分布不均，使得傾斜側(cè)煤料受力大于另一側(cè)，鏈條受力出現(xiàn)波動，兩條鏈同步性較差。

4）刮板輸送機輸送過程存在有效輸送區(qū)域，煤料速度在有效輸送區(qū)域內(nèi)與鏈速基本相同，靠近槽幫處的煤料速度逐漸越低。中板和刮板的平均受力隨著運量增多而增多，但運量過大時，中板的受力會由于煤料無法有效輸送堆積在槽幫處而減小。

5）中板磨損與煤顆粒-中板-刮板（鏈）的接觸相關(guān)，當(dāng)三者形成有效的三體接觸時，中板受力和磨損會出現(xiàn)明顯上升的趨勢。隨著輸運過程中煤顆粒位姿形態(tài)的改變，煤顆粒-中板-刮板（鏈）由三體接觸變?yōu)閮审w接觸，中板受力降低，磨損平穩(wěn)增加。