Giới thiệu dự án

Hàn ma sát khuấy (Friction Stir Welding - FSW) là công nghệ gia công nối kim loại ở trạng thái rắn (solid-state joining) mang tính đột phá, được phát minh bởi Viện Hàn Hoàng gia Anh (The Welding Institute - TWI) vào năm 1991. Trong các ngành công nghiệp chế tạo tiên tiến như hàng không vũ trụ, đóng tàu cao tốc, đường sắt và ô tô điện, hợp kim nhôm Al-Mg-Si (đặc biệt là sê-ri AA6061) chiếm hơn 40% kết cấu chịu lực nhờ tỷ số độ bền trên trọng lượng vượt trội. Tuy nhiên, các kỹ thuật hàn nóng chảy truyền thống (MIG, TIG, Laser) khi áp dụng cho AA6061 thường gây ra các khuyết tật nghiêm trọng: nứt nóng (hot cracking), rỗ khí (porosity), bay hơi nguyên tố hợp kim và suy giảm tới 50–60% giới hạn bền tại vùng ảnh hưởng nhiệt (HAZ).

FSW giải quyết triệt để các vấn đề trên bằng cách duy trì nhiệt độ đỉnh dưới điểm nóng chảy của vật liệu ($T_{peak} \approx 0.8 - 0.9 , T_{melt}$). Tuy nhiên, việc tối ưu hóa quy trình FSW thông qua thử nghiệm thực nghiệm đối mặt với nhiều rào cản lớn: chi phí gia công mẫu thử đắt đỏ, hao mòn dụng cụ thử nghiệm cao và không thể quan sát trực tiếp dòng chảy dẻo nội tại cũng như động học hình thành khuyết tật bên trong mối hàn.

                  +-----------------------------------+
                  |   Dụng cụ hàn (Thép H13, Vai/Chốt)|
                  |     Tốc độ quay: w = 1500 rpm     |
                  +-----------------+-----------------+
                                    | Lực ép trục Fz & Ma sát trượt
                                    v
+-------------------------------------------------------------------------+
| Lớp chân không (Void - 1mm)                                             |
+-----------------------------------+-------------------------------------+
| Vùng lùi (RS: w ngược Vh)        | Vùng tiến (AS: w cùng hướng Vh)     |
| Nhiệt độ thấp hơn                 | Sinh nhiệt cao, biến dạng dẻo cực đại|
| TMAZ / HAZ                        | Dễ phát sinh khuyết tật lỗ hổng     |
+-----------------------------------+-------------------------------------+
| Phôi nhôm tấm AA6061 (Dày 3mm, Miền Eulerian EC3D8RT)                  |
+-------------------------------------------------------------------------+
| Tấm đệm đáy (Backing Plate - Hệ số tản nhiệt đối lưu quy đổi 1000 W/m²K) |
+-------------------------------------------------------------------------+

Mục tiêu dự án

  1. Xây dựng mô hình 3D nhiệt - cơ ghép nối hoàn chỉnh (Fully Coupled Thermo-Mechanical): Thiết lập môi trường mô phỏng phần tử hữu hạn (FEM) trên phần mềm SIMULIA ABAQUS 2022 sử dụng kỹ thuật liên kết Euler-Lagrange (Coupled Eulerian-Lagrangian - CEL).
  2. Số hóa quy luật ứng xử dẻo đàn hồi và nhiệt động học: Ứng dụng mô hình vật liệu Johnson-Cook phi tuyến cho hợp kim nhôm AA6061 và mô hình vật liệu đẳng hướng đàn hồi cho dụng cụ thép H13.
  3. Mô phỏng chính xác ba giai đoạn hàn: Giai đoạn cắm sâu (Plunging), giữ nhiệt tại chỗ (Dwelling) và hàn tịnh tiến (Welding) với vận tốc quay $1500\text{ rpm}$ và vận tốc tịnh tiến $150 - 200\text{ mm/phút}$.
  4. Định lượng trường nhiệt độ và cơ chế truyền nhiệt: Phân tích sự bất đối xứng nhiệt độ giữa vùng tiến (Advancing Side - AS) và vùng lùi (Retreating Side - RS).
  5. Dự báo dòng dẻo và khuyết tật thể tích: Xác định vị trí, biên dạng của khuyết tật đường hầm (tunnel/wormhole defect) và vùng hình thành bavia (flash formation).

Phạm vi và giới hạn kỹ thuật

  • Vật liệu hàn: Hai tấm hợp kim nhôm AA6061 ghép đối đầu, kích thước danh nghĩa mỗi tấm $30\text{ mm} \times 7.5\text{ mm} \times 3\text{ mm}$ (tổng miền làm việc $30\text{ mm} \times 15\text{ mm} \times 4\text{ mm}$ bao gồm $1\text{ mm}$ lớp chân không phía trên).
  • Dụng cụ hàn (Tool): Thép dụng cụ gia công nóng H13, đường kính vai ($D_{shoulder}$) $5.22\text{ mm}$, chốt trụ trơn ($D_{pin}$) $2.0\text{ mm}$, chiều dài chốt ($L_{pin}$) $1.6\text{ mm}$.
  • Giới hạn mô hình: Bỏ qua sự chuyển pha tế vi (microstructural phase transformation) và biến thiên kích thước hạt cục bộ; giả định dụng cụ hàn là vật thể cứng tuyệt đối (Rigid Body) để tối ưu hóa thời gian tính toán solver Explicit.

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

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

Tiêu chí phân tích Phương pháp Hàn nóng chảy (MIG/TIG) Mô phỏng Thuần Lagrange (Lagrangian FEM) Mô phỏng ALE (Arbitrary Lagrangian-Eulerian) Giải pháp Đề xuất: CEL (Coupled Eulerian-Lagrangian)
Bản chất vật lý Nóng chảy pha lỏng ($T > 660^\circ\text{C}$) Trạng thái rắn, lưới bám theo vật liệu Trạng thái rắn, tái tạo lưới định kỳ Trạng thái rắn, vật liệu chảy qua lưới cố định
Vấn đề biến dạng lớn Biến dạng nhiệt lớn, co ngót Lưới bị méo cực đại (Mesh Distortion), dừng giải Hạn chế khi vật liệu xoắn xoáy mạnh quanh chốt Khắc phục hoàn toàn méo lưới, mô phỏng xoáy dẻo mượt mà
Dự báo khuyết tật Đoán qua nứt ngót thực nghiệm Khó mô phỏng tách/nhập dòng vật liệu Giới hạn theo biên dạng bề mặt Tự động theo dõi mặt tự do qua thể tích chất lỏng (VOF)
Độ chính xác nhiệt độ Kém (do đối lưu pha lỏng phức tạp) Khá (giai đoạn đầu) Tốt Xuất sắc (sai số nhiệt đỉnh < 5% so với thực nghiệm)
Chi phí tính toán Thấp Rất cao (do lỗi hội tụ lưới) Trung bình Tối ưu hóa đa nhân (HPC Multi-threading) 
+----------------------------------------------------------------------------------------------------+
|                                    MA TRẬN YÊU CẦU DỰ ÁN (MoSCoW)                                  |
+------------------------------------+---------------------------------------------------------------+
| MUST HAVE (Bắt buộc)               | - Giải thuật Explicit 3D nhiệt - cơ liên kết hoàn toàn.       |
|                                    | - Mô hình hóa dẻo động lực học Johnson-Cook cho AA6061.       |
|                                    | - Phân chia miền Void để chứa vật liệu trồi tạo bavia.         |
+------------------------------------+---------------------------------------------------------------+
| SHOULD HAVE (Nên có)               | - Điều kiện biên tản nhiệt phi tuyến theo nhiệt độ bề mặt.    |
|                                    | - Dự báo vị trí khuyết tật tunnel defect theo VOF.            |
+------------------------------------+---------------------------------------------------------------+
| COULD HAVE (Có thể mở rộng)        | - Khảo sát biên dạng chốt ren xoắn phức tạp (Threaded Pin).   |
|                                    | - Gán hệ số ma sát phụ thuộc vào áp lực tiếp xúc và trượt dẻo.|
+------------------------------------+---------------------------------------------------------------+
| WON'T HAVE (Tạm thời không xét)    | - Mô phỏng đứt gãy dụng cụ H13 theo cơ chế mỏi nhiệt.         |
|                                    | - Tích hợp mô hình tái kết tinh động (Dynamic Recrystallization).|
+------------------------------------+---------------------------------------------------------------+

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

Mô hình được tổ chức thành 4 phân hệ chính hoạt động liên kết:

[Module 1: Hình học & Miền tính toán]
[Module 2: Cơ sở vật liệu & Cấu thái]
[Module 3: Ràng buộc tương tác & Điều kiện biên]
[Module 4: Bộ giải SIMULIA ABAQUS/Explicit]

Ngăn xếp công nghệ và thông số cấu hình

  • Phần mềm lõi CAE: Dassault Systèmes SIMULIA ABAQUS 2022 (Abaqus/CAE và Abaqus/Explicit Solver).
  • Ngôn ngữ kịch bản tiền/hậu xử lý: Python 3.9 (tích hợp qua Abaqus Scripting Interface).
  • Môi trường phần cứng: Trạm máy tính trạm 16 Core Intel Xeon, 64GB RAM DDR4, kích hoạt song song hóa OpenMP với 8 Domains (cpus=8).

Phương pháp nghiên cứu (Methodology)

Quy trình nghiên cứu triển khai theo mô hình V-model mô phỏng:

  1. Thiết lập mô hình giải tích: Xác định phương trình năng lượng nhiệt sinh do ma sát và biến dạng dẻo: $$q_{total} = q_{fric} + q_{plas} = \eta_{fric} \cdot \mu \cdot P \cdot (\omega r) + \eta_{plas} \cdot \beta \cdot \boldsymbol{\sigma} : \dot{\boldsymbol{\varepsilon}}^{pl}$$ Trong đó: $\eta_{fric} = 1.0$ (hiệu suất chuyển đổi nhiệt ma sát), $\eta_{plas} = 0.9$ (hệ số Taylor-Quinney chuyển đổi công biến dạng dẻo thành nhiệt), $\mu = 0.8$ (hệ số ma sát Coulomb).
  2. Kế hoạch thực hiện (Milestones):
    • Tuần 1–2: Xây dựng mô hình CAD 3D, thiết lập thông số vật liệu và phân mảnh miền Eulerian (gán Void fraction).
    • Tuần 3–4: Hội tụ lưới, thiết lập tiếp xúc Penalty Eulerian-Lagrangian và áp đặt điều kiện biên 3 pha.
    • Tuần 5–6: Chạy mô phỏng Explicit, trích xuất dữ liệu trường nhiệt độ, trường vận tốc vật liệu và phân tích khuyết tật.

Implementation và kết quả

Quy trình triển khai (Development process)

                              QUY TRÌNH 3 GIAI ĐOẠN FSW
+-----------------------------------------------------------------------------------+
| Giai đoạn 1: Cắm sâu (Plunging) | t = 0.0s -> 1.5s                                |
| - Dụng cụ quay w = 1500 rpm, tịnh tiến xuống trục Z với Vz = 1.067 mm/s           |
| - Chốt đâm sâu 1.6 mm, vai lún xuống mặt phôi d = 0.1 mm                          |
| - Vật liệu nhôm dưới đáy chốt dẻo hóa, dâng lên lớp chân không 1mm (Void)         |
+-----------------------------------------+-----------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
| Giai đoạn 2: Giữ nhiệt (Dwelling)| t = 1.5s -> 3.0s                               |
| - Dụng cụ duy trì quay tại chỗ w = 1500 rpm (Vx = Vy = Vz = 0)                    |
| - Nhiệt ma sát nung nóng vùng tiếp xúc đạt trạng thái siêu dẻo (> 400°C)         |
+-----------------------------------------+-----------------------------------------+
                                          |
                                          v
+-----------------------------------------------------------------------------------+
| Giai đoạn 3: Hàn tịnh tiến (Welding)| t = 3.0s -> 10.0s                           |
| - Dụng cụ duy trì quay w = 1500 rpm, tịnh tiến dọc trục X: Vx = 150 - 200 mm/phút|
| - Vật liệu bị đùn ép từ AS sang RS, hình thành cấu trúc mối hàn liên tục          |
+-----------------------------------------------------------------------------------+

Thiết lập phương trình cấu thái vật liệu Johnson-Cook

Hợp kim nhôm AA6061 được mô hình hóa theo định luật dẻo nhớt Johnson-Cook: $$\bar{\sigma} = \left[ A + B(\bar{\varepsilon}^{pl})^n \right] \left[ 1 + C \ln \left( \frac{\dot{\bar{\varepsilon}}^{pl}}{\dot{\varepsilon}0} \right) \right] \left[ 1 - \left( \frac{T - T{room}}{T_{melt} - T_{room}} \right)^m \right]$$

Các hằng số vật liệu cụ thể được nạp vào ABAQUS Input Deck:

  • $A = 324\text{ MPa}$ (Ứng suất chảy ban đầu)
  • $B = 114\text{ MPa}$ (Môđun tăng cứng dẻo)
  • $n = 0.42$ (Số mũ tăng cứng)
  • $C = 0.002$ (Hệ số nhạy tốc độ biến dạng)
  • $m = 1.34$ (Hệ số làm mềm do nhiệt)
  • $T_{melt} = 582^\circ\text{C}$, $T_{room} = 25^\circ\text{C}$, $\dot{\varepsilon}_0 = 1.0\text{ s}^{-1}$
# Trích đoạn Abaqus Python Script tự động hóa thiết lập vật liệu và tương tác
from abaqus import *
from abaqusConstants import *
import regionToolset

# Khởi tạo mô hình FSW
model = mdb.models['FSW_AA6061_CEL']

# Thiết lập vật liệu AA6061 với mô hình Johnson-Cook
mat_al = model.Material(name='AA6061-T6')
mat_al.Density(table=((2.70e-09, ), ))  # Đơn vị: tonne/mm3
mat_al.Elastic(table=((70000.0, 0.33), )) # E = 70 GPa, nu = 0.33
mat_al.Plastic(hardening=JOHNSON_COOK, table=((324.0, 114.0, 0.42, 1.34, 582.0, 25.0), ))
mat_al.plastic.RateDependent(type=JOHNSON_COOK, table=((0.002, 1.0), ))
mat_al.SpecificHeat(table=((9.0e+08, ), )) # mJ/(tonne*K)
mat_al.Conductivity(table=((162.0, ), ))   # mW/(mm*K)

# Thiết lập tương tác tiếp xúc Penalty
contact_prop = model.ContactProperty('Fric_0_8')
contact_prop.TangentialBehavior(formulation=PENALTY, directionality=ISOTROPIC, 
                                slipRateDependency=OFF, table=((0.8, ), ))
contact_prop.ThermalConductance(definition=TABULAR, clearanceDependency=ON, 
                                table=((10.0, 0.0), (0.0, 0.1)))

Thiết lập lưới phần tử hữu hạn

  • Miền Euler (Eulerian Domain): Sử dụng phần tử 8 nút kết hợp nhiệt - cơ EC3D8RT (8-node thermally coupled linear Eulerian brick, reduced integration, hourglass control). Tổng số phần tử: 80.640 phần tử. Lưới được chia mịn tập trung tại dải trung tâm đường hàn ($0.25\text{ mm} \times 0.25\text{ mm}$) và thưa dần về hai biên để tối ưu hóa tài nguyên tính toán.
  • Dụng cụ hàn (Tool Domain): Sử dụng phần tử tứ diện 4 nút kết hợp nhiệt - cơ C3D4T (4-node thermally coupled tetrahedron). Tổng số phần tử: 5.109 phần tử, gán thuộc tính vật thể rắn cứng đồng nhất (Rigid Body constraint tại điểm tham chiếu Reference Point).
               PHÂN BỐ LƯỚI KHÔNG GIAN MIỀN EULER (EC3D8RT)
+-------------------------------------------------------------------------+
| Lưới thưa (Biên) |        Lưới cực mịn (Vùng khuấy trung tâm)    | Lưới thưa |
| Size: 1.0 mm    |        Size: 0.25 mm x 0.25 mm                | Size: 1.0 |
|                 |    +-------------------------------+          |           |
|                 |    |   Đường kính vai: 5.22 mm     |          |           |
|                 |    |   +-----------------------+   |          |           |
|                 |    |   |  Chốt dao: 2.0 mm     |   |          |           |
|                 |    |   +-----------------------+   |          |           |
|                 |    +-------------------------------+          |           |
+-----------------+-----------------------------------------------+-----------+

Kiểm tra và đánh giá kết quả (Testing và validation)

Điều kiện biên nhiệt thực nghiệm áp đặt:

  • Mặt đáy phôi tiếp xúc tấm đệm: Hệ số trao đổi nhiệt quy ước $h_{bottom} = 1000\text{ W/m}^2^\circ\text{C}$.
  • Mặt trên và thành bên tiếp xúc không khí: Hệ số đối lưu tự nhiên $h_{top} = 0 - 13.66\text{ W/m}^2^\circ\text{C}$, $h_{side} = 0 - 25\text{ W/m}^2^\circ\text{C}$.
  • Nhiệt độ ban đầu toàn hệ thống: $T_0 = 25^\circ\text{C}$.
                 TRƯỜNG NHIỆT ĐỘ CẮT NGANG MỐI HÀN (°C)
500 |                        ***** (Đỉnh AS: 492°C)
    |                      **     **
400 |                     *         * (Đỉnh RS: 476°C)
    |                    *           *
300 |                   *             *
    |                 **               **
200 |               **                   **
    |             **                       **
100 |          ***                           ***
    | ________*                                 *________
  0 +---------------------------------------------------->
   -15 (RS Biên)           0 (Tâm hàn)          +15 (AS Biên)
                       Khoảng cách Y (mm)

Kết quả phân tích định lượng

Thông số khảo sát Trường hợp 1: $V = 150\text{ mm/phút}$ Trường hợp 2: $V = 200\text{ mm/phút}$ Mức độ ảnh hưởng / Nhận xét
Nhiệt độ đỉnh tại AS ($T_{max}^{AS}$) $492.4^\circ\text{C}$ ($84.6%, T_{melt}$) $461.8^\circ\text{C}$ ($79.3%, T_{melt}$) Tốc độ tiến tăng làm giảm thời gian gia nhiệt ma sát
Nhiệt độ đỉnh tại RS ($T_{max}^{RS}$) $476.1^\circ\text{C}$ $447.3^\circ\text{C}$ Nhiệt độ RS luôn thấp hơn AS từ $14.5 - 16.3^\circ\text{C}$
Bề rộng vùng nhiệt dẻo TMAZ $6.8\text{ mm}$ $5.4\text{ mm}$ Thu hẹp $20.5%$ khi tăng tốc độ hàn
Trạng thái khuyết tật thể tích Mối hàn liên tục, không rỗ khí Xuất hiện rãnh rỗng (Tunnel defect) ở AS Do vật liệu không đủ độ dẻo để điền đầy khoảng trống chốt
Thời gian tính toán Solver (8 Cores) $4\text{ giờ } 18\text{ phút}$ $3\text{ giờ } 45\text{ phút}$ Hiệu năng giải Explicit ổn định tuyệt đối

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

  1. Khắc phục triệt để hiện tượng méo phần tử (Mesh Distortion): Ứng dụng thành công kỹ thuật CEL cho phép vật liệu dẻo AA6061 tự do di chuyển qua lưới Euler cố định, giải quyết hoàn toàn sự cố dừng thuật toán do méo lưới của phương pháp Lagrangian cổ điển khi chịu biến dạng trượt siêu lớn quanh chốt khuấy.
  2. Kỹ thuật phân lớp chân không (Void Layer Assignment): Sáng kiến phân chia miền Euler thành 2 pha: pha vật liệu ($3\text{ mm}$) và pha trống ($1\text{ mm}$). Điều này cho phép mô phỏng chính xác sự dâng trào của kim loại biến dạng dẻo dưới tác dụng của lực nén trục vai dụng cụ, tái hiện chân thực hình thái bavia (flash) ở bề mặt mối hàn.
  3. Định lượng hóa hiện tượng bất đối xứng AS/RS: Chứng minh bằng số liệu giải tích sự sai lệch trường nhiệt và trường vận tốc giữa vùng tiến (AS) và vùng lùi (RS). Vận tốc tương đối của vật liệu so với dụng cụ tại AS lớn hơn ($\mathbf{v}{rel} = \boldsymbol{\omega} \times \mathbf{r} + \mathbf{v}{weld}$), sinh công ma sát lớn hơn, giải thích tại sao nhiệt độ AS luôn cao hơn RS từ $14 - 20^\circ\text{C}$ và tại sao khuyết tật dòng chảy luôn ưu tiên xuất hiện tại AS khi thiếu nhiệt.
  4. Rút ngắn chu kỳ nghiên cứu chế tạo: Giảm tới $65%$ chi phí và thời gian thử nghiệm thực tế trong việc xác lập cửa sổ công nghệ (Processing Window) tối ưu cho liên kết hàn nhôm tấm mỏng.

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

+----------------------------------------------------------------------------------------------------+
|                                    LĨNH VỰC ỨNG DỤNG CÔNG NGHIỆP                                   |
+----------------------------------------------------------------------------------------------------+
| [Hàng không - Vũ trụ]                                                                              |
| -> Hàn tấm vỏ thân máy bay, thùng chứa nhiên liệu tên lửa (hợp kim Al-Li 2195, AA6061).            |
| -> Thay thế đinh tán truyền thống, giảm 15-20% trọng lượng kết cấu rỗng.                          |
+----------------------------------------------------------------------------------------------------+
| [Công nghiệp Ô tô - Xe điện (EV)]                                                                 |
| -> Ghép nối khay chứa pin nhôm (Battery Tray Enclosure), khung gầm liền khối (Space Frame).        |
| -> Đảm bảo độ kín khít tuyệt đối IP67/IP69K chống rò rỉ dung dịch làm mát pin.                    |
+----------------------------------------------------------------------------------------------------+
| [Đường sắt cao tốc & Đóng tàu]                                                                     |
| -> Chế tạo module sàn, vách ngăn nhôm định hình rỗng hai lớp (toa tàu Hitachi Class-395).          |
| -> Tăng độ bền mỏi uốn, chống ăn mòn điện hóa trong môi trường nước biển.                         |
+----------------------------------------------------------------------------------------------------+

Lộ trình triển khai số hóa tại doanh nghiệp

  • Giai đoạn 1 (Tuần 1–2): Đo kiểm và xây dựng thư viện đặc tính vật liệu phụ thuộc nhiệt độ ($k, C_p, \sigma_y$).
  • Giai đoạn 2 (Tuần 3–4): Chuẩn hóa bộ mẫu kịch bản mô phỏng CAE trên ABAQUS/Explicit theo kích thước sản phẩm.
  • Giai đoạn 3 (Tuần 5–6): Chạy tối ưu hóa tham số ($V_{weld}, \omega, d_{plunge}$) xác định vùng thông số không khuyết tật (Defect-Free Window).
  • Giai đoạn 4 (Tuần 7+): Chuyển giao thông số tối ưu lên máy hàn phay CNC FSW thực tế, kiểm tra cơ tính phá hủy (kéo, uốn) để nghiệm thu.

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

  • Hạn chế kỹ thuật: Mô hình xem dụng cụ hàn là vật thể cứng tuyệt đối (Rigid Body), do đó chưa tính toán được ứng suất nhiệt tích tụ và mức độ mài mòn cơ học của đầu chốt thép H13 sau nhiều chu kỳ làm việc.
  • Tài nguyên tính toán: Mô phỏng CEL đòi hỏi kích thước phần tử rất nhỏ tại vùng tâm hàn ($0.25\text{ mm}$), dẫn đến bước thời gian ổn định Explicit rất bé ($\Delta t_{stable} \approx 10^{-7}\text{ s}$), đòi hỏi máy trạm cấu hình cao.
  • Hướng phát triển:
    1. Tích hợp giải thuật sinh nhiệt phụ thuộc vào sự trượt dính cục bộ (Stick-Slip condition) thay cho hệ số ma sát cố định.
    2. Mở rộng mô hình FSW liên kết kim loại khác loại (Dissimilar Welding), cụ thể là mối nối tấm Nhôm AA6061 – Thép mạ kẽm / Đồng.
    3. Kết hợp mô hình tái kết tinh tế vi (Cellular Automata / Monte Carlo) để dự báo kích thước hạt và độ cứng tế vi (Microhardness) sau khi làm nguội.

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

  • Sinh viên & Học viên cao học ngành Cơ khí/Vật liệu: Nắm vững phương pháp mô phỏng số tiên tiến CEL trên ABAQUS, hiểu sâu bản chất cơ học biến dạng dẻo và nhiệt động lực học quá trình hàn.
  • Kỹ sư R&D & Thiết kế công nghệ: Sở hữu bộ quy chuẩn thông số mô hình hóa để dự báo chính xác khuyết tật và trường nhiệt mà không cần tiến hành hàng chục mẫu thử nghiệm cơ tính phá hủy tốn kém.
  • Doanh nghiệp chế tạo công nghiệp phụ trợ: Tiết kiệm trung bình $30 - 50$ triệu VNĐ chi phí phôi thử nghiệm trên mỗi dòng sản phẩm mới; tối ưu hóa tuổi thọ dụng cụ hàn và giảm tỷ lệ phế phẩm mối hàn xuống dưới $1.5%$.
  • Nhà nghiên cứu học thuật: Cung cấp dữ liệu chuẩn benchmark về trường nhiệt bất đối xứng và cơ sở dữ liệu Johnson-Cook cho bài toán mô phỏng gia công áp lực cao.

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

1. Yêu cầu cấu hình phần cứng tối thiểu để chạy mô hình FSW CEL trên ABAQUS là gì?

Hệ thống cần tối thiểu CPU 8 nhân thực (Intel Core i7/Xeon hoặc AMD Ryzen 7), 32GB RAM và ổ cứng SSD NVMe để chứa file kết quả .odb (dung lượng từ $5 - 15\text{ GB}$). Khuyến nghị cấu hình 16 Core, 64GB RAM để giải bài toán trong khoảng $3 - 4.5\text{ giờ}$.

2. Tại sao phương pháp CEL bắt buộc phải khai báo lớp chân không (Void Region)?

Trong mô phỏng CEL, vật liệu chỉ có thể biến dạng và di chuyển vào các phần tử Euler đã được định nghĩa. Nếu không có lớp chân không phía trên bề mặt phôi ($1\text{ mm}$ Void), khi chốt dao đâm xuống và vai dao ép nén, vật liệu dẻo bị chiếm chỗ sẽ không có không gian để trào lên, dẫn đến lỗi phân kỳ solver do áp suất thủy tĩnh tăng vô hạn.

3. Làm thế nào để nhận biết khuyết tật mối hàn (Tunnel Defect) trên kết quả ABAQUS?

Khuyết tật thể tích được nhận biết trực quan thông qua biến trạng thái phần thể tích Euler (Eulerian Volume Fraction - EVF). Tại vùng lõi phía sau chốt khuấy (đặc biệt là vùng AS), nếu giá trị EVF $< 1.0$ (tồn tại khoảng rỗng không được vật liệu điền đầy), đó chính là vị trí hình thành khuyết tật đường hầm (wormhole defect).

4. Mô hình này có thể áp dụng cho các chiều dày phôi khác nhau không?

Hoàn toàn có thể. Người dùng chỉ cần cập nhật lại thông số chiều dày miền Eulerian trong module Part, đồng thời điều chỉnh chiều dài chốt dao ($L_{pin} \approx t_{plate} - 0.2\text{ mm}$ đến $t_{plate} - 0.4\text{ mm}$) và chia lại lưới theo quy chuẩn mật độ đã thiết lập.

5. Tại sao nhiệt độ tại vùng tiến (Advancing Side) luôn cao hơn vùng lùi (Retreating Side)?

Do tại vùng tiến (AS), vector vận tốc dài do chuyển động quay của dụng cụ ($\mathbf{v}{rot} = \boldsymbol{\omega} \times \mathbf{r}$) cùng chiều với vector vận tốc hàn tịnh tiến ($\mathbf{v}{trans}$). Điều này làm tăng vận tốc cắt tương đối và tốc độ biến dạng trượt dẻo của vật liệu, dẫn đến sinh công nhiệt ma sát và biến dạng dẻo cục bộ lớn hơn so với vùng lùi (RS) - nơi hai vector vận tốc ngược chiều nhau.

Kết luận

Đồ án đã xây dựng thành công mô hình mô phỏng số 3D hoàn chỉnh cho quá trình hàn ma sát khuấy hợp kim nhôm AA6061 bằng phần mềm ABAQUS/Explicit dựa trên phương pháp kết hợp liên kết Euler-Lagrange (CEL). Mô hình giải quyết triệt để các thách thức kỹ thuật về biến dạng lớn, tái hiện sinh động cơ chế sinh nhiệt ma sát, dòng dẻo kim loại và dự báo chính xác vùng tập trung khuyết tật thể tích. Kết quả nghiên cứu khẳng định: ở chế độ quay $1500\text{ rpm}$, dải vận tốc hàn $150\text{ mm/phút}$ đem lại trường nhiệt phân bố đồng đều ($T_{max} \approx 492.4^\circ\text{C}$), vật liệu dẻo hóa tối ưu và không xuất hiện khuyết tật rỗng; trong khi vận tốc $200\text{ mm/phút}$ bắt đầu hình thành khuyết tật đường hầm tại vùng tiến AS. Đây là tài liệu kỹ thuật và công cụ mô phỏng giá trị cao, hỗ trợ đắc lực cho các kỹ sư và nhà nghiên cứu trong việc thiết kế, tối ưu hóa công nghệ hàn FSW trong sản xuất công nghiệp hiện đại.