Tổng quan nghiên cứu

Trong kỹ thuật xây dựng hiện đại, kết cấu khung thép nhiều tầng ngày càng được ứng dụng rộng rãi trong các công trình dân dụng và công nghiệp quy mô lớn nhờ khả năng vượt nhịp tiêu chuẩn 9,14 m và tiến độ thi công vượt trội. Tuy nhiên, sau các biến cố địa chấn nghiêm trọng như trận động đất Northridge năm 1994, hiện tượng phá hoại giòn và biến dạng trượt tại các liên kết dầm cột đã chứng minh rằng giả thiết nút khung cứng tuyệt đối trong thiết kế truyền thống không còn phản ánh chính xác sự làm việc thực tế. Việc bỏ qua độ mềm của vùng mắt nút (panel zone) có thể dẫn đến việc đánh giá sai lệch đáng kể về độ cứng tổng thể, khả năng phân tán năng lượng và cơ chế sụp đổ của toàn bộ công trình.

Luận văn tập trung khảo sát chuyên sâu ứng xử của kết cấu khung thép phẳng 16 tầng 4 nhịp thông qua phương pháp phân tích tĩnh phi tuyến đẩy dần (Pushover Analysis). Mục tiêu trọng tâm của nghiên cứu là đánh giá mức độ ảnh hưởng của biến dạng nút khung và sự hình thành khớp dẻo trong dầm đến khả năng chịu tải trọng ngang của công trình. Phạm vi nghiên cứu được thiết lập cho kết cấu khung thép 16 tầng có tổng chiều cao 59,47 m (chiều cao tầng 1 là 4,57 m, các tầng còn lại cao 3,66 m), giả định xây dựng tại khu vực thành phố Pleiku với đỉnh gia tốc nền tham chiếu aGR = 0,0511 g trên nền đất loại C. Ý nghĩa của công trình thể hiện qua việc định lượng các chỉ số an toàn quan trọng như chu kỳ dao động riêng ban đầu T = 0,67 s, sự suy giảm lực cắt đáy và quy luật biến dạng dẻo, mang lại cơ sở khoa học tin cậy cho công tác thiết kế kháng chấn theo TCVN 9386:2012 và tiêu chuẩn quốc tế FEMA 356.

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 vận dụng đồng thời các lý thuyết cơ học kết cấu phi tuyến và mô hình liên kết dầm - cột bán cứng tiên tiến nhằm tái hiện chính xác ứng xử của vật liệu thép khi chịu tải trọng động đất:

  • Mô hình nút khung Gupta và Krawinkler: Khác với mô hình Scissors đơn giản hay mô hình Krawinkler (2000), mô hình Gupta và Krawinkler mô phỏng đầy đủ kích thước hình học vùng mắt nút thông qua hệ thống các thanh cứng biên kết hợp hai lò xo phi tuyến. Mô hình này mô tả chính xác ứng xử cắt của bản bụng cột và sự đóng góp của cánh cột với cường độ cắt dẻo Vy, diện tích cắt hiệu dụng Aeff và sức kháng cắt dẻo Vp.
  • Mô hình suy giảm Modified Ibarra Krawinkler Deterioration: Mô hình sử dụng đường xương sống (backbone curve) gồm ba giai đoạn: giai đoạn đàn hồi ban đầu Ke, giai đoạn chảy dẻo với cường độ Fy và giai đoạn tái bền Ks = 0,03Ke, kết hợp nhánh mềm hóa sau đỉnh để phản ánh sự suy giảm độ cứng và cường độ khi vật liệu chịu biến dạng dẻo lớn.
  • Các khái niệm cốt lõi: Phương pháp phân tích tĩnh phi tuyến đẩy dần (Pushover Analysis), đường cong khả năng (Capacity Curve), biến dạng góc trượt vùng mắt nút (Panel Zone Deformation), và ba mức độ giới hạn hiệu năng kết cấu theo tiêu chuẩn địa chấn gồm mức hoạt động tức thời (IO), mức an toàn sinh mạng (LS) và mức phòng ngừa sụp đổ (CP).

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

Nghiên cứu sử dụng mô hình kết cấu khung thép phẳng 16 tầng 4 nhịp với tổng chiều dài 36,56 m và tổng khối lượng công trình là 263,22 tấn. Các cột biên và cột giữa được chọn lọc từ chuỗi tiết diện W14 (từ W14x120 đến W14x550), trong khi hệ dầm sàn sử dụng các tiết diện chữ I lớn từ W30x116 đến W36x160 với bước nhịp 9,14 m.

Phương pháp phân tích được thực hiện bằng phần mềm mô phỏng phần tử hữu hạn mã nguồn mở OpenSees, đồng thời xây dựng một mô hình kiểm chứng độc lập trên phần mềm SAP2000 v19. Việc lựa chọn công cụ OpenSees xuất phát từ khả năng tùy biến sâu các phần tử lò xo phi tuyến và thuật toán lặp hội tụ mạnh mẽ. Quy trình phân tích đẩy dần được kiểm soát thông qua chuyển vị đỉnh mái dưới tác dụng của hệ tải trọng ngang phân bố dọc theo chiều cao nhà, tính toán song song theo hai tiêu chuẩn: TCVN 9386:2012 (với hệ số gia tốc phổ Sd(T1) = 0,405) và FEMA 356 (với hệ số phân phối k nội suy theo chu kỳ T = 0,67 s và hệ số khối lượng hữu hiệu Cm = 0,9). Toàn bộ dữ liệu phản ứng động lực học và quá trình phân phối nội lực được trích xuất để so sánh đối chiếu giữa 3 trường hợp: nút cứng tuyệt đối, nút biến dạng kết hợp dầm đàn hồi, và nút biến dạng kết hợp dầm khớp dẻo.

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 và phân tích chuyên sâu trên hệ khung thép 16 tầng mang lại các phát hiện kỹ thuật mang tính thực tiễn cao:

  1. Gia tăng chu kỳ dao động và độ mềm tổng thể: Khi xét đến biến dạng vùng mắt nút, chu kỳ dao động riêng ở dạng cơ bản tăng từ T = 0,67 s lên khoảng 0,78 s (tăng 16,4%). Sự linh hoạt của vùng nút làm giảm độ cứng ban đầu của toàn bộ khung thép 16 tầng.
  2. Suy giảm sức kháng cắt đáy và khả năng chịu lực cực hạn: Sức kháng cắt đáy danh nghĩa của kết cấu theo mô hình nút biến dạng giảm khoảng 8,5% đến 11,2% so với mô hình nút cứng truyền thống khi tính toán theo sơ đồ lực của FEMA 356 và TCVN 9386:2012.
  3. Tăng độ trôi tầng nguy hiểm tại các tầng trung gian: Độ trôi tầng cực đại khi có xét đến biến dạng nút tăng từ 18% đến 22% so với trường hợp nút cứng, tập trung chủ yếu từ tầng 7 đến tầng 11, nơi có sự chuyển tiếp giữa các tiết diện cột W14x370 và W14x311.
  4. Tái phân phối và chuyển dịch vị trí khớp dẻo: Khi nút khung bị biến dạng, ứng suất được phân tán đều hơn, làm chậm quá trình hình thành khớp dẻo nguy hiểm tại chân cột tầng 1 (tiết diện W14x550) và tập trung cơ chế tiêu tán năng lượng dẻo vào các liên kết dầm W36x160 ở 6 tầng dưới.
  5. Độ tin cậy cao giữa hai nền tảng tính toán: Tần số dao động riêng ở 3 dạng dao động đầu tiên giữa OpenSees và SAP2000 v19 có mức sai lệch rất nhỏ, chỉ dưới 2,3%, khẳng định tính chính xác của mô hình phần tử hữu hạn xây dựng trong OpenSees.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự suy giảm độ cứng và gia tăng độ trôi tầng là do bản bụng của cột thép W14 phải chịu ứng suất cắt tập trung rất lớn từ các mô men uốn ngược chiều của dầm W36 truyền vào. Biến dạng góc trượt của bản bụng nút khung xảy ra trước khi toàn bộ tiết diện dầm đạt tới trạng thái chảy dẻo hoàn toàn, tạo nên cơ chế liên kết nửa cứng.

Khi so sánh với các công bố khoa học của Krawinkler và các nghiên cứu về kết cấu khung thép chịu mô men (SMRF), kết quả hoàn toàn tương đồng về xu hướng làm mềm của hệ kết cấu. Dữ liệu phân tích có thể được mô tả trực quan qua biểu đồ đường cong khả năng (quan hệ giữa lực cắt đáy Vb và chuyển vị đỉnh mái) và bảng tổng hợp độ trôi 16 tầng. Biểu đồ đường cong khả năng chỉ rõ: đường đặc tuyến của mô hình nút cứng có độ dốc ban đầu lớn nhưng nhanh chóng đạt trạng thái suy giảm, trong khi mô hình nút biến dạng theo Gupta và Krawinkler thể hiện đường cong chuyển tiếp mượt hơn, phản ánh độ dẻo kết cấu tốt hơn trước khi chạm ngưỡng sụp đổ (CP). Ngoài ra, sơ đồ tải trọng ngang theo tiêu chuẩn FEMA 356 cho thấy lực đẩy tác dụng lên các tầng trên đỉnh lớn hơn khoảng 6,5% so với phân bố tam giác của TCVN 9386:2012, mang lại đánh giá an toàn và thận trọng hơn cho công trình cao tầng.

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

Dựa trên các kết quả mô phỏng định lượng từ luận văn, các giải pháp và khuyến nghị kỹ thuật được đề xuất như sau:

  • Tích hợp bắt buộc mô hình biến dạng mắt nút trong tính toán công trình cao tầng: Các đơn vị tư vấn thiết kế kết cấu cần áp dụng mô hình nút Gupta và Krawinkler khi thiết kế kháng chấn cho nhà khung thép từ 10 tầng trở lên. Mục tiêu là kiểm soát chính xác độ trôi tầng thực tế không vượt quá giới hạn an toàn 1,5% đến 2,0% chiều cao tầng theo quy định địa chấn.
  • Gia cường cục bộ bản bụng cột thép chịu lực cắt lớn: Đối với các khung thép sử dụng tiết diện cột W14 chịu tải trọng động đất lớn từ tầng 1 đến tầng 5, chủ đầu tư và kỹ sư thiết kế cần bổ sung sườn gia cường chéo hoặc tấm ốp tăng chiều dày bản bụng (doubler plates). Giải pháp này giúp giảm thiểu 25% đến 30% biến dạng trượt góc tại mắt nút, thực hiện trực tiếp trong giai đoạn thiết kế bản vẽ thi công.
  • Cập nhật quy chuẩn thiết kế kháng chấn quốc gia: Đề xuất Viện Khoa học Công nghệ Xây dựng và các cơ quan ban hành tiêu chuẩn xây dựng bổ sung hướng dẫn chi tiết về phân tích tĩnh phi tuyến đẩy dần có xét đến biến dạng phi tuyến của nút khung trong lộ trình cập nhật TCVN 9386 giai đoạn 2026-2028, tích hợp đầy đủ các hệ số hiệu chỉnh chuyển vị C1, C2, C3 tương thích với tiêu chuẩn quốc tế.
  • Ứng dụng công cụ mô phỏng mã nguồn mở chuyên sâu: Các doanh nghiệp tư vấn và viện nghiên cứu chuyên ngành nên đào tạo đội ngũ kỹ sư sử dụng thành thạo phần mềm OpenSees nhằm phân tích phi tuyến chuyên sâu, kết hợp các công cụ tiền và hậu xử lý trực quan để tối ưu hóa khối lượng kết cấu thép, giảm từ 5% đến 8% chi phí vật liệu mà vẫn bảo đảm hệ số an toàn chịu lực.

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

Nguồn tài liệu và số liệu thực nghiệm số trong luận văn mang lại giá trị ứng dụng thiết thực cho 4 nhóm đối tượng chính:

  • Kỹ sư kết cấu và chuyên gia thiết kế công trình: Sử dụng các bảng thông số tiết diện cột W14, dầm W36 và quy trình mô hình hóa liên kết nửa cứng để tính toán chính xác ứng xử kháng chấn cho các dự án nhà cao tầng và nhà xưởng công nghiệp nhịp lớn 9,14 m.
  • Học viên cao học và nghiên cứu sinh ngành Xây dựng: Tham khảo phương pháp xây dựng thuật toán phân tích tĩnh phi tuyến đẩy dần trên phần mềm OpenSees, cách khai báo các phần tử lò xo xoay và quy trình đối soát kết quả với SAP2000 v19 phục vụ các đề tài nghiên cứu chuyên sâu.
  • Giảng viên các trường đại học khối kỹ thuật: Trích xuất các dữ liệu tính toán mẫu của khung thép 16 tầng 4 nhịp với tổng khối lượng 263,22 tấn làm bài giảng chuyên đề cho các môn học Kết cấu thép nâng cao, Động lực học công trình và Thiết kế công trình chịu tải trọng động đất.
  • Cơ quan thẩm tra và ban quản lý dự án xây dựng: Sử dụng kết quả phân tích độ trôi tầng và cơ chế hình thành khớp dẻo làm tiêu chí đối chứng khi thẩm tra an toàn kết cấu cho các công trình thép cấp I và cấp đặc biệt tại các vùng có nguy cơ địa chấn.

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

Tại sao cần phải xét đến biến dạng vùng mắt nút khi phân tích khung thép cao tầng?

Trong thực tế, vùng mắt nút dầm cột phải chịu lực cắt rất lớn khi công trình chịu tải trọng ngang. Nếu giả định nút cứng tuyệt đối, độ cứng của công trình sẽ bị đánh giá cao hơn thực tế. Đối với khung 16 tầng, biến dạng nút làm tăng độ trôi tầng thêm khoảng 18% và làm suy giảm khả năng chịu lực cực hạn, do đó việc mô hình hóa nút giúp dự báo chính xác nguy cơ mất an toàn.

Mô hình nút khung Gupta và Krawinkler có ưu điểm gì nổi bật so với mô hình Scissors?

Mô hình Scissors xem toàn bộ vùng nút là một khối cứng liên kết với một lò xo xoay tập trung, bỏ qua tương tác cắt giữa cánh và bụng cột nên dẫn đến sai lệch lớn. Mô hình Gupta và Krawinkler mô phỏng đầy đủ kích thước hình học, phân tách rõ ứng xử biến dạng cắt của bản bụng và sự đóng góp của cánh cột, mang lại kết quả mô phỏng sát với thực tế làm việc của liên kết dầm cột W14 và W36.

Phần mềm OpenSees có thế mạnh gì so với SAP2000 v19 trong phân tích phi tuyến?

SAP2000 v19 sở hữu giao diện đồ họa trực quan thuận tiện nhưng bị giới hạn ở một số mô hình vật liệu suy giảm phức tạp. Trong khi đó, OpenSees là nền tảng mã nguồn mở chuyên dụng cho kỹ thuật địa chấn, cung cấp hơn 220 loại vật liệu và 160 loại phần tử, cho phép mô phỏng chính xác các lò xo xoay phi tuyến theo mô hình Modified Ibarra Krawinkler với độ lệch tần số dao động so với SAP2000 dưới 2,3%.

Tiêu chuẩn TCVN 9386:2012 và FEMA 356 khác nhau như thế nào trong phân bố tải trọng ngang?

TCVN 9386:2012 sử dụng phổ gia tốc đàn hồi với gia tốc nền thiết kế xác định theo tọa độ địa lý (như aGR = 0,0511 g tại Pleiku), phân phối lực ngang tăng tuyến tính theo chiều cao. Ngược lại, FEMA 356 sử dụng hệ số phân phối lũy thừa k nội suy từ 1,0 đến 2,0 theo chu kỳ dao động T = 0,67 s, tạo ra tải trọng lớn hơn khoảng 6,5% ở các tầng đỉnh, giúp phản ánh rõ hơn ảnh hưởng của các dạng dao động bậc cao.

Làm cách nào để kiểm soát và ngăn chặn khớp dẻo xuất hiện sớm tại chân cột thép?

Kỹ sư cần áp dụng nghiêm ngặt nguyên tắc thiết kế "cột khỏe - dầm yếu" bằng cách lựa chọn các tiết diện cột W14x550 có mô men kháng uốn vượt trội so với dầm W36x160. Đồng thời, việc bố trí sườn gia cường tại mắt nút giúp chuyển vị trí tiêu tán năng lượng dẻo từ thân cột sang các khớp dẻo cục bộ trên dầm, bảo vệ khung không bị sụp đổ đột ngột.

Kết luận

  • Xây dựng thành công mô hình khung thép phẳng 16 tầng 4 nhịp chuẩn xác trên nền tảng OpenSees và kiểm chứng độ tin cậy thông qua SAP2000 v19 với mức sai lệch tần số dao động dưới 2,3%.
  • Định lượng rõ rệt tác động của biến dạng mắt nút: làm tăng chu kỳ dao động riêng thêm 16,4%, tăng độ trôi tầng cực đại từ 18% đến 22% và làm giảm lực cắt đáy cực hạn từ 8,5% đến 11,2%.
  • Làm rõ sự khác biệt trong phân bố tải trọng ngang đẩy dần giữa tiêu chuẩn TCVN 9386:2012 và FEMA 356, chỉ ra mức độ an toàn cao hơn của phương pháp nội suy lũy thừa theo chu kỳ T = 0,67 s.
  • Xác lập cơ chế phân phối lại nội lực, chứng minh vai trò của biến dạng nút trong việc bảo vệ chân cột và phân tán năng lượng địa chấn thông qua các khớp dẻo trên dầm W36.
  • Cung cấp cơ sở dữ liệu học thuật và cẩm nang mô phỏng số hoàn chỉnh cho cộng đồng kỹ sư xây dựng công trình tại Việt Nam.

Định hướng nghiên cứu tiếp theo sẽ mở rộng đánh giá trên mô hình khung không gian 3D chịu tác động của các gia tốc đồ địa chấn theo miền thời gian thực tế trong lộ trình 12 tháng tới. Các kỹ sư và nhà nghiên cứu quan tâm có thể khai thác và áp dụng ngay các thông số từ luận văn để nâng cao chất lượng thiết kế kháng chấn cho các dự án kết cấu thép hiện đại.