Tổng quan luận án

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

Trong lĩnh vực chế tạo máy hiện đại, áp lực cạnh tranh công nghiệp đòi hỏi các doanh nghiệp phải rút ngắn thời gian và chi phí sản xuất, đồng thời đáp ứng nhu cầu thay đổi nhanh chóng về mẫu mã và kích thước sản phẩm. Công nghệ gia công biến dạng kim loại tấm truyền thống như dập, vuốt, ép sâu đòi hỏi các bộ khuôn kim loại gia công cơ khí chính xác có chi phí chế tạo cao và thời gian phát triển dài, do đó không mang lại hiệu quả kinh tế đối với các sản phẩm đơn chiếc, loạt nhỏ hoặc các mẫu thử nghiệm nhanh trong công nghiệp quốc phòng, hàng không vũ trụ, ô tô và y học.

Công nghệ tạo hình gia tăng (Incremental Sheet Forming - ISF) ra đời nhằm khắc phục hạn chế này bằng cách sử dụng đầu dụng cụ không có lưỡi cắt di chuyển trên máy điều khiển số (CNC) để làm biến dạng phôi tấm từng lớp. Tuy nhiên, phương pháp tạo hình gia tăng đơn điểm (Single Point Incremental Forming - SPIF) bộc lộ nhược điểm lớn về độ chính xác hình học, hiện tượng biến dạng uốn không mong muốn ngoài vùng gia công, hiện tượng phục hồi đàn hồi (springback), và không thể tạo ra các chi tiết có đồng thời cả vùng lồi và vùng lõm.

Phương pháp tạo hình gia tăng đa điểm (Two Point Incremental Forming - TPIF) sử dụng thêm đồ gá chuyên dụng và dưỡng tạo hình cố định bán phần giúp cải thiện độ cứng vững, nâng cao độ chính xác kích thước và mở rộng phạm vi hình học chi tiết. Tại Việt Nam, các công trình trước đây chủ yếu tập trung vào công nghệ SPIF, thiếu các nghiên cứu hệ thống về thiết kế đồ gá chuyên dụng TPIF, xây dựng mô hình thử nghiệm tối ưu hóa góc biến dạng, mô phỏng số dự đoán phá hủy vật liệu và phân tích ảnh hưởng của chế độ bôi trơn trong TPIF. Luận án của tác giả Ma Văn Việt lựa chọn giải quyết khoảng trống này.

Mục tiêu và câu hỏi nghiên cứu

Mục tiêu chung của luận án là nghiên cứu khả năng tạo hình của tấm hợp kim nhôm A1050 H14 có độ dày 1,5 mm bằng công nghệ tạo hình gia tăng đa điểm (TPIF) ở nhiệt độ phòng, khảo sát ảnh hưởng của các thông số công nghệ đến khả năng tạo hình và chất lượng bề mặt sản phẩm đạt được.

Các mục tiêu cụ thể được luận án xác định gồm:

  1. Thiết kế và chế tạo dụng cụ tạo hình, đồ gá, gối đỡ chuyên dụng để thực hiện toàn bộ các thí nghiệm với công nghệ tạo hình TPIF.
  2. Xây dựng mô hình côn nón có đường sinh dạng bậc nhằm khảo sát các góc tường từ 65° đến 85° trong quá trình tạo hình TPIF, giúp tiết kiệm mẫu thí nghiệm và nâng cao hiệu quả khảo sát.
  3. Nghiên cứu tạo hình tấm hợp kim nhôm A1050 H14 có độ dày 1,5 mm bằng công nghệ TPIF; xác định khả năng tạo hình thông qua góc tường sản phẩm $\alpha$ đạt được; khảo sát ảnh hưởng của 4 thông số công nghệ (bước tiến dụng cụ, tốc độ quay trục chính, tốc độ tiến dụng cụ, đường kính dụng cụ) và xây dựng phương trình hồi quy thể hiện mối liên hệ giữa các thông số này với góc biến dạng.
  4. Thực hiện mô phỏng phần tử hữu hạn (FEM) bằng phần mềm Abaqus nhằm dự đoán vết nứt xuất hiện trong quá trình tạo hình TPIF với các thông số tương tự thực nghiệm.
  5. Khảo sát ảnh hưởng của chất bôi trơn và phương pháp bôi trơn thích hợp để đạt chất lượng bề mặt tốt và khả năng biến dạng cao của công nghệ TPIF.

Các câu hỏi nghiên cứu tương ứng:

  • Thiết kế kết cấu đồ gá và dưỡng bán phần nào bảo đảm chuyển động trượt của tấm kim loại theo phương trục Z đồng bộ với chuyển động của dụng cụ tạo hình?
  • Quy luật tác động độc lập và tương tác của bước tiến dụng cụ, vận tốc tiến, tốc độ quay trục chính và đường kính dụng cụ đến góc tạo hình giới hạn $\alpha_{max}$ của tấm nhôm A1050 H14 diễn ra như thế nào?
  • Mô hình phần tử hữu hạn trong Abaqus có dự đoán chính xác vị trí phá hủy, cao độ rách và biên dạng hình học của chi tiết so với thực nghiệm hay không?
  • Loại chất bôi trơn và phương pháp bôi trơn nào làm giảm ma sát tiếp xúc, hạn chế cào xước và mang lại độ nhám bề mặt thấp nhất cho sản phẩm tạo hình TPIF?

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

  • Đối tượng nghiên cứu: Khả năng tạo hình của công nghệ TPIF trên tấm hợp kim nhôm A1050 H14 dày 1,5 mm thông qua 4 thông số công nghệ: tốc độ tiến dụng cụ trong mặt phẳng ($V_{xy}$), bước tiến dụng cụ sau mỗi lớp theo phương đứng ($\Delta z$), đường kính dụng cụ ($D$), tốc độ quay trục chính ($n$); cùng các loại chất bôi trơn và phương pháp bôi trơn tác động đến chất lượng bề mặt và khả năng biến dạng.
  • Phạm vi không gian và thiết bị: Nghiên cứu thực hiện tại Phòng Thí nghiệm trọng điểm Điều khiển số và Kỹ thuật Hệ thống (DCSELAB) - Đại học Quốc gia TP.HCM, trên máy gia công ISF chuyên dụng.
  • Phạm vi vật liệu và nhiệt độ: Thực nghiệm thực hiện hoàn toàn ở nhiệt độ phòng trên hợp kim nhôm A1050 H14, chiều dày ban đầu $t_0 = 1,5$ mm.
  • Phạm vi thông số công nghệ khảo sát:
    • Vận tốc tiến dụng cụ: $V_{xy} = 300, 900, 1500$ mm/phút.
    • Bước tiến dụng cụ theo phương Z: $\Delta z = 0,1; 0,8; 1,5$ mm.
    • Đường kính dụng cụ tạo hình: $D = 6, 12, 18$ mm.
    • Tốc độ quay trục chính: $n = 300, 1050, 1800$ vòng/phút.
  • Phạm vi chất bôi trơn khảo sát: Dầu nhờn công nghiệp (Gear VG 150 EP), mỡ bôi trơn công nghiệp, bột than chì (graphite), hỗn hợp mỡ bôi trơn kết hợp than chì (tỉ lệ 1:1), dầu hướng dương.

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

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

Nguyên lý gia công kim loại tấm bằng công nghệ biến dạng gia tăng (ISF) được Edward Leszak công bố vào năm 1967. Đến năm 1978, Mason thực hiện các nghiên cứu về công nghệ tạo hình gia tăng đơn điểm (SPIF). Tổng kết lịch sử phát triển của công nghệ ISF trong thế kỷ 20 được Emmens (2010) hệ thống hóa toàn diện.

Trong hướng nghiên cứu tạo hình SPIF cho vật liệu dẻo ở nhiệt độ phòng:

  • Franzen (2008, 2009) khảo sát biến dạng tấm nhựa PVC với đường chạy dao xoắn ốc bước 0,5 mm, đường kính dao 10–15 mm, chiều dày tấm 2–3 mm sử dụng bôi trơn dạng soap-emulsion, chứng minh bước tiến $\Delta z$ và đường kính dụng cụ tác động lớn nhất đến góc tạo hình.
  • Crina Radu (2013) nghiên cứu ảnh hưởng của 4 thông số ($D = 3\div 5$ mm, $\Delta z = 0,05\div 0,5$ mm, $V_{xy} = 1500\div 3000$ mm/phút, $n = 500\div 1000$ vòng/phút) đến độ nhám bề mặt tấm thép carbon thấp dày 0,6 mm, ghi nhận $\Delta z$ làm tăng độ nhám mạnh nhất trong khi tăng $D, n, V_{xy}$ giúp giảm độ nhám.
  • Rauch (2009) nghiên cứu ảnh hưởng của đường chạy dao CAM đối với hợp kim nhôm 5086 dày 0,6 mm, chỉ ra rằng quỹ đạo CAM thương mại chưa phù hợp cho ISF và đề xuất xây dựng quỹ đạo dựa trên phản hồi lực tạo hình.

Trong hướng nghiên cứu ISF kết hợp gia nhiệt cho vật liệu khó biến dạng:

  • Ambrogio (2012) dùng hệ thống điện trở gia nhiệt không khí (200–300°C) cho hợp kim magie AZ31, nâng góc tường từ 30° ở nhiệt độ phòng lên 60° ở 250°C.
  • Lê Văn Sỹ (2014, 2017) nghiên cứu gia nhiệt bằng dòng điện một chiều theo định luật Joule-Lenz cho hợp kim nhôm A1010 và magie AZ31, khảo sát $T = 200\div 300^\circ$C, $\Delta z = 0,2\div 1$ mm, $V_{xy} = 1000\div 6000$ mm/phút.
  • Fan (2010) và David Adams (2013) áp dụng gia nhiệt điện trở khảo sát hợp kim magie AZ31 và titan TiA2Mn1,5.
  • Duflou (2007, 2008) ứng dụng gia nhiệt cục bộ bằng chùm tia laser, chứng minh khả năng cải thiện độ chính xác hình học chi tiết so với mô hình CAD.
  • Galdos (2012) dùng dầu lỏng gia nhiệt hợp kim magie AZ31 lên 250°C, kiểm soát sai số vùng kẹp trong phạm vi 2 mm.

Trong hướng nghiên cứu công nghệ tạo hình gia tăng đa điểm (TPIF):

  • Attanasio (2008) so sánh thực nghiệm tạo hình khoang cửa xe ô tô giữa SPIF và TPIF, chứng minh TPIF cải thiện 25–30% độ chính xác hình học và nâng cao chất lượng bề mặt ở các góc lượn.
  • Isabel Bagudanch (2015) khảo sát tạo hình PVC và PC hình chóp cụt sinh tròn bằng TPIF, khẳng định TPIF giảm thiểu độ đàn hồi phục hồi (springback) và tăng độ chính xác đáy chi tiết.
  • Jeswiet (2008) đo lực tạo hình trên tấm hợp kim nhôm AA3003-0 dày 1,21 mm, ghi nhận điểm cực đại của lực tại thời điểm phá hủy ở góc biến dạng giới hạn $\alpha_{max}$.
  • Chenhao Wang (2015) mô phỏng phần tử hữu hạn Abaqus/Explicit kết hợp chương trình con UMAT với mô hình phá hủy liên tục (CDM) cho hợp kim nhôm AA7075 trong TPIF.
  • Perez-Santiago (2016) phân tích lực tạo hình và độ giãn mỏng chiều dày trên mô hình phần tử hữu hạn.
  • Jaber (2019) nghiên cứu biến dạng nhiều giai đoạn trên hợp kim nhôm AA1070 dày 1 mm để đạt góc thành 90° bằng phần mềm ANSYS.
  • Haibo Lu (2016) sử dụng thuật toán điều khiển dự báo mô hình để bù sai số quỹ đạo dao TPIF với dưỡng bán phần trên hợp kim nhôm AA7075-O dày 1,6 mm.
  • Asghari (2016) tối ưu hóa bề mặt đáp ứng (RSM) cho nhôm A1050 tạo hình nón bằng TPIF ($D = 15$ mm, $\alpha = 63^\circ$, $n = 800$ vòng/phút, $\Delta z = 0,2$ mm).
  • Safari (2018) nghiên cứu tạo hình chi tiết có cả phần lồi và lõm trên hợp kim nhôm AA3105 dày 1 mm.
  • Mostafanezhad (2018) áp dụng quy hoạch Box-Behnken và phân tích ANOVA khảo sát lực tạo hình và tỷ lệ chiều dày nhôm A1050.
  • Meier (2009) dùng hai cánh tay robot công nghiệp phối hợp chuyển động biến dạng TPIF.
  • Li (2018) nghiên cứu ảnh hưởng của tấm nội suy đến lực và độ nhám bề mặt composite trong TPIF.

Trong hướng nghiên cứu bôi trơn cho ISF:

  • Diabb (2017) khảo sát dầu thực vật bổ sung hạt nano $\text{SiO}_2$ (0,0125–0,1% khối lượng) trên hợp kim nhôm 6061 trong SPIF.
  • Sornsuwit (2017) kiểm tra phương pháp phun bôi trơn trên SUS304 và mỡ $\text{MoS}_2$ trên TiGr2.
  • Jawale (2018) khảo sát dầu khoáng cho tấm đồng, kết luận bôi trơn cải thiện bề mặt nhưng ít tác động đến góc giới hạn.
  • Shisode (2018) xây dựng mô hình ma sát biên và ma sát thủy động lực dựa trên phương trình Reynolds kết hợp FEM.

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

  • Nguyễn Thanh Nam (2012) khảo sát SPIF trên hợp kim nhôm A1050 H14 với các thông số $\Delta z, V_{xy}, D, n$.
  • Phạm Văn Trung (2014) và Nguyễn Văn Năng (2015) khảo sát tạo hình tấm PVC và composite bằng SPIF.
  • Phan Đình Tuấn (2014) ứng dụng SPIF chế tạo vỏ mô hình xe hơi.
  • Lê Khánh Điền (2016) lập phương trình hồi quy 4 thông số cho nhôm A1050 H14, thép SS330 và SUS304 về góc biến dạng, độ nhám, độ phục hồi và năng suất trong SPIF.
  • Phạm Quốc Tuấn (2015) nghiên cứu mô phỏng SPIF và nhận dạng đường cong giới hạn biến dạng (FLD) khi đứt gãy.
  • Lê Văn Sỹ (2017) nghiên cứu mô phỏng dự đoán khuyết tật tạo hình tấm kim loại ở nhiệt độ cao trong SPIF.
  • Võ Tuyển (2018) và Bùi Anh Phi (2019) nghiên cứu thực nghiệm và mô phỏng lượng phục hồi đàn hồi theo phương hướng kính trên vật liệu titan bằng Hot-SPIF.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết

Các công bố trong nước trước đây hầu như chỉ tập trung vào phương pháp SPIF. Luận án của Ma Văn Việt tập trung vào công nghệ TPIF, giải quyết bốn khoảng trống kỹ thuật và học thuật:

  1. Thiết kế và chế tạo hệ thống đồ gá TPIF hoạt động theo nguyên lý khung kẹp trượt đồng bộ theo trục Z kết hợp dưỡng bán phần cố định.
  2. Thiết kế mô hình mẫu thử nghiệm dạng côn bậc (stepped cone) khảo sát dải góc biến dạng dốc từ 65° đến 85° trong một chu trình đơn, thay thế việc gia công nhiều mẫu hình nón riêng lẻ.
  3. Thiết lập mô hình toán học thực nghiệm bậc 2 (Minitab, ANOVA) xác định mức độ tác động và tương tác của 4 thông số công nghệ ($\Delta z, V_{xy}, D, n$) đến góc tường giới hạn $\alpha$ trên nhôm A1050 H14 dày 1,5 mm.
  4. Đánh giá thực nghiệm toàn diện ảnh hưởng của 5 loại chất bôi trơn và phương pháp bôi trơn đến chất lượng vi cấu trúc bề mặt (SEM) và độ nhám $R_z$ trong TPIF, kết hợp mô phỏng phần tử hữu hạn Abaqus để dự đoán vị trí nứt rách.

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

Khung lý thuyết và khái niệm sử dụng

Luận án vận dụng ba trụ cột lý thuyết chính:

  1. Lý thuyết biến dạng dẻo kim loại: Nghiên cứu trạng thái ứng suất - biến dạng tại một điểm và trạng thái ứng suất khối khi chịu tải trọng cục bộ từ đầu dụng cụ hình cầu; phân tích cơ chế giãn mỏng cục bộ theo định luật bảo toàn thể tích, xác định chiều dày chi tiết sau biến dạng $t_1 = t_0 \cdot \cos(\alpha)$ (định luật Sine) và giới hạn biến dạng trước khi nứt rách.
  2. Lý thuyết quy hoạch thực nghiệm (Design of Experiment - DOE): Sử dụng mô hình quy hoạch bậc 2 (Box-Behnken / quy hoạch xoay) 4 nhân tố trên phần mềm Minitab; phân tích phương sai (ANOVA) để kiểm định ý nghĩa thống kê của các hệ số; sử dụng biểu đồ Pareto để xếp hạng mức độ ảnh hưởng của các thông số; xây dựng bề mặt đáp ứng (Response Surface Methodology - RSM) và phương trình hồi quy đa thức biểu diễn góc biến dạng $\alpha$.
  3. Lý thuyết cơ học vật rắn biến dạng và phương pháp phần tử hữu hạn (FEM): Thiết lập phương trình cân bằng động lực học phi tuyến giải bằng thuật toán tích phân tường minh (Explicit) và tích phân ẩn (Implicit) trong Abaqus. Phân tích điều kiện ổn định của tích phân tường minh dựa trên bước thời gian tới hạn Courant-Friedrichs-Lewy ($\Delta t \le L^{e} / c_{d}$, với $c_d$ là vận tốc sóng âm trong vật liệu).

Phương pháp nghiên cứu thực nghiệm

  • Hệ thống thiết bị thí nghiệm: Máy tạo hình ISF chuyên dùng do DCSELAB chế tạo; đồ gá TPIF chuyên dụng gồm khung trượt dẫn hướng theo trục Z và dưỡng tạo hình bán phần bằng kim loại; đầu dụng cụ dạng viên bi hàn cứng trên cán trụ để giảm mài mòn, với 3 kích thước đường kính $D = 6$ mm, $12$ mm, $18$ mm.
  • Vật liệu phôi: Hợp kim nhôm A1050 H14 tấm, chiều dày danh nghĩa $t_0 = 1,5$ mm. Thành phần hóa học được kiểm định bằng phương pháp quang phổ phát xạ theo tiêu chuẩn BS EN 485-1; cơ tính (giới hạn chảy, độ bền kéo, độ giãn dài) được xác định qua thử kéo trên máy thử vạn năng theo tiêu chuẩn ISO 6361-2:2014.
  • Mô hình chi tiết khảo sát: Mô hình nón cụt sinh bậc (stepped cone) có các góc nghiêng biến thiên liên tục từ 65° đến 85°. Đường chạy dao được lập trình theo quỹ đạo xoắn ốc liên tục (Helical toolpath) trong phần mềm CAM để loại trừ vết tỳ đè khi chuyển lớp.
  • Kỹ thuật đo đạc và thiết bị phân tích:
    • Thước Panme cơ khí Mitutoyo dùng đo độ dày phôi trước và sau biến dạng.
    • Thước đo cao Mitutoyo dùng xác định cao độ vết rách ($h$), từ đó chuyển đổi hình học sang góc biến dạng giới hạn $\alpha$.
    • Máy đo tọa độ ba chiều (CMM) kiểm tra sai lệch biên dạng hình học thực tế so với mô hình CAD.
    • Máy đo độ nhám bề mặt Mitutoyo SJ 301 xác định thông số độ nhám $R_z$.
    • Kính hiển vi điện tử quét (SEM) HITACHI S-4800 chụp ảnh tế vi bề mặt để đánh giá vết xước, cào mòn và kết dính vật liệu.
Thông số công nghệ Ký hiệu Mức thấp (-1) Mức cơ sở (0) Mức cao (+1) Đơn vị
Bước tiến dụng cụ theo trục Z $\Delta z$ 0,1 0,8 1,5 mm
Tốc độ tiến dao mặt phẳng XY $V_{xy}$ 300 900 1500 mm/phút
Đường kính dụng cụ tạo hình $D$ 6 12 18 mm
Tốc độ quay trục chính $n$ 300 1050 1800 vòng/phút

Phương pháp mô phỏng số

Xây dựng mô hình hình học 3D trong môi trường Abaqus/Explicit:

  • Tấm kim loại phôi nhôm A1050 H14 được mô hình hóa là vật thể biến dạng dẻo (deformable), gán các thông số đàn hồi ($E, \nu$) và đường cong ứng suất - biến dạng dẻo thực tế từ thử nghiệm kéo.
  • Dụng cụ tạo hình và dưỡng được xem là vật thể rắn tuyệt đối (analytical/discrete rigid body).
  • Điều kiện biên mô phỏng chính xác trạng thái khung kẹp trượt tự do theo phương Z nhưng hạn chế các bậc tự do quay và tịnh tiến ngang. Quỹ đạo dịch chuyển của đầu bi được nội suy theo thời gian từ mã G-code thực tế.

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

Mở đầu

Phần Mở đầu thiết lập cơ sở khoa học của đề tài: phân tích bối cảnh công nghiệp chuyển dịch sang sản xuất linh hoạt; trình bày tính cấp thiết của việc nghiên cứu công nghệ TPIF; công bố 5 mục tiêu cụ thể, 4 nhiệm vụ nghiên cứu, đối tượng nghiên cứu (nhôm A1050 H14 dày 1,5 mm), phạm vi thông số thực nghiệm, hướng tiếp cận lý thuyết kết hợp thực nghiệm và mô phỏng số; nêu rõ ý nghĩa khoa học và giá trị thực tiễn của công trình.

Chương 1: Tổng quan

Chương 1 hệ thống hóa toàn diện các nghiên cứu về công nghệ tạo hình gia tăng kim loại tấm:

  • Phân loại công nghệ ISF dựa trên cơ chế biến dạng (tiếp xúc đơn mặt SSCIF, tiếp xúc hai mặt DSCIF, tiếp xúc kết hợp); theo chiến lược gia công (một giai đoạn, nhiều giai đoạn); theo nhiệt độ (nhiệt độ thường, nhiệt độ cao); theo thiết bị (máy CNC, robot công nghiệp, máy chuyên dùng ISF).
  • So sánh sự khác biệt bản chất giữa SPIF và TPIF:
Đặc tính so sánh Công nghệ tạo hình SPIF Công nghệ tạo hình TPIF
Kết cấu đồ gá Đơn giản, tấm kẹp cố định trên bàn máy Phức tạp, tấm kẹp trên khung trượt theo trục Z
Sử dụng dưỡng/khuôn Hoàn toàn không sử dụng khuôn/dưỡng Sử dụng dưỡng cố định toàn phần hoặc bán phần
Khả năng tạo hình Chỉ tạo được chi tiết dạng lồi hoặc lõm đơn Tạo được chi tiết phức tạp có cả vùng lồi và lõm
Độ chính xác kích thước Thấp hơn, chịu ảnh hưởng lớn của springback Cao hơn, cải thiện 25–30% so với SPIF
Chi phí và chuẩn bị Chi phí gá đặt thấp, tính linh hoạt rất cao Chi phí gá đặt cao hơn, cần thiết kế dưỡng bán phần
  • Điểm lại các công trình nghiên cứu trong nước (Nguyễn Thanh Nam, Lê Khánh Điền, Lê Văn Sỹ, Phạm Quốc Tuấn, Võ Tuyển, Bùi Anh Phi) và quốc tế (Franzen, Radu, Ambrogio, Attanasio, Bagudanch, Jeswiet, Wang, Lu, Asghari, Mostafanezhad).
  • Tổng quan về chất bôi trơn trong ISF (Diabb, Sornsuwit, Jawale, Shisode), xác lập căn cứ khoa học để lựa chọn 4 thông số công nghệ và khảo sát chế độ bôi trơn cho TPIF.

Chương 2: Cơ sở lý thuyết

Chương 2 trình bày các nền tảng toán học và cơ học làm công cụ cho luận án:

  • Lý thuyết biến dạng kim loại: Thiết lập các biểu thức tính toán ứng suất chính, ứng suất tương đương Von-Mises, biến dạng trượt và trạng thái biến dạng mỏng cục bộ; giải thích cơ chế biến dạng từng lớp dưới tiếp xúc cục bộ của đầu bi.
  • Lý thuyết quy hoạch thực nghiệm: Trình bày phương pháp quy hoạch thực nghiệm bậc hai Box-Behnken; các bước xây dựng ma trận thực nghiệm trực giao; phương pháp phân tích phương sai ANOVA; cách thức đánh giá sự phù hợp của mô hình thông qua hệ số xác định $R^2$ và giá trị p-value; thiết lập các phương trình bề mặt đáp ứng bậc hai.
  • Phương pháp phần tử hữu hạn trong ISF: Trình bày chi tiết toán học của phương pháp tích phân thời gian tường minh (Explicit) và tích phân ẩn (Implicit). Phân tích tiêu chuẩn ổn định bước thời gian giải tích cho bài toán tiếp xúc va chạm và biến dạng dẻo lớn trong Abaqus.

Chương 3: Hệ thống thí nghiệm

Chương 3 trình bày chi tiết về phần cứng thực nghiệm, thiết kế đồ gá và thiết lập thông số:

  • Thiết bị và đồ gá: Giới thiệu máy chuyên dùng ISF của DCSELAB. Bản vẽ và kết cấu đồ gá TPIF gồm khung kẹp tấm phôi có cơ cấu trượt tự do theo chiều thẳng đứng dọc trục Z và hệ thống dưỡng bán phần cố định định vị trên bàn máy.
  • Dụng cụ tạo hình: Chế tạo 3 đầu dụng cụ đường kính $D = 6, 12, 18$ mm từ bi thép chịu lực hàn cứng vào thân trụ, kiểm tra độ tròn và chất lượng bề mặt đầu bi.
  • Mô hình mẫu và quỹ đạo dao: Thiết kế mô hình côn nón sinh bậc khảo sát dải góc từ 65° đến 85°; lập trình quỹ đạo chạy dao dạng đường xoắn ốc (Helical toolpath) để duy trì sự tiếp xúc liên tục và ổn định tải trọng biến dạng.
  • Chất bôi trơn khảo sát: Lựa chọn 5 môi trường bôi trơn gồm mỡ công nghiệp, dầu Gear VG 150 EP, bột than chì, hỗn hợp mỡ và than chì (1:1), dầu hướng dương.
  • Quy trình đo lường: Xây dựng công thức chuyển đổi hình học từ cao độ vết rách ($h$) đo bằng thước đo cao Mitutoyo sang góc biến dạng giới hạn $\alpha$. Giới thiệu quy trình đo biên dạng trên máy CMM, đo độ nhám bằng Mitutoyo SJ 301 và phân tích cấu trúc tế vi bằng SEM HITACHI S-4800.

Chương 4: Mô phỏng quá trình TPIF

Chương 4 tập trung vào xây dựng và kiểm chứng mô hình mô phỏng số bằng phần mềm Abaqus:

  • Mô hình hóa hình học: Thiết lập dưỡng tạo hình và đầu dụng cụ là vật thể cứng tuyệt đối (Rigid Body); tấm kim loại nhôm A1050 H14 là vật thể biến dạng dẻo (Deformable Body).
  • Tạo lưới và điều kiện biên: Thiết lập lưới phần tử vỏ (shell elements) với mật độ lưới mịn tại vùng biến dạng cục bộ; gán điều kiện biên kẹp chặt mép tấm trên khung trượt theo trục Z; cố định toàn bộ các bậc tự do của dưỡng bán phần; định nghĩa tương tác tiếp xúc bề mặt (surface-to-surface contact) có ma sát giữa đầu dao - phôi và phôi - dưỡng.
  • Thông số vật liệu: Khai báo môđun đàn hồi $E$, hệ số Poisson $\nu$, và bảng số liệu ứng suất - biến dạng dẻo thực nghiệm của nhôm A1050 H14.
  • Kết quả mô phỏng: Thực hiện mô phỏng hai trường hợp cụ thể:
    • Phân tích sự phân bố ứng suất tương đương Von-Mises và biến dạng dẻo tích lũy (PEEQ) theo thời gian gia công.
    • Dự đoán chính xác vị trí mỏng hóa cục bộ lớn nhất và thời điểm hình thành vết rách của tấm kim loại.
    • Xác định mặt cắt hình học và cao độ vết rách, cho thấy độ trùng khớp cao giữa kết quả mô phỏng Abaqus và vết nứt thực nghiệm.

Chương 5: Kết quả nghiên cứu

Chương 5 công bố toàn bộ dữ liệu thực nghiệm, phân tích thống kê và đánh giá cơ chế:

  • Đặc tính vật liệu thực nghiệm: Bảng thành phần hóa học phân tích quang phổ của tấm nhôm A1050 H14 đạt tiêu chuẩn BS EN 485-1; kết quả thử kéo mẫu theo tiêu chuẩn ISO 6361-2:2014 xác định cơ tính ban đầu của vật liệu.
  • Kết quả ma trận thực nghiệm và phân tích ANOVA:
    • Thực hiện toàn bộ ma trận thí nghiệm theo quy hoạch bậc 2 với 4 thông số đầu vào.
    • Bảng phân tích phương sai ANOVA cho thấy mô hình hồi quy có ý nghĩa thống kê cao.
    • Biểu đồ Pareto và bảng phần trăm đóng góp xác định thứ tự ảnh hưởng đến góc biến dạng: bước tiến $\Delta z$ và đường kính dao $D$ có tác động chi phối lớn nhất, tiếp theo là vận tốc tiến $V_{xy}$, trong khi tốc độ quay trục chính $n$ có mức độ ảnh hưởng thấp nhất.
    • Thiết lập phương trình hồi quy toán học mô tả góc biến dạng $\alpha$ theo các biến $\Delta z, V_{xy}, D, n$ và các số hạng tương tác bậc hai.
    • Xác định bộ thông số công nghệ hợp lý để đạt được góc tạo hình lớn nhất.
  • Ảnh hưởng của chất bôi trơn và phương pháp bôi trơn:
    • So sánh định lượng giá trị độ nhám bề mặt ngoài $R_z$ của sản phẩm khi sử dụng các chất bôi trơn khác nhau.
    • Phân tích ảnh chụp SEM HITACHI S-4800 chỉ rõ cơ chế bôi trơn: việc sử dụng than chì hoặc hỗn hợp mỡ - than chì tạo lớp màng bôi trơn rắn ngăn cản hiện tượng dính bám kim loại giữa đầu bi và tấm nhôm, làm giảm rõ rệt các vết cào xước vi mô và hạ thấp giá trị độ nhám bề mặt so với gia công khan hoặc dùng dầu bôi trơn thông thường.
  • So sánh mô phỏng và thực nghiệm:
    • Tích hợp dữ liệu đo tọa độ CMM thực tế với biên dạng CAD và biên dạng mô phỏng Abaqus trong hai trường hợp thử nghiệm.
    • So sánh sai số cao độ vết rách và góc biến dạng giữa dự báo số và đo đạc thực tế, khẳng định tính chính xác của mô hình phần tử hữu hạn đã xây dựng.

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. Thiết lập mô hình tương quan toán học cho TPIF: Luận án đã xây dựng thành công phương trình hồi quy thực nghiệm bậc hai mô tả định lượng mối quan hệ giữa 4 thông số công nghệ ($\Delta z, V_{xy}, D, n$) với góc biến dạng giới hạn $\alpha$ của tấm hợp kim nhôm A1050 H14 dày 1,5 mm. Kết quả phân tích ANOVA và biểu đồ Pareto làm sáng tỏ quy luật tác động: bước tiến dụng cụ $\Delta z$ và đường kính dụng cụ $D$ là hai yếu tố chi phối mạnh mẽ nhất đến khả năng tạo hình.
  2. Xây dựng mô hình phần tử hữu hạn Abaqus/Explicit cho TPIF: Luận án giải quyết thành công bài toán mô phỏng số tiếp xúc phi tuyến phức tạp trong điều kiện biến dạng lớn, tiếp xúc thay đổi liên tục và khung kẹp trượt tự do theo phương Z; dự đoán tin cậy quy luật phân bố độ mỏng hóa chiều dày, vùng tập trung biến dạng và cao độ xuất hiện vết nứt rách của phôi.
  3. Làm rõ cơ chế ma sát và bôi trơn trong TPIF: Thông qua phân tích ảnh tế vi SEM, luận án chỉ ra cơ chế tổn hại bề mặt do cào xước và kết dính vật liệu trong quá trình biến dạng gia tăng; chứng minh hiệu quả của màng bôi trơn rắn (than chì, hỗn hợp mỡ - than chì) trong việc giảm ma sát tiếp xúc và cải thiện độ nhám bề mặt chi tiết.

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

  1. Chế tạo hệ thống đồ gá TPIF chuyên dụng: Thiết kế, tính toán độ cứng vững và chế tạo thành công bộ đồ gá TPIF hoạt động ổn định trên máy ISF chuyên dùng, với khung kẹp trượt đồng bộ theo trục Z và dưỡng tạo hình bán phần, phục vụ trực tiếp công tác thực nghiệm và chuyển giao công nghệ.
  2. Sáng tạo mô hình mẫu thử côn bậc: Đề xuất và áp dụng thành công mô hình nón cụt có đường sinh dạng bậc bao phủ dải góc tường $65^\circ\div 85^\circ$, cho phép xác định góc tạo hình giới hạn chỉ trong một lần chạy mẫu, tiết kiệm đáng kể phôi liệu và thời gian thí nghiệm.
  3. Xác lập bộ thông số công nghệ và chế độ bôi trơn hợp lý: Đưa ra bộ thông số tối ưu gồm bước tiến $\Delta z$, vận tốc tiến $V_{xy}$, đường kính $D$, tốc độ quay $n$ cùng tổ hợp chất bôi trơn phù hợp, bảo đảm đạt góc tạo hình lớn nhất và bề mặt chi tiết có độ nhám thấp, cung cấp hướng dẫn kỹ thuật cho các ứng dụng sản xuất thực tế.

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

  • Đề xuất ứng dụng quy trình công nghệ TPIF và bộ đồ gá khung trượt vào sản xuất các sản phẩm nhôm tấm đơn chiếc, loạt nhỏ hoặc tạo mẫu nhanh trong các xưởng gia công cơ khí chính xác tại Việt Nam.
  • Khuyến nghị sử dụng chất bôi trơn than chì hoặc hỗn hợp mỡ - than chì tỉ lệ 1:1 trong tạo hình nhôm tấm bằng TPIF để bảo vệ bề mặt sản phẩm và kéo dài tuổi thọ đầu dụng cụ.
  • Đề xuất tích hợp mô hình mô phỏng số Abaqus vào giai đoạn thiết kế quy trình CAM nhằm dự đoán trước khả năng nứt rách và bù trừ sai số hình học trước khi gia công thực tế.

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

Hạn chế của đề tài

  1. Phạm vi vật liệu và chiều dày: Toàn bộ khảo sát thực nghiệm và mô phỏng chỉ mới thực hiện trên một mác vật liệu duy nhất là hợp kim nhôm A1050 H14 với một cấp chiều dày ban đầu cố định ($t_0 = 1,5$ mm).
  2. Điều kiện nhiệt độ gia công: Nghiên cứu chỉ giới hạn ở nhiệt độ phòng, chưa mở rộng sang tạo hình gia nhiệt (Hot-TPIF) cho các vật liệu khó biến dạng (titan, hợp kim magie, thép không gỉ).
  3. Thời gian tính toán mô phỏng: Do đặc thù của bài toán biến dạng lớn và tiếp xúc thay đổi liên tục theo quỹ đạo xoắn ốc dài, thời gian tính toán của mô hình phần tử hữu hạn Abaqus/Explicit còn lớn.
  4. Vấn đề mòn dụng cụ: Nghiên cứu chưa khảo sát định lượng tốc độ mài mòn của đầu bi dụng cụ trong điều kiện làm việc liên tục dài hạn.

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

  • Mở rộng nghiên cứu công nghệ TPIF trên các vật liệu khác như hợp kim nhôm series 5xxx, 6xxx, hợp kim titan và thép không gỉ.
  • Thiết kế tích hợp hệ thống gia nhiệt (gia nhiệt bằng dòng điện hoặc cảm ứng) vào đồ gá TPIF để gia công các tấm vật liệu khó biến dạng ở nhiệt độ cao.
  • Nghiên cứu phát triển thuật toán tự động bù sai số hình học và độ phục hồi đàn hồi (springback) trực tiếp vào quỹ đạo chạy dao CAM.
  • Áp dụng các kỹ thuật chia lưới thích ứng (adaptive meshing) và tăng tốc khối lượng (mass scaling) hợp lý để tối ưu hóa thời gian tính toán mô phỏng phần tử hữu hạn.

Giá trị tham khảo

Đối tượng thụ hưởng

  • 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í, Cơ kỹ thuật, Gia công áp lực và Công nghệ Chế tạo máy có tài liệu tham khảo chi tiết về lý thuyết biến dạng gia tăng, quy hoạch thực nghiệm và mô phỏng số Abaqus.
  • Kỹ sư công nghệ và thiết kế: Tại các viện nghiên cứu, doanh nghiệp sản xuất ô tô, hàng không vũ trụ, thiết bị y tế và sản phẩm gia dụng cần giải pháp tạo mẫu nhanh hoặc sản xuất loạt nhỏ từ kim loại tấm.
  • Nhà hoạch định công nghệ sản xuất: Cân nhắc ứng dụng công nghệ ISF/TPIF nhằm giảm chi phí đầu tư khuôn mẫu đắt tiền trong các dự án phát triển sản phẩm mới.

Các nội dung đáng tham khảo

  • Hồ sơ thiết kế đồ gá TPIF: Kết cấu cơ khí của khung kẹp trượt dẫn hướng trục Z và dưỡng tạo hình bán phần (Chương 3).
  • Phương pháp mẫu côn bậc: Giải pháp kỹ thuật sử dụng mô hình côn bậc để rút ngắn số lần thử nghiệm khi xác định giới hạn tạo hình $\alpha_{max}$ (Chương 3).
  • Quy trình thiết lập mô hình FEM trong Abaqus/Explicit: Cách khai báo tiếp xúc phi tuyến, điều kiện biên trượt và dự đoán rách tấm kim loại (Chương 4).
  • Dữ liệu thực nghiệm và phương trình hồi quy: Toàn bộ bảng số liệu ma trận thực nghiệm, kết quả phân tích ANOVA và bảng so sánh các chất bôi trơn (Chương 5).

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

1. Tại sao luận án lại sử dụng dưỡng tạo hình cố định bán phần thay vì dưỡng toàn phần trong công nghệ TPIF?

Luận án lựa chọn dưỡng tạo hình cố định bán phần vì phương án này mang lại ưu thế vượt trội về tính linh hoạt và chi phí. Sử dụng dưỡng bán phần giúp người chế tạo không phải gia công toàn bộ biên dạng phức tạp của chi tiết lên dưỡng khuôn như phương pháp toàn phần, từ đó giảm đáng kể thời gian thiết kế, chi phí vật liệu và chế tạo khuôn, trong khi vẫn bảo đảm đủ độ cứng vững để định hình các vùng đáy và hạn chế biến dạng uốn không mong muốn.

2. Mô hình nón cụt sinh bậc mang lại lợi ích gì cho quá trình thực nghiệm so với các mô hình nón truyền thống?

Mô hình nón cụt có đường sinh dạng bậc được thiết kế với các góc tường tăng dần từ 65° đến 85°. Thay vì phải chế tạo và chạy thử hàng loạt mẫu nón cụt có góc tường cố định riêng biệt (mỗi mẫu một góc), mô hình sinh bậc cho phép kiểm tra khả năng biến dạng của vật liệu qua nhiều dải góc dốc khác nhau chỉ trong một lần chạy máy duy nhất. Khi tấm kim loại bị nứt rách ở một độ sâu nhất định, chiều cao vết rách đo được sẽ trực tiếp quy đổi ra góc biến dạng giới hạn $\alpha_{max}$, giúp tiết kiệm tối đa phôi thí nghiệm và thời gian gia công.

3. Trong 4 thông số công nghệ khảo sát, thông số nào tác động mạnh nhất và thông số nào tác động yếu nhất đến góc tạo hình giới hạn?

Kết quả phân tích phương sai ANOVA và biểu đồ Pareto ở Chương 5 chỉ ra rằng bước tiến dụng cụ sau mỗi lớp ($\Delta z$) và đường kính dụng cụ tạo hình ($D$) là hai thông số có mức độ ảnh hưởng lớn nhất và chi phối đến khả năng tạo hình (góc biến dạng $\alpha$). Tốc độ tiến dụng cụ trong mặt phẳng ($V_{xy}$) có ảnh hưởng ở mức độ trung bình, trong khi tốc độ quay trục chính ($n$) có mức độ ảnh hưởng nhỏ nhất đến góc biến dạng giới hạn của tấm hợp kim nhôm A1050 H14.

4. Loại chất bôi trơn nào cho kết quả tốt nhất về chất lượng bề mặt chi tiết và cơ chế tác dụng là gì?

Thực nghiệm và kết quả phân tích ảnh chụp kính hiển vi điện tử quét SEM HITACHI S-4800 chứng minh rằng bột than chì và hỗn hợp mỡ bôi trơn kết hợp than chì (tỉ lệ 1:1) mang lại chất lượng bề mặt tốt nhất với độ nhám bề mặt ngoài $R_z$ thấp nhất. Cơ chế là do các hạt than chì đóng vai trò như một lớp màng bôi trơn rắn chịu áp lực cao tại điểm tiếp xúc hẹp giữa đầu bi và tấm nhôm, ngăn chặn hiệu quả hiện tượng cào xước vi mô và bám dính kim loại vốn rất dễ xảy ra trên nhôm tấm dẻo.

5. Mô hình phần tử hữu hạn Abaqus trong luận án đã kiểm chứng những đại lượng nào so với thực nghiệm?

Mô hình Abaqus/Explicit đã thực hiện mô phỏng hai trường hợp thí nghiệm cụ thể và kiểm chứng đối chiếu với thực tế ở các đại lượng:

  • Biên dạng hình học tổng thể của sản phẩm (so sánh biên dạng CAD, mô phỏng và số liệu đo CMM).
  • Vị trí tập trung biến dạng lớn nhất dẫn đến phá hủy (vị trí rách trên thành chi tiết).
  • Chiều cao/cao độ vết rách ($h$) và góc biến dạng tại thời điểm xuất hiện vết rách. Kết quả cho thấy mô hình mô phỏng số có độ tin cậy cao và phản ánh sát với diễn biến thực nghiệm.

Kết luận

Luận án tiến sĩ của tác giả Ma Văn Việt đã giải quyết một cách có hệ thống bài toán tạo hình kim loại tấm bằng công nghệ biến dạng gia tăng đa điểm (TPIF) trên vật liệu hợp kim nhôm A1050 H14 dày 1,5 mm ở nhiệt độ phòng. Tác giả đã thiết kế, chế tạo thành công bộ đồ gá TPIF chuyên dụng có khung kẹp trượt đồng bộ theo trục Z kết hợp dưỡng bán phần, đồng thời đề xuất giải pháp mẫu côn bậc giúp xác định nhanh góc biến dạng giới hạn từ 65° đến 85°. Bằng sự kết hợp chặt chẽ giữa quy hoạch thực nghiệm bậc hai, phân tích thống kê ANOVA và mô phỏng số trên phần mềm Abaqus, công trình đã xác lập phương trình hồi quy mô tả ảnh hưởng của 4 thông số công nghệ ($\Delta z, V_{xy}, D, n$) đến khả năng tạo hình, dự đoán chính xác vị trí phá hủy và làm sáng tỏ cơ chế bôi trơn bề mặt. Toàn bộ các kết quả đạt được tạo tiền đề kỹ thuật vững chắc cho việc ứng dụng công nghệ TPIF vào sản xuất loạt nhỏ, đơn chiếc và tạo mẫu nhanh trong thực tiễn công nghiệp cơ khí.