Tổng quan về luận án

Sự phát triển vượt bậc của khoa học vật liệu trong thế kỷ XXI đã thúc đẩy sự ra đời của vật liệu có cơ tính biến thiên rỗng (Functionally Graded Porous Material - FGPM), một biến thể tiên tiến của vật liệu FGM truyền thống. Trong kết cấu FGPM, các lỗ rỗng (porosity) phân bố liên tục và biến đổi trơn theo phương chiều dày, cho phép loại bỏ hoàn toàn bề mặt phân chia pha đột ngột – nguyên nhân gốc rễ gây ra hiện tượng tập trung ứng suất và bong tách lớp (delamination) thường gặp ở composite đa lớp truyền thống. Đặc biệt, cấu trúc bọt kim loại tế bào mở (open-cell metal foam) mang lại tỷ lệ độ cứng trên trọng lượng vượt trội; thực nghiệm chứng minh "tấm sandwich dày 16mm với lớp bề mặt dày 1 mm, lõi bọt thép dày 14mm có độ cứng chống uốn bằng tấm thép dày 10mm, trong khi giảm được 35% khối lượng", đồng thời cung cấp khả năng hấp thụ năng lượng va chạm và tiêu tán chấn động tuyệt vời.

Mặc dù tiềm năng ứng dụng trong các ngành công nghiệp mũi nhọn như hàng không vũ trụ (nguyên mẫu tên lửa Ariane 5), công nghiệp ô tô đua, vỏ tàu thủy và kỹ thuật công trình xây dựng là rất lớn, nền tảng lý thuyết cơ học tính toán cho kết cấu tấm FGPM vẫn tồn tại những khoảng trống nghiên cứu (research gap) cốt lõi. Phần lớn các công bố quốc tế trước đây (chẳng hạn như Magnucki và cs. [71], Magnucka-Blandzi [70], Rezaei và Saidi [91]) chủ yếu tập trung vào phân tích tĩnh và dao động tuyến tính, hoặc tiến hành phân tích trên hệ tọa độ gắn với mặt trung bình hình học (geometrical mid-plane). Khi mật độ lỗ rỗng phân bố bất đối xứng qua chiều dày, mặt trung hòa vật lý (physical neutral surface) sẽ bị lệch khỏi mặt trung bình, dẫn đến sự xuất hiện của ma trận liên kết màng - uốn - xoắn ($B_{ij} \neq 0$), khiến hệ phương trình vi phân chủ đạo trở nên cực kỳ cồng kềnh, phức tạp và tiêu tốn tài nguyên tính toán khi giải các bài toán phi tuyến hình học.

Luận án tiến sĩ "Phân tích phi tuyến ứng xử tĩnh và ổn định của tấm bằng vật liệu FGM rỗng" của tác giả Lê Thanh Hải, chuyên ngành Cơ kỹ thuật (Mã số: 9520101) tại Trường Đại học Xây dựng Hà Nội, dưới sự hướng dẫn khoa học của GS. Trần Minh Tú và GS. Lê Xuân Huỳnh, đã giải quyết trọn vẹn bài toán cơ học phi tuyến phức tạp này. Luận án đặt ra và giải quyết 4 câu hỏi nghiên cứu cốt lõi:

  1. RQ1: Làm thế nào để giải trừ hoàn toàn tương tác màng - uốn phức tạp trong phương trình cấu tạo của tấm FGPM có độ rỗng biến thiên bất đối xứng?
  2. RQ2: Ứng xử uốn phi tuyến độ võng lớn của tấm FGPM trên nền đàn hồi chịu tải trọng phân bố ngang phụ thuộc như thế nào vào quy luật phân bố lỗ rỗng và điều kiện biên?
  3. RQ3: Tải trọng tới hạn mất ổn định tĩnh ($N_{cr}$) và đường cong cân bằng sau ổn định (post-buckling) chịu sự chi phối ra sao dưới tác động của trường nén trong mặt phẳng và độ không hoàn hảo hình học ban đầu ($w_0^*$)?
  4. RQ4: Sự tương tác định lượng giữa hệ số nền Pasternak ($K_w, K_s$), tỷ số kích thước ($a/h, b/a$) và hệ số rỗng ($e_0$) diễn ra theo cơ chế cơ học nào?

Hệ thống giả thuyết khoa học tương ứng được xác lập:

  • H1: Việc thiết lập hệ quy chiếu đặt tại mặt trung hòa vật lý ($z_{th} = z_{tb} - C$) sẽ triệt tiêu hoàn toàn thành phần ma trận độ cứng liên kết ($B_{ij} \equiv 0$), đưa hệ phương trình phi tuyến về dạng tách rời tương đương với vật liệu đẳng hướng.
  • H2: Quy luật phân bố lỗ rỗng dạng đối xứng (Dạng 2 - tâm xốp, biên đặc) sẽ duy trì khả năng chịu uốn và độ cứng chống mất ổn định tối ưu hơn đáng kể so với dạng phân bố đều (Dạng 1) và phân bố bất đối xứng (Dạng 3).
  • H3: Thành phần biến dạng phi tuyến von Kármán kết hợp với độ không hoàn hảo ban đầu sẽ làm suy giảm độ cứng tiền tới hạn nhưng tạo ra quỹ đạo sau ổn định ổn định (stable post-buckling path) đối với các điều kiện biên kiềm chế chuyển vị màng.
  • H4: Độ cứng cắt của nền Pasternak ($K_s$) đóng vai trò kiểm soát độ võng và nâng cao tải trọng tới hạn hiệu quả hơn độ cứng uốn Winkler ($K_w$).

Khung lý thuyết của nghiên cứu tích hợp Lý thuyết biến dạng cắt bậc nhất (First-order Shear Deformation Theory - FSDT) với 5 ẩn chuyển vị ($u_0, v_0, w_0, \phi_x, \phi_y$) cho tấm dày ($a/h = 10$) và Lý thuyết tấm cổ điển (Classical Plate Theory - CPT) với 3 ẩn chuyển vị cho tấm mỏng ($a/h = 50$), kết hợp trường biến dạng phi tuyến von Kármán và mô hình nền đàn hồi hai thông số Pasternak ($K_w, K_{sx}, K_{sy}$). Quy mô khảo sát bao quát toàn diện các tỷ số kích thước $a/h \in [10, 50]$, tỷ số cạnh $b/a \in [0.5, 2.0]$, hệ số rỗng $e_0 \in [0, 0.9]$, biên tựa khớp (SSSS), ngàm (CCCC), hỗn hợp (SCSC), cùng hai trạng thái liên kết màng di động (FM) và bất động (IM).

Literature Review và Positioning

Lịch sử nghiên cứu vật liệu cơ tính biến thiên khởi nguồn từ các nhà khoa học Nhật Bản vào năm 1984 (Niino, Koizumi [19]) với mục tiêu ban đầu là tạo ra vật liệu rào cản nhiệt cho cấu trúc hàng không vũ trụ. Trải qua bốn thập kỷ, lý thuyết cơ học tấm vỏ FGM đã được mở rộng mạnh mẽ bởi các học giả hàng đầu thế giới: Reddy [86, 88] đặt nền móng cho các mô hình giải tích và phần tử hữu hạn dựa trên lý thuyết biến dạng cắt bậc cao; Zenkour [132, 134] phát triển các lời giải tĩnh và động lực học giải tích; Liew và cs. [66, 136] tiên phong trong các thuật toán số không lưới; Kiani và cs. [15, 48] đóng góp sâu rộng về cơ nhiệt và độ ổn định kết cấu. Shen [99] trong tài liệu chuyên khảo kinh điển đã hệ thống hóa các kỹ thuật giải tích phi tuyến, dao động và ổn định cho tấm vỏ composite tiên tiến.

Khi cấu trúc vật liệu phát triển sang bọt kim loại xốp (FGPM), Magnucki và cs. [71] là những người đầu tiên đề xuất mô hình toán học cho dầm và tấm FGPM sử dụng hàm chuyển vị bậc ba để khảo sát uốn và ổn định tuyến tính. Tiếp đó, Rezaei và Saidi [91, 92] công bố lời giải chính xác cho dao động tự do của tấm chữ nhật FGPM bão hòa chất lỏng dựa trên lý thuyết biến dạng cắt bậc ba với điều kiện biên Levy. Ebrahimi và cs. [36, 37] khảo sát tần số dao động riêng của vỏ trụ bọt kim loại và tấm FGPM trên nền đàn nhớt. Trong nước, nhóm nghiên cứu của Đào Văn Dũng, Nguyễn Đình Đức [30, 33, 34, 35, 118, 119] đã đóng góp nhiều công trình giải tích xuất sắc về ổn định phi tuyến của tấm FGM và tấm FGM có gân gia cường chịu tải cơ - nhiệt bằng phương pháp hàm ứng suất Airy kết hợp thuật toán Galerkin.

Tuy nhiên, bức tranh học thuật quốc tế tồn tại hai luồng quan điểm và phương pháp luận đối lập:

  1. Trường phái phân tích trên mặt trung bình hình học: Đại đa số nghiên cứu (như Wu và cs. [124], Nam và cs. [78], Cong và cs. [21]) đặt hệ trục quy chiếu tại mặt giữa ($z=0$). Cách tiếp cận này dẫn đến việc tồn tại các số hạng ghép nối màng - uốn ($B_{ij} \neq 0$), khiến hệ phương trình vi phân phi tuyến bậc cao có các hệ số phụ thuộc chặt chẽ vào độ lệch cơ tính, đòi hỏi phải giải hệ phương trình đạo hàm riêng liên kết vô cùng phức tạp hoặc áp dụng kỹ thuật nhiễu loạn hai bước (two-step perturbation technique) với khối lượng tính toán đồ sộ.
  2. Trường phái phân tích trên mặt trung hòa vật lý: Các nghiên cứu tiên phong của Prakash và cs. [81] (dùng FEM) và Kim & Lee [53] (dùng Isogeometric Analysis - IGA) chỉ ra rằng việc chuyển hệ quy chiếu về mặt trung hòa thực sẽ loại bỏ hoàn toàn thành phần ghép nối màng - uốn. Tuy nhiên, các công trình này chỉ dừng lại ở phương pháp số xấp xỉ, chưa từng xây dựng được hệ phương trình giải tích tường minh khép kín cho bài toán uốn phi tuyến và sau ổn định của tấm FGPM.
TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT CƠ HỌC KẾT CẤU TẤM FGM VÀ FGPM
===================================================================================
1984: FGM truyền thống (Niino, Koizumi) -> Chịu nhiệt kim loại/gốm
1990s-2000s: Cơ học tấm FGM tuyến tính & phi tuyến (Reddy, Zenkour, Shen, Đức)
2010s: Vật liệu FGM rỗng - FGPM (Magnucki, Rezaei & Saidi, Ebrahimi)
2021: LUẬN ÁN LÊ THANH HẢI (ĐH Xây dựng Hà Nội)
===================================================================================

Vị thế học thuật của luận án được xác lập vững chắc khi so sánh trực tiếp với các công trình quốc tế điển hình:

  • So với nghiên cứu của Magnucki và cs. [71]: Công trình của Magnucki chỉ xét bài toán uốn tuyến tính đơn thuần của tấm tròn/chữ nhật FGPM và bỏ qua tương tác nền đàn hồi cũng như yếu tố phi tuyến hình học von Kármán. Luận án của Lê Thanh Hải đã mở rộng vượt bậc sang miền ứng xử phi tuyến lớn, thiết lập đường cong sau ổn định và khảo sát chi tiết ảnh hưởng của độ không hoàn hảo hình học ban đầu.
  • So với nghiên cứu của Rezaei và Saidi [91]: Nhóm tác giả Saidi giải bài toán dao động riêng của tấm FGPM dựa trên nghiệm giải tích Levy nhưng bị giới hạn trong điều kiện biên có hai cạnh đối diện tựa khớp và hoàn toàn không đề cập đến bài toán uốn tĩnh phi tuyến hay phá hủy mất ổn định rẽ nhánh. Luận án này giải quyết đa dạng các điều kiện biên phức tạp (SSSS, CCCC, SCSC) dưới cả hai tiếp cận ứng suất và chuyển vị.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã đóng góp một bước tiến lý thuyết cơ bản vào cơ học vật rắn biến dạng bằng việc hoàn thiện mô hình toán học giải tích cho tấm FGPM không hoàn hảo trên nền đàn hồi Pasternak. Đóng góp lý thuyết mang tính cốt lõi bao gồm:

  1. Chứng minh giải tích sự triệt tiêu tương tác màng - uốn: Bằng việc xác định chính xác khoảng cách $C$ từ mặt trung bình đến mặt trung hòa vật lý thông qua điều kiện cân bằng tĩnh học: $$\int_{-h/2}^{h/2} (z_{tb} - C) E(z_{tb}) dz_{tb} = 0 \implies C = \frac{\int_{-h/2}^{h/2} z_{tb} E(z_{tb}) dz_{tb}}{\int_{-h/2}^{h/2} E(z_{tb}) dz_{tb}}$$ Luận án đã chứng minh một cách tường minh rằng trong hệ tọa độ $z_{th} = z_{tb} - C$, các thành phần độ cứng ghép nối: $$B_{ij} = \int_{-h/2 - C}^{h/2 - C} Q_{ij}(z_{th}) z_{th} dz_{th} = \int_{-h/2}^{h/2} (z_{tb} - C) Q_{ij}(z_{tb}) dz_{tb} \equiv 0 \quad (ij = 11, 12, 66)$$ Điều này cho phép tái cấu trúc toàn bộ hệ phương trình đàn hồi phi tuyến, đưa ma trận độ cứng phức tạp của vật liệu không đồng nhất về dạng cấu trúc tương tự vật liệu đẳng hướng đồng nhất, loại bỏ hoàn toàn các số hạng ghép nối giả tạo sinh ra do việc chọn gốc tọa độ hình học.

  2. Hệ thống mệnh đề lý thuyết (Theoretical Propositions):

  • Mệnh đề 1 (P1): Đối với tấm FGPM phân bố lỗ rỗng đối xứng (Dạng 2), vị trí mặt trung hòa trùng với mặt trung bình ($C = 0$). Đối với phân bố bất đối xứng (Dạng 3), $C \neq 0$ và độ lệch tâm tỷ lệ thuận với hệ số rỗng $e_0$.
  • Mệnh đề 2 (P2): Hệ phương trình vi phân cân bằng phi tuyến uốn theo FSDT bao gồm 5 phương trình đạo hàm riêng phi tuyến ghép nối theo 5 ẩn ($u_0, v_0, w_0, \phi_x, \phi_y$), khi chiếu lên mặt trung hòa sẽ tách biệt hoàn toàn thành phần biến dạng uốn - trượt khỏi biến dạng kéo - nén màng thuần túy ở bước thiết lập nội lực.
  • Mệnh đề 3 (P3): Lực tới hạn đàn hồi $N_{cr}$ của tấm FGPM chịu nén hai trục $(\gamma_1, \gamma_2)$ là một hàm giải tích tường minh phụ thuộc vào độ cứng chống uốn $C_{11}$, độ cứng cắt ngang $A_{44}^s$, tỷ số nén và các tham số nền đàn hồi ($K_w, K_s$).
  • Mệnh đề 4 (P4): Quỹ đạo sau ổn định của tấm FGPM có độ không hoàn hảo ban đầu dạng $w_0^*(x, y) = \xi h \sin(\alpha x) \sin(\beta y)$ chuyển đổi trơn tru từ trạng thái uốn sang sau ổn định mà không xảy ra điểm rẽ nhánh đột ngột (bifurcation point), phản ánh đúng bản chất vật lý thực tế.

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

Khung phân tích của luận án tích hợp đồng thời ba cấu trúc lý thuyết nền tảng:

  • Lý thuyết động học chuyển vị: Tích hợp đồng thời mô hình FSDT của Mindlin (kể đến biến dạng trượt ngang với "hệ số hiệu chỉnh cắt $k_s = 5/6$ được sử dụng trong nghiên cứu này") và mô hình CPT của Kirchhoff-Love.
  • Cơ học phi tuyến chuyển vị lớn: Áp dụng quan hệ biến dạng - chuyển vị phi tuyến von Kármán: $$\varepsilon_x^0 = u_{0,x} + \frac{1}{2} w_{0,x}^2 + w_{0,x} w_{0,x}^, \quad \varepsilon_y^0 = v_{0,y} + \frac{1}{2} w_{0,y}^2 + w_{0,y} w_{0,y}^$$ $$\gamma_{xy}^0 = u_{0,y} + v_{0,x} + w_{0,x} w_{0,y} + w_{0,x} w_{0,y}^* + w_{0,y} w_{0,x}^*$$
  • Mô hình tương tác môi trường đàn hồi: Tích hợp mô hình nền Pasternak hai thông số với phản lực nền được mô tả đầy đủ: $$f_e = K_w w_0 - K_{sx} w_{0,xx} - K_{sy} w_{0,yy}$$

Điểm độc đáo vượt trội trong khung phân tích của luận án là việc triển khai song song hai cách tiếp cận giải tích:

  • Tiếp cận theo chuyển vị (Displacement Approach): Biểu diễn các thành phần chuyển vị ($u_0, v_0, w_0, \phi_x, \phi_y$) qua chuỗi lượng giác kép Navier/Fourier thỏa mãn các điều kiện biên động học, sau đó áp dụng thuật toán Bubnov-Galerkin để chiếu trực tiếp hệ phương trình vi phân đạo hàm riêng phi tuyến thành hệ phương trình đại số phi tuyến bậc ba.
  • Tiếp cận theo ứng suất (Stress/Airy Function Approach): Sử dụng hàm ứng suất Airy $F(x, y)$ thỏa mãn phương trình tương thích biến dạng, kết hợp với phương trình cân bằng theo phương vuông góc mặt tấm để rút gọn số ẩn, đưa bài toán ổn định và sau ổn định về biểu thức giải tích tường minh cực kỳ tiện lợi cho tính toán kỹ thuật.
                     KHUNG PHÂN TÍCH GIẢI TÍCH ĐỘC ĐÁO CỦA LUẬN ÁN
                     ===============================================
                                  MÔ HÌNH VẬT LIỆU FGPM
                           (Bọt kim loại: Dạng 1, Dạng 2, Dạng 3)
                           QUY CHIẾU MẶT TRUNG HÒA VẬT LÝ
                               (Xác định C -> Bij = 0)
       LÝ THUYẾT TẤM DÀY (FSDT)                                  LÝ THUYẾT TẤM MỎNG (CPT)
       (5 ẩn: u0, v0, w0, φx, φy; ks=5/6)                         (3 ẩn: u0, v0, w0; γxz=γyz=0)
                       TÍCH HỢP QUAN HỆ PHI TUYẾN VON KÁRMÁN
                       + NỀN PASTERNAK (Kw, Ks) + ĐỘ LỆCH w0*
       TIẾP CẬN CHUYỂN VỊ                                        TIẾP CẬN ỨNG SUẤT (AIRY)
    (Chuỗi kép Navier/Galerkin)                               (Hàm F(x,y) + Phương pháp Galerkin)
    - Phân tích uốn phi tuyến                                 - Phân tích uốn phi tuyến
    - Xác định độ võng w0, nội lực Mx, My, Mxy                - Lập biểu thức tường minh tải tới hạn Ncr
    - Áp dụng biên SSSS, CCCC, SCSC                           - Thiết lập đường cong sau ổn định P - w0
                                THUẬT TOÁN SỐ MATLAB
                     (Kiểm chứng Benchmarks & Khảo sát đa tham số)

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ế trên lập trường triết học thực chứng (Positivism) và chủ nghĩa duy lý diễn dịch (Deductive Rationalism) của cơ học giải tích. Quy trình nghiên cứu kết hợp chặt chẽ giữa toán học giải tích thuần túy (analytical formulation) và phương pháp số thực nghiệm (numerical experimentation).

Thiết kế nghiên cứu phân tầng rõ rệt:

  • Cấp độ mô hình hóa: Xây dựng phương trình biến phân từ Nguyên lý thế năng toàn phần cực tiểu $\delta(U_P + U_F - V) = 0$, trong đó $U_P$ là thế năng biến dạng đàn hồi trong thể tích tấm, $U_F$ là thế năng phản lực nền Pasternak, và $V$ là thế năng của ngoại lực bề mặt và lực nén biên.
  • Cấp độ giải thuật: Thiết lập hệ phương trình vi phân phi tuyến đạo hàm riêng có hệ số là các hàm số cơ tính biến thiên. Áp dụng kỹ thuật trực giao hóa Bubnov-Galerkin trên miền diện tích mặt trung hòa $A = [0, a] \times [0, b]$ để loại bỏ các biến không gian, biến đổi bài toán biên phi tuyến thành hệ phương trình phi tuyến phi thứ nguyên.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được kiểm soát tính chuẩn xác tuyệt đối qua các bước:

  1. Thiết lập quan hệ ứng suất - biến dạng biến thiên: Khảo sát ba mô hình phân bố cơ tính bọt kim loại:
    • Dạng 1 (Phân bố đều): $E(z) = E_1(1 - e_0)$
    • Dạng 2 (Phân bố không đều - đối xứng): $E(z_{th}) = E_1 [1 - e_0 \cos(\frac{\pi z_{th}}{h})]$
    • Dạng 3 (Phân bố không đều - bất đối xứng): $E(z_{th}) = E_1 [1 - e_0 \cos(\frac{\pi z_{th}}{2h} + \frac{\pi}{4})]$ với hệ số rỗng khối lượng tương ứng $e_m = 1 - \sqrt{1 - e_0}$.
  2. Kỹ thuật hàm dạng và điều kiện biên: Lựa chọn các hệ hàm lượng giác và hàm đại số thỏa mãn chính xác các điều kiện biên động học và tĩnh học:
    • Biên tựa khớp bốn cạnh (SSSS): Chuyển vị góc xoay và mô men uốn triệt tiêu tại biên.
    • Biên ngàm bốn cạnh (CCCC): Chuyển vị và độ xoay pháp tuyến bị khóa cứng hoàn toàn ($w_0 = 0, \phi_n = 0$).
    • Biên hỗn hợp (SCSC): Hai cạnh đối diện ngàm, hai cạnh đối diện khớp.
  3. Quy trình kiểm chứng đa giác độ (Triangulation & Convergence):
    • Đánh giá sự hội tụ của chuỗi lượng giác kép với số số hạng gia tăng ($M = N = 1, 3, 5, 7$). Kết quả chỉ ra chuỗi đạt độ hội tụ chính xác tuyệt đối tại $M = N = 3$ với sai số dưới $0.05%$.
    • Kiểm chứng chéo giữa nghiệm theo tiếp cận chuyển vị và tiếp cận ứng suất Airy, đảm bảo tính nhất quán nội tại của mô hình lý thuyết.
    • Đối chiếu chuẩn hóa với các công trình chuẩn quốc tế đã công bố cho tấm đẳng hướng và tấm FGM đặc của Reddy [86], Shen [97], và Timoshenko.

Data và phân tích

Toàn bộ thuật toán giải tích được lập trình hóa hoàn chỉnh trên nền tảng phần mềm kỹ thuật MATLAB. Hệ phương trình đại số phi tuyến sau khi chiếu Galerkin được giải bằng thuật toán lặp Newton-Raphson với điều kiện hội tụ sai số chuẩn bậc hai $|R^{(k)}| < 10^{-6}$.

Bộ thông số cơ tính và hình học chuẩn dùng trong mô phỏng số:

  • Vật liệu nền kim loại: Nhôm, Titan hoặc Thép với $E_1 = 200\text{ GPa}$, hệ số Poisson $\nu = 0.3$.
  • Hệ số mật độ lỗ rỗng: $e_0 \in [0.1, 0.9]$.
  • Kích thước tấm: Tấm dày ($a/h = 10$), tấm mỏng ($a/h = 50$), kích thước mặt bằng vuông ($b/a = 1.0$) và chữ nhật ($b/a = 0.5, 1.5, 2.0$).
  • Tham số nền đàn hồi không thứ nguyên: $K_1 = K_w a^4 / D_{11} \in [0, 200]$, $K_2 = K_s a^2 / D_{11} \in [0, 50]$.

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

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

Qua hàng ngàn lượt khảo sát số chi tiết và phân tích cơ chế cơ học, luận án đã rút ra 5 phát hiện mang tính đột phá khoa học:

  1. Quy luật phân bố lỗ rỗng chi phối tuyệt đối độ cứng kết cấu: Tại cùng một hệ số rỗng $e_0 = 0.5$, tấm FGPM có quy luật phân bố lỗ rỗng đối xứng (Dạng 2 - lõi xốp, hai bề mặt đặc) có độ võng phi tuyến nhỏ nhất và tải trọng tới hạn mất ổn định $N_{cr}$ cao nhất. Cụ thể, tải trọng tới hạn của tấm Dạng 2 cao hơn dạng phân bố đều (Dạng 1) từ $28.4%$ đến $35.6%$, và vượt trội hoàn toàn so với dạng phân bố bất đối xứng (Dạng 3) từ $45.2%$ đến $58.1%$. Cơ chế vật lý được giải thích do dạng Dạng 2 tập trung pha kim loại đặc ở hai biên ngoài xa mặt trung hòa nhất, tối đa hóa mô men quán tính uốn và độ cứng chống uốn $C_{11}$.

  2. Ảnh hưởng phi tuyến hình học von Kármán và hiệu ứng làm cứng kết cấu (Stiffening Effect): Trong bài toán uốn chịu tải trọng phân bố đều $P$, phân tích tuyến tính đánh giá độ võng tăng tuyến tính theo tải trọng và cho kết quả độ võng lớn hơn đáng kể so với phân tích phi tuyến. Khi tải trọng $P$ tăng cao, biến dạng kéo màng trong mặt phẳng xuất hiện do chuyển vị uốn lớn, tạo ra hiệu ứng làm cứng hình học (geometric stiffening). Phân tích tuyến tính truyền thống đánh giá thấp khả năng chịu lực thực tế của tấm và làm sai lệch phân bố nội lực mô men ($M_x, M_y, M_{xy}$) lên tới $40%$ ở vùng chuyển vị lớn.

  3. Tác động vượt trội của độ cứng cắt nền Pasternak so với nền Winkler: Khi khảo sát sự tương tác với môi trường đàn hồi, tham số độ cứng cắt $K_s$ (Pasternak shear parameter) thể hiện vai trò kiềm chế biến dạng và gia cường ổn định mạnh gấp nhiều lần so với hệ số uốn Winkler $K_w$. Khi tăng tham số cắt $K_2$ từ $0$ lên $20$, lực tới hạn $N_{cr}$ của tấm SSSS tăng trung bình $62.3%$, trong khi việc tăng tham số $K_1$ từ $0$ lên $100$ chỉ làm tăng $N_{cr}$ khoảng $18.5%$. Điều này chứng minh sự truyền lực cắt liên tục trong môi trường nền có ý nghĩa quyết định đối với khả năng chịu tải của tấm bọt kim loại.

  4. Độ nhạy cảm khuyết tật hình học ban đầu đối với ứng xử sau ổn định: Độ không hoàn hảo ban đầu ($\xi = w_0^*/h$) làm thay đổi căn bản cơ chế làm việc của tấm FGPM chịu nén. Khi $\xi = 0$ (tấm hoàn hảo), tấm duy trì trạng thái phẳng thuần túy cho đến khi đạt lực tới hạn $N_{cr}$ rồi mới xảy ra hiện tượng rẽ nhánh mất ổn định. Tuy nhiên, khi tồn tại độ lệch ban đầu dù rất nhỏ ($\xi = 0.1, 0.2$), tấm bị uốn ngay từ giai đoạn đầu chịu tải, đường cong tải - độ võng phát triển trơn tru, không có điểm rẽ nhánh gián đoạn. Đáng chú ý, ở giai đoạn sau ổn định sâu, đường cong cân bằng của tấm không hoàn hảo tiệm cận dần về đường cong của tấm hoàn hảo nhờ khả năng chịu lực màng dự trữ.

  5. Phân hóa rõ rệt giữa lý thuyết FSDT và CPT theo tỷ số độ dày: Với tấm dày ($a/h = 10$), biến dạng trượt ngang chiếm tỷ phần đáng kể; lý thuyết tấm cổ điển CPT (bỏ qua trượt ngang) đã đánh giá quá cao độ cứng của tấm, dẫn đến dự báo tải trọng tới hạn $N_{cr}$ cao hơn từ $12.5%$ đến $21.8%$ so với lý thuyết FSDT. Ngược lại, khi tấm mỏng dần ($a/h \ge 50$), kết quả tính toán giữa FSDT và CPT tiệm cận trùng khít hoàn toàn (sai số $< 0.8%$), khẳng định phạm vi áp dụng chuẩn xác của từng mô hình lý thuyết.

            SO SÁNH TẢI TRỌNG TỚI HẠN N_cr THEO QUY LUẬT PHÂN BỐ LỖ RỖNG (e0 = 0.5)

Implications đa chiều

  • Về mặt học thuật và phương pháp luận: Cung cấp bộ công thức giải tích chuẩn xác, đóng vai trò là "lời giải tham chiếu" (benchmark analytical solutions) đáng tin cậy để kiểm tra, đánh giá độ chính xác của các phần mềm phần tử hữu hạn thương mại (Abaqus, Ansys) và các phương pháp số mới (IGA, Meshfree, SBFEM).
  • Về thiết kế kỹ thuật công trình và công nghiệp: Cung cấp giải pháp tối ưu hóa hình học kết cấu. Các kỹ sư thiết kế cấu trúc hàng không và vỏ xe bọc thép nên ưu tiên chế tạo bọt kim loại có gradient mật độ dạng Dạng 2 để giảm thiểu tối đa trọng lượng bản thân mà vẫn tối ưu hóa độ cứng chống uốn và ổn định nén.
  • Về hoàn thiện tiêu chuẩn thiết kế: Cung cấp cơ sở dữ liệu khoa học định lượng để xây dựng các hướng dẫn kỹ thuật và tiêu chuẩn tính toán kết cấu composite tiên tiến chịu tải trọng tới hạn trong điều kiện làm việc thực tế.

Limitations và Future Research

Nhằm duy trì tính khách quan khoa học, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu (boundary conditions):

  1. Giả thiết về hệ số Poisson không đổi: Nghiên cứu giả định hệ số Poisson $\nu$ là hằng số dọc theo chiều dày tấm nhằm đơn giản hóa toán học giải tích. Mặc dù ảnh hưởng của biến thiên $\nu$ là thứ yếu so với mô đun Young $E$, trong vật liệu thực tế với độ rỗng cực lớn ($> 80%$), sự thay đổi cấu trúc tế bào có thể làm $\nu$ biến thiên nhẹ.
  2. Mô hình vật liệu đàn hồi tuyến tính: Luận án giả định vật liệu nền kim loại làm việc trong miền đàn hồi tuyến tính, chưa xét đến hiện tượng chảy dẻo cục bộ của thành vách tế bào bọt kim loại dưới ứng suất tập trung cao.
  3. Môi trường cơ học thuần túy: Chưa tích hợp đồng thời trường nhiệt độ biến thiên (nhiệt phi tuyến) và độ ẩm (hygrothermal environment) – các yếu tố thường xuất hiện trong cấu trúc tàu vũ trụ và động cơ phản lực.
  4. Mô hình nền đàn hồi tuyến tính: Nền Pasternak được giả thiết đàn hồi tuyến tính, chưa xét đến tính phi tuyến hình học của môi trường đất nền hoặc biến dạng nhớt theo thời gian.

Từ những hạn chế trên, 5 hướng nghiên cứu tiếp nối giàu tiềm năng được đề xuất:

  1. Mở rộng bài toán sang phân tích động học phi tuyến, dao động tham số và đáp ứng va chạm nổ (blast resistance) của tấm và vỏ cong FGPM.
  2. Thiết lập mô hình cơ học ghép nối nhiệt - ẩm - cơ phi tuyến (Thermo-hygro-mechanical post-buckling).
  3. Nghiên cứu ứng xử từ biến (viscoelastic) và phá hủy dẻo / tổn thương mỏi của vật liệu bọt kim loại tế bào mở.
  4. Tích hợp mô hình vi cơ học đa tỷ lệ (multiscale homogenization) kết hợp quét ảnh cắt lớp vi tính $\mu$-CT cấu trúc lỗ rỗng thực tế vào mô hình giải tích.
  5. Tiến hành các thử nghiệm thực nghiệm uốn lớn và nén sau ổn định trên máy nén thủy lực vạn năng để đối chứng thực nghiệm vật lý với các đường cong giải tích của luận án.

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

Nghiên cứu tạo ra sức ảnh hưởng sâu rộng trên nhiều bình diện:

  • Tác động học thuật quốc tế: Hệ thống công thức giải tích và phương pháp quy chiếu mặt trung hòa mở ra hướng tiếp cận mới mẻ cho cộng đồng cơ học tính toán, ước tính thu hút trích dẫn cao trong các tạp chí ISI/Scopus hàng đầu nhóm Q1 (như Composite Structures, Thin-Walled Structures, Engineering Structures).
  • Chuyển đổi công nghệ công nghiệp: Cung cấp thuật toán tối ưu hóa cấu trúc cho ngành công nghiệp chế tạo phương tiện giao thông vận tải cao tốc, tàu quân sự và hàng không vũ trụ. Việc ứng dụng giải pháp tấm sandwich bọt kim loại giúp giảm trọng lượng kết cấu tới $35%$, trực tiếp cắt giảm tiêu thụ nhiên liệu và giảm lượng khí thải carbon.
  • Nâng cao năng lực tính toán kỹ thuật: Giúp các viện nghiên cứu và doanh nghiệp công nghệ phát triển các module tính toán chuyên dụng gọn nhẹ chạy trên MATLAB hoặc Python mà không cần đầu tư các hệ thống máy chủ đắt tiền để chạy mô phỏng phần tử hữu hạn phi tuyến nặng nề.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận một quy trình chuẩn mực từ mô hình hóa toán học, biến phân năng lượng, giải tích Galerkin đến lập trình số MATLAB giải quyết bài toán phi tuyến phức tạp.
  • Các nhà khoa học cơ học kết cấu: Sử dụng trực tiếp khung lý thuyết mặt trung hòa vật lý để mở rộng cho các dạng kết cấu phức tạp hơn như vỏ nón, vỏ cầu, tấm sandwich nhiều lớp gia cường nano carbon (CNT/Graphene).
  • Kỹ sư R&D công nghiệp: Vận dụng các biểu đồ phi thứ nguyên và công thức tường minh trong luận án để tra cứu nhanh tải trọng tới hạn, tối ưu hóa kích thước tấm panel sàn, tường giảm chấn và vỏ chịu lực.
  • Cơ quan quản lý và biên soạn tiêu chuẩn: Tham khảo các kết luận định lượng về hệ số an toàn và độ nhạy khuyết tật hình học ban đầu để xây dựng quy chuẩn an toàn kết cấu vật liệu 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ì 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 vi phân phi tuyến chủ đạo cho tấm FGPM trên hệ tọa độ gắn với mặt trung hòa vật lý ($z_{th}$), chứng minh giải tích sự triệt tiêu hoàn toàn của ma trận độ cứng liên kết $B_{ij} = 0$. Công trình đã mở rộng trực tiếp Lý thuyết biến dạng cắt bậc nhất của Mindlin (1951) và Lý thuyết tấm cổ điển Kirchhoff-Love sang môi trường vật liệu bọt kim loại có độ xốp biến thiên liên tục theo chiều dày, giải quyết triệt để sự phức tạp trong ghép nối màng - uốn mà các nghiên cứu quy chiếu mặt trung bình trước đây gặp phải.

2. Đột phá về phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế?

So với nghiên cứu của Magnucki và cs. [71] (chỉ giải tuyến tính cho tấm hoàn hảo) và Shen & Wang [98] (sử dụng phương pháp hàm phạt hai bước trên mặt trung bình), luận án đã:

  • Tích hợp song song cả hai tiếp cận: Tiếp cận chuyển vị (chuỗi lượng giác Navier kép giải bằng Bubnov-Galerkin) và Tiếp cận ứng suất (sử dụng hàm ứng suất Airy $F(x,y)$).
  • Rút ra biểu thức giải tích tường minh cho tải trọng tới hạn $N_{cr}$ và thiết lập toàn diện đường cong cân bằng sau ổn định của tấm FGPM có kể đến đồng thời độ lệch ban đầu $w_0^*$, hiệu ứng nền Pasternak ($K_w, K_s$) và điều kiện biên phức tạp (SSSS, CCCC, SCSC).

3. Phát hiện bất ngờ và phản trực giác nhất từ số liệu là gì?

Phát hiện nổi bật là sự suy giảm tải trọng tới hạn phi tuyến của quy luật phân bố lỗ rỗng bất đối xứng (Dạng 3). Mặc dù cùng có hệ số mật độ lỗ rỗng trung bình $e_0$ và khối lượng tấm hoàn toàn tương đương với Dạng 2, dạng phân bố bất đối xứng Dạng 3 làm mất ổn định kết cấu ở mức tải trọng chỉ bằng khoảng $55%$ so với Dạng 2. Nguyên nhân là do sự lệch tâm của mặt trung hòa làm suy giảm mô men kháng uốn hiệu dụng tại thớ chịu nén cực đại, kích hoạt hiện tượng mất ổn định uốn - nén sớm.

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 minh bạch và đầy đủ. Toàn bộ các bước biến đổi toán học từ trường chuyển vị, biến dạng von Kármán, tích phân nội lực ($A_{ij}, C_{ij}, A_{44}^s$), biểu thức thế năng toàn phần, các hệ số trong ma trận Galerkin đến thuật toán giải phương trình phi tuyến trên MATLAB đều được trình bày chi tiết trong Chương 2, Chương 3 và Chương 4 của luận án. Các bảng số liệu không thứ nguyên chuẩn mực được cung cấp sẵn sàng cho việc tái lập và kiểm chứng độc lập.

5. Lộ trình phát triển nghiên cứu 10 năm tới được phác thảo ra sao?

Lộ trình 10 năm định hình sự phát triển từ giải tích sang công nghệ ứng dụng:

  • Giai đoạn 1-3 năm: Tích hợp trường đa vật lý (Cơ - Nhiệt - Điện - Từ) và động lực học phi tuyến va chạm tốc độ cao.
  • Giai đoạn 4-6 năm: Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) và tối ưu hóa tô pô đa mục tiêu để thiết kế cấu trúc vi mô lỗ rỗng FGPM tối ưu.
  • Giai đoạn 7-10 năm: Chuyển giao công nghệ in 3D kim loại (Additive Manufacturing) chế tạo thực nghiệm tấm FGPM và thương mại hóa các module tính toán kết cấu trong kỹ thuật xây dựng và công nghiệp quốc phòng.

Kết luận

Luận án tiến sĩ của tác giả Lê Thanh Hải là một công trình nghiên cứu khoa học công phu, nghiêm túc và có giá trị học thuật xuất sắc, đóng góp những viên gạch vững chắc cho ngành Cơ kỹ thuật Việt Nam và thế giới:

  1. Xây dựng hoàn chỉnh khung lý thuyết giải tích: Thiết lập hệ thức cơ bản và hệ phương trình chủ đạo phân tích phi tuyến uốn, ổn định và sau ổn định của tấm FGPM trên nền đàn hồi Pasternak theo cả hai mô hình FSDT và CPT.
  2. Giải trừ triệt để liên kết ghép màng - uốn: Chứng minh và ứng dụng thành công hệ quy chiếu mặt trung hòa vật lý ($B_{ij} = 0$), giúp đơn giản hóa cấu trúc toán học và tiết kiệm đáng kể chi phí tính toán số.
  3. Phát triển lời giải giải tích kép: Thiết lập thành công lời giải theo tiếp cận chuyển vị và tiếp cận hàm ứng suất Airy kết hợp thuật toán trực giao hóa Bubnov-Galerkin, rút ra công thức hiển của tải tới hạn $N_{cr}$ và quan hệ tải - độ võng phi tuyến.
  4. Làm sáng tỏ cơ chế tương tác đa tham số: Đánh giá định lượng toàn diện ảnh hưởng của 3 quy luật phân bố lỗ rỗng, hệ số rỗng $e_0$, tỷ số kích thước $a/h, b/a$, hệ số nền Pasternak $K_w, K_s$, điều kiện biên SSSS/CCCC/SCSC và độ không hoàn hảo hình học $w_0^*$.
  5. Định hướng ứng dụng kỹ thuật thực tiễn: Chứng minh ưu thế vượt trội của cấu trúc bọt kim loại phân bố đối xứng (Dạng 2), cung cấp nguồn dữ liệu kỹ thuật tin cậy phục vụ thiết kế, chế tạo các kết cấu nhẹ chịu lực cao trong ngành hàng không vũ trụ, giao thông vận tải và xây dựng công trình hiện đại.