Sự khác biệt căn bản giữa lý thuyết nhiệt đàn hồi và lý thuyết nhiệt đàn dẻo là gì?

Lý thuyết nhiệt đàn hồi chỉ khảo sát ứng xử cơ học trong giới hạn đàn hồi, nơi vật thể có thể hồi phục nguyên vẹn trạng thái hình học ban đầu sau khi tải trọng và nhiệt độ được đưa về mức tham chiếu. Ngược lại, lý thuyết nhiệt đàn dẻo khảo sát toàn diện cả giai đoạn vượt qua giới hạn chảy, nơi các biến dạng dẻo vĩnh viễn không thể đảo ngược bắt đầu tích lũy, đi kèm với sự suy giảm cơ tính vật liệu theo nhiệt độ và hiện tượng tiêu tán năng lượng sinh nhiệt nội tại.

Tại sao quá trình nung nóng và làm nguội không đồng đều trong hàn lại phát sinh ứng suất dư?

Khi hàn, vùng kim loại bị nung nóng cục bộ có xu hướng giãn nở nhiệt nhưng bị kiềm chế bởi vùng kim loại lạnh xung quanh, làm phát sinh ứng suất nén vượt quá giới hạn chảy và gây ra biến dạng dẻo nén. Trong quá trình làm nguội tiếp theo, kim loại co ngót lại nhưng do đã bị biến dạng dẻo trước đó, sự co ngót này bị cản trở ngược lại bởi các vùng lân cận, dẫn đến sự tồn đọng của một trường ứng suất dư kéo có giá trị xấp xỉ giới hạn chảy của vật liệu.

Mô hình nguồn nhiệt Goldak có vai trò như thế nào trong phân tích nhiệt đàn dẻo?

Mô hình nguồn nhiệt khối elip kép của Goldak cho phép mô tả chính xác sự phân bố mật độ thông lượng nhiệt ba chiều của vũng hàn di động, tái hiện chính xác gradient nhiệt độ dốc ở phía trước và thoải ở phía sau mỏ hàn. Đây là thông số đầu vào quyết định tính đúng đắn của trường nhiệt độ phi tuyến trước khi chuyển sang giải bài toán ứng suất biến dạng nhiệt đàn dẻo.

Nhiệt độ giữa các lượt hàn (interpass temperature) ảnh hưởng thế nào đến biến dạng hàn?

Nhiệt độ giữa các lượt hàn quyết định trạng thái nhiệt ban đầu và tốc độ làm nguội của mỗi đường hàn tiếp theo. Nếu nhiệt độ giữa các lượt hàn không được kiểm soát chặt chẽ, sự tích lũy nhiệt độ nền quá cao sẽ làm mở rộng vùng chịu ảnh hưởng nhiệt, giảm giới hạn chảy của kim loại cơ bản và gia tăng mức độ biến dạng cong vênh của toàn bộ kết cấu sau khi nguội hoàn toàn.

Nhiệt đàn dẻo là gì? Cơ sở lý thuyết và mô phỏng số

Nhiệt đàn dẻo (tiếng Anh: thermoelastoplasticity hoặc thermo-elasto-plasticity) là một phân ngành nền tảng của cơ học môi trường liên tục và cơ học vật rắn biến dạng, nghiên cứu quy luật ứng xử cơ học của vật liệu khi chịu đồng thời tác động của tải trọng cơ học bên ngoài và sự biến thiên của trường nhiệt độ không gian thời gian. Trong lý thuyết nhiệt đàn dẻo, quá trình biến dạng của vật thể không chỉ bao gồm các phản ứng đàn hồi tức thời có khả năng hồi phục hoàn toàn, mà còn tích lũy các biến dạng dẻo vĩnh viễn không thể đảo ngược do trượt mạng tinh thể hoặc phá hủy vi mô, kết hợp chặt chẽ với hiện tượng giãn nở hoặc co ngót vì nhiệt. Sự tương tác hai chiều giữa trường nhiệt độ và trường ứng suất biến dạng tạo nên một hệ thống các phương trình vi phân đạo hàm riêng phi tuyến phức tạp, đóng vai trò quyết định trong việc phân tích độ bền, tính toàn vẹn kết cấu và dự đoán ứng suất dư trong các quá trình gia công nhiệt cơ kỹ thuật cao.

Sự hình thành và phát triển của cơ học nhiệt đàn dẻo bắt nguồn từ nhu cầu cấp thiết của các ngành công nghiệp chế tạo máy, hàng không vũ trụ, kỹ thuật hạt nhân và luyện kim, nơi các chi tiết kết cấu phải làm việc trong điều kiện khắc nghiệt với gradient nhiệt độ cực lớn và trạng thái ứng suất đa trục phức tạp. Khác với lý thuyết đàn hồi nhiệt chỉ xét biến dạng trong giới hạn đàn hồi, lý thuyết nhiệt đàn dẻo mở rộng phạm vi khảo sát sang vùng biến dạng dẻo có kể đến sự thay đổi rõ rệt của các hằng số vật liệu theo nhiệt độ, từ đó cung cấp cơ sở khoa học để tối ưu hóa thiết kế và kiểm soát độ an toàn của công trình.

Cơ sở nhiệt động lực học và các phương trình cấu phân

Để xây dựng mô hình toán học nhất quán cho vật liệu nhiệt đàn dẻo, lý thuyết cơ học môi trường liên tục hiện đại dựa trên hai nguyên lý cơ bản của nhiệt động lực học các quá trình không thuận nghịch. Trạng thái nhiệt động của một điểm vật chất được xác định thông qua tập hợp các biến trạng thái quan sát được và các biến nội bộ mô tả cơ chế trượt dẻo vi cấu trúc.

Phân rã tenxơ biến dạng ở phạm vi biến dạng bé

Trong khuôn khổ giả thiết biến dạng bé, tenxơ biến dạng tổng cộng $\boldsymbol{\varepsilon}$ được phân rã cộng tính thành ba thành phần độc lập:

\boldsymbol{\varepsilon} = \boldsymbol{\varepsilon}^e + \boldsymbol{\varepsilon}^p + \boldsymbol{\varepsilon}^{th}

Trong đó:

$\boldsymbol{\varepsilon}^e$ là tenxơ biến dạng đàn hồi có khả năng hồi phục khi dỡ bỏ tải trọng, liên kết trực tiếp với tenxơ ứng suất Cauchy thông qua định luật Hooke mở rộng có mô đun đàn hồi phụ thuộc nhiệt độ.
$\boldsymbol{\varepsilon}^p$ là tenxơ biến dạng dẻo vĩnh viễn, đặc trưng cho sự tích lũy trượt không thuận nghịch của cấu trúc vật liệu khi ứng suất vượt quá ngưỡng chảy dẻo.
$\boldsymbol{\varepsilon}^{th}$ là tenxơ biến dạng do giãn nở nhiệt thuần túy, được xác định thông qua biểu thức:
$\boldsymbol{\varepsilon}^{th} = \alpha (T - T_0) \mathbf{I}$
với $\alpha$ biểu thị hệ số giãn nở nhiệt tuyến tính của vật liệu, $T$ là nhiệt độ tức thời tại điểm khảo sát, $T_0$ là nhiệt độ chuẩn tham chiếu ban đầu, và $\mathbf{I}$ là tenxơ đơn vị bậc hai.

Phân rã tenxơ gradient biến dạng ở phạm vi biến dạng lớn

Khi vật liệu trải qua các biến dạng lớn, chẳng hạn như trong các nguyên công dập nóng, cán kim loại hoặc gia công áp lực cao, giả thiết phân rã cộng tính không còn bảo đảm tính đúng đắn về mặt hình học và cơ học. Đặt nền móng cho cơ học dẻo biến dạng hữu hạn năm 1969, Lee (1969) đã thiết lập khung lý thuyết phân rã nhân của tenxơ gradient biến dạng toàn phần $\mathbf{F}$ thông qua cấu hình trung gian không tải:

\mathbf{F} = \mathbf{F}^e \mathbf{F}^p

Trong cấu trúc phân rã này, $\mathbf{F}^p$ mô tả quá trình biến dạng dẻo không đổi thể tích chuyển từ cấu hình xuất phát sang cấu hình trung gian giải tỏa ứng suất, trong khi $\mathbf{F}^e$ bao hàm cả sự biến dạng đàn hồi hồi phục và sự giãn nở thể tích do nhiệt độ tác động. Dựa trên khung phân rã nhân của Lee, Simo và Miehe (1992) đã phát triển mô hình nhiệt dẻo kết hợp tương hỗ ở biến dạng hữu hạn, bảo đảm tính nhất quán nhiệt động lực học khắt khe theo bất đẳng thức Clausius - Duhem và thiết lập các thuật toán tích phân số ứng suất ổn định cho mô phỏng số.

Tiêu chuẩn chảy dẻo và quy tắc dòng chảy phụ thuộc nhiệt độ

Một đặc trưng chi phối cơ bản của ứng xử nhiệt đàn dẻo là sự phụ thuộc nhạy cảm của các tính chất cơ học vào nhiệt độ. Khi nhiệt độ môi trường hoặc cục bộ tăng cao, dao động nhiệt của các nút mạng tinh thể kim loại gia tăng, làm giảm đáng kể rào cản cản trở sự chuyển động của các lệch mạng. Do đó, giới hạn chảy của hầu hết các kim loại kỹ thuật có xu hướng suy giảm rõ rệt khi nhiệt độ tiệm cận điểm nóng chảy.

Hàm tiêu chuẩn chảy von Mises

Để xác định thời điểm một phân tố vật chất bắt đầu chuyển từ trạng thái biến dạng đàn hồi sang trạng thái dẻo trong điều kiện nhiệt độ biến đổi, lý thuyết nhiệt đàn dẻo sử dụng mặt chảy dẻo mở rộng. Hàm tiêu chuẩn chảy dẻo von Mises phụ thuộc nhiệt độ được biểu diễn dưới dạng:

f(\boldsymbol{\sigma}, \sigma_y, T) = \sigma_{eq} - \sigma_y(T) \le 0

Trong phương trình này, $\sigma_{eq}$ là ứng suất tương đương von Mises, được tính toán từ phần tenxơ ứng suất lệch $\mathbf{s}$ :

\sigma_{eq} = \sqrt{\frac{3}{2} \mathbf{s} : \mathbf{s}}

với $\mathbf{s} = \boldsymbol{\sigma} - \frac{1}{3}\text{tr}(\boldsymbol{\sigma})\mathbf{I}$ , trong đó $\text{tr}(\boldsymbol{\sigma})$ đại diện cho vết của tenxơ ứng suất. Giá trị $\sigma_y(T)$ là giới hạn chảy dẻo của vật liệu, đóng vai trò là một hàm số liên tục đơn điệu giảm theo nhiệt độ $T$ . Khi hàm $f < 0$ , vật liệu phản ứng hoàn toàn đàn hồi; khi $f = 0$ và điều kiện gia tải dẻo được thỏa mãn, biến dạng dẻo bắt đầu phát sinh và tích lũy.

Quy tắc dòng chảy dẻo và tương tác sinh nhiệt

Hướng phát triển của biến dạng dẻo được điều khiển bởi quy tắc dòng chảy liên kết (associative flow rule):

\dot{\boldsymbol{\varepsilon}}^p = \dot{\lambda} \frac{\partial f}{\partial \boldsymbol{\sigma}} = \dot{\lambda} \frac{3 \mathbf{s}}{2 \sigma_{eq}}

trong đó $\dot{\lambda}$ là nhân tử dẻo không âm thỏa mãn các điều kiện bù trừ Kuhn - Tucker. Đáng chú ý, công tiêu tán dẻo do quá trình biến dạng không thuận nghịch sinh ra không hoàn toàn biến mất mà chuyển hóa phần lớn thành nhiệt năng nội sinh, đóng vai trò như một nguồn nhiệt thể tích bên trong vật thể, làm biến thiên cục bộ trường nhiệt độ của chính kết cấu đang xét.

Phương pháp giải số phần tử hữu hạn nhiệt đàn dẻo

Do tính phi tuyến kép bao gồm phi tuyến vật liệu (tính dẻo, tính chất cơ lý phụ thuộc nhiệt độ), phi tuyến hình học (chuyển vị và biến dạng lớn) và liên kết đa trường vật lý (nhiệt - cơ), các phương trình vi phân của cơ học nhiệt đàn dẻo hầu như không thể tìm được nghiệm giải tích dạng đóng cho các kết cấu kỹ thuật thực tế. Việc giải quyết các bài toán này hoàn toàn phụ thuộc vào phương pháp phần tử hữu hạn (Finite Element Method - FEM).

Các khía cạnh tính toán số của phân tích ứng suất nhiệt đàn dẻo kết hợp đã được Argyris và cộng sự (1982) thiết lập một cách toàn diện trong mô phỏng kết cấu và phân tích ứng suất hàn. Trong khuôn khổ tính toán số, quy trình giải thường được phân chia thành hai chiến lược chính:

Phương pháp phân tách tuần tự (Sequentially Coupled Analysis): Trường nhiệt độ được giải trước độc lập thông qua phương trình truyền nhiệt phi tuyến nhằm thu được lịch sử nhiệt độ tại mọi nút lưới theo thời gian. Sau đó, trường nhiệt độ này được đưa vào như một trường tải trọng nhiệt ngoài để giải bài toán cơ học phi tuyến nhiệt đàn dẻo. Chiến lược này đạt hiệu quả tính toán cao và được ứng dụng rộng rãi khi năng lượng nhiệt sinh ra do biến dạng dẻo chiếm tỷ trọng không đáng kể so với nguồn nhiệt ngoài.
Phương pháp kết hợp hoàn toàn (Fully Coupled Analysis): Trường nhiệt độ và trường chuyển vị cơ học được giải đồng thời tại từng bước thời gian bằng thuật toán lặp Newton - Raphson kết hợp ma trận độ cứng tiếp tuyến nhất quán. Chiến lược này bắt buộc phải sử dụng khi tốc độ biến dạng rất cao hoặc hiện tượng sinh nhiệt do ma sát và biến dạng dẻo có ảnh hưởng ngược rõ rệt lên sự phân bố nhiệt độ của hệ thống.

Để mô phỏng chính xác trường nhiệt độ trong các quá trình gia công nhiệt tập trung, Goldak và cộng sự (1984) đã xây dựng mô hình phần tử hữu hạn nguồn nhiệt khối elip kép. Mô hình của Goldak giải quyết triệt để sự mất cân xứng của gradient nhiệt độ giữa phần trước và phần sau của vũng hàn di động, trở thành tiêu chuẩn kinh điển không thể thay thế trong các phần mềm mô phỏng phần tử hữu hạn nhiệt cơ hiện đại.

Ứng dụng thực tiễn trong kỹ thuật công nghệ

Lý thuyết nhiệt đàn dẻo là nền tảng tính toán không thể thiếu trong nhiều lĩnh vực công nghiệp chế tạo tiên tiến:

Dự đoán ứng suất dư và biến dạng trong công nghệ hàn kim loại

Trong quá trình hàn, sự nung nóng cục bộ cường độ cao và quá trình làm nguội nhanh không đồng đều làm phát sinh sự co giãn nhiệt bị kiềm chế, dẫn đến sự xuất hiện của các biến dạng dẻo dư và trường ứng suất dư kéo có độ lớn xấp xỉ giới hạn chảy của kim loại. Ứng suất dư này làm tăng nguy cơ nứt tế vi, biến dạng hình học của kết cấu và đẩy nhanh quá trình ăn mòn mỏi. Nghiên cứu của Han và cộng sự (2024) đã ứng dụng phương pháp phần tử hữu hạn nhiệt đàn dẻo có xét đến nhiệt độ giữa các lượt hàn để dự đoán chính xác sự phân bố ứng suất dư và biến dạng hàn của kết cấu kim loại, chứng minh rằng việc kiểm soát nhiệt độ giữa các lượt hàn là giải pháp cốt lõi để giảm thiểu biến dạng cong vênh của sản phẩm sau chế tạo.

Phân tích phá hủy và tập trung biến dạng ở nhiệt độ cao

Trong các thiết bị làm việc dưới điều kiện chu trình nhiệt cơ khắc nghiệt như tuabin khí, cánh quạt động cơ phản lực hoặc ống dẫn trong lò phản ứng hạt nhân, sự tích lũy biến dạng dẻo chu kỳ kết hợp với hiệu ứng mỏi nhiệt dẫn đến sự hình thành các dải trượt cục bộ và khởi sinh vết nứt. Baruah và Singh (2025) đã xây dựng mô hình tính toán phá hủy nhiệt đàn dẻo sử dụng lý thuyết gradient tổn hại cục bộ và phi cục bộ, mô phỏng chính xác sự tập trung biến dạng dẻo và quá trình phát triển vết nứt nẻ của vật liệu kim loại ở nhiệt độ cao mà không bị phụ thuộc vào kích thước phần tử lưới tính toán.

Công nghệ gia công bồi đắp kim loại (In 3D kim loại)

Tương tự như quá trình hàn nhiều lượt, công nghệ in 3D kim loại bằng chùm tia laser hoặc chùm tia điện tử liên tục làm nóng chảy và đông đặc từng lớp bột kim loại siêu mỏng. Chu kỳ gia nhiệt và làm nguội lặp đi lặp lại hàng nghìn lần tạo ra gradient nhiệt độ cực lớn trong thể tích chi tiết, dễ gây nứt tách lớp hoặc biến dạng làm mất dung sai lắp ghép. Việc áp dụng mô phỏng số nhiệt đàn dẻo cho phép các kỹ sư tối ưu hóa chiến lược quét laser, điều chỉnh tốc độ cấp nhiệt và dự đoán trước trường biến dạng để bù trừ hình học cho phôi in trước khi tiến hành chế tạo thực tế.