仲霜霜, 王琪, 方海峰, 王牧笛
(1.江蘇科技大學機械工程學院,江蘇鎮江 212003;2.江蘇貝爾機械有限公司,江蘇張家港 215600)
我國是世界上農業大國,在農地膜生產和使用上均居世界前列。塑料地膜能夠加快農作物的生長、提高農作物的產量,給農業發展帶來了巨大的經濟效益,但隨之也出現了一系列嚴重的環境污染和生態問題,如何處理覆蓋后的殘膜已經成為人類最關注的問題之一[1]。對農業中所廢棄的覆蓋地膜、大棚種植中所用的塑料薄膜以及日常生活中所用到的塑料袋等等塑料薄膜的回收利用,都離不開撕碎機的撕碎處理。廢舊農地膜撕碎機是一種廢舊農地膜撕碎及砂石秸稈的破碎分離裝置,能夠對廢舊農地膜進行撕碎及細小垃圾的清除,因其破碎能力好,產量高等優點被應用于垃圾預處理階段,對廢膜的回收利用起著重要的作用。滾筒是農地膜撕碎機的主要部件之一,由于在整個工作過程中,滾筒在電動機的驅動下高速轉動,其穩定性是影響整個撕碎機穩定性的重要因素之一,因此,對滾筒進行模態分析,研究其工作時是否容易發生共振現象,對于提高廢舊農地膜撕碎機運行時的穩定性具有一定的意義。
撕碎機的進料倉下部與喂料螺旋連接,喂料螺旋水平安裝在主軸上,主軸由電動機驅動進行高速運轉,螺旋喂料器通過減速機實現和主軸的同軸差速。將廢舊農地膜投入進料倉后,由喂料螺旋將物料輸送到滾筒與外筒的腔體內,螺旋做間歇性的運動,同時還具備自動反轉功能。如果物料中含有大尺寸堅硬的垃圾,當被放入進料倉時,螺旋輸送機會自動反轉并使其進入儲料裝備,保證喂料均勻的同時避免了堵料、卡料。
滾筒由電動機驅動進行高速運轉,滾筒上裝有依照不同螺旋線形式排布的刀體,撕碎機的外筒套接在主軸的外部,外筒內壁上裝有同樣螺旋線排列的刀體,外筒隨喂料螺旋運動狀態做間歇性運動,針對不同濕度和不同結節狀態的物料配備變頻調速功能。筒壁上以螺旋形式排列刀體。刀體末端使用合金刀頭,滾筒上刀體與外筒上刀體之間相互交錯排列,形成對切,能夠配合滾筒上的刀體對進入腔體中的廢膜進行擠壓、剪切、撕碎等作用。以確保物料的徹底撕碎、打散及快速排出。外筒下方配備偏心輪,在轉動的同時做偏心運動,能夠將分離出的泥土、砂石等雜質更徹底地通過筒壁網孔篩出。滾筒及外筒筒身呈錐形形狀,符合廢舊農地膜從進料到出料的過程中體積逐漸變化的過程。筒壁所有刀體可快速自由拆卸,方便客戶在使用過程中快速更換主軸刀體。
被撕碎后的碎膜,將會從出料口被排出,然后落在下面的盤篩上,在盤篩的不停振動作用下,大的農膜被分選出來,進行回收利用,而細小的垃圾雜質將掉落在輸送皮帶上,被運送出去,皮帶采用 V型托輥設計,能夠避免雜質、碎膜的外漏。從而完成整個加工材料的撕碎及雜質的分離過程。其結構如圖1所示。

圖1 廢舊農地膜撕碎機結構圖
首先在SolidWorks中建立滾筒的三維模型。由于滾筒的結構比較復雜,為了能夠簡化計算,在對分析精度不造成影響的前提下,對結構進行適當的簡化,忽略對結果影響不大的小孔、過渡圓角和螺紋孔,簡化和去除一些細小零件[2]。由于零件在焊接時會留有焊縫,因此在建立結構模型時,將各個焊接件作為整體零件來建模。其實體模型如圖2所示。

圖2 滾筒模型
將SolidWorks中處理好的三維滾筒模型通過Parasolid_txt格式導入到ANSYS前處理器中??紤]到滾筒是三維實體結構,對于三維問題,最常用的是實體單元。SOLID187是三維10節點四面體結構單元,這種高階實體單元適用于CAD/CAM中結構復雜、形狀不規則的模型,因此選用該單元。滾筒的材料為Q235鋼,設置其密度為7.85 g/cm3、彈性模量為200 GPa、泊松比為0.3。對滾筒進行自由網格劃分,劃分成四面體單元網格形式,劃分完的完的網格如圖3所示,劃分完后的模型節點總數為72 056,單元數總為34 799。
在滾筒高速旋轉的狀態下,滾筒會產生一定的旋轉預應力,預應力的存在會使滾筒系統的結構剛度產生變化,從而影響滾筒的固有頻率[3-4]。所以對滾筒進行模態分析時,將預應力考慮在內具有實際意義。滾筒工作時的轉速為2000 r/min,計算得其角速度為209 rad/s,設置系統角速度,對其施加速度載荷。根據撕碎機的實際裝配情況,對滾筒內孔徑兩端100 mm處施加徑向零位移約束,對滾筒左右兩個端面施加軸向零位移約束。

圖3 有限元模型
目前,作為一門重要的工程技術,模態分析技術在我國的各個工程領域都得到了廣泛的使用。通過模態分析可以得出結構各階的固有頻率和振動特性,預知結構是否會發生共振現象[5],從而在產品試制之前進行改進、優化,改善產品結構。所以模態分析的應用如今成為動力學分析中必不可少的方法之一[6]。
模態分析理論在技術基礎上,與CAE技術相結合,在理論基礎上,與有限元方法、結構動力學理論相結合[7]。由彈性力學有限元法的理論知識可知N個自由度的系統線性振動微分方程為

式中:F(t)為受到的外界激勵力向量;C為阻尼矩陣;M為質量矩陣;δ¨為加速度向量;K為剛度矩陣;δ˙為速度向量;δ為位移向量。
在對該方程進行求解時,因為模型的固有頻率和振型與外界載荷沒有關系,所以F(t)=0。
由實驗經驗可知,結構的固有頻率基本上不受到阻尼的影響,因此可以不做討論,即C=0,從而得到無阻尼自由振動的運動微分方程為

式中:δ為模型的振型;ω為模型的固有頻率。
只有當稀疏行列式等于零時,此方程才能夠得到非零解,即

對于一個系統來說,求解其固有頻率和振型,其實就是求解矩陣的特征值ω和δ特征向量的值。
運用ANSYS中的模態計算功能,進行滾筒模態頻率的提取。分析中選用Block Lanczos法,一般來說,機械結構的振動是其各階固有振型的線性組合[8],其中高階的振型對機械結構的動力性能影響相對較小,且不容易被激勵,會引起一定的計算誤差[9],所以在分析時只對滾筒的前八階模態進行提取,得到結果如表2所示,其相應的振型如圖4所示。
結合振型動畫與圖4可知每一階固有頻率都有其相應的振型,可看出各部位振動的強弱分布及其在相應模態下的振動響應,滾筒的前8階固有頻率中,除了第一階和第二階固有頻率外,其余相鄰兩階固有頻率之間的值非常相近,并且相互接近的固有頻率對應的振型只在振動方向上有所不同,振型形狀十分相似。滾筒的振動主要是滾筒中間部位的內外振動,類似于膨脹、收縮的運動。由此可見,振動的危險點主要集中在滾筒的中間部位,并且隨著階數的增大,危險點的數目也逐漸增加,滾筒中間的振動就愈加劇烈。

表2 滾筒的固有頻率及振型描述

圖4 滾筒前八階振型
建立了滾筒的三維實體模型和有限元模型,并根據實際工況設置了相應的約束,對滾筒進行了有預應力的模態分析,計算了前八階的固有頻率和振型,結果表明除前兩階固有頻率外,滾筒的振動主要是滾筒中間部位的內外振動,并且隨著階數的增加,滾筒中間的振動就愈加劇烈。滾筒工作頻率約為0~150 Hz,而固有頻率較高,主要集中在200~600 Hz,大于其工作時的頻率,所以在滾筒工作時不會發生共振現象,滿足初步設計要求。該研究為滾筒做諧響應分析、瞬態動力學分析及其他動力學分析提供了依據。
[參考文獻]
[1] 胡宏亮.生物降解地膜的產量效應和降解特性及大田示范研究[D].杭州:浙江大學,2015.
[2] 劉清龍.數控機床結合面特性參數的識別研究及其主結構的動態優化[D].北京:北京工業大學,2007.
[3] 李嘯天,韓振南.基于ANSYS軟件的轉子系統臨界轉速及模態分析[J].機械管理開發,2010(3):184-185.
[4] 王犇,華林.高速旋轉狀態下汽車弧齒錐齒輪的動力學模態分析[J].汽車工程,2011(5):447-451.
[5] 王勇勤,劉峻亦,嚴興春,等.12.5MN雙柱下拉式快鍛液壓機整機模態分析[J].機械設計與制造,2013(2):166-168.
[6] 蘇勛文,王少萍,朱冬梅,等.直升機中減速器諧響應分析與傳感器優化布局[J].北京航空航天大學學報,2011(9):1049-1053.
[7] 梁君,趙登峰.模態分析方法綜述[J].現代制造工程,2006(8):139-141.
[8] 張憲,何洋,鐘江,等.疲勞振動試驗臺的模態與諧響應分析[J].機械設計與制造,2008(4):12-14.
[9] 鄭剛,謝方偉,王和偉,等.盤式制動器關鍵零部件模態分析[J].機械設計與制造,2013(2):172-173.