Tổng quan về luận án

Nghiên cứu về ứng xử dài hạn của kết cấu bê tông cốt thép (BTCT) là một trong những chủ đề trọng tâm của cơ học công trình và kỹ thuật xây dựng hiện đại. Trong thực tế khai thác, các công trình chịu tác động thường trực của các chuyển vị cưỡng bức tại gối tựa liên kết—điển hình như hiện tượng chênh lệch lún giữa các đài móng cọc hoặc sự co ngắn không đồng đều giữa hệ thống cột và vách trong các tòa nhà siêu cao tầng (Gilbert & Ranzi, 2010). Dưới tác động của chuyển vị cưỡng bức duy trì dài hạn, trạng thái ứng suất - biến dạng trong kết cấu biến đổi phức tạp theo thời gian do sự xuất hiện đồng thời của hiện tượng từ biến (creep) và chùng ứng suất (stress relaxation). Luận án tiến sĩ kỹ thuật xây dựng của tác giả Nguyễn Mạnh Hùng với tiêu đề "Nghiên cứu thực nghiệm sự làm việc dài hạn của dầm bê tông cốt thép chịu chuyển vị cưỡng bức gối tựa" (Trường Đại học Xây dựng Hà Nội, 2021) là công trình nghiên cứu tiên phong tại Việt Nam đi sâu vào bản chất cơ học thực nghiệm và giải tích của hiện tượng này.

Khoảng trống nghiên cứu (research gap) trọng yếu mà luận án tập trung giải quyết xuất phát từ việc phần lớn các tiêu chuẩn thiết kế hiện hành (TCVN 5574:2018, ACI 318, Eurocode 2) cũng như các công trình nghiên cứu quốc tế trước đây chỉ xem xét từ biến của bê tông một cách độc lập hoặc gộp chung với biến dạng co ngót và tải trọng sử dụng. Rất ít công trình cô lập và định lượng được mức độ chùng ứng suất thuần túy do chuyển vị cưỡng bức gây ra trên cấu kiện dầm thực tế khi dầm đã chuyển sang trạng thái nứt. Các nghiên cứu quốc tế trước đó như của Jukka Jokela (1990) hay Koichi Maekawa et al. (2006) khảo sát dầm trong điều kiện đồng thời có tải trọng nhiệt và tải trọng sử dụng, dẫn đến sự chồng chéo giữa các thành phần biến dạng, chưa tách biệt được tác động cơ học thuần túy của chuyển vị gối tựa.

Luận án xác lập 3 câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. Câu hỏi nghiên cứu 1 (RQ1): Hiện tượng chùng ứng suất diễn ra theo quy luật thời gian nào trên cấu kiện dầm BTCT chịu chuyển vị cưỡng bức gối tựa khi chịu tác động của nứt và tỷ lệ cốt thép khác nhau?
    • Giả thuyết 1 (H1): Chùng ứng suất làm suy giảm đáng kể phản lực gối tựa và mô men uốn ban đầu, với tốc độ suy giảm lớn nhất tập trung trong giai đoạn đầu (72–100 giờ đầu) và tiệm cận ổn định theo thời gian.
  2. Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để thiết lập mối quan hệ hàm số chính xác giữa hệ số từ biến $\varphi(t, t_0)$ và hệ số già hóa $\chi(t, t_0)$ của bê tông trong điều kiện khí hậu tự nhiên tại Việt Nam?
    • Giả thuyết 2 (H2): Hệ số già $\chi(t, t_0)$ không phải là hằng số cố định mà biến thiên phi tuyến theo thời gian chất tải $t_0$, cấp độ bền bê tông và hệ số từ biến $\varphi(t, t_0)$, đóng vai trò hiệu chỉnh mô đun đàn hồi hiệu quả (AEMM).
  3. Câu hỏi nghiên cứu 3 (RQ3): Mô hình giải tích nào có khả năng dự báo chuẩn xác sự suy giảm phản lực gối tựa $R_B(t)$ có xét đến sự hình thành vùng nứt và độ cứng trung bình dọc trục dầm?
    • Giả thuyết 3 (H3): Việc tích hợp phương pháp mô đun đàn hồi hiệu quả điều chỉnh theo tuổi (AEMM) với giải thuật phân đoạn vùng nứt theo TCVN 5574:2018 sẽ tạo ra mô hình dự báo tương thích với sai số thực nghiệm dưới 10%.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng lý thuyết biến dạng nhớt đàn hồi tuyến tính (linear viscoelasticity) của bê tông, phương trình tích phân Volterra và phương pháp mô đun đàn hồi hiệu quả điều chỉnh theo tuổi (Age-Adjusted Effective Modulus Method - AEMM) do Trost và Bažant phát triển. Về mặt phạm vi, nghiên cứu triển khai chương trình thực nghiệm trên 12 cấu kiện dầm BTCT kích thước $100 \times 150 \times 800\text{ mm}$ với 3 nhóm hàm lượng cốt thép chịu kéo $\mu_1 = 1,163%$ (2Ø10), $\mu_2 = 1,675%$ (2Ø12), $\mu_3 = 2,279%$ (2Ø14) trong điều kiện môi trường tự nhiên tại Phòng thí nghiệm LAS-XD125 (Hà Nội). Đóng góp mang tính đột phá của luận án là việc chế tạo thành công hệ gá thí nghiệm độc lập cho phép loại trừ biến dạng co ngót, biến dạng nhiệt và trọng lượng bản thân, cung cấp bộ cơ sở dữ liệu thực nghiệm chuẩn mực đầu tiên về chùng ứng suất kết cấu BTCT tại Việt Nam.

Literature Review và Positioning

Bản chất biến dạng theo thời gian của bê tông được mô tả tổng quát thông qua phương trình phân rã biến dạng của Gilbert & Ranzi (2010): $$\varepsilon_c(t) = \varepsilon_{in}(t) + \varepsilon_{cr}(t) + \varepsilon_{sh}(t) + \varepsilon_{th}(t)$$ Trong đó $\varepsilon_{in}(t)$ là biến dạng đàn hồi tức thời, $\varepsilon_{cr}(t)$ là biến dạng từ biến, $\varepsilon_{sh}(t)$ là biến dạng co ngót và $\varepsilon_{th}(t)$ là biến dạng nhiệt. Về mặt cấu trúc vi mô, Feldman & Sereda (1968) đã chứng minh từ biến và chùng ứng suất phát sinh từ sự dịch chuyển và trượt của các lớp nước xen kẽ giữa các tấm gel canxi silicat hydrat (C-S-H) trong cấu trúc xốp mao quản của đá xi măng dưới tác dụng của gradient áp suất và ứng suất cơ học.

Trong y văn quốc tế, tổng quan nghiên cứu về ứng xử dài hạn chia thành hai trường phái và nhánh lý thuyết chính:

  1. Nhánh lý thuyết từ biến và mô hình hóa tích phân: Bažant & Baweja (1995) với mô hình B3, ACI 209R và Eurocode 2 (EC2: 2004) tiếp cận từ biến theo hướng hàm từ biến riêng $J(t, t_0)$ và nguyên lý xếp chồng Boltzmann. Ngược lại, phương pháp gần đúng đại số như Mô đun đàn hồi hiệu quả EMM (Effective Modulus Method) và Mô đun hiệu quả điều chỉnh theo tuổi AEMM của Trost (1967) và Bažant (1972) đưa vào hệ số già $\chi(t, t_0)$ nhằm đơn giản hóa việc giải phương trình vi phân tích phân mà vẫn duy trì độ chính xác cao.
  2. Nhánh nghiên cứu thực nghiệm chùng ứng suất trên vật liệu và cấu kiện: Morimoto & Koyanagi (1994) cùng Atrushi (2003) nghiên cứu chùng ứng suất trên mẫu lăng trụ bê tông chịu nén và kéo đơn trục, chỉ ra rằng mức độ chùng ứng suất phụ thuộc chặt chẽ vào tỷ lệ nước/xi măng ($N/X$), tuổi chất tải $t_0$ và trạng thái ứng suất. Tuy nhiên, các thí nghiệm này chỉ dừng lại ở quy mô mẫu vật liệu đẳng chất.

Một cuộc tranh luận khoa học lớn tồn tại giữa các nhà nghiên cứu: Serda & Kristek (1988) cho rằng mối quan hệ giữa ứng suất ban đầu và chùng ứng suất mang tính phi tuyến rõ rệt khi tỷ lệ ứng suất vượt quá $0,4 f_c$, trong khi Gutsch & Rostasy (1994) lập luận rằng trong dải ứng suất làm việc thông thường ($0,3–0,5 f_c$), chùng ứng suất có thể coi là tuyến tính nhưng bị chi phối mạnh mẽ bởi các vết nứt vi mô (microcracks). Hossain & Weiss (2004) bổ sung bằng chứng cho thấy mức độ kiềm chế biến dạng càng cao thì mật độ nứt vi mô càng dày đặc, thúc đẩy tốc độ chùng ứng suất diễn ra nhanh hơn.

Khi so sánh với các nghiên cứu quốc tế trên cấu kiện kết cấu thực tế:

  • Nghiên cứu của Jukka Jokela (1990): Thực hiện thí nghiệm trên 4 mẫu dầm BTCT chịu đồng thời chuyển vị cưỡng bức gối tựa và gradient nhiệt độ. Nghiên cứu này ghi nhận sự suy giảm phản lực gối tựa nhưng không thể cô lập phần suy giảm do chùng ứng suất thuần túy và phần do suy giảm mô đun vì nhiệt độ cao.
  • Nghiên cứu của Koichi Maekawa et al. (2006): Khảo sát dầm liên tục hai nhịp chịu chuyển vị cưỡng bức gối giữa kết hợp tải trọng phân bố sử dụng. Kết quả cho thấy phản lực gối giảm theo thời gian, song tương tác giữa mô men uốn do ngoại tải và nội lực do cưỡng bức biến dạng làm lu mờ ảnh hưởng của tỷ lệ hàm lượng cốt thép đối với quá trình hồi phục độ cứng.

Luận án của NCS. Nguyễn Mạnh Hùng định vị chính xác khoảng trống này bằng cách thiết kế hệ thí nghiệm dầm một đầu ngàm, một đầu khớp, ứng dụng hệ kích nêm khóa cứng chuyển vị tại một thời điểm xác định và thiết lập hệ đo biến dạng khử co ngót - nhiệt độ song song. Cách tiếp cận này giúp cô lập tuyệt đối biến dạng chuyển vị cưỡng bức, từ đó đóng góp một bước tiến thực nghiệm mới cho kho tàng khoa học kết cấu thế giới.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng và hoàn thiện lý thuyết biến dạng nhớt đàn hồi áp dụng cho cấu kiện dầm BTCT làm việc ở trạng thái nứt (Stage II cracked section) chịu điều kiện biên chuyển vị cưỡng bức không đổi:

  1. Mở rộng lý thuyết AEMM cho cấu kiện nứt: Tác giả mở rộng phương pháp mô đun đàn hồi hiệu quả điều chỉnh theo tuổi (AEMM) của Bažant bằng cách đưa hàm độ cứng uốn hiệu quả $B_{eff}(t, t_0)$ vào phương trình phân bố nội lực dài hạn, thay vì chỉ áp dụng cho vật liệu đồng nhất hoặc mặt cắt nguyên đàn hồi: $$E_{eff}(t, t_0) = \frac{E_c(t_0)}{1 + \chi(t, t_0) \varphi(t, t_0)}$$
  2. Thiết lập mối quan hệ giải tích giữa $\chi(t, t_0)$ và $\varphi(t, t_0)$: Dựa trên phương trình tích phân liên hệ giữa hàm từ biến $J(t, t')$ và hàm chùng ứng suất $R(t, t')$: $$\int_{t_0}^t R(\tau, t') \frac{\partial J(t, \tau)}{\partial \tau} d\tau = 1$$ Luận án đã số hóa và giải phương trình tích phân Volterra loại 2 trên môi trường Matlab, từ đó xây dựng biểu đồ và bảng tra hệ số già $\chi(t, t_0)$ theo hệ số từ biến $\varphi(t, t_0)$ ứng với cấp độ bền bê tông từ B20 đến B35 theo tiêu chuẩn Việt Nam.
  3. Mô hình hóa sự dịch chuyển trục trung hòa và phân phối lại ứng suất: Luận án chứng minh rằng dưới tác động của chùng ứng suất, ứng suất kéo trong cốt thép $\sigma_s(t)$ và ứng suất nén trong bê tông $\sigma_b(t)$ suy giảm đồng thời theo thời gian, làm chậm lại quá trình phát triển bề rộng vết nứt mà không làm suy giảm khả năng chịu lực cực hạn của dầm.

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đồng bộ 3 lý thuyết nền tảng: (1) Lý thuyết cơ học chất rắn biến dạng nhớt đàn dẻo, (2) Lý thuyết cấu kiện thanh BTCT chịu uốn có nứt theo mô hình trạng thái giới hạn (TCVN 5574:2018), và (3) Phương pháp hàm chùng ứng suất xấp xỉ ARF (Approximate Relaxation Function) theo hướng dẫn của CEB-FIP Model Code (MC 90-93).

+-------------------------------------------------------------------------------+
|                        KHUNG PHÂN TÍCH TOÀN DIỆN CỦA LUẬN ÁN                   |
+-------------------------------------------------------------------------------+
|  1. ĐIỀU KIỆN BIÊN CƠ HỌC: Chuyển vị cưỡng bức gối tựa cố định Delta(t) = Delta_0 |
+-------------------------------------------------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
|  2. MÔ HÌNH HÓA VẬT LIỆU VISCOELASTIC:                                         |
|     - Hàm từ biến J(t, t0) và Hệ số từ biến phi(t, t0) (Eurocode 2 / TCVN)     |
|     - Giải thuật số giải hàm chùng ứng suất ER(t, t0) qua tích phân Volterra   |
|     - Xác lập hệ số già hóa chi(t, t0) -> Mô đun hiệu quả E_eff(t, t0)          |
+-------------------------------------------------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
|  3. PHÂN TÍCH TIẾT DIỆN VÀ CHIỀU DÀI PHÂN TỬ NỨT:                             |
|     - Phân chia dầm thành 2 phân vùng: Vùng chưa nứt (I_red) & Vùng nứt (D_cr)|
|     - Tính toán chiều dài vùng nứt l_crc(t) và độ cứng trung bình D_tb(t)     |
+-------------------------------------------------------------------------------+
                                        |
                                        v
+-------------------------------------------------------------------------------+
|  4. XÁC LẬP MÔ HÌNH DỰ BÁO SUY GIẢM PHẢN LỰC GỐI TỰA:                        |
|     - Phương trình cân bằng độ võng đàn hồi - nhớt: Delta_B(t) = const        |
|     - Dự báo tỷ lệ suy giảm phản lực: RB(t)/RB(t0) = f(E_eff, mu, D_tb, t)    |
+-------------------------------------------------------------------------------+

Điều kiện biên của mô hình được xác định chặt chẽ: Ứng suất nén trong bê tông tại thời điểm tạo chuyển vị không vượt quá $0,5 f_c$ nhằm bảo toàn tính tuyến tính của biến dạng từ biến; cấu kiện dầm liên kết một đầu ngàm - một đầu khớp chịu chuyển vị cưỡng bức tại gối khớp theo phương thẳng đứng $\Delta_0$.

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

Thiết kế nghiên cứu

Luận án tuân thủ thế giới quan thực chứng (positivism) với phương pháp luận diễn dịch kết hợp chặt chẽ giữa suy diễn giải tích toán học và kiểm định thực nghiệm định lượng. Thiết kế nghiên cứu đa cấp độ (multi-level design) bao gồm:

  • Cấp độ vật liệu: Thử nghiệm xác định cường độ nén mẫu lập phương $150 \times 150 \times 150\text{ mm}$, mô đun đàn hồi $E_c$, hệ số từ biến bê tông $\varphi(t, t_0)$ trên thiết bị nén lò xo chuyên dụng giữ tải không đổi, và thí nghiệm đo chùng ứng suất thuần túy của thanh thép liên kết Ø6.
  • Cấp độ cấu kiện kết cấu: Khảo sát 12 mẫu dầm BTCT chia thành 3 nhóm theo hàm lượng cốt thép dọc chịu kéo ($\mu_1 = 1,163%$, $\mu_2 = 1,675%$, $\mu_3 = 2,279%$), mỗi nhóm gồm 4 dầm thử nghiệm lặp lại để loại trừ sai số ngẫu nhiên. Kích thước hình học tiêu chuẩn: $b \times h \times L = 100 \times 150 \times 800\text{ mm}$.

Quy trình nghiên cứu rigorous

Để đạt được độ chính xác tuyệt đối trong việc đo đạc chùng ứng suất, tác giả đã sáng tạo một hệ thống thí nghiệm phức hợp gồm các cơ cấu chuyên biệt:

  1. Hệ khung trượt và nêm thép khóa cứng: Tạo chuyển vị cưỡng bức $\Delta_0$ tức thời tại đầu gối tựa tự do bằng kích thủy lực, sau đó đóng hệ nêm thép định vị để khóa cứng độ võng tại gối tựa trong suốt thời gian theo dõi dài hạn.
  2. Hệ khung treo và thanh chặn chống xoay ngang: Triệt tiêu hoàn toàn tác động của trọng lượng bản thân dầm lên mô men uốn và ngăn ngừa sự mất ổn định ngoài mặt phẳng uốn.
  3. Mẫu dầm đối chứng khử biến dạng co ngót và nhiệt độ: Bố trí một mẫu dầm có cùng cấp phối và kích thước hình học không chịu chuyển vị cưỡng bức nằm kế bên, trang bị đầy đủ cảm biến đo biến dạng co ngót $\varepsilon_{sh}(t)$ và biến dạng nhiệt $\varepsilon_{th}(t)$. Giá trị biến dạng cơ học thực tế được hiệu chỉnh bằng phương pháp trừ trực tiếp: $$\varepsilon_{mech}(t) = \varepsilon_{total}(t) - \left[ \varepsilon_{sh}(t) + \varepsilon_{th}(t) \right]$$
  4. Hệ thống đo phản lực thông qua thanh thép liên kết Ø6: Phản lực gối tựa được đo gián tiếp và liên tục thông qua cảm biến đo biến dạng (strain gauge) gắn trên thanh treo thép Ø6 đã được hiệu chuẩn chính xác mối quan hệ lực - biến dạng trong miền đàn hồi.
+-------------------------------------------------------------------------------+
|               SƠ ĐỒ BỐ TRÍ THÍ NGHIỆM ĐO CHÙNG ỨNG SUẤT DẦM BTCT              |
+-------------------------------------------------------------------------------+
|                                                                               |
|      [NGÀM CỐ ĐỊNH]                                [GỐI TỰA DI ĐỘNG CƯỠNG BỨC]|
|     +--------------+                              +-----------------------+   |
|     |  Ngàm thép   |==============================| Thanh thép liên kết Ø6|   |
|     |  khóa cứng   |        DẦM BTCT THỬ NGHIỆM   | gắn Strain Gauge đo RB|   |
|     +--------------+  (100 x 150 x 800 mm)        +-----------------------+   |
|            |                                                  |               |
|            |====== [Khung treo dầm triệt tiêu trọng lượng] ===|               |
|            |                                                  v               |
|            |                                      [Kích nêm khóa chuyển vị]   |
|                                                                               |
|     ------------------------------------------------------------------------- |
|      [DẦM ĐỐI CHỨNG KHỬ CO NGÓT & NHIỆT ĐỘ]                                   |
|      (Cùng kích thước, không cưỡng bức, gắn cảm biến đo eps_sh và eps_th)     |
+-------------------------------------------------------------------------------+

Data và phân tích

Số liệu thực nghiệm được thu thập định kỳ theo thang thời gian logarit: đo dày đặc trong 24 giờ đầu (1h, 2h, 4h, 6h, 12h, 24h), giãn cách theo ngày trong tuần đầu và theo chu kỳ tuần/tháng trong các giai đoạn tiếp theo kéo dài hàng trăm ngày.

  • Công cụ phân tích: Xây dựng thuật toán lặp số trên phần mềm Matlab (Phụ lục 8) để tính toán hàm chùng ứng suất $E_R(t, t_0)$ từ dữ liệu từ biến đo được, và thiết lập mô hình tính toán giải tích sự suy giảm phản lực gối tựa bằng phần mềm Maple (Phụ lục 9).
  • Kiểm định độ tin cậy: Tính toán độ tin cậy thông qua hệ số tương quan ($R^2 > 0,92$) giữa đường cong suy giảm phản lực gối tựa lý thuyết và thực nghiệm trên cả 3 nhóm mẫu dầm D1-10, D1-12 và D1-14.

Phát hiện đột phá và implications

Những phát hiện then chốt

Dữ liệu thực nghiệm từ luận án đã chỉ ra các quy luật cơ học đột phá về sự làm việc dài hạn của dầm BTCT chịu chuyển vị cưỡng bức gối tựa:

  1. Mức độ suy giảm phản lực gối tựa sâu sắc: Phản lực gối tựa $R_B(t)$ suy giảm mạnh mẽ theo thời gian do hiện tượng chùng ứng suất. Tỷ lệ suy giảm $R_B(t)/R_B(t_0)$ sau thời gian theo dõi đạt từ 20% đến 35% so với phản lực ban đầu tại thời điểm $t_0$.
  2. Quy luật động học 3 giai đoạn của chùng ứng suất:
    • Giai đoạn 1 (0–72 giờ): Tốc độ chùng ứng suất diễn ra cực nhanh, chiếm khoảng 40–50% tổng lượng suy giảm ứng suất dài hạn.
    • Giai đoạn 2 (3–30 ngày): Tốc độ suy giảm giảm dần rõ rệt theo quy luật hàm mũ.
    • Giai đoạn 3 (sau 30 ngày): Quá trình suy giảm diễn ra chậm chạp và tiệm cận dần tới giá trị ổn định lâu dài ($R_{B,\infty}$).
  3. Ảnh hưởng phi tuyến của hàm lượng cốt thép $\mu$: Nhóm dầm có hàm lượng cốt thép cao ($\mu_3 = 2,279%$, dầm D-14) có tỷ lệ suy giảm phản lực gối tựa $R_B(t)/R_B(t_0)$ thấp hơn so với nhóm dầm có hàm lượng cốt thép thấp ($\mu_1 = 1,163%$, dầm D-10). Nguyên nhân là do cốt thép không bị từ biến, đóng vai trò khung chịu lực giữ ổn định độ cứng của tiết diện nứt, làm giảm tỷ phần đóng góp của từ biến bê tông vào độ chùng chung của kết cấu.
  4. Hiệu ứng giải tỏa ứng suất và kiểm soát nứt: Mặc dù chuyển vị gối tựa $\Delta_0$ được duy trì cố định, ứng suất kéo trong cốt thép tại tiết diện ngàm $\sigma_s(t)$ giảm dần theo thời gian đồng nhịp với sự suy giảm phản lực gối tựa. Hiện tượng này khẳng định chùng ứng suất đóng vai trò tích cực trong việc hạn chế sự mở rộng của các vết nứt uốn dọc theo chiều dài dầm.
Tỷ lệ suy giảm phản lực RB(t)/RB(t0)
1.0 |====================================================
    | \
0.9 |  \  (Nhóm mu3 = 2.279% - Cốt thép cao: Giảm ~22%)
    |   \------------------------------------------------
0.8 |    \
    |     \  (Nhóm mu2 = 1.675% - Cốt thép TB: Giảm ~28%)
0.7 |      \---------------------------------------------
    |       \
0.6 |        \  (Nhóm mu1 = 1.163% - Cốt thép thấp: Giảm ~34%)
    |         \------------------------------------------
0.0 +----------------------------------------------------> Thời gian (t)
    0h   72h     7 ngày     30 ngày       180 ngày

Implications đa chiều

  • Ý nghĩa học thuật: Khẳng định tính chính xác của lý thuyết AEMM khi áp dụng cho bài toán thanh nứt chịu biến dạng cưỡng bức; làm rõ sự tương hỗ hai chiều giữa từ biến và chùng ứng suất trong kết cấu phức hợp bê tông - cốt thép.
  • Ý nghĩa phương pháp luận: Chuẩn hóa quy trình thực nghiệm đo chùng ứng suất kết cấu thực tại Việt Nam, đặc biệt là giải pháp kỹ thuật triệt tiêu ảnh hưởng của co ngót tự do, biến thiên nhiệt độ môi trường và trọng lượng bản thân.
  • Ứng dụng thực tiễn kỹ thuật: Cung cấp cơ sở khoa học để kỹ sư thiết kế tối ưu hóa tính toán kết cấu móng bè - cọc và hệ dầm giằng móng nhà cao tầng. Trong tính toán truyền thống, phản lực do lún lệch thường được coi là không đổi theo thời gian dẫn đến việc bố trí cốt thép quá mức an toàn cục bộ. Việc kể đến hệ số suy giảm chùng ứng suất (0,70–0,80) cho phép tiết kiệm 15–20% hàm lượng cốt thép gối mà vẫn đảm bảo an toàn chịu lực và độ bền lâu.
  • Khuyến nghị chính sách và tiêu chuẩn: Đề xuất bổ sung chỉ dẫn tính toán chùng ứng suất và hệ số già hóa $\chi(t, t_0)$ vào các phiên bản cập nhật tiếp theo của tiêu chuẩn TCVN 5574.

Limitations và Future Research

Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn thừa nhận 4 giới hạn nghiên cứu:

  1. Quy mô kích thước mẫu dầm: Thí nghiệm được thực hiện trên mẫu dầm kích thước nhỏ ($L = 800\text{ mm}$, tiết diện $100 \times 150\text{ mm}$) do giới hạn tải trọng của hệ khung nêm gá đặt trong phòng thí nghiệm; hiệu ứng kích thước (size effect) đối với sự phát triển từ biến và chùng ứng suất cần được kiểm chứng thêm trên dầm nhịp lớn.
  2. Cấp độ bền bê tông đơn lẻ: Thực nghiệm tập trung chủ yếu vào cấp độ bền bê tông B25 ($f_c \approx 30\text{ MPa}$); chưa khảo sát thực nghiệm đối với bê tông cường độ cao (HPC, $f_c > 60\text{ MPa}$) với tỷ lệ $N/X$ thấp và phụ gia khoáng silica fume.
  3. Môi trường tự nhiên không điều hòa: Quá trình đo diễn ra trong điều kiện nhiệt - ẩm tự nhiên của phòng thí nghiệm tại Hà Nội, dẫn đến dao động nhiệt ẩm theo mùa dù đã được khử trừ bằng mẫu đối chứng.
  4. Quy luật chuyển vị cưỡng bức một lần: Nghiên cứu chỉ khảo sát chuyển vị cưỡng bức tức thời đặt một lần $\Delta_0 = \text{const}$, chưa xét đến bài toán chuyển vị gối tựa tăng dần theo thời gian (như quá trình lún cố kết dài hạn của nền đất yếu).

Định hướng nghiên cứu tương lai (Future Agenda):

  • Phát triển mô hình phần tử hữu hạn phi tuyến (Nonlinear FEM) mô phỏng tương tác dầm - nền đất có xét đến tốc độ lún cố kết $\Delta(t)$ biến thiên theo thời gian.
  • Mở rộng chương trình thực nghiệm trên kết cấu dầm BTCT dự ứng lực (căng trước và căng sau) chịu chuyển vị cưỡng bức gối tựa.
  • Nghiên cứu ứng xử chùng ứng suất của bê tông tính năng siêu cao (UHPC) và bê tông geopolymer thân thiện môi trường.
  • Xây dựng mô-đun phần mềm chuyên dụng tích hợp vào các phần mềm kết cấu thương mại (như ETABS, SAP2000) để tự động hóa tính toán suy giảm nội lực do lún lệch.

Tác động và ảnh hưởng

  • Tác động học thuật: Công trình tạo ra một bước nhảy vọt về dữ liệu thực nghiệm tại khu vực Đông Nam Á, nơi có điều kiện khí hậu nhiệt đới ẩm gió mùa đặc thù ảnh hưởng mạnh đến tốc độ phát triển từ biến bê tông. Các công bố khoa học từ luận án trên các tạp chí chuyên ngành uy tín trong nước và hội nghị khoa học đóng vai trò tài liệu tham khảo cốt lõi cho các nghiên cứu sau đại học về độ bền lâu kết cấu.
  • Chuyển đổi ngành công nghiệp xây dựng: Giúp các tổng thầu thiết kế, tư vấn thẩm tra thay đổi tư duy tính toán truyền thống về nội lực do lún lệch đài móng và co ngắn cột - vách trong các dự án nhà chọc trời (trên 40 tầng), trực tiếp giảm chi phí vật liệu thép xây dựng và giảm phát thải carbon ngành xây dựng.
  • Ảnh hưởng chính sách và quy chuẩn: Đóng góp cơ sở thực nghiệm vững chắc cho Viện Khoa học Công nghệ Xây dựng (IBST) và Bộ Xây dựng trong quá trình biên soạn, hoàn thiện các quy chuẩn thiết kế kết cấu bê tông cốt thép tiệm cận với tiêu chuẩn quốc tế Eurocode 2 và fib Model Code 2010.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận thiết kế thí nghiệm độc đáo trong việc khử các thành phần biến dạng ký sinh (co ngót, nhiệt độ, trọng lượng bản thân) và giải thuật số phương trình tích phân Volterra.
  • Các nhà khoa học và giảng viên chuyên ngành công trình: Khai thác tập dữ liệu đo đạc thực nghiệm chi tiết về từ biến và chùng ứng suất theo thời gian phục vụ hiệu chỉnh và kiểm chứng các mô hình toán học giải tích mới.
  • Kỹ sư kết cấu và chuyên gia R&D: Ứng dụng trực tiếp công thức giải tích dự báo tỷ lệ suy giảm phản lực gối tựa $R_B(t)/R_B(t_0)$ vào thiết kế dầm móng và dầm giằng liên kết cột vách chịu lún lệch, nâng cao hiệu quả kinh tế - kỹ thuật của công trình.
  • Cơ quan quản lý và hội đồng biên soạn tiêu chuẩn: Có luận cứ khoa học để định lượng các hệ số suy giảm nội lực dài hạn trong hệ thống tiêu chuẩn quốc gia TCVN.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc mở rộng phương pháp mô đun đàn hồi hiệu quả điều chỉnh theo tuổi (AEMM) của Trost (1967) và Bažant (1972) cho cấu kiện dầm BTCT làm việc ở trạng thái nứt chịu chuyển vị cưỡng bức gối tựa. Luận án đã tích hợp sự biến thiên của chiều dài vùng nứt $l_{crc}(t)$ và độ cứng trung bình $D_{tb}(t)$ dọc trục thanh với hàm chùng ứng suất $E_R(t, t_0)$, chuyển đổi thành công bài toán giải tích phi tuyến phức tạp thành mô hình tính toán đại số tường minh giúp xác định chính xác sự suy giảm phản lực gối tựa $R_B(t)$.

2. Phương pháp thực nghiệm của luận án có gì đổi mới so với ít nhất 2 nghiên cứu quốc tế trước đó?

So với nghiên cứu của Jukka Jokela (1990) (thí nghiệm dầm chịu đồng thời tải trọng nhiệt và chuyển vị cưỡng bức) và Koichi Maekawa et al. (2006) (thí nghiệm dầm 2 nhịp chịu tải trọng sử dụng kết hợp chuyển vị gối), nghiên cứu của NCS. Nguyễn Mạnh Hùng có bước đột phá vượt bậc về kỹ thuật thí nghiệm:

  • Thiết kế hệ khung treo dầm khử triệt để ảnh hưởng của trọng lượng bản thân.
  • Bố trí mẫu dầm đối chứng song song để đo và trừ khử trực tiếp biến dạng co ngót $\varepsilon_{sh}(t)$ và biến dạng nhiệt $\varepsilon_{th}(t)$.
  • Sử dụng hệ kích nêm khóa cứng biến dạng cưỡng bức tức thời $\Delta_0$, tạo ra điều kiện biên chùng ứng suất thuần túy chuẩn mực mà chưa nghiên cứu quốc tế nào trước đó thực hiện được trên cấu kiện dầm nứt.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ là gì?

Phát hiện bất ngờ nhất là tốc độ chùng ứng suất diễn ra cực kỳ nhanh trong 72 giờ đầu tiên, làm suy giảm tới 40–50% tổng lượng phản lực suy giảm của toàn bộ chu kỳ dài hạn, sau đó đường cong suy giảm chuyển pha rõ rệt sang trạng thái ổn định tiệm cận. Dữ liệu thực nghiệm trên nhóm dầm $\mu_1 = 1,163%$ cho thấy tỷ lệ $R_B(t)/R_B(t_0)$ giảm mạnh ngay trong 3 ngày đầu và đạt mức suy giảm tổng cộng xấp xỉ 34% sau thời gian quan trắc, minh chứng cho khả năng tự điều chỉnh và giải tỏa ứng suất nội tại rất mạnh mẽ của bê tông nứt.

4. Luận án có cung cấp quy trình nhân bản thí nghiệm (Replication Protocol) hoàn chỉnh không?

Luận án cung cấp quy trình nhân bản thí nghiệm chi tiết và hoàn chỉnh tuyệt đối:

  • Bản vẽ cấu tạo chi tiết hệ khung trượt tạo chuyển vị, hệ nêm khóa cứng, hệ khung treo và thanh chặn chống xoay ngang (Chương 3).
  • Chi tiết cấp phối bê tông, quy trình bảo dưỡng và thông số thử nghiệm cơ lý vật liệu (cường độ nén, mô đun đàn hồi, biến dạng từ biến trên lò xo chuyên dụng).
  • Thuật toán và mã nguồn chi tiết thiết lập tính lặp trên phần mềm Matlab để tính hàm chùng $E_R(t, t_0)$ (Phụ lục 8) và giải thuật Maple tính toán phản lực gối tựa (Phụ lục 9).

5. Chương trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Chương trình nghiên cứu 10 năm tiếp theo tập trung vào 3 trục chính:

  1. Mở rộng quy mô và vật liệu mới (Năm 1–3): Thử nghiệm chùng ứng suất trên dầm kích thước thật ($L \ge 6\text{ m}$) sử dụng bê tông cường độ siêu cao (UHPC) và bê tông cốt sợi phân tán (FRC).
  2. Tương tác động lực học dầm - nền đất thời gian thực (Năm 4–6): Nghiên cứu chuyển vị cưỡng bức gối tựa biến thiên theo quy luật lún cố kết phi tuyến của đất nền dưới tác động của tải trọng trùng phục và động đất.
  3. Số hóa và tiêu chuẩn hóa (Năm 7–10): Xây dựng bộ tiêu chuẩn thiết kế kháng lún lệch công trình BTCT dựa trên lý thuyết chùng ứng suất và tích hợp gói công cụ phân tích tự động vào các phần mềm thiết kế kết cấu BIM/FEM tiêu chuẩn quốc tế.

Kết luận

  1. Xây dựng thành công mô hình thực nghiệm chuẩn mực: Luận án đã thiết kế và vận hành thành công hệ thống thí nghiệm dầm BTCT một đầu ngàm một đầu khớp chịu chuyển vị cưỡng bức gối tựa, loại trừ hoàn toàn ảnh hưởng của trọng lượng bản thân, co ngót và nhiệt độ môi trường.
  2. Xác lập quy luật suy giảm phản lực gối tựa định lượng: Chứng minh thực nghiệm rằng chùng ứng suất làm giảm từ 20% đến 35% phản lực gối tựa và mô men uốn ban đầu của dầm BTCT, trong đó giai đoạn 72 giờ đầu đóng vai trò chi phối tốc độ suy giảm.
  3. Định lượng ảnh hưởng của hàm lượng cốt thép: Khẳng định hàm lượng cốt thép $\mu$ tỷ lệ nghịch với mức độ suy giảm phản lực gối tựa; cốt thép đóng vai trò khung giằng đàn hồi giữ ổn định độ cứng dài hạn của cấu kiện.
  4. Thiết lập bảng tra hệ số già hóa $\chi(t, t_0)$ cho điều kiện Việt Nam: Giải thành công phương trình tích phân Volterra liên hệ giữa hàm từ biến và hàm chùng ứng suất, cung cấp công cụ tra cứu thực dụng cho kỹ sư thiết kế.
  5. Xây dựng công thức giải tích dự báo chùng ứng suất chính xác: Đề xuất công thức tính toán suy giảm phản lực gối tựa $R_B(t)$ trên nền tảng phương pháp AEMM có xét đến độ cứng vùng nứt trung bình $D_{tb}(t)$, đạt độ tương đồng cao với kết quả thực nghiệm ($R^2 > 0,92$).
  6. Đặt nền móng cho trường phái thiết kế tối ưu kết cấu theo độ bền lâu: Chuyển đổi tư duy thiết kế từ việc xem nội lực cưỡng bức là hằng số sang đại lượng suy giảm theo thời gian, mở ra hướng đi mới tiết kiệm tài nguyên và nâng cao tính an toàn bền vững cho các công trình xây dựng hiện đại.