Tổng quan luận án

Trong kỹ thuật cơ học và xây dựng công trình, việc kiểm soát dao động của các hệ kết cấu quy mô lớn (như tòa nhà cao tầng, tháp truyền hình, công trình biển, cầu nhịp lớn, phương tiện vận tải và thiết bị hàng không) chịu tác động của môi trường ngẫu nhiên đặt ra nhiều thách thức kỹ thuật. Các mô hình động lực học mô tả kết cấu lớn thường đòi hỏi số lượng bậc tự do và biến trạng thái rất cao. Tuy nhiên, trên thực tế, số lượng máy kích động (actuator) và cảm biến đo lường (sensor) luôn bị giới hạn nghiêm ngặt bởi chi phí đầu tư và yêu cầu đảm bảo độ tin cậy vận hành (do việc gia tăng số lượng thiết bị làm tăng xác suất hỏng hóc của toàn hệ thống). Tình trạng này dẫn đến hai bài toán kỹ thuật điển hình: bài toán hạn chế đo (không thể xác định toàn bộ biến trạng thái của hệ) và bài toán hạn chế đặt lực điều khiển (không thể tác động lực điều khiển trực tiếp lên tất cả các bậc tự do).

Mặt khác, các kích động môi trường từ tự nhiên (như gió bão, sóng biển, động đất hay mấp mô mặt đường) mang đặc tính ngẫu nhiên và hầu như không thể đo đạc trực tiếp một cách đầy đủ trong quá trình vận hành. Các thuật toán điều khiển kinh điển—tiêu biểu là bộ điều chỉnh tuyến tính dạng toàn phương (Linear Quadratic Regulator - LQR) và bộ điều khiển tối ưu tuyến tính ngẫu nhiên Gauss (Linear Quadratic Gaussian - LQG)—thường đưa ra giả thiết kích động ngoài là quá trình ồn trắng (white noise). Giả thiết này cho phép xây dựng luật điều khiển thuần phản hồi không phụ thuộc vào kích động ngoài. Tuy nhiên, khi kích động thực tế có quy luật biến thiên khác ồn trắng (đặc biệt là các kích động mang tần số thấp như tải trọng sóng biển), luật điều khiển thuần phản hồi không đạt được hiệu quả tối ưu thực sự. Đồng thời, việc cố gắng nâng cao hiệu quả điều khiển bằng cách tăng độ lớn các phần tử trong ma trận phản hồi sẽ dẫn tới hiện tượng khuyếch đại sai số tính toán và khuyếch đại ảnh hưởng tiêu cực của thời gian trễ.

Từ các khoảng trống và hạn chế kỹ thuật nêu trên, luận án tiến sĩ 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" được thực hiện. Luận án thuộc chuyên ngành Cơ học kỹ thuật / Điều khiển kết cấu công trình, kết hợp phương pháp giải tích hệ thống và mô phỏng số.

Mục tiêu nghiên cứu tổng quát của luận án được xác định rõ:

  1. Phát triển các thuật toán điều khiển kinh điển cho các kết cấu có sự hạn chế đo, bằng cách khai thác tối đa thông tin đo được từ hệ thống để nhận dạng trực tuyến kích động ngoài tác động lên kết cấu.

Đối tượng nghiên cứu:

  • Các kết cấu cơ học và công trình có nhiều bậc tự do, chịu tải trọng ngẫu nhiên từ môi trường và làm việc trong điều kiện bị hạn chế số lượng cảm biến đo đáp ứng.

Phạm vi nghiên cứu:

  • Không gian nghiên cứu: Khảo sát trên các mô hình kết cấu cơ học và công trình kỹ thuật cụ thể đã được công bố trong các tài liệu chuẩn (benchmark), bao gồm mô hình ô tô chuyển động trên đường mấp mô (2 bậc tự do), chuyển động dọc của máy bay, kết cấu dầm thẳng đứng, dầm ngang chịu tải sóng dập dềnh, khung không gian 95 nút 196 thanh (cao 39m, nặng 500 tấn) chịu tải trọng sóng biển, hệ ăng-ten parabol quan sát bầu trời, và kết cấu công trình nhiều tầng lắp đặt hệ thống giảm chấn khối lượng tích cực (Active Mass Damper - AMD).
  • Thời gian nghiên cứu: Tập trung vào phản ứng động lực học của hệ trong miền thời gian liên tục và các bước rời rạc hóa phục vụ tính toán số trực tuyến; khảo sát trong các khoảng thời gian chịu kích động động lực học hữu hạn (từ các chu kỳ dao động ngắn vài giây đến các tiến trình tải trọng sóng kéo dài).

Tổng quan tài liệu và vị trí của luận án

Tổng quan tài liệu trong luận án phân định rõ tiến trình phát triển của lý thuyết điều khiển từ cổ điển đến hiện đại và sự du nhập của các lý thuyết này vào bài toán kết cấu:

  1. Lý thuyết điều khiển cổ điển: Được phát triển trước những năm 1960 dựa trên công cụ biến đổi Laplace cho các phương trình vi phân tuyến tính, sử dụng các phương pháp đồ thị như Nyquist (1932), Bode (1945), và quỹ tích nghiệm Root Locus của Evans (1948). Các phương pháp này trực quan nhưng bị giới hạn ở các hệ một đầu vào - một đầu ra (SISO) và không đưa ra được luật điều khiển tối ưu cho hệ nhiều biến.
  2. Lý thuyết điều khiển hiện đại và điều khiển tối ưu: Hình thành từ thập niên 1960 cùng sự phát triển của máy tính điện tử, chuyển sang sử dụng không gian trạng thái, đại số ma trận và lý thuyết quá trình ngẫu nhiên. Nền tảng tối ưu hóa dựa trên quy hoạch động của Bellman (1957) và nguyên lý cực đại Pontryagin vcs (1967). Vấn đề ước lượng trạng thái được giải quyết qua bộ quan sát của Luenberger (1966) và bộ lọc Kalman-Bucy (Kalman vcs 1961). Các thuật toán điều khiển tối ưu LQR, LQG trên miền thời gian và chuẩn H2 trên miền tần số được phát triển bởi Abdel-Rohman vcs (1978), Kailath (1980), Kwakernaak vcs (1972), Suhardjo vcs (1992).
  3. Các vấn đề đặc thù trong điều khiển kết cấu:
    • Hiện tượng tràn năng lượng (spillover): Khi thiết kế điều khiển cho một số dạng riêng bậc thấp, năng lượng điều khiển có thể khuyếch đại các dạng riêng bậc cao chưa được mô hình hóa, gây mất ổn định kết cấu (Balas 1978, 1980, 1982; Soong vcs 1982; Seto 1998).
    • Tính bền vững (robustness): Xử lý sai lệch giữa mô hình toán học và kết cấu thực tế dựa trên chuẩn H vô cùng và lý thuyết $\mu$-synthesis (Safonov 1980; Chiang 1988; Grimble 1988).
    • Thời gian trễ (time delay): Độ trễ sinh ra từ quá trình tính toán và đáp ứng của máy kích động có thể làm mất ổn định hệ thống (Basharkhah vcs 1984; Pu 1989).
    • Công nghệ máy kích động và cảm biến: Nghiên cứu cảm biến (Lynch vcs 2006); các van thủy lực đóng mở (Feng vcs 1990; Kurata vcs 1994; Sack vcs 1993); vật liệu thông minh như chất lỏng từ biến MR (Carlson vcs 1996; Spencer vcs 1998a), vật liệu áp điện (Crawley 1987; Ikeda 1997), hợp kim nhớ hình SMA (Wayman vcs 1972; Fremond vcs 1996).
    • Ứng dụng thực tế quy mô lớn: Các công trình lắp đặt hệ điều khiển tích cực AMD được tổng kết bởi Soong (1989), Housner vcs (1997), Spencer vcs (1997), Kitamura vcs (2004), Kareem vcs (1999). Các ví dụ tiêu biểu gồm Kyobashi Seiwa Building (Kajima 1991), Yokohama Landmark Tower (Yamazaki vcs 1992), Applause Tower (2001), Nanjing TV Tower (Reinhorn vcs 1998; Cao vcs 1998), Shinjuku Park Tower (Takenaka vcs 1994).
  4. Các giải pháp bổ sung thành phần dẫn tiếp (feedforward):
    • Phương pháp cổ điển triệt tiêu kích động bằng cách đặt điểm không trùng tần số kích động (Ogata 1970; Sievers vcs 1992).
    • Kỹ thuật lọc thích ứng như Filtered-x LMS (Burgess 1981; Elliott vcs 1987; Vipperman vcs 1993) hoặc Higher Harmonic Control (Kewley vcs 1995; Shaw 1980).
    • Mở rộng vectơ trạng thái chứa kích động ồn màu trong không gian trạng thái.
    • Thuật toán nhận dạng kích động ngoài của tác giả N.D. Anh (2000).
  5. Tình hình nghiên cứu trong nước: Các giáo trình lý thuyết điều khiển tự động của Nguyễn Doãn Phước vcs (2000), Phạm Công Ngô (2001) chủ yếu tập trung vào hệ thống điện hoặc dây chuyền sản xuất tự động. Nghiên cứu điều khiển hệ cơ học được thực hiện cho rôbốt, tàu thủy, cáp treo (Dinh Cong Huan vcs 2004; Nguyen Tran Hiep vcs 2004; Do Sanh vcs 2004; Dang Xuan Hoai vcs 2006; Anh vcs 2004a). Các nghiên cứu điều khiển kết cấu công trình chịu tải trọng động lực học được công bố bởi Anh vcs (2004b, 2005), Ninh Quang Hai vcs (2004), và luận án của Nguyễn Chỉ Sáng (2004) về giảm chấn khối lượng.

Khoảng trống nghiên cứu được luận án lựa chọn giải quyết: Phần lớn các công trình nghiên cứu trong nước chưa giải quyết thấu đáo bài toán hạn chế đo biến trạng thái của kết cấu công trình phức tạp. Mặt khác, thuật toán nhận dạng kích động ngoài của N.D. Anh (2000) ban đầu mới chỉ được thiết lập cho điều kiện lý tưởng: đo được toàn bộ các trạng thái và đặt được lực điều khiển lên toàn bộ các trạng thái. Luận án đặt vị trí nghiên cứu vào việc mở rộng và hoàn thiện thuật toán nhận dạng lực, tích hợp thành phần dẫn tiếp thích nghi vào cấu trúc phản hồi kinh điển nhằm xử lý đồng thời bài toán hạn chế đo và hạn chế đặt lực điều khiển.

Cơ sở lý thuyết và phương pháp nghiên cứu

Luận án sử dụng phương pháp luận của lý thuyết điều khiển hiện đại kết hợp với cơ học kết cấu và mô phỏng số kỹ thuật:

  1. Khung lý thuyết và mô hình toán học:

    • Mô hình hóa hệ dao động nhiều bậc tự do: Phương trình vi phân chuyển động cơ học cấp hai: $$M\ddot{d}(t) + D\dot{d}(t) + Kd(t) = L_u u(t) + L_f f(t)$$ được chuyển đổi tương đương về hệ phương trình không gian trạng thái tuyến tính cấp một: $$\dot{x}(t) = Ax(t) + Bu(t) + Hf(t), \quad x(0) = x_0$$ trong đó $x(t) = [d^T(t), \dot{d}^T(t)]^T$ là vectơ trạng thái $n$ chiều; $A$ là ma trận hệ thống cỡ $n \times n$; $u(t)$ là vectơ lực điều khiển $n_u$ chiều; $f(t)$ là vectơ kích động ngoài $n_f$ chiều; $B$ và $H$ là các ma trận định vị lực điều khiển và kích động.
    • Chỉ tiêu chất lượng tối ưu: Độ đo tích phân dạng toàn phương trên miền thời gian vô hạn: $$J = \frac{1}{2}\int_0^\infty [x^T(t) Q x(t) + u^T(t) R u(t)] dt$$ trong đó $Q$ là ma trận trọng số trạng thái đối xứng thực nửa xác định dương ($n \times n$), $R$ là ma trận trọng số điều khiển đối xứng thực xác định dương ($n_u \times n_u$).
    • Lý thuyết biến phân và phép biến đổi Hamilton: Thiết lập hàm Lagrange mở rộng $L$ thông qua vectơ nhân tử Lagrange $\lambda(t)$. Áp dụng phương pháp tách biến liên hợp $\lambda(t) = Px(t) + p(t)$, dẫn đến việc xác định ma trận đối xứng $P$ thông qua phương trình đại số ma trận Riccati (Algebraic Riccati Equation - ARE): $$PA - PBR^{-1}B^T P + A^T P + Q = 0$$ và xác định thành phần dẫn tiếp $p(t)$ qua phương trình vi phân ngược thời gian: $$\dot{p}(t) + (G^T B^T - A^T)p(t) + PHf(t) = 0, \quad p(\infty) = 0$$ với $G = R^{-1}B^T P$ là ma trận phản hồi trạng thái.
    • Nguyên lý tách (Separation Principle): Cho phép tách việc giải bài toán điều khiển phản hồi đầu ra tổng quát thành hai bước độc lập: thiết kế bộ điều khiển phản hồi trạng thái (Chương 2) và thiết kế bộ quan sát trạng thái tối ưu Kalman-Bucy (Chương 3).

  2. Phương pháp nghiên cứu thực nghiệm và công cụ tính toán:

    • Phương pháp mô phỏng số (numerical simulation): Sử dụng phần mềm MATLAB và môi trường SIMULINK để lập trình thuật toán, tích phân số hệ phương trình vi phân có trễ và giải phương trình Riccati bằng phương pháp chéo hóa ma trận Hamilton.
    • Phương pháp Monte Carlo: Dùng để mô tả và mô phỏng các tiến trình tải trọng ngẫu nhiên tác động lên hệ thống.
    • Phương pháp phần tử hữu hạn: Sử dụng phần mềm kết cấu chuyên dụng SAP2000 để phân tích dao động riêng và tải trọng sóng tĩnh/động trên mô hình kết cấu giàn khung không gian, trích xuất 4 dạng riêng đầu tiên phục vụ tính toán điều khiển.

  3. Nguồn tư liệu và dữ liệu kiểm chứng:

    • Toàn bộ thông số vật lý của kết cấu, thiết bị điều khiển và tải trọng kích động được thu thập từ các công trình khoa học chuẩn đã công bố: số liệu mô hình ô tô từ Muller vcs (1997); mô hình chuyển động dọc máy bay từ Blakelock (1991) và Nelson (1998); các thông số kết cấu công trình chịu tải sóng và thiết bị AMD từ Soong vcs (1997) và Spencer vcs (1997).

Nội dung chính theo từng chương

Chương 1: TỔNG QUAN

Chương 1 hệ thống hóa các định nghĩa và nguyên lý hoạt động của hệ thống điều khiển tích cực trong cơ học. Hệ thống bao gồm ba khâu cốt lõi: cảm biến đo đáp ứng/kích động, bộ điều khiển xử lý thuật toán và máy kích động sinh lực. Phân biệt ba phương thức điều khiển: điều khiển mở (dẫn tiếp - feedforward), điều khiển đóng (hồi tiếp - feedback), và điều khiển đóng-mở phối hợp (hồi tiếp - dẫn tiếp - FB-FF). Về bản chất vật lý, điều khiển phản hồi làm thay đổi tham số tương đương (độ cứng và độ cản) của kết cấu, trong khi điều khiển dẫn tiếp triệt tiêu trực tiếp tác động của kích động ngoài.

Về thiết bị sinh lực, chương tổng kết bốn dạng máy kích động chính:

  • Kích động thủy lực: Tạo lực qua áp suất dòng chất lỏng đẩy pittông, lực sinh ra đạt tới vài trăm meganewton (MN).
  • Kích động mô tơ điện: Sinh lực qua chuyển động quay kết hợp vít me bi (ball-screw) biến đổi thành chuyển động tịnh tiến, khả năng sinh lực đạt vài trăm kilonewton (kN).
  • Kích động điện từ: Sử dụng tương tác từ trường giữa hai dòng điện, khả năng sinh lực đạt vài kN.
  • Thiết bị dựa trên vật liệu thông minh: Bộ giảm chấn chất lỏng từ biến MR, vật liệu áp điện, hợp kim nhớ hình SMA.

Chương phân tích ba phương thức cơ học để tác động lực điều khiển vào hệ:

  1. Ngoại lực tác động trực tiếp từ bên ngoài (ví dụ: máy kích động tì vào mố cầu cố định để giảm dao động dầm cầu, hoặc lực đẩy phản lực của máy bay).
  2. Lực sinh ra giữa hai bộ phận chuyển động tương đối bên trong hệ (ví dụ: thanh giằng giữa hai nút khung, cơ cấu góc quay của bánh lái máy bay/tàu thủy).
  3. Lực quán tính sinh ra giữa hệ kết cấu và một khối lượng phụ gắn thêm (hệ AMD chuyển động tịnh tiến hoặc quay).

Các ứng dụng thực tế quy mô lớn được phân tích chi tiết:

  • Kyobashi Seiwa Building (Nhật Bản, 1989): Tòa nhà 11 tầng (cao 33m, nặng 400 tấn), lắp đặt hệ AMD gồm 2 khối lượng phụ (4 tấn giảm dao động ngang, 1 tấn giảm dao động xoắn) vận hành bằng máy kích động thủy lực, làm giảm đáp ứng dao động do động đất ($10\text{ cm/s}^2$) và gió ($20\text{ m/s}$) xuống còn 1/3.
  • Yokohama Landmark Tower (Nhật Bản, 1993): Sử dụng AMD cơ cấu con lắc 3 tầng dài 3m (chuyển động tương đương con lắc 9m), dẫn động bằng mô tơ điện qua vít me bi.
  • Applause Tower (Nhật Bản, 2001): Tận dụng trực tiếp sân bay trực thăng nặng 480 tấn trên nóc làm khối lượng phụ AMD, đặt trên hệ đệm cao su nhiều tầng và điều khiển bởi 5 máy kích động thủy lực.
  • Nanjing TV Tower (Trung Quốc, 1998): Tháp truyền hình cao 340m, lắp đặt AMD dạng vành đai nặng 59 tấn trượt trên đệm cao su Teflon bao quanh vách tháp, điều khiển bằng 3 máy kích động thủy lực hành trình 1.2m bố trí lệch nhau các góc $120^\circ$.
  • Shinjuku Park Tower (Nhật Bản, 1994): Khối lượng phụ nặng 110 tấn chuyển động trên hệ con lăn và ray chữ V, điều khiển bởi động cơ điện với cơ cấu bánh răng. Kết quả đo đạc 5 năm cho thấy gia tốc cực đại giảm từ $6.8\text{ cm/s}^2$ xuống $2.2\text{ cm/s}^2$.

Chương khép lại bằng việc xác định 8 vấn đề cốt lõi của lý thuyết điều khiển hiện đại đối với kết cấu lớn (tính bền vững, hiện tượng spillover, thời gian trễ, tải ngẫu nhiên, tính phi tuyến, hạn chế cảm biến/máy kích động, kỹ thuật rời rạc hóa và tương tác kết cấu - máy kích động) và xác định rõ mục tiêu nghiên cứu của luận án.

Chương 2: BÀI TOÁN ĐIỀU KHIỂN KHÔNG HẠN CHẾ ĐO

Chương 2 giải quyết bài toán với giả thiết đo được đầy đủ vectơ trạng thái $x(t)$ của hệ. Lời giải giải tích tối ưu thực sự từ hàm Hamilton dẫn đến luật điều khiển: $$u(t) = u_{FB}(t) + u_{FF}(t) = -Gx(t) - R^{-1}B^T p(t)$$ Trong đó, thành phần phản hồi $u_{FB}(t)$ được xác định qua ma trận Riccati $P$, còn thành phần dẫn tiếp $u_{FF}(t)$ phụ thuộc vào nghiệm của phương trình vi phân ngược thời gian của $p(t)$. Do đòi hỏi phải biết trước toàn bộ tiến trình tương lai của kích động ngoài ngẫu nhiên $f(t)$, lời giải tối ưu thực sự không thể thực thi trong thực tế. Thuật toán LQR kinh điển buộc phải giả thiết $f(t) = 0$, loại bỏ hoàn toàn $p(t)$ và chỉ sử dụng luật phản hồi $u(t) = -Gx(t)$. Sai số lực điều khiển giữa LQR và tối ưu thực sự là $E_{LQR} = |R^{-1}B^T p(t)|$.

Tính không tối ưu của LQR được minh chứng định lượng trong Ví dụ 2.1 (mô hình ô tô hai bậc tự do chịu mấp mô mặt đường dạng sóng điều hòa với vận tốc $v = 50\text{ km/h}$, biên độ $Y = 50\text{ cm}$, các thông số: $m_1 = 1200\text{ kg}$, $m_2 = 80\text{ kg}$, $c_1 = 4800\text{ Ns/m}$, $k_1 = 3\cdot 10^4\text{ N/m}$, $k_2 = 32\cdot 10^4\text{ N/m}$, trọng số $r = 5\cdot 10^{-5}$). Số liệu trích xuất từ Bảng 2.1 của luận án được tổng hợp dưới đây:

Nửa bước sóng $L$ Phương án đ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.052 2.21 0.00 0.42
Thuật toán LQR 0.024 0.94 0.82 0.28
Tối ưu thực sự 0.021 0.81 0.75 0.22
$L = 6.94\text{ m}$ Không điều khiển 0.038 1.84 0.00 0.18
Thuật toán LQR 0.019 0.76 0.65 0.14
Tối ưu thực sự 0.016 0.63 0.58 0.10
$L = 4.63\text{ m}$ Không điều khiển 0.081 3.65 0.00 0.95
Thuật toán LQR 0.089 4.12 1.45 1.25
Tối ưu thực sự 0.045 2.15 1.10 0.41

Kết quả Bảng 2.1 chỉ ra hiện tượng: ở bước sóng $L = 4.63\text{ m}$, thuật toán LQR làm gia tăng chỉ tiêu $J$ lên $1.25\text{ Js}$, lớn hơn cả trạng thái không điều khiển ($0.95\text{ Js}$). Điều này khẳng định LQR không đảm bảo tính tối ưu khi có kích động ngoài.

Về độ trễ, phương trình chuyển động khi xuất hiện thời gian trễ của lực $\tau_u$ là $\dot{x}(t) = Ax(t) + BGx(t - \tau_u) + Hf(t)$. Tính ổn định được khảo sát qua phương trình đặc trưng siêu việt: $$\det(A - BGe^{-\tau_u \lambda} - \lambda I_n) = 0$$ Khai triển chuỗi Taylor chứng minh rằng khi $\tau_u > 0$, việc tăng độ lớn các phần tử ma trận $G$ sẽ khuyếch đại ảnh hưởng của trễ. Trong Ví dụ 2.2 về điều khiển chuyển động dọc máy bay (ma trận $A$ cấp $4 \times 4$, $B$ cấp $4 \times 2$), khi giảm hệ số trọng số $r$ xuống mức $2\cdot 10^{-6}$ để tăng cường lực phản hồi, thời gian trễ đã làm toàn bộ hệ thống mất ổn định.

Để khắc phục, tác giả đề xuất thuật toán hồi tiếp - dẫn tiếp (FB-FF). Khi kích động có tần số thấp, bỏ qua đạo hàm $\dot{p}(t)$ trong phương trình vi phân để thu được xấp xỉ: $$\bar{p}(t) = -(G^T B^T - A^T)^{-1} P H f(t)$$ Luật điều khiển FB-FF tương ứng là: $$u(t) = -Gx(t) - R^{-1}B^T (G^T B^T - A^T)^{-1} P H f(t)$$ Điều kiện để thuật toán FB-FF vượt trội hơn LQR là đại lượng phi thứ nguyên $I_f < 1$: $$I_f = \frac{\max_t |(G^T B^T - A^T)^{-1} P H \dot{f}(t)|}{\max_t |PHf(t)|} < 1$$ Hiệu quả được kiểm chứng trong Ví dụ 2.3 trên khung không gian 95 nút (cao 39m, nặng 500 tấn, trang bị AMD 20 tấn) chịu tải trọng sóng chu kỳ 11s, chiều cao 15m. Mô hình thu gọn 4 dạng riêng đầu tiên cho thấy FB-FF đạt khả năng giảm gia tốc tương đương LQR nhưng giảm đáng kể lực điều khiển và hành trình dịch chuyển của khối lượng AMD.

Nhằm triển khai FB-FF khi không đo được $f(t)$, tác giả đề xuất thuật toán nhận dạng kích động ngoài trực tuyến từ phương trình trạng thái. Để triệt tiêu vòng lặp đại số (algebraic loop) phát sinh do biến $u_{FF}$ xuất hiện đồng thời ở hai vế, tác giả đưa vào bước trễ nhân tạo $\tau$, hình thành công thức truy hồi: $$u_{FF}(t) = -R^{-1}B^T (G^T B^T - A^T)^{-1} P [\dot{x}(t - \tau) - (A - BG)x(t - \tau) - B u_{FF}(t - \tau)]$$ Quá trình tính toán số được thực hiện phân đoạn qua các khoảng thời gian $\Delta T_k = [ (k-1)\tau, k\tau ]$ với $u_{FF}(t) = 0$ tại thời đoạn đầu tiên. Cuối chương, tác giả đề xuất tiêu chuẩn xác định vị trí đặt lực tối ưu dựa trên việc cực tiểu hóa sai số giữa lực FB-FF và lực tối ưu thực sự.

Chương 3: BÀI TOÁN ĐIỀU KHIỂN KHÔNG HẠN CHẾ ĐẶT LỰC

Chương 3 xem xét bài toán đối ngẫu: hệ thống có đầy đủ máy kích động trên các bậc tự do nhưng bị hạn chế số lượng đầu đo (chỉ thu nhận được vectơ tín hiệu đo đầu ra $y(t) = Cx(t) + v(t)$ với số chiều nhỏ hơn bậc tự do của hệ).

Chương phân tích luật điều khiển dựa trên bộ quan sát tối ưu Kalman-Bucy. Bộ quan sát tái tạo vectơ trạng thái ước lượng $\hat{x}(t)$ từ số liệu đo một phần. Tương tự như LQR, bộ quan sát Kalman-Bucy kinh điển chỉ đạt tính tối ưu khi kích động ngoài là nhiễu trắng ngẫu nhiên; khi kích động mang tính màu hoặc điều hòa, sai số ước lượng trạng thái gia tăng đáng kể.

Để khắc phục, tác giả phát triển thuật toán nhận dạng kích động ngoài dựa trên các biến đo đầu ra:

  • Thiết lập công thức nhận dạng cho hai trường hợp đo đạc cơ học thực tế: đo chuyển động tương đối và đo chuyển động tuyệt đối của kết cấu.
  • Đề xuất phương pháp xác định vị trí đặt đầu đo (cảm biến) tối ưu nhằm cực đại hóa độ chính xác nhận dạng lực kích động và giảm thiểu sai số ước lượng trạng thái của bộ quan sát.
  • Bổ sung thành phần dẫn tiếp được tính toán từ kích động nhận dạng vào luật điều khiển hồi tiếp trạng thái ước lượng.

Kết quả mô phỏng số được thực hiện trên hệ kết cấu dầm thẳng đứng chịu tải trọng sóng và dầm ngang chịu tải sóng dập dềnh. Các đại lượng khảo sát bao gồm:

  • Các giá trị riêng của hệ dầm đứng chịu tải sóng (Bảng 3.1) và hệ dầm ngang (Bảng 3.2).
  • Hiệu quả nhận dạng lực trong trường hợp đo chuyển động tương đối (Bảng 3.3) và đo chuyển động tuyệt đối (Bảng 3.4).
  • Đáp ứng RMS chuyển dịch giữa dầm (Bảng 3.5).
  • Khảo sát sự hội tụ của chuyển dịch đỉnh dầm qua các hệ số trọng số $\mu = 10^{-3}, 10^{-4}, 10^{-5}$ và ảnh hưởng khi có thời gian trễ $0.1\text{ s}$ (Hình 3.3 đến Hình 3.8).
  • Đánh giá phổ phần thực và phần ảo của kích động thực so với kích động xấp xỉ khi bố trí cảm biến tại các vị trí khác nhau: $l/4$, $l/2$, và $3l/4$ (Hình 3.11 đến Hình 3.13).
  • Phân tích chuyển dịch giữa dầm chịu ảnh hưởng của cường độ kích động $I_v$ ($10^{-10}\text{ m}^2\text{/s}$ và $10^{-12}\text{ m}^2\text{/s}$) và tham số thời gian trễ $\tau_v$ ($5\cdot 10^{-4}\text{ s}$ và $25\cdot 10^{-4}\text{ s}$) (Hình 3.15 đến Hình 3.18).

Chương 4: BÀI TOÁN ĐIỀU KHIỂN TỐI ƯU PHẢN HỒI ĐẦU RA

Chương 4 kết hợp các kết quả của Chương 2 và Chương 3 để giải quyết bài toán kỹ thuật tổng quát nhất: kết cấu đồng thời bị hạn chế số lượng cảm biến đo và hạn chế số lượng vị trí đặt máy kích động.

Dựa trên nguyên lý tách (Separation Principle), tác giả phân tích thuật toán điều khiển tối ưu phản hồi đầu ra kinh điển LQG (kết hợp bộ điều khiển trạng thái LQR với bộ lọc tối ưu Kalman-Bucy). Từ các hạn chế cố hữu của LQG trước kích động ngẫu nhiên khác ồn trắng, tác giả đề xuất hai phiên bản điều khiển nâng cao:

  1. Phiên bản điều khiển nâng cao 1: Tích hợp thành phần dẫn tiếp nhận dạng kích động ngoài vào mạch điều khiển phản hồi đầu ra, sử dụng trạng thái ước lượng từ bộ quan sát kết hợp với bước trễ khử vòng lặp đại số.
  2. Phiên bản điều khiển nâng cao 2: Cải tiến cấu trúc ước lượng và nhận dạng đồng thời, tối ưu hóa ma trận khuếch đại phản hồi đầu ra nhằm triệt tiêu dao động hiệu quả hơn dưới tác động của nhiễu đo và độ trễ.

Chương tiến hành đánh giá chi tiết định lượng qua mô phỏng số:

  • Khảo sát ảnh hưởng của ma trận trọng số đến đáp ứng kết cấu (Bảng 4.1).
  • Đánh giá ảnh hưởng của vị trí cảm biến và mức độ nhiễu đo: mô phỏng khi đặt cảm biến tại giữa dầm $d_m = l/2$ với nhiễu đo nhỏ (Bảng 4.2) và nhiễu lớn (Bảng 4.3); khi đặt cảm biến tại vị trí một phần tư dầm $d_m = l/4$ với nhiễu nhỏ (Bảng 4.4) và nhiễu lớn (Bảng 4.5).
  • So sánh hiệu quả kiểm soát dao động qua ba trạng thái: đáp ứng kết cấu khi chưa điều khiển (Bảng 4.6), đáp ứng khi điều khiển bằng thuật toán LQG kinh điển (Bảng 4.7), và đáp ứng khi điều khiển bằng thuật toán nâng cao phiên bản 2 (Bảng 4.8).
  • Các ví dụ kết cấu minh họa gồm có: hệ dầm đứng chịu tải gió được điều khiển bởi AMD, hệ thống điều khiển góc phương vị và hiệu điện thế mô tơ của ăng-ten parabol quan sát bầu trời, kết cấu dầm ngang, và kết cấu công trình 8 tầng chịu tải ngẫu nhiên lắp đặt thiết bị AMD.

Kết quả và những đóng góp mới

  1. Đóng góp mới về mặt lý luận khoa học:

    • Chứng minh về mặt giải tích và định lượng giới hạn của thuật toán phản hồi kinh điển LQR/LQG khi kích động ngoài khác 0; làm rõ cơ chế khuyếch đại thời gian trễ khi gia tăng độ lớn các hệ số trong ma trận phản hồi thông qua phương trình đặc trưng siêu việt.
    • Đề xuất thuật toán điều khiển hồi tiếp - dẫn tiếp (FB-FF) xấp xỉ thành phần dẫn tiếp cho các kích động ngẫu nhiên có tần số thấp, đồng thời thiết lập điều kiện giải tích tường minh dạng phi thứ nguyên ($I_f < 1$) để đảm bảo thuật toán đề xuất đạt độ chính xác cao hơn LQR kinh điển.
    • Phát triển thuật toán nhận dạng trực tuyến kích động ngoài từ đáp ứng đo được của kết cấu; đề xuất giải pháp đưa bước trễ nhân tạo $\tau$ để triệt tiêu vòng lặp đại số trong tính toán số, cho phép tính toán quy nạp theo từng thời đoạn thực thi.
    • Xây dựng tiêu chuẩn xác định vị trí đặt máy kích động tối ưu và vị trí đặt đầu đo tối ưu dựa trên độ đo giảm thiểu sai số so với nghiệm tối ưu thực sự.
    • Thiết lập cấu trúc hai phiên bản điều khiển nâng cao cho bài toán điều khiển phản hồi đầu ra tổng quát (vừa hạn chế đo vừa hạn chế đặt lực) trên cơ sở nguyên lý tách.

  2. Ý nghĩa thực tiễn và giải pháp kiến nghị:

    • Cung cấp công cụ tính toán và thuật toán khả thi để nâng cao hiệu quả giảm chấn cho các công trình lớn (nhà cao tầng, tháp truyền hình, giàn khoan biển) bằng thiết bị AMD mà không đòi hỏi lắp đặt quá nhiều cảm biến hay đo trực tiếp tải trọng môi trường.
    • Kiến nghị việc lựa chọn bước trễ tính toán quy nạp $\tau$ phải đồng bộ với năng lực tính toán của vi xử lý điều khiển và thời gian phản ứng cơ học của máy kích động thủy lực hoặc động cơ điện.
    • Kiến nghị ứng dụng thuật toán FB-FF cho các công trình xây dựng vùng nước nông hoặc công trình ven biển chịu tải trọng sóng có chu kỳ dài (tần số thấp khoảng $0.1\text{ Hz}$), nơi thuật toán thể hiện ưu thế giảm mạnh năng lượng tiêu hao và chuyển vị của khối lượng phụ AMD.

Hạn chế và hướng nghiên cứu tiếp

Các hạn chế được văn bản đề cập:

  1. Thuật toán điều khiển FB-FF đề xuất chỉ chứng minh được tính ưu việt hơn LQR khi kích động ngoài có tần số biến thiên thấp thỏa mãn điều kiện $I_f < 1$. Với các kích động chứa dải tần số cao, hiệu quả xấp xỉ của thành phần dẫn tiếp suy giảm.
  2. Mô hình toán học áp dụng trong luận án hoàn toàn xây dựng trên cơ sở hệ tuyến tính hoặc các hệ phi tuyến đã được tuyến tính hóa quanh vị trí cân bằng tĩnh; chưa bao quát các trạng thái phi tuyến mạnh của vật liệu thông minh (như chất lỏng từ biến MR hoặc hợp kim nhớ hình SMA) hoặc các liên kết kết cấu phi tuyến.
  3. Nghiên cứu khảo sát thời gian trễ chủ yếu dưới dạng tham số trễ cố định $\tau$; việc tìm nghiệm giải tích tổng quát cho phương trình đặc trưng siêu việt chứa trễ vẫn còn là bài toán phức tạp chưa giải quyết trọn vẹn trong khuôn khổ luận án.
  4. Các kết quả đánh giá trong luận án hoàn toàn dựa trên phương pháp mô phỏng số thông qua phần mềm MATLAB/SIMULINK và SAP2000; chưa tiến hành thử nghiệm vật lý trên mô hình thực nghiệm rung chấn ngoài hiện trường.

Hướng nghiên cứu tiếp:

  • Mở rộng thuật toán hồi tiếp - dẫn tiếp kết hợp nhận dạng lực cho các hệ kết cấu phi tuyến hình học và phi tuyến vật liệu.
  • Nghiên cứu cơ chế thích nghi của bước trễ $\tau$ theo sự thay đổi động của tần số kích động và nghiên cứu thử nghiệm thực nghiệm trên mô hình bàn rung (shaking table).

Giá trị tham khảo

Luận án có giá trị tham khảo chuyên môn cho các đối tượng sau:

  • Nghiên cứu sinh, học viên cao học và giảng viên các chuyên ngành Cơ học kỹ thuật, Cơ học công trình, Kỹ thuật xây dựng công trình ngầm/biển, Kỹ thuật điều khiển và tự động hóa: Có thể khai thác làm tài liệu tham khảo chuyên sâu về phương pháp không gian trạng thái, kỹ thuật tối ưu hóa LQR/LQG, lý thuyết bộ lọc Kalman-Bucy và nguyên lý tách trong bài toán dao động công trình.
  • Kỹ sư thiết kế kết cấu công trình chịu tải trọng động lực học môi trường (gió bão, động đất, sóng biển), kỹ sư hàng không và cơ điện tử: Ứng dụng các quy tắc thiết kế hệ thống giảm chấn tích cực AMD, xác định vị trí tối ưu của cảm biến và máy kích động để giảm thiểu chi phí phần cứng.
  • Phần đáng tham khảo nhất: Nội dung Mục 2.4 và 2.5 (kỹ thuật xấp xỉ thành phần dẫn tiếp và thuật toán khử vòng lặp đại số trong nhận dạng lực quy nạp); các bảng dữ liệu mô phỏng đối chuẩn (Bảng 2.1 cho mô hình ô tô; Bảng 4.6 - 4.8 cho hệ phản hồi đầu ra) làm cơ sở benchmark cho các thuật toán điều khiển mới.

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

1. Tại sao thuật toán điều khiển kinh điển LQR không đạt được tính tối ưu thực sự khi kết cấu chịu kích động từ môi trường? Trong điều kiện tối ưu suy ra từ lý thuyết biến phân Lagrange, luật điều khiển tối ưu thực sự bao gồm hai thành phần: thành phần hồi tiếp $u_{FB}(t) = -Gx(t)$ và thành phần dẫn tiếp $u_{FF}(t) = -R^{-1}B^T p(t)$. Vectơ $p(t)$ là nghiệm của phương trình vi phân ngược thời gian phụ thuộc trực tiếp vào kích động ngoài $f(t)$ với điều kiện cuối $p(\infty) = 0$. Thuật toán LQR kinh điển giả thiết $f(t) = 0$ nên triệt tiêu thành phần $p(t)$. Do đó, khi kích động ngoài thực tế khác 0, luật LQR tạo ra sai số $E_{LQR} = |R^{-1}B^T p(t)|$. Trong Ví dụ 2.1, khi mặt đường có bước sóng $L = 4.63\text{ m}$, thuật toán LQR làm chỉ tiêu tích phân $J$ tăng lên $1.25\text{ Js}$, lớn hơn cả mức $0.95\text{ Js}$ của trường hợp hoàn toàn không điều khiển.

2. Vòng lặp đại số (algebraic loop) trong thuật toán nhận dạng kích động ngoài xuất hiện như thế nào và được xử lý bằng cách nào? Khi biểu diễn lực kích động nhận dạng từ phương trình trạng thái $Hf(t) = \dot{x}(t) - (A - BG)x(t) - Bu_{FF}(t)$ và thay thế vào công thức thành phần dẫn tiếp $u_{FF}(t) = -R^{-1}B^T(G^T B^T - A^T)^{-1}PHf(t)$, đại lượng $u_{FF}(t)$ xuất hiện đồng thời ở cả hai vế của phương trình đại số. Hiện tượng này tạo thành một vòng lặp đại số khiến quá trình tính toán số bị tắc nghẽn hoặc không xác định. Tác giả đã giải quyết bằng cách áp dụng một khoảng trễ nhân tạo nhỏ $\tau$ vào các biến đo và lực điều khiển ở vế phải, đưa công thức về dạng quy nạp $u_{FF}(t) = -R^{-1}B^T(G^T B^T - A^T)^{-1}P[\dot{x}(t-\tau) - (A-BG)x(t-\tau) - Bu_{FF}(t-\tau)]$, cho phép thực hiện việc tính toán phân đoạn theo từng bước thời gian $\Delta T_k$.

3. Điều kiện định lượng nào cho phép khẳng định thuật toán điều khiển FB-FF đề xuất vượt trội hơn thuật toán LQR? Thuật toán FB-FF đạt được chuẩn sai số lực $E_{FBFF}$ nhỏ hơn chuẩn sai số của LQR khi chuẩn của sai số ước lượng $e_p(t) = p(t) - \bar{p}(t)$ nhỏ hơn chuẩn của bản thân vectơ $p(t)$. Thông qua việc đổi chiều thời gian và đánh giá chặn trên các tích phân chuẩn ma trận, điều kiện cần thỏa mãn là đại lượng phi thứ nguyên $I_f$ phải nhỏ hơn 1: $$I_f = \frac{\max_t |(G^T B^T - A^T)^{-1} P H \dot{f}(t)|}{\max_t |PHf(t)|} < 1$$ Điều kiện này xác nhận thuật toán FB-FF đạt hiệu quả cao khi kích động ngoài có tốc độ biến thiên chậm (tần số thấp như tải trọng sóng), khiến chuẩn của đạo hàm $\dot{f}(t)$ đủ nhỏ so với chuẩn của $f(t)$.

4. Thời gian trễ của máy kích động ảnh hưởng ra sao đến tính ổn định của hệ điều khiển phản hồi LQR? Khi tồn tại thời gian trễ $\tau_u$ của lực điều khiển, phương trình chuyển động của hệ trở thành $\dot{x}(t) = Ax(t) + BGx(t - \tau_u) + Hf(t)$. Tính ổn định của hệ nghiệm phụ thuộc vào các nghiệm $\lambda$ của phương trình đặc trưng siêu việt $\det(A - BGe^{-\tau_u \lambda} - \lambda I_n) = 0$. Khi không có trễ ($\tau_u = 0$), ma trận $(A - BG)$ luôn có phần thực của các giá trị riêng âm, đảm bảo hệ ổn định tiệm cận. Tuy nhiên, khi $\tau_u > 0$, nếu tăng các hệ số của ma trận phản hồi $G$ (ví dụ bằng cách giảm trọng số $r$ như trong bài toán máy bay ở Ví dụ 2.2 xuống $r = 2\cdot 10^{-6}$), chuỗi số hạng trễ bị khuyếch đại, làm dịch chuyển nghiệm $\lambda$ sang nửa phải mặt phẳng phức, dẫn tới mất ổn định hệ thống.

5. Nguyên lý tách (Separation Principle) được vận dụng như thế nào trong cấu trúc của Chương 4? Nguyên lý tách trong lý thuyết điều khiển ngẫu nhiên khẳng định: đối với hệ tuyến tính chịu chỉ tiêu toàn phương và nhiễu ngẫu nhiên, bài toán điều khiển phản hồi đầu ra tổng quát có thể được tách biệt hoàn toàn thành hai bài toán độc lập: (1) thiết kế bộ điều khiển phản hồi trạng thái tối ưu giả định đo được đầy đủ trạng thái (giải quyết trong Chương 2), và (2) thiết kế bộ quan sát tối ưu tái tạo trạng thái từ các tín hiệu đo đầu ra hạn chế (giải quyết trong Chương 3). Trong Chương 4, tác giả sử dụng trạng thái ước lượng từ bộ quan sát thay thế cho trạng thái thực trong luật điều khiển trạng thái, kết hợp thuật toán nhận dạng lực trực tuyến để hình thành hai phiên bản điều khiển nâng cao cho bài toán hạn chế cả đo lẫn đặt lực.

Kết luận

Luận án "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" đã giải quyết một vấn đề kỹ thuật quan trọng trong kiểm soát dao động kết cấu công trình chịu tải ngẫu nhiên khi số lượng cảm biến và máy kích động bị hạn chế. Bằng việc phân tích giải tích hạn chế của các phương pháp LQR/LQG kinh điển, tác giả đã đề xuất thành công thuật toán điều khiển hồi tiếp - dẫn tiếp (FB-FF) kết hợp với thuật toán nhận dạng kích động ngoài quy nạp có trễ để triệt tiêu vòng lặp đại số. Các giải pháp bổ sung về xác định vị trí tối ưu của cảm biến và máy kích động cùng hai phiên bản điều khiển phản hồi đầu ra nâng cao đã được kiểm chứng nhất quán qua các mô hình mô phỏng số chi tiết (ô tô, máy bay, khung không gian, dầm chịu tải sóng và công trình lắp đặt AMD). Kết quả nghiên cứu cung cấp cơ sở khoa học và công cụ tính toán tin cậy cho việc thiết kế, vận hành các hệ thống giảm chấn tích cực trên các công trình quy mô lớn.