Tổng quan về luận án
Sự suy thoái kết cấu bê tông cốt thép (BTCT) do hiện tượng ăn mòn cốt thép trong môi trường xâm thực mạnh—đặc biệt là vùng ven biển, hải đảo, và các công trình giao thông chịu ion clorua—đặt ra bài toán cấp thiết về mặt kỹ thuật và kinh tế công trình. Thanh polyme cốt sợi thủy tinh (Glass Fiber-Reinforced Polymer - GFRP) với các đặc tính vượt trội như không bị ăn mòn điện hóa, cường độ chịu kéo rất cao (từ 483 MPa đến 1600 MPa) và khối lượng riêng thấp (1,25 đến 2,10 g/cm³, chỉ bằng 1/6 đến 1/4 so với thép 7,90 g/cm³) nổi lên như một giải pháp thay thế tiềm năng. Tuy nhiên, rào cản kỹ thuật lớn nhất của cốt GFRP thuần túy là tính đàn hồi - giòn tuyến tính cho đến khi đứt (không có thềm chảy dẻo) và mô đun đàn hồi thấp ($E_f = 35 \div 51\text{ GPa}$, xấp xỉ 1/4 mô đun cốt thép $E_s = 200\text{ GPa}$). Điều này khiến dầm bê tông cốt GFRP thuần túy có độ võng rất lớn, bề rộng vết nứt mở rộng nhanh và xảy ra phá hoại giòn đột ngột, làm hạn chế nghiêm trọng khả năng đáp ứng trạng thái giới hạn sử dụng (SLS).
Nhằm giải quyết triệt để nghịch lý cơ học trên, giải pháp kết hợp cốt thép và cốt GFRP trong cùng một tiết diện dầm—hình thành hệ cốt hỗn hợp SGFRP (Steel and Glass Fiber Reinforced Polymer)—được đề xuất. Trong hệ kết cấu lai này, cốt thép với mô đun đàn hồi cao ($E_s = 200\text{ GPa}$) được bố trí ở lớp bên trong với chiều dày lớp bê tông bảo vệ lớn nhằm chống ăn mòn, trong khi cốt GFRP chịu lực chính đặt ở lớp ngoài cùng tiếp xúc với môi trường xâm thực.
Khoảng trống nghiên cứu (Research Gap) then chốt hiện nay nằm ở chỗ: hầu hết các nghiên cứu quốc tế (như các tiêu chuẩn ACI 440.1R của Hoa Kỳ, CSA-S806 của Canada, hay FIB Model Code của Châu Âu) chỉ tập trung vào cốt FRP đơn lẻ hoặc tiếp cận dầm hỗn hợp bằng cách quy đổi thô sơ về một hàm lượng cốt tương đương (tương đương thép hoặc tương đương FRP). Cách quy đổi này bỏ sót các dạng phá hoại phức tạp, đặc biệt là dạng phá hoại đứt cốt GFRP sau khi cốt thép đã chảy dẻo nhưng bê tông vùng nén chưa bị nghiền nát. Hơn nữa, toàn bộ hệ thống lý thuyết quốc tế đều xây dựng trên triết lý thiết kế và hệ số an toàn riêng của Bắc Mỹ/Châu Âu, chưa tương thích với hệ thống tiêu chuẩn thiết kế kết cấu bê tông hiện hành của Việt Nam (TCVN 5574:2018) vốn kế thừa phương pháp trạng thái giới hạn của Nga (SNiP).
Luận án thiết lập 4 câu hỏi nghiên cứu và hệ giả thuyết định lượng tương ứng:
- Câu hỏi 1 (RQ1): Cơ chế ứng xử chịu uốn thực tế của dầm bê tông cốt SGFRP diễn ra qua những trạng thái giới hạn nào và làm thế nào để mô hình hóa giải tích quan hệ ứng suất - biến dạng phi tuyến?
- Giả thuyết 1 (H1): Quan hệ ứng suất - biến dạng của hệ cốt hỗn hợp SGFRP có thể mô hình hóa chính xác qua đường song tuyến hai giai đoạn (bilinear), trong đó giai đoạn I đặc trưng cho sự làm việc đồng thời đàn hồi của thép và GFRP, giai đoạn II đặc trưng cho sự tăng cường ứng suất thuần túy của GFRP sau khi thép chảy dẻo.
- Câu hỏi 2 (RQ2): Quy trình giải thuật số lặp từng bước có thể dự đoán chính xác khả năng chịu mô men uốn ($M_u$), quan hệ mô men - độ cong ($M-\phi$) và độ võng ($f$) của dầm theo TCVN 5574:2018 hay không?
- Giả thuyết 2 (H2): Giải thuật số tích phân phân mảnh vùng nén bê tông kết hợp điều kiện tương thích biến dạng phẳng Bernoulli cho phép kiểm soát sai số khả năng chịu uốn dưới 10% so với thực nghiệm.
- Câu hỏi 3 (RQ3): Tỷ lệ diện tích cốt dọc $\omega = A_f / A_s$ (hoặc $m = A_s / A_f$) ảnh hưởng thế nào đến sự chuyển tiếp giữa phá hoại dẻo và phá hoại giòn, cũng như độ võng và bề rộng vết nứt ($a_{crc}$)?
- Giả thuyết 3 (H3): Tồn tại một miền giá trị biên tối ưu của tỷ lệ $\omega = A_f / A_s$ đảm bảo đồng thời tiêu chí cường độ cực hạn (ULS), khả năng kiểm soát nứt - võng (SLS) và chỉ số độ dẻo kết cấu ($DF$).
- Câu hỏi 4 (RQ4): Tính chuẩn xác của các công thức giải tích đề xuất khi đối chứng với thực nghiệm chế tạo dầm tại Việt Nam và cơ sở dữ liệu quốc tế?
- Giả thuyết 4 (H4): Hệ công thức giải tích đóng mở rộng từ TCVN 5574:2018 có độ hội tụ cao và an toàn trong thiết kế thực hành dầm nhịp đơn chịu uốn thuần túy.
Phạm vi nghiên cứu (Scope) tập trung vào dầm bê tông cốt SGFRP tiết diện chữ nhật đặt cốt đơn, chịu tải trọng tĩnh ngắn hạn tác dụng tại hai điểm uốn thuần túy, sử dụng bê tông có cấp độ bền đến B70 (mác 400). Chương trình thực nghiệm được triển khai đồng bộ tại Phòng thí nghiệm LAS-XD125, Trường Đại học Xây dựng Hà Nội, cung cấp bộ dữ liệu đo kiểm đầy đủ từ biến dạng vi mô (strain gauges) đến chuyển vị vĩ mô (LVDTs).
Literature Review và Positioning
Nghiên cứu kết cấu bê tông cốt hỗn hợp thép - FRP trên thế giới đã trải qua hơn hai thập kỷ phát triển với ba trường phái học thuật chính:
Trường phái thực nghiệm tiên phong khởi đầu từ gợi ý của Tan (1997), tiếp nối bởi nghiên cứu kinh điển của Maria Antonietta Aiello và Luciano Ombres (2002) khi thí nghiệm 6 dầm nhịp 2,7 m sử dụng cốt AFRP kết hợp cốt thép ($f_y = 465\text{ MPa}$, cường độ kéo AFRP đạt 1647 MPa, $\omega = 0,88 \div 1,56$). Nhóm tác giả khẳng định việc bổ sung cốt thép làm giảm mạnh độ võng và phân tán vết nứt, song chưa định lượng được ranh giới các dạng phá hoại. Wenjun Qu và cộng sự (2009) thử nghiệm 8 mẫu dầm nhịp 1,8 m cốt thép - GFRP ($E_f = 41\text{ GPa}$), tiếp cận bằng phương pháp quy đổi diện tích cốt hỗn hợp về cốt thép tương đương, nhưng bộc lộ nhược điểm là không nhận diện được pha phá hoại đứt cốt GFRP sau khi thép chảy. Denvid Lau và Hoat Joen Pam (2010) với 12 mẫu dầm nhịp 4,2 m lại quy đổi về hàm lượng cốt FRP tương đương, chứng minh độ dẻo tăng khi bố trí cốt hỗn hợp hàm lượng cao (over-reinforced), song lại bỏ qua dạng phá hoại bê tông vùng nén bị vỡ sớm khi thép chưa kịp chảy. Young-Soo Yoon và cộng sự (2011) khảo sát dầm bê tông cường độ cao 90 MPa với các tổ hợp cốt thép, GFRP và CFRP, đề xuất mô men quán tính hiệu quả $I_e$ nhưng còn hạn chế do mẫu bị phá hoại cắt sớm. Gần đây hơn, WenJie Ge và cộng sự (2011) khảo sát tỷ lệ $\omega = 0,43 \div 0,96$ trên cốt BFRP; Liu Yinghao và Yuan Yong (2012) thử nghiệm bê tông 80 MPa với cốt GFRP cường độ 1301 MPa; Ahmed El Refai và cộng sự (2015) trên 9 mẫu dầm nhịp 3,7 m kết luận dầm hỗn hợp cần hàm lượng cốt dọc vượt 1,4 lần hàm lượng cân bằng; và Xiangjie Ruan và cộng sự (2020) nghiên cứu hiệu ứng bố trí cốt thép và GFRP 1 lớp so với 2 lớp với các tỷ lệ $\omega = 1,0 \div 2,0$.
Trường phái mô phỏng số phần tử hữu hạn (FEA) ghi nhận các công trình đột phá của R. Hawileh (2015) sử dụng ANSYS (phần tử Solid 65, Link 8, mô hình dính bám COMBIN14); Francesco Bencardino và cộng sự (2016) cùng Renyuan Qin, Ao Zhou, Denvid Lau (2017) phát triển mô hình 2D/3D trên ABAQUS (phần tử C3D8R, T3D2) kiểm chứng dữ liệu thực nghiệm với sai số dưới 10%; và Linh Van Hong Bui, Boochai Stitmannaithum, Tamon Ueda (2018) mô phỏng tiêu tán năng lượng, chỉ ra dầm hỗn hợp GFRP-thép đạt 80% năng lượng hấp thụ của dầm BTCT truyền thống khi $\omega < 0,641$.
Trường phái lý thuyết giải tích bao gồm các mô hình của Hau Yan Leung (2004) dựa trên biểu đồ ứng suất chữ nhật ACI; B. Jia và cộng sự (2014) thiết lập điều kiện hàm lượng cốt cân bằng trên 18 mẫu dầm; WenJie Ge và cộng sự (2015) xác lập công thức bề rộng khe nứt cho cốt BFRP; Lei Pang, Wenjun Qu và cộng sự (2015) khảo sát mối quan hệ giữa chỉ số độ dẻo và tỷ lệ diện tích $\omega$; Ali Kheyroddin và Fahimeh Maleki (2017) ứng dụng thuật toán di truyền (Genetic Algorithms - GA) tối ưu hóa số mũ mô men quán tính Branson trên 120 tập dữ liệu; Saruhan Kartal và Ilker Kalkan (2019) phân tích mô men nứt $M_{cr}$ theo Eurocode 2 và ACI; và Yang Yang và cộng sự (2020) xây dựng miền giới hạn tỷ số độ cứng kéo dọc trục $\alpha_{\min} \div \alpha_{\max}$.
Cuộc tranh luận học thuật trung tâm nằm ở hai quan điểm đối nghịch: (1) Nhóm nghiên cứu quy đổi đơn nhất (Qu et al. 2009; Bencardino et al. 2016) giả định biến dạng hai vật liệu cốt là đồng nhất tuyệt đối và xem hệ thống như một vật liệu cốt thép mở rộng; (2) Nhóm nghiên cứu cấu trúc phân tầng (Qin & Lau 2017; El Refai 2015) xem xét sự tách biệt về độ cứng đàn hồi giữa các lớp cốt nhưng lại không đưa ra được mô hình quan hệ ứng suất - biến dạng tương đương tổng quát.
Định vị nghiên cứu của luận án: Khắc phục triệt để hạn chế của phương pháp quy đổi đơn nhất, công trình của Phan Minh Tuấn (2021) tiên phong thiết lập mô hình cơ học cốt hỗn hợp tương đương tích hợp đầy đủ hai giai đoạn đàn hồi - chảy dẻo - tăng bền, đồng thời chuyển hóa toàn bộ hệ phương trình cân bằng nội lực sang hệ quy chuẩn thiết kế theo biến dạng phi tuyến của TCVN 5574:2018.
So sánh với hai nghiên cứu quốc tế tiêu biểu:
- So với Ahmed El Refai và cộng sự (2015): Nghiên cứu của El Refai chỉ quy đổi về hàm lượng cốt FRP tương đương và bố trí cốt dọc một lớp (dễ gây ăn mòn cốt thép), trong khi luận án thiết lập sơ đồ phân lớp tối ưu (thép lớp trong, GFRP lớp ngoài) và giải quyết trọn vẹn bài toán phá hoại nén vỡ bê tông trước hoặc sau khi thép chảy dẻo.
- So với Renyuan Qin, Ao Zhou và Denvid Lau (2017): Nghiên cứu của Qin và Lau dừng lại ở các khuyến nghị định tính trên mô hình số ABAQUS ($\omega = 1,0 \div 2,5$), trong khi luận án đưa ra hệ công thức giải tích đóng tường minh cho các giới hạn $\mu_{\min}, \mu_{\max}, \mu_{\max,e}$, liên kết trực tiếp với độ dẻo $DF$ và kiểm soát nứt $a_{crc}$ theo tiêu chuẩn Việt 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 thuyết dầm Euler-Bernoulli và lý thuyết trạng thái giới hạn phi tuyến của bê tông bằng cách tích hợp mô hình làm việc phi tuyến đồng thời của hai loại vật liệu cốt có đặc tính cơ học trái ngược: thép biến dạng dẻo hoàn hảo và GFRP đàn hồi tuyến tính tuyệt đối.
Khung lý thuyết thiết lập mối quan hệ ứng suất - biến dạng tương đương cho cốt hỗn hợp SGFRP ($\sigma_{sf} - \varepsilon_{sf}$) thông qua phương trình:
$$\sigma_{sf} = \frac{\sigma_s A_s + \sigma_f A_f}{A_s + A_f} = \frac{\sigma_s + \omega \sigma_f}{1 + \omega}$$
Trong đó $\omega = A_f / A_s$. Quan hệ ứng suất - biến dạng này được chia thành hai nhánh tường minh:
- Giai đoạn I ($\varepsilon_{sf} \le \varepsilon_{s,el} = R_s / E_s$): Cả thép và GFRP cùng làm việc đàn hồi. Mô đun đàn hồi tương đương giai đoạn I:
$$E_I = \frac{E_s + \omega E_f}{1 + \omega}$$
Cường độ chảy quy đổi của cốt hỗn hợp:
$$R_{sf,y} = \frac{R_s + \omega \varepsilon_{s,el} E_f}{1 + \omega}$$
- Giai đoạn II ($\varepsilon_{s,el} < \varepsilon_{sf} \le \varepsilon_{sfu} = \varepsilon_{fu}$): Cốt thép đã chảy dẻo ($\sigma_s = R_s$), toàn bộ gia số ứng suất do cốt GFRP gánh nhận. Mô đun đàn hồi quy đổi giai đoạn II:
$$E_{II} = \frac{\omega E_f}{1 + \omega}$$
Cường độ chịu kéo giới hạn quy đổi cực hạn:
$$R_{sf,u} = \frac{R_s + \omega R_f}{1 + \omega}$$
Hệ thống định đề và mệnh đề lý thuyết được xây dựng:
- Mệnh đề 1: Tồn tại ba hình thái phá hoại uốn của dầm SGFRP tương ứng với chiều cao vùng nén quy đổi $\xi = x / h_0$:
- Phá hoại giòn nén ($\xi > \xi_{R,y}$): Bê tông vùng nén bị nghiền nát trước khi cốt thép chảy dẻo và GFRP đứt.
- Phá hoại dẻo ($\xi_{R,e} \le \xi \le \xi_{R,y}$): Cốt thép chảy dẻo trước, dầm thể hiện tính dẻo rõ rệt, biến dạng tăng nhanh trước khi bê tông vùng nén đạt biến dạng cực hạn $\varepsilon_{b2}$.
- Phá hoại giòn kéo ($\xi < \xi_{R,e}$): Cốt thép chảy dẻo và tiếp sau đó cốt GFRP bị đứt đột ngột trước khi bê tông vùng nén bị phá hoại do nén.
Khung phân tích độc đáo
Khung phân tích của luận án tích hợp đồng bộ ba lý thuyết nền tảng: (1) Lý thuyết cơ học biến dạng phi tuyến của bê tông theo TCVN 5574:2018; (2) Lý thuyết đàn dẻo của vật liệu tổ hợp composite-kim loại; (3) Lý thuyết trạng thái giới hạn nứt và võng đàn-nhớt.
Luận án thiết lập các giá trị chiều cao tương đối giới hạn của vùng bê tông chịu nén:
- Giới hạn dưới đảm bảo cốt thép đạt giới hạn chảy ($\xi_{R,y}$):
$$\xi_{R,y} = \frac{\varepsilon_{b2}}{\varepsilon_{b2} + \frac{R_s}{E_s}}$$
- Giới hạn trên hiệu quả đảm bảo cốt GFRP không bị kéo đứt đột ngột ($\xi_{R,e}$):
$$\xi_{R,e} = \frac{\varepsilon_{b2}}{\varepsilon_{b2} + \varepsilon_{fu}}$$
Từ các điều kiện biên này, luận án thiết lập biểu đồ hàm lượng cốt thép tối thiểu $\mu_{\min}$, tối đa $\mu_{\max}$, và tối đa hiệu quả $\mu_{\max,e}$, định hình không gian tham số giải tích hoàn chỉnh cho việc thiết kế cấu kiện chịu uốn.
Điều kiện biên áp dụng: Cấu kiện dầm tiết diện chữ nhật đặt cốt đơn, bê tông cấp độ bền $\le \text{B70}$, làm việc dưới tác dụng của tải trọng tĩnh ngắn hạn trong giai đoạn khai thác tiêu chuẩn.
Phương pháp nghiên cứu tiên tiến
Thiết kế nghiên cứu
Luận án áp dụng thế giới quan thực chứng (Positivism) kết hợp chặt chẽ giữa mô hình hóa giải tích toán học, giải thuật số lặp phi tuyến và kiểm chứng thực nghiệm đa tầng.
Thiết kế nghiên cứu thực nghiệm đa mức độ (multi-level design) được xây dựng với hệ thống mẫu dầm tỷ lệ thực $b \times h = 150 \times 250\text{ mm}$, chiều dài toàn phần $L = 2200\text{ mm}$, nhịp tính toán $L_0 = 1900\text{ mm}$. Sơ đồ gia tải 4 điểm tạo vùng uốn thuần túy dài 500 mm ở giữa nhịp. Chương trình thí nghiệm phân chia thành 5 nhóm dầm đại diện (từ D1 đến D5, mỗi nhóm gồm 2 dầm lặp D-1 và D-2 để kiểm soát tính lặp lại):
- Nhóm D1: Dầm cốt hỗn hợp SGFRP bố trí 2 lớp (lớp trong 2Ф10 thép gai, lớp ngoài 2Ф10 GFRP).
- Nhóm D2: Dầm cốt hỗn hợp SGFRP bố trí 2 lớp (lớp trong 2Ф12 thép gai, lớp ngoài 2Ф10 GFRP).
- Nhóm D3: Dầm cốt hỗn hợp SGFRP bố trí 2 lớp (lớp trong 2Ф14 thép gai, lớp ngoài 2Ф10 GFRP).
- Nhóm D4: Dầm cốt hỗn hợp SGFRP bố trí 1 lớp hỗn hợp cùng cốt thép và GFRP.
- Nhóm D5: Dầm đối chứng cốt GFRP thuần túy (2Ф10 GFRP).
Quy trình nghiên cứu rigorous
Quy trình thu thập dữ liệu tuân thủ các tiêu chuẩn thí nghiệm khắt khe:
- Đặc trưng bê tông: Bê tông mác 400 được kiểm định cường độ chịu nén trên mẫu lập phương $150 \times 150 \times 150\text{ mm}$ ($R_{bm} = 40,2\text{ MPa}$), mẫu lăng trụ $150 \times 150 \times 600\text{ mm}$ xác định mô đun đàn hồi $E_b = 32,5\text{ GPa}$, và mẫu trụ lăng trụ xác định $f'_c$.
- Đặc trưng cốt thép: Thí nghiệm kéo mẫu thép xác định giới hạn chảy thực tế $f_y = 390 \div 420\text{ MPa}$, mô đun $E_s = 2,05 \times 10^5\text{ MPa}$, độ giãn dài cực hạn $\delta > 18%$.
- Đặc trưng cốt GFRP: Thí nghiệm kéo theo tiêu chuẩn ASTM D7205 xác định cường độ kéo đứt $R_{fu} = 980\text{ MPa}$, mô đun đàn hồi $E_f = 46,5\text{ GPa}$, không có biến dạng chảy.
- Hệ thống quan trắc: Sử dụng cảm biến đo biến dạng điện trở (strain gauges dán trực tiếp lên thanh thép, thanh GFRP và thớ nén bê tông trên cùng), hệ thống LVDTs đo độ võng tại gối và giữa nhịp, kính lúp quang học chuyên dụng có độ chính xác 0,01 mm đo bề rộng và theo dõi sự phân bố của khe nứt.
Data và phân tích
Song song với thực nghiệm, luận án phát triển thuật toán số lặp từng bước (Iterative Numerical Algorithm) viết trên nền tảng tính toán khoa học:
- Tiết diện bê tông vùng nén được chia thành $n$ dải phần tử hữu hạn ($h_i = h/n$).
- Với mỗi bước tăng biến dạng thớ nén ngoài cùng $\varepsilon_{b}$, thiết lập quá trình lặp đúng dần để xác định vị trí trục trung hòa $x$ thỏa mãn điều kiện cân bằng lực dọc:
$$\Sigma N = C_b(\varepsilon_b, x) - T_s(\varepsilon_s) - T_f(\varepsilon_f) = 0$$
- Trong đó, ứng suất bê tông được tích phân phi tuyến theo biểu đồ biến dạng của TCVN 5574:2018:
$$\sigma_{bi} = R_b \left[ 1 - \left(1 - \frac{\varepsilon_{bi}}{\varepsilon_{b0}}\right)^2 \right] \quad \text{khi } \varepsilon_{bi} \le \varepsilon_{b0}$$
- Sau khi cân bằng lực hội tụ (dung sai $\le 10^{-4}\text{ kN}$), mô men uốn kháng uốn $M_j$, độ cong $\phi_j = \varepsilon_b / x$, và độ võng $f = S \cdot \phi_j \cdot L^2$ được truy xuất tự động.
Thuật toán số được kiểm chứng chéo (cross-validation) trên cơ sở dữ liệu của 30 dầm thực nghiệm quốc tế đã công bố từ các nhóm tác giả Aiello & Ombres (2002), Qu et al. (2009), Lau & Pam (2010), Yoon et al. (2011), Ge et al. (2011), Hawileh (2015), Ruan et al. (2020), cho hệ số tương quan $R^2 = 0,968$ và độ lệch chuẩn tương đối dưới 6,5%.
Phát hiện đột phá và implications
Những phát hiện then chốt
- Khắc phục hoàn toàn dạng phá hoại giòn của dầm GFRP: Thực nghiệm chứng minh dầm đối chứng GFRP thuần túy (D5) bị phá hoại giòn đột ngột với độ võng tức thời tại thời điểm gãy lớn ($f_u = 38,2\text{ mm}$), vết nứt mở rộng nhanh ($a_{crc} > 1,2\text{ mm}$). Ngược lại, các dầm cốt hỗn hợp SGFRP (D1, D2, D3) đều thể hiện pha chảy dẻo rõ nét; sau khi cốt thép chảy dẻo, dầm tiếp tục chịu tải tăng dần nhờ sự tăng cường ứng suất tuyến tính của cốt GFRP cho đến khi bê tông vùng nén đạt $\varepsilon_{b2} = 0,0035$.
- Cải thiện vượt bậc khả năng kiểm soát võng và nứt (SLS): Dưới cùng cấp tải trọng tiêu chuẩn khai thác ($P = 40\text{ kN}$), độ võng của dầm hỗn hợp SGFRP nhóm D1 giảm 42,5% và nhóm D3 giảm 58,1% so với dầm GFRP thuần túy D5. Bề rộng khe nứt $a_{crc}$ ở dầm SGFRP giảm từ 2,1 đến 2,8 lần so với dầm D5, khoảng cách giữa các khe nứt ngắn lại ($L_{sf}$ giảm 35%), thể hiện sự dính bám vượt trội và khả năng khống chế nứt của mạng cốt hỗn hợp.
- Xác lập khoảng tỷ lệ diện tích cốt dọc tối ưu $\omega = A_f / A_s$: Dữ liệu thực nghiệm và khảo sát tham số số học khẳng định:
- Khi $\omega < 0,5$: Dầm hoạt động gần như dầm BTCT truyền thống, không tận dụng được thế mạnh cường độ cao của GFRP.
- Khi $\omega > 2,0$: Tính dẻo của dầm suy giảm nghiêm trọng ($DF < 2,0$), ứng xử uốn tiến dần về dầm GFRP giòn thuần túy.
- Khoảng tỷ lệ tối ưu được luận án đề xuất là $\omega = A_f / A_s \in [0,6; 1,5]$ (tương ứng $m = A_s / A_f \in [0,67; 1,67]$), đảm bảo chỉ số độ dẻo $DF \ge 2,5$, độ võng và vết nứt thỏa mãn TCVN 5574:2018, đồng thời phát huy tối đa 75–85% cường độ kéo của thanh GFRP.
- Độ chính xác vượt trội của hệ công thức giải tích đề xuất: So sánh giữa mô men uốn giới hạn tính toán $M_{u,LT}$ và thực nghiệm $M_{u,TN}$ trên toàn bộ các dầm D1–D4 cho sai số trung bình chỉ $4,2 \div 7,8%$, mô men nứt $M_{cr}$ sai số dưới 8,5%, chứng minh độ tin cậy tuyệt đối của phương pháp cốt hỗn hợp tương đương.
- Ưu thế kết cấu của phương án bố trí cốt phân lớp (2 lớp): Dầm bố trí 2 lớp (thép lớp trong, GFRP lớp ngoài - nhóm D1) có khả năng kiểm soát vết nứt tốt hơn 18,4% và giữ mô men uốn giới hạn cao hơn 6,2% so với dầm bố trí hỗn hợp trên cùng 1 lớp (nhóm D4) có cùng hàm lượng cốt.
Implications đa chiều
- Đóng góp lý thuyết: Bổ sung vào lý thuyết kết cấu bê tông phi kim mô hình cơ học song tuyến tính tương đương SGFRP, xác lập chuẩn mực mới trong việc phân loại ba trạng thái phá hoại uốn.
- Đổi mới phương pháp luận: Xây dựng quy trình tích hợp giữa giải thuật số lặp phi tuyến và công thức giải tích đóng, có thể chuyển giao trực tiếp vào việc lập trình các phần mềm thiết kế kết cấu chuyên dụng.
- Ứng dụng thực tiễn: Cung cấp cẩm nang tính toán chi tiết cho kỹ sư kết cấu trong việc thiết kế các công trình cảng biển, cầu vượt biển, kè chắn sóng, đê kè ven biển miền Trung và các công trình ngầm xâm thực mặn tại Việt Nam.
- Khuyến nghị chính sách và quy chuẩn: Tạo tiền đề khoa học vững chắc để Bộ Xây dựng và Viện Khoa học Công nghệ Xây dựng (IBST) biên soạn phụ lục tiêu chuẩn tính toán kết cấu bê tông cốt sợi polyme kết hợp cốt thép trong các lần soát xét TCVN 5574 tiếp theo.
Limitations và Future Research
Mặc dù đạt được những kết quả đột phá, luận án thẳng thắn ghi nhận 4 giới hạn nghiên cứu:
- Cấp độ bền bê tông: Nghiên cứu thực nghiệm mới chỉ thực hiện trên bê tông mác 400 (cấp độ bền B30/B35); các cấp độ bền siêu cao (UHPC, B80–B120) chưa được kiểm chứng trực tiếp.
- Chế độ tải trọng: Mới dừng lại ở tải trọng tĩnh ngắn hạn tác dụng đơn kỳ; chưa khảo sát ứng xử mỏi (fatigue), tải trọng động đất chu kỳ (cyclic loading) và biến dạng từ biến/co ngót dài hạn (creep and shrinkage) của cốt SGFRP.
- Hình học cấu kiện: Giới hạn ở dầm tiết diện chữ nhật nhịp đơn đặt cốt đơn; chưa mở rộng cho dầm chữ T, dầm liên tục nhiều nhịp chịu mô men âm tại gối tựa.
- Cơ chế dính bám - trượt: Giả thiết độ dính bám hoàn hảo giữa bê tông và cốt GFRP/thép được áp dụng trong mô hình giải tích; hiện tượng trượt cục bộ tại khe nứt lớn chưa được mô hình hóa riêng biệt.
Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):
- Mở rộng nghiên cứu sự làm việc chịu uốn của dầm liên tục nhiều nhịp cốt hỗn hợp SGFRP xét đến hiện tượng phân phối lại nội lực (moment redistribution).
- Khảo sát độ bền lâu dài, suy giảm cường độ dính bám dưới tác động đồng thời của tải trọng mỏi và môi trường sương muối biển gia tốc.
- Nghiên cứu ứng xử chịu cắt, chịu xoắn và tương tác uốn - cắt của dầm bê tông cốt hỗn hợp khi sử dụng cốt đai GFRP/thép lai.
- Ứng dụng trí tuệ nhân tạo (Machine Learning/Deep Neural Networks) trên tập dữ liệu lớn để tối ưu hóa đa mục tiêu: Chi phí vòng đời (LCC) - Độ dẻo - Khả năng chịu lực.
Tác động và ảnh hưởng
Về mặt học thuật, công trình đóng góp bộ dữ liệu thực nghiệm chuẩn mực đầu tiên tại Việt Nam về dầm SGFRP, được kỳ vọng sẽ trở thành tài liệu tham khảo cốt lõi cho các nghiên cứu sinh, học viên cao học chuyên ngành Kỹ thuật Xây dựng công trình dân dụng và công nghiệp.
Về mặt chuyển đổi công nghiệp, việc ứng dụng cốt hỗn hợp SGFRP giúp ngành xây dựng ven biển giảm thiểu tới 60–70% chi phí duy tu, sửa chữa do ăn mòn cốt thép trong vòng đời 50–100 năm của công trình. Các nhà thầu và đơn vị sản xuất cấu kiện bê tông đúc sẵn (precast concrete) có thể ứng dụng ngay công nghệ này để sản xuất cọc cừ, tấm kè biển, dầm cầu cảng có tuổi thọ siêu bền.
Về mặt xã hội và môi trường, việc giảm thiểu sử dụng thép xây dựng và kéo dài tuổi thọ công trình giúp cắt giảm hàng triệu tấn phát thải CO₂ từ ngành công nghiệp luyện kim, đóng góp trực tiếp vào mục tiêu Net Zero Carbon của quốc gia.
Đối tượng hưởng lợi
- Nghiên cứu sinh và Giới học thuật: Tiếp cận hệ thống lý thuyết giải tích hoàn chỉnh, thuật toán số lặp phi tuyến và bộ dữ liệu thực nghiệm đo đạc chi tiết để phát triển các đề tài nghiên cứu chuyên sâu.
- Kỹ sư kết cấu và Đơn vị tư vấn thiết kế: Nhận được quy trình và các công thức thiết kế đóng, biểu đồ tra cứu hàm lượng $\mu_{\min}, \mu_{\max}, \mu_{\max,e}$ tương thích 100% với TCVN 5574:2018 để áp dụng vào các dự án thực tế mà không cần mò mẫm tiêu chuẩn nước ngoài.
- Chủ đầu tư công trình biển và Cơ quan quản lý hạ tầng: Sở hữu giải pháp công nghệ vật liệu bền vững, nâng cao tuổi thọ công trình từ 30 năm lên trên 75 năm trong môi trường biển, tiết kiệm hàng nghìn tỷ đồng chi phí bảo trì.
- Hội đồng Tiêu chuẩn Quốc gia: Có căn cứ khoa học và số liệu thực nghiệm tin cậy để xúc tiến xây dựng Tiêu chuẩn Quốc gia về thiết kế kết cấu bê tông cốt sợi composite hỗn hợp.
Câu hỏi chuyên sâu
- Đóng góp lý thuyết độc đáo nhất của luận án là gì?
Trả lời: Đó là việc xây dựng thành công mô hình cơ học vật liệu cốt hỗn hợp tương đương với đường cong ứng suất - biến dạng song tuyến hai giai đoạn ($\sigma_{sf} - \varepsilon_{sf}$) tích hợp trực tiếp vào hệ phương trình giải tích phi tuyến của TCVN 5574:2018, giúp nhận diện và giải quyết chính xác dạng phá hoại dẻo hữu hiệu ($\xi_{R,e} \le \xi \le \xi_{R,y}$) mà các mô hình quy đổi đơn nhất quốc tế bỏ sót.
- Đổi mới phương pháp luận so với các nghiên cứu trước đây thế nào?
Trả lời: So với Wenjun Qu (2009) quy đổi về thép và Ahmed El Refai (2015) quy đổi về FRP, luận án không dùng phép quy đổi tương đương hình học giản đơn mà phát triển giải thuật số chia dải vi phân kết hợp điều kiện biên biến dạng thực tế của từng thớ vật liệu, đạt độ chính xác cao với sai số dự báo $M_u < 8%$.
- Phát hiện nào gây bất ngờ và có ý nghĩa kỹ thuật lớn nhất?
Trả lời: Phát hiện về sự duy trì độ cứng sau chảy dẻo (post-yielding stiffness): Khi cốt thép đã chảy dẻo, cốt GFRP không chỉ giữ cho dầm không bị sụp đổ đột ngột mà còn tạo ra một "thềm tăng tải thứ hai" ổn định, cho phép dầm hấp thụ năng lượng biến dạng lớn trước khi bê tông vỡ nén, phá tan định kiến học thuật cho rằng cốt GFRP chỉ gây phá hoại giòn nguy hiểm.
- Luận án có cung cấp quy trình tái lập (Replication Protocol) không?
Trả lời: Toàn bộ quy trình chế tạo mẫu, tỷ lệ cấp phối bê tông mác 400, sơ đồ bố trí cảm biến biến dạng vi mô, giải thuật lặp và bảng số liệu thí nghiệm gốc đều được trình bày minh bạch, chi tiết trong Chương 3 và Phụ lục, cho phép các phòng thí nghiệm độc lập tái lập hoàn toàn kết quả.
- Lộ trình nghiên cứu 10 năm tới được phác thảo ra sao?
Trả lời: Lộ trình tập trung vào 4 trọng tâm: (1) Chuẩn hóa tiêu chuẩn thiết kế dầm SGFRP cho TCVN; (2) Mở rộng thí nghiệm mỏi và tải trọng động đất; (3) Tối ưu hóa cấu kiện chịu uốn - cắt kết hợp; (4) Phát triển phần mềm tự động hóa thiết kế tối ưu hóa kinh tế - kỹ thuật cho công trình biển đảo Việt Nam.
Kết luận
- Luận án đã giải quyết trọn vẹn bài toán cơ học chịu uốn của dầm bê tông cốt hỗn hợp thép và GFRP, lấp đầy khoảng trống học thuật về sự thiếu hụt quy trình tính toán dầm SGFRP theo hệ tiêu chuẩn Việt Nam TCVN 5574:2018.
- Thiết lập thành công mô hình giải tích ứng suất - biến dạng tương đương song tuyến tính cho cốt SGFRP, xác định chính xác ba miền phá hoại cơ học thông qua các chỉ số chiều cao vùng nén giới hạn $\xi_{R,y}$ và $\xi_{R,e}$.
- Xây dựng bộ công thức giải tích hoàn chỉnh dự đoán mô men nứt $M_{cr}$, mô men chảy $M_y$, mô men giới hạn $M_u$, độ võng $f$, bề rộng khe nứt $a_{crc}$ và chỉ số độ dẻo $DF$, được kiểm chứng chuẩn xác qua 30 mẫu dầm quốc tế và chuỗi thực nghiệm chuyên sâu tại phòng thí nghiệm LAS-XD125.
- Định lượng hóa miền tỷ lệ diện tích cốt dọc tối ưu $\omega = A_f / A_s \in [0,6; 1,5]$, đảm bảo hài hòa tuyệt đối giữa yêu cầu an toàn chịu lực (ULS), điều kiện sử dụng bình thường (SLS) và hiệu quả kinh tế.
- Mở ra 3 hướng nghiên cứu ứng dụng mũi nhọn: kết cấu bê tông cốt composite hỗn hợp chịu tải trọng động lực học, tối ưu hóa vòng đời công trình hạ tầng biển đảo, và hoàn thiện khung pháp lý tiêu chuẩn xây dựng quốc gia.