許梁 呂佩師 李文偉
青島海爾滾筒洗衣機有限公司 山東青島 266101
隨著洗衣機向著大容量(幾何尺寸變大,剛度變低)、高轉速及復合式分區洗筒結構的趨勢快速發展,其工作頻率越來越接近于洗衣機的共振頻率(臨界轉速),特別是隨著家電行業競爭越來越激烈,爭相對復合式分區洗筒結構進行創新研究開發,造成復合式筒部件質量加大、共振頻率進一步降低,一旦洗衣機最高轉速接近或達到其共振頻率,內外筒之間的跳動變形會猛然增加,嚴重的會出現內外筒碰磨,出現嚴重振動噪音問題及結構安全失效。因此洗衣機尤其是復合式滾筒的避共振設計變得尤為重要。
洗衣機共振頻率計算,傳統的方式是進行有限元模態分析(Normal Mode),或通過模態錘擊、激振實驗對物理樣機進行頻率實測,但上述方法的不足之處是無法考慮內筒的轉動影響,從而無法考慮由于轉動影響引起的動態陀螺效應對洗衣機共振頻率的改變。
滾筒洗衣機筒總成由內外筒組成,內筒通過軸與軸承與外筒進行連接,內筒可視為高速旋轉的柔性轉子。對于滾筒洗衣機筒部件轉子系統,其內筒會受到陀螺力矩作用,因而需要研究陀螺效應對其的影響。隨著滾筒洗衣機內筒容量的擴大及轉速的提升,作用在內筒及軸上的陀螺力矩更加顯著。對于滾筒洗衣機筒部件來說,內筒轉子的進動頻率隨其自轉轉速變化,特別是當內筒在滿負載狀態(轉動慣量大)及額定速度旋轉時,陀螺力矩作用不可以忽略。若在內外筒系統動力學設計和分析階段忽略陀螺效應的影響,將會導致分析結果嚴重失真,造成嚴重的安全及可靠性問題。
本文將以某款滾筒洗衣機為例,研究考慮陀螺效應對滾筒洗衣機系統臨界轉速的影響,討論以多柔體動力學分析方法進行系統臨界轉速的計算。
對于一般機械系統,通過有限元法可以得到其動力學方程:

式中:[M]為系統質量矩陣;[C]為系統阻尼矩陣;[K]為系統剛度矩陣;{u&&}、{u&}和{u}分別為系統位移、速度和加速度向量;{F}為系統外部激振力向量。
對于轉子系統,考慮陀螺效應后,轉子系統動力學方程變為[4]:

式中:Ω為軸的轉速;[Jp]和[G]分別為極轉動慣量矩陣和陀螺矩陣,為反對稱陣。
陀螺力矩的方向垂直于進動平面,當陀螺力矩為正(正進動)時,會使得自轉軸轉向進動軸,相當于增大了軸系的抗彎剛度,使轉子的臨界轉速提高;當陀螺力矩為負(反進動)時,會使得自轉軸偏離進動軸,相當于降低了軸系的抗彎剛度,使轉子的臨界轉速減小[5]。
文獻[6]在考慮陀螺力矩對軸臨界轉速的影響的前提下,通過影響系數建立了方程,并給出了臨界轉速的解。
對于轉子動力學有限元分析,“陀螺效應”一直是制約其發展的“瓶頸”問題[3],但多篇文獻提到ANSYS有限元分析軟件在進行轉子系統的臨界轉速計算時,可以考慮 “陀螺效應”對臨界轉速的影響。文獻[1],[3],[7]論述了使用ANSYS軟件,對于已建立的有限元模型,在進行模態分析時對軸施加不同的轉速,采用QR Damp模態求解方法計算轉子在各轉速下的各階固有頻率。
本文從多柔體動力學仿真分析入手,計算滾筒洗衣機的臨界轉速,并與不考慮陀螺效應影響的模態分析結果進行比較,考察陀螺效應對滾筒洗衣機系統臨界轉速的影響,并由此提供一種求解滾筒洗衣機系統臨界轉速的計算方法。
本文首先進行不考慮陀螺效應影響的滾筒洗衣機自由模態分析(懸掛約束近似自由),由于質量是關鍵的動力學參數之一,故分析分為空筒和滿載兩種工況,圖1(a)、(b)分別為空筒和滿載的驅動部有限元模型,洗衣機配重以集中質量代替。
略去前六階剛體模態,驅動部的前兩階彈性模態振型均表現為內、外筒的相對彎曲擺動,兩階模態的彎曲擺動方向大致垂直,由于外筒底在圓周方向剛度有差異,所以會出現頻率相近但不相等的兩階模態,本文僅列出第一階彈性模態。圖2、3分別為空筒及滿載工況下,驅動部第一階模態振型圖,其固有頻率分別為57.4Hz,40.3Hz,二者振型一致,但固有頻率差別明顯,可見內筒質量對驅動部固有頻率影響顯著。
以上為不考慮陀螺效應影響的有限元模態分析,為計及陀螺效應,現進行驅動部(內外筒系統)的柔體動力學分析,圖4為建立的同一型號洗衣機的剛柔混合動力學模型,同樣考慮空筒和滿載,分三種工況:(1)空筒&小偏心(規格A);(2)空筒&中等偏心(規格B);(3)滿載&中等偏心(規格B)。
圖5、6、7分別為上述三種工況下內筒前端某節點的位移-轉速曲線,從曲線中可以明顯觀察到共振現象,并讀出考慮陀螺效應影響的驅動部臨界轉速。
將三種工況下的模態分析固有頻率和柔體動力學分析臨界轉速匯列于表1。

圖1 驅動部(內外筒總成)有限元模型

圖2 驅動部(空筒)模態分析(Normal Mode)第一階模態振型圖(57.4Hz)

圖3 驅動部(滿載)模態分析(Normal Mode)第一階模態振型圖(40.3Hz)

圖4 Virtual.lab Motion剛柔混合模型

圖5 內筒前端某節點位移—轉速曲線(空筒&偏心A)

圖6 內筒前端某節點位移—轉速曲線(空筒&偏心B)

圖7 內筒前端某節點位移—轉速曲線(滿載&偏心B)

表1 模態頻率與臨界轉速對比
可以看出,對滾筒洗衣機驅動部而言,陀螺效應使得驅動部臨界轉速降低,因空筒工況下模態頻率所對應的轉速較高,陀螺效應的影響更顯著,所以與模態分析固有頻率相比,臨界轉速降低幅度也更為明顯。
經過以上分析和研究,可以得出以下結論:
(1)由于質量對共振頻率影響明顯,滿載工況下驅動部臨界轉速明顯低于空筒工況,基于安全設計,洗衣機驅動部(內外筒系統)臨界轉速計算應考慮滿載工況;
(2)滾筒洗衣機考慮陀螺效應影響時,與模態分析固有頻率相比,驅動部臨界轉速降低,基于安全設計,洗衣機驅動部臨界轉速計算需考慮陀螺效應影響;
(3)對于大容量、高轉速及復合式分區洗筒結構洗衣機設計,由于質量的增加及陀螺效應影響更需要對其進行柔體動力學分析,以保證結構剛度要求,獲得低振動的最佳用戶體驗。
[1] 師名林,王德忠,張繼革. 基于ANSYS的大型屏蔽電機泵轉子系統建模及動力學分析[J]. 機械科學與技術,2013(32), 3: 358-361.
[2] 隋永楓,呂和詳. 陀螺效應對轉子橫向振動的影響分析[J]. 計算力學學報,2003(20), 6: 711-714.
[3] 劉偉佳. 陀螺效應對轉子臨界轉速的影響[J].吉林師范大學學報,2012,3: 47-49.
[4] 安曉衛,王學永. 渦輪增壓器轉子系統的臨界轉速研究[J]. 機械設計,2015(32), 2: 75-78.
[5] 蔣云帆,廖明夫,劉永泉,等. 同轉/對轉雙轉子系統的動力學特性[J]. 航空動力學報,2013(28), 12: 2771-2780.
[6] 蔡大文. 考慮陀螺力矩時軸的臨界轉速計算[J]. 西安工業學院學報,1995(15), 2: 120-125.
[7] 侯曦,張凱,楊崇倡,等. 高速卷繞機筒管夾頭系統動力學建模與臨界轉速分析[J]. 機械設計與制造,2012, 7: 239-241.
[8] LMS Virtual.lab Online Help.