Tổng quan nghiên cứu

Trong các sự kiện địa chấn lịch sử như trận động đất Northridge năm 1994 hay trận động đất El Centro năm 1940, ứng xử thực tế của kết cấu khung thép bộc lộ rõ tính phi tuyến phức tạp do sự suy giảm độ cứng và tiêu tán năng lượng qua biến dạng dẻo. Tuy nhiên, các phương pháp thiết kế truyền thống theo tiêu chuẩn đàn hồi tuyến tính bậc một thường bỏ qua sự phân bố dẻo dọc chiều dài cấu kiện và phần năng lượng dỡ tải, dẫn đến việc ước tính sai lệch khả năng chịu tải tới hạn từ 15% đến 25%. Vấn đề nghiên cứu trọng tâm của luận văn là phân tích phi tuyến toàn diện khung thép phẳng chịu tải trọng động và tải trọng lặp, có xét đến sự đóng góp của năng lượng dỡ tải nhằm nâng cao độ an toàn cho công trình.

Mục tiêu cụ thể của đề tài là xây dựng một phần tử hữu hạn dầm-cột mới thông qua nguyên lý thế năng toàn phần dừng, đồng thời phát triển chương trình máy tính chuyên dụng có tên SDNASF trên môi trường ngôn ngữ C++ và MATLAB. Nghiên cứu được triển khai trong phạm vi thời gian 10 tháng, từ ngày 18 tháng 08 năm 2014 đến ngày 15 tháng 06 năm 2015 tại Trường Đại học Bách Khoa thuộc Đại học Quốc gia Thành phố Hồ Chí Minh. Ý nghĩa then chốt của công trình thể hiện qua việc cung cấp công cụ phân tích nâng cao có độ chính xác cao, giúp giảm thiểu 10% chi phí vật liệu trong thiết kế trực tiếp theo tiêu chuẩn Eurocode 3 và AISC, đồng thời kiểm soát độ trôi dạt tầng công trình kháng chấn với sai số dưới 3%.

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 lý thuyết thế năng toàn phần dừng kết hợp với phương pháp biến phân Rayleigh-Ritz để thiết lập ma trận độ cứng cát tuyến và tiếp tuyến cho phần tử dầm-cột phẳng gồm 6 bậc tự do. Mô hình vật liệu được xây dựng theo quy luật đàn - dẻo - tái bền tuyến tính, phản ánh chân thực hiệu ứng Bauschinger và quá trình dỡ tải theo mô đun đàn hồi gốc. Ba khái niệm cốt lõi được tích hợp chặt chẽ trong khung lý thuyết gồm: sự dịch chuyển trục trung hòa của lõi đàn hồi khi mặt cắt bị chảy dẻo từng phần, mô hình ứng suất dư chuẩn hóa theo Lehigh và Vogel, cùng ứng xử trễ của liên kết nửa cứng. Để mô phỏng độ cứng xoay phi tuyến tại các nút dầm-cột, luận văn áp dụng mô hình ba thông số Kishi-Chen, mô hình bốn thông số Richard-Abbott và mô hình hàm mũ Chen-Lui với các tham số độ cứng ban đầu và mô men giới hạn.

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

Nguồn dữ liệu kiểm chứng của luận văn được thu thập từ bộ mẫu gồm 8 mô hình kết cấu kinh điển đã được công bố trên các tạp chí khoa học quốc tế uy tín, bao gồm 4 bài toán tĩnh và 4 bài toán động lực học phức tạp. Phương pháp chọn mẫu là chọn mẫu có chủ đích dựa trên các tiêu chuẩn thực nghiệm chuẩn mực, tiêu biểu như khung cổng Vogel với độ xiên ban đầu 1/400, khung Chan-Chui 2 tầng 1 nhịp và khung thép chịu chuỗi dữ liệu gia tốc nền của trận động đất El Centro năm 1940.

Lý do lựa chọn phương pháp phân tích vùng dẻo thay vì phương pháp khớp dẻo tập trung là vì phương pháp này cho phép chia nhỏ mặt cắt ngang thành nhiều thớ vi phân và chia cấu kiện thành nhiều phân đoạn dọc trục. Nhờ đó, chương trình mô phỏng chính xác sự lan truyền dẻo theo cả hai phương mà không cần sử dụng các hàm giả định đơn giản hóa. Hệ phương trình động lực học phi tuyến được thiết lập dựa trên ma trận khối lượng tương thích, ma trận cản Rayleigh và được giải bằng phương pháp tích phân từng bước Newmark kết hợp với thuật toán lặp Newton-Raphson hiệu chỉnh cùng kỹ thuật chiều dài cung. Toàn bộ quy trình mô ph phỏng số được lập trình tự động hóa hoàn toàn trong suốt tiến độ nghiên cứu từ tháng 08 năm 2014 đến tháng 06 năm 2015.

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

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

Thứ nhất, việc bổ sung thành phần năng lượng dỡ tải vào ma trận độ cứng phần tử đã nâng cao độ chính xác khi dự báo quan hệ tải trọng và chuyển vị, giúp hệ số tải giới hạn của khung cổng Vogel đạt mức hội tụ cao với độ sai lệch dưới 2.5% so với kết quả giải tích của Vogel công bố năm 1985.

Thứ hai, trong bài toán phân tích khung 2 tầng chịu kích động địa chấn El Centro năm 1940, khi giới hạn chảy của thép tăng từ 235 MPa lên 260 MPa, biên độ chuyển vị ngang tại nút đỉnh giảm khoảng 18%, đồng thời làm thay đổi rõ rệt hình dáng và diện tích của vòng lặp trễ tiêu tán năng lượng tại mép cánh cấu kiện.

Thứ ba, sự hiện diện của liên kết nửa cứng làm chu kỳ dao động riêng của hệ kết cấu tăng từ 10% đến 25% ở các mode cơ bản, giúp giảm lực cắt đáy công trình tới 15% so với khung có nút liên kết cứng tuyệt đối.

Thứ tư, khi kiểm chứng trên mô hình khung nửa cứng 2 nhịp 6 tầng chịu tải trọng động tập trung 10.44 kN tác dụng trong thời gian 11 giây, chương trình SDNASF cho kết quả phản ứng tức thời về chuyển vị đỉnh trùng khớp hoàn hảo với các nghiên cứu động lực học phi tuyến quốc tế.

Thảo luận kết quả

Khả năng hấp thụ năng lượng của khung thép chịu tải trọng động và tải trọng lặp bắt nguồn từ hai cơ chế chính: năng lượng biến dạng dẻo của vật liệu thép và năng lượng tiêu tán qua vòng trễ ma sát của các liên kết dầm-cột nửa cứng. Khi trục trung hòa của lõi đàn hồi bị dịch chuyển, độ cứng tiếp tuyến của toàn hệ thống suy giảm liên tục theo mức độ dẻo hóa trên từng thớ tiết diện.

Những phản ứng này được thể hiện trực quan thông qua biểu đồ quan hệ tải trọng - chuyển vị đỉnh, biểu đồ vòng lặp trễ mô men - góc xoay khép kín và các bảng so sánh sai số định lượng giữa các phương pháp. Kết quả đối sánh trên các bảng dữ liệu chỉ ra rằng việc kết hợp phương pháp Newmark với thuật toán Newton-Raphson hiệu chỉnh và năng lượng dỡ tải giúp loại bỏ hoàn toàn hiện tượng mất cân bằng lực nội tại, đảm bảo tốc độ hội tụ nhanh hơn 20% so với phương pháp khớp dẻo hiệu chỉnh truyền thống.

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

Thứ nhất, triển khai ứng dụng công cụ phần mềm SDNASF vào quy trình phân tích và thiết kế trực tiếp các dự án nhà cao tầng kết cấu thép tại Việt Nam, nhằm mục tiêu tối ưu hóa từ 8% đến 12% trọng lượng thép tiêu thụ và giảm 30% thời gian chạy mô phỏng phi tuyến trong vòng 12 tháng tới. Đơn vị chủ trì thực hiện là các kỹ sư kết cấu tại các công ty tư vấn thiết kế công trình xây dựng.

Thứ hai, đưa mô hình liên kết nửa cứng ba thông số Kishi-Chen và bốn thông số Richard-Abbott vào hệ thống tiêu chuẩn thiết kế công trình kháng chấn quốc gia trong lộ trình 2025-2030, do các viện nghiên cứu chuyên ngành và Hội đồng Tiêu chuẩn Xây dựng Việt Nam phụ trách.

Thứ ba, tổ chức các khóa tập huấn chuyên môn định kỳ 6 tháng một lần về giải thuật phân tích vùng dẻo và phương pháp tích phân số Newmark cho hơn 200 kỹ sư kết cấu, nâng cao năng lực ứng dụng phân tích phi tuyến trong thực tế.

Thứ tư, mở rộng và phát triển mã nguồn chương trình sang không gian 3 chiều có xét đến hiệu ứng mất ổn định ngoài mặt phẳng và xoắn phi tuyến của cấu kiện thép thành mỏng trong giai đoạn nghiên cứu 2 năm tiếp theo, do nhóm tác giả phối hợp cùng các nghiên cứu sinh chuyên ngành Công trình thực hiện.

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

Nhóm thứ nhất là kỹ sư thiết kế kết cấu công trình dân dụng và công nghiệp, những người cần công cụ phân tích nâng cao để tính toán khung thép kháng chấn với độ chính xác đạt trên 95%, đảm bảo tính tối ưu kinh tế cho các dự án nhà xưởng vượt nhịp lớn và nhà cao tầng.

Nhóm thứ hai là giảng viên, học viên cao học và nghiên cứu sinh chuyên ngành Xây dựng công trình, có thể sử dụng khung lý thuyết thế năng toàn phần dừng và thuật toán phần tử dầm-cột 6 bậc tự do làm tài liệu giảng dạy và nền tảng phát triển các đề tài nghiên cứu chuyên sâu về cơ học kết cấu phi tuyến.

Nhóm thứ ba là các nhà phát triển phần mềm mô phỏng kỹ thuật xây dựng, tiếp cận trực tiếp cấu trúc thuật toán C++ và MATLAB gồm hơn 10 khối module tính toán để tích hợp tính năng phân tích liên kết nửa cứng vào các phần mềm thương mại.

Nhóm thứ tư là các chuyên gia thẩm tra và ban quản lý dự án công trình cấp 1 và cấp đặc biệt, áp dụng các tiêu chí đánh giá tiêu tán năng lượng vòng trễ để nghiệm thu giải pháp kết cấu an toàn tại các vùng có nguy cơ động đất từ cấp 7 đến cấp 9.

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

Câu hỏi 1: Việc xét đến năng lượng dỡ tải trong ma trận độ cứng mang lại lợi ích gì?

Trả lời: Khi công trình chịu tải trọng động đất hoặc tải trọng gió giật, vật liệu thép liên tục trải qua các chu kỳ chất tải và dỡ tải. Tích hợp năng lượng dỡ tải giúp ma trận độ cứng phản ánh chính xác mô đun đàn hồi hồi phục, loại bỏ sai số tích lũy năng lượng và kiểm soát độ lệch chuyển vị dưới 3%.

Câu hỏi 2: Phương pháp vùng dẻo vượt trội hơn phương pháp khớp dẻo ở điểm nào?

Trả lời: Phương pháp vùng dẻo chia cắt tiết diện thành nhiều thớ vi phân và phân đoạn dọc trục, cho phép theo dõi chính xác sự lan truyền dẻo và dịch chuyển trục trung hòa lõi đàn hồi. Trong khi đó, phương pháp khớp dẻo chỉ giả định dẻo cục bộ tại hai đầu mút, dễ gây sai lệch khả năng chịu tải tới 20%.

Câu hỏi 3: Liên kết nửa cứng ảnh hưởng thế nào đến khả năng kháng chấn của khung thép?

Trả lời: Liên kết nửa cứng hoạt động như cơ chế giảm chấn vòng trễ tự nhiên, giúp hấp thụ và tiêu tán phần lớn năng lượng sóng chấn động. Dữ liệu từ trận động đất Northridge năm 1994 chứng minh khung thép có liên kết mềm giảm được 15% lực cắt đáy tác dụng lên móng công trình.

Câu hỏi 4: Thuật toán Newmark kết hợp Newton-Raphson xử lý bài toán phi tuyến ra sao?

Trả lời: Phương pháp Newmark giải quyết phương trình vi phân chuyển động theo từng bước thời gian, trong khi thuật toán Newton-Raphson hiệu chỉnh liên tục cập nhật ma trận độ cứng tiếp tuyến để triệt tiêu lực mất cân bằng. Sự kết hợp này đảm bảo hội tụ ổn định ngay cả khi tiết diện dẻo hóa hơn 50%.

Câu hỏi 5: Chương trình SDNASF đã được kiểm chứng qua các mô hình chuẩn nào?

Trả lời: Phần mềm SDNASF đã được đối chứng tin cậy qua 8 bài toán chuẩn quốc tế, tiêu biểu là khung Vogel 6 tầng 2 nhịp năm 1985 và khung 2 tầng chịu kích động động đất El Centro năm 1940, với đường cong phản ứng hoàn toàn trùng khớp với các công bố thực nghiệm.

Kết luận

  • Xây dựng thành công phần tử hữu hạn dầm-cột phi đàn hồi mới dựa trên nguyên lý thế năng toàn phần dừng có xét đến năng lượng dỡ tải.
  • Tích hợp đồng thời các yếu tố phi tuyến hình học, phi tuyến vật liệu với sự dịch chuyển trục trung hòa của lõi đàn hồi và ứng suất dư chuẩn Lehigh.
  • Mô phỏng chính xác ứng xử vòng trễ phi tuyến của các dạng liên kết nửa cứng dưới tác động của tải trọng lặp và tải trọng động lực học.
  • Thiết lập giải thuật động lực học Newmark kết hợp kỹ thuật lặp Newton-Raphson hiệu chỉnh và phương pháp chiều dài cung đạt tốc độ hội tụ cao.
  • Phát triển hoàn chỉnh chương trình máy tính SDNASF trên nền tảng C++ và MATLAB, kiểm chứng thành công trên 8 mô hình kết cấu kinh điển.

Đóng góp lớn nhất của luận văn là hoàn thiện cơ sở lý thuyết và công cụ tính toán số tin cậy cho phương pháp thiết kế trực tiếp kết cấu thép chịu tải trọng động đất. Kế hoạch tiếp theo trong giai đoạn 2025-2027 là nâng cấp thuật toán cho kết cấu không gian 3 chiều. Các đơn vị tư vấn và kỹ sư công trình hãy ứng dụng ngay giải pháp phân tích nâng cao này vào thực tiễn để nâng cao chất lượng và độ an toàn bền vững cho mọi công trình.