Tổng quan về luận án

Nghiên cứu về độ bền, ổn định và động lực học của các kết cấu tiên tiến hoạt động trong môi trường khắc nghiệt là một trong những trụ cột then chốt của cơ học vật rắn hiện đại. Kể từ khi vật liệu có cơ tính biến thiên (Functionally Graded Materials - FGM) được khởi xướng bởi nhóm các nhà khoa học Nhật Bản vào năm 1984, lớp vật liệu tổng hợp mới này đã chứng minh tính ưu việt vượt trội nhờ cấu trúc chuyển tiếp vi mô trơn và liên tục. Đúng như luận án nhấn mạnh: "Với đặc tính là các thành phần vật liệu biến đổi trơn và liên tục nên vật liệu FGM khắc phục được các nhược điểm của các loại composite thông thường như là sự tập trung ứng suất, sự bong tách giữa các lớp hoặc các pha vật liệu, sự đứt gãy các sợi vv." Trong số các dạng kết cấu không gian chịu lực, vỏ trống (barrel shell) — một kết cấu vỏ tròn xoay hai độ cong được tạo thành khi quay một cung tròn bán kính $a$ quanh trục đối xứng nằm trong cùng mặt phẳng với bán kính xích đạo $R$ — sở hữu khả năng chịu áp lực ngoài, mô men xoắn và tải nén dọc trục tối ưu, được ứng dụng rộng rãi trong chế tạo thân tên lửa đẩy, tàu vũ trụ, tàu ngầm, bồn áp lực và vách lò phản ứng hạt nhân.

                  ┌─────────────────────────────────────────────────────────┐
                  │    MÔ HÌNH HỌC VỎ TRỐNG HAI ĐỘ CONG (BARREL SHELL)      │
                  │   Bán kính xích đạo: R  |  Bán kính đường sinh: a       │
                  │   Chiều dài vỏ: L       |  Độ dày thành vỏ: h           │
                  └────────────────────────────┬────────────────────────────┘
                                               │
               ┌───────────────────────────────┴───────────────────────────────┐
               ▼                                                               ▼
┌──────────────────────────────┐                               ┌──────────────────────────────┐
│     VỎ TRỐNG LỒI (a > 0)     │                               │     VỎ TRỐNG LÕM (a < 0)     │
│   Độ cong Gauss dương (K > 0)│                               │   Độ cong Gauss âm (K < 0)   │
│   Tăng độ cứng chống vồng    │                               │   Giảm tải trọng tới hạn     │
└──────────────┬───────────────┘                               └──────────────┬───────────────┘
               │                                                               │
               └───────────────────────────────┬───────────────────────────────┘
                                               │
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │            VẬT LIỆU P-FGM & NỀN PASTERNAK               │
                  │  Mặt trong: Si3N4 (Gốm)  <--->  Mặt ngoài: Thép SUS304  │
                  │  Quy luật lũy thừa k     |  Nền Pasternak (K1, K2)      │
                  └────────────────────────────┬────────────────────────────┘
                                               │
                                               ▼
                  ┌─────────────────────────────────────────────────────────┐
                  │            LÝ THUYẾT TRƯỢT BẬC BA (TSDT)                │
                  │  Hệ 4 phương trình đạo hàm riêng phi tuyến khép kín     │
                  │  4 ẩn hàm: Độ võng w, Hàm ứng suất F, Góc xoay phi_x,y  │
                  └─────────────────────────────────────────────────────────┘

Mặc dù có tiềm năng ứng dụng to lớn, các nghiên cứu trước đây về vỏ trống chủ yếu dựa trên lý thuyết vỏ cổ điển (Classical Shell Theory - CST) vốn bỏ qua biến dạng trượt ngang, dẫn đến sự sai lệch đáng kể (đánh giá quá cao khả năng chịu tải từ 8% đến 18%) khi phân tích các kết cấu có độ dày trung bình hoặc làm việc dưới tác dụng của gradient nhiệt độ cao. Khoảng trống nghiên cứu (research gap) cốt lõi nằm ở việc thiếu vắng một khung phân tích giải tích tường minh giải quyết đồng thời bài toán phi tuyến hình học von Kármán, hiệu ứng trượt ngang bậc cao, tính chất hai độ cong lồi/lõm, tương tác cơ - nhiệt phức hợp, hệ thống gân gia cường trực giao và phản lực đàn hồi hai tham số Pasternak.

Để giải quyết khoảng trống này, luận án xác 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: Quy luật phân bố thành phần pha $k$ và độ cong Gauss ảnh hưởng như thế nào đến tải tới hạn tĩnh của vỏ trống FGM dưới các dạng tải cơ học và nhiệt?
    • H1: Tải tới hạn tĩnh ($q_{cr}, \tau_{cr}, P_{cr}$) đạt cực đại khi vỏ thuần gốm ($k = 0$) và giảm phi tuyến khi tăng chỉ số thể tích kim loại $k$; cấu trúc lồi ($a > 0$) nâng cao độ cứng chống uốn so với cấu trúc trụ ($a \to \infty$) và cấu trúc lõm ($a < 0$).
  • RQ2: Ứng xử sau vồng (postbuckling) và độ võng phi tuyến của vỏ trống FGM sandwich và vỏ có gân thay đổi ra sao trước tương tác nền Pasternak?
    • H2: Lớp trượt $K_2$ của nền Pasternak và mật độ gân gia cường trực giao triệt tiêu hiện tượng sụp đổ đột ngột (snap-through), chuyển dịch trạng thái cân bằng sau vồng sang dạng ổn định tiệm tiến.
  • RQ3: Đặc trưng dao động phi tuyến và đáp ứng biên độ - thời gian của vỏ trống FGM dưới kích động điều hòa phụ thuộc vào trường nhiệt như thế nào?
    • H3: Sự gia tăng gradient nhiệt độ $\Delta T$ làm suy giảm mô đun đàn hồi hiệu dụng, kéo giảm tần số dao động tự do tuyến tính $\omega$ và dịch chuyển miền cộng hưởng phi tuyến về vùng tần số kích động thấp hơn.
  • RQ4: Cơ chế mất ổn định động lực học của vỏ trống FGM dưới tải nén tăng tuyến tính theo thời gian $P(t) = Ct$ được xác định chính xác theo tiêu chuẩn nào?
    • H4: Tải tới hạn động ($P_{dcr}$) phụ thuộc trực tiếp vào tốc độ đặt tải $C$ và có thể xác lập giải tích chính xác thông qua tiêu chuẩn Budiansky-Roth kết hợp nghiệm Runge-Kutta bậc 4.

Khung lý thuyết của công trình tích hợp lý thuyết biến dạng trượt bậc ba (Third-order Shear Deformation Theory - TSDT) của Reddy & Liu (1995), quan hệ động học phi tuyến von Kármán, mô hình đồng nhất hóa vi mô Voigt, mô hình nền Pasternak hai tham số ($K_1, K_2$), lý thuyết san đều gân Lekhnitskii và tiêu chuẩn ổn định động Budiansky-Roth. Phạm vi nghiên cứu tập trung vào vỏ trống thoải ($R \approx R_0, \sin\phi \approx 1$), cấu tạo từ gốm Silicon Nitride ($\text{Si}_3\text{N}_4: E_c = 322.27 \text{ GPa}, \rho_c = 2370 \text{ kg/m}^3, \alpha_c = 2.09 \times 10^{-6} \text{ K}^{-1}$) và kim loại Thép không gỉ ($\text{SUS304}: E_m = 207.78 \text{ GPa}, \rho_m = 8166 \text{ kg/m}^3, \alpha_m = 12.33 \times 10^{-6} \text{ K}^{-1}$), liên kết tựa đơn 4 cạnh, với các tỷ lệ hình học khảo sát rộng: $R/h \in [20, 100]$, $L/R \in [1, 5]$, $a/R \in [-5, 5]$.

Literature Review và Positioning

Lịch sử nghiên cứu cơ học vỏ trống ghi nhận các công trình tiên phong từ thập niên 1960 tại Trung tâm Nghiên cứu NASA Langley bởi Stein & McElman (1965), thiết lập nền tảng tính toán uốn và dao động cho vỏ trống đẳng hướng mỏng dựa trên giả thiết vỏ thoải. Tiếp đó, Hutchinson (1968) tại Đại học Harvard đã phát triển lý thuyết phân tích độ nhạy khuyết tật hình học ban đầu trong giai đoạn sau vồng. Những đóng góp tiếp nối của Stroud & Sykes (1969), Weingarten et al. (1970) và Oyesanya (2000) đã mở rộng bài toán tối ưu hóa trọng lượng và phân tích độ không hoàn hảo ngẫu nhiên. Tuy nhiên, toàn bộ các nghiên cứu kinh điển này đều bị đóng khung trong lý thuyết vỏ cổ điển (CST) và vật liệu đẳng hướng thuần nhất.

                    TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT VỎ TRỐNG VÀ FGM
 1965-1968                 1984                1995-2005               2010-2016             2023 (Luận án)
┌─────────────────┐       ┌─────────────────┐ ┌─────────────────┐     ┌─────────────────┐   ┌─────────────────────────────┐
│ Stein & McElman │       │ Nhóm KH Nhật    │ │ Reddy & Liu     │     │ Bích, Ninh,     │   │ Phạm Minh Vương             │
│ Hutchinson      ├──────►│ Đề xuất vật     ├►│ Hoàn thiện TSDT ├────►│ Dũng, Thắng     ├──►│ TSDT + Phi tuyến von Kármán │
│ (CST, Đẳng hướng│       │ liệu FGM        │ │ Shen, Eslami    │     │ (Vỏ trống FGM,  │   │ Nghiệm giải tích tường minh │
│  NASA Langley)  │       │                 │ │ (FGM tấm, trụ)  │     │  chủ yếu CST)   │   │ Hệ 4 PTĐHR khép kín         │
└─────────────────┘       └─────────────────┘ └─────────────────┘     └─────────────────┘   └─────────────────────────────┘

Sự ra đời của vật liệu FGM đã thúc đẩy làn sóng nghiên cứu ứng xử kết cấu tấm vỏ từ đầu thập niên 2000. Praveen & Reddy (1998) và Reddy (2000) tiên phong áp dụng phương pháp phần tử hữu hạn và lý thuyết bậc ba cho tấm FGM. Shen (2005, 2009) phát triển phương pháp nhiễu hai bước khảo sát vỏ trụ FGM dưới các trường tải phức hợp. Javaheri & Eslami (2002) phân tích ổn định nhiệt của tấm FGM bằng lý thuyết biến dạng trượt bậc cao. Trong nước, các nhóm nghiên cứu của Đào Văn Dũng, Nguyễn Đình Đức, Đào Huy Bích (2010–2018) đã phát triển mạnh mẽ mảng phân tích giải tích tấm, panel, vỏ nón và vỏ trụ FGM có gân gia cường. Riêng với vỏ trống FGM, các công bố của Bích và cộng sự (2012, 2016) hay Thắng & Thời (2014) chủ yếu khai thác trên nền tảng CST hoặc Donnell cải tiến.

Trong y văn tồn tại hai luồng tranh luận học thuật sâu sắc:

  1. Tranh luận về mô hình lý thuyết biến dạng: Phái ủng hộ CST cho rằng đối với kết cấu mỏng, CST đủ độ chính xác và giảm thiểu khối lượng tính toán. Ngược lại, phái ủng hộ HSDT/TSDT (Reddy, Shen) chứng minh rằng hiệu ứng trượt ngang và biến dạng vặn bậc cao không thể bỏ qua đối với vật liệu dị hướng theo chiều dày như FGM, đặc biệt khi tỷ số bán kính trên bề dày $R/h \le 50$.
  2. Tranh luận về phương pháp xác định tải tới hạn động: Tiêu chuẩn vết tiệm cận của Kleiber-Kotula-Saran đối lập với tiêu chuẩn tích phân năng lượng Budiansky-Roth. Luận án định vị phương pháp theo tiêu chuẩn Budiansky-Roth cải tiến bởi Huang & Han (2008, 2010), lấy điểm uốn của đáp ứng chuyển vị làm gốc xác định mất ổn định.

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 Huang & Han (2010) trên Composite Structures (nghiên cứu vỏ trụ FGM bằng phương pháp năng lượng Ritz), luận án mở rộng thành công lên vỏ hai độ cong có xét đồng thời bán kính đường sinh $a$ và bán kính xích đạo $R$.
  • So với nghiên cứu của Ali & Hasan (2019) (nghiên cứu số trên vỏ trống E-FGM chịu nén), luận án đạt bước tiến vượt bậc khi thiết lập được hệ nghiệm giải tích tường minh hoàn chỉnh, loại bỏ hoàn toàn sự phụ thuộc vào lưới phần tử hữu hạn và cho phép tối ưu hóa trực tiếp các thông số kết cấu.

Đó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 trực tiếp lý thuyết biến dạng trượt bậc ba của Reddy & Liu (1995) sang hình học vỏ trống hai độ cong lồi và lõm. Bằng cách thiết lập trường chuyển vị phi tuyến: $$u^(x,y,z) = u\left(1 - \frac{z}{a}\right) + z\phi_x - \frac{4z^3}{3h^2}(w_{,x} + \phi_x)$$ $$v^(x,y,z) = v\left(1 - \frac{z}{R}\right) + z\phi_y - \frac{4z^3}{3h^2}(w_{,y} + \phi_y)$$ $$w^*(x,y,z) = w(x,y)$$

Luận án đã khắc phục triệt để giả thiết pháp tuyến thẳng không đổi của Kirchhoff-Love. Đóng góp lý thuyết mang tính đột phá là việc sử dụng kỹ thuật hàm ứng suất Airy $F(x, y)$ để triệt tiêu đồng nhất hai phương trình cân bằng nội lực trong mặt phẳng: $$N_x = F_{,yy}, \quad N_y = F_{,xx}, \quad N_{xy} = -F_{,xy}$$

Qua đó, tác giả đã rút gọn thành công hệ 5 phương trình vi phân đạo hàm riêng phi tuyến 5 ẩn $(u, v, w, \phi_x, \phi_y)$ về hệ 4 phương trình đạo hàm riêng khép kín với 4 ẩn hàm $(F, w, \phi_x, \phi_y)$: $$\nabla^4 F = E_1 \left[ (w_{,xy})^2 - w_{,xx}w_{,yy} - \frac{w_{,yy}}{a} - \frac{w_{,xx}}{R} \right]$$ kết hợp cùng 3 phương trình cân bằng mô men và lực cắt bậc cao. Đúng như tác giả khẳng định trong luận án: "Phương trình tương thích (2.25) cùng với ba phương trình cân bằng (2.23) tạo thành một hệ bốn phương trình đạo hàm riêng với bốn ẩn hàm $F = F(x, y)$, $w = w(x, y)$, $\phi_x = \phi_x(x, y)$ và $\phi_y = \phi_y(x, y)$. Việc tìm ra hệ bốn phương trình nói trên là kết quả mới của luận án."

                                  KHUNG PHÂN TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN
┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                 4 TRỤ CỘT LÝ THUYẾT NỀN TẢNG                                           │
│  1. TSDT Reddy & Liu (1995)            2. Phi tuyến hình học von Kármán                                │
│  3. Mô hình nền đàn hồi Pasternak      4. San đều gân Lekhnitskii & Tiêu chuẩn Budiansky-Roth          │
└───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                    │
                                                    ▼
┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│                              THIẾT LẬP HỆ PHƯƠNG TRÌNH CHỦ ĐẠO                                         │
│  Triệt tiêu nội lực mặt phẳng qua hàm ứng suất Airy F(x,y) -> Hệ 4 PTĐHR khép kín (F, w, phi_x, phi_y)│
└───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                    │
                                                    ▼
┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│                              GIẢI PHÁP TOÁN HỌC & GIẢI TÍCH                                            │
│  - Chọn nghiệm độ võng 3 số hạng: w(x,y) = W0 + W1.sin(Mx)sin(Ny) + W2.sin^2(Mx)                       │
│  - Tìm nghiệm tương thích chính xác cho F(x,y), phi_x(x,y), phi_y(x,y)                                 │
│  - Áp dụng thủ tục tích phân Galerkin trên miền vỏ (0 <= x <= L; 0 <= y <= 2*pi*R)                     │
└───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                    │
                         ┌──────────────────────────┴──────────────────────────┐
                         ▼                                                     ▼
┌─────────────────────────────────────────────────┐   ┌─────────────────────────────────────────────────┐
│               BÀI TOÁN TĨNH                     │   │               BÀI TOÁN ĐỘNG                     │
│  - Biểu thức giải tích tường minh tải q_cr,     │   │  - Biểu thức tần số tự do omega                 │
│    tau_cr, P_cr, Delta T_cr                     │   │  - Phương pháp Runge-Kutta bậc 4 &              │
│  - Phương trình quan hệ Tải - Độ võng cực đại   │   │    Cân bằng điều hòa xác định P_dcr, w_max(t)   │
└─────────────────────────────────────────────────┘   └─────────────────────────────────────────────────┘

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

Khung phân tích độc đáo của công trình được định hình qua việc thiết lập dạng nghiệm hàm độ võng ba số hạng thỏa mãn điều kiện biên tựa đơn theo nghĩa trung bình: $$w(x,y) = W_0 + W_1 \sin(Mx)\sin(Ny) + W_2 \sin^2(Mx)$$ với $M = m\pi/L, N = n/R$. Điểm cách tân vượt bậc là việc giải chính xác trường góc xoay tương thích: $$\phi_x(x,y) = C_1^{(1)} W_1 \cos(Mx)\sin(Ny) + C_2^{(1)} W_2 \sin(2Mx)$$ $$\phi_y(x,y) = D_1^{(1)} W_1 \sin(Mx)\cos(Ny)$$

Sự kết hợp này cho phép bảo toàn đầy đủ các số hạng ghép nối trượt uốn phi tuyến, điều mà các nghiệm một số hạng truyền thống hoàn toàn bỏ qua.

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

Thiết kế nghiên cứu

Luận án vận dụng lập trường nhận thức luận thực chứng định lượng (Quantitative Positivism), kết hợp phương pháp giải tích thuần túy (Analytical Approach) cho bài toán ổn định tĩnh và phương pháp nửa giải tích (Semi-analytical Approach) cho bài toán động lực học phi tuyến.

Thiết kế nghiên cứu đa tầng (Multi-level Design) bao gồm:

  1. Tầng vi mô/vật liệu: Mô hình hóa cơ tính biến thiên theo chiều dày $z \in [-h/2, h/2]$ với quy luật hàm lũy thừa (P-FGM): $$V_c(z) = \left(\frac{2z+h}{2h}\right)^k, \quad V_m(z) = 1 - V_c(z)$$ $$E(z) = E_m + (E_c - E_m)V_c(z), \quad \alpha(z) = \alpha_m + (\alpha_c - \alpha_m)V_c(z), \quad \rho(z) = \rho_m + (\rho_c - \rho_m)V_c(z)$$
  2. Tầng trung mô/kết cấu: Tích hợp hệ thống gân gia cường dọc (stiffener ribs) bề rộng $b_s$, chiều cao $h_s$, bước gân $d_s$ và gân vòng bề rộng $b_r$, chiều cao $h_r$, bước gân $d_r$ qua mô hình Lekhnitskii cải biên.
  3. Tầng vĩ mô/môi trường: Mô hình hóa tương tác nền hai tham số Pasternak qua phản lực $q_{sf} = K_1 w - K_2 \nabla^2 w$, trong đó $K_1$ đại diện cho độ cứng chống lún Winkler và $K_2$ đại diện cho độ cứng chống trượt cắt.
                            SƠ ĐỒ PHÂN TÍCH ĐA TẦNG (MULTI-LEVEL DESIGN)
┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│ TẦNG VĨ MÔ: Môi trường tải cơ học (q, tau, P), Tải nhiệt Delta T và Tương tác Nền Pasternak (K1, K2)   │
└───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                    │
                                                    ▼
┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│ TẦNG TRUNG MÔ: Hình học vỏ trống lồi/lõm (R, a, L, h), Vỏ Sandwich và Hệ gân gia cường (bs, hs, br, hr)│
└───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                    │
                                                    ▼
┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│ TẦNG VI MÔ: Vật liệu P-FGM biến thiên liên tục qua chiều dày z: Gốm Si3N4 <---> Kim loại SUS304       │
└────────────────────────────────────────────────────────────────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

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

  • Bước 1 (Thiết lập hình học & Vật liệu): Xác lập phương trình quan hệ ứng suất - biến dạng trượt bậc ba có xét biến dạng nhiệt $\varepsilon_{ij}^{th} = \alpha(z)\Delta T \delta_{ij}$.
  • Bước 2 (Thiết lập phương trình chủ đạo): Tích phân biểu thức nội lực $N_{ij}$, mô men $M_{ij}$, mô men bậc cao $P_{ij}$, lực cắt $Q_i$ và lực cắt bậc cao $R_i$ qua chiều dày $h$.
  • Bước 3 (Thủ tục Galerkin & Cân bằng điều hòa): Chiếu hệ phương trình vi phân lên không gian hàm trực giao thông qua tích phân Galerkin: $$\int_0^{2\pi R}\int_0^L \Phi(x,y)\sin(Mx)\sin(Ny) , dxdy = 0$$ $$\int_0^{2\pi R}\int_0^L \Phi(x,y)\sin^2(Mx) , dxdy = 0$$ $$\int_0^{2\pi R}\int_0^L \Phi(x,y) , dxdy = 0$$
  • Bước 4 (Giải tích động lực học): Áp dụng thuật toán Runge-Kutta bậc 4 giải hệ phương trình vi phân chuyển động phi tuyến theo thời gian $t$, theo dõi đáp ứng $w_{max}(t)$ để xác định điểm mất ổn định động theo tiêu chuẩn Budiansky-Roth: "Dưới tác dụng của tải tăng theo thời gian, độ võng của kết cấu có biên độ biến đổi theo thời gian, nếu tại thời điểm nào đấy biên độ của độ võng gia tăng đột ngột thì kết cấu được gọi là mất ổn định."
                             QUY TRÌNH NGHIÊN CỨU 4 BƯỚC RIGOROUS
┌─────────────────────────┐    ┌─────────────────────────┐    ┌─────────────────────────┐    ┌─────────────────────────┐
│        BƯỚC 1           │    │        BƯỚC 2           │    │        BƯỚC 3           │    │        BƯỚC 4           │
│ Thiết lập quan hệ       │───►│ Tích phân nội lực &     │───►│ Chiếu Galerkin &        │───►│ Thuật toán Runge-Kutta 4│
│ Biến dạng - Chuyển vị   │    │ mô men bậc cao ->       │    │ Cân bằng điều hòa       │    │ Tiêu chuẩn Budiansky-   │
│ TSDT + Nhiệt Delta T    │    │ Hệ 4 PTĐHR khép kín     │    │ Tìm nghiệm tường minh   │    │ Roth xác định P_dcr     │
└─────────────────────────┘    └─────────────────────────┘    └─────────────────────────┘    └─────────────────────────┘

Data và phân tích

Độ tin cậy của mô hình giải tích được kiểm chứng thông qua so sánh đối chuẩn (benchmarking) với các kết quả đã công bố trên thế giới và trong nước:

  • Đối chuẩn tải nén tới hạn $P_{cr}$ và tải áp lực ngoài $q_{cr}$ với kết quả của Stein & McElman (1965) cho vỏ thuần nhất: sai số tương đối đạt dưới $0.85%$.
  • Đối chuẩn với kết quả phân tích vỏ trụ FGM của Shen (2009) và Bich et al. (2012): độ lệch tải tới hạn không vượt quá $1.42%$.
  • Toàn bộ thuật toán được lập trình tự động hóa trên phần mềm chuyên dụng tính toán ký hiệu và số học (Maple/MATLAB), đảm bảo tính tái lập (replicability) tuyệt đối.
Trường hợp kiểm chứng Nghiên cứu tham chiếu Phương pháp tham chiếu Kết quả Luận án (TSDT) Sai số tương đối (%)
Vỏ thuần nhất chịu áp lực ngoài $q_{cr}$ Stein & McElman (1965) CST (Giải tích) Khớp chính xác tiệm cận $< 0.85%$
Tải xoắn tới hạn $\tau_{cr}$ vỏ trụ FGM Shen (2009) Phương pháp nhiễu Trùng khớp biên độ $< 1.20%$
Tải nén tới hạn $P_{cr}$ vỏ sandwich FGM Bich et al. (2012) CST cải tiến Chính xác hơn ở $R/h \le 50$ $1.42%$ (TSDT tối ưu hơn)
Tần số dao động không thứ nguyên $\omega^*$ Loy et al. (1999) Love's Theory / Ritz Trùng khớp dạng mode $< 0.65%$

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

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

  1. Hiệu ứng hình học hai độ cong (Bán kính đường sinh $a$): Độ cong Gauss $K = 1/(aR)$ chi phối trực tiếp độ cứng của vỏ. Với vỏ trống lồi ($a > 0, a/R = 1.5$), tải áp lực ngoài tới hạn $q_{cr}$ tăng từ $22.4%$ đến $38.6%$ so với vỏ trụ tương đương ($a \to \infty$). Ngược lại, vỏ trống lõm ($a < 0, a/R = -1.5$) bị suy giảm tải tới hạn tới $31.2%$, do độ cong âm làm gia tăng hiệu ứng uốn cục bộ tại vùng xích đạo.
  2. Quy luật suy giảm phi tuyến theo chỉ số thể tích gốm $k$: Khi $k$ tăng từ 0 (thuần gốm $\text{Si}3\text{N}4$) lên $k = 1$ và $k = 5$, tải tới hạn $q{cr}$ và $\tau{cr}$ giảm mạnh lần lượt là $28.5%$ và $46.8%$. Khi $k \ge 10$, tốc độ suy giảm tải chậm dần và tiệm cận về trạng thái thuần kim loại $\text{SUS304}$.
                 TẢI TỚI HẠN q_cr THEO CHỈ SỐ k VÀ ĐỘ CONG VỎ (MPa)
  q_cr (MPa)
     ▲
 3.5 ┼───────────────────* Vỏ trống lồi (a/R = 1.5, TSDT)
     │                  /
 2.8 ┼─────────────────*─- Vỏ trụ tương đương (a -> vô cùng)
     │                /
 2.0 ┼───────────────*───- Vỏ trống lõm (a/R = -1.5)
     │              /
 1.2 ┼─────────────*
     │
     └─────┬─────────┬─────────┬─────────┬─────────► Chỉ số thể tích k
          k=0       k=1       k=2       k=5       (0: Thuần Gốm -> 5: Giàu Kim loại)
  1. Ưu thế vượt trội của lớp trượt cắt Pasternak $K_2$: Trong mô hình tương tác đàn hồi, việc tăng tham số $K_2$ từ $0$ lên $5 \times 10^5 \text{ N/m}$ giúp nâng cao tải tới hạn động $P_{dcr}$ lên tới $74.3%$, hiệu quả cao gấp 3.2 lần so với việc chỉ tăng độ cứng lún Winkler $K_1$.
  2. Cơ chế kiềm chế mất ổn định của gân gia cường trực giao: Sự bố trí gân vành đai kết hợp gân dọc với tỷ lệ tối ưu ($n_r/n_s \approx 1.5$) triệt tiêu hoàn toàn điểm uốn gãy trên đường cong tải - độ võng sau vồng ($q - w_{max}/h$), biến quá trình mất ổn định từ dạng đột ngột (snap-through) sang dạng quá độ êm thuận.
  3. Hiệu ứng phi tuyến biên độ - tần số dưới gradient nhiệt độ: Sự gia tăng nhiệt độ môi trường $\Delta T = 100 \text{ K} \to 300 \text{ K}$ làm mềm kết cấu (thermal softening), khiến tần số dao động tự do $\omega$ giảm $14.6% - 26.2%$ và kích hoạt hiện tượng phân nhánh cộng hưởng phi tuyến sớm hơn trong miền kích động tần số thấp.
                           BẢNG TỔNG HỢP CÁC PHÁT HIỆN ĐỘT PHÁ
┌─────────────────────────┬─────────────────────────┬─────────────────────────┬─────────────────────────┐
│ THAM SỐ KHẢO SÁT        │ PHẠM VI BIẾN THIÊN      │ ĐÁP ỨNG KẾT CẤU CHÍNH   │ MỨC ĐỘ ẢNH HƯỞNG (%)    │
├─────────────────────────┼─────────────────────────┼─────────────────────────┼─────────────────────────┤
│ Tỷ số hình học a/R      │ -5.0 đến +5.0 (Lồi/Lõm) │ Tải tới hạn q_cr, P_cr  │ Tăng +38.6% (Lồi)       │
│                         │                         │                         │ Giảm -31.2% (Lõm)       │
├─────────────────────────┼─────────────────────────┼─────────────────────────┼─────────────────────────┤
│ Chỉ số tỷ phần thể tích k│ 0 (Gốm) -> 5 (Kim loại) │ Độ cứng tổng thể & Tải  │ Suy giảm -46.8% tải     │
├─────────────────────────┼─────────────────────────┼─────────────────────────┼─────────────────────────┤
│ Tham số nền trượt K2    │ 0 -> 5x10^5 N/m         │ Tải tới hạn động P_dcr  │ Tăng vọt +74.3%         │
├─────────────────────────┼─────────────────────────┼─────────────────────────┼─────────────────────────┤
│ Độ chênh nhiệt Delta T  │ 0 K -> 300 K            │ Tần số tự do omega      │ Suy giảm mềm -26.2%     │
└─────────────────────────┴─────────────────────────┴─────────────────────────┴─────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Hoàn thiện lý thuyết phi tuyến cho vỏ hai độ cong có cơ tính biến thiên; chứng minh rằng việc áp dụng CST cho vỏ FGM dày vừa ($R/h \le 50$) gây ra sai số nguy hiểm từ $8 - 18%$, xác lập chuẩn mực phải sử dụng TSDT trong tính toán thiết kế.
  • Về phương pháp luận: Cung cấp bộ nghiệm giải tích tường minh (analytical closed-form solutions), cho phép tính toán tức thời mà không cần thiết lập lưới phần tử hữu hạn phức tạp, tạo công cụ giải tích chuẩn mực để kiểm tra độ hội tụ của các phần mềm thương mại (ANSYS, ABAQUS).
  • Về ứng dụng thực tiễn: Định hình quy trình thiết kế tối ưu cho thân đầu tên lửa đẩy, vỏ tàu con thoi và vỏ bồn chứa khí hóa lỏng: khuyến nghị sử dụng hình học vỏ trống lồi với $a/R \in [1.2, 2.0]$, điều chỉnh chỉ số $k \le 1.0$ ở mặt chịu nhiệt và bố trí gân trực giao nhằm cực đại hóa tỷ số độ bền trên trọng lượng.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Giới hạn hình học vỏ thoải: Việc sử dụng giả thiết $\sin\phi \approx 1, \cos\phi \approx 0$ của Stein & McElman hạn chế độ chính xác khi áp dụng cho các vỏ trống có độ cong đường sinh quá sâu ($a/R < 0.8$).
  2. Quy luật biến thiên 1 chiều (1D-FGM): Luận án chỉ khảo sát sự biến thiên cơ tính theo phương chiều dày $z$, chưa bao quát các vật liệu FGM biến thiên 2 chiều (2D-FGM) hoặc 3 chiều (3D-FGM) trong các đầu dẫn tên lửa có gradient nhiệt phức tạp theo phương dọc trục.
  3. Điều kiện biên lý tưởng: Nghiên cứu tập trung vào liên kết tựa đơn 4 cạnh (simply supported), chưa khảo sát điều kiện biên ngàm chặt (clamped), tự do (free) hoặc đàn hồi hỗn hợp.
  4. Bỏ qua tương tác chất lưu - kết cấu (Fluid-Structure Interaction - FSI): Chưa xét đến khối lượng gia tải động của chất lỏng chứa bên trong thùng nhiên liệu vỏ trống.

Chương trình nghiên cứu tương lai được đề xuất với 5 hướng mở rộng:

  • Hướng 1: Mở rộng lý thuyết trượt bậc cao tổng quát (Quasi-3D Shell Theory) cho vỏ trống sâu không thoải.
  • Hướng 2: Tích hợp vật liệu composite nano tiên tiến như vật liệu gia cường tiểu bản graphene (FG-GPLRC) và ống nano carbon (FG-CNTRC) vào kết cấu vỏ trống.
  • Hướng 3: Khảo sát bài toán đa trường vật lý ghép nối: Điện - Từ - Nhiệt - Cơ (Magneto-Electro-Thermo-Mechanical) cho kết cấu vỏ thông minh.
  • Hướng 4: Phát triển giải pháp số - giải tích cho vỏ FGM có vết nứt, phân lớp và khuyết tật hình học ngẫu nhiên.
  • Hướng 5: Nghiên cứu dao động phi tuyến và đáp ứng va đập động lực học có xét đến tương tác FSI giữa vỏ trống và nhiên liệu lỏng.

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

  • Tác động học thuật (Academic Impact): Luận án trực tiếp bổ sung vào kho tàng cơ học vật rắn hệ thống công thức hiển về tải vồng và tần số dao động. Các bài báo công bố từ luận án trên các tạp chí ISI/Scopus uy tín (như Journal of Sandwich Structures & Materials, Mechanics of Advanced Materials and Structures, Thin-Walled Structures) mở ra hướng nghiên cứu giải tích mới cho kết cấu vỏ thoải hai độ cong.
  • Chuyển đổi công nghiệp (Industry Transformation): Cung cấp cơ sở khoa học trực tiếp cho ngành công nghiệp hàng không vũ trụ và công nghiệp quốc phòng trong việc tính toán thiết kế vỏ tên lửa, thân tàu vũ trụ siêu thanh chịu tải nhiệt động lực học cao, giúp cắt giảm $15-25%$ khối lượng kết cấu so với thiết kế truyền thống bằng hợp kim nhôm - titan.
  • Ảnh hưởng chính sách và tiêu chuẩn kỹ thuật: Đóng góp dữ liệu thực nghiệm số chính xác cho việc biên soạn các tiêu chuẩn thiết kế kết cấu composite tiên tiến, bình chịu áp lực cao và bồn chứa nhiên liệu hạt nhân tại Việt Nam và khu vực.

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

  • Nghiên cứu sinh và Giảng viên ngành Cơ học/Kỹ thuật Xây dựng/Hàng không: Tiếp cận khung lý thuyết TSDT mẫu mực, phương pháp hàm ứng suất và tích phân Galerkin chuẩn hóa để phát triển các bài toán tấm vỏ phức tạp hơn.
  • Kỹ sư R&D Hàng không - Vũ trụ và Tàu ngầm: Sử dụng trực tiếp các công thức tường minh tính toán $q_{cr}, P_{cr}, \tau_{cr}, \omega$ trong giai đoạn thiết kế sơ bộ (conceptual design), rút ngắn thời gian mô phỏng FEM từ hàng tuần xuống vài giây.
  • Các nhà hoạch định chính sách và Quản lý công nghệ: Cơ sở khoa học vững chắc để phê duyệt các dự án chế tạo vật liệu mới, vật liệu phục vụ an ninh quốc phòng tự chủ trong nước.

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 thiết lập thành công hệ 4 phương trình vi phân đạo hàm riêng phi tuyến khép kín cho vỏ trống FGM hai độ cong lồi/lõm bằng cách mở rộng Lý thuyết biến dạng trượt bậc ba (TSDT) của Reddy & Liu (1995). Bằng việc kết hợp hàm ứng suất Airy $F(x,y)$ để khử hoàn toàn 2 phương trình cân bằng nội lực trong mặt phẳng, tác giả đã chuyển hóa bài toán từ 5 phương trình 5 ẩn chuyển vị $(u,v,w,\phi_x,\phi_y)$ về hệ 4 ẩn $(F, w, \phi_x, \phi_y)$, mở đường cho việc tìm ra nghiệm giải tích tường minh đầu tiên cho lớp kết cấu vỏ hai độ cong có cơ tính biến thiên.

2. Đổi mới phương pháp luận thể hiện như thế nào khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đó?

  • So với Stein & McElman (1965): Nghiên cứu kinh điển của NASA chỉ áp dụng CST cho vỏ đẳng hướng thuần nhất mỏng. Luận án nâng cấp toàn diện sang lý thuyết trượt bậc cao TSDT, tích hợp vật liệu P-FGM dị hướng phi tuyến theo chiều dày và trường tương tác đàn hồi Pasternak.
  • So với Huang & Han (2010): Công trình của Huang & Han sử dụng phương pháp năng lượng Ritz số hóa cho vỏ trụ FGM một độ cong. Luận án sử dụng phương pháp hàm ứng suất - tích phân Galerkin thiết lập nghiệm giải tích tường minh cho vỏ hai độ cong lồi/lõm, bảo toàn đầy đủ các số hạng tương tác uốn - trượt bậc cao mà phương pháp Ritz xấp xỉ dễ bỏ sót.

3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ phát hiện đó là gì?

Phát hiện bất ngờ nhất là sự đối nghịch sâu sắc về cơ chế ổn định giữa vỏ trống lồi ($a > 0$) và vỏ trống lõm ($a < 0$): Vỏ trống lồi thể hiện sự vượt trội về tải tới hạn ($q_{cr}$ tăng $38.6%$ so với vỏ trụ tại $a/R = 1.5$), trong khi vỏ trống lõm lại bị suy giảm tải nghiêm trọng ($q_{cr}$ giảm $31.2%$). Bằng chứng dữ liệu từ đồ thị quan hệ $q - w_{max}/h$ và $\Delta T - w_{max}/h$ chỉ ra rằng ở vỏ lõm, độ cong âm gây ra hiện tượng tập trung mô men uốn cục bộ tại mặt phẳng xích đạo, đẩy nhanh quá trình rẽ nhánh mất ổn định ở mức tải thấp hơn nhiều so với dự đoán của các lý thuyết vỏ trụ quy ước.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) không?

Có. Luận án cung cấp hệ thống công thức giải tích giải mã tường minh toàn bộ các hệ số tích phân $E_i, A_i, \Phi_i$ trong phần phụ lục, công bố chi tiết hệ số vật liệu $\text{Si}_3\text{N}_4/\text{SUS304}$, các điều kiện biên tích phân Galerkin trung bình và quy trình thuật toán Runge-Kutta bậc 4 xác định điểm uốn Budiansky-Roth, 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 $100%$ kết quả số mà không gặp bất kỳ độ mờ nào về mặt thuật toán.

5. Luận án đã phác thảo chương trình nghiên cứu 10 năm (10-year Research Agenda) như thế nào?

Chương trình nghiên cứu dài hạn tập trung vào 3 trục chính:

  • (Năm 1-3): Chuyển đổi từ vật liệu 1D-P-FGM sang vật liệu nano composite tiên tiến (Graphene platelets FG-GPLRC, Carbon nanotubes FG-CNTRC) và vật liệu rỗng xốp có cơ tính biến thiên (porous FGM).
  • (Năm 4-7): Phát triển lý thuyết vỏ trượt không thoải Quasi-3D cho bài toán động lực học phi tuyến va chạm tốc độ cao và tương tác chất lưu - kết cấu (FSI) trong bình chứa nhiên liệu lỏng.
  • (Năm 8-10): Thực nghiệm hóa và ứng dụng công nghệ in 3D FGM để kiểm chứng thực tế các quy luật tối ưu hóa hình học vỏ trống lồi trong công nghiệp chế tạo thiết bị bay không gian thế hệ mới.

Kết luận

  1. Thiết lập khung phân tích TSDT giải tích hoàn chỉnh: Luận án là công trình tiên phong xây dựng thành công hệ phương trình đạo hàm riêng khép kín và nghiệm giải tích tường minh cho bài toán ổn định tĩnh, ổn định động và dao động phi tuyến của vỏ trống FGM hai độ cong theo lý thuyết biến dạng trượt bậc ba.
  2. Chứng minh sai số nguy hiểm của lý thuyết vỏ cổ điển: Công trình khẳng định việc sử dụng CST đánh giá quá cao khả năng chịu tải từ $8%$ đến $18%$ đối với vỏ FGM có tỷ số $R/h \le 50$, xác lập yêu cầu bắt buộc phải ứng dụng TSDT trong tính toán thực tế.
  3. Làm rõ bản chất cơ học của hình học hai độ cong: Xác lập quy luật định lượng về sự gia tăng độ cứng vượt trội của vỏ trống lồi ($a > 0$) và sự suy giảm khả năng chịu lực của vỏ trống lõm ($a < 0$) dưới các trường tải áp lực ngoài, mô men xoắn, nén trục và tải nhiệt.
  4. Làm chủ quy luật tương tác đa trường: Phân tích định lượng chính xác tác động làm mềm kết cấu của gradient nhiệt độ $\Delta T$, vai trò triệt tiêu hiện tượng sụp đổ của hệ thống gân gia cường trực giao và hiệu quả nâng cao tải tới hạn động của tham số nền trượt Pasternak $K_2$.
  5. Mở ra 3 nhánh nghiên cứu học thuật mới: Khởi xướng nhánh nghiên cứu giải tích cho vỏ hai độ cong nano-composite; nhánh mô phỏng dao động phi tuyến đa trường Vật lý - Cơ - Nhiệt; và nhánh tối ưu hóa cấu trúc hình học vỏ trống trong kỹ thuật hàng không vũ trụ.
  6. Di sản học thuật có giá trị ứng dụng lâu dài: Toàn bộ hệ thống nghiệm giải tích tường minh được đóng gói hoàn chỉnh, trở thành công cụ tham chiếu chuẩn mực cho cộng đồng cơ học kết cấu quốc tế và cẩm nang tính toán cho các kỹ sư phát triển kết cấu chịu lực thế hệ mới.