Tổng quan luận án

Luận án tiến sĩ kỹ thuật của tác giả Nguyễn Hoàng Phương 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), được thực hiện tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh năm 2021 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. Hội đồng đánh giá luận án gồm PGS. Trương Tích Thiện (Phản biện 1), TS. Thái Hoàng Chiến (Phản biện 2) và TS. Phùng Văn Phúc (Phản biện 3).

Tính cấp thiết và khoảng trống nghiên cứu

Trong kỹ thuật hiện đại, các vật liệu mới có cấu trúc phức tạp như vật liệu hỗn hợp gia cường cốt sợi (composite), vật liệu cơ lý biến thiên (Functionally Graded Materials - FGM), vật liệu đa tinh thể và vật liệu có lỗ rỗng được ứng dụng nhằm khai thác các thuộc tính cơ lý về khả năng chịu lực, va đập, chống mài mòn và chịu nhiệt. Vật liệu trong tự nhiên và nhân tạo mang bản chất đa tỉ lệ do kết cấu vĩ mô được cấu tạo từ các vi cấu trúc ở cấp độ nhỏ hơn.

Khi tính toán kết cấu tấm, việc mô phỏng trực tiếp toàn bộ các pha vi mô trên toàn miền kết cấu vĩ mô bằng phương pháp phần tử hữu hạn thông thường dẫn đến sự bùng nổ dữ liệu, quá tải bộ nhớ và tốn kém thời gian tính toán. Ngược lại, các phương pháp giải tích truyền thống như lý thuyết tấm nhiều lớp (laminate plate theory) hoặc quy tắc hỗn hợp thông thường ("the rule of mixture") lại gặp khó khăn khi mô tả chính xác sự phân bố phức tạp của các pha vật liệu theo chiều dày và mặt phẳng tấm. Đối với trạng thái ngoài miền đàn hồi, việc xác định trực tiếp tải trọng giới hạn và miền ứng suất giới hạn (mặt chảy dẻo hữu hiệu) cho vật liệu không đồng nhất bất đẳng hướng chịu tải trọng phức tạp vẫn là thách thức lớn. Do đó, việc xây dựng phương pháp đa tỉ lệ kết hợp đồng nhất hóa và phân tích giới hạn trực tiếp cho kết cấu tấm không đồng nhất là yêu cầu kỹ thuật cấp thiết.

Mục tiêu nghiên cứu

Mục tiêu tổng quát của luận án là thiết lập và phát triển các phương pháp tính toán đa tỉ lệ cho kết cấu tấm không đồng nhất ở cả hai trạng thái: trong miền đàn hồi tuyến tính và ngoài miền đàn hồi dẻo. Mục tiêu cụ thể được phân thành hai nhóm:

  1. Nghiên cứu trong miền đàn hồi: Xác định các thông số đàn hồi hữu hiệu (mô đun đàn hồi, hệ số Poisson, mô đun trượt, mô đun khối, ma trận độ cứng uốn) thông qua kỹ thuật đồng nhất hóa và trung bình thể tích phần tử đại diện ở cấp vi mô. Phương pháp được triển khai cho ba dạng mô hình đại diện gồm: phần tử tấm phẳng hai chiều chịu lực màng, phần tử thể tích không gian ba chiều (3D) và phần tử tấm mỏng chịu uốn.
  2. Nghiên cứu ngoài miền đàn hồi: Xây dựng quy trình xác định miền cường độ hữu hiệu và hàm tiêu chuẩn dẻo hữu hiệu ở cấp vĩ mô bằng cách kết hợp kỹ thuật đồng nhất hóa với lý thuyết phân tích giới hạn động học trên phần tử thể tích đại diện, giải quyết cho vật liệu bất đẳng hướng tuân theo dạng tiêu chuẩn Hill (đối xứng kéo - nén) và dạng tiêu chuẩn Tsai-Wu (bất đối xứng kéo - nén).

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Kết cấu tấm không đồng nhất và các dạng vi cấu trúc đại diện bao gồm vật liệu gia cường cốt sợi (sợi tròn, sợi chữ nhật, sợi ngắn, sợi dài), vật liệu có lỗ rỗng (lỗ đơn tròn, lỗ chữ nhật, hai lỗ, nhiều lỗ phân bố đều và ngẫu nhiên), vật liệu cơ lý biến thiên FGM (hệ kim loại $\text{Al}_2\text{O}_3$ - gốm $\text{ZrO}_2$) và vật liệu đa tinh thể nhôm (Al).
  • Phạm vi nghiên cứu:
    • Về mặt hình học và cơ học: Nghiên cứu xét các bài toán tấm phẳng 2D ứng suất phẳng, kết cấu khối 3D nhiều lớp (cross-ply, angle-ply, xen kẽ) và tấm mỏng chịu uốn theo lý thuyết tấm cổ điển Kirchhoff (CPT) sử dụng phần tử tấm tương thích $C^1$ HCT (Hsieh-Clough-Tocher 12 bậc tự do). Giả thiết chuyển vị và biến dạng bé.
    • Về điều kiện biên vi mô: Áp dụng điều kiện biên tuyến tính và điều kiện biên tuần hoàn. Sử dụng kỹ thuật nén bậc tự do (condensation) để khử các bậc tự do phụ thuộc. Mối liên hệ giữa biến dạng vĩ mô và chuyển vị biên vi mô là bậc nhất (cho tấm phẳng và mô hình 3D) và bậc hai (đối với biến dạng cong trong bài toán tấm uốn).
    • Về mô hình vật liệu ngoài đàn hồi: Các pha vi mô được giả thuyết là vật liệu cứng dẻo lý tưởng (bỏ qua hiện tượng tái bền và làm mềm), tuân theo luật chảy dẻo kết hợp. Vật liệu nền vi mô được xét theo các tiêu chuẩn von Mises, Hill và Tsai-Wu.
    • Về mặt toán học và tính toán: Bài toán phân tích giới hạn cận trên được đưa về bài toán tối ưu hóa nón bậc hai (Second-Order Cone Programming - SOCP) và giải bằng công cụ MOSEK. Mặt dẻo hữu hiệu vĩ mô được xấp xỉ bằng phương pháp bình phương cực tiểu (Least Squares Approximation - LSA).

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước

Luận án đã tổng quan và phân loại có hệ thống các công trình nghiên cứu cơ học vi mô và tính toán đa tỉ lệ trên thế giới:

  • Phương pháp cận và đánh giá giải tích: Khởi đầu từ các nghiên cứu tính toán biên cứng của Voigt [1] và Reuss [2] theo tỉ lệ thể tích; Hashin và Shtrikman [3] phát triển phương pháp biến phân bậc hai; Washizu [9] với nguyên lý biến phân tổng quát cho bài toán phi tuyến; Terada và Kikuchi [10] đề xuất mô phỏng hai cấp độ cho cấu trúc tuần hoàn. Các nghiên cứu mở rộng cho vật liệu đa tinh thể dị hướng [11–14]. Đối với kết cấu tấm chịu uốn, Kolpakov và cs [4, 5] đã hiệu chỉnh cận Voigt - Reuss; Nguyen và cs [15] hiệu chỉnh cận Hashin - Shtrikman. Nhược điểm của nhóm phương pháp này là biên độ dao động giữa cận trên và cận dưới rất lớn khi độ chênh lệch cơ tính giữa các pha cao và không phản ánh được hình học chi tiết của pha vi mô.
  • Phương pháp xấp xỉ chuỗi và biến đổi Fourier: Walker và cs [16] dùng xấp xỉ chuỗi Fourier phân tích ứng xử đàn nhớt dẻo; Fotiu và Nemat [17] áp dụng cho vật liệu tuần hoàn 2D và 3D. Moulinec và Suquet [6, 18] đề xuất phương pháp biến đổi Fourier nhanh (FFT) trên ảnh vi cấu trúc không cần chia lưới. Bary [19, 20] so sánh FEM và FFT trên mô hình bê tông 3D. Eyre và Milton [21], Michel và cs [22] phát triển thuật toán Lagrangian bổ sung để giải quyết sự mất hội tụ khi vật liệu có lỗ rỗng; mở rộng cho kết cấu tấm bởi [23].
  • Phương pháp phần tử hữu hạn miền Voronoi (VCFEM): Ghosh và cs [24, 25] chia miền cấu trúc thành các đa giác Voronoi để xử lý vật liệu hai pha trong miền đàn hồi tuyến tính.
  • Phương pháp phần tử hữu hạn đa tỉ lệ ($\text{FE}^2$): Smit [26], Feyel [8, 27], Miehe và cs [22, 28], Renard và cs [29] phát triển mô hình lặp phi tuyến hai cấp độ. Điểm hạn chế căn bản của $\text{FE}^2$ là bài toán vi mô phải giải tại từng điểm tích phân Gauss của lưới phần tử vĩ mô trong mỗi bước lặp Newton - Raphson, dẫn đến chi phí lưu trữ bộ nhớ và thời gian tính toán quá lớn.
  • Phương pháp đa tỉ lệ ngoài miền đàn hồi và phân tích giới hạn: Meyer [30] đặt nền tảng toán học; Zohdi và Wriggers [31], Ladevèze và Fish [32], Ma và cs [33, 34], Sadowski [35], Li và cs [36], Galvanetto và cs [37] nghiên cứu quan hệ giữa các cấp độ. Về mặt dẻo vĩ mô, các nghiên cứu kinh điển của Hodge [38], Save và Massonnet [39], Zyczkowski [40], Sawczuk [41], Liu và Jiang [42], Yu và cs [43, 44] khảo sát mặt dẻo tương tác và lý thuyết dẻo tuyến tính hóa từng đoạn; Bigoni [45] hiệu chỉnh tiêu chuẩn cho vật liệu giòn. Về cấp vi mô, Kolupaev [46] tổng kết các dạng hàm dẻo; Gurson [25] xây dựng mô hình cho vật liệu có lỗ rỗng dựa trên von Mises; Li và cs [47–50] kết hợp lý thuyết giới hạn động học và đồng nhất hóa cho kim loại gia cường sợi; Schuh và Lund [55, 56] khảo sát chảy dẻo kim loại thủy tinh; ứng dụng trên kết cấu gạch đá xây [51–54].

Khoảng trống luận án lựa chọn giải quyết

Luận án tập trung giải quyết khoảng trống tính toán bằng cách không sử dụng mô hình lặp $\text{FE}^2$ tốn kém, mà xây dựng giải pháp tính toán trực tiếp:

  1. Thiết lập quy trình đồng nhất hóa đàn hồi giải tích - số cho ba dạng phần tử đại diện (tấm phẳng 2D, khối 3D và tấm uốn RPE) có tích hợp kỹ thuật nén bậc tự do cho điều kiện biên tuần hoàn.
  2. Kết hợp lý thuyết phân tích giới hạn động học cận trên với kỹ thuật tối ưu hóa nón bậc hai (SOCP) trên phần tử đại diện vi mô để xác định trực tiếp các điểm ứng suất giới hạn vĩ mô, từ đó dùng phương pháp bình phương cực tiểu thiết lập hàm mặt dẻo hữu hiệu tổng quát (Hill và Tsai-Wu) cho kết cấu tấm không đồng nhất.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết và khái niệm cốt lõi

  • Mô hình vật liệu: Sử dụng mô hình vật liệu cứng dẻo lý tưởng và đàn dẻo lý tưởng. Biến dạng phân tách thành phần đàn hồi và phần dẻo ($d\varepsilon_{ij} = d\varepsilon^e_{ij} + d\varepsilon^p_{ij}$). Định luật Hooke biểu diễn mối quan hệ đàn hồi ($\varepsilon^e_{ij} = H_{ijkl}\sigma_{kl}$). Hệ thống tiêu chuẩn chảy dẻo được khảo sát gồm: Tresca (ứng suất tiếp cực đại), von Mises ($J_2$), Hill (cho vật liệu trực hướng, đối xứng kéo - nén), Mohr - Coulomb và Drucker - Prager (cho vật liệu có góc nội ma sát $\varphi$ và độ dính $c$), và Tsai-Wu (dạng tổng quát bậc hai cho vật liệu bất đẳng hướng, bất đối xứng kéo - nén dạng $\mathbf{\sigma}^T \mathbf{P} \mathbf{\sigma} + \mathbf{\sigma}^T \mathbf{Q} - 1 = 0$).
  • Định đề Drucker và Luật chảy dẻo kết hợp: Đảm bảo tính ổn định vật liệu qua bất đẳng thức Drucker $\oint (\sigma_{ij} - \sigma^0_{ij})d\varepsilon_{ij} \ge 0$. Luật chảy dẻo kết hợp xác định vectơ vận tốc biến dạng dẻo vuông góc với mặt thế dẻo ($\dot{\varepsilon}^p_{ij} = \dot{\lambda} \frac{\partial f}{\partial \sigma_{ij}}, \dot{\lambda} \ge 0$).
  • Lý thuyết đa tỉ lệ và định lý trung bình thể tích: Liên kết trường vĩ mô ($\mathbf{\varepsilon}_M, \mathbf{\sigma}_M$) và trường vi mô ($\mathbf{\varepsilon}_m, \mathbf{\sigma}_m$) trên thể tích phần tử đại diện $\Omega_m$ thông qua định lý trung bình: $$\mathbf{\varepsilon}M = \frac{1}{\Omega_m}\int{\Omega_m}\mathbf{\varepsilon}_m d\Omega_m, \quad \mathbf{\sigma}M = \frac{1}{\Omega_m}\int{\Omega_m}\mathbf{\sigma}m d\Omega_m = \frac{1}{\Omega_m}\sum{i=1}^{N_p} \mathbf{f}_i \mathbf{X}_i$$
  • Lý thuyết phân tích giới hạn cận trên: Đánh giá hệ số tải trọng giới hạn $\lambda^+$ qua nguyên lý công khả dĩ động và năng lượng tiêu tán dẻo nội năng ($W_I = \int_{\Omega} D(\dot{\mathbf{\varepsilon}}^p) d\Omega$) cân bằng với công ngoại lực ($W_E$).
  • Lý thuyết tấm mỏng Kirchhoff (CPT): Giả thiết đường pháp tuyến thẳng và vuông góc với mặt trung bình sau biến dạng. Biến dạng $\mathbf{\varepsilon} = z\mathbf{\kappa}$, trong đó $\mathbf{\kappa} = [-\frac{\partial^2 w}{\partial x^2}, -\frac{\partial^2 w}{\partial y^2}, -2\frac{\partial^2 w}{\partial x \partial y}]^T$.

Phương pháp nghiên cứu và kỹ thuật phân tích

  • Phương pháp phần tử hữu hạn (FEM): Rời rạc hóa trường chuyển vị vi mô bằng các phần tử tam giác 3 nút (T3), phần tử tứ giác 4 nút (Q4), phần tử khối 8 nút (H8) và phần tử tấm HCT tương thích bậc $C^1$ có 12 bậc tự do (độ võng $w$, hai góc xoay $\theta_x, \theta_y$ tại 3 đỉnh và 3 đạo hàm chuẩn tại trung điểm 3 cạnh).
  • Kỹ thuật nén bậc tự do (Condensation): Phân chia hệ bậc tự do thành nhóm độc lập ($u_i$) và phụ thuộc ($u_d$) theo quan hệ biên tuần hoàn $\mathbf{u}_d = \mathbf{C}\mathbf{u}_i$, từ đó khử các bậc tự do phụ thuộc, đảm bảo tính đối xứng và xác định dương của ma trận độ cứng tổng thể.
  • Phương pháp tối ưu hóa nón bậc hai (SOCP): Chuyển đổi hàm năng lượng tiêu tán dẻo phi tuyến và các ràng buộc tương thích, cân bằng, biên tuần hoàn về dạng nón chuẩn $K_q = {x \in \mathbb{R}^n: x_1 \ge \sqrt{\sum_{i=2}^n x_i^2}, x_1 \ge 0}$ và giải bằng bộ giải MOSEK.
  • Phương pháp xấp xỉ bình phương cực tiểu (LSA): Khớp tập hợp các điểm ứng suất giới hạn thu được từ cấp vi mô để xác định các hệ số của ma trận $\mathbf{P}$ và véctơ $\mathbf{Q}$ của mặt dẻo vĩ mô hữu hiệu.

Nội dung chính theo từng chương

Chương 1: Tổng quan

Chương 1 trình bày bối cảnh nghiên cứu cơ học đa tỉ lệ, sự cần thiết của việc xác định thuộc tính hữu hiệu cho vật liệu không đồng nhất cấu tạo nên kết cấu tấm. Tác giả tổng thuật chi tiết 4 nhóm phương pháp trong miền đàn hồi (phương pháp cận, phương pháp chuỗi/FFT, phương pháp miền Voronoi VCFEM, phương pháp $\text{FE}^2$) và các nghiên cứu xác định tiêu chuẩn dẻo vĩ mô - vi mô ngoài miền đàn hồi. Xác định rõ mục tiêu, đối tượng, phạm vi nghiên cứu và cấu trúc 9 chương của luận án.

Chương 2: Lý thuyết nền tảng

Chương 2 thiết lập toàn bộ cơ sở giải tích phục vụ các chương tính toán tiếp theo. Nội dung bao gồm mô hình vật liệu cứng dẻo/đàn dẻo lý tưởng, định đề Drucker, luật chảy dẻo kết hợp và dạng toán học của 6 tiêu chuẩn dẻo (Tresca, von Mises, Hill, Mohr-Coulomb, Drucker-Prager, Tsai-Wu). Trình bày định lý trung bình thể tích, nguyên lý Hill-Mandel cho RVE. Xây dựng lý thuyết phân tích giới hạn cận trên, thiết lập công thức năng lượng tiêu tán dẻo cho vật liệu Hill và Tsai-Wu, định nghĩa bài toán tối ưu hóa hình nón bậc hai (SOCP). Cuối cùng là lý thuyết tấm Kirchhoff và phần tử tấm tương thích $C^1$ HCT 12 bậc tự do.

Chương 3: Bài toán đa tỉ lệ đàn hồi với phần tử đại diện tấm phẳng hai chiều

Chương 3 xây dựng mô hình đồng nhất hóa tấm phẳng 2D chịu trạng thái ứng suất phẳng dưới tác dụng của lực màng. Thiết lập chi tiết điều kiện biên tuyến tính và điều kiện biên tuần hoàn trên 4 cạnh của RVE, giải quyết quan hệ chuyển vị các nút đối xứng và 4 nút góc thông qua kỹ thuật nén bậc tự do. Phân tích 3 mode biến dạng cơ bản ($\varepsilon_{xx}, \varepsilon_{yy}, \gamma_{xy}$) để xác định ma trận độ cứng hữu hiệu $\mathbf{D}M$ và các hằng số đàn hồi hữu hiệu ($E{eff}, \nu_{eff}, G_{eff}, K_{eff}$).

Chương này thực hiện 5 ví dụ số:

  1. Vật liệu có cốt sợi hình chữ nhật: So sánh mô hình cốt sợi ngắn và cốt sợi dài dưới hai điều kiện biên.
  2. Vật liệu có cốt sợi hình tròn: Khảo sát sự thay đổi của mô đun trượt $G_{eff}$ và ma trận vật liệu khi tỉ lệ thể tích sợi $V_f$ biến thiên từ 0.1 đến 0.5.
  3. Vật liệu có lỗ rỗng: Khảo sát ảnh hưởng của thể tích lỗ rỗng tròn đến độ suy giảm các mô đun đàn hồi.
  4. Vật liệu cơ lý biến thiên (FGM): Khảo sát 2 mẫu vi cấu trúc A và B với sự biến thiên liên tục của các thành phần ma trận độ cứng.
  5. Vật liệu đa tinh thể dị hướng: Mô phỏng đa tinh thể nhôm (Al) với sự phân bố ngẫu nhiên góc định hướng hạt $[0^\circ; 90^\circ]$ để xác định mô đun trượt hữu hiệu.

Chương 4: Bài toán đa tỉ lệ đàn hồi với phần tử đại diện 3D

Chương 4 mở rộng mô hình đồng nhất hóa cho phần tử thể tích đại diện không gian 3D sử dụng phần tử khối 8 nút (H8). Thiết lập điều kiện biên tuần hoàn trên 6 mặt ngoài của khối lập phương RVE thông qua các cặp mặt đối diện và 8 đỉnh góc. Kỹ thuật nén bậc tự do được mở rộng để loại bỏ các bậc tự do phụ thuộc trong không gian 3 chiều.

Các bài toán kiểm chứng bao gồm:

  • Mô hình RVE phân bố vật liệu 3 chiều dạng nhiều lớp ngang, nhiều lớp đứng và phân bố xen kẽ; đánh giá sự suy giảm của mô đun dọc trục $D_{11}$ và mô đun khối $K_{eff}$ khi bề dày các lớp giảm dần.
  • Kết cấu tấm composite nhiều lớp chịu uốn: Khảo sát mẫu RVE tấm hướng sợi chữ thập (cross-ply), hướng sợi xiên (angle-ply) và tấm laminate 4 lớp lưới phần tử $12 \times 12 \times 12$ H8.
  • Kết cấu tấm FGM 3D cấu tạo từ kim loại nhôm ($\text{Al}_2\text{O}_3$) và gốm ($\text{ZrO}2$) cho mẫu 1 lớp và mẫu 3 lớp, đánh giá phân bố thống kê mô đun kéo $D{11}$.

Chương 5: Bài toán đa tỉ lệ đàn hồi với phần tử tấm vi mô chịu uốn

Chương 5 phát triển mô hình phần tử tấm đại diện chịu uốn (Representative Plate Element - RPE) dựa trên lý thuyết tấm mỏng Kirchhoff. Mối liên hệ giữa trường biến dạng cong vĩ mô ($\kappa_M = [\kappa_{xx}, \kappa_{yy}, \kappa_{xy}]^T$) và trường chuyển vị góc xoay/độ võng trên biên RPE được thiết lập dưới dạng bậc hai. Sử dụng phần tử tấm HCT 12 bậc tự do để mô hình hóa kết cấu vi mô chịu uốn nhằm đảm bảo tính liên tục $C^1$ của đạo hàm chuyển vị qua các biên phần tử.

Kỹ thuật đồng nhất hóa uốn tính toán trực tiếp ma trận độ cứng uốn hữu hiệu $\mathbf{D}b$ và các thông số $E{eff}, \nu_{eff}, K_{eff}, G_{eff}$. Các ví dụ số gồm:

  • Tấm vi mô khoét lỗ hình vuông với tỉ lệ thể tích lỗ thay đổi, phân tích các mode chuyển vị uốn tương ứng với biến dạng cong $\kappa_{xx}, \kappa_{yy}, \kappa_{xy}$.
  • Tấm mỏng 3 lớp có lỗ tròn tuần hoàn chịu uốn: So sánh trường chuyển vị, độ võng không thứ nguyên tại tâm tấm giữa phương pháp đa tỉ lệ đồng nhất hóa và phương pháp đơn tỉ lệ mô hình hóa chi tiết.

Chương 6: Miền chảy dẻo cho vật liệu theo tiêu chuẩn Hill

Chương 6 tập trung vào bài toán đa tỉ lệ ngoài miền đàn hồi nhằm xác định miền cường độ giới hạn và mặt dẻo hữu hiệu cho vật liệu không đồng nhất đối xứng (cường độ chịu kéo bằng cường độ chịu nén theo từng phương). Áp dụng lý thuyết phân tích giới hạn động học cận trên trên RVE 2D, năng lượng tiêu tán dẻo của vật liệu nền tuân theo tiêu chuẩn Hill trực hướng được chuyển đổi sang hàm theo trường vận tốc biến dạng dẻo $\dot{\mathbf{\varepsilon}}^p$ thông qua ma trận $\mathbf{\Theta} = \mathbf{P}^{-1}$ và khai triển Cholesky $\mathbf{S}$.

Toàn bộ bài toán phân tích giới hạn vi mô được biểu diễn dưới dạng bài toán tối ưu hóa hình nón bậc hai (SOCP) với các biến nón và biến ngoài nón. Sử dụng kỹ thuật bình phương cực tiểu, các hệ số hữu hiệu của mặt dẻo Hill vĩ mô ($P_{11}, P_{22}, P_{12}, P_{33}$) được xác định.

Các trường hợp kiểm chứng:

  1. Tấm kim loại có 1 lỗ hình tròn và 1 lỗ hình chữ nhật chịu kéo/nén: Khảo sát sự thay đổi cường độ dọc trục theo góc xoay $\alpha$ và tỉ lệ kích thước bán kính ($r/a$).
  2. Tấm kim loại gia cường cốt sợi dưới tải trọng dọc trục hai phương.
  3. Vật liệu khoét 2 lỗ rỗng ở các vị trí tương đối khác nhau.
  4. Tấm thép khoét nhiều lỗ phân bố đều tuần hoàn ($l/a = 0.5$) dưới góc tải $\alpha$ biến thiên; xây dựng mặt chảy dẻo 3D trong không gian ứng suất phẳng $(\Sigma_{11}, \Sigma_{22}, \Sigma_{12})$ và biểu đồ phân bố năng lượng tiêu tán dẻo.

Chương 7: Miền chảy dẻo cho vật liệu theo tiêu chuẩn Tsai-Wu

Chương 7 mở rộng bài toán phân tích giới hạn động học cho vật liệu không đồng nhất bất đối xứng (cường độ chịu kéo khác cường độ chịu nén theo mỗi phương) tuân theo tiêu chuẩn dẻo tổng quát Tsai-Wu ($\mathbf{\sigma}^T \mathbf{P} \mathbf{\sigma} + \mathbf{\sigma}^T \mathbf{Q} - 1 = 0$). Thiết lập công thức giải tích cho năng lượng tiêu tán dẻo $D(\dot{\mathbf{\varepsilon}}^p)$ phụ thuộc đồng thời vào ma trận $\mathbf{P}$ và véctơ $\mathbf{Q}$.

Khai triển bài toán phân tích giới hạn cận trên trên RVE vi mô thành bài toán tối ưu hóa ràng buộc nón bậc hai (SOCP) giải bằng MOSEK. Xấp xỉ mặt dẻo hữu hiệu Tsai-Wu ở cấp vĩ mô bằng phương pháp bình phương cực tiểu.

Các ví dụ số:

  1. Vật liệu hỗn hợp gia cường cốt sợi tròn: Khảo sát ảnh hưởng của tỉ lệ thể tích sợi $V_f$ đến hình dạng và sự dịch chuyển tâm của mặt dẻo hữu hiệu 3D.
  2. Vật liệu có lỗ rỗng tròn tuần hoàn ($V_f = 0.1, 0.2, 0.3$): So sánh mặt dẻo hữu hiệu khi vật liệu nền tuân theo von Mises, Hill và Tsai-Wu.
  3. Vật liệu có 16 lỗ rỗng phân bố ngẫu nhiên ($V_f = 0.1, 0.2$): Phân tích trường năng lượng tiêu tán dẻo và tính toán mặt dẻo vĩ mô giới hạn.

Chương 8: Thảo luận

Chương 8 tổng hợp, phân tích và so sánh có hệ thống ưu điểm và hạn chế của từng phương pháp đề xuất trong luận án:

  • Bài toán đa tỉ lệ tấm phẳng 2D đàn hồi: Ưu điểm là thuật toán đơn giản, tốc độ tính toán nhanh, điều kiện biên tuần hoàn nén bậc tự do cho kết quả chính xác; hạn chế là chưa xét được hiệu ứng uốn và hiệu ứng ứng suất theo phương chiều dày.
  • Bài toán đa tỉ lệ kết cấu 3D đàn hồi: Ưu điểm là mô tả đầy đủ trạng thái ứng suất 3 chiều của các lớp vật liệu phức tạp; hạn chế là số bậc tự do lớn, thời gian lập ma trận độ cứng cao.
  • Bài toán đa tỉ lệ tấm chịu uốn: Ưu điểm là tích hợp được phần tử HCT $C^1$ giúp xác định chính xác độ cứng uốn $\mathbf{D}_b$; hạn chế là chỉ áp dụng cho giả thuyết tấm mỏng Kirchhoff, chưa xét biến dạng trượt ngang của tấm dày (Mindlin - Reissner).
  • Mô hình mặt dẻo Hill và Tsai-Wu ngoài đàn hồi: Ưu điểm là giải trực tiếp tải trọng giới hạn bằng SOCP trong thời gian rất ngắn (vài giây), không cần lặp gia tải từng bước; hạn chế là giả thuyết vật liệu cứng dẻo lý tưởng không kể đến hiện tượng tái bền và hóa mềm phá hủy.

Chương 9: Kết luận và kiến nghị

Tổng kết toàn diện các kết quả đạt được về mặt phương pháp luận và thuật toán số cho cả hai nhóm bài toán đàn hồi và thiết kế dẻo. Nêu các kiến nghị khoa học nhằm phát triển nghiên cứu trong tương lai.

Tổng hợp dữ liệu số và bảng số liệu tiêu biểu

Dưới đây là bảng tổng hợp các kết quả số và thông số vật liệu hữu hiệu tiêu biểu được trích xuất từ các chương của luận án:

Hạng mục nghiên cứu Dạng cấu trúc vi mô / RVE Thông số đầu vào của các pha Kết quả hữu hiệu vĩ mô xác định được
Đàn hồi tấm phẳng 2D (Chương 3) Vật liệu gia cường cốt sợi tròn tuần hoàn Vật liệu nền và sợi đẳng hướng; Tỉ lệ thể tích sợi $V_f = 0.1 \div 0.5$ Xác định ma trận hằng số vật liệu $\mathbf{D}M$, các mô đun $E{eff}, \nu_{eff}, G_{eff}, K_{eff}$; điều kiện biên tuần hoàn cho kết quả nằm trong biên hẹp
Đàn hồi tấm phẳng 2D (Chương 3) Vật liệu cơ lý biến thiên (FGM) Mẫu A và B Kim loại ($\text{Al}_2\text{O}_3$) và Gốm ($\text{ZrO}_2$) biến thiên tọa độ Phân bố không gian của các thành phần ma trận hằng số vật liệu và thông số đàn hồi hữu hiệu
Đàn hồi tấm phẳng 2D (Chương 3) Đa tinh thể nhôm (Al) Hạt tinh thể dị hướng, phân bố góc ngẫu nhiên $[0^\circ; 90^\circ]$ Phân bố thống kê của mô đun trượt hữu hiệu $G_{eff}$ của đa tinh thể nhôm
Đàn hồi kết cấu 3D (Chương 4) Tấm 3D nhiều lớp (ngang, đứng, xen kẽ, cross-ply, angle-ply) Các lớp vật liệu trực hướng; Lưới phần tử H8 ($12 \times 12 \times 12$) Ma trận độ cứng 3D; quan hệ suy giảm của mô đun dọc trục $D_{11}$ và mô đun khối $K_{eff}$ theo chiều dày lớp
Đàn hồi tấm chịu uốn (Chương 5) Tấm vi mô khoét lỗ vuông và lỗ tròn (RPE) Phần tử tấm HCT 12 bậc tự do; Tỉ lệ lỗ rỗng $V_f$ biến thiên Ma trận độ cứng uốn $\mathbf{D}_b$; độ võng chuẩn hóa tại tâm tấm biên tựa chu vi khớp với nghiệm mô hình hóa chi tiết
Thiết kế dẻo Hill (Chương 6) Tấm có lỗ tròn, lỗ chữ nhật, 2 lỗ, nhiều lỗ ($l/a = 0.5$) Nền kim loại dẻo theo tiêu chuẩn Hill đối xứng kéo - nén Các hệ số mặt dẻo hữu hiệu $P_{11}, P_{22}, P_{12}, P_{33}$; mặt dẻo 3D $(\Sigma_{11}, \Sigma_{22}, \Sigma_{12})$; trường năng lượng tiêu tán dẻo
Thiết kế dẻo Tsai-Wu (Chương 7) Cốt sợi tròn, lỗ tròn tuần hoàn, 16 lỗ ngẫu nhiên Tiêu chuẩn Tsai-Wu bất đối xứng ($\Sigma_{Yt} \ne \Sigma_{Yc}$); $V_f = 0.1, 0.2$ Xác định ma trận $\mathbf{P}$ và véctơ $\mathbf{Q}$; xây dựng mặt dẻo 3D đóng kín dịch trục; trường tiêu tán năng lượng dẻo

Kết quả và những đóng góp mới

Luận án đã đạt được những đóng góp khoa học và kỹ thuật cụ thể sau:

  1. Về mặt lý luận và phương pháp:

    • Xây dựng thành công khung phân tích đa tỉ lệ đồng nhất hóa cho kết cấu tấm không đồng nhất ở 3 dạng mô hình hình học: tấm phẳng 2D, khối 3D và tấm mỏng chịu uốn RPE.
    • Ứng dụng thành công phần tử tấm tương thích $C^1$ HCT 12 bậc tự do trong bài toán đồng nhất hóa tấm mỏng chịu uốn Kirchhoff, đảm bảo tính liên tục của độ võng và góc xoay trên biên vi mô.
    • Phát triển hoàn chỉnh kỹ thuật nén bậc tự do (condensation) để khử các bậc tự do phụ thuộc trong điều kiện biên tuần hoàn, giúp ma trận độ cứng đối xứng và giảm đáng kể kích thước hệ phương trình đại số.
    • Thiết lập quy trình phân tích giới hạn động học cận trên trên RVE vi mô kết hợp tối ưu hóa hình nón bậc hai (SOCP) cho vật liệu bất đẳng hướng tuân theo tiêu chuẩn Hill và Tsai-Wu, cho phép tìm trực tiếp miền ứng suất giới hạn vĩ mô.

  2. Về mặt ứng dụng và thực tiễn:

    • Cung cấp phương pháp tính toán hiệu quả, giảm thiểu khối lượng bộ nhớ và thời gian tính toán so với phương pháp $\text{FE}^2$ truyền thống, đồng thời vượt qua các hạn chế dao động biên rộng của phương pháp cận Voigt - Reuss và Hashin - Shtrikman.
    • Cung cấp các công thức, bảng thông số đàn hồi hữu hiệu và phương trình mặt dẻo hữu hiệu cho nhiều hệ vật liệu kỹ thuật thực tế (composite sợi, tấm FGM kim loại - gốm, đa tinh thể nhôm, tấm đục lỗ).
    • Đề xuất thuật toán tối ưu hóa nón bậc hai giải bằng MOSEK đạt tốc độ hội tụ cao và ổn định số cho các bài toán cơ học phi tuyến phức tạp.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của luận án

  • Về mô hình vật liệu: Phân tích giới hạn dẻo chỉ áp dụng cho mô hình vật liệu cứng dẻo lý tưởng, chưa xét đến hiện tượng tái bền (hardening) hoặc làm mềm (softening) của vật liệu trong quá trình biến dạng dẻo phát triển.
  • Về hình học và động học: Giả thiết chuyển vị và biến dạng bé; bài toán tấm chịu uốn mới chỉ dừng lại ở lý thuyết tấm mỏng cổ điển Kirchhoff, chưa xét đến biến dạng trượt ngang của tấm dày theo lý thuyết Mindlin - Reissner.
  • Về tương tác vi cấu trúc: Xem sự liên kết giữa các pha vật liệu (như giữa sợi và nền) là liên kết hoàn hảo tuyệt đối, chưa xét đến hiện tượng trượt, nứt tách lớp (debonding) hoặc nứt gãy tế vi cục bộ.

Hướng nghiên cứu tiếp theo

  • Mở rộng phương pháp đa tỉ lệ cho bài toán biến dạng lớn và phi tuyến hình học.
  • Phát triển bài toán đồng nhất hóa tấm chịu uốn cho tấm dày kể đến biến dạng trượt ngang bậc nhất (FSDT) hoặc bậc cao (HSDT).
  • Tích hợp mô hình phá hủy liên kết pha, cơ học phá hủy vật liệu (damage mechanics) và bài toán đa trường vật lý (nhiệt - cơ, áp điện) vào khung tính toán đa tỉ lệ.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các đối tượng và lĩnh vực sau:

  • Nghiên cứu sinh và học viên cao học: Tham khảo phương pháp thiết lập điều kiện biên tuần hoàn, kỹ thuật nén bậc tự do trên RVE/RPE và kỹ thuật lập trình phần tử hữu hạn tấm $C^1$ HCT.
  • Giảng viên và nhà nghiên cứu Cơ học tính toán: Sử dụng khung lý thuyết phân tích giới hạn động học cận trên chuyển đổi sang bài toán tối ưu hóa nón bậc hai (SOCP) để giải quyết các bài toán tải trọng giới hạn trên nền tảng phần mềm MOSEK.
  • Kỹ sư thiết kế kết cấu công trình và vật liệu mới: Ứng dụng các bảng thông số đàn hồi hữu hiệu và phương trình mặt dẻo Hill, Tsai-Wu của vật liệu composite, FGM và kết cấu tấm đục lỗ để tính toán nhanh khả năng chịu lực ở cấp độ kết cấu vĩ mô mà không cần mô hình hóa chi tiết vi cấu trúc.

Câu hỏi thường gặp

1. Luận án giải quyết bài toán đa tỉ lệ trong miền đàn hồi cho những dạng phần tử đại diện nào?

Luận án thực hiện đồng nhất hóa đàn hồi cho 3 dạng phần tử đại diện:

  1. Phần tử đại diện tấm phẳng hai chiều (RVE 2D) chịu trạng thái ứng suất phẳng dưới tác dụng của lực màng (Chương 3).
  2. Phần tử thể tích đại diện không gian ba chiều (RVE 3D) sử dụng phần tử khối H8 cho vật liệu nhiều lớp và FGM (Chương 4).
  3. Phần tử tấm đại diện chịu uốn (RPE) dựa trên lý thuyết tấm mỏng Kirchhoff và phần tử tấm tương thích $C^1$ HCT 12 bậc tự do (Chương 5).

2. Kỹ thuật nén bậc tự do (condensation) đóng vai trò gì trong việc xử lý điều kiện biên tuần hoàn?

Trong điều kiện biên tuần hoàn, chuyển vị của các nút trên các cạnh/mặt đối xứng của RVE có mối liên hệ phụ thuộc lẫn nhau ($\mathbf{u}_d = \mathbf{C}\mathbf{u}_i$). Kỹ thuật nén bậc tự do được sử dụng để khử hoàn toàn các bậc tự do phụ thuộc này ra khỏi hệ phương trình tổng thể, giúp ma trận độ cứng sau khi nén giữ nguyên tính đối xứng, xác định dương, giảm kích thước hệ phương trình và đảm bảo độ ổn định số khi giải.

3. Phương pháp phân tích giới hạn ngoài miền đàn hồi trong luận án được giải quyết bằng thuật toán tối ưu hóa nào?

Luận án sử dụng phương pháp phân tích giới hạn động học cận trên. Năng lượng tiêu tán dẻo của cấu trúc vi mô được thiết lập thông qua luật chảy dẻo kết hợp và biểu diễn dưới dạng hàm theo biến dạng. Bài toán tìm hệ số tải trọng giới hạn nhỏ nhất được chuyển đổi về bài toán tối ưu hóa hình nón bậc hai (SOCP - Second-Order Cone Programming) và giải trực tiếp bằng công cụ MOSEK, giúp tìm nghiệm giới hạn nhanh chóng mà không cần lặp gia tải từng bước.

4. Sự khác biệt cơ bản giữa mặt dẻo hữu hiệu dạng Hill (Chương 6) và dạng Tsai-Wu (Chương 7) là gì?

  • Tiêu chuẩn Hill là dạng tổng quát hóa của tiêu chuẩn von Mises cho vật liệu trực hướng, áp dụng cho vật liệu có cường độ chịu kéo bằng cường độ chịu nén theo mỗi phương ($\Sigma_{Yt} = \Sigma_{Yc}$), do đó mặt dẻo nhận gốc tọa độ làm tâm đối xứng.
  • Tiêu chuẩn Tsai-Wu là tiêu chuẩn bậc hai tổng quát cho vật liệu bất đẳng hướng có khả năng chịu kéo khác khả năng chịu nén theo từng phương ($\Sigma_{Yt} \ne \Sigma_{Yc}$), phương trình chứa cả ma trận $\mathbf{P}$ (bậc hai) và véctơ $\mathbf{Q}$ (bậc nhất), dẫn đến mặt dẻo hữu hiệu 3D bị dịch chuyển tâm khỏi gốc tọa độ.

5. Những loại vật liệu không đồng nhất cụ thể nào đã được mô phỏng kiểm chứng trong luận án?

Luận án đã kiểm chứng trên 5 nhóm cấu trúc vật liệu:

  • Vật liệu composite nền kim loại gia cường cốt sợi (sợi tròn, sợi chữ nhật, sợi ngắn, sợi dài).
  • Vật liệu có lỗ rỗng (lỗ đơn hình tròn, lỗ hình chữ nhật, 2 lỗ rỗng, hệ nhiều lỗ phân bố đều tuần hoàn và hệ 16 lỗ phân bố ngẫu nhiên).
  • Vật liệu cơ lý biến thiên FGM (hệ hợp phần gốm $\text{ZrO}_2$ và kim loại $\text{Al}_2\text{O}_3$).
  • Vật liệu đa tinh thể kim loại (đa tinh thể nhôm Al định hướng ngẫu nhiên).
  • Kết cấu tấm composite nhiều lớp chịu uốn (lớp ngang, lớp đứng, lớp xen kẽ, cross-ply và angle-ply).

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Hoàng Phương đã xây dựng một hệ thống phương pháp tính toán đa tỉ lệ toàn diện cho kết cấu tấm không đồng nhất ở cả hai trạng thái đàn hồi tuyến tính và giới hạn dẻo. Bằng việc kết hợp chặt chẽ giữa kỹ thuật đồng nhất hóa, phương pháp phần tử hữu hạn tấm $C^1$ HCT và lý thuyết phân tích giới hạn tối ưu hóa nón bậc hai (SOCP), công trình đã giải quyết hiệu quả bài toán xác định các hằng số đàn hồi hữu hiệu và thiết lập phương trình mặt dẻo vĩ mô (Hill và Tsai-Wu). Các kết quả số trên nhiều hệ vật liệu phức tạp khẳng định tính chính xác và hiệu quả tính toán của phương pháp đề xuất, đóng góp cơ sở khoa học quan trọng cho lĩnh vực Cơ kỹ thuật và tính toán kết cấu composite.