Tổng quan về luận án

Nghiên cứu về ứng xử cơ học của kết cấu tấm làm từ vật liệu phân lớp chức năng (Functionally Graded Materials - FGM) chịu tải trọng cơ nhiệt là một trong những chủ đề tiên phong của cơ học vật rắn biến dạng và kỹ thuật kết cấu hiện đại. Khái niệm FGM lần đầu tiên được khởi xướng vào giữa thập niên 1980 tại Nhật Bản bởi một nhóm các nhà khoa học vật liệu nhằm giải quyết bài toán lá chắn nhiệt trong ngành hàng không vũ trụ, nơi kết cấu phải làm việc trong điều kiện chênh lệch nhiệt độ khắc nghiệt (mặt ngoài chịu nhiệt độ gốm lên đến $2000\text{ K}$, mặt trong chịu nhiệt độ kim loại khoảng $1000\text{ K}$). Bằng cách thay đổi liên tục thành phần vật liệu theo quy luật gradient không gian, FGM loại bỏ triệt để hiện tượng tập trung ứng suất và suy thoái liên kết tại bề mặt tiếp xúc vốn là nhược điểm chí mạng của vật liệu composite phân lớp truyền thống.

Tuy nhiên, việc mô hình hóa chính xác ứng xử cơ học của tấm FGM đối mặt với các khoảng trống nghiên cứu (research gaps) trọng yếu: (1) Sự không đồng nhất theo chiều dày dẫn đến sự lệch pha giữa mặt trung bình hình học và mặt trung hòa vật lý, gây ra tương tác phức tạp giữa kéo/nén và uốn; (2) Hệ số hiệu chỉnh biến dạng cắt $\kappa_s$ trong lý thuyết biến dạng cắt bậc nhất (First-order Shear Deformation Theory - FSDT) thường bị áp đặt cố định là $5/6$, không phản ánh đúng sự biến thiên liên tục của mô-đun cắt $G(z)$; (3) Các lý thuyết biến dạng cắt bậc cao (Higher-order Shear Deformation Theory - HSDT) và tiếp cận 3 chiều (Quasi-3D) đòi hỏi hàm dạng phức tạp và gặp khó khăn khi xây dựng các mô hình số có tính tương thích $C^1$; (4) Hiện tượng "khóa cắt" (shear locking) làm suy giảm nghiêm trọng độ chính xác của phương pháp phần tử hữu hạn (Finite Element Method - FEM) khi phân tích tấm mỏng; (5) Sự thiếu vắng các đánh giá so sánh hệ thống về ảnh hưởng của các mô hình đồng nhất hóa vi mô (homogenization schemes) đến ứng xử vĩ mô của tấm FGM và tấm sandwich FGM.

Để giải quyết triệt để các khoảng trống trên, luận án xác định hệ thống câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H):

  • RQ1: Làm thế nào để loại bỏ tương tác kéo - uốn nhân tạo và chuẩn hóa hệ số hiệu chỉnh cắt trong lý thuyết FSDT cho tấm FGM?
    H1: Việc tích hợp vị trí mặt trung hòa vật lý $z_0$ kết hợp với độ cứng cắt cải tiến $H$ suy dẫn từ cân bằng năng lượng biến dạng cắt sẽ nâng cao độ chính xác của lời giải giải tích và số mà không làm tăng bậc tự do của bài toán.
  • RQ2: Hàm biến dạng cắt bậc cao mới nào có khả năng mô tả chính xác trường ứng suất cắt ngang mà không cần hệ số hiệu chỉnh cắt, đồng thời mở rộng được cho lý thuyết Quasi-3D có xét biến dạng dọc trục $z$?
    H2: Một hàm biến dạng cắt bậc cao mới kết hợp số hạng phi tuyến của tọa độ chiều dày cho phép thỏa mãn hoàn toàn điều kiện biên ứng suất cắt bằng 0 tại hai bề mặt tự do, phản ánh đúng hiệu ứng biến dạng theo chiều dày của tấm dày.
  • RQ3: Giải pháp nào triệt tiêu hoàn toàn hiện tượng khóa cắt trong phân tích phần tử hữu hạn cho cả bài toán tuyến tính và phi tuyến hình học Von Kármán?
    H3: Kỹ thuật làm trơn biến dạng trên miền (CS), trên cạnh (ES) và trên nút (NS) kết hợp với phần tử nội suy ten-sơ hỗn hợp 3 nút (MITC3) sẽ loại bỏ hiện tượng khóa cắt, tăng tốc độ hội tụ và cho kết quả ổn định ngay cả với lưới phần tử tam giác bất kỳ.
  • RQ4: Mức độ phân tán của các mô hình đồng nhất hóa vi mô tác động như thế nào đến độ võng, lực tới hạn và tần số dao động riêng của kết cấu?
    H4: Sự sai lệch cơ tính giữa các mô hình vi cơ học (Voigt, Reuss, Hashin-Shtrikman, Tamura, SCM, LRVE) có thể tạo ra chênh lệch ứng xử kết cấu lên tới xấp xỉ $48%$, đòi hỏi các giới hạn biên rõ ràng khi lựa chọn mô hình tính toán.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng cơ học môi trường liên tục, nguyên lý biến phân Hamilton, phương trình năng lượng Lagrange, kết hợp lý thuyết nền đàn hồi hai tham số Winkler-Pasternak ($K_w, K_s$). Luận án có dung lượng học thuật đồ sộ gồm 171 trang (chưa kể tài liệu tham khảo và phụ lục), 100 hình vẽ minh họa chi tiết và 60 bảng số liệu thực nghiệm số, được bố cục chặt chẽ thành 7 chương nghiên cứu chuyên sâu, khảo sát toàn diện 4 cấu hình kết cấu tấm: Tấm FGM đơn lớp (Loại A), Tấm sandwich FGM lõi FGM lớp mặt đồng chất (Loại B), Tấm sandwich FGM lõi đồng chất lớp mặt FGM (Loại C) và Tấm composite nhiều lớp có góc nghiêng sợi tùy ý (Loại D).

Literature Review và Positioning

Tổng quan y văn quốc tế cho thấy sự tiến hóa liên tục của các trường phái lý thuyết tấm trong cơ học kết cấu:

  1. Trường phái lý thuyết tấm cổ điển (CPT): Khởi xướng từ các giả thiết Love-Kirchhoff (Timoshenko & Woinowsky-Krieger [20], Reddy [21]), bỏ qua hoàn toàn biến dạng cắt ngang ($\gamma_{xz} = \gamma_{yz} = 0$). Dù được He và cộng sự [24], Chi & Chung [25, 26], Srinivas & Rao [27] áp dụng ban đầu cho FGM, CPT bộc lộ sai số lớn khi áp dụng cho tấm có chiều dày từ trung bình đến lớn ($a/h \le 20$).
  2. Trường phái lý thuyết biến dạng cắt bậc nhất (FSDT): Dựa trên mô hình Reissner-Mindlin, được Praveen & Reddy [28], Batra & Jin [31], Li và cộng sự [30] ứng dụng rộng rãi. Tuy nhiên, nhược điểm cốt lõi là ứng suất cắt ngang giả định phân bố đều qua chiều dày, buộc phải sử dụng hệ số hiệu chỉnh cắt $\kappa_s$. Các nghiên cứu của Nguyen và cộng sự [32, 33, 34, 35] đã nỗ lực đề xuất hệ số hiệu chỉnh nhưng chưa liên kết chặt chẽ với sự dịch chuyển của mặt trung hòa vật lý. Việc xác định mặt trung hòa vật lý $z_0$ được tiếp cận bởi Singha và cộng sự [36, 40], Ma & Lee [37], Zhang và cộng sự [38, 39], Lee và cộng sự [46], chứng minh rằng việc tính từ mặt trung bình hình học sẽ dẫn đến các số hạng tương tác kéo - uốn giả tạo trong ma trận độ cứng $B_{ij}$.
  3. Trường phái lý thuyết biến dạng cắt bậc cao (HSDT) và Quasi-3D: Để loại bỏ hệ số hiệu chỉnh cắt, Reddy [47, 48, 49] phát triển lý thuyết biến dạng cắt bậc 3 (TSDT). Tiếp đó, hàng loạt hàm biến dạng cắt được đề xuất: dạng hyperbolic bởi Soldatos [61], Akavci [63]; dạng hàm sin bởi Touratier [64], Thai & Vo [65]; dạng hàm mũ bởi Karama [66]; dạng lượng giác kết hợp hàm mũ bởi Mantari và cộng sự [67, 68]; đa thức bậc $n$ bởi Xiang & Kang [69]; hàm lượng giác ngược bởi Thai và cộng sự [70] và hyperbolic ngược bởi Grover và cộng sự [71]. Đối với lý thuyết Quasi-3D, Carrera và cộng sự [73], Neves và cộng sự [74, 75], Thai và cộng sự [76, 77, 78], Mantari & Soares [79, 80] đã tích hợp thành phần biến dạng pháp $\varepsilon_z$, phản ánh trạng thái ứng suất 3 chiều thực tế của tấm dày.
  4. Trường phái phương pháp số và xử lý khóa cắt: Trong tính toán phần tử hữu hạn, hiện tượng khóa cắt khi $a/h \to \infty$ được giải quyết qua các kỹ thuật: giả định biến dạng tự nhiên (ANS) [95, 96], biến dạng tự nhiên nâng cao (EAS) [98, 99], rời rạc sự khác biệt cắt (DSG) [100, 101], và nội suy ten-sơ hỗn hợp (MITC) [102-108]. Đồng thời, các kỹ thuật làm trơn biến dạng (S-FEM) trên miền (CS), trên cạnh (ES), trên nút (NS) [109-115] và phương pháp đẳng hình học IGA (Hughes và cộng sự [124], [125-128]) đã mở ra hướng tiếp cận mới vượt trội so với phương pháp không lưới [116-122].

Y văn tồn tại hai cuộc tranh luận học thuật lớn:

  • Tranh luận 1 (Hệ số hiệu chỉnh cắt cố định vs. Động học biến thiên): Một trường phái bảo thủ giữ nguyên giá trị $\kappa_s = 5/6$ từ tấm đồng chất, trong khi trường phái cải tiến chỉ ra rằng trong tấm sandwich FGM phân lớp bất đối xứng, $\kappa_s$ phụ thuộc phi tuyến vào tỷ lệ thể tích $V_c(z)$ và tỷ số mô-đun đàn hồi $E_c/E_m$.
  • Tranh luận 2 (Mô hình vi cơ học đơn giản hóa vs. Mô hình tự nhất quán): Mô hình Voigt tuyến tính (quy tắc trộn thể tích) thường đánh giá quá cao độ cứng, trong khi mô hình Reuss (nghịch đảo) đánh giá quá thấp. Tranh cãi nằm ở việc liệu các mô hình phức tạp như Hashin-Shtrikman [6, 7], Tamura [8, 9], Self-Consistent Method (SCM) [10, 11] hay LRVE (Gasik & Lilius [12, 13]) có thực sự mang lại độ chính xác vượt trội khi đưa vào mô hình phân tích kết cấu vĩ mô hay không.

Luận án định vị nghiên cứu tại điểm giao thoa giữa cơ học giải tích chính xác và phương pháp số bậc cao. Khi so sánh với các công trình quốc tế tiêu biểu như nghiên cứu của Mohammadi và cộng sự [131] về bài toán ổn định và Zenkour [50, 51] về bài toán uốn tĩnh, mô hình của luận án thể hiện sự vượt trội nhờ việc kết hợp đồng thời mặt trung hòa vật lý, hệ số hiệu chỉnh cắt cải tiến phụ thuộc năng lượng và phần tử hữu hạn trơn làm triệt tiêu hoàn toàn khóa cắt.

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

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

Luận án mang lại các đột phá lý thuyết quan trọng trong cơ học kết cấu FGM:

  1. Chuẩn hóa lý thuyết FSDT thông qua mặt trung hòa vật lý ($z_0$): Bằng việc xác định chính xác vị trí mặt trung hòa vật lý theo biểu thức giải tích: $$z_0 = \frac{\int_{-h/2}^{h/2} z E(z) dz}{\int_{-h/2}^{h/2} E(z) dz}$$ toàn bộ các thành phần trong ma trận tương tác kéo - uốn $B_{ij}$ tự động triệt tiêu ($B_{ij} = 0$). Điều này tách rời hoàn toàn phương trình vi phân chủ đạo của chuyển vị trong mặt phẳng ($u, v$) và chuyển vị ngoài mặt phẳng ($w$), đơn giản hóa cấu trúc toán học của hệ phương trình Euler-Lagrange.

  2. Thiết lập biểu thức giải tích cho độ cứng cắt cải tiến ($H$): Thay vì áp đặt hệ số $\kappa_s = 5/6$, luận án xây dựng biểu thức độ cứng cắt $H$ dựa trên nguyên lý cân bằng năng lượng biến dạng cắt trên một đơn vị diện tích: $$\frac{1}{H} = \int_{-h/2}^{h/2} \frac{[R(z)]^2}{G(z)} dz \quad \text{với} \quad R(z) = \int_{-h/2}^{z} C_{11}(\zeta)(b + (\zeta - z_0)d) d\zeta$$ Từ đó, hệ số hiệu chỉnh cắt cải tiến được định nghĩa chính xác theo cấu trúc vật liệu: $$\kappa_s = \frac{H}{\int_{-h/2}^{h/2} C_{55}(z) dz}$$

  3. Đề xuất hàm biến dạng cắt bậc cao mới và mô hình Quasi-3D tổng quát: Luận án xây dựng hàm biến dạng cắt $f(z)$ mới khắc phục nhược điểm của các hàm dạng truyền thống, thỏa mãn điều kiện biên ứng suất tiếp tại hai mặt tự do $\sigma_{xz}(\pm h/2) = \sigma_{yz}(\pm h/2) = 0$. Mô hình Quasi-3D mở rộng trường chuyển vị: $$u_3(x,y,z) = w(x,y) + g(z)\phi_z(x,y) \quad \text{với} \quad g(z) = \frac{df(z)}{dz}$$ cho phép nắm bắt trực tiếp biến dạng co giãn dọc trục $\varepsilon_z$, tạo bước chuyển dịch mô hình (paradigm shift) từ phân tích 2D sang tiệm cận chính xác 3D elasticity.

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

  • Mệnh đề 1: Tấm FGM bất đối xứng quy chiếu về mặt trung hòa vật lý sẽ triệt tiêu hoàn toàn hiệu ứng ghép nối màng - uốn ($B_{ij}=0$), làm cho phản ứng uốn thuần túy không bị nhiễu bởi trường lực dọc trục.
  • Mệnh đề 2: Độ cứng cắt thực tế của tấm sandwich FGM tỷ lệ nghịch với tích phân bình phương hàm ứng suất cắt cảm ứng $R(z)$ chia cho mô-đun cắt cục bộ $G(z)$.
  • Mệnh đề 3: Trường chuyển vị 4 biến tách thành phần uốn ($w_b$) và cắt ($w_s$) bảo toàn tính chính xác của lý thuyết biến dạng cắt bậc cao trong khi giảm số bậc tự do của phần tử từ 5 xuống 4.
  • Mệnh đề 4: Năng lượng biến dạng cắt trong miền làm trơn hội tụ đơn điệu về nghiệm giải tích khi kích thước phần tử tiến tới 0 mà không phát sinh hiện tượng khóa cắt thể tích hoặc khóa cắt mỏng.

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

Khung phân tích của luận án là sự tích hợp đa chiều của 3 lý thuyết cơ học nền tảng: Lý thuyết đàn hồi dị hướng biến thiên liên tục, Lý thuyết biến dạng cắt bậc cao Quasi-3D và Lý thuyết phi tuyến hình học Von Kármán biến dạng nhỏ chuyển vị lớn: $$\varepsilon_{xx} = \frac{\partial u}{\partial x} + \frac{1}{2}\left(\frac{\partial w}{\partial x}\right)^2, \quad \varepsilon_{yy} = \frac{\partial v}{\partial y} + \frac{1}{2}\left(\frac{\partial w}{\partial y}\right)^2, \quad \gamma_{xy} = \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x} + \frac{\partial w}{\partial x}\frac{\partial w}{\partial y}$$

Khung phân tích bao quát toàn diện các điều kiện biên và tương tác nền:

  • Tương tác môi trường nền đàn hồi: Tích hợp thế năng biến dạng nền Winkler ($K_w$) và Pasternak ($K_s$): $$U_F = \frac{1}{2} \int_A \left[ K_w w^2 + K_s \left( \left(\frac{\partial w}{\partial x}\right)^2 + \left(\frac{\partial w}{\partial y}\right)^2 \right) \right] dA$$
  • Không gian biên mở rộng: Phân tích từ liên kết tựa đơn 4 cạnh ($SSSS$), ngàm 4 cạnh ($CCCC$) đến các tổ hợp tựa đơn - ngàm tự do phức tạp ($SCSC, SCSF, CFCF$).

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 (positivist paradigm) và cơ sở nhận thức luận suy diễn chặt chẽ (deductive epistemology) của cơ học giải tích và tính toán số. Thiết kế nghiên cứu phân cấp 3 tầng rõ rệt:

  • Tầng 1 (Cơ học vi mô): Xác định các đặc trưng đàn hồi hữu hiệu $E(z), \rho(z), \nu(z)$ thông qua 6 mô hình đồng nhất hóa: mô hình Voigt, mô hình Reuss, mô hình Hashin-Shtrikman cận trên/cận dưới, mô hình Tamura (với thông số truyền ứng suất - biến dạng $q_T$), mô hình Self-Consistent (SCM) và mô hình LRVE của Gasik & Lilius.
  • Tầng 2 (Cơ học kết cấu biến phân): Thiết lập phương trình cân bằng động học thông qua nguyên lý biến phân Hamilton: $$\int_0^T (\delta U + \delta V - \delta K + \delta U_F) dt = 0$$
  • Tầng 3 (Giải thuật số và giải tích): Triển khai song song lời giải giải tích kép (Chuỗi Navier cho biên $SSSS$ và Chuỗi Ritz trực giao cho biên bất kỳ) đối chiếu trực tiếp với Phương pháp phần tử hữu hạn làm trơn cải tiến.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu số được thiết kế kiểm soát sai số nghiêm ngặt qua kỹ thuật khử khóa cắt và làm trơn biến dạng:

  1. Phần tử tấm MITC3: Sử dụng phần tử tam giác 3 nút với các hàm nội suy Lagrange tiêu chuẩn cho chuyển vị, kết hợp nội suy trực tiếp ten-sơ biến dạng cắt trên các cạnh phần tử để loại bỏ hoàn toàn khóa cắt mỏng: $$\gamma_{xz} = \gamma_{xz}^{MITC}, \quad \gamma_{yz} = \gamma_{yz}^{MITC}$$
  2. Kỹ thuật làm trơn biến dạng (S-FEM): Miền phần tử $\Omega_e$ được chia thành các miền con làm trơn $\Omega_k$. Biến dạng trơn $\tilde{\boldsymbol{\varepsilon}}k$ được tính bằng tích phân đường trên biên $\Gamma_k$: $$\tilde{\boldsymbol{\varepsilon}}k = \frac{1}{A_k} \int{\Omega_k} \boldsymbol{\varepsilon} d\Omega = \frac{1}{A_k} \int{\Gamma_k} \mathbf{n} \mathbf{u} d\Gamma$$ Luận án khảo sát đồng thời 3 mô hình làm trơn:
  • CS-MITC3 (Cell-based Smoothed): Làm trơn trên các miền con bên trong phần tử.
  • ES-MITC3 (Edge-based Smoothed): Làm trơn trên miền liên kết bởi cạnh chung của hai phần tử liền kề.
  • NS-MITC3 (Node-based Smoothed): Làm trơn trên miền bao quanh nút phần tử.
  1. Quy tắc tích phân số Gauss: Đối với phần tử chữ nhật 4 nút cấp thấp, áp dụng kỹ thuật tích phân giảm (reduced integration): ma trận độ cứng cắt $H_{ij}$ được tính bằng 1 điểm Gauss ($\xi=0, \eta=0, w_i=w_j=2$), trong khi các ma trận độ cứng màng và uốn được tích hợp bằng lưới $2 \times 2$ điểm Gauss ($\xi, \eta = \pm 1/\sqrt{3}, w_i=w_j=1$).

Data và phân tích

Toàn bộ hệ thống phương trình và thuật toán được lập trình nguyên bản (in-house scripts) trên nền tảng ngôn ngữ kỹ thuật MATLAB, giải quyết 3 bài toán vật lý trọng tâm:

  • Bài toán tĩnh: Giải hệ phương trình đại số tuyến tính $\mathbf{K} \mathbf{U} = \mathbf{F}$ (hoặc bài toán phi tuyến $\mathbf{K}(\mathbf{U})\mathbf{U} = \mathbf{F}$ bằng giải thuật lặp Newton-Raphson).
  • Bài toán ổn định: Xác định tải trọng tới hạn $N_{cr}$ qua phương trình đặc trưng $\det(\mathbf{K} - N_{cr} \mathbf{K}_g) = 0$.
  • Bài toán dao động tự do: Giải bài toán trị riêng tổng quát $\det(\mathbf{K} - \omega^2 \mathbf{M}) = 0$ để tìm tần số góc $\omega$ và dạng dao động riêng.

Các hệ số vật liệu chuẩn được sử dụng xuyên suốt:

  • Pha kim loại (Aluminum Alloy 1100): $E_m = 70\text{ GPa}$, $\rho_m = 2707\text{ kg/m}^3$, $\nu_m = 0.3$.
  • Pha gốm (Alumina $\text{Al}_2\text{O}_3$): $E_c = 380\text{ GPa}$, $\rho_c = 3800\text{ kg/m}^3$, $\nu_c = 0.3$.
  • Các đại lượng không thứ nguyên chuẩn hóa: $$\bar{w} = \frac{100 h^3 E_c}{a^4 q_0} w\left(\frac{a}{2}, \frac{b}{2}\right), \quad \bar{\sigma}{xx} = \frac{h^2}{a^2 q_0} \sigma{xx}\left(\frac{a}{2}, \frac{b}{2}, \frac{h}{3}\right), \quad \bar{\omega} = \omega h \sqrt{\frac{\rho_0}{E_0}}, \quad \bar{N}{cr} = \frac{N{cr} a^2}{100 E_0 h^3}$$

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

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

  1. Hiệu ứng phân tán của các mô hình đồng nhất hóa vi mô: Kết quả tính toán chỉ ra rằng mô hình Voigt cho giá trị mô-đun Young $E(z)$ lớn nhất, trong khi mô hình Reuss cho giá trị bé nhất. Tại vị trí tỷ lệ thể tích $V(z) = 0.5$, mức độ chênh lệch mô-đun đàn hồi giữa Voigt và Reuss lên tới $48%$. Tương tự, khoảng cách giữa xấp xỉ cận trên và cận dưới của Hashin-Shtrikman là $22%$. Các mô hình Tamura ($q_T = 100\text{ GPa}$), LRVE và SCM cho kết quả rất sát nhau và nằm ở khoảng trung gian tin cậy.

  2. Quy luật biến thiên độ võng theo chỉ số gradient ($p$): Độ võng không thứ nguyên $\bar{u}_3$ tăng đơn điệu khi hệ số đặc trưng vật liệu $p$ tăng từ 0 (gốm thuần túy) đến vô cùng (kim loại thuần túy). Cụ thể, với tấm vuông FGM tựa đơn ($a/h = 10$):

    • Khi $p = 0$ (Gốm $\text{Al}_2\text{O}_3$ thuần): $\bar{u}3 = 0.2960$, ứng suất tiếp cực đại $\bar{\sigma}{xz}^{max} = 0.6226$.
    • Khi $p = 1$: $\bar{u}_3 = 0.5858$ (tăng gần gấp đôi do sự suy giảm của hàm lượng gốm làm mềm kết cấu).
    • Khi $p \to \infty$ (Nhôm thuần): $\bar{u}3 = 1.6083$, $\bar{\sigma}{xz}^{max} = 2.6666$.

  3. Tính chính xác vượt trội của mô hình độ cứng cắt cải tiến ($H$): So sánh với hệ số cố định $\kappa_s = 5/6$, mô hình độ cứng cắt $H$ của luận án phản ánh chính xác sự làm việc của tấm sandwich FGM. Đối với kết cấu tấm sandwich Loại B (lõi cứng FGM), mô hình $\kappa_s = 5/6$ đánh giá thấp độ cứng cắt thực tế, dẫn đến tính toán độ võng lớn hơn thực tế. Ngược lại, đối với tấm sandwich Loại C (lớp mặt FGM, lõi mềm), mô hình $\kappa_s = 5/6$ đánh giá quá cao độ cứng cắt, dẫn đến kết quả độ võng tính toán nhỏ hơn thực tế, tiềm ẩn nguy cơ mất an toàn kết cấu.

  4. Hiệu năng vượt bậc của phần tử hữu hạn trơn ES-MITC3 và CS-MITC3: Khi phân tích bài toán tĩnh và dao động tự do, phần tử ES-MITC3 và CS-MITC3 thể hiện tốc độ hội tụ vượt trội, tiệm cận hoàn toàn với lời giải chính xác 3D Elasticity và triệt tiêu hoàn toàn hiện tượng khóa cắt ngay cả khi tỷ số cạnh trên chiều dày $a/h = 1000$.

  5. Đáp ứng ổn định phi tuyến dưới tải trọng kết hợp: Dưới các trạng thái ứng suất màng khác nhau (nén 1 phương $R_1=1, R_2=0$; nén 2 phương $R_1=1, R_2=1$; nén - kéo kết hợp $R_1=1, R_2=-1$), tải trọng tới hạn $N_{cr}$ đạt giá trị lớn nhất trong trạng thái nén - kéo và giảm mạnh nhất trong trạng thái nén 2 phương, với quy luật phi tuyến hình học Von Kármán làm tăng đáng kể độ cứng sau ổn định (post-buckling stiffness).

Implications đa chiều

  • Về mặt lý thuyết: Đóng góp hệ phương trình giải tích tường minh cho tấm FGM và sandwich FGM quy chiếu về mặt trung hòa vật lý, giải quyết dứt điểm các tranh cãi kéo dài về hệ số hiệu chỉnh cắt trong lý thuyết Reissner-Mindlin.
  • Về mặt phương pháp luận: Thiết lập chuẩn mực mô phỏng số mới với các phần tử tam giác trơn CS/ES-MITC3, cho phép chia lưới tự động trên các miền hình học phức tạp mà không lo ngại suy biến phần tử hay khóa cắt.
  • Về ứng dụng thực tiễn: Cung cấp công cụ tính toán chính xác tuyệt đối cho các kỹ sư thiết kế tấm chắn nhiệt buồng đốt động cơ tên lửa, vỏ chịu áp lực của tàu ngầm hạt nhân, và kết cấu composite nano chịu nhiệt độ cao trong ngành công nghiệp quốc phòng và hàng không vũ trụ.

Limitations và Future Research

Dù đạt được những kết quả đột phá, nghiên cứu tồn tại một số giới hạn học thuật cần được mở rộng:

  1. Giới hạn giả thiết đàn hồi tuyến tính của vật liệu: Nghiên cứu giới hạn ở bài toán đàn hồi thuần túy, chưa xét đến biến dạng dẻo vi mô (micro-plasticity), từ biến (creep) và hiện tượng rạn nứt tế vi tại các pha hạt gốm dưới sốc nhiệt chu kỳ lớn.
  2. Giới hạn hình học kết cấu: Các giải pháp giải tích tập trung vào kết cấu tấm phẳng hình chữ nhật; chưa mở rộng cho kết cấu vỏ cong bậc hai phức tạp (double-curved shells) hoặc kết cấu có lỗ khoét hình học bất kỳ.
  3. Giới hạn liên kết pha: Giả định liên kết hoàn hảo (perfect bonding) giữa các vi hạt và giữa các phân lớp, chưa xét đến hiện tượng trượt tiếp xúc (interfacial debonding) hoặc hư hỏng lớp kết dính.

Chương trình nghiên cứu tiếp nối trong giai đoạn 5-10 năm tới bao gồm:

  • Mở rộng phương pháp đẳng hình học (IGA) kết hợp hàm biến dạng Quasi-3D mới để phân tích kết cấu vỏ FGM không gian phức tạp.
  • Xây dựng mô hình nhiệt - cơ - điện - từ tương tác (Thermo-electro-magneto-mechanical coupling) cho các tấm composite thông minh FGM tích hợp lớp áp điện (piezoelectric layers).
  • Tích hợp mô hình cơ học vi nứt và cơ học phá hủy liên tục (Damage Mechanics) để mô phỏng quá trình suy thoái độ bền mỏi dưới tải trọng nhiệt động lực học.
  • Ứng dụng trí tuệ nhân tạo và giải thuật tối ưu hóa đa mục tiêu (Machine Learning & Genetic Algorithms) để tối ưu hóa hàm phân bố gradient $V_c(z)$ nhằm tối thiểu hóa khối lượng và tối đa hóa tải trọng tới hạn.

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

Luận án tạo ra những tác động khoa học và thực tiễn sâu rộng:

  • Tác động học thuật: Định hình lại phương pháp phân tích tấm FGM trong các tạp chí cơ học hàng đầu quốc tế (như Composite Structures, Computer Methods in Applied Mechanics and Engineering, International Journal of Mechanical Sciences), với tiềm năng trích dẫn ước tính hàng trăm lượt từ các nhóm nghiên cứu cơ học tính toán trên thế giới.
  • Chuyển đổi công nghiệp: Chuyển giao các thuật toán phần tử hữu hạn làm trơn CS/ES-MITC3 vào các phần mềm CAE thương mại và mã nguồn mở, nâng cao năng lực tính toán thiết kế cho ngành cơ khí chính xác, công nghiệp ô tô điện và hàng không vũ trụ.
  • Ảnh hưởng chính sách và tiêu chuẩn: Cung cấp căn cứ khoa học chuẩn xác phục vụ việc xây dựng các bộ tiêu chuẩn tính toán kết cấu vật liệu mới và vật liệu tiên tiến chịu nhiệt tại Việt Nam và khu vực.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một tài liệu tham khảo mẫu mực về quy trình thiết lập phương trình biến phân Hamilton, kỹ thuật khử khóa cắt và phương pháp lập trình số MATLAB hoàn chỉnh từ lý thuyết đến code thực thi.
  • Kỹ sư R&D Hàng không - Vũ trụ và Giao thông: Sở hữu công cụ tính toán chính xác độ võng, ứng suất tập trung và tần số dao động để tối ưu hóa thiết kế tấm cách nhiệt nhẹ, bền, chịu sốc nhiệt cao.
  • Nhà phát triển phần mềm CAE: Ứng dụng trực tiếp cấu trúc ma trận độ cứng phần tử làm trơn CS-MITC3 và ES-MITC3 để nâng cấp lõi tính toán (solver core) cho các bài toán phân tích phi tuyến hình học tấm và vỏ.

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ì?

Đóng góp độc đáo nhất là việc xây dựng thành công mô hình giải tích FSDT tích hợp đồng thời mặt trung hòa vật lý ($z_0$) và độ cứng cắt cải tiến ($H$). Bằng cách định nghĩa lại hệ quy chiếu chuyển vị theo $z_0$, ma trận độ cứng liên kết kéo - uốn $B_{ij}$ tự động triệt tiêu, biến hệ phương trình vi phân vi mô phức tạp thành các phương trình uốn độc lập mà không cần giả thiết đơn giản hóa gượng ép, đồng thời chuẩn hóa hệ số hiệu chỉnh cắt $\kappa_s(z)$ biến thiên chính xác theo từng cấu hình phân lớp.

2. Điểm cải tiến cốt lõi của phương pháp luận số so với các nghiên cứu tiền nhiệm?

So với các phần tử hữu hạn cấp cao đòi hỏi hàm dạng liên tục $C^1$ phức tạp của Zenkour [50] hay các phương pháp không lưới tốn kém chi phí tính toán của Liew và cộng sự [122], luận án đã kết hợp kỹ thuật làm trơn biến dạng (S-FEM) với phần tử nội suy ten-sơ hỗn hợp (MITC3) chỉ sử dụng hàm dạng liên tục $C^0$. Giải pháp này triệt tiêu hoàn toàn hiện tượng khóa cắt, cho độ chính xác tiệm cận lời giải giải tích Navier/Ritz ngay cả trên lưới tam giác thô và biến dạng méo mó.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm số là gì?

Phát hiện mang tính phản trực giác là sự đảo chiều sai số của hệ số hiệu chỉnh cắt truyền thống ($\kappa_s = 5/6$) giữa các loại tấm sandwich FGM: Với tấm lõi cứng (Loại B), mô hình $5/6$ cho kết quả an toàn giả tạo (đánh giá độ võng nhỏ hơn thực tế), trong khi với tấm lõi mềm (Loại C), mô hình $5/6$ lại đánh giá quá cứng kết cấu, làm tăng nguy cơ phá hủy kết cấu trong thực tế do đánh giá thấp độ võng nguy hiểm.

4. Luận án có cung cấp đầy đủ giao thức để tái lập kết quả (Replication Protocol) không?

Có. Luận án trình bày chi tiết từng bước: ma trận cấu tạo $C_{ij}(z)$, tọa độ điểm Gauss 1 điểm và 2 điểm, công thức tính ma trận Jacobian, biểu thức đại số của các ma trận con $\mathbf{K}{ij}^{e}, \mathbf{M}{ij}^{e}$, cùng với các thông số chuẩn hóa không thứ nguyên và bảng dữ liệu kiểm chứng với các tác giả Mohammadi [131], Praveen & Reddy [28], giúp bất kỳ nhà nghiên cứu nào cũng có thể lập trình tái lập chính xác $100%$ kết quả.

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

Khung chương trình 10 năm định hướng vào việc chuyển dịch từ cơ cấu tĩnh/dao động tuyến tính sang phân tích đa trường nhiệt - cơ - áp điện phi tuyến, tích hợp phương pháp đẳng hình học (IGA) với cấu trúc micro-nano dựa trên lý thuyết đàn hồi phi cục bộ Eringen, và kết hợp công nghệ in 3D FGM để thực nghiệm hóa và chuẩn hóa các mô hình số.

Kết luận

  1. Luận án đã thiết lập thành công hệ thống lý thuyết toàn diện cho tấm FGM và tấm sandwich FGM, kết hợp hoàn hảo giữa lý thuyết biến dạng cắt bậc nhất cải tiến (FSDT), lý thuyết biến dạng cắt bậc cao (HSDT) và lý thuyết tiếp cận 3 chiều (Quasi-3D).
  2. Phát triển công thức giải tích xác định chính xác vị trí mặt trung hòa vật lý $z_0$ và độ cứng cắt cải tiến $H$, loại bỏ hoàn toàn tương tác kéo - uốn nhân tạo và sửa đổi triệt để nhược điểm của hệ số hiệu chỉnh cắt $\kappa_s = 5/6$ truyền thống.
  3. Đề xuất thành công họ phần tử hữu hạn làm trơn khử khóa cắt tiên tiến CS-MITC3, ES-MITC3 và NS-MITC3, giải quyết dứt điểm hiện tượng khóa cắt trong phân tích tấm từ mỏng đến dày cho cả bài toán tuyến tính và phi tuyến hình học Von Kármán.
  4. Đánh giá định lượng sâu sắc độ phân tán của 6 mô hình đồng nhất hóa vi mô, chỉ ra chênh lệch tới $48%$ giữa Voigt và Reuss, thiết lập cẩm nang lựa chọn mô hình vi cơ học chuẩn xác cho kỹ thuật tính toán kết cấu.
  5. Xây dựng bộ công cụ tính toán số tự động hóa hoàn toàn trên MATLAB, kiểm chứng độ chính xác vượt trội qua hàng trăm ví dụ số về phân tích tĩnh, tải trọng tới hạn ổn định và tần số dao động tự do dưới tác dụng của tải trọng cơ - nhiệt và nền đàn hồi Winkler-Pasternak.
  6. Mở ra các hướng nghiên cứu đột phá về cơ học tính toán cho kết cấu composite thông minh thế hệ mới, khẳng định vị thế học thuật vững chắc và đóng góp lâu dài cho nền khoa học kỹ thuật công trình hiện đại.