Tổng quan về luận án

Nghiên cứu ổn định tĩnh của các kết cấu vỏ tròn xoay đóng vai trò cốt lõi trong cơ học công trình hiện đại, đặc biệt trong các ngành kỹ thuật hàng không, công nghệ tên lửa, công trình biển và chế tạo bình áp lực cao. Luận án tiến sĩ mang tên "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" của tác giả Phạm Thanh Hiếu (Chuyên ngành: Cơ kỹ thuật, Mã số: 9 52 01 01; Người hướng dẫn: PGS.TS. Đào Như Mai và PGS.TS. Hoàng Văn Tùng, bảo vệ 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, 2022) đại diện cho một công trình nghiên cứu tiên phong, giải quyết đồng thời các bài toán phi tuyến phức tạp về vật liệu tiên tiến và điều kiện biên thực tế.

Trong lịch sử cơ học kết cấu, các phân tích ổn định vỏ kín thường đặt trên các giả thiết biên lý tưởng hóa: hoặc là hoàn toàn tự do di động theo phương tiếp tuyến (freely movable, ứng với độ cứng liên kết $c = 0$), hoặc ngàm cố định tuyệt đối (fully immovable, ứng với độ cứng liên kết $c \to \infty$). Tuy nhiên, trong thực tế kỹ thuật, các kết cấu vỏ trụ và vỏ trống luôn chịu sự ràng buộc đàn hồi hữu hạn từ các kết cấu lân cận ($0 < c < \infty$). Khoảng trống nghiên cứu (research gap) này chưa từng được khảo sát một cách hệ thống và giải tích toàn diện cho các loại vỏ kín chế tạo từ vật liệu cơ tính biến đổi (FGM) và composite gia cường ống nano carbon có cơ tính biến đổi (FG-CNTRC) chịu tải nhiệt - cơ kết hợp, đặc biệt khi tính đến biến dạng trước vồng và hiện tượng hóp (snap-through).

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

  • RQ1: Điều kiện biên đàn hồi tiếp tuyến ($c$) ảnh hưởng như thế nào đến tải tới hạn và đáp ứng sau vồng của vỏ trụ sandwich FGM chịu áp lực ngoài và nhiệt độ?
    • H1: Sự gia tăng độ cứng biên đàn hồi làm tăng đáng kể áp lực tới hạn nhưng làm thay đổi quỹ đạo sau vồng, kích hoạt hiện tượng hóp nhiệt đàn hồi tức thì.
  • RQ2: Quy luật phân bố ống nano carbon (CNT) theo chiều dày và độ cong Gauss ảnh hưởng ra sao đến miền ổn định của vỏ trụ và vỏ trống mỏng FG-CNTRC?
    • H2: Cấu hình phân bố FG-X tối ưu hóa khả năng chịu tải nén và áp lực ngoài vượt trội so với phân bố đều (UD) và phân bố FG-O nhờ tối đa hóa mô men kháng uốn tại các thớ ngoài; độ cong Gauss dương của vỏ trống nâng cao đáng kể độ cứng tổng thể.
  • RQ3: Tính chất vật liệu phụ thuộc nhiệt độ (Temperature-Dependent - T-D) và nền đàn hồi Pasternak tương tác ra sao với tính phi tuyến hình học của vỏ?
    • H3: Bỏ qua tính phụ thuộc nhiệt độ dẫn đến đánh giá quá cao khả năng kháng vồng (sai số từ 15% đến 35%), trong khi nền hai tham số Pasternak kiểm soát hiệu quả biên độ võng sau vồng.
  • RQ4: Làm thế nào để giải quyết bài toán ổn định tuyến tính của vỏ kín FG-CNTRC tương đối dày có xét đến biến dạng trượt ngang mà không gặp bế tắc giải tích của nghiệm ba số hạng truyền thống?
    • H4: Thiết lập dạng nghiệm hàm võng hai số hạng cải tiến kết hợp Lý thuyết biến dạng trượt bậc nhất (FSDT) cho phép xác định chính xác tải tới hạn của vỏ kín dày với độ chính xác tiệm cận mô phỏng số.

Khung lý thuyết của luận án tích hợp Lý thuyết vỏ cổ điển Donnell (CST), Lý thuyết biến dạng trượt bậc nhất (FSDT), quan hệ biến dạng - chuyển vị phi tuyến von Kármán – Donnell, quy tắc hỗn hợp suy rộng bổ sung tham số hiệu quả CNT ($\eta_1, \eta_2, \eta_3$) và phương pháp giải tích Galerkin kết hợp thuật toán lặp phi tuyến.

Về quy mô và phạm vi khảo sát, nghiên cứu phân tích toàn diện trên các dải thông số hình học rộng ($R/h = 30 \div 400$, $L/R = 1 \div 3$), tỷ lệ thể tích CNT ($V_{CNT}^* = 0.11 \div 0.28$), dải nhiệt độ từ nhiệt độ phòng ($T_0 = 300\text{ K}$) đến môi trường nhiệt độ cao ($500\text{ K} - 700\text{ K}$), cùng các cặp vật liệu FGM ($\text{Si}_3\text{N}_4/\text{SUS304}$) và FG-CNTRC trên hai loại nền polymer (PMMA và PmPV) kết hợp cốt gia cường SWCNT (10,10).

+-----------------------------------------------------------------------------------+
|                            KHUNG PHÂN TÍCH TỔNG QUAN                             |
+-----------------------------------------------------------------------------------+
|  VẬT LIỆU TIÊN TIẾN              MÔ HÌNH HÌNH HỌC              ĐIỀU KIỆN BIÊN     |
|  - Sandwich FGM (Loại A & B)     - Vỏ trụ tròn kín (Cylinder)  - Tựa đàn hồi      |
|  - FG-CNTRC (UD, FG-V, X, O)     - Vỏ trống thoải (Barrel)     - 0 <= c <= vô cực |
|  - Cơ tính phụ thuộc T (T-D)     - Nền Pasternak (k1, k2)      - Ràng buộc vĩ tuyến|
+-----------------------------------------------------------------------------------+
                                         │
                                         ▼
+-----------------------------------------------------------------------------------+
|                         CƠ SỞ LÝ THUYẾT & PHƯƠNG PHÁP                             |
|  - Động học phi tuyến von Kármán - Donnell (CST cho vỏ mỏng, FSDT cho vỏ dày)    |
|  - Hàm ứng suất Airy (f) & Hàm độ võng đa số hạng (w0, w1, w2)                    |
|  - Phương pháp tích phân Galerkin & Thuật toán lặp phi tuyến nhiệt đàn hồi       |
+-----------------------------------------------------------------------------------+
                                         │
                                         ▼
+-----------------------------------------------------------------------------------+
|                           KẾT QUẢ ĐỘT PHÁ & ĐÓNG GÓP                              |
|  - Xác lập bản chất chuyển tiếp: Rẽ nhánh (Bifurcation) <-> Cực trị (Snap-through)|
|  - Định lượng chính xác ảnh hưởng của độ cứng tiếp tuyến biên c lên tải tới hạn   |
|  - Cung cấp nghiệm giải tích chuẩn mực cho vỏ trống & vỏ dày FG-CNTRC            |
+-----------------------------------------------------------------------------------+

Literature Review và Positioning

Nghiên cứu về độ ổn định của kết cấu thành mỏng đã trải qua nhiều thập kỷ phát triển với các công trình nền tảng của Donnell, Timoshenko, và đặc biệt là chuyên khảo kinh điển của Brush và Almroth (1975) về ổn định của kết cấu thanh, tấm và vỏ. Khi vật liệu cơ tính biến đổi (FGM) được giới thiệu lần đầu tiên bởi các nhà khoa học vật liệu Nhật Bản (Koizumi, 1997), một làn sóng nghiên cứu ứng xử cơ học của tấm và vỏ FGM đã bùng nổ trên thế giới.

Trong dòng chảy nghiên cứu về ổn định của vỏ FGM, các nghiên cứu giải tích của Shen (2003, 2004, 2007) đã tiên phong sử dụng lý thuyết vỏ cổ điển và 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 bậc cao để khảo sát ứng xử sau vồng của vỏ trụ FGM chịu nén dọc trục, áp lực ngoài và nhiệt độ. Cùng thời kỳ, các công trình của Khazaeinejad và cộng sự (2010), Huang và Han (2008, 2010), Zenkour (2005) và Sofiyev (2010) đã phân tích ổn định tuyến tính và phi tuyến của tấm, vỏ trụ và vỏ nón FGM. Tại Việt Nam, nhóm nghiên cứu của Đào Huy Bích và cộng sự (2011, 2013), Đào Văn Dũng và cộng sự (2013, 2015), Trần Ích Thịnh và cộng sự (2014) đã có những đóng góp xuất sắc về ổn định phi tuyến của vỏ FGM có gân gia cường lệch tâm đặt trên nền đàn hồi.

Bước sang kỷ nguyên vật liệu nano, sau phát hiện mang tính lịch sử về ống nano carbon của Iijima (1991), Shen (2009) đã đề xuất ý tưởng mang tính cách mạng: phân bố CNTs theo quy luật hàm biến thiên qua chiều dày tấm polymer để tạo ra vật liệu FG-CNTRC. Tiếp nối đột phá này, Shen và Xiang (2012, 2014), Ansari và cộng sự (2014), Trang và Tùng (2016, 2018) đã phát triển các mô hình đánh giá ứng xử tĩnh, uốn và sau vồng của dầm, tấm và panel cong FG-CNTRC.

Tuy nhiên, tổng quan y văn quốc tế cho thấy hai cuộc tranh luận học thuật lớn và các khoảng trống tri thức then chốt:

  1. Tranh luận về ảnh hưởng của biến dạng trước vồng và hiện tượng hóp: Các nghiên cứu truyền thống thường bỏ qua biến dạng trước vồng (giả thiết màng thuần túy) để đơn giản hóa hệ phương trình. Tuy nhiên, Brush & Almroth (1975) và Shen (2003) nhấn mạnh rằng đối với các vỏ kín cong hai chiều, biến dạng trước vồng có độ lớn đáng kể. Bỏ qua yếu tố này sẽ làm biến mất điểm rẽ nhánh thực, dẫn đến đánh giá sai lệch nghiêm trọng về tải tới hạn và bỏ qua hiện tượng hóp tức thì (immediate snap-through) hoặc hóp chậm (delayed snap-through).
  2. Tranh luận về tính thực tế của điều kiện biên: Hầu hết các công bố quốc tế (như Shen, 2009; Kiani, 2014) đều giới hạn bài toán ở hai thái cực: tựa di động ($c = 0$) khi chịu tải cơ, hoặc tựa cố định ($c \to \infty$) khi chịu tải nhiệt. Sự phân chia cứng nhắc này không phản ánh đúng tương tác liên kết mềm đàn hồi trong thực tế kỹ thuật. Dù Librescu và cộng sự (2000), Tùng (2014) đã xem xét biên đàn hồi cho tấm và panel hở, nhưng chưa có bất kỳ công trình nào thiết lập giải pháp giải tích toàn diện cho vỏ kín tròn xoay (vỏ trụ và vỏ trống) làm từ FG-CNTRC và sandwich FGM có kể đến biên đàn hồi tiếp tuyến.

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

  • So với công trình của Shen (2014) về vỏ trụ FG-CNTRC chịu tải cơ: Shen chỉ khảo sát biên tựa di động tuyệt đối sử dụng phương pháp tiệm cận bậc cao với nghiệm một số hạng. Luận án của Phạm Thanh Hiếu mở rộng phạm vi bằng cách tích hợp liên kết biên đàn hồi tổng quát ($0 \le c \le \infty$), sử dụng hàm độ võng đa số hạng phản ánh trọn vẹn ba giai đoạn (trước vồng, vồng, sau vồng).
  • So với nghiên cứu của Zenkour (2005) về tấm sandwich FGM: Zenkour sử dụng nghiệm dạng Navier chỉ áp dụng cho kết cấu phẳng tựa bản lề bốn cạnh. Luận án này giải quyết vỏ kín tròn xoay có độ cong kép (vỏ trống), xử lý triệt để điều kiện khép kín chu vi vĩ tuyến $\int_0^L \int_0^{2\pi R} \frac{\partial v}{\partial y} dy dx = 0$.

Vị thế khoa học của luận án được xác lập như một bước tiến quan trọng, lấp đầy hoàn toàn khoảng trống lý thuyết về ổn định phi tuyến của vỏ kín tròn xoay nano composite chịu tương tác đa trường (nhiệt - cơ - biên đàn hồi).

Đó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 và làm phong phú thêm hệ thống lý thuyết ổn định phi tuyến của vỏ mỏng và vỏ dày trung bình thông qua việc phát triển các mô hình giải tích mới:

  • Mở rộng lý thuyết vỏ Donnell và CST: Luận án tích hợp thành công thế năng biến dạng của liên kết biên đàn hồi tiếp tuyến vào hệ phương trình phi tuyến von Kármán – Donnell. Bằng cách thiết lập mối quan hệ giữa hợp lực nén biên bị động $N_{x0}$ và chuyển vị rút ngắn đầu trục trung bình $\bar{\Delta x}$ thông qua độ cứng lò xo tiếp tuyến $c$ ($N_{x0} = c \cdot \bar{\Delta x}$), lý thuyết này xóa bỏ ranh giới cứng nhắc giữa biên di động và cố định, bao quát toàn bộ phổ điều kiện biên vật lý: $$\text{Điều kiện biên: } w = M_x = N_{xy} = 0, \quad N_x = N_{x0} = c \cdot \bar{\Delta x} \quad \text{tại } x = 0, L$$
  • Mô hình hóa lý thuyết vồng cho vỏ trống có độ cong Gauss: Mở rộng phương trình cân bằng cho vỏ thoải có hai bán kính cong ($R_x, R_y$), chứng minh một cách tường minh toán học vai trò của độ cong kinh tuyến đối với việc triệt tiêu sự suy giảm độ cứng sau vồng.
  • Xác lập hệ mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề 1 (P1): Khi độ cứng tiếp tuyến $c$ tăng từ 0 đến $\infty$, cơ chế mất ổn định nhiệt đàn hồi của vỏ trụ mỏng chuyển đổi trạng thái từ mất ổn định kiểu rẽ nhánh êm thuận sang mất ổn định cực trị kèm hiện tượng hóp tức thì với cường độ hóp tỷ lệ thuận với $c$.
    • Mệnh đề 2 (P2): Đối với vỏ FG-CNTRC, hàm phân bố CNTs bất đối xứng qua mặt trung hòa tạo ra sự liên kết uốn - kéo (bending-stretching coupling), dẫn đến sự xuất hiện của độ võng ban đầu ngay khi bắt đầu gia tải nhiệt hoặc nén, làm suy giảm tải tới hạn rẽ nhánh cổ điển.
    • Mệnh đề 3 (P3): Trong phân tích vỏ dày FG-CNTRC sử dụng FSDT, việc sử dụng dạng nghiệm hàm võng hai số hạng ($m, n$) giải quyết triệt để tính bất định của hệ phương trình đại số phi tuyến mà dạng nghiệm ba số hạng gặp phải khi thỏa mãn đồng thời điều kiện biên trượt ngang và biến dạng màng kín.

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

Khung phân tích của luận án được xây dựng dựa trên sự liên kết chặt chẽ của ba trụ cột lý thuyết:

  1. Cơ học vật liệu tiên tiến:
    • Mô hình FGM sandwich loại A (lõi kim loại SUS304, hai lớp mặt FGM giàu gốm $\text{Si}_3\text{N}_4$) và loại B (mặt ngoài gốm, lõi FGM, mặt trong kim loại) với quy luật phân bố hàm mũ $V_M(z) = [ (2z+h)/2h_f ]^N$.
    • Mô hình FG-CNTRC với 4 kiểu phân bố CNTs qua chiều dày (UD, FG-V, FG-O, FG-X), xác định cơ tính hiệu dụng qua quy tắc hỗn hợp mở rộng kết hợp động lực học phân tử: $$E_{11} = \eta_1 V_{CNT} E_{11}^{CNT} + V_m E_m, \quad \frac{\eta_2}{E_{22}} = \frac{V_{CNT}}{E_{22}^{CNT}} + \frac{V_m}{E_m}, \quad \frac{\eta_3}{G_{12}} = \frac{V_{CNT}}{G_{12}^{CNT}} + \frac{V_m}{G_m}$$
  2. Động học phi tuyến và quan hệ cấu tạo: Thiết lập phương trình tương thích biến dạng màng và phương trình cân bằng phi tuyến thông qua hàm ứng suất Airy $f(x,y)$ và hàm độ võng $w(x,y)$: $$\nabla_m^4 f = E_1 \left( w_{,xy}^2 - w_{,xx}w_{,yy} - \frac{w_{,xx}}{R_y} - \frac{w_{,yy}}{R_x} \right)$$ $$D \nabla^4 w + f_{,yy}w_{,xx} - 2f_{,xy}w_{,xy} + f_{,xx}w_{,yy} + \frac{f_{,xx}}{R_y} + \frac{f_{,yy}}{R_x} - q = 0$$
  3. Điều kiện biên khép kín chu vi vĩ tuyến: Tích hợp điều kiện biến dạng liên tục theo phương vòng $\int_0^L \int_0^{2\pi R} \frac{\partial v}{\partial y} dy dx = 0$, đảm bảo tính đóng hình học đặc thù của vỏ tròn xoay so với tấm và panel hở.

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận của chủ nghĩa thực chứng hậu nghiệm (Post-positivism) kết hợp cách tiếp cận cơ học giải tích tất định (deterministic analytical mechanics). Thiết kế nghiên cứu được cấu trúc đa tầng:

[Cấp độ Vi mô (Micro)]       [Cấp độ Cấu trúc (Meso)]        [Cấp độ Kết cấu (Macro)]
Đặc tính CNT & Matrix   ───►  Phân bố FG-X, UD qua h   ───►  Vỏ trụ / Vỏ trống (CST/FSDT)
Quy tắc hỗn hợp + MD         Tính chất hiệu dụng (T-D)       Phương trình Galerkin & Lặp
  • Mức độ vi mô: Tính toán cơ tính của SWCNT (10,10) ($L_t = 9.26\text{ nm}$, $R_t = 0.678\text{ nm}$, $E_{11}^{CNT} = 5.646\text{ TPa}$, $\sigma_b = 63\text{ GPa}$, $\rho_{CNT} = 1300\text{ kg/m}^3$) và nền polymer (PMMA có $E_m = 2.5\text{ GPa}$, $\nu_m = 0.34$; PmPV có $E_m = 2.1\text{ GPa}$, $\nu_m = 0.4$).
  • Mức độ cấu trúc: Tích phân cơ tính qua chiều dày $h$ để xác định các ma trận độ cứng màng ($A_{ij}$), độ cứng uốn ($D_{ij}$), độ cứng liên kết uốn-kéo ($B_{ij}$) và các số hạng nhiệt đàn hồi $(\Phi_0, \Phi_1)$.
  • Mức độ kết cấu: Phân tích ứng xử toàn cục của vỏ kín dưới tải trọng cơ (nén dọc trục $P$, áp lực ngoài $q$) và tải nhiệt tăng đều $\Delta T = T - T_0$.

Quy trình nghiên cứu rigorous

Quy trình giải tích bán số được thực hiện qua các bước chuẩn hóa nghiêm ngặt:

  1. Thiết lập hàm độ võng ba số hạng (Three-term deflection function): $$w(x, y) = W_0 + W_1 \sin(\alpha x) \sin(\beta y) + W_2 \sin^2(\alpha x)$$ Trong đó, $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 với số sóng theo phương trục $m$ và phương vòng $n$ ($\alpha = m\pi/L, \beta = n/R$); $W_2$ là biên độ của số hạng đối xứng trục phi tuyến giai đoạn sau vồng.
  2. Xác định hàm ứng suất Airy $f(x, y)$: Giải chính xác phương trình tương thích biến dạng vi phân bậc 4 để tìm biểu thức giải tích đóng của $f(x, y)$ chứa các số hạng liên kết cơ - nhiệt và tham số biên đàn hồi $c$.
  3. Áp dụng trực giao hóa Galerkin: Chiếu phương trình cân bằng phi tuyến lên không gian hàm thử: $$\int_0^L \int_0^{2\pi R} \left[ \text{Cân bằng phi tuyến} \right] \sin(\alpha x) \sin(\beta y) , dy dx = 0$$ $$\int_0^L \int_0^{2\pi R} \left[ \text{Cân bằng phi tuyến} \right] \sin^2(\alpha x) , dy dx = 0$$
  4. Thuật toán lặp phi tuyến xử lý cơ tính phụ thuộc nhiệt độ (T-D Iterative Algorithm): Vì các mô đun $E(T)$ và hệ số dãn nở $\alpha(T)$ thay đổi liên tục theo nhiệt độ $T$, thuật toán lặp được thiết lập:
    • Gán $T^{(0)} = T_0 = 300\text{ K}$. Tính toán các hệ số vật liệu.
    • Giải phương trình Galerkin tìm tải nhiệt tới hạn $T_{cr}^{(k)}$ hoặc nghiệm sau vồng.
    • Cập nhật cơ tính vật liệu tại $T^{(k)}$. Kiểm tra điều kiện hội tụ $|T^{(k)} - T^{(k-1)}| / T^{(k)} < 10^{-4}$.
    • Nếu chưa hội tụ, lặp lại quy trình cho đến khi đạt độ chính xác tuyệt đối.
                    ┌────────────────────────────────────────┐
                    │ Bắt đầu: Khởi tạo T^(0) = T0 = 300 K  │
                    └───────────────────┬────────────────────┘
                                        │
                                        ▼
                    ┌────────────────────────────────────────┐
                    │ Tính toán cơ tính vật liệu E(T), a(T)  │
                    └───────────────────┬────────────────────┘
                                        │
                                        ▼
                    ┌────────────────────────────────────────┐
                    │ Giải hệ Galerkin tìm tải tới hạn T_cr  │
                    └───────────────────┬────────────────────┘
                                        │
                                        ▼
                           /─────────────────────────\
                          /   |T^(k) - T^(k-1)|/T^(k) \  Chưa hội tụ
                         <            < 10^-4          >───────────────┐
                          \                           /                │
                           \─────────────────────────/                 │
                                        │ Đã hội tụ                    │
                                        ▼                              ▼
                    ┌─────────────────────────┐     ┌───────────────────────────────────┐
                    │ Xuất kết quả tải tới    │     │ Cập nhật cơ tính vật liệu theo    │
                    │ hạn và đường sau vồng   │     │ nhiệt độ mới T^(k) và lặp lại     │
                    └─────────────────────────┘     └───────────────────────────────────┘

Data và phân tích

Tính toán số trị và kiểm chứng kết quả được lập trình tự động hóa trên môi trường tính toán ký hiệu và số học cấp cao (MATLAB/Maple):

  • Kiểm định độ tin cậy (Validation & Triangulation):
    • So sánh tải tới hạn nén dọc trục $F_{cr} = 2\pi R h P_{cr}\text{ (kN)}$ của vỏ trụ FG-CNTRC ($R/h = 100, h = 1\text{ mm}, T = 300\text{ K}$) với kết quả tiệm cận của Shen (2014) cho thấy sai số nhỏ hơn $0.5%$.
    • So sánh áp lực tới hạn $q_{cr}\text{ (MPa)}$ của vỏ trụ FGM $\text{Si}_3\text{N}_4/\text{SUS304}$ ($R/h = 300, T = 300\text{ K}, \Delta T = 200\text{ K}$) với công bố của Huang và Han (2010), độ lệch tuyệt đối nằm trong khoảng $0.12% - 0.85%$, chứng minh tính hội tụ và độ chuẩn xác vượt trội của khung giải tích.

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

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

Các phân tích số trị sâu rộng đã mang lại 5 phát hiện mang tính đột phá về cơ học phi tuyến của vỏ tròn xoay:

  1. Quy luật chi phối của độ cứng ràng buộc tiếp tuyến biên ($c$):
    • Khi độ cứng tiếp tuyến tăng từ trạng thái tự do di động ($c = 0$) lên liên kết đàn hồi ($c = 10^7\text{ N/m}$) và ngàm cố định ($c \to \infty$), áp lực ngoài tới hạn $q_{cr}$ của vỏ trụ sandwich FGM loại A ($R/h = 100, L/R = 2, T = 300\text{ K}$) tăng vọt từ $0.185\text{ MPa}$ lên hơn $0.412\text{ MPa}$ (tăng trên $122%$).
    • Tuy nhiên, độ cứng biên $c$ càng lớn thì ứng xử sau vồng nhiệt đàn hồi càng trở nên bất ổn định dữ dội; vỏ xuất hiện hiện tượng hóp tức thì (immediate snap-through) với bước nhảy độ võng lớn, giải phóng năng lượng biến dạng tích lũy đột ngột.
  2. Tính ưu việt tuyệt đối của kiểu phân bố FG-X:
    • Trong tất cả các dạng phân bố, kiểu phân bố FG-X (CNTs tập trung tối đa ở hai bề mặt ngoài $z = \pm h/2$) luôn mang lại tải tới hạn nén $P_{cr}$ và áp lực ngoài $q_{cr}$ cao nhất. Cụ thể, với cùng tỷ lệ thể tích $V_{CNT}^* = 0.17$, vỏ trụ FG-X có tải tới hạn nén cao hơn $48.6%$ so với phân bố đều (UD) và cao hơn $115.4%$ so với phân bố FG-O (CNTs tập trung ở lớp giữa). Điều này chứng minh quy luật phân bố tối ưu là đưa vật liệu chịu lực cao ra xa mặt trung hòa để cực đại hóa mômen quán tính kháng uốn.
  3. Hiệu ứng gia cường của độ cong Gauss trên vỏ trống:
    • Vỏ trống với tỷ số bán kính cong kinh tuyến trên bán kính vĩ tuyến ($R_x/R_y > 0$) sở hữu khả năng chịu tải cơ kết hợp (nén dọc trục + áp lực ngoài) vượt trội từ $25%$ đến $60%$ so với vỏ trụ tròn thẳng tương đương ($R_x \to \infty$). Độ cong Gauss dương đóng vai trò như một cơ chế gân tự thân hình học (geometric self-stiffening), hạn chế đáng kể xu hướng suy giảm độ cứng sau vồng.
  4. Tác động suy giảm nghiêm trọng của cơ tính phụ thuộc nhiệt độ (T-D):
    • Sự phụ thuộc nhiệt độ của vật liệu làm suy giảm mạnh tải nhiệt tới hạn $\Delta T_{cr}$. Khi tăng nhiệt độ từ $300\text{ K}$ lên $500\text{ K}$, áp lực ngoài tới hạn của vỏ sandwich FGM loại B giảm từ $28%$ đến $42%$ so với khi giả thiết cơ tính không đổi (T-ID). Việc tính toán dựa trên mô hình T-ID sẽ tạo ra sự mất an toàn nghiêm trọng trong thiết kế chịu nhiệt.
  5. Ảnh hưởng tiêu tán của nền đàn hồi Pasternak:
    • Nền Pasternak với độ cứng lớp trượt $K_2$ có hiệu quả kháng vồng và ổn định hóa sau vồng cao gấp $3 - 5$ lần so với mô hình nền Winkler đơn thuần ($K_1$), triệt tiêu hiệu quả các mode vồng bậc cao và ngăn chặn hiện tượng hóp chậm của vỏ trống.
+-----------------------------------------------------------------------------------+
|               SO SÁNH TẢI TỚI HẠN VÀ KHẢ NĂNG KHÁNG VỒNG THEO KIỂU PHÂN BỐ        |
+-----------------------------------------------------------------------------------+
| Kiểu phân bố CNTs  | Vị trí tập trung CNTs | Tải nén P_cr (MPa) | Hiệu quả cơ học |
+--------------------+-----------------------+--------------------+-----------------+
| FG-X               | Hai bề mặt ngoài      | Cao nhất (Max)     | Xuất sắc (100%) |
| UD                 | Phân bố đều           | Trung bình         | Khá (~67%)      |
| FG-V / FG-Lambda   | Một bề mặt (trên/dưới)| Trung bình - Thấp  | Trung bình (~55%)|
| FG-O               | Mặt trung hòa giữa    | Thấp nhất (Min)    | Kém nhất (~46%) |
+--------------------+-----------------------+--------------------+-----------------+

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp một mô hình giải tích tổng quát, khép kín, có khả năng dự báo chính xác hành vi rẽ nhánh và sau vồng của các kết cấu vỏ tròn xoay nano composite mà không cần tiêu tốn tài nguyên tính toán khổng lồ của mô phỏng phần tử hữu hạn 3D phi tuyến.
  • Về mặt phương pháp luận: Chuẩn hóa kỹ thuật nghiệm hai số hạng cho FSDT trên vỏ kín, giải quyết xung đột toán học giữa điều kiện biên trượt và tính tuần hoàn chu vi, mở đường cho các phân tích dao động phi tuyến và động lực học va chạm phức tạp hơn.
  • Về mặt ứng dụng thực tiễn: Định hình trực tiếp quy trình thiết kế bồn chứa nhiên liệu đông lạnh, vỏ thân máy bay siêu âm, ống xả tên lửa và vỏ tàu ngầm chịu áp lực sâu; chỉ rõ cách thiết kế gối đỡ biên đàn hồi có kiểm soát để tránh vùng sụp đổ hóp tức thì.

Limitations và Future Research

Mặc dù đạt được những bước tiến vượt bậc, luận án ghi nhận 4 giới hạn khoa học cần được mở rộng:

  1. Giả thiết vật liệu hoàn hảo: Mô hình giả thiết các ống nano CNTs phân bố hoàn toàn thẳng, liên kết hoàn hảo với nền polymer và không xảy ra hiện tượng vón cục (agglomeration) hay rạn nứt tế vi.
  2. Giả thiết liên kết lớp sandwich: Chưa xét đến hiện tượng bong tách lớp (delamination) do ứng suất cắt tiếp xúc giữa lớp mặt FGM và lõi kim loại dưới tải nhiệt chu kỳ.
  3. Phạm vi hình học của vỏ trống: Giới hạn ở vỏ trống rất thoải ($R_x \gg R_y$), chưa áp dụng trực tiếp cho các vỏ trống có độ cong kinh tuyến sâu hoặc vỏ bán cầu.
  4. Trạng thái biến dạng đàn hồi tuyến tính của vật liệu: Bỏ qua biến dạng dẻo và hiện tượng từ biến (creep) của kim loại ở nhiệt độ siêu cao ($T > 1000\text{ K}$).

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

  • Phát triển mô hình ổn định động lực học phi tuyến và đáp ứng chuyển tiếp (transient dynamic response) của vỏ FG-CNTRC dưới tác động của sóng xung kích nổ và áp lực nhiệt biến thiên theo thời gian.
  • Tích hợp lý thuyết phá hủy liên tục (Damage Mechanics) và hiện tượng vón cục CNTs theo mô hình Eshelby-Mori-Tanaka hai pha.
  • Nghiên cứu tối ưu hóa đa mục tiêu cấu trúc vi mô (topology optimization) nhằm tìm kiếm hàm phân bố CNTs phi tuyến 3 chiều tùy biến theo trường ứng suất cục bộ.
  • Triển khai thực nghiệm chế tạo mẫu in 3D nano composite và đo đạc trường biến dạng thực tế bằng công nghệ quang học tương quan hình ảnh kỹ thuật số (DIC) để kiểm chứng nghiệm giải tích.

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

  • Ảnh hưởng học thuật: Các công trình công bố từ luận án trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus (như Composite Structures, Thin-Walled Structures, Acta Mechanica) tạo nền tảng trích dẫn vững chắc cho cộng đồng cơ học tính toán và vật liệu composite tiên tiến toàn cầu.
  • Chuyển đổi công nghiệp & hàng không vũ trụ: Cung cấp bộ công cụ tính toán giải tích trực tiếp cho các kỹ sư R&D trong ngành hàng không vũ trụ, giúp giảm trọng lượng kết cấu vỏ tên lửa từ $20% - 35%$ nhờ tối ưu hóa vật liệu FG-CNTRC thay thế hợp kim nhôm - titan truyền thống mà vẫn đảm bảo hệ số an toàn ổn định tĩnh.
  • Tác động xã hội và môi trường: Ứng dụng vật liệu composite gia cường nano giúp nâng cao hiệu suất sử dụng nhiên liệu của các phương tiện giao thông vận tải, giảm phát thải khí nhà kính và tăng cường độ bền vận hành an toàn 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 & Nhà nghiên cứu trẻ: Tiếp cận phương pháp giải tích mẫu mực, cách xử lý toán học tinh tế đối với hàm ứng suất Airy, phương pháp Galerkin và thuật toán lặp nhiệt phi tuyến trên vỏ kín.
  • Các giáo sư & Chuyên gia cơ học kết cấu: Sử dụng các nghiệm giải tích đóng làm chuẩn kiểm định (benchmark solutions) chính xác tuyệt đối để đánh giá độ tin cậy của các phần tử hữu hạn mới hoặc các thuật toán không lưới (mesh-free methods).
  • Kỹ sư thiết kế R&D công nghiệp: Ứng dụng trực tiếp các biểu thức toán học và biểu đồ thiết kế để tra cứu tải tới hạn, tối ưu hóa kích thước vỏ sandwich FGM và lựa chọn tỷ lệ thể tích CNTs phù hợp với điều kiện biên của dự án thực tế.
  • Cơ quan quản lý tiêu chuẩn kỹ thuật: Bổ sung cơ sở khoa học để xây dựng và hiệu chỉnh các tiêu chuẩn thiết kế bình bồn chịu áp lực, ống dẫn công nghiệp làm việc trong môi trường nhiệt độ 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ì?

Đó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 thống nhất giải bài toán ổn định phi tuyến của vỏ kín tròn xoay FGM và FG-CNTRC chịu liên kết biên đàn hồi tiếp tuyến tổng quát ($0 \le c \le \infty$). Luận án đã phá vỡ giả thiết nhị phân truyền thống (chỉ xét biên di động tự do hoặc cố định tuyệt đối), chứng minh toán học tường minh cơ chế chuyển pha mất ổn định từ rẽ nhánh êm sang hóp tức thì dưới tác động của độ cứng lò xo biên $c$.

2. Đột phá phương pháp luận so với các nghiên cứu trước đây được thể hiện như thế nào?

So với nghiên cứu của Shen (2014) và Khazaeinejad et al. (2010), đột phá phương pháp luận thể hiện ở hai điểm:

  • Thiết lập hàm độ võng ba số hạng đa mode ($W_0, W_1, W_2$) cho CST, giải quyết đồng thời cả 3 giai đoạn: biến dạng trước vồng, điểm rẽ nhánh và ứng xử sau vồng phi tuyến lớn.
  • Đề xuất nghiệm hàm võng hai số hạng cải tiến trong khuôn khổ FSDT cho vỏ dày trung bình, khắc phục triệt để bế tắc không thể khép kín nghiệm của mô hình ba số hạng khi thỏa mãn điều kiện biên trượt ngang và biến dạng màng khép kín chu vi vĩ tuyến.

3. Phát hiện số trị nào mang tính bất ngờ nhất và bản chất vật lý của nó là gì?

Phát hiện bất ngờ nhất là hiện tượng hóp nhiệt đàn hồi nghịch đảo khi tăng độ cứng biên đàn hồi $c$. Mặc dù $c$ tăng giúp nâng cao đáng kể tải áp lực tới hạn ban đầu, nhưng nó lại biến quỹ đạo sau vồng thành một đường cong suy giảm dốc đứng (unstable postbuckling path), gây ra hiện tượng sụp đổ kết cấu đột ngột không báo trước ở mức tải thấp hơn nhiều so với tải tới hạn rẽ nhánh lý thuyết. Bản chất vật lý là do năng lượng biến dạng nén màng bị giam hãm dọc trục khi biên bị cản trở, giải phóng tức thì thành năng lượng uốn phi tuyến lớn khi vỏ bị lệch khỏi trạng thái cân bằng cơ bản.

Tải P 
  ▲         Đường cân bằng cơ bản (Primary path)
  │              /
P_cr ───────────●  (Điểm rẽ nhánh - Bifurcation point)
  │            / \
  │           /   \  Quỹ đạo sau vồng có hóp tức thì (Snap-through với c -> vô cực)
  │          /     ▼
  │         /       ● (Điểm cực tiểu - Minimum postbuckling load)
  │        /         \
  │       /           \___  Quỹ đạo sau vồng ổn định (với c = 0, tựa di động)
  └──────┴───────────────────► Độ võng W_max

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

Có. Luận án công khai toàn bộ hệ thống phương trình vi phân cơ bản, định nghĩa các ma trận độ cứng tích phân qua chiều dày, công thức giải tích đóng của các hệ số Galerkin, bảng hệ số nhiệt phụ thuộc của vật liệu $\text{Si}_3\text{N}_4$, $\text{SUS304}$, PMMA, PmPV, SWCNT (10,10) và lưu đồ thuật toán lặp nhiệt - cơ. 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ố trị và đường cong cân bằng được trình bày trong luận án.

5. Tầm nhìn chương trình nghiên cứu 10 năm (10-Year Research Agenda) mở ra từ luận án là gì?

Chương trình nghiên cứu dài hạn bao gồm:

  • Mở rộng lý thuyết sang các cấu trúc vỏ nanocomposite thông minh tích hợp vật liệu áp điện (piezoelectric) và từ - nhiệt - đàn hồi (magneto-electro-elastic) để kiểm soát chủ động hiện tượng vồng.
  • Nghiên cứu ổn định động phi tuyến đa trường (Aero-thermo-elastic flutter) của vỏ FG-CNTRC dưới dòng khí siêu âm (hypersonic flow).
  • Phát triển phần mềm tính toán chuyên dụng mã nguồn mở phục vụ thiết kế vỏ mỏng nano composite cho công nghiệp hàng không vũ trụ và năng lượng hydro lỏng.

Kết luận

Luận án tiến sĩ của tác giả Phạm Thanh Hiếu là một công trình nghiên cứu khoa học xuất sắc, mẫu mực về cơ học kết cấu và vật liệu tiên tiến, kết tinh những giá trị học thuật đỉnh cao:

  1. Thiết lập thành công khung lý thuyết giải tích toàn diện giải quyết bài toán ổn định tuyến tính và phi tuyến của vỏ trụ và vỏ trống sandwich FGM và FG-CNTRC có kể đến tính đàn hồi của liên kết biên tiếp tuyến.
  2. Chứng minh và định lượng chính xác quy luật chi phối của độ cứng biên $c$, làm sáng tỏ bản chất chuyển tiếp giữa hiện tượng mất ổn định rẽ nhánh và mất ổn định cực trị/hóp trong vỏ tròn xoay kín.
  3. Khẳng định tính vượt trội của cấu hình phân bố FG-X và độ cong Gauss dương, cung cấp nguyên lý định hướng tối ưu hóa vật liệu và hình học cho các cấu trúc vỏ chịu tải nén và áp lực ngoài khắc nghiệt.
  4. Chuẩn hóa phương pháp giải tích FSDT với nghiệm hai số hạng cải tiến cho vỏ dày trung bình, vượt qua rào cản tính toán tồn tại nhiều năm trong y văn cơ học vỏ kín.
  5. Đánh giá chuẩn xác ảnh hưởng suy giảm của tính chất vật liệu phụ thuộc nhiệt độ, thiết lập quy chuẩn tính toán an toàn cho kết cấu làm việc trong môi trường nhiệt độ cao.
  6. Mở ra 3 nhánh nghiên cứu mới: Ổn định động lực học vỏ nano composite có liên kết biên đàn hồi; Tối ưu hóa vi cấu trúc CNT 3D thích ứng; và Cơ học phá hủy nhiệt đàn hồi trên vỏ mỏng thông minh.

Công trình không chỉ đóng góp to lớn vào kho tàng tri thức cơ học công trình trong nước và quốc tế mà còn khẳng định năng lực nghiên cứu đỉnh cao, độc lập và sáng tạo của các nhà khoa học Việt Nam trên bản đồ học thuật thế giới.