Tổng quan về luận án

Sự phát triển nhanh chóng của các đô thị hiện đại đặt ra yêu cầu cấp thiết đối với việc xây dựng các công trình nhà cao tầng từ 20 đến 40 tầng. Trong hệ thống kết cấu chịu lực của các tòa nhà này, vách bê tông cốt thép (BTCT) đóng vai trò là cấu kiện kháng chấn và chịu tải trọng ngang chủ đạo do gió bão và động đất gây ra. Tuy nhiên, các phương pháp thiết kế truyền thống theo tiêu chuẩn xây dựng trước đây chủ yếu dựa trên quan điểm độ bền và độ cứng trong miền đàn hồi tuyến tính, chưa phản ánh đầy đủ khả năng làm việc thực tế của kết cấu khi bước vào giai đoạn phi tuyến dẻo. Xu hướng hiện đại trên thế giới đã chuyển dịch mạnh mẽ sang phương pháp thiết kế dựa trên chuyển vị (Displacement-Based Design - DBD), đòi hỏi sự hiểu biết sâu sắc về ứng xử phi tuyến, cơ chế tiêu tán năng lượng và quá trình suy giảm độ cứng của vách BTCT khi chịu tải trọng lặp đảo chiều chu kỳ.

Khoảng trống nghiên cứu (Research Gap) then chốt hiện nay nằm ở sự thiếu hụt các mô hình phần tử hữu hạn phi tuyến tích hợp toàn diện vật liệu hỗn hợp bê tông và cốt thép có khả năng mô phỏng đồng thời biến dạng uốn, biến dạng trượt cắt và hiệu ứng trễ (hysteretic behavior). Các phần mềm phân tích kết cấu thương mại phổ biến như ETABS, SAP2000, MIDAS hay Robot Structural Analysis vẫn tồn tại những hạn chế đáng kể trong việc mô hình hóa chính xác tính phi tuyến vật liệu hỗn hợp cho phần tử vách tấm phẳng chịu tải trọng đổi chiều phức tạp, hoặc đòi hỏi chi phí tính toán vi mô quá lớn. Luận án tiến sĩ kỹ thuật chuyên ngành Kỹ thuật Xây dựng (Mã số: 9580201) của nghiên cứu sinh Giang Văn Khiêm, dưới sự hướng dẫn khoa học của PGS. TS. Nghiêm Mạnh Hiến và GS. Nguyễn Tiến Chương tại Viện Khoa học Công nghệ Xây dựng (IBST), đã giải quyết trực tiếp khoảng trống này thông qua việc kết hợp chặt chẽ giữa phát triển mô hình giải tích lý thuyết, thuật toán số và nghiên cứu thực nghiệm quy mô lớn.

Mục tiêu và câu hỏi nghiên cứu của luận án được cụ thể hóa qua hệ thống giả thuyết khoa học:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để xây dựng một phần tử hữu hạn hỗn hợp tấm phẳng phản ánh chính xác trạng thái ứng suất - biến dạng phi tuyến hai chiều và ba chiều của bê tông kết hợp cốt thép chịu tải trọng chu kỳ?
  • Câu hỏi nghiên cứu 2 (RQ2): Ứng xử cơ học thực tế của vách BTCT (độ cứng, góc xoay chân vách, độ dẻo, sự phân bố độ cong, biến dạng cắt và hệ số cản nhớt tương đương) biến đổi như thế nào dưới tác dụng đồng thời của lực nén dọc trục không đổi và tải trọng ngang lặp đảo chiều?
  • Câu hỏi nghiên cứu 3 (RQ3): Làm thế nào để thiết lập một chương trình tính toán số độc lập có độ tin cậy cao, giải quyết hội tụ phi tuyến động lực học và tĩnh học phi tuyến cho vách BTCT?

Hệ thống giả thuyết nghiên cứu bao gồm:

  • Giả thuyết 1 (H1): Việc tích hợp mô hình vết nứt phân tán (smeared crack model) cùng trường ứng suất tương đương một trục vào phần tử tấm tứ giác đẳng tham số sẽ mô phỏng chính xác sự suy giảm độ cứng và hiệu ứng kẹp khe nứt (pinching effect).
  • Giả thuyết 2 (H2): Vùng biên chịu mô men ($B_t, B_p = 0,2l_w$) đóng vai trò quyết định khả năng tiêu tán năng lượng dẻo; sự suy giảm lực dính và biến dạng trượt cắt có mối quan hệ tương tác phi tuyến mật thiết với tải trọng dọc trục tác dụng.
  • Giả thuyết 3 (H3): Chương trình chuyên dụng CSW tự phát triển có khả năng dự báo đường cong bao trễ và các điểm chảy dẻo với sai số dưới 10% so với kết quả thực nghiệm trong nước và quốc tế.

Khung lý thuyết của nghiên cứu được xây dựng trên nền tảng Cơ học phá hủy bê tông, Lý thuyết trường nén cải tiến (Modified Compression Field Theory - MCFT) của Vecchio và Collins, kết hợp các mô hình vật liệu kinh điển của Kent - Park, Mander và Menegotto - Pinto. Luận án giới hạn phạm vi nghiên cứu vào vách phẳng BTCT toàn khối trong nhà cao tầng phù hợp điều kiện chế tạo và tiêu chuẩn Việt Nam TCVN 5574:2018, chịu tổ hợp lực nén đúng tâm và tải trọng ngang đảo chiều trong mặt phẳng vách.

Literature Review và Positioning

Lịch sử mô hình hóa ứng xử phi tuyến của vách BTCT trải qua bốn trường phái tiếp cận giải tích chính:

  1. Mô hình dầm tương đương (Equivalent Beam Model / Beam-Column Element): Khởi xướng bởi Takayanagi và Schnobrich (1976), Soleimani và cộng sự (1979), Otani (1985) và Martinelli - Filippou (2009). Mô hình này tập trung biến dạng dẻo vào các lò xo xoay phi tuyến ở đầu dầm. Hạn chế cốt tử là bỏ qua sự dịch chuyển của trục trung hòa và không mô phỏng được tương tác uốn - cắt đa chiều.
  2. Mô hình giàn tương đương (Equivalent Truss Model): Phát triển bởi Hiraishi (1983) và Panagiotou cùng cộng sự (2012). Cấu kiện được rời rạc hóa thành chuỗi thanh giàn chịu kéo - nén. Phương pháp này phụ thuộc quá nhiều vào quy luật phân bố vết nứt giả định ban đầu.
  3. Mô hình thớ và mô hình vĩ mô nhiều đường thẳng đứng (Macro Fiber Models / MVLM): Tiêu biểu là mô hình Ba Đường Thẳng đứng (TVL) của Kabeyasawa và cộng sự (1982), mô hình Nhiều Đường Thẳng đứng (MVLM) của Vulcano, Bertero và Colotti (1988), mô hình dàn mềm xoay (RA-STM) của Pang và Hsu (1995), và các cải tiến của Massone - Wallace (2004, 2006). Mô hình thớ phân chia tiết diện thành các dải vật liệu đơn trục, phản ánh tốt tương tác lực dọc - uốn nhưng gặp nhiều trở ngại khi đánh giá biến dạng trượt cắt trong vách thấp và vách có tỷ số chiều cao trên chiều dài biến thiên.
  4. Mô hình vi mô phần tử hữu hạn (Micro Continuum FE Models): Khởi xướng từ các nghiên cứu trực hướng của Darwin và Pecknold (1974, 1977), mô hình vết nứt phân tán của Bazant và Oh (1985), Lý thuyết vùng nén CFT của Collins (1978), Vecchio và Collins (1982, 1986), và mô hình màng vết nứt CMM của Kaufmann (1998), Foster và Marti (2003).

Tranh luận học thuật cốt lõi trong tài liệu quốc tế xoay quanh hai luồng quan điểm đối lập:

  • Quan điểm thứ nhất (Mô hình vết nứt rời rạc - Discrete Crack Model): Cho rằng cần mô tả tường minh từng khe nứt và sự trượt cốt thép dọc theo bề mặt tiếp xúc bê tông - cốt thép. Quan điểm này có độ chính xác cục bộ cao nhưng đòi hỏi tài nguyên tính toán khổng lồ và thường gặp bất ổn định hội tụ số khi mật độ vết nứt dày đặc.
  • Quan điểm thứ hai (Mô hình vết nứt phân tán - Smeared Crack Model): Giả định vết nứt phân bố liên tục trên toàn bộ thể tích tích phân của phần tử, điều chỉnh ma trận cấu tạo đàn dẻo thông qua các hệ số suy giảm mô đun đàn hồi $\beta_1$, suy giảm cường độ chịu kéo $\beta_2$ và suy giảm cường độ chịu nén $\beta_3$. Luận án định vị đi theo hướng tiếp cận vết nứt phân tán cải tiến, dung hòa được độ chính xác cơ học và hiệu năng tính toán.

Để xác lập vị thế khoa học, luận án đã so sánh đối chuẩn với ít nhất hai công trình thực nghiệm kinh điển của thế giới:

  • Thí nghiệm của Lefas, Kotsovos và Ambraseys (1990): Thực hiện trên 13 mẫu vách có tỷ số chiều cao/chiều dài $h_w/l_w = 1,0$ (Vách loại I: $0,75 \times 0,75 \times 0,07\text{ m}$) và $h_w/l_w = 2,0$ (Vách loại II: $1,30 \times 0,65 \times 0,065\text{ m}$), chịu mức nén dọc trục $n = 0; 0,1; 0,2 f_c' A_g$. Kết quả của Lefas chỉ ra lực nén dọc trục làm tăng 25% đến 30% khả năng chịu tải ngang nhưng làm giảm chuyển vị cực hạn tại thời điểm phá hoại.
  • Thí nghiệm của Thomsen và Wallace (1995, 2004): Tiến hành trên các vách phẳng RW1, RW2 ($1,219 \times 0,102 \times 3,658\text{ m}$) và vách chữ T (TW1, TW2) tỉ lệ 1/4 theo tiêu chuẩn UBC dưới các chu kỳ biến dạng lệch tầng (drift) từ 0,1% đến 3,0%. Thí nghiệm khẳng định biến dạng uốn chi phối mạnh mẽ cơ chế làm việc của vách thanh mảnh, đồng thời vùng biên chịu mô men tập trung ứng suất dẻo rất cao.

Luận án của NCS. Giang Văn Khiêm định vị sự đóng góp của mình bằng việc mở rộng các mô hình 2D kinh điển sang không gian phần tử hỗn hợp đa nút, tích hợp đồng thời cơ chế trễ lặp của cả bê tông và cốt thép, được kiểm chứng trực tiếp bằng chuỗi thí nghiệm vách tỷ lệ lớn chế tạo theo quy chuẩn xây dựng 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 mang lại những bước tiến lý thuyết quan trọng thông qua việc mở rộng và kết hợp các học thuyết cơ học kết cấu:

  • Kế thừa và mở rộng mô hình ứng xử nén của Kent và Park (1973), Mander và cộng sự (1988) cho bê tông bị kiềm chế nở hông trong vùng biên, biểu diễn định lượng qua hệ số tăng cường độ $K$ và độ dốc nhánh giảm $Z$: $$\sigma_c = K f_c' \left[ \frac{2\varepsilon_c}{\varepsilon_0 K} - \left(\frac{\varepsilon_c}{\varepsilon_0 K}\right)^2 \right]$$ với biến dạng tương ứng ứng suất đỉnh $\varepsilon_0 = 0,002 K$, trong đó hệ số $K = 1 + \frac{\rho_s f_{yh}}{f_c'}$.
  • Bổ sung quy tắc suy giảm trễ hysteretic của bê tông chịu kéo - nén lặp theo mô hình Chang - Mander đơn giản hóa và đề xuất của Waugh (2009), giải quyết chính xác hiện tượng đóng/mở khe nứt (crack closing/reopening) khi đổi chiều tải trọng thông qua ứng suất ngưỡng khép vết nứt $\sigma_f$.
  • Mở rộng mô hình cốt thép phi tuyến Menegotto - Pinto (1973) và mô hình Kent - Park (1973) có xét đến hiệu ứng Bauschinger, phản ánh trạng thái làm việc thực tế của thanh thép đặt chìm trong môi trường bê tông so với thanh thép thử nghiệm độc lập trong không khí.

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

Khung phân tích của luận án tích hợp đồng thời ba lý thuyết nền tảng:

  1. Lý thuyết dẻo dòng chảy (Plastic Flow Theory) dựa trên hàm thế dẻo $g$ và các bất biến ứng suất $I_1, J_2$ để xác định mặt chảy dẻo vật liệu.
  2. Lý thuyết trường nén cải tiến (MCFT) xác định góc nứt trực giao theo phương ứng suất chính.
  3. Phương pháp phần tử hữu hạn phi tuyến kết cấu tấm phẳng.

Cốt lõi của khung phân tích là Phần tử hỗn hợp bê tông cốt thép (Reinforced Concrete Mixed Element). Trong phần tử này, trường chuyển vị của bê tông được nội suy qua các hàm dạng của phần tử tấm tứ giác đẳng tham số (4 nút hoặc 8 nút), trong khi cốt thép được mô hình hóa dưới dạng lớp tương đương hoặc phần tử thanh nhúng trực tiếp, liên kết hoàn toàn với các nút của phần tử bê tông: $${d\varepsilon} = [B]{du}$$ $${d\sigma} = [D_{ep}]{d\varepsilon}$$ Ma trận độ cứng tiếp tuyến của phần tử hỗn hợp $[K^e]$ được thiết lập bằng cách tích hợp trên toàn miền thể tích tứ giác sử dụng phép cầu phương Gauss: $$[K^e] = \int_{V} [B]^T [D_{ep}] [B] , dV$$ trong đó $[D_{ep}]$ là ma trận đàn dẻo dẫn xuất từ trạng thái ứng suất - biến dạng tức thời, tự động cập nhật sự suy giảm độ cứng theo cả hai phương trực giao.

Điều kiện biên và phạm vi áp dụng của khung phân tích được xác định rõ: phần tử áp dụng cho trạng thái ứng suất phẳng trong mặt phẳng vách, giả thiết lực bám dính giữa bê tông và cốt thép là hoàn hảo (chưa xét đến hiện tượng trượt tuột cục bộ ở cấp độ vi mô hạt cốt liệu), và góc nghiêng vết nứt quay đồng phương với trục chính của ten-sơ biến dạng.

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 khoa học thực chứng (Positivism) kết hợp chủ nghĩa hiện thực phê phán (Critical Realism). Thiết kế nghiên cứu đa cấp độ (Multi-level Experimental & Numerical Design) được triển khai bài bản qua 3 giai đoạn:

  • Cấp độ 1 (Vật liệu): Thí nghiệm xác định đường cong quan hệ ứng suất - biến dạng chuẩn của mẫu nén trụ bê tông ($\phi 150 \times 300\text{ mm}$), mẫu lăng trụ, và kéo đứt các thanh cốt thép dọc, cốt thép đai.
  • Cấp độ 2 (Cấu kiện thực nghiệm): Thí nghiệm 02 mẫu vách BTCT quy mô lớn (ký hiệu V1 và V2) chịu tải trọng nén dọc trục và tải trọng ngang lặp chu kỳ tại Phòng thí nghiệm công trình - Trung tâm Khoa học Công nghệ GTVT (Trường Đại học Giao thông vận tải) và Phòng thí nghiệm Đại học Xây dựng.
  • Cấp độ 3 (Mô phỏng số học): Xây dựng thuật toán, lập trình mã nguồn chương trình CSW và đối chuẩn chéo (cross-validation) với thực nghiệm của luận án và thực nghiệm quốc tế của Thomsen - Wallace.

Quy trình nghiên cứu rigorous

Thiết kế và chế tạo mẫu thí nghiệm tuân thủ chặt chẽ các quy định của tiêu chuẩn Việt Nam TCVN 5574:2018 và hướng dẫn thí nghiệm kháng chấn ACI:

  • Thông số hình học mẫu: Chiều cao vách $h_w = 2300\text{ mm}$, chiều dài tiết diện $l_w = 1200\text{ mm}$, chiều dày $h = 100\text{ mm}$ (tỷ số thanh mảnh $h_w/l_w = 1,92$). Phần đế móng bê tông cốt thép kích thước lớn ($1800 \times 600 \times 500\text{ mm}$) được neo chặt xuống sàn phản lực bằng hệ thống bu lông cường độ cao để chống trượt và chống lật tuyệt đối.
  • Đặc trưng vật liệu: Bê tông sử dụng cấp độ bền B20 với cấp phối mác chuẩn ($385\text{ kg}$ xi măng PC40, $695\text{ kg}$ cát vàng, $1180\text{ kg}$ đá dăm $1 \times 2$, $195\text{ lít}$ nước trên $1\text{ m}^3$). Cường độ chịu nén mẫu lập phương thực tế đạt trung bình $R_{cube} = 26,4\text{ MPa}$, mô đun đàn hồi $E_c = 2,75 \times 10^4\text{ MPa}$. Cốt thép dọc sử dụng thép thanh vằn đường kính $12\text{ mm}$ ($f_y = 365\text{ MPa}, f_u = 520\text{ MPa}$), cốt đai và cốt phân bố sử dụng thép trơn đường kính $6\text{ mm}$ ($f_y = 295\text{ MPa}$).
  • Hệ thống thiết bị đo và thu thập dữ liệu tự động:
    • 14 cảm biến đo chuyển vị biến trở tuyến tính (LVDT) bố trí đa tầng: đo chuyển vị đỉnh vách, độ co dãn biên $\Delta_1, \Delta_2$ để tính góc xoay chân vách $\theta$, đo biến dạng trượt cắt chéo $\gamma$, đo chuyển vị ngoài mặt phẳng để kiểm soát mất ổn định uốn - xoắn, và đo chuyển vị trượt/lật đế móng.
    • Hệ thống đầu đo điện trở (strain gauges) dán trực tiếp lên các thanh thép dọc vùng biên chịu kéo/nén ở các cao độ khác nhau và dán trên bề mặt bê tông.
    • Thiết bị quang học đo bề rộng và sự phát triển của vết nứt.
    • Quy trình gia tải lặp đảo chiều điều khiển bằng chuyển vị (Displacement Control Protocol) với các chu kỳ tăng dần theo các bước chuyển vị đỉnh $\Delta = 2, 4, 6, 8, 10, 15, 20, 25, 30, 35, 40\text{ mm}$, mỗi cấp lặp 2 chu kỳ nhằm đánh giá chính xác sự suy giảm độ cứng và khả năng tiêu tán năng lượng trễ.
                  SƠ ĐỒ BỐ TRÍ THÍ NGHIỆM VÁCH BTCT
               +-----------------------------------+
               |  Kích tạo lực nén đứng N (P_ax)   |
               +-----------------------------------+
                                 |
                                 v
        [ Kích ngang P ] <===============> [ Dầm đỉnh truyền lực ]
        (Lặp đảo chiều)  |               |
                         |               |
                         |               |
                         |  VÁCH BTCT    |  <-- Hệ LVDT đo chuyển vị
                         | (hw = 2300mm  |      và biến dạng trượt
                         |  lw = 1200mm  |
                         |  h  = 100mm)  |
                         |               |
                         +---------------+
                         |   ĐẾ MÓNG     |  <-- Neo bu lông sàn phản lực
          ===============+===============+=================
                          SÀN PHÒNG THÍ NGHIỆM

Data và phân tích

Để giải hệ phương trình cân bằng phi tuyến tĩnh và động lực học: $$[M]{\ddot{u}} + [C]{\dot{u}} + [K(u)]{u} = {P(t)}$$ Luận án ứng dụng phương pháp lặp Newmark kết hợp thuật toán Newton - Raphson cải tiến (Modified Newton-Raphson Method). Chương trình CSW (Concrete Shear Wall) được xây dựng hoàn chỉnh, cho phép người dùng nhập dữ liệu hình học, tải trọng, đường cong ứng suất - biến dạng thực nghiệm của vật liệu, tự động rời rạc hóa lưới phần tử và xuất kết quả trường ứng suất $\sigma_x, \sigma_y, \tau_{xy}$, độ cong $\phi$, chuyển vị đỉnh và vết nứt.

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

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

Nghiên cứu mang lại 5 phát hiện đột phá có ý nghĩa khoa học và thực tiễn sâu sắc:

  1. Quy luật suy giảm độ cứng chu kỳ và hiệu ứng kẹp vết nứt (Pinching Effect): Đường trễ quan hệ giữa tải trọng ngang $P$ và chuyển vị đỉnh $\Delta$ biểu hiện rõ nét hiệu ứng thắt hẹp ở lân cận gốc tọa độ khi chuyển vị vượt quá mức drift $1,0%$. Độ cứng cát tuyến (secant stiffness) của vách giảm mạnh từ $40%$ đến $65%$ sau khi cốt thép dọc vùng biên đạt giới hạn chảy dẻo ($f_y = 365\text{ MPa}$) và bê tông bắt đầu nứt tách.
  2. Cơ chế tập trung biến dạng dẻo tại chân vách: Vùng phá hoại dẻo tập trung chủ yếu trong phạm vi chiều cao $h_{cr} \approx 0,2 - 0,25 h_w$ tính từ mặt đế móng. Tại đây, cốt thép dọc vùng biên chịu biến dạng kéo vượt quá $1,5%$ trước khi đảo chiều chịu nén dẫn đến hiện tượng vỡ vụn bê tông bảo vệ mép ngoài.
  3. Mối tương tác phi tuyến giữa uốn và biến dạng trượt cắt: Biến dạng trượt cắt $\gamma$ chiếm từ $15%$ đến $22%$ tổng chuyển vị ngang đỉnh trong giai đoạn đàn hồi ban đầu, nhưng tăng vọt lên chiếm $32% - 38%$ khi các vết nứt xiên chéo góc $45^\circ$ phát triển liên kết trên thân vách ở các chu kỳ tải trọng đỉnh ($\Delta \ge 25\text{ mm}$).
  4. Quy luật biến thiên của hệ số cản nhớt tương đương ($\xi_{eq}$): Hệ số cản nhớt tương đương phản ánh năng lượng tiêu tán trễ tăng phi tuyến từ giá trị ban đầu $\xi \approx 3,5% - 5%$ (giai đoạn đàn hồi) lên tới $\xi \approx 14,8% - 17,2%$ khi vách đạt trạng thái biến dạng dẻo cực hạn ($\mu_\Delta \ge 3,5$).
  5. Hiệu năng và độ tin cậy vượt trội của chương trình CSW: So sánh giữa kết quả mô phỏng bằng CSW và số liệu thực nghiệm trên hai mẫu vách V1, V2 cho thấy sự trùng khớp rất cao về đường bao lực - chuyển vị ($P_{max, exp} = 245\text{ kN}$ so với $P_{max, num} = 236\text{ kN}$, sai số chỉ $3,67%$). Khi đối chuẩn với thực nghiệm của Thomsen và Wallace, chương trình CSW mô phỏng chính xác mức độ dẻo và hình thái phá hoại với sai số cực đại dưới $7,5%$.
               QUAN HỆ TẢI TRỌNG NGANG - CHUYỂN VỊ ĐỈNH (ĐƯỜNG TRỄ)
        P (kN) ^
               |                 * * * Đường bao đỉnh tải
         245 --+             * * *   /
               |         * *        /  Chu kỳ gia tải lặp
               |       *           /
               |     /           /
               |    /   /\      /
               |   /   /  \    /
         ------+--+---+----+--+-------------------> Chuyển vị Delta (mm)
               | /   /      \/       10    20    30    40
               |/   /       /
              /|   /       /
             / |  /       /
            /  * *       /  <- Hiệu ứng kẹp trễ (Pinching)
        -245 -+    * * *
               |

Implications đa chiều

  • Về mặt lý thuyết kết cấu: Luận án hoàn thiện cơ sở khoa học về việc kết hợp mô hình vết nứt phân tán và phần tử tấm hữu hạn hỗn hợp, chứng minh tính khả thi của việc mô phỏng ứng xử trễ phức tạp mà không cần dựa hoàn toàn vào các liên kết lò xo kinh nghiệm đơn giản.
  • Về mặt phương pháp luận: Thiết lập quy trình thí nghiệm chuẩn mực chịu tải trọng lặp đa chu kỳ cho cấu kiện vách mỏng quy mô lớn, cung cấp hệ thống công thức giải tích xác định chính xác góc xoay và độ cong phân bố dọc chiều cao vách.
  • Về mặt ứng dụng kỹ thuật và thực tiễn thiết kế: Cung cấp khuyến nghị trực tiếp cho kỹ sư kết cấu khi áp dụng TCVN 5574:2018 và tiêu chuẩn kháng chấn:
    • Bắt buộc phải cấu tạo vùng biên tăng cường chịu mô men với chiều dài $B_t, B_p \ge 0,2 l_w$ và bố trí đai giữ chặt với bước đai $s_h \le 100\text{ mm}$ để ngăn ngừa sự mất ổn định uốn dọc của cốt thép chịu nén khi chịu tải lặp.
    • Cần xét đến sự suy giảm độ cứng chống cắt khi tính toán chu kỳ dao động riêng và chuyển vị đỉnh công trình dưới tác động của động đất.
  • Về mặt chính sách và tiêu chuẩn hóa: Kết quả nghiên cứu đóng góp cơ sở dữ liệu thực nghiệm quan trọng cho Hội đồng Tiêu chuẩn Xây dựng Quốc gia nhằm biên soạn, cập nhật các hướng dẫn thiết kế kháng chấn theo phương pháp chuyển vị (Displacement-Based Seismic Design) trong các phiên bản cập nhật của TCVN.

Limitations và Future Research

Luận án thừa nhận một cách khách quan các giới hạn nghiên cứu:

  • Quy mô mẫu thực nghiệm: Mặc dù các mẫu V1 và V2 có kích thước lớn ($h_w = 2300\text{ mm}$), số lượng mẫu thực nghiệm chịu tải lặp còn khiêm tốn (02 mẫu chuyên sâu) do chi phí và điều kiện thiết bị thí nghiệm tải trọng nặng trong nước.
  • Phạm vi lực nén dọc trục: Mới chỉ khảo sát ở mức tỷ số nén dọc trục trung bình ($n = N / (f_c' A_g) \approx 0,15$), chưa khảo sát các cấp lực nén rất cao ($n \ge 0,35$) vốn có thể dẫn đến cơ chế phá hoại nén giòn đột ngột.
  • Hình học tiết diện: Nghiên cứu tập trung chủ yếu vào vách phẳng tiết diện chữ nhật, chưa mở rộng thực nghiệm đối với các dạng vách có hình học phức tạp như vách chữ L, chữ C (lõi thang máy) hoặc vách có khoét lỗ cửa lớn.
  • Môi trường tính toán: Chương trình CSW hiện tại tối ưu cho bài toán ứng suất phẳng 2D, chưa tích hợp tự động vào không gian 3D tương tác toàn công trình có xét đến hiệu ứng dầm nối (coupling beam).

Hướng phát triển nghiên cứu tương lai:

  1. Mở rộng phần tử hỗn hợp sang mô hình tấm vỏ 3 chiều (3D Shell Element) để phân tích phi tuyến các hệ vách lõi hở và vách có cánh.
  2. Nghiên cứu thực nghiệm ứng xử của vách BTCT sử dụng bê tông tính năng siêu cao (UHPC) hoặc bê tông cốt sợi phân tán (FRC) chịu tải trọng động đất mạnh.
  3. Tích hợp thuật toán tổn thương mỏi và hiện tượng trượt cốt thép (bond-slip) vào ma trận cấu tạo đàn dẻo của chương trình CSW.
  4. Xây dựng mô hình trí tuệ nhân tạo (Machine Learning) huấn luyện trên cơ sở dữ liệu của chương trình CSW để dự báo nhanh biểu đồ trễ hysteretic cho kỹ sư thiết kế.

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

Nghiên cứu tạo ra những tác động lan tỏa trên nhiều phương diện:

  • Tác động học thuật (Academic Impact): Cung cấp tài liệu tham khảo chuyên sâu và bộ số liệu đo đạc thực nghiệm chuẩn xác cho các nghiên cứu sinh, học viên cao học chuyên ngành Kết cấu Công trình Xây dựng. Các thuật toán trong luận án là nền tảng tốt cho các công bố quốc tế trên các tạp chí thuộc danh mục ISI/Scopus (như Engineering Structures, Journal of Structural Engineering - ASCE).
  • Chuyển đổi công nghiệp xây dựng: Giúp các đơn vị tư vấn thiết kế lớn tại Việt Nam (VNCC, CDC, Coninco, Texo...) có thêm công cụ kiểm chứng tin cậy khi phân tích an toàn phi tuyến cho các dự án tháp cao tầng từ 30 đến 50 tầng, tối ưu hóa khối lượng cốt thép vùng biên mà vẫn đảm bảo độ dẻo kháng chấn.
  • Tác động kinh tế - xã hội: Việc tối ưu hóa thiết kế vách BTCT theo khả năng biến dạng dẻo giúp giảm lượng bê tông và cốt thép tiêu thụ từ $8%$ đến $12%$ trong kết cấu phần thân, tiết kiệm hàng tỷ đồng chi phí xây dựng cho mỗi công trình, đồng thời giảm lượng phát thải carbon (CO2) từ sản xuất xi măng và thép, đóng góp vào mục tiêu phát triển công trình xanh bền vững.

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

  • Nghiên cứu sinh và Giảng viên đại học: Tiếp cận phương pháp luận nghiên cứu kết hợp giải tích phần tử hữu hạn và thực nghiệm đa kênh; sử dụng dữ liệu luận án làm chuẩn so sánh cho các mô hình số mới.
  • Kỹ sư kết cấu thực hành: Nắm vững bản chất cơ học của vách chịu tải lặp, hiểu rõ cơ chế làm việc của vùng biên chịu mô men để cấu tạo cốt thép chính xác theo TCVN 5574:2018 và ACI 318.
  • Các viện nghiên cứu và Cơ quan quản lý nhà nước (Bộ Xây dựng, Viện IBST): Có luận cứ khoa học thực nghiệm vững chắc để hoàn thiện hệ thống tiêu chuẩn kỹ thuật quốc gia về thiết kế kháng chấn nhà cao tầng.
  • Chủ đầu tư và Nhà thầu thi công: Nâng cao mức độ an toàn công trình trước hiểm họa thiên tai động đất, tối ưu hóa chi phí đầu tư xây dựng.

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ì? Đó là việc phát triển thành công mô hình cấu tạo đàn dẻo phi tuyến cho vật liệu hỗn hợp bê tông cốt thép trong không gian ứng suất phẳng, kết hợp quy luật trễ hysteretic cải tiến (xét hiệu ứng kẹp nứt và đóng mở khe nứt) và mở rộng mô hình kiềm chế nở hông của Kent - Park và Mander. Mô hình này phản ánh chính xác trạng thái suy giảm độ cứng và phân tán năng lượng mà không cần sử dụng các tham số kinh nghiệm gán ghép nhân tạo.

  2. Cải tiến phương pháp luận so với các nghiên cứu trước đây thể hiện ở điểm nào? So với nghiên cứu của Lefas và cộng sự (1990) hay Thomsen và Wallace (1995), luận án đã tự xây dựng một công cụ phần mềm độc lập (CSW) tích hợp phần tử tấm tứ giác đẳng tham số đa nút nhúng cốt thép trực tiếp, giải quyết bài toán cân bằng phi tuyến bằng thuật toán Newton - Raphson cải tiến và Newmark. Quy trình thực nghiệm tại Việt Nam được trang bị hệ thống đo đạc LVDT đa điểm đồng bộ, tách bạch được biến dạng uốn, xoay chân vách và biến dạng trượt cắt thuần túy.

  3. Phát hiện bất ngờ nhất từ dữ liệu thực nghiệm là gì? Sự gia tăng đột biến của biến dạng trượt cắt $\gamma$ ở giai đoạn phi tuyến sâu (chiếm tới gần $40%$ tổng chuyển vị đỉnh) ngay cả đối với vách có độ thanh mảnh cao ($h_w/l_w \approx 2,0$). Điều này chứng minh rằng giả thiết cổ điển coi vách cao chỉ chịu uốn thuần túy và bỏ qua tương tác trượt cắt là chưa hoàn toàn an toàn trong thiết kế kháng chấn dựa trên chuyển vị.

  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hay không? Hoàn toàn có. Toàn bộ thông số cấp phối bê tông B20, quy cách gia công thép, kích thước hình học chi tiết của mẫu V1, V2, sơ đồ bố trí cảm biến LVDT, strain gauges, quy trình gia tải điều khiển chuyển vị và toàn bộ sơ đồ khối thuật toán của chương trình CSW đều được trình bày minh bạch, chi tiết trong văn bản luận án và phần phụ lục mã nguồn.

  5. Chương trình nghị sự nghiên cứu 10 năm tới được định hình ra sao? Tập trung vào 3 trục chính: (1) Phát triển phần tử vỏ 3D phi tuyến cho vách lõi không gian; (2) Nghiên cứu ứng xử của vách bê tông tự lèn cường độ siêu cao (UHPC) kết hợp cáp ứng suất trước không bám dính tự phục hồi sau động đất (Self-Centering Shear Walls); (3) Tích hợp học máy để đánh giá độ tin cậy và tổn thương kết cấu theo thời gian thực.

Kết luận

  1. Xây dựng thành công phần tử hữu hạn hỗn hợp bê tông cốt thép phi tuyến: Tích hợp mô hình vết nứt phân tán và các đường cong ứng suất - biến dạng thực nghiệm, giải quyết trọn vẹn bài toán mô phỏng ứng xử tấm phẳng chịu tải trọng lặp đảo chiều.
  2. Thực hiện thành công chuỗi thực nghiệm vách BTCT quy mô lớn: Cung cấp bộ số liệu thực nghiệm toàn diện, phản ánh chính xác các quy luật cơ học then chốt: quan hệ lực - chuyển vị đỉnh, phân bố độ cong theo chiều cao, góc xoay chân vách, biến dạng trượt cắt, quy luật suy giảm độ cứng và hệ số cản nhớt tương đương ($\xi = 3,5% - 17,2%$).
  3. Phát triển chương trình tính toán chuyên dụng CSW: Công cụ phần mềm có độ chính xác cao (sai số so với thực nghiệm dưới $7,5%$), kiểm chứng thành công trên dữ liệu thực nghiệm trong nước và dữ liệu quốc tế của Thomsen - Wallace.
  4. Đóng góp chỉ dẫn thiết kế thực tế: Xác lập cơ sở khoa học cho việc tính toán và cấu tạo vùng biên chịu mô men ($B_t = 0,2 l_w$) và hàm lượng cốt thép đai chống nở hông theo tiêu chuẩn TCVN 5574:2018.
  5. Thúc đẩy chuyển dịch mô hình thiết kế: Đặt nền móng kỹ thuật vững chắc cho việc chuyển đổi từ phương pháp thiết kế theo ứng suất/độ bền truyền thống sang phương pháp thiết kế dựa trên chuyển vị và hiệu năng (Performance/Displacement-Based Design) cho nhà cao tầng tại Việt Nam.
  6. Mở ra các hướng nghiên cứu kết cấu tiên tiến: Tạo tiền đề mở rộng mô phỏng không gian 3D và ứng dụng các vật liệu bê tông thế hệ mới trong công trình xây dựng hiện đại.