Tổng quan nghiên cứu

Trong kỹ thuật xây dựng công trình hiện đại, nguy cơ động đất đang trở thành thách thức lớn đối với độ an toàn của các công trình nhiều tầng tại các đô thị phát triển. Các ghi nhận thực tế từ những trận động đất kinh điển như El Centro với gia tốc nền cực đại đạt khoảng 0,35g hay Kobe với gia tốc vượt mức 0,82g đã chứng minh sức tàn phá nghiêm trọng của sóng địa chấn đối với hệ thống kết cấu hạ tầng. Tuy nhiên, trong quy trình thiết kế truyền thống, kỹ sư thường áp dụng mô hình chân cột ngàm cứng bỏ qua sự mềm của đất nền và cọc (Fixed Base Building - FBB), dẫn đến việc đánh giá sai lệch chu kỳ dao động tự nhiên và phản ứng động lực học của công trình.

Luận văn thạc sĩ chuyên ngành Kỹ thuật 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 - Đại học Quốc gia TP.HCM của tác giả Nguyễn Văn Hiến, dưới sự hướng dẫn của PGS.TS. Chu Quốc Thắng và ThS.NCS. Phạm Nhân Hòa, tập trung giải quyết bài toán cốt lõi: Phân tích sự thay đổi độ cứng động lực học của móng cọc khi chịu tải trọng động đất và đánh giá ảnh hưởng tương tác đất nền - kết cấu (Soil - Structure Interaction - SSI).

Mục tiêu cụ thể của đề tài là xây dựng thuật toán xác định sự suy giảm mô đun cắt của đất theo thời gian thực, tính toán ma trận độ cứng và cản động lực học tương đương cho móng cọc, từ đó phân tích ứng xử động của khung thép 9 tầng và 20 tầng chịu tác động của sóng địa chấn. Nghiên cứu được triển khai thực hiện trong mốc thời gian từ tháng 12 năm 2017 đến tháng 06 năm 2018. Kết quả nghiên cứu mang ý nghĩa thực tiễn to lớn khi giúp giảm thiểu sai số ước tính nội lực và chuyển vị từ 15% đến 35%, cung cấp cơ sở khoa học tin cậy để tối ưu hóa thiết kế kháng chấn cho các công trình cao tầng đặt trên nền móng cọc.

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 xây dựng trên nền tảng lý thuyết tương tác đất nền - kết cấu (SSI) hiện đại, phân tách toàn diện thành hai hiện tượng vật lý chính: tương tác động học (Kinematic Interaction) phát sinh do sự khác biệt về độ cứng giữa kết cấu móng và môi trường đất xung quanh; và tương tác quán tính (Inertial Interaction) sinh ra do quán tính khối lượng của công trình truyền ngược xuống nền đất khi dao động.

Khung lý thuyết của luận văn tích hợp các mô hình tính toán động lực học móng cọc kinh điển của Novak (1974), Gazetas và Dobry (1984), Poulos (1968) cùng mô hình hệ số tương tác của Beredugo và Novak (1972). Theo đó, móng cọc phức tạp được quy đổi tương đương về hệ các phần tử lò xo và giảm chấn 3 bậc tự do gồm: độ cứng ngang $k_x(t)$, độ cứng đứng $k_z(t)$, và độ cứng xoay $k_\theta(t)$ cùng các hệ số cản nhớt tương ứng $c_x(t), c_z(t), c_\theta(t)$. Độ cứng động lực học của cọc không phải là hằng số mà biến thiên phi tuyến theo thời gian dựa trên quan hệ ứng suất cắt $\tau(t)$ và biến dạng trượt $\gamma(t)$, dẫn đến sự suy giảm liên tục của mô đun cắt của đất $G_s(t)$ và biến thiên của hệ số Poisson $\nu(t)$ trong suốt quá trình rung chấn.

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

Nguồn dữ liệu tính toán được thiết lập từ các thông số địa kỹ thuật đất phân lớp thực tế với mô đun đàn hồi $E_s$ biến thiên từ 12 MPa ở lớp đất chứa đài đến 70 MPa ở lớp chịu mũi, khối lượng riêng của đất đạt 1937 kg/m³, và vận tốc truyền sóng cắt $V_s$ từ 79 m/s đến 120 m/s. Cơ sở dữ liệu động đất sử dụng các băng gia tốc nền tiêu chuẩn quốc tế El Centro và Kobe với bước thời gian phân tích chi tiết $\Delta t = 0,02$ giây.

Về cỡ mẫu và phương pháp chọn mẫu, nghiên cứu lựa chọn phương pháp mẫu chuẩn mực điển hình (benchmark case study) gồm 2 mô hình kết cấu khung thép 5 nhịp đại diện cho hai phân khúc công trình quan trọng: công trình trung tầng (9 tầng, tổng khối lượng khoảng 9003 kN) và công trình cao tầng (20 tầng). Cấu hình móng cọc được chọn mẫu theo tiêu chuẩn thực tế với đài cọc $3 \times 3$ gồm 9 cọc bê tông cốt thép M350 có đường kính 0,8 m, mô đun đàn hồi cọc $E_p = 30$ GPa, chiều dài cọc $L_p$ từ 19 m đến 35 m và khoảng cách tim cọc $S = 2,0$ m (tỉ số $S/2R_p = 5$).

Lý do lựa chọn phương pháp phân tích số là việc sử dụng phương pháp phần tử hữu hạn (FEM) kết hợp thuật toán Newmark gia tốc tuyến tính ($\gamma = 0,5; \beta = 1/6$) giải trực tiếp hệ phương trình vi phân chuyển động trên miền thời gian. Phương pháp này đảm bảo tính ổn định số vô điều kiện, phản ánh chính xác tính chất phi tuyến tức thời của hệ lò xo chân cột mà không làm quá tải khối lượng tính toán so với các mô hình 3D toàn phần phức tạp. Toàn bộ thuật toán được tác giả trực tiếp lập trình và thực thi thông qua mã nguồn MATLAB chuyên biệt.

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

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

Nghiên cứu mang lại những kết quả định lượng rõ ràng thông qua việc so sánh đối chứng giữa 3 mô hình: ngàm chân cột truyền thống (FBB), tương tác đất nền tĩnh ($SSI_{sta}$), và tương tác đất nền động biến thiên ($SSI_{dyn}$):

  1. Sự suy giảm mạnh của độ cứng động lực học móng: Dưới tác động rung lắc mạnh của gia tốc nền, mô đun cắt hiệu dụng $G_s(t)$ của nền đất suy giảm mạnh từ 40% đến 65% so với giá trị ban đầu tại các thời điểm gia tốc đạt đỉnh. Kéo theo đó, độ cứng ngang $k_x(t)$ và độ cứng xoay $k_\theta(t)$ của hệ móng cọc giảm từ 30% đến 55%, thể hiện sự mềm hóa rõ rệt của hệ thống liên kết móng - đất.
  2. Hiện tượng kéo dài chu kỳ dao động tự nhiên: Khi xét đến liên kết mềm và biến thiên độ cứng trong mô hình $SSI_{dyn}$, chu kỳ dao động cơ bản của công trình 9 tầng tăng từ 18% đến 24%, và công trình 20 tầng tăng khoảng 12% đến 17% so với mô hình chân ngàm FBB.
  3. Gia tăng chuyển vị đỉnh công trình: Do đài móng bị xoay và trượt tịnh tiến dưới tác động động đất, chuyển vị ngang cực đại tại đỉnh công trình của mô hình $SSI_{dyn}$ tăng từ 15% đến 28% so với mô hình FBB. Điều này cho thấy mô hình truyền thống đã đánh giá thấp mức độ dao động hình học của công trình.
  4. Giảm lực cắt đáy và mô men uốn chân cột: Ngược lại với chuyển vị, phản lực tiếp xúc tại chân kết cấu bao gồm lực cắt đáy và mô men uốn chân cột trong mô hình $SSI_{dyn}$ giảm từ 10% đến 22% so với mô hình FBB. Nguyên nhân chủ yếu do năng lượng sóng dao động đã bị tiêu tán đáng kể thông qua hệ số cản bức xạ sóng và cản nhớt của môi trường đất nền ($c_x(t), c_\theta(t)$).

Thảo luận kết quả

Cơ chế vật lý đằng sau những phát hiện trên nằm ở sự biến dạng trượt phi tuyến của các lớp đất bao quanh thân cọc và mũi cọc. Khi sóng địa chấn lan truyền từ đá gốc lên bề mặt, đất nền bị biến dạng trượt làm giảm độ cứng chống cắt tức thời, khiến công trình dao động với biên độ góc xoay lớn hơn tại cốt không (mặt đài).

Kết quả này hoàn toàn nhất quán với các nghiên cứu lý thuyết của Novak (1974) và Gazetas (1984), đồng thời làm sáng tỏ các điểm nghẽn mà các nghiên cứu trước đây tại Việt Nam chưa giải quyết trọn vẹn (vốn thường xem độ cứng móng cọc là hằng số hoặc chỉ khảo sát nền đất đồng nhất một lớp).

Để phục vụ công tác thiết kế trực quan, toàn bộ tập dữ liệu phản ứng động lực học có thể được biểu diễn thông qua biểu đồ lịch sử thời gian (time-history response) mô tả sự biến thiên của $k(t)$ và $c(t)$ song song với biểu đồ chuyển vị đỉnh theo thời gian. Bảng tổng hợp so sánh độ lệch phần trăm giữa FBB và SSI cung cấp cho kỹ sư cái nhìn định lượng chính xác để hiệu chỉnh hệ số an toàn trong các phần mềm tính toán kết cấu phổ biến.

Thông số phản ứng động Mô hình ngàm FBB Mô hình $SSI_{sta}$ Mô hình $SSI_{dyn}$ Mức độ thay đổi ($SSI_{dyn}$ so với FBB)
Độ cứng móng $k_x, k_\theta$ Vô cùng (Ngàm cứng) Hằng số tĩnh Biến thiên theo $G_s(t)$ Giảm từ 30% đến 55%
Chu kỳ dao động $T_1$ (s) Chu kỳ chuẩn Tăng 8% - 12% Tăng 12% - 24% Tăng đáng kể (công trình mềm hơn)
Chuyển vị đỉnh cực đại (m) Thấp nhất Tăng 10% - 15% Tăng 15% - 28% Tăng rủi ro va chạm khe co giãn
Lực cắt & Mô men đáy (kN, kNm) Lớn nhất (Thiên an toàn) Giảm 5% - 10% Giảm 10% - 22% Giảm tải trọng thiết kế móng

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

Dựa trên các kết luận rút ra từ mô hình tính toán số, luận văn đưa ra 4 nhóm khuyến nghị kỹ thuật có tính ứng dụng cao:

  1. Bắt buộc tích hợp mô hình SSI trong tiêu chuẩn thiết kế kháng chấn: Cơ quan quản lý ngành xây dựng và các hội đồng thẩm định cần yêu cầu bắt buộc áp dụng mô hình phân tích tương tác đất nền - kết cấu đối với các công trình có quy mô từ 9 tầng trở lên hoặc xây dựng trên nền đất yếu (loại D, E theo TCVN 9386:2012). Mục tiêu là kiểm soát độ võng và chuyển vị ngang tương đối giữa các tầng trong giới hạn an toàn từ $h/250$ đến $h/500$.
  2. Cập nhật quy trình xác định độ cứng lò xo móng cọc: Các công ty tư vấn thiết kế kết cấu cần chuyển đổi từ việc gán ngàm cứng tuyệt đối sang việc khai báo hệ lò xo đàn nhớt tương đương 3 phương ($k_x, k_z, k_\theta$) có xét đến sự suy giảm độ cứng theo thời gian. Lộ trình triển khai áp dụng khuyến nghị này có thể thực hiện trong vòng 6 đến 12 tháng thông qua các module liên kết giữa MATLAB và phần mềm ETABS/SAP2000.
  3. Tối ưu hóa bố trí hình học nhóm cọc: Khi thiết kế móng cọc chịu tải trọng động đất, cần duy trì tỉ số khoảng cách giữa các tim cọc $S/2R_p \ge 3$ đến $5$ nhằm hạn chế tối đa hiệu ứng cộng hưởng nhóm cọc và giảm thiểu sự suy giảm đột ngột của hệ số độ cứng nhóm cọc $\alpha_L$.
  4. Ban hành hướng dẫn tính toán thực hành chuyên sâu: Các viện nghiên cứu chuyên ngành cần phối hợp xây dựng sổ tay kỹ thuật chuẩn hóa quy trình phân tích SSI thu gọn, giúp các đơn vị thiết kế tiết kiệm khoảng 20% chi phí vật liệu kết cấu ngầm do không phải thiết kế dôi dư theo các giả định an toàn thái quá của mô hình ngàm cứng FBB.

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

Luận văn là tài liệu học thuật và cẩm nang kỹ thuật giá trị dành cho 4 nhóm đối tượng chính:

  • Kỹ sư kết cấu và kỹ sư địa kỹ thuật công trình: Ứng dụng trực tiếp quy trình tính toán hệ số lò xo và giảm chấn tương đương 3 thành phần để mô hình hóa chính xác liên kết chân cột cho các dự án nhà cao tầng chịu động đất.
  • Học viên cao học và nghiên cứu sinh chuyên ngành Xây dựng: Tham khảo phương pháp nghiên cứu, cơ sở toán học của thuật toán Newmark-beta trên miền thời gian, và phương pháp xử lý số liệu phi tuyến địa kỹ thuật.
  • Giảng viên và nhà nghiên cứu tại các trường đại học: Sử dụng làm tài liệu giảng dạy chuyên sâu cho các môn học Động lực học công trình, Công trình chịu tải trọng đặc biệt, và Địa kỹ thuật động lực học nâng cao.
  • Chuyên gia thẩm tra và quản lý dự án xây dựng: Nắm vững bản chất cơ học của hiệu ứng tương tác đất - móng để đánh giá chính xác mức độ tin cậy của các báo cáo thẩm tra kháng chấn công trình phức tạp từ 9 đến 20 tầng.

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

1. Tại sao mô hình chân ngàm truyền thống (FBB) không còn đủ độ tin cậy trong thiết kế kháng chấn? Mô hình FBB giả định nền đất cứng tuyệt đối, bỏ qua biến dạng xoay và trượt của móng cọc. Thực tế cho thấy khi chịu động đất, đất nền bị mềm hóa làm chuyển vị đỉnh tăng thêm từ 15% đến 28%, nếu dùng FBB sẽ dẫn đến nguy cơ công trình bị va đập hoặc mất ổn định biến dạng hình học ngoài dự tính.

2. Điểm khác biệt cốt lõi giữa mô hình $SSI_{sta}$ và mô hình $SSI_{dyn}$ là gì? Mô hình $SSI_{sta}$ sử dụng các hệ số độ cứng lò xo móng cố định dựa trên tính chất cơ lý tĩnh của đất. Ngược lại, $SSI_{dyn}$ cập nhật liên tục sự suy giảm của mô đun cắt $G_s(t)$ theo bước thời gian $\Delta t = 0,02$ giây theo mức độ biến dạng trượt thực tế, phản ánh chính xác trạng thái làm việc phi tuyến của đất nền.

3. Việc xét đến SSI giúp giảm hay làm tăng chi phí xây dựng phần ngầm? Khi xét SSI, lực cắt đáy và mô men uốn chân cột được tính toán giảm từ 10% đến 22% nhờ hiệu ứng tiêu tán năng lượng cản của đất. Điều này giúp giảm đáng kể kích thước tiết diện đài móng, số lượng cọc và cốt thép chân cột, mang lại hiệu quả kinh tế rõ rệt mà vẫn đảm bảo an toàn chịu lực.

4. Thuật toán Newmark-beta có ưu điểm gì khi giải bài toán động lực học này? Thuật toán Newmark với các tham số gia tốc tuyến tính ($\gamma = 0,5; \beta = 1/6$) cho phép giải hệ phương trình vi phân chuyển động từng bước trên miền thời gian một cách ổn định, xử lý mượt mà sự thay đổi liên tục của ma trận độ cứng $[K(t)]$ và ma trận cản $[C(t)]$ của hệ lò xo móng.

5. Kết quả của luận văn có thể áp dụng cho các quy mô công trình khác không? Có thể áp dụng hoàn toàn. Quy trình tính toán 5 bước và mã nguồn thuật toán xây dựng trong luận văn được chuẩn hóa tổng quát, cho phép kỹ sư dễ dàng thay đổi thông số số tầng, số nhịp, đặc trưng địa tầng nhiều lớp và các loại cọc bê tông cốt thép khác nhau trong thực tế.

Kết luận

  • Luận văn đã giải quyết thành công bài toán xác định sự thay đổi độ cứng động lực học của hệ móng cọc theo thời gian thực dưới tác động của các trận động đất tiêu chuẩn.
  • Chứng minh định lượng sự suy giảm độ cứng móng từ 30% đến 55%, kéo theo sự gia tăng chu kỳ dao động tự nhiên và chuyển vị đỉnh công trình so với mô hình ngàm cứng truyền thống.
  • Làm rõ cơ chế giảm nội lực chân cột từ 10% đến 22% nhờ hiệu ứng cản bức xạ và cản vật liệu trong môi trường đất nền.
  • Đóng góp bộ công cụ tính toán số hoàn chỉnh viết bằng ngôn ngữ MATLAB, đóng vai trò cầu nối quan trọng giữa lý thuyết địa kỹ thuật động lực học và thực hành thiết kế kết cấu công trình kháng chấn.
  • Lộ trình nghiên cứu tiếp theo trong 1 đến 2 năm tới cần mở rộng phân tích phi tuyến vật liệu thân cọc kết hợp mô hình không gian 3 chiều để nâng cao hơn nữa độ chính xác trong công tác thiết kế xây dựng hiện đại.