Tổng quan luận án
Vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM) là loại vật liệu composite tiên tiến được hình thành từ sự kết hợp của hai hay nhiều pha vật liệu (thường là kim loại và gốm), trong đó thành phần và tính chất cơ lý biến đổi liên tục theo không gian (thường là theo chiều dày). Sự biến đổi trơn này giúp loại bỏ mặt phân chia rõ rệt giữa các lớp, từ đó khắc phục nhược điểm tập trung ứng suất tiếp và hiện tượng tách lớp thường gặp ở composite phân lớp truyền thống. Các kết cấu vỏ tròn xoay, vỏ liên hợp và vỏ có gân gia cường chế tạo từ FGM được ứng dụng rộng rãi trong các ngành công nghiệp kỹ thuật cao như hàng không vũ trụ, tên lửa đẩy, thân tàu ngầm, động cơ nhiệt và lò phản ứng hạt nhân. Trong quá trình làm việc, các kết cấu này thường xuyên tiếp xúc với môi trường xung quanh (được mô hình hóa qua nền đàn hồi) và chịu tải trọng động, dẫn đến nguy cơ dao động mạnh, cộng hưởng và phá hủy kết cấu.
Khoảng trống nghiên cứu được xác định: Các công trình trước đây về dao động của vỏ FGM chủ yếu tập trung vào các dạng hình học cơ bản (vỏ trụ, nón cụt đơn lẻ) bằng phương pháp giải tích hoặc phương pháp phần tử hữu hạn (PTHH). Chưa có công trình nào nghiên cứu toàn diện về dao động tự do của các hệ vỏ FGM liên hợp phức tạp (vỏ nón cụt - trụ, nón cụt - vành tròn - trụ, vỏ dạng bậc và vỏ trụ có gân gia cường dạng vành) tương tác với nền đàn hồi Winkler-Pasternak bằng phương pháp Phần tử liên tục (PTLT).
Mục tiêu nghiên cứu của luận án được xác định gồm ba nội dung:
- Tìm lời giải số tin cậy bằng phương pháp phần tử liên tục (PTLT) về tần số dao động tự do cho các kết cấu vỏ FGM đối xứng trục đơn giản và phức tạp (vỏ nón cụt, vỏ trụ, vành tròn, vỏ dạng bậc, vỏ liên hợp) có và không tương tác với nền đàn hồi Winkler-Pasternak; bổ sung và hoàn thiện thư viện phần tử liên tục cho vật liệu FGM.
- Phát triển thuật toán ghép nối phần tử liên tục cho các kết cấu vỏ đối xứng trục phức tạp làm bằng vật liệu FGM, bao gồm: nón bậc, trụ bậc, vành bậc, vỏ liên hợp trụ - nón cụt, nón cụt - vành - trụ, và vỏ trụ có gân gia cường dạng vành.
- Phân tích và khảo sát định lượng ảnh hưởng của điều kiện biên, thông số hình học, thuộc tính vật liệu FGM, cấu hình phân bố vật liệu và các hệ số nền đàn hồi đến tần số và dạng dao động tự do của các kết cấu vỏ FGM.
Đối tượng và phạm vi nghiên cứu:
- Đối tượng nghiên cứu: Kết cấu vỏ tròn xoay đơn (vỏ trụ, nón cụt, vành tròn); kết cấu vỏ liên hợp dạng bậc đối xứng trục (trụ bậc, nón bậc, vành bậc); kết cấu vỏ liên hợp phức tạp (trụ - nón cụt, nón cụt - vành - trụ, vỏ trụ có gân gia cường dạng vành) làm bằng vật liệu FGM có và không tương tác với nền đàn hồi Winkler-Pasternak.
- Phạm vi nghiên cứu: Bài toán dao động tự do trong giới hạn đàn hồi tuyến tính, vật liệu trực hướng, biến dạng bé; sử dụng lý thuyết vỏ biến dạng cắt bậc nhất (First-order Shear Deformation Theory - FSDT) của Reissner-Mindlin có xét đến ảnh hưởng của biến dạng cắt ngang và quán tính quay; nền đàn hồi tuyến tính hai hệ số Pasternak và một hệ số Winkler.
Tổng quan tài liệu và vị trí của luận án
Luận án đã hệ thống hóa và phân tích các nhóm công trình nghiên cứu trong nước và quốc tế liên quan đến dao động của kết cấu vỏ:
-
Nghiên cứu quốc tế về vỏ FGM đơn lẻ và tương tác với nền đàn hồi: Sheng và Wang [10, 13] nghiên cứu giải tích dao động phi tuyến và nhiệt - cơ của vỏ trụ FGM bằng nguyên lý Hamilton và lý thuyết von Karman. Xie và cộng sự [11] dùng phương pháp truyền sóng Haar và lý thuyết FSDT khảo sát dao động vỏ nón và vành FGM. Darabi và cộng sự [12] phân tích cộng hưởng tham số của vỏ trụ FGM. Shen và cộng sự [20, 21] dùng lý thuyết biến dạng cắt bậc cao (HSDT) phân tích trạng thái sau tới hạn của vỏ trụ FGM trên nền đàn hồi. Sheng và Wang [22] khảo sát ảnh hưởng tải nhiệt đến dao động vỏ trụ FGM trên nền đàn hồi tuyến tính. Bagherizadeh và cộng sự [23] nghiên cứu ổn định vỏ trụ FGM trên nền Pasternak theo HSDT. Najafov và cộng sự [24] nghiên cứu dao động xoắn của vỏ trụ FGM trên nền đàn hồi. Nhóm nghiên cứu của Sofiyev cùng các cộng sự Deniz, Kuruoglu, Schnack [14-18, 25-29] đã khảo sát ứng xử dao động và ổn định động của vỏ nón cụt FGM trên nền đàn hồi hai hệ số sử dụng phương pháp Galerkin, lý thuyết vỏ Donnell và thuật toán Runge-Kutta.
-
Nghiên cứu về kết cấu vỏ liên hợp và vỏ dạng bậc: Galletly và Mistry [30] sử dụng sai phân hữu hạn cho vỏ trụ nối kết cấu tròn xoay. Irie và cộng sự [31] dùng phương pháp ma trận truyền khảo sát vỏ nón - trụ và vỏ trụ - tấm cong dựa trên phương trình Flügge. Tavakoli và Singh [32] dùng không gian trạng thái; Redekop [33] dùng phương pháp cầu phương vi phân; Efraim và Eisenberger [34] xây dựng ma trận độ cứng động lực cho vỏ liên hợp nón - trụ đẳng hướng. Qu và cộng sự [38, 39] dùng phương pháp biến phân cho vỏ nón - trụ - vành và nón - trụ - cầu - vành tăng cứng. Ma và cộng sự [40] dùng phương pháp Fourier-Ritz cho vỏ nón - trụ. Patel và cộng sự [41], Kamat [42] dùng PTHH phân tích dao động và mất ổn định của vỏ nón - trụ composite. Kouchakzadeh và Shakouri [43] nghiên cứu vỏ composite nón - nón và nón - trụ theo lý thuyết Donnell và nguyên lý Hamilton.
-
Nghiên cứu tại Việt Nam: Đào Huy Bích và cộng sự [44-52] khảo sát ổn định tĩnh, động phi tuyến và dao động của tấm, vỏ trụ, vỏ thoải FGM có gân gia cường trên nền đàn hồi bằng phương pháp giải tích và kỹ thuật san gân. Đào Văn Dũng và cộng sự [53-56] nghiên cứu ổn định động vỏ trụ FGM có gân và tấm FGM có hệ số Poisson thay đổi theo độ dày. Nguyễn Đình Đức và cộng sự [67-70] nghiên cứu ổn định phi tuyến của vỏ và mảnh trụ FGM chịu nén và tải nhiệt trên nền đàn hồi. Trần Minh Tú, Trần Hữu Quốc và cộng sự [72-77] phân tích tĩnh và động panel trụ, vỏ thoải FGM chịu nhiệt và đặt trên nền đàn hồi theo FSDT. Về phương pháp ma trận độ cứng động lực và PTLT: Đỗ Văn Hiến, Nguyễn Xuân Hùng [78], Nguyễn Tiến Khiêm, Trần Văn Liên [79] áp dụng cho kết cấu dầm; Trần Ích Thịnh, Nguyễn Mạnh Cường [57-62, 94] phát triển PTLT cho vỏ trụ composite chứa lỏng; Vũ Quốc Hiến [95] và Lê Thị Bích Nam [96] mở rộng PTLT cho các kết cấu vỏ composite dạng bậc và vỏ ghép nối.
Khoảng trống luận án lựa chọn giải quyết: Mặc dù các nghiên cứu về FGM rất đa dạng, việc áp dụng phương pháp Phần tử liên tục (Dynamic Stiffness Matrix Method) để giải quyết bài toán dao động tự do của các hệ vỏ liên hợp FGM phức tạp (nón - trụ, nón - vành - trụ, vỏ trụ có gân gia cường dạng vành) đặt trên nền đàn hồi Winkler-Pasternak chưa từng được thực hiện. Luận án lấp đầy khoảng trống này bằng cách xây dựng hệ phương trình PTLT, ma trận độ cứng động lực cấp phần tử và hệ thống, kết hợp phương pháp đường cong đáp ứng để tính toán các tần số dao động riêng.
Cơ sở lý thuyết và phương pháp nghiên cứu
1. Mô hình vật liệu FGM
Vật liệu FGM cấu tạo từ hai thành phần chính là kim loại ($m$) và gốm ($c$). Các đặc tính cơ học gồm mô đun đàn hồi kéo nén $E(z)$, khối lượng riêng $\rho(z)$ và hệ số Poisson $\nu(z)$ biến thiên liên tục theo chiều dày $z$ (với $-h/2 \le z \le h/2$):
$$\begin{aligned}
E(z) &= (E_c - E_m)V_c + E_m \
\rho(z) &= (\rho_c - \rho_m)V_c + \rho_m \
\nu(z) &= (\nu_c - \nu_m)V_c + \nu_m
\end{aligned}$$
Hàm tỉ lệ thể tích pha gốm $V_c$ tuân theo quy luật phân bố tổng quát 4 tham số $(a, b, c, p)$:
$$V_c = \left[ 1 - a\left(\frac{1}{2} + \frac{z}{h}\right) + b\left(\frac{1}{2} + \frac{z}{h}\right)^c \right]^p$$
Khi $p = 0$, vật liệu trở thành gốm đồng nhất ($E = E_c, \rho = \rho_c, \nu = \nu_c$); khi $p \to \infty$, vật liệu chuyển thành kim loại đồng nhất ($E = E_m, \rho = \rho_m, \nu = \nu_m$). Khi $a=1, b=0$, hàm chuyển về dạng quy luật lũy thừa kinh điển P-FGM.
| Vật liệu thành phần |
$E$ ($\text{N/m}^2$) |
$\nu$ |
$\alpha$ ($\text{K}^{-1}$) |
$\kappa$ ($\text{W/mK}$) |
$\rho$ ($\text{kg/m}^3$) |
| Kim loại |
|
|
|
|
|
| Thép không gỉ (SUS 304) |
$207{,}7 \times 10^9$ |
$0{,}3177$ |
$12{,}46 \times 10^{-6}$ |
$15{,}379$ |
$8166$ |
| Nhôm (Al) |
$70{,}0 \times 10^9$ |
$0{,}3000$ |
$23{,}00 \times 10^{-6}$ |
$204{,}000$ |
$2707$ |
| Hợp kim Titan (Ti6Al4V) |
$105{,}7 \times 10^9$ |
$0{,}2980$ |
$6{,}90 \times 10^{-6}$ |
$18{,}100$ |
$4429$ |
| Gốm |
|
|
|
|
|
| Xi-lích nitrit ($\text{Si}_3\text{N}_4$) |
$322{,}2 \times 10^9$ |
$0{,}2400$ |
$3{,}20 \times 10^{-6}$ |
$13{,}723$ |
$2370$ |
| Ôxít Zi-ri-cô-ni ($\text{ZrO}_2$) |
$151{,}0 \times 10^9$ |
$0{,}3000$ |
$10{,}00 \times 10^{-6}$ |
$2{,}090$ |
$3000$ |
| Ôxít nhôm ($\text{Al}_2\text{O}_3$) |
$380{,}0 \times 10^9$ |
$0{,}2600$ |
$7{,}20 \times 10^{-6}$ |
$10{,}400$ |
$3750$ |
2. Mô hình nền đàn hồi Winkler-Pasternak
Tương tác giữa nền đàn hồi và mặt vỏ được mô tả bằng mô hình hai hệ số Pasternak:
$$p_z = k_w w - k_p \nabla^2 w = k_w w - k_p \left( \frac{\partial^2 w}{\partial x^2} + \frac{1}{R^2} \frac{\partial^2 w}{\partial \theta^2} \right)$$
Trong đó $k_w$ ($\text{N/m}^3$) là hệ số độ cứng đàn hồi Winkler mô tả phản lực lò xo tuyến tính vuông góc; $k_p$ ($\text{N/m}$) là hệ số cắt Pasternak mô tả tương tác trượt giữa các phần tử nền. Khi $k_p = 0$, mô hình trở về dạng nền đàn hồi Winkler một hệ số.
3. Lý thuyết vỏ và Phương pháp Phần tử liên tục (PTLT)
Sử dụng lý thuyết biến dạng cắt bậc nhất của Reissner-Mindlin với trường chuyển vị:
$$\begin{aligned}
u(x, \theta, z, t) &= u_0(x, \theta, t) + z \varphi_x(x, \theta, t) \
v(x, \theta, z, t) &= v_0(x, \theta, t) + z \varphi_\theta(x, \theta, t) \
w(x, \theta, z, t) &= w_0(x, \theta, t)
\end{aligned}$$
Hệ số hiệu chỉnh cắt được chọn là $k = 5/6$ (hoặc $f = 5/6$). Các quan hệ nội lực và độ cứng màng $[A]$, độ cứng uốn $[D]$, độ cứng tương tác màng - uốn $[B]$ và độ cứng cắt $[F]$ được tích phân dọc theo chiều dày vỏ $h$.
Quy trình thiết lập PTLT:
- Chọn véc tơ trạng thái kích thước $10 \times 1\$:
$$y = { u_0, v_0, w_0, \varphi_x, \varphi_\theta, N_x, N_{x\theta}, Q_x, M_x, M_{x\theta} }^T$$
- Biểu diễn nghiệm dao động tuần hoàn theo phương chu vi với số sóng chu vi $m$ và tần số góc $\omega$:
$${ u_0, w_0, \varphi_x, N_x, Q_x, M_x } \sim \cos(m\theta) e^{i\omega t}, \quad { v_0, \varphi_\theta, N_{x\theta}, M_{x\theta} } \sim \sin(m\theta) e^{i\omega t}$$
- Đưa hệ phương trình chuyển động và quan hệ nội lực - chuyển vị về dạng phương trình vi phân bậc nhất theo biến tọa độ đường sinh $x$:
$$\frac{d y_m(x, \omega)}{dx} = A_m(x, \omega) y_m(x, \omega)$$
với $A_m(x, \omega)$ là ma trận trạng thái kích thước $10 \times 10$.
- Tích phân phương trình vi phân để xác định ma trận truyền $T(x, \omega)$, từ đó suy ra quan hệ lực nút - chuyển vị nút của phần tử qua ma trận độ cứng động lực:
$$F(\omega) = [K(\omega)] U(\omega)$$
- Kỹ thuật xác định tần số riêng: Sử dụng Phương pháp Đường cong đáp ứng. Tác giả đặt một lực đơn vị điều hòa $F = 1 \cdot e^{i\omega t}$ tại một điểm biên, giải hệ phương trình động lực học $U = [K(\omega)]^{-1} F$ qua một dải tần số $\omega$, và xác định các tần số dao động tự do tại các đỉnh cộng hưởng nơi biên độ chuyển vị $w$ đạt cực đại.
Nội dung chính theo từng chương
Chương 1: Tổng quan về vấn đề nghiên cứu
Chương 1 giới thiệu tổng quan về vật liệu FGM, cấu trúc Maxel, quy luật biến thiên hàm tỉ lệ thể tích 4 tham số và các ứng dụng tiêu biểu của vỏ FGM trong công nghiệp. Trình bày tổng quan tài liệu nghiên cứu trong và ngoài nước về dao động của vỏ tròn xoay, vỏ liên hợp kim loại, composite và FGM có/không có nền đàn hồi. Chương này phân tích chi tiết cơ sở toán học của phương pháp PTLT (ma trận độ cứng động lực), so sánh với các phương pháp giải tích, PTHH, thuật toán Wittrick-Williams và nêu bật lý do lựa chọn phương pháp đường cong đáp ứng để xác định tần số riêng trong luận án.
Chương 2: Nghiên cứu dao động tự do của kết cấu vỏ nón cụt, trụ và vành tròn làm bằng FGM có và không tương tác với nền đàn hồi
Chương 2 thiết lập hệ phương trình vi phân tổng quát cho vỏ nón cụt FGM trên nền đàn hồi Pasternak, từ đó suy biến cho trường hợp vỏ trụ ($\alpha = 0^\circ$) và vành tròn phẳng ($\alpha = 90^\circ$). Ma trận trạng thái $A_m(10 \times 10)$ được thiết lập chi tiết.
Các kết quả tính toán và kiểm chứng:
- Vỏ trụ FGM: So sánh 3 tần số đầu tiên với kết quả công bố cho thấy sự trùng khớp với sai số không đáng kể ở các điều kiện biên C-C (Ngàm-Ngàm), F-C (Tự do-Ngàm), S-S (Tựa-Tựa).
- Vỏ nón cụt FGM: Khảo sát tần số không thứ nguyên $\Omega$ với các bộ tham số vật liệu FGM1I, FGM1II. Đánh giá ảnh hưởng của góc nửa đỉnh nón $\alpha$, tỉ số bán kính $R_2/R_1$, tỉ số $L/R_1$, chiều dày $h/R_1$ và hệ số nền ($k_w, k_p$). Khi $k_w$ tăng từ $0$ lên $5 \times 10^6 \text{ N/m}^3$ và $k_p$ tăng lên $10^5 \text{ N/m}$, tần số dao động riêng cơ bản tăng rõ rệt do độ cứng tổng thể của hệ tăng.
- Vành tròn FGM: Đánh giá 10 tần số đầu tiên với cấu hình FGM2I, FGM2II theo số mũ $p$ và điều kiện biên.
Chương 3: Nghiên cứu dao động tự do của các kết cấu vỏ FGM đối xứng trục dạng bậc có và không tương tác với nền đàn hồi
Chương 3 mở rộng mô hình PTLT cho các kết cấu vỏ tròn xoay có sự thay đổi đột ngột về kích thước hình học hoặc đặc tính vật liệu (kết cấu dạng bậc). Thiết lập điều kiện cân bằng lực và liên tục chuyển vị tại các mặt cắt tiếp giáp giữa các bậc.
Các kết quả tính toán:
- Vỏ nón cụt 4 bậc FGM: Khảo sát cấu hình FGM5I, phân tích ảnh hưởng của tỉ số độ dày giữa các bậc, góc nửa đỉnh $\alpha$, và sự thay đổi vật liệu FGM ở từng bậc riêng biệt.
- Vành tròn 2 bậc và 3 bậc FGM: So sánh 9 tần số đầu tiên, kiểm chứng tính chính xác của thuật toán ghép nối phần tử liên tục nối tiếp.
- Vỏ trụ 2 bậc và 4 bậc: So sánh tần số không thứ nguyên $\Omega = \omega R \sqrt{\rho(1-\nu^2)/E}$ của vỏ trụ kim loại hai bậc với tài liệu tham khảo dưới điều kiện biên C-C và C-F. Khảo sát vỏ trụ FGM 4 bậc trên nền đàn hồi với các giá trị $k_w = 0; 2{,}5 \times 10^6; 5 \times 10^6; 5 \times 10^7 \text{ N/m}^3$ và $k_p = 0; 1 \times 10^5; 2 \times 10^5 \text{ N/m}$, chỉ ra rằng tần số dao động chịu ảnh hưởng đồng thời bởi vị trí bậc, số mũ $p$ và độ cứng nền.
Chương 4: Nghiên cứu dao động tự do của các kết cấu vỏ liên hợp FGM đối xứng trục có và không tương tác với nền đàn hồi
Chương 4 giải quyết các bài toán vỏ liên hợp phức tạp nhất trong luận án:
- Vỏ liên hợp nón - trụ FGM: Thiết lập điều kiện cân bằng và liên tục tại mặt cắt chuyển tiếp giữa phần nón và phần trụ. Khảo sát 10 tần số đầu tiên dưới các điều kiện biên C-C, C-F, F-C với các giá trị số mũ $p$, góc $\alpha$, và hệ số nền $k_w, k_p = 5 \times 10^6$.
- Vỏ liên hợp nón - vành - trụ FGM: Ghép nối ba phần tử liên tục khác loại (nón cụt, vành phẳng, vỏ trụ). Khảo sát ảnh hưởng của chiều dày vỏ, kiểu hàm phân bố FGM3I, FGM3II và thông số nền đàn hồi Pasternak.
- Vỏ trụ FGM có gân gia cường dạng vành: Đề xuất và ứng dụng thuật toán ghép nối song song dạng chữ T cho ma trận độ cứng động lực. Gân gia cường được mô hình hóa chính xác như một phần tử vành tròn đàn hồi liên tục (thay vì làm mịn/san đều như các phương pháp giải tích truyền thống). Kết quả tính toán tần số dao động được kiểm chứng với vỏ trụ kim loại có gân và mở rộng khảo sát cho vỏ FGM có gân trên nền Winkler-Pasternak dưới các điều kiện biên F-C, C-C.
Kết quả và những đóng góp mới
1. Đóng góp về mặt phương pháp và thuật toán
- Xây dựng thành công thuật toán và bộ chương trình tính toán số trong môi trường MATLAB dựa trên phương pháp Phần tử liên tục (ma trận độ cứng động lực) để tính toán tần số dao động tự do cho các kết cấu vỏ tròn xoay FGM đơn lẻ (trụ, nón, vành) và các hệ vỏ liên hợp phức tạp có xét đến tương tác nền đàn hồi Winkler-Pasternak.
- Phát triển và ứng dụng thành công thuật toán ghép nối nối tiếp cho vỏ dạng bậc, vỏ liên hợp nón - trụ, nón - vành - trụ, và đặc biệt là thuật toán ghép nối song song dạng chữ T cho vỏ trụ FGM có gân gia cường dạng vành. Mô hình này mô tả gân như một kết cấu vành đàn hồi thực tế, khắc phục sai số của giả thiết san đều gân.
- Tận dụng ưu điểm của PTLT: số lượng phần tử sử dụng là tối thiểu (mỗi phân đoạn hình học chỉ cần một phần tử duy nhất), không đòi hỏi chia lưới dày đặc, tiết kiệm tài nguyên bộ nhớ và thời gian tính toán, đồng thời tính toán chính xác trong cả miền tần số trung bình và cao.
2. Đóng góp về mặt dữ liệu và phân tích cơ học
- Cung cấp hệ thống dữ liệu số chính xác về tần số dao động tự do và dạng dao động cho nhiều dạng kết cấu vỏ FGM chưa từng được công bố trong các tài liệu trước đây.
- Đánh giá định lượng quy luật biến đổi tần số dao động riêng theo:
- Số mũ phân bố $p$: Khi $p$ tăng, hàm lượng gốm giảm (tăng hàm lượng kim loại), độ cứng tổng thể của kết cấu giảm làm tần số dao động riêng giảm.
- Điều kiện biên: Biên ngàm (C) cho tần số dao động cao nhất, tiếp đến là biên tựa (S) và biên tự do (F).
- Nền đàn hồi: Cả hai hệ số $k_w$ và $k_p$ đều làm tăng tần số dao động riêng của vỏ, trong đó hệ số trượt $k_p$ của mô hình Pasternak có ảnh hưởng đáng kể đến các mode dao động bậc cao.
- Gân gia cường: Bố trí gân dạng vành làm tăng độ cứng cục bộ và tần số dao động riêng của vỏ trụ một cách rõ rệt.
3. Kiến nghị ứng dụng thực tiễn
- Các mô hình và chương trình tính toán của luận án có thể áp dụng trực tiếp trong quy trình tính toán, thiết kế động lực học nhằm tránh hiện tượng cộng hưởng, giảm thiểu rung động và tiếng ồn cho các chi tiết vỏ FGM trong công nghiệp hàng không, tàu thủy, lò phản ứng hạt nhân và thiết bị hóa chất.
Hạn chế và hướng nghiên cứu tiếp
1. Hạn chế của luận án
- Nghiên cứu giới hạn trong phạm vi dao động tuyến tính, biến dạng bé và vật liệu làm việc trong miền đàn hồi.
- Nền đàn hồi Winkler-Pasternak được giả thiết là liên tục và đối xứng trục, chưa xét đến các trường hợp nền đàn hồi phi tuyến, nền chịu tải cục bộ không đều hoặc nền đàn nhớt tiêu tán năng lượng.
- Luận án mới tập trung vào bài toán dao động tự do, chưa phân tích đáp ứng dao động cưỡng bức dưới các dạng tải trọng động ngẫu nhiên hoặc tải va đập theo thời gian thực.
2. Hướng phát triển tiếp theo
- Mở rộng phương pháp PTLT để phân tích dao động của kết cấu vỏ liên hợp FGM có chứa hoặc tiếp xúc với chất lỏng (bài toán tương tác kết cấu - chất lỏng).
- Nghiên cứu dao động và ổn định động của vỏ FGM trong môi trường nhiệt độ cao kết hợp tải cơ học (bài toán nhiệt - cơ ghép nối).
- Khảo sát đáp ứng động lực học phi tuyến hình học của kết cấu vỏ FGM khi chịu tải trọng va đập lớn.
Giá trị tham khảo
Luận án là tài liệu tham khảo chuyên sâu và hữu ích cho:
- Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Cơ học vật rắn biến dạng, Cơ kỹ thuật, Kỹ thuật Cơ khí, Kỹ thuật Hàng không và Xây dựng công trình biển; đặc biệt là chương 2 và chương 4 với hệ thống phương trình vi phân trạng thái $10 \times 10$ và thuật toán ghép nối ma trận độ cứng động lực dạng chữ T.
- Kỹ sư tính toán và chuyên gia thiết kế: Làm việc trong lĩnh vực tính toán kết cấu vỏ mỏng, bình chịu áp lực, thân tên lửa, vỏ động cơ và đường ống dẫn ngầm chế tạo từ vật liệu composite/FGM đặt trên nền đất đá đàn hồi.
- Các nhà phát triển phần mềm cơ học: Tham khảo thuật toán PTLT kết hợp phương pháp đường cong đáp ứng trên môi trường MATLAB để tích hợp mô đun tính toán động lực học kết cấu vỏ tròn xoay với số bậc tự do tối thiểu.
Câu hỏi thường gặp
1. Phương pháp Phần tử liên tục (PTLT) có ưu điểm gì nổi bật so với phương pháp Phần tử hữu hạn (PTHH) truyền thống?
Trong PTHH, trường chuyển vị bên trong phần tử là hàm xấp xỉ tĩnh (đa thức), do đó để thu được tần số dao động chính xác ở miền tần số trung bình và cao thì phải chia lưới rất mịn, dẫn đến số bậc tự do lớn và tăng sai số tích lũy. Ngược lại, phương pháp PTLT (ma trận độ cứng động lực) sử dụng nghiệm chính xác của phương trình vi phân chuyển động động lực học làm hàm dạng. Do đó, mỗi đoạn kết cấu đồng nhất chỉ cần đúng một phần tử liên tục duy nhất (số lượng phần tử là tối thiểu), giải quyết triệt để bài toán lưu trữ, tăng tốc độ tính toán và đảm bảo độ chính xác cao ở cả miền tần số thấp, trung bình và cao.
2. Mô hình nền đàn hồi Pasternak thể hiện sự khác biệt gì so với mô hình Winkler trong bài toán vỏ FGM?
Mô hình Winkler chỉ coi nền như các lò xo tuyến tính độc lập song song, tại đó chuyển vị chỉ xuất hiện cục bộ đúng tại vị trí có lực tác dụng và không có sự liên kết trượt giữa các lò xo lân cận ($p = k_w w$). Mô hình Pasternak bổ sung thêm một lớp trượt cắt với hệ số $k_p$ ($p = k_w w - k_p \nabla^2 w$), mô tả chính xác hơn ứng xử biến dạng liên tục của môi trường nền thực tế xung quanh vỏ.
3. Thuật toán ghép nối dạng chữ T được luận án xây dựng nhằm mục đích gì?
Thuật toán ghép nối dạng chữ T (song song) được xây dựng trong Chương 4 để giải quyết bài toán kết cấu vỏ trụ FGM có gân gia cường dạng vành. Thay vì phải làm mịn/san đều gân vào chiều dày vỏ như các phương pháp giải tích gần đúng, thuật toán chữ T cho phép ghép nối trực tiếp ma trận độ cứng động lực của phần tử vành tròn đàn hồi vào vị trí nút tiếp giáp giữa hai phần tử vỏ trụ liên tục, phản ánh đúng tương tác cơ học thực tế của gân gia cường.
4. Quy luật biến đổi tỉ lệ thể tích $V_c$ ảnh hưởng thế nào đến tần số dao động riêng của kết cấu?
Tỉ lệ thể tích pha gốm $V_c$ quyết định tỉ lệ giữa mô đun đàn hồi $E(z)$ và khối lượng riêng $\rho(z)$. Khi số mũ phân bố $p$ tăng lên, hàm lượng pha gốm ($E_c$ cao) giảm dần và hàm lượng kim loại ($E_m$ thấp hơn) tăng lên, làm giảm độ cứng tương đương của toàn bộ kết cấu vỏ FGM. Kết quả là các tần số dao động riêng của hệ vỏ đều có xu hướng giảm khi $p$ tăng.
5. Phương pháp Đường cong đáp ứng xác định tần số riêng như thế nào trong môi trường MATLAB?
Sau khi thiết lập ma trận độ cứng động lực tổng thể $[K(\omega)]$, một lực kích động điều hòa đơn vị $F = 1 \cdot e^{i\omega t}$ được đặt vào một điểm biên trên kết cấu. Chương trình quét qua một dải tần số $\omega$ với bước tăng nhỏ và giải hệ phương trình $U(\omega) = [K(\omega)]^{-1} F$. Khi tần số kích động $\omega$ tiệm cận tần số dao động tự do của kết cấu, hiện tượng cộng hưởng xuất hiện khiến biên độ chuyển vị $w$ tăng vọt tạo thành các đỉnh nhọn trên đường cong đáp ứng. Tọa độ hoành độ của các đỉnh này chính là các tần số dao động riêng cần tìm.
Kết luận
Luận án của tác giả Lê Quang Vinh đã giải quyết trọn vẹn bài toán dao động tự do của các kết cấu vỏ tròn xoay đơn lẻ, vỏ dạng bậc và vỏ liên hợp phức tạp bằng vật liệu FGM có tương tác với nền đàn hồi Winkler-Pasternak bằng phương pháp Phần tử liên tục. Tác giả đã xây dựng thành công hệ phương trình vi phân trạng thái cấp 10, phát triển các thuật toán ghép nối ma trận độ cứng động lực nối tiếp và dạng chữ T, đồng thời thiết lập chương trình tính toán số hiệu quả trên phần mềm MATLAB dựa trên phương pháp đường cong đáp ứng. Các kết quả và quy luật cơ học được luận án công bố có độ tin cậy cao, bổ sung nguồn dữ liệu chuẩn xác cho thư viện phần tử liên tục và đóng góp cơ sở khoa học quan trọng cho việc tính toán, thiết kế kết cấu FGM trong thực tiễn kỹ thuật.
