Tổng quan về luận án

Nghiên cứu dao động và điều khiển chủ động cho các công trình kỹ thuật quy mô lớn chịu tác động của tải trọng ngẫu nhiên môi trường (như gió bão, động đất, sóng biển) là một trong những thách thức cốt lõi của cơ học kết cấu hiện đại. Luận án tiến sĩ chuyên ngành Cơ học Kỹ thuật với đề tài "Phát triển thuật toán điều khiển tích cực phản hồi cho các kết cấu trong điều kiện đo hạn chế đáp ứng" đặt nền móng giải quyết bài toán điều khiển tối ưu trong điều kiện hạn chế kép: hạn chế đo lường (limited sensing) và hạn chế kích động (limited actuation). Trong thực tế công trình, số lượng cảm biến (sensors) và máy kích động (actuators) luôn nhỏ hơn rất nhiều so với bậc tự do hay số chiều trạng thái thực tế của kết cấu ($n \gg n_u, n \gg n_y$). Luận án nhấn mạnh một thực tế vận hành quan trọng: "Sự hạn chế số lượng của máy kích động và bộ cảm biến không phải chỉ do những yêu cầu về chi phí mà nó còn ảnh hưởng đến độ tin cậy của điều khiển vì số lượng càng tăng thì nguy cơ hỏng hóc càng cao."

Khoảng trống học thuật (research gap) trọng tâm nằm ở chỗ các thuật toán kinh điển như Điều chỉnh Tuyến tính Toàn phương (Linear Quadratic Regulator - LQR) và Điều khiển Tuyến tính Toàn phương Gauss (Linear Quadratic Gaussian - LQG) vốn dựa trên giả thiết kích động ngoài là quá trình ngẫu nhiên ồn trắng (white noise) hoặc bằng 0, dẫn đến luật điều khiển thuần túy phản hồi (pure feedback) không triệt tiêu được tải trọng tiền định hoặc tải trọng tần số thấp. Hơn nữa, việc cố gắng tăng độ lợi phản hồi (gain matrix $G$) để bù đắp năng lượng kích động sẽ trực tiếp gây ra hiện tượng khuyếch đại thời gian trễ ($\tau_u$) và hiện tượng tràn dạng riêng (spillover effect), dẫn đến mất ổn định hệ động lực.

Để giải quyết triệt để vấn đề này, hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để ước lượng trực tuyến thành phần kích động ngoài ngẫu nhiên chưa biết mà không cần lắp đặt thêm cảm biến đo tải trọng trực tiếp?
  • Câu hỏi nghiên cứu 2 (RQ2): Cấu trúc luật điều khiển tích hợp hồi tiếp - dẫn tiếp (Feedback-Feedforward - FB-FF) có thể khắc phục triệt để sai số tối ưu của LQR/LQG trong điều kiện kích động tần số thấp hay không?
  • Câu hỏi nghiên cứu 3 (RQ3): Tiêu chí toán học định lượng nào cho phép xác định vị trí tối ưu của máy kích động và đầu đo nhằm cực tiểu hóa sai số năng lượng điều khiển?
  • Giả thuyết 1 (H1): Thành phần dẫn tiếp $u_{FF}(t)$ xấp xỉ từ phương trình vi phân ngược thời gian kết hợp thuật toán nhận dạng quy nạp trễ thời gian $\Delta\tau$ sẽ thu hẹp sai số lực điều khiển $|e_p(t)|$ so với nghiệm tối ưu thực sự.
  • Giả thuyết 2 (H2): Việc tích hợp thành phần dẫn tiếp nhận dạng trực tuyến sẽ duy trì hiệu quả giảm chấn gia tốc tương đương hoặc vượt trội so với LQR nhưng giảm thiểu từ 30% đến 50% chuyển vị của thiết bị giảm chấn khối lượng tích cực (Active Mass Damper - AMD) và năng lượng sinh lực của actuator.

Khung lý thuyết của luận án được xây dựng trên nền tảng lý thuyết không gian trạng thái (State-Space Formulation), quy hoạch động (Bellman, 1957), nguyên lý cực đại Pontryagin (Pontryagin và cộng sự, 1967), lý thuyết lọc tối ưu Kalman-Bucy (Kalman và Bucy, 1961), và thuật toán nhận dạng ngoại lực trực tuyến (Anh, 2000). Luận án được thực nghiệm mô phỏng trên các mô hình kết cấu phức tạp: từ hệ dao động ô tô 2 bậc tự do chịu mấp mô mặt đường điều hòa, mô hình động lực học bay phi tuyến tuyến tính hóa 4 trạng thái của máy bay, kết cấu dầm đứng/ngang chịu sóng biển, cho đến hệ khung không gian giàn khoan biển cao 39m, nặng 500 tấn (gồm 95 nút, 196 phần tử) lắp đặt AMD 20 tấn dưới tải trọng sóng chu kỳ 11s và chiều cao sóng 15m. Luận án đã khẳng định tính đột phá trong việc nâng cao độ bền vững (robustness), giảm tiêu hao năng lượng và triệt tiêu nguy cơ bất ổn định do trễ pha thời gian.

Literature Review và Positioning

Lý thuyết điều khiển dao động công trình trải qua quá trình tiến hóa từ các phương pháp cổ điển trên miền tần số (biến đổi Laplace, đồ thị Nyquist 1932, biểu đồ Bode 1945, quỹ tích nghiệm Root Locus của Evans 1948) sang lý thuyết điều khiển hiện đại trên miền thời gian dựa trên ma trận trạng thái từ thập niên 1960. Các công trình nền tảng của Kwakernaak và Sivan (1972), Abdel-Rohman và Leipholz (1978), Kailath (1980), cùng Suhardjo và cộng sự (1992) đã thiết lập chuẩn mực cho điều khiển LQR và $H_2$. Tuy nhiên, hai luồng quan điểm đối lập lớn luôn tồn tại trong trường phái điều khiển kết cấu:

Luồng quan điểm thứ nhất (Thuần phản hồi - Pure Feedback): Đại diện bởi các nghiên cứu LQR/LQG truyền thống và bộ quan sát Luenberger (1966). Luồng quan điểm này giả định kích động ngoài là quá trình ngẫu nhiên ồn trắng không thể đo lường, do đó toàn bộ tín hiệu điều khiển chỉ được tổng hợp từ sai lệch trạng thái hồi tiếp ($u = -Gx$). Ưu điểm là thuật toán hoàn toàn độc lập với nguồn kích động, nhưng nhược điểm cốt tử là phản ứng trễ (chỉ tác động sau khi kết cấu đã biến dạng) và hoàn toàn bất lực trong việc triệt tiêu tải trọng kích động tiền định mang tần số thấp.

Luồng quan điểm thứ hai (Điều khiển thích ứng và Dẫn tiếp - Feedforward/Adaptive Control): Đại diện bởi các kỹ thuật lọc thích ứng Filtered-x LMS (Burgess 1981; Elliott và cộng sự, 1987; Vipperman và cộng sự, 1993) hoặc Higher Harmonic Control (Shaw 1980; Kewley và cộng sự, 1995). Nhóm tác giả này đề xuất đưa trực tiếp tín hiệu ngoại lực vào luật điều khiển để đặt các điểm "không" (transmission zeros) triệt tiêu kích động. Tuy nhiên, rào cản thực nghiệm lớn nhất của nhóm này là bắt buộc phải đo đạc trực tiếp kích động đầu vào $f(t)$ — một điều kiện bất khả thi đối với trường sóng biển ngẫu nhiên, áp lực gió gián đoạn hay gia tốc nền địa chấn đa chiều. Khi chuyển sang mô hình hóa ồn màu (colored noise) trong không gian trạng thái mở rộng, sự phụ thuộc quá mức vào hàm mật độ phổ công suất giả định làm suy giảm nghiêm trọng tính thích ứng khi môi trường thay đổi.

                   TIẾN TRÌNH TIẾP CẬN ĐIỀU KHIỂN KẾT CẤU

Vị thế nghiên cứu của luận án được xác lập bằng việc bắc cầu qua hai luồng tư tưởng: sử dụng cấu trúc toán học của lý thuyết điều khiển hiện đại nhưng giải quyết bài toán dẫn tiếp mà hoàn toàn không cần cảm biến đo ngoại lực, thông qua thuật toán nhận dạng kích động trực tuyến thời gian thực kế thừa và phát triển từ Nguyễn Đông Anh (2000).

So sánh với các nghiên cứu và dự án quốc tế kinh điển:

  1. Tòa nhà Kyobashi Seiwa (Tokyo, Nhật Bản, 1989): Công trình cao 11 tầng (33m, nặng 400 tấn) lắp đặt 2 khối lượng AMD thủy lực (4 tấn cho dao động ngang chu kỳ 1.13s và 1 tấn cho dao động xoắn chu kỳ 0.67s) sử dụng thuật toán hồi tiếp vận tốc (Kajima, 1991). Luận án chứng minh rằng việc bổ sung kênh dẫn tiếp nhận dạng lực giúp giảm thiểu xung lực đỉnh của máy kích động thủy lực tốt hơn đáng kể so với hệ thống của Kajima.
  2. Tháp truyền hình Nam Kinh (Nanjing TV Tower, Trung Quốc, 1998): Tháp cao 340m trang bị AMD dạng đai trượt tròn nặng 59 tấn trên đệm cao su Teflon với 3 máy kích động thủy lực góc $120^\circ$ (Reinhorn và cộng sự, 1998; Cao và cộng sự, 1998). Thuật toán của luận án trực tiếp giải quyết hiện tượng trễ điều khiển và tràn mode bậc cao vốn là nguyên nhân gây suy giảm hiệu năng tại tháp Nam Kinh khi chịu tải gió bão giật.
  3. Tháp Shinjuku Park Tower (Nhật Bản, 1994): Tòa nhà sử dụng AMD con lăn trên ray chữ V nặng 110 tấn truyền động mô tơ điện (Takenaka và cộng sự, 1994), đạt mức giảm gia tốc đỉnh từ 6.8 Gal xuống 2.1 Gal. Mô hình nâng cao của luận án cung cấp cơ sở toán học để đạt mức suy giảm tương đương với hành trình con lắc/con lăn ngắn hơn 25-40%.

Đó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 và thách thức các giới hạn cố hữu của phương trình đại số ma trận Riccati trong điều khiển tối ưu LQR kinh điển. Khi thiết lập bài toán cực tiểu hóa phiếm hàm tích phân toàn phương có trọng số: $$J = \frac{1}{2} \int_{0}^{\infty} \left( x^T Q x + u^T R u \right) dt$$ với ràng buộc phương trình trạng thái $\dot{x}(t) = Ax(t) + Bu(t) + Hf(t)$, việc áp dụng phép tính biến phân và hàm Hamilton $H_{mt}$ đã dẫn tới hệ điều kiện tối ưu phân rã nhân tử Lagrange $\lambda(t) = Px(t) + p(t)$.

Luận án đã chứng minh một cách tường minh rằng:

  1. Ma trận $P$ là nghiệm duy nhất xác định dương của phương trình đại số Riccati độc lập với kích động: $$PA - PBR^{-1}B^T P + A^T P + Q = 0$$
  2. Thành phần dẫn tiếp $p(t)$ thỏa mãn phương trình vi phân ngược chiều thời gian với điều kiện biên tại vô cùng: $$\dot{p}(t) + (G^T B^T - A^T)p(t) + PHf(t) = 0, \quad p(\infty) = 0$$ trong đó ma trận độ lợi hồi tiếp $G = R^{-1}B^T P$.

Đóng góp đột phá về mặt lý thuyết của luận án là xác lập điều kiện biên giải tích cho phép xấp xỉ thành phần dẫn tiếp khi kích động có tần số thấp ($\dot{p}(t) \approx 0$): $$\tilde{p}(t) = -(G^T B^T - A^T)^{-1} PHf(t)$$ Sai số lực điều khiển $E_{FBFF} = |R^{-1}B^T (p(t) - \tilde{p}(t))| = |R^{-1}B^T e_p(t)|$ được chứng minh giải tích luôn nhỏ hơn sai số của thuật toán LQR kinh điển $E_{LQR} = |R^{-1}B^T p(t)|$ khi và chỉ khi đại lượng phi thứ nguyên $\kappa_f$ thỏa mãn bất đẳng thức hội tụ: $$\kappa_f = \frac{\max_t |(G^T B^T - A^T)^{-1} P H \dot{f}(t)|}{\max_t |PHf(t)|} < 1$$

Đây là một bước chuyển biến quan trọng (paradigm shift), khẳng định rằng việc tích hợp thành phần dẫn tiếp xấp xỉ không đòi hỏi tri thức tiên nghiệm về toàn bộ lịch sử kích động mà chỉ cần ràng buộc về tỷ số biên độ - đạo hàm tần số thấp của ngoại lực.

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

Khung phân tích của luận án tích hợp đồng thời ba trụ cột lý thuyết lớn:

  • Lý thuyết điều khiển ngẫu nhiên hiện đại (Modern Stochastic Control Theory)
  • Lý thuyết nhận dạng hệ thống thời gian thực (Real-time System Identification)
  • Lý thuyết động lực học công trình và phân tích dạng riêng (Structural Dynamics & Modal Analysis)
                            KHUNG PHÂN TÍCH ĐỘC ĐÁO TÍCH HỢP

Điểm sáng tạo cốt lõi là việc hóa giải hiện tượng "vòng lặp đại số" (algebraic loop). Trong phương trình nhận dạng: $$Hf(t) = \dot{x}(t) - Ax(t) - Bu(t) = \dot{x}(t) - (A - BG)x(t) - Bu_{FF}(t)$$ Đại lượng $u_{FF}(t)$ xuất hiện đồng thời ở cả hai vế (vừa là biến cần tìm, vừa tham gia xác định $Hf(t)$). Luận án giải quyết bằng công thức quy nạp trễ thời gian bước ngắn $\Delta\tau$: $$u_{FF}(t) = -R^{-1}B^T (G^T B^T - A^T)^{-1} P \left[ \dot{x}(t - \Delta\tau) - (A - BG)x(t - \Delta\tau) - Bu_{FF}(t - \Delta\tau) \right]$$ Luận án chứng minh tính đúng đắn dựa trên định lý số gia trung bình Lagrange: $$u_{FF}(t) = -R^{-1}B^T (G^T B^T - A^T)^{-1} P \left[ Hf(t) - \Delta\tau H\dot{f}(\xi) \right], \quad \xi \in (t - \Delta\tau, t)$$ Luận án khẳng định rõ: "Lực điều khiển tính theo biểu thức quy nạp sẽ đủ gần với lực điều khiển trong trường hợp đo được kích động ngoài. Điều này có được cũng là do kích động ngoài không phụ thuộc đáp ứng của hệ nên đại lượng $\dot{f}(\xi)$ không bị khuyếch đại trong quá trình điều khiển."

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

Thiết kế nghiên cứu

Nghiên cứu quán triệt thế giới quan thực chứng (positivism) kết hợp mô hình hóa giải tích chặt chẽ và mô phỏng số ngẫu nhiên (stochastic numerical simulation). Cơ sở phương pháp luận là sự phân rã đa tầng (multi-level design) thông qua Nguyên lý tách (Separation Principle):

  • Tầng 1 (Chương 2): Khảo sát bài toán lý tưởng không hạn chế đo trạng thái ($y(t) = x(t)$) nhưng hạn chế vị trí đặt lực actuator ($B \in \mathbb{R}^{n \times n_u}$).
  • Tầng 2 (Chương 3): Khảo sát bài toán không hạn chế đặt lực nhưng hạn chế đầu đo ($C \in \mathbb{R}^{n_y \times n}$, sử dụng bộ lọc Kalman-Bucy).
  • Tầng 3 (Chương 4): Tổng hợp bài toán thực tế phản hồi đầu ra tổng quát (Output Feedback Problem) chịu tác động đồng thời của nhiễu quá trình $w(t)$ và nhiễu đo lường $v(t)$, từ đó đề xuất 2 phiên bản thuật toán nâng cao.

Phương pháp tiếp cận toàn diện được thể hiện qua mô tả phương trình chuyển động vi phân cấp hai chuyển đổi sang hệ phương trình trạng thái vi phân cấp một: $$\dot{x}(t) = Ax(t) + Bu(t) + Hf(t), \quad x(0) = x_0$$ trong đó: $$x = \begin{bmatrix} d \ \dot{d} \end{bmatrix}, \quad A = \begin{bmatrix} 0 & I \ -M^{-1}K & -M^{-1}D \end{bmatrix}, \quad B = \begin{bmatrix} 0 \ M^{-1}L_u \end{bmatrix}, \quad H = \begin{bmatrix} 0 \ M^{-1}L_f \end{bmatrix}$$

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ các giao thức kiểm định nghiêm ngặt:

  1. Thiết lập mô hình hình học và ma trận phần tử hữu hạn: Sử dụng SAP2000 để trích xuất ma trận khối lượng $M$, ma trận độ cứng $K$, ma trận cản Rayleigh $D = \alpha M + \beta K$.
  2. Thu gọn bậc tự do bằng phân tích dạng riêng (Modal Decomposition): Chuyển dịch vật lý $d(t)$ được biểu diễn qua ma trận mode chính chuẩn $\Phi$: $d(t) = \Phi q(t)$, biến đổi hệ phương trình thành các phương trình tọa độ suy rộng độc lập cho 4 dạng dao động cơ bản đầu tiên nhằm kiểm soát hiện tượng tràn mode (spillover).
  3. Mô phỏng ngẫu nhiên Monte Carlo: Tạo lập trường sóng ngẫu nhiên theo phổ JONSWAP/Pierson-Moskowitz và kích động mặt đường ngẫu nhiên theo tiêu chuẩn ISO 8608.
  4. Phân tích ổn định và trễ pha: Khai triển chuỗi Taylor cho phương trình đặc trưng siêu việt có trễ $\tau_u$: $$\det\left( A - BGe^{-\lambda\tau_u} - \lambda I_n \right) = \det\left( A - BG + \tau_u BG\lambda - \frac{\tau_u^2}{2}BG\lambda^2 + \dots - \lambda I_n \right) = 0$$ để chứng minhเชิง giải tích nguy cơ dịch chuyển nghiệm trị riêng $\lambda$ sang nửa phải mặt phẳng phức khi chỉ số ma trận $G$ quá lớn.

Data và phân tích

Toàn bộ thuật toán được lập trình, giải số và kiểm chứng độc lập trên nền tảng MATLAB và công cụ SIMULINK.

Các thông số động lực học chi tiết trong các bài toán thử nghiệm bao gồm:

  • Hệ thống treo ô tô (Muller và cộng sự, 1997):

    • Khối lượng thân xe $m_1 = 1200\text{ kg}$, khối lượng trục bánh $m_2 = 80\text{ kg}$.
    • Độ cứng lò xo $k_1 = 3\times 10^4\text{ N/m}$, độ cứng lốp $k_2 = 32\times 10^4\text{ N/m}$.
    • Hệ số cản giảm chấn $c_1 = 4800\text{ Ns/m}$.
    • Trọng số tối ưu: $Q = \text{diag}(k_1, k_1+k_2, m_1, m_2)$, $R = r = 5\times 10^{-5}$.
    • Kích động mặt đường sóng điều hòa: $y_s = Y \sin(\frac{\pi v t}{L})$ với biên độ $Y = 0.50\text{ m}$, vận tốc $v = 50\text{ km/h}$, chu kỳ kiểm tra bước sóng: $L = 13.89\text{ m}$, $L = 6.94\text{ m}$, $L = 4.63\text{ m}$.

  • Hệ kết cấu khung không gian ngoài khơi (Offshore Space Frame):

    • Chiều cao công trình $H = 39\text{ m}$, tổng khối lượng kết cấu $M_{total} = 500\text{ tấn}$.
    • Số lượng phần tử hữu hạn: 95 nút, 196 thanh phần tử không gian.
    • Thiết bị AMD: khối lượng $m_d = 20\text{ tấn}$ (chiếm 4% khối lượng kết cấu), lắp đặt tại đỉnh khung.
    • Tải trọng sóng: Chiều cao sóng $H_w = 15\text{ m}$, chu kỳ sóng $T_w = 11\text{ s}$ (tần số thấp $\approx 0.09\text{ Hz}$).

  • Kiểm định độ bền vững (Robustness Checks):

    • Nhiễu đo lường Gauss $v(t)$ với ma trận hiệp phương suất phổ $I_v = 10^{-10}\text{ m}^2/\text{s}$ và $I_v = 10^{-12}\text{ m}^2/\text{s}$.
    • Thời gian trễ máy kích động kiểm tra: $\tau_u = 5\times 10^{-4}\text{ s}$ đến $\tau_u = 25\times 10^{-4}\text{ s}$.

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

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

  1. Hiện tượng suy giảm hiệu năng nghịch lý của LQR dưới kích động chu kỳ: Tại bước sóng mặt đường $L = 4.63\text{ m}$, giá trị chỉ tiêu tích phân toàn phương $J$ của hệ thống điều khiển LQR đạt mức $0.95\text{ Js}$, trong khi trường hợp không điều khiển chỉ tiêu này là $0.41\text{ Js}$. Kết quả thực nghiệm này khẳng định: khi kích động ngoài khác 0, luật LQR thuần phản hồi không những không tối ưu mà còn làm tăng năng lượng dao động do xung đột pha giữa lực điều khiển phản hồi và ngoại lực. Ngược lại, thuật toán điều khiển tối ưu thực sự và thuật toán đề xuất FB-FF duy trì chỉ số $J$ ở mức cực tiểu ổn định.
Bước sóng $L$ Trạng thái điều khiển Biến dạng nhíp RMS (m) Gia tốc thân xe RMS ($\text{m/s}^2$) Lực điều khiển RMS (kN) Chỉ tiêu tối ưu $J$ (Js)
$L = 13.89\text{ m}$ Không điều khiển 0.042 1.85 0.00 0.42
LQR kinh điển 0.021 0.89 1.45 0.22
Tối ưu thực sự 0.018 0.72 1.32 0.18
$L = 6.94\text{ m}$ Không điều khiển 0.065 2.45 0.00 0.68
LQR kinh điển 0.038 1.42 2.10 0.35
Tối ưu thực sự 0.028 1.10 1.85 0.25
$L = 4.63\text{ m}$ Không điều khiển 0.052 2.10 0.00 0.41
LQR kinh điển 0.071 2.85 3.65 0.95
Tối ưu thực sự 0.031 1.15 2.05 0.28

  1. Tiết kiệm biên độ hành trình và công sinh lực của thiết bị AMD: Trong bài toán kết cấu khung giàn khoan 39m chịu tải sóng $H_w = 15\text{ m}$, thuật toán hồi tiếp - dẫn tiếp (FB-FF) kết hợp nhận dạng lực đạt mức giảm gia tốc đỉnh tương đương LQR (giảm trên 65% so với không điều khiển), nhưng giảm được 38.4% chuyển vị tương đối của khối lượng AMD và giảm 32.1% giá trị đỉnh lực điều khiển thủy lực.

  2. Cơ chế khuyếch đại sai số do trễ thời gian trong mô hình máy bay: Khi khảo sát động lực học bay tự động, việc giảm hệ số trọng số điều khiển $r$ từ $2\times 10^{-4}$ xuống $2\times 10^{-6}$ nhằm tăng tốc độ dập tắt độ lệch góc nâng $\theta$ đã dẫn đến sự mất ổn định phân kỳ (divergence instability) khi hệ thống xuất hiện độ trễ $\tau_u = 0.05\text{ s}$. Nghiệm trị riêng $\lambda$ bị đẩy sang bán mặt phẳng bên phải do ma trận $G$ quá lớn. Việc bổ sung thành phần dẫn tiếp cho phép giữ nguyên $r$ ở ngưỡng an toàn nhưng vẫn dập tắt dao động nhanh hơn 45%.

  3. Hiệu quả nhận dạng kích động ngẫu nhiên trực tuyến: Sai số nhận dạng tải sóng $Hf(t)$ từ chuyển vị tương đối đo được tại vị trí $l/2$ đạt độ chính xác trên 92% về cả phần thực và phần ảo trong miền tần số, ngay cả khi dữ liệu đo bị ô nhiễm bởi nhiễu trắng với mật độ phổ $I_v = 10^{-10}\text{ m}^2/\text{s}$.

Implications đa chiều

  • Về mặt lý thuyết: Luận án mở rộng định lý tách trong điều khiển ngẫu nhiên từ bài toán ước lượng trạng thái thuần túy sang bài toán đồng thời ước lượng trạng thái và nhận dạng kích động ngẫu nhiên ngoại cảnh.
  • Về mặt phương pháp luận: Thiết lập thuật toán quy nạp thời gian trễ $\Delta\tau$ như một giải pháp tổng quát để giải quyết triệt để vấn đề vòng lặp đại số trong mô phỏng số hệ động lực liên tục - rời rạc.
  • Về mặt ứng dụng thực tiễn: Cung cấp hướng dẫn tính toán trực tiếp cho kỹ sư thiết kế hệ thống giảm chấn AMD cho nhà cao tầng và công trình biển, cho phép thu nhỏ hành trình pit-tông và giảm công suất trạm nguồn thủy lực, từ đó tiết kiệm chi phí đầu tư thiết bị.
  • Về mặt chính sách và tiêu chuẩn: Đóng góp cơ sở khoa học để hoàn thiện các tiêu chuẩn thiết kế kháng chấn chủ động cho công trình cấp đặc biệt tại Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả lý thuyết và mô phỏng xuất sắc, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Giả thiết tuyến tính hóa hệ thống: Nghiên cứu tập trung vào hệ phương trình vi phân tuyến tính bất biến theo thời gian (LTI). Các phi tuyến hình học lớn và phi tuyến vật liệu của kết cấu khi chịu tải trọng động đất cấp thảm họa chưa được tích hợp trong phương trình Riccati vi phân.
  2. Đặc tính động lực học chi tiết của cơ cấu chấp hành (Actuator Dynamics): Luận án mô hình hóa máy kích động dưới dạng nguồn lực tức thời có trễ thuần túy ($\tau_u$), chưa xét đến hiện tượng tương tác phi tuyến kết cấu - máy kích động thủy lực (Structure-Actuator Interaction - SAI) như bão hòa áp suất (saturation) và bão hòa hành trình.
  3. Phổ tần số của kích động: Thuật toán FB-FF phát huy hiệu quả tối ưu nhất trong miền kích động tần số thấp (sóng biển, gió bão), hiệu quả suy giảm khi áp dụng cho động đất trường gần có dải tần rộng và tính chất xung kích gián đoạn cao.
  4. Kiểm chứng thực nghiệm quy mô lớn: Các kết quả mới dừng lại ở mức mô phỏng số (numerical simulations) trên mô hình phần tử hữu hạn và thông số công trình chuẩn, chưa được kiểm chứng trên bàn rung động học (shake table testing) hoặc mô hình kết cấu vật lý thực tế.

Chương trình nghiên cứu tương lai:

  • Mở rộng thuật toán FB-FF cho hệ phi tuyến bằng phương pháp tuyến tính hóa thống kê tương đương hoặc điều khiển trượt mode (Sliding Mode Control).
  • Tích hợp động lực học phi tuyến của vật liệu thông minh như chất lỏng từ biến (Magnetorheological - MR Damper) vào kênh điều khiển dẫn tiếp.
  • Nghiên cứu thuật toán điều khiển rời rạc thích ứng trên nền tảng vi xử lý DSP thời gian thực phục vụ chế tạo thử nghiệm pilot.

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

Luận án tạo ra tác động học thuật sâu sắc với khả năng trích dẫn cao trong các tạp chí chuyên ngành hàng đầu về cơ học kết cấu và điều khiển tự động (Journal of Structural Engineering, Structural Control and Health Monitoring, Journal of Sound and Vibration). Luận án mở ra hướng chuyển giao công nghệ tiềm năng cho ngành công nghiệp đóng tàu, dầu khí ngoài khơi (giàn khoan fixed jacket và semi-submersible) và xây dựng dân dụng nhà siêu cao tầng tại các đô thị chịu rủi ro bão nhiệt đới và động đất.

Về mặt xã hội và an toàn công trình, các thuật toán nâng cao giúp bảo vệ tính mạng con người, giảm thiểu thiệt hại kinh tế do mỏi kết cấu và gia tăng đáng kể tiện nghi vi khí hậu (human comfort) cho người sử dụng công trình cao tầng bằng cách triệt tiêu gia tốc lắc ngang xuống dưới ngưỡng cảm nhận tiêu chuẩn ISO 2631.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ học / Tự động hóa: Thụ hưởng khung phân tích toán học giải tích chặt chẽ, phương pháp xấp xỉ phương trình vi phân ngược thời gian và giải thuật khử vòng lặp đại số.
  • Kỹ sư R&D Công trình biển & Xây dựng: Ứng dụng trực tiếp quy trình lựa chọn vị trí cảm biến tối ưu $\min |e_p(t)|$ và thuật toán điều khiển nâng cao để tối ưu hóa thiết kế AMD, giảm thể tích phòng máy kỹ thuật trên tầng thượng.
  • Các nhà phát triển công nghệ cảm biến và thiết bị thông minh: Cơ sở tích hợp thuật toán nhận dạng tải trọng trực tuyến vào các hệ thống quan trắc sức khỏe công trình (Structural Health Monitoring - SHM).

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 mở rộng Lý thuyết Điều khiển Tuyến tính Toàn phương (LQR/LQG) sang cấu trúc Hồi tiếp - Dẫn tiếp tích hợp (FB-FF) mà không cần đo trực tiếp kích động đầu vào. Luận án đã thiết lập điều kiện hội tụ giải tích $\kappa_f < 1$ cho phép xấp xỉ nghiệm của phương trình vi phân Riccati ngược thời gian thông qua các chuẩn ma trận của hệ thống, phá vỡ bế tắc lịch sử đòi hỏi phải biết trước toàn bộ lịch sử kích động ngẫu nhiên trong tương lai.

2. Đột phá về phương pháp luận so với các công trình nghiên cứu quốc tế trước đây?

So với các nghiên cứu lọc thích ứng (như Filtered-x LMS của Burgess 1981 đòi hỏi cảm biến đo ngoại lực) và các nghiên cứu LQR kinh điển (Kwakernaak 1972 hoàn toàn bỏ qua dẫn tiếp), luận án đã phát triển phương pháp giải số quy nạp có trễ bước ngắn $\Delta\tau$ để nhận dạng trực tuyến ngoại lực từ dữ liệu đo đáp ứng, loại bỏ hoàn toàn vòng lặp đại số (algebraic loop) mà vẫn đảm bảo tính ổn định tiệm cận của hệ động lực.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu mô phỏng?

Phát hiện phản trực giác nhất là việc áp dụng thuật toán điều khiển tối ưu kinh điển LQR trong điều kiện có kích động ngoài chu kỳ ($L = 4.63\text{ m}$ trong bài toán ô tô) có thể làm tăng chỉ tiêu năng lượng tổng thể $J$ lên hơn 2.3 lần so với khi không điều khiển ($0.95\text{ Js}$ so với $0.41\text{ Js}$). Điều này chứng minh rằng việc áp dụng máy móc các thuật toán điều khiển phản hồi mà không xét đến kênh dẫn tiếp có thể gây nguy hiểm cho kết cấu công trình.

4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) hoàn chỉnh không?

Có. Toàn bộ hệ phương trình trạng thái, các ma trận khối lượng $M$, độ cứng $K$, cản $D$, thông số hình học 95 nút và 196 phần tử khung không gian, cùng các bảng trọng số $Q, R$ và mã nguồn thuật toán giải ma trận Riccati bằng phương pháp chéo hóa (Phụ lục C) đều được chuẩn hóa chi tiết, cho phép tái lập 100% kết quả trên môi trường MATLAB/SIMULINK.

5. Lộ trình phát triển nghiên cứu 10 năm tiếp theo được định hình ra sao?

Lộ trình 10 năm tập trung vào: (1) Mở rộng thuật toán cho hệ kết cấu phi tuyến đa bậc tự do chịu tải chấn động phi dừng; (2) Tích hợp học máy (Physics-Informed Neural Networks - PINN) vào khối nhận dạng ngoại lực siêu nhanh; (3) Thử nghiệm thực tế trên hệ thống bàn rung 6 bậc tự do; (4) Thương mại hóa bộ điều khiển nhúng thời gian thực (Embedded Real-time Controller) cho các công trình giàn khoan dầu khí biển sâu.

Kết luận

Luận án tiến sĩ đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 5 đóng góp cốt lõi:

  1. Phát triển thành công thuật toán điều khiển tích cực Hồi tiếp - Dẫn tiếp (FB-FF) khắc phục triệt để sai số cố hữu của các bộ điều khiển kinh điển LQR và LQG dưới tác động của tải trọng ngẫu nhiên.
  2. Đề xuất giải thuật nhận dạng trực tuyến kích động ngoài từ số liệu đo đáp ứng trễ pha ngắn $\Delta\tau$, giải quyết triệt để bài toán vòng lặp đại số trong điều khiển thời gian thực.
  3. Thiết lập tiêu chí toán học chuẩn xác để xác định vị trí đặt lực (Actuator) và vị trí đặt đầu đo (Sensor) tối ưu dựa trên việc cực tiểu hóa chuẩn sai số năng lượng điều khiển $|e_p(t)|$.
  4. Xây dựng 2 phiên bản thuật toán điều khiển nâng cao cho bài toán phản hồi đầu ra tổng quát trong điều kiện đo lường hạn chế và có nhiễu ngẫu nhiên, chứng minh tính bền vững vượt trội.
  5. Kiểm chứng toàn diện trên các mô hình kết cấu phức tạp thực tế (ô tô, máy bay, dầm chịu sóng, khung không gian giàn khoan 39m lắp đặt AMD 20 tấn), chứng minh khả năng giảm thiểu trên 30% năng lượng điều khiển và hành trình thiết bị.

Công trình tạo ra bước tiến quan trọng trong chuyên ngành Cơ học Kỹ thuật và Kỹ thuật Điều khiển Kết cấu, mở ra ba nhánh nghiên cứu mới: điều khiển bền vững kết cấu biển nước sâu, phát triển thiết bị giảm chấn bán tích cực thông minh (semi-active smart dampers), và tích hợp hệ thống quan trắc sức khỏe công trình tự thích ứng theo thời gian thực. Luận án khẳng định vị thế khoa học vững chắc và giá trị ứng dụng thực tiễn to lớn trong việc bảo vệ an toàn cho các công trình hạ tầng trọng điểm quốc gia trước các hiểm họa thiên tai khốc liệt.