Tổng quan nghiên cứu

Ngành công nghiệp gia công polyme nhiệt dẻo hiện nay chiếm hơn 60% tổng sản lượng sản phẩm nhựa toàn cầu, trong đó công nghệ đùn tấm phẳng và màng mỏng đóng góp khoảng 45% giá trị sản xuất công nghiệp. Tuy nhiên, một trong những thách thức kỹ thuật lớn nhất trong chế tạo khuôn đùn là hiện tượng phân bố dòng chảy không đồng đều tại cửa ra, khiến độ sai lệch chiều dày sản phẩm trong thực tế dao động từ 15% đến 25%, dẫn đến lãng phí nguyên liệu và suy giảm độ bền cơ học.

Vấn đề nghiên cứu trọng tâm của luận văn tập trung vào việc giải quyết bài toán thủy động lực học phi tuyến của chất lỏng đàn nhớt khi chuyển động qua các kênh dẫn phức tạp trong đầu đùn khe phẳng. Mục tiêu cụ thể là xây dựng hệ thống lý thuyết lưu biến học hoàn chỉnh cho polyme nóng chảy, thiết lập phương pháp giải tích thiết kế kênh phân phối dòng chảy có sức cản đồng đều và ứng dụng công nghệ mô phỏng số bằng phần tử hữu hạn để kiểm nghiệm, tối ưu hóa kết cấu hình học của đầu đùn.

Phạm vi nghiên cứu được thực hiện tại Trường Đại học Bách Khoa Hà Nội vào năm 2010, tập trung phân tích hai loại vật liệu nhiệt dẻo phổ biến trong công nghiệp là Polystyrene (PS) và High-Density Polyethylene (HDPE). Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm sai lệch vận tốc dòng chảy tại cửa ra xuống dưới mức 3%, giảm khoảng 20% chi phí thử nghiệm khuôn mẫu thực nghiệm và nâng cao hiệu suất sử dụng vật liệu thêm 12% cho các dây chuyền gia công đùn tấm phẳng.

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

Khung lý thuyết áp dụng

Luận văn xây dựng trên nền tảng cơ học môi trường liên tục kết hợp với lưu biến học chất lỏng phi Newton (Non-Newtonian Fluid Mechanics). Hai mô hình lý thuyết cốt lõi được ứng dụng bao gồm mô hình định luật lũy thừa (Power-Law Model) và mô hình phi tuyến Carreau-Yasuda. Mô hình Carreau-Yasuda cho phép mô tả chính xác sự suy giảm độ nhớt khi tốc độ trượt tăng cao (shear-thinning) trên toàn bộ dải biến thiên từ vùng tốc độ trượt bằng 0 đến vùng trượt vô cùng lớn, đồng thời bổ sung các định luật Herschel-Bulkley và chất lỏng Bingham cho các trường hợp dòng chảy có ứng suất trượt tới hạn.

Hệ thống phương trình bảo toàn cơ bản được thiết lập bao gồm phương trình bảo toàn khối lượng (phương trình liên tục), phương trình bảo toàn động lượng Navier-Stokes mở rộng cho môi trường đàn nhớt và phương trình bảo toàn năng lượng có xét đến hiệu ứng tiêu tán nhiệt do ma sát nhớt. Luận văn chuẩn hóa các khái niệm chuyên ngành then chốt như: độ nhớt cắt trượt ban đầu, gradient vận tốc, nhiệt độ thủy tinh hóa ($T_g$), hiện tượng giãn nở mặt cắt sau khi rời miệng khuôn (die swell hay hiệu ứng Barus) và hiện tượng gãy nứt dòng chảy polyme nóng chảy (melt fracture) khi vượt qua ngưỡng ứng suất trượt giới hạn.

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

Nguồn dữ liệu của đề tài được thu thập từ các thí nghiệm đo đạc đặc tính lưu biến thực nghiệm của nhựa PS và HDPE ở dải nhiệt độ gia công chuẩn từ 180°C đến 220°C. Cỡ mẫu nghiên cứu bao gồm 3 cấu hình hình học kênh tích liệu chính (kênh tích liệu trụ tròn thay đổi chiều sâu, kênh dạng khe chữ nhật và kênh dạng mắc áo - coat-hanger manifold) cùng hơn 50 kịch bản mô phỏng biến thiên các thông số hình học (với chiều sâu miệng tạo hình $t_0$ thay đổi ở các mức 198 mm, 218 mm và 245 mm). Phương pháp chọn mẫu phân tầng có chủ đích được áp dụng nhằm bao phủ toàn bộ các dạng đầu đùn được ứng dụng rộng rãi trong sản xuất công nghiệp tấm nhựa phẳng.

Phương pháp phân tích chủ đạo là sự kết hợp chặt chẽ giữa giải tích toán học giải hệ phương trình vi phân đạo hàm riêng và phương pháp mô phỏng số bằng phần tử hữu hạn (FEM) trên phần mềm chuyên dụng ANSYS Polyflow. Lý do lựa chọn Polyflow là nhờ khả năng xử lý vượt trội các bài toán cơ học chất lưu phi Newton phức tạp, hội tụ nhanh với các chất lỏng có tính đàn nhớt cao và phản ánh chính xác tương tác trường vận tốc - áp suất tại các góc lượn phân dòng với độ chính xác cao hơn khoảng 30% so với các phần mềm mô phỏng CFD thông thường. Toàn bộ tiến trình nghiên cứu được triển khai liên tục và hoàn thiện trong giai đoạn 2009 - 2010.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, kết quả giải tích và mô phỏng số chỉ ra rằng kết cấu kênh tích liệu dạng mắc áo (coat-hanger manifold) cho độ đồng đều phân bố dòng chảy vượt trội nhất. Độ lệch vận tốc của polyme nóng chảy tại cửa ra của đầu đùn mắc áo chỉ dao động trong khoảng 2,4% trên toàn bộ bề rộng tấm, trong khi đầu đùn kênh trụ tròn không hiệu chỉnh có mức sai lệch vận tốc lên đến 18,5%.

Thứ hai, việc tối ưu hóa chiều cao miệng tạo hình theo quy luật giải tích sức cản đồng đều giúp bù trừ triệt để sự tổn thất áp suất dọc theo kênh phân phối. Áp suất phân bố tại cửa vào giảm dần đều từ tâm ra hai biên mà không tạo ra các vùng áp suất cục bộ, giúp tổng độ sụt áp $\Delta p$ trên toàn bộ đầu đùn giảm khoảng 14,8%, từ 3,8 MPa xuống còn 3,24 MPa đối với vật liệu HDPE ở cùng điều kiện nhiệt độ 200°C.

Thứ ba, sự thay đổi tính chất lưu biến theo vật liệu có ảnh hưởng mạnh mẽ đến trường vận tốc. Khi gia tăng lưu lượng nạp từ 150 cm³/s lên 220 cm³/s, vật liệu HDPE có chỉ số Power-Law $n = 0,35$ thể hiện tính giả dẻo cao hơn so với PS ($n = 0,48$), dẫn đến độ nhớt tại vùng thành khuôn giảm mạnh hơn 32%, tạo ra profile vận tốc tại cửa ra có dạng phẳng đều hơn ở khu vực trung tâm.

Thảo luận kết quả

Nguyên nhân chính của hiện tượng mất cân bằng dòng chảy trong các đầu đùn khe phẳng thông thường là do hành trình dòng chảy từ tâm cấp liệu ra hai mép biên dài hơn nhiều so với hành trình đi thẳng ở chính giữa, kết hợp với tính phi tuyến của độ nhớt polyme. Đầu đùn dạng mắc áo khắc phục triệt để nhược điểm này bằng cách thu hẹp tiết diện kênh dẫn dọc theo phương ngang và điều chỉnh chiều sâu khe hở tạo hình, tạo ra một trường kháng lực thủy lực hoàn toàn cân bằng trên mọi đường dòng.

Khi so sánh với các công bố khoa học quốc tế trong cùng lĩnh vực, mô hình tính toán của luận văn đạt độ tương thích cao với độ sai lệch trường vận tốc dưới 4%, đồng thời đề xuất thêm bảng tham số hiệu chỉnh hình học linh hoạt theo từng chỉ số lưu biến cụ thể. Trên phương diện trực quan, các biểu đồ đường cong vận tốc (đặc biệt là đồ thị phân bố vận tốc tại cửa ra ứng với các kích thước $t_0 = 198\text{ mm}$, $218\text{ mm}$, $245\text{ mm}$) và các bảng dữ liệu lưu biến Carreau-Yasuda đã mô tả rõ nét quá trình chuyển đổi trạng thái của dòng chảy từ chế độ trượt chậm sang dòng chảy ổn định. Phát hiện này cung cấp căn cứ khoa học vững chắc giúp các kỹ sư cơ khí chế tạo làm chủ hoàn toàn công nghệ thiết kế khuôn đùn chính xác mà không cần dựa vào phương pháp thử nghiệm thủ công tốn kém.

Đề xuất và khuyến nghị

  1. Chuẩn hóa quy trình thiết kế khuôn đùn bằng công cụ số: Khuyến nghị các viện nghiên cứu và doanh nghiệp cơ khí áp dụng triệt để quy trình tính toán kích thước giải tích kết hợp mô phỏng Polyflow, đặt mục tiêu khống chế sai số chiều dày màng đùn dưới 2% trong vòng 6 tháng kể từ khi tiếp nhận chuyển giao công nghệ.
  2. Chuyển đổi hình học sang kênh tích liệu mắc áo: Đề xuất các nhà máy sản xuất tấm nhựa kỹ thuật thay thế kết cấu kênh chữ nhật cũ bằng kênh dạng mắc áo có biên dạng thay đổi liên tục, hướng tới mục tiêu gia tăng năng suất đùn thêm 15% đến 20% trong giai đoạn 2026 - 2027 do các kỹ sư trưởng bộ phận khuôn mẫu trực tiếp thực hiện.
  3. Kiểm soát nhiệt độ đa vùng trên thân khuôn: Thiết lập hệ thống điều khiển nhiệt độ độc lập với độ chính xác dao động không quá ±0,5°C tại từng phân đoạn dọc thân đầu đùn nhằm loại bỏ triệt để hiện tượng phân bố độ nhớt bất đối xứng, hoàn thành lắp đặt thử nghiệm trong thời hạn 12 tháng tại các cơ sở gia công chất dẻo.
  4. Xây dựng ngân hàng dữ liệu lưu biến vật liệu polyme: Khuyến nghị các phòng thí nghiệm trọng điểm thuộc các trường đại học kỹ thuật tiến hành đo đạc và số hóa thông số lưu biến của các dòng nhựa nguyên sinh và nhựa tái sinh thương dụng tại Việt Nam, nâng tỷ lệ khớp nối giữa mô phỏng và thực tế lên trên 95% trước quý IV năm 2027.

Đối tượng nên tham khảo luận văn

  • Kỹ sư thiết kế khuôn mẫu ngành nhựa: Tiếp cận trực tiếp công thức giải tích tính toán kích thước kênh dẫn và miệng tạo hình để thiết kế các bộ đầu đùn tấm phẳng khổ rộng từ 500 mm đến 1500 mm với độ chính xác cao.
  • Chuyên viên mô phỏng CAE/CFD: Nắm bắt quy trình phân chia lưới phần tử hữu hạn, phương pháp khai báo đặc tính chất lỏng phi Newton và thiết lập điều kiện biên phi tuyến phức tạp trên phần mềm Polyflow để tối ưu hóa thiết kế cơ khí.
  • Giảng viên và học viên cao học ngành Cơ học vật liệu: Sử dụng tài liệu như một giáo trình chuyên khảo mẫu mực về lý thuyết lưu biến học chất đàn nhớt, ứng dụng phương trình Navier-Stokes mở rộng và cơ học môi trường liên tục trong công nghiệp gia công polyme.
  • Giám đốc kỹ thuật và quản lý nhà máy bao bì nhựa: Hiểu rõ bản chất vật lý của các khuyết tật sản phẩm như phồng nở đầu đùn hay sóng bề mặt, từ đó điều chỉnh các thông số vận hành máy đùn (áp suất, nhiệt độ, tốc độ trục vít) nhằm tối ưu hóa chi phí sản xuất.

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

Tại sao kênh tích liệu dạng mắc áo lại ưu việt hơn các dạng kênh dẫn truyền thống?
Kênh dạng mắc áo được thiết kế với tiết diện thu hẹp dần từ tâm ra mép kết hợp góc nghiêng phân dòng hợp lý. Cấu trúc này cân bằng hoàn hảo hành trình dòng chảy và tổn thất áp suất trên mọi tiết diện, giúp vận tốc dòng polyme nóng chảy tại khe ra đạt độ đồng đều cao với sai lệch chỉ khoảng 2,4%.

Phần mềm ANSYS Polyflow đóng vai trò then chốt như thế nào trong nghiên cứu?
Polyflow là phần mềm chuyên biệt sử dụng phương pháp phần tử hữu hạn để mô phỏng dòng chảy đàn nhớt và phi Newton. Phần mềm giúp dự báo chính xác trường phân bố áp suất, gradient vận tốc và sự biến thiên độ nhớt trong kênh dẫn phức tạp, giảm hơn 70% thời gian chế tạo thử nghiệm khuôn mẫu.

Chỉ số Power-Law $n$ ảnh hưởng ra sao đến quá trình đùn polyme?
Chỉ số $n$ biểu thị mức độ phi Newton của vật liệu; khi $n < 1$, vật liệu có tính chất giả dẻo. Chỉ số $n$ càng nhỏ (như HDPE có $n = 0,35$ so với PS có $n = 0,48$) thì độ nhớt giảm càng mạnh tại vùng có tốc độ trượt cao ở thành khuôn, giúp dòng chảy thoát ra dễ dàng và phẳng đều hơn.

Làm cách nào để kiểm soát hiện tượng trương nở đầu đùn (die swell)?
Hiện tượng trương nở đầu đùn xuất hiện do tính đàn hồi tích tụ của các mạch đại phân tử khi bị biến dạng trượt trong khe hẹp. Có thể khắc phục bằng cách kéo dài chiều dài vùng ổn định định hình ở miệng khuôn, nâng nhiệt độ gia công thêm khoảng 5°C đến 10°C hoặc giảm tốc độ trượt qua khe.

Phương pháp tính toán trong luận văn có thể áp dụng cho nhựa tái sinh không?
Hoàn toàn có thể áp dụng được. Tuy nhiên, do nhựa tái sinh thường có độ phân tán phân tử lượng lớn và chứa tạp chất làm biến đổi chỉ số chảy, người thiết kế cần đo đạc lại các thông số lưu biến Carreau-Yasuda thực tế của lô nguyên liệu trước khi đưa vào mô hình mô phỏng số.

Kết luận

  • Luận văn đã hệ thống hóa toàn diện lý thuyết lưu biến học polyme nóng chảy và các mô hình chất lỏng phi Newton phục vụ thiết kế cơ khí chính xác.
  • Thiết lập thành công phương pháp giải tích xác định hình học đầu đùn khe phẳng có sức cản đồng đều, loại bỏ hiện tượng sai lệch lưu lượng dòng chảy.
  • Ứng dụng thành công phần mềm Polyflow để mô phỏng chi tiết trường áp suất và vận tốc của nhựa PS và HDPE trên 3 dạng kênh tích liệu khác nhau.
  • Chứng minh tính ưu việt tuyệt đối của đầu đùn dạng mắc áo với độ đồng đều dòng chảy đạt mức sai lệch dưới 2,4% và giảm sụt áp hơn 14%.
  • Đóng góp chính của công trình là đặt nền móng kết hợp giữa cơ học giải tích và mô phỏng số CAE trong thiết kế khuôn đùn chất lượng cao tại Việt Nam.

Giai đoạn 2026 - 2028 mở ra tiềm năng to lớn trong việc mở rộng mô hình sang bài toán nhiệt - thủy động lực học 3D phi đẳng nhiệt cho các dòng vật liệu polyme phân hủy sinh học thế hệ mới. Các kỹ sư và nhà nghiên cứu quan tâm nên khai thác triệt để các công thức giải tích và quy trình mô phỏng trong luận văn để nâng cao năng lực tự chủ thiết kế khuôn mẫu công nghệ cao.