Tổng quan luận án
Hiện tượng dao động xuất hiện phổ biến trong tự nhiên và nhiều lĩnh vực kỹ thuật như cơ khí, xây dựng, giao thông, mạch điện tử và các hệ vi cơ điện tử (MEMS/NEMS). Trong thực tế, hầu hết các hệ dao động kỹ thuật đều có đặc tính phi tuyến, được mô tả toán học bằng các phương trình vi phân thường hoặc đạo hàm riêng phi tuyến. Khác với dao động tuyến tính, các bài toán dao động phi tuyến phần lớn không tồn tại nghiệm chính xác dạng giải tích, ngoại trừ một số trường hợp rất đặc thù. Các phương pháp số tuy có thể mô phỏng đáp ứng theo thời gian nhưng không xác định được mối quan hệ giải tích tường minh giữa biên độ và tần số – đặc trưng quan trọng nhất của dao động phi tuyến.
Để giải quyết bài toán này, các phương pháp giải tích gần đúng đóng vai trò quan trọng. Trong đó, phương pháp tuyến tính hóa tương đương (Equivalent Linearization method) được sử dụng rộng rãi nhờ tính đơn giản trong thuật toán. Tuy nhiên, phương pháp tuyến tính hóa tương đương truyền thống sử dụng giá trị trung bình cổ điển (lấy trung bình theo thời gian trên một chu kỳ) bộc lộ nhược điểm lớn: khi tính phi tuyến của hệ tăng cao, sai số của nghiệm xấp xỉ tăng nhanh và độ chính xác giảm sút. Nhằm khắc phục hạn chế này, GS. Nguyễn Đông Anh (2015) đã đề xuất phương pháp lấy trung bình mới của hàm tiền định, được gọi là trung bình có trọng số (weighted averaging). Dựa trên cơ sở đó, luận án thực hiện đề tài "Phân tích dao động phi tuyến bằng cách tiếp cận trung bình có trọng số" nhằm thiết lập một khung phân tích giải tích hoàn chỉnh cho các hệ dao động phi tuyến tiền định.
Mục tiêu tổng quát của luận án là phát triển kỹ thuật kết hợp giữa phương pháp tuyến tính hóa tương đương và khái niệm trung bình có trọng số để phân tích đáp ứng của các hệ dao động phi tuyến tự do không cản. Mục tiêu cụ thể được luận án xác định theo hai nội dung chính:
- Xây dựng cơ sở lý thuyết hoàn chỉnh cho bài toán phân tích đáp ứng của các hệ dao động phi tuyến tiền định không cản dựa trên phương pháp tuyến tính hóa tương đương kết hợp với trung bình có trọng số; làm rõ bản chất toán học, tính tổng quát của thuật toán trung bình có trọng số, ưu điểm so với trung bình cổ điển và mối liên hệ với phép biến đổi Laplace.
- Ứng dụng phương pháp đề xuất để khảo sát đáp ứng của các hệ dao động phi tuyến tiền định cụ thể, bao gồm: dao động Duffing từ bậc thấp đến bậc cao, dao động Duffing-điều hòa, dao động phi tuyến mở rộng, dao động phi tuyến với số mũ hữu tỉ, dao động phi tuyến có tính không liên tục, dao động phi tuyến của hệ liên tục gồm dầm micro và dầm nano.
Đối tượng và phạm vi nghiên cứu của luận án được xác định rõ ràng:
- Đối tượng nghiên cứu: Các bài toán dao động phi tuyến tự do không cản của hệ một bậc tự do (dao động Duffing bậc 3, 5, 7, 9, 11, 13; dao động Duffing thế năng dạng giếng đôi; dao động Duffing-điều hòa; dao động phi tuyến mở rộng; dao động có số mũ hữu tỉ; dao động không liên tục) và hệ liên tục (dầm micro tựa trên nền đàn hồi và dầm nano chịu lực tĩnh điện trong MEMS/NEMS).
- Phạm vi nghiên cứu: Giới hạn trong các hệ tiền định, dao động tự do không xét lực cản; mô hình dầm micro được phân tích theo lý thuyết ứng suất cặp sửa đổi (modified couple stress theory); mô hình dầm nano trong MEMS/NEMS được phân tích theo lý thuyết độ dốc biến dạng phi cục bộ (nonlocal strain gradient theory).
Tổng quan tài liệu và vị trí của luận án
Luận án đã hệ thống hóa các hướng tiếp cận giải tích gần đúng kinh điển và hiện đại trong phân tích dao động phi tuyến:
- Phương pháp nhiễu (Perturbation methods): Được tổng kết qua các giáo trình kinh điển [4–9]. Các mở rộng bao gồm phương pháp nhiễu tham số (Parameterized Perturbation) [10], nhiễu đồng luân (Homotopy Perturbation) [11], nhiễu lặp (Iteration Perturbation) [12] do J.H. He đề xuất, và phương pháp nhiễu đồng luân sửa đổi của Beléndez và cộng sự [13] cho hệ có số mũ hữu tỉ. Hạn chế của phương pháp nhiễu là quy trình tìm nghiệm phức tạp khi hệ chứa số hạng phi tuyến hữu tỉ hoặc vô tỉ.
- Phương pháp cân bằng điều hòa (Harmonic Balance method): Được Mickens và cộng sự [14–19] ứng dụng cho hệ một bậc tự do và hệ có số mũ hữu tỉ; Hu và Tang [20] áp dụng cho dao động Duffing-điều hòa; Razzak [21] tìm nghiệm cho hệ có lực phục hồi dạng hữu tỉ/vô tỉ; Beléndez và cộng sự [22, 23] khảo sát dao động bảo toàn; Stupnicka [24] nghiên cứu cộng hưởng tham số nhiều bậc tự do; Hosen và cộng sự [25] khảo sát hệ phi tuyến mạnh; Wu và cộng sự [26] kết hợp phương pháp Newton và cân bằng điều hòa cho hệ Duffing giếng đôi. Nhược điểm là nghiệm bậc một có độ chính xác thấp, trong khi các xấp xỉ bậc cao đòi hỏi giải hệ phương trình đại số phi tuyến phức tạp.
- Phương pháp khai triển tham số (Parameter Expansion method): Gồm phương pháp tham số bookkeeping [27] và Lindstedt-Poincare sửa đổi [28–30] của J.H. He; Ozis và Yildirm [31] áp dụng cho dao động có số mũ hữu tỉ. Phương pháp này chỉ phù hợp với hàm phi tuyến dạng đa thức.
- Các phương pháp năng lượng: Phương pháp cân bằng năng lượng (Energy Balance Method - EBM) do J.H. He đề xuất [32], phương pháp biến phân tổng quát (Variational Approach) [33] và phương pháp Hamiltonian [34]. Phương pháp này được Younesian và cộng sự [35] áp dụng cho dao động Duffing tổng quát; Momeni và cộng sự [36] cho dao động Duffing-điều hòa; Saadatnia và cộng sự [37], Shou [38], Yildirim và cộng sự [39, 40] cho các hệ phi tuyến gián đoạn, hữu tỉ và vô tỉ.
- Nghiên cứu dao động của dầm micro/nano và MEMS/NEMS: Emam và Nayfeh [44], Fallah và Aghdam [45, 46], Şimşek [47, 65, 66, 74], Gunda và cộng sự [48], Azrar và cộng sự [49], Sedighi và Reza [50]. Về lý thuyết đàn hồi phụ thuộc kích thước, luận án điểm lại lý thuyết ứng suất cặp [54–57], lý thuyết độ dốc biến dạng thứ nhất [58], lý thuyết đàn hồi phi cục bộ [59, 60], lý thuyết ứng suất cặp sửa đổi của Yang và cộng sự [61], lý thuyết độ dốc biến dạng của Mindlin [62], Aifantis [63], và lý thuyết độ dốc biến dạng phi cục bộ của Lim và cộng sự [64]. Các ứng dụng trong MEMS/NEMS có các công trình của Fu và cộng sự [78], Qian và cộng sự [79], Sadeghzadeh và Kabiri [80], Luo và Wang [81], Younis [82], Chaterjee và Pohit [83], Krylov [84], Ma và cộng sự [85], Fu và Zhang [86], Duan và Rach [87], Abdi và cộng sự [88], Miandoab và cộng sự [89], Sedighi [90], Esfahani và cộng sự [91].
- Phương pháp tuyến tính hóa tương đương và các cải tiến: Caughey [1] đề xuất cho hệ ngẫu nhiên từ ý tưởng của Krylov-Bogoliubov [6–9]; Iwan và cộng sự [92, 93], Atalik và Utku [94], Park và cộng sự [95], Su và cộng sự [96], Younespour và cộng sự [97]. Các cải tiến bao gồm: tuyến tính hóa tương đương điều chỉnh (Anh và Di Paola [98], Elishakoff và cộng sự [99]); tiêu chuẩn cực tiểu sai số thế năng (Zhang và cộng sự [100]); tiêu chuẩn tuyến tính hóa từng phần (Elishakoff và Cai [101]); tuyến tính hóa tương đương có đuôi (Fujimura và Kiureghian [102–104]); tiêu chuẩn đối ngẫu (Nguyễn Đông Anh [105, 106]); tiêu chuẩn đối ngẫu có trọng số (Nguyễn Đông Anh và Nguyễn Ngọc Linh [107–109]).
- Tình hình nghiên cứu trong nước: Nhóm GS. Nguyễn Đông Anh [98, 105–109, 116–118]; nhóm GS. Nguyễn Văn Khang [119–121] nghiên cứu dao động với đạo hàm cấp phân số; nhóm GS. Đào Huy Bích và GS. Nguyễn Đăng Bích [122–125] nghiên cứu phương pháp lặp hệ số thay đổi và xấp xỉ nối tiếp; Nguyễn Thái Chung và Phạm Lê Bình [126] nghiên cứu dầm nứt bằng phần tử hữu hạn; các luận án tiến sĩ gần đây của Nguyễn Ngọc Linh [109], Lưu Xuân Hùng [127], Dương Ngọc Hảo [128], Nguyễn Như Hiếu [129].
Từ tổng quan tài liệu, luận án xác định ba khoảng trống nghiên cứu then chốt:
- Phương pháp tuyến tính hóa tương đương chủ yếu phát triển cho các hệ dao động phi tuyến ngẫu nhiên; chưa có công trình nghiên cứu hệ thống áp dụng cho hệ dao động tiền định dựa trên trung bình có trọng số.
- Trong bài toán dao động của dầm micro trên nền đàn hồi ba lớp theo lý thuyết ứng suất cặp sửa đổi của Şimşek [65], ảnh hưởng của tỉ số độ cứng chống uốn (flexural rigidity ratio) và tỉ số độ mảnh (slenderness ratio) vẫn chưa được làm rõ.
- Trong mô hình MEMS/NEMS dầm micro/nano đặt giữa hai bản tích điện [78–80], các nghiên cứu trước đây dùng lý thuyết đàn hồi cổ điển nên chưa quan sát được hiệu ứng phụ thuộc kích thước. Bài toán dao động của mô hình này dựa trên lý thuyết độ dốc biến dạng phi cục bộ vẫn chưa được thực hiện.
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án kết hợp phương pháp giải tích và phương pháp số để phân tích và kiểm chứng đáp ứng của các bài toán dao động phi tuyến.
1. Khung lý thuyết tuyến tính hóa tương đương cho hệ tiền định
Xét phương trình vi phân phi tuyến một bậc tự do dạng tiền định:
$$\ddot{X} + g(X, \dot{X}) = F(t)$$
Phương trình tuyến tính hóa tương đương được thiết lập dưới dạng:
$$\ddot{X} + \alpha \dot{X} + \beta X = F(t)$$
Theo tiêu chuẩn độ lệch trung bình bình phương (tiêu chuẩn kinh điển của Caughey [1]):
$$\langle e^2(X, \dot{X}) \rangle = \langle [g(X, \dot{X}) - \alpha \dot{X} - \beta X]^2 \rangle \to \text{Min}_{\alpha, \beta}$$
Từ điều kiện cực trị $\frac{\partial \langle e^2 \rangle}{\partial \alpha} = 0$ và $\frac{\partial \langle e^2 \rangle}{\partial \beta} = 0$, các hệ số $\alpha$ và $\beta$ được xác định thông qua các toán tử lấy trung bình của các hàm dịch chuyển và vận tốc:
$$\alpha = \frac{\langle g(X, \dot{X})\dot{X} \rangle \langle X^2 \rangle - \langle g(X, \dot{X})X \rangle \langle X\dot{X} \rangle}{\langle X^2 \rangle \langle \dot{X}^2 \rangle - \langle X\dot{X} \rangle^2}, \quad \beta = \frac{\langle g(X, \dot{X})X \rangle \langle \dot{X}^2 \rangle - \langle g(X, \dot{X})\dot{X} \rangle \langle X\dot{X} \rangle}{\langle X^2 \rangle \langle \dot{X}^2 \rangle - \langle X\dot{X} \rangle^2}$$
2. Định nghĩa và tính chất của trung bình có trọng số
Trung bình cổ điển trên một chu kỳ $T = 2\pi/\omega$ được xác định bởi:
$$\langle f(t) \rangle_c = \frac{1}{T} \int_0^T f(t) dt$$
Hạn chế của trung bình cổ điển là $\langle \cos(nt) \rangle_c = \langle \sin(nt) \rangle_c = 0$ ($n \in \mathbb{N}^*$), dẫn đến việc mất thông tin của hàm điều hòa sau khi lấy trung bình.
Để khắc phục, định nghĩa trung bình có trọng số của hàm $x(t)$ theo GS. Nguyễn Đông Anh [2] được đưa vào:
$$\langle x(t) \rangle_w = \int_0^\infty h(t) x(t) dt$$
Hàm trọng số $h(t)$ được lựa chọn 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^{-st}$:
$$h(t) = s^2 2t e^{-st}, \quad s > 0$$
Trong đó $s$ là tham số điều chỉnh. Hàm $h(t)$ đạt cực đại tại $t_{\max} = 1/s$. Khi gán $t_{\max} = T/n = 2\pi/(n\omega)$, ta thu được mối liên hệ $s = n\omega/(2\pi)$. Với $n=1$, $h(t)$ đạt cực đại sau một chu kỳ; $n=2$ sau nửa chu kỳ; $n=4$ sau 1/4 chu kỳ; và khi $s \to 0$ ($n=0$), trung bình có trọng số trở về trung bình cổ điển.
Các tính chất quan trọng của trung bình có trọng số:
- Tính tuyến tính: $\langle C_1 x_1(t) + C_2 x_2(t) \rangle_w = C_1 \langle x_1(t) \rangle_w + C_2 \langle x_2(t) \rangle_w$.
- Liên hệ với phép biến đổi Laplace: Giá trị trung bình có trọng số là ảnh của hàm $s^2 t x(t)$ qua biến đổi Laplace:
$$\langle x(t) \rangle_w = s^2 \int_0^\infty 2t e^{-st} x(t) dt = -s^2 \frac{d}{ds} \mathcal{L}{x(t)}(s)$$
Bảng giá trị trung bình có trọng số của một số hàm điều hòa cơ bản được luận án thiết lập:
| STT |
Hàm điều hòa |
Biểu thức trung bình có trọng số $\langle \cdot \rangle_w$ |
| 1 |
$\cos(nt)$ |
$\frac{s^2(s^2 - n^2)}{(s^2 + n^2)^2}$ |
| 2 |
$\sin(nt)$ |
$\frac{2sn s^2}{(s^2 + n^2)^2}$ |
| 3 |
$\cos^2(t)$ |
$\frac{s^4 + 2s^2 + 8}{(s^2 + 4)^2}$ |
| 4 |
$\sin^2(t)$ |
$\frac{6s^2 + 8}{(s^2 + 4)^2}$ |
| 5 |
$\cos^4(t)$ |
$\frac{s^8 + 28s^6 + 248s^4 + 416s^2 + 1536}{(s^2 + 4)^2(s^2 + 16)^2}$ |
| 6 |
$\sin^4(t)$ |
$\frac{24(5s^4 + 60s^2 + 64)}{(s^2 + 4)^2(s^2 + 16)^2}$ |
| 7 |
$\cos^6(t)$ |
$\frac{s^{12} + 94s^{10} + 3168s^8 + 45712s^6 + 282496s^4 + 440064s^2 + 1658880}{(s^2 + 4)^2(s^2 + 16)^2(s^2 + 36)^2}$ |
| 8 |
$\sin^6(t)$ |
$\frac{720(7s^6 + 2805s^4 + 2352s^2 + 2304)}{(s^2 + 4)^2(s^2 + 16)^2(s^2 + 36)^2}$ |
3. Phương pháp số và công cụ kiểm chứng
Kết quả giải tích của luận án được kiểm chứng độc lập bằng phương pháp số Runge-Kutta bậc 4, nghiệm giải tích chính xác (tích phân elliptic) và các kết quả công bố của các tác giả khác (phương pháp EBM, MAFF, Min-Max). Các công cụ phần mềm hỗ trợ tính toán gồm Maple và Matlab.
Nội dung chính theo từng chương
Chương 1: Tổng quan
Chương 1 trình bày tổng quan về các hiện tượng dao động phi tuyến trong tự nhiên và kỹ thuật từ quy mô vĩ mô đến vi mô. Chương này điểm lại có hệ thống các phương pháp giải tích gần đúng (nhiễu, cân bằng điều hòa, khai triển tham số, cân bằng năng lượng), phương pháp tuyến tính hóa tương đương và các hướng phát triển, phương pháp trung bình có trọng số cùng các ứng dụng trong thống kê, xử lý tín hiệu và lý thuyết mờ. Đồng thời, chương tổng kết tình hình nghiên cứu dao động dầm micro/nano trên thế giới và các hướng nghiên cứu dao động phi tuyến trong nước. Từ đó, chương xác định mục tiêu và 2 định hướng nghiên cứu chính của luận án.
Chương 2: Phương pháp tuyến tính hóa tương đương cho hệ dao động tiền định và trung bình có trọng số
Chương 2 tập trung xây dựng cơ sở toán học cho phương pháp đề xuất. Tác giả trình bày quy trình tuyến tính hóa hệ tiền định dựa trên tiêu chuẩn cực tiểu sai số trung bình bình phương, phân tích nguyên nhân gây sai số của trung bình cổ điển khi triệt tiêu các thành phần điều hòa, và định nghĩa hình thức toán học của trung bình có trọng số với hàm trọng số $h(t) = s^2 2t e^{-st}$. Chương chứng minh tính chất tuyến tính, tính chất hội tụ qua chuỗi Fourier và mối liên hệ bản chất với phép biến đổi Laplace, đồng thời cung cấp bảng giải tích tường minh các tích phân trung bình có trọng số cho các hàm điều hòa cơ bản bậc 1 đến bậc 6.
Chương 3: Dao động phi tuyến của hệ một bậc tự do
Chương 3 áp dụng phương pháp tuyến tính hóa tương đương kết hợp trung bình có trọng số để khảo sát các hệ dao động phi tuyến một bậc tự do không cản:
- Dao động Duffing tổng quát:
Phương trình vi phân:
$$\ddot{X} + \alpha_1 X + \alpha_3 X^3 + \alpha_5 X^5 + \dots + \alpha_{2n+1} X^{2n+1} = 0, \quad X(0) = A, , \dot{X}(0) = 0$$
Phương trình tuyến tính tương đương $\ddot{X} + \omega^2 X = 0$ với nghiệm $X = A\cos(\omega t)$ dẫn đến tần số xấp xỉ:
$$\omega^2 = \alpha_1 + \alpha_3 \frac{\langle X^4 \rangle_w}{\langle X^2 \rangle_w} + \alpha_5 \frac{\langle X^6 \rangle_w}{\langle X^2 \rangle_w} + \dots + \alpha_{2n+1} \frac{\langle X^{2n+2} \rangle_w}{\langle X^2 \rangle_w}$$
- Duffing bậc 3 ($n=1$): Tần số giải tích thu được phụ thuộc vào $s$ và biên độ $A$. Khi chọn $s=2$, sai số tương đối của luận án duy trì ở mức 1.5% khi $A$ tăng, trong khi sai số của phương pháp cân bằng năng lượng EBM (Younesian [35]) lên tới 2.2%. Khi $\alpha_1=10, \alpha_3=10, A=1$, nghiệm chính xác là $\omega = 4.1672$, tần số luận án đạt độ chính xác cao tại vùng $s \in [0.5, 2]$.
- Duffing bậc 5 ($n=2$): Với $s=2$, tần số xấp xỉ là $\omega^2 = \alpha_1 + 0.7208 \alpha_3 A^2 + 0.5401 \alpha_5 A^4$. Khi $\alpha_1=1, \alpha_3=10, \alpha_5=100, A=10$, tần số chính xác là 751.6951; sai số luận án là 1.52% so với 2.26% của EBM.
- Duffing bậc 7 ($n=3$): Tần số xấp xỉ với $s=2$ là $\omega^2 = \alpha_1 + 0.7208 \alpha_3 A^2 + 0.5401 \alpha_5 A^4 + 0.4132 \alpha_7 A^6$.
- Duffing bậc cao ($n=4, 5, 6$ tương ứng bậc 9, 11, 13): Nghiệm giải tích của luận án được kiểm chứng với nghiệm số Runge-Kutta bậc 4, cho thấy mức độ bám sát cao của đường cong pha và dịch chuyển theo thời gian.
- Dao động phi tuyến mở rộng (Generalized nonlinear oscillator):
Phương trình tổng quát:
$$\ddot{u} + \alpha u + \beta u^m + \frac{\mu \ddot{u} u^n + \delta u^p + \gamma u^{m+p}}{1 + u^p} = 0, \quad u(0)=A, , \dot{u}(0)=0$$
Luận án hệ thống hóa các trường hợp riêng của phương trình này:
| STT |
Trường hợp dao động |
Tham số tương ứng |
Tác giả & Tài liệu trích dẫn |
| 1 |
Duffing bậc 3 |
$m=3, \mu=\delta=\gamma=0$ |
Bayat và cộng sự [43] |
| 2 |
Duffing bậc 5 |
$m=5, \mu=\delta=\gamma=0$ |
Sedighi và Reza [50] |
| 3 |
Duffing tổng quát |
$m=2k+1, \mu=\delta=\gamma=0$ |
Younesian và cộng sự [35] |
| 4 |
Dao động trong vật lý plasma |
$\alpha=0, \beta=0, \gamma=1, m=1$ |
Mickens [19] |
| 5 |
Động lực học hệ thống treo |
$\alpha=0, \beta=1, m=3/2, \mu=\delta=\gamma=0$ |
Zhu & Ishitoby [130], Cveticanin [131] |
| 6 |
Hệ động lực phẳng |
$\alpha=0, \beta=1, m=1/3, \mu=\delta=\gamma=0$ |
Mickens [16–18], Ozis [31], Beléndez [13] |
| 7 |
Dao động Duffing-điều hòa |
$\alpha=0, \beta=1, m=3, \mu=1, \delta=1, \gamma=1, n=1, p=2$ |
Momeni và cộng sự [36], Fan [132] |
| 8 |
Dao động có số mũ hữu tỉ |
$\alpha=a, \beta=b, m=3, \mu=1, \delta=0, \gamma=c$ |
J.H. He [33] |
- Dao động Duffing-điều hòa: Phương trình $\ddot{u} + \frac{\ddot{u}u + u + u^3}{1+u^2} = 0$, tần số xấp xỉ xác định dạng đóng $\omega = \sqrt{\frac{1 + 0.7208 A^2}{1 + 0.7208 A^2}} \dots$, kiểm chứng bằng phương pháp Runge-Kutta bậc 4 và phương pháp MAFF (Fan [132]).
- Dao động Duffing thế năng giếng đôi: Phương trình $\ddot{u} - u + u^3 = 0$ với thế năng $U(u) = -u^2/2 + u^4/4$, có 3 điểm cân bằng ($u=0$ không ổn định, $u=\pm 1$ ổn định). Với biên độ $A \ge \sqrt{2}$, dao động đối xứng qua các vị trí cân bằng, tần số xấp xỉ thu được là $\omega = \sqrt{-1 + 0.7208 A^2}$.
Chương 4: Dao động phi tuyến của dầm micro và nano
Chương 4 mở rộng phương pháp cho hệ liên tục thông qua hai bài toán cụ thể:
- Dao động phi tuyến của dầm micro tựa trên nền đàn hồi ba lớp theo lý thuyết ứng suất cặp sửa đổi: giải quyết ảnh hưởng của tham số chiều dài vật liệu, tham số nền đàn hồi, đồng thời làm rõ vai trò của hai tham số mới là tỉ số độ cứng chống uốn (flexural rigidity ratio) và tỉ số độ mảnh (slenderness ratio).
- Dao động phi tuyến của dầm nano chịu lực tĩnh điện trong hệ vi cơ điện tử MEMS/NEMS theo lý thuyết độ dốc biến dạng phi cục bộ: phân tích tương tác phi tuyến giữa lực tĩnh điện, hiệu ứng phi cục bộ và gradient biến dạng lên tần số dao động phi tuyến của hệ.
Kết quả và những đóng góp mới
Luận án đã đạt được các kết quả và đóng góp mới sau:
1. Đóng góp về mặt lý luận
- Thiết lập hoàn chỉnh khung phương pháp phân tích giải tích gần đúng dựa trên sự kết hợp giữa phương pháp tuyến tính hóa tương đương và trung bình có trọng số cho các hệ dao động phi tuyến tiền định không cản.
- Làm sáng tỏ bản chất toán học của thuật toán trung bình có trọng số với hàm trọng số $h(t) = s^2 2t e^{-st}$, thiết lập mối liên hệ trực tiếp với phép biến đổi Laplace, và xây dựng bảng công thức tích phân giải tích cho các hàm điều hòa cơ bản.
- Chứng minh rằng trung bình có trọng số bảo toàn được đặc tính tuần hoàn và thông tin biên độ - tần số trong nửa chu kỳ đầu tiên, khắc phục hiện tượng triệt tiêu thông tin của phương pháp trung bình cổ điển.
2. Đóng góp về mặt thực tiễn
- Cung cấp các biểu thức giải tích tường minh dạng đóng cho quan hệ biên độ - tần số của nhiều hệ dao động phi tuyến một bậc tự do phức tạp (Duffing bậc 3 đến bậc 13, Duffing thế năng giếng đôi, Duffing-điều hòa, hệ có số mũ hữu tỉ và phân thức). Kết quả cho thấy độ chính xác cao hơn phương pháp cân bằng năng lượng (EBM) khi biên độ dao động lớn.
- Giải quyết hoàn chỉnh bài toán dao động của dầm micro trên nền đàn hồi ba lớp theo lý thuyết ứng suất cặp sửa đổi, bổ sung khảo sát định lượng ảnh hưởng của tỉ số độ cứng chống uốn và tỉ số độ mảnh.
- Xây dựng lời giải giải tích cho bài toán dao động dầm nano chịu tác dụng của lực tĩnh điện trong MEMS/NEMS theo lý thuyết độ dốc biến dạng phi cục bộ, mô tả chính xác hiệu ứng phụ thuộc kích thước hình học.
3. Kiến nghị và giải pháp đề xuất
- Kiến nghị ứng dụng quy trình tuyến tính hóa tương đương với trung bình có trọng số vào tính toán thiết kế sơ bộ các phần tử cảm biến và cơ cấu chấp hành vi mô (MEMS/NEMS) nhằm giảm tải chi phí tính toán mô phỏng số.
- Đề xuất mở rộng nghiên cứu sang các hệ dao động có cản phi tuyến, hệ chịu kích động điều hòa cưỡng bức hoặc kích động ngẫu nhiên.
Hạn chế và hướng nghiên cứu tiếp
- Hạn chế: Các khảo sát trong luận án chỉ giới hạn ở các hệ dao động tự do không cản dạng tiền định. Việc xác định giá trị tối ưu của tham số điều chỉnh $s$ ($s=2$) được thực hiện dựa trên khảo sát thực nghiệm số so sánh với nghiệm chính xác hoặc Runge-Kutta cho từng lớp phương trình, chưa xây dựng được tiêu chuẩn tối ưu hóa giải tích thuần túy để chọn $s$ tự động.
- Hướng nghiên cứu tiếp:
- Mở rộng phương pháp cho các hệ dao động phi tuyến có cản (cản nhớt, cản phi tuyến, cản Coulomb).
- Khảo sát hệ dao động chịu kích động ngoài tuần hoàn, kích động tham số và kích động ngẫu nhiên.
- Mở rộng cho hệ dao động nhiều bậc tự do và các kết cấu dầm/tấm kích thước micro/nano có hình học phức tạp.
Giá trị tham khảo
- Đối tượng thụ hưởng: Luận án là tài liệu tham khảo cho các nghiên cứu sinh, học viên cao học, giảng viên và kỹ sư chuyên ngành Cơ kỹ thuật, Cơ học vật rắn biến dạng, Toán ứng dụng và Kỹ thuật vi hệ thống (MEMS/NEMS).
- Nội dung giá trị:
- Bảng tích phân trung bình có trọng số qua biến đổi Laplace ở Chương 2 là công cụ tính toán thuận tiện cho các bài toán xấp xỉ giải tích.
- Quy trình xử lý các phương trình phi tuyến chứa số hạng hữu tỉ, vô tỉ ở Chương 3 và kỹ thuật mô hình hóa dầm liên tục kích thước micro/nano ở Chương 4 cung cấp khung tham chiếu rõ ràng cho các nghiên cứu cơ học tính toán kết cấu vi mô.
Câu hỏi thường gặp
1. Cơ sở toán học của việc lựa chọn hàm trọng số $h(t) = s^2 2t e^{-st}$ là gì?
Hàm $h(t)$ được xây dựng từ 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^{-st}$. Cấu trúc này đảm bảo $h(t)$ đạt cực trị tại $t_{\max} = 1/s$. Khi điều chỉnh tham số $s = n\omega/(2\pi)$, hàm trọng số sẽ tập trung trọng số vào các khoảng thời gian xác định (sau 1 chu kỳ, nửa chu kỳ, hoặc 1/4 chu kỳ), giúp giữ lại thông tin tuần hoàn của hàm điều hòa mà không bị triệt tiêu về 0 như khi lấy tích phân trung bình cổ điển.
2. Mối liên hệ giữa trung bình có trọng số và phép biến đổi Laplace được biểu diễn như thế nào?
Giá trị trung bình có trọng số của hàm tuần hoàn $x(t)$ là ảnh của hàm $s^2 t x(t)$ qua phép biến đổi Laplace: $\langle x(t) \rangle_w = s^2 \int_0^\infty 2t e^{-st} x(t) dt = -s^2 \frac{d}{ds} \mathcal{L}{x(t)}(s)$. Nhờ mối liên hệ này, việc tính toán giá trị trung bình của các hàm điều hòa bậc cao $\cos^n(t)$, $\sin^n(t)$ được chuyển thành các phép đạo hàm đại số của hàm Laplace.
3. Vì sao phương pháp trong luận án cho kết quả chính xác hơn phương pháp cân bằng năng lượng (EBM) khi biên độ dao động $A$ lớn?
Trong các hệ dao động Duffing bậc 3, 5, 7, khi biên độ $A$ tăng cao, tính phi tuyến của hệ trở nên mạnh mẽ. Phương pháp trung bình có trọng số với tham số điều chỉnh $s=2$ giúp phân bổ trọng số hợp lý vào nửa chu kỳ dao động đầu tiên, làm giảm độ lệch cực đại và giữ lại đặc trưng phi tuyến của hệ tốt hơn so với tiêu chuẩn tích phân biến phân của EBM, từ đó duy trì sai số tương đối ở mức khoảng 1.5% (so với 2.2% - 2.26% của EBM).
4. Luận án đã giải quyết những khoảng trống cụ thể nào trong bài toán dầm micro và nano ở Chương 4?
Luận án đã giải quyết hai bài toán chưa được khảo sát đầy đủ trước đó: (1) Khảo sát ảnh hưởng của tỉ số độ cứng chống uốn và tỉ số độ mảnh đến đáp ứng dao động của dầm micro trên nền đàn hồi ba lớp theo lý thuyết ứng suất cặp sửa đổi; (2) Thiết lập và phân tích bài toán dao động phi tuyến của dầm nano đặt giữa hai bản tích điện trong hệ MEMS/NEMS dựa trên lý thuyết độ dốc biến dạng phi cục bộ nhằm quan sát hiệu ứng kích thước hình học.
5. Phương trình dao động phi tuyến mở rộng trong luận án bao quát những mô hình kỹ thuật cụ thể nào?
Phương trình dao động phi tuyến mở rộng trong luận án bao quát 8 mô hình vật lý và kỹ thuật khác nhau, bao gồm: dao động Duffing bậc 3, Duffing bậc 5, Duffing bậc tổng quát, dao động trong vật lý plasma, mô hình động lực học hệ thống treo ô tô, hệ động lực phẳng, dao động Duffing-điều hòa và dao động có số mũ hữu tỉ.
Kết luận
Luận án đã phát triển thành công phương pháp tuyến tính hóa tương đương kết hợp với trung bình có trọng số để phân tích đáp ứng của các hệ dao động phi tuyến tiền định không cản. Bản chất toán học của thuật toán trung bình có trọng số và mối liên hệ với phép biến đổi Laplace đã được chứng minh và lập bảng công thức giải tích hoàn chỉnh. Độ chính xác của phương pháp được khẳng định qua việc giải quyết hàng loạt bài toán dao động một bậc tự do từ Duffing bậc cao, Duffing-điều hòa, giếng đôi đến các hệ liên tục dầm micro và nano trong MEMS/NEMS. Các kết quả giải tích dạng đóng thu được cung cấp công cụ tính toán hiệu quả cho việc nghiên cứu và thiết kế các hệ dao động phi tuyến trong kỹ thuật.
