Tổng quan về luận án

Nghiên cứu tác động của tải trọng gió lên kết cấu công trình cao tầng là một trong những bài toán then chốt của cơ học công trình và kỹ thuật gió hiện đại. Tại Việt Nam, đặc biệt là dải duyên hải Miền Trung, các công trình cao tầng thường xuyên phải đối mặt với các cơn bão nhiệt đới có cường độ mạnh kèm theo các xung vận tốc biến đổi phức tạp theo không gian và thời gian. Luận án tiến sĩ kỹ thuật chuyên ngành Cơ kỹ thuật (Mã số: 62 52 01 01) của nghiên cứu sinh Bùi Thiên Lam, dưới sự hướng dẫn khoa học của PGS.TS. Phan Quang Minh và PGS.TS. Lê Cung tại Đại học Đà Nẵng (2018), mang tên "Nghiên cứu xác định tải trọng gió lên nhà cao tầng có kết cấu khung giằng", đại diện cho một bước tiến tiên phong trong việc chuẩn hóa và tối ưu hóa quy trình tính toán động lực học công trình chịu gió bão tại Việt Nam.

Luận án khẳng định tính cấp thiết kinh tế – kỹ thuật cốt lõi: "Chi phí cho hệ kết cấu để chịu tải trọng thẳng đứng (trọng lượng bản thân và hoạt tải sử dụng) tỉ lệ bậc nhất với chiều cao nhà. Trong khi đó chi phí cho hệ kết cấu để chịu tải trọng ngang tăng nhanh hơn nhiều khi chiều cao nhà tăng lên" (Trang 1). Trong khi các tiêu chuẩn quốc tế tiên tiến như ASCE/SEI 7-10 (Hoa Kỳ), Eurocode EN 1991-1-4 (Châu Âu) hay AS/NZS 1170.2 (Úc/New Zealand) đã ứng dụng rộng rãi phương pháp Hệ số gió giật (Gust Loading Factor – GLF) do Alan Garnett Davenport phát triển từ năm 1967, tiêu chuẩn Việt Nam hiện hành TCVN 2737:1995 (biên soạn dựa trên tiêu chuẩn Liên Xô SNiP 2.01.07-85*) vẫn áp dụng phương pháp tách biệt thành phần gió tĩnh và thành phần gió động thông qua phân tích dạng dao động.

Khoảng trống nghiên cứu (Research Gap) then chốt được luận án xác định cụ thể: Công thức gần đúng xác định thành phần gió động trong Điều 5.5 của TCVN 2737:1995 mặc định dạng dao động riêng thứ nhất của công trình biến thiên tuyến tính theo chiều cao $\phi_1(z) \approx z/H$. Giả thiết này chỉ phù hợp với kết cấu biến dạng cắt thuần túy (khung thuần). Đối với hệ kết cấu khung giằng (Dual Frame-Shear Wall System) – giải pháp kết cấu phổ biến nhất cho nhà cao tầng hiện đại – sự tương tác uốn của vách cứng và biến dạng cắt của khung tạo nên đường cong chuyển vị uốn-cắt phi tuyến dạng chữ S. Việc áp đặt giả thiết bậc nhất dẫn đến sai số rất lớn về phân bố tải trọng gió theo chiều cao và lực cắt đáy công trình.

Nghiên cứu đặt ra 3 câu hỏi nghiên cứu và hệ thống giả thuyết khoa học tương ứng:

  • Câu hỏi nghiên cứu 1 (RQ1): Mức độ sai số định lượng của công thức gần đúng trong TCVN 2737:1995 khi áp dụng cho hệ khung giằng có độ cứng tương đối $\alpha H$ khác nhau là bao nhiêu?
    Giả thuyết 1 (H1): Công thức gần đúng TCVN 2737:1995 gây ra sai số vượt quá $30%$ tại các tầng trên cao và đánh giá sai lệch lực cắt đáy đáng kể khi thông số $\alpha H > 2.0$.
  • Câu hỏi nghiên cứu 2 (RQ2): Làm thế nào để thiết lập một hàm phân bố giải tích mới thay thế hệ số tuyến tính $z/H$ mà vẫn giữ được cấu trúc tính toán đơn giản, phản ánh chính xác tương tác khung – vách?
    Giả thuyết 2 (H2): Tích hợp nghiệm phương trình vi phân tương tác môi trường liên tục (Continuum Method) vào hàm tải trọng động $K_1(\alpha H, z/H)$ sẽ kiểm soát sai số lực cắt đáy dưới $5%$.
  • Câu hỏi nghiên cứu 3 (RQ3): Có thể xây dựng một công thức tính hệ số gió giật $G$ (GLF) tường minh, tương thích với cơ sở dữ liệu khí hậu bão và thông số khí động của TCVN 2737:1995 cho nhà cao tầng tại vùng duyên hải như Thành phố Đà Nẵng hay không?
    Giả thuyết 3 (H3): Hệ số gió giật $G$ quy đổi cho vùng gió II-B (áp lực $W_0 = 95\text{ daN/m}^2$, chu kỳ lặp 20 năm, địa hình B và C) có thể biểu diễn qua dạng hàm đại số phụ thuộc trực tiếp vào số tầng $n$ hoặc chiều cao $H$, đơn giản hóa triệt để quy trình thiết kế thực hành.

Phạm vi nghiên cứu bao quát các công trình bê tông cốt thép cao từ 15 đến 35 tầng ($H = 60\text{m} \div 120\text{m}$), mặt bằng đối xứng, bề rộng đón gió $B \le 50\text{m}$, tần số dao động cơ bản $f_1$ lớn hơn tần số giới hạn $f_L$ (công trình thuộc diện chỉ cần xét dạng dao động thứ nhất), khảo sát chuyên sâu trên điều kiện khí quyển và địa hình thực tế tại Thành phố Đà Nẵng.

Literature Review và Positioning

Lý thuyết kỹ thuật gió và động lực học công trình chịu tải trọng gió đã trải qua hơn nửa thế kỷ phát triển với hai trường phái chính:

Trường phái phân tích phổ dao động ngẫu nhiên và Hệ số gió giật (GLF) bắt đầu từ công trình kinh điển của Davenport (1961, 1967), sau đó được chuẩn hóa và mở rộng bởi Simiu & Scanlan (1996), Solari (1993), và Zhou, Kijewski & Kareem (2002). Davenport đã phân rã phản ứng động của công trình thành thành phần phản ứng nền (Background response $B$) do các xoáy khí quyển quy mô lớn tần số thấp gây ra và thành phần phản ứng cộng hưởng (Resonant response $R$) do các xoáy gần tần số dao động riêng của công trình gây ra. Hướng tiếp cận này trở thành nền tảng của ASCE/SEI 7-10 và Eurocode EN 1991-1-4.

Trường phái phương pháp lực quán tính tương đương dựa trên phân tích dạng dao động riêng biệt của Liên Xô (SNiP 2.01.07-85*), được chuyển giao trực tiếp vào TCVN 2737:1995 tại Việt Nam. Theo phương pháp này, thành phần gió động tại tầng $k$ được tính bằng tích của thành phần gió tĩnh, hệ số xung áp lực gió $\zeta$, hệ số tương quan không gian $\nu$, và thừa số dạng dao động $\xi_k = \frac{y_k \sum M_j y_j}{\sum M_j y_j^2}$. Để đơn giản hóa, TCVN 2737:1995 đưa ra công thức gần đúng (1.50): $$W_{dk} = \frac{k}{n} \cdot m_k \cdot \nu \cdot \xi \cdot W_{tk}$$ với giả thiết chuyển vị ngang $y_k \approx \frac{z_k}{H}$.

                      ┌────────────────────────────────────────┐
                      │    HỆ THỐNG PHÂN TÍCH GIÓ CÔNG TRÌNH   │
                      └──────────────────┬─────────────────────┘
                                         │
                 ┌───────────────────────┴───────────────────────┐
                 ▼                                               ▼
   ┌───────────────────────────┐                   ┌───────────────────────────┐
   │ TRƯỜNG PHÁI QUỐC TẾ (GLF) │                   │  TRƯỜNG PHÁI TCVN / SNiP  │
   │  Davenport (1967), ASCE   │                   │   TCVN 2737, SNiP 2.01    │
   │   EN 1991, Solari (1993)  │                   │  Tách biệt Gió Tĩnh/Động  │
   └─────────────┬─────────────┘                   └─────────────┬─────────────┘
                 │                                               │
                 │    ┌─────────────────────────────────────┐    │
                 └───►│      KHOẢNG TRỐNG NGHIÊN CỨU        │◄───┘
                      │  Giả thiết tuyến tính y(z) = z/H    │
                      │  gây sai số lớn trong Khung-Giằng   │
                      │  Chưa có GLF tương thích TCVN 2737  │
                      └──────────────────┬──────────────────┘
                                         │
                                         ▼
                      ┌─────────────────────────────────────┐
                      │   ĐÓNG GÓP ĐỘT PHÁ CỦA LUẬN ÁN      │
                      │ 1. Hàm phân bố phi tuyến K1(αH, z/H)│
                      │ 2. Công thức GLF đóng gói cho VN-IIB│
                      └─────────────────────────────────────┘

Tranh luận học thuật cốt lõi tồn tại giữa giả thiết đơn giản hóa dạng dao động bậc nhất và ứng xử thực tế của hệ kết cấu khung giằng. Các nghiên cứu cơ học kết cấu của Chitty (1947), Rosman (1964), và Stafford Smith & Coull (1991) đã chứng minh rằng kết cấu khung giằng chịu tải trọng ngang bị chi phối bởi phương trình vi phân bậc bốn: $$\frac{d^4 y}{dz^4} - \alpha^2 \frac{d^2 y}{dz^2} = \frac{q(z)}{EI_w}$$ trong đó $\alpha = \sqrt{\frac{GA_f}{EI_w}}$ biểu thị tỷ số giữa độ cứng chống trượt của khung ($GA_f$) và độ cứng chống uốn của vách ($EI_w$). Chuyển vị của hệ khung giằng kết hợp là sự chồng hình của đường cong uốn (vách) và đường cong cắt (khung), khiến cho dạng dao động thực tế $\phi_1(z)$ cong gập mạnh ở phần dưới và duỗi thẳng ở phần trên, khác biệt hoàn toàn với đường thẳng $y = z/H$.

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  1. Nghiên cứu của Boggs & Peterka (1989) và Holmes (2001) trong tiêu chuẩn ASCE 7 chỉ ra rằng việc giả định sai dạng dao động có thể làm sai lệch lực quán tính tương đương từ $15%$ đến $35%$, nhưng ASCE 7 giải quyết thông qua hệ số dạng dao động tích hợp trong phổ năng lượng gió.
  2. Tiêu chuẩn Eurocode EN 1991-1-4:2005 sử dụng hệ số số mũ dạng dao động $\zeta_1$ phụ thuộc vào sơ đồ kết cấu ($\zeta_1 = 1.0$ cho khung thuần, $\zeta_1 = 1.5$ cho khung giằng, $\zeta_1 = 2.0$ cho vách cứng thuần túy).

Tuy nhiên, các tiêu chuẩn này chưa được chuyển đổi tương thích với tập tham số khí tượng và phương pháp tách lực tĩnh/động của tiêu chuẩn xây dựng Việt Nam. Luận án của Bùi Thiên Lam định vị chính xác tại điểm giao thoa này: vừa hiệu chỉnh sai số dạng dao động cho phương pháp truyền thống TCVN 2737, vừa xây dựng cầu nối giải tích đưa phương pháp GLF Davenport vào ứng dụng thực hành tại Việt Nam.

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

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

Luận án đã mở rộng và tích hợp thành công hai trụ cột lý thuyết lớn của cơ học công trình: Lý thuyết môi trường liên tục phân tích kết cấu khung giằng (Continuum Theory for Coupled Frame-Wall Structures của Chitty, Rosman, Stafford Smith) và Lý thuyết dao động ngẫu nhiên kỹ thuật gió (Random Vibration Theory & Spectral Wind Engineering của Davenport).

                      ┌────────────────────────────────────────┐
                      │      KHUNG PHÂN TÍCH LÝ THUYẾT         │
                      └──────────────────┬─────────────────────┘
                                         │
        ┌────────────────────────────────┴────────────────────────────────┐
        ▼                                                                 ▼
┌───────────────────────────────┐                       ┌───────────────────────────────┐
│     LÝ THUYẾT MÔI TRƯỜNG      │                       │     LÝ THUYẾT DAO ĐỘNG        │
│          LIÊN TỤC             │                       │           NGẪU NHIÊN          │
│ (Chitty, Rosman, Smith-Coull) │                       │     (Davenport, Liepmann)     │
│   d⁴y/dz⁴ - α²d²y/dz² = q/EI  │                       │   S_w(n), GLF, G = 1 + g·I·R  │
└───────────────┬───────────────┘                       └───────────────┬───────────────┘
                │                                                       │
                └───────────────────────┬───────────────────────────────┘
                                        ▼
                      ┌───────────────────────────────────┐
                      │    HÀM PHÂN BỐ HIỆU CHỈNH K1      │
                      │  TẢI TRỌNG ĐỘNG DỌC CHIỀU CAO     │
                      └─────────────────┬─────────────────┘
                                        │
                                        ▼
                      ┌───────────────────────────────────┐
                      │    CÔNG THỨC HỆ SỐ GIÓ GIẬT G     │
                      │  CHUẨN HÓA VÙNG GIÓ II-B ĐÀ NẴNG  │
                      └───────────────────────────────────┘

Mô hình lý thuyết thiết lập mối quan hệ hàm số giữa phản ứng động lực học của công trình và thông số không thứ nguyên $\alpha H$:

  • Mệnh đề 1 (Proposition 1): Khi $\alpha H \le 2.0$, kết cấu làm việc chủ yếu theo sơ đồ vách cứng chịu uốn; dạng dao động bậc nhất có độ cong lớn tại gốc, hàm phân bố tải trọng động tập trung mạnh về phía đỉnh nhà.
  • Mệnh đề 2 (Proposition 2): Khi $\alpha H \in [2.0, 6.0]$, kết cấu thể hiện rõ nét tương tác hỗn hợp khung – vách; lực tương tác cắt giữa khung và vách đạt cực đại ở khoảng $(0.4 \div 0.6)H$, làm phẳng dạng dao động ở phần giữa thân nhà.
  • Mệnh đề 3 (Proposition 3): Thành phần gió động không biến thiên tuyến tính theo chiều cao mà biến thiên theo hàm số $K_1(\alpha H, \bar{z})$, trong đó $\bar{z} = z/H$.

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

Bản chất của đóng góp lý thuyết thể hiện qua việc thiết lập biểu thức tường minh cho hàm hiệu chỉnh $K_1$: $$K_1(\alpha H, \bar{z}) = \frac{y(\bar{z}) \int_0^1 \bar{m}(\bar{z}) y(\bar{z}) d\bar{z}}{\int_0^1 \bar{m}(\bar{z}) y^2(\bar{z}) d\bar{z}}$$

Bằng việc giải hệ phương trình vi phân biên cho dầm công-xôn tương đương chịu tải trọng ngang phân bố đều và tam giác, luận án đưa ra công thức giải tích thu gọn cho hệ số $K_1$: $$K_1(\alpha H, \bar{z}) = \left[ C_1 \cdot \bar{z} + C_2 \cdot \bar{z}^2 + C_3 \cdot (1 - \cosh(\alpha H \bar{z})) \right] \cdot \Phi(\alpha H)$$ trong đó $C_1, C_2, C_3$ là các hằng số tích phân phụ thuộc vào điều kiện biên tại ngàm chân cột và đỉnh mái tự do.

Điều kiện biên (Boundary conditions) của khung phân tích được xác định rõ ràng:

  • Công trình nhà cao tầng bằng bê tông cốt thép toàn khối, mặt bằng đối xứng tuyệt đối hoặc gần đối xứng (độ lệch tâm $e/B < 0.05$ để loại trừ dao động xoắn không gian).
  • Chiều cao công trình từ 20 đến 35 tầng ($60\text{m} \le H \le 120\text{m}$).
  • Khối lượng và độ cứng phân bố đều hoặc biến thiên tuyến tính nhẹ theo chiều cao.
  • Tần số dao động tự nhiên thứ nhất thỏa mãn $f_1 > f_L$ theo quy định của TCVN 2737:1995.

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ủ thế giới quan Thực chứng (Positivism) kết hợp Chủ nghĩa thực tế phản biện (Critical Realism) trong kỹ thuật công trình, đặt trọng tâm vào việc lượng hóa chính xác các hiện tượng cơ học thông qua giải tích toán học, mô phỏng số phần tử hữu hạn (FEM) và kiểm chứng thực nghiệm trên các công trình thực tế.

Thiết kế nghiên cứu đa cấp độ (Multi-level design):

  • Cấp độ 1 (Analytical Modeling): Thiết lập mô hình giải tích môi trường liên tục dầm tương đương cho hệ khung giằng và giải hệ phương trình vi phân chuyển vị bằng phương pháp giải tích chính xác.
  • Cấp độ 2 (Numerical Simulation): Xây dựng ma trận mô hình phần tử hữu hạn 3D không gian trên phần mềm SAP2000 và ETABS với 24 cấu hình kết cấu điển hình để trích xuất tần số riêng ($f_i$), chu kỳ dao động ($T_i$), và các vector dạng dao động (${\phi_i}$).
  • Cấp độ 3 (Empirical Case Verification): Kiểm chứng mô hình trên 2 công trình cao tầng biểu tượng thực tế tại Thành phố Đà Nẵng.
Cấp độ Phương pháp tiếp cận Công cụ / Phần mềm Đối tượng khảo sát Mục tiêu đầu ra
Cấp 1 Giải tích liên tục MATLAB, Maple, Excel Dầm công-xôn tương đương Hàm giải tích $K_1(\alpha H, \bar{z})$ và hệ số $G$
Cấp 2 Số học phần tử hữu hạn ETABS, SAP2000 v18 24 mô hình 20, 25, 30, 35 tầng Dạng dao động thực, chuyển vị đỉnh, lực cắt đáy
Cấp 3 Kiểm chứng công trình Báo cáo thiết kế kỹ thuật Đà Nẵng Plaza & Cục Hải quan ĐN Đánh giá sai số thực tế, chuẩn hóa ứng dụng

Quy trình nghiên cứu và Triangulation

Quy trình nghiên cứu kiểm định chéo tam giác (Triangulation) nghiêm ngặt giữa 4 nguồn dữ liệu:

  1. Kết quả tính toán giải tích lý thuyết giải từ phương trình tương tác khung – vách.
  2. Kết quả phân tích động lực học trực tiếp bằng phương pháp phần tử hữu hạn theo tiêu chuẩn TCVN 2737:1995 (chuẩn so sánh $W_{dk}^{TCVN}$).
  3. Kết quả tính toán theo công thức gần đúng hiện hành của Điều 5.5 TCVN 2737:1995.
  4. Kết quả tính toán từ công thức đề xuất mới của luận án ($W_{dk}^{DX}$ và phương pháp $GLF$).
                  ┌──────────────────────────────────────────────┐
                  │          TAM GIÁC KIỂM CHỨNG KẾT QUẢ         │
                  └──────────────────────┬───────────────────────┘
                                         │
                 ┌───────────────────────┼───────────────────────┐
                 ▼                       ▼                       ▼
   ┌───────────────────────────┐ ┌───────────────┐ ┌───────────────────────────┐
   │ 1. Giải tích lý thuyết    │ │ 2. Mô hình số │ │ 3. Công thức đề xuất mới  │
   │    Phương trình vi phân   │ │    FEM ETABS  │ │    Hàm K1(αH) & GLF       │
   │    chuyển vị khung-vách   │ │    SAP2000    │ │    độ lệch < 3.8%         │
   └─────────────┬─────────────┘ └───────┬───────┘ └─────────────┬─────────────┘
                 │                       │                       │
                 └───────────────────────┼───────────────────────┘
                                         ▼
                         ┌───────────────────────────────┐
                         │   4. KIỂM CHỨNG THỰC TẾ       │
                         │   Đà Nẵng Plaza (20 tầng)     │
                         │   Cục Hải quan TP. ĐN (22 t.) │
                         └───────────────────────────────┘

Dữ liệu và phân tích tham số

Bộ mẫu nghiên cứu tham số được thiết kế chi tiết bao gồm 4 nhóm mặt bằng kết cấu điển hình (Dạng 1, 2, 3, 4) và các phân nhóm mở rộng (Dạng 1a, 2a, 3a, 4, 5, 6):

  • Số tầng: Khảo sát tại các mốc 20 tầng ($H = 68.0\text{m}$), 25 tầng ($H = 85.0\text{m}$), 30 tầng ($H = 102.0\text{m}$), và 35 tầng ($H = 119.0\text{m}$). Chiều cao tầng điển hình $h = 3.4\text{m}$.
  • Chiều dày vách cứng ($t_w$): Biến thiên từ $200\text{mm}, 250\text{mm}, 300\text{mm}$ đến $350\text{mm}$ để quét toàn bộ dải biến thiên của thông số độ cứng tương đối $\alpha H \in [0.5, 6.0]$.
  • Đặc trưng tải trọng gió: Vùng gió II-B (Thành phố Đà Nẵng) với áp lực gió cơ sở $W_0 = 95\text{ daN/m}^2$, hệ số độ cao tính toán theo 3 dạng địa hình A, B, C theo TCVN 2737:1995.

Phân tích định lượng độ tin cậy được thực hiện qua việc tính toán sai số tương đối phần trăm tại từng tầng $i$ và sai số lực cắt đáy tổng thể $Q_0$: $$\Delta W_{di} = \frac{W_{di}^{DX} - W_{di}^{TCVN}}{W_{di}^{TCVN}} \times 100% ; \quad \Delta Q_0 = \frac{Q_0^{DX} - Q_0^{TCVN}}{Q_0^{TCVN}} \times 100%$$

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

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

Luận án đã chứng minh bằng số liệu thực nghiệm mô phỏng 5 phát hiện mang tính bước ngoặt:

┌───────────────────────────────────────────────────────────────────────────────────┐
│               SO SÁNH SAI SỐ LỰC CẮT ĐÁY TẢI TRỌNG GIÓ ĐỘNG                       │
├─────────────────────────┬───────────────────────────────┬─────────────────────────┤
│    Dạng kết cấu / Tầng  │ Công thức gần đúng TCVN 2737  │   Công thức đề xuất     │
├─────────────────────────┼───────────────────────────────┼─────────────────────────┤
│ Nhà 20 tầng (Dạng 1)    │ Sai số: -14.82% đến +18.65%   │ Sai số: -1.24% đến +2.15%│
│ Nhà 30 tầng (Dạng 1)    │ Sai số: -22.40% đến +28.90%   │ Sai số: -1.85% đến +2.80%│
│ Nhà 20 tầng (Dạng 2-4)  │ Sai số cục bộ tầng mái: >35%  │ Sai số cục bộ: < 4.50%  │
│ Nhà 30 tầng (Dạng 2-4)  │ Sai số cục bộ tầng mái: >48%  │ Sai số cục bộ: < 5.00%  │
│ Đà Nẵng Plaza (20 tầng) │ Sai số lực cắt đáy: -12.60%   │ Sai số lực cắt đáy: +1.82%│
│ Cục Hải quan ĐN (22 t.) │ Sai số lực cắt đáy: -15.40%   │ Sai số lực cắt đáy: +2.14%│
└─────────────────────────┴───────────────────────────────┴─────────────────────────┘
  1. Vạch rõ giới hạn sai số nghiêm trọng của TCVN 2737:1995: Luận án chứng minh rằng công thức gần đúng hiện hành (1.50) trong TCVN 2737:1995 gây ra sai số cục bộ từ $25%$ đến $48%$ tại các tầng trên cao và sai số lực cắt đáy từ $12%$ đến $23%$ đối với các công trình có vách cứng tham gia chịu lực lớn ($\alpha H < 2.5$).
  2. Quy luật phi tuyến của hàm phân bố thành phần động $K_1$: Dạng phân bố của thành phần gió động phụ thuộc phi tuyến tính vào tham số $\alpha H$. Khi vách càng dày (tính uốn tăng, $\alpha H$ nhỏ), tải trọng gió động dịch chuyển mạnh lên phía đỉnh công trình. Khi khung chiếm ưu thế ($\alpha H > 6.0$), quy luật phân bố mới tiệm cận về đường tuyến tính bậc nhất.
  3. Hiệu chỉnh giải tích thành công cho TCVN 2737: Công thức cải tiến đề xuất tính thành phần gió động dựa trên hàm $K_1(\alpha H, \bar{z})$ đã triệt tiêu sai số có hệ thống, đưa sai số lực cắt đáy trên tất cả 24 mô hình khảo sát 20–35 tầng xuống dưới mức $\pm 2.8%$.
  4. Thiết lập công thức giải tích hệ số gió giật $G$ cho Việt Nam: Luận án đã đóng gói thành công phương pháp GLF Davenport thành công thức đơn giản hóa áp dụng trực tiếp cho vùng gió II-B:
    • Đối với công trình có $\alpha H \le 2.0$: $$G = 1 + \frac{\xi \cdot \varepsilon \cdot \nu}{\bar{q}_0} \cdot \left[ 0.75 + 0.25 \left(\frac{z}{H}\right)^2 \right]$$
    • Đối với công trình có $\alpha H = 2.0 \div 6.0$: $$G = 1 + \frac{\xi \cdot \varepsilon \cdot \nu}{\bar{q}_0} \cdot \left[ 0.60 + 0.40 \left(\frac{z}{H}\right) \right]$$
  5. Kiểm chứng thành công trên công trình thực tế: Tại tòa nhà Đà Nẵng Plaza (20 tầng) và Cục Hải quan TP. Đà Nẵng (22 tầng), công thức đề xuất đạt độ tương đồng vượt trội với kết quả phân tích động lực học chi tiết: chênh lệch lực cắt đáy chỉ lần lượt là $+1.82%$ và $+2.14%$, trong khi công thức TCVN 2737 cũ đánh giá non tải trọng tương ứng $-12.60%$ và $-15.40%$.

Implications đa chiều

  • Về mặt học thuật: Cung cấp mô hình giải tích hợp nhất giữa biến dạng uốn – cắt và dao động ngẫu nhiên do gió giật, lấp đầy khoảng trống chuyển giao công nghệ giữa tiêu chuẩn SNiP và ASCE/Eurocode.
  • Về mặt kỹ thuật thực hành: Cho phép kỹ sư kết cấu tính toán nhanh chóng, chính xác tải trọng gió động và hệ số gió giật mà không bắt buộc phải chạy các bài toán phổ phản ứng dao động phức tạp trên phần mềm thương mại.
  • Về mặt chính sách tiêu chuẩn: Đóng góp cơ sở khoa học và công thức toán học trực tiếp cho Ban soạn thảo Tiêu chuẩn Quốc gia trong quá trình biên soạn, sửa đổi và nâng cấp tiêu chuẩn thay thế TCVN 2737:1995.

Limitations và Future Research

Nghiên cứu có các giới hạn biên cần được nhìn nhận khách quan:

  • Giới hạn hình học và đối xứng: Công thức chỉ áp dụng cho công trình có mặt bằng đối xứng hoặc gần đối xứng, bỏ qua tương tác xoắn – uốn không gian (3D coupled torsional-lateral vibration).
  • Giới hạn chiều cao: Phạm vi kiểm chứng tối ưu cho công trình từ 15 đến 35 tầng ($H \le 120\text{m}$). Với các công trình siêu cao tầng ($H > 150\text{m}$), ảnh hưởng của các dạng dao động bậc cao (bậc 2, bậc 3) và dao động xoáy ngang luồng gió (Cross-wind vortex shedding) trở nên áp đảo.
  • Giới hạn dữ liệu khí tượng: Các tham số xung vận tốc gió và phổ mật độ Davenport được hiệu chuẩn tập trung cho vùng gió II-B (Đà Nẵng và ven biển Miền Trung).

Chương trình nghiên cứu 10 năm tiếp theo (Future Research Agenda):

  1. Mở rộng mô hình giải tích cho kết cấu khung – vách – lõi bất đối xứng chịu dao động xoắn ghép.
  2. Ứng dụng mô hình phần tử hữu hạn nâng cao kết hợp thí nghiệm hầm gió (Wind Tunnel Testing) sử dụng công nghệ cân lực đa thành phần tần số cao (HFFB) và mảng cảm biến áp lực đồng bộ (SPAN).
  3. Nghiên cứu ảnh hưởng của hiện tượng tương tác khí đàn hồi (Aeroelasticity) và hiệu ứng che chắn, giao thoa khí động học (Interference effects) trong cụm nhà cao tầng đô thị.
  4. Xây dựng bản đồ hệ số gió giật số hóa chi tiết cho toàn bộ các vùng gió I, II, III, IV, V trên lãnh thổ Việt Nam.

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

Luận án tạo ra tác động sâu rộng trên cả 4 phương diện:

  • Tác động học thuật (Academic Impact): Tạo tiền đề cho các công trình nghiên cứu tiếp theo tại Việt Nam về động lực học công trình ngẫu nhiên, ước tính thu hút hàng trăm trích dẫn trong các tạp chí chuyên ngành kết cấu, xây dựng và cơ học kỹ thuật.
  • Chuyển đổi ngành công nghiệp xây dựng (Industry Transformation): Nâng cao chất lượng thiết kế kết cấu nhà cao tầng tại các công ty tư vấn thiết kế hàng đầu. Giúp tối ưu hóa hàm lượng cốt thép và tiết diện vách cứng, giảm từ $8%$ đến $12%$ chi phí vật liệu kết cấu chịu lực ngang không cần thiết do thiên lệch sai số tính toán cũ.
  • Ảnh hưởng chính sách tiêu chuẩn (Policy Influence): Cung cấp luận cứ cốt lõi cho Bộ Xây dựng và Viện Khoa học Công nghệ Xây dựng (IBST) trong việc hiện đại hóa quy chuẩn xây dựng Việt Nam.
  • Lợi ích an toàn xã hội (Societal Benefits): Đảm bảo an toàn chịu lực tuyệt đối cho hàng triệu người dân sinh sống và làm việc trong các tòa nhà cao tầng dọc dải bờ biển Miền Trung thường xuyên chịu bão mạnh.

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

  • Nghiên cứu sinh và Giảng viên ngành Kỹ thuật Xây dựng / Cơ học: Tiếp cận phương pháp tích hợp dầm liên tục tương đương và cơ chế lan truyền năng lượng gió giật Davenport.
  • Kỹ sư kết cấu thực hành: Sở hữu công cụ tính tay và bảng biểu tra cứu nhanh hệ số $K_1$ và hệ số gió giật $G$, tăng tốc độ thiết kế sơ bộ và thẩm tra kết cấu.
  • Các Chủ đầu tư bất động sản: Tối ưu hóa chi phí đầu tư xây dựng phần thô kết cấu nhà cao tầng mà vẫn đảm bảo an toàn tuyệt đối trước gió bão.
  • Cơ quan quản lý nhà nước và Hội đồng nghiệm thu: Có căn cứ khoa học chuẩn xác để đánh giá, kiểm tra an toàn kết cấu công trình cao tầng chịu tải trọng gió bão.

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ì?

Đóng góp độc đáo nhất là việc thiết lập thành công hàm phân bố giải tích phi tuyến $K_1(\alpha H, \bar{z})$ dựa trên việc giải phương trình vi phân tương tác khung – vách môi trường liên tục (Stafford Smith & Coull), giải quyết triệt để sự mâu thuẫn giữa giả thiết dạng dao động tuyến tính trong TCVN 2737:1995 và bản chất biến dạng uốn – cắt phi tuyến của kết cấu thực tế.

2. Đổi mới phương pháp luận so với các nghiên cứu trước đây?

So với phương pháp gần đúng của TCVN 2737:1995 (vốn giả định $y = z/H$) và phương pháp số thuần túy (đòi hỏi mô hình hóa phức tạp và phân tích phổ dao động riêng), luận án đã kết hợp phương pháp giải tích giải chính xác phương trình dầm tương đương với kỹ thuật hồi quy tham số trên mô hình số phần tử hữu hạn 3D (ETABS/SAP2000), tạo ra công thức thực hành vừa đạt độ chính xác tương đương phân tích động lực học chi tiết (sai số $< 3.8%$), vừa duy trì dạng công thức đại số tường minh, dễ áp dụng.

3. Phát hiện gây ngạc nhiên nhất về mặt dữ liệu là gì?

Phát hiện gây ngạc nhiên nhất là công thức hiện hành của Điều 5.5 TCVN 2737:1995 không chỉ gây sai số lớn mà còn đánh giá non tải trọng (underestimate) lực cắt đáy từ $12.6%$ đến $15.4%$ ở các công trình thực tế có vách cứng lớn (như Đà Nẵng Plaza và Cục Hải quan TP. Đà Nẵng). Điều này đồng nghĩa với việc nhiều công trình thực tế thiết kế theo công thức gần đúng cũ có thể đang thiếu an toàn chịu lực trước tải trọng gió bão cực hạn.

4. Quy trình tái lập nghiên cứu (Replication Protocol) có được cung cấp đầy đủ không?

Luận án cung cấp tường minh toàn bộ hệ phương trình vi phân, các bước biến đổi toán học, bảng số liệu tọa độ tiết diện cột/dầm/vách của 24 mô hình kiểm chuẩn, dữ liệu khí động học vùng gió II-B, và hồ sơ thiết kế chi tiết của 2 công trình thực nghiệm kiểm chứng, đảm bảo khả năng tái lập và kiểm chứng độc lập hoàn toàn.

5. Kế hoạch nghiên cứu 10 năm được vạch ra như thế nào?

Kế hoạch 10 năm tập trung vào 3 hướng chính: (1) Mở rộng công thức cho kết cấu lõi vách bất đối xứng có xét dao động xoắn ghép 3D; (2) Thực nghiệm hầm gió đối chứng tương tác khí đàn hồi cho nhà siêu cao tầng $H > 150\text{m}$; (3) Xây dựng phần mềm tích hợp tự động tính toán hệ số gió giật GLF theo bản đồ số hóa bão biển Việt Nam.

Kết luận

  1. Lượng hóa chính xác sai số của TCVN 2737:1995: Chỉ rõ công thức gần đúng cũ gây sai số cục bộ lên đến $48%$ trên các tầng cao và đánh giá thiếu an toàn lực cắt đáy từ $12%$ đến $15%$ đối với hệ kết cấu khung giằng.
  2. Đề xuất hàm phân bố động lực học $K_1$ mới: Phát triển thành công hàm giải tích $K_1(\alpha H, \bar{z})$ phản ánh trung thực bản chất tương tác uốn – cắt, giảm sai số lực cắt đáy xuống dưới $2.8%$.
  3. Hiện đại hóa phương pháp GLF Davenport tại Việt Nam: Thiết lập hệ thống công thức giải tích tính hệ số gió giật $G$ tương thích hoàn hảo với cơ sở dữ liệu khí quyển và tiêu chuẩn xây dựng Việt Nam cho vùng gió II-B.
  4. Kiểm chứng thực nghiệm đa cấp độ thành công: Khảo sát toàn diện trên 24 mô hình tham số và 2 công trình thực tế (Đà Nẵng Plaza 20 tầng, Cục Hải quan TP. Đà Nẵng 22 tầng) với độ tin cậy vượt trội (độ lệch $< 3.8%$).
  5. Đóng góp nền tảng cho tiêu chuẩn quốc gia: Cung cấp cơ sở khoa học trực tiếp phục vụ biên soạn các thế hệ Tiêu chuẩn Tải trọng và Tác động mới thay thế TCVN 2737:1995.
  6. Mở ra các hướng nghiên cứu chuyển tiếp: Khởi xướng các nghiên cứu sâu rộng về dao động xoắn không gian, tương tác khí đàn hồi và hầm gió khí động học cho nhà siêu cao tầng tại Việt Nam.