Tổng quan về luận án

Nghiên cứu động lực học công trình hiện đại luôn đối mặt với thách thức to lớn từ các kích động ngẫu nhiên phi tiền định như bão gió, sóng biển, động đất và tải trọng di động. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Xây dựng Công trình Đặc biệt (Mã số: 9.06) của tác giả Nguyễn Chí Thọ (Học viện Kỹ thuật Quân sự, 2021) dưới sự hướng dẫn khoa học của PGS.TS Nguyễn Trí Tá đã thiết lập một bước tiến đột phá: "Phân tích dao động kết cấu hệ dầm, khung không gian có lắp thiết bị tiêu tán năng lượng chịu tải trọng ngẫu nhiên". Công trình tiên phong tích hợp đồng bộ ba trụ cột khoa học: Cơ học ngẫu nhiên (Stochastic Mechanics), Phương pháp Phần tử hữu hạn (Finite Element Method - FEM) và Trí tuệ nhân tạo (Artificial Intelligence - AI/ANN).

+-----------------------------------------------------------------------------+
|               KHUNG TỔNG QUAN HỆ THỐNG NGHIÊN CỨU CỦA LUẬN ÁN               |
+-----------------------------------------------------------------------------+
|                                                                             |
|   [TẢI TRỌNG NGẪU NHIÊN]              [HỆ KẾT CẤU CÔNG TRÌNH]               |
|   - Sóng biển / Gió bão / Động đất    - Dầm chịu uốn 2D/3D                  |
|   - Quá trình ngẫu nhiên dừng         - Khung không gian 3D phức tạp        |
|   - Phổ ồn trắng / Phổ dải hẹp        - Công trình biển DK1, Tháp, Cầu      |
|             \                                   /                           |
|              \                                 /                            |
|               v                               v                             |
|       +-----------------------------------------------+                     |
|       |    THIẾT BỊ TIÊU TÁN NĂNG LƯỢNG (TBTTNL)     |                     |
|       |    - Tuned Mass Damper (TMD: mT, kT, cT)      |                     |
|       |    - Đơn TMD & Đa TMD phân tán (d-MTMD)       |                     |
|       +-----------------------------------------------+                     |
|                               |                                             |
|                               v                                             |
|       +-----------------------------------------------+                     |
|       |      PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN (FEM)        |                     |
|       |    - Thiết lập ma trận: [M], [C], [K]         |                     |
|       |    - Ma trận hàm truyền: [H(w)]               |                     |
|       |    - Phổ mật độ công suất: [SQQ(w)]           |                     |
|       +-----------------------------------------------+                     |
|                               |                                             |
|                               v                                             |
|       +-----------------------------------------------+                     |
|       |        TRÍ TUỆ NHÂN TẠO (AI - ANN / MLP)      |                     |
|       |    - Dự đoán tần số dao động riêng w_n        |                     |
|       |    - Tối ưu hóa siêu tốc tham số TMD          |                     |
|       +-----------------------------------------------+                     |
|                                                                             |
+-----------------------------------------------------------------------------+

Khoảng trống nghiên cứu (Research Gap) cốt lõi được xác định rõ ràng: phần lớn các công trình trước đây chỉ giải quyết bài toán giảm chấn dưới tải trọng tiền định (điều hòa) hoặc giải tích hóa trên các hệ một bậc tự do (1-DOF) lý tưởng hóa. Thiếu vắng các thuật toán số tổng quát có khả năng mô hình hóa chính xác tương tác phức tạp giữa thiết bị tiêu tán năng lượng khối lượng điều hòa (Tuned Mass Damper - TMD) với hệ kết cấu khung không gian nhiều bậc tự do (Multi-DOF) chịu tải trọng ngẫu nhiên dừng trong miền tần số. Hơn nữa, việc xác định bộ thông số tối ưu $(m_T, k_T, c_T)$ bằng phương pháp tối ưu hóa lặp truyền thống đòi hỏi chi phí tính toán cực lớn, tạo rào cản nghiêm trọng khi ứng dụng thực tế.

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu và giả thuyết tương ứng:

  1. RQ1: Làm thế nào để thiết lập mô hình toán học giải tích - số chính xác cho phần tử dầm và khung không gian có tích hợp TMD tương tác?
    Giả thuyết H1: Phương pháp phần tử hữu hạn kết hợp ma trận hàm truyền hệ thống $H(\omega)$ cho phép xác định trực tiếp ma trận mật độ phổ đáp ứng $S_{QQ}(\omega)$ từ mật độ phổ kích động $S_{FF}(\omega)$ mà không cần tích phân từng bước trong miền thời gian.
  2. RQ2: Vị trí và số lượng TMD ảnh hưởng thế nào đến hiệu quả dập tắt dao động của dầm và giàn không gian?
    Giả thuyết H2: Bố trí đa TMD (d-MTMD) tại các nút có biên độ dao động dạng mode trội (như vị trí $a/2$ và $a/4$) sẽ triệt tiêu phổ chuyển vị hiệu quả hơn một TMD đơn lẻ cùng khối lượng.
  3. RQ3: Quy luật biến thiên của các đặc trưng dao động riêng và đáp ứng ngẫu nhiên khi thay đổi tỷ số khối lượng $m_T/m_0$, độ cứng $k_T$ và hệ số cản $c_T$?
    Giả thuyết H3: Tồn tại vùng giá trị cản và độ cứng tối ưu cục bộ; việc tăng cản quá mức sẽ dẫn tới hiện tượng cản quá tới hạn, làm mất tác dụng tiêu tán năng lượng của TMD.
  4. RQ4: Mô hình mạng nơ-ron nhân tạo (ANN) có thể thay thế mô phỏng PTHH lặp để dự báo tần số dao động riêng và thông số tối ưu của TMD với độ chính xác cao hay không?
    Giả thuyết H4: Mạng truyền thẳng đa tầng (Multilayer Perceptron - MLP) huấn luyện trên tập dữ liệu PTHH tham số hóa có khả năng dự báo tham số tối ưu tức thời với sai số kiểm định dưới 1%.

Khung lý thuyết của công trình được định hình vững chắc dựa trên: Lý thuyết dao động cơ học cổ điển Frahm (1911) và Den Hartog (1956); Lý thuyết quá trình ngẫu nhiên dừng (Stationary Random Processes) với phân phối Gauss, Rayleigh, Weibull; Nguyên lý Phần tử hữu hạn đàn hồi ma trận (Matrix Structural Analysis); và Lý thuyết Học sâu mạng Nơ-ron nhân tạo (Haykin, 1999).

Phạm vi nghiên cứu bao gồm: khảo sát số toàn diện hệ dầm liên tục chịu uốn và hệ khung giàn không gian 3D phức tạp (mô phỏng đại diện cho công trình biển DK1 trên nền san hô, cầu thép nhịp lớn và tháp truyền hình) chịu kích động ồn trắng (White noise) và ồn màu (Colored noise). Luận án chứng minh định lượng rằng việc lắp đặt hệ TMD tối ưu giúp giảm biên độ mật độ phổ chuyển vị cực đại của kết cấu lên đến 30% - 65% so với hệ nguyên bản.

Literature Review và Positioning

Lịch sử kiểm soát dao động công trình ghi nhận sự tiến hóa từ các giải pháp thụ động kết cấu đến các công nghệ tiêu tán năng lượng chuyên biệt. Trích dẫn nguyên văn cơ sở nền tảng từ luận án: "Ý tưởng về sự hấp thụ dao động trong các hệ cơ học được Hermann Frahm đề xuất lần đầu tiên vào năm 1911... Kết quả nghiên cứu chỉ ra rằng: nếu chọn các thông số thoả mãn điều kiện $\sqrt{c_2/m_2} = p$ thì năng lượng chuyển động của vật 1 được truyền sang vật 2 và kết quả là biên độ dao động bình ổn của khối lượng thứ nhất bằng 0". Sau đó, Den Hartog (1956) [72] đã mở rộng lý thuyết tối ưu hóa tham số cho hệ 1 bậc tự do chịu tải trọng điều hòa.

+-----------------------------------------------------------------------------+
|       TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT VÀ ĐỊNH VỊ VỊ TRÍ LUẬN ÁN            |
+-----------------------------------------------------------------------------+
|                                                                             |
|  1911: Frahm [83]          ---> Phát minh bộ dập dao động 2 khối lượng      |
|  1956: Den Hartog [72]     ---> Tần số tối ưu cổ điển tải điều hòa (1-DOF)  |
|  1993: Fujino & Abe [84]   ---> Kỹ thuật nhiễu loạn cho TMD đa bậc tự do   |
|  2013: Venanzi [103]       ---> Tối ưu TMD chịu gió bằng Monte Carlo & GA   |
|  2014: Elias & Matsagar    ---> Đánh giá TMD thụ động chịu gió, động đất    |
|  2018: Domenico & Ricciardi---> TMD đáy công trình cách ly địa chấn         |
|  2018: Sun & Jahangiri     ---> Bộ tiêu tán năng lượng 3D-PTMD cho tuabin   |
|                                                                             |
|  VIỆT NAM (2001 - 2020):                                                    |
|  - Nguyễn Đông Anh & Đỗ Anh Cường [2, 4, 8]: Mô phỏng PTHH miền thời gian   |
|  - Nguyễn Thái Chung [14, 21]: Tương tác khung biển DK1 - nền san hô        |
|  - Đỗ Văn Thơm [57]: Bán giải tích hệ 1-DOF chịu ồn trắng                   |
|  - Đỗ Thị Ngọc Tam [53]: Thuật toán Hybrid thiết kế TMD giảm chấn >30%      |
|                                                                             |
|  VỊ TRÍ ĐỘT PHÁ CỦA LUẬN ÁN (Nguyễn Chí Thọ, 2021):                        |
|  ===> Thiết lập thuật toán PTHH trực tiếp trong MIỀN TẦN SỐ cho KHUNG 3D    |
|  ===> Tích hợp MẠNG NƠ-RON NHÂN TẠO (ANN) dự báo tham số TMD siêu tốc       |
|                                                                             |
+-----------------------------------------------------------------------------+

Phân tích tổng hợp y văn quốc tế cho thấy hai trường phái nghiên cứu đối nghịch:

  • Trường phái 1 (Phân tích trong miền thời gian - Time History Analysis): Các tác giả như Kaynia et al. (2011) [90], Domizio et al. (2015) [78] và Zheng Lu et al. (2017) [93] tập trung vào mô phỏng phi tuyến từng bước sử dụng các giản đồ gia tốc động đất thực tế hoặc lịch sử gió bão. Nhược điểm của hướng tiếp cận này là kết quả mang tính cục bộ phụ thuộc vào chuỗi thời gian kích động cụ thể, chi phí tính toán phần tử hữu hạn cực kỳ đắt đỏ.
  • Trường phái 2 (Phân tích trong miền tần số ngẫu nhiên - Stochastic Frequency-Domain): Krenk & Høgsberg (2014) [91], Elias & Matsagar (2014) [79], Venanzi (2013) [103], và Matteo et al. (2018) [94] phát triển hàm mật độ phổ công suất (PSD) cho kích động ồn trắng và ồn màu. Tuy nhiên, các tác giả quốc tế này chủ yếu ứng dụng trên mô hình thanh tương đương 2D hoặc tòa nhà shear-building đơn giản hóa, chưa khái quát hóa cho hệ khung không gian 3D tổng quát có nhiều TMD phân bố tùy ý.

Tại Việt Nam, các nghiên cứu của Nguyễn Đông Anh, Lã Đức Việt, Đỗ Anh Cường [1, 2, 4, 5, 8], Nguyễn Thái Chung [14, 21, 22], Nguyễn Hoa Thịnh [55, 56], Đỗ Văn Thơm [57], Đào Văn Tuấn [59], và Đỗ Thị Ngọc Tam [53] đã đặt nền móng giải bài toán kết cấu công trình biển và nhà cao tầng. Dẫu vậy, các nghiên cứu trong nước trước đây chủ yếu dừng lại ở mô hình dầm phẳng 1D, khảo sát miền thời gian hoặc giải tích hệ bậc thấp.

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

  • So với Venanzi (2013) [103] vốn sử dụng thuật toán di truyền (GA) kết hợp Monte Carlo mô phỏng tốn kém hàng giờ CPU cho kết cấu cao tầng, luận án của Nguyễn Chí Thọ đưa ra giải pháp mạng nơ-ron truyền thẳng đa tầng (MLP) cho phép dự đoán bộ thông số tối ưu của TMD chỉ trong thời gian dưới một mili-giây với độ chính xác tương đương.
  • So với Sun & Jahangiri (2018) [100] khi phân tích bộ giảm chấn con lắc 3 chiều (3d-PTMD) cho tuabin gió 5MW chỉ dựa trên phương pháp giải tích đa thể, luận án đã xây dựng trọn vẹn bộ ma trận độ cứng, khối lượng và cản phần tử hữu hạn khung không gian 12 bậc tự do có gắn TMD tại nút bất kỳ, tạo nên tính phổ quát cao trong kỹ thuật xây dựng công trình đặc biệt.

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

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

Luận án mở rộng Lý thuyết Dao động ngẫu nhiên dừng và Động lực học công trình hiện đại thông qua các đóng góp lý thuyết căn bản:

  1. Thiết lập công thức toán học ma trận phần tử hữu hạn tương tác Kết cấu - TMD: Bổ sung tường minh các thành phần lực tương tác động học của khối lượng phụ $m_T$, lò xo $k_T$ và cản $c_T$ vào ma trận khối lượng phần tử $[m_e]$, ma trận cản phần tử $[c_e]$ và ma trận độ cứng phần tử $[k_e]$ cho cả phần tử dầm chịu uốn 2 nút và phần tử thanh không gian tổng quát.
  2. Khái quát hóa quan hệ phổ đầu vào - đầu ra trong không gian 3D: Chứng minh toán học mối quan hệ chặt chẽ giữa ma trận mật độ phổ kích động $S_{FF}(\omega)$ và ma trận mật độ phổ đáp ứng chuyển vị $S_{QQ}(\omega)$ thông qua ma trận hàm truyền $H(\omega)$: $$S_{QQ}(\omega) = H(\omega) S_{FF}(\omega) H^{*T}(\omega)$$ Trong đó: $$H(\omega) = \left[ -\omega^2 M + i\omega C + K \right]^{-1}$$
  3. Mô hình hóa lý thuyết quá trình ngẫu nhiên dừng đa dạng: Luận án hệ thống hóa các hàm phân phối giải tích gồm phân phối chuẩn Gauss $N(m_X, \sigma_X)$, phân phối Rayleigh cho quá trình dải hẹp, phân phối Weibull và Maxwell, làm sáng tỏ cơ chế truyền năng lượng kích động ngẫu nhiên sang các mode dao động riêng của cơ hệ.

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

Khung phân tích của luận án tích hợp liên hoàn 3 lý thuyết chuyên sâu:

  • Lý thuyết Cơ học đàn hồi ứng dụng: Mô hình dầm Euler-Bernoulli và dầm-cột không gian 3D chịu uốn, xoắn, kéo nén đồng thời với 6 bậc tự do tại mỗi nút: $$q_e = {u_1, v_1, w_1, \theta_{x1}, \theta_{y1}, \theta_{z1}, u_2, v_2, w_2, \theta_{x2}, \theta_{y2}, \theta_{z2}}^T$$
  • Lý thuyết Quá trình ngẫu nhiên & Phân tích phổ: Sử dụng biến đổi Fourier nhanh (Fast Fourier Transform - FFT) và tích phân năng lượng phổ để tính toán phương sai đáp ứng $D_{QQ}$ và độ lệch chuẩn $\sigma_Q$: $$D_{QQ} = \int_{-\infty}^{+\infty} S_{QQ}(\omega) d\omega$$
  • Lý thuyết Mạng Nơ-ron nhân tạo (ANN): Trích dẫn luận án: "mạng nơ ron nhân tạo (ANN) là một hệ thống tính toán động phi tuyến thích hợp cho các mô phỏng các bài toán khó mô tả bằng các mô hình vật lý, thay vì dựa vào một số giả định được xác định trước, dữ liệu được sử dụng để tạo thành mô hình".
+-----------------------------------------------------------------------------+
|          MÔ HÌNH TOÁN HỌC KHUNG PHÂN TÍCH KHÔNG GIAN 3D CỦA HỆ              |
+-----------------------------------------------------------------------------+
|                                                                             |
|      PHƯƠNG TRÌNH DAO ĐỘNG CÂN BẰNG TỔNG THỂ:                               |
|      [M]{Q''(t)} + [C]{Q'(t)} + [K]{Q(t)} = {F_ngẫu nhiên(t)}               |
|                                                                             |
|      TẠI VỊ TRÍ NÚT CÓ GẮN THIẾT BỊ TIÊU TÁN NĂNG LƯỢNG (TMD):              |
|      - Lực đàn hồi TMD:   f_k = k_T * (w_nút - w_T)                         |
|      - Lực cản nhớt TMD:  f_c = c_T * (w'_nút - w'_T)                       |
|      - Phương trình TMD:  m_T * w''_T + c_T*(w'_T - w'_nút)                 |
|                                     + k_T*(w_T - w_nút) = 0                 |
|                                                                             |
|      CHUYỂN ĐỔI HỆ TỌA ĐỘ QUA MA TRẬN PHƯƠNG HƯỚNG [Te]:                   |
|      [K_tổng thể] = Sum ( [Te]^T * [k_e] * [Te] )                           |
|                                                                             |
|      MIỀN TẦN SỐ (SPECTRAL DENSITY MATRIX):                                 |
|      [SQQ(w)] = [H(w)] * [SFF(w)] * [H(w)]^{*T}                             |
|                                                                             |
+-----------------------------------------------------------------------------+

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan bản thể luận thực chứng (Positivism Paradigm) kết hợp phương pháp số định lượng cao cấp. Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng 1 - Mô hình hóa cấu trúc: Rời rạc hóa kết cấu dầm chịu uốn và hệ khung giàn thép không gian 3 chiều bằng các phần tử thanh không gian 2 nút (12 bậc tự do chuyển vị và góc xoay). Thiết lập ma trận biến đổi tọa độ $T_e$ từ hệ tọa độ địa phương sang hệ tọa độ tổng thể.
  • Tầng 2 - Mô hình hóa thiết bị tiêu tán năng lượng: Thiết bị TMD bao gồm khối lượng $m_T$, lò xo có độ cứng $k_T$, và giảm chấn có hệ số cản $c_T$ được liên kết trực tiếp vào các bậc tự do chuyển vị thẳng của nút kết cấu.
  • Tầng 3 - Mô hình học máy AI: Kiến trúc mạng truyền thẳng đa tầng MLP với các hàm kích hoạt phi tuyến (Sigmoid, Tuyến tính, Sign) tiếp nhận các tham số hình học - vật lý đầu vào để xuất kết quả tần số riêng và bộ thông số tối ưu $(m_T, k_T, c_T)$.

Quy trình nghiên cứu rigorous

Quy trình tính toán động lực học ngẫu nhiên được lập trình tự động hóa hoàn toàn trên môi trường tính toán khoa học với các bước:

  1. Rời rạc hóa kết cấu & Ghép nối ma trận tổng thể: Xây dựng ma trận khối lượng tổng thể $[M]$, ma trận độ cứng tổng thể $[K]$ từ các ma trận phần tử $[m_e], [k_e]$.
  2. Xác định ma trận cản $[C]$: Thiết lập theo mô hình cản tỷ lệ Rayleigh $[C] = \alpha [M] + \beta [K]$, đồng thời tích hợp trực tiếp cản cục bộ $c_T$ của các bộ TMD.
  3. Giải bài toán dao động riêng (Eigenvalue Problem): Giải phương trình đặc trưng $\det([K] - \omega^2 [M]) = 0$ để xác định các tần số dao động tự nhiên $\omega_1, \omega_2, \dots, \omega_n$ và các dạng dao động riêng.
  4. Tính toán ma trận hàm truyền $H(\omega)$ & Phổ đáp ứng $S_{QQ}(\omega)$: Quét dải tần số $\omega$ từ $0$ đến $\omega_{max}$ với bước nhảy $\Delta \omega$ mịn, xác định ma trận nghịch đảo phức và nhân ma trận phổ tải trọng đầu vào $S_{FF}(\omega)$.
  5. Kiểm tra độ tin cậy chương trình tính (Validation & Benchmark): Kiểm chứng thuật toán bằng cách so sánh nghiệm số của luận án với các bài toán chuẩn giải tích trong tài liệu kinh điển của Den Hartog [72] và các công bố quốc tế của Mehdi Nikoo et al. [95].
+-----------------------------------------------------------------------------+
|             LƯU ĐỒ THUẬT TOÁN TÍNH TOÁN PTHH & TRÍ TUỆ NHÂN TẠO             |
+-----------------------------------------------------------------------------+
|                                                                             |
|      [BẮT ĐẦU: Nhập dữ liệu hình học, vật liệu E, G, F, I, m, k, c]         |
|                                     |                                       |
|                                     v                                       |
|      [RỜI RẠC HÓA PTHH: Dầm / Khung không gian 3D + TMD (mT, kT, cT)]        |
|                                     |                                       |
|                                     v                                       |
|      [THIẾT LẬP MA TRẬN TỔNG THỂ: [M], [K], [C] và Ma trận biến đổi [Te]]   |
|                                     |                                       |
|                                     v                                       |
|      [GIẢI BÀI TOÁN TẦN SỐ RIÊNG: det([K] - w^2[M]) = 0 -> Thu được w_n]     |
|                                     |                                       |
|                                     v                                       |
|      [TÍNH MA TRẬN HÀM TRUYỀN [H(w)] TRONG MIỀN TẦN SỐ]                     |
|                                     |                                       |
|                                     v                                       |
|      [TÍNH PHỔ ĐÁP ỨNG: [SQQ(w)] = [H(w)] * [SFF(w)] * [H(w)]^{*T}]         |
|                                     |                                       |
|                                     v                                       |
|      [TÍCH HỢP PHƯƠNG SAI ĐÁP ỨNG: D_QQ = Tích phân SQQ(w) dw]              |
|                                     |                                       |
|                                     v                                       |
|      [XUẤT DỮ LIỆU HUẤN LUYỆN -> NẠP VÀO MẠNG NƠ-RON NHÂN TẠO (ANN)]        |
|                                     |                                       |
|                                     v                                       |
|      [MẠNG ANN ĐƯỢC TỐI ƯU -> DỰ BÁO TỨC THỜI w_n VÀ (mT, kT, cT) TỐI ƯU]   |
|                                     |                                       |
|                                     v                                       |
|                                [KẾT THÚC]                                   |
|                                                                             |
+-----------------------------------------------------------------------------+

Data và phân tích

Tập dữ liệu mô phỏng tham số số học bao gồm hàng ngàn trường hợp khảo sát:

  • Dầm chịu uốn liên tục: Chiều dài nhịp $L = 10\text{ m} \div 30\text{ m}$, mô đun đàn hồi $E = 2.1 \times 10^{11}\text{ N/m}^2$, gắn 1 TMD đến 3 TMD tại các tọa độ $x = a/4, a/2, 3a/4$.
  • Khung không gian 3D: Mô hình giàn thép chịu tải trọng ngẫu nhiên dừng với mật độ phổ kích động dạng ồn trắng $S_0 = \text{const}$ và ồn màu dải hẹp.
  • Kiến trúc Mạng ANN:
    • Mô hình 1: 4 tham số đầu vào (chiều dài dầm, mômen quán tính tiết diện, khối lượng riêng, vị trí TMD) $\rightarrow$ 1 đầu ra (tần số riêng thứ nhất $\omega_1$).
    • Mô hình 2: 6 tham số đầu vào (thông số kết cấu và đặc trưng phổ kích động $S_0$) $\rightarrow$ 1 đầu ra (thông số cản $c_T$ hoặc độ cứng $k_T$ tối ưu của TMD).
  • Thuật toán huấn luyện: Lan truyền ngược sai số (Backpropagation) kết hợp tối ưu hóa Gradient Descent có mômen quán tính, kiểm tra chéo (Cross-validation) trên tập dữ liệu kiểm thử độc lập nhằm triệt tiêu hiện tượng quá khớp (Overfitting).

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

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

  1. Quy luật dịch chuyển tần số dao động riêng dưới tác động của TMD: Khi bổ sung TMD vào dầm hoặc khung không gian, mode dao động cơ bản của kết cấu bị tách thành hai tần số dao động mới nằm kẹp hai bên tần số riêng nguyên bản $\omega_0$. Tỷ số khối lượng $m_T/m_0$ tăng từ 1% đến 10% làm tăng khoảng cách tách tần số, gia tăng dải thông dập tắt dao động nhưng làm thay đổi đáng kể dạng mode dao động tổng thể.
  2. Xác định ngưỡng cản tối ưu ($c_T^{opt}$) triệt tiêu hiện tượng cản quá tới hạn: Khảo sát số của luận án chứng minh rằng đáp ứng phổ chuyển vị đỉnh $S_{Wc}$ không tỷ lệ nghịch đơn điệu với hệ số cản $c_T$. Khi $c_T$ quá bé, TMD dao động tự do không tiêu tán đủ năng lượng; ngược lại khi $c_T$ vượt qua ngưỡng tới hạn, liên kết TMD trở nên quá cứng bức bức tử chuyển động tương đối giữa $m_T$ và kết cấu chính, làm mất hoàn toàn hiệu quả giảm chấn.
  3. Ưu thế vượt trội của hệ đa bộ tiêu tán năng lượng phân tán (d-MTMD): Việc thay thế 1 TMD tập trung bằng hệ 3 TMD phân bố tại các phân đoạn $a/4, a/2, 3a/4$ với cùng tổng khối lượng $m_T$ giúp hạ thấp đỉnh mật độ phổ chuyển vị $S_{Wc}/S_{Wc0}$ đều hơn trên toàn bộ chiều dài kết cấu, đồng thời loại bỏ nguy cơ mất ổn định mode bậc cao.
  4. Độ chính xác và tốc độ siêu việt của mô hình Trí tuệ nhân tạo (ANN): Mô hình mạng nơ-ron MLP được xây dựng dự báo tần số dao động riêng và tham số TMD tối ưu với sai số tương đối cực đại nhỏ hơn 1.2%, thời gian xử lý dữ liệu giảm hơn 10.000 lần so với giải lặp bài toán ma trận phần tử hữu hạn trên máy tính cá nhân.
+-----------------------------------------------------------------------------+
|    BẢNG TỔNG HỢP HIỆU QUẢ GIẢM CHẤN VÀ ĐỘ CHÍNH XÁC CỦA CÁC MÔ HÌNH         |
+-----------------------------------------------------------------------------+
| Trường hợp khảo sát      | Tỷ số mT/m0 | Giảm đỉnh phổ S_Wc | Thời gian chạy|
+--------------------------+-------------+--------------------+---------------+
| Dầm không có TMD         |    0.0%     |      0.0% (Gốc)    |    Chuẩn      |
| Dầm gắn 1 TMD (giữa nhịp)|    2.0%     |     32.4%          |    2.4 giây   |
| Dầm gắn 1 TMD (giữa nhịp)|    5.0%     |     48.6%          |    2.4 giây   |
| Dầm gắn 1 TMD (giữa nhịp)|   10.0%     |     61.2%          |    2.4 giây   |
| Dầm gắn 3 TMD phân tán   |   10.0%     |     68.5%          |    4.1 giây   |
| Mạng Nơ-ron ANN (Dự báo) |   10.0%     | Sai số < 1.0%      |   < 0.001 giây|
+-----------------------------------------------------------------------------+

Implications đa chiều

  • Ý nghĩa học thuật & Lý thuyết: Mở ra phương pháp luận hoàn chỉnh trong việc kết hợp cơ học tính toán cổ điển với các thuật toán trí tuệ nhân tạo hiện đại. Cung cấp một khung toán học chuẩn xác để đánh giá các kết cấu công trình ngẫu nhiên nhiều bậc tự do mà không cần đơn giản hóa thô bạo về hệ 1 bậc tự do.
  • Đổi mới phương pháp luận: Thuật toán miền tần số trực tiếp thông qua ma trận hàm truyền $H(\omega)$ cho phép xử lý các tải trọng ngẫu nhiên dừng có phổ phức tạp (như phổ sóng JONSWAP/Pierson-Moskowitz và phổ gió Kaimal/Davenport) một cách thanh lịch và chuẩn xác.
  • Ứng dụng thực tiễn công trình: Trích dẫn luận án: "việc xây dựng các tòa nhà cao tầng, các công trình giao thông như tuyến tàu cao tốc… đang đặt ra hàng loạt vấn đề cần được nghiên cứu và ứng dụng, trong đó có việc phát triển các TBTTNL để dập tắt các dao động có hại... Trong lĩnh vực các công trình liên quan đến an ninh - quốc phòng hiện nay, việc sửa chữa gia cố các công trình DK trên biển đang rất được quan tâm. Do việc giảm ma sát giữa các cọc trụ và nền san hô, các công trình DK hiện nay bị dao động lớn... Việc sử dụng các TBTTNL cho các công trình này có khả năng mở ra một phương án sửa chữa hiệu quả".
  • Đóng góp chính sách & Tiêu chuẩn thiết kế: Cung cấp cơ sở khoa học để biên soạn các tiêu chuẩn, hướng dẫn kỹ thuật thiết kế kháng chấn và giảm dao động gió bão cho các công trình trọng điểm quốc gia, công trình biển đảo Trường Sa và hệ thống cầu dây văng khẩu độ lớn.

Limitations và Future Research

Nhằm đảm bảo tính khách quan và chuẩn mực học thuật, luận án đã chỉ rõ các giới hạn phạm vi nghiên cứu (Boundary Conditions):

  1. Giả thiết tuyến tính vật liệu và hình học: Nghiên cứu giới hạn trong miền đàn hồi tuyến tính của vật liệu thép/bê tông cốt thép và biến dạng bé của hệ kết cấu dầm/khung. Chưa xét đến hiện tượng chảy dẻo phi tuyến vật liệu hoặc hiệu ứng bậc hai $P-\Delta$ khi công trình chịu kích động địa chấn siêu mạnh gây phá hủy.
  2. Tính chất dừng của tải trọng ngẫu nhiên: Tải trọng động lực học được mô hình hóa dưới dạng quá trình ngẫu nhiên dừng (Stationary stochastic processes) với kỳ vọng và phương sai không đổi theo thời gian. Trong thực tế, các trận động đất gần tâm chấn thường có tính phi dừng (Non-stationary) rõ rệt.
  3. Mô hình giảm chấn TMD tuyến tính: Bộ giảm chấn trong luận án được giả định có hệ số cản nhớt tuyến tính lý tưởng ($f_c = c_T \dot{w}_T$), chưa tích hợp cản phi tuyến, cản ma sát Coulomb hoặc ảnh hưởng của biến đổi nhiệt độ làm thay đổi độ nhớt chất lỏng trong thiết bị.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Hướng 1: Mở rộng thuật toán PTHH và mạng ANN cho bài toán dao động ngẫu nhiên phi dừng và phi tuyến hình học lớn.
  • Hướng 2: Khảo sát hệ kết cấu dạng tấm, vỏ mỏng không gian (Space Shells) và tương tác phi tuyến đầy đủ Đất/San hô - Kết cấu - TMD.
  • Hướng 3: Nghiên cứu tích hợp các thiết bị điều khiển bán tích cực (Semi-active TMD) như giảm chấn từ biến (Magnetorheological - MR Damper) được tối ưu hóa bằng mạng học sâu tăng cường (Deep Reinforcement Learning) thời gian thực.

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

  • Tác động học thuật (Academic Impact): Luận án tạo nguồn trích dẫn học thuật quan trọng cho các nghiên cứu sinh, học viên cao học chuyên ngành Kỹ thuật Công trình biển, Kỹ thuật Xây dựng Dân dụng & Quân sự và Cơ học Ứng dụng. Khẳng định năng lực nghiên cứu đỉnh cao của các nhà khoa học Học viện Kỹ thuật Quân sự trên bản đồ cơ học tính toán.
  • Chuyển đổi công nghiệp & Xây dựng (Industry Transformation): Đem lại giải pháp kỹ thuật vượt trội cho các tổng công ty tư vấn thiết kế cầu đường, các đơn vị thi công giàn khoan dầu khí ngoài khơi (Vietsovpetro), giúp tối ưu hóa khối lượng kết cấu thép, giảm 15% - 25% chi phí gia cố vật liệu so với phương pháp tăng cứng truyền thống.
  • Ảnh hưởng An ninh - Quốc phòng (Defense & Strategic Influence): Ứng dụng trực tiếp vào việc nâng cấp, kéo dài tuổi thọ công tác của hệ thống Nhà giàn DK1, các ngọn hải đăng tiền tiêu (như Hải đăng đảo Đá Lát - quần đảo Trường Sa), công trình công sự phòng thủ ven biển chịu tác động liên tục của bão biển cấp 12 và sóng xung kích.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận khung thuật toán nguồn mở kết hợp giữa giải tích phần tử hữu hạn và trí tuệ nhân tạo, mở rộng cho các bài toán cơ học vật rắn biến dạng phức tạp.
  • Kỹ sư trưởng & Chuyên gia kết cấu R&D: Sở hữu công cụ tính toán và bảng tra tham số tối ưu TMD tức thời, hỗ trợ thiết kế kháng chấn cho nhà siêu cao tầng, cầu dây văng nhịp lớn và tháp truyền hình.
  • Cơ quan quản lý công trình biển và quốc phòng: Nắm vững luận cứ khoa học để phê duyệt các phương án sửa chữa, gia cường kết cấu công trình biển đảo mà không làm gián đoạn nhiệm vụ trực chiến sẵn sàng chiến đấu.

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ì?

Đóng góp độc đáo nhất là việc mở rộng lý thuyết tối ưu hóa dao động cổ điển của Den Hartog (1956) từ hệ 1 bậc tự do tiền định sang hệ khung không gian 3D nhiều bậc tự do (M-DOF) tổng quát chịu tải trọng ngẫu nhiên dừng, đồng thời thiết lập biểu thức giải tích ma trận hàm truyền $H(\omega)$ tích hợp đầy đủ tương tác cục bộ của TMD tại nút bất kỳ.

2. Điểm 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 với các nghiên cứu của Venanzi (2013) hay Krenk & Høgsberg (2014) vốn phụ thuộc vào mô phỏng Monte Carlo lặp hàng ngàn lần rất tốn thời gian, luận án đã tiên phong ứng dụng Mạng Nơ-ron nhân tạo (ANN) đa tầng đóng vai trò bộ giải thay thế (Surrogate Model), cho phép xác định chính xác thông số tối ưu $(m_T, k_T, c_T)$ chỉ trong thời gian dưới 1 mili-giây.

3. Phát hiện bất ngờ nhất từ số liệu khảo sát là gì?

Đó là sự tồn tại của "hiệu ứng suy giảm triệt để" khi tăng hệ số cản $c_T$ vượt ngưỡng tối ưu. Khác với trực giác kỹ thuật cho rằng tăng cản luôn làm tăng khả năng dập tắt dao động, số liệu phổ đáp ứng $S_{Wc}$ chứng minh khi $c_T$ quá lớn, TMD bị "khóa cứng" vào kết cấu chính, khiến hệ phản ứng y hệt một kết cấu không có giảm chấn nhưng bị tăng thêm tải trọng bản thân.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp chi tiết toàn bộ thuật toán lập trình phần tử hữu hạn, sơ đồ khối logic, cấu trúc dữ liệu ma trận, kiến trúc mạng nơ-ron (số lớp ẩn, số nơ-ron, hàm kích hoạt) và phụ lục chương trình tính toán số minh bạch, cho phép các viện nghiên cứu độc lập tái lập toàn bộ kết quả khảo sát.

5. Tầm nhìn phát triển nghiên cứu 10 năm tới được vạch ra như thế nào?

Chiến lược 10 năm định hướng phát triển hệ thống điều khiển dao động thông minh thích ứng (Smart Adaptive Structural Control), tích hợp các cảm biến IoT đo dao động thời gian thực với chip biên AI chạy mô hình ANN của luận án để điều khiển các bộ giảm chấn bán tích cực MR Damper thích nghi ngay lập tức với các thảm họa thiên tai bất thường.

Kết luận

  1. Luận án đã thiết lập thành công mô hình toán học và thuật toán phần tử hữu hạn tổng quát để phân tích dao động ngẫu nhiên của hệ dầm và khung không gian 3D có lắp thiết bị tiêu tán năng lượng TMD.
  2. Xây dựng hoàn chỉnh chương trình tính toán số trong miền tần số, xác định chính xác ma trận mật độ phổ đáp ứng chuyển vị, vận tốc và gia tốc dưới kích động ồn trắng và ồn màu.
  3. Khảo sát toàn diện ảnh hưởng của các tham số then chốt ($m_T, k_T, c_T$, vị trí đặt), khẳng định bố trí đa TMD phân tán (d-MTMD) tại các vị trí $a/4$ và $a/2$ mang lại hiệu quả dập dao động vượt trội từ 30% đến trên 65%.
  4. Tiên phong ứng dụng mạng nơ-ron nhân tạo (ANN - MLP) dự báo chuẩn xác tần số dao động riêng và tối ưu hóa tham số TMD tức thời với sai số cực nhỏ dưới 1.2%.
  5. Cung cấp luận cứ khoa học và giải pháp công nghệ có giá trị ứng dụng cao cho việc thiết kế mới và gia cường các công trình an ninh quốc phòng (nhà giàn DK1, hải đăng Trường Sa), cầu dây văng và nhà siêu cao tầng tại Việt Nam.