Tổng quan về luận án

Dao động phi tuyến là một hiện tượng phổ biến và giữ vai trò cốt lõi trong nhiều hệ thống kỹ thuật hiện đại, trải dài từ các công trình vĩ mô như cầu treo, nhà cao tầng, phương tiện giao thông cho đến các cấu trúc vi mô và nano trong hệ thống vi cơ điện tử (MEMS/NEMS) hay kính hiển vi lực nguyên tử (AFM). Trong khi dao động tuyến tính chỉ là sự lý tưởng hóa trong phạm vi chuyển động vô cùng bé, hầu hết các hệ vật lý thực tế đều vận hành theo quy luật phi tuyến phức tạp được mô tả bởi các phương trình vi phân thường hoặc đạo hàm riêng phi tuyến. Ngoại trừ một số rất ít trường hợp đặc biệt có nghiệm giải tích chính xác (dưới dạng tích phân elliptic), việc phân tích động lực học phi tuyến phụ thuộc lớn vào phương pháp số và các phương pháp giải tích xấp xỉ. Tuy nhiên, các phương pháp số thuần túy không thể thiết lập được biểu thức tường minh của mối quan hệ biên độ – tần số ($\omega - A$), vốn là đặc trưng động lực học then chốt nhất của hệ phi tuyến.

Khoảng trống nghiên cứu then chốt xuất phát từ giới hạn nội tại của phương pháp tuyến tính hóa tương đương kinh điển được Caughey (1963) phát triển dựa trên kỹ thuật Krylov-Bogoliubov. Dù phương pháp tuyến tính hóa tương đương có ưu điểm nổi bật là quy trình đại số đơn giản, nó sử dụng toán tử trung bình thời gian cổ điển trên một chu kỳ: $$\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ộ thông tin của các số hạng điều hòa bậc cao do tính chất đối xứng: $$\langle \cos(nt) \rangle_c = \langle \sin(nt) \rangle_c = 0 \quad (\forall n \in \mathbb{N}^*)$$ Hệ quả là khi bậc phi tuyến tăng cao hoặc khi hệ chứa các hàm phi tuyến phức tạp (dạng phân số hữu tỉ, căn thức vô tỉ, thế năng giếng đôi hay tính không liên tục), phương pháp tuyến tính hóa cổ điển dẫn đến sai số rất lớn, thậm chí không thể chấp nhận được. Dù đã có nhiều nỗ lực cải tiến như phương pháp tuyến tính hóa điều chỉnh (Anh & Di Paola, 1995), tiêu chuẩn cực tiểu sai số thế năng (Zhang et al., 1991), tiêu chuẩn tuyến tính hóa từng phần (Elishakoff & Cai, 1993), hay tiêu chuẩn đối ngẫu (Nguyen Dong Anh, 2010), các cách tiếp cận này chủ yếu giải quyết hệ ngẫu nhiên hoặc giới hạn ở đa thức bậc thấp. Luận án tiến sĩ của Đặng Văn Hiếu (2021) với tiêu đề "Phân tích dao động phi tuyến bằng cách tiếp cận trung bình có trọng số" đã tạo ra bước đột phá lý thuyết khi kết hợp phương pháp tuyến tính hóa tương đương với toán tử trung bình có trọng số do GS. Nguyễn Đông Anh (2015) khởi xướng, thiết lập hệ phương pháp giải tích hoàn chỉnh cho các hệ dao động tiền định không cản từ một bậc tự do đến hệ liên tục phụ thuộc kích thước.

Luận án giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 4 giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để xây dựng thuật toán giải tích tổng quát xác định tần số phi tuyến cho các hệ dao động tiền định có tính phi tuyến tùy ý mà không làm mất thông tin điều hòa bậc cao?
    • H1: Hàm trọng số tổ hợp dạng tích của hàm tăng tuyến tính và hàm suy giảm hàm mũ cho phép bảo toàn các đặc trưng chu kỳ và liên hệ trực tiếp với biến đổi Laplace.
  • RQ2: Mức độ chính xác và phạm vi hội tụ của phương pháp trung bình có trọng số đối với các hệ dao động phi tuyến mạnh bậc cao (lên đến bậc 13), thế năng giếng đôi và hàm phi tuyến không khả vi/không liên tục đạt được như thế nào khi so sánh với nghiệm số Runge-Kutta bậc 4?
    • H2: Tồn tại tham số điều chỉnh tối ưu $s = n/2$ giúp triệt tiêu sai số phương trình và đưa sai số tương đối của tần số về dưới 1-3%.
  • RQ3: Hiệu ứng kích thước hình học, tỉ số độ cứng chống uốn và tỉ số độ mảnh chi phối như thế nào đến tần số dao động phi tuyến của dầm vi mô (microbeam) tựa trên nền đàn hồi phi tuyến ba tham số Pasternak-Winkler?
    • H3: Lý thuyết ứng suất cặp sửa đổi (MCST) kết hợp với trung bình có trọng số sẽ phản ánh chính xác hiệu ứng tăng cứng (stiffening effect) khi kích thước mặt cắt tiệm cận tham số độ dài vật liệu $l$.
  • RQ4: Sự cạnh tranh giữa hiệu ứng phi cục bộ (làm mềm kết cấu) và gradient biến dạng (làm cứng kết cấu) tác động như thế nào đến tần số phi tuyến và hiện tượng mất ổn định tĩnh điện (pull-in instability) trong dầm nano NEMS?
    • H4: Lý thuyết độ dốc biến dạng phi cục bộ (NSGT) tích hợp trong mô hình giải tích sẽ phân lập chính xác ranh giới ổn định động lực học dưới tác dụng đồng thời của lực nén dọc trục và điện thế tĩnh điện phi tuyến.

Khung lý thuyết của luận án tích hợp: Lý thuyết tuyến tính hóa tương đương (Equivalent Linearization), Giải tích biến đổi Laplace (Laplace Transformation), Lý thuyết ứng suất cặp sửa đổi (Modified Couple Stress Theory - Yang et al., 2002) và Lý thuyết độ dốc biến dạng phi cục bộ (Nonlocal Strain Gradient Theory - Lim et al., 2015). Phạm vi nghiên cứu bao quát các hệ 1 bậc tự do điển hình (Duffing bậc 3, 5, 7, 9, 11, 13; Duffing-điều hòa; thế năng giếng đôi; số mũ hữu tỉ; hệ phi tuyến từng đoạn không liên tục) và các hệ liên tục (dầm micro trên nền đàn hồi phi tuyến $K_L, K_P, K_{NL}$, dầm nano MEMS/NEMS chịu kích thích tĩnh điện $V_0$ và lực nén $P_0$).

Literature Review và Positioning

Nghiên cứu dao động phi tuyến trong nửa thế kỷ qua tập trung vào bốn dòng phương pháp giải tích gần đúng chính:

  1. Dòng phương pháp nhiễu (Perturbation Methods): Khởi xướng từ các khai triển tiệm cận cổ điển của Nayfeh (1973, 1995) dành cho hệ phi tuyến yếu. Để vượt qua rào cản tham số bé, J.H. He đã phát triển phương pháp nhiễu tham số (1999), nhiễu đồng luân (HPM - 2000) và nhiễu lặp (2001). Beléndez et al. (2009, 2011) đã cải tiến HPM để tìm nghiệm bậc cao cho hệ có số mũ hữu tỉ. Tuy nhiên, phương pháp nhiễu bộc lộ hạn chế lớn khi giải tích các hệ có lực phục hồi vô tỉ hoặc thế năng đa giếng do cấu trúc đại số tiệm cận bùng nổ phức tạp.
  2. Dòng phương pháp cân bằng điều hòa (Harmonic Balance Method - HBM): Phát triển mạnh mẽ bởi Mickens (1984, 1996, 2010) cho các hệ có số mũ hữu tỉ và phân số; Hu & Tang (2006) cho hệ Duffing-điều hòa; Beléndez et al. (2007) cho hệ có lực phục hồi vô tỉ; Hosen et al. (2012) và Razzak (2016) cho các nghiệm xấp xỉ bậc cao. Dù HBM xấp xỉ bậc nhất đơn giản nhưng cho sai số lớn, còn HBM bậc cao đòi hỏi giải hệ phương trình phi tuyến siêu việt nhiều ẩn cực kỳ phức tạp.
  3. Dòng phương pháp năng lượng và biến phân (Energy & Variational Approaches): Khởi xướng bởi J.H. He với phương pháp cân bằng năng lượng (EBM - 2002), phương pháp biến phân tổng quát (VA - 2007) và phương pháp Hamiltonian (2010). Các phương pháp này được Younesian et al. (2010) áp dụng cho dao động Duffing tổng quát, Shou (2009) cho hệ phi tuyến không liên tục, Sedighi et al. (2012) và Yildirim et al. (2012) cho các phương trình dao động bậc cao. Nhược điểm chính của EBM là phụ thuộc vào việc lựa chọn điểm định vị (collocation point), trong khi VA bậc cao đòi hỏi các phép tích phân năng lượng phức tạp.
  4. Dòng phương pháp tuyến tính hóa tương đương (Equivalent Linearization - EQL): Đặt nền móng bởi Caughey (1963), Iwan (1973), Atalik & Utku (1976). Tại Việt Nam và trên thế giới, GS. Nguyễn Đông Anh cùng các cộng sự đã tạo ra nhiều đột phá: tuyến tính hóa điều chỉnh (Anh & Di Paola, 1995), tiêu chuẩn đối ngẫu (Anh, 2010), tiêu chuẩn đối ngẫu có trọng số (Anh & Linh, 2013). Tuy nhiên, hầu hết các công trình này đều hướng vào dao động ngẫu nhiên (random vibrations).

Về cơ học vi mô và nano, các nghiên cứu của Şimşek (2010, 2015), Reddy (2011), Akgöz & Civalek (2013), Fu et al. (2008), Qian et al. (2009) đã áp dụng lý thuyết ứng suất cặp sửa đổi (MCST) hoặc lý thuyết độ dốc biến dạng phi cục bộ (NSGT - Lim et al., 2015). Tuy nhiên, các nghiên cứu trước đây chưa khảo sát tường minh ảnh hưởng đồng thời của tỉ số độ cứng chống uốn và tỉ số độ mảnh trên nền 3 thông số, cũng như thiếu một công cụ giải tích dạng đóng cho bài toán NEMS chịu trường lực tĩnh điện phi tuyến mạnh. Luận án định vị chính xác vào khoảng trống này, cung cấp lời giải giải tích dạng đóng có độ chính xác cao tương đương nghiệm số.

Đó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 (1963) bằng việc thay thế tiên đề trung bình chu kỳ cổ điển bằng toán tử trung bình có trọng số của Nguyễn Đông Anh (2015). Trích dẫn trực tiếp từ luận án:

"Phương pháp tuyến tính hóa tương đương với trung bình cổ điển thường có những nhược điểm đó là khi tính phi tuyến của bài toán tăng sẽ dẫn đến kết quả của phương pháp này thường không chính xác, đôi khi không thể chấp nhận được."

Về mặt bản chất nhận thức luận (epistemological stance), luận án thực hiện bước chuyển dịch mô hình (paradigm shift) từ việc xem xét giá trị trung bình phẳng đồng đều trên miền thời gian tuần hoàn $[0, T]$ sang mô hình trung bình tập trung có trọng số suy giảm thời gian (memory-weighted decay). Luận án chứng minh định lý toán học:

"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 qua phép biến đổi Laplace."

Phương trình toán tử cốt lõi được xác lập: $$\langle x(t) \rangle_w = \int_0^\infty s^2 \omega^2 t e^{-s\omega t} x(t) dt = s^2 \left. \left[ -\frac{d}{dp} \mathcal{L}{x(t)} \right] \right|_{p = s\omega}$$

Nhờ tính chất này, toàn bộ các tích phân trung bình phức tạp của hàm điều hòa được chuyển hóa thành các hàm đại số hữu tỉ theo tham số điều chỉnh $s$: $$\langle \cos(nt) \rangle_w = s^2 \frac{s^2 - n^2}{(s^2 + n^2)^2}$$ $$\langle \cos^2(t) \rangle_w = \frac{s^4 + 2s^2 + 8}{(s^2 + 4)^2}$$ $$\langle \cos^4(t) \rangle_w = \frac{248s^4 + 416s^2 + 1536 - 28s^6 + s^8}{(s^2 + 4)^2 (s^2 + 16)^2}$$

Trích dẫn từ phân tích toán học của tác giả:

"Hàm hệ số trọng số $h(t)$ điều chỉnh giá trị của các hàm $\cos(t)$ và $\cos^2(t)$, giữ lại chính các đặc tính tuần hoàn của các hàm $\cos(t)$ và $\cos^2(t)$, cũng như làm ngưng tụ giá trị của các hàm này trong chu kỳ đầu tiên, giữ lại phần trọng trong nửa chu kỳ đầu tiên, làm giảm các giá trị khác biệt lớn nhất và nhỏ nhất cũng như điều chỉnh các hàm theo thời gian."

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

Khung phân tích của luận án kết hợp đa tầng 4 trụ cột lý thuyết:

  1. Lý thuyết tuyến tính hóa tương đương tiền định: Thiết lập phương trình sai số $e(X, \dot{X}) = g(X, \dot{X}) - 2\alpha\dot{X} - \omega^2 X$ và cực tiểu hóa $\langle e^2 \rangle_w \to \text{Min}$.
  2. Giải tích biến đổi Laplace cho hàm trọng số: Thiết lập biểu thức giải tích đóng cho các hàm phi tuyến đa thức bậc lẻ từ bậc 3 đến bậc 13, hàm phi tuyến phân thức, vô tỉ và không liên tục.
  3. Lý thuyết ứng suất cặp sửa đổi (MCST): Mô tả biến dạng dầm vi mô với mật độ năng lượng biến dạng $U_s = \frac{1}{2} \int_V (\sigma_{ij}\varepsilon_{ij} + m_{ij}\chi_{ij}) dV$, đưa vào tham số độ dài vật liệu $l$ nhằm nắm bắt hiệu ứng kích thước micro.
  4. Lý thuyết độ dốc biến dạng phi cục bộ (NSGT): Tổ hợp phương trình cấu trúc $\left(1 - e_0^2 a^2 \nabla^2\right) \sigma_{xx} = E\left(1 - l^2 \nabla^2\right) \varepsilon_{xx}$, giải quyết đồng thời tương tác tầm xa nguyên tử ($e_0a$) và gradient biến dạng mạng tinh thể ($l$).
Hệ thống khảo sát Phương trình phi tuyến gốc Toán tử tuyến tính hóa có trọng số
Duffing bậc cao ($2n+1$) $\ddot{x} + \omega_0^2 x + \sum_{k=1}^n \alpha_{2k+1} x^{2k+1} = 0$ $\omega_{eq}^2 = \omega_0^2 + \sum_{k=1}^n \alpha_{2k+1} A^{2k} \langle \cos^{2k+2}(t) \rangle_w / \langle \cos^2(t) \rangle_w$
Duffing-Điều hòa $\ddot{x} + \frac{\alpha x^m}{1 + \beta x^n} = 0$ Biến đổi Taylor kết hợp trọng số qua chuỗi Fourier rút gọn
Dầm Micro trên nền đàn hồi $EI_{eq} \frac{\partial^4 w}{\partial x^4} + \rho A \frac{\partial^2 w}{\partial t^2} + K_L w - K_P \frac{\partial^2 w}{\partial x^2} + K_{NL} w^3 = 0$ Tuyến tính hóa tương đương không gian - thời gian qua hàm dạng Galerkin $\phi(x)$
Dầm Nano NEMS tĩnh điện $\mathcal{L}_{NSGT}(w) + \rho A \ddot{w} = \frac{\epsilon_0 b V_0^2}{2(g_0 - w)^2} - \frac{\epsilon_0 b V_0^2}{2(g_0 + w)^2}$ Khai triển lực tĩnh điện phi tuyến + Tuyến tính hóa trọng số

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ủ chặt chẽ thế giới quan thực chứng (positivism) và cơ học môi trường liên tục giải tích (analytical continuum mechanics). Thiết kế nghiên cứu là sự kết hợp song hành giữa suy diễn giải tích toán học (analytical deduction) và kiểm chứng mô phỏng số học (numerical validation).

Quy trình nghiên cứu rigorous

Quy trình giải tích được chuẩn hóa theo 5 bước nghiêm ngặt:

  1. Thiết lập phương trình sai số: Phương trình vi phân phi tuyến tổng quát $\ddot{X} + g(X, \dot{X}) = 0$ được thay thế bởi hệ tuyến tính $\ddot{X} + \omega_{eq}^2 X = 0$. Sai số phương trình xác định bởi $e(X) = g(X, 0) - \omega_{eq}^2 X$.
  2. Cực tiểu hóa sai số bình phương trung bình có trọng số: Điều kiện $\frac{\partial \langle e^2 \rangle_w}{\partial (\omega_{eq}^2)} = 0$ dẫn đến: $$\omega_{eq}^2 = \frac{\langle X g(X) \rangle_w}{\langle X^2 \rangle_w}$$
  3. Giải tích hóa tích phân qua biến đổi Laplace: Với nghiệm xấp xỉ dạng $X(t) = A \cos(\omega t)$, các đại lượng $\langle X g(X) \rangle_w$ và $\langle X^2 \rangle_w$ được giải tích hóa hoàn toàn thông qua bảng tích phân biến đổi Laplace cho các hàm $\cos^k(\omega t)$.
  4. Hiệu chỉnh tham số $s$: Tác giả xác lập liên hệ vật lý giữa đỉnh cực đại của hàm trọng số $t_{max} = \frac{1}{s\omega}$ với chu kỳ dao động $T = \frac{2\pi}{\omega}$. Bằng cách chọn $t_{max} = \frac{T}{n}$, thu được $s = \frac{n}{2\pi}$. Khi $s \to 0$, nghiệm tiệm cận trung bình cổ điển; khi $s = 1$ hoặc $s = 2$, hàm trọng số tập trung vào nửa chu kỳ đầu, triệt tiêu sự méo dạng phi tuyến.
  5. Kiểm chứng độc lập: Toàn bộ quan hệ tần số - biên độ thu được đều được so chuẩn đối chứng với:
    • Nghiệm tích phân elip chính xác (Exact Elliptic Solutions) đối với dao động Duffing bậc 3, 5, 7 và thế năng giếng đôi.
    • Nghiệm số trực tiếp Runge-Kutta bậc 4 giải trên phần mềm tính toán kỹ thuật số chuyên dụng.
    • Các công bố quốc tế sử dụng phương pháp cân bằng năng lượng (EBM), phương pháp biến phân (VA), phương pháp nhiễu đồng luân (HPM).

Data và phân tích

Luận án thực hiện khảo sát tham số trên tập dữ liệu mô phỏng toàn diện:

  • Hệ Duffing bậc cao: Khảo sát các hệ số phi tuyến từ yếu ($\alpha_i = 1$) đến cực mạnh ($\alpha_i = 100$) và bậc phi tuyến từ bậc 3, 5, 7, 9, 11 đến 13 với dải biên độ ban đầu $A \in [0.1, 10]$.
  • Dao động với thế năng giếng đôi: Khảo sát hai miền biên độ đặc thù $1 \le A \le \sqrt{2}$ (chuyển động quanh vị trí cân bằng lệch) và $A \ge \sqrt{2}$ (chuyển động vượt qua rào thế năng).
  • Dầm vi mô: Khảo sát tham số tỉ lệ độ dài vật liệu $\theta = h/l \in [1, 10]$, tỉ số độ cứng chống uốn $\alpha = E_c/E_m$, tỉ số độ mảnh $S = L/h \in [20, 60]$, hệ số nền Winkler $K_L \in [0, 100]$, Pasternak $K_P \in [0, 100]$, và nền phi tuyến $K_{NL} \in [0, 100]$.
  • Dầm nano NEMS: Khảo sát tham số phi cục bộ $\mu = (e_0a/L)^2 \in [0, 0.1]$, tham số gradient biến dạng $\eta = l/L \in [0, 0.1]$, điện thế tĩnh điện không thứ nguyên $V \in [0, V_{cr}]$ và lực nén dọc trục $P \in [0, P_{cr}]$.

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

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

  1. Độ chính xác vượt trội đối với hệ phi tuyến bậc cao: Đối với dao động Duffing bậc 7, bậc 9, bậc 11 và bậc 13, phương pháp trung bình có trọng số cho tần số giải tích xấp xỉ gần như trùng khít với nghiệm số Runge-Kutta bậc 4 (sai số tương đối $< 0.5%$). Trong khi phương pháp tuyến tính hóa cổ điển cho sai số tăng vọt theo biên độ $A$, phương pháp đề xuất duy trì sai số ổn định ngay cả khi hệ số phi tuyến $\alpha_{2k+1} = 100$.
  2. Giải quyết triệt để bài toán thế năng giếng đôi và bất liên tục: Luận án chứng minh phương pháp trung bình có trọng số mô tả chính xác dao động trong thế năng giếng đôi $\ddot{x} - x + x^3 = 0$ ở cả hai trường hợp $1 \le A \le \sqrt{2}$ và $A \ge \sqrt{2}$. Đối với hệ có lực phục hồi không liên tục (piecewise linear/discontinuous): $$\ddot{x} + x + \epsilon \cdot \text{sgn}(x) = 0$$ nghiệm giải tích thu được bám sát quỹ đạo pha chính xác mà không gặp hiện tượng phân kỳ toán học như các phương pháp nhiễu cổ điển.
  3. Hiệu ứng tăng cứng mạnh mẽ của dầm vi mô theo MCST: Tần số dao động phi tuyến $\omega_{NL}$ tăng đột biến khi chiều dày dầm $h$ giảm dần về kích thước tham số độ dài vật liệu $l$ (tức $\theta = h/l \to 1$). Tỉ số độ cứng chống uốn và tham số nền cắt Pasternak $K_P$ đóng vai trò chi phối lớn hơn nền Winkler $K_L$ trong việc triệt tiêu độ võng phi tuyến.
  4. Quy luật cạnh tranh kép phi cục bộ - gradient biến dạng trong NEMS: Trong dầm nano tĩnh điện, tham số phi cục bộ $e_0a$ làm mềm dầm (giảm tần số tự nhiên và giảm điện thế mất ổn định $V_{cr}$), trong khi tham số độ dài vật liệu $l$ làm cứng dầm (tăng tần số và tăng $V_{cr}$). Khi hai tham số này cân bằng ($e_0a = l$), ứng xử của dầm nano quay trở về mô hình cổ điển nhưng ở trạng thái ứng suất nội vi mô triệt tiêu.

Implications đa chiều

  • Về mặt lý thuyết cơ học: Luận án hoàn thiện hóa trường phái giải tích tuyến tính hóa tương đương của Việt Nam do GS. Nguyễn Đông Anh khởi xướng, nâng tầm ứng dụng từ miền ngẫu nhiên sang hệ tiền định phi tuyến tùy ý và hệ liên tục vi mô/nano.
  • Về mặt phương pháp luận: Cung cấp một công cụ giải tích dạng đóng thay thế cho các thuật toán HBM và VA phức tạp, giúp các nhà nghiên cứu tính toán nhanh tần số mà không cần giải hệ đại số phi tuyến bậc cao.
  • Về mặt kỹ thuật và ứng dụng công nghiệp:
    • Thiết kế cảm biến MEMS/NEMS: Cho phép dự báo chính xác ngưỡng điện thế kích hoạt tới hạn ($V_{cr}$) để tránh hiện tượng dính bám điện cực (pull-in failure).
    • Chế tạo cảm biến cộng hưởng vi mô: Tối ưu hóa tần số làm việc dựa trên việc điều chỉnh kích thước hình học $h$ so với tham số vật liệu $l$.
    • Thiết kế móng công trình và nền đường cao tốc: Mô hình dầm trên nền 3 tham số $K_L, K_P, K_{NL}$ cung cấp công thức kiểm tra dao động đàn - dẻo phi tuyến dưới tải trọng động.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 3 hạn chế chính:

  1. Giới hạn ở hệ dao động tự do không cản (undamped free vibration): Chưa mở rộng sang các hệ có cản phi tuyến phức tạp (cản Rayleigh, cản tỷ lệ với lũy thừa vận tốc $|\dot{x}|\dot{x}$) hoặc dao động cưỡng bức tuần hoàn có hiện tượng cộng hưởng phụ và cộng hưởng hỗn loạn.
  2. Xấp xỉ không gian đơn mode (single-mode discretization): Sử dụng phương pháp Galerkin 1 mode để đưa dầm micro/nano về hệ 1 bậc tự do tương đương. Dù chính xác với dao động cơ bản, mô hình này chưa phản ánh được hiện tượng truyền năng lượng giữa các mode (internal resonance).
  3. Mô hình vật liệu đồng nhất đẳng hướng: Chưa xét đến ảnh hưởng của trường nhiệt độ ngẫu nhiên, độ ẩm môi trường hoặc vật liệu áp điện thông minh (piezoelectric/flexoelectric coupling).

Chương trình nghiên cứu tương lai định hướng:

  • Hướng 1: Mở rộng phương pháp trung bình có trọng số cho hệ dao động cưỡng bức phi tuyến có cản, xác định đường cong cộng hưởng biên độ - tần số và hiện tượng bước nhảy (jump phenomenon).
  • Hướng 2: Phát triển mô hình liên tục đa mode (multi-mode formulation) cho tấm vi mô/nano (micro/nanoplates) và vỏ nano (nanoshells).
  • Hướng 3: Ứng dụng giải tích dao động trong cảm biến sinh học nano (nanobiosensors) phát hiện khối lượng phân tử virus/kháng thể dựa trên độ dịch chuyển tần số phi tuyến.

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

  • Tác động học thuật: Các kết quả của luận án đã được công bố trên các tạp chí khoa học chuyên ngành uy tín trong nước và quốc tế. Phương pháp trung bình có trọng số mở ra một nhánh tiếp cận giải tích mới, tạo tiềm năng trích dẫn cao trong cộng đồng cơ học tính toán phi tuyến (Nonlinear Dynamics, International Journal of Non-Linear Mechanics, Composite Structures).
  • Chuyển giao công nghệ vi cơ điện tử: Đóng góp trực tiếp vào cẩm nang thiết kế linh kiện MEMS/NEMS cho các trung tâm R&D công nghệ cao, giảm thiểu chi phí thử nghiệm thực nghiệm đắt đỏ nhờ các công thức dự báo giải tích chính xác.
  • Phát triển phần mềm tính toán: Thuật toán giải tích hóa qua biến đổi Laplace dễ dàng lập trình tích hợp vào các gói phần mềm kỹ thuật như MATLAB, Mathematica hay Maple để phân tích dao động tự động.

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

  • Nghiên cứu sinh & Học viên cao học ngành Cơ học, Vật lý kỹ thuật: Tiếp cận phương pháp toán học giải tích mới, khắc phục các giới hạn cổ điển, mở rộng đề tài cho các kết cấu phức tạp.
  • Các nhà khoa học đầu ngành: Sở hữu một công cụ mạnh mẽ để đối chuẩn (benchmark) với các thuật toán số (FEM, DQM) và các phương pháp giải tích năng lượng khác.
  • Kỹ sư R&D trong ngành Vi điện tử & Bán dẫn: Sử dụng các công thức giải tích đóng để tính toán kích thước tới hạn, điện áp vận hành an toàn và tần số cộng hưởng của vi dầm trong cảm biến gia tốc, con quay vi cơ và công tắc vô tuyến RF-MEMS.
  • Kỹ sư Xây dựng - Cơ khí kết cấu: Áp dụng mô hình dầm trên nền đàn hồi 3 tham số để tính toán kết cấu ngầm, móng cọc chịu tải trọng động lớ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 thiết lập cơ sở toán học giải tích liên kết giữa phương pháp tuyến tính hóa tương đương tiền định và phép biến đổi Laplace thông qua hàm trọng số tổ hợp $h(t) = s^2 \omega^2 t e^{-s\omega t}$. Luận án đã mở rộng căn bản Lý thuyết tuyến tính hóa tương đương kinh điển của Caughey (1963) và Krylov-Bogoliubov, khắc phục hoàn toàn nhược điểm triệt tiêu số hạng điều hòa bậc cao ($\langle \cos(nt) \rangle_c = 0$) vốn tồn tại hơn 50 năm qua.

2. Phương pháp của luận án có điểm gì ưu việt hơn so với các nghiên cứu quốc tế đã công bố?

So với phương pháp Cân bằng điều hòa biến đổi (Modified HBM) của Beléndez et al. và Mickens đòi hỏi giải hệ phương trình đại số phi tuyến siêu việt phức tạp, phương pháp của luận án đưa ra công thức dạng đóng tường minh thông qua các biểu thức đại số hữu tỉ đơn giản. So với phương pháp Biến phân (VA) và Cân bằng năng lượng (EBM) của J.H. He phụ thuộc vào việc chọn điểm định vị và tích phân phức vi mô, phương pháp trung bình có trọng số tự động tối ưu hóa qua tham số $s = n/2\pi$, đạt độ chính xác tương đương hoặc cao hơn với khối lượng tính toán giảm hơn 70%.

3. Phát hiện nào trong dữ liệu gây bất ngờ nhất và có bằng chứng định lượng cụ thể ra sao?

Phát hiện bất ngờ nhất là tính ổn định tuyệt đối của sai số khi bậc phi tuyến tăng từ 3 lên đến 13. Trong khi các phương pháp giải tích cổ điển bộc lộ sai số bùng nổ khi số mũ phi tuyến tăng cao, phương pháp trung bình có trọng số duy trì sai số tương đối so với nghiệm số Runge-Kutta ở mức dưới 1% cho các hệ Duffing bậc 7, bậc 9, bậc 11 và bậc 13 với hệ số phi tuyến cực đại $\alpha_i = 100$.

4. Luận án có cung cấp quy trình tái lập (replication protocol) đầy đủ không?

Có. Luận án trình bày chi tiết bảng tra tích phân trung bình có trọng số cho toàn bộ các hàm lượng giác bậc chẵn và bậc lẻ $\cos^k(t), \sin^k(t)$ (Bảng 2.1 trong văn bản gốc), các bước biến đổi trung gian của phương trình sai số vi phân, các biểu thức thế năng giếng đôi và quy trình giải thuật số RK4 đối chứng, cho phép bất kỳ nhà nghiên cứu nào cũng có thể lập trình tái lập chính xác 100% kết quả.

5. Lộ trình phát triển học thuật 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm tập trung vào: (1) Mở rộng lý thuyết tuyến tính hóa trọng số cho dao động tham số và dao động ngẫu nhiên phi tuyến phi Gauss; (2) Tích hợp hiệu ứng nhiệt - điện - từ trong dầm và tấm nano 2D (Graphene, MXene); (3) Thiết kế vật liệu cơ tính biến thiên xốp (porous FGMs) và cấu trúc metamaterials giảm chấn động lực học phi tuyến.

Kết luận

  1. Luận án đã phát triển thành công một hệ phương pháp giải tích hoàn chỉnh, kết hợp phương pháp tuyến tính hóa tương đương với toán tử trung bình có trọng số cho các hệ dao động phi tuyến tiền định.
  2. Thiết lập bản chất toán học chặt chẽ chứng minh trung bình có trọng số là ảnh Laplace của hàm thời gian nhân tử, bảo toàn trọn vẹn thông tin điều hòa bậc cao.
  3. Giải quyết thành công hàng loạt bài toán dao động phi tuyến phức tạp 1 bậc tự do: Duffing bậc 3 đến bậc 13, Duffing-điều hòa, thế năng giếng đôi, số mũ hữu tỉ và hệ có lực phục hồi bất liên tục với sai số nhỏ hơn 1-3% so với nghiệm chính xác.
  4. Xây dựng mô hình giải tích mới phân tích dao động phi tuyến của dầm vi mô trên nền đàn hồi 3 tham số Pasternak-Winkler theo lý thuyết ứng suất cặp sửa đổi, làm sáng tỏ vai trò của tỉ số độ cứng chống uốn và độ mảnh.
  5. Giải quyết trọn vẹn bài toán dao động phi tuyến của dầm nano NEMS dưới tác dụng của lực tĩnh điện phi tuyến theo lý thuyết độ dốc biến dạng phi cục bộ, phân lập cơ chế cạnh tranh đa tỉ lệ giữa hiệu ứng làm mềm phi cục bộ và hiệu ứng làm cứng gradient biến dạng.
  6. Mở ra 3 nhánh nghiên cứu đột phá mới về giải tích dao động phi tuyến trong công nghệ vi cơ điện tử, vật liệu nano tiên tiến và cơ học công trình hiện đại.