Chênh nhiệt độ chưa đủ để xác định ứng suất
Một thành thiết bị nóng không đều có thể vừa giãn dài, vừa cong, vừa mang ứng suất tự cân bằng. Lấy EαΔT với ΔT là chênh nhiệt độ hai mặt sẽ bỏ qua điều kiện biến dạng. Một dải được phép cong khác một thành vỏ khép kín hoặc một chi tiết bị giữ phẳng, dù các nhiệt kế trên hai mặt chỉ cùng giá trị.
Để nhìn rõ cơ chế, bài này xét một dải chữ nhật dài, đàn hồi tuyến tính và đồng nhất. Nó là mô hình minh họa của một dải vật liệu, không phải nghiệm ứng suất vòng của toàn ống. Phép cắt tưởng tượng thành ống thành dải đã bỏ sự liên tục quanh chu vi; phải trả lại liên kết đó khi tính thiết bị thực.
Biến dạng nhiệt và biến dạng cơ học
Gọi z là tọa độ qua chiều dày t, từ −t/2 tới t/2; x dọc dải. Theo giả thiết tiết diện phẳng, biến dạng dọc có dạng ε_x(z)=ε_0+κz. Quy ước κ là hệ số gradient của biến dạng dọc; dấu độ cong hình học có thể khác tùy quy ước chuyển vị. E là môđun Young, Pa; α là hệ số giãn nở, K⁻¹; ΔT(z) là nhiệt độ tăng so với trạng thái không ứng suất.
σ_x(z)=E[ε_0+κz−αΔT(z)].
Đây là ứng suất một chiều của mô hình dải, chưa bao gồm tương tác Poisson của một trạng thái ứng suất hai hoặc ba chiều. DSPE phân biệt biến dạng tổng với phần cơ học và phần nhiệt, đồng thời tách ảnh hưởng giãn dọc khỏi chuyển vị do gradient ngang tiết diện. [1] Không được lấy toàn bộ ε_x nhân E khi đã có giãn nở nhiệt.
Với dải tự do về lực dọc và mômen, cân bằng đòi hỏi N=∫σ_x dA=0 và M=∫zσ_x dA=0. Với bề rộng không đổi, tính chất không đổi và gốc z ở giữa chiều dày:
ε_0=(α/t)∫ΔT(z)dz.
κ=(12α/t³)∫zΔT(z)dz.
Các tích phân chạy từ −t/2 tới t/2. Điều kiện tự do là hai tích phân lực bằng không; nó không bắt buộc σ_x bằng không tại mọi điểm.
Ba phần của phân bố nhiệt độ
Đặt ξ=2z/t, nằm trong khoảng −1 tới1, và viết một phân bố giả định:
ΔT(ξ)=a_0+a_1ξ+a_2ξ².
a_0, a_1, a_2 có đơn vị K. Trung bình nhiệt độ tăng là a_0+a_2/3. Hai phương trình cân bằng cho ε_0=α(a_0+a_2/3) và κ=2αa_1/t. Thay lại công thức ứng suất:
σ_x(ξ)=Eαa_2(1/3−ξ²).
Phần nhiệt độ đều được giải phóng bằng giãn dài. Phần tuyến tính được giải phóng bằng độ cong. Phần bậc hai không khớp với một trường biến dạng tuyến tính qua tiết diện nên còn ứng suất. Với a_2=0 và các điều kiện tự do đã nêu, mô hình cho ứng suất bằng không dù nhiệt độ hai mặt khác nhau.
Ví dụ kiểm chứng của COMSOL cho thấy gradient nhiệt qua tiết diện gây uốn một dầm với các điều kiện liên kết được khai báo cụ thể. [2] Kết luận ấy giúp hiểu cơ chế, nhưng không cho phép mặc định mọi mô hình “không có lực ngoài” đều không ứng suất. Cấu trúc bị giữ, hình học vỏ hoặc phân bố nhiệt phi tuyến có thể thay đổi đáp án.
Ví dụ số với một phân bố phi tuyến
Giả định b=0,100 m, t=0,020 m, L=1,50 m, E=190 GPa, α=12×10⁻⁶ K⁻¹; chọn a_0=40 K, a_1=30 K và a_2=60 K. Đây là ảnh chụp một trường nhiệt giả định tại một thời điểm, không phải nghiệm dẫn nhiệt ổn định đã chứng minh cho cấu tạo mái hoặc thành bình nào.
Nhiệt độ tăng tại mặt ξ=−1 là70 K, giữa dải40 K và mặt ξ=1 là130 K. Trung bình qua chiều dày là60 K. Biến dạng đường giữa bằng0,000720, tương ứng tăng chiều dài1,080 mm; κ=0,036000 m⁻¹.
Nếu chỉ gá hai điểm đầu để tạo đường dây cung và không cản giãn dọc hay góc quay, độ võng lớn nhất theo xấp xỉ góc nhỏ có độ lớn |κ|L²/8=10,125 mm. Giá trị |κ|L/2=0,027 rad ở đầu giúp kiểm mức độ nhỏ của góc trong ví dụ. Không gán phép tính dây cung này cho một ống liên tục qua nhiều gối.
Ứng suất tại giữa chiều dày là +45,600 MPa, còn cả hai mặt là −91,200 MPa, với kéo dương. Hai mặt có nhiệt độ khác nhau nhưng ứng suất bằng nhau ở mô hình tự do vì phần chênh tuyến tính đã đi vào độ cong. Hai điểm đổi dấu nằm tại |z|=t/(2√3)=5,773503 mm.
Tích phân của 1/3−ξ² trên −1 tới1 bằng0; tích phân của ξ(1/3−ξ²) cũng bằng0. Đó là kiểm độc lập cho N và M. Chỉ đọc ứng suất mặt rồi nhân diện tích sẽ tạo ra một lực nén giả, vì đã bỏ phần kéo ở vùng giữa.
Giữ phẳng tạo thêm mômen
Xét một trường hợp khác: vẫn cho dải giãn dài tự do, nhưng hệ gá bắt κ=0. Điều kiện N=0 giữ ε_0=α60, trong khi ứng suất trở thành Eα(a_2/3−a_1ξ−a_2ξ²). Tại ξ=−1,0,1, ứng suất lần lượt −22,800; +45,600; −159,600 MPa. Cực đại kéo không còn ở giữa, mà tại ξ=−a_1/(2a_2)=−0,25, bằng54,150 MPa.
Mômen quán tính I=bt³/12=6,666667×10⁻⁸ m⁴. Mômen bên trong cần để ngăn độ cong nhiệt là M=−EIκ_th=−456,000 N·m theo quy ước đã chọn. Các liên kết phải cung cấp tác động tương thích với mômen đó; “giữ phẳng” không phải một điều kiện miễn phí về lực.
Nếu đồng thời giữ cả ε_0=0, ứng suất bằng −EαΔT. Khi ấy lực dọc N=−273,600 kN và mômen vẫn −456,000 N·m. Ứng suất tại mặt nóng là −296,400 MPa. Kết quả này là dự đoán đàn hồi có điều kiện; chưa có giới hạn chảy hoặc cơ sở thiết kế để kết luận vật liệu thực còn đàn hồi.
Các độ nhạy dễ bị diễn giải nhầm
Giữ phân bố theo ξ như cũ nhưng giảm chiều dày một nửa, độ cong tự do tăng gấp đôi. Ứng suất tự cân bằng do a_2 trong mô hình vẫn như cũ, còn mômen giữ phẳng giảm còn một phần tư vì I tỷ lệ t³ và κ tỷ lệ t⁻¹. Nếu thay vào đó giữ gradient K/m không đổi, hệ số a_1 phải đổi theo chiều dày và độ cong không còn tăng theo quy luật vừa nêu.
Giảm riêng a_2 từ60 xuống30 K làm ứng suất tự cân bằng giảm một nửa, nhưng không đổi κ vì a_1 giữ nguyên. Giảm E làm các ứng suất và mômen đàn hồi giảm tương ứng, trong khi độ cong tự do với α giữ nguyên không đổi. Đây là ba phép thay đổi riêng biệt; không nên gộp thành phát biểu chung rằng thành mỏng hơn hay vật liệu mềm hơn luôn làm giảm mọi tác động nhiệt.
Trở lại thành thiết bị thật
Thành ống còn ràng buộc quanh chu vi, điều kiện đầu, áp suất và chỗ chuyển tiếp độ dày. Với tấm rộng, tương tác theo hai phương cùng hình học phi tuyến có thể quan trọng. Bower phân tích riêng tấm tròn tự do chịu gradient xuyên chiều dày và chỉ ra độ cong có thể gắn với các trạng thái cân bằng phi tuyến. [3] Kết quả dải một chiều không bao phủ bài toán đó.
Một phân tích thiết bị cần trường nhiệt theo thời gian, nhiệt độ không ứng suất, tính chất theo nhiệt độ và điều kiện liên kết. Mô hình dải giúp kiểm dấu, đơn vị và phân biệt thành phần giãn, uốn, tự cân bằng. Sau kiểm ấy, phải giải bài toán vỏ hoặc khối phù hợp trước khi dùng ứng suất để đánh giá mỏi nhiệt hay khả năng làm việc.
Nguồn tham khảo
[1] 1.3 Principles of thermal deformation — Dutch Society for Precision Engineering. Truy cập ngày 09/11/2026.
[2] Thermally Loaded Beam — COMSOL Multiphysics 6.4 — COMSOL. Truy cập ngày 09/11/2026.
[3] Applied Mechanics of Solids — Section 10.7.5, Bending of a circular plate caused by a through-thickness temperature gradient — Allan F. Bower, Brown University. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một thành thiết bị nóng không đều có thể vừa giãn dài, vừa cong, vừa mang ứng suất tự cân bằng. Lấy EαΔT với ΔT là chênh nhiệt độ hai mặt sẽ bỏ qua điều kiện biến dạng. Một dải được phép cong khác một thành vỏ khép kín hoặc một chi tiết bị giữ phẳng, dù các nhiệt kế trên hai mặt chỉ cùng giá trị.
Để nhìn rõ cơ chế, bài này xét một dải chữ nhật dài, đàn hồi tuyến tính và đồng nhất. Nó là mô hình minh họa của một dải vật liệu, không phải nghiệm ứng suất vòng của toàn ống. Phép cắt tưởng tượng thành ống thành dải đã bỏ sự liên tục quanh chu vi; phải trả lại liên kết đó khi tính thiết bị thực.
Biến dạng nhiệt và biến dạng cơ học
Gọi z là tọa độ qua chiều dày t, từ −t/2 tới t/2; x dọc dải. Theo giả thiết tiết diện phẳng, biến dạng dọc có dạng ε_x(z)=ε_0+κz. Quy ước κ là hệ số gradient của biến dạng dọc; dấu độ cong hình học có thể khác tùy quy ước chuyển vị. E là môđun Young, Pa; α là hệ số giãn nở, K⁻¹; ΔT(z) là nhiệt độ tăng so với trạng thái không ứng suất.
σ_x(z)=E[ε_0+κz−αΔT(z)].
Đây là ứng suất một chiều của mô hình dải, chưa bao gồm tương tác Poisson của một trạng thái ứng suất hai hoặc ba chiều. DSPE phân biệt biến dạng tổng với phần cơ học và phần nhiệt, đồng thời tách ảnh hưởng giãn dọc khỏi chuyển vị do gradient ngang tiết diện. [1] Không được lấy toàn bộ ε_x nhân E khi đã có giãn nở nhiệt.
Với dải tự do về lực dọc và mômen, cân bằng đòi hỏi N=∫σ_x dA=0 và M=∫zσ_x dA=0. Với bề rộng không đổi, tính chất không đổi và gốc z ở giữa chiều dày:
ε_0=(α/t)∫ΔT(z)dz.
κ=(12α/t³)∫zΔT(z)dz.
Các tích phân chạy từ −t/2 tới t/2. Điều kiện tự do là hai tích phân lực bằng không; nó không bắt buộc σ_x bằng không tại mọi điểm.
Ba phần của phân bố nhiệt độ
Đặt ξ=2z/t, nằm trong khoảng −1 tới1, và viết một phân bố giả định:
ΔT(ξ)=a_0+a_1ξ+a_2ξ².
a_0, a_1, a_2 có đơn vị K. Trung bình nhiệt độ tăng là a_0+a_2/3. Hai phương trình cân bằng cho ε_0=α(a_0+a_2/3) và κ=2αa_1/t. Thay lại công thức ứng suất:
σ_x(ξ)=Eαa_2(1/3−ξ²).
Phần nhiệt độ đều được giải phóng bằng giãn dài. Phần tuyến tính được giải phóng bằng độ cong. Phần bậc hai không khớp với một trường biến dạng tuyến tính qua tiết diện nên còn ứng suất. Với a_2=0 và các điều kiện tự do đã nêu, mô hình cho ứng suất bằng không dù nhiệt độ hai mặt khác nhau.
Ví dụ kiểm chứng của COMSOL cho thấy gradient nhiệt qua tiết diện gây uốn một dầm với các điều kiện liên kết được khai báo cụ thể. [2] Kết luận ấy giúp hiểu cơ chế, nhưng không cho phép mặc định mọi mô hình “không có lực ngoài” đều không ứng suất. Cấu trúc bị giữ, hình học vỏ hoặc phân bố nhiệt phi tuyến có thể thay đổi đáp án.
Ví dụ số với một phân bố phi tuyến
Giả định b=0,100 m, t=0,020 m, L=1,50 m, E=190 GPa, α=12×10⁻⁶ K⁻¹; chọn a_0=40 K, a_1=30 K và a_2=60 K. Đây là ảnh chụp một trường nhiệt giả định tại một thời điểm, không phải nghiệm dẫn nhiệt ổn định đã chứng minh cho cấu tạo mái hoặc thành bình nào.
Nhiệt độ tăng tại mặt ξ=−1 là70 K, giữa dải40 K và mặt ξ=1 là130 K. Trung bình qua chiều dày là60 K. Biến dạng đường giữa bằng0,000720, tương ứng tăng chiều dài1,080 mm; κ=0,036000 m⁻¹.
Nếu chỉ gá hai điểm đầu để tạo đường dây cung và không cản giãn dọc hay góc quay, độ võng lớn nhất theo xấp xỉ góc nhỏ có độ lớn |κ|L²/8=10,125 mm. Giá trị |κ|L/2=0,027 rad ở đầu giúp kiểm mức độ nhỏ của góc trong ví dụ. Không gán phép tính dây cung này cho một ống liên tục qua nhiều gối.
Ứng suất tại giữa chiều dày là +45,600 MPa, còn cả hai mặt là −91,200 MPa, với kéo dương. Hai mặt có nhiệt độ khác nhau nhưng ứng suất bằng nhau ở mô hình tự do vì phần chênh tuyến tính đã đi vào độ cong. Hai điểm đổi dấu nằm tại |z|=t/(2√3)=5,773503 mm.
Tích phân của 1/3−ξ² trên −1 tới1 bằng0; tích phân của ξ(1/3−ξ²) cũng bằng0. Đó là kiểm độc lập cho N và M. Chỉ đọc ứng suất mặt rồi nhân diện tích sẽ tạo ra một lực nén giả, vì đã bỏ phần kéo ở vùng giữa.
Giữ phẳng tạo thêm mômen
Xét một trường hợp khác: vẫn cho dải giãn dài tự do, nhưng hệ gá bắt κ=0. Điều kiện N=0 giữ ε_0=α60, trong khi ứng suất trở thành Eα(a_2/3−a_1ξ−a_2ξ²). Tại ξ=−1,0,1, ứng suất lần lượt −22,800; +45,600; −159,600 MPa. Cực đại kéo không còn ở giữa, mà tại ξ=−a_1/(2a_2)=−0,25, bằng54,150 MPa.
Mômen quán tính I=bt³/12=6,666667×10⁻⁸ m⁴. Mômen bên trong cần để ngăn độ cong nhiệt là M=−EIκ_th=−456,000 N·m theo quy ước đã chọn. Các liên kết phải cung cấp tác động tương thích với mômen đó; “giữ phẳng” không phải một điều kiện miễn phí về lực.
Nếu đồng thời giữ cả ε_0=0, ứng suất bằng −EαΔT. Khi ấy lực dọc N=−273,600 kN và mômen vẫn −456,000 N·m. Ứng suất tại mặt nóng là −296,400 MPa. Kết quả này là dự đoán đàn hồi có điều kiện; chưa có giới hạn chảy hoặc cơ sở thiết kế để kết luận vật liệu thực còn đàn hồi.
Các độ nhạy dễ bị diễn giải nhầm
Giữ phân bố theo ξ như cũ nhưng giảm chiều dày một nửa, độ cong tự do tăng gấp đôi. Ứng suất tự cân bằng do a_2 trong mô hình vẫn như cũ, còn mômen giữ phẳng giảm còn một phần tư vì I tỷ lệ t³ và κ tỷ lệ t⁻¹. Nếu thay vào đó giữ gradient K/m không đổi, hệ số a_1 phải đổi theo chiều dày và độ cong không còn tăng theo quy luật vừa nêu.
Giảm riêng a_2 từ60 xuống30 K làm ứng suất tự cân bằng giảm một nửa, nhưng không đổi κ vì a_1 giữ nguyên. Giảm E làm các ứng suất và mômen đàn hồi giảm tương ứng, trong khi độ cong tự do với α giữ nguyên không đổi. Đây là ba phép thay đổi riêng biệt; không nên gộp thành phát biểu chung rằng thành mỏng hơn hay vật liệu mềm hơn luôn làm giảm mọi tác động nhiệt.
Trở lại thành thiết bị thật
Thành ống còn ràng buộc quanh chu vi, điều kiện đầu, áp suất và chỗ chuyển tiếp độ dày. Với tấm rộng, tương tác theo hai phương cùng hình học phi tuyến có thể quan trọng. Bower phân tích riêng tấm tròn tự do chịu gradient xuyên chiều dày và chỉ ra độ cong có thể gắn với các trạng thái cân bằng phi tuyến. [3] Kết quả dải một chiều không bao phủ bài toán đó.
Một phân tích thiết bị cần trường nhiệt theo thời gian, nhiệt độ không ứng suất, tính chất theo nhiệt độ và điều kiện liên kết. Mô hình dải giúp kiểm dấu, đơn vị và phân biệt thành phần giãn, uốn, tự cân bằng. Sau kiểm ấy, phải giải bài toán vỏ hoặc khối phù hợp trước khi dùng ứng suất để đánh giá mỏi nhiệt hay khả năng làm việc.
Nguồn tham khảo
[1] 1.3 Principles of thermal deformation — Dutch Society for Precision Engineering. Truy cập ngày 09/11/2026.
[2] Thermally Loaded Beam — COMSOL Multiphysics 6.4 — COMSOL. Truy cập ngày 09/11/2026.
[3] Applied Mechanics of Solids — Section 10.7.5, Bending of a circular plate caused by a through-thickness temperature gradient — Allan F. Bower, Brown University. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
