Tổng quan về luận án

Trong kỹ thuật xây dựng công trình hiện đại, thiết kế kháng chấn cho nhà nhiều tầng đang chuyển dịch mạnh mẽ từ trường phái tiếp cận thuần túy theo khả năng chịu tải trọng sang tiếp cận dựa trên năng lượng và kiểm soát tính năng công trình. Bản chất phản ứng của kết cấu công trình dưới tác động của tải trọng động đất là bài toán cân bằng và phân tán năng lượng: $\text{Năng lượng truyền vào} = \text{Năng lượng đàn hồi} + \text{Năng lượng động năng} + \text{Năng lượng cản} + \text{Năng lượng tiêu tán trễ}$. Để công trình không sụp đổ khi xảy ra động đất mạnh, hệ thống kết cấu phải sở hữu cơ chế hấp thu và tiêu tán năng lượng hiệu quả, duy trì khả năng chịu tải và hạn chế tối đa hư hại vĩnh viễn ở các cấu kiện chịu lực chính.

Tường chịu cắt bê tông cốt thép (BTCT) là giải pháp chịu tải trọng ngang phổ biến nhất trong các tòa nhà cao tầng nhờ độ cứng lớn và khả năng chống cắt vượt trội. Khi quy mô công trình mở rộng hoặc khi xuất hiện các dãy lỗ thông tầng phục vụ yêu cầu kiến trúc (như cửa đi, cửa sổ, hành lang), hệ thống sẽ hình thành kết cấu tường kép (coupled shear walls). Các vách tường đơn được liên kết với nhau thông qua hệ dầm nối (coupling beams). Theo tiêu chuẩn thiết kế kháng chấn Việt Nam TCVN 9386:2012, kết cấu tường kép thể hiện đặc tính tiêu tán năng lượng và khả năng ứng xử vượt trội hơn so với tường đơn chịu cắt thông thường: "kết cấu tường kép có hệ số ứng xử cao hơn kết cấu tường chịu cắt thông thường, có nghĩa là về mặt phản ứng với tác động của tường kép tốt hơn. Cụ thể, khi thiết kế nhà nhiều tầng, nhiều nhịp theo cấp dẻo kết cấu trung bình, hệ số ứng xử của hệ tường kép và tường thông thường tương ứng là 3,9 và 3,0. Với cấp dẻo kết cấu cao, hệ số ứng xử là 5,85 và 5,2 tương ứng với hệ tường kép và hệ không thuộc hệ tường kép".

Tuy nhiên, dầm nối BTCT truyền thống thường chịu lực cắt rất lớn và chuyển vị góc xoay đảo chiều liên tục trong động đất, dẫn đến nguy cơ phá hoại giòn do cắt nứt chéo hoặc phá hoại trượt ở đầu dầm, gây suy giảm đột ngột độ cứng toàn công trình. Dù các giải pháp bố trí cốt thép chéo hay cốt thép hình thoi đã được nghiên cứu, việc thi công cốt thép phức tạp trong các vách mỏng gây cản trở lớn trên thực địa. Nhằm giải quyết triệt để vấn đề này, giải pháp sử dụng kết cấu tường kép lai (hybrid coupled shear walls) kết hợp tường bê tông cốt thép với dầm nối bằng thép có gắn thiết bị tiêu tán năng lượng (TTNL) đã ra đời. Thiết bị TTNL đóng vai trò như các "cầu chì kết cấu" (structural fuses), tập trung toàn bộ biến dạng dẻo và năng lượng trễ tại thiết bị, giữ cho tường bê tông và dầm thép chính luôn làm việc trong giới hạn đàn hồi.

Khoảng trống nghiên cứu then chốt mà luận án tiến sĩ kỹ thuật "Nghiên cứu tính toán kết cấu tường kép có gắn thiết bị tiêu tán năng lượng chịu tác động động đất" (chuyên ngành Cơ học vật rắn, mã số: 9440107, thực hiện bởi NCS. Phạm Thu Hiền dưới sự hướng dẫn của GS. TS. Nguyễn Tiến Chương và TS. Nguyễn Hải Quang tại Trường Đại học Thủy lợi, năm 2023) tập trung giải quyết là: Các mô hình số hiện đại trên thế giới khi phân tích dầm nối có gắn thiết bị TTNL phi tuyến đều phải chia nhỏ cấu kiện thành tổ hợp phức tạp gồm nhiều phần tử thanh, phần tử dầm cứng và nhiều phần tử liên kết phi tuyến có chiều dài bằng 0 (zero-length elements), làm tăng vọt bậc tự do, gây khó khăn lớn cho việc tính toán động lực học phi tuyến lịch sử thời gian và chưa đánh giá được toàn diện ảnh hưởng của các tham số kỹ thuật, số lượng cũng như vị trí lắp đặt thiết bị theo chiều cao công trình.

Luận án xác lập 3 câu hỏi nghiên cứu và các giả thuyết khoa học cụ thể:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập một mô hình tính toán giải tích - phần tử hữu hạn thu gọn cho dầm nối có gắn thiết bị TTNL phi tuyến trễ thành một siêu phần tử đơn nhất mà vẫn bảo toàn đầy đủ ứng xử cơ học phi tuyến phức tạp?
    • Giả thuyết 1 (H1): Toàn bộ dầm nối thép tích hợp thiết bị TTNL ở giữa có thể được mô hình hóa bằng một phần tử thanh tương đương duy nhất thông qua việc thiết lập ma trận độ cứng cục bộ tích hợp quy luật quan hệ lực cắt - biến dạng trượt phi tuyến trễ và chuyển đổi tọa độ ma trận độ cứng tổng thể.
  2. Câu hỏi nghiên cứu 2 (RQ2): Ảnh hưởng định lượng của các đặc trưng kỹ thuật của thiết bị TTNL (lực chảy dẻo $Q_y$, độ cứng đàn hồi $k_1$, độ cứng sau chảy $k_2$) và độ cứng hình học dầm thép đến phản ứng giảm chấn của hệ tường kép là gì?
    • Giả thuyết 2 (H2): Tồn tại một ngưỡng lực chảy dẻo $Q_y$ và độ cứng thiết bị tối ưu; việc vượt quá hoặc thấp hơn ngưỡng này sẽ làm suy giảm đáng kể hiệu quả tiêu tán năng lượng của hệ kết cấu.
  3. Câu hỏi nghiên cứu 3 (RQ3): Số lượng và vị trí bố trí thiết bị TTNL theo chiều cao công trình ảnh hưởng như thế nào đến khả năng kiểm soát chuyển vị đỉnh, gia tốc đỉnh và chuyển vị ngang tương đối giữa các tầng dưới các phổ gia tốc động đất khác nhau?
    • Giả thuyết 3 (H3): Việc phân bổ thiết bị tập trung tại các tầng có lực cắt và góc xoay dầm nối lớn nhất (dải tầng trung tâm) sẽ mang lại hiệu quả giảm chấn tương đương việc lắp đặt trên 100% các tầng, giúp tiết kiệm đáng kể chi phí kinh tế.

Phạm vi nghiên cứu của luận án bao gồm việc mô hình hóa lý thuyết, xây dựng chương trình số chuyên dụng TTKGC (Tính toán khung gắn thiết bị giảm chấn), kiểm chứng với phần mềm tiêu chuẩn quốc tế SAP2000 và phân tích phản ứng động lực học phi tuyến của hệ tường kép 10 tầng và 30 tầng dưới tác động của 4 phổ gia tốc nền động đất thực tế (KOBE, TCU078, IND, CAST360_DBE).

Literature Review và Positioning

Lịch sử phát triển các giải pháp tăng cường tính năng kháng chấn cho kết cấu tường kép trải qua 4 giai đoạn tiến hóa lý thuyết và công nghệ rõ rệt:

  1. Giai đoạn nghiên cứu cấu tạo cốt thép trong dầm nối bê tông cốt thép: Khởi đầu từ các nghiên cứu kinh điển của Paulay và cộng sự (1969), các dạng phá hoại dầm nối BTCT chịu uốn - cắt chu kỳ đã được làm rõ, bao gồm phá hoại kéo theo đường chéo và phá hoại trượt do cắt. Để khắc phục hiện tượng phá hoại giòn này, Paulay và Binney (1974) đã đề xuất giải pháp bố trí cốt thép chéo trong dầm nối. Mặc dù nâng cao đáng kể độ dẻo và khả năng tiêu tán năng lượng, nhược điểm cốt tử của cốt thép chéo là sự xung đột cốt thép nghiêm trọng, gây khó khăn cho công tác đổ bê tông khi chiều dày tường $t < 250\text{ mm}$. Đến các nghiên cứu thực nghiệm của Galano và Vignoli, tác giả đã thí nghiệm 15 mẫu dầm nối có tỷ lệ nhịp trên chiều cao $L/h = 1,5$ với các cấu tạo cốt thép khác nhau (cốt thép thông thường, cốt thép chéo không đai, cốt thép chéo có đai và cốt thép bố trí hình thoi). Dữ liệu thực nghiệm của Galano và Vignoli chứng minh: "dầm nối có cấu tạo cốt thép hình thoi đạt giá trị độ dẻo cao nhất. Tuy nhiên, dầm nối bố trí cốt thép hình thoi có độ bền thấp hơn dầm khác có cùng hàm lượng cốt thép (giảm xấp xỉ 17 (%))... Về khả năng tiêu tán năng lượng, dầm nối với cốt thép hình thoi cho thấy khả năng tiêu tán năng lượng nhiều hơn 9,6 (%) so với năng lượng được tiêu tán bởi các mẫu thí nghiệm có cấu tạo cốt thép đường chéo". Mặc dù vậy, cấu tạo này vẫn không giải quyết được vấn đề tập trung hư hại và biến dạng dư vĩnh viễn trên thân dầm sau động đất.

  2. Giai đoạn dầm nối thép và dầm nối liên hợp thép - bê tông: Nhằm tăng tính dẻo và thuận tiện hóa thi công, Harries và cộng sự đã tiên phong nghiên cứu kết cấu tường kép lai liên hợp tường BTCT và dầm nối bằng thép hình. Tiếp đó, Gong và Shahrooz đã đề xuất dầm nối liên hợp thép - bê tông bọc ngoài, chứng minh rằng sự kết hợp giữa bản bụng thép chữ I và bê tông bọc xung quanh giúp ngăn chặn hiện tượng oằn cục bộ của bản bụng thép, đồng thời tăng cường độ cứng liên kết của hệ tường kép.

  3. Giai đoạn phát triển thiết bị giảm chấn và cầu chì năng lượng: Các giải pháp giảm chấn thụ động, chủ động và bán chủ động đã được phát triển rộng rãi. Trong khi phương pháp cách chấn đáy (Base Isolation) bộc lộ nhược điểm là làm kết cấu quá mềm theo phương ngang, không hiệu quả trước tải trọng gió đối với công trình siêu cao tầng, thì thiết bị tiêu tán năng lượng gắn trực tiếp tại dầm nối thể hiện ưu thế vượt trội. Các dạng thiết bị giảm chấn thụ động bao gồm thiết bị cản kim loại (Metallic Dampers - MD), thiết bị cản ma sát (Friction Dampers - FD), thiết bị cản nhớt (Viscous Dampers - VD), thiết bị cản khối lượng điều chỉnh (Tuned Mass Dampers - TMD). Dẫn chứng về TMD, Villaverde đã nghiên cứu và chỉ ra rằng: "thiết bị TMD gắn tại đỉnh đã giảm được tới 84 (%) năng lượng động đất tác dụng vào kết cấu, điều này cũng có nghĩa là việc giảm năng lượng động đất đáng kể chỉ có thể xảy ra khi đặt thiết bị TMD ở trên đỉnh của công trình", tiêu biểu như quả cầu thép 660 tấn tại tháp Taipei 101. Đối với dầm nối, Fortney đã phát triển dầm nối thép có khớp nối trung tâm (fuse) liên kết trượt bằng bu lông; Ahn và cộng sự tích hợp thiết bị cản ma sát kết hợp vật liệu NAO và cản kim loại; Montgomery và cộng sự đề xuất dầm nối vật liệu đàn nhớt; S. Oh và cộng sự kết hợp cao su cản cao và thép chữ U; H. Jiang và cộng sự kết hợp cản kim loại và cản nhớt.

  4. Giai đoạn mô phỏng số phi tuyến kết cấu tường kép có gắn thiết bị TTNL: Các nghiên cứu quốc tế gần đây đã ứng dụng phần mềm số để phân tích ứng xử động học của hệ kết cấu này, song đều tồn tại hạn chế lớn về mô hình hóa:

    • Cong Chen và cộng sự (2017): Sử dụng OpenSees mô phỏng dầm nối gắn thiết bị TTNL bằng mô hình giàn song song (parallel truss), đòi hỏi tới 5 phần tử giàn và 10 phần tử liên kết cứng cho một nút dầm.
    • Yun Chen và cộng sự (2019): Sử dụng ABAQUS mô hình hóa phần tử khối 3D C3D8R 8 nút cho cả tường, dầm và thiết bị với mô hình tái bền đẳng hướng, dẫn đến khối lượng tính toán khổng lồ.
    • Molham Salameh và cộng sự (2022): Mô hình hóa dầm nối gắn thiết bị TTNL trên OpenSees bằng cách sử dụng 1 phần tử dầm Euler-Bernoulli kết hợp 2 phần tử liên kết zero-length ở mỗi vị trí liên kết. Tổng cộng để mô phỏng một dầm nối cần tới 3 phần tử thanh và 4 phần tử liên kết zero-length.
    • A. Kheyroddin và cộng sự (2022): Mô hình hóa dầm nối gồm 3 phần, trong đó thiết bị TTNL là 1 phần tử liên kết zero-length phi tuyến cắt, kết nối với các đoạn dầm bằng các liên kết phụ, đòi hỏi 2 phần tử thanh và 3 phần tử liên kết zero-length.
                    ┌────────────────────────────────────────────────────────┐
                    │     TIẾN TRÌNH NGHIÊN CỨU DẦM NỐI TƯỜNG KÉP KHÁNG CHẤN │
                    └───────────────────────────┬────────────────────────────┘
                                                │
         ┌──────────────────────────────────────┴──────────────────────────────────────┐
         ▼                                                                             ▼
┌─────────────────────────────────┐                                           ┌─────────────────────────────────┐
│     HƯỚNG TIẾP CẬN KẾT CẤU      │                                           │      HƯỚNG MÔ HÌNH HÓA SỐ       │
├─────────────────────────────────┤                                           ├─────────────────────────────────┤
│ • Paulay & Binney (1974):       │                                           │ • Yun Chen et al. (2019):       │
│   Cốt thép chéo (khó thi công)  │                                           │   Phần tử khối 3D C3D8R (rất nặng)│
│ • Galano & Vignoli:             │                                           │ • Cong Chen et al. (2017):      │
│   Thép hình thoi (+9.6% NL, -17% B)│                                        │   5 thanh giàn + 10 LK cứng     │
│ • Harries et al., Gong et al.:  │                                           │ • Salameh et al. (2022):        │
│   Dầm thép liên hợp             │                                           │   3 phần tử thanh + 4 LK zero-L │
│ • Fortney, Ahn, Oh, Jiang:      │                                           │ • Kheyroddin et al. (2022):     │
│   Gắn thiết bị TTNL / Khớp nối  │                                           │   2 phần tử thanh + 3 LK zero-L │
└────────────────┬────────────────┘                                           └────────────────┬────────────────┘
                 │                                                                             │
                 └──────────────────────────────────────┬──────────────────────────────────────┘
                                                        ▼
                                   ┌────────────────────────────────────────┐
                                   │       ĐỘT PHÁ CỦA LUẬN ÁN (2023)       │
                                   ├────────────────────────────────────────┤
                                   │ • Đóng gói Siêu phần tử dầm nối duy nhất│
                                   │ • Ma trận độ cứng giải tích phi tuyến   │
                                   │ • Chương trình TTKGC chuyên dụng        │
                                   └────────────────────────────────────────┘

Positioning khoa học của luận án: Luận án định vị tại điểm giao thoa giữa cơ học công trình phi tuyến và kỹ thuật phần tử hữu hạn nâng cao. Khắc phục sự cồng kềnh của các mô hình đa phần tử trong OpenSees của Salameh (2022) và Kheyroddin (2022), luận án đề xuất một ma trận độ cứng giải tích thu gọn cho một siêu phần tử dầm nối tích hợp trực tiếp thiết bị TTNL, cho phép phân tích phi tuyến lịch sử thời gian nhanh chóng, chính xác và mở đường cho việc khảo sát tối ưu hóa đa tham số trên các công trình nhiều tầng thực tế.

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

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

Luận án tạo ra bước tiến quan trọng trong lý thuyết tính toán kết cấu công trình chịu động đất thông qua các đóng góp cốt lõi:

  1. Mở rộng lý thuyết phần tử dầm tương đương có xét vùng cứng vô cùng và khớp dẻo trượt phi tuyến: Luận án đã phát triển công thức ma trận độ cứng cho phần tử thanh gồm hai đoạn cứng ở hai đầu (đoạn dầm chôn sâu vào vách bê tông và đoạn dầm thép liên kết) và một vùng biến dạng trượt phi tuyến trễ ở giữa.
  2. Thiết lập mô hình giải tích cho quy luật quan hệ lực cắt - biến dạng trượt ($Q - \Delta$) của thiết bị TTNL đa trạng thái: Mô hình toán học mô tả chính xác 6 trạng thái động học của vòng lặp trễ song tuyến tính (Bilinear Hysteresis Model):
    • Trạng thái 1: Lực cắt $Q$ tăng trong miền đàn hồi với độ cứng ban đầu $k_1$.
    • Trạng thái 2: Lực cắt $Q$ vượt qua lực chảy dẻo $Q_y$, chuyển sang giai đoạn dẻo với độ cứng sau chảy $k_2$.
    • Trạng thái 3: Quá trình dỡ tải (lực cắt $Q$ giảm) với độ cứng đàn hồi phục hồi $k_1$.
    • Trạng thái 4: Quá trình dỡ tải tiếp tục chuyển sang giai đoạn biến dạng dẻo âm với độ cứng $k_2$.
    • Trạng thái 5: Quá trình đảo chiều biến dạng chuyển tiếp từ độ cứng $k_2$ sang $k_1$.
    • Trạng thái 6: Quá trình đảo chiều biến dạng chuyển tiếp từ $k_1$ sang $k_2$.

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

Khung phân tích của luận án là sự tích hợp chặt chẽ giữa 3 lý thuyết nền tảng:

  • Lý thuyết động lực học công trình phi tuyến: Thiết lập phương trình vi phân chuyển động tổng quát của hệ nhiều bậc tự do chịu kích động gia tốc nền động đất: $$[M]{\ddot{u}(t)} + [C]{\dot{u}(t)} + [K(t)]{u(t)} = -[M]{I}\ddot{u}_g(t)$$ Trong đó $[M]$ là ma trận khối lượng tập trung, $[C]$ là ma trận cản Rayleigh ($[C] = \alpha_m [M] + \beta_k [K]$), $[K(t)]$ là ma trận độ cứng tức thời thay đổi theo từng bước thời gian phụ thuộc trạng thái biến dạng của thiết bị TTNL, và $\ddot{u}_g(t)$ là gia tốc nền động đất.
  • Lý thuyết năng lượng bổ sung (Complementary Energy Theorem) và Nguyên lý công ảo: Sử dụng để vi phân thiết lập chính xác các hệ số trong ma trận độ cứng địa phương $[K_r]$ của siêu phần tử dầm nối có thiết bị TTNL.
  • Phương pháp chuyển đổi tọa độ ma trận phần tử hữu hạn: Chuyển đổi ma trận độ cứng từ hệ tọa độ cục bộ $[K_r]$ sang hệ tọa độ tổng thể $[K_c]$ có xét đến góc nghiêng phần tử và các bậc tự do liên kết tại nút: $$[K_c] = [T]^T [K_r] [T]$$ Phương pháp này loại bỏ hoàn toàn nhu cầu sử dụng các phần tử liên kết zero-length trung gian, giảm số lượng phương trình đại số cần giải trong mỗi vòng lặp phi tuyến.

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

Thiết kế nghiên cứu

Luận án tuân thủ lập trường nhận thức luận thực chứng (Positivism), kết hợp suy diễn toán học cơ học giải tích với mô phỏng số phần tử hữu hạn nâng cao và đối soát kiểm chứng chéo đa cấp (Multilevel Cross-Verification).

                                  QUY TRÌNH NGHIÊN CỨU VÀ ĐỐI CHỨNG ĐA CẤP
  ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 1. XÂY DỰNG LÝ THUYẾT VÀ MA TRẬN ĐỘ CỨNG SIÊU PHẦN TỬ                                                  │
  │    Thiết lập ma trận độ cứng [Kr] của dầm nối có thiết bị TTNL -> Chuyển đổi tọa độ [Kc]               │
  └───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                      │
                                                      ▼
  ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 2. XÂY DỰNG THUẬT TOÁN VÀ PHẦN MỀM CHUYÊN DỤNG TTKGC                                                  │
  │    Tích hợp phương pháp số Newmark-beta, cập nhật ma trận [K(t)] theo 6 trạng thái trễ                │
  └───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                      │
                                                      ▼
  ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 3. QUY TRÌNH KIỂM CHỨNG ĐA CẤP (TRIANGULATION & BENCHMARKING)                                          │
  │  ┌───────────────────────────────┐ ┌───────────────────────────────┐ ┌───────────────────────────────┐ │
  │  │ Cấp 1: Kiểm chứng giải tích   │ │ Cấp 2: Kiểm chứng hội tụ lưới │ │ Cấp 3: Kiểm chứng SAP2000     │ │
  │  │ Khung 1 tầng chịu tải điều hòa│ │ Mô hình tấm 14.400 phần tử    │ │ Hệ tường kép 10 tầng gắn TTNL │ │
  │  │ Sai số chuyển vị: ~0.00%      │ │ (15.040 nút), sai lệch 0.9%   │ │ Khớp hoàn hảo lịch sử thời gian│ │
  │  └───────────────────────────────┘ └───────────────────────────────┘ └───────────────────────────────┘ │
  └───────────────────────────────────────────────────┬────────────────────────────────────────────────────┘
                                                      │
                                                      ▼
  ┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
  │ 4. KHẢO SÁT ỨNG XỬ ĐỘNG LỰC HỌC CÔNG TRÌNH CAO TẦNG (30 TẦNG)                                          │
  │    Phân tích 4 sóng động đất (KOBE, TCU078, IND, CAST360_DBE) x Khảo sát Qy, k1, k2 x Khảo sát vị trí │
  └────────────────────────────────────────────────────────────────────────────────────────────────────────┘

Mô hình tường kép được phân tích dựa trên mô hình khung tương đương (Equivalent Frame Model). Để xác lập tính hợp lý của mô hình khung tương đương so với mô hình tấm 2D (Shell Element Model), luận án đã thực hiện nghiên cứu độ hội tụ lưới trên hệ tường kép 10 tầng, chiều cao tầng 3,5 m, chiều dày tường $t = 0,25\text{ m}$, mô đun đàn hồi bê tông $E = 25.000\text{ MPa}$, hệ số Poisson $\nu = 0,2$ chịu tải trọng ngang $200\text{ kN}$ tại các tầng:

  • Khi chia lưới phần tử tấm với số lượng phần tử trên một tầng tăng dần từ 5, 13, 60, 90, 360 đến 1440 phần tử, chuyển vị ngang tại đỉnh tường lần lượt đạt $6,0547\text{ cm}$; $6,2300\text{ cm}$; $6,5335\text{ cm}$; $6,6006\text{ cm}$; $6,7330\text{ cm}$ và $6,7933\text{ cm}$.
  • Kết quả cho thấy khi tăng từ 360 lên 1440 phần tử/tầng (tương ứng toàn công trình có 15.040 nút và 14.400 phần tử tấm), độ chênh lệch chuyển vị đỉnh chỉ là $0,9%$, khẳng định trạng thái hội tụ nghiệm của mô hình số. Mô hình khung tương đương rút gọn thành công hệ thống từ 15.040 nút xuống còn 22 nút và 30 phần tử thanh, giảm thời gian tính toán hàng trăm lần mà vẫn đảm bảo sai số cơ học nằm trong phạm vi kỹ thuật cho phép.

Quy trình nghiên cứu rigorous và thuật toán phần mềm TTKGC

Tác giả đã lập trình hoàn chỉnh phần mềm TTKGC dựa trên thuật toán tích phân từng bước Newmark kết hợp vòng lặp điều chỉnh độ cứng tiếp tuyến (Tangent Stiffness Method) để giải hệ phương trình vi phân chuyển động phi tuyến.

Quy trình kiểm chứng được tiến hành 3 bước nghiêm ngặt:

  1. Kiểm chứng mức độ cấu kiện cơ sở: So sánh kết quả tính toán chu kỳ dao động riêng và chuyển vị đỉnh của khung 1 tầng có dầm nối gắn thiết bị TTNL giữa chương trình TTKGC và lời giải giải tích chính xác, cho thấy sự trùng khít hoàn toàn về mặt số học.
  2. Kiểm chứng phản ứng động lực học với SAP2000 trên kết cấu 10 tầng: Mô hình hóa hệ tường kép 10 tầng gắn thiết bị TTNL tại tầng 1 và từ tầng 3 đến tầng 5, chịu tải trọng động đất CAST360_DBE. Kết quả so sánh giữa TTKGC và SAP2000 về các giá trị lớn nhất của chuyển vị đỉnh ($d_{\max}$), gia tốc đỉnh ($a_{\max}$), lực cắt đáy ($Q_{\text{base}}$), mô men chân vách ($M_{\text{base}}$) và chuyển vị ngang tương đối theo tầng thể hiện mức độ tương đồng gần như tuyệt đối, trong khi thời gian phân tích của TTKGC vượt trội hơn hẳn.
  3. Phân tích kết cấu quy mô lớn: Khảo sát toàn diện hệ tường kép 30 tầng chịu tác động của 4 trận động đất có đặc trưng phổ khác biệt sâu sắc: Kobe (1995), Chi-Chi TCU078 (1999), India (IND) và CAST360_DBE.

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

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

  1. Hiệu năng tiêu tán năng lượng và giảm chấn vượt trội: Khi lắp đặt thiết bị TTNL trên tất cả các tầng (từ tầng 1 đến tầng 30), phản ứng động lực học của công trình giảm đáng kể trên toàn bộ các chỉ tiêu so với kết cấu không gắn thiết bị. Cụ thể, biên độ dao động chuyển vị đỉnh giảm từ $35%$ đến hơn $50%$, gia tốc đỉnh giảm rõ rệt và chuyển vị ngang tương đối giữa các tầng (drift ratio) được kiểm soát chặt chẽ dưới ngưỡng an toàn danh định, loại bỏ nguy cơ va đập và hư hỏng cấu kiện phi kết cấu.
┌──────────────────────────────────────────────────────────────────────────────────────────────────────┐
│ BIỂU ĐỒ SO SÁNH HIỆU QUẢ GIẢM CHẤN KHI GẮN THIẾT BỊ TTNL TRÊN TƯỜNG KÉP 30 TẦNG                     │
├──────────────────────────────────────────────────────────────────────────────────────────────────────┤
│ Chuyển vị đỉnh (cm)                                                                                  │
│ Không TTNL : ████████████████████████████████████████ 100% (Gốc)                                     │
│ Gắn 10 tầng: ████████████████████████ 62% (Giảm 38%)                                                │
│ Gắn 15 tầng: ██████████████████ 48% (Giảm 52%)                                                      │
│ Gắn 30 tầng: ████████████████ 45% (Giảm 55%)                                                        │
├──────────────────────────────────────────────────────────────────────────────────────────────────────┤
│ Ứng xử cấu kiện                                                                                      │
│ Dầm nối BTCT cổ điển: Phá hoại giòn nứt chéo, bê tông vỡ vụn, biến dạng dư vĩnh viễn                 │
│ Giải pháp Luận án   : Dầm thép & Vách BTCT đàn hồi 100%, Thiết bị TTNL hấp thu trọn năng lượng trễ  │
└──────────────────────────────────────────────────────────────────────────────────────────────────────┘

  1. Quy luật phi tuyến về lực chảy dẻo ($Q_y$) của thiết bị: Khảo sát tham số lực chảy qua chuỗi thiết bị từ TB1 đến TB10 cho thấy hiệu quả giảm chấn không tăng tuyến tính theo lực chảy:

    • Nếu $Q_y$ quá nhỏ: Thiết bị chảy dẻo quá sớm tại biên độ kích động thấp, diện tích vòng trễ nhỏ, khả năng hấp thu năng lượng động đất hạn chế.
    • Nếu $Q_y$ quá lớn: Thiết bị không bị kích hoạt biến dạng dẻo trong suốt trận động đất, làm việc thuần túy đàn hồi, dẫn đến kết cấu hoạt động như một hệ dầm cứng nguyên khối, lực truyền trực tiếp vào tường gây tăng mô men chân vách và gia tốc đỉnh.
    • Luận án xác định được miền giá trị lực chảy tối ưu $Q_{y,\text{opt}}$ đảm bảo diện tích vòng lặp trễ $W_D = \oint Q , d\Delta$ đạt cực đại.

  2. Phát hiện mang tính đột phá về vị trí và số lượng bố trí thiết bị tối ưu: Khảo sát các phương án bố trí thiết bị tại các dải tầng khác nhau: $1/5$ số tầng (6 tầng), $1/3$ số tầng (10 tầng), $1/2$ số tầng (15 tầng), $2/3$ số tầng (20 tầng), $5/6$ số tầng (25 tầng) và toàn bộ 30 tầng:

    • Trong kết cấu tường kép, biến dạng trượt và lực cắt của dầm nối phân bố không đều theo chiều cao, đạt giá trị cực đại ở khoảng $1/3$ đến $2/3$ chiều cao công trình (từ tầng 10 đến tầng 20 đối với nhà 30 tầng).
    • Kết quả phân tích chỉ ra rằng: Phương án bố trí thiết bị TTNL tập trung tại $1/2$ đến $2/3$ số tầng ở dải giữa công trình mang lại hiệu quả giảm chuyển vị đỉnh và lực cắt đáy đạt tới $88% - 94%$ so với phương án gắn thiết bị trên toàn bộ 30 tầng. Phát hiện này mang ý nghĩa thực tiễn to lớn, cho phép cắt giảm từ $33%$ đến $50%$ số lượng thiết bị lắp đặt, tối ưu hóa triệt để chi phí đầu tư xây dựng.

  3. Bảo vệ toàn vẹn tường BTCT và khả năng tự phục hồi: Dưới các trận động đất mạnh như Kobe hay Chi-Chi, vòng lặp trễ của thiết bị mở rộng tối đa, chứng minh cơ chế tiêu tán năng lượng biến dạng dẻo thuần túy diễn ra trong thiết bị. Toàn bộ thân dầm thép và vách bê tông cốt thép hoàn toàn không xuất hiện biến dạng dẻo, bảo toàn 100% khả năng chịu lực đàn hồi, cho phép công trình nhanh chóng đưa vào sử dụng lại sau động đất chỉ bằng thao tác kiểm tra hoặc thay thế chốt năng lượng tại dầm nối.

Implications đa chiều

  • Về mặt học thuật: Cung cấp phương pháp luận và công cụ toán học mô hình hóa giải tích thu gọn cho các cấu kiện phi tuyến trễ phức tạp, đóng góp vào lý thuyết kết cấu kháng chấn hiện đại.
  • Về mặt phương pháp luận: Mở ra khả năng ứng dụng mô hình siêu phần tử tương đương cho các loại kết cấu lai khác như khung - giằng giảm chấn, vách liên hợp có lỗ mở bất đối xứng.
  • Về mặt thực tiễn kỹ thuật: Cung cấp cho các kỹ sư tư vấn thiết kế kết cấu cẩm nang định lượng rõ ràng về việc lựa chọn thông số kỹ thuật thiết bị ($Q_y, k_1, k_2$) và sơ đồ phân bổ tối ưu theo tầng, tránh việc lắp đặt dàn trải gây lãng phí kinh phí.

Limitations và Future Research

Mặc dù đạt được những đóng góp khoa học đột phá, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu nội tại:

  1. Giới hạn mô hình phẳng 2D: Phân tích hiện tại được thực hiện trong mặt phẳng khung 2D, chưa xét đến hiệu ứng xoắn không gian 3D khi công trình có mặt bằng bất đối xứng hoặc chịu tác động đồng thời của 3 thành phần gia tốc động đất không gian ($X, Y, Z$).
  2. Giả thiết ứng xử đàn hồi của tường bê tông: Mô hình giả định tường BTCT luôn làm việc trong miền đàn hồi nhờ sự bảo vệ của thiết bị TTNL. Trong kịch bản động đất siêu cấp vượt tần suất thiết kế cực hạn, tường BTCT có thể xuất hiện vết nứt dẻo cục bộ ở chân vách chưa được mô hình hóa phi tuyến phi đàn hồi trong thuật toán hiện tại.
  3. Điều kiện biên ngàm cứng tại chân móng: Luận án chưa xét đến hiệu ứng tương tác động lực học giữa đất nền và kết cấu (Soil-Structure Interaction - SSI), vốn có thể làm thay đổi chu kỳ dao động riêng và gia tăng chuyển vị xoay ở chân móng trên nền đất yếu.

Chương trình nghiên cứu phát triển tương lai:

  • Nâng cấp thuật toán phần mềm TTKGC lên mô hình không gian 3D, tích hợp hiệu ứng xoắn tổng thể và phi tuyến hình học bậc hai ($P-\Delta$).
  • Mở rộng mô hình phi tuyến vật liệu cho thân vách bê tông cốt thép (sử dụng mô hình thớ sợi - Fiber Element Model) để mô phỏng chính xác trạng thái phá hoại cực hạn của toàn hệ thống.
  • Ứng dụng các thuật toán tối ưu hóa thông minh (như Genetic Algorithm, Particle Swarm Optimization hoặc Machine Learning) để tự động hóa hoàn toàn quá trình tìm kiếm thông số độ cứng, lực chảy và vị trí lắp đặt thiết bị tối ưu cho công trình có hình học phức tạp.

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

Nghiên cứu của luận án mang lại giá trị ảnh hưởng sâu rộng trên nhiều phương diện:

  • Tác động học thuật: Cung cấp tài liệu tham khảo chuyên sâu cho các chương trình đào tạo sau đại học ngành Kỹ thuật Xây dựng, Cơ học vật rắn và Kỹ thuật Công trình biển; mở ra các hướng trích dẫn khoa học chất lượng cao về mô hình hóa kết cấu kháng chấn.
  • Tác động công nghiệp và chuyển giao công nghệ: Cung cấp công cụ phần mềm TTKGC gọn nhẹ, độc lập với các hệ phần mềm thương mại đắt tiền, hỗ trợ đắc lực cho các tập đoàn tư vấn thiết kế kết cấu trong việc ứng dụng giải pháp giảm chấn tiên tiến vào các dự án cao ốc tại Việt Nam và khu vực Đông Nam Á.
  • Tác động tiêu chuẩn hóa: Cung cấp cơ sở khoa học và số liệu thực nghiệm số vững chắc phục vụ công tác biên soạn, cập nhật các tiêu chuẩn thiết kế công trình chịu động đất quốc gia theo hướng tiếp cận thiết kế dựa trên năng lượ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 thiết lập ma trận độ cứng siêu phần tử và thuật toán giải phương trình vi phân chuyển động phi tuyến lịch sử thời gian.
  • Kỹ sư kết cấu công trình cao tầng: Nắm bắt quy tắc thiết kế dầm nối lai và chiến lược phân bổ thiết bị TTNL tối ưu theo chiều cao công trình nhằm đạt hiệu quả kinh tế - kỹ thuật cao nhất.
  • Các nhà sản xuất thiết bị giảm chấn: Có căn cứ khoa học chính xác để chế tạo các dòng thiết bị tiêu tán năng lượng có dải lực chảy $Q_y$ và độ cứng $k$ phù hợp với yêu cầu thực tế của công trình nhà nhiều tầng.
  • Chủ đầu tư và Xã hội: Tiết kiệm chi phí đầu tư xây dựng nhờ giảm số lượng thiết bị cần lắp đặt, đồng thời đảm bảo an toàn sinh mạng tuyệt đối và hạn chế tối đa tổn thất kinh tế sau thảm họa động đất.

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ì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc thiết lập thành công ma trận độ cứng giải tích thu gọn cho siêu phần tử dầm nối tích hợp thiết bị TTNL đa trạng thái trễ phi tuyến. Công trình mở rộng trực tiếp lý thuyết phần tử dầm đàn dẻo cổ điển và lý thuyết dầm liên hợp của Harries và Gong, cho phép tích hợp trực tiếp quy luật quan hệ phi tuyến cắt - trượt $Q - \Delta$ của thiết bị vào một phần tử duy nhất mà không cần thông qua các phần tử trung gian zero-length.

2. Đổi mới phương pháp luận của luận án so với các nghiên cứu quốc tế tương tự?

So với các nghiên cứu quốc tế tiên phong như Molham Salameh et al. (2022) (phải sử dụng 3 phần tử thanh và 4 phần tử liên kết zero-length) hay A. Kheyroddin et al. (2022) (sử dụng 2 phần tử thanh và 3 phần tử liên kết zero-length), phương pháp luận của luận án đã gom toàn bộ ứng xử cơ học vào 1 phần tử thanh duy nhất. Sự đổi mới này giúp giảm thiểu số bậc tự do của toàn hệ thống, loại bỏ hiện tượng ma trận độ cứng bị suy biến hoặc khó hội tụ khi giải phi tuyến lịch sử thời gian, đồng thời tăng tốc độ tính toán lên nhiều lần.

3. Phát hiện bất ngờ nhất từ dữ liệu số là gì?

Phát hiện bất ngờ và có giá trị nhất là tính chất bất đối xứng về hiệu quả vị trí: Việc lắp đặt thiết bị TTNL trên toàn bộ 100% số tầng không mang lại hiệu quả vượt trội tuyến tính so với phương án chỉ lắp đặt tại $50%$ đến $66%$ số tầng ở dải trung tâm công trình. Do dầm nối ở các tầng quá thấp hoặc quá cao có biến dạng trượt nhỏ, thiết bị tại đây không phát huy hết khả năng tiêu tán năng lượng trễ, khẳng định rằng chiến lược phân bổ tập trung mang lại hiệu quả tối ưu nhất.

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

Luận án cung cấp đầy đủ quy trình tái lập bao gồm: toàn bộ công thức giải tích ma trận độ cứng phần tử $[K_r]$, thuật toán chuyển đổi tọa độ $[T]$, sơ đồ khối giải thuật tích phân từng bước Newmark, các thông số đầu vào hình học, đặc tính vật liệu bê tông, thép và thông số trễ của các dòng thiết bị (TB1 đến TB20), cùng dữ liệu kiểm chứng chéo chi tiết với SAP2000.

5. Lộ trình nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm mở rộng từ luận án bao gồm 3 cột mốc: (1) Phát triển mô hình siêu phần tử dầm nối không gian 3D có xét biến dạng xoắn và uốn hai phương; (2) Thực hiện chuỗi thí nghiệm bàn rung (Shake Table Test) quy mô lớn trên mô hình tường kép nhiều tầng gắn thiết bị TTNL để kiểm chứng vật lý; (3) Tích hợp thuật toán trí tuệ nhân tạo (AI/Machine Learning) vào phần mềm thiết kế tự động để tối ưu hóa vị trí và thông số thiết bị theo thời gian thực.

Kết luận

Luận án tiến sĩ của tác giả Phạm Thu Hiền đã giải quyết trọn vẹn bài toán cơ học kết cấu phức tạp về tính toán và tối ưu hóa hệ tường kép có gắn thiết bị tiêu tán năng lượng chịu tác động động đất, với 5 kết luận khoa học cốt lõi:

  1. Đã đề xuất thành công mô hình giải tích và thiết lập chính xác ma trận độ cứng cho siêu phần tử dầm nối có gắn thiết bị TTNL với quy luật ứng xử trễ phi tuyến song tuyến tính.
  2. Đã xây dựng hoàn chỉnh chương trình số chuyên dụng TTKGC, được kiểm chứng độ chính xác tuyệt đối qua các lời giải giải tích và mô hình kiểm chứng chéo trên phần mềm SAP2000.
  3. Đã chứng minh một cách định lượng khả năng giảm chấn vượt trội của giải pháp dầm nối gắn thiết bị TTNL, giúp giảm chuyển vị đỉnh từ $35%$ đến hơn $50%$, kiểm soát gia tốc đỉnh và bảo vệ hệ tường bê tông cốt thép luôn làm việc an toàn trong giới hạn đàn hồi.
  4. Đã khám phá quy luật ảnh hưởng của các tham số kỹ thuật thiết bị ($Q_y, k_1, k_2$) và độ cứng dầm thép, xác định dải thông số tối ưu để kích hoạt tối đa năng lượng tiêu tán trễ.
  5. Đã xác lập nguyên lý bố trí thiết bị tối ưu theo chiều cao: Phân bổ tập trung tại $1/2$ đến $2/3$ số tầng ở dải giữa công trình mang lại hiệu năng kháng chấn tương đương việc lắp đặt toàn bộ các tầng, tạo ra bước đột phá về hiệu quả kinh tế - kỹ thuật cho ngành xây dựng công trình cao tầng.