Tổng quan luận án

Gia công tinh và xử lý bề mặt chi tiết cơ khí bằng môi trường hạt mài chuyển động thể rời trong các thiết bị gia công ly tâm - roto (ЦРО - Центробежно-роторная обработка) là một phương pháp gia công cơ khí nhằm cải thiện chất lượng bề mặt, loại bỏ ba-via và nâng cao độ bền mỏi của chi tiết. Quá trình gia công này diễn ra bên trong buồng công tác kín, nơi hỗn hợp gồm hạt mài, chi tiết gia công và dung dịch công tác chịu tác động của lực ly tâm, lực ma sát và trường vận tốc xoáy phức tạp. Khoảng trống nghiên cứu được luận án giải quyết là việc thiếu hụt các mô hình toán học - thủy động lực học tích hợp mô phỏng trường vận tốc không gian hai chiều (2D Swirl Flow) với mô hình thực nghiệm đa biến phản ánh sự phụ thuộc của cường độ mài mòn vào tập hợp đầy đủ các thông số công nghệ, động học và cơ tính của vật liệu.

Mục tiêu nghiên cứu của luận án được xác định gồm các nội dung sau:

  1. Xây dựng và thiết lập hệ phương trình thủy động lực học mô tả dòng chảy xoáy (Swirl Flow) của môi trường hạt mài trong buồng làm việc của thiết bị ly tâm - roto.
  2. Xác định các thành phần vận tốc dòng chảy (vận tốc xuyên tâm $u$, vận tốc vòng $v$, vận tốc dọc trục $w$) và áp suất $p$ trong môi trường làm việc thông qua công cụ mô phỏng phần tử hữu hạn.
  3. Thiết lập mô hình thực nghiệm tính toán cường độ mài mòn $j_a$ (%/h) của bề mặt chi tiết gia công theo các thông số chế độ công nghệ và đặc tính vật liệu.
  4. Phát triển quy trình xử lý dữ liệu thực nghiệm bằng phương pháp trích xuất ảnh số và tích phân số nhằm xác định thể tích $V_i$ và phân bố năng lượng $E_i$ tại các vùng không gian trong buồng chứa.

Đối tượng nghiên cứu của đề tài là quá trình gia công chi tiết cơ khí trong môi trường hạt mài chuyển động xoáy; trường thủy động lực học của hỗn hợp công tác bên trong buồng hình trụ của thiết bị gia công ly tâm - roto; cơ chế mài mòn và tương tác cơ học giữa hạt mài và bề mặt chi tiết.

Phạm vi nghiên cứu bao gồm thiết bị gia công ly tâm - roto với buồng chứa hình trụ; mô hình dòng xoáy 2D trong phần mềm Comsol Multiphysics; vật liệu hạt mài với kích thước xác định $d$; chi tiết cơ khí có độ cứng $HB$ và khối lượng $G$; các điều kiện biên ma sát trượt và thông số tỷ lệ điền đầy trong không gian làm việc.

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

Tổng quan các công trình nghiên cứu về lĩnh vực gia công trong môi trường thể rời cho thấy các phương pháp gia công thể tích truyền thống (như gia công rung động, gia công trong thùng quay, gia công dòng hạt mài áp lực) đã được nghiên cứu từ nhiều khía cạnh động học và thực nghiệm. Tuy nhiên, đối với thiết bị gia công ly tâm - roto (ЦРО), chuyển động của môi trường hạt mài có tính chất bất đối xứng và tạo thành dòng xoáy phức tạp với gradient vận tốc lớn giữa lớp đáy roto và thành buồng đứng yên.

Các nghiên cứu trước đây chủ yếu dừng lại ở việc khảo sát thực nghiệm đơn biến hoặc xây dựng mô hình giả định dòng chảy một chiều đơn giản hóa, chưa phản ánh đầy đủ bản chất không gian của các thành phần vận tốc trong môi trường nhiều pha (pha rắn của hạt mài, pha lỏng của dung dịch công tác và pha chi tiết gia công). Luận án định vị trọng tâm giải quyết bài toán mô hình hóa trường dòng chảy xoáy thông qua module chuyên dụng Swirl Flow trên nền tảng Comsol Multiphysics, kết hợp đồng thời phương pháp thực nghiệm quang học, xử lý ma trận ảnh màu và thuật toán hồi quy đa thức trong môi trường MathCAD để xây dựng mối quan hệ hàm số giữa cường độ mài mòn và các thông số công nghệ.

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

Luận án sử dụng khung phân tích kết hợp giữa cơ học chất lỏng chuyển động xoáy, lý thuyết truyền nhiệt - động lực học hạt mài, mô phỏng số bằng phương pháp phần tử hữu hạn và phương pháp quy hoạch thực nghiệm.

1. Cơ sở lý thuyết thủy động lực học dòng xoáy (Swirl Flow)

Chuyển động của môi trường hạt mài trong buồng làm việc của thiết bị ly tâm - roto được mô tả trong hệ tọa độ trụ $(r, \theta, z)$. Trường vận tốc của dòng môi trường bao gồm ba thành phần không gian:

  • $u$: Thành phần vận tốc xuyên tâm (radial velocity component, $m/s$).
  • $v$: Thành phần vận tốc vòng hay vận tốc chu vi (circumferential/azimuthal velocity component, $m/s$).
  • $w$: Thành phần vận tốc dọc trục (axial velocity component, $m/s$).

Hệ phương trình vi phân mô tả chuyển động của môi trường dòng xoáy được thiết lập với các lực thể tích: thành phần lực xuyên tâm $F_r = 0$, thành phần lực vòng $F_\theta = 0$, và thành phần lực dọc trục $F_z = g(\rho - \rho_0)$, trong đó $g$ là gia tốc trọng trường, $\rho$ là khối lượng riêng của môi trường và $\rho_0$ là khối lượng riêng chuẩn. Phương trình vi phân theo phương dọc trục $z$ có dạng:

$$\rho \left( u \frac{\partial w}{\partial r} + w \frac{\partial w}{\partial z} \right) = -\frac{\partial p}{\partial z} + \mu \left( \frac{\partial^2 w}{\partial r^2} + \frac{1}{r} \frac{\partial w}{\partial r} + \frac{\partial^2 w}{\partial z^2} \right) + F_z$$

Hệ thống cho phép module Swirl Flow trong Comsol Multiphysics xác định vận tốc góc và phân bố vận tốc ba chiều $(u, v, w)$ trên mặt cắt không gian hai chiều $(r, z)$ có tính đối xứng trục.

2. Điều kiện biên trong mô phỏng số

Điều kiện biên áp dụng cho bề mặt hình trụ (ký hiệu là mặt 2) giới hạn thể tích môi trường làm việc được xác định như sau:

  • Chuyển vị theo phương bán kính bằng 0: $R_r = 0$.
  • Các phản lực tiếp tuyến $T_x$ và $T_y$ được xác định tự động trong quá trình mô phỏng lặp.
  • Phản lực theo phương thẳng đứng $T_z$ tỷ lệ thuận với phản lực pháp tuyến theo phương hướng tâm $T_r$ thông qua định luật ma sát Coulomb: $$T_z = k_{fr} \cdot T_r$$ trong đó hệ số ma sát giữa môi trường hạt mài và thành buồng công tác được xác định là $k_{fr} = 0,3$.

3. Phương pháp xử lý dữ liệu ảnh và giải tích số

Để xác nhận và hiệu chuẩn mô hình mô phỏng với dòng chuyển động thực tế của hạt mài, luận án áp dụng quy trình xử lý ảnh và giải tích số:

  • Thu thập dữ liệu ảnh: Quá trình chuyển động của các điểm đánh dấu (dấu chữ thập màu) trên tiết diện môi trường được ghi nhận và lưu dưới định dạng bitmap (.bmp).
  • Trích xuất ma trận màu: Ma trận tương ứng với kênh màu đỏ được trích xuất và xử lý trong môi trường tính toán MathCAD nhằm tìm tọa độ điểm giao nhau của các dấu chữ thập.
  • Biến đổi tọa độ: Chuyển đổi các tọa độ ảnh điểm thành hệ tọa độ không gian thực $(r, z)$ trên tiết diện của buồng chứa thiết bị ly tâm - roto.
  • Xấp xỉ hàm số: Xây dựng các đường mức và hàm xấp xỉ trường vận tốc thông qua cặp toán tử hồi quy và nội suy regress và interp trong MathCAD: $$vs := \text{regress}(rlevel, zlevel, npoly)$$
  • Tích phân số thể tích và năng lượng: Thể tích $V_i$ của vùng không gian thứ $i$ trong buồng công tác được tính toán bằng phương pháp tích phân số theo công thức: $$V_i = \int \int \int 2\pi r , dz dr d\theta = \int_{r_i}^{r_{i+1}} \int_{z_{min}}^{level(i)} 2\pi r , dz dr$$ Năng lượng $E_i$ tích lũy trong từng phân vùng $i$ được xác định thông qua tích phân số: $$E_i = \int_{r_i}^{r_{i+1}} \int_{z_{min}}^{level(i)} 2\pi r \cdot level(i) , dz dr$$

4. Mô hình thực nghiệm xác định cường độ mài mòn

Cường độ mài mòn $j_a$ (%/h), biểu thị phần trăm khối lượng chi tiết bị tách bóc trên một đơn vị thời gian, được mô tả dưới dạng phương trình lũy thừa đa biến:

$$j_a = C \cdot \omega^{b_1} \cdot p^{b_2} \cdot d^{b_3} \cdot K_v^{b_4} \cdot K_w^{b_5} \cdot G^{b_6} \cdot HB^{b_7}$$

Các thông số trong mô hình được tổng hợp tại bảng dưới đây:

Ký hiệu Tên thông số / Biến số Đơn vị đo Ý nghĩa cơ học / Công nghệ
$j_a$ Cường độ mài mòn bề mặt %/h Mức độ suy giảm khối lượng chi tiết theo thời gian
$\omega$ Vận tốc góc / Tốc độ quay roto $s^{-1}$ Cung cấp động năng chính cho dòng hạt mài
$p$ Mức độ điền đầy hạt mài Không thứ nguyên Tỷ lệ thể tích môi trường hạt mài trong buồng
$d$ Kích thước hạt mài mm Độ hạt của môi trường công tác
$K_v$ Tỷ lệ thể tích chi tiết / chất độn Không thứ nguyên Tương quan thể tích giữa chi tiết và hạt mài
$K_w$ Tỷ lệ thể tích pha lỏng / nạp liệu Không thứ nguyên Lượng dung dịch công tác cung cấp vào hệ
$G$ Khối lượng chi tiết gia công g (hoặc kg) Tải trọng tĩnh của phôi trong quá trình đảo trộn
$HB$ Độ cứng Brinell của chi tiết Không thứ nguyên Khả năng chống lại biến dạng dẻo bề mặt
$C$ Hệ số thực nghiệm Không thứ nguyên Hằng số tỷ lệ phụ thuộc cặp vật liệu
$b_1 \dots b_7$ Các số mũ thực nghiệm Không thứ nguyên Mức độ ảnh hưởng của từng thông số đến $j_a$

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

Chương 1: Tổng quan về gia công chi tiết cơ khí trong môi trường hạt mài thể rời và thiết bị ly tâm - roto

Chương này trình bày cơ sở lý luận về các phương pháp gia công tinh bề mặt bằng hạt mài tự do. Tác giả phân tích kết cấu buồng công tác của thiết bị ly tâm - roto (ЦРО), bao gồm roto quay ở đáy và thân vỏ hình trụ đứng yên. Phân tích trạng thái chuyển động của dòng hạt mài dưới tác dụng của lực ly tâm làm môi trường dâng lên dọc theo thành buồng và cuộn ngược trở lại tâm roto, tạo thành dòng xáo trộn dạng xuyến khép kín. Chương này nêu bật sự cần thiết của việc xây dựng mô hình toán mô tả định lượng trường vận tốc và cường độ tác động cơ học của dòng môi trường lên chi tiết gia công.

Chương 2: Mô hình hóa toán học và mô phỏng số dòng chảy xoáy của môi trường hạt mài

Chương 2 tập trung thiết lập hệ phương trình vi phân thủy động lực học mô tả dòng chảy xoáy (Swirl Flow) bên trong buồng hình trụ. Tác giả sử dụng phần mềm Comsol Multiphysics để giải số trường vận tốc 3 thành phần $(u, v, w)$ và áp suất $p$ trên mặt cắt 2D đối xứng trục.

Các kết quả mô phỏng số đạt được trong chương:

  • Thiết lập điều kiện biên hình học tại thành trụ giới hạn với chuyển vị hướng kính $R_r = 0$.
  • Ứng dụng điều kiện biên ma sát trượt với hệ số ma sát $k_{fr} = 0,3$, xác định phản lực thẳng đứng $T_z = 0,3 \cdot T_r$.
  • Tính toán phân bố các thành phần vận tốc: vận tốc tiếp tuyến $v$ đạt giá trị cực đại tại vùng tiếp giáp mép ngoài của roto quay với tốc độ $\omega$ ($s^{-1}$), trong khi vận tốc dọc trục $w$ và vận tốc hướng tâm $u$ quyết định chu kỳ hoàn lưu của hạt mài trong buồng chứa.

Chương 3: Phương pháp thực nghiệm quang học, xử lý ma trận ảnh và tích phân số

Chương 3 trình bày kỹ thuật đo đạc thực nghiệm và thuật toán xử lý dữ liệu dòng hạt mài:

  • Sử dụng phương pháp chụp ảnh quang học ghi lại chuyển động của các điểm dấu chữ thập trên tiết diện công tác, xuất dữ liệu sang tập tin ảnh định dạng .bmp.
  • Trích xuất ma trận giá trị màu đỏ trong phần mềm MathCAD để tính toán tọa độ giao điểm và ánh xạ về hệ tọa độ không gian $(r, z)$ của buồng chứa.
  • Sử dụng toán tử regress và interp trong biểu thức vs := regress(rlevel, zlevel, npoly) để xây dựng bản đồ đường đồng mức mật độ dòng chảy (contour map of flow density) trên mặt cắt bán kính của thiết bị ЦРО.
  • Ứng dụng các công thức tích phân số để tính toán thể tích vùng $V_i = \iint 2\pi r , dz dr$ và năng lượng $E_i = \iint 2\pi r \cdot level(i) , dz dr$, làm cơ sở định lượng sự phân bố mật độ công suất gia công theo các phân vùng không gian bên trong buồng chứa.

Chương 4: Xây dựng hàm thực nghiệm tính toán cường độ mài mòn và tối ưu hóa chế độ gia công

Chương 4 tổng hợp các kết quả thực nghiệm để xác định các hệ số $C$ và số mũ $b_1, b_2, b_3, b_4, b_5, b_6, b_7$ trong phương trình tính toán cường độ mài mòn $j_a$. Phân tích mức độ ảnh hưởng của từng thông số công nghệ:

  • Vận tốc góc roto $\omega$ và mức độ điền đầy $p$ là các yếu tố ảnh hưởng trực tiếp đến động năng va đập và tần số tiếp xúc giữa hạt mài và bề mặt chi tiết.
  • Kích thước hạt $d$ và tỷ lệ thể tích $K_v, K_w$ ảnh hưởng đến mật độ va đập và tính linh động của môi trường hạt mài thể rời.
  • Khối lượng $G$ và độ cứng $HB$ của chi tiết quyết định cường độ biến dạng dẻo và tốc độ bóc tách kim loại trên bề mặt gia công.

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

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

  • Xây dựng mô hình toán - cơ học dòng chảy xoáy (Swirl Flow) 2D có tính đối xứng trục trong Comsol Multiphysics, mô tả chính xác trường phân bố ba thành phần vận tốc $(u, v, w)$ và áp suất $p$ của môi trường hạt mài trong thiết bị ly tâm - roto.
  • Thiết lập phương pháp tích phân số xác định thể tích $V_i$ và phân bố năng lượng $E_i$ trên tiết diện làm việc, kết hợp toán tử hồi quy đa thức regress và interp trên nền tảng MathCAD.

2. Đóng góp mới về mặt thực tiễn và công nghệ

  • Xây dựng công thức giải tích thực nghiệm đa thông số xác định cường độ mài mòn $j_a$ (%/h) của bề mặt chi tiết cơ khí, tích hợp đầy đủ các biến số công nghệ $(\omega, p, d, K_v, K_w)$ và đặc tính cơ học của phôi $(G, HB)$.
  • Phát triển giải pháp xử lý ảnh ma trận màu từ ảnh bitmap (.bmp) để số hóa trường vận tốc thực nghiệm của môi trường hạt mài.
  • Các kết quả của luận án được công bố trong 23 công trình khoa học, bao gồm 4 bài báo thuộc cơ sở dữ liệu quốc tế Web of Science / Scopus và 5 bài báo thuộc danh mục tạp chí khoa học chuyên ngành được công nhận bởi Hội đồng Chứng chỉ Nhà nước Liên bang Nga (ВАК РФ).

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

  • Hạn chế: Mô hình thủy động lực học Swirl Flow được giải trên mặt cắt 2D đối xứng trục, chưa xét đến các hiệu ứng bất đối xứng cục bộ khi chi tiết gia công có kích thước hình học lớn hoặc biên dạng bất quy tắc; hệ số ma sát tiếp xúc giữa môi trường và thành buồng được cố định ở mức $k_{fr} = 0,3$; chưa xét đến sự biến đổi độ nhớt dung dịch công tác theo nhiệt độ phát sinh trong quá trình gia công dài hạn.
  • Hướng nghiên cứu tiếp: Mở rộng mô hình mô phỏng 3D phi đối xứng; khảo sát động học dòng hạt mài với các loại chi tiết có hình học phức tạp (bánh răng, cánh tuabin); nghiên cứu ảnh hưởng của áp suất thủy tĩnh và thành phần hóa học của dung dịch công tác đến hệ số bóc tách vật liệu.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu cho các đối tượng sau:

  • Nghiên cứu sinh và học viên cao học ngành Kỹ thuật Cơ khí, Chế tạo máy: Tham khảo phương pháp thiết lập điều kiện biên và ứng dụng module Swirl Flow trong Comsol Multiphysics để giải các bài toán chuyển động môi trường nhiều pha dạng xoáy.
  • Kỹ sư công nghệ và cán bộ nghiên cứu tại các cơ sở chế tạo máy: Ứng dụng công thức tính toán $j_a$ để lựa chọn chế độ công nghệ tối ưu (tốc độ quay $\omega$, độ điền đầy $p$, tỷ lệ nạp liệu $K_v, K_w$, kích thước hạt $d$) phù hợp với từng loại vật liệu có độ cứng $HB$ và khối lượng $G$ xác định.
  • Chuyên gia tính toán giải tích số: Ứng dụng quy trình xử lý ma trận ảnh thực nghiệm .bmp và kỹ thuật tích phân số thể tích - năng lượng trong MathCAD.

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

1. Cường độ mài mòn $j_a$ trong luận án được xác định phụ thuộc vào những thông số nào?

Cường độ mài mòn $j_a$ (%/h) được xác định thông qua hàm số thực nghiệm phụ thuộc vào 7 biến số chính: vận tốc góc roto $\omega$ ($s^{-1}$), mức độ điền đầy hạt mài trong buồng làm việc $p$, kích thước hạt mài $d$ (mm), tỷ lệ thể tích chi tiết/chất độn $K_v$, tỷ lệ thể tích pha lỏng/nạp liệu $K_w$, khối lượng chi tiết $G$, độ cứng vật liệu chi tiết $HB$, cùng hệ số tỷ lệ $C$ và các số mũ thực nghiệm $b_1 \dots b_7$.

2. Điều kiện biên tại bề mặt hình trụ giới hạn buồng làm việc được thiết lập như thế nào trong mô phỏng?

Tại bề mặt hình trụ (ký hiệu mặt 2), chuyển vị theo phương bán kính được đặt bằng 0 ($R_r = 0$). Các phản lực tiếp tuyến $T_x$ và $T_y$ được xác định qua quá trình mô phỏng, trong khi phản lực theo phương thẳng đứng $T_z$ được tính theo định luật ma sát tỷ lệ thuận với phản lực pháp tuyến hướng tâm: $T_z = k_{fr} \cdot T_r$ với hệ số ma sát $k_{fr} = 0,3$.

3. Phần mềm Comsol Multiphysics đóng vai trò gì trong phương pháp nghiên cứu của luận án?

Comsol Multiphysics sử dụng module Swirl Flow (Dòng chảy xoáy) để giải hệ phương trình thủy động lực học trên mặt cắt 2D, cho phép xác định vận tốc góc, áp suất $p$ và cả ba thành phần vận tốc không gian: vận tốc xuyên tâm $u$, vận tốc vòng $v$ và vận tốc dọc trục $w$.

4. Dữ liệu thực nghiệm quang học được xử lý như thế nào trong MathCAD?

Dữ liệu chuyển động được lưu dưới dạng ảnh bitmap (.bmp). Ma trận kênh màu đỏ được trích xuất trong MathCAD để tìm tọa độ giao điểm của các dấu chữ thập và chuyển đổi sang hệ tọa độ không gian $(r, z)$. Sau đó, các toán tử regress và interp (vs := regress(rlevel, zlevel, npoly)) được dùng để xây dựng đường đồng mức mật độ dòng chảy và thực hiện tích phân số xác định thể tích $V_i$ cùng phân bố năng lượng $E_i$.

5. Kết quả nghiên cứu của luận án đã được công bố với khối lượng học thuật như thế nào?

Tác giả đã công bố 23 công trình khoa học từ kết quả của luận án, trong đó có 4 bài báo thuộc hệ thống Web of Science và Scopus, cùng 5 bài báo trên các tạp chí khoa học thuộc danh mục VAK của Liên bang Nga (ВАК РФ).

Kết luận

Luận án tiến sĩ đã giải quyết bài toán nâng cao hiệu quả gia công chi tiết cơ khí trong môi trường hạt mài chuyển động xoáy trên thiết bị ly tâm - roto thông qua sự kết hợp giữa mô phỏng số thủy động lực học và nghiên cứu thực nghiệm. Tác giả đã thiết lập thành công mô hình 2D Swirl Flow trong Comsol Multiphysics, phát triển thuật toán xử lý ma trận ảnh và tích phân số trong MathCAD để xác định trường vận tốc, thể tích $V_i$ và phân bố năng lượng $E_i$. Đồng thời, luận án đã xây dựng mô hình thực nghiệm tính toán cường độ mài mòn $j_a$ theo tập hợp các thông số công nghệ và cơ tính của vật liệu, cung cấp cơ sở khoa học tin cậy cho việc tính toán thiết kế và tối ưu hóa quá trình gia công tinh bề mặt trong ngành chế tạo máy.