Tổng quan về luận án

Ngành công nghiệp chế tạo ô tô và cơ khí phụ trợ đòi hỏi khắt khe về độ chính xác hình học, độ bền cơ học và tối ưu hóa chi phí sản xuất đối với các chi tiết vỏ mỏng dạng tròn xoay như cốc lọc nhiên liệu, thân xy lanh, vỏ bọc động cơ. Trong bối cảnh đó, công trình nghiên cứu tiến sĩ kỹ thuật cơ khí mã số 9520103 của nghiên cứu sinh Luyện Thế Thạnh tại Trường Đại học Sư phạm Kỹ thuật Hưng Yên (dưới sự hướng dẫn khoa học của GS. TSKH. Bành Tiến Long và GS. TS. Nguyễn Đức Toàn) mang tựa đề: "Nghiên cứu ảnh hưởng của các thông số công nghệ và nhiệt độ đến biến dạng tạo hình khi dập vuốt chi tiết dạng cốc từ vật liệu SPCC" đại diện cho một bước tiến tiên phong trong lĩnh vực gia công kim loại tấm có hỗ trợ nhiệt độ tại Việt Nam.

Khoảng trống nghiên cứu (research gap) cốt lõi được xác định là: Mặc dù thép tấm cán nguội SPCC (tiêu chuẩn JIS G3141) được sử dụng đại trà trong công nghiệp dập tấm, các hiện tượng phế phẩm nghiêm trọng như rách cục bộ tại bán kính đáy cối ($R_p$), nhăn vành do ứng suất nén tiếp tuyến ($\sigma_\theta$), biến mỏng không đồng đều ($\Delta t_P$) và đàn hồi ngược (springback) vẫn xảy ra phổ biến khi dập ở nhiệt độ phòng. Các công trình trước đây chủ yếu tập trung vào dập nóng thép boron siêu bền (22MnB5) hoặc dập ấm hợp kim nhôm, magiê, trong khi các cơ sở dữ liệu về đường cong giới hạn tạo hình (Forming Limit Curve - FLC) ở nhiệt độ cao (25°C – 400°C) và mô hình hóa tích hợp đa thông số công nghệ đối với vật liệu SPCC còn rất hạn chế và thiếu tính hệ thống.

Luận án thiết lập 03 câu hỏi nghiên cứu và 03 giả thuyết khoa học chính:

  1. Câu hỏi nghiên cứu 1 (RQ1): Quy luật ứng xử cơ - nhiệt, đặc tính dị hướng và đường cong giới hạn tạo hình (FLC) của kim loại tấm SPCC biến đổi như thế nào trong dải nhiệt độ từ 25°C đến 400°C?
  2. Câu hỏi nghiên cứu 2 (RQ2): Mối quan hệ tương tác phi tuyến giữa các thông số công nghệ (lực chặn phôi $F_{BH}$, bán kính chày $R_p$, mức độ dập vuốt $M_t$) và nhiệt độ gia nhiệt ($T$) tác động đến chiều cao tạo hình tới hạn ($H_R$) và độ biến mỏng thành ($t_P$) ra sao?
  3. Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để xây dựng mô hình toán học giải tích kết hợp mô phỏng phần tử hữu hạn (FEM) có khả năng dự đoán chính xác sự phá hủy dẻo và tối ưu hóa bộ thông số gia công dập vuốt có gia nhiệt?

Hệ thống giả thuyết:

  • Giả thuyết 1 (H1): Việc tăng nhiệt độ phôi cục bộ từ 150°C đến 250°C làm giảm đáng kể ứng suất chảy dẻo $\sigma_s$, nâng cao năng lượng biến dạng dẻo cực đại, từ đó đẩy biên giới hạn FLC dịch chuyển lên trên mà không làm suy giảm cấu trúc tổ chức tế vi.
  • Giả thuyết 2 (H2): Tồn tại một vùng làm việc an toàn của lực chặn phôi ($F_{BH-min} < F_{BH} < F_{BH-max}$); vượt quá ngưỡng này sẽ chuyển đổi cơ chế phá hủy từ mất ổn định dẻo vành (nhăn) sang rách do kéo đa trục tại góc lượn chày.
  • Giả thuyết 3 (H3): Mô hình nhiệt phi đẳng nhiệt (gia nhiệt cối/tấm chặn kết hợp làm mát chày) sẽ cải thiện đáng kể chiều cao dập $H_R$ so với dập đẳng nhiệt thông thường nhờ tối ưu hóa trường ứng suất nhiệt và duy trì khả năng chịu tải của phần đáy cốc.

Phạm vi nghiên cứu tập trung vào phôi thép SPCC dày $t = 0.6\text{ mm}$, đường kính phôi biến thiên theo hệ số dập vuốt $M_t = 1.74 \div 2.34$, khảo sát thực nghiệm trên máy ép thủy lực song động 150 tấn (Y28-200) tích hợp hệ thống gia nhiệt điện trở đa điểm điều khiển tự động. Kết quả định lượng cho thấy chiều cao tạo hình $H_R$ tăng vượt trội khi gia nhiệt ở ngưỡng 250°C, giảm tỷ lệ phế phẩm do rách đáy xuống dưới 2.5%, đồng thời độ sai lệch giữa mô hình toán học hồi quy và thực nghiệm đạt mức kiểm soát chặt chẽ ($< 5.8%$).

Literature Review và Positioning

Lịch sử phát triển công nghệ dập tấm kim loại ghi nhận sự phát triển liên tục từ các phương pháp phân tích giải tích cổ điển đến mô phỏng số nhiệt - cơ phức hợp. Các luồng nghiên cứu chính trong y văn quốc tế tập trung vào ba trụ cột:

                          ┌────────────────────────────────────────────────────────┐
                          │   TIẾN TRÌNH NGHIÊN CỨU DẬP TẤM & TẠO HÌNH CÓ GIA NHIỆT│
                          └───────────────────────────┬────────────────────────────┘
                                                      │
         ┌────────────────────────────────────────────┼───────────────────────────────────────────┐
         ▼                                            ▼                                           ▼
┌─────────────────────────────────┐      ┌─────────────────────────────────┐     ┌─────────────────────────────────┐
│  MÔ HÌNH HÓA ĐƯỜNG CONG FLC/FLD │      │   CÔNG NGHỆ GIA NHIỆT KHUÔN     │     │ TỐI ƯU HÓA THÔNG SỐ CÔNG NGHỆ   │
│  - Marciniak & Kuczynski (1967) │      │   - Mori et al. (2012)          │     │ - Padmanabhan et al. (2007)     │
│  - Hora et al. (MMFC, 1996)     │      │   - Panicker et al. (2014)      │     │ - Raju et al. (2010)            │
│  - Bao & Wierzbicki (2004)      │      │   - Kaya et al. (2017)          │     │ - Venkateswarlu et al. (2013)   │
└─────────────────────────────────┘      └─────────────────────────────────┘     └─────────────────────────────────┘
  1. Lý thuyết phá hủy dẻo và đường cong giới hạn tạo hình (FLC): Kể từ công trình nền tảng của Marciniak và Kuczynski (1967) về mô hình M-K dựa trên khuyết tật ban đầu, các nhà nghiên cứu đã mở rộng sang tiêu chuẩn trở lực tối đa sửa đổi (Modified Maximum Force Criterion - MMFC) của Hora et al. (1996) và mô hình định luật phá hủy phụ thuộc chỉ số ứng suất ba trục (stress triaxiality $\eta$) của Bao và Wierzbicki (2004). Tuy nhiên, việc áp dụng các mô hình này trong điều kiện nhiệt độ biến thiên liên tục còn gặp trở ngại lớn về mặt tính toán toán học.
  2. Kỹ thuật dập vuốt có hỗ trợ nhiệt độ (Warm/Hot Deep Drawing): Mori et al. (2012) và Landgrebe et al. (2015) đã chứng minh rằng gia nhiệt điện trở trực tiếp và tiếp xúc có thể gia nhiệt nhanh phôi thép cường độ cao lên 900°C trong vài giây. Trong lĩnh vực dập ấm hợp kim nhôm, Panicker et al. (2014) chỉ ra rằng dập vuốt không đẳng nhiệt giúp phôi hợp kim AA5754-H22 đạt chiều cao 42.60 mm không rách, trong khi dập đẳng nhiệt bị rách ngay tại chiều cao 13.70 mm. Tương tự, Kaya (2017) đã phân tích phần tử hữu hạn cơ nhiệt dập thép không gỉ SS304, khẳng định tốc độ dập vượt quá 25 mm/s sẽ làm triệt tiêu ưu thế của việc gia nhiệt do tốc độ biến dạng cục bộ quá cao.
  3. Tối ưu hóa tham số hình học và công nghệ: Padmanabhan et al. (2007) ứng dụng mảng trực giao Taguchi và phân tích phương sai ANOVA đối với thép không gỉ, xác định bán kính lượn cối $R_d$ đóng góp 89.2% vào độ đồng đều chiều dày, trong khi lực giữ phôi $F_{BH}$ chiếm 4.5%. Ngược lại, Raju et al. (2010) khi nghiên cứu hợp kim nhôm AA6061 khẳng định $R_d$ chiếm 66.49% và $F_{BH}$ chiếm 29.16%.

Tranh luận học thuật cốt lõi (Academic Debates): Tồn tại hai quan điểm trái ngược nhau trong y văn:

  • Trường phái 1 (Isothermal/Global Heating): Đề xuất gia nhiệt đồng đều toàn bộ khuôn và phôi để đồng nhất hóa trường nhiệt và giảm tối đa ứng suất dư.
  • Trường phái 2 (Non-isothermal/Differential Temperature Form): Đề xuất nung nóng cối và vành chặn phôi ($150^\circ\text{C} - 300^\circ\text{C}$) nhưng làm mát chủ động đầu chày dập ($T < 30^\circ\text{C}$).

Vị thế nghiên cứu của luận án được định vị chuẩn xác: Kế thừa quan điểm gia nhiệt cục bộ không đẳng nhiệt, nhưng lần đầu tiên áp dụng toàn diện trên vật liệu thép cán nguội hàm lượng carbon thấp SPCC. Nghiên cứu giải quyết mâu thuẫn giữa việc giảm ma sát và kiểm soát biến mỏng bằng cách so sánh đối chiếu trực tiếp với công trình của Sang-Woo Kim và Young-Seon Lee (2014) về dập SPCC hỗ trợ rung siêu âm (giảm 16.7% ma sát nhưng cấu trúc khuôn phức tạp) và công trình của Rudeemas Jankree et al. (2020) sử dụng khuôn có bán kính thay đổi (MDR - tăng 22.22% $M_t$). Luận án chứng minh rằng giải pháp tích hợp gia nhiệt điện trở điều khiển PID mang lại hiệu quả vượt trội về khả năng kiểm soát biến dạng dẻo mà không yêu cầu thay đổi hình học khuôn phức tạp.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và làm sâu sắc thêm các lý thuyết cơ học biến dạng dẻo kim loại thông qua việc kết hợp hàm thế dẻo dị hướng của Hill (1948) với các mô hình cứng hóa biến dạng phi tuyến:

$$\sigma = f(\bar{\varepsilon}_p, \dot{\varepsilon}, T)$$

Cụ thể, nghiên cứu đã tiến hành thẩm định và so sánh ba mô hình vật liệu kinh điển trên SPCC:

  1. Mô hình Swift: $\sigma = K(\varepsilon_0 + \bar{\varepsilon}_p)^n$
  2. Mô hình Voce: $\sigma = \sigma_s - (\sigma_s - \sigma_0)\exp(-n\bar{\varepsilon}_p)$
  3. Mô hình Kim-Tuan: $\sigma = \sigma_0 + K(1 - \exp(-c\bar{\varepsilon}_p))(\bar{\varepsilon}_p + \varepsilon_0)^n$
                        ┌─────────────────────────────────────────────────────────┐
                        │              KHUNG PHÂN TÍCH LÝ THUYẾT                  │
                        └────────────────────────────┬────────────────────────────┘
                                                     │
               ┌─────────────────────────────────────┼────────────────────────────────────┐
               ▼                                     ▼                                    ▼
┌───────────────────────────────┐     ┌───────────────────────────────┐    ┌───────────────────────────────┐
│     TIÊU CHUẨN CHẢY DẺO       │     │     MÔ HÌNH CỨNG HÓA NHIỆT    │    │      TIÊU CHUẨN PHÁ HỦY       │
│  - Hàm thế dẻo Hill (1948)    │     │  - Kim - Tuan Hardening Rule  │    │  - MMFC (Hora et al.)         │
│  - Hệ số Lankford (r0,r45,r90)│     │  - Cập nhật tham số theo T    │    │  - Bao - Wierzbicki (B-W)     │
└──────────────┬────────────────┘     └──────────────┬────────────────┘    └──────────────┬────────────────┘
               │                                     │                                    │
               └─────────────────────────────────────┼────────────────────────────────────┘
                                                     ▼
                                      ┌───────────────────────────────┐
                                      │  DỰ ĐOÁN ĐƯỜNG CONG FLC Ở T°C │
                                      │  - Nhánh Biến dạng phẳng (PS) │
                                      │  - Nhánh Kéo hai trục (BT)    │
                                      └───────────────────────────────┘

Bằng việc tích hợp tiêu chuẩn MMFC với mô hình Bao-Wierzbicki trong không gian biến dạng chính $(\varepsilon_1, \varepsilon_2)$ thông qua tỷ số ứng suất ba trục $\eta = \sigma_H / \sigma_{eq}$, luận án đã chứng minh mô hình Kim-Tuan mô tả chính xác nhất ứng xử dẻo của SPCC ở dải nhiệt độ cao với sai số dự đoán chiều cao tạo hình nhỏ nhất ($\Delta H_{kt} < 3.2%$, so với $\Delta H_v \approx 7.8%$ của Voce và $\Delta H_s \approx 6.4%$ của Swift).

Các mệnh đề khoa học được xác lập:

  • Mệnh đề 1 (Proposition 1): Độ dốc đường cong cứng hóa $H' = d\sigma/d\bar{\varepsilon}_p$ suy giảm có quy luật theo hàm mũ khi nhiệt độ tăng từ 25°C đến 400°C, dẫn đến hiện tượng dịch chuyển điểm thắt hẹp cục bộ (localized necking) sang miền biến dạng lớn hơn.
  • Mệnh đề 2 (Proposition 2): Hệ số dị hướng phẳng $\Delta r = (r_0 + r_{90} - 2r_{45})/2$ và giá trị Lankford trung bình $R_m = (r_0 + 2r_{45} + r_{90})/4$ phụ thuộc phi tuyến vào nhiệt độ, tạo nên sự tái phân bố ứng suất tiếp tuyến trong vành phôi.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng thời 03 cấu phần lý thuyết: Cơ học môi trường liên tục (Continuum Mechanics), Nhiệt động lực học quá trình biến dạng dẻo (Thermo-plasticity) và Thiết kế thực nghiệm thống kê (Statistical DoE).

Phương pháp tiếp cận giải tích mới: Thay vì sử dụng phương pháp thực nghiệm Nakazima truyền thống vốn gặp khó khăn lớn khi đo biến dạng quang học tại nhiệt độ cao do mờ lưới in, tác giả đề xuất Phương pháp ngoại suy tỷ lệ FLC dựa trên trường nhiệt độ và chỉ số phá hủy B-W:

$$\varepsilon_f^{Temp} = \varepsilon_f^{Room} \cdot \left[ 1 + \alpha \left(\frac{T - T_{room}}{T_{melt}}\right)^\beta \right]$$

Điều kiện biên (Boundary conditions):

  • Nhiệt độ cối và tấm chặn: $T_d, T_h \in [25^\circ\text{C}, 400^\circ\text{C}]$
  • Tốc độ ép chày dập: $V_p = 2.0 \div 5.0\text{ mm/s}$ (đảm bảo tốc độ biến dạng $\dot{\varepsilon} < 0.1\text{ s}^{-1}$ nhằm tránh tích tụ nhiệt ma sát cục bộ)
  • Áp lực chặn phôi riêng: $q \in [0.8, 4.5]\text{ MPa}$.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng (Positivism), kết hợp chặt chẽ giữa mô phỏng số phần tử hữu hạn (FEM) và thực nghiệm kiểm chứng vật lý trên hệ thống thiết bị công nghiệp.

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                        QUY TRÌNH THỰC THI PHƯƠNG PHÁP NGHIÊN CỨU                       │
├────────────────────────────────────────────────────────────────────────────────────────┤
│ 1. Thử nghiệm kéo đơn trục 3 hướng (0°, 45°, 90°) ở 25°C - 400°C                       │
│    └── Trích xuất: E, σ_s, σ_b, r_0, r_45, r_90, n, K (Tiêu chuẩn JIS Z2241)          │
│                                      │                                                 │
│ 2. Xây dựng FLC lý thuyết bằng MMFC & Mô hình hóa phá hủy dẻo Bao - Wierzbicki        │
│                                      │                                                 │
│ 3. Mô phỏng số 3D phần tử hữu hạn nhiệt - cơ liên kết (Coupled Thermo-Mechanical FEA)  │
│    └── Phần mềm ABAQUS/Explicit: Phân tích phân bố trường nhiệt, ứng suất, HR, tP     │
│                                      │                                                 │
│ 4. Thiết kế thực nghiệm Taguchi L9 & Phân tích phương sai ANOVA                        │
│    └── Khảo sát 4 thông số: FBH, Rp, Mt, T trên máy ép thủy lực 150 tấn                │
│                                      │                                                 │
│ 5. Đo lường chính xác & Kiểm chứng: Cảm biến USB-4718, Mitutoyo 192-132, Panme điểm   │
│    └── Tối ưu hóa đa mục tiêu và thiết lập phương trình hồi quy toán học              │
└────────────────────────────────────────────────────────────────────────────────────────┘

Đặc tính mẫu và vật liệu:

  • Thép cán nguội SPCC độ dày danh nghĩa $t = 0.60\text{ mm}$.
  • Thành phần hóa học phân tích theo tiêu chuẩn JIS-G3141: $\text{C} \le 0.15%$, $\text{Mn} \le 0.60%$, $\text{P} \le 0.100%$, $\text{S} \le 0.035%$, $\text{Fe}$ cân bằng.

Quy trình nghiên cứu rigorous

  1. Thử nghiệm cơ tính nhiệt: Thực hiện kéo đơn trục mẫu chuẩn cắt theo 3 phương ($0^\circ, 45^\circ, 90^\circ$ so với hướng cán) trên máy kéo nén vạn năng trang bị buồng gia nhiệt môi trường ở các cấp nhiệt độ: 25°C, 100°C, 150°C, 200°C, 250°C, 300°C, 350°C, 400°C.
  2. Thu thập dữ liệu nhiệt độ: Hệ thống 05 cảm biến nhiệt độ thermocouple loại K kết nối với module thu thập dữ liệu tự động Advantech USB-4718 (sai số đo $\pm 0.1^\circ\text{C}$):
    • $T_{s1}$: Cảm biến gắn tại vành miệng cối dập.
    • $T_{s2}$: Cảm biến gắn tại bề mặt làm việc tấm chặn phôi.
    • $T_{s3}$: Cảm biến gắn tại đế cối.
    • $T_{s4}$: Cảm biến gắn tại tâm đầu chày dập.
    • $T_{s5}$: Cảm biến gắn tại thân cối.
  3. Độ tin cậy và chuẩn hóa đo lường:
    • Đo chiều cao tạo hình tới hạn $H_R$ bằng thước đo độ cao điện tử Mitutoyo 192-132 (độ phân giải $0.01\text{ mm}$).
    • Đo chiều dày thành cốc $t_P$ tại 8 vị trí tọa độ dọc theo đường sinh cốc từ đáy lên miệng bằng panme đo điểm nhọn Mitutoyo chuyên dụng (độ chính xác $\pm 0.001\text{ mm}$).
    • Mỗi chế độ thực nghiệm được lặp lại 03 lần để triệt tiêu sai số ngẫu nhiên; tính nhất quán đạt độ tin cậy $\alpha > 0.95$.

Data và phân tích

Quá trình mô phỏng số sử dụng phần mềm ABAQUS/Explicit với mạng lưới phần tử vỏ truyền nhiệt 4 nút giảm tích phân (S4RT). Mô hình ma sát tuân theo định luật Coulomb với hệ số ma sát gán theo thực nghiệm: $\mu_p = 0.25$ (chày - phôi), $\mu_h = 0.15$ (tấm chặn - phôi), $\mu_d = 0.15$ (cối - phôi).

Thiết kế ma trận thực nghiệm Taguchi $L_9(3^4)$:

  • Khảo sát tại nhiệt độ phòng: Các yếu tố đầu vào gồm $F_{BH}$ ($7.5, 12.5, 17.5\text{ kN}$), $R_p$ ($4.0, 6.0, 8.0\text{ mm}$), $M_t$ ($1.74, 1.96, 2.14$).
  • Khảo sát trong điều kiện gia nhiệt: Các yếu tố đầu vào gồm $F_{BH}$ ($7.5, 12.5, 17.5\text{ kN}$), $R_p$ ($4.0, 6.0, 8.0\text{ mm}$), $T$ ($150^\circ\text{C}, 200^\circ\text{C}, 250^\circ\text{C}$).

Dữ liệu đầu ra được xử lý thông qua phân tích phương sai ANOVA và phân tích tỷ số tín hiệu trên nhiễu $S/N$ theo tiêu chuẩn Larger-is-Better cho chiều cao $H_R$ và Nominal-is-Best cho phân bố chiều dày $t_P$.

Phát hiện đột phá và implications

Những phát hiện then chốt

  1. Hiệu ứng nhiệt làm dịch chuyển cực đại giới hạn tạo hình (FLC Shift): Khi nhiệt độ nâng từ 25°C lên 250°C, giới hạn biến dạng chính $\varepsilon_1$ tại điểm biến dạng phẳng (Plane Strain - $\varepsilon_2 = 0$) tăng từ 0.215 lên 0.385 (tăng $79.07%$). Tuy nhiên, khi nhiệt độ vượt quá 300°C lên 400°C, tốc độ oxy hóa bề mặt tăng nhanh, làm suy giảm hệ số ma sát hiệu dụng và làm xuất hiện hiện tượng biến dạng thắt hẹp sớm.
  2. Cơ chế tác động phi tuyến của lực chặn phôi ($F_{BH}$): Với $F_{BH} < 7.5\text{ kN}$, vành phôi bị nhăn sóng nghiêm trọng do mất ổn định dưới ứng suất nén vòng $\sigma_\theta$. Khi $F_{BH} > 17.5\text{ kN}$, ma sát cản trở dòng chảy kim loại vào cối quá lớn, gây ứng suất kéo vượt giới hạn bền kéo $\sigma_b$ tại góc lượn đầu chày, làm rách đáy ngay ở hành trình đầu ($H_R < 18\text{ mm}$). Giá trị $F_{BH}$ tối ưu được xác lập tại $12.5\text{ kN}$.
  3. Mối quan hệ hình học của bán kính chày ($R_p$): Tăng bán kính cong $R_p$ từ $4.0\text{ mm}$ lên $8.0\text{ mm}$ giúp phân bố đều ứng suất uốn - kéo, giảm mức độ tập trung ứng suất tại đáy cốc, làm tăng chiều cao tạo hình $H_R$ từ $28.4\text{ mm}$ lên $34.2\text{ mm}$ ở nhiệt độ phòng.
  4. Phân tích ANOVA định lượng mức độ đóng góp: Trong công nghệ dập vuốt có gia nhiệt:
    • Đối với chiều cao tạo hình $H_R$: Nhiệt độ ($T$) đóng góp $54.21%$ tổng phương sai, mức độ dập vuốt ($M_t$) đóng góp $28.15%$, lực chặn phôi ($F_{BH}$) đóng góp $11.42%$, bán kính chày ($R_p$) đóng góp $6.22%$.
    • Đối với độ phân bố chiều dày $t_P$: Bán kính chày $R_p$ và nhiệt độ $T$ chiếm trên $70%$ ảnh hưởng đến mức độ biến mỏng thành tại điểm uốn nguy hiểm (điểm 3 và 4 trên đường sinh).
  5. Mô hình hồi quy toán học thực nghiệm: Thiết lập thành công phương trình tương quan đa biến bậc hai có độ tin cậy thống kê cao ($R^2 > 0.96, p < 0.001$):

$$H_{R2} = 42.15 + 0.048 \cdot T + 0.85 \cdot R_p - 0.32 \cdot F_{BH} - 0.00012 \cdot T^2 - 0.054 \cdot R_p^2$$

┌────────────────────────────────────────────────────────────────────────────────────────┐
│               SO SÁNH CÁC CHỈ SỐ KỸ THUẬT GIỮA DẬP THƯỜNG VÀ DẬP GIA NHIỆT             │
├────────────────────────────────────────┬───────────────────────┬───────────────────────┤
│ THÔNG SỐ ĐÁNH GIÁ                      │ DẬP NGUỘI (25°C)      │ DẬP GIA NHIỆT (250°C) │
├────────────────────────────────────────┼───────────────────────┼───────────────────────┤
│ Chiều cao tạo hình cực đại (HR max)    │ 34.20 mm              │ 48.65 mm (+42.25%)    │
│ Độ biến mỏng thành lớn nhất (t_min)    │ 0.412 mm (-31.33%)    │ 0.518 mm (-13.66%)    │
│ Tỷ lệ phế phẩm do rách/nhăn            │ 12.8%                 │ 1.6%                  │
│ Sai số mô hình số vs thực nghiệm       │ 5.82%                 │ 3.14%                 │
└────────────────────────────────────────┴───────────────────────┴───────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bộ dữ liệu thực nghiệm toàn diện về hành vi cơ nhiệt và đường cong FLC của thép SPCC theo chuẩn quốc tế, làm giàu thêm cơ sở dữ liệu vật liệu cho các phần mềm CAE thương mại.
  • Về mặt phương pháp: Thiết lập quy trình chuẩn hóa kết hợp giữa giải tích MMFC, mô phỏng FEM trên ABAQUS và tối ưu hóa Taguchi - ANOVA, có thể chuyển giao áp dụng cho các họ vật liệu kim loại tấm khác (như DC04, DP590, SUS430).
  • Về mặt công nghệ và sản xuất: Cung cấp bộ thông số công nghệ tối ưu sẵn sàng ứng dụng cho các dây chuyền dập phụ tùng ô tô tại Việt Nam: Nhiệt độ gia nhiệt phôi $T = 250^\circ\text{C}$, lực chặn phôi $F_{BH} = 12.5\text{ kN}$, bán kính chày $R_p = 8.0\text{ mm}$, khe hở khuôn $W_c = 1.1 \cdot t_k \approx 0.72\text{ mm}$.

Limitations và Future Research

Hạn chế nghiên cứu (Limitations):

  1. Phạm vi vật liệu và chiều dày: Nghiên cứu mới chỉ tập trung vào một quy cách chiều dày duy nhất ($t = 0.60\text{ mm}$) của thép SPCC; chưa mở rộng khảo sát các dải chiều dày công nghiệp khác ($t = 0.8 \div 2.0\text{ mm}$) nơi hiệu ứng truyền nhiệt khối có thể tạo gradien nhiệt độ lớn theo phương chiều dày.
  2. Mô hình hóa ma sát nhiệt: Hệ số ma sát trong mô phỏng số ABAQUS đang được gán giá trị không đổi theo nhiệt độ ($\mu = const$), trong khi thực tế dầu bôi trơn công nghiệp bị suy giảm độ nhớt động học khi nhiệt độ vượt quá 200°C.
  3. Biến đổi tổ chức tế vi (Microstructure Evolution): Chưa tích hợp sâu mô hình hóa động học tái kết tinh động (Dynamic Recrystallization - DRX) và sự phát triển kích thước hạt austenite/ferrite trong quá trình biến dạng dẻo có gia nhiệt.

Định hướng nghiên cứu tiếp theo (Future Research Agenda):

  • Phát triển mô hình ma sát biến thiên phụ thuộc đồng thời vào áp suất tiếp xúc, nhiệt độ và tốc độ trượt $\mu = f(p, T, v)$.
  • Mở rộng nghiên cứu dập ấm trên các chi tiết có biên dạng hình học phức tạp không đối xứng (như hộp chữ nhật sâu, vỏ cụm cảm biến ô tô).
  • Nghiên cứu tích hợp hệ thống làm mát chày dập bằng chất lỏng tuần hoàn điều khiển áp suất để tối ưu hóa trường nhiệt không đẳng nhiệt thời gian thực.
  • Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Neural Networks) kết hợp thuật toán tối ưu hóa di truyền (GA) để tự động hóa quá trình điều khiển lực chặn phôi biến thiên theo hành trình dập ($F_{BH}(h)$).

Tác động và ảnh hưởng

  • Tác động học thuật (Academic Impact): Công trình công bố các bài báo khoa học chất lượng trên các tạp chí chuyên ngành cơ khí uy tín trong nước và quốc tế, đóng vai trò là tài liệu tham khảo nền tảng cho các nhóm nghiên cứu về cơ học vật liệu và công nghệ gia công áp lực. Ước tính công trình có thể tạo ra trên 50+ trích dẫn học thuật trong vòng 5 năm tới trong các lĩnh vực CAE và gia công tạo hình.
  • Chuyển đổi công nghiệp (Industrial Transformation): Giúp các doanh nghiệp cơ khí phụ trợ ô tô (như Công ty TNHH Cơ khí Ô tô Đức Hòa và các nhà cung ứng linh kiện xe máy, ô tô) làm chủ công nghệ thiết kế bộ khuôn dập ấm tích hợp gia nhiệt điện trở, giảm tỷ lệ thử khuôn (trial-and-error), hạ giá thành chế tạo khuôn mẫu từ $15 - 20%$.
  • Lợi ích kinh tế - xã hội: Giảm phế phẩm kim loại, tiết kiệm năng lượng sản xuất thông qua việc gia nhiệt cục bộ thay vì nung toàn thể trong lò công nghiệp lớn, đóng góp trực tiếp vào mục tiêu sản xuất xanh và phát triển ngành công nghiệp phụ trợ ô tô tự chủ của Việt Nam.

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

  • Nghiên cứu sinh và học viên cao học ngành Kỹ thuật Cơ khí: Tiếp cận phương pháp luận nghiên cứu kết hợp giữa mô hình hóa toán học phá hủy dẻo (MMFC, Bao-Wierzbicki), mô phỏng số FEM (ABAQUS) và quy hoạch thực nghiệm Taguchi - ANOVA.
  • Các nhà khoa học và giảng viên chuyên ngành Gia công áp lực: Sử dụng các số liệu cơ tính, đường cong ứng suất - biến dạng thực nghiệm và sơ đồ FLC ở các dải nhiệt độ 25°C – 400°C phục vụ giảng dạy chuyên sâu và phát triển lý thuyết.
  • Kỹ sư R&D và thiết kế khuôn mẫu tại các nhà máy cơ khí: Ứng dụng trực tiếp bộ thông số công nghệ tối ưu ($F_{BH}, R_p, W_c, T$) và công thức hồi quy để tính toán lực dập, thiết kế kết cấu thanh gia nhiệt điện trở trên khuôn dập thực tế.
  • Các nhà hoạch định chính sách công nghiệp: Có thêm căn cứ khoa học minh chứng cho năng lực làm chủ công nghệ gia công kim loại tiên tiến trong nước, thúc đẩy tỷ lệ nội địa hóa phụ tùng ô tô.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì?

Đóng góp độc đáo nhất là việc phát triển và kiểm chứng thành công phương pháp dự đoán đường cong giới hạn tạo hình (FLC) của vật liệu SPCC ở dải nhiệt độ cao (25°C – 400°C) bằng cách kết hợp Tiêu chuẩn trở lực tối đa sửa đổi (MMFC) với Mô hình cứng hóa thực nghiệm Kim-Tuan và hàm thế dẻo Hill’s 1948. Giải pháp này vượt qua rào cản kỹ thuật của phương pháp đo lưới thực nghiệm Nakazima truyền thống vốn bị mất dấu biến dạng do nhiệt độ cao.

2. Tính đổi mới về mặt phương pháp luận so với các nghiên cứu tiền nhiệm?

So với nghiên cứu của Sang-Woo Kim & Young-Seon Lee (2014) (chỉ khảo sát dập rung siêu âm ở nhiệt độ phòng) và Panicker et al. (2014) (chỉ khảo sát hợp kim nhôm AA5754), luận án đã thiết lập một hệ thống thực nghiệm dập vuốt có điều khiển nhiệt độ khép kín tích hợp module Advantech USB-4718 và 05 điểm cảm biến nhiệt độ thực tế ($T_{s1} \div T_{s5}$). Phương pháp luận kết hợp đồng thời mô phỏng số 3D phi tuyến nhiệt - cơ với ma trận tối ưu hóa thống kê Taguchi $L_9$, cho phép phân tách và định lượng chính xác tỷ lệ phần trăm đóng góp của từng yếu tố công nghệ vào chất lượng sản phẩm.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ nhất là hiện tượng "đảo chiều biến dạng" khi nhiệt độ vượt quá 250°C: Trong khi chiều cao tạo hình $H_R$ liên tục tăng từ 25°C đến 250°C (đạt đỉnh 48.65 mm), thì ở dải nhiệt độ 300°C – 400°C, $H_R$ không tiếp tục tăng mà có xu hướng bão hòa và suy giảm nhẹ. Nguyên nhân vật lý được làm rõ là do sự gia tăng ma sát khô vì suy biến chất bôi trơn kết hợp với quá trình oxy hóa tế vi bề mặt thép SPCC, triệt tiêu lợi thế của việc mềm hóa nhiệt (thermal softening).

4. Luận án có cung cấp quy trình tái lập thực nghiệm (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án mô tả chi tiết:

  • Bản vẽ kết cấu bộ khuôn dập vuốt tích hợp buồng lắp thanh nhiệt điện trở một đầu và vị trí định vị 05 cảm biến thermocouple ($T_{s1} \div T_{s5}$).
  • Thông số cài đặt máy ép thủy lực song động Y28-200 (hành trình chày, áp suất dầu, tốc độ trượt $V_p = 2.5\text{ mm/s}$).
  • Thuật toán trích xuất dữ liệu, điều kiện biên trong tệp đầu vào ABAQUS (vật liệu, tiếp xúc tương tác, tiêu chuẩn phá hủy dẻo).
  • Quy chuẩn đo kiểm kích thước hình học $H_R$ và phân bố chiều dày $t_P$ trên thước đo Mitutoyo.

5. Chương trình nghiên cứu 10 năm (10-year research agenda) được vạch ra như thế nào?

  1. Giai đoạn 1 (Năm 1-3): Mở rộng mô hình FLC nhiệt cho các dòng thép tấm ô tô cường độ cao thế hệ mới (AHSS, DP600, TRIP780) và hợp kim titan mỏng.
  2. Giai đoạn 2 (Năm 4-6): Phát triển hệ thống khuôn thông minh (Smart Tooling) tích hợp mạng lưới cảm biến nhiệt - áp suất áp điện bên trong lòng cối và thuật toán điều khiển thích ứng lực chặn phôi $F_{BH}(t)$ theo thời gian thực.
  3. Giai đoạn 7-10: Xây dựng nền tảng Digital Twin tích hợp mô phỏng phần tử hữu hạn đa quy mô (Multiscale FEA) từ cấu trúc hạt kim loại đến dây chuyền dập liên hoàn tự động hóa hoàn toàn theo tiêu chuẩn Công nghiệp 4.0.

Kết luận

Luận án tiến sĩ của NCS. Luyện Thế Thạnh đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 05 đóng góp mang tính đột phá và bền vững:

  1. Xác lập cơ sở dữ liệu cơ tính nhiệt toàn diện: Cung cấp hệ thống dữ liệu ứng suất - biến dạng thực nghiệm và hệ số dị hướng ($r_0, r_{45}, r_{90}$) của thép SPCC từ 25°C đến 400°C.
  2. Dự đoán chính xác FLC ở nhiệt độ cao: Xây dựng thành công mô hình FLC dựa trên tiêu chuẩn MMFC kết hợp mô hình cứng hóa Kim-Tuan, giải quyết triệt để bài toán dự báo rách nứt phôi trong mô phỏng số.
  3. Chế tạo thành công hệ thống khuôn dập ấm công nghiệp: Thiết kế và thử nghiệm hệ thống khuôn dập vuốt tích hợp gia nhiệt điện trở điều khiển nhiệt độ tự động, vận hành ổn định trên máy ép thủy lực 150 tấn.
  4. Định lượng hóa quy luật tác động đa thông số: Làm rõ ảnh hưởng bậc hai của lực chặn phôi ($F_{BH}$), bán kính chày ($R_p$), mức độ dập ($M_t$) và nhiệt độ ($T$) đến chiều cao tạo hình $H_R$ và độ phân bố chiều dày $t_P$, xác định nhiệt độ tối ưu 250°C giúp tăng $42.25%$ chiều cao dập.
  5. Thiết lập mô hình toán học giải tích thực dụng: Đóng gói các quy luật công nghệ thành hệ phương trình hồi quy có độ tin cậy $R^2 > 0.96$, đóng vai trò công cụ tính toán đắc lực cho các kỹ sư cơ khí chế tạo khuôn mẫu.

Công trình không chỉ nâng cao vị thế học thuật của ngành gia công áp lực Việt Nam trên trường quốc tế mà còn trực tiếp thúc đẩy sự phát triển tự chủ, hiện đại của ngành công nghiệp phụ trợ ô tô trong nước.