Tổng quan về luận án
Nghiên cứu ổn định tĩnh và ứng xử sau tới hạn của các kết cấu thành mỏng và dày vừa làm bằng vật liệu cơ tính biến thiên (Functionally Graded Materials - FGM) có gân gia cường là một bài toán then chốt trong cơ học vật rắn biến dạng hiện đại. Luận án tiến sĩ cơ học của tác giả Nguyễn Thị Nga (2018), bảo vệ tại Trường Đại học Khoa học Tự nhiên – Đại học Quốc gia Hà Nội dưới sự hướng dẫn của cố GS. Đào Văn Dũng và PGS. Vũ Đỗ Long, mang tiê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" (Mã số: 62440107). Công trình này đại diện cho bước phát triển tiên phong trong việc xây dựng khung lý thuyết giải tích giải quyết bài toán mất ổn định tĩnh phi tuyến (non-linear static buckling and post-buckling) của các kết cấu FGM phức tạp (tấm chữ nhật, vỏ trụ tròn sandwich, vỏ nón cụt sandwich) có gân gia cường lệch tâm FGM (ES-FGM) đặt trên nền đàn hồi Pasternak chịu tải trọng cơ học, tải nhiệt và tải cơ - nhiệt kết hợp.
+------------------------------------------------------------------------------------+
| KẾT CẤU ES-FGM |
| (Tấm chữ nhật / Vỏ trụ tròn Sandwich / Vỏ nón cụt Sandwich + Gân gia cường FGM) |
+------------------------------------------------------------------------------------+
|
v
+------------------------------------------------------------------------------------+
| MÔ HÌNH LÝ THUYẾT & VẬT LIỆU |
| - Quy luật thể tích P-FGM / S-FGM (phụ thuộc T-D & không phụ thuộc nhiệt độ T-ID) |
| - Lý thuyết biến dạng trượt bậc nhất (FSDT) & bậc ba (TSDT) + Hình học Von Karman |
| - Kỹ thuật san đều tác dụng gân Lekhnitskii + Nền đàn hồi Pasternak (K1, K2) |
+------------------------------------------------------------------------------------+
|
v
+------------------------------------------------------------------------------------+
| PHƯƠNG PHÁP GIẢI TÍCH GALERKIN |
| - Thiết lập hàm ứng suất Airy f(x, y) & phương trình tương thích biến dạng |
| - Xây dựng nghiệm hiển giải tích cho tải tới hạn (Fxcr, Pcr, qcr, ΔTcr) |
| - Thiết lập phương trình phi tuyến đường cong tải - độ võng sau tới hạn |
+------------------------------------------------------------------------------------+
Research Gap trong y văn cơ học kết cấu
Trước thời điểm công trình này được công bố, các nghiên cứu quốc tế và trong nước về ổn định của kết cấu FGM bộc lộ các khoảng trống học thuật (research gaps) mang tính hệ thống:
- Giới hạn của lý thuyết tấm/vỏ cổ điển (CPT/CST): Phần lớn các công trình của Eslami và cộng sự (2005, 2007), Sofiyev (2011), Shen (2009) chủ yếu dựa trên giả thuyết Kirchhoff-Love (CPT/CST), vốn bỏ qua biến dạng trượt ngang ($\gamma_{xz} = \gamma_{yz} = 0$). Cách tiếp cận này dẫn tới sự sai lệch nghiêm trọng và đánh giá quá cao (overestimate) tải trọng tới hạn khi áp dụng cho các kết cấu dày vừa và kết cấu sandwich chịu gradient nhiệt độ cao.
- Khoảng trống về gân gia cường FGM: Các nghiên cứu có xét đến gân gia cường lệch tâm (như Đào Huy Bích và cộng sự, 2011; Nguyễn Đình Đức và Phạm Hồng Công, 2013) hầu hết giả định gân làm bằng vật liệu thuần nhất (homogeneous stiffeners - gân thuần gốm hoặc thuần kim loại). Giả định này gây ra sự bất liên tục về ứng suất tiếp xúc và gradient nhiệt tại bề mặt tiếp giáp giữa tấm/vỏ và chân gân, chưa phản ánh đúng bản chất của công nghệ chế tạo vật liệu biến tính liên tục.
- Thiếu vắng nghiệm giải tích dạng tường minh cho kết cấu sandwich vỏ nón và vỏ trụ phức tạp: Chưa có một khung giải tích hoàn chỉnh nào tích hợp đồng thời: lý thuyết biến dạng trượt bậc ba (TSDT), tính phi tuyến hình học Von Kármán, kỹ thuật san đều tác dụng gân Lekhnitskii, mô hình nền đàn hồi hai thông số Pasternak và sự phụ thuộc phi tuyến của cơ tính vào nhiệt độ (T-D) cho kết cấu sandwich FGM có gân gia cường FGM.
Câu hỏi nghiên cứu và hệ thống giả thuyết
Luận án tập trung giải quyết 4 câu hỏi nghiên cứu cốt lõi:
- RQ1: Bản chất của tương tác giữa biến dạng trượt ngang bậc cao và gradient nhiệt độ phi tuyến ảnh hưởng như thế nào đến điểm rẽ nhánh ổn định và đáp ứng sau tới hạn của tấm và vỏ ES-FGM?
- RQ2: Việc sử dụng hệ gân gia cường biến tính FGM so với gân thuần nhất tạo ra sự thay đổi định lượng như thế nào đối với khả năng chịu tải tới hạn cơ học và nhiệt?
- RQ3: Mô hình cấu trúc sandwich (lõi FGM phủ kim loại/gốm hoặc lõi thuần nhất phủ FGM) cùng các thông số tỷ lệ chiều dày lõi ($h_{co}/h_{FG}$) điều biến ứng xử mất ổn định của vỏ trụ và vỏ nón cụt ra sao dưới các điều kiện biên khác nhau?
- RQ4: Mức độ nhạy cảm của tải tới hạn trước độ không hoàn hảo hình học ban đầu ($\xi$) và tương tác phản lực nền đàn hồi Pasternak ($K_1, K_2$) được mô tả chính xác bằng biểu thức giải tích nào?
Hệ thống giả thuyết khoa học được kiểm chứng bao gồm:
- H1: Lý thuyết TSDT cho kết quả tải nén và tải nhiệt tới hạn thấp hơn đáng kể (từ 5% đến 18%) so với CPT/CST đối với kết cấu có tỷ số $h/b > 0.05$ hoặc $h/R > 0.01$, phản ánh chính xác hiệu ứng suy giảm độ cứng chống cắt ngang.
- H2: Gân gia cường FGM với sự biến thiên cơ tính trơn liên tục loại bỏ hiện tượng tập trung ứng suất tại liên kết, đồng thời nâng cao tải tới hạn của kết cấu lên từ 1.2 đến 2.5 lần so với kết cấu không gân và vượt trội rõ rệt so với gân thuần nhất cùng khối lượng.
- H3: Độ không hoàn hảo hình học ban đầu ($\xi \in [0, 0.3]$) làm suy giảm đường cong cân bằng sau tới hạn và làm mất tính rẽ nhánh đột ngột, chuyển kết cấu sang trạng thái biến dạng uốn - nén phi tuyến liên tục.
Khung lý thuyết và phạm vi nghiên cứu
Nghiên cứu được xây dựng trên nền tảng vững chắc của:
- Lý thuyết đàn hồi phi tuyến và động học Von Kármán;
- Lý thuyết biến dạng trượt bậc nhất FSDT (Reissner-Mindlin) và bậc ba TSDT (Reddy);
- Tiêu chuẩn ổn định rẽ nhánh Euler - Poincaré;
- Kỹ thuật san đều tác dụng gân Lekhnitskii cải tiến cho vật liệu FGM;
- Mô hình tương tác cơ học nền đàn hồi Pasternak hai tham số ($K_1$ mô đun Winkler, $K_2$ lớp cắt trượt).
Phạm vi nghiên cứu bao gồm: tấm FGM chữ nhật kích thước $a \times b \times h$; vỏ trụ tròn sandwich ES-FGM bán kính $R$, chiều dài $L$; và vỏ nón cụt sandwich ES-FGM góc nửa đỉnh $\alpha$, bán kính đáy bé $r_1$, bán kính đáy lớn $r_2$. Hệ vật liệu khảo sát là gốm kim loại ($ZrO_2/Ti-6Al-4V$, $Si_3N_4/SUS304$, $Al_2O_3/Al$) với các hàm phân bố tỷ phần thể tích lũy thừa (P-FGM) và Sigmoid (S-FGM), có xét đến tính chất độc lập nhiệt độ (T-ID) và phụ thuộc nhiệt độ (T-D) thông qua hàm nhiệt động phi tuyến $P(T) = P_0(P_{-1}T^{-1} + 1 + P_1T + P_2T^2 + P_3T^3)$.
Literature Review và Positioning
Tổng quan các dòng nghiên cứu lớn trong cơ học FGM
Lịch sử nghiên cứu ứng xử mất ổn định của kết cấu composite biến tính FGM kể từ khi vật liệu này được phát minh tại Viện Sendai (Nhật Bản, 1984) được chia thành ba trường phái giải thuật chính:
TIẾN TRÌNH NGHIÊN CỨU ỔN ĐỊNH KẾT CẤU FGM
|
+--------------------------------------------+--------------------------------------------+
| | |
v v v
TRƯỜNG PHÁI GIẢI TÍCH CỔ ĐIỂN TRƯỜNG PHÁI BIẾN DẠNG TRƯỢT CAO CẤP TRƯỜNG PHÁI KẾT CẤU GIA CƯỜNG GÂN
- Eslami & cs (2005, 2007) [CPT] - Reddy (2000), Zenkour (2005) [TSDT] - Najafizadeh & cs (2009) [Gân FGM]
- Sofiyev (2011, 2017) [CST/Galerkin] - Shen (2009, 2014) [Perturbation/TSDT] - Đào Huy Bích & cs (2011) [Gân thuần nhất]
- Brush & Almroth (1975), Seide (1962) - Hoàng Văn Tùng (2011, 2014) [FSDT/TSDT] - Nguyễn Đình Đức (2015, 2016) [Gân thuần/TSDT]
| | |
+--------------------------------------------+--------------------------------------------+
|
v
LUẬN ÁN NGUYỄN THỊ NGA (2018)
- Tích hợp TSDT/FSDT + Phi tuyến Von Karman + Gân FGM lệch tâm
- Vỏ trụ/nón cụt Sandwich FGM + Nền Pasternak + Tính chất T-D
- Nghiệm hiển giải tích Galerkin đóng kín cho tải tới hạn & sau tới hạn
- Dòng nghiên cứu tấm và vỏ FGM đẳng hướng không gân: Eslami cùng các cộng sự (2005, 2007) mở đầu với việc phân tích ổn định tuyến tính của tấm chữ nhật và vỏ nón FGM chịu tải cơ - nhiệt dựa trên CPT và lý thuyết động học Sander. Shariat và Eslami (2007) mở rộng sang lý thuyết biến dạng trượt bậc cao nhưng chỉ giới hạn cho tấm phẳng đơn lớp, không gân và không nền đàn hồi. Shen (2009, 2014) phát triển phương pháp giải tích nhiễu (perturbation technique) khảo sát đáp ứng sau tới hạn của vỏ trụ và panel trụ FGM có xét đến tính chất vật liệu phụ thuộc nhiệt độ $P(T)$ và lớp áp điện, tuy nhiên cách tiếp cận này chỉ áp dụng cho kết cấu không có gân gia cường.
- Dòng nghiên cứu tấm và vỏ sandwich FGM: Zenkour (2005, 2006) cùng Zenkour và Alghamdi (2008) chứng minh việc chèn lớp lõi đồng nhất kẹp giữa hai lớp phủ FGM (hoặc lõi FGM kẹp giữa hai lớp kim loại/gốm) tối ưu hóa khả năng chống uốn mà hầu như không tăng đáng kể trọng lượng bản thân. Sofiyev và cộng sự (2011, 2017) tập trung giải quyết bài toán dao động và mất ổn định của vỏ nón cụt và vỏ trụ sandwich FGM chịu nén dọc trục và áp lực ngoài bằng lý thuyết CST và FSDT.
- Dòng nghiên cứu kết cấu FGM có gân gia cường: Nhận thức được hiệu quả tối ưu hóa kết cấu của gân, Najafizadeh và cộng sự (2009) lần đầu tiên đề xuất mô hình gân FGM cho vỏ trụ dựa trên CST. Tại Việt Nam, Đào Huy Bích và cộng sự (2011) tiên phong phát triển phương pháp giải tích cho tấm và vỏ thoải ES-FGM dựa trên CST với gân thuần nhất. Tiếp đó, Đào Văn Dũng và cộng sự (2013, 2017), Nguyễn Đình Đức và cộng sự (2015, 2016) mở rộng sang lý thuyết FSDT/TSDT cho tấm và vỏ trụ mỏng, nhưng đa số vẫn xử lý gân như một thành phần thuần nhất (gân kim loại hoặc gân gốm).
Tranh luận học thuật và định vị nghiên cứu
Một cuộc tranh luận học thuật lớn tồn tại giữa hai quan điểm:
- Quan điểm 1 (Simplification School - CPT/CST): Cho rằng đối với kết cấu mỏng ($h/a < 0.02$), ảnh hưởng của biến dạng trượt ngang là không đáng kể, việc áp dụng CPT giúp đơn giản hóa hệ phương trình vi phân đạo hàm riêng phi tuyến từ 4 phương trình bậc cao xuống 2 phương trình, tiết kiệm thời gian tính toán và dễ tìm nghiệm giải tích.
- Quan điểm 2 (Higher-order Shear School - FSDT/TSDT): Chứng minh rằng trong môi trường tải nhiệt độ tăng cao kết hợp với vật liệu có gradient cơ tính lớn theo chiều dày như FGM, sự tập trung ứng suất cắt ngang $\tau_{xz}, \tau_{yz}$ và hiệu ứng vồng phi tuyến là rất lớn. Việc bỏ qua biến dạng trượt sẽ dẫn đến việc dự báo sai lệch trạng thái mất ổn định sau tới hạn, đặc biệt là hiện tượng nhảy snap-through.
Định vị của luận án: Luận án của NCS. Nguyễn Thị Nga định vị chính xác tại giao điểm phức tạp nhất: vừa mở rộng tính chặt chẽ của lý thuyết TSDT cho các kết cấu có gân gia cường FGM lệch tâm, vừa vượt qua rào cản toán học để thiết lập các công thức giải tích đóng kín (closed-form solutions) thay vì dựa hoàn toàn vào phương pháp số xấp xỉ (FEM).
So sánh đối chuẩn quốc tế (International Benchmarking)
- So sánh với nghiên cứu kinh điển của Shariat và Eslami (2007) (International Journal of Mechanical Sciences): Khi suy biến mô hình tấm FSDT/TSDT của luận án về trường hợp tấm FGM không gân, không nền, kết quả tải nén tới hạn không thứ nguyên $F_x^* = F_x a^2/(\pi^2 D_m)$ đạt độ hội tụ và sai số tuyệt đối dưới 0.35%.
- So sánh với nghiên cứu của Sofiyev (2011) (Composite Structures) và Seide (1962): Đối với bài toán vỏ nón cụt đẳng hướng và vỏ nón sandwich FGM chịu tải nén dọc trục và áp lực ngoài, biểu thức giải tích của luận án cho kết quả tải tới hạn $P_{cr}$ và $q_{cr}$ trùng khớp hoàn hảo với các kết quả chuẩn mực của Seide (1962), Brush và Almroth (1975) khi suy biến về vỏ thuần nhất, và sai lệch dưới 1.2% so với nghiệm bán giải tích của Sofiyev (2011).
Đóng góp lý thuyết và khung phân tích
+-------------------------------------------------------------------------------------------------------------------------+
| KHUNG ĐÓNG GÓP LÝ THUYẾT CỦA LUẬN ÁN |
+-------------------------------------------------------------------------------------------------------------------------+
| 1. MỞ RỘNG LÝ THUYẾT BIẾN DẠNG TRƯỢT BẬC BA (TSDT) CHO KẾT CẤU ES-FGM |
| - Mô tả chính xác trường chuyển vị: u(x,y,z) = u0 + z*phi_x - c1*z^3*(phi_x + w,x) |
| - Triệt tiêu hoàn toàn ứng suất cắt ngang trên hai bề mặt tự do (z = ±h/2) mà không cần hệ số hiệu chỉnh trượt ks |
+-------------------------------------------------------------------------------------------------------------------------+
| 2. MÔ HÌNH HÓA TOÀN DIỆN HỆ GÂN GIA CƯỜNG FGM LỆCH TÂM (ES-FGM) |
| - Đảm bảo tính liên tục hoàn hảo về mô đun đàn hồi: Ep(-h/2) = Esx(-h/2) = Esy(-h/2) = Em |
| - Thiết lập ma trận độ cứng tương đương tích hợp độ lệch tâm của gân dọc và gân vòng/ngang |
+-------------------------------------------------------------------------------------------------------------------------+
| 3. THIẾT LẬP HỆ PHƯƠNG TRÌNH PHI TUYẾN 4 BIẾN VỚI HÀM ỨNG SUẤT AIRY & GALERKIN |
| - Khử chuyển vị trong mặt phẳng, giữ lại 4 hàm ẩn: w(x, y), phi_x(x, y), phi_y(x, y), f(x, y) |
| - Xây dựng biểu thức nghiệm giải tích hiển xác định trực tiếp tải tới hạn và hàm liên hệ tải - độ võng sau tới hạn |
+-------------------------------------------------------------------------------------------------------------------------+
Đóng góp cho lý thuyết cơ học kết cấu
-
Hoàn thiện lý thuyết biến dạng trượt bậc cao cho kết cấu có gân FGM: Luận án mở rộng thành công lý thuyết TSDT của Reddy cho tấm và vỏ sandwich có gân gia cường FGM. Trường chuyển vị được biểu diễn chính xác:
$$u(x, y, z) = u_0(x, y) + z\phi_x(x, y) - c_1 z^3 \left(\phi_x + \frac{\partial w}{\partial x}\right)$$
$$v(x, y, z) = v_0(x, y) + z\phi_y(x, y) - c_1 z^3 \left(\phi_y + \frac{\partial w}{\partial y}\right)$$
$$w(x, y, z) = w(x, y)$$
với $c_1 = 4/(3h^2)$. Cấu trúc này triệt tiêu hoàn toàn ứng suất tiếp tuyến $\tau_{xz}, \tau_{yz}$ tại hai bề mặt tự do $z = \pm h/2$, loại bỏ sự phụ thuộc vào hệ số điều chỉnh trượt $k_s$ vốn gây tranh cãi trong lý thuyết FSDT.
-
Xây dựng quy luật cơ tính gân biến thiên liên tục (Continuous FGM Stiffener Law): Khác với các công bố trước đây coi gân là thuần gốm hoặc thuần kim loại, luận án thiết lập hàm biến thiên mô đun đàn hồi và hệ số giãn nở nhiệt của hệ gân dọc ($s$) và gân ngang/vòng ($r$):
$$E_{sx}(z) = E_m + (E_c - E_m)\left(\frac{2z + h}{2h_1}\right)^{k_2}, \quad -\frac{h}{2} - h_1 \le z \le -\frac{h}{2}$$
$$E_{sy}(z) = E_m + (E_c - E_m)\left(\frac{2z + h}{2h_2}\right)^{k_3}, \quad -\frac{h}{2} - h_2 \le z \le -\frac{h}{2}$$
Công thức này bảo toàn điều kiện liên tục tuyệt đối về vật liệu tại mặt tiếp giáp $z = -h/2$ ($E_p(-h/2) = E_{sx}(-h/2) = E_{sy}(-h/2) = E_m$), loại bỏ kỳ dị toán học và nứt gãy ứng suất tại liên kết chân gân.
Khung phân tích độc đáo và tích hợp lý thuyết
Khung phân tích của luận án tích hợp thành công 4 trụ cột cơ học lý thuyết:
- Động học phi tuyến Von Kármán cho biến dạng hữu hạn (finite displacement);
- Kỹ thuật san đều Lekhnitskii mở rộng tích hợp độ cứng uốn, kéo và xoắn của gân FGM lệch tâm;
- Phương pháp hàm ứng suất Airy $f(x, y)$ thỏa mãn chính xác phương trình cân bằng trong mặt phẳng;
- Thuật toán tích phân trực giao Galerkin chiếu hệ phương trình vi phân phi tuyến bậc cao lên không gian nghiệm lượng giác.
Hệ phương trình chủ đạo phi tuyến 4 ẩn $(w, \phi_x, \phi_y, f)$ cho tấm ES-FGM không hoàn hảo trên nền Pasternak được thiết lập dưới dạng tổng quát:
$$B_{21}^* f_{,xxxx} + (B_{11}^* + B_{22}^* + 2B_{66}^) f_{,xxyy} + B_{12}^ f_{,yyyy} - D_{11}^* \phi_{x,xxx} - (D_{12}^* + 2D_{66}^) \phi_{y,xxy} - (D_{21}^ + 2D_{66}^) \phi_{x,xyy} - D_{22}^ \phi_{y,yyy} - \mathcal{L}(w + w^, f) + K_1 w - K_2 \nabla^2 w = 0$$
$$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} + \frac{1}{2}\mathcal{L}(w, w + 2w^) = 0$$
trong đó toán tử vi phân phi tuyến $\mathcal{L}(p, q) = p_{,xx}q_{,yy} - 2p_{,xy}q_{,xy} + p_{,yy}q_{,xx}$, và các hệ số $A_{ij}^, B_{ij}^, D_{ij}^$ là các đại lượng tích phân phức hợp phản ánh cơ tính hiệu dụng của tấm, gân và trường nhiệt độ.
Phương pháp nghiên cứu tiên tiến
+----------------------------------------------------------------------------------------------------+
| QUY TRÌNH PHƯƠNG PHÁP NGHIÊN CỨU |
+----------------------------------------------------------------------------------------------------+
| BƯỚC 1: MÔ HÌNH HÓA ĐỘNG HỌC & CẤU TRÚC VẬT LIỆU |
| - Thiết lập trường biến dạng phi tuyến Von Karman & cấu hàm P-FGM/S-FGM phụ thuộc nhiệt độ T-D |
| - San đều hệ gân dọc/vòng FGM theo kỹ thuật Lekhnitskii cải tiến |
+----------------------------------------------------------------------------------------------------+
|
v
+----------------------------------------------------------------------------------------------------+
| BƯỚC 2: THIẾT LẬP HỆ PHƯƠNG TRÌNH VI PHÂN CHỦ ĐẠO |
| - Đưa vào hàm ứng suất Airy f(x, y) để triệt tiêu phương trình cân bằng trong mặt phẳng |
| - Biểu diễn hệ phương trình vi phân phi tuyến 4 ẩn: {w, phi_x, phi_y, f} |
+----------------------------------------------------------------------------------------------------+
|
v
+----------------------------------------------------------------------------------------------------+
| BƯỚC 3: GIẢI TÍCH GALERKIN & THIẾT LẬP NGHIỆM HIỂN |
| - Lựa chọn hệ hàm thử lượng giác thỏa mãn điều kiện biên tựa đơn |
| - Thực hiện tích phân Galerkin đóng kín, suy dẫn trực tiếp công thức tải tới hạn & sau tới hạn |
+----------------------------------------------------------------------------------------------------+
|
v
+----------------------------------------------------------------------------------------------------+
| BƯỚC 4: THỰC THI THUẬT TOÁN SỐ & ĐỐI SÁNH ĐỘ TIN CẬY |
| - Lập trình tính toán tự động bằng Maple / MATLAB |
| - Đối chuẩn đa tầng (CPT vs FSDT vs TSDT, so sánh Shariat-Eslami, Sofiyev, Brush-Almroth) |
+----------------------------------------------------------------------------------------------------+
Triết lý và thiết kế nghiên cứu
Nghiên cứu tuân thủ chặt chẽ triết lý diễn dịch giải tích thực chứng (deductive analytical positivism). Nghiên cứu không sử dụng phương pháp phần tử hữu hạn (FEM) thương mại như một "hộp đen" mà xây dựng phương pháp giải tích tường minh (analytical closed-form formulation) kết hợp với thuật toán lập trình biểu tượng (symbolic computation) trên môi trường Maple và MATLAB.
Thiết kế nghiên cứu đa tầng bao gồm:
- Tầng cấu trúc vi mô (Micro-mechanics level): Tính toán cơ tính hiệu dụng $E(z, T), \alpha(z, T), \nu$ theo quy tắc hỗn hợp biến tính phụ thuộc nhiệt độ.
- Tầng mặt cắt tương đương (Equivalent cross-section level): Tích phân số học theo chiều dày $z$ tích hợp độ cứng của tấm và hệ gân trực giao san đều.
- Tầng kết cấu vĩ mô (Macro-structural level): Thiết lập bài toán biên phi tuyến dưới các điều kiện biên tựa đơn tự do (movable edges), tựa đơn cố định (immovable edges) và biên hỗn hợp.
Thuật toán giải tích và quy trình Galerkin
Hệ hàm thử cho độ võng $w(x, y)$, độ không hoàn hảo ban đầu $w^(x, y)$, các góc xoay $\phi_x, \phi_y$ và hàm ứng suất $f(x, y)$ được lựa chọn nghiêm ngặt thỏa mãn điều kiện biê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(x, y) = \Phi_{10} \cos(\alpha x) \sin(\beta y) + \Phi_{11} \sin(2\alpha x)$$
$$\phi_y(x, y) = \Phi_{20} \sin(\alpha x) \cos(\beta y) + \Phi_{21} \sin(2\beta y)$$
$$f(x, y) = f_1 \cos(2\alpha x) + f_2 \cos(2\beta y) + f_3 \sin(\alpha x) \sin(\beta y) - \frac{1}{2}N_{x0}y^2 - \frac{1}{2}N_{y0}x^2$$
trong đó $\alpha = m\pi/a$, $\beta = n\pi/b$, với $(m, n)$ là số nửa sóng biến dạng.
Bằng cách áp dụng phép chiếu trực giao Galerkin:
$$\int_0^a \int_0^b \Phi_i(x, y) \cdot \mathcal{E}(w, \phi_x, \phi_y, f) , dx dy = 0$$
luận án thu được hệ phương trình đại số phi tuyến dạng tường minh kết nối trực tiếp biên độ võng không thứ nguyên $\bar{W} = W/h$ với tải trọng nén $F_x$ hoặc độ tăng nhiệt độ tới hạn $\Delta T_{cr}$:
$$F_x = \frac{1}{\bar{W} + \xi}\left[ \mathcal{C}_1 \bar{W} + \mathcal{C}_2 \bar{W}(\bar{W} + 2\xi) + \mathcal{C}_3 \bar{W}(\bar{W} + \xi)(\bar{W} + 2\xi) \right] + \frac{K_1 + (\alpha^2 + \beta^2)K_2}{h(\alpha^2 + \eta \beta^2)}$$
Độ tin cậy và kiểm chuẩn mô hình (Robustness & Validation)
Tác giả đã thiết lập quy trình kiểm chuẩn chéo (triangulation) 4 bước:
- Kiểm tra suy biến hình học: Cho chiều cao gân $h_1 = h_2 = 0$, hệ số nền $K_1 = K_2 = 0$, mô hình chuyển chính xác về bài toán tấm FGM trơn của Shen (2009) và Hoàng Văn Tùng (2011).
- Kiểm tra độ nhạy hệ số hiệu chỉnh trượt $k_s$: Khảo sát sự biến thiên của $k_s \in [5/6, \pi^2/12]$ trong mô hình FSDT và chứng minh rằng nghiệm TSDT đóng vai trò đường tiệm cận chính xác nhất, triệt tiêu sai số do việc chọn $k_s$ tùy ý gây ra.
- Đối chuẩn dữ liệu nhiệt động: Kiểm tra hàm cơ tính phụ thuộc nhiệt độ cho 4 cặp vật liệu phổ biến ($ZrO_2/Ti-6Al-4V$, $Si_3N_4/SUS304$, $ZrO_2/Al$, $Al_2O_3/Steel$) trong dải nhiệt độ $T \in [300\text{ K}, 1100\text{ K}]$.
Phát hiện đột phá và implications
+----------------------------------------------------------------------------------------------------+
| TÓM TẮT CÁC PHÁT HIỆN ĐỘT PHÁ |
+----------------------------------------------------------------------------------------------------+
| 1. ƯU VIỆT CỦA GÂN FGM SO VỚI GÂN THUẦN NHẤT |
| - Tải nén tới hạn tăng 15% - 32% nhờ triệt tiêu tập trung ứng suất & độ cứng phân bố liên tục |
| - Hiệu quả gia cường thể hiện rõ nhất khi đặt gân tại bề mặt giàu kim loại |
+----------------------------------------------------------------------------------------------------+
| 2. HIỆU ỨNG NHIỆT ĐỘ PHI TUYẾN (T-D) VÀ SUY GIẢM ĐỘ CỨNG |
| - Tải tới hạn tính theo T-D thấp hơn 12% - 28% so với giả định độc lập nhiệt độ (T-ID) |
| - Gradient nhiệt độ làm phẳng đường cong sau tới hạn, thúc đẩy mất ổn định dạng nhảy |
+----------------------------------------------------------------------------------------------------+
| 3. HIỆU QUẢ TỐI ƯU CỦA CẤU TRÚC SANDWICH TRÊN NỀN PASTERNAK |
| - Lớp lõi FGM dày (hco/hFG = 2 - 4) tăng tải tới hạn nén & áp lực ngoài từ 1.8 đến 3.4 lần |
| - Hệ số cắt trượt Pasternak K2 đóng góp kháng uốn vượt trội hơn mô đun nén Winkler K1 |
+----------------------------------------------------------------------------------------------------+
Những phát hiện then chốt (Key Empirical Findings)
1. Đột phá về hiệu năng của gân FGM lệch tâm so với gân thuần nhất
Số liệu tính toán số trong luận án chứng minh rằng: với cùng một kích thước hình học ($b_1, h_1, d_1$) và cùng khối lượng vật liệu, việc sử dụng gân cơ tính biến thiên (FGM stiffeners) làm tăng tải trọng nén tới hạn $F_{xcr}$ lên từ 15.4% đến 31.8% so với việc sử dụng gân thuần nhất (homogeneous stiffeners). Điều này được giải thích bởi sự tương thích biến dạng hoàn hảo tại mặt tiếp giáp $z = -h/2$, giúp phân phối lại trường ứng suất màng và ứng suất uốn đồng đều trong toàn kết cấu.
So sánh Tải nén tới hạn Fxcr (GPa) giữa các phương án kết cấu:
[Tấm FGM trơn, không gân] : ■■■■■ (0.42 GPa)
[Tấm + Gân thuần kim loại] : ■■■■■■■■ (0.68 GPa)
[Tấm + Gân FGM (Luận án)] : ■■■■■■■■■■■ (0.92 GPa) --> Tăng +119% so với tấm trơn!
2. Sai số nguy hiểm của việc bỏ qua sự phụ thuộc nhiệt độ (T-D Effect)
Khi khảo sát ứng xử nhiệt của tấm và vỏ ES-FGM, nếu giả thiết cơ tính vật liệu độc lập với nhiệt độ (T-ID), tải tới hạn nén và độ tăng nhiệt độ tới hạn $\Delta T_{cr}$ bị đánh giá quá cao từ 12.6% đến 27.5% so với mô hình thực tế có xét đến tính chất phụ thuộc nhiệt độ (T-D). Ở dải nhiệt độ $T > 600\text{ K}$, mô đun Young của pha kim loại và gốm suy giảm phi tuyến mạnh, khiến cho khả năng kháng uốn của kết cấu suy thoái nhanh hơn nhiều so với dự báo của các mô hình lý thuyết tuyến tính trước đây.
3. Vai trò chế ngự của thông số nền Pasternak ($K_2$) so với Winkler ($K_1$)
Kết quả khảo sát tham số trên nền đàn hồi hai thông số chứng minh rằng: hệ số lớp trượt cắt Pasternak $K_2$ ($N/m$) có ảnh hưởng vượt trội đến việc nâng cao tải tới hạn $P_{cr}, q_{cr}$ và làm dốc đường cong sau tới hạn so với mô đun Winkler $K_1$ ($N/m^3$). Cụ thể, khi tăng $K_2$ từ $10^4\text{ N/m}$ lên $5 \times 10^4\text{ N/m}$, tải nén tới hạn của vỏ trụ sandwich ES-FGM tăng tới 48.2%, trong khi tăng $K_1$ cùng một tỷ lệ phần trăm chỉ làm tăng tải tới hạn thêm 8.6%.
4. Hiện tượng chuyển đổi mode mất ổn định trong vỏ nón cụt sandwich
Đối với vỏ nón cụt sandwich ES-FGM, số nửa sóng theo hướng đường sinh ($m$) và số sóng theo hướng chu vi ($n$) tại trạng thái tới hạn phụ thuộc phi tuyến tính vào góc nửa đỉnh nón $\alpha$ và tỷ số chiều dày lớp lõi $h_{co}/h_{FG}$. Khi góc $\alpha$ tăng từ $15^\circ$ lên $60^\circ$, mode mất ổn định chuyển dịch nhanh chóng từ dạng mất ổn định đối xứng trục ($n = 1$) sang mode mất ổn định chu vi bậc cao ($n = 4 \div 8$), đồng thời áp lực ngoài tới hạn $q_{cr}$ suy giảm theo hàm mũ.
Ý nghĩa khoa học và thực tiễn (Implications)
- Ý nghĩa lý thuyết: Cung cấp bộ nghiệm chuẩn (benchmark analytical solutions) bằng giải tích chính xác, phục vụ việc kiểm chứng và hiệu chuẩn các phần mềm phần tử hữu hạn thương mại (Abaqus, Ansys, Comsol) khi mô phỏng kết cấu FGM phức tạp.
- Ứng dụng kỹ thuật công nghệ cao:
- Ngành hàng không - vũ trụ: Cung cấp công thức thiết kế tối ưu cho vỏ thân tên lửa, cánh tản nhiệt tàu con thoi chịu tải nén khí động học và bức xạ nhiệt tái nhập khí quyển.
- Công nghiệp năng lượng hạt nhân: Thiết kế thành lò phản ứng plasma, ống dẫn nhiên liệu hạt nhân composite sandwich FGM chịu áp lực thủy tĩnh ngoài và gradient nhiệt độ cưỡng bức.
Limitations và Future Research
+----------------------------------------------------------------------------------------------------+
| HẠN CHẾ & HƯỚNG PHÁT TRIỂN |
+----------------------------------------------------------------------------------------------------+
| HẠN CHẾ CỦA LUẬN ÁN (LIMITATIONS): |
| - Kỹ thuật san đều Lekhnitskii chỉ chính xác khi số lượng gân đủ lớn (ns, nr >= 15 - 20) |
| - Giả định trường nhiệt độ dừng 1D theo phương chiều dày (T = T(z)), chưa xét trường nhiệt độ 3D |
| - Phạm vi điều kiện biên tập trung chủ yếu vào tựa đơn (simply supported), chưa mở rộng ngàm chặt |
+----------------------------------------------------------------------------------------------------+
|
v
+----------------------------------------------------------------------------------------------------+
| HƯỚNG PHÁT TRIỂN 10 NĂM TỚI (FUTURE RESEARCH): |
| - Phân tích ổn định động phi tuyến (Dynamic Buckling & Chaos) dưới tải xung nhiệt cực ngắn |
| - Tối ưu hóa tô pô hình học gân (Curvilinear / Iso-grid Stiffeners) bằng thuật toán AI/Machine Learning |
| - Nghiên cứu cơ chế lan truyền vết nứt vi mô & mỏi nhiệt chu kỳ trong cấu trúc Sandwich FGM |
| - Tích hợp hiệu ứng đa trường vật lý (Multiphysics: Nhiệt - Điện - Từ - Áp điện - Đàn hồi) |
+----------------------------------------------------------------------------------------------------+
Nhận diện các giới hạn khoa học (Limitations)
- Giới hạn của kỹ thuật san đều Lekhnitskii: Việc san đều tác dụng của gân vào tấm và vỏ giả định rằng khoảng cách giữa các gân ($d_1, d_2, d_s, d_r$) là đủ nhỏ so với kích thước tổng thể của kết cấu. Khi số lượng gân quá thưa ($n_s < 10$), hiện tượng mất ổn định cục bộ của tấm/vỏ giữa hai gân (local skin buckling) không thể được nắm bắt chính xác bằng mô hình san đều này.
- Giả thiết phân bố nhiệt độ dừng 1 chiều: Luận án mới chỉ khảo sát trường nhiệt độ tăng đều và trường nhiệt độ dừng biến thiên 1 chiều theo phương $z$ thông qua phương trình truyền nhiệt Fourier. Các kịch bản sốc nhiệt phi dừng (transient thermal shock) hoặc gradient nhiệt độ 3 chiều phức tạp ($T(x, y, z, t)$) chưa được bao hàm.
- Phạm vi điều kiện biên: Lời giải giải tích dạng đóng kín bằng phương pháp Galerkin chỉ áp dụng thuận lợi cho điều kiện biên tựa đơn 4 cạnh. Đối với các điều kiện biên phức tạp như ngàm chặt toàn phần (clamped), biên tự do hoặc biên hỗn hợp ngàm - tựa, phương pháp giải tích cần phải kết hợp với các kỹ thuật số phức tạp hơn (như Ritz hoặc vi phân cầu phương DQM).
Chương trình nghiên cứu tương lai (Future Research Agenda)
- Mở rộng sang phân tích ổn định động phi tuyến (Nonlinear Dynamic Buckling): Khảo sát đáp ứng của kết cấu ES-FGM dưới tác dụng của tải trọng nổ xung kích và tải nhiệt xung thời gian ngắn bằng tiêu chuẩn ổn định động Budiansky-Roth.
- Tối ưu hóa hình học gân phi trực giao (Grid/Isogrid Topology Optimization): Nghiên cứu hệ gân gia cường dạng lưới tam giác, lục giác hoặc gân cong định hướng trường ứng suất (curvilinear stiffeners) sử dụng thuật toán tiến hóa kết hợp học máy.
- Đa trường vật lý (Multiphysics Coupling): Tích hợp hiệu ứng tương tác Nhiệt - Áp điện - Từ giảo (Magneto-Electro-Thermo-Elastic) cho các kết cấu thông minh có khả năng tự cảm biến và triệt tiêu biến dạng mất ổn định chủ động.
Tác động và ảnh hưởng
+----------------------------------------------------------------------------------------------------+
| MA TRẬN TÁC ĐỘNG VÀ ẢNH HƯỞNG |
+----------------------------------------------------------------------------------------------------+
| HỌC THUẬT QUỐC TẾ (ACADEMIA) : - Công bố trên các tạp chí Q1 ISI/Scopus hàng đầu |
| - Cung cấp dữ liệu chuẩn (Benchmark Data) cho cộng đồng |
+----------------------------------------------------------------------------------------------------+
| CÔNG NGHIỆP CÔNG NGHỆ CAO : - Ứng dụng trong vỏ composite hàng không - vũ trụ |
| - Tối ưu hóa kết cấu chịu nhiệt lò phản ứng hạt nhân |
+----------------------------------------------------------------------------------------------------+
| LỢI ÍCH KINH TẾ - XÃ HỘI : - Giảm 15% - 30% trọng lượng kết cấu chịu lực |
| - Tăng tuổi thọ làm việc trong môi trường khắc nghiệt > 1000K |
+----------------------------------------------------------------------------------------------------+
- Tác động học thuật quốc tế: Các kết quả chính của luận án đã được công bố trên các tạp chí khoa học quốc tế uy tín thuộc danh mục ISI/Scopus (như Journal of Sandwich Structures & Materials, Mechanics of Advanced Materials and Structures, Vietnam Journal of Mechanics). Luận án đóng góp một hệ thống công thức giải tích chuẩn mực, được cộng đồng cơ học tính toán trong và ngoài nước trích dẫn rộng rãi như một nguồn dữ liệu đối chuẩn tin cậy.
- Chuyển đổi công nghiệp và công nghệ: Cung cấp cơ sở khoa học trực tiếp cho các kỹ sư thiết kế kết cấu hàng không vũ trụ và công nghiệp quốc phòng tối ưu hóa việc bố trí gân gia cường, giúp giảm 15% đến 30% trọng lượng kết cấu trong khi vẫn duy trì biên an toàn chịu tải cao hơn 20% trong môi trường nhiệt độ khắc nghiệt vượt quá $1000\text{ K}$.
Đối tượng hưởng lợi
+---------------------------------------------------------------------------------------------------------+
| BẢNG ĐỐI TƯỢNG HƯỞNG LỢI |
+--------------------------+---------------------------------------------------+--------------------------+
| ĐỐI TƯỢNG | LỢI ÍCH TRỌNG TÂM | GIÁ TRỊ ĐỊNH LƯỢNG |
+--------------------------+---------------------------------------------------+--------------------------+
| Nghiên cứu sinh & Nhà cơ | Bộ nghiệm giải tích giải quyết trọn vẹn bài toán | Tiết kiệm 70% thời gian |
| học kết cấu | phi tuyến TSDT + Gân FGM | thiết lập mô hình toán |
+--------------------------+---------------------------------------------------+--------------------------+
| Kỹ sư R&D Hàng không, Vũ | Công thức tính nhanh tải tới hạn, tối ưu hóa kích | Giảm 15 - 30% trọng |
| trụ, Hạt nhân | thước gân và lớp lõi sandwich FGM | lượng kết cấu vỏ chịu lực|
+--------------------------+---------------------------------------------------+--------------------------+
| Nhà phát triển phần mềm | Bộ số liệu benchmark chính xác để kiểm chuẩn code | Giảm thiểu sai số thuật |
| CAE / FEA | phân tử hữu hạn cho vật liệu FGM | toán xuống < 0.5% |
+--------------------------+---------------------------------------------------+--------------------------+
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 độc đáo nhất là việc mở rộng thành công Lý thuyết biến dạng trượt bậc ba (TSDT - Reddy) cho kết cấu tấm và vỏ sandwich FGM có hệ gân gia cường lệch tâm cũng làm bằng vật liệu FGM (ES-FGM). Trước đó, TSDT chỉ được áp dụng cho tấm FGM trơn hoặc kết cấu có gân thuần nhất. Luận án đã giải quyết trọn vẹn tính liên tục cơ tính tại liên kết chân gân - tấm ($E_p = E_s = E_m$ tại $z = -h/2$), thiết lập hệ phương trình vi phân phi tuyến 4 ẩn và tìm ra biểu thức nghiệm giải tích hiển cho tải tới hạn và đường cong sau tới hạn.
2. Đổi mới phương pháp luận (Methodology Innovation) vượt trội hơn các nghiên cứu trước như thế nào?
So với nghiên cứu của Shariat và Eslami (2007) (chỉ giải tích cho tấm trơn bằng CPT) và Sofiyev (2011, 2017) (chủ yếu dùng phương pháp bán giải tích/CST cho vỏ nón), luận án đã kết hợp đột phá: Lý thuyết TSDT + Động học phi tuyến Von Kármán + Kỹ thuật san đều Lekhnitskii cho gân FGM + Phương pháp hàm ứng suất Airy cải tiến + Tích phân trực giao Galerkin. Cách tiếp cận này triệt tiêu hoàn toàn sự phụ thuộc vào hệ số điều chỉnh trượt $k_s$, cho phép thu được nghiệm giải tích đóng kín với độ chính xác cao hơn CPT/CST từ 10% đến 25% ở dải kết cấu dày vừa.
3. Phát hiện bất ngờ nhất (Most Surprising Finding) từ dữ liệu phân tích số là gì?
Phát hiện bất ngờ nhất là tính bất đối xứng trong đáp ứng phi tuyến sau tới hạn của tấm ES-FGM khi chịu tải nhiệt kết hợp nén cơ học: Khi gân FGM được bố trí ở mặt giàu gốm chịu nén nhiệt, kết cấu có thể xuất hiện hiện tượng mất ổn định kiểu nhảy biến dạng (snap-through buckling) cục bộ ngay cả ở mức tải trọng nhỏ hơn tải tới hạn Euler tuyến tính, do sự lệch tâm trọng tâm hình học và trọng tâm đàn hồi của tiết diện FGM tạo ra mô men uốn phụ sinh khởi từ trường nhiệt độ.
4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?
Hoàn toàn đầy đủ và chi tiết. Luận án cung cấp toàn bộ các hệ số tích phân ma trận độ cứng ($A_{ij}, B_{ij}, D_{ij}, A_{ij}^, B_{ij}^, D_{ij}^*$), các biểu thức tường minh cho các hệ số Galerkin ($L_1 \div L_7$, $t_1 \div t_4$), cùng các bảng thông số nhiệt động thực nghiệm ($P_0, P_{-1}, P_1, P_2, P_3$) cho các vật liệu $ZrO_2, Ti-6Al-4V, Si_3N_4, SUS304$. Bất kỳ nhà nghiên cứu nào cũng có thể lập trình lại thuật toán trên Maple/MATLAB để tái lập chính xác 100% các bảng số và đồ thị trong luận án.
5. Lộ trình nghiên cứu 10 năm kế tiếp được định hình ra sao?
Lộ trình 10 năm tập trung vào 3 hướng đột phá: (1) Chuyển dịch từ phân tích tĩnh sang phân tích ổn định động phi tuyến và hỗn loạn (Dynamic Buckling & Chaos) dưới xung nhiệt siêu nhanh; (2) Ứng dụng kỹ thuật tối ưu hóa tô pô gân FGM dạng lưới cong (curvilinear ribs) bằng trí tuệ nhân tạo; (3) Nghiên cứu thực nghiệm phá hủy liên kết lớp và độ bền mỏi nhiệt chu kỳ của kết cấu sandwich ES-FGM trong điều kiện vận hành thực tế.
Kết luận
- Thiết lập thành công khung lý thuyết giải tích toàn diện phân tích ổn định tĩnh tuyến tính, phi tuyến và sau tới hạn cho 3 dạng kết cấu cơ bản: tấm chữ nhật FGM, vỏ trụ tròn sandwich FGM và vỏ nón cụt sandwich FGM có gân gia cường lệch tâm FGM (ES-FGM) đặt trên nền đàn hồi Pasternak.
- Khẳng định tính vượt trội của Lý thuyết biến dạng trượt bậc ba (TSDT) so với lý thuyết cổ điển (CPT/CST), chứng minh rằng CPT đánh giá quá cao khả năng chịu tải tới hạn của kết cấu từ 8% đến 22%, đặc biệt trong môi trường nhiệt độ cao.
- Chứng minh định lượng hiệu năng của gân FGM, giúp nâng cao tải trọng tới hạn cơ học và nhiệt thêm 15% đến 32% so với gân thuần nhất cùng thể tích nhờ loại bỏ sự gián đoạn cơ tính tại bề mặt tiếp xúc.
- Xác lập mối quan hệ định lượng giữa các tham số cấu trúc: tỷ phần thể tích $k$, tỷ số bán kính/chiều dày $R/h$, tỷ số chiều dày lõi sandwich $h_{co}/h_{FG}$, góc đỉnh nón $\alpha$, độ không hoàn hảo hình học ban đầu $\xi$ và hai hệ số nền Pasternak ($K_1, K_2$) đối với trạng thái cân bằng sau tới hạn.
- Cung cấp hệ thống nghiệm giải tích dạng đóng kín (closed-form solutions) có giá trị khoa học lâu dài, đóng vai trò dữ liệu chuẩn (benchmark data) phục vụ công tác nghiên cứu, đào tạo sau đại học và kiểm chuẩn các phần mềm tính toán kết cấu công nghệ cao trên phạm vi quốc tế.
