Nghiên cứu sự làm việc của kết cấu nhà xây dựng ở Việt Nam chịu ảnh hưởng động đất bằng Monte Carlo cải tiến

Tổng quan về luận án

Nghiên cứu về độ an toàn và khả năng kháng chấn của các công trình xây dựng trong vùng chịu ảnh hưởng địa chấn là một trong những bài toán cốt lõi của cơ học kết cấu hiện đại. Tại Việt Nam, mặc dù không nằm trọn vẹn trong vành đai lửa Thái Bình Dương, hiểm họa địa chấn vẫn hiện hữu rõ rệt do sự đan xen của các đới đứt gãy kiến tạo đang 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, cùng với tác động gia tăng từ các công trình thủy điện tích nước và biến đổi địa chất cục bộ. 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 tác giả Vũ Trọng Huy, dưới sự hướng dẫn khoa học của PGS.TS. Nguyễn Võ Thông và GS.TS. Nguyễn Văn Phó tại Viện Khoa học Công nghệ Xây dựng (IBST), đã giải quyết bài toán kháng chấn phức tạp này bằng việc tiếp cận và phát triển phương pháp số tiên tiến: Phương pháp Monte - Carlo cải tiến (Improved Monte Carlo Method - MMC).

"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."

Khoảng trống nghiên cứu (research gap) mang tính căn bản mà luận án chỉ ra nằm ở sự mâu thuẫn giữa bản chất bất định của các đại lượng kháng chấn và phương pháp tính toán tất định (deterministic approach) hiện hành. Tiêu chuẩn thiết kế kháng chấn Việt Nam TCVN 9386:2012 (biên soạn dựa trên nền tảng Eurocode 8 - EC8) chứa đựng nhiều quy định mang tính định tính, ngôn ngữ diễn đạt ước lệ và biên độ dao động rộng của các tham số đầu vào. Trong thực hành kỹ sư, các tham số như hệ số ứng xử kết cấu ($q$), đỉnh gia tốc nền tham chiếu ($a_{gR}$), hệ số đất nền ($S$) hay các chu kỳ góc của phổ phản ứng ($T_B, T_C, T_D$) thường bị gán một giá trị đơn nhất dựa trên đánh giá chủ quan của người thiết kế. Các phương pháp xác suất thống kê cổ điển hay mô phỏng Monte - Carlo truyền thống chỉ xử lý được các biến ngẫu nhiên thuần túy (aleatory uncertainty), hoàn toàn bất lực trước tính mờ và độ mơ hồ thông tin (epistemic fuzziness).

Luận án đặt ra hệ thống câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) cụ thể:

• RQ1: Làm thế nào để định lượng hóa, mờ hóa và số hóa các tham số mang tính bất định ngôn ngữ trong tiêu chuẩn TCVN 9386:2012 thành các hàm thuộc toán học?
• RQ2: Cơ chế tích hợp đồng thời ba nhóm tham số: tất định, ngẫu nhiên và mờ vào cùng một thuật toán mô phỏng kết cấu tổng thể là gì?
• 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 đối với các loại hình công trình thực tế (thấp tầng, cao tầng, đều đặn và không đều đặn) trong các vùng địa chấn khác nhau tại Việt Nam là bao nhiêu?
• H1: Việc số hóa các đại lượng mờ thành hàm thuộc $\mu(x)$ kết hợp biến đổi định lượng tương đương cho phép loại bỏ sự thiên vị chủ quan của người thiết kế, tạo lập miền phản ứng phổ tiệm cận với bản chất ngẫu nhiên mờ của động đất.
• H2: Thuật toán Monte - Carlo cải tiến (MMC) xử lý mẫu toàn thể có gán trọng số khả dĩ sẽ thiết lập được bao nội lực và chuyển vị an toàn hơn, phản ánh chính xác các vùng ứng xử phi tuyến mà phương pháp tất định bỏ sót.

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 Lofti A. Zadeh, 1965), Lý thuyết độ tin cậy kết cấu (Structural Reliability Theory), Động lực học công trình phi tuyến (Nonlinear Structural Dynamics) và Phương pháp số Newmark tích phân từng bước. Phạm vi khảo nghiệm bao trùm cả nhà thấp tầng (công trình Nhà văn hóa, mô hình khung 4 tầng, 6 tầng, 8 tầng) và nhà cao tầng (tổ hợp chung cư CT1, CT2, CT3 quy mô 20 tầng, 30 tầng, 40 tầng) tại các địa bàn địa chấn trọng điểm: Mường Lay, Tuần Giáo, Ba Vì, Hà Nội (Gia Lâm, Ba Đình), Hải Phòng (Ngô Quyền), Đà Nẵng (Ngũ Hành Sơn) và TP. Hồ Chí Minh (Quận 10).

Literature Review và Positioning

Lịch sử phát triển của kỹ thuật kháng chấn toàn cầu ghi nhận các bước nhảy vọt về nhận thức luận. Khởi đầu từ khái niệm "độ chấn" đơn giản của Giáo sư Riki Sano (1915) sau động đất San Francisco (1906) và Kanto (1924) với việc lượng hóa lực quán tính $F = R \cdot W$, lý thuyết phổ phản ứng đàn hồi chính thức được thiết lập bởi Maurice Anthony Biot (1932) và hoàn thiện bởi George W. Housner (1959). Tiếp đó, phương pháp tích phân số từng bước của Nathan M. Newmark (1959) đã mở đường cho việc giải phương trình vi phân dao động của hệ nhiều bậc tự do (MDOF) chịu gia tốc nền theo lịch sử thời gian $u_g(t), \dot{u}_g(t), \ddot{u}_g(t)$.

Trong việc thiết lập tương quan giữa gia tốc cực đại trung bình $a_{max}$ với Magnitude ($M$) và cự ly chấn tiêu ($R$), mô hình kinh nghiệm của Donovan (1973) được ứng dụng rộng rãi: $$a_{max} = \frac{1320 e^{0.58M}}{(R+25)^{1.52}} \quad (\text{cm/s}^2)$$ Tuy nhiên, khi tiến gần tâm chấn ($R \le 15 - 20\text{ km}$), tương quan giải tích này phân tán dữ liệu rất lớn, đòi hỏi phải tiếp cận qua mô hình xác suất và lý thuyết mờ.

Trong trường phái mô phỏng độ tin cậy địa chấn quốc tế, cuộc tranh luận khoa học diễn ra gay gắt giữa hai hướng tiếp cận chính:

1. Trường phái xác suất ngẫu nhiên thuần túy: Đại diện bởi Siu-Kui Au và James L. Beck (2001, 2003) với kỹ thuật mô phỏng tập con (Subset Simulation) dựa trên thuật toán sửa đổi Metropolis-Hastings (Markov Chain Monte Carlo - MCMC) để ước tính xác suất phá hoại nhỏ trong không gian nhiều chiều; Carlo Del Gaudio và cộng sự (2017) với mô hình phân tích đường cong tổn thương địa chấn (fragility curves) khu vực đô thị; Gian Felice Giaccu và Luca Caracoglia (2018) với phương pháp xấp xỉ ngẫu nhiên phân lớp (Layered Stochastic Approximation Monte Carlo - LSAMC). Hạn chế của trường phái này là giả định mọi biến số bất định đều tuân theo các hàm mật độ xác suất đã biết ($PDF$), bỏ qua tính thiếu hụt thông tin và sự mơ hồ ngữ nghĩa trong các tiêu chuẩn kỹ thuật.
2. Trường phái số học mờ và ngẫu nhiên mờ (Fuzzy Randomness): Khởi xướng mạnh mẽ bởi Bernd Möller và Michael Beer (2004) qua công trình kinh điển "Fuzzy Randomness: Uncertainty in Civil Engineering and Computational Mechanics", cùng với Schuëller (2004), Xinggang Bai và Shorab Asgarpoor (2004), David Moens và Dirk Vandepitte (2005). Möller và Beer đề xuất phương pháp tối ưu hóa mức $\alpha$ ($\alpha\text{-level optimization}$) để tìm kiếm hàm thuộc kết quả bằng cách giải hàng loạt bài toán biên tất định. Tuy nhiên, rào cản lớn nhất của phương pháp này là khối lượng tính toán bùng nổ cấp số nhân khi số bậc tự do lớn, đồng thời ép buộc các tham số tất định và ngẫu nhiên thực nghiệm phải chuyển đổi ngược về dạng mờ, gây khó khăn cho kỹ thuật công trình.

Tại Việt Nam, các nghiên cứu của Lê Xuân Huỳnh và Nguyễn Hùng Tuấn (2016), Nguyễn Đình Xân (2005), Lê Công Duy (2014), Nguyễn Thanh Hưng (2017) đã tiếp cận phương pháp phần tử hữu hạn mờ (Fuzzy FEM) và độ tin cậy mờ cho các cấu kiện đơn lẻ (dầm bê tông cốt thép ăn mòn, khung phẳng 2D chịu tải tĩnh hoặc động đơn giản) theo mức tối ưu $\alpha$. Nhóm nghiên cứu của Nguyễn Văn Phó, Chu Thanh Bình, Vũ Trọng Huy (2006, 2013, 2017, 2019, 2021) đã phát triển mô hình giao thoa và đặt nền móng cho thuật toán Monte - Carlo cải tiến.

Luận án của Vũ Trọng Huy đã định vị chính xác khoảng trống chưa ai giải quyết: Ứng dụng phương pháp Monte - Carlo cải tiến vào bài toán kết cấu không gian 3D tổng thể của công trình nhà chịu động đất theo tiêu chuẩn TCVN 9386:2012, trực tiếp số hóa 6 nhóm tham số mờ quy định trong tiêu chuẩn và liên kết tính toán với phần mềm kỹ thuật ETABS v18.

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

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

Luận án đã mở rộng Lý thuyết tập mờ của Zadeh và Lý thuyết ngẫu nhiên mờ của Möller & Beer vào việc diễn giải các điều khoản kỹ thuật của tiêu chuẩn xây dựng quốc gia. Luận án chứng minh một cách tường minh rằng:

"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."

Bản chất của sự chuyển dịch nhận thức luận (epistemological paradigm shift) nằm ở việc thay thế mô hình đánh giá "đúng/sai" tất định bằng mô hình "độ thuộc/khả dĩ". Mô hình lý thuyết đề xuất bao gồm hệ thống các mệnh đề khoa học:

• Mệnh đề 1 (Đặc trưng hóa 3 nhóm biến): Không gian đầu vào của bài toán kháng chấn kết cấu $\Omega$ được cấu thành từ tập hợp tích hợp $\Omega = {X_{det}} \cup {X_{ran}} \cup {X_{fuz}}$, trong đó $X_{det}$ là véc-tơ tham số tất định (kích thước tiết diện hình học), $X_{ran}$ là véc-tơ biến ngẫu nhiên (mô đun đàn hồi $E$, tĩnh tải $G$, hoạt tải $Q$ với hàm phân bố xác suất $f(x)$), và $X_{fuz}$ là véc-tơ các biến mờ trong tiêu chuẩn kháng chấn ($a_g, q, S, T_B, T_C, T_D, \eta$).
• Mệnh đề 2 (Tính tương đương định lượng): Một tham số mờ cấp một $P$ với hàm thuộc $\mu_P(x)$ có thể được rời rạc hóa thành các biến mờ thứ cấp $P_i$ có trọng số $w_i$, cho phép thiết lập không gian mẫu thử toàn thể kết hợp với các biến ngẫu nhiên mà không làm biến dạng bản chất vật lý của hệ kết cấu.
• Mệnh đề 3 (Giải mờ bao đáp ứng): Nội lực ($M, N, V$) và chuyển vị đỉnh/chuyển vị tầng ($d_x, d_y, \Delta/h$) của hệ kết cấu không gian là tập hợp mờ của các phản ứng thành phần, được giải mờ qua phương pháp trọng tâm có trọng số, tạo nên đường bao thiết kế an toàn tối ưu.

+-----------------------------------------------------------------------------------+
|               KHUNG PHÂN TÍCH TỔNG HỢP MONTE - CARLO CẢI TIẾN (MMC)               |
+-----------------------------------------------------------------------------------+
|  [Nhóm 1: Tất định (Det)]   |  [Nhóm 2: Ngẫu nhiên (Ran)] | [Nhóm 3: Biến Mờ (Fuz)]|
|  - Tiết diện hình học       |  - Mô đun đàn hồi E         | - Đỉnh gia tốc nền agR|
|  - Cường độ danh định fck   |  - Tĩnh tải (G), Hoạt tải(Q)| - Hệ số ứng xử q      |
|                             |    (Hàm phân bố Gauss, Weib)| - Tham số nền: S,TB,TC|
+-----------------------------+-----------------------------+-----------------------+
                                             |
                                             v
                      +---------------------------------------------+
                      |   SỐ HÓA & BIẾN ĐỔI ĐỊNH LƯỢNG TƯƠNG ĐƯƠNG  |
                      |   - Xây dựng hàm thuộc mu(x) (tam giác/thang)|
                      |   - Phân rã biến mờ thứ cấp & gán trọng số w|
                      +---------------------------------------------+
                                             |
                                             v
                      +---------------------------------------------+
                      |   THUẬT TOÁN MONTE - CARLO CẢI TIẾN (MMC)   |
                      |   - Thiết lập tổ hợp phổ khả dĩ (Sd(T)_k)   |
                      |   - Tích hợp mô phỏng số ngẫu nhiên         |
                      +---------------------------------------------+
                                             |
                                             v
                      +---------------------------------------------+
                      |   PHÂN TÍCH KẾT CẤU TRÊN ETABS v18          |
                      |   - Phổ phản ứng đa trường hợp              |
                      |   - Phân tích phi tuyến lịch sử thời gian   |
                      |   - Mô hình hóa tường chèn (Gap Element)    |
                      +---------------------------------------------+
                                             |
                                             v
                      +---------------------------------------------+
                      |   GIẢI MỜ ĐÁP ỨNG ĐẦU RA (NỘI LỰC, CHUYỂN VỊ|
                      |   - Bao nội lực phản lực chân cột (N, V)    |
                      |   - Chuyển vị ngang và độ lệch tầng (dx, dy)|
                      |   - Đánh giá chênh lệch so với TCVN 9386    |
                      +---------------------------------------------+

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

Luận án đã phân loại và số hóa thành công 6 nhóm tham số mờ cốt lõi trong tiêu chuẩn TCVN 9386:2012:

1. Tham số đỉnh gia tốc nền tham chiếu ($a_{gR}$): Tiêu chuẩn quy định $a_{gR}$ theo các đường đẳng trị và cho phép nội suy tuyến tính. Giá trị thực tế nằm trong khoảng $[a_{gR1}, a_{gR2}]$. Luận án số hóa bằng hàm thuộc tam giác hoặc hình thang, xác định mức độ tin tưởng tại tâm vùng và độ dốc suy giảm ra biên.
2. Hệ số ứng xử kết cấu ($q$): Phản ánh khả năng tiêu tán năng lượng dẻo ngoài miền đàn hồi. Thay vì chọn giá trị cố định $q = q_0 \cdot k_w$, luận án mô hình hóa $q$ dưới dạng số mờ biến thiên trong khoảng $[a, b]$ tùy thuộc vào độ dẻo kết cấu (DCM, DCH), tính đều đặn mặt bằng/mặt đứng và sự tham gia của vách cứng.
3. Các tham số phổ đất nền ($S, T_B, T_C, T_D$): Nhận dạng địa tầng qua vận tốc sóng cắt $v_{s,30}$ hoặc số nhát đập xuyên tiêu chuẩn $N_{SPT}$ có ranh giới mờ giữa các loại đất (A, B, C, D, E). Hàm thuộc của $S, T_B, T_C$ được thiết lập dưới dạng hình thang mờ đa chiều.
4. Hệ số tầm quan trọng công trình ($\gamma_I$): Phân cấp công trình từ cấp I ($\gamma_I = 1.25$), cấp II ($\gamma_I = 1.0$) đến cấp III ($\gamma_I = 0.75$). Luận án xây dựng hàm thuộc hình chuông và hình chữ nhật cho các công trình đa công năng.
5. Hệ số điều chỉnh độ cản ($\eta$): Biểu diễn sự biến thiên của tỷ số cản nhớt $\xi$ khác $5%$ thông qua công thức $\eta = \sqrt{10 / (5 + \xi)} \ge 0.55$, với $\xi$ là biến mờ phụ thuộc vật liệu hỗn hợp bê tông - cốt thép - tường xây.
6. Mô hình tường chèn tương đương: Mô hình hóa thanh chống chịu nén tương đương kết hợp phần tử phi tuyến tiếp xúc (Gap element) nhằm giải quyết sự không chắc chắn về độ cứng đóng góp của khối xây.

Điều kiện biên lý thuyết (boundary conditions) được xác lập: các biến mờ được giả thiết độc lập trong không gian xác suất, phổ thiết kế $S_d(T)$ tuân thủ giới hạn cận dưới $S_d(T) \ge 0.2 a_g$, và chu kỳ dao động công trình $T \le 4.0\text{ s}$.

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

Thiết kế nghiên cứu

Nghiên cứu được kiến tạo trên lập trường thực chứng thực tế (Pragmatic Critical Realism), dung hòa giữa các quy luật cơ học kết cấu giải tích chính xác và sự bất định cố hữu của hiện tượng tự nhiên. Thiết kế nghiên cứu là sự kết hợp đa tầng (multi-level simulation design):

• Cấp độ vật liệu/cấu kiện: Khảo sát phi tuyến vật liệu bê tông cốt thép, biến dạng dẻo cục bộ tại các vị trí khớp dẻo (plastic hinges).
• Cấp độ hệ kết cấu tổng thể: Mô hình hóa không gian 3D đầy đủ 6 bậc tự do tại mỗi nút, khảo sát tương tác giữa khung, vách chịu cắt và khối xây tường chèn.
• Cấp độ ngoại lực địa chấn: Chuyển hóa phổ phản ứng chuẩn tắc TCVN 9386 thành một họ các hàm phổ khả dĩ $S_d(T)_k$ kèm theo trọng số xác suất mờ.

Quy trình nghiên cứu rigorous

Quy trình tính toán theo phương pháp Monte - Carlo cải tiến (MMC) được triển khai qua các bước chuẩn hóa nghiêm ngặt:

[BƯỚC 1: XÁC ĐỊNH KHÔNG GIAN ĐẦU VÀO]
- Tất định: Kích thước hình học, tải trọng bản thân kết cấu.
- Ngẫu nhiên: E, tải trọng sử dụng (hàm Gauss, N(mu, sigma^2)).
- Mờ: agR, q, S, TB, TC, TD (thiết lập hàm thuộc mu(x)).
          |
          v
[BƯỚC 2: MỜ HÓA VÀ RỜI RẠC HÓA THAM SỐ]
- Chia khoảng biến thiên của các biến mờ thành m mức rời rạc x_i.
- Tính mức độ thuộc mu(x_i) và chuẩn hóa thành trọng số w_i = mu(x_i) / SUM(mu(x_k)).
- Tổ hợp các trường hợp phổ thiết kế khả dĩ k = 1..K với trọng số W_k = PROD(w_i).
          |
          v
[BƯỚC 3: MÔ PHỎNG SỐ VÀ PHÂN TÍCH KẾT CẤU]
- Khởi tạo số giả ngẫu nhiên cho các biến ngẫu nhiên theo phân bố xác suất.
- Nhập ma trận phổ phản ứng khả dĩ Sd(T)_k vào phần mềm ETABS v18.
- Chạy phân tích động đàn hồi (Response Spectrum Analysis) và phi tuyến (Time-History Analysis).
          |
          v
[BƯỚC 4: GIẢI MỜ VÀ TỔNG HỢP ĐÁP ỨNG]
- Thu nhận tập kết quả phản lực đáy Fb_k, chuyển vị d_k, nội lực M_k, V_k.
- Giải mờ theo phương pháp trung bình trọng số: R_final = SUM(W_k * R_k).
- Thiết lập đường bao ứng xử so sánh trực tiếp với phương pháp tất định chuẩn.

Độ tin cậy và tính hợp thức (validity & reliability) của thuật toán được bảo đảm thông qua:

• Triangulation phương pháp: Kiểm chứng chéo giữa phân tích phổ phản ứng đàn hồi (RSA), phân tích tĩnh phi tuyến (Pushover) và phân tích động phi tuyến theo lịch sử thời gian (Nonlinear Time History Analysis - NTHA) sử dụng giản đồ gia tốc El Centro chuẩn hóa.
• Kiểm soát hội tụ: Số lượng mẫu mô phỏng ngẫu nhiên mờ được kiểm định với dung sai sai số tương đối $\varepsilon \le 1.5%$ ở độ tin cậy $95%$.

Data và phân tích

Ma trận dữ liệu thực nghiệm mô phỏng bao gồm 2 nhóm công trình lớn:

Mô hình tường chèn được khai báo bằng thanh chống xiên chịu nén tương đương với bề rộng tính toán $b_w$ và phần tử Gap chỉ chịu nén, không chịu kéo, triệt tiêu liên kết khi phát sinh chuyển vị tách rời khung.

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

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

1. Sự đánh giá thấp nguy hiểm của phương pháp tất định trong vùng động đất mạnh: Tại các khu vực có gia tốc nền cao như Mường Lay ($a_g \approx 0.218g$) và Tuần Giáo ($a_g \approx 0.187g$), phương pháp tất định theo TCVN 9386:2012 cho kết quả phản lực chân công trình và mô men uốn đáy nhỏ hơn đáng kể so với phương pháp Monte - Carlo cải tiến (MMC). Tỷ lệ sai lệch phản lực chân cột đạt từ $14.2%$ đến $28.6%$. Điều này chứng minh việc chọn một giá trị tất định $q$ cố định (ví dụ $q = 3.9$) đã vô tình triệt tiêu các kịch bản bất lợi khi cấu kiện không đạt độ dẻo lý tưởng.

2. Khuếch đại độ lệch chuyển vị đối với kết cấu không đều đặn mặt bằng: Đối với công trình 8 tầng không đều đặn, sự biến thiên mờ của tâm cứng và độ lệch tâm ngẫu nhiên $e_{ai} = \pm 0.05 L_i$ kết hợp với phổ mờ MMC dẫn đến chuyển vị ngang đỉnh tầng ($d_x, d_y$) tăng vọt $18.7% - 31.4%$ so với tính toán tất định. Hiện tượng xoắn tổng thể công trình bị trầm trọng hóa do tương tác giữa mode dao động xoắn và tịnh tiến.

3. Tác động hai mặt của khối xây tường chèn trong phân tích phi tuyến: Khi phân tích lịch sử thời gian có xét phần tử Gap cho tường chèn, độ cứng tổng thể ban đầu của công trình tăng lên làm giảm chu kỳ dao động riêng cơ bản ($T_1$ giảm từ $0.82\text{ s}$ xuống $0.54\text{ s}$ ở nhà văn hóa thấp tầng), dẫn đến lực cắt đáy tăng vọt $35.2%$. Tuy nhiên, khi xảy ra phá hoại giòn ở khối xây, sự thay đổi độ cứng đột ngột gây ra hiện tượng tập trung ứng suất cắt cực lớn tại các đầu cột tầng trệt (hiệu ứng tầng mềm).

4. Độ nhạy của nhà cao tầng (20T, 30T, 40T) đối với tham số đất nền mờ: Tại khu vực nền đất yếu Gia Lâm (Hà Nội) và Ngô Quyền (Hải Phòng) (đất loại D, E với hệ số $S = 1.4 - 1.6$), chuyển vị ngang tương đối giữa các tầng ($d_r$) theo phương pháp MMC vượt giới hạn cho phép của TCVN 9386:2012 tới $16.8%$ ở các tầng từ 15 đến 28 của tòa nhà CT3 (40 tầng). Phương pháp tất định đã bỏ qua dải bất định của chu kỳ $T_C, T_D$, dẫn tới việc tính thiếu cốt thép vách kháng chấn.

+-----------------------------------------------------------------------------------+
|     BIỂU ĐỒ SO SÁNH PHẢN LỰC ĐÁY VÀ CHUYỂN VỊ GIỮA TCVN 9386 VÀ MMC              |
+-----------------------------------------------------------------------------------+
|  Lực cắt đáy Fb (kN)                                                              |
|  4500 |                                                 [MMC: 4120 kN]            |
|  4000 |                                              +-----------------+          |
|  3500 |                      [TCVN: 3250 kN]         |                 |          |
|  3000 |                   +-----------------+        |  Chênh lệch     |          |
|  2500 |                   |                 |        |  +26.77%        |          |
|  2000 |                   |                 |        | (Vùng Mường Lay)|          |
|       +-------------------+-----------------+--------+-----------------+          |
|                            Phương pháp TCVN           Phương pháp MMC             |
+-----------------------------------------------------------------------------------+

Implications đa chiều

• Về mặt học thuật và lý thuyết: Xác lập khung phương pháp luận kết hợp ngẫu nhiên - mờ hoàn chỉnh cho cơ học công trình tại Việt Nam, giải quyết triệt để bài toán thiếu dữ liệu quan trắc địa chấn lịch sử thông qua xử lý số hóa thông tin mờ.
• Về mặt thực tiễn thiết kế: Cung cấp quy trình tích hợp thuật toán MMC trực tiếp vào phần mềm ETABS v18 thông qua tổ hợp tải trọng phổ có trọng số, giúp các kỹ sư tư vấn thiết kế tạo lập được đường bao nội lực an toàn thực chất mà không cần viết lại toàn bộ phần mềm phần tử hữu hạn.
• Về mặt chính sách và quy chuẩn: Đóng vai trò cơ sở khoa học then chốt kiến nghị Ban soạn thảo tiêu chuẩn quốc gia sửa đổi TCVN 9386:2012 theo hướng: chuyển đổi việc tra cứu tham số cứng nhắc sang định dạng dải tin cậy hoặc bổ sung hệ số an toàn mờ cho các công trình có cấp tầm quan trọng cao ($\gamma_I > 1.0$).

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật:

• Giả thiết tính độc lập của các biến mờ: Do hạn chế về dữ liệu thực nghiệm địa chất sâu, các biến mờ ($a_g, q, S$) được giả thiết là độc lập thống kê. Trên thực tế, giữa loại đất nền ($S$) và chu kỳ góc ($T_C$) tồn tại mối tương quan phi tuyến chặt chẽ.
• Số lượng giản đồ sóng địa chấn lịch sử: Chưa có đủ cơ sở dữ liệu gia tốc đồ mạnh ghi nhận tại các đới đứt gãy nội địa Việt Nam, việc phân tích phi tuyến thời gian vẫn phải dựa vào giản đồ chuẩn hóa quốc tế (El Centro 1940).
• Mô hình tương tác đất - kết cấu (SSI): Công trình ngàm cứng tại mặt móng, chưa xét đến tính mềm của cọc và sự tiêu tán năng lượng qua nền đất (Soil-Structure Interaction).

Chương trình nghiên cứu 5 - 10 năm tới bao gồm:

1. Phát triển mô hình hàm liên hợp mờ (Fuzzy Copula Theory) để mô tả cấu trúc tương quan phụ thuộc giữa các tham số kháng chấn.
2. Mở rộng thuật toán MMC cho bài toán phân tích sụp đổ động lực học phi tuyến 3D (Incremental Dynamic Collapse Analysis) dưới tác động của sóng địa chấn trường gần có hiệu ứng xung vận tốc (near-fault pulse-like ground motions).
3. Ứng dụng trí tuệ nhân tạo và học máy (Machine Learning Surrogate Models) để tăng tốc độ mô phỏng số MMC cho các siêu đô thị phức hợp.

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

• Tác động học thuật: Luận án mở ra hướng nghiên cứu đột phá về "Độ tin cậy mờ công trình" tại Việt Nam, dự kiến thu hút hàng chục trích dẫn khoa học trên các tạp chí uy tín thuộc danh mục ISI/Scopus (như Engineering Structures, Soil Dynamics and Earthquake Engineering, Tạp chí Khoa học Công nghệ Xây dựng - IBST).
• Chuyển đổi ngành công nghiệp xây dựng: Giúp các tổng công ty tư vấn thiết kế lớn (VNCC, CDC, Texo, Coninco) tối ưu hóa việc bố trí cốt thép và hình học kết cấu vách, giảm thiểu chi phí khắc phục hư hỏng sau động đất lên tới hàng tỷ đồng cho mỗi dự án chung cư cao tầng.
• An sinh xã hội và giảm nhẹ thiên tai: Nâng cao ngưỡng an toàn sinh mạng cho hàng triệu cư dân sinh sống tại các đô thị nhạy cảm địa chấn (Hà Nội, Điện Biên, Lai Châu), bảo vệ tính toàn vẹn của các công trình công cộng trọng điểm (bệnh viện, trường học, nhà văn hóa).

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

• Nghiên cứu sinh và học viên cao học: Tiếp cận phương pháp luận mô phỏng số MMC hoàn chỉnh, mã nguồn thuật toán xử lý biến mờ và quy trình tích hợp phần mềm thương mại.
• Giảng viên và nhà khoa học chuyên ngành Kết cấu - Công trình: Tài liệu tham khảo chuẩn mực phục vụ giảng dạy môn Động lực học công trình nâng cao, Độ tin cậy kết cấu và Thiết kế công trình chịu tải trọng đặc biệt.
• Kỹ sư kết cấu thực hành (R&D Engineering): Nắm vững kỹ thuật lập phổ thiết kế bao an toàn và phương pháp gán phần tử Gap mô phỏng tường chèn phi tuyến.
• Cơ quan quản lý nhà nước (Bộ Xây dựng, Viện KHCN Xây dựng): Dữ liệu thực chứng đắt giá phục vụ lộ trình cập nhật, biên soạn Quy chuẩn Kỹ thuật Quốc gia và Tiêu chuẩn Thiết kế Kháng chấn thế hệ mới.

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 định lượng tương đương giữa các đại lượng mờ ngôn ngữ trong tiêu chuẩn xây dựng với không gian biến số ngẫu nhiên. Luận án đã mở rộng Lý thuyết tập mờ của Lofti Zadeh (1965) và Lý thuyết ngẫu nhiên mờ của Möller & Beer (2004) từ phạm vi cơ học lý thuyết thuần túy sang ứng dụng trực tiếp vào tiêu chuẩn kỹ thuật TCVN 9386:2012, chứng minh rằng phổ thiết kế kháng chấn thực chất là một hàm mờ đa chiều có thể giải mờ qua trọng số khả dĩ.

2. Sự đổi mới về phương pháp luận của luận án khi so sánh với các nghiên cứu quốc tế trước đây?

So với phương pháp mô phỏng tập con (Subset Simulation) của Siu-Kui Au và James L. Beck (2001, 2003) vốn chỉ xử lý biến ngẫu nhiên, và phương pháp tối ưu hóa mức $\alpha$ của Bernd Möller và Michael Beer (2004) vốn làm bùng nổ khối lượng tính toán và biến đổi gượng ép mọi biến thành mờ, thuật toán MMC của luận án duy trì nguyên vẹn bản chất 3 nhóm biến (Tất định - Ngẫu nhiên - Mờ). Thuật toán xử lý trên mẫu toàn thể có gán trọng số và liên kết trực tiếp với phần mềm thương mại ETABS v18, giảm $80%$ thời gian tính toán so với phương pháp tối ưu $\alpha$ mà vẫn đảm bảo độ chính xác kỹ thuật.

3. Phát hiện bất ngờ nhất (counter-intuitive finding) có số liệu chứng minh là gì?

Phát hiện bất ngờ nhất là: Việc xây chèn tường gạch vào khung BTCT nhà thấp tầng tưởng chừng làm tăng độ an toàn nhưng thực tế lại làm tăng nguy cơ phá hoại giòn tầng trệt khi chịu động đất. Số liệu chứng minh: chu kỳ $T_1$ giảm từ $0.82\text{ s}$ xuống $0.54\text{ s}$ khiến gia tốc phổ $S_d(T)$ rơi vào đỉnh cộng hưởng, lực cắt đáy $F_b$ tăng vọt $35.2%$, gây biến dạng trượt cắt phá hủy đầu cột trước khi dầm kịp phát triển khớp dẻo tiêu tán năng lượng.

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

Luận án cung cấp quy trình nhân bản 4 bước chi tiết:

1. Bảng số hóa hàm thuộc 6 tham số tiêu chuẩn ($\mu(a_g), \mu(q), \mu(S), \dots$).
2. Thuật toán tạo số giả ngẫu nhiên và phân rã biến mờ thứ cấp trên mẫu toàn thể.
3. Bảng tra trọng số $W_k$ nhập liệu trực tiếp các trường hợp phổ vào Load Cases của ETABS v18.
4. Công thức giải mờ nội lực đầu ra bằng giải thuật tích phân trọng tâm.

5. Kế hoạch nghiên cứu 10 năm tiếp theo được định hình như thế nào?

Lộ trình 10 năm định hình 3 mục tiêu chiến lược:

• Xây dựng bản đồ nguy cơ địa chấn mờ động (Dynamic Fuzzy Seismic Hazard Map) cho toàn lãnh thổ Việt Nam tích hợp dữ liệu GIS.
• Hoàn thiện gói phần mềm chuyên dụng tích hợp thuật toán MMC tự động hóa hoàn toàn với các phần mềm kết cấu BIM (Revit, ETABS, SAP2000).
• Mở rộng tiêu chuẩn hóa phương pháp thiết kế dựa trên hiệu năng mờ (Fuzzy Performance-Based Seismic Design - FPBSD) cho các công trình siêu cao tầng và công trình ngầm đô thị.

Kết luận

1. Số hóa toàn diện tiêu chuẩn TCVN 9386:2012: Lần đầu tiên tại Việt Nam, 6 nhóm tham số mờ cốt lõi trong tiêu chuẩn thiết kế kháng chấn đã được nhận diện, phân loại và toán học hóa thành các hàm thuộc chính xác.
2. Hoàn thiện phương pháp Monte - Carlo cải tiến (MMC): Xây dựng thành công thuật toán mô phỏng kết cấu không gian 3D xử lý đồng thời 3 nhóm biến: tất định, ngẫu nhiên và mờ.
3. Định lượng mức độ sai lệch của phương pháp tất định: Chỉ ra sai số thiếu an toàn từ $14.2%$ đến $28.6%$ về phản lực chân công trình và hơn $30%$ về chuyển vị xoắn ở các vùng địa chấn mạnh khi áp dụng phương pháp tất định truyền thống.
4. Làm rõ ứng xử phi tuyến của kết cấu có tường chèn: Xác lập mô hình phần tử Gap mô phỏng chính xác hiện tượng suy giảm độ cứng và tập trung ứng suất cắt phá hoại.
5. Mở ra 3 nhánh nghiên cứu học thuật mới: Độ tin cậy ngẫu nhiên mờ trong kết cấu công trình; Tối ưu hóa kháng chấn dựa trên lý thuyết khả dĩ; Tích hợp thuật toán số hóa mờ vào nền tảng thiết kế tự động hóa 4.0.
6. Giá trị thực tiễn và di sản khoa học: Đóng góp luận cứ khoa học vững chắc phục vụ công tác sửa đổi tiêu chuẩn xây dựng quốc gia, bảo vệ an toàn sinh mạng cộng đồng và nâng cao khả năng chống chịu động đất của cơ sở hạ tầng Việt Nam.