Luận án tiến sĩ: Nghiên cứu thực nghiệm dầm BTCT chịu chuyển vị cưỡng bức gối tựa

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 bài toán phức tạp và then chốt nhất của cơ học công trình hiện đại. Bản chất của bê tông là vật liệu đàn dẻo nhớt (viscoelastic-plastic), trong đó tổng biến dạng $\varepsilon_c(t)$ theo thời gian chịu sự chi phối đồng thời của nhiều thành phần: biến dạng tức thời đàn hồi $\varepsilon_{in}(t)$, biến dạng từ biến $\varepsilon_{cr}(t)$, biến dạng co ngót $\varepsilon_{sh}(t)$ và biến dạng nhiệt $\varepsilon_{th}(t)$, được mô tả qua phương trình tương thích biến dạng tổng quát của Gilbert và Ranzi (2010):

$$\varepsilon_c(t) = \varepsilon_{in}(t) + \varepsilon_{cr}(t) + \varepsilon_{sh}(t) + \varepsilon_{th}(t)$$

Trong các công trình thực tế, kết cấu dầm BTCT thường xuyên chịu tác động của chuyển vị cưỡng bức gối tựa được duy trì lâu dài. Các kịch bản điển hình bao gồm: dầm móng liên kết các đài cọc chịu chuyển vị do lún lệch địa kỹ thuật, hoặc dầm giằng liên kết giữa hệ cột và vách cứng trong nhà cao tầng chịu chuyển vị cưỡng bức do hiện tượng co ngắn không đều (differential axial shortening) giữa các cấu kiện thẳng đứng. Dưới tác động của chuyển vị cưỡng bức duy trì cố định, hiện tượng chùng ứng suất (stress relaxation) xảy ra đồng thời với từ biến, làm suy giảm nội lực và phản lực gối tựa theo thời gian.

Khoảng trống nghiên cứu (research gap) cốt lõi hiện nay là: trong khi hiện tượng từ biến chịu tải trọng tĩnh đã được quy chuẩn hóa trong các tiêu chuẩn thiết kế như TCVN 5574:2018, ACI 318:2014 và Eurocode 2 (EC2:2004), thì hiện tượng chùng ứng suất trên cấu kiện dầm BTCT chịu uốn lại chưa có mô hình dự báo hoàn chỉnh. Đa số các nghiên cứu quốc tế trước đây (Jukka Jokela, 1990; Koichi Maekawa et al., 2006) khảo sát chùng ứng suất trong điều kiện chịu tải trọng hỗn hợp (tải trọng sử dụng kết hợp biến dạng co ngót), dẫn đến sự chồng lấn các hiệu ứng cơ học và không thể cô lập để định lượng chính xác quy luật chùng ứng suất thuần túy do chuyển vị cưỡng bức gối tựa.

Luận án tiến sĩ kỹ thuật xây dựng (mã số 9580201) của nghiên cứu sinh Nguyễn Mạnh Hùng (2021) tại Trường Đại học Xây dựng, dưới sự hướng dẫn của PGS.TS. Nguyễn Trung Hiếu và GS.TS. Ngô Thế Phong, tập trung giải quyết 3 câu hỏi nghiên cứu (Research Questions - RQ) và kiểm chứng 2 giả thuyết khoa học (Hypotheses - H):

• RQ1: Quy luật suy giảm phản lực gối tựa và phân phối lại ứng suất trong dầm BTCT chịu chuyển vị cưỡng bức diễn ra như thế nào theo thời gian dưới các hàm lượng cốt thép khác nhau?
• RQ2: Mối quan hệ tương hỗ 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 nhiệt đới ẩm được xác lập ra sao để phục vụ tính toán dài hạn?
• RQ3: Khung giải tích nào cho phép dự báo chính xác tốc độ suy giảm phản lực gối tựa $R_B(t)$ có xét đến sự suy giảm độ cứng và phân vùng nứt của dầm BTCT?
• H1: Hiện tượng 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 dài hạn (từ 25% đến 45%), đóng vai trò có lợi trong việc kiểm soát ứng suất kéo trong cốt thép và hạn chế bề rộng khe nứt bê tông.
• H2: Hàm lượng cốt thép dọc chịu kéo $\mu$ chi phối trực tiếp đến độ cứng tương đương của tiết diện nứt $I_{red}$, làm thay đổi tốc độ và giá trị chùng ứng suất cuối cùng của dầm.

Phạm vi thực nghiệm của luận án bao gồm 12 mẫu dầm BTCT kích thước $100 \times 150 \times 800\text{ mm}$, chia thành 3 nhóm hàm lượng cốt thép: $\mu_1 = 1{,}163%$ ($2\Phi10$), $\mu_2 = 1{,}675%$ ($2\Phi12$), và $\mu_3 = 2{,}279%$ ($2\Phi14$), được chế tạo từ bê tông nặng cấp độ bền B25–B30 ($f_c \approx 30\text{ MPa}$), thí nghiệm tại Phòng thí nghiệm LAS-XD125 trong điều kiện môi trường tự nhiên tại Hà Nội.

Literature Review và Positioning

Nghiên cứu về cơ chế vi mô của từ biến và chùng ứng suất khởi nguồn từ mô hình cấu trúc lực dính xi măng của Feldman và Sereda (1968), giải thích rằng biến dạng phụ thuộc thời gian bắt nguồn từ sự trượt nhớt của các tấm gel canxi silicat hydrat (C-S-H) và sự khuếch tán của các phân tử nước xen kẽ giữa các lớp mao quản dưới tác dụng của gradient áp suất. Gilbert và Ranzi (2010) cùng Alexander và Beushausen (2009) đã hệ thống hóa từ biến thành hai thành phần chính: từ biến cơ bản (basic creep) xảy ra khi không có trao đổi ẩm và từ biến khô (drying creep) do mất nước ra môi trường bên ngoài.

Trong lĩnh vực chùng ứng suất vật liệu, các tranh luận học thuật lớn đã diễn ra xung quanh tính tuyến tính và các yếu tố ảnh hưởng:

1. Tranh luận về mức ứng suất ban đầu: Morimoto và Koyanagi (1994) cho rằng mối quan hệ giữa ứng suất ban đầu và mức độ chùng ứng suất là tuyến tính khi tỷ số ứng suất/cường độ đạt tới 80%. Ngược lại, Gutsch và Rostasy (1994) cùng Serda và Kristek (1988) chứng minh rằng khi ứng suất vượt quá 11,8 MPa hoặc đạt trên 50% cường độ chịu nén đặc trưng $0{,}5f'_c$, các vết nứt vi mô (microcracks) xuất hiện trong ma trận vữa làm mất độ cứng cục bộ, dẫn đến quá trình chùng ứng suất phi tuyến tính rõ rệt.
2. Tranh luận về ảnh hưởng của nhiệt độ: Gutsch và Rostasy (1994) nhận định nhiệt độ vượt quá $20^\circ\text{C}$ sẽ đẩy nhanh tốc độ chùng ứng suất. Tuy nhiên, Morimoto và Koyanagi (1994) lại quan sát thấy ảnh hưởng không đáng kể trong dải nhiệt độ lên đến $60^\circ\text{C}$, trong khi các nghiên cứu của Atrushi (2003) trên mẫu bê tông chứa 5% silica fume ở các mức nhiệt độ $5^\circ\text{C}, 13^\circ\text{C}, 20^\circ\text{C}, 45^\circ\text{C}$ lại chỉ ra sự biến thiên phi hệ thống do tác động chồng lấn của co ngót tự sinh (autogenous shrinkage).

                        TIẾN TRÌNH LÝ THUYẾT & THỰC NGHIỆM CHÙNG ỨNG SUẤT
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ Feldman & Sereda (1968)      Morimoto & Koyanagi (1994)       Gilbert & Ranzi (2010)   │
│ Vi cấu trúc Gel C-S-H   ───► Thí nghiệm chùng kéo/nén   ───►  Mô hình AEMM & phân tích │
│ trượt giữa các lớp nước      trên mẫu vật liệu vữa            ứng xử dài hạn dầm/cột   │
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ CÁC NGHIÊN CỨU KẾT CẤU QUỐC TẾ:                                                        │
│ • Jukka Jokela (1990): Dầm BTCT chịu tải nhiệt kết hợp chuyển vị gối tựa               │
│ • Koichi Maekawa et al. (2006): Dầm liên tục 2 nhịp chịu tải trọng & co ngót          │
│ • Ghali et al. (1969): Chùng ứng suất dầm liên tục nhưng chưa cô lập biến dạng ký sinh │
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │ (Research Gap: Chưa cô lập chuyển vị thuần)
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│ ĐỘT PHÁ CỦA LUẬN ÁN (Nguyễn Mạnh Hùng, 2021):                                          │
│ Thiết kế hệ khung treo khử trọng lượng bản thân + Khung đo biến dạng co ngót đối chứng │
│ ──► Thiết lập công thức giải tích chuẩn xác dự báo RB(t) và hàm quan hệ phi(t)-chi(t)  │
└────────────────────────────────────────────────────────────────────────────────────────┘

Khi định vị với các công trình quốc tế trên cấu kiện thực:

• So sánh với Jukka Jokela (1990): Jokela khảo sát 04 mẫu dầm BTCT chịu đồng thời tải trọng nhiệt và chuyển vị gối tựa, ghi nhận sự suy giảm phản lực gối nhưng biến dạng nhiệt lớn làm lu mờ động học chùng ứng suất do vật liệu.
• So sánh với Koichi Maekawa et al. (2006): Maekawa thực hiện thí nghiệm trên dầm 2 nhịp chịu chuyển vị gối giữa kết hợp tải trọng sử dụng. Mô hình này phản ánh tổng ứng xử làm việc nhưng không thể bóc tách độc lập phần suy giảm lực do chùng ứng suất thuần túy và phần do võng uốn từ biến gây ra bởi ngoại tải.

Luận án của NCS. Nguyễn Mạnh Hùng đã tạo bước tiến quan trọng khi thiết kế giải pháp thực nghiệm độc lập hoàn toàn, triệt tiêu ảnh hưởng của tải trọng bản thân và biến dạng co ngót/nhiệt, lấp đầy khoảng trống dữ liệu thực nghiệm về chùng ứng suất kết cấu tại khu vực Đông Nam Á.

Đó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 nền tảng lý thuyết đàn dẻo nhớt của Bazant (1972) và hướng dẫn CEB-FIP Model Code (MC 90-93: 2010) thông qua việc phát triển 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) cho cấu kiện dầm BTCT chịu uốn có xét đến sự nứt cục bộ:

1. Thiết lập phương trình liên hệ Từ biến - Chùng ứng suất: Dựa trên phương trình tích phân Volterra loại 2, biến dạng dài hạn và ứng suất được liên kết thông qua hàm từ biến $J(t, t_0)$ và hàm chùng ứng suất $R(t, t_0)$:

$$\varepsilon_c(t) = \sigma(t_0) J(t, t_0) + \int_{t_0}^t J(t, \tau) \frac{\partial \sigma(\tau)}{\partial \tau} d\tau$$

$$\sigma(t) = \varepsilon(t_0) R(t, t_0) + \int_{t_0}^t R(t, \tau) \frac{\partial \varepsilon(\tau)}{\partial \tau} d\tau$$

2. Mô hình hóa quan hệ giữa Hệ số từ biến $\varphi(t, t_0)$ và Hệ số già hóa $\chi(t, t_0)$: Luận án đã xây dựng quy luật biến thiên của hệ số già hóa $\chi$ theo thời gian và cấp độ bền bê tông ($f_c = 30\text{ MPa}$). Hệ số $\chi(t, t_0)$ đóng vai trò là trọng số điều chỉnh mức độ suy giảm mô đun đàn hồi hiệu quả $E_{eff}$:

$$E_{eff}(t, t_0) = \frac{E_c(t_0)}{1 + \chi(t, t_0) \varphi(t, t_0)}$$

3. Lý thuyết phân phối lại nội lực có lợi: Khẳng định chuyển vị cưỡng bức gối tựa tạo ra trạng thái ứng suất ban đầu, nhưng nhờ hiện tượng chùng ứng suất, ứng suất nén trong bê tông giảm dần, kéo theo sự giảm độ cong $1/r$ và sự hạ thấp ứng suất kéo trong cốt thép $\sigma_s(t)$, giúp kết cấu tự điều tiết và gia tăng độ bền lâu.

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

Khung phân tích của luận án tích hợp 3 trụ cột lý thuyết: (1) Cơ học kết cấu đàn dẻo nhớt, (2) Lý thuyết tiết diện quy đổi nứt theo TCVN 5574:2018 và Eurocode 2, và (3) Phương pháp xấp xỉ số giải phương trình tích phân:

                  KHUNG PHÂN TÍCH GIẢI TÍCH & THỰC NGHIỆM CỦA LUẬN ÁN
┌────────────────────────────────────────────────────────────────────────────────────────┐
│                                 DỮ LIỆU ĐẦU VÀO VẬT LIỆU                               │
│  • Cấp độ bền bê tông B25-B30 (fc, Ec)         • Cốt thép CB300-V (Es, fy, mu1, mu2, mu3)│
│  • Hàm từ biến J(t, t0) & hệ số phi(t, t0)     • Chuyển vị cưỡng bức cố định Delta0    │
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│                               MÔ PHỎNG GIẢI TÍCH (MATLAB/MAPLE)                         │
│  • Phân vùng dầm: Chiều dài vùng nứt lcr và vùng không nứt (l - lcr)                   │
│  • Độ cứng tiết diện biến thiên: D1 = Ec*Ired (chưa nứt) ──► D2 = Eeff*Ired (đã nứt)  │
│  • Giải tích phân hàm chùng ER(t, t0) qua AEMM có hiệu chỉnh hệ số già hóa chi(t, t0)  │
└──────────────────────────────────────────┬─────────────────────────────────────────────┘
                                           │
                                           ▼
┌────────────────────────────────────────────────────────────────────────────────────────┐
│                                DỰ BÁO SUY GIẢM PHẢN LỰC GỐI                             │
│       RB(t) / RB(t0) = f [ Eeff(t, t0), Ired(t), lcr(t), phi(t, t0), chi(t, t0) ]      │
│  ──► Kiểm chứng với dữ liệu thực nghiệm hệ khung treo độc lập tại LAS-XD125 (Hà Nội)   │
└────────────────────────────────────────────────────────────────────────────────────────┘

Điều kiện biên áp dụng cho mô hình giải tích được giới hạn chặt chẽ:

• Mức ứng suất ban đầu $\sigma_c \le 0{,}5 f_c$ nhằm đảm bảo nằm trong miền từ biến tuyến tính.
• Liên kết: Một đầu ngàm tuyệt đối, một đầu gối tựa di động chịu chuyển vị cưỡng bức $\Delta_0$ thẳng đứng không đổi theo thời gian.
• Tiết diện chữ nhật đặt cốt thép đơn hoặc cốt thép kép, có xét đến tỷ số quy đổi đàn hồi $\alpha_e = E_s / E_c$.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ thế giới quan thực chứng (positivism) với phương pháp thực nghiệm cơ học cấp độ phòng thí nghiệm kết hợp mô hình hóa giải tích. Thiết kế nghiên cứu phân cấp đa mức (multi-level design):

• Cấp độ vật liệu: Xác định các chỉ tiêu cơ lý thực tế bao gồm cường độ nén khối lập phương $R$, cường độ lăng trụ $R_b$, mô đun đàn hồi tức thời $E_c$, đường cong chùng ứng suất của thanh thép $\Phi6$, và thí nghiệm đo hệ số từ biến lăng trụ $\varphi(t, t_0)$ bằng hệ thống khung ép lò xo chuyên dụng.
• Cấp độ cấu kiện: Khảo sát sự suy giảm phản lực $R_B(t)$ trên hệ dầm BTCT thực tế kích thước $b \times h \times L = 100 \times 150 \times 800\text{ mm}$, chiều dài nhịp tính toán $L_0 = 700\text{ mm}$.

                 THIẾT KẾ MÔ HÌNH THỰC NGHIỆM ĐỘC LẬP TẠI LAS-XD125
                 
         Ngàm cố định                                      Gối trượt tạo cưỡng bức
        ┌────────────┐                                          ┌────────────┐
        │ █  █  █  █ │══════════════════════════════════════════│ Load Cell  │
        │ █ (Khung   │      ▲                              ▲    │ & Khung    │
        │ █   ngàm)  │      │ Dầm BTCT: 100x150x800 mm     │    │ Nêm thép   │──► Khóa cứng Delta0
        └────────────┘      │                              │    └────────────┘
              │             └── Hệ treo đối trọng khử G ───┘           │
              │                                                        │
              ▼                                                        ▼
   ┌──────────────────────┐                                 ┌──────────────────────┐
   │ Cặp mẫu dầm D1..D4   │                                 │ Cặp mẫu đối chứng    │
   │ Đo phản lực RB(t)    │                                 │ Đo co ngót eps_sh &  │
   │ & Biến dạng eps_c, s │                                 │ biến dạng nhiệt th   │
   └──────────────────────┘                                 └──────────────────────┘

Ma trận mẫu dầm thí nghiệm gồm 12 cấu kiện chính:

1. Nhóm 1 (D1-10, D2-10, D3-10, D4-10): Đặt $2\Phi10$ vùng kéo ($\mu_1 = 1{,}163%$).
2. Nhóm 2 (D1-12, D2-12, D3-12, D4-12): Đặt $2\Phi12$ vùng kéo ($\mu_2 = 1{,}675%$).
3. Nhóm 3 (D1-14, D2-14, D3-14, D4-14): Đặt $2\Phi14$ vùng kéo ($\mu_3 = 2{,}279%$).

Quy trình nghiên cứu rigorous

Nhằm đạt độ tin cậy và tính hợp lệ cấu trúc (construct validity) cao nhất, tác giả đã chế tạo hệ thống thiết bị thí nghiệm chuyên biệt giải quyết triệt để 3 thách thức lớn:

1. Loại bỏ ảnh hưởng của trọng lượng bản thân: Dầm được treo trên hệ khung thép định hình có cơ cấu đối trọng cân bằng (Hình 3.8 và Hình 3.9), triệt tiêu hoàn toàn mô men uốn do tĩnh tải bản thân gây ra.
2. Khóa cứng chuyển vị cưỡng bức: Sử dụng hệ khung trượt kết hợp kích thủy lực tạo chuyển vị ban đầu $\Delta_0$, sau đó khóa cố định tuyệt đối bằng hệ nêm thép cơ học (Hình 3.7). Chuyển vị được giám sát liên tục bằng đồng hồ so điện tử độ chính xác $0{,}001\text{ mm}$.
3. Bóc tách biến dạng co ngót và nhiệt độ: Thiết lập song song một nhóm tổ mẫu dầm đối chứng không chịu lực cùng cấp phối và điều kiện bảo dưỡng (Hình 3.11, Phụ lục 4) để đo trực tiếp $\varepsilon_{sh}(t)$ và $\varepsilon_{th}(t)$, làm cơ sở khử trừ số liệu.

Hệ thống đo lường sử dụng:

• Đầu đo lực (Load Cell) kỹ thuật số kiểm soát phản lực gối tựa $R_B(t)$.
• Tensometer điện trở (Strain Gauges) dán trên bề mặt bê tông vùng nén ($\varepsilon_c$) và hàn trực tiếp trên thanh cốt thép dọc vùng kéo ($\varepsilon_s$).
• Cầu đo biến dạng tĩnh chuyên dụng thu thập số liệu tự động theo chu kỳ thời gian thực (1 giờ, 3 giờ, 6 giờ, 12 giờ, 24 giờ, 3 ngày, 7 ngày, 28 ngày, 90 ngày, 180 ngày).

Data và phân tích

• Đặc trưng vật liệu: Bê tông sử dụng xi măng PC40, đá dăm $D_{max} = 15\text{ mm}$, cát vàng sông Lô; cốt thép CB300-V có mô đun đàn hồi $E_s = 2{,}0 \times 10^5\text{ MPa}$.
• Kỹ thuật tính toán nâng cao: Luận án thiết lập thuật toán tính lặp trên nền tảng Matlab (Phụ lục 8) để giải hệ phương trình xác định mô đun chùng $E_R(t, t_0)$ và sử dụng phần mềm Maple (Phụ lục 9) để giải phương trình vi phân xác định sự suy giảm phản lực gối tựa $R_B(t)$ theo thời gian.
• Độ tin cậy (Reliability): Kết quả thí nghiệm thu được từ 4 mẫu dầm trong mỗi tổ mẫu cho thấy độ lệch chuẩn cực tiểu ($CV < 5%$), hệ số xác định tương quan giữa mô hình lý thuyết và thực nghiệm đạt $R^2 > 0{,}94$.

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

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

1. Quy luật chùng ứng suất 3 giai đoạn của dầm BTCT: Quá trình suy giảm phản lực gối tựa $R_B(t)$ diễn ra theo 3 giai đoạn động học rõ rệt:
   • Giai đoạn 1 (0 – 72 giờ đầu): Tốc độ chùng ứng suất diễn ra mãnh liệt nhất, giá trị phản lực $R_B(t)$ suy giảm từ 15% đến 22% so với giá trị ban đầu $R_B(t_0)$.
   • Giai đoạn 2 (3 ngày – 90 ngày): Tốc độ suy giảm chậm dần, đạt xấp xỉ 50% tổng lượng chùng ứng suất dài hạn sau 2 đến 3 tháng.
   • Giai đoạn 3 (sau 90 ngày đến 360 ngày): Đường cong phản lực tiệm cận giá trị ổn định, mức suy giảm phản lực tích lũy đạt từ 32% đến 44% tùy thuộc hàm lượng cốt thép.
2. Ảnh hưởng của hàm lượng cốt thép $\mu$: Nhóm dầm có hàm lượng cốt thép thấp ($\mu_1 = 1{,}163%$) có tỷ lệ suy giảm phản lực $R_B(t)/R_B(t_0)$ lớn hơn nhóm dầm có hàm lượng cốt thép cao ($\mu_3 = 2{,}279%$). Cốt thép với mô đun đàn hồi $E_s$ không đổi đóng vai trò là khung chịu lực đàn hồi, kìm hãm biến dạng từ biến của bê tông và hạn chế mức độ chùng ứng suất tổng thể của cấu kiện.
3. Quy luật chuyển dịch ứng suất từ Bê tông sang Cốt thép: Khi phản lực gối tựa $R_B(t)$ giảm, mô men uốn tại tiết diện ngàm giảm theo, dẫn đến biến dạng nén $\varepsilon_c(t)$ của bê tông vùng nén giảm rõ rệt. Đồng thời, biến dạng kéo của cốt thép $\varepsilon_s(t)$ ổn định và không phát triển thêm vết nứt mới sau 28 ngày chịu chuyển vị cưỡng bức.
4. Kiểm chứng hàm già hóa $\chi(t, t_0)$: Thực nghiệm chứng minh giá trị hệ số già hóa của bê tông trong điều kiện miền Bắc Việt Nam dao động trong khoảng $0{,}65 \le \chi(t, t_0) \le 0{,}85$, phù hợp cao với mô hình tính toán AEMM điều chỉnh.

┌────────────────────────────────────────────────────────────────────────────────────────┐
│                ĐỘNG HỌC SUY GIẢM PHẢN LỰC GỐI TỰA RB(t)/RB(t0) THEO THỜI GIAN          │
│                                                                                        │
│  RB(t)/RB(t0)                                                                          │
│   1.00 ──┐                                                                             │
│   0.90   │\  Giai đoạn 1: 0 - 72h (Suy giảm cực nhanh: 15 - 22%)                       │
│   0.80   │ \                                                                           │
│   0.70   │  \─── Giai đoạn 2: 3 - 90 ngày (Tốc độ giảm dần, đạt ~50% lượng chùng)     │
│   0.60   │      \────────────────────────── Giai đoạn 3: > 90 ngày (Tiệm cận ổn định) │
│   0.50   │                                 ─────────────────────────────────           │
│   0.00 ──┴─────┬──────────┬──────────┬──────────┬──────────┬──────────┬──────────►      │
│         t0   3 ngày    28 ngày    90 ngày    180 ngày   270 ngày   360 ngày (Thời gian)│
│                                                                                        │
│  [Nhóm mu1 = 1.163%]: Suy giảm mạnh nhất (~40-44%) do ít cốt thép kiềm chế đàn hồi    │
│  [Nhóm mu3 = 2.279%]: Suy giảm thấp nhất (~32-35%) do hàm lượng thép cao giữ độ cứng   │
└────────────────────────────────────────────────────────────────────────────────────────┘

Implications đa chiều

• Đóng góp cho lý thuyết kết cấu: Bổ sung cơ sở khoa học hoàn chỉnh về bài toán tương tác đàn dẻo nhớt trong dầm BTCT nứt chịu uốn. Chứng minh rằng việc bỏ qua hiện tượng chùng ứng suất khi tính toán dầm chịu lún lệch hoặc co ngắn cột sẽ dẫn đến đánh giá quá cao nội lực dài hạn (overestimation) từ 30% đến 45%.
• Ứng dụng thực tiễn kỹ thuật (Practical Engineering): Cho phép các kỹ sư thiết kế tối ưu hóa tiết diện và bố trí cốt thép cho hệ dầm móng công trình và dầm nối (coupling beam/link beam) trong các tòa nhà siêu cao tầng, tránh việc gia cường cốt thép lãng phí khi tính toán theo sơ đồ đàn hồi thuần túy.
• Khuyến nghị chính sách và tiêu chuẩn: Cung cấp cơ sở thực nghiệm vững chắc để ban soạn thảo tiêu chuẩn quốc gia xem xét bổ sung phụ lục tính toán chùng ứng suất và hệ số già hóa $\chi$ vào các phiên bản cập nhật tiếp theo của TCVN 5574.

Limitations và Future Research

Mặc dù đạt được những bước tiến đột phá, luận án có một số giới hạn biên (boundary conditions):

1. Phạm vi vật liệu: Nghiên cứu chỉ tập trung vào bê tông nặng thông thường cấp B25–B30; chưa khảo sát trên bê tông cường độ cao (HPC, $f_c > 60\text{ MPa}$), bê tông hạt mịn hoặc bê tông cốt sợi phân tán (FRC).
2. Kịch bản chuyển vị: Thí nghiệm áp đặt chuyển vị cưỡng bức tức thời tại thời điểm $t_0$ và giữ không đổi $\Delta(t) = \text{const}$. Trong thực tế địa kỹ thuật, lún lệch của móng hoặc co ngắn cột thường diễn ra tiệm tiến theo tiến độ thi công tải trọng tầng nhà $\Delta(t) = f(t)$.
3. Trạng thái ứng suất: Nghiên cứu tập trung vào cấu kiện dầm chịu uốn đơn thuần một nhịp tĩnh định; chưa mở rộng cho hệ kết cấu không gian 3D phức tạp có xét đến liên kết sàn toàn khối.

Chương trình nghiên cứu tiếp theo (Future Research Agenda):

• Phát triển mô hình giải tích cho chuyển vị cưỡng bức tăng dần theo thời gian $\Delta(t)$ tương thích với tiến độ xây dựng nhà nhiều tầng.
• Mở rộng nghiên cứu thực nghiệm trên dầm bê tông cốt thép dự ứng lực (prestressed concrete) và bê tông geopolymer thân thiện môi trường.
• Ứng dụng thuật toán học máy (Machine Learning) và mạng nơ-ron để nhận dạng tham số hàm chùng ứng suất từ dữ liệu quan trắc kết cấu thực tế (Structural Health Monitoring - SHM).

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

• Tác động học thuật (Academic Impact): Tạo lập bộ dữ liệu thực nghiệm chuẩn (benchmark dataset) gồm 12 mẫu dầm được theo dõi dài hạn với độ tin cậy cao, mở ra hướng nghiên cứu mới về cơ học viscoelasticity cho các nghiên cứu sinh và học giả trong nước và quốc tế.
• Chuyển đổi ngành công nghiệp xây dựng: Giúp các tổng thầu, đơn vị tư vấn thiết kế kết cấu (như VNCC, CDC, Coninco) giảm chi phí vật liệu cốt thép từ 10% đến 15% trong các cấu kiện chịu lún lệch hoặc co ngắn cưỡng bức, đồng thời nâng cao tuổi thọ công trình thông qua kiểm soát nứt chính xác.
• Lợi ích xã hội và bền vững: Giảm lượng tiêu thụ xi măng và thép xây dựng, trực tiếp cắt giảm phát thải khí nhà kính $\text{CO}_2$, đóng góp vào mục tiêu phát triển công trình xanh và bền vững.

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

                               CÁC NHÓM ĐỐI TƯỢNG HƯỞNG LỢI
┌──────────────────────────────────────┬──────────────────────────────────────┐
│       NGHIÊN CỨU SINH & HỌC GIẢ      │        KỸ SƯ THIẾT KẾ KẾT CẤU        │
│ • Tiếp cận phương pháp luận cô lập   │ • Sở hữu công thức tính toán AEMM    │
│   biến dạng mẫu treo độc lập.        │   dự báo chính xác phản lực RB(t).   │
│ • Sử dụng bộ số liệu thực nghiệm làm │ • Tối ưu hóa cốt thép dầm móng và    │
│   benchmark kiểm chứng mô hình FEM.  │   dầm giằng cột-vách nhà cao tầng.   │
├──────────────────────────────────────┼──────────────────────────────────────┤
│      CHỦ ĐẦU TƯ & DOANH NGHIỆP       │        BAN SOẠN THẢO TIÊU CHUẨN      │
│ • Tiết kiệm 10-15% khối lượng thép   │ • Dữ liệu khoa học thực tế tại VN để │
│   ở các cấu kiện dầm chuyển vị.      │   bổ sung hướng dẫn tính chùng       │
│ • Hạn chế chi phí xử lý nứt nẻ móng. │   ứng suất vào TCVN 5574.            │
└──────────────────────────────────────┴──────────────────────────────────────┘

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ì?

Luận án đã mở rộng thành công lý thuyết Mô đun đàn hồi hiệu quả điều chỉnh theo tuổi (AEMM) của Bazant (1972) và CEB-FIP cho bài toán dầm BTCT nứt chịu chuyển vị cưỡng bức gối tựa. Đóng góp độc bản là việc thiết lập giải pháp đóng (closed-form solution) xác định hàm chùng ứng suất tương đương $E_R(t, t_0)$ gắn liền với quy luật già hóa $\chi(t, t_0)$ của bê tông trong điều kiện nhiệt đới ẩm Việt Nam.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây?

So với Jukka Jokela (1990) và Koichi Maekawa et al. (2006), luận án đã thiết kế thành công hệ khung treo đối trọng độc lập khử trọng lượng bản thân kết hợp hệ mẫu dầm đối chứng song song khử co ngót và biến dạng nhiệt. Giải pháp này cho phép đo đạc và lượng hóa chính xác 100% động học chùng ứng suất do chuyển vị cưỡng bức thuần túy mà không bị nhiễu bởi ngoại tải hay biến dạng môi trường.

3. Phát hiện thực nghiệm nào gây bất ngờ nhất?

Phát hiện về tốc độ suy giảm phản lực gối tựa cực lớn trong 72 giờ đầu tiên (đạt tới 15–22% tổng phản lực ban đầu) và sự suy giảm này đạt tới hơn 50% tổng lượng chùng chỉ sau 2 đến 3 tháng. Điều này chứng minh quá trình phân phối lại ứng suất diễn ra ngay tức thì sau khi xuất hiện chuyển vị cưỡng bức, đóng vai trò giải tỏa nguy cơ nứt vỡ sớm cho bê tông.

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

Có. Quy trình thí nghiệm được mô tả chi tiết từ bản vẽ chế tạo khung nêm thép khóa chuyển vị, hệ thống treo cân bằng, sơ đồ bố trí tensometer, lò xo giữ tải đo từ biến, cho đến toàn bộ mã nguồn tính lặp trên Matlab (Phụ lục 8) và bảng tính tích phân giải tích trên Maple (Phụ lục 9).

5. Định hướng nghiên cứu 10 năm tới được phác thảo thế nào?

Mở rộng nghiên cứu bài toán chùng ứng suất sang hệ kết cấu không gian 3 chiều chịu chuyển vị biến thiên theo thời gian $\Delta(t)$, tích hợp ảnh hưởng của cốt sợi phi kim (FRP) và xây dựng hệ thống giám sát sức khỏe công trình (SHM) dựa trên nền tảng kỹ thuật số (Digital Twin).

Kết luận

1. Thiết lập thành công mô hình thực nghiệm chuẩn xác: Xây dựng hệ thống thí nghiệm dầm BTCT chịu chuyển vị cưỡng bức gối tựa có khả năng khử hoàn toàn ảnh hưởng của trọng lượng bản thân, co ngót và nhiệt độ, đảm bảo tính khách quan tuyệt đối của dữ liệu chùng ứng suất.
2. Khám phá quy luật động học 3 giai đoạn: Chứng minh phản lực gối tựa suy giảm từ 32% đến 44% sau 1 năm, trong đó hơn 50% lượng suy giảm tập trung trong 2-3 tháng đầu tiên.
3. Lượng hóa ảnh hưởng của hàm lượng cốt thép: Xác định rõ cơ chế cốt thép kìm hãm chùng ứng suất; hàm lượng cốt thép tăng từ $\mu_1 = 1{,}163%$ lên $\mu_3 = 2{,}279%$ làm giảm mức độ suy giảm phản lực tương đối của dầm.
4. Hoàn thiện công thức giải tích AEMM cải tiến: Xây dựng phương pháp tính toán dự báo chùng ứng suất và hệ số già hóa $\chi(t, t_0)$ trên phần mềm Maple và Matlab, cho kết quả tương đồng cao với thực nghiệm ($R^2 > 0{,}94$).
5. Mở ra 3 nhánh nghiên cứu mới: (1) Động lực học chùng ứng suất dưới chuyển vị biến thiên $\Delta(t)$, (2) Chùng ứng suất trong bê tông cường độ cao và bê tông tính năng siêu cao (UHPC), (3) Tối ưu hóa thiết kế dầm giằng nhà siêu cao tầng theo trạng thái giới hạn sử dụng dài hạn.
6. Giá trị thực tiễn và chuyển giao: Cung cấp cơ sở khoa học định lượng giúp các kỹ sư kết cấu tự tin tối ưu hóa thiết kế dầm móng và dầm liên kết cột - vách, mang lại hiệu quả kinh tế kỹ thuật cao và đóng góp trực tiếp vào sự phát triển của hệ thống tiêu chuẩn xây dựng Việt Nam.