Tổng quan nghiên cứu

Trong kỹ thuật kết cấu hiện đại, thép mềm và kim loại chịu lực chiếm hơn 70% tổng khối lượng vật liệu xây dựng cầu đường, nhà cao tầng và các thiết bị cơ khí công nghiệp. Các loại thép kết cấu thông dụng thường có giới hạn tỷ lệ đàn hồi trong khoảng từ 200 MPa đến 280 MPa, trong khi thép cường độ cao có thể đạt trên 550 MPa. Khi công trình chịu tải trọng lớn vượt qua ngưỡng đàn hồi ban đầu, vật liệu sẽ chuyển sang trạng thái chảy dẻo với độ giãn dài lớn gấp 10 đến 15 lần so với giai đoạn đàn hồi. Nếu kỹ sư chỉ áp dụng các mô hình tuyến tính truyền thống, sai số dự báo khả năng chịu lực có thể vượt quá 25%, dẫn đến lãng phí vật tư hoặc gây mất an toàn nghiêm trọng cho công trình.

Luận văn tập trung giải quyết bài toán phân tích phi tuyến vật liệu đối với kết cấu dầm chịu tải, nhằm khắc phục triệt để các hạn chế của phương pháp giải tích cổ điển vốn không thể giải được các phương trình vi phân phi tuyến phức tạp. Mục tiêu cốt lõi của đề tài là xây dựng công thức phần tử dầm Bernoulli hai nút, thiết lập ma trận độ cứng tiếp tuyến, phát triển thuật toán cập nhật ứng suất và lập trình hệ thống tính toán số hoàn chỉnh trên nền tảng MATLAB. Nghiên cứu được triển khai tại Viện Cơ học thuộc Viện Khoa học và Công nghệ Việt Nam phối hợp cùng Trường Đại học Công nghệ thuộc Đại học Quốc gia Hà Nội. Kết quả nghiên cứu mang lại ý nghĩa kinh tế và kỹ thuật to lớn, cho phép tối ưu hóa kích thước mặt cắt dầm từ 10% đến 15%, đồng thời bảo đảm hệ số an toàn kết cấu luôn đạt trên 1,5 lần dưới các điều kiện tải trọng tập trung và mô-men uốn khắc nghiệt.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng hệ thống lý thuyết cơ học vật thể rắn biến dạng phi tuyến, trọng tâm là lý thuyết dẻo và lý thuyết dầm Bernoulli. Khi vật liệu làm việc trong miền đàn hồi, mối quan hệ giữa ứng suất và biến dạng tuân thủ định luật Hooke với mô-đun đàn hồi Young của thép xây dựng đạt mức xấp xỉ 210.000 MPa, trong khi biến dạng đàn hồi cực hạn duy trì ở mức dưới 5%. Khi tải trọng ngoài vượt quá giới hạn chảy, ứng xử của kết cấu được mô tả qua mô hình vật liệu đàn-dẻo lưỡng tuyến tính với mô-đun tiếp tuyến đặc trưng cho giai đoạn tái bền.

Ba trụ cột lý thuyết dẻo được tích hợp trong mô hình gồm có: tiêu chuẩn chảy dẻo von Mises để xác định thời điểm bắt đầu chảy dẻo trong không gian ứng suất, luật chảy Prandtl-Reuss thiết lập mối quan hệ giữa gia số biến dạng dẻo với ten-xơ ứng suất lệch, và hai quy luật tái bền gồm tái bền đẳng hướng cùng tái bền động học nhằm phản ánh chính xác hiệu ứng Bauschinger khi kết cấu chịu chu trình chất tải và cất tải đảo chiều.

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

Phương pháp nghiên cứu chủ đạo là phương pháp phần tử hữu hạn phi tuyến kết hợp thuật toán số giải tích gia số. Nguồn dữ liệu kiểm chứng dựa trên các thông số thực nghiệm kéo mẫu kim loại theo tiêu chuẩn quốc tế ASTM với đường kính mẫu đo chuẩn 12,8 mm và chiều dài đo hiệu dụng 50,8 mm. Cỡ mẫu mô phỏng bao gồm 3 mô hình kết cấu dầm điển hình: dầm tựa đơn chịu lực tập trung tại điểm giữa, dầm công-xôn chịu lực tập trung tại đầu tự do, và dầm công-xôn chịu mô-men uốn tại đầu mút.

Phương pháp chọn mẫu điển hình được áp dụng nhằm đại diện trọn vẹn cho các trạng thái chịu uốn thuần túy và uốn kết hợp lực cắt trong thực tế xây dựng. Lý do lựa chọn phương pháp phần tử hữu hạn là do khả năng xử lý bài toán tiêu tán năng lượng khi nguyên lý thế năng dừng không còn khả dụng, đòi hỏi phải áp dụng nguyên lý công ảo để rời rạc hóa kết cấu. Ma trận độ cứng tiếp tuyến của phần tử dầm được tính toán chính xác thông qua phép tích phân số Gauss hai chiều với lưới điểm gồm 4 điểm theo trục dọc và 7 điểm theo chiều cao tiết diện (tổng cộng 28 điểm Gauss trên mỗi phần tử). Quy trình giải hệ phương trình cân bằng phi tuyến sử dụng thuật toán lặp gia số Newton-Raphson với tiêu chuẩn hội tụ dựa trên chuẩn sai số lực dư đạt độ chính xác khắt khe là 0,0001 (tức 10 mũ trừ 4).

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng số và phân tích kết cấu dầm đàn-dẻo đã đem lại nhiều phát hiện định lượng quan trọng về mặt cơ học:

Thứ nhất, đồ thị phản ứng lực và chuyển vị cho thấy khi tải trọng ngoài vượt qua điểm chảy ban đầu, độ võng tại điểm giữa của dầm tựa đơn gia tăng nhanh chóng, đạt mức tăng hơn 180% so với giả định đàn hồi tuyến tính thuần túy ở cùng một mức tải trọng lớn.

Thứ hai, tham số mô-đun tiếp tuyến đóng vai trò quyết định đến độ cứng dư và khả năng chịu tải sau đàn hồi của dầm. Khi mô-đun tiếp tuyến đạt giá trị 0,3 lần mô-đun đàn hồi Young, độ võng ở đầu tự do của dầm công-xôn giảm khoảng 35% so với trường hợp vật liệu có mô-đun tiếp tuyến chỉ bằng 0,1 lần mô-đun đàn hồi Young dưới cùng một giá trị lực tác dụng.

Thứ ba, thuật toán cập nhật ứng suất ba trường hợp (đã chảy dẻo trước đó, giữ nguyên trạng thái đàn hồi, và chuyển tiếp từ đàn hồi sang dẻo) đạt độ ổn định tuyệt đối. Quy trình lặp Newton-Raphson chỉ cần từ 4 đến 6 bước lặp tại mỗi bước tăng tải để chuẩn lực dư suy giảm hơn 99,99%, đáp ứng chuẩn dung sai 0,0001.

Thảo luận kết quả

Sự suy giảm độ cứng tổng thể của dầm khi tải trọng tăng được giải thích bởi sự lan rộng của vùng chảy dẻo từ các thớ ngoài cùng của tiết diện (nơi có khoảng cách lớn nhất tới trục trung hòa) dần vào thớ trung tâm. Dữ liệu tính toán từ chương trình số được trình bày trực quan thông qua các bảng thống kê ứng suất, biến dạng tại điểm Gauss thứ nhất (như Bảng 3.1, Bảng 3.2 và Bảng 3.3) cùng các đường cong quan hệ tải trọng và độ võng. Khi so sánh với các nghiên cứu giải tích đàn-dẻo cổ điển, kết quả từ chương trình số cho thấy sự trùng khớp với sai lệch dưới 2% trong miền đàn hồi và phản ánh mượt mà bước chuyển pha phi tuyến. Việc kiểm soát chính xác 28 điểm tích phân Gauss trên mỗi phần tử giúp chương trình bắt trọn quá trình hình thành khớp dẻo cục bộ, loại bỏ hoàn toàn hiện tượng khóa số hoặc phân kỳ thuật toán thường gặp trong các mô hình đơn giản hóa.

Đề xuất và khuyến nghị

Nhằm chuyển giao hiệu quả các đóng góp học thuật của luận văn vào thực tiễn thiết kế cơ học và công trình, các đề xuất cụ thể được kiến nghị như sau:

Thứ nhất, tích hợp mô hình phần tử dầm đàn-dẻo Bernoulli vào phần mềm thiết kế kết cấu công trình. Bộ phận nghiên cứu và phát triển tại các tổng công ty tư vấn thiết kế xây dựng cần triển khai ứng dụng thuật toán trong vòng 6 tháng tới, đặt mục tiêu nâng cao độ chính xác khi tính toán khả năng chịu tải dự trữ của khung dầm thép lên thêm 20%.

Thứ hai, mở rộng thuật toán sang bài toán phi tuyến hình học kết hợp phi tuyến vật liệu. Nhóm nghiên cứu chuyên sâu tại Viện Cơ học nên tiếp tục hoàn thiện mã nguồn trong thời gian 12 tháng, nhằm kiểm soát các hiện tượng mất ổn định tổng thể và biến dạng lớn với sai số thuật toán khống chế dưới 1,5%.

Thứ ba, chuẩn hóa quy trình thí nghiệm xác định mô-đun tiếp tuyến vật liệu. Các phòng thí nghiệm cơ học vật liệu tại các viện và trường đại học kỹ thuật cần chuẩn hóa phép thử kéo nén chu kỳ theo tiêu chuẩn ASTM trong quý 3 năm 2026, giúp cung cấp thông số đầu vào chính xác với độ tin cậy trên 95% cho các phân tích phần tử hữu hạn.

Thứ tư, đổi mới chương trình đào tạo sau đại học ngành Cơ học và Xây dựng. Các trường đại học kỹ thuật nên đưa chuyên đề thuật toán phần tử hữu hạn phi tuyến và mô hình hóa dẻo vào chương trình giảng dạy thạc sĩ từ năm học 2026-2027, phấn đấu tăng tỷ lệ sinh viên có khả năng tự phát triển chương trình tính toán số lên 30%.

Đối tượng nên tham khảo luận văn

Luận văn là tài liệu tham khảo giá trị cao cho 4 nhóm đối tượng chính:

  1. Kỹ sư thiết kế kết cấu công trình dân dụng, cầu đường và công trình biển: Tài liệu cung cấp cơ sở tính toán chính xác để đánh giá khả năng chịu tải của dầm thép sau khi vượt giới hạn chảy đàn hồi, giúp tối ưu hóa vật liệu trong các dự án quy mô lớn.
  2. Học viên cao học và nghiên cứu sinh chuyên ngành Cơ học vật thể rắn và Cơ kỹ thuật: Nghiên cứu đóng vai trò như một cẩm nang phương pháp luận hoàn chỉnh từ việc thiết lập phương trình công ảo, xây dựng ma trận độ cứng tiếp tuyến đến thuật toán Newton-Raphson.
  3. Kỹ sư phát triển phần mềm mô phỏng kỹ thuật CAE và FEA: Các lập trình viên có thể kế thừa và thương mại hóa thuật toán cập nhật ứng suất cùng cấu trúc hàm xử lý tích phân Gauss 28 điểm viết bằng ngôn ngữ MATLAB.
  4. Giảng viên và nhà nghiên cứu tại các viện cơ học và trường đại học kỹ thuật: Tài liệu là nguồn học liệu mẫu mực để biên soạn bài giảng môn học Phương pháp phần tử hữu hạn nâng cao và Cơ học biến dạng dẻo.

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

  1. Tại sao cần phân tích đàn-dẻo thay vì chỉ tính toán đàn hồi tuyến tính truyền thống? Phân tích đàn hồi tuyến tính bỏ qua khả năng chịu lực sau khi vật liệu đạt giới hạn chảy từ 200 MPa đến 280 MPa. Mô hình đàn-dẻo giúp khai thác tối đa vùng biến dạng dẻo có kiểm soát, nâng cao hiệu quả sử dụng vật liệu từ 10% đến 15% mà vẫn bảo đảm công trình an toàn tuyệt đối.

  2. Tiêu chuẩn von Mises và luật chảy Prandtl-Reuss đóng vai trò gì trong mô hình dầm? Tiêu chuẩn von Mises xác định chính xác thời điểm ứng suất hữu hiệu đạt giới hạn chảy ban đầu. Kết hợp với luật chảy Prandtl-Reuss, mô hình tính toán chuẩn xác tỷ lệ gia số biến dạng dẻo, giúp thuật toán cập nhật đúng trạng thái cơ học tại từng điểm Gauss của dầm.

  3. Thuật toán Newton-Raphson giải quyết tính phi tuyến của kết cấu như thế nào? Thuật toán chia nhỏ tải trọng tác dụng thành các bước gia số và tuyến tính hóa phương trình cân bằng phi tuyến quanh điểm hiện tại. Thông qua ma trận độ cứng tiếp tuyến, thuật toán lặp từ 4 đến 6 bước để triệt tiêu lực dư với độ chính xác đạt chuẩn sai số 0,0001.

  4. Vì sao cần sử dụng lưới tích phân 28 điểm Gauss trên mỗi phần tử dầm? Tiết diện dầm phẳng chịu uốn có sự biến thiên ứng suất phức tạp theo cả chiều dài và chiều cao. Việc bố trí 4 điểm Gauss theo phương dọc và 7 điểm theo chiều dày cho phép theo dõi sát sao sự hình thành từng lớp dẻo cục bộ với độ chính xác trên 99%.

  5. Hiệu ứng Bauschinger được xử lý ra sao trong quy luật tái bền của nghiên cứu? Nghiên cứu ứng dụng luật tái bền động học để mô tả hiện tượng tâm mặt dẻo dịch chuyển trong không gian ứng suất. Nhờ đó, giới hạn chảy khi nén ngược chiều sau khi kéo dẻo được tái hiện chính xác, hạn chế tối đa sai số khi công trình chịu tải trọng đổi chiều.

Kết luận

  • Xây dựng thành công công thức phần tử hữu hạn dầm Bernoulli hai nút tích hợp đầy đủ ứng xử đàn-dẻo phi tuyến vật liệu.
  • Thiết lập hoàn chỉnh ma trận độ cứng tiếp tuyến đối xứng và thuật toán lặp tăng dần Newton-Raphson với tiêu chuẩn hội tụ 0,0001.
  • Phát triển trọn bộ mã nguồn tính toán số trên nền tảng MATLAB với cơ chế cập nhật ứng suất tự động tại 28 điểm tích phân Gauss.
  • Kiểm chứng độ chính xác qua 3 bài toán dầm tựa đơn và dầm công-xôn, chứng minh tính vượt trội trong việc mô phỏng khớp dẻo so với lý thuyết cổ điển.
  • Đóng góp giải pháp khoa học giúp tiết kiệm từ 10% đến 15% vật liệu kết cấu và bảo đảm hệ số an toàn công trình trên 1,5 lần.

Về lộ trình phát triển tiếp theo, nhóm tác giả định hướng hoàn thiện mô hình dầm không gian 3 chiều và tích hợp bài toán động lực học phi tuyến trong vòng 12 tháng tới. Các kỹ sư và nhà nghiên cứu được khuyến khích áp dụng rộng rãi thuật toán này vào thực tiễn thiết kế nhằm nâng cao chất lượng và độ an toàn bền vững cho các công trình kỹ thuật.