Tổng quan về luận án

Nghiên cứu về cơ học kết cấu vỏ kín đóng vai trò then chốt trong kỹ thuật hàng không vũ trụ, công nghiệp chế tạo tên lửa, công trình biển và các bồn bể chịu áp lực cao. Luận án tiến sĩ kỹ thuật của tác giả Phạm Thanh Hiếu (người hướng dẫn: PGS.TS Hoàng Văn Tùng và PGS.TS Đào Như Mai, thực hiện tại Viện Cơ học thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam) mang tiêu đề "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" là một công trình tiên phong giải quyết bài toán ứng xử cơ - nhiệt phức tạp của kết cấu vỏ thoải xoay tròn.

Khoảng trống học thuật cốt lõi (research gap) mà luận án giải quyết nằm ở ba điểm nghẽn chính trong cơ học kết cấu tính toán hiện đại:

  1. Thiếu mô hình phân tích giải tích phi tuyến chính xác cho vỏ kín chịu tải nhiệt và cơ - nhiệt kết hợp có xét đến điều kiện biên đàn hồi thực tế (partially movable edges), khắc phục giả thiết lý tưởng hóa cứng nhắc (hoàn toàn tự do $c=0$ hoặc ngàm cứng $c=\infty$) thường thấy trong các nghiên cứu kinh điển của Brush & Almroth (1975) hay Shen (2009).
  2. Sự khan hiếm nghiêm trọng các phân tích giải tích về ứng xử vồng (buckling) và sau vồng (postbuckling) của vỏ trống (barrel shell) có độ cong Gauss biến thiên làm từ vật liệu composite gia cường ống nano các-bon có cơ tính biến đổi (FG-CNTRC).
  3. Hạn chế toán học của dạng nghiệm ba số hạng truyền thống khi áp dụng lý thuyết biến dạng trượt bậc nhất (FSDT) cho bài toán ổn định tuyến tính của vỏ kín tương đối dày.

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

  • Câu hỏi nghiên cứu 1 (RQ1): Tính đàn hồi trong liên kết biên tiếp tuyến tác động như thế nào đến tải tới hạn và đường cân bằng sau vồng của vỏ trụ sandwich FGM và vỏ FG-CNTRC dưới tải áp lực ngoài và nhiệt độ tăng đều?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố CNT (UD, FG-$\Lambda$, FG-V, FG-O, FG-X) kết hợp với độ cong Gauss của vỏ trống tối ưu hóa khả năng chịu tải cơ - nhiệt kết hợp ra sao?
  • Giả thuyết 1 (H1): Điều kiện biên đàn hồi tiếp tuyến trung gian ($0 < c < \infty$) tạo ra sự dịch chuyển đáng kể của điểm rẽ nhánh và thay đổi cường độ hóp (snap-through) so với các biên lý tưởng.
  • Giả thuyết 2 (H2): Dạng nghiệm độ võng hai số hạng trong khuôn khổ FSDT cho phép giải quyết chính xác bài toán ổn định tuyến tính của vỏ kín FG-CNTRC tương đối dày mà không bị suy biến điều kiện biên đóng kín theo phương vĩ tuyến.

Phạm vi nghiên cứu bao quát các tỷ số hình học từ mỏng đến tương đối dày ($R/h = 30 \div 400$, $L/R = 1 \div 3$, $L^2/Rh = 300$), nhiệt độ môi trường từ $300\text{ K}$ đến $500\text{ K}$, với các tham số độn hạt nano SWCNT (10,10) có mô đun đàn hồi $E > 1\text{ TPa}$ và tỷ lệ thể tích $V_{CNT}^* = 0.11 \div 0.28$. Ý nghĩa đột phá của luận án là cung cấp hệ thống công thức giải tích tường minh, đóng vai trò làm chuẩn kiểm định (benchmark) tin cậy cho các thuật toán số (FEM, Meshfree, DQM).

Literature Review và Positioning

Lịch sử nghiên cứu vật liệu tiên tiến ghi nhận bước chuyển dịch từ vật liệu cơ tính biến đổi truyền thống FGM (cấu tạo từ ceramic như $\text{Si}_3\text{N}_4$ và kim loại như SUS304 nhằm khắc phục sự tập trung ứng suất và tách lớp của composite cổ điển) sang composite gia cường ống nano các-bon FG-CNTRC được Shen khởi xướng năm 2009. Ống nano các-bon (CNTs) từ phát hiện của Iijima (1991) đã chứng minh các đặc tính vượt bậc: khối lượng riêng siêu nhẹ xấp xỉ $1300\text{ kg/m}^3$ (thấp hơn nhiều so với nhôm $2700\text{ kg/m}^3$ hay thép $7800\text{ kg/m}^3$), độ bền kéo đo được lên tới $63\text{ GPa}$ (gấp 10–100 lần thép cường độ cao), độ ổn định nhiệt đạt $2800^\circ\text{C}$ trong chân không và độ dẫn nhiệt gấp đôi kim cương.

Tổng quan tài liệu cho thấy hai luồng quan điểm và tranh luận học thuật sâu sắc:

  1. Tranh luận về mô hình hóa cơ tính tương đương: Một trường phái (Shen, 2009; Trang & Tung, 2018) ủng hộ quy tắc hỗn hợp suy rộng (extended rule of mixtures) với các tham số hiệu quả $\eta_1, \eta_2, \eta_3$ hiệu chỉnh bằng mô phỏng động lực học phân tử (MD); trong khi trường phái khác (Macías et al., 2019) ưu tiên mô hình vi cơ học Eshelby-Mori-Tanaka. Luận án định vị tiếp cận theo quy tắc hỗn hợp suy rộng phụ thuộc nhiệt độ nhằm bảo đảm tính tương thích giải tích chính xác.
  2. Tranh luận về ảnh hưởng của biến dạng trước vồng (pre-buckling deformation): Các công trình cổ điển thường bỏ qua biến dạng trước vồng để đơn giản hóa phương trình trạng thái màng. Tuy nhiên, đối với vỏ kín chịu tải nén hoặc áp lực ngoài, Brush & Almroth (1975) và Hoàng Văn Tùng (2014) khẳng định việc bỏ qua biến dạng trước vồng gây sai số nghiêm trọng trong việc dự đoán điểm rẽ nhánh và hiện tượng hóp tức thì (immediate snap-through).

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

  • So sánh với Shen (2011, 2014): Shen sử dụng lý thuyết HSDT kết hợp phương pháp khai triển tiệm cận hai bước (two-step perturbation technique) để khảo sát vỏ trụ FG-CNTRC nhưng chỉ giới hạn ở điều kiện biên tựa di động hoàn toàn (movable) hoặc tựa cố định (immovable), bỏ qua tính đàn hồi biên thực tế và không giải quyết kết cấu vỏ trống.
  • So sánh với Zenkour (2005) và Sofiyev et al. (2017): Các tác giả này tập trung vào bài toán tuyến tính và dao động tự do của tấm/vỏ FGM sandwich bằng FSDT và nghiệm chuỗi Navier, chưa thiết lập hệ nghiệm giải tích phi tuyến đa số hạng kết hợp phương pháp Galerkin có xét đến hàm ứng suất cho vỏ có liên kết lò xo biên tiếp tuyến.

Luận án của Phạm Thanh Hiếu đã định vị chuẩn xác điểm giao thoa giữa cơ học vật liệu nano và phi tuyến hình học kết cấu, mở rộng lý thuyết cho cả hai dạng hình học: vỏ trụ tròn và vỏ trống mỏng/dày với liên kết biên tiếp tuyến đàn hồi $N_{x0} = c \Delta_x$.

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

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

Luận án mở rộng nền tảng của Lý thuyết vỏ Donnell (Donnell Shell Theory) và quan hệ phi tuyến động học von Kármán – Donnell bằng việc tích hợp các liên kết đàn hồi biên tiếp tuyến và sự phụ thuộc nhiệt độ của vật liệu thành phần FGM và FG-CNTRC:

[Vật liệu tiên tiến: FGM Sandwich / FG-CNTRC (SWCNT 10,10)]
[Mô hình biến dạng: CST (Vỏ mỏng) / FSDT (Vỏ dày)]
[Phi tuyến hình học von Kármán - Donnell + Nền đàn hồi Pasternak]
[Điều kiện biên đàn hồi tiếp tuyến: Nx0 = c · Δx (0 ≤ c ≤ ∞)]
[Hệ phương trình vi phân phi tuyến: Cân bằng (w) & Tương thích (f)]
[Phương pháp Galerkin & Thuật toán lặp nhiệt độ]
[Phát hiện: Tải tới hạn, Quỹ đạo sau vồng, Hiện tượng hóp, Miền ổn định]

Mô hình lý thuyết thiết lập 4 mệnh đề toán - cơ học then chốt:

  • Mệnh đề 1 (Độ cứng tiếp tuyến đàn hồi): Tham số độ cứng biên tiếp tuyến không thứ nguyên $\bar{c} = c L / (E_m h)$ chi phối trực tiếp lực nén bị động $N_{x0}$ sinh ra do hiệu ứng giãn nở nhiệt hoặc giãn nở ngang Poisson, quyết định sự dịch chuyển từ mất ổn định kiểu rẽ nhánh (bifurcation-type) sang mất ổn định kiểu giới hạn (limit-type).
  • Mệnh đề 2 (Cấu hình sandwich FGM): Vỏ sandwich Loại A (lõi kim loại, hai mặt FGM giàu gốm đối xứng) có tính đối xứng cơ học cao hơn, triệt tiêu sự ghép nối uốn - kéo nội tại ($B_{ij} = 0$), mang lại khả năng kháng vồng nhiệt vượt trội so với vỏ Loại B (mặt ngoài gốm, lõi FGM, mặt trong kim loại bất đối xứng).
  • Mệnh đề 3 (Quy luật phân bố FG-X): Phân bố CNT dạng FG-X (tập trung mật độ cao ở hai bề mặt ngoài $z = \pm h/2$) tạo ra mô men quán tính uốn cục bộ cực đại, nâng cao độ cứng uốn hiệu dụng $D$ lên mức tối ưu so với các dạng UD, FG-O, FG-V và FG-$\Lambda$.
  • Mệnh đề 4 (Nghiệm FSDT hai số hạng cho vỏ dày): Việc áp dụng dạng hàm độ võng hai số hạng kết hợp hàm góc xoay trong FSDT giải quyết trọn vẹn sự tương thích biến dạng cắt ngang mà vẫn bảo đảm điều kiện khép kín theo chu vi vĩ tuyến $\int_0^{2\pi R} \int_0^L \varepsilon_{y0} dx dy = 0$.

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

Khung phân tích tích hợp đồng thời 3 hệ lý thuyết:

  1. Lý thuyết vật liệu vi cơ học có cơ tính biến đổi: Mô hình hóa hàm biến thiên tỷ lệ thể tích $V_M(z) = [ (2z+h)/(2h_f) ]^N$ cho FGM và $V_{CNT}(z) = 2(2|z|/h)V_{CNT}^*$ cho FG-X CNTRC.
  2. Lý thuyết vỏ phi tuyến von Kármán – Donnell: Thiết lập hệ phương trình đạo hàm riêng phi tuyến bậc bốn liên kết giữa hàm độ võng $w(x,y)$ và hàm ứng suất Airy $f(x,y)$: $$D \nabla^4 w - f_{,yy} w_{,xx} + 2 f_{,xy} w_{,xy} - f_{,xx} w_{,yy} - \frac{f_{,xx}}{R} - q = 0$$ $$\nabla^4 f + E_1 \left( w_{,xy}^2 - w_{,xx} w_{,yy} - \frac{w_{,xx}}{R} \right) = 0$$
  3. Lý thuyết tương tác kết cấu - nền đàn hồi hai tham số Pasternak: Mô hình hóa môi trường đàn hồi bao quanh bằng hệ số độ cứng Winkler $k_1$ và độ cứng lớp trượt $k_2$ (chuyển đổi sang tham số không thứ nguyên $K_1 = k_1 R^2 / D$, $K_2 = k_2 / D$).

Điều kiện biên tiếp xúc bao quát: Giả thiết vỏ hoàn hảo, tiếp xúc liên tục toàn phần với nền Pasternak, bán kính cong kinh tuyến của vỏ trống rất thoải ($R_x \gg R_y = R$), và vật liệu làm việc hoàn toàn trong miền đàn hồi tuyến tính Hooke.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ triết lý diễn dịch toán học (mathematical deductive paradigm) thuộc trường phái thực chứng thực nghiệm kết hợp lý thuyết giải tích cơ học chính xác (analytical mechanics).

Thiết kế nghiên cứu đa tầng bao gồm:

  • Tầng 1 - Cơ học vi mô (Micro-scale): Tính toán cơ tính nhiệt - đàn hồi hiệu dụng phụ thuộc nhiệt độ $T$ của pha thành phần: ceramic ($\text{Si}3\text{N}4$), kim loại (SUS304), polymer (PmPV, PMMA) và ống nano SWCNT (10,10) có chiều dài $L_t = 9.26\text{ nm}$, bán kính $r_t = 0.678\text{ nm}$, $E{11}^{CNT} = 5.6466\text{ TPa}$, $E{22}^{CNT} = 7.0800\text{ TPa}$, $G_{12}^{CNT} = 1.9445\text{ TPa}$.
  • Tầng 2 - Cơ học vĩ mô kết cấu mỏng (Macro-scale Thin Shells): Sử dụng Lý thuyết vỏ cổ điển (CST) khảo sát vỏ sandwich FGM và vỏ FG-CNTRC mỏng ($R/h \ge 80$).
  • Tầng 3 - Cơ học kết cấu dày có kể đến trượt ngang (Macro-scale Thick Shells): Áp dụng Lý thuyết biến dạng trượt bậc nhất (FSDT) cho dải kích thước $R/h = 30 \div 50$.

Quy trình nghiên cứu rigorous

Quy trình biến đổi toán học và giải tích nghiệm được chuẩn hóa nghiêm ngặt:

  1. Thiết lập dạng nghiệm hàm độ võng ba số hạng: $$w(x,y) = W_0 + W_1 \sin(\alpha_m x) \sin(\beta_n y) + W_2 \sin(2\alpha_m x)$$ trong đó $\alpha_m = m\pi / L$, $\beta_n = n / R$; $W_0$ biểu thị độ võng đều giai đoạn trước vồng, $W_1$ là biên độ vồng tuyến tính, $W_2$ là biên độ thành phần đối xứng trục phi tuyến sau vồng.
  2. Nghiệm giải tích hàm ứng suất Airy: $$f(x,y) = A_1 \cos(2\alpha_m x) + A_2 \cos(2\beta_n y) + A_3 \sin(\alpha_m x)\sin(\beta_n y) + A_4 \sin(3\alpha_m x)\sin(\beta_n y) - \frac{\sigma_{0y} h x^2}{2} - \frac{N_{x0} y^2}{2}$$ Các hệ số $A_1 \div A_4$ được xác định tường minh từ phương trình tương thích biến dạng.
  3. Thực thi thuật toán trực giao Galerkin: $$\int_0^L \int_0^{2\pi R} \Phi(w, f, q, T) \sin(\alpha_m x) \sin(\beta_n y) dy dx = 0$$ dẫn xuất hệ phương trình phi tuyến đại số bậc ba xác định tải tới hạn $P_{cr}, q_{cr}, T_{cr}$ và đường cong cân bằng sau vồng.
  4. Kiểm định độ tin cậy và giá trị hội tụ (Validation & Reliability): So sánh đối chuẩn (cross-validation) đa nguồn với các công bố uy tín quốc tế:
  • So sánh áp lực ngoài tới hạn $q_{cr}$ (MPa) của vỏ trụ FGM $\text{Si}_3\text{N}_4/\text{SUS304}$ ($R/h=300, T_0=300\text{ K}, \Delta T=200\text{ K}, m=1$): Sai số tuyệt đối so với kết quả của Shen (2014) dưới $0.85%$.
  • So sánh lực nén tới hạn $F_{cr} = 2\pi R h P_{cr}$ (kN) của vỏ trụ CNTRC ($R/h=100, h=1\text{ mm}, T=300\text{ K}$): Trùng khớp hoàn toàn với nghiệm số phần tử hữu hạn của Macías et al. (2019) và nghiệm giải tích của Shen & Xiang (2014).

Data và phân tích

Mô hình hóa số liệu được lập trình và thực thi bằng công cụ tính toán biểu tượng hóa học thuật (Maple/MATLAB).

Các tập tham số chuẩn xác định cấu hình phân tích:

  • Tỷ lệ kích thước: $R/h \in [30, 400]$, $L/R \in [1, 3]$, tham số vỏ hình học $L^2/(Rh) = 300$.
  • Phân bố nhiệt độ: Nhiệt độ phòng ban đầu $T_0 = 300\text{ K}$; miền khảo sát $\Delta T = 0 \div 500\text{ K}$; xét hai kịch bản vật liệu phụ thuộc nhiệt độ (T-D) và không phụ thuộc nhiệt độ (T-ID).
  • Chỉ số thể tích gốm/kim loại $N \in [0, 5]$; tỷ lệ chiều dày lớp mặt sandwich $h_f/h \in [0.1, 0.4]$.
  • Tỷ lệ thể tích CNT: $V_{CNT}^* \in {0.11, 0.14, 0.17, 0.28}$.
  • Độ cứng biên đàn hồi: Khảo sát liên tục từ $\bar{c} = 0$ (tự do) qua các giá trị hữu hạn $\bar{c} = 10^{-2}, 10^{-1}, 1, 10^1, 10^2$ đến $\bar{c} \to \infty$ (ngàm cứng tuyệt đối).

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

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

Luận án phát hiện 5 quy luật cơ học đột phá có ý nghĩa nền tảng:

           ÁP LỰC NGOÀI / LỰC NÉN DỌC TRỤC
[Kiểu phân bố CNT: FG-X]         [Mức độ liên kết biên c]
  Tăng 35% - 48% tải vồng          c tăng -> Tải vồng tăng
  Kháng uốn mặt ngoài tối đa      Chuyển đổi hình thái hóp
             [Môi trường nhiệt độ cao T]
          Vật liệu suy giảm mô đun E(T)
          Tải tới hạn giảm 18% - 32%
          Hiện tượng hóp tức thì gia tăng
  1. Hiệu ứng gia cường vượt trội của cấu hình FG-X: Trong mọi trường hợp tải (nén dọc trục, áp lực ngoài, cơ - nhiệt kết hợp), kiểu phân bố FG-X luôn đem lại tải tới hạn $P_{cr}$ và $q_{cr}$ cao nhất, vượt trội từ $35%$ đến $48%$ so với phân bố đều (UD) và vượt hơn $65%$ so với phân bố FG-O. Nguyên nhân do mật độ CNT tập trung cao nhất tại hai biên bề mặt $z = \pm h/2$, tối đa hóa mô men kháng uốn của mặt cắt.
  2. Quy luật chi phối của liên kết biên đàn hồi tiếp tuyến ($\bar{c}$): Khi $\bar{c}$ tăng từ $0$ lên $10^2$, lực tới hạn nén dọc trục $P_{cr}$ tăng mạnh nhưng đối với tải nhiệt thuần túy, $\bar{c}$ tăng lại làm giảm mạnh nhiệt độ tới hạn $T_{cr}$. Sự ngăn cản dịch chuyển tiếp tuyến khi nhiệt độ tăng sinh ra lực nén màng bị động $N_{x0}$ cực lớn, đẩy vỏ vào trạng thái mất ổn định nhiệt sớm hơn.
  3. Ảnh hưởng phi tuyến của độ cong Gauss trên vỏ trống: Vỏ trống có tỷ số bán kính $R/a > 0$ (độ cong Gauss dương) thể hiện khả năng chịu áp lực ngoài cao hơn vỏ trụ tròn tương đương ($R/a = 0$) từ $15%$ đến $28%$. Độ cong kinh tuyến làm tăng độ cứng chống méo hình học, giảm biên độ độ võng sau vồng.
  4. Hiện tượng suy giảm độ ổn định do phụ thuộc nhiệt độ (T-D effect): Tính toán theo mô hình vật liệu phụ thuộc nhiệt độ (T-D) cho giá trị tải tới hạn thấp hơn từ $18%$ đến $32%$ so với mô hình bỏ qua ảnh hưởng nhiệt độ (T-ID). Việc bỏ qua tính chất T-D dẫn đến sự đánh giá quá cao biên độ an toàn của kết cấu trong môi trường nhiệt.
  5. Cơ chế hóp tức thì (immediate snap-through) trong miền sau vồng: Khi vỏ trụ chịu tải cơ kết hợp trong môi trường nhiệt độ cao ($T = 500\text{ K}$), đường cân bằng sau vồng xuất hiện bước nhảy biến dạng đột ngột (drop-down snap-through) với cường độ hóp gia tăng tỷ lệ thuận với tỷ số bán kính trên chiều dày $R/h$.

Implications đa chiều

  • Đóng góp lý thuyết: Mở rộng thành công lý thuyết phi tuyến Donnell cho hệ vỏ kín dị hướng nano composite FG-CNTRC với liên kết biên tổng quát, lấp đầy khoảng trống tồn tại nhiều thập kỷ giữa lý thuyết biên lý tưởng và thực tế kỹ thuật.
  • Đổi mới phương pháp luận: Xây dựng quy trình tích hợp giữa phương pháp giải tích đa số hạng Galerkin và thuật toán lặp xác định điểm cân bằng nhiệt, giải quyết triệt để bài toán phi tuyến phụ thuộc tham số phi cơ học $E(T), \alpha(T)$.
  • Ứng dụng thực tiễn: Cung cấp biểu đồ thiết kế tối ưu hóa tỷ lệ phân bố thể tích $V_{CNT}^*$ và chiều dày lớp mặt $h_f/h$ cho các kỹ sư chế tạo vỏ bọc động cơ tên lửa, thân tàu ngầm mini và bồn chứa khí hóa lỏng chịu nhiệt - áp suất cao.
  • Khuyến nghị chính sách và tiêu chuẩn kỹ thuật: Đề xuất cập nhật tiêu chuẩn thiết kế vỏ mỏng chịu áp lực trong ngành giao thông vận tải và hàng không vũ trụ, bắt buộc đưa tham số độ cứng đàn hồi liên kết gối đỡ $c$ và mô hình nhiệt - đàn hồi T-D vào quy chuẩn kiểm tra an toàn chịu tải.

Limitations và Future Research

Luận án thẳng thắn chỉ ra 4 giới hạn nghiên cứu (boundary conditions):

  1. Giới hạn hình học và khuyết tật ban đầu: Mô hình giải tích giả thiết vỏ có dạng hình học hoàn hảo (perfect geometry), chưa xét đến ảnh hưởng của độ méo hình học ban đầu (initial geometric imperfections) vốn có thể làm giảm tải tới hạn thực nghiệm.
  2. Giới hạn trường biến dạng phi đàn hồi: Toàn bộ phân tích giả định vật liệu làm việc trong miền đàn hồi tuyến tính vi mô, chưa tính đến các hiện tượng dẻo hóa của nền polymer, rão nhiệt (thermal creep) hoặc sự trượt/tách dính cục bộ giữa bề mặt ống CNT và pha nền matrix.
  3. Giới hạn dạng vỏ trống: Nghiên cứu tập trung vào vỏ trống rất thoải ($R_x \gg R_y$), chưa bao quát các dạng vỏ có độ cong hai chiều bất đối xứng phức tạp hoặc vỏ nón cụt có góc côn lớn.
  4. Giới hạn điều kiện nhiệt: Chỉ khảo sát trường nhiệt độ tăng đều và truyền nhiệt đối xứng trục một chiều qua chiều dày, chưa mô hình hóa trường nhiệt độ biến đổi ngẫu nhiên dạng 3D hoặc sốc nhiệt động lực học (thermal shock).

Chương trình nghiên cứu phát triển tương lai (Future Research Agenda):

  • Phát triển mô hình giải tích phi tuyến cho vỏ FGM/FG-CNTRC có gân gia cường nghiêng (stiffened shells) chịu tải sốc nhiệt động lực học và va chạm cơ học tốc độ cao.
  • Tích hợp lý thuyết biến dạng trượt bậc cao HSDT kết hợp mô hình phần tử hữu hạn đẳng hình học (IGA) để phân tích dao động phi tuyến và đáp ứng chuyển tiếp của vỏ composite có khuyết tật hình học ngẫu nhiên.
  • Nghiên cứu cơ chế phá hủy đa cấp độ (multiscale damage mechanics) kết hợp giữa phương pháp động lực học phân tử (MD) và mô hình liên tục để dự đoán hiện tượng tách lớp nano và đứt gãy CNT trong pha nền.

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

  • Tác động học thuật (Academic Impact): Luận án đã công bố nhiều công trình trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (Q1/Q2) chuyên ngành cơ học kết cấu (Composite Structures, Thin-Walled Structures, Applied Mathematical Modelling). Ước tính các công thức giải tích chuẩn của luận án sẽ thu hút lượng trích dẫn cao từ cộng đồng nghiên cứu cơ học tính toán toàn cầu.
  • Chuyển đổi công nghiệp (Industry Transformation): Định hình tiêu chuẩn thiết kế mới cho ngành vật liệu composite tiên tiến, trực tiếp phục vụ các tập đoàn công nghiệp quốc phòng, hàng không vũ trụ và đóng tàu trong việc chế tạo các cấu kiện vỏ nhẹ, siêu bền và chịu nhiệt cao.
  • Ý nghĩa xã hội và môi trường: Việc tối ưu hóa vật liệu FG-CNTRC giúp giảm khối lượng kết cấu lên đến $25% - 30%$ so với hợp kim truyền thống, giúp tiết kiệm đáng kể nhiên liệu vận hành trong ngành vận tải hàng không và giảm phát thải khí nhà kính.

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

  • Nghiên cứu sinh và học viên cao học ngành Cơ học/Kỹ thuật công trình: Tiếp cận khung phương pháp luận toán học chặt chẽ về lý thuyết vỏ phi tuyến và kỹ thuật giải phương trình đạo hàm riêng phi tuyến bằng phương pháp Galerkin.
  • Các nhà khoa học và chuyên gia cơ học kết cấu: Sử dụng hệ thống nghiệm giải tích tường minh làm cơ sở đối chuẩn (benchmark solutions) để kiểm tra độ tin cậy của các phần mềm mô phỏng thương mại (Abaqus, Ansys) và các thuật toán số mới.
  • Kỹ sư R&D trong ngành hàng không vũ trụ và dầu khí: Ứng dụng trực tiếp các biểu đồ miền ổn định cơ - nhiệt để lựa chọn cấu hình vật liệu sandwich và phân bố CNT tối ưu cho bồn chứa áp lực cao.
  • Cơ quan quản lý và thiết lập quy chuẩn kỹ thuật: Có cơ sở khoa học định lượng để hoàn thiện các tiêu chuẩn an toàn cho kết cấu composite làm việc trong môi trường khắc nghiệt.

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?

Trả lời: Đóng góp độc đáo nhất là việc thiết lập thành công mô hình giải tích phi tuyến xác định ứng xử sau vồng của vỏ trụ sandwich FGM và vỏ kín FG-CNTRC có xét đến tính đàn hồi của liên kết biên tiếp tuyến ($N_{x0} = c \Delta_x$). Luận án đã mở rộng Lý thuyết vỏ mỏng Donnell và lý thuyết phi tuyến von Kármán, vượt qua ranh giới của các giả thiết biên lý tưởng hóa (tựa tự do hoàn toàn hoặc ngàm cứng tuyệt đối) tồn tại trong các công trình kinh điển của Shen (2009) và Brush & Almroth (1975).

2. Đổi mới phương pháp luận nổi bật khi so sánh với ít nhất 2 nghiên cứu quốc tế?

Trả lời:

  • So với nghiên cứu của Shen & Xiang (2014) sử dụng phương pháp khai triển tiệm cận hai bước (perturbation method) vốn phức tạp và nhạy cảm với việc chọn tham số bé, luận án sử dụng nghiệm hàm độ võng ba số hạng kết hợp trực giao Galerkin giúp thu được nghiệm tường minh trực tiếp cho toàn bộ ba giai đoạn (trước vồng, vồng và sau vồng).
  • So với nghiên cứu phần tử hữu hạn số của Macías et al. (2019), luận án giải quyết trọn vẹn sự phụ thuộc liên tục của cơ tính vào nhiệt độ $E(T), \alpha(T)$ thông qua thuật toán lặp giải tích chính xác, loại bỏ hoàn toàn sai số rời rạc hóa lưới không gian.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì?

Trả lời: Kết quả bất ngờ nhất là tác động trái ngược của độ cứng biên đàn hồi $\bar{c}$ đối với hai loại tải trọng: Khi $\bar{c}$ tăng từ $0$ lên $10^2$, khả năng kháng vồng do nén trục tăng đáng kể, nhưng nhiệt độ tới hạn $T_{cr}$ lại sụt giảm nghiêm trọng (giảm tới hơn $40%$). Điều này chứng minh rằng việc gia cố ngàm biên quá cứng trong môi trường nhiệt độ cao sẽ phản tác dụng, gây mất ổn định nhiệt sớm do tích tụ ứng suất nén nhiệt nội tại cực lớn.

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

Trả lời: Luận án cung cấp đầy đủ và chi tiết toàn bộ quy trình tái lập: Từ biểu thức giải tích các tích phân Galerkin, hệ số ma trận trong phụ lục toán học, hệ phương trình xác định ứng suất màng $\sigma_{0y}$ và lực nén $N_{x0}$, đến bảng thông số vật liệu phụ thuộc nhiệt độ của $\text{Si}_3\text{N}_4$, SUS304, PmPV, PMMA và cấu trúc hạt SWCNT (10,10). Mọi nhà nghiên cứu độc lập đều có thể lập trình tái lập chính xác các kết quả số.

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

Trả lời: Lộ trình 10 năm tập trung vào 3 hướng: (1) Mở rộng bài toán sang miền mất ổn định động lực học và dao động phi tuyến dưới tác dụng của sóng xung kích và sốc nhiệt; (2) Phân tích cơ chế suy giảm độ bền do mỏi nhiệt và rão vật liệu trong môi trường bức xạ nhiệt khắc nghiệt; (3) Tích hợp học máy (Physics-Informed Neural Networks - PINNs) để tối ưu hóa đa mục tiêu cấu trúc vi mô phân bố CNT và vật liệu FGM xốp (porous materials).

Kết luận

Luận án tiến sĩ của tác giả Phạm Thanh Hiếu đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 6 đóng góp mang tính bước ngoặt:

  1. Xây dựng thành công khung phân tích giải tích phi tuyến hoàn chỉnh cho vỏ trụ sandwich FGM và vỏ FG-CNTRC chịu tải cơ, nhiệt và cơ - nhiệt kết hợp có xét đến liên kết biên tiếp tuyến đàn hồi.
  2. Chứng minh định lượng cấu hình phân bố ống nano các-bon FG-X mang lại khả năng kháng vồng và chịu tải sau vồng tối ưu nhất trong tất cả các dạng phân bố nanocomposite.
  3. Làm sáng tỏ vai trò hai mặt của độ cứng liên kết biên $\bar{c}$ và xác lập ranh giới chuyển dịch giữa mất ổn định rẽ nhánh và mất ổn định giới hạn.
  4. Đề xuất dạng nghiệm hai số hạng đột phá trong khuôn khổ FSDT cho vỏ kín FG-CNTRC tương đối dày, khắc phục hạn chế lịch sử của dạng nghiệm ba số hạng.
  5. Cung cấp bộ dữ liệu chuẩn mực có độ chính xác cao về tải tới hạn và quỹ đạo sau vồng cho vỏ trống có độ cong Gauss biến thiên.
  6. Thiết lập quy trình tính toán lặp giải tích chính xác phản ánh sự suy giảm cơ tính theo nhiệt độ của vật liệu composite nano.

Công trình tạo tiền đề khai mở các hướng nghiên cứu chuyên sâu về cơ học cấu trúc nano đa tỷ lệ, đóng góp trực tiếp vào sự phát triển của khoa học cơ học tính toán và công nghệ chế tạo vật liệu tiên tiến tại Việt Nam cũng như trên trường quốc tế.