Tổng quan về luận án

Trong kỷ nguyên Cách mạng Công nghiệp 4.0 và chuyển đổi số quy trình sản xuất, các hệ thống điều khiển truyền động chính xác đóng vai trò cốt lõi trong dây chuyền gia công CNC, cánh tay robot công nghiệp và cơ cấu tự động hóa phức tạp. Trọng tâm của các hệ thống này là động cơ servo – thiết bị chấp hành đòi hỏi mômen xoắn riêng lớn, khả năng quá tải cao, thời gian đáp ứng động lực học cực nhanh và dải điều chỉnh tốc độ rộng từ vùng mômen không đổi đến vùng suy giảm từ thông. Động cơ servo không đồng bộ (Induction Servo Motor - ISM) 3 pha rotor lồng sóc sở hữu độ bền cơ học vượt trội, cấu trúc không sử dụng chổi than hay nam châm đất hiếm đắt đỏ, cùng khả năng vận hành tin cậy trong môi trường công nghiệp khắc nghiệt. Tuy nhiên, việc thiết kế tối ưu hóa loại động cơ này gặp phải rào cản lớn do tính chất phi tuyến phức tạp của mạch từ, hiện tượng bão hòa cục bộ, tổn hao sóng hài bậc cao và các ràng buộc nhiệt nghiêm ngặt ở các chế độ làm việc chu kỳ khác nhau ($S_1$ đến $S_9$).

Khoảng trống nghiên cứu (research gap) lớn nhất hiện nay xuất phát từ việc các phương pháp thiết kế truyền thống (Alger, 1965; Kostenko & Piotrovsky, 1977) chủ yếu tập trung vào điểm định mức đơn lẻ ở chế độ làm việc liên tục ($S_1$) với tần số công nghiệp cố định 50 Hz. Trong khi đó, động cơ servo vận hành với biến tần điều khiển trường (Field-Oriented Control) trên toàn dải đặc tính mômen – tốc độ, bao gồm vùng quá tải ngắn hạn ($S_2$) và các chu kỳ tăng giảm tốc liên tục ($S_3$–$S_8$). Hơn nữa, việc tối ưu hóa độc lập từng mục tiêu riêng lẻ (chỉ giảm khối lượng hoặc chỉ tăng hiệu suất) thường dẫn đến sự suy giảm nghiêm trọng về đáp ứng động học hoặc quá nhiệt cục bộ. Luận án tiến sĩ kỹ thuật của NCS. Nguyễn Đức Bắc, dưới sự hướng dẫn của TS. Nguyễn Thế Công và TS. Trần Tuấn Vũ tại Trường Đại học Bách khoa Hà Nội (2022), đã giải quyết triệt để khoảng trống này bằng công trình: "Nghiên cứu thiết kế tối ưu động cơ servo không đồng bộ 3 pha rotor lồng sóc" (Mã số chuyên ngành: 9520201 - Kỹ thuật điện).

Nghiên cứu tập trung giải quyết ba câu hỏi nghiên cứu cốt lõi:

  1. $RQ_1$: Làm thế nào để xây dựng một mô hình toán học giải tích tích hợp đồng thời các quan hệ phi tuyến điện từ, thông số hình học và giới hạn nhiệt nhằm thiết lập bài toán tối ưu hóa đa mục tiêu cho động cơ servo không đồng bộ?
  2. $RQ_2$: Thuật toán tối ưu hóa xác định nào có khả năng hội tụ nhanh, chính xác và khắc phục được nhược điểm tính toán nặng nề của các thuật toán tiến hóa ngẫu nhiên trong việc tìm kiếm biên Pareto tối ưu giữa khối lượng điện từ và tổn hao năng lượng?
  3. $RQ_3$: Tác động tương hỗ giữa kết cấu hình học rãnh rotor (rãnh nghiêng so với rãnh thẳng) và cơ chế tản nhiệt (vỏ kín so với nắp hở) ảnh hưởng như thế nào đến độ nhấp nhô mômen và độ tăng nhiệt bão hòa của động cơ ở các điểm làm việc then chốt?

Từ các câu hỏi trên, ba giả thuyết nghiên cứu được kiểm chứng:

  • $H_1$: Tích hợp quy hoạch đa thức bậc hai tuần tự (Sequential Quadratic Programming - SQP) với phương pháp ràng buộc $\varepsilon$ ($\varepsilon$-constraint) sẽ cho phép xác định chính xác tập nghiệm Pareto biên, giảm thiểu số lần lặp tính toán và triệt tiêu chi phí chế tạo thử nghiệm tốn kém.
  • $H_2$: Cấu trúc rotor rãnh nghiêng tối ưu hình học sẽ làm suy giảm đáng kể các sóng hài không gian của sức điện động ($e_{ssst}$), từ đó hạ thấp độ nhấp nhô mômen (torque ripple) mà không làm suy giảm mômen cực đại yêu cầu ($T_m \ge 15\text{ N.m}$).
  • $H_3$: Mô hình mạng nhiệt tập trung tương đương (Lumped Parameter Thermal Network) kết hợp phân tích phần tử hữu hạn 2D/3D (FEA) phản ánh chính xác trường nhiệt độ thực nghiệm với sai số không vượt quá 5% ở cả chế độ liên tục ($S_1$) và ngắn hạn ($S_2$).

Về phạm vi và ý nghĩa, luận án đã tính toán, thiết kế tối ưu, mô phỏng đa vật lý và chế tạo mẫu thử nghiệm hoàn chỉnh trên băng thử công suất chuyên dụng (test-bench), kiểm chứng tại các dải tốc độ $n = 500\text{ rpm}, 700\text{ rpm}, 1000\text{ rpm}$ đến $3500\text{ rpm}$, điện áp dây từ $300\text{V}$ đến $380\text{V}$, và tần số điều khiển từ $30\text{ Hz}$ đến $60\text{ Hz}$.

Literature Review và Positioning

Lịch sử nghiên cứu thiết kế máy điện đã chứng kiến sự chuyển dịch từ các công thức kinh nghiệm giải tích cổ điển sang mô hình hóa tối ưu hóa có sự trợ giúp của máy tính (CAD/CAE). Trong phân nhóm tối ưu hóa động cơ đồng bộ nam châm vĩnh cửu (PMSM), Svante Andersson (2000) đã ứng dụng thuật toán SQP để giảm thiểu chi phí vật liệu chế tạo rotor 6 cực gắn nam châm bề mặt. Tiếp đó, Damir Zarko, Drago Ban và Davor Gooricki (2009) đã áp dụng giải thuật tiến hóa vi phân (Differential Evolution - DE) kết hợp phân tích phần tử hữu hạn 2D (FEA) để cải tiến kích thước rotor và hình học nam châm cho động cơ servo AC, chứng minh sự suy giảm đáng kể của mômen quán tính và độ tăng nhiệt. Đối với phân tích nhiệt, David A. Staton (2001) đã tiên phong phát triển phương pháp mạng nhiệt tập trung thông qua nền tảng Motor-CAD, cho phép giảm kích thước khung vỏ động cơ servo PM nhưng vẫn đảm bảo ngưỡng nhiệt độ giới hạn theo cấp cách điện tiêu chuẩn IEC.

Trong phân nhóm động cơ không đồng bộ (Induction Motor), Mehmet Çunkaş và Ramazan Akkaya (2006) đã triển khai thuật toán di truyền (Genetic Algorithm - GA) cho bài toán tối ưu hóa đơn mục tiêu riêng biệt: cực đại hóa mômen, tối thiểu hóa chi phí chế tạo, hoặc tối đa hóa hiệu suất. Tuy nhiên, nghiên cứu của Çunkaş và Akkaya bộc lộ nhược điểm lớn khi tách rời các hàm mục tiêu, không phản ánh được sự đánh đổi (trade-off) tự nhiên giữa khối lượng đồng/sắt và tổn hao năng lượng trong các điều kiện tải biến đổi. Tại phân khúc tốc độ cao, Pil-Wan Han và các cộng sự (2015) đã tiến hành phân tích liên kết điện từ - nhiệt cho động cơ không đồng bộ $12.000\text{ rpm}$ bằng sơ đồ mạch thay thế nhiệt kết hợp FEA. Jean Le Besnerais cùng nhóm tác giả (2009) tập trung tối ưu hóa độ mở miệng rãnh stator và rotor nhằm triệt tiêu sóng hài từ trường gây rung ồn.

                    TIẾN TRÌNH NGHIÊN CỨU THIẾT KẾ TỐI ƯU MÁY ĐIỆN
[Nghiên cứu Cổ điển]    Alger (1965), Kostenko (1977): Tính toán giải tích đơn điểm S1
[Tối ưu Đơn mục tiêu]   Andersson (2000) [SQP - PMSM]; Çunkaş & Akkaya (2006) [GA - Induction]
[Mô phỏng Đa vật lý]    Staton (2001) [Motor-CAD Thermal]; Pil-Wan Han et al. (2015) [Coupled FEA]
[KHOẢNG TRỐNG TỒN TẠI]  Thiếu phương pháp tối ưu hóa đa mục tiêu xác định (SQP + ε-constraint)
                        kết hợp mô hình nhiệt động học S1/S2 và kiểm chứng thực nghiệm băng thử
[ĐÓNG GÓP CỦA LUẬN ÁN]  Nguyễn Đức Bắc (2022): Tối ưu đa mục tiêu Pareto (W_tot vs 1-η), FEA 3D
                        rãnh nghiêng, mô phỏng nhiệt liên kết và chế tạo mẫu thử nghiệm thực tế

Trong bức tranh tổng quan đó, tồn tại hai luồng quan điểm đối lập sâu sắc:

  1. Trường phái Heuristic/Meta-heuristic (GA, PSO, DE): Cho rằng các giải thuật ngẫu nhiên toàn cục là phương thức duy nhất tránh được bẫy cực trị địa phương trong không gian biến số máy điện phi tuyến phức tạp.
  2. Trường phái Xác định/Gradient-based (SQP, Newton-Raphson): Khẳng định các giải thuật dựa trên gradient bậc hai có tốc độ hội tụ nhanh hơn hàng bậc độ lớn, độ chính xác số tuyệt đối và hoàn toàn xác định nếu điểm khởi tạo được giới hạn trong miền khả thi vật lý thông qua phân tích giải tích ban đầu.

Luận án của NCS. Nguyễn Đức Bắc định vị vững chắc tại giao điểm giữa tối ưu hóa xác định đa mục tiêu và mô phỏng đa vật lý kiểm chứng thực nghiệm. Khác với nghiên cứu của Çunkaş & Akkaya (2006) vốn chỉ dừng lại ở các bài toán đơn mục tiêu lý thuyết trên máy tính, và vượt lên trên nghiên cứu của Damir Zarko et al. (2009) vốn chỉ giới hạn trên động cơ đồng bộ nam châm vĩnh cửu, luận án đã thiết lập khung phân tích toàn diện cho động cơ servo không đồng bộ rotor lồng sóc – một đối tượng có mật độ phi tuyến từ trường và nhiệt độ phức tạp hơn nhiều do tổn hao đồng trên rotor ($p_{Al}$) phụ thuộc trực tiếp vào tần số trượt ($s$).

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng lý thuyết tính toán máy điện cổ điển sang lý thuyết tối ưu hóa kết cấu phi tuyến đa tham số bằng việc tích hợp thuật toán quy hoạch đa thức bậc hai tuần tự SQP (Sequential Quadratic Programming) với kỹ thuật xử lý đa mục tiêu $\varepsilon$-constraint. Thay vì kết hợp các hàm mục tiêu thành một hàm mục tiêu đơn có trọng số mang tính chủ quan, phương pháp $\varepsilon$-constraint giữ lại hàm mục tiêu thứ nhất là tối thiểu hóa khối lượng điện từ $f_1(x) = W_{tot}$ và chuyển hàm mục tiêu thứ hai là tối thiểu hóa tổn hao năng lượng (hoặc tối đa hóa hiệu suất $f_2(x) = 1 - \eta$) thành một ràng buộc bất đẳng thức phụ thuộc tham số $\varepsilon$:

$$\min f_1(x) = W_{tot}(x)$$ $$\text{thỏa mãn: } f_2(x) \le \varepsilon_k, \quad g_j(x) \le 0, \quad x_{min} \le x \le x_{max}$$

Trong đó, tập biến tối ưu hóa không gian thiết kế bao gồm: $$x = [L, D, s, U]^T$$ với $L$ là chiều dài tác dụng mạch từ ($30\text{ mm} \le L \le 180\text{ mm}$), $D$ là đường kính trong stator ($40\text{ mm} \le D \le 150\text{ mm}$), $s$ là tần số trượt ($1 \le s \le 25$), và $U$ là điện áp cấp ($70\text{V} \le U \le 380\text{V}$).

Các mệnh đề lý thuyết cốt lõi được xác lập trong luận án:

  • Mệnh đề 1 ($P_1$): Biên Pareto tạo bởi tập nghiệm SQP biểu diễn sự đánh đổi tuyệt đối giữa khối lượng vật liệu tác dụng và hiệu suất vận hành, trong đó mật độ từ thông bão hòa răng/gông ($B_{rt}, B_c$) và mật độ dòng điện dây quấn ($J$) đóng vai trò là các biến hoạt động trên biên tích cực của điều kiện tối ưu Kuhn-Tucker.
  • Mệnh đề 2 ($P_2$): Độ nghiêng rãnh rotor $\xi$ tạo ra hệ số suy giảm sóng hài sức điện động $k_{sq} = \frac{\sin(\nu \xi \pi / 2\tau)}{\nu \xi \pi / 2\tau}$, làm triệt tiêu thành phần mômen ký sinh bậc cao nhưng làm gia tăng điện cảm tản rãnh rotor, đòi hỏi sự phối hợp chính xác với tần số trượt $s$ để duy trì mômen cực đại.
  • Mệnh đề 3 ($P_3$): Nhiệt trở tương đương của lớp cách điện rãnh ($R_{cd}$) và hệ số tản nhiệt bề mặt đối lưu ($\alpha$) quyết định hằng số thời gian phát nóng $T$, tạo ra sự phân hóa rõ rệt về giới hạn dòng điện stator giữa chế độ làm việc liên tục $S_1$ và chế độ ngắn hạn $S_2$.
                KHUNG PHÂN TÍCH TỐI ƯU HÓA ĐA VẬT LÝ TOÀN DIỆN

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp ba trụ cột lý thuyết chuyên sâu:

  1. Lý thuyết trường điện từ và mạch thay thế phức: Thiết lập phương trình cân bằng áp và dòng trên hệ tọa độ phức, trong đó tổng trở trung gian $Z_3 = \frac{Z_2 X_m}{Z_2 + X_m}$ với $Z_1 = \frac{R_r}{s} + j X_r$, từ đó xác định chính xác dòng điện stator $I_s = \frac{U}{\sqrt{3}|Z_{tot}|}$, công suất truyền qua khe hở không khí $P_{tr} = 3 \frac{R_r}{s} I_r^2$, và mômen điện từ thực $T_{m_real} = \frac{P_{out}}{\Omega_m}$.
  2. Lý thuyết tản nhiệt và dẫn nhiệt phi tuyến: Sử dụng mô hình nhiệt đồng nhất kết hợp định luật Fourier và Newton về làm mát bề mặt, tính toán độ tăng nhiệt bão hòa $\theta_\infty = \frac{Q}{\alpha S_c}$ và hằng số thời gian phát nhiệt $T = \frac{C}{\alpha S_c}$, có xét đến biến thiên hệ số tản nhiệt bề mặt $\alpha = \alpha_0 (1 + k_0 \sqrt{v})$.
  3. Lý thuyết tối ưu hóa toán học phi tuyến có ràng buộc: Ứng dụng điều kiện dừng chặt chẽ với các sai số hội tụ hàm mục tiêu, biến số và hàm ràng buộc đạt mức $10^{-5}$.

Điều kiện biên (Boundary conditions) được xác định rõ ràng: Cấp cách điện F với giới hạn tăng nhiệt độ cuộn dây $\Delta \theta \le 105^\circ\text{C}$ (nhiệt độ môi trường chuẩn $40^\circ\text{C}$); giới hạn cảm ứng từ khe hở không khí $B_g \in [0.65\text{T}, 0.85\text{T}]$; mật độ dòng điện stator $J \le 9\text{ A/mm}^2$ ở chế độ làm mát tự nhiên.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan khoa học thực chứng (Positivism) với lập trường tri thức luận khách quan chặt chẽ. Thiết kế nghiên cứu được tổ chức theo quy trình 4 giai đoạn khép kín:

[Giai đoạn 1: Mô hình toán giải tích & Thuật toán SQP (MATLAB)]
[Giai đoạn 2: Mô phỏng số phần tử hữu hạn 2D/3D (FEA Electromagnetic & Thermal)]
[Giai đoạn 3: Thiết kế chế tạo mô hình vật lý nguyên mẫu (Physical Prototype)]
[Giai đoạn 4: Đo kiểm thực nghiệm trên hệ thống Dyno Test-bench & Đối chuẩn]

Quy trình nghiên cứu rigorous

Quy trình tối ưu hóa giải tích được lập trình tự động trên nền tảng MATLAB thông qua module fmincon với thuật toán active-set. Cấu hình giải thuật thiết lập như sau:

"options = optimset('Algorithm','active-set','Display','iter','TolX',1e-5,'TolCon',1e-5,'TolFun',1e-5,'MaxFunEvals',1000,'FunValCheck','on');" — (Trích xuất từ mã nguồn tối ưu hóa giải tích của Luận án, Phụ lục D, tr. 109)

Mô hình tính toán xác định chính xác tổn hao và trọng lượng từng phần tử vật liệu:

"Wtot=1.01(Wcus+Wt+Wc+Wri+Wcur+Wcue); %(Kg)- Total wt"* "Tm_real = Pout/(Nm2pi/60); func =Wtot;" — (Trích xuất từ mã nguồn thuật toán đánh giá hiệu năng của Luận án, Phụ lục C, tr. 107-108)

Sau khi xác định được tập nghiệm tối ưu trên đường cong Pareto, cấu trúc hình học được chuyển sang phần mềm mô phỏng phần tử hữu hạn FEA để phân tích trường từ 2D và 3D. Bước cực $\tau$, hệ số Carter $k_g$, và góc nghiêng rãnh rotor $\xi$ được phân tích chi tiết nhằm kiểm tra hiện tượng bão hòa cục bộ tại đỉnh răng stator ($B_{st}$) và răng rotor ($B_{rt}$).

Trong mô phỏng nhiệt, hai trạng thái vận hành chính được khảo sát:

  1. Chế độ làm việc liên tục $S_1$: Động cơ vận hành kéo dài đến khi nhiệt độ toàn bộ các phần tử đạt trạng thái cân bằng nhiệt bão hòa $\theta_\infty$.
  2. Chế độ làm việc ngắn hạn $S_2$: Động cơ vận hành phát mômen cực đại trong khoảng thời gian ngắn định trước ($t_p = 10 - 30\text{ phút}$), sau đó dừng nghỉ hoàn toàn để nhiệt độ hạ về bằng nhiệt độ môi trường.

Data và phân tích

Mẫu thử nghiệm thực tế được chế tạo và lắp đặt trên hệ thống băng thử công suất chuyên dụng (test-bench) bao gồm: động cơ thử nghiệm, cảm biến mômen xoắn - tốc độ quay độ chính xác cao, máy phát tải phụ tải điều khiển được, hệ thống biến tần công nghiệp và khối thu thập dữ liệu nhiệt độ đa kênh gắn trực tiếp vào đầu cuộn dây, lõi sắt stator và thân vỏ.

Bảng dữ liệu thực nghiệm đối chuẩn giữa thiết kế tối ưu, mô phỏng phần tử hữu hạn và đo kiểm thực tế được thể hiện chi tiết tại các chế độ vận hành:

Điểm làm việc / Chế độ Điện áp $U$ (V) Tần số $f$ (Hz) Tốc độ $n$ (rpm) Mômen tính toán / FEA (N.m) Mômen đo thực nghiệm (N.m) Nhiệt độ cuộn dây mô phỏng ($^\circ\text{C}$) Nhiệt độ cuộn dây thực nghiệm ($^\circ\text{C}$)
Điểm 1 (Ngắn hạn $S_2$) $340$ $35.0$ $500$ $15.42$ $15.10$ $88.5$ $90.2$
Điểm 2 (Ngắn hạn $S_2$) $360$ $32.0$ $500$ $16.80$ $16.35$ $94.2$ $96.8$
Điểm 3 (Liên tục $S_1$) $350$ $42.0$ $700$ $8.65$ $8.50$ $72.4$ $74.1$
Điểm 4 (Liên tục $S_1$) $380$ $46.0$ $700$ $10.20$ $9.95$ $78.6$ $81.3$
Điểm 5 (Định mức $S_1$) $380$ $54.4$ $1000$ $7.15$ $7.02$ $68.0$ $70.5$

Sự hội tụ của thuật toán SQP diễn ra nhanh chóng sau 23-35 bước lặp (iterations), vượt trội hoàn toàn so với hàng nghìn thế hệ của giải thuật di truyền GA. Sai số tương đối giữa mô hình tính toán giải tích tối ưu, mô phỏng FEA và số liệu đo thực nghiệm kiểm chứng đều nằm trong biên độ hẹp ($< 3.2%$ đối với mômen và $< 4.5%$ đối với nhiệt độ bão hòa), khẳng định độ tin cậy và giá trị khoa học vững chắc của phương pháp.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Xác lập tập nghiệm tối ưu Pareto chính xác bằng thuật toán SQP kết hợp $\varepsilon$-constraint: Khác với các tiếp cận giải thuật ngẫu nhiên, phương pháp SQP xác định chính xác đường biên Pareto thể hiện sự đánh đổi giữa khối lượng điện từ $W_{tot}$ và hiệu suất $\eta$. Khi giảm khối lượng động cơ từ $14.5\text{ kg}$ xuống $10.2\text{ kg}$, hiệu suất động cơ giảm tương ứng từ $86.5%$ xuống $81.2%$, cung cấp cơ sở dữ liệu định lượng cho kỹ sư lựa chọn thiết kế tối ưu tùy theo yêu cầu cụ thể của từng ứng dụng truyền động.
  2. Hiệu ứng triệt tiêu sóng hài và giảm độ nhấp nhô mômen của rotor rãnh nghiêng: Kết quả phân tích phần tử hữu hạn 3D và đo kiểm thực nghiệm chỉ ra rằng, việc nghiêng rãnh rotor đi một bước răng stator đã làm giảm biên độ nhấp nhô mômen từ $18.4%$ (đối với rotor rãnh thẳng) xuống chỉ còn $4.2%$ (đối với rotor rãnh nghiêng), đồng thời loại bỏ gần như hoàn toàn sóng hài răng bậc cao trong sức điện động cảm ứng $e_{ssst}$.
  3. Cơ chế tản nhiệt và vai trò của cấu trúc nắp hở ở chế độ liên tục: Nghiên cứu phát hiện rằng ở chế độ làm việc liên tục $S_1$ với phương thức làm mát tự nhiên, thiết kế vỏ kín khiến nhiệt độ bão hòa cuộn dây vượt ngưỡng an toàn ($> 125^\circ\text{C}$). Việc chuyển đổi sang thiết kế nắp hở có lưu thông không khí cục bộ làm tăng hệ số tản nhiệt đối lưu bề mặt $\alpha$ từ $14.5\text{ W/m}^2\text{K}$ lên $23.8\text{ W/m}^2\text{K}$, giúp hạ nhiệt độ cuộn dây bão hòa xuống dưới $85^\circ\text{C}$ ở tốc độ $1000\text{ rpm}$.
  4. Kiểm chứng khả năng chịu quá tải mômen ở chế độ làm việc ngắn hạn $S_2$: Dữ liệu thực nghiệm xác nhận tại điểm làm việc $n = 500\text{ rpm}, U = 360\text{V}, f = 32\text{ Hz}$, động cơ phát sinh mômen cực đại $16.35\text{ N.m}$ (vượt yêu cầu thiết kế $15\text{ N.m}$) trong thời gian $15\text{ phút}$ mà nhiệt độ dây quấn chỉ đạt $96.8^\circ\text{C}$, nằm hoàn toàn trong giới hạn cho phép của cấp cách điện F ($155^\circ\text{C}$).
                             SO SÁNH ĐẶC TÍNH KẾT CẤU

Implications đa chiều

  • Về mặt học thuật (Theoretical Advances): Nghiên cứu đã làm phong phú lý thuyết thiết kế máy điện hiện đại bằng việc mô hình hóa thành công tương tác đa vật lý (điện từ - nhiệt) trong môi trường tối ưu hóa phi tuyến có ràng buộc, cung cấp phương pháp luận chuẩn xác cho các bài toán tối ưu máy điện xoay chiều.
  • Về mặt phương pháp luận (Methodological Innovations): Quy trình tích hợp từ giải tích MATLAB $\rightarrow$ FEA 2D/3D $\rightarrow$ Mô phỏng nhiệt $\rightarrow$ Thực nghiệm Test-bench thiết lập một tiêu chuẩn mới cho nghiên cứu chế tạo máy điện đặc chủng, có thể chuyển giao áp dụng cho động cơ đồng bộ từ trở (SynRM) hoặc động cơ bước công suất lớn.
  • Về ứng dụng thực tiễn (Practical Applications): Cung cấp tài liệu thiết kế hoàn chỉnh, thông số hình học dập lá thép silic và sơ đồ quấn dây cho các doanh nghiệp sản xuất động cơ điện trong nước, giảm phụ thuộc vào động cơ servo nhập khẩu đắt tiền từ Đức hay Nhật Bản.
  • Về chính sách và tự chủ công nghệ (Policy Recommendations): Đóng góp trực tiếp vào chiến lược phát triển sản phẩm công nghiệp hỗ trợ công nghệ cao của quốc gia, hiện thực hóa mục tiêu tự chủ thiết bị cơ điện tử trong dây chuyền sản xuất thông minh.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, nghiên cứu thẳng thắn thừa nhận 4 hạn chế mang tính điều kiện biên:

  1. Mô hình hóa điện từ 2D giai đoạn đầu chưa bao quát toàn diện hiệu ứng tản đầu nối 3D: Việc tính toán điện kháng tản đầu nối cuộn dây stator ($X_0$) và vành ngắn mạch rotor ($R_e, X_e$) trong mô hình giải tích tối ưu ban đầu vẫn dựa trên các công thức kinh nghiệm giải tích, đòi hỏi phải hiệu chỉnh bằng mô phỏng FEA 3D chi tiết.
  2. Giả thiết về tính đẳng hướng của vật liệu cách điện trong mô hình nhiệt tập trung: Mô hình nhiệt giả định hệ số dẫn nhiệt của lớp cách điện rãnh ($\lambda_c$) là đồng nhất theo mọi phương, chưa phản ánh hết sự bất đối xứng dẫn nhiệt tế vi do các bọt khí sót lại trong quá trình tẩm vec-ni cách điện.
  3. Chưa xét đến tổn hao sóng hài do điều biến độ rộng xung (PWM) của biến tần: Mô hình tính toán và mô phỏng hiện tại sử dụng nguồn cấp hình sin thuần túy, chưa tích hợp đầy đủ tổn hao sắt phụ do sóng hài bậc cao sinh ra từ tần số đóng cắt của khóa bán dẫn IGBT/SiC.
  4. Phạm vi kiểm chứng thực nghiệm tập trung vào dải công suất vừa và nhỏ: Mẫu thử nghiệm được chế tạo tối ưu cho mức mômen $15\text{ N.m}$ ở dải tốc độ dưới $3500\text{ rpm}$, chưa mở rộng khảo sát các động cơ servo công nghiệp siêu tốc ($> 10.000\text{ rpm}$).

Chương trình nghiên cứu tiếp theo (Future Research Agenda) mở ra 5 định hướng cụ thể:

  1. Xây dựng mô hình tối ưu hóa đồng thời liên kết điện từ - nhiệt - điều khiển (Co-simulation optimization) xét đến thuật toán điều khiển FOC và tần số đóng cắt PWM.
  2. Ứng dụng vật liệu từ tính tiên tiến như thép vô định hình (Amorphous) hoặc tôn silic tổn hao siêu thấp xử lý bằng công nghệ cắt laser không ứng suất dư.
  3. Nghiên cứu tích hợp hệ thống làm mát bằng chất lỏng trực tiếp bên trong vỏ động cơ cho các ứng dụng servo môi trường kín yêu cầu chuẩn IP67/IP68.
  4. Mở rộng bài toán tối ưu hóa đa mục tiêu sang 3-4 hàm mục tiêu đồng thời: Khối lượng, Hiệu suất, Mômen quán tính rotor ($J$), và Chi phí vật liệu tổng thể.
  5. Nghiên cứu công nghệ tự động hóa khâu lồng dây stator và đúc áp lực rotor nhôm/đồng độ dẫn điện cao trên dây chuyền công nghiệp quy mô lớn.

Tác động và ảnh hưởng

  • Tác động học thuật (Academic Impact): Đặt nền móng học thuật vững chắc cho chuyên ngành Kỹ thuật điện tại Việt Nam trong lĩnh vực thiết kế động cơ servo; mở ra hàng loạt đề tài nghiên cứu chuyên sâu cho học viên cao học và nghiên cứu sinh tiến sĩ. Dự kiến công trình sẽ tạo ra nguồn trích dẫn học thuật quan trọng trong các nghiên cứu về máy điện đặc chủng và điều khiển chuyển động chính xác.
  • Chuyển đổi công nghiệp (Industry Transformation): Tạo bước đột phá kỹ thuật giúp các đơn vị sản xuất máy điện trong nước (như HEM, VIHEM, ABB Việt Nam) nâng cao năng lực công nghệ từ chế tạo động cơ công nghiệp tiêu chuẩn sang làm chủ hoàn toàn quy trình thiết kế, chế tạo động cơ servo chất lượng cao phục vụ các nhà máy chế tạo máy CNC và hệ thống cánh tay robot.
  • Tác động kinh tế - xã hội (Societal Benefits): Giảm chi phí đầu tư hệ thống tự động hóa cho các doanh nghiệp cơ khí chế tạo vừa và nhỏ từ $30 - 45%$ so với nhập khẩu nguyên chiếc từ các thương hiệu phương Tây hay Nhật Bản; nâng cao hàm lượng khoa học công nghệ nội địa trong các sản phẩm công nghiệp mang thương hiệu Việt Nam.

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

                              CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công mô hình toán học giải tích phi tuyến đa tham số tích hợp đồng thời trường điện từ và trường nhiệt, sau đó giải quyết bài toán bằng thuật toán SQP kết hợp phương pháp ràng buộc $\varepsilon$-constraint. Công trình đã mở rộng lý thuyết tính toán thiết kế máy điện kinh điển của Alger và Kostenko sang miền tối ưu hóa đa mục tiêu thời gian thực, cho phép xác lập chính xác phân bố nghiệm tối ưu Pareto giữa khối lượng điện từ $W_{tot}$ và tổn hao năng lượng $1-\eta$, giải quyết triệt để sự mất cân đối giữa kích thước cơ khí và độ tăng nhiệt trong động cơ servo.

2. Điểm mới về phương pháp luận nghiên cứu khi so sánh với các công trình quốc tế tiêu biểu?

So với nghiên cứu của Svante Andersson (2000) vốn chỉ giải bài toán đơn mục tiêu trên động cơ đồng bộ PMSM, và nghiên cứu của Mehmet Çunkaş & Ramazan Akkaya (2006) sử dụng thuật toán di truyền GA cho động cơ không đồng bộ với thời gian tính toán lớn và nghiệm ngẫu nhiên, luận án đã:

  • Ứng dụng thành công thuật toán xác định SQP đạt tốc độ hội tụ bậc hai siêu nhanh ($< 35$ bước lặp) với sai số dừng chặt chẽ $10^{-5}$.
  • Xây dựng quy trình xác thực 4 cấp khép kín: Mô hình giải tích $\rightarrow$ Tối ưu hóa MATLAB $\rightarrow$ Mô phỏng phần tử hữu hạn FEA 2D/3D & Mạng nhiệt $\rightarrow$ Băng thử Test-bench thực tế, đảm bảo tính tin cậy tuyệt đối của kết quả nghiên cứu.

3. Phát hiện bất ngờ nhất trong quá trình nghiên cứu thực nghiệm là gì?

Phát hiện bất ngờ nhất là độ tăng nhiệt của động cơ ở chế độ làm việc ngắn hạn $S_2$ với tải mômen cực đại $16.35\text{ N.m}$ ở tốc độ thấp ($500\text{ rpm}$) diễn ra chậm hơn đáng kể so với dự đoán của mô hình nhiệt tuyến tính cổ điển. Nhờ nhiệt dung riêng của khối lõi sắt stator và rotor hấp thụ phần lớn nhiệt lượng ban đầu ($Q$), nhiệt độ cuộn dây trong $15\text{ phút}$ đầu chỉ đạt $96.8^\circ\text{C}$, thấp hơn nhiều so với nhiệt độ bão hòa liên tục ($> 135^\circ\text{C}$). Điều này mở ra khả năng khai thác mômen vượt định mức lên tới $160 - 200%$ trong các chu kỳ làm việc ngắn mà không cần gia tăng kích thước khung vỏ động cơ.

4. Luận án có cung cấp đầy đủ quy trình và thông số để tái lập nghiên cứu (Replication Protocol) không?

Luận án cung cấp chi tiết toàn bộ mã nguồn thuật toán tối ưu hóa trong môi trường MATLAB (Phụ lục C và D), đầy đủ hệ phương trình giải tích, thông số hình học chi tiết của stator/rotor (đường kính ngoài, đường kính trong, kích thước rãnh, bề rộng răng, khe hở không khí), đặc tính vật liệu thép silic và dây quấn đồng, cùng cấu hình cài đặt thử nghiệm trên băng thử dyno (Phụ lục E). Bất kỳ nhóm nghiên cứu nào cũng có thể tái lập hoàn toàn quy trình tính toán, mô phỏng và chế tạo thực tế dựa trên dữ liệu công bố.

5. Định hướng lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình 10 năm được định hướng qua 3 giai đoạn chiến lược:

  • Giai đoạn 1 (1–3 năm): Hoàn thiện mô hình tối ưu hóa co-simulation đa vật lý liên kết trực tiếp với thuật toán điều khiển vector không gian (SVPWM-FOC), xét đến tổn hao sóng hài biến tần.
  • Giai đoạn 2 (4–6 năm): Nghiên cứu chế tạo động cơ servo không đồng bộ công suất lớn ($15 - 55\text{ kW}$) tích hợp áo làm mát bằng dung dịch phục vụ xe điện và máy công cụ hạng nặng.
  • Giai đoạn 3 (7–10 năm): Ứng dụng trí tuệ nhân tạo (AI-driven surrogate modeling) và công nghệ in 3D kim loại để tự động hóa hoàn toàn quy trình thiết kế và chế tạo các thế hệ động cơ servo thông minh tự chẩn đoán nhiệt - cơ.

Kết luận

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Đức Bắc là một công trình khoa học công phu, nghiêm túc và có giá trị học thuật cũng như giá trị thực tiễn rất cao trong lĩnh vực Kỹ thuật điện và Tự động hóa truyền động. Toàn bộ các mục tiêu nghiên cứu đề ra ban đầu đã được hoàn thành trọn vẹn với các đóng góp cốt lõi được đúc kết như sau:

  1. Xây dựng thành công hệ thống mô hình toán học giải tích đa vật lý tích hợp trọn vẹn quan hệ phi tuyến điện từ, mạch từ bão hòa, tổn hao công suất và tản nhiệt đối lưu/dẫn nhiệt cho động cơ servo không đồng bộ 3 pha rotor lồng sóc.
  2. Phát triển và hiện thực hóa thuật toán tối ưu hóa đa mục tiêu xác định SQP kết hợp phương pháp ràng buộc $\varepsilon$-constraint, giải quyết thành công bài toán tìm tập nghiệm Pareto tối ưu giữa khối lượng vật liệu điện từ ($W_{tot}$) và hiệu suất vận hành ($\eta$) với tốc độ hội tụ vượt trội.
  3. Chứng minh và định lượng hóa bằng mô phỏng phần tử hữu hạn 3D vai trò quyết định của cấu trúc rotor rãnh nghiêng trong việc triệt tiêu sóng hài sức điện động cảm ứng, hạ độ nhấp nhô mômen xoắn xuống mức $4.2%$, đảm bảo tính năng vận hành êm ái cho hệ truyền động servo.
  4. Xác lập mô hình trường nhiệt chính xác ở các chế độ làm việc $S_1$ và $S_2$, làm rõ cơ chế phát nhiệt và đề xuất giải pháp kết cấu nắp hở giúp hạ nhiệt độ cuộn dây bão hòa xuống dưới $85^\circ\text{C}$, đảm bảo độ bền nhiệt vượt cấp cách điện F.
  5. Thiết kế, chế tạo nguyên mẫu vật lý hoàn chỉnh và thử nghiệm thực chứng thành công trên hệ thống test-bench, khẳng định sự trùng khớp tuyệt vời giữa lý thuyết tính toán, mô phỏng số và dữ liệu vận hành thực tế với sai số dưới $5%$.
  6. Mở ra một phương pháp luận thiết kế chuẩn mực, tạo tiền đề công nghệ vững chắc thúc đẩy ngành sản xuất máy điện và thiết bị tự động hóa trong nước tiến lên nấc thang tự chủ công nghệ cao trong kỷ nguyên chuyển đổi số công nghiệp.