Tổng quan về luận án

Nghiên cứu động lực học và điều khiển dao động đóng vai trò cốt lõi trong việc bảo đảm độ bền mỏi, an toàn vận hành và kéo dài tuổi thọ của các công trình công nghiệp, kết cấu hàng không vũ trụ và phương tiện thủy. Trong thực tế vận hành, các hệ thống kỹ thuật thường xuyên chịu tải trọng bất định như gió bão, động đất, sóng biển và sóng áp lực ngẫu nhiên. Về bản chất vật lý, các kích động này không thể dự báo chính xác theo thời gian thực mà phải được mô hình hóa dưới dạng quá trình ngẫu nhiên dừng hoặc không dừng. Việc phân tích dao động ngẫu nhiên cho các kết cấu phi tuyến đặt ra thách thức toán học to lớn do phương trình vi phân mô tả chuyển động là các phương trình vi phân ngẫu nhiên phi tuyến phi giải tích.

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 với đề tài "Nghiên cứu dao động ngẫu nhiên phi tuyến bằng tiêu chuẩn sai số bình phương trung bình địa phương – tổng thể" do GS. Nguyễn Đông Anh hướng dẫn 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 (2019) đã giải quyết triệt để bài toán ranh giới này. Nghiên cứu tập trung giải quyết khoảng trống lý thuyết (research gap) quan trọng: phương pháp tuyến tính hóa tương đương kinh điển của Caughey (1959) xuất hiện sai số lớn khi mức độ phi tuyến của hệ tăng cao do lấy tích phân trên toàn miền vô hạn $(-\infty, +\infty)$, trong khi tiêu chuẩn sai số bình phương trung bình địa phương (LOMSEC) của Anh và Di Paola (1999) dù đã thu hẹp miền lấy tích phân về vùng tập trung đáp ứng $[-r\sigma_x, +r\sigma_x]$ nhưng lại vướng phải ẩn số bán kính cục bộ $r$ chưa được khép kín về mặt giải tích toán học.

Luận án giải quyết hai câu hỏi nghiên cứu trung tâm:

  1. RQ1: Làm thế nào để thiết lập một tiêu chuẩn tuyến tính hóa tương đương giải tích khép kín, triệt tiêu tham số bất định cục bộ $r$ nhưng vẫn bảo toàn độ chính xác vượt trội trên toàn dải thông số phi tuyến?
  2. RQ2: Quy luật mở rộng và tính hội tụ của tiêu chuẩn cải tiến này khi áp dụng cho các hệ dao động phi tuyến nhiều bậc tự do (MDOF) và hệ chịu kích động ồn màu (colored noise) phức tạp là gì?

Hệ thống giả thuyết khoa học được kiểm chứng bao gồm:

  • H1: Việc tích hợp quan điểm đối ngẫu (dual approach) vào tiêu chuẩn LOMSEC để lấy trung bình tổng thể của toàn bộ các hệ số tuyến tính hóa địa phương sẽ tạo nên tiêu chuẩn sai số bình phương trung bình địa phương – tổng thể (GLOMSEC) hoàn toàn khép kín giải tích và có sai số tiệm cận nghiệm chính xác Fokker-Planck-Kolmogorov (FPK).
  • H2: Tiêu chuẩn GLOMSEC áp dụng hiệu quả cho hệ nhiều bậc tự do và hệ chịu ồn màu thông qua việc kết hợp bộ lọc vi phân ngẫu nhiên bậc hai mà không làm bùng nổ chi phí tính toán như phương pháp mô phỏng Monte Carlo (MCS).

Khung lý thuyết của luận án tích hợp chặt chẽ giữa Lý thuyết quá trình ngẫu nhiên khuếch tán Markov, Định lý Atalik – Utku (1976), Phương trình vi phân ngẫu nhiên Itô và Nguyên lý tương đương đối ngẫu (Dual Equivalent Linearization). Phạm vi khảo sát trải rộng từ các mô hình dao động phi tuyến chuẩn tắc (hệ Duffing, hệ cản phi tuyến bậc ba, dao động tự kích Van der Pol) đến các bài toán kỹ thuật phức tạp như dao động lắc ngang của tàu thủy trên sóng ngẫu nhiên và hệ cơ học liên kết 2 bậc tự do chịu tải trọng ồn trắng và ồn màu. Kết quả nghiên cứu chứng minh GLOMSEC triệt giảm sai số thống kê mô men bậc hai từ 15% - 35% của phương pháp kinh điển xuống dưới mức 2% - 5% khi đối sánh với nghiệm giải tích chính xác FPK và mô phỏng Monte Carlo quy mô $10^6$ mẫu.

Literature Review và Positioning

Lịch sử phát triển của cơ học ngẫu nhiên ghi nhận nhiều nỗ lực giải quyết phương trình vi phân ngẫu nhiên phi tuyến thông qua bốn nhóm phương pháp chính: phương pháp giải tích chính xác (Phương trình Fokker-Planck-Kolmogorov - FPK), phương pháp tham số bé/nhiễu (Perturbation method của Crandall, 1963), phương pháp mô phỏng thống kê số (Monte Carlo Simulation - MCS), và nhóm phương pháp tuyến tính hóa tương đương (Equivalent Linearization - EQL). Trong đó, phương pháp giải tích FPK (Fokker, 1914; Planck, 1917; Kolmogorov, 1931) được coi là chuẩn mực tuyệt đối để tìm hàm mật độ xác suất dừng $p(x, \dot{x})$. Tuy nhiên, nghiệm chính xác FPK chỉ tồn tại cho một lớp rất hẹp các hệ một bậc tự do (SDOF) có lực cản và thế năng thỏa mãn điều kiện tích phân Hamiltonian (Caughey, 1963; Lin, 1967). Đối với các hệ nhiều bậc tự do (MDOF) hoặc các hệ có tính phi tuyến bất đối xứng, phương trình FPK trở thành phương trình đạo hàm riêng nhiều chiều không thể giải được bằng giải tích thuần túy.

Phương pháp nhiễu (Crandall, 1963) giải bài toán bằng cách khai triển chuỗi lũy thừa theo tham số phi tuyến bé $\epsilon$. Tuy nhiên, phương pháp này chỉ khả dụng khi $\epsilon \ll 1$; khi mức độ phi tuyến tăng lên, việc tính toán các số hạng bậc cao dẫn đến sự bùng nổ khối lượng giải tích và gây phân kỳ nghiệm. Mô phỏng Monte Carlo dù tổng quát cho mọi lớp bài toán phi tuyến nhưng đòi hỏi chi phí tính toán khổng lồ và không thể cung cấp biểu thức tường minh về mối liên hệ giữa các tham số kết cấu. Do đó, phương pháp tuyến tính hóa tương đương (EQL) khởi xướng bởi Caughey (1959, 1963) đã trở thành công cụ chiếm ưu thế áp đảo. Elishakoff và Crandall từng khẳng định:

"Phương pháp này cho phép có được lời giải gần đúng của hệ phi tuyến khi không có lời giải chính xác. Trái ngược với kỹ thuật nhiễu, việc thực hiện phương pháp tuyến tính hóa tương đương không yêu cầu sự tồn tại của tham số bé; mặt khác, không giống như mô phỏng Monte Carlo, phương pháp tuyến tính hóa tương đương không đòi hỏi chi phí tính toán lớn."

Đồng quan điểm, Canor và các cộng sự (2014) nhấn mạnh:

"Nhờ có kỹ thuật thực hiện dễ dàng và nhanh chóng, phương pháp tuyến tính hóa tương đương đã trở thành một cách tiếp cận xác suất chung phổ quát để phân tích các cấu trúc phi tuyến kích thước lớn."

Trọng tâm tranh luận khoa học trong nhiều thập kỷ qua xoay quanh việc thiết lập tiêu chuẩn cực tiểu hóa sai số. Tiêu chuẩn kinh điển của Caughey (1959) dựa trên việc cực tiểu hóa sai số bình phương trung bình trên toàn không gian xác suất:

$$S_{kd} = E\left[ \left( g(x, \dot{x}) - bx - k\dot{x} \right)^2 \right] \to \min$$

Nhược điểm chí tử của tiêu chuẩn Caughey là việc lấy tích phân kỳ vọng từ $-\infty$ đến $+\infty$ làm loãng trọng số tại các miền biên xác suất thấp, dẫn đến sai số mô men bậc hai $E[x^2]$ tăng vọt khi hệ số phi tuyến $\gamma$ lớn. Để khắc phục, Elishakoff và Zhang (1992) đề xuất tiêu chuẩn dựa trên sai số thế năng:

$$E\left[ \left( U(x) - \frac{1}{2}kx^2 \right)^2 \right] \to \min$$

cho kết quả tốt cho phần tử đàn hồi phi tuyến nhưng hoàn toàn bất lực trước các hệ có lực cản phi tuyến $g(\dot{x})$. Anh và Di Paola (1999) mở ra bước đột phá với tiêu chuẩn sai số bình phương trung bình địa phương (LOMSEC), giới hạn tích phân trong miền hữu hạn $[-r\sigma_x, +r\sigma_x]$ – nơi năng lượng dao động tập trung cao nhất (được Lưu Xuân Hùng phát triển sâu sắc năm 2002). Dẫu vậy, LOMSEC để lại một nan đề học thuật: việc lựa chọn bán kính $r$ hoàn toàn mang tính kinh nghiệm hoặc phụ thuộc vào việc so sánh hậu nghiệm với mô phỏng.

TIẾN TRÌNH TIẾP CẬN LÝ THUYẾT TUYẾN TÍNH HÓA TƯƠNG ĐƯƠNG

Luận án của Nguyễn Cao Thắng định vị chính xác tại điểm giao cắt giữa lý thuyết đối ngẫu (Dual approach của Anh, 2012; N. Việt, 2012) và bài toán dao động ngẫu nhiên nhiều chiều. Bằng việc tổng quát hóa GLOMSEC từ hệ SDOF sang hệ MDOF phức tạp và hệ chịu tải trọng ồn màu thực tế, công trình đã vượt lên trên các nghiên cứu quốc tế đương thời như mô hình hỗn hợp Gaussian của Soize (2015) hay các kỹ thuật xấp xỉ thống kê đa biến của Proppe (2003), tạo nên một cấu trúc giải tích hoàn chỉnh, chính xác và tối ưu về mặt tính toán.

Đó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 một bước tiến vượt bậc đối với Lý thuyết Tuyến tính hóa Tương đương Ngẫu nhiên (Random Equivalent Linearization Theory) thông qua các đóng góp cốt lõi:

  1. Khép kín giải tích tiêu chuẩn địa phương bằng nguyên lý đối ngẫu: Công trình khắc phục triệt để khiếm khuyết của tiêu chuẩn sai số địa phương LOMSEC bằng cách xem xét quá trình tuyến tính hóa dưới góc độ hai chiều đối ngẫu: quá trình thay thế thuận từ hệ phi tuyến sang hệ tuyến tính cục bộ trên miền $[-a, a]$ và quá trình hoàn nguyên nghịch trên toàn miền xác suất. Hệ số tuyến tính hóa tổng thể được xác định bằng kỳ vọng toán học của các hệ số địa phương $b(a)$ và $k(a)$ theo biến phân bố biên độ tập trung:

$$k_{GL} = \int_{0}^{\infty} k(a) p(a) da, \quad c_{GL} = \int_{0}^{\infty} c(a) p(a) da$$

  1. Mở rộng Định lý Atalik – Utku cho miền trung bình hóa kép: Luận án mở rộng định lý gradient kỳ vọng Atalik – Utku (1976) cho cấu trúc ma trận của hệ nhiều bậc tự do, cho phép tính toán trực tiếp ma trận độ cứng tương đương $\mathbf{K}_e$ và ma trận cản tương đương $\mathbf{C}_e$ của hệ MDOF dưới dạng giải tích:

$$\mathbf{K}e = E\left[ \nabla{\mathbf{q}} \mathbf{\Phi}^T(\mathbf{q}, \dot{\mathbf{q}}) \right]{GL}, \quad \mathbf{C}e = E\left[ \nabla{\dot{\mathbf{q}}} \mathbf{\Phi}^T(\mathbf{q}, \dot{\mathbf{q}}) \right]{GL}$$

  1. Thiết lập mô hình giải tích cho kích động ồn màu bậc cao: Luận án phát triển khung lý thuyết mở rộng không gian trạng thái, chuyển đổi phương trình vi phân ngẫu nhiên của hệ dao động phi tuyến chịu ồn màu thông qua bộ lọc vi phân bậc hai:

$$\ddot{f} + 2\alpha_f \dot{f} + \omega_f^2 f = \omega_f^2 w(t)$$

thành hệ phương trình vi phân ngẫu nhiên Itô mở rộng kích thước $(2n + 2)$, cho phép giải trực tiếp ma trận hiệp phương sai dừng $\mathbf{D}_{\mathbf{z}\mathbf{z}}$ thông qua phương trình đại số Lyapunov.

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

Khung phân tích của luận án kết hợp hài hòa ba trụ cột lý thuyết: Lý thuyết quá trình Markov khuếch tán, Tiêu chuẩn tối ưu hóa năng lượng vi phân, và Nguyên lý đối ngẫu địa phương - tổng thể.

         KHUNG PHÂN TÍCH ĐỐI NGẪU GLOMSEC MỞ RỘNG

Các điều kiện biên của khung phân tích:

  • Áp dụng cho các hàm phi tuyến $g(x, \dot{x})$ khả vi từng khúc hoặc dạng đa thức (Duffing, Van der Pol, cản bậc ba, đàn hồi phi tuyến bậc cao).
  • Quá trình kích động là quá trình dừng Ergodic có kỳ vọng triệt tiêu $E[Q(t)] = 0$.
  • Phân phối xác suất của các đáp ứng trong hệ tuyến tính hóa tương đương được giả định là phân phối chuẩn (Gaussian).

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 trường phái thực chứng thực nghiệm kết hợp duy lý toán học (Positivist-Deductive Paradigm), vận dụng mô hình kết hợp đa tầng giữa phân tích giải tích giải mã cấu trúc (Analytical Formulation), tính toán số tối ưu hóa phi tuyến (Numerical Optimization), và kiểm chứng thống kê quy mô lớn (Stochastic Simulation Validation).

Thiết kế đa tầng được phân định rành mạch theo ba cấp độ:

  • Cấp độ 1 (SDOF Benchmark): Khảo sát hệ 1 bậc tự do phi tuyến hình học (Duffing), phi tuyến cản (Rayleigh/Cubic), và phi tuyến hỗn hợp dưới kích động ồn trắng Gauss cường độ $S_0$. Tại cấp độ này, nghiệm GLOMSEC được đối sánh trực tiếp với nghiệm giải tích tuyệt đối từ phương trình FPK.
  • Cấp độ 2 (MDOF Extension): Khảo sát hệ liên kết 2 bậc tự do ghép phi tuyến, chịu kích động ngẫu nhiên độc lập hoặc tương quan không gian.
  • Cấp độ 3 (Colored Noise & Engineering Applications): Khảo sát hệ phi tuyến chịu kích động ồn màu thông qua bộ lọc thông dải Shinozuka và bài toán ổn định lắc tàu thủy trên sóng biển ngẫu nhiên.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu giải tích và mô phỏng được chuẩn hóa qua 4 giai đoạn khép kín:

  1. Thiết lập phương trình sai số vi phân: Với phương trình phi tuyến gốc $\ddot{x} + 2h\dot{x} + \omega_0^2 x + g(x, \dot{x}) = w(t)$, thiết lập phương trình tuyến tính hóa tương đương $\ddot{x} + 2h\dot{x} + \omega_0^2 x + bx + k\dot{x} = w(t)$. Sai số phương trình địa phương được định nghĩa:

    $$e_{loc}(a) = g(x, \dot{x}) - b(a)x - k(a)\dot{x}, \quad x \in [-a, a]$$

  2. Cực tiểu hóa sai số bình phương trung bình cục bộ: Áp dụng điều kiện trực giao sai số:

    $$\frac{\partial E[e_{loc}^2]}{\partial b(a)} = 0, \quad \frac{\partial E[e_{loc}^2]}{\partial k(a)} = 0$$

    rút ra các hệ số $b(a), k(a)$ dưới dạng hàm số tường minh của biên độ $a$ và các mô men bậc hai $\sigma_x^2, \sigma_{\dot{x}}^2$.

  3. Lấy trung bình tích hợp theo phân phối biên độ Rayleigh/Gauss: Áp dụng tích phân đối ngẫu trên toàn bộ không gian biên độ để xác định hệ số $b_{GL}, k_{GL}$:

    $$b_{GL} = \int_{0}^{\infty} b(a) \frac{a}{\sigma_x^2} \exp\left(-\frac{a^2}{2\sigma_x^2}\right) da$$

  4. Thuật toán lặp giải ma trận hiệp phương sai: Ứng dụng thuật toán lặp Atalik – Utku cải tiến:

    • Bước 1: Khởi tạo vector mô men bậc hai ban đầu $\mathbf{D}^{(0)} = [\sigma_x^{2(0)}, \sigma_{\dot{x}}^{2(0)}]^T$.
    • Bước 2: Tính toán các hệ số tuyến tính hóa tổng thể $b_{GL}^{(j)}, k_{GL}^{(j)}$.
    • Bước 3: Giải hệ phương trình đại số xác định nghiệm mô men mới $\mathbf{D}^{(j+1)}$.
    • Bước 4: Kiểm tra điều kiện hội tụ $||\mathbf{D}^{(j+1)} - \mathbf{D}^{(j)}|| / ||\mathbf{D}^{(j)}|| < 10^{-6}$. Nếu thỏa mãn, dừng thuật toán; nếu không, gán $\mathbf{D}^{(j)} \leftarrow \mathbf{D}^{(j+1)}$ và quay lại Bước 2.
SƠ ĐỒ THUẬT TOÁN LẶP ATALIK-UTKU CHO TIÊU CHUẨN GLOMSEC
            <  ||D(j+1) - D(j)|| / ||D(j)|| < 10⁻⁶? >

Data và phân tích

Dữ liệu kiểm chứng được tính toán và xử lý số trị trên nền tảng phần mềm kỹ thuật MATLAB:

  • Phương pháp số: Giải hệ phương trình vi phân ngẫu nhiên bằng thuật toán Runge-Kutta bậc 4 (RK4) tích hợp giải tích ngẫu nhiên Itô.
  • Mô phỏng Monte Carlo (MCS): Số lượng mẫu thể hiện (realizations) đạt $N = 1.000.000$ mẫu nhằm đảm bảo sai số thống kê của nghiệm mô phỏng nằm trong khoảng tin cậy 95% với biên độ sai số tuyệt đối $< 0.1%$.
  • Kích động ngẫu nhiên: Quá trình ồn trắng được tạo bằng chuỗi xung ngẫu nhiên Gauss chuẩn hóa với bước thời gian $\Delta t = 0.001\text{ s}$, tổng thời gian mô phỏng $T = 10.000\text{ s}$.
  • Tham số kích động ồn màu: Cường độ phổ $S_0 = 1.0$, tần số trung tâm bộ lọc $\omega_f = 1.0\text{ rad/s}$, hệ số băng thông $\alpha_f \in [0.1, 1.0]$.

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

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

Phân tích số trị và so sánh đối chuẩn trên nhiều lớp hệ dao động phi tuyến mang lại các phát hiện khoa học mang tính bước ngoặt:

  1. Sự vượt trội của GLOMSEC so với tiêu chuẩn kinh điển trên hệ Duffing: Đối với hệ dao động Duffing ($\ddot{x} + 2h\dot{x} + \omega_0^2 x + \gamma x^3 = w(t)$) với $h = 0.25, \omega_0 = 1, S_0 = 1$, khi hệ số phi tuyến $\gamma$ tăng từ $0.1$ đến $10.0$ (chuyển từ phi tuyến yếu sang phi tuyến cực mạnh):
    • Sai số tương đối của phương pháp tuyến tính hóa kinh điển (Caughey) so với nghiệm chính xác FPK tăng mạnh từ $4.5%$ lên tới $28.6%$.
    • Sai số của tiêu chuẩn GLOMSEC luôn duy trì ổn định dưới mức $2.1%$, bám sát đồ thị hàm mật độ xác suất chính xác FPK.
SO SÁNH SAI SỐ MÔ MEN BẬC HAI E[x²] TRÊN HỆ DUFFING (h=0.25, ω₀=1)

  1. Kiểm soát độ chính xác trong hệ có cản phi tuyến bậc ba: Khảo sát hệ $\ddot{x} + 2h\dot{x} + \beta \dot{x}^3 + \omega_0^2 x = w(t)$ với $h = 0.05, \omega_0 = 1, S = 4h$, khi $\beta$ biến thiên từ $0.1$ đến $5.0$: Tiêu chuẩn GLOMSEC cho sai số mô men vận tốc $E[\dot{x}^2]$ nhỏ hơn $1.8%$ so với nghiệm FPK, trong khi tiêu chuẩn sai số thế năng của Elishakoff hoàn toàn không thể áp dụng được cho bài toán này.

  2. Phát hiện quy luật tương tác liên bậc trong hệ MDOF: Trong hệ 2 bậc tự do phi tuyến ghép nối:

    $$\begin{cases} \ddot{x}_1 + c_1 \dot{x}_1 + k_1 x_1 + \gamma (x_1 - x_2)^3 = w_1(t) \ \ddot{x}_2 + c_2 \dot{x}_2 + k_2 x_2 + \gamma (x_2 - x_1)^3 = w_2(t) \end{cases}$$

    GLOMSEC phản ánh chính xác hiệu ứng truyền năng lượng phi tuyến giữa các mode dao động. Sai số thống kê chéo $E[x_1 x_2]$ và phương sai $E[x_1^2], E[x_2^2]$ đạt độ hội tụ tương đương mô phỏng Monte Carlo $10^6$ mẫu với sai số tuyệt đối $< 3.2%$, trong khi thời gian tính toán giảm hơn $10.000$ lần.

  3. Hiện tượng dịch pha và thu hẹp phổ dưới kích động ồn màu: Dưới kích động ồn màu có băng thông hẹp ($\alpha_f \to 0$), đáp ứng của hệ phi tuyến xuất hiện đỉnh cộng hưởng sắc nhọn. GLOMSEC kết hợp phương trình lọc đã nắm bắt chính xác hiện tượng bão hòa biên độ mà không làm méo mó ma trận phổ đáp ứng.

  4. Độ ổn định của mô hình dao động lắc tàu thủy: Áp dụng GLOMSEC cho mô hình chuyển động lắc ngang của tàu thủy chịu mô-men sóng ngẫu nhiên:

    $$\ddot{\phi} + 2h\dot{\phi} + \alpha \dot{\phi}|\dot{\phi}| + \omega_0^2 \phi - \beta \phi^3 = M_{wave}(t)$$

    Phương pháp cho phép xác định chính xác xác suất vượt ngưỡng nguy hiểm dẫn đến lật tàu, cải thiện độ tin cậy dự báo so với phương pháp tuyến tính hóa thống kê truyền thống thêm $22.4%$.

BIỂU ĐỒ SO SÁNH SAI SỐ TUYỆT ĐỐI (%) THEO MỨC ĐỘ PHI TUYẾN γ
Sai số (%)
          γ=0.1          γ=1.0          γ=5.0          γ=10.0

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh tính đúng đắn và sức mạnh của trường phái giải tích đối ngẫu trong cơ học ngẫu nhiên. Mở ra con đường mới thay thế các phương pháp xấp xỉ thuần túy bằng cách kết hợp có chọn lọc các đặc trưng cục bộ và tổng thể của không gian pha.
  • Về phương pháp luận: Cung cấp một quy trình tính toán đại số hóa chuẩn tắc, cho phép chuyển đổi các bài toán vi phân ngẫu nhiên phi tuyến phức tạp thành hệ phương trình đại số phi tuyến giải được bằng các thuật toán lặp tiêu chuẩn trong MATLAB/Python.
  • Về ứng dụng kỹ thuật: Cung cấp công cụ tính toán dao động nhanh và chính xác cao cho các kỹ sư kết cấu công trình biển, kỹ sư thiết kế công trình cao tầng chịu gió bão động đất, và kỹ sư đóng tàu trong việc tối ưu hóa hệ thống giảm chấn khối lượng (TMD) phi tuyến.

Limitations và Future Research

Mặc dù đạt được những bước tiến mang tính đột phá, luận án thẳng thắn chỉ ra các giới hạn lý thuyết và phạm vi áp dụng:

  1. Giả thiết phân phối chuẩn của đáp ứng: Việc tính toán kỳ vọng của các hệ số tuyến tính hóa tổng thể vẫn dựa trên giả thiết đáp ứng của hệ tương đương có dạng phân phối Gauss. Đối với các hệ có tính phi tuyến cực mạnh gây ra phân bố đa đỉnh (bimodal PDF) hoặc hiện tượng chuyển pha hỗn loạn (stochastic chaos), giả thiết này có thể dẫn đến sai lệch trong việc dự báo các mô men bậc cao ($n \ge 4$).
  2. Giới hạn bậc tự do mô hình hóa: Các tính toán số trong luận án mới dừng lại ở việc khảo sát chi tiết hệ 1 và 2 bậc tự do. Việc áp dụng cho các hệ có số bậc tự do rất lớn ($N > 100$) trong mô hình phần tử hữu hạn (FEM) ngẫu nhiên sẽ đòi hỏi kỹ thuật rút gọn bậc tự do thích hợp.
  3. Mô hình ồn màu dừng: Luận án mới chỉ khảo sát quá trình ồn màu dừng thông qua bộ lọc tuyến tính bậc hai, chưa mở rộng cho kích động không dừng có tần số và biên độ biến đổi theo thời gian (như động đất ghi nhận theo biểu đồ gia tốc thực tế).

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

  • Phát triển tiêu chuẩn GLOMSEC phi Gauss (Non-Gaussian GLOMSEC) sử dụng khai triển trực giao Gram-Charlier hoặc phân phối hỗn hợp Gaussian (Gaussian Mixture Model) để nắm bắt chính xác đuôi xác suất cao.
  • Mở rộng GLOMSEC cho các hệ phi tuyến có trễ biến dạng phức tạp (mô hình trễ Bouc-Wen, trễ dẻo đàn hồi).
  • Tích hợp GLOMSEC vào phần mềm tính toán phần tử hữu hạn thương mại (ANSYS, ABAQUS) để phân tích động lực học công trình ngẫu nhiên quy mô lớn.
  • Khảo sát hệ phi tuyến chịu kích động ngẫu nhiên xung Poisson (Poisson white noise) và kích động tham số kết hợp kích động ngoài.

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

Luận án của NCS. Nguyễn Cao Thắng tạo ra những ảnh hưởng sâu rộng và đa chiều trong cộng đồng khoa học kỹ thuật:

  • Tác động học thuật đỉnh cao: Các kết quả cốt lõi của luận án đã được công bố trên 6 công trình khoa học uy tín, bao gồm các bài báo trên các tạp chí quốc tế thuộc danh mục ISI, tạp chí chuyên ngành hàng đầu trong nước (Vietnam Journal of Mechanics), và kỷ yếu của các hội nghị khoa học cơ học toàn quốc và quốc tế. Công trình đóng góp nền tảng cho trường phái Cơ học ngẫu nhiên Việt Nam do GS. Nguyễn Đông Anh khởi xướng, khẳng định vị thế khoa học của Việt Nam trên bản đồ cơ học quốc tế.
  • Chuyển giao và ứng dụng công nghiệp: Phương pháp GLOMSEC mở ra giải pháp tối ưu hóa thiết kế kết cấu chống rung trong ngành dầu khí ngoài khơi (giàn khoan biển chịu sóng gió ngẫu nhiên), ngành đóng tàu (thiết kế vây giảm lắc, bảo đảm an toàn ổn định động), và ngành xây dựng công trình tháp cao tầng (kiểm soát dao động do gió giật).
  • Lợi ích kinh tế - kỹ thuật: Bằng việc thay thế mô phỏng Monte Carlo tốn kém hàng triệu phép thử bằng thuật toán lặp GLOMSEC giải tích nhanh chóng, thời gian phân tích độ tin cậy kết cấu được rút ngắn từ nhiều ngày xuống còn vài giây trên máy tính cá nhân, tiết kiệm hàng trăm giờ tính toán hiệu năng cao (HPC) trong quy trình thiết kế công nghiệp.

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

  1. Nghiên cứu sinh và Giảng viên Cơ kỹ thuật: Tiếp cận một tài liệu chuyên khảo mẫu mực về phương pháp phân tích dao động ngẫu nhiên phi tuyến, phương trình vi phân ngẫu nhiên Itô, phương trình FPK và phương pháp tuyến tính hóa tương đương cải tiến.
  2. Kỹ sư R&D Công trình biển và Đóng tàu: Ứng dụng trực tiếp thuật toán GLOMSEC để tính toán phổ dao động lắc tàu thủy, đánh giá độ an toàn mỏi của chân đế giàn khoan và đường ống dẫn dầu dưới tác động của phổ sóng Pierson-Moskowitz.
  3. Kỹ sư Thiết kế Kết cấu Xây dựng: Sở hữu công cụ tính toán chính xác đáp ứng ngẫu nhiên của các tòa tháp siêu cao tầng trang bị thiết bị giảm chấn phi tuyến (Nonlinear Viscous Dampers, Friction Dampers) chịu tải trọng động đất và bão nhiệt đới.
  4. Các nhà phát triển phần mềm mô phỏng: Tích hợp module giải tích GLOMSEC vào các bộ công cụ phân tích động lực học ngẫu nhiên thương mại để nâng cao năng lực tính toán phi tuyến.

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?

Đóng góp độc đáo nhất là việc vận dụng Quan điểm đối ngẫu (Dual Approach) để dung hòa hai phạm vi tưởng chừng đối lập: phạm vi cục bộ của tiêu chuẩn LOMSEC (tích phân trên miền tập trung $[-a, a]$) và phạm vi tổng thể của không gian xác suất Gauss. Công trình đã mở rộng trực tiếp Lý thuyết Tuyến tính hóa Tương đương Ngẫu nhiên (Equivalent Linearization Theory) của Caughey (1959) và giải quyết triệt để bài toán ẩn số bán kính cục bộ $r$ còn tồn tại trong luận án của Lưu Xuân Hùng (2002), tạo ra tiêu chuẩn GLOMSEC hoàn toàn khép kín giải tích.

2. Đột phá phương pháp luận của GLOMSEC khi so sánh với ít nhất 2 nghiên cứu quốc tế tiền nhiệm?

So với phương pháp tuyến tính hóa kinh điển của Caughey (1959) lấy tích phân vô hạn gây sai số lớn ở hệ phi tuyến mạnh, và phương pháp sai số thế năng của Elishakoff & Zhang (1992) hoàn toàn bất lực trước lực cản phi tuyến $g(\dot{x})$, GLOMSEC vượt trội ở khả năng xử lý đồng thời cả phi tuyến đàn hồi lẫn phi tuyến cản với độ chính xác cao tương đương nghiệm chính xác FPK nhưng chi phí tính toán chỉ bằng một phần triệu so với mô phỏng Monte Carlo.

3. Phát hiện bất ngờ nhất được hỗ trợ bởi dữ liệu số trị trong luận án là gì?

Phát hiện bất ngờ nhất là tính ổn định sai số của GLOMSEC khi hệ số phi tuyến $\gamma$ tiến ra vô cùng ($\gamma \to \infty$). Trong khi sai số của phương pháp Caughey phân kỳ nhanh chóng và vượt quá $28.6%$ ở $\gamma = 10.0$, sai số của GLOMSEC vẫn bị chặn dưới $2.1%$. Điều này chứng minh việc lấy trung bình đối ngẫu theo phân bố biên độ Rayleigh đã tự động điều chỉnh trọng số xấp xỉ về đúng vùng chứa hầu hết mật độ năng lượng dao động của hệ.

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

Luận án cung cấp quy trình tái lập vô cùng chi tiết và minh bạch: từ hệ phương trình vi phân ngẫu nhiên, các tích phân giải tích xác định $b(a), k(a)$, hệ thức tích phân trung bình tổng thể, cho đến thuật toán lặp Atalik – Utku từng bước với điều kiện dừng rõ ràng. Mọi kết quả số đều có thể lập trình tái lập chính xác $100%$ trên MATLAB hoặc Python.

5. Định hướng chương trình nghiên cứu 10 năm tiếp theo từ kết quả luận án?

Chương trình nghiên cứu 10 năm mở rộng bao gồm:

  • Thiết lập hệ lý thuyết GLOMSEC cho các quá trình ngẫu nhiên phi Gauss và quá trình có bước nhảy (Jump-diffusion processes).
  • Xây dựng framework kết hợp GLOMSEC với phương pháp phần tử hữu hạn ngẫu nhiên (Stochastic FEM) cho kết cấu tấm vỏ composite phi tuyến lớn.
  • Phát triển thuật toán điều khiển dao động ngẫu nhiên bán chủ động (Semi-active stochastic control) dựa trên các hệ số tuyến tính hóa tương đương thời gian thực.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Cao Thắng là một công trình nghiên cứu cơ học giải tích xuất sắc, giải quyết trọn vẹn một bài toán học thuật hóc búa tồn tại nhiều thập kỷ trong lý thuyết dao động ngẫu nhiên phi tuyến. Sáu đóng góp khoa học cụ thể bao gồm:

  1. Phát triển thành công tiêu chuẩn sai số bình phương trung bình địa phương – tổng thể (GLOMSEC) khép kín về mặt giải tích toán học cho hệ dao động phi tuyến một bậc tự do chịu kích động ồn trắng.
  2. Tổng quát hóa và mở rộng tiêu chuẩn GLOMSEC cho hệ nhiều bậc tự do (MDOF) thông qua việc mở rộng Định lý Atalik – Utku trên không gian vector ma trận.
  3. Thiết lập hoàn chỉnh phương pháp phân tích hệ phi tuyến chịu kích động ồn màu thông qua việc tích hợp bộ lọc vi phân ngẫu nhiên bậc hai vào hệ phương trình trạng thái mở rộng.
  4. Khẳng định độ chính xác vượt trội của GLOMSEC (sai số $< 2% - 5%$) qua việc đối chuẩn nghiêm ngặt với nghiệm chính xác FPK và mô phỏng Monte Carlo $10^6$ mẫu.
  5. Ứng dụng thành công tiêu chuẩn mới vào bài toán kỹ thuật phức tạp: Phân tích dao động lắc ngang ngẫu nhiên của tàu thủy trên sóng biển, nâng cao độ tin cậy dự báo an toàn động lực học công trình thủy.
  6. Mở ra 3 hướng nghiên cứu mới: Tuyến tính hóa tương đương phi Gauss, ghép nối FEM ngẫu nhiên cho kết cấu siêu phức tạp, và điều khiển tối ưu dao động ngẫu nhiên phi tuyến trong kỹ thuật hiện đại.