Tổng quan nghiên cứu

Trong bối cảnh phát triển cơ sở hạ tầng hiện đại, các công trình kết cấu dạng tháp giàn thanh mảnh như tháp truyền hình, cột viễn thông và cột điện cao thế với chiều cao từ 40 mét đến trên 100 mét ngày càng xuất hiện phổ biến. Tuy nhiên, các công trình này đặc biệt nhạy cảm với tác động của tải trọng gió động, bao gồm hiện tượng tự kích khí đàn hồi và rung lắc ngẫu nhiên. Vấn đề nghiên cứu trọng tâm đặt ra là đánh giá định lượng sự suy giảm khả năng chịu lực khí động học khi kết cấu thép xuất hiện hư hại cục bộ như nứt ngầm hoặc mất liên kết bu lông.

Mục tiêu cụ thể của công trình là thiết lập hệ phương trình vi phân chuyển động không gian ba chiều chịu đồng thời hiệu ứng tự kích phân kỳ (Flutter) và gió rối ngẫu nhiên (Buffeting), từ đó xây dựng tiêu chí phát hiện tổn thương thông qua phân tích quỹ đạo dao động đỉnh tháp. Nghiên cứu được thực hiện tại Hà Nội với mô hình tháp giàn không gian cao 40 mét, kích thước đáy 2,5 mét x 2,5 mét. Ý nghĩa khoa học và thực tiễn của đề tài thể hiện qua việc chứng minh kết cấu có hư hỏng cột chịu lực sẽ bị sụt giảm tới 37,7% vận tốc gió tới hạn, khiến một công trình vốn được thiết kế chịu bão cấp 14 có thể bị sụp đổ hoàn toàn ở cấp gió 10 với vận tốc chỉ 28 mét/giây. Kết quả này cung cấp cơ sở quan trọng cho công tác cảnh báo sớm và bảo trì kết cấu thanh mảnh.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu ứng dụng lý thuyết dao động tự kích khí đàn hồi của Scanlan (1978) kết hợp cùng mô hình phân tích miền thời gian tích phân chập của Chen (2000) nhằm mô tả trọn vẹn lực khí động học tác động lên từng phân đoạn thanh giàn. Vận tốc gió dao động ngẫu nhiên được mô phỏng liên tục dựa trên mô hình hàm mật độ phổ năng lượng Van der Hoven (1957).

Khung lý thuyết vận hành dựa trên các khái niệm nền tảng bao gồm: hiện tượng mất ổn định Flutter với lực tự kích sinh ra do sự tương tác tương đối giữa bề mặt kết cấu và dòng khí; hiện tượng Buffeting do các xoáy khí rối gây rung lắc cưỡng bức; ma trận cản nhớt Rayleigh biểu diễn qua 2 tần số dao động riêng đầu tiên; và chỉ số ma trận chuyển đổi Boolean liên kết bậc tự do địa phương vào hệ tọa độ tổng thể. Mỗi nút trong mô hình không gian ba chiều được thiết lập với 6 bậc tự do độc lập, tương ứng 12 bậc tự do cho mỗi phần tử thanh dầm không gian.

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

Nghiên cứu sử dụng mô hình giàn thép không gian ba chiều cao 40 mét gồm 4 cột chính và hệ thanh giằng chéo phân bố theo các tầng. Các thông số vật liệu được thiết lập với mô đun đàn hồi Young E bằng 1,3 x 10^11 N/m2, khối lượng riêng của thép đạt 7849,0474 kg/m3, diện tích mặt cắt thanh cột là 0,02 m2 và thanh giằng là 0,01 m2. Mô hình chọn mẫu mô phỏng hư hỏng tại vị trí cao độ 20 mét với hai kịch bản độc lập: mất liên kết thanh cột số 134 và mất liên kết thanh giằng số 387, khảo sát mức độ sâu vết nứt từ 0% đến 100%.

Phương pháp phần tử hữu hạn dạng thanh dầm 3D được lựa chọn nhờ khả năng mô tả chính xác đồng thời độ cứng uốn, xoắn Saint-Venant và lực dọc trục. Để giải phương trình động lực học phi tuyến theo bước thời gian, tác giả áp dụng thuật toán tích phân số Newmark với hệ số gamma bằng 0,5 và beta bằng 0,25. Tác giả tự lập trình phần mềm chuyên dụng Wind Effects để tính toán số liệu và đối chuẩn trực tiếp với phần mềm thương mại SAP2000. Kết quả kiểm chứng 5 tần số dao động riêng đầu tiên cho thấy sai số hoàn toàn nằm trong mức cho phép, cụ thể: dạng 1 (1,1489 Hz) đạt sai số 0%, dạng 2 (1,1565 Hz) đạt sai số 0%, dạng 3 (1,8077 Hz) sai số 5%, dạng 4 (5,1260 Hz) sai số 4% và dạng 5 (6,4784 Hz) sai số 8%.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, hư hỏng kết cấu làm suy giảm đột ngột vận tốc gió tới hạn Flutter. Đối với tháp nguyên vẹn, vận tốc gió tới hạn đạt 45 mét/giây (tương đương bão cấp 14 theo thang Beaufort). Khi thanh giằng số 387 bị đứt gãy liên kết, vận tốc tới hạn giảm 17,7% xuống còn 37 mét/giây (tương ứng cấp 11). Nghiêm trọng hơn, khi thanh cột số 134 bị mất liên kết, vận tốc tới hạn giảm mạnh tới 37,7%, chỉ còn 28 mét/giây (tương ứng cấp 10).

Thứ hai, mức độ nhạy cảm của phản ứng động phụ thuộc vào độ sâu vết nứt. Khảo sát mô phỏng cho thấy các vết nứt bề mặt có độ sâu dưới 80% chỉ gây biến đổi nhỏ về tần số riêng, nhưng khi độ sâu vết nứt đạt từ 90% đến 100% thì biên độ dao động đỉnh tháp tăng vọt và làm thay đổi toàn diện quỹ đạo chuyển động.

Thứ ba, quỹ đạo chuyển động mặt phẳng Oxy của 4 nút đỉnh cột (nút 90, 180, 270, 360) dưới tác dụng của tải trọng gió Buffeting mang tính định danh hư hỏng rất cao. Ở trạng thái nguyên vẹn, quỹ đạo của 4 nút đỉnh hoàn toàn đồng dạng ở mọi góc gió từ 0 độ đến 45 độ. Khi thanh cột số 134 thuộc cột số 2 bị hỏng, nút đỉnh số 180 có biên độ chuyển vị lớn nhất; hai nút 90 và 270 liên kết trực tiếp qua giằng có biên độ trung bình; còn nút 360 không liên kết trực tiếp có biên độ dao động nhỏ nhất.

Thảo luận kết quả

Các kết quả phản ứng động lực học có thể được minh họa trực quan thông qua đồ thị hình chiếu quỹ đạo dao động đỉnh tháp trên hệ trục tọa độ hai chiều Oxy kết hợp bảng so sánh biến thiên biên độ chuyển vị theo từng góc hướng gió beta. Khi góc hướng gió thay đổi từ 0 độ đến 30 độ, quỹ đạo chuyển động của nút đỉnh chuyển đổi dần từ dạng đường dẹt kéo dài sang dạng elip mở rộng.

Nguyên nhân cốt lõi của sự bất đối xứng này xuất phát từ vai trò chịu lực chính của các cột đứng. Khi một thanh cột bị hư hại, độ cứng chống uốn tổng thể của nhánh cột đó suy giảm cục bộ, khiến năng lượng kích động từ dòng gió chuyển hóa mạnh thành biên độ dao động uốn và xoắn tại chính cột bị tổn thương. Phát hiện này hoàn toàn tương thích với các nguyên lý khí động lực học của Scanlan nhưng được phát triển mở rộng thành công cho kết cấu thanh mảnh không gian 3D. Điều này chứng minh rằng việc theo dõi hình chiếu quỹ đạo nút đỉnh là một công cụ chẩn đoán không phá hủy cực kỳ tiềm năng để xác định nhánh cột chịu lực đang bị hư hại.

Đề xuất và khuyến nghị

Thứ nhất, triển khai lắp đặt hệ thống cảm biến gia tốc và đo chuyển vị đa trục tại 4 góc đỉnh tháp đối với các công trình tháp truyền hình và trạm phát sóng viễn thông trọng điểm. Mục tiêu nhằm giám sát độ lệch quỹ đạo nút đỉnh với độ chính xác trên 95%, thời gian hoàn thành trong vòng 6 tháng do các đơn vị quản lý hạ tầng vận hành trực tiếp.

Thứ hai, cập nhật và điều chỉnh quy chuẩn tính toán tải trọng gió động cho công trình dạng thanh mảnh cao trên 30 mét. Cơ quan soạn thảo tiêu chuẩn cần bổ sung hệ số suy giảm an toàn khí động học tối thiểu 30% đến 40% để dự phòng kịch bản xuất hiện hư hỏng mỏi cục bộ trong quá trình khai thác, hoàn thiện văn bản trong vòng 12 tháng.

Thứ ba, thiết lập quy trình kiểm tra định kỳ 6 tháng một lần đối với hệ thống bu lông liên kết và mối hàn của các thanh cột chịu lực tại khoảng cao độ từ 15 mét đến 25 mét. Đây là khu vực chịu mô men uốn và chuyển vị động lớn nhất, giúp loại trừ 80% nguy cơ gãy đổ kết cấu trong mùa mưa bão.

Thứ tư, nâng cấp mã nguồn phần mềm Wind Effects lên phiên bản tự động hóa, tích hợp thuật toán học máy nhằm nhận dạng tự động vị trí vết nứt dựa trên dữ liệu quỹ đạo thời gian thực trong vòng dưới 15 giây xử lý, giao cho nhóm chuyên gia cơ kỹ thuật phát triển trong thời hạn 9 tháng.

Đối tượng nên tham khảo luận văn

Thứ nhất, kỹ sư thiết kế kết cấu công trình xây dựng và giao thông có thể ứng dụng trực tiếp thuật toán tích phân Newmark kết hợp lực khí đàn hồi để tính toán tối ưu hóa mặt cắt giàn thép, đảm bảo độ an toàn chống hiện tượng Flutter cho công trình cao trên 40 mét.

Thứ hai, các đơn vị quản lý, vận hành mạng lưới viễn thông và truyền tải điện lực sử dụng phương pháp phân tích quỹ đạo 4 nút đỉnh tháp làm tài liệu hướng dẫn kỹ thuật để đánh giá nhanh sức khỏe công trình sau các đợt bão lớn mà không cần dừng khai thác.

Thứ ba, nghiên cứu sinh và học viên cao học thuộc khối ngành Cơ học công trình, Kỹ thuật xây dựng công trình ngầm và mỏ, Cơ kỹ thuật sử dụng luận văn làm tài liệu tham khảo chuyên sâu về mô hình phần tử hữu hạn 3D Frame và kỹ thuật lập trình mô phỏng phổ gió ngẫu nhiên.

Thứ tư, các cơ quan kiểm định chất lượng công trình xây dựng cấp tỉnh và thành phố tham khảo số liệu để xây dựng định mức đánh giá rủi ro mỏi và độ suy thoái liên kết kim loại theo niên hạn sử dụng.

Câu hỏi thường gặp

Hiện tượng Flutter gây nguy hiểm như thế nào đối với công trình giàn cao tầng? Flutter là hiện tượng tự kích khí động học phân kỳ sinh ra do sự tương tác giữa kết cấu và dòng gió lớn. Khi đạt vận tốc tới hạn, biên độ dao động uốn và xoắn tăng nhanh không kiểm soát, dẫn đến phá hủy và sụp đổ công trình trong thời gian ngắn mà không có dấu hiệu dẻo báo trước.

Vết nứt cục bộ làm suy giảm khả năng chịu tải trọng gió của tháp như thế nào? Khi thanh cột ở cao độ 20 mét bị nứt gãy hoàn toàn, độ cứng của tháp giảm cục bộ khiến vận tốc gió tới hạn Flutter giảm từ 45 mét/giây xuống 28 mét/giây, tức giảm 37,7%. Một công trình an toàn ở gió cấp 14 có thể sụp đổ ngay khi gặp gió cấp 10.

Phương pháp quỹ đạo nút đỉnh có định vị được chính xác vị trí thanh bị gãy không? Phương pháp quỹ đạo nút đỉnh phân tích trên mặt phẳng Oxy xác định chính xác nhánh cột bị hư hại dựa trên nút có biên độ dao động lớn nhất trong 4 đỉnh cột. Tuy nhiên, phương pháp hiện tại chỉ khoanh vùng được cột bị hỏng chứ chưa chỉ ra chính xác số thứ tự thanh trên độ cao.

Phần mềm Wind Effects của tác giả có độ tin cậy ra sao khi so sánh với SAP2000? Wind Effects đạt độ tin cậy rất cao khi kiểm chứng cùng SAP2000 trên cùng mô hình giàn 3D. Sai số tần số dao động riêng ở 2 dạng đầu tiên bằng đúng 0%, dạng 3 lệch 5%, dạng 4 lệch 4% và dạng 5 lệch 8%, chứng minh thuật toán lập trình hoàn toàn chuẩn xác.

Sự khác biệt chính giữa lực gió Flutter và lực gió Buffeting là gì? Lực Flutter là lực tự kích sinh ra từ chính chuyển động dao động tương đối của kết cấu cắt qua dòng khí tĩnh, có thể gây mất ổn định phân kỳ. Trong khi đó, lực Buffeting là lực cưỡng bức ngẫu nhiên sinh ra do tính chất rối và biến đổi vận tốc của dòng khí tự nhiên.

Kết luận

  • Thiết lập thành công mô hình phần tử hữu hạn 3D Frame phân tích đồng thời tác động của hiện tượng tự kích Flutter và gió rối Buffeting trên kết cấu tháp giàn cao 40 mét.
  • Chứng minh định lượng mức độ sụt giảm nghiêm trọng của vận tốc gió tới hạn Flutter từ 45 mét/giây xuống còn 28 mét/giây (giảm 37,7%) khi xuất hiện hư hỏng ở thanh cột chịu lực.
  • Khám phá quy luật bất đối xứng trong quỹ đạo chuyển động 4 nút đỉnh tháp, chứng minh nút trên đỉnh cột bị hỏng luôn có biên độ dao động lớn nhất so với các nút còn lại.
  • Phát triển độc lập phần mềm tính toán chuyên dụng Wind Effects cho kết quả tương thích cao với phần mềm chuẩn quốc tế SAP2000 với sai số tần số dao động dưới 8%.
  • Đề xuất giải pháp giám sát sức khỏe kết cấu (SHM) chi phí thấp và hiệu quả cao thông qua việc quan trắc hình chiếu chuyển vị đỉnh tháp.

Về lộ trình nghiên cứu tiếp theo, trong giai đoạn 12 tháng tới, đề tài cần tiếp tục mở rộng mô hình hóa phi tuyến hình học và vật liệu để xác định chính xác tọa độ cao độ của vị trí hư hỏng. Các đơn vị nghiên cứu và doanh nghiệp xây dựng hãy chủ động ứng dụng mô hình phân tích khí động học tiên tiến này vào quy trình thiết kế và kiểm định để bảo vệ an toàn tối đa cho công trình tháp thanh mảnh trước hiểm họa thiên tai bão gió.