Tổng quan nghiên cứu

Trong kỹ thuật xây dựng hiện đại, vật liệu thép đóng vai trò then chốt với mô đun đàn hồi danh định đạt khoảng 200.000 MPa cùng giới hạn chảy dẻo xấp xỉ 248 MPa, tạo nên khả năng chịu tải và độ dai vượt trội. Tuy nhiên, khi kết cấu chịu tác động của tải trọng tĩnh lớn hoặc tải trọng ngang cực hạn, thép thể hiện rõ tính chất phi tuyến vật liệu phức tạp với các vùng biến dạng dẻo cục bộ. Phương pháp phân tích phần tử hữu hạn truyền thống dựa trên chuyển vị thường bộc lộ hai hạn chế lớn: đòi hỏi phải chia nhỏ cấu kiện thành 10 đến 20 phần tử để đạt độ hội tụ chính xác và chỉ áp dụng thuận lợi cho các tiết diện lăng trụ đồng nhất, gây lãng phí bộ nhớ đệm và kéo dài thời gian tính toán.

Nhằm giải quyết triệt để vấn đề này, luận văn thạc sĩ chuyên ngành Xây dựng công trình dân dụng và công nghiệp tại Trường Đại học Bách Khoa thuộc Đại học Quốc gia Thành phố Hồ Chí Minh (thực hiện từ tháng 7 năm 2012 đến tháng 6 năm 2013) đã tập trung nghiên cứu mô hình phân tích khớp dẻo cho khung thép phẳng bằng phương pháp lực. Mục tiêu trọng tâm của đề tài là xây dựng hoàn chỉnh ma trận độ mềm phần tử dầm-cột thông qua hàm nội suy lực chính xác, tích hợp ảnh hưởng của ứng suất dư theo tiêu chuẩn châu Âu và xét đến chiều dài khớp dẻo thực tế. Nghiên cứu ứng dụng kỹ thuật tích phân số Gauss-Radau cải tiến kết hợp giải thuật lặp phi tuyến Newton-Raphson trên môi trường lập trình MATLAB. Đóng góp nổi bật của luận văn là cho phép mô phỏng chính xác toàn bộ ứng xử uốn dẻo của một cấu kiện dầm hoặc cột chỉ bằng đúng 1 phần tử duy nhất, giúp tối ưu hóa hơn 50% thời gian xử lý và kiểm soát sai số tải trọng tới hạn dưới mức 1,5% so với các phần mềm thương mại chuyên dụng.

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

Khung lý thuyết áp dụng

Luận văn thiết lập nền tảng tính toán dựa trên hai trụ cột lý thuyết chính: lý thuyết phần tử hữu hạn dựa trên phương pháp lực và mô hình khớp dẻo thớ phân tích phi tuyến vật liệu. Khác với phương pháp chuyển vị giả định trường chuyển vị gần đúng, phương pháp lực vận dụng nguyên lý lực ảo để thiết lập quan hệ trực tiếp giữa lực hai đầu phần tử và nội lực mặt cắt thông qua ma trận hàm nội suy lực. Hàm nội suy này là nghiệm chính xác thỏa mãn phương trình vi phân cân bằng tĩnh học dọc theo chiều dài thanh, trong đó biểu đồ lực dọc có dạng hằng số và biểu đồ moment uốn biến thiên tuyến tính.

Để mô tả sự biến dạng dẻo từng phần, mô hình phân tử thớ được áp dụng cho tiết diện thép hình chữ I. Toàn bộ mặt cắt được chia thành 66 thớ thành phần, bao gồm 9 đoạn theo phương ngang và 3 lớp theo chiều dày cho mỗi bản cánh, kết hợp cùng 12 thớ dọc theo chiều cao bản bụng. Mỗi thớ tuân theo quan hệ ứng suất - biến dạng song tuyến tính có xét đến hệ số tái bền vật liệu và sự phân bố của ứng suất dư theo mô hình tiêu chuẩn ECCS.

Bên cạnh đó, khái niệm chiều dài khớp dẻo được định lượng thông qua hệ số chiều dài dao động từ 0,02 đến 0,08 lần chiều dài nhịp. Thay vì sử dụng tích phân điểm giữa Gauss-Legendre hay tích phân điểm đầu Gauss-Radau thông thường vốn gây sai số tại vị trí moment cực đại hoặc tập trung biến dạng giả tạo, tác giả áp dụng kỹ thuật tích phân hai điểm Gauss-Radau cải tiến. Kỹ thuật này mở rộng phạm vi tích phân lên 4 lần chiều dài khớp dẻo, đặt điểm lấy tích phân chính xác tại hai đầu mút cấu kiện, giúp nắm bắt trọn vẹn giá trị nội lực lớn nhất và phản ánh đúng hàm biến dạng bậc hai trong vùng chảy dẻo.

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

Nghiên cứu sử dụng nguồn dữ liệu kiểm chứng bao gồm 4 mô hình kết cấu thép kinh điển từ các công bố quốc tế uy tín: dầm công xôn tiết diện W21x50 có chiều dài 10,8 m chịu tải tập trung đầu tự do, cột thép liên kết một đầu ngàm một đầu khớp chịu uốn nén đồng thời, khung phẳng 1 tầng 2 nhịp và khung cổng Vogel 1 tầng chịu tải trọng phức hợp. Phương pháp chọn mẫu tập trung vào các dạng kết cấu chịu lực điển hình nhằm khảo sát đa dạng trạng thái ứng suất từ uốn thuần túy đến tương tác uốn - nén dọc trục phức tạp.

Lý do lựa chọn phương pháp phân tích lực kết hợp tích phân Gauss-Radau cải tiến xuất phát từ khả năng loại bỏ hoàn toàn hiện tượng khóa lưới và giảm thiểu tối đa số lượng bậc tự do của toàn hệ thống. Quy trình tính toán được lập trình tự động hóa trên nền tảng MATLAB với thuật giải lặp phi tuyến Newton-Raphson hai cấp độ: vòng lặp ngoài giải quyết cân bằng lực không cân bằng tại các nút kết cấu và vòng lặp con lặp qua từng thớ mặt cắt để cập nhật ma trận độ mềm tiếp tuyến. Toàn bộ chương trình tính toán và kiểm chứng được hoàn thành sau 12 tháng nghiên cứu chuyên sâu từ ngày 02 tháng 07 năm 2012 đến ngày 21 tháng 06 năm 2013.

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

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

Quá trình phân tích số trên 4 bài toán kiểm chứng đã mang lại các phát hiện cốt lõi:

Thứ nhất, tính chính xác vượt trội của mô hình 1 phần tử trên mỗi cấu kiện. Kết quả tính toán hệ số tải trọng tới hạn đối với dầm công xôn W21x50 và tải trọng giới hạn trên khung cổng Vogel cho thấy độ lệch số học chỉ dao động từ 0,8% đến 1,4% khi đối chiếu với mô hình chi tiết sử dụng hàng chục phần tử trong phần mềm ABAQUS và OpenSees.

Thứ hai, tác động rõ nét của ứng suất dư lên ứng xử chịu lực. Khi đưa mô hình ứng suất dư ECCS vào tính toán trên 66 thớ tiết diện dầm W21x50 dài 10,8 m, khả năng chịu tải tới hạn của dầm bị suy giảm từ 4,2% đến 7,8% so với trường hợp vật liệu lý tưởng không có ứng suất dư. Biến dạng dẻo xuất hiện sớm hơn tại các mép cánh dầm làm giảm nhanh độ cứng chống uốn tổng thể.

Thứ ba, sự ổn định của kỹ thuật tích phân Gauss-Radau cải tiến theo các hệ số chiều dài khớp dẻo. Khi khảo sát hệ số chiều dài khớp dẻo thay đổi ở các mức 0,02; 0,04 và 0,08, đường cong quan hệ lực - chuyển vị vẫn duy trì độ mượt mà, loại bỏ triệt để hiện tượng phân kỳ toán học và phản ánh sát thực tế vùng lan truyền dẻo ở các đầu ngàm liên kết.

Thứ tư, hiệu quả vượt bậc về tài nguyên tính toán. Việc giảm số lượng phần tử từ 10 xuống còn 1 phần tử cho mỗi cấu kiện giúp giảm hơn 65% dung lượng ma trận độ cứng tổng thể, đồng thời rút ngắn thời gian lặp hội tụ hơn 50% trong các bài toán phân tích phi tuyến khung nhiều tầng nhiều nhịp.

Thảo luận kết quả

Các phát hiện trên khẳng định tính ưu việt của việc xây dựng hàm nội suy lực dựa trên phương trình cân bằng nội lực chính xác thay vì trường chuyển vị ước lượng. Kết quả nghiên cứu có thể được biểu diễn một cách trực quan thông qua bảng so sánh tải trọng tới hạn và các đồ thị quan hệ lực - chuyển vị đỉnh dầm.

Trên đồ thị biểu diễn quan hệ hệ số tải và chuyển vị đầu mút, đường cong phản ứng kết cấu thể hiện rõ 3 giai đoạn làm việc: giai đoạn đàn hồi tuyến tính ban đầu, giai đoạn chuyển tiếp đàn - dẻo khi các thớ ngoài cùng của bản cánh bắt đầu chảy dẻo, và giai đoạn hình thành khớp dẻo hoàn toàn tại vùng ngàm dẫn đến độ võng tăng nhanh không kiểm soát. Khi so sánh với các công trình nghiên cứu của Scott và Fenves năm 2006 trên cấu kiện bê tông cốt thép, giải thuật trong luận văn đã mở rộng thành công cho kết cấu thép hình có xét đến tính dị hướng của ứng suất dư cán nóng, tạo nên sự tương thích cao với dữ liệu thực nghiệm quốc tế.

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

Nhằm phát huy giá trị học thuật và tính ứng dụng thực tiễn của đề tài, các giải pháp và lộ trình triển khai cụ thể được đề xuất như sau:

Thứ nhất, tích hợp thuật toán phân tích khớp dẻo dựa trên phương pháp lực vào các phần mềm tính toán kết cấu thương mại trong nước nhằm rút ngắn 40% thời gian phân tích phi tuyến tĩnh cho các dự án nhà cao tầng và nhà xưởng nhịp lớn trong giai đoạn 2026-2028, do các kỹ sư kết cấu và viện nghiên cứu chuyên ngành chủ trì thực hiện.

Thứ hai, chuẩn hóa quy trình chia lưới thớ mặt cắt với mật độ tối thiểu 66 thớ trên tiết diện chữ I và bắt buộc xét đến ứng suất dư ECCS trong công tác thẩm tra an toàn chịu lực công trình để nâng cao độ tin cậy kết cấu lên trên 98%, triển khai thực hiện trong vòng 6 tháng tới bởi các công ty tư vấn thiết kế và giám định xây dựng.

Thứ ba, mở rộng mô hình toán học từ khung phẳng 2D sang hệ khung không gian 3D chịu xoắn và phân tích phi tuyến động lực học dưới tác động của động đất, hướng đến mục tiêu hoàn thiện gói thuật toán mở rộng trong thời gian 12 đến 24 tháng do các nhóm nghiên cứu học thuật tại các trường đại học khối kỹ thuật phụ trách.

Thứ tư, biên soạn tài liệu hướng dẫn và chuyển giao mã nguồn chương trình MATLAB vào chương trình đào tạo sau đại học ngành Kỹ thuật Xây dựng để gia tăng 30% tỷ lệ ứng dụng tính toán phi tuyến trong các đề tài thạc sĩ và kỹ sư từ năm học 2026-2027, dưới sự chỉ đạo của các hội đồng khoa học chuyên ngành.

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

Luận văn mang lại giá trị chuyên môn sâu sắc cho 4 nhóm đối tượng chính:

Nhóm thứ nhất là các kỹ sư thiết kế kết cấu công trình xây dựng dân dụng và công nghiệp. Luận văn cung cấp phương pháp đánh giá chính xác khả năng chịu tải sau đàn hồi của khung thép, giúp tối ưu hóa kích thước tiết diện dầm cột và tiết kiệm từ 5% đến 10% khối lượng thép sử dụng trong các công trình nhà công nghiệp nhịp lớn.

Nhóm thứ hai là các học viên cao học, nghiên cứu sinh và nhà khoa học chuyên ngành Xây dựng. Tài liệu này đóng vai trò là tài liệu tham khảo nền tảng về lý thuyết phương pháp lực, tích phân số Gauss-Radau cải tiến và thuật giải Newton-Raphson, hỗ trợ nghiên cứu mở rộng sang các bài toán phi tuyến hình học bậc cao hoặc kết cấu liên hợp thép - bê tông.

Nhóm thứ ba là các chuyên gia phát triển phần mềm mô phỏng kỹ thuật và giải pháp tính toán kết cấu. Cấu trúc chương trình MATLAB chi tiết trong luận văn cung cấp mã nguồn mẫu để xây dựng các module giải phi tuyến vật liệu hiệu năng cao cho các ứng dụng tính toán thương mại.

Nhóm thứ tư là giảng viên và sinh viên các trường đại học kỹ thuật công trình. Luận văn là học liệu giá trị phục vụ giảng dạy các môn học chuyên sâu như Cơ học kết cấu nâng cao, Phương pháp phần tử hữu hạn và Phân tích phi tuyến công trình.

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

Phương pháp lực có ưu thế vượt trội nào so với phương pháp chuyển vị truyền thống khi phân tích khung thép? Phương pháp lực sử dụng hàm nội suy lực thỏa mãn chính xác các phương trình cân bằng tĩnh học dọc theo phần tử. Nhờ đó, mô hình chỉ cần sử dụng đúng 1 phần tử cho mỗi cấu kiện dầm hoặc cột thay vì phải chia nhỏ từ 10 đến 20 phần tử như phương pháp chuyển vị, giúp giảm hơn 50% thời gian tính toán mà vẫn đảm bảo sai số dưới 1,5%.

Mô hình phần tử thớ được thiết lập và tính toán như thế nào trên tiết diện dầm thép chữ I? Tiết diện dầm chữ I được chia thành 66 thớ nhỏ, gồm 9 thớ theo chiều rộng và 3 thớ theo chiều dày cho mỗi bản cánh cùng 12 thớ ở bản bụng. Mỗi thớ được gán đặc tính ứng suất - biến dạng song tuyến tính và ứng suất dư riêng biệt, cho phép theo dõi chính xác quá trình chảy dẻo cục bộ tại từng tọa độ mặt cắt.

Vì sao kỹ thuật tích phân Gauss-Radau cải tiến lại được lựa chọn thay vì các phương pháp tích phân khác? Gauss-Radau cải tiến đặt điểm lấy tích phân chính xác tại hai đầu mút cấu kiện nơi có moment uốn cực đại, đồng thời mở rộng miền tích phân lên 4 lần chiều dài khớp dẻo. Phương pháp này chỉ cần 2 điểm tính toán cho mỗi khớp dẻo mà vẫn mô tả chính xác hàm biến dạng bậc hai, triệt tiêu hiện tượng tập trung biến dạng giả tạo.

Ứng suất dư ảnh hưởng như thế nào đến khả năng chịu lực cực hạn của dầm thép? Theo mô hình ứng suất dư ECCS, ứng suất dư làm các thớ mép cánh đạt đến giới hạn chảy dẻo 248 MPa sớm hơn bình thường. Điều này dẫn đến sự suy giảm độ cứng chống uốn sớm, làm giảm từ 4,2% đến 7,8% tải trọng tới hạn của dầm thép W21x50 dài 10,8 m so với trường hợp không xét ứng suất dư.

Độ tin cậy của chương trình tính toán viết bằng MATLAB đã được kiểm chứng qua các mô hình nào? Độ chính xác của chương trình đã được kiểm chứng qua 4 bài toán kết cấu kinh điển gồm dầm công xôn W21x50, cột ngàm khớp uốn nén, khung phẳng 1 tầng 2 nhịp và khung cổng Vogel 1 tầng. Mọi kết quả tải trọng tới hạn đều trùng khớp chặt chẽ với dữ liệu từ các phần mềm chuyên dụng như OpenSees và ABAQUS.

Kết luận

Luận văn đã giải quyết trọn vẹn và chuẩn xác bài toán phân tích phi tuyến vật liệu cho khung thép phẳng với những đóng góp khoa học nổi bật:

  • Thiết lập thành công mô hình phần tử dầm-cột dựa trên phương pháp lực, tích hợp trọn vẹn ảnh hưởng của ứng suất dư ECCS và chiều dài khớp dẻo thực tế.
  • Ứng dụng xuất sắc kỹ thuật tích phân hai điểm Gauss-Radau cải tiến, cho phép mô phỏng chính xác ứng xử chảy dẻo của toàn bộ cấu kiện chỉ với 1 phần tử duy nhất.
  • Lập trình hoàn chỉnh hệ thống tính toán trên ngôn ngữ MATLAB sử dụng thuật giải lặp phi tuyến Newton-Raphson hai cấp độ đạt độ ổn định và hội tụ cao.
  • Kiểm chứng độ tin cậy tuyệt đối qua 4 bài toán kết cấu chuẩn mực với sai số tải trọng tới hạn kiểm soát dưới mức 1,5% so với OpenSees và ABAQUS.
  • Tối ưu hóa hiệu năng tính toán vượt bậc thông qua việc giảm hơn 65% dung lượng bộ nhớ ma trận và tiết kiệm trên 50% thời gian xử lý dữ liệu.

Định hướng tiếp theo của nghiên cứu là tiếp tục phát triển mô hình cho kết cấu khung không gian 3D và bài toán phân tích động lực học phi tuyến chịu động đất trong giai đoạn 2026-2027. Các kỹ sư và nhà nghiên cứu kết cấu được khuyến khích khai thác ngay phương pháp lực và mã nguồn thuật toán để nâng cao độ chính xác và hiệu quả kinh tế cho các dự án xây dựng hiện đại.