Tổng quan nghiên cứu

Thảm họa động đất luôn là một trong những hiểm họa tự nhiên khốc liệt nhất, đe dọa trực tiếp đến tính mạng con người và sự ổn định của cơ sở hạ tầng. Trong thực tế tính toán thiết kế kháng chấn hiện nay, hơn 80% kỹ sư vẫn áp dụng giả thiết chân cột ngàm lý tưởng tại mặt nền (Fixed Base Building - FBB). Giả định này đã bỏ qua hoàn toàn sự tương tác phức tạp giữa móng cọc và môi trường địa chất xung quanh (Soil - Structure Interaction - SSI). Tại các đô thị hiện đại, hơn 70% công trình nhà cao tầng được xây dựng trên hệ móng cọc nằm trong các tầng địa chất có độ biến dạng động phức tạp, khiến mô hình ngàm cứng truyền thống không còn phản ánh chính xác ứng xử chịu lực thực tế.

Mục tiêu trọng tâm của nghiên cứu là xây dựng mô hình toán học và giải thuật toàn diện nhằm phân tích đáp ứng động lực học của kết cấu khung không gian ba chiều (3D) có xét đến tương tác móng cọc - đất nền dưới tác dụng của tải trọng động đất. Luận văn tập trung giải quyết các nhiệm vụ: xác định ma trận độ cứng và hệ số cản động lực học 6 bậc tự do của hệ móng cọc, thiết lập phương trình vi phân chuyển động hoàn chỉnh cho khung không gian 3D, và đánh giá định lượng sự sai lệch về chuyển vị đỉnh, gia tốc, lực cắt đáy cùng mô men uốn chân cột giữa mô hình FBB và SSI. Phạm vi nghiên cứu khảo sát trên mô hình kết cấu thép 9 tầng điển hình chịu sóng gia tốc động đất El Centro, được thực hiện trong giai đoạn từ tháng 8/2017 đến tháng 6/2018 tại Thành phố Hồ Chí Minh. Kết quả nghiên cứu mang lại ý nghĩa then chốt, giúp giảm thiểu rủi ro đánh giá sai lệch chuyển vị lên đến hơn 30%, tối ưu hóa giải pháp kháng chấn và nâng cao hệ số an toàn công trình.

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 được phát triển dựa trên sự tích hợp giữa lý thuyết động lực học kết cấu công trình và cơ học đất động học. Hai lý thuyết nền tảng chi phối bài toán tương tác đất nền - kết cấu bao gồm: lý thuyết tương tác động học (Kinematic Interaction) sinh ra do sự khác biệt về độ cứng giữa kết cấu móng và đất nền xung quanh, và lý thuyết tương tác quán tính (Inertial Interaction) xuất hiện từ chênh lệch khối lượng dao động giữa kết cấu bên trên và nền đất.

Mô hình nghiên cứu kế thừa và hoàn thiện các lý thuyết tính toán độ cứng động học của Milos Novak (1974), George Gazetas (1984) và Ramana (2011). Hệ thống khái niệm chủ đạo bao gồm: hệ số độ cứng động lực học 6 thành phần (3 thành phần tịnh tiến, 2 thành phần xoay quanh trục ngang và 1 thành phần xoắn quanh trục đứng), hệ số giảm chấn động lực học, và ma trận cản Rayleigh tổng thể. Trong mô hình khung không gian 3D, mỗi nút kết cấu sở hữu 6 bậc tự do, ma trận độ cứng phần tử dầm cột có kích thước 12x12, kết hợp cùng ma trận phần tử lò xo móng 12x12 chứa đầy đủ các thành phần độ cứng tịnh tiến k0x, k0y, k0z, độ cứng xoay kx, ky và độ cứng xoắn k(T). Hệ số cản tỉ lệ Rayleigh được xác định thông qua 2 tần số dao động riêng đầu tiên với tỉ số cản vật liệu được kiểm soát chặt chẽ trong khoảng từ 2% đến 5%.

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

Nghiên cứu sử dụng mô hình kết cấu thép 9 tầng không gian làm mô hình đối chứng (benchmark model). Hệ kết cấu bên trên gồm hệ dầm cột thép liên hợp, bên dưới liên kết với hệ móng cọc bê tông cốt thép có chiều dài 17.8m, bán kính quy đổi 0.226m, cọc xuyên qua môi trường địa chất đồng nhất và không đồng nhất. Cỡ mẫu khảo sát bao gồm 4 kịch bản bố trí cọc (nhóm 4 cọc 2x2 và nhóm 9 cọc 3x3 với khoảng cách tim cọc biến thiên từ 3 đến 5 lần bán kính cọc) kết hợp cùng 4 phân loại đất nền từ đất mềm đến đất rất cứng có mô đun đàn hồi biến thiên từ 10 MPa đến 160 MPa. Phương pháp chọn mẫu tập trung vào các dạng kết cấu và địa tầng phổ biến nhất tại các vùng đô thị chịu rủi ro động đất.

Phương pháp phân tích chính là phương pháp tích phân số Newmark trên miền thời gian (Time-History Analysis) với bước tích phân thời gian cực nhỏ dt = 0.02 giây nhằm đảm bảo tính hội tụ và độ chính xác cao khi giải hệ phương trình vi phân chuyển động phi tuyến. Toàn bộ thuật giải được lập trình tự động hóa trên môi trường MATLAB và được đối soát kiểm chứng chu kỳ dao động riêng với phần mềm ETABS, ghi nhận mức sai số tương đối luôn duy trì dưới 2%. Tiến độ nghiên cứu được triển khai liên tục và hoàn thiện trong 10 tháng, bảo đảm độ tin cậy khoa học cao.

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

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

Phân tích số liệu mô phỏng đã chỉ ra sự khác biệt mang tính bản chất giữa mô hình kết cấu ngàm cứng truyền thống (FBB) và mô hình có xét đến tương tác móng - đất nền (SSI) qua các phát hiện nổi bật:

Thứ nhất, chu kỳ dao động riêng của toàn hệ thống khi xét SSI luôn lớn hơn đáng kể so với mô hình FBB. Cụ thể, chu kỳ dao động mode 1 của khung thép 9 tầng tăng từ 1.12 giây ở mô hình FBB lên 1.35 giây ở mô hình SSI (tăng 20.5%), minh chứng rõ nét cho thấy sự xuất hiện của hệ lò xo móng làm tăng tính linh hoạt và độ mềm tổng thể của công trình.

Thứ hai, chuyển vị ngang đỉnh công trình tầng 9 gia tăng rõ rệt khi xét đến SSI. Dưới tác động của trận động đất El Centro, biên độ chuyển vị đỉnh lớn nhất theo phương trục X tăng từ 18% đến 32% tùy thuộc vào độ cứng của nền đất, trong đó môi trường đất mềm ghi nhận mức tăng chuyển vị cực đại lên tới 35.4% so với mô hình ngàm chân cột.

Thứ ba, lực cắt đáy và mô men uốn tại chân cột tầng 1 có xu hướng suy giảm từ 12% đến 28% ở một số vị trí cột giữa do hiệu ứng giảm chấn bức xạ năng lượng vào lòng đất. Tuy nhiên, sự phân bố lại nội lực khiến một số cột biên chịu mô men xoắn gia tăng từ 15% đến 22% khi động đất tác dụng theo phương xiên XY.

Thứ tư, cách bố trí nhóm cọc (khoảng cách cọc S/2Rp = 3 so với S/2Rp = 5) làm thay đổi độ cứng chống xoắn k(T) của đài móng hơn 40%, chứng minh mật độ bố trí cọc chi phối trực tiếp đến khả năng chống xoắn toàn cục của kết cấu không gian.

Thảo luận kết quả

Nguyên nhân vật lý sâu xa của sự gia tăng chu kỳ dao động và chuyển vị đỉnh nằm ở độ lún và độ xoay đàn hồi tại đài móng khi chịu lực kích động chu kỳ. Khi móng bị xoay, góc xoay tích lũy từ chân cột khuếch đại chuyển vị ngang tại các tầng cao. Hiện tượng lực cắt đáy giảm ở mô hình SSI có thể giải thích do chu kỳ dao động dài hơn đã đẩy phản ứng của công trình ra khỏi vùng đỉnh phổ phản ứng gia tốc của sóng địa chấn El Centro.

Kết quả này hoàn toàn tương đồng với các công bố của Novak (1974) và Ramana (2011) trên mô hình lý thuyết đơn giản, đồng thời mở rộng và làm sáng tỏ ứng xử xoắn 3D mà các nghiên cứu khung phẳng của Đặng Văn Út (2016) hay Lương Minh Sang (2014) chưa bao quát được. Về mặt trình bày dữ liệu, các kết quả này được thể hiện trực quan qua hệ thống biểu đồ lịch sử thời gian (Time-History Response) so sánh từng giây chuyển vị đỉnh, biểu đồ bao biến dạng trôi tầng, và các bảng đối chiếu giá trị hiệu dụng căn quân phương (Root Mean Square - RMS) với độ tin cậy toán học cao, cho thấy rõ mức độ phân kỳ giữa hai phương pháp tính toán.

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

Nhằm nâng cao chất lượng thiết kế kháng chấn và thu hẹp khoảng cách giữa mô hình tính toán lý thuyết với ứng xử công trình thực tế, các khuyến nghị cụ thể được đề xuất như sau:

Một là, tích hợp bắt buộc ma trận lò xo độ cứng 6 bậc tự do (gồm độ cứng tịnh tiến, xoay và xoắn) vào các phần mềm thiết kế thương mại như ETABS hoặc SAP2000 khi tính toán công trình từ 8 tầng trở lên trên nền đất yếu. Mục tiêu giảm sai số dự báo chuyển vị ngang xuống dưới 5%, lộ trình triển khai áp dụng tại các đơn vị tư vấn thiết kế trong vòng 12 đến 24 tháng tới do các kỹ sư kết cấu chủ trì.

Hai là, bổ sung và chuẩn hóa các chỉ dẫn tính toán tương tác móng cọc - đất nền vào Quy chuẩn và Tiêu chuẩn Thiết kế Công trình Chịu Động Đất quốc gia. Mục tiêu hoàn thiện hệ số điều chỉnh phản ứng động học cho từng phân loại đất nền (đặc biệt là đất loại D và E), thực hiện trong giai đoạn 2026 - 2028 dưới sự phối hợp của Viện Khoa học Công nghệ Xây dựng và các trường đại học chuyên ngành.

Ba là, tối ưu hóa khoảng cách bố trí cọc trong đài móng với tỷ lệ khoảng cách tim cọc đạt từ 3 đến 5 lần đường kính cọc. Giải pháp này giúp tăng hiệu quả độ cứng xoắn của nhóm cọc lên trên 35%, đồng thời tiết kiệm 10% đến 15% khối lượng bê tông đài móng, được thực hiện bởi các kỹ sư địa kỹ thuật ngay từ giai đoạn thiết kế cơ sở.

Bốn là, đẩy mạnh công tác quan trắc dao động động lực học thực tế (Full-scale Dynamic Testing) cho 100% các công trình cao tầng cấp đặc biệt nhằm hiệu chỉnh mô hình số. Thời gian đo đạc định kỳ 6 tháng một lần sau khi hoàn thành xây dựng do các trung tâm kiểm định chất lượng công trình thực hiện.

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

Nội dung và giải thuật trong luận văn mang lại giá trị thực tiễn và học thuật to lớn cho 4 nhóm đối tượng chính:

Thứ nhất, Kỹ sư thiết kế kết cấu công trình dân dụng và công nghiệp. Luận văn cung cấp phương pháp luận chi tiết và bộ công thức toán học tường minh để xác định độ cứng lò xo móng 3D, giúp kỹ sư dễ dàng đưa vào mô hình tính toán để kiểm tra độ trôi tầng và chuyển vị đỉnh chính xác, giảm thiểu nguy cơ nứt vỡ kết cấu bao che.

Thứ hai, Học viên cao học, nghiên cứu sinh và giảng viên ngành Kỹ thuật Xây dựng. Đề tài là tài liệu tham khảo mẫu mực về quy trình giải bài toán động lực học phi tuyến trên miền thời gian, cung cấp toàn bộ mã nguồn chương trình MATLAB để phát triển các hướng nghiên cứu chuyên sâu về tương tác kết cấu - địa kỹ thuật.

Thứ ba, Đơn vị thẩm tra và tư vấn giám sát xây dựng. Dữ liệu từ luận văn là cơ sở tin cậy giúp các chuyên gia thẩm tra phát hiện các rủi ro tiềm ẩn khi nhà thầu chỉ sử dụng mô hình ngàm cứng FBB cho các dự án xây dựng trên nền đất yếu có nguy cơ khuếch đại chuyển vị ngang.

Thứ tư, Cơ quan quản lý nhà nước và các viện nghiên cứu chuyên ngành xây dựng. Báo cáo nghiên cứu đóng góp cơ sở dữ liệu thực nghiệm số quan trọng, hỗ trợ đắc lực cho việc biên soạn, sửa đổi và hoàn thiện các tiêu chuẩn thiết kế kháng chấn tại Việt Nam.

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

Mô hình ngàm chân cột truyền thống (FBB) gây ra những sai lệch nghiêm trọng nào khi tính toán động đất cho nhà cao tầng? Mô hình FBB giả định nền đất cứng tuyệt đối và triệt tiêu toàn bộ chuyển vị tại đáy cột. Trong thực tế, đất nền luôn bị biến dạng động, làm chu kỳ dao động riêng của công trình tăng từ 15% đến 25%. Bỏ qua SSI khiến kỹ sư đánh giá thấp chuyển vị đỉnh thực tế tới hơn 30%, dẫn đến nguy cơ va đập giữa các công trình liền kề.

Độ cứng xoắn của móng cọc k(T) có vai trò như thế nào trong mô hình khung không gian 3D? Độ cứng xoắn k(T) kiểm soát dao động xoắn ngoài mặt phẳng của công trình khi chịu tải trọng động đất tác dụng theo phương bất kỳ hoặc phương xiên XY. Nếu bỏ qua thành phần này, độ cứng chống xoắn tổng thể bị sai lệch hơn 40%, làm tăng nội lực xoắn cục bộ tại các cột góc từ 10% đến 20%.

Đặc tính của nền đất ảnh hưởng như thế nào đến mức độ tương tác giữa móng và kết cấu? Độ cứng của nền đất tỉ lệ thuận với độ cứng động lực học của hệ lò xo móng. Nền đất càng mềm (mô đun cắt thấp), ảnh hưởng của SSI càng rõ rệt, khiến chu kỳ dao động tăng trên 30% và chuyển vị ngang tăng hơn 25% so với nền đất cứng, đòi hỏi phải tăng cường tiết diện hoặc gia cố móng.

Tại sao phương pháp tích phân Newmark trên MATLAB lại vượt trội trong việc phân tích SSI? Thuật toán Newmark giải trực tiếp phương trình vi phân chuyển động từng bước trên miền thời gian với bước chia dt = 0.02 giây, cho phép cập nhật đồng thời ma trận độ cứng động lực học và ma trận cản Rayleigh của cả kết cấu lẫn đất nền mà không làm mất đi các mode dao động bậc cao như các phương pháp giải gần đúng.

Khi nào bắt buộc phải áp dụng mô hình phân tích có xét đến tương tác móng - đất nền (SSI)? Phân tích SSI là bắt buộc đối với các công trình có chiều cao trên 8 tầng, công trình có tỷ lệ chiều cao trên bề rộng lớn hơn 3, hoặc xây dựng trên nền đất yếu thuộc phân loại D, E với đỉnh gia tốc nền địa chấn lớn hơn 0.1g theo quy chuẩn an toàn kháng chấn.

Kết luận

Tổng kết những đóng góp cốt lõi và định hướng tương lai của đề tài:

  • Xây dựng thành công thuật toán và mô hình toán học giải bài toán tương tác móng cọc - đất nền (SSI) cho khung không gian 3D đầy đủ 6 thành phần độ cứng động lực học.
  • Khẳng định mô hình SSI làm tăng chu kỳ dao động riêng hơn 20% và tăng chuyển vị đỉnh từ 18% đến 32% so với mô hình ngàm FBB dưới tác động của động đất El Centro.
  • Phát hiện sự suy giảm lực cắt chân cột tầng 1 từ 12% đến 28% kết hợp với sự phân bố lại mô men uốn tại các cột góc khi xét động đất tác dụng theo phương xiên XY.
  • Chứng minh vai trò then chốt của khoảng cách bố trí nhóm cọc (S/2Rp = 3 đến 5), quyết định trực tiếp tới hơn 40% độ cứng chống xoắn của móng đài.
  • Đề xuất quy trình ứng dụng thực tiễn và chuyển giao thuật toán Newmark trên MATLAB phục vụ công tác thiết kế kháng chấn chính xác trong công nghiệp xây dựng.

Đóng góp khoa học lớn nhất của luận văn là đã hoàn thiện bài toán phân tích động lực học kết cấu 3D có xét đến đầy đủ độ cứng xoay và xoắn của móng cọc trong môi trường đất nền nhiều lớp. Trong lộ trình 12 tháng tới, hướng nghiên cứu tiếp theo sẽ mở rộng khảo sát ứng xử phi tuyến hình học và phi tuyến vật liệu của đất nền, kết hợp hiệu ứng cản nhớt phi tuyến nhằm nâng cao hơn nữa độ chuẩn xác. Các kỹ sư kết cấu, nhà nghiên cứu và đơn vị tư vấn hãy áp dụng ngay phương pháp luận tính toán SSI vào các dự án thực tế để bảo đảm sự an toàn bền vững cho mọi công trình trước thảm họa thiên tai.