Tổng quan về luận án

Nghiên cứu về ứng xử cơ học của các kết cấu thành mỏng làm bằng vật liệu tiên tiến là một trong những trụ cột then chốt của cơ kỹ thuật hiện đại. Luận án tiến sĩ kỹ thuật cơ khí và cơ kỹ thuật của tác giả Phạm Văn Hoàn (mã số: 9520101), thực hiện dưới sự hướng dẫn khoa học của PGS. Lê Khả Hòa và PGS. Đào Như Mai tại Viện Hàn lâm Khoa học và Công nghệ Việt Nam (2024), mang tiêu đề: "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". Công trình đặt trọng tâm vào phân tích giải tích phi tuyến trạng thái ổn định tĩnh và ứng xử sau mất ổn định (post-buckling) của các cấu trúc panel trụ và vỏ trụ làm bằng vật liệu xốp có cơ tính biến đổi (Functionally Graded Porous - FGP), kết cấu sandwich FGP ba lớp có gân gia cường FGM (gân dọc, gân vòng, gân trực giao, gân xiên) đặt trên nền đàn hồi hai tham số Pasternak chịu tổ hợp tải trọng cơ - nhiệt phức tạp.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế: mặc dù các kết cấu vật liệu có cơ tính biến đổi truyền thống (Functionally Graded Materials - FGM) đã được khảo sát rộng rãi bởi Shen (2009) [25] hay Reddy & Starnes (1993) [52], nhưng các công trình tiếp cận giải tích đối với vật liệu xốp biến đổi (FGP) dạng panel trụ và vỏ trụ sandwich ba lớp chịu tải nén dọc trục, áp lực ngoài, tải xoắn kết hợp gradient nhiệt độ trong điều kiện có kể đến độ không hoàn hảo hình học ban đầu ($w^*$) và gân xiên gia cường vẫn còn rất hạn chế.

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âu hỏi nghiên cứu 1 (RQ1): Quy luật phân bố lỗ rỗng 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) và hệ số độ rỗng $e_0 \in [0, 1)$ chi phối như thế nào đến khả năng chịu lực tới hạn và quỹ đạo cân bằng sau mất ổn định của kết cấu?
  • Câu hỏi nghiên cứu 2 (RQ2): Tác động tương hỗ giữa cấu trúc sandwich FGP ba lớp (lớp lõi xốp $h_{\text{core}}$ kẹp giữa hai lớp mặt FGM $h_{\text{FG}}/2$) và mật độ/hướng gân FGM gia cường (gân trực giao, gân xoắn $n_p$ góc $\psi$) biến đổi ra sao dưới các trường nhiệt độ $\Delta T$?
  • Câu hỏi nghiên cứu 3 (RQ3): Mức độ ảnh hưởng của tương tác nền đàn hồi Pasternak ($K_1, K_2$) và điều kiện biên (bốn cạnh tựa đơn - SSSS so với hai cạnh cong tựa đơn, hai cạnh thẳng ngàm - SSCC) đến tải trọng tới hạn trên và tải trọng tới hạn dưới?
  • Giả thuyết khoa học (H1 & H2): Luận án giả thuyết rằng mô hình phân bố xốp đối xứng sẽ bảo toàn độ cứng uốn hiệu dụng cao nhất so với phân bố xốp đồng đều khi cùng tổng thể tích độ rỗng; đồng thời, việc tích hợp gân xiên FGM tối ưu hóa góc nghiêng $\psi$ sẽ triệt tiêu hiệu ứng suy giảm độ cứng do nhiệt độ sinh ra dưới tác dụng của mô men xoắn thuần túy.

Khung lý thuyết của luận án được định hình từ sự giao thoa giữa lý thuyết biến dạng phi tuyến hình học von Kármán, lý thuyết vỏ Donnell mỏng, phương pháp san đều tác dụng gân Leckhnitsky (smeared stiffeners technique) và mô hình nền hai tham số Pasternak. Luận án mang lại đóng góp đột phá khi thiết lập được hệ phương trình vi phân phi tuyến tương thích chuyển vị - hàm ứng suất, giải trực tiếp bằng phương pháp Galerkin để thu được nghiệm giải tích dạng tường minh (closed-form solutions) cho tải tới hạn trên ($r_0^{\text{upper}}, P^*, q^{\text{cr}}, \tau^{\text{upper}}$) và tải tới hạn dưới ($r_0^{\text{lower}}, P_{\text{lower}}, \tau^{\text{lower}}$). Toàn bộ thuật toán được lập trình trên nền tảng MATLAB với phạm vi khảo sát tham số hình học rộng ($a/b, a/R, R/h, L/R$) và kiểm chứng chuẩn hóa với các công bố quốc tế uy tín.

Literature Review và Positioning

Lịch sử phát triển của cơ học kết cấu composite biến đổi cơ tính chứng kiến bước nhảy vọt từ FGM hai thành phần kinh điển do các nhà khoa học Nhật Bản khởi xướng năm 1984 [1] sang vật liệu xốp biến đổi FGP tiên tiến nhằm tối ưu tỷ lệ độ cứng trên khối lượng và khả năng tiêu tán năng lượng, cách âm, cách nhiệt (Chen et al., 2016 [68]; Chen, Kitipornchai & Yang, 2016 [69]).

                               TIẾN TRÌNH LITERATURE VÀ ĐỊNH VỊ LUẬN ÁN
                               
 1984: Khởi xướng FGM              1996 - 2003: Phân tích FGM           2016 - Nay: Kỷ nguyên FGP
 (Vật liệu gốm - kim loại)         (Shen, Jaunky, Turvey)               (Chen, Keleshteri, Zhenhuan)

Trong bức tranh tổng quan học thuật, các dòng nghiên cứu lớn tập trung vào phân tích ổn định kết cấu dạng vỏ mỏng FGM có thể phân loại qua hai xu hướng phương pháp luận đối nghịch:

  1. Trường phái số trị và bán giải tích (FEM & Semi-analytical): Kweon & Hong (1994) [41] dùng phần tử vỏ suy biến 8 nút Lagrange; Zghal et al. (2018) [31] dùng FEM phân tích vỏ gia cường ống nano carbon; García-Macías et al. (2017) [32, 33] áp dụng mô hình Eshelby-Mori-Tanaka; trong khi Yang et al. (2004) [43] và Alijani & Aghdam (2012) [44] áp dụng phương pháp vi phân cầu phương và Kantorovich mở rộng. Ưu điểm của trường phái này là xử lý linh hoạt biên phức tạp nhưng hạn chế lớn là chi phí tính toán cao và không cung cấp được các công thức tham số hiển để phục vụ tối ưu hóa thiết kế kết cấu nhanh.
  2. Trường phái giải tích phi tuyến giải nghiệm hàm đóng (Analytical Closed-form): Đại diện tiêu biểu là Shen & Liew (2003, 2004) [23, 24], Shen (2003, 2009) [25, 53] sử dụng phương pháp nhiễu bậc cao kết hợp lý thuyết biến dạng cắt bậc cao (HSDT); Hoàng Văn Tùng và Nguyễn Đình Đức (2010, 2013) [34, 48] sử dụng lý thuyết biến dạng cắt bậc nhất kết hợp Galerkin. Nhóm các nhà cơ học Việt Nam như Đào Văn Dũng, Lê Khả Hòa, Vũ Hoài Nam, Hoàng Thị Thiêm (2014–2021) [50, 54, 62–67, 81] đã phát triển hệ phương trình giải tích màng Donnell và von Kármán cho vỏ FGM có gân gia cường trực giao và xiên.

Tranh luận học thuật cốt lõi diễn ra giữa việc mô hình hóa gân gia cường: phương pháp gân rời rạc cục bộ (discrete stiffener approach) phản ánh chính xác trường ứng suất tại chân gân nhưng phức tạp hóa bài toán vi phân, trong khi kỹ thuật san đều gân Leckhnitsky (smeared stiffener technique) cho phép tính toán toàn cục ổn định kết cấu với độ chính xác cao đối với mật độ gân dày. Về mặt vật liệu xốp, các công trình của Keleshteri & Jelovica (2020) [70] chỉ giới hạn ở dao động phi tuyến tấm; Zhenhuan et al. (2020) [76] phân tích vỏ trụ composite gia cường graphene xốp dưới tải nén thuần túy; Dong et al. (2019) [79] nghiên cứu vỏ quay nanocompozit; Belica & Magnucki (2020) [80] phân tích ổn định vỏ bọt kim loại theo nguyên lý thế năng toàn phần nhưng bỏ qua gân xiên gia cường và gradient nhiệt độ.

Luận án của NCS. Phạm Văn Hoàn định vị chính xác tại điểm giao thoa chưa được khai phá: xây dựng khung giải tích phi tuyến toàn diện cho vỏ trụ và panel trụ sandwich FGP ba lớp có gân gia cường đa hướng (trực giao, xiên) chịu đa tải trọng (nén trục, áp lực ngoài, xoắn thuần túy, nhiệt độ đồng đều $\Delta T$) trên nền Pasternak ($K_1, K_2$). Độ chính xác của mô hình giải tích được đối chuẩn thành công với hai nghiên cứu quốc tế kinh điển:

  • Đối chuẩn với Jaunky & Knight (1999) [39] và Shen (2003) [25] về tải tới hạn nén dọc trục $N_{\text{cr}}$ của panel trụ đẳng hướng ($a/b=0.5, a/R=0.2, b/h=500$): kết quả của luận án đạt $N_{\text{cr}} = 0.0384$, hoàn toàn trùng khớp với giá trị $0.0384$ của Shen [25] và mô hình FEM/HSDT của Jaunky [39] (sai số < 0.01%).
  • Đối chuẩn với Turvey (1996) [38] và Shen (2003) [25] về tải không thứ nguyên $N_x / (E R)$: với thông số $a/b=0.2, a/R=0.1, b/h=500$, kết quả luận án đạt $0.0553 \times 10^{-5}$ ứng với mode sóng tới hạn $(m,n)=(4,2)$, tái khẳng định tính chuẩn xác và độ tin cậy tuyệt đối của thuật toán giải tích.

Đó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 biên giới của cơ học vật rắn biến dạng phi tuyến qua các đóng góp nền tảng:

  • Mở rộng lý thuyết màng Donnell và quan hệ biến dạng phi tuyến von Kármán: Thiết lập thành công mô hình cấu nạp ứng suất - biến dạng tổng quát cho vật liệu composite xốp FGP dị hướng cục bộ theo chiều dày, tích hợp hàm phân bố độ rỗng phi tuyến vào ten-xơ độ cứng uốn ($D_{ij}$) và độ cứng kéo nén ($A_{ij}$).
  • Khung khái niệm và hệ phương trình vi phân liên kết: Thiết lập hệ phương trình vi phân phi tuyến bậc bốn đối với hai trường biến số độc lập: độ võng uốn $w(x,y)$ và hàm ứng suất mặt giữa $f(x,y)$ theo dạng: $$\mathcal{A}3 \nabla^4 f - \frac{1}{R} f{,xx} + \mathcal{A}_4 \nabla^4 w + \mathcal{L}(w+w^, f) = 0$$ $$\mathcal{A}1 \nabla^4 w + \frac{1}{R} w{,xx} - \mathcal{A}_2 \nabla^4 f + \mathcal{M}(w+w^, w) = 0$$
  • Hệ thống mệnh đề khoa học (Theoretical Propositions):
    1. Mệnh đề 1: Độ rỗng phân bố đối xứng (Loại a) bảo toàn thế năng biến dạng uốn cực đại, dẫn đến khả năng chống mất ổn định vượt trội so với phân bố bất đối xứng (Loại b, c) và phân bố đồng đều (Loại d).
    2. Mệnh đề 2: Sự suy giảm độ cứng nhiệt độ ($\Delta T$) tác động phi tuyến lên tải tới hạn dưới ($r_0^{\text{lower}}, \tau^{\text{lower}}$) mạnh hơn đáng kể so với tải tới hạn trên ($r_0^{\text{upper}}, \tau^{\text{upper}}$).
    3. Mệnh đề 3: Bố trí gân xiên xoắn ốc ($n_p$) có khả năng triệt tiêu mode mất ổn định xoắn bất đối xứng, biến đổi quỹ đạo rẽ nhánh cân bằng từ dạng nhảy vọt (snap-through) sang dạng suy biến êm.

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

Khung phân tích của luận án là sự hợp nhất của ba trụ cột giải tích:

  1. Quy luật biến đổi vật liệu sandwich xốp FGP:
    • Lớp mặt FGM ($z \in [h_{\text{core}}/2, h/2]$ và $z \in [-h/2, -h_{\text{core}}/2]$) biến thiên theo hàm lũy thừa thể tích $V_c(z) = (2z/h)^k$ kết hợp tổ hợp vật liệu gốm - kim loại ($\text{ZrO}_2/\text{Ti-6Al-4V}, \text{Si}_3\text{N}_4/\text{SUS304}$).
    • Lớp lõi xốp FGP ($z \in [-h_{\text{core}}/2, h_{\text{core}}/2]$) tuân theo bốn mô hình độ xốp:
      • Phân bố đối xứng (Loại a): $E(z) = E_m \left[1 - e_0 \cos\left(\frac{\pi z}{h_{\text{core}}}\right)\right]$
      • Phân bố không đối xứng (Loại b, c): $E(z) = E_m \left[1 - e_0 \cos\left(\frac{\pi z}{2h} \pm \frac{\pi}{4}\right)\right]$
      • Phân bố đồng đều (Loại d): $E(z) = E_m [1 - e_0 \alpha]$
  2. Kỹ thuật san đều gân Leckhnitsky: Tích phân liên tục hóa độ cứng của hệ gân dọc ($s$), gân vòng ($r$) và gân xiên ($p$) với các kích thước hình học $(h_s, b_s, d_s), (h_r, b_r, d_r), (h_p, b_p, d_p)$ vào hệ ma trận nội lực tổng quát của vỏ.
  3. Phương pháp xấp xỉ liên tục Galerkin: Sử dụng các hàm dạng lượng giác thỏa mãn chính xác các điều kiện biên động học và lực học (SSSS và SSCC), biến đổi hệ phương trình đạo hàm riêng phi tuyến thành hệ phương trình đại số phi tuyến bậc ba, cho phép giải tích hóa toàn bộ quan hệ tải trọng - độ võng $(r_0 - W/h, q - W_{\text{max}}/h, \tau - \psi)$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo thế giới quan thực chứng duy lý (Positivist / Rationalist Mechanics paradigm), dựa trên logic diễn dịch toán lý cơ học môi trường liên tục. Luận án tích hợp mô hình hóa đa tỉ lệ:

  • Tỉ lệ vi mô/vật liệu: Mô hình hóa vi cấu trúc lỗ rỗng phân bố liên tục và gradient cơ tính FGM qua chiều dày.
  • Tỉ lệ cấu trúc: Mô hình hóa tương tác giữa panel/vỏ cong, hệ thống gân gia cường bố trí cách đều và lớp nền đàn hồi Pasternak.
  • Tỉ lệ môi trường: Tích hợp trường tải trọng nén, áp lực, lực xoắn cắt cùng trường nhiệt độ biến thiên đồng đều $\Delta T$.
               QUY TRÌNH PHÂN TÍCH GIẢI TÍCH VÀ MÔ PHỎNG SỐ

Quy trình nghiên cứu rigorous

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

  1. Thiết lập quan hệ biến dạng - chuyển vị: Kể đến số hạng phi tuyến bậc hai von Kármán: $$\epsilon_x^0 = u_{,x} + \frac{1}{2} w_{,x}^2 + w_{,x} w_{,x}^, \quad \epsilon_y^0 = v_{,y} - \frac{w}{R} + \frac{1}{2} w_{,y}^2 + w_{,y} w_{,y}^, \quad \gamma_{xy}^0 = u_{,y} + v_{,x} + w_{,x} w_{,y} + w_{,x} w_{,y}^* + w_{,y} w_{,x}^*$$
  2. Thiết lập phương trình tương thích biến dạng: Khử các thành phần chuyển vị trong mặt phẳng $u, v$ để tạo phương trình liên hệ giữa $\epsilon_x^0, \epsilon_y^0, \gamma_{xy}^0$ với độ võng $w$ và độ võng ban đầu $w^* = \xi h \sin(m\pi x/a)\sin(n\pi y/b)$.
  3. Cân bằng lực và mô men: Tích hợp phương trình Donnell: $$N_{x,x} + N_{xy,y} = 0, \quad N_{xy,x} + N_{y,y} = 0$$ $$M_{x,xx} + 2M_{xy,xy} + M_{y,yy} + \frac{N_y}{R} + N_x (w_{,xx} + w_{,xx}^) + 2N_{xy} (w_{,xy} + w_{,xy}^) + N_y (w_{,yy} + w_{,yy}^*) - K_1 w + K_2 \nabla^2 w = 0$$
  4. Trực giao hóa Galerkin: Nhân phương trình với hàm trọng số tương ứng và lấy tích phân trên toàn bộ miền kết cấu $\int_0^a \int_0^b \Omega(x,y) dx dy = 0$.

Data và phân tích

Dữ liệu tính toán được thiết lập trên các hệ vật liệu kỹ thuật tiêu chuẩn cao cấp:

  • Hệ Zirconia / Ti-6Al-4V: $E_c = 168\text{ GPa}, \nu_c = 0.3, E_m = 105.7\text{ GPa}, \nu_m = 0.3$.
  • Hệ Silicon Nitride / Stainless Steel ($\text{Si}_3\text{N}_4/\text{SUS304}$): $E_c = 348.43\text{ GPa}, E_m = 201.04\text{ GPa}$.
  • Miền khảo sát tham số: Hệ số rỗng $e_0 \in [0, 0.8]$; tỷ lệ chiều dày lõi/vỏ $h_{\text{core}}/h_{\text{FG}} \in [0.2, 5.0]$; chỉ số thể tích FGM $k \in [0, 10]$; hệ số nền đàn hồi Winkler $K_1 \in [0, 10^8]\text{ N/m}^3$, Pasternak $K_2 \in [0, 10^6]\text{ N/m}$; gradient nhiệt độ $\Delta T \in [0, 300]\text{ K}$.
  • Cực tiểu hóa số sóng mất ổn định: Mọi giá trị tải tới hạn đều được quét tối ưu hóa qua không gian số sóng $(m, n) \in [1, 20] \times [1, 20]$ bằng thuật toán lặp trên MATLAB để đảm bảo tìm đúng mode tới hạn cực tiểu.

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

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

Luận án đã rút ra 5 phát hiện mang tính quy luật cơ học đột phá:

STT Phát hiện cơ học cốt lõi Dẫn chứng số liệu định lượng cụ thể từ luận án Ý nghĩa cơ học kết cấu
1 Ưu thế tuyệt đối của quy luật phân bố độ xốp đối xứng (Loại a) Với $e_0 = 0.5$, tải tới hạn trên đạt $r_0^{\text{upper}} = 293.7449\text{ MN}$ (Loại a), vượt trội hoàn toàn so với Loại b & c ($267.2059\text{ MN}$, cao hơn 9.93%) và Loại d ($260.6720\text{ MN}$, cao hơn 12.69%). Đặt lỗ rỗng tại trục trung hòa giúp duy trì mômen quán tính uốn của 2 lớp mặt ngoài.
2 Mức suy giảm phi tuyến của tải tới hạn theo hệ số rỗng $e_0$ Khi $e_0$ tăng từ $0 \to 0.5 \to 0.8$ ở mô hình đối xứng, $r_0^{\text{upper}}$ giảm từ $393.9360\text{ MN} \to 293.7449\text{ MN} \to 212.7248\text{ MN}$ (suy giảm 46.0%). Cung cấp định lượng chính xác cho đánh đổi giữa giảm tự trọng và độ cứng chịu lực.
3 Tác động kích kháng hình học ($a/R$ và $a/b$) Tại tỷ lệ $a/b = 2$, khi tăng độ cong $a/R$ từ $0.2 \to 0.9$, tải tới hạn $r_0^{\text{upper}}$ tăng vọt từ $132.9329\text{ MN} \to 400.3292\text{ MN}$ (tăng trưởng 201.15% hay gấp hơn 3 lần). Hiệu ứng màng do độ cong hình học áp chế biến dạng uốn cục bộ của panel mỏng.
4 Cơ chế suy thoái do gradient nhiệt độ $\Delta T$ Nhiệt độ tăng $\Delta T$ làm dịch chuyển toàn bộ đường cong cân bằng sau tới hạn $(p - W_{\text{max}}/h)$ xuống dưới, làm hạ thấp điểm tới hạn dưới $P_{\text{lower}}$ nhanh hơn $P_{\text{upper}}$. Ứng suất nén nhiệt phụ sinh làm suy yếu liên kết tinh thể và thúc đẩy mất ổn định sớm.
5 Hiệu ứng gia cường của gân xiên FGM ($n_p, \psi$) chịu tải xoắn Hệ gân xiên làm tăng tải xoắn tới hạn $\tau^{\text{upper}}$ lên gấp 1.8 – 2.5 lần so với vỏ trơn; hiệu quả gia cường đạt tối ưu tại góc nghiêng $\psi \approx 45^\circ$. Hướng gân xiên trùng với phương ứng suất kéo chính sinh ra bởi mômen xoắn. 
    TẢI TỚI HẠN TRÊN r₀ᵘᵖᵖᵉʳ (MN) THEO CÁC DẠNG PHÂN BỐ XỐP (e₀ = 0.5)
    
    Loại a (Đối xứng)    : ▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 293.74 MN
    Loại b, c (Bất đối)  : ▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 267.21 MN  (-9.9%)
    Loại d (Đồng đều)    : ▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 260.67 MN    (-12.7%)
    Không xốp (e₀ = 0)   : ▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇▇ 393.94 MN

Hiện tượng phản trực giác (counter-intuitive) được làm sáng tỏ: mặc dù cấu trúc xốp đồng đều (Loại d) có cùng mật độ rỗng trung bình như cấu trúc xốp đối xứng (Loại a), nhưng khả năng chịu nén của Loại d lại thấp nhất trong tất cả các cấu hình. Nguyên nhân là do ở Loại a, vật liệu đặc (độ rỗng bằng 0) tập trung tối đa tại hai bề mặt ngoài ($z = \pm h/2$), nơi chịu ứng suất uốn lớn nhất theo phân bố Euler-Bernoulli, do đó tối ưu hóa độ cứng chống uốn hiệu dụng mà không làm tăng khối lượng tổng thể. Ngoài ra, sự xuất hiện của độ không hoàn hảo hình học ban đầu $\xi$ làm triệt tiêu hoàn toàn điểm rẽ nhánh lý thuyết (bifurcation point), biến quá trình mất ổn định thành dạng uốn phi tuyến tiệm tiến ngay từ khi bắt đầu gia tải.

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp mô hình giải tích mẫu mực kết hợp đồng thời 4 yếu tố: FGP sandwich + Gân xiên/trực giao FGM + Nền Pasternak + Môi trường nhiệt độ, bổ sung hoàn thiện cho kho tàng lý thuyết vỏ phi tuyến thế giới.
  • Về mặt phương pháp luận: Bộ công thức hiển cho phép thay thế các chương trình mô phỏng FEM phi tuyến phức tạp trong giai đoạn thiết kế sơ bộ, rút ngắn thời gian tính toán từ nhiều giờ xuống vài phần nghìn giây.
  • Về ứng dụng thực tiễn: Định hình trực tiếp quy trình thiết kế cấu kiện vỏ tên lửa, thân máy bay siêu âm, vỏ tàu ngầm và kết cấu chịu áp lực làm việc trong điều kiện nhiệt độ cao, nơi yêu cầu tối ưu hóa tự trọng nghiêm ngặt.
  • Khuyến nghị kỹ thuật: Ưu tiên chế tạo lõi xốp theo phân bố đối xứng (Loại a) với hệ số $e_0 \le 0.5$; bố trí gân xiên FGM với góc nghiêng $45^\circ$ cho các bộ phận chịu tải xoắn chính; và chủ động tăng cường hệ số nền $K_2$ (lớp cắt Pasternak) để gia tăng tải tới hạn sau mất ổn định.

Limitations và Future Research

Luận án thừa nhận một số giới hạn học thuật và điều kiện biên thực nghiệm:

  1. Giới hạn lý thuyết vỏ Donnell: Phù hợp tối ưu cho vỏ mỏng và độ cong thoải ($h/R \ll 1$); đối với các kết cấu vỏ dày, hiệu ứng biến dạng trượt ngang qua chiều dày cần được mô tả bằng lý thuyết biến dạng cắt bậc cao (HSDT) hoặc lý thuyết đàn hồi 3D Biot.
  2. Giả thiết liên kết lớp hoàn hảo: Bỏ qua hiện tượng tách lớp (delamination) cục bộ tại mặt tiếp giáp giữa lớp phủ FGM và lõi xốp FGP dưới tải trọng nhiệt chu kỳ.
  3. Phạm vi tải trọng tĩnh: Luận án tập trung vào phân tích ổn định tĩnh, chưa khảo sát mất ổn định động dưới tác dụng của sóng xung kích (blast loading) hoặc tải trọng cơ học thay đổi theo thời gian.
  4. Mô hình nhiệt độ đồng đều: Giả thiết trường nhiệt độ phân bố đều $\Delta T$ qua chiều dày; thực tế các kết cấu FGM chịu gradient nhiệt độ phi tuyến dạng dừng hoặc không dừng.
                    CHƯƠNG TRÌNH NGHIÊN CỨU 10 NĂM TỚI

Chương trình nghiên cứu phát triển tiếp theo bao gồm:

  • Mở rộng phân tích ổn định động và đáp ứng quá độ của vỏ FGP sandwich dưới tác dụng của tải va đập và sóng nổ.
  • Khảo sát mất ổn định nhiệt đàn hồi phi tuyến với gradient nhiệt độ thay đổi phức tạp theo quy luật truyền nhiệt Fourier và phi-Fourier.
  • Phát triển mô hình phần tử hữu hạn nâng cao kết hợp cơ học phá hủy để mô phỏng sự lan truyền vết nứt và tách lớp trong lõi xốp.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning) và thuật toán tối ưu hóa đa mục tiêu để thiết kế tối ưu vi cấu trúc lỗ rỗng và góc nghiêng gân xiên.

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

  • Tác động học thuật (Academic Impact): Luận án đóng góp hệ phương trình và các biểu thức giải tích chuẩn mực, đóng vai trò làm dữ liệu chuẩn (benchmark data) tin cậy để các nhà nghiên cứu quốc tế kiểm chứng các thuật toán số mới (FEM, XFEM, IGA, Meshless). Ước tính các công bố từ luận án có khả năng thu hút trích dẫn cao trong các tạp chí ISI/Scopus hàng đầu nhóm Q1 ngành Cơ học và Kỹ thuật Kết cấu.
  • Chuyển đổi công nghiệp (Industry Transformation): Đặt nền tảng tính toán vững chắc cho các ngành công nghiệp công nghệ cao:
    • Hàng không vũ trụ: Chế tạo vỏ động cơ phản lực, thân tên lửa đẩy và cánh tàu con thoi bằng hợp kim xốp titan/gốm.
    • Quốc phòng và hàng hải: Thiết kế thân vỏ tàu ngầm chịu áp lực sâu, tấm chắn chống đạn tiêu thấu năng lượng va chạm.
    • Xây dựng & năng lượng: Ứng dụng bọt nhôm FGP làm panel cách âm, cách nhiệt chống cháy cho các công trình siêu cao tầng và nhà máy điện hạt nhân.
  • Lợi ích xã hội và môi trường: Việc tối ưu hóa cấu trúc lỗ rỗng giúp giảm tới 30–45% trọng lượng kết cấu nhưng vẫn đảm bảo độ an toàn chịu lực, từ đó tiết kiệm nhiên liệu vận tải, giảm phát thải khí nhà kính và tối ưu hóa chi phí nguyên vật liệu quý hiếm.

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

                              BẢN ĐỒ ĐỐI TƯỢNG HƯỞNG LỢI
  1. Nghiên cứu sinh và Giảng viên ngành Cơ học: Khai thác bộ khung lý thuyết giải tích mẫu mực, hệ phương pháp Galerkin giải tích phi tuyến và mã nguồn MATLAB làm tài liệu giảng dạy chuyên đề cao học và tiến sĩ.
  2. Kỹ sư thiết kế R&D Hàng không - Vũ trụ - Cơ khí: Sử dụng trực tiếp hệ biểu thức đóng để tính toán nhanh tải trọng tới hạn trên/dưới trong pha thiết kế ý niệm (conceptual design), cắt giảm 80% thời gian chạy mô phỏng số.
  3. Cơ quan quản lý và Viện nghiên cứu Tiêu chuẩn hóa: Tiếp cận cơ sở khoa học định lượng để biên soạn các quy chuẩn kỹ thuật quốc gia về thử nghiệm và ứng dụng vật liệu composite xốp biến đổi trong quốc phòng và dân dụng.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc thiết lập thành công hệ phương trình phi tuyến giải tích giải theo độ võng - hàm ứng suất cho vỏ trụ sandwich ba lớp FGP tích hợp hệ gân FGM gia cường xiên xoắn ốc trên nền Pasternak trong môi trường nhiệt. Công trình đã mở rộng trực tiếp lý thuyết vỏ Donnell và quan hệ biến dạng von Kármán truyền thống sang lớp vật liệu xốp biến đổi phi tuyến 3D, vượt qua các giới hạn của mô hình màng đẳng hướng thuần túy.

2. Phương pháp nghiên cứu có tính đột phá gì so với các nghiên cứu quốc tế trước đây?

So với phương pháp phần tử hữu hạn thuần số của Jaunky & Knight (1999) [39] hay phương pháp nhiễu phức tạp của Shen (2003, 2009) [25, 53], luận án sử dụng phương pháp giải tích Galerkin trực giao để rút ra biểu thức hiển đại số tường minh cho cả 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}}$). Cách tiếp cận này loại bỏ hoàn toàn sự phụ thuộc vào ma trận phần tử lớn, cho phép phân tích tham số liên tục với độ chính xác cao và thời gian xử lý tức thời.

3. Phát hiện nào bất ngờ nhất về mặt số liệu thực nghiệm cơ học?

Phát hiện bất ngờ nhất là mô hình phân bố độ xốp đồng đều (Loại d) có khả năng chịu tải nén tới hạn kém nhất ($260.6720\text{ MN}$), thấp hơn mô hình xốp đối xứng Loại a ($293.7449\text{ MN}$) tới 12.69% và thấp hơn phân bố bất đối xứng Loại b, c ($267.2059\text{ MN}$) dù cùng chia sẻ tổng thể tích lỗ rỗng. Điều này khẳng định quy luật phân bố vị trí lỗ rỗng theo chiều dày đóng vai trò quyết định hơn nhiều so với tổng hàm lượng rỗng của kết cấu.

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

Hoàn toàn đầy đủ. Luận án cung cấp chi tiết:

  • Toàn bộ các hệ số giải tích $A_1, A_2, A_3, A_4$ và các tích phân độ cứng $D_{ij}$.
  • Công thức hiển của các hệ số Galerkin $S_1, S_2, S_3, S_4, S_5$ và $C_1, C_2, C_3, C_4$.
  • Cấu trúc thuật toán MATLAB quét cực tiểu hóa số sóng $(m, n)$ và vẽ đường cong quan hệ tải - độ võng.

5. Chương trình nghiên cứu 10 năm tiếp theo được vạch ra như thế nào?

Luận án định hướng lộ trình nghiên cứu 10 năm tập trung vào: (1) Ổn định động phi tuyến dưới tải nổ và va chạm vận tốc cao; (2) Truyền nhiệt phi-Fourier kết hợp ứng suất nhiệt quá độ; (3) Mô hình hóa tách lớp giao diện sandwich và lan truyền vết nứt 3D; (4) Tối ưu hóa cấu trúc FGP bằng trí tuệ nhân tạo và công nghệ in 3D kim loại.

Kết luận

Luận án tiến sĩ của NCS. Phạm Văn Hoàn là một công trình nghiên cứu khoa học xuất sắc, hoàn chỉnh và có tính tiên phong cao trong lĩnh vực Cơ kỹ thuật, được đúc kết qua 6 đóng góp cốt lõi:

  1. 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ụ FGP, sandwich FGP ba lớp chịu tải nén dọc trục, áp lực ngoài, tải xoắn và tải nhiệt độ trên nền đàn hồi Pasternak.
  2. Thiết lập hệ biểu thức hiển tường minh xác định chính xác tải trọng tới hạn trên, tải trọng tới hạn dưới và phương trình đường cong cân bằng sau mất ổn định.
  3. Phát hiện và chứng minh quy luật tối ưu của phân bố xốp đối xứng (Loại a), mang lại khả năng chịu lực vượt trội 9.9% – 12.7% so với các dạng phân bố khác.
  4. Làm sáng tỏ vai trò kháng xoắn của hệ gân xiên FGM, chứng minh hiệu quả nâng cao tải tới hạn xoắn lên tới 2.5 lần tại góc bố trí tối ưu $\psi = 45^\circ$.
  5. Đánh giá toàn diện ảnh hưởng của nền đàn hồi Pasternak và điều kiện biên (SSSS, SSCC), chỉ ra vai trò nâng cao độ ổn định sau tới hạn của lớp cắt $K_2$.
  6. Cung cấp bộ dữ liệu chuẩn hóa và công cụ tính toán MATLAB tin cậy, mở ra hướng phát triển ứng dụng vật liệu composite xốp biến đổi trong các ngành công nghiệp hàng không, quốc phòng và xây dựng hiện đại trên phạm vi quốc tế.