Tổng quan về luận án

Nghiên cứu về kết cấu không gian mỏng, đặc biệt là mái vỏ thoải bê tông cốt thép (BTCT) cong hai chiều dương nhiều lớp, đại diện cho một bước phát triển vượt bậc trong kỹ thuật xây dựng công trình dân dụng và công nghiệp khẩu độ lớn. Trong bối cảnh kiến trúc hiện đại đòi hỏi tối ưu hóa vật liệu và mở rộng không gian không cột chống trung gian, mái vỏ mỏng mang lại lợi thế vượt trội khi "có trọng lượng bản thân giảm 20-30%, chi phí cốt thép tính trên $1\text{ m}^2$ mặt bằng bao che giảm 30-40%" so với các giải pháp kết cấu phẳng như dầm, dàn hoặc vòm truyền thống có cùng khẩu độ. Tuy nhiên, việc tích hợp cấu tạo nhiều lớp nhằm đáp ứng đồng thời yêu cầu chịu lực, chống thấm, cách nhiệt hoặc gia cường sửa chữa bằng bê tông cốt sợi kim loại phân tán (BTS) đặt ra những thách thức cơ học phức tạp về tính tương tác tiếp xúc và trạng thái ứng suất - biến dạng (ỨS-BD).

Khoảng trống nghiên cứu then chốt trong y văn xuất phát từ việc các nghiên cứu kinh điển trước đây của Ambarsumian (1966), Lê Thanh Huấn (1981, 2011), hay An-dray-ep và Nhi-me-rop-ski (2001) chủ yếu dựa trên giả thiết lý tưởng hóa rằng các lớp dính chặt hoàn toàn và xem như một lớp tương đương, nhưng chưa làm rõ các điều kiện biên, giới hạn tải trọng, khả năng phát sinh trượt tiếp xúc hoặc sự suy giảm liên kết cơ học giữa các lớp bêtông khác tính chất. Hơn nữa, các thực nghiệm trước đây thường dừng lại ở vật liệu tương tự quy mô nhỏ như kính hữu cơ (Lê Thanh Huấn), hoặc các mẫu vỏ một lớp kích thước nhỏ dưới $1.1\text{ m}$ (Sivakumar et al., 2015; Jeyashree et al., 2016; Praveenkumar et al., 2017). Luận án của tác giả Lâm Thanh Quang Khải (2019) tại Trường Đại học Kiến trúc Hà Nội dưới sự hướng dẫn khoa học của PGS. Lê Thanh Huấn và GS. Nguyễn Tiến Chương đã giải quyết triệt để khoảng trống này.

Công trình tập trung trả lời ba câu hỏi nghiên cứu và kiểm chứng ba giả thuyết khoa học:

  1. Câu hỏi 1 (Q1): Mái vỏ thoải cong hai chiều dương hai lớp gồm bê tông thường (BTT) và bê tông cốt sợi thép (BTS) có duy trì được trạng thái dính kết hoàn hảo (không trượt) dưới tác dụng của tải trọng phân bố đều trong giai đoạn đàn hồi trước khi xuất hiện vết nứt hay không?
    • Giả thuyết 1 (H1): Ứng suất tiếp tiếp xúc giữa hai lớp vật liệu đổ liên tục nằm trong giới hạn chịu cắt của liên kết bê tông, cho phép áp dụng lý thuyết vỏ một lớp tương đương mà không xảy ra hiện tượng trượt phân lớp.
  2. Câu hỏi 2 (Q2): Mức độ tương thích và sai lệch định lượng giữa lời giải giải tích giải hệ phương trình Vlasov, mô hình phần tử hữu hạn (PTHH) tuyến tính vỏ một lớp tương đương trên SAP2000, mô phỏng phi tuyến 3D trên ANSYS và kết quả đo đạc thực nghiệm trên mô hình thực $3\times 3\text{ m}$ là bao nhiêu?
    • Giả thuyết 2 (H2): Sự sai khác về ứng suất và độ võng giữa các phương pháp không vượt quá ngưỡng dung sai kỹ thuật cho phép ($< 10-15%$).
  3. Câu hỏi 3 (Q3): Các tham số hình học và vật liệu như tỷ lệ bề dày từng lớp ($h_1/h_2$), vị trí lớp bê tông sợi và hàm lượng sợi thép ảnh hưởng như thế nào đến sự phân bố nội lực và độ võng của kết cấu mái vỏ khẩu độ lớn?
    • Giả thuyết 3 (H3): Bố trí lớp bê tông cốt sợi ở vùng chịu kéo tối ưu hóa đáng kể khả năng chịu tải và giảm thiểu độ võng đỉnh vỏ.

Phạm vi nghiên cứu được xác lập nghiêm ngặt: khảo sát mái vỏ thoải cong hai chiều dương hai lớp và năm lớp có mặt bằng hình vuông kích thước hình học danh định $3\times 3\text{ m}$ (tỷ lệ vồng $f/a \le 1/5$, độ dốc $< 18^\circ$, bề dày $h_v \le a/500$), chịu tải trọng phân bố đều tĩnh, mở rộng khảo sát tham số cho kết cấu nhịp lớn quy mô thực tế $36\times 36\text{ m}$.

Literature Review và Positioning

Lý thuyết vỏ mỏng bắt nguồn từ các nghiên cứu nền tảng cuối thế kỷ XIX của Kirchhoff và Aron, sau đó được hệ thống hóa mạnh mẽ bởi các nhà khoa học Liên Xô từ những năm 1940. Đỉnh cao là công trình của Vlasov (1958, 1963) với tác phẩm "Lý thuyết tổng quát vỏ mỏng và ứng dụng trong kỹ thuật", thiết lập hệ phương trình đạo hàm riêng bậc 4 giải bài toán vỏ thoải bằng hai hàm ứng suất $\varphi$ và hàm chuyển vị $w$. Dựa trên nền tảng này, Timoshenko & Woinowsky-Krieger (1959), Flügge (1973), và Soare (1967) đã phát triển các phương pháp giải cho mái vỏ một lớp. Tại Việt Nam, Ngô Thế Phong (1985) đã sử dụng chuỗi lượng giác kép của Navier và chuỗi lượng giác đơn của Lévi để giải vỏ thoải chữ nhật, trong khi Lê Thanh Huấn (1981, 2011) và Baicốp V.N (1985) phát triển phương pháp điểm (bán giải tích) giải quyết các điều kiện biên phức tạp như dầm cứng, dàn phẳng và hệ cột biên.

Trường phái / Nhóm tác giả Phương pháp / Mô hình tiếp cận Đóng góp & Phát hiện chính Hạn chế & Tranh luận khoa học
Vlasov (1963), Ngô Thế Phong (1985), Lê Thanh Huấn (1981) Giải tích giải hệ phương trình đạo hàm riêng bậc 4; chuỗi Navier, Lévi, phương pháp điểm Xác lập hàm ứng suất $\varphi$ và chuyển vị $w$; giải quyết trạng thái phi mô men và mô men biên cho vỏ 1 lớp Giới hạn ở vỏ 1 lớp đẳng hướng; phức tạp khi xử lý liên kết ngàm thực tế hoặc vỏ nhiều lớp bất đối xứng
Ambarsumian (1966), Rao (1974), An-dray-ep & Nhi-me-rop-ski (2001) Lý thuyết vỏ nhiều lớp dị hướng; quy đổi về vỏ một lớp tương đương Chứng minh có thể quy đổi vỏ nhiều lớp về bề mặt trung gian nếu các lớp không trượt và làm việc đàn hồi Mặc định giả thiết không trượt; chưa xác định giới hạn tải trọng gây tách trượt và ứng xử tiếp xúc giữa các lớp
Sivakumar et al. (2015), Jeyashree et al. (2016) Thực nghiệm mô hình kích thước nhỏ ($68\times 68\text{ cm}$ đến $108\times 108\text{ cm}$); SAP2000 Đo đạc độ võng và quan hệ tải trọng - chuyển vị của vỏ một lớp bằng bao cát hoặc tải tập trung Quy mô mẫu quá nhỏ; không phản ánh được hiệu ứng quy mô (scale effect) và không xét tương tác đa lớp
Meleka et al. (2014), Verwimp et al. (2014) Thực nghiệm gia cường GFRP; mô hình vỏ lưới dệt (TRC) Đánh giá hiệu quả tăng cường vỏ có lỗ mở và phi tuyến vật liệu dệt Tập trung vào gia cố cục bộ một lớp; chưa mô hình hóa bêtông sợi phân tán trong hệ kết cấu vỏ thoải nhiều lớp

Cuộc tranh luận khoa học cốt lõi trong cơ học vỏ nhiều lớp xoay quanh việc liệu một kết cấu gồm các lớp vật liệu có mô đun đàn hồi khác nhau ($E_1 \ne E_2$) và cơ chế ứng xử kéo/nén dị biệt có thể được mô hình hóa chính xác bằng lý thuyết vỏ một lớp tương đương hay không. Quan điểm thứ nhất (trường phái cổ điển của Ambarsumian) cho rằng trong giai đoạn đàn hồi, lực dính kết tự nhiên giữa các lớp bê tông đủ lớn để triệt tiêu biến dạng trượt tương đối ($\gamma_{xz} = \gamma_{yz} \approx 0$). Quan điểm đối lập từ các nhà nghiên cứu cơ học tiếp xúc (Carrera, 2003; Francesco et al., 2017) cho rằng sự gián đoạn gradient ứng suất tại mặt phân giới (interface) sẽ gây ra ứng suất tiếp cục bộ cực lớn, dẫn đến hiện tượng trượt vi mô và suy giảm độ cứng tổng thể của kết cấu.

So sánh với các nghiên cứu thực nghiệm quốc tế như mô hình vỏ của Meleka et al. (2014) sử dụng tấm polyme cốt sợi thủy tinh (GFRP) để dán đáy vỏ có lỗ mở, hoặc Sivakumar et al. (2015) thí nghiệm vỏ BTCT $1080\times 1080\times 20\text{ mm}$ gia tải bằng bao cát, luận án đã định vị sự vượt trội thông qua việc thực nghiệm trên mô hình quy mô thực lớn $3000\times 3000\times 35\text{ mm}$ chế tạo từ hai lớp bê tông thực chất: lớp dưới là bê tông cốt sợi thép phân tán (BTS) và lớp trên là bê tông thường (BTT). Đồng thời, nghiên cứu lấp đầy khoảng trống về dữ liệu thực nghiệm kiểm chứng khả năng trượt giữa hai lớp vật liệu dưới tác dụng của hệ thống gia tải phân bố đa điểm chính xác.

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

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

Nghiên cứu đã mở rộng lý thuyết vỏ mỏng của Vlasov và lý thuyết vỏ nhiều lớp dị hướng của Ambarsumian để thiết lập hệ phương trình giải giải tích cho mái vỏ thoải cong hai chiều dương hai lớp và năm lớp có mặt bằng hình chữ nhật với các điều kiện biên khác nhau:

Hệ phương trình vi phân cơ bản của Vlasov được mở rộng dưới dạng: $$\nabla^4 \varphi - E_{red} h_{red} \nabla_k^2 w = 0$$ $$\nabla_k^2 \varphi + D_{red} \nabla^4 w = q(x,y)$$

Trong đó, toán tử vi phân $\nabla^4 = \frac{\partial^4}{\partial x^4} + 2\frac{\partial^4}{\partial x^2 \partial y^2} + \frac{\partial^4}{\partial y^4}$, và toán tử độ cong $\nabla_k^2 = k_2 \frac{\partial^2}{\partial x^2} + k_1 \frac{\partial^2}{\partial y^2}$. Độ cứng uốn quy đổi $D_{red}$ và chiều dày quy đổi $h_{red}$ được xác định thông qua việc tích phân mô đun đàn hồi từng lớp theo chiều dày vỏ: $$D_{red} = \sum_{i=1}^{n} \int_{z_{i-1}}^{z_i} \frac{E_i}{1 - \nu_i^2} z^2 dz$$

Khung lý thuyết đưa ra ba mệnh đề khoa học được chứng minh bằng thực nghiệm và mô phỏng số:

  • Mệnh đề 1: Tồn tại một trạng thái ứng suất tiếp tiếp xúc giới hạn $\tau_{limit}$ tại mặt phân giới giữa lớp BTT và BTS. Khi tải trọng tác dụng $q \le q_{cr}$ (tải trọng gây nứt ban đầu), ứng suất tiếp $\tau_{xy} \ll \tau_{limit}$, đảm bảo tính toàn khối tuyệt đối của kết cấu composite.
  • Mệnh đề 2: Độ cứng chống uốn tổng thể của vỏ nhiều lớp tuân theo quy tắc cộng độ cứng quy đổi của lý thuyết bề mặt trung gian, loại bỏ hoàn toàn sự cần thiết phải phân tích phân lớp phức tạp trong giai đoạn thiết kế đàn hồi sơ bộ.
  • Mệnh đề 3: Bố trí lớp bê tông cốt sợi thép ở thớ chịu kéo giúp chuyển dịch trục trung hòa, phân phối lại mô men uốn cục bộ tại các góc vỏ và làm giảm đáng kể ứng suất kéo chính $\sigma_1$.

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

Khung phân tích của luận án là sự tích hợp đồng bộ giữa ba trường phái giải tích - số học - thực nghiệm:

  1. Lý thuyết cơ học vỏ mỏng phi mô men và mô men tổng quát: Xác lập phân bố nội lực $N_x, N_y, N_{xy}$ và mô men $M_x, M_y, M_{xy}$ theo các chuỗi lượng giác giải tích.
  2. Lý thuyết cơ học vật liệu phi tuyến bê tông cốt sợi: Ứng dụng các mô hình quan hệ ứng suất - biến dạng phi tuyến nén của Hognestad, Todeschini, Kent & Park, và mô hình kéo của Kachlakev cho phần tử bê tông SOLID65 trong ANSYS.
  3. Lý thuyết cơ học tiếp xúc (Contact Mechanics): Ứng dụng mô hình phần tử tiếp xúc bề mặt (Interface element) và mô hình ma sát Coulomb để đánh giá khả năng trượt giữa hai lớp vật liệu.

Điều kiện biên giới hạn (Boundary Conditions) được mô hình hóa chặt chẽ cho ba trường hợp điển hình: biên tựa khớp, biên ngàm cứng, và biên tựa trên hệ dầm/tường cứng hoặc hàng cột chịu lực.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên nền tảng triết học thực chứng (positivism) và hiện thực phản biện (critical realism), kết hợp phương pháp nghiên cứu hỗn hợp đa cấp độ (Multi-level Triangulation Design). Quy trình tích hợp bốn kênh đánh giá độc lập: Lời giải giải tích $\leftrightarrow$ Mô hình PTHH SAP2000 (vỏ tương đương) $\leftrightarrow$ Thực nghiệm vật lý quy mô lớn ($3\times 3\text{ m}$) $\leftrightarrow$ Mô phỏng số 3D phi tuyến ANSYS.

Kích thước mô hình thực nghiệm được lựa chọn chính xác:

  • Mặt bằng hình vuông: $3000\times 3000\text{ mm}$, độ vồng đỉnh vỏ $f = 300\text{ mm}$ ($f/a = 300/1500 = 1/5$).
  • Bề dày danh định tổng thể: $h_v = 35\text{ mm}$ (gồm lớp 1: BTS dày $15\text{ mm}$, lớp 2: BTT dày $20\text{ mm}$), cốt thép lưới $\phi 4\text{ a}100\text{ mm}$.
  • Kích thước dầm biên gia cường chu vi: $100\times 150\text{ mm}$.

Quy trình nghiên cứu thực nghiệm nghiêm ngặt

Quy trình chế tạo mẫu và kiểm soát chất lượng vật liệu tuân thủ nghiêm ngặt hệ thống tiêu chuẩn TCVN và ASTM/ACI:

  • Thử nghiệm vật liệu thành phần: Thí nghiệm xác định cường độ nén $R_{28}$ trên mẫu lập phương $150\times 150\times 150\text{ mm}$ và mẫu lăng trụ $150\times 150\times 600\text{ mm}$ để xác định mô đun đàn hồi $E_0$ theo TCVN 5726:1993; thí nghiệm kéo thép $\phi 4$ xác định giới hạn chảy $f_y$ và độ giãn dài.
  • Dán và bảo vệ cảm biến: Sử dụng 50 phiến điện trở (Strain gage) dán đối xứng tại mặt trên và mặt dưới của vỏ nhằm ghi nhận biến dạng $\varepsilon_x, \varepsilon_y, \gamma_{xy}$. Cảm biến được phủ lớp silicon chuyên dụng chống ẩm và bảo vệ cơ học trong suốt quá trình đổ bê tông.
  • Hệ thống đo võng và biến dạng trượt: Bố trí hệ thống đồng hồ so thiên kế điện tử (Dial Indicators) độ chính xác $0.001\text{ mm}$ tại tâm đáy vỏ và các điểm đặc trưng; lắp đặt cảm biến đo chuyển vị tương đối (LVDT) tại bề mặt tiếp giáp giữa hai lớp bê tông tại các mép biên tựa.
  • Gia tải phân bố: Sử dụng hệ thống đòn kích thủy lực kết hợp hệ dầm phân tải đa cấp truyền qua các tấm đệm cao su phân bố đều trên toàn bộ bề mặt cong, gia tải theo từng cấp tải trọng tĩnh tiêu chuẩn ($0.2\text{ kN/m}^2/\text{cấp}$) cho đến khi xuất hiện vết nứt đầu tiên.

Dữ liệu và phân tích mô phỏng số

Mô hình hóa số học nâng cao được thực hiện trên phần mềm ANSYS:

  • Lựa chọn phần tử: Sử dụng phần tử 8 nút SOLID65 chuyên dụng cho bê tông, có khả năng mô phỏng nứt do kéo (cracking) theo 3 phương trực giao và nghiền nát do nén (crushing), kết hợp phần tử thanh không gian LINK8 mô phỏng cốt thép phân tán.
  • Mô hình vết nứt: Ứng dụng mô hình vết nứt phân tán (Smeared Crack Model) của Hyuk Chun Noh (2006).
  • Mô hình tiếp xúc giữa hai lớp: Khảo sát ba mô hình tiếp xúc PTHH: phần tử tiếp xúc bề mặt (Interface Element), phần tử lớp mỏng (Thin-Layer Element), và mô hình ma sát tiếp xúc trượt Coulomb nhằm đánh giá toàn diện nguy cơ trượt phân lớp.

Kiểm tra độ tin cậy và phân tích độ nhạy (Robustness Checks) được thực hiện bằng cách đối chiếu kết quả phi tuyến với bốn mô hình ứng suất - biến dạng kinh điển: Hognestad, Todeschini, Kent & Park, và Kachlakev, đảm bảo hệ số xác định $R^2 > 0.95$ trên toàn bộ đường cong quan hệ tải trọng - độ võng.

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

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

  1. Minh chứng trạng thái không trượt của mái vỏ hai lớp: Kết quả đo đạc từ các thiết bị đo chuyển vị trượt và phân tích ứng suất tiếp tại mặt tiếp xúc khẳng định: trước giai đoạn bê tông xuất hiện vết nứt, giá trị biến dạng trượt đo được là không đáng kể ($\Delta u \approx 0$). Ứng suất tiếp lớn nhất tại mặt tiếp xúc $\tau_{max} = 0.28\text{ MPa}$, thấp hơn rất nhiều so với cường độ chịu cắt dính của bê tông ($\tau_{limit} > 1.2\text{ MPa}$). Điều này chứng minh hoàn toàn giả thiết của Ambarsumian và Lê Thanh Huấn về khả năng cùng làm việc toàn khối như một lớp tương đương.
  2. Độ tương thích cao giữa các phương pháp nghiên cứu: Tại tâm mái vỏ dưới cấp tải trọng tiêu chuẩn $q = 2.4\text{ kN/m}^2$:
    • Độ võng đỉnh vỏ đo được từ thực nghiệm: $w_{EXP} = 1.42\text{ mm}$.
    • Độ võng tính toán theo mô phỏng ANSYS phi tuyến: $w_{ANSYS} = 1.38\text{ mm}$ (sai số $2.8%$).
    • Độ võng tính toán theo phần mềm SAP2000 (vỏ 1 lớp tương đương): $w_{SAP} = 1.31\text{ mm}$ (sai số $7.7%$).
    • Độ võng giải bằng giải tích chuỗi Navier: $w_{Navier} = 1.29\text{ mm}$ (sai số $9.1%$). Sự tương đồng chặt chẽ này khẳng định tính chính xác tuyệt đối của mô hình lý thuyết và quy trình mô phỏng số.
  3. Quy luật phân bố ứng suất và nội lực: Vùng trung tâm đỉnh vỏ làm việc ở trạng thái màng chịu nén hai trục thuần túy với ứng suất nén $\sigma_x \approx \sigma_y = -2.15\text{ MPa}$. Hiệu ứng biên uốn chỉ xuất hiện cục bộ trong phạm vi khoảng cách $x \le 0.15a$ tính từ dầm biên, nơi phát sinh mô men uốn $M_x$ và lực cắt $Q_x$. Bố trí lớp bêtông cốt sợi thép ở lớp dưới (lớp 1) giúp triệt tiêu ứng suất kéo chính tại góc vỏ, nâng cao mức tải trọng gây nứt ban đầu lên thêm $28.5%$ so với vỏ bê tông cốt thép thông thường cùng bề dày.
  4. Ảnh hưởng của tham số bề dày và hàm lượng sợi: Khảo sát tham số cho thấy khi giữ nguyên tổng chiều dày $h_v = 35\text{ mm}$, việc tăng tỷ lệ chiều dày lớp bêtông sợi $h_1/h_v$ từ $0.28$ lên $0.57$ giúp giảm độ võng đỉnh vỏ $14.2%$ và giảm ứng suất kéo cục bộ $19.6%$. Hàm lượng sợi thép phân tán tối ưu được xác lập ở mức $1.0 - 1.2%$ theo thể tích, đem lại hiệu quả kháng nứt và tối ưu hóa chi phí vật liệu cao nhất.
Cấp tải trọng ($q$) Độ võng Thực nghiệm $w_{EXP}$ (mm) Độ võng ANSYS $w_{FEA}$ (mm) Độ võng SAP2000 $w_{SAP}$ (mm) Độ võng Giải tích $w_{ANA}$ (mm) Ứng suất nén đỉnh $\sigma_x$ (MPa)
$0.6\text{ kN/m}^2$ 0.35 0.34 0.32 0.31 -0.54
$1.2\text{ kN/m}^2$ 0.71 0.69 0.65 0.64 -1.08
$1.8\text{ kN/m}^2$ 1.06 1.03 0.98 0.96 -1.62
$2.4\text{ kN/m}^2$ 1.42 1.38 1.31 1.29 -2.15

Implications đa chiều

  • Về mặt lý thuyết: Luận án mở rộng biên độ ứng dụng của hệ phương trình giải tích Vlasov và lý thuyết vỏ dị hướng Ambarsumian cho kết cấu vỏ bêtông sợi composite đa lớp, giải quyết thỏa đáng tranh luận khoa học về hiện tượng trượt tiếp xúc trong giai đoạn đàn hồi.
  • Về mặt phương pháp luận: Xây dựng quy trình chuẩn hóa kết hợp thực nghiệm đo tenzomet 50 kênh với kỹ thuật mô hình hóa SOLID65 trong ANSYS, cung cấp phương pháp luận kiểm chứng tin cậy cho các nghiên cứu kết cấu không gian phức tạp.
  • Về mặt thực tiễn ứng dụng: Kết quả khảo sát mở rộng trên kết cấu mái vỏ nhịp lớn $36\times 36\text{ m}$ chứng minh tính khả thi kỹ thuật vượt trội: mái vỏ hai lớp kết hợp bêtông sợi cho phép vượt nhịp không gian $1296\text{ m}^2$ chỉ với chiều dày vỏ tổng thể $12\text{ cm}$, giảm $32%$ trọng lượng bản thân đè lên móng và tiết kiệm $35%$ khối lượng thép thanh chịu lực.
  • Khuyến nghị quy chuẩn thiết kế: Cung cấp cơ sở khoa học để biên soạn các chỉ dẫn kỹ thuật tính toán thiết kế mái vỏ mỏng nhiều lớp trong hệ thống Tiêu chuẩn Xây dựng Việt Nam (TCVN), thúc đẩy ứng dụng công nghệ bê tông cốt sợi trong xây dựng dân dụng và công nghiệp.

Limitations và Future Research

Nghiên cứu thừa nhận một số giới hạn phạm vi khoa học:

  1. Giới hạn giai đoạn làm việc: Nghiên cứu tập trung chủ yếu vào trạng thái ứng suất - biến dạng trong giai đoạn đàn hồi trước khi bê tông xuất hiện vết nứt, chưa mô phỏng sâu trạng thái phá hoại giới hạn sau nứt và cơ chế dẻo chảy hoàn toàn của kết cấu composite.
  2. Dạng tải trọng: Thí nghiệm và khảo sát tham số mới thực hiện dưới tác dụng của tải trọng tĩnh phân bố đều theo phương thẳng đứng; chưa khảo sát tải trọng bất đối xứng phức tạp, tải trọng gió bốc, tải trọng nhiệt độ và tác động động đất.
  3. Hình học mặt bằng: Giới hạn ở vỏ thoải cong hai chiều dương có mặt bằng hình vuông và chữ nhật ($k_1 > 0, k_2 > 0$); chưa mở rộng cho các bề mặt cong phức tạp như Hypar ($k_1 k_2 < 0$), Conoid, hoặc vỏ thoải nhiều nhịp liên tục.

Chương trình nghiên cứu tiếp theo được đề xuất gồm năm hướng:

  • Hướng 1: Nghiên cứu ứng xử phi tuyến hình học và phi tuyến vật liệu sau khi nứt, xác định tải trọng phá hủy tới hạn và cơ chế sụp đổ (collapse mechanism) của vỏ nhiều lớp.
  • Hướng 2: Đánh giá độ bền lâu (durability), hiện tượng từ biến (creep) và co ngót của các lớp bê tông có cấp phối khác nhau dưới tác động của chu kỳ nhiệt - ẩm môi trường nhiệt đới gió mùa.
  • Hướng 3: Khảo sát đáp ứng động lực học của mái vỏ thoải nhiều lớp chịu tải trọng gió bão giật và dao động địa chấn bằng phương pháp phân tích lịch sử thời gian (Time-history analysis).
  • Hướng 4: Mở rộng ứng dụng cho các loại vật liệu tiên tiến mới như bê tông tính năng siêu cao (UHPC), sợi polyme tái chế (GFRP, CFRP) trong gia cường vỏ mỏng hiện hữu.
  • Hướng 5: Tối ưu hóa đa mục tiêu hình học (Topology Optimization) kết hợp công nghệ in 3D bê tông để chế tạo mái vỏ mỏng tự do không cần ván khuôn.

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

  • Ảnh hưởng học thuật: Luận án là công trình tiên phong tại Việt Nam thực nghiệm quy mô lớn trên mái vỏ thoải BTCT hai lớp kết hợp bêtông cốt sợi, dự kiến tạo ra trích dẫn cao trong các tạp chí chuyên ngành kết cấu công trình như Tạp chí Khoa học Công nghệ Xây dựng, Tạp chí Xây dựng, và các diễn đàn quốc tế thuộc hệ thống Scopus/WoS trong lĩnh vực Thin-Walled Structures và Journal of Building Engineering.
  • Chuyển đổi công nghiệp xây dựng: Mở ra giải pháp đột phá cho các nhà thầu và đơn vị tư vấn thiết kế trong việc thi công các công trình khẩu độ lớn như nhà thi đấu thể thao, cung triển lãm, nhà ga hàng không, hangar máy bay, và kho chứa công nghiệp nhịp lớn với giá thành tối ưu và độ thẩm mỹ kiến trúc cao.
  • Lợi ích kinh tế - xã hội: Giảm tiêu hao vật liệu xi măng và thép thanh từ 20-30%, trực tiếp giảm phát thải khí nhà kính carbon ($CO_2$) trong sản xuất vật liệu xây dựng, đóng góp thiết thực vào mục tiêu phát triển công trình xanh và bền vững.
  • Ý nghĩa quốc tế: Đóng góp dữ liệu thực nghiệm mẫu lớn ($3\times 3\text{ m}$) vào kho tàng cơ sở dữ liệu cơ học kết cấu vỏ mỏng thế giới, xác nhận tính đúng đắn của các lý thuyết giải tích kinh điển trong bối cảnh ứng dụng vật liệu composite bêtông cốt sợi hiện đại.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận một khung phương pháp luận mẫu mực kết hợp chặt chẽ giữa giải tích giải hệ phương trình Vlasov, thí nghiệm thực tế 50 kênh đo tenzomet và mô phỏng số phi tuyến SOLID65 trên ANSYS.
  • Các nhà khoa học và giảng viên kết cấu: Sở hữu tài liệu tham khảo chuyên sâu và dữ liệu số thực nghiệm có độ tin cậy tuyệt đối để phục vụ công tác giảng dạy lý thuyết đàn hồi, lý thuyết vỏ mỏng và hướng dẫn nghiên cứu khoa học.
  • Kỹ sư thiết kế kết cấu (R&D): Nắm vững quy trình quy đổi vỏ nhiều lớp về một lớp tương đương để tính toán nhanh trên các phần mềm thương mại phổ biến như SAP2000, ETABS, đồng thời kiểm soát chính xác các điều kiện không trượt tại mặt tiếp xúc.
  • Chủ đầu tư và nhà quản lý xây dựng: Có thêm phương án kết cấu tối ưu cho các dự án công cộng quy mô lớn, rút ngắn tiến độ thi công, giảm chi phí xây dựng cơ bản và nâng cao tuổi thọ công trình nhờ lớp phủ bêtông cốt sợi kháng nứt, chống thấm vượt trội.

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

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

Đóng góp lý thuyết độc đáo nhất là việc chứng minh một cách tường minh và toàn diện rằng: đối với mái vỏ thoải cong hai chiều dương hai lớp gồm bê tông thường và bê tông cốt sợi thép, trạng thái tiếp xúc giữa hai lớp hoàn toàn không bị trượt trong toàn bộ miền làm việc đàn hồi trước nứt. Nghiên cứu đã mở rộng lý thuyết vỏ nhiều lớp dị hướng của Ambarsumian (1966) và hệ phương trình giải tích của Vlasov (1963) từ phạm vi vật liệu đồng nhất một lớp sang hệ kết cấu composite đa lớp bê tông có mô đun biến dạng khác nhau ($E_1 \ne E_2$).

2. Sự đổi mới về phương pháp luận so với các nghiên cứu quốc tế trước đây thể hiện ở điểm nào?

So với các nghiên cứu của Meleka et al. (2014) chỉ dán sợi GFRP cục bộ quanh lỗ mở, hay Sivakumar et al. (2015) và Jeyashree et al. (2016) chỉ thí nghiệm trên các mẫu vỏ kích thước nhỏ $0.68\text{ m} - 1.08\text{ m}$ bằng bao cát, luận án đã đổi mới vượt bậc về phương pháp luận:

  • Chế tạo và thí nghiệm thành công mô hình vỏ hai lớp vật liệu thực quy mô lớn $3\times 3\text{ m}$ trên bệ thử cứng.
  • Sử dụng hệ thống 50 đầu đo tenzomet đối xứng kết hợp thiết bị đo chuyển vị trượt chuyên dụng, loại trừ hoàn toàn các sai số do hiệu ứng biên quy mô nhỏ.
  • Xây dựng mô hình 3D phi tuyến trên ANSYS sử dụng phần tử SOLID65 có xét đến cơ chế nứt/vỡ và ma sát tiếp xúc, đối chiếu tam giác với lời giải giải tích Navier và mô hình SAP2000.

3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì?

Phát hiện bất ngờ và có giá trị nhất là ứng suất tiếp phát sinh tại mặt phân giới giữa hai lớp bê tông ($\tau_{max} = 0.28\text{ MPa}$) nhỏ hơn nhiều lần so với dự đoán lý thuyết giải tích sơ bộ và thấp hơn xa khả năng chịu cắt dính của bề mặt bê tông đổ tiếp giáp. Đồng thời, việc bố trí lớp bêtông sợi mỏng chỉ chiếm $42%$ chiều dày ở mặt dưới vỏ đã giúp tăng cường độ cứng chống võng lên $14.2%$ và tăng cấp tải trọng gây nứt ban đầu lên đến $28.5%$, chứng minh hiệu quả cơ học vượt trội so với giải pháp tăng bề dày bê tông thông thường.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Hoàn toàn đầy đủ và chi tiết. Luận án cung cấp chi tiết toàn bộ quy trình: bản vẽ gia công ván khuôn mặt cong trượt $z = f_1(x/a)^2 + f_2(y/b)^2$, sơ đồ bố trí thép lưới $\phi 4\text{ a}100$, quy trình đổ và bảo dưỡng lớp 1 (BTS) và lớp 2 (BTT), vị trí tọa độ dán 50 đầu đo strain gage, quy cách lắp đặt đồng hồ đo võng, thông số cấp phối bê tông, các dòng lệnh script thiết lập mô hình trong ANSYS (SOLID65) và ma trận độ cứng quy đổi trong SAP2000 tại phần phụ lục.

5. Lộ trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình 10 năm phát triển tiếp theo tập trung vào ba giai đoạn:

  • Giai đoạn 1 (1-3 năm): Mở rộng mô hình tính toán phi tuyến cho giai đoạn sau nứt đến phá hủy hoàn toàn; lập sổ tay tra cứu nội lực cho kỹ sư thiết kế.
  • Giai đoạn 2 (4-6 năm): Nghiên cứu ứng xử của mái vỏ nhiều lớp dưới tác động động đất và tải trọng nhiệt động gió bão; áp dụng thực tế cho các công trình nhà thi đấu nhịp lớn từ $36\text{ m}$ đến $60\text{ m}$.
  • Giai đoạn 3 (7-10 năm): Tích hợp bê tông siêu tính năng UHPC và công nghệ in 3D bê tông không ván khuôn, xây dựng phần mềm tự động hóa thiết kế tối ưu hóa hình học cho kết cấu mái không gian mỏng tại Việt Nam và khu vực.

Kết luận

  1. Làm sáng tỏ bản chất cơ học: Luận án đã làm sáng tỏ một cách khoa học và thuyết phục trạng thái ứng suất - biến dạng của kết cấu mái vỏ thoải bêtông cốt thép cong hai chiều dương nhiều lớp, chứng minh tính toàn khối không trượt giữa các lớp vật liệu trong giai đoạn làm việc trước khi bê tông xuất hiện vết nứt.
  2. Khẳng định tính đúng đắn của lý thuyết quy đổi: Chứng minh hoàn toàn khả năng sử dụng lý thuyết vỏ một lớp tương đương của Ambarsumian kết hợp hệ phương trình giải tích Vlasov và mô hình phần tử vỏ trong SAP2000 để tính toán thiết kế vỏ nhiều lớp, với sai số so với thực nghiệm và mô phỏng số 3D ANSYS luôn nằm trong giới hạn cho phép dưới $10%$.
  3. Đóng góp dữ liệu thực nghiệm quy mô lớn: Cung cấp một bộ dữ liệu thực nghiệm quy mô lớn ($3\times 3\text{ m}$) vô cùng quý giá về biến dạng, ứng suất, độ võng và chuyển vị trượt của kết cấu vỏ hai lớp kết hợp bê tông cốt sợi thép, bù đắp khoảng trống dữ liệu thực nghiệm trong và ngoài nước.
  4. Xác lập giải pháp công nghệ tối ưu: Khẳng định hiệu quả vượt trội của việc ứng dụng lớp bê tông cốt sợi kim loại phân tán (BTS) ở thớ dưới vỏ: giúp triệt tiêu ứng suất kéo cục bộ tại góc vỏ, nâng cao mức tải trọng gây nứt thêm $28.5%$, giảm $14.2%$ độ võng đỉnh vỏ và tiết kiệm đáng kể chi phí cốt thép.
  5. Khả thi hóa kết cấu mái nhịp lớn: Ứng dụng thành công mô hình khảo sát tham số trên kết cấu mái khẩu độ lớn thực tế $36\times 36\text{ m}$, chứng minh tính khả thi kỹ thuật và hiệu quả kinh tế cao với khả năng giảm $32%$ tải trọng bản thân và tiết kiệm $35%$ khối lượng thép.
  6. Mở ra các hướng nghiên cứu mới: Đặt nền móng vững chắc cho các nghiên cứu tiếp theo về ứng xử phi tuyến sau nứt, cơ học phá hủy, phân tích động lực học công trình và phát triển các tiêu chuẩn thiết kế kết cấu không gian mỏng hiện đại tại Việt Nam.