Khi phân tích tải trọng giới hạn bằng ANSYS, giá trị ứng suất Mises vượt quá giới hạn dẻo

Chủ bài dùng phương pháp giới hạn để tiến hành phân tích phần tử hữu hạn cho kết cấu vòng gia cố à? Tôi thường trực tiếp sử dụng ống dẫn dày hơn để tiến hành phân tích thông thường. Loại có vòng gia cố thì chỉ cần sử dụng phương pháp tính toán thông thường là được. Nếu vượt quá tiêu chuẩn, nên tăng độ dày của ống nối.

Đang sử dụng vật liệu đàn hồi lý tưởng phải không?

Đây là bài viết cũ từ bốn năm trước rồi. Gần đây tôi lại gặp phải vấn đề tương tự, nên xin chia sẻ suy nghĩ của mình để mọi người cùng thảo luận và tham khảo. Rất mong nhận được những ý kiến đóng góp và chỉ trích từ mọi người. Sau khi thiết lập giá trị giới hạn chảy trong Ansys, về mặt lý thuyết, khi ứng suất Mises tại điểm đó (có vẻ như Ansys không hỗ trợ nguyên tắc chảy dựa trên độ bền cấp ba) đạt đến giá trị giới hạn chảy, mối quan hệ giữa ứng suất và biến dạng tại khu vực đó sẽ thay đổi. Hình thức cụ thể của sự thay đổi này phụ thuộc vào mô hình mối quan hệ đặc tính vật liệu được sử dụng. Tôi sử dụng phương pháp tăng cường tuyến tính kép, với giá trị tangens được đặt thành 0. Trong điều kiện này, sau khi khối vật liệu bắt đầu chịu lực, sẽ không xảy ra hiện tượng cứng hóa; áp suất sẽ không tiếp tục tăng theo mức độ biến dạng, mà sẽ luôn giữ nguyên ở mức của đường biến dạng giới hạn. Trong quá trình tự thực hiện, tôi cũng gặp tình trạng giá trị ứng suất Mises vượt quá mức giới hạn chịu đựng, với mức chênh lệch khoảng 30%. Điều này là do trong phép tính phần tử hữu hạn, giá trị được tính trực tiếp chính là giá trị tại các điểm tích phân bên trong phần tử; sau đó, giá trị này được mở rộng ra bên ngoài (có thể so sánh với việc nội suy) để thu được dữ liệu của các nút. Giá trị ứng suất được hiển thị trong biểu đồ mây chính là giá trị tại các điểm nút. Sai số xảy ra trong quá trình dự báo có thể khiến ứng suất tại các điểm nút lớn hơn giá trị ứng suất tại các điểm tích phân ở hai bên. Theo quan điểm cá nhân, việc ước lượng các giá trị tại các điểm tích phân và sai số phát sinh từ việc ước lượng này là điều không thể tránh khỏi. Chúng ta có thể giảm bớt sai số đó bằng cách tăng độ chi tiết của mạng lưới, từ đó có thêm nhiều điểm tích phân hơn trong khu vực cần xem xét, và khoảng cách giữa các điểm này cũng sẽ ngắn hơn. Trong quá trình tính toán của bản thân tôi, phương pháp này cũng rất hiệu quả. Khi mạng lưới được phân chia chi tiết hơn, giá trị ứng suất tối đa vẫn vượt quá giá trị giới hạn chịu lực, nhưng độ lệch tương đối chỉ khoảng 1,5% mà thôi.

Nếu bạn tự thêm các vùng vào thang đo biểu diễn dạng đám mây, và thiết lập một vùng hẹp gần đường giới hạn chịu lực, thì trong bản đồ đám mây mới, các vùng có ứng suất vượt quá giá trị giới hạn sẽ nằm ở những chỗ có chất lượng mạng lưới kém hoặc ở rìa vùng chịu lực. Ở những nơi này, mức độ thay đổi của ứng suất khá lớn, do đó sai số khi tính toán cũng sẽ cao. Ansys trước tiên tính ứng suất biến dạng, sau đó mới tính ứng suất. Sau khi thực hiện các thao tác nêu trên, người ta sẽ nhận thấy rằng vùng mà ứng suất đạt đến giá trị giới hạn chịu lực và vùng xuất hiện biến dạng dẻo gần như trùng khớp với nhau. Điều này cũng cho thấy rằng các giá trị ứng suất lớn hơn giá trị giới hạn chịu lực được hiển thị trên bản đồ ứng suất thực chất là do sai số trong quá trình tính toán; những giá trị này không ảnh hưởng đến mức độ biến dạng của kết cấu. Dù sao thì việc tính toán biến dạng trước, sau đó mới tính toán ứng suất cũng là một quy trình hợp lý. Nói cách khác, bỏ qua những sai số phát sinh từ chính quá trình mô phỏng (do tải trọng, các điều kiện biên, việc đơn giản hóa mô hình và cấu trúc lưới), mặc dù ở một số vị trí trong kết quả tính toán, giá trị ứng suất lớn hơn giá trị giới hạn chịu đựng của vật liệu, nhưng những sai số cục bộ này không ảnh hưởng đến độ chính xác tổng thể của kết quả mô phỏng. Kết quả ứng suất của kết cấu dưới các điều kiện tính toán là đáng tin cậy.

Tiếp theo là vấn đề làm chi tiết hóa mạng lưới. Do có những hạn chế về hiệu năng của máy tính dùng để xử lý công việc, tôi đã sử dụng phương pháp mô hình con trong quá trình tính toán chi tiết. Dựa trên kết quả tính toán từ mô hình chính, tôi liên tục lấy ra các phần cụ thể trong cấu trúc cần được xem xét để tiến hành tính toán. Quả thật, càng có mạng lưới các mô hình con chi tiết hơn, thì độ lệch tương đối giữa ứng suất tối đa và đường uốn cũng càng nhỏ. Tuy nhiên, phương pháp mô hình con cũng cần chú ý tránh các vùng tập trung ứng suất, nhằm giúp dữ liệu định nghĩa các điều kiện biên của mô hình con trở nên chính xác nhất có thể. Vì vậy, mô hình con không nên quá nhỏ. Hơn nữa, dựa vào mối quan hệ giữa các giá trị ứng suất và biến dạng nêu trên, có thể kết luận rằng trong các phép tính thực tế, không cần thiết phải làm cho mạng lưới điểm dữ liệu trở nên chi tiết đến mức giá trị ứng suất tại các điểm đó bằng hoàn toàn giá trị giới hạn uốn. Chỉ cần các giá trị này gần bằng nhau thôi, thì vẫn có thể thu được kết quả về biến dạng và ứng suất với độ chính xác khá cao. Dù sao thì sai số trong quá trình mô phỏng cũng chắc chắn sẽ tồn tại; mức độ chi tiết của các mô hình con và lưới mô hình phụ thuộc vào mức độ chấp nhận sai số của bạn.

Trước đó tôi quên nói rồi, bằng cách so sánh các bản đồ ứng suất và biến dạng dẻo, có thể nhận thấy rằng ở những vị trí mà ứng suất lớn hơn nhiều so với giá trị giới hạn chảy, thì khả năng xảy ra biến dạng dẻo lại rất thấp. Theo hiểu biết cá nhân của tôi, khi Ansys tự đánh giá các kết quả như giá trị chịu uốn, có lẽ họ sử dụng các giá trị tại các điểm tích phân; còn những gì được hiển thị trong biểu đồ là các giá trị tại các nút. Vì vậy, bản đồ ứng suất có thể chứa sai số, nhưng bản đồ ứng biến thì nên tương đối chính xác; điều này không mâu thuẫn gì cả. Tất cả những nội dung được đề cập ở trên, bao gồm sai số do việc ước lượng các điểm tích điểm, phương pháp sử dụng các mô hình con, việc xác định giới hạn chịu uốn, v.v., đều có thể dễ dàng tìm thấy. Nếu có điều gì sai sót, xin mọi người hãy vui lòng chỉ ra để tôi có thể sửa chữa.

Bài viết này cuối cùng được chỉnh sửa bởi Ou Lan vào ngày 14-12-2020 lúc 18:04. Cấu trúc áp suất bên trong ống tròn khá đơn giản, mạng lưới các thành phần cấu tạo có trật tự, và sự phân bố ứng suất và biến dạng cũng khá lý tưởng; vì vậy, khó có thể xảy ra những vấn đề mà tôi đã đề cập trước đây. Bạn hãy thử kéo giãn đơn thuần xem, bạn sẽ nhận thấy rằng mức độ uốn cong của mẫu thí nghiệm cũng khá phù hợp với giá trị lý tưởng.