Luận án Dao động của cọc đối với công trình biển

Chương 1 TẢI TRỌNG TÁC ĐỘNG 1. Giơi thiệu : Tất cả những công trình thực tế luôn chịu sự tác động của tải trọng động, phụ thuộc vào thời gian làm gia tăng đáp ứng động . Vấn đề quan trọng nhất đối với những công trình biển là tải trọng động do sóng gây ra . 2. Phương trình Morison : Để tính toán tải trọng do sóng lên kết cấu cứng phải thừa nhận giả thuyết tải trọng của sóng là hàm tuyến tính của tổng lực cản và lực quán tính khi dòng chảy xuyên qua công trình . Hợp lực đối với chiều dài vi phân ds của thanh hình trụ nằm trong chất lỏng phương trình Morison cho bởi : Trong đó : F – Hợp lực tác dụng lên hình trụ h - Chiều cao mực nước tức thời r - Trọng lượng riêng của nước U – Vận tốc của phần tử nước - Gia tốc của phần tử nước, vuông góc với trục của phần tử kết cấu . D – Chiều rộng hoặc đường kính của mặt cắt hình trụ A – Diện tích mặt cắt ngang của mặt cắt hình trụ ds – Chiều dài vi phân của hình trụ theo phương đứng CD – Hệ số cản Cm - Hệ số quán tính Để tìm được hợp lực F cần phải xác định các đặc trưng về vận tốc, gia tốc của phần tử nước, các hệ số cản CD và hệ số quán tính Cm . 3 . Tải trọng sóng : Tải trọng sóng được tính toán từ phương trình chuyển động của Morison , nếu giả thuyết rằng chuyển động không ảnh hưởng đến bản thân kết cấu . Điều này có nghĩa là đặc trưng kích thước của kết cấu không vượt quá 0.2 lần chiều dài sóng . Đối với những kết cấu lớn , tham số của sóng đối với kết cấu phải kể đến lý thuyết nhiễu xạ của sóng . Khi : D/L>1 điều kiện phản xạ hầu như hoàn toàn D/L>0.2 , nhiễu xạ bắt đầu gia tăng D/L<0.2 , Công thức Morison mới có ý nghĩa . D/W>0.2 , lực quán tính chiếm ưu thế D/W<0.2 , lực cản chiếm ưu thế Trong đó : D – Chiều rộng hoặc đường kính cấu kiện L – Chiều dài sóng W – Chiều rộng quỹ đạo hạt nước cho bởi Trong đó : H – chiều cao sóng d – chiều sâu nước Một số giả thuyết khi sử dụng phương trình Morison : 1 – Vận tốc và gia tốc tức thời theo lý thuyết sóng tuyến tính và kích thước của kết cấu không ảnh hưởng đặc trưng của sóng . Giới hạn kích thước của kết cấu để sử dụng phương trình Morison là : D/L£0.2 Ở đây : D – Chiều rộng các thành phần kết cấu L – Chiều dài sóng Chiều dài của sóng được xác định từ các đặc trưng của sóng như chiều cao sóng H, chu kỳ sóng T và chiều sâu nước d. 2 – Hệ số CD và Cm xác định từ thí nghiệm . Thành phần lực cản là do lưu chất tác dụng lên công trình và lực cản được xác định từ dòng chảy đều . Hệ số cản phụ thuộc vào hệ số Reynold . Trong thực hành giá trị của Reynold được lấy là giá trị trung bình và được dùng tính toán tại mọi điểm của sóng . Hệ số Cm được lấy tuỳ thuộc vào hình dạng của kết cấu . Giá trị của Cm tra bảng Bảng 1.1 Hệ số cản của một số kết cấu thông dụng Hình dạng mặt cắt CD Hình dạng mặt cắt CD OR 2.0 1.9 r b 0.6 1.3 r b 0.5 1.3 2.0 1.3 b 1.5 0.5 Bảng 1.2 Hệ số quán tính của một số kết cấu thông dụng Dạng mặt cắt Cm 2.0 2.5 2.5 b A=D2 1.6 2.3 2.2 3 – Dạng chuẩn của phương trình Morison giả thuyết rằng kết cấu là cứng khi lực tác dụng . Tuy nhiên nếu kết cấu có đáp ứng động hoặc có một phần nổi chuyển động kích thích với vận tốc Ub, và gia tốc đối với vận tốc và gia tốc của phần tử nước . Trong trường hợp này dạng động học của phương trình có thể viết : Trong đó : Ub- Vận tốc gia tăng do mặt cắt của kết cấu - gia tốc tương ứng của mặt cắt kết cấu M – Khối lượng của mặt cắt kết cấu 4 – Phương trình Morison sử dụng giá trị CD cho lực dọc trục kết cấu và chỉ áp dụng cho những cấu kiện có lực ma sát nhỏ . 4 . Lý thuyết sóng tuyến tính Các phương trình chủ đạo và điều kiện biên được tuyến tính hoá bằng cách dùng 3 giả thiết : i) Vận tốc u và đường mặt sóng h là các giá trị nhỏ . ii ) Dòng chảy hai chiều . iii) Độ sâu không thay đổi . Phương trình mặt sóng có dạng : h(t) = asin(kx-wt) hoặc trong đó : H : chiều cao sóng , H = 2* a k: số sóng , k = 2p /L w : tần số sóng , w = 2p/T T : Chu kỳ sóng . Tốc độ sóng cho bởi : d: chiều sâu nước và Nếu thì Khi c= L/T và co = Lo/T , thì L/Lo=tanh(2p/L) Ở đây trong đó T : tính bằng giây(s) L : Tính bằng (m) Vận tốc và gia tốc theo phương đứng và phương ngang tại một điểm theo thời gian t cho bởi: Trong vùng nước sâu d/L³ 0.5, vận tốc trở thành Aùp suất dưới mặt nước cho bởi phương trình sau : Tổng năng lượng của sóng trên một đơn vị chiều rộng đỉnh sóng Lực trên một đơn vị chiều dài đỉnh sóng : 5 .Lý thuyết sóng phi tuyến : Khi chiều cao sóng tương đối lớn, không thể bỏ qua các số hạng phi tuyến trong lý thuyết sóng tuyến tính (sóng có biên độ nhỏ ). 5.1. Lý thuyết sóng Stokes : Lý thuyết sóng Stokes được áp dụng trong vùng nước sâu Phương trình mặt sóng có dạng : Với : Và : Khi đó : Với : Có thể viết h(t) dưới dạng sau : Phương trình vận tốc hạt nước cho bởi : Trong vùng nước sâu : Ở đây S’ quĩ đạo hạt nước phái trên đường mực nước S’ = y , với S= d phía trên đáy biển . 5.2 . Lý thuyết sóng đơn : Có dạng của các sóng cnoidal là hàm tuần hoàn nhưng có xu hướng trở thành một sóng không tuần hoàn khi có một đỉnh sóng k ® 1. Sóng giới hạn này được gọi là sóng đơn . d Đáy biển H h y x Sơ đồ định nghĩa của một sóng đơn Các kết quả chính của lời giải bậc 1 của lý thuyết sóng đơn Trong đó : H – Là chiều cao của sóng h - Cao trình mặt nước d – Chiều sâu nước c – Vận tốc sóng cho bỡi Vận tốc và biên độ của nước cho bỡi 5.3. Lý thuyết sóng Cnoidal: Lý thuyết sóng Cnoidal được áp dụng khi 1/50 <d/L<1/10 . Lý thuyết này quan tâm đến số hạng bậc hai,do đó cho kết quả chính xác hơn . Chi tiết lập thành bảng ở phần sau 5.4. Bảng công thức tính toán các đặc trưng của vận tốc và gia tốc sử dụng trong phương trình Morison : Tính toán cho sóng hình sin và sóng Cnoidal . BẢNG TÍNH TOÁN LỰC SÓNG Qui ước c = Pha vận tốc cg = Vận tốc nhóm E = Mật độ năng lượng S = Tổng các thành phần (=S(m)) Ef = Năng lượng của dòng Fm = Động lượng của ứng suất tán xạ Fp = Aùp lực của ứng suất tán xạ Fw = Tổng ứng suất tán xạ (= Fm+Fp) g = Gia tốc trọng trường d= chiều sâu nước từ mực nước trung bình tới đáy biển Dh = Độ hạ của sóng H = Chiều cao sóng k = Số sóng (=2p/L) K = Số hạng thứ nhất của một elliptic hoàn chỉnh K(=m) ) L = Chiều dài sóng m = Thông số của hàm elliptic m1 = Hệ số elliptic bổ sung (=1 – m) MEL = Mức năng lượng MWL = Mực nước trung bình (y= 0 tại MWL) p = Aùp lực sóng p+ = Aùp lực giới hạn (= p+rgz) T = Chu kỳ sóng u = Vận tốc hạt nước theo phương ngang U = Hệ số Stokes (=HL2/d3) v = Vận tốc hạt nước theo phương đứng a = Biên độ theo phương ngang hạt nước b = Biên độ hạt nước theo phương đứng h = Chiều cao đường mặt thoáng r = Trong lượng riêng của nước SÓNG HÌNH SIN 1 Các thông số: Cho d và T, xác định L,c Tính Tính d/Lo tìm tanh kd Xác định L =Lo tanh kd , kiểm tra * Xác định c = L/T (Cho giá trị của d,L: tìm d/Lo từ d/L trong bảng 1 và do đó xác định T) 2 Vùng nước nông : Cho Ha và T(hoặc La ) ở chiều sâu nước da 1 Xác định hb và Lb ở chiều sâu nước db Xác định L0 từ 1 Và Ha/Ho từ bảng 1 và tìm Ho Xác định Lb từ và tìm Hb , kiểm tra* *Tính U º HL2/d3 Nếu U<15 :Aùp dụng lý thuyết sóng hình Sin (ST) Nếu U>15 và d/ Lo < 0.1 : Lý thuyết sóng Cnoidal (CT) cho kết quả tin cậy hơn Nếu U>15 và d/Lo >0.1 : CT vô nghĩa và ST cũng không tin cậy BẢNG CÔNG THỨC Số hạng Đơn vị( SI) Biểu thức chung Nước sâu Nước nông (Trong khoảng ) 0<<: 0.50<<: 0<<0.050 (Trong khoảng ) 0<<: 0.50<<: 0.50<<0.050 c (m/s) cg (m/s) 0.5(c(1+G) 0.5co c E (J/m2) Ef (W/m) Fm (N/m) Fp (N/m) 0 Fw (N/m) G 0 1 Dh (m) 0 pmax (N/m) P+max-rgy P+max,0-rgy P+max (N/m) + umax (m/s) + vmax (m/s) + a (m) + b (m) + Ghi chú : + Dấu cộng sử dụng cho trường hợp y< d/2 SÓNG CNOIDAL 1 2 Các thông số 1.1 Cho d ,H và T, xác định L,c Tính H/d và T (g/d)0.5 Tìm L/d từ bảng 3 và từ đó tìm được L Xác định c = L/T 1.2 Cho d, H và L , xác định c,T Kiểm tra * tính UºHL2/d3 Tìm A từ bảng 2 Xác định c = (gd(1+AH/d))0.5và T=L/c Vùng nước nông Cho Ha và T(hoặc La ) ở chiều sâu nước da Xác định Hb và Lb ở chiều sâu nước db 1 Xác định L0 sử dụng 1.1 , hoặc T và Lo sử dụng Tính Ua và tìm Ba từ bảng 2 Tính Ho=4Ha(BaLa/Lo)0.5 Tính db/Lo và Ho/Lo Kiểm tra * : Xác định Hb/Ho từ bảng 4 và từ Hb xác định được Lb dùng 1.1 * Nếu d/Lo>0.10(hoặc d/L>0.13) lý thuyết Cnoidal ít có ý nghĩa và lý thuyết sóng hình Sin được áp dụng cho trường hợp này . CÔNG THỨC CƠ BẢN Số hạng Đơn vị A B’ c (m/s) Ef (W/m) rgH2Bc Fm (N/m) 2rgdDh Fp (N/m) rgdDh Fw (N/m) 3rgdDh Dh (m) P (N/m2) p+ - rgy p+max (N/m2) p+min (N/m2) umax (m/s) umin (m/s) U h (m) hmin+Hcn2(0,m) hmax (m) hmin+H hmin (m) q II - Tác động do dòng chảy : Trong phần này tập trung ở hai dạng dao động cơ bản phụ thuộc vào hình dạng của công trình và chế độ chảy rối mà đáp ứng động phụ thuộc chính vào dòng chảy . Lực do dòng chảy tác động lên hình trụ tròn trên một đơn vị chiều dài : Trong đó : Cm và CD là hệ số cản và hệ số nâng tương ứng . V, tương ứng là vận tốc và gia tốc của dòng chảy . 1 . Dòng chảy xung quanh hình trụ : Khi dòng chảy chuyển động bao quanh hình trụ,ở phía sau hình trụ sẽ xuất hiện một dòng xoáy . Dòng xoáy này tạo ra lực dao động theo phương vuông góc dòng chảy gọi là lực nâng FL, có tần số fL. Ngoài ra theo phương của dòng chảy có lực cản FD. Xoáy cuộn xung quanh hình trụ tròn phụ thuộc vào số Reynold,Re cho bởi : Trong đó : V- Vận tốc dòng chảy D - Kích thước mặt cắt ngang n - Hệ số nhớt . Xoáy cuộn phụ thuộc vào sự phân bố vận tốc V của dòng chảy . Và để tiện lợi sử dụng đại lượng không thứ nguyên gọi là số Strouhal, cho bỡi : Lực nâng và lực cản đặc trưng bỡi hệ số CL’,CD và CD’ Hệ số dao động của lực nâng Hệ số lực cản trung bình Hệ số lực cản Trong đó : r - Khối lượng riêng của nước A – Diện tích của mặt cắt 2 .Đáp ứng của kết cấu do xoáy cuộn : Hiện tượng cộng hưởng xảy ra khi tầng số dao động của lực nâng fL bằng với tầng số dao động tự nhiên N ,khi đó fL = N và hệ số Strouhal = 0.2 Þ fLD/V=0.2 Þ Theo phương của chuyển động fD’D/V=0.4 Þ Thông số để xác định độ lớn biên độ khi xảy ra xoáy cuộn Hệ số độ giảm LOGA của kết cấu d Hệ số giảm vận tốc V/ND Hệ số khối lượng Trong đó khối lượng hiệu quả trên một đơn vị chiều dài,xác định bỡi : Ở đây : m- Khối lượng trên một đơn vị chiều dài (gồm phần ngập trong nước và trên mặt nước) bao gồm khối lượng nước thêm vào và khối lượng nước trong lỗ rỗng L’ – Chiều cao từ vị trí ngàm của cọc đến đầu trên của kết cấu y(x) – Dạng hàm dạng x – Tính từ chiều sâu ngàm d – Chiều sâu nước ,tính từ chiều dài ngàm của cọc 3 . Qui luật dao động do xoáy cuộn gây ra : Dao động có thể theo phương của dòng chảy hoặc theo phương vuông góc với dòng chảy, phụ thuộc vào tỷ số . Dao động theo phương của dòng chảy xảy ra khi giá trị nhỏ hơn theo phương vuông góc với dòng chảy . 3.1 Dao động theo phương dòng chảy : Dao động theo phương của dòng chảy trong khu vực không ổn định có thể gia tăng khi , phụ thuộc vào tỷ số và đường cong ổn định của King –1970 . Dao động theo phương của dòng chảy khu vực mất ổn định thứ hai bắt đầu xảy ra khi . Điều kiện không có kích thích trong dạng gốc dãi cung là Tùy thuộc vào sự giảm biên độ cho cả hai khu vực mất ổn định thể hiện trên hình 1.3 của King- 1974 3.2 . Dao động theo phương ngang : Dạng dao động dãi cung sẽ bị khống chế đối với mặt cắt dạng hình tròn khi thoã mãn điều kiện Thể hiện trên hình 1.4 tùy thuộc vào 4 . Dao động không liên kết trong nhóm : Theo kết quả thực hiện của King (1975) dao động theo phương ngang dòng chảy ở hạ lưu hình trụ sẽ cản trở khi thoã mãn điều kện : 2md/rD2 >30 Sau đây là một số biểu đồ về đáp ứng của cọc trong dòng triều được trích từ sách DYNAMIC OF MARINE STRUCTRES –AIT – Asian instite of Technology –M.G.Hallam BScPhD,N.J.Heaf Beng PhD, L.R. Wootton BScPhD MICE MRAeS- 1977 Dao động Khu vực mất ổn định thứ hai Khu vực mất ổn định thứ nhất 2.5 2.0 1.5 1.0 0 0.5 1.0 1.5 2.0 2.0 1.5 1.0 Không dao động V/ND Hình 1.2 : Khu vực mất ổn định của dòng chảy là hàm của hệ số cản và khối lượng theo phương dòng chảy 2.0 1.5 1.0 Khu vực mất ổn định thứ hai giá trị cao Khu vực mất ổn định thứ nhất giá trị thấp Theo phương chuyển động 0 0.5 1.0 1.5 2.0 0 0.15 0.10 0.05 h Hình 1.3 : Biên độ đáp ứng với giá trị nhỏ hơn 3 Giá trị gần 1.9 và cao hơn 2.5 Biên độ Đường kính 2.5 2.0 1.5 1.0 0.5 0 2 4 6 8 10 12 14 16 18 Hình 1.4 : Biên độ đáp ứng là một hàm của khối lượng và hệ số cản theo phương vuông góc dòng chảy 104 105 106 107 6 5 4 3 2 0 Hình 1.5 : Dòng chảy bắt đầu dao động theo phương dòng chảy và vuông góc dòng chảy trong khu vực mất ổn định thứ hai Hệ số Renold Re Dao động bắt đầu theo phương vuông góc dòng chảy Đỉnh của dao động theo phương vuông góc dòng chảy Khu vực mất ổn định thứ hai Bắt đầu dao động theo phương dòng chảy