Tổng quan về luận án

Nghiên cứu về độ an toàn sinh mạng và độ bền vững cơ học của các công trình biển quảng cáo tấm lớn (large-panel billboards/hoardings) đặt ngoài trời là một bài toán kỹ thuật gió (Wind Engineering) mang tính thời sự cấp thiết. Tại Việt Nam, hệ thống pano quảng cáo độc lập dạng một cột đỡ phân bố dày đặc dọc theo các trục đường cao tốc và khu vực ngoại vi đô thị. Với diện tích đón gió đặc trưng từ $40\text{ m}^2$ đến hơn $200\text{ m}^2$, chiều dày kết cấu rất nhỏ so với kích thước mặt bảng ($b/c$ dao động từ $1,0$ đến $5,0$), công trình mang tính thanh mảnh cao và chịu rủi ro rất lớn trước tác động của bão nhiệt đới. Thực tế ghi nhận hàng loạt sự cố thảm khốc: vụ gãy đổ cột bảng quảng cáo đường kính $2,0\text{ m}$ tại TP. Thủ Đức ngày 19/7/2021, cùng nhiều sự cố phá hoại liên kết chân cột, xé rách tấm bọc, và vặn xoắn kết cấu khung giàn tại các tỉnh ven biển miền Trung và miền Bắc khi bão đổ bộ.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ sự bất cập và thiếu đồng nhất trong các tiêu chuẩn thiết kế: tiêu chuẩn cũ TCVN 2737:1995 áp dụng hệ số khí động tĩnh quy ước khá thấp ($c_x = 1,2$) và giả định tải trọng phân bố đều với lực đặt đúng tâm hình học ($e_y = 0$), hoàn toàn bỏ qua hiệu ứng xoắn khí động học ($M_z$). Ngược lại, các tiêu chuẩn quốc tế tiên tiến như EN 1991-1-4:2005 ($C_f = 1,8; e_y = 0,25b$), ASCE/SEI 7-16 ($C_f = 1,7 \div 1,8; e_y = 0,2b$), và AS/NZS 1170.2:2011 quy định hệ số lực cao hơn đáng kể và bắt buộc xét đến độ lệch tâm ngang. Mặc dù tiêu chuẩn mới TCVN 2737:2023 đã nâng hệ số khí động lên $c_x = 1,6 \div 1,7$ và $e_y = 0,25b$, việc thiếu hụt dữ liệu thực nghiệm khí động học trong điều kiện lớp biên khí quyển thực tế tại Việt Nam đặt ra yêu cầu phải kiểm chứng độc lập và chuẩn hóa phương pháp tính.

Luận án của nghiên cứu sinh Nguyễn Lệ Thủy tại Viện Khoa học Công nghệ Xây dựng (IBST) dưới sự hướng dẫn của PGS. TS Nguyễn Hồng Sơn và TS. Vũ Thành Trung đã giải quyết bài toán này thông qua hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  1. Câu hỏi nghiên cứu 1 (RQ1): Sự khác biệt định lượng về lực cắt đáy ($F_x$), mô men uốn lật ($M_y$), và mô men xoắn đáy ($M_z$) khi tính toán tải trọng gió lên pano tấm lớn giữa TCVN 2737:1995, TCVN 2737:2023, EN 1991-1-4, và ASCE/SEI 7-16 là bao nhiêu?
  2. Câu hỏi nghiên cứu 2 (RQ2): Các tham số hình học pano ($b/c, c/h$), tỷ lệ mô hình, mô hình dòng chảy rối (RANS, SST $k-\omega$, RSM), và hàm tường (wall functions) trong mô phỏng số CFD (ANSYS Fluent) phản ánh trường áp suất và hệ số lực gió như thế nào?
  3. Câu hỏi nghiên cứu 3 (RQ3): Dưới tác động của dòng gió rối mô phỏng địa hình Việt Nam trong ống thổi khí động, quy luật biến thiên của hệ số lực pháp tuyến ($C_{f,x}$), hệ số lực ngang ($C_{f,y}$), và độ lệch tâm tương đối ($e_y/b$) theo góc hướng gió ($\theta = 0^\circ \div 90^\circ$) diễn ra như thế nào?

Giả thuyết nghiên cứu: Tải trọng gió thực tế gây ra đồng thời lực cản chính diện cực đại ở góc đón gió trực diện ($\theta = 0^\circ \div 15^\circ$) và mô men xoắn nguy hiểm nhất ở các góc đón gió xiên ($\theta = 30^\circ \div 60^\circ$), đòi hỏi phải hiệu chỉnh đồng thời hệ số lực $C_{f,x}$ và tỷ số lệch tâm $e_y/b$ cho các phân loại pano phổ biến ($3\times6\text{ m}, 4\times8\text{ m}, 5\times14\text{ m}, 6\times18\text{ m}$).

Phạm vi nghiên cứu bao quát kết cấu pano bảng quảng cáo hai mặt phẳng song song bọc kín đặt trên một cột đỡ ống thép độc lập ngoài trời, chịu tác động của 5 vùng áp lực gió theo QCVN 02:2022/BXD với vận tốc gió bão cực hạn đạt cấp 13–14 ($37 \div 46\text{ m/s}$) và gió giật cấp 16–17 ($61 \div 67\text{ m/s}$).

Literature Review và Positioning

Nghiên cứu khí động học công trình bảng quảng cáo trên thế giới phân bổ theo hai dòng chính: đo áp suất bề mặt tĩnh/động đa điểm (pressure tap measurement) và đo cân bằng lực tần số cao (High-Frequency Force Balance - HFFB).

Dòng nghiên cứu lý thuyết và thực nghiệm kinh điển của Letchford (2001) trên các tấm tường và bảng hiệu đã chứng minh hệ số lực cản $C_f$ không cố định mà biến thiên phức tạp theo tỷ số cạnh $b/c$ và tỷ số độ cao $c/h$, qua đó đề xuất công thức bán thực nghiệm: $$C_f = 1,45 + 0,5 \times 0,7 \log_{10}(b/c) \times (c/h - 0,5)$$ với sai số $\pm 0,1%$ trong miền $0,2 < b/c < 5,0$. Tuy nhiên, dữ liệu của Letchford chỉ ra tính bất thường: khi $b/c = 1,0$, hệ số khí động giảm khi $c/h > 0,67$ nhưng lại tăng khi $c/h < 0,3$, cho thấy cấu trúc dòng chảy cuộn xoáy mép dưới pano khi đặt sát mặt đất có sự chuyển pha phức tạp.

Trong khi đó, Delong Zuo và cộng sự (2012, 2014) thực hiện cả thí nghiệm hiện trường tỷ lệ 1:1 tại Lubbock (Texas) và thí nghiệm hầm gió tỷ lệ 1:50 trên 4 dạng pano (tấm đơn, tấm đôi song song, tấm chữ V, và dạng hộp). Kết quả công bố hệ số lực cản chính diện lần lượt là $C_x = 1,41$ (tấm đơn), $1,35$ (tấm đôi), $1,44$ (chữ V), và $1,23$ (hộp); cùng hệ số xoắn $C_e$ tương ứng là $0,19; 0,20; 0,21$ và $0,12$. Phát hiện của Delong Zuo khẳng định mô men xoắn thực tế đo được vượt xa các quy định trong các tiêu chuẩn thiết kế hiện hành. Douglas và cộng sự (2014) phân tích 470 bản ghi hầm gió và hiện trường, chứng minh giá trị cực trị của tải trọng gió tuân theo hàm phân phối xác suất Fisher-Tippett Loại I (Gumbel), đồng thời chỉ ra hệ số lực trung bình thu được từ mô hình thu nhỏ trong hầm gió có xu hướng cao hơn từ 5–10% so với mô hình thực địa do hiệu ứng số Reynolds ($Re$) và vi mô trường dòng.

Đối lập với quan điểm của Delong Zuo, nghiên cứu của Wang và cộng sự (2016) trên mô hình tỷ lệ 1:20 cho các pano hai mặt và ba mặt dạng hộp lại cho rằng mô men xoắn tổng thể có biên độ thấp hơn quy định tiêu chuẩn và đóng góp không đáng kể vào tổng ứng suất phá hoại cột đỡ khi góc đón gió thay đổi. Tồn tại một cuộc tranh luận học thuật sâu sắc giữa trường phái xem mô men xoắn khí động là yếu tố thứ yếu (Wang et al., 2016) và trường phái khẳng định mô men xoắn kết hợp với độ lệch tâm ngang $e_y = (0,20 \div 0,25)b$ là nguyên nhân phá hủy bu lông chân cột chính (Warnitchai, 2009; Meyer et al., 2017; Samith Buddika et al., 2022).

Nghiên cứu của Warnitchai (2009) tại Viện Công nghệ Châu Á (AIT) bằng thiết bị HFFB trên mô hình gỗ pano hai mặt xác định độ lệch tâm ngang $e = 0,1b$ tại góc $\theta = 0^\circ$ và vọt lên $e = 0,2b$ tại $\theta = 45^\circ$. Gần đây, Meyer và cộng sự (2017) thí nghiệm 11 mẫu pano tỷ lệ 1:3 tại Đại học Quốc tế Florida (FIU) với các tỷ số $b/c = 1, 3, 5$ và $d/c = 0,1; 0,4; 0,7$, chỉ ra hệ số xoắn $C_{M,z}$ đạt từ $0,02$ đến $0,15$, tương ứng với tỷ số độ lệch tâm $r/b$ dao động từ $0,02$ đến $0,13$. Samith Buddika và cộng sự (2022) mô phỏng số 09 mô hình pano hai mặt bằng mô hình rối SST $k-\omega$, thiết lập phương trình hồi quy hệ số lực: $$C_{Ft} = 1,787 - 0,04AR - 0,228CR - 1,246\cos\theta + 1,214\cos^2\theta \quad (\text{với } 0^\circ < \theta < 45^\circ)$$ và kiến nghị độ lệch tâm $e_x = 0,18b$ đến $0,21b$.

Vị thế của luận án được xác lập rõ ràng: Đây là công trình đầu tiên tại Việt Nam tích hợp toàn diện giữa phân tích giải tích đối chuẩn 4 hệ tiêu chuẩn, khảo sát tham số khí động số vi mô (CFD ANSYS Fluent với UDF profile gió bão), và thực nghiệm hầm gió lớp biên khí quyển tỷ lệ lớn trên cảm biến đa trục 6-DOF HFFB JR3. Luận án giải quyết triệt để bất đồng học thuật về hệ số lực cản và độ lệch tâm cho dạng pano hai mặt bọc kín trong điều kiện địa hình và bão nhiệt đới đặc thù của 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 tạo ra bước tiến quan trọng trong lý thuyết khí động học công trình bằng việc mở rộng và tinh chỉnh lý thuyết lớp biên khí quyển (Atmospheric Boundary Layer - ABL) và lý thuyết giá trị cực trị (Extreme Value Theory) áp dụng cho kết cấu cản gió dạng tấm phẳng đặt trên cao:

  1. Mở rộng lý thuyết phân bố cực trị Gumbel (Fisher-Tippett Type I) trong chuyển đổi chu kỳ lặp và thời gian lấy mẫu: Thiết lập hệ phương trình chuyển đổi tương thích giữa vận tốc gió giật 3 giây ($V_{3s}$ theo ASCE 7-16, TCVN 2737:2023) và vận tốc gió trung bình 10 phút / 1 giờ ($V_{10min}, V_{1hr}$ theo EN 1991-1-4, ISO 4354:2009): $$U_R = u - a \ln\left[-\ln\left(1 - \frac{1}{R}\right)\right] \approx u + a \ln(R)$$ Chứng minh tính phi tuyến của hệ số chuyển đổi vận tốc khi tính toán phản ứng động lực học của kết cấu pano.

  2. Bổ khuyết mô hình trường dòng xoáy phân ly và tái bám (Separation and Reattachment Flow Model): Chứng minh đối với pano tấm mỏng ($d/c < 0,1$), hiện tượng xoáy tách rời ở mép đón gió tạo ra trường áp suất âm cực đại không đối xứng ở mặt sau khi gió thổi xiên ($\theta = 15^\circ \div 45^\circ$), làm dịch chuyển tâm áp lực (Center of Pressure - CoP) ra xa trục đối xứng hình học, sinh ra hiệu ứng ngẫu lực xoắn $M_z$ mà lý thuyết tĩnh học cổ điển không thể mô tả.

                  DÒNG GIÓ TỚI (Vận tốc U, Góc tới θ)

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

Khung phân tích của luận án tích hợp đồng thời 3 trụ cột lý thuyết: (i) Cơ học chất lưu tính toán (Computational Fluid Dynamics); (ii) Lý thuyết tương tác gió - công trình (Bluff Body Aerodynamics); và (iii) Cơ học kết cấu phi tuyến chịu tải trọng ngẫu nhiên.

Phương pháp luận phân tích: Thiết lập chu trình kiểm tra chéo (cross-validation matrix) 3 cấp độ:

  • Cấp 1: Phân tích giải tích chuẩn tắc theo các tiêu chuẩn TCVN 2737:1995, TCVN 2737:2023, EN 1991-1-4, ASCE/SEI 7-16.
  • Cấp 2: Thử nghiệm số (Numerical Experimentation) bằng phương pháp thể tích hữu hạn (FVM) khảo sát các biến thiên $b/c$ từ $1,0$ đến $5,0$, kiểm soát độ nhạy lưới qua các chỉ số chất lượng trực giao (Orthogonal Quality $> 0,85$) và độ lệch lưới (Skewness $< 0,25$).
  • Cấp 3: Thực nghiệm vật lý trong hầm gió đóng kín (Closed-loop wind tunnel) với hệ thống mô phỏng độ rối nhân tạo (spires, barrier wall, roughness cubes) tạo lập chính xác profile vận tốc gió dạng lũy thừa theo số mũ Hellman $m = 0,27$ (dạng địa hình B).

Điều kiện biên xác lập: Áp dụng cho kết cấu bảng quảng cáo độc lập đặt ngoài trời có tỷ số cạnh $b/c = 1,0 \div 5,0$; tỷ số độ cao $c/h = 0,2 \div 0,7$; chiều dày tấm $d/c \le 0,1$; không xét đến hiện tượng biến dạng động khí đàn hồi tức thời (mô hình cứng - rigid model).

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 hiện thực phản biện (Critical Realism), sử dụng phương pháp hỗn hợp định lượng bậc cao (Advanced Quantitative Multi-Method Design) với sự phối hợp chặt chẽ giữa mô phỏng số và thực nghiệm vật lý hầm gió.

       THIẾT KẾ NGHIÊN CỨU THỰC NGHIỆM VÀ SỐ HÓA ĐA CẤP ĐỘ

  [PHÂN TÍCH TIÊU CHUẨN]         [MÔ PHỎNG SỐ CFD]             [THỰC NGHIỆM HẦM GIÓ]
  - TCVN 2737:1995               - ANSYS Fluent (FVM)          - Ống thổi khí động IBST
  - TCVN 2737:2023               - SST k-ω, RSM, RANS          - Bàn xoay 360°, Spire, Roughness
  - EN 1991-1-4                  - Lưới hexahedral/tetrahedral - Cảm biến 6-DOF JR3 (1000 Hz)
  - ASCE/SEI 7-16                - UDF Profile gió bão         - 6 Mẫu thử nghiệm (M1 - M6)
                       [ĐỐI CHUẨN VÀ TỐI ƯU HÓA THIẾT KẾ]
                       - Hệ số lực chuẩn hóa: Cf,x = 1,6 ÷ 1,7
                       - Hệ số lệch tâm: ey/b = 0,25
                       - Bảng tính tự động hóa & Hướng dẫn kỹ thuật

Quy trình nghiên cứu rigorous

Giai đoạn 1: Mô phỏng số CFD trên phần mềm ANSYS Fluent:

  • Không gian tính toán (Computational Domain): Thiết lập hầm gió ảo 3D với kích thước tiêu chuẩn nhằm loại bỏ hiệu ứng nghẽn dòng (Blockage ratio $< 3%$). Miền tính toán có khoảng cách từ lối vào (Inlet) đến mặt pano là $5H$, từ mặt pano đến lối ra (Outlet) là $10H \div 15H$, chiều rộng $10B$ và chiều cao $5H$ (với $H, B$ là chiều cao và chiều rộng tổng thể mô hình).
  • Lập trình hàm người dùng UDF (User-Defined Function): Định nghĩa chính xác profile vận tốc gió trung bình $V(z)$ và phân bố động năng rối $k(z)$, tỷ lệ tiêu tán rối $\varepsilon(z)$ / $\omega(z)$ theo quy luật lớp biên khí quyển thực tế.
  • Khảo sát mô hình dòng rối: So sánh nghiêm ngặt giữa các mô hình $k-\varepsilon$ (Standard, RNG, Realizable), Shear-Stress Transport (SST) $k-\omega$, Reynolds Stress Model (RSM), và đánh giá hàm tường tiêu chuẩn (Standard Wall Functions) so với hàm tường tăng cường (Enhanced Wall Treatment) với giá trị $y^+ \approx 1$.

Giai đoạn 2: Thiết lập thực nghiệm trong ống thổi khí động IBST:

  • Cấu trúc hầm gió: Ống thổi khí động tuần hoàn dạng kín theo phương đứng, kích thước phòng thí nghiệm $8,0\text{ m} \times 35,0\text{ m} \times 9,0\text{ m}$. Đoạn thử nghiệm nhánh dưới có tiết diện $2,0\text{ m} \times 2,0\text{ m}$ tích hợp bàn xoay đường kính $4,0\text{ m}$ (vận tốc dòng $0 \div 50\text{ m/s}$).
  • Thiết bị đo lường lực 6 thành phần JR3 (High-Frequency Force Balance - Model 30E12A4): Gắn tại gốc chân cột của mô hình, dải đo tiêu chuẩn $F_X, F_Y = \pm 200\text{ N}$; $F_Z = \pm 400\text{ N}$; $M_X, M_Y, M_Z = \pm 16\text{ Nm}$. Tần số lấy mẫu đạt $1000\text{ Hz}$ (1000 giá trị/giây), độ phân giải kỹ thuật số lực là $0,025\text{ N}$ và mô men là $0,002\text{ Nm}$, cấp chính xác danh nghĩa trên tất cả các trục đạt $0,25%$.
  • Thiết lập mô hình vật lý: Chế tạo 06 mẫu thử nghiệm thu nhỏ (ký hiệu M1 đến M6) bằng vật liệu cứng (hợp kim nhôm và mica chuyên dụng) với các tỷ lệ hình học $b/c$ và chiều cao đặt pano $h$ biến thiên đại diện cho toàn bộ các dạng pano thực tế. Thiết lập hàng rào chắn (barrier), thanh tạo xoáy nhọn (spires), và thảm cục tạo nhám (roughness cubes) bằng gỗ để tạo dựng trường gió địa hình B có profile vận tốc và phổ mật độ công suất rối (PSD) đồng dạng hoàn toàn với tự nhiên.

Data và phân tích

Ma trận thực nghiệm và mô phỏng được thực hiện đồng bộ trên 6 mẫu pano cứng đại diện:

Mẫu thử Kích thước pano nguyên mẫu $(b \times c)$ [m] Chiều cao đỉnh pano nguyên mẫu $h$ [m] Tỷ số cạnh $b/c$ Tỷ lệ mô hình thu nhỏ Vận tốc gió thí nghiệm $U_{ref}$ [m/s] Tần số quét HFFB [Hz]
M1 $6,0 \times 3,0$ $12,0$ $2,0$ 1:10 / 1:30 15,0 – 35,0 1000
M2 $8,0 \times 4,0$ $14,0$ $2,0$ 1:10 / 1:30 15,0 – 35,0 1000
M3 $14,0 \times 5,0$ $16,0$ $2,8$ 1:10 / 1:30 15,0 – 35,0 1000
M4 $18,0 \times 6,0$ $18,0$ $3,0$ 1:10 / 1:30 15,0 – 35,0 1000
M5 $18,0 \times 6,0$ $21,0$ $3,0$ 1:10 / 1:30 15,0 – 35,0 1000
M6 $15,0 \times 3,0$ $12,0$ $5,0$ 1:10 / 1:30 15,0 – 35,0 1000

Dữ liệu chuỗi thời gian 6 thành phần phản lực đáy $[F_X(t), F_Y(t), F_Z(t), M_X(t), M_Y(t), M_Z(t)]$ được xử lý thống kê loại trừ nhiễu bằng thuật toán lọc thông thấp (low-pass filter), tính toán giá trị trung bình (mean), độ lệch chuẩn (RMS), giá trị đỉnh (peak values) với hệ số đỉnh $g = 3,5$, từ đó trích xuất hệ số lực khí động: $$C_{f,x} = \frac{F_X}{0,5 \rho U_{ref}^2 A_s}; \quad C_{f,y} = \frac{F_Y}{0,5 \rho U_{ref}^2 A_s}; \quad \frac{e_y}{b} = \frac{M_Z}{F_X \cdot b}$$

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

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

  1. Sự thiếu hụt nguy hiểm của TCVN 2737:1995: Tính toán đối chuẩn trên nguyên mẫu pano tiêu chuẩn kích thước $14\times5\text{ m}$ đặt ở độ cao đỉnh $16\text{ m}$ chỉ ra: lực cắt đáy $F_x$ theo TCVN 2737:1995 thấp hơn $24,6%$ so với TCVN 2737:2023, thấp hơn $31,2%$ so với EN 1991-1-4, và thấp hơn $28,5%$ so với ASCE/SEI 7-16. Đặc biệt, mô men lật $M_y$ và mô men xoắn $M_z$ theo tiêu chuẩn cũ bị đánh giá thấp từ $35%$ đến $100%$ (do tiêu chuẩn cũ mặc định $M_z = 0$). Điều này lý giải trực tiếp nguyên nhân hàng loạt vụ sụp đổ, gãy chân cột bu lông trong thực tế khi gặp bão cấp 11–12.

  2. Quy luật phụ thuộc góc hướng gió của hệ số lực pháp tuyến và lực ngang: Kết quả thực nghiệm hầm gió trên cảm biến JR3 xác nhận lực cản trực diện $F_x$ và hệ số $C_{f,x}$ đạt cực đại tại góc gió vuông góc $\theta = 0^\circ \div 15^\circ$ với giá trị thực nghiệm trung bình $C_{f,x} = 1,62 \div 1,68$. Khi góc đón gió tăng dần lên $\theta = 45^\circ \div 60^\circ$, hệ số $C_{f,x}$ giảm dần xuống $1,21 \div 1,01$, nhưng lực ngang $F_y$ và hệ số $C_{f,y}$ tăng vọt từ $0,00$ lên giá trị đỉnh $0,151 \div 0,287$, tạo ra trạng thái chịu lực uốn hai phương kết hợp xoắn cực kỳ bất lợi cho cột thép tròn.

  3. Độ lệch tâm khí động thực nghiệm $e_y/b$: Dữ liệu thực nghiệm chứng minh hợp lực gió không bao giờ đặt tại tâm hình học của pano ngay cả khi gió thổi vuông góc ($\theta = 0^\circ$, độ lệch tâm tồn tại do tính bất đối xứng tức thời của dòng rối đạt $e_y/b = 0,05 \div 0,08$). Khi góc hướng gió dịch chuyển sang góc xiên $\theta = 30^\circ \div 52,5^\circ$, độ lệch tâm tương đối $e_y/b$ tăng vọt và đạt cực đại trong khoảng $0,20 \div 0,25$. Phát hiện này bác bỏ hoàn toàn giả định đặt lực đúng tâm và cung cấp bằng chứng thực nghiệm chuẩn xác khẳng định quy định $e_y = 0,25b$ trong TCVN 2737:2023 và EN 1991-1-4 là hoàn toàn phù hợp với bản chất vật lý của dòng khí động.

  4. Tương quan phi tuyến giữa tỷ số cạnh $b/c$ và hệ số khí động: Khi tỷ số cạnh $b/c$ tăng từ $1,0$ lên $2,8$ và $5,0$ (pano trải dài theo phương ngang), hệ số khí động $C_{f,x}$ có xu hướng tăng từ $1,42$ lên $1,68$, do diện tích cản dòng lớn làm giảm hiệu ứng rò rỉ khí động học (3D flow relief) qua các mép biên đứng, tạo vùng áp suất hút chân không rất mạnh và đồng nhất ở mặt sau pano.

  5. Hiệu chuẩn mô phỏng số CFD: Mô hình dòng rối SST $k-\omega$ kết hợp lưới phân tầng đa lớp cho kết quả sai số dưới $4,8%$ so với thực nghiệm hầm gió, vượt trội hoàn toàn so với mô hình chuẩn $k-\varepsilon$ (sai số $> 14,2%$ do đánh giá thấp vùng xoáy bóc tách mép).

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bộ cơ sở dữ liệu thực nghiệm khí động học chuẩn xác cao về trường lực 6 bậc tự do cho kết cấu dạng tấm mỏng đặt trên cao trong dòng rối khí quyển, đóng góp vào lý thuyết khí động học công trình xây dựng vùng bão nhiệt đới Tây Bắc Thái Bình Dương.
  • Về phương pháp luận: Chuẩn hóa quy trình tích hợp mô phỏng số CFD (với UDF biên dòng gió) và thí nghiệm cân bằng lực tần số cao HFFB JR3 cho các phòng thí nghiệm kỹ thuật gió tại Việt Nam.
  • Về thực tiễn thiết kế: Cung cấp bảng tra và công thức tính toán tường minh hệ số lực $C_{f,x} = 1,6 \div 1,7$ kết hợp độ lệch tâm bắt buộc $e_y = 0,25b$, giúp các kỹ sư kết cấu tính toán chính xác kích thước tiết diện thân cột, chiều dày bản mã chân cột, đường kính và chiều sâu neo bu lông cường độ cao, ngăn chặn nguy cơ gãy đổ kết cấu.
  • Về chính sách và tiêu chuẩn: Là luận cứ khoa học trực tiếp phục vụ hoàn thiện các chỉ dẫn kỹ thuật áp dụng Tiêu chuẩn Quốc gia TCVN 2737:2023 và quy chuẩn QCVN 02:2022/BXD của Bộ Xây dựng.

Limitations và Future Research

  1. Mô hình hóa kết cấu cứng (Rigid Model Limitation): Thí nghiệm hầm gió và mô phỏng CFD trong luận án chỉ dừng lại ở mô hình cứng tuyệt đối, đo lực phản lực đáy mà chưa mô phỏng được chuyển vị động khí đàn hồi (Aeroelastic model). Do đó, các hiện tượng dao động tự kích thích nguy hiểm như mất ổn định khí đàn hồi (Galloping) hay hiện tượng tương tác rung giật xoáy (Vortex-Induced Vibration - VIV) ở vận tốc gió siêu cao chưa được định lượng chi tiết.
  2. Phạm vi cấu hình hình học: Luận án tập trung chuyên sâu vào loại pano phẳng hai mặt song song bọc kín đặt trên một cột đơn. Các cấu hình bảng quảng cáo chữ V, lăng trụ 3 mặt, bảng có khoét lỗ thông gió, hoặc bảng quảng cáo đặt trên mái nhà nhiều tầng chưa được khảo sát thực nghiệm toàn diện trong khuôn khổ này.
  3. Đặc trưng dòng gió phức hợp: Thí nghiệm thiết lập theo dòng gió tầng biên đồng nhất địa hình B; chưa xét đến tác động của địa hình phức tạp cục bộ (đồi dốc, khe núi) hoặc hiệu ứng chắn/hút gió (interference/shielding effects) khi có các công trình cao tầng lân cận.

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

  • Thiết lập hệ thống thí nghiệm mô hình khí đàn hồi (aeroelastic model) tích hợp đo cảm biến gia tốc và đo dịch chuyển quang học nhằm xác định chính xác hệ số suy giảm khí động (aerodynamic damping) và dự báo phản ứng cộng hưởng gió giật.
  • Ứng dụng mô phỏng dòng xoáy lớn LES (Large Eddy Simulation) và mô phỏng thích ứng quy mô SAS (Scale-Adaptive Simulation) trên hệ thống tính toán hiệu năng cao (HPC) để bắt trọn cấu trúc xoáy chuyển tiếp tức thời.
  • Mở rộng ngân hàng dữ liệu khí động cho bảng quảng cáo ba mặt tam giác và hệ thống pano lắp đặt trên mái công trình dân dụng.

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

  • Học thuật: Đặt nền móng cho các nghiên cứu chuyên sâu về khí động học công trình biển báo tại các trường đại học khối kỹ thuật; dự kiến tạo ra các công bố quốc tế uy tín trong danh mục ISI/Scopus (Q1/Q2) thuộc lĩnh vực kết cấu và kỹ thuật gió.
  • Ngành công nghiệp xây dựng và quảng cáo ngoài trời: Định hình lại quy trình tính toán thiết kế kết cấu thép quảng cáo cho hơn 500 doanh nghiệp xây lắp và quảng cáo trên toàn quốc, giúp tối ưu hóa chi phí vật liệu thép (giảm $10 \div 15%$ trọng lượng kết cấu không cần thiết nhờ phân bố tải trọng chính xác, đồng thời tăng cường đúng mức độ an toàn liên kết chân cột).
  • Quản lý nhà nước và xã hội: Cung cấp công cụ thẩm định khoa học cho Sở Xây dựng và cơ quan cấp phép quy hoạch đô thị các địa phương, trực tiếp bảo vệ an toàn tính mạng người tham gia giao thông và phương tiện lưu thông dưới chân các đại lộ trong mùa mưa bão.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận tài liệu chuẩn mực về phương pháp tích hợp mô phỏng CFD bậc cao và kỹ thuật đo đạc hầm gió HFFB 6 thành phần lực.
  • Giảng viên và chuyên gia nghiên cứu kết cấu: Cơ sở dữ liệu đối chuẩn thực nghiệm có độ tin cậy cao để giảng dạy chuyên đề tải trọng gió và kỹ thuật công trình ven biển.
  • Kỹ sư thiết kế kết cấu (Structural Engineers): Sở hữu trực tiếp quy trình, thuật toán và bảng tính Excel/UDF chuẩn hóa để áp dụng ngay vào thiết kế hàng ngày theo TCVN 2737:2023.
  • Cơ quan quản lý quy chuẩn xây dựng (Cục Giám định, Viện IBST): Luận cứ thực nghiệm vững chắc để bảo vệ và cập nhật các điều khoản kỹ thuật trong hệ thống tiêu chuẩn quốc gia.

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 mô hình hóa và định lượng hóa thành công bản chất phi đối xứng của trường áp suất khí động lên pano mỏng đặt trên cao trong dòng rối lớp biên khí quyển, mở rộng Lý thuyết Khí động học vật thể cản gió (Bluff Body Aerodynamics). Luận án chứng minh rằng ngay cả khi luồng gió thổi trực diện vuông góc ($\theta = 0^\circ$), sự bất đối xứng tức thời của cấu trúc xoáy phân ly đã sinh ra độ lệch tâm $e_y = (0,05 \div 0,08)b$; và khi góc gió lệch $\theta = 30^\circ \div 52,5^\circ$, hiệu ứng chênh áp hút chân không cực đại ở mặt sau đẩy độ lệch tâm lên $e_y = (0,20 \div 0,25)b$, tạo ra mô men xoắn $M_z$ áp đảo.

2. Đột phá về phương pháp luận so với các nghiên cứu tiền nhiệm quốc tế?

So với các nghiên cứu trước đây (như Letchford (2001) chỉ tập trung đo áp suất trên tấm tường đơn giản; Warnitchai (2009) chỉ đo mô hình gỗ với dải tỷ số cạnh hẹp; Delong Zuo (2014) chủ yếu đo hầm gió hở với cảm biến giới hạn), đột phá của luận án nằm ở quy trình thực nghiệm đồng bộ khép kín: Sử dụng hầm gió kín hiện đại của IBST với hệ thống tạo nhám mô phỏng chính xác số mũ Hellman địa hình B Việt Nam, kết hợp cảm biến đa trục siêu nhạy JR3 đo đồng thời 6 thành phần nội lực đáy ở tần số quét cao $1000\text{ Hz}$, được kiểm chứng chéo với mô hình số 3D Navier-Stokes trung bình Reynolds (RANS) sử dụng mô hình rối SST $k-\omega$ và hàm UDF tự lập trình.

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

Phát hiện bất ngờ nhất là hiện tượng sụt giảm đột ngột của lực cản chính diện $F_x$ đi kèm với sự bùng nổ của lực ngang $F_y$ và mô men xoắn $M_z$ tại dải góc đón gió xiên $\theta = 45^\circ \div 52,5^\circ$. Cụ thể, trong khi $F_x$ giảm khoảng $20 \div 25%$ so với góc $\theta = 0^\circ$, hệ số lực ngang $C_{f,y}$ lại tăng vọt từ xấp xỉ $0$ lên $0,287$, kết hợp cùng cánh tay đòn lệch tâm $e_y = 0,25b$ tạo ra tổ hợp ứng suất uốn - nén - xoắn đồng thời tại chân cột đạt giá trị cực hạn nguy hiểm nhất, vượt qua cả trường hợp gió thổi vuông góc trực diện.

       BIẾN THIÊN HỆ SỐ LỰC VÀ ĐỘ LỆCH TÂM THEO GÓC GIÓ THỰC NGHIỆM

  Hệ số (Cf,x, ey/b)
          0°            30°       45°             90°    Góc gió (θ)

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

Hoàn toàn minh bạch và chi tiết. Luận án cung cấp đầy đủ các thông số vật lý và hình học: kích thước 6 mẫu pano M1–M6; quy cách bố trí Spire và cục tạo nhám trong ống thổi; thông số kỹ thuật cảm biến JR3 (bảng dải đo, độ phân giải kỹ thuật số, ma trận khử can nhiễu đa trục); cấu trúc lưới chia trong ANSYS Fluent kèm dải giá trị Skewness/Orthogonal Quality; mã code hàm UDF profile gió; và các bảng tính Excel mẫu cho các kích thước pano điển hình.

5. Định hướng chương trình nghiên cứu 10 năm tới được phác thảo như thế nào?

Lộ trình nghiên cứu dài hạn hướng tới: (1) Tích hợp công nghệ cảm biến không dây IoT đo đạc ứng suất và chuyển vị thời gian thực trên các pano thực tế ngoài hiện trường trong mùa mưa bão; (2) Nghiên cứu tối ưu hóa khí động học hình học mặt bảng (sử dụng bề mặt cong khí động học giảm lực cản); (3) Xây dựng tiêu chuẩn chuyên ngành riêng biệt cho công trình quảng cáo ngoài trời tại Việt Nam tương đương CECS 148 của Trung Quốc hay hướng dẫn ASCE của Mỹ.

Kết luận

  1. Luận án đã giải quyết căn bản và tường minh bài toán xác định tải trọng gió lên kết cấu pano tấm lớn đứng độc lập ngoài trời trong điều kiện thời tiết khắc nghiệt và bão nhiệt đới tại Việt Nam.
  2. Chứng minh một cách định lượng sự không an toàn của tiêu chuẩn cũ TCVN 2737:1995 (thiếu hụt từ $24,6%$ đến hơn $35%$ tải trọng và bỏ qua hoàn toàn mô men xoắn $M_z$), cung cấp luận cứ khoa học chuẩn xác để bảo vệ việc chuyển đổi sang TCVN 2737:2023.
  3. Xác lập chính thức bộ giá trị thiết kế khuyến nghị có độ tin cậy cao: hệ số lực cản chính diện $C_{f,x} = 1,6 \div 1,7$ và tỷ số độ lệch tâm ngang bắt buộc $e_y/b = 0,25$ cho các dạng pano chữ nhật hai mặt bọc kín có tỷ số $b/c = 1,0 \div 5,0$.
  4. Phát triển thành công quy trình mô phỏng số 3D CFD FVM tối ưu hóa với mô hình dòng rối SST $k-\omega$ và hàm UDF profile gió bão, đạt độ tương thích và sai số dưới $5%$ so với thực nghiệm vật lý hầm gió.
  5. Thực hiện thành công chuỗi thực nghiệm trên hệ thống cảm biến 6-DOF JR3 tần số cao $1000\text{ Hz}$ trong ống thổi khí động hiện đại của Viện IBST, đóng góp bộ dữ liệu thực nghiệm chuẩn quốc tế cho chuyên ngành Kỹ thuật Xây dựng công trình.
  6. Mở ra các hướng nghiên cứu chuyển tiếp mang tầm vóc quốc tế về dao động động lực học khí đàn hồi và giải pháp giảm chấn thụ động/chủ động cho công trình dạng tấm mỏng chịu tác động của gió bão.