Tổng quan về luận án

Nghiên cứu ứng xử động lực học và khả năng chịu lực của kết cấu công trình xây dựng chịu tải trọng động đất là một bài toán then chốt trong kỹ thuật kết cấu hiện đại. Bối cảnh tự nhiên tại Việt Nam ghi nhận sự tồn tại của nhiều hệ thống đứt gãy kiến tạo còn hoạt động như đới đứt gãy Điện Biên – Mường Lay, Sông Mã – Tuần Giáo – Lai Châu, Cao Bằng – Tiên Yên. Đồng thời, sự phát triển mạnh mẽ của hệ thống các hồ chứa thủy điện quy mô lớn kết hợp với tác động phức tạp của biến đổi khí hậu đã gia tăng đáng kể hiện tượng động đất kích thích và độ nhạy cảm kháng chấn của các vùng dân cư đô thị. Trong thực tiễn thiết kế kháng chấn theo tiêu chuẩn TCVN 9386:2012 (Thiết kế công trình chịu động đất, biên soạn dựa trên nền tảng Tiêu chuẩn Châu Âu Eurocode 8), quy trình tính toán phổ biến hiện hành chủ yếu dựa trên cách tiếp cận tất định (deterministic approach) thông qua phổ thiết kế đàn hồi $S_d(T)$ đã được chiết giảm bởi hệ số ứng xử $q$.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế rằng các tham số động lực học địa chấn và phản ứng kết cấu mang bản chất bất định kép: tính bất định ngẫu nhiên (aleatory randomness) do cơ chế giải phóng năng lượng đứt gãy địa chấn và tính bất định nhận thức hay tính mờ (epistemic uncertainty/fuzziness) do sự thiếu hụt nghiêm trọng các trạm quan trắc gia tốc đồ mạnh tại Việt Nam. Như tác giả luận án đã trích dẫn trực tiếp: "Việt Nam tuy không nằm trong vành đai lửa động đất của thế giới, nhưng vẫn bị ảnh hưởng do tồn tại nhiều đứt gãy như đứt gãy Ðiện Biên - Mường Lay, đứt gãy Sông Mã - Tuần Giáo - Lai Châu…. Đặc biệt là do biến đổi khí hậu và việc xây dựng các công trình thủy điện đang làm gia tăng mức độ nhạy cảm với động đất ở Việt Nam, nhiều trận động đất khá mạnh đã xảy ra". Việc áp dụng các giá trị đơn lẻ mang tính chủ quan của người thiết kế khi nội suy gia tốc nền hay lựa chọn hệ số ứng xử $q$ từ các mô tả định tính trong tiêu chuẩn đã tạo ra những sai lệch tiềm ẩn lớn, làm giảm độ an toàn thực tế của công trình.

Luận án tập trung giải quyết 4 câu hỏi nghiên cứu (Research Questions - RQ) và 4 giả thuyết khoa học (Hypotheses - H) tương ứng:

  • RQ1: Làm thế nào để nhận diện, phân loại và số hóa các quy định mang tính định tính, mờ hóa trong tiêu chuẩn kháng chấn TCVN 9386:2012 thành các đại lượng toán học khả vi?
    • H1: Các tham số mô tả đặc trưng động đất và phản ứng kết cấu trong tiêu chuẩn có thể được số hóa chính xác thông qua việc xây dựng các hàm thuộc (membership functions) $\mu(x)$ và hàm thuộc biến đổi $\mu^*(x)$ trên cơ sở kỹ thuật thực nghiệm.
  • RQ2: Cơ chế tích hợp đồng thời ba loại tham số (tất định, ngẫu nhiên và mờ) trong bài toán phân tích phản ứng kết cấu công trình chịu động đất được thiết lập như thế nào?
    • H2: Thuật toán Monte Carlo cải tiến (MMC) kết hợp biến đổi tương đương định lượng giữa biến mờ và biến ngẫu nhiên cho phép giải quyết triệt để bài toán phân tích động lực học kết cấu trong điều kiện thiếu dữ liệu quan trắc lịch sử.
  • RQ3: Mức độ sai lệch về nội lực và chuyển vị giữa phương pháp tất định hiện hành và phương pháp Monte Carlo cải tiến trên các dạng công trình nhà thấp tầng và cao tầng ở các vùng động đất khác nhau tại Việt Nam là bao nhiêu?
    • H3: Phương pháp tất định hiện hành có xu hướng đánh giá thấp phản lực chân và chuyển vị đỉnh từ 12% đến hơn 30% ở các vùng có nguy cơ động đất từ trung bình đến cao, đặc biệt đối với kết cấu không đều đặn.
  • RQ4: Tương tác phi tuyến của các bộ phận phi kết cấu (như tường chèn) làm thay đổi đáp ứng động lực học của kết cấu khi chịu kích động phổ Monte Carlo cải tiến như thế nào?
    • H4: Việc mô hình hóa phi tuyến có xét đến phần tử tiếp xúc (Gap element) của tường chèn làm tăng đáng kể độ cứng tổng thể ban đầu nhưng làm tập trung ứng suất cắt cục bộ nghiêm trọng tại các chân cột biên.

Khung lý thuyết của công trình được định hình trên sự giao thoa giữa: Lý thuyết tập mờ (Fuzzy Set Theory của Lotfi A. Zadeh, Berd Möller & Michael Beer), Lý thuyết xác suất và mô phỏng ngẫu nhiên (Monte Carlo Simulation của Metropolis & Ulam), cùng Lý thuyết phổ phản ứng động lực học công trình (Response Spectrum Theory của Biot, Housner và Newmark). Đóng góp đột phá của luận án là đã lượng hóa thành công các biên độ bất định trong TCVN 9386:2012, cung cấp một thuật toán phân tích kết cấu có khả năng số hóa hoàn toàn các kịch bản động đất khả dĩ, kiểm soát rủi ro phá hoại dẻo và biến dạng phi tuyến. Phạm vi nghiên cứu bao quát các công trình nhà bê tông cốt thép thấp tầng (Nhà văn hóa 1 tầng vượt nhịp, nhà 4 tầng, 6 tầng, 8 tầng dạng đều đặn và không đều đặn) và nhà cao tầng (các tổ hợp chung cư CT1, CT2, CT3 với quy mô 20, 30 và 40 tầng) phân bố tại các vùng động đất trọng điểm của Việt Nam gồm Mường Lay, Tuần Giáo, Ba Vì, Gia Lâm (Hà Nội), Quận 10 (TP. Hồ Chí Minh), Ngũ Hành Sơn (Đà Nẵng), và Ngô Quyền (Hải Phòng).

Literature Review và Positioning

Lịch sử phát triển các phương pháp đánh giá an toàn kháng chấn ghi nhận sự tiến hóa từ các phương pháp tĩnh lực ngang tương đương (Sano, 1915) đến lý thuyết phổ phản ứng đàn hồi (Biot, 1932; Housner, 1959; Newmark, 1973). Công thức tương quan gia tốc cực đại $a_{max}$ kinh điển của Donovan (1973): $$a_{max} = \frac{1320 e^{0.58M}}{(R + 25)^{1.52}}$$ đã bộc lộ nhiều hạn chế lớn tại vùng gần chấn tiêu ($R \le 15 - 20\text{ km}$), nơi dao động địa chấn thể hiện tính phi tuyến và bất định nghiêm trọng.

Trên trường quốc tế, hai trường phái học thuật lớn đã hình thành nhằm xử lý các bất định trong động lực học công trình:

  1. Trường phái Xác suất - Ngẫu nhiên thuần túy (Probabilistic Simulation Paradigm): Siu-Kui Au & James L. Beck (2001, 2003) phát triển kỹ thuật mô phỏng tập con (Subset Simulation) dựa trên thuật toán chuỗi Markov Monte Carlo (MCMC) và thuật toán Modified Metropolis-Hastings để đánh giá xác suất phá hoại nhỏ trong không gian nhiều chiều của rủi ro địa chấn. Tiếp đó, Carlo Del Gaudio, Paolo Ricci, Gerardo M. Verderame, Gaetano Manfredi (2017) đã mô phỏng Monte Carlo để xây dựng đường cong tổn thương địa chấn quy mô đô thị cho các tòa nhà bê tông cốt thép chịu ảnh hưởng của trận động đất L'Aquila. Gian Felice Giaccu & Luca Caracoglia (2018) đưa ra phương pháp Monte Carlo ngẫu nhiên phân lớp (Layered Stochastic Approximation - LSAMC). Các nghiên cứu này mặc dù giải quyết sâu sắc tính ngẫu nhiên nhưng đòi hỏi một cơ sở dữ liệu thống kê cực lớn về các trận động đất lịch sử – điều kiện mà các quốc gia như Việt Nam không thể đáp ứng.
  2. Trường phái Mờ và Mờ - Ngẫu nhiên kết hợp (Fuzzy & Fuzzy-Randomness Paradigm): Berd Möller & Michael Beer (2004) trong công trình kinh điển Fuzzy Randomness: Uncertainty in Civil Engineering and Computational Mechanics đã chứng minh rằng các tham số kỹ thuật không chắc chắn không chỉ có tính ngẫu nhiên mà còn mang tính mờ do thiếu hụt thông tin. Möller và Beer đã đề xuất phương pháp tối ưu hóa mức $\alpha$ (Alpha-level optimization) để giải bài toán cơ học kết cấu mờ. Tiếp cận này được mở rộng bởi D. Moens & D. Vandepitte (2005) trong phân tích phần tử hữu hạn mờ (Fuzzy Finite Element Analysis).

Tại Việt Nam, các công trình của Lê Xuân Huỳnh & Nguyễn Hùng Tuấn (2016) về phương pháp phần tử hữu hạn tham số mờ, Nguyễn Đình Xân (2005), Chu Thanh Bình (2014), Lê Công Duy (2014), Nguyễn Thanh Hưng (2017) đã bước đầu đặt nền móng cho việc ứng dụng lý thuyết tập mờ và độ tin cậy trong cơ học xây dựng. Tuy nhiên, phần lớn các nghiên cứu trước đây chỉ giải quyết bài toán cơ học lý thuyết thuần túy trên các kết cấu thanh đơn lẻ và giả định toàn bộ tham số đầu vào đều là biến mờ, khiến khối lượng tính toán bùng nổ theo cấp số nhân và không thể áp dụng vào quy trình thiết kế công trình tổng thể theo tiêu chuẩn hiện hành.

                    ┌─────────────────────────────────────────────────────────┐
                    │      CÁC TIẾP CẬN BẤT ĐỊNH TRONG THIẾT KẾ KHÁNG CHẤN    │
                    └────────────────────────────┬────────────────────────────┘
                                                 │
            ┌────────────────────────────────────┴───────────────────────────────────┐
            ▼                                                                        ▼
┌──────────────────────────────────────┐                   ┌──────────────────────────────────────┐
│ TRƯỜNG PHÁI XÁC SUẤT THUẦN TÚY       │                   │ TRƯỜNG PHÁI MỜ - NGẪU NHIÊN          │
│ (Au & Beck 2003; Del Gaudio 2017)    │                   │ (Möller & Beer 2004; Zadeh 1965)     │
├──────────────────────────────────────┤                   ├──────────────────────────────────────┤
│ - Giả định biết rõ hàm mật độ xác suất│                   │ - Xử lý thông tin thiếu hụt/bất định │
│ - Đòi hỏi kho dữ liệu gia tốc kế lớn │                   │ - Tối ưu mức Alpha phức tạp          │
│ - Không xử lý được ngôn ngữ định tính│                   │ - Chưa gắn với tiêu chuẩn thực hành  │
└──────────────────┬───────────────────┘                   └──────────────────┬───────────────────┘
                   │                                                          │
                   └─────────────────────────────┬────────────────────────────┘
                                                 │
                                                 ▼
                    ┌─────────────────────────────────────────────────────────┐
                    │     VỊ TRÍ ĐỘT PHÁ CỦA LUẬN ÁN (VŨ TRỌNG HUY, 2021)     │
                    ├─────────────────────────────────────────────────────────┤
                    │ • Tích hợp 3 loại biến: Tất định + Ngẫu nhiên + Mờ      │
                    │ • Số hóa 6 nhóm tham số mờ trong TCVN 9386:2012         │
                    │ • Phương pháp Monte Carlo cải tiến (MMC) xử lý trọng số │
                    │ • Tương thích hoàn toàn với quy trình mô phỏng ETABS v18│
                    └─────────────────────────────────────────────────────────┘

Luận án của Vũ Trọng Huy định vị chính xác tại điểm nghẽn của hai dòng nghiên cứu: xây dựng cầu nối phương pháp luận giữa lý thuyết mờ ngẫu nhiên phức tạp với tiêu chuẩn thực hành kỹ thuật TCVN 9386:2012, thông qua phương pháp Monte Carlo cải tiến có khả năng xử lý đồng thời cả 3 dạng biến: tất định, ngẫu nhiên và mờ. So sánh với các tiêu chuẩn quốc tế như ASCE 7-10/IBC 2012 của Mỹ, BSL 2000 của Nhật Bản, SP 14.13330 của Nga và GB 50011 của Trung Quốc, nghiên cứu này chứng minh rằng việc số hóa các điều khoản ngôn ngữ định tính thành các hàm mờ khả dĩ là phương pháp tối ưu nhất để loại bỏ tính chủ quan trong thiết kế kháng chấn ở các quốc gia đang phát triển.

Đó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 một bước tiến lý thuyết đáng kể trong việc mở rộng cơ học kết cấu mờ (Fuzzy Structural Mechanics) sang lĩnh vực kỹ thuật kháng chấn thực hành. Tác giả luận án khẳng định rõ nét trong luận điểm nghiên cứu: "Công cụ xác xuất thống kê không thể giải các bài toán có biến mờ như bài toán tìm đáp ứng của kết cấu công trình nhà chịu động đất, nên việc nghiên cứu ứng dụng phương pháp Monte - Carlo cải tiến để giải các bài toán có đồng thời các loại biến: tiền định, ngẫu nhiên và mờ là cần thiết, có ý nghĩa khoa học".

                               ┌────────────────────────────────────────┐
                               │   PHÂN LOẠI 3 THÀNH PHẦN THAM SỐ       │
                               │   TRONG PHÂN TÍCH KHÁNG CHẤN CÔNG TRÌNH│
                               └──────────────────┬─────────────────────┘
                                                  │
         ┌────────────────────────────────────────┼────────────────────────────────────────┐
         ▼                                        ▼                                        ▼
┌─────────────────────────┐              ┌─────────────────────────┐              ┌─────────────────────────┐
│   THAM SỐ TẤT ĐỊNH      │              │   THAM SỐ NGẪU NHIÊN    │              │      THAM SỐ MỜ         │
├─────────────────────────┤              ├─────────────────────────┤              ├─────────────────────────┤
│ - Kích thước hình học   │              │ - Mô đun đàn hồi ($E$)   │              │ - Gia tốc nền ($a_{gR}$) │
│   tiết diện ($b, h, t$) │              │ - Tải trọng tĩnh ($G$)  │              │ - Hệ số ứng xử ($q$)    │
│ - Sơ đồ hình học khung  │              │ - Hoạt tải sử dụng ($Q$)│              │ - Hệ số nền ($S$)       │
│ - Cường độ danh định    │              │ - Tải trọng gió tĩnh    │              │ - Chu kỳ $T_B, T_C, T_D$│
└─────────────────────────┘              └─────────────────────────┘              └─────────────────────────┘

Mô hình lý thuyết mới đề xuất phân rã bài toán kết cấu chịu động đất thành một hệ thống toán học thống nhất với 3 thành phần tham số:

  1. Tham số tất định ($X_D$): Kích thước hình học của cấu kiện dầm, cột, vách, sàn ($b, h, L$), điều kiện biên liên kết cứng hoặc khớp.
  2. Tham số ngẫu nhiên ($X_R$): Mô đun đàn hồi của bê tông cốt thép $E$, dung trọng vật liệu, tĩnh tải hoàn thiện $G$, hoạt tải sử dụng $Q$, tuân theo các phân phối xác suất liên tục (Gaussian, Lognormal).
  3. Tham số mờ ($X_F$): Đỉnh gia tốc nền tham chiếu $a_{gR}$, hệ số tầm quan trọng công trình $\gamma_I$, hệ số ứng xử kết cấu $q$, hệ số nền $S$, các chu kỳ đặc trưng của phổ phản ứng $T_B, T_C, T_D$, được biểu diễn bằng các tập mờ với các hàm thuộc $\mu(x) \in [0, 1]$.

Mô hình lý thuyết thiết lập quan hệ chuyển đổi tương đương định lượng giữa biến mờ thứ cấp $P$ và biến ngẫu nhiên $E$ thông qua việc chuẩn hóa diện tích hàm thuộc: $$\mu^*(x) = \frac{\mu(x)}{\int_{-\infty}^{+\infty} \mu(x) dx}$$ từ đó chuyển bài toán phân tích mờ phức tạp về bài toán mô phỏng thống kê có trọng số tích hợp, khắc phục được sự bùng nổ tính toán của các phương pháp giải mờ truyền thống.

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

Khung phân tích của luận án tích hợp liền mạch ba trụ cột khoa học: Lý thuyết tập mờ, Phổ phản ứng gia tốc thiết kế $S_d(T)$ của TCVN 9386:2012/Eurocode 8, và Phương pháp Monte Carlo cải tiến (MMC).

                                  ┌──────────────────────────────────────────────┐
                                  │ BƯỚC 1: SỐ HÓA THAM SỐ MỜ THEO TCVN 9386     │
                                  │ (Xác định hàm thuộc tam giác, hình thang...) │
                                  └──────────────────────┬───────────────────────┘
                                                         │
                                                         ▼
                                  ┌──────────────────────────────────────────────┐
                                  │ BƯỚC 2: BIẾN ĐỔI TƯƠNG ĐƯƠNG ĐỊNH LƯỢNG      │
                                  │ (Chuẩn hóa hàm thuộc gốc thành hàm mật độ)   │
                                  └──────────────────────┬───────────────────────┘
                                                         │
                                                         ▼
                                  ┌──────────────────────────────────────────────┐
                                  │ BƯỚC 3: PHÁT SINH MẪU TOÀN THỂ MONTE CARLO   │
                                  │ (Tổ hợp các kịch bản phổ khả dĩ Sk(T))       │
                                  └──────────────────────┬───────────────────────┘
                                                         │
                                                         ▼
                                  ┌──────────────────────────────────────────────┐
                                  │ BƯỚC 4: GIẢI BÀI TOÁN KẾT CẤU TRÊN ETABS V18 │
                                  │ (Phân tích phổ đàn hồi & phi tuyến theo time)│
                                  └──────────────────────┬───────────────────────┘
                                                         │
                                                         ▼
                                  ┌──────────────────────────────────────────────┐
                                  │ BƯỚC 5: GIẢI MỜ NỘI LỰC & CHUYỂN VỊ          │
                                  │ (Tổng hợp thống kê có trọng số xác suất)     │
                                  └──────────────────────────────────────────────┘

Sáu loại đại lượng mờ cốt lõi trong tiêu chuẩn được phân loại và số hóa thành các dạng hàm thuộc hình học:

  • Đại lượng mờ loại 1 (Khoảng giá trị đối xứng/hình thang): Áp dụng cho hệ số nền $S$ và các chu kỳ $T_B, T_C, T_D$ của nền đất loại A, B, C, D, E.
  • Đại lượng mờ loại 2 (Giá trị tin tưởng duy nhất hoặc khoảng tin tưởng hữu hạn): Áp dụng cho đỉnh gia tốc nền tham chiếu $a_{gR}$ theo đường đẳng trị tại các địa danh cụ thể (ví dụ: quận Ba Đình $a_{gR} \in [0.0984g; 0.1158g]$ với giá trị nội suy $x_0 = 0.1064g$).
  • Đại lượng mờ loại 3 và loại 4 (Hàm thuộc bán vô hạn và hình chuông): Áp dụng cho hệ số ứng xử $q$ khi xét đến tính đều đặn trên mặt bằng và mặt đứng ($q = q_0 \cdot k_w \ge 1.5$).

Quy trình giải mờ được thực hiện thông qua việc gán trọng số khả dĩ $w_k$ cho từng trường hợp phổ thiết kế thứ $k$: $$w_k = \prod_{j=1}^{m} \mu^*j(x{j,k})$$ $$\bar{R} = \frac{\sum_{k=1}^{N} w_k \cdot R_k}{\sum_{k=1}^{N} w_k}$$ trong đó $R_k$ là đáp ứng kết cấu (lực cắt đáy $F_b$, mô men uốn $M$, chuyển vị tầng $d_r$) ứng với kịch bản phổ thứ $k$.

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 lập trường thực chứng kỹ thuật kết hợp thực dụng luận tính toán (Computational Engineering Pragmatism). Thiết kế nghiên cứu triển khai theo mô hình đa tầng bậc (multi-level design), từ cấp độ ứng xử vật liệu, phần tử cấu kiện (dầm, cột, vách bê tông cốt thép, thanh chống tương đương của tường chèn) đến toàn bộ không gian 3D của hệ kết cấu chịu lực.

┌────────────────────────────────────────────────────────────────────────────────────────────────────────┐
│                                    THIẾT KẾ ĐỐI CHỨNG NGHIÊN CỨU                                       │
├──────────────────────────────────┬──────────────────────────────────┬──────────────────────────────────┤
│ PHƯƠNG PHÁP TẤT ĐỊNH HIỆN HÀNH   │ PHƯƠNG PHÁP MONTE CARLO CẢI TIẾN │ PHÂN TÍCH ĐỘNG PHI TUYẾN         │
│ (Deterministic TCVN 9386:2012)   │ (Modified Monte Carlo - MMC)     │ (Nonlinear Time-History + Gap)   │
├──────────────────────────────────┼──────────────────────────────────┼──────────────────────────────────┤
│ - Sử dụng giá trị danh định      │ - Số hóa 6 nhóm tham số mờ       │ - Dữ liệu sóng động đất thực tế  │
│ - 1 kịch bản phổ thiết kế duy nhất│ - Tổ hợp $N$ kịch bản phổ khả dĩ │ - Khớp dẻo phân tán (Plasticity) │
│ - Đánh giá tuyến tính đàn hồi    │ - Giải mờ thống kê có trọng số   │ - Mô hình phần tử tiếp xúc Gap   │
└──────────────────────────────────┴──────────────────────────────────┴──────────────────────────────────┘

Quy trình nghiên cứu rigorous

Quy trình thực nghiệm số và kiểm định được tiến hành qua các giao thức nghiêm ngặt:

  1. Giao thức lấy mẫu và phát sinh số giả ngẫu nhiên: Sử dụng thuật toán sinh số ngẫu nhiên trên đoạn $[0, 1]$ được chuẩn hóa, quét toàn bộ không gian biến thiên của các tham số mờ theo phương pháp mẫu toàn thể (Full domain sampling).
  2. Quy trình lập mô hình phần tử hữu hạn: Toàn bộ kết cấu được mô hình hóa không gian 3D trên phần mềm chuyên dụng ETABS v18, tuân thủ các điều kiện biên về khối lượng tham gia dao động ($> 90%$), độ lệch tâm ngẫu nhiên tầng $e_{ai} = \pm 0.05 L_i$.
  3. Tam giác đạc phương pháp (Methodological Triangulation): Kết quả phân tích được đối chứng chéo giữa ba phương pháp độc lập:
    • Phân tích phổ phản ứng đàn hồi theo phương pháp tất định của TCVN 9386:2012.
    • Phân tích phổ phản ứng bằng phương pháp Monte Carlo cải tiến (MMC).
    • Phân tích động lực học phi tuyến theo lịch sử thời gian (Nonlinear Dynamic Time-History Analysis) có xét đến hiệu ứng tương tác của tường gạch xây thông qua mô hình thanh chống chéo tương đương và phần tử tiếp xúc Gap.

Data và phân tích

Tập dữ liệu khảo sát bao gồm hai nhóm công trình đại diện lớn:

Nhóm công trình Mô hình chi tiết Vị trí địa lý / Gia tốc nền tham chiếu $a_{gR}$ Loại đất nền / Tham số $S, T_B, T_C, T_D$ Hệ số ứng xử $q$
Nhà thấp tầng - Nhà văn hóa (1 tầng, nhịp lớn)
- Khung 4 tầng, 6 tầng, 8 tầng (đều đặn và không đều đặn)		 - Mường Lay ($a_g = 0.150g$)
- Tuần Giáo ($a_g = 0.135g$)
- Ba Vì ($a_g = 0.086g$)		 Nền loại C ($S=1.15, T_B=0.20s, T_C=0.60s, T_D=2.0s$) và Nền loại D ($S=1.35$) $q \in [1.5; 3.9]$ (Cấp dẻo DCL, DCM)
Nhà cao tầng - Tổ hợp CT1 (20 tầng)
- Tổ hợp CT2 (20 tầng)
- Tổ hợp CT3 (20, 30, 40 tầng)		 - Gia Lâm - Hà Nội ($a_g = 0.106g$)
- Q.10 - TP.HCM ($a_g = 0.065g$)
- Ngũ Hành Sơn - ĐN ($a_g = 0.054g$)
- Ngô Quyền - HP ($a_g = 0.082g$)		 Nền loại D ($S=1.35, T_B=0.20s, T_C=0.80s, T_D=2.0s$) và Nền loại E ($S=1.40$) $q \in [2.5; 3.6]$ (Hệ kết cấu hỗn hợp vách - khung)

Phương pháp kiểm tra tính vững (robustness checks) được thực hiện bằng cách thay đổi hình dạng hàm thuộc (từ tam giác sang hình thang và Gauss) để đánh giá độ nhạy của phản lực chân và chuyển vị đỉnh, bảo đảm kết quả hội tụ với sai số thống kê dưới 2.5%.

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

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

Nghiên cứu mang lại 4 phát hiện đột phá với đầy đủ bằng chứng định lượng từ kết quả mô phỏng số:

                            TỶ LỆ SAI LỆCH PHẢN LỰC CHÂN CÔNG TRÌNH (MMC vs TCVN)
  Địa điểm / Mô hình              [% Sai lệch phản lực chân]
  ─────────────────────────────────────────────────────────────────────────────
  Nhà văn hóa (Mường Lay)         ██████████████████████████ 25.8%
  Nhà văn hóa (Tuần Giáo)         █████████████████████ 21.4%
  Nhà văn hóa (Ba Vì)             ████████████ 12.8%
  Khung 8T không đều đặn          ████████████████████████████████ 32.4%
  Chung cư CT1 20T (Gia Lâm)      ████████████████████ 20.3%
  Chung cư CT3 40T (Gia Lâm)      ████████████████████████ 24.1%
  ─────────────────────────────────────────────────────────────────────────────
  1. Sự đánh giá thấp hệ thống của phương pháp tất định: Tại các vùng có cường độ địa chấn mạnh như Mường Lay ($a_g = 0.150g$) và Tuần Giáo ($a_g = 0.135g$), phản lực chân cột và lực cắt đáy tính toán theo phương pháp Monte Carlo cải tiến lớn hơn phương pháp tất định hiện hành từ $21.4%$ đến $25.8%$. Tại Ba Vì ($a_g = 0.086g$), mức chênh lệch đạt $12.8%$. Điều này chỉ ra rằng việc bỏ qua tính mờ của tham số gia tốc nền và hệ số nền $S$ trong TCVN 9386:2012 đã dẫn đến việc thiết kế thiếu an toàn đối với các công trình công cộng tập trung đông người.
  2. Khuếch đại phản ứng phi tuyến trên kết cấu không đều đặn: Đối với kết cấu nhà thấp tầng có sự gián đoạn độ cứng trên chiều cao (không đều đặn), độ lệch phản lực chân giữa hai phương pháp tăng vọt lên mức $32.4%$ ở công trình 8 tầng, so với chỉ $18.6%$ ở công trình 8 tầng đều đặn. Sự biến thiên mờ của hệ số ứng xử $q$ làm lộ rõ các vị trí tập trung biến dạng dẻo mà phương pháp tất định không thể phát hiện.
  3. Hiệu ứng chu kỳ dao động bậc cao ở công trình cao tầng (CT3 20T, 30T, 40T): Khi chiều cao công trình tăng từ 20 tầng ($H = 65\text{ m}$) lên 40 tầng ($H = 130\text{ m}$), tỷ lệ lệch chuyển vị đỉnh ngang giữa MMC và TCVN tăng từ $15.2%$ lên $24.1%$ tại khu vực Gia Lâm – Hà Nội. Phổ Monte Carlo cải tiến kích hoạt mạnh mẽ các dạng dao động bậc cao ($T_2, T_3$), làm gia tăng lực cắt ở các tầng phía trên đỉnh tháp vượt quá $22.5%$ so với dự báo của tiêu chuẩn hiện hành.
  4. Ảnh hưởng hai mặt của tường chèn trong phân tích phi tuyến: Kết quả phân tích lịch sử thời gian có xét phần tử Gap cho thấy: tường chèn gạch làm tăng độ cứng tổng thể của công trình, kéo chu kỳ dao động cơ bản $T_1$ giảm từ $0.82\text{s}$ xuống $0.48\text{s}$, làm giảm chuyển vị đỉnh công trình khoảng $35% - 42%$. Tuy nhiên, tường chèn làm tăng lực cắt đáy cục bộ lên $1.45$ lần và gây ra hiện tượng phá hoại cắt giòn tại các đầu cột ngắn (short-column effect).

Implications đa chiều

  • Về mặt lý thuyết: Luận án chứng minh sự cần thiết phải chuyển dịch mô hình tính toán kháng chấn từ tất định sang mô hình bất định mờ – ngẫu nhiên, hoàn thiện lý thuyết động lực học công trình trong điều kiện dữ liệu mờ.
  • Về phương pháp luận: Thiết lập quy trình chuẩn hóa số hóa tham số mờ, mở ra khả năng tích hợp linh hoạt thuật toán Monte Carlo cải tiến vào các phần mềm kết cấu thương mại (ETABS, SAP2000) mà không cần can thiệp vào mã nguồn lõi.
  • Về thực tiễn thiết kế: Cung cấp bộ hệ số hiệu chỉnh an toàn cho kỹ sư khi lựa chọn $q$ và $a_{gR}$ theo TCVN 9386:2012, bảo đảm kích thước tiết diện và hàm lượng cốt thép đủ khả năng tiêu tán năng lượng trong các trận động đất vượt cấp thiết kế.
  • Về chính sách và tiêu chuẩn: Cung cấp cơ sở khoa học định lượng vững chắc cho Viện Khoa học Công nghệ Xây dựng (IBST) và Bộ Xây dựng trong quá trình soát xét, bổ sung và cập nhật tiêu chuẩn quốc gia TCVN 9386 trong giai đoạn tới.

Limitations và Future Research

Nghiên cứu thẳng thắn thừa nhận 4 giới hạn khoa học cụ thể:

  1. Giả thiết độc lập thống kê: Các biến mờ đầu vào ($a_{gR}, S, q, T_B, T_C$) được giả định độc lập với nhau, trong khi trên thực tế có thể tồn tại mối tương quan phụ thuộc phức tạp giữa loại nền đất $S$ và các chu kỳ đặc trưng $T_B, T_C$.
  2. Tương tác đất – kết cấu (Soil-Structure Interaction - SSI): Mô hình nghiên cứu giả định ngàm cứng tại chân cột móng, chưa xét đến hiệu ứng suy giảm độ cứng và cản bức xạ của nền đất mềm dưới móng bè/móng cọc.
  3. Phạm vi vật liệu: Nghiên cứu chỉ tập trung vào hệ kết cấu bê tông cốt thép toàn khối, chưa mở rộng cho các hệ kết cấu thép tiền chế, kết cấu liên hợp thép – bê tông hoặc kết cấu gỗ hiện đại.
  4. Cơ sở dữ liệu gia tốc đồ: Số lượng băng ghi gia tốc động đất mạnh thực tế tại Việt Nam còn hạn chế, các phân tích phi tuyến thời gian vẫn phải sử dụng một phần dữ liệu chuẩn hóa quốc tế (như sóng El Centro 1940).

Chương trình nghiên cứu tiếp theo (Future Research Agenda) mở ra 4 hướng đi chiến lược:

  • Xây dựng mô hình Copula mờ nhiều chiều (Fuzzy Multidimensional Copula) để mô tả chính xác sự tương quan phụ thuộc giữa các tham số địa chấn mờ.
  • Ứng dụng trí tuệ nhân tạo (Mạng nơ-ron sâu/Physics-Informed Neural Networks - PINNs) để tự động hóa quá trình tối ưu hóa việc sinh số giả ngẫu nhiên có trọng số trong thuật toán MMC.
  • Mở rộng phân tích độ tin cậy mờ cho các công trình sử dụng công nghệ giảm chấn hiện đại (Passive/Active Dampers) và cách chấn đáy (Base Isolation).
  • Nghiên cứu thực nghiệm trên bàn rung động học (Shaking table test) cho mô hình khung có xét tường chèn chịu kích động động đất theo các kịch bản phổ mờ.

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

  • Ảnh hưởng học thuật (Academic Impact): Luận án tạo lập một dấu mốc mới trong việc ứng dụng toán học mờ vào kỹ thuật công trình tại Việt Nam. Các bài báo công bố trên Tạp chí Xây dựng (2019, 2021) và các kỷ yếu hội nghị chuyên ngành đã nhận được sự quan tâm và trích dẫn sâu rộng từ cộng đồng nghiên cứu độ tin cậy kết cấu.
  • Chuyển đổi công nghiệp xây dựng (Industry Transformation): Thay đổi tư duy thiết kế kháng chấn của các tổng công ty tư vấn thiết kế lớn (VNCC, CDC, Coninco). Cung cấp giải pháp đánh giá rủi ro kết cấu tin cậy hơn cho các dự án chung cư cao tầng, bệnh viện, trường học tại Hà Nội, TP. Hồ Chí Minh và khu vực Tây Bắc.
  • Tác động chính sách công (Policy Influence): Là tài liệu tham khảo cốt lõi cho các hội đồng khoa học chuyên ngành thuộc Bộ Xây dựng trong việc hoàn thiện hệ thống tiêu chuẩn kỹ thuật quốc gia về an toàn công trình trước thiên tai động đất.
  • Lợi ích xã hội (Societal Benefits): Góp phần bảo vệ tính mạng con người và giảm thiểu thiệt hại hàng nghìn tỷ đồng cho cơ sở hạ tầng đô thị khi xảy ra các biến cố động đất cực đoan trong tương lai.

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

                               ┌────────────────────────────────────────────────────────┐
                               │                 ĐỐI TƯỢNG HƯỞNG LỢI                    │
                               └──────────────────────────┬─────────────────────────────┘
                                                          │
         ┌────────────────────────┬───────────────────────┴───────────────┬────────────────────────┐
         ▼                        ▼                                       ▼                        ▼
┌──────────────────┐    ┌──────────────────┐                    ┌──────────────────┐    ┌──────────────────┐
│ NGHIÊN CỨU SINH  │    │ KỸ SƯ KẾT CẤU &  │                    │ CƠ QUAN QUẢN LÝ  │    │ CHỦ ĐẦU TƯ &     │
│ & NHÀ KHOA HỌC   │    │ CÔNG TY TƯ VẤN   │                    │ & SOẠN THẢO TCVN │    │ CỘNG ĐỒNG XÃ HỘI │
├──────────────────┤    ├──────────────────┤                    ├──────────────────┤    ├──────────────────┤
│ Kế thừa khung    │    │ Tiếp cận quy     │                    │ Sở hữu bằng chứng│    │ Công trình an    │
│ toán học mờ &    │    │ trình mô phỏng   │                    │ định lượng để    │    │ toàn, bền vững;  │
│ thuật toán MMC   │    │ ETABS chính xác, │                    │ soát xét tiêu    │    │ giảm thiểu rủi   │
│ cho đề tài mới   │    │ loại bỏ rủi ro   │                    │ chuẩn kháng chấn │    │ ro sinh mạng     │
└──────────────────┘    └──────────────────┘                    └──────────────────┘    └──────────────────┘
  1. Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phân tích toán học hoàn chỉnh kết hợp giữa lý thuyết tập mờ và mô phỏng ngẫu nhiên; mở rộng nghiên cứu sang các loại tải trọng bất định khác như gió bão cực đoan, sóng thần, va chạm nổ.
  2. Kỹ sư kết cấu và Doanh nghiệp tư vấn thiết kế: Nắm vững công cụ số hóa để đánh giá định lượng độ an toàn thực tế của công trình; tránh được các sai số chủ quan nguy hiểm khi tra cứu các bảng biểu định tính trong TCVN 9386:2012.
  3. Cơ quan quản lý Nhà nước và Ban soạn thảo tiêu chuẩn: Sở hữu cơ sở dữ liệu thực nghiệm số vững chắc để điều chỉnh phân vùng gia tốc nền $a_{gR}$ và biên độ hệ số ứng xử $q$ trong các bản cập nhật tiêu chuẩn tiếp theo.
  4. Chủ đầu tư và Cộng đồng xã hội: Bảo đảm tối ưu hóa chi phí vật liệu xây dựng song hành với việc nâng cao năng lực thích ứng và khả năng chống chịu thảm họa của công trình đô thị.

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 lý thuyết độc đáo nhất là việc thiết lập mô hình chuyển đổi tương đương định lượng giữa biến mờ thứ cấp và biến ngẫu nhiên, tích hợp đồng thời 3 loại biến (tất định, ngẫu nhiên, mờ) trong bài toán động lực học công trình. Công trình đã mở rộng trực tiếp Lý thuyết tập mờ của Lotfi A. Zadeh (1965) và Lý thuyết ngẫu nhiên mờ (Fuzzy Randomness) của Möller & Beer (2004) vào quy trình thiết kế kháng chấn thực hành theo TCVN 9386:2012.

2. Phương pháp Monte Carlo cải tiến (MMC) trong luận án có điểm gì vượt trội so với các nghiên cứu tiền nhiệm?

So với phương pháp Subset Simulation của Au & Beck (2003) (chỉ xử lý biến ngẫu nhiên qua chuỗi Markov) và phương pháp phần tử hữu hạn mờ của Möller & Beer (2004) hay Lê Xuân Huỳnh (2016) (đòi hỏi giải tối ưu hóa mức $\alpha$ phi tuyến phức tạp), phương pháp MMC của luận án sử dụng kỹ thuật sinh mẫu toàn thể kết hợp gán trọng số xác suất tích hợp $\prod \mu^*(x)$. Phương pháp này không làm bùng nổ số lượng phép tính, cho phép tận dụng trực tiếp các phần tử hữu hạn đàn hồi/phi tuyến thương mại (như ETABS v18) để giải mờ nội lực cho toàn bộ hệ kết cấu không gian thực tế.

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

Phát hiện bất ngờ nhất là mức độ thiếu an toàn danh định nghiêm trọng của tiêu chuẩn hiện hành đối với kết cấu không đều đặn trên mặt đứng. Độ lệch phản lực chân giữa MMC và phương pháp tất định của TCVN 9386:2012 lên tới $32.4%$ ở công trình 8 tầng không đều đặn, cao gấp $1.74$ lần so với công trình đều đặn ($18.6%$). Điều này chứng minh rằng quy định chọn hệ số ứng xử $q_0$ trơn trong tiêu chuẩn đã đánh giá thấp đáng kể sự tập trung ứng suất dẻo tại các tầng gián đoạn độ cứng.

4. Luận án có cung cấp giao thức tái lập (replication protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án trình bày chi tiết từng bước: từ bảng số giả ngẫu nhiên trên đoạn $[0, 1]$, bảng trọng số giải mờ cho biến $P$ và $E$, công thức giải tích xác định hàm thuộc $\mu(x)$ của 6 nhóm tham số mờ, đến sơ đồ khối tích hợp và các bước nhập dữ liệu phổ khả dĩ trực tiếp vào bảng định nghĩa phổ phản ứng (Response Spectrum Functions) và tổ hợp tải trọng trọng số (Load Combinations) trong ETABS v18.

5. Chương trình nghị sự nghiên cứu 10 năm được vạch ra như thế nào?

Chương trình nghiên cứu 10 năm tập trung vào 3 trụ cột: (1) Mở rộng thuật toán MMC cho bài toán tương tác phi tuyến 3D đất – móng – kết cấu (SSI); (2) Tích hợp học máy (Machine Learning) và dữ liệu lớn (Big Data) quan trắc thời gian thực từ mạng lưới cảm biến IoT trên công trình để tự động cập nhật hàm thuộc $\mu(x)$; (3) Xây dựng bộ tiêu chuẩn thiết kế kháng chấn dựa trên hiệu năng mờ (Fuzzy Performance-Based Seismic Design Code) thế hệ mới cho Việt Nam.

Kết luận

  1. Luận án đã giải quyết thành công bài toán cốt lõi về nhận diện, phân loại và số hóa hoàn chỉnh 6 nhóm tham số mờ trong tiêu chuẩn thiết kế kháng chấn TCVN 9386:2012, bao gồm gia tốc nền $a_{gR}$, hệ số tầm quan trọng $\gamma_I$, hệ số nền $S$, các chu kỳ giới hạn $T_B, T_C, T_D$ và hệ số ứng xử $q$.
  2. Hoàn thiện và ứng dụng thành công thuật toán Monte Carlo cải tiến (MMC) có khả năng tích hợp đồng thời ba nhóm tham số (tất định, ngẫu nhiên và mờ), giải quyết triệt để bài toán phân tích động lực học kết cấu trong bối cảnh thiếu hụt dữ liệu địa chấn lịch sử tại Việt Nam.
  3. Cung cấp bằng chứng thực nghiệm số toàn diện trên các công trình nhà thấp tầng và cao tầng (20, 30, 40 tầng) tại các đô thị trọng điểm, chỉ ra mức sai lệch nội lực và chuyển vị từ $12.8%$ đến $32.4%$ giữa phương pháp MMC và phương pháp tất định hiện hành, chứng minh sự cần thiết phải kiểm soát biên độ bất định trong thiết kế kháng chấn.
  4. Làm sáng tỏ ứng xử phi tuyến phức tạp của kết cấu có xét tương tác tường chèn thông qua mô hình phần tử tiếp xúc Gap, chỉ ra nguy cơ hình thành hiệu ứng cột ngắn và phá hoại cắt giòn cục bộ dưới tác động của các kịch bản động đất mờ.
  5. Tạo ra bước chuyển dịch mô hình nhận thức (paradigm shift) từ thiết kế kháng chấn tất định truyền thống sang phân tích bất định mờ – ngẫu nhiên tiên tiến, mở ra 3 nhánh nghiên cứu chuyên sâu về cơ học mờ, trí tuệ nhân tạo trong động lực học công trình và thiết kế kháng chấn dựa trên độ tin cậy hiệu năng.
  6. Thiết lập một công trình khoa học mẫu mực có giá trị kế thừa lâu dài, đóng góp trực tiếp vào sự nghiệp xây dựng tiêu chuẩn kỹ thuật quốc gia và bảo đảm an toàn bền vững cho các công trình xây dựng trước hiểm họa thiên tai động đất tại Việt Nam.