Tổng quan luận án

Công nghệ bồi đắp vật liệu (Additive Manufacturing - AM), đặc biệt là công nghệ lắng đọng nóng chảy (Fused Deposition Modeling - FDM), chiếm thị phần lớn trong sản xuất công nghiệp nhờ ưu điểm sử dụng vật liệu thông dụng, kết cấu thiết bị đơn giản, chi phí đầu tư và bảo trì thấp do không cần nguồn năng lượng laser hay tia tử ngoại (UV). Tuy nhiên, khi ứng dụng FDM để chế tạo các sản phẩm kỹ thuật cơ khí, chất lượng sản phẩm vẫn còn bộc lộ nhiều hạn chế về độ chính xác kích thước và cơ tính (độ bền kéo). Các chỉ tiêu này phụ thuộc chặt chẽ vào việc lựa chọn và kiểm soát các thông số công nghệ trong quá trình tạo hình.

Khoảng trống nghiên cứu được luận án xác định xuất phát từ thực tế: phần lớn các nghiên cứu trước đây chủ yếu tập trung vào các loại vật liệu polymer nguyên sinh (như PLA, ABS, PEEK) hoặc chỉ khảo sát riêng lẻ một vài thông số mà chưa xem xét đầy đủ tác động đồng thời của cụm thông số công nghệ trọng yếu. Đặc biệt, vật liệu composite nền polymer pha trộn hạt kim loại (nhựa PLA-đồng) có tiềm năng ứng dụng cơ khí cao nhưng chưa được nghiên cứu toàn diện về mặt mô phỏng lưu nhiệt động học trong đầu đùn, quy hoạch thực nghiệm tối ưu hóa đa mục tiêu và mô hình hóa dự báo chất lượng.

Luận án tiến hành giải quyết các mục tiêu nghiên cứu cụ thể sau:

  1. Xác định ảnh hưởng của năm thông số công nghệ quan trọng (chiều dày lớp, góc đường đùn, tốc độ đầu đùn, mật độ điền đầy và nhiệt độ buồng tạo sản phẩm) đến độ chính xác kích thước và độ bền kéo của sản phẩm FDM trong điều kiện gia công cụ thể.
  2. Xây dựng mô hình thực nghiệm biểu diễn mối quan hệ định lượng giữa các thông số công nghệ đầu vào với độ chính xác kích thước (theo các phương X, Y, Z) và độ bền kéo.
  3. Xác định các bộ giá trị thông số công nghệ tối ưu cho từng mục tiêu đơn lẻ và tối ưu hóa đa mục tiêu nhằm đồng thời nâng cao độ chính xác kích thước và độ bền kéo của sản phẩm FDM.
  4. Xây dựng quy trình hỗ trợ quá trình kiểm định và hiệu chỉnh các thiết bị FDM được nghiên cứu, chế tạo tại Việt Nam.
  5. Đề xuất mô hình mạng nơ-ron nhân tạo (ANN) dự đoán độ chính xác kích thước và độ bền kéo dựa trên các thông số đầu vào với sai số thấp nhằm giảm thời gian và chi phí thực nghiệm.

Đối tượng và phạm vi nghiên cứu của luận án được xác định rõ:

  • Đối tượng nghiên cứu: Quá trình đùn ép tạo hình, độ chính xác kích thước theo ba phương tọa độ không gian ($X, Y, Z$) và độ bền kéo (chịu tải theo phương vuông góc với phương liên kết các lớp) của chi tiết chế tạo bằng công nghệ FDM.
  • Vật liệu nghiên cứu: Vật liệu composite polymer kết hợp kim loại PLA-đồng (do hãng Formfutura sản xuất).
  • Phạm vi thông số khảo sát: Khảo sát 5 thông số công nghệ biến thiên gồm: chiều dày lớp ($t = 0{,}1 \div 0{,}3\text{ mm}$), góc đường đùn ($\theta = 0^\circ \div 90^\circ$), tốc độ đầu đùn ($V_p = 40 \div 60\text{ mm/s}$), mật độ điền đầy ($D = 20% \div 60%$) và nhiệt độ buồng tạo sản phẩm ($T_b = 50^\circ\text{C} \div 60^\circ\text{C}$). Các thông số được cố định gồm: nhiệt độ đầu đùn ($190^\circ\text{C}$), bề rộng đường đùn ($0{,}4\text{ mm}$), số đường bao ngoài ($3\text{ đường}$), kiểu điền đầy dạng đường thẳng song song (Rectilinear) và hướng chế tạo theo phương mặt phẳng.
  • Phạm vi thiết bị: Thiết bị FDM công nghiệp FDM-Vina2015 do Trường Đại học Bách khoa – ĐHQG-HCM chế tạo, máy kiểm tra cơ tính vạn năng DELTALAB và thước cặp điện tử Mitutoyo.

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

Luận án đã phân tích và tổng hợp có hệ thống các công trình nghiên cứu trong và ngoài nước liên quan đến các yếu tố ảnh hưởng và phương pháp tối ưu hóa chất lượng sản phẩm FDM:

Các hướng nghiên cứu quốc tế

  • Ảnh hưởng đến độ chính xác kích thước: Oguzhan Emre Akbas (2019, Thổ Nhĩ Kỳ) phân tích ảnh hưởng của nhiệt độ đầu đùn và tốc độ cấp liệu đến sai số kích thước mẫu PLA và ABS; Beibit Zharylkassyn và cộng sự (2021, Kazakhstan) kết luận chiều dày lớp $0{,}1 \div 0{,}2\text{ mm}$ phù hợp cho ABS/PLA và hướng chế tạo $0^\circ, 90^\circ$ cho độ chính xác cao nhất; D. Murthy (2018, Ấn Độ) khảo sát chiều dày lớp, góc đường đùn và hướng chế tạo trên mẫu ABS và PLA; A. Mohanty và cộng sự (2021, Ấn Độ) đánh giá 5 thông số và xác định hướng chế tạo là nhân tố chi phối lớn nhất đến độ chính xác kích thước; K. Manolakos (2020, Ai Cập) chỉ ra nhiệt độ đầu đùn là thông số ảnh hưởng mạnh nhất đến dung sai hình học.
  • Ảnh hưởng đến cơ tính sản phẩm: Fuda Ning và cộng sự (2016, Đại học Kỹ thuật Texas, Hoa Kỳ) thử nghiệm mẫu kéo theo tiêu chuẩn ASTM D638-10, kết luận góc đùn $0^\circ$ và $90^\circ$, tốc độ $25\text{ mm/s}$, nhiệt độ đùn $220^\circ\text{C}$ và chiều dày lớp $0{,}15\text{ mm}$ mang lại độ bền kéo cực đại; Shilpesh R. Rajpurohit và Harshit K. Dave (2018, Ấn Độ) xác định góc đùn $0^\circ$ và giảm chiều dày lớp giúp tăng độ bền kéo của PLA (ASTM D638 loại I); Gebisa và Lemu (2019, Na Uy) sử dụng máy Fortus 450Mc và tiêu chuẩn ASTM D3039, chỉ ra góc đường đùn là nhân tố then chốt; Deng và cộng sự (2018, Trung Quốc) khảo sát vật liệu PEEK với độ bền kéo tối ưu tại tốc độ $60\text{ mm/s}$, chiều dày lớp $0{,}2\text{ mm}$, nhiệt độ đùn $370^\circ\text{C}$ và mật độ điền đầy $40%$; Dinesh Yadav và cộng sự (2019, Ấn Độ) nghiên cứu vật liệu PETG, ABS và hỗn hợp ABS/PETG theo ASTM D638-IV; M. Heidari-Rarani và cộng sự (2020, Iran) tối ưu hóa cơ tính PLA với mật độ điền đầy $80%$, tốc độ $40\text{ mm/s}$, chiều dày lớp $0{,}1\text{ mm}$.
  • Tối ưu hóa đa mục tiêu và vật liệu composite: J. Chacón và cộng sự (2017, 2019, Tây Ban Nha) phân tích vật liệu composite CFRTPCs và PLA bằng kính hiển vi điện tử quét SEM và ANOVA; Lu Wang và cộng sự (2018, Hoa Kỳ) so sánh mẫu PP giữa FDM và ép phun; Omar Ahmed Mohamed và cộng sự (2016, Úc) nghiên cứu mẫu composite PC/ABS pha kim loại và ứng dụng thiết kế tối ưu Q-optimal; Shuheng Wang và cộng sự (2020, Trung Quốc) thử nghiệm tiêu chuẩn ISO 527-2:2012; Wang và cộng sự (2019, Đại học Sơn Đông, Trung Quốc) ứng dụng mô phỏng phần tử hữu hạn (FEA) khảo sát dòng chảy lưu biến; Carmita Camposeco Negrete (2020, Mexico) tối ưu 5 mục tiêu trên máy Stratasys F270; S. Jiang và cộng sự (2018, Trung Quốc) khảo sát nhiệt độ đùn lên tới $370^\circ\text{C}$ trên máy Makerpi M2030; Tambrallimath (2021, Ấn Độ) nghiên cứu ABS gia cố nano carbon; Nikzad và cộng sự (Đại học Kỹ thuật Swinburne, Úc) mô phỏng dòng chảy 2D/3D của hỗn hợp ABS-sắt; Sa'ude (Malaysia) và Omar Ahmed Mohamed (2016) phát triển vật liệu nhựa nhiệt dẻo pha bột kim loại ABS-Iron và ABS-Copper.
  • Nhiệt độ buồng tạo sản phẩm: Kou và cộng sự (Đài Loan) chỉ ra nhiệt độ bàn máy và nhiệt độ buồng ảnh hưởng lần lượt $54{,}18%$ và $39{,}87%$ đến biến dạng cong vênh của ABS; Wichniarek và cộng sự (Ba Lan) khảo sát cơ tính ABS; Casavola và cộng sự (Ý) chứng minh buồng nhiệt của Stratasys giúp giảm ứng suất dư so với MakerBot; Wang và cộng sự (Trung Quốc) lập mô hình co rút vật liệu; Wu và cộng sự (Trung Quốc) xác định buồng nhiệt $130^\circ\text{C}$ hạn chế biến dạng của PEEK.

Các nghiên cứu trong nước

Đoàn Yên Thế (2019, Đại học Thủy lợi) ứng dụng phương pháp Taguchi xác định thông số tối ưu giảm sai số kích thước mẫu PLA; Phan Quang Thế (2018, Đại học Kỹ thuật Thái Nguyên) khảo sát khe hở đường đùn đến độ bền chi tiết nhựa PC trên máy Stratasys Fortus 400MC; Đề tài cấp nhà nước do PGS.TS Đặng Văn Nghìn chủ nhiệm ứng dụng tạo mẫu nhanh trong công nghiệp và y tế (chế tạo mảnh ghép sọ não cho 25 bệnh nhân); Đề tài do PGS.TS Thái Thị Thu Hà chủ nhiệm chế tạo máy FDM công nghiệp FDM-Vina2015; GS.TS Nguyễn Thanh Nam (2019) nghiên cứu máy in 3D theo mô-đun; PGS.TS Bùi Trọng Hiếu (2020) nghiên cứu bộ điều khiển nhiệt độ buồng máy in 3D.

Kỹ thuật và phương pháp tối ưu Độ chính xác Tối ưu đa mục tiêu Yêu cầu mô hình toán Thời gian tính toán Phương pháp tiếp cận Mô hình tuyến tính / phi tuyến Phạm vi dữ liệu Khả năng tương tác giữa các biến Độ khó thực hiện
Giải thuật di truyền (GA) Cao Có thể Không Rất dài Đơn giản Cả hai Cao Không thể Rất khó
Mạng nơ-ron nhân tạo (ANN) Rất cao Có thể Có Rất dài Thông qua mô phỏng Cả hai Cao Không thể Trung bình
Thiết kế giai thừa Trung bình Không thể Có Ngắn Đơn giản Tuyến tính Trung bình Có thể Dễ
Một nhân tố tại một thời điểm (OFAT) Thấp Không thể Không Rất dài Đơn giản Tuyến tính Cao Không thể Dễ
Lôgic mờ Cao Có thể Có Dài Thông qua mô phỏng Cả hai Cao Không thể Rất khó
Phương pháp Taguchi Thấp Không thể Không Ngắn Đơn giản Tuyến tính Trung bình Có thể Khó
Phân tích quan hệ xám Trung bình Có thể Không Ngắn Đơn giản Tuyến tính Trung bình Có thể Khó
Bề mặt đáp ứng (RSM) Rất cao Có thể Có Ngắn Thông qua mô phỏng Cả hai Thấp Có thể Trung bình
Xử lý dữ liệu nhóm (GMDH) Cao Có thể Có Trung bình Thông qua mô phỏng Cả hai Cao Không thể Trung bình
Tối ưu bầy đàn (PSO) Rất cao Có thể Có Rất dài Thông qua mô phỏng Cả hai Cao Không thể Rất khó

Vị trí và khoảng trống giải quyết của luận án

Luận án kết hợp mô phỏng thủy nhiệt động lực học trên Ansys Fluent nhằm xác định vận tốc cấp liệu an toàn dựa trên giới hạn uốn dọc Euler của sợi vật liệu, sau đó sử dụng quy hoạch thực nghiệm theo mô hình phối hợp tâm hướng mặt (FCCCD) kết hợp phân tích phương sai (ANOVA), thuật toán tối ưu hóa bề mặt đáp ứng (RSM) và mô hình học máy mạng nơ-ron nhân tạo (ANN) trên vật liệu mới PLA-đồng.

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

Khung lý thuyết và mô hình dòng chảy

Quá trình tạo hình FDM dựa trên nguyên lý đùn ép lớp vật liệu nóng chảy qua đầu đùn có vòi phun đường kính nhỏ ($d = 0{,}4\text{ mm}$). Khi gia nhiệt đến trạng thái nóng chảy dẻo ($190^\circ\text{C}$), vật liệu composite PLA-đồng biểu hiện tính chất của chất lỏng phi Newton. Độ nhớt biểu kiến $\eta$ của dòng chảy được mô tả qua định luật hàm số mũ (Power-law):

$$\eta = K \cdot \dot{\gamma}^{n-1}$$

Trong đó:

  • $\eta$: Độ nhớt của lưu chất ($\text{Pa}\cdot\text{s}$)
  • $K$: Hệ số đồng nhất dòng chảy ($3000\text{ kg}\cdot\text{s}^{n-2}/\text{m}$)
  • $\dot{\gamma}$: Tốc độ trượt ($\text{s}^{-1}$)
  • $n$: Chỉ số luật lũy thừa / hệ số công suất ($n = 0{,}77$)

Sợi vật liệu ở trạng thái rắn được đẩy vào đầu đùn bởi cặp con lăn dẫn động. Để tránh hiện tượng sợi nhựa bị uốn cong, mất ổn định hoặc gãy trước khi đi vào vùng hóa lỏng, áp lực đẩy cấp liệu phải nhỏ hơn lực tới hạn nén uốn dọc Euler $P_{cr}$:

$$P_{cr} = \frac{\pi^2 E I}{(\mu L)^2}$$

Trong đó $E$ là mô đun đàn hồi của vật liệu, $I$ là mômen quán tính tiết diện sợi dây, $L$ là chiều dài tự do của đoạn sợi chịu nén và $\mu$ là hệ số liên kết chiều dài.

Thông số nhiệt - cơ của vật liệu PLA-đồng Giá trị định lượng Đơn vị đo
Khối lượng riêng ($\rho$) $4000$ $\text{kg/m}^3$
Nhiệt dung riêng ($C_p$) $706{,}424$ $\text{J/kg}\cdot\text{K}$
Hệ số dẫn nhiệt ($k$) $0{,}381$ $\text{W/m}\cdot\text{K}$
Nhiệt độ nóng chảy ($T_m$) $433$ $\text{K}$
Hệ số định luật công suất ($K$) $3000$ $\text{kg}\cdot\text{s}^{n-2}/\text{m}$
Hệ số công suất ($n$) $0{,}77$ Không thứ nguyên
Mô đun đàn hồi ($E$) $93{,}77$ $\text{MPa}$

Phương pháp nghiên cứu và công cụ xử lý

  • Phương pháp mô phỏng số: Sử dụng phần mềm Ansys Fluent (phiên bản 2011) tích hợp công cụ chia lưới Ansys Meshing và hình học DesignModeler / SpaceClaim để giải hệ phương trình bảo toàn khối lượng, bảo toàn động lượng Navier-Stokes và phương trình năng lượng truyền nhiệt của dòng chảy phi Newton trong lòng vòi phun.
  • Quy hoạch thực nghiệm: Ứng dụng phương pháp thiết kế tổng hợp trung tâm hướng mặt (Face-Centered Central Composite Design - FCCCD) thuộc phương pháp bề mặt đáp ứng (RSM) trên phần mềm Minitab. Cấu trúc FCCCD gồm các điểm nhân tố bậc $2^k$, các điểm trục (với khoảng cách trục $\alpha = 1{,}0$) và các điểm tâm lặp nhằm đánh giá mô hình hồi quy phi tuyến bậc hai:

$$Y = \beta_0 + \sum_{i=1}^{k} \beta_i X_i + \sum_{i=1}^{k} \beta_{ii} X_i^2 + \sum_{i<j} \beta_{ij} X_i X_j + \varepsilon$$

  • Phân tích thống kê: Sử dụng phân tích phương sai (ANOVA) ở mức ý nghĩa $\alpha = 0{,}05$ (độ tin cậy $95%$) để đánh giá hệ số xác định ($R^2$, $R^2_{\text{adj}}$), giá trị kiểm định chuẩn F ($F$-value) và giá trị xác suất $p$-value của từng thông số đơn lẻ cũng như các tương tác chéo.
  • Mô hình học máy dự báo: Xây dựng mạng nơ-ron truyền thẳng nhiều lớp (Artificial Neural Network - ANN) với thuật toán lan truyền ngược sai số (Backpropagation) gồm 1 lớp vào (5 nút), 1 lớp ẩn và 1 lớp ra (tương ứng với từng chỉ tiêu sai số $X, Y, Z$ hoặc độ bền kéo). Đánh giá mô hình mạng qua sai số toàn phương trung bình (Mean Squared Error - MSE) và hệ số tương quan ($R$).
  • Nguồn dữ liệu và thiết bị đo kiểm: Mẫu kiểm tra độ bền kéo được gia công theo tiêu chuẩn và thử nghiệm trên thiết bị kiểm tra cơ tính DELTALAB; kích thước hình học theo ba trục được đo bằng thước cặp điện tử Mitutoyo có độ chính xác $0{,}01\text{ mm}$.

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

Chương 1: Tổng quan

Chương 1 hệ thống hóa bức tranh toàn cảnh về công nghệ AM và FDM, chỉ ra tốc độ tăng trưởng kinh tế và tỷ trọng sử dụng trong công nghiệp cơ khí, khuôn mẫu và chế tạo chi tiết hoàn thiện. Chương này khảo sát chi tiết hơn 50 công trình nghiên cứu tiêu biểu trên thế giới và tại Việt Nam, phân loại rõ các nhóm ảnh hưởng chính: chiều dày lớp, góc quét sợi, mật độ điền đầy, tốc độ đầu đùn, nhiệt độ đùn và nhiệt độ buồng máy. Trên cơ sở đối chiếu ưu nhược điểm của các phương pháp thực nghiệm (Taguchi, GA, ANN, RSM, OFAT), chương 1 xác lập 5 mục tiêu cụ thể, phạm vi nghiên cứu thực nghiệm trên vật liệu composite PLA-đồng, phương pháp tiếp cận mô phỏng kết hợp quy hoạch thực nghiệm và ý nghĩa khoa học - thực tiễn của đề tài.

Chương 2: Ảnh hưởng của thông số công nghệ đến chất lượng sản phẩm FDM

Chương 2 định nghĩa bản chất vật lý và cơ chế ảnh hưởng của 9 thông số công nghệ trong quy trình FDM đến độ chính xác và cơ tính:

  1. Chiều dày lớp ($t$): Quyết định diện tích tiếp xúc giữa các lớp và hiệu ứng bậc thang trên bề mặt; mức thực nghiệm chọn: thấp ($-1$) là $0{,}1\text{ mm}$, trung bình ($0$) là $0{,}2\text{ mm}$, cao ($+1$) là $0{,}3\text{ mm}$.
  2. Góc đường đùn ($\theta$): Hướng xếp sợi nhựa so với trục tọa độ, quyết định khả năng truyền lực kéo; mức thực nghiệm chọn: thấp ($-1$) là $0^\circ$, trung bình ($0$) là $45^\circ$, cao ($+1$) là $90^\circ$.
  3. Tốc độ đầu đùn ($V_p$): Tốc độ dịch chuyển vòi phun, tương thích với vận tốc cấp liệu; mức thực nghiệm chọn: thấp ($-1$) là $40\text{ mm/s}$, trung bình ($0$) là $50\text{ mm/s}$, cao ($+1$) là $60\text{ mm/s}$.
  4. Mật độ điền đầy ($D$): Tỷ lệ phần trăm thể tích điền đầy bên trong chi tiết; mức thực nghiệm chọn: thấp ($-1$) là $20%$, trung bình ($0$) là $40%$, cao ($+1$) là $60%$.
  5. Nhiệt độ buồng tạo sản phẩm ($T_b$): Môi trường nhiệt xung quanh giúp giảm gradient nhiệt độ, hạn chế nguội đột ngột và co ngót; mức thực nghiệm chọn: thấp ($-1$) là $50^\circ\text{C}$, trung bình ($0$) là $55^\circ\text{C}$, cao ($+1$) là $60^\circ\text{C}$.
  6. Hướng chế tạo: Cố định theo phương mặt phẳng (flat build orientation) nhằm giảm số lượng lớp đắp, tối ưu hóa thời gian và diện tích tựa nền.
  7. Kiểu điền đầy: Cố định dạng Rectilinear (các đường song song đan chéo) giúp tiết kiệm vật liệu và thời gian chế tạo.
  8. Nhiệt độ đầu đùn: Cố định ở $190^\circ\text{C}$ đảm bảo vật liệu PLA-đồng nóng chảy hoàn toàn mà không bị thoái hóa nhiệt hay rỗ khí.
  9. Bề rộng đường đùn và số đường bao ngoài: Cố định bề rộng đùn là $0{,}4\text{ mm}$ và số đường bao ngoài là $3\text{ đường}$ để ổn định biên dạng ngoài.

Chương 3: Mô phỏng quá trình đùn nhựa

Chương 3 xây dựng bài toán mô phỏng phần tử hữu hạn dòng lưu chất truyền nhiệt trong cụm đầu đùn bằng Ansys Fluent 2011. Các bước thực hiện bao gồm thiết lập mô hình 3D hình học vòi phun, chia lưới phần tử, định nghĩa thông số nhiệt - vật lý của PLA-đồng và giải phương trình Navier-Stokes cho dòng phi Newton Power-law.

Kết quả mô phỏng làm rõ:

  • Trường phân bố vận tốc, tốc độ trượt và ứng suất cắt lớn nhất tập trung tại sát thành vòi phun.
  • Mối tương quan giữa vận tốc cấp liệu ($V_e$) với chiều dài vùng nóng chảy hoàn toàn và độ giảm áp suất trong ống đùn.
  • Xác định ngưỡng áp suất đùn cực đại để áp lực đẩy không vượt quá giới hạn bền nén mất ổn định Euler ($P_{cr}$) của sợi PLA-đồng.
  • Kết quả mô phỏng làm cơ sở khoa học để xác lập khoảng vận tốc đầu đùn an toàn và hợp lý đưa vào thực nghiệm là $40 \div 60\text{ mm/s}$.

Chương 4: Thực nghiệm và phân tích

Chương 4 triển khai ma trận thực nghiệm FCCCD với 5 thông số ($t, \theta, V_p, D, T_b$) trên máy in FDM-Vina2015. Các mẫu kiểm tra sai số kích thước ($\Delta X, \Delta Y, \Delta Z$) và mẫu kéo ASTM D638 được chế tạo và đo lường đồng bộ.

Thông số công nghệ Ký hiệu Đơn vị Mức thấp (-1) Mức trung bình (0) Mức cao (+1)
Chiều dày lớp $t$ $\text{mm}$ $0{,}1$ $0{,}2$ $0{,}3$
Góc đường đùn $\theta$ Độ ($^\circ$) $0$ $45$ $90$
Tốc độ đầu đùn $V_p$ $\text{mm/s}$ $40$ $50$ $60$
Mật độ điền đầy $D$ $%$ $20$ $40$ $60$
Nhiệt độ buồng tạo sản phẩm $T_b$ $^\circ\text{C}$ $50$ $55$ $60$

Phân tích phương sai ANOVA chỉ ra:

  • Sai số kích thước theo phương $X$ và $Y$ chịu ảnh hưởng mạnh của chiều dày lớp, nhiệt độ buồng và tương tác giữa tốc độ đùn với mật độ điền đầy.
  • Sai số kích thước theo phương $Z$ phụ thuộc chủ yếu vào chiều dày lớp $t$ do tích lũy số lượng lớp đắp.
  • Độ bền kéo chịu chi phối lớn nhất bởi mật độ điền đầy $D$, tiếp theo là góc đường đùn $\theta$ và nhiệt độ buồng $T_b$.
  • Xây dựng mô hình ANN với 1 lớp ẩn cho thấy khả năng xấp xỉ dữ liệu phi tuyến rất cao; hệ số tương quan $R$ của các mô hình dự báo kích thước và độ bền kéo đều đạt mức cao, sai số MSE ở mức nhỏ, khẳng định tính tương thích tốt giữa mô hình thực nghiệm RSM và mạng ANN.

Chương 5: Tối ưu thông số công nghệ

Chương 5 ứng dụng thuật toán tối ưu hóa hàm mục tiêu đáp ứng (Desirability Function) trên phần mềm Minitab để xác định các bộ thông số công nghệ:

  1. Bộ thông số tối ưu cho độ chính xác kích thước: Cân bằng sai lệch $\Delta X, \Delta Y, \Delta Z$ nhỏ nhất.
  2. Bộ thông số tối ưu cho độ bền kéo: Tối đa hóa giá trị ứng suất kéo cực đại của sản phẩm.
  3. Bộ thông số tối ưu đa mục tiêu: Dung hòa đồng thời sai số kích thước tối thiểu và độ bền kéo tối đa.

Tác giả đã tiến hành chế tạo lại các mẫu thử độc lập từ các bộ thông số tối ưu tìm được để kiểm chứng. Kết quả kiểm nghiệm thực tế cho thấy sai số giữa giá trị dự đoán lý thuyết và giá trị đo thực nghiệm nằm trong giới hạn cho phép, xác thực độ tin cậy của mô hình tối ưu.

Chương 6: Kết luận – Đề xuất

Chương 6 tổng hợp các kết quả khoa học đạt được, khẳng định việc làm chủ công nghệ kiểm soát thông số FDM trên vật liệu PLA-đồng và đưa ra các đề xuất nghiên cứu mở rộng trong tương lai.

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

  • Thiết lập thành công mô hình toán học giải tích kết hợp mô phỏng số (Ansys Fluent) dòng chảy phi Newton Power-law của vật liệu composite PLA-đồng trong đầu đùn FDM, xác định được quan hệ giữa tốc độ cấp liệu và giới hạn bền uốn dọc Euler của sợi cấp liệu.
  • Xây dựng hệ phương trình hồi quy bề mặt đáp ứng bậc hai (RSM/FCCCD) biểu diễn mối quan hệ định lượng giữa 5 thông số công nghệ trọng yếu ($t, \theta, V_p, D, T_b$) với sai số kích thước không gian ($\Delta X, \Delta Y, \Delta Z$) và độ bền kéo của chi tiết FDM bằng vật liệu PLA-đồng.
  • Xây dựng và kiểm chứng mô hình học máy mạng nơ-ron nhân tạo (ANN) có khả năng dự đoán chính xác chất lượng sản phẩm đầu ra, làm phong phú thêm cơ sở lý thuyết tối ưu hóa công nghệ bồi đắp vật liệu composite nền polymer.

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

  • Cung cấp các bộ thông số công nghệ tối ưu hóa đơn mục tiêu và đa mục tiêu đã qua kiểm chứng thực nghiệm, hỗ trợ trực tiếp cho các cơ sở sản xuất trong việc cài đặt thông số gia công chi tiết composite PLA-đồng đạt cơ tính và độ chính xác cao.
  • Xây dựng quy trình thực nghiệm và kiểm soát thông số công nghệ chuẩn mực, đóng vai trò là tài liệu kỹ thuật phục vụ việc hiệu chuẩn, kiểm định và đánh giá tính năng làm việc của các dòng máy in 3D FDM công nghiệp được nghiên cứu chế tạo trong nước (như FDM-Vina2015).

Kiến nghị và giải pháp đề xuất

  • Áp dụng các bộ thông số công nghệ tối ưu tương ứng với từng yêu cầu làm việc cụ thể của chi tiết: ưu tiên lắp ráp chính xác hoặc ưu tiên khả năng chịu tải.
  • Sử dụng phần mềm dự báo dựa trên mạng ANN đã huấn luyện để giảm thiểu thời gian thử nghiệm và lãng phí phôi liệu khi đưa vật liệu mới vào sản xuất.

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

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

  • Nghiên cứu tập trung khảo sát trên một loại vật liệu composite cụ thể là PLA-đồng của hãng Formfutura.
  • Cố định hướng chế tạo theo phương mặt phẳng và một số thông số phụ (nhiệt độ đầu đùn $190^\circ\text{C}$, bề rộng đùn $0{,}4\text{ mm}$, kiểu điền đầy Rectilinear) nhằm tập trung vào 5 thông số biến thiên chính.
  • Kiểm tra cơ tính mới tập trung vào độ bền kéo đơn trục vuông góc với phương liên kết lớp, chưa khảo sát các trạng thái tải trọng phức tạp khác như uốn, nén, va đập hoặc độ bền mỏi.

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

  • Mở rộng phạm vi nghiên cứu sang các loại vật liệu composite nền nhựa pha trộn sợi carbon, bột sắt hoặc nhựa chịu nhiệt cao (PEEK, ULTEM).
  • Khảo sát các hướng đặt mẫu nghiêng không gian và các kiểu điền đầy phức tạp (tam giác, tổ ong, lượn sóng) kết hợp phân tích vi cấu trúc đứt gãy bằng kính hiển vi điện tử SEM.
  • Nghiên cứu ứng dụng các thuật toán tối ưu hóa thông minh khác như giải thuật di truyền (GA) hoặc tối ưu bầy đàn (PSO) kết hợp điều khiển nhiệt độ buồng động học trong thời gian thực.
  • Mở rộng phương pháp luận tối ưu hóa sang các công nghệ AM khác như SLA, DLP, SLM.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho:

  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Kỹ thuật Cơ khí, Chế tạo máy, Vật liệu: Tham khảo phương pháp luận quy hoạch thực nghiệm FCCCD, phân tích phương sai ANOVA, mô phỏng lưu nhiệt động học dòng chảy phi Newton bằng Ansys Fluent và kỹ thuật huấn luyện mạng ANN dự báo chất lượng chế tạo.
  • Kỹ sư thiết kế, vận hành công nghệ in 3D trong công nghiệp: Ứng dụng trực tiếp các bộ thông số tối ưu để gia công chi tiết máy, đồ gá, chi tiết composite PLA-đồng đòi hỏi độ chính xác kích thước và độ bền kéo.
  • Các nhóm nghiên cứu và doanh nghiệp chế tạo thiết bị cơ khí chính xác: Sử dụng quy trình của luận án làm tài liệu tham khảo cho công tác kiểm định, đánh giá và hiệu chuẩn thiết bị FDM công nghiệp sản xuất tại Việt Nam.

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

1. Vì sao luận án lựa chọn vật liệu composite PLA-đồng thay vì các vật liệu truyền thống như PLA hay ABS nguyên sinh?

Vật liệu composite PLA-đồng (nền nhựa PLA pha trộn hạt đồng) có ưu điểm vượt trội về khả năng dẫn nhiệt, độ ổn định kích thước, ít cong vênh, cải thiện dòng chảy và độ bám dính giữa các lớp, đồng thời rất thuận lợi cho các công đoạn hậu xử lý cơ khí (cắt, mài, đánh bóng). Đây là vật liệu có tiềm năng ứng dụng cao trong chế tạo chi tiết công nghiệp nhưng chưa được nghiên cứu toàn diện về quy luật công nghệ FDM.

2. Mô hình dòng chảy nào được sử dụng để mô phỏng vật liệu PLA-đồng trong đầu đùn?

Luận án sử dụng mô hình chất lỏng phi Newton theo định luật công suất Power-law ($\eta = K \cdot \dot{\gamma}^{n-1}$) với hệ số đồng nhất $K = 3000\text{ kg}\cdot\text{s}^{n-2}/\text{m}$ và chỉ số công suất $n = 0{,}77$, kết hợp các phương trình bảo toàn khối lượng, động lượng và năng lượng trong phần mềm Ansys Fluent 2011.

3. Năm thông số công nghệ biến thiên được lựa chọn trong quy hoạch thực nghiệm FCCCD là gì?

Năm thông số gồm: (1) Chiều dày lớp $t$ ($0{,}1 \div 0{,}3\text{ mm}$); (2) Góc đường đùn $\theta$ ($0^\circ \div 90^\circ$); (3) Tốc độ đầu đùn $V_p$ ($40 \div 60\text{ mm/s}$); (4) Mật độ điền đầy $D$ ($20% \div 60%$); (5) Nhiệt độ buồng tạo sản phẩm $T_b$ ($50^\circ\text{C} \div 60^\circ\text{C}$).

4. Tốc độ cấp liệu an toàn của đầu đùn được kiểm soát dựa trên điều kiện cơ học nào?

Tốc độ cấp liệu được kiểm soát sao cho áp suất đùn ngược tác động lên đoạn sợi nhựa chưa nóng chảy không vượt quá tải trọng nén tới hạn gây mất ổn định uốn dọc theo công thức Euler ($P_{cr} = \frac{\pi^2 E I}{(\mu L)^2}$), ngăn ngừa hiện tượng cong gãy sợi vật liệu trước khi đi vào buồng hóa lỏng.

5. Kết quả so sánh giữa hai phương pháp tối ưu hóa RSM và mô hình ANN trong luận án như thế nào?

Cả hai phương pháp đều cho kết quả dự báo rất sát với thực nghiệm. Phương pháp RSM cung cấp phương trình hồi quy rõ ràng và xác định các bộ thông số tối ưu toàn cục thông qua hàm đáp ứng, trong khi mạng nơ-ron nhân tạo ANN có khả năng học và dự báo các quan hệ phi tuyến phức tạp với sai số bình phương trung bình (MSE) rất nhỏ và hệ số tương quan $R$ ở mức cao.

Kết luận

Luận án của tác giả Huỳnh Hữu Nghị đã giải quyết trọn vẹn bài toán nâng cao chất lượng sản phẩm công nghệ FDM trên vật liệu composite PLA-đồng thông qua sự kết hợp chặt chẽ giữa mô phỏng thủy nhiệt Ansys Fluent, quy hoạch thực nghiệm FCCCD, phân tích phương sai ANOVA, tối ưu hóa bề mặt đáp ứng RSM và mô hình dự báo học máy ANN. Công trình đã xác định chính xác quy luật ảnh hưởng của 5 thông số công nghệ then chốt đến độ chính xác kích thước 3D và độ bền kéo của sản phẩm, đồng thời xác lập các bộ thông số tối ưu phục vụ chế tạo thực tế. Kết quả nghiên cứu không chỉ đóng góp cơ sở khoa học giá trị cho ngành kỹ thuật cơ khí chế tạo bồi đắp mà còn mang ý nghĩa thực tiễn quan trọng trong việc làm chủ quy trình kiểm định và vận hành thiết bị FDM công nghiệp tại Việt Nam.