Tổng quan về luận án

Nghiên cứu động lực học và dao động ngẫu nhiên phi tuyến đóng vai trò cốt lõi trong kỹ thuật công trình, hàng hải, hàng không vũ trụ và cơ học ứng dụng. Khi các kết cấu kỹ thuật chịu tác động của các tải trọng ngẫu nhiên trong tự nhiên như sóng biển tác động lên tàu thủy và công trình biển, xung động gió giật lên nhà cao tầng, cầu dây văng, hay rung chấn động đất ngẫu nhiên, các đáp ứng động lực học của hệ thống thường thể hiện tính phi tuyến rõ rệt. Về bản chất, các quá trình ngẫu nhiên này không thể dự báo chính xác theo thời gian tiền định mà phải phân tích thông qua các đặc trưng thống kê và lý thuyết xác suất. Việc giải chính xác phương trình vi phân ngẫu nhiên phi tuyến thông qua phương trình Fokker-Planck-Kolmogorov (FPK) bị giới hạn nghiêm ngặt trong một số rất ít hệ một bậc tự do (SDOF) có tính chất đặc biệt. Do đó, phát triển các phương pháp giải tích xấp xỉ có độ chính xác cao và khả năng mở rộng cho hệ nhiều bậc tự do (MDOF) là đòi hỏi cấp thiết của cơ học ngẫu nhiên hiện đại.

Khoảng trống nghiên cứu then chốt (Research Gap) xuất phát từ hạn chế cố hữu của phương pháp tuyến tính hóa tương đương kinh điển (Classical Equivalent Linearization - CEL) do Caughey (1959, 1963) khởi xướng. CEL áp dụng tiêu chuẩn cực tiểu hóa sai số bình phương trung bình trên toàn bộ không gian biến xác suất $(-\infty, +\infty)$, dẫn đến sai số tăng vọt khi hệ số phi tuyến mạnh. Nhằm khắc phục điều này, tiêu chuẩn sai số bình phương trung bình địa phương (Local Mean Square Error Criterion - LOMSEC) được Anh và Di Paola (1993) đề xuất và Lưu Xuân Hùng (2002) phát triển, đã thu hẹp miền lấy tích phân vào khoảng tập trung năng lượng hữu hạn $[-r_x, +r_x]$. Tuy nhiên, LOMSEC tồn tại một nhược điểm lớn: tham số biên hữu hạn $r_x$ là một ẩn số chưa thể xác định giải tích một cách tự động và đóng kín (closed-form), buộc phải dựa vào kinh nghiệm hoặc tinh chỉnh bán thực nghiệm.

Luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (Mã số: 9 52 01 01) của nghiên cứu sinh Nguyễn Cao Thắng dưới sự hướng dẫn khoa học của GS. Nguyễn Đông Anh và TS. Lưu Xuân Hùng tại Học viện Khoa học và Công nghệ - Viện Hàn lâm Khoa học và Công nghệ Việt Nam đã giải quyết triệt để bài toán này bằng cách phát triển tiêu chuẩn sai số bình phương trung bình địa phương – tổng thể (Global-Local Mean Square Error Criterion - GLOMSEC).

Các câu hỏi nghiên cứu và giả thuyết khoa học được xác định cụ thể:

  • $RQ_1$: Làm thế nào để khép kín về mặt toán học việc xác định hệ số tuyến tính hóa địa phương của LOMSEC mà không cần giả định tùy tiện về miền tích phân hữu hạn?
  • $RQ_2$: Quy luật tích hợp đối ngẫu giữa trường trung bình địa phương và phân bố xác suất tổng thể được thiết lập như thế nào cho hệ dao động ngẫu nhiên phi tuyến nhiều bậc tự do (MDOF)?
  • $RQ_3$: Tiêu chuẩn mới đáp ứng như thế nào đối với các dạng kích động phức tạp như ồn màu (colored noise) và các hệ động lực học kỹ thuật thực tế (hệ Duffing, Van der Pol, dao động lắc ngang tàu thủy)?
  • $H_1$: Việc lấy trung bình thống kê tổng thể (global expectation) của các hệ số tuyến tính hóa địa phương (local coefficients) theo toàn bộ không gian trạng thái sẽ tự động tối ưu hóa trọng số năng lượng, loại bỏ hoàn toàn tính bất định của tham số $r_x$.
  • $H_2$: Tiêu chuẩn GLOMSEC duy trì độ chính xác tiệm cận nghiệm FPK chính xác và mô phỏng Monte Carlo (MCS) ngay cả trong miền phi tuyến mạnh, vượt trội hoàn toàn so với tiêu chuẩn kinh điển của Caughey.

Khung lý thuyết của luận án được xây dựng dựa trên sự giao thoa giữa Lý thuyết phương trình vi phân ngẫu nhiên (Stochastic Differential Equations), Quá trình khuếch tán Markov (Markov Diffusion Processes), Phương trình Fokker-Planck-Kolmogorov (FPK), và Lý thuyết xấp xỉ đối ngẫu trong cơ học phi tuyến (Dual Approach in Nonlinear Mechanics). Về mặt quy mô và ý nghĩa, luận án đã phân tích toàn diện trên 5 mô hình SDOF tiêu biểu, hệ 2-DOF ghép phi tuyến mạnh, và hệ chịu kích động ồn màu bậc hai; được kiểm chứng định lượng độc lập thông qua phương pháp mô phỏng Monte Carlo với hàng triệu mẫu dữ liệu số.

Literature Review và Positioning

Lịch sử phát triển của lý thuyết dao động ngẫu nhiên phi tuyến ghi nhận các nhánh phương pháp chính: phương pháp giải tích chính xác FPK, phương pháp nhiễu tham số bé (Perturbation Methods), phương pháp trung bình hóa ngẫu nhiên (Stochastic Averaging Method), phương pháp phi tuyến hóa tương đương (Equivalent Non-linearization), và phương pháp tuyến tính hóa tương đương (Equivalent Linearization - EEL).

Trong nhánh tuyến tính hóa tương đương, Caughey (1959, 1963) và sau đó là Atalik & Utku (1976) đã đặt nền móng với tiêu chuẩn cực tiểu hóa kỳ vọng sai số bình phương $E[e^2] \to \min$, chứng minh định lý gradient kỳ vọng trên véc-tơ Gaussian: $E[\mathbf{z}\mathbf{\Phi}^T(\mathbf{z})] = E[\mathbf{z}\mathbf{z}^T]E[\nabla\mathbf{\Phi}^T(\mathbf{z})]$. Tuy nhiên, Crandall (1963, 2000) cùng Elishakoff (1992, 2000) đã chỉ rõ sự suy giảm nghiêm trọng về độ chính xác của tiêu chuẩn kinh điển khi tham số phi tuyến tăng cao hoặc khi hệ xuất hiện tính chất phi đối xứng và tự kích.

Một cuộc tranh luận khoa học lớn trong nhiều thập kỷ xoay quanh hai trường phái tiếp cận:

  • Trường phái tiếp cận toàn cục (Global Approach): Tiêu biểu là Caughey (1959) và Atalik & Utku (1976), duy trì miền tích phân $(-\infty, +\infty)$ nhằm bảo toàn tính giải tích và tính tổng quát cho hệ nhiều chiều, nhưng đánh đổi bằng việc làm phẳng (smearing) sai số ở các vùng tập trung năng lượng cao.
  • Trường phái tiếp cận cải tiến năng lượng và biến đổi cục bộ (Energy & Local Truncation Approach): Elishakoff & Zhang (1992) đề xuất tiêu chuẩn sai số thế năng cực tiểu $E[(U(x) - \frac{1}{2}kx^2)^2] \to \min$, giải quyết tốt phi tuyến đàn hồi nhưng bất lực trước phi tuyến cản; Anh & Di Paola (1993, 1995) đề xuất tiêu chuẩn tuyến tính hóa có điều chỉnh bậc cao và tiêu chuẩn LOMSEC trên miền hữu hạn $[-r_x, +r_x]$, tập trung vào vùng mật độ xác suất cao nhất. Tuy nhiên, rào cản lớn nhất của LOMSEC là không thể khép kín nghiệm giải tích do tham số $r_x$ chưa có biểu thức xác định tường minh.

Gần đây, các nghiên cứu quốc tế đã mở rộng tuyến tính hóa tương đương theo các hướng mới: Canor et al. (2016) áp dụng tuyến tính hóa cấu trúc kích thước lớn, Silva-González et al. (2017) nghiên cứu hệ đàn dẻo chịu động đất, Jalali (2015) nghiên cứu hệ trễ phi tuyến, và Su et al. (2018) phát triển quy trình cho kích động không dừng. Tại Việt Nam, trường phái nghiên cứu của GS. Nguyễn Đông Anh đã phát triển quan điểm tiếp cận đối ngẫu (Dual Approach) bao gồm quá trình thay thế lượt đi và lượt về (Nguyễn Đông Anh 2000; Nguyễn Ngọc Linh 2010; Nguyễn Như Hiếu 2012; Nguyễn Minh Triết 2015; N. Việt 2012).

Luận án của Nguyễn Cao Thắng định vị chính xác tại điểm nghẽn của LOMSEC: kế thừa quan điểm đối ngẫu để thiết lập mối liên kết toán học hữu cơ giữa xấp xỉ địa phương và kỳ vọng tổng thể. Luận án mở rộng thành công GLOMSEC từ hệ SDOF của N. Việt (2012) sang hệ nhiều bậc tự do MDOF và hệ chịu kích động ồn màu dải hẹp phức tạp, xác lập một bước tiến vượt bậc so với các công trình quốc tế cùng thời kỳ.

Đó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 trực tiếp lý thuyết tuyến tính hóa tương đương của Caughey (1959), lý thuyết sai số địa phương của Anh & Di Paola (1993, 2002) và khung đối ngẫu của Nguyễn Đông Anh (2000, 2012) thông qua các đóng góp cụ thể:

  1. Khép kín giải tích tiêu chuẩn địa phương: Chứng minh rằng hệ số tuyến tính hóa thực chất là một biến ngẫu nhiên phụ thuộc vào miền tập trung đáp ứng địa phương $r$. Bằng cách lấy kỳ vọng toán tổng thể trên toàn bộ phân bố mật độ xác suất (PDF), hệ số tuyến tính hóa tổng thể $b_{GL}$ và $k_{GL}$ được xác lập hoàn toàn tự động mà không cần bất kỳ tham số tinh chỉnh thực nghiệm nào.
  2. Mô hình lý thuyết tổng quát cho hệ phi tuyến nhiều bậc tự do: Thiết lập hệ phương trình ma trận tuyến tính hóa cho hệ phương trình vi phân: $$\mathbf{M}\ddot{\mathbf{x}} + \mathbf{C}\dot{\mathbf{x}} + \mathbf{K}\mathbf{x} + \mathbf{\Phi}(\mathbf{x}, \dot{\mathbf{x}}, \ddot{\mathbf{x}}) = \mathbf{Q}(t)$$ thành hệ tuyến tính hóa tương đương: $$(\mathbf{M} + \mathbf{M}_e)\ddot{\mathbf{q}} + (\mathbf{C} + \mathbf{C}_e)\dot{\mathbf{q}} + (\mathbf{K} + \mathbf{K}_e)\mathbf{q} = \mathbf{Q}(t)$$ trong đó ma trận $\mathbf{M}_e, \mathbf{C}_e, \mathbf{K}_e$ được đóng kín thông qua toán tử tích phân địa phương - tổng thể.
  3. Định lý chuyển tiếp tiệm cận: Chứng minh về mặt giải tích rằng khi miền tích phân mở rộng dần ra vô hạn ($r \to \infty$), nghiệm của GLOMSEC hội tụ trơn về nghiệm tuyến tính hóa kinh điển Caughey, xác nhận CEL chỉ là một trường hợp suy biến đặc biệt của GLOMSEC.
+-------------------------------------------------------------+
|               Hệ Phi Tuyến Gốc (SDOF / MDOF)                |
|       M x'' + C x' + K x + Phi(x, x', x'') = Q(t)          |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
|    1. Phân Tích Miền Tập Trung Đáp Ứng Địa Phương (Local)   |
|         Miền tích phân hữu hạn [-rx, +rx], [-rv, +rv]        |
|    Cực tiểu sai số LOMSEC: E_L[e^2] -> b_L(r), k_L(r)       |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
|      2. Tiếp Cận Đối Ngẫu & Khép Kín Tổng Thể (Global)      |
|    Hệ số GLOMSEC = Kỳ vọng toàn cục của hệ số địa phương    |
|      b_GL = E_G[b_L(r)], k_GL = E_G[k_L(r)] qua PDF p(r)    |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
|           3. Hệ Tuyến Tính Hóa Tương Đương Đóng Kín         |
|         (M + Me) q'' + (C + Ce) q' + (K + Ke) q = Q(t)      |
+-------------------------------------------------------------+
                              |
                              v
+-------------------------------------------------------------+
|    4. Giải Thuật Toán Lặp Ma Trận & Xuất Moment Bậc Hai     |
|   Tính E[x^2], E[x'^2], E[x1 x2] -> Đạt độ chính xác cao     |
+-------------------------------------------------------------+

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

Khung phân tích của luận án tích hợp liền mạch ba trụ cột lý thuyết:

  • Lý thuyết vi phân ngẫu nhiên Itô và Fokker-Planck-Kolmogorov: Sử dụng hàm thế năng, hàm Hamilton $H(x, \dot{x}) = \frac{1}{2}\dot{x}^2 + \int_0^x g(u)du$ để xác định cấu trúc hàm mật độ xác suất chính xác $p(x, \dot{x}) = C \exp\left(-\frac{2}{\sigma^2}\int_0^H f(v)dv\right)$ làm chuẩn đối sánh.
  • Lý thuyết không gian xác suất Hilbert và chuẩn sai số Euclid: Đánh giá sai số trên không gian Hilbert $L^2$ của các biến ngẫu nhiên qua chuẩn $e^2 = \mathbf{e}^T\mathbf{e}$.
  • Khung biến đổi tương đương đối ngẫu: Xử lý bản chất xấp xỉ hai chiều của quá trình tuyến tính hóa, cân bằng giữa mất mát năng lượng do tuyến tính hóa và sự bù trừ năng lượng phân bố.

Điều kiện biên và phạm vi hiệu lực của khung phân tích: Áp dụng cho các quá trình ngẫu nhiên dừng (stationary) và ergodic; lực kích động chuẩn Gaussian ồn trắng hoặc ồn màu thu được từ bộ lọc tuyến tính bậc hai; các hàm phi tuyến đa thức hoặc hàm phi tuyến liên tục trơn của dịch chuyển và vận tốc.

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

Thiết kế nghiên cứu

Luận án tuân thủ nghiêm ngặt quan điểm nhận thức luận thực chứng định lượng (Positivist Quantitative Paradigm) trong Cơ học tính toán và Động lực học ngẫu nhiên. Thiết kế nghiên cứu kết hợp đa tầng phương pháp (Methodological Triangulation):

  1. Nghiên cứu giải tích thuần túy (Analytical Derivation): Thiết lập các biểu thức tường minh cho hệ số tuyến tính hóa địa phương - tổng thể.
  2. Nghiên cứu chuẩn FPK chính xác (Exact Benchmark Solution): Áp dụng cho các hệ SDOF có nghiệm FPK giải tích để làm thước đo sai số tuyệt đối.
  3. Nghiên cứu mô phỏng số thống kê (Numerical Monte Carlo Simulation - MCS): Đóng vai trò là phương pháp kiểm chứng độc lập cho các hệ phức tạp, hệ nhiều bậc tự do MDOF và hệ chịu ồn màu không có nghiệm FPK chính xác.

Thiết kế đa tầng được phân chia rõ rệt:

  • Tầng 1: Hệ dao động ngẫu nhiên SDOF phi tuyến cản bậc ba, hệ Van der Pol, hệ Duffing đàn hồi phi tuyến, hệ Duffing kết hợp cản phi tuyến, dao động lắc ngang tàu thủy chịu sóng biển.
  • Tầng 2: Hệ dao động ngẫu nhiên phi tuyến 2-DOF ghép liên hợp.
  • Tầng 3: Hệ dao động ngẫu nhiên chịu kích động ồn màu dải hẹp thông qua phương trình lọc vi phân bậc hai.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu số và giải tích được tiêu chuẩn hóa theo các bước chặt chẽ:

+-----------------------------------------------------------------------------------+
| BƯỚC 1: Xây dựng mô hình toán phi tuyến hệ SDOF/MDOF                              |
|         Thiết lập phương trình vi phân ngẫu nhiên: M x'' + C x' + K x + g = Q(t)  |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
| BƯỚC 2: Thiết lập toán tử sai số địa phương LOMSEC                                 |
|         Xác định e_L = g(x, x') - (b_L x' + k_L x) trên [-r_x, +r_x]              |
|         Cực tiểu hóa sai số bình phương trung bình: Min S_L -> b_L(r), k_L(r)     |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
| BƯỚC 3: Áp dụng toán tử kỳ vọng đối ngẫu GLOMSEC                                  |
|         Tính tích phân tổng thể: b_GL = \int b_L(r) p(r) dr                       |
|                                  k_GL = \int k_L(r) p(r) dr                       |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
| BƯỚC 4: Thiết lập thuật toán lặp Atalik-Utku cải tiến                             |
|         a) Khởi tạo ma trận hiệp phương sai ban đầu [Dx]^(0)                      |
|         b) Cập nhật các hệ số b_GL, k_GL                                          |
|         c) Giải hệ phương trình Lyapunov xác định mô men bậc hai [Dx]^(k+1)       |
|         d) Kiểm tra hội tụ: || [Dx]^(k+1) - [Dx]^(k) || < 10^(-6)                 |
+-----------------------------------------------------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
| BƯỚC 5: Kiểm chứng chéo & Đánh giá sai số                                         |
|         So sánh trực tiếp: Sai số tuyệt đối & Tương đối so với FPK / Monte Carlo  |
+-----------------------------------------------------------------------------------+

Giao thức mô phỏng Monte Carlo được thực hiện với độ tin cậy thống kê cao: tích phân số phương trình vi phân ngẫu nhiên bằng thuật toán Runge-Kutta bậc 4, số lượng mẫu ngẫu nhiên $N = 10^5 \div 10^6$ thể hiện (realizations), bước thời gian tích phân $\Delta t = 0.001 \div 0.01\text{ s}$, loại bỏ giai đoạn dao động chuyển tiếp ban đầu (transient phase) để đảm bảo quá trình đạt trạng thái dừng Ergodic ổn định.

Data và phân tích

Hệ thống tính toán giải tích số và mô phỏng được lập trình toàn diện trên môi trường MATLAB. Quá trình xử lý dữ liệu tập trung vào các đặc trưng thống kê bậc hai của đáp ứng:

  • Mô men bậc nhất: $\mu_x = E[x] = 0$ (quá trình dừng kỳ vọng không).
  • Mô men bậc hai dịch chuyển: $E[x^2] = \sigma_x^2 = \int_{-\infty}^{+\infty} x^2 p_1(x) dx$.
  • Mô men bậc hai vận tốc: $E[\dot{x}^2] = \sigma_{\dot{x}}^2 = \int_{-\infty}^{+\infty} \dot{x}^2 p_2(\dot{x}) d\dot{x}$.
  • Hiệp phương sai tương quan chéo: $D_{12} = E[(x_1 - \mu_1)(x_2 - \mu_2)]$.

Đối với kích động ồn màu dải hẹp $f(t)$, bộ lọc vi phân tuyến tính bậc hai được mở rộng vào không gian trạng thái: $$\ddot{f} + 2\alpha_f \dot{f} + \omega_f^2 f = \omega_f^2 w(t)$$ với mật độ phổ công suất xác định: $$S_f(\omega) = |H(i\omega)|^2 S_0 = \frac{\omega_f^4 S_0}{(\omega_f^2 - \omega^2)^2 + 4\alpha_f^2 \omega^2}$$ Các kiểm tra độ bền vững (Robustness Checks) được thực hiện bằng cách quét liên tục các tham số phi tuyến từ rất bé ($\gamma, \beta = 0.01$) đến miền phi tuyến cực mạnh ($\gamma, \beta = 5.0 \div 10.0$), đồng thời thay đổi cường độ kích động ồn trắng $S_0$ và tần số lọc $\omega_f$.

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

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

  1. Vượt trội hoàn toàn về độ chính xác so với phương pháp kinh điển (CEL): Trong tất cả các mô hình khảo sát, sai số tương đối của mô men bậc hai $E[x^2]$ tính theo GLOMSEC luôn nhỏ hơn đáng kể so với CEL khi so sánh với nghiệm chính xác FPK hoặc mô phỏng Monte Carlo. Đối với hệ Duffing có mức độ phi tuyến đàn hồi cao ($\beta = 1.0 \div 5.0$), sai số của CEL thường vượt quá $15% - 30%$, trong khi sai số của GLOMSEC được duy trì dưới mức $1% - 4%$.
Mô hình dao động ngẫu nhiên Tham số phi tuyến tiêu biểu Sai số CEL (%) Sai số GLOMSEC (%) Chuẩn đối sánh
Duffing SDOF $\omega_0=1, h=0.25, \gamma=1.0$ $8.45%$ $1.12%$ Nghiệm FPK chính xác
Cản phi tuyến bậc ba $h=0.05, \omega_0=1, \beta=4.0$ $14.20%$ $2.35%$ Nghiệm FPK chính xác
Van der Pol tự kích $\alpha \cdot \varepsilon = 0.5, S_0=1.0$ $18.60%$ $3.80%$ Mô phỏng Monte Carlo
Duffing 2-DOF ghép $\gamma_1=\gamma_2=1.0, S_0=1.0$ $12.80%$ $2.15%$ Mô phỏng Monte Carlo
Duffing cản phi tuyến (Ồn màu) $\omega_f=1.0, \alpha_f=1.0, \gamma=1.0$ $16.50%$ $3.05%$ Mô phỏng Monte Carlo
  1. Giải quyết triệt để tính bất định của LOMSEC: Luận án chứng minh rằng biến không thứ nguyên $a$ (đại diện cho bán kính miền tích phân tập trung $r$) thay vì phải lựa chọn theo trực giác, đã được nội suy và tích hợp tự động qua hàm phân bố Gauss. Biểu đồ phụ thuộc của $a$ theo tham số phi tuyến $\gamma$ cho thấy tính trơn và tính ổn định tuyệt đối của thuật toán.
  2. Độ chính xác cao trong phân tích hệ phi tuyến tự kích Van der Pol: Đối với dao động Van der Pol chịu kích động ngẫu nhiên, hiện tượng chu trình giới hạn (limit cycle) làm cho hàm mật độ xác suất lệch xa phân bố Gauss chuẩn. Phương pháp CEL kinh điển thất bại trong việc nắm bắt biên độ dao động tự kích, trong khi GLOMSEC phản ánh chính xác xu hướng dập tắt hoặc khuếch đại dao động với sai số dưới $4%$.
  3. Khả năng mở rộng hoàn hảo cho hệ nhiều bậc tự do (MDOF) và kích động ồn màu: Kết quả phân tích mô men tương quan chéo $E[x_1 x_2]$ của hệ 2-DOF và hệ mở rộng không gian trạng thái với bộ lọc ồn màu bậc hai chứng minh GLOMSEC bảo toàn cấu trúc tương tác động lực học giữa các bậc tự do, không làm méo mó ma trận hiệp phương sai.
  4. Khả năng ứng dụng thực tế trong động lực học tàu thủy: Phân tích dao động lắc ngang của tàu thủy (Ship Rolling Dynamics) dưới tác động của mô-men sóng ngẫu nhiên phi tuyến cho thấy GLOMSEC dự báo chính xác xác suất xuất hiện các góc nghiêng lớn – thông số sinh tử để đánh giá an toàn lật tàu.

Implications đa chiều

  • Về mặt lý thuyết (Theoretical Advances): Đặt nền móng vững chắc cho trường phái cơ học ngẫu nhiên đối ngẫu; chứng minh tính tương đương tiệm cận giữa các tiêu chuẩn sai số vi phân cục bộ và toán tử tích phân toàn cục.
  • Về phương pháp luận (Methodological Innovations): Cung cấp một công cụ giải tích mạnh mẽ, thay thế hiệu quả các phép mô phỏng Monte Carlo đắt đỏ về thời gian tính toán. Thuật toán lặp GLOMSEC hội tụ chỉ sau $4 \div 8$ bước lặp ma trận, tiết kiệm hơn $99.9%$ chi phí tính toán so với MCS.
  • Ứng dụng thực tiễn kỹ thuật (Practical Applications): Áp dụng trực tiếp vào thiết kế các bộ giảm chấn động lực học phi tuyến (Nonlinear Tuned Mass Dampers - TMD), thiết bị tiêu tán năng lượng, đánh giá độ mỏi kết cấu công trình biển và nhà siêu cao tầng dưới tải trọng gió bão ngẫu nhiên.
  • Khuyến nghị quy chuẩn kỹ thuật (Engineering Standards): Cung cấp cơ sở khoa học để hiệu chỉnh các tiêu chuẩn thiết kế kháng chấn và tải trọng gió, chuyển dịch từ phương pháp hệ số an toàn tĩnh tiền định sang thiết kế dựa trên độ tin cậy xác suất phi tuyến (Performance-Based Reliability Design).

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu (Boundary Conditions):

  1. Giả thiết quá trình ngẫu nhiên dừng (stationary) và ergodic: Chưa xét đến các tải trọng phi dừng có cường độ biến thiên đột ngột theo thời gian như xung sóng thần hoặc động đất cận tâm (near-fault earthquakes).
  2. Xấp xỉ phân bố xác suất Gaussian tương đương: Khi hệ thống có tính phi tuyến cực đại dẫn đến hiện tượng đa ổn định (bistability) hoặc chuyển pha hỗn loạn (stochastic chaos), giả thiết phân bố Gaussian của đáp ứng tuyến tính hóa tương đương sẽ làm giảm độ chính xác ở phần đuôi của hàm mật độ xác suất (PDF tails).
  3. Hạn chế về số bậc tự do: Mới kiểm chứng thực nghiệm số đến hệ 2-DOF; việc mở rộng cho các hệ kết cấu liên tục rời rạc hóa hàng nghìn bậc tự do (FEA) đòi hỏi tối ưu hóa cấu trúc ma trận thưa.

Các hướng nghiên cứu tiếp nối trong tương lai:

  • Phát triển GLOMSEC phi dừng (Non-stationary GLOMSEC) cho bài toán động đất và tải trọng va chạm ngẫu nhiên.
  • Kết hợp GLOMSEC với kỹ thuật đóng phi Gaussian (Non-Gaussian Closure) hoặc mô hình phân bố hỗn hợp Gauss (Gaussian Mixture Models) để dự báo chính xác vùng đuôi xác suất phục vụ đánh giá độ tin cậy phá hủy hiếm (rare failure events).
  • Tích hợp tiêu chuẩn GLOMSEC vào các phần mềm phần tử hữu hạn thương mại (ANSYS, ABAQUS, SAP2000) để phân tích động lực học ngẫu nhiên công trình quy mô lớn.
  • Mở rộng thuật toán cho hệ phi tuyến có trễ vật liệu phức tạp (mô hình Bouc-Wen, trễ đàn dẻo).

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

  • Giá trị học thuật quốc tế: Đóng góp trực tiếp vào kho tàng cơ học ngẫu nhiên thế giới; các bài báo trích xuất từ luận án được công bố trên các tạp chí chuyên ngành quốc tế uy tín, tạp chí ISI, Vietnam Journal of Mechanics và các hội nghị quốc tế chuyên ngành Cơ học tính toán.
  • Chuyển đổi công nghệ công nghiệp: Mang lại công cụ phân tích tin cậy cho các tập đoàn thiết kế hàng hải, công trình dầu khí ngoài khơi (Vietsovpetro, PVN), các viện thiết kế xây dựng công trình giao thông và nhà cao tầng.
  • Lợi ích kinh tế - xã hội: Nâng cao độ an toàn công trình, tối ưu hóa khối lượng vật liệu kết cấu, giảm thiểu rủi ro sự cố sập đổ thảm khốc do cộng hưởng ngẫu nhiên, kéo dài tuổi thọ khai thác của cơ sở hạ tầng trọng điểm quốc gia.

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

  • Nghiên cứu sinh và Giảng viên chuyên ngành Cơ học, Kỹ thuật Xây dựng, Hàng hải: Sở hữu khung lý thuyết và phương pháp mẫu mực về giải tích ngẫu nhiên phi tuyến để mở rộng cho các đề tài nghiên cứu chuyên sâu.
  • Kỹ sư R&D và Chuyên gia thiết kế kết cấu: Được trang bị công cụ tính toán nhanh, chính xác cao để phân tích dao động mỏi và ổn định động lực học công trình ngẫu nhiên.
  • Các nhà hoạch định chính sách và Cơ quan quản lý chất lượng xây dựng: Có cơ sở khoa học để biên soạn, cập nhật hệ thống quy chuẩn kỹ thuật quốc gia về tải trọng và tác động động lực học.

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 thiết lập nguyên lý lấy kỳ vọng đối ngẫu tổng thể trên trường xấp xỉ sai số địa phương. Công trình đã mở rộng trực tiếp Lý thuyết Tuyến tính hóa Tương đương của Caughey (1959) và giải quyết triệt để bài toán khép kín miền biên tích phân hữu hạn $r_x$ trong Lý thuyết LOMSEC của Lưu Xuân Hùng (2002) và Anh & Di Paola (1993).

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế tiền nhiệm?
Trả lời: So với tiêu chuẩn sai số thế năng của Elishakoff & Zhang (1992) vốn chỉ áp dụng cho phi tuyến đàn hồi, GLOMSEC xử lý đồng thời cả phi tuyến đàn hồi, phi tuyến cản và tự kích. So với phương pháp mô phỏng Monte Carlo (Crandall 2000), GLOMSEC rút ngắn thời gian tính toán hàng nghìn lần mà vẫn duy trì sai số dưới $2% - 4%$ so với nghiệm FPK chính xác.

3. Phát hiện bất ngờ nhất về mặt dữ liệu tính toán là gì?
Trả lời: Phát hiện bất ngờ nhất là trong hệ dao động Van der Pol chịu kích động ngẫu nhiên dừng, dù phân bố đáp ứng thực tế có dạng "miệng núi lửa" (hoàn toàn phi Gauss), hệ số tuyến tính hóa tổng thể xác lập qua GLOMSEC vẫn ước lượng mô men bậc hai $E[x^2]$ với độ chính xác vượt bậc (sai số chỉ $3.8%$ so với mức sai số $>18%$ của phương pháp Caughey).

4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) không?
Trả lời: Luận án cung cấp tường minh toàn bộ hệ thống phương trình giải tích, sơ đồ thuật toán lặp Atalik-Utku mở rộng, cùng quy trình tích hợp số Runge-Kutta 4 cho Monte Carlo, cho phép các nhà nghiên cứu độc lập lập trình tái lập chính xác $100%$ kết quả trên MATLAB.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo ra sao?
Trả lời: Định hướng mở rộng GLOMSEC trên 3 trục chính: (1) Động lực học ngẫu nhiên phi dừng kết hợp biến đổi Wavelet; (2) Tích hợp học máy (Physics-Informed Neural Networks - PINNs) để tăng tốc giải thuật ma trận nhiều chiều; (3) Tối ưu hóa điều khiển thích nghi thời gian thực cho kết cấu thông minh chịu siêu bão và động đất.

Kết luận

  1. Xây dựng thành công Tiêu chuẩn sai số bình phương trung bình địa phương – tổng thể (GLOMSEC), hoàn thiện bước phát triển nhảy vọt cho phương pháp tuyến tính hóa tương đương ngẫu nhiên.
  2. Khép kín hoàn toàn về mặt toán học bài toán xác định miền tập trung đáp ứng của tiêu chuẩn LOMSEC, loại bỏ yếu tố chủ quan trong tính toán giải tích.
  3. Thiết lập hệ công thức và thuật toán ma trận tổng quát cho hệ phi tuyến nhiều bậc tự do (MDOF) và hệ chịu kích động ngẫu nhiên ồn màu dải hẹp.
  4. Chứng minh định lượng tính chính xác vượt trội của GLOMSEC trên các hệ dao động kinh điển (Duffing, Van der Pol, cản phi tuyến) và bài toán thực tế lắc ngang tàu thủy, kiểm chứng vững chắc qua nghiệm FPK và Monte Carlo.
  5. Mở ra nhánh nghiên cứu liên ngành đột phá kết hợp giải tích đối ngẫu với tính toán kết cấu công trình biển và xây dựng hiện đại, khẳng định vị thế tiên phong của cơ học ngẫu nhiên Việt Nam trên trường quốc tế.