Tổng quan về luận án

Sự phát triển đột phá của khoa học vật liệu tiên tiến trong những thập kỷ gần đây đã thúc đẩy sự ra đời của các hệ kết cấu tấm và vỏ mỏng nanocomposite có khả năng làm việc trong môi trường cơ học và nhiệt học khắc nghiệt. Luận án tiến sĩ chuyên ngành Cơ kỹ thuật (Mã số: 952 01 01) của nghiên cứu sinh Lê Thị Như Trang với đề tài "Ổn định nhiệt đàn hồi của tấm và vỏ thoải composite gia cường carbon nanotube", dưới sự hướng dẫn khoa học của PGS.TS Hoàng Văn Tùng và PGS.TS Nguyễn Đình Kiên 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, đại diện cho một công trình nghiên cứu hàn lâm tiên phong, giải quyết toàn diện bài toán ứng xử phi tuyến và mất ổn định nhiệt đàn hồi của kết cấu nanocomposite thế hệ mới.

                    +-------------------------------------------------------------+
                    |  ỐNG NANO CARBON (SWCNT 10,10: E > 1 TPa, Độ bền = 63 GPa)  |
                    +-------------------------------------------------------------+
                                                   |
                                                   v
                    +-------------------------------------------------------------+
                    |  MÔ HÌNH VẬT LIỆU NANOCOMPOSITE CƠ TÍNH BIẾN ĐỔI (FG-CNTRC) |
                    |  Phân bố: UD, FG-V, FG-O, FG-X, FG-Λ | Quy tắc hỗn hợp suy rộng|
                    +-------------------------------------------------------------+
                                                   |
                                                   v
                    +-------------------------------------------------------------+
                    |  KHUNG LÝ THUYẾT ĐA CẤP (CPT / FSDT / HSDT + von Kármán)    |
                    |  Biên đàn hồi tiếp tuyến (k1, k2) + Nền Pasternak (K1, K2)   |
                    +-------------------------------------------------------------+
                                                   |
                         +-------------------------+-------------------------+
                         |                                                   |
                         v                                                   v
+--------------------------------------------------+ +--------------------------------------------------+
|          MẤT ỔN ĐỊNH CƠ HỌC & NHIỆT ĐÀN HỒI       | |            QUÁ TRÌNH ỨNG XỬ SAU VỒNG             |
| - Mất ổn định kiểu rẽ nhánh (Bifurcation)        | | - Trạng thái cân bằng thứ cấp sau vồng          |
| - Mất ổn định kiểu cực trị / giới hạn (Limit)    | | - Hiện tượng hóp kết cấu (Snap-through)         |
| - Hiện tượng mất ổn định ngay từ đầu             | | - Đáp ứng tải - độ võng phi tuyến giải tích     |
+--------------------------------------------------+ +--------------------------------------------------+

Bối cảnh khoa học của nghiên cứu khởi nguồn từ việc khai thác các đặc tính cơ lý siêu việt của ống nano carbon (CNTs). Kể từ công trình bước ngoặt năm 1991 của Iijima, CNTs đã được chứng minh là vật liệu gia cường lý tưởng nhờ sở hữu khối lượng riêng thấp xấp xỉ $1300\text{ kg/m}^3$, mô đun đàn hồi vượt trên $1\text{ TPa}$ và "độ bền chịu kéo cao nhất đã đo được đối với một CNT là 63 GPa, tức là cao gấp khoảng 10-100 lần so với độ bền chịu kéo của thép cường độ cao" [1, 3-5]. Kế thừa ý tưởng vật liệu có cơ tính biến đổi (Functionally Graded Materials - FGM), Shen (2009) [13] đã đề xuất mô hình composite gia cường ống nano carbon có cơ tính biến thiên (FG-CNTRC), trong đó tỷ lệ thể tích CNTs được tối ưu hóa qua chiều dày theo các hàm toán học tuyến tính.

Tuy nhiên, một khoảng trống nghiên cứu (research gap) cốt lõi tồn tại trong y văn quốc tế là phần lớn các công trình trước đây chỉ xem xét các điều kiện biên lý tưởng hóa (hoàn toàn tự do dịch chuyển hoặc hoàn toàn ngàm cứng tuyệt đối) và giả định trường nhiệt đồng nhất. Thực tế kỹ thuật đòi hỏi phải giải quyết bài toán với điều kiện biên liên kết đàn hồi nửa cứng theo phương tiếp tuyến (tangential elastic restraint), trường nhiệt độ biến thiên không đều trong mặt phẳng (dạng sin, dạng tuyến tính), và sự tương tác phức tạp của các trường ứng suất trước cơ - nhiệt. Luận án đã thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Mức độ ràng buộc đàn hồi dịch chuyển tiếp tuyến ở biên chi phối như thế nào đến sự chuyển dịch giữa các dạng mất ổn định (rẽ nhánh, cực trị, hay không tồn tại điểm rẽ nhánh) của tấm và panel vỏ thoải FG-CNTRC?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật phân bố CNTs (UD, FG-V, FG-O, FG-X, $\text{FG-}\Lambda$) và tỷ lệ thể tích tổng $V^*_{CNT}$ tác động định lượng ra sao đến khả năng chịu tải tới hạn và đường cân bằng sau vồng nhiệt đàn hồi?
  • Câu hỏi nghiên cứu 3 (RQ3): Sự phụ thuộc nhiệt độ của vật liệu (Temperature-Dependent - T-D) kết hợp với các trường nhiệt không đều và tải nén/áp lực ngoài trước ảnh hưởng như thế nào đến độ võng và độ ổn định tĩnh phi tuyến?
  • Giả thuyết 1 (H1): Cấu hình phân bố FG-X mang lại độ cứng uốn và độ ổn định nhiệt đàn hồi vượt trội nhất do tập trung mật độ CNTs cao tại các bề mặt biên chịu ứng suất uốn lớn nhất.
  • Giả thuyết 2 (H2): Sự gia tăng độ cứng ràng buộc biên tiếp tuyến sẽ tăng cường khả năng chịu tải cơ học nhưng làm suy giảm nhiệt độ tới hạn rẽ nhánh do thúc đẩy ứng suất nén nhiệt nội tại phát triển nhanh hơn.
  • Giả thuyết 3 (H3): Tồn tại các ngưỡng độ cong và ứng suất cơ - nhiệt trước khiến panel vỏ thoải mất hoàn toàn trạng thái cân bằng màng ban đầu, dẫn đến mất ổn định ngay từ thời điểm bắt đầu gia tải.

Khung lý thuyết của luận án tích hợp chặt chẽ: Lý thuyết tấm vỏ cổ điển (CPT), Lý thuyết biến dạng trượt bậc nhất (FSDT - Reissner-Mindlin), Lý thuyết biến dạng trượt bậc cao (HSDT - Reddy), Lý thuyết phi tuyến hình học von Kármán, Thuyết vỏ thoải Donnell, và Mô hình nền đàn hồi hai tham số Pasternak. Luận án khảo sát phạm vi toàn diện trên các mẫu vật liệu nanocomposite nền polymer (PMMA, PmPV) gia cường ống nano carbon đơn vách SWCNT (10,10) với các tỷ lệ thể tích $V^*_{CNT} \in {0.11, 0.14, 0.17, 0.28}$, dải nhiệt độ $T \in [300\text{ K}, 500\text{ K}]$, trên các dạng kết cấu tấm chữ nhật, panel trụ mỏng/dày và panel hai độ cong.

Literature Review và Positioning

Nghiên cứu về ổn định kết cấu composite tiên tiến đã trải qua các bước tiến hóa quan trọng trong ba thập kỷ qua. Bắt đầu từ các công trình nền tảng về tính chất nano của Thostenson et al. (2001, 2005) [3, 4] và Coleman et al. (2006) [5], tiềm năng gia cường của CNTs đã được chứng minh với nhận định kinh điển: "chỉ cần bổ sung 1% khối lượng (khoảng 0.5% thể tích) CNTs vào nền polystyrene có thể làm cho mô đun đàn hồi của composite tăng khoảng 36 - 42% và độ bền chịu kéo của composite tăng khoảng 25%". Sau khi Shen (2009) [13] phát triển mô hình FG-CNTRC cho phân tích uốn phi tuyến tấm, hàng loạt các hướng nghiên cứu quốc tế đã mở rộng sang phân tích tĩnh, dao động và ổn định tuyến tính.

Luồng nghiên cứu chính Tác giả tiêu biểu & Năm Phương pháp tiếp cận Giới hạn nghiên cứu chính
Dao động & Phân tích tĩnh 3D/Quasi-3D Alibeigloo et al. (2011) [14-16]; Pourasghar & Chen (2016) [17]; Fazzolari (2014) [25] Giải tích 3D Elasticity, DQM, Quasi-3D Ritz Giới hạn ở tải cơ học vuông góc hoặc dao động tuyến tính, chưa khảo sát mất ổn định phi tuyến sau vồng.
Phần tử hữu hạn & Mesh-free phi tuyến Liew, Zhang et al. (2014-2015) [20-22, 58, 59]; Tornabene et al. (2016) [23, 24] IMLS-Ritz, Mesh-free kp-Ritz, GDQ Chi phí tính toán cao, chủ yếu áp dụng biên lý tưởng, khó trích xuất biểu thức giải tích dạng đóng.
Ổn định nhiệt đàn hồi tuyến tính Kiani, Mirzaei et al. (2013-2015) [41-43, 92-94]; Torabi et al. (2013) [96] Chebyshev-Ritz, DQM, DSC Chỉ xác định tải nhiệt tới hạn tại điểm rẽ nhánh, bỏ qua lộ trình cân bằng phi tuyến sau vồng.
Ổn định phi tuyến tấm vỏ FG-CNTRC Shen & Zhang (2012) [100]; Shen & Xiang (2014) [103]; Tùng (2016) [102] Phương pháp tiệm cận hai bước, Giải tích Galerkin Chưa khảo sát đồng thời biên đàn hồi tiếp tuyến cho panel cong và trường nhiệt phân bố không đều trong mặt phẳng.

Trong y văn học thuật, tồn tại một cuộc tranh luận sâu sắc giữa hai quan điểm đối nghịch về cơ chế mất ổn định nhiệt đàn hồi của panel vỏ cong:

  1. Quan điểm cổ điển truyền thống: Giả định rằng panel vỏ hoàn hảo chịu tải nhiệt đồng nhất luôn trải qua trạng thái màng trước vồng ổn định và mất ổn định thông qua một điểm rẽ nhánh (bifurcation buckling) xác định, tương tự như tấm phẳng chịu nén.
  2. Quan điểm hiện đại phi tuyến: Cho rằng do độ cong hình học ban đầu kết hợp với tính bất đối xứng cấu trúc và điều kiện ràng buộc biên, trạng thái màng bị phá vỡ ngay khi gia nhiệt, khiến kết cấu bị uốn võng tức thì theo hướng lồi hoặc lõm mà không hề tồn tại điểm rẽ nhánh thực sự.

Luận án định vị chính xác vị trí khoa học của mình bằng cách hòa giải và làm sáng tỏ cuộc tranh luận này. So sánh trực tiếp với hai nghiên cứu quốc tế cốt lõi:

  • So với công trình của Shen & Zhang (2012) [100] (nghiên cứu tấm FG-CNTRC chịu nhiệt đều và truyền nhiệt phi tuyến qua chiều dày bằng phương pháp tiệm cận với biên ngàm cố định): Luận án mở rộng toàn diện sang điều kiện biên đàn hồi tiếp tuyến tổng quát ($k_1, k_2$), đồng thời giải quyết bài toán truyền nhiệt dạng hàm sin và hàm tuyến tính trong mặt phẳng kết cấu.
  • So với nghiên cứu của Librescu et al. (2005) [106] (khảo sát ảnh hưởng của liên kết đàn hồi biên tiếp tuyến trên vỏ composite cốt sợi truyền thống): Luận án nâng tầm lý thuyết khi tích hợp vi cơ học nanocomposite FG-CNTRC dị hướng cao, tính phụ thuộc nhiệt độ của polymer matrix, và các mô hình cắt trượt bậc cao HSDT.

Đó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à thách thức các giới hạn của các lý thuyết cơ học kết cấu kinh điển thông qua 4 đóng góp lý thuyết cốt lõi:

  1. Mở rộng Quy tắc hỗn hợp suy rộng (Extended Rule of Mixtures): Tích hợp các tham số hiệu quả CNT ($\eta_1, \eta_2, \eta_3$) được hiệu chuẩn từ mô phỏng Động lực học phân tử (Molecular Dynamics - MD) vào các phương trình cấu tạo của vỏ dị hướng: $$E_{11} = \eta_1 V_{CNT} E_{11}^{CNT} + V_m E^m$$ $$E_{22} = \frac{\eta_2 E_{22}^{CNT} E^m}{V_m E_{22}^{CNT} + \eta_2 V_{CNT} E^m}$$ $$G_{12} = \frac{\eta_3 G_{12}^{CNT} G^m}{V_m G_{12}^{CNT} + \eta_3 V_{CNT} G^m}$$
  2. Mở rộng Lý thuyết Vỏ thoải Donnell phi tuyến: Thiết lập hệ phương trình vi phân đạo hàm riêng phi tuyến ghép cặp hoàn toàn giữa hàm ứng suất Airy ($\Phi$) và hàm độ võng ($w$) kết hợp góc xoay trượt ($\psi_x, \psi_y$), mô tả chính xác tương tác giữa lực màng và mô men uốn trong không gian biến dạng lớn von Kármán.
  3. Mô hình hóa Biên đàn hồi tiếp tuyến tổng quát: Lượng hóa các hệ số độ cứng lò xo tiếp tuyến không thứ nguyên ($\chi_1, \chi_2$) và độ cứng liên kết ($k_1, k_2$), chuyển hóa các điều kiện biên từ hai trạng thái cực đoan (Freely Movable - FM và Fully Immovable - TR) thành một hàm biến thiên liên tục, phản ánh đúng bản chất vật lý của gối tựa công trình.
  4. Chứng minh bước chuyển Paradigm (Paradigm Shift) trong cơ chế mất ổn định: Đưa ra bằng chứng giải tích chứng minh rằng sự tồn tại của điểm rẽ nhánh (bifurcation load) hay điểm cực trị giới hạn (limit point) không phải là thuộc tính cố định của hình học vỏ, mà là hàm phụ thuộc trực tiếp vào tỷ số độ cong/chiều dày ($R/h$) và độ cứng chống dịch chuyển tiếp tuyến biên ($\chi$).
           +-------------------------------------------------------------------+
           |            HỆ PHƯƠNG TRÌNH ĐIỀU KHIỂN CỐT LÕI (FSDT / HSDT)        |
           +-------------------------------------------------------------------+
                                              |
                   +--------------------------+--------------------------+
                   |                                                     |
                   v                                                     v
+------------------------------------+                 +------------------------------------+
|   PHƯƠNG TRÌNH TƯƠNG THÍCH BIẾN DẠNG|                 |     PHƯƠNG TRÌNH CÂN BẰNG PHI TUYẾN |
|   d^4(Phi) = L(w, w) + L(w, w*)    |                 |   D*d^4(w) + L(Phi, w) + q_ext     |
|   + Tác động giãn nở nhiệt Gradient |                 |   - P_elast(K1, K2) = 0            |
+------------------------------------+                 +------------------------------------+
                   |                                                     |
                   +--------------------------+--------------------------+
                                              |
                                              v
           +-------------------------------------------------------------------+
           |     PHÉP CHIẾU GALERKIN & THUẬT TOÁN LẶP PHỤ THUỘC NHIỆT ĐỘ (T-D)  |
           |      w(x, y) = W * sin(alpha*x) * sin(beta*y) -> F(Load, W, T) = 0|
           +-------------------------------------------------------------------+

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

Khung phân tích của luận án kết hợp đồng thời 4 trụ cột lý thuyết: Lý thuyết cấu tạo FG-CNTRC, Động học biến dạng trượt HSDT Reddy, Mô hình tiếp xúc nền Pasternak ($K_1, K_2$), và Cơ học nhiệt đàn hồi phi tuyến. Khung phân tích giải tích được xây dựng trên hệ điều kiện biên tựa bản lề đàn hồi:

  • Tại các cạnh $x = 0, a$: $w = 0$, $M_x = 0$, $N_{xy} = 0$, và $\frac{1}{b}\int_0^b N_x dy + k_1 \frac{1}{b}\int_0^b u|_{x=0,a} dy = 0$.
  • Tại các cạnh $y = 0, b$: $w = 0$, $M_y = 0$, $N_{xy} = 0$, và $\frac{1}{a}\int_0^a N_y dx + k_2 \frac{1}{a}\int_0^a v|_{y=0,b} dx = 0$.

Điều kiện biên này loại bỏ hoàn toàn các sai số làm mềm nhân tạo (artificial softening) thường gặp trong phương pháp số, cho phép kiểm soát chặt chẽ điều kiện biên từ tự do ($\chi = 0$) đến ngàm cứng ($\chi \to \infty$).

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 theo hệ hình thực chứng (Positivism) và chủ nghĩa duy lý diễn dịch trong cơ học môi trường liên tục. Thiết kế nghiên cứu kết hợp phương pháp giải tích hình thức (closed-form analytical formulation) và thuật toán lặp số học (numerical iteration) đối với các đặc trưng phụ thuộc nhiệt độ. Thiết kế đa tầng bao gồm:

  • Cấp độ nano: Trích xuất các tham số đàn hồi của SWCNT (10,10) ($r = 0.678\text{ nm}$, $h = 0.067\text{ nm}$, mô đun đàn hồi $E_{11}^{CNT} = 5.6466\text{ TPa}$, $E_{22}^{CNT} = 7.0800\text{ TPa}$, $G_{12}^{CNT} = 1.9445\text{ TPa}$, hệ số Poisson $\nu_{12}^{CNT} = 0.175$, hệ số giãn nở nhiệt $\alpha_{11}^{CNT} = 3.4584 \times 10^{-6}\text{ K}^{-1}$).
  • Cấp độ vi mô: Tính toán các hệ số hiệu quả $\eta_j$ ($j=1,2,3$) tương ứng với các mức tỷ lệ thể tích $V^*_{CNT}$ thông qua quy tắc hỗn hợp suy rộng.
  • Cấp độ vĩ mô: Thiết lập và giải các phương trình chủ đạo cho tấm và panel vỏ thoải dưới các điều kiện tải trọng phức hợp.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được triển khai theo các bước chuẩn mực nghiêm ngặt:

  1. Thiết lập quan hệ ứng suất - biến dạng: Tích hợp tensor ứng suất phi tuyến von Kármán và biến dạng nhiệt: $${\epsilon} = {\epsilon^0} + z{\kappa} + z^3{\kappa^*} - {\alpha}\Delta T$$
  2. Thiết lập hàm ứng suất Airy và hàm chuyển vị: Biểu diễn trường chuyển vị qua chuỗi hàm lượng giác trực giao Fourier thỏa mãn chính xác các điều kiện biên động học.
  3. Áp dụng tích phân Galerkin: Chiếu các phương trình đạo hàm riêng phi tuyến lên không gian hàm dạng để triệt tiêu số dư, thiết lập phương trình đại số phi tuyến bậc ba liên hệ trực tiếp giữa Tải trọng - Độ võng ($P - W$) hoặc Nhiệt độ - Độ võng ($T - W$).
  4. Tam giác đạc phương pháp (Method Triangulation): So sánh và kiểm chứng chéo kết quả giữa ba cấp lý thuyết CPT, FSDT và HSDT.
  5. Kiểm chứng độ tin cậy ngoại vi (External Validation): Đối chiếu độ chính xác với các kết quả đã công bố của Shen [13], Kiani [101], và Liew et al. [58] trên các trường hợp suy biến chuẩn (độ lệch sai số dưới $0.5%$).
                                QUY TRÌNH NGHIÊN CỨU VÀ TÍNH TOÁN
                                
[Dữ liệu MD cấp Nano] ---> [Quy tắc hỗn hợp suy rộng] ---> [Mô hình hóa FG-CNTRC (UD/FG)]
                                                                     |
                                                                     v
[Hệ phương trình FSDT/HSDT] <--- [Biên đàn hồi (k1, k2)] <--- [Động học von Kármán-Donnell]
          |
          v
[Phương pháp Galerkin] ---> [Hệ phương trình phi tuyến đại số P-W / T-W]
                                         |
               +-------------------------+-------------------------+
               | (Trường hợp T-ID)                                 | (Trường hợp T-D phụ thuộc T)
               v                                                   v
     [Nghiệm giải tích dạng đóng]                         [Thuật toán lặp Picard lặp nhiệt]
               |                                                   |
               +-------------------------+-------------------------+
                                         |
                                         v
                         [ĐƯỜNG CÂN BẰNG VÀ TẢI TỚI HẠN]

Data và phân tích

Dữ liệu tính toán được thiết lập qua hàng trăm kịch bản số trên phần mềm tính toán ký hiệu toán học chuyên dụng (MATLAB/Maple). Các tham số kiểm định bao gồm:

  • Tỷ lệ kích thước cạnh: $a/b \in [0.5, 2.0]$; Tỷ lệ độ mỏng: $a/h \in [10, 100]$.
  • Tỷ số độ cong: $R/a \in [5, 50]$ và $R_x/R_y$.
  • Tham số nền Pasternak: Độ cứng lớp Winkler $K_1 \in [0, 300]$ và độ cứng lớp trượt Pasternak $K_2 \in [0, 50]$.
  • Kiểm tra tính vững (Robustness checks): Đánh giá độ nhạy của kết quả trước các mức độ không hoàn hảo hình dáng ban đầu $W^* \in [0.01, 0.3]$.

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 khoa học đột phá có giá trị thực nghiệm và lý thuyết sâu sắc:

  1. Quy luật ưu việt tuyệt đối của cấu hình FG-X: Trong tất cả các dạng phân bố, kiểu FG-X (CNTs tập trung cực đại tại hai bề mặt ngoài) luôn tạo ra tải trọng tới hạn cơ học và nhiệt độ tới hạn rẽ nhánh cao nhất. Cụ thể, tấm FG-X có tải tới hạn nén cao hơn từ $40%$ đến $65%$ so với phân bố đều (UD) và vượt trội gấp $2.5$ lần so với phân bố FG-O (vốn tập trung CNTs tại mặt trung hòa).
  2. Tác động hai mặt đối nghịch của liên kết đàn hồi biên tiếp tuyến:
    • Đối với tải nén cơ học, việc tăng độ cứng lò xo biên ($k_1, k_2$ hay $\chi_1, \chi_2$) làm tăng đáng kể khả năng chịu tải sau vồng nhờ khai thác hiệu ứng màng kéo căng (tensile membrane action).
    • Đối với tải nhiệt, việc tăng độ cứng liên kết biên lại làm suy giảm nghiêm trọng nhiệt độ tới hạn rẽ nhánh $T_{cr}$, do kết cấu bị cưỡng bức biến dạng nhiệt mạnh hơn, kích hoạt sớm trạng thái mất ổn định.
  3. Sự triệt tiêu điểm rẽ nhánh ở panel vỏ cong chịu nhiệt: Khi panel trụ hoặc panel hai độ cong có độ cong đáng kể ($R/a \le 20$), kết cấu không xuất hiện hiện tượng rẽ nhánh vồng cổ điển mà biến dạng uốn phi tuyến phát triển đơn điệu ngay từ khi bắt đầu gia nhiệt, uốn cong về phía lồi hoặc phía lõm tùy thuộc vào điều kiện ứng suất trước.
  4. Hiệu ứng triệt tiêu độ bền do ứng suất trước cơ - nhiệt kết hợp: Khi tấm và panel chịu áp lực nén trước dọc trục hoặc áp lực mặt ngoài tồn tại trước, khả năng chống chịu nhiệt độ tới hạn bị suy giảm phi tuyến; một lượng tải nén trước bằng $30%$ tải tới hạn có thể làm giảm nhiệt độ tới hạn của kết cấu lên đến $50%$.
  5. Ưu thế áp đảo của tham số nền Pasternak: Lớp cắt trượt Pasternak ($K_2$) thể hiện hiệu quả kiềm chế độ võng sau vồng và nâng cao tải tới hạn cao gấp $3-5$ lần so với lớp lò xo Winkler đơn thuần ($K_1$).
                         ĐÁP ỨNG TẢI TRỌNG - ĐỘ VÕNG SAU VỒNG (POSTBUCKLING)
   Tải P, Nhiệt T
        ^
        |                                     / (FG-X: Khả năng chịu tải cao nhất)
        |                                    /
        |         Điểm rẽ nhánh             /-- (UD: Phân bố đều)
 P_cr --+-----------------x               /
 (T_cr) |                / \             /---- (FG-O: Khả năng chịu tải thấp nhất)
        |  Trạng thái   /   \           /
        |  màng        /     \         /
        |  trước vồng /       \       /   (Đường sau vồng phi tuyến dạng dốc lên)
        |            /         \     /
        +-----------+-----------+---+--------------------------------------------> Độ võng W
       O|           |           |
        |           |           +-- [Hiện tượng Hóp vỏ (Snap-through) ở panel cong]
        |           +-- [Đường rẽ nhánh đối xứng của tấm phẳng hoàn hảo]

Implications đa chiều

  • Về mặt lý thuyết: Hoàn thiện lý thuyết phi tuyến vỏ thoải dị hướng nanocomposite, cung cấp công cụ giải tích chuẩn xác định rõ ranh giới giữa mất ổn định rẽ nhánh và mất ổn định cực trị.
  • Về mặt phương pháp luận: Thiết lập quy trình giải tích - bán giải tích tin cậy cao, đóng vai trò là "chuẩn đối sánh" (benchmark solution) cho việc phát triển và kiểm chứng các thuật toán số, phần tử hữu hạn phi tuyến thương mại (như ANSYS, ABAQUS).
  • Về mặt kỹ thuật thực tiễn: Cung cấp hướng dẫn kỹ thuật cốt lõi cho các kỹ sư kết cấu hàng không vũ trụ: tối ưu hóa vị trí gia cường CNTs theo mẫu FG-X giúp giảm $20-30%$ trọng lượng kết cấu vỏ máy bay siêu thanh mà vẫn bảo toàn độ ổn định nhiệt đàn hồi.
  • Về tiêu chuẩn và chính sách: Đề xuất bổ sung hệ số suy giảm ổn định do liên kết biên đàn hồi vào các quy chuẩn thiết kế kết cấu composite chịu nhiệt độ cao trong ngành công nghiệp quốc phòng và năng lượng.

Limitations và Future Research

Nhằm duy trì tính khách quan và chuẩn mực học thuật, luận án chỉ rõ 4 giới hạn nghiên cứu:

  1. Phạm vi đàn hồi tuyến tính của vật liệu: Nghiên cứu giả định các thành phần vật liệu làm việc hoàn toàn trong miền biến dạng đàn hồi, chưa tính đến hiện tượng dẻo hóa của polymer hay từ biến (creep) khi nhiệt độ tiến sát nhiệt độ chuyển pha thủy tinh ($T_g$).
  2. Giả thiết liên kết pha lý tưởng: Giả định ống nano carbon liên kết hoàn hảo với nền polymer, bỏ qua hiện tượng trượt phân tử tại bề mặt phân chia pha, hiện tượng kết cụm (agglomeration) và độ uốn lượn cục bộ (waviness) của CNTs.
  3. Giới hạn hình học vỏ thoải: Mô hình toán học Donnell bị giới hạn ở các panel cong có độ cong nhỏ ($R \gg a, b, h$), không áp dụng trực tiếp cho các vỏ trụ tròn kín hoàn toàn hoặc vỏ có độ cong sâu.
  4. Trạng thái tải tĩnh: Chưa khảo sát đáp ứng mất ổn định động học (dynamic buckling) dưới tác dụng của sóng xung kích nhiệt hoặc tải nén va đập theo thời gian.

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

  • Mở rộng phân tích ổn định động lực học và tương tác khí - nhiệt - đàn hồi (aerothermoelastic flutter) của panel FG-CNTRC trong dòng khí siêu âm.
  • Phát triển mô hình đa quy mô tích hợp hiện tượng kết cụm CNTs dựa trên lý thuyết Eshelby-Mori-Tanaka kết hợp cơ học phá hủy liên kết phân tử.
  • Nghiên cứu kết cấu tấm vỏ FG-CNTRC thông minh tích hợp các lớp áp điện (piezoelectric) hoặc vật liệu nhớ hình (SMA) nhằm chủ động điều khiển và dập tắt độ võng sau vồng do nhiệt.
  • Triển khai thực nghiệm chế tạo và thử nghiệm nén nhiệt đa trục để kiểm chứng trực tiếp các đường cong sau vồng giải tích.

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

  • Tác động học thuật: Luận án mở ra hướng nghiên cứu mới về ổn định phi tuyến kết cấu nano có liên kết biên thực tế, với tiềm năng trích dẫn cao trong các tạp chí cơ học hàng đầu (như Composite Structures, Thin-Walled Structures, Aerospace Science and Technology).
  • Chuyển đổi công nghiệp: Ứng dụng trực tiếp trong thiết kế vỏ vệ tinh không gian, thân tên lửa đẩy, cánh tà máy bay không người lái (UAV) tốc độ cao, và tấm chắn nhiệt trong lò phản ứng hạt nhân.
  • Lợi ích xã hội và môi trường: Việc tối ưu hóa vật liệu nanocomposite giúp tạo ra các phương tiện vận tải nhẹ hơn, giảm thiểu tiêu hao nhiên liệu hóa thạch và cắt giảm lượng khí thải carbon quy mô toàn cầu.

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

  • Nghiên cứu sinh & Nhà khoa học: Tiếp cận khung giải tích mẫu mực để phát triển các bài toán kết cấu phức tạp hơn (vỏ nón, vỏ cầu, kết cấu xếp lớp chức năng).
  • Giảng viên & Học giả cao cấp: Sử dụng làm tài liệu tham khảo chuyên sâu cho các chương trình đào tạo sau đại học ngành Cơ học vật rắn biến dạng và Kỹ thuật hàng không vũ trụ.
  • Kỹ sư R&D công nghiệp: Ứng dụng công thức giải tích dạng đóng để tính toán nhanh tải trọng phá hủy và biên độ an toàn nhiệt mà không cần chạy các mô hình phần tử hữu hạn đắt đỏ.
  • Cơ quan quản lý & Cơ quan lập quy chuẩn: Khai thác dữ liệu định lượng để xây dựng bộ quy chuẩn thiết kế độ bền kết cấu vật liệu mới tại Việt Nam.

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 hệ phương trình phi tuyến giải tích giải quyết đồng thời ba yếu tố: mô hình vi cơ học FG-CNTRC với các tham số hiệu quả $\eta_j$, điều kiện biên đàn hồi tiếp tuyến nửa cứng ($k_1, k_2$), và trường nhiệt độ không đều dạng hàm trong mặt phẳng. Công trình đã mở rộng trực tiếp Lý thuyết Biến dạng trượt bậc cao HSDT của Reddy và Thuyết vỏ thoải Donnell sang miền bài toán vật liệu nanocomposite chịu trường vật lý phức hợp.

2. Đột phá về phương pháp luận so với các nghiên cứu tiền nhiệm? Trả lời: Khác với cách tiếp cận số thuần túy (IMLS-Ritz của Liew et al., DQM của Ansari et al.) vốn phụ thuộc vào lưới và khó khảo sát tính liên tục của tham số, hoặc cách tiếp cận tiệm cận hai bước của Shen vốn áp dụng cho biên ngàm cứng; luận án sử dụng phương pháp đặt bài toán theo hàm ứng suất Airy kết hợp phép chiếu trực giao Galerkin và thuật toán lặp Picard cho đặc trưng T-D. Phương pháp này cung cấp nghiệm giải tích tường minh, bảo toàn tính liên tục vật lý và đạt độ chính xác tương đương nghiệm giải tích giải tích chính xác.

3. Phát hiện bất ngờ nhất mang tính trực giác ngược (counter-intuitive) là gì? Trả lời: Phát hiện về "nghịch lý biên đàn hồi trong trường nhiệt": Khi tăng độ cứng ràng buộc biên tiếp tuyến ($k_1, k_2$), khả năng chịu tải nén cơ học tăng mạnh, nhưng khả năng chịu nhiệt tới hạn rẽ nhánh lại suy giảm nghiêm trọng. Điều này cảnh báo các nhà thiết kế không được phép tăng cứng biên một cách cơ học khi kết cấu làm việc trong môi trường nhiệt độ cao.

4. Quy trình lặp lại nghiên cứu (Replication Protocol) có được đảm bảo không? Trả lời: Hoàn toàn đảm bảo. Luận án cung cấp đầy đủ hệ phương trình cấu tạo vi mô, toàn bộ các bảng hệ số vật liệu phụ thuộc nhiệt độ của SWCNT (10,10) và polymer (PmPV, PMMA), tường minh các tích phân Galerkin và quy trình thuật toán lặp, cho phép bất kỳ nhà nghiên cứu nào cũng có thể lập trình lại chính xác trên các phần mềm toán học.

5. Tầm nhìn nghiên cứu 10 năm tiếp theo từ bệ phóng của luận án? Trả lời: Phát triển hệ thống phần mềm tính toán chuyên dụng phân tích độ bền - ổn định nhiệt đàn hồi đa trường vật lý (Cơ - Nhiệt - Điện - Từ) cho các kết cấu nanocomposite đa thang đo, hướng tới tự động hóa thiết kế tối ưu hóa topo cho ngành công nghiệp hàng không vũ trụ Việt Nam.

Kết luận

Luận án tiến sĩ của NCS Lê Thị Như Trang đã giải quyết trọn vẹn và xuất sắc các mục tiêu khoa học đề ra, đóng góp những giá trị nền tảng cho cơ học kết cấu nanocomposite:

  1. Xác lập thành công mô hình giải tích toàn diện cho bài toán ổn định nhiệt đàn hồi phi tuyến của tấm, panel trụ và panel hai độ cong FG-CNTRC trên nền đàn hồi Pasternak.
  2. Khám phá và lượng hóa chi tiết ảnh hưởng của liên kết biên đàn hồi tiếp tuyến, giải quyết triệt để sự thiếu nhất quán trong y văn quốc tế về cơ chế mất ổn định nhiệt đàn hồi của panel vỏ cong.
  3. Chứng minh định lượng tính ưu việt vượt trội của cấu hình phân bố ống nano carbon dạng FG-X trong việc tối ưu hóa khả năng kháng vồng cơ - nhiệt.
  4. Xây dựng thuật toán lặp chính xác cao xử lý tính phi tuyến vật liệu phụ thuộc nhiệt độ kết hợp trường nhiệt phân bố không đều.
  5. Cung cấp hệ thống dữ liệu số chuẩn mực, đóng vai trò tham chiếu quốc tế cho các nghiên cứu lý thuyết và ứng dụng thực tiễn trong thiết kế kết cấu hàng không vũ trụ thế hệ mới.