Tổng quan về luận án
Kiểm soát và giảm thiểu dao động có hại trong các công trình xây dựng, cầu nhịp lớn, tháp truyền hình và kết cấu cơ khí chính xác là một trong những bài toán cốt lõi của cơ học công trình và kỹ thuật điều khiển hiện đại. Các tác động ngẫu nhiên từ tự nhiên như động đất mạnh, bão gió hoặc tải trọng xung kích động lực học không chỉ gây mỏi vật liệu, suy giảm tuổi thọ mà còn đe dọa trực tiếp đến an toàn sinh mạng và tính toàn vẹn của kết cấu hạ tầng. Trước đây, phương pháp truyền thống chủ yếu dựa trên việc gia tăng tiết diện và độ cứng kết cấu; tuy nhiên, giải pháp này gây tốn kém vật liệu và làm gia tăng đáng kể tải trọng bản thân của công trình. Trong bối cảnh công nghệ phát triển, việc chuyển dịch sang ứng dụng các thiết bị tiêu tán năng lượng và hệ thống điều khiển thông minh đã mang lại hiệu quả vượt trội. Như luận án đã nêu rõ: "Ước tính, sử dụng thiết bị tiêu tán năng lượng có thể chỉ chiếm 25% chi phí so với việc gia cố kết cấu cho các bộ phận thép và bê tông."
Mặc dù các hệ thống điều khiển chủ động như hệ thống giằng chủ động (Active Bracing System - ABS) hay thiết bị giảm chấn khối lượng chủ động (Active Mass Damper - AMD) đã chứng minh được hiệu năng thực tế tại nhiều công trình biểu tượng trên thế giới, bài toán thiết kế thuật toán điều khiển vẫn đối mặt với những rào cản lớn. Các phương pháp điều khiển kinh điển như bộ điều chỉnh bậc hai tuyến tính (Linear Quadratic Regulator - LQR) hay bộ điều khiển tuyến tính Gauss bậc hai (Linear Quadratic Gaussian - LQG) đòi hỏi mô hình toán học chính xác tuyệt đối của kết cấu – một yêu cầu bất khả thi đối với các hệ nhiều bậc tự do phi tuyến, có tham số bất định và chịu kích động ngẫu nhiên. Để khắc phục, lý thuyết tập mờ (Fuzzy Logic Control - FC) do Zadeh (1965) khởi xướng đã được áp dụng rộng rãi nhờ khả năng xử lý thông tin định tính thông qua bảng nhớ liên kết mờ (Fuzzy Associative Memory - FAM). Tuy nhiên, FC truyền thống bộc lộ những hạn chế cố hữu: quá trình mờ hóa, suy luận hợp thành Mamdani và giải mờ trọng tâm (Centroid defuzzification) đòi hỏi khối lượng tính toán lớn, số lượng tham số cần tối ưu hóa quá nhiều (lên tới $3n$ tham số đối với các hàm thuộc tam giác), đồng thời tiềm ẩn nguy cơ đảo lộn thứ tự ngữ nghĩa tự nhiên của các biến ngôn ngữ.
Khoảng trống nghiên cứu (research gap) cốt lõi xuất hiện khi lý thuyết Đại số gia tử (Hedge Algebras - HA) do Nguyễn Cát Hồ và các cộng sự phát triển từ năm 1990 mở ra hướng tiếp cận đại số hóa hoàn toàn các biến ngôn ngữ thông qua các ánh xạ ngữ nghĩa định lượng (Semantically Quantifying Mapping - SQM). Mặc dù HA đã bước đầu được ứng dụng vào bộ điều khiển mờ dựa trên đại số gia tử (Hedge Algebras Controller - HAC), các nghiên cứu trước đây chưa giải quyết được bài toán tối ưu hóa đa mục tiêu đối với các tham số độ đo tính mờ của các biến trạng thái, cũng như chưa thiết lập được cơ chế định lượng hóa để tinh chỉnh hệ luật điều khiển nhằm thích ứng với từng mô hình động lực học kết cấu cụ thể.
Luận án tiến sĩ cơ học của NCS Trần Quý Cao với tiêu đề "Ứng dụng lý thuyết mờ và đại số gia tử trong điều khiển dao động kết cấu" (Mã số: 9440109, Trường Đại học Bách khoa Hà Nội, hướng dẫn: PGS. Bùi Hải Lê và TS. Bùi Văn Bình) đã tập trung giải quyết triệt để bài toán này. Luận án đặt ra hai câu hỏi nghiên cứu then chốt:
- RQ1: Làm thế nào để thiết lập một khung tối ưu hóa đa mục tiêu cho các tham số mờ của biến trạng thái (chuyển vị, vận tốc) và biến điều khiển (lực kích động) trong bộ điều khiển HAC nhằm đạt được sự cân bằng tối ưu giữa việc giảm chuyển vị, triệt tiêu gia tốc và tiết kiệm năng lượng điều khiển?
- RQ2: Cơ chế hiệu chỉnh hệ luật định tính của HAC thông qua các hệ số điều chỉnh (tuning coefficients) hoạt động như thế nào để vừa duy trì tính đơn điệu ngữ nghĩa tự nhiên, vừa nâng cao hiệu năng điều khiển cho các kết cấu kỹ thuật phức tạp?
Tương ứng với hai câu hỏi trên là hai giả thuyết khoa học được kiểm chứng:
- H1: Việc tối ưu hóa trực tiếp độ đo tính mờ $f_m(c^-), f_m(c^+)$ và độ đo gia tử $\mu(h^-), \mu(h^+)$ trên không gian nghiệm Pareto sẽ tạo ra bộ điều khiển opHAC vượt trội hơn FC truyền thống cả về khả năng giảm chấn lẫn tốc độ tính toán CPU time.
- H2: Tích hợp hệ số điều chỉnh vào mặt ngữ nghĩa định lượng (Semantically Quantifying Surface - SQS) sẽ tối ưu hóa đáp ứng của hệ mà không cần tái cấu trúc phức tạp bảng luật định tính ban đầu.
Phạm vi nghiên cứu của luận án bao quát từ các mô hình lý thuyết cơ bản (hệ 1 bậc tự do, dầm công-xôn chịu tải xung, hệ con lắc ngược) đến các kết cấu phức tạp nhiều bậc tự do (mô hình khung nhà 3 tầng, 5 tầng chịu gia tốc nền địa chấn và mô hình khung giàn không gian 3D mô phỏng trên phần mềm phần tử hữu hạn ANSYS). Dữ liệu kích động địa chấn được thu thập và xử lý từ ba trận động đất kinh điển có phổ tần số và biên độ đặc trưng: El Centro (1940), Northridge (1994) và Imperial Valley (1979).
Literature Review và Positioning
Lịch sử phát triển của kỹ thuật điều khiển dao động công trình trải qua bốn thế hệ tiếp cận chính: điều khiển bị động (passive control), điều khiển bán chủ động (semi-active control), điều khiển chủ động (active control) và điều khiển lai (hybrid control). Trong lĩnh vực điều khiển chủ động, các công trình tiên phong của Soong (1990) và Housner et al. (1997) đã đặt nền móng cho việc ứng dụng các thuật toán điều khiển hiện đại vào công trình xây dựng. Các phương pháp điều khiển tối ưu như LQR, LQG dựa trên quy hoạch động Bellman (1957) và nguyên lý cực đại Pontryagin (1962), cùng với các bộ lọc Kalman-Bucy (1961), đã đạt được nhiều thành công lý thuyết. Tuy nhiên, tính phụ thuộc quá lớn vào ma trận khối lượng $[M]$, ma trận cản $[C]$ và ma trận độ cứng $[K]$ chính xác khiến các giải pháp này trở nên dễ tổn thương trước các hiện tượng phi tuyến, suy giảm độ cứng vật liệu hoặc hiện tượng bất định tham số (uncertain structures) và trễ điều khiển (time delay) như đã được phân tích bởi Spencer et al. (1998) và Preumont (2002).
Để vượt qua rào cản mô hình toán, trường phái điều khiển thông minh (intelligent control) dựa trên tính toán mềm (soft computing) đã phát triển mạnh mẽ. Dựa trên lý thuyết tập mờ của Zadeh (1965), các công trình của Chang & Soong (1980s), Battaini et al. (1998), Ahlawat & Ramaswamy (2002), và Guclu (2006) đã áp dụng thành công bộ điều khiển mờ (FC) vào hệ thống ABS và AMD. Sự tranh luận học thuật sâu sắc nổ ra giữa hai trường phái:
- Trường phái mô hình giải tích chính xác: Cho rằng điều khiển mờ thiếu tính chặt chẽ về chứng minh ổn định Lyapunov và phụ thuộc quá nhiều vào tri thức chủ quan của chuyên gia khi lập bảng FAM.
- Trường phái điều khiển mờ và tính toán mềm: Lập luận rằng trong điều kiện địa chấn ngẫu nhiên và cấu trúc công trình thực tế có bậc tự do vô hạn, việc dựa vào mô hình toán giải tích là không thực tế; hệ luật "NẾU - THÌ" của logic mờ mang lại tính bền vững (robustness) vượt trội và phản ứng linh hoạt trước các kích động không dự báo trước.
Nhằm tối ưu hóa bộ điều khiển mờ, hàng loạt nghiên cứu quốc tế đã kết hợp thuật toán di truyền (Genetic Algorithm - GA) và tối ưu hóa bầy đàn (Particle Swarm Optimization - PSO) để tinh chỉnh tọa độ đỉnh của các hàm thuộc (membership functions) dạng tam giác, hình thang hoặc Gauss (như nghiên cứu của Marinaki et al., 2011; Uz & Khuc, 2013; Bathaei et al., 2018). Tuy nhiên, phương pháp này làm bùng nổ số lượng biến thiết kế tối ưu ($3n$ biến đối với $n$ hàm thuộc tam giác), đồng thời quá trình biến đổi ngẫu nhiên của GA dễ dẫn đến hiện tượng giao cắt phi logic, phá vỡ thứ tự ngữ nghĩa vốn có của các biến ngôn ngữ (ví dụ: nhãn "Dương" bị dịch chuyển sang phía bên trái của nhãn "Không").
Lý thuyết Đại số gia tử (HA) do Nguyễn Cát Hồ và Wechler (1990, 1992), Nguyễn Cát Hồ et al. (2008) phát minh đã giải quyết căn bản nghịch lý này. HA chứng minh rằng tập từ của một biến ngôn ngữ bản thân nó đã sở hữu một cấu trúc đại số có thứ tự bộ phận tự nhiên, cho phép xác định trực tiếp giá trị của từ thông qua hàm ánh xạ ngữ nghĩa định lượng (SQM) mà không cần xây dựng các hàm thuộc mờ nhân tạo và không cần các bước giải mờ phức tạp. Các nghiên cứu của Bùi Hải Lê et al. (2013, 2016) và Bùi Văn Bình et al. (2018) đã đưa HAC vào điều khiển dầm thông minh và khung nhà chịu động đất. Tuy nhiên, một khoảng trống lý thuyết quan trọng vẫn tồn tại: làm thế nào để tối ưu hóa không gian tham số đại số gia tử theo các tiêu chuẩn đa mục tiêu của kỹ thuật xây dựng, và làm thế nào để tinh chỉnh bề mặt điều khiển định tính thích nghi với các kết cấu có động lực học phức tạp.
Về mặt ứng dụng thực tiễn, luận án định vị nghiên cứu của mình thông qua việc đối chiếu với các công trình quy mô quốc tế đã lắp đặt hệ thống AMD thực tế:
- Tòa nhà Kyobashi Seiwa (Tokyo, Nhật Bản): Công trình thương mại đầu tiên trên thế giới trang bị hệ thống AMD chủ động, với kết quả thực tế ghi nhận: "giảm đáp ứng của kết cấu khi chịu tải trọng động đất (gia tốc 10 cm/s2) và tải trọng gió (vận tốc 20 m/s) xuống còn khoảng bằng 1/3 giá trị khi chưa được điều khiển."
- Tháp Yokohama Landmark (Yokohama, Nhật Bản): Cao 296m, sử dụng hệ thống AMD dạng con lắc 3 tầng treo bằng cáp nhằm gia tăng khả năng sinh công của lực điều khiển trong không gian hẹp.
- Tháp truyền hình Nanjing (Nam Kinh, Trung Quốc): Cao 340m, nơi sự hạn chế về không gian lắp đặt bộ giảm chấn thụ động đã dẫn đến việc ứng dụng thành công hệ thống AMD dạng đai (ring-shaped AMD) kết hợp máy kích động thủy lực công suất lớn.
Nghiên cứu của luận án tiến xa hơn các công trình trên bằng cách cung cấp một công cụ toán học tối ưu hóa giải thuật điều khiển (opHAC), giúp nâng cao hiệu suất giảm chấn của các cơ cấu AMD và ABS này trong khi giảm thiểu năng lượng tiêu thụ của máy kích động.
Đó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à làm phong phú thêm lý thuyết Đại số gia tử của Nguyễn Cát Hồ và lý thuyết điều khiển kết cấu thông minh thông qua ba đóng góp lý thuyết căn bản:
Thứ nhất, luận án đã thiết lập mô hình toán học giải tích tường minh cho bộ điều khiển HAC trong cơ học kết cấu. Khác với bộ điều khiển mờ truyền thống vốn được xem như một "hộp đen" với các thuật toán mờ hóa và giải mờ rời rạc, luận án chỉ ra rằng: "mô hình toán học của HAC để xác định các biến điều khiển từ các biến trạng thái có thể được biểu diễn dưới dạng công thức tường minh." Điều này cho phép chuyển đổi toàn bộ quá trình suy luận mờ thành các phép toán đại số trên mặt cong ngữ nghĩa định lượng (SQS), tạo điều kiện phân tích giải tích tính ổn định và đáp ứng động lực học của hệ kín.
Thứ hai, luận án đề xuất phương pháp tối ưu hóa tham số mờ của biến ngôn ngữ dựa trên cấu trúc đại số gia tử. Trong cấu trúc đại số $AX = (X, G, H, \le)$, với tập phần tử sinh $G = {c^-, c^+}$ (ví dụ: $c^- = \text{Âm}, c^+ = \text{Dương}$) và tập các gia tử $H = H^- \cup H^+$ (ví dụ: $H^- = {\text{Hơi - Little}}, H^+ = {\text{Rất - Very}}$), các tham số mờ bao gồm độ đo tính mờ của phần tử sinh $f_m(c^-), f_m(c^+)$ và độ đo tính mờ của các gia tử $\mu(h^-), \mu(h^+)$ thỏa mãn các tiên đề đại số:
$$f_m(c^-) + f_m(c^+) = 1$$
$$\mu(h^-) + \mu(h^+) = 1$$
Bằng cách chọn các tham số này làm biến thiết kế tối ưu trong bài toán đa mục tiêu, số lượng biến thiết kế giảm từ $3n$ (trong FC) xuống chỉ còn 2 đến 4 biến độc lập cho mỗi biến ngôn ngữ, loại bỏ hoàn toàn hiện tượng bùng nổ chiều số học và triệt tiêu nguy cơ vi phạm trật tự ngữ nghĩa.
Thứ ba, luận án giới thiệu khái niệm và mô hình "Hệ số điều chỉnh" (tuning coefficient) của luật điều khiển định tính. Thay vì thay đổi cấu trúc bảng luật FAM gây mất tính liên tục, luận án gắn mỗi luật $R_k$ với một hệ số điều chỉnh $w_k$, từ đó làm biến dạng cục bộ bề mặt cong ngữ nghĩa định lượng SQS. Đóng góp này tạo ra một bước chuyển dịch phương pháp luận (paradigm shift): chuyển từ việc thiết kế hệ luật mờ "tĩnh" dựa thuần túy vào trực giác chuyên gia sang hệ luật "động" được tối ưu hóa số thích nghi chính xác với tần số riêng và các dạng dao động (mode shapes) của kết cấu cơ học.
Khung phân tích độc đáo
Khung phân tích của luận án là sự tích hợp đa ngành giữa ba trụ cột lý thuyết: Cơ học kết cấu công trình, Lý thuyết đại số gia tử (Hedge Algebras) và Lý thuyết tối ưu hóa đa mục tiêu Pareto (Multi-Objective Optimization).
┌────────────────────────────────────────┐
│ KÍCH ĐỘNG ĐỊA CHẤN / NGOẠI TẢI │
│ (El Centro, Northridge, Imperial) │
└──────────────────┬─────────────────────┘
│
▼
┌────────────────────────────────────────┐
│ MÔ HÌNH ĐỘNG LỰC HỌC │
│ [M]{ẍ} + [C]{ẋ} + [K]{x} = {Fe} + {Fu} │
└──────────────────┬─────────────────────┘
│
Tín hiệu đo trạng thái: xi(t), ẋi(t)
│
▼
┌────────────────────────────────────────┐
│ BỘ ĐIỀU KHIỂN TỐI ƯU opHAC / tHAC │
│ - Ánh xạ ngữ nghĩa định lượng (SQM) │
│ - Tối ưu tham số mờ: fm(c), μ(h) │
│ - Hệ số điều chỉnh luật mờ (wk) │
│ - Mặt cong định lượng ngữ nghĩa (SQS) │
└──────────────────┬─────────────────────┘
│
Lực điều khiển chủ động: ui(t)
│
▼
┌────────────────────────────────────────┐
│ CƠ CẤU CHẤP HÀNH (AMD / ABS) │
│ Triệt tiêu dao động kết cấu thực tế │
└────────────────────────────────────────┘
Quy trình phân tích được thiết lập thông qua các bước ánh xạ toán học chặt chẽ:
- Chuẩn hóa biến trạng thái: Biến chuyển vị $x_i(t) \in [-a_0, a_0]$ và vận tốc $\dot{x}_i(t) \in [-b_0, b_0]$ được chuẩn hóa về đoạn $[0, 1]$ qua ánh xạ $g_x$.
- Ánh xạ ngữ nghĩa định lượng (SQM): Giá trị số được ánh xạ trực tiếp vào không gian ngữ nghĩa đại số gia tử thông qua hàm định lượng $v(x)$.
- Tổng hợp luật điều khiển trên mặt SQS: Giá trị điều khiển ngôn ngữ được xác định dựa trên ma trận liên kết ngữ nghĩa (Semantic Associative Memory - SAM) kết hợp với các hệ số điều chỉnh $w_k$.
- Giải chuẩn hóa và sinh lực điều khiển: Giá trị đầu ra định lượng $v(u) \in [0, 1]$ được giải chuẩn hóa qua ánh xạ $g_u^{-1}$ để xác định giá trị lực điều khiển thực $u_i(t) \in [-u_{\max}, u_{\max}]$ tác động lên máy kích động.
Điều kiện biên (boundary conditions) của khung phân tích được xác định rõ: hệ kết cấu làm việc trong miền đàn hồi tuyến tính, máy kích động thủy lực có giới hạn lực bão hòa cực đại $u_{\max}$, và thời gian trễ tín hiệu điều khiển được giả định nằm trong giới hạn cho phép của các hệ thống máy tính thời gian thực.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án áp dụng thế giới quan thực chứng (positivism) kết hợp với phương pháp luận cơ học tính toán (computational mechanics). Thiết kế nghiên cứu là sự kết hợp chặt chẽ giữa mô hình hóa toán học giải tích, phát triển thuật toán số trên môi trường MATLAB và kiểm chứng nâng cao bằng mô phỏng số phần tử hữu hạn trên phần mềm chuyên dụng ANSYS.
Hệ thống phương trình vi phân dao động của kết cấu tòa nhà $n$ tầng chịu gia tốc nền địa chấn $\ddot{x}_0(t)$ trang bị hệ thống giằng chủ động (ABS) hoặc thiết bị giảm chấn khối lượng chủ động (AMD) được thiết lập dưới dạng ma trận:
$$[M]{\ddot{x}(t)} + [C]{\dot{x}(t)} + [K]{x(t)} = -[M]{\delta}\ddot{x}_0(t) + [\Gamma]{u(t)}$$
Trong đó:
- $[M], [C], [K]$ lần lượt là các ma trận khối lượng, ma trận cản nhớt Rayleigh và ma trận độ cứng của hệ kết cấu bậc $n \times n$.
- ${x(t)}, {\dot{x}(t)}, {\ddot{x}(t)}$ là các véc-tơ chuyển vị tương đối, vận tốc tương đối và gia tốc tương đối của các tầng sàn.
- ${\delta}_{n \times 1} = [1, 1, \dots, 1]^T$ là véc-tơ ảnh hưởng của gia tốc nền địa chấn.
- ${u(t)}_{r \times 1}$ là véc-tơ lực điều khiển sinh ra bởi $r$ máy kích động ($r \le n$).
- $[\Gamma]_{n \times r}$ là ma trận định vị vị trí lắp đặt của các cơ cấu chấp hành trên các tầng công trình.
Đối với kết cấu lắp đặt AMD trên tầng thượng (tầng thứ $n$), hệ thống bổ sung thêm 1 bậc tự do của khối lượng giảm chấn $m_d$, nâng tổng số bậc tự do của toàn hệ thống lên $(n + 1)$ với tọa độ chuyển vị tương đối của thiết bị là $x_d(t) - x_n(t)$.
Quy trình nghiên cứu rigorous
Quy trình tính toán động lực học công trình và mô phỏng điều khiển hồi tiếp vòng kín (closed-loop feedback control) được thực hiện nghiêm ngặt theo phương pháp số Newmark-$\beta$ (với các tham số tích phân không điều kiện ổn định $\gamma = 0.5$ và $\beta = 0.25$). Kích thước bước thời gian tích phân $\Delta t$ được lựa chọn bằng $0.005\text{ s}$ đến $0.02\text{ s}$, đảm bảo nhỏ hơn nhiều so với chu kỳ dao động riêng thấp nhất của kết cấu, triệt tiêu hoàn toàn sai số phân kỳ số học.
Giao thức thu thập và xử lý dữ liệu địa chấn sử dụng ba bộ bản ghi gia tốc nền chuẩn quốc tế:
- Trận động đất El Centro (1940): Thành phần NS, đỉnh gia tốc nền PGA = $0.348g$ ($3.417\text{ m/s}^2$), thời gian kích động $T = 30\text{ s}$, đại diện cho kích động dải tần rộng với năng lượng xung kích lớn.
- Trận động đất Northridge (1994): PGA = $0.843g$, đại diện cho địa chấn trường gần (near-fault) có vận tốc xung kích cực cao.
- Trận động đất Imperial Valley (1979): PGA = $0.315g$, đại diện cho sóng địa chấn có chu kỳ dài và độ duy trì kéo dài.
Quy trình đánh giá chất lượng điều khiển được định lượng hóa thông qua ba hàm mục tiêu độc lập không thứ nguyên:
- Chỉ tiêu chuyển vị tầng đáy ($F_1$):
$$F_1 = \frac{\max_{t} |d_1(t)|}{d_{1,\max}^{\text{uncontrolled}}}$$
- Chỉ tiêu gia tốc tuyệt đối tầng đỉnh ($F_2$):
$$F_2 = \frac{\max_{t} |\ddot{x}{a,n}(t)|}{\ddot{x}{a,n,\max}^{\text{uncontrolled}}}$$
- Chỉ tiêu năng lượng lực điều khiển trung bình ($F_3$):
$$F_3 = \frac{\frac{1}{T}\int_0^T |u(t)| dt}{u_{\max}}$$
Data và phân tích
Để tìm kiếm tập nghiệm tối ưu cho bộ điều khiển opHAC, luận án ứng dụng Thuật toán Di truyền Đa mục tiêu (Multi-Objective Genetic Algorithm - MOGA) và Tối ưu hóa bầy đàn (PSO). Không gian tìm kiếm được giới hạn chặt chẽ bởi các ràng buộc đại số:
$$f_m(c^-) \in [0.1, 0.9], \quad \mu(h^-) \in [0.1, 0.9]$$
$$w_k \in [0.1, 2.0] \quad (\text{đối với } k = 1, \dots, m)$$
Tập giải pháp Pareto được xây dựng đồng thời trên các không gian hai mục tiêu ($F_1 - F_2$, $F_1 - F_3$, $F_2 - F_3$) và ba mục tiêu ($F_1 - F_2 - F_3$). Dựa vào đường biên Pareto, người kỹ sư thiết kế có thể chủ động lựa chọn cấu hình bộ điều khiển tối ưu tùy thuộc vào yêu cầu công trình: ưu tiên bảo vệ kết cấu chống sụp đổ (cực tiểu hóa $F_1$), ưu tiên bảo vệ tiện nghi người sử dụng và thiết bị thứ cấp nhạy cảm (cực tiểu hóa $F_2$), hay ưu tiên giới hạn công suất máy kích động (cực tiểu hóa $F_3$).
Độ tin cậy của thuật toán HAC còn được kiểm chứng chéo thông qua mô hình phần tử hữu hạn khung không gian 3D trên ANSYS, với các phần tử dầm (BEAM), tấm vỏ (SHELL) và khối đặc (SOLID), nhằm đánh giá sự phân bố của trường ứng suất tương đương von Mises lớn nhất ($\sigma_{\text{eqv}}$) khi có và không có điều khiển chủ động.
Phát hiện đột phá và implications
Những phát hiện then chốt
Các kết quả mô phỏng số và phân tích so sánh chi tiết trong chương 4 của luận án đã đem lại 5 phát hiện đột phá:
-
Hiệu năng triệt tiêu dao động vượt trội của opHAC: Trên mô hình kết cấu nhà 3 tầng chịu tải trọng động đất El Centro, bộ điều khiển tối ưu hóa tham số mờ opHAC đã làm giảm chuyển vị tương đối cực đại của tầng 1 tới $62.4%$ và giảm gia tốc tuyệt đối cực đại của tầng 3 tới $58.7%$ so với trường hợp không điều khiển. So với bộ điều khiển mờ truyền thống FC được thiết kế theo kinh nghiệm, opHAC cải thiện thêm từ $12%$ đến $18%$ hiệu quả giảm chấn trong khi lực điều khiển yêu cầu giảm $15.3%$.
-
Rút ngắn đột phá thời gian tính toán (CPU Time): Một phát hiện có ý nghĩa then chốt cho ứng dụng điều khiển thời gian thực (real-time control) là bộ điều khiển HAC và opHAC giảm từ $40%$ đến $60%$ thời gian xử lý CPU so với FC truyền thống. Bảng so sánh định lượng thời gian tính toán trong luận án chỉ ra rằng do loại bỏ hoàn toàn các vòng lặp tích phân số giải mờ trọng tâm phức tạp, thuật toán HAC cho phép xuất tín hiệu điều khiển gần như tức thời tại mỗi bước trích mẫu $\Delta t$.
-
Tính bền vững cao trước sự biến đổi của phổ địa chấn: Bộ điều khiển opHAC được huấn luyện tối ưu trên dữ liệu động đất El Centro (ký hiệu opHAC_El) khi đem thử nghiệm trực tiếp trên các trận động đất hoàn toàn khác biệt về đặc tính phổ như Northridge (opHAC_No) và Imperial Valley (opHAC_Im) vẫn duy trì hiệu suất giảm chấn đỉnh cao: giảm hơn $50%$ biên độ dao động mà không hề rơi vào trạng thái mất ổn định hay bão hòa máy kích động.
-
Hiệu quả của hệ số điều chỉnh luật trên kết cấu liên tục: Đối với mô hình dầm công-xôn chịu tải trọng xung lực và độ võng ban đầu, việc tối ưu hóa hệ số điều chỉnh $w_k$ (bộ điều khiển tHAC) đã giúp đưa độ võng đầu tự do của dầm về vị trí cân bằng nhanh hơn $45%$ so với bộ điều khiển HAC thông thường, đồng thời triệt tiêu hoàn toàn hiện tượng dao động dư (residual vibration).
-
Giảm ứng suất tập trung trong khung không gian 3D ANSYS: Kết quả mô phỏng đồng thời (co-simulation) trên ANSYS cho thấy lực điều khiển từ opHAC không chỉ làm giảm chuyển vị đỉnh khung không gian từ $38.2\text{ mm}$ xuống còn $14.6\text{ mm}$, mà còn làm giảm ứng suất tương đương von Mises cực đại tại các nút liên kết chân cột từ $245\text{ MPa}$ xuống dưới $130\text{ MPa}$, ngăn ngừa kết cấu đi vào vùng biến dạng dẻo phá hủy.
Implications đa chiều
- Về mặt học thuật và lý thuyết: Nghiên cứu đã chứng minh tính ưu việt vững chắc của lý thuyết Đại số gia tử trong việc giải quyết các bài toán cơ học kỹ thuật phi tuyến; mở ra cầu nối phương pháp luận giữa đại số trừu tượng, logic mờ và động lực học công trình.
- Về đổi mới phương pháp luận: Khung tối ưu hóa đa mục tiêu dựa trên tham số độ đo tính mờ có thể chuyển giao và áp dụng trực tiếp cho nhiều lĩnh vực kỹ thuật khác như điều khiển hệ thống treo chủ động ô tô, tay máy robot công nghiệp, hệ thống ổn định cánh máy bay hoặc điều khiển dao động cáp cầu dây văng.
- Về ứng dụng thực tiễn: Cung cấp cho các kỹ sư xây dựng và cơ khí một giải pháp điều khiển chủ động khả thi cao, thuật toán tinh gọn, dễ dàng nhúng trực tiếp vào các vi điều khiển DSP hoặc hệ thống PLC công nghiệp mà không đòi hỏi phần cứng tính toán đắt tiền.
- Về chính sách và tiêu chuẩn xây dựng: Đóng góp luận cứ khoa học để hoàn thiện các tiêu chuẩn thiết kế kháng chấn cho công trình cao tầng, đề xuất lộ trình tích hợp hệ thống giảm chấn chủ động thông minh vào các quy chuẩn xây dựng quốc gia tại các vùng có nguy cơ động đất.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, luận án cũng thẳng thắn chỉ ra 4 giới hạn nghiên cứu (boundary conditions):
- Phạm vi mô hình tuyến tính: Các mô hình kết cấu khảo sát chủ yếu giả định vật liệu làm việc trong miền đàn hồi tuyến tính, chưa xét đến hiện tượng khớp dẻo, nứt vỡ bê tông hoặc phi tuyến hình học khi kết cấu chịu biến dạng cực lớn.
- Hiện tượng bão hòa máy kích động: Dù đã có các hàm mục tiêu khống chế lực điều khiển, luận án chưa tích hợp sâu bài toán điều khiển trượt thích nghi để xử lý triệt để hiện tượng bão hòa phi tuyến máy kích động (actuator saturation) và độ trễ pha do đường ống thủy lực.
- Giới hạn kiểm chứng thực nghiệm: Toàn bộ các kết quả của luận án được thực hiện thông qua mô phỏng số học và phần tử hữu hạn có độ tin cậy cao, chưa được thử nghiệm trên mô hình bàn rung địa chấn (shaking table test) ở quy mô thực tế do hạn chế về kinh phí trang thiết bị.
- Loại hình điều khiển: Luận án tập trung chủ yếu vào điều khiển chủ động (active control), chưa mở rộng chi tiết sang điều khiển bán chủ động (semi-active) với thiết bị giảm chấn chất lỏng từ biến (Magnetorheological - MR damper).
Từ các giới hạn trên, 5 hướng nghiên cứu phát triển tiếp theo (future research agenda) được đề xuất:
- Phát triển bộ điều khiển Nơ-ron - Mờ dựa trên Đại số gia tử (Neuro-Fuzzy HA Controller) có khả năng tự học trực tuyến (online learning) để thích nghi với các kết cấu có tham số thay đổi theo thời gian do suy thoái mỏi.
- Nghiên cứu tích hợp HAC cho các hệ thống điều khiển bán chủ động sử dụng giảm chấn từ biến MR damper nhằm giảm thiểu tối đa nguồn năng lượng tiêu thụ bên ngoài.
- Mở rộng thuật toán cho bài toán điều khiển bền vững kháng trễ ($H_\infty$ HA control) có xét đến độ trễ thời gian trong mạng lưới cảm biến không dây (wireless sensor networks).
- Khảo sát động lực học phi tuyến hoàn toàn của kết cấu bê tông cốt thép chịu tải trọng động đất vượt cấp kết hợp thuật toán tối ưu hóa thích nghi.
- Tiến hành các thử nghiệm bàn rung thực tế trên mô hình kết cấu nhiều tầng thu nhỏ để xác thực toàn diện thuật toán opHAC trước khi ứng dụng thương mại.
Tác động và ảnh hưởng
Nghiên cứu của luận án tạo ra những tác động sâu rộng trên nhiều phương diện:
- Tác động học thuật (Academic Impact): Tạo tiền đề vững chắc cho việc ứng dụng Đại số gia tử vào cơ học tính toán; dự kiến tạo ra các trích dẫn học thuật quan trọng trong các tạp chí chuyên ngành hàng đầu về Structural Control and Health Monitoring, Journal of Sound and Vibration, và Engineering Structures.
- Chuyển đổi công nghiệp (Industry Transformation): Ngành công nghiệp xây dựng nhà cao tầng, cầu đường và cơ khí giao thông được cung cấp giải pháp giảm chấn tiên tiến, giúp giảm kích thước cấu kiện chịu lực, tiết kiệm hàng triệu USD chi phí vật liệu xây dựng cho mỗi dự án quy mô lớn.
- Lợi ích xã hội (Societal Benefits): Nâng cao mức độ an toàn công trình trước các thảm họa thiên tai động đất và bão lớn, bảo vệ tính mạng con người, giảm thiểu thiệt hại kinh tế do sụp đổ hạ tầng và đảm bảo hoạt động liên tục của các công trình trọng điểm quốc gia (bệnh viện, tháp truyền hình, trung tâm điều hành).
- Ý nghĩa quốc tế (Global Relevance): Khẳng định tính tiên phong và năng lực hội nhập khoa học đỉnh cao của các nhà nghiên cứu Việt Nam trong việc phát triển và ứng dụng lý thuyết toán học Đại số gia tử thuần Việt vào giải quyết các bài toán kỹ thuật phức tạp của thế giới.
Đối tượng hưởng lợi
- Nghiên cứu sinh và học viên cao học ngành Cơ học, Tự động hóa, Xây dựng: Tiếp cận được phương pháp luận nghiên cứu kết hợp giữa toán học hiện đại, cơ học kết cấu và trí tuệ nhân tạo; kế thừa các mã nguồn thuật toán tối ưu opHAC trên MATLAB và quy trình mô phỏng trên ANSYS.
- Các nhà khoa học và chuyên gia đầu ngành: Có thêm bằng chứng thực nghiệm số vững chắc về tính ưu việt của Đại số gia tử so với logic mờ kinh điển, làm nền tảng phát triển các lý thuyết điều khiển thông minh mới.
- Bộ phận R&D của các tập đoàn xây dựng và chế tạo thiết bị giảm chấn: Ứng dụng thuật toán opHAC để nâng cấp phần mềm điều khiển nhúng cho các dòng sản phẩm giảm chấn khối lượng AMD và giằng ABS thế hệ mới.
- Các cơ quan quản lý nhà nước và nhà hoạch định chính sách xây dựng: Sử dụng kết quả nghiên cứu làm căn cứ khoa học để xây dựng, cập nhật các tiêu chuẩn thiết kế kháng chấn quốc gia, khuyến khích ứng dụng công nghệ điều khiển thông minh trong quy hoạch đô thị hiện đạ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?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng Lý thuyết Đại số gia tử (Hedge Algebras) của Nguyễn Cát Hồ sang bài toán tối ưu hóa đa mục tiêu động lực học công trình. Luận án đã đại số hóa toàn bộ quá trình suy luận mờ thành mô hình toán học giải tích tường minh thông qua ánh xạ ngữ nghĩa định lượng SQM, đồng thời phát triển phương pháp tối ưu hóa tham số độ đo tính mờ $f_m(c^-), \mu(h^-)$ trên tập Pareto, loại bỏ hoàn toàn việc sử dụng hàm thuộc mờ nhân tạo của Zadeh.
2. Điểm đổi mới về phương pháp luận khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?
Trả lời: So với nghiên cứu của Marinaki et al. (2011) và Uz & Khuc (2013) – vốn sử dụng GA/PSO để tối ưu hóa $3n$ tọa độ đỉnh hàm thuộc mờ tam giác gây bùng nổ số chiều và tiềm ẩn nguy cơ đảo lộn thứ tự ngữ nghĩa – phương pháp của luận án chỉ sử dụng từ 2 đến 4 biến thiết kế đại số độc lập, đảm bảo tính đơn điệu ngữ nghĩa tự nhiên 100%. So với nghiên cứu của Guclu (2006) sử dụng bảng luật cố định, luận án đã đổi mới bằng cách tích hợp hệ số điều chỉnh luật $w_k$, cho phép biến dạng linh hoạt mặt SQS thích ứng với động lực học kết cấu.
3. Phát hiện bất ngờ nhất trong quá trình tính toán mô phỏng là gì?
Trả lời: Phát hiện bất ngờ nhất là việc giảm thiểu đáng kể thời gian tính toán CPU time (từ $40%$ đến $60%$) của bộ điều khiển opHAC so với FC truyền thống mà không hề làm suy giảm độ chính xác điều khiển. Điều này chứng minh rằng việc đại số hóa logic mờ không những không làm mất mát thông tin mà còn tối ưu hóa tài nguyên phần cứng vượt bậc.
4. Luận án có cung cấp giao thức tái lập (replication protocol) đầy đủ không?
Trả lời: Luận án cung cấp đầy đủ các phương trình vi phân chuyển động, ma trận độ cứng, khối lượng, cản của các hệ 1-DOF, 3-DOF, 5-DOF, thông số dầm công-xôn, dữ liệu kích động địa chấn chuẩn hóa (El Centro, Northridge, Imperial Valley), cùng thuật toán Newmark-$\beta$ và các công thức ánh xạ ngữ nghĩa định lượng SQM chi tiết, cho phép các nhà nghiên cứu độc lập tái lập chính xác toàn bộ kết quả mô phỏng trên MATLAB và ANSYS.
5. Lộ trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?
Trả lời: Lộ trình 10 năm bao gồm: Giai đoạn 1 (1-3 năm): Phát triển bộ điều khiển mờ - nơ ron tự thích nghi HA và tích hợp vào giảm chấn bán chủ động MR damper; Giai đoạn 2 (4-6 năm): Thực hiện các thí nghiệm bàn rung quy mô lớn và giải quyết bài toán bù trễ điều khiển mạng cảm biến không dây; Giai đoạn 3 (7-10 năm): Thương mại hóa hệ thống nhúng vi xử lý opHAC và tích hợp vào các tiêu chuẩn thiết kế nhà thông minh kháng chấn toàn cầu.
Kết luận
- Luận án đã giải quyết thành công bài toán điều khiển chủ động dao động kết cấu bằng cách kết hợp nhuần nhuyễn Lý thuyết mờ, Đại số gia tử và Kỹ thuật tối ưu hóa đa mục tiêu.
- Đã đề xuất phương án tối ưu hóa đa mục tiêu đột phá với biến thiết kế là các tham số độ đo tính mờ $f_m(c^-)$ và $\mu(h^-)$ của các biến trạng thái và biến điều khiển, giúp thu nhỏ không gian tìm kiếm và đảm bảo tính nhất quán ngữ nghĩa tuyệt đối.
- Đã sáng tạo ra phương pháp sử dụng hệ số điều chỉnh luật $w_k$ để tối ưu hóa hệ luật định tính của HAC, giúp bề mặt điều khiển thích nghi hoàn hảo với các dạng kết cấu khác nhau từ hệ rời rạc nhiều tầng đến hệ liên tục dầm công-xôn và khung không gian 3D.
- Chứng minh bằng số liệu định lượng rằng bộ điều khiển opHAC làm giảm tới hơn $60%$ chuyển vị đỉnh và gia tốc công trình dưới tác động của các trận động đất dữ dội (El Centro, Northridge, Imperial Valley), đồng thời giảm từ $40%$ đến $60%$ thời gian tính toán CPU so với bộ điều khiển mờ truyền thống.
- Mở ra 3 nhánh nghiên cứu chuyên sâu mới: Điều khiển bán chủ động thông minh với giảm chấn MR dựa trên HA, Điều khiển thích nghi mờ - nơ ron HA cho kết cấu phi tuyến, và Thiết kế kháng chấn công trình tối ưu trên không gian Pareto.
- Toàn bộ kết quả của luận án tạo nên một bước tiến quan trọng về mặt học thuật và thực tiễn, đóng góp cơ sở khoa học vững chắc để nâng cao độ an toàn và khả năng chống chịu động đất của các công trình xây dựng trong kỷ nguyên công trình thông minh.