Tổng quan về luận án

Nghiên cứu về cơ học vật liệu và tính toán ổn định phi tuyến kết cấu vỏ mỏng đóng vai trò trọng yếu trong kỹ thuật hàng không vũ trụ, công trình ngầm, công nghệ hạt nhân và chế tạo bồn bể chịu áp lực cao. Luận án tiến sĩ chuyên ngành Cơ kỹ thuật (Mã số: 9 52 01 01) với đề tài "Phân tích ổn định tĩnh của vỏ trụ và vỏ trống làm từ FGM và FG-CNTRC có kể đến tính đàn hồi của liên kết biên" do nghiên cứu sinh Phạm Thanh Hiếu thực hiện dưới sự hướng dẫn khoa học của PGS.TS. Hoàng Văn Tùng và PGS.TS. Đào Như Mai 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 đã giải quyết thấu đáo khoảng trống nghiên cứu then chốt: mô hình hóa giải tích ứng xử vồng (buckling) và sau vồng (postbuckling) của các cấu trúc vỏ tròn xoay kín làm từ vật liệu composite tiên tiến chịu tải trọng phức hợp kết hợp điều kiện biên đàn hồi thực tế.

Khoảng trống học thuật (research gap) cốt lõi xuất phát từ thực trạng các công bố quốc tế kinh điển (như Shen [15, 33–38], Brush và Almroth [22]) chủ yếu quy ước điều kiện biên lý tưởng hóa: tự do dịch chuyển tiếp tuyến hoàn toàn (movable) khi chịu tải cơ học hoặc bị chặn cứng hoàn toàn (immovable) khi chịu tải nhiệt. Thực tế kỹ thuật đòi hỏi phải mô hình hóa liên kết đàn hồi biên bán di động (partially movable elastic constraints). Hơn nữa, các phân tích trước đây đa phần tập trung vào kết cấu hở (tấm, panel phẳng hoặc cong) mà thiếu vắng các giải pháp giải tích tường minh cho kết cấu vỏ tròn xoay khép kín theo phương vĩ tuyến (closed circumference) như vỏ trụ tròn (cylindrical shells) và vỏ trống thoải (barrel shells) làm từ vật liệu cơ tính biến đổi FGM (Functionally Graded Materials) và composite gia cường ống nano carbon cơ tính biến đổi FG-CNTRC (Functionally Graded Carbon Nanotube-Reinforced Composites).

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được xác lập chặt chẽ:

  1. Câu hỏi nghiên cứu 1 (RQ1): Mức độ ràng buộc đàn hồi tiếp tuyến tại biên ảnh hưởng như thế nào đến tải tới hạn rẽ nhánh và đường cân bằng sau vồng của vỏ trụ sandwich FGM cấu hình loại A (lõi thuần nhất, mặt FGM) và loại B (lõi FGM, mặt thuần nhất)?
  2. Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố nano carbon (UD, FG-V, FG-A, FG-O, FG-X) tương tác ra sao với độ cong Gauss, nền đàn hồi hai tham số Pasternak và sự phụ thuộc nhiệt độ ($T\text{-D}$) trong việc chi phối hiện tượng hóp (snap-through) của vỏ trống FG-CNTRC?
  3. Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để khắc phục hạn chế toán học của dạng nghiệm độ võng ba số hạng khi mở rộng từ Lý thuyết vỏ cổ điển (CST) sang Lý thuyết biến dạng trượt bậc nhất (FSDT) đối với vỏ trụ và vỏ trống tương đối dày?
  • Giả thuyết khoa học 1 (H1): Tăng độ cứng ràng buộc đàn hồi tiếp tuyến sẽ làm tăng phi tuyến độ cứng tổng thể của vỏ, nâng cao tải tới hạn nhưng đồng thời làm tăng độ nhạy và cường độ hóp (snap-through intensity) trong giai đoạn sau vồng nhiệt đàn hồi.
  • Giả thuyết khoa học 2 (H2): Cấu hình phân bố sợi nano FG-X tối ưu hóa mô men quán tính uốn tại các lớp biên mặt ngoài, mang lại khả năng kháng vồng vượt trội nhất so với phân bố đều UD và các dạng phân bố FG khác dưới tác dụng của trường tải trọng cơ - nhiệt kết hợp.

Khung lý thuyết tổng quát tích hợp Lý thuyết vỏ mỏng Donnell với quan hệ phi tuyến hình học von Kármán – Donnell, Lý thuyết biến dạng trượt bậc nhất FSDT, Quy tắc hỗn hợp suy rộng hiệu chỉnh bởi Động lực học phân tử (MD) và Tiêu chuẩn tĩnh về ổn định cân bằng lân cận. Nghiên cứu mang lại đột phá phương pháp luận khi thiết lập thành công nghiệm giải tích dạng tường minh, kiểm chứng định lượng trên dải thông số nhiệt độ từ 300 K đến 600 K, áp lực ngoài từ dải kPa đến MPa, lực nén dọc trục hàng trăm kN, tạo lập cơ sở dữ liệu benchmark cho thiết kế kết cấu vỏ thế hệ mới.

Literature Review và Positioning

Bối cảnh y văn quốc tế ghi nhận sự tiến hóa của hai dòng vật liệu nanocomposite và functionally graded materials:

                                                   [HIỆN DIỆN NGHIÊN CỨU: Luận án NCS Phạm Thanh Hiếu]
                                                   • Vỏ trụ tròn & Vỏ trống FGM / FG-CNTRC
                                                   • Liên kết biên đàn hồi tiếp tuyến (Partially movable)
                                                   • Tải trọng kết hợp Cơ - Nhiệt & Nền Pasternak
                                                   • Nghiệm giải tích 3 số hạng (CST) & 2 số hạng (FSDT)

Nghiên cứu FGM bắt đầu từ các công bố tiên phong của Koizumi [1] và Suresh et al. [2], tập trung vào composite kim loại - gốm biến thiên liên tục tỷ lệ thể tích để triệt tiêu tập trung ứng suất tiếp xúc lớp. Đối với vỏ FGM, Shahsiah và Eslami [16], Wu cùng ccs [17] giải bài toán vồng tuyến tính nhiệt theo CST; Khazaeinejad và ccs [18], Sofiyev và ccs [28–30] phát triển phương pháp giải tích dựa trên FSDT; Shen [33–38] ứng dụng lý thuyết biến dạng trượt bậc cao HSDT kết hợp phương pháp tiệm cận giải tích hai bước và thuật toán lặp xác định ứng xử sau vồng vỏ FGM chịu nén dọc trục, áp lực ngoài và trường nhiệt. Về vật liệu gia cường nano carbon, kể từ phát hiện của Iijima [3, 4] về SWCNT và MWCNT với các đặc tính phi thường (mô đun đàn hồi vượt 1 TPa, độ bền kéo thực nghiệm đạt 63 GPa, độ ổn định nhiệt trong chân không lên tới 2800 °C), Shen [15] đã đưa ra bước ngoặt mang tính nền tảng: đề xuất khái niệm FG-CNTRC với tỷ lệ thể tích CNTs biến thiên hàm tuyến tính dọc chiều dày kết cấu tấm.

Cuộc tranh luận học thuật sâu sắc tồn tại giữa hai quan điểm tiếp cận mô hình hóa:

  1. Quan điểm vi cơ học liên tục cổ điển (Classical micromechanics view): Đại diện bởi Macías và ccs [83, 84] sử dụng mô hình Eshelby-Mori-Tanaka, cho rằng tương tác pha nền polymer - hạt độn nano carbon có thể mô phỏng thuần túy qua trường tenxơ ứng suất bao hàm đẳng trị mà không cần tính đến hiệu ứng kích thước nano phụ thuộc phương pháp tổng hợp.
  2. Quan điểm quy tắc hỗn hợp suy rộng kết hợp động lực học phân tử (Extended rule of mixtures with MD parameters): Khởi xướng bởi Shen [15, 119], chứng minh rằng mô hình Eshelby-Mori-Tanaka truyền thống đánh giá quá cao độ cứng của nanocomposite khi không xét đến liên kết liên bề mặt lỏng lẻo, hiện tượng trượt nano và gợn sóng ống. Quy tắc hỗn hợp suy rộng bổ sung các hệ số hiệu quả $\eta_j$ ($j=1,2,3$) được chuẩn hóa trực tiếp từ mô phỏng động lực học phân tử (MD) phản ánh chính xác cơ tính dị hướng thực tế.

Positioning của luận án định vị tại điểm giao thoa chưa từng được giải quyết: mở rộng mô hình FG-CNTRC của Shen sang kết cấu vỏ kín hai chiều cong (vỏ trụ và vỏ trống) đặt trên nền đàn hồi Pasternak hai tham số ($k_1, k_2$), kết hợp xử lý đồng thời tính đàn hồi của liên kết biên tiếp tuyến và đặc tính vật liệu phụ thuộc nhiệt độ ($T\text{-D}$).

So sánh với các nghiên cứu quốc tế điển hình:

  • So với công trình của Shen [121] về vỏ trụ FG-CNTRC chịu nhiệt và Jam & Kiani [104] về vỏ nón FG-CNTRC: Các tác giả chỉ xem xét điều kiện biên chặn cứng hoàn toàn ($u=0$ tại biên), dẫn đến việc đánh giá thiên lệch tải nhiệt tới hạn. Luận án tổng quát hóa với tham số độ cứng biên tiếp tuyến không thứ nguyên $c_t \in [0, \infty)$, bao trùm trạng thái tự do ($c_t=0$), ngàm đàn hồi ($0 < c_t < \infty$) và chặn cứng tuyệt đối ($c_t \to \infty$).
  • So với các phân tích của Sofiyev et al. [61, 62, 110] và Ansari et al. [96, 106]: Vốn chủ yếu dừng lại ở bài toán vồng tuyến tính đơn số hạng hoặc sử dụng phương pháp số lưới phần tử hữu hạn/GDQM tốn kém chi phí tính toán, luận án xây dựng hệ nghiệm giải tích đa số hạng tường minh, thiết lập chính xác điểm phân nhánh rẽ nhánh (bifurcation point) và quỹ đạo cân bằng phi tuyến sau vồng.

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

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

Luận án đóng góp mở rộng nền tảng cơ học phi tuyến kết cấu composite qua 4 trụ cột lý thuyết:

  1. Mở rộng lý thuyết ổn định vỏ mỏng Donnell và quan hệ biến dạng - chuyển vị von Kármán – Donnell: Thiết lập hệ phương trình vi phân phi tuyến tương thích và cân bằng biểu diễn hoàn toàn qua hàm ứng suất Airy ($F$) và hàm độ võng ($W$), tích hợp biến dạng màng trước vồng (prebuckling deformation $W_0$), biến dạng vồng tuyến tính ($W_1$) và số hạng biến dạng phi tuyến đối xứng trục ($W_2$): $$W(x, y) = W_0 + W_1 \sin(\alpha x)\sin(\beta y) + W_2 \sin^2(\alpha x)$$ trong đó $\alpha = m\pi/L$, $\beta = n/R$, $m$ là số nửa sóng dọc trục, $n$ là số sóng theo chu vi.
  2. Mở rộng trường ứng xử nhiệt đàn hồi sang cấu trúc vỏ trống hai độ cong (Barrel Shells): Khác với vỏ trụ ($R_x = \infty$), vỏ trống sở hữu độ cong kinh tuyến $R_x$ và độ cong vĩ tuyến $R_\theta$ đồng thời, tạo ra độ cong Gauss dương. Luận án mở rộng phương trình Donnell cho vỏ trống thoải, chứng minh định lượng cơ chế tương tác giữa độ cong Gauss và ứng suất nén trước vồng nhiệt.
  3. Đặc tả mô hình vật liệu Sandwich FGM phân lớp: Xác lập giải tích biểu thức cơ tính biến thiên phi đối xứng (Loại B: Ceramic/FGM/SUS304) và đối xứng (Loại A: FGM/SUS304/FGM) theo quy luật hàm lũy thừa: $$V_m(z) = \left(\frac{2z + h - 2h_f}{2h - 4h_f}\right)^k, \quad V_c(z) = 1 - V_m(z)$$ chỉ rõ sự chuyển dịch của mặt trung hòa hình học so với mặt trung hòa vật lý gây ra hiệu ứng ghép uốn - kéo (bending-stretching coupling).
  4. Chuẩn hóa mô hình nanocomposite dị hướng FG-CNTRC: Cụ thể hóa 4 kiểu phân bố hàm gradient theo phương chiều dày $h$:
  • Phân bố đều (UD): $V_{cnt}(z) = V_{cnt}^*$
  • Phân bố chữ V (FG-V): $V_{cnt}(z) = \left(1 + \frac{2z}{h}\right)V_{cnt}^*$
  • Phân bố chữ O (FG-O): $V_{cnt}(z) = 2\left(1 - \frac{2|z|}{h}\right)V_{cnt}^*$
  • Phân bố chữ X (FG-X): $V_{cnt}(z) = 4\left(\frac{|z|}{h}\right)V_{cnt}^$ với $V_{cnt}^ = \frac{w_{cnt}}{w_{cnt} + (\rho_{cnt}/\rho_m) - (\rho_{cnt}/\rho_m)w_{cnt}}$, trong đó $w_{cnt}$ là tỷ lệ phần trăm khối lượng nano carbon, $\rho_{cnt}$ và $\rho_m$ là khối lượng riêng tương ứng của hạt độn CNT và vật liệu nền polymer (PMMA hoặc PmPV).

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

Khung phân tích của luận án là sự hợp nhất của ba lý thuyết cơ học kết cấu: Lý thuyết vỏ cổ điển Donnell (CST), Lý thuyết biến dạng trượt bậc nhất (FSDT), và Mô hình tương tác nền đàn hồi Pasternak hai thông số với áp lực tiếp xúc $q_e = k_1 W - k_2 \nabla^2 W$.

Điểm sáng tạo phương pháp luận đột phá nằm ở việc dẫn xuất phương trình điều kiện biên tiếp tuyến đàn hồi trung bình: $$\frac{c_t E_0 h}{L} \int_0^{2\pi R} \int_0^L \frac{\partial u}{\partial x} dx dy + \int_0^{2\pi R} N_x \Big|_{x=0, L} dy = 0$$ Cho phép liên kết biến thiên lực dọc trục bị động sinh ra trong quá trình biến dạng nhiệt màng trực tiếp với độ cứng tiếp tuyến lò xo biên, khắc phục triệt để giả định phi thực tế của các mô hình biên cứng truyền thống.

Điều kiện biên điều chỉnh:

  • Giả thiết vỏ mỏng và vừa dày hoàn hảo, biến dạng đàn hồi tuyến tính tuân theo định luật Hooke tổng quát.
  • Bỏ qua hiện tượng bong tách lớp (delamination) trong kết cấu sandwich FGM.
  • Ống nano carbon SWCNT phân bố định hướng lý tưởng, không bị uốn gấp khúc cục bộ.
  • Bán kính cong kinh tuyến $R_x \gg R_\theta$ đối với vỏ trống (giả thiết vỏ thoải).

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

Thiết kế nghiên cứu

Nghiên cứu vận dụng thế giới quan duy thực phản biện (critical realism) kết hợp trường phái thực chứng định lượng (positivism) trong cơ học vật rắn biến dạng. Thiết kế nghiên cứu đa tầng (multi-level research design) được cấu trúc chặt chẽ qua ba cấp độ mô hình hóa:

  1. Cấp độ vi mô (Micro-scale): Mô hình hóa đặc trưng đàn hồi dị hướng của nanocomposite đơn lớp thông qua hệ số hiệu quả $\eta_1, \eta_2, \eta_3$ phụ thuộc nhiệt độ, xác định mô đun đàn hồi dọc $E_{11}$, mô đun đàn hồi ngang $E_{22}$, và mô đun trượt $G_{12}$.
  2. Cấp độ trung mô (Meso-scale): Tích phân trường cơ tính biến thiên liên tục qua chiều dày kết cấu $h$ nhằm thiết lập ma trận độ cứng màng ($A_{ij}$), độ cứng ghép ($B_{ij}$), độ cứng uốn ($D_{ij}$) và các ma trận hiệu chỉnh trượt ($A_{44}^s, A_{55}^s$).
  3. Cấp độ vĩ mô (Macro-scale): Thiết lập phương trình vi phân đạo hàm riêng phi tuyến bậc 4 chi phối trường ứng suất và chuyển vị của toàn thể kết cấu vỏ trụ tròn và vỏ trống.

Quy trình nghiên cứu rigorous

Quy trình giải tích bán số (semi-analytical protocol) trải qua 5 bước nghiêm ngặt:

[BƯỚC 1: XÁC LẬP HỆ PHƯƠNG TRÌNH ĐIỀU KHIỂN]
  - CST: Hệ 2 phương trình vi phân phi tuyến cho F và W (von Kármán - Donnell).
  - FSDT: Hệ 5 phương trình đạo hàm riêng chứa W, F, và góc xoay psi_x, psi_theta.
[BƯỚC 2: THIẾT LẬP DẠNG NGHIỆM THỎA MÃN ĐIỀU KIỆN BIÊN]
  - CST: W(x,y) = W0 + W1.sin(alpha.x)sin(beta.y) + W2.sin^2(alpha.x)
  - FSDT: Dạng nghiệm 2 số hạng khử triệt để kỳ dị biên cho vỏ kín.
  - Giải giải tích tường minh phương trình tương thích để tìm hàm ứng suất Airy F.
[BƯỚC 3: CHIẾU PHI TUYẾN PHƯƠNG PHÁP GALERKIN]
  - Tích phân trực giao: int_0^{2pi R} int_0^L [Phương trình cân bằng] . delta W dx dy = 0
  - Đại số hóa hệ phương trình vi phân thành phương trình đại số phi tuyến bậc 3.
[BƯỚC 4: THUẬT TOÁN LẶP NHIỆT ĐÀN HỒI (CHO VẬT LIỆU T-D)]
  - Khởi tạo nhiệt độ T_0 = 300 K -> Tính toán cơ tính E(T), alpha(T).
  - Tính tải nhiệt tới hạn T_cr^(i) -> Cập nhật lại cơ tính tại T_cr^(i) -> Lặp hội tụ (|T^(i+1) - T^(i)| < 10^-4).
[BƯỚC 5: XÁC MINH (VERIFICATION) & PHÂN TÍCH ĐỘ NHẠY ĐA THAM SỐ]
  - Đối chuẩn dữ liệu với nghiệm của Shen, Sofiyev, Brush & Almroth.
  - Khảo sát ảnh hưởng của c_t, k_1, k_2, V*_cnt, tỷ số R/h, L/R.

Độ tin cậy (reliability) và tính hợp thức (validity) được đảm bảo thông qua kiểm định chéo (triangulation):

  • Tính hợp thức cấu trúc (Construct validity): Mô hình hóa đầy đủ phương trình liên tục và tenxơ biến dạng phi tuyến.
  • Tính hợp thức nội tại (Internal validity): Khử hoàn toàn sai số làm tròn số học bằng cách dẫn xuất biểu thức giải tích tường minh trước khi đưa vào tính toán số.
  • Tính hợp thức bên ngoài (External validity): So sánh đối chuẩn kết quả với các công bố uy tín quốc tế trong các trường hợp suy biến biên kinh điển.

Data và phân tích

Bộ thông số cơ lý và hình học kiểm định:

  • Vật liệu FGM: Kết hợp Kim loại (Thép không gỉ SUS304: $E_m(T) = 201.04 \times 10^9 - 0.0543T\text{ Pa}$, $\alpha_m(T) = 12.33 \times 10^{-6} + 8.08 \times 10^{-9}T\text{ K}^{-1}$) và Gốm (Silicon Nitride $\text{Si}_3\text{N}_4$: $E_c(T) = 348.43 \times 10^9 - 0.0307T\text{ Pa}$, $\alpha_c(T) = 5.87 \times 10^{-6} + 9.09 \times 10^{-10}T\text{ K}^{-1}$).
  • Vật liệu FG-CNTRC: Pha nền polymer PMMA ($E_m = (3.52 - 0.0034T)\text{ GPa}$, $\nu_m = 0.34$, $\alpha_m = 45(1 + 0.0005\Delta T)\times 10^{-6}\text{ K}^{-1}$) và PmPV. Pha gia cường SWCNT (10,10) ($r = 0.678\text{ nm}$, $h = 0.067\text{ nm}$, $L = 100\text{ nm}$, $E_{11}^{cnt} = 5.6466\text{ TPa}$, $E_{22}^{cnt} = 7.0805\text{ TPa}$, $G_{12}^{cnt} = 1.9445\text{ TPa}$, $\alpha_{11}^{cnt} = 3.4584 \times 10^{-6}\text{ K}^{-1}$). Tỷ lệ thể tích nano carbon khảo sát: $V_{cnt}^* = 0.11, 0.14, 0.17, 0.28$.
  • Phần mềm và môi trường tính toán: Toàn bộ biểu thức giải tích và thuật toán lặp hội tụ được lập trình trên nền tảng tính toán ký hiệu Wolfram Mathematica và xử lý dữ liệu ma trận bằng MATLAB.

Kiểm chứng độ hội tụ và đối chuẩn kết quả (Robustness checks): Kết quả đối chuẩn áp lực ngoài tới hạn $P_{cr}\text{ (MPa)}$ của vỏ trụ FGM ($\text{Si}_3\text{N}_4/\text{SUS304}$) tại $T = 300\text{ K}$ cho thấy sự trùng khớp tuyệt đối với nghiên cứu của Shen [34]:

  • Với $R/h = 100, L/R = 2, k = 0$ (thuần gốm): Luận án đạt $P_{cr} = 3.421\text{ MPa}$ (Shen [34] công bố $3.420\text{ MPa}$, sai số $< 0.03%$).
  • Với $k = 1$: Luận án đạt $P_{cr} = 2.654\text{ MPa}$ (Shen [34] công bố $2.652\text{ MPa}$, sai số $< 0.07%$).
  • Với $k \to \infty$ (thuần kim loại): Luận án đạt $P_{cr} = 1.968\text{ MPa}$ (Shen [34] công bố $1.967\text{ MPa}$).

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

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

Năm phát hiện then chốt được minh chứng bằng dữ liệu định lượng:

  1. Cơ chế phân bố nano carbon FG-X mang lại khả năng kháng vồng tối ưu: Bằng chứng dữ liệu từ Chương 3 và 4 chỉ rõ: Vỏ trụ FG-CNTRC với phân bố FG-X luôn đạt tải nén dọc trục và áp lực ngoài tới hạn cao nhất trong mọi dải nhiệt độ. Khi tỷ lệ thể tích $V_{cnt}^*$ tăng từ $0.11$ lên $0.28$, tải nén tới hạn tăng tuyến tính trên $120%$. Ngược lại, phân bố FG-O cho khả năng chịu tải thấp nhất do tập trung sợi ở mặt trung hòa – vùng có ứng suất uốn cực tiểu.
  2. Tác động ức chế hai mặt của độ cứng liên kết biên tiếp tuyến ($c_t$): Đối với tải cơ học, tăng $c_t$ giúp nâng cao khả năng mang tải của vỏ. Tuy nhiên, đối với bài toán tải nhiệt, độ cứng liên kết biên càng lớn ($c_t \to \infty$ - biên tựa cố định) làm sinh ra lực nén trước vồng $N_x^{thermal} = \int E(T)\alpha(T)\Delta T dz$ cực lớn trong mặt phẳng, dẫn đến tải nhiệt tới hạn $\Delta T_{cr}$ tụt giảm nghiêm trọng (từ trạng thái không bị vồng nhiệt khi $c_t=0$ xuống còn $\Delta T_{cr} \approx 65 - 180\text{ K}$ khi bị chặn cứng).
  3. Hiện tượng hóp tức thì (Immediate snap-through) trong ứng xử nhiệt phi tuyến: Khác với ứng xử vồng rẽ nhánh êm thuận của tấm phẳng, vỏ trụ và vỏ trống mỏng chịu nhiệt độ tăng đều dưới liên kết biên đàn hồi xuất hiện quỹ đạo cân bằng sau vồng có độ dốc âm sâu ($dP/dW_{max} < 0$). Sau khi đạt điểm rẽ nhánh tới hạn, kết cấu mất ổn định đột ngột và nảy sang trạng thái biến dạng lớn với cường độ hóp tỷ lệ thuận với tham số $c_t$ và tỷ số $R/h$.
  4. Hiệu ứng hình học vượt trội của vỏ trống (Barrel Shell) so với vỏ trụ tròn: Nhờ có bán kính cong kinh tuyến $R_x$, vỏ trống thoải tích lũy năng lượng biến dạng màng hiệu quả hơn, nâng tải tới hạn nén kết hợp áp lực lên $18 - 32%$ so với vỏ trụ có cùng đường kính danh nghĩa $R_0$ và chiều dài $L$.
  5. Vai trò của tham số cắt nền Pasternak ($k_2$): Lớp cắt trượt Pasternak $k_2$ có tác dụng liên kết các lò xo Winkler độc lập, phân bố đều độ võng cục bộ và làm tăng đáng kể độ cứng sau vồng, triệt tiêu gần như hoàn toàn hiện tượng hóp nguy hiểm ở miền tải trọng cao.

Implications đa chiều

  • Về mặt lý thuyết (Theoretical advances): Đóng góp trực tiếp vào lý thuyết vỏ phi tuyến của Donnell và lý thuyết biến dạng trượt FSDT qua việc tích hợp thành công điều kiện biên đàn hồi bán di động tổng quát, làm sáng tỏ bản chất cơ học của sự tương tác ghép uốn - nén màng trong nanocomposite có cấu trúc dị hướng.
  • Về mặt phương pháp luận (Methodological innovations): Thiết lập quy trình giải tích bán số đóng (closed-form semi-analytical approach) sử dụng dạng nghiệm hai số hạng cho FSDT vỏ kín, cho phép tính toán tức thời (real-time prediction) với độ chính xác tương đương mô phỏng phần tử hữu hạn phi tuyến phi mô men (non-linear FEM) nhưng giảm thiểu $99%$ thời gian xử lý vi tính.
  • Về ứng dụng thực tiễn kỹ thuật (Practical engineering applications): Cung cấp các biểu đồ thiết kế (design charts) và công thức giải tích chuẩn xác để kỹ sư kết cấu lựa chọn cấu hình vật liệu tối ưu: kiến nghị sử dụng phân bố FG-X cho vỏ chịu áp lực ngoài và cấu hình sandwich Loại A (lõi kim loại dẻo, vỏ FGM gốm ngoài) cho các bồn chứa nhiên liệu tên lửa, vỏ tàu lặn biển sâu và ống dẫn nhiệt lò phản ứng.
  • Về chính sách và tiêu chuẩn công nghiệp (Policy and standard pathways): Cung cấp cơ sở khoa học để xây dựng các tiêu chuẩn kỹ thuật thiết kế kết cấu vỏ mỏng vật liệu mới trong ngành công nghiệp quốc phòng và hàng không vũ trụ tại Việt Nam và khu vực.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án chỉ rõ các giới hạn biên nội tại:

  1. Giả thiết hình học hoàn hảo (Perfect geometric configuration): Mô hình giải tích chưa tích hợp trường độ lệch hình học ngẫu nhiên ban đầu (initial geometric imperfections), vốn có thể làm giảm tải tới hạn thực nghiệm của vỏ mỏng xuống $20-40%$.
  2. Mô hình tiếp xúc lý tưởng giữa vỏ và nền: Giả thiết bề mặt vỏ tiếp xúc liên tục $100%$ với nền đàn hồi Pasternak (no-separation assumption), chưa xét đến hiện tượng tách lớp cục bộ khi biên độ võng sau vồng vượt quá giới hạn đàn hồi của nền.
  3. Trạng thái biến dạng đàn hồi thuần túy: Bỏ qua ứng xử dẻo, từ biến nhiệt (creep) của pha nền polymer và sự phá hủy liên kết liên bề mặt giữa ống nano carbon và chất nền ở nhiệt độ vượt quá 600 K.

Chương trình nghiên cứu 10 năm tiếp theo được định hình:

  • Hướng 1: Phát triển mô hình giải tích ngẫu nhiên (stochastic analytical models) đánh giá ảnh hưởng của khuyết tật hình học và phân bố ngẫu nhiên của nano carbon.
  • Hướng 2: Mở rộng nghiên cứu sang bài toán ổn định động học phi tuyến (dynamic buckling) dưới tác dụng của sóng xung kích nổ, tải nhiệt biến thiên theo thời gian sử dụng tiêu chuẩn Budiansky-Roth.
  • Hướng 3: Khảo sát kết cấu vỏ sandwich FGM và FG-CNTRC có gân gia cường lệch tâm (stiffeners) xoắn ốc hoặc gân đẳng hướng (isogrid).
  • Hướng 4: Nghiên cứu thực nghiệm (experimental validation) chế tạo mẫu vỏ nanocomposite bằng công nghệ in 3D và thử nghiệm nén thủy tĩnh trong buồng nhiệt.

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

                                  KỸ THUẬT HÀNG KHÔNG VŨ TRỤ
                             (Vỏ tên lửa, thân máy bay siêu âm)
(Vỏ tàu ngầm, ống dẫn đáy biển)      CỦA NGHIÊN CỨU         (Bồn chịu áp, vỏ lò phản ứng)
                                 HỆ VI CƠ ĐIỆN TỬ (MEMS/NEMS)
                               (Cảm biến nano, màng vi cơ đàn hồi)

Nghiên cứu tạo lập tầm ảnh hưởng đa ngành sâu rộng:

  • Tác động học thuật (Academic impact): Mở ra hướng tiếp cận giải tích chuẩn tắc cho bài toán vỏ kín phi tuyến đàn hồi biên. Với chất lượng khoa học cao, các kết quả của luận án đã được công bố trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (Q1/Q2) trong ngành Cơ học kết cấu (như Thin-Walled Structures, Composite Structures, Aerospace Science and Technology), dự kiến thu hút chỉ số trích dẫn cao từ cộng đồng cơ học tính toán quốc tế.
  • Chuyển đổi công nghiệp (Industry transformation): Trực tiếp hỗ trợ các ngành công nghiệp mũi nhọn:
    1. Hàng không vũ trụ: Tối ưu hóa khối lượng cấu trúc thân tên lửa đẩy và vệ tinh thông qua vật liệu siêu nhẹ FG-CNTRC.
    2. Công nghiệp dầu khí và quốc phòng: Thiết kế bồn chứa khí hóa lỏng áp lực cao, thân tàu ngầm mini chịu áp lực nước biển sâu.
    3. Năng lượng hạt nhân: Gia cường ống trao đổi nhiệt làm việc trong môi trường bức xạ và gradient nhiệt độ khắc nghiệt.
  • Lợi ích xã hội và môi trường: Giảm trọng lượng phương tiện vận tải giúp tiết kiệm nhiên liệu hóa thạch, cắt giảm phát thải khí nhà kính $\text{CO}_2$, thúc đẩy phát triển vật liệu xanh và bền vững.

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

  • Nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật, Cơ học vật rắn: Tiếp cận phương pháp luận giải tích mẫu mực, nắm bắt kỹ thuật xử lý phi tuyến hình học von Kármán và phương pháp Galerkin giải tích đa số hạng.
  • Các nhà khoa học và chuyên gia Cơ học tính toán: Sử dụng hệ thống công thức giải tích đóng tường minh và số liệu bảng biểu làm bộ dữ liệu chuẩn (benchmarks) để kiểm định độ chính xác của các thuật toán số mới (FEM, XFEM, IGA, Meshless).
  • Kỹ sư R&D trong ngành công nghiệp vật liệu mới: Khai thác các quy luật phân bố nano carbon tối ưu (FG-X) để ứng dụng trực tiếp vào quy trình đúc, kéo sợi và chế tạo sản phẩm composite tiên tiến.
  • Cơ quan quản lý và thiết lập quy chuẩn kỹ thuật: Tham khảo cơ sở dữ liệu khoa học để xây dựng tiêu chuẩn an toàn kết cấu mỏng làm việc trong môi trường chịu nhiệt và chịu áp.

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 mở rộng thành công Lý thuyết vỏ Donnell phi tuyến và Lý thuyết biến dạng trượt bậc nhất (FSDT) cho các kết cấu vỏ tròn xoay khép kín (vỏ trụ và vỏ trống) thông qua việc tích hợp định lượng điều kiện biên đàn hồi tiếp tuyến bán di động (partially movable tangential boundary condition). Luận án đã thiết lập hệ nghiệm giải tích đa số hạng biểu diễn đồng thời ba trạng thái: biến dạng màng trước vồng ($W_0$), dạng vồng rẽ nhánh cơ bản ($W_1$), và dạng võng đối xứng trục phi tuyến sau vồng ($W_2$). Đóng góp này khắc phục hoàn toàn sự chia cắt phi thực tế giữa hai giả định biên kinh điển (chặn cứng immovable và thả tự do movable) tồn tại suốt nhiều thập kỷ trong các công trình của Shen [33–38] hay Brush & Almroth [22].

2. Điểm mới về phương pháp luận nghiên cứu so với các công bố quốc tế tiền nhiệm?

So sánh với nghiên cứu của Shen [15, 121] sử dụng phương pháp tiệm cận hai bước (two-step perturbation method) vốn đòi hỏi khai triển chuỗi theo tham số bé và phụ thuộc nặng nề vào điều kiện biên đơn giản, phương pháp của luận án kết hợp biến đổi hàm ứng suất Airy với phương pháp chiếu Galerkin toàn cục tích phân trên toàn bộ miền vỏ kín. Đặc biệt, đối với vỏ tương đối dày theo FSDT, luận án đề xuất dạng nghiệm hai số hạng đột phá, giải quyết triệt để vấn đề xung đột bậc tự do và kỳ dị biên mà dạng nghiệm ba số hạng gặp phải khi mô hình hóa vỏ kín có xét đến biến dạng trượt ngang ($A_{44}^s, A_{55}^s$).

3. Phát hiện bất ngờ nhất (counter-intuitive) có sự ủng hộ của dữ liệu là gì?

Phát hiện bất ngờ nhất nằm ở tác động nghịch đảo của độ cứng liên kết biên tiếp tuyến đối với bài toán ổn định nhiệt: Trong khi việc tăng độ cứng lò xo biên ($c_t$) giúp tăng khả năng chịu tải cơ học (áp lực ngoài và nén dọc trục), thì đối với tải nhiệt độ, tăng $c_t$ lại làm giảm mạnh tải nhiệt tới hạn $\Delta T_{cr}$ và kích hoạt hiện tượng hóp tức thì (immediate snap-through) với mức độ giải phóng năng lượng uốn cực kỳ dữ dội. Dữ liệu tính toán chứng minh rằng khi $c_t$ tăng từ 0 lên vô cùng, vỏ trụ FGM chuyển từ trạng thái hoàn toàn không bị mất ổn định nhiệt sang trạng thái bị vồng sập ở nhiệt độ chỉ vài trăm Kelvin.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp đầy đủ:

  • Tường minh toàn bộ hệ số ma trận độ cứng vật liệu ($A_{ij}, B_{ij}, D_{ij}, A_{44}^s, A_{55}^s$).
  • Công thức giải tích chi tiết của các tích phân Galerkin và hệ phương trình đại số phi tuyến.
  • Bảng thông số phụ thuộc nhiệt độ $T\text{-D}$ của vật liệu nền polymer, thép không gỉ SUS304, gốm $\text{Si}_3\text{N}_4$, và ống nano carbon SWCNT (10,10).
  • Thuật toán lặp sai phân nhiệt hồi quy được chuẩn hóa tường minh từng bước logic.

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình 10 năm tập trung chuyển dịch từ mô hình hóa giải tích đơn lẻ sang hệ thống kỹ thuật phức hợp:

  • Giai đoạn 1 (Năm 1–3): Mở rộng giải tích ổn định động lực học phi tuyến vỏ FG-CNTRC chịu sóng nổ và va chạm thủy động lực học.
  • Giai đoạn 2 (Năm 4–6): Tích hợp tính không hoàn hảo hình học ngẫu nhiên (random geometric imperfections) và mô hình tổn thương liên kết nano (interfacial debonding).
  • Giai đoạn 3 (Năm 7–10): Thực nghiệm chế tạo vỏ trống nanocomposite bằng công nghệ in 3D sợi liên tục và xây dựng bộ tiêu chuẩn thiết kế quốc gia cho kết cấu composite chịu áp lực cao.

Kết luận

Luận án tiến sĩ của NCS Phạm Thanh Hiếu là một công trình khoa học công phu, nghiêm túc, đạt chuẩn mực học thuật quốc tế xuất sắc trong lĩnh vực Cơ học vật rắn và Cơ kỹ thuật.

Năm đóng góp cốt lõi kết tinh giá trị của công trình:

  1. Tiên phong thiết lập thành công mô hình giải tích ổn định phi tuyến cho vỏ trụ tròn sandwich FGM (loại A và B) và vỏ trống FG-CNTRC có kể đến tính đàn hồi tiếp tuyến của liên kết biên.
  2. Làm sáng tỏ cơ chế cơ học lượng hóa ưu thế vượt trội của kiểu phân bố sợi FG-X, giúp tăng khả năng chịu tải nén và áp lực ngoài lên đến $68%$ so với phân bố đều.
  3. Khám phá và giải thích bản chất hiện tượng hóp (snap-through) dưới tác động đồng thời của trường nhiệt độ phân bố đối xứng trục và liên kết biên bán di động.
  4. Đột phá phương pháp luận với dạng nghiệm hai số hạng cho lý thuyết FSDT vỏ kín, mở ra giải pháp giải tích chuẩn xác cho kết cấu vỏ composite dày trung bình.
  5. Xây dựng hệ thống công thức giải tích đóng và dữ liệu số benchmark tin cậy, phục vụ trực tiếp cho công tác nghiên cứu khoa học và thiết kế kỹ thuật ứng dụng trong các ngành công nghiệp công nghệ cao.