Giới thiệu dự án

Trong bối cảnh nền công nghiệp 4.0 đang thúc đẩy quá trình tự động hóa và hiện đại hóa sản xuất toàn cầu, ngành gia công cơ khí chính xác đóng vai trò then chốt trong chuỗi cung ứng công nghiệp phụ trợ. Theo các khảo sát công nghiệp chế tạo, việc ứng dụng máy công cụ điều khiển số (Computer Numerical Control - CNC) giúp nâng cao năng suất từ 35% đến 50% và giảm tỷ lệ phế phẩm xuống dưới 1.5% so với phương pháp gia công cơ khí truyền thống. Tuy nhiên, các hệ thống máy phay cơ khí cũ hoặc máy điều khiển vòng hở (open-loop) thường gặp hạn chế nghiêm trọng về độ cứng vững động học, sai số do ma sát trượt và hiện tượng mất bước khi gia công tải nặng, dẫn đến giảm độ chính xác kích cỡ và độ nhám bề mặt chi tiết gia công.

Đồ án "Thiết kế hệ thống điều khiển CNC cho bàn máy phay 3 trục" do sinh viên Nguyễn Anh Tuấn (MSSV: 20184654, Lớp CĐT 02 – K63, Trường Cơ khí - Đại học Bách Khoa Hà Nội) thực hiện dưới sự hướng dẫn của TS. Lê Đức Độ, mã học phần ME5512 (Mã đề VCK02-9), giải quyết trực tiếp bài toán thiết kế hệ thống điều khiển truyền động vị trí vòng kín cho bàn máy phay CNC 3 trục chịu tải trọng lớn.

+-------------------------------------------------------------------------+
|                        MỤC TIÊU CỐT LÕI CỦA DỰ ÁN                       |
+-------------------------------------------------------------------------+
| 1. Xây dựng mô hình toán học động học/động lực học bàn máy X-Y          |
| 2. Tính toán độ cứng vững tương đương (K = 182x10^3 N/m)                |
| 3. Thiết kế bộ điều khiển PID tối ưu hóa theo Ziegler-Nichols 1         |
| 4. Khảo sát ổn định hệ thống qua tiêu chuẩn Routh, Nyquist, Bode        |
| 5. Tích hợp phần cứng PLC S7-200, Servo 400W và Encoder vòng kín        |
+-------------------------------------------------------------------------+

Dự án đặt ra các chỉ tiêu kỹ thuật khắt khe: gia công phôi thép hợp kim JIS S45C (Grade 4040) có khối lượng tối đa $M_p = 700\text{ kg}$, vận tốc chạy dao không tải $V_1 = 18\text{ m/ph}$, vận tốc cắt gọt có tải $V_2 = 12\text{ m/ph}$, gia tốc cực đại $a = 5\text{ m/s}^2$, chế độ cắt phay mặt đầu 6 lưỡi cắt đường kính $D = 80\text{ mm}$, chiều sâu cắt $t = 1.2\text{ mm}$, lượng tiến dao $F = 900\text{ mm/ph}$ và vận tốc cắt $v = 100\text{ m/ph}$.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trước khi tiến hành thiết kế, việc đánh giá các giải pháp truyền động và điều khiển máy phay trên thị trường được tiến hành nhằm xác định cấu trúc tối ưu:

Tiêu chí so sánh Hệ thống vòng hở (Stepper Motor) Hệ thống vòng kín công nghiệp (CNC Closed-Loop) Giải pháp đề xuất của đồ án
Cơ cấu phản hồi Không có cảm biến phản hồi vị trí Thước quang tuyến tính hoặc Encoder tuyệt đối Encoder Omron tương đối gắn trục Servo + Vít me bi
Độ chính xác Thấp ($0.05 - 0.1\text{ mm}$), dễ mất bước Rất cao ($< 0.005\text{ mm}$) Cao ($< 0.01\text{ mm}$), kiểm soát phản hồi kín
Khả năng bù tải Kém, sụt bước khi lực cắt $F_m$ biến thiên Tự động thích nghi thời gian thực Khâu tích phân (Ki) triệt tiêu sai số xác lập
Chi phí triển khai Rất thấp Rất cao (yêu cầu bộ điều khiển CNC chuyên dụng) Vừa phải, module hóa cao (PLC S7-200 + Driver rời)
Bảo trì, sửa chữa Đơn giản Phức tạp, phụ thuộc nhà sản xuất Dễ dàng chẩn đoán và thay thế linh kiện

Yêu cầu người dùng được phân loại theo mô hình MoSCoW:

  • Must Have: Khả năng chịu tải phôi $700\text{ kg}$, cơ chế dừng khẩn cấp và giới hạn hành trình bằng công tắc hành trình cơ điện, chu trình về gốc máy (Home Return).
  • Should Have: Thời gian xác lập $T_s < 0.2\text{ s}$, độ quá điều chỉnh $\Delta h < 5%$.
  • Could Have: Giao diện điều khiển liên kết màn hình HMI.
  • Won't Have: Bù sai số giãn nở nhiệt của trục chính trong phạm vi đồ án môn học.

Thiết kế hệ thống

Kiến trúc điều khiển được phân tầng rõ ràng từ tầng giao tiếp, tầng điều khiển logic đến tầng chấp hành cơ điện:

[ G-code / CAD/CAM / Input Switch ]
[ Bộ điều khiển logic khả trình PLC Siemens S7-200 (CPU 224/226) ]
    (Xung PTO/PWM & Hướng)
                                                                [ Bàn máy X/Y + Phôi 700kg ]

Technology Stack và thông số linh kiện được chuẩn hóa:

  • Bộ điều khiển lập trình: Siemens S7-200 CPU (24 Digital Inputs, 16 Digital Outputs, 6 bộ đếm tốc độ cao HSC 32-bit, 2 bộ phát xung tốc độ cao PTO/PWM tần số $20\text{ kHz}$).
  • Động cơ chấp hành: AC Servo Motor model BPTA004GBA6D1, công suất định mức $400\text{ W}$, điện áp $220\text{ VAC}$, dòng định mức $2.5\text{ A}$, mô-men xoắn định mức $1.27\text{ Nm}$ (mô-men xoắn cực đại $3.82\text{ Nm}$), tốc độ quay định mức $3000\text{ rpm}$ (max $4500\text{ rpm}$), quán tính rotor $J_m = 0.28\text{ kg}\cdot\text{cm}^2$.
  • Cảm biến vị trí: Incremental Encoder Omron E6B2-CWZ6C (3 pha A, B, Z, ngõ ra NPN Open Collector, tần số đáp ứng $100\text{ kHz}$, cấp bảo vệ IEC 60529 IP50).
  • Thiết bị an toàn: Công tắc hành trình Omron Z-15GQ-B (tiếp điểm SPDT 15A 125VAC, độ bền cơ học 10 triệu chu kỳ, tốc độ tác động $1\text{ mm/s} \div 1\text{ m/s}$) và nút ấn công nghiệp LA38-11BN.

Mô hình toán học của hệ thống truyền động bàn máy được thiết lập thông qua phương trình vi phân cân bằng lực:

$$m \ddot{x}(t) + B \dot{x}(t) + K x(t) = F_{dc}(t) - F_{ms}(t)$$

Trong đó:

  • $m$: Tổng khối lượng chuyển động ($m_X = M_p + M_x = 700 + 140 = 840\text{ kg}$; $m_Y = M_p + M_x + M_y = 700 + 140 + 250 = 1090\text{ kg}$).
  • $l = 10\text{ mm} = 0.01\text{ m}$: Bước vít me bi.
  • Lực kéo động cơ: $F_{dc}(t) = K \frac{l}{2\pi} \theta(t)$.
  • Lực ma sát: $F_{ms}(t) = f m g \dot{x}(t)$ với hệ số ma sát trượt ray dẫn hướng $f = 0.1$.
  • Độ cứng vững tổng hợp $K$ được tính toán từ chuỗi độ cứng nối tiếp:

$$\frac{1}{K} = \frac{1}{k_c} + \frac{1}{k_b} + \frac{1}{k_g + k_s} + \frac{1}{k_{ch}} + \frac{1}{k_n}$$

Với các thông số thực tế: độ cứng nối trục đàn hồi Daidoseimitu $k_c = 205.9 \times 10^4\text{ N/m}$, ổ bi đỡ chặn SKF $k_b = 120 \times 10^6\text{ N/m}$, ray dẫn hướng PMI $k_g = 3 \times 10^9\text{ N/m}$, vít me bi $k_s = 23.86 \times 10^7\text{ N/m}$, càng gá $k_{ch} = 2 \times 10^5\text{ N/m}$, đai ốc bi $k_n = 108.33 \times 10^7\text{ N/m}$. Giá trị độ cứng tương đương tính được là $K = 182 \times 10^3\text{ N/m}$.

Hệ số cản nhớt tương đương với hệ số tắt dần $\xi = 0.35$:

$$B = 2 \xi \sqrt{m K}$$

Hàm truyền đạt vòng hở của bàn máy trục X qua phép biến đổi Laplace:

$$G_X(s) = \frac{X(s)}{\theta(s)} = \frac{K \frac{l}{2\pi}}{m s^2 + (B + f m g) s + K} = \frac{289.7}{840 s^2 + 8964 s + 182000}$$

Tương tự đối với bàn máy trục Y:

$$G_Y(s) = \frac{289.7}{1090 s^2 + 10949 s + 182000}$$

Methodology

Quy trình phát triển đồ án được tổ chức theo phương pháp luận Waterfall-V mô phỏng theo chuẩn kỹ thuật công nghiệp trong vòng 15 tuần:

  • Tuần 1 - 4: Thu thập tài liệu, phân tích động học, tính toán lực cắt $F_m$ và lựa chọn cơ cấu vít me bi, ray dẫn hướng, khớp nối.
  • Tuần 5 - 8: Xây dựng mô hình toán học giải tích, tính toán độ cứng tổng hợp, xác định hàm truyền đạt động cơ - bàn máy.
  • Tuần 9 - 11: Khảo sát ổn định (Routh, Nyquist, Bode) và thiết kế thông số bộ điều khiển PID bằng công cụ MATLAB/Simulink R2022a.
  • Tuần 12 - 14: Lập trình điều khiển logic PLC S7-200 (STEP 7-Micro/WIN v4.0 SP9), xây dựng thuật toán quét gốc Home và liên động an toàn.
  • Tuần 15: Đánh giá tổng hợp chỉ tiêu chất lượng, hoàn thiện thuyết minh và bản vẽ mạch điện điều khiển.

Implementation và kết quả

Development process

Quá trình tổng hợp bộ điều khiển PID được thực hiện theo phương pháp thực nghiệm Ziegler-Nichols thứ nhất kết hợp thuật toán tối ưu hóa đáp ứng trên nền tảng MATLAB/Simulink. Hàm truyền đạt của bộ điều khiển PID có cấu trúc:

$$R(s) = K_p \left( 1 + \frac{1}{T_i s} + T_d s \right) = K_p + \frac{K_i}{s} + K_d s$$

% Script tính toán và mô phỏng hệ thống điều khiển bàn X
clear; clc;
num = [289.7];
den = [840 8964 182000];
Gx = tf(num, den);

% Kiểm tra tính ổn định hệ hở
poles = pole(Gx);
disp('Cực của hệ thống hở:');
disp(poles);

% Thiết kế bộ điều khiển PID tối ưu
Kp = 1215.9832;
Ki = 9262.6158;
Kd = 39.5397;
C = pid(Kp, Ki, Kd);

% Hệ thống hồi tiếp âm đơn vị có PID
Sys_closed = feedback(C * Gx, 1);

% Phân tích đáp ứng bước nhảy
figure(1);
step(Sys_closed);
grid on;
title('Đáp ứng bước nhảy bàn X với bộ điều khiển PID');
stepinfo(Sys_closed)

Chương trình điều khiển logic chạy trên PLC Siemens S7-200 quản lý quy trình chạy máy và thuật toán quét gốc Home (Origin Return):

// Thuật toán quét gốc HOME trên PLC Siemens S7-200 (Structured Text / Statement List)
NETWORK 1: Kích hoạt chu trình Home
LD     I0.0          // Nút nhấn Start Home
AN     M0.1          // Cờ hoàn thành Home
S      M0.0, 1       // Bật cờ thực hiện chạy về gốc

NETWORK 2: Điều khiển phát xung lùi tìm công tắc hành trình Limit X-
LD     M0.0
AN     I0.2          // Công tắc hành trình Limit X- chưa tác động
CALL   PTO0_RUN, 5000 // Phát xung tần số 5kHz quay ngược chiều

NETWORK 3: Giảm tốc và dừng khi gặp cờ giới hạn, tìm xung Z của Encoder
LD     I0.2          // Đã chạm Limit X-
R      M0.0, 1
S      M0.2, 1       // Bật cờ tìm xung Zero pha Z
CALL   PTO0_RUN, 500  // Giảm tốc độ còn 500Hz

NETWORK 4: Xác lập tọa độ gốc máy
LD     M0.2
A      I0.5          // Tín hiệu pha Z từ Encoder E6B2-CWZ6C
R      M0.2, 1
S      M0.1, 1       // Xác nhận hoàn tất Home X
MOVD   0, VD100      // Reset thanh ghi vị trí hiện tại về 0

Testing và validation

Tính ổn định của hệ thống được kiểm chứng qua 4 tiêu chuẩn giải tích và miền tần số:

+-------------------------------------------------------------------------+
|                  KẾT QUẢ KHẢO SÁT ỔN ĐỊNH HỆ THỐNG TRỤC X               |
+-------------------------------------------------------------------------+
| 1. Tiêu chuẩn Routh: Đa thức mẫu số A(s) có toàn bộ hệ số dương        |
|    Nghiệm cực s1,2 = -5.336 ± 13.718j nằm hoàn toàn bên trái trục ảo    |
| 2. Tiêu chuẩn Nyquist: Đường cong Nyquist của Gh(s) không bao điểm -1+0j|
| 3. Tiêu chuẩn Bode: Biên độ và góc pha thỏa mãn -pi < phi < 0           |
| 4. Đáp ứng Dirac: Biên độ đạt cực đại 0.0147 tại t=0.0863s (< 0.5)      |
+-------------------------------------------------------------------------+

Kết quả đạt được

Việc bổ sung bộ điều khiển PID đã cải thiện vượt bậc các chỉ tiêu chất lượng động học của hệ thống truyền động bàn máy:

Chỉ tiêu chất lượng Hệ thống ban đầu (Chưa có PID) Hệ thống sau khi có PID Mức độ cải thiện (%)
Thời gian tăng ($T_r$) $0.0958\text{ s}$ $0.0713\text{ s}$ Nhanh hơn 25.57%
Thời gian xác lập ($T_s$) $0.7370\text{ s}$ $0.1300\text{ s}$ Giảm 82.36%
Độ quá điều chỉnh ($\Delta h$) $29.40%$ $1.00%$ Giảm 96.60% (triệt tiêu dao động)
Biên độ cực đại (Peak Amplitude) $0.00206$ (chưa bù độ lợi) $1.01000$ (bám chính xác giá trị đặt) Chuẩn hóa đáp ứng vị trí 100%
Sai số xác lập ($e_{xl}$) Lớn do ảnh hưởng lực ma sát $F_{ms}$ Xấp xỉ $0\text{ mm}$ Triệt tiêu hoàn toàn nhờ khâu $K_i$

Các thông số điều khiển PID tối ưu thu được cho trục X:

  • Hệ số khuếch đại tỷ lệ: $K_p = 1215.9832$
  • Hệ số tích phân: $K_i = 9262.6158$
  • Hệ số vi phân: $K_d = 39.5397$

Đổi mới và đóng góp

  1. Mô hình hóa độ cứng tương đương đa thành phần: Khác với các mô hình đơn giản hóa chỉ xét độ cứng vít me, đồ án đã giải tích hóa toàn diện chuỗi độ cứng bao gồm khớp nối đàn hồi Daidoseimitu ($k_c$), ổ lăn đỡ chặn SKF ($k_b$), ray dẫn hướng PMI ($k_g$), trục vít me bi ($k_s$), gá đỡ ($k_{ch}$) và đai ốc bi ($k_n$), đạt giá trị tổng hợp chính xác $K = 182 \times 10^3\text{ N/m}$.
  2. Triệt tiêu hoàn toàn sai số ma sát tĩnh và ma sát trượt: Tích hợp thành phần lực ma sát $F_{ms} = f m g \dot{x}$ trực tiếp vào hàm truyền đạt miền phức $s$, giúp khâu tích phân $K_i$ bù đắp tức thời lượng sụt áp mô-men khi máy bắt đầu chuyển động.
  3. Giảm thiểu thời gian quá độ tới 82.36%: Tối ưu hóa bộ thông số PID theo phương pháp Ziegler-Nichols thứ nhất giúp hệ thống bàn máy nặng gần 1 tấn đạt trạng thái ổn định chỉ trong $0.13\text{ s}$, đáp ứng hoàn hảo yêu cầu nội suy đường cong biên dạng phức tạp trong gia công CNC.

Ứng dụng thực tế và triển khai

Giải pháp trong đồ án có khả năng ứng dụng rộng rãi vào các lĩnh vực chế tạo và cải tạo máy công cụ:

  • Nâng cấp (Retrofit) máy công cụ cũ: Ứng dụng gói điều khiển PLC S7-200 kết hợp AC Servo và Encoder Omron để chuyển đổi các máy phay cơ truyền thống (vạn năng) thành máy phay CNC 3 trục bán tự động với chi phí chỉ bằng 25-30% so với mua máy mới.
  • Dây chuyền gia công chi tiết ô tô, xe máy: Chịu tải phôi tĩnh lên tới $700\text{ kg}$, phù hợp gia công các chi tiết vỏ hộp số, khuôn dập, đồ gá kích thước lớn từ vật liệu thép S45C, gang xám.
  • Yêu cầu môi trường triển khai: Tủ điều khiển đạt chuẩn IP54 chống bụi kim loại và dầu làm nguội, hệ thống tiếp địa độc lập cho Driver Servo để lọc nhiễu sóng hài tần số cao tác động lên tín hiệu Encoder.
+-------------------------------------------------------------------------+
|                  ƯỚC TÍNH CHI PHÍ VÀ LỢI ÍCH KINH TẾ (ROI)              |
+-------------------------------------------------------------------------+
| • Chi phí phần cứng điều khiển (PLC, 3 bộ Servo 400W, Encoder, cảm biến)|
|   Dao động: 45.000.000 - 60.000.000 VNĐ / cụm bàn máy 3 trục            |
| • Năng suất gia công tăng: +40% nhờ tốc độ chạy không tải 18 m/ph       |
| • Thời gian hoàn vốn (ROI): Ước tính 6 - 9 tháng hoạt động 2 ca/ngày     |
+-------------------------------------------------------------------------+

Hạn chế và hướng phát triển

  • Hạn chế kỹ thuật hiện tại:
    • Mô hình toán học chưa tính đến sự phi tuyến hóa của hiện tượng khe hở cơ khí (backlash) khi đảo chiều chuyển động vít me bi sau thời gian dài vận hành mòn cơ học.
    • Bộ điều khiển PID tuyến tính với hệ số cố định chưa tự thích nghi khi khối lượng phôi $M_p$ thay đổi đột ngột từ $0\text{ kg}$ lên $700\text{ kg}$.
  • Hướng phát triển tiếp theo:
    • Nghiên cứu ứng dụng bộ điều khiển mờ thích nghi (Adaptive Fuzzy-PID) hoặc điều khiển trượt (Sliding Mode Control) để tự động căn chỉnh thông số theo tải trọng thực.
    • Mở rộng giao thức truyền thông công nghiệp thời gian thực (như Modbus RTU hoặc PROFINET) kết nối đồng bộ 3 trục X-Y-Z và trục chính Spindle với phần mềm máy tính trung tâm (Mach3/LinuxCNC).

Đối tượng hưởng lợi

  • Sinh viên Cơ điện tử và Kỹ thuật Cơ khí: Cung cấp tài liệu mẫu mực về phương pháp mô hình hóa cơ điện, tính toán độ cứng cơ học và thiết kế hệ thống điều khiển tự động vòng kín từ lý thuyết đến mô phỏng.
  • Kỹ sư R&D và Chế tạo máy: Tham khảo bộ thông số phần cứng thực tế (Servo BPTA004GBA6D1, PLC S7-200, Encoder E6B2-CWZ6C) và thuật toán giải tích để rút ngắn thời gian thiết kế máy CNC công nghiệp.
  • Doanh nghiệp vừa và nhỏ (SMEs): Tiếp cận giải pháp tự động hóa chi phí thấp, tối ưu hóa năng lực sản xuất từ các hệ thống máy gia công sẵn có.

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

  1. Yêu cầu kỹ thuật phần cứng tối thiểu để triển khai hệ thống là gì? Hệ thống yêu cầu nguồn cấp 3 pha 220V/380V ổn định, bộ điều khiển PLC hỗ trợ tối thiểu 2 kênh phát xung tốc độ cao PTO tần số từ $20\text{ kHz}$ trở lên, các kênh ngắt tốc độ cao HSC đọc tín hiệu Encoder vi sai hoặc NPN, và động cơ Servo có mô-men xoắn giữ tối thiểu $1.27\text{ Nm}$.

  2. Hệ thống xử lý bài toán quá tải trọng khi phôi vượt quá 700kg như thế nào? Khi khối lượng phôi vượt định mức, quán tính hệ thống tăng làm thay đổi hàm truyền đạt $G(s)$, dẫn đến tăng thời gian xác lập $T_s$. Driver Servo tích hợp cảnh báo quá dòng/quá tải (Overload Alarm) sẽ tự động kích hoạt ngắt an toàn và dừng phát xung PLC để bảo vệ cơ cấu vít me bi và cuộn dây động cơ.

  3. Làm thế nào để tích hợp hệ thống điều khiển này với phần mềm CAD/CAM? Mã G-code xuất ra từ phần mềm CAM (như Mastercam, Fusion 360) được bộ vi xử lý trung tâm biên dịch thành tập lệnh vị trí và tọa độ, sau đó truyền thông qua cổng RS485/PPI xuống khối nhớ dữ liệu của PLC S7-200 để thực thi phát chuỗi xung điều khiển từng trục.

  4. Bảo trì và hiệu chuẩn hệ thống định kỳ bao gồm những công việc gì? Cần tra mỡ bôi trơn chuyên dụng định kỳ cho trục vít me bi và ray dẫn hướng PMI sau mỗi 500 giờ làm việc, kiểm tra độ căng của khớp nối đàn hồi Daidoseimitu và kiểm tra xung lệch Zero pha Z của Encoder để đảm bảo độ chính xác vị trí bàn máy.

  5. Thời gian thu hồi vốn (ROI) khi nâng cấp máy phay bằng giải pháp này là bao lâu? Với chi phí đầu tư phần cứng và lắp đặt khoảng 50 - 60 triệu VNĐ, hệ thống giúp tiết kiệm 70% thời gian gá đặt và tăng tốc độ gia công, thời gian thu hồi vốn trung bình ước tính từ 6 đến 9 tháng tùy thuộc vào công suất vận hành của xưởng sản xuất.

Kết luận

Đồ án môn học "Thiết kế hệ thống điều khiển CNC cho bàn máy phay 3 trục" đã giải quyết triệt để và đồng bộ từ bài toán cơ sở động lực học cơ khí, mô hình hóa giải tích độ cứng tương đương ($K = 182 \times 10^3\text{ N/m}$), khảo sát tính ổn định hệ thống qua các tiêu chuẩn Routh, Nyquist, Bode, đến việc thiết kế bộ điều khiển PID tối ưu trên MATLAB/Simulink và tích hợp phần cứng công nghiệp PLC S7-200, AC Servo 400W. Kết quả nghiên cứu đạt độ tin cậy cao với thời gian xác lập cực nhanh $0.13\text{ s}$, độ vọt lố triệt tiêu về $1%$, tạo nền tảng ứng dụng vững chắc cho việc phát triển và chế tạo các hệ máy CNC chuyên dụng trong công nghiệp cơ khí Việt Nam.