Tổng quan về luận án

Nghiên cứu của tác giả Vũ Thành Long (2022) tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam, 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ập trung giải quyết bài toán: "Ổn định của vỏ FGM có lỗ rỗng và panel sandwich FG-CNTRC với các cạnh biên chịu liên kết đàn hồi" thuộc chuyên ngành Cơ kỹ thuật (Mã số: 9 52 01 01).

Trong kỹ thuật kết cấu hiện đại, đặc biệt là hàng không vũ trụ, công trình biển và bình chứa cao áp, vật liệu composite cơ tính biến đổi (Functionally Graded Material - FGM) và composite gia cường ống nano carbon có cơ tính biến đổi (Functionally Graded Carbon Nanotube Reinforced Composite - FG-CNTRC) đóng vai trò then chốt nhờ khả năng chịu nhiệt và tỷ số độ cứng trên trọng lượng vượt trội. Tuy nhiên, quá trình chế tạo thực tế không thể tránh khỏi các khuyết tật vi mô dạng lỗ rỗng (porosities). Đồng thời, điều kiện biên cơ học trong thực tế kỹ thuật hiếm khi là ngàm cứng tuyệt đối hay tựa tự do lý tưởng, mà mang bản chất liên kết đàn hồi bán cứng (elastically restrained edges).

Khoảng trống học thuật cốt lõi (research gap) nằm ở chỗ: các công bố quốc tế kinh điển thường giả thiết cấu trúc vật liệu hoàn hảo ($\alpha = 0$), bỏ qua sự tương tác phức tạp giữa độ rỗng phân bố đều/không đều và độ cứng đàn hồi của liên kết biên lên ứng xử mất ổn định tuyến tính lẫn phi tuyến sau vồng của các dạng vỏ phức tạp như vỏ trống (toroidal shell), vỏ trụ tương đối dày, vỏ cầu thoải (spherical cap) và kết cấu panel sandwich FG-CNTRC đa lớp.

Luận án đặt ra và giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học tương ứng (Hypotheses - H):

  • RQ1 & H1: Ảnh hưởng của độ rỗng vi mô (phân bố đều Even và không đều Uneven) cùng liên kết đàn hồi tiếp tuyến $\chi$ làm thay đổi căn bản tải tới hạn tuyến tính của vỏ trống dày FGM; giả thuyết rằng độ rỗng làm suy giảm mạnh độ cứng chống uốn nhưng có thể bù đắp bằng việc tối ưu hóa độ cong Gauss và nền đàn hồi Pasternak.
  • RQ2 & H2: Ứng xử phi tuyến sau vồng của vỏ trụ tròn FGM tương đối dày có độ võng màng đáng kể trước khi vồng; giả thuyết rằng việc sử dụng nghiệm dạng hai số hạng sẽ khắc phục sai số lớn của nghiệm đơn số hạng truyền thống.
  • RQ3 & H3: Ổn định phi tuyến của chỏm cầu thoải và tấm tròn FGM có lỗ rỗng chịu áp lực ngoài và nhiệt độ cao trên nền đàn hồi phi tuyến ba tham số; giả thuyết rằng hiện tượng hóp (snap-through) chịu chi phối mạnh bởi tỷ số độ cong $a/R$ và tham số nền phi tuyến $k_3$.
  • RQ4 & H4: Hiệu năng mang tải sau vồng của hai mô hình panel sandwich FG-CNTRC (Loại A: lõi đẳng hướng - bề mặt FG-CNTRC; Loại B: lõi FG-CNTRC - bề mặt đẳng hướng); giả thuyết rằng mô hình Loại A tối ưu hóa khả năng chống vồng nhiệt và nén dọc trục hơn mô hình Loại B.

Khung lý thuyết tích hợp bao gồm: Lý thuyết biến dạng trượt bậc cao (HSDT) của Reddy & Liu, Lý thuyết biến dạng trượt bậc nhất (FSDT) của Mindlin-Reissner, Lý thuyết vỏ Donnell-Stein-McElman, Quy tắc hỗn hợp cải tiến mở rộng cho vật liệu có lỗ rỗng, Mô hình vi cơ học FG-CNTRC, và Phương pháp giải tích Galerkin phi tuyến kết hợp thuật toán lặp xử lý tính chất vật liệu phụ thuộc nhiệt độ (T-D). Phạm vi khảo sát trải rộng từ vỏ trống $R/h = 80 - 100$, $L^2/Rh = 120 - 400$, tấm tròn $a/h = 50$, panel sandwich $h_f/h = 0.1 - 0.2$, tỷ lệ thể tích CNT $V_{CNT}^* = 0.11 - 0.28$, dải nhiệt độ $T = 300\text{ K} - 1000\text{ K}$.

Literature Review và Positioning

Bức tranh tổng quan học thuật về cơ học kết cấu tấm vỏ composite phân nhánh thành ba dòng nghiên cứu chính:

  1. Dòng lý thuyết tấm vỏ và ổn định tuyến tính/phi tuyến: Khởi nguồn từ lý thuyết vỏ cổ điển (CST) của Donnell, Brush & Almroth (1975) với tiêu chuẩn cân bằng lân cận (adjacent equilibrium criterion). Shahsiah & Eslami (2004), Wu và cộng sự (2005) nghiên cứu vỏ trụ FGM dựa trên CST. Tuy nhiên, CST bỏ qua hoàn toàn biến dạng trượt ngang, dẫn đến ước tính quá cao tải tới hạn đối với kết cấu dày. Reddy & Liu (1985) phát triển HSDT khắc phục triệt để nhược điểm này bằng cách thỏa mãn điều kiện ứng suất tiếp bằng không tại hai bề mặt tự do mà không cần hệ số hiệu chỉnh trượt $K_s$. Shen (2009, 2014) phát triển phương pháp nghiệm nhiễu tiệm cận phân tích sau vồng của vỏ FGM và FG-CNTRC.
  2. Dòng nghiên cứu vật liệu FGM có lỗ rỗng: Wattanasakulpong & Ungbhakorn (2014) tiên phong mô hình hóa ảnh hưởng của lỗ rỗng lên dao động dầm FGM. Tiếp đó, Cong và cộng sự (2018), Gupta & Talha (2018) khảo sát tấm phẳng FGM rỗng. Tuy nhiên, các cấu trúc vỏ kín FGM (vỏ trụ, vỏ trống) và vỏ cầu thoải có lỗ rỗng hoàn toàn bị bỏ ngỏ.
  3. Dòng nghiên cứu kết cấu FG-CNTRC và liên kết biên đàn hồi: Dựa trên đề xuất của Shen (2009) về các dạng phân bố CNT (UD, FG-O, FG-X, FG-V), hàng loạt nghiên cứu của Liew và cộng sự (2014), Kiani (2016), Sofiyev và cộng sự (2017) đã phân tích tấm và panel đơn lớp. Về điều kiện biên, Librescu và cộng sự (1989), Zhang và cộng sự (2015), Hoàng Văn Tùng và cộng sự (2015-2021) đã chứng minh rằng việc giả thiết cạnh biên tự do dịch chuyển (movable) hoặc cố định tuyệt đối (immovable) tạo ra khoảng cách lớn so với ứng xử thực tế của liên kết bu-lông, đinh tán hay gối tựa đàn hồi.

Tranh luận học thuật cốt lõi diễn ra giữa quan điểm sử dụng nghiệm đơn số hạng giải tích truyền thống (tiết kiệm tính toán nhưng bỏ qua trạng thái màng trước vồng của vỏ kín) và nghiệm đa số hạng/nghiệm số phần tử hữu hạn (chính xác hơn nhưng phức tạp và tốn kém tài nguyên). Luận án định vị chính xác khoảng trống này bằng cách thiết lập phương pháp bán giải tích tối ưu, kết hợp nghiệm độ võng hai số hạng và hàm ứng suất Airy trên miền biên liên kết đàn hồi.

Khi so sánh với hai nghiên cứu quốc tế tiêu biểu:

  • So với công trình của Wattanasakulpong & Ungbhakorn (2014) trên dầm FGM, luận án mở rộng quy tắc hỗn hợp cải tiến sang kết cấu vỏ hai độ cong phức tạp, đồng thời khảo sát đầy đủ hai dạng phân bố lỗ rỗng đều (Even) và không đều (Uneven) tương tác với trường nhiệt độ phi tuyến.
  • So với nghiên cứu của Shen & Xiang (2014) về panel FG-CNTRC trên nền đàn hồi Pasternak, luận án phát triển thêm mô hình kết cấu sandwich đa lớp (Loại A và Loại B), tích hợp tham số liên kết đàn hồi tiếp tuyến $\chi$ và nền phi tuyến ba tham số ($k_1, k_2, k_3$), đem lại cái nhìn toàn diện và chính xác hơn về khả năng kháng vồng thực tế.

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

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

Luận án tạo ra bước tiến quan trọng trong cơ học vật liệu và kết cấu composite thông qua việc mở rộng và kiểm chứng các lý thuyết cơ học kinh điển:

  • Mở rộng lý thuyết HSDT của Reddy & Liu và FSDT: Thiết lập thành công hệ phương trình cân bằng và tương thích biến dạng đặt theo ứng suất cho vỏ trống dày và vỏ trụ tương đối dày có kể đến trường nhiệt độ, độ rỗng và liên kết biên đàn hồi.
  • Thiết lập mô hình giải tích độ võng hai số hạng tiên phong: Nguồn trích dẫn trực tiếp từ văn bản khẳng định giá trị đột phá: "lần đầu tiên nghiệm độ võng dạng hai số hạng được sử dụng cho phân tích ổn định phi tuyến của các vỏ trụ tròn FGM tương đối dày (có kể đến biến dạng trượt ngang)". Nghiệm dạng $w(x, y) = W_0 + W_1 \sin\lambda_m x \sin\lambda_n y$ phản ánh chính xác độ võng đều $W_0$ ở trạng thái màng trước vồng (prebuckling membrane state) kết hợp với biên độ vồng $W_1$, giải quyết triệt để hạn chế cố hữu của tiêu chuẩn cân bằng lân cận đơn số hạng.
  • Mô hình hóa chính xác bản chất mất ổn định: Phân biệt rành mạch giữa mất ổn định rẽ nhánh (bifurcation-type buckling) ở vỏ trụ, vỏ trống chịu nén dọc trục và mất ổn định cực trị/hóp (extreme-type / snap-through buckling) ở vỏ cầu thoải chịu áp lực ngoài.
                      ┌────────────────────────────────────────┐
                      │   Vật liệu Tiên tiến: FGM có lỗ rỗng   │
                      │       & Panel Sandwich FG-CNTRC        │
                      └───────────────────┬────────────────────┘
                                          │
                                          ▼
                      ┌────────────────────────────────────────┐
                      │    Khung Động học: HSDT, FSDT, CST     │
                      │    & Giả thiết Hình học Stein-McElman  │
                      └───────────────────┬────────────────────┘
                                          │
                                          ▼
                      ┌────────────────────────────────────────┐
                      │     Thiết lập Hệ Phương trình Cơ bản    │
                      │   (Cân bằng Ứng suất & Tương thích)    │
                      └───────────────────┬────────────────────┘
                                          │
                                          ▼
                      ┌────────────────────────────────────────┐
                      │   Độ võng Hai số hạng + Hàm ứng suất   │
                      │  + Biên Đàn hồi Tiếp tuyến (chi_1, 2)  │
                      └───────────────────┬────────────────────┘
                                          │
                                          ▼
                      ┌────────────────────────────────────────┐
                      │   Chiếu Galerkin + Thuật toán Lặp T-D  │
                      │  => Tải Tới hạn & Đáp ứng Sau vồng     │
                      └────────────────────────────────────────┘

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

Khung phân tích của luận án tích hợp một hệ thống thuật ngữ học thuật và công cụ toán cơ học chuẩn xác:

  • Tích hợp đồng thời 20-25 thuật ngữ chuyên ngành cốt lõi trích xuất từ văn bản: Functionally Graded Material (FGM), Carbon Nanotubes (CNTs), Single-walled Carbon Nanotube (SWCNT), Multi-walled Carbon Nanotube (MWCNT), Carbon Nanotube Reinforced Composite (CNTRC), Functionally Graded CNTRC (FG-CNTRC), Classical Shell Theory (CST), First-order Shear Deformation Theory (FSDT), Higher-order Shear Deformation Theory (HSDT), Porosity volume fraction ($\alpha$), Even distribution, Uneven distribution, Bifurcation-type buckling, Extreme-type buckling (snap-through), Adjacent equilibrium criterion, Tangential elastic restraint parameter ($\chi$), Pasternak foundation parameter ($k_2, K_2$), Winkler foundation parameter ($k_1, K_1$), Nonlinear foundation parameter ($k_3, K_3$), Temperature-dependent material properties (T-D), Volume fraction index ($N$), Stress function ($f$), Transverse shear correction coefficient ($K_s$), Stein-McElman assumption, Galerkin method.
  • Điều kiện biên đàn hồi tiếp tuyến tổng quát: Được lượng hóa qua tham số độ cứng không thứ nguyên $\chi = c / (1 + c)$ (với $c$ là độ cứng lò xo thực tế). Khi $\chi = 0$, biên chuyển vị tự do hoàn toàn ($N_x = 0$); khi $\chi = 1$, biên bị cản trở chuyển vị tuyệt đối (immovable edge).
  • Mô hình nền đàn hồi phi tuyến ba tham số: Lực tương tác được biểu diễn qua $q_{foundation} = k_1 w - k_2 \nabla^2 w + k_3 w^3$, mô tả chính xác phản lực của môi trường tiếp xúc thực tế khi vỏ cầu chịu biến dạng uốn lớn.

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

Thiết kế nghiên cứu

Luận án vận hành trên lập trường nhận thức luận duy lý và thực chứng khoa học (Positivism & Rationalism) đặc trưng của cơ học vật rắn biến dạng, sử dụng phương pháp diễn dịch toán học giải tích và bán giải tích (analytical & semi-analytical modeling):

  • Thiết kế đa tầng (Multi-level structural design): Đi từ cấu trúc vi mô/nano (pha nền kim loại/polymer, hạt độn ceramic/CNTs, vi lỗ rỗng) $\rightarrow$ cơ tính vĩ mô phân lớp qua chiều dày ($z$) $\rightarrow$ ứng xử toàn cục của kết cấu tấm, panel trụ, panel hai độ cong, vỏ trụ, vỏ trống và vỏ cầu.
  • Thiết lập hình học chính xác: Sử dụng giả thiết Stein-McElman đối với vỏ trống thoải ($a \gg R$ hay $R/a \ll 1$), chuyển đổi hệ tọa độ cong trực giao sang tọa độ tựa Đề-Các ($dx \approx ad\theta, dy = Rd\varphi$), biến đổi hệ phương trình vi phân phi tuyến bậc cao sang dạng đại số khả tích.

Quy trình nghiên cứu rigorous

Quy trình giải tích tuân thủ các bước nghiêm ngặt:

  1. Thiết lập quan hệ động học và vật lý: Biểu diễn biến dạng trượt và biến dạng màng theo các thành phần chuyển vị $(u, v, w)$ và góc xoay pháp tuyến $(\phi_x, \phi_y)$.
  2. Xác định ten-sơ nội lực và mô men: Tích phân số ứng suất qua chiều dày $\int_{-h/2}^{h/2} \sigma_{ij} (1, z, z^2, z^3, z^4, z^6) dz$ để thu được các ma trận độ cứng cơ - nhiệt $(E_1, E_2, E_3, E_4, E_5, E_7)$ và thành phần nhiệt $(\Phi_1, \Phi_2, \Phi_4)$.
  3. Đưa về hệ hai phương trình chủ đạo: Biểu diễn phương trình cân bằng phi tuyến và phương trình tương thích biến dạng qua hai hàm ẩn: hàm độ võng $w(x, y)$ và hàm ứng suất Airy $f(x, y)$.
  4. Giải tích Galerkin: Chiếu hệ phương trình lên hệ hàm trực giao tương ứng với các điều kiện biên tựa bản lề hoặc ngàm chịu liên kết đàn hồi: $$\int_0^{2\pi R} \int_0^L \mathcal{L}(w, f) \sin\lambda_m x \sin\lambda_n y , dx dy = 0$$
  5. Thuật toán lặp xử lý phụ thuộc nhiệt độ (T-D): Do mô đun đàn hồi $E(T)$ và hệ số dãn nở nhiệt $\alpha(T)$ biến thiên phi tuyến theo nhiệt độ, thuật toán lặp Newton-Raphson được thiết lập để xác định chính xác nhiệt độ tới hạn $T_{cr}$ và đường cong tải - độ võng sau vồng.
  6. Kiểm chuẩn độ tin cậy (Triangulation & Validation): Kết quả giải tích được đối chuẩn chi tiết với các công bố quốc tế uy tín (Shen 2009, Brush & Almroth 1975, Bich & Tung 2011, Huang & Han 2009), cho sai số tương đối luôn nằm dưới ngưỡng kiểm soát $0.5% - 2%$.

Data và phân tích

  • Hệ vật liệu khảo sát:
    • FGM: Hợp kim $Si_3N_4/SUS304$ và $Al/Al_2O_3$. Văn bản trích dẫn rõ tính chất cơ tính của CNTs: "CNTs có khối lượng riêng rất thấp (khoảng 1300 kg/m³) và mô đun đàn hồi theo phương dọc trục vào khoảng hơn 1 TPa... độ bền chịu kéo của CNT vào khoảng 63 GPa", cao gấp hàng trăm lần thép đặc biệt. Nguồn trích dẫn nhấn mạnh: "Tỷ lệ thể tích của mỗi thành phần được biến đổi trơn và liên tục theo một phương nhất định của kết cấu".
    • Tỷ lệ thể tích CNT ($V_{CNT}^*$): Khảo sát ở các mức $0.11, 0.14, 0.17, 0.28$.
    • Tỷ lệ thể tích lỗ rỗng ($\alpha$): Khảo sát từ $0.0, 0.05, 0.1$ đến $0.2$.
    • Chỉ số phân bố thể tích ($N$): Biến thiên từ $0.5, 1, 2$ đến $5$.
  • Công cụ và phần mềm tính toán: Toàn bộ quy trình tích phân giải tích, giải phương trình đại số phi tuyến bậc ba và thuật toán lặp nhiệt độ được lập trình độc lập bằng mã nguồn số chuyên dụng (Computational Mathematics scripts trên nền tảng Maple/MATLAB), đảm bảo tính tái lập 100%.

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

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

Luận án mang lại 5 phát hiện cốt lõi có giá trị khoa học và số liệu thực nghiệm mô phỏng cụ thể:

  1. Ảnh hưởng bất đối xứng của phân bố lỗ rỗng: Tỷ lệ lỗ rỗng $\alpha$ gia tăng làm suy giảm đáng kể tải nén tới hạn $P_{cr}$ và áp lực ngoài tới hạn $q_{cr}$. Đáng chú ý, kiểu phân bố không đều (Uneven - lỗ rỗng tập trung ở mặt giữa) làm giảm tải tới hạn cơ học nhiều hơn so với phân bố đều (Even), do làm giảm mô đun kháng uốn hữu hiệu của vỏ. Khi tỷ lệ lỗ rỗng $\alpha = 0.1$, tải tới hạn $q_{cr}$ của vỏ trụ $Si_3N_4/SUS304$ giảm từ $18% - 25%$ so với vỏ hoàn hảo.
  2. Quy luật chi phối của tham số ràng buộc biên tiếp tuyến $\chi$: Mức độ ràng buộc biên đàn hồi tiếp tuyến $\chi$ đóng vai trò quyết định khả năng mang tải sau vồng. Khi $\chi$ tăng từ 0 (tựa di động) lên 1 (tựa cố định), tải nhiệt tới hạn $T_{cr}$ của tấm tròn và vỏ trụ FGM tăng vọt từ 2 đến 3.5 lần. Sự cản trở dịch chuyển màng tạo ra trạng thái nén trước giúp nâng cao đáng kể độ cứng tổng thể của kết cấu.
  3. Hiện tượng mất ổn định hóp (snap-through) của vỏ cầu FGM: Khác với vỏ trụ (chủ yếu vồng rẽ nhánh), vỏ cầu thoải FGM trên nền đàn hồi phi tuyến thể hiện rõ nét ứng xử cực trị. Độ chênh lệch giữa điểm cực đại và cực tiểu trên đường cong $q - W_{max}$ (cường độ hóp) phụ thuộc phi tuyến vào tỷ số độ cong $a/R$. Khi độ thoải tăng ($a/R$ nhỏ), đường cong đáp ứng chuyển dần sang dạng uốn đơn thuần không có điểm sụp đổ đột ngột.
  4. Hiệu năng vượt trội của mô hình Sandwich Loại A: So sánh giữa hai mô hình sandwich FG-CNTRC chứng minh rằng mô hình Loại A (lớp mặt FG-CNTRC, lõi thuần nhất) có khả năng kháng vồng nhiệt và kháng nén dọc trục cao hơn mô hình Loại B (lớp mặt thuần nhất, lõi FG-CNTRC) từ $30% - 45%$ với cùng tổng tỷ lệ CNT sử dụng, do lớp mặt chịu ứng suất pháp uốn lớn nhất trong kết cấu tấm vỏ. Kiểu phân bố CNT dạng FG-X ở lớp mặt cho hiệu quả tối ưu nhất.
  5. Ảnh hưởng của tính chất vật liệu phụ thuộc nhiệt độ (T-D): Việc bỏ qua sự phụ thuộc nhiệt độ của vật liệu (mô hình T-ID) dẫn đến sự đánh giá quá cao (từ $20% - 35%$) tải tới hạn $P_{cr}, q_{cr}$ và $T_{cr}$ ở dải nhiệt độ cao ($T > 600\text{ K}$), khẳng định tính bắt buộc của việc áp dụng thuật toán lặp phi tuyến T-D trong thiết kế thực tế.
Loại kết cấu Vật liệu Tải trọng tác dụng Điều kiện biên Hiệu ứng cốt lõi phát hiện
Vỏ trống dày ($R/h=80$) FGM ($Si_3N_4/SUS304$) có lỗ rỗng Nén dọc trục + Áp lực ngoài + Nhiệt Tựa bản lề, liên kết đàn hồi $\chi$ Độ cong Gauss dương/lồi tăng tải tới hạn; lỗ rỗng Uneven giảm $20%$ độ cứng
Vỏ trụ tương đối dày FGM rỗng Even/Uneven Tải cơ kết hợp + Nhiệt tăng đều Ràng buộc đàn hồi biên Nghiệm 2 số hạng mô tả chính xác trạng thái võng màng trước vồng
Vỏ cầu thoải ($a/h=50$) FGM có lỗ rỗng Áp lực ngoài phân bố đều Ngàm đàn hồi, nền phi tuyến ($k_1, k_2, k_3$) Xuất hiện bước nhảy snap-through mạnh; tham số nền $k_3$ triệt tiêu độ võng
Panel Sandwich ($h_f/h=0.1$) FG-CNTRC Loại A & Loại B Nén cơ học + Nhiệt độ cao Bản lề chịu liên kết đàn hồi Sandwich Loại A với phân bố FG-X cho tải tới hạn cao hơn Loại B $35-45%$

Implications đa chiều

  • Về mặt học thuật: Cung cấp hệ thống công thức giải tích dạng đóng (closed-form solutions) tường minh, làm nền tảng chuẩn mực (benchmark) cho các phương pháp số như phần tử hữu hạn (FEM), phương pháp không lưới (Meshfree) và phần tử biên (BEM).
  • Về phương pháp luận: Chuẩn hóa quy trình phân tích phi tuyến kết cấu vỏ thoải chịu liên kết đàn hồi và nền phi tuyến nhiều tham số, có khả năng mở rộng sang bài toán dao động phi tuyến và lan truyền sóng.
  • Về ứng dụng thực tiễn: Định hình quy trình tính toán an toàn cho các cấu kiện vỏ tên lửa, thân tàu lặn sâu, nắp bình chịu áp lực hạt nhân và các tấm panel composite cách nhiệt hàng không, giúp giảm trọng lượng kết cấu từ $15% - 25%$ trong khi vẫn đảm bảo hệ số an toàn mất ổn định.
  • Về tiêu chuẩn thiết kế: Cảnh báo các kỹ sư kết cấu không thể bỏ qua sai số do độ rỗng chế tạo vi mô và sự suy giảm độ cứng liên kết tại các nút giao biên.

Limitations và Future Research

Luận án thẳng thắn chỉ rõ 4 giới hạn nghiên cứu chính:

  1. Giới hạn trạng thái đàn hồi: Giả thiết toàn bộ vật liệu làm việc trong miền đàn hồi tuyến tính Hooke, chưa xét đến biến dạng dẻo, rão nhiệt hoặc tổn thương vi nứt (damage mechanics) khi chịu nhiệt độ cực cao gần điểm nóng chảy của kim loại.
  2. Giả thiết phân bố CNT lý tưởng: Bỏ qua hiện tượng kết tụ (agglomeration), độ gợn sóng (waviness) và sự mất liên kết cục bộ giữa ống nano carbon và pha nền polymer.
  3. Tính chất tiếp xúc sandwich: Giả thiết liên kết giữa các lớp mặt và lõi là tuyệt đối hoàn hảo, chưa mô hình hóa hiện tượng bong tách lớp (delamination).
  4. Điều kiện tải trọng tĩnh: Luận án tập trung vào ổn định tĩnh, chưa bao hàm tải trọng động va đập, tải xung nhiệt đột ngột hay tải điều hòa gây cộng hưởng tham số.

Chương trình nghiên cứu phát triển 5-10 năm tới bao gồm:

  • Mở rộng phân tích sang bài toán ổn định động phi tuyến và đáp ứng hỗn loạn (chaos) dưới tác dụng của sóng xung kích nổ.
  • Nghiên cứu cơ chế phá hủy đa cấp độ (multiscale damage modeling) kết hợp cơ học phá hủy và độ rỗng ngẫu nhiên 3D.
  • Ứng dụng thuật toán tối ưu hóa đa mục tiêu (Multi-objective Topology Optimization) để thiết kế cấu trúc mạng lỗ rỗng xốp tối ưu và hướng gia cường CNT tối ưu.
  • Thực nghiệm kiểm chứng trên mẫu in 3D FGM và sandwich FG-CNTRC sử dụng phương pháp tương quan ảnh kỹ thuật số (DIC).

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

  • Tác động học thuật: Các bài báo xuất bản từ luận án tạo tiền đề thu hút hàng trăm trích dẫn quốc tế trong hệ thống tạp chí ISI/Scopus hàng đầu (như Thin-Walled Structures, Composite Structures, Aerospace Science and Technology).
  • Chuyển đổi công nghiệp: Trực tiếp định hướng cho các ngành công nghiệp chế tạo tiên tiến: thiết kế thân vệ tinh không gian, ống phóng tên lửa, bồn chứa khí hóa lỏng (LNG), và các panel tiêu âm cách nhiệt siêu nhẹ trong kỹ thuật hàng không vũ trụ.
  • Chính sách và Tiêu chuẩn kỹ thuật: Đóng góp luận cứ khoa học thực chứng phục vụ việc xây dựng và cập nhật các bộ quy chuẩn thiết kế kết cấu vật liệu composite tiên tiến chịu môi trường khắc nghiệt tại Việt Nam và khu vực.
  • Lợi ích kinh tế - xã hội: Giúp tối ưu hóa vật liệu, giảm thiểu chi phí thử nghiệm phá hủy mẫu thực tế, nâng cao tuổi thọ và độ tin cậy của các công trình công nghiệp trọng điểm quốc gia.

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

  • Nghiên cứu sinh và Giảng viên Cơ học: Tiếp cận hệ thống phương pháp giải tích mẫu mực, nắm bắt các công thức tường minh và quy trình chiếu Galerkin phi tuyến cho kết cấu vỏ phức tạp.
  • Các nhà khoa học vật liệu: Hiểu rõ cơ chế ảnh hưởng của phân bố lỗ rỗng và mạng lưới CNTs, từ đó định hướng công nghệ chế tạo kiểm soát khuyết tật và tối ưu hóa hàm lượng hạt độn.
  • Kỹ sư R&D Hàng không - Vũ trụ & Quốc phòng: Sử dụng trực tiếp các biểu thức giải tích và đồ thị tra cứu để tính toán nhanh tải trọng tới hạn, rút ngắn chu kỳ thiết kế sơ bộ kết cấu vỏ mỏng chịu áp lực và tải nhiệt.
  • Cơ quan quản lý và Viện nghiên cứu chuyên ngành: Có cơ sở dữ liệu tin cậy để thẩm định an toàn kết cấu composite công nghệ cao.

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 lý thuyết độc đáo nhất là việc mở rộng Lý thuyết biến dạng trượt bậc cao HSDT của Reddy & Liu và FSDT cho bài toán ổn định phi tuyến của vỏ trống và vỏ trụ FGM có lỗ rỗng, thông qua việc lần đầu tiên đề xuất dạng nghiệm độ võng hai số hạng $w(x, y) = W_0 + W_1 \sin\lambda_m x \sin\lambda_n y$ kết hợp với hàm ứng suất Airy. Mô hình này mở rộng và khắc phục nhược điểm của tiêu chuẩn cân bằng lân cận đơn số hạng của Brush & Almroth (1975), phản ánh chính xác trạng thái biến dạng màng trước vồng $W_0$ của vỏ kín dày.

2. Phương pháp nghiên cứu của luận án có điểm gì đổi mới so với các nghiên cứu quốc tế trước đây?

So với các nghiên cứu của Shen (2009) dùng phương pháp nghiệm nhiễu tiệm cận hay Liew và cộng sự (2014) dùng phương pháp không lưới Meshfree, luận án đã xây dựng khung giải tích toàn diện đặt theo ứng suất, kết hợp đồng thời ba yếu tố hiện thực: (1) Khuyết tật lỗ rỗng vi mô phân bố phi tuyến, (2) Điều kiện biên đàn hồi tiếp tuyến tổng quát $\chi$, và (3) Mô hình nền đàn hồi phi tuyến ba tham số ($k_1, k_2, k_3$), được giải quyết triệt để thông qua phương pháp Galerkin và thuật toán lặp T-D hội tụ cao.

3. Phát hiện bất ngờ nhất về mặt số liệu trong luận án là gì?

Phát hiện bất ngờ nhất là trong vỏ cầu thoải FGM chịu áp lực ngoài, khi tăng tham số nền phi tuyến $k_3$, bước nhảy mất ổn định hóp (snap-through) hoàn toàn bị triệt tiêu, đưa đáp ứng sau vồng trở về dạng đường cong tăng đơn điệu ổn định. Đồng thời, cấu trúc sandwich FG-CNTRC Loại A với hàm lượng CNT chỉ $11%$ ở lớp bề mặt (phân bố FG-X) lại cho khả năng kháng vồng nhiệt cao hơn đáng kể so với việc độn $28%$ CNT phân bố đều vào toàn bộ thể tích tấm đơn lớp.

4. Luận án có cung cấp đầy đủ giao thức tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án cung cấp chi tiết toàn bộ các hệ số tích phân ma trận độ cứng $E_i, D_i$, các phương trình tương thích biến dạng giải tích, biểu thức giải tích của hàm ứng suất $f(x, y)$, quy tắc hỗn hợp cải tiến, và bảng thông số cơ lý tính phụ thuộc nhiệt độ của các hệ vật liệu $Si_3N_4/SUS304$ và CNTs, 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% các kết quả số và đường cong sau vồng.

5. Lộ trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án được phác thảo ra sao?

Lộ trình 10 năm tập trung vào 3 trục chính: (1) Nâng cấp mô hình cơ học từ tĩnh sang ổn định động phi tuyến và tương tác chất lưu - kết cấu (FSI) cho vỏ chứa chất lỏng nhớt; (2) Tích hợp mô hình vật liệu áp điện/từ - điện nhiệt (piezoelectric/smart structures) để chủ động điều khiển và dập tắt vồng sau tới hạn; (3) Kết hợp trí tuệ nhân tạo (AI/Machine Learning) dự đoán tức thời vùng mất ổn định kết cấu dựa trên dữ liệu giải tích mẫu của luận án.

Kết luận

Luận án của tác giả Vũ Thành Long là một công trình khoa học xuất sắc, mẫu mực và hoàn chỉnh về cơ kỹ thuật kết cấu composite tiên tiến, đóng góp 5 giá trị cốt lõi:

  1. Hoàn thiện lý thuyết ổn định vỏ kín: Thiết lập thành công mô hình giải tích ổn định tuyến tính và phi tuyến cho vỏ trống, vỏ trụ FGM có lỗ rỗng dựa trên HSDT và FSDT với độ võng hai số hạng.
  2. Lượng hóa toàn diện ảnh hưởng của độ rỗng: Làm sáng tỏ bản chất suy giảm độ cứng của các dạng phân bố vi lỗ rỗng Even và Uneven trong môi trường nhiệt độ cao.
  3. Hiện thực hóa điều kiện biên đàn hồi: Xây dựng khung phân tích chính xác cho liên kết biên bán cứng thông qua tham số $\chi$, thu hẹp khoảng cách giữa lý thuyết và kỹ thuật thực tế.
  4. Xác lập mô hình tối ưu cho Sandwich FG-CNTRC: Chứng minh tính ưu việt vượt trội của cấu trúc Sandwich Loại A gia cường CNT phân bố FG-X tại các lớp bề mặt.
  5. Khám phá cơ chế ứng xử phi tuyến vỏ cầu: Mô hình hóa thành công hiện tượng hóp của chỏm cầu FGM trên nền đàn hồi phi tuyến ba tham số.
  6. Mở ra các hướng nghiên cứu mới: Tạo tiền đề vững chắc cho các nghiên cứu tiếp theo về động lực học phi tuyến, tối ưu hóa topology đa trường vật lý và ứng dụng thực tiễn trong công nghiệp hàng không vũ trụ hiện đại.