Tổng quan về luận án

Nghiên cứu về các hệ thống động lực phi tuyến luôn là tâm điểm của cơ học vật rắn biến dạng, vật lý kỹ thuật và khoa học cấu trúc hiện đại. Luận án tiến sĩ "Phân tích dao động phi tuyến bằng cách tiếp cận trung bình có trọng số" của tác giả Đặng Văn Hiếu (Chuyên ngành: Cơ học vật rắn, Mã số: 9440107; thực hiệ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 dưới sự hướng dẫn của PGS. Ninh Quang Hải và TS. Dương Thế Hùng) 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 giải tích trong phân tích các hệ dao động phi tuyến mạnh tiền định mà không làm bùng nổ độ phức tạp tính toán.

Trong cơ học phi tuyến, phương pháp tuyến tính hóa tương đương (Equivalent Linearization - EQL) do Caughey (1963) phát triển trên nền tảng của Krylov và Bogoliubov (1937) được ứng dụng rộng rãi nhờ cấu trúc thuật toán tường minh. Tuy nhiên, một research gap mang tính căn bản của EQL cổ điển là việc sử dụng toán tử trung bình thời gian cổ điển trên một chu kỳ dao động $T = 2\pi/\omega$:

$$\langle f(t) \rangle_c = \frac{1}{T} \int_0^T f(t) dt$$

Toán tử này làm triệt tiêu toàn bộ giá trị trung bình của các hàm điều hòa cơ bản: $\langle \cos(n\omega t) \rangle_c = \langle \sin(n\omega t) \rangle_c = 0$ với mọi $n \in \mathbb{N}^*$. Khi các hệ dao động xuất hiện tính phi tuyến mạnh (bậc cao, số mũ hữu tỉ, thế năng giếng đôi hoặc có bước nhảy gián đoạn), toán tử trung bình cổ điển làm thất thoát nghiêm trọng thông tin cấu trúc pha và biên độ, dẫn đến sai số của nghiệm xấp xỉ bậc một vượt quá $15% - 30%$, thậm chí sai lệch hoàn toàn về mặt định tính.

Để giải quyết khoảng trống học thuật này, luận án xây dựng hệ phương pháp luận dựa trên 4 câu hỏi nghiên cứu và 4 giả thuyết khoa học tương ứng:

  • RQ1: Làm thế nào để thiết lập một toán tử lấy trung bình mới bảo toàn được thông tin dạng sóng của các hàm điều hòa bậc cao trong quá trình tuyến tính hóa tương đương?
    • H1: Toán tử trung bình có trọng số sử dụng hàm phân bố tắt dần theo thời gian $h(t) = s^2 \omega^2 t e^{-s\omega t}$ sẽ giữ lại các đặc trưng động lực học chủ đạo trong nửa chu kỳ đầu, triệt tiêu sự mất mát thông tin của các số hạng điều hòa.
  • RQ2: Bản chất toán học của toán tử trung bình có trọng số liên hệ như thế nào với các phép biến đổi tích phân kinh điển?
    • H2: Giá trị trung bình có trọng số của hàm tuần hoàn thực chất là một đẳng cấu trực tiếp với ảnh của hàm gốc qua phép biến đổi Laplace, cho phép đại số hóa hoàn toàn quá trình giải tích.
  • RQ3: Mức độ chính xác và phạm vi ứng dụng của phương pháp EQL kết hợp trung bình có trọng số đối với các hệ phi tuyến 1 bậc tự do có độ phức tạp cao (Duffing bậc 3 đến bậc 13, dao động giếng đôi, số mũ hữu tỉ, hệ gián đoạn) đạt ngưỡng nào so với nghiệm số Runge-Kutta?
    • H3: Phương pháp đề xuất cho phép thiết lập biểu thức giải tích dạng đóng (closed-form analytical solutions) của tần số phi tuyến $\omega_{NL}$ với sai số tương đối giảm xuống dưới $1% - 2%$ trên toàn dải biên độ ban đầu $A$.
  • RQ4: Hiệu ứng phụ thuộc kích thước (size-dependent effects) của các cấu trúc micro/nano (dầm micro trên nền đàn hồi 3 tham số, dầm nano MEMS/NEMS chịu điện trường tĩnh điện) được biểu diễn chính xác như thế nào qua khung phân tích giải tích mới?
    • H4: Việc tích hợp EQL trung bình có trọng số với Lý thuyết ứng suất cặp sửa đổi (Modified Couple Stress Theory - MCST) và Lý thuyết độ dốc biến dạng phi cục bộ (Nonlocal Strain Gradient Theory - NSGT) sẽ làm sáng tỏ sự cạnh tranh giữa hiệu ứng làm cứng (stiffness hardening) và làm mềm (stiffness softening) của kết cấu vi mô/nano.

Khung lý thuyết của luận án tích hợp chặt chẽ giữa: (1) Cơ sở tuyến tính hóa tương đương của Caughey (1963); (2) Lý thuyết trung bình có trọng số do GS. Nguyễn Đông Anh đề xuất (2015); (3) Lý thuyết ứng suất cặp sửa đổi của Yang và cộng sự (2002); (4) Lý thuyết độ dốc biến dạng phi cục bộ của Lim và cộng sự (2015). Đóng góp mang tính đột phá của công trình là thiết lập thành công quy trình giải tích dạng đóng cho 5 lớp dao động phi tuyến 1 bậc tự do điển hình và 2 bài toán hệ liên tục phức tạp ở quy mô micro/nano, nâng cao độ chính xác định lượng vượt trội so với các phương pháp giải tích đương đại (HPM, EBM, Variational Approach).

Literature Review và Positioning

Bản đồ tiến trình nghiên cứu dao động phi tuyến thế giới ghi nhận nhiều trường phái giải tích tiệm cận và biến phân:

                            TIẾN TRÌNH LÝ THUYẾT KHÁCH QUAN
                            
  [1937] Krylov - Bogoliubov: Tuyến tính hóa tương đương tiền định
  [1963] Caughey: Tuyến tính hóa tương đương ngẫu nhiên (Toán tử cổ điển ⟨·⟩_c)
  [1991-1995] Các biến thể EQL    [2000-2010] Phương pháp năng    [2002-2015] Cơ học Micro/Nano
  - Zhang (1991): Thế năng          lượng & Nhiễu biến phân       - Yang et al. (2002): MCST
  - Elishakoff (1993): Từng phần  - He J.H. (2002): EBM           - Lim et al. (2015): NSGT
  - Anh & Di Paola (1995):        - Beléndez (2009): HPM          - Şimşek (2010): Variational
    Regulated EQL                 - Sedighi (2012): Hamiltonian     Beam Dynamics
                 [2015] GS. Nguyễn Đông Anh: Toán tử Trung bình có trọng số ⟨·⟩_w
                 [2021] Luận án Đặng Văn Hiếu: KHUNG GIẢI TÍCH ĐÓNG TOÀN DIỆN
                 (EQL-WA + Laplace Transform + MCST/NSGT cho Hệ 1-DOF & Micro/Nano)
  1. Trường phái Nhiễu tiệm cận và Nhiễu đồng luân: Xuất phát từ lý thuyết nhiễu cổ điển (Nayfeh & Mook, 1979), Ji-Huan He (2000, 2006) đã mở rộng với Phương pháp nhiễu đồng luân (Homotopy Perturbation Method - HPM) và Phương pháp tham số hóa (Parameterized Perturbation Method). Beléndez và cộng sự (2009, 2011) đã phát triển các xấp xỉ bậc cao cho dao động phi tuyến mạnh. Tuy nhiên, hạn chế cố hữu của trường phái này là quy trình tính toán đại số cực kỳ nặng nề khi hàm phục hồi có chứa số mũ hữu tỉ hoặc phân thức phức tạp.
  2. Trường phái Cân bằng điều hòa và Phương pháp năng lượng: Phương pháp cân bằng điều hòa (Harmonic Balance Method - HBM) được Ronald E. Mickens (2001, 2005, 2010) phát triển mạnh mẽ cho các hệ có số mũ phân số. Bên cạnh đó, Ji-Huan He (2002, 2007, 2010) đề xuất Phương pháp cân bằng năng lượng (Energy Balance Method - EBM), Phương pháp biến phân (Variational Approach - VA) và Phương pháp tiếp cận Hamiltonian. Dù mang lại kết quả tốt ở xấp xỉ bậc nhất, các phương pháp này phụ thuộc vào việc lựa chọn điểm định vị (location point) hoặc đòi hỏi giải hệ phương trình phi tuyến siêu việt phức tạp khi xét đến các bậc cao (Saadatnia et al., 2015; Yildirim et al., 2012).
  3. Trường phái Tuyến tính hóa tương đương và các cải tiến: Nhằm khắc phục hạn chế sai số của Caughey (1963), nhiều cải tiến đã ra đời như: Tiêu chuẩn cực tiểu sai số thế năng (Zhang et al., 1991), Tuyến tính hóa từng phần (Elishakoff & Cai, 1993), Tuyến tính hóa điều chỉnh (Anh & Di Paola, 1995; Elishakoff et al., 1996), Tuyến tính hóa có đuôi (Fujimura & Der Kiureghian, 2007), Tiêu chuẩn đối ngẫu và đối ngẫu có trọng số (Nguyen Dong Anh, 2010; Nguyen Dong Anh & Nguyen Ngoc Linh, 2013). Dù vậy, phần lớn các công trình này tập trung vào kích động ngẫu nhiên, để ngỏ việc hoàn thiện lời giải giải tích dạng đóng cho các hệ tiền định phi tuyến phức tạp.

Về phương diện cơ học cấu trúc vi mô, các công trình quốc tế của Şimşek (2010, 2016), Reddy (2011), Akgöz & Civalek (2013), Sedighi (2014) và Fu et al. (2011) đã ứng dụng các lý thuyết đàn hồi phi cổ điển để mô hình hóa dầm micro/nano. Tuy nhiên, phần lớn các công trình sử dụng phương pháp số hoặc biến phân đơn thuần, chưa khảo sát toàn diện ảnh hưởng liên đới giữa tỉ số độ cứng chống uốn, tỉ số độ mảnh và nền đàn hồi phi tuyến 3 thông số Winkler-Pasternak-$K_{NL}$, cũng như hiện tượng mất ổn định tĩnh điện (pull-in instability) kết hợp lực dọc trục dưới tác động đồng thời của hiệu ứng phi cục bộ và độ dốc biến dạng.

Luận án của Đặng Văn Hiếu định vị chính xác tại giao điểm này: kế thừa toán tử trung bình có trọng số của GS. Nguyễn Đông Anh (2015), phát triển thành một phương pháp giải tích hoàn chỉnh, kết nối với biến đổi Laplace để giải quyết triệt để bài toán dao động phi tuyến tiền định từ 1 bậc tự do đến các cấu trúc micro/nano.

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

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

Luận án đã mở rộng lý thuyết tuyến tính hóa tương đương của Krylov-Bogoliubov và Caughey bằng việc chính thức xác lập cơ sở toán học vững chắc cho toán tử trung bình có trọng số trong miền thời gian thực:

"Giá trị trung bình có trọng số của một hàm phụ thuộc thời gian $x(t)$ được xác định bởi: $$\langle x(t) \rangle_w = \int_0^\infty h(t) x(t) dt$$ với hàm hệ số trọng số tổ hợp: $$h(t) = s^2 \omega^2 t e^{-s\omega t}, \quad s > 0$$ trong đó $s$ là tham số điều chỉnh." (Trích Chương 2, Luận án)

       Hàm trọng số h(t)
          0       t_max = 1/(sω)
                  (= T / 2πs)

Bản chất của hàm trọng số $h(t)$ là tích của hàm tăng cơ bản $O(t) = t$ và hàm giảm cơ bản $P(t) = e^{-s\omega t}$. Hàm đạt cực đại tại thời điểm:

$$t_{\max} = \frac{1}{s\omega} = \frac{T}{2\pi s}$$

Khi lựa chọn $s = n/2$ ($n \in \mathbb{N}^*$), điểm tập trung trọng số sẽ rơi vào các phân vị chu kỳ cụ thể: $n=1 \Rightarrow t_{\max} = T/\pi$; $n=2 \Rightarrow t_{\max} = T/2$ (nửa chu kỳ); $n=4 \Rightarrow t_{\max} = T/4$ (một phần tư chu kỳ).

Đặc biệt, luận án chứng minh định lý về tính bao hàm và bảo toàn:

  • Khi $s \to 0$, toán tử trung bình có trọng số tiệm cận chính xác về toán tử trung bình cổ điển: $\lim_{s \to 0} \langle x(t) \rangle_w = \langle x(t) \rangle_c$.
  • Khi $h(t) = f_X(x)$ (hàm mật độ xác suất), toán tử chuyển hóa thành kỳ vọng toán học kinh điển: $\langle x(t) \rangle_w = E(X) = \int_{-\infty}^\infty x f_X(x) dx$.

Đóng góp lý thuyết mang tính bước ngoặt là việc phát hiện và chứng minh liên hệ toán học giữa trung bình có trọng số và biến đổi Laplace:

"Giá trị trung bình có trọng số của hàm tuần hoàn $x(t)$ thực chất là ảnh của hàm $s^2 t x(t)$ qua phép biến đổi Laplace: $$\langle x(t) \rangle_w = \int_0^\infty s^2 \omega^2 t e^{-s\omega t} x(t) dt = s^2 \int_0^\infty \tau e^{-s\tau} x(\tau) d\tau = s^2 \bar{X}(s)$$ với $\tau = \omega t$." (Trích Chương 2, Luận án)

Khám phá này chuyển hóa toàn bộ các tích phân phức tạp của hàm lượng giác bậc cao thành các phép đại số hữu tỉ đơn giản thông qua bảng Laplace:

$$\langle \cos(n\omega t) \rangle_w = \frac{s^2(s^2 - n^2)}{(s^2 + n^2)^2}, \quad \langle \sin(n\omega t) \rangle_w = \frac{2s^3 n}{(s^2 + n^2)^2}$$

$$\langle \cos^2(\omega t) \rangle_w = \frac{s^4 + 2s^2 + 8}{(s^2 + 4)^2}, \quad \langle \cos^4(\omega t) \rangle_w = \frac{s^8 + 28s^6 + 248s^4 + 416s^2 + 1536}{(s^2 + 4)^2 (s^2 + 16)^2}$$

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

Khung phân tích của luận án kết hợp 3 trụ cột: (1) Cơ học giải tích biến phân; (2) Toán tử trung bình có trọng số Laplace hóa; (3) Cơ học môi trường liên tục bậc cao phụ thuộc kích thước.

                      KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN
                      
     [HỆ RỜI RẠC 1 BẬC TỰ DO]                 [HỆ LIÊN TỤC MICRO / NANO]
     - Duffing bậc 3, 5, 7, ..., 13           - Dầm Micro: MCST (Yang et al.)
     - Duffing - Điều hòa                     - Dầm Nano: NSGT (Lim et al.)
     - Thế năng giếng đôi                     - Chiếu Galerkin qua hàm dạng W(x)
     - Số mũ hữu tỉ, Gián đoạn                - Rời rạc hóa về phương trình q̈(t)

Điều kiện biên và phạm vi áp dụng (boundary conditions): Phương pháp áp dụng tối ưu cho các hệ dao động bảo toàn tiền định không cản, biên độ dao động từ nhỏ đến rất lớn ($A \in (0, 100]$), tham số điều chỉnh trọng số $s \in [0.5, 2.0]$.

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ủ thế giới quan thực chứng (positivism) và cơ học duy lý tiền định. Thiết kế nghiên cứu kết hợp chặt chẽ giữa suy diễn giải tích giải pháp đóng (closed-form analytical deduction) và kiểm chứng số vi phân chính xác cao (high-precision numerical benchmarking).

Cấp độ phân tích Đối tượng nghiên cứu Mô hình cơ học / Phương trình toán Công cụ phân tích
Cấp độ 1: 1-DOF Dao động Duffing bậc 3, 5, 7, 9, 11, 13 $\ddot{X} + \sum_{k=1}^n \alpha_{2k-1} X^{2k-1} = 0$ EQL + Trung bình có trọng số
Cấp độ 1: 1-DOF Dao động Duffing - Điều hòa $\ddot{X} + \frac{\alpha X + \beta X^3}{1 + \gamma X^2} = 0$ EQL + Trung bình có trọng số
Cấp độ 1: 1-DOF Dao động thế năng giếng đôi $\ddot{X} - X + X^3 = 0$ ($A \ge \sqrt{2}$ và $1 < A < \sqrt{2}$) EQL + Tích phân Jacobi Elliptic
Cấp độ 1: 1-DOF Hệ số mũ hữu tỉ & Gián đoạn $\ddot{X} + c_1 X^{1/3} + c_2 X^{5/3} = 0$; $\ddot{X} + \text{sgn}(X) = 0$ EQL + Khai triển chuỗi trọng số
Cấp độ 2: Micro Dầm Micro trên nền 3 tham số MCST + Euler-Bernoulli + Pasternak Galerkin + EQL Trọng số
Cấp độ 2: Nano Dầm Nano MEMS/NEMS tĩnh điện NSGT + Lực tĩnh điện phi tuyến + Lực nén $P_0$ Phép biến phân + EQL Trọng số

Quy trình nghiên cứu rigorous

Quy trình giải tích được chuẩn hóa qua 5 giai đoạn nghiêm ngặt:

  1. Thiết lập mô hình vật lý: Xây dựng phương trình vi phân chuyển động từ nguyên lý Hamilton và biến phân năng lượng.
  2. Thiết lập quan hệ tuyến tính hóa: Giả thiết dạng nghiệm bậc nhất $X(t) = A\cos(\omega t)$ (với $\dot{X}(0)=0, X(0)=A$) và đưa vào phương trình sai số $e(X, \dot{X}) = g(X, \dot{X}) - k_{eq} X$.
  3. Cực tiểu hóa sai số bình phương trung bình có trọng số: $$\frac{\partial \langle e^2 \rangle_w}{\partial k_{eq}} = 0 \implies k_{eq} = \frac{\langle X g(X, \dot{X}) \rangle_w}{\langle X^2 \rangle_w}$$
  4. Đại số hóa tích phân Laplace: Chuyển đổi các số hạng $\langle X^{2k} \rangle_w$ sang miền Laplace $s$, triệt tiêu sự phụ thuộc thời gian và thu được biểu thức tường minh của tần số dao động phi tuyến $\omega_{NL} = \sqrt{k_{eq}}$.
  5. Kiểm chứng độ tin cậy tam giác (Triangulation): So sánh đối chiếu trực tiếp giữa:
    • Nghiệm giải tích đề xuất (EQL-WA).
    • Nghiệm giải tích cổ điển (EQL cổ điển, HPM, EBM, Variational Approach).
    • Nghiệm số chính xác cao giải bằng thuật toán Runge-Kutta bậc 4/5 (RK45) trong môi trường MATLAB.
       QUY TRÌNH KIỂM CHỨNG TAM GIÁC (TRIANGULATION PROTOCOL)

Data và phân tích

Bộ thông số mô phỏng đại diện trong luận án bao gồm:

  • Dao động Duffing bậc cao: Bậc 7 ($\alpha_1=1, \alpha_3=10, \alpha_5=10, \alpha_7=10$), Bậc 9, Bậc 11, Bậc 13 với dải tham số $\alpha_i \in [1, 100]$ và biên độ $A \in [0.1, 10]$.
  • Dầm Micro: Tỉ số độ mảnh $S = L/r \in [20, 100]$, tham số tỉ lệ chiều dài vật liệu $\alpha = l/h \in [0, 1]$, độ cứng chống uốn $\theta = 6$, thông số nền đàn hồi không thứ nguyên Winkler $K_L \in [0, 100]$, Pasternak $K_P \in [0, 100]$, phi tuyến $K_{NL} \in [0, 100]$.
  • Dầm Nano MEMS/NEMS: Tham số phi cục bộ $e_0 a \in [0, 2]\text{ nm}$, tham số chiều dài vật liệu $l \in [0, 1]\text{ nm}$, lực nén dọc trục $P_0/P_{cr} \in [0, 0.8]$, điện áp tĩnh điện một chiều $V_0/V_{cr} \in [0, 0.9]$.

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

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

        SO SÁNH SAI SỐ TƯƠNG ĐỐI (%) GIỮA CÁC PHƯƠNG PHÁP GIẢI TÍCH
        (Hệ Duffing bậc cao & Dao động giếng đôi phi tuyến mạnh)
        
  Sai số (%)
    0.1          1.0                 10.0  Biên độ dao động (A)
  1. Khắc phục triệt để sai số của EQL cổ điển trong hệ Duffing bậc cao: Với hệ Duffing bậc 7, khi biên độ $A$ tăng từ $0.1$ đến $10$, sai số tương đối của phương pháp EQL cổ điển vượt quá $20%$, trong khi phương pháp EQL kết hợp trung bình có trọng số với $s=1$ duy trì sai số dưới $0.85%$ so với nghiệm Runge-Kutta. Với Duffing bậc 13 ($\alpha_1=1, \alpha_3=10, \alpha_5=10, \alpha_7=10, \alpha_9=5, \alpha_{11}=10, \alpha_{13}=20$), tần số giải tích thu được khớp chính xác với nghiệm số ở 4 chữ số thập phân có nghĩa.
  2. Giải quyết bài toán phi tuyến thế năng giếng đôi: Trong miền biên độ $A \ge \sqrt{2}$ và đặc biệt là miền chuyển tiếp nhạy cảm $1 < A < \sqrt{2}$, nơi phương pháp tuyến tính hóa cổ điển hoàn toàn thất bại do rẽ nhánh nghiệm (bifurcation), nghiệm giải tích của luận án bám sát nghiệm chính xác elliptic với sai số cực đại chỉ $1.42%$.
  3. Quy luật chi phối của hiệu ứng kích thước lên dầm Micro: Khi tham số tỉ lệ chiều dài vật liệu $\alpha = l/h$ tăng từ $0$ lên $1.0$, tần số dao động riêng phi tuyến của dầm micro tăng hơn $180%$ đối với cả hai điều kiện biên (tựa bản lề và ngàm hai đầu). Tỉ số tần số $\omega_{NL}/\omega_L$ giảm khi độ cứng nền Pasternak ($K_P$) tăng, chứng minh vai trò cản trở chuyển vị phi tuyến của lớp cắt đàn hồi.
  4. Sự triệt tiêu và cạnh tranh kích thước trong cấu trúc Nano NEMS: Luận án phát hiện cơ chế tương tác đối nghịch: tham số phi cục bộ $e_0 a$ gây ra hiệu ứng làm mềm kết cấu (softening effect - giảm tần số tự nhiên), trong khi tham số độ dốc biến dạng $l$ gây ra hiệu ứng làm cứng (hardening effect - tăng tần số tự nhiên). Khi điện thế tác dụng $V_0$ tiệm cận điện thế tới hạn $V_{cr}$ (pull-in voltage), tần số phi tuyến sụt giảm đột ngột về 0, báo hiệu sự mất ổn định động học.
       CƠ CHẾ CẠNH TRANH ĐỘ CỨNG Ở CẤP ĐỘ NANO TRONG MEMS/NEMS
       
         Tham số phi cục bộ (e_0 a)           Tham số chiều dài vật liệu (l)
         [Hiệu ứng phi cục bộ vi mô]          [Hiệu ứng độ dốc biến dạng]
             LÀM MỀM KẾT CẤU                       LÀM CỨNG KẾT CỨNG
            (Stiffness Softening)                 (Stiffness Hardening)
          Giảm tần số riêng (ω_NL ↓)            Tăng tần số riêng (ω_NL ↑)
                  ĐIỀU KHIỂN ĐIỆN ÁP TĨNH ĐIỆN & ĐIỂM SẬP PULL-IN
                  - Dự báo chính xác điện áp tới hạn V_cr
                  - Tối ưu hóa dải thông cảm biến cộng hưởng NEMS

Implications đa chiều

  • Về mặt lý thuyết: Thiết lập chuẩn tắc mới cho việc áp dụng phương pháp tuyến tính hóa tương đương trong cơ học phi tuyến tiền định; mở rộng biên giới của lý thuyết đàn hồi phi cổ điển khi kết hợp với nền đàn hồi phi tuyến.
  • Về phương pháp luận: Cung cấp quy trình đại số hóa tính toán giải tích thông qua bảng biến đổi Laplace của hàm trọng số, mở ra khả năng tự động hóa giải tích trên các phần mềm đại số máy tính (Maple, Mathematica).
  • Về kỹ thuật ứng dụng: Cung cấp công thức dạng đóng giúp các kỹ sư thiết kế cảm biến vi cơ điện tử MEMS/NEMS, đầu dò kính hiển vi lực nguyên tử (AFM) có thể tính toán tức thời tần số cộng hưởng và ngưỡng mất ổn định pull-in mà không cần chạy mô phỏng phần tử hữu hạn (FEM) phi tuyến đắt đỏ.

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn học thuật:

  1. Phạm vi hệ động lực: Nghiên cứu chỉ tập trung vào các hệ dao động tự do không cản tiền định ($F(t)=0, c=0$). Các hiện tượng phi tuyến phức tạp như chu trình giới hạn (limit cycles) của hệ tự chấn (Van der Pol), dao động hỗn loạn (chaos) dưới kích động điều hòa đa tần chưa được bao hàm.
  2. Kỹ thuật rời rạc hóa: Trong bài toán dầm liên tục, luận án sử dụng phép chiếu Galerkin một mode đơn lẻ ($N=1$). Xấp xỉ này có độ chính xác cao cho mode cơ bản nhưng chưa mô tả được hiện tượng truyền năng lượng giữa các mode (internal resonance).
  3. Quy tắc xác định tham số $s$: Việc lựa chọn tham số điều chỉnh $s$ ($s=1, s=0.5$ hoặc $s=n/2$) chủ yếu dựa trên các đặc tính cực trị hình học của hàm $h(t)$ và kiểm chứng thực nghiệm số, chưa xây dựng được thuật toán biến phân tối ưu hóa tự động cho $s$ đối với từng dạng phi tuyến bất kỳ.

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  • Mở rộng phương pháp EQL trung bình có trọng số cho hệ dao động ngẫu nhiên chịu kích động ồn màu (colored noise) và kích động tham số.
  • Phân tích tương tác đa mode và cộng hưởng nội trong ống nano carbon truyền dẫn chất lỏng (CNT conveying fluid).
  • Tích hợp cản phi tuyến fractional-order (đạo hàm cấp phân số) vào khung trung bình có trọng số.
  • Phát triển phần mềm nguồn mở tự động hóa việc xuất biểu thức tần số giải tích cho các hệ cơ điện tử đa vật thể.

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

                               CÂY TÁC ĐỘNG HỌC THUẬT VÀ CÔNG NGHỆ
                               
                                   [ ĐỔI MỚI CƠ HỌC TÍNH TOÁN ]
 [HỌC THUẬT QUỐC TẾ]                   [CÔNG NGHỆ BÁN DẪN MEMS]                [GIẢNG DẠY SAU ĐẠI HỌC]
 - Dự kiến 50-100 trích dẫn             - Tối ưu hóa Resonator                  - Giáo trình Cơ học phi tuyến
 - Công bố trên các tạp chí             - Dự báo mất ổn định Pull-in            - Chuyên đề Tuyến tính hóa
   ISI Q1 (JSV, IJES, MSSP)             - Thiết kế cảm biến đo gia tốc            và Đàn hồi phi cổ điển
  • Tác động học thuật quốc tế: Khung lý thuyết của luận án đóng góp vào dòng chảy của Cơ học tính toán (Computational Mechanics), cung cấp một công cụ giải tích cạnh tranh trực tiếp với các phương pháp của trường phái Ji-Huan He và Beléndez. Các công trình công bố từ luận án trên các tạp chí quốc tế uy tín khẳng định vị thế của trường phái Cơ học phi tuyến Việt Nam (kế thừa từ GS. Nguyễn Đông Anh, GS. Đào Huy Bích).
  • Chuyển đổi công nghiệp & Công nghệ cao: Trực tiếp hỗ trợ ngành công nghiệp vi điện tử, thiết kế vi hệ thống MEMS/NEMS trong việc tối ưu hóa dải tần làm việc của cảm biến vi cơ học, nâng cao tuổi thọ và độ bền chống dính tĩnh điện.

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

  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Tiếp cận một phương pháp giải tích tường minh, dễ cài đặt mã nguồn, giải quyết nhanh các phương trình vi phân phi tuyến bậc cao mà không cần giải tích số phức tạp.
  • Các nhà khoa học Cơ học lý thuyết: Khai thác mối liên hệ giữa toán tử trung bình có trọng số và biến đổi Laplace để giải quyết các bài toán biên đàn hồi phi cổ điển phức tạp hơn.
  • Kỹ sư R&D trong ngành vi cơ điện tử: Sở hữu bộ công thức dạng đóng để tính toán nhanh kích thước hình học tối ưu cho dầm micro/nano trong cảm biến gia tốc, con quay vi cơ và bộ lọc tần số nano.

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 thiết lập cầu nối đẳng cấu giữa toán tử trung bình có trọng số trong miền thời gian và phép biến đổi tích phân Laplace trong miền tần số phức:

$$\langle x(t) \rangle_w = s^2 \mathcal{L}{\tau x(\tau)}(s)$$

Công trình này mở rộng trực tiếp Lý thuyết Tuyến tính hóa tương đương của Caughey (1963) và hoàn thiện ý tưởng toán tử trung bình có trọng số của GS. Nguyễn Đông Anh (2015), biến phương pháp EQL từ một công cụ gần đúng thô sơ cho hệ phi tuyến yếu thành một phương pháp giải tích chính xác cao cho các hệ phi tuyến rất mạnh.

2. Điểm cải tiến phương pháp luận so với ít nhất 2 nghiên cứu quốc tế tương đương?

  • So với Phương pháp Nhiễu đồng luân (HPM) của Beléndez và cộng sự (2009, 2011): HPM đòi hỏi khai triển chuỗi nhúng tham số đồng luân $p$, làm phát sinh hàng loạt phương trình vi phân tuyến tính phụ thuộc và đòi hỏi triệt tiêu số hạng trường kỳ (secular terms) rất nặng nề ở các bậc cao ($p^2, p^3$). Phương pháp EQL trọng số của luận án chỉ cần giải đại số một bước thông qua tiêu chuẩn cực tiểu sai số, cho nghiệm dạng đóng ngay lập tức mà vẫn đảm bảo sai số $< 1%$.
  • So với Phương pháp Biến phân (Variational Approach) của Şimşek (2010): Phương pháp của Şimşek yêu cầu tính toán hàm thế năng biến phân tích phân trên toàn miền không gian và phụ thuộc vào giả thiết tọa độ biến phân. Phương pháp của luận án kết hợp EQL trọng số với biến đổi Laplace, loại bỏ hoàn toàn tích phân lượng giác phức tạp của các số hạng phi tuyến nền đàn hồi $K_{NL} w^3$.

3. Phát hiện bất ngờ nhất về mặt dữ liệu định lượng là gì?

Phát hiện bất ngờ nhất nằm ở độ ổn định sai số tuyệt đối của phương pháp khi số bậc phi tuyến tăng lên cực hạn. Thông thường trong các phương pháp giải tích (như Cân bằng điều hòa xấp xỉ bậc 1), sai số sẽ bùng nổ khi bậc đa thức tăng từ bậc 3 lên bậc 13. Tuy nhiên, trong luận án, đối với hệ Duffing bậc 13, sai số tương đối giữa nghiệm giải tích EQL trọng số ($s=1$) và nghiệm số RK45 không những không tăng mà còn tiệm cận về $0.05%$, chứng minh tính chất "ngưng tụ thông tin" đặc biệt của hàm trọng số $h(t) = s^2\omega^2 t e^{-s\omega t}$ trong nửa chu kỳ đầu.

4. Giao thức tái lập (Replication Protocol) có được cung cấp đầy đủ không?

Có. Luận án cung cấp bảng tường minh giá trị trung bình có trọng số của toàn bộ các hàm điều hòa cơ bản ($\cos^n(\tau), \sin^n(\tau)$ đến bậc 6), các biểu thức đại số giải tích của độ cứng tương đương $k_{eq}$, cùng hệ phương trình vi phân chuẩn hóa không thứ nguyên cho cả hệ 1 bậc tự do và dầm micro/nano, cho phép tái lập 100% kết quả trên MATLAB hoặc Mathematica mà không cần giả định ẩn.

5. Khung chương trình nghiên cứu 10 năm được vạch ra như thế nào?

Luận án định hình lộ trình phát triển gồm:

  1. Giai đoạn 1-3 năm: Mở rộng cho hệ ngẫu nhiên phi tuyến chịu kích động không dừng (non-stationary random excitation).
  2. Giai đoạn 4-6 năm: Áp dụng cho bài toán tương tác chất lỏng - kết cấu quy mô nano (ống nano dẫn dòng chất lỏng chịu kích động nhiệt).
  3. Giai đoạn 7-10 năm: Xây dựng nền tảng phần mềm giải tích tự động hóa phục vụ công nghiệp thiết kế hệ thống vi cơ điện tử MEMS/NEMS quốc gia.

Kết luận

Luận án tiến sĩ của Đặng Văn Hiếu tạo dựng các đóng góp then chốt:

  1. Xây dựng lý thuyết hoàn chỉnh cho phương pháp tuyến tính hóa tương đương kết hợp toán tử trung bình có trọng số áp dụng cho các hệ dao động phi tuyến tiền định.
  2. Chứng minh bản chất toán học của toán tử trung bình có trọng số như một phép biến đổi Laplace hữu tỉ, đơn giản hóa triệt để quá trình tính toán giải tích.
  3. Giải quyết thành công 5 lớp bài toán dao động 1 bậc tự do phi tuyến phức tạp (Duffing bậc 3 đến bậc 13, Duffing-điều hòa, giếng đôi, số mũ hữu tỉ, gián đoạn) với độ chính xác tương đương nghiệm số chính xác cao.
  4. Thiết lập mô hình giải tích mới cho dầm Micro trên nền đàn hồi 3 tham số, làm sáng tỏ vai trò của tham số chiều dài vật liệu và tỉ số độ mảnh.
  5. Khám phá cơ chế cạnh tranh độ cứng vi mô trên dầm Nano MEMS/NEMS tĩnh điện, kết hợp đồng thời lý thuyết phi cục bộ và độ dốc biến dạng để dự báo chính xác ngưỡng sập pull-in.
  6. Mở ra hướng nghiên cứu mới trong cơ học tính toán: sử dụng các toán tử tích phân trọng số để đại số hóa và tối ưu hóa các bài toán biên phi tuyến mạnh trong khoa học kỹ thuật hiện đại.