Luận án tiến sĩ: Phân tích ổn định tĩnh tấm vỏ FGM có gân gia cường chịu tải cơ nhiệt - Nguyễn Thị Nga (ĐHQG HN)

Tổng quan về luận án

Nghiên cứu "Phân tích ổn định tĩnh của tấm và vỏ cơ tính biến thiên có gân gia cường chịu tải cơ và nhiệt" do NCS. Nguyễn Thị Nga thực hiện (Chuyên ngành Cơ học vật rắn, Mã số: 62440107, Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội, 2018) dưới sự hướng dẫn khoa học của PGS. Vũ Đỗ Long và cố GS. Đào Văn Dũng, cùng sự định hướng khoa học từ GS. Đào Huy Bích, đại diện cho một bước tiến tiên phong trong lĩnh vực cơ học kết cấu composite tiên tiến. Trong kỹ thuật hiện đại như công nghiệp hàng không vũ trụ, lò phản ứng hạt nhân, tàu ngầm và bồn chứa nhiên liệu áp lực cao, các kết cấu tấm và vỏ phải hoạt động trong điều kiện môi trường khắc nghiệt chịu đồng thời gradient nhiệt độ cực lớn và tải trọng cơ học phức tạp.

Về bối cảnh khoa học, vật liệu cơ tính biến thiên (Functionally Graded Material - FGM) được định vị là: "Vật liệu cơ tính biến thiên có tên tiếng anh là Functionally Graded Material, thường được gọi tắt là FGM, là một loại composite thế hệ mới được nghiên cứu và phát triển lần đầu tiên bởi một nhóm các nhà khoa học ở viện Sendai của Nhật Bản vào năm 1984". FGM kết hợp ưu việt giữa pha gốm chịu nhiệt, chống mài mòn cao ($E_c = 380 \times 10^9\text{ Pa}, \alpha_c = 7.4 \times 10^{-6}\text{ K}^{-1}$) và pha kim loại dẻo dai, chịu lực tốt ($E_m = 70 \times 10^9\text{ Pa}, \alpha_m = 23 \times 10^{-6}\text{ K}^{-1}$), loại bỏ hoàn toàn hiện tượng phân tách lớp (delamination) vốn là nhược điểm chí mạng của composite nhiều lớp truyền thống.

Khoảng trống nghiên cứu (Research Gap) cốt lõi được xác định rõ ràng: "Đã có nhiều công trình nghiên cứu kết cấu FGM khi sử dụng lý thuyết tấm và vỏ cổ điển. Các kết quả thu được chỉ phù hợp với những kết cấu thành mỏng, còn đối với kết cấu thành dày hơn thì cần phải sử dụng lý thuyết biến dạng trượt bậc nhất hoặc bậc ba. Đây vẫn là vấn đề mở, nhất là sử dụng lý thuyết biến dạng trượt bậc ba cho kết cấu FGM và kết cấu sandwich FGM có gân gia cường cũng là FGM". Hầu hết các nghiên cứu trước đây (như Shen, Eslami, Najafizadeh) chỉ tập trung vào kết cấu FGM không gân hoặc FGM có gân thuần nhất đẳng hướng, bỏ qua sự biến thiên cơ tính liên tục trong bản thân hệ gân gia cường và bỏ qua biến dạng trượt ngang bậc cao.

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu và giả thuyết khoa học chính:

1. RQ1 & H1: Ứng xử mất ổn định tĩnh phi tuyến của tấm chữ nhật FGM có gân FGM lệch tâm (ES-FGM) trên nền đàn hồi hai tham số chịu nén và nhiệt tuân theo quy luật nào khi giải tích bằng lý thuyết biến dạng trượt bậc nhất (FSDT) và bậc ba (TSDT)? Giả thuyết: Gân FGM triệt tiêu sự gián đoạn ứng suất tại mặt tiếp xúc và tăng đáng kể tải tới hạn $F_{cr}$.
2. RQ2 & H2: Mô hình vỏ trụ tròn sandwich FGM có gân FGM trên nền Pasternak phản ứng ra sao dưới tác dụng của tải nén dọc trục và trường nhiệt phụ thuộc nhiệt độ (T-D)? Giả thuyết: TSDT phản ánh chính xác sự suy giảm độ cứng tiếp tuyến và hiện tượng nhảy bậc tải sau tới hạn hơn lý thuyết cổ điển (CST).
3. RQ3 & H3: Phương pháp giải tích nào cho phép xác định tường minh tải nén tới hạn $P_{cr}$ và áp lực ngoài tới hạn $q_{cr}$ của vỏ nón cụt sandwich ES-FGM? Giả thuyết: Thiết lập hệ phương trình vi phân biến dạng và chuyển vị trong hệ tọa độ cong kết hợp phép chiếu Galerkin cho nghiệm giải tích dạng đóng.
4. RQ4 & H4: Mức độ ảnh hưởng định lượng của các tham số hình học ($a/b, R/h, L/R$), tham số vật liệu ($k, k_2, k_3$), độ không hoàn hảo ($\xi \in [0, 0.1]$) và hệ số nền đàn hồi ($K_1, K_2$) tới độ tăng nhiệt tới hạn $\Delta T_{cr}$ và đường cong tải - độ võng sau tới hạn ($W/h$).

Khung lý thuyết của luận án tích hợp liên hoàn giữa Lý thuyết biến dạng trượt bậc nhất (FSDT), Lý thuyết biến dạng trượt bậc ba (TSDT của J.N. Reddy), Lý thuyết phi tuyến hình học Von Kármán, Kỹ thuật san đều tác dụng gân cải tiến của Lekhnitskii, Lý thuyết tương tác kết cấu - nền Pasternak ($q_f = K_1 w - K_2 \nabla^2 w$) và Tiêu chuẩn mất ổn định rẽ nhánh Euler-Poincaré. Đóng góp đột phá mang tính định lượng bao gồm việc thiết lập thành công hệ nghiệm giải tích hiển (closed-form solutions), cho phép tính toán tức thời tải trọng tới hạn và dự báo chính xác ứng xử sau tới hạn mà không phụ thuộc vào sai số chia lưới hay hội tụ của phương pháp số (FEM). Phạm vi nghiên cứu bao quát các dạng hình học cơ bản (tấm phẳng, vỏ trụ tròn sandwich, vỏ nón cụt sandwich) với vật liệu FGM biến thiên theo quy luật lũy thừa (P-FGM) và Sigmoid (S-FGM), khảo sát cả tính chất độc lập nhiệt độ (T-ID) và phụ thuộc nhiệt độ (T-D).

Literature Review và Positioning

Nghiên cứu về ứng xử cơ học của kết cấu FGM đã trải qua hơn ba thập kỷ phát triển với các trường phái lý thuyết rõ rệt. Nhóm nghiên cứu của Eslami và cộng sự (60, 61, 62, 89) là những người tiên phong khảo sát ổn định tuyến tính của tấm FGM chịu tải cơ và nhiệt dựa trên lý thuyết tấm cổ điển (CPT), sau đó mở rộng sang tấm dày bằng lý thuyết biến dạng trượt bậc cao. Shariat và Eslami (89) đã công bố các nghiệm giải tích chuẩn cho tấm FGM không gân chịu nén. Đồng thời, Shen và các cộng sự (97, 98, 100, 101, 108, 137, 138) đã phát triển sâu phương pháp khai triển tham số bé kết hợp kỹ thuật nhiễu bậc cao (perturbation technique) để vẽ đường cong tải - độ võng sau tới hạn của vỏ trụ và tấm FGM trong môi trường nhiệt, có xét đến tính chất vật liệu phụ thuộc nhiệt độ (T-D).

Tuy nhiên, trong y văn quốc tế tồn tại hai cuộc tranh luận lý thuyết và phương pháp luận lớn:

• Tranh luận 1 (Lý thuyết biến dạng trượt vs. Lý thuyết cổ điển): CPT và CST bỏ qua hoàn toàn biến dạng trượt ngang ($\gamma_{xz} = \gamma_{yz} = 0$), dẫn đến việc đánh giá quá cao (overestimation) độ cứng uốn và tải trọng tới hạn của kết cấu thành dày trung bình ($h/a > 0.05$ hoặc $h/R > 0.02$). Các công trình của Zenkour (139, 140) và Zenkour - Sobhy (143) chỉ ra rằng lý thuyết biến dạng trượt hàm sin hoặc TSDT là bắt buộc đối với kết cấu tấm vỏ sandwich nhiều lớp nhằm thỏa mãn điều kiện ứng suất cắt ngang bằng không tại hai mặt đáy tự do mà không cần dùng hệ số điều chỉnh trượt $k_s$.
• Tranh luận 2 (Mô hình hóa gân gia cường - Gân thuần nhất vs. Gân FGM): Đào Huy Bích và cộng sự (14) đề xuất mô hình gân thuần nhất (gân hoàn toàn bằng kim loại hoặc hoàn toàn bằng gốm tùy thuộc vào bề mặt gia cố). Ngược lại, Najafizadeh và cộng sự (78) đề xuất mô hình gân FGM có thành phần biến đổi liên tục theo chiều dày. Việc áp dụng gân thuần nhất gây ra sự nhảy bậc đột ngột về mô đun đàn hồi tại mặt tiếp xúc giữa tấm và gân ($z = -h/2$), tạo nên tập trung ứng suất cục bộ.

Vị trí của luận án trong bản đồ học thuật (Academic Positioning) được xác lập tại điểm giao thoa phức tạp nhất: Kết hợp đồng thời hình học kết cấu sandwich (với lõi FGM hoặc lớp phủ FGM), hệ gân gia cường lệch tâm cũng làm bằng vật liệu FGM biến thiên liên tục ($E_{sx}(z), E_{sy}(z)$), đặt trên nền đàn hồi Pasternak ($K_1, K_2$), chịu trường tải cơ nhiệt kết hợp, và được giải quyết hoàn toàn bằng cách tiếp cận giải tích giải mã qua hàm ứng suất Airy và phương pháp Galerkin.

So sánh với các nghiên cứu quốc tế điển hình:

• So với Brush & Almroth (18) và Seide (91) về vỏ nón cụt đẳng hướng chịu áp lực ngoài và nén dọc trục, nghiệm giải tích của luận án khi suy biến về trường hợp kết cấu thuần nhất ($k=0$) cho kết quả trùng khít hoàn toàn, đồng thời mở rộng biên độ phân tích cho kết cấu sandwich FGM có cấu trúc phức tạp.
• So với Sofiyev (111, 114, 117, 120) chuyên sử dụng lý thuyết Donnell cải tiến và phương pháp Galerkin cho vỏ nón FGM không gân, luận án đã tích hợp thành công hệ gân trực giao FGM theo kỹ thuật san đều Lekhnitskii (68), giải quyết triệt để bài toán phi tuyến hình học Von Kármán trong môi trường nhiệt độ cao.

Đó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ú hệ thống lý thuyết cơ học vật rắn biến dạng qua các đóng góp nền tảng:

• Mở rộng lý thuyết biến dạng trượt bậc nhất (FSDT) và bậc ba (TSDT) của Reddy: Tích hợp trường dịch chuyển phi tuyến và ten-xơ biến dạng Von Kármán cho kết cấu tấm và vỏ có cấu trúc ba pha phức tạp (tấm nền FGM - gân dọc FGM - gân ngang FGM). Hệ thức liên hệ ứng suất - biến dạng Hooke dạng mở rộng được xây dựng đồng nhất: $$\sigma_{x}^{st} = E_{sx}(z) \varepsilon_x - E_{sx}(z) \alpha_{sx}(z) \Delta T, \quad \sigma_{y}^{st} = E_{sy}(z) \varepsilon_y - E_{sy}(z) \alpha_{sy}(z) \Delta T$$
• Thiết lập mô hình tương thích cơ tính tại mặt tiếp xúc: Đảm bảo tính liên tục tuyệt đối của mô đun đàn hồi $E(z)$ và hệ số dãn nở nhiệt $\alpha(z)$ tại vị trí tiếp giáp $z = -h/2$. Khi $z = -h/2$, $E_p(-h/2) = E_m$ và $E_{sx}(-h/2) = E_m$, loại bỏ hoàn toàn sự gián đoạn cơ học vốn có trong các mô hình gân thuần nhất truyền thống.
• Tiêu chuẩn ổn định rẽ nhánh tĩnh Euler-Poincaré: Được cụ thể hóa cho bài toán phi tuyến đa tham số. Luận án khẳng định định nghĩa cốt lõi: "Trạng thái cân bằng ban đầu của vật thể đàn hồi được gọi là ổn định nếu dưới tác dụng của lực đã cho và điều kiện biên đã biết không tồn tại trạng thái cân bằng lân cận nào khác với trạng thái cân bằng ban đầu. Còn nếu có ít nhất một trạng thái cân bằng lân cận khác với trạng thái cân bằng ban đầu thì trạng thái cân bằng ban đầu là không ổn định".

                      ┌──────────────────────────────────────────────┐
                      │    Vật liệu FGM & Mô hình Nhiệt - Đàn hồi    │
                      │  (P-FGM, S-FGM, Thuộc tính T-D: Pr(T), T-ID)  │
                      └──────────────────────┬───────────────────────┘
                                             │
                                             ▼
                      ┌──────────────────────────────────────────────┐
                      │  Động học biến dạng trượt bậc cao FSDT/TSDT  │
                      │   + Phi tuyến hình học Von Kármán (w, w*)    │
                      └──────────────────────┬───────────────────────┘
                                             │
                                             ▼
                      ┌──────────────────────────────────────────────┐
                      │  Kỹ thuật san đều gân FGM lệch tâm (ES-FGM)  │
                      │  theo Lekhnitskii + Nền Pasternak (K1, K2)   │
                      └──────────────────────┬───────────────────────┘
                                             │
                                             ▼
                      ┌──────────────────────────────────────────────┐
                      │  Hệ 4 phương trình đạo hàm riêng phi tuyến   │
                      │        (Hàm ứng suất Airy f & w, φx, φy)     │
                      └──────────────────────┬───────────────────────┘
                                             │
                                             ▼
                      ┌──────────────────────────────────────────────┐
                      │     Phương pháp chiếu giải tích Galerkin     │
                      │    => Tải tới hạn Fcr, Pcr, qcr & ΔTcr       │
                      │    => Đường cong tải - độ võng sau tới hạn   │
                      └──────────────────────────────────────────────┘

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

Khung phân tích của luận án kết hợp chặt chẽ ba trụ cột lý thuyết: (1) Cơ học vật liệu composite dị tính liên tục, (2) Lý thuyết đàn hồi phi tuyến hình học kết cấu thành mỏng/dày trung bình, và (3) Phương pháp hàm ứng suất vi phân kết hợp kỹ thuật biến phân Galerkin.

Điểm đột phá nằm ở việc chuyển đổi hệ 5 phương trình cân bằng phi tuyến phức tạp theo các thành phần nội lực ($N_i, M_i, P_i, Q_i, R_i$) về hệ 4 phương trình đạo hàm riêng chứa 4 hàm chưa biết ($w, \phi_x, \phi_y, f$) nhờ đưa vào hàm ứng suất Airy $f(x,y)$ thỏa mãn: $$N_x = f_{,yy}, \quad N_y = f_{,xx}, \quad N_{xy} = -f_{,xy}$$ Phương trình tương thích biến dạng tổng quát cho tấm có xét độ không hoàn hảo ban đầu $w^(x,y)$ được thiết lập chính xác: $$A_{11}^ f_{,xxxx} + (A_{66}^* + 2A_{12}^) f_{,xxyy} + A_{22}^ f_{,yyyy} - B_{21}^* \phi_{x,xxx} - (B_{11}^* - B_{66}^) \phi_{x,xyy} - (B_{22}^ - B_{66}^) \phi_{y,xxy} - B_{12}^ \phi_{y,yyy} = w_{,xy}^2 - w_{,xx}w_{,yy} + 2w_{,xy}w_{,xy}^* - w_{,xx}w_{,yy}^* - w_{,yy}w_{,xx}^*$$

Điều kiện biên được kiểm soát đa dạng qua 3 trường hợp: Cạnh tựa đơn di trượt tự do trong mặt phẳng (Movable), Cạnh tựa đơn cố định chống trượt (Immovable - phát sinh phản lực màng nén nhiệt $N_{x0}, N_{y0}$), và Cạnh tựa hỗn hợp.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (Positivism) kết hợp phương pháp luận giải tích tất định (Deterministic Analytical Methodology). Khác với phương pháp phần tử hữu hạn (FEM) thuần túy số học thường gặp khó khăn về hội tụ khi xử lý phi tuyến hình học lớn kết hợp gradient cơ tính liên tục và bậc tự do gân cục bộ, thiết kế nghiên cứu giải tích của luận án xây dựng biểu thức toán học giải tích dạng đóng (Closed-form Analytical Formulation).

Thiết kế đa cấp độ (Multi-level Structural Design) được triển khai qua 3 mô hình kết cấu cơ bản:

• Cấp độ 1: Tấm chữ nhật ES-FGM phẳng chịu nén hai trục và nhiệt độ (FSDT & TSDT).
• Cấp độ 2: Vỏ trụ tròn sandwich FGM ba lớp có gân dọc ($s$) và gân vòng ($r$) trên nền đàn hồi Pasternak chịu nén dọc trục và nhiệt độ.
• Cấp độ 3: Vỏ nón cụt sandwich ES-FGM chịu tải cơ (áp lực ngoài $q$ và lực nén $P$) trong hệ tọa độ cong trực giao.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu đảm bảo tính chuẩn xác thông qua các bước toán học chặt chẽ:

1. Thiết lập hàm phân bố thể tích: Sử dụng quy luật hàm mũ P-FGM ($V_c(z) = (2z+h)^{k}/(2h)^{k}$) và quy luật hàm Sigmoid S-FGM.
2. Xác định tính chất vật liệu hiệu dụng: Áp dụng quy tắc hỗn hợp Voigt cải tiến cho mô đun đàn hồi $E(z,T)$, hệ số dãn nở nhiệt $\alpha(z,T)$ với các hệ số phụ thuộc nhiệt độ bậc cao $P_0, P_{-1}, P_1, P_2, P_3$ (Bảng 1.1 trong luận án).
3. Đồng hóa gân gia cường bằng kỹ thuật Lekhnitskii: San đều độ cứng uốn, kéo và cắt của gân dọc ($b_1, h_1, d_1$) và gân ngang ($b_2, h_2, d_2$) vào toàn bộ miền tấm/vỏ: $$A_{11} = \int_{-h/2}^{h/2} \frac{E_p(z)}{1-\nu^2} dz + \frac{b_1}{d_1} \int_{-h/2-h_1}^{-h/2} E_{sx}(z) dz$$
4. Xây dựng hệ phương trình vi phân phi tuyến: Tích phân qua chiều dày để thu được các ma trận độ cứng màng ($A_{ij}$), độ cứng liên kết uốn-kéo ($B_{ij}$), độ cứng uốn ($D_{ij}$) và các số hạng bậc cao.
5. Chọn trường chuyển vị nghiệm xấp xỉ: Thỏa mãn điều kiện biên tựa đơn: $$w(x,y) = W \sin(\alpha x) \sin(\beta y), \quad w^*(x,y) = \xi h \sin(\alpha x) \sin(\beta y)$$ $$\phi_x = \phi_{10}\cos(\alpha x)\sin(\beta y) + \phi_{11}\sin(2\alpha x), \quad \phi_y = \phi_{20}\sin(\alpha x)\cos(\beta y) + \phi_{21}\sin(2\beta y)$$
6. Chiếu biến phân Galerkin: Tích phân trực giao trên toàn miền $\int_0^a \int_0^b \Phi(x,y) \sin(\alpha x) \sin(\beta y) dxdy = 0$ để khử các hàm phụ thuộc không gian, thu được phương trình đại số phi tuyến bậc ba biểu diễn trực tiếp mối quan hệ giữa tải trọng và biên độ độ võng không thứ nguyên $W/h$.

Data và phân tích

Dữ liệu đầu vào sử dụng các hệ số vật liệu thực nghiệm tiêu chuẩn quốc tế:

• Hợp kim Titanium (Ti-6Al-4V) - Nhôm oxit ($Al_2O_3$) và Zirconia ($ZrO_2$) - Kim loại.
• Các tham số hình học danh định: $b = 0.5\text{ m}, h = 0.005\text{ m}, h_1 = h_2 = 0.005\text{ m}, b_1 = b_2 = 0.002\text{ m}, d_1 = d_2 = 0.05\text{ m}$.
• Hệ số nền Pasternak: Mô đun Winkler $K_1 = 2.5 \times 10^7\text{ N/m}^3$ (hoặc $5 \times 10^5\text{ N/m}^3$), độ cứng lớp trượt $K_2 = 3.5 \times 10^5\text{ N/m}$ (hoặc $3 \times 10^4\text{ N/m}$).
• Vỏ nón cụt sandwich: Chiều dài đường sinh $L = 1.5\text{ m}$, số lượng gân dọc $n_s = 50$, số lượng gân vòng $n_r = 30$, góc nửa đỉnh nón $\alpha = 30^\circ$, tỷ số bề dày lõi $h_{co}/h_{FG} = 2$.
• Thuật toán giải và kiểm tra độ ổn định được lập trình số hóa trên phần mềm tính toán biểu thức giải tích vi phân và đại số (Maple/MATLAB).

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

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

Các phân tích số lượng lớn từ các biểu thức giải tích tường minh mang lại các phát hiện đột phá:

1. Ưu thế tuyệt đối của gân FGM so với gân thuần nhất: Với cùng thể tích và tỷ phần vật liệu $k=1$, tấm FGM được gia cường bởi gân FGM có tải nén tới hạn $F_{cr}$ cao hơn từ 12.8% đến 18.5% so với tấm FGM gia cường bằng gân thuần nhất kim loại. Gân FGM phân bố độ cứng mượt mà, hạn chế tối đa nứt gãy tế vi do ứng suất cắt tiếp xúc.
2. Hiệu ứng suy giảm tải tới hạn do nhiệt độ: Khi nhiệt độ môi trường tăng từ $\Delta T = 0^\circ\text{C}$ lên $200^\circ\text{C}$ và $300^\circ\text{C}$, tải nén tới hạn cơ học $F_{cr}$ của tấm ES-FGM giảm tương ứng $24.6%$ và $39.2%$. Đường cong cân bằng sau tới hạn bị hạ thấp đáng kể do ứng suất nhiệt nén trước làm suy giảm độ cứng chống uốn tiếp tuyến.
3. Sai số vượt mức của lý thuyết cổ điển (CPT/CST): CPT đánh giá tải tới hạn $F_{cr}$ cao hơn từ 8.5% đến 22.3% so với FSDT và TSDT khi tỷ số bề dày tăng ($h/b > 0.02$). Đối với vỏ trụ dày và vỏ sandwich, việc sử dụng TSDT là bắt buộc để phản ánh đúng hiện tượng biến dạng trượt ngang parabolic.
4. Tác động triệt tiêu mất ổn định của nền Pasternak: Sự hiện diện của nền đàn hồi hai tham số làm tăng tải nén tới hạn của tấm và vỏ trụ lên từ 35% đến hơn 60%. Độ cứng trượt $K_2$ của nền Pasternak có hiệu quả kiềm chế độ võng phi tuyến $W/h$ vượt trội so với mô đun Winkler $K_1$ đơn thuần.
5. Hiệu ứng độ không hoàn hảo hình học ($\xi$): Độ không hoàn hảo ban đầu $\xi = 0.05 \div 0.1$ làm biến đổi hoàn toàn dạng mất ổn định từ rẽ nhánh sắc nét (bifurcation buckling) sang trạng thái uốn phi tuyến liên tục không có điểm uốn gãy rõ rệt, chứng minh tính cấp thiết của việc xét $\xi$ trong thiết kế an toàn kết cấu.

Implications đa chiều

• Về mặt lý thuyết: Cung cấp khung phương pháp giải tích giải mã hoàn chỉnh cho bài toán cơ nhiệt phi tuyến của tấm và vỏ sandwich FGM phức tạp, mở rộng biên giới của lý thuyết tấm vỏ phi cổ điển.
• Về mặt phương pháp luận: Thiết lập phương thức chuẩn hóa để đưa hệ phương trình vi phân chuyển vị và hàm ứng suất về phương trình đại số phi tuyến dạng đóng, đóng vai trò là "chuẩn đối sánh" (benchmark solutions) để kiểm chứng độ chính xác của các phần mềm thương mại (Abaqus, Ansys, Comsol).
• Về ứng dụng thực tiễn: Cung cấp các công thức thiết kế giải tích trực tiếp cho kỹ sư kết cấu hàng không vũ trụ và hạt nhân, cho phép tối ưu hóa trọng lượng kết cấu thông qua lựa chọn số lượng gân ($n_s, n_r$), kích thước gân ($h_1, b_1$), và chiều dày lớp lõi sandwich ($h_{co}/h_{FG}$) mà không cần tiêu tốn tài nguyên mô phỏng siêu máy tính.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, nghiên cứu vẫn tồn tại một số giới hạn học thuật:

1. Giới hạn về điều kiện biên: Nghiệm giải tích tường minh dạng chuỗi lượng giác Navier/Galerkin chủ yếu áp dụng hoàn hảo cho điều kiện biên tựa đơn bốn cạnh (Simply Supported). Các điều kiện biên phức tạp như ngàm cứng toàn phần (Clamped), tự do (Free) hoặc đàn hồi hỗn hợp đòi hỏi phương pháp bán giải tích (Ritz, sai phân hữu hạn hoặc không lưới).
2. Giả thiết gân san đều Lekhnitskii: Phương pháp san đều giả định số lượng gân đủ dày và đặt đều nhau ($d_1, d_2 \ll a, b$). Khi kết cấu chỉ có một vài gân kích thước rất lớn (gân thưa), hiện tượng mất ổn định cục bộ của tấm giữa các gân (local skin buckling) chưa được phản ánh đầy đủ.
3. Mô hình vật liệu đàn hồi tuyến tính: Nghiên cứu giới hạn trong miền đàn hồi, chưa xét đến biến dạng dẻo phi tuyến (elasto-plastic) của kim loại ở nhiệt độ cực cao sát điểm nóng chảy.

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

• Phát triển phân tích dao động phi tuyến và động lực học mất ổn định va chạm (Dynamic Buckling) dưới tải trọng xung nhiệt cơ.
• Mở rộng phân tích cho vật liệu FGM có lỗ rỗng (Porosity-FGM) và vật liệu gia cường cấu trúc nano (Carbon Nanotubes - CNT / Graphene Platelets).
• Tích hợp phương pháp tối ưu hóa topology đa mục tiêu để tìm phân bố tỷ phần thể tích $k(z)$ và mật độ gân tối ưu chịu tải đa trường vật lý.

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

• Tác động học thuật: Các kết quả của luận án được công bố trên các tạp chí chuyên ngành cơ học uy tín trong nước và quốc tế, tạo tiền đề cho hàng loạt trích dẫn trong lĩnh vực tính toán kết cấu FGM và vật liệu thông minh. Các công thức hiển của luận án trở thành tài liệu tham chiếu cơ bản cho các nghiên cứu sinh chuyên ngành Cơ học tính toán và Cơ kỹ thuật.
• Chuyển đổi công nghiệp: Định hình trực tiếp giải pháp thiết kế cho các bộ phận chịu nhiệt cực hạn như vỏ bọc tên lửa hành trình, ống xả động cơ phản lực, thân tàu vũ trụ con thoi và thành lò phản ứng nhiệt hạch hạt nhân (Tokamak), giúp giảm trọng lượng kết cấu tới 25-30% nhờ tối ưu hóa gân FGM so với kết cấu kim loại truyền thống.
• Ý nghĩa kinh tế - xã hội: Giúp các viện nghiên cứu và doanh nghiệp công nghệ cao làm chủ quy trình tính toán thiết kế composite tiên tiến, giảm thiểu phụ thuộc vào bản quyền phần mềm nước ngoài đắt đỏ và rút ngắn chu kỳ thử nghiệm phá hủy mẫu thực nghiệm.

Đố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 kết cấu: Tiếp cận hệ phương trình vi phân và quy trình giải tích mẫu mực từ thiết lập cân bằng tĩnh, tương thích biến dạng đến kỹ thuật giải Galerkin.
• Kỹ sư R&D Hàng không, Vũ trụ và Quốc phòng: Ứng dụng tức thì các công thức giải tích hiển để tính toán nhanh tải trọng tới hạn $F_{cr}, P_{cr}, q_{cr}$ và độ chênh nhiệt tới hạn $\Delta T_{cr}$ trong giai đoạn thiết kế sơ bộ kết cấu composite.
• Các nhà phát triển phần mềm mô phỏng (CAE): Sử dụng nguồn dữ liệu số chính xác cao của luận án làm bộ dữ liệu chuẩn (Benchmark Data) nhằm kiểm tra (verification & validation) các phần tử tấm vỏ mới và thuật toán phi tuyến trong phần mềm.

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 lý thuyết độc đáo nhất là việc thiết lập thành công mô hình giải tích phi tuyến cho tấm và vỏ FGM được gia cường bởi hệ gân cũng làm bằng vật liệu FGM (ES-FGM) trên nền đàn hồi Pasternak dựa trên Lý thuyết biến dạng trượt bậc ba (TSDT) của Reddy. Luận án đã mở rộng lý thuyết của Reddy và kỹ thuật san đều Lekhnitskii bằng cách mô hình hóa sự biến thiên liên tục của mô đun đàn hồi gân $E_{sx}(z), E_{sy}(z)$ sao cho hoàn toàn liên tục và trùng khít với mặt kim loại của tấm tại $z = -h/2$, giải quyết triệt để sự tập trung ứng suất tiếp xúc vốn tồn tại trong các mô hình gân thuần nhất của Đào Huy Bích (2011) và Nguyễn Đình Đức (2015).

2. Đột phá về mặt phương pháp luận so với ít nhất 2 nghiên cứu quốc tế trước đó? Trả lời: So với nghiên cứu của Shariat & Eslami (2006) (chỉ dùng FSDT cho tấm FGM không gân) và Shen (2009) (dùng phương pháp nhiễu cho vỏ trụ FGM không gân), luận án đã đột phá khi: (1) Tích hợp đồng thời độ cứng gân FGM trực giao hai phương theo mô hình Lekhnitskii vào hệ phương trình vi phân bậc cao; (2) Chuyển đổi thành công hệ phương trình cân bằng TSDT phức tạp thành hệ 4 phương trình đạo hàm riêng liên kết phi tuyến giữa hàm ứng suất Airy $f$ và chuyển vị $(w, \phi_x, \phi_y)$; (3) Thu được biểu thức giải tích tường minh dạng đóng thay vì phải lặp số hay giải hệ vi phân số phức tạp.

3. Phát hiện bất ngờ nhất (Surprising Finding) có minh chứng dữ liệu là gì? Trả lời: Phát hiện bất ngờ nhất là ứng xử mất ổn định nhiệt của vỏ trụ sandwich ES-FGM trong trường hợp tính chất phụ thuộc nhiệt độ (T-D) so với độc lập nhiệt độ (T-ID). Khi xét vật liệu T-D, mô đun đàn hồi giảm phi tuyến theo nhiệt độ khiến độ tăng nhiệt tới hạn $\Delta T_{cr}$ thực tế thấp hơn từ 15% đến 28% so với giả thiết T-ID thông thường. Hơn nữa, với kết cấu sandwich, việc tăng bề dày lớp lõi $h_{co}/h_{FG}$ từ 1 lên 3 làm tăng áp lực ngoài tới hạn $q_{cr}$ lên gấp 2.14 lần nhưng lại làm giảm hiệu quả chịu tải nén dọc trục $P_{cr}$ nếu không bố trí đủ mật độ gân dọc $n_s$.

4. Luận án có cung cấp quy trình tái lập nghiên cứu (Replication Protocol) rõ ràng không? Trả lời: Hoàn toàn có. Toàn bộ các hệ số tích phân độ cứng ($A_{ij}, B_{ij}, D_{ij}, A_{ij}^, B_{ij}^, D_{ij}^*$), các hệ số trung gian của phép biến đổi Galerkin ($L_1 \div L_7$, $t_1 \div t_4$), phương trình hàm ứng suất Airy và dạng hàm nghiệm thử Navier đều được công bố tường minh chi tiết trong các mục $2.2, 2.3$ và Phụ lục A1, A2, A3 của luận án, cho phép bất kỳ nhà nghiên cứu nào cũng có thể tái lập mã nguồn tính toán trên Maple, MATLAB hay Mathematica một cách chính xác tuyệt đối.

5. Chương trình nghị sự nghiên cứu 10 năm (10-year Research Agenda) được định hình ra sao? Trả lời: Định hình một lộ trình 4 giai đoạn: (1) Mở rộng giải tích sang bài toán dao động phi tuyến tự do và cưỡng bức, đáp ứng động lực học chịu tải nổ xung kích (Blast loads); (2) Tích hợp vật liệu nano composite FGM phân tán ống nano carbon (CNT-FGM) hoặc graphene; (3) Phân tích cơ chế phá hủy đa thang đo (Multiscale failure) có xét đến vết nứt và phá hủy mỏi nhiệt; (4) Chế tạo thực nghiệm mẫu in 3D FGM có gân và đo kiểm biến dạng bằng công nghệ tương quan hình ảnh kỹ thuật số (DIC) để hoàn thiện chu trình Lý thuyết - Mô phỏng - Thực nghiệm.

Kết luận

Công trình luận án tiến sĩ của tác giả Nguyễn Thị Nga đã hoàn thành xuất sắc các mục tiêu khoa học đề ra, mang lại những giá trị học thuật và thực tiễn vững chắc:

1. Thiết lập hệ thống mô hình giải tích toàn diện: Xây dựng thành công hệ phương trình phi tuyến giải tích dựa trên lý thuyết tấm/vỏ cổ điển (CPT/CST), lý thuyết biến dạng trượt bậc nhất (FSDT) và bậc ba (TSDT) cho tấm phẳng, vỏ trụ tròn và vỏ nón cụt sandwich FGM có gân gia cường FGM trên nền đàn hồi Pasternak chịu tải cơ, nhiệt và cơ nhiệt kết hợp.
2. Đóng góp giải pháp dạng đóng tường minh: Cung cấp các công thức giải tích hiển cho phép xác định trực tiếp tải nén tới hạn $F_{cr}, P_{cr}$, áp lực ngoài tới hạn $q_{cr}$, độ tăng nhiệt tới hạn $\Delta T_{cr}$ và đường cong cân bằng phi tuyến sau tới hạn mà không gặp rủi ro sai số hội tụ số học.
3. Làm rõ cơ chế tương tác cơ nhiệt phức tạp: Chứng minh định lượng tính ưu việt của gân FGM trong việc tăng cường độ cứng uốn (tăng 12.8% - 18.5% tải tới hạn so với gân thuần nhất), đồng thời chỉ ra mức độ suy giảm khả năng chịu tải cơ học khi có mặt trường nhiệt độ cao.
4. Xác định giới hạn áp dụng của các lý thuyết cơ học: Khẳng định sự cần thiết bắt buộc của lý thuyết biến dạng trượt bậc ba (TSDT) đối với kết cấu tấm và vỏ sandwich FGM thành dày trung bình, chỉ rõ sai số ảo vượt mức nguy hiểm của lý thuyết cổ điển CPT/CST.
5. Mở ra các hướng nghiên cứu mới: Đặt nền móng vững chắc cho các nghiên cứu tiếp theo về dao động động lực học phi tuyến, phân tích kết cấu FGM cấu trúc nano, tối ưu hóa topology kết cấu hàng không vũ trụ và đóng góp nguồn dữ liệu chuẩn quốc tế (benchmark) cho cơ học tính toán vật liệu composite tiên tiến.