Sử dụng thiết bị điều chỉnh khối lượng để hạn chế chuyển vị ngang của kết cấu

Bài viết Sử dụng thiết bị điều chỉnh khối lượng để hạn chế chuyển vị ngang của kết cấu trình bày đề cập đến hai vấn đề chính: hạn chế chuyển vị ngang của kết cấu có chiều cao lớn nhờ sử dụng một thiết bị điều chỉnh khối lượng,... Mời các bạn cùng tham khảo.

KHẢO SÁT - THIẾT KẾ XÂY DỰNG SỬ DỤNG THIẾT BỊ ĐIỀU CHỈNH KHỐI LƯỢNG ĐỂ HẠN CHẾ CHUYỂN VỊ NGANG CỦA KẾT CẤU ThS PHÙNG NGỌC DŨNG Trường Đại học Kiến trúc Hà Nội Tóm tắt: Đối với kết cấu có chiều cao lớn chịu tải trọng động gió hay động đất, để đảm bảo chuyển vị ngang nằm giới hạn cho phép kết cấu phải có kích thước hình học định tạo độ cứng đủ lớn hạn chế chuyển vị ngang Việc tăng kích thước tiết diện dẫn đến nhiều hệ lụy như: chi phí tăng lên, biện pháp thi cơng phức tạp hơn,… Ngồi ra, tác dụng tải trọng động lên kết cấu cao tầng ống khói, tháp hay nhà nhiều tầng độ cản thân kết cấu (damping) có ảnh hưởng quan trọng tới ứng suất chuyển vị ngang kết cấu Bài báo đề cập đến hai vấn đề chính: hạn chế chuyển vị ngang kết cấu có chiều cao lớn nhờ sử dụng thiết bị điều chỉnh khối lượng (Tuned Mass Damper - TMD); tầm quan trọng độ cản kết cấu tới phản ứng chúng chịu tải trọng động Giới thiệu a Kết cấu nghiên cứu Kết cấu ống khói thép, thân hình trụ tròn Tiết diện ngang hình vành khun khơng đổi suốt chiều cao Ống khói xem conson ngàm móng (cốt 0.00) Về mặt kết cấu, ống khói mơ dầm Bernoulli, bỏ qua ảnh hưởng lực dọc, lực cắt Vật liệu thép xem làm việc giới hạn đàn hồi Các tính chất ống khói cho bảng sau: Bảng Các tính chất ống khói thép Chiều cao (m) Đường kính ngồi (m) Đường kính (m) Chiều dày (m) Mơ đun đàn hồi E (GPa) Khối lượng riêng  (kg/m ) Mơ men qn tính (m ) Mơ men chống uốn (m ) 38 3.6 3.585 0.015 200 7800 2.18491 1.21384 b Lực tác dụng Tải trọng tác dụng lên kết cấu gió Tải gió bao gồm hai thành phần tĩnh động Thành phần tĩnh xem lực tĩnh không thay đổi theo thời gian Việc thiết kế hay phân tích kết cấu tác dụng tải tĩnh đơn giản, nên không đề cập khuôn khổ báo Thành phần động mơ hình theo số phương pháp khác [4], [5] Khi tác dụng vào kết cấu mỏng cao ống khói, chế gió cuộn (vortex shedding) thường chiếm ưu [5] Nó mô tả hàm sine thời gian vận tốc [5] Vận tốc gió thay đổi khoảng (0:100) km/h, tức (0:27.8) m/s Hàm tải trọng gió thể sau: p(t )  p0 D sin(2 ns t ) (1) [4], [5] Trong đó: p0 – áp lực động lực - p0   air CdU (2) - air = 1.2kg/m ; Cd = 0.5; (air - tỷ trọng khơng khí Cd - hệ số); U – vận tốc trung bình gió; ns – tần số gió lốc - ns  SU (3) [5]; với S = 0.4; D – D kích thước vật chắn gió c Thiết bị điều chỉnh khối lượng - Tuned Mass Damper (TMD) TMD thiết bị gắn vào đỉnh ống khói nhằm hạn chế chuyển vị chịu tải trọng động Các tính chất TMD bao gồm: độ cứng (k), khối lượng (m) độ cản (damping-c) Trong thực tế giá trị xác định xác sau tính tốn phản ứng kết cấu với có mặt khơng có mặt TMD Tuy nhiên, khuôn khổ báo này, giá trị khối lượng tỷ số cản giả sử trước, sau dựa vào việc thay đổi giá trị độ cứng TMD, ta xác định giá trị độ cứng hợp lý TMD Tạp chí KHCN Xây dựng - số 3/2011 KHẢO SÁT - THIẾT KẾ XÂY DỰNG cho chuyển vị ngang TMD chuyển vị ngang toàn ống khói nhỏ giá trị cho phép (tương đương với việc tính tốn theo trạng thái giới hạn 2) Giá trị khối lượng m TMD thường nhỏ 300kg, hệ số cản damping () khoảng 5%-10% [1-3] d Các yêu cầu kiểm tra Dưới tác dụng tải trọng gió, phản ứng kết cấu (như chuyển vị, vận tốc, gia tốc, nội lực, ) thay đổi theo thời gian phụ thuộc vào vận tốc gió Các tiêu chí kiểm tra báo: - Xác định chuyển vị lớn ống khói với có mặt khơng có mặt TMD; - Kiểm tra kết cấu với tính chất khác TMD, với điều kiện ràng buộc là: chuyển vị ngang TMD phải nhỏ 20cm chuyển vị ngang kết cấu phải nhỏ 4cm; - Kiểm tra ứng suất cực đại đáy hai trường hợp: không kể đến độ cản kết cấu có kể đến độ cản kết cấu Ứng suất cực đại chân ống khói phải nhỏ 50MPa Các phương pháp giải bước giải Do tải trọng gió thay đổi theo thời gian vận tốc dạng hàm hình sin, phản ứng kết cấu động học Để xác định phản ứng điểm mặt cắt kết cấu, sử dụng số phương pháp sau: a Phương pháp Xác định trực tiếp chuyển vị từ phương trình dao động: ống khói xem dầm conson Bernoulli nên phương trình chuyển động phân tố vơ bé (khơng có cản) là:   2 ( x, t )      x, t    EI    p  x, t  (4) [1 - 3] t x  x  Trong đó:  - khối lượng riêng dầm;  - chuyển vị dầm; x – hướng trục dầm; EI – độ cứng uốn dầm; t – thời gian p – tải trọng phân bố phụ thuộc vào thời gian tọa độ Rõ ràng rằng, việc giải trực tiếp phương trình dao động công việc dễ dàng dầm có khối lượng phân bố nên số bậc tự dầm vô hạn b Phương pháp Xác định phản ứng kết cấu nhờ sử dụng phân tích phương thức (Modal Analysis): chuyển vị dầm chia thành hai thành phần độc lập  ( x, t )   ( x) (t ) (5); tích hàm dạng (x) (phụ thuộc vào tọa độ hình học vị trí xác định chuyển vị) tung độ giá trị chuyển vị (phụ thuộc vào thời gian -  (t)) Phương trình dao động tự dầm chia thành hai phương trình có hai biến số độc lập (khơng tính đến độ cản hai phương trình này):   ( x)   2 (t )   ( x )    2 (t )  (6)[1],[2] x EI t Giải hai phương trình (5) (6)Error! Reference source not found., ta xác định tính chất động học dầm (các chu kỳ, tần số, dạng dao động riêng) Nhờ vào tính trực giao dạng dao động, dầm với số bậc tự vô hạn chuyển thành dầm với bậc tự dạng dao động Độ cứng, khối lượng hay độ cản damping lực tác dụng chiếu lên dạng dao động riêng để xác định độ cứng chuẩn hóa, khối lượng chuẩn hóa, lực chuẩn hóa Việc giải phương trình tung độ để xác định giá trị  (t) trở nên dễ dàng lúc dạng dao động, hệ trở thành bậc tự Kết hợp tất dạng dao động (hoặc kết hợp vài dạng bản), ta xác định phản ứng tổng cộng kết cấu [1 - 3] c Phương pháp Phương pháp phần tử hữu hạn kết hợp với việc phân tích phương thức phương pháp số: ống khói chia thành nhiều phần tử rời rạc, phần tử xem phần tử dầm Bernoulli Mỗi phần tử có bậc tự ràng buộc chuyển vị nút Như vậy, từ dầm có số bậc tự vơ hạn chuyển thành dầm có số bậc tự hữu hạn Dựa vào phương pháp phân tích phương thức, ta tính tốn hệ có Tạp chí KHCN Xây dựng - số 3/2011 KHẢO SÁT - THIẾT KẾ XÂY DỰNG số bậc tự hữu hạn số dạng dao động kết cấu Đối với dạng, phương pháp số (phương pháp Newmark [1-3]) sử dụng để giải phương trình vi phân dao động hệ bậc tự do, từ xác định tung độ  k(t) dạng dao động thứ k Sau đó, nhờ việc tổ hợp tất dạng (hoặc tổ hợp số dạng đầu tiên) ta xác định phản ứng tổng cộng kết cấu d Phương pháp Phương pháp phần tử hữu hạn kết hợp với việc áp dụng thẳng phương pháp số để giải hệ nhiều bậc tự do: nhờ việc sửa đổi giải thuật Newmark ([1 - 3]) áp dụng cho hệ nhiều bậc tự (dạng ma trận), phản ứng hệ xác định trực tiếp từ phương trình dao động hệ nhiều bậc tự do: M v  C v  K v  P (6) Trong đó: M – ma trận khối lượng; C – ma trận cản; K – ma trận độ cứng; P - vector ngoại lực tác dụng nút; v – vector chuyển vị nút phần tử Chuyển vị ngang ống khói với có mặt hay khơng có mặt TMD xác định theo giải thuật Newmark thay đổi [1 - 3] e Trình tự tính tốn Hai phương pháp số sử dụng để tính tốn kết cấu Các bước để giải tổng kết bảng Bảng Các phương pháp sử dụng Phương Các pháp số bước Sơ đồ công việc cần thực bước Kết cấu khơng có TMD Chia ống khói thành hữu hạn đoạn, xác định đặc trưng phần tử: ma trận độ cứng, khối lượng, cản Nối tất ma trận phần tử để có ma trận tổng thể kết cấu: độ cứng K, khối lượng M Áp dụng điều kiện biên vào ma trận kết cấu ta có ma trận cuối kết cấu: ma trận độ cứng K1, ma trận khối lượng M1 Xác định trị riêng i, dạng dao động i kết cấu giải phương trình K1   M  Xác định tính chất động học số dạng liên quan tới chuyển vị ngang ống khói (trong trường hợp này, dạng có liên quan tới chuyển vị ngang) Xác định ma trận cản hệ với giả sử xác định theo ma trận cản Rayleight [1] xác định theo K1 M1: C1=K1+M1 Các hệ số   xác định với giả sử hệ số cản hai dạng dao động 1% Xác định ma trận lực tác dụng P(U,t) nút phần tử từ lực phân bố p(u,t) Chiếu ma trận độ cứng, khối lượng, cản kết cấu ma trận lực tác dụng nút lên bốn dạng dao động đầu tiên, ta có ma trận độ cứng chuẩn hóa (  chuẩn hóa (  T T K 1 ), khối lượng M 1 ), cản chuẩn hóa (  C1 ) ma trận lực chuẩn hóa Các ma trận T ma trận đường chéo Đối với dạng, sử dụng phương pháp số (phương pháp Newmark cho hệ bậc tự do) để giải phương trình chuyển động thời điểm t(i) (theo miền thời gian) Kết đạt tung độ (t) tất nút kết cấu Xác định phản ứng kết cấu dạng Xác định phản ứng tổng cộng kết cấu tất dạng dao động gây Kiểm tra điều kiện yêu cầu kết luận 1-4 Giống phương pháp Giải trực tiếp phương trình dao động nhờ sử dụng phương pháp gần Newmark cho hệ nhiều bậc tự (dạng ma trận) để xác định phản ứng tổng cộng kết cấu Kiểm tra điều kiện yêu cầu kết luận Tạp chí KHCN Xây dựng - số 3/2011 KHẢO SÁT - THIẾT KẾ XÂY DỰNG Kết cấu có TMD Phương pháp khơng áp dụng cho kết cấu với TMD ma trận cản kết cấu ma trận đường chéo trường hợp khơng có TMD Do vậy, phương trình tung độ kết cấu ( M t   C t   K  t   PU , t  , [M ], [K*] – ma trận đường chéo, [C ] – khơng * * * * * phải ma trận đường chéo) Giống phương pháp 3, kết cấu khơng có TMD 1-4 Xác định tính chất TMD: k – độ cứng, c – giá trị cản, m – khối lượng Xây dựng lại tính chất động học kết cấu nhờ việc nối K1, M1, C1 trước với tính chất động học k, m, c TMD để đạt ma trận độ cứng, khối lượng độ cản mở rộng tổng thể K2, M2 C2 Giải trực tiếp phương trình dao động dùng giải thuật Newmark áp dụng cho hệ nhiều bậc tự để xác định phản ứng tổng cộng kết cấu với TMD Kiểm tra điều kiện yêu cầu kết luận vấn đề Kết 3.1 Các tính chất kết cấu Tồn ống khói chia thành 16 phần tử hữu hạn rời rạc Như giả thiết, phần tử gồm có hai nút với bậc tự (2 chuyển vị thẳng chuyển vị xoay nút) phần tử phần tử dầm Bernoulli Như vậy, tổng cộng số bậc tự toàn kết cấu 51 3.1.1.Kết cấu không chưa gắn thiết bị TMD a Các ma trận độ cứng khối lượng Các ma trận độ cứng, khối lượng phần tử dầm xác định theo số phương pháp chuyển vị ảo Sau xác định ma trận phần tử, ma trận độ cứng khối lượng tổng thể kết cấu xác định cách nối ma trận phần tử Với 16 phần tử hữu hạn, ma trận khối lượng độ cứng kết cấu (K M) có kích thước 51x51 Do ống khói xem ngàm móng nên ba bậc tự phần tử thứ (sát chân ống khói) khơng Sau áp dụng điều kiện biên, ma trận độ cứng khối lượng kết cấu trở thành ma trận có kích thước 48x48 (K1 M1) b Các chu kỳ, tần số dạng dao động kết cấu Đặt chuyển vị ngang ống khói (x,t) theo:  ( x, t )   ( x) (t ) , để xác định chu kỳ, tần số dạng dao động ta cần xác định trị riêng hệ phương trình:  K   M ( x)  0 T=2/ (7) Trong đó: K – ma trận độ cứng kết cấu; K K1; M – ma trận khối lượng kết cấu, M M1;  - tần số dao động riêng kết cấu; T – chu kỳ dao động riêng kết cấu (x) – dạng dao động riêng ứng với tần số  chu kỳ T Dễ dàng nhận thấy rằng, để phương trình (8) thỏa mãn với điều kiện kết cấu có dao động tác dụng lực tức ((x) 0), hệ phải có nghiệm khơng tầm thường, đó, định     thức ma trận K   M phải không Tức là: det K   M  K   M  Thơng qua ngơn ngữ MATLAB với lệnh eig, ta xác định 48 trị riêng i Sắp xếp 1 < 2 < 3… T1=2/1> T2=2/2> T3=2/3…, ta xác định tần số chu kỳ dao động riêng kết cấu Tuy nhiên, quan tâm tới chuyển vị ngang ống khói nên ta cần loại bớt giá trị riêng liên quan tới chuyển vị xoay chuyển vị dọc trục kết cấu dạng đầu tiên: dạng dao động đầu tiên, liên quan tới chuyển vị ngang ống khói Giá trị chu kỳ tần số dạng dao động sau loại bỏ dạng dao động xoay dọc trục cho bảng Sau xác định tần số dạng dao động đầu tiên, thay giá trị tần số i Tạp chí KHCN Xây dựng - số 3/2011 KHẢO SÁT - THIẾT KẾ XÂY DỰNG vào phương trình (8) ta xác định vectơ riêng tương ứng với tần số Hình cho thấy dạng dao động ngang conson ống khói Các dạng dao động chuẩn hóa với giá trị chuyển vị đỉnh conson (bậc tự số 46) đơn vị Bảng Chu kỳ tần số dao động tự nhiên bốn dạng dao động ngang Dạng Tần số(Hz) Chu kỳ (s) 1/2/3/4 1 = 15.6604/2 =98.1428/ =274.8128/4 =538.5823 T1 = 0.4/ T2 =0.06/ T3 =0.02/ T4 =0.01 c Ma trận cản Hình Bốn dạng dao động ngang chuẩn hóa Hình Chuyển vị ngang đỉnh ống khói bốn dạng dao động U =50 km/h Như đề cập bước bảng 2, ma trận cản kết cấu xem ma trận cản Rayleight [1], [2] Ma trận cản kết cấu C1, xác định từ ma trận khối lượng độ cứng: C1=M1+K1 (8) Các hệ số   xác định nhờ giả sử hệ số cản  (daping ratio) hai dạng dao động 1% Chúng ta có tần số dao động riêng hai dạng dao động 1 2 Hệ phương trình dao động kết cấu chiếu lên tất dạng dao động là: M t   C  t   K  * * * t   P * U , t  (9) Ở dạng dao động thứ i, giá trị cản chuẩn hóa Ci*(generalised damping) phương trình dao động dạng thứ i: M i*i t   C i*i t   K i* i t   Pi * (11) xác định theo công thức: C i*   Ti C1 i  (1210) Thay (9) vào (12) ta có:   C i*   Ti C1 i   M i*   K i* (13) * * Ci  i  * Từ (11), ta chia hai vế cho Mi , đặt   K i ;  i  (4) ta có phương trình   i * * 2 i M i 2 i Mi xác định tung độ dao động theo thời gian dạng thứ i là: P* (1511) i t   2 i  i i t    i2 i t   i * Mi Từ (14), với giá trị i hai dạng 1%, ta dễ dàng xác định hai hệ số  : Tạp chí KHCN Xây dựng - số 3/2011 KHẢO SÁT - THIẾT KẾ XÂY DỰNG   2    2  1     0.01    0.2701           2     0.01   0.0001757   2 (16) Thay (16) vào (9) ta xác định ma trận cản kết cấu Giá trị i hệ số cản damping dạng dao động thứ i, thường từ 1-2% cho kết cấu thép 3.1.2 Kết cấu với có mặt TMD Khi ống khói gắn thêm thiết bị TMD đỉnh, tính chất động học sẵn có kết cấu bị thay đổi Hệ kết cấu từ 51 bậc tự trở thành 52 bậc tự do, sau gắn điều kiện biên, hệ cũ có 48 bậc tự trở thành 49 bậc tự Bậc tự thứ 49 chuyển vị ngang TMD Các tính chất học hệ với có mặt TMD xác định với quan niệm hệ cũ (chưa gắn TMD) hệ bậc tự với độ cứng K1, khối lượng M1, độ cản C1 nối với TMD có độ cứng k, khối lượng m, độ cản c Với quan niệm trên, phương trình dao động hệ viết thành:      M   u1  C1  c c   u1   K  k    m      c c      k   u    2 u     k   u1   P1(t )     k      u  2 (17) Trong đó: u1 – chuyển vị quy ước hệ cũ khơng có TMD; u2 – chuyển vị TMD; P1(t) – lực quy ước tác dụng lên hệ cũ Từ phương trình (17), ta dễ dàng xác định tính chất học cho hệ mới: ma trận độ cứng: C1  c c   M  ; ma trận cản  K1  k k  ; ma trận khối lượng: C2   K2   M2    c  k   k  m  c Cần nhấn mạnh rằng, TMD gắn vào đỉnh ống khói, khơng làm thay đổi tính chất học tồn kết cấu mà làm thay đổi đặc tính học phần tử (phần tử 16) gắn với TMD Do vậy, giá trị độ cứng K1+k độ cản C1+c phương trình (5) hệ tính tốn tổng độ cứng k, độ cản c TMD với độ cứng K1 độ cản C1 vào bậc tự số 46 (bậc tự tương ứng với chuyển vị ngang đỉnh) thuộc phần tử thứ 16 3.2 Xác định lực tác dụng lên hệ a Đối với phương pháp Theo phương pháp này, lực phân bố chuyển thành lực tập trung nút phần tử, lực tập trung sau chiếu lên bốn dạng dao động gọi lực chuẩn hóa dạng dao dộng Giá trị lực chuẩn hóa tác dụng lên kết cấu dạng dao động xác định theo bước sau: nodej - Chuyển tải phân bố p(t) thành lực tập trung Pc (t) = p(t)*l nút phần tử (trong đó: l – chiều dài phần tử) Lực tập trung đỉnh ống khói ½ giá trị nút khác; - Chiếu lực tập trung lên dạng dao động riêng để xác định lực tác dụng dạng Ví dụ dạng thứ i: NDOFs Pi  t     inodej ( x) Pcnodej (t ) (1812) j 1 (NDOFs – tổng số bậc tự hệ; nod - nút) b Đối với phương pháp Tạp chí KHCN Xây dựng - số 3/2011 KHẢO SÁT - THIẾT KẾ XÂY DỰNG Theo phương pháp này, lực phân bố gió chuyển thành lực tập trung nút nodej phần tử Các giá trị lực tập trung tính theo Pc (t) = p(t)*l nút (l – chiều dài phần tử) Lực tập trung đỉnh ½ so với nút khác Cần lưu ý rằng, lực phân bố gió chuyển thành lực tập trung tương ứng với bậc tự có chuyển vị ngang hệ Các lực tập trung tác dụng theo chuyển vị khác không 3.3 Các tham số phương pháp Newmark Các tham số phương pháp bao gồm hệ số alpha, delta, bước thời gian DT, số lượng bước thời gian Nstep [1 – 3] Chúng lựa chọn cho phản ứng kết cấu ổn định đảm bảo độ ổn định phương pháp số Newmark Do đó, tổng thời gian tác dụng lực phải lớn T1/2 tỷ số DT/T1 phải nhỏ 0.551 [1], [2] Giá trị tham số chọn là: alpha=0.25, delta=0.5; DT=0.02s Nstep =2500 3.4 Phản ứng kết cấu khơng có TMD 3.4.1.Chuyển vị ngang đỉnh ống khói a Đối với phương pháp Để xác định phản ứng tổng cộng kết cấu, ta chiếu ma trận độ cứng, khối lượng độ cản lên dạng dao động để xác định ma trận độ cứng, khối lượng ma trận cản chuẩn hóa Sau thực phép chiếu lên dạng dao động, phương trình dao động trở thành phương trình số (9) Các giá trị * * * [M ], [K ] [C ] ma trận đường chéo Các giá trị đường chéo ma trận chuẩn hóa xác định cho dạng thứ i sau: M i*   i 1xn M 1 i nx1 ; C i*   i 1 xn C1 i nx1 ; K i*   i 1xn K1 i nx1 T T T Rõ ràng, hệ phương trình số (3) hệ gồm phương trình độc lập, khơng phụ thuộc vào Ta xác định vectơ tung độ {(t)}i dạng thứ i thuật giải số Newmark cho hệ bậc tự Ngôn ngữ MATLAB lựa chọn để lập giải thuật Tuy nhiên cần nhấn mạnh rằng, tải trọng gió thay đổi theo thời gian vận tốc gió U, để xác định chuyển vị ngang kết cấu, giá trị vận tốc giữ nguyên, sau cho vận tốc gió thay đổi (U trở thành vectơ) để xác định giá trị vận tốc gió nguy hiểm Sau xác định giá trị {(t)}i, chuyển vị kết cấu xác định theo  ( x, t )   ( x) (t ) tổ hợp tất dạng Hình thể giá trị chuyển vị đỉnh ống khói thay đổi theo thời gian bốn dạng dao động đầu tương ứng với vận tốc gió U = 50km/h=13.89m/s Từ hình ta thấy rằng, đóng góp dạng dao động vào chuyển vị ngang hệ chiếm ưu Ba dạng dao động lại ảnh hưởng tới phản ứng hệ Hình cho thấy chuyển vị ngang lớn đỉnh ống khói với giá trị vận tốc gió U thay đổi (từ 0km/h đến 100 km/h) Hình thể giá trị chuyển vị đỉnh ống khói bốn dạng dao động tương ứng với vận tốc U = 79.92km/h=22.2m/s Như vậy, chuyển vị lớn đỉnh ống khói 0.2515m tương ứng với giá trị vận tốc gió U = 22.2m/s (79.92 km/h) Tại U=22.2m/s, tần số lực gió tác dụng   2 SU  15.49856 , D đó, tần số dạng dao động 1 = 15.6604 Có thể nhận giá trị vận tốc gió U=22.2m/s, tần số lực tác dụng gần tần số dao động riêng nguy hiểm hệ, tượng cộng hưởng xảy Tại giá trị vận tốc này, kết cấu nằm dải cộng hưởng Tạp chí KHCN Xây dựng - số 3/2011 KHẢO SÁT - THIẾT KẾ XÂY DỰNG Hình Các chuyển vị ngang lớn đỉnh với giá trị vận tốc U khác Hình Các chuyển vị ngang đỉnh bốn dạng dao động với U = 79.92 km/h b Đối với phương pháp Hình Giá trị chuyển vị ống khói đỉnh theo thời gian với U =50 km/h – Phương pháp Hình Các chuyển vị ngang đỉnh ống khói với vận tốc gió U =79.92 km/h – Phương pháp Theo phương pháp này, phản ứng kết cấu xác định trực tiếp từ phương trình dao động dạng ma trận nhờ áp dụng trực tiếp giải thuật Newmark cho hệ nhiều bậc tự Phương trình dao động kết cấu là:     M 1 u  C1 u   K1 u  P (t ) (1913) Trong đó: Các giá trị M1, C1, K1 P(t) xác định Hình thể chuyển vị đỉnh ống khói theo thời gian với giá trị vận tốc gió U = 50 km/h = 13.88889m/s Các giá trị chuyển vị đỉnh ống khói với U = 50 km/2 phương pháp gần giống phương pháp Hình cho thấy chuyển vị đỉnh ống khói theo thời gian với vận tốc gió U = 79.92 km/h=22.2m/s Tương tự phương pháp 3, chuyển vị lớn đỉnh ống khói 0.2515m tương ứng với giá trị vận tốc gió U = 22.2m/s (79.92 km/h) Tại U=22.2m/s, tần số lực gió tác dụng   2 SU  15.49856 , đó, tần số dạng dao động 1 = 15.6604 Có thể dễ dàng D nhận rằng, giá trị vận tốc gió U=22.2m/s, tần số lực tác dụng gần tần số dao động riêng nguy hiểm hệ, tượng cộng hưởng xảy 3.5 Phản ứng kết cấu với TMD Bảng Chuyển vị ngang TMD đỉnh ống khói với giá trị độ cứng, khối lượng khác TMD Tạp chí KHCN Xây dựng - số 3/2011 KHẢO SÁT - THIẾT KẾ XÂY DỰNG Khối lượng m – kg 300 400 280 800 Các tính chất TMD Hệ số 1/2 Độ cứng k =(k/m) tần cản  – EI số TMD (%) 14412 6.931 28825 9.802 43237 12.005 57649 13.86 72061 15.499 79267 16.255 108090 18.98 96082 15.499 67257 15.499 192160 15.499 Vận tốc gió cộng hưởng (m/s) Chuyển vị lớn đỉnh ống khói khơng có TMD (m) 22.2 0.2515 22.2 22.2 22.2 0.2515 0.2515 0.2515 Chuyển vị lớn đỉnh ống khói có TMD (m) 0.1865 0.1278 0.0764 0.0384 0.0240 0.0241 0.0551 0.0206 0.0247 0.0174 Chuyển vị lớn TMD (m) 0.0467 0.0829 0.1070 0.1187 0.1490 0.1471 0.1550 0.1156 0.1572 0.0644 Khi kết cấu có thêm thiết bị TMD, ma trận độ cứng, khối lượng độ cản bị thay đổi Nếu dùng phương pháp để xác định phản ứng kết cấu, chuyển ma trận hệ thành ma trận phương trình phương thức, ma trận khối lượng độ cứng chuẩn hóa ma trận đường chéo Tuy nhiên, ma trận cản hệ không ma trận đường chéo có mặt hệ số cản TMD Phương trình dao động hệ phân tích phương thức là:  M   C    K    P (t )  ( - tung độ) (2014)   Trong đó: ma trận độ cứng chuẩn hóa: K *  T49 x 49 K 249 x 49 49 x 49 - ma trận đường chéo; ma trận khối T lượng chuẩn hóa: M *  49 x 49 M 249 x 49 49 x 49 - ma trận đường chéo; ma trận độ cản chuẩn hóa: T C *  49 x 49 C 249 x 49 49 x 49 ma trận đường chéo; [P(t)] - ma trận lực tập trung nút     Dễ dàng nhận thấy, phương trình (20) hệ phương trình lệ thuộc lẫn sử dụng phương pháp Newmark cho hệ bậc tự để giải trường hợp Do vậy, phương pháp không áp dụng cho trường hợp hệ có TMD Vậy phương pháp lựa chọn để giải hệ phương trình phụ thuộc Hệ phương trình dao động hệ có TMD là:  M 2          u  C u  K u   P(t ) (2115) Các giá trị M2, C2, K2 P(t) xác định Phần tính tốn hệ khơng có TMD cho thấy gia tốc gió 22.2 m/s xảy tượng cộng hưởng hệ dao động đỉnh ống khói 0.2515m Giá trị vận tốc gió U=22.2m/s, sử dụng để đánh giá ảnh hưởng TMD lên phản ứng hệ Như đề cập phần 1.3, để xác định tính chất TMD, ta cần phải thử với nhiều giá trị khác khối lượng, độ cứng độ cản TMD Tuy nhiên, để đơn giản hóa q trình tính tốn, giá trị khối lượng độ cản giả thiết khơng đổi, sau độ cứng TMD xác định cho giá trị chuyển vị ngang TMD đỉnh ống khói nằm giới hạn cho phép Bảng thể giá trị chuyển vị ngang TMD (bậc tự thứ 49) chuyển vị ngang đỉnh ống khói với giá trị khác độ cứng khối lượng TMD Rõ ràng rằng, với m=300kg, tần số TMD tăng dần nhỏ tần số hệ tải trọng gió chuyển vị ngang đỉnh ống khói giảm dần chuyển vị ngang TMD tăng dần lên Ngược lại, tần số TMD lớn tần số hệ tải trọng gió chuyển vị ngang hệ tăng dần lên chuyển vị ngang TMD tăng dần lên Khi tần số TMD SU đạt 15.499, với tần số hệ tải trọng gió   2  15.49856 , chuyển vị ngang đỉnh ống khói D nhỏ (0.024m=2.4cm