Tổng quan về luận án

Nghiên cứu của tác giả Phạm Văn Hoàn (2024) với đề tài "Phân tích ổn định phi tuyến panel trụ và vỏ trụ làm bằng vật liệu FGP chịu tải cơ trong môi trường nhiệt" thuộc chuyên ngành Cơ kỹ thuật (mã số: 9520101), thực hiện tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn của PGS.TS. Lê Khả Hòa và PGS.TS. Đào Như Mai, là một công trình khoa học chuyên sâu và tiên phong trong lĩnh vực cơ học vật liệu composite tiên tiến. Trong bối cảnh ngành kỹ thuật hàng không vũ trụ, công trình ngầm, năng lượng nguyên tử và công nghệ quốc phòng đòi hỏi các kết cấu dạng vỏ mỏng vừa có trọng lượng siêu nhẹ, vừa sở hữu khả năng hấp thụ năng lượng và chịu được môi trường nhiệt cơ khắc nghiệt, vật liệu xốp có cơ tính biến đổi (Functionally Graded Porous material - FGP) đã trở thành đối tượng nghiên cứu then chốt của cơ học hiện đại.

Khoảng trống học thuật (Research Gap) trước nghiên cứu này rất rõ rệt: phần lớn các công trình công bố quốc tế trước đây chủ yếu tập trung vào vật liệu FGM truyền thống hoặc tấm phẳng FGP (Shen, 2009; Reddy & Starnes, 1993). Các khảo sát về ứng xử mất ổn định phi tuyến tĩnh và động của vỏ trụ và panel trụ FGP ba lớp (sandwich) có kể đến đầy đủ các yếu tố thực tế như: độ không hoàn hảo hình học ban đầu ($w^*$), hệ gân gia cường đa hướng (gân dọc, gân vòng, gân trực giao, gân xiên xoắn ốc), tương tác phức hợp với nền đàn hồi hai tham số Pasternak ($K_1, K_2$), và chịu tải trọng phức hợp (nén dọc trục $r_0$, áp lực ngoài $q$, tải xoắn thuần túy $\tau$, cùng trường nhiệt độ biến thiên $\Delta T$) hầu như chưa có lời giải giải tích tường minh khép kín.

Luận án đã giải quyết hệ thống câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  • Câu hỏi nghiên cứu 1 (RQ1): Quy luật phân bố lỗ xốp theo chiều dày (đối xứng, không đối xứng loại 1, không đối xứng loại 2, đồng đều) chi phối như thế nào đến trường ứng suất màng và ứng suất uốn của kết cấu vỏ?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để thiết lập được biểu thức giải tích dạng hiển (explicit closed-form formulations) xác định tải trọng tới hạn trên ($r_{0,\text{upper}}, \tau_{\text{upper}}, q_{\text{cr}}$) và tải trọng tới hạn dưới ($r_{0,\text{lower}}, \tau_{\text{lower}}$) cũng như đường cong cân bằng sau mất ổn định (post-buckling path)?
  • Câu hỏi nghiên cứu 3 (RQ3): Tác động gia cường của gân xoắn/gân xiên và phản lực nền Pasternak cải thiện bao nhiêu phần trăm khả năng chịu tải xoắn và nén khi có mặt độ không hoàn hảo hình học?
  • Giả thuyết khoa học 1 (H1): Phân bố xốp đối xứng với mật độ xốp cực đại ở mặt trung hòa sẽ cho khả năng chịu tải nén và uốn vượt trội so với phân bố xốp không đối xứng hoặc đồng đều do duy trì được độ cứng uốn ở hai lớp biên.
  • Giả thuyết khoa học 2 (H2): Bố trí gân xiên (oblique stiffeners) kết hợp với lớp phủ mặt FGM sẽ triệt tiêu hiệu quả biến dạng trượt xoắn và dịch chuyển điểm rẽ nhánh mất ổn định sang miền tải trọng cao hơn.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp giữa lý thuyết vỏ thoải Donnell, giả thuyết phi tuyến hình học von-Kármán, lý thuyết biến dạng cắt bậc nhất (FSDT), phương pháp san đều tác dụng của gân Leckhnitsky và kỹ thuật chiếu trực giao Galerkin. Đóng góp đột phá của công trình được định lượng cụ thể: thiết lập hệ phương trình vi phân đạo hàm riêng phi tuyến bậc cao theo hàm ứng suất $f(x,y)$ và độ võng $w(x,y)$, xây dựng thành công thuật toán tính toán trên MATLAB cho phép dự báo chính xác các điểm rẽ nhánh và ứng xử snap-through với sai số dưới $1.5%$ so với các công bố kinh điển quốc tế. Phạm vi nghiên cứu bao trùm toàn diện các cấu hình hình học từ tấm phẳng đến panel trụ và vỏ trụ kín, mở ra cơ sở lý thuyết vững chắc cho thiết kế kết cấu vỏ mỏng chịu nhiệt thế hệ mới.

Literature Review và Positioning

Lịch sử phát triển của cơ học kết cấu composite biến tính bắt đầu từ khái niệm FGM của các nhà khoa học vật liệu Nhật Bản năm 1984 (Koizumi, 1984). Tiếp nối dòng nghiên cứu này, Shen và Liew (2004, 2009) đã áp dụng lý thuyết biến dạng cắt bậc cao (HSDT) và phương pháp nhiễu để khảo sát ổn định của panel trụ FGM nhiều lớp chịu tác động cơ - nhiệt. Tại Việt Nam, trường phái nghiên cứu của GS. Nguyễn Đình Đức và cộng sự (2011, 2017) cùng nhóm nghiên cứu của PGS. Đào Văn Dũng, PGS. Lê Khả Hòa (2012, 2015, 2017) đã phát triển mạnh mẽ phương pháp giải tích cho vỏ FGM có gân gia cường trên nền đàn hồi.

Tuy nhiên, sự chuyển dịch sang vật liệu xốp có cơ tính biến đổi (FGP) - dạng vật liệu có lỗ rỗng phân bố liên tục theo hàm toán học (Chen et al., 2016; Kitipornchai et al., 2017) - đã đặt ra những thách thức cơ học hoàn toàn mới. Trong y văn tồn tại hai luồng quan điểm và tranh luận học thuật sâu sắc:

  • Tranh luận 1 (Về tính hợp lệ của mô hình giải tích vỏ mỏng): Một số tác giả (như Zghal et al., 2018; Kweon & Hong, 1994) cho rằng khi xuất hiện các lỗ xốp phân bố ngẫu nhiên hoặc bất đối xứng, biến dạng cắt ngang trở nên thống trị và bắt buộc phải dùng phương pháp phần tử hữu hạn phi tuyến 3D (Nonlinear 3D-FEM) để tránh đánh giá quá cao độ cứng kết cấu. Ngược lại, trường phái giải tích (Shen, 2009; Đào Văn Dũng & Lê Khả Hòa, 2015) chứng minh rằng lý thuyết vỏ Donnell kết hợp với quan hệ biến dạng von-Kármán hoàn toàn kiểm soát được độ chính xác đối với vỏ mỏng và vừa, đồng thời cung cấp bản chất vật lý tường minh mà phương pháp số thuần túy bị giới hạn bởi thời gian tính toán và độ hội tụ.
  • Tranh luận 2 (Về cơ chế mất ổn định): Quan điểm truyền thống xem mất ổn định thuần túy là hiện tượng rẽ nhánh Euler (bifurcation buckling) xảy ra ở trạng thái màng không có biến dạng uốn trước tới hạn. Tuy nhiên, các nghiên cứu thực nghiệm và lý thuyết của Belica & Magnucki (2020) chỉ ra rằng với vật liệu xốp FGP và kết cấu có độ không hoàn hảo hình học ban đầu $w^*$, mất ổn định xảy ra theo kiểu cực trị (limit-point buckling / snap-through) với sự xuất hiện đồng thời của tải tới hạn trên ($r_{0,\text{upper}}$) và tải tới hạn dưới ($r_{0,\text{lower}}$).

Vị trí học thuật của luận án Phạm Văn Hoàn được xác lập chính xác tại giao điểm của các khoảng trống nêu trên. Bằng việc định vị mô hình kết cấu sandwich FGP 3 lớp (lớp lõi xốp đối xứng ở giữa hai lớp phủ FGM thuần nhất) và hệ gân FGM gia cường 4 dạng (dọc, vòng, trực giao, xiên), luận án đã vượt lên trên các giới hạn của nghiên cứu trước.

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

  1. So với nghiên cứu của Jaunky & Knight (1999) trên tạp chí Composite Structures sử dụng phương pháp Rayleigh–Ritz cho panel FGM nhiều lớp, mô hình giải tích của luận án cho kết quả tải tới hạn nén dọc trục tương đồng chặt chẽ (độ lệch tuyệt đối $< 1.2%$), đồng thời mở rộng được sang trường hợp vật liệu xốp FGP có hệ số rỗng $e_0$ thay đổi và tích hợp phản lực nền Pasternak.
  2. So với công trình của Shen (2009) trên tạp chí Thin-Walled Structures khảo sát mất ổn định nhiệt của panel FGM, giải pháp giải tích của luận án đã bao quát thêm ứng xử sau mất ổn định (post-buckling response) của vỏ sandwich FGP chịu đồng thời tải xoắn $\tau$ và áp lực ngoài $q$, đưa ra công thức tường minh cho điểm cực trị mà mô hình của Shen chưa giải quyết trọn vẹn cho cấu trúc xốp.

Đó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à hoàn thiện lý thuyết mất ổn định phi tuyến của hệ thanh vỏ mỏng (Thin-walled shell buckling theory) bằng việc tích hợp các quy luật cơ tính biến đổi của vật liệu xốp FGP vào hệ phương trình vi phân liên hợp. Cụ thể:

  • Mở rộng định luật Hooke tổng quát cho môi trường xốp liên tục không đẳng hướng với mô đun đàn hồi Young $E(z)$ và hệ số Poisson $\nu(z)$ biến thiên phi tuyến theo chiều dày $z \in [-h/2, h/2]$: $$\text{Loại a (Đối xứng): } E(z) = E_m \left[ 1 - e_0 \cos\left(\frac{\pi z}{h}\right) \right]$$ $$\text{Loại b, c (Bất đối xứng): } E(z) = E_m \left[ 1 - e_0 \cos\left(\frac{\pi z}{2h} \pm \frac{\pi}{4}\right) \right]$$ $$\text{Loại d (Đồng đều): } E(z) = E_m (1 - e_0 \alpha)$$
  • Xây dựng mô hình toán học giải thích sự dịch chuyển trục trung hòa hình học và trục trung hòa đàn hồi trong kết cấu panel sandwich FGP ba lớp ($h = h_{\text{core}} + h_{\text{FG}}$), chứng minh rằng khi $E_2 = \int_{-h/2}^{h/2} z E_{\text{sh}}(z) dz = 0$ (đối với phân bố xốp đối xứng), hiện tượng ghép nối uốn - kéo (bending-extension coupling) bị triệt tiêu, giúp tối ưu hóa khả năng ổn định đàn hồi của kết cấu.
  • Xác lập hệ mệnh đề lý thuyết (Theoretical Propositions):
    • Mệnh đề 1 (P1): Tải trọng tới hạn trên $r_{0,\text{upper}}$ đạt cực đại toàn cục khi hàm phân bố độ rỗng tập trung giá trị lớn nhất tại mặt $z=0$ và giảm dần về 0 tại các biên $z = \pm h/2$.
    • Mệnh đề 2 (P2): Đường cong cân bằng tải - độ võng sau tới hạn của vỏ không hoàn hảo ($w^* \neq 0$) chuyển dịch tiệm cận trơn về nhánh cân bằng của vỏ hoàn hảo khi biên độ chuyển vị $W/h \gg 1$.
+-----------------------------------------------------------------------------+
|                             KHUNG PHÂN TÍCH ĐỘC ĐÁO                         |
+-----------------------------------------------------------------------------+
|                                                                             |
|   +--------------------------+           +------------------------------+   |
|   | 1. CƠ HỌC VẬT LIỆU FGP   |           | 2. HÌNH HỌC & ĐỘ HOÀN HẢO    |   |
|   | - Hàm phân bố rỗng e0    |           | - Độ cong vỏ Donnell         |   |
|   | - Cấu trúc Sandwich 3 lớp|           | - Phi tuyến von-Kármán       |   |
|   | - Quy luật lũy thừa FGM  |           | - Độ lệch ban đầu w*(x,y)    |   |
|   +------------+-------------+           +--------------+---------------+   |
|                |                                        |                   |
|                +-------------------+--------------------+                   |
|                                    |                                        |
|                                    v                                        |
|                     +-------------------------------+                       |
|                     | 3. MÔ HÌNH HỆ THỐNG PHỨC HỢP  |                       |
|                     | - San đều gân Leckhnitsky     |                       |
|                     | - Nền Pasternak (K1, K2)      |                       |
|                     | - Trường nhiệt đồng đều DeltaT|                       |
|                     +---------------+---------------+                       |
|                                     |                                       |
|                                     v                                       |
|               +-------------------------------------------+                 |
|               | 4. THIẾT LẬP HỆ ĐẠO HÀM RIÊNG PHI TUYẾN   |                 |
|               |    A3*Laplace^4(f) + Laplace_R(w) = ...   |                 |
|               |    A4*Laplace^4(w) - Laplace_R(f) = ...   |                 |
|               +---------------------+---------------------+                 |
|                                     |                                       |
|                                     v                                       |
|               +-------------------------------------------+                 |
|               | 5. TRỰC GIAO HÓA GALERKIN & THUẬT TOÁN    |                 |
|               |    - Biểu thức hiển: r0_upper, r0_lower   |                 |
|               |    - Quỹ đạo tải - độ võng (MATLAB)       |                 |
|               +-------------------------------------------+                 |
+-----------------------------------------------------------------------------+

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

Khung phân tích của luận án kết hợp chặt chẽ 4 lý thuyết nền tảng:

  1. Lý thuyết vỏ mỏng Donnell tích hợp biến dạng von-Kármán: Mô tả chính xác mối quan hệ phi tuyến giữa trường chuyển vị $(u, v, w)$ và trường biến dạng $(\varepsilon_x^0, \varepsilon_y^0, \gamma_{xy}^0, k_x, k_y, k_{xy})$.
  2. Kỹ thuật san đều gân Leckhnitsky: Quy đổi độ cứng rời rạc của các gân gia cường FGM (kích thước $h_s, b_s, d_s$ cho gân dọc; $h_r, b_r, d_r$ cho gân vòng; $h_p, b_p, d_p, \theta$ cho gân xiên) vào các ma trận độ cứng màng và độ cứng uốn tương đương của vỏ trụ.
  3. Mô hình tương tác nền Pasternak 2 tham số: Phản ánh chính xác cả phản lực đàn hồi chống lún ($K_1$) và tương tác cắt giữa các phần tử nền ($K_2$), khắc phục hoàn toàn nhược điểm phi vật lý của mô hình Winkler cổ điển.
  4. Phương pháp hàm ứng suất Airy kết hợp phép chiếu trực giao Galerkin: Chuyển đổi hệ hai phương trình vi phân đạo hàm riêng phi tuyến cấp 4 thành hệ phương trình đại số phi tuyến dạng: $$r_0 = \frac{S_1 W + S_2 W^2 + S_3 W^3 + \dots}{S_5 (W + \xi h) \dots}$$ Cho phép tính đạo hàm bậc nhất $dr_0/dW = 0$ để trích xuất trực tiếp giá trị độ võng tới hạn $W_{\text{th}}$ và tải trọng tới hạn dưới $r_{0,\text{lower}}$ mà không cần sử dụng các thuật toán lặp dò nghiệm Newton-Raphson phức tạp.

Điều kiện biên được thiết lập rõ ràng cho hai trường hợp điển hình: bốn cạnh tựa đơn ($SSSS$) và hai cạnh cong tựa đơn kết hợp hai cạnh thẳng ngàm ($SSCC$), xác định ranh giới ứng dụng chính xác cho các cấu kiện thực tế trong thân tàu và vỏ động cơ.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (Positivism) với lập luận diễn dịch toán học (Deductive mathematical modeling) chuẩn mực của cơ học giải tích. Thiết kế nghiên cứu kết hợp giữa mô hình hóa cơ học giải tích đa quy mô (Multi-scale modeling) và kiểm chứng số thông qua thuật toán máy tính. Cấu trúc thiết kế được phân chia làm 3 cấp độ:

  • Cấp độ vật liệu (Micro-meso level): Mô hình hóa cấu trúc xốp vi mô và thành phần gốm/kim loại ($ZrO_2/Ti-6Al-4V$, $Si_3N_4/SUS304$, $Al_2O_3/Al$) thông qua các hàm phân bố độ rỗng $e_0$ và tỉ phần thể tích $k$.
  • Cấp độ cấu kiện (Component level): Thiết lập phương trình cơ bản cho panel trụ và vỏ trụ kín sandwich ba lớp gia cường gân.
  • Cấp độ hệ thống chịu tải đa trường (Multiphysics system level): Đặt cấu kiện vào không gian tương tác cơ - nhiệt - nền đàn hồi với các tải trọng nén $r_0$, áp lực $q$, xoắn $\tau$ và gradient nhiệt $\Delta T$.

Quy trình nghiên cứu rigorous

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

  1. Thiết lập quan hệ hình học và ứng suất: Xác định các thành phần lực dãn $(N_x, N_y, N_{xy})$ và mô-men uốn $(M_x, M_y, M_{xy})$ thông qua việc lấy tích phân đại số giải tích theo chiều dày từ $-h/2$ đến $h/2$.
  2. Thiết lập phương trình tương thích biến dạng và phương trình cân bằng: Biểu diễn hệ phương trình Donnell thông qua toán tử vi phân bậc 4 ($\nabla^4$) và toán tử độ cong ($1/R$).
  3. Lựa chọn hàm dạng nghiệm thỏa mãn điều kiện biên:
    • Với biên tựa đơn ($SSSS$): $$w(x,y) = W \sin\left(\frac{m\pi x}{a}\right) \sin\left(\frac{n\pi y}{b}\right)$$ $$w^*(x,y) = \xi h \sin\left(\frac{m\pi x}{a}\right) \sin\left(\frac{n\pi y}{b}\right)$$
    • Với biên ngàm-tựa ($SSCC$): $$w(x,y) = W \sin\left(\frac{m\pi x}{a}\right) \left[ 1 - \cos\left(\frac{2n\pi y}{b}\right) \right]$$
  4. Trực giao hóa Galerkin: Nhân phương trình vi phân với hàm trọng số tương ứng và lấy tích phân trên toàn bộ miền diện tích $\Omega = [0, a] \times [0, b]$ để loại bỏ các biến không gian, chuyển về dạng đại số phi tuyến.
  5. Cực tiểu hóa số sóng (Mode minimization): Tải trọng tới hạn trên và dưới được xác định bằng cách quét và tìm giá trị cực tiểu của $r_0(m,n)$ trên miền số sóng không gian $(m, n \in \mathbb{N}^*)$.
+-----------------------------------------------------------------------------+
|                           QUY TRÌNH NGHIÊN CỨU RIGOROUS                     |
+-----------------------------------------------------------------------------+
|                                                                             |
|  [Định luật Hooke FGP]                                                      |
|           |                                                                 |
|           v                                                                 |
|  [Tích phân ứng suất qua chiều dày h]                                       |
|           |                                                                 |
|           v                                                                 |
|  [Ma trận độ cứng tương đương Dij (kể cả san đều gân Leckhnitsky)]          |
|           |                                                                 |
|           v                                                                 |
|  [Hệ 2 phương trình vi phân Donnell (w, f)]                                 |
|           |                                                                 |
|           v                                                                 |
|  [Chọn nghiệm Fourier phù hợp SSSS / SSCC]                                  |
|           |                                                                 |
|           v                                                                 |
|  [Tích phân Galerkin trên miền phẳng Omega: Integral (Res * phi) dxdy = 0]  |
|           |                                                                 |
|           v                                                                 |
|  [Biểu thức hiển r0(W, m, n, e0, k, K1, K2, DeltaT)]                        |
|           |                                                                 |
|           v                                                                 |
|  [Giải tích cực trị: dr0/dW = 0 --> r0_upper & r0_lower]                    |
|           |                                                                 |
|           v                                                                 |
|  [Tối ưu hóa tìm mode nguy hiểm nhất: min_{m,n} r0_upper]                   |
|           |                                                                 |
|           v                                                                 |
|  [Mô phỏng số trên MATLAB & Đối chuẩn kiểm chứng]                           |
|                                                                             |
+-----------------------------------------------------------------------------+

Data và phân tích

Toàn bộ hệ thống công thức giải tích được lập trình hóa trên phần mềm MATLAB. Thuật toán tự động hóa việc tính toán các ma trận độ cứng $D_{ij}$, xác định nghiệm của các đa thức bậc cao và vẽ đường cong cân bằng tĩnh. Độ tin cậy và tính hợp lệ (construct & external validity) của phương pháp được kiểm chứng tam giác (Triangulation) thông qua việc so sánh đối chiếu đa chiều:

  • Kiểm chứng 1 (Bảng 2.1): So sánh tải tới hạn trên $N_{\text{cr}} = r_{0,\text{upper}}/(bh)$ của panel trụ đẳng hướng với kết quả của Shen (2009) theo lý thuyết Donnell và Jaunky (1999) theo HSDT/FEM. Ở cấu hình $a/b = 1.0, a/R = 0.2, b/h = 100$, kết quả của luận án đạt $N_{\text{cr}} = 3.9683 \times 10^{-4}$, hoàn toàn trùng khớp với giá trị $3.9683 \times 10^{-4}$ của Shen và độ sai lệch $< 0.8%$ so với FEM của Jaunky.
  • Kiểm chứng 2 (Bảng 2.2): So sánh tải tới hạn nén không thứ nguyên $N_x/(ER)$ với công bố của Turvey (1995) và Shen (2009). Với $a/b = 0.5, a/R = 0.5, b/h = 200$, kết quả luận án cho $N_x/(ER) = 0.0553 \times 10^{-5}$ ứng với mode $(m=4, n=2)$, đồng nhất hoàn toàn với công bố quốc tế.
  • Kiểm định độ bền vững (Robustness checks): Khi cho bán kính $R \rightarrow \infty$, phương trình tự động thoái hóa về bài toán ổn định của tấm phẳng FGP; khi hệ số rỗng $e_0 \rightarrow 0$, kết quả trở về đúng nghiệm của panel FGM đặc; khi độ không hoàn hảo $\xi \rightarrow 0$, phương trình mô tả chính xác trạng thái rẽ nhánh cổ điển.

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

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

Phân tích số liệu tham số sâu rộng trong luận án đã đưa ra 5 phát hiện mang tính đột phá với bằng chứng định lượng chính xác:

  1. Ưu thế tuyệt đối của mô hình phân bố xốp đối xứng (Loại a): Trong 4 dạng phân bố rỗng được khảo sát, kết cấu có phân bố xốp đối xứng luôn đạt tải trọng tới hạn cao nhất, tiếp theo là hai dạng bất đối xứng (Loại b và Loại c có khả năng chịu lực tương đương nhau), và thấp nhất là phân bố xốp đồng đều (Loại d). Cụ thể, với panel trụ xốp FG có hệ số rỗng $e_0 = 0.5$, tải trọng tới hạn trên đạt:

    • Loại a (Đối xứng): $r_{0,\text{upper}} = 293.7449\text{ MN}$
    • Loại b, c (Không đối xứng): $r_{0,\text{upper}} = 267.2059\text{ MN}$ (thấp hơn $9.03%$)
    • Loại d (Đồng đều): $r_{0,\text{upper}} = 260.6729\text{ MN}$ (thấp hơn $11.26%$) Đường cong quan hệ tải - độ võng ($r_0 - W/h$) của Loại a luôn nằm ở vị trí cao nhất trên đồ thị cân bằng.

  2. Quy luật suy giảm độ cứng phi tuyến theo độ rỗng $e_0$: Hệ số độ rỗng $e_0$ tỷ lệ nghịch mạnh mẽ với khả năng chịu tải. Khi $e_0$ tăng từ $0.0$ (vật liệu đặc hoàn toàn) lên $0.5$, tải trọng tới hạn trên $r_{0,\text{upper}}$ của panel phân bố đối xứng giảm đột ngột từ $393.9360\text{ MN}$ xuống $212.9157\text{ MN}$ (suy giảm đến $45.95%$). Điều này chứng minh rằng việc tăng độ rỗng để giảm trọng lượng phải được đánh đổi bằng sự suy giảm đáng kể độ cứng chống mất ổn định.

  3. Ảnh hưởng phi tuyến mạnh của tham số hình học ($a/b, a/R$): Tải trọng nén tới hạn trên giảm mạnh khi tỷ số chiều dài/chiều rộng ($a/b$) tăng (panel càng dài càng yếu), nhưng lại tăng vọt khi tỷ số độ cong ($a/R$) tăng. Cụ thể, khi giữ nguyên $a/b = 2.0$, việc tăng tỷ số $a/R$ từ $0.2$ lên $0.8$ đã làm tăng tải tới hạn trên từ $r_{0,\text{upper}} = 132.84\text{ MN}$ lên $401.76\text{ MN}$ (tăng trưởng vượt bậc trên $302%$).

  4. Hiện tượng mất ổn định sớm do độ không hoàn hảo hình học ban đầu ($w^*$): Khác với kết cấu hoàn hảo chỉ xuất hiện độ võng khi tải trọng đạt ngưỡng rẽ nhánh ($r_0 = r_{0,\text{upper}}$), kết cấu có độ không hoàn hảo ($\xi = 0.1, 0.3, 0.5$) xuất hiện chuyển vị võng ngay từ khi bắt đầu gia tải ($r_0 > 0$), và đường cong cân bằng chuyển dịch trơn qua điểm cực trị mà không có hiện tượng nhảy bước rẽ nhánh gián đoạn.

  5. Hiệu ứng gia cường của gân xiên (oblique) và nền Pasternak trong môi trường nhiệt: Dưới tác dụng của mô-men xoắn $\tau$, việc bố trí gân xiên theo góc xoắn tối ưu $\theta$ làm tăng tải trọng xoắn tới hạn $\tau_{\text{upper}}$ lên gấp $2.5 - 3.2$ lần so với vỏ không gân. Đồng thời, nền Pasternak với độ cứng chống cắt $K_2$ đóng vai trò cốt lõi trong việc nâng cao đường cong hậu mất ổn định, bù đắp trực tiếp sự suy giảm khả năng chịu lực do ứng suất nén nhiệt sinh ra bởi độ tăng nhiệt độ $\Delta T$.

+-----------------------------------------------------------------------------+
|              BẢNG TỔNG HỢP SO SÁNH SỐ LIỆU ĐỘT PHÁ TỪ LUẬN ÁN               |
+-----------------------------------------------------------------------------+
| Thông số khảo sát       | Cấu hình / Giá trị          | Tải tới hạn r0_upper|
+-------------------------+-----------------------------+---------------------+
| Mô hình phân bố xốp     | Loại a (Đối xứng)           | 293.7449 MN         |
| (e0 = 0.5)              | Loại b, c (Không đối xứng)  | 267.2059 MN         |
|                         | Loại d (Đồng đều)           | 260.6729 MN         |
+-------------------------+-----------------------------+---------------------+
| Ảnh hưởng độ rỗng e0    | e0 = 0.0 (Kim loại đặc)     | 393.9360 MN (100%)  |
| (Loại a)                | e0 = 0.5 (Xốp 50%)          | 212.9157 MN (54.05%)|
+-------------------------+-----------------------------+---------------------+
| Tỷ số độ cong a/R       | a/R = 0.2 (Độ cong nhỏ)     | 132.84 MN           |
| (a/b = 2.0)             | a/R = 0.8 (Độ cong lớn)     | 401.76 MN (+302.4%) |
+-------------------------+-----------------------------+---------------------+
| Điều kiện biên          | SSSS (4 cạnh tựa đơn)       | Mức tham chiếu      |
|                         | SSCC (2 cong tựa, 2 ngàm)   | Tăng ~22 - 35%      |
+-------------------------+-----------------------------+---------------------+

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp hệ thống nghiệm giải tích giải quyết trọn vẹn bài toán cơ học phi tuyến vật liệu xốp cấu trúc sandwich, đóng góp một chương mới vào lý thuyết ổn định vỏ mỏng composite.
  • Về mặt phương pháp luận: Xây dựng công cụ giải tích - số kết hợp, cho phép phân tích tham số nhanh chóng, vượt qua rào cản tính toán chi phí cao của các phần tử hữu hạn phi tuyến thương mại (như ABAQUS, ANSYS).
  • Về mặt kỹ thuật thực tiễn: Định hướng cho các kỹ sư kết cấu lựa chọn công nghệ chế tạo bọt kim loại rỗng (bột kim loại, in 3D laser melting) sao cho đạt được phân bố xốp đối xứng với mật độ tập trung ở lõi, tối ưu hóa tỷ lệ $h_{\text{core}}/h_{\text{FG}}$ để đạt độ bền uốn tối đa với khối lượng kết cấu nhỏ nhất.
  • Về mặt chính sách và tiêu chuẩn: Cung cấp cơ sở dữ liệu tính toán tin cậy làm tiền đề cho việc xây dựng các tiêu chuẩn thiết kế vỏ mỏng vật liệu mới trong ngành công nghiệp quốc phòng và xây dựng công trình biển tại Việt Nam.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án vẫn tồn tại một số giới hạn học thuật cần được nhìn nhận khách quan:

  1. Giả thiết môi trường đàn hồi tuyến tính: Nghiên cứu giả định vật liệu làm việc hoàn toàn trong miền đàn hồi, bỏ qua hiện tượng dẻo hóa cục bộ xung quanh các vi lỗ xốp khi chịu ứng suất tập trung cao.
  2. Xấp xỉ lý thuyết vỏ Donnell: Lý thuyết Donnell chỉ cho độ chính xác cao đối với vỏ mỏng và độ cong vừa phải ($h/R \le 1/50$), chưa xét đến hiệu ứng biến dạng trượt ngang bậc cao xuyên suốt chiều dày đối với các kết cấu panel sandwich rất dày.
  3. Mô hình nhiệt độ đồng đều: Luận án mới khảo sát trường hợp tăng nhiệt độ đồng đều $\Delta T$, chưa giải phương trình truyền nhiệt phi tuyến 3D theo thời gian với gradient nhiệt phi tuyến xuyên tâm và dọc trục.
  4. Giả thiết liên kết lớp dính bám hoàn hảo: Bỏ qua hiện tượng tách lớp (delamination) giữa lớp phủ FGM và lớp lõi xốp FGP dưới tải trọng trượt cắt lớn.

Chương trình nghiên cứu tiếp theo (Future Research Agenda) cần tập trung vào 4 hướng chiến lược:

  • Hướng 1: Mở rộng nghiên cứu sang bài toán ổn định động lực học phi tuyến (Nonlinear dynamic buckling) và đáp ứng dao động tự do/cưỡng bức dưới tải trọng xung va đập (blast/impact loading).
  • Hướng 2: Tích hợp mô hình phá hủy đàn dẻo (Elasto-plastic damage mechanics) và cơ học đứt gãy để đánh giá sự phát triển vết nứt từ các vi lỗ rỗng trong môi trường nhiệt độ cao.
  • Hướng 3: Nghiên cứu bài toán tối ưu hóa đa mục tiêu (Multi-objective topology optimization) đồng thời tỷ trọng, khả năng tiêu âm cách nhiệt và độ cứng chống mất ổn định của vỏ trụ FGP.
  • Hướng 4: Mở rộng mô hình giải tích cho vỏ composite lai nano (hybrid nanocomposite) gia cường sợi graphene/carbon nanotube trong nền xốp FG.

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

Luận án của NCS. Phạm Văn Hoàn mang lại giá trị học thuật và ứng dụng sâu rộng:

  • Tác động học thuật (Academic Impact): Các bài báo trích xuất từ luận án được công bố trên các tạp chí khoa học uy tín trong nước và quốc tế, đóng góp tài liệu tham khảo chất lượng cao cho cộng đồng cơ học tính toán vật liệu composite. Dự báo công trình sẽ thu hút trích dẫn lớn từ các nghiên cứu sinh và chuyên gia cơ học kết cấu vỏ mỏng.
  • Chuyển đổi công nghiệp (Industrial Transformation): Trực tiếp phục vụ các ngành công nghiệp công nghệ cao:
    • Hàng không - Vũ trụ: Thiết kế thân tên lửa đẩy, vỏ khoang tàu vũ trụ siêu nhẹ chịu tải nén và ma sát nhiệt khi tái nhập khí quyển.
    • Công nghiệp tàu thủy và công trình ngầm: Chế tạo vỏ chịu áp lực cao, panel tiêu âm giảm chấn cho tàu ngầm và kết cấu giàn khoan ngoài khơi.
    • Xây dựng dân dụng & công nghiệp: Chế tạo các tấm panel sandwich bọt nhôm cách âm, cách nhiệt, chống cháy cho các tòa nhà chọc trời và vách ngăn chống ồn đường cao tốc.
  • Lợi ích xã hội: Giúp giảm tiêu hao nhiên liệu trong giao thông vận tải nhờ tối ưu hóa giảm trọng lượng kết cấu, đồng thời nâng cao mức độ an toàn chống sụp đổ công trình dưới điều kiện cháy nổ hoặc va chạm cơ học.

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

  1. Nghiên cứu sinh và Giảng viên Cơ kỹ thuật: Tiếp cận phương pháp thiết lập lời giải giải tích tường minh cho hệ phương trình vi phân phi tuyến bậc cao; tận dụng mã nguồn MATLAB để phát triển tiếp các bài toán vỏ thoải, vỏ nón và vỏ cầu.
  2. Các nhà nghiên cứu vật liệu tiên tiến: Nắm bắt được quy luật ảnh hưởng của vi cấu trúc xốp đến ứng xử cơ học vĩ mô để định hướng quy trình công nghệ chế tạo bọt kim loại (như phương pháp luyện kim bột, thiêu kết sợi).
  3. Kỹ sư R&D trong ngành Quốc phòng và Hàng không: Sử dụng trực tiếp các công thức hiển ($r_{0,\text{upper}}, r_{0,\text{lower}}, \tau_{\text{upper}}, q_{\text{cr}}$) để tính toán sơ bộ kết cấu vỏ bọc khí tài quân sự mà không cần phụ thuộc vào phần mềm thương mại đắt tiền.
  4. Các cơ quan quản lý và thiết kế công trình: Áp dụng các kết quả nghiên cứu vào việc thẩm định độ ổn định của các bồn bể chứa áp lực, tháp xi-lô bằng vật liệu composite mới.

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ì? Đó là việc thiết lập thành công hệ phương trình vi phân phi tuyến liên hợp kín và nghiệm giải tích dạng hiển xác định chính xác tải trọng tới hạn trên và tải trọng tới hạn dưới cho panel và vỏ trụ sandwich FGP ba lớp có xét đến độ không hoàn hảo hình học và tương tác nền Pasternak. Luận án đã làm sáng tỏ bản chất cơ học của hiện tượng chuyển dịch trục trung hòa do gradient độ rỗng tạo ra.

2. Đột phá về phương pháp luận so với các nghiên cứu trước đây thể hiện ở điểm nào? So với Shen (2009) (dùng phương pháp nhiễu cho FGM đặc) và Jaunky & Knight (1999) (dùng số trị Rayleigh-Ritz), luận án đã kết hợp nhuần nhuyễn kỹ thuật san đều gân Leckhnitsky với phép chiếu trực giao Galerkin trên không gian nghiệm hàm lượng giác 2 chiều, cho phép giải trực tiếp bài toán phi tuyến lớn von-Kármán mà không làm mất đi các thành phần ghép nối phi tuyến quan trọng.

3. Phát hiện nào bất ngờ nhất từ kết quả phân tích số? Phát hiện bất ngờ nhất là mô hình phân bố xốp đồng đều (Loại d) - vốn là mô hình thường được giả định trong các tính toán kỹ thuật đơn giản vì tính thuận tiện - lại là mô hình cho khả năng chịu lực nén kém nhất trong tất cả các dạng phân bố rỗng ($r_{0,\text{upper}}$ thấp hơn tới $11.26%$ so với phân bố đối xứng Loại a tại $e_0 = 0.5$). Điều này cảnh báo các kỹ sư không được đơn giản hóa cấu trúc xốp trong thiết kế chịu lực an toàn.

4. Luận án có cung cấp đầy đủ giao thức để tái lập kết quả (Replication Protocol) không? Có. Toàn bộ các hệ số tích phân độ cứng $D_{ij}$, các hệ số vi phân $A_1, A_2, A_3, A_4$, các hàm trọng số Galerkin $S_1 \dots S_5$, cùng thông số cơ lý của vật liệu ($ZrO_2$, $Ti-6Al-4V$, thép, bọt nhôm) và thuật toán MATLAB đều được trình bày minh bạch, chi tiết từng bước trong các chương 2, 3, 4 và phần phụ lục, đảm bảo khả năng tái lập kết quả độc lập đạt $100%$.

5. Lộ trình phát triển nghiên cứu 10 năm từ luận án này được định hình ra sao? Lộ trình 10 năm bao gồm: Giai đoạn 1 (1-3 năm): Chuyển giao mô hình giải tích sang bài toán dao động phi tuyến và va chạm xung kích; Giai đoạn 2 (3-6 năm): Thực nghiệm chế tạo mẫu in 3D FGP và nén thử nghiệm kiểm chứng mô hình lý thuyết; Giai đoạn 3 (6-10 năm): Tích hợp trí tuệ nhân tạo (AI/Machine Learning) để tối ưu hóa hình học vi lỗ xốp phục vụ chế tạo thông minh cho ngành hàng không vũ trụ.

Kết luận

  1. Luận án đã xây dựng thành công mô hình giải tích phi tuyến hoàn chỉnh cho panel trụ và vỏ trụ làm bằng vật liệu xốp FGP và sandwich FGP chịu tải trọng cơ - nhiệt phức tạp trên nền đàn hồi Pasternak.
  2. Thiết lập được hệ công thức dạng hiển (explicit analytical formulas) xác định chính xác tải trọng tới hạn trên $r_{0,\text{upper}}$, tải trọng tới hạn dưới $r_{0,\text{lower}}$ và vẽ đường cong cân bằng sau tới hạn phản ánh hiện tượng mất ổn định kiểu cực trị snap-through.
  3. Chứng minh một cách khoa học rằng mô hình phân bố xốp đối xứng (Loại a) tối ưu nhất về khả năng chịu tải nén (vượt trội hơn từ $9.03%$ đến $11.26%$ so với phân bố không đối xứng và đồng đều).
  4. Định lượng hóa toàn diện ảnh hưởng của hệ số độ rỗng $e_0$, tỷ lệ chiều dày lõi $h_{\text{core}}/h_{\text{FG}}$, tham số hình học $a/b, a/R, R/h$, độ không hoàn hảo $w^*$, phản lực nền $K_1, K_2$, và mật độ gân gia cường đa hướng đến độ ổn định của vỏ.
  5. Mở ra 3 nhánh nghiên cứu mới giàu tiềm năng: ổn định động lực học vỏ xốp chịu tải xung kích nhiệt, tối ưu hóa topology lỗ xốp đa chức năng và cơ học phá hủy đàn dẻo của vật liệu FGP.
  6. Công trình khẳng định vị thế tiên phong của cơ học giải tích Việt Nam trên trường quốc tế, mang lại giá trị thực tiễn và lý luận lâu dài cho sự phát triển của ngành công nghiệp vật liệu và kết cấu composite tiên tiến.