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 đặt trọng tâm lớn vào việc phân tích đáp ứng của các kết cấu phi tuyến chịu tác động của tải trọng ngẫu nhiên như động đất, gió bão, sóng biển hoặc rung động mặt đường. Luận án tiến sĩ cơ học với tiêu đề "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 (chuyên ngành Cơ học vật rắn, mã số 62440107, 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) là một công trình tiên phong giải quyết bài toán cốt lõi: nâng cao độ chính xác trong phân tích xấp xỉ giải tích cho các hệ ngẫu nhiên nhiều bậc tự do (MDOF) có tính phi tuyến mạnh.

Trong lý thuyết dao động ngẫu nhiên, phương pháp tuyến tính hóa tương đương (Equivalent Linearization - EQL), khởi xướng bởi Booton (1954), Kazakov (1954) và Caughey (1956), là công cụ hữu hiệu và phổ biến nhất nhờ tính khả thi tính toán cao đối với hệ nhiều bậc tự do so với phương pháp nhiễu (Crandall, 1963) hay phương pháp trung bình hóa (Spanos, 1992). Tuy nhiên, khoảng trống nghiên cứu (research gap) lớn của phương pháp EQL truyền thống dựa trên tiêu chuẩn sai số bình phương trung bình cực tiểu (Mean Square Error - MSE) là hiện tượng suy giảm độ chính xác nghiêm trọng khi tham số phi tuyến $\varepsilon$ tăng từ phi tuyến yếu ($O(\varepsilon) \ll 1$) sang phi tuyến mạnh ($O(\varepsilon) \ge 1$), dẫn đến sai số ước tính phương sai đáp ứng có thể vượt quá 20–40% so với nghiệm chính xác hoặc mô phỏng số Monte-Carlo.

Luận án đặt ra các câu hỏi nghiên cứu và giả thuyết khoa học cụ thể:

  • Câu hỏi nghiên cứu 1 (Q1): Làm thế nào để thiết lập một tiêu chuẩn tuyến tính hóa giải tích mới có khả năng bù trừ sai số thay thế một chiều của phương pháp EQL kinh điển cho hệ nhiều bậc tự do?
  • Câu hỏi nghiên cứu 2 (Q2): Cấu trúc đại số của hệ đóng các phương trình xác định ma trận hệ số tuyến tính hóa tương đương ($M_e, C_e, K_e$) và ma trận thay thế ngược $D$ được hình thành như thế nào khi áp dụng tiêu chuẩn đối ngẫu?
  • Câu hỏi nghiên cứu 3 (Q3): Tiêu chuẩn đối ngẫu thể hiện ưu thế định lượng ra sao khi mở rộng cho các hệ liên tục phức tạp như dầm Euler-Bernoulli phi tuyến hình học mang khối lượng tập trung dưới kích động ồn trắng Gaussian?

Hệ thống giả thuyết nghiên cứu bao gồm:

  • Giả thuyết 1 (H1): Việc đưa vào bước thay thế ngược (đối ngẫu) hệ tuyến tính hóa bằng một hệ phi tuyến cùng lớp với hệ gốc, thông qua ma trận hệ số $D$ và phiếm hàm sai số kết hợp $E[e^T e] + \gamma E[\psi^T \psi] \to \min$, sẽ thiết lập trạng thái cân bằng năng lượng và giảm triệt để độ lệch thống kê bậc hai.
  • Giả thuyết 2 (H2): Tồn tại một hệ đại số tuyến tính đóng kín bậc $3n^2$ đối với các phần tử $k_{ij}^e, c_{ij}^e, d_{ij}$, cho phép kết hợp trực tiếp với ma trận hàm đáp ứng tần số $H(\omega)$ và ma trận mật độ phổ $S_Q(\omega)$ để hội tụ nghiệm dừng Gaussian.
  • Giả thuyết 3 (H3): Sai số tương đối của đáp ứng bình phương trung bình ($E[x_1^2], E[x_2^2], E[w_1^2]$) tính theo tiêu chuẩn đối ngẫu sẽ giảm từ 3 đến 5 lần so với phương pháp tuyến tính hóa kinh điển và vượt trội hơn phương pháp tiêu chuẩn năng lượng (Energy Criterion) của Zhang và đồng nghiệp (1990).

Phạm vi nghiên cứu (Scope) của luận án bao gồm các hệ rời rạc 2 bậc tự do (2-DOF) có độ cứng phi tuyến bậc ba (Duffing), bậc năm, hệ ghép cản và độ cứng phi tuyến, cùng hệ liên tục là dầm Euler-Bernoulli tựa giản đơn chịu kích động ngẫu nhiên dừng, xét từ mô hình đơn mode, 2 mode đến 5 mode ($M=5$), có xét đến ảnh hưởng của độ mảnh $L/h$, tỉ số khối lượng tập trung $M_0/m_0$ và vị trí đặt tải $d_0/L$.

Literature Review và Positioning

Lịch sử phát triển của các phương pháp phân tích động lực học ngẫu nhiên ghi nhận sự đóng góp mang tính nền tảng của phương trình Fokker-Planck-Kolmogorov (FPK) (Kolmogorov, 1931; Andronov, 1933; Caughey, 1963). Tuy nhiên, việc tìm nghiệm giải tích chính xác của phương trình FPK chỉ giới hạn trong một số lớp hệ hẹp thỏa mãn điều kiện thế năng dừng (Stationary Potential) hoặc cân bằng chi tiết (Detailed Balance) (Caughey & Dienes, 1961; Cai & Lin, 1988; Wen-Yao & Tong, 1992; Crandall, 2006).

Để khắc phục rào cản tính toán trên hệ MDOF, các trường phái tiếp cận xấp xỉ đã phân hóa rõ rệt:

  1. Trường phái Tuyến tính hóa thống kê kinh điển (Classical EQL): Atalik & Utku (1976) đã tổng quát hóa tiêu chuẩn sai số bình phương trung bình cho hệ MDOF bằng cách sử dụng toán tử đạo hàm riêng ma trận kỳ vọng gradient $\nabla \Phi(q, \dot{q})$, biến đổi phương trình phi tuyến thành hệ đại số ma trận tương đương: $$M_e = E\left[\frac{\partial \Phi}{\partial \ddot{q}}\right], \quad C_e = E\left[\frac{\partial \Phi}{\partial \dot{q}}\right], \quad K_e = E\left[\frac{\partial \Phi}{\partial q}\right]$$ Roberts & Spanos (1990) và Socha (2008) đã hệ thống hóa hoàn chỉnh phương pháp này. Điểm hạn chế chí tử là EQL cổ điển chỉ tối thiểu hóa phần dư phương trình vi phân $e = \Phi(q, \dot{q}) - C_e \dot{q} - K_e q$, dẫn đến việc đánh giá quá cao độ cứng tương đương và làm sai lệch đáng kể hàm mật độ xác suất (PDF) ở vùng đuôi phân bố.

  2. Trường phái Tiêu chuẩn Năng lượng (Energy-based EQL): Zhang, Elishakoff & Zhang (1990) và Elishakoff & Colombi (1993) đề xuất cực tiểu hóa sai số bình phương trung bình của thế năng $E[(V(q) - V_e(q))^2] \to \min$ hoặc hàm tiêu tán năng lượng $E[(W(\dot{q}) - W_e(\dot{q}))^2] \to \min$. Mặc dù cải thiện độ chính xác cho một số hệ 1-DOF, phương pháp năng lượng bộc lộ tranh luận lớn (Debates): không thể áp dụng đồng thời khi hàm phi tuyến chứa cả tọa độ vị trí và vận tốc ghép kênh, đồng thời việc tính kỳ vọng của đa thức thế năng bậc cao gặp khó khăn giải tích lớn trên hệ MDOF.

  3. Trường phái Tuyến tính hóa phi tham số (Non-parametric Linearization): Fujimura & Der Kiureghian (2007), Broccardo & Der Kiureghian (2015) sử dụng phương pháp bám điểm thiết kế trong không gian chuẩn Gaussian rời rạc hóa. Phương pháp này mô phỏng tốt phân bố non-Gaussian nhưng đòi hỏi tài nguyên tính toán lặp phi tuyến phức tạp, tiệm cận mô phỏng số hơn là công thức giải tích đóng.

Định vị của luận án: Luận án của Nguyễn Như Hiếu kế thừa tư tưởng điều chỉnh tuyến tính hóa của Nguyễn Đông Anh & Di Paola (1995), Elishakoff, Andrimasy & Dolley (2009), và bước phát triển ban đầu của tiêu chuẩn đối ngẫu cho hệ 1-DOF (Anh, 2010; Anh và cộng sự, 2012a). Luận án định vị một bước tiến đột phá: mở rộng hoàn chỉnh tiêu chuẩn đối ngẫu từ hệ 1-DOF sang hệ MDOF rời rạc và hệ liên tục chịu kích động ngẫu nhiên, thiết lập công thức giải tích tường minh giải quyết triệt để sự mất cân bằng cấu trúc trong phép tuyến tính hóa tương đương.

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

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

Luận án thách thức và mở rộng trực tiếp lý thuyết tuyến tính hóa tương đương của Caughey (1956) và Atalik & Utku (1976). Về mặt nhận thức luận, luận án chứng minh rằng việc tuyến tính hóa một chiều thuần túy ($Phi \to \text{Tuyến tính}$) gây ra sự thiên lệch hệ thống (systematic bias) trong việc ước lượng các mô-men thống kê bậc hai.

Mô hình lý thuyết đối ngẫu đề xuất một cơ chế bù trừ kép:

  • Bước thuận (Primal step): Thay thế vector lực phi tuyến $\Phi(q, \dot{q})$ bằng vector tuyến tính $C_e \dot{q} + K_e q$, phát sinh vector sai số: $$e = \Phi(q, \dot{q}) - C_e \dot{q} - K_e q$$
  • Bước nghịch/đối ngẫu (Dual step): Thay thế ngược vector tuyến tính $C_e \dot{q} + K_e q$ bằng một cấu trúc phi tuyến tương đương $D\Phi(q, \dot{q})$ cùng lớp với hệ gốc (với $D = [d_{ij}]_{n \times n}$ là ma trận hệ số thay thế trở lại), phát sinh vector sai số đối ngẫu: $$\psi = C_e \dot{q} + K_e q - D\Phi(q, \dot{q})$$

Phiếm hàm tối ưu hóa kết hợp đồng thời hai bước thay thế được thiết lập: $$J(C_e, K_e, D) = E[e^T e] + \gamma E[\psi^T \psi] \to \min_{c_{ij}^e, k_{ij}^e, d_{ij}}$$ trong đó $\gamma > 0$ là tham số điều chỉnh đối ngẫu (chọn $\gamma = 1$ cho trạng thái cân bằng chuẩn).

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

Khung phân tích của luận án tích hợp chặt chẽ 3 lý thuyết nền tảng: Lý thuyết quá trình ngẫu nhiên dừng (quá trình Gaussian, Wiener, đạo hàm tổng quát ồn trắng), Lý thuyết dao động cơ học phi tuyến nhiều bậc tự do, và Lý thuyết phổ đáp ứng tần số liên tục.

Từ điều kiện cực trị đạo hàm riêng của phiếm hàm: $$\frac{\partial J}{\partial c_{ij}^e} = 0, \quad \frac{\partial J}{\partial k_{ij}^e} = 0, \quad \frac{\partial J}{\partial d_{ij}} = 0 \quad (\forall i, j = 1, \dots, n)$$

Luận án đã thiết lập nên hệ phương trình ma trận đóng gồm $3n^2$ ẩn số: $$\begin{bmatrix} (1+\gamma) E[\dot{q}\dot{q}^T] & (1+\gamma) E[\dot{q}q^T] & -E[\dot{q}\Phi^T] - \gamma E[\dot{q}\Phi^T] \ (1+\gamma) E[q\dot{q}^T] & (1+\gamma) E[qq^T] & -E[q\Phi^T] - \gamma E[q\Phi^T] \ -\gamma E[\Phi\dot{q}^T] & -\gamma E[\Phi q^T] & \gamma E[\Phi\Phi^T] \end{bmatrix} \begin{bmatrix} (C_i^e)^T \ (K_i^e)^T \ (D_i)^T \end{bmatrix} = \begin{bmatrix} E[\dot{q}\Phi_i] \ E[q\Phi_i] \ 0 \end{bmatrix}$$

Điều kiện biên lý thuyết (Boundary conditions) được xác định rõ: phi tuyến phải là các hàm liên tục khả vi từng khúc, kích động đầu vào là quá trình ngẫu nhiên dừng Gaussian có kỳ vọng triệt tiêu $E[Q(t)] = 0$, ma trận mật độ phổ $S_Q(\omega)$ liên tục xác định dương, và ma trận hiệp phương sai đáp ứng $R_u$ không suy biến.

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

Thiết kế nghiên cứu

Luận án vận hành trên lập trường nhận thức luận thực chứng thực nghiệm kết hợp duy lý giải tích (Analytical Positivism). Thiết kế nghiên cứu là sự kết hợp chặt chẽ giữa 3 tầng phương pháp:

  1. Tầng nghiệm chính xác giải tích: Sử dụng phương trình vi phân đạo hàm riêng Parabolic Fokker-Planck-Kolmogorov (FPK) ở trạng thái dừng với điều kiện dòng xác suất triệt tiêu: $$G_i(x) = A_i(x)f(x) - \frac{1}{2} \sum_{j=1}^n \frac{\partial}{\partial x_j} [B_{ij}(x)f(x)] = 0$$ Nghiệm mật độ xác suất chính xác thu được dưới dạng: $$f(x) = \frac{1}{b(x)} \exp\left[ -2 \int_{a_0}^x \sum_{i=1}^n \frac{A_i(x)}{b(x)} dx_i \right]$$
  2. Tầng nghiệm xấp xỉ giải tích đối ngẫu: Thiết lập ma trận hàm đáp ứng biên độ – tần số $H(\omega)$: $$H(\omega) = \left[ -\omega^2 M + i\omega (C + C_e) + (K + K_e) \right]^{-1}$$ kết hợp tích phân Duhamel xác định ma trận mô-men bậc hai: $$E[qq^T] = \int_{-\infty}^{\infty} H(\omega) S_Q(\omega) H^T(-\omega) d\omega, \quad E[\dot{q}\dot{q}^T] = \int_{-\infty}^{\infty} \omega^2 H(\omega) S_Q(\omega) H^T(-\omega) d\omega$$
  3. Tầng mô phỏng số đối chứng (Monte-Carlo Simulation - MCS): Xây dựng thuật toán sinh mẫu kích động ngẫu nhiên dừng, giải tích phân từng bước phương trình vi phân ngẫu nhiên phi tuyến trên miền thời gian lớn để ước lượng kỳ vọng thực nghiệm với độ tin cậy $95%$.

Quy trình nghiên cứu và cấu trúc thuật toán

Quy trình giải bài toán tuyến tính hóa đối ngẫu được thực thi qua vòng lặp tự phối hợp (Self-consistent Iteration Scheme):

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

Những phát hiện then chốt với minh chứng định lượng

Dữ liệu tính toán số và mô phỏng từ Chương 4 của luận án cung cấp các bằng chứng thực nghiệm và giải tích vững chắc:

Mô hình khảo sát Tham số phi tuyến / Độ mảnh Nghiệm chính xác (CX) / MCS Tuyến tính hóa thông thường (TT) Tuyến tính hóa năng lượng (NL) Tiêu chuẩn đối ngẫu (DN) Mức cải thiện sai số của DN so với TT
Hệ 2-DOF phi tuyến bậc 3 ($x_1$) $\beta_1 = 1.0, \beta_2 = 0.0$ $E[x_1^2] = 0.3633$ $0.4185$ (Sai số $+15.19%$) $0.3921$ (Sai số $+7.92%$) $0.3702$ (Sai số $+1.89%$) Giảm sai số 8.0 lần
Hệ 2-DOF phi tuyến bậc 3 ($x_2$) $\beta_1 = 1.0, \beta_2 = 0.0$ $E[x_2^2] = 0.4421$ $0.4950$ (Sai số $+11.96%$) $0.4710$ (Sai số $+6.53%$) $0.4498$ (Sai số $+1.74%$) Giảm sai số 6.9 lần
Hệ 2-DOF phi tuyến bậc 5 ($x_1$) $\mu_1 = 2.0, \mu_2 = 0.0$ $E[x_1^2] = 0.2845$ $0.3520$ (Sai số $+23.72%$) $0.3210$ (Sai số $+12.83%$) $0.2930$ (Sai số $+2.98%$) Giảm sai số 7.9 lần
Dầm đơn mode ($w_1$) $1/R = 100, S_0 = 1$ $E[w_1^2] = 0.0412\text{ m}^2$ $0.0518\text{ m}^2$ ($+25.72%$) $0.0475\text{ m}^2$ ($+15.29%$) $0.0425\text{ m}^2$ ($+3.15%$) Giảm sai số 8.1 lần
Dầm 2 mode ($w_1$) $1/R = 200, S_0 = 1$ $E[w_1^2] = 0.0245\text{ m}^2$ $0.0321\text{ m}^2$ ($+31.02%$) $0.0289\text{ m}^2$ ($+17.95%$) $0.0254\text{ m}^2$ ($+3.67%$) Giảm sai số 8.4 lần
  1. Hiện tượng triệt tiêu sai số phi tuyến cao bậc: Khi bậc phi tuyến tăng từ bậc 3 lên bậc 5 và cường độ kích động ngẫu nhiên $S_0$ tăng cao, sai số của phương pháp EQL thông thường bùng nổ vượt ngưỡng $25-35%$. Tiêu chuẩn đối ngẫu duy trì độ lệch ổn định dưới $3.5%$, chứng minh tính bền vững giải tích (analytical robustness).
  2. Hiệu ứng phân phối năng lượng đa mode trên kết cấu dầm: Khảo sát dầm Euler-Bernoulli với $M=5$ mode dao động chỉ ra rằng mode dao động cơ bản ($m=1$) chiếm ưu thế tuyệt đối ($>85%$) trong tổng đáp ứng bình phương trung bình. Tỉ số đáp ứng giữa mode 1 và các mode cao ($m=2, 3, 4, 5$) tăng phi tuyến theo tham số độ mảnh $L/h$, khẳng định tiêu chuẩn đối ngẫu nắm bắt chính xác tương tác phi tuyến liên mode (inter-mode coupling).
  3. Ảnh hưởng của khối lượng tập trung và vị trí đặt: Khi tỉ số khối lượng $M_0/m_0$ tăng từ 1 lên 5 và vị trí đặt khối lượng $d_0/L$ dịch chuyển về tâm dầm ($d_0/L = 0.5$), đáp ứng dao động $E[w_1^2]$ suy giảm mạnh do quán tính cục bộ, đồng thời làm tăng tính phi tuyến hình học của dầm. Phương pháp đối ngẫu mô tả trùng khít đường cong biến thiên của mô phỏng Monte-Carlo trong khi EQL thông thường đánh giá thấp biên độ uốn.

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Mở rộng biên giới của cơ học ngẫu nhiên (Stochastic Mechanics), cung cấp một công cụ giải tích chuẩn mực mới thay thế cho quy trình tuyến tính hóa cổ điển Atalik-Utku vốn đã tồn tại hơn 4 thập kỷ mà chưa có sự đột phá căn bản về cấu trúc đại số.
  • Về mặt kỹ thuật công trình: Cho phép các kỹ sư kết cấu tính toán chính xác ứng suất và chuyển vị ngẫu nhiên của các công trình siêu cao tầng, giàn khoan ngoài khơi, cầu treo dây võng chịu tải trọng sóng và gió bão mà không cần tiêu tốn hàng trăm giờ chạy mô phỏng Monte-Carlo trên siêu máy tính.
  • Về thiết kế kháng chấn và giảm chấn: Đặt cơ sở tính toán tối ưu cho các thiết bị giảm chấn phi tuyến (Nonlinear Tuned Mass Dampers - TMD, Bouc-Wen hysteretic dampers) bảo vệ công trình trước thảm họa động đất.

Limitations và Future Research

Luận án nhận diện một cách khách quan các giới hạn nội tại:

  1. Giới hạn về mô hình tải kích động: Các nghiệm giải tích và công thức đạo hàm kỳ vọng được xây dựng trên giả thiết kích động ngẫu nhiên là quá trình dừng Gaussian ồn trắng (Gaussian White Noise). Trong thực tế, tải trọng gió và động đất mang bản chất ồn màu (Colored Noise) hoặc quá trình không dừng (Non-stationary Process) có tần số thay đổi theo thời gian.
  2. Giới hạn về tính chất vật liệu: Mô hình phi tuyến khảo sát tập trung vào các hàm giải tích liên tục (đa thức bậc ba, bậc năm, cản nhớt phi tuyến). Các hệ có tính chất trễ phi đàn hồi mạnh (Hysteresis) như mô hình Bouc-Wen có suy giảm độ cứng và trượt (pinching) chưa được bao hàm trực tiếp trong hệ đóng của luận án.
  3. Giới hạn số bậc tự do của nghiệm giải tích chính xác: Nghiệm chính xác từ phương trình FPK thế năng dừng chỉ thiết lập được cho hệ 2-DOF, việc đối chứng ở hệ nhiều bậc tự do cao hơn ($n > 2$) hoàn toàn phải phụ thuộc vào mô phỏng số Monte-Carlo.

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

  • Mở rộng tiêu chuẩn đối ngẫu cho hệ phi tuyến chịu kích động ồn màu (sử dụng phương pháp lọc định hình trạng thái - State-space Shaping Filter).
  • Phát triển công thức đối ngẫu cho các kết cấu thông minh chịu kích động không dừng (Wavelet-based EQL).
  • Tích hợp tiêu chuẩn đối ngẫu vào phần mềm phần tử hữu hạn thương mại (ANSYS, ABAQUS, SAP2000) qua các module lập trình UMAT/VUMAT để mô phỏng vỏ mỏng composite và kết cấu không gian phức tạp.

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

  • Tác động học thuật quốc tế: Công trình đặt nền móng cho chuỗi ấn phẩm chuyên sâu trên các tạp chí hàng đầu danh mục ISI/Scopus Q1 trong lĩnh vực cơ học ứng dụng và kỹ thuật công trình (như Journal of Sound and Vibration, Nonlinear Dynamics, Probabilistic Engineering Mechanics).
  • Chuyển giao công nghệ R&D: Ứng dụng trực tiếp trong các đề tài nghiên cứu cơ bản cấp quốc gia do Quỹ Phát triển Khoa học và Công nghệ Quốc gia (NAFOSTED) tài trợ, nâng cao năng lực tính toán cơ học kết cấu tại Viện Cơ học – Viện Hàn lâm Khoa học và Công nghệ Việt Nam.
  • Tối ưu hóa chi phí xây dựng: Giảm biên độ an toàn giả tạo (over-conservative safety factors) do sai số của phương pháp tính cũ gây ra, giúp tiết kiệm $8-12%$ khối lượng vật liệu thép và bê tông trong các công trình kết cấu nhạy cảm với rung động ngẫu nhiên.

Đố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, Giao thông: Tiếp cận hệ thống phương pháp luận giải tích chặt chẽ, mở ra các đề tài luận án mới về dao động ngẫu nhiên phi tuyến.
  • Các giáo sư, chuyên gia cơ học lý thuyết: Sử dụng hệ phương trình đóng đối ngẫu $3n^2$ như một công cụ chuẩn mực để so sánh và kiểm chuẩn các phương pháp xấp xỉ số mới.
  • Kỹ sư trưởng R&D tại các tập đoàn tư vấn thiết kế kết cấu: Sở hữu thuật toán giải tích tốc độ cao để đánh giá độ tin cậy và tuổi thọ mỏi của công trình cầu đường, công trình biển và nhà cao tầng.
  • Cơ quan quản lý và xây dựng quy chuẩn xây dựng: Bổ sung cơ sở khoa học để cập nhật các tiêu chuẩn thiết kế kháng chấn quốc gia tiệm cận tiêu chuẩn quốc tế (Eurocode 8, ASCE 7).

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 phát hiện và khắc phục sự mất cân bằng cấu trúc của phương pháp tuyến tính hóa tương đương thông thường bằng việc thiết lập tiêu chuẩn tối ưu hóa kết hợp hai bước thuận - nghịch (Dual Criterion): $$J = E[e^T e] + \gamma E[\psi^T \psi] \to \min$$ Luận án đã mở rộng thành công lý thuyết tuyến tính hóa thống kê của Caughey (1956) và Atalik & Utku (1976) sang hệ nhiều bậc tự do và hệ liên tục dầm Euler-Bernoulli phi tuyến, chuyển hóa bài toán tối ưu phi tuyến thành hệ đại số tuyến tính đóng kín giải được tường minh.

2. Sự đổi mới về mặt phương pháp luận khi so sánh với các nghiên cứu quốc tế?

So với phương pháp EQL kinh điển (Atalik & Utku, 1976) vốn có sai số tăng vọt khi hệ phi tuyến mạnh, và phương pháp tiêu chuẩn năng lượng (Zhang và cộng sự, 1990) vốn bị bế tắc khi hàm phi tuyến chứa đồng thời cả vị trí và vận tốc, phương pháp đối ngẫu của luận án tổng quát hóa được toàn bộ các dạng phi tuyến ghép kênh $\Phi(q, \dot{q})$, đồng thời hạ thấp sai số đáp ứng bình phương trung bình xuống dưới $3.5%$ mà vẫn duy trì dạng đóng giải tích thuần túy.

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

Phát hiện bất ngờ nhất là tính ổn định phi thường của sai số xấp xỉ khi bậc phi tuyến tăng từ bậc 3 lên bậc 5. Trong khi sai số của EQL thông thường tăng vọt từ $15.19%$ lên $23.72%$ đối với hệ 2-DOF, sai số của tiêu chuẩn đối ngẫu chỉ biến thiên nhẹ từ $1.89%$ lên $2.98%$. Điều này chứng minh rằng ma trận thay thế ngược $D$ đã triệt tiêu hiệu quả các số hạng nhiễu loạn bậc cao của phân bố xác suất non-Gaussian.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án cung cấp tường minh toàn bộ các ma trận đại số, thuật toán lặp số, công thức giải tích tích phân Gauss-Hermite để tính kỳ vọng hàm gradient, công thức giải tích ma trận truyền tần số $H(\omega)$ và tích phân phổ Duhamel, cho phép các nhà nghiên cứu lập trình tái lập chính xác $100%$ kết quả trên phần mềm MATLAB, Mathematica hoặc Python.

5. Tầm nhìn chương trình nghiên cứu 10 năm tiếp theo được định hình thế nào?

Chương trình 10 năm hướng tới: (1) Mở rộng tiêu chuẩn đối ngẫu cho hệ động lực học ngẫu nhiên không dừng và ồn phân đoạn (Fractional Brownian Motion); (2) Ứng dụng giải tích độ tin cậy động lực học (Dynamic Reliability) và bài toán vượt ngưỡng an toàn đầu tiên (First-passage Failure); (3) Tích hợp học máy giải tích (Physics-Informed Neural Networks - PINNs) kết hợp tiêu chuẩn đối ngẫu để tối ưu hóa kết cấu tự thích ứng.

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 các mục tiêu nghiên cứu đặt ra, tạo nên dấu ấn học thuật quan trọng trong chuyên ngành Cơ học vật rắn:

  1. Xây dựng thành công Tiêu chuẩn đối ngẫu cho hệ MDOF: Thiết lập hệ đóng gồm $3n^2$ phương trình đại số xác định chính xác các ma trận hệ số tuyến tính hóa tương đương $K_e, C_e$ và ma trận thay thế ngược $D$.
  2. Khảo sát toàn diện các hệ động lực phi tuyến tiêu biểu: Giải quyết trọn vẹn bài toán đáp ứng ngẫu nhiên cho hệ 2-DOF phi tuyến bậc ba, bậc năm, hệ có cản và độ cứng phi tuyến, đạt độ chính xác tiệm cận nghiệm chính xác FPK.
  3. Tiên phong mở rộng cho hệ liên tục dầm Euler-Bernoulli: Phân tích thành công dao động ngẫu nhiên của dầm phi tuyến hình học mang khối lượng tập trung từ đơn mode đến đa mode ($M=5$), làm sáng tỏ cơ chế truyền năng lượng giữa các mode dao động.
  4. Chứng minh định lượng sự vượt trội của phương pháp: Dữ liệu tính toán chỉ ra tiêu chuẩn đối ngẫu giúp giảm sai số tương đối từ 3 đến 8 lần so với phương pháp tuyến tính hóa kinh điển và vượt trội hơn phương pháp năng lượng.
  5. Mở ra 3 nhánh nghiên cứu mới trong tương lai: Tạo tiền đề phát triển các phương pháp giải tích cho hệ chịu kích động ồn màu phi Gaussian, hệ trễ Bouc-Wen phức tạp, và công nghệ điều khiển giảm chấn bán tích cực cho công trình hiện đại.
  6. Giá trị thực tiễn và di sản học thuật: Đóng góp trực tiếp vào kho tàng cơ học tính toán Việt Nam, cung cấp công cụ đắc lực phục vụ thiết kế an toàn và tối ưu hóa kết cấu công trình chịu tải trọng thiên tai ngẫu nhiên.