Tổng quan nghiên cứu

Đánh giá an toàn kết cấu công trình chịu động đất là một bài toán kỹ thuật trọng yếu trong bối cảnh các đô thị hiện đại ngày càng phát triển nhiều công trình cao tầng. Trong thực tế tính toán kết cấu kháng chấn, phương pháp phân tích phản ứng phi tuyến theo miền thời gian (NL-RHA) được xem là thước đo chuẩn xác nhất về ứng xử thực tế của công trình khi chịu các rung chấn mạnh với xác suất xảy ra 10% trong 50 năm (chu kỳ lặp khoảng 475 năm). Tuy nhiên, phương pháp NL-RHA đòi hỏi khối lượng tài nguyên tính toán rất lớn, mô hình hóa phức tạp và thời gian xử lý kéo dài, gây khó khăn lớn khi áp dụng trực tiếp vào các dự án thiết kế thông thường.

Để giải quyết thách thức này, các kỹ sư thường dùng phương pháp phân tích tĩnh phi tuyến đẩy dần tiêu chuẩn (SPA) theo các tiêu chuẩn quốc tế như FEMA-356 và ATC-40. Mặc dù vậy, phương pháp SPA chỉ phù hợp với các công trình thấp tầng, và bộc lộ sai số lớn lên đến hơn 35% đối với các công trình nhiều tầng do bỏ qua sự đóng góp của các dạng dao động bậc cao.

Nghiên cứu của học viên Nguyễn Thị Thảo Hiền tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh tập trung làm rõ và ứng dụng phương pháp phân tích tĩnh đẩy dần có xét đến các dạng dao động cao (MPA) kết hợp với mô hình phần tử dầm - cột có khớp thớ (Fiber Model). Nghiên cứu được triển khai trên 3 mô hình khung bê tông cốt thép đại diện gồm 6 tầng, 12 tầng và 20 tầng trong khoảng thời gian từ tháng 8 năm 2013 đến tháng 6 năm 2014. Mục tiêu cốt lõi là đánh giá độ chính xác, độ trôi tầng và mức độ hư hại của kết cấu từ phương pháp MPA và SPA khi đối chiếu với kết quả chuẩn xác từ NL-RHA, cung cấp giải pháp tính toán tối ưu giúp giảm thiểu thời gian mô phỏng mà vẫn đảm bảo độ tin cậy trên 85% cho các công trình cao tầng.

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 xây dựng trên nền tảng động lực học công trình phi tuyến và lý thuyết biến dạng dẻo của vật liệu bê tông cốt thép. Hai khung lý thuyết chính được áp dụng xuyên suốt gồm:

Thứ nhất là lý thuyết phân tích đẩy dần theo từng dạng dao động (Modal Pushover Analysis - MPA) do Chopra và Goel đề xuất năm 2001. Khung lý thuyết này phân tích kết cấu nhiều bậc tự do (MDOF) không đàn hồi thông qua việc chuyển đổi phản ứng của từng mode dao động sang hệ một bậc tự do (SDOF) không đàn hồi tương đương, sau đó tích phân phản ứng đỉnh theo quy tắc căn bậc hai của tổng bình phương (SRSS).

Thứ hai là mô hình phi tuyến vật liệu với mô hình ứng suất - biến dạng của Mander (1984) cho bê tông chịu nén có xét đến hiệu ứng ép ngang của cốt đai (confined concrete). Mô hình này xác định cường độ nén dọc trục và hệ số nén ngang hiệu quả ke, giúp phản ánh chính xác khả năng tiêu tán năng lượng và gia tăng độ dẻo của lõi bê tông. Cốt thép được mô phỏng theo mô hình song tuyến tính có xét đến hiện tượng cứng hóa biến dạng với hệ số giảm độ cứng sau chảy dẻo bằng 0.03 (tương đương 3%).

Các khái niệm chính trong đề tài bao gồm: mô hình phần tử dầm - cột có khớp thớ (Fiber Model) chia tiết diện thành nhiều thớ vi mô để theo dõi quá trình lan truyền dẻo; độ trôi tầng cực đại (Maximum Inter-story Drift Ratio); chuyển vị mục tiêu tại đỉnh mái (Target Roof Displacement); và chu kỳ dao động tự nhiên của hệ kết cấu.

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

Mẫu nghiên cứu gồm 3 mô hình khung bê tông cốt thép phẳng đối xứng với chiều cao lần lượt là 6 tầng (công trình thấp tầng), 12 tầng (công trình trung tầng) và 20 tầng (công trình cao tầng). Mỗi nhịp khung có chiều dài 6 mét và chiều cao mỗi tầng là 3.6 mét, được thiết kế theo nguyên tắc cột mạnh - dầm yếu. Dữ liệu động đất đầu vào sử dụng tổ hợp 5 gia tốc đồ dao động nền của các trận động đất mạnh tại Los Angeles thuộc dự án SAC, có xác suất vượt quá 10% trong 50 năm.

Phương pháp phân tích được thực hiện thông qua quy trình 9 bước kết hợp hai công cụ mô phỏng chuyên sâu: phần mềm Drain-2DX được sử dụng để xây dựng mô hình khớp thớ (Type 15) và thực hiện phân tích đẩy dần phi tuyến; phần mềm MATLAB R2011b đảm nhiệm việc tính toán phản ứng của hệ SDOF không đàn hồi và giải thuật tích phân thời gian. Toàn bộ nghiên cứu được thực hiện trong thời gian 10 tháng. Phương pháp phân tích bằng mô hình khớp thớ được lựa chọn vì tính chính xác vượt trội trong việc nắm bắt tương tác lực dọc - mô men uốn (P-M-M) và sự suy giảm độ cứng cục bộ mà không đòi hỏi chi phí giải tích quá cao như mô hình phần tử hữu hạn 3D dạng khối.

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 số liệu trên 3 hệ khung bê tông cốt thép đã đưa ra 4 phát hiện quan trọng:

Đối với khung 6 tầng, dạng dao động đầu tiên (Mode 1) chiếm hơn 82% tỷ lệ khối lượng tham gia dao động. Kết quả phân tích cho thấy cả hai phương pháp SPA và MPA đều đưa ra kết quả chuyển vị đỉnh và độ trôi tầng rất gần với phương pháp chính xác NL-RHA, với sai số trung bình chỉ dao động từ 5.2% đến 8.5%.

Đối với khung 12 tầng và 20 tầng, đóng góp của các dạng dao động bậc cao (Mode 2 và Mode 3) tăng mạnh, chiếm từ 25% đến 35% tổng phản ứng lực cắt và chuyển vị. Phương pháp đẩy dần tiêu chuẩn SPA đánh giá thấp độ trôi tầng cực đại ở các tầng phía trên từ 28% đến 42%, tạo ra nguy cơ mất an toàn nghiêm trọng trong thiết kế kháng chấn.

Phương pháp MPA khi kết hợp 3 dạng dao động đầu tiên bằng quy tắc SRSS đã cải thiện rõ rệt độ chính xác, kéo giảm độ sai lệch so với NL-RHA xuống dưới mức 12% đối với khung 12 tầng và dưới 16.5% đối với khung 20 tầng.

Về biến dạng dẻo của vật liệu, mô hình khớp thớ chỉ ra rằng các khớp dẻo xuất hiện đầu tiên tại các đầu dầm ở khoảng 60% chiều cao công trình trước khi phát triển dần xuống chân cột tầng trệt. Độ dẻo uốn của các cấu kiện có bố trí thép đai theo mô hình Mander đạt giá trị lớn hơn 3.5 lần so với cấu kiện không có đai gia cường.

Thảo luận kết quả

Sự khác biệt rõ rệt giữa phương pháp SPA và MPA ở các khung 12 và 20 tầng bắt nguồn từ việc phân bố lực tĩnh của SPA không thể phản ánh sự thay đổi quán tính liên tục khi các tầng trên bị kích động bởi các mode dao động bậc cao. Khi các kết quả được biểu diễn trực quan trên đồ thị phân bố độ trôi tầng theo chiều cao và đường cong quan hệ lực cắt đáy - chuyển vị đỉnh (Pushover Curve), các đỉnh biến dạng của khung 20 tầng luôn có xu hướng phình to ở khu vực từ tầng 12 đến tầng 16. Phương pháp MPA mô tả chính xác hình dạng đường cong uốn lượn này, trong khi đường cong của SPA có dạng tuyến tính đơn giản làm mất đi các vị trí nguy hiểm cục bộ.

Khi so sánh với các nghiên cứu kinh điển trên thế giới của Chopra và Chintanapakdee trên kết cấu khung thép, nghiên cứu này chứng minh rằng quy luật chi phối của các dạng dao động bậc cao vẫn giữ nguyên bản chất đối với kết cấu bê tông cốt thép. Việc sử dụng mô hình khớp thớ đã nâng cao độ nhạy cảm của mô hình với hiện tượng nứt bê tông và chảy dẻo cốt thép, giúp thu hẹp khoảng cách giữa mô phỏng số và ứng xử thực nghiệm mà các mô hình khớp dẻo tập trung truyền thống không làm được.

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

Dựa trên các kết quả định lượng đạt được, đề tài đưa ra 4 nhóm khuyến nghị mang tính ứng dụng cao cho ngành xây dựng công trình:

Thứ nhất, chuẩn hóa việc áp dụng phương pháp phân tích MPA vào các tiêu chuẩn thiết kế kháng chấn quốc gia cho các công trình bê tông cốt thép có chiều cao từ 10 tầng trở lên (hoặc chu kỳ dao động cơ bản lớn hơn 1.0 giây). Cơ quan quản lý chuyên ngành và các viện nghiên cứu xây dựng cần sớm ban hành tài liệu hướng dẫn kỹ thuật trong thời gian 12 đến 24 tháng tới nhằm khống chế sai số thiết kế dưới 15%.

Thứ hai, tích hợp mô hình khớp thớ vào quy trình thiết kế của các công ty tư vấn kỹ thuật kết cấu thay cho mô hình khớp dẻo đơn giản. Việc này giúp kiểm soát chính xác mức độ dẻo yêu cầu lớn hơn 3.0 tại các nút khung dầm - cột, đảm bảo mục tiêu bảo vệ công trình không sụp đổ dưới các trận động đất cấp 7 và cấp 8.

Thứ ba, thiết lập quy chuẩn kiểm tra bắt buộc về độ trôi tầng giới hạn, khống chế độ trôi tầng không vượt quá 1.5% đối với trạng thái sử dụng bình thường và 2.0% đối với trạng thái giới hạn phá hoại cho 100% công trình nhà cao tầng tại các khu vực có gia tốc nền thiết kế lớn hơn 0.08g.

Thứ tư, mở rộng nghiên cứu sang không gian 3 chiều có xét đến hiệu ứng xoắn và tính bất đối xứng mặt bằng. Các nhóm nghiên cứu tại các trường đại học kỹ thuật nên triển khai thử nghiệm với tập hợp từ 10 đến 20 trận động đất thực tế để hoàn thiện phần mềm tự động hóa phân tích phi tuyến trong lộ trình 3 năm tới.

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

Nội dung và dữ liệu thực nghiệm của luận văn mang lại giá trị chuyên môn thiết thực cho 4 nhóm đối tượng:

Kỹ sư thiết kế kết cấu công trình cao tầng: Có thể áp dụng ngay quy trình 9 bước của phương pháp MPA để tính toán kiểm tra động đất cho các tòa nhà từ 10 đến 25 tầng, giúp tiết kiệm hơn 70% thời gian chạy mô phỏng so với phương pháp NL-RHA mà vẫn đảm bảo độ tin cậy kết cấu.

Chuyên gia thẩm tra và ban quản lý dự án xây dựng: Sử dụng các tiêu chí về độ trôi tầng và mô hình khớp thớ trong luận văn làm cơ sở kỹ thuật để đánh giá mức độ an toàn kháng chấn của các dự án chung cư, văn phòng quy mô lớn.

Học viên cao học và nghiên cứu sinh ngành Xây dựng Dân dụng: Tham khảo phương pháp thiết lập ma trận độ cứng, kỹ thuật mô hình hóa vật liệu phi tuyến Mander và quy trình kết nối dữ liệu giữa Drain-2DX với MATLAB.

Các đơn vị biên soạn tiêu chuẩn và phần mềm xây dựng: Khai thác hệ thống số liệu so sánh định lượng giữa SPA, MPA và NL-RHA để làm dữ liệu kiểm chuẩn (benchmark) khi phát triển các module tính toán kháng chấn mới.

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

Tại sao phương pháp phân tích đẩy dần chuẩn SPA lại cho sai số lớn đối với công trình trên 10 tầng? Phương pháp SPA chỉ áp dụng mô hình phân bố lực quán tính cố định dựa trên dạng dao động cơ bản đầu tiên. Khi công trình cao trên 10 tầng (khoảng trên 30 mét), chu kỳ dao động kéo dài khiến các mode dao động bậc 2 và bậc 3 tham gia mạnh vào phản ứng công trình, gây ra lực cắt tầng và độ trôi tầng lớn ở nửa trên tòa nhà mà SPA không thể dự báo được.

Mô hình khớp thớ Fiber Model có điểm gì vượt trội so với mô hình khớp dẻo tập trung? Mô hình khớp dẻo tập trung giả định toàn bộ biến dạng dẻo chỉ gom lại tại một điểm với các quan hệ lực - biến dạng gán sẵn. Ngược lại, mô hình khớp thớ chia mặt cắt thành nhiều thớ bê tông và cốt thép riêng biệt, tự động cập nhật quan hệ ứng suất - biến dạng phi tuyến và tương tác lực dọc - uốn theo từng bước tải, mang lại độ chính xác cao hơn 25% trong phân tích phá hoại.

Quy tắc tổ hợp căn bậc hai của tổng bình phương SRSS trong phương pháp MPA hoạt động như thế nào? Quy tắc SRSS tính toán phản ứng tổng thể của kết cấu bằng cách lấy căn bậc hai của tổng bình phương các giá trị cực đại thu được từ từng mode dao động riêng biệt. Trong nghiên cứu này, việc kết hợp phản ứng của 3 dạng dao động đầu tiên bằng SRSS giúp triệt tiêu sai số dao động cùng pha và cho kết quả sát với phân tích miền thời gian thực tế.

Dữ liệu động đất nào được sử dụng để kiểm chứng độ tin cậy trong nghiên cứu? Nghiên cứu sử dụng 5 gia tốc đồ dao động nền từ các trận động đất thực tế đã được chuẩn hóa theo chương trình SAC tại khu vực Los Angeles. Các trận động đất này có xác suất xuất hiện 10% trong vòng 50 năm, tương ứng với chu kỳ lặp 475 năm, đại diện cho mức động đất thiết kế tiêu chuẩn của các công trình công cộng.

Các kỹ sư có thể ứng dụng kết quả nghiên cứu vào phần mềm thương mại hiện nay như thế nào? Các kỹ sư hoàn toàn có thể áp dụng nguyên lý 9 bước của phương pháp MPA vào các phần mềm thương mại phổ biến như ETABS hoặc SAP2000 bằng cách khai báo các mô hình tải trọng tỉ lệ với từng dạng dao động riêng biệt, thực hiện phân tích pushover cho từng mode và kết hợp kết quả chuyển vị bằng bảng tính Excel.

Kết luận

  • Luận văn đã chứng minh tính hiệu quả vượt trội của phương pháp phân tích tĩnh phi tuyến MPA kết hợp mô hình khớp thớ trong việc đánh giá an toàn kháng chấn cho khung bê tông cốt thép.
  • Phương pháp MPA giúp kiểm soát sai số dưới 15% so với phương pháp chuẩn NL-RHA trên các hệ khung 6, 12 và 20 tầng, trong khi phương pháp đẩy dần tiêu chuẩn SPA có thể sai lệch tới hơn 40%.
  • Xác định rõ vai trò quyết định của các dạng dao động bậc cao đối với phản ứng độ trôi tầng ở các công trình có chiều cao trên 12 tầng chịu tác động của động đất.
  • Mô hình vật liệu Mander cho bê tông có đai gia cường đã mô phỏng thành công khả năng tăng độ dẻo lên hơn 3.5 lần, giúp giảm thiểu nguy cơ phá hoại giòn của cấu kiện chịu lực.
  • Đề xuất quy trình tính toán 9 bước hoàn chỉnh, tạo tiền đề khoa học vững chắc để phát triển các công cụ phân tích kháng chấn 3D tự động hóa trong giai đoạn 12 đến 36 tháng tới.

Các đơn vị tư vấn và kỹ sư kết cấu nên chủ động áp dụng quy trình phân tích MPA và mô hình khớp thớ vào thực tiễn thiết kế công trình cao tầng nhằm nâng cao độ an toàn và tối ưu hóa chi phí kháng chấn cho công trình.