Tổng quan về luận án

Trong kỹ thuật cơ khí, xây dựng công trình và động lực học máy hiện đại, các kết cấu dạng thanh và trục chịu tải trọng động phức tạp thường xuất hiện các dao động cưỡng bức nguy hiểm. Hiện tượng phá hủy do cộng hưởng xảy ra khi tần số kích động từ môi trường hoặc nguồn dẫn động trùng khớp với tần số dao động riêng của hệ thống. Nhằm giải quyết triệt để bài toán này, luận án tiến sĩ chuyên ngành Cơ học Kỹ thuật / Kỹ thuật Cơ khí với đề tài "Ứng dụng lý thuyết điều khiển trong tối ưu tần số riêng và khối lượng của kết cấu thanh" đã tiên phong ứng dụng Nguyên lý Cực đại Pontryagin (Pontryagin's Maximum Principle - PMP) để tái cấu trúc hình học, kiểm soát dải tần số riêng và khối lượng của hệ kết cấu đàn hồi 1D.

Khoảng trống nghiên cứu (research gap) cốt lõi được luận án chỉ ra là: các nghiên cứu trước đây chủ yếu dừng lại ở bài toán tối ưu đơn mục tiêu (cực đại tần số riêng thứ nhất hoặc cực tiểu hóa thể tích riêng lẻ), hoặc các bài toán đa mục tiêu đơn giản dạng $Z = (1 - k_1)\frac{\omega_1}{\omega_{01}} + k_1\frac{W}{W_0} \to \min$ (Trần Đức Trung & Bùi Hải Lê, 2009; Minh-Quy Le et al., 2012). Cấu trúc hàm mục tiêu kinh điển này hoàn toàn bỏ sót các kịch bản quan trọng khác trong kỹ thuật như đồng thời cực tiểu hóa tần số và cực đại hóa khối lượng, hoặc mở rộng dải tần giữa các bậc dao động bậc cao ($\omega_1, \omega_2, \dots, \omega_n$). Luận án giải quyết hai câu hỏi nghiên cứu trung tâm:

  1. RQ1: Làm thế nào để thiết lập một hàm đa mục tiêu tổng quát và xác định tường minh hệ số tỷ lệ liên hợp $k$ trong hàm Hamilton nhằm điều khiển đồng thời nhiều tần số riêng và khối lượng kết cấu?
  2. RQ2: Quy luật phân bố hình học tối ưu tương đương giữa các kết cấu có chiều dài và điều kiện biên khác nhau được hình thành như thế nào dựa trên dạng dao động riêng?

Nghiên cứu đặt ra hai giả thuyết khoa học:

  • H1: Tồn tại một hệ số tỷ lệ giải tích $k$ liên kết trực tiếp biến trạng thái dao động với biến trạng thái liên hợp trong hàm Hamilton cho bài toán đa mục tiêu tổng quát, giúp thuật toán PMP hội tụ xác định mà không phụ thuộc vào phép dò tìm số ngẫu nhiên.
  • H2: Cấu hình phân bố vật liệu tối ưu của kết cấu thanh chịu dao động dọc và trục chịu xoắn có mối tương quan định tính trực tiếp với các điểm nút và điểm bụng của dạng dao động riêng tương ứng.

Khung lý thuyết của nghiên cứu tích hợp Lý thuyết Điều khiển Tối ưu Pontryagin (1956), Phương pháp Ma trận Truyền (Transfer Matrix Method - TMM), và Phương trình Động lực học Thanh - Trục đàn hồi Saint-Venant / Euler-Bernoulli. Phạm vi khảo sát tập trung vào các kết cấu thanh chịu dao động dọc, trục tròn chịu dao động xoắn, và dầm chịu uốn với vật liệu đẳng hướng, đàn hồi tuyến tính dưới 4 cấu hình điều kiện biên chuẩn hóa: Ngàm - Tự do, Tự do - Tự do, Ngàm - Ngàm, và Tự do - Ngàm. Kết quả mô phỏng số thực hiện trên các mô hình rời rạc hóa từ 6 bậc, 21 nút, 41 nút đến 81 nút, chứng minh khả năng dịch chuyển dải tần số riêng lên tới 15% - 35% trong khi tiết giảm khối lượng tổng thể từ 8% - 20%.

Literature Review và Positioning

Lịch sử ứng dụng lý thuyết điều khiển tối ưu trong cơ học kết cấu ghi nhận bước ngoặt khi Viện sĩ Lev Semionovich Pontryagin công bố Nguyên lý Cực đại vào năm 1956 (Pontryagin, 1956, 1958). Sau đó, Grinev và Filipov (1975), Grinev (1975) đã đặt nền móng đưa PMP vào cơ học vật rắn biến dạng để giải quyết bài toán tối ưu lực tới hạn và tìm hình dáng thanh có khối lượng cực tiểu với tần số cho trước. Tại Việt Nam, trường phái nghiên cứu của Trần Đức Trung và Vũ Thế Bình (1977), Trần Đức Trung và Nguyễn Dũng (1978, 1979) đã phát triển việc ứng dụng PMP để tối ưu hóa thanh chịu nén đúng tâm, tối ưu dao động dọc của thanh ngàm - tự do có tính đến khối lượng tập trung.

Giai đoạn thập niên 1980 - 2010 chứng kiến sự phát triển mạnh mẽ của các nghiên cứu quốc tế:

  • Szymczak (1983, 1984) áp dụng PMP thiết kế tối ưu dầm thành mỏng tiết diện chữ I chịu xoắn và uốn có xét lực dọc trục.
  • Atanackovic và Simic (1999) sử dụng PMP khảo sát sự phân nhánh và cấu hình tối ưu của cột Pflüger chống mất ổn định đàn hồi.
  • Glavardanov và Atanackovic (2001), Atanackovic (2001) tối ưu hóa trục quay đàn hồi chịu tải tập trung để chống mất ổn định với thể tích cực tiểu.
  • Atanackovic và Braun (2005), Braun (2008) kết hợp công thức Rayleigh và PMP giải bài toán biên phi tuyến bậc 4 và bậc 5 cho thanh quay tròn, chứng minh tiết diện đầu tự do có giá trị hữu hạn, khác biệt với nghiệm suy biến bằng 0 của lý thuyết Bernoulli - Euler cổ điển.
  • Atanackovic và cộng sự (2006, 2010, 2012) mở rộng PMP cho dầm trên nền đàn hồi Winkler và cấu trúc micro/nano dựa trên lý thuyết đàn hồi phi cục bộ của Eringen.

Trong bức tranh tổng quan đó, tồn tại hai luồng quan điểm đối nghịch trong tối ưu hóa kết cấu: (1) Nhóm phương pháp tối ưu toán học thuần túy (Mathematical Programming / Quy hoạch phi tuyến) tìm kiếm trong Không gian Tham số (Parametric Space - PS) đòi hỏi khối lượng tính toán khổng lồ và dễ rơi vào cực trị địa phương; (2) Nhóm phương pháp biến phân và nguyên lý cực trị (PMP, Bellman) giải bài toán trong Không gian Trạng thái (State Space - SS). Luận án định vị nghiên cứu của mình vào nhóm thứ hai nhưng khắc phục triệt để nhược điểm của các công trình trước đó (Trần Đức Trung & Bùi Hải Lê, 2009, 2010; Minh-Quy Le et al., 2012).

Bằng cách thiết lập phiếm hàm mục tiêu Maier mở rộng và trích xuất tường minh hệ số tỷ lệ $k$, luận án định vị như một công trình bước ngoặt, mở rộng từ bài toán tối ưu trị riêng đơn lẻ sang bài toán tối ưu đồng thời dải tần số đa bậc $(\omega_1, \omega_2)$ kết hợp kiểm soát khối lượng $W$, xây dựng hoàn chỉnh Tập giải pháp khả thi (Feasible Region) và Mặt trận Pareto (Pareto Front).

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

graph TD
    A["Hệ phương trình vi phân trạng thái dao động"] --> B["Hệ phương trình vi phân liên hợp & Hàm Hamilton H"]
    C["Phiếm hàm mục tiêu Maier tổng quát (Tần số & Khối lượng)"] --> B
    B --> D["Xác định tường minh hệ số tỷ lệ liên hợp k"]
    D --> E["Thiết lập Điều kiện cần tối ưu (NOC) theo PMP"]
    E --> F["Phương pháp Ma trận truyền TMM & Thuật toán lặp số"]
    F --> G["Tập Pareto đa mục tiêu & Tập giải pháp khả thi"]
    F --> H["Quy luật Cấu hình tối ưu tương đương & Dự báo định tính"]

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

Luận án tạo ra bước đột phá khi mở rộng và hoàn thiện lý thuyết điều khiển tối ưu cấu trúc kết cấu thông qua các đóng góp lý thuyết cụ thể:

  1. Xây dựng hệ hàm đa mục tiêu tổng quát: Khắc phục hạn chế của các hàm mục tiêu cục bộ trước đây, luận án chuẩn hóa không gian tối ưu cho 4 tổ hợp trạng thái mục tiêu: $$\begin{cases} \text{TH1: } \omega_i \to \max, W \to \min \ \text{TH2: } \omega_i \to \min, W \to \max \ \text{TH3: } \omega_i \to \max, W \to \max \ \text{TH4: } \omega_i \to \min, W \to \min \end{cases}$$ Đồng thời mở rộng cho bài toán tối ưu hóa khoảng cách giữa hai tần số riêng liên tiếp: $F = k_1 \frac{\omega_1}{\omega_{01}} - k_2 \frac{\omega_2}{\omega_{02}} \to \max(\min)$.
  2. Xác định tường minh hệ số tỷ lệ liên hợp $k$: Luận án chứng minh mối quan hệ tỷ lệ giữa biến trạng thái dao động ($u, \phi, N, M$) và biến trạng thái liên hợp ($p_u, p_\phi, p_N, p_M$) trong hàm Hamilton $H$. Nhờ phát hiện này, việc giải hệ phương trình vi phân liên hợp phức tạp được quy về việc giải hệ phương trình trạng thái dao động tự do với hệ số tỷ lệ xác định giải tích, giảm bậc tự do tính toán xuống 50%.
  3. Mô hình hóa biến mở rộng trong Không gian Trạng thái: Tần số riêng $\omega$ và khối lượng $W$ được chuyển đổi thành các biến trạng thái mở rộng tuân theo phương trình $\frac{d\omega}{dx} = 0$ và $\frac{dW}{dx} = f_W(U, x)$, biến bài toán biên thành bài toán điều khiển trạng thái cuối dạng Maier: $$G = \sum_{i=1}^m b_i q_i(x_f)$$

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

Khung phân tích của luận án là sự hợp nhất của ba trụ cột lý thuyết: (1) Nguyên lý Cực đại Pontryagin, (2) Phương pháp Ma trận Truyền (Transfer Matrix Method), và (3) Lý thuyết Đa mục tiêu Pareto.

Hàm Hamilton tổng quát được thiết lập dưới dạng: $$H(p, q, U, x) = \sum_{i=1}^n p_i f_i(q, U, x)$$ Biến điều khiển (biến thiết kế) $U(x)$ là diện tích mặt cắt ngang $A(x)$ hoặc đường kính trục $d(x)$, bị chặn trong miền đóng $U(x) \in [U_{\min}, U_{\max}]$. Điều kiện cần tối ưu (Optimal Necessary Condition - NOC) đòi hỏi hàm Hamilton đạt cực trị đối với mọi $U \in [U_{\min}, U_{\max}]$: $$H(p^, q^, U^, x) = \max_{U \in G_u} H(p^, q^*, U, x)$$

Điều kiện biên áp dụng cho hệ liên hợp được xác định qua phương trình hoành: $$\sum_{i=1}^n p_i(x_f) \delta q_i(x_f) - \sum_{i=1}^n p_i(x_0) \delta q_i(x_0) + \delta F = 0$$

Biên điều kiện xác định tính hợp lệ của mô hình gồm: vật liệu làm việc trong miền đàn hồi tuyến tính, biến dạng nhỏ, mặt cắt ngang biến thiên bậc (stepped bar/shaft) với chiều dài mỗi đoạn $L_e$ đủ lớn để ứng suất trượt cục bộ Saint-Venant triệt tiêu nhanh chóng.

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

graph LR
    subgraph "Mô hình hóa Toán học"
        M1["Hệ PTVP trạng thái"] --> M2["Ma trận truyền TMM"]
        M2 --> M3["Hàm Hamilton & Biến liên hợp"]
    end
    subgraph "Giải thuật Số (MATLAB)"
        A1["Khởi tạo U_0"] --> A2["Giải TMM tìm omega_i, q_i"]
        A2 --> A3["Tính p_i qua hệ số k"]
        A3 --> A4["Cực trị hóa hàm Hamilton H"]
        A4 --> A5{"Hội tụ |Delta F| < epsilon?"}
        A5 -- No --> A1
        A5 -- Yes --> A6["Xuất tập Pareto & Cấu hình"]
    end
    subgraph "Kiểm chứng Độc lập (ANSYS)"
        V1["Mô hình 3D Solid"] --> V2["Phân tích Modal & Harmonic"]
        V2 --> V3["Đối chuẩn tần số & Dạng dao động"]
    end
    M3 --> A1
    A6 --> V1

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng định lượng (Positivism), dựa trên nền tảng toán học giải tích kết hợp mô phỏng số chính xác cao. Thiết kế nghiên cứu đa tầng (Multi-level design) được cấu trúc qua 4 cấp độ:

  • Cấp độ 1: Thiết lập phương trình vi phân trạng thái liên tục từ định luật nhiệt động lực học và nguyên lý Ostrogradsky - Hamilton.
  • Cấp độ 2: Rời rạc hóa cấu trúc thành $n$ phân tố bậc hữu hạn kết hợp ma trận truyền giải tích cục bộ.
  • Cấp độ 3: Tối ưu hóa điều khiển lặp theo PMP bằng chương trình tự phát triển trên môi trường MATLAB.
  • Cấp độ 4: Kiểm nghiệm chéo và đối chuẩn (Benchmarking) bằng phương pháp phần tử hữu hạn phi tuyến trên phần mềm công nghiệp ANSYS.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu số tuân theo giao thức kiểm soát sai số nghiêm ngặt:

  1. Rời rạc hóa ma trận truyền: Với mỗi phân đoạn $e$ có chiều dài $L_e$, ma trận truyền $[H_e]$ cho dao động xoắn có dạng: $$[H_e] = \begin{bmatrix} \cos\left(\frac{\omega}{a} L_e\right) & \frac{a}{G J_p} \sin\left(\frac{\omega}{a} L_e\right) \ -\frac{G J_p}{a} \sin\left(\frac{\omega}{a} L_e\right) & \cos\left(\frac{\omega}{a} L_e\right) \end{bmatrix}$$ với $a = \sqrt{G/\rho}$. Tương tự cho dao động dọc với $a = \sqrt{E/\rho}$ và độ cứng $EA$.
  2. Khử sai số đầu cuối: Thiết lập ma trận truyền toàn hệ $[H] = \prod_{e=1}^{n-1} [H_e]$. Điều kiện biên ở hai đầu (Ngàm - Tự do, Tự do - Tự do, Ngàm - Ngàm) dẫn đến phương trình đại số xác định tần số riêng $\det([H_{BC}]) = 0$.
  3. Thuật toán lặp cải tiến theo hàm Hamilton: Khởi tạo phân bố điều khiển $U_0(x) \in [U_{\min}, U_{\max}]$, giải bài toán trị riêng tìm $\omega_i^{(k)}$ và dạng riêng $q_i^{(k)}$, xác định các biến liên hợp $p_i^{(k)}$ thông qua hệ số $k$, cập nhật $U^{(k+1)}(x)$ từ điều kiện cực trị của $H$. Quá trình dừng khi $|\omega_i^{(k+1)} - \omega_i^{(k)}| / \omega_i^{(k)} < 10^{-5}$ và $|W^{(k+1)} - W^{(k)}| / W^{(k)} < 10^{-5}$.

Data và phân tích

Dữ liệu mô phỏng được thực hiện trên tập mẫu tham số gồm các trục và thanh có kích thước kỹ thuật tiêu chuẩn:

  • Mô hình kiểm tra độ tin cậy: Trục 6 đoạn, 7 nút ($n=7$), mở rộng kiểm tra tính độc lập của lưới với các mô hình 21 nút, 41 nút, và 81 nút.
  • Thông số vật liệu: Thép kết cấu chuẩn có $E = 2.1 \times 10^{11} \text{ N/m}^2$, $G = 8.1 \times 10^{10} \text{ N/m}^2$, khối lượng riêng $\rho = 7850 \text{ kg/m}^3$.
  • Biến thiết kế bị chặn: Cận dưới $d_{\min} = 0.02 \text{ m}$, cận trên $d_{\max} = 0.08 \text{ m}$.
  • Kiểm tra tính vững chắc (Robustness Checks): So sánh nghiệm tối ưu PMP với kết quả phân tích phổ (Spectral Method), phần tử hữu hạn ANSYS (phần tử SOLID186), và phương pháp On-Off (Bang-bang control) cho bài toán dầm uốn. Sai số tương đối giữa mô phỏng PMP trên MATLAB và kiểm chứng ANSYS luôn duy trì dưới mức $0.85%$, chứng minh độ tin cậy tuyệt đối của giải thuật.

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

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

Trích xuất trực tiếp từ các dữ liệu và bảng biểu số học trong luận án, các phát hiện đột phá bao gồm:

"Khi tần số dao động cưỡng bức trùng hoặc gần với tần số dao động riêng sẽ xảy ra hiện tượng cộng hưởng làm dao động tăng lên rất mạnh. Do vậy kiểm soát và điều khiển được quá trình dao động, đặc biệt là hiện tượng cộng hưởng là yêu cầu rất cần thiết."

  1. Khả năng tái cấu trúc mở rộng dải tần vượt trội: Trong bài toán tối ưu đa mục tiêu đồng thời hai tần số riêng đầu tiên ($\omega_1, \omega_2$) của trục chịu xoắn ngàm - tự do (Bảng 4.1 và Bảng 4.2), thuật toán PMP cho phép tăng tần số riêng thứ nhất $\omega_1$ lên $28.6%$ trong khi dịch chuyển $\omega_2$ theo hướng mong muốn để triệt tiêu hoàn toàn dải cộng hưởng làm việc của động cơ.
  2. Quy luật cấu hình tối ưu định tính dựa trên dạng dao động riêng: Luận án phát hiện rằng cấu hình hình học tối ưu của thanh và trục phân bố hoàn toàn trùng khớp với phân bố năng lượng biến dạng của các dạng riêng dao động. Tại các vị trí nút dao động (biên độ dịch chuyển bằng 0, nội lực cực đại), tiết diện tự động đạt cận trên $d_{\max}$; ngược lại tại các vị trí bụng dao động (biên độ cực đại, nội lực bằng 0), tiết diện co lại về cận dưới $d_{\min}$. Phát hiện này cho phép kỹ sư phác thảo ngay cấu hình tối ưu định tính mà không cần chạy thuật toán số.
  3. Quy luật cấu hình tối ưu tương đương: Khảo sát trên các trục có chiều dài $L$ và điều kiện biên khác nhau (Ngàm - Tự do chiều dài $L$ tương đương với một nửa trục Ngàm - Ngàm chiều dài $2L$), luận án chứng minh sự tương đương toán học về quy luật biến thiên tiết diện không thứ nguyên dọc theo tọa độ chuẩn hóa $\xi = x/L$.
  4. Giải pháp điều khiển On-Off (Bang-bang) cho độ cứng dầm chịu uốn: Đối với dầm chịu uốn tĩnh và uốn động, điều khiển tối ưu theo PMP chuyển đổi trạng thái dạng nhảy bậc (switching function) giữa hai cận $U_{\min}$ và $U_{\max}$, chứng minh dạng dầm 2 bậc tối ưu có độ võng tại đầu tự do nhỏ hơn $34.2%$ so với dầm tiết diện không đổi có cùng khối lượng.
  5. Đánh đổi đa mục tiêu trên Tập Pareto: Luận án định lượng hóa mức độ trade-off giữa tần số và khối lượng thông qua các hệ số không thứ nguyên $R_1, R_2, R_{12}, R_W$, cung cấp một tập nghiệm khả thi liên tục cho phép người thiết kế chủ động lựa chọn phương án kết cấu tùy theo ưu tiên tải trọng hay thể tích.

Implications đa chiều

  • Về mặt học thuật: Mở rộng phạm vi ứng dụng của Nguyên lý Cực đại Pontryagin sang lĩnh vực động lực học phi tuyến và tối ưu hóa hệ đa bậc tự do, cung cấp công cụ toán học giải quyết bài toán biên phức tạp trong cơ học vật rắn.
  • Về mặt phương pháp luận: Sự kết hợp giữa Ma trận truyền (TMM) và PMP tạo ra một khung tính toán hiệu năng cao, giảm chi phí thời gian tính toán hàng chục lần so với các thuật toán tiến hóa (Genetic Algorithms - GA, Particle Swarm Optimization - PSO) khi tìm kiếm cấu hình liên tục.
  • Về mặt ứng dụng thực tiễn: Cung cấp giải pháp thiết kế trực tiếp cho trục truyền động turbine, trục động cơ máy đùn ép nhựa, máy lốc tôn, máy đột dập công nghiệp nhằm triệt tiêu dao động tự do và tăng tuổi thọ làm việc của ổ bi đỡ.
  • Về mặt thiết kế chính sách & tiêu chuẩn kỹ thuật: Cung cấp căn cứ khoa học để hoàn thiện các quy chuẩn thiết kế giảm chấn thụ động và tiêu chuẩn kiểm định an toàn động lực học cho trục máy quay tốc độ cao trong công nghiệp chế tạo máy.

Limitations và Future Research

Mặc dù đạt được những đột phá quan trọng, nghiên cứu ghi nhận một số giới hạn học thuật cần được mở rộng:

  1. Giới hạn miền tuyến tính: Mô hình nghiên cứu dừng lại ở dao động tự do tuyến tính, chưa xét đến các hiệu ứng phi tuyến hình học khi góc xoay lớn hoặc biến dạng dọc trục vượt quá giới hạn đàn hồi Hooke.
  2. Vật liệu đồng nhất đẳng hướng: Các phương trình trạng thái được thiết lập cho vật liệu kim loại đồng chất, chưa mở rộng cho vật liệu composite phân lớp Functionally Graded Materials (FGM) hoặc vật liệu thông minh có tính chất biến thiên liên tục.
  3. Bỏ qua lực cản môi trường: Khảo sát chưa tính đến ma sát trong của vật liệu (damping nội) và lực cản nhớt khí động học bên ngoài trong quá trình dao động tự do.
  4. Ràng buộc công nghệ gia công: Tiết diện tối ưu dạng bước nhảy bậc (stepped profile) tại các điểm chuyển tiếp xuất hiện sự tập trung ứng suất, cần được làm mịn (fillet rounding) trong sản xuất cơ khí thực tế.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda) được đề xuất:

  • Hướng 1: Mở rộng thuật toán PMP cho bài toán tối ưu động lực học kết cấu phi tuyến chịu kích động ngẫu nhiên (stochastic vibration).
  • Hướng 2: Tích hợp giảm chấn khối lượng điều chỉnh (Tuned Mass Damper - TMD) chủ động và bán chủ động kết hợp tối ưu hóa hình học dầm/trục.
  • Hướng 3: Ứng dụng PMP cho kết cấu dầm/tấm nano dựa trên lý thuyết đàn hồi phi cục bộ Eringen có xét đến hiệu ứng bề mặt.
  • Hướng 4: Phát triển thuật toán tự động làm mịn biên và kết nối trực tiếp với công nghệ in 3D kim loại (Additive Manufacturing).
  • Hướng 5: Nghiên cứu tối ưu hóa đa trường vật lý (Thermo-mechanical-vibration coupling) cho cánh tuabin làm việc trong môi trường nhiệt độ cao.

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

Luận án tạo ra dấu ấn học thuật sâu sắc với 08 công trình khoa học công bố trên các diễn đàn uy tín, bao gồm:

  • 02 bài báo quốc tế trên tạp chí ISI/SCI uy tín: Meccanica (Nhà xuất bản Springer, Q1/Q2 về Cơ học Kỹ thuật).
  • 01 bài báo trên Tạp chí Cơ học - Viện Hàn lâm Khoa học và Công nghệ Việt Nam.
  • 03 bài báo trên Tạp chí Khoa học và Công nghệ - Viện Hàn lâm Khoa học và Công nghệ Việt Nam.
  • 02 bài báo trong Kỷ yếu Hội nghị Khoa học Cơ học Kỹ thuật Toàn quốc và Quốc tế.

Về mặt công nghiệp và kinh tế xã hội, giải pháp tối ưu hóa giúp giảm từ 10% - 18% lượng vật liệu kim loại chế tạo trục trong các dây chuyền máy công cụ hạng nặng, đồng thời kéo dài chu kỳ bảo dưỡng ổ đỡ lên 1.5 - 2 lần nhờ triệt tiêu biên độ rung động cộng hưởng.

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

graph TD
    User["Các nhóm thụ hưởng giá trị nghiên cứu"]
    User --> D["Nghiên cứu sinh & Nhà khoa học trẻ"]
    User --> S["Giảng viên & Chuyên gia Cơ học"]
    User --> I["Kỹ sư R&D & Doanh nghiệp Cơ khí"]
    User --> P["Cơ quan Quản lý & Tổ chức Tiêu chuẩn"]

    D --- D1["Kế thừa khung toán học PMP & Ma trận truyền TMM"]
    S --- S1["Giáo trình giảng dạy Động lực học & Điều khiển nâng cao"]
    I --- I1["Quy trình thiết kế trục máy tối ưu giảm chấn"]
    P --- P1["Cập nhật tiêu chuẩn định mức rung động kỹ thuật"]
  • Nghiên cứu sinh và Giới nghiên cứu học thuật: Kế thừa khung lý thuyết giải tích số PMP kết hợp Ma trận truyền để mở rộng cho các bài toán cơ học kết cấu phức tạp hơn.
  • Giảng viên và Chuyên gia Động lực học: Tài liệu tham khảo chuẩn mực cho các học phần cao học chuyên sâu về Điều khiển Tối ưu Hệ cơ học, Tối ưu hóa Kết cấu.
  • Kỹ sư R&D và Nhà sản xuất Chế tạo máy: Ứng dụng quy luật phân bố tiết diện định tính và bộ công cụ MATLAB/ANSYS để thiết kế tối ưu hệ thống trục truyền động, cánh khuấy, rotor tuabin đạt tiêu chuẩn chống rung cao cấp.
  • Cơ quan Quản lý Tiêu chuẩn Kỹ thuật: Cơ sở dữ liệu khoa học tin cậy để xây dựng định mức kỹ thuật về giới hạn dải tần an toàn cho thiết bị công nghiệp.

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 độc đáo nhất là việc thiết lập Hàm đa mục tiêu tổng quát tích hợp đầy đủ 4 trường hợp cực trị tần số riêng và khối lượng ($\omega_i \to \max/\min, W \to \max/\min$), đồng thời xác định tường minh hệ số tỷ lệ liên hợp $k$ trong hàm Hamilton. Luận án đã mở rộng trực tiếp Lý thuyết Điều khiển Tối ưu Pontryagin (1956) và khắc phục khoảng trống nghiên cứu của Trần Đức Trung & Bùi Hải Lê (2009), biến việc giải hệ phương trình vi phân liên hợp phức tạp thành bài toán giải tích xác định, cho phép tìm nghiệm tối ưu toàn cục (Global Optimum) trong không gian trạng thái.

2. Sự đổi mới về phương pháp luận so với các nghiên cứu quốc tế trước đây?

So sánh với các nghiên cứu quốc tế:

  • So với Atanackovic & Simic (1999) và Glavardanov & Atanackovic (2001) vốn sử dụng phương pháp Ritz hoặc tích phân số trực tiếp giải phương trình vi phân phi tuyến bậc cao dễ gặp bất ổn định số, luận án đã tích hợp thành công Phương pháp Ma trận Truyền (TMM) vào thuật toán lặp Hamilton.
  • So với Braun (2008) và Minh-Quy Le et al. (2012) vốn chỉ giải bài toán đơn mục tiêu hoặc ràng buộc thể tích cố định, phương pháp luận của luận án cho phép quét toàn bộ Không gian Trạng thái mở rộng để kiến tạo Tập giải pháp khả thi (Feasible Region) và Mặt trận Pareto hoàn chỉnh đa mục tiêu.

3. Phát hiện bất ngờ và phản trực giác nhất thu được từ dữ liệu số là gì?

Phát hiện bất ngờ nhất là: Cấu hình tối ưu của kết cấu không phụ thuộc vào giá trị tải trọng kích động mà phụ thuộc hoàn toàn vào dạng dao động riêng (eigenmode) của hệ. Cụ thể, việc tăng độ cứng động của trục không đòi hỏi phải tăng thể tích toàn bộ trục, mà chỉ cần bổ sung vật liệu chính xác tại các điểm nút dao động (nơi mô men và lực cắt cực đại) và đồng thời giảm mạnh tiết diện tại các điểm bụng dao động. Kết quả này giúp giảm tổng khối lượng kết cấu tới $15.4%$ nhưng tần số riêng dao động lại tăng thêm $22.8%$, phá vỡ định kiến trực giác thông thường rằng muốn tăng độ cứng thì phải tăng khối lượng.

4. Luận án có cung cấp quy trình tái lặp (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp quy trình tái lập vô cùng chi tiết và minh bạch tại Phụ lục A, B, C, D:

  • Toàn bộ sơ đồ thuật toán dòng (flowchart) giải bài toán điều khiển tối ưu đa mục tiêu PMP.
  • Hệ phương trình ma trận truyền chuẩn hóa cho từng dạng liên kết (2.15 - 2.26).
  • Quy trình đối chuẩn trên ANSYS với việc chỉ định rõ ràng loại phần tử (SOLID186), điều kiện hội tụ lưới và thiết lập bài toán Modal / Harmonic Response. Các nhà nghiên cứu độc lập có thể lập trình lại hoàn toàn trên MATLAB dựa trên các công thức giải tích được công bố.

5. Lộ trình nghiên cứu 10 năm được phác thảo cụ thể như thế nào?

Lộ trình 10 năm được cấu trúc qua 3 giai đoạn:

  • Giai đoạn 1 (Năm 1 - 3): Mở rộng thuật toán PMP cho dao động uốn - xoắn kết hợp của dầm thành mỏng phi đối xứng và kết cấu FGM.
  • Giai đoạn 2 (Năm 4 - 6): Phát triển hệ thống điều khiển dao động chủ động (Active Vibration Control) tích hợp vật liệu áp điện (Piezoelectric) và bộ giảm chấn TMD điều khiển thông minh.
  • Giai đoạn 3 (Năm 7 - 10): Xây dựng nền tảng phần mềm tối ưu hóa hình học tự động kết nối API với các phần mềm CAD/CAM/CAE và tích hợp công nghệ in 3D kim loại phục vụ chế tạo các kết cấu có biên dạng tối ưu phức tạp trong công nghiệp hàng không vũ trụ.

Kết luận

Luận án tiến sĩ "Ứng dụng lý thuyết điều khiển trong tối ưu tần số riêng và khối lượng của kết cấu thanh" là một công trình nghiên cứu khoa học công phu, chuẩn mực và mang tính đột phá cao trong lĩnh vực Cơ học Kỹ thuật và Điều khiển Tối ưu Kết cấu. Năm đóng góp cốt lõi của luận án bao gồm:

  1. Xây dựng thành công hệ hàm đa mục tiêu tổng quát cho bài toán tối ưu đồng thời dải tần số riêng và khối lượng kết cấu thanh, trục.
  2. Chứng minh giải tích và xác định tường minh hệ số tỷ lệ liên hợp $k$, tạo bước nhảy vọt về hiệu năng giải thuật PMP.
  3. Thiết lập bộ công cụ giải tích số kết hợp Phương pháp Ma trận Truyền và thuật toán lặp Hamilton trên MATLAB, được kiểm chứng chuẩn xác bằng ANSYS.
  4. Phát hiện và chứng minh quy luật cấu hình tối ưu định tính và quy luật cấu hình tối ưu tương đương giữa các hệ kết cấu có điều kiện biên và chiều dài khác nhau.
  5. Kiến tạo hoàn chỉnh Tập giải pháp khả thi và Mặt trận Pareto đa mục tiêu, định lượng hóa sự đánh đổi giữa tần số dao động và khối lượng vật liệu.

Công trình không chỉ làm phong phú thêm kho tàng lý thuyết điều khiển tối ưu cơ học mà còn mở ra các hướng nghiên cứu hiện đại về kiểm soát dao động và chẩn đoán kết cấu, khẳng định vị thế khoa học của nghiên cứu cơ học Việt Nam trên trường quốc tế với các bài báo công bố trên tạp chí ISI Meccanica danh giá.