Tổng quan nghiên cứu

Trong các hệ thống máy công nghiệp hiện đại như tua bin khí, máy nén khí ly tâm, quạt công nghiệp và động cơ điện công suất lớn, tốc độ quay thường đạt từ 3.000 đến trên 30.000 vòng/phút. Hiện tượng mất cân bằng khối lượng của chi tiết quay (rotor) là nguyên nhân hàng đầu gây ra hơn 75% các sự cố rung động và tiếng ồn trong công nghiệp cơ khí. Khi trục quay với vận tốc góc cao, lực ly tâm quán tính tỷ lệ thuận với bình phương tốc độ quay, khiến một độ lệch tâm chỉ vài micromet cũng có thể sinh ra lực động học lên tới hàng chục kilonewton truyền xuống gối đỡ và móng máy. Tình trạng này đẩy nhanh quá trình mòn rão cổ trục, phá hủy màng dầu bôi trơn trong ổ trượt và gây mất ổn định động học toàn hệ thống.

Mục tiêu cốt lõi của nghiên cứu là xây dựng cơ sở khoa học, hệ thống hóa các giải pháp giải tích và thực nghiệm nhằm tối ưu hóa quy trình cân bằng động chi tiết quay cứng (dạng phẳng và dạng dài) trên cả hệ tuyến tính và hệ phi tuyến. Nghiên cứu được hoàn thiện tại Viện Cơ học thuộc Viện Khoa học và Công nghệ Việt Nam kết hợp với Trường Đại học Công nghệ, Đại học Quốc gia Hà Nội. Ý nghĩa thực tiễn của công trình thể hiện qua việc giảm biên độ dao động dư của hệ thống máy từ mức 0,87992 mm xuống còn 0,0513 mm (giảm tới 94,17%), giúp kéo dài tuổi thọ gối đỡ từ 30% đến 50%, đồng thời cắt giảm hơn 40% chi phí và thời gian dừng máy thông qua các kỹ thuật cân bằng động tại hiện trường.

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

Khung lý thuyết áp dụng

Nghiên cứu được thiết lập dựa trên hệ thống lý thuyết dao động cơ học cổ điển và hiện đại, bao gồm:

  • Phương trình vi phân chuyển động Lagrange loại II: Sử dụng để mô hình hóa động học rotor dài nhiều bậc tự do trên gối đàn hồi có cản, với việc xác định động năng, thế năng và các hàm tiêu tán năng lượng suy rộng.
  • Nguyên lý chồng chất nghiệm trong dao động tuyến tính: Áp dụng cho các hệ rotor mà đáp ứng dao động tỷ lệ thuận tuyệt đối với biên độ lực kích động ly tâm.
  • Phương pháp giải tích phi tuyến và phương pháp tham số bé: Dùng để xấp xỉ khai triển Taylor hàm phi tuyến của lực đàn hồi gối đỡ và xác định các biên tần số cộng hưởng cận kề.

Các khái niệm cơ bản được chuẩn hóa theo tiêu chuẩn quốc tế ISO 1940-1 và tiêu chuẩn Việt Nam TCVN 6373:1998:

  • Mất cân bằng rotor: Biểu diễn qua tích số giữa khối lượng lệch tâm và bán kính quay với công thức $U = ur$.
  • Khối lượng cân bằng triệt tiêu: Giá trị khối lượng đối trọng cần bổ sung hoặc khấu trừ sao cho tổng mất cân bằng đưa về không.
  • Phân loại rotor theo tỷ số hình học: Tiêu chuẩn phân định giữa rotor ngắn (dạng đĩa mỏng, tỷ số chiều dài trên đường kính nhỏ hơn 0,5 hoạt động dưới 900 vòng/phút chỉ cần 1 mặt phẳng cân bằng) và rotor dài (tỷ số chiều dài trên đường kính lớn hơn 1,0 đòi hỏi cân bằng trên 2 mặt phẳng độc lập).
  • Cấp độ cân bằng cho phép: Cấp G2.5 cho tua bin hơi nước hoạt động ở tốc độ 5.100 vòng/phút (độ lệch tâm riêng cho phép 4 micromet đối với rotor 11.000 kg, tương ứng mất cân bằng dư cho phép 44 mm.kg) và cấp G6.3 cho máy ly tâm công nghiệp 300 vòng/phút (độ lệch tâm riêng cho phép 200 micromet đối với rotor 2.000 kg, tương ứng lượng mất cân bằng dư 400 mm.kg).

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

  • Nguồn dữ liệu và đối tượng mô hình: Nghiên cứu sử dụng dữ liệu thực nghiệm dao động từ hệ rotor cứng tiêu chuẩn có khối lượng 2.000 kg, tần số riêng 5 Hz, kết hợp phân tích kiểm chứng trên các rotor công nghiệp tải trọng lớn 11.000 kg.
  • Phương pháp chọn mẫu: Lựa chọn mẫu nghiên cứu có chủ đích theo cấu trúc hình học điển hình của máy quay công nghiệp, chia thành 2 nhóm cấu trúc rotor phẳng (bánh đà, đĩa mài, cánh quạt gió) và rotor dài (trục cán, ru-lô con lăn, rotor máy phát điện).
  • Phương pháp phân tích: Kết hợp mô phỏng số Runge-Kutta, phân tích biến đổi Fourier nhanh (FFT) trên miền tần số dao động từ 3 Hz đến 34 Hz, và phương pháp ma trận số phức để giải hệ phương trình hệ số ảnh hưởng 2 mặt phẳng.
  • Lý do lựa chọn: Phương pháp ma trận và phân tích phổ FFT cho phép bóc tách chính xác các hài bậc cao và nhận diện rõ biên độ dao động tại đúng tần số quay, từ đó xử lý triệt để hiện tượng lệch pha mà không bị nhiễu bởi các dao động nền của xưởng sản xuất.
  • Kế hoạch nghiên cứu: Toàn bộ quá trình khảo sát lý thuyết, lập trình mô phỏng và đo đạc thực nghiệm được thực hiện xuyên suốt trong khung thời gian 12 tháng.

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

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

Thứ nhất, đối với rotor phẳng tuyến tính, các phương pháp cân bằng không cần đo pha (như phương pháp 3 lần thử, phương pháp 2 lần thử tạo góc 90 độ, và phương pháp 4 lần chạy máy gắn tải thử cách đều 120 độ) đạt độ chính xác tương đương phương pháp đo pha truyền thống. Khi kiểm chứng tại tần số làm việc 3 Hz, biên độ rung ban đầu $A_0 = 0,87992\text{ mm}$ đã giảm xuống còn $A_{cb1} = 0,12242\text{ mm}$ ngay sau lần hiệu chỉnh đối trọng đầu tiên (khối lượng thử 0,7 kg, đối trọng cân bằng 1,17 kg tại góc 1,248 radian). Sau lần tinh chỉnh thứ hai với tải thử 0,5 kg, biên độ dao động chỉ còn $A_{cb2} = 0,0513\text{ mm}$, đạt tỷ lệ triệt tiêu rung động lên đến 94,17%.

Thứ hai, nghiên cứu đã chỉ ra sự biến dạng phổ phức tạp của rotor trong môi trường phi tuyến (hệ số phi tuyến bậc ba $b_3 = 30$). Trong dải tần số cộng hưởng nguy hiểm từ 3,75 Hz đến 25 Hz (tại tần số kích thích 24 Hz), đồ thị phổ FFT xuất hiện đồng thời hai đỉnh dao động lớn tại 8 Hz và 24 Hz, khiến các công thức cân bằng tuyến tính hoàn toàn mất tác dụng. Tuy nhiên, khi vận hành ngoài miền cộng hưởng tại tần số cao 34 Hz, biên độ dao động giảm mạnh từ 1,21738 mm xuống 0,09898 mm sau cân bằng (đạt hiệu suất giảm rung 91,87%).

Thứ ba, đối với rotor dài trên hai gối đỡ, phương pháp ma trận hệ số ảnh hưởng số phức đã giải quyết dứt điểm sự tương tác dao động chéo giữa hai mặt phẳng. Trong bài toán rotor đối xứng chịu dao động ban đầu 25 micromet lệch pha 210 độ tại gối I và 60 micromet lệch pha 310 độ tại gối II, giải ma trận nghịch đảo đã xác định chính xác đối trọng cần gắn là 1,85 kg tại góc âm 8,86 độ trên mặt phẳng I và 2,03 kg tại góc âm 25,36 độ trên mặt phẳng II.

Thảo luận kết quả

Các kết quả thực nghiệm và mô phỏng khẳng định tính ưu việt của việc kết hợp đo phổ biên độ và góc pha. Nguyên nhân chính khiến phương pháp cân bằng truyền thống thất bại trên hệ phi tuyến là do hiện tượng nhảy biên độ và sự sinh tần số phụ, khiến lực ly tâm không còn tỷ lệ bậc nhất với độ dời.

Dữ liệu nghiên cứu được biểu diễn trực quan thông qua đồ thị véc-tơ cực (đồ thị Nyquist) thể hiện quỹ đạo dịch chuyển pha của các thành phần mất cân bằng tĩnh và mất cân bằng động ngẫu lực. Các bảng đối sánh phổ dao động trước và sau cân bằng minh chứng rằng việc lựa chọn vùng tần số cân bằng thấp hơn tần số riêng (như mức 3 Hz) mang lại độ an toàn kỹ thuật cao hơn nhiều so với việc cố gắng cân bằng ở tốc độ cận kề vùng cộng hưởng. So với các công bố của những chuyên gia quốc tế trong ngành động lực học máy, mô hình tính toán trong luận văn này đã đơn giản hóa đáng kể quy trình thao tác ngoài hiện trường nhờ tối ưu hóa thuật toán tính hệ số ảnh hưởng $a_{ij}$.

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

  1. Chuẩn hóa quy trình cân bằng động tại hiện trường: Các doanh nghiệp sản xuất và bảo dưỡng cơ khí cần áp dụng đồng bộ quy trình cân bằng 3 lần thử hoặc 4 lần chạy máy cách đều 120 độ cho các quạt công nghiệp và máy bơm ly tâm, nhằm đưa mức rung động vận hành về dưới ngưỡng tiêu chuẩn 0,05 mm/s; thời gian triển khai áp dụng thí điểm trong vòng 3 tháng đầu năm.
  2. Ứng dụng phần mềm tự động hóa tính toán ma trận hệ số ảnh hưởng: Phòng kỹ thuật cơ điện tại các nhà máy cần tích hợp thuật toán xử lý số phức cho rotor dài 2 mặt phẳng lên máy tính công nghiệp, giúp rút ngắn thời gian tính toán đối trọng xuống dưới 30 phút cho mỗi ca hiệu chỉnh; hoàn thành cài đặt trong vòng 6 tháng.
  3. Thiết lập dải tần số an toàn khi cân bằng hệ phi tuyến: Đội ngũ kỹ sư vận hành phải thực hiện đo quét tần số riêng trước khi cân bằng, chủ động tránh xa vùng tần số cộng hưởng phi tuyến từ 3,75 Hz đến 25 Hz và ưu tiên tiến hành cân bằng ở tần số thấp dưới 3 Hz để đảm bảo độ chính xác trên 90%; áp dụng định kỳ trong quy trình bảo dưỡng hàng quý.
  4. Đầu tư nâng cấp thiết bị đo phổ FFT và kiểm định theo tiêu chuẩn quốc gia: Ban quản lý nhà máy cần trang bị cảm biến gia tốc đa kênh kết hợp máy phân tích phổ động năng nhằm kiểm soát chất lượng cân bằng đạt chuẩn G2.5 cho tua bin và G6.3 cho các máy quay thông thường theo TCVN 6373:1998; lộ trình đầu tư thực hiện trong 12 tháng.

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

  • Kỹ sư cơ điện và chuyên viên bảo trì nhà máy: Tiếp cận quy trình cân bằng động tại chỗ chuẩn xác, giúp xử lý nhanh các sự cố rung lắc của quạt lò nung, máy nghiền, bơm cao áp mà không phải tháo dỡ thiết bị vận chuyển về xưởng.
  • Học viên cao học và nghiên cứu sinh ngành Cơ học kỹ thuật: Khai thác mô hình toán học giải tích Lagrange loại II, các phương pháp xấp xỉ dao động phi tuyến và thuật toán ma trận phục vụ phát triển đề tài nghiên cứu chuyên sâu.
  • Các đơn vị chế tạo máy và thiết bị quay công nghiệp: Ứng dụng bộ công thức tính toán hệ số ảnh hưởng $a_{ij}$ để lập trình phần mềm điều khiển cho các dòng máy cân bằng động chuyên dụng sản xuất hàng loạt.
  • Cán bộ quản lý an toàn và kiểm định thiết bị công nghiệp: Sử dụng hệ thống chỉ tiêu rung động và bảng cấp độ mất cân bằng theo tiêu chuẩn ISO 1940-1 làm căn cứ nghiệm thu kỹ thuật và đánh giá an toàn máy móc định kỳ.

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

  1. Khi nào nên cân bằng rotor trên 1 mặt phẳng và khi nào bắt buộc phải cân bằng trên 2 mặt phẳng? Rotor có tỷ số chiều dài trên đường kính nhỏ hơn 0,5 và tốc độ dưới 900 đến 1.800 vòng/phút (như cánh quạt mỏng, đĩa mài) chỉ cần cân bằng trên 1 mặt phẳng. Rotor có tỷ số chiều dài trên đường kính lớn hơn 1,0 (như trục dài, ru-lô) bắt buộc phải cân bằng trên 2 mặt phẳng để triệt tiêu đồng thời cả mất cân bằng tĩnh lẫn mất cân bằng ngẫu lực.

  2. Cân bằng động tại hiện trường có ưu điểm gì vượt trội so với cân bằng trên máy chuyên dùng tại xưởng? Cân bằng tại hiện trường giúp tiết kiệm hơn 50% thời gian dừng máy do không cần tháo rời rotor, giảm chi phí vận chuyển cơ giới, đồng thời đảm bảo rotor được cân chỉnh chính xác ngay trong điều kiện lắp ráp thực tế cùng hệ thống gối đỡ, móng máy và tải trọng vận hành.

  3. Công thức nào giúp lựa chọn khối lượng thử an toàn khi chưa biết trước lượng mất cân bằng? Khối lượng thử $m_t$ được chọn qua công thức kinh nghiệm $m_t = (Mg / n^2 r) \cdot 100 \cdot k$ (với hệ số $k$ từ 0,2 đến 0,5 tùy khối lượng máy), hoặc chọn tải thử bằng khoảng 5 lần lượng mất cân bằng dư cho phép theo cấp chất lượng ISO 1940-1 nhằm đảm bảo đủ độ nhạy tín hiệu mà không gây quá tải cho ổ đỡ.

  4. Tại sao không được áp dụng các phương pháp cân bằng tuyến tính trong vùng cộng hưởng của hệ phi tuyến? Trong vùng tần số cộng hưởng từ 3,75 Hz đến 25 Hz của hệ phi tuyến, dao động xuất hiện nhiều đỉnh phổ phụ (như 8 Hz kết hợp với 24 Hz) và biên độ dao động không còn tỷ lệ bậc nhất với lực ly tâm, khiến các phép tính véc-tơ tuyến tính bị sai lệch hoàn toàn.

  5. Phương pháp ma trận hệ số ảnh hưởng đem lại lợi ích gì cho các cơ sở sản xuất hàng loạt? Đối với các rotor cùng chủng loại lắp trên cùng hệ gối đỡ, ma trận hệ số ảnh hưởng $a_{ij}$ giữ giá trị không đổi. Do đó, cơ sở sản xuất chỉ cần đo dao động ban đầu một lần duy nhất là tính toán được ngay vị trí và khối lượng đối trọng mà không cần chạy thử gắn tải lặp lại.

Kết luận

  • Luận văn đã giải quyết toàn diện bài toán cân bằng động cho cả chi tiết quay phẳng và chi tiết quay dài, thiết lập hệ thống phương pháp thực nghiệm không cần đo pha và có đo pha đạt độ chuẩn xác cao.
  • Đã xây dựng thành công mô hình toán học giải tích dao động rotor dài dựa trên phương trình Lagrange loại II và làm chủ thuật toán ma trận hệ số ảnh hưởng số phức cho bài toán cân bằng hai mặt phẳng.
  • Đã khám phá bản chất động học của hệ rotor phi tuyến, chỉ ra dải tần số cộng hưởng nguy hiểm và chứng minh hiệu quả của phương pháp cân bằng đa bước ở tần số thấp với độ giảm rung đạt 94,17%.
  • Đóng góp giải pháp công nghệ cân bằng tại chỗ có tính ứng dụng cao, hỗ trợ đắc lực cho ngành cơ khí chế tạo và bảo trì công nghiệp tại Việt Nam.
  • Các đơn vị sản xuất và viện nghiên cứu nên đẩy mạnh chuyển giao thuật toán này vào các thiết bị đo kiểm hiện đại trong vòng 6 đến 12 tháng tới để tối ưu hóa năng suất và an toàn vận hành.