Tổng quan luận án

Luận án tiến sĩ "Nghiên cứu các phương pháp đa tỉ lệ kết cấu tấm không đồng nhất", chuyên ngành Cơ kỹ thuật (mã số 62520101) do nghiên cứu sinh Nguyễn Hoàng Phương thực hiện tại Trường Đại học Sư phạm Kỹ thuật Thành phố Hồ Chí Minh 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 (bảo vệ năm 2021), tập trung giải quyết các bài toán tính toán cơ học cho kết cấu tấm có cấu trúc vi mô phức tạp bằng phương pháp đa tỉ lệ trong cả hai miền ứng xử: đàn hồi tuyến tính và phi tuyến ngoài miền đàn hồi (phân tích giới hạn).

Tính cấp thiết của đề tài bắt nguồn từ sự phát triển của các loại vật liệu mới có cấu trúc không đồng nhất như vật liệu composite gia cường sợi, vật liệu có lỗ rỗng, vật liệu cơ tính biến thiên (Functionally Graded Materials - FGM) và vật liệu đa tinh thể. Khi phân tích cơ học các kết cấu này, việc mô hình hóa trực tiếp toàn bộ chi tiết cấu trúc vi mô bằng phương pháp phần tử hữu hạn truyền thống gặp trở ngại lớn do số lượng bậc tự do quá lớn, dẫn đến hiện tượng quá tải bộ nhớ và chi phí tính toán cao. Ngược lại, các lý thuyết giải tích cổ điển như lý thuyết tấm nhiều lớp (laminate plate theory) hay quy tắc pha trộn (rule of mixture) gặp khó khăn khi mô tả sự phức tạp về hình học và phân bố pha vi mô theo chiều dày tấm. Khoảng trống nghiên cứu được luận án xác định là sự cần thiết phải phát triển một khung tính toán đa tỉ lệ thống nhất, có khả năng chuyển đổi dữ liệu biến dạng từ cấp độ vĩ mô xuống phần tử đại diện vi mô, áp đặt chính xác các điều kiện biên động học, và đồng nhất hóa để tìm ra ma trận hằng số đàn hồi hữu hiệu cũng như miền ứng suất giới hạn (hàm chảy dẻo hữu hiệu) cho kết cấu tấm.

Mục tiêu nghiên cứu của luận án bao gồm hai nội dung chính:

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

Đối tượng nghiên cứu là kết cấu tấm không đồng nhất với các dạng cấu trúc vi mô: composite cốt sợi (ngắn, dài, tròn), vật liệu có lỗ rỗng (lỗ tròn, vuông, lỗ kép, nhiều lỗ tuần hoàn và ngẫu nhiên), vật liệu FGM kim loại - gốm ($Al_2O_3$ và $ZrO_2$), và vật liệu đa tinh thể nhôm ($Al$).

Phạm vi nghiên cứu của luận án được giới hạn:

  • Về không gian hình học: Xét phần tử thể tích đại diện (RVE) hai chiều cho tấm phẳng chịu kéo/nén trong mặt phẳng, phần tử đại diện không gian 3D, và phần tử tấm đại diện (RPE) cho tấm mỏng chịu uốn tuân theo giả thuyết Kirchhoff.
  • Về ứng xử vật liệu: Trong miền đàn hồi, các pha thành phần vi mô được giả định là đồng nhất và đẳng hướng, ứng xử theo định luật Hooke; ngoài miền đàn hồi, vật liệu nền vi mô được giả thiết tuân theo mô hình cứng dẻo lý tưởng (von Mises, Hill hoặc Tsai-Wu) với giả thiết biến dạng bé và áp dụng luật chảy dẻo kết hợp.

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

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

Trong miền đàn hồi:

  • Phương pháp cận biến phân: Bắt đầu từ các nghiên cứu biên cứng của Voigt (1889) [1] và Reuss (1929) [2] dựa trên tỉ lệ thể tích pha, sau đó được mở rộng bởi Hashin và Shtrikman (1963) [3] dựa trên công thức biến phân. Phương pháp biến phân tổng quát của Washizu (1968) [9] được Terada và Kikuchi (2001) [10] áp dụng cho bài toán mô phỏng hai cấp độ với vi cấu trúc tuần hoàn. Các nghiên cứu xác định mô đun đàn hồi cho vật liệu đa tinh thể dị hướng được phát triển bởi các tác giả [11–14]. Ứng dụng cho kết cấu tấm được Kolpakov và cs. (1991, 1999) [4, 5] phát triển qua việc hiệu chỉnh cận Voigt - Reuss, và Nguyen và cs. (2008) [15] qua hiệu chỉnh cận Hashin - Shtrikman. Hạn chế của phương pháp cận là biên độ dao động giữa cận trên và cận dưới rất lớn khi hai pha vật liệu có độ chênh lệch cơ tính cao.
  • Phương pháp xấp xỉ chuỗi và biến đổi Fourier: Walker và cs. (1991) [16] dùng chuỗi Fourier giải bài toán composite nhớt dẻo tuần hoàn. Fotiu và Nemat-Nasser (1996) [17] ước lượng tính chất vật liệu composite tuần hoàn 3D và 2D. Moulinec và Suquet (1998) [6, 18] đề xuất phương pháp biến đổi Fourier nhanh (FFT) sử dụng trực tiếp ảnh vi cấu trúc, tránh việc chia lưới phần tử. Bary (2008) [19, 20] so sánh FFT và FEM cho bê tông. Eyre và Milton (1999) [21] cùng Michel và cs. (1999, 2001) [22, 23] cải thiện thuật toán FFT cho vật liệu có lỗ rỗng và tấm tuần hoàn.
  • Phương pháp phần tử Voronoi (VCFEM): Ghosh và Mukhopadhyay (1991, 1993) [24, 25] phát triển kỹ thuật chia lưới đa giác Voronoi để mô phỏng vật liệu hai pha dạng hạt nhân không đồng nhất.
  • Phương pháp đa tỉ lệ phần tử hữu hạn ($FE^2$): Bắt đầu từ ý tưởng của Renard và cs. (1995) [29], Feyel (1999, 2000) [8, 27] hoàn thiện mô hình $FE^2$ giải bằng thuật toán Newton-Raphson. Miehe và cs. (1999, 2002) [22, 28] cùng Smit (1998) [26] phát triển khung tính toán liên kết vi mô - vĩ mô, tuy nhiên phương pháp này đòi hỏi lưu trữ và lặp lại tại từng điểm Gauss ở cấp độ toàn cục.

Ngoài miền đàn hồi:

  • Tiêu chuẩn dẻo vĩ mô và tương tác: Các nghiên cứu của Hodge (1959) [38], Save và Massonnet (1972) [39], Zyczkowski (1981) [40], Sawczuk (1989) [41] xây dựng mặt dẻo tương tác cho các cấu kiện chịu lực phức hợp. Liu và Jiang (1998) [42] khảo sát sự phát triển vùng dẻo của tấm tròn. Yu và cs. (1998, 2006) [43, 44] áp dụng lý thuyết dẻo tuyến tính hóa từng đoạn cho kết cấu xây dựng và đất gia cường [51–54]. Bigoni (2000) [45] hiệu chỉnh tiêu chuẩn dẻo cho vật liệu giòn và đất. Kolupaev (2018) [46] tổng kết các hàm tiêu chuẩn dẻo.
  • Kỹ thuật phân tích giới hạn đa tỉ lệ: Mô hình Gurson (1977) [25] cho kim loại có lỗ rỗng nền von Mises. Li và cs. (2000–2006) [47–50] kết hợp lý thuyết phân tích giới hạn động học và đồng nhất hóa để tìm cường độ giới hạn cho composite kim loại gia cường sợi. Schuh và Lund (2003, 2004) [55, 56] nghiên cứu tiêu chuẩn dẻo cho kim loại thủy tinh.

Khoảng trống luận án lựa chọn giải quyết: Luận án tập trung xây dựng mô hình phần tử hữu hạn hai cấp độ kết hợp điều kiện biên tuần hoàn và tuyến tính thông qua kỹ thuật nén bậc tự do (condensation technique) cho các cấu trúc tấm (phẳng, 3D, tấm mỏng chịu uốn). Ở miền ngoài đàn hồi, luận án xây dựng bài toán phân tích giới hạn cận trên cho phần tử RVE dưới dạng bài toán tối ưu hóa hình nón bậc hai (Second-Order Cone Programming - SOCP), giải bằng thuật toán điểm trong của công cụ Mosek, từ đó xác định miền cường độ hữu hiệu dạng Hill và Tsai-Wu thông qua kỹ thuật xấp xỉ bình phương cực tiểu.

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

Cơ sở lý thuyết

  1. Mô hình vật liệu và tiêu chuẩn chảy dẻo:
  • Mô hình cứng dẻo lý tưởng von Mises (1913): Bỏ qua biến dạng đàn hồi ($d\epsilon_{ij} = d\epsilon_{ij}^p$), phân tích trực tiếp trạng thái sụp đổ dẻo.
  • Mô hình đàn dẻo lý tưởng Prandtl (1928): Tổng biến dạng gồm đàn hồi tuân theo định luật Hooke và biến dạng dẻo: $d\epsilon_{ij} = d\epsilon_{ij}^e + d\epsilon_{ij}^p$.
  • Định đề Drucker: Đảm bảo tính ổn định của vật liệu trong chu trình tăng tải và dỡ tải khép kín: $\oint (\sigma_{ij} - \sigma_{ij}^0) d\epsilon_{ij} \ge 0$.
  • Luật chảy dẻo kết hợp: Vectơ vận tốc biến dạng dẻo $\dot{\epsilon}{ij}^p$ trực giao với mặt chảy dẻo: $\dot{\epsilon}{ij}^p = \dot{\lambda} \frac{\partial f}{\partial \sigma_{ij}}$ với $\dot{\lambda} \ge 0$.
  • Các tiêu chuẩn dẻo áp dụng: Tiêu chuẩn von Mises, tiêu chuẩn Tresca, tiêu chuẩn Hill (cho vật liệu bất đẳng hướng đối xứng kéo - nén: $\Sigma_{Ytx} = \Sigma_{Yty} \neq \Sigma_{Ycx} = \Sigma_{Ycy}$), tiêu chuẩn Mohr-Coulomb, tiêu chuẩn Drucker-Prager, và tiêu chuẩn tổng quát Tsai-Wu cho vật liệu bất đẳng hướng bất đối xứng kéo - nén ($\Sigma_{Ytx} \neq \Sigma_{Yty} \neq \Sigma_{Ycx} \neq \Sigma_{Ycy}$) với phương trình $\Psi_T(\sigma) = \sigma^T P \sigma + \sigma^T Q - 1 = 0$.
  1. Lý thuyết đa tỉ lệ và định lý trung bình thể tích:
  • Phần tử thể tích đại diện (RVE - $\Omega_m$): Mô tả đầy đủ vi cấu trúc vật liệu nhưng kích thước vi mô $l_{\text{vi mô}} \ll l_{\text{vĩ mô}}$ cho phép bỏ qua lực thể tích trong phương trình cân bằng vi mô $\nabla \sigma_m = 0$.
  • Định lý trung bình thể tích (Hill 1963, Nemat-Nasser 1999): $$\epsilon_M = \frac{1}{V_m} \int_{\Omega_m} \epsilon_m dV_m$$ $$\sigma_M = \frac{1}{V_m} \int_{\Omega_m} \sigma_m dV_m = \frac{1}{V_m} \sum_{i=1}^{N_p} f_i X_i$$ Trong đó $f_i$ và $X_i$ lần lượt là lực nút và tọa độ của nút $i$ trên biên $\Gamma_m$.
  1. Lý thuyết phân tích giới hạn cận trên (Kinematic Limit Analysis):
  • Hệ số tải trọng giới hạn cận trên $\lambda^+$ được xác định qua bài toán tối ưu hóa cực tiểu hóa công suất tiêu tán dẻo nội năng với điều kiện chuẩn hóa công ngoại lực $W_{\text{ext}} = 1$: $$\lambda^+ = \min \int_{\Omega} \sigma_{ij} \epsilon_{ij} dV \quad \text{với} \quad W_{\text{ext}}(u) = 1, \quad \epsilon_{ij} = \frac{1}{2}\left(\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right)$$
  • Năng lượng tiêu tán dẻo cho vật liệu Tsai-Wu được thiết lập theo biến dạng dẻo: $$D(\epsilon^p) = \frac{1}{2\mu} (\epsilon^p)^T P^{-1} \epsilon^p - \frac{1}{2} (\epsilon^p)^T P^{-1} Q$$
  1. Lý thuyết tấm mỏng Kirchhoff và phần tử tấm chịu uốn Hsieh-Clough-Tocher (HCT):
  • Quan hệ biến dạng - chuyển vị uốn: $\epsilon = z \kappa$, với $\kappa = - \nabla^2 u = L u$.
  • Phần tử tam giác tương thích $C^1$ HCT gồm 3 tam giác con bên trong; 12 bậc tự do ban đầu (độ võng $w_i$, hai góc xoay $\theta_{xi}, \theta_{yi}$ tại 3 đỉnh và góc xoay pháp tuyến $\theta_k$ tại 3 trung điểm cạnh) được khử chuyển vị nút trọng tâm thông qua điều kiện liên tục, sau đó giả thiết độ dốc pháp tuyến thay đổi tuyến tính dọc cạnh để khử 3 bậc tự do tại trung điểm cạnh, rút gọn phần tử về 9 bậc tự do.

Phương pháp nghiên cứu và công cụ tính toán

  • Rời rạc hóa phần tử hữu hạn (FEM): Sử dụng các loại phần tử phần tử tam giác 3 nút (T3), tứ giác 4 nút (Q4), phần tử khối 8 nút (H8), phần tử tam giác làm trơn trên cạnh (ES-T3), và phần tử tấm HCT.
  • Điều kiện biên: Áp đặt điều kiện biên tuyến tính (LBC) và điều kiện biên tuần hoàn (PBC). Ràng buộc điều kiện biên tuần hoàn giữa các cặp nút đối xứng trên biên được khử bằng phương pháp nén bậc tự do (condensation technique).
  • Thuật toán tối ưu hóa: Khai triển phi tuyến về dạng quy hoạch nón bậc hai (SOCP) và giải bằng phần mềm tối ưu Mosek.
  • Thuật toán xấp xỉ mặt dẻo: Phương pháp bình phương cực tiểu (LSA) khớp các hệ số ma trận $P$ và $Q$ từ tập nghiệm ứng suất vĩ mô giới hạn.

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

Chương 1: Tổng quan

Chương 1 giới thiệu tính cấp thiết của bài toán tính toán cơ học đa tỉ lệ cho kết cấu tấm không đồng nhất. Tác giả tổng quan có hệ thống các công trình nghiên cứu trong nước và quốc tế về các phương pháp đa tỉ lệ trong miền đàn hồi (phương pháp cận Voigt, Reuss, Hashin-Shtrikman, phương pháp FFT, phương pháp miền Voronoi VCFEM, phương pháp đa tỉ lệ phần tử hữu hạn $FE^2$) và ngoài miền đàn hồi (các tiêu chuẩn dẻo vĩ mô, lý thuyết phân tích giới hạn động học kết hợp đồng nhất hóa). Từ đó, chương xác lập mục tiêu nghiên cứu, phạm vi khảo sát và cấu trúc tổng thể gồm 9 chương của luận án.

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

Chương 2 trình bày đầy đủ các cơ sở lý thuyết toán - cơ học làm nền tảng cho toàn bộ luận án, bao gồm:

  • Các mô hình ứng xử vật liệu: cứng dẻo lý tưởng von Mises, đàn dẻo lý tưởng Prandtl; định đề Drucker về tính ổn định của vật liệu; luật chảy dẻo kết hợp.
  • Biểu thức giải tích của các tiêu chuẩn chảy dẻo: von Mises, Tresca, Hill, Mohr-Coulomb, Drucker-Prager, Tsai-Wu.
  • Lý thuyết đa tỉ lệ, phần tử thể tích đại diện RVE, nguyên lý Hill-Mandel và định lý trung bình thể tích xác định quan hệ ứng suất - biến dạng giữa cấp vi mô và vĩ mô.
  • Lý thuyết phân tích giới hạn (định lý cận trên và cận dưới), năng lượng tiêu tán dẻo của vật liệu Hill và Tsai-Wu; định nghĩa và cấu trúc của bài toán tối ưu hóa hình nón bậc hai (SOCP) giải bằng công cụ Mosek.
  • Lý thuyết tấm mỏng Kirchhoff chịu uốn và công thức phần tử tấm tam giác tương thích $C^1$ HCT 12 bậc tự do rút gọn về 9 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 quy trình tính toán đa tỉ lệ cho kết cấu tấm phẳng 2D chịu tải trọng trong mặt phẳng. Các công thức áp đặt điều kiện biên tuyến tính (LBC) và tuần hoàn (PBC) cho RVE được thiết lập chi tiết. Ma trận độ cứng hữu hiệu được xác định qua kỹ thuật nén bậc tự do và trung bình thể tích. Tác giả thực hiện 5 ví dụ số khảo sát các vi cấu trúc khác nhau:

  1. Composite cốt sợi chữ nhật ngắn và dài: So sánh lưới phần tử T3 và Q4 dưới hai điều kiện biên.
  2. Composite cốt sợi hình tròn: Đánh giá ảnh hưởng của tỉ lệ thể tích sợi $V_f$ đến các hằng số đàn hồi hữu hiệu ($E_{eff}, \nu_{eff}, G_{eff}, K_{eff}$) sử dụng phần tử T3 và phần tử làm trơn cạnh ES-T3.
  3. Vật liệu có lỗ rỗng: Khảo sát sự suy giảm của ma trận hằng số đàn hồi khi tỉ lệ thể tích lỗ rỗng thay đổi.
  4. Vật liệu cơ tính biến thiên FGM (mẫu A và mẫu B): Đánh giá phân bố không gian của các thành phần ma trận hằng số vật liệu.
  5. Vật liệu đa tinh thể nhôm ($Al$): Phân tích phân bố ngẫu nhiên góc định hướng hạt $[0^\circ, 90^\circ]$ và xác định mô đun trượt hữu hiệu $G_{eff}$.
Loại cấu trúc vi mô RVE 2D Loại phần tử hữu hạn Điều kiện biên khảo sát Đại lượng đàn hồi hữu hiệu trích xuất
Cốt sợi ngắn / cốt sợi dài T3, Q4 Tuần hoàn (PBC), Tuyến tính (LBC) Mô đun đàn hồi hữu hiệu, ma trận hằng số vật liệu
Cốt sợi gia cường hình tròn T3, ES-T3 Tuần hoàn (PBC), Tuyến tính (LBC) Mô đun kháng cắt $G_{eff}$, hằng số vật liệu theo $V_f$
Vật liệu có lỗ rỗng tròn T3 Tuần hoàn (PBC), Tuyến tính (LBC) $E_{eff}, \nu_{eff}, K_{eff}, G_{eff}$ theo thể tích lỗ
Cơ lý biến thiên FGM (mẫu A, B) T3 Tuần hoàn (PBC) Phân bố thành phần ma trận hằng số vật liệu
Đa tinh thể nhôm ($Al$) T3 Tuần hoàn (PBC) Mô đun đàn hồi kháng trượt hữu hiệu dị hướng

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 đa tỉ lệ đàn hồi cho không gian ba chiều bằng phần tử RVE 3D. Mối quan hệ giữa biến dạng vĩ mô và chuyển vị cưỡng bức trên biên RVE 3D được thiết lập ở bậc nhất. Tác giả áp dụng phần tử khối 8 nút (H8) để giải quyết các bài toán vi cấu trúc 3D phức tạp:

  • Vật liệu cấu trúc lớp tuần hoàn: Khảo sát 3 dạng phân bố gồm lớp nằm ngang, lớp đứng và phân bố xen kẽ. Xác định sự biến thiên của mô đun dọc trục $D_{11}$ và mô đun khối $K_{eff}$ khi bề dày mỗi lớp giảm dần.
  • Tấm composite nhiều lớp (laminate plate): Phân tích tấm 4 lớp có hướng sợi chữ thập (cross-ply) và hướng sợi xiên (angle-ply) trên hệ lưới $12 \times 12 \times 12$ phần tử H8.
  • Tấm 3D chịu uốn có cơ lý biến thiên (FGM): Mô hình hóa tấm FGM cấu tạo từ kim loại nhôm oxit ($Al_2O_3$) và gốm Zirconia ($ZrO_2$) dạng 1 lớp và 3 lớp, trích xuất ma trận hằng số vật liệu hữu hiệu 3D.

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 nghiên cứu bài toán đa tỉ lệ cho kết cấu tấm mỏng chịu uốn (lực vuông góc với mặt phẳng tấm) bằng cách xây dựng phần tử tấm đại diện vi mô (RPE - Representative Plate Element). Mối liên hệ động học giữa độ cong vĩ mô ($\kappa_{xx}, \kappa_{yy}, 2\kappa_{xy}$) và chuyển vị cưỡng bức trên biên RPE là quan hệ bậc hai. Chuyển vị tấm được rời rạc hóa bằng phần tử tam giác tương thích $C^1$ HCT. Các ví dụ số gồm:

  • Tấm mỏng vi mô khoét lỗ hình vuông: Tính toán độ võng và sự biến thiên của các mô đun $E_{eff}, \nu_{eff}, K_{eff}, G_{eff}$ theo kích thước lỗ.
  • Tấm nhiều lớp có lỗ tròn: Khảo sát ảnh hưởng của thể tích lỗ rỗng $V_f$ đến độ võng không thứ nguyên tại tâm tấm và biên tựa chu vi, so sánh trực tiếp kết quả giữa phương pháp đa tỉ lệ đồng nhất hóa và mô hình đơn tỉ lệ (single-scale).

Chương 6: Miền cường độ hữu hiệu cho vật liệu theo tiêu chuẩn Hill

Chương 6 trình bày phương pháp xác định miền cường độ hữu hiệu (mặt chảy dẻo hữu hiệu) cho vật liệu không đồng nhất đối xứng kéo - nén bằng cách kết hợp lý thuyết phân tích giới hạn động học và kỹ thuật đồng nhất hóa. Bài toán phân tích giới hạn cận trên của phần tử RVE được chuyển đổi thành bài toán tối ưu hóa hình nón bậc hai (SOCP) với biến dạng dẻo và hàm tiêu tán dẻo theo tiêu chuẩn Hill. Các bài toán số được khảo sát bao gồm:

  • Tấm có lỗ hình tròn và hình chữ nhật chịu kéo/nén: Thiết lập miền cường độ vĩ mô, khảo sát sự thay đổi của cường độ kéo dọc trục theo góc phương tải $\alpha$ và tỉ lệ bán kính $r/a$; mô tả cơ cấu phá hoại dẻo và dựng mặt chảy dẻo 3D ($\Sigma_{11}, \Sigma_{22}, \Sigma_{12}$).
  • Tấm kim loại gia cường cốt sợi: Xác định mặt chảy dẻo và cơ cấu phá hoại dưới tải trọng dọc trục hai phương cho vật liệu mẫu A và mẫu B.
  • Tấm khoét hai lỗ và nhiều lỗ tuần hoàn chịu kéo/nén: Đánh giá sự tương tác vị trí lỗ rỗng đến khả năng chịu lực giới hạn vĩ mô và thời gian tính toán của thuật toán SOCP.

Chương 7: Miền cường độ hữu hiệu cho vật liệu theo tiêu chuẩn Tsai-Wu

Chương 7 mở rộng mô hình phân tích giới hạn đa tỉ lệ cho vật liệu bất đẳng hướng bất đối xứng kéo - nén tuân theo tiêu chuẩn Tsai-Wu ($\Sigma_{Ytx} \neq \Sigma_{Yty} \neq \Sigma_{Ycx} \neq \Sigma_{Ycy}$). Năng lượng tiêu tán dẻo được xây dựng tường minh theo biến dạng dẻo và đưa về dạng quy hoạch nón bậc hai SOCP. Khảo sát số gồm 3 trường hợp:

  1. Composite gia cường cốt sợi tròn: Xác định mặt dẻo hữu hiệu và cơ cấu phá hoại dưới tải trọng đơn trục và đa trục, đánh giá ảnh hưởng của thể tích sợi $V_f$.
  2. Vật liệu có lỗ rỗng tròn tuần hoàn: So sánh mặt dẻo hữu hiệu khi vật liệu nền tuân theo các tiêu chuẩn khác nhau (von Mises, Hill, Tsai-Wu).
  3. Vật liệu có 16 lỗ rỗng phân bố ngẫu nhiên ($V_f = 0.1$): Xác định năng lượng tiêu tán dẻo, phân bố ứng suất giới hạn vĩ mô và xây dựng miền cường độ hữu hiệu.

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

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

  • Bài toán tấm phẳng 2D và 3D đàn hồi: Ưu điểm là tính toán chính xác hằng số đàn hồi hữu hiệu với điều kiện biên tuần hoàn khử bậc tự do; hạn chế là giả thiết biến dạng bé và chi phí tăng khi chia lưới 3D mịn.
  • Bài toán tấm mỏng chịu uốn: Ưu điểm là sử dụng phần tử HCT 9 bậc tự do đảm bảo tính tương thích $C^1$ độ dốc pháp tuyến; hạn chế là chỉ áp dụng cho lý thuyết tấm mỏng Kirchhoff, chưa xét biến dạng trượt ngang của tấm dày.
  • Phân tích giới hạn theo tiêu chuẩn Hill và Tsai-Wu: Ưu điểm là giải trực tiếp tải trọng giới hạn và cơ cấu phá hoại thông qua bài toán tối ưu SOCP giải bằng Mosek có tốc độ hội tụ cao; hạn chế là giả thiết vật liệu cứng dẻo lý tưởng bỏ qua giai đoạn tái bền/mềm hóa.

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

Chương 9 tổng kết toàn bộ các kết quả đạt được về mặt phương pháp tính đa tỉ lệ đàn hồi và thiết kế dẻo ngoài miền đàn hồi, đồng thời đề xuất các định hướng nghiên cứu mở rộng tiếp theo.

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

Những đóng góp mới về mặt khoa học và lý luận

  1. Phát triển phương pháp đa tỉ lệ tính toán hằng số đàn hồi hữu hiệu cho kết cấu tấm phẳng 2D, kết cấu 3D và tấm mỏng chịu uốn: Luận án đã xây dựng thuật toán áp đặt điều kiện biên tuần hoàn (PBC) và tuyến tính (LBC) thông qua kỹ thuật nén bậc tự do (condensation technique), loại bỏ hoàn toàn các bậc tự do phụ thuộc trên biên của phần tử đại diện RVE/RPE, nâng cao độ chính xác khi xác định ma trận hằng số đàn hồi hữu hiệu.
  2. Xây dựng mô hình đa tỉ lệ cho tấm mỏng chịu uốn sử dụng phần tử tương thích $C^1$ HCT: Thiết lập mối liên hệ bậc hai giữa trường độ cong vĩ mô và chuyển vị cưỡng bức biên vi mô, ứng dụng phần tử HCT 9 bậc tự do giúp đảm bảo tính tương thích chuyển vị và đạo hàm góc xoay pháp tuyến giữa các phần tử con.
  3. Thiết lập khung phân tích giới hạn đa tỉ lệ dạng tối ưu hóa hình nón bậc hai (SOCP) cho vật liệu bất đẳng hướng: Chuyển đổi thành công bài toán phân tích giới hạn cận trên của phần tử RVE vi mô tuân theo tiêu chuẩn dẻo Hill và Tsai-Wu thành bài toán tối ưu hóa nón bậc hai SOCP. Phương pháp cho phép tìm trực tiếp các điểm ứng suất giới hạn vĩ mô mà không cần lặp tải từng bước phi tuyến đàn dẻo.
  4. Xây dựng mặt chảy dẻo hữu hiệu tổng quát cho vật liệu composite và vật liệu có lỗ rỗng: Kết hợp phương pháp bình phương cực tiểu (LSA) với tập nghiệm ứng suất giới hạn vi mô để xác định các hệ số của ma trận $P$ và $Q$, cung cấp phương trình mặt dẻo hữu hiệu 3D cho vật liệu có cấu trúc lỗ rỗng phức tạp (tuần hoàn và ngẫu nhiên 16 lỗ).

Bảng tổng hợp các thuật toán và mô hình số trong luận án

Bài toán đa tỉ lệ Cấp độ hình học Phần tử sử dụng Điều kiện biên Hàm mục tiêu / Kỹ thuật đồng nhất Tiêu chuẩn dẻo vi mô
Đàn hồi tấm phẳng 2D T3, Q4, ES-T3 PBC, LBC Định lý trung bình thể tích, nén bậc tự do Đàn hồi tuyến tính Hooke
Đàn hồi không gian 3D H8 PBC Đồng nhất hóa thể tích ma trận độ cứng 3D Đàn hồi tuyến tính Hooke
Đàn hồi tấm chịu uốn 2D (Uốn) HCT (9 DOFs) PBC, LBC Đồng nhất hóa độ cong - mô men tấm Kirchhoff Đàn hồi tuyến tính Hooke
Thiết kế dẻo Hill 2D Tam giác PBC Cực tiểu hóa tiêu tán dẻo (SOCP qua Mosek) Tiêu chuẩn Hill, von Mises
Thiết kế dẻo Tsai-Wu 2D Tam giác PBC Cực tiểu hóa tiêu tán dẻo (SOCP qua Mosek) Tiêu chuẩn Tsai-Wu, Mohr-Coulomb

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

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

  • Giả thiết cơ học vi mô: Ở bài toán phân tích giới hạn ngoài miền đàn hồi, vật liệu nền vi mô chỉ mới được khảo sát dưới mô hình cứng dẻo lý tưởng, chưa xét đến các hiện tượng phi tuyến vật liệu phức tạp như biến cứng (tái bền), hóa mềm dẻo hoặc tổn thương tích lũy.
  • Giới hạn lý thuyết tấm: Bài toán tấm chịu uốn mới chỉ dừng lại ở lý thuyết tấm mỏng Kirchhoff (CPT) với biến dạng bé, chưa xét đến ảnh hưởng của biến dạng trượt ngang trong tấm dày (lý thuyết Mindlin - Reissner) và phi tuyến hình học (biến dạng lớn).
  • Ràng buộc hình học RVE: Các phần tử đại diện vi mô ngoài miền đàn hồi chủ yếu được xây dựng trên miền hình vuông đơn vị hai chiều, chưa mở rộng phân tích giới hạn cho phần tử đại diện không gian 3D chịu trạng thái ứng suất khối phức hợp.

Hướng phát triển tiếp theo:

  • Mở rộng bài toán phân tích giới hạn đa tỉ lệ cho mô hình vật liệu đàn dẻo có xét đến hiện tượng tái bền và phá hủy cục bộ tại cấp độ vi mô.
  • Phát triển phần tử tấm đại diện chịu uốn cho kết cấu tấm dày nhiều lớp có xét đến biến dạng cắt bậc cao.
  • Triển khai bài toán thiết kế dẻo cho phần tử RVE không gian ba chiều (3D) với phân bố pha bất kỳ và phân bố lỗ rỗng ngẫu nhiên trong không gian.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu trong lĩnh vực Cơ học tính toán, Cơ kỹ thuật và Khoa học vật liệu:

  • Đối với nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật, Cơ học kết cấu: Luận án cung cấp chi tiết quy trình toán học và thuật toán áp đặt điều kiện biên tuần hoàn, kỹ thuật nén bậc tự do và cách thức lập trình phần tử tấm HCT tương thích $C^1$.
  • Đối với các nhà nghiên cứu cơ học tính toán: Cung cấp phương pháp chuyển đổi bài toán phân tích giới hạn cơ học phi tuyến thành bài toán tối ưu hóa hình nón bậc hai (SOCP), giải hiệu quả bằng công cụ Mosek với thời gian tính toán tối ưu.
  • Đối với kỹ sư thiết kế kết cấu composite và vật liệu tiên tiến: Cung cấp bộ công thức giải tích và quy trình xác định thông số đàn hồi hữu hiệu ($E_{eff}, \nu_{eff}, G_{eff}, K_{eff}$) cùng các phương trình mặt dẻo hữu hiệu phục vụ tính toán trực tiếp kết cấu ở cấp độ vĩ mô mà không cần mô phỏng vi cấu trúc phức tạp.

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

1. Kỹ thuật nén bậc tự do (condensation technique) được áp dụng như thế nào khi xử lý điều kiện biên tuần hoàn của RVE?

Khi áp đặt điều kiện biên tuần hoàn (PBC), chuyển vị của các cặp nút đối xứng trên biên phần tử đại diện vi mô bị ràng buộc lẫn nhau và phụ thuộc vào chuyển vị của các nút góc (vốn được xác định từ biến dạng vĩ mô). Luận án sử dụng kỹ thuật nén bậc tự do để khử toàn bộ các bậc tự do phụ thuộc này khỏi hệ phương trình đại số tổng thể, giúp hệ phương trình chỉ còn chứa các bậc tự do độc lập, đảm bảo tính đối xứng, xác định dương của ma trận độ cứng và ổn định số khi giải.

2. Phần tử tấm HCT được cấu tạo và rút gọn bậc tự do ra sao trong bài toán tấm mỏng chịu uốn?

Phần tử HCT (Hsieh-Clough-Tocher) là phần tử tam giác tương thích $C^1$ gồm 3 phần tử tam giác con bên trong kết nối tại trọng tâm tấm. Ban đầu phần tử có 12 bậc tự do (3 chuyển vị võng và 6 góc xoay tại 3 đỉnh, cùng 3 góc xoay pháp tuyến tại trung điểm 3 cạnh ngoài). Ba bậc tự do tại nút trọng tâm bên trong được khử bằng điều kiện liên tục góc xoay giữa các tam giác con. Sau đó, bằng giả thuyết đạo hàm góc xoay pháp tuyến biến thiên tuyến tính dọc theo từng cạnh biên, 3 bậc tự do tại trung điểm cạnh được biểu diễn qua góc xoay tại các đỉnh, rút gọn phần tử HCT về 9 bậc tự do độc lập.

3. Làm thế nào bài toán phân tích giới hạn cận trên cho RVE vi mô được chuyển về bài toán tối ưu hóa hình nón bậc hai (SOCP)?

Trong phương pháp cận trên, hệ số tải trọng giới hạn được định nghĩa là giá trị cực tiểu của năng lượng tiêu tán dẻo dưới các ràng buộc tương thích động học, điều kiện biên tuần hoàn và điều kiện chuẩn hóa công ngoại lực ($W_{\text{ext}} = 1$). Bằng cách sử dụng luật chảy dẻo kết hợp, hàm tiêu tán dẻo (vốn có dạng căn bậc hai của dạng toàn phương theo biến dạng đối với tiêu chuẩn Hill và Tsai-Wu) được biểu diễn tuyến tính thông qua việc đặt các biến phụ hình nón. Bài toán tối ưu hóa phi tuyến được chuyển hoàn toàn về dạng quy hoạch nón bậc hai chuẩn (SOCP), giải trực tiếp và hội tụ nhanh bằng phần mềm Mosek.

4. Sự khác biệt cơ bản giữa mặt chảy dẻo hữu hiệu theo tiêu chuẩn Hill và tiêu chuẩn Tsai-Wu trong luận án là gì?

Tiêu chuẩn Hill là dạng mở rộng của tiêu chuẩn von Mises cho vật liệu bất đẳng hướng nhưng đối xứng kéo - nén, tức là cường độ chịu kéo và chịu nén theo mỗi phương tọa độ là bằng nhau ($\Sigma_{Ytx} = \Sigma_{Yty}$ và $\Sigma_{Ycx} = \Sigma_{Ycy}$). Trong khi đó, tiêu chuẩn Tsai-Wu là tiêu chuẩn tổng quát bậc hai áp dụng cho vật liệu bất đẳng hướng bất đối xứng kéo - nén, trong đó khả năng chịu kéo và chịu nén theo từng phương là hoàn toàn độc lập và khác nhau ($\Sigma_{Ytx} \neq \Sigma_{Yty} \neq \Sigma_{Ycx} \neq \Sigma_{Ycy}$), cho phép mô tả chính xác ứng xử phá hoại của các vật liệu có ứng xử giòn, composite nền polyme hoặc vật liệu có phân bố pha bất đối xứng.

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 khung phương pháp đa tỉ lệ hoàn chỉnh cho kết cấu tấm không đồng nhất trong cả miền ứng xử đàn hồi và phân tích giới hạn dẻo. Bằng việc kết hợp lý thuyết đồng nhất hóa thể tích, điều kiện biên tuần hoàn nén bậc tự do, phần tử tấm tương thích $C^1$ HCT và giải thuật tối ưu hóa nón bậc hai SOCP thông qua công cụ Mosek, công trình đã xác định chính xác các hằng số đàn hồi hữu hiệu và thiết lập các hàm miền cường độ hữu hiệu dạng Hill và Tsai-Wu cho nhiều dạng vi cấu trúc phức tạp. Kết quả nghiên cứu cung cấp cơ sở khoa học và công cụ tính toán số tin cậy cho việc mô phỏng và thiết kế các kết cấu tấm vật liệu mới trong kỹ thuật công trình và cơ học kết cấu.