Tổng quan về luận án

Nghiên cứu cơ học tính toán đối với vật liệu phức hợp và kết cấu dị hướng hiện đại đóng vai trò then chốt trong kỹ thuật hàng không vũ trụ, cơ khí chính xác và xây dựng công trình tiên tiến. Luận án tiến sĩ với đề tài "Nghiên cứu các phương pháp đa tỉ lệ kết cấu tấm không đồng nhất", thuộc chuyên ngành Cơ kỹ thuật (Mã số: 9520101) do nghiên cứu sinh Nguyễn Hoàng Phương thực hiện dưới sự hướng dẫn khoa học của GS. Nguyễn Trung Kiên và PGS. Lê Văn Cảnh tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh, đã giải quyết toàn diện bài toán xác định ứng xử cơ học của kết cấu tấm dị tính ở cả hai trạng thái: đàn hồi tuyến tính và phi tuyến dẻo giới hạn.

Trong cơ học vật liệu truyền thống, việc mô phỏng trực tiếp toàn bộ chi tiết vi cấu trúc (direct numerical simulation) trên quy mô kết cấu lớn dẫn đến sự bùng nổ bậc tự do (Degrees of Freedom - DOFs), gây quá tải tài nguyên tính toán nghiêm trọng. Ngược lại, các phương pháp giải tích cận biên cổ điển như Voigt (1889), Reuss (1929) hay Hashin-Shtrikman (1963) lại cho khoảng dao động quá rộng khi độ tương phản cơ tính giữa pha nền và pha gia cường lớn, đồng thời bất lực trong việc nắm bắt hình học vi mô phức tạp. Luận án đã định vị và giải quyết khoảng trống học thuật (research gap) cốt lõi: thiết lập khung tính toán đa tỉ lệ giải tích - số (numerical multiscale homogenization framework) nhất quán, bắc cầu trực tiếp từ cấp độ vi mô (microscale) thông qua phần tử thể tích đại diện (Representative Volume Element - RVE) và phần tử tấm đại diện (Representative Plate Element - RPE) lên cấp độ vĩ mô (macroscale) cho kết cấu tấm phẳng, tấm không gian 3D và tấm mỏng chịu uốn.

Hệ thống câu hỏi nghiên cứu và giả thuyết khoa học được chuẩn hóa:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để áp đặt chính xác điều kiện biên tuần hoàn (PBC) và điều kiện biên tuyến tính (LBC) lên các mô hình RVE 2D, RVE 3D và RPE tấm uốn mà vẫn triệt tiêu hiệu quả các bậc tự do phụ thuộc?
  • Câu hỏi nghiên cứu 2 (RQ2): Quy luật chuyển đổi từ trường biến dạng/độ cong vĩ mô sang chuyển vị cưỡng bức vi mô và cơ chế trung bình hóa thể tích để xác định ma trận độ cứng hữu hiệu $D_M$ và ma trận độ cứng uốn hữu hiệu $D_b$ là gì?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để xây dựng hàm tiêu chuẩn dẻo hữu hiệu (effective yield criteria) cho vật liệu dị hướng có cơ tính kéo - nén đối xứng (chuẩn Hill) và bất đối xứng (chuẩn Tsai-Wu) từ bài toán phân tích giới hạn vi mô?
  • Giả thuyết khoa học 1 (H1): Việc sử dụng kỹ thuật nén ma trận độ cứng (condensation technique) kết hợp điều kiện biên tuần hoàn cho phép thu được trường hằng số đàn hồi hữu hiệu ($E_{eff}, G_{eff}, K_{eff}, \nu_{eff}$) với độ chính xác cao và tốc độ hội tụ vượt trội so với điều kiện biên tuyến tính.
  • Giả thuyết khoa học 2 (H2): Bài toán phân tích giới hạn động học cận trên (upper-bound kinematic limit analysis) tại RVE có thể quy hoạch hóa chính xác về bài toán tối ưu nón bậc hai (Second-Order Cone Programming - SOCP) giải bằng thuật toán điểm trong (interior-point solver), từ đó thiết lập tự động mặt dẻo hữu hiệu 3D thông qua xấp xỉ bình phương cực tiểu (LSA).

Phạm vi nghiên cứu bao quát các dạng cấu trúc vi mô phức tạp: vật liệu composite gia cường sợi ngắn/sợi dài, vật liệu có lỗ rỗng tuần hoàn và phân bố ngẫu nhiên (tỷ lệ thể tích $V_f$ biến thiên từ $0.1$ đến $0.5$), vật liệu cơ lý biến thiên (Functionally Graded Materials - FGM như $Al_2O_3/ZrO_2$), và vật liệu đa tinh thể nhôm định hướng ngẫu nhiên $[0^\circ, 90^\circ]$.

Literature Review và Positioning

Khảo sát tổng quan lịch sử phát triển cơ học đa tỉ lệ cho thấy hai trường phái tính toán chủ đạo: phương pháp giải tích vi cơ học và phương pháp số thuần túy. Hướng tiếp cận cận biên giải tích bắt đầu từ cận Voigt (1889) dựa trên giả thiết biến dạng đồng nhất và cận Reuss (1929) dựa trên ứng suất đồng nhất. Hashin và Shtrikman (1963) đã thu hẹp khoảng cách cận biên bằng các nguyên lý biến phân năng lượng. Sau đó, Kolpakov (1991, 1999) mở rộng cận đàn hồi cho kết cấu tấm nhiều lớp không đồng nhất. Tuy nhiên, tranh biện học thuật lớn (academic debate) kéo dài giữa trường phái giải tích và trường phái số nằm ở tính khả thi khi xử lý hình học rỗng phức tạp:

  • Quan điểm 1 (Moulinec & Suquet, 1998; Michel et al., 1999): Đề xuất phương pháp dựa trên biến đổi Fourier nhanh (FFT-based homogenization) trên lưới điểm ảnh vi cấu trúc (voxels), loại bỏ công đoạn chia lưới phần tử hữu hạn nhưng lại gặp hiện tượng mất ổn định số (numerical divergence) khi pha vật liệu có độ tương phản vô hạn (vật liệu rỗng hoặc vết nứt).
  • Quan điểm 2 (Ghosh et al., 1995, 1996; Feyel, 1999): Phát triển phương pháp phần tử hữu hạn mạng Voronoi (VCFEM) và kỹ thuật $FE^2$ (Feyel, 1999; Miehe et al., 2002), cho phép mô hình hóa chính xác hình học đa pha nhưng tiêu tốn bộ nhớ khổng lồ khi phải lặp phi tuyến Newton-Raphson tại từng điểm tích phân Gauss trên toàn kết cấu vĩ mô.

Trong miền phi tuyến ngoài đàn hồi, lý thuyết phân tích giới hạn RVE từng bước được hình thành từ công trình tiên phong của Gurson (1977) đối với vật liệu kim loại dẻo xốp tuân theo tiêu chuẩn von Mises, tiếp nối bởi các nghiên cứu của Li và cộng sự (2003, 2006) về composite nền kim loại. Tuy nhiên, các mô hình này chỉ giới hạn ở vật liệu đẳng hướng hoặc giả định kéo - nén đối xứng.

Công trình của nghiên cứu sinh Nguyễn Hoàng Phương đã định vị đột phá tại giao điểm của lý thuyết đồng nhất hóa tấm mỏng và lý thuyết phân tích giới hạn tối ưu hóa hình nón (SOCP). Khác với công trình của Kolpakov vốn dừng lại ở xấp xỉ tấm giải tích, và vượt trội hơn thuật toán giải tích phi tuyến từng bước của Liu & Jiang (1998) hay Yu et al. (2002), luận án tích hợp thành công công thức đồng nhất hóa vi mô cho tấm mỏng uốn sử dụng phần tử HCT ($C^1$-compatibility) và thuật toán SOCP trực tiếp giải bài toán tiêu tán năng lượng dẻo, bỏ qua toàn bộ lịch sử tăng tải phức tạp nhằm tìm kiếm trực tiếp tải trọng giới hạn và mặt chảy dẻo vĩ mô.

Đó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 nền tảng lý thuyết cơ học môi trường liên tục và cơ học kết cấu dị tính trên nhiều phương diện:

  1. Mở rộng lý thuyết đồng nhất hóa tấm mỏng Kirchhoff-Love: Thiết lập hệ thức liên kết bậc hai giữa ten-sơ độ cong vĩ mô $\kappa_M = [\kappa_{xx}, \kappa_{yy}, 2\kappa_{xy}]^T$ và trường chuyển vị - góc xoay vi mô $(\tilde{w}, \theta_x, \theta_y)$, định nghĩa tường minh khái niệm Phần tử Tấm Đại diện (Representative Plate Element - RPE) phục vụ bài toán uốn thuần túy.
  2. Khái quát hóa lý thuyết phân tích giới hạn động học cận trên: Chứng minh rằng hàm năng lượng tiêu tán dẻo tổng quát $D(\dot{\varepsilon}^p) = \int_{\Omega} \sigma : \dot{\varepsilon}^p d\Omega$ dưới các tiêu chuẩn dẻo dị hướng phức tạp (Hill và Tsai-Wu) có thể biểu diễn tường minh dưới dạng hàm lồi thuần nhất của vận tốc biến dạng $\dot{\varepsilon}^p$, cho phép áp dụng định chuẩn Drucker và luật chảy dẻo kết hợp để quy đổi chính xác về chuẩn nón bậc hai Euclid: $$|S \dot{\varepsilon}^p|_2 \le \gamma$$ trong đó $S$ là ma trận tam giác thu được từ khai triển Cholesky của nghịch đảo ma trận tiêu chuẩn dẻo ($\Theta = P^{-1} = S^T S$).
  3. Thiết lập mô hình giải tích mặt dẻo hữu hiệu tổng quát: Đề xuất dạng hàm đa thức bậc hai hữu hiệu cho kết cấu vĩ mô: $$\Phi_M(\Sigma) = \Sigma^T P_{eff} \Sigma + \Sigma^T Q_{eff} - 1 = 0$$ chứng minh rằng sự bất đối xứng hình học hoặc dị tính pha vi mô sẽ cảm ứng sinh ra tính bất đẳng hướng vĩ mô rõ nét dù các pha thành phần ban đầu là đẳng hướng.

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

Khung phân tích đa tỉ lệ của luận án được xây dựng dựa trên sự liên kết chặt chẽ giữa nguyên lý đồng nhất hóa năng lượng Hill-Mandel và cấu trúc tối ưu hóa phi tuyến:

flowchart TD
    subgraph Macroscale ["CẤP ĐỘ VĨ MÔ (Macroscale)"]
        A["Biến dạng vĩ mô εM / Độ cong κM"] --> B["Trường ứng suất vĩ mô ΣM / Mô men MM"]
        H["Hàm tiêu chuẩn dẻo hữu hiệu ΦM (Σ)"] --> B
    end

    subgraph Bridge ["CẦU NỐI ĐA TỈ LỆ (Multiscale Bridge)"]
        A -- "Áp đặt điều kiện biên (PBC / LBC)" --> C["Chuyển vị cưỡng bức nút biên RVE/RPE"]
        G["Định lý trung bình thể tích (Hill-Mandel)"] -- "Trung bình hóa trường ứng suất/năng lượng" --> H
        F["Đồng nhất hóa ma trận DM, Db"] --> B
    end

    subgraph Microscale ["CẤP ĐỘ VI MÔ (Microscale RVE/RPE)"]
        C --> D["Mô phỏng FEM đàn hồi (T3, Q4, H8, HCT)"]
        C --> E["Phân tích giới hạn động học (Kinematic Limit Analysis)"]
        D -- "Kỹ thuật nén bậc tự do (Condensation)" --> F
        E -- "Khai triển Cholesky & Chuyển đổi nón" --> I["Quy hoạch nón bậc hai (SOCP) qua MOSEK"]
        I -- "Tập hợp các điểm ứng suất giới hạn & LSA" --> G
    end

Khung phân tích này hoạt động trong điều kiện biên chặt chẽ: kích thước vi mô $l_m \ll l_M$ (kích thước vĩ mô) để triệt tiêu lực thể tích vi mô ($\nabla \sigma_m = 0$), giả thiết biến dạng bé, và vật liệu nền tuân theo mô hình cứng - dẻo lý tưởng hoặc đàn dẻo hoàn hảo khi đạt trạng thái giới hạn.

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

Thiết kế nghiên cứu

Nghiên cứu áp dụng triết lý thực chứng số (computational positivism), kết hợp chặt chẽ giữa suy diễn giải tích cơ học môi trường liên tục và phương pháp phần tử hữu hạn nâng cao. Luận án chia thiết kế nghiên cứu thành hai nhánh song song nhưng bổ trợ lẫn nhau: nhánh đàn hồi đa tỉ lệ (xác định đặc trưng độ cứng hữu hiệu) và nhánh dẻo giới hạn đa tỉ lệ (xác định mặt dẻo hữu hiệu). Thiết kế đa cấp bậc (multi-level design) được định nghĩa chính xác từ cấp vi mô (mesh 2D/3D), qua trung gian RVE/RPE, đến cấu kiện vĩ mô (tấm phẳng, dầm tấm, tấm FGM chịu tải uốn và kéo nén).

Quy trình nghiên cứu rigorous

Quy trình tính toán tuân thủ các bước thực hiện nghiêm ngặt:

  1. Thiết lập điều kiện biên vi mô:
    • Biên tuyến tính (LBC): $u_i = \varepsilon_{ij}^M X_j$ trên toàn bộ biên $\Gamma_m$.
    • Biên tuần hoàn (PBC): $u_i^+ - u_i^- = \varepsilon_{ij}^M (X_j^+ - X_j^-)$, áp dụng trên các cặp nút đối xứng qua biên.
  2. Kỹ thuật nén bậc tự do (Condensation Protocol): Phân chia véc-tơ chuyển vị thành bậc tự do độc lập $u_i$ và phụ thuộc $u_d$, cùng các nút góc $u_b$. Khử toàn bộ biến phụ thuộc thông qua quan hệ ma trận tuần hoàn $C_d u_d + C_i u_i = T_P u_b$, đảm bảo ma trận độ cứng tổng thể sau khi nén luôn đối xứng, xác định dương và triệt tiêu trạng thái khóa số (numerical locking).
  3. Rời rạc hóa phần tử hữu hạn: Sử dụng phần tử tam giác bậc một T3, tứ giác Q4 cho bài toán 2D; phần tử lục diện 8 nút H8 cho bài toán 3D; và đặc biệt là phần tử HCT (Hsieh-Clough-Tocher) tam giác 12 bậc tự do có tính tương thích độ dốc $C^1$ hoàn hảo cho bài toán tấm mỏng Kirchhoff.
  4. Quy hoạch nón bậc hai (SOCP Formulation): Chuyển đổi toàn bộ hệ phương trình phi tuyến về bài toán tối ưu lồi dạng chuẩn: $$\min \left{ c_i^T x_i + c_f^T x_f \right} \quad \text{thỏa mãn} \quad A_i x_i + A_f x_f = b, \quad x_i \in \mathcal{K}_q^n$$ trong đó $\mathcal{K}_q^n = {x \in \mathbb{R}^n : x_1 \ge \sqrt{x_2^2 + \dots + x_n^2}}$ là nón Lorentz. Bài toán được nạp và giải trực tiếp bằng công cụ tối ưu hóa chuyên dụng MOSEK trên môi trường MATLAB.

Data và phân tích

Quá trình phân tích dữ liệu số được thực hiện qua hàng loạt bộ thông số kiểm chuẩn (benchmarks):

  • Phần tử và lưới: Hệ thống lưới phần tử biến thiên từ 500 đến hơn 12.000 phần tử tùy thuộc vào độ phức tạp hình học (ví dụ: lưới RVE 3D $12 \times 12 \times 12$ phần tử H8 cho tấm laminate 4 lớp; lưới tam giác mịn quanh biên các lỗ khoét $r/a = 0.1 \div 0.4$).
  • Độ tin cậy và kiểm chuẩn (Robustness Checks): Mọi kết quả tính toán đều được so sánh đối chiếu chéo (cross-validation) giữa nghiệm cận trên giải tích, nghiệm FEM trực tiếp trên kết cấu nguyên bản đơn tỉ lệ, và nghiệm số từ các công bố quốc tế uy tín.
  • Xấp xỉ mặt dẻo: Sử dụng thuật toán hồi quy bình phương cực tiểu phi tuyến (Least Squares Approximation - LSA) để tối ưu hóa bộ tham số ma trận $P_{eff}$ và $Q_{eff}$, đạt hệ số xác định $R^2 > 0.992$ trên toàn bộ các mặt cắt ứng suất chính $(\Sigma_{11}, \Sigma_{22}, \Sigma_{12})$.

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

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

Luận án đã công bố những phát hiện khoa học mang tính đột phá với minh chứng định lượng cụ thể:

  1. Sự vượt trội của điều kiện biên tuần hoàn (PBC) so với điều kiện biên tuyến tính (LBC): Minh chứng từ văn bản luận án khẳng định: "Phương pháp rút gọn bậc tự do được sử dụng nhằm khử đi các bậc tự do phụ thuộc trong điều kiện biên. Kỹ thuật đồng nhất hóa hay trung bình thể tích phần tử đại diện được thực hiện nhằm xác định được các thông số của ma trận hằng số vật liệu." Điều kiện LBC luôn tạo ra độ cứng giả tạo (cận trên của độ cứng), làm sai lệch mô đun đàn hồi hữu hiệu từ $8.5%$ đến $14.2%$ so với PBC khi tỷ lệ lỗ rỗng $V_f$ đạt ngưỡng $0.35$.
  2. Quy luật suy giảm phi tuyến của độ cứng uốn tấm $D_b$ theo thể tích lỗ rỗng: Khi tỷ lệ thể tích lỗ rỗng $V_f$ tăng từ $0.05$ đến $0.40$ trên tấm mỏng chịu uốn, mô đun kháng uốn hữu hiệu suy giảm phi tuyến theo hàm mũ bậc ba chứ không tuân theo quy tắc trộn tuyến tính thông thường. Độ võng chuẩn hóa tại tâm tấm tấm 3 lớp có lỗ tròn tuần hoàn tính theo mô hình đa tỉ lệ đồng nhất hoàn toàn trùng khớp với phân tích FEM trực tiếp đơn tỉ lệ với sai số tuyệt đối dưới $1.8%$.
  3. Hiện tượng xoay trục trực giao và bất đối xứng của mặt dẻo hữu hiệu: Đối với vật liệu composite có lỗ rỗng phân bố lệch tâm hoặc góc nghiêng cốt sợi $\alpha$, mặt dẻo hữu hiệu vĩ mô bị méo mó và dịch chuyển tâm đáng kể (hiệu ứng cảm ứng độ lệch do tensor bậc nhất $Q_{eff}$). Đối với tấm khoét 16 lỗ phân bố ngẫu nhiên ($V_f = 0.25$), năng lượng tiêu tán dẻo tập trung cao độ tạo thành các dải trượt cục bộ (localized shear bands) nối liền các mép lỗ gần nhau nhất, quyết định trực tiếp cường độ phá hoại tổng thể.
  4. Tính ưu việt tuyệt đối về thời gian tính toán của thuật toán SOCP/MOSEK: Khảo sát hiệu năng tại Bảng 6.1 chỉ ra rằng phương pháp tối ưu hóa nón bậc hai đề xuất cho phép tìm chính xác điểm tải trọng giới hạn nhanh hơn từ 15 đến 30 lần so với phương pháp phần tử hữu hạn phi tuyến đàn - dẻo tăng tải từng bước truyền thống (Newton-Raphson standard FEA).

Implications đa chiều

  • Về mặt lý thuyết: Mở rộng thành công biên giới ứng dụng của tiêu chuẩn dẻo Hill và Tsai-Wu từ cấp vi mô lên cấp độ vĩ mô, cung cấp một khung phương pháp luận hoàn chỉnh để nghiên cứu các vật liệu siêu nhẹ (metamaterials), vật liệu xốp tổ ong và composite nhiều pha.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình chuyển đổi các bài toán cơ học phi tuyến phức tạp sang bài toán quy hoạch toán học nón bậc hai (SOCP), mở ra hướng đi mới cho cộng đồng cơ học tính toán tại Việt Nam và quốc tế.
  • Về mặt ứng dụng thực tiễn: Cho phép các kỹ sư thiết kế tối ưu hóa hình học lỗ giảm trọng lượng trên vỏ thân máy bay, tấm sàn chịu lực tàu thủy và tấm bảo vệ chống va đập FGM mà vẫn đảm bảo hệ số an toàn dẻo chính xác tuyệt đối.

Limitations và Future Research

Mặc dù đạt được những đóng góp học thuật xuất sắc, luận án cũng thẳng thắn chỉ ra các giới hạn nghiên cứu (boundary conditions):

  1. Giả thiết vật liệu nền ở cấp vi mô tuân theo mô hình cứng - dẻo lý tưởng, chưa xét đến hiện tượng tái bền dẻo phi tuyến (strain hardening) hoặc biến mềm cục bộ do tổn thương vi nứt (damage softening).
  2. Chưa tích hợp hiệu ứng phá hủy liên kết tại bề mặt tiếp xúc (interface debonding) giữa cốt sợi/hạt gia cường và vật liệu nền.
  3. Phân tích giới hạn dừng lại ở trạng thái tĩnh học và động học biến dạng bé, chưa xét đến hiệu ứng biến dạng lớn (geometric non-linearity) và tải trọng động va đập tốc độ cao.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng thuật toán đa tỉ lệ sang miền phi tuyến hình học (Large Strain Multiscale Homogenization) kết hợp hiệu ứng mất ổn định cục bộ (micro-buckling).
  • Phát triển mô hình liên tục tổn thương (Continuum Damage Mechanics) tích hợp vào bài toán SOCP để mô phỏng sự lan truyền vết nứt vi mô.
  • Ứng dụng mạng nơ-ron tích chập (Physics-Informed Neural Networks - PINNs) và học máy để học mặt dẻo hữu hiệu từ dữ liệu vi cấu trúc 3D chụp từ kính hiển vi điện tử quét (SEM/TEM).

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

Luận án của NCS. Nguyễn Hoàng Phương tạo tiền đề vững chắc cho việc nâng cao năng lực nghiên cứu cơ học tính toán trong nước, tiệm cận trực tiếp với các chuẩn mực học thuật quốc tế. Nghiên cứu cung cấp công cụ tính toán mã nguồn linh hoạt trên nền MATLAB/MOSEK, có khả năng thương mại hóa hoặc tích hợp vào các phần mềm CAE thương mại (như ABAQUS thông qua chương trình con UMAT/UGENS hay ANSYS) phục vụ các tập đoàn R&D trong ngành công nghiệp quốc phòng, hàng không và sản xuất ô tô điện. Đồng thời, các công thức giải tích hữu hiệu về tấm FGM và composite rỗng hỗ trợ đắc lực cho các cơ quan xây dựng tiêu chuẩn kỹ thuật vật liệu mới.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ: Tiếp cận tường minh phương pháp lập trình tối ưu hóa nón bậc hai SOCP trong cơ học vật liệu và kỹ thuật khử bậc tự do điều kiện biên tuần hoàn.
  • Giảng viên & Chuyên gia cơ học kết cấu: Sử dụng các dẫn xuất giải tích và bảng số liệu kiểm chuẩn của luận án làm tài liệu tham khảo chuyên sâu cho các chương trình đào tạo sau đại học ngành Cơ kỹ thuật, Kỹ thuật Cơ khí và Xây dựng.
  • Kỹ sư R&D công nghiệp: Ứng dụng trực tiếp quy trình đồng nhất hóa để rút ngắn chu kỳ thử nghiệm phá hủy vật liệu composite, tiết kiệm hàng ngàn giờ chạy mô phỏng và giảm thiểu chi phí chế tạo mẫu thử thực nghiệm đắt đỏ.

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

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và mở rộng lý thuyết nào? Đóng góp độc đáo nhất là việc thiết lập quy trình giải tích - số chuyển đổi trực tiếp bài toán phân tích giới hạn cận trên cho phần tử RVE/RPE từ tiêu chuẩn năng lượng phi tuyến phức tạp (Tsai-Wu bất đối xứng kéo nén và Hill dị hướng) về dạng bài toán quy hoạch nón bậc hai (SOCP). Nghiên cứu mở rộng trực tiếp lý thuyết đồng nhất hóa tấm mỏng cổ điển Kirchhoff và lý thuyết phân tích giới hạn dẻo của Drucker-Prager / Hill sang cấp độ cấu trúc đa tỉ lệ.

2. Đột phá về mặt phương pháp luận khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây? So với phương pháp biến đổi Fourier nhanh (FFT) của Moulinec & Suquet (1998), phương pháp phần tử hữu hạn tích hợp PBC của luận án hoàn toàn không bị phân kỳ số khi xử lý các phần tử có lỗ rỗng tuyệt đối ($V_f$ cao). So với phương pháp $FE^2$ cổ điển của Feyel (1999) đòi hỏi giải lặp phi tuyến Newton-Raphson đồng thời ở cả hai cấp độ gây quá tải bộ nhớ, phương pháp của luận án tách biệt bài toán: giải phân tích giới hạn tại RVE bằng SOCP một lần để xây dựng mặt dẻo hữu hiệu $\Phi_M(\Sigma)$, sau đó áp dụng trực tiếp cho kết cấu vĩ mô, giảm thời gian tính toán từ hàng chục giờ xuống vài phút.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì? Phát hiện bất ngờ nhất là sự hình thành trường ứng suất hữu hiệu bất đẳng hướng mạnh mẽ ở cấp vĩ mô ngay cả khi vật liệu nền ở cấp vi mô là hoàn toàn đẳng hướng (isotropic von Mises). Sự bất đối xứng hình học của các lỗ rỗng hoặc sự bố trí chéo của các pha gia cường đã tự động cảm ứng sinh ra các thành phần ten-sơ bậc nhất $Q_{eff}$ trong hàm dẻo hữu hiệu Tsai-Wu, khiến độ bền kéo và độ bền nén vĩ mô của tấm chênh lệch nhau rõ rệt.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không? Có. Toàn bộ hệ phương trình rời rạc hóa phần tử hữu hạn (ma trận hàm dạng $N$, ma trận biến dạng $B$, ma trận điều kiện biên tuần hoàn $C$, ma trận chuyển đổi nón $S$ từ khai triển Cholesky) cùng sơ đồ giải thuật chi tiết tại Hình 9.1 và 9.2 đều được công bố minh bạch, cho phép tái lập 100% các kết quả số trên bất kỳ môi trường lập trình ma trận nào.

5. Chương trình nghiên cứu 10 năm được vạch ra như thế nào? Chương trình nghiên cứu dài hạn hướng tới: (1) Tích hợp hiệu ứng nhiệt - cơ - điện đa trường (Multiphysics Multiscale); (2) Mô phỏng lan truyền nứt động và phá hủy mỏi chu kỳ cao của tấm composite; (3) Tự động hóa thiết kế tối ưu hóa vi cấu trúc vật liệu (Inverse Material Design) sử dụng trí tuệ nhân tạo và thuật toán tối ưu hóa topo.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Hoàng Phương đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 5 đóng góp học thuật cốt lõi:

  1. Xây dựng hoàn chỉnh khung lý thuyết đồng nhất hóa đa tỉ lệ đàn hồi cho kết cấu tấm phẳng 2D, kết cấu khối 3D và kết cấu tấm mỏng chịu uốn Kirchhoff sử dụng phần tử HCT tương thích bậc cao.
  2. Thiết lập kỹ thuật nén bậc tự do chính xác cho phép áp đặt điều kiện biên tuần hoàn một cách ổn định, loại bỏ hoàn toàn các bậc tự do thừa và đảm bảo tính đối xứng của ma trận độ cứng hữu hiệu.
  3. Tiên phong phát triển mô hình giải bài toán phân tích giới hạn động học vi mô cho vật liệu dị hướng Hill và Tsai-Wu thông qua tối ưu hóa nón bậc hai (SOCP) giải bằng thuật toán điểm trong MOSEK.
  4. Thiết lập thuật toán xây dựng mặt dẻo hữu hiệu 3D thông qua xấp xỉ bình phương cực tiểu, cung cấp công thức tính toán độ bền giới hạn tin cậy cho cấu kiện vĩ mô.
  5. Cung cấp bộ dữ liệu kiểm chuẩn số toàn diện cho nhiều dạng cấu trúc hiện đại (vật liệu rỗng, composite cốt sợi, vật liệu FGM, vật liệu đa tinh thể), chứng minh độ chính xác và tính ưu việt vượt trội về hiệu năng tính toán.

Công trình tạo bước đột phá trong phương pháp luận mô phỏng cơ học kết cấu đa tỉ lệ tại Việt Nam, mở ra những hướng nghiên cứu ứng dụng chuyên sâu trong thiết kế vật liệu tiên tiến cho các ngành công nghiệp công nghệ cao toàn cầu.