Tổng quan luận án

Tấm composite sandwich lõi lượn sóng, đặc biệt là carton lượn sóng, là dạng kết cấu được ứng dụng rộng rãi trong kỹ thuật bao bì, công nghiệp đóng tàu, vận tải và xây dựng nhờ sở hữu tỷ lệ độ bền trên khối lượng cao, khả năng hấp thu xung lực tốt và khả năng tái chế vượt trội. Trong quy trình thiết kế công nghiệp, phương pháp phần tử hữu hạn (FEM) là công cụ tiêu chuẩn để dự đoán ứng xử kết cấu dưới các điều kiện tải trọng khác nhau. Tuy nhiên, việc mô hình hóa trực tiếp cấu trúc hình học 3D chi tiết của tấm carton lõi lượn sóng gặp phải trở ngại lớn do đòi hỏi khối lượng phần tử rất lớn, thời gian chuẩn bị mô hình kéo dài và chi phí tính toán CPU cao.

Khoảng trống nghiên cứu được tác giả xác định xuất phát từ thực trạng phần lớn các nghiên cứu về phương pháp đồng nhất hóa tấm carton lượn sóng hiện nay mới dừng lại ở phạm vi ứng xử đàn hồi tuyến tính. Các mô hình này chưa phản ánh đầy đủ các đặc tính cơ học phi tuyến của vật liệu giấy trong thực tế, bao gồm biến dạng dẻo, ảnh hưởng của lực cắt ngang và phản ứng của kết cấu dưới các trạng thái tải trọng phức tạp hoặc tải trọng động. Một số nghiên cứu bước đầu đề xuất mô hình đàn dẻo nhưng vẫn đơn giản hóa ứng xử vật liệu hoặc chưa được kiểm chứng đầy đủ trong các điều kiện tải trọng kết hợp và tải trọng va đập. Thêm vào đó, việc xác định chính xác các tham số vật liệu đầu vào phi tuyến cho mô hình số cũng gặp khó khăn do tính chất dị hướng phức tạp của từng lớp giấy cấu thành.

Mục tiêu nghiên cứu của luận án được xác định cụ thể như sau:

  1. Phát triển mô hình phần tử hữu hạn 2D tương đương đàn hồi – dẻo (mô hình H-2D) thay thế cho mô hình 3D của tấm carton lõi lượn sóng nhằm giảm thiểu thời gian tính toán và thời gian xây dựng mô hình khi phân tích ứng xử cơ học dưới các điều kiện tải trọng khác nhau.
  2. Xây dựng quy trình nhận dạng ngược sử dụng thuật toán tối ưu hóa để xác định chính xác các tham số trong mô hình ứng xử cơ học đàn dẻo của vật liệu giấy cấu thành.
  3. Ứng dụng mô hình 2D tương đương để phân tích ứng xử động học, xác định độ bền và xây dựng đường cong phá hủy (Damage Boundary Curve - DBC) của hộp carton khi chịu rung động và va đập trong quá trình vận chuyển.

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

  • Đối tượng nghiên cứu: Tấm carton lõi lượn sóng và kết cấu hộp bao bì chế tạo từ carton lượn sóng.
  • Phạm vi nghiên cứu về nội dung: Tập trung vào mô hình giải tích, mô hình hóa số bằng phương pháp phần tử hữu hạn trên phần mềm Abaqus, xây dựng chương trình con UGENS và VUMAT, kết hợp thuật toán tối ưu hóa trong phần mềm ModeFrontier; thực nghiệm cơ học kéo mẫu giấy tách lớp trong phòng thí nghiệm để phục vụ bài toán nhận dạng ngược.
  • Phạm vi không gian và thời gian: Nghiên cứu được thực hiện tại Trường Đại học Kỹ thuật Công nghiệp – Đại học Thái Nguyên, hoàn thành và bảo vệ năm 2025.

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

Luận án đã hệ thống hóa các hướng nghiên cứu trong nước và quốc tế liên quan đến cấu trúc sandwich, cơ học vật liệu giấy và các kỹ thuật đồng nhất hóa:

Các hướng nghiên cứu về đồng nhất hóa tấm sandwich và carton lượn sóng:

  • Phương pháp giải tích và lý thuyết tấm nhiều lớp: Aboura và cộng sự [4] xây dựng mô hình ứng xử đàn hồi cho carton sóng dựa trên lý thuyết tấm nhiều lớp cổ điển, phân tích sự ảnh hưởng của chiều dày tấm tới mô đun EMD và các độ cứng uốn D11, D22, cũng như ảnh hưởng của bước sóng tới ECD. Carlsson [30] và Nordstrand [31, 32] sử dụng phương pháp phân tích kết hợp nguyên lý năng lượng để xác định các đặc tính cơ học hiệu dụng, làm rõ vai trò của độ cứng cắt ngang trên mặt cắt MD và tải trọng tới hạn khi tấm chịu nén [33, 34]. Luo và cộng sự [29] phát triển phương pháp tính toán độ cứng uốn; Abbès và Guo [36], Dương và cộng sự [37] cải tiến công thức tính toán độ cứng xoắn.
  • Phương pháp đồng nhất hóa dựa trên phần tử hữu hạn và năng lượng tương đương: Biancolini [22, 42] đề xuất kỹ thuật đồng nhất hóa số sử dụng mô hình phần tử thể tích đại diện (RVE/VER) kết hợp điều kiện biên tuần hoàn để xác định tensor độ cứng màng, uốn và mở rộng cho độ cứng cắt ngang. Talbi và cộng sự [19] áp dụng lý thuyết độ cứng tương đương để quy đổi kết cấu nhiều lớp thành tấm đồng nhất tương đương, kiểm chứng qua các trạng thái kéo, nén, uốn, xoắn đàn hồi. Young-Jo Cheon [41] kết hợp phương pháp năng lượng biến dạng tương đương (DEEM) và lý thuyết tấm cổ điển để lập ma trận độ cứng tấm lượn sóng chịu kéo và uốn. Damian Mrówczyński và Tomasz Garbowski [39, 43] phát triển mô hình DEEM cho tấm 5 lớp và nghiên cứu ảnh hưởng của các khuyết tật hình học, đường gấp nếp và lỗ đục tới sự suy giảm độ cứng đàn hồi [11]. Huimin Li [44] áp dụng mô hình đồng nhất hóa cho tấm sandwich lõi lượn sóng hình thang 2D.
  • Hướng nghiên cứu kết hợp trong nước: Các tác giả trong nước đã sử dụng mô hình đồng nhất hóa kết hợp với mô hình Halpin–Tsai mở rộng để tính độ cứng tấm composite dạng sóng [22], hoặc kết hợp với lý thuyết vỏ Donnell và phi tuyến hình học von Karman cho cấu trúc xoắn ốc dạng sóng [23], và tích hợp trường nhiệt vào phân tích ứng xử của vỏ hình xuyến sandwich FG-GRC lõi sóng [24].

Các hướng nghiên cứu về ứng xử vật liệu giấy và mô hình đàn dẻo:

  • N. Stenberg [48] khảo sát ứng xử đàn dẻo theo chiều dày của giấy; P. Mäkelä [49] đề xuất mô hình đàn dẻo trực giao. Nhóm nghiên cứu [50] phát triển mô hình đàn dẻo phi đẳng hướng của giấy và bìa giấy, kết hợp đàn hồi tuyến tính và mặt chảy dẻo phụ thuộc vào các hướng MD, CD và ZD. V. Luong và cộng sự [47, 51] sử dụng phương pháp phần tử hữu hạn đánh giá khả năng chịu va chạm của hộp carton và xây dựng mô hình đồng nhất hóa mô tả ứng xử sau giai đoạn đàn hồi cho tấm carton sóng.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết: Phần lớn các nghiên cứu đã công bố giới hạn ở ứng xử đàn hồi tuyến tính hoặc chỉ xét biến dạng dẻo trong các điều kiện tải trọng đơn giản. Các mô hình hiện có thiếu sự kiểm chứng toàn diện dưới các chế độ tải phức tạp như tải trọng va chạm, tải trọng chu kỳ lặp (chất tải - dỡ tải) và rung động động học. Hơn nữa, việc xác định đồng thời bộ thông số đàn hồi và dẻo dị hướng của các lớp giấy cấu thành chưa có quy trình nhận dạng ngược tự động hóa. Luận án giải quyết khoảng trống này bằng cách xây dựng một mô hình 2D đồng nhất hóa hoàn chỉnh có khả năng mô tả liên tục từ giai đoạn đàn hồi tuyến tính đến giai đoạn biến dạng dẻo phi tuyến, tích hợp thuật toán nhận dạng ngược tham số vật liệu và kiểm chứng qua các bài toán kết cấu thực tế.

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

Cơ sở lý thuyết và khung phân tích

Luận án vận dụng hệ thống lý thuyết cơ học kết cấu tấm vỏ và cơ học vật liệu dị hướng:

  • Lý thuyết tấm mỏng Kirchhoff–Love: Áp dụng giả thiết đoạn pháp tuyến giữ thẳng và vuông góc với mặt trung bình sau biến dạng, bỏ qua biến dạng cắt ngang ($\gamma_{xz} = \gamma_{yz} = 0$). Trường chuyển vị được xác định theo công thức $u(x,y,z) = -z \frac{\partial w_0}{\partial x}$, $v(x,y,z) = -z \frac{\partial w_0}{\partial y}$, $w(x,y,z) = w_0(x,y)$, với độ cứng uốn $D = \frac{E h^3}{12(1-\nu^2)}$.
  • Lý thuyết tấm Mindlin (First-order Shear Deformation Theory - FSDT): Kể đến ảnh hưởng của biến dạng cắt ngang thông qua các góc xoay độc lập của pháp tuyến $\beta_x$ và $\beta_y$: $u = u_0 + z\beta_x$, $v = v_0 + z\beta_y$, $w = w_0$. Trường biến dạng phân rã thành thành phần màng và thành phần uốn: $\boldsymbol{\varepsilon} = \boldsymbol{\varepsilon}_m + z\boldsymbol{\kappa}$.
  • Lý thuyết tấm nhiều lớp (CLPT mở rộng): Thiết lập quan hệ giữa nội lực (lực màng $\mathbf{N}$, mô men uốn - xoắn $\mathbf{M}$, lực cắt ngang $\mathbf{T}$) và biến dạng tổng quát thông qua hệ ma trận độ cứng: $$\begin{bmatrix} \mathbf{N} \ \mathbf{M} \ \mathbf{T} \end{bmatrix} = \begin{bmatrix} \mathbf{A} & \mathbf{B} & \mathbf{0} \ \mathbf{B} & \mathbf{D} & \mathbf{0} \ \mathbf{0} & \mathbf{0} & \mathbf{F} \end{bmatrix} \begin{bmatrix} \boldsymbol{\varepsilon}_m \ \boldsymbol{\kappa} \ \boldsymbol{\gamma} \end{bmatrix}$$ Trong đó $\mathbf{A}$ là ma trận độ cứng màng, $\mathbf{B}$ là ma trận tương tác màng - uốn, $\mathbf{D}$ là ma trận độ cứng uốn - xoắn, và $\mathbf{F}$ là ma trận độ cứng cắt ngang.
  • Mô hình ứng xử vật liệu đàn dẻo: Giấy và bìa carton là vật liệu trực hướng có sự khác biệt rõ rệt theo ba phương chính: phương máy chạy (Machine Direction - MD), phương ngang máy (Cross Direction - CD) và phương chiều dày (Thickness Direction - ZD). Ứng xử dẻo được mô hình hóa bằng mô hình dẻo tương đương đẳng hướng (Isotropic Plastic Equivalent - IPE) kết hợp với thuật toán hồi quy xuyên tâm (Radial Return) để cập nhật trạng thái ứng suất tại các điểm tích phân.

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

  • Phương pháp giải tích và số: Thiết lập ô cơ sở đại diện tuần hoàn (VER) với chiều dài bằng bước sóng $P$ của rãnh sóng, hình học đường sóng trung bình mô tả theo hàm sin $h(x) = \frac{h_c}{2} \sin\left(\frac{2\pi x}{P}\right)$. Chia VER thành các lát cắt vi phân $dx$, thực hiện tích phân theo chiều dày từng lát cắt rồi đồng nhất hóa trung bình theo phương $x$ để suy ra các ma trận độ cứng tương đương.
  • Lập trình chương trình con (User Subroutines): Luận án xây dựng hai chương trình con bằng ngôn ngữ Fortran và Python tích hợp trong phần mềm Abaqus:
    • UGENS: Sử dụng để định nghĩa độ ứng xử cơ học phi tuyến của mặt cắt vỏ tổng quát, gán ma trận độ cứng và quản lý tích phân qua chiều dày cho mô hình phần tử hữu hạn 2D tương đương.
    • VUMAT: Sử dụng trong giải thuật phân tích động tường minh (Abaqus/Explicit) nhằm mô tả quan hệ ứng suất – biến dạng đàn hồi - dẻo phi tuyến của vật liệu giấy dị hướng.
  • Phương pháp thực nghiệm và nhận dạng ngược: Tiến hành thí nghiệm kéo tĩnh các mẫu giấy thành phần (tách từ tấm carton sóng loại B) theo các hướng MD, CD và $45^\circ$ trên máy kéo nén vạn năng. Xây dựng mô hình mô phỏng số tương ứng của bài toán thử kéo trong Abaqus, kết nối với phần mềm ModeFrontier để thực hiện tối ưu hóa đa mục tiêu, tự động hiệu chỉnh và xác định bộ thông số đàn hồi và dẻo của mô hình IPE sao cho sai lệch giữa đường cong lực - biến dạng mô phỏng và thực nghiệm là nhỏ nhất.

Nguồn tư liệu và dữ liệu

  • Dữ liệu hình học các dạng sóng chuẩn (A, B, C, E, F) theo tiêu chuẩn kỹ thuật bao bì quốc tế.
  • Số liệu thực nghiệm kéo mẫu giấy lớp mặt (liner) và lớp lõi (flute) thực hiện trực tiếp trong phòng thí nghiệm.
  • Dữ liệu mô phỏng phần tử hữu hạn từ các mô hình 3D chi tiết, mô hình tấm 2D đồng nhất hóa đàn hồi - dẻo (H-2D) và mô hình đàn hồi tương đương (E-2D).
  • Dữ liệu thực nghiệm uốn 3 điểm và thử nghiệm rung động hộp carton từ các công bố thực nghiệm đối chứng.

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

CHƯƠNG 1. NGHIÊN CỨU TỔNG QUAN

Chương 1 trình bày tổng quan cấu tạo, phân loại và đặc tính cơ học của kết cấu tấm sandwich và tấm carton lõi lượn sóng. Tấm carton sóng gồm các lớp giấy phẳng bên ngoài (liner) làm từ bột giấy kraft gỗ mềm định lượng 125 đến 440 $\text{g/m}^2$, và lớp lõi lượn sóng (flute) định lượng 80 đến 180 $\text{g/m}^2$. Quá trình sản xuất công nghiệp định hướng 70–80% sợi xenlulo theo phương MD, 20% theo phương CD và 10% theo phương ZD, tạo ra tính trực hướng rõ rệt cho vật liệu.

Bảng tổng hợp đặc tính hình học các loại rãnh sóng phổ biến được trình bày trong luận án:

Loại rãnh sóng Chiều cao rãnh $h_c$ (mm) Bước sóng $P$ (mm) Số rãnh trên 1 mét dài (rãnh/m)
A 4,8 8,0 – 9,5 $108 \pm 10$
B 3,2 5,5 – 6,5 $154 \pm 10$
C 4,0 6,8 – 7,9 $128 \pm 10$
E 1,6 3,0 – 3,5 $295 \pm 13$
F 0,8 1,9 – 2,6 $420 \pm 10$

Chương này cũng phân tích dữ liệu quy mô thị trường bao bì tại Việt Nam giai đoạn 2015–2020 và cấu trúc vật liệu bao bì tại Cộng hòa Pháp năm 2016 (Hội đồng Bao bì Quốc gia CNE), chứng minh tỷ lệ tái chế của bao bì giấy carton vượt mốc 90%, cao nhất trong các nhóm vật liệu đóng gói. Phân tích tổng quan chỉ rõ các phương pháp đồng nhất hóa phổ biến như CLPT, FSDT, DEEM và AH, từ đó xác lập tính cấp thiết của việc phát triển mô hình 2D đàn dẻo thay thế mô hình 3D.

CHƯƠNG 2. XÂY DỰNG MÔ HÌNH PHẦN TỬ HỮU HẠN TƯƠNG ĐƯƠNG

Chương 2 thiết lập cơ sở toán học và giải thuật tính toán cho mô hình 2D tương đương. Phân loại hình học tấm dựa trên tỷ số chiều dày $h$ trên chiều dài nhỏ nhất $L_{\min}$: màng mỏng ($h/L_{\min} \le 1/80$), tấm mỏng ($1/80 < h/L_{\min} \le 1/5$), và tấm dày ($h/L_{\min} > 1/5$).

Nội dung xây dựng mô hình đồng nhất hóa bao gồm:

  1. Đồng nhất hóa đàn hồi: Xác định ô đơn vị tuần hoàn (VER) có chiều dài bằng bước sóng $P$. Tọa độ bề mặt sóng được mô tả qua hàm $h(x) = \frac{h_c}{2} \sin\left(\frac{2\pi x}{P}\right)$. Cắt VER thành các phân tố $dx$, biến đổi ma trận độ cứng cục bộ của từng lớp nghiêng về hệ tọa độ toàn cục $(x, y, z)$, sau đó tích phân theo chiều dày để tính các ma trận $\mathbf{A}(x)$, $\mathbf{B}(x)$, $\mathbf{D}(x)$, $\mathbf{F}(x)$. Độ cứng trung bình của toàn tấm được xác lập bằng phép tích phân trên chu kỳ sóng: $$\bar{\mathbf{A}} = \frac{1}{P}\int_0^P \mathbf{A}(x)dx; \quad \bar{\mathbf{B}} = \frac{1}{P}\int_0^P \mathbf{B}(x)dx; \quad \bar{\mathbf{D}} = \frac{1}{P}\int_0^P \mathbf{D}(x)dx$$
  2. Đồng nhất hóa đàn hồi – dẻo: Tích hợp quan hệ ứng suất – biến dạng phi tuyến của vật liệu giấy dựa trên tiêu chuẩn chảy dẻo IPE. Xây dựng thuật toán hiệu chỉnh ứng suất Radial Return để cập nhật tensor ứng suất và biến dạng dẻo tích lũy tại mỗi bước gia tải.
  3. Cấu trúc chương trình con: Thiết lập mã nguồn Fortran cho VUMAT và UGENS. Trong đó, UGENS định nghĩa cách gán ma trận độ cứng tương đương tại từng điểm tích phân qua chiều dày phần tử vỏ, trong khi VUMAT tính toán ứng xử phi tuyến vật liệu trong phân tích động học Abaqus/Explicit.

CHƯƠNG 3. KIỂM CHỨNG MÔ HÌNH ĐỒNG NHẤT HÓA

Chương 3 triển khai kiểm chứng độ tin cậy của mô hình đề xuất qua hai giai đoạn: nhận dạng ngược tham số vật liệu và so sánh mô phỏng giữa mô hình 3D chi tiết với mô hình 2D tương đương (H-2D).

Quy trình nhận dạng ngược:

  • Mẫu thí nghiệm kéo được chế tạo theo tiêu chuẩn từ các lớp giấy thành phần tách ra từ tấm carton sóng B (gồm giấy mặt 1, giấy mặt 3 và giấy sóng 2). Thử nghiệm kéo đơn trục theo các hướng MD, CD và $45^\circ$ cung cấp các đường cong quan hệ lực – độ giãn dài thực nghiệm.
  • Xây dựng bài toán tối ưu hóa nhận dạng ngược trong ModeFrontier kết hợp Abaqus. Thuật toán tối ưu hóa tự động duyệt qua các biến thiết kế thuộc tính đàn hồi và dẻo của mô hình IPE. Kết quả đạt sự hội tụ tại vòng lặp tính toán số 168, khớp nối chính xác giữa đường cong số và đường cong thực nghiệm cho cả ba lớp giấy.

Kiểm chứng mô hình đàn hồi và đàn dẻo:

  • So sánh các thông số độ cứng giữa mô hình 3D và mô hình H-2D trong giai đoạn đàn hồi: Độ cứng kéo theo phương MD, độ cứng kéo theo phương CD, độ cứng uốn quanh trục $x$, độ cứng uốn quanh trục $y$, độ cứng cắt trong mặt phẳng $xy$ và độ cứng xoắn theo hai mặt cắt MD, CD. Độ sai lệch giữa hai mô hình duy trì ở mức rất nhỏ.
  • Kiểm chứng đàn hồi – dẻo qua các bài toán chịu tải phức tạp:
    1. Tấm phẳng chịu kéo theo phương MD và CD;
    2. Tấm có lỗ khoét tập trung ứng suất chịu kéo: Mô hình H-2D mô tả chính xác mức độ tập trung ứng suất và trường biến dạng dẻo quanh mép lỗ tương đương mô hình 3D;
    3. Tấm chịu nén dọc trục theo phương CD;
    4. Phân tích dao động tự do: Xác định tần số tự nhiên và dạng dao động (Mode 1, Mode 2, Mode 3) cho thấy sự tương thích giữa hai mô hình;
    5. Đánh giá thời gian tính toán: Mô hình H-2D tiết kiệm đáng kể thời gian CPU và tài nguyên bộ nhớ so với mô hình 3D chi tiết.

CHƯƠNG 4. ỨNG DỤNG MÔ HÌNH ĐỒNG NHẤT HÓA

Chương 4 áp dụng mô hình 2D đồng nhất hóa vào các bài toán kỹ thuật thực tế mô phỏng điều kiện làm việc của tấm và hộp bao bì carton:

  1. Dao động của tấm khi chịu uốn và va chạm: Phân tích phản hồi gia tốc của tấm và dầm chịu va chạm xung lực theo hướng MD và hướng CD. So sánh gia tốc kích thích và gia tốc phản hồi giữa mô hình 3D và H-2D cho kết quả trùng khớp cao về cả biên độ và pha dao động.
  2. Mô phỏng chu kỳ chất tải – dỡ tải: Thực hiện gia tải tuần hoàn với hai trường hợp: biên độ tăng dần và biên độ ổn định. Kết quả khẳng định mô hình H-2D mô tả chính xác hiện tượng suy giảm độ cứng và sự tích lũy biến dạng dẻo dư thừa sau mỗi chu kỳ dỡ tải.
  3. Thí nghiệm uốn ba điểm: Mô phỏng thí nghiệm uốn 3 điểm tấm carton lượn sóng, đối chiếu trực tiếp đường cong lực – độ võng và trạng thái ứng suất uốn với kết quả thực nghiệm.
  4. Ứng dụng phân tích rung động hộp carton hoàn chỉnh:
    • Mô hình hóa hộp carton chứa khối hàng tải trọng $M = 16,9\text{ kg}$ chịu rung động ngẫu nhiên và va đập mô phỏng quá trình vận chuyển.
    • So sánh ba mô hình: mô hình 3D chi tiết, mô hình đồng nhất hóa đàn dẻo (H-2D) và mô hình đàn hồi tương đương (E-2D).
    • Xây dựng đường cong biên phá hủy (Damage Boundary Curve - DBC) biểu diễn mối quan hệ giữa gia tốc tới hạn và độ thay đổi vận tốc tới hạn. Kết quả từ mô hình H-2D bám sát đường cong DBC thực nghiệm, chứng minh tính vượt trội so với mô hình đàn hồi E-2D thuần túy.

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

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

  1. Luận án đã phát triển thành công mô hình phần tử hữu hạn 2D tương đương đàn hồi – dẻo (mô hình H-2D) cho kết cấu tấm composite sandwich lõi lượn sóng. Mô hình khắc phục được hạn chế của các phương pháp đồng nhất hóa truyền thống vốn chỉ áp dụng cho vùng đàn hồi tuyến tính.
  2. Xây dựng và tích hợp thành công hai chương trình con UGENS và VUMAT bằng ngôn ngữ Fortran trong môi trường Abaqus, cho phép mô phỏng chính xác ứng xử cơ học phi tuyến, dị hướng của vật liệu giấy và hiện tượng biến dạng dẻo tích lũy dưới các chế độ tải trọng tĩnh và động.
  3. Thiết lập hoàn chỉnh quy trình nhận dạng ngược tự động kết hợp phần mềm tối ưu hóa ModeFrontier và Abaqus để xác định đồng thời các tham số đàn hồi và dẻo của từng lớp giấy cấu thành từ dữ liệu thực nghiệm kéo đơn trục.

Đóng góp mới về mặt thực tiễn và kỹ thuật

  1. Cung cấp giải pháp mô phỏng số rút ngắn thời gian chuẩn bị mô hình FEM và tiết kiệm thời gian tính toán CPU so với mô hình 3D chi tiết, trong khi vẫn duy trì độ chính xác cao về trường ứng suất, biến dạng và phản hồi động học.
  2. Ứng dụng thành công mô hình 2D tương đương trong phân tích độ bền của hộp bao bì carton thực tế tải trọng $16,9\text{ kg}$, xây dựng chính xác đường cong biên phá hủy (DBC) chịu va chạm và rung động trong vận tải, hỗ trợ trực tiếp cho công tác tối ưu hóa thiết kế bao bì bảo vệ sản phẩm.

Kiến nghị và đề xuất kỹ thuật

  • Áp dụng quy trình mô hình hóa tương đương 2D và chương trình con UGENS/VUMAT vào các phần mềm CAE thương mại trong quy trình thiết kế công nghiệp để giảm thiểu chi phí thử nghiệm thực tế phá hủy.
  • Khuyến nghị sử dụng mô hình nhận dạng ngược trong ModeFrontier để chuẩn hóa cơ sở dữ liệu vật liệu giấy công nghiệp tại các nhà máy sản xuất bao bì carton.

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

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

  • Các dữ liệu thực nghiệm kiểm chứng và nhận dạng ngược trong luận án mới chỉ tập trung triển khai trên kết cấu carton lõi lượn sóng đơn loại B (3 lớp), chưa mở rộng đầy đủ sang các dạng sóng phức hợp nhiều lớp (sóng 5 lớp, 7 lớp).
  • Mô hình hiện tại chưa xét đến ảnh hưởng của các yếu tố môi trường bên ngoài như nhiệt độ và độ ẩm tương đối của không khí – những yếu tố có tác động đáng kể đến liên kết hydro trong sợi xenlulo của giấy.
  • Cơ chế phá hủy do tách lớp (delamination) hoặc bong tróc liên kết keo dán giữa lớp giấy sóng và lớp giấy đáy chưa được tích hợp hoàn chỉnh trong mô hình đàn dẻo tổng thể.

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

  • Mở rộng thuật toán đồng nhất hóa đàn dẻo cho các kết cấu tấm sandwich đa lớp sóng (sóng đôi BC, AB hoặc sóng ba) và các cấu trúc lõi phức tạp khác.
  • Nghiên cứu tích hợp trường nhiệt – ẩm vào luật ứng xử vật liệu để mô phỏng sự suy giảm độ bền bao bì trong điều kiện môi trường nhiệt đới ẩm.
  • Bổ sung mô hình phá hủy bề mặt tiếp xúc (cohesive zone model) nhằm mô tả hiện tượng tách lớp keo dán khi tấm chịu tải trọng bóc tách hoặc uốn trượt lớn.

Giá trị tham khảo

Luận án mang lại giá trị tham khảo thiết thực cho nhiều nhóm đối tượng:

  • Nghiên cứu sinh, học viên cao học và nhà nghiên cứu cơ học kết cấu: Cung cấp tài liệu tham khảo về phương pháp đồng nhất hóa giải tích và số, kỹ thuật lập trình chương trình con UGENS, VUMAT trong phần mềm Abaqus, và phương pháp luận thiết lập bài toán nhận dạng ngược thông số vật liệu.
  • Kỹ sư thiết kế kết cấu, kỹ sư bao bì công nghiệp: Ứng dụng trực tiếp mô hình H-2D để tính toán nhanh khả năng chịu lực, độ bền nén hộp (BCT), khả năng chịu va đập và xây dựng đường cong DBC của bao bì vận chuyển mà không cần đầu tư hệ thống máy tính hiệu năng cao để chạy mô hình 3D.
  • Giảng viên và các viện nghiên cứu cơ khí: Sử dụng làm tài liệu giảng dạy chuyên đề cho các học phần Cơ học vật liệu composite, Phương pháp phần tử hữu hạn nâng cao và Tối ưu hóa kết cấu cơ khí.

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

1. Điểm khác biệt mấu chốt giữa mô hình H-2D của luận án và các mô hình đồng nhất hóa trước đây là gì?

Các mô hình đồng nhất hóa truyền thống (như CLPT, DEEM cổ điển) chủ yếu dừng lại ở giai đoạn đàn hồi tuyến tính, coi vật liệu có độ cứng không đổi. Mô hình H-2D của luận án là mô hình 2D tương đương đàn hồi – dẻo, tích hợp thuật toán Radial Return và tiêu chuẩn chảy dẻo IPE thông qua chương trình con UGENS và VUMAT, cho phép mô tả chính xác sự suy giảm độ cứng, biến dạng dẻo tích lũy và ứng xử phi tuyến sau giới hạn đàn hồi dưới các tải trọng tĩnh, chu kỳ và va chạm động học.

2. Thuật toán nhận dạng ngược trong ModeFrontier được vận hành như thế nào?

Quy trình nhận dạng ngược kết hợp phần mềm ModeFrontier và Abaqus để xác định bộ thông số vật liệu tối ưu. Đầu vào là dữ liệu thử kéo thực nghiệm của các mẫu giấy thành phần theo các hướng MD, CD và $45^\circ$. ModeFrontier tự động điều khiển Abaqus chạy mô phỏng kéo với các biến thiết kế thuộc tính vật liệu trong một dải giá trị xác định, sau đó so sánh đường cong lực – biến dạng mô phỏng với thực nghiệm qua hàm mục tiêu giảm thiểu sai số. Quy trình hội tụ tối ưu tại vòng lặp thứ 168.

3. Tại sao giấy và carton lượn sóng lại có tính chất cơ học trực hướng dị biệt giữa hai phương MD và CD?

Tính chất trực hướng xuất phát từ quá trình chế tạo giấy trên máy làm giấy: khoảng 70–80% sợi xenlulo được định hướng sắp xếp theo phương máy chạy (MD), trong khi chỉ có khoảng 20% định hướng theo phương ngang máy (CD) và 10% theo phương chiều dày (ZD). Do đó, độ bền kéo, độ cứng và khả năng chịu lực theo phương MD luôn cao hơn đáng kể so với phương CD.

4. Vai trò cụ thể của hai chương trình con UGENS và VUMAT trong luận án là gì?

UGENS là chương trình con cho phép người dùng định nghĩa ứng xử cơ học phi tuyến của mặt cắt vỏ tổng quát trong Abaqus, quản lý các điểm tích phân số qua độ dày và ma trận độ cứng tương đương của tấm. VUMAT là chương trình con vật liệu dùng trong giải thuật phân tích động tường minh (Abaqus/Explicit), thực hiện cập nhật tensor ứng suất và biến dạng dẻo theo từng bước thời gian cho vật liệu dị hướng khi phân tích các bài toán va chạm, phá hủy và tải trọng động.

5. Kết quả ứng dụng mô hình đồng nhất hóa cho bài toán hộp carton chịu rung động đạt độ chính xác ra sao?

Khi áp dụng cho hộp carton chứa tải trọng $16,9\text{ kg}$ chịu tải trọng rung động và va đập mô phỏng quá trình vận chuyển, mô hình đàn dẻo H-2D cho kết quả xây dựng đường cong phá hủy (Damage Boundary Curve - DBC) bám sát chặt chẽ đường cong đo đạc từ thực nghiệm, vượt trội hoàn toàn so với mô hình thuần đàn hồi E-2D, đồng thời tiết kiệm đáng kể thời gian tính toán so với mô hình 3D chi tiết.

Kết luận

Luận án của nghiên cứu sinh Đào Liên Tiến đã giải quyết trọn vẹn bài toán xây dựng mô hình phần tử hữu hạn 2D tương đương đàn hồi – dẻo (H-2D) thay thế cho mô hình 3D của tấm carton lõi lượn sóng. Thông qua việc kết hợp cơ sở lý thuyết tấm vỏ nhiều lớp, các chương trình con UGENS/VUMAT trong Abaqus và quy trình nhận dạng ngược tham số vật liệu trong ModeFrontier, nghiên cứu đã tạo ra công cụ tính toán cơ học tin cậy cho kết cấu sandwich lõi sóng. Các kết quả kiểm chứng qua bài toán tĩnh, bài toán va chạm chu kỳ và bài toán phân tích rung động hộp bao bì $16,9\text{ kg}$ đã khẳng định tính chính xác và hiệu quả giảm tải thời gian tính toán của mô hình. Đây là công trình nghiên cứu hoàn chỉnh, đóng góp cả về phương pháp luận khoa học cơ học vật liệu lẫn ứng dụng thực tiễn trong ngành kỹ thuật bao bì công nghiệp.