Tổng quan nghiên cứu

Trong ngành công nghiệp sản xuất thiết bị gia dụng hiện đại, máy giặt lồng ngang chiếm hơn 65% thị phần tại các quốc gia phát triển nhờ hiệu suất giặt vượt trội, khả năng tiết kiệm nước gấp 3 lần so với máy giặt lồng đứng và tốc độ vắt ly tâm đạt từ 1.000 đến 1.400 vòng/phút. Mặc dù mang lại hiệu quả giặt vắt cao, máy giặt lồng ngang phải đối mặt với thách thức kỹ thuật nghiêm trọng về rung động và tiếng ồn sinh ra do sự phân bố không đồng đều của khối lượng vật phẩm giặt trong lồng quay. Lực quán tính ly tâm từ khối lượng mất cân bằng khi quay ở vận tốc cao có thể khiến biên độ dao động của lồng chứa vượt ngưỡng an toàn 5 mm, dẫn đến hiện tượng va chạm mạnh vào thành vỏ máy, gây hư hỏng cơ khí và làm máy giặt tự dịch chuyển trên sàn nhà.

Đề tài luận văn thạc sĩ "Mô hình hóa rung động lồng giặt của máy giặt lồng ngang" do học viên Nguyễn Thị Bích Ngọc thực hiện dưới sự hướng dẫn khoa học của PGS. Ngô Như Khoa tại Trường Đại học Kỹ thuật Công nghiệp – Đại học Thái Nguyên (năm 2019) tập trung giải quyết bài toán động lực học then chốt này. Mục tiêu cốt lõi của nghiên cứu là xây dựng mô hình toán học 2D mô tả chính xác chuyển động rung động của cụm lồng giặt; kiểm chứng và đánh giá độ tin cậy của mô hình thông qua thực nghiệm trên máy giặt thương mại LG WD-8990TDS; đồng thời khảo sát định lượng mức độ ảnh hưởng của độ cứng lò xo, góc nghiêng và đặc tính giảm chấn đến đáp ứng dao động của toàn hệ thống.

Nghiên cứu được hoàn thành trong giai đoạn 2018 - 2019, đóng góp trực tiếp vào đề tài nghiên cứu khoa học cấp Bộ và công bố 02 báo cáo khoa học tại Hội nghị Quốc tế ICERA 2018. Kết quả nghiên cứu cung cấp cơ sở dữ liệu phân tích động lực học chuẩn xác với khối lượng lồng chứa 32,5 kg và tải lệch tâm 0,62 kg, mở ra phương án tối ưu hóa hệ thống treo giúp giảm thiểu biên độ rung động từ 15% đến 25% mà không làm tăng chi phí sản xuất thiết bị.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Luận văn vận dụng lý thuyết cơ học chuyển động nhiều vật (multibody dynamics), nguyên lý D'Alembert cho cơ hệ phẳng và lý thuyết tiêu tán năng lượng trong hệ dao động chịu cản ma sát. Để mô tả chuyển động của lồng giặt gắn trên hệ thống treo gồm 2 lò xo và 3 giảm chấn, tác giả xây dựng 2 mô hình toán học 2 bậc tự do (chuyển vị ngang $x$ và chuyển vị đứng $y$):

  1. Mô hình 2D tuyến tính (Mô hình 1): Giả thiết lực cản ma sát của bộ giảm chấn tỷ lệ thuận bậc nhất với vận tốc chuyển dịch ($F = C \cdot v$), sử dụng hệ số cản nhớt không đổi $C = 230\text{ N}\cdot\text{s/m}$ và độ cứng lò xo giãn dài $K = 5880\text{ N/m}$.
  2. Mô hình 2D phi tuyến (Mô hình 2): Kể đến bản chất làm việc thực tế của giảm chấn thủy lực ma sát nửa ướt, trong đó lực cản ma sát là hàm phi tuyến bậc cao của vận tốc dạng đa thức $F = C_0 v^2 + C_1 v + C_2$, với các hệ số thực nghiệm tiêu biểu như $C_0 = 47,9$.

Các khái niệm chuyên ngành nòng cốt được thiết lập bao gồm: lực kích thích quán tính ly tâm $F(t) = m \omega^2 r$, góc lệch lò xo treo $\theta_S = 10^\circ$, góc lệch giảm chấn $\theta_D = 33^\circ$, khối lượng quy đổi toàn phần $M = 32,5\text{ kg}$, và hàm năng lượng hao tán dập tắt rung động của hệ thống treo.

Phương pháp nghiên cứu

Nghiên cứu áp dụng phương pháp tích hợp giữa mô hình hóa giải tích, giải số trên phần mềm mô phỏng và kiểm chứng thực nghiệm trực tiếp:

  • Thiết lập mẫu và thông số thí nghiệm: Đối tượng khảo sát là máy giặt LG WD-8990TDS với tổng khối lượng lồng chứa và đối trọng $M = 32,5\text{ kg}$, khối lượng mất cân bằng giả lập $m = 0,62\text{ kg}$ đặt tại bán kính quay $r = 0,2\text{ m}$, và dải tốc độ khảo sát thực nghiệm tại các mốc 400 vòng/phút, 611 vòng/phút và 764,77 vòng/phút.
  • Phương pháp chọn mẫu và bố trí đo kiểm: Tác giả xây dựng 2 hệ thống thử nghiệm độc lập: (1) Băng thử chuyên dụng xác định đường đặc tính lực - vận tốc (F-V) của ống giảm chấn qua cơ cấu tay quay con trượt và động cơ điều tốc; (2) Hệ khung đo động lực học trang bị 5 cảm biến tải trọng Loadcell (đo phản lực tại 2 điểm treo lò xo và 3 gối đỡ giảm chấn) cùng 2 cảm biến dịch chuyển LVDT (đo chuyển vị hai chiều $x, y$) kết nối máy tính qua card thu thập dữ liệu và phần mềm điều khiển LabVIEW.
  • Lý do lựa chọn phương pháp phân tích: MATLAB/Simulink được lựa chọn để giải hệ phương trình vi phân phi tuyến bậc hai bằng phương pháp tích phân số Runge-Kutta vì khả năng xử lý nhanh các hàm điều kiện trễ ma sát signum, kết hợp với phương pháp hồi quy dữ liệu thực nghiệm F-V để trích xuất chính xác bộ thông số cản $C_0, C_1, C_2$ trong chu trình đo từ 0 đến 60 giây.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng số học và đối chứng thực nghiệm đã làm sáng tỏ 4 phát hiện quan trọng:

  1. Sự vượt trội của mô hình giảm chấn phi tuyến: Tại tốc độ quay vắt $N = 609 - 611\text{ vòng/phút}$, Mô hình 2 (phi tuyến) cho kết quả biên độ dịch chuyển phương ngang $x$ trùng khớp với số liệu thực nghiệm LVDT với sai số dưới 6,5%, trong khi Mô hình 1 (tuyến tính) tạo ra sai số lớn lên tới 18,2% do tuyến tính hóa lực cản.
  2. Đặc tính dao động tại vùng cộng hưởng: Tại tốc độ chuyển tiếp $N = 611\text{ vòng/phút}$, biên độ dao động của lồng giặt đạt giá trị cực đại với chuyển vị đứng $y$ xấp xỉ 4,8 mm và chuyển vị ngang $x$ đạt 3,2 mm. Khi máy giặt tăng tốc vượt qua vùng cộng hưởng lên trạng thái quay ổn định $N = 764,77\text{ vòng/phút}$, biên độ dao động giảm mạnh hơn 42%, giữ ổn định quanh mức 1,8 mm.
  3. Vùng góc nghiêng dập tắt dao động tối ưu: Năng lượng tiêu tán của hệ thống treo biến thiên mạnh theo góc nghiêng của ống giảm chấn; tại dải góc $\theta_D = 33^\circ \pm 3^\circ$, bộ 3 giảm chấn hấp thụ hơn 78% tổng động năng rung động phát sinh, giúp hạn chế hiện tượng xô lệch phương ngang.
  4. Hiệu quả phân bố lực của cấu hình 3 giảm chấn: Bố trí 2 giảm chấn bên trái và 1 giảm chấn bên phải giúp cân bằng triệt để mô-men xoắn phụ do động cơ truyền động gây ra, giảm biên độ quỹ đạo chuyển vị tâm lồng hơn 26% so với kết cấu 2 giảm chấn truyền thống ở tốc độ vắt 800 vòng/phút.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự sai khác giữa tính toán lý thuyết cổ điển và thực tế vận hành là tính chất trễ phi tuyến (hysteresis loop) của dầu giảm chấn thủy lực khi đổi chiều chuyển dịch với gia tốc cao. Dưới tác động của lực kích thích biến thiên tuần hoàn $F_x(t) = m \omega^2 r \cos(\omega t)$ và $F_y(t) = m \omega^2 r \sin(\omega t)$, lực cản xuất hiện thành phần ma sát phi tuyến $C_0 v^2$ đóng vai trò ghìm giữ chuyển vị tại thời điểm đảo chiều vận tốc.

So sánh với kết quả nghiên cứu của Pınar Boyraz (2013) ở tốc độ 900 vòng/phút và mô hình nhiều vật của T. Argentini (2018), mô hình 2D phi tuyến trong luận văn thể hiện độ tin cậy cao với sai số biên độ tổng hợp chỉ xấp xỉ 5,3%. Dữ liệu nghiên cứu được biểu diễn trực quan qua các đồ thị hồi quy đặc tính lực - vận tốc F-V, đồ thị miền thời gian trích xuất từ 5 kênh Loadcell, và biểu đồ pha mô tả quỹ đạo chuyển vị hình elip khép kín của tâm trục lồng giặt ở các mốc 600 vòng/phút và 800 vòng/phút. Phân tích bảng dữ liệu phản lực liên kết khẳng định lực động truyền xuống chân đế giảm rõ rệt khi các thông số góc nghiêng được hiệu chỉnh tối ưu.

Đề xuất và khuyến nghị

Dựa trên các kết luận định lượng từ mô hình toán và thực nghiệm, luận văn đưa ra 4 khuyến nghị kỹ thuật có tính ứng dụng cao:

  1. Hiệu chỉnh góc lắp đặt cụm hệ thống treo: Doanh nghiệp sản xuất cần cố định góc nghiêng làm việc của bộ giảm chấn $\theta_D = 33^\circ$ và góc nghiêng lò xo treo $\theta_S = 10^\circ \pm 2^\circ$. Target metric: Giảm từ 18% đến 25% biên độ chuyển vị ngang của lồng giặt trong suốt chu trình vắt. Timeline: 3 đến 6 tháng trong khâu thiết kế kết cấu khung vỏ mới. Chủ thể thực hiện: Kỹ sư phòng R&D cơ khí tại các nhà máy sản xuất thiết bị gia dụng.
  2. Ứng dụng thuật toán cản phi tuyến vào phần mềm thiết kế: Thay thế hoàn toàn mô hình cản tuyến tính bằng hàm cản phi tuyến bậc hai ($C_0 = 45 - 50$, $C_1 = 180 - 220$) khi tính toán mô phỏng động lực học ảo. Target metric: Nâng cao độ chính xác dự báo phản lực tác dụng lên vỏ máy đạt trên 92%. Timeline: 6 tháng cập nhật quy trình tính toán CAE/NVH. Chủ thể: Kỹ sư mô phỏng động lực học công nghiệp.
  3. Chuẩn hóa cấu hình 3 giảm chấn bất đối xứng: Duy trì kết cấu bố trí 2 ống giảm chấn tại phía chịu tải trọng động cơ và 1 ống giảm chấn tại phía đối diện để triệt tiêu mô-men xoắn cục bộ. Target metric: Giảm lực truyền động xuống nền sàn dưới 120 N ở chế độ vắt cao tốc 1.200 vòng/phút. Timeline: 1 năm áp dụng trên dây chuyền sản xuất hàng loạt. Chủ thể: Bộ phận quản lý chất lượng và thử nghiệm độ bền sản phẩm.
  4. Tích hợp chiến lược điều khiển tăng tốc thông minh: Lập trình thuật toán điều khiển biến tần động cơ vượt nhanh qua dải vận tốc cộng hưởng nguy hiểm từ 580 đến 640 vòng/phút trong khoảng thời gian dưới 3,5 giây. Target metric: Cắt giảm hơn 60% thời gian lưu trú tại đỉnh cộng hưởng, ngăn chặn triệt để hiện tượng va đập thùng giặt. Timeline: 9 tháng phát triển và nạp vi mã điều khiển động cơ. Chủ thể: Kỹ sư lập trình hệ thống nhúng và điều khiển tự động.

Đối tượng nên tham khảo luận văn

  1. Học viên cao học và sinh viên ngành Kỹ thuật Cơ khí, Cơ điện tử: Tiếp cận trọn vẹn phương pháp luận thiết lập phương trình vi phân chuyển động bằng nguyên lý D'Alembert, kỹ năng xây dựng sơ đồ khối mô phỏng trên MATLAB/Simulink và kỹ thuật thu thập dữ liệu cảm biến đa kênh thông qua phần mềm LabVIEW.
  2. Kỹ sư R&D thiết kế máy giặt và thiết bị quay công nghiệp: Ứng dụng trực tiếp bộ thông số độ cứng lò xo $K = 5880\text{ N/m}$ và hàm cản phi tuyến để tối ưu hóa kết cấu máy thực tế, rút ngắn 30% thời gian chế tạo và thử nghiệm các cụm nguyên mẫu cơ khí.
  3. Giảng viên và nhà nghiên cứu trong lĩnh vực Động lực học & Dao động máy: Sử dụng các số liệu thực nghiệm đối chứng làm tài liệu giảng dạy chuyên sâu, đồng thời phát triển các hướng nghiên cứu nâng cao về hệ thống treo bán chủ động sử dụng giảm chấn từ biến (MR damper).
  4. Doanh nghiệp sản xuất linh kiện giảm xóc và phụ tùng gia dụng: Tham khảo quy trình và sơ đồ đồ gá băng thử nghiệm lực - vận tốc F-V để kiểm chuẩn chất lượng xuất xưởng cho các dòng giảm chấn thủy lực với độ tin cậy đạt trên 95%.

Câu hỏi thường gặp

1. Tại sao máy giặt lồng ngang lại xuất hiện rung lắc mạnh nhất ở dải tốc độ từ 600 đến 650 vòng/phút? Dải tốc độ từ 600 đến 650 vòng/phút tương ứng với tần số kích thích dao động từ 10 Hz đến 10,8 Hz, trùng khớp với tần số dao động riêng của hệ thống treo (gồm khối lượng lồng $M = 32,5\text{ kg}$ và 2 lò xo $K = 5880\text{ N/m}$). Hiện tượng cộng hưởng cơ học xuất hiện làm biên độ dao động tăng vọt lên mức 4,8 mm trước khi giảm dần khi vòng quay đạt mức ổn định cao hơn.

2. Sự khác biệt cốt lõi giữa mô hình cản tuyến tính và mô hình cản phi tuyến trong luận văn là gì? Mô hình 1 coi lực cản là hàm bậc nhất $F = C \cdot v$ với $C = 230\text{ N}\cdot\text{s/m}$, trong khi Mô hình 2 mô tả chính xác lực cản thực tế qua hàm phi tuyến bậc hai $F = C_0 v^2 + C_1 v + C_2$ với $C_0 = 47,9$. Mô hình phi tuyến giúp hạ thấp sai số dự báo chuyển vị từ 18,2% xuống dưới 6,5% khi so sánh với dữ liệu đo thực tế.

3. Việc bố trí 3 bộ giảm chấn bất đối xứng mang lại ưu thế kỹ thuật gì? Cấu hình 3 giảm chấn (2 chiếc bên trái, 1 chiếc bên phải) giúp triệt tiêu mô-men quán tính xoắn sinh ra do cụm động cơ điện được lắp lệch về một bên thân máy giặt. Cách bố trí này giúp tăng khả năng hấp thụ năng lượng rung động lên hơn 26% so với kết cấu 2 giảm chấn đặt đối xứng.

4. Khối lượng mất cân bằng 0,62 kg ảnh hưởng như thế nào đến lực quán tính ly tâm? Khối lượng mất cân bằng $m = 0,62\text{ kg}$ quay ở bán kính $r = 0,2\text{ m}$ tại vận tốc 764,77 vòng/phút tạo ra lực kích thích ly tâm $F(t) = m \omega^2 r$ lên tới 79,5 N. Lực này đảo chiều liên tục theo chu kỳ quay, là nguyên nhân trực tiếp gây ra rung động cưỡng bức và truyền rung chấn nguy hiểm xuống sàn nhà.

5. Dữ liệu thực nghiệm của luận văn được kiểm chứng qua các thiết bị đo lường chuẩn xác nào? Hệ thống thử nghiệm sử dụng 5 cảm biến lực Loadcell đo phản lực tại các điểm liên kết lò xo và chân giảm chấn, kết hợp 2 cảm biến LVDT đo chuyển vị hai chiều $x-y$, toàn bộ tín hiệu được số hóa qua card thu nhận dữ liệu và xử lý trên nền tảng LabVIEW với độ chính xác cao trên 95%.

Kết luận

  • Xây dựng thành công hệ phương trình vi phân mô tả chuyển động rung động 2D của lồng giặt máy giặt lồng ngang có xét đến tính chất phi tuyến của giảm chấn.
  • Chứng minh mô hình giảm chấn phi tuyến với hệ số $C_0 = 47,9$ đạt độ chính xác cao vượt trội với sai số thực nghiệm dưới 6,5% tại tốc độ quay cộng hưởng 611 vòng/phút.
  • Xác định góc nghiêng làm việc tối ưu của giảm chấn $\theta_D = 33^\circ$ cùng kết cấu 3 ống giảm xóc giúp tiêu tán hơn 78% động năng rung lắc của lồng chứa.
  • Hoàn thành bộ cơ sở dữ liệu thực nghiệm toàn diện trên máy giặt LG WD-8990TDS với khối lượng $32,5\text{ kg}$ và tải lệch tâm $0,62\text{ kg}$.
  • Kế hoạch phát triển trong 6 đến 12 tháng tới: Mở rộng sang mô hình không gian 3D sáu bậc tự do và thử nghiệm giải pháp giảm chấn từ biến bán chủ động MR.

Hãy khai thác và ứng dụng ngay các giải pháp mô hình hóa cùng bộ thông số tối ưu từ luận văn này vào quy trình R&D để nâng cao chất lượng vận hành và triệt tiêu tiếng ồn cho các dòng máy giặt thế hệ mới!