zunia.vn

Tuyển sinh 2024 dành cho Gen-Z

zunia.vn

» Luận Văn - Báo Cáo

» Báo cáo khoa học

NGHIÊN CỨU PHẢN ỨNG ĐỘNG LỰC CẦU DÂY VĂNG DƯỚI TÁC DỤNG CỦA TẢI TRỌNG DI ĐỘNG

Chia sẻ: Phạm Đức Linh | Ngày: | Loại File: PDF | Số trang:6

Báo xấu

223
lượt xem 86
download

Download Vui lòng tải xuống để xem tài liệu đầy đủ

Trong phần nghiên cứu này, thể hiện cái nhìn tổng quan về động lực học và đưa ra một phương pháp phân tích đơn giản để đánh giá phản ứng động của CDV dưới tác dụng của tải trọng di động. Cây cầu được lý tưởng hóa theo mô hình Bernoulli-Euler trên gối đàn hồi với độ cứng gối thay đổi.

Chủ đề:

Bình luận(0) Đăng nhập để gửi bình luận!

Lưu

Nội dung Text: NGHIÊN CỨU PHẢN ỨNG ĐỘNG LỰC CẦU DÂY VĂNG DƯỚI TÁC DỤNG CỦA TẢI TRỌNG DI ĐỘNG

Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008 NGHIÊN CỨU PHẢN ỨNG ĐỘNG LỰC CẦU DÂY VĂNG DƯỚI TÁC DỤNG CỦA TẢI TRỌNG DI ĐỘNG A STUDY RESPONSE OF CABLE-STAYED BRIDGES TO MOVING VEHICLES SVTH : PHAN HOÀNG NAM Lớp 03X3C – Trường Đại học Bách khoa, ĐHĐN GVHD: GVC. Th.S LÊ VĂN LẠC Khoa XDCĐ – Trường Đại học Bách khoa, ĐHĐN TÓM TẮT Trong phần nghiên cứu này, thể hiện cái nhìn tổng quan về động lực học và đưa ra một phương pháp phân tích đơn giản để đánh giá phản ứng động của CDV dưới tác dụng của tải trọng di động. Cây cầu được lý tưởng hóa theo mô hình Bern oulli-Euler trên gối đàn hồi với độ cứng gối thay đổi. Để giải phương trình dao động của cầu, tác giả sử dụng phương pháp sai phân hữu hạn và kỹ thuật xếp chồng các mode dao động. ABSTRACT This part of the thesis presents a state-of-the-art review and a simplified analysis method for evaluating the dynamic response of cable-stayed bridges to moving vehicles. The bridge is idealized as a Bernoulli-Euler beam on elastic supports with varying support stiffness. To solve the equation of motion of the bridge, the finite difference method and the mode superposition technique are used. 1. Mở đầu Do tính thẩm mỹ, ứng dụng hiệu quả các vật liệu xây dựng và những ưu điểm đáng kể khác mà cầu dây văng (CDV) đã được xây dựng phổ biến trong những thập niên gần đây. Loại cầu này đang có xu hướng tăng khả năng vượt nhịp, sử dụng các loại vật liệu có cường độ cao để giảm trọng lượng đến mức thấp nhất. Khi đó kết cấu trở nên thanh mảnh, gọn nhẹ hơn và rất nhạy cảm với các tác động của tải trọng. Kể từ vụ sụp cầu thảm khốc của cây cầu Tacoma Narrows đầu tiên năm 1940, người ta đã chú ý nhiều hơn đến đặc tính động lực học của cầu treo. Trong suốt 68 năm qua, nhiều nghiên cứu thực nghiệm và lý thuyết đã được tiến hành nhằm hiểu rõ hơn về các yếu tố khác nhau ảnh hưởng đến phản ứng của loại công trình này đối với gió, động đất và xe cộ lưu thông. Để nghiên cứu những ảnh hưởng về động lực học do xe cộ lưu thông trên cầu, các kĩ sư trên thế giới dựa vào hệ số khuếch đại động (DAF) được xác định trong các tiêu chuẩn thiết kế. Những hệ số này thường là hàm của tần số dao động riêng của kết cấu nhịp hoặc chiều dài phạm vi ảnh hưởng của tải trọng. Đây là phương pháp thường được dùng cho mục đích thiết kế có thể cho ra các thiết kế vững chắc và đắt tiền đối với một số cây cầu nhưng lại đánh giá thấp các ảnh hưởng động lực học đối với các cây cầu khác. Hơn nữa, khi so w (t) v(t) sánh các tiêu chuẩn thiết kế ở mỗi quốc gia, các hệ số DAF dao 2 động trong một phạm vi lớn. m2 Vì vậy, để kiểm tra khả năng chịu đựng của các cây cầu đang tồn tại đối với sự lưu thông nặng nề hơn và để thiết kế các cây cầu mới hợp lý hơn thì cần phải sử dụng các kỹ thuật phân w1(t) ks cs tích cải tiến có xét đến tất cả các thông số quan trọng ảnh hưởng đến phản ứng động lực học của cầu như: sự tương tác cầu - xe m1 cộ, độ gồ ghề của bề mặt đường, hệ số cản của cầu, … Hçnh 2-1: Mä hçnh hoïa xe 2. Nội dung 2.1. Mô hình hóa xe 186
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008 Xe nặng với một số bộ phận cơ bản như: máy kéo, rơmooc, hệ thống kéo, … có thể được mô hình hóa như một chất điểm có khối lượng đặt trên hệ lò xo và hệ giảm chấn (hình 2-1). Gọi F(t) là lực liên hệ giữa xe và cầu, được xác định là dương khi hướng xuống dưới cầu, phương trình vi phân chuyển động được thiết lập như sau: d 2w dw dw (m1  m2 ) g  m1 2 1  ks ( w2  w1 )  cs ( 2  1 )  F (t )  0 (2.1) dt dt dt d 2 w2 dw dw m2  ks ( w2  w1 )  cs ( 2  1 )  0 (2.2) 2 dt dt dt Với: w1(t), w2(t) - Chuyển vị theo phương đứng của vật m1, m2 theo thời gian. ks - Độ cứng của lò xo ; cs – Hệ số cản nhớt ; g - Gia tốc trọng trường. Từ (2.1) và (2.2) ta có thể xác định được lực liên hệ F(t) giữa xe và cầu như sau: d 2 w1 d 2 w2 F (t )  (m1  m2 ) g  m1  m2 (2.3) dt 2 dt 2 Giả thiết rằng xe luôn luôn tiếp xúc với cầu khi đó lực F(t)>0, và biến dạng giữa hệ vật và bề mặt của cầu không đáng kể, ta có: y y r w1 (t )  y( x(t ), t )  r ( x(t )) ; w1 (t )  v   v  (2.4a-b) x t x 2 y 2 y y  2 y  2 r 2 r 1 (t )  2 v 2  2v  a  2  2 v  a w (2.4c) x x.t x t x x   Với: w1 (t ) và w1 (t ) - Vận tốc và gia tốc theo phương đứng của khối vật. v và a - Vận tốc và gia tốc theo phương dọc cầu. y(x,t), r(x) - Hàm chuyển vị và hàm thể hiện quy luật bề mặt của cầu. 2.2. Mô hình hóa kết cấu cầu y(x,t) v(t) m2 Li, Ai, Ei Hi m1 x Ig, Eg, mg xv(t) s Ls=99.5m Ls=99.5m L=204m v(t) m2 y(x,t) m1 kx Ig, Eg, mg x Hình 2-2: Mô hình hóa kết cấu CDV 3 nhịp 99,5m+204m+99,5m 2.2.1. Các giả thiết chính - Hệ thống nhiều cáp với khoảng cách các điểm neo cáp nhỏ so với chiều dài của cầu. - Bỏ qua khối lượng cáp, hệ số cản của cầu. - Các cáp được lý tưởng hóa như các lò xo liên tục phân bố theo phương dọc cầu. - Lực căng của cáp dưới tác dụng của tĩnh tải được điều chỉnh sao cho chuyển vị tại nút bằng 0. - Chỉ xét mômen uốn trong mặt phẳng, bỏ qua mômen xoắn do tải lệch tâm trên mặt cầu. - Khi xe qua cầu, chuyển vị và vận tốc di chuyển theo phương đứng của xe được xem bằng 0. 2.2.2. Thiết lập phương trình vi phân chuyển động Theo [1], phương trình vi phân chuyển động cho sự dao động của cầu tại các mặt cắt của dầm chủ (được lý tưởng hóa theo mô hình Euler-Bernoulli trên các gối đàn hồi): 187
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008  4 y ( x, t )  2 y ( x, t )   .( x  xv ).F (t )  k ( x). y ( x, t )  mg . (2.5) Eg .I g . x 4 x 2 Với:  - Hàm Delta-Dirac ; Eg - Mô đun đàn hồi của dầm chủ Ig - Mômen quán tính của mặt cắt ngang dầm chủ mg - Khối lượng trên 1 đơn vị chiều dài dầm ; k(x) - Độ cứng của lò xo Các điều kiện biên: (2.6a-d, 2.7a-b)  2 y (0, t )  2 y ( L, t ) y ( x, 0) 0  0 ; y( L, t )  0 ;  0 ; y(x,0) = 0 ; y(0, t )  0 ; t x x 2 2 2.2.3. Độ cứng của lò xo Theo [7], độ cứng của lò xo lý tưởng hóa từ dây cáp thứ i, được xác định: E . A .sin 2 i ki  i i (2.8) Li Biểu thị ứng suất cho phép của cáp là  a , tĩnh tải và hoạt tải trên 1 đơn vị chiều dài là qg và qq diện tích mặt cắt ngang của cáp i được tính như sau: (qg  qq ).s Ai  (2.9)  a .sin i Tính phi tuyến hình học của cáp do sự thay đổi của độ chùng và hình dạng dưới các lực căng khác nhau, được xét đến gần đúng bằng cách giới thiệu môđun đàn hồi tương đương [5].  L Ec E ( x)  Ei  0  x   (2.10)  c .x  2 22 1 .Ec 12 03 Với: Ec, c - Môđun đàn hồi, khối lượng riêng của vật liêu chế tạo cáp. 0 - Ứng suất kéo ban đầu trong cáp. Lực trong cáp gây ra bởi hoạt tải xe nhỏ so với lực do tĩnh tải nên 0 được thay bằng g ứng suất trong cáp do tĩnh tải gây ra. qg g  a (2.11) (qg  qq ) Thay thế (2.9) vào (2.8) phương trình sau đây được thiết lập cho độ cứng của lò xo trên một đơn vị chiều dài của cáp. E ( x).(qg  qq )  L 1 k ( x)  0  x   . (2.12)  a .H x2 2 1 ( ) H 2.3. Phân tích phản ứng động 2.3.1. Dao động tự do Phương trình vi phân cho dao động tự do của mô hình cầu:  4 y( x, t )  2 y( x, t )  k ( x).y( x, t )  m g . 0 E g .I g . (3.1) x 4 x 2 Hàm chuyển vị y(x,t) có thể được biểu diễn bởi 2 hàm, hàm tọa độ không gian z(x) và hàm tọa độ thời gian (t). Khi đó dao động của cầu ứng với mode dao động thứ i và chuyển vị cầu tại các vị trí khác nhau ứng với thời gian t được diễn đạt như sau: yi ( x, t )  zi ( x)i (t )  zi ( x)(ai cos it  bi sin it ) (3.2) Với: ωi là tần số góc dao động thứ i của dao động tự do của cầu. 188
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008 Thay thế phương trình (3.2) vào phương trình (3.1) ta được: d 4 y ( x, t )  k ( x).z i ( x)  m g .i2 z i ( x)  0 (3.3) E g .I g . 4 dx Phân chia mô hình cầu thành n đoạn giống nhau với chiều dài mỗi đoạn là h. Lấy (n-1) phần chưa biết thay thế vào phương trình (3.3), theo phương pháp sai phân hữu hạn thì phương trình trên được thiết lập như sau: AZ = λZ (3.4)  c1 4 1 0  4 c 4 1    2  1 4 c3 4 1    A  (3.5a-f)   1 4 c3 4 1    1 4 c2 4   0 1 4 c1   n 1n 1   1  z1,1  z1,n 1  0     ; Z   z1 , z2 , z3 ,..., zn 1          zn 1,1  zn1,n1  n 1  ( n 1n 1) 0     ( n 1n 1) h4 k ( x 2,3 ) 4 k ( x1 ) 4 i  mg . i2 ; c1  5  h ; c 2,3  6  h Eg I g Eg I g Eg I g 2.3.2. Phân tích động Sử dụng kỹ thuật xếp chồng các mode dao động, ta tìm nghiệm phương trình (2.5) dưới dạng: y ( x, t )   i 1 zi ( x)i (t )  (3.6) Thay (3.6) vào (2.5) ta được: ..   ( Eg I g ziIV i  k ( x) zii  mg zi  i )   ( x  xv ).F (3.7) i 1 F .. Từ đó suy ra:  i  i2i   zi ( xv ) (3.8) mg Nghiệm của phương trình (3.8): . Fz ( x ) ( Fzi ( xv )  0ii mg )  2 i   i2 v  cosi t  0i sin i t (3.9) i mg i mg i 2 . Với 0i và  0i là những giá trị ban đầu cho đừng đoạn, y0 và y0 là chuyển vị theo  phương đứng và vectơ vận tốc của cầu ở đoạn trước đó, ta có: . . 0i  ziT y0 ;  0i  ziT y0 (3.10a-b)  Với y0 và y0 là chuyển vị theo phương đứng và vectơ vận tốc của cầu ở đoạn trước đó. Thay vào phương trình (3.6) ta có được vectơ chuyển vị theo phương đứng của cầu tại mỗi đoạn:   . Fzi ( xv ) ( Fzi ( xv )  0ii mg )  0i 2 s y   zi   2 sin i t  cosi t   (3.11)  i mg  i2 mg i i 1   189
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008 Với s là số mode dao động, s  n-1. Vectơ vận tốc và vectơ gia tốc của cầu được xác định như sau:  ( Fzi ( xv )  0ii2 mg )  s .    zi   sin i t   0i cos i t  y (3.12)   i mg   i 1  ( Fzi ( xv )  0ii2 mg )  s .    zi   cosi t   0i i sin i t  y (3.13)   mg   i 1 Kỹ thuật xếp chồng các mode dao động cũng được sử dụng để xác định vectơ mômen uốn trong dầm tại mỗi đoạn. s m   Eg I g y ''   Eg I g  zi ''i (3.14) i 1 Sử dụng ngôn ngữ lập trình MATLAB để phân tích phản ứng động, sự tương tác động giữa xe và cầu, các bước lặp, kiểm tra, tính toán được thể hiện ở thuật toán bên dưới. Điều kiện hội tụ, trong kết quả bài toán này chọn Tolerance = 5.10-6:   max y j  y j 1  Tolerance  5.106 (3.15)  max y j Với yj – Vectơ chuyển vị của cầu tại đoạn đang xét. 2.3.3. Phân tích tĩnh Véc tơ chuyển vị tĩnh được xác định như sau: 0     (m1  m2 ) gh 1  3 Pj  1   j yj   st A Pj ; (3.16a-b)   Eg I g  0   -1 Với: A là nghịch đảo của ma trận A được xác định từ biểu thức (3.5a). j - Điểm chịu tác động của tải trọng xe. y sjt - Vectơ chuyển vị của cầu trong trường hợp tải trọng xe tác dụng vào điểm thứ j. Mômen uốn trong dầm ở từng đoạn: Eg I g mistj   2 ( yist 1, j  2 yistj  yist1, j ) (3.17)   , , h 2.4. Ứng dụng phân tích tĩnh và động CDV Mỹ Thanh – Tỉnh Sóc Trăng 2.4.1. Số liệu đầu vào Xe: m1 =3000 (kg), m2 =31700 (kg), cs = 8,6.104 (Ns/m), ks = 9,12.106 (N/m). Cầu: Ec = 1,97.1011 (N/m2), Eg = 0,38.1011 (N/m2), Ig = 1,48 (m4), a = 750.106 N/m2, qg = 17.104 (N/m), qq = 6.104 (N/m), c = 2,5.104 (N/m3), mg = 2,71.105 (N/m). Phân chia mô hình cầu thành 403 đoạn, chiều dài mỗi đoạn h = 1m. Xét 10 mode dao động đầu tiên. 2.4.2. Kết quả Tần số dao động tự do ứng với các mode dao động: Mode dao động 1 2 3 4 5 6 7 8 9 10 F (Hz) 0,36 0,47 0,60 0,76 0,99 1,29 1,65 2,08 2,56 3,11 190
Tuyển tập Báo cáo “Hội nghị Sinh viên Nghiên cứu Khoa học” lần thứ 6 Đại học Đà Nẵng - 2008 Biểu đồ chuyển vị tại vị trí giữa 4.0 nhịp giữa DAFd = 1,062 (v = 30m/s ; Dynamic x = 201,5m). Vehicle displacement (mm) 0 Biểu đồ mômen uốn tại vị trí Static giữa nhịp giữa DAFm = 1,094 (v = -4.0 30m/s ; x = 201,5m). -8.0 3. Kết luận và kiến nghị Đối với các nghiên cứu sơ bộ, sử -12.0 dụng các mô hình CDV đơn giản để xác định tính khả thi của các mẫu thiết -16.0 kế, có thể áp dụng phương pháp trình bày trong đề tài vì nó đơn giản và đủ -20.0 0.2 0.4 0.6 0.8 0 1 chính xác cho việc phân tích phản ứng Vehicle position (xv/L) động. Nên dùng kỹ thuật xếp chồng 1000 các mode dao động đối với việc phân Dynamic 0 tích như trên đặc biệt khi phân tích Bending moment (kNm) Static các mô hình cầu lớn, nhiều cấp độ. -1000 Với hầu hết các trường hợp, có thể đưa ra được các kết quả đủ chính xác -2000 bằng cách sử dụng 25 đến 30 dạng dao động đầu tiên. -3000 Phản ứng động của cầu chỉ thể hiện chính xác khi có xét đến độ gồ -4000 ghề của mặt cầu, tương tác cầu - xe -5000 1 cộ, hệ số cản của cầu và dao động của 0.2 0.4 0.6 0.8 0 cáp. Vehicle position (xv/L) Trong đề tài, chỉ trình bày tổng quan nghiên cứu đầu tiên về vấn đề tải trọng động của CDV với việc sử dụng mô hình, kỹ thuật phân tích đơn giản. Phương pháp phần tử hữu hạn phi tuyến mới chính là nghiên cứu tổng quát nhất về vấn đề này. TÀI LIỆU THAM KHẢO Nguyễn Văn Khang (1998), Dao động kỹ thuật , Nxb Khoa học và Kỹ thuật, Hà Nội. [1] Nguyễn Hoài Sơn, Đỗ Thanh Việt, Bùi Xuân Lâm (2002), Ứng dụng Matlab trong tính [2] toán kỹ thuật, Nxb Đại học Quốc gia, Tp Hồ Chí Minh. Lê Đình Tâm, Phạm Duy Hòa (2000), Cầu dây văng, Nxb Khoa học kỹ thuật, Hà Nội. [3] Nguyễn Xuân Toản (2008), Nghiên cứu xây dựng phần mềm phân tích tương tác động lực [4] học giữa cầu dây văng và đoàn tải trọng di động mô hình 2 khối lượng, Luận án Tiến sỹ Kỹ thuật, Đà Nẵng. Mai Lựu (2003), Phân tích tĩnh và động trong cầu dây văng, Luận văn Thạc sĩ, Tp Hồ [5] Chí Minh. [6] Raid Karoumi (1998), Response of Cable-Stayed and Suspension Bridges to Moving Vehicles, Doctoral Thesis, Stockholm. [7] Troitsky M.S (1988), Cable-Stayed Bridges, 2nd ed, BSP Professional Books, London. 191

CÓ THỂ BẠN MUỐN DOWNLOAD

LV.11: Bộ Luận Văn Tốt Nghiệp Chuyên Ngành Tài Chính Ngân Hàng

172 tài liệu

842 lượt tải

THÔNG TIN

TRỢ GIÚP

HỖ TRỢ KHÁCH HÀNG

Theo dõi chúng tôi

Chịu trách nhiệm nội dung:

Nguyễn Công Hà - Giám đốc Công ty TNHH TÀI LIỆU TRỰC TUYẾN VI NA

LIÊN HỆ

Địa chỉ: P402, 54A Nơ Trang Long, Phường 14, Q.Bình Thạnh, TP.HCM

Hotline: 093 303 0098

Email: support@tailieu.vn

Giấy phép Mạng Xã Hội số: 670/GP-BTTTT cấp ngày 30/11/2015 Copyright © 2022-2032 TaiLieu.VN. All rights reserved.