張義民, 張 睿, 朱麗莎, 趙春雨
(東北大學 機械工程與自動化學院,遼寧 沈陽 110819)
采煤機是煤礦綜采的關鍵裝備,搖臂作為采煤機截割荷載傳遞及傳動件支承的主要部件,受到截割荷載及搖臂傳動件內激勵的作用,是采煤機可靠性的主要薄弱環節。搖臂的振動和變形直接影響傳動件受載及可靠性,有必要了解其動態特性。國外對采煤機搖臂研究較少,國內主要通過數學建?;蛴邢拊浖?,利用動力學軟件分析,存在缺乏實驗驗證的問題。具體表現在固有特性缺乏實驗驗證、截割荷載利用近似公式模擬、仿真模型精確度不高等[1-11]。
依托國家973課題“深部危險煤層無人采掘裝備關鍵基礎研究”,以MG500/1180-WD滾筒采煤機為研究對象,通過有限元模型對搖臂固有特性進行分析,以固有頻率為判定指標,用實驗數據對有限元模型進行驗證。重點研究了搖臂固有頻率隨設計變量的變化規律。根據傳動系統嚙合頻率對搖臂進行諧響應分析,得到搖臂應力分布及薄弱環節。
本文以MG500/1180-WD型采煤機的搖臂為研究對象,該采煤機搖臂由一個三級直齒輪減速系統和兩級行星減速系統以及殼體組成。利用Pro/E建立MG500/1180-WD搖臂實體模型并導入ANSYS。實體建模忽略微小倒角、圓角,采用自由網格劃分使網格適應搖臂的復雜結構。圖1為搖臂有限元模型。
有限元模型單元為Solid185,材料為ZG25MnNi,屈服強度270 MPa,密度7 800 kg/m3,將材料參數、單元參數輸入ANSYS,并將搖臂上下鉸耳與牽引部接觸面設置固定約束,提取前4階模態振型。ANSYS計算的前8階固有頻率如表1所示。圖2為搖臂固有振型。1階振型為行星頭沿Z軸彎曲;2階振型為行星頭沿Y軸彎曲;3階振型為搖臂中部和行星頭沿Y軸彎曲;4階振型為電機殼體繞Y軸扭轉。
階次12345678頻率/Hz368528685801963104512361508
(a) 1階
(b) 2階
(c) 3階
(d) 4階
由模態分析可知,搖臂主要振動模式為電機殼體及行星頭彎曲扭轉振動。電機殼體和行星頭是搖臂的薄弱環節,相對振動大。
搖臂固有頻率取決于剛度及質量,搖臂質量和剛度與壁厚及肋板尺寸密切相關。選取5個設計變量分析固有頻率隨設計變量的變化趨勢,即:搖臂中部壁厚(老塘側),搖臂中部壁厚(煤壁側),電機殼體壁厚,行星頭壁厚,肋板高度。通過修改實體模型尺寸并在ANSYS中計算每一設計尺寸對應的固有頻率,共調用有限元模型165次。設計變量對各階固有頻率影響規律如圖3所示,圖中曲線從下至上依次為1~8階固有頻率。
(a) 肋板高度
(b) 電機殼體壁厚
(c) 老塘側壁厚
(d) 煤壁側壁厚
(e) 行星頭壁厚
(1) 觀察圖3,設計尺寸對搖臂固有頻率影響程度不同。老塘側壁厚對固有頻率影響最大,肋板高度和電機殼體厚度次之,煤壁側壁厚及行星頭壁厚變化對搖臂固有頻率影響很小。
(2) 觀察圖3(c),在老塘側壁厚87 mm處,除f7外,其他固有頻率數值上均發生變化。其中f1~f4及f8增加,f5~f6減小,2階固有頻率改變最大,達到141.73 Hz。說明老塘側壁厚變化對除f7之外的所有固有頻率均有明顯影響,改變老塘側壁厚可以明顯改變搖臂固有頻率。
(3) 觀察圖3(a)、圖3(b),隨著肋板高度和電機殼體壁厚的增加,f3,f7,f8逐漸增加,固有頻率變化值最大為50.34 Hz,遠小于老塘側最大變化值141.73 Hz。其余固有頻率變化值在5 Hz以內,說明肋板高度和電機殼體壁厚分段影響搖臂固有頻率,但影響程度遠小于老塘側壁厚。
依托張家口國家能源采掘裝備研發實驗中心,進行采煤機搖臂模態實驗,獲取搖臂固有頻率,驗證有限元模態分析結果。
合理選取激振點及拾振點,測得激振點與拾振點的振動加速度,運用數字信號處理技術得到頻響函數,進而識別搖臂模態參數。
根據頻響函數的定義,頻響函數矩陣為
(1)
式中:H(ω)為頻響函數矩陣,頻響函數矩陣表達式為
(2)
實驗模態分析還可通過功率譜密度得到頻響函數,表達式為
(3)
式中:Gxx(ω)為激勵自功率譜密度;Gxy(ω)為互功率譜密度。
頻響函數矩陣任意一行(列)均包含系統所有模態參數,當激勵和響應相關度較好時,得到的數據可信度高,對綜合頻響函數利用峰值法即可得到固有頻率[12-13]。
實驗平臺依托國家能源采掘裝備研發實驗中心綜采工作面,實驗平臺組成如圖4所示。
(a)實驗設備(b)實驗現場
圖4 采煤機實驗平臺
Fig.4 Shearer experiment platform
實驗平臺由采煤機、加速度傳感器、數據采集器、計算機、力錘及連接線組成(如圖4所示)。傳感器分為三向及單向傳感器,左右搖臂各兩個三向傳感器及兩個單向傳感器。
采用單點激振、多點拾振脈沖(Single Input/Multi-Output,SIMO)激勵法進行模態試驗。拾振點即傳感器安裝點,傳感器安裝位置如圖5所示。傳感器安裝位置包括直齒輪系統、惰輪系統及行星輪系統,選取靠近軸承座的位置安裝。
圖5 傳感器安裝位置
激振點選在滾筒截齒??紤]采煤機工作時截齒所受荷載方向,將激振方向選為截齒所受三向力方向。為增加信噪比,每個方向敲擊3次,共進行35次測試。滾筒敲擊位置如圖6所示。
圖6 滾筒敲擊位置
搖臂的約束采用原裝約束,邊界條件與工作狀態一致。將得到的激勵、響應數據通過模態測試軟件進行SIMO分析。
得到搖臂綜合頻響函數如圖7所示。對系統綜合頻響函數進行峰值分析得到低階固有頻率如表2所示。
圖7 頻響函數曲線
階次12345678頻率/Hz350525675828976112512751423
將結果與理論計算結果對比,前8階固有頻率相對誤差在10%以內,說明理論計算結果是合理的。
諧響應分析的載荷包括幅值及頻率范圍。載荷幅值通過重載截割試驗及理論分析獲得;載荷頻率范圍由搖臂傳動系統嚙合頻率計算得出。
重載截割試驗煤壁以煤粉為本體,以水泥、水及一些特骨料為粘合劑,按一定時間間隔分層澆筑而成,煤壁模擬出了煤的層理和節理。煤巖體主要參數:硬度f3,密度1 359 kg/m3,彈性模量259.23 MPa,抗壓強度2.33 MPa。采煤機割煤工況:割煤截深600 mm,啟動后將牽引速度由0增加至3 m/min截割煤壁2 min后停機,重載截割試驗現場如圖8所示。
采煤機搖臂諧響應載荷幅值的確定方法,目前沒有相關文獻報道,其難點主要在于采煤機搖臂動力學模型復雜,截割載荷是隨機載荷,載荷及動力學模型需要進行理論假設和簡化,求解出的振動響應準確度得不到保證。針對這種情況,本文采用重載截割實驗振動數據,直接得到搖臂軸承座附近振動加速度,運用數值積分方法得到振動速度;同時對理論模型進行求解,以測點處三向振動速度均方根值(Root Mean Square,RMS)為判定指標,將實驗結果與理論計算進行對比并設定允差ε(本文設定為10%),當RMS值滿足
(4)
表明理論計算與實驗測試RMS值吻合較好,則計算結束,否則應修改截割隨機荷載模型,直到得到的RMS值接近為止。圖9為重載截割試驗數據(篇幅所限,僅列出行星級X向)。將振動加速度進行數值積分得到振動速度并求出RMS值;以項目組前期對采煤機截割荷載的研究成果為基礎,將單齒截割荷載視作頻域Gamma分布并編制載荷譜,通過滾筒受力分析將截齒三向力合成得到滾筒三向力及力矩,代入項目組所建立的搖臂有限元模型中得到理論計算振動速度并求出RMS值。各測點振動速度RMS值統計表如表3所示。
最大允差εmax為6.74%,在合理區間內。通過計算得到軸承座內圈接觸剛度,與接觸面三向相對位移代入搖臂有限元模型,最終得到搖臂諧響應載荷幅值為140 kN。
表3 各測點振動速度RMS值統計表
圖8 重載截割試驗現場
圖9 重載截割測試數據
根據MG500/1180-WD型滾筒采煤機傳動系統參數,利用直齒輪及行星輪傳動嚙合頻率計算公式,得到采煤機搖臂齒輪傳動系統嚙合頻率如表4所示。
表4 傳動系統嚙合頻率
得到載荷的幅值和頻率范圍后,對搖臂進行諧響應分析。根據傳動系統嚙合頻率及其倍頻的取值范圍,諧響應激勵頻率范圍確定為200~1 400 Hz,在搖臂傳動系統軸承座孔處施加激勵,選取200個頻率點進行分析[14-15]。
圖10為頻域響應。由圖可知,激勵頻率為325 Hz、716 Hz及1 380 Hz時,搖臂動態響應較大。進一步分析可知,325 Hz為行星1級2倍頻,716 Hz為直齒1級2倍頻附近,1 380 Hz為直齒2級2倍頻及8階固頻附近。利用ANSYS求得激振頻率為325 Hz、716 Hz及1 380 Hz時搖臂響應。圖11為應力云圖,從圖中可知,搖臂應力在716 Hz時最大,應力值30.9 MPa;應力最大處為電機殼體與搖臂中部交界處;搖臂中部模態振動及諧響應應力均遠小于行星頭及電機殼體,其中應力最大值20.5 MPa。
(a) 電機殼體
(b) 搖臂中部
(c) 行星頭
(a) 325 Hz應力云圖
(b) 716 Hz應力云圖
(c) 1 380 Hz應力云圖
由諧響應分析結果可知,325 Hz、716 Hz及1 380 Hz是激勵頻率中容易導致搖臂振動響應偏大的成分,從應力云圖看,1 380 Hz時響應明顯小于325 Hz、716 Hz,因此與325 Hz和716 Hz接近的1階、3階固有振型易被激發。據此可判斷搖臂的主要振動形態為1階、3階振型模態振動為主要特征的彈性振動,即搖臂行星頭及中部的彎曲振動。采煤機搖臂可視作鉸接在牽引部上的懸臂梁結構,因此行星頭及搖臂中部的彎曲振動形態會導致搖臂與牽引部鉸接處徑向以及搖臂中部豎直方向振動響應偏大。將理論分析結果與振動試驗測得的結果進行比較(如圖12所示)可發現,搖臂中部豎直方向(圖中X向)振動加速度大于Y、Z向;鉸接點處徑向(圖中X、Y向)振動加速度大于Z向(軸向)。實驗結果與理論推斷吻合良好。
(a) 搖臂中部
(b) 鉸接點
(1) 通過有限元模態分析及實驗模態分析,獲得MG500/1180-WD采煤機搖臂固有特性,得到前8階固有頻率及振型。
(2) 通過分析固有頻率隨設計變量的變化規律,得出搖臂中部老塘側壁厚是影響搖臂固有特性的關鍵因素,肋板高度和電機殼體壁厚分段影響搖臂固有頻率,但影響程度遠小于老塘側壁厚,煤壁側壁厚及行星頭壁厚的變化對搖臂固有頻率幾乎沒有影響。
(3) 利用ANSYS對MG500/1180-WD搖臂進行諧響應分析,結果表明:搖臂在行星1級2倍嚙頻、直齒1級2倍嚙頻、直齒2級2倍嚙頻及8階固頻激勵下動態應力較大。
(4) 采煤機搖臂形成了以1階、3階振型模態振動為主要特征的彈性振動并導致搖臂鉸接耳處產生應力集中。分析結果可為采煤機搖臂結構優化設計及可靠性分析提供參考。
參 考 文 獻
[1] ZAGHBANI I,SONGMENE V.Estimation of machine-tool dynamic parameters during machining operation through operational modal analysis[J].International Journal of Machine Tools and Manufacture,2009,49(12/13):947-957.
[2] 劉長釗,秦大同,廖映華.采煤機截割部機電傳動系統動力學特性分析[J].機械工程學報,2016,52(7):15-22.
LIU Changzhao,QIN Datong, LIAO Yinghua. Dynamic characteristics analysis of mechanical and electrical transmission system of cutting part of shear[J].Journal of Mechanical Engineering,2016,52(7):15-22.
[3] 趙麗娟,田震.薄煤層采煤機振動特性研究[J].振動與沖擊,2015,34(1):195-199.
ZHAO Lijuan,TIAN Zhen.Vibration characteristics of thin coal seam shearer[J].Journal of Vibration and Shock,2015,34(1):195-199.
[4] 申寶宏,郭玉輝.我國綜合機械化采煤技術裝備發展現狀與趨勢[J]. 煤炭科學技術,2012,40(2):1-3.
SHEN Baohong,GUO Yuhui. Development status and tendency of technology and equipment for fully mechanized coal mining in China[J]. Coal Science and Technology,2012,40(2):1-3.
[5] 羅晨旭.滾筒采煤機開采含煤巖界面煤層截割特性研究[D]. 徐州:中國礦業大學,2015.
[6] 趙麗娟,馬連偉.薄煤層采煤機可靠性分析與疲勞壽命預測[J]. 煤炭學報,2013,38(7):1288-1292.
ZHAO Lijuan,MA Lianwei.Thin seam shearer reliability analysis and fatigue life prediction[J].Journal of China Coal Society,2013,38(7):1288-1292.
[7] 張義民,黃婧,朱麗莎,等. 采煤機搖臂傳動系統可靠性穩健優化設計[J].煤炭學報,2015,40(11):2540-2545.
ZHANG Yimin,HUANG Jing,ZHU Lisha.Robustreliability improvement design of rocker arm drive system of coal mining machine[J].Journal of China Coal Society,2015,40(11):2540-2545.
[8] 張世洪.我國綜采采煤機技術的創新研究[J]. 煤炭學報,2010,35(11):1898-1902.
ZHANG Shihong.Study on the innovation of fully mechanized coal shearer technology in China[J].Journal of China Coal Society,2010,35(11):1898-1902.
[9] 李曉豁,李婷,焦麗,等.滾筒采煤機截割載荷的模擬系統開發及其模擬[J].煤炭學報,2016,41(2):503-507.
LI Xiaohuo,LI Ting,JIAO Li. Development of cutting load simulation system and its simulation study on drum shearers[J].Journal of China Coal Society,2016,41(2):503-507.
[10] 趙麗娟,劉旭楠,馬聯偉.基于經濟截割的采煤機運動學參數優化研究[J]. 煤炭學報,2013,38(8):1490-1495.
ZHAO Lijuan,LIU Xu’nan,MA Lianwei. Optimization research on shearer’s kinematic parameters based on economical cutting[J].Journal of China Coal Society,2013,38(8):1490-1495.
[11] 楊陽,馬鵬程.采煤機差動2K-H行星輪系固有特性分析[J].煤炭學報,2016,41(2):509-512.
YANG Yang,MA Pengcheng. Analysis of shearer 2K-H differential planetary gear train inherent characteristics[J].Journal of China Coal Society,2016,41(2):509-512.
[12] 肖正明,秦大同,武文輝.土壓平衡盾構機行星減速器動態特性測試與分析[J].振動與沖擊,2012,31(9):156-159.
XIAO Zhengming,QIN Datong,WU Wenhui, et al. Test and analysis of the dynamic performance of the earth pressure balance shield machine[J].Journal of Vibration and Shock,2012,31(9):156-159.
[13] 張力.機械振動實驗與分析[M].北京:清華大學出版社,2013:78-86.
[14] 常樂浩,賀朝霞,劉更,等.集平行軸齒輪傳動系統動力學建模的有限單元法[J].振動與沖擊,2016,35(20):48-53.
CHANG Lehao,HE Zhaoxia,LIU Geng, et al. Dynamic modeling of parallel shaft gear transmissions using finite element method[J].Journal of Vibration and Shock,2016,35(20):48-53.
[15] 向玲,賈秩,李媛媛,等.內外激勵作用下多自由度齒輪系統的非線性動力學特性[J].振動與沖擊,2016,35(13):154-159.
XIANG Ling,JIA Zhi,LI Yuanyuan. Non-linear dynamic features of a gear system with multi-DOF subjected to internal and external excitation[J].Journal of Vibration and Shock,2016,35(13):154-159.