Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Trong lĩnh vực xây dựng công trình dân dụng và công nghiệp, kết cấu bê tông cốt thép ứng lực trước (BTCT ỨLT) được ứng dụng rộng rãi cho các hệ kết cấu dầm, sàn vượt khẩu độ lớn (dầm nhịp từ 10m đến 30m, sàn nhịp từ 7m đến 15m). Giải pháp này giúp giảm chiều dày kết cấu, giảm chiều cao tầng (tiết kiệm 30–70 cm/tầng), giảm 15–30% khối lượng bê tông, 30–80% khối lượng cốt thép và rút ngắn thời gian thi công từ 3–5 ngày cho mỗi sàn so với bê tông cốt thép thông thường.

Tuy nhiên, tại Việt Nam, thực tế ghi nhận một số công trình sàn BTCT ỨLT xuất hiện vết nứt sau 5–6 năm đưa vào khai thác hoặc xuất hiện ngay trong quá trình thi công sau khi tạo ứng lực trước. Nguyên nhân chủ yếu xuất phát từ sự sai khác giữa các thông số giả định tiền định trong giai đoạn thiết kế với các thông số thực tế trong quá trình thi công và khai thác (như hệ số ma sát giữa cáp và ống lồng, hệ số dao động đường cong cáp, độ tụt neo, cường độ và mô đun đàn hồi thực tế của vật liệu, độ ẩm, nhiệt độ môi trường, hiện tượng chùng ứng suất và sai số độ võng quỹ đạo cáp). Các sai lệch này dẫn đến việc tính toán sai ứng suất hữu hiệu trong tao cáp, khiến lượng cáp bố trí không phù hợp, làm ứng suất kéo/nén trong bê tông vượt giới hạn cho phép và hình thành vết nứt.

Khoảng trống nghiên cứu được tác giả xác định là: trong nước chưa có công trình nào nghiên cứu một cách toàn diện việc thiết lập hàm công năng tổng quát tích hợp các yếu tố tổn hao ứng suất ngẫu nhiên và nhận dạng các phân bố xác suất thực nghiệm của vật liệu, thi công tại Việt Nam để đánh giá độ tin cậy (xác suất an toàn) theo tiêu chí kiểm soát nứt của kết cấu sàn BTCT ỨLT căng sau có bám dính.

Mục tiêu nghiên cứu

  1. Nghiên cứu tổng quan và lựa chọn tiêu chuẩn tính toán thiết kế sàn BTCT ỨLT phù hợp với điều kiện Việt Nam, xây dựng quy trình tính toán ứng suất bê tông cho hệ sàn nhiều nhịp.
  2. Thiết lập hàm công năng giới hạn theo tiêu chí kiểm soát vết nứt của sàn BTCT ỨLT căng sau có bám dính, có xét đến các thành phần tổn hao ứng suất tức thời và dài hạn.
  3. Nhận dạng quy luật phân bố xác suất của các biến ngẫu nhiên liên quan đến vật liệu, hình học và thi công dựa trên số liệu thực nghiệm thu thập tại các công trình thực tế ở Việt Nam.
  4. Xây dựng chương trình tính toán độ tin cậy bằng phương pháp mô phỏng Monte Carlo; khảo sát định lượng mức độ ảnh hưởng của các thành phần tổn hao ứng suất đến xác suất an toàn của sàn; lập bảng tra quan hệ kỹ thuật – kinh tế giữa chiều dày sàn, số lượng tao cáp, độ tin cậy và giá thành.

Đối tượng và phạm vi nghiên cứu

  • Đối tượng nghiên cứu: Độ tin cậy (xác suất an toàn chống nứt) của kết cấu sàn phẳng bê tông cốt thép ứng lực trước căng sau có bám dính (CSBD).
  • Phạm vi không gian: Khảo sát các công trình nhà cao tầng dân dụng (văn phòng và chung cư) tại Việt Nam với kết cấu sàn 3 nhịp có các khẩu độ điển hình gồm $7{,}5\text{ m} \times 7{,}5\text{ m}$, $9{,}0\text{ m} \times 9{,}0\text{ m}$ và $12{,}0\text{ m} \times 12{,}0\text{ m}$.
  • Phạm vi tiêu chuẩn và thời gian: Nghiên cứu đối chiếu giữa các hệ thống tiêu chuẩn thiết kế gồm AS 3600-2009 (Úc), ACI 318-08 (Mỹ), BS EN 1992 (Anh/Châu Âu) và TCVN 5574:2012; số liệu thực nghiệm thi công và vật liệu được tổng hợp từ các dự án thực tế tại Việt Nam đến năm 2014.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu ngoài nước

  • Dự đoán tổn hao ứng suất và quan trắc thực nghiệm: Scanlon (1995) công bố nghiên cứu về dự đoán tổn hao ứng suất trong dầm cầu BTCT ỨLT, chỉ ra ảnh hưởng của các biến ngẫu nhiên (mô đun đàn hồi, từ biến, co ngót, nhiệt độ thủy hóa bê tông) đến tổn hao dài hạn và đề xuất giảm thiểu bất định bằng đo đạc tổn hao sớm. Chunguang Lan, Zhi Zhou, Jinping Ou (2011, 2012) nghiên cứu đo ứng suất cáp ỨLT bằng cảm biến quang học Brillouin.
  • Độ tin cậy kết cấu ỨLT chịu ăn mòn và suy thoái: Nessim công bố nghiên cứu về độ tin cậy của tấm bê tông ỨLT cáp không bám dính chịu ăn mòn do ẩm. Sigit Darmawan và Mark G. Stewart nghiên cứu phân tích độ tin cậy dầm cầu BTCT ỨLT bị ăn mòn cốt thép và cáp do xâm nhập ion clorua theo trạng thái giới hạn uốn và giới hạn sử dụng. Pillai và cộng sự phân tích ảnh hưởng ăn mòn cáp trong ống lồng khi vữa bơm không lấp kín khoảng trống đến độ tin cậy cầu BTCT ỨLT bằng mô phỏng Monte Carlo.
  • Độ tin cậy trong điều kiện đặc biệt và công trình hạt nhân: Christopher D. Eamon và Elin Jensen (2012) nghiên cứu độ tin cậy dầm BTCT ỨLT chịu tác động hỏa hoạn, chứng minh độ tin cậy suy giảm nhanh theo thời gian tiếp xúc lửa. Fabio Biondini, Franco Bontempi, Dan M. Frangopol, Pier Giorgio Malerba (2004) khảo sát độ tin cậy cầu BTCT ỨLT đúc hẫng bằng Monte Carlo theo bài toán nứt và biến dạng. M. Pandey (1997) và Patrik Andersson, Martin Hansson, Sven Thelandersson (2008) nghiên cứu ảnh hưởng tổn hao ứng suất dài hạn (co ngót, từ biến, chùng ứng suất) đến độ tin cậy vỏ lò phản ứng hạt nhân tại Thụy Điển.
  • Tối ưu hóa thiết kế dựa trên độ tin cậy: Dan M. Frangopol nghiên cứu kết hợp tối ưu hóa và độ tin cậy cho dầm ỨLT. Samer Barakat, Khaldoon Bani-Hani, Mohammed Q. Taha (2004) nghiên cứu tối ưu hóa đa mục tiêu kết cấu BTCT ỨLT (giảm thiểu chi phí, tối đa hóa chỉ số độ tin cậy độ bền uốn và độ bền kéo). Vincent nghiên cứu dự đoán biến dạng sàn BTCT ỨLT nhà cao tầng tại Đại học Adelaide bằng Monte Carlo. Thamarie Jayasinghe (2011) tại Đại học RMIT nghiên cứu dự báo biến dạng dầm sàn ỨLT căng sau theo tiêu chuẩn AS 3600-2009, tổng hợp ảnh hưởng từ biến, co ngót và cốt thép chịu nén theo Washa & Fluck (1952) và Stevens (1972). Ed Cross và cộng sự (Austress Freyssinet) phân tích bài toán kỹ thuật - kinh tế giữa chiều dày sàn, số lượng cáp và rung động sàn.

Các hướng nghiên cứu trong nước

  • Độ tin cậy công trình xây dựng và kết cấu: Nguyễn Xuân Chính (2003) nghiên cứu phương pháp đánh giá độ tin cậy theo chỉ số tin cậy $\beta$ và xác suất không an toàn $P_f$ qua 3 ví dụ (khung bê tông 3 tầng 2 nhịp, khán đài 12.000 chỗ, dầm consol). Nghiêm Quang Hà, Nguyễn Quang Huy (2003) nghiên cứu an toàn dàn lưới không gian. Nguyễn Đình Xân, Nguyễn Thạc Vũ (2003) nghiên cứu độ tin cậy kết cấu hỗn hợp cung thể thao. Lê Kiều (2005) nghiên cứu đảm bảo tuổi thọ công trình từ khâu thiết kế. Nguyễn Tuấn Anh (2006) đánh giá tuổi thọ kết cấu thép bị ăn mòn. Nguyễn Văn Phó, Nguyễn Xuân Chính, Tạ Thanh Vân (2006) nghiên cứu phương pháp đánh giá độ tin cậy công trình. Lê Xuân Huỳnh, Lê Công Duy (2006) nghiên cứu độ tin cậy mờ kết cấu khung. Phạm Đức Cương (2007) tính toán cấu kiện BTCT cơ bản theo độ tin cậy. Nguyễn Hùng Tuấn, Lê Xuân Huỳnh (2009) đánh giá độ tin cậy nhà chung cư theo mô hình mờ. Bùi Đức Năng (2010) bảo vệ luận án tiến sĩ về tính xác suất không hỏng của hệ thanh bằng mô phỏng Monte Carlo có xét biến ngẫu nhiên vật liệu và hình học.
  • Độ tin cậy công trình chuyên ngành giao thông, thủy lợi, biển: Nguyễn Văn Huân, Phùng Vĩnh An (2003) nghiên cứu cống dưới đê; Mai Châu Anh, Lê Ngọc Thạch (2006) nghiên cứu vỏ thoải chữ nhật; Đỗ Ngọc Tú (2006) nghiên cứu cột điện chịu bão; Đỗ Văn Đệ (2007, 2010) nghiên cứu công trình thủy và dàn chịu tải trọng sóng; Nguyễn Hoàng Long, Phạm Cao Thăng (2008) ứng dụng Monte Carlo tính độ tin cậy mặt đường cứng ô tô; Nguyễn Vi và cộng sự (2009, 2010, 2012) nghiên cứu độ tin cậy công trình cảng, mái dốc, tường chắn cứng và tường cừ một tầng neo; Bùi Đức Chính (2010) đánh giá độ tin cậy cầu BTCT đang khai thác; Phạm Khắc Hùng (2010) nghiên cứu độ tin cậy kết cấu công trình biển cố định bằng thép.
  • Tối ưu hóa thiết kế và tổn hao ứng suất kết cấu ỨLT: Lê Xuân Huỳnh, Hoàng Xuân Long (2002) nghiên cứu quan hệ chi phí và độ tin cậy; Nguyễn Trọng Hà (2009) phân tích độ nhạy độ tin cậy; Nguyễn Hữu Lộc, Đinh Lê Cao Kỳ (2009) tối ưu hóa kết cấu dựa trên độ tin cậy bằng giải thuật di truyền; Hoàng Quang Nhu (2007, 2008) nghiên cứu tổn hao ứng suất do từ biến và co ngót; Nguyễn Xuân Khang, Nguyễn Văn Thuyên (2007) nghiên cứu kiểm soát lực kéo và độ giãn dài; Nguyễn Thanh Hưng (2009) khảo sát độ tin cậy sàn phẳng bê tông dự ứng lực theo bài toán ứng suất tức thời và dài hạn nhưng sử dụng các giá trị độ lệch giả định và chưa thiết lập hàm công năng tổng quát có xét đầy đủ thành phần tổn hao ứng suất.

Vị trí và khoảng trống luận án giải quyết

Luận án kế thừa cơ sở lý thuyết độ tin cậy kết cấu và phương pháp mô phỏng xác suất, khắc phục hạn chế của các nghiên cứu trước bằng cách:

  1. Xây dựng hàm công năng trạng thái giới hạn kiểm soát vết nứt hoàn chỉnh cho sàn BTCT ỨLT, tích hợp trực tiếp các công thức giải tích tính toán tổn hao ứng suất tức thời (ma sát, tụt neo) và dài hạn (co ngót, từ biến bê tông, chùng ứng suất cáp) theo tiêu chuẩn AS 3600-2009.
  2. Thu thập dữ liệu thực tế tại hiện trường thi công ở Việt Nam để nhận dạng quy luật phân bố xác suất của các biến ngẫu nhiên bằng phương pháp phi tham số (Kernel) và xấp xỉ hàm giải tích, thay vì áp đặt phân bố chuẩn giả định.
  3. Thiết lập mối quan hệ định lượng giữa tổn hao ứng suất, các sai số thi công và mức độ an toàn chống nứt của kết cấu, đồng thời xây dựng công cụ tra cứu tối ưu hóa kỹ thuật - kinh tế phục vụ trực tiếp công tác thiết kế thực hành.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khái niệm và lý thuyết độ tin cậy kết cấu

  • Định nghĩa độ tin cậy: Khả năng của kết cấu hoàn thành các tiêu chuẩn thiết kế và yêu cầu công năng quy định trong suốt tuổi thọ thiết kế (theo Palle Thoft-Christensen, Yoshisada Murotsu; Eurocode 0; ISO 2394; Tiêu chuẩn độ tin cậy kết cấu xây dựng và nền của Nga được Bộ Xây dựng thông qua).
  • Hàm công năng (Performance Function): Xác định trạng thái làm việc của kết cấu thông qua hàm: $$Z = g(X_1, X_2, \dots, X_n)$$ Trong đó: $Z > 0$ là miền an toàn; $Z < 0$ là miền không an toàn (hư hại/nứt); $Z = 0$ là mặt giới hạn. Xác suất an toàn được định nghĩa: $$P_s = P(Z > 0) = 1 - P_f$$ Với $P_f$ là xác suất không an toàn ($P_f = P(Z \le 0)$).

Cơ sở tính toán tổn hao ứng suất và trạng thái ứng suất

  • Tổn hao ứng suất tức thời:
    • Tổn hao do ma sát giữa cáp và ống lồng ($\Delta \sigma_{fr}$), phụ thuộc hệ số ma sát $\mu$ và hệ số dao động đường cong $k$.
    • Tổn hao do biến dạng đàn hồi của bê tông ($\Delta \sigma_{el}$).
    • Tổn hao do độ tụt neo ($\Delta L$).
  • Tổn hao ứng suất dài hạn:
    • Tổn hao do co ngót của bê tông ($\Delta \sigma_{sh}$), phụ thuộc thời gian dưỡng hộ ẩm $t_s$ và độ ẩm môi trường.
    • Tổn hao do từ biến của bê tông ($\Delta \sigma_{cr}$), phụ thuộc tuổi bê tông khi chịu tải $t_0$ và hệ số từ biến.
    • Tổn hao do chùng ứng suất của cốt thép ($\Delta \sigma_{re}$), phụ thuộc nhiệt độ môi trường $T$, chùng cơ bản $R_b$ và tỷ số mức ứng suất kéo ban đầu $P_j/f_p$.
  • Mô hình tính toán cơ học sàn: Sử dụng phương pháp phân phối mô men / khung tương đương để xác định nội lực dải cột, dải giữa nhịp và tính toán ứng suất thớ dưới/thớ trên của bê tông ở giai đoạn sử dụng dài lâu.

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

  1. Phương pháp phân tích lý thuyết và so sánh tiêu chuẩn: Phân tích, so sánh các công thức giải tích tính tổn hao ứng suất và kiểm tra ứng suất nứt giữa AS 3600-2009, ACI 318-08, BS EN 1992 và TCVN 5574:2012 để lựa chọn tiêu chuẩn thiết kế phù hợp.
  2. Phương pháp thống kê phi tham số (Kernel Density Estimation): Ước lượng hàm mật độ xác suất thực nghiệm từ tập dữ liệu mẫu quan trắc tại hiện trường: $$\hat{f}h(x) = \frac{1}{nh} \sum{i=1}^{n} K\left(\frac{x - X_i}{h}\right)$$ Với $K(u)$ là hàm nhân (Gaussian/Epanechnikov) và $h$ là độ rộng cửa sổ (bandwidth).
  3. Phương pháp bình phương tối thiểu: Xấp xỉ hàm mật độ thực nghiệm tìm được từ phương pháp Kernel sang các dạng hàm phân bố xác suất lý thuyết liên tục (Chuẩn - Normal, Mũ - Exponential, Gamma, Lognormal) và kiểm định thống kê.
  4. Phương pháp xác suất - hình học Hasofer-Lind và mô phỏng Monte Carlo: Sử dụng phương pháp Hasofer-Lind trong không gian biến chuẩn hóa và thuật toán mô phỏng ngẫu nhiên Monte Carlo với số lần lặp lớn để tính xác suất an toàn $P_s$, khảo sát sự hội tụ của nghiệm.

Nguồn dữ liệu thực nghiệm

Số liệu phục vụ nhận dạng biến ngẫu nhiên được thu thập trực tiếp từ các dự án thực tế do Công ty Cổ phần Đầu tư và Công nghệ Xây dựng IBST thi công, bao gồm:

  • Tiết diện ngang của tao cáp ứng lực ($A_p$).
  • Mô đun đàn hồi của vật liệu làm cáp ($E_p$).
  • Độ tụt nêm neo tại đầu neo ($\Delta L$).
  • Sai số đặc trưng hình học thi công / độ võng của tuyến cáp ($L$).

Nội dung chính theo từng chương

Chương I: Tổng quan

Chương I trình bày tổng quan về kết cấu BTCT ỨLT, phân loại công nghệ và tình hình ứng dụng tại Việt Nam; điểm luận các công trình nghiên cứu trong và ngoài nước về độ tin cậy kết cấu; phân tích các vấn đề tồn tại trong thiết kế, thi công dẫn đến hư hại sàn.

So sánh công nghệ căng sau có bám dính (CSBD) và không bám dính

Luận án tiến hành so sánh toàn diện hai giải pháp thi công căng sau:

Tiêu chí Căng sau có bám dính (CSBD) Căng sau không bám dính
Cách tạo ỨLT Kết hợp lực bám dính giữa cáp với bê tông (qua vữa bơm xi măng áp lực cao lấp kín ống ghen tôn mạ kẽm/nhựa) và phản lực ép mặt đầu neo. Chỉ sử dụng phản lực ép mặt đầu neo truyền lực vào bê tông; không có lực bám dính dọc thân cáp.
Ứng xử khi hỏa hoạn Cáp được bảo vệ bởi lớp vữa xi măng cường độ cao trong ống lồng, ứng suất chỉ suy giảm khi toàn bộ tiết diện bê tông bị phá hoại do nhiệt. Cáp giãn tự do giữa hai đầu neo khi nhiệt độ tăng; lực căng suy giảm rất nhanh, tiềm ẩn nguy cơ sụp đổ kết cấu.
Mức độ an toàn đầu neo Đầu neo được bịt vữa xi măng cường độ cao; nếu đầu neo gặp sự cố, cáp vẫn làm việc tốt nhờ lực bám dính dọc chiều dài. Đầu neo bọc mỡ trung tính; nếu đầu neo hỏng, toàn bộ tao cáp mất hoàn toàn ứng lực dọc suốt chiều dài.
Khả năng xử lý sự cố Khó xử lý đứt cáp sau khi đã bơm vữa, nhưng trong quá trình kéo căng nếu đứt/tụt cáp có thể rút và thay thế trong ống lồng. Khi xảy ra đứt cáp trong quá trình thi công hoặc sử dụng, việc khắc phục rất khó khăn do phải đục phá sàn dọc tuyến cáp.
Bố trí và thi công Cáp gom thành bó đặt trong ống lồng; khoảng cách giữa các bó cáp thưa, dễ tránh các lỗ mở kỹ thuật; dễ tạo profile đường cong. Cáp bố trí đơn lẻ; mật độ cáp dày đặc, khó bố trí tránh lỗ mở kỹ thuật; thi công kéo dài và khó kiểm soát dung sai profile.

Luận điểm rút ra: Sàn BTCT ỨLT căng sau có bám dính sở hữu nhiều ưu điểm vượt trội về an toàn chịu lực lâu dài, chống cháy và kiểm soát nứt nên được ưu tiên áp dụng cho các công trình cao tầng lớn (như Tòa tháp BIDV, Tòa tháp Keangnam, Usilk City, The Pride, Siêu thị Mê Linh Plaza).

Chương II: Nghiên cứu lựa chọn tiêu chuẩn tính toán và tính ứng suất trong bê tông cho trường hợp thiết kế tổng quát

Chương II đi sâu phân tích lý thuyết tính toán tổn hao ứng suất theo các tiêu chuẩn quốc tế và Việt Nam; xây dựng ví dụ so sánh định lượng; lựa chọn tiêu chuẩn AS 3600-2009 làm chuẩn tính toán; thiết lập thuật toán xác định nội lực và ứng suất bê tông sàn nhiều nhịp.

+-------------------------------------------------------------------------+
|                  QUY TRÌNH THIẾT KẾ & KIỂM SOÁT ỨNG SUẤT                |
|                           (Theo AS 3600-2009)                           |
+-------------------------------------------------------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       | 1. Lựa chọn sơ bộ: Cường độ bê tông (f'c), chiều dày sàn  |
       |    (Ds), loại cáp, ứng suất căng ban đầu (sigma_pj)       |
       +-----------------------------------------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       | 2. Phân tích nội lực: Phương pháp dải thiết kế, phân phối |
       |    mô men ngoại tải (M0 -> M_gối, M_nhịp -> Dải cột/giữa) |
       +-----------------------------------------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       | 3. Tính toán tổn hao ứng suất:                           |
       |    - Tức thời: Ma sát (mu, k), tụt neo (Delta L)          |
       |    - Dài hạn: Co ngót (eps_cs), từ biến (phi_cc), chùng   |
       +-----------------------------------------------------------+
                                     |
                                     v
       +-----------------------------------------------------------+
       | 4. Kiểm tra ứng suất bê tông giai đoạn sử dụng lâu dài:   |
       |    sigma_bt = (M_ngoại tải - M_bal) / W - P_eff / A_c     |
       |    Điều kiện kiểm soát nứt: sigma_bt <= [sigma_kéo]       |
       +-----------------------------------------------------------+

So sánh tổn hao ứng suất giữa các tiêu chuẩn

  • Tổn hao do ma sát và tụt neo: So sánh theo AS 3600, ACI 318, BS EN 1992 cho thấy mức độ tổn hao phụ thuộc chặt chẽ vào hệ số ma sát $\mu$ và hệ số dao động đường cong $k$.
  • Tổn hao do co ngót bê tông: Tính toán theo ACI 318-08 và AS 3600-2009 với các giá trị thời gian dưỡng hộ ẩm $t_s$ khác nhau.
  • Tổn hao do từ biến bê tông: So sánh giữa AS 3600-2009 và BS EN 1992 với các giá trị tuổi bê tông khi chịu tải $t_0$ (khảo sát tại mốc $t_0 = 180\text{ ngày}$).
  • Tổn hao do chùng ứng suất cốt thép: Đối chiếu TCVN 5574:2012, ACI 318-08, AS 3600-2009 và BS EN 1992 theo các tỷ số $f_{si}/f_{pu}$, $P_j/f_p$, chùng cơ bản $R_b$ và nhiệt độ $T$.

Thuật toán tính ứng suất sàn n nhịp

Thiết lập bài toán giải tích phân phối mô men từ mô men tổng $M_0$ thành mô men tại gối $M_N$ và mô men tại nhịp $M_M$, sau đó phân bổ lên dải cột (Column Strip) và dải giữa nhịp (Middle Strip). Xác định lực căng hiệu dụng còn lại trong cáp sau khi trừ tổn hao tức thời và tổn hao dài hạn, từ đó tính toán chính xác ứng suất kéo/nén thớ mặt ngoài bê tông ở giai đoạn làm việc lâu dài.

Chương III: Xây dựng hàm công năng theo tiêu chí kiểm soát vết nứt sàn bê tông cốt thép ứng lực trước và nhận dạng các biến ngẫu nhiên qua các số liệu thực nghiệm

Chương III thiết lập hàm trạng thái giới hạn kiểm soát nứt cho hệ sàn 3 nhịp; trình bày phương pháp Kernel ước lượng hàm mật độ phi tham số; thực hiện nhận dạng thống kê 4 biến ngẫu nhiên cốt lõi từ số liệu công trường tại Việt Nam.

Thiết lập hàm công năng kiểm soát vết nứt

Hàm công năng được xây dựng dựa trên điều kiện ứng suất kéo lớn nhất sinh ra ở thớ bê tông chịu kéo không vượt quá giới hạn kéo cho phép của bê tông theo AS 3600-2009: $$Z = g(X) = [\sigma_{kt}] - \sigma_{bt}(A_p, E_p, \Delta L, L, P_{pj}, T, RH, R_b, \dots)$$ Trong đó:

  • $[\sigma_{kt}] = 0{,}6\sqrt{f'_c}$ (hoặc theo cấp độ kiểm soát nứt của tiêu chuẩn).
  • $\sigma_{bt}$ là ứng suất tính toán thực tế tại thớ kéo dưới tác dụng của ngoại tải và lực nén hiệu dụng $P_{eff}$ do cáp truyền vào bê tông (đã trừ toàn bộ tổn hao ứng suất ngẫu nhiên).

Quy trình nhận dạng biến ngẫu nhiên bằng phương pháp Kernel

  1. Lập bảng tần suất thực nghiệm từ tập số liệu đo đạc hiện trường.
  2. Dùng hàm nhân Gaussian xác định đường cong mật độ phân bố xác suất thực nghiệm liên tục $\hat{f}(x)$.
  3. Áp dụng phương pháp bình phương tối thiểu để làm mịn và khớp đường cong thực nghiệm với các dạng hàm phân bố lý thuyết chuẩn tắc: Chuẩn (Normal), Mũ (Exponential), Gamma.
       Số liệu thực nghiệm (IBST)
   [Ap, Ep, Tụt neo Delta L, Sai số L]
                  |
                  v
       +----------------------+
       |  Tổ chức đồ tần suất |
       +----------------------+
                  |
                  v
       +------------------------------------------+
       | Phương pháp Kernel (Hàm nhân Gaussian)   |
       |  Ước lượng hàm mật độ thực nghiệm liên tục|
       +------------------------------------------+
                  |
                  v
       +------------------------------------------+
       | Khớp hàm bằng Bình phương tối thiểu     |
       |  (Normal, Exponential, Gamma, Lognormal) |
       +------------------------------------------+
                  |
                  v
       Kỳ vọng (E) & Độ lệch chuẩn (Std)
          phục vụ Mô phỏng Monte Carlo

Kết quả nhận dạng các biến ngẫu nhiên thi công tại Việt Nam

Dựa trên số liệu thống kê thu thập từ IBST:

  • Tiết diện ngang cáp ($A_p$): Có kỳ vọng thực tế sát với giá trị danh định của nhà sản xuất, độ tán xạ tuân theo phân bố Chuẩn/Lognormal.
  • Mô đun đàn hồi cáp ($E_p$): Xác định giá trị kỳ vọng và độ lệch chuẩn đặc trưng cho các lô thép nhập khẩu thi công tại Việt Nam.
  • Độ tụt neo ($\Delta L$): Khảo sát biến động độ trượt nêm neo khi khóa cáp, phân bố thực nghiệm phản ánh tay nghề công nhân và chất lượng thiết bị kích kéo.
  • Sai số đặc trưng thi công quỹ đạo cáp ($L$): Biến thiên độ võng và tọa độ điểm đặt cáp so với bản vẽ thiết kế hình học.

Chương IV: Xây dựng chương trình tính độ tin cậy và khảo sát ảnh hưởng của tổn hao ứng suất đến độ tin cậy sàn bê tông cốt thép ứng lực trước căng sau có bám dính

Chương IV xây dựng phần mềm tính toán độ tin cậy kết hợp thuật toán Monte Carlo và hàm công năng cơ học; khảo sát định lượng 6 trường hợp thiết kế cụ thể; đánh giá độ nhạy của các tham số tổn hao; lập bảng tra kinh tế - kỹ thuật phục vụ thiết kế tối ưu.

Mô hình khảo sát thực nghiệm số

Khảo sát 6 mô hình sàn 3 nhịp phẳng chịu tải trọng văn phòng và chung cư:

  1. Sàn văn phòng nhịp $7{,}5\text{ m} \times 7{,}5\text{ m} \times 7{,}5\text{ m}$: Chiều dày $D_s = 200\text{ mm}$, số tao cáp $N = 15$.
  2. Sàn chung cư nhịp $7{,}5\text{ m} \times 7{,}5\text{ m} \times 7{,}5\text{ m}$: Chiều dày $D_s = 210\text{ mm}$, số tao cáp $N = 19$.
  3. Sàn văn phòng nhịp $9{,}0\text{ m} \times 9{,}0\text{ m} \times 9{,}0\text{ m}$: Chiều dày $D_s = 230\text{ mm}$, số tao cáp $N = 24$.
  4. Sàn chung cư nhịp $9{,}0\text{ m} \times 9{,}0\text{ m} \times 9{,}0\text{ m}$: Chiều dày $D_s = 240\text{ mm}$, số tao cáp $N = 30$.
  5. Sàn văn phòng nhịp $12{,}0\text{ m} \times 12{,}0\text{ m} \times 12{,}0\text{ m}$: Chiều dày $D_s = 330\text{ mm}$, số tao cáp $N = 50$.
  6. Sàn chung cư nhịp $12{,}0\text{ m} \times 12{,}0\text{ m} \times 12{,}0\text{ m}$: Chiều dày $D_s = 340\text{ mm}$, số tao cáp $N = 60$.

Khảo sát ảnh hưởng của các yếu tố biến động đến xác suất an toàn chống nứt ($P_s$)

  1. Ảnh hưởng của lực kéo căng cáp ban đầu ($P_{pj}$): Khảo sát tại 3 mức $95%$, $100%$ và $105%$ ứng suất cho phép. Kết quả cho thấy khi lực căng cáp tăng, xác suất an toàn $P_s$ tăng mạnh do lực nén trước triệt tiêu ứng suất kéo trong bê tông.
  2. Ảnh hưởng của nhiệt độ môi trường ($T$): Nhiệt độ môi trường tăng làm gia tăng tốc độ chùng ứng suất của thép và tổn hao từ biến/co ngót của bê tông, làm giảm lực căng hiệu dụng và làm giảm xác suất an toàn $P_s$.
  3. Ảnh hưởng của độ ẩm môi trường ($RH$): Độ ẩm môi trường giảm làm tăng mạnh độ co ngót của bê tông ($\varepsilon_{cs}$), dẫn đến gia tăng tổn hao ứng suất dài hạn và làm suy giảm độ tin cậy chống nứt.
  4. Ảnh hưởng của chùng ứng suất cơ bản ($R_b$): Thép có độ chùng thấp (Low Relaxation) duy trì độ tin cậy cao hơn đáng kể so với thép có độ chùng thông thường.
  5. Ảnh hưởng của sai số độ võng quỹ đạo cáp trong thi công: Sai lệch vị trí đặt cáp tại mặt cắt giữa nhịp và gối làm thay đổi cánh tay đòn nội lực, gây sụt giảm xác suất an toàn $P_s$ nhanh chóng nếu dung sai thi công vượt quá giới hạn cho phép.

Đánh giá mối quan hệ kỹ thuật - kinh tế ($N - D_s - P_s$ - Giá thành)

Luận án xây dựng bảng tra và biểu đồ liên hệ giữa chiều dày sàn $D_s$, số lượng tao cáp $N$, xác suất an toàn $P_s$ và chỉ số giá thành ($10^3\text{ VNĐ/m}^2$) cho hệ sàn nhịp $12\text{ m}$ (nhà văn phòng và chung cư):

Nhóm thông số Phương án giảm chiều dày sàn ($D_s$ nhỏ) Phương án tăng chiều dày sàn ($D_s$ lớn)
Đặc điểm kỹ thuật Tăng số lượng tao cáp $N$ để bù độ võng và ứng suất; độ mảnh kết cấu cao, tăng dao động rung. Giảm số lượng tao cáp $N$; tăng khối lượng bê tông bản sàn và tĩnh tải bản thân.
Độ tin cậy ($P_s$) Nhạy cảm cao với sai số thi công và tổn hao dài hạn; xác suất an toàn cần được kiểm soát chặt. Xác suất an toàn ổn định hơn, độ võng nhỏ.
Giá thành ($10^3\text{ VNĐ/m}^2$) Tối ưu về mặt không gian kiến trúc nhưng chi phí cáp và neo tăng cao. Tiết kiệm chi phí cáp nhưng chi phí bê tông và kết cấu cột/móng chịu tải tăng.

Kết quả và những đóng góp mới

Đóng góp về mặt lý luận khoa học

  1. Thiết lập thành công hàm công năng trạng thái giới hạn kiểm soát vết nứt hoàn chỉnh cho kết cấu sàn BTCT ỨLT căng sau có bám dính nhiều nhịp, tích hợp các công thức giải tích tính toán tổn hao ứng suất tức thời và dài hạn theo tiêu chuẩn AS 3600-2009.
  2. Ứng dụng phương pháp thống kê phi tham số Kernel kết hợp bình phương tối thiểu để nhận dạng chính xác quy luật phân bố xác suất thực tế của các biến ngẫu nhiên liên quan đến vật liệu và thi công tại Việt Nam ($A_p, E_p, \Delta L, L$), khắc phục việc áp đặt phân bố chuẩn thuần túy lý thuyết.
  3. Ứng dụng thuật toán mô phỏng ngẫu nhiên Monte Carlo để định lượng xác suất an toàn chống nứt của sàn BTCT ỨLT, hình thành phương pháp luận đánh giá an toàn kết cấu dựa trên lý thuyết độ tin cậy trong điều kiện các thông số đầu vào có tính bất định cao.

Đóng góp về mặt thực tiễn

  1. Cung cấp bộ số liệu thống kê thực nghiệm có giá trị về đặc trưng cơ lý của vật liệu cáp ứng lực và sai số thi công thu thập tại các công trình thực tế ở Việt Nam.
  2. Xây dựng chương trình tính toán tự động độ tin cậy của sàn BTCT ỨLT, cho phép kỹ sư đánh giá định lượng rủi ro xuất hiện vết nứt ngay từ giai đoạn thiết kế.
  3. Thiết lập các biểu đồ độ nhạy và bảng tra quan hệ giữa chiều dày sàn $D_s$, số lượng tao cáp $N$, xác suất an toàn $P_s$ và giá thành đơn vị ($10^3\text{ VNĐ/m}^2$), giúp các đơn vị tư vấn thiết kế và chủ đầu tư lựa chọn phương án kết cấu tối ưu về kỹ thuật và kinh tế.
  4. Đưa ra các chỉ dẫn kỹ thuật cụ thể trong thi công: kiểm soát lực căng kéo cáp, kiểm soát độ tụt neo và độ chính xác quỹ đạo cáp nhằm hạn chế tối đa nguy cơ nứt sàn trong quá trình khai thác lâu dài.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế

  • Số liệu thực nghiệm dùng để nhận dạng các biến ngẫu nhiên ($A_p, E_p, \Delta L, L$) được thu thập chủ yếu từ các dự án do một đơn vị chuyên ngành thi công (IBST), chưa bao phủ toàn bộ các vùng miền khí hậu và các nhà thầu khác nhau tại Việt Nam.
  • Hàm công năng kiểm soát nứt mới chỉ tập trung khảo sát hệ kết cấu sàn phẳng nhiều nhịp đối xứng chịu tải trọng tĩnh phân bố đều trong nhà cao tầng dân dụng, chưa xét tới các dạng sàn phức tạp (sàn có dầm nấm, sàn có lỗ mở lớn, sàn chịu tải trọng động đất hoặc rung động máy móc).
  • Chưa xét đến sự tương tác đồng thời giữa tổn hao ứng suất do từ biến/co ngót với quá trình xâm thực ăn mòn cốt thép trong môi trường biển mặn dài hạn.

Hướng nghiên cứu tiếp

  • Mở rộng quy mô điều tra và cập nhật cơ sở dữ liệu ngẫu nhiên về vật liệu, thiết bị kích kéo và tay nghề thi công trên phạm vi toàn quốc.
  • Xây dựng hàm công năng độ tin cậy đa mục tiêu, kết hợp đồng thời trạng thái giới hạn kiểm soát nứt, trạng thái giới hạn độ võng dài hạn và trạng thái giới hạn cực hạn chịu uốn/cắt của sàn.
  • Tích hợp mô hình tính toán độ tin cậy vào các phần mềm thiết kế kết cấu chuyên dụng và mô hình thông tin công trình (BIM) nhằm tự động hóa quy trình thiết kế tối ưu hóa độ tin cậy - chi phí vòng đời công trình.

Giá trị tham khảo

  • Đối với nghiên cứu sinh và học viên cao học: Cung cấp tài liệu tham khảo về phương pháp luận áp dụng lý thuyết độ tin cậy vào kết cấu BTCT ỨLT, kỹ thuật ước lượng phi tham số Kernel và thuật toán mô phỏng Monte Carlo trong kỹ thuật công trình xây dựng.
  • Đối với giảng viên và các viện nghiên cứu: Là tài liệu chuyên khảo trong giảng dạy các môn học Kết cấu bê tông cốt thép nâng cao, Độ tin cậy công trình và Tối ưu hóa kết cấu; làm cơ sở dữ liệu phục vụ biên soạn, hoàn thiện các tiêu chuẩn thiết kế và thi công nghiệm thu kết cấu ỨLT tại Việt Nam (như các dự thảo TC 03, TC 04, TC 51-05).
  • Đối với kỹ sư thiết kế và tư vấn giám sát: Giúp nắm vững bản chất các thành phần tổn hao ứng suất theo tiêu chuẩn AS 3600-2009; sử dụng bảng tra $N - D_s - P_s$ để lựa chọn giải pháp thiết kế hợp lý; nhận thức rõ tầm quan trọng của việc kiểm soát sai số thi công (độ tụt neo, profile cáp, lực căng kéo) tại công trường nhằm ngăn ngừa sự cố nứt sàn.

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

1. Vì sao sàn BTCT ỨLT căng sau có bám dính (CSBD) được ưu tiên ứng dụng hơn loại không bám dính trong nhà cao tầng?

Theo phân tích ở Chương I, sàn căng sau có bám dính có cáp được bọc trong ống ghen và bơm vữa xi măng áp lực cao lấp kín sau khi kéo căng. Điều này mang lại 3 ưu thế vượt trội: (1) Khả năng chịu hỏa hoạn cao vì cáp được lớp vữa bảo vệ; (2) Khi đầu neo gặp sự cố, cáp vẫn duy trì lực căng nhờ bám dính dọc thân; (3) Cáp gom thành bó thưa, dễ bố trí tránh các lỗ mở kỹ thuật và dễ kiểm soát hình học quỹ đạo cáp hơn so với cáp không bám dính phân tán.

2. Tiêu chuẩn thiết kế nào được luận án lựa chọn làm cơ sở tính toán chính và lý do?

Luận án lựa chọn tiêu chuẩn Úc AS 3600-2009 làm cơ sở tính toán chính (kết hợp đối chiếu ACI 318-08, BS EN 1992 và TCVN 5574:2012). Tiêu chuẩn này chỉ rõ trình tự thiết kế giải tích chi tiết, đưa ra các công thức tính toán tổn hao ứng suất do co ngót, từ biến, chùng ứng suất và quy định rõ ràng phương pháp dải thiết kế kiểm tra điều kiện ứng suất kiểm soát vết nứt trong bê tông ở giai đoạn sử dụng dài lâu.

3. Phương pháp Kernel được sử dụng trong luận án nhằm mục đích gì?

Phương pháp Kernel (ước lượng mật độ nhân) là công cụ thống kê phi tham số được tác giả dùng để xây dựng hàm mật độ xác suất thực nghiệm liên tục từ các tập số liệu đo đạc thực tế tại hiện trường ($A_p, E_p, \Delta L, L$). Phương pháp này không áp đặt phân bố giả định trước, giúp phản ánh trung thực quy luật biến động ngẫu nhiên của vật liệu và chất lượng thi công tại Việt Nam, sau đó mới dùng phương pháp bình phương tối thiểu để khớp sang các dạng hàm giải tích lý thuyết.

4. Kết quả khảo sát cho thấy những yếu tố thi công và môi trường nào ảnh hưởng mạnh nhất đến xác suất an toàn chống nứt ($P_s$)?

Kết quả ở Chương IV chỉ ra rằng:

  • Lực kéo căng ban đầu ($P_{pj}$): Nếu lực kéo đạt $105%$ ứng suất cho phép, xác suất an toàn tăng đáng kể; nếu kéo thiếu lực ($95%$), nguy cơ nứt sàn tăng cao.
  • Sai số quỹ đạo cáp ($L$) và độ tụt neo ($\Delta L$): Sai lệch chiều cao đặt cáp tại gối và nhịp làm suy giảm nhanh chóng ứng suất nén hiệu dụng.
  • Độ ẩm ($RH$) và nhiệt độ ($T$): Môi trường khô hanh (độ ẩm thấp) làm tăng co ngót, nhiệt độ cao làm tăng chùng ứng suất cáp, dẫn đến gia tăng tổn hao dài hạn và làm giảm xác suất an toàn $P_s$.

5. Bảng tra kỹ thuật - kinh tế ở Chương IV giúp giải quyết bài toán gì cho kỹ sư thiết kế?

Bảng tra cung cấp mối quan hệ giữa chiều dày sàn $D_s$, số lượng tao cáp $N$, xác suất an toàn $P_s$ và giá thành ($10^3\text{ VNĐ/m}^2$). Nó giúp kỹ sư không chỉ thiết kế sàn đạt yêu cầu chịu lực đơn thuần mà còn cân đối được bài toán kinh tế: chọn chiều dày sàn tối ưu để vừa đảm bảo độ tin cậy chống nứt theo tiêu chuẩn, vừa giảm thiểu chi phí vật liệu (cáp, neo, bê tông) và thuận lợi cho thi công.

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Chí Hiếu đã giải quyết trọn vẹn vấn đề đánh giá ảnh hưởng của các thành phần tổn hao ứng suất đến độ tin cậy của sàn bê tông cốt thép ứng lực trước căng sau có bám dính tại Việt Nam. Bằng việc kết hợp phương pháp dải thiết kế theo AS 3600-2009, ước lượng phi tham số Kernel từ dữ liệu hiện trường và mô phỏng xác suất Monte Carlo, tác giả đã thiết lập được hàm công năng kiểm soát nứt và lượng hóa tác động của các sai số thi công, điều kiện môi trường đến an toàn kết cấu. Các kết quả nghiên cứu, chương trình tính toán và bảng tra quan hệ kỹ thuật - kinh tế cung cấp cơ sở khoa học và công cụ thực tiễn giá trị cho công tác thiết kế, thi công và biên soạn tiêu chuẩn xây dựng nhà cao tầng.