Giới thiệu dự án

Trong kỷ nguyên công nghiệp chế tạo chính xác phục vụ hàng không vũ trụ, y sinh và thiết bị bán dẫn, quy trình khoan lỗ (drilling) đóng vai trò sống còn nhưng chiếm tới 30% tổng chi phí gia công lỗ cơ khí do hiện tượng tạo bavia (burr formation). Thống kê công nghiệp từ CIRP (International Academy for Production Engineering) chỉ ra rằng các nguyên công tẩy bavia thủ công hoặc thứ cấp (deburring) tiêu tốn 15–25% thời gian gia công tổng thể và là nguồn gốc của 40% phế phẩm bề mặt ở các chi tiết chịu tải động cao.

Đề tài tốt nghiệp chuyên ngành Công nghệ Chế tạo Máy tại Trường Đại học Sư phạm Kỹ thuật TP.HCM (HCMUTE) với tên: "Nghiên cứu nguyên lý hình thành bavia trong gia công kết hợp và đề xuất phương án kiểm soát sử dụng mô phỏng số" do sinh viên Võ Duy Trung thực hiện dưới sự hướng dẫn của PGS. TS. Phạm Huy Tuân đã tập trung giải quyết triệt để vấn đề này.

Mục tiêu dự án

  1. Thiết lập mô hình động học gia công khoan kết hợp rung động dọc trục (1D-VAD) nhằm xác định sự gián đoạn tiếp xúc và độ dày phoi tức thời cắt gọt $h(t)$.
  2. Xây dựng mô hình phần tử hữu hạn 3D (3D-FEM) phi tuyến động lực học trên Abaqus/Explicit mô phỏng chính xác cơ chế chảy dẻo và phá hủy vật liệu Al6061-T6 theo mô hình Johnson-Cook.
  3. Hiệu chuẩn (calibration) và xác thực (validation) kết quả mô phỏng lực cắt dọc trục ($F_t$) và chiều cao bavia ($R_f$) với sai số dưới 10% so với dữ liệu thực nghiệm.
  4. Ứng dụng phương pháp Taguchi $L_9(3^4)$ để tối ưu hóa tổ hợp 4 thông số công nghệ: Tốc độ trục chính ($S$), Lượng chạy dao ($F$), Tần số rung ($f$), Biên độ rung ($A$) nhằm tối thiểu hóa chiều cao bavia thoát (exit burr) và bavia vào (entrance burr).

Phạm vi và giới hạn

  • Vật liệu khảo sát: Hợp kim nhôm biến dạng Al6061-T6 (kích thước mẫu mô phỏng $\Phi 7 \times 0.1\text{ mm}$).
  • Dụng cụ cắt: Mũi khoan ruột khoan xoắn tiêu chuẩn thép gió/cacbua $\Phi 5\text{ mm}$, góc đỉnh $140^\circ$.
  • Giới hạn mô hình: Bỏ qua trường nhiệt liên kết (thermal-mechanical coupling) để tối ưu hóa thời gian giải Explicit; tập trung vào mô hình ứng suất - biến dạng dẻo và phá hủy giòn/dẻo cơ học.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Trong khoan truyền thống (Conventional Drilling - CD), khi mũi khoan tiến gần mặt đáy phôi, độ cứng vững của lớp kim loại đáy giảm nhanh. Ứng suất cắt và lực đẩy dọc trục ($F_t$) biến dạng uốn dẻo lớp kim loại này thành bavia hình vương miện (crown burr) hoặc bavia uốn (rollover burr).

Tiêu chí so sánh Khoan truyền thống (CD) Khoan hỗ trợ siêu âm cao tần (UVAD) Khoan rung động tần số thấp (LFVAD)
Cơ chế cắt Cắt liên tục, ma sát trượt lớn Cắt va đập siêu âm ($f > 20\text{ kHz}$) Cắt ngắt quãng tuần hoàn ($f < 1\text{ kHz}$)
Kiểm soát bavia Kém, phụ thuộc lớn vào độ mòn dao Triệt tiêu bavia tốt ($>65%$) Giảm bavia hiệu quả ($50-60%$)
Chi phí đầu tư bộ rung Thấp ($0\text{ USD}$) Rất cao ($>15{,}000\text{ USD}$, bộ PZT siêu âm) Trung bình ($1{,}500 - 3{,}000\text{ USD}$)
Độ ổn định hệ thống Cao Dễ quá nhiệt bộ biến năng Cao, dễ tích hợp trên máy CNC

Bảng ma trận yêu cầu hệ thống (MoSCoW)

  • Must have: Tái hiện chính xác hình học mũi khoan xoắn $\Phi 5\text{ mm}$, thiết lập chuẩn cơ tính Johnson-Cook cho Al6061-T6, giải thuật Explicit không phân kỳ.
  • Should have: Tối ưu lưới chuyển tiếp tại vùng tạo bavia (kích thước phần tử $L_e = 0.01\text{ mm}$), áp dụng Mass Scaling kiểm soát bước thời gian.
  • Could have: Tích hợp thuật ngữ lực cày (plowing force model) để nâng cao độ chính xác lực dọc trục.
  • Won't have (lần này): Khảo sát mài mòn dụng cụ cắt theo thời gian thực và trường nhiệt độ phi đẳng nhiệt.

Thiết kế hệ thống mô phỏng

Mô hình CAE hoàn chỉnh được xây dựng dựa trên ngăn xếp công nghệ:

  • Bộ tiền/hậu xử lý: Dassault Systèmes SIMULIA Abaqus/CAE v2023.
  • Solver: Abaqus/Explicit (tích hợp sai phân trung tâm - Central Difference Rule).
  • Phần mềm đo đạc & hiệu chuẩn hình ảnh: NIH ImageJ v1.54f.
  • Xử lý số liệu & Tối ưu hóa: Python v3.10 (thư viện NumPy, SciPy, Matplotlib) và Minitab v21.

Mô hình vật liệu Johnson-Cook cho Al6061-T6

Ứng suất chảy dẻo tương đương $\sigma_{eq}$ được định nghĩa: $$\sigma_{eq} = \left[ A + B(\varepsilon_p)^n \right] \left[ 1 + C \ln\left(\frac{\dot{\varepsilon}_p}{\dot{\varepsilon}_0}\right) \right] \left[ 1 - \left(T^*\right)^m \right]$$

Mô hình phá hủy dẻo (Ductile Damage) xác định độ biến dạng khi gãy $\varepsilon_f$: $$\varepsilon_f = \left[ d_1 + d_2 \exp\left(d_3 \frac{\sigma_m}{\sigma_{eq}}\right) \right] \left[ 1 + d_4 \ln\left(\frac{\dot{\varepsilon}_p}{\dot{\varepsilon}_0}\right) \right] \left[ 1 + d_5 T^* \right]$$

Thông số vật liệu Al6061-T6 chuẩn hóa trong đề tài:

  • Mật độ $\rho = 2700\text{ kg/m}^3$, Môđun đàn hồi $E = 68.9\text{ GPa}$, Hệ số Poisson $\nu = 0.33$.
  • $A = 324.1\text{ MPa}$, $B = 113.8\text{ MPa}$, $n = 0.42$, $C = 0.002$, $m = 1.34$.
  • Tham số phá hủy: $d_1 = -0.77$, $d_2 = 1.45$, $d_3 = -0.47$, $d_4 = 0$, $d_5 = 1.6$.

Phương pháp nghiên cứu và quy hoạch thực nghiệm

Sử dụng phương pháp quy hoạch thực nghiệm Taguchi với mảng trực quan $L_9(3^4)$ nhằm đánh giá mức độ ảnh hưởng của 4 yếu tố độc lập ở 3 mức:

Yếu tố công nghệ Ký hiệu Mức 1 (Thấp) Mức 2 (Trung bình) Mức 3 (Cao)
Vận tốc quay trục chính (rpm) $S$ 1000 2500 4000
Lượng chạy dao (mm/vòng) $F$ 0.02 0.05 0.08
Tần số rung (Hz) $f$ 100 500 1000
Biên độ rung ($\mu\text{m}$) $A$ 5 15 25

Tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio - $S/N$) theo đặc tính "Càng nhỏ càng tốt" (Smaller-is-Better) cho chiều cao bavia ($y_i$): $$\eta = -10 \log_{10} \left( \frac{1}{n} \sum_{i=1}^n y_i^2 \right)$$

Triển khai thực hiện và kết quả

Quy trình tính toán động học và mã nguồn mô phỏng

Trong quá trình khoan VAD, quỹ đạo dao động của lưỡi cắt theo phương trục $Z$ là tổng hợp của chuyển động tiến dao đều và dao động điều hòa: $$z(t) = A \sin(2\pi f t) + F \cdot n \cdot t$$ Vận tốc tức thời của mũi cắt: $$\dot{z}(t) = 2\pi f A \cos(2\pi f t) + F \cdot n$$

Khi $\dot{z}(t) < 0$, mũi khoan lùi về phía sau so với bề mặt gia công, tạo ra chu kỳ "tách - cắt" liên tục giúp giải phóng phoi và hạ thấp biến dạng uốn dẻo tại đáy lỗ.

# abaqus_vad_amplitude.py
# Script tạo bảng biên độ dao động điều hòa cho Reference Point mũi khoan
import numpy as np

def generate_vad_amplitude_curve(amplitude_um, frequency_hz, feed_rate_mm_rev, rpm, total_time, dt=1e-6):
    """
    Tạo dữ liệu biên độ dịch chuyển Z(t) = A*sin(2*pi*f*t) + F*(rpm/60)*t
    áp dụng cho Abaqus/Explicit Amplitude Table
    """
    A = amplitude_um * 1e-3  # Chuyển đổi sang mm
    f = frequency_hz          # Hz
    v_feed = feed_rate_mm_rev * (rpm / 60.0) # mm/s
    
    time_steps = np.arange(0, total_time, dt)
    z_displacement = A * np.sin(2 * np.pi * f * time_steps) + v_feed * time_steps
    
    with open("vad_disp_amplitude.amp", "w") as f_out:
        f_out.write("*AMPLITUDE, NAME=VAD_Z_MOTION, DEFINITION=TABULAR\n")
        for t, z in zip(time_steps, z_displacement):
            f_out.write(f"{t:.8e}, {z:.8e}\n")
            
    print(f"Generated {len(time_steps)} VAD kinematic steps successfully.")

# Tham số mô phỏng: A=15um, f=1000Hz, F=0.02mm/rev, S=4000rpm, T=0.005s
generate_vad_amplitude_curve(15.0, 1000.0, 0.02, 4000.0, 0.005)

Thuật toán sai phân trung tâm và kiểm soát bước thời gian

Abaqus/Explicit giải hệ phương trình cân bằng động tại nút: $$M \ddot{u}{(t)} = P{(t)} - I_{(t)} \implies \ddot{u}{(t)} = M^{-1} \left( P{(t)} - I_{(t)} \right)$$

Bước thời gian ổn định cực đại $\Delta t_{stable}$ được giới hạn bởi kích thước phần tử nhỏ nhất $L_e$ và tốc độ sóng giãn nở trong vật liệu $c_d$: $$\Delta t_{stable} \le \frac{L_e}{c_d} = L_e \sqrt{\frac{\rho}{E}} = 0.01\text{ mm} \times \sqrt{\frac{2700\text{ kg/m}^3}{68.9 \times 10^9\text{ Pa}}} \approx 1.98 \times 10^{-9}\text{ s}$$

Bằng việc áp dụng Mass Scaling cục bộ với hệ số tỷ lệ ổn định, bước thời gian giải được kiểm soát ở mức $\Delta t = 2.5 \times 10^{-8}\text{ s}$, đảm bảo tỷ số biến dạng không vượt quá tốc độ truyền sóng âm và hạn chế triệt để lỗi "excessively distorted elements".

Kiểm nghiệm và xác thực mô hình

Kết quả mô phỏng số được hiệu chuẩn chéo với dữ liệu thực nghiệm khoan hợp kim Al6061-T6 từ nghiên cứu của Chang et al. và Lauderbaugh:

Đại lượng kiểm nghiệm Điều kiện cắt ($S$ - $F$ - $f$ - $A$) Thực nghiệm (Exp) Mô phỏng số (FEM) Sai số đối đối (%)
Lực đẩy $F_t$ (CD) $1000\text{ rpm} - 0.05\text{ mm/vòng}$ $112.5\text{ N}$ $119.8\text{ N}$ $+6.49%$
Lực đẩy $F_t$ (VAD) $1000\text{ rpm} - 0.05\text{ mm/v} - 10\text{ kHz} - 10\mu\text{m}$ $94.5\text{ N}$ $101.2\text{ N}$ $+7.09%$
Chiều cao bavia đáy $R_f$ (CD) $4000\text{ rpm} - 0.02\text{ mm/vòng}$ $88.4\ \mu\text{m}$ $94.6\ \mu\text{m}$ $+7.01%$
Chiều cao bavia đáy $R_f$ (VAD) $4000\text{ rpm} - 0.02\text{ mm/v} - 1\text{ kHz} - 25\mu\text{m}$ $36.2\ \mu\text{m}$ $39.1\ \mu\text{m}$ $+8.01%$

Nhận xét: Độ sai lệch giữa mô phỏng và thực nghiệm luôn duy trì dưới $8.5%$, nằm trong ngưỡng tin cậy kỹ thuật cho phép ($<10%$).

Thrust Force Ft (N)
         Khoan thông thường (CD)       Khoan rung động (VAD)
         [Giảm ~16% lực cắt dọc trục nhờ rung động bổ trợ]

Kết quả đạt được và phân tích Taguchi

Kết quả mô phỏng 9 thí nghiệm theo mảng trực giao Taguchi $L_9$:

Thí nghiệm $S$ (rpm) $F$ (mm/vòng) $f$ (Hz) $A$ ($\mu\text{m}$) Chiều cao bavia thoát $H_{exit}$ ($\mu\text{m}$) $S/N$ Ratio bavia thoát (dB) Chiều cao bavia vào $H_{ent}$ ($\mu\text{m}$)
1 1000 0.02 100 5 52.3 -34.37 24.1
2 1000 0.05 500 15 68.7 -36.74 31.5
3 1000 0.08 1000 25 74.2 -37.41 38.2
4 2500 0.02 500 25 41.5 -32.36 19.8
5 2500 0.05 1000 5 55.4 -34.87 27.6
6 2500 0.08 100 15 89.1 -38.99 44.3
7 4000 0.02 1000 15 38.4 -31.68 17.2
8 4000 0.05 100 25 61.2 -35.73 29.4
9 4000 0.08 500 5 82.6 -38.34 41.8

Phân tích mức độ ảnh hưởng của các yếu tố (ANOVA Delta Values)

  • Chiều cao bavia thoát ($H_{exit}$): Yếu tố ảnh hưởng lớn nhất là Lượng chạy dao $F$ (Delta $S/N = 5.34\text{ dB}$, chiếm $54.2%$ tỷ trọng), tiếp theo là Tần số rung $f$ ($22.1%$), Biên độ rung $A$ ($14.5%$) và cuối cùng là Tốc độ trục chính $S$ ($9.2%$).
  • Tổ hợp thông số tối ưu để triệt tiêu bavia: $S_3 F_1 f_3 A_3$ ($S = 4000\text{ rpm}$, $F = 0.02\text{ mm/vòng}$, $f = 1000\text{ Hz}$, $A = 25\ \mu\text{m}$). Tại chế độ này, bavia thoát giảm từ $94.6\ \mu\text{m}$ (CD) xuống còn $38.4\ \mu\text{m}$ (VAD), tương đương mức giảm $59.41%$.

Đổi mới và đóng góp

  1. Ứng dụng thành công phương pháp tích phân hiện Explicit phi tuyến kết hợp giải thuật Mass Scaling cục bộ trên Abaqus để giải quyết bài toán biến dạng cực lớn và tách phoi động học trong gia công lai 1D-VAD mà không bị phân kỳ phần tử.
  2. Xác lập định lượng cơ chế giảm bavia: Rung động dọc trục tạo ra hiệu ứng tách dụng cụ ngắt quãng (intermittent cutting), làm giảm nhiệt ma sát và lực đẩy dọc trục $F_t$ trung bình đi 16%, đồng thời cưỡng bức phoi bị bẻ gãy sớm tại điểm kỳ dị tập trung ứng suất đáy lỗ, ngăn chặn sự phát triển của nếp gấp uốn dẻo (plastic hinge rollover).
  3. Cung cấp bộ thông số công nghệ chuẩn mực: Xác lập quy hoạch Taguchi định lượng hóa mối quan hệ giữa biên độ/tần số dao động và chiều cao bavia trên nhôm Al6061-T6, cung cấp cẩm nang thực nghiệm cho các kỹ sư chế tạo máy.
Giải pháp công nghệ Chiều cao bavia thoát trung bình Lực đẩy dọc trục $F_t$ Khả năng tự bẻ phoi Nhu cầu khử bavia phụ trợ
Khoan truyền thống (CD) $85 - 120\ \mu\text{m}$ (Gốc) $110 - 130\text{ N}$ Kém (Phoi dải dài) Bắt buộc (Tốn 20% chi phí)
Khoan phủ dao Nano (TiAlN) $70 - 95\ \mu\text{m}$ (Giảm 18%) $100 - 115\text{ N}$ Trung bình Vẫn cần xử lý mép
Mô hình VAD tối ưu (Đề tài) $35 - 45\ \mu\text{m}$ (Giảm 59.4%) $85 - 100\text{ N}$ Xuất sắc (Phoi vụn rời) Khử bỏ hoặc giảm 80% công đoạn

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng công nghiệp

  • Chế tạo cấu kiện vỏ và khung thân máy bay: Gia công hàng triệu lỗ tán đinh trên các tấm hợp kim nhôm Al6061-T6/Al7075. Việc giảm bavia giúp triệt tiêu hoàn toàn rủi ro tập trung ứng suất gây nứt mỏi (fatigue cracking) quanh lỗ đinh tán.
  • Sản xuất khối động cơ và hộp số ô tô: Khoan các đường dẫn dầu bôi trơn vi mô sâu bên trong thân máy, nơi các dụng cụ cạo bavia cơ học truyền thống không thể tiếp cận.
  • Gia công chi tiết cấy ghép y sinh (Medical Implants): Đảm bảo độ nhẵn thành lỗ và không có bavia bám dính gây viêm nhiễm mô sinh học.
                    QUY TRÌNH TRIỂN KHAI VAD TRÊN MÁY PHAY CNC

Phân tích hiệu quả kinh tế và ROI

  • Giảm chi phí vận hành: Khử bỏ công đoạn làm sạch bavia thủ công giúp một xưởng gia công 500,000 lỗ/năm tiết kiệm khoảng $18,000 - $25,000\text{ USD}$ chi phí nhân công và dụng cụ mài.
  • Tăng tuổi thọ dụng cụ cắt: Lực cắt giảm 16% và nhiệt độ vùng cắt giảm giúp kéo dài tuổi thọ mũi khoan thêm $35 - 45%$.
  • Thời gian hoàn vốn (ROI): Chi phí chế tạo đồ gá rung động tần số thấp khoảng $2,500\text{ USD}$, thời gian thu hồi vốn ước tính chỉ từ 4 đến 6 tháng.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  1. Giả thiết đẳng nhiệt: Chưa xét đến trường nhiệt cắt sinh ra do ma sát tại gờ mũi khoan, vốn có thể làm mềm cục bộ vật liệu nhôm ở vận tốc trục chính rất cao ($S > 6000\text{ rpm}$).
  2. Kích thước mẫu mô phỏng thu nhỏ: Do giới hạn tài nguyên trạm tính toán (Workstation RAM/CPU), chiều dày phôi mô phỏng giới hạn ở mức $0.1\text{ mm}$ thay vì chiều dày phôi công nghiệp thực tế ($2 - 5\text{ mm}$).

Hướng phát triển tiếp theo

  • Tích hợp mô hình nhiệt - cơ liên kết hoàn toàn (Coupled Temp-Displacement Explicit) để khảo sát ảnh hưởng của nhiệt độ cắt đến sự chuyển biến tổ chức tế vi (microstructure evolution) tại mép lỗ.
  • Mở rộng sang gia công các vật liệu khó gia công: Hợp kim Titan Ti-6Al-4V và thép không gỉ Austenitic AISI 316L.
  • Phát triển thuật toán học máy (Machine Learning): Sử dụng mạng nơ-ron nhân tạo (ANN) kết hợp thuật toán di truyền (GA) để tự động dự báo và điều khiển thích nghi biên độ rung thời gian thực.

Đối tượng hưởng lợi

  • Sinh viên & Học viên cao học Cơ khí: Nắm vững phương pháp luận xây dựng mô hình CAE động học phi tuyến, kỹ thuật chia lưới phần tử khử suy biến (C3D8R) và quy trình tối ưu Taguchi.
  • Kỹ sư R&D & Lập trình CAM/CNC: Sở hữu cơ sở dữ liệu chế độ cắt tối ưu cho gia công kết hợp rung động, ứng dụng trực tiếp vào tối ưu hóa năng suất nhà xưởng.
  • Doanh nghiệp gia công cơ khí chính xác: Tiết kiệm hàng trăm giờ lao động khử bavia, nâng cao tỷ lệ sản phẩm đạt chuẩn chất lượng xuất khẩu (chỉ số Cpk tăng $>1.33$).
  • Nhà nghiên cứu khoa học: Cung cấp mô hình tham chiếu được kiểm chứng thực nghiệm về cơ học biến dạng và phá hủy vật liệu theo mô hình Johnson-Cook trong môi trường Abaqus/Explicit.

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

1. Yêu cầu cấu hình phần cứng để chạy mô hình mô phỏng số này là gì?

Cấu hình tối thiểu: CPU Intel Core i7/Xeon 8 nhân 16 luồng, 32GB RAM DDR4/DDR5, ổ cứng NVMe SSD với dung lượng trống tối thiểu 50GB để lưu trữ file kết quả .odb. Cấu hình khuyến nghị: Máy trạm Workstation Dual Intel Xeon Gold, 64GB RAM ECC, card đồ họa chuyên dụng NVIDIA Quadro/RTX Ada để tối ưu hóa thời gian tính toán mô hình Explicit (khoảng 3–5 giờ cho một bài toán với 150,000 phần tử).

2. Tại sao tần số rung cao lại giúp giảm chiều cao bavia thoát?

Tần số rung cao ($f \ge 1000\text{ Hz}$) làm tăng số lần tách dụng cụ khỏi phôi trong một vòng quay của mũi khoan. Hiện tượng này chia nhỏ độ dày phoi tức thời thành các đoạn phoi vi mô, tạo ra trạng thái tải động chu kỳ làm tập trung biến dạng dẻo cục bộ và kích hoạt sớm các vết nứt tế vi theo tiêu chí phá hủy Johnson-Cook, khiến phoi đứt gãy trước khi kịp uốn cong thành bavia lớn.

3. Phương pháp này có thể áp dụng cho các vật liệu composite sợi carbon (CFRP) không?

Có thể áp dụng, tuy nhiên cơ chế tạo khuyết tật ở CFRP là tách lớp (delamination) và đứt gãy sợi chứ không phải biến dạng dẻo uốn bavia như kim loại. Khi đó, mô hình phần tử hữu hạn cần chuyển sang sử dụng tiêu chí phá hủy Hashin Damage hoặc phần tử Cohesive (CZM) thay cho mô hình Johnson-Cook.

4. Giải pháp chống lỗi biến dạng phần tử quá mức (Distorted Element Error) trong Abaqus/Explicit?

Cần áp dụng đồng thời: (1) Sử dụng phần tử lục giác tích phân giảm C3D8R với tính năng kiểm soát hiệu ứng phập phồng (Hourglass Control: Enhanced); (2) Bật tính năng xóa phần tử hư hỏng (Element Deletion) khi chỉ số biến dạng $D \ge 1$; (3) Tinh chỉnh bước thời gian giải thông qua hệ số tỷ lệ khối lượng (Mass Scaling Factor) hợp lý.

5. Chi phí đầu tư bộ gá rung động VAD trên máy CNC hiện hữu là bao nhiêu?

Chi phí chế tạo một đầu gá rung động tần số thấp (Low-Frequency VAD Tool Holder) dao động từ $1,500$ đến $3,000\text{ USD}$ (bao gồm bộ truyền động áp điện PZT, vòng tiếp điện không dây hoặc cổ góp điện, và mạch nguồn khuếch đại tín hiệu), thấp hơn nhiều so với hệ thống siêu âm công nghiệp chuyên dụng ($>15{,}000\text{ USD}$).

Kết luận

Đồ án tốt nghiệp của sinh viên Võ Duy Trung dưới sự hướng dẫn của PGS. TS. Phạm Huy Tuân đã giải quyết thành công bài toán cốt lõi trong gia công cơ khí chính xác: Làm sáng tỏ cơ chế hình thành và đề xuất giải pháp triệt tiêu bavia trong khoan kết hợp rung động (VAD) thông qua mô phỏng phần tử hữu hạn Abaqus/Explicit.

Với mô hình số đạt độ chính xác cao (sai số $<8.5%$ so với thực nghiệm), nghiên cứu chứng minh công nghệ VAD giúp giảm $16%$ lực đẩy dọc trục và giảm tới $59.41%$ chiều cao bavia thoát ở chế độ cắt tối ưu ($S = 4000\text{ rpm}, F = 0.02\text{ mm/vòng}, f = 1000\text{ Hz}, A = 25\ \mu\text{m}$). Kết quả này mở ra hướng ứng dụng đột phá cho các nhà máy gia công hàng không và thiết bị y tế, giảm thiểu chi phí hoàn thiện bề mặt và nâng cao chất lượng sản phẩm công nghiệp.