王 慧 宋宇寧
遼寧工程技術(shù)大學(xué)機械工程學(xué)院,阜新,123000
實現(xiàn)采煤機滾筒的自動調(diào)高對實現(xiàn)采煤工作面無人化具有重要意義。與此同時,在實現(xiàn)工作面煤層和巖層界面的自動識別后,可以通過自動調(diào)高滾筒避開巖石,避免滾筒截齒與巖石碰撞,延長采煤機的使用壽命[1]。
文獻[2-4]對采煤機自動調(diào)高控制及其關(guān)鍵技術(shù)進行了研究,提出了一種基于模變結(jié)構(gòu)控制的自動調(diào)高控制策略,并搭建了調(diào)高系統(tǒng)模擬加載實驗臺進行實驗研究。文獻[5]研究了一種基于地球物理探測、位姿測量技術(shù)以及地學(xué)信息建模多信息融合的自動調(diào)高方法,并以山西某煤礦為例對所研究的調(diào)高方法進行了應(yīng)用,應(yīng)用結(jié)果表明,該調(diào)高方法的精度能夠滿足現(xiàn)場需求。文獻[6]基于雙坐標系的采煤機截割路徑控制策略來處理采煤機截割路徑不平整問題,研究了基于微分理論的搖臂傾角控制算法以實現(xiàn)滾筒高度變化的連續(xù)控制。文獻[7]對采煤機滾筒高度變化與液壓缸活塞桿位移的關(guān)系進行了數(shù)學(xué)建模,建立了折算到油缸上的負載和滾筒搖臂質(zhì)量的函數(shù)關(guān)系,通過仿真方法對調(diào)高機構(gòu)的運動特性進行研究。
目前專家學(xué)者針對采煤機調(diào)高過程的控制策略進行深入研究,期望快速、穩(wěn)定、精準地控制采煤機滾筒高度。本文針對不同牽引速度和不同煤巖截割阻抗情況下,使用常規(guī)PID和模糊PID控制器對采煤機調(diào)高過程的響應(yīng)進行分析研究。
如圖1所示,采煤機調(diào)高機構(gòu)主要由調(diào)高油缸、搖臂、電液比例換向閥、油泵、電機、安全閥、過濾器等組成。通過電液比例換向閥控制液壓缸的伸出與縮回,從而調(diào)節(jié)搖臂的高度。
圖1 采煤機調(diào)高機構(gòu)基本組成Fig.1 Basic components of shearer height adjustment mechanism
圖2 調(diào)高機構(gòu)運動軌跡示意圖Fig.2 Schematic diagram of motion track of height adjusting mechanism
根據(jù)滾筒采煤機調(diào)高機構(gòu)的實際工作情況,對調(diào)高機構(gòu)運動學(xué)進行分析。為了便于分析,首先分析采煤機靜止情況下,調(diào)高液壓缸伸出時的截齒運動軌跡,此時的調(diào)高機構(gòu)運動軌跡如圖2所示。圖2中,將坐標系OXY的原點選取在搖臂的回轉(zhuǎn)軸心上,隨動坐標系O1X1Y1表示螺旋滾筒回轉(zhuǎn)中心的起始位置,調(diào)高液壓缸以速度vu伸出,同時,螺旋滾筒以角速度ω逆時針回轉(zhuǎn),經(jīng)過時間t后,隨動坐標系O1X1Y1變換到O2X2Y2。
根據(jù)上述描述,容易得到經(jīng)過時間t后,滾筒運動的位置方程:
(1)
式中,l1為搖臂的回轉(zhuǎn)中心到螺旋滾筒回轉(zhuǎn)中心的距離,m;ω2為搖臂的擺動的角速度,rad/s;θ為搖臂與水平面的夾角,rad。
根據(jù)截割部的調(diào)高原理,可以知道,搖臂的擺動角速度ω2與調(diào)高液壓缸的伸出速度vu有關(guān),根據(jù)三角形原理可以得到
(2)
(3)
式中,θ1為螺旋滾筒位于初始位置時,小搖臂與搖臂回轉(zhuǎn)中心和調(diào)高液壓缸安裝軸心連線的夾角,(°);θ2為經(jīng)過t時間后,小搖臂與搖臂回轉(zhuǎn)中心和調(diào)高液壓缸安裝軸心連線的夾角,(°);l2為小搖臂的長度,m;l3為搖臂回轉(zhuǎn)中心和調(diào)高液壓缸安裝軸心連線長度,m;l0為螺旋滾筒位于起始位置時,調(diào)高液壓缸的長度,m。
根據(jù)式(2)、式(3)便可得到搖臂擺動角速度:
(4)
之后,分析采煤機行走情況下,調(diào)高液壓缸伸出時,滾筒的運動軌跡:
(5)
式中,vq為采煤機的牽引速度。
采煤機調(diào)高機構(gòu)油缸是非對稱油缸,假設(shè)換向閥為理想滑閥,調(diào)高機構(gòu)中,閥控非對稱油缸進行調(diào)高的原理如圖3所示[8-9]。
圖3 閥控非對稱油缸調(diào)高原理圖Fig.3 Valve control asymmetric cylinder height adjustment schematic
下面根據(jù)調(diào)高系統(tǒng)原理圖建立調(diào)高系統(tǒng)非對稱閥控缸的數(shù)學(xué)模型:
(1)建立油缸流量連續(xù)性方程。無桿腔流量連續(xù)性方程為
(6)
有桿腔流量連續(xù)性方程為
(7)
式中,y為油缸活塞位移量;K為液壓油的彈性模量;Q1、Q2分別為油缸的進出油腔體積流量;A1、A2分別為油缸的無桿腔和有桿腔的作用面積;Cic、Cec分別為油缸的內(nèi)外泄漏系數(shù);p1、p2分別為油缸的無桿腔和有桿腔的液壓油壓力;V1、V2分別為油缸的無桿腔和有桿腔的容積。
設(shè)負載壓力pL和流量QL分別為
(8)
則有
(9)
式中,Vt為油缸無桿腔初始容積。
(2)建立油缸和負載的力平衡方程:
(10)
式中,F(xiàn)g為油缸輸出力;m為活塞質(zhì)量;B為油缸和負載的等效阻尼;k為負載等效剛度;FL為負載擾動力。
分析調(diào)高機構(gòu)的工作原理及負載可知,負載擾動力FL由五部分組成:
FL=Fe1+Fe2-Fe3-Fe4+Fe5
(11)
式中,F(xiàn)Y(t)為滾筒推進阻力;FZ(t)為滾筒截割阻力;G1為滾筒重力;G2為搖臂重力。
(3)建立滑閥壓力流量方程:
QL=Kqxv-KcpL
(12)
式中,Kq為滑閥流量增益;Kc為滑閥流量壓力系數(shù)[10]。
(4)調(diào)高機構(gòu)的數(shù)學(xué)模型建立。將上述得到的油缸流量連續(xù)性方程、油缸和負載的力平衡方程以及滑閥壓力流量方程進行拉氏變換:
(13)
ps=p1+p2/m2
Ctc=(1+m)Cic/(1+m3)+Cec/(1+m2)
式中,ps為供油壓力。
可獲得調(diào)高機構(gòu)數(shù)學(xué)模型:
(14)
其中,ωh為液壓固有頻率;ξh為液壓阻尼比;A為煤巖截割阻抗均值,如果小于180 kN/m,則為軟煤,如果在180~240 kN/m范圍內(nèi),則為中硬度煤,如果在240~360 kN/m范圍內(nèi),則為硬度煤。
將上述調(diào)高機構(gòu)數(shù)學(xué)模型作為調(diào)高機構(gòu)控制器的控制對象。
滾筒截割負載主要載荷的形成如圖4所示。采煤機滾筒在截割煤巖過程中,受到側(cè)向力FX(t)、推進阻力FY(t)、截割阻力FZ(t)及產(chǎn)生的截割阻力矩M(t)。側(cè)向力FX(t)不是調(diào)高機構(gòu)的主要負載,因此本文不予以考慮。
圖4 滾筒截割負載主要載荷的形成Fig.4 Formation of main load of the drum cutting load
推進阻力、截割阻力及產(chǎn)生的截割阻力矩分別為[11-12]
(15)
式中,φi為第i個截齒與z軸的夾角;N為參與截割的截齒數(shù)量;Dc為滾筒直徑;FZi為第i個截齒的截割阻力;FYi為第i個截齒的牽引阻力,F(xiàn)Yi=(0.5~0.8)FZ0+100KySdσy;Sd為截齒磨損面積;f為阻力系數(shù);Ky為單向抗壓強度比值;σy為煤巖單軸抗壓強度。
可以看出,F(xiàn)Zi和FYi均用第i個截齒的截割阻力均值FZ0表示,如果同時考慮采煤機牽引速度和煤巖截割阻抗,則FZ0為
(16)
式中,Kf為截齒前刃面形狀系數(shù);Kp為截齒配置系數(shù);Ka為截角影響系數(shù);Km為煤巖體裸露系數(shù);h為切屑厚度[13-14]。
使用常規(guī)PID和模糊PID這兩種易于工程實現(xiàn)的控制器,對不同滾筒截割載荷和牽引速度情況下的調(diào)高機構(gòu)進行控制。使用MATLAB編寫滾筒截割阻力程序,獲取調(diào)高機構(gòu)負載變化規(guī)律,將數(shù)據(jù)導(dǎo)入仿真模型。
本文將MG400/930-WD型采煤機作為研究對象。截割煤巖阻力系數(shù)f=0.40;單向抗壓強度比值Ky=1.0;通過整定后的PID參數(shù)Kp=2.8,Ki=0.25,Kd=15.8。
截割阻抗A分別選取為100 kN/m(軟煤巖),200 kN/m(中硬度煤巖)和300 kN/m(硬度煤巖),牽引速度分別設(shè)定為2 m/s、3 m/s和4 m/s。
通過仿真得到不同情況(牽引速度分設(shè)定為2 m/s,截割阻抗A分別選取為100 kN/m、200 kN/m和300 kN/m,以及截割阻抗A取200 kN/m,牽引速度分別設(shè)定為2 m/s、3 m/s和4 m/s,滾筒受到的Y方向和Z方向載荷以及Z方向的轉(zhuǎn)矩的變化情況如圖5、圖6所示。仿真結(jié)果得到的規(guī)律是,隨著截割阻抗以及牽引速度的增大,滾筒的截割載荷增大,因此本文將研究不同情況下,不同截割載荷對調(diào)高機構(gòu)控制精度的影響。
(a)A=100 kN/m
(b)A=200 kN/m
(c)A=300 kN/m
圖5不同截割阻抗情況下的滾筒載荷(vq=2 m/s)
Fig.5Drumloadwith2m/sanddifferentcuttingimpedance
使用AMESim仿真軟件建立采煤機調(diào)高機構(gòu)的仿真模型,如圖7所示。
AMESim仿真軟件能夠根據(jù)采煤機調(diào)高機構(gòu)的實際物理模型建立仿真模型,滾筒載荷以及調(diào)高機構(gòu)的控制器參數(shù)通過聯(lián)合仿真接口由MATLAB軟件獲得。使用正弦波和方波作為調(diào)高機構(gòu)控制器的設(shè)定信號,首先將牽引速度設(shè)定為2 m/s,截割阻抗A分別選取100 kN/m、200 kN/m和300 kN/m,之后設(shè)定截割阻抗A為200 kN/m,牽引速度分別設(shè)定為2 m/s、3 m/s和4 m/s,研究各種情況時,常規(guī)PID控制器和模糊PID控制器作用下,調(diào)高機構(gòu)的滾筒位移響應(yīng)如圖8~圖12所示。
(a)vq=2 m/s
(b)vq=3 m/s
(c)vq=4 m/s
圖6不同牽引速度下的滾筒載荷(A=200 kN/m)
Fig.6Drumloadwithdifferenttractionspeedand200kN/mcuttingimpedance
圖7 采煤機調(diào)高機構(gòu)的仿真模型Fig.7 Simulation model of shearer height adjusting mechanism
仿真結(jié)果表明,牽引速度為2 m/s,截割阻抗A分別選取為100 kN/m、200 kN/m和300 kN/m,以及截割阻抗A取200 kN/m,牽引速度分別設(shè)定為3 m/s和4m/s這樣五種情況下,使用模糊PID控制器對調(diào)高機構(gòu)進行控制時,調(diào)高機構(gòu)輸出的滾筒位置響應(yīng)對設(shè)定曲線具有較好的跟蹤性能,誤差較小,并且其控制效果受截割阻抗和牽引速度的變化影響不大,能夠?qū)Σ煌r有較好的適應(yīng)性。
(a)A=100 kN/m
(b)A=200 kN/m
(c)A=300 kN/m
圖8不同截割阻抗情況下的滾筒調(diào)高位移響應(yīng)(vq=2 m/s)
Fig.8Drumtransferresponsewith2m/sanddifferentcuttingimpedance
為驗證AMESim仿真模型得到結(jié)果的精度,使用MATLAB對采煤機調(diào)高機構(gòu)動態(tài)特性響應(yīng)進行數(shù)值計算。由于PID控制器作用下,調(diào)高機構(gòu)動態(tài)響應(yīng)誤差更加明顯,因此在此僅列出PID控制器作用下,調(diào)高誤差極值的數(shù)據(jù)。數(shù)值計算得到結(jié)果與AMESim仿真模型得到結(jié)果對比見表1、表2。
(a)vq=2 m/s
(b)vq=3 m/s
(c)vq=4 m/s
圖9不同牽引速度下的滾筒調(diào)高位移響應(yīng)(A=200 kN/m)
Fig.9Drumtransferresponsewithdifferenttractionspeedand200kN/mcuttingimpedance
圖10A=200kN/m、vq=3m/s滾筒調(diào)高位移響應(yīng)
Fig.10Drumtransferresponsewith2m/stractionspeedand200kN/mcuttingimpedance
由表1、表2中的數(shù)據(jù)可知,數(shù)值計算結(jié)果略大于仿真結(jié)果果,誤差最大值為12.89%,小于15%,滿足精度要求。
表1 輸入方波信號的仿真結(jié)果與數(shù)值計算結(jié)果
表2 輸入正弦信號的仿真結(jié)果與數(shù)值計算結(jié)果
(a)A=100 kN/m
(b)A=200 kN/m
(c)A=300 kN/m圖11 不同截割阻抗情況下的滾筒調(diào)高位移響應(yīng)誤差(vq=2 m/s)Fig.11 Drum transfer response error with 2 m/s and different cutting impedance
(a)vq=2 m/s
(b)vq=3 m/s
(c)vq=4 m/s圖12 不同牽引速度下的滾筒調(diào)高位移響應(yīng)誤差(A=200 kN/m)Fig.12 Drum transfer response error with different traction speed and 200 kN/m cutting impedance
使用常規(guī)PID控制器對調(diào)高機構(gòu)進行控制時,調(diào)高機構(gòu)輸出的滾筒位置響應(yīng)對設(shè)定曲線的跟蹤性能較差,誤差較大,其控制效果受截割阻抗和牽引速度的變化影響較大。仿真結(jié)果表明,隨著截割阻抗以及牽引速度的增大,滾筒的截割載荷增大,因此調(diào)高機構(gòu)的擾動負載增大,使得控制精度下降,這說明,由于常規(guī)PID控制器的參數(shù)整定是針對某一特定工況進行的,因此常規(guī)PID控制器對不同工況的適應(yīng)性較弱。
(1)不同牽引速度和不同煤巖截割阻抗會影響滾筒截割負載,從而影響調(diào)高機構(gòu)的控制效果。
(2)使用模糊PID控制器對調(diào)高機構(gòu)進行控制相比常規(guī)PID,具有更好的跟蹤性能,更低的誤差。
(3)常規(guī)PID控制器控制效果受截割阻抗和牽引速度的變化影響較大,模糊PID控制器控制效果受截割阻抗和牽引速度的變化影響相對較小。
[1] 劉春生. 滾筒式采煤機理論設(shè)計基礎(chǔ)[M]. 徐州:中國礦業(yè)大學(xué)出版社,2003.
LIU Chunsheng. Theoretical Basis of Roller Shearer [M]. Xuzhou: China University of Mining and Technology Press, 2003.
[2] 蘇秀平,朱樓坤,李威. 采煤機調(diào)高系統(tǒng)模擬加載實驗臺的研制[J].機械設(shè)計, 2015, 32(2): 83-86.
SU Xiuping, ZHU Loukun, LI Wei. Development of Loading Experiment Table for Shearer-drum-height Mechanism[J]. Journal of Machine Design, 2015, 32(2): 83-86.
[3] 蘇秀平.采煤機自動調(diào)高控制及其關(guān)鍵技術(shù)研究[D].徐州:中國礦業(yè)大學(xué),2013.
SU Xiuping. Study on Key Technologies of Auto-height Adjustment for Shearer[D]. Xuzhou:China University of Mining and Technology,2013.
[4] 蘇秀平,李威,樊啟高.采煤機滾筒調(diào)高滑模變結(jié)構(gòu)控制策略[J]. 煤炭學(xué)報, 2012, 37(12): 2107-2111.
SU Xiuping, LI Wei, FAN Qigao. A Shearer Drum Height Adjusting Strategy Using Sliding-mode Variable Structure Control[J]. Journal of China Coal Society, 2012, 37(12): 2107-2111.
[5] 劉鵬,孟磊,王勃,等. 基于位姿測量與煤層DEM的采煤機滾筒自動調(diào)高方法[J]. 煤炭學(xué)報, 2015, 40(2): 470-475.
LIU Peng, MENG Lei, WANG Bo, et al. An Automatic Height Adjustment Method for Shearer Drums Based on Pose Measurement and Coal Seam DEM[J]. Journal of China Coal Society, 2015, 40(2): 470-475.
[6] 周信,王忠賓,譚超,等. 基于雙坐標系的采煤機截割路徑平整性控制方法[J]. 煤炭學(xué)報, 2014, 39(3):574-579.
ZHOU Xin, WANG Zhongbin, TAN Chao, et al. A Smoothness Controlling Method for the Cutting Path of the Shearer Based on the Double-coordinators[J]. Journal of China Coal Society, 2014, 39(3): 574-579.
[7] 張春雨.采煤機滾筒調(diào)高電液比例控制系統(tǒng)的研究[D]. 太原:太原理工大學(xué),2014.
ZHANG Chunyu. Research on Shearer Roller's Position Adjustment System That Based on Electro-hydraulic Proportional Control[D]. Taiyuan:Taiyuan University of Technology,2014.
[8] 張俊梅,范迅,趙雪松. 采煤機自動調(diào)高控制系統(tǒng)研究[J]. 中國礦業(yè)大學(xué)學(xué)報, 2002, 31(4): 415-418.
ZHANG Junmei, FAN Xun, ZHAO Xuesong. Automatic Horizon Control System of Coal Mining Machine[J]. Journal of China University of Mining & Technology, 2002, 31(4): 415-418.
[9] YAO B, BU F, CHID G T C. Nonlinear Adaptive Robust Control of Electro-hydraulic Servo Systems with Discontinuous Projections[C]//Proceedings of IEEE Conf. Decision and Control. Tampa, 1998: 265-270.
[10] 劉春生, 楊秋, 李春華. 采煤機滾筒記憶程控截割的模糊控制系統(tǒng)仿真[J]. 煤炭學(xué)報, 2008, 33(7): 822-825.
LIU Chunsheng, YANG Qiu, LI Chunhua. Simulation of Shearer Drum Cutting with Memory Program Controlling by Fuzzy Control[J]. Journal of China Coal Society, 2008, 33(7): 822-825.
[11] 曾慶良,張海忠,王成龍,等. 采煤機調(diào)高系統(tǒng)的機液協(xié)同仿真分析[J]. 煤炭科學(xué)技術(shù), 2015, 43(1): 86-90.
ZENG Qingliang, ZHANG Haizhong, WANG Chenglong, et al. Analysis on Mechanic and Hydraulic Collaborated Simulation of Height Adjusting System in Coal Shearer [J]. Coal Science and Technology, 2015, 43(1): 86-90.
[12] FAN Qigao, LI Wei, WANG Yuqiao, et al. Control Strategy for an Intelligent Shearer Height Adjusting System[J]. Mining Science and Technology, 2010, 20(6): 908-912.
[13] 張修榮. 采煤機滾筒自適應(yīng)液壓調(diào)高系統(tǒng)研究[D]. 西安:西安科技大學(xué),2009.
ZHANG Xiurong. Research of Adaptive Hydraulic Shearer Drum System[D]. Xi’an: Xi’an University of Science and Technology, 2009.
[14] 王忠賓,徐志鵬,董曉軍. 基于人工免疫和記憶切割的采煤機滾筒自適應(yīng)調(diào)高[J]. 煤炭學(xué)報, 2009, 34(10): 1405-1409.
WANG Zhongbin, XU Zhipeng, DONG Xiaojun. Self-adaptive Adjustment Height of the Drum in the Shearer Based on Artificial Immune and Memory Cutting[J]. Journal of China Coal Society, 2009, 34(10): 1405-1409.