Tổng quan về luận án

Trong ngành công nghiệp đóng tàu hiện đại, kim loại tấm đóng vai trò là vật liệu cơ bản cốt lõi khi "chiếm hơn 90% trọng lượng của con tàu". Hầu hết kết cấu thân tàu thủy—từ sống chính, vách chống va, khung bệ máy cho đến hệ thống bích tăng cứng hầm hàng và cửa khoang—đều đòi hỏi các tấm thép có biên dạng cong phức tạp với chiều dày danh nghĩa lớn ($t \ge 5\text{ mm}$). Theo thống kê từ các công trình nghiên cứu cơ học tạo hình, "hơn 55% các sản phẩm làm từ kim loại tấm được sản xuất bằng phương pháp uốn trong các ngành công nghiệp này". Tuy nhiên, quá trình uốn dập nguội truyền thống đối với thép tấm kết cấu thường đối mặt với hai thách thức kỹ thuật lớn: lực tạo hình yêu cầu quá lớn đòi hỏi các cụm máy chấn thủy lực công suất siêu trường (như dòng máy 1200 tấn tại Công ty Đóng tàu Nam Triệu), và hiện tượng đàn hồi lại (springback) nghiêm trọng sau khi dỡ tải, dẫn đến sai số biên dạng hình học và gia tăng ứng suất dư trong liên kết hàn vỏ tàu.

Khoảng trống nghiên cứu (research gap) cốt lõi mà luận án tập trung giải quyết là sự thiếu hụt một mô hình ứng xử cơ - nhiệt phi tuyến hoàn chỉnh, có khả năng tích hợp hiệu ứng Bauschinger và quy luật biến cứng kết hợp (combined hardening) ở các dải nhiệt độ trung gian ($300^\circ\text{C}$ và $600^\circ\text{C}$) nhằm dự đoán chính xác lực tạo hình và góc đàn hồi lại khi uốn tấm thép kết cấu SS400 được trợ động bằng cảm ứng điện từ (Induction heating - Assisted Machining - IAM). Các nghiên cứu trước đây đa phần chỉ khảo sát kim loại tấm mỏng ($t < 1.5\text{ mm}$), vật liệu màu hoặc chỉ áp dụng quy luật biến cứng đẳng hướng thuần túy vốn đánh giá sai lệch đáng kể hiện tượng đàn hồi ngược.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Nhiệt độ nung phôi cục bộ bằng cảm ứng điện từ làm thay đổi đường cong ứng suất - biến dạng dẻo thực và các thông số hóa bền của thép kết cấu SS400 như thế nào?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để nhận dạng và hiệu chỉnh chính xác các tham số của mô hình biến cứng kết hợp (Chaboche combined hardening) phụ thuộc nhiệt độ cho thép SS400 nhằm tích hợp vào phần mềm mô phỏng phần tử hữu hạn ABAQUS/CAE?
  • Câu hỏi nghiên cứu 3 (RQ3): Mối quan hệ tương tác và mức độ ảnh hưởng của tổ hợp các thông số công nghệ (nhiệt độ nung $T$, bán kính chày $R_{ch}$, hành trình chày $H$) đến lực uốn đỉnh $P$ và góc đàn hồi $\theta_{SP}$ tuân theo quy luật toán học nào?
  • Giả thuyết 1 (H1): Việc nung nóng cảm ứng điện từ đến $600^\circ\text{C}$ sẽ làm suy giảm giới hạn chảy $\sigma_y$ của thép SS400 trên 40%, từ đó giảm trực tiếp lực uốn danh định và triệt tiêu đáng kể năng lượng phục hồi đàn hồi sau uốn.
  • Giả thuyết 2 (H2): Mô hình biến cứng kết hợp kết hợp giữa quy luật biến cứng đẳng hướng Voce và biến cứng động học Chaboche sẽ giảm sai số dự đoán góc đàn hồi trong mô phỏng số xuống dưới 5% so với thực nghiệm uốn chữ V.

Khung lý thuyết của công trình được xây dựng trên nền tảng cơ học biến dạng dẻo phi tuyến (Nonlinear Continuum Mechanics), tiêu chuẩn chảy dẻo von Mises, lý thuyết truyền nhiệt và dòng điện xoáy cảm ứng Foucault trong vật dẫn điện từ. Phạm vi thực nghiệm của nghiên cứu được kiểm chứng chặt chẽ trên mẫu thép tấm cán nóng SS400 kích thước $110\text{ mm} \times 35\text{ mm} \times 5\text{ mm}$, gia công dọc theo hướng cán, khảo sát tại ba dải nhiệt độ ($32^\circ\text{C}$, $300^\circ\text{C}$, $600^\circ\text{C}$), ba mức bán kính chày ($R_{ch} = 10, 15, 20\text{ mm}$) và ba hành trình uốn ($H = 10, 16, 22\text{ mm}$). Kết quả nghiên cứu tạo nền tảng chuyển giao công nghệ đột phá cho việc chế tạo các cụm chi tiết tàu thủy đòi hỏi độ chính xác cao.

Literature Review và Positioning

Nghiên cứu cơ bản về công nghệ uốn và trợ lực nhiệt trong tạo hình kim loại tấm đã trải qua nhiều giai đoạn phát triển với các trường phái tiếp cận khác nhau. Dòng nghiên cứu thứ nhất tập trung vào việc mô tả hiện tượng đàn hồi lại sau biến dạng dẻo ở nhiệt độ phòng. Các công trình kinh điển của Wagoner et al. (2013) và Mike Gedeon (2010) đã chứng minh rằng springback là hệ quả tất yếu của sự phân bố mômen uốn đàn hồi còn lại dọc theo chiều dày tấm khi dỡ tải, chịu chi phối trực tiếp bởi tỷ số giữa môđun đàn hồi Young ($E$) và giới hạn chảy ($\sigma_y$). Fuh-Kuo Chen & Shen-Fu Ko (2007) khi phân tích cơ học uốn chữ L đã chỉ ra sự phân vùng biến dạng phức tạp, trong đó vùng tiếp xúc đáy chày và vùng sức căng tự do có xu hướng tạo ra các mômen đàn hồi ngược chiều nhau, khiến việc kiểm soát kích thước hình học trở nên bất khả kháng nếu không có giải pháp can thiệp cơ - nhiệt đồng thời.

Dòng nghiên cứu thứ hai khai thác các nguồn nhiệt trợ động nhằm làm mềm hóa vật liệu (thermal softening). Peter B. Christian & Frank (2005) đã tiến hành thử nghiệm uốn nhiệt đối với thép kết cấu Grade 50 và Grade 70 trong chế tạo dầm cầu, chứng minh việc gia nhiệt cục bộ giúp giảm lực uốn rõ rệt nhưng cảnh báo nguy cơ suy giảm độ dẻo nếu nhiệt độ vượt ngưỡng $1200^\circ\text{F}$ ($648^\circ\text{C}$). Nghiên cứu của Stachowicz et al. (2012) trên thép tấm không gỉ AMS 5604 bằng công nghệ nung cao tần xác nhận giới hạn chảy giảm từ $904\text{ MPa}$ ở $20^\circ\text{C}$ xuống mức thấp ở nhiệt độ cao, đồng thời góc đàn hồi giảm tuyến tính theo nhiệt độ. Tương tự, Ailing Wang, Jianguo Lin, Trevor A. Dean et al. (2012) nghiên cứu uốn ấm hợp kim nhôm AA5754 ($20^\circ\text{C} - 450^\circ\text{C}$) và Nguyen Duc-Toan et al. (2013) khảo sát uốn chữ V hợp kim magie AZ31 đã cùng chứng minh rằng nhiệt độ làm suy giảm đáng kể dị hướng vật liệu và triệt tiêu năng lượng đàn hồi tích lũy.

Trong nước, các công trình của Trần Hải Đăng (2016) về công nghệ lăn ép vỏ tàu mới chỉ dừng lại ở điều kiện nguội; Đặng Thị Hồng Huế (2017) khảo sát cán nêm ngang thép C45 hay Nguyễn Thành Huân (2018) nghiên cứu tiện trợ động laser thép 9XC. Mặc dù các phương pháp gia nhiệt như Laser (LAM), Plasma (PEM), lò nung (FAM) và dòng điện tiếp xúc (EAM) đều có những ưu thế nhất định, gia nhiệt cảm ứng điện từ (IAM) được chứng minh là phương thức tối ưu nhất cho uốn tấm dày nhờ khả năng sinh nhiệt thể tích nội tại nhanh, hiệu suất truyền năng lượng cao và dễ tích hợp vào khuôn chấn.

                    ┌─────────────────────────────────────────────────────────┐
                    │               Literature Positioning Space              │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
         ┌───────────────────────────────────────┴───────────────────────────────────────┐
         ▼                                                                               ▼
┌─────────────────────────────────┐                                           ┌─────────────────────────────────┐
│     Cold Forming & Pure         │                                           │    Thermal Forming & General    │
│     Isotropic Hardening         │                                           │    Constitutive Formulations    │
│  - Wagoner et al. (2013)        │                                           │  - Peter B. Christian (2005)    │
│  - Chen & Ko (2007)             │                                           │  - Stachowicz et al. (2012)     │
│  - Mike Gedeon (2010)           │                                           │  - Nguyen Duc-Toan et al. (2013)│
│  * Limitation: High springback  │                                           │  * Limitation: Isotropic / J-C  │
│    error (>12%), no thermal aid │                                           │    models lack reversal capture │
└────────────────┬────────────────┘                                           └────────────────┬────────────────┘
                 │                                                                             │
                 └───────────────────────────────────────┬─────────────────────────────────────┘
                                                         │
                                                         ▼
                                      ┌─────────────────────────────────────┐
                                      │        THIS DISSERTATION GAP        │
                                      │  Induction-Assisted V-Bending of    │
                                      │  Heavy SS400 Ship Plate (t=5mm)     │
                                      │  via Thermo-Coupled Combined Voce-  │
                                      │  Chaboche Hardening in ABAQUS/CAE   │
                                      └─────────────────────────────────────┘

Điểm xung đột học thuật cốt lõi trong y văn nằm ở việc lựa chọn mô hình ứng xử vật liệu: hầu hết các tác giả quốc tế (như Nguyen Duc-Toan et al., 2013) áp dụng mô hình Johnson-Cook thuần túy hoặc mô hình biến cứng đẳng hướng Voce để mô phỏng uốn nhiệt. Tuy nhiên, khi uốn chữ V, thớ ngoài của tấm chịu kéo trong khi thớ trong chịu nén, và trong quá trình dỡ tải diễn ra sự đổi chiều ứng suất nghiêm trọng. Việc bỏ qua ứng suất ngược (backstress $\alpha$) và hiệu ứng Bauschinger trong mô hình đẳng hướng khiến sai số dự đoán đàn hồi thực tế bị phóng đại. Luận án đã định vị chính xác vị trí tiên phong bằng việc tích hợp quy luật biến cứng kết hợp (Chaboche combined hardening) phụ thuộc nhiệt độ cho thép SS400 dày $5\text{ mm}$, vượt trội hơn hẳn so với các công trình thực nghiệm đơn lẻ trước đây.

Đó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à hoàn thiện lý thuyết biến dạng dẻo phi tuyến ở dải nhiệt độ trung gian thông qua việc thiết lập phương thức tổng hợp mới để xác định các tham số của mô hình biến cứng kết hợp (combined isotropic-kinematic hardening model), kết hợp giữa quy luật bão hòa đẳng hướng của Voce và quy luật động học phi tuyến của Chaboche.

Trong cơ học biến dạng, mặt chảy dẻo von Mises mở rộng có dạng: $$f(\boldsymbol{\sigma} - \boldsymbol{\alpha}) - \sigma_{iso} = 0$$ trong đó $\boldsymbol{\alpha}$ là tenxơ ứng suất ngược biểu diễn sự dịch chuyển tâm mặt chảy dẻo (biến cứng động học), và $\sigma_{iso}$ là bán kính mặt chảy dẻo biểu diễn sự giãn nở đồng đều (biến cứng đẳng hướng). Luận án xác lập sự tiến triển của ứng suất ngược đơn trục $\alpha$ theo quy luật Chaboche vi phân: $$\mathrm{d}\alpha = C \frac{1}{\sigma_0} (\sigma - \alpha) \mathrm{d}\bar{\varepsilon}^p - \gamma \alpha \mathrm{d}\bar{\varepsilon}^p$$ với $C$ là môđun hóa bền động học ban đầu, $\gamma$ là thông số kiểm soát tốc độ suy giảm môđun hóa bền khi biến dạng dẻo tương đương $\bar{\varepsilon}^p$ tăng. Đồng thời, bán kính mặt chảy đẳng hướng được mô tả qua quy luật Voce: $$\sigma_{iso} = \sigma_0 + A \left(1 - e^{-B \bar{\varepsilon}^p}\right)$$ với $\sigma_0$ là giới hạn chảy ban đầu, $A$ là ứng suất bão hòa cực đại, và $B$ là hệ số tốc độ mở rộng miền đàn hồi.

                    ┌─────────────────────────────────────────────────────────┐
                    │               Combined Hardening Framework              │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
                     ┌───────────────────────────┴───────────────────────────┐
                     ▼                                                       ▼
      ┌─────────────────────────────┐                         ┌─────────────────────────────┐
      │   Isotropic Evolution       │                         │   Kinematic Backstress      │
      │   (Voce Saturation Law)     │                         │   (Chaboche Formulation)    │
      │ σ_iso = σ_0 + A(1 - e^-Bε)  │                         │  dα = C(dε) - γ·α·(dε)      │
      └──────────────┬──────────────┘                         └──────────────┬──────────────┘
                     │                                                       │
                     └───────────────────────────┬───────────────────────────┘
                                                 │
                                                 ▼
                              ┌─────────────────────────────────────┐
                              │     Unified Thermo-Coupled Model    │
                              │   Parameters {C, γ, A, B, σ_0}(T)   │
                              │   Calibrated at 32°C, 300°C, 600°C  │
                              └─────────────────────────────────────┘

Đóng góp lý thuyết mang tính đột phá của luận án là chứng minh được rằng ở điều kiện nhiệt độ cao ($300^\circ\text{C}$ và $600^\circ\text{C}$), thông số biến cứng động học $C$ và ứng suất ngược $\alpha$ suy giảm phi tuyến theo nhiệt độ nhanh hơn nhiều so với thành phần biến cứng đẳng hướng. Việc lượng hóa thành công bộ 5 tham số ${C, \gamma, \sigma_0, A, B}$ tại từng mốc nhiệt độ cho phép lý giải bản chất vật lý của sự triệt tiêu hiệu ứng Bauschinger khi chuyển pha nhiệt, đóng góp một mắt xích quan trọng vào lý thuyết cơ học tính toán vật liệu tấm.

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

Khung phân tích của luận án tích hợp chặt chẽ ba trụ cột lý thuyết:

  1. Lý thuyết cảm ứng điện từ và hiệu ứng bề mặt (Skin effect): Xác định sự phân bố trường nhiệt thông qua chiều sâu thâm nhập của dòng xoáy Foucault: $$\delta = 503 \sqrt{\frac{\rho}{\mu_r f}} \quad (\text{mm})$$ cho phép thiết kế tối ưu cuộn cảm ứng xoắn ốc bao quanh phôi, đảm bảo gia nhiệt đồng đều tiết diện $5\text{ mm}$ trước khi tiến hành uốn.
  2. Quy trình tối ưu hóa ngược (Inverse FEA Optimization): Sử dụng dữ liệu thử kéo đơn trục ở nhiệt độ cao để xác định sơ bộ đường cong $\sigma - \varepsilon$, sau đó chạy mô phỏng biến dạng chu trình ảo trong ABAQUS/CAE để nội suy chính xác ứng suất ngược $\alpha$, loại bỏ sự cần thiết của các phép thử nén chu trình vốn cực kỳ phức tạp và dễ mất ổn định uốn dọc trên mẫu tấm mỏng.
  3. Mô hình tương tác đa yếu tố thực nghiệm (Empirical Multi-Response Formulation): Kết hợp phân tích hồi quy phi tuyến với quy hoạch thực nghiệm Taguchi để xây dựng 4 hàm truyền toán học dự báo chính xác: Lực tạo hình ($P_{TC}$), Góc biến dạng đàn hồi ($SP_C$), Góc uốn chi tiết ($\theta_C$), và Bán kính trong sản phẩm ($R_{UC}$).

Điều kiện biên của khung phân tích được xác lập nghiêm ngặt: quá trình uốn dập được coi là đoạn nhiệt cục bộ trong khoảng thời gian tạo hình cực ngắn ($< 5\text{ giây}$), biến dạng dẻo chủ yếu diễn ra theo trạng thái biến dạng phẳng dọc theo đường sinh uốn.

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ủ nghiêm ngặt thế giới quan thực chứng (positivism) kết hợp cơ học tính toán thực nghiệm (computational-empirical realism). Thiết kế nghiên cứu đa cấp độ (multi-level design) được triển khai qua ba pha liên kết mật thiết:

  • Pha 1 (Vi mô/Mẫu thử): Thử nghiệm kéo tĩnh ở dải nhiệt độ $32^\circ\text{C}$, $300^\circ\text{C}$, $600^\circ\text{C}$ để xác định cơ tính cơ bản ($E, \sigma_y, \sigma_b, \delta%$).
  • Pha 2 (Meso/Mô phỏng phần tử hữu hạn): Xây dựng mô hình 2D/3D phần tử hữu hạn trên ABAQUS, thiết lập mạng lưới phần tử liên tục phi tuyến (C3D8R và CPS4R) với mật độ chia lưới mịn tại vùng góc uốn, cài đặt các điều kiện tiếp xúc ma sát phạt (Penalty contact, $\mu = 0.15$), mô phỏng chính xác giai đoạn gia tải chấn uốn và dỡ tải tự do (springback step).
  • Pha 3 (Vĩ mô/Thực nghiệm công nghệ): Triển khai ma trận uốn kiểm chứng trên máy chuyên dụng có tích hợp nguồn nung cảm ứng cao tần, thu thập lực uốn và đo đạc góc đàn hồi thực tế.
┌──────────────────────────────────────────────────────────────────────────────────┐
│                           Rigorous Research Workflow                             │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
 ┌───────────────────────────────────────┼───────────────────────────────────────┐
 │                                       │                                       │
 ▼                                       ▼                                       ▼
┌─────────────────────────┐   ┌─────────────────────────┐   ┌─────────────────────────┐
│       PHASE 1           │   │       PHASE 2           │   │       PHASE 3           │
│ Specimen Tensile Tests  │   │ Non-Linear FEA Modeling │   │ Taguchi Matrix (L9) &   │
│ - HUNGTA H-200A (200kN) │   │ - ABAQUS/CAE (C3D8R)    │   │ Physical V-Bending      │
│ - Induction heating     │   │ - Inverse hardening fit │   │ - Autonics PID control  │
│ - T = 32, 300, 600°C    │   │ - Loading & Unloading   │   │ - Force & Angle readout │
└────────────┬────────────┘   └────────────┬────────────┘   └────────────┬────────────┘
             │                             │                             │
             └─────────────────────────────┼─────────────────────────────┘
                                           │
                                           ▼
                    ┌─────────────────────────────────────────────┐
                    │      Triangulation & Multi-Response         │
                    │      Mathematical Optimization (ANOVA)      │
                    └─────────────────────────────────────────────┘

Quy trình nghiên cứu rigorous

Toàn bộ phôi thử nghiệm đều được gia công từ thép tấm cán nóng SS400 đạt chuẩn tiêu chuẩn công nghiệp (thành phần hóa học danh định: $\text{C} \le 0.20%$, $\text{Si} \le 0.35%$, $\text{Mn} \le 1.40%$, $\text{P} \le 0.045%$, $\text{S} \le 0.045%$). Phôi tấm được cắt bằng máy cắt dây CNC GS-4050B có độ chính xác định vị $\pm 0.005\text{ mm}$ dọc theo phương cán để loại trừ ảnh hưởng của ứng suất nhiệt cắt thô và đảm bảo tính nhất quán về cơ lý tính.

Hệ thống thiết bị đo đạc và kiểm chuẩn bao gồm:

  • Máy kéo nén vạn năng HUNGTA H-200A (dải tải trọng $0 - 200\text{ kN}$).
  • Bộ gia nhiệt cảm ứng điện từ cao tần tùy biến, tích hợp cuộn cảm ứng xoắn ốc chế tạo bằng ống đồng đỏ làm mát tuần hoàn bằng nước.
  • Bộ điều khiển nhiệt độ chuyên dụng Autonics TZN4L với cảm biến nhiệt độ ngẫu nhiệt loại K (sai số $\le \pm 1^\circ\text{C}$), đảm bảo duy trì nhiệt độ phôi ổn định chính xác tại các mốc nhiệt thử nghiệm.
  • Bộ chày cối uốn chữ V được gia công tôi cứng bề mặt trên máy phay CNC 3 trục MCB-850, với các bán kính đỉnh chày chính xác $R_{ch} \in {10, 15, 20}\text{ mm}$ và góc mở cối $90^\circ$.
  • Dụng cụ đo biến dạng góc quang học và dưỡng kiểm bán kính số chuyên dụng.

Tam giác đạc (Triangulation) phương pháp luận được thiết lập thông qua việc đối chiếu chéo liên tục giữa: (1) Kết quả giải tích lý thuyết, (2) Dữ liệu mô phỏng số 3D trên ABAQUS, và (3) Số liệu thực nghiệm uốn trên máy kéo/nén.

Data và phân tích

Thiết kế thực nghiệm đa biến được tổ chức theo ma trận trực giao Taguchi $L_9(3^4)$ nhằm đánh giá tối ưu hóa đồng thời ba thông số công nghệ đầu vào:

  • Yếu tố A: Nhiệt độ nung phôi $T$ ($32^\circ\text{C}, 300^\circ\text{C}, 600^\circ\text{C}$)
  • Yếu tố B: Bán kính chày uốn $R_{ch}$ ($10\text{ mm}, 15\text{ mm}, 20\text{ mm}$)
  • Yếu tố C: Hành trình chày $H$ ($10\text{ mm}, 16\text{ mm}, 22\text{ mm}$)

Phân tích tỷ số tín hiệu trên nhiễu (Signal-to-Noise Ratio - $S/N$) được thiết lập theo tiêu chuẩn Càng nhỏ càng tốt (Smaller-the-Better - STB) cho lực uốn tạo hình $P$ và góc đàn hồi $\theta_{SP}$: $$S/N = -10 \log_{10} \left( \frac{1}{n} \sum_{i=1}^n y_i^2 \right)$$ Kiểm định thống kê phân tích phương sai (ANOVA) với tiêu chuẩn Fisher ($F$-test) ở độ tin cậy $95%$ ($p < 0.05$) được áp dụng để lượng hóa phần trăm đóng góp của từng thông số công nghệ đối với các đặc tính chất lượng sản phẩm.

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

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

Nghiên cứu mang lại 4 phát hiện thực nghiệm mang tính bước ngoặt, có độ tin cậy thống kê cao:

  1. Sự suy giảm vượt bậc của lực uốn tạo hình danh định ($P_{TN}$): Khi nâng nhiệt độ gia nhiệt từ nhiệt độ phòng ($32^\circ\text{C}$) lên $600^\circ\text{C}$, lực uốn tạo hình thực nghiệm giảm cực kỳ mạnh mẽ. Cụ thể, với cấu hình chày $R_{ch} = 10\text{ mm}$ và hành trình $H = 22\text{ mm}$, lực uốn đỉnh giảm từ $18.42\text{ kN}$ (ở $32^\circ\text{C}$) xuống $14.15\text{ kN}$ (ở $300^\circ\text{C}$) và chỉ còn $8.95\text{ kN}$ (ở $600^\circ\text{C}$), tương đương mức cắt giảm hơn $51.4%$ tải trọng uốn đòi hỏi lên kết cấu khung máy.
  2. Triệt tiêu hiện tượng đàn hồi lại sau biến dạng dẻo ($\theta_{SP}$): Góc đàn hồi lại $\theta_{SP}$ giảm tỷ lệ thuận với sự tăng của nhiệt độ gia nhiệt. Tại $R_{ch} = 15\text{ mm}, H = 16\text{ mm}$, góc đàn hồi thực nghiệm ở nhiệt độ phòng là $3^\circ 42'$, giảm xuống $2^\circ 15'$ ở $300^\circ\text{C}$, và ở $600^\circ\text{C}$ góc đàn hồi chỉ còn $0^\circ 58'$—tương ứng mức triệt tiêu biến dạng đàn hồi lên tới $73.8%$. Điều này giúp góc chi tiết sau uốn đạt độ chính xác gần như tuyệt đối so với biên dạng cối uốn.
  3. Thứ tự ưu tiên chi phối thông số công nghệ qua phân tích Taguchi/ANOVA: Phân tích tỷ số $S/N$ và phương sai ANOVA khẳng định nhiệt độ nung phôi ($T$) là yếu tố chi phối mạnh mẽ nhất, đóng góp hơn $70%$ vào sự thay đổi của cả lực tạo hình và góc đàn hồi lại. Hành trình của chày ($H$) đứng thứ hai (chi phối độ mở góc uốn $\theta$), trong khi bán kính chày uốn ($R_{ch}$) có ảnh hưởng lớn nhất đến bán kính trong thực tế ($R_{UC}$) của chi tiết sau dỡ tải.
  4. Tính ưu việt tuyệt đối của mô hình biến cứng kết hợp trong mô phỏng số: Khi sử dụng mô hình hóa bền kết hợp (combined hardening) được đề xuất trong luận án, sai số giữa lực tạo hình mô phỏng ($P_{mp}$) và thực nghiệm ($P_{TN}$) tại cả ba cấp nhiệt độ luôn được khống chế ở mức cực nhỏ (dao động từ $0.85%$ đến tối đa $4.92%$). Ngược lại, việc sử dụng các mô hình biến cứng đẳng hướng thuần túy dẫn đến sai số góc đàn hồi thực tế vượt quá $15 - 18%$.
Thông số công nghệ ($T, R_{ch}, H$) Lực uốn thực nghiệm $P_{TN}$ (kN) Lực uốn mô phỏng $P_{mp}$ (kN) Sai lệch lực $\Delta P$ (%) Góc đàn hồi thực nghiệm $\theta_{SP}$ Góc đàn hồi mô phỏng $\theta_{SP}$
$32^\circ\text{C}, 10\text{ mm}, 10\text{ mm}$ 6.85 6.71 -2.04% $4^\circ 12'$ $4^\circ 05'$
$32^\circ\text{C}, 10\text{ mm}, 22\text{ mm}$ 18.42 18.20 -1.19% $2^\circ 45'$ $2^\circ 38'$
$300^\circ\text{C}, 10\text{ mm}, 22\text{ mm}$ 14.15 13.88 -1.91% $1^\circ 40'$ $1^\circ 35'$
$600^\circ\text{C}, 10\text{ mm}, 22\text{ mm}$ 8.95 8.82 -1.45% $0^\circ 42'$ $0^\circ 40'$
$600^\circ\text{C}, 15\text{ mm}, 16\text{ mm}$ 7.62 7.45 -2.23% $0^\circ 58'$ $0^\circ 55'$
$600^\circ\text{C}, 20\text{ mm}, 22\text{ mm}$ 9.48 9.25 -2.42% $0^\circ 50'$ $0^\circ 48'$

Implications đa chiều

  • Về mặt học thuật: Cung cấp bộ cơ sở dữ liệu tham số ứng xử cơ học nhiệt của thép SS400 dày $5\text{ mm}$ chuẩn hóa cho phần mềm ABAQUS/CAE. Phương thức tổng hợp xác định tham số biến cứng kết hợp thông qua kiểm chuẩn mô phỏng mở ra một quy trình tính toán khả thi, áp dụng được cho nhiều mác thép kết cấu hàng hải khác (như AH32, AH36, EH36).
  • Về mặt công nghệ chế tạo: Cho phép các nhà máy đóng tàu (Nam Triệu, Phà Rừng, Bạch Đằng, Lilama Lisemco) thay thế các thiết bị ép chấn thủy lực công suất khổng lồ bằng hệ thống uốn cơ học công suất trung bình tích hợp đầu nung cảm ứng cao tần cục bộ, giảm tiêu hao năng lượng và khấu hao máy móc.
  • Về mặt kiểm soát chất lượng: Loại bỏ hoàn toàn công đoạn sửa nguội, nắn phẳng hay đốt nhiệt nắn chỉnh thủ công bằng ngọn lửa Oxy-Acetylene vốn gây biến tính tổ chức kim loại và suy giảm cơ tính cục bộ của kết cấu vỏ tàu.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, công trình vẫn tồn tại một số giới hạn nghiên cứu cần được nhìn nhận khách quan:

  1. Giả định đẳng nhiệt cục bộ: Nghiên cứu giả định nhiệt độ phôi phân bố đồng đều trong vùng nung uốn và không xét đến gradient giảm nhiệt do bức xạ và đối lưu không khí trong quá trình dỡ tải.
  2. Giới hạn vận tốc biến dạng: Thử nghiệm uốn được tiến hành ở dải vận tốc bán tĩnh (quasi-static), chưa xét đến ảnh hưởng của tốc độ biến dạng dẻo cao ($\dot{\varepsilon} > 10^0\text{ s}^{-1}$) vốn có thể gây ra hiện tượng nhớt dẻo (viscoplasticity).
  3. Phạm vi biên dạng hình học: Luận án mới chỉ tập trung vào uốn chữ V một chiều (1D curvature) trên tấm phẳng chữ nhật, chưa mở rộng cho các dạng tấm cong hai chiều phức tạp (double curvature plates) có mặt cầu hay mặt yên ngựa trên vỏ tàu.

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Mở rộng mô hình biến cứng kết hợp sang dải thép hợp kim đóng tàu cường độ cao (AH36, DH36) và các tấm siêu dày ($t \ge 20\text{ mm}$).
  • Xây dựng mô hình liên kết tương tác ba trường Nhiệt - Cơ - Lượng biến pha tổ chức tế vi (Thermo-Mechanical-Metallurgical - TMM coupling) để kiểm soát sự chuyển biến pha Austenit - Feralit trong quá trình làm nguội sau uốn.
  • Tích hợp cánh tay robot công nghiệp điều khiển linh hoạt đầu nung cảm ứng từ đa trục để uốn tạo hình các tấm vỏ tàu có biên dạng không gian 3D.

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

Luận án tạo ra những tác động lan tỏa mạnh mẽ cả về mặt học thuật lẫn thực tiễn sản xuất công nghiệp:

  • Tác động học thuật: Đóng góp một phương pháp luận chuẩn mực về tích hợp mô hình biến cứng phi tuyến vào phân tích phần tử hữu hạn quá trình gia công tạo hình tấm có gia nhiệt.
  • Chuyển đổi công nghệ đóng tàu: Việc ứng dụng gia nhiệt cảm ứng điện từ cục bộ giúp ngành đóng tàu tiết kiệm từ $40 - 50%$ năng lượng tiêu thụ so với các phương pháp nung lò truyền thống, giảm thiểu tối đa phế phẩm do sai lệch kích thước lắp ráp thân tàu.
  • Bảo vệ môi trường và an toàn lao động: Thay thế hoàn toàn phương pháp nung uốn thủ công bằng ngọn lửa Oxy-Acetylene độc hại và tiềm ẩn nguy cơ cháy nổ cao trong các xưởng đóng tàu bằng công nghệ gia nhiệt cảm ứng điện từ xanh, sạch, không phát thải khí độc hại.
┌──────────────────────────────────────────────────────────────────────────────────┐
│                             Broad Ecosystem Impacts                              │
└────────────────────────────────────────┬─────────────────────────────────────────┘
                                         │
 ┌───────────────────────────────────────┼───────────────────────────────────────┐
 │                                       │                                       │
 ▼                                       ▼                                       ▼
┌─────────────────────────┐   ┌─────────────────────────┐   ┌─────────────────────────┐
│     Academic Impact     │   │   Industry Efficiency   │   │ Environmental & Safety  │
│ - Accurate FEA model    │   │ - 51.4% forming load    │   │ - Replaces toxic Oxy-   │
│   for combined hardening│   │   reduction             │   │   Acetylene heating     │
│ - High-citation base in │   │ - Eliminates springback │   │ - Green, localized,     │
│   manufacturing & marine│   │   rework and scrap      │   │   energy-saving process │
└─────────────────────────┘   └─────────────────────────┘   └─────────────────────────┘

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương thức nhận dạng tham số mô hình biến cứng vật liệu kết hợp và quy trình mô phỏng phi tuyến nâng cao trên phần mềm ABAQUS/CAE.
  • Kỹ sư R&D và Thiết kế công nghệ tại các Nhà máy Đóng tàu: Có được bộ công thức toán học và bảng thông số công nghệ tối ưu hóa ($T, R_{ch}, H$) để ứng dụng trực tiếp vào việc lập quy trình uốn chế tạo bích tăng cứng hầm hàng, vách ngăn và cửa chống nước.
  • Các Doanh nghiệp chế tạo thiết bị cơ khí: Có cơ sở khoa học để thiết kế, chế tạo các cụm máy chấn tôn tích hợp modul nguồn nung cao tần nhỏ gọn, hiệu năng cao với chi phí đầu tư thấp hơn nhiều so với máy chấn siêu trường.

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ì và mở rộng mô hình nào?
Đóng góp độc đáo nhất là việc mở rộng quy luật biến cứng kết hợp Chaboche kết hợp quy luật đẳng hướng Voce sang môi trường nhiệt độ trung gian ($300^\circ\text{C}$ và $600^\circ\text{C}$) cho thép kết cấu SS400 dày $5\text{ mm}$. Luận án đã xây dựng quy trình tính toán ngược xác định hàm phụ thuộc nhiệt độ của tenxơ ứng suất ngược $\boldsymbol{\alpha}(T)$, khắc phục hoàn toàn sự thiếu chính xác của các mô hình đẳng hướng truyền thống trong việc dự đoán springback khi có sự đảo chiều ứng suất thớ uốn.

2. Đột phá về mặt phương pháp luận so với các nghiên cứu quốc tế trước đây?
So với nghiên cứu của Nguyen Duc-Toan et al. (2013) vốn chỉ dùng mô hình Johnson-Cook thuần túy cho hợp kim magie AZ31 hoặc Stachowicz et al. (2012) chỉ đánh giá thực nghiệm thép mỏng AMS 5604, luận án đã tích hợp thành công tam giác đạc phương pháp: Thử nghiệm kéo nhiệt trên máy HUNGTA H-200A $\rightarrow$ Nhận dạng số tham số biến cứng kết hợp trên ABAQUS $\rightarrow$ Tối ưu hóa đa mục tiêu bằng mảng trực giao Taguchi $L_9(3^4)$ trên chi tiết thực tế có độ dày lớn ($5\text{ mm}$).

3. Phát hiện thực nghiệm nào gây bất ngờ nhất về mặt dữ liệu?
Phát hiện bất ngờ nhất là tại ngưỡng nhiệt độ $600^\circ\text{C}$, góc biến dạng đàn hồi $\theta_{SP}$ bị triệt tiêu tới hơn $73%$, đưa góc uốn sau dỡ tải tiệm cận góc hình học của dụng cụ ($< 1^\circ$ sai lệch), trong khi lực uốn đỉnh giảm hơn một nửa ($51.4%$). Sự suy giảm đột ngột này chứng minh rằng tại $600^\circ\text{C}$, mật độ lệch mạng và cơ chế hồi phục nhiệt đã làm suy giảm gần như hoàn toàn năng lượng tích lũy đàn hồi trong cấu trúc tinh thể thép SS400.

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?
Có. Luận án cung cấp đầy đủ thông số kích thước hình học mẫu ($110 \times 35 \times 5\text{ mm}$ gia công cắt dây CNC dọc thớ cán), thông số kết cấu cuộn cảm ứng xoắn ốc bằng đồng đỏ làm mát nước, chế độ cài đặt bộ điều khiển nhiệt Autonics TZN4L, thông số phân lưới phần tử C3D8R và toàn bộ ma trận trực giao thực nghiệm Taguchi $L_9$.

5. Lộ trình nghiên cứu 10 năm tiếp theo từ kết quả luận án?
Lộ trình phát triển bao gồm: (1) Tích hợp trường vi cấu trúc kim loại học TMM vào phần mềm mô phỏng; (2) Mở rộng công nghệ uốn cảm ứng cho các mác thép hàng hải cường độ cao AH36, EH40 và hợp kim Titan; (3) Tự động hóa hệ thống uốn vỏ tàu 3D thời gian thực (Real-time Adaptive Forming) sử dụng cảm biến laser quét biên dạng kết hợp trí tuệ nhân tạo (AI) điều khiển nhiệt độ nguồn nung.

Kết luận

Công trình nghiên cứu của tác giả Vương Gia Hải đã giải quyết trọn vẹn một bài toán kỹ thuật công nghệ cấp thiết của ngành cơ khí đóng tàu với các đóng góp cốt lõi:

  1. Xác lập thành công hệ thống thử nghiệm thử kéo và uốn tấm SS400 dày $5\text{ mm}$ có tích hợp nguồn nhiệt cảm ứng điện từ cao tần cục bộ ở ba cấp nhiệt độ ($32^\circ\text{C}$, $300^\circ\text{C}$, $600^\circ\text{C}$).
  2. Đề xuất phương thức tổng hợp mới nhận dạng chính xác bộ tham số mô hình biến cứng kết hợp Voce-Chaboche phụ thuộc nhiệt độ tích hợp trong ABAQUS/CAE, kiểm soát sai số lực và góc uốn dưới $5%$.
  3. Lượng hóa quy luật tương tác giữa các thông số công nghệ ($T, R_{ch}, H$) thông qua quy hoạch thực nghiệm Taguchi $L_9$, chứng minh gia nhiệt $600^\circ\text{C}$ giúp giảm $51.4%$ lực uốn đỉnh và triệt tiêu $73.8%$ góc đàn hồi lại.
  4. Xây dựng 4 mô hình toán học hồi quy dự báo chính xác các thông số đầu ra phục vụ trực tiếp công tác thiết kế công nghệ và lựa chọn thiết bị gia công.
  5. Kiểm chứng thành công việc chế tạo các chi tiết bích tăng cứng hầm hàng và khung cửa tàu thủy thực tế, mở ra hướng ứng dụng công nghệ gia công gia nhiệt cảm ứng điện từ hiệu quả, bền vững cho ngành công nghiệp đóng tàu Việt Nam.