Tổng quan nghiên cứu
Trong xu hướng phát triển của ngành kỹ thuật cơ điện tử hiện đại, nhu cầu tối ưu hóa khối lượng và nâng cao tốc độ vận hành của các cơ cấu truyền động chính xác ngày càng trở nên cấp thiết. Các cấu trúc thanh mềm công xôn (flexible cantilever beam) được ứng dụng rộng rãi từ các cánh tay robot công nghiệp nhẹ, anten vệ tinh không gian cho đến các cảm biến vi cơ điện tử MEMS. Tuy nhiên, việc giảm khối lượng kết cấu thường đi kèm với hiện tượng suy giảm độ cứng vững, khiến hệ thống dễ bị kích động dao động dư nghiêm trọng khi bệ gá di chuyển tịnh tiến ở gia tốc lớn. Dao động dư này làm suy giảm độ chính xác định vị và kéo dài thời gian ổn định của cơ cấu.
Nghiên cứu của tác giả Vũ Nguyễn Trí Giang dưới sự hướng dẫn của Phó Giáo sư Nguyễn Quốc Chí tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh đã tập trung giải quyết bài toán mô hình hóa động lực học và điều khiển triệt tiêu dao động cho thanh dầm mềm công xôn gắn trên bệ trượt tịnh tiến. Mục tiêu trọng tâm là xây dựng mô hình toán học rời rạc có kể đến cơ chế cản dao động, thiết kế bộ điều khiển nạp xung trước (Input Shaping) đa mode để triệt tiêu dao động dư, đồng thời mở rộng mô hình hóa hệ thống trong không gian ba chiều (3D).
Được thực hiện trong giai đoạn từ tháng 08 năm 2019 đến tháng 06 năm 2020, đề tài đã tiến hành kiểm chứng mô hình lý thuyết và thuật toán điều khiển trên hệ thống thực nghiệm bao gồm động cơ tuyến tính Yokogawa, bộ điều khiển UMAC và camera tốc độ cao Keyence. Kết quả nghiên cứu mang ý nghĩa thực tiễn lớn khi giúp giảm biên độ dao động dư của thanh dầm từ mức chưa điều khiển xuống dưới 0,3% với bộ định hình hai mode ZVDD, đồng thời rút ngắn thời gian đáp ứng ổn định của hệ thống xuống còn khoảng 0,9 đến 1,2 giây.
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 ứng dụng lý thuyết dầm Euler-Bernoulli kết hợp với Nguyên lý Hamilton mở rộng để thiết lập phương trình vi phân chuyển động liên tục của thanh dầm mềm chuyển động tịnh tiến. Dầm được giả thiết có tiết diện không biến dạng trong mặt phẳng và vuông góc với trục trung hòa sau biến dạng. Do phương trình đạo hàm riêng (PDE) mô tả hệ thống có bậc vô hạn gây khó khăn cho việc phân tích, phương pháp rời rạc hóa Galerkin được áp dụng để chuyển đổi hệ PDE thành hệ phương trình vi phân thường (ODE) hữu hạn bậc.
Nhằm mô tả chính xác sự suy hao năng lượng trong thực tế, mô hình cản tỷ lệ Rayleigh ($C = a_0 M + a_1 K$) được tích hợp vào hệ phương trình động lực học không cản. Về mặt điều khiển, lý thuyết định hình xung (Input Shaping) của Singer và Smith được sử dụng như một giải pháp điều khiển vòng hở tiến tiến (feedforward). Bằng cách tích chập chuỗi xung lực kích thích với tín hiệu điều khiển ban đầu dựa trên tần số dao động riêng và tỷ số cản của hệ thống, các sóng dao động do các xung tạo ra sẽ tự triệt tiêu lẫn nhau. Các dạng bộ định hình xung được khảo sát bao gồm Zero Vibration (ZV), Zero Vibration Derivative (ZVD), Zero Vibration Derivative-Derivative (ZVDD) cùng các cấu trúc nạp xung đa mode (Multi-mode input shapers). Đối với không gian 3 chiều, giả thiết thanh không dãn (inextensional beam) khi tỷ số chiều dài trên bán kính quán tính đạt trên 100 được áp dụng cùng phép biến đổi tọa độ chính để tách rời 3 phương dao động uốn và dọc trục.
Phương pháp nghiên cứu
Hệ thống thực nghiệm được thiết lập với mẫu dầm bằng thép không gỉ có chiều dài $L = 0,9\text{ m}$, chiều rộng $h = 30\text{ mm}$, chiều dày $b = 1\text{ mm}$, mật độ khối lượng theo chiều dài $\rho = 0,235\text{ kg/m}$ và mô đun đàn hồi $E = 200 \times 10^9\text{ N/m}^2$, gắn trên bệ gá có khối lượng $m = 2,5\text{ kg}$. Việc lựa chọn kích thước mẫu này nhằm đảm bảo tỷ số độ mảnh lớn, phù hợp hoàn toàn với giả thiết dầm Euler-Bernoulli.
Để xác định chính xác các thông số dao động phục vụ thiết kế bộ điều khiển, phương pháp kích thích quét tần số (sine sweep) từ 0 Hz đến 8 Hz được thực hiện thông qua động cơ tuyến tính Yokogawa điều khiển bởi card chuyển động NI PCI và phần mềm LabVIEW. Dữ liệu chuyển vị đầu thanh dầm được thu thập trực tiếp bằng camera quang học tốc độ cao Keyence với tần số lấy mẫu cao. Quá trình xử lý tín hiệu và nhận dạng thông số modal được tiến hành thông qua hai hướng tiếp cận: phương pháp suy giảm logarit (Logarithmic Decrement) kết hợp biến đổi Fourier nhanh (FFT), và phương pháp phân tích modal thực nghiệm (Experimental Modal Analysis - EMA) sử dụng hàm modalfrf và modalfit với thuật toán xấp xỉ phân thức hữu tỷ bình phương bé nhất trong phần mềm MATLAB. Toàn bộ quy trình nghiên cứu thực nghiệm và mô phỏng được thực hiện liên tục trong thời gian 10 tháng.
Kết quả nghiên cứu và thảo luận
Những phát hiện chính
Quá trình phân tích phổ phản hồi tần số (FRF) và thực nghiệm nhận dạng thông số hệ thống đã chỉ ra hai tần số dao động tự nhiên chủ đạo của thanh dầm dẹt lần lượt là $f_1 \approx 1,34\text{ Hz}$ ở mode thứ nhất và $f_2 \approx 8,4\text{ Hz}$ ở mode thứ hai. Tỷ số cản tương ứng của kết cấu được xác định ở mức xấp xỉ $\zeta_1 = 0,0071$ và $\zeta_2 = 0,0074$. Nghiên cứu phát hiện rằng dạng kích thích quét tần số hàm mũ (Exponential Sine Sweep với tốc độ tăng trưởng 10 đến 15) mang lại độ chính xác nhận dạng modal vượt trội so với quét tần số tuyến tính thông thường.
Khi áp dụng phương pháp điều khiển định hình xung đơn mode cho mode dao động thứ nhất, bộ định hình ZV chỉ triệt tiêu dao động dư về mức khoảng 5,0% so với biên độ ban đầu khi chưa điều khiển, trong khi bộ định hình ZVD và ZVDD duy trì mức dao động dư khoảng 4,5%. Nguyên nhân là do năng lượng dao động của mode thứ hai vẫn còn tồn tại trong hệ thống. Tuy nhiên, khi chuyển sang áp dụng bộ định hình xung đa mode (Two-mode shapers) tích hợp đồng thời cả hai tần số $1,34\text{ Hz}$ và $8,4\text{ Hz}$, hiệu quả triệt tiêu dao động tăng vọt rõ rệt. Bộ định hình Two-mode ZV đưa dao động dư về mức 1,7% với thời gian tăng tốc đạt 0,90 giây. Đáng chú ý nhất, bộ định hình Two-mode ZVD và Two-mode ZVDD đã triệt tiêu gần như hoàn toàn rung động, đưa dao động dư về mức chỉ 0,3%, đồng thời hạ thấp độ vọt lố xuống dưới 0,02%.
Thảo luận kết quả
Hiệu quả triệt tiêu dao động khi nâng cấp từ bộ điều khiển đơn mode sang đa mode đạt mức cải thiện hơn 50% về độ chính xác biên độ. Điều này chứng minh rằng đối với các cấu trúc dầm mềm có độ đàn hồi cao, việc bỏ qua ảnh hưởng của các mode dao động bậc cao sẽ khiến hệ thống không thể đạt được độ chính xác định vị dưới ngưỡng sai số 2%.
Trong quá trình đánh giá hiệu năng, một sự đánh đổi rõ ràng giữa tính bền vững (robustness) và thời gian trễ đáp ứng đã được lượng hóa. Bộ định hình Two-mode ZV có thời gian tăng tốc nhanh nhất (0,90 giây) nhưng nhạy cảm hơn với sai lệch thông số mô hình. Ngược lại, bộ định hình Two-mode ZVDD đạt độ ổn định tuyệt đối trước biến thiên tần số và loại bỏ hoàn toàn hiện tượng rung lắc với mức dư 0,3%, nhưng đòi hỏi thời gian xác lập lâu hơn, đạt 1,20 giây (tăng thêm 33,3% thời gian so với ZV). Các đồ thị đáp ứng miền thời gian và bảng so sánh các chỉ số hiệu năng (thời gian tăng, độ vọt lố, phần trăm dao động dư) được xây dựng nhất quán, khẳng định tính chính xác của mô hình lý thuyết Galerkin kết hợp giảm chấn Rayleigh so với đáp ứng đo đạc thực tế từ camera tốc độ cao.
Đề xuất và khuyến nghị
Dựa trên các kết quả mô phỏng và thực nghiệm đạt được, nghiên cứu đưa ra 4 khuyến nghị và giải pháp kỹ thuật cụ thể nhằm ứng dụng và phát triển công nghệ điều khiển cơ cấu thanh mềm:
- Triển khai ứng dụng bộ điều khiển nạp xung đa mode (Two-mode ZV/ZVD) vào hệ thống truyền động của các cánh tay robot công nghiệp nhẹ và máy gắp đặt linh kiện điện tử tốc độ cao. Giải pháp này giúp giảm biên độ rung lắc của khâu chấp hành xuống dưới 1,0% mà không đòi hỏi lắp đặt thêm cảm biến phản hồi đắt tiền, hướng tới áp dụng thử nghiệm trong vòng 6 đến 12 tháng tới bởi các kỹ sư R&D cơ điện tử.
- Tích hợp thuật toán nhận dạng thông số trực tuyến (Online Modal Identification) sử dụng tín hiệu quét tần số hàm mũ kết hợp bộ lọc thích nghi. Mục tiêu là tự động bù sai lệch tần số dao động riêng khi cơ cấu mang tải trọng thay đổi từ 0,5 kg đến 5 kg, hoàn thiện trong giai đoạn 1 năm tiếp theo bởi các nhóm nghiên cứu phần mềm điều khiển nhúng.
- Mở rộng kiểm chứng thực nghiệm mô hình dầm mềm 3D cho các cơ cấu không gian phức tạp. Cần trang bị hệ thống cảm biến quang học ba chiều hoặc cảm biến đo biến dạng (strain gauge) đa điểm để đo đạc đồng thời dao động uốn hai phương và dao động dọc trục, triển khai trong khung thời gian 2024–2026 tại các phòng thí nghiệm trọng điểm.
- Chuẩn hóa quy trình ước lượng hệ số cản Rayleigh trong các phần mềm thiết kế và mô phỏng động lực học cơ cấu máy. Các doanh nghiệp chế tạo máy chính xác cần áp dụng phương pháp phân tích modal thực nghiệm EMA để thu thập dữ liệu cản thực tế, nâng cao độ tin cậy của mô hình số lên trên 95% trước khi chế tạo mẫu thử cơ khí.
Đối tượng nên tham khảo luận văn
Nội dung và phương pháp luận của luận văn mang lại giá trị học thuật và ứng dụng thực tiễn sâu sắc cho 4 nhóm đối tượng sau:
- Học viên cao học và nghiên cứu sinh chuyên ngành Cơ điện tử, Cơ kỹ thuật và Kỹ thuật Điều khiển - Tự động hóa: Cung cấp tài liệu tham khảo chuẩn mực về phương pháp rời rạc hóa Galerkin, lý thuyết dầm Euler-Bernoulli và kỹ thuật phân tích modal thực nghiệm trên phần mềm MATLAB.
- Kỹ sư R&D thiết kế robot và thiết bị tự động hóa chính xác: Nắm bắt giải pháp điều khiển feedforward bằng phương pháp định hình xung để triệt tiêu rung động cho các cơ cấu chuyển động nhanh mà không cần gia tăng độ dày kết cấu hay chi phí phần cứng cảm biến.
- Chuyên gia phân tích động lực học kết cấu trong ngành hàng không - vũ trụ: Vận dụng mô hình toán học dầm mềm 3 chiều có tính đến hiện tượng không dãn và phép biến đổi tách kênh dao động để tính toán kết cấu cánh máy bay, tấm pin năng lượng mặt trời vệ tinh.
- Giảng viên và nhà nghiên cứu tại các trường đại học kỹ thuật: Sử dụng làm bài giảng chuyên đề hoặc phát triển các đề tài nghiên cứu chuyên sâu về điều khiển bền vững cho hệ thống tham số phân bố (Distributed Parameter Systems).
Câu hỏi thường gặp
Phương pháp Input Shaping có ưu điểm gì vượt trội so với điều khiển phản hồi kín truyền thống trong bài toán dầm mềm?
Input Shaping là kỹ thuật điều khiển vòng hở (feedforward) can thiệp trực tiếp vào quỹ đạo lệnh phát động. Phương pháp này loại bỏ hoàn toàn nhu cầu lắp đặt thêm cảm biến đo biến dạng hoặc cơ cấu chấp hành áp điện đắt tiền dọc thân dầm, đồng thời triệt tiêu nguy cơ mất ổn định vòng kín (closed-loop spillover instability) do các mode dao động bậc cao chưa được mô hình hóa gây ra.
Tại sao bộ định hình xung đơn mode không thể đưa dao động dư về mức dưới 2%?
Trong kết cấu dầm mềm dài 0,9 m, năng lượng dao động phân bố trên nhiều mode tần số. Bộ định hình đơn mode chỉ triệt tiêu năng lượng tại tần số riêng thứ nhất ($1,34\text{ Hz}$), trong khi mode dao động thứ hai ($8,4\text{ Hz}$) vẫn bị kích động trong quá trình bệ gá tăng tốc, khiến mức dao động dư còn lại duy trì ở khoảng 4,5% đến 5,0%.
Tại sao phương pháp quét tần số hàm mũ lại tối ưu hơn quét tần số tuyến tính khi nhận dạng modal?
Phương pháp quét tần số hàm mũ (Exponential Sine Sweep) phân bổ thời gian lưu lại ở dải tần số thấp lâu hơn so với quét tuyến tính. Điều này cung cấp đủ năng lượng kích thích cho mode dao động riêng đầu tiên ($1,34\text{ Hz}$) vốn nằm sát giới hạn dưới của dải đo, giúp hàm truyền phản hồi tần số FRF thu được có biên độ đỉnh sắc nét và thuật toán khớp hàm đạt độ chính xác cao.
Mô hình giảm chấn Rayleigh ảnh hưởng như thế nào đến độ tin cậy của mô phỏng số?
Mô hình Rayleigh ($C = a_0 M + a_1 K$) kết hợp tuyến tính giữa ma trận khối lượng và ma trận độ cứng, cho phép bảo toàn tính trực giao của các hàm dạng mode. Điều này giúp tách rời hệ phương trình vi phân nhiều bậc tự do mà vẫn phản ánh chính xác tốc độ suy giảm dao động thực tế với hệ số cản $\zeta \approx 0,0075$, đảm bảo độ tương đồng cao giữa mô phỏng và thực nghiệm.
Làm cách nào để tách rời phương trình chuyển động của thanh dầm mềm trong không gian 3 chiều?
Bằng cách đặt gốc tọa độ tại tâm diện tích mặt cắt ngang và thực hiện phép quay hệ trục tọa độ một góc xác định để trùng với các trục uốn chính của tiết diện dầm, tích quán tính ly tâm sẽ bị triệt tiêu bằng 0. Khi kết hợp với giả thiết dầm không dãn, phương trình dao động 3D phức tạp được phân rã thành ba phương trình độc lập theo hai phương uốn ngang và một phương dọc trục.
Kết luận
- Luận văn đã xây dựng thành công mô hình động lực học rời rạc có cản Rayleigh cho cơ cấu thanh dầm mềm công xôn 1D và mô hình hóa hệ thống dầm 3D tách kênh trong không gian.
- Nhận dạng chính xác hai tần số dao động riêng của kết cấu thực nghiệm là $1,34\text{ Hz}$ và $8,4\text{ Hz}$ cùng tỷ số cản tương ứng $\zeta \approx 0,0075$ thông qua kỹ thuật quét tần số hàm mũ và phân tích modal thực nghiệm.
- Chứng minh tính ưu việt của bộ định hình xung đa mode Two-mode ZVDD khi triệt tiêu dao động dư xuống mức ấn tượng 0,3% và độ vọt lố dưới 0,02% so với trạng thái không điều khiển.
- Kiểm chứng thực nghiệm trên bệ thử động cơ tuyến tính Yokogawa và camera Keyence cho thấy sự trùng khớp chặt chẽ giữa tính toán số và hành vi thực tế của cơ cấu.
- Kế hoạch tiếp theo tập trung vào việc hiện thực hóa thuật toán tự chỉnh định trên phần cứng nhúng và mở rộng thực nghiệm cho hệ thống thanh dầm không gian 3D.
Quý độc giả, nhà nghiên cứu và kỹ sư quan tâm có thể khai thác chi tiết các phương trình toán học, thuật toán mô phỏng MATLAB và bản vẽ thiết kế thực nghiệm trong toàn văn luận văn thạc sĩ để ứng dụng trực tiếp vào các dự án cơ điện tử và tự động hóa chính xác.