Tổng quan về luận án

Nghiên cứu động lực học công trình và cơ học vật rắn hiện đại thường xuyên phải đối mặt với các bài toán dao động phi tuyến chịu tác động của tải trọng ngẫu nhiên như động đất, sóng biển, gió bão và rung động giao thông. Luận án tiến sĩ chuyên ngành Cơ học vật rắn (mã số 62440107) với đề tài "Tiêu chuẩn đối ngẫu trong phương pháp tuyến tính hóa tương đương cho hệ phi tuyến nhiều bậc tự do chịu kích động ngẫu nhiên" của nghiên cứu sinh Nguyễn Như Hiếu (Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, dưới sự hướng dẫn khoa học của GS. Nguyễn Đông Anh và PGS. Ninh Quang Hải) đại diện cho một bước đột phá lý thuyết trong việc giải quyết các giới hạn cố hữu của cơ học ngẫu nhiên phi tuyến.

Vấn đề cốt lõi (research gap) mà luận án giải quyết xuất phát từ hạn chế căn bản của phương pháp tuyến tính hóa tương đương (Equivalent Linearization - EQL) kinh điển bắt nguồn từ các nghiên cứu tiên phong của Booton (1954), Kazakov (1954) và Caughey (1956). Phương pháp truyền thống dựa trên tiêu chuẩn sai số bình phương trung bình chỉ mang lại kết quả xấp xỉ chấp nhận được đối với các hệ phi tuyến yếu ($O(\varepsilon)$ với $\varepsilon \ll 1$). Khi hệ có tính phi tuyến mạnh ($O(1)$) hoặc có sự tương tác phức tạp giữa nhiều bậc tự do (MDOF), sai số tính toán của phương pháp kinh điển tăng vọt, dẫn đến những dự báo sai lệch về độ an toàn và độ tin cậy kết cấu.

Luận án tập trung giải quyết 3 câu hỏi nghiên cứu và kiểm chứng 3 giả thuyết khoa học:

  1. RQ1: Làm thế nào để xây dựng một tiêu chuẩn tuyến tính hóa mới có khả năng bù trừ sai số khi thay thế phần tử phi tuyến mạnh bằng phần tử tuyến tính tương đương trong hệ nhiều bậc tự do?
    • Giả thuyết H1: Việc thiết lập bước thay thế ngược (đối ngẫu) từ hệ tuyến tính hóa về một hệ phi tuyến cùng lớp sẽ tạo ra cơ chế tự cân bằng sai số thống kê.
  2. RQ2: Hệ phương trình xác định các ma trận hệ số tuyến tính hóa theo tiêu chuẩn mới có thể đóng kín giải tích và giải được đối với hệ nhiều bậc tự do hay không?
    • Giả thuyết H2: Kết hợp gradient ngẫu nhiên của hàm phi tuyến với kỹ thuật ma trận hàm đáp ứng tần số và mật độ phổ cho phép thiết lập hệ phương trình ma trận đóng kín khả nghịch.
  3. RQ3: Tiêu chuẩn đối ngẫu có thể áp dụng chính xác cho các hệ liên tục phức tạp như dầm Euler-Bernoulli chịu kích động ngẫu nhiên dải rộng và dầm mang khối lượng tập trung trên nền đàn hồi hay không?
    • Giả thuyết H3: Biểu diễn biên độ mode dao động thông qua tiêu chuẩn đối ngẫu sẽ giảm thiểu sai số mô men bậc hai xuống dưới mức sai số của các phương pháp năng lượng hiện hành.

Khung lý thuyết của luận án được định vị trên nền tảng cơ học giải tích ngẫu nhiên, lý thuyết quá trình ngẫu nhiên Markov/Wiener và phương trình đạo hàm riêng parabolic Fokker-Planck-Kolmogorov (FPK). Đóng góp mang tính định lượng của công trình thể hiện ở việc giảm thiểu sai số đáp ứng bình phương trung bình $\mathbb{E}[x_1^2], \mathbb{E}[x_2^2]$ từ mức 20%–45% của phương pháp kinh điển xuống dưới 5%–10% so với nghiệm chính xác FPK và mô phỏng Monte-Carlo, ngay cả khi tham số phi tuyến $\beta_1, \beta_2$ tăng mạnh. Phạm vi nghiên cứu bao quát các hệ rời rạc 2-DOF phi tuyến bậc ba, bậc năm, cản phi tuyến kết hợp độ cứng phi tuyến, và hệ liên tục dầm Euler-Bernoulli mở rộng đến 5 mode dao động ($M=5$).

Literature Review và Positioning

Lịch sử phát triển của phương pháp tuyến tính hóa ngẫu nhiên được định hình qua nhiều thập kỷ với các dòng nghiên cứu chủ đạo:

Dòng nghiên cứu cổ điển khởi đầu từ Booton (1954), Kazakov (1954) và Caughey (1956, 1963). Caughey đã chính thức hóa khái niệm tuyến tính hóa tương đương cho các hệ dao động ngẫu nhiên kích động ồn trắng Gaussian dựa trên tiêu chuẩn cực tiểu hóa sai số bình phương trung bình $\mathbb{E}[e^T e] \to \min$. Atalik và Utku (1976) đã tổng quát hóa phương pháp cho hệ nhiều bậc tự do bằng việc áp dụng định lý đạo hàm kỳ vọng cho các biến ngẫu nhiên Gaussian: $\mathbb{E}[z \Phi^T(z)] = \mathbb{E}[z z^T] \mathbb{E}[\nabla \Phi^T(z)]$, đặt nền móng tính toán cho ma trận hệ số tuyến tính hóa $C_e, K_e$. Tuy nhiên, Roberts và Spanos (1990) cũng như Socha (2008) đã chỉ ra rằng điểm yếu chí mạng của dòng phương pháp này là việc đánh mất các đặc trưng phi tuyến mạnh do phép xấp xỉ một chiều mang tính cục bộ.

Dòng nghiên cứu tiếp cận theo tiêu chuẩn năng lượng xuất hiện nhằm khắc phục khiếm khuyết trên. Zhang et al. (1990) đề xuất tiêu chuẩn cực tiểu hóa bình phương trung bình thế năng $\mathbb{E}[(V(q) - V_e(q))^2] \to \min$, trong khi Elishakoff và Colombi (1993) mở rộng tiêu chuẩn dựa trên hàm hao tán năng lượng $\mathbb{E}[(s(\dot{x})[W(\dot{x}) - W_e(\dot{x})])^2] \to \min$. Song song đó, các nghiên cứu về tuyến tính hóa phi tham số (Nonparametric Linearization) của Fujimura và Der Kiureghian (2007), Luca và Der Kiureghian (2010), Reza và Mohsen (2014, 2016), Broccardo và Armen (2015) đã nỗ lực khớp các điểm thiết kế trong không gian biến chuẩn, nhưng phải trả giá bằng chi phí tính toán số rất cao.

Cuộc tranh luận khoa học sâu sắc tồn tại giữa hai trường phái:

  • Trường phái thứ nhất (Tiêu chuẩn năng lượng/hình học): Cho rằng cần bảo toàn các đại lượng vô hướng toàn cục như thế năng $V(x)$ hoặc năng lượng hao tán để duy trì tính ổn định của đáp ứng biên độ lớn (Zhang et al., 1990; Elishakoff & Colombi, 1993).
  • Trường phái thứ hai (Tiêu chuẩn điều chỉnh sai số vi phân): Cho rằng năng lượng không phản ánh đầy đủ phân bố không gian và tính ghép kênh giữa các bậc tự do, do đó cần điều chỉnh trực tiếp phương trình chuyển động thông qua các tham số phụ trợ (Nguyễn Đông Anh & Di Paola, 1995; Elishakoff et al., 2009).

Luận án của Nguyễn Như Hiếu định vị tại giao điểm của việc giải quyết bài toán phi tuyến mạnh cho hệ MDOF bằng cách mở rộng ý tưởng "tiêu chuẩn đối ngẫu" vốn mới chỉ được nghiên cứu sơ khởi trên hệ 1-DOF (Nguyễn Đông Anh, 2010; Anh et al., 2012a, 2012b). So với hai công trình quốc tế tiêu biểu là phương pháp tiệm cận phổ của Spanos và Agathoklis (2013) trên mô hình trễ Bouc-Wen và kỹ thuật hỗn hợp Gaussian của Wang và Song (2016), phương pháp tiêu chuẩn đối ngẫu của luận án vượt trội về sự tinh gọn trong biểu thức giải tích đóng kín, không đòi hỏi thuật toán tối ưu hóa lặp phức tạp mà vẫn đạt độ chính xác tương đương nghiệm FPK chính xác.

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

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

Luận án tạo ra bước tiến quan trọng mở rộng lý thuyết dao động ngẫu nhiên và cơ học kết cấu phi tuyến:

Mở rộng lý thuyết tuyến tính hóa tương đương của Caughey – Atalik – Utku: Thay vì chỉ thực hiện phép thay thế thuận một chiều từ phương trình phi tuyến gốc: $$M\ddot{x} + C\dot{x} + Kx + \Phi(x, \dot{x}) = Q(t)$$ sang hệ tuyến tính hóa tương đương: $$M\ddot{q} + (C + C_e)\dot{q} + (K + K_e)q = Q(t)$$ với sai số vi phân thuận $e = \Phi(q, \dot{q}) - C_e \dot{q} - K_e q$, luận án đề xuất bước thay thế đối ngẫu thứ hai. Trong bước này, tổ hợp tuyến tính $C_e \dot{q} + K_e q$ lại được xấp xỉ ngược về dạng phi tuyến $D\Phi(q, \dot{q})$, tạo ra sai số đối ngẫu: $$\psi = C_e \dot{q} + K_e q - D\Phi(q, \dot{q})$$ trong đó $D = [d_{ij}]_{n \times n}$ là ma trận hệ số thay thế trở lại.

Xác lập hàm mục tiêu tối ưu hóa đối ngẫu: Phiếm hàm sai số tổng quát được thiết lập dưới dạng: $$\mathbb{E}[e^T e] + \alpha \mathbb{E}[\psi^T \psi] \to \min_{C_e, K_e, D}$$ với $\alpha > 0$ là tham số điều chỉnh trọng số đối ngẫu (chọn $\alpha = 1$ cho nghiệm cân bằng chuẩn).

+-------------------------------------------------------------------------------+
|                       KHUNG PHÂN TÍCH ĐỐI NGẪU CHO HỆ MDOF                    |
+-------------------------------------------------------------------------------+
|  Hệ phi tuyến gốc:          Phi(q, dq)                                        |
|                                |                                              |
|  [Thay thế thuận]             | Sai số thuận: e = Phi - Ce*dq - Ke*q         |
|                                v                                              |
|  Hệ tuyến tính hóa:         Ce*dq + Ke*q                                      |
|                                |                                              |
|  [Thay thế đối ngẫu]          | Sai số đối ngẫu: psi = Ce*dq + Ke*q - D*Phi  |
|                                v                                              |
|  Tái cấu trúc phi tuyến:    D * Phi(q, dq)                                    |
+-------------------------------------------------------------------------------+
|  TIÊU CHUẨN ĐỐI NGẪU: E[e^T * e] + alpha * E[psi^T * psi] --> min_{Ce, Ke, D} |
+-------------------------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp nhuần nhuyễn 3 trụ cột lý thuyết: (1) Lý thuyết phương trình vi phân ngẫu nhiên Itô/Stratonovich, (2) Lý thuyết phổ và biến đổi Fourier trong miền tần số, (3) Lý thuyết dao động mode đàn hồi Euler-Bernoulli.

Hệ thống phương trình ma trận xác định hệ số tuyến tính hóa đối ngẫu được thiết lập từ điều kiện triệt tiêu đạo hàm riêng của phiếm hàm mục tiêu $\frac{\partial \varepsilon_i}{\partial c_{ije}} = 0, \frac{\partial \varepsilon_i}{\partial k_{ije}} = 0, \frac{\partial \varepsilon_i}{\partial d_{ij}} = 0$:

$$\begin{bmatrix} (1+\alpha)\mathbb{E}[\dot{q}\dot{q}^T] & (1+\alpha)\mathbb{E}[\dot{q}q^T] & -\alpha\mathbb{E}[\dot{q}\Phi^T] \ (1+\alpha)\mathbb{E}[q\dot{q}^T] & (1+\alpha)\mathbb{E}[qq^T] & -\alpha\mathbb{E}[q\Phi^T] \ -\alpha\mathbb{E}[\Phi\dot{q}^T] & -\alpha\mathbb{E}[\Phi q^T] & \alpha\mathbb{E}[\Phi\Phi^T] \end{bmatrix} \begin{bmatrix} (C_i^e)^T \ (K_i^e)^T \ (D_i)^T \end{bmatrix} = \begin{bmatrix} \mathbb{E}[\dot{q}\Phi_i] \ \mathbb{E}[q\Phi_i] \ 0 \end{bmatrix}$$

Khung phân tích xác định rõ các điều kiện biên (boundary conditions): các phần tử phi tuyến $\Phi(x, \dot{x})$ phải là các hàm liên tục, khả vi từng đoạn, kích động đầu vào $Q(t)$ là quá trình ngẫu nhiên dừng Gaussian có kỳ vọng triệt tiêu $\mathbb{E}[Q(t)] = 0$, và ma trận tương quan $\mathbb{E}[qq^T]$ không suy biến.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo quan điểm bản thể luận hiện thực phản biện (Critical Realism) và nhận thức luận thực chứng định lượng suy diễn (Quantitative Positivist Paradigm). Thiết kế nghiên cứu kết hợp chặt chẽ giữa 3 phân hệ phương pháp luận:

  1. Phân hệ Giải tích lý thuyết (Analytical formulation): Khai triển giải tích ma trận đáp ứng tần số $H(\omega)$ và tích phân phổ Duhamel để đóng kín các hệ số phi tuyến.
  2. Phân hệ Nghiệm chính xác Fokker-Planck-Kolmogorov (FPK Benchmark): Thiết lập phương trình đạo hàm riêng mô tả hàm mật độ xác suất dừng $f(x, \dot{x})$ với thế năng dừng triệt tiêu dòng xác suất $G(x) = 0$ để làm chuẩn tuyệt đối (ground truth).
  3. Phân hệ Mô phỏng số Monte-Carlo (Numerical MCS): Tạo lập 10.000 đến 100.000 mẫu quỹ đạo kích động ồn trắng Gaussian để kiểm chứng các hệ không tồn tại nghiệm FPK chính xác (như hệ cản phi tuyến kết hợp hoặc dầm nhiều mode).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ tính chuẩn mực học thuật cao:

  • Xác định ma trận đáp ứng tần số: $$H(\omega) = [-\omega^2 M + i\omega(C + C_e) + (K + K_e)]^{-1}$$
  • Thiết lập tích phân ma trận mật độ phổ: Biểu diễn các mô men bậc hai của chuyển vị và vận tốc thông qua ma trận mật độ phổ kích động $S_Q(\omega)$: $$\mathbb{E}[qq^T] = \int_{-\infty}^{\infty} H(\omega) S_Q(\omega) H^{*T}(\omega) d\omega$$ $$\mathbb{E}[\dot{q}\dot{q}^T] = \int_{-\infty}^{\infty} \omega^2 H(\omega) S_Q(\omega) H^{*T}(\omega) d\omega$$
  • Thiết lập chu trình lặp tự phối (Self-consistent Iteration): Giải hệ phương trình phi tuyến đóng kín kết hợp giữa các phần tử $k_{ije}, c_{ije}, d_{ij}$ và các mô men $\mathbb{E}[q_i q_j], \mathbb{E}[\dot{q}_i \dot{q}_j]$ bằng thuật toán giải tích tích phân thặng dư Cauchy và thuật toán lặp hội tụ cao với sai số dung sai $10^{-6}$.
+-------------------------------------------------------------------------------+
|                       SƠ ĐỒ THUẬT TOÁN ĐÓNG KÍN VÀ KIỂM CHỨNG                 |
+-------------------------------------------------------------------------------+
|  1. Thiết lập ma trận Ke, Ce, D từ Tiêu chuẩn đối ngẫu                        |
|                                |                                              |
|                                v                                              |
|  2. Tính Ma trận hàm truyền H(omega) = [-w^2*M + i*w*(C+Ce) + (K+Ke)]^(-1)   |
|                                |                                              |
|                                v                                              |
|  3. Tính mô men bậc hai: E[q*q^T] = Tích phân [ H(w) * S_Q(w) * H*(w)^T dw ] |
|                                |                                              |
|                                v                                              |
|  4. Lặp hội tụ hệ đóng kín (Dung sai sai số < 10^(-6))                        |
|                                |                                              |
|                                v                                              |
|  5. So sánh đối chứng chuẩn:                                                  |
|     - Nghiệm chính xác FPK (Hệ thế năng dừng)                                 |
|     - Mô phỏng Monte-Carlo MCS (10^4 - 10^5 samples)                          |
|     - Tuyến tính hóa kinh điển (TT) & Tuyến tính hóa năng lượng (NL)          |
+-------------------------------------------------------------------------------+

Data và phân tích

Luận án tiến hành kiểm thử số trị trên các hệ kết cấu đa dạng với bộ tham số chuẩn hóa:

  • Hệ 2-DOF phi tuyến bậc ba: $m_1=m_2=1, c_1=c_2=0.2, k_1=k_2=1$, lực phi tuyến $\Phi_1 = \beta_1 x_1^3 + \beta_2 (x_1 - x_2)^3, \Phi_2 = \beta_2 (x_2 - x_1)^3$, khảo sát tham số phi tuyến $\beta_1, \beta_2$ biến thiên từ $0$ đến $5.0$.
  • Hệ dầm Euler-Bernoulli liên tục tựa giản đơn: Chiều dài dầm $L$, chiều dày $h$, độ mảnh kết cấu $L/h \in [50, 400]$, bán kính quán tính $R$, hệ số cản nhớt $c$, hệ số nền đàn hồi $\gamma \in [0, 50]$, mật độ phổ kích động $S_0$.
  • Hệ dầm mang khối lượng tập trung: Khối lượng tập trung $M_0$ đặt tại tọa độ $d_0/L \in [0.1, 0.9]$, tỉ số khối lượng $M_0/m_0$ khảo sát từ $1.0$ đến $10.0$.
  • Công cụ tính toán: Xây dựng mã nguồn giải tích và mô phỏng số chuyên sâu trên môi trường MATLAB/Mathematica, thực hiện kiểm tra độ nhạy (sensitivity analysis) với các bước lưới tích phân thích ứng.

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

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

  1. Khả năng triệt tiêu sai số phi tuyến vượt bậc của Tiêu chuẩn đối ngẫu (DN): Trong hệ 2-DOF có độ cứng phi tuyến bậc ba ($\beta_1 = 1, \beta_2 = 2.0$), phương pháp đối ngẫu (DN) mang lại đáp ứng $\mathbb{E}[x_1^2] = 0.4418$, gần như trùng khớp hoàn hảo với nghiệm chính xác FPK ($\mathbb{E}[x_1^2]_{CX} = 0.4491$, sai số chỉ $-1.62%$). Trong khi đó, phương pháp thông thường (TT) cho giá trị $0.5186$ (sai số $+15.48%$) và phương pháp năng lượng (NL) cho giá trị $0.3872$ (sai số $-13.78%$). Khi tính phi tuyến tăng mạnh lên $\beta_1 = 5.0$, sai số của DN vẫn được khống chế dưới $3.5%$, chứng minh tính bền vững của ma trận đối ngẫu $D$.

  2. Quy luật chi phối tuyệt đối của Mode thứ nhất trong đáp ứng dầm liên tục: Khi phân tích dao động dầm tựa giản đơn mở rộng đến 5 mode dao động ($M=5$), tỷ số đáp ứng bình phương trung bình giữa mode 1 và các mode bậc cao là áp đảo (Bảng 4.11 từ nguồn văn bản): tỷ số $\mathbb{E}[w_1^2]/\mathbb{E}[w_2^2] = 81.0$, $\mathbb{E}[w_1^2]/\mathbb{E}[w_3^3] = 729.0$, $\mathbb{E}[w_1^2]/\mathbb{E}[w_4^2] = 4096.0$, và $\mathbb{E}[w_1^2]/\mathbb{E}[w_5^2] = 15625.0$. Phát hiện này khẳng định tính khả thi của việc sử dụng xấp xỉ đơn mode hoặc hai mode trong tính toán thiết kế kỹ thuật thực hành mà không làm mất đi độ chính xác cơ học cốt lõi.

  3. Hiện tượng phi tuyến hóa đáp ứng phụ thuộc vào độ mảnh $L/h$ và nền đàn hồi: Với dầm có độ mảnh cao ($L/h = 200$), tỷ số đáp ứng bình phương trung bình giữa hệ phi tuyến và hệ tuyến tính giảm mạnh xuống xấp xỉ $0.25 - 0.30$. Tác dụng gia cường của nền đàn hồi $\gamma$ làm tăng độ cứng tổng thể của hệ, triệt tiêu đáng kể biên độ dao động ngẫu nhiên, đưa hàm mật độ xác suất $f(w_1)$ co cụm nhanh quanh vị trí cân bằng.

  4. Hiệu ứng vị trí và tỷ số khối lượng tập trung ($M_0/m_0, d_0/L$): Khi khối lượng tập trung dịch chuyển về phía giữa nhịp ($d_0/L \to 0.5$) hoặc tỷ số $M_0/m_0$ tăng từ 1 lên 5, hiệu ứng quán tính cục bộ làm thay đổi cấu trúc mode dao động, khiến sự chênh lệch giữa phương pháp kinh điển và phương pháp đối ngẫu càng giãn rộng, khẳng định tính bắt buộc phải dùng tiêu chuẩn đối ngẫu khi tính toán các kết cấu dầm cầu có phương tiện dừng đỗ.

Implications đa chiều

  • Về mặt lý thuyết: Hoàn thiện lý thuyết tuyến tính hóa ngẫu nhiên thông qua việc cung cấp cấu trúc đại số đối ngẫu 3 thành phần ($C_e, K_e, D$), giải quyết bài toán bế tắc kéo dài hơn 6 thập kỷ về tính phi tuyến mạnh của Caughey (1956).
  • Về mặt phương pháp luận: Mở ra phương pháp tiếp cận đóng kín cho các bài toán tương tác phức tạp (kết cấu - nền đàn hồi, kết cấu mang khối lượng di động) mà không cần phụ thuộc vào mô phỏng Monte-Carlo tốn kém thời gian.
  • Về ứng dụng thực tiễn: Ứng dụng trực tiếp trong thiết kế kết cấu công trình biển chịu sóng ngẫu nhiên, nhà cao tầng chịu bão, dầm cầu chịu tải trọng động đất, và công nghệ điều khiển giảm chấn bán tích cực (TMD phi tuyến, thiết bị cách ly đáy).
  • Về khuyến nghị chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học để cập nhật các tiêu chuẩn thiết kế kháng chấn công trình (TCVN, Eurocode 8, ASCE 7), cho phép giảm hệ số an toàn dư thừa do sai số tính toán của các mô hình tuyến tính hóa cũ.

Limitations và Future Research

Mặc dù đạt được những đột phá vượt bậc, luận án cũng ghi nhận các giới hạn nghiên cứu mang tính biên:

  1. Giới hạn về dạng kích động: Kích động đầu vào chủ yếu được giả định là quá trình ngẫu nhiên dừng ồn trắng Gaussian hoặc ồn màu dải rộng. Các kích động non-Gaussian xung kích (Poissonian pulse) hoặc kích động không dừng (non-stationary earthquake records) chưa được bao quát hoàn toàn.
  2. Giới hạn về số bậc tự do giải tích: Các công thức nghiệm giải tích tường minh chi tiết được triển khai tối ưu cho hệ 2-DOF và hệ liên tục dầm quy về $M \le 5$ mode. Với các hệ siêu nhiều bậc tự do ($n > 100$ trong phần tử hữu hạn FEM), hệ phương trình ma trận đòi hỏi các thuật toán xử lý ma trận thưa.
  3. Đặc trưng hàm phi tuyến: Chưa giải quyết triệt để các hệ có tính chất phi tuyến gián đoạn mạnh như va đập (piecewise linear/impact) hoặc hệ có hiện tượng trễ suy giảm độ cứng và trễ trượt phức tạp (degrading hysteretic systems).

Định hướng nghiên cứu tương lai (Future Agenda):

  • Mở rộng tiêu chuẩn đối ngẫu cho hệ phi tuyến chịu đồng thời kích động tham số ngẫu nhiên (parametric excitation) và kích động ngoài ngẫu nhiên.
  • Phát triển thuật toán tuyến tính hóa đối ngẫu tích hợp vào phần mềm phần tử hữu hạn (FEM/OpenSees/ANSYS) để phân tích kết cấu không gian phức tạp.
  • Khảo sát tối ưu hóa tham số điều chỉnh $\alpha(\beta)$ như một hàm phụ thuộc vào mức độ phi tuyến thay vì cố định $\alpha = 1$.
  • Ứng dụng tiêu chuẩn đối ngẫu trong đánh giá độ tin cậy kết cấu phụ thuộc thời gian (Time-variant Structural Reliability) và bài toán vượt ngưỡng an toàn (First-passage failure).

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

  • Ảnh hưởng học thuật: Luận án đóng góp trực tiếp vào các chương trình nghiên cứu cơ bản quốc gia (tài trợ bởi Quỹ NAFOSTED) do GS. Nguyễn Đông Anh chủ trì, với các công bố quốc tế uy tín, tạo tiềm năng trích dẫn cao trong cộng đồng cơ học kết cấu ngẫu nhiên quốc tế (Journal of Sound and Vibration, Probabilistic Engineering Mechanics, Nonlinear Dynamics).
  • Chuyển đổi công nghiệp: Cung cấp công cụ tính toán chính xác cho các viện thiết kế xây dựng, dầu khí, giao thông trong việc tối ưu hóa kích thước dầm, trụ cầu cảng và tháp gió ngoài khơi, giúp tiết kiệm vật liệu xây dựng từ 5%–12% nhờ loại bỏ sai số ước tính biên độ dao động cực đại.
  • Lợi ích xã hội: Nâng cao độ an toàn của các công trình công cộng trọng điểm trước thiên tai ngẫu nhiên (động đất cấp hiểm họa, bão lớn), giảm thiểu rủi ro sự cố công trình và bảo vệ tính mạng con ngườ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/Xây dựng/Toán ứng dụng: Tiếp cận một phương pháp luận giải tích ngẫu nhiên chuẩn mực, cung cấp nền tảng giải quyết các đề tài luận án về hệ động lực phi tuyến phức tạp.
  • Các nhà khoa học cao cấp và giảng viên đại học: Bổ sung một lý thuyết cải tiến xuất sắc vào giáo trình giảng dạy chuyên đề Dao động ngẫu nhiên và Động lực học công trình nâng cao.
  • Kỹ sư R&D và chuyên gia kết cấu công trình: Sở hữu công thức xấp xỉ nhanh, chính xác để tính toán đáp ứng dao động và thiết kế hệ thống giảm chấn cho kết cấu dầm, khung, tấm vỏ.
  • Các cơ quan biên soạn tiêu chuẩn quy chuẩn xây dựng: Có dữ liệu khoa học tin cậy để hiệu chỉnh các công thức tính toán tải trọng động ngẫu nhiên trong tiêu chuẩn quốc gia.

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à nó mở rộng lý thuyết nào? Trả lời: Đóng góp độc đáo nhất là việc thiết lập Tiêu chuẩn đối ngẫu (Dual Criterion) thông qua việc tích hợp bước thay thế ngược với ma trận $D$ và hàm $D\Phi(q, \dot{q})$. Công trình đã mở rộng trực tiếp lý thuyết tuyến tính hóa tương đương kinh điển của Caughey (1956) và Atalik - Utku (1976), chuyển đổi cấu trúc tối ưu hóa đơn lẻ thành cấu trúc tối ưu hóa bù trừ sai số kép, giải quyết trọn vẹn điểm nghẽn phi tuyến mạnh trong hệ MDOF.

  2. Điểm 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ế tiền nhiệm? Trả lời:

    • So với Phương pháp tuyến tính hóa năng lượng của Zhang et al. (1990): Phương pháp đối ngẫu của luận án không bị giới hạn ở các hệ thuần thế năng vị trí mà áp dụng tổng quát cho cả hệ có cản phi tuyến $\Phi(x, \dot{x})$ và giữ được dạng phương trình vi phân vector cục bộ thay vì ép về một hàm thế năng vô hướng $V(x)$.
    • So với Phương pháp tuyến tính hóa phi tham số của Fujimura & Der Kiureghian (2007): Phương pháp đối ngẫu cung cấp hệ phương trình đóng kín giải tích thuần túy, loại bỏ hoàn toàn nhu cầu rời rạc hóa điểm thiết kế và tối ưu phi tuyến số đa biến, giúp giảm thời gian tính toán hàng trăm lần.

  3. Phát hiện gây bất ngờ nhất có số liệu thực nghiệm/mô phỏng chứng minh? Trả lời: Phát hiện bất ngờ nhất là tính chất tự cân bằng sai số của tham số $\alpha = 1$: trong khi phương pháp thông thường (TT) luôn ước lượng thừa (overestimate) phương sai đáp ứng $\mathbb{E}[x_1^2]$ lên $+15.48%$ và phương pháp năng lượng (NL) luôn ước lượng thiếu (underestimate) xuống $-13.78%$, thì tiêu chuẩn đối ngẫu (DN) tạo ra sự triệt tiêu tương hỗ giữa hai xu hướng sai số này, kéo sai số về mức lý tưởng $-1.62%$ tại $\beta_1 = 1, \beta_2 = 2.0$ (Bảng 4.6).

  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) rõ ràng không? Trả lời: Hoàn toàn có. Luận án cung cấp hệ thống công thức giải tích tường minh từ ma trận khối lượng $M$, ma trận cản tương đương $C_e$, ma trận độ cứng tương đương $K_e$, ma trận đáp ứng tần số $H(\omega)$, cho đến quy trình lấy tích phân phổ năng lượng xác định mô men $\mathbb{E}[q_i q_j]$. Bất kỳ nhà nghiên cứu nào cũng có thể lập trình tái lập chính xác 100% kết quả các bảng số liệu (từ Bảng 4.1 đến Bảng 4.14) trên các phần mềm toán học tiêu chuẩn.

  5. Chương trình nghiên cứu 10 năm (10-Year Research Agenda) được vạch ra như thế nào? Trả lời: Chương trình 10 năm tập trung vào 3 giai đoạn: (1) Mở rộng tiêu chuẩn đối ngẫu cho hệ kết cấu liên tục 2D/3D (tấm, vỏ thoái hóa phi tuyến von Kármán), (2) Phát triển lý thuyết đối ngẫu cho quá trình ngẫu nhiên non-Gaussian và quá trình phi dừng, (3) Tích hợp thuật toán đối ngẫu vào lõi bộ giải phần mềm tính toán kết cấu thương mại phục vụ chẩn đoán sức khỏe công trình theo thời gian thực (Real-time Structural Health Monitoring).

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Như Hiếu đã hoàn thành xuất sắc mục tiêu nghiên cứu với 5 đóng góp cốt lõi:

  1. Đề xuất thành công Tiêu chuẩn đối ngẫu mới trong phương pháp tuyến tính hóa tương đương cho hệ ngẫu nhiên nhiều bậc tự do phi tuyến mạnh.
  2. Xây dựng trọn vẹn hệ phương trình ma trận đóng kín xác định các ma trận hệ số tuyến tính hóa $C_e, K_e$ và ma trận thay thế ngược $D$.
  3. Chứng minh tính ưu việt vượt trội của phương pháp đối ngẫu thông qua việc giảm sai số đáp ứng bình phương trung bình xuống dưới 5% so với nghiệm chính xác Fokker-Planck và mô phỏng Monte-Carlo.
  4. Mở rộng thành công phương pháp cho bài toán dao động ngẫu nhiên của hệ liên tục dầm Euler-Bernoulli tựa giản đơn và dầm mang khối lượng tập trung trên nền đàn hồi.
  5. Xác lập quy luật hội tụ mode dao động của dầm liên tục, chứng minh vai trò chi phối tuyệt đối của mode thứ nhất trong tính toán thiết kế kỹ thuật thực tiễn.

Công trình tạo ra bước tiến mẫu mực về nhận thức luận cơ học tính toán ngẫu nhiên, mở ra 3 hướng nghiên cứu ứng dụng đột phá và xác lập vị thế tiên phong của cơ học giải tích Việt Nam trên trường quốc tế.